Equation of State Overview

Introduction

Equations of State (EoS) provide a rigorous thermodynamic framework for predicting fluid phase behavior from composition. Unlike empirical PVT correlations that treat oil as a "black box," EoS methods model individual hydrocarbon components and their interactions.

EoS modeling is essential when:

  • Compositional effects matter — gas injection, condensate reservoirs, volatile oils
  • Phase boundaries are needed — bubble/dew point curves, phase envelopes
  • Multiple phases coexist — vapor-liquid equilibrium (VLE), liquid-liquid equilibrium (LLE)
  • Process simulation — separator optimization, pipeline transport
  • PVT lab tuning — matching experimental data with thermodynamic models

When to Use EoS vs. Black Oil

Loading diagram...
Scenario Recommended Approach
Primary depletion, black oil PVT correlations
Gas condensate reservoir EoS (PR or SRK)
Volatile oil (GOR > 1000) EoS preferred
Gas injection (CO2, N2) EoS required
Compositional gradient EoS required
Near-critical fluids EoS required
Simple screening PVT correlations

Cubic Equations of State

General Form

All cubic EoS can be written as:

P=RTVba(T)(V+δ1b)(V+δ2b)P = \frac{RT}{V - b} - \frac{a(T)}{(V + \delta_1 b)(V + \delta_2 b)}

Where:

  • a(T)a(T) = attraction parameter (temperature-dependent)
  • bb = co-volume parameter (molecular size)
  • δ1,δ2\delta_1, \delta_2 = EoS-specific constants

Available Models

EoS Year δ1\delta_1 δ2\delta_2 Strength
Peng-Robinson (PR) 1976 1+21 + \sqrt{2} 121 - \sqrt{2} Best for liquid density
Soave-Redlich-Kwong (SRK) 1972 1 0 Better for vapor properties

Best Practice: Peng-Robinson is the industry standard for petroleum applications. Use SRK when comparing with legacy studies or when vapor-phase accuracy is more important.


Key Concepts

Component Properties

Each component in an EoS model requires:

Property Symbol Description
Critical temperature TcT_c Temperature above which no liquid phase exists
Critical pressure PcP_c Pressure at the critical point
Acentric factor ω\omega Measure of molecular non-sphericity
Molecular weight MwM_w For density calculations

For well-defined components (C1-C10, CO2, N2, H2S), these are known from databases. For heavy fractions (C7+), they must be estimated using characterization correlations.

Binary Interaction Parameters ($k_

Binary interaction parameters correct the geometric mean mixing rule for unlike molecular interactions:

aij=aiaj(1kij)a_{ij} = \sqrt{a_i \cdot a_j} \cdot (1 - k_{ij})

Pair Type Typical kijk_{ij} Notes
Hydrocarbon-hydrocarbon 0.0 - 0.02 Often set to zero
N2-hydrocarbon 0.02 - 0.12 Increases with carbon number
CO2-hydrocarbon 0.10 - 0.15 Important for CO2 injection
H2S-hydrocarbon 0.05 - 0.10 Moderate interaction

Alpha Functions

The temperature dependence of the attraction parameter is modeled through alpha functions:

a(T)=acα(Tr,ω)a(T) = a_c \cdot \alpha(T_r, \omega)

Alpha Function Form Best For
Soave (1972) (1+m(1Tr))2(1 + m(1 - \sqrt{T_r}))^2 General hydrocarbons
Graboski-Daubert (1978) Modified Soave coefficients Improved accuracy

EoS Workflow

Complete Phase Behavior Analysis

Loading diagram...

Available Calculations

Component Data

Calculation Description
Component properties TcT_c, PcP_c, ω\omega, MwM_w from database
SCN properties Single Carbon Number group properties
Component list Available components in database

Flash Calculations

Calculation Description
PT Flash (PR) Phase split at given P, T using Peng-Robinson
PT Flash (SRK) Phase split at given P, T using SRK
Bubble point Pressure/temperature where first gas appears
Dew point Pressure/temperature where first liquid appears

📖 Full Documentation: Flash Calculations

Phase Properties

Property PR SRK Description
Z-factor Yes Yes Compressibility factor per phase
Density Yes Yes Mass density per phase
Molar volume Yes Yes Volume per mole
Fugacity Yes Yes Thermodynamic fugacity
K-values Yes Yes Equilibrium ratios Ki=yi/xiK_i = y_i/x_i

Phase Envelope

Calculation Description
Envelope (PR) Complete phase boundary using Peng-Robinson
Envelope (SRK) Complete phase boundary using SRK

📖 Full Documentation: Phase Envelope


C7+ Characterization

Heavy fractions (C7+) are not pure components and must be characterized:

  1. Property estimation — estimate TcT_c, PcP_c, ω\omega from MwM_w and γ\gamma (specific gravity)
  2. Splitting — divide C7+ into pseudo-components using distribution functions
  3. Lumping — reduce component count for computational efficiency

Available characterization methods include Kesler-Lee, Twu, and Riazi-Daubert correlations.

📖 Full Documentation: C7+ Characterization


Viscosity Models

EoS provides volumetric properties but not transport properties. Separate viscosity models are needed:

Model Type Best For
Lee-Gonzalez-Eakin Gas viscosity Natural gas systems
Lorentz-Bray-Clark (LBC) Corresponding states General hydrocarbon mixtures
Pedersen Corresponding states Heavy oil, wide composition range

📖 Full Documentation: Viscosity Models


Best Practices

Model Selection

  1. Start with Peng-Robinson — industry standard, best liquid density
  2. Use SRK only when comparing with legacy data or for vapor-dominated systems
  3. Always specify kijk_{ij} for non-hydrocarbon components (CO2, N2, H2S)

Validation

  1. Match saturation pressure — tune kijk_{ij} and C7+ properties to match measured PbP_b or PdewP_{dew}
  2. Check density — compare calculated liquid density with measured values
  3. Verify phase envelope — ensure critical point and cricondenbar are physically reasonable
  4. Cross-check with PVT — black oil properties from EoS should agree with PVT correlations

Common Pitfalls

Issue Cause Solution
Wrong saturation pressure Poor C7+ characterization Tune MwM_w and γ\gamma of heaviest fraction
Liquid density off Volume translation needed Apply Peneloux correction
Flash doesn't converge Near critical point Use stability analysis first
Too many components Slow computation Lump to 7-12 pseudo-components

EoS Details

Supporting Topics


References

  1. Peng, D.Y. and Robinson, D.B. (1976). "A New Two-Constant Equation of State." Industrial & Engineering Chemistry Fundamentals, 15(1), 59-64.

  2. Soave, G. (1972). "Equilibrium Constants from a Modified Redlich-Kwong Equation of State." Chemical Engineering Science, 27(6), 1197-1203.

  3. Whitson, C.H. and Brule, M.R. (2000). Phase Behavior. SPE Monograph Vol. 20.

  4. Pedersen, K.S. and Christensen, P.L. (2007). Phase Behavior of Petroleum Reservoir Fluids. CRC Press.

  5. Michelsen, M.L. and Mollerup, J.M. (2007). Thermodynamic Models: Fundamentals and Computational Aspects, 2nd Edition. Tie-Line Publications.

An unhandled error has occurred. Reload X