Multiple ray tracing in 3D layered media using the traveltime gradient method
Zhang Tingting1,2, Zhang Dong1, Qiu Da2
1. School of Physics and Technology, Wuhan University, Wuhan, Hubei 430072, China;
2. School of Information Engineering, Hubei University for Nationalities, Enshi, Hubei 445000, China
Abstract：The purpose of this paper is to present a multistage maximum taveltime gradient ray tracing method (MTG) for tracing multiple seismic arrivals in a 3D layer media. First a reflected or refracted wavefront can be calculated in the order of the multiple transmissions by partitioning layers into separate computational domains in which the dynamic successive partitioning scheme is applied. The next step is to track the multiple arrivals along the negative traveltime gradient using the traveltime information we have obtained. The successive partitioning of cell interfaces is used to improve the efficiency of the multiple wavefront traveltime field calculation,hence,we can find a more accurate secondary source point with a minimum traveltime. The multiply-arriving ray path can be traced from the receiver to the source along the negative traveltime gradient. In the area of velocity interface we adopted a local ray tracing strategy to improve the precision of multiply arrival ray paths. Several examples (including 3D Marmousi model and layered media) and error analysis indicate that the multistage MTG method retains the basic characteristics of the first arrival MTG method with high computational accuracy, algorithm simplicity, and numerical stability,and the result demonstrate the viability of the multistage MTG in heterogeneous media.
张婷婷, 张东, 邱达. 三维层状介质中基于走时梯度的多次波射线追踪[J]. 石油地球物理勘探, 2014, 49(6): 1097-1105.
Zhang Tingting, Zhang Dong, Qiu Da. Multiple ray tracing in 3D layered media using the traveltime gradient method. OGP, 2014, 49(6): 1097-1105.
Huang Guojiao, Bai Chaoying. Simultaneous inversion with multiple traveltimes within 2-D complex layered media.Chinese Journal of Geophysics, 2010,53(12):2972-2981.
[2]
James W D, Satish C and Timothy A. Three-dimensional tomographic inversion of combined reflection and refraction seismic traveltime data. Geophysical Journal International. 2003,152(1):79-93.
[3]
Rawlinsona N, Pozgaya S, Fishwickb S. Seismic tomography: A window into deep Earth. Physics of the Earth and Planetary Interiors, 2010,178(3):101-135.
Li Peng, Liu Yike, Chang Xu et al. Progress on the multiple problems. Progress in Geophysics, 2006,21(3):888-897.
[5]
Vidale J E. Finite-difference calculations of travel-times.Bulletin of the Seismological of America, 1988,78(6):2062-2076.
[6]
Vidale J E. Finite-difference calculations of travel-times in three dimensions. Geophysics, 1990,55(5):521-526.
[7]
Hole J A and Zelt B C. 3-D finite-difference reflection travel times. Geophysical Journal International, 1995,121(2):427-434.
[8]
Afnimar and Koketsu K. Finite difference traveltime calculation for head waves travelling along an irregular interface.Geophysical Journal International, 2000,143(3):729-734.
[9]
Podvin P and Lecomte I. Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools. Geophysical Journal International, 1991,105(1):271-284.
[10]
Qin F, Luo Y, Olsen K B, Cai W et al. Finite-difference solution of the eikonal equation along expanding wavefronts. Geophysics, 1992,57(3):478-487.
[11]
Nakanishi I and Yamaguchi K. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional islandarc structure. Physics of the Earth and Planetary Interiors, 1986,34(2):195-201.
[12]
Moser T J. Shortest path calculation of seismic rays. Geophysics, 1991,56(1):59-67.
[13]
Fischer R and Lees J M. Shortest path ray tracing with sparse graphs.Geophysics,1993,58(7):987-996.
[14]
Cheng N and House L. Minimum traveltime calculations in 3-D graph theory. Geophysics, 1996,61(6):1895-1898.
[15]
Asakawa E and Kanawa T. Seismic ray tracing using linear traveltime interpolation. Geophysical Prospecting, 1993, 41(1): 99-111.
Nie Jianxin, Yang Huizhu. Quadratic/linear travel time interpolation of seismic ray-tracing. Journal of Tsinghua University (Science and Technology Edition). 2003,43(11):1495-1498.
Zhang Dong, Xie Baolian,Yang Yan et al. A ray tracing method based on improved linear traveltime interpolation. Chinese Journal of Geophysics, 2009,52(1):200-205.
[18]
Cardarelli E, Cerreto A. Ray tracing in elliptical anisotropic media using the linear traveltime inter polation (LTI) method applied to traveltime seismic tomography. Geophysical Prospecting,2002,50(1):55-72.
Sun Zhangqing, Sun Jianguo, Han Fuxing. The comparion of three schemes for computing seismic wave traveltimes in complex topographical conditions. Chinese Journal of Geophysics, 2012, 55(2):560-568.
[20]
Sethian J A. A fast marching level set method for monotonically advancing fronts. Proceedings of the National Academy of Sciences, 1996,93(4):1591-1595.
[21]
Sethian J A and Popovici A M. 3-D traveltime computation using the fast marching method.Geophysics, 1999,64(2):516-523.
[22]
Alkhalifah T and Fomel S. Implementing the fast marching eikonal solver: spherical versus Cartesian coordinates. Geophysical Prospecting, 2001,49(2):165-178.
[23]
Popovici A M and Sethian J A. 3-D imaging using higher order fast marching traveltimes. Geophysics, 2002,67(2):604-609.
[24]
Moser T J. Shortest path calculation of seismic rays. Geophysics, 1991,56(1):59-67.
[25]
Podvin P, Lecomte I. Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools. Geophysical Journal International, 1991,105(1):271-284.
[26]
Hole J A and Zelt B C. 3-D finite-difference reflection traveltimes. Geophysical Journal International, 1995,121(2):427-434.
[27]
Riahi M A, Juhlin C. 3-D interpretation of reflected arrival times by finite-difference techniques. Geophysics, 1994,59(5):844-849.
Wang Tongkui, Zhang Meigen, Li Xiaofan et al. Interface points as secondary sources method raytracing PS converted waves. Progress in Geophysics, 2007,22(1):165-170.
Zhao Gaishan,Hao Shouling,Yang Erhao et al. Seis-mic ray tracing algorithm based on the linear traveltime interpolation. GPP,1998,37(2):14-24.
[31]
Rawlinson N and Sambridge M. Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 2004,69(5):1338-1350.
[32]
De Kool M, Rawlinson N, Sambridge M. A praetieal gridbased method for tracking mutiple refraction and reflection Phases in three-dimensional heterogeneous media. Geophysical Journal International, 2006,167(1):253-270.
Tang Xiaoping, Bai Chaoying. Multiple ray tracing in 3-D layered media with the shortest path method. Chinese Journal of Geophysics, 2009, 52(10):2635-2643.
[35]
Huang Yueqin, Zhang Jianzhong and Liu Qinghuo.Three-dimensional GPR ray tracing based on wavefront expansion with irregular cells. IEEE Transactions on Geoscience and Remote Sensing, 2011,49(2):679-687.
[36]
Zhang Dong, Zhang Tingting et al. A new 3-D ray tracing method based on LTI using successive partitioning of cell interfaces and traveltime gradients. Journal of Applied Geophysics. 2013,92:20-29.
Zhang Dong, Fu Xiangru, Yang Yan et al. 3-D seismic ray tracing algorithm based on LTI and partition of grid interface. Chinese Journal of Geophysics, 2009,52(9):2370-2376.
