[发明专利]自由表面多次波的压制方法和装置

说明书

本发明提供了一种自由表面多次波的压制方法和装置，其中，该方法包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取反射部分中的第一自由表面多次波；提取非反射部分中的第二自由表面多次波；根据第一单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。在本发明实施例中，所得到的压制自由表面多次波之后的第二单炮地震数据精度更高。

技术领域

本发明涉及地质勘探技术领域，特别涉及一种自由表面多次波的压制方法和装置。

背景技术

在地震资料处理过程中，多次波通常被视为一种相干噪声，特别是针对海上采集的地震数据，多次波发育尤为严重，处理过程中，如果不能采取有效的手段对多次波予以消除，最终的成像结果往往会造成成果数据出现地质假象，甚至会降低地震资料信噪比以及分辨率。多次波根据其产生机理，其种类非常复杂，针对海洋资料而言，能量最强的往往是自由表面多次波。

通常，通过采用多次波压制方案，可以有效的消除自由表面反射率对地震数据的影响。目前，常规多次波压制方法通常只需要提供震源处地震资料的波形特征及自由表面反射率，并不需要提取详细的地下地质信息，而随着地震勘探技术的不断进步，自由表面多次波压制方法也逐渐演变成了适应性的反演处理问题，即：首先，估算震源与自由表面的反射率；其次，将不含多次的数据视为反演剩余量。但是此方法受计算成本及精度的限制，无法应用于三维地震数据。

发明内容

本发明提供了一种自由表面多次波的压制方法和装置，以解决现有技术无法满足三维地震数据压制多次波的问题。

本发明实施例提供了一种自由表面多次波的压制方法，可以包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取所述反射部分中的第一自由表面多次波；提取所述非反射部分中的第二自由表面多次波；根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。

在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射波中的第一自由表面多次波。

在一个实施例中，按照以下公式提取所述反射波中的第一自由表面多次波：

其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。

在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。

在一个实施例中，按照以下公式提取所述反射部分中的第一自由表面多次波：

其中，表示第n次迭代后得到的第一自由表面多次波，表示所述反射部分中第 n次迭代去除自由表面多次波后的波形数据，D_r表示所述反射部分，D_s表示所述非反射部分。