The plus-minus method
It often is difficult to use first breaks to estimate the intercept time and velocities for the weathering layer and bedrock. This is primarily because the base of weathering typically is an undulating surface, which makes traveltime plots difficult to interpret. Traveltime plots also are affected by severe elevation changes. Additionally, a typical field cable layout does not provide a sufficient number of channels inside the crossover distance xc (Figure 3.4-11a) for a reliable estimate of the weathering velocity or thickness. In most cases, vw cannot be measured and a reasonable value is assumed for it.
 formulated a method to indirectly estimate intercept time and bedrock velocity. The method still requires picking the first breaks. However, it does not require interpreting the traveltime profile (Figure 3.4-11a). (Interpretation means drawing the linear segments for the direct and refracted waves.) Figure 3.4-12a shows three raypaths associated with shot-receiver pairs AD, DG, and AG. The basis of Hage-doorn’s method involves computing two time values, the plus and minus times given by
Figure 3.4-11 (a) Geometry for refracted arrivals. Here, vw = weathering velocity, vb = bedrock velocity, zw = depth to the refractor equivalent to the base of the weathering layer, θc = critical angle, and xc = crossover distance. The direct wave arrival has a slope equal to 1/vw and the refracted wave arrival has a slope equal to 1/vb. (b) A shot record that exhibits the direct wave and the refracted wave depicted in (a). (c) Geometry for a dipping refractor with forward traveltime profile associated with the direct wave and refracted wave arrivals, and (d) with both forward and reverse traveltime profiles. See text for details.
Figure 3.4-12 (a) Geometry for the plus-minus method. (b) Geometry for the generalized reciprocal method. Here, zw is the depth to the refractor at the surface station where the plus-minus times as for (a) and intercept times as for (b) are to be estimated, vw is the weathering velocity, and θc is the critical angle of refraction.
The times given on the right side of these equations are the measured (picked) values from the first breaks for the three raypaths shown in Figure 3.4-12a. From the raypath configuration, we find that (Section C.7)
Rewrite equation (41a)
and note that the plus time t+ in equation (48a) is identical to the intercept time ti in equation (48b). Hence, instead of measuring ti directly from the shot record, Hagedoorn’s method suggests estimating ti from the first break picks given by the right-hand terms in equation (47a).
By applying algebra (Section C.7), we find that minus time t− is related to bedrock velocity vb by
where x is the source-receiver separation AD.
Thus, Hagedoorn’s plus-minus method involves:
- Picking the first breaks,
- Computing the plus-minus times, t− and t+ (equations 47a and 47b),
- Deriving from the plus-minus times the intercept time ti (equations 48a and 48b) and bedrock velocity vb (equation 48c),
- Assuming a value for weathering velocity vw,
- Computing the depth zD to bedrock below station D (Figure 3.4-12a) from equation (41a), and
- Computing the shot-receiver static shift ΔτD at that station by
where ES and ED are the surface and datum elevations at station D (Figure 3.4-12a). If there is a shot at station D, ΔτD represents the shot static, and if there is a receiver at station D, it represents the receiver static. Again, uphole and elevation corrections are needed before making the plus-minus statics corrections.
- First breaks
- Field statics corrections
- Flat refractor
- Dipping refractor
- The generalized reciprocal method
- The least-squares method
- Processing sequence for statics corrections
- Model experiments
- Field data examples
- Topics in moveout and statics corrections
- Hagedoorn (1959), Hagedoorn, J. G., 1959, The plus-minus method of interpreting seismic refraction sections: Geophys. Prosp., 7, 158–182.