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

摘要

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

主权项

1.一种虚拟震源二维波前构建地震波走时计算方法，其特征在于，包括以下步骤：步骤1：读入相关参数文件、速度模型，其中，所述参数文件包含速度模型的网格点数、网格间距、震源位置、追踪步长、射线角度范围、射线采样间隔；步骤2：追踪射线，并在追踪过程中插入新射线以此来保证射线的覆盖率。追踪射线的本质是使用龙格库塔法求解运动学射线追踪方程组，如下式所示：其中，xi表示位置分量，pi表示慢度分量，τ表示地震波走时，v表示地震波速度。步骤3：通过计算获取的射线上离散点的空间位置信息将计算空间划分为多个波前四边形；步骤4：通过矢量乘积法判断网格点与波前四边形的位置关系，找到波前四边形所包含的网格点；步骤5：通过波前四边形顶点相关信息计算出各个点所对应的虚拟震源位置，并基于这些虚拟震源位置计算出波前四边形所包含网格点的地震波走时；步骤6：完成所有网格节点的走时计算，并输出最终的走时计算结果。