Traveltime calculation based on linear interpolation in hybrid meshes for rugged topographical conditions
Wang Qi1, Zhu Pan1, Ye Pei2, Li Qin3, Li Qingchun1
1. College of Geology Engineering and Geomatics, Chang'an University, Xi'an, Shaanxi 710054, China; 2. China Aero Geophysical Survey and Remote Sensing Center, Land and Resources, Beijing 100083, China; 3. College of Geological and Environmental Sciences, Xi'an University of Science and Technology, Xi'an, Shaanxi 710054, China
Abstract:In order to improve the accuracy and efficiency of traveltimes calculation under complex structures with rugged topographical conditions,we propose a new LTI method based on hybrid meshes.In rectangle meshes and irregular quadrilateral meshes,we deduce new formulas for local travetime calculation in hybrid meshes based on formulas only for local traveltime calculation in rectangle meshes.The formulas are proved by the tests.Then,we promote this new method into multi-stage LTI method in hybrid meshes.Numerical test results show that the proposed method can achieve accurate traveltimes and raypaths in the sharply rugged topography and complex near-surface.
王琦, 朱盼, 叶佩, 李勤, 李庆春. 起伏地表地震波旅行时混合网格线性插值射线追踪计算方法[J]. 石油地球物理勘探, 2018, 53(1): 35-46.
Wang Qi, Zhu Pan, Ye Pei, Li Qin, Li Qingchun. Traveltime calculation based on linear interpolation in hybrid meshes for rugged topographical conditions. Oil Geophysical Prospecting, 2018, 53(1): 35-46.
Asakawa E and Kawanaka T.Seismic ray tracing using linear traveltime interpolation.Geophysical Prospecting,1993,41(4):99-111.
[2]
赵改善,郝守玲.基于旅行时线性插值的地震射线追踪算法.石油物探,1998,32(2):14-24.Zhao Gaishan,Hao Shouling.Seismic ray tracing algorithm based on the linear traveltime interpolation.GPP,1998,32(2):14-24.
[3]
聂建新,杨慧珠.地震波旅行时二次/线性联合插值法.清华大学学报,2003,43(11):1495-1498.Nie Jianxin,Yang Huizhu.Quadrarical linear travel time interpolation of seismic ray-tracing.Journal of Tsinghua University,2003,43(11):1495-1498.
[4]
张赛民,周竹生,陈灵君等.对旅行时进行抛物型插值的地震射线追踪方法.地球物理学进展,2007,22(1):43-48.Zhang Saimin,Zhou Zhusheng,Chen Lingjun et al.Seismic ray tracing method of applying parabolic interpolation to travel time.Progress in Geophysics,2007,22(1):43-48.
[5]
张东,谢宝莲,杨艳等.一种改进的线性旅行时插值射线追踪算法.地球物理学报,2009,52(1):200-205.Zhang Dong,Xie Baolian,Yang Yan et al.A ray tra-cing method based on improved linear traveltime interpolation.Chinese Journal of Geophysics,2009,52(1):200-205.
[6]
Rawlinson N,Sambridge M.Multiple reflection and transmission phases in complex layered media using multistage fast marching method.Geophysics,2004,69(5):1338-1350.
[7]
De Kool M,Rawlinson N,Sambridge M.A practical grid based method for tracking multiple refraction and reflection phases in three-dimensional heterogeneous media.Geophysical Journal International,2006,167(1):253-270.
[8]
唐小平,白超英.最短路径算法下三维层状介质中多次波追踪.地球物理学报,2009,52(10):2635-2643.Tang Xiaoping,Bai Chaoying.Multiple ray tracing within 3-d layered media with the shortest path me-thod.Chinese Journal of Geophysics,2009,52(10):2635-2643.
[9]
赵瑞,白超英.复杂层状模型中多次波快速追踪——一种基于非规则网格的最短路径算法.地震学报,2010,32(4):433-444.Zhao Rui,Bai Chaoying.Fast multiple ray tracing within complex layered media:the shortest path method based on irregular grid cells.Acta Seismologica Sinica,2010,32(4):433-444.
[10]
岳玉波.起伏地表下初值射线追踪的实现.勘探地球物理进展,2007,30(5):388-391.Yue Yubo.Realization of initial value ray tracing for rugged topography.Progress in Exploration Geophy-sics,2007,30(5):388-391.
[11]
孙章庆,孙建国.波前快速推进法起伏地表地震波走时计算.勘探地球物理进展,2007,30(5):92-97.Sun Zhangqing,Sun Jianguo.Traveltime computation using fast marching method from rugged topography.Progress in Exploration Geophysics,2007,30(5):92-97.
[12]
孙章庆,孙建国,韩复兴.复杂地表条件下基于线性插值和窄带技术的地震波走时计算.地球物理学报,2009,52(11):2846-2853.Sun Zhangqing,Sun Jianguo and Han Fuxing.Traveltimes computation using linear interpolation and narrow band technique under complex topographical conditions.Chinese Journal of Geophysics,2009,52(11):2846-2853.
[13]
Lelièvre P G,Farquharson C G and Hurich C A.Computing first-arrival seismic traveltimes on unstructured 3-D tetrahedral grids using the fast marching method.Geophysical Journal International,2011,184(2):885-896.
[14]
Sun Jianguo,Sun Zhangqing and Han Fuxing.A finite difference scheme for solving the eikonal equation including surface to topography.Geophysics,2011,76(4):53-56.
[15]
Lan Haiqing and Zhang Zhongjie.Topography-dependent eikonal equation and its solver for calculating first-arrival traveltimes with an irregular surface.Geophysical Journal International,2013,193(2):1010-1026.
[16]
赵后越,张美根.起伏地表条件下各向异性地震波最短路径射线追踪.地球物理学报,2014,57(9):2910-2917.Zhao Houyue,Zhang Meigen.Tracing seismic shortest path rays in anisotropic medium with rolling surface.Chinese Journal of Geophysics,2014,57(9):2910-2917.
[17]
孙章庆.起伏地表条件下的地震波走时与射线路径计算[学位论文].吉林长春:吉林大学,2011.Sun Zhangqing.The Seismic Traveltimes and Ray Path Computation under Undulating Earth's Surface Condition[D].Jilin University,Changchun,Jilin,2011.
[18]
李永博,李庆春,吴琼等.快速行进法射线追踪提高旅行时计算精度和效率的改进措施.石油地球物理勘探,2016,51(3):467-473.Li Yongbo,Li Qingchun,Wu Qiong et al.Improved fast marching method for higher calculation accuracy and efficiency of traveltime.OGP,2016,51(3):467-473.
[19]
王珺,李永庆.遗传算法和正交时频原子相结合的地震记录快速匹配追踪.石油地球物理勘探,2016,51(5):881-888,893.Wang Jun,Li Yongqing.Seismic trace fast matching pursuit based on genetic algorithm and orthogonal time-frequency atom.OGP,2016,51(5):881-888,893.
[20]
秦宁,王延光,杨晓东等.非水平地表高斯束叠前深度偏移及山前带应用实例.石油地球物理勘探,2017,52(1):81-86.Qi Ning,Wang Yanguang,Yang Xiaodong et al.Gaussian bean prestack depth migration for undula-ting-surface area in piedmont zone.OGP,2017,52(1):81-86.
[21]
黄国娇,白超英.二维复杂层状介质中地震多波旅行时联合反演成像.地球物理学报,2010,53(12):2972-2981.Huang Guojiao,Bai Chaoying.Simultaneous inversion with multiple traveltimes within 2-d complex layered media.Chinese Journal of Geophysics,2010,53(12):2972-2981.
[22]
桑运云,孙晓军.起伏地表下基于抛物插值的最短路径射线追踪.石油物探,2014,53(2):1441-1448.Sang Yunyun,Sun Xiaojun.Shortest path raytracing based on parabolic traveltime interpolation in irregular topography.GPP,2014,53(2):1441-1448.
[23]
侯爵,张忠杰,兰海强等.起伏地表下地震波传播数值模拟方法研究进展.地球物理学进展,2014,29(2):488-497.Hou Jue,Zhang Zhongjie,Lan Haiqiang et al.Progress in numerical simulation of seismic wave propagation under an undulating surface.Progress in Geophysics,2014,29(2):488-497.
[24]
孙章庆,孙建国,韩复兴.基于线性插值和窄带技术的走时计算方法.石油地球物理勘探,2009,44(4):436-441.Sun Zhangqing,Sun Jianguo and Han Fuxing.Travel-time computation based on linear interpolation and narrow-band technique.OGP,2009,44(4):436-441.
[25]
叶佩,李庆春.旅行时线性插值射线追踪提高计算精度和效率的改进方法.吉林大学学报(地球科学版),2013,43(1):291-298.Ye Pei,Li Qingchun.Improvements of linear traveltime interpolation ray tracing for the accuracy and efficiency.Journal of Jilin University (Earth Science Edition),2013,43(1):291-298.
[26]
宋御杰.基于混合网格剖分的起伏地表旅行时计算方法研究[学位论文].陕西西安:长安大学,2015.Song Yujie.The Travel-time Computation Method Study Under Undulating Earth's Surface Based on Hybrid mesh[D].Chang'an University,Xi'an,Shaanxi,2015.
[27]
Kim S.3-D eikonal solvers:First-arrival traveltimes.Geophysics,2002,67(4):1225-1231.