穿越富水?dāng)鄬悠扑閹淼浪綑C(jī)理分析與預(yù)防

　　摘 要 穿越山嶺隧道常遇到斷層破碎帶。由于缺乏斷層復(fù)雜巖體穩(wěn)定性機(jī)理的研究與理解，隧道施工曾導(dǎo)致多起大變形甚至塌方事故。本文使用理論分析、離散元數(shù)值模擬方法，結(jié)合現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果，對(duì)江西某高速鐵路隧道過(guò)斷層施工案例進(jìn)行綜合分析。在考慮斷層破碎帶地下水影響的前提下，開(kāi)展富水?dāng)鄬悠扑閹┕み^(guò)程中圍巖失穩(wěn)塌方的內(nèi)在機(jī)理研究。研究得出：隧道穿越富水?dāng)鄬悠扑閹r(shí)，巖體結(jié)構(gòu)面發(fā)育，且結(jié)構(gòu)面間多充填碎屑或黏土充填物。隨著施工擾動(dòng)引發(fā)結(jié)構(gòu)面發(fā)生錯(cuò)動(dòng)，流動(dòng)的地下水進(jìn)一步對(duì)結(jié)構(gòu)面形成沖刷，從而降低圍巖的完整性與力學(xué)強(qiáng)度。斷層破碎帶抗剪強(qiáng)度越低，圍巖松動(dòng)范圍越大，塌方概率與范圍越大。根據(jù)“防水排水，強(qiáng)化圍巖”的思路，本文提出了穿越斷層破碎帶防塌方的措施：設(shè)置迂回導(dǎo)坑、地層預(yù)加固、增加支護(hù)強(qiáng)度、優(yōu)化施工工法及加強(qiáng)監(jiān)控量測(cè)，可為此類(lèi)隧道工程的安全施工提供參考。

　　本文源自楊建輝; 沈愷; 周杰; 薛亞?wèn)|， 工程地質(zhì)學(xué)報(bào) 發(fā)表時(shí)間：2021-05-13《工程地質(zhì)學(xué)報(bào)》是我國(guó)工程地質(zhì)學(xué)科綜合性的高級(jí)學(xué)術(shù)期刊。《工程地質(zhì)學(xué)報(bào)》辦刊宗旨是加強(qiáng)學(xué)術(shù)交流，促進(jìn)工程地質(zhì)科學(xué)的了理論、應(yīng)用和技術(shù)的發(fā)展，使工程地質(zhì)學(xué)科更好地為國(guó)民經(jīng)濟(jì)建設(shè)服務(wù)。主管單位：中國(guó)科學(xué)院。

　　關(guān)鍵詞 鐵路隧道;塌方;離散元;機(jī)理分析;斷層破碎帶

　　0 引 言

　　近年來(lái)，隨鐵路隧道建設(shè)規(guī)模不斷增大，隧道施工遭遇軟弱圍巖、巖溶、斷層破碎帶等不良地質(zhì)的情況越來(lái)越多見(jiàn)，致使突水涌泥、隧道塌方等安全事故頻發(fā)(吳學(xué)智等，2016;魏雪斐等，2019)。大斷面隧道不利于承受荷載，穿越斷層破碎帶時(shí)易發(fā)生塌方事故，因此，對(duì)大斷面隧道塌方機(jī)理開(kāi)展研究具有重要的工程意義。

　　目前，國(guó)內(nèi)外學(xué)者將隧道塌方影響因素歸納為地質(zhì)條件、地下水影響、降雨影響及設(shè)計(jì)施工影響(Mete Kun et al.，2013;豐明海等，2014;尚彥軍等，2018;楊成忠等，2018)，采用模型試驗(yàn)和數(shù)值分析等研究方法研究隧道塌方機(jī)理。數(shù)值分析方法種類(lèi)繁多，目前適用于隧道塌方研究的主要有有限差分法、有限單元法、離散元法及非連續(xù)變形法。吳強(qiáng)(2009)等通過(guò) FLAC3D 研究隧道塌方形成過(guò)程中圍巖變形特征和初支受力情況。劉利生等(2019)通過(guò) FLAC3D 和現(xiàn)場(chǎng)監(jiān)測(cè)研究了隧道塌方的力學(xué)機(jī)制。徐前衛(wèi)等(2016)利用數(shù)值模擬和室內(nèi)試驗(yàn)相結(jié)合的方法，認(rèn)為斷層破碎帶會(huì)降低隧道圍巖穩(wěn)定性。Zhong et al.(2020)通過(guò) Abaqus 軟件模擬了斷層破碎帶走滑運(yùn)動(dòng)及斷層傾角對(duì)隧道的影響。基于連續(xù)介質(zhì)力學(xué)的有限差分法和有限元法在模擬隧道塌方過(guò)程存在局限性，而離散元法由于在模擬巖石斷裂過(guò)程中不需要對(duì)單元重新劃分，能更好的模擬非連續(xù)介質(zhì)大變形運(yùn)動(dòng)(薛亞?wèn)|等，2020)。汪成兵(2008)等通過(guò) PFC2D 模擬了隧道塌方全過(guò)程，提出大斷面隧道高跨比越小，隧道塌方程度越嚴(yán)重。 Huang et al.(2020)運(yùn)用 PFC2D 對(duì)比了斷層比鄰隧道和斷層穿越隧道兩種情況，認(rèn)為斷層穿越隧道時(shí)，對(duì)隧道影響更為顯著。高峰等(2018)采用 UDEC 離散元法分別模擬了不同埋深、不同巖石結(jié)構(gòu)和不同圍巖級(jí)別隧道的塌方過(guò)程，認(rèn)為圍巖等級(jí)越高，塌方范圍越大。Vazaios et al.(2019)通過(guò)有限離散元法討論了天然節(jié)理的空間分布對(duì)開(kāi)挖時(shí)圍巖穩(wěn)定性的影響。

　　在塌方處治措施研究中，張學(xué)文(2018)采用真空降水、雙層大管棚等措施解決了富水粉細(xì)砂塌體的圍巖二次變形、開(kāi)挖困難等問(wèn)題。李志厚等(2008)采用了洞內(nèi)+洞外綜合處治措施，以期從根本上預(yù)防塌方事故的發(fā)生。Zhang et al.(2019)通過(guò)注漿、排水等措施，對(duì)坍塌地層進(jìn)行快速修復(fù)，并對(duì)坍塌段進(jìn)行了現(xiàn)場(chǎng)監(jiān)測(cè)，結(jié)果表明這些措施能有效地穩(wěn)定塌方段地層。Wang et al.(2019)采取了改善圍巖應(yīng)力狀態(tài)、提高支護(hù)承載力及調(diào)整開(kāi)挖方式的方法處理塌方段施工。

　　目前，模擬隧道塌方過(guò)程的離散元分析以二維模型為主，不能直觀、系統(tǒng)地研究穿越斷層破碎帶時(shí)隧道塌方空間規(guī)律。本文針對(duì)江西省某鐵路隧道塌方事故，對(duì)塌方影響因素進(jìn)行系統(tǒng)分析，并結(jié)合現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)研究塌方前隧道圍巖變形特征，同時(shí)采用三維接觸算法離散元軟件 MatDEM 模擬穿越斷層破碎帶時(shí)塌方全過(guò)程，在塌方機(jī)理的研究基礎(chǔ)上，提出了針對(duì)性預(yù)防措施。

　　1 工程概況

　　江西省某鐵路隧道為雙線隧道，全長(zhǎng) 10 240.225m，最大埋深約 580m。隧道地質(zhì)構(gòu)造及水文地質(zhì)條件復(fù)雜，以變質(zhì)砂巖、花崗巖、砂巖、石英砂巖等地層為主，圍巖級(jí)別為Ⅱ~Ⅴ級(jí)，共發(fā)育 11 處斷裂帶。為便于表述，取掌子面里程為 K0+000，從洞內(nèi)向洞外排定相對(duì)空間位置。

　　1.1 F8-1 斷層破碎帶地質(zhì)概況

　　如圖 1 所示，隧道穿越 F8-1 斷層破碎帶長(zhǎng)度約 65m，起訖里程 K0+012~+077，產(chǎn)狀 66°∠40°。 F8 斷層破碎帶為壓性斷裂，受旋扭作用影響，不同部位斷裂面傾向不一，性質(zhì)互異。由圖 2 可見(jiàn)，掌子面揭示巖性為變質(zhì)砂巖，灰黑色，局部褐黃色，全風(fēng)化夾強(qiáng)風(fēng)化，節(jié)理裂隙發(fā)育，巖體破碎，結(jié)構(gòu)面間夾雜泥沙，為Ⅴ級(jí)圍巖。根據(jù)地勘報(bào)告顯示，斷層破碎帶富含承壓水，洞身水壓達(dá) 0.5MPa，最大涌水量為 5998.04m3 /d，為強(qiáng)富水區(qū)。

　　1.2 施工工況

　　塌方段 K0+057~K0+74 施工采用三臺(tái)階法，開(kāi)挖跨度 14.96m，開(kāi)挖高度 12.64m，斷面面積 158.25m2，屬大斷面隧道。預(yù)支護(hù)采用φ89 洞身長(zhǎng)管棚注漿，初期支護(hù)全環(huán)采用 HW175 型鋼，縱向間距 0.6m，采用φ22 連接筋連接;拱墻設(shè)置φ6 鋼筋網(wǎng)，網(wǎng)格間距 20×20cm;拱部采用φ22 組合中空錨桿，長(zhǎng)度 4m，縱向和橫向間距均為 1.5m;噴 C25 混凝土，厚度 28cm;邊墻采用φ22 砂漿錨桿，長(zhǎng)度 4m，縱向和橫向間距分別為 1.2m 和 1.0m，隧道斷面與支護(hù)參數(shù)如圖 3 所示。

　　1.3 塌方過(guò)程

　　2019 年 10 月 11 日上臺(tái)階開(kāi)挖至 K0+077 時(shí)拱腳位置發(fā)生涌水，并在 40min 內(nèi)淹沒(méi)至仰拱面，初始涌水量約 460m3 /h，之后涌水量逐步減小，最終穩(wěn)定在 25m3 /h。14 日，揭示掌子面圍巖為變質(zhì)砂巖，涌水量穩(wěn)定在 50～70m3 /h，水質(zhì)渾濁，見(jiàn)圖 4。10 月 31 日，對(duì) K0+063.6 進(jìn)行開(kāi)挖。12 月 6 日，K0+061 拱頂至左側(cè)邊墻初支面上出現(xiàn)一條環(huán)向細(xì)裂紋。為防止初支進(jìn)一步變形開(kāi)裂，在 K0+060~065 中臺(tái)階邊墻位置采用錨桿鉆機(jī)打設(shè)φ89 管進(jìn)行初支加固，鋼管間距為 60cm，左右各打 8 根，打入長(zhǎng)度 5m，鋼管采用 L 筋與初支鋼架焊接，并用槽鋼與各根鋼管縱向焊接為同一結(jié)構(gòu)整體。

　　12 月 21 日對(duì) K0+072~K0+077 處仰拱基底進(jìn)行注漿加固，加固深度 5~9 m。12 月 27 日，在對(duì) K0+000 上臺(tái)階初支進(jìn)行局部補(bǔ)噴時(shí)， K0+063.6 處拱頂及左側(cè)邊墻發(fā)生坍塌(見(jiàn)圖 5)。塌方段埋深約 110.21m，如圖 6 所示，拱頂塌方范圍約為 15m，坍體體積約 2000m3，預(yù)估塌腔體為帽型。

　　2 塌方前隧道圍巖變形監(jiān)測(cè)

　　選取了塌方段前后 7 個(gè)監(jiān)測(cè)斷面(K0+072、 K0+067、K0+062、K0+057、K0+052、K0+047 及 K0+042)，時(shí)間從 K0+072 初支施作完畢(2019 年 10 月 29 日)至塌方事故發(fā)生當(dāng)天(2019 年 12 月 27 日)，共計(jì) 60d，監(jiān)測(cè)頻率為 1 次/ d。圍巖變形時(shí)程曲線見(jiàn)圖 7。

　　由圖 7 可見(jiàn)，截止塌方當(dāng)天，7 個(gè)監(jiān)測(cè)斷面均未表現(xiàn)出收斂的趨勢(shì)，曲線以線性形式增長(zhǎng)，說(shuō)明隧道圍巖受富水?dāng)鄬悠扑閹в绊懀瑖鷰r長(zhǎng)時(shí)間仍較難達(dá)到穩(wěn)定。拱頂下沉最大值 57.37mm 出現(xiàn)在 K0+072 處，拱頂下沉速率最大值 3.50mm/d 出現(xiàn)在 K0+062 處，邊墻收斂最大值 54.96mm 出現(xiàn)在 K0+062 處，邊墻收斂速率最大值 3.86mm/d 出現(xiàn)在 K0+052 處，表明塌方與圍巖變形存在相關(guān)性，收斂量大且收斂速率快的段落塌方風(fēng)險(xiǎn)較高。K0+072 處距離涌水點(diǎn)位置最近，拱頂下沉量最大，但卻并未塌方，塌方發(fā)生在 K0+63.6 處，說(shuō)明塌方事故產(chǎn)生的原因不僅與斷層破碎帶、地下水等地質(zhì)因素有關(guān)，也與施工、支護(hù)等因素有關(guān)。K0+062 處初期支護(hù)相較于其他斷面支護(hù)薄弱，存在施工不到位的可能，加之上臺(tái)階 K0+57.8 處開(kāi)挖爆破振動(dòng)及機(jī)械振動(dòng)的影響，誘發(fā)了本次塌方事故。

　　3 MatDEM 離散元模擬

　　MatDEM 是采用 GPU 矩陣計(jì)算法和三維接觸算法的離散元軟件(Liu 等，2017 朱晨光等，2019)，適用于模擬基于離散介質(zhì)力學(xué)法的非連續(xù)材料在靜載或動(dòng)載作用下的力學(xué)響應(yīng)過(guò)程。本研究通過(guò) MatDEM 構(gòu)建三維隧道施工模型利用顆粒間彈簧的法向剛度、切向剛度及抗剪強(qiáng)度來(lái)模擬結(jié)構(gòu)面，進(jìn)而模擬出隧道塌方全過(guò)程，分析塌方機(jī)理。

　　3.1 模型參數(shù)及 3D 建模

　　根據(jù)圖 1 隧道塌方段，結(jié)合地質(zhì)勘測(cè)資料，運(yùn)用 MatDEM 離散元軟件建立 110m×60m×55m 的全尺寸隧道模型，隧道模型見(jiàn)圖 8。模型中，顆粒半徑為 0.6m，總共 267794 個(gè)顆粒，四周邊界采用固定約束，上邊界采用應(yīng)力約束，應(yīng)力取 1MPa。支護(hù)方案基于文獻(xiàn)(谷拴成等，2018)換算為支護(hù)與圍巖等效錨固體并根據(jù)實(shí)際塌方情況進(jìn)行了一定修正，圍巖與錨固體宏觀參數(shù)由室內(nèi)試驗(yàn)及根據(jù)設(shè)計(jì)資料確定，微觀參數(shù)通過(guò)轉(zhuǎn)換公式(劉春等，2019，見(jiàn)附錄 A)由宏觀參數(shù)計(jì)算得出，具體參數(shù)見(jiàn)表 1，根據(jù)材料性質(zhì)，模擬時(shí)間步確定為 5.0×10-5 s。

　　3.2 塌方過(guò)程模擬

　　隧道開(kāi)挖前后圍巖的應(yīng)力狀態(tài)見(jiàn)圖 9，可以直觀得出：隧道開(kāi)挖前，斷層破碎帶與地層處于穩(wěn)定狀態(tài)，隧道開(kāi)挖后，圍巖應(yīng)力二次重分布，斷層破碎帶及隧道周邊出現(xiàn)了明顯的應(yīng)力集中。開(kāi)挖后相較于開(kāi)挖前，斷層破碎帶邊界所受剪應(yīng)力增大，隧道圍巖壓應(yīng)力減小，而隧道圍巖處于臨空狀態(tài)，導(dǎo)致巖塊之間結(jié)構(gòu)面法向力減小。由圖 10 可見(jiàn)，當(dāng)隧道開(kāi)挖至斷層破碎帶處時(shí)，拱頂和邊墻產(chǎn)生了應(yīng)力集中，且左側(cè)邊墻應(yīng)力大于右側(cè)，與隧道在拱頂與左側(cè)邊墻處發(fā)生塌方的事實(shí)吻合。

　　由圖 11(a)可見(jiàn)，隧道開(kāi)挖后，斷層破碎帶內(nèi)的軟弱巖體產(chǎn)生了大范圍松動(dòng)，圍巖位移最大值 1.12m 位于掌子面前方 30m 處，和現(xiàn)場(chǎng)塌方位置接近。上臺(tái)階和中臺(tái)階產(chǎn)生了一定的隆起，因此在掌子面開(kāi)挖后應(yīng)盡早施作初支，且閉合成環(huán)。由圖 11 (b)可見(jiàn)，拱肩上部巖體位移較大，拱腳處位移較小，是隧道塌方后圍巖應(yīng)力釋放導(dǎo)致。當(dāng)顆粒間彈簧作用力為 0 時(shí)，可視作塌方，經(jīng)計(jì)算塌方量為 2640m3，與現(xiàn)場(chǎng)坍塌體體積基本吻合。從隧道橫剖面和縱剖面變形圖可以看出，圍巖最大松動(dòng)范圍在 13~17m，這與預(yù)估塌腔范圍相符，并且可以觀察到塌落拱的大致形狀。

　　3.3 富水?dāng)鄬訋M分析

　　基于現(xiàn)場(chǎng)調(diào)查與理論分析，地下水對(duì)斷層帶的影響主要表現(xiàn)為斷層破碎帶內(nèi)巖體結(jié)構(gòu)面抗剪強(qiáng)度的弱化，考慮兩種極端抗剪強(qiáng)度(0.1 倍，2 倍)條件下，計(jì)算圍巖最大位移量和塌方量，以研究地下水對(duì)隧道開(kāi)挖帶來(lái)的影響。由圖 13 可見(jiàn)，圍巖松動(dòng)范圍受結(jié)構(gòu)面抗剪強(qiáng)度影響，取 0.1 倍抗剪強(qiáng)度時(shí)，圍巖最大松動(dòng)范圍在 20 ~23m，塌方量為 13500m3，圍巖最大位移 2.51m;采用實(shí)際抗剪強(qiáng)度時(shí)，圍巖最大松動(dòng)范圍在 13~17m，塌方量 2640 m3，圍巖最大位移 1.12m;取 2 倍抗剪強(qiáng)度時(shí)，圍巖最大松動(dòng)范圍在 4~7m，塌方量為 16.7m3，圍巖最大位移 0.80m。表明富水?dāng)鄬悠扑閹е薪Y(jié)構(gòu)面抗剪強(qiáng)度低，圍巖松動(dòng)范圍大，出現(xiàn)圍巖失穩(wěn)塌方的概率與規(guī)模越大。

　　4 塌方機(jī)理分析

　　斷層破碎帶節(jié)理裂隙發(fā)育，巖體強(qiáng)度主要受結(jié)構(gòu)面影響。結(jié)構(gòu)面組數(shù)多，軟、硬結(jié)構(gòu)面混雜且排列無(wú)序?qū)е聰鄬悠扑閹r體強(qiáng)度低，受地下水軟化、泥化和潤(rùn)滑作用的進(jìn)一步影響易發(fā)生塌方事故。地質(zhì)構(gòu)造及地下水因素是隧道塌方的主要因素及客觀因素，支護(hù)及施工措施不到位是隧道塌方的次要因素及主觀因素。通過(guò) MatDEM 數(shù)值模擬計(jì)算結(jié)果，以及現(xiàn)場(chǎng)勘測(cè)的實(shí)際情況，對(duì)隧道塌方機(jī)理展開(kāi)分析。

　　4.1 斷層破碎帶影響

　　斷層破碎帶結(jié)構(gòu)面的性質(zhì)主要由法向變形、剪切變形及抗剪強(qiáng)度三方面組成。隧道開(kāi)挖后，受斷層破碎帶影響的淺部巖體出現(xiàn)應(yīng)力集中現(xiàn)象，巖體產(chǎn)生兩種破壞形式：巖體沿結(jié)構(gòu)面法向張拉破壞、巖體沿結(jié)構(gòu)面切向剪切破壞。為研究結(jié)構(gòu)面法向變形能力和剪切變形能力，Goodman 提出了法向變形剛度