Fast algorithm and numerical simulation for ray-tracing in 3D structure

Abstract

Beginning with the method of whole path iterative ray-tracing and according to the positive definiteness of the coefficient matrix of the systems of linear equations, a symmetry block tridiagonal matrix was decomposed into the product of block bidiagonal triangular matrix and its transpose by means of Cholesky decomposition. Then an algorithm for solving systems of block bidiagonal triangular linear equations was given, which is not necessary to treat with the zero elements out of banded systems. A fast algorithm for solving the systems of symmetry block tridiagonal linear equations was deduced, which can quicken the speed of ray-tracing. Finally, the simulation based on this algorithm for ray-tracing in three dimensional media was carried out. Meanwhile, the segmentally-iterative ray-tracing method and banded method for solving the systems of block tridiagonal linear equations were compared in the same model mentioned above. The convergence condition was assumed that the L-2 norm summation for m_{k, 1} and m_{k, 2} in the whole ray path was limited in 10⁻⁶. And the calculating speeds of these methods were compared. The results show that the calculating speed of this algorithm is faster than that of conventional method and the calculated results are accurate enough. In addition, its precision can be controlled according to the requirement of ray-tracing

Keywords

ray-tracing / seismic migration / Cholesky decomposition / tomography

References

[1] AkiK., ChristofferssonA., HusebyeE. S.. Determination of the three-dimensional seismic structure of the lithosphere [J]. J Geophys Res, 1977, 82(2): 277-296

[2] DuveneckE.. 3D tomographic velocity model estimation with kinematic wavefield attributes [J]. Geophysical Prospecting, 2004, 52(6): 535-545

[3] IvanssonS.. Seismic borehole tomography theory and computational methods [J]. Proc IEEE, 1986, 74(2): 328-338

[4] SchusterG. T.. Fermat’s interferometric principle for target-oriented traveltime tomography [J]. Geophysics, 2005, 70(4): U47-50

[5] YangZ.-x., ZhangX.-k., YangJ., ZhaoJ.-ren.. Joint inversion of 3-D crustal structure and velocity distribution [J]. Geophysical Progress, 1997, 12(1): 41-52

[6] ZhangZ.-j., QinY.-l., ChenY., ZhangC.-y., SunS.-x., ZhaoB., LiuY.-f., LiuZ.-kuan.. Reconstruction of the semblance section for the crust and mantle reflection structure by wide-angle reflection seismic data [J]. Chinese J Geophys, 2004, 47(3): 469-474

[7] ZhangZ.-j., TengJ.-w., ZhangL.-b., YangD.-hui.. The modified Snell law for seismic ray in anisotropic media [J]. Chinese Science Bulletin, 1995, 40(20): 1724-1728

[8] XuT., XuG.-m., GaoE.-g., ZhuL.-b., JiangX.-yi.. Block modeling and shooting ray in complex 3-D media [J]. Chinese J Geophys, 2004, 47(3): 1118-1126

[9] XuT., XuG.-m., GaoE.-g., LiY.-c., JiangX.-y., LuoK.-yun.. Block modeling and segmentally iterative ray tracing in complex 3D media [J]. Geophysics, 2006, 71(3): T41-51

[10] XuT., XuG.-m., GaoE.-g., JiangX.-yi.. Refraction ray-tracing in complex medium [J]. Oil Geophysical Prospecting, 2004, 39(6): 690-693

[11] XuT., XuG.-M., GaoE.-G., JiangX.-y., LuoK.-yun.. 3-D shooting ray-tracing sub-triangle method [J]. Oil Geophysical Prospecting, 2005, 40(4): 391-399

[12] GjystdalH., IversenE., LecomteI., KaschwichT., DrottningA., MispelJ.. Improved applicability of ray tracing in seismic acquisition, imaging, and interpretation [J]. Geophysics, 2007, 72(5): SM261-271

[13] BuskeS., KästnerU.. Efficient and accurate computation of seismic traveltimes and amplitudes [J]. Geophysical Prospecting, 2004, 52(4): 313-322

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

[15] 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

[16] GaoE.-g., XuG.-m., JiangX.-y., LuoK.-y., LiuT.-q., XieD., ShiJ.-liang.. Iterative ray-tracing method segment by segment under 3-D construction [J]. Oil Geophysical Prospecting, 2002, 37(1): 11-16

[17] GaoE.-g., UkH., TenJ.-wen.. Fast ray-tracing method in 3-D structure and its proof of positive definiteness [J]. J Cent South Univ Technol, 2007, 14(1): 100-103

[18] GolubG. H., van LoanC. F.Matrix computation [M], 1983, Baltimore, Johns Hopkins University Press