随机噪声的压制c语言,一种地震数据随机噪声压制方法及装置的制造方法

一种地震数据随机噪声压制方法及装置的制造方法

【技术领域】

[0001] 本发明涉及地震勘探领域，尤其涉及一种地震数据随机噪声压制方法及装置。

【背景技术】

[0002] 按噪声在地震剖面上出现的特征，将噪声分为规则噪声和不规则噪声。规则噪声 主要是指有一定主频和一定视速度的噪声，如面波、交流电干扰、声波、浅层折射等。不规则 噪声即随机噪声，主要指没有固定频率和固定传播方向的波，在地震数据中形成杂乱无章 的背景。随机噪声在地震记录上表现为频率很宽，无一定视速度，因而很难利用随机噪声同 有效波之间在频谱上的差异或传播方向上的差异对其进行压制。由于地表条件的复杂性， 地震资料中常含有许多随机噪声，如微震、背景干扰等。这些噪声分布很广，严重影响了地 震资料的信噪比。

[0003] 地震勘探领域中，地震数据随机随机噪声压制方法的研究一直受到相关研究人员 的广泛关注。目前地震数据随机噪声压制方法有傅氏变换滤波法、独立成分分析法和二维 小波变换法等。傅氏变换滤波法需要假定信号是平稳的，并且难以表征信号的局部特点，当 信号非平稳时，该方法难以达到理想的效果;独立成分分析法要求地震信号和随机噪声在 统计上相互独立，但该方法假设数据道是对相同道的多次观测，而实际地震道数据并不满 足这种假设，指示一种很粗略的近似，因此其对随机噪声的压制能力并不理想。二维小波变 换法在地震数据处理中有广泛应用，但由于其基函数是各项同性的，在处理地震数据时变 换系数的局部极大值只能反映出这个小波系数出现的位置是"过"边缘的，而无法表达"沿" 边缘的信息，因此二维小波变换法对图像边缘数据处理并不理想。

[0004] 随着技术的发展，勘探目标逐渐向深处和地表复杂地区转移，地震数据处理要求 也越来越高，因此，如何有效的压制噪声，从而提高地震资料的信噪比是地震勘探领域中亟 需解决的问题。

【发明内容】

[0005] 本申请实施例提供了一种地震数据随机噪声压制方法及装置，以保持有效信号， 压制噪声，从而提高地震资料的信噪比。

[0006] 为达到上述目的，一方面，本申请实施例提供了一种地震数据随机噪声压制方法， 包括以下步骤：

[0007] 对获取的地震数据进行互补经验态分解CEEMD，获得本征模函数頂F分量序列；

[0008] 确定所述IMF分量序列中含有随机噪声的IMF分量和不含有随机噪声的IMF分量；

[0009] 对所述含有随机噪声的MF分量进行曲波最优迭代阈值去噪处理，获得处理后的 IMF分量；

[0010] 利用所述处理后的IMF分量以及所述不含有随机噪声的IMF分量进行重构获得去 除随机噪声的地震数据。

[0011] 本申请实施例的地震数据随机噪声压制方法，所述对含有随机噪声的IMF分量进 行曲波最优迭代阈值去噪处理，获得处理后的頂F分量，包括：

[0012] 根据所述含有随机噪声的IMF分量中所含随机噪声的不同，对所述含有随机噪声 的頂F分量选取不同的阈值进行曲波最优迭代阈值去噪处理，获得处理后的頂F分量。

[0013] 本申请实施例的地震数据随机噪声压制方法，所述确定IMF分量序列中含有噪声 的頂F分量和不含有噪声的頂F分量，包括以下步骤：

[0014] 计算所述地震数据和所述頂F分量序列之间的互相关系数，获得互相关系数序列；

[0015] 确定所述互相关系数序列中的极小值；

[0016] 根据所述极小值对应的IMF分量位置确定含有随机噪声的IMF分量和不含有随机 噪声的MF分量。

[0017] 本申请实施例的地震数据随机噪声压制方法，所述计算地震数据和所述MF分量 之间的互相关系数，具体公式为：

[0019] 其中，R(j)为互相关系数，Rdata(0)为地震数据自相关后零时刻所对应的最大峰值 点(仍为第j个IMF分离的自相关后零时刻所对应的最大峰值点，为原始数 据与第j个頂F分量互相关后零时刻所对应的值。

[0020] 本申请实施例的地震数据随机噪声压制方法，所述利用处理后的MF分量以及所 述不含有随机噪声的MF分量进行重构获得去除随机噪声的地震数据的计算公式为：

[0022] C~m为处理后的頂F分量;Cn为不含有随机噪声的MF分量，其中η为CEEMD分解得到 的IMF分量的个数;Χ~为去除噪声后的地震数据;k为极小值对应的IMF分量位置。

[0023]另一方面，本申请实施例还提供一种地震数据随机噪声压制装置，包括：

[0024] CEEMD分解单元，用于对获取的地震数据进行互补经验态分解CEEMD，获得本征模 函数IMF分量序列；

[0025] IMF分量分析单元，用于确定所述IMF分量序列中含有随机噪声的IMF分量和不含 有随机噪声的IMF分量；

[0026] 去噪处理单元，用于对所述含有随机噪声的IMF分量进行曲波最优迭代阈值去噪 处理，获得处理后的MF分量；

[0027]重构单元，用于利用所述处理后的頂F分量以及所述不含有随机噪声的頂F分量进 行重构获得去除随机噪声的地震数据。

[0028] 本申请实施例的地震数据随机噪声压制装置，所述对含有随机噪声的IMF分量进 行曲波最优迭代阈值去噪处理，获得处理后的頂F分量，包括：

[0029] 根据所述含有随机噪声的IMF分量中所含随机噪声的不同，对所述含有随机噪声 的頂F分量选取不同的阈值进行曲波最优迭代阈值去噪处理，获得处理后的頂F分量。

[0030] 本申请实施例的地震数据随机噪声压制装置，所述頂F分量分析单元包括：

[0031] 互相关系数计算子单元，用于计算所述地震数据和所述MF分量序列之间的互相 关系数，获得互相关系数序列；

[0032]极小值确定子单元，用于确定所述互相关系数序列中的极小值；

[0033]噪声分析子单元，用于根据所述极小值对应的IMF分量位置确定含有随机噪声的 IMF分量和不含有随机噪声的IMF分量。

[0034]本申请实施例的地震数据随机噪声压制装置，所述计算地震数据和所述IMF分量 之间的互相关系数，具体公式为：

[0036]其中，R(j)为互相关系数，Rdata(0)为地震数据自相关后零时刻所对应的最大峰值 点，Aa/f, (〇)为第j个IMF分离的自相关后零时刻所对应的最大峰值点，爲_￡^ (〇)为原始数 据与第j个頂F分量互相关后零时刻所对应的值。

[0037]本申请实施例的地震数据随机噪声压制装置，所述利用处理后的MF分量以及所 述不含有随机噪声的MF分量进行重构获得去除随机噪声的地震数据的计算公式为：

[0039] C~m为处理后的頂F分量;Cn为不含有随机噪声的MF分量，其中η为CEEMD分解得到 的IMF分量的个数;Χ~为去除噪声后的地震数据;k为极小值对应的IMF分量位置。

[0040] 本申请实施例中，通过对地震数据进行CEEMD分解获得MF分量序列，然后通过地 震数据和頂F分量之间的互相关系数确定含有随机噪声的頂F分量和不含有随机噪声的頂F 分量;再根据含有随机噪声的IMF分量中所含随机噪声的不同，对所述含有随机噪声的IMF 分量选取不同的