[发明专利]用于海洋地震资料的二次震源衰减方法

权利要求书

1.一种用于海洋地震资料的二次震源衰减方法，其特征在于，包括：

步骤1：在包含二次震源水面直达波的时窗范围内，根据设定的振幅阈值获取单炮单缆地震数据中第一个强振幅的波峰，确定所述单炮单缆地震数据中二次震源的水面直达波的最小起始时间；

步骤2：在所述包含所述二次震源水面直达波的时窗范围内，输出水面直达波模型并去除所述水面直达波模型中二次震源的水面直达波；其中，所述输出水面直达波模型包括：利用水速动校正方法拉平所述单炮单缆地震数据中的水面直达波，通过多个相邻道叠加衰减随机噪音，然后输出所述水面直达波模型；

步骤3：根据所述二次震源的水面直达波的最小起始时间，在所述包含二次震源水面直达波的时窗范围内，截取出包含二次震源的地震数据，对所述包含二次震源的地震数据通过曲波变换方法重构出二次震源反射信号模型；其中，所述通过曲波变换方法重构出二次震源反射信号模型包括：对截取出的所述包含二次震源地震数据进行多尺度多方向二维曲波变换，得到曲波域系数；在曲波域采用局部阈值去噪方法，针对二维曲波变换后的每个尺度、每个方向的二次震源曲波系数分别选取一个阈值算子并进行阈值处理，得到每个尺度、每个方向下的曲波系数；分别针对每个尺度、每个方向下的所述曲波系数进行逆变换，重构出所述二次震源反射信号模型；

步骤4：将步骤2去除所述二次震源的水面直达波后的水面直达波模型减去步骤3重构的所述二次震源反射信号模型，得到消除二次震源后的地震数据。

2.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，所述根据设定的振幅阈值获取所述单炮单缆地震数据中第一个强振幅的波峰包括：选取所述单炮单缆地震数据中距离炮点最近的地震道，给定所述包含二次震源水面直达波的时窗范围，通过测试设定所述振幅阈值，拾取所述第一个强振幅的波峰。

3.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，在所述步骤2中，利用自适应相减方法去除所述水面直达波模型中二次震源的水面直达波。

4.根据权利要求3所述的用于海洋地震资料的二次震源衰减方法，其特征在于，还包括：选择所述单炮单缆地震数据中近偏移距的地震道集并重复所述步骤2，直至所述近偏移距的地震道集的二次震源的水面直达波的能量被完全消除。

5.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，所述截取出包含二次震源的地震数据包括：从所述单炮单缆地震数据中截取出所述二次震源的水面直达波的最小起始时间以后的地震数据。

6.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，所述阈值处理包括：比较每个尺度、每个方向的所述二次震源曲波系数与所述阈值算子，去除小于所述阈值算子的所述二次震源据曲波系数。

7.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，还包括：调整所述阈值算子的值，控制二次震源反射信号的去除量。

8.根据权利要求1所述的用于海洋地震资料的二次震源衰减方法，其特征在于，所述步骤4包括：将重构的所述二次震源反射信号模型与去除所述二次震源的水面直达波后的水面直达波模型进行自适应匹配后，用去除所述二次震源的水面直达波后的水面直达波模型减去所述二次震源反射信号模型得到所述消除二次震源的地震数据。