[发明专利]一种虚拟震源二维波前构建地震波走时计算方法

说明书

本发明公开了一种虚拟震源波前构建地震波走时计算方法，包括以下步骤：读入相关参数、速度模型；从震源沿着不同方向使用龙格库塔法追踪射线，并在射线追踪过程中依据相关判断条件插入新射线；通过相邻射线与相邻波前面上的离散点把计算空间划分为多个波前四边形；寻找出波前四边形内部网格点；通过虚拟震源法计算得出波前四边形内部网格点走时；完成所有波前四边形内部网格点走时计算。本发明通过采用虚拟震源法计算波前四边形内部网格点走时，提高了网格点走时计算精度，实现了一种高精度波前构建地震波走时计算方法。

技术领域

本发明涉及地震波走时计算领域，特别是一种二维波前构建地震波走时计 算方法。

背景技术

《吉林大学学报》(地球科学版)2008年第2期公开了韩复兴等“基于二维 三次卷积插值算法的波前构建射线追踪”，介绍了一种改进的波前构建地震波走 时计算方法，将二维三次卷积插值方法应用到了网格点走时计算中，以此提高 了获得的地震波走时计算精度以及算法的计算效率。并且通过均匀介质对基于 二维三次卷积插值算法的波前构建方法进行了误差分析，分析结果得到了比较 好的效果。

《计算物理》2008年第2期公开了韩复兴等“不同插值算法在波前构建射 线追踪中的应用与对比”，对比分析了邻近域插值法、双线性插值发、分片性插 值法以及二维三次卷积插值法在波前构建射线追踪中的应用效果，模型计算结 果反映出二维三次卷积插值方法相对于其它三种方法可以得到更加准确的射线 路径。

《地球物理学进展》2009年第5期公开了韩复兴等“波前构建法中网格点 相对定位及属性计算研究”，介绍了如何使用矢量叉乘法判断矩形网格节点与非 规则波前四边形的相对位置关系，并根据网格点与非规则波前四边形的位置关 系，给出了相对应的插值方法计算网格点的属性信息。并且通过均匀模型和计 算实例对提出的方法进行了验证，取得了良好的计算结果。

通过以上例子可以看出，现有二维波前构建地震波走时计算方法在一定程 度上能够提升计算精度，但是其中的插值方法并没有考虑地震波的传播规律， 提升的计算精度也有限。

发明内容

本发明所要解决的技术问题是提供一种虚拟震源波前构建地震波走时计算 方法，考虑了地震波在介质中的传播规律，在计算波前四边形内网格节点走时 过程中，使用了一种虚拟震源计算方法以代替原有的二维三次卷积插值方法， 提高了波前构建地震波走时计算方法的计算精度以及稳定性。

为解决上述技术问题，本发明采用的技术方案是：

一种高精度波前构建地震波走时计算方法，包括以下步骤：

步骤1：读入相关参数文件、速度模型，其中，所述参数文件包含速度模型 的网格点数、网格间距、震源位置、追踪步长、射线角度范围、射线采样间隔；

步骤2：追踪射线，并在追踪过程中插入新射线以此来保证射线的覆盖率。 追踪射线的本质是使用龙格库塔法求解运动学射线追踪方程组，如下式所示：

其中，x_i表示位置分量，p_i表示慢度分量，τ表示地震波走时，v表示地震波 速度。

步骤3：通过计算获取的射线上离散点的空间位置信息将计算空间划分为多 个波前四边形；

步骤4：通过矢量乘积法判断网格点与波前四边形的位置关系，找到波前四 边形所包含的网格点；

步骤5：通过波前四边形顶点相关信息计算出各个点所对应的虚拟震源位置， 并基于这些虚拟震源位置计算出波前四边形所包含网格点的地震波走时；

步骤6：完成所有网格节点的走时计算，并输出最终的走时计算结果。