拉裂-墜落式巖質崩塌失穩過程聲信號特征模擬試驗研究

廖之戀,陳世萬,楊 貴,余 琪,李本云

(1. 貴州大學 資源與環境工程學院,貴州 貴陽 550025; 2. 貴州大學 喀斯特地質資源與環境教育部重點實驗室,貴州 貴陽 550025;3. 貴州省交通規劃勘察設計研究院股份有限公司,貴州 貴陽 550081)

0 引言

我國西南巖溶集中發育,面積約55萬km2,巖溶山地地形切割劇烈,巖體崩塌十分發育。以貴州省為例,據2018年地質災害排查結果,全省范圍內共有3681個崩塌災害隱患點[1];調查發現,超過80%崩塌災害點主體為碳酸鹽巖,巖溶構造發育;而底座多為軟弱碎屑巖,呈典型“上硬下軟”的特征。差異風化強烈,下部軟巖逐步形成凹巖腔最終導致上覆危巖體臨空,在降雨侵蝕作用或地下采礦等工程作用下在危巖體內形成主控裂隙[2-3],主控裂隙經歷起裂、擴展、直至貫通失穩形成崩塌[4-8]。崩塌具有突發、快速、致災性強的特性,極大地威脅人民生命財產安全[9-10]。

巖溶山區巖體崩塌點多面廣、失穩過程迅速,使有效監測預警極具挑戰。聲發射/微震監測精度高,具有連續、實時動態監測的特性[11-12],對巖石內部微破裂演化過程十分敏感。聲發射(AE)技術可以很好地揭示微裂紋的演化規律,由聲發射特征參數進行綜合分析的裂紋擴展和損傷演化的階段性特征與巖石變形破壞過程具有較好的一致性[13]。KAISER[14]最早在金屬材料中發現聲發射現象,隨后國內外對巖石類材料的聲發射特征進行了研究[15-20],并將聲發射技術用于礦山、隧道、邊坡等工程穩定性的監測[21-24]。近年來,一些研究證明了聲發射技術監測邊坡失穩破壞具有潛在優勢[25-26],CODEGLIA等[27]開展了凍融作用下危巖體崩塌的監測,通過巖質邊坡上長周期內振鈴計數的分析,發現地下水及積雪的變化會影響聲發射信號的強弱。DIXON等[28]利用有源波導的聲發射設備監測了英國鐵路軌道邊坡的變形。魏學勇等[29]利用聲發射技術監測了三峽庫區萬州危巖體變形和裂縫擴展。

微震技術與聲發射技術相似,也是利用巖石破壞過程中的振動信號來研究巖體內部的裂縫擴展情況,微震監測所采集的振動頻率較低[30]。微震可有效定位數百米外巖體崩塌崩落產生的信號,同時通過識別信號特征參數可準確識別巖體破裂產生的微震信號。微震監測是分析巖質邊坡穩定性和巖體損傷的有效手段,微震監測技術已開始應用于崩塌災害研究,如:崩塌落石監測定位[31-32]、巖體破裂微震信號特征識別[33-34]、落石運動撞擊解體信號分析[35-36]、落石活動性[37]、崩塌過程反演[38],這些研究證實了微震可有效定位數百米外巖體崩塌崩落產生的信號,同時通過識別信號特征參數可準確識別巖體破裂產生的微震信號。此外,一些研究發現巖體崩塌失穩前也會出現明顯微震前兆信號。SENFAUTE等[39]根據陡崖危巖崩落前15 h內微震信號頻率演化特征將崩塌孕育過程劃分為:裂紋啟裂段、裂紋生長發展段和裂紋貫通破壞段;LIANG等[40]應用巖質邊坡破壞前落石撞擊、彈跳和滾動累計微震能量對巖體崩塌進行預警并成功提前3 min預警一次較大規模崩塌。

總體上看,巖體崩塌監測預警極具挑戰性[41-42]。本文根據巖體崩塌孕育發展力學過程,開展分級蠕變加載和模擬自然漸變2種工況下的拉裂型崩塌相似模擬室內試驗,采用聲發射系統和無線微震傳感器研究崩塌模型失穩過程聲信號特征,研究成果為發展巖體崩塌監測預警方法和技術提供支撐。

1 模型設計及試驗方法

1.1 相似材料試樣配置

1.1.1 脆性相似材料

相似模擬試驗中,相似常數主要有幾何相似常數、重度相似常數和應力相似常數。一般來說,模型的尺寸與原型越相近,則試驗效果越好。本試驗以巖溶山區典型巖體崩塌主巖灰巖為研究對象,主要物理力學參數為單軸抗壓強度σc、楊氏模量E、泊松比μ、密度ρ以及幾何尺寸L。根據相似準則可得到相似常數之間的關系為:

(1)

(2)

CE=CC=Cσc

(3)

Cσ=αL×Cρ

(4)

式中: 下標H、M分別為原型和模型;αL為幾何相似常數;Cσ為應力相似常數;CE為彈性模量相似常數;Cρ為密度相似常數。

根據式(1)～式(4),結合灰巖性能計算選取Cρ=1.7,由于室內試驗儀器及場地限制選取幾何相似常數αL為50,則其他力學相似常數CE,Cc,Cσc均為85,并根據相似比確定灰巖相似材料力學參數范圍,如表1所示。

表1 灰巖物理力學參數及對應相似材料物理力學參數Table 1 Physical and mechanical parameters of limestone and corresponding similar materials

獲得脆性相似材料,確保模型在破壞過程中產生與巖石破壞相似的聲信號是本研究的前提。采用石英砂為骨料,石膏、水泥為膠結材料配制相似材料[43],配制過程中加入3%硼砂作為緩凝劑防止石膏在制樣過程中凝結過快。以骨料與膠結材料比例及膠結成分比例為影響因素,配置脆性相似材料。配比材料中石英砂粒徑40～80目,石膏粉280目,水泥325目。

對不同配比的相似材料試樣進行了單軸壓縮條件下的聲發射試驗如圖1所示。由試驗結果可知,砂膠比5∶1,膠結物中水泥∶石膏為3∶7的材料在力學曲線和聲發射活動特征方面與脆性巖石具有較好的相似性如圖2所示,其表現為:在初始加載孔隙壓密段,孔隙閉合過程使應力-應變曲線呈下凹,聲發射活動較強;在彈性壓縮階段,聲發射活動保持穩定;而在峰前裂隙加速破裂階段,應力-應變曲線斜率降低,彈性模量減小,聲發射出現激增;峰后應力突降,表現出顯著脆性特征。試驗中崩塌主體采用此配比的相似材料。

圖1 部分不同配比試樣聲發射效果圖Fig. 1 Acoustic emission effect diagram of partial samples with different proportions

圖2 灰巖單軸壓縮條件下的聲發射效果圖Fig. 2 Acoustic emission effect diagram of limestone under uniaxial compression condition

1.1.2 軟弱基座相似材料

軟弱基座相似材料需具有遇水強度顯著軟化的特性,通過加入黃豆面親水軟化材料獲得親水軟化特性[44]。軟弱基座相似材料的配比為砂膠比5∶1,即石英砂與水泥加高嶺土比例為5∶1,并添加1%的黃豆面。

1.2 試驗模型與試驗設備

設計巖層結構為上硬下軟型,上部為灰巖,采用材料為脆性相似材料;下部軟弱層為泥巖,采用材料為軟弱基座相似材料,如圖3所示。上覆危巖體通過預制裂紋模擬崩塌主控結構面,下部軟弱層持續注水軟化模擬凹巖腔形成過程。

圖3 試驗示意圖Fig. 3 Diagram of test表2 試驗設備及主要參數設置Table 2 Test equipment and main parameter setting試驗設備型號性能參數及設置高速攝像機千眼狼高速攝像機采集頻率1000 Hz動態應變采集儀DH3818最大采集頻率0.2 Hz測量量程19999 με最高分辨率1 με聲發射監測系統PCI-2E5.40傳感器中心頻率150 kHz聲發射門檻值35 dB峰值定義時間50 μs撞擊定義時間100 μs撞擊閉鎖時間1000 μs微震監測系統新型無線微震監測系統Wi-Fi傳輸方式采樣頻率0.5 kHz

本次試驗所用設備有高速攝像機、動態應變采集儀、無線微震監測系統、聲發射監測系統等。試樣整體尺寸為800 mm×200 mm×800 mm,軟弱基座尺寸為450 mm×200 mm×450 mm。在試樣預制裂紋處粘貼應變片,通過動態應變采集儀可獲得預制裂紋處應變演化過程。采用高速攝像機記錄巖體破裂發展過程,在攝像區域布置密集測點陣列,再將采集到的圖片進行數字圖像相關法(digital image correlation, DIC)獲得破裂過程應變場演化。采用專門無線微震監測系統采集破裂發展過程的低頻聲信號,微震傳感器內置無線傳輸模塊,由無線路由器傳輸數據。采用美國PAC物理聲學公司聲發射監測系統(PCI-2E5.40)采集高頻聲信號。試驗設備及主要參數如表2所示。

1.3 試驗方法

1.3.1 預制裂紋長度設計

研究表明拉裂型巖體崩塌主控結構面演化過程為Ⅰ-Ⅱ型復合裂紋問題[45],主控結構面受力情況主要由彎矩作用和垂直于主控結構面的拉應力Ⅰ型裂紋(圖4(a))和平行于主控結構面的剪應力Ⅱ型裂紋(圖4(b))兩部分組成。

圖4 拉裂型崩塌簡化斷裂力學模型Fig. 4 Simplified fracture mechanics model for crack rock collapse

Ⅰ型裂紋應力強度因子計算公式為:

(5)

(6)

(7)

式中:F為常量系數;σmax為主控結構面受到的最大拉應力;β為貫通率;l為凹巖腔深度;H為上覆危巖體高度;γ為危巖體的重度。

Ⅱ型裂紋應力強度因子計算公式為:

(8)

(9)

復合裂紋斷裂問題可通過最大周向正應力理論的復合應力強度因子Ke來判斷崩塌的穩定性。

(10)

式中θ為破裂角,可由式(11)計算得到:

(11)

式中,k=(KⅠ/KⅡ)2, (KⅡ≠0)。

將計算得到的Ⅰ型和Ⅱ型應力強度因子代入式(11),即可得到破裂角θ,再由θ即可得到復合應力強度因子Ke。復合斷裂韌度判別為:Ke≥KⅠC,即當Ke-KⅠC≥0時,將發生失穩破壞。

本次試驗中巖體崩塌模型各參數為:H=0.35 m,γ=1800 kN/m3,l=0.45 m,β取0.1,0.2,……,0.9。本試驗所涉及的貫通率β為預制裂隙長度與上覆危巖體高度之比[46],經計算貫通率為90%時,在自重條件下可發生破壞。本文試驗中,分級蠕變加載條件下貫通率取80%,模擬自然漸變條件下貫通率取90%。

1.3.2 加載方式

巖體崩塌孕育失穩過程受控于在自然漸變荷載作用下主控結構面的發展。自然狀態下主控結構面的受力方式類似于分級蠕變加載[47],每一次的擾動后(如洪水沖刷基座等)崩塌應力場重分布,達到新的穩定狀態并持續一段時間,擾動和穩定時間段可視為分級加載蠕變;同時分析拉裂巖質崩塌發育力學過程,部分崩塌發育過程可視為基座逐步軟化引起主控裂紋處荷載發生改變,直至崩塌失穩。本文通過概化兩類崩塌發育力學演化特征,開展了分級蠕變加載和基座軟化逐步演化2種類型的室內相似模擬試驗。

1)分級蠕變加載

通過在試樣斷裂自由端上部進行分級蠕變加載,模擬崩塌發育多次外荷載改變引起崩塌失穩破壞過程?；诰€彈性斷裂力學理論,根據試樣I型裂紋斷裂韌性,設計分級加載的荷載為500 N和1000 N,試樣編號為CL-5。

2)自然漸變軟化荷載

為模擬基座漸變軟化引起拉裂型巖質崩塌失穩破壞過程,開展僅通過注水軟化基座,模擬崩塌發育過程。在基座不斷軟化,剛度降低的自然漸變過程中,研究漸進破壞全過程聲發射和微震信號演化規律,試樣編號為DL-2。

2 試驗結果分析

2.1 應變及應變場演化

2.1.1 裂縫應變分析

分析動態應變采集儀獲得預制裂紋處應變演化過程。圖5(a)為分級加載過程中裂縫產生的應變,圖5(b)為圖5(a)中紅色圈定區域的局部放大圖。按分級加載時間可將全過程分為3個階段,前2個階段(0～4500s)為加載階段,應變響應極其微弱,應變增長速度緩慢;第三階段初(第5502.04 s)應變值僅為0.035%,直至破壞前0.04 s時應變出現突增,由0.035%升至5.9%。

圖5 各試樣應變-時間曲線圖Fig. 5 Strain-time curve of each specimen

圖5(c)為自然漸變軟化過程中裂縫產生的應變,圖5(d)為圖5(c)中紅色圈定區域的局部放大圖。由圖5(c)可知,自然漸變軟化荷載下DL-2試樣前期應變響應微弱;直至破壞前0.05 s應變瞬間增大而發生快速失穩破壞?？梢?模擬自然漸變條件與分級蠕變加載條件的變形監測結果類似,均未出現明顯失穩前兆信息。

2.1.2 數字散斑應變場分析

CL-5試樣裂紋從起裂擴展到貫穿失穩全過程,應變場有明顯改變并集中在裂紋處,如圖6所示。以0.001 s為時間間隔,從裂紋起裂到貫穿失穩全過程最大應變分別為1%和8%,崩塌主控裂紋尖端變形增長率迅速增大,由起裂擴展至貫通僅用了0.03 s。DL-2試樣由起裂到貫通僅0.00612 s。2個試樣的破壞速度極快,應變響應在瞬間完成,變形監測難以提供可靠破裂前兆信息。

2.2 高頻信號時頻特征演化規律研究

本試驗中聲發射采集的信號頻率均為高頻信號,從時域和頻域兩方面分析聲發射采集的高頻信號。

2.2.1 聲發射信號時域參數分析

聲發射事件的幅值變化直觀反映了巖樣在加載過程中巖石所釋放能量變化;累計撞擊數為超過門檻并使某一通道獲取數據的任何信號的累計,反映聲發射活動的總量和頻度[48]。

圖7(a)為分級蠕變加載下聲發射信號分布,CL-5試樣在初始注水軟化階段幾乎沒有信號產生,基本處于穩定階段。在第一級荷載下導致主控結構面附近出現應力集中,產生聲發射信號,此階段聲發射幅值較小,集中分布在30～55 dB之間。經第二次加載之后,信號持續增加,但整體幅值相比于第一次加載并未明顯增大;直至破壞前25.6 s,累計撞擊數曲線出現突增,大于60 dB的信號數量劇增?？梢妿r體崩塌孕育過程中聲發射信號響應較為靈敏,距離破壞前較長時間出現聲發射信號突增。

圖7 各試樣聲發射幅值和累計撞擊數演化Fig. 7 Evolution of acoustic emission amplitude and cumulative impact number of each specimen

圖7(b)為模擬自然漸變軟化荷載下聲發射信號分布,根據試驗過程可將全過程分為凹巖腔軟化階段Ⅰ、加速損傷階段Ⅱ、失穩破壞階段Ⅲ。在凹巖腔軟化階段(0～4000 s)末期集中產生較多信號,信號幅值較低,主要為小尺度破裂。隨著不斷注水軟化,在加速損傷階段(4000～5292.1 s)產生較多高幅值信號,于5292.1 s時累計撞擊數曲線出現陡增。相比分級加載,模擬自然漸變條件下的聲發射信號表現出明顯的階段性特征,可用于確定所處發展階段。

2)r值演化規律

r值為聲發射累計撞擊數與聲發射累計能量的比值,用來反映試件加載過程中內部能量的集中度[49-51]。r值持續減小并保持較低的數值,表明聲發射數量少且能量高,處于大破裂孕育階段,其計算公式為:

r=∑N∕∑E

(12)

式中: ∑N為某相同時間間隔的AE累計撞擊數;∑E為某相同時間間隔的AE累計能量。

圖8(a)為分級蠕變加載CL-5試樣的r值演化過程,圖8(b)為圖8(a)中紅色圈定區域的局部放大圖。在每級加載之后均出現降低,表明荷載作用對試樣造成了損傷,產生小尺度破裂。破壞前25.6 s時r值出現突增,突增后緊隨突降,顯示破裂從小破裂急劇轉為大破裂占主導;對比聲發射撞擊數,r值在破裂前出現急劇改變,顯示r值是更為穩健的前兆參數。

圖8 各試樣r值隨時間變化圖Fig. 8 Plot of r value of each specimen with time

圖8(c)為模擬自然漸變荷載下DL-2試樣r值演化過程,圖8(d)為圖8(c)中紅色圈定區域的局部放大圖。在凹巖腔軟化階段的聲發射信號均為小能量信號。進入加速損傷階段r值突增后開始降低,表明在此階段開始產生大尺度破裂。進入失穩破壞階段r值又出現一次突增,之后r值持續降低(距破壞前6.4 s),并持續保持至試樣破壞。

3)b值演化過程

例 ：原文：“yet,as Sir Pitt Crawley would say with perfect justice in his elegant way,‘Rotten!be hanged-it produces me a good fifteen hundred a year.”（Thackeray，2001，p.57）

b值是表征試樣宏微觀裂隙擴展尺度的有效參數,其動態變化特征具有直接的物理意義,揭示了巖石內部應力的變化及裂紋擴展規律[52-53]。計算巖石聲發射b值的公式為[54-55]:

lgN=a-b(A∕20)

(13)

式中:b為表征聲發射活動水平的物理量;A為聲發射幅值;N為震級間隔下的聲發射統計累積頻次;a為常數。

圖9(a)是分級蠕變加載試驗過程中b值演化過程,圖9(b)為圖9(a)中紅色圈定區域的局部放大圖。隨著載荷的逐漸增大,聲發射小事件所占的比例較多b值有所增加,說明此時試樣內部以小尺度微破裂為主,巖石處于相對穩定階段。在試樣破壞前,巖石內部大尺度微裂紋所占比例增加并出現群集現象,微裂紋空間分布開始從無序分布向有序分布轉變,b值開始較快速下降(距破壞前22.25 s)。圖9(c)是模擬自然漸變加載試驗過程中b值演化過程,圖9(d)為圖9(c)中紅色圈定區域的局部放大圖。隨著破裂的發展,DL-2試樣b值同樣逐漸增大,進入加速損傷階段,b值開始較快速下降(距破壞前1.52 s)。

圖9 各試樣b值演化規律Fig. 9 Evolution pattern of b value for each specimen

2.2.2 聲發射信號頻域參數分析

1)聲發射信號主頻

圖10 (a)～(c)為分級蠕變加載試樣破壞過程中高頻聲信號的主頻演化,圖10(c)為圖10(b)中黑色圈定區域的局部放大圖。由圖10 (a)～(c)可知,首次加載后分級蠕變加載試樣主頻主要分布在400 kHz左右且以高頻低幅值為主。第二次加載后(4500～5494.5 s)持續產生低頻大幅值信號,數秒后信號數量突增且此時產生的信號頻段較寬,分布在0～500 kHz之間。直至破壞前25.55 s,信號由高頻向低頻轉化,大尺度破裂產生。由此可見,聲發射主頻帶變寬,由高頻向低頻迅速轉換是巖石破裂可靠的前兆信息[56]。

圖10 主頻及幅值隨時間變化規律圖Fig. 10 Main frequency and amplitude variation pattern with time

圖10(d)～(f)為模擬自然漸變荷載試樣破壞過程中高頻聲信號的主頻演化,圖10(f)為圖10(e)中黑色圈定區域的局部放大圖。由圖10(d)～(f)可知,在凹巖腔軟化階段(0～4000 s)產生的信號多為低頻高幅值;在加速損傷階段(4000～5292.4 s)信號持續呈現低頻特征。臨近破壞時的信號在20～150 kHz之間均有分布且在距破壞前12.08 s出現更多低頻信號。

2)基于AF/RA破裂類型判別

圖11 基于RA與AF的破裂機制判定Fig. 11 Crack mechanism analysis based on RA and AF

在聲發射波形參數指標中,RA(上升時間/振幅)與AF(平均頻率)值可用于對破裂機制進行定性分析。一般來說,低AF值、高RA值表示剪切裂紋的產生或發育,高AF值、低RA值則表示張拉裂紋的產生或發育[57-60],如圖11所示。

圖12(a)～(d)為分級蠕變加載不同時段內的RA與AF值分布,CL-5試樣經第一次加載后,微裂紋總量稀疏,RA值分布在0～10 kHz之間,AF值主要分布30～80 kHz之間。經第二次加載后,RA逐漸增大,AF值逐漸減小,表明剪切裂紋開始增加。臨近破壞時,RA值增大,AF值降低,剪切破裂成分迅速增加,AF-RA值呈現明顯條狀分布。

圖12 不同時段內的RA與AF值分布 Fig. 12 Distribution of RA and AF in different periods

圖12(e)～(h)為模擬自然漸變荷載下不同時段內的RA與AF值分布,DL-2試樣前期RA分布在0～50 ms/V之間,AF值分布在10～60 kHz之間。在失穩破壞階段RA值增大,AF值減小,AF-RA分布集中在剪切破裂區域。

2.3 低頻信號分析

2.3.1 微震信號全局波形分析

DL-2試樣的全過程波形如圖13所示,將微震信號幅值做歸一化處理之后選定閾值y=0.55,分析試驗全過程微震信號幅值演化特征。圖13(b)為DL-2試樣全過程微震信號累計過閾值數出現3次明顯的突增,相對于變形監測,微震信號對加載和破壞十分靈敏;突增后信號保持平穩,表明破裂呈階段發育;而破壞前出現信號突增,裂紋加劇破壞。

圖13 微震信號累計過閾值數Fig. 13 Cumulative number of microseismic signals over the threshold

2.3.2 微震信號主頻分析

微震監測系統采集到的信號為連續波形,采用滑動時窗法(窗口大小為500個/s)截取微震信號,計算每段信號的主頻,DL-2試樣的微震信號主頻演化如圖14所示。

圖14 自然漸變荷載主頻演化規律Fig. 14 Main frequency evolution pattern of natural gradual loading

由圖14可知,DL-2試樣主頻具有呈連續密集條帶狀分布規律,DL-2試樣主頻可分為:2～18 Hz,44～92 Hz,144～446 Hz。按主頻分布特征可將微震信號分為A(高頻)、B(中頻)、C(低頻)三類信號,A類信號數量遠多于B類和C類信號。由圖15可知,三類信號的時頻域波形均有較為明顯的差異。根據時域波形可將A類和B類信號定義為連續信號,C類信號為突變信號。按頻域進行分類,A類信號為單峰結構,B類和C類信號為多峰結構。

2.3.3 主頻變異系數

采用滑動法計算主頻變異系數,主頻變異系數為主頻標準差和均值的比值,計算公式為:

(14)

(15)

DL-2試樣三類信號的主頻變異系數計算結果如圖16所示。由圖16(a)可知,DL-2試樣A類信號的主頻變異系數較為密集,總體上都呈現先升后降的趨勢,在失穩前較長時間內均出現降低趨勢。由圖16(b)可知,B類信號變化趨勢表現為波動持續上升。C類信號過少,不宜用于分析趨勢。綜合分析,A類信號失穩前均出現顯著降低且提前時間較長,可通過高頻的A類信號主頻變異系數的突降來提取微震信號前兆點。

3 討論

通過對分級蠕變和模擬自然漸變2種工況下的拉裂型崩塌失穩過程高頻段和低頻段聲信號進行了系統的分析,研究了模型試驗崩塌孕育全過程聲信號特征(累計撞擊數、r值、主頻、主頻變異系數、AF/RA分布、b值)。試驗結果表明,相對于傳統的位移監測,采用聲發射系統和無線微震傳感器監測裂隙發展過程聲信號表現出明顯的前兆信息,表3總結了各類參數在破裂失穩前發生突變時可提供的預警時間。根據變形監測得到的預警時間不到0.1 s,難以達到預警效果;而各類聲信號參數可以提供數秒至數十秒的預警時間,對于更大規模的巖體崩塌,可提供更長的預警時間。

分級蠕變加載下的預警時間比模擬自然漸變的預警時間長,分析其原因應為分級蠕變加載較模擬自然漸變擾動大,分級蠕變加載加速危巖體損傷。統計參數(b值)的提前量少于物理參數(累計撞擊數、r值、主頻、頻率質心、主頻變異系數、AF/RA分布),預警時間較為保守,分析其原因應為統計參數計算的是某一時段內的數據,而物理參數則是精確的某一刻數據的變化;而統計參數為無量綱的表征破裂本質轉換的參數,更具可靠性。

應該指出,本文研究中未考慮與時間相關的加載速率(分級加載中的時間間隔、自然漸變加載中基座軟化速率)等因素,研究的量綱中未包含時間,因此不能通過本次試驗結果獲得真實巖體崩塌微震預警時間,在后續研究中將進一步考慮荷載施加速率與自然巖體崩塌的相似性。

4 結論

根據巖體崩塌孕育發展力學過程,開展分級蠕變和自然漸變2種工況下的拉裂型崩塌相似模擬試驗,采用聲發射系統和無線微震傳感器研究崩塌模型失穩過程聲信號特征(累計撞擊數、r值、主頻、主頻變異系數、AF/RA分布、b值),得到以下結論:

1)位移和應變場均僅在主破裂發展瞬間產生突變,難以提供指示失穩破壞的前兆信息。

2)根據2種工況下相似模擬試驗聲發射和微震信號的特征,崩塌發展過程可分為漸變階段、主控面加速損傷階段和失穩破壞階段。模擬自然漸變荷載下主控裂紋貫穿失穩更為突然。

3)分析2種工況下巖體崩塌相似試驗的高頻聲發射信號時域和頻域參數特征,可為巖體崩塌失穩提供有效的失穩前兆指示。聲發射時域表現為累計撞擊數的突增、r值的持續降低、b值突降。頻域為主頻帶變寬,臨近破壞時信號由高頻迅速向低頻轉化;RA-AF在加速損傷階段和失穩破壞階段出現轉化,可有效指示破壞過程中破裂主控機制的轉化。

4)采用滑動窗口法分析了試驗全過程的微震波形信號,根據主頻特征可清晰區分高頻、中頻、低頻3類信號,其中低頻信號在加速損傷階段出現突降,可通過低頻信號的變異系數來確定微震信號的前兆信息。