Skip to content
Brent$117.6Gas TTF$1,192USD453.87 ₸
MunaiHub

Equations: Gas properties

The equations of the “Gas properties” lessons in “Reservoir engineering”, with their conditions and units, and what each symbol means in them.

Equation of state of an ideal gas
pV=nRTpV = nRT

SI (any consistent system of units)

Holds when: The molecules have no volume of their own and do not attract one another. This is a good approximation at low pressure and at a temperature well above the critical one (the highest temperature at which pressure can still turn the gas into a liquid). Pressure and temperature are absolute.

Ideal and real gas: the gas deviation factor (Z) →
Equation of state of a real gas
pV=ZnRT,Z=VrealVidealpV = ZnRT,\qquad Z = \frac{V_{\text{real}}}{V_{\text{ideal}}}

any consistent system of units

Holds when: ZZ is not a constant. It depends on the composition of the gas, the pressure and the temperature, so it must be found for each state. The gas is single-phase (no liquid drops out in the reservoir or in the pipe).

Ideal and real gas: the gas deviation factor (Z) →
Reduced temperature and pressure
Tr=TTc,pr=ppcT_r = \frac{T}{T_c},\qquad p_r = \frac{p}{p_c}

dimensionless ratios; temperature only in kelvins or degrees Rankine

Holds when: The law works best for simple non-polar molecules (those with evenly spread electric charge: methane, nitrogen, argon). For large, asymmetric and polar molecules (heavy hydrocarbons, hydrogen sulfide, water) the deviations are larger.

Ideal and real gas: the gas deviation factor (Z) →
Pseudo-critical properties of a mixture (Kay's rule)
Tpc=∑iyi Tc,i,ppc=∑iyi pc,iT_{pc} = \sum_i y_i\,T_{c,i},\qquad p_{pc} = \sum_i y_i\,p_{c,i}

K and MPa (or °R and psia)

Holds when: The composition of the gas is fully known, ∑yi=1\sum y_i = 1. For the C7+ fraction (the heavy end of the mixture, treated as one component) the critical properties are estimated from correlations that use its molar mass and density. There are few non-hydrocarbon components. Hydrogen sulfide and carbon dioxide need a correction (below).

Ideal and real gas: the gas deviation factor (Z) →
Pseudo-reduced temperature and pressure
Tpr=TTpc,ppr=pppcT_{pr} = \frac{T}{T_{pc}},\qquad p_{pr} = \frac{p}{p_{pc}}

dimensionless ratios

Holds when: TT and pp are absolute and in the same units as TpcT_{pc} and ppcp_{pc}.

Ideal and real gas: the gas deviation factor (Z) →
Pseudo-critical properties of natural gas from specific gravity (Standing's approximation)
Tpc=168+325 γg−12.5 γg2,ppc=677+15 γg−37.5 γg2T_{pc} = 168 + 325\,\gamma_g - 12.5\,\gamma_g^2,\qquad p_{pc} = 677 + 15\,\gamma_g - 37.5\,\gamma_g^2

oilfield: TpcT_{pc} in degrees Rankine, ppcp_{pc} in psia

Holds when: Natural gases that are dry and lean (with little heavy hydrocarbon) and have no more than a few per cent of nitrogen, carbon dioxide and hydrogen sulfide. This is an approximation of a graphical correlation. For gas-condensate mixtures (gases that drop liquid when the pressure falls) another curve is used (advanced section). For wet gases (gases rich in heavier hydrocarbons), later correlations are more accurate. When there is an appreciable amount of sour gases, we account for them separately, with a correction, as in the next section.

Ideal and real gas: the gas deviation factor (Z) →
Wichert–Aziz correction to the pseudo-critical properties
Tpc′=Tpc−ε,ppc′=ppc Tpc′Tpc+B (1−B) εT'_{pc} = T_{pc} - \varepsilon,\qquad p'_{pc} = \frac{p_{pc}\,T'_{pc}}{T_{pc} + B\,(1 - B)\,\varepsilon}

temperatures and ε\varepsilon in the same units (K or °R), pressure in any

Holds when: The correction was built from experiments with gases that contained up to about 54 mol% (mole per cent) of carbon dioxide and up to 74 mol% of hydrogen sulfide. If the gas contains a lot of nitrogen or heavy hydrocarbons, an accurate calculation needs a laboratory analysis or an equation of state.

Ideal and real gas: the gas deviation factor (Z) →
Wichert–Aziz correction term
ε=120 (A0.9−A1.6)+15 (B0.5−B4)\varepsilon = 120\,\left(A^{0.9} - A^{1.6}\right) + 15\,\left(B^{0.5} - B^{4}\right)

ε\varepsilon in degrees Rankine; for kelvins, divide the result by 1.8

Holds when: The fractions are mole fractions, written as numbers between 0 and 1, not as percentages.

Ideal and real gas: the gas deviation factor (Z) →
Gas deviation factor of a van der Waals gas
Z=VmVm−b−aR T VmZ = \frac{V_m}{V_m - b} - \frac{a}{R\,T\,V_m}

SI or any consistent units

Holds when: This is a qualitative model. It is not accurate enough for natural gas in numbers, so in practice more advanced equations of state are used (Cubic equations of state).

Ideal and real gas: the gas deviation factor (Z) →
Dranchuk–Abou-Kassem equation
Z=1+(A1+A2Tpr+A3Tpr3+A4Tpr4+A5Tpr5)ρr+(A6+A7Tpr+A8Tpr2)ρr2−A9(A7Tpr+A8Tpr2)ρr5+A10(1+A11ρr2)ρr2Tpr3 e−A11ρr2Z = 1 + \left(A_1 + \frac{A_2}{T_{pr}} + \frac{A_3}{T_{pr}^3} + \frac{A_4}{T_{pr}^4} + \frac{A_5}{T_{pr}^5}\right)\rho_r + \left(A_6 + \frac{A_7}{T_{pr}} + \frac{A_8}{T_{pr}^2}\right)\rho_r^2 - A_9\left(\frac{A_7}{T_{pr}} + \frac{A_8}{T_{pr}^2}\right)\rho_r^5 + A_{10}\left(1 + A_{11}\rho_r^2\right)\frac{\rho_r^2}{T_{pr}^3}\,\e^{-A_{11}\rho_r^2}

dimensionless

Holds when: It applies for 1.0<Tpr≤3.01.0 < T_{pr} \le 3.0 and 0.2≤ppr<300.2 \le p_{pr} < 30. For TprT_{pr} close to 1 and ppr>1p_{pr} > 1 the accuracy is noticeably worse. It approximates the same generalised chart, so it has the same limits: sour gases need corrected pseudo-critical properties.

Ideal and real gas: the gas deviation factor (Z) →
Pseudo-critical properties of a gas-condensate mixture from specific gravity (Standing's approximation)
Tpc=187+330 γg−71.5 γg2,ppc=706−51.7 γg−11.1 γg2T_{pc} = 187 + 330\,\gamma_g - 71.5\,\gamma_g^2,\qquad p_{pc} = 706 - 51.7\,\gamma_g - 11.1\,\gamma_g^2

oilfield: °R and psia (divide by 1.8 for K, multiply by 0.006895 for MPa)

Holds when: The mixture in the reservoir is single-phase (the pressure is above the dew point pressure, the pressure at which the first drops of liquid form). There are few non-hydrocarbon components. For rich mixtures it is more accurate to calculate from the composition, with C7+ pseudo-components (lumped groups of heavy components), or to use an equation of state.

Ideal and real gas: the gas deviation factor (Z) →
Gas formation volume factor
Bg=VresVsc=pscp⋅Z TTscB_g = \frac{V_{\text{res}}}{V_{\text{sc}}} = \frac{p_{\text{sc}}}{p}\cdot\frac{Z\,T}{T_{\text{sc}}}

absolute pressures and temperatures in consistent units

Holds when: The gas is single-phase in the reservoir and remains a gas all the way to standard conditions. For wet gas (gas rich in heavier hydrocarbons) and gas condensate (gas that drops a light liquid), part of the gas turns into liquid at the surface, and the conversion is made from the composition (Types of reservoir fluids). Zsc=1Z_{\text{sc}} = 1. For natural gas at standard conditions, ZZ differs from 1 by only fractions of a per cent.

Gas formation volume factor, compressibility and viscosity →
Gas density in the reservoir
ρg=ρg,scBg,ρg,sc=1.204 γg\rho_g = \frac{\rho_{g,\text{sc}}}{B_g},\qquad \rho_{g,\text{sc}} = 1.204\,\gamma_g

kg/m³; BgB_g is for the same standard conditions as ρg,sc\rho_{g,\text{sc}}

Holds when: The gas is single-phase both at the surface and in the reservoir (all its mass remains gas).

Gas formation volume factor, compressibility and viscosity →
Isothermal compressibility of a gas
cg=−1V(∂V∂p)T=1p−1Z(∂Z∂p)Tc_g = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T = \frac{1}{p} - \frac{1}{Z}\left(\frac{\partial Z}{\partial p}\right)_T

any consistent units: MPa and MPa⁻¹ or psi and psi⁻¹

Holds when: The temperature is constant: we treat the reservoir as isothermal. The gas is single-phase. The second form comes from substituting V=ZnRT/pV = ZnRT/p (the derivation is in the advanced section).

Gas formation volume factor, compressibility and viscosity →
Pseudo-reduced compressibility of a gas
cpr=cg ppc=1ppr−1Z(∂Z∂ppr)Tprc_{pr} = c_g\,p_{pc} = \frac{1}{p_{pr}} - \frac{1}{Z}\left(\frac{\partial Z}{\partial p_{pr}}\right)_{T_{pr}}

dimensionless

Holds when: We take the slope from the same generalised chart (or its approximation) that gives ZZ.

Gas formation volume factor, compressibility and viscosity →
Viscosity of natural gas (Lee–Gonzalez–Eakin correlation)
μg=10−4 K exp⁡(Xρg Y),K=(9.4+0.02 M) T1.5209+19 M+T,X=3.5+986T+0.01 M,Y=2.4−0.2 X\begin{gathered}\mu_g = 10^{-4}\,K\,\exp\left(X\rho_g^{\,Y}\right),\qquad K = \frac{(9.4 + 0.02\,M)\,T^{1.5}}{209 + 19\,M + T},\\ X = 3.5 + \frac{986}{T} + 0.01\,M,\qquad Y = 2.4 - 0.2\,X\end{gathered}

empirical: TT in °R, ρg\rho_g in g/cm³, μg\mu_g in mPa·s

Holds when: It was built from measurements on natural gases at temperatures from about 40 to 170 °C and pressures up to 55 MPa. The error is a few per cent. For sour gases, gases with a lot of nitrogen, and gas condensates the accuracy is lower. Calculate the gas density with a ZZ that includes the correction for the sour components.

Gas formation volume factor, compressibility and viscosity →
Real gas pseudo-pressure
m(p)=2∫p0pp′μg(p′) Z(p′)  dp′m(p) = 2\int_{p_0}^{p}\frac{p'}{\mu_g(p')\,Z(p')}\,\dd p'

consistent; in oilfield units psi²/cP

Holds when: The reservoir is isothermal, and the gas is single-phase.

Gas formation volume factor, compressibility and viscosity →
Viscosity of a gas mixture at low pressure
μmix=∑iμi yiMi∑iyiMi\mu_{\text{mix}} = \frac{\sum_i \mu_i\,y_i\sqrt{M_i}}{\sum_i y_i\sqrt{M_i}}

viscosities in the same units

Holds when: The pressure is low (of the order of atmospheric) and the gas is close to ideal. For reservoir pressures, multiply the result by the ratio μ/μ1\mu/\mu_1 from a chart, or recalculate it from the density.

Gas formation volume factor, compressibility and viscosity →
Van der Waals equation
p=RTVm−b−aVm2p = \frac{RT}{V_m - b} - \frac{a}{V_m^2}

metric: MPa, K, m³/kmol (any consistent system)

Holds when: aa and bb are constants for a substance and do not depend on temperature. This is the main simplification of the equation. Pressure and temperature are absolute. If a=b=0a = b = 0, we get the ideal gas equation pVm=RTpV_m = RT.

Cubic equations of state →
Parameters of the van der Waals equation from the critical constants
a=2764 R2Tc2pc,b=18 RTcpc,Zc=pcVcRTc=38a = \frac{27}{64}\,\frac{R^2T_c^2}{p_c},\qquad b = \frac{1}{8}\,\frac{RT_c}{p_c},\qquad Z_c = \frac{p_cV_c}{RT_c} = \frac{3}{8}

any consistent units; with R=0.008314R = 0.008314 MPa·m³/(kmol·K), aa comes out in MPa·m⁶/kmol² and bb in m³/kmol

Holds when: At the critical point the isotherm has a horizontal inflection (it flattens and changes its bend): ∂p/∂Vm=∂2p/∂Vm2=0\partial p/\partial V_m = \partial^2 p/\partial V_m^2 = 0 at T=TcT = T_c. The parameters are chosen so that the equation passes exactly through the critical temperature and pressure. The critical volume then comes out as the equation gives it.

Cubic equations of state →
General form of two-parameter cubic equations
p=RTVm−b−ac α(T)(Vm+δ1b)(Vm+δ2b),ac=ΩaR2Tc2pc,b=ΩbRTcpcp = \frac{RT}{V_m - b} - \frac{a_c\,\alpha(T)}{(V_m + \delta_1 b)(V_m + \delta_2 b)},\qquad a_c = \Omega_a\frac{R^2T_c^2}{p_c},\quad b = \Omega_b\frac{RT_c}{p_c}

any consistent system; Ωa\Omega_a, Ωb\Omega_b, δ1\delta_1, δ2\delta_2 are dimensionless

Holds when: A substance is described by three constants: TcT_c, pcp_c and (for SRK and PR) the acentric factor ω\omega. The equation passes exactly through the critical temperature and pressure of each substance, but not through its critical volume. bb does not depend on temperature. All the dependence on temperature is in α(T)\alpha(T).

Cubic equations of state →
Peng–Robinson equation
p=RTVm−b−ac α(T)Vm2+2bVm−b2,ac=0.45724 R2Tc2pc,b=0.07780 RTcpcp = \frac{RT}{V_m - b} - \frac{a_c\,\alpha(T)}{V_m^2 + 2bV_m - b^2},\qquad a_c = 0.45724\,\frac{R^2T_c^2}{p_c},\quad b = 0.07780\,\frac{RT_c}{p_c}

metric: MPa, K, m³/kmol (the constants 0.45724 and 0.07780 are dimensionless and suit any system)

Holds when: The denominator is Vm2+2bVm−b2=(Vm+2.414 b)(Vm−0.414 b)V_m^2 + 2bV_m - b^2 = (V_m + 2.414\,b)(V_m - 0.414\,b), that is, δ1=1+2\delta_1 = 1 + \sqrt2 and δ2=1−2\delta_2 = 1 - \sqrt2. For a mixture, take aa and bb from the mixing rules (section on mixtures).

Cubic equations of state →
Temperature function of the attraction in the SRK and PR equations
α(T)=[1+m(1−T/Tc)]2\alpha(T) = \left[1 + m\left(1 - \sqrt{T/T_c}\right)\right]^2

dimensionless; temperatures absolute

Holds when: The dependence m(ω)m(\omega) was fitted to the vapour pressures of hydrocarbons, so each equation has its own. For heavy substances with ω>0.49\omega > 0.49, the PR equation uses a refined dependence (advanced section).

Cubic equations of state →
Acentric factor
ω=−lg⁡pspc∣T=0.7 Tc−1\omega = -\lg\left.\frac{p_s}{p_c}\right|_{T = 0.7\,T_c} - 1

dimensionless; pressures absolute, in the same units

Holds when: It is defined from the vapour pressure curve of the liquid. For carbon dioxide the temperature 0.7 Tc0.7\,T_c is below the triple point (where solid, liquid and gas coexist), so the curve is extended beyond it. For C7+ fractions and pseudo-components, ω\omega is estimated from correlations.

Cubic equations of state →
Peng–Robinson equation in terms of the Z-factor
Z3−(1−B) Z2+(A−3B2−2B) Z−(AB−B2−B3)=0Z^3 - (1 - B)\,Z^2 + (A - 3B^2 - 2B)\,Z - (AB - B^2 - B^3) = 0

dimensionless form; pp and TT absolute

Holds when: AA and BB are calculated at the same absolute pressure and temperature for which ZZ is sought. For a mixture, aa and bb are calculated from the composition of the phase whose properties we want.

Cubic equations of state →
Van der Waals mixing rules
a=∑i∑jxixj (1−kij)aiaj,b=∑ixibia = \sum_i\sum_j x_ix_j\,(1 - k_{ij})\sqrt{a_ia_j},\qquad b = \sum_i x_ib_i

as for the parameters of the components

Holds when: The mixture is treated as one “average” substance (the van der Waals mixing rules). The kijk_{ij} do not depend on pressure, and in the simplest version not on temperature either. Use the composition of the phase whose properties you are calculating. The gas and the liquid of the same mixture have different aa and bb.

Cubic equations of state →
Fugacity and the condition of phase equilibrium
fi=φi xi p,fi L=fi V  for each componentf_i = \varphi_i\,x_i\,p,\qquad f_i^{\,L} = f_i^{\,V}\ \ \text{for each component}

pressure in any units; φi\varphi_i is dimensionless

Holds when: The phases are at the same pressure and temperature. The fugacities of each component are equal in all phases. This is the condition of equilibrium, and it is equivalent to equal chemical potentials (the quantities that drive a component to move from one phase to another). For an ideal gas φi=1\varphi_i = 1, and the fugacity equals the partial pressure yipy_ip.

Cubic equations of state →
Fugacity coefficient of a mixture component from the Peng–Robinson equation
ln⁡φi=bib(Z−1)−ln⁡(Z−B)−A22 B(2∑jxjaija−bib)ln⁡Z+(1+2)BZ+(1−2)B\ln\varphi_i = \frac{b_i}{b}(Z - 1) - \ln(Z - B) - \frac{A}{2\sqrt2\,B}\left(\frac{2\sum_j x_ja_{ij}}{a} - \frac{b_i}{b}\right)\ln\frac{Z + (1 + \sqrt2)B}{Z + (1 - \sqrt2)B}

dimensionless; natural logarithms

Holds when: Van der Waals mixing rules with the coefficients kijk_{ij}. ZZ, AA and BB are calculated for the composition and the root of this particular phase.

Cubic equations of state →
Volume translation (Péneloux correction)
Vm=VmEOS−c,c=∑ixiciV_m = V_m^{\text{EOS}} - c,\qquad c = \sum_i x_ic_i

m³/kmol

Holds when: The cic_i do not depend on pressure, and in the simplest version not on temperature either. The cic_i are fitted to the measured density of a liquid (the pure component, stock-tank oil, which is oil at surface conditions, or reservoir oil) or estimated from correlations. For c>0c > 0 the volume decreases and the density rises.

Cubic equations of state →
Depressed cubic equation
t3+P t+Q=0,P=c1−c223,Q=2c2327−c2c13+c0t^3 + P\,t + Q = 0,\qquad P = c_1 - \frac{c_2^2}{3},\quad Q = \frac{2c_2^3}{27} - \frac{c_2c_1}{3} + c_0

dimensionless

Holds when: The coefficients are real numbers. The sign of D=Q2/4+P3/27D = Q^2/4 + P^3/27 decides how many real roots there are.

Cubic equations of state →
Fugacity coefficient from an equation of state
ln⁡φi=∫V∞[1RT(∂p∂ni)T,V,nj≠i−1V] dV−ln⁡Z\ln\varphi_i = \int_V^{\infty}\left[\frac{1}{RT}\left(\frac{\partial p}{\partial n_i}\right)_{T,V,n_{j\ne i}} - \frac{1}{V}\right]\dd V - \ln Z

dimensionless; RR, TT, VV, pp in consistent units

Holds when: We integrate at constant temperature and constant amounts of the other components. At infinite volume the phase is an ideal gas.

Cubic equations of state →
Refined dependence m(ω) of the Peng–Robinson equation for heavy components
m=0.379642+1.48503 ω−0.164423 ω2+0.016666 ω3m = 0.379642 + 1.48503\,\omega - 0.164423\,\omega^2 + 0.016666\,\omega^3

dimensionless

Holds when: The authors of the equation proposed it in 1978 for substances with ω>0.49\omega > 0.49. For smaller ω\omega the formula of 1976 is used.

Cubic equations of state →

Symbols

SymbolMeaning and context
(∂Z/∂p)T\left(\partial Z/\partial p\right)_T
aa
AA
A1…A11A_1 \ldots A_{11}
aca_c
ac, ba_c,\ b
ai, bia_i,\ b_i
  • parameters of the component at temperature TT: ai=ac,i αi(T)a_i = a_{c,i}\,\alpha_i(T), MPa·m⁶/kmol², m³/kmolCubic equations of state
aija_{ij}
  • (1−kij)aiaj(1 - k_{ij})\sqrt{a_ia_j} — the attraction of a pair of components, MPa·m⁶/kmol²Cubic equations of state
a, ba,\ b
  • parameters of the substance at temperature TT (for a mixture, from the mixing rules), MPa·m⁶/kmol², m³/kmolCubic equations of state
a, b, A, Ba,\ b,\ A,\ B
  • parameters of the phase from the mixing rules and their dimensionless forms, MPa·m⁶/kmol², m³/kmol; AA and BB are dimensionlessCubic equations of state
bb
BB
BgB_g
bib_i
cc
c2, c1, c0c_2,\ c_1,\ c_0
cgc_g
cic_i
cprc_{pr}
fif_i
kijk_{ij}
  • binary interaction coefficient; kii=0k_{ii} = 0, kij=kjik_{ij} = k_{ji}, —Cubic equations of state
L, VL,\ V
lg⁡\lg
mm
MM
m(p)m(p)
nn
nin_i
pp
p0p_0
ppcp_{pc}
ppc, pprp_{pc},\ p_{pr}
psp_s
psc, Tscp_{\text{sc}},\ T_{\text{sc}}
p, Tp,\ T
P, QP,\ Q
p, V, n, R, Tp,\ V,\ n,\ R,\ T
p, T, Vm, Rp,\ T,\ V_m,\ R
p, T, Vm, R, bp,\ T,\ V_m,\ R,\ b
p(T,V,n1,n2,…)p(T, V, n_1, n_2, \ldots)
  • the equation of state written for nn kilomoles of the mixture: it contains n2a=∑i∑jninjaijn^2a = \sum_i\sum_j n_in_ja_{ij} and nb=∑inibinb = \sum_i n_ib_i, MPaCubic equations of state
RR
tt
TT
Tc, pcT_c,\ p_c
Tc,i, pc,iT_{c,i},\ p_{c,i}
TpcT_{pc}
Tpc, ppcT_{pc},\ p_{pc}
Tpr, pprT_{pr},\ p_{pr}
Tr, prT_r,\ p_r
T, pT,\ p
Tpc′, ppc′T'_{pc},\ p'_{pc}
T/TcT/T_c
VV
VcV_c
VidealV_{\text{ideal}}
VmV_m
VmEOSV_m^{\text{EOS}}
VrealV_{\text{real}}
xix_i
yiy_i
yi, Miy_i,\ M_i
ZZ
ZcZ_c
α(T)\alpha(T)
γg\gamma_g
δ1, δ2\delta_1,\ \delta_2
ε\varepsilon
μg\mu_g
μg, Z\mu_g,\ Z
μi\mu_i
ρg\rho_g
ρg,sc\rho_{g,\text{sc}}
ρr\rho_r
φi\varphi_i
  • fugacity coefficient of the component in the phase (not to be confused with porosity); it tends to 1 as p→0p \to 0, —Cubic equations of state
ω\omega
Ωa, Ωb\Omega_a,\ \Omega_b