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

说明书

本发明提供了基于曲波变换和聚类的地震数据去噪方法和装置，包括：获取第一地震数据，将第一地震数据进行曲波变换，得到第一曲波变换系数；对第一曲波变换系数进行去噪处理，得到第二曲波变换系数；对第二曲波变换系数进行聚类处理，得到第三曲波变换系数；将第三曲波变换系数进行逆曲波变换，得到第二地震数据，可以有效去噪，并保证有效信号的完整性。

技术领域

本发明涉及信号处理技术领域，尤其是涉及基于曲波变换和聚类的地震数据去噪方法和装置。

背景技术

在地震勘探中，野外采集的地震数据不可避免的记录了各种噪音干扰，降低了地震数据的信噪比，使各种叠前处理技术受到一定的影响。为了提高最终地震剖面的信噪比和分辨率，需要对地震数据进行去噪处理。

目前，对地震数据进行去噪处理的方法是通过傅里叶变换和小波变换等方法进行处理，但是，上述处理方法在对地震数据进行去噪处理时，不能保证信号的完整性。

发明内容

有鉴于此，本发明的目的在于提供基于曲波变换和聚类的地震数据去噪方法和装置，可以有效去噪，并保证有效信号的完整性。

第一方面，本发明实施例提供了基于曲波变换和聚类的地震数据去噪方法，所述方法包括：

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

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

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

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

结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所述将所述第一地震数据进行曲波变换，得到第一曲波变换系数包括：

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

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

结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述对所述第一曲波变换系数进行去噪处理，得到第二曲波变换系数包括：

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

结合第一方面的第二种可能的实施方式，本发明实施例提供了第一方面的第三种可能的实施方式，其中，所述去除所述第一曲波变换系数中绝对值小于预设阈值的系数，保留所述第一曲波变换系数中绝对值大于预设阈值的系数，从而得到所述第二曲波变换系数包括：

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

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

结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，所述对所述第二曲波变换系数进行聚类处理，得到第三曲波变换系数包括：

将所述第二曲波变换系数中密度相连的核心点和边界点生成聚类簇；

去除所述第二曲波变换系数中孤立的点；

恢复所述聚类簇中所述核心点邻域内点的数值，从而得到所述第三曲波变换系数。

结合第一方面，本发明实施例提供了第一方面的第五种可能的实施方式，其中，所述将所述第三曲波变换系数进行逆曲波变换，得到第二地震数据包括：