徐宏斌 李庶林， 陳際經(jīng)

1）中國福建廈門361005廈門大學(xué)建筑與土木工程學(xué)院

2）中國湖南郴州423037湖南柿竹園有色金屬有限責(zé)任公司

基于小波變換的大尺度巖體結(jié)構(gòu)微震監(jiān)測信號(hào)去噪方法研究＊

徐宏斌1）李庶林1），陳際經(jīng)2）

1）中國福建廈門361005廈門大學(xué)建筑與土木工程學(xué)院

2）中國湖南郴州423037湖南柿竹園有色金屬有限責(zé)任公司

為將小波去噪方法應(yīng)用于大尺度巖體結(jié)構(gòu)微震監(jiān)測信號(hào)的去噪研究，首先在MATLAB環(huán)境下進(jìn)行仿真，驗(yàn)證了使用Symlet6小波進(jìn)行小波去噪的可行性；利用4種自適應(yīng)閾值規(guī)則對含噪信號(hào)進(jìn)行去噪對比，結(jié)果表明4種閾值去噪后的信號(hào)在均方差較小的情況下都極大地提高了信號(hào)的信噪比，有效地去除了噪聲，對不同的含噪信號(hào)，無偏似然原則閾值去噪后的信噪比最高，同時(shí)均方差也最小，在去噪時(shí)顯得更為有效；以柿竹園全數(shù)字多通道微震監(jiān)測系統(tǒng)為背景，將MATLAB仿真結(jié)論應(yīng)用于現(xiàn)場微震信號(hào)的去噪研究，結(jié)果表明小波閾值去噪特別適合大尺度巖體微震信號(hào)這一類非穩(wěn)定信號(hào)的去噪分析，既可以對低信噪比的微震信號(hào)提取出有效信號(hào)，也可以對頻率覆蓋范圍廣的微震信號(hào)在各尺度上提取并重構(gòu)出有效信號(hào)，實(shí)現(xiàn)了對微震真實(shí)信號(hào)和噪聲信號(hào)的有效分離.

小波去噪 大尺度巖體 微震技術(shù) MATLAB仿真

引言

微震監(jiān)測技術(shù)作為一種先進(jìn)的實(shí)時(shí)監(jiān)測技術(shù)，在地下工程安全監(jiān)測和防災(zāi)減災(zāi)監(jiān)測方面具有較好的應(yīng)用價(jià)值和前景（李庶林，2009）.自2003年本文作者在凡口鉛鋅礦建立起了我國礦山地壓監(jiān)測的第一套全數(shù)字型64通道微震監(jiān)測系統(tǒng)以來（李庶林等，2005），全數(shù)字多通道微震監(jiān)測系統(tǒng)開始在國內(nèi)得到了較多的應(yīng)用，尤其是在礦山等大尺度巖體結(jié)構(gòu)的穩(wěn)定性和地壓監(jiān)測方面，為礦山的安全生產(chǎn)提供了有力的技術(shù)支撐.微震監(jiān)測技術(shù)是一門典型的理論落后于實(shí)踐的技術(shù)，該技術(shù)在應(yīng)用領(lǐng)域還有許多亟待解決的應(yīng)用理論和技術(shù)問題.通常，微震系統(tǒng)的工作環(huán)境復(fù)雜，采集的信號(hào)包含著各式各樣的噪聲，因此對微震監(jiān)測系統(tǒng)采集到的信號(hào)進(jìn)行去噪和提純就是其中一個(gè)重要問題，它是準(zhǔn)確計(jì)算信號(hào)特征參量以及精確定位微震事件的前提.

傳統(tǒng)的信號(hào)處理方法包括低通、高通和帶通濾波，這些方法對信號(hào)的處理是比較粗略的，在噪音環(huán)境比較簡單的情況下是適用的.但是對于噪聲環(huán)境復(fù)雜的微震信號(hào)而言，傳統(tǒng)方法就顯現(xiàn)出了局限性，既不能有效地去除各頻段的噪聲，也不能在去除噪聲的同時(shí)有效地保留信號(hào)的細(xì)節(jié)信息.小波變換在時(shí)域和頻域都有很好的局部化性質(zhì)，具有多分辨率分析的特點(diǎn)，在非穩(wěn)定信號(hào)的去噪分析方面具有重要的應(yīng)用.Mallat和Hwahg（1992）提出的利用小波變換模極大值原理進(jìn)行信號(hào)去噪的方法是小波去噪中最經(jīng)典的方法.Xu等（1994）提出了一種基于空域相關(guān)性的噪聲去除方法.Donoho和Johnstone（1994）提出了信號(hào)去噪的軟閾值方法和硬閾值方法.同年，Coifman和Donoho（1994）提出了平移不變小波去噪.Johnstone和Silverman（1997）給出一種相關(guān)噪聲去除的小波閾值估計(jì)器.Chang等（2000）提出一種針對圖像的空域自適應(yīng)小波閾值去噪方法.Zhang和Salari（2005）提出了基于神經(jīng)網(wǎng)絡(luò)的圖像去噪算法.近年來，基于小波的去噪方法在國內(nèi)信號(hào)處理和圖像處理中得到了較為廣泛地研究（潘泉等，2005）.已經(jīng)開始有部分學(xué)者將小波去噪技術(shù)應(yīng)用于巖石類材料的微震、聲發(fā)射研究中.趙奎等（2006）通過確定室內(nèi)巖石聲發(fā)射實(shí)驗(yàn)Kaiser點(diǎn)信號(hào)的頻率范圍，提出了基于小波變換的Kaiser點(diǎn)信號(hào)的信噪分離方法.王更峰（2006）對室內(nèi)巖石聲發(fā)射信號(hào)進(jìn)行小波消噪，并提出了基于能量閾值法的小波包消噪.金解放等（2007）總結(jié)了適合于室內(nèi)巖石聲發(fā)射信號(hào)分析和處理的小波基應(yīng)具有的特點(diǎn)，并在理論上歸納出Daubechies小波族中的Db4小波基能滿足巖石聲發(fā)射信號(hào)處理的要求.許大為等（2008）對小波去噪方法在礦山微震信號(hào)的去噪研究中進(jìn)行了初步的應(yīng)用研究，指出小波去噪方法可應(yīng)用于微震信號(hào)的預(yù)處理.

目前，將小波去噪方法應(yīng)用于巖石類材料的微震、聲發(fā)射技術(shù)的研究較少，研究范圍也多限于室內(nèi)實(shí)驗(yàn)，尤其是基于全數(shù)字多通道微震監(jiān)測系統(tǒng)的現(xiàn)場實(shí)際應(yīng)用研究尚未看到.鑒于此，本文依托柿竹園礦微震監(jiān)測系統(tǒng)（袁節(jié)平等，2010），在仿真實(shí)驗(yàn)的基礎(chǔ)上，將小波去噪技術(shù)應(yīng)用于現(xiàn)場微震監(jiān)測信號(hào)的去噪研究.

1 小波變換基本理論

1.1 連續(xù)小波變換

連續(xù)小波變換（continuous wavelet transform，簡寫為CWT）（潘泉等，2005）定義為

其中函數(shù)系

1.2 離散小波變換

對連續(xù)小波變換的定義式中，將其尺度參數(shù)a和平移參數(shù)b進(jìn)行如下的離散化采樣：

則對應(yīng)的離散小波為

離散小波變換的定義為

1.3 Mallat算法

在應(yīng)用這些Mallat快速算法做時(shí)頻分析時(shí)，只需要分析信號(hào)的有關(guān)數(shù)據(jù)和二尺度方程的傳遞系數(shù)｛h（k）｝和｛g（k）｝，并不涉及尺度函數(shù)φ（t）的具體表達(dá)式.

分解算法要實(shí)現(xiàn)的目標(biāo)是：在｛φ（t－k）｝是標(biāo)準(zhǔn)正交基條件下，已知｛h（k）｝和｛g（k）｝，求出｛cj，k｝（cj－1，k的近似信號(hào)）和｛dj，k｝（cj－1，k的細(xì)節(jié)信號(hào)）.其具體的計(jì)算公式為

對于公式（6）和公式（7），可以用更為簡單的圖形表示.記Aj＝｛cj，k｝，Dj＝｛dj，k｝表示公式（6）和公式（7）所表達(dá)的運(yùn)算，則Mallat分解算法可用圖1展示分解算法結(jié)構(gòu)示意圖.

信號(hào)經(jīng)過分解處理后，便要進(jìn)行信號(hào)的重構(gòu)，就要用到重構(gòu)算法.重構(gòu)算法作為分解算法的逆過程，其計(jì)算公式為

2 小波閾值去噪原理

目前比較常用的去噪方法有3種：模極大值去噪法、空域相關(guān)去噪法和小波閾值去噪法.本文的研究選用的是應(yīng)用最廣泛的小波閾值去噪法.

令基本的噪聲模型表示為

式中，e（n）為噪聲，σ為噪聲強(qiáng)度.

小波分解過程如圖2所示.一般來說，一維信號(hào)的閾值消噪過程可分為以下3個(gè)步驟進(jìn)行：

1）計(jì)算含噪信號(hào)的正交小波變換.選擇一個(gè)正交小波并確定一個(gè)小波分解的層次N，然后對信號(hào)x（n）進(jìn)行N層小波分解.

2）對小波系數(shù)進(jìn)行非線性閾值處理.為保持信號(hào)的整體形狀不變，保留所有的低頻系數(shù)，對第1到第N層的每一層高頻系數(shù)，采用硬閾值或軟閾值方法進(jìn)行量化處理.

3）進(jìn)行逆小波變換.根據(jù)小波分解的第N層的低頻系數(shù)和經(jīng)過量化處理后的第1層到第N層的高頻系數(shù)，做逆小波變換進(jìn)行重構(gòu)，得到恢復(fù)的原始信號(hào)的估計(jì)值.

3 閾值的量化和閾值選取

3.1 閾值的量化規(guī)則

硬閾值處理是把信號(hào)小波變換系數(shù)的絕對值與閾值比較，小于或等于閾值的小波系數(shù)變?yōu)榱?，大于閾值的點(diǎn)保持不變.軟閾值處理則是把比較后小波系數(shù)大于閾值的點(diǎn)變?yōu)樵擖c(diǎn)與閾值的差值.用公式表示為

硬閾值

軟閾值

3.2 閾值的選擇

小波去噪有4種自適應(yīng)閾值原則可供選擇（潘泉等，2005）：無偏似然估計(jì)原則、固定閾值原則、啟發(fā)式閾值原則和極值閾值原則.

1）無偏似然估計(jì)（sure）原則.無偏似然估計(jì)原則是一種基于stein無偏似然估計(jì)（二次方程）原理的自適應(yīng)閾值選擇.對一個(gè)給定的閾值t，先得到它的似然估計(jì)，再將似然t最小化，就得到了所選的閾值，它是一種軟閾值估計(jì)器.

2）固定閾值原則.固定閾值計(jì)算公式如下：

式中，n為信號(hào)x（i）的長度，T為閾值.

3）啟發(fā)式閾值原則.啟發(fā)式閾值原則是無偏似然估計(jì)原則與固定閾值原則的折中.如果信噪比很小，即按無偏似然估計(jì)處理；如果信噪比較大，則采用固定閾值原則處理.

4）極值閾值原則.極值閾值原則采用極大極小原則選擇閾值，它產(chǎn)生一個(gè)最小均方誤差的極值，而不是沒有誤差.計(jì)算公式如下：

4 小波閾值去噪的仿真實(shí)驗(yàn)分析

本節(jié)進(jìn)行小波閾值去噪的仿真實(shí)驗(yàn)分析，針對小波基的選取以及去噪閾值的不同對去噪效果的影響進(jìn)行對比分析，尋求可用于微震信號(hào)去噪分析的小波及去噪閾值.

為了更加精確以及量化評價(jià)去噪效果，定義信噪比（SNR）和均方差（RMSE）兩個(gè)參數(shù).設(shè)原信號(hào)為x（n），去噪后的信號(hào)為x（n）′，則信噪比的定義為

原始信號(hào)和去噪后信號(hào)的均方根誤差定義為

信噪比越高，均方根誤差越小，去噪信號(hào)就越接近原始信號(hào)，去噪效果就越好.

4.1 小波基的選取

金解放等（2007）指出Db4小波基能滿足巖石聲發(fā)射信號(hào)處理的要求.本文選用經(jīng)典的symlet小波對幾種不同類型的含噪信號(hào)進(jìn)行小波去噪，來檢驗(yàn)使用symlet小波進(jìn)行小波去噪的可行性.為了便于比較去噪結(jié)果，小波分解采用symlet6小波對信號(hào)進(jìn)行5層分解，采用啟發(fā)式閾值進(jìn)行去噪，閾值量化為軟閾值.

1）信噪比為5的沖擊（bump）信號(hào).含噪沖擊信號(hào)的原信號(hào)及去噪后信號(hào)如圖3所示.

2）信噪比為5的塊狀（block）信號(hào).含噪塊狀信號(hào)的原信號(hào)及去噪后信號(hào)如圖4所示.

3）信噪比為5的寬正弦（heavysin）信號(hào).含噪寬正弦信號(hào)的原信號(hào)及去噪后信號(hào)如圖5所示.各種信號(hào)去噪后的信噪比和均方根如表1所示.

表1 去噪后信號(hào)的信噪比和均方根Table 1 SNR and RMSE of denoised signal

從表1可以看出，選用symlet小波去噪可以對不同類型的含噪信號(hào)進(jìn)行有效地去噪，在均方差較小的情況下，極大地提高了信號(hào)的信噪比，較好地重現(xiàn)了含噪前的真實(shí)信號(hào)，可以用于信號(hào)的去噪應(yīng)用.

4.2 不同去噪閾值的去噪結(jié)果分析

為了便于比較不同閾值的消噪結(jié)果，選用4.1節(jié)中的symlet6小波，分解層數(shù)為5層，分別以無偏似然估計(jì)閾值、固定閾值、啟發(fā)式閾值和極值閾值對該節(jié)中信噪比為5的沖擊信號(hào)、塊狀信號(hào)以及寬正弦信號(hào)進(jìn)行消噪，閾值量化為軟閾值.限于篇幅，只給出沖擊信號(hào)的原信號(hào)及去噪后信號(hào)，如圖6所示.

含噪沖擊信號(hào)、塊狀信號(hào)和寬正弦信號(hào)4種閾值去噪后信號(hào)的信噪比和均方根分別如表2—4所示.

表2 沖擊信號(hào)去噪后的信噪比和均方根Table 2 SNR and RMSE of denoised bump signal

表3 塊狀信號(hào)去噪后的信噪比和均方根Table 3SNR and RMSE of denoised blocks signal

從表2—4中可以看出4種閾值去噪的均方差都較小，都極大地提高了信號(hào)的信噪比，有效地去除了噪聲.經(jīng)比較可以發(fā)現(xiàn)，對不同的含噪信號(hào)，無偏似然原則閾值降噪后的信噪比最高，同時(shí)均方根誤差也最小，在去噪時(shí)顯得更為有效，既能有效地去除了噪聲，同時(shí)也能保留信號(hào)里的高頻信息.因此，下文中利用小波閾值去噪對現(xiàn)場微震信號(hào)進(jìn)行去噪時(shí)將選用無偏似然原則閾值.

表4 寬正弦信號(hào)去噪后的信噪比和均方根Table 4 SNR and RMSE of denoised heavysin signal

5 大尺度巖體結(jié)構(gòu)微震監(jiān)測信號(hào)的小波閾值去噪

5.1 柿竹園礦微震監(jiān)測系統(tǒng)介紹

柿竹園多金屬礦位于湖南省郴州市蘇仙區(qū)境內(nèi)，自1987年采用分段鑿巖階段礦房法開采以來，因種種原因井下的采空區(qū)（礦房）未能進(jìn)行處理，截止2002年，井下留下占礦段約60%的礦柱礦量和近300萬m3的巨大采空區(qū)群，累計(jì)頂板暴露面積達(dá)3萬m2，連續(xù)頂板暴露面積近1萬m2，15m厚的連續(xù)條帶礦柱多處垮塌.由于采空區(qū)量大集中，暴露時(shí)間長，大爆破振動(dòng)影響頻繁和受破碎帶與巖體結(jié)構(gòu)面的影響，以及應(yīng)力集中導(dǎo)致礦柱破損等，殘留大尺度礦柱結(jié)構(gòu)穩(wěn)定性安全問題凸顯.

針對常規(guī)的地壓監(jiān)測缺點(diǎn)，柿竹園多金屬礦于2008年引進(jìn)國際先進(jìn)的全數(shù)字型多通道微震監(jiān)測系統(tǒng)，可進(jìn)行全范圍立體空間實(shí)時(shí)監(jiān)測.柿竹園礦多通道全數(shù)字型微震監(jiān)測系統(tǒng)為加拿大ESG公司第四代全數(shù)字型監(jiān)測系統(tǒng)，共計(jì)5臺(tái)Paladin數(shù)據(jù)采集儀（30個(gè)通道），攜帶30個(gè)單軸加速度傳感器.該系統(tǒng)在組成上由地表數(shù)據(jù)處理站、井下微震數(shù)據(jù)采集儀、井下傳感器三大部分組成.

5.2 微震監(jiān)測信號(hào)的小波閾值去噪

柿竹園微震系統(tǒng)自建立以來，采集了大量的微震波形信號(hào)，現(xiàn)已建成波形數(shù)據(jù)庫，為微震技術(shù)及理論的進(jìn)一步研究提供了寶貴的數(shù)據(jù).經(jīng)過長期的現(xiàn)場研究，已形成系統(tǒng)和有效的微震信號(hào)辨識(shí)方法，包括經(jīng)驗(yàn)分析法、理論分析法和綜合分析法等.但這只是信號(hào)處理的初步，辨識(shí)后的各種信號(hào)仍然含有各種隨機(jī)噪聲.噪聲對信號(hào)參數(shù)計(jì)算以及微震事件定位會(huì)產(chǎn)生不同程度的影響，微震信號(hào)還需進(jìn)一步的信號(hào)去噪研究.

本部分將選用現(xiàn)場采集到的幾種典型的微震信號(hào)，利用4.1節(jié)仿真實(shí)驗(yàn)的結(jié)論，采用symlet6小波及無偏似然原則閾值進(jìn)行去噪研究，閾值量化為軟閾值，小波分解層數(shù)為5層.需要指出的是，本部分的所有信號(hào)都采用傳統(tǒng)去噪方法中的帶通濾波器進(jìn)行了初步的去噪，獲得包含有效信號(hào)頻率段50—5 000Hz范圍內(nèi)的信號(hào).當(dāng)然，初步帶通去噪后信號(hào)中依然還包含有其它噪聲.

1）可定位巖體破裂事件.巖體破裂信號(hào)是指巖石或巖體在應(yīng)力的作用下產(chǎn)生破裂時(shí)產(chǎn)生的彈性波信號(hào)，是微震監(jiān)測的主要對象.包括巖體受力變形、產(chǎn)生裂紋、裂紋擴(kuò)展、產(chǎn)生裂縫、裂縫擴(kuò)展等過程產(chǎn)生的微震信號(hào)，也包括巖體錯(cuò)動(dòng)、片幫、崩落、冒頂、巖爆、大爆破后的余震、滑坡以及大規(guī)模巖體失穩(wěn)導(dǎo)致的礦震等微震信號(hào).可定位巖體破裂事件是指震源激起彈性波的能量較大，可以被多個(gè)傳感器感應(yīng)、接收并產(chǎn)生觸發(fā)，可由系統(tǒng)定位軟件計(jì)算出具體位置的事件.

圖7a所示的巖體破裂事件信號(hào)，是在柿竹園礦一次大爆破后由于礦體應(yīng)力重分布引起的能量較大的巖體破裂信號(hào)；圖7b為小波閾值去噪后的波形圖；圖7c為巖體破裂事件信號(hào)的小波分解各尺度波形圖，可以看到，各尺度都包含著明顯的有效微震信號(hào)，如果選用常用的低通或帶通濾波器將會(huì)遺失大量的有效信息；圖7d為各尺度細(xì)節(jié)小波系數(shù)作用閾值的過程，在每個(gè)尺度進(jìn)行閾值去噪，去除噪聲，保留有效信號(hào).

巖體破裂及爆破事件信號(hào)去噪后的信噪比和均方根如表5所示.

表5 巖石破裂及爆破去噪后信號(hào)的信噪比和均方根Table 5 SNR and RMSE of denoised rock mass fracture signal and denoised blast signal

圖8給出了可定位巖體破裂信號(hào)中的噪聲及其頻譜圖.由圖8b可以看到，隨機(jī)噪聲頻率分布范圍較廣，從幾十Hz到3 000Hz.如果為了去除噪聲，選擇3 000—5 000Hz的帶通濾波器，則單純對頻域的濾波有“一刀切”的缺陷，對帶通之外的頻譜不加區(qū)分的去除的同時(shí)也過濾掉了很多有效信號(hào).

2）爆破信號(hào).爆破信號(hào)是系統(tǒng)監(jiān)測到的主要事件之一.當(dāng)研究爆破對巖體穩(wěn)定的影響或爆破波在巖體中的傳播衰減較規(guī)律時(shí)，它是系統(tǒng)研究的對象，爆破信號(hào)能量大，一般都為可定為事件.

圖9a所示的爆破事件信號(hào)，是柿竹園礦某一次爆破的信號(hào)波形；圖9b為小波閾值去噪后的波形圖.限于篇幅，未給出小波分解各尺度波形圖以及各尺度細(xì)節(jié)小波系數(shù)作用閾值的過程.

圖10給出了爆破信號(hào)中的噪聲及其頻譜圖.由圖10b可以看到，隨機(jī)噪聲頻率分布范圍從幾十Hz到5000Hz，涵蓋了非常廣的頻率段.通過以上的分析可以發(fā)現(xiàn)，對于能量較大的定位事件，降噪后的信號(hào)消除了含噪信號(hào)中覆蓋的較強(qiáng)的隨機(jī)噪聲和毛刺現(xiàn)象，并且與原信號(hào)具有同等的光滑性；降噪后信號(hào)的信噪比和均方差均達(dá)到理想的數(shù)值，能有效地區(qū)分信號(hào)中的突變部分和噪聲，實(shí)現(xiàn)了在各個(gè)尺度上進(jìn)行去噪，較好地重現(xiàn)了含噪前的真實(shí)信號(hào).

3）不可定位巖體破裂事件.不可定位巖體破裂事件在這里指的是單通道巖體破裂事件或者事件只有被2—3個(gè)傳感器觸發(fā)的巖體破裂事件，事件能量較可定位事件小.不可定位巖體破裂事件的有效辨識(shí)對微震日常監(jiān)測具有極其重要的意義，它直接反應(yīng)各個(gè)傳感器臨近區(qū)域巖體穩(wěn)定性和地壓活動(dòng)性的相對程度.并且，由于是能量較小的破裂，得到時(shí)域上事件數(shù)或相對能量值的演變特征后，結(jié)合現(xiàn)場考察情況進(jìn)行分析總結(jié)，并尋求相關(guān)理論的合理解釋，這對研究和評價(jià)巖體失穩(wěn)過程具有重要的意義.

圖11a所示的微震監(jiān)測系統(tǒng)采集到的一段波形數(shù)據(jù)，從波形圖中可以初步判斷該段波形圖包含著巖體破裂事件，但由于能量較小，被噪聲強(qiáng)烈的干擾，巖體破裂事件的有效波形信號(hào)不明顯，其起跳和衰減的過程都不清晰；圖11b是小波閾值去噪后的信號(hào)，可以清楚地看到，該段波形包含3個(gè)能量較小的巖體破裂事件，小波閾值去噪實(shí)現(xiàn)了在低信噪比情況下有效地提取真實(shí)的巖體破裂信號(hào).

通過以上的分析可以看到，小波閾值去噪既可以對低信噪比的微震信號(hào)提取有效信號(hào)，也可以對頻率覆蓋范圍廣的微震信號(hào)在各尺度上提取有效信號(hào)，實(shí)現(xiàn)了對微震信號(hào)真實(shí)信號(hào)與噪聲信號(hào)的有效分離.小波閾值去噪特別適合于大尺度巖體結(jié)構(gòu)微震信號(hào)這一類非穩(wěn)定信號(hào)的去噪分析，是微震信號(hào)進(jìn)行有效提純和進(jìn)一步信號(hào)分析的有效工具.

6 討論與結(jié)論

大尺度巖體結(jié)構(gòu)的穩(wěn)定性分析是一個(gè)復(fù)雜的研究內(nèi)容，微震監(jiān)測技術(shù)是研究大尺度巖體結(jié)構(gòu)穩(wěn)定性的一個(gè)有力的手段.全數(shù)字多通道微震監(jiān)測系統(tǒng)的建立，實(shí)現(xiàn)了對工程現(xiàn)場大尺度巖體結(jié)構(gòu)的實(shí)時(shí)全過程監(jiān)測，解決了傳統(tǒng)聲發(fā)射技術(shù)現(xiàn)場工程應(yīng)用數(shù)據(jù)不完整性的問題.本文采用小波閾值去噪理論實(shí)現(xiàn)了微震信號(hào)的信噪分離.基于小波的信號(hào)分析，如微震信號(hào)特征參量的準(zhǔn)確計(jì)算、微震事件的精確定位等極具研究意義.可以預(yù)見的是，如果采用更為先進(jìn)的小波理論進(jìn)行微震信號(hào)的信號(hào)分析，將會(huì)有力地推進(jìn)微震監(jiān)測技術(shù)的發(fā)展.

本文得到以下結(jié)論：

1）通過Matlab仿真實(shí)驗(yàn)，symlet6小波滿足小波閾值去噪的要求，驗(yàn)證了使用symlet6小波進(jìn)行小波去噪的可行性.

2）4種閾值去噪的均方差都較小，都極大地提高了信號(hào)的信噪比，有效地去除了噪聲.對不同的含噪信號(hào)，無偏似然原則閾值降噪后的信噪比最高，同時(shí)均方根誤差也最小，在去噪時(shí)顯得更為有效.

3）小波閾值去噪特別適合于大尺度巖體微震信號(hào)這一類非穩(wěn)定信號(hào)的去噪分析，既可以對低信噪比的微震信號(hào)提取有效信號(hào)，也可以對頻率覆蓋范圍廣的微震信號(hào)在各尺度上提取有效信號(hào)，實(shí)現(xiàn)了對微震信號(hào)真實(shí)信號(hào)與噪聲信號(hào)的有效分離.

金解放，趙奎，王曉軍，趙康.2007.巖石聲發(fā)射信號(hào)處理小波基選擇的研究［J］.礦業(yè)研究與開發(fā)，27（2）：12－15.

李庶林.2009.試論微震監(jiān)測技術(shù)在地下工程中的應(yīng)用［J］.地下空間與工程學(xué)報(bào)，5（1）：122－128.

李庶林，尹賢剛，鄭文達(dá)，Trifu C.2005.凡口鉛鋅礦多通道微震監(jiān)測系統(tǒng)及其應(yīng)用研究［J］.巖石力學(xué)與工程學(xué)報(bào)，24（12）：2048－2053.

潘泉，張磊，孟晉麗，張洪才.2005.小波濾波方法及應(yīng)用［M］.北京：清華大學(xué)出版社：3－5.

王更峰.2006.基于小波分析的巖石聲發(fā)射信號(hào)處理技術(shù)［J］.礦業(yè)工程，4（5）：11－14.

許大為，潘一山，李國臻，王巖.2008.基于小波變換的礦山微震信號(hào)濾波方法研究［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，19（3）：74－77.

袁節(jié)平，胡靜云，周愛民，林峰，陳際經(jīng)，李庶林.2010.柿竹園多通道微震監(jiān)測系統(tǒng)的建立及其應(yīng)用［J］.礦業(yè)研究與開發(fā)，30（6）：12－14.

趙奎，鄧飛，金解放.2006.巖石聲發(fā)射Kaiser點(diǎn)信號(hào)的小波分析及其應(yīng)用初步研究［J］.巖石力學(xué)與工程學(xué)報(bào)，25（2）：3854－3858.

Coifman R R，Donoho D L.1994.Translation－invariant de－noising［M］∥Wavelets and Statistics.Springer Lecture Notes in Statistics 103.New York：Springer－Verlag：125－150.

Chang S G，Yu B，Vetterli M.2000.Adaptive wavelet thresholding for image denoising and compression［J］.IEEE Transation on Image Processing，9（9）：1532－1546.

Donoho D L，Johnstone I M.1994.Ideal sPatial adaptation via wavelet shrinkage［J］.Biometrika，81：425－455.

Mallat S，Hwang W L.1992.Singularity detection and processing with wavelet［J］.IEEE Transation on Information Theory，38（2）：617－643.

Johnstone L M，Silverman B W.1997.Wavelet threshold estimators for data with correlated noise［J］.J R Statist Soc，59（B）：319－351.

Zhang S，Salari E.2005.Image denoising using a neural network based non－linear filter in wavelet domain［J］.IEEE Transation on Image Processing，5（2）：989－992.

Xu Y，Weaver J，Healy M.1994.Wavelet transform domain filters：A sPatially selective noise filtration teehnique［J］.IEEE Transation on Image Processing，3（6）：747－758.

徐宏斌 現(xiàn)在華潤置地（福州）有限公司工作.2011年廈門大學(xué)建筑與土木工程學(xué)院土木工程系巖土工程專業(yè)畢業(yè)，獲工學(xué)碩士學(xué)位.曾主要從事聲發(fā)射與微地震監(jiān)測技術(shù)等方面的研究工作.

注：孟令媛、李宇彤、鄭建常、蔣海昆、王秀英、宋美琴等的簡介分別見本刊：Vol.33，No.4；Vol.30，No.4；Vol.29，No.4；Vol.17，No.4；Vol.31，No.6；Vol.33，No.5.

A study on method of signal denoising based on wavelet transform for micro－seismicity monitoring in large－scale rockmass structures

Xu Hongbin1）Li Shulin1），Chen Jijing2）
1）School of Architecture and Civil Engineering of Xiamen University，F(xiàn)ujan，Xiamen361005，China
2）Hunan Shizhuyuan Nonferrous Metals Company Limited，Hunan，Chenzhou423037，China

This paper applied wavelet denoising method to monitoring microseismicity in large－scale rockmass structure.The feasibility of using symlet6in wavelet denoising was validated with MATLAB simulation.Then four types of adaptive threshold rules for wavelet denoising are used to denoise three noisy signals.The result shows that the noise in signals can be filtered effectively with the four threshold rules and the Rigrsure threshold for wavelet denoising is more effective with the least mean square deviation and highest signal to noise ratio.Based on the multi－channel digital microseism monitoring system in Shizhuyuan mine，this paper applied wavelet denoising method to three different microseismic signals with the result of MATLAB simulation.The results show that the true microseismic signals can be recovered from the noisy signals by removing noise at every wavelet scale，even though noisy signals have low signal to noise ratio or include wide frequency range.The wavelet threshold denoising is suited especially to the denoising of microseismic monitoring signals in largescale rockmass structures.

wavelet de－noising；large－scale rockmass；microseismic technology；MATLAB simulation

10.3969／j.issn.0253－3782.2012.01.008

P315.9

A

徐宏斌，李庶林，陳際經(jīng).2012.基于小波變換的大尺度巖體結(jié)構(gòu)微震監(jiān)測信號(hào)去噪方法研究.地震學(xué)報(bào)，34（1）：85－96.

Xu Hongbin，Li Shulin，Chen Jijing.2012.A study on method of signal denoising based on wavelet transform for microseismicity monitoring in large－scale rockmass structures.Acta Seismologica Sinica，34（1）：85－96.

國家自然科學(xué)基金（10572122）資助.

2011－03－10收到初稿，2011－09－06決定采用修改稿.

e－mail：shulin.li＠163.com