基于vmd的地震信号随机噪声压制方法技术

本发明专利技术公开了一种基于vmd的地震信号随机噪声压制方法，该方法包括以下步骤：选择一个时间窗口将原始含噪二维地震信号作傅里叶变换，变换至f‑x域；对每一个频率切片数据进行模态分解；将分解得到的BIMF分量组合得到滤波后的信号；将信号作傅里叶逆变换，变换回t‑x域；重复对下一个时间窗口进行同上操作，地震记录全部处理完毕后，即得到最终的二维地震去噪结果。同时，本发明专利技术中将一维变分模态分解拓展至二维复数形式，可用于三维地震去噪处理。基于VMD的地震信号随机噪声压制方法具有优异的噪声压制、幅值保持性能的同时，还具备较高的计算效率，可满足高维大尺度地震数据的处理要求。

Random noise suppression method of seismic signal based on VMD

The invention discloses a method for suppressing random noise of seismic signals based on vmd. The method comprises the following steps: selecting a time window to transform the original two-dimensional noisy seismic signals into the f_x domain by Fourier transform; modal decomposition of each frequency slice data; combining the decomposed BIMF components to filter. The signal after wave is transformed back to t_x domain by inverse Fourier transform, and the next time window is repeatedly operated as above. After all the seismic records are processed, the final two-dimensional seismic denoising results are obtained. Meanwhile, in the present invention, one-dimensional variational mode decomposition is extended to two-dimensional complex form, which can be used for three-dimensional seismic denoising processing. VMD-based seismic signal random noise suppression method has excellent noise suppression, amplitude preservation performance, but also has high computational efficiency, can meet the requirements of high-dimensional large-scale seismic data processing.

【技术实现步骤摘要】
基于vmd的地震信号随机噪声压制方法
本专利技术涉及地球物理勘探的地震资料处理
，尤其涉及地震信号中随机噪声压制和有效能量保持处理技术。技术背景地震信号随机噪声衰减是地震数据处理领域的热点和难点问题。从观测数据中有效去除随机噪声干扰，提高信号信噪比和分辨率，是正反演计算和地质解释的前提。目前，地震信号随机噪声压制方法主要有：基于滤波理论的方法、基于小波域变换方法、基于矩阵理论的变换方法和基于信号分解理论方法等。近年来，基于信号分解的经验模态分解(EmpiricalModeDecomposition，EMD)方法，由于在处理非平稳及非线性数据上具有很高的信噪比，已成为信号分解领域研究的热点。完备集合经验模态分解(CompleteEnsembleEMD，CEEMD)避免了EMD方法中模态混叠的问题，有效改善集合经验模态分解EEMD方法中加入随机白噪声对原始信号造成污染的缺点，但CEEMD方法中依然采用的是递归迭代式筛选分解过程，非平稳地震信号极值点插值和包络计算过程依然耗时过长，在处理多维度和多尺度地震数据时存在一定限制。为解决这些问题，进一步提升信号分解的精度，本专利技术提出基于变分模态分解(VariationalModalDecomposition，VMD)的地震信号随机噪声压制处理方法。通过将模态分解过程转变为变分模型的最优解迭代求解过程，以每个分量的估计聚集带宽之和最小为约束，双重对偶上升使得每个分量的频率中心及带宽在频域内进行自适应分离。与EMD的递归“筛选”模式不同，VMD将信号分解转化非递归、变分模型泛函极值求解问题，由变分模型优化目标函数求解过程可知，其实质是多个维纳滤波组的自适应拓展，具有优良的噪声鲁棒性，对于非平稳的时序信号分解更为适用。
技术实现思路
本专利技术的目的在于提供一种基于VMD的地震信号随机噪声压制方法，以解决地震勘探任务中采集的地震数据被随机噪声干扰的问题，该方法具有优异的噪声压制、幅值保持性能的同时，还具备较高的计算效率，可满足高维大尺度地震数据的处理要求。本专利技术的技术方案为：一种基于vmd的地震信号随机噪声压制方法，其技术方案步骤包括：S1：地震数据的时域信号转换为频域信号：选择一个时间窗口将原始含噪二维地震信号d(x,t)作傅里叶变换，变换至f-x域；S2：信号的频段分解：对每一个频段数据进行变分模态分解VMD(VariationalModeDecomposition)，生成带限固有模态函数(Band-LimitedIntrinsicModeFunction)分量，简称BIMF分量；S3：信号的组合：将VMD分解得到的BIMF分量进行组合，得到新的滤波信号；S4：地震数据的频域信号转换为时域信号：将S3步骤得的滤波信号作傅里叶逆变换，变换回t-x域；S5：迭代过程：重复对下一个时间窗口进行同上操作，地震记录全部处理完毕后，即得到最终的二维地震去噪结果。所述S2步骤的VMD分解过程可转换至变分泛函最优求解过程，其分解包括：通过在频率域中估计固有模态函数BIMF(Band-LimitedIntrinsicModeFunction)分量频率带宽目标函数，其数学表达式如下：其中，K为预设分解尺度个数，t是时间变量，uk即为VMD分解后具备带限性质的BIMF分量，ωk为对应模态的频率中心，是关于时间t的导数，δ(t)为Dirac(狄利克雷)冲击函数，*为卷积符号，f为原始的频率域实值信号，是平方L2范数。的意义是通过Hilbert变换，将每个模态函数uk变为解析信号，使实值信号uk转变为复值，以获得uk的单边频谱。式(1)中的变分问题使得每个BIMF分量的频谱带宽在其中心频率ωk附近，且要求BIMF的带宽具备稀疏性。通过将VMD方法的变分优化框架拓展至复数空间，因此专利要求2中的等式改写为如下形式：采用增广Lagrange函数来计算其最优解。通过引入二次惩罚因子α和拉格朗日乘法算子λ，将约束性变分问题转换为如下的非约束性变分形式：其中，二次惩罚因子α是控制数据保真度的均衡参数，用于平衡变分正则项和二次约束项，在含噪声情形时可保证信号重构精度；Lagrange乘子λ(x)可以保证模型约束条件的严格性；然后，通过采用交替方向乘子算法ADMM(AlternateDirectionMethodofMultipliers)求解式，具体步骤如下：MM_1：初始化λ1，n＝0；MM_2：n＝n+1，执行主循环；MM_3：Fork＝1:K-1，执行第一个内循环更新uk：MM_4：k＝K，结束第一个内循环；:MM_5：Fork＝1:K-1，执行第二个内循环更新ωk：:MM_6：k＝K，结束第二个内循环；MM_7：对于所有ωk＞0，双重对偶上升，更新λ：其中，τ表示噪声容限参数；在噪声压制任务中(而不是信号的重构)，将更新参数τ＝0，以得到更好的去噪效果。MM_8：给定判定精度ε＞0，重复步骤2)-7)，直至满足迭代停止条件：结束迭代，即得到K个带限BIMF分量。ADMM的求解过程包含VMD的模态更新和频率中心更新。其中，ωk频率中心的更新由对应模态的能量谱重心得到，uk模态更新对应于1/αω2的Wiener滤波器结构，α为白噪声方差，1/ω2表示信号的能量谱为低通形式。参数α控制着Wiener滤波器的宽度，本专利技术称为保真度均衡参数。增大α值，Wiener滤波器宽度变窄，可以滤除更多的噪声，但也使得其包含更少的真实峰信息，同时算法趋于发散不收敛的几率增加，反之亦然。BIMF分量的个数K值的取值是VMD算法的关键问题，不同的模态个数会对分解结果产生影响，从而影响最终解的评价。本专利技术通过计算BIMF分量瞬时频率均值变化的极大值确定最佳的模态分解个数。本专利技术的优点：VMD方法的分解过程可转换至变分泛函最优求解过程，与EMD类方法的迭代筛选模式不同，信号分解过程被转移到变分框架，以每个带限BIMF分量的估计带宽之和最小为约束，通过增广Lagrange目标函数将变分问题由约束性变为非约束性，采用交替方向乘子ADMM算法寻求变分泛函的最优解达到信号自适应分解的目的。ADMM中频率中心及带宽交替更新对偶上升，使得两者同时达到最优趋势，一次性得到所有BIMF分量，具有更高的时间效率，每个模态分量在频谱上均具有带限特性，实现信号频带的自适应剖分。BIMF分量的更新过程中具备维纳滤波特性，因而VMD可看作是维纳滤波的多重化和自适应阶推广，即分解出的各组分能量已进行了类似维纳滤波操作，该过程具有坚实的理论基础。一维VMD可拓展至频率域二维复数的广义形式。基于VMD的地震信号随机噪声压制处理方法具有优异的噪声压制、幅值保持性能的同时，还具备较高的计算效率，可满足高维大尺度地震数据的处理要求。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作进一步介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。附图1是一种基于VMD的地震信号随机噪声压制处理方法实现流程图；附图2是各组分信号及VMD分解BIMF分量；附图3是附图2合成信号的频谱分布与VMD分解后本文档来自技高网...

【技术保护点】
1.一种基于VMD的地震信号随机噪声压制方法，其步骤包括：S1：选择一个时间窗口将原始含噪二维地震信号d(x,t)作傅里叶变换，变换至f‑x域；S2：对每一个频段数据进行变分模态分解，生成带限固有模态函数(BIMF)分量；S3：将VMD分解得到的BIMF分量进行组合，生成滤波后的信号；S4：将信号作傅里叶逆变换，变换回t‑x域；S5：重复对下一个时间窗口进行同上操作，地震记录全部处理完毕后，即得到最终的二维地震去噪结果。

【技术特征摘要】
1.一种基于VMD的地震信号随机噪声压制方法，其步骤包括：S1：选择一个时间窗口将原始含噪二维地震信号d(x,t)作傅里叶变换，变换至f-x域；S2：对每一个频段数据进行变分模态分解，生成带限固有模态函数(BIMF)分量；S3：将VMD分解得到的BIMF分量进行组合，生成滤波后的信号；S4：将信号作傅里叶逆变换，变换回t-x域；S5：重复对下一个时间窗口进行同上操作，地震记录全部处理完毕后，即得到最终的二维地震去噪结果。2.根据权利要求1所述的基于VMD的地震信号随机噪声压制方法，其特征在于：所述S2步骤的VMD分解过程包括：首先通过在频率域中估计固有模态函数BIMF分量频率带宽目标函数，其数学表达式如下：其中，K为预设分解尺度个数，t是时间变量，uk即为VMD分解后具备带限性质的BIMF分量，ωk为对应模态的频率中心，是关于时间t的导数，δ(t)为Dirac(狄利克雷)冲击函数，*为卷积符号，f为原始的频率域实值信号，是平方L2范数；的意义是通过Hilbert变换，将每个模态函数uk变为解析信号，使实值信号uk转变为复值，以获得uk的单边频谱；式(1)中的变分问题使得每个BIMF分量的频谱带宽在其中心频率ωk附近，且要求BIMF的带宽具备稀疏性；然后，通过将VMD方法的变分优化框架拓展至复数空间，因此专利要求2中的等式改写为如下形式：采用增广Lagrange函数来计算其最优解；通过引入二次惩罚因子α和拉格朗日乘法算子λ，将约束性变分问题转换为如下的非约束性变分形式：其中，二次惩罚因子α是控制数据保真度的均衡参数，用于平衡变分正则项和二次约束项，在含噪声情形时可保证信号重构精度；Lagrange乘子λ(x)可以保证模型约束条件的严格性；最后，通过采用交替方向乘子算法ADMM求解式，具体步骤如下：MM_1：初始化n＝0；MM_2：n＝n+1，执行主循环；MM_3：Fork＝1:K-1，执行第一个内循环更新uk：MM_4：k＝K，结束第一个内循环；:MM_5...

【专利技术属性】
技术研发人员：方江雄，温志平，顾华奇，刘军，刘雄涛，
申请(专利权)人：东华理工大学，
类型：发明
国别省市：江西,36