陶連金， 殷福鸞， 安 韶

(北京工業(yè)大學(xué) 城市與工程安全減災(zāi)省部共建教育部重點實驗室， 北京 100124)

0 引 言

國內(nèi)外學(xué)者對隧道穿越斷層方面做了較為詳細的研究。張偉喜等[2]針對活動正斷層黏滑錯動系統(tǒng)計算分析了整體式隧道襯砌結(jié)構(gòu)縱向應(yīng)變、主應(yīng)力、裂縫和圍巖-初期支護附加接觸力等分布規(guī)律。趙寶平[3]針對不同斷層錯距分析逆斷層錯動作用下隧道襯砌受力與變形規(guī)律得出襯砌破壞形式。熊煒等[4]研究了斷層錯動量、斷層傾角、隧道埋深以及隧道與斷層的交角等因素對正斷層公路山嶺隧道襯砌受力變形的影響。趙伯明等[5]對跨斷層隧道在斷層發(fā)生錯動時的響應(yīng)機理進行數(shù)值模擬，對比分析在不同斷層傾角、斷層寬度及圍巖條件下，隧道的位移響應(yīng)、結(jié)構(gòu)變形和襯砌內(nèi)力變化規(guī)律，以及隧道結(jié)構(gòu)的塑性區(qū)開展情況。孫飛等[6-7]通過數(shù)值模擬、變形監(jiān)測和模型實驗方法，分析了斷層錯動下，正斷層塑性區(qū)分布規(guī)律、裂縫分布特征及隧道結(jié)構(gòu)在斷層剪切錯動下的破壞范圍及破壞形態(tài)。焦鵬飛等[8]采用FLAC3D有限差分軟件,建立數(shù)值模型,分析由于地震引發(fā)的逆斷層錯動作用對正交穿越斷層隧道的影響,并揭示其影響機理。陳海亮等[9]以現(xiàn)場監(jiān)測與數(shù)值模擬相結(jié)合的手段，對斷層寬度20 m、傾角70°,以及斷層寬度8 m、傾角65°的跨斷層隧道變形及破壞機理進行分析，探索了斷層錯動對跨斷層隧道的影響及破壞特征。從前人研究可以看出，深埋山嶺隧道襯砌受力變形方面已取得不少成果，但淺埋地鐵隧道跨活動斷裂方面的研究比較少。筆者在前人對斷層調(diào)查和研究基礎(chǔ)上，采用有限元軟件ABAQUS，考慮了斷層錯動距離和隧道埋深兩個影響因素，模擬分析斷層錯動作用下，兩種不同因素條件下隧道襯砌結(jié)構(gòu)的受力與變形，以期為相關(guān)工程提供一定的指導(dǎo)。

1 工程概況

烏魯木齊地鐵2號線某隧道區(qū)間大角度通過西山斷層北支，見圖1。該斷層形成于中更新世中晚期，最新活動時間為晚更新世晚期，屬晚更新世活動斷層，走向N45°～70°E，傾向N，傾角44°～ 83°，具逆沖性質(zhì)。由于人類活動(如抽取地下水)和全新世構(gòu)造應(yīng)力的影響，目前，活動斷層每年的垂直速率達到1～5 mm/a。地質(zhì)剖面見圖2。

圖1 西山斷層?xùn)|段南北支示意

圖2 地質(zhì)剖面

2 模型建立

2.1 計算方案

取斷層上下盤各100 m建立計算模型，模型包括斷層破碎帶、上盤、下盤、上覆土層和襯砌5個部分，均采用實體單元，其中斷層破碎帶為煤層。模型長×寬×高分別為252 m×68 m×60 m，斷層傾角為45°，隧道直徑為6.7 m，初期支護和二次襯砌厚度分別為350和250 mm。計算模型幾何尺寸和網(wǎng)格劃分見圖3(以隧道埋深12 m為例)，土體和襯砌網(wǎng)格劃分均采用C3D8R單元。

圖3 模型尺寸與網(wǎng)格劃分

根據(jù)實際工程，隧道結(jié)構(gòu)建于上覆土層中，將該次計算斷層錯距取為5.0、10.0、16.5、20.0 cm，隧道埋深取為12、10和8 m，計算方案見表1。

表1 計算方案

2.2 計算參數(shù)

數(shù)值模擬中圍巖和斷層破碎帶采用彈塑性模型，摩爾庫倫屈服準(zhǔn)則，襯砌為復(fù)合式襯砌，采用彈性模型模擬，外襯與土體采用摩擦接觸，摩擦系數(shù)為0.4[10]。模擬分為3步：初始地應(yīng)力平衡—隧道開挖施加襯砌—施加斷層錯動位移。在前兩個分析步中，模型底部及側(cè)邊界施加法向約束，上部邊界為自由邊界，在最后一個分析步中，由于逆斷層錯動為上盤相對下盤作向上運動，因此為實現(xiàn)該錯動過程，釋放上盤底部及側(cè)邊界法向約束，并在相應(yīng)的位置上施加位移來模擬斷層錯動，位移加載邊界條件見圖4。

圖4 位移加載模擬

圍巖力學(xué)參數(shù)參考所依托隧道工程地質(zhì)勘測資料，計算參數(shù)如表2。

表2 計算參數(shù)

3 計算結(jié)果與分析

3.1 不同斷層錯距對襯砌受力變形的影響

3.1.1 位移分析

兩組患者VAS評分結(jié)果見表3，兩組患者VAS評分均隨時間延長而顯著下降，不同時間點間差異有統(tǒng)計學(xué)意義（P＜0.05）。但相應(yīng)時間點早期組與晚期組VAS評分差異均無統(tǒng)計學(xué)意義（P＞0.05））。末次隨訪所有患者結(jié)核中毒癥狀及局部疼痛均消失，結(jié)核病灶均無局部復(fù)發(fā)。

圖5是隧道埋深12 m，逆斷層錯動16.5 cm時，放大50倍的地層豎向位移云圖和襯砌的豎向位移云圖。從圖5可見，上盤受錯動位移的作用，使上盤土體產(chǎn)生向上的位移16.5 cm，隧道襯砌沿著縱向發(fā)生了“S”狀彎曲，襯砌產(chǎn)生向上的位移16.5 cm。

圖5 豎向位移云圖

圖6是逆斷層錯動不同錯距時襯砌拱頂處的豎向位移曲線。圖6中，以模型下盤邊界處為0點，往右側(cè)指向上盤方向為正，破碎帶兩側(cè)邊界分別在縱向坐標(biāo)100和152 m處，斷層錯動發(fā)生在縱向坐標(biāo)152 m處，縱坐標(biāo)的豎向位移s，以豎直向上為正，向下為負。隨著縱向距離dy的減小，襯砌豎向位移逐漸減小。隨著斷層錯動距離的增加，襯砌豎向位移逐漸增大。

圖6 拱頂豎向位移曲線

襯砌的豎向位移分為3個階段，整體上升段(一階段)、差異變形段(二階段)、漸趨平緩段(三階段)。斷層錯動作用下，襯砌豎向位移整體上升段發(fā)生在上盤區(qū)段；差異變形段發(fā)生在斷層錯動面附近一定范圍，在該范圍內(nèi)，襯砌豎向位移變化差異明顯；漸趨平緩段發(fā)生在斷層下盤，襯砌豎向位移逐漸趨于平緩。

錯動5 cm時，襯砌漸趨平緩段范圍沿縱向0～112 m，差異變形段112～168 m，襯砌整體上升段168～252 m；錯動10 cm時，襯砌漸趨平緩段范圍沿縱向0～108 m，差異變形段108～168 m，襯砌整體上升段168～252 m；錯動16.5 cm時，襯砌漸趨平緩段范圍沿縱向0～112 m，差異變形段112～168 m，襯砌整體上升段168～252 m；錯動20 cm時，襯砌漸趨平緩段范圍沿縱向0～110 m，差異變形段110～168 m，襯砌整體上升段168～252 m。總體來看，襯砌變形的影響范圍(差異變形段)幾乎不受斷層錯距的影響。

3.1.2 結(jié)構(gòu)內(nèi)力分析

圖7為斷層錯動作用下襯砌內(nèi)力分布曲線。

從圖7可以看出，襯砌彎矩沿縱向距離呈“S”形變化，且正負彎矩交界點在斷層錯動面處，兩側(cè)大致呈現(xiàn)反對稱分布，襯砌的彎曲變形集中在該范圍。當(dāng)斷層錯動5 cm時，正負彎矩最大值分別出現(xiàn)在沿縱向140、160 m處，分別為-126.3、152.7 MN·m；斷層錯動10.0 cm時，正負彎矩最大值分別出現(xiàn)在136、160 m處，分別為-229.5、285.8 MN·m；斷層錯動16.5 cm時，正負彎矩最大值分別出現(xiàn)在136、160 m處，分別為-352.5、428 MN·m；斷層錯動20.0 cm時，正負彎矩最大值分別出現(xiàn)在132、160 m處，分別為-411.8、494.8 MN·m?？梢钥闯?，隨著斷層錯距的增加，襯砌彎矩M逐漸增大但變化趨勢相同。

當(dāng)斷層錯動5.0、10.0、16.5、20.0 cm時，剪力FQ最大值分別為24.2、40.8、55.3、60.8 MN，軸力FN最大值分別為-90、-165.7、-211.1、-218.2 MN。襯砌剪力和軸力在斷層錯動面處達到最大，向兩側(cè)逐漸減小，隨著斷層錯距的增大，襯砌軸力和剪力逐漸增大。

3.1.3 襯砌縱向應(yīng)力分析

縱向應(yīng)力采用建立圓上三節(jié)點柱局部坐標(biāo)系提取，局部坐標(biāo)系示意見圖8。

圖8 局部坐標(biāo)系

圖9為隧道埋深12 m，斷層錯動16.5 cm下，襯砌拱頂、拱底縱向應(yīng)力分布云圖，其中拉應(yīng)力為正，壓應(yīng)力為負。

圖9 襯砌縱向應(yīng)力云圖

圖10為不同斷層錯動距離下，襯砌拱頂和拱底縱向應(yīng)力變化曲線。由圖10可見,隨著斷層錯動距離的增加，襯砌拱頂和拱底縱向應(yīng)力逐漸增大。

圖10 不同斷層錯距時襯砌的縱向應(yīng)力變化曲線

由圖10可見，拱頂和拱底的壓應(yīng)力均大于拉應(yīng)力，拉、壓應(yīng)力變化點出現(xiàn)在斷層錯動面處，襯砌的拉、壓應(yīng)力集中分布在斷層錯動面兩側(cè)，向兩側(cè)逐漸減小，由于周圍圍巖的作用，遠離斷層破碎帶的區(qū)域，襯砌拱頂和拱底均呈現(xiàn)一定的壓應(yīng)力。斷層錯動5.0、10.0、16.5、20.0 cm時，襯砌拱頂壓應(yīng)力最大值出現(xiàn)在沿縱向140、136、136和136 m處，分別為-14.8、-28.2、-40.2、-44.6 MPa；拱底壓應(yīng)力最大值均出現(xiàn)在沿縱向168 m處，分別為-17.6、-32.4、-45.0、-49.0 MPa；襯砌拱頂拉應(yīng)力最大值均出現(xiàn)在沿縱向160 m處，分別為2.8、5.3、11.5、15.8 MPa；拱底拉應(yīng)力最大值出現(xiàn)在沿縱向148、144、140、140 m處，分別為0.7、1.3、7.7、8.0 MPa。可見隨著斷層錯距的增加，襯砌縱向應(yīng)力逐漸增大。

3.2 不同隧道埋深對襯砌受力變形的影響

3.2.1 結(jié)構(gòu)內(nèi)力分析

圖11為斷層錯動16.5 cm作用下不同埋深隧道襯砌沿軸向的內(nèi)力分布曲線。由圖11可以看出，隧道埋深12、10、8 m時，襯砌彎矩最大值出現(xiàn)在沿縱向160、164、164 m處,分別為428.0、376.1、340.4 MN·m。隨著隧道埋深的減小，襯砌彎矩最大值逐漸減小，但總體變化趨勢相同。隧道埋深12、10、8 m時，襯砌剪力最大值分別為55.3、41.4、33.3 MN，軸力最大值分別為-211.1、-176.7、-156.1 MN，襯砌剪力和軸力均在斷層錯動面處達到最大值，隨著隧道埋深的減小，襯砌剪力和軸力逐漸減小。

圖11 不同隧道埋深時襯砌的內(nèi)力分布

綜上，隨著隧道遠離斷層，襯砌對斷層錯動的響應(yīng)逐步減弱。遠離斷層帶一定距離后，襯砌內(nèi)力變化趨于平緩，收斂于一個穩(wěn)定的較小值。

3.2.2 襯砌縱向應(yīng)力分析

圖12為不同埋深條件下襯砌縱向應(yīng)力變化曲線。

圖12 不同隧道埋深時襯砌的縱向應(yīng)力變化曲線

由圖12可見，埋深12、10、8 m時，襯砌拱頂壓應(yīng)力最大值均出現(xiàn)在沿縱向136 m處，分別為-40.2、-37.6、-35.0 MPa；拱底壓應(yīng)力最大值均出現(xiàn)在沿縱向168 m處，分別為-45.0、-37.1、-33.0 MPa；襯砌拱頂拉應(yīng)力最大值均出現(xiàn)在沿縱向160 m處，分別為11.5、11.0、10.6 MPa；拱底拉應(yīng)力最大值均出現(xiàn)在沿縱向140 m處，分別為7.7、7.5、7.4 MPa。隨著隧道埋深的減小，襯砌拉壓應(yīng)力最大值逐漸減小，但襯砌縱向應(yīng)力總體變化趨勢不變。

4 結(jié) 論

(1)襯砌拱頂豎向位移在斷層錯動面附近呈三區(qū)段“S”形變化，襯砌的影響范圍(差異變形段)幾乎不受錯距的影響。

(2)隨著斷層錯距的增加，地層釋放的能量增加，襯砌豎向位移、內(nèi)力和縱向應(yīng)力均逐漸增大。

(3)隨著隧道埋深的減小，隧道底部距斷層錯動面距離越大，土層耗散能量越多，襯砌內(nèi)力和縱向應(yīng)力最大值逐漸減小。

(4)襯砌的剪力和軸力在斷層錯動面處取得最大值，彎矩和縱向應(yīng)力最大值分布在斷層錯動面兩側(cè)。襯砌在斷層錯動面及其附近容易發(fā)生張拉擠壓型破壞和剪斷破壞，應(yīng)重點設(shè)防。