朵生君，李國(guó)良，張 景

(1. 陜西省鐵道及地下交通工程重點(diǎn)實(shí)驗(yàn)室，西安 710043；2. 中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，西安 710043)

根據(jù)中國(guó)國(guó)家鐵路集團(tuán)有限公司的總體安排及規(guī)劃，川藏鐵路東部起始于四川省成都市，向西途徑雅安、康定、昌都、林芝，西止于西藏拉薩，新建線(xiàn)路總長(zhǎng)度大約為1 567.37 km[1]，川藏鐵路沿線(xiàn)地震活動(dòng)強(qiáng)烈、地震頻發(fā)，是影響工程安全的主要地質(zhì)災(zāi)害，自有記載以來(lái)，至少發(fā)生了7級(jí)以上地震22次以上[2].國(guó)內(nèi)外諸多案例表明，穿越活動(dòng)斷裂帶隧道極其容易產(chǎn)生嚴(yán)重的破壞[3-5].活動(dòng)斷層的運(yùn)動(dòng)分為粘滑運(yùn)動(dòng)和蠕滑運(yùn)動(dòng)，活動(dòng)斷層從粘結(jié)到瞬間發(fā)生滑移的這一過(guò)程被稱(chēng)為粘滑，與其相對(duì)應(yīng)的另外一種相對(duì)運(yùn)動(dòng)較為緩慢的方式叫蠕滑，二者區(qū)別就在于蠕滑沒(méi)有從粘結(jié)到瞬間發(fā)生滑移，引起地層發(fā)生瞬時(shí)的脆性斷裂，而是以較穩(wěn)定的速率蠕動(dòng)，粘滑運(yùn)動(dòng)導(dǎo)致的脆性斷裂會(huì)在瞬間釋放巨大能量，引起地震.而川藏線(xiàn)沿線(xiàn)的有諸多活動(dòng)斷裂帶，根據(jù)目前的勘探結(jié)果顯示，多數(shù)斷裂兼具走滑和逆沖的活動(dòng)性質(zhì)[6]，且具有明顯的年平均位移速率，雖然蠕滑導(dǎo)致的兩盤(pán)相對(duì)位移較小，但隨著時(shí)間不斷累積，對(duì)隧道結(jié)構(gòu)的安全性仍有威脅.

國(guó)內(nèi)外許多學(xué)者對(duì)蠕滑斷裂帶對(duì)隧道的影響規(guī)律展開(kāi)了研究，蠕滑活動(dòng)對(duì)隧道安全的影響與斷層傾角、活動(dòng)斷裂帶寬度以及活動(dòng)斷裂帶的運(yùn)動(dòng)方式有關(guān)[7-9]，國(guó)內(nèi)外規(guī)范針對(duì)活動(dòng)斷裂的處理原則，首先考慮避讓?zhuān)瑹o(wú)法避讓時(shí)，應(yīng)選擇大角度穿越[10].而受限于地理?xiàng)l件，通麥隧道只能以小角度穿越F48嘉黎-迫龍藏布斷裂，這給隧道的設(shè)計(jì)和施工都帶來(lái)了巨大挑戰(zhàn)，因此，基于通麥隧道的初步設(shè)計(jì)資料，提出了兩種穿越活動(dòng)斷層的抗錯(cuò)斷措施，通過(guò)數(shù)值模擬計(jì)算二者的抗錯(cuò)斷效果，從而為隧道的抗錯(cuò)斷措施研究以及相似工程的設(shè)計(jì)提供參考.

1 計(jì)算模型及參數(shù)

1.1 工程概況

通麥隧道位于西藏自治區(qū)南部東段波密縣通麥鎮(zhèn)，整體上沿G318國(guó)道北側(cè)展布，與國(guó)道相距最遠(yuǎn)處約3.7 km.進(jìn)口為茶隆隆巴曲左岸，道路交通不便，出口為通麥小集鎮(zhèn).隧道洞身段不良地質(zhì)類(lèi)型較多，分別為巖爆、大變形、高地溫、活動(dòng)斷裂、高瓦斯氣體等.最新地質(zhì)資料顯示，洞身強(qiáng)烈?guī)r爆段800 m，50～60 ℃中高溫段2 580 m，二級(jí)大變形段135 m，高瓦斯氣體855 m以及近出口段穿越415 m長(zhǎng)活動(dòng)斷裂帶，隧道出口端與通麥車(chē)站相連，該區(qū)域地形地貌受到區(qū)域地質(zhì)構(gòu)造的控制，山脈總體走向與構(gòu)造跡線(xiàn)一致，總體地勢(shì)北高南低，地形起伏較大.通麥隧道地理位置如圖1所示.

圖1 通麥隧道地理位置圖示Fig.1 Location of Tongmai tunnel

隧道所在區(qū)域地處念青唐古拉山東段，為青藏高原東南部隆升及河流切割最為強(qiáng)烈的地區(qū)之一，Ⅱ類(lèi)場(chǎng)地條件下通麥隧道沿線(xiàn)50年超越概率10%地震動(dòng)加速度峰值0.20 g，地震基本烈度為Ⅷ度，在通麥隧道出口段DK1166+395～DK1167+810(415 m)以小角度通過(guò)F48嘉黎-迫龍藏布斷裂，斷裂帶寬度為415 m，中國(guó)地震局地質(zhì)研究所的GPS最新觀(guān)測(cè)結(jié)果顯示嘉黎斷裂帶現(xiàn)今活動(dòng)以右旋走滑活動(dòng)為主兼有逆沖活動(dòng)性質(zhì)，右旋走滑活動(dòng)速率約2～3 mm/a，斷層破碎帶及影響帶內(nèi)以混雜巖、碎裂巖、斷層泥礫、蝕變巖等為主.通麥隧道出口段平面布置及縱斷面圖如圖2所示.

圖2 通麥隧道穿越活動(dòng)斷裂段平面布置及縱斷面圖Fig.2 Layout and longitudinal section of Tongmai Tunnel crossing active fault

結(jié)合工程實(shí)際，川藏鐵路穿越活動(dòng)斷裂帶隧道襯砌初步設(shè)計(jì)遵循“預(yù)留空間、優(yōu)化斷面、節(jié)段設(shè)計(jì)、運(yùn)營(yíng)監(jiān)測(cè)”的設(shè)計(jì)原則.活動(dòng)斷裂帶采用活動(dòng)斷裂襯砌，并適當(dāng)往兩端延伸，活動(dòng)斷裂襯砌原則上采用近圓形襯砌結(jié)構(gòu)，根據(jù)各活動(dòng)斷裂性質(zhì)考慮在襯砌節(jié)段之間設(shè)置變形縫，斷面內(nèi)凈空預(yù)留補(bǔ)強(qiáng)空間，并采用加大預(yù)留變形量、加強(qiáng)初期支護(hù)及超前支護(hù)、破碎段周邊注漿、二襯加強(qiáng)等措施.具體到通麥隧道，主要考慮設(shè)置變形縫及在初支和二襯之間設(shè)置柔性緩沖層(0.15 m)，以及鉸接設(shè)計(jì)(節(jié)段長(zhǎng)度10 m，變形縫寬度0.15 m)，加強(qiáng)隧道對(duì)蠕滑變形的適應(yīng)性.穿越活動(dòng)斷裂帶襯砌設(shè)計(jì)示意圖如圖3所示.

圖3 穿越活動(dòng)斷裂帶襯砌設(shè)計(jì)示意圖(單位：m)Fig.3 Lining design of tunnels crossing the active fault zone (unit:m)

1.2 模型建立及計(jì)算參數(shù)的選擇

根據(jù)通麥隧道平面及縱斷面布置圖，建立三維計(jì)算模型，隧道處于半無(wú)限空間中而斷層蠕滑錯(cuò)動(dòng)是一個(gè)大范圍的地層運(yùn)動(dòng)，且斷層錯(cuò)動(dòng)涉及水平方向及垂直方向上的運(yùn)動(dòng)，因此蠕滑錯(cuò)動(dòng)量并非直接作用于隧道上而是通過(guò)地層傳遞到隧道結(jié)構(gòu)上，為了減小邊界效應(yīng)，模型尺寸取2 500 m×500 m×500 m，斷層傾角為70°，斷層寬度為415 m，上盤(pán)與斷層、下盤(pán)與斷層以及隧道襯砌與圍巖之間采用庫(kù)倫摩擦接觸，參照已有研究成果，摩擦系數(shù)取0.8，建立計(jì)算模型如圖4所示.

圖4 計(jì)算模型示意圖Fig.4 Calculation model

根據(jù)川藏鐵路穿越活動(dòng)斷裂帶襯砌初步設(shè)計(jì)資料，初支采用C30噴射混凝土，二襯采用C35鋼筋混凝土，斷裂帶范圍內(nèi)圍巖主要以壓碎巖為主，隧道洞身通過(guò)地層主要為古-中元古界念青唐古拉巖群，巖性以片麻巖為主，局部分布有花崗巖、閃長(zhǎng)巖，云母片巖、綠簾石巖等透鏡體；斷裂帶內(nèi)發(fā)育壓碎巖等構(gòu)造巖.在隧道進(jìn)出口斜坡及洞身溝谷內(nèi)分布有第四系地層，巖性主要為碎石類(lèi)土，涉及冰磧、洪積、坡積、沖積等多種成因，地質(zhì)鉆孔圖如圖5所示.根據(jù)勘測(cè)資料結(jié)合含水率測(cè)定、粘聚力測(cè)定、內(nèi)摩擦角測(cè)定等室內(nèi)土工試驗(yàn)結(jié)果，得到隧道圍巖物理力學(xué)參數(shù)如表1所列，圍巖采用Mohr-Coulomb本構(gòu)，重點(diǎn)討論設(shè)置變形縫及柔性緩沖層對(duì)隧道抗錯(cuò)斷的提升效果.初支、二襯、變形縫填充材料及柔性緩沖層材料采用線(xiàn)彈性本構(gòu)，計(jì)算參數(shù)如表2所列.

圖5 地質(zhì)鉆孔圖Fig.5 Geological drilling core

表1 圍巖計(jì)算參數(shù)表Tab.1 Calculation parameters of the ground layers

表2 隧道結(jié)構(gòu)計(jì)算參數(shù)表Tab.2 Calculation parameters of tunnel structure

1.3 斷層蠕滑的模擬

斷層蠕滑產(chǎn)生的累計(jì)永久位移是導(dǎo)致隧道結(jié)構(gòu)破壞的主因，因此在計(jì)算蠕滑作用下隧道的響應(yīng)時(shí)，將其視為擬靜力問(wèn)題，在計(jì)算時(shí)，首先平衡地應(yīng)力，初始應(yīng)力由隧道埋深及上覆土體確定，然后進(jìn)行隧道的開(kāi)挖、支護(hù)過(guò)程模擬.斷層活動(dòng)是一個(gè)大范圍的運(yùn)動(dòng)，因此在模型底部和側(cè)面施加彈簧，模擬周?chē)貙雍湍Ｐ蛥^(qū)域地層的相互影響關(guān)系[13]，如圖6(a)所示，基床系數(shù)按照圍巖分級(jí)參照規(guī)范[14]中要求進(jìn)行取值.模擬斷層蠕滑時(shí)，固定模型左右兩側(cè)邊界，在上盤(pán)側(cè)面施加沿著斷層與上盤(pán)交接面方向的水平位錯(cuò)，如圖6(b)所示，該斷層的百年蠕滑量為0.2～0.3 m，計(jì)算時(shí)取0.3 m作為施加的蠕滑量.

圖6 模型邊界及位錯(cuò)施加示意圖Fig.6 Model boundary and dislocation application

1.4 工況設(shè)置

主要目的為驗(yàn)證所提出的設(shè)置變形縫以及設(shè)置柔性緩沖層對(duì)隧道抗蠕滑錯(cuò)斷的提升效果，因此主要設(shè)置了4個(gè)工況進(jìn)行對(duì)比分析，分別為無(wú)措施工況、設(shè)置變形縫工況、設(shè)置柔性緩沖層工況以及同時(shí)設(shè)置變形縫及柔性緩沖層工況，考慮到一般的混凝土二襯臺(tái)車(chē)長(zhǎng)度為6～12 m，因此本次計(jì)算每隔10 m設(shè)置一道寬度為15 cm的變形縫，柔性緩沖層厚度取0.3 m，各工況設(shè)置如表3所列.

表3 計(jì)算工況表Fig.3 Calculaiton cases

2 計(jì)算結(jié)果分析

2.1 隧道變形受力分析

以工況1為例，隧道沿軸向變形云圖如圖7所示，沿隧道軸向結(jié)構(gòu)三向位移如圖8所示.可以看出，隧道變形是追隨著地層變形的，處于破碎帶與上下盤(pán)交接面附近的隧道受斷層錯(cuò)動(dòng)影響最大，由不同方向位移統(tǒng)計(jì)可知，最大橫向位移為8.03 cm，最大豎向位移為5.93 cm，最大縱向位移為1.97 cm，橫向和豎向位移較大，縱向變形較小，在斷層破碎帶和上下盤(pán)交界面處，位移會(huì)出現(xiàn)明增加，這是因?yàn)樯舷卤P(pán)與斷層破碎帶地層交接面在此處，在此范圍內(nèi)地層位錯(cuò)很大.

圖7 隧道軸向變形云圖Fig.7 Deformation nephogram of tunnel in axial direction

圖8 不同位移統(tǒng)計(jì)Fig.8 Displacement of tunnel in different cases

工況1隧道結(jié)構(gòu)最小主應(yīng)力、最大主應(yīng)力云圖如圖9～10所示.隧道最小主應(yīng)力云圖顯示，在此處會(huì)產(chǎn)生明顯的受壓區(qū)域，應(yīng)力集中顯著，最小主應(yīng)力為-7.731 MPa，處于斷層破碎帶范圍內(nèi)的隧道最小主應(yīng)力絕對(duì)值比處于上下盤(pán)范圍隧道的大，說(shuō)明處于斷層破碎帶內(nèi)的隧道受到圍巖的壓力比處于上下盤(pán)范圍的隧道；而隧道最大主應(yīng)力云圖顯示處于斷層破碎帶中的隧道在拱肩-拱腰部位會(huì)出現(xiàn)受拉區(qū)，最大主應(yīng)力為1.648 MPa.

圖9 最小主應(yīng)力云圖Fig.9 Minimum principal stress nephogram

圖10 最大主應(yīng)力云圖Fig.10 Maximum principal stress nephogram

2.2 不同抗錯(cuò)斷措施效果分析

從已有工程來(lái)看，隧道穿越活動(dòng)斷裂帶最常用的應(yīng)對(duì)方法是設(shè)置柔性連接，讓隧道能夠更好的適應(yīng)地層的變形，本節(jié)重點(diǎn)討論設(shè)置變形縫及柔性緩沖墊層能否有效提高隧道的抗錯(cuò)斷效果.提取隧道拱頂節(jié)點(diǎn)各個(gè)方向位移沿縱向變化如圖11所示，可以看到，未設(shè)置變形縫時(shí)候，隧道各個(gè)方向位移變化趨勢(shì)沿軸向基本為連續(xù)變化，而當(dāng)設(shè)置變形縫時(shí)，襯砌沿軸向的變形不再是連續(xù)變化，在變形縫處可能會(huì)出現(xiàn)明顯階梯狀突變，兩段襯砌之間由變形縫連接，變形縫的設(shè)置減小了隧道的整體剛度，設(shè)置柔性緩沖層會(huì)吸納部分地層變形，將地層的部分變形抵御在襯砌之外.

圖11 沿隧道各方向位移統(tǒng)計(jì)Fig.11 Deformatio of tunnel in all cases

統(tǒng)計(jì)各工況不同方向位移極值如表4所列，各工況下同一方向的位移差距不大，這也說(shuō)明無(wú)法依靠設(shè)置變形縫和柔性緩沖材料來(lái)抵御地層的變形，斷層為水平方向錯(cuò)動(dòng)，因此橫向位移大于另外兩個(gè)方向的位移，因斷層與X、Y、Z平面均不平行，因此在豎向和縱向均有較大的位移.各工況下襯砌位移極值大小排序?yàn)椋汗r2>工況4>工況1>工況3，從原理上來(lái)看，變形縫的設(shè)置讓隧道變?yōu)楠?dú)立的節(jié)段，更容易追隨地層的變形，而設(shè)置柔性緩沖層會(huì)吸收部分變形，導(dǎo)致反應(yīng)在襯砌上的位移稍小.

表4 位移統(tǒng)計(jì)表Tab.4 Deformation of tunnel in all cases

提取各工況下，隧道的最大、最小主應(yīng)力繪制圖表如圖12所示，當(dāng)未設(shè)置變形縫和柔性緩沖材料時(shí)，隧道襯砌的最大，最小主應(yīng)力分別為1.648 MPa、-7.731 MPa，而工況2(設(shè)置變形縫)，工況3(設(shè)置柔性緩沖材料)，工況4(同時(shí)設(shè)置變形縫和柔性緩沖材料)的最大分別為1.364 MPa、1.618 MPa、1.318 MPa，最小主應(yīng)力分別為-4.033 MPa、-6.828 MPa、-3.943 MPa，與未設(shè)置抗錯(cuò)斷措施相比，工況2的最大、最小主應(yīng)力分減少了17.23%、47.83%，工況3分別減少了1.82%、11.68%，工況4分別減少了20.02%、48.99%.從計(jì)算結(jié)果可以看出，變形縫及柔性緩沖材料的設(shè)置均可以有效的減少襯砌的內(nèi)力，對(duì)比工況2及工況3計(jì)算結(jié)果，可以看出，變形縫的設(shè)置對(duì)襯砌內(nèi)力減小幅度明顯大于工況3，抗錯(cuò)斷效果顯著.

圖12 最大最小主應(yīng)力統(tǒng)計(jì)圖Fig.12 Maximum and minimum stress in all cases

3 結(jié)論

依托川藏鐵路通麥隧道，詳細(xì)的介紹了通麥隧道的地理地質(zhì)條件，以及其所穿越的F48嘉黎-迫龍藏布斷裂的情況，對(duì)隧道穿越活動(dòng)斷層時(shí)設(shè)置變形縫及設(shè)置柔性緩沖墊層的抗錯(cuò)斷效果進(jìn)行了數(shù)值仿真，得到了以下結(jié)論：

1) 蠕滑是一個(gè)長(zhǎng)期緩慢的過(guò)程，在活動(dòng)斷層與上下盤(pán)交界面附近的襯砌易產(chǎn)生位移突變，且此范圍內(nèi)會(huì)產(chǎn)生明顯的應(yīng)力集中現(xiàn)象，是穿越活動(dòng)斷層的薄弱環(huán)節(jié)，應(yīng)加強(qiáng)該部位的設(shè)計(jì).

2) 地層運(yùn)動(dòng)是一個(gè)大范圍運(yùn)動(dòng)，無(wú)法依靠?jī)H設(shè)置變形縫或設(shè)施隔振層而明顯減少活動(dòng)斷層蠕滑所造成的襯砌位移，從數(shù)值上來(lái)看，隧道位移極值排序?yàn)椋汗r2>工況4>工況1>工況3，這是因?yàn)樽冃慰p的設(shè)置減小了隧道的整體剛度，位移增大，而設(shè)置柔性緩沖層會(huì)吸納部分地層變形.

3) 對(duì)比襯砌最大、最小主應(yīng)力，設(shè)置變形縫和設(shè)置柔性緩沖墊層均可以減小襯砌應(yīng)力，而設(shè)置變形縫的抗錯(cuò)斷效果高于設(shè)置柔性緩沖墊層.