[发明专利]一种复杂地形条件下的克希霍夫型波束偏移方法

摘要

本发明公开了一种复杂地形条件下的克希霍夫型波束偏移方法，包括以下步骤：读入相关参数文件、偏移速度模型及地震记录；将单炮地震记录分为若干不同的以窗中心为核心的时间域数据体，并对每一个数据体进行局部平面波分解；从炮点沿着不同方向追踪射线，并计算储存每条射线对应分段式射线束范围内网格节点的属性信息；从窗中心沿着不同方向追踪射线，并计算储存每条射线对应分段式射线束范围内网格节点的属性信息；选取炮点和窗中心相应的射线束对进行成像计算；累加所有射线束对的成像结果，得到偏移成像结果。本发明增加了射线束在模型浅层的覆盖率，提高了克希霍夫型波束偏移方法在复杂地形模型中浅层构造的成像效果。

主权项

1.一种复杂地形条件下的克希霍夫型波束偏移方法，其特征在于，包括以下步骤：步骤1：读入相关参数文件、偏移速度模型以及地震记录，其中，所述参数文件包含模型的网格点数、网格间距、参考频率、最大频率、初始波束宽度、地震记录采样点数、每一炮的道数和道间距；步骤2：将单炮地震记录按照距离分为若干不同的以窗中心为核心的时间域数据体，并对每一个数据体进行局部平面波分解；其中，V_avg为偏移速度场的平局值，f_min为参考频率，f_max为选取的频率最大值；步骤3：从炮点沿着不同方向追踪射线，并计算储存每条射线对应分段式射线束范围内网格节点的属性信息，所述属性信息包括走时和振幅；从炮点沿着不同方向发射射线，使用龙哥库塔法求解运动学射线追踪方程组得到中心射线信息；其中，x_i表示空间位置，p_i表示慢度，τ表示地震波走时，v表示离散点处的速度值；得出中心射线上离散点信息后，通过榜轴近似方法求取射线束覆盖范围内网格节点的属性信息，所述属性信息包括走时和振幅；通过分段式方法计算射线束的宽度：其中，Δa为相邻射线的角度差，σ为速度沿着射线路径的积分；步骤4：从窗中心沿着不同方向追踪射线，并计算储存每条射线对应分段式射线束范围内网格节点的属性信息，所述属性信息包括走时和振幅；同样，射线束的宽度通过步骤3中分段式方法进行计算；对窗中心对应的地震数据进行局部倾斜叠加，得到局部平面波分解结果；步骤5：选取炮点和窗中心相应的射线束对进行成像计算：其中，I_s(x)表示点x处的成像值，p_s表示炮点发出射线的慢度值，p_bc表示窗中心发出射线的慢度值，A表示权函数，D_s表示局部平面波分解结果；步骤6：累加所有射线束对的成像结果，得到最终的偏移成像结果。