繆海賓

（1.中煤科工集團沈陽研究院有限公司，遼寧 撫順 113122；2.煤礦安全技術(shù)國家重點實驗室，遼寧 撫順 113122）

2013年初，撫順西露天礦南幫整體邊坡發(fā)生大變形，礦坑底部底鼓現(xiàn)象愈發(fā)突顯，坡體上部地表裂縫總長度超過了2.7km，變形體前緣與后緣之間的距離達到1.5km，整個南幫大變形體面積約2.9 km2，初步估算滑體體積超過1億m3，變形體南北向最大水平變形量達150m，垂直向最大變形量達80 m，監(jiān)測數(shù)據(jù)顯示其日變形量在13mm/d以上，最大時達到近200mm/d，變形體整體突顯出了蠕變-大變形特性[1-3]。

軟巖賦存于高陡露天邊坡弱層環(huán)境下，其強度不足以承受露天礦坑底開挖引起的應力集中而產(chǎn)生塑性剪切滑移，開挖應力釋放及調(diào)整后圍巖將形成峰前彈性區(qū)、峰后塑性軟化區(qū)和殘余流動塑性區(qū)，隨著時間不斷增加，處于峰前彈性狀態(tài)的圍巖將由黏彈性變形向黏塑性變形轉(zhuǎn)化，并伴隨黏塑性應變軟化現(xiàn)象，致使峰前彈性區(qū)圍巖發(fā)生非線性蠕變-大變形現(xiàn)象，不僅如此，明顯的非線性黏塑性流變特性也在高陡露天邊坡弱層軟巖峰后塑性區(qū)表現(xiàn)出來，其結(jié)果使得圍巖變形程度增大、變形釋放時間延長，大大地制約了治理工程的建設(shè)工期，進而嚴重影響露天礦的安全生產(chǎn)和經(jīng)濟效益。因此，在研究高陡露天邊坡弱層軟巖蠕變大變形特性的基礎(chǔ)上，提出軟巖蠕變大變形本構(gòu)模型，揭示高陡露天邊坡弱層軟巖大變形非線性蠕變規(guī)律，探尋高陡露天礦邊坡軟巖蠕變-大變形問題的實用理論是露天開采工程中亟需開展的基礎(chǔ)性工作[4-11]。

1 實驗方案

實驗主要設(shè)備為英國GDS高精度軟巖流變儀，該儀器設(shè)備可進行常規(guī)單軸、三軸壓縮和蠕變等實驗。該實驗系統(tǒng)主要包含250kN電機驅(qū)動數(shù)字荷載架、32MPa壓力/體積控制系統(tǒng)、局部應變傳感器、以及多功能測試模塊，可開展軟巖排水或不排水三軸壓縮實驗（UU或CU法）、單軸壓縮及流變實驗和三軸壓縮流變實驗等，具有高測試精度，可以滿足軟巖單軸及三軸流變大變形實驗的要求。試樣上安裝有3個位移傳感器，其中2個測定局部軸向位移，另1個測定徑向位移。

通過設(shè)定不同的軸壓σ1和圍壓σ3，分別開展常規(guī)三軸壓縮試驗、峰前壓縮蠕變實驗、峰后壓縮蠕變實驗。

1）通過常規(guī)三軸壓縮實驗，得到試樣破壞全過程的σ-ε曲線，分析σ-ε曲線演化過程及特征，進而確定軟巖試樣的物理力學參數(shù)。

2）通過對不同圍壓作用下的軟巖開展三軸峰前壓縮蠕變實驗，獲取軟巖蠕變曲線，研究軟巖巖樣軸向、徑向的蠕變應變隨時間的變化規(guī)律，不同圍壓條件下峰前蠕變試驗應力施壓梯度見表1。

表1 不同圍壓條件下峰前蠕變試驗應力施壓梯度Table1 Stress loading gradient under different confining pressures in pre-peak creep test

3）考慮到深部軟巖達到峰后塑性階段時，圍巖的變形擴展仍具有明顯的時效特性。對不同圍壓條件下軟巖峰后應變軟化初段、中段及末段分別進行軟巖峰后壓縮蠕變實驗，獲得軟巖試樣軸向、徑向的蠕變應變隨時間的變化規(guī)律。采取分級加載進行峰后壓縮蠕變試驗，不同圍壓條件下峰后蠕變試驗應力加載梯度見表2。

表2 不同圍壓條件下峰后蠕變試驗應力加載梯度Table2 Stress loading gradient under different confining pressures in post-peak creep test

2 實驗結(jié)果

2.1 軟巖常規(guī)壓縮試驗結(jié)果

采用非線性最小二乘法，得到的弱層軟巖三軸壓縮實驗全過程σ-ε曲線如圖1。

圖1中分別列舉了圍壓為0、5、8、10MPa時的σ-ε曲線，經(jīng)過研究σ-ε曲線的特征可知，當圍壓處于較低狀態(tài)時，弱層軟巖總體呈現(xiàn)為彈脆塑性特征，具有較大的強度衰減幅度，隨著圍壓不斷增大，軟巖呈現(xiàn)出較好的延性特征，當圍壓達到10MPa時，軟巖表現(xiàn)出了典型的彈塑性特征，通過對軟巖的力學參數(shù)進行非線性擬合分析得到的擬合曲線如圖2。由圖2可知，通過實驗數(shù)據(jù)擬合獲得的Mohr-Coulomb強度曲線的擬合精度較高。擬合得到的軟巖常規(guī)力學參數(shù)見表3。現(xiàn)場的實際應用中需考慮巖石剪脹因素影響，由于獲得摩爾庫倫強度參數(shù)的難度較低，因此對弱層軟巖峰后本構(gòu)模型開展了基于摩爾庫倫強度準則的研究。

圖1 軟巖應力-應變曲線Fig.1 Stress-strain curves of soft rock

圖2 軟巖強度擬合曲線Fig.2 Strength fitting curves of soft rock

表3 軟巖強度擬合參數(shù)Table3 Strength fitting parameters of soft rock

2.2 軟巖峰值前壓縮蠕變結(jié)果

弱層軟巖峰值前蠕變曲線如圖3。

通過研究圖3曲線特征可知，在不斷增大圍壓的過程中，弱層軟巖的蠕應變也逐漸增大，由低圍壓狀態(tài)下的2.5 %增大到14.5 %，彰顯出弱層軟巖在高應力作用下具有典型的大變形特征。

圖3 不同圍壓條件下軟巖峰前蠕變及擬合曲線Fig.3 Pre-peak creep and fitting curves of soft rock under different confining pressures

通過分析軟巖的蠕變參數(shù)可得，軟巖的蠕變特性與圍巖的作用密不可分，具有極大相關(guān)性，圍壓增大，軟巖的蠕應變量逐漸遞增，將軟巖的黏滯系數(shù)與圍巖的關(guān)系進行擬合分析，得到的不同圍壓條件下軟巖峰前蠕變參數(shù)變化規(guī)律如圖4。

由圖4可知，圍壓增大，K體（Kelvin體）與M體（Maxwell體）的黏滯系數(shù)隨之卻逐漸減小，其中圍壓對M體的影響最大，當圍壓達到3MPa以上，K體與M體的黏滯系數(shù)受到圍巖的影響基本平穩(wěn)，說明在實際工程中，應力釋放較大的巖體黏滯系數(shù)要遠大于圍巖內(nèi)部的巖體。分別對K體與M體曲線進行指數(shù)擬合研究，得到的軟巖峰前蠕變K體黏滯系數(shù)η1K和M體黏滯系數(shù)η1M隨圍壓σ3變化的指數(shù)方程為：

圖4 不同圍壓條件下軟巖峰前蠕變參數(shù)變化規(guī)律Fig.4 Variation law of pre-peak creep parameters of soft rock under different confining pressures

2.3 軟巖峰值后壓縮蠕變結(jié)果

通過峰值后蠕變實驗獲得的弱層軟巖峰值后蠕變曲線如圖5。

研究圖3和圖5曲線特征可得，峰值前、后的蠕變曲線特征具有較大差異，主要體現(xiàn)在蠕應變方面，通過觀察蠕變曲線可知，軟巖峰值后的蠕應變較小，在峰值后的前期和后期，其蠕應變的范圍在0.1 %～0.6 %，說明軟巖峰值后的蠕變特性極不穩(wěn)定，這與軟巖峰值后的應力路徑和軟巖的破裂特征有關(guān)。

圖5 不同圍壓條件下軟巖峰后蠕變及擬合曲線Fig.5 Post-peak creep and fitting curves of soft rock under different confining pressures

通過分析軟巖的峰值后蠕變參數(shù)，發(fā)現(xiàn)同峰值前的規(guī)律類似，峰值后軟巖的蠕變特性與圍巖的相關(guān)性同樣較大，即隨著圍壓的增大，軟巖的蠕應變量逐漸遞增，將軟巖的黏滯系數(shù)與圍巖的關(guān)系進行擬合分析，得到的不同圍壓條件下軟巖峰后蠕變參數(shù)變化規(guī)律如圖6。

由圖6可知，隨著圍壓的增大，K體與M體的黏滯系數(shù)隨之減小，圍壓對K體的影響最大，而M體的黏滯系數(shù)受到圍巖的影響較小，說明在實際工程中，隨著邊坡巖體的開挖，弱層軟巖會由彈性階段向塑性階段轉(zhuǎn)化，而進入塑性軟化階段后，軟巖仍然具有蠕變特性，導致弱層結(jié)構(gòu)強度持續(xù)弱化的主要原因是軟巖特性發(fā)生了改變。分別對K體與M體曲線進行指數(shù)擬合研究，得到的軟巖峰后蠕變K體黏滯系數(shù)η2K和M體黏滯系數(shù)η2M隨圍壓σ3變化的指數(shù)方程為：

圖6 不同圍壓條件下軟巖峰后蠕變參數(shù)變化規(guī)律Fig.6 Variation law of post-peak creep parameters of soft rock under different confining pressures

3 軟巖BNSS蠕變損傷模型

通過開展軟巖壓縮蠕變實驗并分析結(jié)果，在充分研究Burgers模型的基礎(chǔ)上[12-16]，將非線性M-C塑性元件引入其中，并實現(xiàn)Burgers模型與非線性摩爾庫倫的應變軟化S-S塑性元件的串聯(lián)，建立了BNSS蠕變損傷模型，軟巖BNSS模型如圖7。

圖7 軟巖BNSS模型Fig.7 BNSS model of soft rock

衰減蠕變和穩(wěn)定蠕變可以合理的在Burgers蠕變模型能描述出來，其σ-ε偏量本構(gòu)關(guān)系為：

式中：EK為K體彈性模量；EM為M體彈性模量；ηK為K體黏滯系數(shù)；ηM為M體黏滯系數(shù)；eij為偏應變率；Sij為偏應力張量。

在σ1和σ3作用下，推導出基于Burgers模型的巖石軸向蠕變應變ε（t）為：

式中：GK為K體剪切模量；GM為K體剪切模量；K為體積模量；t為時間。

式中：v為泊松比。

將S-S塑性非線性元件引入其中，假定其服從摩爾庫倫準則，即屈服函數(shù)f=0，則剪切屈服fs和拉伸屈服ft在主軸應力空間的公式分別為：

式中：φ為巖石的內(nèi)摩擦角；c為巖石的黏聚力；σt為抗拉強度；σ1、σ3分別為最大、最小主應力（拉為正）。

假設(shè)c0、φ0分別為巖石初始黏聚力和內(nèi)摩擦角，在蠕變過程中，c0、φ0的蠕變損傷方程應為：

式中：D為損傷變量；t為時間。

將式（9）代入式（8），得到應變軟化S-S塑性非線性元件與時間有關(guān)的剪切屈服fs準則為：

BNSS蠕變損傷模型中總偏應變率eij為：

式中：eijK為K體偏應變；eijM為M體偏應變；eijP為黏塑性偏應變。

對Kelvin體，有：

對Maxwell體，有：

應變軟化S-S非線性體的偏應變率為：

式中：g為勢函數(shù)；evolP為黏塑性體應變；λ為僅在塑性流動階段非零的參數(shù)，通過fs=0確定；σij為應力分量；δij為應力張量；σ1、σ2、σ3分別為最大主應力、中間主應力和最小主應力。

勢函數(shù)g形式如下：

式 中：Nψ為 塑 性 勢 函 數(shù)，Nψ=（1+sinΨ）/（1-sinΨ）；Ψ為巖石剪脹角。

4 結(jié) 語

1）通過軟巖三軸壓縮蠕變實驗，得到的其常規(guī)應力-應變曲線表現(xiàn)出了明顯的彈脆塑性向彈塑性轉(zhuǎn)化的趨勢，且該現(xiàn)象較為明顯。通過擬合得到軟巖峰值前及峰值后殘余強度指標黏聚力和內(nèi)摩擦角。

2）通過分析軟巖的蠕變曲線規(guī)律，軟巖的蠕應變對圍壓的敏感度較高，在低圍壓、高應力作用下具有典型的大變形特征；其峰值后的蠕變特性不穩(wěn)定，易受外界因素的干擾。

3）軟巖峰值前的黏滯系數(shù)受圍壓的影響大，隨著圍壓的增大，K體與M體的黏滯系數(shù)逐漸減小，圍壓達到3MPa后二者趨于平穩(wěn)；圍壓對峰值后的黏滯系數(shù)影響較小，但對K體黏滯系數(shù)影響的要較M體的大，弱層軟巖會由彈性階段向塑性階段轉(zhuǎn)化。

4）通過實驗獲得的軟巖蠕變-大變形試驗規(guī)律及建立的本構(gòu)模型，可為后期蠕變-大變形高陡邊坡破壞機理及防治措施的研究提供實驗依據(jù)及理論支撐。