C7+ Characterization
Overview
Heavy hydrocarbon fractions (C7+) are not single components but complex mixtures of hundreds of molecules. Since equations of state require specific critical properties (, , ) for each component, the C7+ fraction must be characterized — broken into pseudo-components with estimated properties.
C7+ characterization is often the most important step in EoS modeling because:
- C7+ typically represents 20-60% of reservoir fluid by mole
- Phase behavior predictions are highly sensitive to heavy fraction properties
- Small changes in C7+ characterization can shift bubble/dew points by hundreds of psi
Property Estimation
Input Data
The minimum data needed for C7+ characterization:
| Property | Symbol | Source |
|---|---|---|
| Molecular weight | GC analysis or PVT report | |
| Specific gravity | Measured or estimated | |
| Normal boiling point | Measured or estimated from and |
Correlation Methods
Kesler-Lee (1976)
Estimates , , from and :
Best for: General petroleum fractions, widely validated.
Twu (1984)
Uses a perturbation approach starting from n-alkane properties:
- Calculate n-alkane properties at the same
- Apply perturbation corrections based on deviation from n-alkane
Best for: Heavy fractions, better extrapolation to high molecular weights.
Riazi-Daubert (1987)
Simple two-parameter correlations:
Where represents , , or , and , , are correlation-specific constants.
Best for: Quick estimates, widely used in refining applications.
Comparison
| Method | Strength | Weakness |
|---|---|---|
| Kesler-Lee | Well-validated, industry standard | Less accurate for very heavy fractions |
| Twu | Best for C20+ components | More complex computation |
| Riazi-Daubert | Simple, quick | Less accurate for reservoir fluids |
Splitting Methods
Purpose
A single C7+ pseudo-component is insufficient for accurate phase behavior. The C7+ fraction should be split into multiple pseudo-components that represent the molecular weight distribution.
Whitson Gamma Distribution
The molar distribution function follows a three-parameter gamma distribution:
Where:
- = minimum molecular weight (typically 84 for C7)
- = shape parameter (controls distribution shape)
| Value | Distribution Shape | Typical Fluid |
|---|---|---|
| Exponential (most common) | Gas condensates | |
| Skewed, heavy tail | Heavy oils | |
| More symmetric | Volatile oils |
Katz Splitting
Assumes exponential decay of mole fractions with carbon number:
Where is fitted to match the measured C7+ mole fraction and molecular weight. This is simpler than the Whitson method but less flexible.
Number of Pseudo-Components
| Application | Components | Notes |
|---|---|---|
| Screening | 1-3 | Quick estimates |
| Black oil reservoir simulation | 3-5 | Standard practice |
| Gas condensate | 5-8 | Need detail in light end |
| Compositional simulation | 7-12 | Balance accuracy and speed |
| PVT matching | 10-20 | Maximum detail for tuning |
Lumping
Purpose
After splitting, the total component count may be too large for efficient simulation. Lumping groups similar components to reduce the count while preserving phase behavior.
Guidelines
- Group by volatility — components with similar K-values
- Preserve C1, CO2, N2, H2S as individual components
- Keep C2-C6 as pure components or small groups
- Lump C7+ pseudo-components into 2-5 groups
Typical Lumped Composition
| Group | Components | Purpose |
|---|---|---|
| N2 | N2 | Separate — very different K-values |
| CO2 | CO2 | Separate — non-hydrocarbon interactions |
| C1 | CH4 | Lightest hydrocarbon |
| C2-C3 | C2H6 + C3H8 | Light intermediates |
| C4-C6 | iC4 + nC4 + iC5 + nC5 + C6 | Heavy intermediates |
| C7-C12 | Pseudo-components | Light heavy end |
| C13-C19 | Pseudo-components | Medium heavy end |
| C20+ | Pseudo-components | Heaviest fraction |
Tuning Strategy
Matching Measured Data
After initial characterization, the EoS model must be tuned to match laboratory PVT data:
| Data Type | Tuning Parameter | Sensitivity |
|---|---|---|
| Saturation pressure (, ) | C7+ , , or | High |
| Liquid density | Volume translation | Direct |
| GOR / CGR | C7+ splitting parameters | Moderate |
| Phase volumes | Multiple parameters | Combined |
Recommended Tuning Order
- Match saturation pressure by adjusting C7+ molecular weight ($\pm$5-10%)
- Match liquid density by tuning volume translation
- Match gas-oil ratio by adjusting C7+ splitting
- Fine-tune for non-hydrocarbon components
Related Topics
- EoS Overview — Role of characterization in EoS workflow
- Peng-Robinson EoS — How component properties feed into EoS
- Flash Calculations — Sensitivity of flash to characterization
- Phase Envelope — Impact on predicted phase boundaries
References
Whitson, C.H. (1983). "Characterizing Hydrocarbon Plus Fractions." SPE Journal, 23(4), 683-694. SPE-12233-PA.
Kesler, M.G. and Lee, B.I. (1976). "Improve Prediction of Enthalpy of Fractions." Hydrocarbon Processing, 55(3), 153-158.
Twu, C.H. (1984). "An Internally Consistent Correlation for Predicting the Critical Properties and Molecular Weights of Petroleum and Coal-Tar Liquids." Fluid Phase Equilibria, 16(2), 137-150.
Riazi, M.R. and Daubert, T.E. (1987). "Characterization Parameters for Petroleum Fractions." Industrial & Engineering Chemistry Research, 26(4), 755-759.
Pedersen, K.S., Christensen, P.L., and Shaikh, J.A. (2015). Phase Behavior of Petroleum Reservoir Fluids, 2nd Edition. CRC Press.