Reflection Coefficient
The reflection coefficient in seismic exploration is a dimensionless quantity (ranging from -1 to +1) that describes what fraction of an incident seismic wave's amplitude is reflected at an interface between two formations with different acoustic impedances — calculated as R = (Z₂ - Z₁) / (Z₂ + Z₁) where Z₁ is the acoustic impedance (density × P-wave velocity) of the upper layer and Z₂ is the acoustic impedance of the lower layer — with the reflection coefficient determining the strength and polarity of seismic reflections observed on seismic records and serving as the fundamental link between seismic amplitude data and subsurface lithology, porosity, and fluid content that geophysicists use to interpret seismic data for hydrocarbon exploration and reservoir characterization.
Key Takeaways
- Acoustic impedance (Z = ρ × Vp, the product of formation density and P-wave velocity) is the formation property that controls seismic reflection — a large acoustic impedance contrast between two formations produces a strong reflection (high-amplitude event on the seismic section); a small contrast produces a weak reflection; zero contrast produces no reflection (the seismic wave passes through without reflection); the magnitude and polarity of the reflection coefficient determine whether the seismic reflection appears as a peak (positive R, high-impedance layer below low-impedance) or trough (negative R, low-impedance layer below high-impedance) on a zero-phase seismic section processed with standard polarity conventions.
- Fluid substitution effects on reflection coefficient are the basis of amplitude-versus-offset (AVO) analysis for direct hydrocarbon indication — replacing brine with gas in a porous sandstone reduces both density (gas is less dense than brine) and P-wave velocity (gas compresses more easily than brine), reducing the acoustic impedance of the sand and changing the reflection coefficient at the top of the gas sand; this acoustic impedance reduction produces a negative reflection coefficient at the top of the gas sand if the overlying shale has higher impedance, creating the "bright spot" high-amplitude reflection anomaly associated with gas sands in young clastic sections; interpreting this AVO response correctly requires understanding how the reflection coefficient changes with angle of incidence (since P-wave to S-wave conversion varies with angle, captured by the Zoeppritz equations or their Shuey approximations).
- The seismic wavelet convolution model describes a seismic trace as the mathematical convolution of the earth's reflection coefficient series (the sequence of R values at each formation boundary from shallow to deep) with the seismic source wavelet plus noise — this model is the foundation of seismic deconvolution, which attempts to remove the wavelet effect from the seismic trace to recover a closer approximation of the reflection coefficient series (the earth's impedance contrasts); the reflection coefficient series derived from well logs (by computing Z at each depth point and calculating R at each boundary) provides the calibration tie between the seismic trace at the well location and the actual formation boundaries intersected by the well.
- Normal moveout (NMO) correction in seismic processing removes the time delay between near-offset and far-offset seismic traces caused by the varying travel path length for different source-receiver offsets, making reflections from the same subsurface boundary line up correctly across the receiver array before stacking; after NMO and stacking, the stacked reflection amplitude is approximately proportional to the zero-offset reflection coefficient modified by the average AVO effect across the offset range, and amplitude anomalies in the stacked seismic data are interpreted as reflection coefficient variations caused by changes in lithology or fluid content along the reflector.
- Synthetic seismogram generation uses the reflection coefficient series computed from sonic log (velocity) and density log (density) to produce a simulated seismic trace at the well location that can be compared directly with the actual seismic data — the well-to-seismic tie quality (the visual similarity between the synthetic seismogram and the nearby seismic trace) validates the time-depth conversion function, confirms the seismic polarity convention, and identifies which seismic events correspond to specific formation boundaries seen in the wireline logs; a good synthetic-to-seismic tie with cross-correlation coefficient above 0.7 is required for reliable seismic interpretation calibrated to formation tops in the well.
Fast Facts
The Zoeppritz equations (published by Karl Zoeppritz in 1919) are the exact mathematical description of seismic reflection and transmission coefficients at an interface for all angles of incidence, accounting for the conversion of P-waves to S-waves and vice versa; for most practical seismic exploration purposes at angles less than 30 to 35 degrees, the simpler Shuey two-term and three-term approximations to the Zoeppritz equations provide sufficient accuracy for AVO analysis and are computationally far simpler; the Shuey approximation expresses the reflection coefficient as a function of angle as R(θ) ≈ R₀ + G × sin²θ + F × sin²θ × tan²θ, where R₀ is the intercept (near-angle reflection coefficient), G is the gradient (sensitivity to fluid and porosity), and F is the far-angle term; AVO crossplot analysis of R₀ versus G is the primary tool for classifying gas sand responses and distinguishing lithology from fluid effects in amplitude anomaly interpretation.
What Is the Reflection Coefficient?
Every seismic wave traveling through the earth is partially reflected and partially transmitted at every boundary where rock properties change. The fraction that is reflected — the reflection coefficient — depends entirely on the acoustic impedance contrast across that boundary. High contrast means a strong reflection; low contrast means most of the energy passes through without reflection; zero contrast means the wave continues as if no boundary exists.
The reflection coefficient's importance in seismic exploration stems from this direct connection to formation properties. When a seismic wave from a surface source travels downward and encounters the boundary between a shale and a sandstone, the reflection it generates carries the acoustic impedance contrast as its amplitude. If the sandstone contains gas instead of brine, its acoustic impedance is lower — and the reflection coefficient and seismic amplitude at that boundary change in a predictable way. By measuring seismic amplitude carefully, geophysicists can infer not just where boundaries are (the traditional use of seismic) but what kind of fluid fills the pores of the rock below the boundary.
This amplitude-based reservoir characterization capability — amplitude versus offset analysis, seismic inversion for acoustic impedance, direct hydrocarbon indicators — is built entirely on the physics of the reflection coefficient and its sensitivity to formation density, velocity, and pore fluid. The reflection coefficient is therefore not just a processing parameter but the fundamental link between seismic measurements at surface and the rock and fluid properties that petroleum geologists and reservoir engineers need to understand the subsurface.
Reflection Coefficient in Seismic Interpretation
Seismic impedance inversion converts a stacked seismic volume (amplitude vs. time) into an acoustic impedance volume (Z vs. depth or time) by deconvolving the reflection coefficient series from the estimated seismic wavelet and integrating to recover absolute impedance values calibrated by well log impedance constraints — model-based inversion, sparse-spike inversion, and simultaneous P-impedance/S-impedance inversion are the principal methods, each producing a different resolution and geological smoothness in the output impedance volume; the resulting impedance volume allows geologists and reservoir engineers to map reservoir quality, net-to-gross, and fluid content variations spatially between wells, extending the point measurements from wells across the inter-well space using the seismic data's lateral continuity.
AVO analysis uses the angle-dependent reflection coefficient behavior to separate fluid effects from lithology effects in seismic amplitude anomalies — gas sands typically show Class II or Class III AVO response where the reflection amplitude increases in magnitude with offset (increasing angle from near-vertical incidence to 35 to 45 degree incidence at far offsets), while brine-saturated sands and carbonates show flat or decreasing amplitude with offset; identifying AVO anomalies on pre-stack seismic gathers calibrated against known gas discoveries is the standard workflow for DHI (direct hydrocarbon indicator) analysis in clastic basins and supports well placement decisions in exploration and appraisal drilling programs.
Reflection Coefficient Across International Jurisdictions
Canada (AER / WCSB): WCSB seismic exploration for Montney, Duvernay, and Cardium tight plays uses amplitude and AVO analysis calibrated against well data to identify sweet spots within tight formations before drilling — reflection coefficient analysis from 3D seismic data helps identify high-porosity, high-organic-content intervals in the Duvernay shale play (where reflection coefficient correlates with TOC through density effects) and high-permeability zones in the Cardium tight oil play; AER requires that seismic data used in resource assessments supporting well license applications be properly processed and interpreted, with well-to-seismic ties documented in the geological assessment.
United States (API / BSEE): Gulf of Mexico deepwater seismic interpretation relies heavily on amplitude anomaly analysis using reflection coefficient principles — the "bright spots" that indicate gas sands in the shallow Pleistocene section above the salt are among the clearest reflection coefficient expressions of hydrocarbon presence in seismic data anywhere in the world; BSEE's offshore well permitting process requires geological documentation of the exploration or delineation target that typically includes seismic interpretation referencing amplitude anomalies and reflection coefficient analysis to justify the well's prospectivity. The US Geological Survey uses reflection coefficient analysis in regional basin analysis for resource assessment reports that inform federal leasing decisions.
Norway (Sodir / NORSOK): NCS exploration relies on high-quality 3D seismic amplitude and AVO analysis for hydrocarbon trap identification in the North Sea and Norwegian Sea — Sodir's exploration data portal provides access to extensive reprocessed 3D seismic datasets from the NCS that operators use for reflection coefficient analysis in exploration license applications; NCS seismic interpretation has contributed significantly to global advances in seismic amplitude interpretation methodology, particularly through the development of AVO interpretation techniques at Statoil (now Equinor) and its academic partners in the late 1980s and 1990s.
Middle East (Saudi Aramco): Saudi Aramco uses 3D seismic interpretation including reflection coefficient and impedance analysis to map Arab Formation reservoir heterogeneity and fluid contacts across the Ghawar field — the enormous areal extent of Ghawar (approximately 280 km long) makes seismic-based reservoir characterization between wells essential for managing the waterflood front advancing through the Arab D reservoir; Aramco's geophysics research group has developed proprietary methods for integrating reflection coefficient and AVO analysis with rock physics modeling specific to Arab Formation carbonate reservoir properties that improve impedance-to-porosity and fluid-fill predictions compared to generic methods applied without formation-specific calibration.