[发明专利]一种平面波叠前深度偏移方法

说明书

技术领域

本发明涉及石油勘探反射波地震数据处理过程中的叠前偏移成像技术，具体的说是一种在不降低成像质量的前提下，提高叠前深度偏移处理速度的基于稀疏频率和大步长波场延拓的平面波叠前深度偏移方法。

背景技术

在石油勘探反射波地震数据处理过程中对复杂地质构造成像时，波动方程叠前深度偏移方法具有许多积分法无法比拟的优点，如没有高频近似、理论完整、计算合理、容易处理振幅问题、对复杂构造成像精度高等。其缺点是计算效率低。为了提高计算效率，Taner提出用独立炮记录抽出的共接收点道集，在地表合成平面波震源记录的方法。Berkhout进一步提出面炮记录叠前深度偏移方法，该方法大大减小了叠前数据量，其实质是一种基于平面波合成的偏移方法。该方法相对于单炮偏移计算效率成倍提高，但是计算量依然较大，在生产中广泛应用的成本较高。

发明内容

本发明是提供一种基于稀疏频率和大步长波场延拓，可将平面波叠前深度偏移的计算效率成倍提高的平面波叠前深度偏移方法。

本发明具体步骤包括：

1)野外采集的地震数据，通过常规预处理流程，得到统一基准面上的叠前炮集数据；

2)根据地震资料的复杂程度确定稀疏频率采样系数α和延拓大步长NΔz，所述Δz为通常延拓步长，N为延拓大步长相对于通常延拓步长的倍数；

步骤2所述的稀疏频率采样系数α是为重采样后的频率个数为重采样前的一半至三分之一：

步骤2所述的N为5至10。

步骤2所述的N在地震资料复杂时为5。

3)对频率进行稀疏采样；

利用步骤2中给出的采样系数α，按如下公式对频率进行重采样：

f_m＝α^-mf₁，m＝0，1，…，M

其中f₁为偏移的最小频率，f_n为偏移的最大频率，M为重采样后的频率个数，

4)对重采样后得到的每一个频率按照平面波叠前深度偏移方法以NΔz的大步长进行波场延拓，得到下一个延拓步的波场；

5)对相邻两个大延拓步之间的波场进行线性相位移插值，得到这两个延拓步中间每一个以小延拓步长Δz进行延拓的波场；

深度位置z₁到深度z₁+NΔz之间以小步长step1＝Δz延拓得到的波场可由如下公式计算：

P(x,z1+nΔz,ω)=(N-nN)Pd(x,z1+nΔz,ω)+(nN)Pu(x,z1+nΔz,ω)]]>