孟令超 徐榮超 王安明 黃志全 袁廣祥
(華北水利水電大學(xué), 地球科學(xué)與工程學(xué)院,鄭州 450046,中國(guó))
高應(yīng)力條件下巖體開挖卸荷導(dǎo)致圍巖發(fā)生的脆性破壞現(xiàn)象(劈裂破壞、巖爆等)促使巖石力學(xué)工作者們逐步認(rèn)識(shí)到揭示巖石脆性破壞機(jī)制對(duì)于工程安全穩(wěn)定性研究的重要性。巖石漸進(jìn)的脆性破壞過程包含了裂紋的起裂、擴(kuò)展及貫通直至宏觀破裂面的形成,不同的階段所對(duì)應(yīng)的應(yīng)力水平及變形特性不同,這決定了巖石應(yīng)力-應(yīng)變曲線形態(tài)上的差異性。因而,巖石應(yīng)力門檻值的研究對(duì)于認(rèn)識(shí)巖石脆性破壞過程及機(jī)制具有重要意義。國(guó)外學(xué)者率先系統(tǒng)開展了Lac du Bonnet 花崗巖起裂及損傷強(qiáng)度特性的室內(nèi)及現(xiàn)場(chǎng)試驗(yàn)(Martin, 1993, 1997; Martin et al., 1994, 1997; Cai et al.,2004)。張曉平等(2011)基于丹巴二云英片巖單軸試驗(yàn)結(jié)果,研究了其應(yīng)力門檻值之間的相互關(guān)系,并發(fā)現(xiàn)由于片理面的存在,巖石的裂紋擴(kuò)展過程存在顯著的各向異性特征。梁昌玉等(2012)系統(tǒng)研究了不同加載應(yīng)變率的條件下巖石應(yīng)力門檻值的變化規(guī)律,結(jié)果表明特征應(yīng)力門檻值對(duì)加載應(yīng)變率不敏感。黃達(dá)等(2012)研究了9種加載應(yīng)變率對(duì)粗晶大理巖應(yīng)力-應(yīng)變曲線、特征應(yīng)力的影響特性,發(fā)現(xiàn)起裂及擴(kuò)容應(yīng)力隨加載速率的增大而提高。周輝等(2015)系統(tǒng)總結(jié)了花崗巖及大理巖應(yīng)力門檻值在圍壓影響下的變化規(guī)律,發(fā)現(xiàn)兩種巖石應(yīng)力門檻值隨圍壓呈現(xiàn)近似線性相關(guān)關(guān)系。曹洋兵等(2020)提出一種能全面反映花崗巖變形破壞全過程的脆性評(píng)價(jià)新指標(biāo),分析了不同含水率花崗巖的應(yīng)力門檻值變化規(guī)律。
由于巖石內(nèi)部裂紋的擴(kuò)展行為無法被直接觀測(cè),作為一種無損檢測(cè)手段,聲發(fā)射技術(shù)在巖石力學(xué)試驗(yàn)中,尤其是裂紋擴(kuò)展方面逐步發(fā)揮重要作用。常用的聲發(fā)射分析參數(shù)包括:撞擊、計(jì)數(shù)、能量及b值等。Eberhardt et al.(1998, 1999)利用聲發(fā)射測(cè)試方法,研究了Lac du Bonnet花崗巖在單軸壓縮條件下的起裂強(qiáng)度和損傷強(qiáng)度特性。楊永杰等(2006)采用最大Lyapunov指數(shù)計(jì)算了煤樣單軸試驗(yàn)條件下的能量計(jì)數(shù)率,提出了預(yù)測(cè)試樣破裂時(shí)間的方法。曹樹剛等(2009)對(duì)比分析了煤巖樣單軸及三軸條件下的聲發(fā)射特征。劉寧等(2012)采用裂紋體積應(yīng)變拐點(diǎn)及聲發(fā)射測(cè)試分析了圍巖下大理巖的強(qiáng)度特征。周輝等(2014b)系統(tǒng)探討了聲發(fā)射參數(shù)在應(yīng)力門檻值確定應(yīng)用上的優(yōu)缺點(diǎn)。楊百存等(2017)為研究巖石在不穩(wěn)定擴(kuò)展階段的破裂行為,建立了b值與weibull分布函數(shù)形狀參數(shù)m之間的定量關(guān)系。韓偉歌等(2017)開展了圍壓條件下砂巖的三軸試驗(yàn),對(duì)砂巖試樣破壞過程中裂紋分布情況進(jìn)行了定位研究。王春來等(2018)通過對(duì)比分析不同硬巖的聲發(fā)射b值變化特征。王林均等(2019)對(duì)花崗巖和砂巖兩類巖石在單軸壓縮過程中的聲發(fā)射特性(計(jì)數(shù)、能量、AF值、RA值及b值)進(jìn)行了詳細(xì)分析。
巖石的脆性指的是其抵抗塑性變形的能力(Morley, 1944; Nejati et al.,2014; Akinbinu, 2017),是巖石礦物成分、顆粒膠結(jié)強(qiáng)弱及結(jié)構(gòu)特征的綜合反映,不同脆性的巖石其變形破壞機(jī)制不同,因而聲發(fā)射特性也不同。已有研究,主要針對(duì)不同巖性巖石強(qiáng)度特征定量確定,及其在壓縮破壞過程中聲發(fā)射現(xiàn)象規(guī)律總結(jié),而針對(duì)脆性對(duì)巖石聲發(fā)射特征的影響規(guī)律并沒有系統(tǒng)開展。本研究采用兩種類型硬巖花崗巖及大理巖,開展室內(nèi)單軸壓縮試驗(yàn),同步進(jìn)行聲發(fā)射特征測(cè)試; 結(jié)合試驗(yàn)結(jié)果及試樣破壞現(xiàn)象評(píng)價(jià)兩者的脆性及變形特性的差異; 基于聲發(fā)射測(cè)試結(jié)果,分析脆性對(duì)巖石聲發(fā)射特征參數(shù)(能量及b值)的影響規(guī)律。研究結(jié)果對(duì)于硬巖脆性破壞過程及機(jī)制的認(rèn)識(shí)具有一定指導(dǎo)意義。
試驗(yàn)采用兩種不同巖性的硬巖,花崗巖和大理巖,天然密度分別為2.61 g·cm-3和2.85 g·cm-3,巖樣為直徑50 mm×高度100 mm的標(biāo)準(zhǔn)圓柱試樣,試驗(yàn)前經(jīng)量測(cè),試樣尺寸及精度符合相關(guān)試驗(yàn)規(guī)范的規(guī)定。單軸壓縮試驗(yàn)所用儀器為RMT-150C型試驗(yàn)機(jī),采用軸向位移控制的方式進(jìn)行加載,加載速率設(shè)定為0.12 mm·min-1。聲發(fā)射系統(tǒng)為美國(guó)PAC公司生產(chǎn)的8通道聲發(fā)射監(jiān)測(cè)系統(tǒng),聲發(fā)射探頭型號(hào)為PICO。為有效過濾噪音信號(hào)對(duì)聲發(fā)射結(jié)果的影響,信號(hào)門檻值設(shè)置為40 dB,采樣頻率為1 MHz。試驗(yàn)系統(tǒng)如圖1所示。
圖1 試驗(yàn)系統(tǒng)Fig.1 Testing systema.加載及聲發(fā)射測(cè)試系統(tǒng); b.變形及聲發(fā)射傳感器布置
圖2所示為單軸壓縮下花崗巖和大理巖典型的應(yīng)力-應(yīng)變曲線。從圖中曲線形態(tài)可知,兩種巖石均為典型的脆性巖石,與大理巖相比,花崗巖峰值強(qiáng)度后應(yīng)力跌落速率更快、曲線更陡且峰值強(qiáng)度前未表現(xiàn)出明顯的屈服現(xiàn)象,因而其脆性特性表現(xiàn)得更為顯著。
圖2 花崗巖和大理巖典型的應(yīng)力-應(yīng)變曲線Fig.2 Typical curves of stress with strain for granite and marble
巖石單軸壓縮過程中,依據(jù)應(yīng)力-應(yīng)變曲線形態(tài),可分為5個(gè)不同的階段:裂隙壓密、線彈性、裂紋穩(wěn)定擴(kuò)展、裂紋非穩(wěn)定擴(kuò)展及峰后破壞階段。其中,由裂紋穩(wěn)定擴(kuò)展向裂紋非穩(wěn)定擴(kuò)展轉(zhuǎn)變時(shí),意味著巖石內(nèi)部微裂紋開始相互連接、貫通,產(chǎn)生宏觀的破裂面。該轉(zhuǎn)折點(diǎn)所對(duì)應(yīng)的應(yīng)力值稱為巖石的損傷強(qiáng)度σcd,在該應(yīng)力水平下,即使軸向壓縮荷載不再增大,巖石內(nèi)部裂紋也會(huì)繼續(xù)擴(kuò)展,最終導(dǎo)致巖石宏觀破壞的發(fā)生,該應(yīng)力值也被稱之為巖石的長(zhǎng)期強(qiáng)度。由此可見,損傷強(qiáng)度的分析對(duì)于認(rèn)識(shí)巖石脆性破壞具有重要意義。損傷強(qiáng)度可依據(jù)巖石的體積應(yīng)變曲線或者聲發(fā)射特性曲線來確定,本文采用應(yīng)變測(cè)量法確定,即體積應(yīng)變的拐點(diǎn)所對(duì)應(yīng)的軸向應(yīng)力即為巖石的損傷強(qiáng)度,如圖3所示。圖4所示為兩種巖石損傷強(qiáng)度與峰值強(qiáng)度比值大小,花崗巖的σcd/σp之比介于67.60%~74.51%之間,平均為70.57%,而大理巖的σcd/σp之比介于43.93%~61.46%,平均為52.04%。由此可知,花崗巖的損傷強(qiáng)度與峰值強(qiáng)度應(yīng)力比遠(yuǎn)大于大理巖。顯然這是由巖石本身的礦物成分、膠結(jié)強(qiáng)度及構(gòu)造等因素導(dǎo)致的。
圖3 損傷應(yīng)力σcd 取值方法示意圖Fig.3 Sketch of method to obtain the damage strength σcd
圖4 試樣破壞形態(tài)Fig.4 Sample failure morphology of granite and marblea.花崗巖; b.大理巖
表1所列為花崗巖和大理巖試樣主要力學(xué)參數(shù)統(tǒng)計(jì)值?;◢弾r和大理巖平均抗壓強(qiáng)度分別為198.63 MPa和120.49 MPa,彈性模量分別為55.44 GPa和44.93 GPa,而由于兩者脆性較強(qiáng),其殘余強(qiáng)度都為0 MPa。
表1 試樣主要力學(xué)參數(shù)統(tǒng)計(jì)表Table1 Main mechanical parameters of rock samples
目前,巖石脆性定量評(píng)價(jià)方法眾多,如基于巖石礦物成分、硬度及應(yīng)力-應(yīng)變曲線形態(tài)等。由于巖石應(yīng)力-應(yīng)變曲線具有直觀性且容易獲取,因此,該方法被廣泛采納。周輝等(2014a,2014b)基于巖石峰后應(yīng)力跌落程度及跌落速率提出了一個(gè)表征巖石脆性大小的定量評(píng)價(jià)指標(biāo),該指標(biāo)表達(dá)式為:
(1)
式中:σp和σr分別為峰值強(qiáng)度及殘余強(qiáng)度;kde(DE)表征巖石峰后應(yīng)力跌落速率,其中,D為峰值強(qiáng)度點(diǎn),E為殘余強(qiáng)度起始點(diǎn)。需要說明的是,該脆性指標(biāo)的具體物理意義及計(jì)算方法在原文中有詳述,在此不再贅述。
利用該脆性評(píng)價(jià)指標(biāo),兩種巖石脆性指標(biāo)計(jì)算結(jié)果如表1中所列。根據(jù)表中脆性大小評(píng)價(jià)結(jié)果可知,花崗巖脆性均值為0.517大于大理巖的0.483,這與根據(jù)圖2中兩者應(yīng)力-應(yīng)變曲線形態(tài)所體現(xiàn)出的脆性強(qiáng)弱定性分析結(jié)果一致。
此外,對(duì)比圖4中兩種巖石典型的破壞形態(tài),可以看出:盡管兩者在單軸壓縮條件下均表現(xiàn)出典型的豎向劈裂破壞,但花崗巖試樣內(nèi)裂紋擴(kuò)展更充分、試樣被眾多密集分布的豎向裂紋分裂成板、片狀,而且發(fā)生了明顯的巖塊、巖片彈射現(xiàn)象,這與巖爆現(xiàn)象類似; 而大理巖試樣則是被豎向宏觀破裂面分裂為幾個(gè)巖塊,從試驗(yàn)機(jī)上殘留的粉末狀巖屑可知,其破壞過程中剪切機(jī)制發(fā)揮了作用,除此之外,未產(chǎn)生巖塊、巖片的彈射現(xiàn)象。由此可見,從試樣破壞現(xiàn)象及形態(tài)上來看,花崗巖脆性程度大于大理巖,這進(jìn)一步驗(yàn)證了上述基于應(yīng)力-應(yīng)變曲線形態(tài)和脆性指標(biāo)計(jì)算評(píng)價(jià)結(jié)果的正確性。需要指出的是,本試驗(yàn)采用的花崗巖及大理巖與周輝等(2014b)采用的不是同種巖石,因而對(duì)于兩種巖石的脆性大小評(píng)價(jià)結(jié)果不一致。
巖石變形特性的分析,對(duì)于認(rèn)識(shí)其脆性破壞過程及現(xiàn)場(chǎng)圍巖變形的預(yù)測(cè)具有重要意義。
直接對(duì)比兩種巖石的損傷應(yīng)變(損傷強(qiáng)度對(duì)應(yīng)的軸向應(yīng)變)、峰值應(yīng)變等參數(shù)的絕對(duì)大小,對(duì)于揭示其破壞過程的內(nèi)在機(jī)制意義不明確。實(shí)際上,當(dāng)軸向應(yīng)力超過損傷強(qiáng)度時(shí),巖石的環(huán)向變形增長(zhǎng)速率急劇增加,意味著大量裂紋開始擴(kuò)展與貫通。因此,對(duì)于兩種脆性巖石變形特性的分析,本文以損傷強(qiáng)度至峰值強(qiáng)度即裂紋不穩(wěn)定擴(kuò)展階段為重點(diǎn)。表2給出了兩種巖石的軸向應(yīng)變比(損傷強(qiáng)度所對(duì)應(yīng)的軸向應(yīng)變與峰值應(yīng)變比值)和環(huán)向應(yīng)變比(損傷強(qiáng)度所對(duì)應(yīng)的環(huán)向應(yīng)變與峰值強(qiáng)度對(duì)應(yīng)的環(huán)向應(yīng)變)。以花崗巖試樣1-1為例,當(dāng)壓縮荷載由損傷強(qiáng)度增大至峰值強(qiáng)度時(shí)(即軸向應(yīng)變由74.75%峰值增至峰值大小時(shí)),試樣的軸向變形增幅為25.12%,而環(huán)向變形增幅卻達(dá)到68.03%; 同樣地對(duì)于大理巖試樣2-2,軸向變形增幅為40.82%,而環(huán)向變形增幅卻達(dá)到86.89%。由此可見,脆性巖石環(huán)向變形主要發(fā)生在裂紋不穩(wěn)定擴(kuò)展階段,環(huán)向變形增長(zhǎng)幅值遠(yuǎn)大于軸向變形。
表2 試樣應(yīng)變比統(tǒng)計(jì)表Table2 Strain ration of the sample
進(jìn)一步,將巖樣由損傷強(qiáng)度至峰值強(qiáng)度階段的軸向應(yīng)變?cè)隽繛闄M坐標(biāo),以對(duì)應(yīng)的環(huán)向應(yīng)變?cè)隽繛榭v坐標(biāo)繪圖,如圖5所示。圖5中,擬合直線的斜率的大小定量表征了環(huán)向應(yīng)變相對(duì)于軸向應(yīng)變?cè)鲩L(zhǎng)速率的相對(duì)大小,即斜率越大則在產(chǎn)生相同軸向應(yīng)變?cè)隽康那闆r下,環(huán)向應(yīng)變?cè)鲩L(zhǎng)幅值越大。結(jié)合2.2節(jié)中脆性大小的評(píng)價(jià)結(jié)果(花崗巖大于大理巖)可知,巖石脆性程度越大,在裂紋不穩(wěn)定擴(kuò)展階段,在產(chǎn)生相同的軸向壓縮變形的情況下,環(huán)向變形量越大。上述結(jié)論,可由圖4中兩種巖石的破壞現(xiàn)象得到驗(yàn)證:花崗巖由損傷強(qiáng)度增至峰值時(shí),大量裂紋迅速擴(kuò)展、貫通,同時(shí)裂隙面之間的張開、分離變形增大,導(dǎo)致試樣環(huán)向變形在短時(shí)間內(nèi)急劇增大,而軸向荷載作用下裂隙面切割形成的巖板、巖片發(fā)生劇烈的彈射而脫離試樣母體,進(jìn)而造成試樣整體的承載能力迅速下降,表現(xiàn)在應(yīng)力-應(yīng)變曲線上則是峰后軸向應(yīng)力急劇下降。
圖5 兩種巖石的環(huán)向應(yīng)變?cè)隽颗c軸向應(yīng)變?cè)隽筷P(guān)系Fig.5 Relationship between axial strain increment and circumferential strain increment of two different rocks
有效地利用聲發(fā)射信號(hào)所攜帶的信息,極大提高了人們對(duì)巖石破裂機(jī)制的認(rèn)識(shí)水平。本節(jié)將對(duì)花崗巖和大理巖在單軸壓縮過程中的聲發(fā)射能量演化特性以及聲發(fā)射b值進(jìn)行系統(tǒng)分析。
聲發(fā)射撞擊(AE hit)指的是系統(tǒng)識(shí)別的有效破裂信號(hào),而聲發(fā)射能量(AE energy)則是該信號(hào)強(qiáng)度強(qiáng)弱的定量表征參數(shù)。圖6為花崗巖及大理巖聲發(fā)射能量、累計(jì)能量隨加載時(shí)間的變化曲線(花崗巖以1-1試樣為例,大理巖以試樣2-1為例),為反映單軸加載過程中不同變形階段的聲發(fā)射特性,圖中同時(shí)給出了軸向應(yīng)力隨時(shí)間的變化曲線。由圖中曲線可以看出,在不同的變形階段,巖石的聲發(fā)射能量特征不同:裂隙壓密階段,聲發(fā)射能量始終處于低水平,其大小一般不超過100,這一階段巖石內(nèi)原生裂隙的閉合及裂隙面摩擦?xí)a(chǎn)生一些聲發(fā)射信號(hào),但由于不是真實(shí)的巖石破裂,因而聲發(fā)射撞擊所攜帶的能量值很低; 在線彈性及裂紋穩(wěn)定擴(kuò)展階段,聲發(fā)射能量隨著加載時(shí)間的增加而呈現(xiàn)穩(wěn)定增長(zhǎng)的變化趨勢(shì),累計(jì)能量變化曲線近似線性的增長(zhǎng)規(guī)律也充分說明了這一點(diǎn); 進(jìn)入裂紋不穩(wěn)定擴(kuò)展階段之后,兩種巖石的聲發(fā)射能量表現(xiàn)出不同的變化規(guī)律,對(duì)于花崗巖而言,高能級(jí)的聲發(fā)射信號(hào)急劇增加,說明大尺度的微破裂事件持續(xù)發(fā)生,而大理巖則恰好相反,其聲發(fā)射能量則表現(xiàn)出持續(xù)降低的趨勢(shì),即微破裂事件逐步增多; 峰值強(qiáng)度后,花崗巖有個(gè)別高能級(jí)信號(hào),而大理巖則不同,隨著其軸向應(yīng)力的幾次跌落現(xiàn)象,聲發(fā)射能量也隨之突增,同時(shí)累計(jì)能量曲線也持續(xù)快速增長(zhǎng),最后,隨著軸向應(yīng)力的持續(xù)跌落,聲發(fā)射能量急劇下降,而累計(jì)能量曲線出現(xiàn)停滯增長(zhǎng)的平臺(tái)。
圖6 聲發(fā)射能量及累計(jì)能量隨加載時(shí)間變化曲線Fig.6 Change curves of AE energy and cumulative AE energy with loading timea.花崗巖; b.大理巖
結(jié)合2.2節(jié)中兩種巖石的脆性評(píng)價(jià)結(jié)果,脆性對(duì)巖石聲發(fā)射能量演化特征具有重要影響:強(qiáng)脆性的花崗巖在裂紋不穩(wěn)定擴(kuò)展階段持續(xù)出現(xiàn)高能級(jí)的聲發(fā)射信號(hào),而弱脆性的大理巖則表現(xiàn)出能量持續(xù)降低的變化趨勢(shì)。
聲發(fā)射b值(AE b-value)來源于地震學(xué)的研究,b值的物理意義在于其能夠表征巖石破裂過程中大尺度破裂與小尺度破裂之間的比例關(guān)系,b值越小則大尺度破裂事件所占比例越大。
圖7所示為兩種脆性巖石單軸加載全過程中,縱坐標(biāo)為撞擊數(shù)目,聲發(fā)射撞擊幅值的分布規(guī)律(花崗巖以1-1試樣為例,大理巖以試樣2-1為例)。從圖中可觀察到,整個(gè)巖石破壞過程中大理巖聲發(fā)射撞擊總數(shù)遠(yuǎn)多于花崗巖,而高幅值的撞擊數(shù)目(大于70 dB)卻比花崗巖少得多。由此可知,對(duì)于脆性較弱的大理巖,其變形破壞過中盡管產(chǎn)生了數(shù)目更多的破裂,但多屬于尺度較小的微破裂,而脆性更強(qiáng)的花崗巖則恰好與之相反。
圖7 聲發(fā)射幅值分布Fig.7 Distribution of AE amplitude
以Adb/20為橫坐標(biāo),以log10(N)為縱坐標(biāo),將各幅值區(qū)間內(nèi)的撞擊數(shù)目繪圖,如圖8所示(AdB為聲發(fā)射幅值,N為撞擊數(shù))。由線性擬合結(jié)果可知,花崗巖b值為0.498小于大理巖的b值1.134。結(jié)合聲發(fā)射b值物理意義可知,單軸壓縮下花崗巖破壞過程中大尺度的破裂事件比例較大,而大理巖則是小尺度的微破裂事件所占比例較大,這與上述幅值分布特征分析的結(jié)果相一致。
圖8 聲發(fā)射b值擬合曲線Fig.8 Fitting curves of AE b-value
本文采用花崗巖及大理巖兩種不同巖性的巖石,開展了單軸壓縮及聲發(fā)射測(cè)試試驗(yàn),獲取了兩種巖石的強(qiáng)度及變形特性,并對(duì)其脆性大小進(jìn)行了定量評(píng)價(jià),最后分析了單軸壓縮過程中兩種巖石聲發(fā)射特性,結(jié)合聲發(fā)射b值計(jì)算結(jié)果及其物理意義,對(duì)比了兩種巖石破裂機(jī)制的差異性。主要結(jié)論如下:
(1)本次試驗(yàn)所采用的兩種巖石,花崗巖的σcd/σp之比介于0.676~0.745之間,平均為0.706,而大理巖的σcd/σp之比介于0.4393~0.615,平均為0.52。
(2)基于脆性定量評(píng)價(jià)指標(biāo),結(jié)合試樣宏觀破壞現(xiàn)象及形態(tài),對(duì)花崗巖及大理巖進(jìn)行了脆性程度大小的比較,結(jié)果表明,本次試驗(yàn)采用的花崗巖其脆性大于大理巖。
(3)基于花崗巖及大理巖在裂紋不穩(wěn)定擴(kuò)展階段的軸向及環(huán)向變形特性分析發(fā)現(xiàn):巖石脆性程度越大,在裂紋不穩(wěn)定擴(kuò)展階段,在產(chǎn)生相同的軸向壓縮變形的情況下,環(huán)向變形量越大,這與試樣所表現(xiàn)出的宏觀破壞現(xiàn)象一致。
(4)強(qiáng)脆性的花崗巖在裂紋不穩(wěn)定擴(kuò)展階段持續(xù)出現(xiàn)高能級(jí)的聲發(fā)射信號(hào),而弱脆性的大理巖則表現(xiàn)出能量持續(xù)降低的變化趨勢(shì); 大理巖在整個(gè)變形破壞過程中聲發(fā)射撞擊總數(shù)目遠(yuǎn)多于脆性更強(qiáng)的花崗巖; 花崗巖破壞過程中大尺度的破裂事件所占比例較大,而大理巖則是小尺度的微破裂事件所占比例較大。