基于微震監(jiān)測(cè)試驗(yàn)的微震信號(hào)類別識(shí)別降噪方法研究

劉躍成，康躍明，趙智輝，周 創(chuàng)，周 水

(1.云南交投集團(tuán)投資有限公司，云南 昆明 675800； 2.中煤科工集團(tuán)重慶研究院有限公司，重慶 400039)

高埋深隧道以“高地應(yīng)力，高溫，高滲流性”，等顯著性特點(diǎn)，對(duì)其隧道開挖施工安全性產(chǎn)生較大影響[1-2]。如果圍巖狀態(tài)掌握不準(zhǔn)確，造成應(yīng)力積累和突然釋放，造成巖爆發(fā)生，不僅容易造成經(jīng)濟(jì)的損失，而且會(huì)使整個(gè)隧道陷入困境，甚至造成人員傷亡[3]。因此，開展對(duì)微震監(jiān)測(cè)信號(hào)識(shí)別研究，不僅可以通過微震監(jiān)測(cè)手段對(duì)圍巖狀態(tài)有更清晰地認(rèn)識(shí)，而且能為隧道開挖圍巖應(yīng)力理論和支護(hù)技術(shù)方面提供更可靠、更全面的資料。

目前，對(duì)于微震監(jiān)測(cè)技術(shù)的研究很多，主要是關(guān)于微震信號(hào)發(fā)生源的定位、能量識(shí)別以及巖爆預(yù)測(cè)。Ma 等學(xué)者[4]采用數(shù)值分析和微震監(jiān)測(cè)試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析，提出了判斷圍巖穩(wěn)定性的綜合預(yù)測(cè)方法。文獻(xiàn)[5]通過建立微震信號(hào)與地下工程開挖之間的關(guān)系，分析圍巖變形特征與位置。Young 等學(xué)者[6-7]證明了采用微震監(jiān)測(cè)技術(shù)來判斷圍巖開挖損傷特征的可行性，并提出了更好處理和解釋微震數(shù)據(jù)的非常規(guī)方法。Cai 等[8]采用大數(shù)據(jù)對(duì)微震監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析，提出了定量描述了巖體損傷程度的方法。Li Yue等[9]結(jié)通過整理和歸類西北地區(qū)沙漠地震資料，采用信號(hào)特征結(jié)構(gòu)原理對(duì)其噪聲去噪，改進(jìn)了去噪聲法。賈寶新等[10]采用TFA-DC方法解決了雙震相定位中S波初至到時(shí)拾取困難的問題，并且通過雙震相目標(biāo)函數(shù)，計(jì)算震出源坐標(biāo)。竇林名等[11]通過數(shù)值模擬、理論分析及現(xiàn)場(chǎng)實(shí)測(cè)的方法，分析了震動(dòng)波在巖層中的傳播歸類和礦震震動(dòng)波傳播過程中質(zhì)點(diǎn)峰值速度和震動(dòng)波能量的衰減規(guī)律。李邵軍等[12]采用微震監(jiān)測(cè)系統(tǒng)對(duì)中國錦屏地下工程開挖過程中巖層變化進(jìn)行了微震監(jiān)測(cè)，分析了圍巖松動(dòng)程度和危害以及針對(duì)巖爆可能性進(jìn)行了分析。

因此，基于微震監(jiān)測(cè)試驗(yàn)研究信號(hào)類別對(duì)采集信號(hào)的特征影響規(guī)律，提出通過波形信號(hào)特征[13]、體變勢(shì)大小、視體積以及視應(yīng)力大小來對(duì)所監(jiān)測(cè)的微震信號(hào)進(jìn)行歸類劃分的手段；首先針對(duì)監(jiān)測(cè)的波形特征進(jìn)行歸類分析，并與實(shí)際現(xiàn)場(chǎng)施工場(chǎng)景所可能采集的微震信號(hào)類別進(jìn)行對(duì)比分析；然后針對(duì)體變勢(shì)大小、視體積以及視應(yīng)力大小綜合特征進(jìn)行展開分析，最終通過現(xiàn)場(chǎng)實(shí)際情況來進(jìn)行驗(yàn)證。

1 微震監(jiān)測(cè)原理及系統(tǒng)介紹

擬試驗(yàn)場(chǎng)地位于鳳慶縣鳳山鎮(zhèn)安石村—勐佑鎮(zhèn)中河村的安石隧道，該隧道區(qū)屬中山地貌，地形起伏較大，最大埋深約453.11 m，為分離式特長隧道，雙向長度分別為5 338 m和5 263 m。隧址區(qū)下伏基巖前半段為燕山期侵入花崗巖，區(qū)域地質(zhì)較穩(wěn)定,隧道圍巖中殘存的構(gòu)造應(yīng)力在隧道開挖中易使堅(jiān)硬巖石產(chǎn)生巖爆。石隧道施工中采用上下臺(tái)階爆破開挖方式?，F(xiàn)場(chǎng)隧道開挖施工中由微震監(jiān)測(cè)系統(tǒng)對(duì)其巖層內(nèi)部及隧道周圍信號(hào)進(jìn)行監(jiān)測(cè)。

1.1 微震監(jiān)測(cè)系統(tǒng)組成

微震監(jiān)測(cè)系統(tǒng)主要由3部分組成[14-15]：傳感器、數(shù)據(jù)采集系統(tǒng)及用于數(shù)據(jù)實(shí)時(shí)監(jiān)測(cè)與處理的微震服務(wù)器，如圖1所示。傳感器接收到與巖體破裂有關(guān)的波形信號(hào)(模擬電信號(hào))以后，由連接電纜將信號(hào)傳輸至數(shù)據(jù)采集儀，采集儀把接收到的模擬信號(hào)A/D轉(zhuǎn)換為數(shù)字信號(hào)，并經(jīng)過光纜通過網(wǎng)絡(luò)交換機(jī)將信號(hào)傳輸?shù)轿挥谒淼纼?nèi)二襯下方的微震監(jiān)測(cè)與處理服務(wù)器上，此過程如圖2所示。

圖1 微震監(jiān)測(cè)原理示意Fig.1 Principle of microseismic monitoring

圖2 微震監(jiān)測(cè)系統(tǒng)網(wǎng)絡(luò)拓?fù)銯ig.2 Network topology of microseismicmonitoring system

1.2 微震監(jiān)測(cè)原理

20世紀(jì)30年代后期，文獻(xiàn)[14-15]在深埋硬巖礦中開展聲波研究時(shí)發(fā)現(xiàn)了微震現(xiàn)象。高應(yīng)力作用下的巖體受到開挖卸荷等擾動(dòng)時(shí)，伴隨著裂縫的產(chǎn)生及擴(kuò)展，巖體內(nèi)集聚的彈性應(yīng)變能以應(yīng)力波的形式釋放并傳播，即為微震事件。微震活動(dòng)包含了巖體變形破壞的豐富信息，通過布設(shè)的微震監(jiān)測(cè)系統(tǒng)采集、識(shí)別并定位微震信號(hào)，可獲得微裂隙位置及巖體損傷程度，由此可以對(duì)巖體失穩(wěn)災(zāi)害進(jìn)行預(yù)警，從而確保圍巖安全穩(wěn)定。與傳統(tǒng)監(jiān)測(cè)技術(shù)相比，微震監(jiān)測(cè)技術(shù)具有遠(yuǎn)距離、動(dòng)態(tài)、三維、整體、實(shí)時(shí)監(jiān)測(cè)等特點(diǎn)，其最大優(yōu)點(diǎn)是能夠定位出圍巖潛在破壞區(qū)域的位置。

微震監(jiān)測(cè)技術(shù)是一種被動(dòng)監(jiān)測(cè)方法，在隧道中的應(yīng)用主要包括以下3個(gè)步驟：①在潛在圍巖失穩(wěn)區(qū)附近安裝傳感器和數(shù)據(jù)采集設(shè)備，傳感器臺(tái)陣的布設(shè)以其能夠空間覆蓋關(guān)注區(qū)域?yàn)樽罴?②檢測(cè)和處理微震事件，包括噪聲濾除、微震信號(hào)識(shí)別與初至拾取及微震事件定位;③根據(jù)現(xiàn)場(chǎng)地質(zhì)條件及工程施工活動(dòng)分析并解釋微震數(shù)據(jù)，根據(jù)微震活動(dòng)有關(guān)特征預(yù)警潛在的圍巖失穩(wěn)區(qū)域，確保工程安全。

2 微震監(jiān)測(cè)方案設(shè)計(jì)

2.1 設(shè)計(jì)基本原則

微震監(jiān)測(cè)不同于常規(guī)應(yīng)力、變形監(jiān)測(cè)，微震監(jiān)測(cè)范圍是一個(gè)空間的體，微震監(jiān)測(cè)方案設(shè)計(jì)時(shí)應(yīng)考慮微震信號(hào)分析與工程災(zāi)害預(yù)警所必需的條件。因此，微震監(jiān)測(cè)方案設(shè)計(jì)時(shí)應(yīng)遵循以下原則：①滿足“微震源定位需要有4個(gè)傳感器監(jiān)測(cè)到有效信號(hào)”的必要條件；②考慮微震源定位原理，盡可能避免微震源定位的“盲區(qū)”；③考慮微震監(jiān)測(cè)的特殊性及傳感器的類型與性能，根據(jù)巖體結(jié)構(gòu)特點(diǎn)，在確保監(jiān)測(cè)范圍內(nèi)對(duì)巖石破裂事件定位及時(shí)、準(zhǔn)確的前提下，盡可能擴(kuò)大監(jiān)測(cè)范圍；④盡量避免環(huán)境噪聲對(duì)微震信號(hào)的影響；⑤考慮現(xiàn)場(chǎng)布置的可操作性、設(shè)備銜接和走線的方便，確保線路安全，確保監(jiān)測(cè)數(shù)據(jù)連續(xù)、準(zhǔn)確；⑥在條件允許的情況下，盡可能確保監(jiān)測(cè)對(duì)象在傳感器陣列之內(nèi);⑦各傳感器應(yīng)協(xié)同作業(yè)原則：充分利用傳感器協(xié)同工作的特點(diǎn)，提高監(jiān)測(cè)系統(tǒng)自我容錯(cuò)能力，當(dāng)某區(qū)域傳感器工作不正常時(shí)，其他區(qū)域的傳感器仍可保障該區(qū)域的基本監(jiān)測(cè)。

2.2 傳感器安裝設(shè)計(jì)

傳感器的布設(shè)以其能夠空間覆蓋監(jiān)測(cè)區(qū)域?yàn)樽罴?。然而，在深埋隧道中進(jìn)行微震監(jiān)測(cè)時(shí)，傳感器通常只能布設(shè)于工作面之后，而微震事件主要集中在隧道工作面附近，因此，傳感器陣列很難或不可能實(shí)現(xiàn)對(duì)監(jiān)測(cè)區(qū)域的三維覆蓋。在這種情況下，傳感器應(yīng)盡可能靠近監(jiān)測(cè)區(qū)域。

微震傳感器現(xiàn)場(chǎng)安裝的質(zhì)量，很大程度上影響著微震監(jiān)測(cè)結(jié)果的準(zhǔn)確性。①應(yīng)確保傳感器與巖體的有效耦合，提高微震信號(hào)的信噪比，以便后續(xù)數(shù)據(jù)的處理及解譯；②應(yīng)確保傳感器安裝到設(shè)計(jì)位置，保證傳感器坐標(biāo)位置的準(zhǔn)確性，提高微震事件的定位精度。

安石隧道中微震傳感器布設(shè)于一個(gè)斷面，每個(gè)斷面分別于不同位置處(隧道左側(cè)和隧道右側(cè))布設(shè)6個(gè)傳感器，具體布設(shè)如圖3所示。

圖3 傳感器布設(shè)Fig.3 Sensor layout diagram

傳感器現(xiàn)場(chǎng)安裝方式主要分為孔內(nèi)安裝和孔外安裝2種方法，此次采用孔內(nèi)安裝方法?？變?nèi)安裝時(shí)先將傳感器和回收裝置連接，主要用于傳感器的回收利用，并將其置于3 m深的孔底，而后用回收桿將其擰緊，確保傳感器與巖體充分耦合。

3 監(jiān)測(cè)信號(hào)規(guī)律研究

3.1 微震監(jiān)測(cè)過程中爆破信號(hào)

通過此次實(shí)驗(yàn)結(jié)果，可得知爆破信號(hào)的能量較大，觸發(fā)傳感器數(shù)量較多，部分傳感器接收的信號(hào)呈滿量程狀態(tài)，信號(hào)波形被削弱，波形呈現(xiàn)雙峰形。爆破信號(hào)衰減較快，尾波不發(fā)育，持續(xù)時(shí)間約為240～380 ms，頻譜最大頻率分布區(qū)間為0～800 Hz，多數(shù)分布在0～10 Hz。微震監(jiān)測(cè)系統(tǒng)在2021年6月2日—2021年7月30日監(jiān)測(cè)到的現(xiàn)場(chǎng)爆破采集事件體變視統(tǒng)計(jì)如圖4所示。研究得知爆破信號(hào)的體變勢(shì)分布在0～510 m3,且大多集中在10～90 m3，體變勢(shì)分布區(qū)間較為集中。爆破采集事件視體積統(tǒng)計(jì)如圖5所示。

圖4 爆破信號(hào)體變勢(shì)統(tǒng)計(jì)Fig.4 Statistical chart of blasting signal volume variation

圖5 爆破信號(hào)視體積統(tǒng)計(jì)Fig.5 Visual volume statistics of blasting signal

通過分析研究可知，爆破信號(hào)的視體積數(shù)量級(jí)最大，分布在0～45 000 000 m3，且多數(shù)分布在1 000～7 000 000 m3，最為分散。爆破采集事件視應(yīng)力統(tǒng)計(jì)如圖6所示。研究得知爆破信號(hào)視應(yīng)力數(shù)量級(jí)較大，分布在0～650 kPa，多數(shù)在50～450 kPa，較為分散。

圖6 爆破信號(hào)視應(yīng)力統(tǒng)計(jì)Fig.6 Visual stress statistics of blasting signal

3.2 微震監(jiān)測(cè)過程中打炮眼信號(hào)

通過分析研究，打炮眼信號(hào)持續(xù)時(shí)間為270～380 ms，頻譜最大頻率分布區(qū)間為190～230 Hz區(qū)間，多數(shù)分布在190～210 Hz。打炮眼采集事件體變勢(shì)統(tǒng)計(jì)如圖7所示。通過分析研究，得知打炮眼信號(hào)的體變勢(shì)大多集中在0～0.10 m3，數(shù)值最少，體變勢(shì)分布區(qū)間較為集中。

圖7 打炮眼信號(hào)體變勢(shì)統(tǒng)計(jì)Fig.7 Statistics of the signal volume of blasting boreholes

現(xiàn)場(chǎng)打炮眼采集事件視體積統(tǒng)計(jì)如圖8所示。通過分析研究，得知爆破信號(hào)的視體積數(shù)量級(jí)較大，分布在0～300 000 m3，且多數(shù)分布在0～100 000 m3，較為分散。

圖8 打炮眼信號(hào)視體積統(tǒng)計(jì)Fig.8 Signal apparent volume statistics of blasting boreholes

打炮眼采集事件視應(yīng)力統(tǒng)計(jì)如圖9所示。通過分析研究，得知打炮眼信號(hào)的視應(yīng)力數(shù)量級(jí)較大，分布區(qū)間為0～35 kPa，多數(shù)分布在2.5～20 kPa，較為集中。

圖9 打炮眼信號(hào)視應(yīng)力統(tǒng)計(jì)Fig.9 Statistics of apparent stress of blasting boreholes signal

3.3 電氣干擾信號(hào)波形分析

通過分析研究，打炮眼信號(hào)持續(xù)時(shí)間為270～330 ms，頻譜最大頻率分布區(qū)間為0～55 Hz，多數(shù)分布在48～55 Hz，極少數(shù)分布在10 Hz以下。

現(xiàn)場(chǎng)電氣干擾采集事件體變視統(tǒng)計(jì)如圖10所示。通過分析研究，得知電氣干擾信號(hào)的體變勢(shì)分布在0～0.17 m3,且大多集中在0～0.5 m3，數(shù)值最少，體變勢(shì)分布區(qū)間較為集中。

圖10 電氣干擾信號(hào)體變勢(shì)統(tǒng)計(jì)Fig.10 Statistical chart of electrical interference signal variant

現(xiàn)場(chǎng)電氣干擾采集事件視體積統(tǒng)計(jì)如圖11所示。通過分析研究，得知電氣干擾的視體積數(shù)量級(jí)較大(0～6 500 000 m3)，且多數(shù)分布在0～60 000 m3，較為分散。

圖11 電氣干擾信號(hào)視體積統(tǒng)計(jì)Fig.11 Statistics of electrical interference signal apparent volume

現(xiàn)場(chǎng)電氣干擾采集事件視應(yīng)力統(tǒng)計(jì)圖如圖12所示。通過分析研究，得知電氣干擾信號(hào)的視應(yīng)力數(shù)量級(jí)較大，為0～80 kPa，多數(shù)分布在0～20 kPa，較為集中。

圖12 電氣干擾信號(hào)視應(yīng)力統(tǒng)計(jì)Fig.12 Statistics of electrical interference signal apparent stress

3.4 巖石破裂微震監(jiān)測(cè)信號(hào)

通過分析研究，巖石破裂持續(xù)時(shí)間為300～430 ms，頻譜最大頻率分布區(qū)間為0～200 Hz，多數(shù)分布在50～60 Hz和180～200 Hz。

監(jiān)測(cè)到的現(xiàn)場(chǎng)巖石破裂采集事件體變勢(shì)統(tǒng)計(jì)如圖13所示。通過分析研究，得知巖石破裂產(chǎn)生的微震信號(hào)的體變勢(shì)分布在0～110m3,且大多集中在0～10 m3，數(shù)值較大，體變勢(shì)分布區(qū)間較為集中。

圖13 微震信號(hào)體變勢(shì)統(tǒng)計(jì)Fig.13 Statistical chart of microseismic signal volume variability

監(jiān)測(cè)到的現(xiàn)場(chǎng)巖石破裂采集事件視應(yīng)力統(tǒng)計(jì)如圖14所示。通過分析研究，得知巖石破裂產(chǎn)生的微震信號(hào)的視體積數(shù)量級(jí)較大，分布在0～17 500 000 m3，且多數(shù)分布在0～25 000 m3，較為集中。

圖14 微震信號(hào)視體積統(tǒng)計(jì)Fig.14 Statistics of apparent volume of microseismic signals

監(jiān)測(cè)到的現(xiàn)場(chǎng)爆破采集事件視應(yīng)力統(tǒng)計(jì)如圖15所示。

圖15 微震信號(hào)視應(yīng)力統(tǒng)計(jì)Fig.15 Statistics of apparent stress of microseismic signals

通過分析研究，得知巖石破裂產(chǎn)生的微震信號(hào)的視應(yīng)力數(shù)量級(jí)較大，分布區(qū)間為0～5 500 kPa，多數(shù)分布在1～5 kPa，較為集中。

4 結(jié)論

(1)隨著隧道的開挖，隧道周邊及巖層內(nèi)部發(fā)生變化產(chǎn)生微震信號(hào)，其中包括施工過程的干擾信號(hào)(爆破信號(hào)、打炮眼信號(hào)、電氣信號(hào))和巖石破裂信號(hào)。因此，提出了通過波形信號(hào)特征、體變勢(shì)大小、視體積以及視應(yīng)力大小來對(duì)所監(jiān)測(cè)的微震信號(hào)進(jìn)行歸類劃分的方法。

(2)爆破信號(hào)體變勢(shì)的數(shù)值最大，分布在0～510 m3,且大多集中在10～90 m3。打炮眼信號(hào)的體變勢(shì)大多集中在0～0.10 m3，數(shù)值最少，電氣干擾信號(hào)的體變勢(shì)分布在0～0.17m3,且大多集中在0～0.5 m3，數(shù)值較少。巖石破裂產(chǎn)生的微震信號(hào)的體變勢(shì)分布在0～110 m3,且大多集中在0～10 m3，數(shù)值較大。

(3)爆破信號(hào)頻譜最大頻率多數(shù)分布在0～10 Hz區(qū)間。打炮眼頻譜最大頻率分布區(qū)間為190～230 Hz，多數(shù)分布在190～210 Hz。極少數(shù)分布在10 Hz以下。現(xiàn)場(chǎng)電氣干擾信號(hào)頻譜最大頻率多數(shù)分布在48～55 Hz?，F(xiàn)場(chǎng)巖石破裂產(chǎn)生的微震信號(hào)頻譜最大頻率分布區(qū)間為0～200 Hz，多數(shù)分布在50～60 Hz和180～200 Hz。

(4)爆破信號(hào)的視體積數(shù)量級(jí)最大，多數(shù)分布在1 000～7 000 000 m3，最為分散；炮眼視體積數(shù)量級(jí)較大，分布在0～300 000 m3，且多數(shù)分布在0～100 000 m3，較為分散。電氣干擾的視體積數(shù)量級(jí)較大，分布在0～6500000 m3，且多數(shù)分布在0～60 000 m3，較為分散。巖石破裂產(chǎn)生的微震信號(hào)的視體積數(shù)量級(jí)較小，多數(shù)分布在0～25 000 m3，較為集中。

(5)爆破信號(hào)的視應(yīng)力數(shù)量級(jí)最大，多數(shù)分布在50～450 kPa，較為分散。打炮眼、電氣干擾及巖石破裂產(chǎn)生的視應(yīng)力信號(hào)較為集中。打炮眼信號(hào)的視應(yīng)力多數(shù)分布在2.5～20.0 kPa。電氣干擾信號(hào)的視應(yīng)力多數(shù)分布在0～20 kPa。巖石破裂產(chǎn)生的微震信號(hào)的視應(yīng)力多數(shù)分布在1～5 kPa。

首次考慮信號(hào)源對(duì)產(chǎn)生的波形信號(hào)特征、體變勢(shì)大小、視體積以及視應(yīng)力大小差異性特征，考慮以上因素特征進(jìn)行歸類與現(xiàn)場(chǎng)實(shí)測(cè)信號(hào)類型吻合良好，對(duì)今后的微震信號(hào)識(shí)別有很好的指導(dǎo)意義。