首页|基于快速扫描法计算射线路径的走时层析成像方法/

基于快速扫描法计算射线路径的走时层析成像方法/

扫码查看
本文提出了一种梯度下降走时层析成像方法.本方法与目前流行的伴随状态法同为基于程函方程的走时层析成像方法,两者都克服了基于两点射线追踪的传统走时层析成像方法的低效问题,能够同时计算一震源对应的全部检波器对梯度的贡献.但本文所述方法避免了伴随状态法中复杂数学过程的引入,梯度由对弗雷歇导数的直接求解得到,因此更容易推导实现.对弗雷歇导数求解的关键是基于快速扫描法计算一个射线路径项:首先计算走时场,随后基于走时场计算波前面的传播方向,进一步基于波前面的传播方向计算射线路径项.每个步骤均可应用并行计算技术进行加速.因此,本文所述方法在保有与伴随状态法同样的高效、存储成本低和鲁棒性好的特点的同时,能够显式地计算出射线路径,推导和实现过程更加简单直观.
Travel time tomography by ray tracing using the fast sweeping method
This paper presents a gradient-descent travel time tomography method for solving the acoustic-type velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equation,enabling simultaneous calculation of contributions from all common-source receivers to the gradient.This overcomes the inefficiency inherent in conventional travel time tomography methods,which rely on a two-point ray tracing process.By directly calculating Fréchet derivatives,our method avoids the complex derivation processes associated with the adjoint-state method.The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term.Consequently,compared to the adjoint-state method,the proposed method can explicitly obtain the ray paths,resulting in a more concise and intuitive derivation process.Furthermore,our method retains the benefits of the adjoint-state method,such as speed,low memory usage,and robustness.This paper focuses on elucidating the principles and algorithms for calculating the ray-path term based on the fast sweeping method.The algorithms could be further speeded up by using parallel computational techniques.Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths,regardless of the complexity of the model and recording geometry.

tomographyvelocity model buildingray tracingfast-sweeping method

唐怀谷、谢仁军、吴怡、周长所、袁俊亮、秦玮

展开 >

中海油研究总院,北京 100028

层析成像 速度建模 射线追踪 快速扫描法

2024

应用地球物理(英文版)
中国地球物理学会

应用地球物理(英文版)

影响因子:1.01
ISSN:1672-7975
年,卷(期):2024.21(4)