Gathering Network Optimization - Multiphase Flow Modeling, Separator Pressure Selection, and Compression Design

Gathering Network Optimization - Multiphase Flow Modeling, Separator Pressure Selection, and Compression Design

The gathering network is the collection of flowlines, trunk lines, separators, compressors, and pumps that transport produced fluids from individual wellheads to the central processing facility. It is the hydraulic link between the reservoir and the surface plant, and its design determines how much of the reservoir's productive capacity actually reaches the facility as saleable product. A gathering network that is over-designed - with oversized flowlines and excessive separation capacity - represents capital that could have been deployed elsewhere. A network that is under-designed creates backpressure on the wells that reduces production below the reservoir's deliverability, increases the risk of liquid loading in gas wells, creates slugging problems in multiphase flowlines, and limits the field's plateau production rate. The optimization of gathering network design is therefore a direct production optimization problem: for a given reservoir deliverability profile and a given capital budget, what combination of flowline sizes, separator pressures, and compression capacity maximizes the net present value of the field over its producing life? This guide covers the quantitative engineering framework that answers this question.

1. Multiphase Flow in Gathering Flowlines

1.1 Flow Regime Identification - The Foundation of Multiphase Design

Unlike single-phase flow where pressure drop depends only on velocity and friction, multiphase flow behavior depends critically on the flow regime - the spatial distribution of gas and liquid phases in the pipe. Two pipes with identical gas and liquid flow rates but different inclination angles or pipe diameters can exhibit completely different flow regimes with dramatically different pressure drop characteristics and operational behavior:

Flow Regime	Description	Typical Conditions	Operational Consequence
Bubble flow	Small gas bubbles dispersed in continuous liquid phase. Gas fraction (void fraction) typically <25%.	Low GOR, high liquid rate, high pressure (gas compressed into small bubbles)	Most stable flow regime. Predictable pressure drop. Liquid phase controls hydraulics. Favorable for production.
Slug flow	Alternating large gas bubbles (Taylor bubbles) and liquid slugs. Gas void fraction 25-75%. Most common regime in production flowlines.	Moderate GOR, intermediate velocities. Very common in horizontal and slightly inclined flowlines.	Highly variable pressure at separator inlet. Slug catcher required. Can cause liquid carryover in gas separators and gas breakthrough in liquid separators.
Stratified flow	Gas and liquid flow separately with a smooth or wavy interface. Gas on top, liquid on bottom. Only possible in horizontal or slightly downward-inclined pipes at low velocities.	Low velocity, horizontal pipe, moderate GOR	Stable but risk of corrosion at bottom of pipe (liquid contains CO2/H2S). Transition to slug flow with any velocity increase.
Annular flow	Gas core flows at high velocity in center of pipe. Liquid film flows on pipe wall. Gas void fraction >75%.	High GOR, high gas velocity, low liquid rate. Common in gas wells and gas transmission lines.	High erosion risk at elbows (liquid droplets in gas core). As rate declines, transition to slug flow signals onset of liquid loading.
Mist flow	Liquid dispersed as fine droplets in continuous gas phase. Extreme case of annular flow where liquid film is stripped from wall.	Very high gas velocity (>50 ft/sec), very low liquid holdup	Maximum erosion risk. Droplets impact pipe wall at high velocity at bends. Requires erosion-resistant materials or velocity reduction.

1.2 Beggs-Brill Pressure Drop Correlation - The Industry Workhorse

The Beggs-Brill correlation is the most widely used method for calculating pressure drop in multiphase flowlines. It was developed empirically from laboratory data covering horizontal, vertical, and inclined pipes and remains the default method in PIPESIM and most other multiphase flow simulators for flowline design:

Beggs-Brill method - key parameters:

Input variables:
q_L = liquid flow rate (bbl/day): oil + water at pipeline conditions
q_g = gas flow rate (MMscf/day) at pipeline conditions
D = pipe inside diameter (inches)
theta = pipe inclination angle (degrees from horizontal, positive = uphill)
rho_L = liquid density (lb/ft3), rho_g = gas density (lb/ft3)
mu_L = liquid viscosity (cp), mu_g = gas viscosity (cp)
sigma_L = liquid-gas surface tension (dynes/cm)

Step 1 - No-slip liquid holdup and mixture velocity:
lambda_L = q_L / (q_L + q_g) (input liquid fraction, no slip)
V_m = V_SL + V_SG (mixture superficial velocity, ft/sec)
V_SL = q_L / (2.448 x D^2) (liquid superficial velocity, ft/sec, with q_L in bbl/day)
V_SG = q_g x 4.706e5 / (D^2 x P_avg x z) (gas superficial velocity, ft/sec)

Step 2 - Froude number and flow regime:
N_FR = V_m^2 / (32.174 x D/12) (Froude number, dimensionless)

Step 3 - Liquid holdup correction for actual slip:
H_L (actual) = H_L(horizontal) x phi (inclination correction)

Example calculation - 8" flowline, 10 km horizontal, oil well:
Production: q_L = 3,500 bbl/day (2,800 oil + 700 water), GOR = 400 scf/STB
q_g_surface = 3,500 x 400 = 1,400,000 scf/day = 1.4 MMscf/day
Pipeline conditions: P_avg = 250 psia, T = 80°F, z = 0.92
D = 7.981" (8" schedule 40), rho_oil = 52 lb/ft3, rho_water = 64 lb/ft3
rho_L = (2800 x 52 + 700 x 64)/(2800+700) = (145,600 + 44,800)/3,500 = 54.4 lb/ft3

V_SL = 3,500 / (2.448 x 7.981^2) = 3,500 / 155.97 = 22.44 ft/day = 0.260 ft/sec
q_g_pipeline = 1,400,000 x 14.7/250 x (540/520) x 0.92 = 1,400,000 x 0.0589 x 1.038 x 0.92
= 1,400,000 x 0.0562 = 78,680 scf/min → actually in ft3/sec:
q_g = 1,400,000 scf/day / 86,400 sec/day x (14.7/250) x (540/520) / 0.92 = 16.20 x 0.0588 x 1.038 / 0.92
= 16.20 x 0.06637 = 1.075 ft3/sec = 1.075 / (pi/4 x (7.981/12)^2) = 1.075/0.3472 = 3.10 ft/sec V_SG

V_m = 0.260 + 3.10 = 3.36 ft/sec mixture velocity
lambda_L = 0.260 / 3.36 = 0.0774 (7.74% no-slip liquid fraction)

N_FR = 3.36^2 / (32.174 x 7.981/12) = 11.29 / 21.44 = 0.527

At N_FR = 0.527 and lambda_L = 0.0774: Flow regime = Slug flow (from Beggs-Brill regime map)

Actual liquid holdup H_L ≈ 0.18-0.25 for these conditions (gas slips relative to liquid)
Mixture density = H_L x rho_L + (1-H_L) x rho_g = 0.22 x 54.4 + 0.78 x 1.85 = 11.97 + 1.44 = 13.4 lb/ft3

Total pressure drop over 10 km (32,808 ft):
dP_friction ≈ f x rho_m x V_m^2 x L / (2 x D_m) [Darcy-Weisbach adapted for multiphase]
f = 0.020 (Moody friction factor at Re ~500,000), D_m = 7.981/12 = 0.665 ft
dP = 0.020 x 13.4 x 3.36^2 x 32,808 / (2 x 0.665) = 0.020 x 13.4 x 11.29 x 32,808 / 1.330
= 0.020 x 13.4 x 11.29 x 24,668 = 74,745 lb/ft2 = 519 psi total pressure drop

(Note: This is a simplified friction-dominated calculation. PIPESIM would iterate with full Beggs-Brill including holdup effects across multiple segments.)

2. Separator Pressure Optimization

2.1 Optimal First-Stage Separator Pressure

In a multi-stage separation train (typically three stages), the pressure at each separator determines how much gas is released at that stage and at what composition. The separator pressures affect the stock tank oil volume (API gravity), the gas-oil ratio at each stage, and the total liquid recovery from the gas. The optimal set of separator pressures maximizes the combined value of the oil and gas products:

Effect of first-stage separator pressure on liquid recovery:
High first-stage pressure (800-1,200 psia):
- More intermediate components (C3-C6) remain in liquid phase
- Higher stock tank oil volume (higher Bo)
- Lower GOR at first stage (less gas released)
- Gas at first stage has lower heating value (lighter composition)

Low first-stage pressure (200-400 psia):
- Intermediate components flash to gas
- Lower stock tank liquid volume
- Higher GOR at first stage
- Gas is richer (more C3-C6 NGL content)

Three-stage separation pressure optimization example:
Reservoir fluid: GOR = 850 scf/STB, C7+ = 28 mol%, API = 38°
Wellhead pressure: 2,200 psia

Optimize P1 (first stage), P2 (second stage), P3 = 50 psia (stock tank)

Flash calculation results (from equation of state simulation):

Case A: P1 = 1,200 psia, P2 = 200 psia:
Stock tank oil = 1.000 STB/STB (baseline)
Total gas = 850 scf/STB (baseline)
Oil API = 38.2°

Case B: P1 = 800 psia, P2 = 150 psia:
Stock tank oil = 1.012 STB/STB (+1.2% more oil)
Total gas = 843 scf/STB (-0.8% less gas)
Oil API = 38.7°

Case C: P1 = 500 psia, P2 = 100 psia:
Stock tank oil = 1.018 STB/STB (+1.8% more oil)
Total gas = 835 scf/STB (-1.8% less gas)
Oil API = 39.3°

Case D: P1 = 300 psia, P2 = 75 psia:
Stock tank oil = 1.015 STB/STB (+1.5% more oil - diminishing returns)
Total gas = 838 scf/STB (-1.4% less gas)
Oil API = 39.8°

Revenue comparison per STB wellhead production (at oil = $65/STB, gas = $3.50/Mscf):
Case A: 1.000 x $65 + 850/1,000 x $3.50 = $65.00 + $2.975 = $67.975/STB
Case B: 1.012 x $65 + 843/1,000 x $3.50 = $65.78 + $2.951 = $68.731/STB (+$0.756 vs Case A)
Case C: 1.018 x $65 + 835/1,000 x $3.50 = $66.17 + $2.923 = $69.093/STB (+$1.118 vs Case A)
Case D: 1.015 x $65 + 838/1,000 x $3.50 = $65.975 + $2.933 = $68.908/STB (+$0.933 vs Case A)

Case C (P1=500 psia) maximizes revenue at $69.093/STB
For a 10,000 STB/day field: Case C generates $1,118 x 10,000 = $11,180/day = $4.08 million/year additional revenue vs the default high-pressure case.

3. Gathering Network Pressure Management

3.1 Backpressure Effect on Well Production Rate

Every psi of gathering system pressure at the wellhead is backpressure that reduces the pressure differential available for production from the reservoir. For wells on natural flow, gathering pressure is directly subtracted from the available drawdown. For wells on artificial lift, gathering pressure increases the total dynamic head requirement and therefore the lift energy cost:

Production rate sensitivity to wellhead pressure (backpressure):
Using inflow performance relationship (IPR) with Vogel equation for solution gas drive reservoir:
q/q_max = 1 - 0.20 x (Pwf/Pr) - 0.80 x (Pwf/Pr)^2

Where q_max = AOF (Absolute Open Flow), Pr = reservoir pressure

Pwf = Pwh + dP_tubing + dP_perforations
dP_tubing (vertical lift) = rho_fluid x 0.052 x TVD - (small correction for friction)
For a 6,000 ft TVD well producing 800 bbl/day at 45% WC: dP_tubing ≈ 1,200 psi

Current wellhead pressure: Pwh = 150 psi
Pwf = 150 + 1,200 = 1,350 psi

Pr = 2,800 psi, q_max = 3,200 bbl/day
q = 3,200 x [1 - 0.20 x (1,350/2,800) - 0.80 x (1,350/2,800)^2]
= 3,200 x [1 - 0.20 x 0.482 - 0.80 x 0.2323]
= 3,200 x [1 - 0.0964 - 0.1859] = 3,200 x 0.7177 = 2,297 bbl/day current production

Effect of reducing wellhead pressure from 150 to 75 psi (gathering system upgrade):
New Pwf = 75 + 1,200 = 1,275 psi
q_new = 3,200 x [1 - 0.20 x (1,275/2,800) - 0.80 x (1,275/2,800)^2]
= 3,200 x [1 - 0.20 x 0.4554 - 0.80 x 0.2073]
= 3,200 x [1 - 0.0911 - 0.1659] = 3,200 x 0.7430 = 2,378 bbl/day new production

Production increase = 2,378 - 2,297 = 81 bbl/day per well from 75 psi wellhead pressure reduction

For a 25-well field: 25 x 81 = 2,025 bbl/day additional production
Annual revenue = 2,025 x 365 x $60 = $44.3 million/year additional revenue from the gathering pressure reduction.

The capital cost to reduce gathering pressure by 75 psi across a 25-well field (larger trunk line, additional compression) must be compared to this $44.3M/year revenue benefit to justify the investment.

3.2 Slugging in Gathering Networks - Diagnosis and Mitigation

Terrain-induced slugging occurs when liquid accumulates at low points (valleys) in the flowline topography during periods of low flow rate. When sufficient gas pressure builds behind the liquid accumulation, the liquid slug is expelled as a large pulse that can exceed the separator inlet capacity, cause liquid overflow, and trip production equipment on high-level shutdowns:

Slug volume estimation (severe slugging in hilly terrain):
V_slug (bbls) ≈ Volume of pipeline low-point segment x liquid holdup fraction

For a 6" flowline low-point segment 2 km long with holdup H_L = 0.45:
V_pipe = pi/4 x (6.065/12)^2 x (2,000 x 3.281) ft = 0.2006 x 6,562 = 1,316 ft3 = 234 bbls
V_slug = 234 x 0.45 = 105 bbls per slug event

Slug frequency (severe slugging cycle time):
T_slug = V_slug / (q_liquid at pipeline conditions)
q_L = 800 bbl/day = 33.3 bbl/hr
T_slug = 105/33.3 = 3.15 hours per slug cycle

Peak slug arrival rate at separator:
Slug arrives over approximately 10-20 minutes (rapid expulsion)
Peak rate = 105 bbls / (15 min / 60 min/hr) = 420 bbl/hr = 10,080 bbl/day peak slug rate

Average design rate: 800 bbl/day
Peak slug rate: 10,080 bbl/day → 12.6x average rate

This peak rate requires the separator to have either:
1. A liquid slug catcher rated for 10,000+ bbl/day inlet surge, OR
2. Terrain-following flowline rerouting to eliminate the low point, OR
3. Gas lift injection at the low point to continuously aerate the liquid and prevent accumulation, OR
4. Active slug control valve at the wellhead that throttles gas flow to prevent slug buildup (requires real-time monitoring)

4. Compression System Design

4.1 Compression Ratio and Stage Design

Gathering system compression is required when reservoir pressure has declined to the point where natural flow can no longer overcome the gathering system pressure plus the wellbore hydraulic losses. The compression design must account for the full life-of-field compression requirement, including the declining suction pressure as the reservoir depletes and the increasing water-gas ratio that affects the gas composition and dew point entering the compressor:

Compression power calculation (polytropic compression):
W_poly (kW) = (n/(n-1)) x q_g x P1 x [(P2/P1)^((n-1)/n) - 1] / (eta_poly x eta_mech)

Where:
n = polytropic index (typically 1.2-1.4 for natural gas)
q_g = actual inlet gas flow rate (m3/s at suction conditions)
P1 = suction pressure (Pa), P2 = discharge pressure (Pa)
eta_poly = polytropic efficiency (0.75-0.82 for centrifugal, 0.80-0.88 for reciprocating)
eta_mech = mechanical efficiency (0.95-0.98)

Practical field formula (in oilfield units):
HP = 0.0857 x q_g (MMscfd) x T1 (°R) / eta_overall x [(P2/P1)^0.2857 - 1]

Example: Gas gathering compression
q_g = 12 MMscf/day, T1 = 80°F = 540°R
P1 = 80 psia (suction, low pressure gathering), P2 = 600 psia (discharge to sales pipeline)
Compression ratio = 600/80 = 7.5
eta_overall = 0.78 (centrifugal with gas cooler)

For compression ratio > 4, two stages are recommended to limit discharge temperature:
Optimal inter-stage pressure = sqrt(P1 x P2) = sqrt(80 x 600) = sqrt(48,000) = 219 psia inter-stage pressure
Stage 1: 80 → 219 psia (ratio = 2.74)
Stage 2: 219 → 600 psia (ratio = 2.74)

HP per stage = 0.0857 x 12 x 540 / 0.78 x [(2.74)^0.2857 - 1]
= 0.0857 x 12 x 540 / 0.78 x [1.3145 - 1]
= 0.0857 x 12 x 540 / 0.78 x 0.3145
= 554.9 / 0.78 x 0.3145 = 711.4 x 0.3145 = 223.7 HP per stage
Total for 2 stages: 447 HP = 333 kW total compression power

Discharge temperature check (to prevent lubricant degradation and valve damage):
T2 = T1 x (P2/P1)^((n-1)/n) = 540 x (2.74)^(0.4/1.4) = 540 x (2.74)^0.286 = 540 x 1.315 = 710°R = 250°F
Maximum allowable: 250-300°F for most compressor designs → at limit, gas cooler between stages is mandatory.

4.2 Life-of-Field Compression Staging Strategy

Field Life Phase	Reservoir Pressure	Required Compression	Strategy
Early production (0-5 years)	High (>2,500 psi)	None or minimal	Natural flow to gathering system. Wellhead chokes control rate. Compression facilities on standby or not yet installed.
Mid-life (5-12 years)	Moderate (1,200-2,500 psi)	Low-pressure gathering compression	First compression stage installed. Reduces wellhead backpressure from 150 to 50 psi, extending plateau production rate by 2-4 years.
Late production (12-20 years)	Low (<1,200 psi)	Multi-stage compression	Second and possibly third compression stage. Wellhead pressure may be as low as 20-30 psia. Artificial lift on most wells. Gas-liquid separation at wellhead to reduce liquid loading on compressors.
Abandonment phase (>20 years)	Very low (<400 psi)	Maximum compression ratio	Economic limit reached when compression operating cost exceeds revenue from produced gas. P&A decision triggered by compression economics rather than reservoir depletion.

5. Network Simulation and Optimization - PIPESIM Workflow

5.1 Building the Network Model

A full gathering network simulation links the inflow performance of each well (IPR curves from reservoir data) to the surface hydraulics of the flowlines, manifolds, separators, and compression system. The simulation finds the operating point where all wells simultaneously satisfy both their IPR and the gathering system hydraulics:

Network simulation workflow (PIPESIM or equivalent):

Step 1: Input well data
- IPR for each well: Vogel or Jones equation with current reservoir pressure
- Well completion data: tubing size, depth, artificial lift specification
- Produced fluid composition: GOR, water cut, oil gravity, gas SG

Step 2: Build surface network
- Flowline geometry: length, diameter, elevation profile, insulation
- Manifold connections: which wells connect to which trunk lines
- Separator: pressure, temperature, efficiency
- Compression: suction pressure, discharge pressure, capacity

Step 3: Specify boundary conditions
- Delivery pressure at custody transfer point (fixed)
- Ambient temperature (for heat transfer calculations)

Step 4: Run nodal analysis
The simulator iterates until pressures and flow rates satisfy mass and energy balance simultaneously at every node in the network.

Step 5: Sensitivity analysis
- Vary separator pressure: find optimal P1 for maximum revenue
- Vary flowline diameter: find economic optimum (production gain vs capital cost)
- Vary compression suction pressure: find pressure that maximizes field NPV

Key output metrics:
- Production rate by well and by total field
- Wellhead pressure for each well
- Pressure drop across each flowline segment
- Liquid holdup and flow regime in each segment
- Slug volume at separator inlet
- Compression power requirement
- Dew point and hydrate risk assessment for each segment

Conclusion

The separator pressure optimization in this article - Case C at P1 = 500 psia generating $69.093/STB versus Case A at P1 = 1,200 psia generating $67.975/STB - demonstrates that separator pressure is not an arbitrary operating choice. The $1.118/STB difference across a 10,000 STB/day field generates $4.08 million additional annual revenue from nothing more than operating the separator at the thermodynamically optimal pressure for that specific fluid composition. This optimization must be performed using equation-of-state flash calculations that account for the actual reservoir fluid composition, because the optimal pressure is specific to the C3-C6 content of the produced fluid - a gas condensate with high C4-C6 fraction has a very different optimal separator pressure than a black oil with high C7+.

The slug volume calculation - 105 bbl slugs arriving at peak rates of 10,080 bbl/day (12.6x the average 800 bbl/day design rate) - illustrates why separator sizing based on average production rate is inadequate for terrain-influenced gathering networks. A separator sized for 800 bbl/day would be overwhelmed in 6 minutes by a 105-bbl slug. The choice between slug catcher installation, flowline rerouting, and active slug control is an engineering-economic decision that requires knowledge of the slug volume, the slug frequency, and the cost of each mitigation option. The slug volume calculation is the starting point of that analysis, and it requires knowledge of the terrain profile of the flowline - which is why accurate survey data for the gathering system topography is as important as the production rate data when designing the separation and slug management infrastructure.

Want to access our gathering network design toolkit with multiphase pressure drop calculator, separator pressure optimizer, slug volume estimator, and compression power calculator, or discuss gathering system design for a specific field? Join our Telegram group for production engineering and flow assurance discussions, or visit our YouTube channel for step-by-step tutorials on multiphase flow modeling and gathering network optimization.