Annular Flow: Multiphase Flow Regime, Liquid Loading, and Gas Well Design
Annular flow is a multiphase flow regime in which the lighter fluid (typically gas) occupies the central core of a pipe or wellbore at high velocity while the heavier fluid (typically liquid) exists simultaneously as a continuous film coating the inner pipe wall and as dispersed droplets entrained in the gas core. The liquid film is kept in motion by the interfacial drag of the high-velocity gas core, which exerts a shear stress on the liquid surface and drives liquid upward in a vertical well or along the pipe wall in a horizontal line. Annular flow is distinguished from other multiphase regimes by the continuity of both the gas and the liquid phases: unlike slug flow (where liquid and gas alternate as distinct plugs occupying the full pipe cross-section) or bubble flow (where small bubbles are dispersed in a continuous liquid phase), in annular flow the gas core is continuous from the source to the outlet and the liquid film is continuous around the pipe circumference. The transition to annular flow from slug flow occurs when gas velocity is high enough relative to liquid velocity that the Kelvin-Helmholtz instability that generates slug waves is suppressed and the liquid is forced to the pipe wall. In oil and gas production, annular flow is the dominant flow regime in vertical wellbores producing at high gas-to-liquid ratios (above approximately 1,000 to 5,000 standard cubic metres per cubic metre of liquid), in gas-condensate wells above the critical liquid loading velocity, and in gas gathering pipelines downstream of high-GOR wellheads. Understanding annular flow is essential for predicting wellbore pressure gradients, calculating bottom-hole flowing pressure from wellhead measurement, designing artificial lift systems, diagnosing liquid loading in declining gas wells, optimizing velocity string completions, and designing surface separator inlet conditions that account for the liquid film distribution arriving in the separator feed line.
Key Takeaways
- Flow regime transitions and the flow pattern map: The transition from slug flow to annular flow in a vertical pipe is mapped on a flow pattern diagram (also called a flow regime map) that plots superficial gas velocity (volumetric gas flow rate divided by total pipe cross-sectional area) versus superficial liquid velocity on logarithmic axes, with boundaries between flow regime regions determined empirically from laboratory and field measurements. The Taitel-Dukler and Barnea models are the most widely used mechanistic flow pattern maps for vertical and inclined pipes in petroleum engineering, defining the annular flow region as occurring when the superficial gas velocity exceeds approximately 10 to 25 metres per second at typical production wellbore pressures and temperatures. In the transition zone between slug and annular flow, an unstable churn flow regime exists where the liquid film is periodically disrupted by gas slugs, creating erratic surface flow rates and separator level fluctuations that complicate production measurement. WCSB tight gas wells in the Horseshoe Canyon formation at 700 to 1,200 metres depth typically operate in the annular flow regime early in their production life (high-rate period) but transition through churn and slug flow regimes as reservoir pressure declines and gas velocity decreases.
- Turner critical velocity and liquid loading: The Turner critical velocity (named after Hubert Turner's 1969 paper establishing the model) is the minimum gas velocity required to lift liquid droplets entrained in the gas core against gravity in a vertical wellbore, preventing their fallback and accumulation at the well bottom. Below the Turner critical velocity, liquid droplets entrained in the gas core are not efficiently transported to surface: they fall back through the rising gas and accumulate in the wellbore, increasing bottom-hole pressure, reducing gas production rate, and eventually killing the well if the accumulated liquid column grows tall enough to exceed the reservoir driving pressure. The Turner velocity model gives VT = 5.62 x ((σ x (ρL minus ρG) / ρG2)0.25, where sigma is gas-liquid surface tension in dynes/cm, rho_L and rho_G are liquid and gas densities in lbm/ft3. For a typical Cretaceous shallow gas well in Alberta producing sweet gas at 1,500 kPa wellhead pressure with condensate at 10 barrels per MMscf, Turner velocity is approximately 4.5 to 6.5 metres per second in a 3.5-inch tubing, corresponding to a minimum production rate of approximately 150 to 250 Mscf/d to sustain annular flow and prevent liquid loading.
- Annular flow pressure gradient and wellbore multiphase modeling: The pressure gradient in a producing well operating in the annular flow regime differs significantly from both single-phase gas flow and slug flow, because the liquid holdup (fraction of the pipe cross-section occupied by liquid at any given point) affects both the mixture density and the friction factor applied to the flowing stream. In annular flow, liquid holdup is low (typically 2 to 15 percent by volume) but concentrated in the film at the pipe wall where friction is highest, producing wall friction pressure losses that are intermediate between single-phase gas and slug flow values. The Hagedorn-Brown correlation and the Beggs-Brill correlation are the two most widely used empirical multiphase pressure gradient models for annular flow in production wellbores, both requiring input of wellbore geometry (tubing OD and ID), gas and liquid rates, PVT properties at each depth increment, and inclination angle. Bottom-hole flowing pressure (BHFP) calculated from wellhead pressure and production rate using these correlations is used to construct an inflow performance relationship that, combined with the tubing performance curve, identifies the natural flow rate at which reservoir inflow and wellbore deliverability balance for the given completion design.
- Velocity string completions to sustain annular flow: When a declining gas well's production rate falls below the Turner critical velocity in the production tubing string, the operator can restore annular flow and prevent liquid loading by installing a smaller-diameter velocity string inside the existing tubing. The smaller velocity string has a reduced cross-sectional area that increases gas velocity for the same volumetric flow rate: if Turner velocity is 5 m/s in the original 3.5-inch tubing at 200 Mscf/d, installing a 1.9-inch velocity string reduces the flow area by 70 percent and increases gas velocity to 16 m/s at the same rate, far above the Turner threshold and sufficient to re-establish annular flow and unload the accumulated liquid column. Velocity string installation in Alberta shallow Horseshoe Canyon and Pekisko gas wells typically costs CAD 28,000 to CAD 55,000 for coiled tubing installation of 1,900 metres of 1-1/4 or 1-1/2 inch velocity tubing, with production response observable within 24 to 72 hours of installation as the accumulated liquid is swept to surface in the annular flow regime restored by the higher gas velocity. A successful velocity string installation can extend the producing life of a liquid-loaded gas well by 3 to 8 years compared to abandonment at the natural liquid loading threshold.
- Annular flow in horizontal and near-horizontal wellbores: In horizontal and deviated wellbores, the annular flow regime takes on a different character than in vertical wells because gravity acts perpendicular to the flow direction rather than opposing it. In horizontal pipes, gravity segregates the liquid film to the bottom of the pipe rather than distributing it uniformly around the circumference, creating a stratified-annular or wavy-annular regime where the liquid phase is concentrated at the pipe bottom and the gas phase occupies the upper portion of the cross-section. This asymmetric liquid distribution affects pressure gradient calculations, increases the risk of liquid accumulation at low points in directional profiles, and complicates the design of artificial lift for horizontal well sections. In Montney and Duvernay horizontal wells producing above the dew point at high gas rates, annular flow exists in the vertical portion of the wellbore and transitions through stratified-wavy flow in the horizontal section, creating a complex flow assurance challenge at the heel of the well where the vertical and horizontal sections meet. Downhole flow meters placed in horizontal Montney wells confirm that individual perforation clusters contribute gas flow in proportion to their reservoir contact, with inter-cluster liquid re-distribution in the stratified annular flow regime creating commingled flow profiles that challenge individual stage production allocation.
Liquid Loading Diagnosis and Remediation in Alberta Shallow Gas Wells
Liquid loading is the most common artificial lift problem in shallow gas wells in the WCSB and represents the primary production optimization challenge for the approximately 65,000 shallow gas wells in Alberta, many of which produce from Horseshoe Canyon coal seams and associated Cretaceous sands at depths of 300 to 1,200 metres. Liquid loading develops gradually as reservoir pressure declines over the first 3 to 8 years of production and gas velocity in the tubing falls below the Turner critical velocity for the produced water and condensate. The early warning signs are visible in the wellhead surface data: increasing wellhead flowing pressure as liquids accumulate against the bottomhole pressure, erratic production rates with periodic large liquid slugs at the separator, declining instantaneous gas rate while cumulative production flattens, and eventually complete well death when the accumulated liquid head equals the reservoir pressure.
Diagnosing liquid loading versus other production decline mechanisms (reservoir depletion, formation damage, equipment failure) requires careful analysis of the wellhead flowing pressure and production rate history combined with a nodal analysis of the well's current flow regime. The key diagnostic question is whether the well's current gas velocity exceeds the Turner critical velocity: if the gas rate at current wellhead conditions places the well in the annular flow region of the flow pattern map, then the decline is probably reservoir depletion rather than liquid loading. If the well's operating point has moved below the Turner velocity and into the slug or churn flow region, liquid loading is confirmed and artificial lift or wellbore unloading is required. Field confirmation of liquid loading uses a flow gradient survey with a dead-weight pressure recorder run on slickline, which shows the pressure gradient increasing below the point of liquid accumulation compared to a pure gas gradient, confirming a liquid column in the wellbore.
Remediation options for liquid-loaded Horseshoe Canyon wells progress from low-cost to high-cost interventions in order of severity. Intermittent cycling (shutting in the well to allow pressure buildup, then opening rapidly to surge the accumulated liquid to surface in a high-velocity annular flow event) costs nothing except the deferred production during shut-in periods of 4 to 12 hours and can sustain production for 1 to 3 years beyond the natural loading threshold if the reservoir has sufficient pressure to sustain the surge velocity. Chemical treatment with foaming agents (surfactant injection at 100 to 500 mL per day down the casing-tubing annulus, falling into the wellbore and foaming the accumulated liquid into a lighter mixture that gas velocity can transport at a lower Turner velocity) costs CAD 800 to CAD 2,400 per year in chemical and adds 6 to 18 months of producing life per treated well. Velocity string installation at CAD 28,000 to CAD 55,000 per well is the next cost increment, followed by plunger lift installation (CAD 5,000 to CAD 15,000 for the plunger, lubricator, and controller) for wells with sufficient casing-to-tubing pressure differential to cycle a mechanical plunger that sweeps liquid slugs to surface on each upstroke.
In the Hanna area gas fields of southeast Alberta, where Horseshoe Canyon and Belly River gas wells at 400 to 900 metres depth are operated by small independents and royalty trusts, a systematic liquid loading surveillance program covering 420 wells identified 87 wells producing below their Turner velocity by applying monthly wellhead pressure and gas rate data to automated Turner velocity calculations in a production database. Of these 87 wells, 34 were treated with intermittent cycling, 28 with chemical foaming, 19 with velocity string installation, and 6 with plunger lift, at a total remediation cost of CAD 1.85 million. Total incremental production recovered from the 87 wells over the 24 months following remediation was approximately 280 MMscf of natural gas, valued at CAD 1.12 per GJ for ATC shallow gas reference price, generating CAD 5.6 million in gross revenue against the CAD 1.85 million remediation investment and confirming that systematic liquid loading remediation targeting wells operating below their annular flow Turner velocity threshold is one of the highest-return production optimization activities available to shallow gas operators in the WCSB.