Normal-Moveout Correction: CMP Stacking, Velocity Analysis, and 3D Seismic Processing for WCSB Exploration

Normal-moveout correction (NMO) is a seismic data processing step that compensates for the geometric time delay introduced when a downgoing acoustic wave reflects off a subsurface interface and returns to receivers located at progressively greater horizontal distances from the source. For a flat reflector beneath horizontally layered ground, the two-way travel time at offset x is approximated by the hyperbolic equation t(x) = sqrt(t0^2 + (x/v)^2), where t0 is the zero-offset two-way time and v is the root-mean-square (RMS) velocity to the reflector. NMO correction subtracts the offset-dependent excess time from each trace so that, after correction, all reflections from the same subsurface point arrive at the same time across the common midpoint (CMP) gather. The corrected traces are then stacked, summing constructively at the reflector and averaging out random noise, which boosts signal-to-noise ratio dramatically. The entire workflow assumes hyperbolic moveout, which is exact only for a single flat layer with constant velocity; in real WCSB stratigraphy with dipping layers, lateral velocity variations, and anisotropic shales, the assumption is approximate, and modern processors use higher-order moveout corrections (fourth-order or non-hyperbolic NMO) or prestack depth migration to handle the residuals. NMO correction is the workhorse of conventional 2D and 3D seismic processing and was the foundational step that made the common-depth-point (CDP) method economically viable when introduced by Harry Mayne at Petty-Ray Geophysical in the late 1950s. Velocity analysis, the iterative selection of the RMS velocity function that best flattens the gathers, is performed on velocity spectra (semblance panels) typically computed at 250 to 500 m intervals across a survey. Picked velocities are then interpolated and used to apply the NMO correction across the full volume. In the Western Canadian Sedimentary Basin, NMO correction is applied to every commercial seismic dataset, from shallow Mannville exploration grids at Cenovus to deep Slave Point and Leduc carbonate plays at Rainbow Lake. Particular care is needed for the Duvernay and Montney intervals where rapid lateral velocity changes (related to porosity and TOC variability) and significant anisotropy (Thomsen epsilon values of 0.10 to 0.25) introduce residual moveout that requires anisotropic NMO or full-azimuth processing. NMO stretch, the unavoidable distortion of waveforms at far offsets, is muted out before stack to prevent the corrupted samples from polluting the stacked trace. The correction integrates closely with seismic attribute analysis and AVO analysis, both of which rely on accurate NMO-corrected prestack data.

Velocity Analysis Workflow

The processor builds a velocity function by examining semblance panels at supergather locations spaced 250 to 500 m apart across the survey. At each location, the semblance spectrum displays coherence as a function of two-way time and RMS velocity, with high-amplitude peaks indicating the velocity that best flattens that reflector. The interpreter picks velocities along major events: a Mannville Channel sand at 1,200 ms two-way time might have RMS velocity of 2,400 m/s, while a deeper Beaverhill Lake reef at 2,400 ms might pick at 4,100 m/s. Picks are quality-controlled against the geological model, interpolated laterally, and applied as the NMO velocity field. The entire process can be iterated to refine the residual moveout on key target horizons.

NMO Stretch, Muting, and Higher-Order Corrections

NMO stretch is the unavoidable consequence of applying a time shift that varies along the waveform. A reflection 40 ms wide at zero offset may stretch to 50 to 60 ms at a 3,000 m offset, lowering the dominant frequency from 40 Hz to perhaps 27 Hz. Far-offset traces beyond about 30% stretch are muted to preserve resolution. For deep targets and long offsets typical in WCSB Duvernay surveys, fourth-order non-hyperbolic NMO using the eta parameter is applied to recover far-offset signal that would otherwise be muted, improving AVO analysis on prestack gathers and supporting better lithology and fluid discrimination.

Related Terms

NMO correction sits in the middle of a long seismic processing chain. Common midpoint is the gather geometry that NMO correction operates on, while seismic velocity (specifically RMS velocity) is the input the correction requires. Stacking is the constructive summing step that follows NMO and produces the signal-to-noise improvement, and migration is the downstream process that repositions reflectors to their true subsurface locations. Together these terms describe the sequence that turns raw field records into the interpretable seismic image used by geophysicists across the WCSB.

WCSB Field Scenario: Montney 3D Reprocessing at Karr

A WCSB operator commissions a reprocessing job over a 280 km² Karr-area Montney 3D survey originally processed in 2014 with isotropic NMO. The new flow at a Calgary seismic processing house, costing approximately CAD 425,000, applies anisotropic NMO with picked eta values of 0.15 to 0.22 across the Montney interval, full-azimuth gather conditioning, and reverse-time migration. The reprocessed AVO product reveals previously hidden brittleness and porosity anomalies in the upper Montney C interval at depths of 2,950 to 3,100 m subsurface.

The operator uses the new attribute volume to high-grade four horizontal well locations on its 2026 development program, each costing roughly CAD 9.5 million drill-complete-tie-in. The reprocessing pays out in less than 90 days through improved well placement that adds approximately 12% to expected ultimate recovery per well, a return that is now the standard business case for anisotropic NMO reprocessing across overpressured unconventional plays in northeast BC and northwest Alberta.