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