[发明专利]多次波干扰的压制方法、存储介质和计算机设备

说明书

本发明提供一种多次波干扰的压制方法、存储介质和计算机设备。该方法包括：获取原始地震数据，通过预处理得到共反射点道集数据；利用逆散射级数法得到层间多次波模型数据；在F‑K域，利用基于层间多次波模型数据设计的自适应滤波器对共反射点道集数据进行滤波，将共反射点道集数据分离成含有部分一次波的多次波数据和未包含多次波的一次波数据，通过F‑K反变换变换到时间空间域；在时间空间域，利用层间多次波模型数据对含有部分一次波的多次波数据执行曲波域匹配相减，得到剩余的一次波数据；将剩余的一次波数据与未被多次波干扰的一次波数据进行相加融合，得到压制层间多次波干扰后的地震数据。本实施例能够较好的压制层间多次波的干扰。

技术领域

本发明涉及地球物理勘探技术领域，尤其涉及一种多次波干扰的压制方法、存储介质和计算机设备。

背景技术

多次波干扰普遍存在于地震勘探中，其中，海上和陆上某些地区尤为严重。在海洋地震勘探中，由于海水表面强反射层的存在，反射回来的地震波很大一部分会进行又一次的表层反射，因此，在海洋数据中发育着丰富的自由表面多次波。而陆上地震勘探相对于海上地震勘探而言，除了发育自由表面多次波外，还发育着大量的层间多次波。无论是在海洋勘探还是陆地勘探中，多次波对地震资料都有多方面的影响。例如，多次波会降低地震资料的信噪比进而影响最终的成像精度；当目的层反射较弱时，多次波还会影响地质认识，从而影响勘探部署。因此，正确的识别和有效压制多次波成为地震勘探中的一个重要问题。

目前，压制多次波的方法主要包括两大类：滤波法和基于波动方程的预测法。滤波法主要利用了多次波的周期性和可分离性，常用的方法有2D SRME(Surface RelatedMultiple Elimination)、FK变换、拉东变换、预测反褶积等。但是这些算法受地震速度影响较大，近道去除不彻底，一般在CRP道集(共反射点道集)上，不能满足RTM(逆时偏移)成像需要。波动方程预测法利用波动方程或数据驱动方法预测多次波，并将预测得到的多次波通过自适应相减的方法从地震记录中减去以达到压制多次波的目的。该方法代表了多次波压制方法的发展趋势，主要包括反馈迭代法和逆散射级数法。无论是滤波法还是波动方程预测法，在层间多次波的压制中都有所应用。Jakubowicz(1998)方法长期以来一直是层间多次波的主要预测算法，由于其需要对不同的多次波产生的位置进行分层，使得整个预测过程对于实际生产是不可行的。其他方法，如Weglein等人使用从Lippmann-Schwinger方程导出的逆散射级数和Verschuur等人基于层位近似的方法也在层间多次波衰减的背景下进行了探索。此外，van der Neut和Wapenaar曾尝试使用Marchenko算法策略来压制层间多次波。

在地震数据处理过程中，层间多次波往往难以被完全压制，因为在某些情况下层间多次波对一次波有强烈的干涉作用，进而影响了油气藏解释。表面多次波压制已经成为了标准处理流程，而层间多次波的压制技术仍然不成熟，原因如下：首先，由于层间多次波的预测过程大多是数据驱动的，其典型实现本质上比表面多次波计算成本高出大约两到三个数量级，因此在大规模的三维数据上应用计算量令人望而却步，使得现有的算法通常很难应用于实际生产；第二，数据驱动的方法通常对数据的质量要求非常高，如果预处理的效果不理想，最终层间多次波的压制效果也很难满足生产需求；第三，从原始数据中难以减去层间多次波，因为多次波的振幅通常较弱，而且与一次波强烈的干涉在一起，在减法运算的过程中算法选择不当很容易造成一次波的损害或者进而引入新的噪声。

发明内容

本发明的主要目的是提供一种多次波干扰的压制方法、存储介质和计算机设备，以有效压制地震数据中的层间多次波干扰。