[发明专利]一种多源混合地震数据表层相关多次波压制方法

说明书

本发明涉及一种多源混合地震数据表层相关多次波压制方法，其步骤：对于给定的时空域地震数据，按炮点、检波点和偏移距信息进行排列；对排列后地震数据每一道进行傅里叶变换，得到频率域的切片数据；利用最小二乘方法，求取混合算子的逆；取频率域的切片数据与混合算子的逆相乘，得到预测的多次波；对频率域预测得到的多次波进行傅里叶逆变换，得到时空域预测的多次波记录；利用最小能量准则，将得到的多次波记录从原始给定的多源混合地震数据中自适应减去。本发明充分利用多源地震数据的信息实现多次波压制，从而有效提高地震数据深层信噪比，提高后期处理方法精度。

技术领域

本发明涉及一种地震数据多次波压制方法，特别是关于一种基于反馈迭代理论的多源混合地震数据表层相关多次波压制方法。

背景技术

随着新的硬件和软件技术的引入，地震采集技术实现了长久以来的一些设想，带来了新的采集理念，地震资料的质量和数量也有了很大的变化。多震源同时激发采集系统就是代表这一趋势的新技术，其与传统地震方法最大的区别在于多个震源同时激发，检波器连续接收信号，从时间域的间断、逐炮接收，到连续不断地采集资料，这是多震源同时激发采集技术在数据采集方面的重大变革，极大地提高了数据的获取量和采集效率。

多源混合地震采集技术通过降低相邻炮激发的等待时间来降低采集成本，或者通过增加相同勘探时间内的震源数目来提高采集数据的质量，可以较好的提高地震数据采集的效率。该方法采集得到的地震记录是由两个或多个常规采集地震记录组成。如果直接对多震源地震数据进行处理，可以明显地提高处理效率，特别是可以极大地改善叠前地震数据偏移成像和反演的效率。但是，多源地震数据和常规采集地震数据类似，也存在严重的多次波问题，多次波的存在对于多源地震数据的偏移成像和反演是一个严重的干扰，因此多次波压制对于多源地震数据直接处理的算法极为重要。

发明内容

针对上述问题，本发明的目的是提供一种多源混合地震数据表层相关多次波压制方法，其充分利用多源地震数据的信息实现多次波压制，从而有效提高地震数据深层信噪比，提高后期处理方法精度。

为实现上述目的，本发明采取以下技术方案：一种多源混合地震数据表层相关多次波压制方法，其包括以下步骤：1)对于给定的时空域地震数据，按炮点、检波点和偏移距信息进行排列；2)对排列后地震数据每一道进行傅里叶变换，得到频率域的切片数据P_bl；3)利用最小二乘方法，求取混合算子的逆；4)取频率域的切片数据P_bl与混合算子的逆相乘，得到预测的多次波M_bl；5)对频率域预测得到的多次波进行傅里叶逆变换，得到时空域预测的多次波记录；6)利用最小能量准则，将得到的多次波记录从原始给定的多源混合地震数据中自适应减去；7)根据深层多次波压制的效果，判断是否仍有多次波存在，若仍有，则返回步骤2)；反之，则输出多次波压制后的多源混合地震记录。

进一步，所述步骤3)中，混合算子的逆(Γ_bl)-¹求取方法为：设多源混合地震数据混合算子为Γ_bl，采用最小二乘理论对其进行求逆处理，该算子的最小二乘逆为：

式中，上标H代表矩阵的共轭转置，μ为稳定求逆过程的阻尼因子，I为单位矩阵。

进一步，所述步骤4)中，预测的多次波M_bl为：

M_bl＝P_0,bl(Γ_bl)^-1P_bl

其中，P_0,bl为多源混合地震数据一次波记录。