# Interpretation of four-shot refraction data

Series Geophysical References Series Problems in Exploration Seismology and their Solutions Lloyd P. Geldart and Robert E. Sheriff 14 497 - 503 http://dx.doi.org/10.1190/1.9781560801733 ISBN 9781560801153 SEG Online Store

## Problem 14.3

Find the velocities and depths of refractors in Figure 14.3b.

### Background

The four-shot method uses a refraction spread with a source at each end to provide “short-shot” data and two sources offset inline beyond the critical distance (see problem 4.18) to provide “long-shot” data that permit accurate measurement of the refraction velocities. A refraction event appearing on both profiles can be identified by the slopes of the ${\displaystyle t}$-${\displaystyle x}$ curves. Geophones at the shot points are offset when shots at ${\displaystyle A}$ or ${\displaystyle B}$ are used.

### Solution

Figure 14.3b shows east-west reversed profiles with sources ${\displaystyle A}$ and ${\displaystyle B}$ providing the short-shot (SS) data and sources ${\displaystyle A^{*}}$ and ${\displaystyle B^{*}}$ (not shown) that are offset inline to provide the long-shot (LS) data. We have labeled the different segments from ${\displaystyle C}$ to ${\displaystyle J}$ and fitted each with a straight line whose slope gives the apparent velocity. The SS segments were used to get intercept times as shown in the figure. Note that fitting straight lines to the data assumes that the refractors are planar, clearly an approximation.

The data suggest that this is a three-layer situation with dip down to the east. We note that the LS data segments ${\displaystyle (I,J)}$, ${\displaystyle (E,F)}$, and ${\displaystyle (C,D)}$ show a time-offset between geophones 6 and 7, as if indicating a fault downthrown to the east. However, evidences for features at depth (such as a fault cutting the refractor) should be displaced away from the source (because the head-wave path from the refractor to a geophone is inclined, as in Figure 14.3a). However profiles ${\displaystyle I,J}$ and ${\displaystyle E,F}$ show the same anomaly at the same location even though shot from opposite directions. Because our geophones spacing is only 20 m, this suggests that they are caused by something shallower than 20 m, and the most probable cause is a near-surface anomaly—perhaps a statics error in near-surface velocity or elevation corrections. We could make an empirical surface correction, but have not done so.

Figure 14.3a.  A feature on a refractor would be observed at ${\displaystyle D}$ from source ${\displaystyle A}$ but at ${\displaystyle C}$ from source ${\displaystyle B}$.
Figure 14.3b.  Time-distance curves for four-shot data. Open squares and triangles represent short-shot data, solid squares and triangles are long-shot data. (after Milson, 1989, 173)

The SS data show a shallow refraction event (segments ${\displaystyle C}$ and ${\displaystyle G}$) that probably is a refraction from the top of an intermediate layer. The measured velocities are ${\displaystyle V_{1}=700\ {\rm {m/s}}}$, ${\displaystyle V_{A2}=1030\ {\rm {m/s}}}$, ${\displaystyle V_{B2}=1470\ {\rm {m/s}}}$ (note that neither ${\displaystyle C}$ nor ${\displaystyle G}$ is well determined, both being based on two points only), with intercept times ${\displaystyle t_{A1}=25\ {\rm {ms}}}$, ${\displaystyle t_{B1}=77\ {\rm {ms}}}$. Using equation (4.24f) we get

{\displaystyle {\begin{aligned}V_{2}&=2(1/1030+1/1470)^{-1}=1210\ {\rm {m/s}};\\\theta _{c}&=\sin ^{-1}(700/1210)=35.3^{\circ };\\h_{A}&=V_{1}t_{d1}/2\cos \theta _{c}=700\times 0.025/2\times 0.816=11\ {\rm {m}};\\h_{B}&=700\times 0.077/2\times 0,816=33\ {\rm {m}};\\\xi _{1}&=\tan ^{-1}[(33-11)/220]=5.7^{\circ }.\end{aligned}}}

Because of the small dip, we can take ${\displaystyle h_{A}}$ and ${\displaystyle h_{B}}$ as vertical depths.

Comparison of velocities shows that segments ${\displaystyle E}$ and ${\displaystyle F}$ are the LS equivalent of SS segment ${\displaystyle D}$, while ${\displaystyle I}$ and ${\displaystyle J}$ are equivalent to ${\displaystyle H}$. To interpret these events we shall use the LS velocities and the SS intercept times.

Based on our previous discussion we ignore the displacements of segments ${\displaystyle E}$ and ${\displaystyle F}$ and ${\displaystyle I}$ and ${\displaystyle J}$ and draw the best-fit straight lines as shown in Figure 14.3b, obtaining ${\displaystyle V_{3d}=1860\ {\rm {m/s}}}$, ${\displaystyle V_{3u}=5240\ {\rm {m/s}}}$. We use these values plus the intercept times ${\displaystyle t_{A}=40\ {\rm {ms}}}$ and ${\displaystyle t_{B}=165\ {\rm {ms}}}$ to get the depth and dip.

We use the method of problem 4.24b to get ${\displaystyle \theta _{c2}}$ and ${\displaystyle \xi _{2}}$:

{\displaystyle {\begin{aligned}V_{2}/V_{3d}=\sin(\theta _{c2}+\xi _{2})=1210/1860,\quad (\theta _{c2}+\xi _{2})=40.6^{\circ };\\V_{2}/V_{3u}=\sin \theta _{c2}-\xi _{2})=1210/5240,\quad (\theta _{c2}-\xi _{2})=13.4^{\circ };\end{aligned}}}

so ${\displaystyle \theta _{c2}=27.0^{\circ }}$ , ${\displaystyle \xi _{2}=13.6^{\circ }}$, ${\displaystyle V_{3}=1210/\sin 27.0^{\circ },=2670\ {\rm {m/s}}}$.

Now ${\displaystyle \sin \theta _{1}/700=\sin \theta _{c2}/1210=1/2670}$, ${\displaystyle \theta _{1}=15.2^{\circ }}$, ${\displaystyle \theta _{c2}=26.9^{\circ }}$ (the fact that this value of ${\displaystyle \theta _{c}}$ agrees so closely with the above value is because both values depend upon the velocities).

The SS intercept times of 40 and 165 ms give the following depths:

{\displaystyle {\begin{aligned}{\mbox{at }}A:0.040&=2\times 11\cos 15.2^{\circ }/700+2\ h_{A2}\cos 26.9^{\circ }/1210\\&=0.030+0.00147\ h_{A2},\quad h_{A2}=6.8\ {\rm {m}};\\{\mbox{vertical depth}}&=11+7/\cos 13.6^{\circ }=18\ {\rm {m}}\\{\mbox{at }}B:0.165&=2\times 33\cos 15.2^{\circ }/700+2h_{B2}\cos 26.9^{\circ }/1210\\&=0.091+0.00147\ h_{B2},\quad h_{B2}=50.3\ {\rm {m}};\\{\mbox{vertical depth}}&=33+50.3/\cos 13.6^{\circ }=84\ {\rm {m}};\\\tan \xi &=(84-18)/220=16.7^{\circ }.\end{aligned}}}

Note that this value of ${\displaystyle \xi }$ is based on intercept times plus velocities, whereas the value ${\displaystyle 13.6^{\circ }}$ is based on the velocities only.