ChemProCal
Knowledge Base + New Calculation Login Sign Up
ChemProCal • Thermodynamics • Gas Well Deliverability & IPR Curve
THERMODYNAMICS

Gas Well Deliverability & IPR Curve

Generate complete gas well Inflow Performance Relationship (IPR) curves and calculate Absolute Open Flow Potential (AOFP) using Rawlins-Schellhardt and Houpeurt LIT equations.



Deliverability Formulation

Reservoir Pressure

0.5 - 1.0

Well Test Point (Flowing Conditions)

Sm³/d

Operating Point Evaluation (Optional)

Sm³/d

📐 Engineering Formula & Specification Sheet

Gas Well Deliverability & IPR Curve: Sizing Methodology & Engineering Basis

Calculate gas well AOFP deliverability and generate Inflow Performance Relationship (IPR) curves using Rawlins-Schellhardt backpressure and Houpeurt LIT models.

👨‍🔬 Reviewed by: Process Design Engineer (B.Tech Chemical Engineering, 5+ Yrs Exp)
🗓️ Standards: ASME / API / GPSA / ISA
⏱️ Format: Quick Reference Spec Sheet

1. Rawlins & Schellhardt Backpressure Formulation

Derived by USBM (1935), the empirical gas deliverability equation accounts for high-velocity non-Darcy turbulence near the wellbore:

\[ q_g = C \cdot (P_R^2 - P_{wf}^2)^n \]

Where \(P_R\) is average reservoir pressure, \(P_{wf}\) is flowing bottomhole pressure, \(n\) is the turbulence exponent (\(0.5 \le n \le 1.0\)), and \(C\) is the stabilized deliverability coefficient.

2. Absolute Open Flow Potential (AOFP)

The theoretical maximum gas flow rate when flowing bottomhole pressure is pulled down to atmospheric pressure (\(P_{wf} = 0 \text{ barg}\)):

\[ \text{AOFP} = C \cdot (P_R^2 - P_{atm}^2)^n \]

📚 ChemProCal Knowledge Base

Need Complete Mathematical Derivations & Industrial Case Studies?

Explore our in-depth engineering textbook guide: Worked Example: Gas Well Deliverability & Inflow Performance Relationship (IPR) Modeling. Features step-by-step manual hand calculations, thermodynamic proofs, and commercial simulator benchmarks.

Read Full Theory Guide →
Peer Review & Engineering Accuracy

This specification sheet and its mathematical formulations adhere to internationally recognized process codes (ASME, API, GPSA, TEMA, and ISA). Sizing parameters and default design allowances are cross-referenced with empirical industry design practices.

About This Tool

The Gas Well Deliverability & Inflow Performance Relationship (IPR) Calculator determines the stabilized production capacity of natural gas reservoirs and generates complete bottomhole pressure-drawdown curves. Supporting both the empirical Rawlins-Schellhardt backpressure equation and the analytical Houpeurt Laminar-Inertial-Turbulent (LIT) model, this engine computes Absolute Open Flow Potential (AOFP) and evaluates gas well deliverability for field development planning and nodal analysis.

In petroleum and natural gas reservoir engineering, the Inflow Performance Relationship (IPR) describes the fluid flow from the reservoir drainage boundary to the wellbore face. Unlike liquid oil wells where Darcy's law produces a linear productivity index, gas wells suffer severe non-Darcy inertial turbulence and gas expansion at high velocities near the wellbore. This tool provides production engineers with verified calculations to calibrate multi-rate well tests and match reservoir deliverability against tubing intake curves.

1. Rawlins & Schellhardt Backpressure Formulation

Derived by Rawlins and Schellhardt in 1935 (USBM Monograph 7), the empirical backpressure deliverability equation remains the oil and gas industry's most widely recognized method for rating gas wells and reporting deliverability to government regulatory bodies:

\[ q_g = C \cdot \left( P_R^2 - P_{wf}^2 \right)^n \]

Where:

  • \(q_g\): Stabilized natural gas flow rate in \(\text{Sm}^3\text{/day}\) or \(\text{Mscfd}\).
  • \(P_R\): Average static reservoir shut-in pressure in \(\text{bar(a)}\) or \(\text{psia}\).
  • \(P_{wf}\): Flowing bottomhole pressure at the sandface in \(\text{bar(a)}\) or \(\text{psia}\).
  • \(C\): Stabilized well deliverability performance coefficient, accounting for rock permeability, reservoir thickness, gas viscosity, drainage radius, and mechanical skin.
  • \(n\): Backpressure turbulence exponent (dimensionless), ranging strictly from 0.5 to 1.0.

Physical Meaning of the Turbulence Exponent (\(n\)):

\(n = 1.0\) (Darcy Flow)

Purely laminar Darcy viscous flow. Negligible inertial resistance. Common only in very low permeability or tight gas formations at low production rates.

\(0.65 \le n \le 0.85\) (Typical)

Mixed transition regime. Darcy flow deep in the reservoir transitioning into non-Darcy inertial flow as gas accelerates converging toward the wellbore.

\(n = 0.5\) (Turbulent)

Completely turbulent non-Darcy flow. High-velocity inertial forces dominate near-wellbore perforations and gravel packs. Maximum possible flow resistance.

2. Houpeurt LIT Formulation: Separating Laminar from Inertial Skin

While the Rawlins-Schellhardt equation is empirical, Houpeurt (1959) derived the rigorous analytical Laminar-Inertial-Turbulent (LIT) equation directly from the Forchheimer equation for high-velocity fluid flow in porous media:

\[ P_R^2 - P_{wf}^2 = a \cdot q_g + b \cdot q_g^2 \]

Where:

  • Laminar Coefficient (\(a\)): Represents classical Darcy viscous pressure drop through the porous rock matrix:

    \[ a = \frac{1422 \cdot \bar{\mu}_g \cdot \bar{Z} \cdot T}{k \cdot h} \left[ \ln\left(\frac{r_e}{r_w}\right) - 0.75 + s \right] \]

    Notice that the mechanical damage skin factor (\(s\)) is contained entirely within the laminar term \(a\).
  • Non-Darcy Inertial Coefficient (\(b\)): Represents velocity-squared inertial drag and rate-dependent skin (\(D \cdot q_g\)) caused by microscopic flow convergence through pore throats:

    \[ b = \frac{1422 \cdot \bar{\mu}_g \cdot \bar{Z} \cdot T}{k \cdot h} \cdot D \]

    Where \(D\) is the non-Darcy flow factor proportional to the Geertsma / Jones / Cooke rock turbulence coefficient (\(\beta\)).

Solving Deliverability for Flow Rate:

Because Houpeurt's equation is quadratic in flow rate \(q_g\), deliverability at any specified flowing bottomhole pressure \(P_{wf}\) is solved analytically using the quadratic formula:

\[ q_g = \frac{-a + \sqrt{a^2 + 4 \cdot b \cdot (P_R^2 - P_{wf}^2)}}{2b} \]

3. Absolute Open Flow Potential (AOFP) & Well Nodal Analysis

The Absolute Open Flow Potential (AOFP) is the theoretical maximum flow rate a gas well could deliver if the sandface flowing bottomhole pressure could be reduced completely to atmospheric pressure (\(P_{wf} = 0 \text{ barg}\)):

\[ \text{AOFP} = C \cdot \left( P_R^2 - P_{atm}^2 \right)^n \approx C \cdot (P_R^2)^n \]

AOFP is the universal regulatory benchmark used by government bodies to assign well allowable production quotas. However, no gas well ever produces at its AOFP in practice because:

  1. Hydrostatic fluid column and frictional pressure losses in the production tubing string require a significant backpressure (\(P_{wf} \gg 0\)) to lift gas to surface.
  2. Producing at excessive drawdown induces formation sand failure, fines migration, casing collapse, and catastrophic water or gas coning.
  3. Industry best practice typically chokes gas wells to produce at 50% to 65% of AOFP for long-term reservoir recovery optimization.

Frequently Asked Questions

What is the difference between a flow-after-flow test and a modified isochronal test?

In a flow-after-flow test, the well is produced at four successive flow rates, each run until stabilized reservoir pressure equilibrium is reached. In low-permeability reservoirs, reaching true stabilization at each step can require weeks, wasting large volumes of flared gas. The modified isochronal test solves this by producing the well at four equal short time periods (e.g. 4 hours each) separated by shut-in periods of equal duration, followed by one single extended flow test to stabilization. This enables accurate determination of turbulence exponent \(n\) and stabilized coefficient \(C\) in tight gas formations within a fraction of the testing time.

Why should gas wells rarely produce at more than 60% to 70% of their AOFP?

Operating a gas well at maximum drawdown (> 70% of AOFP) creates extreme near-wellbore pressure gradients and flow velocities. This high velocity leads to severe non-Darcy inertial pressure drops (the \(b \cdot q^2\) effect), causing premature reservoir energy dissipation. Furthermore, excessive drawdown triggers formation sand collapse, perforations erosion, and causes bottom-water to cone rapidly upwards into the perforated interval, permanently drowning the well.

When should reservoir engineers use Pseudo-Pressure m(P) instead of Pressure-Squared P²?

The pressure-squared formulation (\(P^2\)) assumes that the product of gas viscosity and compressibility factor (\(\mu_g \cdot Z\)) is constant. This assumption is highly accurate for low to moderate reservoir pressures (\(P < 2,000 \text{ psia} / 140 \text{ bar}\)). However, in high-pressure, high-temperature (HPHT) reservoirs (> 3,000 psia), \(\mu_g \cdot Z\) varies dramatically with pressure. In deep HPHT reservoirs, engineers must integrate Al-Hussainy and Ramey's Real Gas Pseudo-Pressure \(m(P) = 2 \int (P / \mu_g Z) dP\) to maintain mathematical rigor.

What causes the deliverability coefficient C to decline over time during production?

The performance coefficient \(C\) is directly proportional to reservoir permeability-thickness (\(k \cdot h\)) and inversely proportional to gas viscosity and drainage radius \(\ln(r_e / r_w)\). Over years of production, \(C\) declines primarily due to: (1) permeability impairment (skin damage) from scale, condensate dropout, or fines migration near the wellbore, (2) two-phase liquid loading where water or condensate builds up in the pore space, reducing relative permeability to gas (\(k_{rg}\)), and (3) shrinking reservoir boundary drainage radius as offset infill wells are drilled.

Tool Units

Input / Output Units

Configure global unit system preset or individual parameter dimensions.