[发明专利]共炮检距域Beam叠前时间偏移成像方法、系统、介质及应用

权利要求书

1.一种共炮检距域Beam叠前时间偏移成像方法，其特征在于，所述Beam叠前时间偏移成像方法包括：

完成读入共炮检距域数据并划分子集，实际实验中先采用正演模拟出的二维SEGOverthrust高密度地震资料进行实验验证有效性，再用某地区实际三维地震资料进行实验；在读入地震数据之后，根据Beam中心间隔将共炮检距域记录划分为不同的子集，完成最初的数据输入工作；

进行不同子集平面波的计算，对于不同地震数据的子集，根据给定的射线参数范围和间隔，通过引入二阶校正项计算得到对应不同参数的局部平面波；

进行地震波旅行时的计算，在不同的Beam中心位置，计算对应每个成像点的地震波双程旅行时和水平射线参数，然后提取对应的振幅值根据Beam叠前时间偏移成像公式进行累加得到成像结果；

完成Beam叠前时间偏移成像，对地震资料中的每一个的共炮检距域地震记录都进行上述计算过程，直到所有共炮检距域记录计算完成；

所述Beam叠前时间偏移成像方法读入地震数据的单炮地震记录，二维SEG Overthrust模型的网格为640×375，水平采样间距为15m，垂直采样间距为7.5m；利用有限差分法正演了200炮单炮记录，每炮201道，炮间距和道间距均为20m，将炮记录抽取成共炮检距域地震记录；

所述Beam叠前时间偏移成像方法根据Beam中心间隔将单炮记录划分为不同的子集，二维SEG Overthrust模型数据中选择的中心间隔为240m；

Beam中心间隔△L是Beam叠前时间偏移中的重要参数，△L的选取同地震波的平均速度V_avg和主频f₀有关：

ΔL＝2V_avg/f₀；

所述共炮检距域Beam叠前时间偏移成像方法对于不同的数据子集，根据给定的射线参数范围和间隔，通过引入二阶校正项得到对应不同射线参数的局部平面波，局部平面波用下面的公式进行计算：

式中t表示截距时间，Φ表示常数因子，Δr’表示接收点到束中心的空间距离，h代表炮检距，r_m代表中心点位置，L代表Beam中心位置，表示共炮检距域地震记录，是用来提高平面波合成精度的二阶走时校正项；

所述Beam叠前时间偏移成像方法在不同的Beam中心位置，计算对应每个成像点的地震波双程走时和水平射线参数，已知Beam中心对应的震源和检波点到成像点的走时为t_s+t_d，近似计算L附近的震源和检波点到成像点的走时为：

t′_sd≈t_s+t_d+p_m·Δr′

Beam中心处射线参数的水平分量为：

每个成像点的地震波双程走时公式为：

式中t_s和t_d分别为震源和接收点到成像点的单程走时，V_rms为成像点处的均方根速度。

2.如权利要求1所述的共炮检距域Beam叠前时间偏移成像方法，其特征在于，所述Beam叠前时间偏移成像方法提取对应位置振幅值累加得到成像结果，Beam叠前时间偏移成像公式为：

其中，r表示成像点空间位置，A^*(L+h)A^*(L-h)为偏移加权函数，提取对应的振幅值根据上式累加得到成像结果

3.如权利要求2所述的共炮检距域Beam叠前时间偏移成像方法，其特征在于，所述Beam叠前时间偏移成像方法对于地震资料中每一个共炮检距记录重复计算过程，直到所有地震记录计算完成。

4.一种计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器执行时，使得所述处理器执行如权利要求1～3任意一项所述共炮检距域Beam叠前时间偏移成像方法。