next up previous print clean
Next: Nonhyperbolic curves Up: CURVED WAVEFRONTS Previous: Root-mean-square velocity

Layered media

From the assumption that experimental data can be fit to hyperbolas (each with a different velocity and each with a different apex $ \tau $)let us next see how we can fit an earth model of layers, each with a constant velocity. Consider the horizontal reflector overlain by a stratified interval velocity v(z) shown in Figure 10.

 
stratrms
Figure 10
Raypath diagram for normal moveout in a stratified earth.

stratrms
view burn build edit restore

The separation between the source and geophone, also called the offset, is 2h and the total travel time is t. Travel times are not be precisely hyperbolic, but it is common practice to find the best fitting hyperbolas, thus finding the function $V^2(\tau)$.

 
 \begin{displaymath}
t^2 \eq \tau^2 + \frac{4h^2}{V^2(\tau)}\end{displaymath} (24)
where $ \tau $ is the zero-offset two-way traveltime.

An example of using equation (24) to stretch t into $ \tau $is shown in Figure 11. (The programs that find the required $V(\tau )$ and do the stretching are coming up in chapter [*].)

 
nmogath
nmogath
Figure 11
If you are lucky and get a good velocity, when you do NMO, everything turns out flat. Shown with and without mute.


view burn build edit restore

Equation (21) shows that $V(\tau )$is the ``root-mean-square'' or ``RMS'' velocity defined by an average of v2 over the layers. Expressing it for a small number of layers we get  
 \begin{displaymath}
V^2(\tau) \eq \frac{1}{\tau}\ \sum_i v^2_i \Delta\tau_i\end{displaymath} (25)
where the zero-offset traveltime $ \tau $ is a sum over the layers:  
 \begin{displaymath}
\tau \eq \sum_i \ \Delta\tau_i\end{displaymath} (26)
The two-way vertical travel time $\tau_i$in the ith layer is related to the thickness $\Delta z_i$ and the velocity vi by  
 \begin{displaymath}
\Delta\tau_i \eq \frac{2\ \Delta z_i}{v_i} \ \ .\end{displaymath} (27)

Next we examine an important practical calculation, getting interval velocities from measured RMS velocities: Define in layer i, the interval velocity vi and the two-way vertical travel time $\Delta\tau_i$.Define the RMS velocity of a reflection from the bottom of the i-th layer to be Vi. Equation (25) tells us that for reflections from the bottom of the first, second, and third layers we have
      \begin{eqnarray}
V_1^2 &=& {v_1^2\Delta\tau_1 
 \over \Delta\tau_1 }
\\ 
V_2^2 &...
 ...3^2\Delta\tau_3 
 \over \Delta\tau_1 + \Delta\tau_2 +\Delta\tau_3}\end{eqnarray} (28)
(29)
(30)

Normally it is easy to measure the times of the three hyperbola tops, $\Delta\tau_1$, $\Delta\tau_1 + \Delta\tau_2$and $\Delta\tau_1 + \Delta\tau_2 +\Delta\tau_3$.Using methods in chapter [*] we can measure the RMS velocities V2 and V3. With these we can solve for the interval velocity v3 in the third layer. Rearrange (30) and (29) to get

      \begin{eqnarray}
 (\Delta\tau_1 + \Delta\tau_2 +\Delta\tau_3)
V_3^2 &=& v_1^2\De...
 ..._1 + \Delta\tau_2)
V_2^2 &=& v_1^2\Delta\tau_1+ v_2^2\Delta\tau_2 \end{eqnarray} (31)
(32)
and subtract getting the squared interval velocity v32  
 \begin{displaymath}
v_3^2 \eq {
 (\Delta\tau_1 + \Delta\tau_2 +\Delta\tau_3) V_3^2 -
 (\Delta\tau_1 + \Delta\tau_2 ) V_2^2
 \over
 \Delta\tau_3}\end{displaymath} (33)
For any real earth model we would not like an imaginary velocity which is what could happen if the squared velocity in (33) happened to be negative. You see that this means that the RMS velocity we estimate for the third layer cannot be too much smaller than the one we estimate for the second layer.


next up previous print clean
Next: Nonhyperbolic curves Up: CURVED WAVEFRONTS Previous: Root-mean-square velocity
Stanford Exploration Project
12/26/2000