📑 Table of Contents (Tap to view sections)
1. Problem Statement & Industrial Context
A rich natural gas gathering pipeline delivers $15\text{ MMSCFD}$ ($17,700\text{ Sm}^3/\text{h}$) from an offshore production platform to a mainland processing terminal. The raw gas composition determined by on-line gas chromatography is:
- Methane ($\text{CH}_4$): $82.0\text{ mol\%}$
- Ethane ($\text{C}_2\text{H}_6$): $7.5\text{ mol\%}$
- Propane ($\text{C}_3\text{H}_8$): $4.5\text{ mol\%}$
- i-Butane ($\text{i-C}4\text{H}$): $1.2\text{ mol\%}$
- n-Butane ($\text{n-C}4\text{H}$): $1.8\text{ mol\%}$
- n-Pentane ($\text{C}5\text{H}$): $0.8\text{ mol\%}$
- Hexane-Plus ($\text{C}_6+$ as $n\text{-C}_6$): $0.7\text{ mol\%}$
- Nitrogen ($\text{N}_2$): $1.0\text{ mol\%}$
- Carbon Dioxide ($\text{CO}_2$): $0.5\text{ mol\%}$
Current Operating Conditions
- Platform Departure: $P_1 = 75.0\text{ barg}$ ($76.01\text{ bar abs}$), $T_1 = 35.0\text{ }^\circ\text{C}$
- Subsea River Crossing: $P_2 = 48.0\text{ barg}$ ($49.01\text{ bar abs}$), $T_2 = 18.0\text{ }^\circ\text{C}$ due to seabed ambient cooling and pipe frictional pressure drop.
Engineering Tasks
- Calculate the mixture pseudo-critical temperature ($T_c$) and pseudo-critical pressure ($P_c$).
- Determine the cricondenbar ($P_{ccb}$) and cricondentherm ($T_{cct}$).
- Evaluate the fluid phase state at both platform departure and the subsea crossing to verify whether the gas enters the retrograde condensation envelope.
- Quantify the expected liquid dropout volume and specify mitigation strategies.
2. Step-by-Step Mathematical Derivation
Step 1: Mixture Pseudo-Critical Properties (Kay's Rule)
The pseudo-critical temperature ($T_{pc}$) and pressure ($P_{pc}$) are evaluated as molar averages:
Using standard physical property constants: * Methane: $T_c = 190.56\text{ K}$, $P_c = 45.99\text{ bar}$ * Ethane: $T_c = 305.32\text{ K}$, $P_c = 48.72\text{ bar}$ * Propane: $T_c = 369.83\text{ K}$, $P_c = 42.48\text{ bar}$ * i-Butane: $T_c = 407.85\text{ K}$, $P_c = 36.40\text{ bar}$ * n-Butane: $T_c = 425.12\text{ K}$, $P_c = 37.96\text{ bar}$ * n-Pentane: $T_c = 469.70\text{ K}$, $P_c = 33.70\text{ bar}$ * n-Hexane: $T_c = 507.60\text{ K}$, $P_c = 30.25\text{ bar}$ * Nitrogen: $T_c = 126.20\text{ K}$, $P_c = 33.90\text{ bar}$ * Carbon Dioxide: $T_c = 304.13\text{ K}$, $P_c = 73.77\text{ bar}$
Computing the summation:
Similarly for pseudo-critical pressure:
The true thermodynamic critical point of the multicomponent mixture is shifted from the molar Kay's average due to the large disparity in component molecular sizes and acentric factors:
Step 2: Evaluation of Cricondentherm and Cricondenbar
Applying the Peng-Robinson EOS multicomponent bubble and dew point tracking algorithm: * Cricondenbar ($P_{ccb}$): The maximum two-phase pressure is determined to be: * Cricondentherm ($T_{cct}$): The maximum hydrocarbon dew point temperature occurs at an intermediate pressure of $38.0\text{ bar abs}$ ($37.0\text{ barg}$):
The retrograde condensation zone is bounded between: and extends from the dew point curve up to $P_{ccb} = 88.5\text{ bar}$.
Step 3: Phase State Classification at Operating Points
Operating Point 1 (Platform Departure: $76.01\text{ bar abs}$, $35.0\text{ }^\circ\text{C}$):
- $T_{op1} = 35.0\text{ }^\circ\text{C} > T_{cct} = 24.6\text{ }^\circ\text{C}$.
- Because the operating temperature exceeds the cricondentherm, no liquid can exist at any pressure.
- Phase Classification: Single-Phase Superheated Gas / Dense Phase. Margin above dew point = $+10.4\text{ }^\circ\text{C}$. Safe.
Operating Point 2 (Subsea River Crossing: $49.01\text{ bar abs}$, $18.0\text{ }^\circ\text{C}$):
- At $P = 49.01\text{ bar abs}$, we calculate the hydrocarbon dew point temperature ($T_{dp}$):
- The actual flowing gas temperature is $T_2 = 18.0\text{ }^\circ\text{C}$.
- Comparing temperatures:
- The operating point lies inside the two-phase retrograde condensation envelope!
- Phase Classification: Two-Phase Retrograde Condensation (Liquid Dropout Occurring). Margin: $-4.4\text{ }^\circ\text{C}$ below dew point.
Step 4: Hydrocarbon Liquid Dropout Volume
Performing an isothermal PR-EOS flash at $49.01\text{ bar abs}$ and $18.0\text{ }^\circ\text{C}$: * Liquid Molar Fraction ($L/F$): $1.85\text{ mol\%}$ of the stream condenses into liquid. * Liquid Density ($\rho_L$): $628.5\text{ kg/m}^3$ (stabilized pentanes, hexanes, and dissolved butanes). * Gas Flow Rate: $15.0\text{ MMSCFD} = 17,700\text{ Sm}^3/\text{h} = 12,240\text{ kg/h}$. * Liquid Mass Condensation Rate: * Volumetric Liquid Dropout:
Over $100\text{ bbl/day}$ of valuable hydrocarbon condensate drops out inside the subsea pipeline, creating severe liquid holdup, increased pipeline friction loss, and massive slugging risk during pipeline pigging.
3. Engineering Recommendations & Flow Assurance Actions
- Upstream Dew Point Control (DPCU): Install a Joule-Thomson refrigeration or mechanical chilling unit at the offshore platform to extract pentanes and heavier components ($\text{C}_5+$), lowering the gas cricondentherm to $-5\text{ }^\circ\text{C}$.
- Subsea Pipeline Insulation: If an upstream chilling unit cannot be retrofitted, insulate the subsea river crossing section (polyurethane coating) to maintain flowing temperatures above $25\text{ }^\circ\text{C}$.
- Downstream Slug Catcher Sizing: Ensure the receiving terminal slug catcher is sized to handle transient pigging surges of at least $50\text{ m}^3$ ($315\text{ bbl}$) of accumulated condensate.
Natural Gas Mixture Phase Envelope Explorer
Apply this methodology directly in the ChemProCal calculator.
Open Calculator →References
- Peng, D.Y. and Robinson, D.B. - A New Two-Constant Equation of State (Ind. Eng. Chem. Fundam., 15(1), 59-64, 1976) - Paper [Link ↗]
- GPSA Engineering Data Book, Section 25: Phase Equilibria (14th Edition, Gas Processors Suppliers Association) - Standard [Link ↗]
- Campbell, J.M. - Gas Conditioning and Processing, Vol 1: The Basic Principles (9th Edition, Campbell Petroleum Series, 2014) - Book
Live Gas Pipeline Velocity & Momentum Flux Estimator
Adjust parameters below to test the methodology equations in real time before running full simulations: