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 (TcT_c, PcP_c, ω\omega) 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 MwM_w GC analysis or PVT report
Specific gravity γ\gamma Measured or estimated
Normal boiling point TbT_b Measured or estimated from MwM_w and γ\gamma

Correlation Methods

Kesler-Lee (1976)

Estimates TcT_c, PcP_c, ω\omega from TbT_b and γ\gamma:

Tc=341.7+811γ+(0.4244+0.1174γ)Tb+(0.46693.2623γ)×105TbT_c = 341.7 + 811\gamma + (0.4244 + 0.1174\gamma)T_b + \frac{(0.4669 - 3.2623\gamma) \times 10^5}{T_b}

Best for: General petroleum fractions, widely validated.

Twu (1984)

Uses a perturbation approach starting from n-alkane properties:

  1. Calculate n-alkane properties at the same TbT_b
  2. Apply perturbation corrections based on γ\gamma deviation from n-alkane

Best for: Heavy fractions, better extrapolation to high molecular weights.

Riazi-Daubert (1987)

Simple two-parameter correlations:

θ=aMwbγc\theta = aM_w^b \gamma^c

Where θ\theta represents TcT_c, PcP_c, or TbT_b, and aa, bb, cc 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:

p(M)=(Mη)α1βαΓ(α)exp(Mηβ)p(M) = \frac{(M - \eta)^{\alpha-1}}{\beta^\alpha \Gamma(\alpha)} \exp\left(-\frac{M - \eta}{\beta}\right)

Where:

  • η\eta = minimum molecular weight (typically 84 for C7)
  • α\alpha = shape parameter (controls distribution shape)
  • β=(Mw,C7+η)/α\beta = (M_{w,C7+} - \eta) / \alpha
α\alpha Value Distribution Shape Typical Fluid
α=1\alpha = 1 Exponential (most common) Gas condensates
α<1\alpha < 1 Skewed, heavy tail Heavy oils
α>1\alpha > 1 More symmetric Volatile oils

Katz Splitting

Assumes exponential decay of mole fractions with carbon number:

zn=z7exp[A(n7)]z_n = z_7 \exp[-A(n - 7)]

Where AA 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

  1. Group by volatility — components with similar K-values
  2. Preserve C1, CO2, N2, H2S as individual components
  3. Keep C2-C6 as pure components or small groups
  4. 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 (PbP_b, PdewP_{dew}) C7+ MwM_w, γ\gamma, or kijk_{ij} High
Liquid density Volume translation cic_i Direct
GOR / CGR C7+ splitting parameters Moderate
Phase volumes Multiple parameters Combined
  1. Match saturation pressure by adjusting C7+ molecular weight ($\pm$5-10%)
  2. Match liquid density by tuning volume translation
  3. Match gas-oil ratio by adjusting C7+ splitting
  4. Fine-tune kijk_{ij} for non-hydrocarbon components


References

  1. Whitson, C.H. (1983). "Characterizing Hydrocarbon Plus Fractions." SPE Journal, 23(4), 683-694. SPE-12233-PA.

  2. Kesler, M.G. and Lee, B.I. (1976). "Improve Prediction of Enthalpy of Fractions." Hydrocarbon Processing, 55(3), 153-158.

  3. 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.

  4. Riazi, M.R. and Daubert, T.E. (1987). "Characterization Parameters for Petroleum Fractions." Industrial & Engineering Chemistry Research, 26(4), 755-759.

  5. Pedersen, K.S., Christensen, P.L., and Shaikh, J.A. (2015). Phase Behavior of Petroleum Reservoir Fluids, 2nd Edition. CRC Press.

An unhandled error has occurred. Reload X