next up previous print clean
Next: FUTURE WORK Up: Rickett, et al.: STANFORD Previous: Synthetic Example

TOMOGRAPHY IN TIME VS. TOMOGRAPHY IN DEPTH

Because of the velocity-depth coupling, non-linear tomographic velocity estimation in the depth coordinates often does not converge effectively. To avoid some of the drawbacks of traditional tomography Biondi et al. proposed a tomographic velocity estimation in the two-ways traveltime coordinates ($\tau$,$\xi$). The method is based on a transformation of the eikonal equation from the depth coordinates (z,x) into the time coordinates ($\tau$,$\xi$), according the following transformation ov variables:
\begin{eqnarray}
\tau(z,x) &=& \int_0^z {2\over V(z',x)} dz' \nonumber \\ \xi(z,x) &=& x .\end{eqnarray}
(70)

 
panel
panel
Figure 8
Velocity model after one, two, three, and four non-linear iterations.
view burn build edit restore

Starting from the eikonal equation in the time coordinates, we can remap the tomography goals of equation ([*]) into the time coordinates. The operators $\bold{L_s^i}$, $\bold{H}$, $\bold{G}$ are different in the time domain from the equivalent depth-domain operators; and also $\bold A^{-1}$ has a slightly different orientation, because it now operates on a velocity function defined in ($\tau$,$\xi$).

Some of the advantages of formulating the tomography problem in ($\tau$,$\xi$) domain can be seen in Figure [*]. In this new domain reflector movement is minimal, significantly decreasing the velocity-depth connectivity problem that we saw in (z,x) space. The initial reflector position is closer to the correct reflector position. As a result, the orientation of the steering filters more accurately follows true reflector dip.

 
movement
movement
Figure 9
Left, depth velocity model with the correct reflector position (dashed) and the estimated reflector position using the initial guess at the velocity model (solid). Right, ($\tau$,$\xi$) velocity model with the same reflectors superimposed. Note how reflector movement is significantly less than in the ($\tau$,$\xi$) case.
view burn build edit restore

To analyze the advantages of working in ($\tau$,$\xi$) space we performed two non-linear iterations of tomography and remapped the resulting velocity model back into (z,x) space (Figure [*]). Comparing the results of tomography in the time domain (Figure [*]), with the results of tomography in the depth domain (Figure [*]), we can notice how the doublet, seen in the depth tomography case, has significantly decreased and how the dome shape is much more prominent.


next up previous print clean
Next: FUTURE WORK Up: Rickett, et al.: STANFORD Previous: Synthetic Example
Stanford Exploration Project
7/5/1998