石油地球物理勘探2025,Vol.60Issue(5):1160-1167,8.DOI:10.13810/j.cnki.issn.1000-7210.20240489
应用三阶WENO格式和Lax-Friedrichs算法求解因式分解程函方程
Factored eikonal equation solver based on third-order WENO scheme and Lax-Friedrichs fast sweeping method
摘要
Abstract
Seismic wave traveltime computation is a critical foundational technology for seismic migration,to-mography and subsurface static correction,etc.The third-order weighted essentially non-oscillatory(WENO)finite difference scheme and the Lax-Friedrichs method are used to solve the factored eikonal equation to im-prove seismic wave traveltime accuracy.First,the general eikonal equation is factorized to avoid singularities near the source point.Then,the factored eikonal equation is rewritten in the Lax-Friedrichs traveltime iterative scheme,and the third-order WENO difference scheme is employed to replace the partial differential items of traveltime along the horizontal and vertical directions in the iterative scheme.Finally,the fast-sweeping itera-tive algorithm is utilized to solve the factored iterative scheme of traveltime.Calculation accuracy of seismic wave traveltime can be improved effectively due to the adoption of the factorization equation and high-order WENO difference scheme.The numerical solutions of the homogeneous medium,and 2D and 3D constant gra-dient velocity models show that compared with the conventional eikonal equations,the proposed traveltime cal-culation method significantly improves the calculation accuracy of seismic wave travel time.关键词
旅行时计算/三阶WENO格式/Lax-Friedrichs算法/因式分解/程函方程Key words
traveltime computation/third-order WENO scheme/Lax-Friedrichs method/factorization method/eikonal equation分类
天文与地球科学引用本文复制引用
黄光南,康婕,尹玉茹,陶学佳,李瑞华,吕安琪..应用三阶WENO格式和Lax-Friedrichs算法求解因式分解程函方程[J].石油地球物理勘探,2025,60(5):1160-1167,8.基金项目
本项研究受国家自然科学基金项目"一般各向异性介质地震qP波,qSV波和qSH波的有限差分走时计算方法研究"(42064007)和教育部春晖计划项目"基于弱各向异性假设的TI介质高精度地震qP波走时层析成像方法研究"(HZKY20220386)联合资助. (42064007)