A linear traveltime perturbation interpolation me-thod for seismic ray tracing
Li Tongyu1,3, Zhang Jianzhong1,2,3
1. Key Laboratory of Submarine Geosciences and Prospecting Techniques, Ministry of Education, Qingdao, Shandong 266100, China; 2. Laboratory of Marine Mineral Resources Evaluation and Detection Techniques, National Laboratory for Marine Science and Technology, Qingdao, Shandong 266061, China; 3. College of Marine Geosciences, Ocean University of China, Qingdao, Shandong 266100, China
Abstract:The linear traveltime interpolation (LTI) method assumes that traveltime is linear along a cell boundary,and the traveltime at any point on a cell boundary can be computed by interpolating the known traveltimes at adjacent grid nodes.In reality,however,the traveltime is not linear along a cell boundary,and the linear assumption may result in big calculation errors especially in large grid cells.To solve this problem,we calculate the traveltime with a linear traveltime perturbation interpolation (LTPI) method.In this method,the traveltime at a grid node is decomposed into two parts:a traveltime perturbation and a reference traveltime propagating in an equivalent homogeneous medium.The traveltime perturbation is much less than the reference traveltime.We assume that the traveltime perturbation is linear along the cell boundary,while the reference traveltime is nonlinear.We discretize the complex media with undulant topography and interfaces into a series of irregular cells.We also derive the LTPI formula for 2D irregular cells.Based on the LTPI,we propose a ray tracing method.Numerical experiments indicate that the LTPI is more accurate and stable than the LTI in computing wavefront traveltimes and ray paths.The LTPI performs better than the LTI in complex media.Under a certain accuracy requirement,the LTPI is more efficient than the LTI.
Zhang J,Shi T K,Zhao Y et al.Static corrections in mountainous areas using Fresnel-wavepath tomography.Journal of Applied Geophysics,2014,111(12):242-249.
[2]
Zhang Jianzhong,Liu Han,Zou Zhihui et al.Velocity modeling and inversion techniques for locating micro-seismic events in unconventional reservoirs.Journal of Earth Science,2015,26(4):495-501.
[3]
Julian B R,Gubbins D.Three-dimensional seismic ray tracing.Journal of Geophysics Research,1977,43(1):95-114.
[4]
Thurber C H,Ellsworth W L.Rapid solution of ray tracing problems in heterogeneous media.Bulletin of the Seismological Society of America,1980,70(4):1137-1148.
[5]
Vidale J E.Finite-difference calculation of traveltimes.Bulletin of the Seismological Society of America,1988,78(6):2062-2076.
[6]
张霖斌,许云,乌达巴拉等.地震反射波波前法射线追踪.石油地球物理勘探,1995,30(1):50-55.Zhang Linbin,Xu Yun,Wudabala et al.Ray tracing for seismic reflection wave.OGP,1995,30(1):50-55.
[7]
谢海兵,马在田,王燕琨.一种新程函方程差分法层析成像.石油地球物理勘探,1997,32(5):623-627.Xie Haibing,Ma Zaitian and Wang Yankun.Seismic tomography based on new finite difference scheme of eikonal equation.OGP,1997,32(5):623-627.
[8]
李永博,李庆春,吴琼等.快速行进法射线追踪提高旅行时计算精度和效率的改进措施.石油地球物理勘探,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.
[9]
Nakanishi I,Yamaguchi K.A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure.Earth Planets & Space,1986,34(2):195-201.
[10]
Moser T J.Shortest path calculation of seismic rays.Geophysics,1991,56(1):59-67.
[11]
王辉,常旭.基于图形结构的三维射线追踪方法.地球物理学报,2000,43(4):534-541.Wang Hui,Chang Xu.3-D ray tracing method based on graphic structure.Chinese Journal of Geophysics,2000,43(4):534-541.
[12]
张美根,程冰洁,李小凡等.一种最短路径射线追踪的快速算法.地球物理学报,2006,49(5):1467-1474.Zhang Meigen,Cheng Bingjie,Li Xiaofan et al.A fast algorithm of shortest path ray tracing.Chinese Journal of Geophysics,2006,49(5):1467-1474.
[13]
桑运云,李振春,张凯.抛物旅行时插值最短路径射线追踪.石油地球物理勘探,2013,48(3):403-409.Sang Yunyun,Li Zhenchun and Zhang Kai.Shortest path ray tracing based on parabolic travel-time interpolation.OGP,2013,48(3):403-409.
[14]
Asakawa E,Kawanaka T.Seismic ray tracing using linear traveltime interpolation.Geophysical Prospecting,1993,41(1):99-111.
[15]
张建中,陈世军,余大祥.最短路径射线追踪方法及改进.地球物理学进展,2003,18(1):146-150.Zhang Jianzhong,Chen Shijun,Yu Daxiang.Improvement of shortest path ray tracing method.Progress in Geophysics,2003,18(1):146-150.
[16]
张建中,陈世军.动态网络最短路径射线追踪.地球物理学报,2004,47(5):899-904.Zhang Jianzhong,Chen Shijun.A method of shortest path ray tracing with dynamic networks.Chinese Journal of Geophysics,2004,47(5):899-904.
[17]
黄月琴,张建中.基于波前传播时间插值的三维射线追踪算法.声学学报,2008,33(1):21-27.Huang Yueqin,Zhang Jianzhong.A three-dimensional sound ray tracing algorithm based on wavefront travel-times interpolation.Chinese Journal of Acoustics,2008,33(1):21-27.
[18]
Zhang J,Huang Y,Song L et al.Fast and accurate 3-D ray tracing using bilinear traveltime interpolation and the wave front group marching.Geophysical Journal International,2011,184(3):1327-1340.
[19]
Huang Y,Zhang J,Liu Q H.Three-dimensional GPRray tracing based on wavefront expansion with irregular cells.IEEE Transactions on Geoscience & Remote Sensing,2011,49(2):679-687.
[20]
李培明,梅胜全,马青坡.一种改进的双线性插值射线追踪方法.石油地球物理勘探,2013,48(4):553-558.Li Peiming,Mei Shengquan and Ma Qingpo.An improved bilinear interpolation travel-time ray-tracing method.OGP,2013,48(4):553-558.
[21]
黄翼坚,朱光明,敦敏.直达波旅行时线性插值算法.石油地球物理勘探,2010,45(2):225-229.Huang Yijian,Zhu Guangming and Guo Min.Direct wave travel linear interpolation algorithm.OGP,2010,45(2):225-229.
[22]
张东,张婷婷,乔友锋等.三维旅行时场B样条插值射线追踪方法.石油地球物理勘探,2013,48(4):559-566.Zhang Dong,Zhang Tingting,Qiao Youfeng et al.A 3-D ray tracing method based on B-spline traveltime interpolation.OGP,2013,48(4):559-566.
[23]
张婷婷,张东,邱达.三维层状介质中基于走时梯度的多次波射线追踪.石油地球物理勘探,2014,49(6):1097-1105.Zhang Tingting,Zhang Dong and Qiu Da.Multiple ray tracing in 3D layered media using the traveltime gradient method.OGP,2014,49(6):1097-1105.
[24]
张婷婷,邱达,张东.一种改进的三维旅行时梯度射线追踪方法.石油地球物理勘探,2016,51(5):916-923.Zhang Tingting,Qiu Da and Zhang Dong.A modified 3D traveltime gradient ray tracing method.OGP,2016,51(5):916-923.
[25]
Zhang J,Shi J,Song L P et al.Linear traveltime perturbation interpolation:a novel method to compute 3-D traveltimes.Geophysical Journal International,2015,203(1):548-552.
[26]
Shi T K,Zhang J Z,Huang Z L et al.A layer-stripping method for 3D near-surface velocity model building using seismic first-arrival times.Journal of Earth Science,2015,26(4):502-507.