Angle of Approach: Seismic Refraction, Head Waves, and Near-Surface Statics
The angle of approach is the acute angle formed between an incoming seismic ray and the reflecting or refracting interface at the precise point of intersection, measured from the interface plane itself rather than from the perpendicular to that interface. It is the geometric complement of the angle of incidence: because the interface plane and its normal are perpendicular, the two angles always sum to exactly 90 degrees, so angle of approach (θa) equals 90 degrees minus the angle of incidence (θi). A ray traveling vertically straight down onto a horizontal reflector has an angle of incidence of zero and an angle of approach of 90 degrees, meaning it strikes the interface perpendicular to it; a ray grazing nearly parallel to the interface has an angle of incidence close to 90 degrees and an angle of approach close to zero degrees. The term appears in three distinct but mathematically unified contexts within applied petroleum geophysics. In seismic reflection work, angle of approach is simply a geometric descriptor of ray orientation at any reflector and is rarely used independently because the angle of incidence is the conventional parameter in Zoeppritz-based amplitude analysis. In seismic refraction surveying, however, the angle of approach takes on specific operational importance: it is the emergent angle at which a critically refracted head wave arrives at a surface geophone, and this emergent angle directly encodes the seismic velocity of the refracting layer through Snell's Law, forming the basis for near-surface velocity models and the datum static corrections applied before reflection data are stacked. In acoustic borehole logging with full-waveform monopole or dipole sonic tools, the angle at which the compressional head wave impinges on the receiver array after propagating along the formation wall is also described as an angle of approach and governs the tool's transmitter-receiver spacing and frequency design. All three usages are consistent with Snell's Law and with the critical-angle condition at which refracted energy travels exactly along the interface, producing an angle of incidence equal to arcsin(V1/V2) and a corresponding angle of approach at the surface of 90 degrees minus that critical value.
Key Takeaways
- Geometric complement of angle of incidence: The angle of approach and the angle of incidence are always complementary, summing to 90 degrees at any interface regardless of that interface's dip or the complexity of the overlying velocity field. The distinction matters for communication rather than physics: geophysicists working in refraction interpretation describe arrival geometry using the angle of approach because it is the angle visually observable at the surface receiver, while reflection seismologists default to the angle of incidence because Snell's Law and the Zoeppritz equations are formulated relative to the normal. In depth-imaging workflows that integrate refraction-derived statics with reflection-derived velocity models, understanding both conventions prevents sign errors and unit confusion when comparing deliverables across contractors or software packages. Knowing that a head wave arriving at 68 degrees from horizontal represents an angle of incidence of 22 degrees allows the interpreter to directly verify consistency with the refractor velocity model without additional conversion steps.
- Head wave emergence and refractor velocity determination: In seismic refraction surveying, a critically refracted head wave travels along the top of a high-velocity refractor at that layer's P-wave velocity and continuously radiates energy back to the surface at the critical angle of approach. By measuring the inverse slope of the refraction travel-time curve (time versus offset), the interpreter recovers the refractor velocity V2 directly. Applying the complement relationship converts the measured emergent angle of approach to the angle of incidence, and Snell's Law (sin θc = V1/V2) then solves for V2 given V1 from the direct-wave slope. In the Western Canada Sedimentary Basin, near-surface refractors are typically the base of a glacial till or the top of a Paleozoic carbonate, with weathering-zone velocities of 400 to 1,200 m/s and sub-weathering velocities of 1,800 to 2,400 m/s, producing emergent angles of approach between 57 and 72 degrees at typical offsets of 200 to 600 metres.
- Near-surface static corrections from refraction geometry: Datum static corrections are time shifts applied to every seismic trace before stacking to remove the effect of irregular topography and the slow weathering zone and to project all receivers mathematically onto a smooth floating datum. The refraction component of this workflow uses first-break travel times, which are controlled by the head wave angle of approach at each receiver, to map the thickness and velocity of the weathering layer across the entire 2D line or 3D program. Delay-time methods, generalized reciprocal methods, and first-break tomography all invert these travel times with the angle of approach as the critical geometric observable. Errors in the refraction model translate directly into residual static anomalies that create false structural relief on horizon maps. In the Pembina Cardium fairway of west-central Alberta, near-surface static corrections routinely reach 40 to 120 milliseconds, representing potential structural misplacement of 60 to 180 metres in two-way time if the refraction model is incorrect.
- Acoustic borehole logging and formation head wave geometry: Full-waveform monopole sonic logging tools measure compressional and shear head waves that propagate from a transmitter along the borehole wall and arrive at receiver transducers spaced 0.5 and 1.0 metres above the transmitter. The head wave departs the transmitter at the critical angle into the formation, travels along the formation wall at the formation P-wave or S-wave velocity, and re-enters the borehole fluid at the reciprocal critical angle, impinging on the receiver array at an angle of approach that is determined by the formation velocity relative to the borehole fluid velocity. Higher formation velocities relative to fluid velocity produce a larger critical angle of incidence and thus a smaller angle of approach at the receiver. Dipole shear tools, which generate flexural-wave modes at 1 to 8 kHz rather than the 10 to 30 kHz used by monopole compressional tools, are required when the formation shear velocity is slower than the borehole fluid compressional velocity, a condition called slow formation, because no shear head wave can then exist at the borehole wall.
- Survey design and near-offset muting constraints: In 3D seismic reflection programs, refraction head waves recorded at short source-receiver offsets produce first arrivals that overprint and interfere with shallow reflection events on common-midpoint gathers. A top-mute function is applied to remove these refraction-contaminated near-offset traces before normal-moveout correction and stacking. The crossover distance at which the head wave first-break overtakes the direct wave, and therefore the minimum safe offset for reflection data, is determined by the angle of approach of the head wave and the two-way time to the shallowest target reflector. In the Montney fairway of northeast British Columbia, where targets lie at 2,500 to 3,500 metres depth and near-surface velocities average 1,900 m/s beneath thin till cover, crossover offsets typically fall between 800 and 1,400 metres, constraining the inner mute boundary on 3D gathers and limiting the minimum near-offset available for AVO analysis of the Montney A and B intervals.
Refraction Survey Processing and Static Correction Workflows
A seismic refraction survey acquires first-break travel times from a surface source to an array of geophones extending several hundred metres to several kilometres from the shot point. Processing begins by picking the first-break onset on each trace to within one millisecond accuracy, a task performed interactively in workstation software by identifying the inflection point where the trace leaves the noise floor and begins a coherent upswing. The time versus offset plot of first breaks reveals two or more linear segments with different slopes: the steep-slope near-field segment represents the direct wave traveling through the slow weathering zone, while the shallow-slope far-field segment represents the head wave that traveled ahead through the fast sub-weathering refractor. The inverse of the far-field slope equals the refractor velocity V2, and the intercept time of the far-field segment encodes the depth to the refractor at the mid-point of the survey. Converting the measured angle of approach of the far-field arrivals to the angle of incidence and applying Snell's Law cross-checks this velocity estimate independently.
Modern refraction processing has moved beyond simple two-layer slope-intercept analysis toward generalized reciprocal methods and seismic refraction tomography. The generalized reciprocal method separates individual receiver delay times from the overall velocity model, allowing the refractor to be mapped as a continuously undulating surface rather than a planar layer. This is important in glaciated terrains across the WCSB plains, where buried valleys filled with soft clay can accelerate the weathering zone depth from 10 metres at the valley margins to 80 metres at the valley axis over horizontal distances of 200 metres, creating static anomalies of 40 to 90 milliseconds concentrated in a narrow corridor. Seismic refraction tomography inverts all first-break times simultaneously using iterative ray-tracing through a starting velocity model, updating the model at each iteration until predicted and observed travel times converge. Cell sizes of five metres vertically and 20 metres horizontally resolve fine-scale velocity heterogeneity that the reciprocal method would smear into a smooth model.
The final deliverable from refraction processing is a gridded near-surface velocity and thickness model used by the datum static workflow to calculate a scalar time shift at each surface receiver location. These datum statics are applied to every trace recorded at that location before the traces are sorted into common-midpoint gathers and subjected to normal-moveout correction, deconvolution, and migration. Residual static corrections, estimated iteratively from the consistency of reflection moveout within CMP gathers, address the short-wavelength component of the static anomaly not captured by refraction modeling. However, residual statics can only correct anomalies smaller than one-quarter of the dominant wavelength of the reflection data; larger anomalies must be captured correctly by the refraction model or they persist as structural artifacts in the final migrated image.
In the Peace River Arch area of north-central Alberta, thick Quaternary drift deposits overlie tight Keg River and Slave Point carbonate reservoirs at depths of 1,400 to 2,000 metres. Refraction surveys in this area routinely map two sub-surface velocity boundaries: the base of active weathering at 10 to 25 metres depth with velocities of 450 to 750 m/s, and the base of the till at 20 to 65 metres depth with velocities of 1,600 to 2,000 m/s. A 160-channel 2D refraction line at 10-metre receiver spacing and maximum shot offset of 1,600 metres requires 10 to 14 days of field acquisition at a total cost of CAD 90,000 to CAD 130,000 including crew, recording equipment, explosive charges, and environmental access fees. The resulting static corrections, validated at existing well ties, reduce the standard deviation of the interpreted top-Keg River two-way time from 14 milliseconds to 4 milliseconds across the program area, improving structural closure mapping accuracy by 25 to 50 metres and enabling drill location ranking that would otherwise be unreliable without the refraction-derived static model.
Fast Facts
Angles of approach in WCSB near-surface refraction surveys typically range from 55 to 80 degrees above horizontal, reflecting the strong velocity contrast between glacial till (400 to 1,500 m/s) and underlying Paleozoic carbonates or crystalline basement (1,800 to 3,500 m/s). First breaks usable for angle-of-approach measurement appear at offsets starting 100 to 300 metres from the shot point, depending on weathering-zone thickness, which ranges from 5 to 100 metres across the Alberta plains. Full-waveform sonic tools designed around head-wave geometry measure compressional slowness values from 140 microseconds per metre in tight Montney dolomite to 380 microseconds per metre in organic-rich shale. Datum static corrections derived from refraction surveys in heavily glaciated terrain, such as the Horseshoe Canyon coalbed methane play south of Drumheller, routinely exceed 80 milliseconds and represent the single largest source of structural uncertainty in shallow gas exploration.