存在局部破碎帶條件下海底隧道滲透穩(wěn)定性分析

王建新，鄭穎人，包林海

(1．中國地震局地殼應(yīng)力研究所，北京 100085;2．后勤工程學院，重慶 400041)

0 引言

目前，世界上有多條海底隧道正在籌劃中，如日韓海底隧道、直布羅陀海峽海底隧道、意大利墨西拿海峽海底隧道等，已建成的如日本青函隧道、英法海峽隧道、丹麥大海峽隧道、日本東京灣渡海公路隧道等。隧道技術(shù)在我國己經(jīng)得到了突飛猛進的發(fā)展，應(yīng)用前景廣闊。近些年來，水下隧道以其獨特的優(yōu)勢得到了專家的青睞，目前跨江河的水下交通隧道己有很多成功的工程范例，許多穿江(河)隧道己形成多重環(huán)線交通。山東青島的青黃海底公路隧道以及渤海海峽、長江口、杭州灣、珠江口、伶仃洋及瓊州海峽等地區(qū)建造海隧道的規(guī)劃也在醞釀之中，在此基礎(chǔ)上以上海為中心的連接周邊港口、城市、島嶼的高速公路網(wǎng)和環(huán)太平洋海岸的南北公路干線亦可能建成。這些交通線網(wǎng)上的各種隧道工程，在海底、大江或大河下修建的交通隧道，尤其是海底隧道，不同于穿江隧道，有自身的特殊性，海底隧道的設(shè)計和施工，都面臨很多關(guān)鍵技術(shù)問題亟待解決，如安全系數(shù)與可靠性指標［1］。

繼我國第一條海底隧道——廈門海底隧道正式施工后，青島的青黃海底公路隧道也已開工建設(shè)，初步地質(zhì)勘測顯示，大部分圍巖情況良好，多為微風化花崗巖和火成巖，無大的斷層和軟弱結(jié)構(gòu)面通過，巖體完整性好，地質(zhì)條件良好。但是隧道通過線路內(nèi)局部存在斷層破碎帶，需要在設(shè)計、施工中重點考慮。本文以該隧道為研究背景，計算分析海水沿破碎帶入滲后隧道圍巖滲流場的分布，應(yīng)用滲流分析有限元及強度折減理論，分別進行堵水和不堵水2種工況條件下，隧道的滲透穩(wěn)定性流固耦合分析;并討論破碎帶條件下海水入滲在隧道設(shè)計支護中的考慮。

1 巖體滲透穩(wěn)定性計算分析理論

1．1 滲流計算數(shù)學模型

模擬巖土體內(nèi)滲流問題基本數(shù)學模型如下［2］:

式中:Kx，Ky，Kz，μ 分別為含水層的滲透系數(shù)和貯水率;t為時間;Ω為滲流問題區(qū)域;Γ1為水頭邊界;Γ2為流量邊界;n為外法線;H0，H1分別為0時刻和t時刻的水頭。上述滲流微分方程一般采用伽遼金加權(quán)余量求解法。在實際計算中，多采用8節(jié)點等參元對空間進行離散，采用隱式差分格式對時間進行離散。在此基礎(chǔ)上，編寫了飽和滲流計算程序。

1．2 考慮滲透作用下的強度折減法

考慮滲透作用的強度折減法基本原理為:采用滲流有限元程序計算出的海水入滲條件下，含破碎帶巖體在堵水和不堵水條件下圍巖的滲流場分布，并求解其梯度分布。在已知梯度條件下，根據(jù)有限元理論可方便得到由單元滲透體積力等效的節(jié)點荷載，具體計算公式如下:

式中:{FS}為滲透體積力引起的等效節(jié)點荷載;γw為水的容重;［N］為形函數(shù);H為滲流場水頭函數(shù)。

強度折減法在隧道中的應(yīng)用目前已被巖土界逐步認可，其基本思想就是在理想彈塑性有限元計算中將隧道巖體抗剪切強度參數(shù)(黏聚力和內(nèi)摩擦角)逐漸降低直到其達到極限破壞狀態(tài)為止，此時程序可以自動根據(jù)其彈塑性有限元計算結(jié)果得到隧道的破裂面，同時得到隧道的強度儲備的安全系數(shù)ω。隧道工程中的強度折減形式與邊坡穩(wěn)定分析的傳統(tǒng)極限平衡條分法安全系數(shù)定義形式是一致的，基本公式如下:

式中:ω 為安全系數(shù);c，φ，c'，φ'分別為巖體初始和折減后的黏聚力和摩擦角。極限狀態(tài)下的ω為隧道穩(wěn)定系數(shù)［3-5］。

基于上述基本理論，本文采用ANSYS軟件進行強度折減后的隧道穩(wěn)定性計算。此外，應(yīng)用Fortran語言對ANSYS軟件進行了二次開發(fā)，編制了滲透穩(wěn)定性計算部分的程序。本文即采用該程序來計算考慮海水滲透作用下的滲流應(yīng)力耦合穩(wěn)定性計算。

2 隧道計算模型參數(shù)及工況

青黃海底隧道進口(團島端)在團島路小瞿塘峽路交匯口附近出洞，后沿四川路向北延伸;其另一端在薛家島，在島上北莊村和后岔灣村之間出洞［6］。隧道全長6．17 km，埋深線路平直，灣口隧址大部分圍巖情況良好，多為微風化花崗巖和火成巖，無大的斷層和軟弱結(jié)構(gòu)面通過，巖體完整性好，地質(zhì)條件良好，不良地質(zhì)地段總長小于100 m，海水最大深度40 m。隧道數(shù)值模型如圖1所示，模型計算范圍兩側(cè)取4倍隧道寬度，上下取2倍隧道高度。依據(jù)地質(zhì)勘測報告，破碎帶近似取40°傾角，堵水加固厚度為3 m。

圖1 海底隧道計算模型Fig．1 Numerical model of subsea tunnel

為了比較破碎帶對隧道圍巖滲流場的影響，文中首先應(yīng)用滲流有限元計算方法計算了完整圍巖條件下的滲流場分布，在此基礎(chǔ)上，對海水沿破碎帶入滲條件下隧道圍巖的滲流場進行了模擬，隨后計算堵水條件下隧道圍巖滲流場分布特征。滲流有限元模擬中，取海底平面為第一邊界條件，即海水深度的水頭邊界;模型兩側(cè)為不透水邊界;底部為下滲流出梯度邊界。考慮到隧道圍巖滲流特征的不均勻性，取隧道開挖面節(jié)點為溢出點，通過程序內(nèi)溢出邊界迭代搜索技術(shù)，能更好地模擬隧道圍巖的滲流狀態(tài)特征。

經(jīng)過數(shù)值計算研究得出，隧道圍巖處于極限狀態(tài)時圍巖發(fā)生塑性應(yīng)變突變時的情況就是圍巖發(fā)生流動破壞的情況，此時，恰好計算不收斂，因而可依據(jù)塑性應(yīng)變和位移突變來確定潛在破壞面及對應(yīng)的安全系數(shù)。利用強度折減有限元法求得安全系數(shù)與潛在滑動面，不僅可以評價隧道的穩(wěn)定性和設(shè)計的合理性，還可以對支護參數(shù)和施工工藝提出改進建議。在考慮滲透荷載的數(shù)值模擬計算中，為更好地比較計算結(jié)果，應(yīng)用有限元強度折減法對隧道完整圍巖無滲透作用條件下的穩(wěn)定性進行計算分析。在此基礎(chǔ)上，本文對海水入滲荷載對圍巖穩(wěn)定的影響進行了模擬，并計算了堵水和不堵水2種工況。具體工況見表1。

表1 計算工況列表Table 1 List of calculation cases

在滲透穩(wěn)定性計算中，上邊界按照海水壓力面荷載施加，模型兩側(cè)為水平位移約束，模型底部固定約束。巖體物理力學參數(shù)見表2。根據(jù)勘測報告提供的巖塊強度，參考有關(guān)資料，計算時實際巖體強度大約為巖塊強度的1/6。

表2 巖體物理力學參數(shù)Table 2 Physical and mechanical parameters of rock mass

3 計算結(jié)果分析

3．1 隧道滲流場計算結(jié)果分析

各工況滲流場水頭和水力梯度分布情況見圖2—4。由圖2可以看出:無破碎帶存在，低滲透性條件下，隧道圍巖滲流場分布均勻，且水頭不高。海水滲透梯度分布規(guī)律，靠近隧道開挖面隧道梯度分布密集，但梯度最大值為4．3，小于隧道圍巖的容許最大水力梯度40。隧道圍巖入滲流量較小，為0．32 m3。

由圖3可以看到:滲透性較大的破碎帶的存在，使海水沿著破碎帶入滲到隧道圍巖周圍。破碎帶起到了泄水通道的作用，隧道圍巖的水頭降低，尤其在破碎帶周圍，水頭降低非常明顯;梯度值也較小，隧道圍巖周圍最大梯度值為2．6，小于巖體的容許最大水力梯度40。但是，破碎帶的存在，使隧道圍巖滲流量明顯加大，達到3．6 m3/d。隧道圍巖破碎帶的存在造成滲流的加大對隧道的施工安全穩(wěn)定性影響很大，這不僅會加大施工的難度，而且在采取措施不及時的情況下，極有可能造成隧道突水的災(zāi)害性后果。在海底隧道工程中，隧道突水是一個必須嚴格防范和注意的問題。

圖2 完整圍巖海底隧道滲流場水頭與梯度分布Fig．2 Distribution of water head and hydraulic gradient in seepage field of subsea tunnel in intact rock mass

圖3 存在破碎帶海底隧道滲流場水頭與梯度分布Fig．3 Distribution of water head and hydraulic gradient in seepage field of subsea tunnel in rock mass with fracture zones

按照設(shè)計，破碎帶堵水圈厚度為3m，由圖4知:堵水后隧道內(nèi)圍巖水頭明顯上升，特別在堵水部位，水頭值較高，水力梯度也較大，達到15，但依然小于堵水圈巖體和隧道圍巖的容許水力梯度40。堵水處理后隧道內(nèi)滲流量明顯減小，為0．36 m3/d，幾乎和完整圍巖滲流量大小相似。低滲透性的堵水圈巖體使破碎帶下滲的海水聚集在堵水部位，造成了高水頭和較高的梯度，可見需要提高堵水部位巖體的加固強度來預(yù)防由此可能造成的突水災(zāi)害。

圖4 堵水加固后海底隧道圍巖滲流場水頭與梯度分布Fig．4 Distribution of water head and hydraulic gradient in seepage field of subsea tunnel after water-sealing consolidation

上述計算結(jié)果基本正確反應(yīng)了海水沿破碎帶入滲后隧道圍巖范圍滲流場的分布特征。根據(jù)本文的滲流有限元計算方法，可以較為精確地描述隧道滲流場的分布和演化，并通過各工況對比分析，來反應(yīng)不同地質(zhì)條件以及封堵治理措施下，隧道圍巖的滲流場情況，為隧道的設(shè)計和災(zāi)害治理提供重要的參考。

3．2 考慮滲透荷載的穩(wěn)定性分析

應(yīng)用Fortran語言對ANSYS軟件進行了二次開發(fā)，編制了滲透穩(wěn)定性計算部分的程序。本文即采用該程序來考慮海水滲透作用下的滲流應(yīng)力耦合穩(wěn)定性計算［7-9］?？紤]滲透力作用下，完整圍巖隧道的安全系數(shù)為7．21;破碎帶的存在使海底隧道整體安全性大大降低，安全系數(shù)降為2．20;堵水加固后隧道整體安全系數(shù)為5．6。各個工況下隧道圍巖在極限狀態(tài)下的塑性區(qū)和潛在的破裂面如圖5—7所示。

由圖5—7可以看出:完整圍巖條件下，隧道圍巖潛在破裂面發(fā)生在兩側(cè)，呈對稱的圓弧狀，極限狀態(tài)下的最大的塑性應(yīng)變值為0．003;由于局部破碎帶的存在，海底隧道圍巖塑性區(qū)范圍明顯擴大，主要集中在隧道兩側(cè)以及破碎帶內(nèi)，此條件下隧道圍巖內(nèi)最大的塑性應(yīng)變位于破碎帶內(nèi)，最大值為0．017(見圖6);對破碎帶進行堵水加固后，隧道圍巖周圍塑性區(qū)范圍減小，最大塑性應(yīng)變位于堵水部位與破碎帶結(jié)合的部位，這主要是由于堵水后在堵水圈外的破碎帶內(nèi)產(chǎn)生了較大水力梯度的緣故。由圖7可以看出:加固堵水后，極限狀態(tài)下隧道的潛在破裂面位于堵水部位，所以，在堵水加固的同時，要注意提高堵水圈的支護強度，以免造成隧道突水災(zāi)害。

圖5 完整圍巖海底隧道塑性區(qū)和潛在破裂面分布圖Fig．5 Distribution of plasticized zone and potential fracture plane of subsea tunnel in intact rock mass

4 滲透壓力在隧道支護設(shè)計中的考慮

對海底隧道而言，圍巖滲水對海底隧道的最大危害就是造成突發(fā)的頂板透水以及隨之而來的頂板塌落事故，如挪威Oslofjord海底隧道。所以在對海底隧道通過線路地質(zhì)條件勘測準確以后，如何確定作用在隧道圍巖上的海水壓力，是海底隧道設(shè)計、施工中需要重點考慮的問題。海底隧道設(shè)計中，作用在隧道襯砌上的水壓力究竟該取多大，工程界的意見尚不統(tǒng)一。隧道襯砌設(shè)計中考慮全水頭水壓力的情況基本上有2種:1)土質(zhì)隧道以及上覆巖層較為破碎且滲透性較強的隧道，隧道襯砌周邊很快達到全水頭;2)采用盾構(gòu)法施工的隧道，施工工藝要求隧道設(shè)計為不排水。第2種情況，即使隧道圍巖完整性較好，在隧道修建初期，雖然地下水沿裂隙滲入水量很少，但經(jīng)過若干年以后，作用在隧道襯砌的外滲透壓力也會達到全水頭壓力，這種狀況已被重慶市某長江水下隧道的試驗所證實，假如隧道設(shè)計為不排水，襯砌上必須考慮全水頭壓力。

圖6 存在破碎帶海底隧道塑性區(qū)和潛在破裂面分布圖Fig．6 Distribution of plasticized zone and potential fracture plane of subsea tunnel in rock mass with fracture zones

圖7 堵水加固后海底隧道塑性區(qū)和潛在破裂面分布圖Fig．7 Distribution of plasticized zone and potential fracture plane of subsea tunnel after water-sealing consolidation

顯然，如果隧道上覆巖層完整性非常好，滲透性不大，圍巖內(nèi)裂隙分布范圍很小，海水只是沿圍巖裂隙少量流出，按挪威海底隧道規(guī)范，允許滲流量為300 L/(km·min)。這種情況下，沒有必要進行堵水，隧道的整體穩(wěn)定性問題不大，計算分析可以忽略滲透壓力的作用，在此地段圍巖構(gòu)筑薄型襯砌就可滿足隧道穩(wěn)定性要求，甚至只進行混凝土噴面就可以保持隧道穩(wěn)定，在這方面挪威的經(jīng)驗值得借鑒［10-11］。

上述計算結(jié)果顯示，海底隧道的設(shè)計與開挖施工可采用“排堵結(jié)合”的措施來降低水的滲流。尤其在不良地段，海水滲流主要通過局部破碎帶、斷層及部分裂隙侵入隧道，因而必須對破碎帶進行注漿堵水，降低海水的入滲。漏水是必然的、正常的，沒有必要完全堵水，完全堵水會造成堵水部位水力梯度的急劇增大，對堵水部位非常不利;而且適量的排水對隧道襯砌有卸壓的作用，有利于隧道穩(wěn)定。通過對封堵加固后隧道的滲透穩(wěn)定分析，加固堵水后，極限狀態(tài)下隧道的潛在破裂面位于堵水部位，所以，在堵水加固的同時，要注意提高堵水圈的支護強度，在不良地段內(nèi)適當加厚襯砌，并增設(shè)鋼拱架或長錨索等設(shè)施，以免造成隧道突水災(zāi)害。

5 結(jié)論

基于滲流有限元計算方法，模擬了海水沿破碎帶入滲條件下，海底隧道圍巖范圍的滲流場分布特征;應(yīng)用Fortran語言對ANSYS軟件進行了二次開發(fā)，編制了滲透穩(wěn)定性計算部分的程序。采用該程序來考慮海水滲透作用下的滲流應(yīng)力耦合穩(wěn)定性計算。

1)不堵水情況下，海水沿著破碎帶入滲，使隧道圍巖內(nèi)水頭降低，水力梯度變小，但是隧道內(nèi)會產(chǎn)生較大的滲流量;堵水加固后，隧道內(nèi)滲流量明顯減小，但是在破碎帶封堵部位，產(chǎn)生較高的水力梯度，對隧道的穩(wěn)定性不利，尤其要提高封堵部位的支護強度。

2)考慮滲透壓力影響條件下，破碎帶的存在使隧道的安全穩(wěn)定性系數(shù)大大降低，塑性區(qū)范圍擴大，塑性應(yīng)變值增大;堵水加固后，隧道的整體穩(wěn)定性提高，但是由于堵水部位較大的水力梯度造成較大的滲透壓力，使該部位成為潛在的破壞點，應(yīng)該在堵水加固的同時提高支護強度。

3)對海底隧道而言，如果破碎帶范圍較小，可采用“排堵結(jié)合”的措施，沒有必要對破碎帶全封閉，這樣容易造成破碎帶較大的水力梯度;在堵水加固的同時，要注意提高堵水圈的支護強度，在不良地段內(nèi)適當加厚襯砌，并增設(shè)鋼拱架或長錨索等設(shè)施，以免造成隧道突水災(zāi)害。

［1］ DAHLO T S，NILSEN B．Stability and rock cover of hard rock subsea tunnels［J］．Tunneling and Underground Space Technology，1994，9(2):151-158．

［2］ 榮冠，張偉，周創(chuàng)兵．降雨入滲條件下邊坡巖體飽和非飽和滲流計算［J］．巖土力學，2005，26(10):1545-1550．(RONG Guan，ZHANG Wei，ZHOU Chuangbing．Numerical analysis of saturated unsaturated seepage problem of rock slope under rainfall infiltration［J］．Rock and Soil Mechanics，2005，26(10):1545-1550．(in Chinese))

［3］ 李廷春，李術(shù)才，陳衛(wèi)忠，等．廈門海底隧道的流固耦合分析［J］．巖土工程學報，2004，26(3):397-401．(LI Tingchun，LI Shucai，CHEN Weizhong，et al．Coupled fluidmechanical analysis of Xiamen subsea tunnel［J］．Chinese Journal of Geotechnical Engineering，2004，26(3):397-401．(in Chinese))

［4］ 張勝軍，孫云志，蔣小娟，等．廈門海底隧道工程巖體滲透穩(wěn)定性試驗研究［J］．人民長江，2005，36(3):22-23．(ZHANG Shengjun，SUN Yunzhi，JIANG Xiaojuan，et al．Test and study on the rock mass seepage stability of eastroute sea bottom tunnel works in Xiamen city［J］．Yangtze River，2005，36(3):22-23．(in Chinese))

［5］ 鄭穎人，趙尚毅．巖土工程極限分析有限元法及其應(yīng)用［J］．土木工程學報，2005，38(1):91-98，104．(ZHENG Yingren，ZHAO Shangyi．Limit state finite element method for geotechnical engineering analysis and its applications［J］．China Civil Engineering Journal，2005，38(1):91-98，104．(in Chinese))

［6］ 張明聚，郜新軍，郭衍敬．海底隧道突水分析及其在翔安隧道中的應(yīng)用［J］．北京工業(yè)大學學報，2007，33(3):273-277．(ZHANG Mingju，GAO Xinjun，GUO Yanjing．Analysis of water inrush in undersea tunnel and its application in Xiang'an tunnel［J］．Journal of Beijing University of Technology，2007，33(3):273-277．(in Chinese))

［7］ 劉才華，陳從新，付少蘭．剪應(yīng)力作用下巖體裂隙滲流特性研究［J］．巖石力學與工程學報，2003，22(10):1651-1655．(LIU Caihua，CHEN Congxin，F(xiàn)U Shaolan．Study on seepage characteristics of a single rock fracture under shear stresses［J］．Chinese Journal of Rock Mechanics and Engineering，2003，22(10):1651-1655．(in Chinese))

［8］ 張黎明，鄭穎人，王在泉，等．有限元強度折減法在公路隧道中的應(yīng)用探討［J］．巖土力學，2007，28(1):97-101，106．(ZHANG Liming，ZHENG Yingren，WANG Zaiquan，et al．Application of strength reduction finite element method to road tunnels［J］．Rock and Soil Mechanics，2007，28(1):97-101，106．(in Chinese))

［9］ 陳衛(wèi)忠，楊建平，楊家?guī)X，等．裂隙巖體應(yīng)力滲流耦合模型在壓力隧洞工程中的應(yīng)用［J］．巖石力學與工程學報，2006，25(12):2384-2391．(CHEN Weizhong，YANG Jianping，YANG Jialing，et al．Hydro-mechanical coupled model of jointed rock mass and its application to pressure tunnels［J］．Chinese Journal of Rock Mechanics and Engineering，2006，25(12):2384-2391．(in Chinese))

［10］ 孫鈞．海底隧道工程設(shè)計施工若干關(guān)鍵技術(shù)的商榷［J］．巖石力學與工程學報，2006，25(8):1513-1521．(SUN Jun．Discussion on some key technical issues for design and construction of undersea tunnels［J］．Chinese Journal of Rock Mechanics and Engineering，2006，25(8):1513-1521．(in Chinese))

［11］ Eisenstein Z D．Large undersea tunnels and the progress of tunneling technology［J］．Tunneling and Underground Space Technology，1994，19(3):283-292．