Surface Wave: Love Waves, Rayleigh Waves, and Ground Roll Suppression in Land Seismic
A surface wave is a seismic wave that propagates along the interface between two media rather than through the interior of a single medium, with most of its elastic energy concentrated within a depth equal to about one wavelength of the boundary. In oil and gas seismic exploration, the two most operationally relevant surface wave types are Love waves and Rayleigh waves, both named after the British mathematicians who first described them: A.E.H. Love (Love waves, 1911) and Lord Rayleigh (Rayleigh waves, 1885). Rayleigh waves travel along the Earth-to-air interface with retrograde elliptical particle motion in the vertical plane containing the propagation direction, while Love waves are horizontally polarized shear waves trapped between a low-velocity surface layer and a higher-velocity substrate. Rayleigh waves are dispersive in layered media, meaning different frequencies travel at different velocities, and this dispersion is the basis for Multichannel Analysis of Surface Waves (MASW) used for near-surface shear-wave velocity profiling in geotechnical and seismic exploration applications. In WCSB land seismic acquisition, surface waves manifest as "ground roll", a slow, high-amplitude, low-frequency wavetrain that dominates near-source traces at velocities of 200 to 1,200 m/s (656 to 3,937 ft/s) compared to body wave velocities of 3,000 to 6,000 m/s. Ground roll overlies and obscures useful reflection signals from deeper targets including the Montney, Duvernay, Cardium, and Mannville formations, and its suppression is a primary objective of seismic processing flows including frequency-wavenumber (FK) filtering, radial-trace filtering, surface wave deconvolution, and adaptive subtraction. Acquisition crews also use field techniques such as vibroseis sweeps biased to higher frequencies (8 to 96 Hz instead of 2 to 96 Hz) and geophone arrays designed to attenuate ground roll through destructive interference. Surface wave behaviour is also fundamental to earthquake seismology and the assessment of induced seismicity associated with hydraulic fracturing in the Duvernay and Montney plays, where regulators including the AER (Subsurface Order 2) and BCER track magnitudes and epicentre locations using surface and body wave arrivals at regional broadband stations. The discipline connects to seismic wave taxonomy, ground roll suppression, and induced seismicity monitoring across the WCSB unconventional fairways.
Key Takeaways
- Rayleigh Versus Love Waves: Rayleigh waves involve coupled P and SV motion with retrograde elliptical particle paths in the vertical plane and require only a free surface to exist. Love waves are pure horizontally polarized shear (SH) motion and require a low-velocity surface layer over a higher-velocity substrate to be trapped. Both are dispersive, meaning velocity depends on frequency, with low-frequency components sampling deeper than high-frequency components, a property exploited in MASW surveys for near-surface characterization to depths of 30 to 100 m below ground level.
- Ground Roll Characteristics in WCSB Land Seismic: Ground roll velocities in WCSB plains and foothills typically range from 300 to 900 m/s with dominant frequencies of 6 to 18 Hz, while target reflections from the Montney or Duvernay arrive at 4,000 to 4,500 m/s with frequencies of 20 to 60 Hz. The order-of-magnitude velocity and frequency separation enables FK and FX filtering to attenuate ground roll by 20 to 40 decibels without affecting the reflection signal, recovering otherwise obscured target events from the deep section.
- MASW for Near-Surface Velocity Profiling: Multichannel Analysis of Surface Waves uses the dispersion of Rayleigh waves to invert for shear-wave velocity (Vs) versus depth in the upper 30 to 100 m. This is critical for static corrections in WCSB land seismic, where weathering layers vary from 5 to 40 m thick and have Vs of 150 to 800 m/s. MASW surveys cost CAD 8,000 to CAD 25,000 for a 200 m line and dramatically improve static solutions in muskeg or glacial till near-surface conditions common across northern Alberta and northeastern BC.
- Vibroseis and Geophone Array Design: Vibroseis trucks generate seismic energy with controlled sweep frequencies. Operators commonly use sweeps starting at 6 to 10 Hz rather than 2 to 4 Hz to minimize ground roll generation, accepting reduced low-frequency reflection content as a tradeoff. Geophone arrays of six to twelve phones spaced over 6 to 24 m create directional sensitivity that attenuates surface waves through destructive interference, leveraging the slow apparent velocity of ground roll across the array compared to near-vertical reflections from depth.
- Induced Seismicity Monitoring: Surface wave arrivals at regional broadband stations are routinely used to locate and characterize induced earthquakes triggered by Duvernay or Montney hydraulic fracturing. The AER and BCER traffic-light protocols use local magnitudes derived from local seismic networks, with surface wave magnitudes (Ms) reported for events exceeding local magnitude 4.0. AER Subsurface Order 2 requires real-time monitoring and operational shutdown when local magnitude exceeds 4.0 within 5 km of active operations in red-light zones.
FK Filtering Workflow in Montney 3D Seismic
A typical Montney 3D processing flow attacks ground roll in the early stages. After initial trace edits and amplitude recovery, processors transform shot gathers from the time-offset (t-x) domain to the frequency-wavenumber (f-k) domain via 2D FFT. Ground roll energy appears as a coherent linear streak at low frequency and slow apparent velocity, distinct from the higher-velocity reflection hyperbolas. A surgical f-k filter zeroes the ground roll quadrant before inverse FFT returns the data to t-x. Typical FK filter parameters reject apparent velocities of 200 to 1,000 m/s and frequencies of 4 to 20 Hz, removing roughly 30 dB of ground roll while preserving target reflection bandwidth from the Doig phosphate marker through the Montney C-zone.
Surface Wave Magnitude in Duvernay Induced Seismicity
The January 2016 local magnitude 4.8 induced earthquake near Fox Creek, Alberta, triggered by hydraulic fracturing in the Duvernay, was characterized using both local body waves and regional surface waves at distances of 50 to 500 km. Surface wave magnitudes (Ms) derived from 20-second-period Rayleigh waves at stations including the Canadian National Seismograph Network confirmed the event size and supported the AER's traffic-light response. The operator's drilling and fracking program was paused, and the event prompted Subsurface Order 2 (later refined in 2018) imposing real-time monitoring and shutdown thresholds across the entire Duvernay play in Alberta.
Fast Facts
Lord Rayleigh predicted the existence of the waves now bearing his name in 1885 in a single page of his "Theory of Sound" volume, decades before geophones could record them. A.E.H. Love followed in 1911 with his prediction of horizontally polarized surface waves, which he derived to explain why earthquake records showed two distinct horizontal motions arriving at different times. Both predictions preceded modern seismology by decades, and yet every WCSB 3D seismic survey today still designs its acquisition geometry and processing flow to mitigate the very waves these two mathematicians predicted from pure elastic theory.
Related Terms
Surface waves sit within a broader seismic wave taxonomy described in related glossary entries. Seismic wave is the overarching category that includes body waves (P and S), surface waves (Rayleigh and Love), and guided waves. P-wave and S-wave are the compressional and shear body wave counterparts that carry useful reflection energy through the subsurface. Ground roll is the operational term for the Rayleigh wave noise that dominates near-source traces in land seismic and is the primary target of FK and radial filtering during processing.
WCSB Field Scenario: Ground Roll Suppression on a Duvernay 3D Survey Near Kaybob
An operator acquired a 180 square kilometre 3D seismic survey over a Duvernay block near Kaybob, Alberta, in winter 2024 using vibroseis trucks on 50 by 50 m source spacing with 100 m by 200 m receiver geometry. Initial shot gathers showed dominant ground roll at 8 to 14 Hz with apparent velocities of 350 to 650 m/s, obscuring Duvernay reflections at 1.6 to 1.8 seconds two-way time. Standard f-k filtering reduced ground roll by 28 dB but introduced minor frequency artifacts on the Duvernay reflectors. The processor switched to a hybrid radial-trace and adaptive subtraction flow that delivered 36 dB ground roll suppression with cleaner preservation of the Duvernay signal envelope.
Total processing cost for the upgraded flow ran CAD 240,000 over the 180 square kilometre survey, versus CAD 165,000 for the standard flow, but the resulting Duvernay reflector quality supported tighter horizontal well placement and reduced lateral redrills. Final inversion of the dataset yielded acoustic impedance and density attributes with 12 percent improved resolution over the standard processing flow, justifying the incremental processing spend across the eight-well first-phase development program.