An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation  

作　　者：Bao-Ping Qiao Qing-Qing Li Wei-Guang He Dan Zhao Qu-Bo Wu 

机构地区：[1]CNNC Beijing Research Institute of Uranium Geology,Beijing 100029,China [2]National Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China),Qingdao 266580,Shangdong,China [3]SINOPEC Geophysical Research Institute,Nanjing,211103,Jiangsu,China

出　　处：《Petroleum Science》2024年第1期195-205,共11页石油科学（英文版）

基　　金：The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04;National Natural Science Foundation of China(42074132,42274154).

摘　　要：3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

关 键 词：3D eikonal equation Accurate traveltimes Global fast sweeping 3D inhomogeneous media Adaptive finite-difference method 

分 类 号：P631.4[天文地球—地质矿产勘探]