Velocity Correction: NMO, Statics, Migration Velocity, and Seismic Imaging in the WCSB
A velocity correction in seismic processing is any adjustment applied to recorded data so that reflectors appear at their correct two-way travel time, true depth, and accurate lateral position, compensating for the fact that seismic waves travel at velocities ranging from roughly 1,500 m/s in shallow water-saturated sediments to over 6,500 m/s in deep Devonian carbonates and Precambrian basement of the Western Canadian Sedimentary Basin. Several distinct velocity corrections are routinely applied during the processing flow: normal moveout (NMO) correction, which removes the offset-dependent travel time delay between source and receiver pairs at non-zero offsets so that reflections from the same subsurface point align horizontally across a common-midpoint (CMP) gather before stacking; statics corrections, which adjust for irregular surface elevation, weathered layer thickness, and near-surface velocity variations that would otherwise scramble traveltimes shot-by-shot in the rolling Foothills, Peace River Arch, and Plains terrain; dip moveout (DMO) and pre-stack migration velocity corrections, which honour reflector dip and lateral velocity variations during the imaging step; and depth conversion velocity functions, which transform the final stacked or migrated time section into a true depth section that can be tied to well control. All velocity corrections require assumptions about the seismic velocity distribution in the rocks and sediments through which the waves have travelled, and these assumptions are themselves the largest single source of imaging uncertainty in modern WCSB seismic interpretation. Velocities used in correction can be derived from sonic well logs (such as a DT log calibrated through a vertical seismic profile or check shot survey), from semblance velocity analysis on CMP gathers (the classic Dix-equation interval velocity approach), from full-waveform inversion run on the recorded data, or from tomographic refraction inversion of first-break picks. Each of these has accuracy limits typically ±2 to ±8 percent of true interval velocity, which translates into structural depth uncertainty of ±25 to ±150 m at typical Montney and Duvernay target depths of 2,500 to 3,800 m below ground level in west-central Alberta. The most common operator mistake is to over-trust a single velocity model. Best practice on a major WCSB exploration prospect today is to build velocity models that honour all available constraints, including sonic logs from offset wells, VSP check shots from the most recent vertical pilot hole, stacking velocities from the seismic itself, and refraction tomography for the near surface, then perform sensitivity testing on the depth conversion to see how the structural mapping changes under reasonable velocity perturbations. Velocity corrections also have a major impact on amplitude versus offset (AVO) analysis used for fluid prediction: an incorrect NMO velocity will distort the AVO gradient and lead to false-positive or false-negative direct hydrocarbon indicators, particularly in WCSB clastic plays such as the Glauconite, Sparky, Belly River, and Viking where AVO is the primary fluid-discrimination tool. Cost implications are material. A poor velocity model contributing to a missed structural closure at the AER Mississippian Pekisko or Devonian Leduc target level can mean a CAD 8 to 15 million dry hole in the Foothills or a missed reserves booking under NI 51-101.
Key Takeaways
- NMO correction fundamentals: Normal moveout correction removes the offset-dependent travel-time delay across a CMP gather, allowing reflections from the same subsurface point to align horizontally before stacking. The hyperbolic NMO equation requires an accurate stacking velocity Vrms at every reflector; incorrect velocity values produce residual moveout that smears the stack and destroys the AVO gradient relied on for fluid prediction.
- Statics correction in WCSB terrain: Statics adjust for irregular surface elevation and weathered-layer velocity variations and are critical in rolling Foothills, glaciated Peace River Arch, and Athabasca topography. Refraction statics from first-break picks plus surface-consistent residual statics together typically improve stack signal-to-noise by 6 to 12 dB and resolve subtle Devonian and Mississippian structural features otherwise masked by near-surface noise.
- Velocity sources and accuracy: Velocity models are built from sonic logs (DT), VSP check shots, stacking velocity semblance, refraction tomography, and increasingly full-waveform inversion. Each has accuracy of ±2 to ±8 percent of true interval velocity, translating into structural depth uncertainty of ±25 to ±150 m at typical 2,500 to 3,800 m WCSB target depths.
- Depth conversion ambiguity: Time-to-depth conversion uses an average velocity model and is the final step before reserves volumetrics. Sensitivity testing on alternative velocity functions is required best practice for any NI 51-101 reserves booking, because a 5 percent velocity error at 3,500 m depth shifts the mapped structural closure depth by 175 m and can move a prospect across the spill point.
- AVO and fluid prediction impact: Incorrect NMO velocity distorts the AVO gradient and intercept used for direct hydrocarbon indication in WCSB clastic plays (Glauconite, Sparky, Belly River, Viking). Industry studies report 15 to 30 percent of AVO false positives in published case histories are traceable to velocity analysis errors rather than rock physics ambiguity.
Velocity Picking Workflow on Modern WCSB 3D Data
A typical Calgary processing shop velocity workflow begins with brute-stack velocities semblance-picked on 100 m by 100 m supergather grids, followed by surface-consistent refraction statics, then dense interval velocity picking on a 250 m grid using semblance plus dynamic correction QC overlays. A residual moveout pass on the migrated gathers fine-tunes the model. Tomographic update inversion is then applied for areas of strong lateral velocity variation such as the Bashaw reef complex or the Crowsnest dolomite belt. Final velocity volumes are exported as both interval (Vint) and RMS (Vrms) cubes for migration and depth conversion respectively, with the entire workflow taking 4 to 8 weeks for a 200 km2 3D survey at a cost of CAD 320,000 to CAD 540,000.
Depth Conversion and Reserves Volumetrics
Once a final time-migrated 3D volume is produced, structural mapping at the Montney, Duvernay, Cardium, or Viking target reflector level is converted to depth using a velocity model anchored to offset well sonic logs and corrected for systematic time-depth mistie. Operators report multiple depth-conversion scenarios in their NI 51-101 reserves filings, with the chosen base case explicitly justified against alternate velocity functions. Reserves auditors at GLJ Petroleum Consultants and McDaniel Associates routinely require depth-uncertainty sensitivity analyses showing how the contingent and probable categories respond to velocity perturbation, and have rejected booking submissions where this analysis is absent or unrealistic.
Fast Facts
The Dix equation, used since 1955 to convert stacking velocities to interval velocities for depth conversion, was published by Conrad Dix while working at California Research Corporation (later part of Chevron) and is still cited in every modern WCSB processing report. Dix demonstrated mathematically that if you know the RMS velocity at the top and bottom of a layer, the interval velocity through that layer is the square root of the difference of squared RMS velocities divided by the difference of their travel times, a deceptively simple formula that underpins virtually every velocity correction step in seismic data processing today.
Related Terms
Velocity corrections are fundamental to normal moveout processing and to seismic migration, both of which require accurate interval velocity fields. The velocity model itself is typically anchored by a vertical seismic profile or check shot survey, and the final time-to-depth transformation feeds directly into structural mapping and reserves volumetric calculations under NI 51-101 standards. Errors in any of these steps propagate through the entire interpretation chain.
Pieridae Foothills Velocity Model Update, 2026
A 2026 Pieridae Energy reprocessing of legacy 1998 and 2007 Alberta Foothills 3D seismic volumes covering 412 km2 near Coleman included a comprehensive velocity model update driven by 14 newly acquired offset VSPs and 7 long-offset refraction lines. The reprocessing cost CAD 1.8 million through a Calgary processing contractor over a 14-week schedule and reduced average structural depth uncertainty on Mississippian Banff and Pekisko reflectors from ±180 m to ±55 m as measured against subsequent well control. Three of the previously high-graded prospects shifted significantly in mapped position once the new velocity model was applied to time-to-depth conversion.
One prospect previously mapped as a 12 km2 four-way structural closure resolved to a 3 km2 closure once the corrected velocity field repositioned the Banff peak by 145 m laterally and 38 m deeper. The reduced volumetric reserve estimate (P50 from 95 to 22 Bcfe) caused Pieridae to drop the prospect from the 2026 drilling programme, avoiding an estimated CAD 11 million dry-hole cost on a target that had previously been ranked second in the Foothills inventory.