Angle of Incidence: Definition, Snell's Law, and AVO Analysis

The angle of incidence is the acute angle formed between an incoming seismic ray and the normal (perpendicular) to a reflecting or refracting interface at the point of intersection. It is conventionally measured in degrees from the normal, not from the interface itself; the complement of the angle of incidence is the angle of approach, which is measured from the interface plane. Understanding the angle of incidence is central to nearly every quantitative seismic discipline, from basic ray-path geometry and refraction surveying to amplitude-versus-offset (AVO) analysis, full-waveform inversion, and acoustic borehole logging. When a seismic wave strikes an interface at any angle other than exactly perpendicular, the relationship between incoming and outgoing energy is determined by the angle of incidence through Snell's Law and the full set of Zoeppritz equations, and the balance between reflected and transmitted energy shifts dramatically as the angle increases toward and beyond the critical angle.

Key Takeaways

• The angle of incidence (theta-i) is measured between the incoming ray and the normal to the interface; it equals zero for a wave traveling straight down and 90 degrees for a wave traveling along the interface.
• Snell's Law governs both reflection (angle of reflection equals angle of incidence) and refraction (sin theta-i / V1 = sin theta-t / V2), directly linking the angle of incidence to the velocity contrast across the interface.
• The critical angle (theta-c = arcsin(V1/V2)) is reached when the transmitted ray travels exactly along the interface; beyond this angle, total internal reflection occurs and all energy is reflected, a phenomenon exploited in refraction seismic surveys.
• AVO analysis relies on the systematic variation of reflection amplitude with angle of incidence to infer rock properties such as porosity, fluid content, and acoustic impedance contrast; near, mid, and far angle stacks (roughly 0-15 deg, 15-30 deg, 30-45 deg) are routinely extracted in seismic processing.
• In acoustic borehole logging, the angle of incidence at the borehole wall determines whether formation compressional and shear head waves are excited, directly controlling the measurement of formation slowness in acoustic log tools.

Fundamental Definition and Geometry

In classical wave physics and seismic ray theory, an interface is any surface across which elastic properties change: a bedding plane separating shale from sand, the top of a salt body, the boundary between gas-saturated and brine-saturated rock, or the wall of a borehole. When a ray traveling through medium 1 reaches such an interface, it simultaneously generates a reflected ray (returning into medium 1) and a transmitted (refracted) ray continuing into medium 2. The angles of all three rays are referenced to the normal to the interface, an imaginary line perpendicular to the interface at the point of contact. The angle between the incident ray and this normal is theta-i (the angle of incidence), and by Snell's Law the reflected ray leaves at the same angle on the other side of the normal (theta-reflected = theta-i). The transmitted ray bends toward or away from the normal according to the velocity ratio.

In vertical seismic profiling (VSP), where geophones are deployed in a wellbore and a surface source generates waves, the angle of incidence at each reflector changes as a function of source offset; analyzing these angle-dependent amplitudes is the cornerstone of VSP-based AVO and impedance inversion. In conventional surface seismic reflection surveys, each recorded trace corresponds to a reflection from a specific source-receiver offset, and that offset maps to a specific angle of incidence at each reflector depth through the velocity model of the overburden. Converting from offset to angle is therefore an essential step in any AVO workflow.

The relationship between angle and offset is not linear; it depends on the velocity field through the Dix (NMO) conversion. For a flat reflector at depth z in a homogeneous medium with velocity V, the angle of incidence theta at offset x is:

tan(theta) = x / (2z),   or equivalently   sin(theta) = x / sqrt(x^2 + 4z^2)

In layered media, the apparent angle must be computed using the ray-parameter approach through the full interval velocity model. This offset-to-angle conversion is performed during seismic migration or post-migration in the angle-gather domain, producing angle gathers in which each trace represents reflections at a specific incidence angle rather than a specific offset.

Snell's Law and Wave Propagation at Interfaces

Snell's Law describes the relationship between the angles and velocities on both sides of a planar interface for any wave type. For P-to-P transmission (the most common case in conventional seismic reflection surveys):

sin(theta-i) / V1 = sin(theta-t) / V2

where V1 is the P-wave velocity in the incident medium and V2 is the P-wave velocity in the transmitted medium. Rearranging: theta-t = arcsin(V2 / V1 * sin(theta-i)). If V2 > V1 (a velocity increase downward, which is the normal case for increasing depth and compaction), the transmitted ray bends away from the normal (theta-t > theta-i). If V2 < V1 (a velocity decrease, as occurs at a gas-charged sand or at an unconsolidated sediment below a hard carbonate), the transmitted ray bends toward the normal (theta-t < theta-i). The critical condition is reached when theta-t equals 90 degrees, meaning the transmitted ray travels exactly along the interface. Setting sin(theta-t) = 1 in Snell's Law gives the critical angle:

theta-c = arcsin(V1 / V2)

Note that a critical angle exists only if V2 > V1. For incidence angles beyond theta-c, the wave is totally internally reflected: all incident energy returns to medium 1, none is transmitted. In refraction seismic surveys, the wave that travels along the interface at the critical angle is called a head wave (or refracted wave or first-break wave); it continuously re-radiates energy back into medium 1 at the critical angle, and these arrivals, which precede direct and reflected waves at large source-receiver offsets, are used to map shallow velocity structure. For example, if V1 = 1,800 m/s (5,900 ft/s) in a near-surface layer and V2 = 3,500 m/s (11,480 ft/s) in a deeper consolidated formation, the critical angle is arcsin(1,800/3,500) = 30.9 degrees. At a receiver offset of 200 m from a 50 m deep refractor, the angle of incidence is approximately 63.4 degrees, well past critical, so the refracted first-break arrival will be observed.

At elastic interfaces, Snell's Law must be applied simultaneously to P-to-P, P-to-S, S-to-P, and S-to-S wave conversions, using the appropriate velocities:

sin(theta-P1) / VP1 = sin(theta-S1) / VS1 = sin(theta-P2) / VP2 = sin(theta-S2) / VS2 = ray parameter p

The single conserved quantity is the ray parameter (also called horizontal slowness) p, which is preserved across all boundaries along a ray path. Because VS is always less than VP for the same rock, a P-wave converting to S at an interface will always produce an S-wave transmitted at a smaller angle than the P-wave component; the P-to-S critical angle is also larger (or may not exist) compared to P-to-P.

AVO Analysis and the Zoeppritz Equations

Amplitude variation with offset (AVO) analysis, also termed amplitude variation with angle (AVA), is the single most commercially important application of angle-of-incidence physics in the modern oil and gas industry. The Zoeppritz equations (1919), later simplified by Bortfeld (1961), Aki and Richards (1980), and Shuey (1985), give the exact reflection and transmission coefficients for P and S waves as a function of angle of incidence and the elastic properties on both sides of an interface. For a P-to-P reflection, the Shuey two-term linearized approximation is:

R(theta) = R0 + G * sin^2(theta)

where R0 is the zero-angle (normal-incidence) reflection coefficient, G is the AVO gradient, and theta is the angle of incidence. R0 is controlled mainly by the acoustic impedance contrast (product of density and velocity), while G is controlled by the contrast in Poisson's ratio (equivalently, the VP/VS ratio). A gas-saturated sand, which has a lower Poisson's ratio than the surrounding shale, typically shows a negative R0 and a strongly negative G, meaning the reflection amplitude becomes increasingly negative with increasing angle. This is the classic Class III or "bright spot" AVO response that has been used to identify gas sands since the mid-1970s. A brine-saturated sand may show a similar R0 but a weaker or positive G, allowing AVO analysis to discriminate gas from brine in many cases.

In seismic processing, AVO analysis is performed on angle stacks. Processors extract near-angle stacks (0-15 degrees), mid-angle stacks (15-30 degrees), and far-angle stacks (30-45 degrees) by sorting migrated gathers into angle bins and stacking. The far stack has the most AVO sensitivity but also the most contamination from wide-angle noise, residual NMO stretching, and amplitude instabilities. A mute function based on maximum incidence angle is applied before stacking to exclude traces beyond a cutoff angle (commonly 45-50 degrees) where wide-angle reflections, mode conversions, and head-wave energy degrade the signal. In depth-migrated angle gathers, the offset-to-angle conversion is performed explicitly using the migration velocity model, producing angle gathers that are more accurate than those from NMO-based offset-to-angle conversion.

The three-term Shuey approximation adds a curvature term:

R(theta) = R0 + G * sin^2(theta) + F * (tan^2(theta) - sin^2(theta))

where F is related to the contrast in P-wave velocity and contributes mainly at large angles (beyond 30 degrees). The curvature term improves fitting at wide angles but is also more sensitive to noise, so it is typically used only when high-quality far-angle data are available. See also: AVO, acoustic impedance.