[发明专利]一种基于克希霍夫积分反偏移的层间多次波预测方法

说明书

一种基于克希霍夫积分反偏移的层间多次波预测方法，属于地震多次波压制领域。该方法首先基于输入的叠前时间偏移剖面，应用逆散射级数法预测出仅含有层间多次波的叠前时间偏移剖面，然后在此基础上利用时域克希霍夫反偏移算法最终获得层间多次波炮集记录。模型数据的实验结果显示，基于叠前时间偏移剖面的层间多次波预测有效减少了多次波预测的计算量，且通过克希霍夫反偏移创建层间多次波炮集记录其显著降低了实际问题的复杂程度，使层间多次波预测能够用于实际地震资料处理，是一种较为理想的层间多次波预测方法。

技术领域

本发明属于地震资料数据处理领域，是一种既可减少计算量又可保证精度的基于克希霍夫积分反偏移的层间多次波预测方法。

背景技术

由于地层当中存在强反射界面，地震波在界面之间经过多次反射形成层间多次波，若不加以去除将严重影响成像剖面的质量，甚至会误导后期地震资料的解释，因此预测与衰减层间多次波是地震资料数据处理的关键环节。

多次波压制主要是根据多次波与一次波性质上的差异，应用一定的数学处理方法分离或预测出多次波，并在此基础上衰减多次波。当前压制多次波的方法主要有两种：一种是滤波法，其主要是依据多次波与一次波视速度或周期性的差异，利用数学变换将多次波分离，达到压制多次波的效果，但是在多次波与一次波视速度相近的情况下，该方法无法达到理想的压制效果；另一种常用的多次波压制方法是预测减去法，该方法是依据多次波与一次波的关系预测出多次波，并在相应的炮集域或CMP域中减去相应的预测结果，以达到压制多次波的目的。

常规的多次波压制方法对于层间多次波无法起到良好的效果，其原因在于层间多次波视速度与一次波相近，并且周期性差异也不大，因此用常规的滤波方法难以识别与剔除层间多次波，而采用预测减去法计算量大，成本高，限制了其在实际资料处理中的应用。

发明内容

本发明要解决的技术问题在于提供一种基于克希霍夫积分反偏移的层间多次波预测方法，所述方法保证在较小的计算量下对多次波实现最有效的压制。总共分为两步实现对层间多次波的预测，第一步在成像剖面域预测一次波对应的多次波，第二步利用克希霍夫积分反偏移方法将其变换到炮集域。

本发明采取以下技术方案：

一种基于克希霍夫积分反偏移的层间多次波预测方法，其特征在于它具体包括以下步骤：

(1)给定一次波叠前时间剖面IMT，给定偏移速度v(x,t)，利用式(1)对一次波叠前时间偏移记录进行“时——深”转换，得到叠前深度偏移剖面IMD；

式中，IMD表示叠前深度偏移剖面，IMT表示叠前时间偏移剖面，v(x,t)为地震波速，x表示横向距离，t表示时间，z表示深度；

(2)输入叠前深度偏移剖面IMD，利用式(2)得到波数域的层间多次波预测结果IMK_m；

式中，x表示横向距离，k_z为波数，z′₁、z′₂、z′₃为深度，IMK_m为波数域的层间多次波预测结果；

(3)利用式(3)对IMK_m作反傅里叶变换，得到深度域层间多次波剖面IMD_m；

式中x表示横向距离，z表示深度，k_z为波数

(4)利用式(4)对深度域层间多次波剖面IMD_m做“深——时”转换，最终得到层间多次波的叠前时间偏移剖面IMT_m；