基于高通滤波的频率-空间域经验模态分解压制高频噪声
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
高频噪声压制是高分辨率地震数据处理中提高信噪比的关键性问题.本文针对f-x(频率-空间)反褶积空间预测滤波器无法处理非平稳、非线性信号的缺点,提出了一种基于高通滤波的频率-空间域经验模态分解(Empirical Mode Decomposition in the frequency-space domain,f-xEMD)压制地震剖面中高频噪声的方法.该方法采用全域高通滤波从原始数据中分离出含有部分有效信号的高频数据,将其变换到f-x域,然后在滑动的短窗口内提取每一个频率的空变数据序列进行EMD分解得到高频复本征模态函数(Intrinsic Mode Function,IMF)IMF1,将所有频率的IMF1序列反Fourier变换到时间域得到噪声剖面,将其与原始数据相减,达到高频噪声压制的目的.该方法可克服传统EMD分解方法中的模态混叠现象,保护陡倾角反射同相轴;压制后的噪声剖面中不包含有效信号能量,地震剖面的信噪比得到了提高.模拟数据和实际数据处理结果充分证明了该方法的有效性.
It is a key issue that high frequency noise is suppressed in high resolution seismic data processing.In order to overcome some shortcomings that f-x(frequency-spatial)deconvolution spatial prediction filter can't cope with the nonstationary and nonlinear data,this paper presents a method that high frequency noise in the seismic profile is suppressed by Empirical Mode Decomposition in the frequency-spatial domain(f-xEMD) based on high-pass filter.This method can separate the high frequency data in which include some effective signal from the raw seismic data by the whole region high-pass filter and transform it into frequency-spatial domain.For every frequency,the complex high frequency intrinsic mode function(IMF1) is computed from the spatial-varying data in the sliding time window by f-xEMD,then inverse transform all the complex IMF1 back to the time-spatial domain for gaining the noise section.Signal-to-noise enhancement can be achieved by subtracting noise section from raw data.This method can overcome the shortcoming of mode mixing in the conventional EMD method and protect the steeply dipping reflect event.The noise section doesn't actually contain the effective signal energy so that this method can restrain the high frequency noise and enhance the signal-to-noise of seismic profile.The processing result of simulated data and actual data indicate that the high frequency noise suppression method by f-xEMD based on high-pass filter is very effective for high resolution seismic data processing.
引文
[1]W ang Y H.A stable and efficient approach of inverse Qfiltering.Geophysics,2002,67(2):657-663.
    [2]Hargreaves,N D,Calvert,A J.Inverse Q filtering by Fouriertransform.Geophysics,1991,56(4):519-527.
    [3]Zhang J F,Wapenaar K.Wavefield extrapolation and prestackdepth migration in anelastic inhomogeneous media.GeophysicalProspecting,2002,50(6):629-643.
    [4]Canales L.Random noise reduction.54th Annual InternationalMeeting,SEG,Expanded Abstracts,1984:525-527.
    [5]Galbraith M.Random noise attenuation by f-x prediction:Atutorial.61st SEG,Annual International Meeting,ExpandedAbstracts,1991:1428-1431.
    [6]Harris P E,White R E.Improving the performance of f-xprediction filtering at low signal-to-noise ratios.GeophysicalProspecting,1997,45(2):269-302.
    [7]陈凯.基于经验模式分解的去噪方法.石油地球物理勘探,2009,44(5):603-608.Chen K.A new denoising method based on Empirical ModeDecomposition(EMD).Oil Geophysical Prospecting(inChinese),2009,44(5):603-608.
    [8]Bekara M,van der Baan M.Random and coherent noiseattenuation by empirical mode decomposition.Geophysics,2009,74(5):V89-V98.
    [9]Huang N E,Shen Z,Long S R,et al.The empirical modedecomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis.Proceedings of the RoyalSociety A:Mathematical,Physical and Engineering Sciences,1998,454(1971):903-995.
    [10]Huang N E,Zheng S,Long S R.A new view of nonlinearwater waves:the Hilbert spectrum.Annu.Rev.Fluid Mech.,1999,31:417-457.
    [11]Huang N E,Wu M L C,Long S R,et al.A confidence limitfor the empirical mode decomposition and Hilbert spectralanalysis.Proc.R.Soc.Lond,SERIESA,2003,459(2073):2317-2345.
    [12]Wu Z,Huang N E.A study of the characteristics of whitenoise using the empirical mode decomposition method.Proc.R.Soc.Lond,SERIESA,2004,460(2046):1587-1611.
    [13]Flandrin P,Rilling G,Goncalves P.Empirical mode decompositionas a filter bank.IEEE Signal Processing Letters,2004,11(2):112-114.
    [14]Rilling G,Flandrin P.Sampling effects on the empirical modedecomposition.Advanced in Adaptive Data Analysis,2009,1(1):43-59.
    [15]Battista B M,Knapp C,McGee T,et al.Application of theempirical mode decompositionand Hilbert-Huang transformto seismic reflection data.Geophysics,2007,72(2):H29-H37.
    [16]Nunes J C,Guyot S,Deléchelle E.Texture analysis based onlocal analysis of the bidimensional empirical modedecomposition.Machine Vision and Applications,2005,16(3):177-188.
    [17]Linderhed A.2Dempirical mode decompositions in the spiritof image compression.Society of Photo-Optical InstrumentationEngineers(SPIE)Conference Series,2002,4738(1):1-8.
    [18]邓拥军,王伟,钱成春等.EMD方法及Hilbert变换中边界问题的处理.科学通报,2001,46(3):257-263.Deng Y J,Wang W,Qian C C,et al.Comment andmodification on EMD and Hilbert transform method.ChineseScience Bulletin(in Chinese),2001,46(3):257-263.
    [19]Magrin-Chagnolleau I,Baraniuk R G.Empirical modedecomposition based time-frequency attributes.69thAnnualInternational Meeting,SEG,Expanded Abstracts,1999:1949-1953.
    [20]Hassan H,Peirce J W.Empirical mode decomposition(EMD)of potential field data:Airborne gravity data as an example.75th Annual International Meeting,SEG,ExpandedAbstracts,2005,24:704-707.
    [21]Sinclair S,Pegram G G S.Empirical mode decomposition in 2-D space and time:a tool for space-time rainfall analysis andnow casting.Hydrology and Earth System Sciences,2005,9(3):127-137.
    [22]胡劲松,杨世锡,吴昭同等.基于经验模态分解的旋转机械振动信号滤波技术研究.振动、测试与诊断,2003,23(2):96-98.Hu J S,Yang S X,Wu Z T.The study of EMD-based filteringtechnique for vibration signals in rotating machinery.Journalof Vibration Measurement&Diagnosis(in Chinese),2003,23(2):96-98.
    [23]岳焕印,郭华东,韩春明等.经验模态分解技术在SAR干涉图滤波中的应用.高技术通讯,2001,11(12):37-40.Yue H Y,Guo H D,Hang C M,et al.A SAR interferogramfilter based on the Empirical Mode Decomposition method.High Technology Letters(in Chinese),2001,11(12):37-40.
    [24]黄朝军,徐彤,吴健等.利用EMD方法提取地磁Ap指数周期分量.地球物理学进展,2009,24(6):1968-1974.Huang C J,Xu T,Wu J,et al.Extracting periodic componentsof geomagnetic Apindex using EMD method.Progress inGeophys.(in Chinese),2009,24(6):1968-1974.
    [25]宋海斌,拜阳,董崇志等.南海东北部内波特征——经验模态分解方法应用初探.地球物理学报,2010,53(2):393-400.Song H B,Bai Y,Dong C Z,et al.A preliminary study ofapplication of Empirical Mode Decomposition method inunderstanding the features of internal waves in thenortheastern South China Sea.Chinese J.Geophys.(inChinese),2010,53(2):393-400.
    [26]杨培杰,印兴耀,张广智.希尔伯特-黄变换地震信号时频分析与属性提取.地球物理学进展,2007,22(5):1585-1990.Yang P J,Yin X Y,Zhang G Z.Seismic signal time-frequencyanalysis and attributes extraction based on HHT.Progress inGeophys.(in Chinese),2007,22(5):1585-1990.
    [27]陈林,宋海滨.基于经验模态分解的地震瞬时属性提取.地球物理学进展,2008,23(4):1179-1185.Chen L,Song H B.Seismic instantaneous attribute extractionbased on empirical mode decomposition.Progress in Geophys.(in Chinese),2008,23(4):1179-1185.
    [28]钱昌松,刘代志,刘志刚等.基于递归高通滤波的经验模态分解及其在地震信号分析中的应用.地球物理学报,2010,53(5):1215-1225.Qian C S,Liu D Z,Liu Z G,et al.EMD based on recursivehigh-pass filter and its application on seismic signal analysis.Chinese J.Geophys.(in Chinese),2010,53(5):1215-1225.
    [29]宗涛,孟鸿鹰,贾玉兰等.小波包F-X域前后向预测去噪.地球物理学报,1998,41(S1):337-346.Zong T,Meng H Y,Jia Y L,et al.Denoising method based onwavelet packet and forward-backward prediction in F-Xdomain.Chinese J.Geophys.(in Chinese),1998,41(S1):337-346.
    [30]章珂,刘贵忠,邹大文等.二进小波变换方法的地震信号分时分频去噪处理.地球物理学报,1996,39(2):265-271.Zhang K,Liu G Z,Zou D W,et al.Seismic data time-frequency domain denoising based on dyadic wavelettransform.Chinese J.Geophys.(in Chinese),1996,39(2):265-271.
    [31]刘红伟,刘洪,邹振等.地震叠前逆时偏移中的去噪与存储.地球物理学报,2010,53(9):2171-2180.Liu H W,Liu H,Zou Z,et al.The problems of denoise andstorage in seismic reverse time migration.Chinese J.Geophys.(in Chinese),2010,53(9):2171-2180.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心