Multiples raytracing in 3D anisotropic TI media with a tetrahedron cell model
Li Xingwang1, Bai Chaoying1,2, Li Xiaoling3
1. College of Geology Engineering and Geomatics, Chang'an University, Xi'an, Shannxi 710054, China;
2. Institute of Computing Geophysics, Chang'an University, Xi'an, Shannxi 710054, China;
3. Patent Examination Cooperation Center of the Patent Office, SIPO, Zhengzhou, Henan 450000, China
Abstract:The current algorithms can only trace first arrivals and primary reflections (or conversions),but cannot well trace multiples.To overcome this problem,with a tetrahedron cell model we extend the function of multistage irregular shortest-path method into 3D anisotropic TTI media for tracking multiple arrivals composed of any kind of combinations of transmissions,conversions and reflections.The basic principle is that the 7 parameters (5 elastic parameters and 2 titled angles of symmetry axis) of TTI media sampled at the primary nodes are used to calculate the group velocity of three-mode (qP,qSV,and qSH) waves.The group velocity values at the secondary nodes are obtained by a linear interpolation function via the defined values at the primary nodes of the cell.Finally we conduct grid (or cell) wavefront expansion to trace multiple arrivals from one layer to the next.The calculation in the uniform anisotropic media indicates that the numerical result is coincided with the analytic solution,which verifies the accuracy of the methodology.Several simulation results show that the proposed method is able to efficiently approximate the undulated topography and velocity discontinuity,and suitable for any combination of multiple arrivals tracking in TTI media with any titled angle of symmetry axis and not subject to the strength (no weak anisotropic assumption) of anisotropic media.
李兴旺, 白超英, 李晓玲. 四面体单元剖分下三维各向异性TI介质中多次波射线追踪[J]. 石油地球物理勘探, 2017, 52(1): 48-55.
Li Xingwang, Bai Chaoying, Li Xiaoling. Multiples raytracing in 3D anisotropic TI media with a tetrahedron cell model. OGP, 2017, 52(1): 48-55.
erveny V.Seismic rays and rays intensities in inhomogeneous anisotropic media.Geophysical Journal of Royal Astronomic Society,1972,29(1):1-13.
[2]
Babich V M.Ray method of calculating the intensity of wavefronts in the case of a heterogeneous,anisotropic,elastic medium.Geophysical Journal International,1994,118(2):379-383.
[3]
Cardarelli E,Cerreto A.Ray tracing in elliptical anisotropic media using the linear traveltime interpolation method applied to traveltime seismic tomography.Geophysical Prospecting,2002,50(1):55-72.
[4]
Qian J,Symes W W.Finite-difference quasi-P traveltimes for anisotropic media.Geophysics,2002,67(1):147-155.
Ma Detang,Zhu Guangming.Calculating traveltime of the first arrival in transversely isotropic medium.OGP,2006,41(1):26-31.
[6]
Grechka V Y,McMechan G A.3-D tow-point ray tracing for heterogeneous,weakly transversely isotropic media.Geophysics,1996,61(6):1883-1894.
[7]
Psencik I,Farra V.First-order ray tracing for qP waves in inhomogeneous weakly anisotropic Media.Geophysics,2005,70(6):D65-D75.
[8]
erveny V,Firbas P.Numerical modelling and inversion of traveltimes of seismic body waves in inhomogeneous anisotropic media.Geophysical Journal Royal Astronomical Society,1984,76(1):41-51.
[9]
Gajewski D,Psencik I.Computation of high-frequency seismic wavefields in 3-D laterally inhomogeneous anisotropic media.Geophysical Journal Royal Astronomical Society,1987,91(2):383-411.
[10]
Shearer P M,Chapman C H.Ray tracing in anisotropic media with a linear gradient.Geophysical Journal International,1988,94(3):575-580.
[11]
Farm V.First-order ray tracing for qS waves inhomogeneous weakly anisotropic media.Geophysical Journal International,2005,161(2):309-324.
Zhang Zhongjie and He Qiaoden.The study of refracted wave and its traveltime curve in anisotropic medium.Journal of Xi'an Petroleum Institute,1992,7(3):1-6.
Zhao Aihua,Zhang Meigen,Ding Zhifeng.Seismic traveltime computation for transversely isotropic media.Chinese Journal of Geophysics,2006,49(6):1762-1769.
[14]
Zhou B,Greenhalgh S.On the computation of elastic wave group velocities for a general anisotropic medium.Journal of Geophysics and Engineering,2004,1(3):205-215.
[15]
Zhou B,Greenhalgh S.Shortest path'ray tracing for most general 2D/3D anisotropic media.Journal of Geophysics and Engineering,2005,2(1):54-63.
[16]
Zhou B,Greenhalgh S.Raypath and traveltime computations for 2D transversely isotropic media with dipping symmetry axes.Exploration Geophysics,2006,37(2):150-159.
Ma Detang,Zhu Guangming,Fan Ting'en.Search algorithm of first arrival traveltime in 2D TTI medium.OGP,2011,46(5):710-714.
[18]
Bai C Y,Huang G J,Li X L et al.Ray tracing of multiply transmitted/reflected/converted waves in 2D/3D layered anisotropic TTI media and application to crosswell traveltime tomography.Geophysical Journal International,2013,195(2):1068-1087.
Li Xiaoling,Bai Chaoying,Hu Guangyi.Multiple ray tracing in an undulated layered TI media.OGP,2013,48(6):924-931.
[20]
Bai C Y,Li X L,Tang X P.Seismic wavefront evolution of multiply reflected,transmitted,and converted phases in 2D/3D triangular cell model.Journal of Seismology,2011,15(4):637-652.
[21]
Bai C Y,Li X,Wang Q L,Peng J B.Multiple arrival tracking within irregular triangular or tetrahedral cell model.Journal of Geophysics and Engineering,2012,9(1):29-38.
[22]
Daley P,Hron F.Reflection and transmission coeffi-cients for transversely isotropic media.Bulletin of the Seismological Society of America,1977,67(3):661-675.
[23]
Berryman J G.Long-wave elastic anisotropy in transversely isotropic media.Geophysics,1979,44(5):896-917.
[24]
Dwyer R A.Higher-dimensional voronoi diagram in linear expected time.Discrete and Computational Geometry,1991,6(4):64-67.
[25]
Radovitzky R,Ortiz M.Tetrahedral mesh generation based on node insertion in crystal lattice arrangements and advancing-front-Delaunay triangulation.Computer Methods in Applied Mechanics and Engineering,2000,187(3-4):543-569.
Guo Jiyuan,Gong Junfang.Algorithms of producing tetrahedral network from three dimensional dispersed data.Journal of China University of Geosciences (Earth Science),2002,27(3):271-273.
Kong Juanhua,Zheng Jiangbin.An algorithm of generating unstructured tetrahedrons from 3D discrete points.Computer Engineering and Science,2009,31(1):35-37.