[发明专利]基于曲波变换和聚类的地震数据去噪方法和装置

权利要求书

1.一种基于曲波变换和聚类的地震数据去噪方法，其特征在于，所述方法包括：

获取第一地震数据，将所述第一地震数据进行曲波变换，得到第一曲波变换系数；

对所述第一曲波变换系数进行去噪处理，得到第二曲波变换系数；

对所述第二曲波变换系数进行聚类处理，得到第三曲波变换系数；

所述对所述第二曲波变换系数进行聚类处理，得到第三曲波变换系数包括：将所述第二曲波变换系数中密度相连的核心点和边界点生成聚类簇，去除所述第二曲波变换系数中孤立的点，恢复所述聚类簇中所述核心点邻域内点的数值，从而得到所述第三曲波变换系数；

将所述第三曲波变换系数进行逆曲波变换，得到第二地震数据。

2.根据权利要求1所述的基于曲波变换和聚类的地震数据去噪方法，其特征在于，所述将所述第一地震数据进行曲波变换，得到第一曲波变换系数包括：

根据下式计算所述第一曲波变换系数：

其中，j,l,k分别为尺度、方向和位移，为变换的基函数，C(j,l,k)为所述第一曲波变换系数，f(x)为所述第一地震数据，x为所述第一地震数据的时间和空间参数。

3.根据权利要求1所述的基于曲波变换和聚类的地震数据去噪方法，其特征在于，所述对所述第一曲波变换系数进行去噪处理，得到第二曲波变换系数包括：

去除所述第一曲波变换系数中绝对值小于预设阈值的系数，保留所述第一曲波变换系数中绝对值大于预设阈值的系数，从而得到所述第二曲波变换系数。

4.根据权利要求3所述的基于曲波变换和聚类的地震数据去噪方法，其特征在于，所述去除所述第一曲波变换系数中绝对值小于预设阈值的系数，保留所述第一曲波变换系数中绝对值大于预设阈值的系数，从而得到所述第二曲波变换系数包括：

根据下式计算所述第二曲波变换系数：

其中，为所述第二曲波变换系数，τ_j为尺度j下的所述预设阈值，l为方向，k为位移。

5.根据权利要求1所述的基于曲波变换和聚类的地震数据去噪方法，其特征在于，所述将所述第三曲波变换系数进行逆曲波变换，得到第二地震数据包括：

根据下式计算所述第二地震数据：

其中，f₁(x)为所述第二地震数据，为所述第三曲波变换系数，为变换的基函数，j,l,k分别为尺度、方向和位移，x为所述第一地震数据的时间和空间参数。

6.一种基于曲波变换和聚类的地震数据去噪装置，其特征在于，所述装置包括：

曲波变换单元，用于获取第一地震数据，将所述第一地震数据进行曲波变换，得到第一曲波变换系数；

去噪处理单元，用于对所述第一曲波变换系数进行去噪处理，得到第二曲波变换系数；

聚类处理单元，用于对所述第二曲波变换系数进行聚类处理，得到第三曲波变换系数；

逆曲波变换单元，用于将所述第三曲波变换系数进行逆曲波变换，得到第二地震数据；

所述聚类处理单元包括：

将所述第二曲波变换系数中密度相连的核心点和边界点生成聚类簇，去除所述第二曲波变换系数中孤立的点，恢复所述聚类簇中所述核心点邻域内点的数值，从而得到所述第三曲波变换系数。

7.根据权利要求6所述的基于曲波变换和聚类的地震数据去噪装置，其特征在于，所述曲波变换单元包括：

根据下式计算所述第一曲波变换系数：

其中，j,l,k分别为尺度、方向和位移，为变换的基函数，C(j,l,k)为所述第一曲波变换系数，f(x)为所述第一地震数据，x为所述第一地震数据的时间和空间参数。

8.根据权利要求6所述的基于曲波变换和聚类的地震数据去噪装置，其特征在于，所述去噪处理单元包括：

去除所述第一曲波变换系数中绝对值小于预设阈值的系数，保留所述第一曲波变换系数中绝对值大于预设阈值的系数，从而得到所述第二曲波变换系数。

9.根据权利要求8所述的基于曲波变换和聚类的地震数据去噪装置，其特征在于，所述去噪处理单元还包括：

根据下式计算所述第二曲波变换系数：

其中，为所述第二曲波变换系数，τ_j为尺度j下的所述预设阈值，l为方向，k为位移。