一种改进的三维旅行时梯度射线追踪方法

最大旅行时梯度射线追踪(The Maximum Traveltime Gradient Ray Tracing,MTG)法利用三次B样条插值计算旅行时场的最大旅行时梯度,从而追踪初至波射线路径。但由于三次B样条插值是连续光滑的插值函数,在射线通过速度突变区域时,计算次级源点的旅行时和梯度会产生一定的误差。在分析MTG算法计算误差来源的基础上,提出了一种改进的B样条/线性联合插值的三维射线追踪(The Modified Maximum Traveltime Gradient Ray Tracing,简称MMTG)算法,即在速度均匀变化的区域采用B样条插值计算旅行时梯度确定下一个次级源点,在速度变化剧烈的区域利用线性插值法在网格界面上找出旅行时最小的入射点作为次级源点,从而保证射线在速度分界面上折射产生的不连续性。数值实验结果表明,MMTG算法在保持MTG算法优点的同时进一步提高了射线路径的计算精度,能够适应更复杂的速度介质模型。

基于Chebyshev多项式的三维射线追踪

采用快速匹配追踪(FMM)算法计算旅行时场;根据函数逼近和正交多项式的定义,利用Chebyshev正交多项式逼近旅行时场,从而使误差在最小二乘意义下达到最小,起到一定的平滑作用;然后直接对Chebyshev插值函数施以微分运算计算旅行时场的梯度,并最终得到射线路径。数值实验表明该方法提高了射线路径计算的精度和效率。