專利名稱：基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法
技術(shù)領(lǐng)域：
本發(fā)明涉及信號消除噪音的方法，特別是一種地震信號的消噪的方法。
二背景技術(shù)：
地震信號被廣泛應(yīng)用于地質(zhì)勘探，地震預(yù)報等應(yīng)用領(lǐng)域。由于地震信號采集環(huán)境 的不同和采集儀器中存在各種干擾，所采集到的地震波形中通常包含了很多噪聲信號，嚴(yán) 重影響了對地震信號的進一步應(yīng)用。要提高地震資料的信噪比，就需要去除地震波形中的 隨機噪聲和想干噪聲。目前常用的用于消除信號噪聲的方法有基于IIR的濾波方法、FIR 濾波方法、基于小波閾值濾波方法等。但用上述方法時需要反復(fù)修改參數(shù)或基函數(shù)，并需要 不斷觀察消噪的效果是否理想。
三
發(fā)明內(nèi)容
本發(fā)明的目的是克服上述不足問題，提供一種地震信號的消噪的方法。該方法采 用自適應(yīng)消噪方法，拋棄了傳統(tǒng)閾值消噪方法中以單個信號點為閾值處理對象的思想，根 據(jù)經(jīng)驗周期模態(tài)單元濾波的原理，該方法不需要先驗基函數(shù)和不斷改變和修改參數(shù)。方法 簡單、可行。 本發(fā)明為實現(xiàn)上述目的所采用的技術(shù)方案是基于經(jīng)驗周期模態(tài)單元濾波的地震 信號消噪方法，首先采集地震波形圖和頻譜圖，利用MATLAB軟件中的imread命令進行讀取 處理，初始預(yù)定信號模擬單元的經(jīng)驗周期范圍，確定兩個經(jīng)驗周期的閾值；然后根據(jù)確定的 閾值和模態(tài)單元的處理系數(shù)確定去除噪聲后的波形形狀，最后利用MATLAB軟件中的plot 命令處理結(jié)果顯示出來。 所述確定閾值和模態(tài)單元的方法為先將信號模態(tài)單元的經(jīng)驗周期范圍初始化為 [thpi—， thpi+]，把imfi中第k個模態(tài)單元記為111,，其中/附乂-{附/， A^7，2,…A), i = 1， 2， . . . ， I。
的經(jīng)驗周期閾值處理方法為 <formula>formula see original document page 3</formula>其中 <formula>formula see original document page 3</formula>
所述所選取的作為初始信號模擬單元的經(jīng)驗周期的范圍為6TS《thpi—《16TS， 25TS《thpi+《45Ts，其中Ts為采樣周期。 本發(fā)明的特點該方法拋棄了傳統(tǒng)閾值消噪方法中以單個信號點為閾值處理對象 的思想，根據(jù)經(jīng)驗周期模態(tài)單元濾波的原理，對地震波形圖進行處理。該方法不需要先驗基 函數(shù)和不斷改變和修改參數(shù)。方法簡單、可行。具體實施例方式
下面結(jié)合具體實施方式
對本發(fā)明作進一步詳細(xì)說明，但不限于具體實施例。 基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法，具體的實施步驟為 1、采集地震波形圖和頻譜圖，利用MATLAB軟件中的imread命令進行讀取處理。 2、初始預(yù)定信號模擬單元的經(jīng)驗周期范圍，確定兩個經(jīng)驗周期的閾值。初始的信
號模擬單元的經(jīng)驗周期范圍的上限和下限分別為25TS《thpi+《45TS，6TS《thpi—16TS，
其中L為采樣周期。 3、先將信號模態(tài)單元的經(jīng)驗周期范圍初始化為[thpi—，thpi+]，把imfi中第k個模 態(tài)單元記為1^，其中!'附/-{附/， &=7，2,...,}， i = 1，2， ， I。 的經(jīng)驗周期閾值處理 方法為m/ (Ap/V/^廣),//2/7廣)w,，閾值處理后的z'm乂-(m, (^pf,^w'+) ， &=J，2W。
權(quán)利要求
基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法，其特征是首先采集地震波形圖和頻譜圖，利用MATLAB軟件中的imread命令進行讀取處理，初始預(yù)定信號模擬單元的經(jīng)驗周期范圍，確定兩個經(jīng)驗周期的閾值；然后根據(jù)確定的閾值和模態(tài)單元的處理系數(shù)確定去除噪聲后的波形形狀，最后利用MATLAB軟件中的plot命令處理結(jié)果顯示出來。
2. 根據(jù)權(quán)利要求1所述的基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法，其特征 是所述確定閾值和模態(tài)單元的方法為先將信號模態(tài)單元的經(jīng)驗周期范圍初始化為[thpi—， thpi+]，把imfi中第k個模態(tài)單 元記為111,，其中/附/-{附/， &=7,2，...人}， i = 1，2， ， I。 的經(jīng)驗周期閾值處理方法閾值處理后的/附乂-{附/ ， ^:=7,2，...^}
3.根據(jù)權(quán)利要求l所述的基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法，其 特征是所選取的作為初始信號模擬單元的經(jīng)驗周期的范圍為6K thpi—《16TS， 25TS《thpi+《45Ts，其中Ts為采樣周期。
全文摘要
本發(fā)明涉及信號消除噪音的方法。基于經(jīng)驗周期模態(tài)單元濾波的地震信號消噪方法，采用自適應(yīng)消噪方法，拋棄于傳統(tǒng)閾值消噪方法中以單個信號點為閾值處理對象的思想，根據(jù)經(jīng)驗周期模態(tài)單元濾波的原理，該方法不需要先驗基函數(shù)和不斷改變和修改參數(shù)。方法簡單、可行。
文檔編號G01V1/36GK101738634SQ20081022860
公開日2010年6月16日 申請日期2008年11月5日 優(yōu)先權(quán)日2008年11月5日
發(fā)明者鄒萍 申請人:鄒萍