趙 奎 宋 林 曾 鵬 謝文健 龔 囪 伍文凱 劉 宇

(1.江西理工大學資源與環(huán)境工程學院,江西 贛州 341000;2.江西省礦業(yè)工程重點實驗室,江西 贛州 341000;3.江西省地質(zhì)局生態(tài)地質(zhì)大隊,江西 南昌 330001)

聲發(fā)射監(jiān)測是礦山常用的無損監(jiān)測手段,在礦巖體、礦柱及充填體的穩(wěn)定性監(jiān)測方面均發(fā)揮了重要作用,聲發(fā)射監(jiān)測過程一般要求將傳感器盡可能貼近監(jiān)測對象,從而防止聲發(fā)射信號過度衰減[1]。但實際采礦工程中,常常會出現(xiàn)聲發(fā)射傳感器無法直接安放在監(jiān)測對象表面的情形。以采用充填采礦法的礦山為例,聲發(fā)射傳感器有時須安放在隔條礦柱或充填體側(cè),達到監(jiān)測目的。因此,研究通過充填體介質(zhì)監(jiān)測巖石聲發(fā)射信號和通過巖石介質(zhì)監(jiān)測充填體聲發(fā)射信號特征具有一定的實際意義。

巖石類材料內(nèi)部損傷演化過程中往往伴隨聲發(fā)射現(xiàn)象[2-3]。學者們對巖石類材料損傷破壞過程中的聲發(fā)射信號進行研究,探究其破壞失穩(wěn)聲發(fā)射前兆特征,并取得了豐富的成果。張省軍等[4]研究了單軸壓縮下粗?；◢弾r、大理巖和砂巖等巖石破壞過程的聲發(fā)射信號,提出可根據(jù)累計聲發(fā)射數(shù)快速增加、聲發(fā)射能量快速釋放、聲發(fā)射b值下降到最低值作為巖石破裂失穩(wěn)的前兆信息;張艷博等[5]通過巷道巖爆模擬試驗,對該過程中的聲發(fā)射信號進行聚類分析及神經(jīng)網(wǎng)絡(luò)識別,并建立了巖爆實時預(yù)警方法;曾鵬等[6]對粗砂巖單軸壓縮聲發(fā)射信號頻段占比隨應(yīng)力的變化特征進行了分析,通過特征頻段占比及聲發(fā)射幅值關(guān)聯(lián)維數(shù)演化可更準確地對巖石破壞前兆進行判識和預(yù)測;宋朝陽等[7]對不同粒徑的弱膠結(jié)砂巖進行單軸壓縮聲發(fā)射試驗,將聲發(fā)射信號主頻頻段增加作為其破壞前兆判識信息;張志博等[8]采用時—空維度聚類分析法分析了循環(huán)加卸載條件下煤樣破壞過程的聲發(fā)射信號,將無聲發(fā)射信號區(qū)域縮小、能量劇烈釋放作為煤樣失穩(wěn)破壞前兆判識信息;CHEN等[9]對不同節(jié)理角度的砂巖試件進行剪切試驗,研究了不同節(jié)理角度條件下砂巖剪切破壞過程中的裂紋演化規(guī)律及其聲發(fā)射信號特征;梁學健等[10]開展了不同灰砂比充填體單軸壓縮試驗,提出中、高頻信號明顯增加這一特征可作為充填體破壞前兆判識信息;趙奎等[11]基于顆粒流軟件建立了花崗巖等效晶體模型,對不同粒徑大小的花崗巖進行單軸壓縮聲發(fā)射模擬試驗,研究了晶粒大小對聲發(fā)射相對平靜期的影響規(guī)律,為相對平靜期現(xiàn)象應(yīng)用于巖石失穩(wěn)破壞的監(jiān)測預(yù)警提供了一定的依據(jù)。

綜上分析,學者們在不同巖性巖石及類巖石材料破壞過程的聲發(fā)射信號特征及其預(yù)警方面做了大量研究,并取得了豐碩成果,但關(guān)于通過充填體介質(zhì)監(jiān)測巖石聲發(fā)射信號或通過巖石介質(zhì)監(jiān)測充填體聲發(fā)射信號特征方面的研究涉及較少。為此,本研究在室內(nèi)開展砂巖及充填體—砂巖組合體單軸壓縮聲發(fā)射試驗,分別在充填體—砂巖組合體兩側(cè)各安放同類傳感器,研究兩者之間聲發(fā)射信號的差異性,并基于BP神經(jīng)網(wǎng)絡(luò)對充填體—砂巖組合體的聲發(fā)射信號類型進行識別分析,探究充填體—砂巖組合體破壞前兆判識信息,以期為礦山開采過程中的圍巖與充填體組合礦柱破壞失穩(wěn)監(jiān)測及預(yù)警提供依據(jù)。

1 試驗簡介

1.1 試件制備

試驗參照充填采礦法回采過程(圖1)分別制備砂巖試件與充填體—砂巖組合體試件,具體步驟如下:

(1)砂巖試件制備。試驗所用砂巖均取自同一塊結(jié)構(gòu)較完整的砂巖巖塊,巖塊取自江西贛南地區(qū),顏色為磚紅色,結(jié)構(gòu)致密,遇水無泥化現(xiàn)象。參照國際巖石力學學會(ISRM)建議的尺寸和精度,共制備6個標準圓柱體砂巖試件,尺寸為?50 mm×100 mm,編號為H1～H6。

(2)充填體—砂巖組合體試件制備。將制備的砂巖試件H4～H6對半切開,得到6塊半圓柱體砂巖試件,并將其置于標準型圓柱模具中,采用灰砂比1∶4、質(zhì)量分數(shù)74%的充填料漿進行澆筑,待其初凝后進行刮模、脫模處理,然后置于YH-35B型標準恒溫恒濕養(yǎng)護箱進行養(yǎng)護。基于實際開采過程中,礦山有時在充填體構(gòu)筑完成1～2周后,進行下一階段開采工作,故本次試驗充填體—砂巖組合體試件養(yǎng)護齡期為14 d。養(yǎng)護完成后,最終得到3個合格的充填體—砂巖組合體試件,編號為Z1～Z3(圖2)。

圖2 砂巖及充填體—砂巖組合體試件Fig.2 Specimens of sandstone and filling body-sandstone combination

本次試驗對砂巖、充填體—砂巖組合體試件進行單軸壓縮聲發(fā)射試驗,采集相應(yīng)的力學試驗數(shù)據(jù)及聲發(fā)射原始數(shù)據(jù)。

1.2 試驗設(shè)備及方法

本次試驗采用中國科學院武漢巖土力學研究所研發(fā)的RMT-150C巖石力學加載試驗系統(tǒng),試驗全程采用位移控制方式,加載速率為0.002 mm/s。

試驗聲發(fā)射采集系統(tǒng)采用PCI-Ⅱ型聲發(fā)射儀,AE傳感器選用2個型號相同的UT-1000傳感器,前置放大增益40 dB,采樣頻率為1MSPS,門檻值為40 dB。試驗過程中,將AE傳感器安裝于試件側(cè)面居中位置,安放位置事先用比例尺標記好點位并輕微打磨,然后在試件與傳感器接觸位置涂抹凡士林耦合劑使二者充分耦合接觸,并用膠帶將傳感器固定。為了對比分析,在充填體—砂巖組合體試件砂巖側(cè)與充填體側(cè)各布置1個聲發(fā)射傳感器,傳感器安裝于試件砂巖側(cè)(充填體側(cè))居中位置。試驗設(shè)備與試件安裝連接如圖3所示。

圖3 試驗設(shè)備與試件安裝連接示意Fig.3 Schematic of installation and connection of test equipment and test piece

1.3 單軸抗壓強度

充填體—砂巖組合體試件的時間—應(yīng)力曲線(圖4)、單軸抗壓強度與砂巖試件存在明顯差異。在應(yīng)力水平較低時,砂巖試件曲線表現(xiàn)為下凹彎曲狀態(tài),組合體試件曲線呈微傾斜緩慢上升;當加載至峰值應(yīng)力附近時,各試件均表現(xiàn)為上凸屈服狀態(tài)。同時,砂巖試件的單軸抗壓強度明顯高于組合體試件,砂巖試件單軸抗壓強度平均值為52.91 MPa,而組合體試件的平均單軸抗壓強度為20.55 MPa,僅為砂巖試件的38.85%。

圖4 單軸壓縮下不同類型試件時間—應(yīng)力曲線Fig.4 Time-stress curves of different types of specimens under uniaxial compression

2 聲發(fā)射信號差異性分析

2.1 聲發(fā)射振鈴計數(shù)

聲發(fā)射振鈴計數(shù)是聲發(fā)射信號超過某設(shè)定閾值的振蕩次數(shù),可反映聲發(fā)射信號強度與頻度,常用于評價聲發(fā)射活動性程度強弱[13]。充填體—砂巖組合體試件破壞過程中砂巖側(cè)與充填體側(cè)及砂巖試件破壞過程的聲發(fā)射振鈴計數(shù)與峰前應(yīng)力百分比關(guān)系如圖5所示。

圖5 聲發(fā)射振鈴計數(shù)與應(yīng)力百分比關(guān)系(試件H 1、Z1)Fig.5 Relationship between acoustic emission ring count and stress percentage (specimen H 1 and Z1)

加載初期,組合體試件砂巖部分壓密過程中與充填體部分摩擦、擠壓,造成損傷,試件兩側(cè)均監(jiān)測到大量的聲發(fā)射信號產(chǎn)生,砂巖試件僅有少量聲發(fā)射信號產(chǎn)生;當加載至峰值應(yīng)力的80%左右時,砂巖試件及組合體試件砂巖側(cè)聲發(fā)射振鈴計數(shù)均出現(xiàn)明顯的“激增”現(xiàn)象,而組合體試件充填體側(cè)聲發(fā)射振鈴計數(shù)“激增”現(xiàn)象相對不明顯;同時,臨近峰值應(yīng)力,砂巖試件有少量振鈴計數(shù)值較高的聲發(fā)射信號產(chǎn)生,而組合體試件砂巖側(cè)與充填體側(cè)峰值破壞前聲發(fā)射信號振鈴計數(shù)值均處于較低水平,與砂巖試件存在明顯差異。從整體上看,峰值應(yīng)力前砂巖試件中聲發(fā)射信號振鈴計數(shù)整體呈“平穩(wěn)—激增”變化趨勢,組合體試件砂巖側(cè)聲發(fā)射信號振鈴計數(shù)整體呈“W”狀變化,組合體試件充填體側(cè)聲發(fā)射信號振鈴計數(shù)整體呈“V”狀變化,相互之間存在明顯差異。

2.2 聲發(fā)射信號主頻

對各傳感器所采集到的聲發(fā)射信號原始波形文件進行快速傅里葉變換,得到各試件破壞過程的聲發(fā)射信號主頻。根據(jù)聲發(fā)射采樣率(1 MSPS),得到其奈奎斯特(Nyquist)頻率為500 kHz,用MATLAB自帶小波包將聲發(fā)射信號分解到第4層,分別提取了16個頻段,各頻段長度為31.25 kHz,分別為0～31.25 kHz,31.25～62.5 kHz,…,468.75～500 kHz[6]。分別得到砂巖、組合體(砂巖側(cè)與充填體側(cè)信號)破壞過程的聲發(fā)射信號主頻、應(yīng)力隨時間的變化關(guān)系如圖6所示。

圖6 聲發(fā)射信號主頻與應(yīng)力百分比關(guān)系(試件H 1、Z1)Fig.6 Relationship between dominant frequency of acoustic emission signal and stress percentage (specimen H 1 and Z1)

加載初期,組合體試件砂巖側(cè)與充填體側(cè)聲發(fā)射信號主頻主要集中在31.25～62.5 kHz頻段,同時存在少部分0～31.25 kHz與93.75～125 kHz頻段的聲發(fā)射信號,而砂巖試件中的聲發(fā)射信號主頻幾乎全部集中在31.25～62.5 kHz頻段。隨著應(yīng)力水平不斷增加,砂巖試件聲發(fā)射信號主頻逐漸由31.25～62.5 kHz 1個頻段向31.25～62.5 kHz、62.5～93.75 kHz及93.75～125 kHz 3個頻段演化,組合體試件砂巖側(cè)與充填體側(cè)聲發(fā)射信號主頻逐漸由0～31.25 kHz、31.25～62.5 kHz與93.75～125 kHz 3個頻段向31.25～62.5 kHz及93.75～125 kHz 2個頻段演化。峰值應(yīng)力前,砂巖試件聲發(fā)射信號主頻由原先3個頻段向6個頻段向演化,并產(chǎn)生少量250～281.25 kHz及281.25～312.5 kHz高頻段聲發(fā)射信號。組合體試件砂巖側(cè)聲發(fā)射信號主頻由原先2個頻段向0～31.25 kHz、31.25～62.5 kHz、62.5～93.75 kHz與93.75～125 kHz 4個頻段演化;充填體側(cè)聲發(fā)射信號主頻由原先2個頻段向0～31.25 kHz、31.25～62.5 kHz與93.75～125 kHz 3個頻段演化,且93.75～125 kHz頻段聲發(fā)射信號量明顯少于試件砂巖側(cè);同時能夠發(fā)現(xiàn),峰值應(yīng)力前,組合體試件砂巖側(cè)與充填體側(cè)均未監(jiān)測到高頻段聲發(fā)射信號產(chǎn)生,與砂巖試件存在明顯差異。

3 基于BP神經(jīng)網(wǎng)絡(luò)的聲發(fā)射信號類型識別及分析

3.1 聲發(fā)射信號類型識別

BP神經(jīng)網(wǎng)絡(luò)具有信號向前傳遞,誤差反向傳播的特點,也被稱為誤差反向傳播神經(jīng)網(wǎng)絡(luò),在信號識別等領(lǐng)域有著十分廣泛的應(yīng)用[5,14-15]。因此,本研究將大量已知的砂巖和充填體聲發(fā)射信號建立樣本庫,利用提取的特征參數(shù)對BP神經(jīng)網(wǎng)絡(luò)進行訓(xùn)練,誤差達到合理范圍后使用BP神經(jīng)網(wǎng)絡(luò)對組合體試件峰值破壞前的聲發(fā)射信號進行識別,進一步探究充填體—砂巖組合體破壞前兆判識信息,具體結(jié)果如圖7所示。

圖7 聲發(fā)射信號類型與應(yīng)力百分比關(guān)系(試件Z1)Fig.7 Relationship between acoustic emission signal type and stress percentage (specimen Z1)

由圖7可知:組合體試件破壞過程中的聲發(fā)射信號主要為砂巖類聲發(fā)射信號,充填體類聲發(fā)射信號較少。加載初期,組合體試件聲發(fā)射信號大部分為砂巖類聲發(fā)射信號,充填體類聲發(fā)射信號較少,說明此時組合體試件由砂巖部分起主要承載作用;超過峰值應(yīng)力的40%后,充填體類聲發(fā)射信號逐漸減少,試件聲發(fā)射信號幾乎全為砂巖類聲發(fā)射信號;達到峰值應(yīng)力的80%后,充填體類聲發(fā)射信號數(shù)逐漸增加,表明此時試件砂巖部分內(nèi)部微裂紋萌生擴展達到一定程度,承載能力降低,試件充填體部分承載作用增強[16];超過峰值應(yīng)力的90%后,試件砂巖類聲發(fā)射信號與充填體類聲發(fā)射信號數(shù)量均出現(xiàn)“突增”現(xiàn)象,說明此時試件內(nèi)部微裂紋不斷擴展、貫通,聲發(fā)射活動十分劇烈,預(yù)示試件即將破壞,可作為組合體試件破壞前兆特征之一。

3.2 不同類型聲發(fā)射信號占比分析

為進一步研究峰值應(yīng)力前組合體試件聲發(fā)射信號類型規(guī)律,對不同應(yīng)力百分比區(qū)間內(nèi)砂巖類聲發(fā)射信號占比進行了統(tǒng)計[6,17-18],結(jié)果如圖8所示。

圖8 不同應(yīng)力百分比區(qū)間內(nèi)砂巖類聲發(fā)射信號占比關(guān)系Fig.8 Proportion of acoustic emission signals of sandstone in different stress percentage intervals

峰值應(yīng)力前,組合體試件砂巖側(cè)聲發(fā)射信號中砂巖類聲發(fā)射信號占比較大,總體呈“降—升—降”變化趨勢。加載初期,砂巖側(cè)聲發(fā)射信號大部分為砂巖類聲發(fā)射信號,說明此時組合體試件由砂巖部分起主要承載作用;持續(xù)加載過程,試件逐漸被壓密,試件充填體部分逐漸開始起承載作用并產(chǎn)生少量充填體類聲發(fā)射信號,砂巖類聲發(fā)射信號占比逐漸下降;應(yīng)力持續(xù)增大,組合體試件砂巖部分內(nèi)部微裂紋逐漸萌生擴展,產(chǎn)生大量砂巖類聲發(fā)射信號,此階段試件充填體部分依然處于被壓密狀態(tài),只有少量充填體類聲發(fā)射信號產(chǎn)生,砂巖類聲發(fā)射信號占比逐漸上升;當加載至峰值應(yīng)力的70%～80%時,試件砂巖部分已產(chǎn)生較大程度損傷,充填體部分承載作用迅速增強[16],內(nèi)部微裂紋不斷萌生擴展,產(chǎn)生大量充填體類聲發(fā)射信號,砂巖類聲發(fā)射信號占比逐漸下降,預(yù)示試件即將破壞。

4 結(jié) 論

(1)充填體—砂巖組合體試件的時間—應(yīng)力曲線特征與砂巖試件存在明顯差異。組合體試件的平均單軸抗壓強度明顯小于砂巖試件,為砂巖的38.85%。

(2)對比分析充填體—砂巖組合體與砂巖試件破壞前的聲發(fā)射信號特征發(fā)現(xiàn),砂巖試件臨近破壞前有少量250.0～281.25 kHz及281.25～312.5kHz頻段的較高頻信號出現(xiàn),而組合體試件兩側(cè)均未出現(xiàn)此頻段的聲發(fā)射信號。且組合體試件中所采集到的聲發(fā)射信號振鈴計數(shù)量整體少于砂巖試件,差異明顯。

(3)峰值應(yīng)力前,組合體試件中砂巖一側(cè)的聲發(fā)射振鈴計數(shù)出現(xiàn)明顯的“激增”現(xiàn)象,同時產(chǎn)生大量93.75～125kHz頻段的聲發(fā)射信號;而充填體一側(cè)的聲發(fā)射振鈴計數(shù)“激增”現(xiàn)象相對不明顯,且只有極少量93.75～125 kHz頻段的聲發(fā)射信號。

(4)采用BP神經(jīng)網(wǎng)絡(luò)對組合體試件兩側(cè)采集到的信號進行識別,計算了不同應(yīng)力百分比區(qū)間內(nèi)的砂巖類信號的占比。組合體試件破壞過程中產(chǎn)生的聲發(fā)射信號以砂巖類信號為主;峰值應(yīng)力前,砂巖類信號的占比逐漸下降,充填體類信號的占比逐漸上升,這一現(xiàn)象可作為組合體試件破壞的前兆特征。