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

权利要求书

1.一种自由表面多次波的压制方法，其特征在于，包括：

将所获取的第一单炮地震数据划分为反射部分以及非反射部分；

提取所述反射部分中的第一自由表面多次波；其中，按照以下公式提取所述反射部分中的第一自由表面多次波：

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

提取所述非反射部分中的第二自由表面多次波；其中，按照以下公式提取所述非反射部分中的第二自由表面多次波：

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

根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。

2.如权利要求1所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：

获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；

获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；

将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射部分中的第一自由表面多次波。

3.如权利要求2所述的方法，其特征在于，按照以下公式提取所述反射部分中的第一自由表面多次波：

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

4.如权利要求2所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：

获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；

计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；

计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；

计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。

5.如权利要求1所述的方法，其特征在于，所述提取所述非反射部分中的第二自由表面多次波包括：

获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；

获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；

将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射部分中的第二自由表面多次波。

6.如权利要求5所述的方法，其特征在于，按照以下公式提取所述非反射部分中的第二自由表面多次波：

其中，表示第n次迭代后得到的第二自由表面多次波，表示第n次迭代后得到的第五自由表面多次波，表示第n次迭代后得到的第六自由表面多次波。

7.如权利要求5所述的方法，其特征在于，所述提取所述非反射部分中的第二自由表面多次波包括：

获取所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；

计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第五自由表面多次波；

计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第六自由表面多次波；

计算所述第五自由表面多次波以及所述第六自由表面多次波的和，将所得到的和作为所述非反射部分中的第二自由表面多次波。