How do you calculate gas density at operating conditions?

Real gas density is calculated as ρ = (P × MW) / (Z × R × T), where P is absolute pressure, MW is molecular weight, Z is the compressibility factor, R is the gas constant (10.7316 psia·ft³/lbmol·°R), and T is absolute temperature in °R.

When do you need to use a Z-factor instead of the ideal gas law?

The ideal gas law is accurate (less than 1% error) when pressure is below 100 psia and temperature exceeds 100°F. Above 500 psia, Z-factor correction is required as errors exceed 5%. Near the hydrocarbon dew point, real gas equations must always be used.

What is the typical Z-factor range for natural gas transmission pipelines?

For transmission pipelines operating at 800–1000 psia and 50–80°F, the typical Z-factor range is 0.85–0.90. At these conditions, real gas is about 10–15% denser than the ideal gas prediction.

How do you estimate pseudo-critical properties when gas composition is unknown?

The Sutton correlation (1985) estimates pseudo-critical properties from specific gravity: T_pc = 169.2 + 349.5×SG − 74.0×SG² (°R) and P_pc = 756.8 − 131.0×SG − 3.6×SG² (psia). It is valid for SG = 0.57 to 1.68 with ±1% accuracy for most sweet natural gases.

What is the specific heat ratio (k) and why does it matter for compressors?

The specific heat ratio k = Cp/Cv (also called the isentropic exponent) governs compression work, discharge-temperature rise, efficiency conversion, and the speed of sound c = √(kZRT/MW). Natural gas is typically k = 1.26–1.32 (methane 1.31, N2 1.40, CO2 1.29, rich gas 1.15–1.25). Using ideal-gas k at high pressure introduces significant errors — calculate real-gas Cp and Cv from an equation of state.

How are pseudo-critical properties corrected for sour gas (H2S, CO2)?

Use the Wichert-Aziz correction: ε = 120(A^0.9 − A^1.6) + 15(B^0.5 − B^4), where A = mol fraction CO2 + H2S and B = mol fraction H2S. Then T_pc' = T_pc − ε and P_pc' = P_pc × T_pc'/[T_pc + B(1−B)ε]. The correction is applied before computing reduced properties and the Z-factor.

Gas Density Fundamentals | Engineering Guide

1. Overview & Applications

Gas density is the mass per unit volume of a gas at specified pressure and temperature. Accurate density calculations are fundamental to:

Pipeline hydraulics

Flow calculations

Reynolds number, pressure drop, erosional velocity all depend on ρ.

Metering

Mass flow rate

Orifice, turbine, and ultrasonic meters require density for mass flow.

Equipment sizing

Compressors, separators

Compressor power and separator sizing use gas density.

Custody transfer

Contract volumes

Standard volume calculations require accurate density at base conditions.

Key Concepts

Density (ρ): Mass per unit volume, typically lb/ft³ or kg/m³
Specific gravity (SG): Gas density relative to air density (dimensionless)
Molecular weight (MW): Average MW of gas mixture, lb/lbmol or g/mol
Compressibility (Z): Deviation from ideal gas behavior (Z = 1 for ideal gas)

Why density matters: A 5% error in density translates directly to 5% error in mass flow rate, which impacts custody transfer revenue, pipeline capacity calculations, and equipment sizing.

2. Ideal Gas Law

The ideal gas law assumes no intermolecular forces and that gas molecules occupy negligible volume. Valid at low pressures and high temperatures relative to the critical point.

Fundamental Equation

Ideal Gas Law: PV = nRT Where: P = Absolute pressure (psia or kPa abs) V = Volume (ft³ or m³) n = Number of moles (lbmol or kmol) R = Universal gas constant = 10.7316 psia·ft³/lbmol·°R (US) = 8.314 kJ/kmol·K (SI) T = Absolute temperature (°R or K)

Ideal gas molecular model showing point mass assumption with negligible volume, no intermolecular forces, elastic collisions only, and random molecular motion with PV=nRT equation and validity conditions listed — Ideal gas assumptions: molecules as point masses with no intermolecular forces, elastic collisions, and random motion. Valid at P < 100 psia when gas is far from condensation.

Density Form

Gas Density (Ideal): ρ = (P × MW) / (R × T) Where: ρ = Gas density (lb/ft³ or kg/m³) MW = Molecular weight (lb/lbmol or kg/kmol) P = Absolute pressure (psia or kPa abs) R = 10.7316 psia·ft³/lbmol·°R or 8.314 kJ/kmol·K T = Absolute temperature °R = °F + 459.67 K = °C + 273.15

Specific Gravity Form

Using Specific Gravity: SG = MW_gas / MW_air = MW_gas / 28.9647 Standard air density (14.73 psia, 60°F): ρ_air = 0.0765 lb/ft³ Quick estimation formula: ρ = 2.70 × SG × P / T Where: P = Absolute pressure (psia) T = Absolute temperature (°R) ρ = Gas density (lb/ft³) Note: This is an ideal gas approximation.

When Ideal Gas Law Applies

Condition	Ideal Gas Accuracy	Recommendation
P < 100 psia, T > 100°F	< 1% error	Use ideal gas law
100 < P < 500 psia	1–5% error	Consider Z-factor correction
P > 500 psia	> 5% error	Must use Z-factor or EOS
Near dew point	Highly inaccurate	Must use real gas equations

Example Calculation

Calculate density of methane (MW = 16.04) at 100 psia and 80°F:

T = 80 + 459.67 = 539.67 °R P = 100 psia ρ = (100 × 16.04) / (10.73 × 539.67) ρ = 1604 / 5790.5 ρ = 0.277 lb/ft³ Alternatively using SG: SG = 16.04 / 28.97 = 0.554 ρ = 2.70 × 0.554 × 100 / 539.67 ρ = 0.277 lb/ft³ ✓

3. Real Gas & Compressibility Factor

Real gases deviate from ideal behavior due to intermolecular attractive/repulsive forces and finite molecular volume. The compressibility factor Z quantifies this deviation.

Real Gas Equation

Real Gas Law: PV = Z n R T Where Z = compressibility factor (dimensionless) For density: ρ = (P × MW) / (Z × R × T) Or: ρ = ρ_ideal / Z Since Z < 1 typically, real gas density is higher than ideal prediction.

Standing-Katz Chart Method

The Standing-Katz correlation (1942) is the industry-standard graphical method for determining Z-factor from reduced pressure and temperature:

Reduced Properties: P_r = P / P_pc (pseudo-reduced pressure) T_r = T / T_pc (pseudo-reduced temperature) Where: P_pc = pseudo-critical pressure (psia) T_pc = pseudo-critical temperature (°R) For natural gas mixtures (Kay's mixing rule): P_pc = Σ(y_i × P_ci) T_pc = Σ(y_i × T_ci) Z = f(P_r, T_r) from Standing-Katz chart

Standing-Katz compressibility factor chart showing Z versus pseudo-reduced pressure P_r from 0 to 8 with curves for pseudo-reduced temperatures T_r from 1.05 to 3.0, ideal gas line at Z=1.0, and annotations showing maximum deviation from ideal gas and pipeline conditions — Standing-Katz compressibility chart: Z-factor vs. pseudo-reduced pressure at various T_r values. Lower T_r curves show deeper dips (more deviation from ideal gas behavior).

Sutton Correlation for P_pc and T_pc

When full gas composition is unavailable, use specific gravity to estimate pseudo-critical properties:

Sutton Correlation (1985): T_pc = 169.2 + 349.5×SG - 74.0×SG² (°R) P_pc = 756.8 - 131.0×SG - 3.6×SG² (psia) Valid for SG = 0.57 to 1.68 (sweet gas) Accuracy: ±1% for most natural gases Example (SG = 0.65): T_pc = 169.2 + 349.5(0.65) - 74.0(0.65)² = 365°R P_pc = 756.8 - 131.0(0.65) - 3.6(0.65)² = 670 psia

Dranchuk-Abou-Kassem (DAK) Correlation

The DAK correlation provides an explicit equation to calculate Z-factor, widely used in pipeline simulation software:

DAK Z-Factor Correlation: ρ_r = 0.27 × P_r / (Z × T_r) Z = 1 + (A₁ + A₂/T_r + A₃/T_r³ + A₄/T_r⁴ + A₅/T_r⁵) × ρ_r + (A₆ + A₇/T_r + A₈/T_r²) × ρ_r² - A₉(A₇/T_r + A₈/T_r²) × ρ_r⁵ + A₁₀(1 + A₁₁ρ_r²)(ρ_r²/T_r³) × exp(-A₁₁ρ_r²) Requires iterative solution (Newton-Raphson). Valid range: 0.2 < P_r < 30, 1.0 < T_r < 3.0 Accuracy: ±0.5% for natural gas

Typical Z-Factor Values

Application	Pressure (psia)	Temperature (°F)	Typical Z
Low-pressure distribution	50–100	40–100	0.98–1.00
Transmission pipeline	800–1000	50–80	0.85–0.90
High-pressure storage	2000–3000	60–100	0.75–0.85
Compressor suction	400–600	80–120	0.90–0.95
Compressor discharge	1200–1500	120–180	0.80–0.88

Impact of Z-factor: At 1000 psia and Z = 0.88, real gas is 14% denser than ideal gas prediction (ρ_real = ρ_ideal / 0.88 = 1.14 × ρ_ideal). Ignoring Z causes significant errors in flow rate and equipment sizing.

Z-factor versus pressure chart for natural gas (SG=0.65) showing curves at 60°F, 100°F, and 150°F with shaded regions indicating Ideal Gas OK zone below 100 psia, Consider Z zone 100-500 psia, and Z Required zone above 500 psia, with typical operating points marked — Z-factor vs. pressure for natural gas: Higher temperature → Z closer to 1.0. Shaded zones show when ideal gas is acceptable vs. when Z-factor correction is required.

Z-Factor Variations with Gas Composition

Dry gas (lean): High methane content → higher T_c relative to MW → Z closer to 1.0
Wet gas (rich): Higher ethane+ content → lower T_c → lower Z (more deviation)
Acid gas (CO₂/H₂S): High CO₂ or H₂S → significantly affects P_c and T_c → custom correlations needed
Nitrogen dilution: High N₂ → raises P_c → affects Z, especially at high pressure

4. Equations of State

Equations of state (EOS) relate pressure, volume, and temperature through thermodynamic models. More accurate than Z-factor charts for complex mixtures and extreme conditions.

Peng-Robinson EOS

Peng-Robinson Equation (1976): P = RT/(V - b) - a·α(T) / [V(V+b) + b(V-b)] Where: a, b = substance-dependent constants α(T) = temperature function V = molar volume Widely used in oil/gas industry for: - Phase equilibrium (VLE) - Density near critical point - Hydrocarbon mixtures with heavy components

AGA-8 Detail Method

The AGA-8 (American Gas Association Report No. 8) detail characterization method is the industry standard for custody transfer and high-accuracy applications.

AGA-8 Detail Method: Z = 1 + B/V + C/V² + D/V³ + E/V⁴ + F/V⁵ + G/V⁶ Where B, C, D, E, F, G are virial coefficients that depend on: - Gas composition (21 components) - Temperature - Density (solved iteratively) Accuracy: ±0.1% for natural gas mixtures Required inputs: - Composition (C1 through C10+, N2, CO2, H2S) - Pressure - Temperature

GERG-2008 EOS

GERG-2008 (Groupe Européen de Recherches Gazières) is the reference equation for natural gas, similar to AGA-8 but with extended range and 21 components.

Method	Accuracy	Application	Complexity
Ideal Gas Law	± 1–10%	Low P, preliminary calcs	Very simple
Z-factor (Standing-Katz)	± 1–2%	Pipeline, general design	Simple (chart or correlation)
CNGA (Dranchuk)	± 0.5–1%	Pipeline, compressor calcs	Moderate (iterative)
Peng-Robinson	± 0.5–2%	Phase behavior, VLE	Complex (requires solver)
AGA-8 Detail	± 0.1%	Custody transfer, metering	Very complex (iterative)
GERG-2008	± 0.05–0.1%	Reference standard, LNG	Very complex

When to Use Each Method

Ideal gas: Screening studies, low-pressure (<100 psia), academic problems
Z-factor (Standing-Katz or CNGA): Pipeline design, compressor sizing, general engineering
Peng-Robinson: Multiphase flow, gas processing, near critical conditions
AGA-8 or GERG: Custody transfer, revenue metering, high-accuracy requirements

Method selection: For pipeline hydraulics and equipment sizing, Z-factor methods (Standing-Katz or CNGA) provide sufficient accuracy. For custody transfer and revenue metering, AGA-8 is the contractual standard.

5. Practical Applications

Pipeline Flow Rate Calculations

Gas density directly affects mass flow rate and Reynolds number:

Mass Flow Rate: ṁ = ρ × Q Where: ṁ = Mass flow rate (lb/hr or kg/s) ρ = Gas density (lb/ft³ or kg/m³) Q = Volumetric flow rate (ft³/hr or m³/s) Reynolds Number: Re = ρ V D / μ Where: V = Gas velocity (ft/s) D = Pipe diameter (ft) μ = Dynamic viscosity (lb/ft·s) Re determines flow regime (laminar vs turbulent) and friction factor.

Orifice Meter Flow Calculation

Orifice meters measure differential pressure to infer flow rate:

Orifice Flow Equation (AGA Report 3): Q = C × E × Y × d² × √(ΔP / ρ) Where: Q = Volumetric flow rate at flowing conditions (ft³/hr) C = Discharge coefficient E = Velocity of approach factor Y = Expansion factor d = Orifice diameter (in) ΔP = Differential pressure (in H2O or psi) ρ = Gas density at flowing P/T (lb/ft³) A 1% error in ρ causes 0.5% error in calculated flow rate.

Compressor Power Calculation

Compressor power depends on inlet density and compression ratio:

Adiabatic Compression Power: HP = (Q × ρ₁ × R × T₁ × Z_avg / (MW × 33000)) × [(k/(k-1)) × ((P₂/P₁)^((k-1)/k) - 1)] Where: Q = Inlet volumetric flow (ft³/min) ρ₁ = Inlet density (lb/ft³) k = Specific heat ratio (Cp/Cv ≈ 1.27 for natural gas) P₂/P₁ = Compression ratio Inlet density affects both suction volume handling and power requirement.

Standard Volume Conversion

Convert actual flow to standard conditions for contracts and custody transfer:

Standard Volume: Q_std = Q_actual × (P_actual / P_std) × (T_std / T_actual) × (Z_std / Z_actual) Common standard conditions: - US: 14.73 psia, 60°F (Z_std ≈ 1.0) - ISO: 101.325 kPa abs, 15°C Example: Q_actual = 10 MMscfd at 800 psia, 80°F, Z = 0.88 Q_std = 10 × (800/14.73) × (520/540) × (1.0/0.88) Q_std = 10 × 54.31 × 0.963 × 1.136 Q_std = 594 MMscfd flowing equivalent

Erosional Velocity Check

API RP 14E provides erosional velocity limit to prevent pipe erosion:

Erosional Velocity (API RP 14E): V_erosion = C / √ρ Where: V_erosion = Maximum safe velocity (ft/s) C = Empirical constant (100 for non-corrosive, 125 for clean gas) ρ = Gas density (lb/ft³) For ρ = 0.3 lb/ft³: V_erosion = 100 / √0.3 = 183 ft/s Higher density → lower allowable velocity → larger pipe diameter required.

Common Pitfalls

Using gauge pressure instead of absolute: Always add atmospheric pressure (14.7 psia at sea level)
Mixing °F and °R: Gas law requires absolute temperature (°R = °F + 459.67)
Ignoring Z-factor at high pressure: Z < 1.0 makes gas denser than ideal prediction
Assuming constant density in long pipelines: Density varies with P/T profile along line
Using air MW (28.97) for natural gas: Natural gas MW typically 16–22 depending on composition
Using Z = 1.0 for custody transfer: Unacceptable for revenue metering—use AGA-8

Bar chart showing percent error in gas density when using ideal gas law versus real gas at pressures from 100 to 2000 psia, with color gradient from green (acceptable) to red (unacceptable), tolerance lines at 2% and 5%, and formula showing error equals (1/Z - 1) times 100% — Density error from ignoring Z-factor: 2% at 100 psia (acceptable) to 32% at 2000 psia (unacceptable). At 1000 psia, 16% error directly impacts mass flow and revenue calculations.

6. Specific Heat Ratio (k)

Density is only one of the gas properties needed for equipment design. The specific heat ratio (also called the isentropic exponent or gamma) governs compression work, discharge-temperature rise, and the speed of sound — it is essential for compressor head and power calculations.

Definition

Specific Heat Ratio: k = C_p / C_v Where: C_p = specific heat at constant pressure C_v = specific heat at constant volume

Where k Appears in Compressor Calculations

Head calculation: the (k-1)/k term in the polytropic/adiabatic head equation
Temperature rise: discharge temperature prediction
Efficiency conversion: between polytropic and isentropic efficiency
Speed of sound: c = √(k·Z·R·T / MW), which limits impeller tip speed (Mach number)

Typical k-Values

Gas	k-Value	Notes
Methane (CH₄)	1.31	Primary natural gas component
Natural Gas	1.26–1.32	Depends on composition
Nitrogen (N₂)	1.40	Diatomic gas
Carbon Dioxide (CO₂)	1.29	Triatomic, more complex
Hydrogen (H₂)	1.41	Light diatomic gas
Rich Gas (C₂+)	1.15–1.25	Heavier components lower k

Temperature, Pressure & Composition Effects

Unlike ideal gases where k is constant, real-gas k varies with conditions:

Temperature: k generally decreases with increasing temperature
Pressure: k changes significantly near critical conditions
Composition: heavier components decrease k

Don't use ideal-gas k at high pressure: Using ideal-gas k values at high pressures introduces significant errors. Always calculate real-gas C_p and C_v from an equation of state for accurate head and discharge-temperature results. For quick screening of natural gas, k ≈ 1.28 is a reasonable assumption.

7. Gas Properties for Compressor Calculations

Compressor sizing is highly sensitive to gas-property accuracy. A 5% error in Z-factor can produce a 5–10% error in head, potentially undersizing or oversizing the machine; an incorrect k gives the wrong discharge temperature. This section collects the properties and corrections most relevant to compression.

Molecular Weight Impact on Compressor Performance

For a gas mixture, the apparent molecular weight is the mole-fraction weighted average (Kay's rule): MW = Σ(y_i × MW_i), and SG = MW / 28.966.

Head requirement: head is inversely proportional to MW for the same pressure ratio
Number of stages: lighter gases require more head (more stages) per unit pressure ratio
Impeller tip speed: limited by MW (Mach-number concerns)
Power: MW cancels from the power equation at the same volumetric flow and pressure ratio; power depends primarily on k-value, not MW directly

Heavy vs. light gas: A compressor designed for natural gas (MW ≈ 18) achieves a much lower pressure ratio on hydrogen (MW = 2) at the same head, since lighter gases require far more head to achieve the same compression. Conversely, on propane (MW = 44), the same head produces a higher pressure ratio.

Z-Factor at Multiple Conditions (Iteration)

Centrifugal compressor calculations need Z at more than one point:

Z₁: at suction conditions (inlet density)
Z₂: at discharge conditions (discharge density)
Z_avg: average for head calculation = (Z₁ + Z₂) / 2

Iteration required: Since Z₂ depends on discharge temperature, which depends on head, which depends on Z_avg, an iterative solution is needed. Start with an assumed Z₂, calculate T₂, recalculate Z₂, and repeat until converged.

Wichert-Aziz Acid-Gas Correction

For sour gas containing H₂S or CO₂, the pseudo-critical properties from Kay's rule (or Sutton, §3) must be corrected before computing reduced properties and Z:

Wichert-Aziz Correction: ε = 120(A^0.9 − A^1.6) + 15(B^0.5 − B^4) T_pc' = T_pc − ε P_pc' = P_pc × T_pc' / [ T_pc + B(1 − B)ε ] Where: A = mol fraction CO₂ + mol fraction H₂S B = mol fraction H₂S ε = pseudo-critical temperature adjustment (°R)

Soave-Redlich-Kwong (SRK) Equation of State

SRK complements the Peng-Robinson EOS covered in §4. It is widely used in process simulation for vapor-phase property calculations:

Soave-Redlich-Kwong (SRK): P = RT/(V − b) − a(T) / [ V(V + b) ] Good for vapor-phase calculations and widely used in process simulators. Peng-Robinson (§4) generally gives better liquid-density predictions; SRK is often preferred for gas-phase compressor work.

A properly tuned EOS outputs all the properties a compressor calculation needs: Z, density (ρ), enthalpy (H), entropy (S), heat capacities (C_p, C_v, k), speed of sound (c), and the Joule-Thomson coefficient.

Worked Example — Multi-Component Natural Gas Density

Given: Sweet natural gas (mole fractions): CH₄ = 0.85, C₂H₆ = 0.10, C₃H₈ = 0.05. Operating at P = 1000 psia, T = 100°F (559.67 °R). Step 1 — Molecular weight (Kay's rule): MW = 0.85(16.043) + 0.10(30.07) + 0.05(44.097) = 18.85 lb/lbmol Step 2 — Specific gravity: SG = 18.85 / 28.966 = 0.651 Step 3 — Pseudo-critical properties (Kay's rule): T_pc = 0.85(343.37) + 0.10(549.92) + 0.05(665.68) = 380.1 °R P_pc = 0.85(667.8) + 0.10(708.3) + 0.05(616.3) = 669.3 psia (No H₂S/CO₂, so Wichert-Aziz correction is skipped.) Step 4 — Reduced properties: T_r = 559.67 / 380.1 = 1.472 ; P_r = 1000 / 669.3 = 1.494 Step 5 — Compressibility factor: Iterating the DAK correlation (or reading Standing-Katz) at these reduced conditions gives Z ≈ 0.854. Step 6 — Gas density: ρ = P·MW / (Z·R·T) = (1000 × 18.85) / (0.854 × 10.7316 × 559.67) ρ = 3.68 lbm/ft³ Real-gas correction: the ideal-gas density would be 3.14 lbm/ft³ — the Z-factor recovers the ~17% mass that an ideal-gas assumption would have lost in compressor sizing.

References

GPSA Engineering Data Book, Section 23 — Physical Properties
API 617 — Axial and Centrifugal Compressors
Campbell, J.M. — Gas Conditioning and Processing
Standing, M.B. and Katz, D.L. — Density of Natural Gases (1942)
Dranchuk, P.M. and Abou-Kassem, J.H. — Calculation of Z Factors
Wichert, E. and Aziz, K. — Calculate Z's for Sour Gases (1972)

Natural Gas Density: GPSA Engineering Fundamentals

Low P, High T

High P or Low T

0.7–1.0

Contents

Use this guide when you need to: