Next Article in Journal
Energy Security in Danger? A Comparative Analysis of Oil and Copper Supply
Previous Article in Journal
Integrated Estimation of a Cyber-Physical System’s Sustainability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Prediction Models for the Physical Properties in Fire-Flooding Exhaust Reinjection Process

1
College of Petroleum Engineering, Xi’an Shiyou University, Xi’an 710065, China
2
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University), Southwest Petroleum University, Chengdu 610500, China
3
Research Institute of Engineering Technology, Xinjiang Petroleum Corporation, China National Petroleum Corporation, Karamay 834000, China
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(2), 562; https://doi.org/10.3390/en15020562
Submission received: 25 October 2021 / Revised: 1 January 2022 / Accepted: 2 January 2022 / Published: 13 January 2022

Abstract

:
The reinjection of the fire-flooding exhaust is a novel disposal process for handling the exhaust produced by the in-situ combustion technology. For reasonable process design and safe operation, it is of great significance to select an optimum property calculation method for the fire-flooding exhaust. However, due to the compositional particularity and the wide range of operating parameters during reinjection, the state equations in predicting the exhaust properties over the wide range of operating parameters have not been studied clearly yet. Hence, this paper investigates the applicability of several commonly-used equations of state, including the Soave–Redlich–Kwong equation, Peng–Robinson equation, Lee–Kesler–Plocker equation, Benedict–Webb–Rubin–Starling equation, and GERG-2008 equations. Employing Aspen Plus software, the gas densities, compressibility factors, volumetric coefficients, and dew points for five exhaust compositions are calculated. In comparison with the experimental data comprehensively, the result indicates that the Soave–Redlich–Kwong equation shows the highest precision over a wide range of temperature and pressure. The mean absolute percentage error for the above four parameters is 3.84%, 5.17%, 5.53%, and 4.33%, respectively. This study provides a reference for the accurate calculation of the physical properties of fire-flooding exhausts when designing and managing a reinjection system of fire-flooding exhaust.

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], CO2 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 CO2 and H2S in the fire-flooding exhaust is lower. However, the latter contains much N2 and other trace components. Compared with flue gas, being equivalent in the content of N2 and CO2, the fire-flooding exhaust has light hydrocarbons, H2O, and trace amount of O2, H2S, H2, 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], CO2 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 CO2 by Pereira et al. [30]. Modekurti et al. [24] conducted a dynamic modeling of a multistage CO2 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 CO2 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 CO2+C6F6, 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.

2. Equations of State to Be Evaluated

2.1. Soave–Redlich–Kwong EoS

Redlich and Kwong [37] proposed a modified cubic EoS with two parameters, which obtained satisfactory results at the pressure above the critical temperature. A modification of this equation was developed by Soave [38], thus the SRK equation was obtained. It was the first EoS to be associated with the acentric factor and could be applied to all non-polar substances [39], particularly for multicomponent mixtures. The equation is given below:
p = R T v b a T v v + b a T = α T 0.42747 R 2 T c 2 / p c b = 0.08664 R T c / p c α T = 1 + m 1 T / T c 0.5 2 m = 0.480 + 1.5741 ω 0.1715 ω 2
Where p, R, T, and v is the pressure, gas constant, temperature, and molar volume respectively, a and b are parameters of the EoS, only related to physical properties, Tc is the critical temperature, and ω is the acentric factor.
For mixtures, mixing rules are expressed as
a = i j x i y j a i a j 0.5 1 k i j b = i y i b i T c m = i x i T c i p c m = i x i p c i
where xi and xj is the molar fractions of components i and j, Tcm and pcm is the critical temperature and critical pressure of the mixture, respectively, ai and bi are the pure component parameters obtained from Equation (1), and kij is the binary interaction parameter of the SRK equation.

2.2. Peng–Robinson EoS

The SRK equation is not accurate in predicting liquid density. For hydrocarbon components (except for the methane), the predicted liquid density is generally smaller than the experimental data. To further improve the prediction accuracy of physical parameters and phase equilibrium, Peng and Robinson [40] developed a two-constant EoS based on the Soave’s model. It showed great advantages in predicting liquid phase density. The equation is given below as
p = R T v b a T v v + b + b v b a T = 0.45724 α T R 2 T c 2 / p c b T = 0.07780 R T c / p c α T = 1 + k 1 T / T c 0.5 2 k = 0.37464 + 1.5422 ω 0.26992 ω 2
where the notation declarations and mixing rules are consistent with the SRK equation.

2.3. Lee–Kesler–Plocker EoS

The LKP equation is a modified LK equation proposed by Plocker, Knapp, and Prausnitz [41]. It retained the prototype of the LK equation but modified the mixing rule. The LKP equation is applicable to the liquid phase, especially with good accuracy in the calculation of the physical properties of hydrocarbons. It was deemed to be the best method for forecasting the gas thermophysical properties [42]. The equation is expressed as
Z = Z 0 + ω ω r Z r Z 0 Z = p r v r T r = 1 + B v r + C v r 2 + D v r 5 + c 4 T r 3 v r 2 β + γ v r 2 exp γ v r B = b 1 b 2 T r b 3 T r 2 b 4 T r 3 C = c 1 c 2 T r + c 3 T r 3 D = d 1 + d 2 T r
where Z is the compression factor, 0 and r represents the simple and reference fluid parameters, respectively, and b1~b4, c1~c4, d1, d2, β, and γ are model constants.
For mixtures, mixing rules are expressed as
T c m v c m 0.25 = i j x i x j v c i j 0.25 T c i j v c m = i j x i x j v c i j ω m = i x i ω i T c i j = T c i T c j 0.5 k i j v c i j = 1 8 v c i 1 / 3 + v c j 1 / 3 3
where kij represents the binary interaction parameter of the LKP equation.

2.4. Benedict–Webb–Rubin–Starling EoS

To expand the application range and improve the accuracy at high pressure and low temperature, Benedict–Webb–Rubin conducted an 8-parameters method which was applied in a gas–liquid system in 1940, called BWR equation [43]. Compared with the SRK and PR method, the equation performed high accuracy in predicting the thermophysical properties of light hydrocarbons and their mixtures. However, it had poor adaptability to the mixtures with more non-hydrocarbon gases, heavier hydrocarbon components, and lower temperatures (Tr < 0.6). Combining the thermodynamic data of C1~C8 alkanes, Staring and Powers [44], Starling and Han [45], updated the BWR model to an 11-constant equation, namely the BWRS equation, which showed high prediction accuracy for the physical properties of hydrocarbon, CO2, H2S, and N2. The BWRS equation is given as
p = ρ R T + B 0 R T A 0 C 0 T 2 + D 0 T 3 E 0 T 4 ρ 2 + b R T a d T ρ 3 + α a + d T ρ 6 + c ρ 3 T 2 1 + γ ρ 2 exp γ ρ 2
where A0, B0, C0, D0, E0, γ, a, b, c, d, and α are the parameters of BWRS equation of which form are shown in reference [44].

2.5. GERG-2008 EoS

As an international standard named ISO20765-2 [46], the GERG-2008 equation is the extension of the GERG-2004 model [47]. This equation covers 21 natural gas components, including CH4, N2, CO2, C2H6, C3H8, C4H10, I-C4H10, C5H12, I-C5H12, C6H14, C7H16, C8H18, C9H20, C10H22, H2, O2, CO, H2O, H2S, He, and Ar. Being explicit in the Helmholtz free energy as a function of density, temperature, and composition, the GERG-2008 equation was constructed by fitting a large number of experimental data to obtain the regression coefficients. Furthermore, it has been improving and upgrading with more experimental data nowadays [48]. Hence, it showed excellent predictive performance in the gas phase, liquid phase, supercritical region, and vapor–liquid equilibrium states for mixtures of the components mentioned above, and was regarded as the most comprehensive calculation method in predicting natural gas thermophysical properties. However, for the fire-flooding exhaust, the precision has not been assessed and needs to be studied. The equation is expressed as
a ρ , T , x ¯ = a o ρ , T , x ¯ + a r ρ , T , x ¯
Using the dimensionless form of Helmholtz free energy α = a/(RT), the following equation is deduced (note that αo does not depend on δ and τ of the mixture but ρ and T)
α δ , τ , x ¯ = α o ρ , T , x ¯ + α r δ , τ , x ¯ δ = ρ / ρ r , τ = T / T r
The dimensionless form of the Helmholtz free energy for the ideal-gas mixture αo is given by
α o ρ , T , x ¯ = i = 1 N x i α o i o ρ , T + ln x i
The residual part of the Helmholtz free energy is expressed as
α r δ , τ , x ¯ = i = 1 N x i α o i r δ , τ + i = 1 N 1 j = i + 1 N x i x j F i j α i j r δ , τ
For mixtures, the mixture density and mixture temperature are given below as
1 ρ r x ¯ = i = 1 N x i 2 1 ρ c i + i = 1 N 1 j = i + 1 N 2 x i x j β v , i j γ v , i j · x i + x j β v , i j 2 x i + x j · 1 8 1 ρ c i 1 / 3 + 1 ρ c j 1 / 3 3 T r x ¯ = i = 1 N x i 2 T c i + i = 1 N 1 j = i + 1 N 2 x i x j β T , i j γ T , i j · x i + x j β T , i j 2 x i + x j T c i · T c j 0.5
The detailed expressions involved in the above equations can be seen in reference [47].
The above property calculation methods have been embedded in the Aspen Plus software which can be used as a platform to analyze the accuracy of each EoS. In Aspen Plus, the above property methods are named as SRK, PENG-ROB, LK-PLOCK, BWRS, and GERG2008. In this paper, for the SRK equation, the PR equation, the LKP equation, and the BWRS equation, the binary interaction parameters of Knapp et al. [49] were used. The precision of the five models mentioned is analyzed by comparing the calculation results with the experimental data.

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
C C = cov x cal , x exp D x cal D x exp
M S E = 1 N i = 1 N x i exp x i cal 2
M A E = 1 N i = 1 N x i exp x i cal
M A P E = 1 N i = 1 N x i cal x i exp x i exp × 100 %
R M S E = 1 N i = 1 N x i cal x i exp 2
N R M S E = 1 x ¯ exp 1 N i = 1 N x i cal x i exp 2
M A R N E = 1 N i = 1 N x i cal x i exp max x i c a l × 100 %
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.

Author Contributions

Conceptualization, Y.W.; Methodology, Y.W. and M.C.; Validation, M.C., D.H. and S.F.; Formal analysis, Y.W. and M.C.; Resources, L.C.; Writing—original draft preparation, Y.W. and M.C.; Writing—review and editing, Y.W.; Supervision, L.C. and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Science Basic Research Program of Shaanxi (Program No. 2021JQ-593), and also supported by Open Fund (PLN2020-16) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

EoSequation of state
PRPeng–Robinson
SRKSoave–Redlich–Kwong
LKPLee–Kesler–Plocker
BWRSBenedict–Webb–Rubin–Starling
PR-CPAPeng–Robinson cubic plus association
PC-SAFTperturbed chain-statistical associating fluid theory
NRTLnon-random two liquid
NISTNational Institute of Standards and Technology
PVTpressure, volume, and temperature
LNGliquefied natural gas
BOGboil-off gas
HPHThigh pressure and high temperature
CCcorrelation coefficient
MSEmean squared error
MAEmean absolute error
MAPEmean absolute percentage error
RMSEroot mean square error
NRMSEnormalized root mean squared error
MARNEmean absolute range normalized error

References

  1. Le Van, S.; Chon, B.H. Chemical Flooding in Heavy-Oil Reservoirs: From Technical Investigation to Optimization Using Response Surface Methodology. Energies 2016, 9, 711. [Google Scholar] [CrossRef] [Green Version]
  2. Alvarado, V.; Manrique, E. Enhanced Oil Recovery: An Update Review. Energies 2010, 3, 1529–1575. [Google Scholar] [CrossRef]
  3. Kar, T.; Hascakir, B. Effect of solvent type on emulsion formation in steam and solvent-steam flooding processes for heavy oil recovery. Colloids Surf. A Physicochem. Eng. Asp. 2021, 611, 125783. [Google Scholar] [CrossRef]
  4. Yang, J.; Hou, J.; Qu, M.; Liang, T.; Wen, Y. Experimental study the flow behaviors and mechanisms of nitrogen and foam assisted nitrogen gas flooding in 2-D visualized fractured-vuggy model. J. Pet. Sci. Eng. 2020, 194, 107501. [Google Scholar] [CrossRef]
  5. Thorne, R.J.; Sundseth, K.; Bouman, E.; Czarnowska, L.; Mathisen, A.; Skagestad, R.; Stanek, W.; Pacyna, J.M.; Pacyna, E.G. Technical and environmental viability of a European CO2 EOR system. Int. J. Greenh. Gas Control 2020, 92, 102857. [Google Scholar] [CrossRef]
  6. Guan, W.; Xi, C.; Chen, Y.; Zhang, X.; Muhetar; Liang, J.; Huang, J.; Wu, J. Fire-flooding technologies in post-steam-injected heavy oil reservoirs. Pet. Explor. Dev. 2011, 38, 452–463. [Google Scholar] [CrossRef]
  7. Wang, T.; Xiu, J.; Cui, Q.; Huang, L.; Ma, Y.; Yu, L. Research advances in numerical simulation of microbial enhanced oil recovery. J. Cent. South Univ. (Sci. Technol.) 2019, 50, 1474–1484. [Google Scholar] [CrossRef]
  8. Yuan, S.; Jiang, H.; Yang, F.; Shi, Y.; Bai, Y.; Du, K. Research on characteristics of fire flooding zones based on core analysis. J. Pet. Sci. Eng. 2018, 170, 607–610. [Google Scholar] [CrossRef]
  9. Mendoza, A.; Perozo, H.; Oliveros, D.A.; Reyes, N. Study for improving the in-situ combustion of Venezuela extra-heavy crude oil using unconventional additive. In Proceedings of the World Heavy Oil Congress 2011, Edmonton, AB, Canada, 14–17 March 2011. [Google Scholar]
  10. Feng, T.; Wu, Y.; Hu, C.; Gong, Y.; Yang, S. Technology of fire flooding process control: Application to fire flooding in a deep and thick heavy oil reservoir in China. In Proceedings of the SPE Asia Pacific Oil and Gas Conference and Exhibition, Brisbane, Australia, 23–25 October 2018. SPE-191897-MS. [Google Scholar] [CrossRef]
  11. Zhao, R.; Li, B.; Yang, J.; Gao, B.; Zhang, C.; Fu, J.Y.; Wang, R. Temperature prediction via reaction heat calculation of burned pseudo-components during in-situ combustion. Fuel 2018, 214, 264–271. [Google Scholar] [CrossRef]
  12. Wang, S.; Tang, L.; Carroll, J.J. A state of the art of acid gas injection and its application prospects in China. Nat. Gas Ind. 2013, 33, 105–111. [Google Scholar] [CrossRef]
  13. Dong, M.; Huang, S. Flue Gas Injection for Heavy Oil Recovery. J. Can. Pet. Technol. 2002, 41, 44–50. [Google Scholar] [CrossRef]
  14. Chen, D.; Zhang, Z.; Yang, J.; Ma, X.; Li, P.; Hao, X.; Guan, G. Process simulation and energy saving analysis of CO2 capture by chemical absorption method based on self-heat recuperation. J. Chem. Ind. Eng. 2019, 70, 2938. [Google Scholar] [CrossRef]
  15. Hussain, A.; Hägg, M.-B. A feasibility study of CO2 capture from flue gas by a facilitated transport membrane. J. Membr. Sci. 2010, 359, 140–148. [Google Scholar] [CrossRef]
  16. Aouini, I.; Ledoux, A.; Estel, L.; Mary, S.; Evrard, P.; Valognes, B. Experimental Study of Carbon Dioxide Capture from Synthetic Industrial Incinerator Flue Gas with a Pilot and Laboratory Measurements. Procedia Eng. 2012, 42, 704–720. [Google Scholar] [CrossRef] [Green Version]
  17. Qin, H. Analysis and control measures of tail gas in the development process of the in situ combustion. Oil Gas Field Surf. Eng. 2016, 35, 102–104. [Google Scholar] [CrossRef]
  18. Wang, Y.; Wang, S.; Chen, L.; Chen, L.; Wang, S. Numerical calculation and analysis on phase behavior of fireflood exhaust from HQ-1 well in Xinjiang oilfield. IOP Conf. Ser. Earth Environ. Sci. 2019, 242, 022024. [Google Scholar] [CrossRef]
  19. Askari, O.; Vien, K.; Wang, Z.; Sirio, M.; Metghalchi, H. Exhaust gas recirculation effects on flame structure and laminar burning speeds of H2/CO/air flames at high pressures and temperatures. Appl. Energy 2016, 179, 451–462. [Google Scholar] [CrossRef]
  20. Kang, M.S.; Hwang, J. Effect of SO2 Concentration on NOx Removal Efficiency in NaOH-Based Wet Scrubbing. J. Korean Soc. Atmos. Environ. 2018, 34, 659–667. [Google Scholar] [CrossRef]
  21. Majchrzak-Kucęba, I.; Wawrzyńczak, D.; Zdeb, J.; Smółka, W.; Zajchowski, A. Treatment of Flue Gas in a CO2 Capture Pilot Plant for a Commercial CFB Boiler. Energies 2021, 14, 2458. [Google Scholar] [CrossRef]
  22. Razmara, N.; Kirch, A.; Meneghini, J.; Miranda, C. Efficient CH4/CO2 Gas Mixture Separation through Nanoporous Graphene Membrane Designs. Energies 2021, 14, 2488. [Google Scholar] [CrossRef]
  23. Li, Q.; Liu, X.; Du, L.; Bai, B.; Fang, Z.; Jing, M.; Li, X. Economics of Acid Gas Injection with Comparison to Sulfur Recovery in China. Energy Procedia 2013, 37, 2505–2510. [Google Scholar] [CrossRef] [Green Version]
  24. Li, S.; Li, Z.; Sun, X. Effect of flue gas and n -hexane on heavy oil properties in steam flooding process. Fuel 2017, 187, 84–93. [Google Scholar] [CrossRef]
  25. Wang, Y.; Wei, N.; Wan, D.; Wang, S.; Yuan, Z. Numerical Simulation for Preheating New Submarine Hot Oil Pipelines. Energies 2019, 12, 3518. [Google Scholar] [CrossRef] [Green Version]
  26. Modekurti, S.; Eslick, J.; Omell, B.; Bhattacharyya, D.; Miller, D.; Zitney, S.E. Design, dynamic modeling, and control of a multistage CO2 compression system. Int. J. Greenh. Gas Control 2017, 62, 31–45. [Google Scholar] [CrossRef] [Green Version]
  27. Yuan, Z.; Cui, M.; Song, R.; Xie, Y. Evaluation of prediction models for the physical parameters in natural gas liquefaction processes. J. Nat. Gas Sci. Eng. 2015, 27, 876–886. [Google Scholar] [CrossRef]
  28. Zhou, J.; Chen, J.; Gao, P.; Gao, H. Wellbore flow model of acid gas reinjection. Nat. Gas Ind. 2005, 8, 76–78. [Google Scholar] [CrossRef]
  29. Wang, T.; El Ahmar, E.; Coquelet, C.; Kontogeorgis, G.M. Improvement of the PR-CPA equation of state for modelling of acid gases solubilities in aqueous alkanolamine solutions. Fluid Phase Equilib. 2018, 471, 74–87. [Google Scholar] [CrossRef] [Green Version]
  30. Pereira, V.J.; Regueira, V.B.; Costa, G.M.N.; De Melo, S.A.B.V. Modeling the Saturation Pressure of Systems Containing Crude Oils and CO2 Using the SRK Equation of State. J. Chem. Eng. Data 2019, 64, 2134–2142. [Google Scholar] [CrossRef]
  31. Adom, E.; Islam, S.Z.; Ji, X. Modelling of boil-off gas in LNG tanks: A case study. Int. J. Eng. Technol. 2010, 2, 292–296. [Google Scholar]
  32. Wan, P.; Shi, P.; Dai, Y. Simulation and optimization technology of the hydrocarbons removal unit of acid gas removal. Comput. Appl. Chem. 2020, 3, 201–208. [Google Scholar] [CrossRef]
  33. Yan, W.; Varzandeh, F.; Stenby, E.H. PVT modeling of reservoir fluids using PC-SAFT EoS and Soave-BWR EoS. Fluid Phase Equilib. 2015, 386, 96–124. [Google Scholar] [CrossRef]
  34. Bertini, M.; Fiaschi, D.; Manfrida, G.; Niknam, P.H.; Talluri, L. Evaluation of the property methods for pure and mixture of CO2 for power cycles analysis. Energy Convers. Manag. 2021, 245, 114568. [Google Scholar] [CrossRef]
  35. Di Marcoberardino, G.; Morosini, E.; Manzolini, G. Preliminary investigation of the influence of equations of state on the performance of CO2 + C6F6 as innovative working fluid in transcritical cycles. Energy 2022, 238, 121815. [Google Scholar] [CrossRef]
  36. Varzandeh, F.; Stenby, E.H.; Yan, W. Comparison of GERG-2008 and simpler EoS models in calculation of phase equilibrium and physical properties of natural gas related systems. Fluid Phase Equilib. 2017, 434, 21–43. [Google Scholar] [CrossRef] [Green Version]
  37. Redlich, O.; Kwong, J.N.S. On the Thermodynamics of Solutions. V. An Equation of State. Fugacities of Gaseous Solutions. Chem. Rev. 1949, 44, 233–244. [Google Scholar] [CrossRef] [PubMed]
  38. Soave, G. Equilibrium constants from a modified Redlich-Kwong equation of state. Chem. Eng. Sci. 1972, 27, 1197–1203. [Google Scholar] [CrossRef]
  39. Zhong, J.; Zhao, X.; Li, B.; Lin, X.; Geng, Z. Calculation of material properties of SF6/CF4 gas mixture based on SRK state equation. High Voltage Appar. 2016, 52, 48–53. [Google Scholar] [CrossRef]
  40. Peng, D.-Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  41. Plocker, U.; Knapp, H.; Prausnitz, J. Calculation of High-Pressure Vapor-Liquid Equilibria from a Corresponding-States Correlation with Emphasis on Asymmetric Mixtures. Ind. Eng. Chem. Process. Des. Dev. 1978, 17, 324–332. [Google Scholar] [CrossRef]
  42. Liu, L.; Liu, B.; Li, S.; Qi, J. Calculation of LNG Thermodynamic Properties. Contemp. Chem. Ind. 2010, 39, 88–90. [Google Scholar] [CrossRef]
  43. Benedict, M.; Webb, G.B.; Rubin, L.C. An Empirical Equation for Thermodynamic Properties of Light Hydrocarbons and Their Mixtures, I. Methane, Ethane, Propane and n-Butane. J. Chem. Phys. 1940, 8, 334–345. [Google Scholar] [CrossRef]
  44. Starling, K.E.; Powers, J.E. Enthalpy of Mixtures by Modified BWR equation. Ind. Eng. Chem. Fundam. 1970, 9, 531–537. [Google Scholar] [CrossRef]
  45. Starling, K.E.; Han, M.S. Thermo data refined for LPG-14 mixtures. Hydrocarbon Process. 1972, 51, 129–132. [Google Scholar]
  46. Farzaneh-Gord, M.; Mohseni-Gharyehsafa, B.; Toikka, A.; Zvereva, I. Sensitivity of natural gas flow measurement to AGA8 or GERG2008 equation of state utilization. J. Nat. Gas Sci. Eng. 2018, 57, 305–321. [Google Scholar] [CrossRef]
  47. Kunz, O.; Wagner, W. The GERG-2008 Wide-Range Equation of State for Natural Gases and Other Mixtures: An Expansion of GERG-2004. J. Chem. Eng. Data 2012, 57, 3032–3091. [Google Scholar] [CrossRef]
  48. Gong, K.; Wang, Z.; Jia, Y. Generalization of thermodynamic research of LNG and other mixtures storage and transportation process. Sci. Technol. Eng. 2013, 13, 121–131, CNKI:SUN:KXJS.0.2013-35-023. [Google Scholar]
  49. Renon, H. Vapor-Liquid Equilibria for Mixtures of Low Boiling Substances, Vol. VI: H. Knapp, R. Döring, L. Oellrich, V. Plöcker and J.M. Prausnitz, 1982, 9890 pp. ISBN 3-921 567-38-6, DECHEMA, Chemistry Data Series, Frankfurt am Main. Fluid Phase Equilib. 1985, 21, 174–175. [Google Scholar] [CrossRef]
  50. Panapakidis, I.P.; Dagoumas, A.S. Day-ahead natural gas demand forecasting based on the combination of wavelet transform and ANFIS/genetic algorithm/neural network model. Energy 2017, 118, 231–245. [Google Scholar] [CrossRef]
  51. Wei, N.; Li, C.; Peng, X.; Zeng, F.; Lu, X. Conventional models and artificial intelligence-based models for energy consumption forecasting: A review. J. Pet. Sci. Eng. 2019, 181, 106187. [Google Scholar] [CrossRef]
  52. Sedaghat, M.; Kiomarsiyan, A. Applying MLP-ANN as a novel and accurate method to estimate gas density. Pet. Sci. Technol. 2019, 37, 2128–2133. [Google Scholar] [CrossRef]
  53. Panah, H.S. Modeling binary vapor–liquid equilibrium data containing perfluorocarbons using the Peng–Robinson and the PC-SAFT equations of state. Int. J. Refrig. 2018, 85, 13–26. [Google Scholar] [CrossRef]
  54. Bisotti, F.; Galeazzi, A.; Galatioto, L.; Masserdotti, F.; Bigi, A.; Gritti, P.; Manenti, F. Implementing robust thermodynamic model for reliable bubble/dew problem solution in cryogenic distillation of air separation units. Int. J. Thermofluids 2021, 10, 100083. [Google Scholar] [CrossRef]
Figure 1. Temperature and pressure change during the proposed reinjection process of the fire-flooding exhaust.
Figure 1. Temperature and pressure change during the proposed reinjection process of the fire-flooding exhaust.
Energies 15 00562 g001
Figure 2. Experimental and calculation results of dew points for Sample 1.
Figure 2. Experimental and calculation results of dew points for Sample 1.
Energies 15 00562 g002
Figure 3. Experimental and calculation results of dew points for Sample 2.
Figure 3. Experimental and calculation results of dew points for Sample 2.
Energies 15 00562 g003
Figure 4. Experimental and calculation results of dew points for Sample 3.
Figure 4. Experimental and calculation results of dew points for Sample 3.
Energies 15 00562 g004
Figure 5. Experimental and calculation results of dew points for Sample 4.
Figure 5. Experimental and calculation results of dew points for Sample 4.
Energies 15 00562 g005
Figure 6. Experimental and calculation results of dew points for Sample 5.
Figure 6. Experimental and calculation results of dew points for Sample 5.
Energies 15 00562 g006
Table 1. Comparison of the components of acid gas, flue gas, and fire-flooding exhaust (unit: mol %).
Table 1. Comparison of the components of acid gas, flue gas, and fire-flooding exhaust (unit: mol %).
GasCO2N2O2H2H2OH2SC1C2OthersNotes
Acid gas No.186 10 Acheson(Canada) [12]
Acid gas No.265 35 Lisbon (America) [12]
Acid gas No.3100 Salah(Algeria) [12]
Flue gas No.115.684.4 Combustion gas in a thermal power plant [13]
Flue gas No.214.280.85 Reference [13]
Flue gas No.318.4175.284.71 1.6 Reference [14]
Flue gas No.410805 5 Reference [15]
Flue gas No.55–949–578–12 saturation Reference [16]
Fire-flooding exhaust No.115.2764.851.14 17.61 Liaohe Oilfield, China [17]
Fire-flooding exhaust No.213.2973.60.35 11.870.36C3+: 0.53
Fire-flooding exhaust No.313.5676.760.970.492.580.14.380.38C3+: 1.01, CO: 0.15Xinjiang Oilfield, China [18]
Table 2. Compositions of five fire-flooding exhaust samples.
Table 2. Compositions of five fire-flooding exhaust samples.
Sample NameComposition (mol%)Analysis Method
CH4C2H6C3H8C4+N2CO2O2H2H2S
Sample 11.90.140.070.0878.1317.592.07-0.02PVT property analysis
Sample 23.080.220.110.07480.415.20.250.530.088
Sample 34.490.290.130.1677.7613.933.27-0.062
Sample 45.610.360.160.279.2213.981.15-0.088
Sample 50.970.0950.040.03481.8916.973.5-0.011
Table 3. Equipment parameters used in this study.
Table 3. Equipment parameters used in this study.
Technical IndicatorsValue
Working pressure0–100 MPa
Working temperatureFrom ambient temperature to 473.15 K
Pump capacity500 mL
Pump precision0.0001 mL
Table 4. Experimental data and calculated results of the gas density for Sample 1.
Table 4. Experimental data and calculated results of the gas density for Sample 1.
Experimental DataCalculated Results
SRKPRLKPBWRSGERG-2008
T (K)p (MPa)ρ (kg/m3)ρ (kg/m3)
303.1524.040280.818287.1567303.4508296.5895285.2223295.9078
22.020260.102265.6166280.0517274.2284264.7271273.7772
20.000238.235245.3392258.0877253.0746245.2787252.7865
18.047217.519222.6385233.5732229.3109223.3094229.1533
16.027194.501199.0326208.1689204.5558200.2309204.488
14.007171.483175.4618182.899179.8466176.9487179.8341
11.919147.315151.2386157.043154.5195152.7901154.5381
9.899121.995125.6173129.8411127.8607127.0213127.8929
8.01397.826101.4431104.3418102.8719102.5618102.9049
5.99372.50675.3370677.0261776.0928476.0689776.11882
3.97348.33850.1958850.9859550.513350.5658250.52896
333.1524.040242.839252.8345265.6839259.5303251.8671259.5332
22.020224.425234.3723245.7615240.4981234.2006240.5679
20.067207.161215.8794225.8603221.3727216.3539221.4786
18.047188.747196.0967204.6345200.8674197.0882200.9854
16.027169.182175.68182.8025179.6819177.015179.7921
14.007149.616154.6753160.4283157.8923156.1701157.9821
11.987127.749133.1463137.5982135.5974134.6212135.6611
9.966105.882111.1749114.4223112.9163112.4723112.9546
7.94685.16688.8586991.0321189.9819789.8630490.00016
5.99363.29967.0642668.360267.7045867.7319867.7106
3.90640.28143.6512744.2278243.9203543.975543.92078
363.1524.040218.670227.4674237.9229232.4172227.0283232.7277
22.088202.558211.2601220.5309215.8245211.4446216.1178
20.067186.445193.9769202.0344198.0834194.6852198.3424
17.980169.182175.5751182.4016179.1564176.6794179.3675
16.027153.069157.8771163.5837160.933159.2082161.0933
14.007134.655139.1037143.6979141.5992140.5207141.7069
11.987115.090119.8922123.4378121.8317121.2511121.8923
9.96695.524100.2885102.87101.7001101.4656101.7242
7.94674.80880.3476682.0741781.2846681.2533481.28522
5.99356.39460.8124961.8425461.3638161.4152861.35442
4.51241.43245.8635646.4680146.182446.2422446.17241
393.1524.040197.954207.4401216.1139211.2069207.3155211.6463
22.088182.992192.5665200.2326196.0667192.9535196.4434
20.067168.031176.7583183.399179.9353177.5635180.2405
18.047153.069160.5257166.1658163.3382161.627163.5694
16.027136.957143.878148.5513146.294145.148146.4537
14.007120.844126.8295130.5811128.8292128.1415128.9249
11.987104.731109.4004112.2884110.9786110.6362111.0221
10.03487.46892.215594.3390393.397793.2821393.40498
8.01370.20574.119775.5428974.9228874.9400374.90869
5.99351.79055.739656.5754556.2150956.279556.19466
4.71440.28143.9716944.5046944.2755644.3377744.25744
Table 5. MAPEs of experimental and calculated results in the gas density prediction.
Table 5. MAPEs of experimental and calculated results in the gas density prediction.
MAPESRKPRLKPBWRSGERG-2008Number of Experimental Data Sets
Sample 14.62%8.13%6.50%5.32%6.54%76
Sample 22.92%6.21%4.62%3.46%4.66%77
Sample 33.91%7.33%5.61%4.46%5.68%75
Sample 44.95%2.14%3.50%4.59%3.45%76
Sample 52.82%0.80%1.28%2.39%1.23%74
Overall MAPE3.84%4.92%4.30%4.04%4.31%-
Table 6. Experimental data and calculated results of the compressibility factors for Sample 1.
Table 6. Experimental data and calculated results of the compressibility factors for Sample 1.
Experimental DataCalculated Results
SRKPRLKPBWRSGERG-2008
T (K)p (MPa)Z-FactorZ-Factor
303.1523.9371.0681.0232770.96853290.99074891.0298850.9929652
21.9241.0581.0097670.95782030.97808671.0130050.9796722
19.9111.0460.99781640.9487610.96743120.9977470.9684879
17.8971.0370.98763240.94160310.95902990.98442970.9596604
15.8841.0280.97942670.93661470.95311640.97341180.9534156
13.9821.0220.97367770.93414570.94998960.96545920.9500521
11.8571.0160.9697460.93425270.94945190.95978420.9493242
9.8431.0160.96859810.93740250.95185590.95786990.95161
8.0541.0170.96975940.94286330.95632270.95918820.9560161
6.0401.0230.97354730.95213140.96383590.96416020.9635069
4.0271.0260.97992690.9647710.97378390.9727740.9734824
333.1523.8261.1191.0533051.0025811.026151.0570011.026104
21.8121.1071.0409270.9929161.0144691.041361.01415
19.9111.0921.0302470.9849081.0047521.0277521.00426
17.8971.0801.0201250.9777610.9959941.014790.995405
15.9961.0701.0117910.9724140.9892821.0041240.988676
13.9821.0611.004380.9684030.9839410.9947450.983383
12.0811.0560.9988140.9663260.9806440.9879190.980178
9.8431.0530.9941670.9662080.9790050.9826930.978681
8.0541.0500.9919990.9680620.9794550.9808680.97925
6.0401.0520.9912630.9723470.9818140.9814480.981725
4.0271.0520.9923570.979060.9860850.9848680.986072
363.1523.8261.1461.0741541.0271681.0512831.0758961.049874
21.9241.1331.0629641.0184571.0405091.0617921.039098
20.0221.1181.0524761.0105541.030671.0485851.029326
18.0091.1061.0422051.0031621.0213581.0357261.020153
15.9961.0941.0328670.996881.0132731.0241991.012268
13.9821.0861.024540.9918221.0065021.014191.00574
11.9691.0801.0172980.9881081.001121.005891.000624
9.9551.0741.0112140.9858590.997190.9994830.996956
7.9421.0701.0063530.9851930.9947580.9951380.994752
5.9281.0641.0027710.9862230.9938510.9929890.994007
4.2511.0621.0007950.9884510.9942650.9929420.994485
393.1524.0491.1761.0893541.0456241.0699271.0900221.067705
21.9241.1601.0772031.0361351.0579811.0748141.055969
20.0221.1451.0668711.0282921.0480421.0619761.046273
18.0091.1331.0565361.0207221.0383551.0492951.036894
15.9961.1241.0468741.0139831.0295991.0376871.028482
13.9821.1131.0379381.0081551.0218421.0272951.021088
11.9691.1041.0297791.0033221.0151491.0182711.014754
9.9551.0971.0224510.9995751.009581.0107641.009513
7.9421.0911.0160040.9970021.0051881.0049151.005386
6.0401.0851.0107660.9957331.002161.0010271.002521
4.8101.0801.0078350.9955471.0007930.99941.001202
Table 7. MAPEs of experimental and calculated results in the compressibility factor prediction.
Table 7. MAPEs of experimental and calculated results in the compressibility factor prediction.
MAPESRKPRLKPBWRSGERG-2008Number of Experimental Data Sets
Sample 15.72%8.78%7.39%6.34%7.42%76
Sample 29.48%12.35%11.00%9.94%11.03%76
Sample 33.78%6.86%5.34%4.32%5.41%75
Sample 41.92%4.94%3.44%2.45%3.49%76
Sample 54.93%1.51%3.11%4.27%3.06%74
Overall MAPE5.17%6.89%6.06%5.46%6.08%-
Table 8. Experimental data and calculated results of the volumetric coefficient for Sample 1.
Table 8. Experimental data and calculated results of the volumetric coefficient for Sample 1.
Experimental DataCalculated Results
SRKPRLKPBWRSGERG-2008
T (K)p (MPa)Volumetric CoefficientVolumetric Coefficient
303.1524.1420.0047000.00444960.00421070.00430810.00447980.004318
22.0180.0050360.00481050.00456250.00465940.00482660.0046671
20.1060.0054560.00520810.00495080.00504890.00520930.0050546
18.0530.0060430.00573910.00547040.00557210.00572180.0055759
16.0000.0067150.00641980.0061380.00624640.00638130.0062485
14.0180.0076380.00728220.00698610.00710460.0072210.0071051
12.0350.0088970.00844850.00813630.00826910.00836270.0082681
9.9820.0107430.01017170.00984080.00999320.01005930.0099907
8.0710.0133450.01259570.01224570.01242070.01245830.0124167
6.0180.0178780.01696030.01658840.01679190.01679720.0167862
4.0350.0269420.02545510.02506070.02529520.0252690.0252873
333.1524.1420.0052880.0050360.0047920.0049060.0050560.004906
22.0890.0057910.0054380.0051850.0052990.0054420.005297
20.1770.0062110.0058910.0056290.0057440.0058790.005741
18.1240.0068820.0064910.006220.0063370.0064590.006333
16.0710.0077220.0072550.0069720.0070930.0072010.007089
14.0890.0087290.0082160.007920.0080480.0081380.008043
12.0350.0101560.009560.0092490.0093860.0094550.009382
9.9820.0121700.0114760.011150.0112990.0113440.011295
8.0710.0151080.014160.0138180.0139810.0140010.013978
6.0890.0201440.0187560.0183960.0185760.018570.018574
4.0350.0303840.028330.0279490.028150.0281160.02815
363.1524.1420.0059590.0055970.0053510.0054780.0056090.005471
22.0890.0064630.0060490.0057940.0059210.0060430.005913
20.1060.0069660.0065760.0063140.006440.0065530.006431
18.0530.0077220.0072510.0069790.0071060.0072060.007097
16.1420.0085610.008040.0077590.0078870.0079740.007879
14.0180.0098200.009180.0088860.0090180.0090870.009011
12.0350.0113310.0106170.0103110.0104470.0104980.010442
10.0530.0135130.0126350.0123160.0124580.0124880.012455
8.0000.0168710.0157990.0154650.0156160.0156220.015615
6.0180.0224100.0209290.0205790.020740.0207230.020743
4.3890.0306360.0286360.0282720.0284430.0284060.028449
393.1524.1420.0067150.0061380.0058910.0060280.0061420.006016
22.0890.0072180.0066360.0063820.0065170.0066220.006505
20.1060.0078060.0072170.0069550.0070890.0071840.007077
18.1240.0085610.007930.007660.0077930.0078760.007782
16.1420.0094840.0088230.0085440.0086770.0087470.008667
14.0890.0108270.010020.0097310.0098640.0099180.009857
12.0350.0125060.0116350.0113350.0114690.0115050.011465
10.0530.0149400.0138310.013520.0136560.0136730.013655
8.0710.0185490.0171210.0167960.0169360.0169330.016939
6.0890.0245080.0225720.0222340.0223790.0223540.022387
4.8850.0306360.0280530.0277070.0278550.0278160.027866
Table 9. MAPEs of experimental and calculated results in the volumetric coefficient prediction.
Table 9. MAPEs of experimental and calculated results in the volumetric coefficient prediction.
MAPESRKPRLKPBWRSGERG-2008Number of Experimental Data Sets
Sample 16.22%9.28%7.89%6.83%7.92%76
Sample 29.28%12.12%10.78%9.74%10.81%77
Sample 34.08%7.11%5.61%4.59%5.67%76
Sample 42.92%5.97%4.49%3.41%4.53%76
Sample 55.16%1.74%3.34%4.49%3.29%74
Overall MAPE5.53%7.24%6.42%5.81%6.44%-
Table 10. MAPEs of experimental and calculated results in the dew point prediction.
Table 10. MAPEs of experimental and calculated results in the dew point prediction.
MAPESRKPRLKPBWRSGERG-2008Number of Experimental Data Sets
Sample 12.81%2.84%1.70%8.20%3.78%43
Sample 21.56%2.44%3.80%6.85%3.02%43
Sample 37.31%8.84%3.17%9.84%7.59%43
Sample 43.73%4.51%4.72%8.68%3.71%44
Sample 56.22%8.52%5.15%8.91%9.65%44
Overall MAPE4.33%5.43%3.71%8.50%5.55%-
Table 11. Summary of the precisions of the five physical property calculation methods.
Table 11. Summary of the precisions of the five physical property calculation methods.
Overall MAPESRKPRLKPBWRSGERG-2008
Gas density3.84%4.92%4.30%4.04%4.31%
Compressibility factor5.17%6.89%6.06%5.46%6.08%
Volumetric coefficient5.53%7.24%6.42%5.81%6.44%
Dew point4.33%5.43%3.71%8.50%5.55%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Chang, M.; Chen, L.; Wang, S.; Fan, S.; Hua, D. Evaluation of Prediction Models for the Physical Properties in Fire-Flooding Exhaust Reinjection Process. Energies 2022, 15, 562. https://doi.org/10.3390/en15020562

AMA Style

Wang Y, Chang M, Chen L, Wang S, Fan S, Hua D. Evaluation of Prediction Models for the Physical Properties in Fire-Flooding Exhaust Reinjection Process. Energies. 2022; 15(2):562. https://doi.org/10.3390/en15020562

Chicago/Turabian Style

Wang, Yong, Mingliang Chang, Long Chen, Shouxi Wang, Shihao Fan, and Dongyang Hua. 2022. "Evaluation of Prediction Models for the Physical Properties in Fire-Flooding Exhaust Reinjection Process" Energies 15, no. 2: 562. https://doi.org/10.3390/en15020562

APA Style

Wang, Y., Chang, M., Chen, L., Wang, S., Fan, S., & Hua, D. (2022). Evaluation of Prediction Models for the Physical Properties in Fire-Flooding Exhaust Reinjection Process. Energies, 15(2), 562. https://doi.org/10.3390/en15020562

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop