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将信号分解转化非递归、变分模型泛函极值求解问题,由变分模型优化目标 ...
【技术保护点】
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
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。