A robust approach to time-to-depth conversion and interval velocity estimation from time migration in the presence of lateral velocity variations |

vz0
(Top) a synthetic model and
(middle) Dix velocity converted to depth. Both overlaid with image
rays. (Bottom) the model perturbation for testing linearization.
Figure 12. |
---|

Figure 12 shows a synthetic model borrowed from Cameron et al. (2008). The Dix inversion recovers the shallow part of the model but deteriorates quickly as geometrical spreading of image rays grows in the deeper section.

As a simple verification for the linearization process, we add a small positive velocity perturbation at location km to the synthetic model. Comparisons between the exact and linearly predicted attributes are illustrated in Figures 13, 14 and 15. In accordance with forward modeling, where we solve firstly , then , and finally , the linearization (see Appendix B) is carried out following the same sequence. First, Figure 13 justifies our upwind finite-differences implementation of the linearized eikonal equation. The positive perturbation in in Figure 12 causes to decrease in a narrow downwind region. Next, the area affected by the perturbation in Figure 14 is wider than that in Figure 13. It also has both positive and negative amplitudes. These phenomenon are physical because image rays should bend in opposite directions in response to the perturbation. Finally, effects in cost in Figure 15 show alternating polarities and are broader in width compared to that of and . They indicate a complicated dependency of on . Note the good agreements in both shape and magnitude between exact and linearly predicted quantities in all three steps.

Because there is no analytical formula for Dix velocity in this model, we compute by tracing image rays numerically in the exact model . Also, based on Figure 12, there is no in-flow boundary other than . Therefore, we do not need to extend the domain as in the preceding examples. We use the Dix-inverted model as the prior model and run the inversion. It turns out that the first linearization update is sufficient for achieving the desired global minimum as shown in Figure 16.

pdt
(Top) exact and (bottom)
linearly predicted by equation B-4.
Figure 13. |
---|

pdx
(Top) exact and (bottom)
linearly predicted by equation B-3.
Figure 14. |
---|

diffcost
(Top) exact and (bottom)
linearly predicted by equation B-5.
Figure 15. |
---|

bgrad
(Top) the exact
and (middle) the computed of the first linearization step.
(Bottom) the inverted interval velocity model. Compare with
Figure 12.
Figure 16. |
---|

A robust approach to time-to-depth conversion and interval velocity estimation from time migration in the presence of lateral velocity variations |

2015-03-25