Gas Well Nodal Analysis & Tubing Performance Sizer (VLP vs. IPR)

In gas production engineering, Nodal Analysis divides the complete fluid flow path into two hydraulic components at a selected evaluation node—conventionally the bottomhole perforations at total vertical depth ($TVD$): 1. **Inflow Performance Relationship (IPR):** Represents fluid delivery capacity from the reservoir drainage boundary to the wellbore sandface. For gas wells, the classic Rawlins and Schellhardt (1936) backpressure formulation models laminar-turbulent reservoir flow: $$Q_g = C \cdot (P_R^2 - P_{wf}^2)^n$$ where $C$ is the stabilized deliverability coefficient ($\text{MMSCFD/psi}^{2n}$), $n$ is the turbulence exponent ($0.5 \le n \le 1.0$), $P_R$ is static average reservoir pressure, and $P_{wf}$ is flowing bottomhole pressure. The theoretical maximum deliverability at $P_{wf} = 0$ is the Absolute Open Flow Potential ($AOFP = C \cdot P_R^{2n}$). 2. **Vertical Lift Performance (VLP):** Represents the hydraulic pressure required to lift fluid against gravity and tubing friction from the sandface to the surface separator: $$P_{wf} = P_{th} + \Delta P_{\text{hydrostatic}} + \Delta P_{\text{friction}}$$ where $P_{th}$ is flowing tubing head pressure. The physical operating point occurs at the exact mathematical intersection where $P_{wf,\text{IPR}}(Q) = P_{wf,\text{VLP}}(Q)$. If no intersection exists, the reservoir energy is insufficient to overcome hydrostatic and frictional losses, and the well cannot flow naturally.

Even in nominally 'dry' gas wells, retrograde condensate dropout and formation water condensation generate a multi-phase flow regime inside the production tubing. The Beggs & Brill (1973) correlation determines the in-situ liquid holdup ($H_L$) and two-phase friction factor across inclination angles: 1. **Input Liquid Volume Fraction (No-Slip Holdup, $\lambda_L$):** $$\lambda_L = \frac{q_L}{q_L + q_g}$$ 2. **Flow Pattern Classification:** The mixture Froude number ($Fr_m = v_m^2 / (g D)$) and $\lambda_L$ determine whether the local flow pattern is Segregated, Intermittent (Slug), or Distributed (Mist). 3. **In-Situ Liquid Holdup Correction:** Horizontal holdup $H_L(0)$ is adjusted for upward vertical inclination ($\theta = 90^\circ$): $$H_L(\theta) = H_L(0) \cdot \left[1 + C \cdot (\sin(1.8 \theta) - 0.333 \sin^3(1.8 \theta))\right]$$ 4. **Effective Multiphase Density & Friction:** $$\rho_m = \rho_L H_L + \rho_g (1 - H_L)$$ $$\left(\frac{dP}{dz}\right)_{\text{total}} = \rho_m g + \frac{2 f_{tp} \rho_n v_m^2}{D_{ti}}$$ At low gas velocities, liquid slips backward, driving up $H_L$ and hydrostatic pressure loss. At high gas velocities, friction dominates.

Liquid loading is the progressive accumulation of co-produced water and hydrocarbon condensate in the wellbore, which occurs when superficial gas velocity falls below the terminal settling velocity of the largest liquid droplet. Turner et al. (1969), refined by Coleman et al. (1991) for low-pressure gas wells, modeled droplet suspension based on a balance between upward aerodynamic drag and downward gravitational force at the critical droplet Weber number ($We = 30$): $$v_{\text{crit}} = C_{\text{lift}} \cdot \left[\frac{\sigma \cdot (\rho_L - \rho_g)}{\rho_g^2}\right]^{0.25}$$ where: * $v_{\text{crit}}$ is critical terminal superficial gas velocity ($\text{m/s}$ or $\text{ft/s}$) * $\sigma$ is liquid interfacial tension ($60\text{ dynes/cm}$ for water, $20\text{ dynes/cm}$ for condensate) * $\rho_L, \rho_g$ are liquid and gas densities ($\text{kg/m}^3$) * $C_{\text{lift}}$ is $1.593$ in SI units ($3.06$ for Turner; $2.55$ for Coleman in US units) Converting $v_{\text{crit}}$ to gas volumetric rate at flowing conditions defines the Turner Critical Gas Flow Rate ($Q_{\text{crit}}$): $$Q_{\text{crit}} = 3.06 \cdot \frac{A_{\text{tubing}} \cdot P_{th}}{Z \cdot T_{th}} \cdot \left[\frac{\sigma \cdot (\rho_L - \rho_g)}{\rho_g^2}\right]^{0.25}$$ If actual production $Q_{\text{opt}} < Q_{\text{crit}}$, the well is actively loading liquids and requires artificial lift intervention (velocity string, plunger lift, foaming, or compression).