critical rate

Critical rate in petroleum engineering refers to the minimum flow rate at which the flow regime in a pipe, wellbore, or porous medium transitions from laminar (Darcy) flow to turbulent (non-Darcy) flow, a transition that occurs at a Reynolds number above approximately 2,100 in pipe flow (fully turbulent above Reynolds number 4,000) and at Forchheimer velocity above the critical interstitial velocity in porous media, with the onset of turbulence increasing pressure losses beyond what Darcy's law predicts and introducing an additional rate-dependent pressure drop term proportional to the square of the flow velocity; in the Western Canada Sedimentary Basin, the critical rate concept is applied in two distinct contexts: turbulent flow onset in high-rate WCSB Cardium, Viking, Montney, and Deep Basin gas wells (where non-Darcy flow in the near-wellbore critical matrix and within hydraulic fracture proppant packs imposes a rate-dependent turbulence skin that reduces gas well deliverability below Darcy-only predictions), and the Turner critical rate for liquid loading in declining WCSB gas wells (the minimum gas production rate at which gas velocity in the tubing string is sufficient to continuously lift produced water and condensate droplets to surface, below which liquid accumulation causes wellbore hydraulic loading that kills production). The Forchheimer equation extends Darcy's law to include turbulent pressure drop in porous media: the pressure gradient equals the sum of mu times v divided by k (Darcy viscous term) plus beta times rho times v squared (inertial or turbulent term), where mu is gas viscosity, v is interstitial velocity, k is permeability, beta is the Forchheimer turbulence coefficient (inversely proportional to permeability to the 1.5 power in most WCSB reservoir rock correlations), and rho is gas density; in a WCSB Cardium gas well producing at 50,000 to 200,000 m3/d through perforations in a 5 mD sandstone, the inertial term in the Forchheimer equation near the wellbore accounts for 20 to 50 percent of total near-wellbore pressure drop, equivalent to a positive non-Darcy turbulence skin of 2 to 8 that reduces gas productivity below Darcy theory. The rate-dependent nature of turbulence skin in WCSB gas wells is identified through multi-rate gas well tests (back-pressure tests or isochronal tests) in which the apparent skin varies with flow rate, with the turbulence coefficient D extracted from plotting apparent skin S' versus rate q (the slope of S' versus q equals D, the turbulence factor in units of 1/(103 m3/d)); AER Directive 040 (Pressure and Deliverability Testing) requires WCSB gas wells to conduct multi-rate deliverability tests that capture the turbulence contribution to the inflow performance relationship.

• Non-Darcy turbulent flow in WCSB high-rate gas well near-wellbore completions: Non-Darcy flow in the near-wellbore region of WCSB high-rate gas wells arises because the convergent radial geometry forces a progressively higher interstitial gas velocity as flow converges from the drainage area toward the wellbore, with the velocity proportional to 1/r; for a WCSB Cardium gas well producing 100,000 m3/d through a 0.1 m wellbore in a 3 m net pay with porosity 0.18, the interstitial velocity in the 0.1 to 0.2 m near-wellbore zone exceeds 15 m/s, well above the critical velocity for turbulent onset in the 3 to 5 mD WCSB Cardium sandstone (turbulent onset occurs at interstitial velocities of 1 to 3 m/s in medium-permeability sandstone). Perforation skin in WCSB perforated completions adds a concentrated turbulence effect: gas accelerating through perforation tunnel constrictions (6 to 12 mm diameter, 0.3 to 0.6 m length) at high flow rates generates localized turbulent pressure drop that can add 5 to 20 MPa of perforation friction over the 0.3 m perforation length, overwhelming the Darcy radial flow pressure drop in the reservoir at rates above 100,000 m3/d; WCSB operators minimize perforation turbulence by maximizing shot density (16 shots per meter versus 4 shots per meter standard), using large entrance hole diameter charges (16 to 20 mm versus 8 to 10 mm), and selecting perforation phasing and orientation to maximize the flow area through each perforation cluster.
• Proppant pack turbulence and non-Darcy flow in WCSB hydraulically fractured gas wells: Hydraulic fractures in WCSB Montney, Cardium, and Viking gas wells create high-permeability flow conduits (propped fracture permeability 50,000 to 500,000 mD for 40-mesh sand at 20 to 35 MPa closure) that dramatically increase gas deliverability over unfractured completions, but also concentrate flow through a narrow proppant pack where velocities are high enough to generate significant inertial pressure drop at WCSB gas production rates above 30,000 to 100,000 m3/d per stage. The Forchheimer beta coefficient for proppant packs is lower than for reservoir rock (typically 50,000 to 200,000 m-1 for 40-mesh sand versus 500,000 to 5,000,000 m-1 for 5 mD reservoir rock), but the much higher flow velocities in the propped fracture near the wellbore still produce non-Darcy pressure losses of 0.5 to 3 MPa per stage at WCSB Montney production rates of 50,000 to 200,000 m3/d. In WCSB Montney multistage horizontal well completion optimization, the non-Darcy fracture pressure loss is reduced by maximizing proppant pack width (using larger proppant volumes and slickwater treatments that create wider fractures) and by increasing the number of perforation clusters per stage (reducing the flow rate per cluster and therefore reducing velocity and turbulence within each fracture).
• Turner critical rate for liquid loading in WCSB declining shallow gas wells: The Turner critical rate (Turner, Hubbard, and Dukler, 1969) defines the minimum gas production rate needed to continuously lift liquid droplets (formation water or condensate) against gravity in a gas well tubing string, below which liquid accumulates in the wellbore and progressively kills the well; the critical velocity Vc = 6.6 times the fourth root of surface tension times (liquid density minus gas density) divided by gas density squared, expressed in m/s, and the corresponding critical flow rate Qc is Vc times the tubing cross-sectional area times the gas compressibility factor. In WCSB shallow gas wells (Medicine Hat, Horseshoe Canyon, Belly River) at late-life wellhead pressures of 0.5 to 2 MPa producing water and condensate through 60 mm (2-3/8 inch) or 73 mm (2-7/8 inch) tubing, the Turner critical rate ranges from 5,000 to 25,000 m3/d; when production falls below this critical rate, liquid loading begins within days to weeks, and the well requires intervention (plunger lift installation, velocity string running, or compression addition) to maintain economic rates. AER Directive 023 (Mandatory Liquid Unloading Reporting) requires WCSB operators to report wells with documented liquid loading and the remediation method applied, as chronic liquid loading in WCSB shallow gas pools is an economic and conservation concern where unloaded reserves may be stranded if the well kills before depletion.
• Multi-rate deliverability testing and LIT analysis for WCSB gas well turbulence quantification: Multi-rate deliverability testing of WCSB gas wells uses the laminar-inertial-turbulent (LIT) analysis method (Jones, Blount, and Glaze, 1976) to separately quantify the Darcy (laminar) and non-Darcy (inertial-turbulent) contributions to the wellbore pressure drop at each test rate: the stabilized flowing bottomhole pressure at each rate is used to calculate the pressure-squared drawdown divided by rate (P-squared form) or the pseudopressure drawdown divided by rate (pseudopressure form), which plots linearly versus rate for non-Darcy flow with slope equal to the turbulence parameter D and intercept equal to the Darcy deliverability coefficient. In WCSB Cardium and Viking gas pools where AER Directive 040 requires stabilized deliverability tests before pool production allocation, the LIT deliverability analysis gives the absolute open flow potential (AOFP) under both Darcy-only and turbulence-inclusive assumptions; the turbulence-inclusive AOFP (typically 20 to 40 percent lower than the Darcy-only value in high-permeability WCSB gas wells) is the operationally relevant figure for determining well productivity and production rate allocation under the pool production allocation scheme.
• Critical rate in WCSB steam injection and thermal production operations: In WCSB thermal heavy oil and oil sands production (CSS in Cold Lake and Lloydminster, SAGD in Athabasca and Cold Lake), the critical rate concept applies to the minimum steam injection rate required to maintain the steam chamber growth rate in SAGD operations and to the minimum bitumen production rate needed to prevent the SAGD production well from flooding with hot water that cools the wellbore and reduces bitumen mobility. In WCSB CSS operations, the critical injection rate is the minimum steam injection rate per cycle sufficient to heat the near-wellbore critical matrix to a temperature above the viscosity threshold for gravity drainage (bitumen viscosity below approximately 500 mPa-s at reservoir temperature), which for Athabasca bitumen (viscosity 1,000,000 mPa-s at 10 degrees Celsius) requires heating to 140 to 160 degrees Celsius through injection of steam at 8 to 12 MPa over 15 to 30 days; below the critical injection rate, the steam chamber does not reach sufficient areal extent to produce economic bitumen volumes during the subsequent production phase, reducing the steam-to-oil ratio and the economic viability of the CSS cycle.

Turner Critical Rate Remediation in WCSB Medicine Hat Gas Well

A WCSB Medicine Hat shallow gas operator observed production decline from 18,000 m3/d to 4,200 m3/d over 6 months on a well producing through 73 mm (2-7/8 inch) tubing at 1.1 MPa wellhead pressure, with erratic surface pressure indicating liquid slugging. Turner critical rate calculation for the well conditions (wellhead pressure 1.1 MPa, gas gravity 0.60, condensate surface tension 20 mN/m, condensate density 720 kg/m3) gave a critical rate of 12,500 m3/d; the well was producing well below this threshold. A 38 mm (1-1/2 inch) velocity string was run inside the existing tubing to reduce the flow area and increase gas velocity above the Turner critical rate at the measured 4,200 m3/d production rate. Post-velocity-string production stabilized at 6,800 m3/d (above the new critical rate of 3,200 m3/d in the smaller tubing), recovering 2,600 m3/d of lost production attributable to liquid loading. The velocity string installation cost $38,000 and returned payout in 4 months at current Medicine Hat gas prices.

Related Terms

Non-Darcy flow in WCSB high-rate gas wells introduces turbulent pressure losses proportional to velocity squared; the Forchheimer equation adds the inertial beta term to Darcy's law, producing rate-dependent turbulence skin in near-wellbore and proppant pack flow regimes. Liquid loading in WCSB declining gas wells begins below the Turner critical rate; water and condensate accumulate in the tubing string and progressively kill production, requiring plunger lift or velocity string intervention. Deliverability test under AER Directive 040 uses multi-rate LIT analysis to separate Darcy and non-Darcy pressure drop components and determine the turbulence-inclusive AOFP for WCSB gas well production allocation. Plunger lift is the most common WCSB shallow gas well remediation for production below the Turner critical rate; a free-traveling plunger cycles between the wellbore bottom and surface to periodically unload liquid accumulations. Inflow performance relationships for WCSB gas wells must incorporate the non-Darcy turbulence skin to correctly predict well deliverability at rates above the critical turbulent onset rate; Darcy-only IPR curves overestimate deliverability by 20-50% in high-permeability WCSB gas completions.