王金貴，張 蘇，胡維喜，張佳慶

(1. 福州大學(xué)環(huán)境與資源學(xué)院，福建 福州 350116； 2. 中國礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院，北京 100083； 3. 福州大學(xué)紫金礦業(yè)學(xué)院，福建 福州 350116； 4. 國網(wǎng)安徽省電力公司電力科學(xué)研究院，安徽 合肥 230601)

0 引言

歷經(jīng)長時間、高強(qiáng)度的開采，我國大部分礦井進(jìn)入了深部開采時期. 隨之而來的是地應(yīng)力和煤層瓦斯壓力明顯增大，地質(zhì)條件也日益復(fù)雜，煤巖動力災(zāi)害日益嚴(yán)重[1-4]. 動力災(zāi)害實質(zhì)是煤巖體原有應(yīng)力平衡狀態(tài)被打破所致，已有研究[3-11]表明，該過程中往往伴隨著電磁及微震的異常波動，且信號特征與煤巖變形破裂過程密切相關(guān).

任學(xué)坤等[5]通過對預(yù)制裂紋巖板進(jìn)行單軸壓縮實驗，發(fā)現(xiàn)預(yù)制裂紋巖板電位和電磁輻射隨載荷的變化而變化. Kong等[6]研究高溫處理后的砂巖變形破裂過程中的電磁信號特征，發(fā)現(xiàn)電磁信號與應(yīng)力值正相關(guān)，且經(jīng)過不同溫度處理的試件伴隨的電磁信號存在差異. 姚精明等[7]通過煤巖單軸壓縮實驗，研究發(fā)現(xiàn)該過程中電磁信號的多重分形譜寬度ΔDq與試樣所受的應(yīng)力水平密切相關(guān). Song等[8]研究了煤巖受載破裂過程中電磁信號頻譜的變化特征，發(fā)現(xiàn)隨著所受應(yīng)力的增加，電磁信號主頻從0～60 kHz下降至0～20 kHz. 李成武等[9]通過利用霍普金森沖擊壓桿對煤體施加沖擊載荷，測試分析了煤體沖擊破裂過程中的電磁特征. 對于煤巖破裂微震信號特征的研究，朱權(quán)潔等[10]利用煤與瓦斯突出模擬實驗系統(tǒng)及微震監(jiān)測系統(tǒng)，研究發(fā)現(xiàn)突出過程中的微震信號時頻特征具有階段性和差異性. 雷文杰等[11]研究了頂板、煤壁破裂過程伴隨的微震信號響應(yīng)特征，提出利用微震事件及其主頻帶能量趨勢進(jìn)行采動應(yīng)力轉(zhuǎn)移、調(diào)整的預(yù)測. 郭超等[12]利用一種隨信號自適應(yīng)改變參數(shù)的信號濾波器，進(jìn)行礦山實測震動信號的預(yù)處理，以消除低頻基線漂移及工頻干擾的影響. 陸菜平等[13]分析了沖擊地壓、煤層泄壓爆破和頂板泄壓爆破3種過程中伴隨的微震信號頻譜特征，發(fā)現(xiàn)它們的頻域都低于200 Hz，但主頻差異較大. 綜上所述，煤體破裂電磁及微震信號已被廣泛應(yīng)用于礦井動力災(zāi)害的預(yù)測，并取得較好效果，但現(xiàn)有研究較少考慮采掘作用的動力擾動對煤體損傷的影響.

通過統(tǒng)計礦井歷年的沖擊礦壓案例，梁冰等[14]發(fā)現(xiàn)93.3%的沖擊事件是被動力擾動誘發(fā)的. 蔡成功等[15]統(tǒng)計了我國有明確作業(yè)方式記錄的8 669次突出事例，發(fā)現(xiàn)由放炮、打鉆等作業(yè)方式擾動誘發(fā)的突出有8 362次，占比高達(dá)96. 5%. 掘進(jìn)爆破、放炮落煤及鉆探作業(yè)等是煤炭采掘過程中常見的動力擾動形式，這些擾動打破原有的應(yīng)力平衡狀態(tài)，引起載荷遷移和應(yīng)力集中現(xiàn)象； 應(yīng)力波在煤巖體中傳播造成裂隙的產(chǎn)生與發(fā)展，使其損傷變形，內(nèi)摩擦力減小，承載能力降低，最終造成煤巖體失穩(wěn)而觸發(fā)煤巖動力災(zāi)害. 特別在大采深時期，巨大的地應(yīng)力使煤巖體積聚了很大的能量，其自身已是臨界穩(wěn)定狀態(tài)，顯著增加了動力擾動觸發(fā)煤巖動力災(zāi)害的可能性. 因此，掌握煤巖體在采掘擾動作用下破壞失穩(wěn)時的電磁及微震信號特征，對于災(zāi)害預(yù)警中區(qū)分動力擾動下的正常和異常信號具有重要意義.

去噪是分析信號特征前的關(guān)鍵環(huán)節(jié)，雷文杰等[11]、李學(xué)龍等[16]、楊楨等[17]利用小波變換對微震信號去噪，分析其頻譜特征. 小波變換雖具有較強(qiáng)時頻分析能力，但應(yīng)用中小波基函數(shù)和分解層數(shù)的選擇對去噪成效影響很大，導(dǎo)致去噪結(jié)果不穩(wěn)定且有效性低[18]. Huang等[19]提出的經(jīng)驗?zāi)B(tài)分解算法可很好地應(yīng)對隨機(jī)非平穩(wěn)信號，該算法不需選擇基函數(shù)及分解層次，而是經(jīng)迭代自適應(yīng)，把原信號分解成多尺度的內(nèi)蘊(yùn)模態(tài)函數(shù)分量，該方法已被成功用于處理振動信號、電磁信號和診斷機(jī)械故障等[20-22].

本研究利用ZDKT-1型煤巖動力災(zāi)害監(jiān)測系統(tǒng)對平煤十礦己15-21180機(jī)巷掘進(jìn)工作面進(jìn)行動力災(zāi)害危險性實時監(jiān)測，分析現(xiàn)場煤巖體在強(qiáng)烈的采掘動力擾動下的電磁及微震信號特征，探究該破裂過程中電磁和微震信號的耦合關(guān)系.

1 測試現(xiàn)場及設(shè)備概況

1.1 測試現(xiàn)場簡介

爆破測試地點為平煤集團(tuán)十礦己15-21180機(jī)巷掘進(jìn)工作面，該礦主采煤層(丁組、戊組和己組)均發(fā)生過不同強(qiáng)度的突出，煤巖動力災(zāi)害嚴(yán)重威脅著該礦的生產(chǎn)作業(yè). 測試工作面開采的己15煤層厚度在1.6～2.3 m之間，煤層傾角8°～20°，最大煤層瓦斯壓力為2.4 MPa，煤層瓦斯含量13～20 m3/t.

1.2 監(jiān)測設(shè)備及布置簡介

采用中國礦業(yè)大學(xué)(北京)生產(chǎn)的ZDKT-1型煤巖動力災(zāi)害監(jiān)測系統(tǒng)(見圖1)對前述掘進(jìn)工作面進(jìn)行動力災(zāi)害危險性實時監(jiān)測，該系統(tǒng)通過配套的電磁接收天線及微震傳感器采集相應(yīng)數(shù)據(jù). 為了避免監(jiān)測系統(tǒng)對施工的影響，同時也為降低作業(yè)對設(shè)備的干擾，在保證振動、電磁傳感器有效監(jiān)測范圍的基礎(chǔ)上，將兩類傳感器安置在距掘進(jìn)面迎頭20 m處，詳見圖2.

電磁接收天線的磁棒直徑為100 mm，1 000匝線圈繞在該磁棒上，內(nèi)阻為21.8 Ω，電感為11.39 mH，電容為13.61 μF. 接收天線的有效接收頻段為0～2.2 kHz，在該頻段的信號能量級在10-7V2以上. 微震傳感器采用COLIBRYS公司生產(chǎn)的SF1500MEMS電容性加速度計，頻響范圍為0～1.5 kHz，對弱信號可提高到5 kHz. 動態(tài)范圍為120 dB，可受沖擊限制在1.5 kg以內(nèi).

圖1 ZDKT-1型煤巖動力災(zāi)害實驗?zāi)M系統(tǒng)Fig.1 ZDKT-1 type of coal and rock dynamic disasters experiment simulation system

圖2 平煤十礦掘進(jìn)工作面測點布置圖Fig.2 Layout of measuring points on excavating working face of Pingmei ten mine

2 信號去噪方法

2.1 現(xiàn)場信號噪聲特征

為了便于選擇信號去噪方法，需先分析現(xiàn)場電磁及微震信號背景噪聲的特征.

2.1.1 現(xiàn)場電磁噪聲

電磁信號監(jiān)測中的干擾源分成兩類：采集系統(tǒng)的內(nèi)部噪聲和外部噪聲，前者主要是元器件的散粒噪聲以及熱噪聲，后者主要是50 Hz工頻干擾及其諧波干擾. 圖3為截取的一段現(xiàn)場電磁背景噪聲及其頻譜圖. 從圖3中可以看出，現(xiàn)場電磁背景噪聲的幅值在1.5 mV以下，持續(xù)波動且無明顯變化. 從頻譜圖可以看出電磁背景噪聲在620、780和1 420 Hz有明顯的變化，其中主頻在780 Hz，這與實驗室背景噪聲基本為50 Hz的倍頻規(guī)律不同. 這是由于電力行業(yè)的頻率標(biāo)準(zhǔn)是50 Hz，頻率在48～52 Hz都屬于正常范圍，實驗室使用的為城市用電，較為規(guī)范； 而煤礦掘進(jìn)機(jī)使用的為1.14 kV三相電，煤礦井下高壓用電頻率不能保證為50 Hz.

圖3 現(xiàn)場電磁背景噪聲及其頻譜Fig.3 Electromagnetic signal background noise in the mine

2.1.2 現(xiàn)場微震噪聲

采集的微震信號中混雜了一些幅值穩(wěn)定的噪聲，這些噪聲可能來自地層的規(guī)律性振動及采集系統(tǒng)自帶噪聲. 圖4為截取的5 s現(xiàn)場微震噪聲及其快速傅里葉變換頻譜圖，由圖4可以看出，現(xiàn)場微震信號背景噪聲在0.5 mV左右持續(xù)波動，且無明顯變化. 從頻譜圖可以看出微震背景噪聲在120、500、650、780和1 420 Hz有明顯的變化，其中主頻在1 420 Hz.

圖4 現(xiàn)場微震噪聲及其頻譜Fig.4 Micro-seismic signal background noise in the mine

2.2 信號去噪方法

通過前述現(xiàn)場噪聲特征的分析，發(fā)現(xiàn)現(xiàn)場電磁及微震背景噪聲均為具有多個主頻的寬頻擾動，結(jié)合當(dāng)前常用去噪方法的優(yōu)缺點，選擇希爾伯特-黃變換(Hilbert-Huang transform，HHT)對信號進(jìn)行處理，其核心是借助整體經(jīng)驗?zāi)Ｊ椒纸夥?ensemble empirical mode decomposition，EEMD)分解重構(gòu)信號，主要有兩個作用：去除疊加波和使波形更加對稱[23].

EEMD的基本原理[24-25]是：先找出原信號X(t)的所有極值點，分別連接各極大、極小值點得到上包絡(luò)線Xmax(t)和下包絡(luò)線Xmin(t)，計算上、下包絡(luò)線的平均值線m1(t)，用X(t)減m1(t)得到h1(t). 再將h1(t)按前述步驟處理，即得到h11(t)=h1(t)-m11(t)，m11(t)為h1(t)的上、下包絡(luò)線均值, 如下式所示，重復(fù)k次，若h1k(t)和h1(k-1)(t)滿足SD小于0.25，則終止運(yùn)算.

2.3 信號能量表達(dá)

電磁和微震信號能量值是衡量煤巖損傷程度的重要指標(biāo)，該值可通過巴什瓦(Parseval)定理(式(2))求取，即利用前述EEMD的分解結(jié)果，求除噪聲模量以外的各模量能量值之和[25].

其中：i表示模量數(shù)；m是采樣點數(shù)；Ei，j(tj)為模量i在節(jié)點j處的頻帶能量值;xj，k為重構(gòu)信號fi，j(tj)的采樣點幅值.

2.4 信號去噪效果

利用前述去噪方法對監(jiān)測信號進(jìn)行去噪處理，限于篇幅，只列舉一組去噪前后的信號圖進(jìn)行對比分析，如圖5所示. 從圖5可以看出，通過前述方法，原信號中背景噪聲得到明顯減弱，且有效地保留了信號特征，原信號幅值未出現(xiàn)明顯下降.

圖5 電磁原始信號及其降噪后信號Fig.5 Electromagnetic signals before and after denoising

3 電磁和微震信號關(guān)聯(lián)性測試結(jié)果

利用前述監(jiān)測系統(tǒng)，監(jiān)測2012年9月22日至2012年9月26日期間該掘進(jìn)工作面10次爆破作業(yè)的電磁和微震信號，觀察發(fā)現(xiàn)在爆破作業(yè)后炮孔附近的煤體均發(fā)生了失穩(wěn)破壞，且在這些爆破掘進(jìn)作業(yè)中均存在明顯的電磁和微震信號.

3.1 電磁及微震信號

利用前述EEMD去噪方法，對現(xiàn)場的電磁及微震信號進(jìn)行降噪處理. 圖6為去噪后的電磁信號與微震信號對應(yīng)圖(左側(cè)為電磁，右側(cè)為微震). 分析各次作業(yè)過程中的電磁和微震信號出現(xiàn)時間點，發(fā)現(xiàn)兩類信號具有時間上的同步性，即煤體破裂損傷過程中產(chǎn)生的微震和電磁異常幾乎是同步的.

分析信號特征發(fā)現(xiàn)，現(xiàn)場爆破過程中的電磁信號存在明顯脈沖特性，單邊幅值在14 mV左右； 而微震信號與電磁信號相比，波形存在差異，微震信號波形為更具對稱性的簡諧振動，單邊幅值也明顯大于電磁信號，達(dá)到100 mV左右. 總體而言，兩類信號的幅值要明顯大于室內(nèi)煤體破裂實驗，這主要是由于煤巖破裂的速率及體量不在一個量級導(dǎo)致的.

圖6 去噪后電磁及微震信號對應(yīng)圖Fig.6 Electromagnetic signals and micro-seismic signals after denoising

3.2 電磁與微震信號關(guān)聯(lián)性

前述分析發(fā)現(xiàn)，同次作業(yè)過程的電磁和微震信號具有時間上的同步性. 為了分析兩類信號的單邊幅值、能量等參數(shù)的相關(guān)性，對去噪后信號的幅值及能量進(jìn)行統(tǒng)計，并分別進(jìn)行相關(guān)性分析(見圖7).

從圖7可以看出，現(xiàn)場兩類信號幅值的線性擬合相關(guān)性系數(shù)(0.756)要小于能量線性擬合關(guān)系(0.966)，分析認(rèn)為現(xiàn)場出現(xiàn)這種現(xiàn)象的原因是由于兩類信號在波形及頻率上的差異，微震信號比電磁信號更易出現(xiàn)尖峰型的波形，但這并不影響兩類信號在能量上符號較好的線性擬合關(guān)系.

圖7 煤巷爆破掘進(jìn)電磁和微震振幅及能量關(guān)系Fig.7 Amplitude and power relations of electromagnetic and micro-seismic of blasting excavation site

4 結(jié)語

以動力擾動下煤體破裂電磁、微震信號特征及關(guān)聯(lián)性為研究對象，基于煤巖動力災(zāi)害監(jiān)測系統(tǒng)，開展了爆破擾動煤體電磁和微震信號的現(xiàn)場試驗. 主要得到以下結(jié)論：1) 受供電頻率波動的影響，煤礦井下電磁和微震背景噪聲為具有多個主頻的寬頻擾動，且不嚴(yán)格符合50 Hz的倍頻； 2) 現(xiàn)場爆破過程中的電磁信號存在明顯脈沖特性，單邊幅值在14 mV左右； 而微震信號波形為更具對稱性的簡諧振動，單邊幅值達(dá)到100 mV左右； 3) 煤體破裂電磁和微震兩類信號具有時間上的同步性； 4) 電磁和微震信號能量符合較好的線性關(guān)系，擬合相關(guān)性系數(shù)達(dá)0.966，因兩類信號在波形及頻率上的差異，微震信號比電磁信號更易出現(xiàn)尖峰型的波形，導(dǎo)致兩類信號幅值的線性擬合相關(guān)性系數(shù)較低(0.756). 研究成果對礦井煤巖動力災(zāi)害的電磁、微震預(yù)測工作具有重要意義.

[1] 姜福興，史先鋒，王存文，等. 高應(yīng)力區(qū)分層開采沖擊地壓事故發(fā)生機(jī)理研究[J]. 巖土工程學(xué)報，2015，37(6)：1 123-1 131.

[2] 楊偉利，姜福興，楊鵬，等. 特厚煤層沖擊地壓重復(fù)發(fā)生的機(jī)理研究[J]. 巖土工程學(xué)報，2015，37(11)：2 045-2 050.

[3] 張宏偉，朱峰，韓軍，等. 沖擊地壓的地質(zhì)動力條件與監(jiān)測預(yù)測方法[J]. 煤炭學(xué)報，2016，41(3)：545-551.

[4] 何學(xué)秋，竇林名，牟宗龍，等． 煤巖沖擊動力災(zāi)害連續(xù)監(jiān)測預(yù)警理論與技術(shù)[J]. 煤炭學(xué)報，2014，39(8)：1 485-1 491.

[5] 任學(xué)坤，王恩元，李忠輝． 預(yù)制裂紋巖板破壞電位與電磁輻射特征的實驗研究[J]． 中國礦業(yè)大學(xué)學(xué)報，2016，45(3)：440-446．

[6] KONG B, WANG E Y, LI Z H,etal. Electromagnetic radiation characteristics and mechanical properties of deformed and fractured sandstone after high temperature treatment[J]. Engineering Geology，2016(209)：82-92．

[7] 姚精明，董文山，閆永業(yè)，等． 受載煤巖體電磁輻射動態(tài)多重分形特征[J]． 煤炭學(xué)報，2016，41(6)：1 429-1 433．

[8] SONG D Z, WANG E Y, SONG X Y,etal. Changes in frequency of electromagnetic radiation from loaded coal rock[J]. Rock Mech Rock Eng，2016(49)：291-302.

[9] 李成武，楊 威，徐曉萌，等. 煤體沖擊破壞超低頻/極低頻電磁異常特征分析[J]. 煤炭學(xué)報，2014，39(10)：2 014-2 021．

[10] 朱權(quán)潔，李青松，李紹泉，等. 煤與瓦斯突出試驗的微震動態(tài)響應(yīng)與特征分析[J]. 巖石力學(xué)與工程學(xué)報，2015，34(增刊2)：3 813-3 821.

[11] 雷文杰，王兆豐，韓真理. 采動效應(yīng)微震特征及其變化規(guī)律[J]. 巖石力學(xué)與工程學(xué)報，2015，34(4)：795-803.

[12] 郭超，陳鎣. 礦山震動信號預(yù)處理及其對時頻的影響[J]. 煤炭學(xué)報，2016，41(7)：1 859-1 865．

[13] LU C P, DOU L M, LIU B,etal. Microseismic low-frequency precursor effect of bursting failure of coal and rock[J]. Journal of Applied Geophysics, 2012,79(79): 55-63.

[14] 梁冰，汪北方，李剛，等. 動力擾動誘發(fā)巷幫煤體沖擊礦壓機(jī)制研究[J]. 安全與環(huán)境學(xué)報，2014，14(2)：13-18．

[15] 蔡成功，王佑安. 煤與瓦斯突出一般規(guī)律定性定量分析研究[J]. 中國安全科學(xué)學(xué)報，2004，14(6)：109-112．

[16] 李學(xué)龍，李忠輝，王恩元，等. 礦山微震信號干擾特征及去噪方法研究[J]. 中國礦業(yè)大學(xué)學(xué)報，2015，44(5)：788-822．

[17] 楊楨，李艷，李鑫，等． 改進(jìn)小波變換的煤巖電磁輻射信號去噪方法[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(自然科學(xué)版)，2015，34(3)：410-413．

[18] ZHANG L, BAO P, WU X. Multiscale LMMSE-based image de-noising with optimal wavelet selection[J]. IEEE Transaction on Circuits System Video Technology，2005，15(4)：469- 481.

[19] HUANG N E, SHEN Z, LONG S R,etal. The empirical mode de-composition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proceedings of the Royal Society of London，1998，454(1971)：903- 995.

[20] 王金貴，張?zhí)K. 基于頻域約束獨立成分分析的經(jīng)驗?zāi)B(tài)分解去噪方法[J]. 煤炭學(xué)報, 2017, 42(3): 621-629.

[21] WU FJ, QU L S. An improved method for restraining the end effect in empirical mode decomposition and its applications to the fault diagnosis of large rotating machinery[J]. Journal of Sound and Vibration. 2008，314(3/4/5) ：586-602．

[22] 李成武，董利輝，王啟飛，等． 煤巖微弱電磁信號的噪聲源識別及去噪方法[J]. 煤炭學(xué)報，2016，41(8)：1 933-1 940．

[23] WU Z H，HUANG N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. World Scientific，2009，1(1)：1-41.

[24] 李夕兵，凌同華，張義平. 爆破震動信號分析理論與技術(shù)[M]. 北京：科學(xué)出版社，2009：82-84.

[25] 張?zhí)K，王金貴，艾立昊． 水泥材料靜爆過程中的微震信號特征[J]. 振動、測試與診斷，2016，36(6)：1 078-1 084．