常新科, 吳順川*, 程海勇 , 王勇兵, 夏志遠 , 汪小東

(1. 昆明理工大學(xué)國土資源工程學(xué)院, 昆明 650093; 2. 云南馳宏鋅鍺股份有限公司會澤礦業(yè)分公司, 曲靖 654200; 3. 長沙有色冶金設(shè)計研究院有限公司, 長沙 410019)

在深部地下工程中,大多數(shù)隧道、巷道及硐室會承受較大的拉應(yīng)力或結(jié)構(gòu)損傷進而發(fā)生壓拉失穩(wěn)破壞,巖體的穩(wěn)定性監(jiān)測是保障深部巖體工程安全建設(shè)的必要手段[1-3]。聲發(fā)射(acoustic emission,AE)是巖石內(nèi)部結(jié)構(gòu)發(fā)生變化釋放的瞬態(tài)彈性波,AE信號特征與巖石的力學(xué)性質(zhì)及破裂機制具有密切關(guān)系[4]。因此,深入分析巖石在壓拉荷載下破壞的AE特性,有助于科學(xué)認識其演化規(guī)律,對實際工程穩(wěn)定性研究和監(jiān)測預(yù)警更具指導(dǎo)意義。

學(xué)者們從AE事件率[5-6]、能率[7-8]、b值[9-10]等時域特征參數(shù)方面入手,展開了對巖石破壞過程的一系列研究,取得了豐盛的成果。但這些AE時域特征參數(shù)只是對信號單一特征的統(tǒng)計描述,僅能提供有限的破裂信息。為了彌補研究上的不足,近年來關(guān)于AE頻域頻譜信息的研究熱度逐年增長,獲得了巖石受力狀態(tài)、結(jié)構(gòu)、物理力學(xué)特性等豐富信息。胡英等[11]從白云巖破壞的AE頻域特征出發(fā),分析了巖石破壞的前兆信息。Wang等[12]分析了細砂巖破壞的AE頻域信息,認為其破壞模式與波形主頻及幅值相關(guān)。Zhang等[13]進行了不同類別巖石的單軸壓縮AE試驗,認為AE高頻對應(yīng)剪切破壞,低頻對應(yīng)拉伸破壞。這些研究成果證實了從AE頻域頻譜信息角度,能夠獲得巖石的損傷機制。但主頻僅代表AE頻譜中最高幅值對應(yīng)的頻率,對于巖石這種復(fù)雜的非均勻材質(zhì)而言,只考慮主頻顯然是不合適的,容易忽略AE信號的全局性特征,時域內(nèi)的信息仍值得關(guān)注。Shiotani等[14]通過RA值(上升時間與幅值的比值)表征巖石破壞模式,表明高RA值對應(yīng)剪切破壞,低RA值對應(yīng)張拉破壞?；炷敛牧贤鶎F值(振鈴計數(shù)與持續(xù)時間的比值)的降低,視為破壞模式由拉伸向剪切的轉(zhuǎn)變[15]。Du等[16]進行了多組不同力學(xué)路徑的AE試驗,表明RA-AF值具有對巖石破裂模式進行分類的性質(zhì)。吳順川等[17]分析了隧道內(nèi)的AE監(jiān)測信號,表明RA-AF值具有表征巖石破壞的可行性。因此,RA-AF值可以作為從AE信號時域特征直接獲取巖石破裂模式的重要參數(shù)。

盡管對于巖石破壞過程中AE損傷表征已有廣泛研究,但其關(guān)注點主要局限在時域?qū)用婊蝾l域?qū)用鎯?nèi)的定性研究上,對二者耦合的定量研究仍然較少。由于深部巖石的復(fù)雜性,單一AE指標(biāo)不足以起到定量表征作用。如果能聯(lián)合AE頻-時域信息深入了解巖石微觀破壞機制,將有利于認識深部巖石從起裂發(fā)展至宏觀破裂的裂紋演變過程及破壞特性,與巖石內(nèi)部的微觀結(jié)構(gòu)變化起到緊密聯(lián)系。

據(jù)此,現(xiàn)對深部灰?guī)r進行巴西圓盤劈裂AE實驗,分析AE實測全波形特征,聯(lián)合頻-時域內(nèi)的主頻及RA-AF信息,探討深部巖石在張拉破壞下的AE特性及損傷表征,在一定程度上促進深部巖體工程破壞分析方法由定性分析向定量分析的轉(zhuǎn)變。

1 實驗設(shè)計

1.1 灰?guī)r試件制備

實驗所用巖石為云南某地下鋅鍺礦大塘組深部灰?guī)r,現(xiàn)場取芯時均按同一方向鉆取(井口標(biāo)高2 308 m,取芯孔標(biāo)高為1 104 m)。按照國際巖石力學(xué)與巖石工程學(xué)會(International Society for Rock Mechanics,ISRM)建議,將現(xiàn)場取回的巖芯加工成直徑50 mm ×厚25 mm的巴西圓盤試樣,編號BS,加工精度滿足試驗要求,取芯現(xiàn)場及巴西圓盤試樣,如圖1所示。

圖1 取芯現(xiàn)場及巴西圓盤試樣Fig.1 Coring site and specimen

1.2 實驗方案及參數(shù)設(shè)置

巴西圓盤劈裂AE試驗采用線荷載加載方式[18],所需試驗力由YAW-300型微機控制電液伺服壓力機提供,預(yù)加荷載為200 N,加載速率為0.1 MPa/s。AE信號通過DS5-16B/C型16通道AE監(jiān)測系統(tǒng)采集,采集參數(shù)基于前期開展的大量AE實驗研究,按表1設(shè)置。

在試件表面固定了6個涂抹硅脂的PAC-Nano30型AE傳感器,其中兩個布設(shè)在試件正面劈裂中心兩側(cè),剩余4個對稱布設(shè)在試樣圓環(huán)面上,分別與水平線呈±30°。AE傳感器布設(shè)方案及室內(nèi)試驗現(xiàn)場,如圖2所示。

表1 AE系統(tǒng)采集參數(shù)設(shè)置Table 1 AE system acquisition parameter setting

圖2 AE傳感器布設(shè)方案及試驗現(xiàn)場Fig.2 AE sensor installation scheme and test site

2 試驗結(jié)果分析

2.1 AE實測全波形特征

巖石在微觀尺度上的裂紋萌生、擴展和聚合變形,是巖石因應(yīng)力作用發(fā)生破壞的重要機制,這種與結(jié)構(gòu)變化機制直接相關(guān)的源,稱為AE源。AE源的本征性貯存在向外釋放的瞬態(tài)彈性波內(nèi),其幅值變化與巖石損傷強度直接關(guān)聯(lián)。而巖石損傷量的活躍程度則可通過統(tǒng)計單位時間內(nèi)的波形事件率進行測量。

通過對巖石破壞全過程的AE實測波形進行適當(dāng)分析,可以為巖石穩(wěn)定性監(jiān)測提供更多細節(jié)[19]?；?guī)r的巴西圓盤劈裂AE實驗結(jié)果表明,6個性能相同的AE傳感器采集到的波形存在差異,經(jīng)分析認為主要與AE源的距離、方向、傳播路徑及衰減特性相關(guān)。為此,將全部傳感器共同采集到同一信號進行耦合,選擇到達時間最早的波形作為實測AE波形,繪制了灰?guī)r在巴西圓盤劈裂條件下的AE實測全波形圖,如圖3所示。AE實測全波形圖中的不同顏色波形,代表了不同的觸發(fā)時刻。

分析圖3可知,通過AE實測全波形特征,可以將灰?guī)r在巴西圓盤劈裂條件下的損傷過程劃分為3個階段(階段Ⅰ～Ⅲ),依次對應(yīng)巖石微損傷萌生階段,裂紋擴展階段和結(jié)構(gòu)破壞階段。為了更清晰地表述階段Ⅰ～Ⅲ內(nèi)的波形特征,將各階段內(nèi)的AE極值波形(幅值最低和幅值最高的波形)進行了繪制,如圖4所示。

圖3 實測AE全波形圖Fig.3 AE full waveform

圖4 不同階段的AE極值波形Fig.4 AE extreme waveforms at different stages

圖4(a)和圖4(b)為分別階段Ⅰ內(nèi)幅值最低(No.6)和最高(No.32)的波形,圖4(c)和圖4(d)分別階段Ⅱ內(nèi)幅值最低(No.59)和最高(No.85)的波形,圖4(e)和圖4(f)分別為階段Ⅲ內(nèi)的幅值最低(No.277)和最高(No.792)的波形,其中波形的藍色區(qū)域為AE事件的持續(xù)時間,橘色區(qū)域為撞擊定義時間(hit definition time,HDT),黃色區(qū)域為撞擊閉鎖時間(hit lockout time,HLT)。HDT確保了AE系統(tǒng)采集到的信號僅為一次撞擊,HLT則避免了單個撞擊信號中包含其他撞擊信號,消除了反射波的干擾,保障了撞擊的完整性。不同階段的AE極值波形幅值比,如表2所示。

表2 不同階段的AE極值波形幅值比Table 2 AE polar waveform amplitude ratio at different stages

分析圖4和表2可以看出,低幅值波形多為連續(xù)型信號,幅值梯度變化小,高幅值波形則多為突發(fā)型信號,幅值梯度衰減快。隨試驗荷載的增加,低幅值信號的持續(xù)時間段明顯加長,極值波形幅值比顯著增加。這些貫穿試驗全過程低幅值波形,主要與微觀層面上的位錯、晶體孿生、晶界面移動相關(guān),相反,逐級遞增的高幅值波形,則主要與宏觀層面上的礦物顆粒、結(jié)構(gòu)面的滑移和分離相關(guān)。

階段Ⅰ,灰?guī)rAE實測波形較為離散,相鄰信號間隔較大,且幅值較低。在外部荷載作用下AE源釋放的AE波形較少,巖石微損傷活躍度較低,使得AE事件率維持在較低水平,AE累積事件數(shù)曲線表現(xiàn)出緩慢增長特征。這是由于灰?guī)r在裂紋擴展前期,發(fā)生的變化主要是內(nèi)部應(yīng)力狀態(tài)的重新分布,巖石結(jié)構(gòu)仍處于穩(wěn)定狀態(tài),巖石抵抗外力破壞的能力仍然較強。

階段Ⅱ,灰?guī)rAE實測波形表現(xiàn)出連貫密集的分布特征,但幅值仍然較低,AE事件率較前一階段有所增加,AE累積事件數(shù)曲線增長速率出現(xiàn)明顯提高。這表明巖石內(nèi)部局部區(qū)域已經(jīng)形成應(yīng)力集中,此前的微損傷演化成內(nèi)部裂紋,并在應(yīng)力集中區(qū)域得到擴展,使得巖石內(nèi)部結(jié)構(gòu)出現(xiàn)變化,損傷量迅速積累疊加。

階段Ⅲ,灰?guī)rAE實測波形表現(xiàn)出更加連貫密集的特征,幅值相比前一階段顯著增高,AE事件率進一步增加,特別是在試驗59 s左右出現(xiàn)了極高值,造成AE累積事件數(shù)曲線豎直向上攀升。這主要是因為巖石內(nèi)部結(jié)構(gòu)發(fā)生的劇烈變化,巖石抵抗外力破壞的能力消失,AE源的活躍度和強度大幅增加,損傷裂紋向宏觀破裂面聚攏,巖石發(fā)生斷裂。值得注意的是,在試驗?zāi)┪?AE事件率及累積事件曲線不再大幅增加,出現(xiàn)了1～2 s的“平靜期”,經(jīng)分析認為,由于巖石宏觀裂隙地形成,阻礙了瞬態(tài)彈性波地傳播。

2.2 AE信號主頻演化特征

通過快速傅里葉變換可以將信號從時域角度轉(zhuǎn)化到頻域角度,獲得實測波形的二維頻譜信息。二維頻譜圖中最大幅值對應(yīng)的頻率即為主頻,是AE頻域分析中表征波形特性的一個重要參數(shù)。圖4中極值波形(No.6、No.32、No.59、No.85、No.277和No.792)的二維頻譜圖,如圖5所示。

分析圖5可知,根據(jù)Nyquist采樣定律,AE采樣頻率為2.5 MHz,可以獲得0～1 250 kHz的頻域信息?；?guī)r的AE波形二維頻譜主要特征集中在0～500 kHz內(nèi),其中No.85和No.792為高頻信號,其余為低頻信號,灰?guī)r在階段Ⅰ～Ⅱ內(nèi)的AE頻譜多為較窄的單峰結(jié)構(gòu),進入階段Ⅲ后,頻譜增寬并逐漸向雙峰或多峰結(jié)構(gòu)轉(zhuǎn)化。這主要與巖石的損傷量及損傷程度密切相關(guān)。在階段Ⅰ～Ⅱ,巖石發(fā)生的變化主要是微損傷的萌生和裂紋的擴展,其損傷形式較單一。進入階段Ⅲ后,試驗臨近斷裂,其內(nèi)部結(jié)構(gòu)發(fā)生破壞,損傷形式增加,使得頻率成分逐漸復(fù)雜。

為了獲得更全面的分析結(jié)果,計算了實驗獲得的全部波形的AE主頻,同時結(jié)合試樣的時間-應(yīng)力曲線,繪制了灰?guī)r斷裂全過程的AE主頻演化特征圖,如圖6所示。

分析圖6可知,灰?guī)rAE主頻演化階段與實測AE全波形演化階段一致,可分為相同的階段Ⅰ、階段Ⅱ和階段Ⅲ,對應(yīng)巖石微損傷萌生階段,裂紋擴展階段,結(jié)構(gòu)破壞階段。按照主頻分布及聚類特征,又可分為4個頻帶,頻帶1(1～15 kHz)、頻帶2(33～55 kHz)、頻帶3(142～155 kHz)和頻帶4(307～358 kHz),其中頻帶1和頻帶2屬于低頻帶,頻帶3和頻點4屬于高頻帶。對于巖石AE而言,不同類型的AE源對應(yīng)著不同的破壞模式,不同的破壞模式產(chǎn)生不同頻率的信號[20]。

以BS-1為例,可以看出頻帶1貫穿試驗全過程,對應(yīng)著微觀層面上的位錯、晶體孿生,頻帶2在階段Ⅰ僅有零星幾個,在階段Ⅱ增加,在階段Ⅲ聚集,對應(yīng)著微觀層面上的晶界面移動,頻帶3僅在階段Ⅲ出現(xiàn),對應(yīng)著宏觀層面上的礦物顆粒和結(jié)構(gòu)面的分離,頻帶4主要分布在階段Ⅱ和階段Ⅲ內(nèi),對應(yīng)著宏觀層面上礦物顆粒和結(jié)構(gòu)面的滑移。臨近巖石斷裂前,頻帶的增加意味著破壞模式的復(fù)雜化。

圖5 AE波形二維頻譜圖Fig.5 2D spectrogram of AE waveform

圖6 AE主頻演化特征圖Fig.6 AE dominant frequency evolutionary feature

2.3 AE信號RA-AF分布特性

RA值越小時波形梯度越大,表明脈沖越窄,但RA值僅能刻畫上升時間段內(nèi)波形,對信號整體特征關(guān)注不足。AF值補充了描述波形的全局特征,在持續(xù)時間內(nèi)信號越過的門限的次數(shù)越多,表明頻率越高。

拉伸破壞信號具有窄脈沖,高頻率的特性,剪切破壞信號則具有寬脈沖,低頻率的特性[21],如圖8(a)所示。據(jù)此,研究者們根據(jù)RA值與AF值之間的關(guān)系對AE信號進行分類,將具有低RA值、高AF值特征信號視為剪切事件,將高RA值、低AF值特征的信號視為張拉事件,如圖8(b)所示。

不同AE主頻帶的RA-AF散點密度分布圖,如圖9所示。RA-AF的分布范圍,如表3所示。

分析圖9可知,RA值核心密度區(qū)域為0～8 ms/V,AF值核心密度區(qū)域為0～500 kHz,為便于分析,將RA-AF值的核心密度區(qū)域進行了放大觀察。可以看出,絕大多數(shù)的RA-AF散點貼近AF軸,表明試樣的破壞以拉伸破壞為主。其中頻帶1的RA-AF值分布范圍較廣,頻帶3的信號數(shù)量較少,且不具有明顯的聚類特征,在分析其破壞模式時應(yīng)予剔除。頻帶2和頻帶4的RA-AF值主要聚集在核心密度區(qū)域,且二者僅出現(xiàn)在臨近巖石完全破壞前,因此能夠作為表征巖石的破壞模式的優(yōu)勢頻帶。

圖7 AE特征參數(shù)示意圖Fig.7 Graphical representations of AE characteristic parameters

表3 不同AE主頻帶的RA-AF分布差異Table 3 Differences in RA-AF distribution for different AE dominant bands

圖9 不同AE主頻帶的RA-AF散點分布Fig.9 RA-AF scatter distribution for different AE dominant frequency bands

現(xiàn)階段研究均證實了拉伸破壞對應(yīng)較低的RA值和較高的AF值,剪切破壞對應(yīng)較低的AF值和較高的RA值,但對于RA-AF分布的劃分,尚未形成統(tǒng)一的定論。日本JC MS-Ⅲ B5706混凝土建筑規(guī)范[22]提出,將AF/RA定義為k,將AF/RA>k的信號視為張拉破壞,將AF/RA

通過圖10可以看出,基于RA-AF區(qū)分張拉破壞與剪切破壞的分界線,可能不是一條固定的斜率的直線。這主要是因為巖石作為一種非均質(zhì)的材質(zhì),其損傷特性與外部荷載、原位環(huán)境,礦物顆粒及細觀結(jié)構(gòu)特征等相關(guān),在破壞過程中包含了原生缺陷的壓密,微孔隙及微裂紋的萌生、擴張、延伸,宏觀裂隙的貫通等,表現(xiàn)出階段性特性。因此,在進行巖石AE的破壞模式分析時,不能將全部的RA-AF進行批量處理,更不能將RA-AF分界線簡單定義為一條確定斜率的直線?，F(xiàn)通過結(jié)合AE主頻優(yōu)勢頻帶特征與RA-AF分布密度劃分的分界線,具有一定實際物理意義,增加了AE信號損傷表征能力,在一定程度上促進了深部巖體工程破壞預(yù)測方法由定性分析向定量分析的轉(zhuǎn)變,分析結(jié)果可為深地巖石的破壞模式的研究提供參考。

圖10 優(yōu)勢頻帶的RA-AF散點密度特征Fig.10 RA-AF scatter density characteristics of dominant frequency band

3 結(jié)論

本研究對深部灰?guī)r進行了巴西圓盤劈裂AE試驗,通過對試驗全過程的實測AE波形進行頻-時域聯(lián)合分析,得到以下結(jié)論。

(1)通過實測AE全波形特征,可以將灰?guī)r在巴西圓盤劈裂條件下的損傷過程分為階段Ⅰ～Ⅲ,分別對應(yīng)微損傷萌生階段,裂紋擴展階段,結(jié)構(gòu)破壞階段。在實測AE波形特征中,低幅值主要表現(xiàn)為連續(xù)型特征,高幅值波形則主要表現(xiàn)為突發(fā)型特征。

(2)AE頻譜結(jié)構(gòu)與灰?guī)r損傷量及損傷程度密切相關(guān)。當(dāng)試樣損傷量較小時,AE信號頻譜較窄,表現(xiàn)出單峰結(jié)構(gòu),損傷量較大時,AE信號頻譜變寬,表現(xiàn)出雙峰或多峰結(jié)構(gòu)。

(3)AE主頻具有聚類分布特征,可以分為低頻帶信號與高頻帶信號。低頻帶信號貫穿試驗全過程,主要與微觀層面上的位錯、晶體孿生、晶界面移動相關(guān)。高頻信號僅在巖石臨近完全斷裂前出現(xiàn),主要與宏觀層面上的礦物顆粒、結(jié)構(gòu)面的滑移和分離相關(guān)。

(4)巖石的破壞過程具有階段性特性,基于AE主頻劃分的RA-AF分界線能有效增加AE信號表征灰?guī)r的破裂模式及裂紋分類的能力,為損傷表征定量分析提供重要的理論依據(jù)。