Gas Well Liquid Loading & Critical Velocity

The mathematical physics of liquid loading was pioneered by Turner, Hubbard, and Dukler (1969). They established that in vertical gas wells producing at low liquid-gas ratios (< 100 bbl/MMSCF), liquid is transported primarily as entrained liquid droplets suspended in a continuous upward gas core.

Derivation of the Force Balance

A falling spherical liquid droplet in an upward-flowing gas stream experiences two opposing vertical forces:

1. Downward Gravity Force (corrected for buoyancy):

   \[ F_g = \frac{\pi}{6} d_p^3 (\rho_L - \rho_g) g \]

2. Upward Aerodynamic Drag Force:

   \[ F_d = C_D \cdot \left(\frac{\pi d_p^2}{4}\right) \cdot \frac{1}{2} \rho_g v_g^2 \]

At the critical terminal velocity (\(v_{crit}\)), aerodynamic drag exactly balances gravitational weight (\(F_d = F_g\)), causing the droplet to remain stationary in the gas stream:

\[ v_{crit} = \sqrt{\frac{4 g d_p (\rho_L - \rho_g)}{3 C_D \rho_g}} \]

The Hinze Critical Weber Number & Droplet Shatter Limit

The droplet diameter (\(d_p\)) cannot grow indefinitely; if a droplet becomes too large, aerodynamic shear forces overcome the liquid surface tension (\(\sigma\)), shattering it into smaller stable droplets. Hinze (1955) demonstrated that droplet breakup is governed by the critical Weber number (\(We_{crit}\)):

\[ We_{crit} = \frac{\rho_g v^2 d_{max}}{\sigma} \approx 30 \]

Solving for maximum droplet diameter \(d_{max} = \frac{30 \sigma}{\rho_g v^2}\) and substituting back into the terminal velocity equation with an empirical drag coefficient \(C_D \approx 0.44\) yields the universal Turner droplet lift equation:

\[ v_{crit} = K \cdot \left[ \frac{\sigma \cdot (\rho_L - \rho_g)}{\rho_g^2} \right]^{0.25} \]

While Turner's formulation is universal, the empirical coefficient \(K\) has been refined based on extensive laboratory and field production data:

Critical Unloading Gas Flow Rate (\(Q_{crit}\))

The critical gas flow rate required to unload a well with tubing internal cross-sectional area \(A_{pipe}\) is:

\[ Q_{crit} = v_{crit} \cdot A_{pipe} \cdot \left(\frac{P_{wh}}{P_{sc}}\right) \left(\frac{T_{sc}}{T_{wh}}\right) \frac{1}{Z_{wh}} \quad [\text{Sm}^3\text{/day or MMSCFD}] \]

When actual gas flow rate falls below \(Q_{crit}\), the well enters the liquid loading death cycle: fallback droplets coalesce at the bottom → hydrostatic head increases → bottomhole flowing pressure \(P_{wf}\) rises → reservoir drawdown \((P_R - P_{wf})\) shrinks → production drops further → the well ceases flowing permanently.

Velocity String Retrofit Engineering

Because critical velocity \(v_{crit}\) is independent of pipe diameter, reducing the internal diameter of the production conduit reduces the cross-sectional flow area \(A_{pipe}\), thereby increasing in-situ gas velocity for the same volumetric production rate. Sizing a smaller velocity string (e.g. installing a 1.50" or 1.75" coiled tubing string inside existing 2.875" or 3.50" tubing) restores gas velocity above \(v_{crit}\), permanently unloading the well without requiring surface power.

Alternative Artificial Lift Techniques for Loaded Gas Wells:

• Plunger Lift: A free-traveling mechanical piston cycled between the well bottom and wellhead. The piston provides a physical solid seal that prevents gas slippage through the liquid column, lifting 100% of accumulated fluid slugs to surface using the well's own shut-in reservoir energy.
• Foaming Surfactant Injection: Capillary injection of chemical foaming agents into the rat-hole. Surfactants lower interfacial surface tension (\(\sigma\)) from 60 to ~20 dynes/cm and trap liquids in low-density foam, reducing critical velocity by over 50%.
• Downhole Gas Compressors & Pumping: Progressive cavity pumps (PCP) or electrical submersible pumps (ESP) installed to mechanically pump liquids up the annulus while gas flows up the tubing.