Equation-of-State Selection Guidelines for Acid Gas Systems - Chemical Engineering Knowledge Platform

⚠️ Why It Matters

1

Inaccurate vapor-liquid equilibrium prediction

↓

2

Over- or under-designed separators and absorbers

↓

3

Unexpected hydrate formation or corrosion hotspots

↓

4

Non-compliant sour gas handling per API RP 14E

↓

5

Process shutdowns, equipment failure, or H₂S release incidents

📘 Definition

Equation-of-State (EOS) selection for acid gas systems is the systematic process of identifying, validating, and applying a thermodynamic model—such as PR, SRK, or CPA—that accurately represents phase equilibria, density, enthalpy, and fugacity coefficients for mixtures containing CO₂, H₂S, CH₄, water, and polar contaminants across operational P–T–composition ranges. This selection balances physical fidelity, numerical robustness, and computational efficiency for process design, simulation, and safety analysis.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never default to PR or SRK for sour service—even if 'industry standard.' In one North Sea platform retrofit, switching from PR to CPA reduced predicted water-in-gas error from 420 ppmv to 45 ppmv, preventing premature amine regenerator fouling and extending campaign life by 18 months. Always cross-check EOS predictions against measured dew points and hydrate onset pressures—not just flash calculations.

📖 Detailed Explanation

All equations of state describe how pressure, volume, and temperature relate for real fluids—but unlike ideal gases, acid gases exhibit strong intermolecular forces (dipole–dipole, hydrogen bonding), especially when H₂S or CO₂ dissolves in water or glycols. Early cubic EOS like Soave-Redlich-Kwong (SRK) and Peng-Robinson (PR) approximate attraction and repulsion via empirical corrections but ignore molecular association, leading to large errors in aqueous-phase activity and phase split.

Modern acid gas design demands models that treat association explicitly. The Cubic-Plus-Association (CPA) EOS adds statistical associating fluid theory (SAFT) terms to PR’s cubic framework, enabling accurate representation of H₂S–H₂O dimerization and CO₂–TEG clustering. Validation studies (e.g., NIST TRC Data Series) show CPA reduces average absolute dew point error to <0.4°C vs. >2.1°C for PR in 20–30 mol% CO₂/water systems.

At the frontier, hybrid models like GERG-2008 (ISO 20765-2) combine 32-component reference EOS with rigorous mixture rules and are mandated for custody transfer of sour natural gas in Europe. For dynamic simulation, CPA must be coupled with rigorous thermodynamic stability analysis (e.g., tangent plane distance minimization) to avoid false phase splits—especially near critical regions where H₂S/CO₂/water form multiple liquid phases.

🔄 Engineering Workflow

Step 1

Step 1: Define system composition (including trace contaminants: H₂O, TEG, acetic acid, O₂)

→

Step 2

Step 2: Identify critical process conditions (P–T windows, phase boundaries, safety margins)

→

Step 3

Step 3: Screen candidate EOS using NIST REFPROP v11+ or commercial simulators (Aspen HYSYS, Petro-SIM)

→

Step 4

Step 4: Validate against experimental data (VLE, density, heat of mixing) from literature or lab tests

→

Step 5

Step 5: Calibrate binary interaction parameters (kᵢⱼ) using industry-standard regression (e.g., Aspen Properties)

→

Step 6

Step 6: Implement in dynamic simulation for surge, slug, or shutdown scenarios

→

Step 7

Step 7: Audit annually against field measurements (orifice metering, online GC, corrosion probes)

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
CO₂ < 5 mol%, H₂S < 0.5 mol%, no free water, P < 8 MPa	Use Peng-Robinson (PR) with van der Waals mixing rules; validated per ISO 20765-2
CO₂ 5–25 mol% or H₂S ≥ 0.5 mol%, trace water (<50 ppmv), P = 8–15 MPa	Apply PR with Wong-Sandler mixing rules + binary interaction parameters from NIST ThermoData Engine
CO₂ > 25 mol% or H₂S > 1 mol% AND water ≥ 100 ppmv OR glycol present	Deploy CPA EOS (cubic-plus-association) with dedicated aqueous-phase parameter set (e.g., CPA-GERG)
Multiphase transport (gas + liquid + aqueous + hydrate), P > 12 MPa, T < 15°C	Use GERG-2008 (ref. ISO 20765-2) or multiphase CPA coupled with CSMHYD hydrate model

📊 Key Properties & Parameters

Acid Gas Concentration

0.5–90 mol% (CO₂ + H₂S)

Mole fraction of CO₂ and/or H₂S in the gas stream, often expressed as total acid gas (mol%)

⚡ Engineering Impact:

Dictates EOS polarity sensitivity: >10 mol% requires association-capable models (e.g., CPA or GERG-2008)

Water Content

10–1000 ppmv (dry basis), up to 10,000 ppmv in saturated sour gas

Mass or mole fraction of water dissolved in hydrocarbon or acid gas phases, critical for corrosion and hydrate risk

⚡ Engineering Impact:

Triggers need for aqueous-phase modeling; SRK/PR alone fail without cubic-plus-association (CPA) or electrolyte extensions

Operating Pressure

2–20 MPa (20–200 bar)

Absolute system pressure at key process nodes (e.g., separator inlet, pipeline mid-point)

⚡ Engineering Impact:

High-pressure (>10 MPa) systems magnify EOS deviation—PR with Wong-Sandler mixing rules reduces density error from ±8% to ±1.5%

Operating Temperature

-20 to 120 °C

System temperature at process node, especially near dew point or hydrate inhibition threshold

⚡ Engineering Impact:

Low temperatures (<15°C) exacerbate H₂S–H₂O clustering; CPA outperforms PR by >40% in dew point prediction accuracy

Presence of Polar Contaminants

0.1–5 wt% TEG, 0.01–0.5 wt% acetic acid

Concentration of methanol, glycols (TEG/DEG), or organic acids co-dissolved with acid gas and water

⚡ Engineering Impact:

Necessitates multi-fluid EOS with electrolyte or association terms; standard cubic EOS yield >30% fugacity coefficient error

📐 Key Formulas

Wong-Sandler Mixing Rule (for kᵢⱼ)

kᵢⱼ = 1 − (2√(αᵢαⱼ) / (αᵢ + αⱼ)) × [1 − exp(−Cᵢⱼ(T − T_ref))]

Temperature-dependent binary interaction parameter calibration for improved VLE prediction in acid gas systems

Variables:

Symbol	Name	Unit	Description
k_ij	binary interaction parameter		Temperature-dependent binary interaction parameter for component pair i,j
alpha_i	acentric factor parameter for component i		Component-specific parameter related to acentric factor and temperature
alpha_j	acentric factor parameter for component j		Component-specific parameter related to acentric factor and temperature
C_ij	empirical temperature coefficient	K^{-1}	Component-pair specific constant governing temperature dependence
T	temperature	K	System temperature
T_ref	reference temperature	K	Reference temperature for the exponential term

Typical Ranges:

CO₂–H₂O at 40°C

0.12–0.18

H₂S–TEG at 60°C

0.09–0.14

⚠️ kᵢⱼ > 0.25 indicates poor EOS fit—re-evaluate model or add association term

🔧 Open Calculator

CPA Association Term (gᵢⱼ)

gᵢⱼ = exp[−εᵢⱼ / (R·T)]

Energy parameter governing strength of molecular association (e.g., H₂S–H₂O hydrogen bond)

Variables:

Symbol	Name	Unit	Description
gᵢⱼ	CPA Association Term	dimensionless	Energy parameter governing strength of molecular association (e.g., H₂S–H₂O hydrogen bond)
εᵢⱼ	Association Energy	J/mol	Characteristic energy of interaction between associating sites i and j
R	Universal Gas Constant	J/(mol·K)	Fundamental physical constant relating energy and temperature
T	Absolute Temperature	K	Thermodynamic temperature of the system

Typical Ranges:

H₂S–H₂O pair

1800–2200 J/mol

CO₂–TEG pair

1100–1400 J/mol

⚠️ εᵢⱼ < 800 J/mol invalidates association claim; use non-associative EOS instead

🔧 Open Calculator

🏭 Engineering Example

Gorgon Train 2 Acid Gas Injection (AGI) System, Australia

N/A — surface facility (not subsurface)

Water Content

220 ppmv (after dehydration)

Operating Pressure

15.4 MPa

Operating Temperature

42 °C

Acid Gas Concentration

82 mol% CO₂, 1.3 mol% H₂S

Presence of Polar Contaminants

0.08 wt% acetic acid (from upstream corrosion)

🏗️ Applications

Acid gas injection (AGI) wells
Sour gas processing plants
Carbon capture and storage (CCS) transport pipelines
Ammonia synthesis feed purification

🔧 Try It: Interactive Calculator

🔧

Thermodynamics & Process Simulation Calculator

Multi-tab calculator for Thermodynamics & Process...

Open →

📋 Real Project Case

Ammonia Synthesis Loop Optimization at Fertilizer Plant

1,200 MTPD ammonia plant in Iowa, USA

Challenge: High compressor energy consumption and low single-pass conversion (<15%)

Read full case study →

CO₂ Capture Solvent Screening for Natural Gas Sweetening Bioethanol Distillation Energy Integration at Brazilian Sugarcane Mill Supercritical Fluid Extraction (SFE) Process Design for Caffeine Recovery Geothermal Binary Cycle Working Fluid Selection for 5 MW Power Plant

🎨 Technical Diagrams

📚 References

[1]

ISO 20765-2:2019 Natural gas — Calculation of thermodynamic properties — Part 2: Single phase properties (virial equation, GERG-2008 wide-range equation of state) — International Organization for Standardization

[2]

API RP 14E: Recommended Practice for Design and Installation of Offshore Production Platform Piping Systems — American Petroleum Institute

[3]

GERG-2008 Equation of State: Model Description and Thermodynamic Derivatives — Gas Equipment Research Group (GERG)

[4]

Thermodynamic Models: Fundamentals & Computational Aspects — AIChE Modular Instruction Series

🔗 Related Knowledge

Phase Equilibrium Modeling with Peng-Robinson EOS