1. Introduction
Heavy oil is abundant but difficult to exploit. To improve the production of heavy oil reservoirs, some enhanced oil recovery (EOR) methods [
1,
2] have been adopted, such as steam flooding [
3], high-pressure nitrogen flooding [
4], CO
2 flooding [
5], fire flooding [
6], microbial technology [
7], and so on. Among them, taking the technical advantages of efficient oil displacement, low unit thermal cost, and wide applicability, the fire flooding technology is one of the most effective and promising ways to improve the thermal recovery of heavy-oil resources [
8]. By combusting the air underground, the heat and mixed gas generated push the heavy oil from the gas injection well to the production well to realize short-distance displacement of heavy oil [
9].
However, the combustion of heavy oil produces a large amount of exhaust gas, named fire-flooding exhaust, which contains both valuable and harmful components, such as methane, ethane, carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, oxygen, etc. [
10,
11].
Table 1 compares the composition of fire-flooding exhausts with other greenhouse gasses referenced from published literature [
12,
13,
14,
15,
16,
17,
18], including acid gas and flue gas. It can be seen that the composition of fire-flooding exhaust is considerably different from the other gases listed. Compared with the acid gas, the content of CO
2 and H
2S in the fire-flooding exhaust is lower. However, the latter contains much N
2 and other trace components. Compared with flue gas, being equivalent in the content of N
2 and CO
2, the fire-flooding exhaust has light hydrocarbons, H
2O, and trace amount of O
2, H
2S, H
2, and CO. It is conceivable that, with increasingly strict controls on environmental protection globally, how to efficiently and safely handle the exhaust gas which is produced by the fire flooding will become a key factor in restricting its worldwide application.
Some technologies—including burning [
19], desulfurization [
20], pressure swing adsorption [
21], CO
2 separation [
22], etc.—have been applied to effectively treat the harmful components in the fire-flooding exhaust. Unfortunately, the valuable hydrocarbon resource in it is wasted meanwhile. To improve the situation, like the acid gas injection procedure [
23], a novel treatment strategy called the ‘fire-flooding exhaust reinjection process’ was proposed in an oilfield in northwest China [
18]. In the proposed procedure, the fire-flooding exhaust will be dehydrated, compressed, and injected back into another proper oil reservoir, achieving zero emissions of harmful gases. Meanwhile, by adopting the principle of flue gas flooding [
24], the novel process will make full use of the fire-flooding resource to enhance the oil recovery of the selected oil layers. Hence, it can be regarded as a promising method to realize multiple goals simultaneously, including waste handling, carbon capture and utilization, heavy oil enhancement, etc.
In the process, the exhaust needs to undergo dehydration, compression, pipeline transportation, cooling, and reinjection. Accordingly, the pressure will gradually increase to be higher than 20 MPa from the wellhead condition, with the temperature fluctuating over 160 K meanwhile, as shown in
Figure 1 [
18]. Together with the unique composition of fire-flooding exhausts, getting the physical properties accurately during the whole operation is a major problem to be solved for reinjection engineering.
Simulation is the most economical and intuitive research method for oil and gas processing and transporting systems [
25]. For the reinjection project of the fire-flooding exhaust, it can realize the prediction of the gas physical properties under different operating conditions, thus contributes to the precise design and safe operation. There is a number of chemical simulation software nowadays—such as Aspen HYSYS, Aspen Plus, Pro II, UniSim Design, etc.—which includes a variety of physical property calculation methods to achieve the above simulation goals. The first and most important step in numerical simulation is to select a precise physical property calculation method, which should provide a good approximation of the properties over the wide process ranges of pressure, temperature and composition [
26]. Whether the prediction equation is suitable or will not directly affect the reliability of the calculation results [
27].
Some well-known equations of state (EoS)—such as the Peng–Robinson (PR) equation, the Soave–Redlich–Kwong (SRK) equation, the Lee–Kesler–Plocker (LKP) equation, the Benedict–Webb–Rubin–Starling (BWRS) equation, etc.—have been used or modified to study the natural gas treatment, acid gas treatment, reinjection, and other energy fields. Zhou et al. [
28] combined PR and Patel-Teja equation to predict the densities and viscosities of the acid gas during an acid gas reinjection process. Wang et al. [
29] developed the Peng–Robinson-cubic-plus-association (PR-CPA) EoS with a pseudo-chemical reaction approach to calculate the solubility of acid gases in aqueous alkanolamines solutions. Based on the SRK equation, Hussain [
15] established a HYSYS model to perform a feasibility analysis of recovering carbon dioxide from flue gas. The SRK equation was also used to calculate the saturation pressure of the crude oil and its mixtures with pure and impure CO
2 by Pereira et al. [
30]. Modekurti et al. [
24] conducted a dynamic modeling of a multistage CO
2 compression system, in which the LKP method was used in the property prediction for the compression section. Adom et al. [
31] applied the BWRS model to predict the compressibility factor, enthalpy, and heat leakage at various pressures for determining the factors which affect the BOG in typical LNG tanks of different capacities. Combing the PR and the eNRTL activity coefficient equation, Wan et al. [
32] conducted an acid gas-chemical solvent model to simulate the desulfurization and decarbonization process of a natural gas treatment plant.
Research results reveal that the accuracy of calculation methods for physical property forecasting closely relates to the composition and operation condition studied. A comparative study of different prediction models is helpful to master their applicability in different scenarios, which benefits in model selection. W. Yan et al. [
33] made a comparative study of non-cubic models (PC-SAFT and BWRS) and cubic models (SRK and PR) in several important aspects related to PVT modeling of reservoir fluids. The result shows that the non-cubic models are clearly advantageous in density calculation of pure components and some HP and HT fluids. Compared with the NIST results, M. Bertini et al. [
34] evaluated 10 property methods for pure and mixture of CO
2 for power cycles, including the cubic-type, virial-type, Helmholtz-type, and SAFT-type. G.D. Marcoherardino et al. [
35] conducted an analysis of the performance of five EoSs on CO
2+C
6F
6, revealing that although all equations converge at low temperature, the big difference and non-convergence occurs at higher temperatures. The region close to the critical point is the most difficult to be modeled. Varzandeh [
36] compared the GERG-2008 with other cubic and non-cubic EoS in calculation of the phase equilibrium and physical properties of natural gas related systems. A similar study was carried out by Yuan et al. [
27]. After conducting a comprehensive comparison and analysis among the SRK, PR, and GERG-2008 equations, the GERG-2008 equation was recommended as the basis for the property calculation in the natural gas liquefaction processes. Many similar studies exist that are not enumerated here.
However, according to the papers published, there is little literature studying the application of the above EoSs in predicting the physical properties of the fire-flooding exhaust. Considering the compositional singularity of the fire-flooding exhaust and the wide range of operating parameters accompanied with the reinjection procedure, the accuracy of property calculation methods is unclear yet. Since every model has its applicable boundary, the purpose of this research is to investigate the application of some commonly-used EoSs in predicting the physical properties of the fire-flooding exhaust within a wide range of pressure and temperature. The EoSs evaluated in this work include the PR, the SRK, the LKP, the BWRS, and the GERG-2008 equations. Using a HPHT physical property analyzer and the Aspen Plus software, the physical properties were detected and predicted, including the gas densities, compressibility factors, volumetric coefficients, and dew point temperatures for five exhaust compositions which were collected from a real oilfield. A comprehensive comparison and analysis between the experimental data and the calculated value for the five equations was then conducted. This work is helpful to select a precise property calculation method when simulating and designing a reinjection system which is used to reduce the emission of the harmful gas produced by the fire flooding process. It can be applied in the simulation of other waste treatment processes, too.
3. Accuracy Analysis of the Calculation Models for the Physical Properties in the Fire-Flooding Exhaust Reinjection Process
In order to find out the most precise physical equation among the five EoSs for the fire-flooding exhaust reinjection process, a comparison between the experimental data and the calculated value under the same temperature and pressure for different physical properties is needed. In this study, the properties of five representative fire-flooding exhausts collected from an oilfield in the northwest China were detected. The gas compositions are shown in
Table 2.
The experimental data, including the densities, the compressibility factors, the volumetric coefficients and the dew point temperatures for the mentioned components above, is obtained through a high pressure (HP) and high temperature (HT) physical property analyzer of which parameters are shown in
Table 3. This equipment is particularly used to detect the physical properties and phase states of formation fluids. Combining the field operating parameters and the experimental condition, the tested pressure of 4–25 MPa was set, while the tested temperature contained 303.15 K, 333.15 K, 363.15 K, and 393.15 K. After the experimental testing, a total of 1351 valid experimental data was obtained, which was used for the model evaluation and analysis.
Some indicators are necessary for evaluating the performance of the prediction models. The common used includes the correlation coefficient (
CC), the mean squared error (
MSE), the mean absolute error (
MAE), the mean absolute percentage error (
MAPE), the root mean square error (
RMSE), the normalized root mean squared error (
NRMSE), and the mean absolute range normalized error (
MARNE) [
50]. These indicators are defined as
where
xcal is the calculated value,
xexp is the experimental value,
N is the dataset size.
For different goals and data sets, the significance of every indicator—such as
MAE,
MSE, and
RMSE—is not the same [
51]. According to the similar literature published [
30,
52], in our work, the
MAPE is employed as the main criterion to evaluate the accuracy of different prediction methods.
In the following paragraphs,
Section 3.1 Section 3.2,
Section 3.3 and
Section 3.4 are the detailed comparison results between the calculated and the experimental data. Tables 4, 6, and 8 lists the comparison of gas densities, compressibility factors and volumetric coefficients for Sample 1 as an example, respectively. The comparison of dew point temperature for the five samples is illustrated from
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6. Tables 5, 7, 9, and 10 shows the prediction accuracy for each model in predicting the four properties, respectively. An overall summary was conducted in the Results.
3.1. Accuracy Analysis of the Prediction Models for the Gas Density
The density is the most important and fundamental property for process design. Using the Aspen Plus software, the densities of the collected fire-flooding exhausts were calculated by employing the five prediction models under the same conditions as the experiments. The total number of the experimental data for the five samples is 378.
Table 4 shows the experimental and the calculation results of Sample 1 as an example. By comparing with the experimental data,
Table 5 lists the prediction accuracies for each model. The comparison result shows that, the SRK equation performs the minimum overall
MAPE in calculating gas density, followed by the BWRS equation. Although the PR equation shows the highest prediction accuracy for Sample 4 and 5, its overall prediction accuracy is the lowest among the five equations. In comparison, the calculation accuracy of the SRK equation is maintained at a relatively high level for the five sample compositions as a whole. Therefore, it can be considered that the SRK equation is the suitable model for predicting the density of fire-flooding exhausts during the reinjection process.
3.2. Accuracy Analysis of the Prediction Models for the Compressibility Factor
Defined as the ratio of the actual gas volume to the ideal value, the compressibility factor, also called Z-factor, is a dimensionless parameter representing the degree to which a real gas deviates from an ideal gas. It depends on the working condition and fluid composition. Considering that the PVT relationship can be calculated when the compressibility factor is given, it is significant for prediction other gas properties.
Using the Aspen Plus software, the compressibility factors of the collected fire-flooding exhausts were predicted by employing the five prediction models under the reinjection condition. The total number of the experimental data for the five samples is 377.
Table 6 shows the experimental and the calculation results of the compressibility factors for Sample 1 as an example. By comparing with experimental data,
Table 7 lists the prediction accuracies for each model. The comparison result shows that, in the prediction of compressibility factor, the SRK equation performs the minimum overall
MAPE, followed by the BWRS equation, while that of the LKP, GERG-2008, and PR equation stands the highest. Although the PR equation shows the highest prediction accuracy for Sample 5, that of SRK equation is maintained at a relatively high level for the five sample compositions as a whole in comparison. Therefore, it can be considered that the SRK equation is the suitable model for predicting the compressibility factor of fire-flooding exhausts during the reinjection process.
3.3. Accuracy Analysis of the Prediction Models for the Volumetric Coefficient
The volumetric coefficient refers to the gas volume under the specified condition for the unit volume of gas at the standard condition. It is a function of the gas compressibility factor. Using the Aspen Plus software and the computational relation, the volumetric coefficients of the fire-flooding exhaust samples were calculated by employing the above five prediction models under the detection conditions. The total number of the experimental data for the five samples is 379.
Table 8 shows the experimental and the calculation results of Sample 1 as an example. By comparing with the experimental data,
Table 9 lists the prediction accuracies for each model. The comparison results show that the model precision is greatly affected by the gas composition. Among the evaluated methods, the SRK equation performs the highest accuracy, followed by the BWRS equation, while the PR, LKP, and GERG-2008 equations bring out bigger errors. Therefore, it can be considered that the SRK equation is suitable for predicting the volumetric coefficient of fire-flooding exhausts during the reinjection process.
3.4. Accuracy Analysis of the Prediction Models for Dew Point
Forecasting dew points involves the calculation of both gas and liquid properties, including the compressibility factor, fugacity, etc. The calculation accuracy reflects the ability of each EoS in predicting phase equilibrium parameters. Using the Aspen Plus software, the dew point lines of the fire-flooding exhaust samples were calculated by employing the above five models. The total number of the experimental data for the five samples is 217.
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6 show the experimental and predicted value of dew point temperatures under different pressures for the five samples, respectively. It can be seen that, except for the BWRS equation, the incontinuity of dew point lines occurs within some pressure range when the other four equations were adopted. According to the literature published, it stands for the non-convergence or non-solution at the specified conditions for solving the liquid-vapor equilibrium [
35,
53,
54]. The case happens for both cubic and non-cubic equations, especially near the critical zone. Improvements can be done from the aspects of model description, mathematical algorithm, etc.
Meanwhile, for every model, as the pressure increases, the deviation between the experimental and calculated value increases gradually. With the non-convergence points excluded,
Table 10 lists the prediction accuracy of the calculation models. It can be found that, the LKP equation performs the highest accuracy with a
MAPE of 3.71%, followed by SRK equation of which
MAPE is 4.33%. Errors brought by other equations are beyond 5%. Although the prediction accuracy of the dew points through the LKP equation is slightly higher than the SRK equation, the latter has stronger computational convergence than the former over the whole pressure range. Hence, considering the prediction accuracy and computational convergence, the SRK equation is suitable for predicting the dew points of fire-flooding exhaust during the reinjection process.
4. Results
Comparing all the predicted data with the experimental data, the overall summary of the calculation accuracy for each EoS is shown in
Table 11.
It can be found that, for calculating the different physical properties, the precision of the evaluated equations is uneven. The SRK equation performs well in predicting the gas densities, but the calculation of compressibility factors and volume coefficients has a notable deviation, exceeding 5%. The PR equation shows good accuracy in terms of the gas density prediction, whereas the errors of the compressibility factors and volume coefficients computation happen to be around 7%. The LKP equation has the highest precision in forecasting the dew points, but it is not accurate enough in calculating the compressibility factors and volume coefficients with the MAPE of beyond 6%. The BWRS equation has good precision in forecasting the gas densities, but the biggest error of 8.50% is observed when it is used to conduct the dew point calculation. The GERG-2008 equation also achieves good performance in the gas density prediction. However, a big error of more than 6% exists for computing the compressibility factors and the volume coefficients. Hence, it is revealed that although the GERG-2008 equation is considered as the international standard model for the natural gas, it is not entirely suitable for the fire-flooding exhaust.
Above all, for the five samples of the fire-flooding exhaust tested, the SRK equation takes the minimum MAPE in predicting the properties mentioned. In terms of dew point prediction, although the accuracy of the LKP equation is slightly higher than that of SRK equation, the latter has stronger computational convergence than the former over the whole pressure range. Taking both the prediction accuracy and computational convergence into consideration, the SRK equation is recommended as the basic prediction method for the physical properties during the reinjection of the fire-flooding exhaust.
5. Conclusions
Accurate prediction of the physical properties with the change of working conditions is of significant importance for the design and simulation of the reinjection processes of fire-flooding exhaust. In this paper, five property calculation methods—including PR, SRK, LKP, BWRS, and GERG-2008 EoS—were evaluated in predicting the properties of the fire-flooding exhaust under the conditions of reinjection process. By employing the Aspen Plus software, four properties of five representative fire-flooding exhaust compositions which collected from an oilfield in northwest China were calculated, including the gas densities, compressibility factors, volumetric coefficients, and dew point temperatures. Through a comprehensive comparison and analysis with 1351 sets of experimental detecting data, it is revealed that the SRK equation shows the best prediction performance among the five studied equations, thus it can be recommended as the basis for the property prediction in the reinjection process of the fire-flooding exhaust.
The result in this paper is helpful to select a precise property calculation method for the fire-flooding exhaust under the wide range of temperature and temperature. Thus the reinjection system can be simulated and designed in higher accuracy, which will reduce the harmful gas emission made by the in-situ combustion in a novel and efficient way. Although the compositions were collected from China, it is believed that this work can also provide reference for the study of other similar gases globally. By applying the SRK equation, the high-accuracy prediction on the physical properties and phase transition can be conducted. This is the prerequisite to carry out theoretic and applied research for the related gas. Considering the components in this study contain hydrocarbons, nitrogen, carbon dioxide, hydrogen sulfide, etc., with the increasingly higher standards on the environmental protection, the results can be also applied in the waste treatment processes.
Since more detailed components of other countries or regions have not been found out from the published literature unfortunately, our work was based on the samples collected from Chinese oilfields. Although it is believed that the samples used in this study are representative, whether the result is suitable for other oilfields can be further evaluated by the researchers based on more local component data.
Meanwhile, the binary interaction parameters of Knapp et al. were used for the SRK equation, the PR equation, the LKP equation, and the BWRS equation. The prediction accuracy can be improved through the regression of binary parameters using more experimental data.
Finally, the experimental data in this paper were obtained without a water component. There may be trace amounts of water in the fire-flooding exhaust during the reinjection process. Considering the influence of water on gas physical properties, the prediction accuracy of cubic EoSs, virial EoSs, and GERG-2008 can be further studied based on the research results of this paper.