Amplitude Variation with Offset
Amplitude variation with offset (AVO) is the change in seismic reflection amplitude as a function of source-receiver distance (offset) or, equivalently, as a function of the angle of incidence of the seismic ray at the reflecting interface. The physical basis of AVO is the Zoeppritz equations (1919), which describe how the partitioning of energy between reflected and transmitted compressional (P-wave) and shear (S-wave) waves at an interface depends on the incidence angle and on the contrasts in compressional velocity (Vp), shear velocity (Vs), and density (ρ) across the interface. At normal incidence (zero offset), the reflection amplitude depends only on the compressional acoustic impedance contrast (Vp × ρ); as the angle of incidence increases, the amplitude increasingly reflects both the Vp contrast and the Vs contrast across the interface, because shear wave conversion becomes significant. Since gas saturation in a porous rock strongly reduces Vp (lowering pore fluid bulk modulus from approximately 2.2 GPa for brine to approximately 0.02 GPa for gas) but has minimal effect on Vs (which is insensitive to pore fluid under Gassmann's theory because the shear modulus of a rock is controlled by the grain framework, not the fluid), gas-saturated sands have anomalously low Vp/Vs ratios (approximately 1.5) compared to brine-saturated sands (approximately 1.7 to 2.0) and shales (approximately 1.8 to 2.2). This Vp/Vs contrast is amplified in the AVO response (the angle-dependent amplitude change), making AVO the primary seismic technique for distinguishing gas saturation from lithological changes in WCSB Viking, Cardium, and Deep Basin tight gas exploration. The Shuey (1985) two-term approximation expresses the angle-dependent reflection coefficient as R(θ) ≈ R₀ + G × sin²θ, where R₀ is the normal-incidence reflectivity (AVO intercept, proportional to acoustic impedance contrast) and G is the AVO gradient (proportional to the Vs/Vp contrast change across the interface); AVO classes I through IV are defined by the sign and magnitude of R₀ and G, with Class III (negative R₀, negative G) being the most common gas sand signature in WCSB shallow clastic plays, and Class IIp (near-zero R₀, negative G causing polarity reversal with offset) being diagnostically valuable for reservoirs at intermediate impedance contrast where normal-incidence amplitude alone provides weak discrimination between gas and brine.
Key Takeaways
- AVO is classified into four standard classes based on the sign and magnitude of the normal-incidence reflectivity (intercept R₀) and the AVO gradient (G), with each class corresponding to a distinct geological and fluid scenario: Class I (positive R₀, negative G) indicates a reservoir with higher impedance than the surrounding shale (hard reflection that dims with offset), Class II (near-zero R₀) is a polarity-sensitive crossover case, Class III (negative R₀, negative G) is the diagnostic bright-spot gas signature in low-impedance sand below shale, and Class IV (negative R₀, positive G) occurs in unusual geological settings where gas sand amplitude decreases with offset despite being a soft reflector: Class I AVO occurs where the reservoir sand (even when gas-saturated) has higher acoustic impedance than the encasing shale, which is typical of deeply buried tight sandstones with significant quartz cementation (Deep Basin Cadomin and Spirit River formations in northwest Alberta, depths 2,500 to 4,500 m). In this case, gas saturation causes the impedance to decrease from the brine-saturated value but not below the shale impedance, so the reflection remains positive-polarity (hard) at normal incidence but decreases with offset as the Vs constraint reduces the gradient. Class III AVO (negative R₀, negative G) is the most prevalent gas sand signature in WCSB shallow exploration (Viking, Cardium, Mannville at 500 to 1,500 m depth), where low confining stress, high porosity, and high gas saturation combine to produce gas-sand impedance well below the shale impedance, generating bright spots that increase in amplitude with offset. The AVO class is not a fixed property of a formation but depends on the specific Vp-Vs-density relationships of both the reservoir and the surrounding shale at the burial depth and compaction state of the target.
- The AVO intercept-gradient crossplot is the primary diagnostic tool for identifying gas sands in pre-stack seismic data, with the background trend (shale-shale reflections plotting near the origin with a negative slope defined by the local mudrock line) providing a reference against which hydrocarbon anomalies plot as outlier clusters in the quadrant corresponding to the expected AVO class for the target formation, enabling semi-automatic identification of anomalous reflectors from hundreds of millions of seismic data points: On an intercept (R₀) versus gradient (G) crossplot, elastic reflections from shale-shale interfaces follow a linear trend (the background or wet trend) with a negative slope reflecting the normal Vp/Vs ratio in the shale sequence. Gas sand reflectors at the top of a Class III gas reservoir plot in the third quadrant (negative R₀, negative G), displaced from the background trend by the anomalously low Vs/Vp ratio of the gas sand relative to the surrounding shale. Brine-saturated sands with the same lithology plot near the background trend on the crossplot, providing the discrimination between gas and brine that normal-incidence amplitude alone cannot provide when the acoustic impedance contrast is similar for both fluid states. In WCSB Viking play evaluation, intercept-gradient crossplot analysis of 3D pre-stack seismic data routinely identifies 15 to 30 anomalous reflectors per 100 km² of 3D seismic coverage, of which 40 to 60% are subsequently ranked as medium to high confidence DHI prospects by AVO class, amplitude conformance to structure, and well calibration consistency.
- Gassmann's equations (1951) provide the rock physics foundation for AVO by predicting the bulk modulus (and therefore Vp) of a porous rock saturated with different fluids, enabling quantitative fluid substitution modelling that calculates the expected AVO response for gas saturation, oil saturation, and brine saturation from well log data, validating that the observed seismic AVO anomaly is consistent with the expected fluid contrast before committing to a drilling decision: Gassmann's relation: K_sat = K_dry + (1 - K_dry/K_mineral)² / (φ/K_fluid + (1-φ)/K_mineral - K_dry/K_mineral²), where K_sat is the saturated-rock bulk modulus, K_dry the dry-frame modulus, K_mineral the mineral grain modulus, K_fluid the pore fluid modulus, and φ the porosity. By substituting K_fluid values for brine (2.2 GPa), oil (0.7 to 1.5 GPa depending on API gravity), and gas (0.02 to 0.10 GPa depending on pressure and composition), the geophysicist can calculate the Vp for each fluid scenario from the same K_dry measured on dry core plugs (or estimated from log Vp and the Biot coefficient), then compute the expected AVO response for each scenario and compare to the observed seismic AVO. For WCSB Viking gas sands, Gassmann modelling using core-plug K_dry values typically predicts Class III AVO with gradient G = -0.04 to -0.09 for gas saturation above 50%, while brine substitution predicts G = -0.01 to +0.02, providing a clear discrimination between gas and brine fluid states on the AVO gradient measurement.
- Azimuthal AVO (AVAZ) analysis uses the variation of AVO gradient as a function of seismic acquisition azimuth to detect seismic anisotropy caused by aligned natural fractures or in-situ stress anisotropy, providing fracture orientation and intensity information from 3D seismic data that guides horizontal well landing orientation and perforation cluster spacing decisions in naturally fractured Montney and Duvernay reservoirs in Alberta: In a naturally fractured reservoir, aligned vertical fractures create seismic velocity anisotropy: P-wave velocity propagating perpendicular to the fractures is lower than velocity propagating parallel to the fractures (because fractures open in the perpendicular direction and reduce the elastic stiffness). This velocity anisotropy translates into azimuth-dependent AVO responses: the AVO gradient is larger (more negative for a soft reflector) in the azimuth perpendicular to the fractures than in the azimuth parallel to the fractures. Fitting the AVO gradient variation with azimuth (using a cos²θ azimuth function) yields the fracture orientation (from the direction of maximum gradient) and the fracture intensity (from the magnitude of the azimuthal gradient variation). In the Montney Formation of northeast British Columbia (Dawson Creek, Pouce Coupe, Encana and ConocoPhillips acreage), AVAZ analysis of wide-azimuth 3D seismic has identified NW-SE to NNW-SSE fracture orientations (consistent with the regional maximum horizontal stress direction) that are used to orient horizontal well laterals perpendicular to the fractures and to select perforation cluster spacings that maximise interaction with the natural fracture network, improving stimulated reservoir volume (SRV) and first-year cumulative production by an estimated 10 to 25% over wells landed and completed without AVAZ guidance.
- AVO analysis requires careful data quality control and processing to preserve the angle-dependent amplitude variations that carry the fluid information, including offset-balancing of near and far trace amplitudes, removal of amplitude distortions from acquisition geometry (source and receiver coupling variations), residual normal moveout correction (to ensure that reflections are properly aligned before AVO extraction), and overburden Q-compensation (to correct for the frequency-dependent attenuation that distorts amplitude differently at near and far offsets), with any uncorrected error in these steps potentially creating false AVO anomalies that are mistaken for hydrocarbon indicators: The most common AVO artefact in WCSB 3D seismic data is near-surface coupling variations: the amplitude of seismic energy recorded at receivers in soft near-surface materials (glacial lake clays, peat bogs, river valley alluvium) is attenuated differently than at receivers in hard near-surface materials (compact till, limestone outcrop), creating receiver-term amplitude variations that appear in the recorded data as offset-dependent amplitude changes (because different offsets use different receiver locations). Uncorrected, these near-surface coupling variations produce false AVO gradients that can be mistaken for lithological or fluid signals. Amplitude-preserving surface-consistent amplitude corrections (SCAC) remove the source and receiver coupling terms from the data in a least-squares fit using all reflections from multiple horizons, enabling AVO gradient extraction with accuracy sufficient to discriminate Class III gas sands from Class IIp background in WCSB shallow gas exploration. High-quality AVO results in the WCSB require raw data fold greater than 60 (typically 80 to 120 for DHI-grade data quality) and signal-to-noise ratio above 8 dB at the target level, setting the minimum 3D seismic survey design specification for AVO-capable datasets.
AVO Classes and WCSB Formation Examples
The four standard AVO classes were first formally defined by Rutherford and Williams (1989) in Geophysics and have since been applied globally to classify gas sand reflections. In the WCSB, the depth and compaction state of the target formation largely determine which AVO class is observed. Viking Formation gas sands (500 to 900 m depth, porosity 22 to 30%, unconsolidated to weakly consolidated) are consistently Class III: at shallow depth, the low confining stress means the dry-frame bulk modulus is small, so gas saturation strongly reduces the total bulk modulus and Vp, driving gas-sand impedance well below the shale impedance. The Cadomin and Spirit River tight gas sands of the Deep Basin (2,500 to 4,500 m depth, porosity 8 to 14%, heavily cemented) are Class I: at depth, cementation increases the dry-frame bulk modulus, so even gas saturation cannot reduce the gas-sand impedance below the shale impedance, and the AVO gradient is negative but the intercept remains positive. The Montney siltstone of northeast BC (1,500 to 3,500 m depth, porosity 4 to 12%, mixed dolomite-quartz-clay mineralogy) shows variable AVO class from I to II depending on the specific Montney member and depth, making AVO analysis a powerful discriminator for sweet spot identification when calibrated to the local Montney rock physics.