楊彥亮，李錫波，毛洪錄

(1.聊城市公路工程總公司第三工程處，山東聊城 252000;2.山東省交通規(guī)劃設(shè)計院，山東濟南 250031;3.山東大學土建與水利學院，山東濟南 250061)

0 引言

隨著我國交通事業(yè)的迅猛發(fā)展和國家西部大開發(fā)戰(zhàn)略的持續(xù)推進，復雜地形、地質(zhì)條件下的深挖高填工程日趨增多，尤其是長大、深埋的隧道工程。深埋隧道中因高地應力的存在，易引起巖爆、流變、斷層軟巖擠入大變形等典型災害，為工程勘測、設(shè)計及施工帶來了嚴峻的挑戰(zhàn)［1］，因此亟需對高地應力條件下隧道的施工響應進行模擬、分析。因其顯著的物理非線性、結(jié)構(gòu)非線性和接觸邊界非線性等幾何和力學特征，傳統(tǒng)數(shù)學手段遇到非常大的困難，因此現(xiàn)階段主要采用數(shù)值分析和模型試驗等手段對其進行研究。趙旭峰等［1］使用ANSYS軟件對深埋軟巖隧洞施工力學形態(tài)和變形時空效應進行了三維非線性黏彈性數(shù)值模擬分析，并基于與現(xiàn)場實測數(shù)據(jù)的對比分析得到了一些定性的認識。張志強等［2］和陶波等［3］使用FLAC對高地應力條件下軟弱圍巖隧道的變形規(guī)律進行了數(shù)值模擬分析。朱維申等［4］對高地應力條件下地下洞室開挖過程中的力學及變形規(guī)律進行了模型試驗和數(shù)值模擬分析。劉招偉等［5］針對烏鞘嶺隧道高地應力段出現(xiàn)的初期支護大變形問題，通過模型試驗對不同支護條件下的隧道變形和支護結(jié)構(gòu)內(nèi)力進行了監(jiān)測、分析。

關(guān)于不同開挖方法對隧洞施工過程力學和變形響應的研究分析，現(xiàn)階段以數(shù)值手段為主。段慧玲等［6］對不同圍巖等級條件下大跨度隧道施工方法進行了經(jīng)驗和數(shù)值比選，并給出了相應合理的開挖方法。湯勁松等［7］結(jié)合實際工程，對破碎圍巖下大跨扁平隧道的開挖和支護方式進行了數(shù)值和理論分析。然而，現(xiàn)有成果對下導洞超前開挖法的相關(guān)研究和經(jīng)驗介紹較少，特別是對高地應力條件下隧道開挖方法如何選取、不同開挖方式的影響規(guī)律和適應條件的研究更為少見。

某隧道由于埋深較大，具有顯著的高地應力特征，原設(shè)計中深埋段采用全斷面法，但施工過程中多次出現(xiàn)圍巖剝落掉塊、支護結(jié)構(gòu)扭曲破壞等現(xiàn)象，因此擬改進開挖方式。由于隧道圍巖巖性為砂巖，不適宜采用分部過多的CD法、CRD法及預留核心土環(huán)形導坑等開挖方法，經(jīng)初步比選擬采用下導洞超前開挖法或上下臺階法，為論證擬采取方案的可行性，本文利用有限差分軟件FLAC3D針對高地應力條件，對不同開挖方式的隧道斷面位移場、塑性區(qū)及應力場的分布特征進行數(shù)值模擬和對比分析，并據(jù)此推薦了最優(yōu)施工方法。

1 工程概況

隧道采用分幅式，全長約7 500 m，隧道最大埋深約678 m，屬特長深埋隧道。隧道屬構(gòu)造剝蝕、侵蝕中低山深切溝谷地貌區(qū)，地形切割較深，溝谷斜坡地貌發(fā)育，峰頂呈渾圓狀，丘脊寬緩，山嶺穿越區(qū)地形切割較強，山間溝谷發(fā)育，溝谷峽窄。施工過程中開挖揭露圍巖發(fā)育有褶皺等小型構(gòu)造，圍巖軟硬交替，且具有高地應力特點，施工中多處出現(xiàn)偏幫、擠壓掉快，圍巖沿錨桿層狀剝落致使支護扭曲破壞等現(xiàn)象，如圖1所示。

圖1 現(xiàn)場支護結(jié)構(gòu)扭曲及擠壓掉塊Fig.1 Phenomena of support structure distortion and localized rock fall

2 隧道施工數(shù)值模擬

本文選取三維快速拉格朗日有限差分軟件FLAC3D進行數(shù)值模擬計算，其在材料的彈塑性分析、大變形分析以及模擬施工過程等領(lǐng)域有其獨到的優(yōu)點。針對所研究內(nèi)容，根據(jù)工程地質(zhì)、構(gòu)造特點，并結(jié)合類似數(shù)值計算經(jīng)驗［8－11］選取并確定了數(shù)值計算模型。

2.1 計算模型

根據(jù)圣維南原理和相關(guān)經(jīng)驗、結(jié)論可知隧道開挖影響范圍為洞徑的3～5倍，對數(shù)值模型進行了適當?shù)哪Ｐ秃喕?。模型計算范?水平方向x軸長度取120 m，豎直方向y軸取120 m，縱向z軸沿隧道軸線方向取90 m。

模型左、右、前、后和下部邊界均施加法向約束，模型上部施加邊界應力，其等效地應力由σz=γh確定(γ為上覆巖層的平均加權(quán)容重;h為上覆巖層總厚度［11］，計算取 γ =24 kN/m3，h=620 m，得 σz=14.88 MPa)，水平應力結(jié)合隧道地質(zhì)資料及相關(guān)應力數(shù)據(jù)資料按垂直地應力的1.2倍取值，計算模型如圖2所示。

圖2 隧道計算模型Fig.2 Computation model of tunnel

2.2 力學參數(shù)及邊界條件

計算中屈服準則采用Mohr-Coulomb準則，并考慮巖體的受拉屈服、彈塑性變形及大變形。圍巖材料模型采用 Mohr-Coulomb理想彈塑性模型，開挖采用FLAC3D中的Null模型，復合式支護中的初期支護采用FLAC3D中的Shell結(jié)構(gòu)單元。

初期支護采用10 cm厚的C20噴射混凝土，φ8鋼筋網(wǎng)片。對于支護中采用的錨桿(φ22砂漿錨桿，長3 m，梅花型布置，間距1.2 m)，計算過程中根據(jù)作用等效原則來考慮，即提高圍巖的黏聚力和摩擦角來替代錨桿的作用，由于摩擦角改變較小，不予考慮，錨固圍巖體的黏聚力

式中:C0為未加錨桿時圍巖的黏聚力;C為加錨桿時圍巖的黏聚力;τ為錨桿最大抗剪應力;Sm為錨桿的面積;a，b分別為錨桿的縱、橫向間距;η為經(jīng)驗系數(shù)，可取2～5。

鋼拱架的作用可采用等效的方法來考慮，即將鋼拱架彈性模量折算給噴射混凝土［13］，按式(2)進行折算。

式中:E為折算后混凝土彈性模量;E0為原混凝土彈性模量;Sg為鋼拱架截面積;Eg為鋼材彈性模量;Sc為混凝土截面積。

根據(jù)該隧道工程地質(zhì)、水文地質(zhì)及相關(guān)地質(zhì)勘測資料，圍巖物理力學參數(shù)如表1所示。

表1 圍巖物理力學參數(shù)表Table 1 Physical and mechanical parameters of surrounding rock

2.3 開挖方案

為模擬實際施工開挖順序，沿z軸方向以3 m長度為一個開挖循環(huán)進行模擬;采用如下開挖方案:分15步開挖，共開挖45 m(自Z=0開挖至Z=45)，下導洞超前開挖法先行開挖下導洞，而后開挖上弧部，臺階法先行開挖上臺階，而后開挖下臺階。以Z=18斷面的位移場、塑性區(qū)、應力場為研究對象，分析判斷圍巖的穩(wěn)定性。

3 數(shù)值模擬結(jié)果分析

3.1 位移場分析

隧道水平收斂和拱頂下沉量測是監(jiān)控量測的主要內(nèi)容，也是隧道圍巖應力狀態(tài)變化的最直觀反映［14］。表2為采用不同開挖方式時隧道周邊位移的計算值，采用下導洞超前開挖法施工時，隧道拱頂下沉、水平收斂和拱肩位移(此處拱肩位移指上部4/5隧道內(nèi)斷面高度處測點在xy平面內(nèi)的合位移)最大，采用全斷面法時，拱頂下沉最小，采用臺階法時，水平收斂和拱肩位移最小。

表2 隧道周邊位移計算值Table 2 Calculated displacement of surrounding rock mm

巖體工程開挖變形具有很強的時空效應，隧道開挖后受開挖面的約束影響，圍巖應力逐步釋放，直到開挖面空間約束效應完全消失，圍巖應力才得以全部釋放［1，15－16］，因此采用不同的施工方式，隧道圍巖應力釋放過程不同，從而使巖體卸荷速率不同?，F(xiàn)以水平收斂值為研究對象，分析開挖方式對卸荷速率的影響。水平收斂值隨隧道開挖的變化曲線如圖3所示。

圖3 水平收斂值變化曲線Fig.3 Variation curves of horizontal convergence

分析圖3可知:水平收斂值增大主要發(fā)生在開挖Z=12～27時，單步開挖使水平收斂值產(chǎn)生最大增幅的是第7開挖步，即自Z=18～21。由全斷面開挖法、下導洞超前開挖法和上下臺階法計算得:水平收斂變形速率最大值依次為:5.61，4.06，4.07 mm/d，但下導洞超前開挖法和上下臺階法因斷面錯次開挖，其同一斷面的水平收斂速率可顯著減小。由此可見，采用下導洞超前開挖法和臺階法施工時，可充分利用開挖面的約束作用，降低隧道周邊位移變形速率，保證隧道施工安全。

3.2 塑性區(qū)分析

分別采用全斷面法、下導洞超前開挖法、臺階法開挖時，塑性區(qū)分布如圖4(a)、(b)、(c)所示。分析圖4可知:采用下導洞超前開挖時，隧道施工造成圍巖的多次應力擾動，產(chǎn)生的塑性區(qū)范圍最大，拱腰及其上部塑性區(qū)范圍約2.7 m，拱腳處塑性區(qū)范圍約1.8 m，拱底塑性區(qū)范圍最大，約3.6 m，且拱底0.9 m范圍內(nèi)曾發(fā)生拉伸破壞(下導洞開挖時發(fā)生)，不利于隧道圍巖的穩(wěn)定，施工中應盡量避免;采用全斷面法和臺階法時產(chǎn)生的塑性區(qū)大小基本一致，但塑性區(qū)集中區(qū)域不同，采用全斷面法開挖時，塑性區(qū)主要集中在拱頂和拱底，塑性區(qū)范圍約2.7 m，采用臺階法開挖時，塑性區(qū)主要集中在拱底，塑性區(qū)范圍約3.6 m。因此，從塑性區(qū)范圍及破壞機制2方面考慮，隧道施工時應優(yōu)先采用臺階法。

圖4 圍巖塑性區(qū)分布圖Fig.4 Plasticized zone of surrounding rock

3.3 應力場分析

分別采用全斷面開挖法、下導洞超前開挖法和臺階法開挖時，第一主應力分布如圖4(a)、(b)、(c)所示。分析圖4可知:采用全斷面法開挖時，應力集中主要發(fā)生在拱頂和拱底，而采用下導洞超前開挖法和臺階法開挖時，應力集中主要發(fā)生在隧道左右兩側(cè)拱腰處;分別采用上述3種方法開挖時，隧道周邊圍巖第一主應力的最大值依次為 30.007，26.805，27.128 MPa，且高應力分布范圍依次為全斷面法＞下導洞超前開挖法＞臺階法，如表3所示。

表3 第一主應力(SIG1)高應力區(qū)分布范圍Table 3 High stress area of first principal stress m

4 工程應用

根據(jù)計算分析結(jié)果，結(jié)合工程特點，最終選用上下臺階法進行隧道開挖，并現(xiàn)場選取YK112+750斷面對隧道開挖斷面的圍巖壓力、噴射混凝土的應力、支架的應力、位移以及水平收斂值和拱頂沉降進行了監(jiān)測，測點布設(shè)如圖5所示，初期支護監(jiān)測結(jié)果如圖6所示。

圖5 測點布設(shè)圖(單位:m)Fig.5 Layout of monitoring points(m)

圖6 初始支護壓力、應力分布圖Fig.6 Distribution of primary support pressure and stress

監(jiān)測周期為60 d，在監(jiān)測期間內(nèi)隧道最終水平收斂值為25.02 mm，最大收斂速率為4.26 mm/d;拱頂最終下沉量為15.68 mm，最大下沉速率為5.20 mm/d;監(jiān)測后期收斂速率小于0.2 mm/d，下沉速率小于0.1 mm/d，收斂時間曲線和拱頂下沉時間曲線均有明顯收斂趨勢。實測拱頂下沉量達計算值的93.6%，實測水平收斂值達計算值的90.0%，滿足規(guī)范穩(wěn)定要求。

5 結(jié)論與建議

針對工程實際特點，應用有限差分法對高地應力條件下深埋隧道開挖時的力學和變形響應特性進行了模擬分析，模擬結(jié)果與工程監(jiān)測數(shù)據(jù)吻合度較高。根據(jù)所建立數(shù)值模型，對不同開挖方式進行對比分析可知:綜合考慮隧道周邊位移值、水平收斂變化曲線、塑性區(qū)分布及高應力區(qū)分布等關(guān)鍵控制因素，上下臺階法較全斷面法和下導洞超前開挖法具有明顯的優(yōu)勢，更適合于該類工程，且工程應用效果較好;同時，在處理好下導洞施工時拱底小范圍內(nèi)出現(xiàn)的拉伸破壞的前提下，下導洞超前開挖法具有一定的可行性，但仍需進一步驗證。

［1］ 趙旭峰，王春苗，孔祥利.深部軟巖隧道施工性態(tài)時空效應分析［J］.巖石力學與工程學報，2007，26(2):404－409.(ZHAO Xufeng，WANG Chunmiao，KONG Xiangli.Analysis of time-space effects of construction behavior of deep soft rock tunnel［J］.Chinese Journal of Rock Mechanics and Engineering，2007，26(2):404－409.(in Chinese))

［2］ 張志強，關(guān)寶樹.軟弱圍巖隧道在高地應力條件下的變形規(guī)律研究［J］.巖土工程學報，2000，22(6):696－700.(ZHANG Zhiqiang，GUAN Baoshu.Research on the deformation rule of weak rock mass tunnel under high initial geostress［J］.Chinese Journal of Geotechnical Engineering，2000，22(6):696－700.(in Chinese))

［3］ 陶波，伍法權(quán)，郭啟良，等.高地應力環(huán)境下烏鞘嶺深埋長隧道軟弱圍巖流變規(guī)律實測與數(shù)值分析研究［J］.巖石力學與工程學報，2006，25(9):1828－1834.(TAO Bo，WU Faquan，GUO Qiliang，et al.Research on rheology rule of deep-buried long Wuqiaoling tunnel under high crustal stress by monitoring and numerical analysis［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(9):1828－1834.(in Chinese))

［4］ 朱維申，李勇，張磊，等.高地應力條件下洞群穩(wěn)定性的地質(zhì)力學模型試驗研究［J］.巖石力學與工程學報，2008，27(7):1308 － 1314.(ZHU Weishen，LI Yong，ZHANG Lei，et al.Geomechanical model test on stability of cavern group under high seostress［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27(7):1308 －1314.(in Chinese))

［5］ 劉招偉，王明勝，方俊波.高地應力大變形隧道支護系統(tǒng)的試驗研究［J］.土木工程學報，2010，43(5):111－116.(LIU Zhaowei，WANG Mingsheng，F(xiàn)ANG Junbo.Insitu study of reinforcement systems for tunnels under high geo-stress and large deformation［J］.China Civil Engineering Journal，2010，43(5):111 －116.(in Chinese))

［6］ 段慧玲，張林.大跨度公路隧道合理開挖方法對比研究［J］.土木工程學報，2009，42(9):114－119.(DUAN Huiling， ZHANG Lin. Comparative study ofrational excavation methods for large-span highway tunnels［J］.China Civil Engineering Journal，2009，42(9):114 － 119.(in Chinese))

［7］ 湯勁松，劉松玉，童立元，等.破碎巖體淺埋大跨公路隧道開挖方案對比研究［J］.巖土力學，2007，28(S):469－ 473.(TANG Jinsong，LIU Songyu，TONG Liyuan，et al.Comparative study on tunneling scheme of shallow overburden and large-span highway tunnel in cracked rock mass［J］.Rock and Soil Mechanics，2007，28(S):469 －473.(in Chinese))

［8］ 佘鍵，何川.軟弱圍巖段隧道施工過程中圍巖位移的三維彈塑性數(shù)值模擬［J］.巖石力學與工程學報，2006，25(3):623 － 629.(SHE Jian，HE Chuan.3D elastoplastic numerical simulation of surrounding rock displacement in soft surrounding rock section during construction process［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25(3):623－629.(in Chinese))

［9］ 邵祥澤，潘志存，張培森.高地應力巷道圍巖的蠕變數(shù)值模擬［J］.采礦與安全工程學報，2006，23(2):245－248.(SHAO Xiangze，PAN Zhicun，ZHANG Peisen.Numerical simulation of creep for surrounding rock under high stress［J］.Journal of Mining＆ Safety Engineering，2006，23(2):245－248.(in Chinese))

［10］ 譚代明，漆泰岳，莫陽春.側(cè)部巖溶隧道圍巖穩(wěn)定性數(shù)值分析與研究［J］.巖石力學與工程學報，2009，28(S2):3497 － 3503.(TAN Daiming，QI Taiyue，MO Yangchun.Numerical analysis and research on surrounding rock stability of lateral karst cave tunnel［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28(S2):3497－3503.(in Chinese))

［11］ 吳波，劉維寧，高波，等.地鐵分岔隧道施工性態(tài)的三維數(shù)值模擬與分析［J］.巖石力學與工程學報，2004，23(18):3081 －3086.(WU Bo，LIU Weining，GAO Bo，et al.3D simulation and analysis on construction behavior of forked metro tunnel［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23(18):3081－3086.(in Chinese))

［12］ 朱維申，李術(shù)才，陳衛(wèi)忠.節(jié)理巖體破壞機理和錨固效應及工程應用［M］.北京:科學出版社，2002.

［13］ 蔣樹屏，李建軍.公路隧道前置式洞口工法的三維數(shù)值分析［J］.巖土工程學報，2007，29(4):484－489.(JIANG Shuping，LI Jianjun.3D simulation and analysis of highway tunnels by use of advanced entrance construction method［J］.Chinese Journal of Geotechnical Engineering，2007，29(4):484－489.(in Chinese))

［14］ 金豐年，王波，蔣美榮.雙連拱隧道穿越破碎山體圍巖穩(wěn)定性數(shù)值分析［J］.巖土力學，2008，28(S):227－231.(JIN Fengnian，WANG Bo，JIANG Meirong.Numerical analysis of stability of wall-rock of double-arch tunnel through fractured mountain body［J］.Rock and Soil Mechanics，2008，28(S):227－231.(in Chinese))

［15］ 黃潤秋，黃達.高地應力條件下卸荷速率對錦屏大理巖力學特性影響規(guī)律試驗研究［J］.巖石力學與工程學報，2010，29(1):21－33.(HUANG Runqiu，HUANG Da.Experimental research on affection laws of unloading rates on mechanical properties of Jinping marble under high geostress［J］.Chinese Journal of Rock Mechanics and Engineering，2010，29(1):21－33.(in Chinese))

［16］ 黃潤秋.巖石高邊坡發(fā)育的動力過程及其穩(wěn)定性控制［J］.巖石力學與工程學報，2008，27(8):1525－1544.(HUANG Runqiu.Geodynamical process and stability control of high rock slope development［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27(8):1525－1544.(in Chinese))