Journal of Jilin University(Earth Science Edition) ›› 2019, Vol. 49 ›› Issue (4): 1169-1181.doi: 10.13278/j.cnki.jjuese.20190076

Previous Articles     Next Articles

Ray Tracing and Kinematic Characteristics of Different Types of Seismic Waves in Complex Seabed

Sun Zhangqing, Wang Dengke, Han Fuxing   

  1. College of GeoExploration Science and Technology, Jilin University, Changchun 130026, China
  • Received:2019-04-08 Online:2019-07-26 Published:2019-07-26
  • Supported by:
    Supported by National Science Foundation of China (41404085), Training Program for Excellent Young Teachers in Jilin University (419080500337) and Project of Jilin Oilfield Group Company (JS2018-W-32-JZ-10-20)

Abstract: In order to establish a flexible, stable and accurate ray tracing method for complex seabed geological conditions such as rugged seabed, steep slope seabed, submarine volcano, and to analyze the kinematic characteristics of different types of seismic waves, we put forward a fast marching upwind interpolation method through taking the advantages of many kinds of algorithms together. Firstly, a complex seabed model is built by using hybrid grids. Secondly, through incorporating the idea of upwind differences, the linear interpolation scheme is adapted to construct high precision,unconditionally stable, and flexible local computational formulas for travel-time and ray-path. Then, by using these formulas and the flexibility of the fast marching method, the travel-time information of various types of seismic waves in the whole model are calculated, and by using reverse tracing, the ray-path is calculated. Finally, the accuracy of this method is analyzed. It is found that the computational accuracy is relatively high by using this method for travel-time and ray-path, and the computational accuracy of reflected waves is much higher than that of incident waves. Furthermore, based on some computational examples, the kinematic characteristics of seismic waves can be obtained, for example, refracted waves are abundant in the first-arrival wave, reflected wave caused by steep structures can be received at large offsets, and the ray paths of various waves are complicated in rugged seabed, etc. Finally, based on the above kinematic characteristics, we suggest to use the travel-time of direct wave for replacing the one of first-arrival wave and increase the length of acquisition array to improve imaging quality of steep dip structures.

Key words: complex seabed conditions, different types of seismic waves, travel-time computation, ray tracing, kinematics characteristic analysis

CLC Number: 

  • P631.4
[1] Valente L S S, Santos H B, Costa J C, et al. Time-to-Depth Conversion and Velocity Estimation by Image-Wavefront Propagation[J]. Geophysics, 2017, 82(6):U75-U85.
[2] Huang Xingguo, Sun Hui, Sun Jianguo. Born Modeling for Heterogeneous Media Using the Gaussian Beam Summation Based Green's Function[J]. Journal of Applied Geophysics, 2016, 131:191-201.
[3] Julian B R, Gubbins D. Three-Dimensional Seismic Ray Tracing[J]. Geophysics, 1977, 43:95-113.
[4] Moser T J. Shortest Path Calculation of Seismic Rays[J]. Geophysics, 1991, 56(1):59-67.
[5] Vidale J. Finite-Difference Calculation of Travel Times[J]. Bull Seism Soc Am, 1988, 78(6):2062-2076.
[6] Hao Q, Alkhalifah T. An Acoustic Eikonal Equation for Attenuating Transversely Isotropic Media with a Vertical Symmetry Axis Acoustic Attenuating VTI Media[J].Geophysics, 2017, 82(1):C9-C20.
[7] Waheed U B, Alkhalifah T. A Fast Sweeping Algorithm for Accurate Solution of the Tilted Transversely Isotropic Eikonal Equation Using Factorization[J]. Geophysics, 2017, 82(6):WB1-WB8.
[8] Hu J, Cao J, Wang H, et al. 3D Travel-Time Computation for Quasi-P-Wave in Orthorhombic Media Using Dynamic Programming[J]. Geophysics, 2017, 83(1):C27-C35.
[9] Asakawa E, Kawanaka T. Seismic Ray Tracing Using Linear Traveltime Interpolation[J]. Geophysical Prospecting, 1993, 41(1):99-111.
[10] Schneider W A J, Ranzinger K A, Balch A H, et al. A Dynamic Programming Approach to First Arrival Traveltime Computation in Media with Arbitrarily Distributed Velocities[J]. Geophysics, 1992, 57(1):39-50.
[11] Han S, Zhang W, Zhang J. Calculating QP-Wave Traveltimes in 2-D TTI Media by High-Order Fast Sweeping Methods with a Numerical Quartic Equation Solver[J]. Geophysical Journal International, 2017, 210(3):1560-1569.
[12] Vinje V, Iversen E, Gjoystdal H. Traveltime and Amplitude Estimation Using Wavefront Construction[J]. Geophysics, 1993, 58(8):1157-1166.
[13] Leidenfrost A, Ettrich N, Gajewski D, et al. Comparison of Six Different Methods for Calculating Traveltimes[J]. Geophysical Prospecting, 1999, 47:269-297.
[14] 徐涛,徐果明,高尔根,等.三维复杂介质的块状建模和试射射线追踪[J].地球物理学报,2004,47(6):1118-1126. Xu Tao, Xu Guoming, Gao Ergen, et al. Block Modeling and Shooting Ray Tracing in Complex 3-D Media[J]. Chinese Journal of Geophysics,2004, 47(6):1118-1126.
[15] Sun Jianguo, Sun Zhangqing, Han Fuxing. A Finite Difference Scheme for Solving the Eikonal Equation Including Surface Topography[J]. Geophysics, 2011, 76(4):T53-T63.
[16] Zhou B, Greenhalgh S. Shortest Path' Ray Tracing for Most General 2D/3D Anisotropic Media[J]. Geophysics, 2005, 2:54-63.
[17] Sethian J A, Popovici A M. 3-D Traveltime Computation Using the Fast Marching Method[J]. Geophysics, 1999, 64:516-523.
[18] 熊金良,李国发,王濮,等.地震波走时计算的逆风差分算法[J].物探化探计算技术,2004, 26(4):295-298. Xiong Jinliang, Li Guofa, Wang Pu, et al. Upwind Finite-Difference Solution of Eikonal Equation[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2004, 26(4):295-298.
[19] Rawlinson N, Sambridge M. Multiple Reflection and Transmission Phases in Complex Layered Media Using a Multistage Fast Marching Method[J]. Geophysics, 2004, 69(5):1338-1350.
[20] Li Zhenchun, Liu Yulian, Zhang Jianlei, et al. Finite-Difference Calculation of Traveltimes Based on Rectangular Grid[J]. Acta Seismologica Sinica, 2004, 26(6):644-650.
[21] Liu Zhenwu, Sun Yonghe, Sun Jianguo, et al. A Method for Computing the 2D Seismic Traveltime Including the Surface Topography by Combining the Linear Interpolation and the Narrow Band Technique[C]//International Geophysical Conference & Exposition. Beijing:CPS/SEG, 2009:130.
[22] Wang Huazhong, Kuang Bin, Ma Zaitian, et al. 3-D Traveltime Calculation in Arbitrary Velocity Distribution with Dynamic Programming Approach[J]. Chinese Journal of Geophysics, 2001, 44(Sup. 1):179-189.
[23] 孙章庆,孙建国,韩复兴.复杂地表条件下基于线性插值和窄带技术的地震波走时计算[J].地球物理学报,2009,52(11):2846-2853. Sun Zhangqing, Sun Jianguo, Han Fuxing. Traveltimes Computation Using Linear Interpolation and Narrow Band Technique Under Complex Topographical Conditions[J]. Chinese Journal of Geophysics, 2009, 52(11):2846-2853.
[24] 孙章庆.起伏地表条件下的地震波走时与射线路径计算[D].长春:吉林大学, 2011. Sun Zhangqing. The Seismic Traveltimes and Raypath Computation Under Undulating Earth's Surface Condition[D]. Changchun:Jilin University, 2011.
[25] 高亮,李幼铭,陈旭荣,等.地震射线辛几何算法初探[J]. 地球物理学报, 2000, 43(3):402-410. Gao Liang, Li Youming, Chen Xurong, et al. An Attempt to Seismic Ray Tracing with Symplectic Algorithm[J]. Chinese Journal of Geophysics, 2000, 43(3):402-410.
[26] Bancroft J C, Du Xiang. The Computation of Traveltimes When Assuming Locally Circular or Spherica Wavefronts[C]//76th Annual International Meeting. New Orleans:SEG, 2006:2231-2235.
[27] Kononov A, Gisolf D, Verschuur E. Application of Neural Networks to Travel-Times Computation[C]//77th Annual International Meeting. San Antonio:SEG, 2007:1785-1789.
[28] Bóna A, Slawinski M A, Smith P. Ray Tracing by Simulated Annealing:Bending Method[J]. Geophysics, 2009, 74(2):T25-T32.
[29] Huang X, Greenhalgh S. Linearized Formulations and Approximate Solutions for the Complex Eikonal Equation in Orthorhombic Media and Applications of Complex Seismic Traveltime[J]. Geophysics, 2018, 83(3):C115-C136.
[1] Ye Pei,Li Qingchun. Improvements of Linear Traveltime Interpolation Ray Tracing for the Accuracy and Efficiency [J]. Journal of Jilin University(Earth Science Edition), 2013, 43(1): 291-298.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] ZHAO Hong-guang,SUN Jing-gui, CHEN Jun-qiang,ZHAO Jun-kang, YAO Feng-liang,DUAN Zhan. The Genesis and Evolution of Orebearing Fluids of the Xiaoxinancha Goldbearing Porphyry Copper Deposit in Yanbian Area: H,O,C,S,Pb Isotope Tracing[J]. J4, 2005, 35(05): 601 -606 .
[2] CHI Bao-ming, YI Shu-ping, LI Zhi-jun, ZHOU Yan-zhang. Water Resources’ Artificial Regulation in Dalian Area[J]. J4, 2005, 35(05): 632 -635 .
[3] WANG Ming-chang,YING Shen,LI Lin. Multiscale Representation of Spatial Information on Voronoi Images[J]. J4, 2005, 35(04): 539 -0542 .
[4] LAI Ya-wen,WANG Lin-gen, XIAO Guo-shi,GAN Shu-cai, DUAN Guo-zheng. Study on the Occurrence States of the Platinum Group Elements in Basaltfrom Western Guizhou Province and Their Geological Significance[J]. J4, 2005, 35(05): 607 -610 .
[5] DU Chong, LIN Nian-feng, TANG Jie, LI Zhao-yang. RS-GIS Integration Research on Land Use Dynamic Changes in Western Plain of Jilin Province[J]. J4, 2005, 35(03): 366 -0372 .
[6] WANG Yu-sheng, MA Xiao-fan,WANG Yu-Liang, CHENG Hai-jun, LV Zhi-yong. A Study on Treatment of Municipal Wastewater by Hybrid Biological Reactor[J]. J4, 2005, 35(03): 373 -0377 .
[7] YANG De-ming, HE Zhong-hua, HUANG Ying-cong, ZHAO Liang, DAI Lin-na. Metamorphism Characteristics of Songduo Group in Menba Township Mozhugongka County, Tibet and the Discussion on Its Age[J]. J4, 2005, 35(04): 430 -0435 .
[8] ZHANG Jun,PANG Xiong-qi,JIANG Zhen-xue,SHAO Zhen-jun,CHEN Dong-xia,WANG Liang. Quantificational Prospection of Petroliferous Quality of Lithologic Reservoir,Dongying Depression[J]. J4, 2005, 35(06): 732 -0737 .
[9] XIAO Chang-lai,LIANG Xiu-juan, CUI Jian-ming, LAN Ying-ying,ZHANG Jun, LI Shu-lan, LIANG Rui-qi, ZHENG Ce. Whole Curve Matching Method for Aquifer Parameters Determination[J]. J4, 2005, 35(06): 751 -0755 .
[10] WANG Qiang, ZAN Shu-qin, JIN Li-yong,CHEN Jun. A New Oospecies,Dictyoolithus Gongzhulingensis,from theEarly Cretaceous Quantou Formation in the Central Jilin Province[J]. J4, 2006, 36(02): 153 -0157 .