OilAuthority.com

Snell's Law: Definition, Seismic Refraction, and AVO Applications

What Is Snell's Law?

Snell's law states that the ray parameter p = sin(θ)/V is conserved as a seismic wave crosses an interface between two formations with different velocities, governing how P-waves and S-waves refract, reflect, and mode-convert at subsurface boundaries — making it the foundational equation for seismic ray tracing, velocity model building, refraction statics corrections, and AVO analysis across all oil and gas exploration programmes.

Key Takeaways

  • The conservation of the ray parameter p = sin(θ)/V at each interface means that a wave going from a slow formation into a faster one bends away from vertical (increases angle), while a wave entering a slower formation bends toward vertical.
  • The critical angle (θc = arcsin(V1/V2) for V2 > V1) is where the refracted wave travels along the interface, generating head waves used as first arrivals in refraction surveys for near-surface velocity model building.
  • At non-normal incidence, P-waves generate both reflected P-waves and converted S-waves; the angle of the converted wave is set by requiring p = sin(θP)/VP = sin(θS)/VS simultaneously, producing a shallower conversion angle because VS ≈ VP/2.
  • AVO analysis uses Snell's law geometry through the Zoeppritz equations to predict how reflection amplitude varies with offset, enabling gas sand identification in the Permian Basin (US), Montney (Canada), North Field (Qatar), and Johan Sverdrup (Norway).
  • Refraction seismic surveys applying Snell's law to first-break arrivals are used for statics corrections and near-surface velocity models in land seismic programmes worldwide, including AER-regulated 3D surveys in Alberta.

How Snell's Law Works

For a seismic wave crossing an interface between a layer with velocity V1 and a layer with velocity V2, Snell's law states: sin(θ1)/V1 = sin(θ2)/V2, where θ1 is the angle of incidence and θ2 is the angle of refraction, both measured from the normal to the interface. The conserved quantity p = sin(θ)/V is called the ray parameter. In a layered earth where velocity increases with depth, the ray parameter is constant throughout the entire ray path, causing rays to curve progressively away from vertical as they enter faster formations — the basis for the curved raypath geometry in refraction surveys and in turning-wave tomography used to build long-wavelength velocity models.

Mode conversion at an interface follows Snell's law simultaneously for all generated wave types. An incident P-wave generates reflected P, transmitted P, reflected S (converted), and transmitted S (converted) at the same ray parameter p. The converted S-wave travels at the S-wave velocity VS, so its angle is arcsin(p × VS), which is roughly half the P-wave angle because VS ≈ VP/2. This geometry enables long-offset converted-wave (C-wave) acquisition where the mode-converted S-wave leg covers a wider range of depths than the P-wave leg from the same offset.

Snell's Law Across International Applications

In Canada, refraction seismic surveys applying Snell's law are routinely acquired ahead of or jointly with reflection surveys in Alberta and British Columbia to establish the near-surface velocity model required for static corrections; AER Directive 082 governs seismic operations including acquisition geometry that determines how well the refraction geometry samples the critical angle arrivals. The Montney and Duvernay 3D programmes all use refraction-derived statics as a foundation for reflection processing.

In the United States, BOEM-permitted seismic programmes on the Outer Continental Shelf use refraction analysis for water-bottom velocity characterisation in deepwater surveys. Onshore, the Permian Basin, Eagle Ford, and Marcellus shale programmes use full-waveform inversion (FWI) — a computationally intensive extension of Snell's law ray-tracing — to build high-resolution near-surface velocity models. In Norway, Sodir's data requirements include processed seismic velocity models; the wide-azimuth surveys over the Johan Sverdrup Field use Kirchhoff migration governed by Snell's law ray-tracing in an anisotropic velocity model. In the Middle East, refraction surveys over the Arabian Shield establish the overburden velocity model for deep Khuff and Arab formation targets in Saudi Aramco and ADNOC concessions. Australia's NOPSEMA-regulated surveys in the Carnarvon and Cooper basins use refraction first-break tomography as standard pre-processing.

Fast Facts

Snell's law was formulated in the context of light refraction by Dutch mathematician Willebrord Snel in 1621 — over three centuries before it was applied to seismic exploration — yet remains the governing equation for every seismic migration algorithm and velocity model used in oil and gas exploration today.

Snell's Law in AVO and Migration

AVO (amplitude variation with offset) analysis depends directly on Snell's law through the Zoeppritz equations, which compute the amplitude of each wave type generated at an interface as a function of incidence angle, and the two-term Shuey approximation commonly used in AVO attribute extraction. The incidence angle at the reservoir for a given source-receiver offset is computed by ray-tracing through the overburden velocity model using Snell's law at every interface. An error in the velocity model produces incorrect incidence angles and incorrect AVO attributes, which can cause a gas sand to be misidentified as a wet sand — the most consequential interpretation error in seismic exploration.

Kirchhoff depth migration, the most widely used seismic migration algorithm for complex structures, implements Snell's law by computing diffraction traveltimes through the velocity model using ray tracing. Phase-shift migration in the frequency-wavenumber domain is a wave-equation method that is equivalent to applying Snell's law continuously through a laterally varying velocity field. Both methods require an accurate velocity model — typically built by iterative tomography that minimises residual moveout, itself a measure of how well the velocity model satisfies Snell's law at every recorded offset.

Tip: In shallow refraction surveys for statics corrections, the critical angle geometry requires that receiver offsets extend to at least twice the depth of the refractor being mapped; insufficient offset means the head wave never arrives as a first break and the refractor velocity is not sampled, resulting in an incomplete statics solution and residual long-wavelength velocity errors that degrade reflection stack quality.

Snell's law is also known as:

  • Snell-Descartes law — the French term honouring both Snel and René Descartes, who independently derived the refraction relationship
  • Law of refraction — the generic descriptive term used in introductory physics and seismology texts
  • Ray parameter conservation — the form of the law expressed in terms of the invariant p = sin(θ)/V, used in advanced seismological texts

Related terms: P-wave, S-wave, AVO, Kirchhoff migration, VSP, Fermat's principle

Frequently Asked Questions

What is Snell's law in seismic exploration?

In seismic exploration, Snell's law states that the ray parameter (sin θ / V) is conserved as a seismic wave crosses an interface between formations of different velocities. It governs how seismic rays bend when entering faster or slower rock layers, underpins refraction statics corrections, controls AVO incidence angle geometry, and is the basis for all seismic migration algorithms used in subsurface imaging.

What is the critical angle in seismic refraction?

The critical angle is the angle of incidence at which the refracted wave travels along the interface rather than transmitting into the lower formation. It is given by θc = arcsin(V1/V2), where V1 is the upper layer velocity and V2 is the lower (faster) velocity. Waves at the critical angle generate head waves — first arrivals recorded at long offsets in refraction surveys used to map near-surface velocity structure.

How does Snell's law relate to AVO analysis?

AVO analysis uses the Zoeppritz equations, which require the angle of incidence at the reservoir for each source-receiver offset. These angles are computed by ray-tracing through the overburden velocity model using Snell's law at every interface. An inaccurate velocity model produces wrong incidence angles and distorted AVO attributes, potentially misidentifying the pore fluid in the reservoir.

Why Snell's Law Matters in Oil and Gas

Snell's law is the physical foundation of every seismic workflow used in oil and gas exploration and development — from the simple refraction statics correction that removes near-surface noise from a land seismic record to the complex full-waveform inversion that builds the velocity model for a deepwater pre-salt image. Every AVO-driven drill decision, every depth-migrated structural map, and every 4D reservoir monitoring result depends on the correct application of Snell's law. It is, in practical terms, the equation that connects surface seismic observations to the rock properties and fluid contacts below the drill bit.

Browse the Oil and Gas Glossary