S变换时频滤波与其它滤波方法的比较.docVIP

S变换时频滤波与其它滤波方法的比较 第29卷第3期 2007年9月 西北地震 NORTHWESTERNSEISMOLOGICALJOuRNAL Vo1.29No.3 Sept.,2007 S变换时频滤波与其它滤波方法的比较 赵淑红,王璇 (长安大学地质工程及测绘工程学院,陕西西安710054) 摘要:在提高信噪比时通常根据信号和噪音的差异选取不同的滤波方法.由于地震记录的复杂 性,有时有效波与干扰波在视速度或频率上没有明显差异,为此本文提出一种利用S变换时频滤波 的方法,并比较了该方法相对其它滤波方法在时频滤波方面的优势. 关键词:S变换;时频滤波;短时傅立叶变换;F—K滤波 中图分类号:P315.3l文献标识码:A文章编号:i000—0844(2007)03—0224～06 ComparisonbetweentheTime.frequencyFilteringofSTransform MethodandOtherFilterMethod ZHAOShu—hong,WANGXuan (Schoolo/GeologicalEngineeringandSurveyingEngineering.ChannUniversity,Xn710054.China) Abstract:Inusua1differentfiltermethodsareselectedaccordingtothedifferencesbetweeneffec— tivesignalandnoiseforincreasingthesignal—noiseratio.Butforcomplicationofseismicdata, sometimetheeffectivesignalandnoisehavenocleardifferenceinapparentvelocityandfrequen— cy.ForthisproblemamethodusingStransformtofulfillthetime—frequencyfilteringispresen— tedinthispaper,andtheadvantageofStransformintime—frequencyfilteringiscomparedwith otherfiltermethodsthroughexample. Keywords:Stransform;Time-frequencyfiltering;Short—timeFouriertransform:F—Kfiltering 0前言 在提高地震波信噪比方面,人们根据信号与噪 音在某些方面的差异(例如视速度,频率等)来确定 选取什么样滤波方法.常用的去噪方法有:一维滤 波,F—K域滤波,r—P变换滤波,SVD滤波,小波 变换去噪等.但是各种去噪方法都有其使用条件, 例如一维频率滤波可以把与有效信号频带不同的噪 音消除掉;利用F—K域滤波可以把与有效信号视 速度不同的干扰滤掉;SVD滤波假定有效波与干扰 波在空间相干性上有明显差别等.只有地震记录满 足某种去噪方法所需要的条件,才有可能取得良好 的效果.实际地震记录非常复杂,本文提出的利用 S变换时频滤波的方法可用在有效波与干扰波在视 速度或者频率等方面都没有明显的差异时.并且通 过实例比较S变换在时频滤波方面的优势. 1方法原理 1.1S变换 S变换首先由Stockwellc等人提出,它是以 Morlet小波为基本小波的连续小波变换的延伸. 在S变换中基本小波是由简谐波与高斯函数的乘积 构成的.与连续小波变换不同,基本小波中的简谐 波在时间域仅作伸缩变换,而高斯函数则进行伸缩 和平移.而在连续小波变换中,简谐波与高斯函数 进行同样的伸缩和平移. S变换可以认为是连续小波变换(CWT)的相 位校正.函数^()的S变换表示为 厂∞ S(t,厂)一I^(r)?(r一)?exp(--i27c)dr J—.∞ (1) 这里 收稿日期:2006—10—07 作者简介:赵淑红(1973一),女(汉族),陕西大荔人,博士,讲师,目前从事地球物理的教学与研究工作 第3期赵淑红等:S变换时频滤波与其它滤波方法的比较225 ,一exp(一);口~/27【,z(2) 厂南 sj,?唧(一) exp(一i27【/)dr (3) 用w()表示训()的傅立叶变换 wc,一exp(一等.)㈥ 的单位和_厂的单位相同. s变换的窗函数必须满足条件:1w(r—t, 厂)dr一1,因此,S变换和傅立叶变换满足关系式 Is(t,厂)dr=Ih()exp(--27【i)×I训(r—t,厂)drdt—∞.. 一 Ih()exp(--2~d)dt—H(.厂) (5) 同时,式(5)可以确保S变换是可逆的. 根据短时傅立叶变换在时间域和频