Regular Article
An improved method for calculating critical temperatures and critical pressures in natural gas mixtures with up to nC_{11} hydrocarbons
School of Chemical Sciences and Engineering, Yachay Tech University, San Miguel de Urcuquí 100650, Ecuador
^{*} Corresponding author: mricaurte@yachaytech.edu.ec
Received:
24
January
2019
Accepted:
15
April
2019
This study suggests an improvement to the empirical model proposed by Peng (1986, Can. J. Chem. Eng. 64, 827–830) to calculate critical temperatures and critical pressures in natural gas mixtures. It aims to extend its application to natural gas mixtures containing hydrocarbons compounds up to undecane (nC_{11}). This work focuses on establishing new matrices of coefficients A_{ij} by obtaining new binary interactions between heavy compounds and the rest of compounds present in natural gas mixtures. The analysis considered more than 300 natural gas mixtures. Different comparisons were made between calculated critical properties, and referenced critical properties. Mean absolute errors <1.00% for the critical temperatures, and <2.70% for critical pressures were obtained. These low average deviations demonstrate the accuracy of this study, and could be considered as an easytouse engineering tool for estimating critical properties in natural gas mixtures, applicable to lean gas, rich gas, gas condensate, and Natural Gas Liquids (NGL).
© M. Ricaurte et al., published by IFP Energies nouvelles, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Natural gas is a mixture of hydrocarbon compounds, which are in a gaseous state under most ambient conditions. This mixture is mainly composed of paraffinic hydrocarbons molecules, from methane (C_{1}) to undecane (nC_{11}), and some impurities, such as carbon dioxide (CO_{2}), hydrogen sulfide (H_{2}S), nitrogen (N_{2}), and water vapor (H_{2}O) (Gas Processors Suppliers Association, 1998). Given the industrial importance of natural gas as an energy vector in the 21st century, it is necessary to count with thermodynamicbased models to calculate or estimate physical and chemical properties of interest in the oil and gas industry (Mazyan et al., 2016; Nasrifar and Rahmanian, 2018). Among their industrial applications are: incrementing oil recovery factors from reservoirs (Sayed and AlMuntasheri, 2016), production and operational analyses for transport of natural gas in onshore and offshore environments (Campos et al., 2017), design of surface facilities for processing natural gas (Yang et al., 2018), and underground natural gas storage studies (Demirel et al., 2017).
Traditionally, composition of natural gas mixtures is reported up to heptane (nC_{7}), even though in the natural gas industry, some of these mixtures (rich gas) could have a significant content of paraffins larger than nC_{7.} New developments on chromatographic techniques to characterize natural gas mixtures, it is now common practice to determine, with a high degree of accuracy, the hydrocarbon content up to undecane (nC_{11}) (Faramawy et al., 2016; International Organization for Standardization, 2000).
There are rigorous and empirical methods to calculate critical properties of VaporLiquid Equilibrium (VLE) in multicomponent mixtures, applied to natural gas samples (Huang and Guo, 1995). Rigorous methods are based on classical thermodynamic theories, and on the concept of Helmholtz free energy minimization (Bell and Jäger, 2017), which require simultaneous resolution of nonlinear equation systems. Several authors have proposed different strategies and algorithms to solve these equation systems (Heidemann and Khalil, 1980; Michelsen, 1980; Peng and Robinson, 1977). On the other hand, empirical methods are mathematical correlations used to calculate directly critical properties (Ahmed, 2016). Peng (1986) suggested an empirical method to calculate critical temperature (T_{c}) and critical pressure (P_{c}) in natural gas mixtures, which was based on the concept of excess function (Eq. (1)):(1)where ϕ is the desired critical property (T_{c} or P_{c}), x_{i} and x_{j} molar composition of the studied mixture, ϕ_{i} critical properties of individual components, C an arbitrary constant, and A_{ij} fitted parameter for component pair i–j. Despite its easytouse implementation, it is limited to paraffins raging from methane through nC_{7} (Bahadori, 2014), so its application to rich gas mixtures is not recommended (Huang and Guo, 1995).
This paper suggests an enhanced approach to the method proposed by Peng (1986). It takes into account, the effect of heavy paraffins (namely, nC_{7}, nC_{8}, nC_{9}, nC_{10,} and nC_{11}) in calculating critical properties (P_{c} and T_{c}) in natural gas mixtures. This improved approach mainly focuses on optimizing the excess function (Eq. (1)), by proposing new matrices of coefficients A_{ij}, which consider the interaction between heavy paraffins and the rest of components present in natural gas. The significance of this study is that it extends the applicability of the method of Peng (1986), to calculate critical properties of natural gas mixtures containing considerable amounts of heavy hydrocarbons, i.e., rich gas, gas condensates, and Natural Gas Liquids (NGL) (Faramawy et al., 2016; Jin and Lim, 2018).
2 Methods
Figure 1 shows a schematic representation of the methodology developed in this study. It is divided into two sections: (i) excess function optimization, and (ii) validation of calculated critical properties. Samples used during the excess function optimization section (Data Set 1) were different from the ones used for the validation section (Data Set 2).
Fig. 1 Schematic representation of the methodology developed in this study. 
In both sections, the criteria to classify the natural gas is based on the content of (Faramawy et al., 2016). Natural gas containing ≤10 mol% of is considered “lean gas”, and natural gas containing >10 mol% of is defined as “rich gas”.

Excess function optimization section: critical temperatures and critical pressures were obtained from Data Set 1 using commercial chemical process simulators (e.g., SimSci Pro/II, Aspen Hysys, and PVTSim Nova), with Peng and Robinson (1976) Equation of State (EoS), as the suitable thermodynamic model. Details on the accuracy and confidence of thermodynamic packages used by chemical process simulators can be found in a recent book by Foo et al. (2017). The arithmetic average of critical properties obtained from the simulators was taken for each sample (P _{c} & T _{c} Data 1), and used as “critical properties reference”.

The critical properties are fitted using the Peng (1986) correlation by defining new matrices coefficients A _{ ij,} in order to obtain an enhanced excess function. This optimization analysis was performed using the MatLab Optimization Toolbox (MathWorks, 2018).

Validation section: this section was performed using a set of natural gas samples different than those used in the excess function optimization section_{,} thus assuring that the enhanced excess function presented in this study, is not dependent of the dataset used to determine the matrices coefficient A _{ ij }. Critical pressures and temperatures obtained from the enhanced excess function were compared against critical properties calculated from three different sources:

Comparison A: it considers the average of the results obtained from chemical process simulators.

Comparison B: it considers results using the original proposal of Peng (1986).

Comparison C: it takes into consideration some models available in the literature to predict vaporliquid equilibrium properties, including critical points (JustoGarcía et al., 2008; Nasrifar et al., 2005; GuerreroZárate et al., 2017; Jiang and Prausnitz, 2000). These models use experimental critical temperatures and pressures to verify accuracy. The comparison C shows the applicability and robustness of the current study against available models.

2.1 Natural gas mixtures
Table 1 shows number of natural gas mixtures used in this study. It includes both natural gas mixtures to determine matrices coefficients A _{ ij } (199 samples), and natural gas mixtures to validate critical properties calculated (108 samples).
Number of natural gas samples considered in this study.
The total number of natural gas mixtures considered was 307, of which 65.5% were rich gas samples from associated and nonassociated gas fields (Faramawy et al., 2016). To demonstrate the molar composition amplitude of the samples considered in this study, all natural gas samples were classified taking into account: (a) C_{1} content in lean gas and rich gas, (b) content of nonhydrocarbon compounds and (c) the content of heavy compounds (Fig. 2).
Fig. 2 Characterization of natural gas samples: a) C_{1} content in lean gas and rich gas, b) nonhydrocarbon compounds, and c) heavy compounds content. (Note: black symbols represent the Data Set 1 used in excess function optimization and red ones denote the Data Set 2 for validation section). 
2.2 Calculation method used to obtain matrices coefficients A _{ ij }
Matrices coefficients were optimized using the leastsquares method of estimated errors for each critical temperature and pressure analyzed. Equation (1) governing the behavior of the critical constants is the same as that used by Peng (1986).
The new matrices of coefficients A _{ ij } are squares matrices of order 16. The order of a matrix corresponds to its number of compounds in natural gas mixtures: paraffinic hydrocarbons (from C_{1} to nC_{11}), and nonhydrocarbon compounds (N_{2}, CO_{2,} and H_{2}S). The form of matrix A _{ ij } is given by equation (2):(2)The coefficients a _{ ij } of matrix A _{ ij } can be considered as the binary relationship between the different hydrocarbon molecules in a natural gas mixture. Binary relationships between same hydrocarbon molecules are assumed to be 1, therefore, unit values are assigned to the elements in the main diagonal of each matrix (Peng, 1986).
Matrices coefficients A _{ ij } were calculated using the LevenbergMarquardt method (Marquardt, 1963), which was implemented using MatLab Optimization Toolbox. The function to be minimized (Eq. (3)) was established from the quadratic difference between the value calculated using equation (1) (f(x _{ i })) and the average value obtained with simulators (y _{ i }) for a given number (N) of natural gas mixtures:(3)The average absolute percentage error is determined using equation (4), and considering y _{ i } as the correct value:(4)
3 Results
3.1 Determination of coefficients matrices Aij and calculation of critical properties
Tables 2 and 3 show new matrices coefficients A _{ ij } established to calculate critical properties of natural gas mixtures using equation (1). The first subscript (i) corresponds to rows and the second subscript (j) to columns on the new matrices. The constant C was assumed as 2500 K for calculation of critical temperatures, and 250 MPa for calculation of critical pressures.
Coefficients (A _{ ij }) for the calculation of critical temperatures.
Coefficients (A _{ ij }) for the calculation of the critical pressures.
Figure 3 shows results comparison between results obtained from analyzed mixtures to determine coefficients A _{ ij }.
Fig. 3 Critical property calculations: a) critical temperatures, b) critical pressures, and c) Average Absolute Deviation (AAD). 
Figure 3, graphs (a) and (b), shows that calculated critical properties behave similarly to the reference values. These results deviate from the expected values by less than 10% in both cases. Critical temperature results presented less statistical dispersion. Critical pressure results presented greater dispersion, only four out of 199 natural gas mixtures showed results with deviations greater than 10%.
Mean absolute percentage errors were <1.00% for critical temperatures results, and <2.70% for the critical pressure (Fig. 3c). These average errors are acceptable when compared to mean errors reported in the literature (Sloan and Koh, 2008). These small dispersion values confirm the reliability of all selected mixtures, and a high level of statistical confidence to the applied method. With respect to the type of sample, a higher dispersion of the results was found in lean gas samples.
3.2 Validation
Validation was performed using natural gas mixtures other than those used to calculate the coefficients A _{ ij } (Sect. 3.1). Using these natural gas mixtures, confirm the applicability and robustness of the proposed modification to the correlation by Peng (1986) in this study.
3.2.1 Comparison A
Figure 4 shows the comparison of results considering average errors from results obtained using simulators.
Fig. 4 Critical property calculations – Comparison A: a) critical temperatures, b) critical pressures, and c) Average Absolute Deviation (AAD). 
Average errors of 0.90% for the critical temperature and 2.61% for the critical pressure were obtained on all natural gas mixtures. These errors are within the same error range obtained in Section 3.1. Lean gas mixtures have the highest error dispersion and average absolute deviation error.
3.2.2 Comparison B
Comparison B (Tab. 4) is made limiting contents of nC_{7} in natural gas mixtures. Two types of mixtures were considered: (1) natural gas mixtures containing hydrocarbons molecules lighter than or equal to nC_{7} (≤nC_{7}), and (2) natural gas mixtures containing hydrocarbons molecules up to nC_{11}.
Critical property calculations – Comparison B.
Due to the fact that the semiempirical model developed by Peng (1986), considers only mixtures up to nC_{7,} for mixtures type (2) components heavier than nC_{7} and components equal or lighter than nC_{11} were considered as one component, and grouped as the heavy fraction of the mixture, assuming this fraction as nC_{7} plus () (Huang and Guo, 1995).
The improvement made to the semiempirical model of Peng (1986), for the calculation critical properties in natural gas mixtures is noticed by reducing significantly a significant reduction of the average errors of 28.38% for the critical temperature and 66.70% for the critical pressure, applicable to natural gas mixtures with significant amounts of both light and heavy compounds.
3.2.3 Comparison C
Comparison C (Tab. 5) considers other suggested models to calculate critical properties: (i) using simulators with the PengRobinson EoS and (ii) using experimental data reported by different authors (Jiang and Prausnitz, 2000; JustoGarcía et al., 2008; Nasrifar et al., 2005; GuerreroZárate et al., 2017).
Critical property calculations – Comparison C.
In general terms, average errors obtained using the adjustment of the coefficient matrix A _{ ij } performed in this study are comparable to those average errors obtained by other authors, using as reference the values from the simulators and experimental data. This proves the power of coefficient matrices A _{ ij } developed in this study to calculate critical properties of natural gas mixtures, which is applicable to natural gas mixtures containing significant amounts of heavy hydrocarbons components, (up to nC_{11}).
4 Conclusion
This study proposes a modification to the semiempirical model advanced by Peng (1986) for the calculation of the critical properties in natural gas mixtures, aiming to extend its application to rich natural gas streams, i.e., natural gas mixtures containing significant amounts of heavy hydrocarbons (up to nC_{11}). The proposed modification focuses on establishing new matrices of coefficients A _{ ij } used in the calculation. Different comparisons were made between the experimental and reference values of the critical properties, which show the robustness of the improvement of the semiempirical model as an easytouse engineering tool for estimating the critical properties in natural gas mixtures, applicable to both lean and rich gas mixtures. The average errors obtained are less than 1.00% and 2.70% for calculating the critical temperature and pressure, respectively.
Acknowledgments
The authors would like to owe great thanks to Arlán Briceño and JeanPiero Suárez, at Universidad Central de Venezuela for their technical advices.
References
 Ahmed T. (2016) Equations of state and PVT analysis: Applications for improved reservoir modeling, 2nd ed., Gulf Professional Publishing – Elsevier Inc., Cambridge, USA. [Google Scholar]
 Atilhan M., Ejaz S., Zhou J., Cristancho D., Mantilla I., Holste J., Hall K. (2010) Characterization of deepwater natural gas samples. Part 1: 78% methane mixture with heavy components, J. Chem. Eng. Data 55, 4907–4911. [Google Scholar]
 Ávila S., Blanco S., Velasco I., Rauzy E., Otín S. (2002) Thermodynamic properties of synthetic natural gases. 1. Dewpoint curves of synthetic natural gases and their mixtures with water and methanol. Measurement and correlation, Ind. Eng. Chem. Res. 41, 3714–3721. [Google Scholar]
 Ávila S., Benito A., Berro C., Blanco S., Otín S., Velasco I. (2006) Dewpoint curves of natural gas. Measurement and modeling, Ind. Eng. Chem. Res. 45, 5179–5184. [Google Scholar]
 Bahadori A. (2014) Natural gas processing: Technology and engineering design, Gulf Professional Publishing – Elsevier Inc., Waltham, USA. [Google Scholar]
 Bell I., Jäger A. (2017) Calculation of critical points from Helmholtzenergyexplicit mixture models, Fluid Phase Equilib. 433, 159–173. [Google Scholar]
 Campos A., da Silva N., Pereira M., Vasconcelos M. (2017) A review of Brazilian natural gas industry: Challenges and strategies, Renew. Sust. Energ. Rev. 75, 1207–1216. [CrossRef] [Google Scholar]
 Demirel N., Demirel T., Deveci M., Vardar G. (2017) Location selection for underground natural gas storage using Choquet integral, J. Nat. Gas Sci. Eng. 45, 368–379. [Google Scholar]
 Faramawy S., Zaki T., Sakr A. (2016) Natural gas origin, composition, and processing: A review, J. Nat. Gas Sci. Eng. 34, 34–54. [Google Scholar]
 FarzanehGord M., Niazmand A., DeymiDashtebayaz M., Rahbari H. (2015) Effects of natural gas compositions on CNG (compressed natural gas) reciprocating compressors performance, Energy 90, 1, 1152–1162. [CrossRef] [Google Scholar]
 Foo D., Chemmangattuvalappil N., Ng D., Elyas R., Chen C., Elms R., Lee H., Chien I., Chong S., Chong C. (2017) Chemical engineering process simulation, 1st ed., Elsevier, Amsterdam, Netherlands. [Google Scholar]
 Gas Processors Suppliers Association (1998) Engineering data book, 11th ed., GPSA, Tulsa, Okla. [Google Scholar]
 GuerreroZárate D., EstradaBaltazar A., IglesiasSilva A. (2017) Calculation of critical points for natural gas mixtures with the GERG2008 equation of state, Fluid Phase Equilib. 437, 69–82. [Google Scholar]
 Heidemann R., Khalil M. (1980) The calculation of critical points, AlChE J. 26, 2, 769–779. [CrossRef] [Google Scholar]
 Huang Q., Guo T. (1995) Prediction of the critical points of natural gas mixtures by rigorous and semiempirical methods, J. Pet. Sci. Eng. 13, 3–4, 233–245. [Google Scholar]
 International Organization for Standardization (2000) Natural gas – Determination of composition with defined uncertainty by gas chromatography – Part 4: Determination of nitrogen, carbon dioxide and C_{1} to C_{5} and C_{6}+ hydrocarbons for a laboratory and online measuring system using two columns (ISO 69744:2000). [Google Scholar]
 Jaeschke M., Humpherys A. (1991) The GERG databank of high accuracy compressibility factor measurements, Fortschr.Ber. VDI Reihe 6, 251, 58. [Google Scholar]
 Jiang J., Prausnitz J. (2000) Critical temperatures and pressures for hydrocarbon mixtures from an equation of state with renormalizationgroup theory corrections, Fluid Phase Equilib. 169, 127–147. [Google Scholar]
 Jin C., Lim Y. (2018) Economic evaluation of NGL recovery process schemes for lean feed compositions, Chem. Eng. Res. Des. 129, 297–305. [Google Scholar]
 JustoGarcía D., GarcíaSánchez F., ÁguilaHernández J., EustaquioRincón R. (2008) Application of the simulated annealing technique to the calculation of critical points of multicomponent mixtures with cubic equations of state, Fluid Phase Equilib. 264, 164–173. [Google Scholar]
 Kunz O., Wagner W. (2012) The GERG2008 widerange equation of state for natural gases and other mixtures: An expansion of GERG2004, J. Chem. Eng. Data 57, 3032–3091. [Google Scholar]
 Louli V., Pappa G., Boukouvalas C., Skouras S., Solbraa E., Christensen K., Voutsas E. (2012) Measurement and prediction of dew point curves of natural gas mixtures, Fluid Phase Equilib. 334, 2012, 1–9. [Google Scholar]
 Marquardt D. (1963) An algorithm for leastsquares estimation of nonlinear parameters, J. Soc. Ind. Appl. Math. 11, 2, 431–441. [CrossRef] [MathSciNet] [Google Scholar]
 MathWorks (2018) Optimization toolbox™ user’s guide. Version 8.2 (R2018b), The MathWorks Inc, Natick, USA. [Google Scholar]
 Mazyan W., Ahmadi A., Ahmed H., Hoorfar M. (2016) Market and technology assessment of natural gas processing: A review, J. Nat. Gas Sci. Eng. 30, 487–514. [Google Scholar]
 Michelsen M. (1980) Calculation of phase envelopes and critical points for multicomponent mixtures, Fluid Phase Equilib. 4, 1–2, 1–10. [Google Scholar]
 Nasir Q., Sabil K., Lau K. (2015) Measurement of isothermal (vapor + liquid) equilibria, (VLE) for binary (CH_{4} + CO_{2}) from T = (240.35 to 293.15) K and CO_{2} rich synthetic natural gas systems from T = (248.15 to 279.15) K, J. Nat. Gas Sci. Eng. 27, 1, 158–167. [Google Scholar]
 Nasrifar K., Bolland O., Moshfeghian M. (2005) Predicting natural gas dew points from 15 equations of state, Energy Fuels 19, 561–572. [Google Scholar]
 Nasrifar K., Rahmanian N. (2018) Equations of state with group contribution binary interaction parameters for calculation of twophase envelopes for synthetic and real natural gas mixtures with heavy fractions, Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles 73, 7. [CrossRef] [Google Scholar]
 Peng D. (1986) An empirical method for calculating vapor–liquid critical points of multicomponent mixtures, Can. J. Chem. Eng. 64, 827–830. [Google Scholar]
 Peng D., Robinson D. (1976) A new twoconstant equation of state, Ind. Eng. Chem. Fundam. 15, 1, 59–64. [CrossRef] [Google Scholar]
 Peng D., Robinson D. (1977) A rigorous method for predicting the critical properties of multicomponent systems from an equation of state, AlChE J. 23, 2, 137–144. [CrossRef] [Google Scholar]
 Sayed M., AlMuntasheri G. (2016) Mitigation of the effects of condensate banking: A critical review, SPE Prod. Oper. 31, 2, 85–102. [Google Scholar]
 Sloan E., Koh C. (2008) Clathrate hydrates of natural gases, 3rd ed., CRC Press, Taylor & Francis, New York, NY. [Google Scholar]
 Yang Y., Li J., Wang S., Wen C. (2018) Gasliquid twophase flow behavior in terraininclined pipelines for gathering transport system of wet natural gas, Int. J. Press. Vessels Piping 162, 52–58. [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Schematic representation of the methodology developed in this study. 

In the text 
Fig. 2 Characterization of natural gas samples: a) C_{1} content in lean gas and rich gas, b) nonhydrocarbon compounds, and c) heavy compounds content. (Note: black symbols represent the Data Set 1 used in excess function optimization and red ones denote the Data Set 2 for validation section). 

In the text 
Fig. 3 Critical property calculations: a) critical temperatures, b) critical pressures, and c) Average Absolute Deviation (AAD). 

In the text 
Fig. 4 Critical property calculations – Comparison A: a) critical temperatures, b) critical pressures, and c) Average Absolute Deviation (AAD). 

In the text 