一种基于改进程函方程的快速扫描地震波走时计算方法

本发明公开了一种基于改进程函方程的快速扫描地震波走时计算方法，涉及地震波走时计算领域，首先读入速度模型以及相关参数信息；然后进行震源初始化，设定慢度固定项，计算走时固定项，依据震源点位置，设置走时扰动项；接下来进行Gauss‑Seidel迭代计算，在交替的4个方向上，使用Gauss‑Seidel型迭代法求解因式分解的程函方程，并通过因果条件选择合适的走时扰动项，最后，将满足终止条件的走时计算结果输出，本发明提供了一种基于因式分解程函方程的高精度地震波走时计算方法，良好地解决了震源奇异性问题，保证有较高精度的同时，算法也具有较高的计算效率，提高后续地震勘探数据处理的准确度。