卢江波, 方志. 2014: 一种线性走时插值射线追踪改进算法. 地震学报, 36(6): 1089-1100. DOI: 10.3969/j.issn.0253-3782.2014.06.010
引用本文: 卢江波, 方志. 2014: 一种线性走时插值射线追踪改进算法. 地震学报, 36(6): 1089-1100. DOI: 10.3969/j.issn.0253-3782.2014.06.010
Lu Jiangbo, Fang Zhi. 2014: An improved ray-tracing algorithm based on linear travel-time interpolation. Acta Seismologica Sinica, 36(6): 1089-1100. DOI: 10.3969/j.issn.0253-3782.2014.06.010
Citation: Lu Jiangbo, Fang Zhi. 2014: An improved ray-tracing algorithm based on linear travel-time interpolation. Acta Seismologica Sinica, 36(6): 1089-1100. DOI: 10.3969/j.issn.0253-3782.2014.06.010

一种线性走时插值射线追踪改进算法

An improved ray-tracing algorithm based on linear travel-time interpolation

  • 摘要: 针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍.

     

    Abstract: In order to solver for the problem that the original LTI algorithm could not trace the reverse propagation ray, several linear travel-time interpolation (LTI for short) improved algorithms, such as extension-compaction LTI algorithm, loop computation LTI algorithm, the shortest path ray tracing algorithm with dynamic networks, have been presented, but the computational efficiency of these algorithms are low. After analyzing these improved algorithms, this paper presented a new improved shortest path ray tracing algorithm with dynamic networks. According to the law of wave propagation and the basic equation of LTI, a large number of redundancy node calculation are excluded, and the traditional binary heap sort algorithm was used to manage node of wavefront array. The numerical examples show that, the improved algorithm presented in this paper has the highest computational efficiency among all of improved algorithms; its calculation efficiency is about 4.5—30 times of the shortest path ray tracing algorithm with dynamic networks, and about 2—6.5 times of the original LTI algorithm, and about 3.5—15 times of the shortest path ray tracing algorithm with dynamic networks when the traditional binary heap sort algorithm is also used.

     

/

返回文章
返回