[发明专利]一种压制地震噪声的方法、装置及系统

权利要求书

1.一种压制地震噪声的方法，其特征在于，包括：

将原始地震数据进行曲波变换处理，获得曲波系数；

对所述曲波系数求取梯度获得局部波向量；

对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数；

对所述阈值处理后的曲波系数进行逆曲波变换，获得压制噪声后的地震数据。

2.根据权利要求1所述的压制地震噪声的方法，其特征在于，所述方法还包括：

获取绝对值大于等于第一预设阈值的曲波系数；

相应的，所述对所述曲波系数求取梯度获得局部波向量包括：对绝对值大于等于第一预设阈值的曲波系数求取梯度，获得局部波向量。

3.根据权利要求1所述的压制地震噪声的方法，其特征在于，所述对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数，包括：

对所述局部波向量的大小进行第一阈值处理，获得第一阈值处理后的局部波向量；

根据预先存储的局部波向量与曲波系数的索引关系，获得第一阈值处理后的局部波向量对应的曲波系数作为所述阈值处理后的曲波系数。

4.根据权利要求1所述的压制地震噪声的方法，其特征在于，所述对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数，包括：

根据对地震噪声与有效地震波的局部波向量数据特征差异性的分析结果，对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数。

5.根据权利要求1-4任一项所述的压制地震噪声的方法，其特征在于，所述对所述曲波系数求取梯度获得局部波向量，包括：

根据如下公式获得不同位置处的局部波向量：

其中，W_d(a,θ,b)表示地震数据d的曲波系数，a表示曲波系数的尺度参数，θ表示曲波系数的方向参数，b表示曲波系数的位置参数，表示在位置b上对曲波系数做梯度，i表示虚数单位，v_d(a,θ,b)表示位置b上的局部波向量。

6.根据权利要求1所述的压制地震噪声的方法，其特征在于，所述方法还包括：

根据所述局部波向量获取不同位置处的波数，将相同波数对应的曲波系数平方值进行叠加，获得同步压缩能量，对所述同步压缩能量进行第二阈值处理；

相应的，对所述局部波向量的大小进行第一阈值处理之后，所述阈值处理后的曲波系数包括：

选取所述第一阈值处理和第二阈值处理后的曲波系数作为所述阈值处理后的曲波系数。

7.根据权利要求6所述的压制地震噪声的方法，其特征在于，所述根据所述局部波向量获取不同位置处的波数，将相同波数对应的曲波系数平方值进行叠加，获得同步压缩能量，包括：

其中，T_d(v,b)表示同步压缩能量，表示局部波矢量的实部，δ表示delta函数，v表示不同位置处的波数。

8.一种压制地震噪声的装置，其特征在于，包括：

曲波变换模块，用于将原始地震数据进行曲波变换处理，获得曲波系数；

波向量计算模块，用于对所述曲波系数求取梯度获得局部波向量；

阈值处理模块，用于对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数；

逆曲波变换模块，用于对所述阈值处理后的曲波系数进行逆曲波变换，获得压制噪声后的地震数据。

9.根据权利要求8所述的压制地震噪声的装置，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

将原始地震数据进行曲波变换处理，获得曲波系数；

对所述曲波系数求取梯度获得局部波向量；

对所述局部波向量的大小进行第一阈值处理，获得阈值处理后的曲波系数；

对所述阈值处理后的曲波系数进行逆曲波变换，获得压制噪声后的地震数据。

10.一种压制地震噪声的系统，其特征在于，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1-7中任意一项所述方法的步骤。