Body Waves in Seismic Exploration: P-Wave Velocity, S-Wave Anisotropy, and WCSB Formation Characterization
A body wave is a seismic wave that propagates through the interior of a solid or fluid medium — as opposed to surface waves that travel along an interface — carrying elastic energy outward from the seismic source through all formation layers in the subsurface. Two fundamentally different body wave types exist, each governed by different physics and encoding different information about the formations they traverse: P-waves (compressional or longitudinal waves) in which particle motion is parallel to the propagation direction (the rock alternately compresses and expands in the direction of wave travel), and S-waves (shear or transverse waves) in which particle motion is perpendicular to the propagation direction (the rock shears side-to-side while the wave moves forward). P-waves travel faster than S-waves in almost all geological materials — in typical WCSB formations, Vp (P-wave velocity) ranges from 2,000-2,800 m/s in shallow unconsolidated sands to 5,000-6,200 m/s in tight Montney siltstone and Devonian carbonates at depth, while Vs (S-wave velocity) for the same lithologies ranges from approximately 55-65% of Vp. The Vp/Vs ratio (or its equivalent, Poisson's ratio) is a uniquely powerful rock physics diagnostic: Vp/Vs is lower in gas-saturated sands than in water-saturated sands of the same porosity and lithology (because gas has negligible shear modulus, reducing bulk modulus and Vp without affecting Vs), a principle that underlies AVO (amplitude versus offset) analysis on WCSB 3D seismic datasets used to distinguish gas-charged Montney and Cardium zones from brine-saturated equivalents. S-waves cannot propagate through liquids (liquid has no shear modulus), meaning that while P-waves pass through both solid rock and interstitial fluid, S-waves respond only to the rock frame — a difference exploited in crosswell S-wave surveys and multicomponent geophones to separate lithology and fluid effects that P-wave data alone cannot disambiguate. Vertical seismic profiles (VSPs) conducted in WCSB Montney and Duvernay wells using 3-component borehole receivers record both P and S body waves generated by the surface seismic source, allowing direct measurement of interval velocity (and therefore depth conversion accuracy) that cannot be obtained from surface seismic data alone, and generating synthetic seismograms calibrated to the sonic log that are the primary tool for tying surface seismic reflections to specific formation depths in the well.
Key Takeaways
- Vp/Vs ratio as a gas indicator in WCSB seismic: The Gassmann equation quantifies how fluid substitution changes seismic velocity: replacing brine with gas in a reservoir reduces the bulk modulus (K_fluid drops from ~2.2 GPa for brine to ~0.02 GPa for gas), reducing Vp while leaving Vs nearly unchanged (because Vs depends on the shear modulus of the rock frame, not the fluid). A typical Montney tight siltstone at 3,000 m has Vp 5,200 m/s and Vs 3,100 m/s when brine-saturated (Vp/Vs = 1.68), dropping to Vp 4,900 m/s and Vs 3,090 m/s when gas-saturated (Vp/Vs = 1.59). The 0.09 reduction in Vp/Vs identifies gas presence and drives AVO Class IIp anomaly patterns used by WCSB explorationists to rank Montney and Duvernay prospects before committing exploration well capital.
- S-wave splitting and fracture characterization: In fractured WCSB reservoirs (Montney SRV, Devonian reef systems), S-waves traveling through the fractured rock split into two components: a fast S-wave polarized parallel to the fracture orientation and a slow S-wave polarized perpendicular to fractures. The time delay between fast and slow S-waves (the splitting delay, measured in milliseconds per kilometre of path length) is proportional to fracture density, and the fast S-wave polarization direction indicates the dominant fracture strike — which is approximately equivalent to the maximum horizontal stress orientation in most WCSB wells. Multicomponent surface seismic and VSP surveys use S-wave splitting to map natural fracture networks that guide horizontal well azimuth selection and hydraulic fracture design in Montney completions.
- P-wave normal moveout and velocity analysis in WCSB seismic processing: In a 3D seismic gather, P-waves reflected from a flat reflector display a characteristic hyperbolic moveout pattern (NMO — normal moveout) that allows the interval velocity above the reflector to be estimated from the curvature of the hyperbola. WCSB Montney 3D seismic datasets typically use bin sizes of 15-25 m and maximum offsets of 3,000-5,000 m, providing a wide-angle reflection dataset suitable for precise NMO velocity picking that yields interval velocities accurate to plus or minus 2-3% for depth conversion of the Montney horizon. Velocity errors of 2% propagate to depth errors of 60 m at 3,000 m TVD — significant enough to affect well placement in stacked Montney A/B/C landing targets spaced 50-100 m apart vertically.
- Body wave attenuation as a hydrocarbon indicator: P-waves attenuate (lose amplitude) more rapidly in gas-saturated reservoirs than in water-saturated formations because gas-brine capillary squirt flow at the grain-to-grain contacts dissipates seismic energy at the frequencies used in reflection seismic (10-100 Hz). This attenuation contrast produces the seismic attribute called "gas chimney" in 3D seismic volumes — a zone of anomalously low P-wave amplitude and high attenuation (low Q factor) above a gas accumulation or gas migration pathway. WCSB explorationists use gas chimney attributes from 3D seismic to identify leaky faults, gas migration paths from deeper Montney into shallower zones, and potential gas pockets hazardous to shallow drilling.
- VSP body wave recording in WCSB wells: A zero-offset VSP places a borehole receiver array at successive depth stations (typically every 15 m from surface to TD), recording downgoing and upgoing P and S body waves from a surface source. The downgoing P-wave arrival time at each receiver gives the interval transit time (reciprocal of velocity) more accurately than the sonic log alone (which averages over a 30-50 cm tool aperture and may be affected by borehole enlargement). For WCSB Duvernay and Montney wells, VSP interval velocities are used to calibrate the seismic velocity model to within 15-20 m depth accuracy versus the 40-60 m residual uncertainty in pure seismic velocity analysis — an improvement that can make the difference between landing a horizontal well in the sweet spot of a 30-m Montney pay interval versus missing it.
AVO Analysis: Montney Gas Sand vs Brine-Saturated Equivalent at Groundbirch
A WCSB explorationist uses 3D seismic AVO (amplitude versus offset) analysis to evaluate a Montney A zone gas prospect at Groundbirch. On the angle gathers (stacked into near, mid, and far angle substacks at 5-15°, 15-30°, and 30-45° respectively), the Montney A reflection shows a Class IIp AVO response: negative near-angle reflection coefficient (weak negative amplitude on the near stack) becoming increasingly negative with offset (strong negative amplitude on the far stack) — consistent with a gas-saturated sand above a brine-saturated equivalent. Gradient-intercept crossplot confirms the prospect falls in the gas-sand quadrant (negative intercept, negative gradient) with a separation of 0.08 Vp/Vs units from the brine trend established by calibration against 14 offset wells with known fluid contacts. The prospect is drilled as Montney A horizontal well: gas productive at 4.2 MMcf/day initial rate, confirming the AVO interpretation. The brine-wet alternative scenario (which would have shown a near-zero to positive AVO response) is confirmed absent on all 23 Montney A wells drilled in the Groundbirch area during the subsequent 3-year development program.
Fast Facts
The distinction between P-waves and S-waves was first established theoretically by Siméon-Denis Poisson in 1829, who derived from the equations of elasticity that two distinct wave modes must exist in a solid: one involving volumetric compression (what we now call the P-wave) and one involving shear deformation (the S-wave). Poisson also derived that the ratio Vp/Vs in a homogeneous isotropic solid must be at least sqrt(2) — approximately 1.414 — when Poisson's ratio (the ratio of lateral to longitudinal strain) is 0.25, the value for many consolidated sedimentary rocks. The practical application of Vp/Vs ratios for fluid discrimination was not exploited until AVO analysis became computationally practical in the 1980s, when digital 3D seismic acquisition and workstation-based processing enabled the amplitude-versus-offset analysis that exploits the fluid sensitivity of Vp/Vs on a routine exploration basis.
Related Terms
The borehole seismic applications that record body waves in WCSB wells — VSP surveys, crosswell tomography, and downhole microseismic — are covered under borehole seismic data, including the specific processing steps used to separate downgoing and upgoing waves on the borehole receiver array. The Vp/Vs rock physics framework that allows body wave velocity ratios to distinguish gas-saturated from brine-saturated formations — and that underlies most WCSB 3D seismic exploration for Montney and Duvernay gas — is the same framework applied in the formation evaluation crossplot methods described under bivariate analysis, where acoustic impedance versus Vp/Vs crossplots are the primary seismic reservoir characterization tool for ranking pre-drill prospect fluid predictions before committing exploration well capital.