Acoustic Impedance Section: Definition, Inversion, and Seismic

An acoustic impedance section is a two-dimensional or three-dimensional seismic data volume that has been mathematically transformed, through a process called seismic inversion, from its native form as a record of reflection amplitudes into a spatial representation of acoustic impedance values. Rather than showing where seismic waves bounced back to the surface, an acoustic impedance section shows the absolute physical property of each layer, enabling geoscientists to directly compare seismic data with well-log measurements and to identify lithology changes, fluid effects, and reservoir quality across a survey area. Sonic and density logs from nearby wellbores are used to calibrate the inversion and to supply the low-frequency content that seismic data alone cannot recover.

Key Takeaways

• Acoustic impedance (Z) is the product of compressional-wave velocity (Vp) and bulk density (rho), measured in kg/m²/s (rayl); typical values range from approximately 4 × 10&sup6; rayl for soft marine sediments to more than 20 × 10&sup6; rayl for tight carbonates.
• The reflection coefficient at a boundary is R = (Z2 - Z1) / (Z2 + Z1); seismic inversion reverses this relationship to recover Z1 and Z2 from observed reflectivity R.
• Seismic data is bandlimited: low frequencies below roughly 8 Hz must be supplied from well logs, while high frequencies above 80-100 Hz are attenuated by the earth and cannot be recovered. This makes log calibration mandatory.
• Different inversion algorithms, including recursive inversion, model-based inversion, sparse-spike inversion, and geostatistical inversion, each offer different trade-offs among resolution, noise sensitivity, and computational cost.
• Simultaneous inversion extends the acoustic impedance approach to jointly solve for P-impedance (Ip), S-impedance (Is), and density (rho), enabling rock-physics crossplots that discriminate lithology from fluid effects.

How Acoustic Impedance and Seismic Reflectivity Are Related

Every seismic reflection that a geophysicist sees on a conventional seismic section originates at an interface where acoustic impedance changes. The reflection coefficient R at a planar boundary between layer 1 (impedance Z1) and layer 2 (impedance Z2) is given by the expression R = (Z2 - Z1) / (Z2 + Z1). When Z2 is greater than Z1, R is positive and the reflected wavelet has the same polarity as the incident wave. When Z2 is less than Z1, R is negative and the wavelet is reversed. What a seismic section actually records is the convolution of the earth's reflectivity series with a seismic wavelet, plus noise. The wavelet blurs adjacent reflections together and imposes the bandlimited frequency content described above, making direct lithological or fluid interpretation difficult from raw amplitudes alone.

Seismic inversion addresses this problem by deconvolving the wavelet effect, recovering a broadband estimate of the impedance contrast series, and then integrating that series upward from a known reference level to produce an absolute impedance volume. The result is a dataset whose display is directly analogous to a wireline log curve repeated laterally across the entire seismic survey. A geoscientist can place a synthetic log computed from the inversion at any point and compare it immediately with measured sonic and density data from a nearby well, providing a powerful QC workflow. Gas sands, which have low velocity and low density, yield characteristically low acoustic impedance values, typically below 5 × 10&sup6; rayl in Tertiary clastic basins, and appear as prominent low-AI anomalies on the inverted section. Tight, water-saturated sands or carbonates yield higher impedance values and appear bright on a reversed-polarity impedance display.

The Model-Based Inversion Workflow

The most widely applied class of seismic inversion in commercial exploration is model-based inversion, also known by proprietary algorithm names such as BLIMP (Band-Limited Impedance from Model Parameters). The workflow proceeds in three major stages. First, the interpreter performs a well tie, extracting or statistically estimating the dominant seismic wavelet by cross-correlating a synthetic seismogram (computed from acoustic log and density log data) with the actual seismic trace at the well location. The quality of the well tie, typically expressed as a cross-correlation coefficient, sets the upper bound on inversion quality; ties below 0.7 are generally considered insufficient for reliable quantitative inversion. Second, a low-frequency model is constructed by interpolating smoothed impedance logs between all available wells, guided by seismic horizons that control the lateral geometry of each stratigraphic interval. This model provides the sub-8 Hz frequency content that the seismic data cannot contain. Third, the inversion optimization runs iteratively, perturbing the impedance model until the synthetic seismogram computed from the model matches the actual seismic data within a user-defined misfit tolerance, while simultaneously penalizing large departures from the low-frequency model through a regularization term.

Sparse-spike inversion, such as the CSSI (Constrained Sparse-Spike Inversion) algorithm, takes a different philosophical approach. Rather than continuously updating a smooth background model, it searches for the minimum number of large impedance contrasts that can reproduce the observed seismic data. This approach yields sharper, blocky impedance profiles that more closely resemble the abrupt layer boundaries seen in well logs. It is particularly effective in thin-bed environments where model-based approaches tend to produce smoothly varying impedance profiles that smear distinct reservoir layers together. Geostatistical inversion methods, which combine co-simulation algorithms such as Sequential Gaussian Simulation (SGS) with seismic data conditioning, generate multiple equally probable realizations of the subsurface impedance volume. This enables uncertainty quantification and probabilistic resource estimation, which is increasingly demanded by operators and regulatory bodies when booking reserves.

Calibration with Sonic and Density Logs

The acoustic impedance section derives its quantitative value entirely from its calibration to measured well data. The acoustic log (sonic log) measures the compressional-wave travel time in microseconds per foot or microseconds per metre through the formation adjacent to the borehole, from which interval velocity is computed. The density log measures bulk density in g/cm³ using a gamma-gamma backscatter tool. Multiplying these two curves sample by sample produces the impedance log Z(depth) that anchors the inversion. In wells where the density log is absent or of poor quality, empirical rock-physics relationships such as Gardner's equation (rho = a × Vp^b, where typical constants are a = 0.31 in SI units, b = 0.25) may be used to estimate density from velocity, though this introduces additional uncertainty.

Before inversion, the impedance log must be blocked to the seismic resolution scale, typically 10-30 m depending on dominant frequency and depth, because sub-resolution log heterogeneity cannot be recovered from seismic data and will introduce spurious misfit if included in the wavelet estimation. The blocked log is then used in the well tie, the low-frequency model construction, and post-inversion QC. A blind well test, in which one well is withheld from the inversion and its impedance log is compared against the inverted volume at that location, is the industry standard for validating that the inversion has genuinely transferred information from the seismic data rather than simply interpolating between the calibration wells.