Fast ray tracing method in 3-D structure and its proof of positive definiteness

Abstract

Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symmetry block tridiagonal matrix equation set was deduced. Further, the positive definiteness of coefficient matrix was discussed, and the positive definiteness was accurately proved in a mathematical way. It assured that the algorithm was well-posed. Associated with iterative method, the solution to ray tracing can be got through step-by-step linearized iteration of the nonlinear problem. An algorithm of the whole path iterative ray tracing method in three-dimensional velocity structure was obtained. This method shows a clear and simple as well as explicit computation formula, which makes ray tracing computation easily applicable in practice. The correction vector is obtained through finding the solution to the positive definite block tridiagonal equation set, which ensures the method is robust convergence. This study offers a new kind of feasible and efficient ray tracing method for three dimensional seismic migration and tomography. Meanwhile, it also provides the prerequisite guarantee to design a fast algorithm.

Keywords

Fermat’s principle / ray tracing / ray path / positivity / seismic migration / tomography

References

[1] KaralF. C., KellerJ. B.. Elastic wave propagation in homogeneous and inhomogeneous media[J]. J Accoust Soc Am, 1959, 31(6): 694-705

[2] ČervenýV., MolotkovI. A., PšenčíkI.Ray Method in Seismology[M], 1977, Prague, Charles University Press

[3] FarraV.. Bending method revisited: a Hamiltonian approach[J]. Geophys J Int, 1992, 109(1): 138-150

[4] GaoE.-g., XuG.-m., ZhaoYi.. Segmentally-iterative ray tracing method for any interface[J]. Oil Geophysical Prospecting, 1998, 33(1): 54-60

[5] GaoE.-g., XuG.-ming.. A new kind of step by step iterative ray-tracing method[J]. Acta Geophysica Sinica, 1996, 39(S): 302-308

[6] GaoE.-g., XuG.-m., LiG.-p., et al.. A new total iterative ray-tracing method in random interface[J]. Acta Custica, 2002, 27(3): 282-287

[7] SambridgeM. S., KennettB. L. N.. Boundary value ray tracing in heterogeneous medium: a simple and versatile algorithm[J]. Geophys J Int, 1990, 101(1): 157-168

[8] UmJ., ThurberC. H.. A fast algorithm for two-point seismic ray tracing[J]. Bull Seis Soc Am, 1987, 77(3): 972-986

[9] XuG.-m., WeiS., GaoE.-g., et al.. Block model-building and ray-tracing in 2-D complicated medium[J]. Oil Geophysical Prospecting, 2001, 36(2): 213-219

[10] ChiuS. K. L., KanasewichE. R., PhadkeS.. Three-dimensional determination of structure and velocity by seismic tomography[J]. Geophysics, 1986, 51(8): 1559-1571

[11] FarraV.. Ray tracing in complex media[J]. J Appl Geophys, 1993, 30(1): 55-73

[12] LiuH., MengF.-l., LiY.-ming.. The interface grid method for seeking global minimum travel-time and the correspondent ray path[J]. Acta Geophysica Sinica, 1995, 38(6): 823-832

[13] VidaleJ. E.. Finite different calculation of travel-time in three dimensions[J]. Geophysics, 1990, 55(5): 521-526

[14] MaZ.-m., LiY.-da.. Two-step ray tracing method[J]. Acta Geophysica Sinica, 1991, 34(4): 501-508

[15] YangC.-c., LengC.-b., LiY.-ming.. Fast and accurate ray tracing in 3-D media[J]. Acta Geophysica Sinica, 1997, 40(3): 414-420

[16] YangW.-c., LiY.-m., LiS.-x., et al.Applied seismic Tomography[M], 1993, Beijing, Geological Press

[17] ZhangL.-b., YaoZ.-x., JiChen.. Finite difference calculation of seismic first-arrival travel-time[J]. Progress in Geophysics, 1996, 11(4): 47-52

[18] ZhouZ.-s., ZhangS.-m., ChenL.-jun.. Seismic ray tracing calculation based on parabolic travel-time interpolation[J]. J Cent South Univ Technol, 2004, 11(2): 199-205