Differential Weathering Correction

Key Takeaways

• The weathering layer (low-velocity zone, LVZ) is characterized by its dramatically lower seismic velocity relative to the underlying bedrock: unconsolidated dry sand and gravel near the surface has P-wave velocities of 200 to 500 m/s, saturated alluvium ranges from 600 to 1,500 m/s, and partially weathered rock transitions from 800 to 2,000 m/s to the unweathered bedrock velocity of 2,000 to 6,000 m/s; the thickness of the weathering layer is highly variable, ranging from less than 1 meter on exposed bedrock outcrops (desert pavement, glacially scoured shield) to more than 100 meters in deeply weathered tropical terrains, alluvial valleys, sand dunes (where the loose eolian sand has very low velocity throughout the dune thickness), or karst terrain (where collapse of carbonate cave systems creates highly irregular LVZ geometry); the key quantity in differential weathering correction is not the absolute thickness or velocity of the LVZ, but rather the lateral variation in the two-way traveltime through the LVZ from source to the weathering base (the replacement velocity datum) and from the weathering base to the receiver, because this lateral variation is what introduces the spurious time shifts in the reflection data that need to be removed.
• Refraction statics, derived from the analysis of first-arrival traveltimes in seismic data, are the most widely used method for differential weathering correction in land seismic processing: the first arrivals at short receiver offsets (less than the critical distance for refraction) are direct waves traveling near the surface through the low-velocity weathering layer, while at longer offsets the first arrivals are refracted waves that have traveled down to the top of the high-velocity bedrock, along the bedrock surface at the bedrock velocity, and back up to the receiver; by fitting a two-layer refraction model to these first-arrival traveltime curves, the weathering layer thickness and velocity can be estimated at each shot and receiver location across the entire seismic spread; the refraction static correction for each shot and receiver is then the difference between the actual traveltime through the LVZ and the traveltime that would have been observed if the LVZ had the uniform replacement velocity (typically 2,000 to 3,000 m/s) to the datum elevation; refraction statics work well when the base of the LVZ is a distinct refracting boundary with sufficient velocity contrast, but fail when the LVZ velocity gradient is continuous (no distinct refractor), when the refraction velocity is less than the LVZ velocity (hidden layers), or when the LVZ thickness exceeds the maximum refraction offset recorded in the data.
• Uphole surveys provide the most direct measurement of weathering layer properties by firing shots at depth in shallow boreholes (typically 20 to 100 meters deep) and recording the seismic arrivals at surface receivers; the traveltime from the shot at depth through the weathering layer to the surface receivers, combined with knowledge of the shot depth and the depth of the LVZ base (determined by the depth where the velocity transitions from LVZ values to bedrock values), directly gives the LVZ velocity and thickness at the uphole location; uphole surveys are acquired at selected intervals across a seismic survey area (typically every 1 to 5 km for regional 2D surveys, every 0.5 to 2 km for 3D surveys) and the results are interpolated between uphole locations to provide a continuous weathering model for the static correction calculation; uphole data provide a ground truth reference that refraction statics can be calibrated to, and are particularly valuable in complex near-surface environments (permafrost, thick sand dunes, rough karst terrain) where refraction first-arrival picking and inversion are difficult or ambiguous.
• Residual statics corrections are an iterative refinement step applied after the initial field statics (elevation and refraction corrections) to remove the remaining short-wavelength static errors caused by near-surface heterogeneities too complex to be fully resolved by the refraction model; residual statics are computed from the moveout-corrected reflection data by measuring the time differences between traces from the same common midpoint that should have identical reflection times if the NMO correction and field statics were perfect, and attributing these residual time differences to short-period shot and receiver statics that were not captured by the smooth refraction model; the residual static correction at each shot and receiver location is estimated by maximizing the cross-correlation of traces sharing the same shot or receiver (the surface-consistent residual static decomposition), and the corrections are typically limited to one or two wavelengths at the dominant frequency to prevent the algorithm from removing real structural relief and mistaking it for near-surface static error; residual statics are particularly important in areas with irregular near-surface conditions such as sand dunes, glacial till, river crossings, road cuts, or agricultural areas with depth-variable irrigation-induced saturation changes.
• Permafrost terrain presents some of the most severe differential weathering correction challenges in the petroleum industry because the active layer (the seasonally thawed zone above the permafrost table) has dramatically lower velocity than the frozen permafrost below, and its thickness varies seasonally (expanding downward in summer and contracting in winter) and spatially (deeper under lakes, shallower under elevated terrain, absent in taliks where heat flow maintains unfrozen conditions); Arctic seismic surveys must be acquired consistently during winter (when the active layer is frozen and part of the high-velocity permafrost column) or consistently during summer (when the active layer is fully thawed) to minimize temporal static variations; even so, permafrost thickness itself varies widely from tens of meters in sub-Arctic lowlands to hundreds of meters in continuous permafrost zones, creating long-wavelength differential weathering corrections of tens to hundreds of milliseconds that fundamentally limit the structural accuracy of reflection data; near-surface characterization for permafrost static correction uses a combination of uphole surveys, ground-penetrating radar (GPR), and refraction seismic to map permafrost thickness and active layer properties, with the resulting model providing the static corrections that make deep exploration reflection images interpretable.