Bubble Flow in Vertical Production Tubing: Dispersed Gas Bubble Transport in Continuous Liquid Phase and Artificial Lift Optimization in WCSB Oil Wells
Bubble flow in multiphase wellbore hydraulics is the flow pattern that exists in a vertical or near-vertical production tubing string when gas is present in relatively small volumes as discrete, dispersed bubbles rising through a continuous liquid phase, with the bubbles distributed across the pipe cross-section and travelling upward at a velocity that is the sum of the bulk mixture velocity and the bubble slip velocity relative to the liquid, producing a flow regime in which the gas phase is entirely entrained and does not form a continuous gas channel along the tubing wall or coalesce into large gas slugs. The bubble flow regime occupies the low-gas-fraction end of the multiphase flow map for vertical pipes: it is bounded on the lower side by single-phase liquid flow (gas volume fraction approaching zero) and on the upper side by slug flow (where bubble coalescence creates intermittent large Taylor bubbles that occupy nearly the full pipe cross-section and alternate with liquid slugs). The gas volume fraction (void fraction) in the bubble flow regime for vertical production tubing typically ranges from approximately 0.01 to 0.25 depending on the flow pattern transition criteria applied, with the Duns-Ros, Griffith-Wallis, and Hagedorn-Brown correlations each defining the bubble-to-slug transition at somewhat different void fractions based on the tubing diameter, fluid properties, and superficial velocities. In WCSB oil production, the bubble flow regime is most commonly encountered in the upper portions of the production tubing string in medium-API-gravity wells producing close to or slightly below the bubble point pressure, where the gas that has evolved from solution in the reservoir arrives at the tubing intake as a minor gas phase co-mingled with the produced oil, and in the injection side of gas-lift installations where the injection gas volume rate has been set low enough that the gas is dispersed as small bubbles in the produced liquid column rather than forming slugs that would cause cyclic flow and pressure oscillations. The critical practical importance of correctly identifying bubble flow conditions in WCSB tubing design is that pressure traverse calculations in the bubble flow regime must account for gas slip (the tendency of small gas bubbles to rise faster than the liquid, increasing the actual in-situ gas void fraction above the no-slip value at a given gas-liquid ratio) to avoid underestimating the hydrostatic component of the flowing bottomhole pressure, which in turn affects artificial lift design, production rate optimization, and nodal analysis predictions for WCSB Cardium, Viking, and Devonian oil producers.
Key Takeaways
- Flow pattern transition criteria from bubble to slug flow in vertical WCSB production tubing and the conditions favoring each regime: The bubble-to-slug flow transition in vertical tubing occurs when bubbles coalesce faster than they are dispersed by turbulent breakup, typically when the gas void fraction exceeds a critical value (approximately 0.20-0.25 in the Griffith-Wallis model, 0.25 in the Duns-Ros model) or when the superficial gas velocity exceeds a diameter-dependent threshold. For a WCSB Cardium oil well producing through 60.3 mm (2.375-inch) production tubing with 30 m3/m3 GOR oil at 1,200 kPa wellhead pressure, the transition to slug flow typically occurs at gas-liquid ratios above approximately 15-25 standard m3/m3 at tubing flowing conditions, depending on the mixture velocity and fluid properties. Wells operating in bubble flow are hydraulically smoother and more amenable to steady-state inflow-outflow analysis; wells in slug flow exhibit intermittent pressure and flow-rate fluctuations that complicate nodal analysis and can cause artificial lift instability. Designing a gas-lift installation to keep injection near the bubble-to-slug transition (maximizing gas void fraction without triggering slug instability) maximizes lift efficiency while preserving flow stability.
- Bubble slip velocity and its effect on in-situ void fraction and hydrostatic pressure gradient in WCSB wellbore pressure traverse calculations: In bubble flow, individual gas bubbles rise faster than the bulk liquid-gas mixture because buoyancy drives them upward relative to the denser surrounding liquid. The slip velocity of small spherical bubbles in vertical pipes is described by the Harmathy equation: Vslip = 1.53 × (g × sigma × (rho_L - rho_G) / rho_L squared)^0.25, where sigma is the surface tension, g is gravitational acceleration, and rho are phase densities. For WCSB light oil (API 38, sigma approximately 0.020 N/m, rho_L approximately 800 kg/m3) at reservoir conditions, Vslip is approximately 0.15-0.25 m/s for small bubbles. This slip means the in-situ gas void fraction is higher than the no-slip (homogeneous) void fraction at any given gas-liquid ratio: at a superficial gas velocity of 0.3 m/s and superficial liquid velocity of 1.0 m/s, the no-slip void fraction is 0.23 but the in-situ void fraction with slip is approximately 0.29-0.33 depending on the drift-flux model parameters. Underestimating void fraction by 0.10 in the bubble flow regime underestimates the average mixture density by approximately 60-80 kg/m3, which translates to a flowing BHP underestimate of approximately 200-280 kPa per 1,000 m of tubing, an error large enough to incorrectly size an artificial lift installation or misdiagnose a Cardium well's actual reservoir pressure.
- Duns-Ros and Hagedorn-Brown pressure traverse methods for bubble flow in WCSB oil well nodal analysis and their accuracy limitations: Two of the most widely applied multiphase pressure traverse correlations for WCSB vertical oil well calculations are the Duns-Ros (1963) method, which uses dimensionless groups (liquid velocity number, gas velocity number, diameter number, and liquid viscosity number) calibrated against a range of pipe sizes and fluid properties, and the Hagedorn-Brown (1965) method, which uses empirical correlations for the liquid hold-up (1 minus void fraction) derived from field measurements on a small-diameter test well at various GOR and flow rate conditions. In the bubble flow regime specifically, Hagedorn-Brown tends to overpredict liquid holdup (underpredict void fraction) by 5-15% for WCSB medium-gravity oils compared to mechanistic drift-flux models, while Duns-Ros performs better at matching void fraction when properly evaluated in the correct flow regime. For detailed artificial lift design in WCSB Cardium and Viking infill horizontal wells where the tubing string transitions through multiple flow patterns from perforations to surface, mechanistic models (Hasan-Kabir, Ansari et al.) that explicitly solve the drift-flux equations in each flow regime are preferred over empirical correlations, particularly when GOR and water cut vary over the producing life of the well.
- Gas-lift design optimization using bubble flow regime targeting for WCSB medium-gravity oil production: In gas-lift design for WCSB Cardium or Viking medium-gravity oil wells (API 32-40, producing GOR 50-200 m3/m3, water cut 20-60%), the injected gas volume rate is selected to maximize the reduction in tubing hydrostatic gradient while avoiding slug flow instability. Operating in the upper range of the bubble flow regime (void fraction 0.15-0.22) provides a favourable balance: the in-situ mixture density is reduced by 20-35% from single-phase liquid density, meaningfully lowering the required flowing bottomhole pressure for a given surface production rate, while the dispersed bubble pattern prevents the intermittent pressure surges and production rate fluctuations characteristic of slug flow. Gas-lift mandrel spacing calculations for WCSB wells must account for the bubble-to-slug transition depth, which varies with temperature, pressure, and GOR along the tubing string: the injection point should be placed above the depth where slug flow would spontaneously initiate at the target injection rate, to ensure that gas disperses as bubbles throughout the column above the injection point rather than forming slugs that would reduce lift efficiency and cause surface separator surging.
- Bubble flow identification in WCSB production well monitoring using pressure and acoustic sensors and its implications for well performance diagnosis: Field identification of the flow regime in a WCSB production well's tubing string is valuable for diagnosing operational problems and optimizing artificial lift performance. Downhole pressure gauges recording flowing bottomhole pressure can detect the characteristic gradients associated with bubble flow: the pressure-depth gradient in bubble flow is nearly linear and intermediate between the single-phase liquid gradient and the single-phase gas gradient, while slug flow produces a stepped or oscillatory gradient profile that can be detected by fast-sampling (1-second) downhole gauges as slug-induced pressure waves passing the gauge depth. Acoustic flow meters at the surface (based on Doppler or cross-correlation of acoustic signals) and casing-annulus pressure oscillation monitoring provide indirect indicators of the flow regime in the tubing string: steady low-amplitude oscillations suggest bubble flow, while periodic high-amplitude pulses at 30-120 second intervals suggest slug flow. For WCSB Cardium gas-lift wells showing declining production efficiency over time, transitioning from bubble to slug flow as reservoir pressure declines and GOR increases is a common cause of artificial lift underperformance that can be corrected by reducing injection gas rate to move back into the stable bubble flow window.
Gas-Lift Rate Optimization to Maintain Bubble Flow in a Declining WCSB Cardium Producer
A Pembina Cardium gas-lift well (API 37 oil, 2.875-inch tubing, 1,600 m depth to perforations) is producing 85 m3/d oil with 40% water cut, GOR 95 m3/m3, at a wellhead injection rate of 18,000 standard m3/d lift gas. Surface separator shows intermittent 45-second surging cycles consistent with slug flow. A nodal analysis review confirms: at 18,000 m3/d injection, the calculated void fraction at mid-tubing (800 m) is 0.28, above the Griffith-Wallis slug transition of 0.25 for this tubing diameter. Reducing injection to 13,500 m3/d brings the mid-tubing void fraction to 0.21 (bubble flow). Post-adjustment: surging eliminated, production rate increases from 85 to 94 m3/d oil as the steadier inflow pressure profile allows the reservoir to drain at higher average productivity. Lift efficiency improves from 9.0 to 11.4 m3 oil per 1,000 m3 injection gas, reducing compression fuel costs by CAD 8,400 per month.
Fast Facts
The systematic classification of vertical multiphase flow into bubble, slug, churn, and annular regimes was developed primarily in the 1950s and 1960s through experimental work in transparent test pipes using air-water and air-oil mixtures. The petroleum industry adapted these laboratory flow maps to actual wellbore conditions (higher pressures, field-measured fluid properties, inclined and deviated tubing) through correlations developed by Duns and Ros in 1963 and Hagedorn and Brown in 1965, which remain in commercial nodal analysis software used for WCSB well design today.
Related Terms
The slug flow regime that develops when bubble coalescence exceeds dispersal, creating intermittent large gas Taylor bubbles alternating with liquid slugs in vertical production tubing, including its pressure oscillation signature and artificial lift instability implications for WCSB gas-lift and plunger-lift systems, is described under slug flow. The nodal analysis technique for matching inflow performance to tubing outflow performance, including the construction of tubing performance curves using multiphase pressure traverse correlations in bubble and slug flow regimes and their application to gas-lift rate optimization in WCSB Cardium oil wells, is described under nodal analysis. The gas-lift artificial lift method that injects compressed gas into the production tubing to reduce the hydrostatic head, including valve spacing, injection rate selection to target the bubble flow regime, and compressor sizing for WCSB medium-gravity oil fields, is described under gas lift.