Bond Number: When Gravity Beats Capillary Pressure in WCSB Waterflood and EOR Operations
The Bond number (Bo) is a dimensionless parameter in fluid mechanics and pore-scale physics that quantifies the relative importance of gravitational forces versus interfacial tension (capillary) forces acting on a fluid interface in a porous medium. The Bond number is defined as Bo = (Δρ × g × L²) / σ, where Δρ is the density difference between the two fluid phases (kg/m³), g is gravitational acceleration (9.81 m/s²), L is a characteristic length scale (typically the pore radius or tube diameter, in metres), and σ is the interfacial tension between the phases (N/m or mJ/m²). A Bond number much less than 1 (Bo ≪ 1) indicates capillary dominance: interfacial tension forces are large relative to buoyancy, and the non-wetting phase is trapped as residual saturation in pore throats regardless of the density contrast. A Bond number much greater than 1 (Bo ≫ 1) indicates gravitational dominance: buoyancy forces are large enough to drive phase separation and fluid migration against capillary resistance. In reservoir engineering practice, the Bond number is most directly useful in EOR (enhanced oil recovery) design, where the goal is to mobilize residual oil that waterflooding has left trapped in pore throats by capillary snap-off. Surfactant flooding reduces σ from the typical oil-water value of 20-35 mN/m to ultra-low values of 0.001-0.01 mN/m, increasing the Bond number by 3-4 orders of magnitude and allowing buoyancy and viscous forces to displace oil that would otherwise remain as residual saturation in the pore space. For a typical WCSB Cardium sandstone pore radius of 5 micrometres with crude oil density 850 kg/m³, water density 1,010 kg/m³, and σ = 25 mN/m, the Bond number is approximately 1.5 × 10⁻⁷ — deeply in the capillary-dominated regime, confirming that residual oil after waterflooding will remain trapped regardless of further water injection without chemical EOR intervention. Surfactant injection reducing σ to 0.005 mN/m raises the Bond number to 7.5 × 10⁻⁴ — still capillary-dominated, but approaching the critical Bond number for oil mobilization (Bo_critical ≈ 10⁻⁴ to 10⁻³ depending on pore geometry), which explains why surfactant EOR in tight carbonate and sandstone formations requires ultra-low interfacial tension to achieve significant incremental oil recovery above conventional waterflooding.
Key Takeaways
- Bond number in CO2 gravity override during WCSB EOR: In CO2 miscible or immiscible flooding of WCSB Devonian carbonate and Cardium sandstone reservoirs, CO2 (density 400-700 kg/m³ at reservoir conditions) is significantly less dense than reservoir brine (1,000-1,050 kg/m³), generating a large Δρ that drives CO2 to migrate upward through the reservoir by buoyancy — the "gravity override" problem that reduces CO2 sweep efficiency. The Bond number for CO2 in a macro-scale fracture (L = 1 cm) is approximately 0.2-0.6 — approaching gravitational dominance, explaining why CO2 gravity override is severe in fractured reservoirs but moderate in tight matrix rock. Foam injection reduces CO2 effective mobility and partially counteracts gravity override by creating a high-viscosity, density-modified CO2 foam bank that overrides less severely than free CO2.
- Capillary number vs Bond number in recovery factor analysis: The Bond number is often confused with the capillary number (Ca = μv/σ, where μ is fluid viscosity and v is Darcy velocity) but addresses a different force competition: Bond number is gravity versus capillary tension; capillary number is viscous versus capillary tension. Both must be considered in EOR design because residual oil trapping is controlled by the "total trapping number" Nt = sqrt(Ca² + 2Ca×Bo×cos(θ) + Bo²) — the vector sum of viscous and gravitational forces acting on a pore throat relative to capillary force. For a WCSB Cardium waterflood at typical injection velocity (1 m/day Darcy velocity), Ca ≈ 5 × 10⁻⁶ and Bo ≈ 10⁻⁷, so Nt ≈ Ca (viscous forces dominate over gravity at the pore scale), and the only effective EOR strategy is reducing σ by chemical means rather than relying on gravity.
- Bond number in SAGD gravity drainage physics: SAGD (steam-assisted gravity drainage) is a thermal recovery process where the Bond number is deliberately exploited at the macroscopic scale (L = reservoir thickness, 15-40 m in Athabasca oil sands). With Δρ ≈ 450 kg/m³ (bitumen vs steam condensate), L = 20 m, and σ reduced to approximately 0.5-2 mN/m at SAGD temperatures (250-300°C), Bo is in the range of 20-200 — strongly gravitational-dominance regime. This is why SAGD works: gravity drainage of mobilized bitumen is the dominant transport mechanism at the steam chamber interface, not viscous pressure-driven flow. At reservoir temperatures, the Bond number confirms that bitumen will drain by gravity regardless of the specific petrophysical heterogeneity of the oil sand matrix, as long as the steam chamber provides sufficient thermal mobilization to reduce bitumen viscosity below 50 mPa-s — the viscosity threshold below which gravity drainage rates become economically meaningful.
- Bond number effects on waterflood residual oil saturation (Sor): The relationship between Bond number and residual oil saturation for WCSB Cardium and Viking sandstones has been measured in core flooding experiments: at typical waterflood Bond numbers (Bo ≈ 10⁻⁷), Sor is approximately 0.20-0.35 (20-35% of pore volume remains as trapped oil after waterflood). Raising Bo to 10⁻³ by surfactant injection (ultra-low IFT) reduces Sor to approximately 0.02-0.05. This Sor reduction is the incremental recovery target for WCSB chemical EOR projects: if the flooded oil volume in a Cardium pool is 50 million barrels at Sor 0.25, a surfactant flood that reduces Sor to 0.05 could recover an additional 20 million barrels — the economic case that drives interest in chemical EOR pilot programs in mature WCSB Cardium and Viking pools despite the high surfactant cost (CAD 8,000-20,000/tonne) and injection logistics challenges.
- Critical Bond number and pore geometry effects: The critical Bond number at which gravity forces just overcome capillary trapping (the onset of buoyancy-driven oil mobilization) depends strongly on pore geometry: in ideal cylindrical pore throats, Bo_critical ≈ 4 (cos θ / (1+cos θ))² × (r_throat/r_pore)², where θ is the contact angle and r ratios are the throat-to-pore size ratio. In irregular, angular pore geometries (common in WCSB Devonian carbonates), the critical Bond number is higher because angular pores retain oil in corners by capillary action even when buoyancy would displace oil from the main pore body. This geometry dependence means that Bo_critical must be measured from laboratory drainage-imbibition experiments on representative core rather than calculated from simplified geometry assumptions — a reason why WCSB EOR pilot design always includes a core flooding program before field-scale implementation.
Surfactant EOR Design: Bond Number Analysis for a Cardium Waterflood Pool
An operator evaluating surfactant enhanced recovery on a mature Cardium pool at Pembina (current Sor 0.24, pore radius 4 µm average, oil-water IFT 28 mN/m, oil density 835 kg/m³, formation water 1,018 kg/m³) calculates the baseline Bond number: Bo = (183 × 9.81 × (4×10⁻⁶)²) / 0.028 = 1.03 × 10⁻⁷. The EOR target is to achieve Bo_critical of 5 × 10⁻⁴ to begin oil mobilization from trapped ganglia. Required IFT reduction: σ_target = (183 × 9.81 × (4×10⁻⁶)²) / (5×10⁻⁴) = 5.7 × 10⁻⁶ N/m = 0.0057 mN/m. Ultra-low IFT surfactant (sulfonated betaine blend) at 0.5 wt% concentration achieves σ = 0.004-0.008 mN/m in laboratory spindle tests with Cardium brine and Pembina crude — just reaching the target. Laboratory core flood at surfactant conditions recovers an additional 11% of pore volume (OOIP) beyond waterflood Sor, equivalent to approximately 8.5 million barrels incremental recovery at field scale. Economic threshold: surfactant slug volume required = 0.3 pore volumes at 0.5% concentration = 2.8 million kg of surfactant at CAD 15,000/tonne = CAD 42 million surfactant cost. At CAD 65/bbl incremental oil price and 8.5 million barrels, revenue = CAD 552 million. EOR project NPV positive at CAD 3.50/GJ gas cost (for waterflood injection pump energy) and 10% discount rate. Pilot design approved for 3-well pattern in the most favourable Cardium zone.
Fast Facts
The Bond number is named after Wilfrid Noel Bond, a British surface chemist who in the 1920s studied the behavior of bubbles and drops under the combined influence of surface tension and gravity — work that led to the Bond number as a formalization of the observation that large drops fall while small drops float, because in large drops (large L) gravitational forces overcome surface tension, while in small drops (small L) the reverse is true. Bond's work predates petroleum reservoir engineering by decades, but his dimensionless number became essential in petroleum science in the 1970s when laboratory studies by Larson, Davis, and Scriven at the University of Minnesota established the quantitative relationship between Bond number and residual oil saturation that underpins modern surfactant EOR design — work that forms the theoretical basis for every chemical enhanced recovery project proposed for mature WCSB waterfloods today.
Related Terms
The Bond number analysis for gravity drainage in SAGD thermal recovery connects directly to the steam generation system covered under boiler: the SAGD steam quality and injection rate determine the temperature (and therefore the oil viscosity reduction) at the steam chamber boundary, with higher steam quality delivering more latent heat and achieving lower bitumen viscosity that enables gravity drainage at lower Bond numbers for a given SAGD reservoir geometry. The capillary trapping that the Bond number quantifies at the pore scale is directly related to the pore size distribution and wettability measured by the NMR tools whose polarization physics are described under Boltzmann probability distribution — NMR T2 relaxation time distributions provide the pore size distribution that is the input for calculating the pore-scale Bond number distribution in heterogeneous WCSB reservoir rocks.