Zoeppritz Equations
The Zoeppritz equations are a set of mathematical relationships derived in 1919 by German geophysicist Karl Zoeppritz that describe how seismic wave energy is partitioned at a boundary between two rock layers with different elastic properties. When a seismic wave reaches such a boundary, a portion of its energy is reflected back toward the surface, a portion is transmitted into the lower layer, and energy may be converted between compressional (P-wave) and shear (S-wave) modes. The Zoeppritz equations calculate the amplitude ratios of each reflected and transmitted wave as a function of the angle of incidence at the boundary. These equations form the theoretical foundation of Amplitude Versus Offset (AVO) analysis, one of the most powerful seismic methods for detecting hydrocarbons and characterizing reservoir fluids before drilling.
Physical Basis and Input Parameters
At any interface between two rock units, the boundary conditions require continuity of both displacement and stress across the contact. Satisfying these conditions for an incident P-wave produces four output wave modes: reflected P-wave, reflected S-wave, transmitted P-wave, and transmitted S-wave. The Zoeppritz equations express the amplitude ratios (reflection and transmission coefficients) of each output mode as functions of six input parameters: the P-wave velocities (Vp1 and Vp2) and S-wave velocities (Vs1 and Vs2) on each side of the interface, and the bulk densities (rho1 and rho2) on each side. At zero angle of incidence (normal incidence), the P-wave reflection coefficient simplifies to the familiar expression (Z2 minus Z1) / (Z2 plus Z1), where Z = rho x Vp is acoustic impedance. As the angle increases, the interaction between all six parameters determines how reflectivity changes with offset, and this variation contains information about the S-wave velocity and density contrast that is invisible at zero offset.
AVO Analysis and Linearized Approximations
The full Zoeppritz equations are transcendental and must be solved numerically, which made them impractical for routine interpretation before digital computing. In the 1980s, several linearized approximations were developed for small-to-moderate angles of incidence. The Aki-Richards approximation expresses the P-wave reflection coefficient as a linear combination of relative contrasts in P-wave velocity, S-wave velocity, and density. The Shuey two-term approximation further simplifies this to R(theta) = R0 + G x sin^2(theta), where R0 is the zero-offset reflectivity (also called the intercept) and G is the AVO gradient. The intercept primarily responds to acoustic impedance contrast (controlled by Vp and density), while the gradient responds strongly to Vp/Vs ratio changes, which in turn are sensitive to fluid content. Gas-saturated sandstones typically have lower Vp/Vs ratios than brine-saturated sandstones, causing predictable AVO behavior that distinguishes them from wet sands and shales.
AVO Classification and Gas Sand Detection
AVO analysis classifies anomalies into four standard types based on the signs and magnitudes of the intercept and gradient. Class I anomalies have positive intercept (hard kick) and positive gradient; the reflection weakens with offset. Class II anomalies have near-zero intercept with either positive or negative gradient and can flip polarity with offset. Class III anomalies, the most common gas sand signature in Tertiary basins, have negative intercept and negative gradient, meaning reflectivity increases in magnitude (brightens) with offset. Class IV anomalies have negative intercept and positive gradient, where amplitude decreases with offset despite a bright zero-offset response. The AVO crossplot of intercept versus gradient is a standard interpretation tool: brine sands and shale-on-shale reflections cluster near the background trend, while gas sands plot away from this trend in characteristic quadrants depending on their class. Direct Hydrocarbon Indicators (DHIs) such as bright spots, flat spots, and dim spots on seismic data gain interpretive confidence when supported by AVO analysis consistent with the predicted Zoeppritz behavior for a given fluid substitution scenario.
Connection to Seismic Inversion and Rock Physics
The Zoeppritz equations link the observable seismic reflection response to the subsurface elastic properties (Vp, Vs, density) that control it. Seismic inversion uses this relationship in reverse: given the observed angle-dependent reflectivity from near, mid, and far offset stacks, inversion algorithms recover the elastic property contrasts at each reflector. Simultaneous pre-stack inversion produces volumes of P-impedance (Vp x rho), S-impedance (Vs x rho), and density, which can then be projected into lithology and fluid probability through rock physics templates calibrated to well log data. Gassmann fluid substitution is used to predict how the Zoeppritz-derived reflection response would change if brine were replaced by gas or oil, enabling prediction of DHI behavior at undrilled locations. Lambda-Mu-Rho (LMR) attributes derived from pre-stack inversion are a further transformation of the Zoeppritz-based elastic parameters that can separate fluid effects (lambda) from frame stiffness effects (mu) with greater clarity than Vp/Vs ratios alone.
Limitations and Practical Considerations
The Zoeppritz equations assume planar, isotropic, homogeneous layers separated by a flat interface. Real reservoirs violate all of these assumptions to varying degrees. Anisotropy, particularly vertical transverse isotropy (VTI) caused by laminated shales or horizontal fractures, alters the angle-dependent reflectivity in ways not captured by the isotropic Zoeppritz formulation. Curved interfaces, thin beds below seismic resolution, and wavelet tuning effects all complicate AVO analysis. Noise amplification is a significant practical problem because the gradient term in the Shuey approximation is estimated from the difference between near and far offsets, and both are noisy; gradient volumes therefore have lower signal-to-noise ratios than zero-offset stacks. Overburden anisotropy also introduces apparent AVO effects unrelated to the target reflector. Despite these limitations, AVO analysis grounded in Zoeppritz theory has a well-documented track record of improving drilling success rates when applied carefully with well control and quality-controlled offset gathers.
Key Takeaways
- The Zoeppritz equations describe how seismic reflection and transmission coefficients vary with angle of incidence at a rock boundary, using P-wave velocity, S-wave velocity, and density on both sides as inputs.
- The Shuey two-term linearization reduces the equations to an intercept (R0) and gradient (G), where the gradient is sensitive to Vp/Vs ratio changes caused by fluid substitution, enabling gas sand detection from surface seismic data.
- AVO Class III anomalies (negative intercept, negative gradient) are the classic bright-spot gas sand signature in Tertiary basins and represent the most widely used Zoeppritz-based DHI indicator.
- Pre-stack seismic inversion uses Zoeppritz-derived reflectivity as the forward model to recover elastic property volumes (P-impedance, S-impedance, density) that feed rock physics and reservoir characterization workflows.