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
| 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:
Where:
- = attraction parameter (temperature-dependent)
- = co-volume parameter (molecular size)
- = EoS-specific constants
Available Models
| EoS | Year | Strength | ||
|---|---|---|---|---|
| Peng-Robinson (PR) | 1976 | 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 | Temperature above which no liquid phase exists | |
| Critical pressure | Pressure at the critical point | |
| Acentric factor | Measure of molecular non-sphericity | |
| Molecular weight | 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:
| Pair Type | Typical | 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:
| Alpha Function | Form | Best For |
|---|---|---|
| Soave (1972) | General hydrocarbons | |
| Graboski-Daubert (1978) | Modified Soave coefficients | Improved accuracy |
EoS Workflow
Complete Phase Behavior Analysis
Available Calculations
Component Data
| Calculation | Description |
|---|---|
| Component properties | , , , 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 |
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:
- Property estimation — estimate , , from and (specific gravity)
- Splitting — divide C7+ into pseudo-components using distribution functions
- 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
- Start with Peng-Robinson — industry standard, best liquid density
- Use SRK only when comparing with legacy data or for vapor-dominated systems
- Always specify for non-hydrocarbon components (CO2, N2, H2S)
Validation
- Match saturation pressure — tune and C7+ properties to match measured or
- Check density — compare calculated liquid density with measured values
- Verify phase envelope — ensure critical point and cricondenbar are physically reasonable
- 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 and 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 |
Related Documentation
EoS Details
- Peng-Robinson EoS — PR76 equation, mixing rules, parameters
- Flash Calculations — PT flash, Rachford-Rice, K-values
- Phase Envelope — Bubble/dew point curves, cricondenbar
- C7+ Characterization — Splitting, property estimation
- Viscosity Models — Lee-Gonzalez-Eakin, LBC, Pedersen
Supporting Topics
- PVT Overview — Empirical correlations for comparison
- VFP Overview — Multiphase flow requiring fluid properties
References
Peng, D.Y. and Robinson, D.B. (1976). "A New Two-Constant Equation of State." Industrial & Engineering Chemistry Fundamentals, 15(1), 59-64.
Soave, G. (1972). "Equilibrium Constants from a Modified Redlich-Kwong Equation of State." Chemical Engineering Science, 27(6), 1197-1203.
Whitson, C.H. and Brule, M.R. (2000). Phase Behavior. SPE Monograph Vol. 20.
Pedersen, K.S. and Christensen, P.L. (2007). Phase Behavior of Petroleum Reservoir Fluids. CRC Press.
Michelsen, M.L. and Mollerup, J.M. (2007). Thermodynamic Models: Fundamentals and Computational Aspects, 2nd Edition. Tie-Line Publications.