Pipe Flow Overview

Introduction

Pipe flow calculations determine pressure changes along wellbores, flowlines, and pipelines. These calculations are essential for:

  • Tubing performance (VLP) — pressure profile from bottomhole to wellhead
  • Flowline sizing — pressure loss in surface lines
  • Artificial lift design — gas lift, ESP, rod pump optimization
  • Production system analysis — nodal analysis, system optimization

Flow Types

Single-Phase Flow

When only one phase (liquid or gas) flows through the pipe:

Flow Type Examples Correlations Used
Liquid Water injection, dead oil Fanning equation
Gas Dry gas wells, gas pipelines Compressible flow equations

📖 Documentation: Single-Phase Pipe Flow

Multiphase Flow

When oil and gas (and possibly water) flow together, the flow becomes much more complex:

  • Slippage — gas moves faster than liquid
  • Flow patterns — bubble, slug, churn, annular, mist
  • Liquid holdup — fraction of pipe occupied by liquid

Empirical correlations are required because theoretical solutions are impractical.


Correlation Selection Guide

Decision Framework

Loading diagram...

Correlation Comparison

Correlation Inclination Flow Type Best Application
Beggs-Brill Any angle Oil + gas General purpose, horizontal/inclined
Hagedorn-Brown Vertical only Oil + gas Vertical oil wells, tubing
Gray Any angle Gas + liquid Gas wells with condensate/water
Fanning (liquid) Any angle Liquid only Water injection, dead oil
Compressible Any angle Gas only Dry gas wells, gas pipelines

Available Functions

Single-Phase Liquid (Incompressible)

Function Description
ReynoldsNumberLiquid Reynolds number for liquid flow
FrictionPressureDropLiquid Fanning friction pressure loss
PotentialEnergyPressureDropLiquid Elevation pressure change
InletPipePressureLiquid Calculate inlet from outlet
OutletPipePressureLiquid Calculate outlet from inlet

📖 Documentation: Single-Phase Pipe Flow


Single-Phase Gas (Compressible)

Function Description
ReynoldsNumberGas Reynolds number for gas flow
InletPipePressureGas Calculate inlet from outlet
OutletPipePressureGas Calculate outlet from inlet

📖 Documentation: Single-Phase Pipe Flow


Multiphase: Beggs-Brill (1973)

The most versatile correlation — applicable at any pipe inclination from horizontal to vertical, for both upward and downward flow.

Function Description
PressureGradientBeggsBrill Pressure gradient, psi/ft
InletPressureBeggsBrill Calculate inlet from outlet
OutletPressureBeggsBrill Calculate outlet from inlet

Key Features:

  • Flow pattern prediction (segregated, intermittent, distributed)
  • Liquid holdup correlation for each pattern
  • Inclination correction factor
  • Friction factor modification for multiphase

Best For: Horizontal flowlines, inclined wells, general-purpose calculations


Multiphase: Hagedorn-Brown (1965)

Developed specifically for vertical upward flow in oil wells. Includes Griffith bubble flow modification.

Function Description
PressureGradientHarBrown Pressure gradient, psi/ft
InletPressureHarBrown Calculate inlet from outlet
OutletPressureHarBrown Calculate outlet from inlet

Key Features:

  • Empirical holdup correlations from extensive test data
  • Griffith modification for bubble flow at low gas rates
  • No-slip density calculation

Best For: Vertical oil producers, tubing performance curves

Limitations: Only valid for vertical (90°) upward flow


Multiphase: Gray (1974)

Developed for gas wells producing liquids (condensate or water). Part of API 14B.

Function Description
PressureGradientGray Pressure gradient, psi/ft
InletPressureGray Calculate inlet from outlet
OutletPressureGray Calculate outlet from inlet

Key Features:

  • Designed for high gas-liquid ratios (GLR > 5000 scf/STB)
  • Accounts for liquid loading effects
  • Applicable at any inclination

Best For: Gas wells with condensate, wet gas wells, high-GLR producers


Input Parameters

Common Parameters

Parameter Symbol Units Description
Liquid rate QLQ_L bbl/d Total liquid (oil + water)
Gas rate QgQ_g mmscf/d Gas at standard conditions
Liquid density ρL\rho_L lb/ft³ At flowing conditions
Gas specific gravity γg\gamma_g - Air = 1.0
Liquid viscosity μL\mu_L cP At flowing conditions
Gas viscosity μg\mu_g cP At flowing conditions
Z-factor ZZ - Gas compressibility factor
IFT σ\sigma dynes/cm Gas-liquid interfacial tension

Pipe Parameters

Parameter Symbol Units Typical Range
Inner diameter dd in 2 - 12 (tubing), 4 - 24 (flowlines)
Length LL ft 100 - 20,000
Roughness ε/d\varepsilon/d - 0.0001 - 0.01
Angle θ\theta degrees -90 to +90

Angle Convention

Angle Description Example
Horizontal Surface flowline
+90° Vertical upward Producer wellbore
-90° Vertical downward Injector wellbore
+45° Inclined upward Deviated producer

Workflow: Tubing Performance Curve (VLP)

Step 1: Gather Data

  • Tubing: ID, length, roughness, well trajectory
  • Fluids: PVT data or correlations
  • Conditions: wellhead pressure, temperature profile

Step 2: Select Correlation

Well Type Recommended Correlation
Vertical oil well Hagedorn-Brown
Deviated oil well Beggs-Brill
Gas well with liquids Gray
Horizontal flowline Beggs-Brill

Step 3: Calculate Pressure Profile

For inlet pressure calculation (bottom to top):

P_inlet = InletPressureBeggsBrill(Ql, Rho_l, Ul, Qg, SGgas, IFT,
                                   pipeID, pipeLength, pipeRoughness,
                                   pipeAngle, P_wellhead, T_avg)

Step 4: Generate VLP Curve

For multiple flow rates:

For each rate Q:
  P_bhf = InletPressure(..., Q, ..., P_whp)
  Plot (Q, P_bhf)

Step 5: Find Operating Point

Intersect VLP curve with IPR curve to determine:

  • Operating flow rate
  • Required bottomhole flowing pressure

Calculation Tips

Pressure Segmentation

For better accuracy with long pipes or large pressure changes:

  1. Divide pipe into segments (e.g., 500 ft each)
  2. Calculate properties at average conditions for each segment
  3. March from known pressure to unknown

Temperature Effects

  • Use average temperature for short pipes
  • Use temperature profile for deep wells
  • Gas properties (Z, μg\mu_g) are temperature-sensitive

Critical Flow Check

Near critical flow conditions (approaching sonic velocity), standard correlations may be inaccurate. Check if:

vmixture<0.5vsonicv_{mixture} < 0.5 \cdot v_{sonic}


Troubleshooting

Problem Likely Cause Solution
Unrealistic pressure drop Wrong correlation for inclination Match correlation to geometry
Negative outlet pressure Flow rate too high for pipe Check pipe size, reduce rate
Results differ from field PVT data mismatch Verify fluid properties
Convergence issues Extreme conditions Use smaller segments

Detailed Correlations

Fluid Properties

Well Performance


References

  1. Beggs, H.D. and Brill, J.P. (1973). "A Study of Two-Phase Flow in Inclined Pipes." Journal of Petroleum Technology, May 1973, pp. 607-617. SPE-4007-PA.

  2. Hagedorn, A.R. and Brown, K.E. (1965). "Experimental Study of Pressure Gradients Occurring During Continuous Two-Phase Flow in Small-Diameter Vertical Conduits." Journal of Petroleum Technology, April 1965, pp. 475-484. SPE-940-PA.

  3. Gray, H.E. (1974). "Vertical Flow Correlation in Gas Wells." User's Manual for API 14B, Appendix B.

  4. Brill, J.P. and Mukherjee, H. (1999). Multiphase Flow in Wells. SPE Monograph Vol. 17.

  5. Brown, K.E. (1984). The Technology of Artificial Lift Methods, Vol. 1. PennWell Books.

An unhandled error has occurred. Reload X