尤哲敏，陳建平，李 濤，徐 方

(中國地質(zhì)大學(武漢)工程學院，湖北武漢 430074)

深埋特長隧道穿越中—高地應力區(qū)，而高地應力意味著深埋隧道在硬巖段的開挖過程中，將會出現(xiàn)不同程度的巖爆災害［1］，影響施工安全和開挖進度，因此日益成為地下工程研究的熱點和重點。

巖爆是在高應力條件下，由于地下工程掘進造成圍巖應力重分布，從而引起硬脆性巖體在短時間內(nèi)產(chǎn)生的儲備彈性應變能突然快速釋放，并發(fā)生巖塊剝落、彈射甚至拋擲［2］。

對巖爆的預測目前國內(nèi)外還沒有成熟的理論與方法［3-4］。大多是針對具體工程采用單一或多種判據(jù)，判據(jù)的適用性還有待商榷［5-7］。

大坪山隧道是湖北省谷竹高速公路的一條分離式特長深埋隧道，尚未開挖段存在中—強烈?guī)r爆的可能，因此，對于大坪山隧道硬巖段的巖爆預測亟需研究。本文采用三維有限元回歸分析法，對大坪山隧道硬巖段開挖過程中巖爆發(fā)生的可能性進行預測。

1 工程概況

大坪山隧道位于襄陽市谷城縣紫金鎮(zhèn)與寺坪鎮(zhèn)境內(nèi)，左右洞長分別為8 263 m和8 242 m，隧道軸線呈NE—SW向展布。

隧址區(qū)屬構(gòu)造剝蝕侵蝕低山—低中山區(qū)，地形起伏較大，最大埋深約900 m。該區(qū)出露的地層主要為志留系頁巖、奧陶系生物碎屑灰?guī)r夾少量炭質(zhì)頁巖、粉砂質(zhì)頁巖;下寒武統(tǒng)泥質(zhì)條帶灰?guī)r、炭質(zhì)灰?guī)r局部含頁巖夾層;上震旦統(tǒng)—下寒武統(tǒng)白云巖夾灰?guī)r。大地構(gòu)造上位于揚子準地臺北緣的青峰臺褶束，走向近EW向。青峰斷裂帶與隧道軸線小角度相交，其形成的數(shù)條逆沖斷層、破碎帶及韌性剪切帶見圖1，使隧道圍巖局部變得較破碎，地下水發(fā)育，整體穩(wěn)定性差。

圖1 大坪山隧道軸線地質(zhì)剖面示意

2 工程區(qū)地應力場反演分析

初始地應力的分布是巖爆發(fā)生與否的關鍵性影響因素，地應力的準確評價直接影響到巖爆預測的可靠性［8］。

地應力回歸分析法是利用數(shù)值分析和多元回歸分析，將計算區(qū)域的地應力場視為自重應力場和水平x，y向擠壓和xy面剪切等構(gòu)造應力場的線性疊加，通過分別模擬4個子應力場，組合成計算地應力場［9］。通過三維應力場反演和實測資料擬合，得出隧址區(qū)的初始地應力場特征［10］，為隧道各斷面的巖爆預測提供依據(jù)。

2.1 初始地應力場的現(xiàn)場實測結(jié)果分析

工程區(qū)采用水壓致裂法進行了2個鉆孔的地應力測量，地應力分布具有如下特點［11］:①最大水平主應力σH＞最小水平主應力σh＞上覆巖體的自重應力σV，σV=γH;其中γ為上覆巖體的重度，H為上覆巖體的埋深;②水平主應力隨深度的增加呈線性增大;③最大主應力方向為NE15°～18°，表明測孔附近地應力以NE向擠壓為主;④各測孔側(cè)壓系數(shù)均＞1，表明工程區(qū)地應力以構(gòu)造應力為主。

2.2 初始地應力場的三維有限元反演分析

根據(jù)大坪山區(qū)域地形地貌特點及巖層力學特性，選取長×寬=8 375 m×3 000 m，底部高程為 －1 000 m，以隧址區(qū)為主的計算區(qū)域，主要巖體的物理力學參數(shù)見表1。采用線彈性本構(gòu)模型，對四個子應力場分別施加位移邊界條件進行求解［12］。

表1 巖體力學參數(shù)

將鉆孔實測應力和回歸計算值進行對比，如圖2所示?？芍?在測試深度范圍內(nèi)，兩者大部分較接近;鉆孔附近地應力場以構(gòu)造應力場為主。表明所采用的計算模型是可行的，適合隧址區(qū)的整體地應力場反演。

圖2 鉆孔ZK11地應力實測值與計算值對比

2.3 大坪山隧址區(qū)的地應力場特征

圖3給出了大坪山隧道軸線縱剖面的最大水平主應力和最小水平主應力等值線分布圖?？芍?σH，σh隨深度呈現(xiàn)近似線性增長的趨勢;應力的分布與地形地貌關系密切，在淺部，受地形地貌影響較大，地形變化較大時，其應力等值線相對較密;隨著埋深增大，其影響逐漸變?nèi)酢＿@與實際地應力場規(guī)律一致，說明大坪山隧道縱剖面處計算的地應力場是合理的。此外，應力等值線受地質(zhì)構(gòu)造的影響也較大［13］，在斷層及其影響帶中，應力等值線出現(xiàn)突變現(xiàn)象。

圖3 大坪山隧道軸線縱剖面應力分布(單位:MPa)

根據(jù)分析結(jié)果，隧道軸線縱剖面處屬于中等應力水平，應注意有發(fā)生中等巖爆的可能。

3 巖爆判據(jù)及預測

對于大坪山特長深埋隧道的巖爆預測研究，選擇了4種國內(nèi)外代表性的判據(jù)，即 Russenes判據(jù)、Turchaninov判據(jù)、Rb/σ1判據(jù)及陶振宇判據(jù)，如表2所示(選取原則:最大限度利用巖石的物理力學參數(shù)及有限元分析結(jié)果)。分別對隧道各個洞段依次分析，在多種判據(jù)判斷的基礎上，綜合評價了各個洞段發(fā)生巖爆的可能性和強度［14］。隧洞截面應力采用式(1)～式(3)計算。

表2 巖爆判據(jù)

式中，σmax為橫截面最大初始應力，α為最大水平主應力與隧洞走向的夾角。

預測時各洞段的最大主應力σ1和拱壁切向應力σθ采用前面有限元分析的數(shù)據(jù)，對四種判據(jù)預測結(jié)果綜合評價得出沿線巖爆發(fā)生的可能性及強度，并將已開挖段巖爆發(fā)生情況與預測結(jié)果進行對比，見表3～表7。

表3 Russenes判據(jù)巖爆預測結(jié)果

表4 Turchaninov判據(jù)巖爆預測結(jié)果

表5 Rb/σ1判據(jù)巖爆預測結(jié)果

表6 陶振宇判據(jù)巖爆預測結(jié)果

表7 大坪山隧道各斷面巖爆綜合預測結(jié)果及與實際情況對比

根據(jù)室內(nèi)巖石物理力學試驗結(jié)果，取灰?guī)r和白云巖天然巖樣的單軸抗壓強度Rc=100.5 MPa，飽和單軸抗壓強度Rb=63.0 MPa(僅考慮硬巖段的灰?guī)r和白云巖)。

由巖爆綜合預測結(jié)果可知，大坪山隧道的硬巖段巖爆發(fā)生情況:除了3 301 m長度范圍內(nèi)不發(fā)生巖爆外，弱巖爆區(qū)長2 900 m，中等巖爆區(qū)長600 m，強烈?guī)r爆區(qū)長3 00 m(埋深900 m處完整堅硬白云巖段)，隧道開挖過程中應特別注意中等—強烈?guī)r爆發(fā)生區(qū)段，及時做好支護。

大坪山隧道出口段左右洞尚處于頁巖段，所以不考慮發(fā)生巖爆。進口段施工進入巖爆區(qū)段后未發(fā)生巖爆，這與沿線出現(xiàn)多條斷層破碎帶、隧道埋深不大、地下水相對較發(fā)育等有關。但隨著掘進的深入，隧道出口段將逐漸進入硬巖段，進口段埋深也隨之增大，洞內(nèi)出現(xiàn)巖爆的可能性增大。因此，在隧道開挖施工過程中應加強監(jiān)測、開展超前預報工作，并結(jié)合工程建設采取巖爆防治措施，避免災害發(fā)生。

4 結(jié)語

1)通過三維有限元回歸分析，對大坪山深埋特長隧道兩個鉆孔的實測應力進行分析，計算值與實測值對比結(jié)果表明所選模型是可行的，可用來進行隧址區(qū)的整體地應力場的反演。

2)采用4種國內(nèi)外常用的巖爆判據(jù)對隧道各斷面進行巖爆預測，綜合評價各斷面發(fā)生巖爆的可能性及強度，預測有600 m長為中等巖爆區(qū)，300 m長為強烈?guī)r爆區(qū)。但以上預測是根據(jù)設計資料進行的，還需要結(jié)合施工階段巖爆的發(fā)生情況進行現(xiàn)場檢驗。

本文巖爆預測僅采用了理論法中的4種判據(jù)，未涉及能量法及現(xiàn)場實測法，且對巖爆發(fā)生的影響因素未做分析，還需結(jié)合現(xiàn)場監(jiān)測和測試，分析影響巖爆的因素，以更準確判斷巖爆發(fā)生的可能性。

［1］張永雙，熊探宇，杜宇本，等.高黎貢山深埋隧道地應力特征及巖爆模擬試驗［J］.巖石力學與工程學報，2009，28(11):2286-2294.

［2］康勇，李曉紅，王青海，等.隧道地應力測試及巖爆預測研究［J］.巖土力學，2005，26(6):959-963.

［3］張鏡劍，傅冰駿.巖爆及其判據(jù)和防治［J］.巖石力學與工程學報，2008，27(10):2034-2042.

［4］肖清華，秦陽，代承勇，等.高地應力下大斷面隧洞施工巖爆工程對策［J］.鐵道建筑，2012(10):76-78.

［5］吳德興，楊健.蒼嶺特長公路隧道巖爆預測和工程對策［J］.巖石力學與工程學報，2005，24(21):3965-3971.

［6］陶振宇，潘別桐.巖石力學原理與方法［M］.武漢:中國地質(zhì)大學出版社，1991.

［7］張倬元，王士天，王蘭生.工程地質(zhì)分析原理［M］.北京:地質(zhì)出版社，1994.

［8］朱宏銳.錦屏引水隧洞巖爆特征及其影響因素分析［J］.鐵道建筑，2009(5):71-73.

［9］徐士良，崔振東.秦嶺公路隧道2號豎井地應力與巖爆分析［J］.工程地質(zhì)學報，2010，18(3):407-412.

［10］許博，謝和平，涂揚舉.瀑布溝水電站地下廠房開挖過程中巖爆應力狀態(tài)的數(shù)值模擬［J］.巖石力學與工程學報，2007，26(增1):2894-2900.

［11］中國科學院武漢巖土力學研究所.湖北省谷竹高速公路大坪山隧道ZK10(補)、ZK11(補)孔地應力測試報告［R］.武漢:中國科學院武漢巖土力學研究所，2010.

［12］張奇華，鐘作武，龔壁新.施加邊界位移產(chǎn)生純剪應力及反分析應用［J］.長江科學院院報，2000，17(2):34-36.

［13］KONIETZKY H，KAMP L，HAMMER H，et al.Numerical modelling of in situ stress conditions as an aid in route selection for rail tunnels in complex geological formations in South Germany［J］.Computers and Geotechnics，2001，28(6):495-516.

［14］呂慶，孫紅月，尚岳全，等.深埋特長公路隧道巖爆預測綜合研究［J］.巖石力學與工程學報，2005，24(16):2982-2988.