基于能量特征的不同深度煤巖脆塑性評(píng)價(jià)方法

張嵐斌， 劉 洋， 賈哲強(qiáng)， 張朝鵬， 艾 婷， 張澤天

（1.四川大學(xué)a.水利水電學(xué)院；b.深地科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室，成都610065；2.中鐵二院工程集團(tuán)有限責(zé)任公司，成都610031）

0 引 言

隨著地球淺部資源日趨枯竭，礦產(chǎn)資源的開發(fā)深度不斷增加，深部資源開采過程的可行性與工程安全性等關(guān)鍵技術(shù)問題亟需相應(yīng)的理論來支撐保障。研究表明，深部巖體的物性與賦存環(huán)境不同于淺部巖體，適用于淺部巖體變形破壞的傳統(tǒng)巖石力學(xué)理論應(yīng)用于深部巖體時(shí)存在爭議［1-2］，因此需要結(jié)合深部巖體物性與賦存環(huán)境，進(jìn)一步探索深部巖體變形破壞過程中的力學(xué)行為響應(yīng)規(guī)律，發(fā)展適用于深部巖體力學(xué)行為的相關(guān)理論。

對(duì)深部巖體力學(xué)性質(zhì)的研究表明，開展巖石脆性分析對(duì)于深部礦產(chǎn)資源開采具有重要意義［3］。對(duì)于深部煤炭開采而言，煤巖脆性是煤巖力學(xué)性質(zhì)的重要指標(biāo)，也是預(yù)裂爆破抽采瓦斯效果的重要參考，對(duì)預(yù)測(cè)巖爆及巷道圍巖穩(wěn)定性分析具有重要作用。目前學(xué)者們從不同角度對(duì)脆性進(jìn)行了定義［4-5］，并提出了多種表征巖石脆性的評(píng)價(jià)方法［6-7］。其中，有學(xué)者從能量角度進(jìn)行了分析。由熱力學(xué)定律可知，能量轉(zhuǎn)化是物質(zhì)物理過程的本質(zhì)特征，物質(zhì)破壞是能量驅(qū)動(dòng)下的一種狀態(tài)失穩(wěn)現(xiàn)象［8］。研究表明，巖石的脆塑性破壞與能量積聚、釋放現(xiàn)象密切相關(guān)。因此，從能量角度分析不同賦存深度巖石的脆塑性特征能夠更好地揭示巖石脆延性破壞的力學(xué)特征、煤巖彈塑性強(qiáng)化和彈脆性弱化現(xiàn)象。根據(jù)巖石變形破壞全應(yīng)力-應(yīng)變曲線所反映的應(yīng)變能演化規(guī)律，Baron［9］將脆性指數(shù)定義為可釋放彈性能與全過程輸入的總能量的比值；Kidybinski［10］則將脆性指數(shù)定義為峰前耗散能與峰值應(yīng)力處儲(chǔ)存的可釋放彈性能的比值；Tarasov等［11］將峰值應(yīng)力后破裂能量與釋放的彈性能的比值作為脆性評(píng)價(jià)指數(shù)。上述評(píng)價(jià)指標(biāo)均只考慮了峰前或峰后的能量演化。隨后，張軍等［12］在Tarasov的基礎(chǔ)上，對(duì)評(píng)價(jià)標(biāo)準(zhǔn)進(jìn)行了修正，加入了峰前耗散能的影響，但該評(píng)價(jià)指數(shù)依賴于屈服點(diǎn)的選取。

不同賦存深度巖體力學(xué)行為特征不同，脆塑性特征也隨深度變化呈現(xiàn)規(guī)律性。一方面深部煤巖由于高賦存應(yīng)力的影響，峰前的塑性特征與峰后延性特征均會(huì)增強(qiáng)，因此，針對(duì)于不同賦存深度的煤巖有必要同時(shí)考慮巖樣峰前及峰后變形破壞過程中的能量演化規(guī)律來進(jìn)行脆塑性評(píng)估；另一方面不同深度煤巖應(yīng)力-應(yīng)變曲線的屈服段不明顯，依賴于屈服點(diǎn)選取的脆性評(píng)價(jià)方法在分析不同深度煤巖時(shí)適用性較差。因此，為了提高不同賦存深度煤巖脆塑性評(píng)價(jià)指標(biāo)在實(shí)際工程中的適用性，本文根據(jù)能量特征的物理內(nèi)涵，即彈性能反映了巖石的儲(chǔ)能能力，耗散能反映了巖石的塑性變形特征，從峰前和峰后輸入能量的分配入手，對(duì)峰前峰后相應(yīng)的能量特征進(jìn)行處理得到新的脆性度指標(biāo)計(jì)算公式，并結(jié)合不同埋深煤巖的三軸壓縮試驗(yàn)結(jié)果對(duì)不同賦存深度煤巖脆塑性進(jìn)行分析。

1 能量特征的煤巖脆性度指標(biāo)

1.1 特征能量求解

煤巖的應(yīng)力-應(yīng)變曲線反映了巖石內(nèi)部能量狀態(tài)的轉(zhuǎn)變。煤巖試樣在變形破壞全過程中，與外界能量交換貫穿始終，巖石變形破壞過程是峰值應(yīng)力前不斷吸收外界能量，而在峰值應(yīng)力后不斷釋放能量過程［13］。總體來說，是能量輸入、積聚、耗散、釋放相互平衡和動(dòng)態(tài)演化過程。

圖1中對(duì)煤巖內(nèi)部積聚、耗散和釋放能量數(shù)值的計(jì)算進(jìn)行了說明，提出了5個(gè)能量特征，即Uec、Uer、Ud1、Ur、Ud2，并針對(duì)它們的選取意義進(jìn)行分析：①煤巖結(jié)構(gòu)本身具有一定的儲(chǔ)能能力，峰值應(yīng)力處積聚的彈性能Uec（儲(chǔ)能極限）和殘余應(yīng)力處積聚的彈性能Uer（殘余彈性能）綜合反映煤巖峰值強(qiáng)度、變形特征、殘余強(qiáng)度和變形特征。②考慮到煤巖內(nèi)部的微裂隙和缺陷是不斷演化的，能量的耗散作為煤巖變形破壞的基本特征，使得煤巖的強(qiáng)度和儲(chǔ)能能力降低，導(dǎo)致最終煤巖破壞。因此，煤巖能量的耗散影響煤巖的內(nèi)部損傷發(fā)育及強(qiáng)度劣化，煤樣峰值應(yīng)力處的能量耗散Ud1為此處外界輸入的總能量與儲(chǔ)存的彈性能的差值，綜合反映了煤巖在破壞前結(jié)構(gòu)的損傷劣化程度。③在煤樣達(dá)到峰值應(yīng)力σtc后，煤樣儲(chǔ)存的應(yīng)變能迅速釋放，軸向應(yīng)力迅速降低，一些能量轉(zhuǎn)化為了其他形式的能量，如裂紋擴(kuò)展的表面能、結(jié)構(gòu)的震動(dòng)、煤巖碎片的動(dòng)能等。這一過程中，煤樣裂隙貫通，結(jié)構(gòu)趨于破碎，如果繼續(xù)加載，由于煤巖破碎后沿?cái)嗔驯砻婊瑒?dòng)，最后逐漸穩(wěn)定到最小值，即為殘余應(yīng)力σr。峰后的破壞模式取決于峰后能量的釋放和耗散特征。對(duì)于脆性破壞，峰后外力做功幾乎為零，以峰前積聚能量的釋放為主，因此，峰值應(yīng)力后煤巖所釋放的彈性能Ur反映了煤巖峰后的脆性破壞特征；對(duì)于延性破壞，儲(chǔ)存的彈性能逐漸釋放，同時(shí)外力繼續(xù)做功，持續(xù)產(chǎn)生能量耗散，應(yīng)力跌落過程外力作用產(chǎn)生的耗散能Ud2反映了外界驅(qū)動(dòng)煤巖峰后產(chǎn)生持續(xù)塑性破壞額外輸入的能量。

圖1 特征能量計(jì)算方法［14］

煤巖變形破壞過程中積聚的可釋放彈性能在三軸條件下表示為［14］：

式中：E0（MPa）為三軸應(yīng)力狀態(tài)下的初始彈性模量；ν為泊松比；σ1（MPa）為施加在試樣上的軸向應(yīng)力；σ2（MPa）、σ3（MPa）為施加在試樣上的圍壓；ε1e為軸向彈性應(yīng)變；ε2e、ε3e為環(huán)向彈性應(yīng)變。其中，煤巖在三軸應(yīng)力狀態(tài)下的彈性模量與泊松比根據(jù)推薦方法計(jì)算［15］：

峰值應(yīng)力處積聚的彈性能Uec和殘余應(yīng)力處積聚的Uer可表示為：

外力做功輸入的能量可以表示為：

在煤樣峰值應(yīng)力處的能量耗散Ud1為此處外界輸入的總能量與儲(chǔ)存的彈性能的差值，即：

式中：Wc為峰值應(yīng)力前外力向煤巖試樣輸入的能量；Uec為峰值應(yīng)力處煤巖試樣儲(chǔ)存的可釋放彈性能；εtc和ε3tc分別為峰值應(yīng)力對(duì)應(yīng)的軸向和環(huán)向應(yīng)變。

應(yīng)力跌落過程產(chǎn)生的耗散能Ud2可表示為

式中：εr和ε3r分別為σr對(duì)應(yīng)的軸向應(yīng)變和環(huán)向應(yīng)變。

峰值應(yīng)力后煤巖所釋放的彈性能Ur可表示為

1.2 脆性度指標(biāo)的求解

在峰前階段，外界能量輸入一部分作為彈性能儲(chǔ)存，峰值應(yīng)力處即Uec，而同時(shí)產(chǎn)生一部分能量耗散，峰值應(yīng)力處即Ud1。對(duì)于理想脆性巖石來說，Ud1為零，外界輸入能量全部轉(zhuǎn)化為彈性能儲(chǔ)存起來，因此定義峰前脆性度指標(biāo)如下［10］：

在峰后階段，驅(qū)動(dòng)煤樣破壞過程除了積聚在煤樣內(nèi)部的彈性能，也有外界能量的輸入。在該過程中，耗散的能量即能量輸入用于塑性做功而額外消耗的能量為Ud2，而應(yīng)力到達(dá)殘余應(yīng)力時(shí)釋放的彈性能為Ur，對(duì)于理想脆性巖石來說，Ud2為零，巖石的峰后破壞具有自我維持的特點(diǎn)，僅靠彈性能的釋放來驅(qū)動(dòng)，因此定義峰后脆性度指標(biāo)如下［11］：

可以看出，Bpre和Bpost均為0～+∞，且當(dāng)Bpre為1時(shí)，表明峰前輸入的能量一半積聚為彈性能，另一半轉(zhuǎn)化為耗散能；同理，當(dāng)Bpost為1時(shí)，表明峰后破壞的能量驅(qū)動(dòng)，一半依賴于彈性能的釋放，另一半依賴于外界能量的輸入。Bpre=1與Bpost=1可分別作為峰前、峰后脆塑性的轉(zhuǎn)變點(diǎn)，但是由于這樣得到的脆性評(píng)價(jià)方法中脆性和塑性的范圍不對(duì)等，脆性范圍為0～1，而塑性范圍則為1～+∞，需要進(jìn)行范圍的統(tǒng)一，進(jìn)行反正切變換有：

兩種模型脆塑性范圍劃分相同，可將兩者綜合考慮，發(fā)展能量型綜合脆性評(píng)價(jià)指標(biāo)如下：

圖2給出了脆性指標(biāo)的變化情況，這樣的指標(biāo)僅依賴峰前和峰后能量特征的獲取，更適合不同賦存深度煤巖的脆性分析，并可根據(jù)不同的目的有側(cè)重性地參考峰前、峰后或綜合指標(biāo)，使其在煤巖脆性分析中更具靈活性和適應(yīng)性。

圖2 考慮峰前峰后能量特征的脆性指標(biāo)變化

2 新脆性評(píng)價(jià)方法的試驗(yàn)分析

本文新脆性評(píng)價(jià)方法基于Zhang等［16］的力學(xué)試驗(yàn)數(shù)據(jù)進(jìn)行計(jì)算分析：實(shí)驗(yàn)所用煤樣采自于河南省平煤礦區(qū)己15煤層300、600、700、850和1 050 m 5個(gè)不同深度采面，施加相應(yīng)深度圍壓σ3=σH（實(shí)測(cè)水平最大應(yīng)力）開展三軸壓縮試驗(yàn)。該試驗(yàn)綜合考慮不同賦存深度煤巖的物性特征和賦存應(yīng)力環(huán)境特征，從實(shí)驗(yàn)上獲得了不同賦存深度煤巖全應(yīng)力應(yīng)變曲線，如圖3所示。

圖3 不同賦存深度煤巖典型應(yīng)力應(yīng)變曲線［16］

2.1 不同賦存深度煤樣的能量特征分析

根據(jù)上述特征能量計(jì)算方法，進(jìn)行煤樣三軸壓縮試驗(yàn)自初始狀態(tài)（未施加圍壓）至試樣破壞的全過程能量特征計(jì)算，得到不同賦存深度煤巖變形破壞過程中的能量特征。如圖4（a）所示，煤巖破壞前產(chǎn)生的耗散能隨賦存深度增加呈增長趨勢(shì)，其平均值由300 m的59 kJ/m3增長至1 050 m的198 kJ/m3，表明隨著賦存深度增加，煤巖破壞前產(chǎn)生了更多的損傷，用于裂紋發(fā)育、擴(kuò)展以及基質(zhì)塑性變形的能量更多。如圖4（b）所示，煤巖破壞后額外產(chǎn)生的耗散能也同樣呈現(xiàn)隨賦存深度增加而增長的趨勢(shì)，其平均值由300 m的222 kJ/m3增長至1 050 m的546 kJ/m3，表明在煤樣破壞后，仍有較多能量用于煤樣進(jìn)一步的損傷和持續(xù)的塑性變形；且隨著賦存深度增加，峰后應(yīng)變軟化現(xiàn)象呈弱化趨勢(shì)，后續(xù)能量的輸入持續(xù)轉(zhuǎn)化為耗散能驅(qū)動(dòng)破壞，且呈現(xiàn)出更強(qiáng)的延性特征。如圖4（c）所示，煤樣在峰值應(yīng)力處儲(chǔ)存的可釋放彈性能隨賦存深度增加，呈現(xiàn)先增加后趨于平穩(wěn)的趨勢(shì)，300 m試樣平均值僅407 kJ/m3，而1 050 m時(shí)試樣平均值達(dá)到了875 kJ/m3。峰前階段，從初始狀態(tài)至相應(yīng)賦存深度的靜水應(yīng)力狀態(tài)時(shí)，高賦存應(yīng)力的煤樣會(huì)儲(chǔ)存較多的彈性能，繼續(xù)增加軸向壓力至峰值應(yīng)力，隨著深度增加，盡管深部煤巖峰前耗散更多能量，剛度下降更多，結(jié)構(gòu)劣化更嚴(yán)重，但由于對(duì)應(yīng)更大的高賦存應(yīng)力環(huán)境影響，很大一部分吸收的能量仍是轉(zhuǎn)換為彈性能儲(chǔ)存起來。如圖4（d）所示，煤樣在峰值應(yīng)力后，儲(chǔ)存的彈性能釋放，應(yīng)力跌落，逐漸穩(wěn)定在殘余應(yīng)力，此時(shí)由于圍壓的限制，煤樣仍具備儲(chǔ)存能量的能力。隨著賦存深度增加，峰值應(yīng)力后殘余的彈性能隨之增加，1 050 m試樣平均值786 kJ/m3，是300 m試樣平均值的3.14倍。煤樣破壞之后，深部煤巖仍儲(chǔ)存了較多的彈性能，有可能繼續(xù)釋放。如圖4（e）所示，煤樣由峰值應(yīng)力跌落至殘余應(yīng)力的過程中，一方面外力繼續(xù)做功，另一方面彈性能大量釋放，驅(qū)動(dòng)煤樣破壞、損傷迅速發(fā)育直至宏觀裂紋的貫通。由此可見，煤巖破壞時(shí)釋放的彈性能隨著賦存深度增加，有略微下降的趨勢(shì)，平均值從300 m的156 kJ/m3下降至1 050 m的89 kJ/m3。

圖4 不同賦存深度煤巖變形破壞的能量特征

2.2 煤巖脆塑性轉(zhuǎn)變的深度特征分析

將不同賦存深度煤巖變形破壞的能量特征代入，得到不同賦存深度煤巖的脆塑性特征，如圖5所示。

對(duì)于峰前脆性度Bpre見圖5（a），不同賦存深度煤巖Bpre均低于1，屬于脆性范疇，并隨賦存深度逐漸增加，用于峰前裂紋擴(kuò)展等損傷和塑性變形的能量遠(yuǎn)小于積聚的彈性能，導(dǎo)致表現(xiàn)出的脆性特征更強(qiáng)，表明就峰前階段而言，不同賦存深度依然以脆性為主，外界輸入的能量更多以可釋放彈性能的形式儲(chǔ)存起來；對(duì)于峰后脆性度Bpost而言，如圖5（b）所示，與峰前情況不同，不同賦存深度煤巖Bpost幾乎均大于1，屬于塑性范疇，且隨著深度逐漸增大，驅(qū)動(dòng)煤樣破壞過程的輸入能與煤樣系統(tǒng)釋放的彈性能的差值增大，表明深部煤巖需要更多的能量輸入來驅(qū)動(dòng)煤樣的破壞，彈性能釋放相對(duì)減少，這樣導(dǎo)致深部煤巖峰后應(yīng)力跌落不明顯，延性特征更為顯著，1 050 m煤巖試樣趨近于理想彈塑性的情況；對(duì)于同時(shí)考慮峰前峰后的能量型脆性指標(biāo)Benergy而言，在綜合峰前峰后脆性特征情況后，如圖5（c）所示，不同賦存深度煤巖的Benergy依然隨深度呈現(xiàn)增長趨勢(shì)，其300、600、700、850 m平均值分別為0.68、0.81、0.94、0.89，均屬于脆性范圍。隨著深度進(jìn)一步增加，煤樣塑性特征更加明顯，至1 050 m時(shí)達(dá)到了1.03。由此可見，平煤脆塑性轉(zhuǎn)變的特征深度大致在700～1 000 m范圍，而該范圍也對(duì)應(yīng)著平煤現(xiàn)場(chǎng)煤炭開采工程的深部開采范圍，深部開采條件下巷道大變形、支護(hù)失效等情況時(shí)有發(fā)生，由此推斷，上述現(xiàn)象與其脆性特征減弱，塑性特征增強(qiáng)關(guān)聯(lián)明顯。綜上，對(duì)不同賦存深度煤巖而言，綜合能量型脆性指標(biāo)更為合理。

圖5 不同賦存深度煤巖脆塑性特征

進(jìn)一步分析，峰前和峰后脆性度出現(xiàn)較大差距的原因在于：平煤礦區(qū)煤樣結(jié)構(gòu)致密性差，微裂紋和孔隙結(jié)構(gòu)分布較多，淺部煤巖在低圍壓作用下，裂隙孔隙閉合，屈服后裂隙迅速擴(kuò)展達(dá)到峰值應(yīng)力，峰前線彈性特征顯著，屈服段不明顯，破壞以脆性斷裂為主，峰前耗散的能量遠(yuǎn)小于積聚的彈性能。而峰后，煤樣受到圍壓限制，裂紋擴(kuò)展貫通受限，更多能量消耗在破壞的煤巖顆粒間的摩擦錯(cuò)動(dòng)，斷裂面的摩擦滑移上，但其仍然具有一定承載能力，能量釋放不徹底，導(dǎo)致峰后更趨于塑性；而對(duì)于深部煤巖，賦存應(yīng)力大，峰前階段積聚了大量的彈性能，由于圍壓作用，煤樣各個(gè)微元變形趨于統(tǒng)一，塑性變形相對(duì)于淺部巖石低圍壓時(shí)更為明顯，但就其自身能量分配而言，耗散的能量依然低于積聚的能量，因此峰前塑性特征相對(duì)較少。而峰后則由于高圍壓限制，試樣依然具備較高承載能力，釋放的彈性能不足以驅(qū)動(dòng)煤樣完全破壞，外界輸入更多能量驅(qū)動(dòng)煤樣產(chǎn)生斷裂面塑性滑移等塑性變形，由此消耗的能量遠(yuǎn)大于應(yīng)力跌落時(shí)釋放的能量，因此峰后呈現(xiàn)出更明顯的塑性特征。而對(duì)于全過程能量型脆性指標(biāo)來說，由于綜合考慮了峰前和峰后的脆塑性特征，得到了更為合理的不同賦存深度煤樣脆塑性演化規(guī)律。

3 結(jié) 論

基于平煤不同賦存深度煤巖三軸壓縮試驗(yàn)結(jié)果，采用新的不同賦存深度煤巖能量脆性度指標(biāo)，揭示了不同賦存深度煤巖能量特征差異，深入分析了煤巖脆延性的深淺差異，主要結(jié)論如下：

（1）煤巖破壞前后產(chǎn)生的耗散能、峰值應(yīng)力處儲(chǔ)存的可釋放彈性能、峰值應(yīng)力后殘余的彈性能均隨深度增加呈現(xiàn)增長趨勢(shì)，而煤樣破壞時(shí)釋放的彈性能隨深度增加呈現(xiàn)減小趨勢(shì)。由此可見，深部煤巖儲(chǔ)能能力更強(qiáng)但損傷更為嚴(yán)重，破壞后釋放的能量相對(duì)更少，塑延性特征顯著，破壞突發(fā)性下降。

（2）提出的能量型新脆性度指標(biāo)Benergy更適合平煤不同賦存深度煤巖的脆塑性分析。

（3）根據(jù)新指標(biāo)，不同賦存深度煤巖峰前以脆性為主，峰后以塑性為主，隨深度增加，煤巖峰前的塑性特征和峰后的延性特征均更為明顯。

（4）綜合考慮峰前、峰后情況的能量型脆性度新指標(biāo)Benergy隨賦存深度增加呈現(xiàn)增長趨勢(shì)，深部煤巖綜合塑性更為明顯，推測(cè)在700～1 000 m范圍內(nèi)發(fā)生脆-塑性特征轉(zhuǎn)化。

本文致力于從能量角度建立不同賦存深度煤巖的脆塑性評(píng)價(jià)指標(biāo)，針對(duì)不同賦存深度煤巖脆塑性評(píng)價(jià)有較好的適用性，但該指標(biāo)對(duì)于其他處于不同賦存深度的巖石是否適用還需進(jìn)一步驗(yàn)證。