李廷春，呂連勛, 2，段會(huì)玲，陳偉, 3

?

深埋隧道穿越富水破碎帶圍巖突水機(jī)理

李廷春1，呂連勛1, 2，段會(huì)玲1，陳偉1, 3

(1. 山東科技大學(xué)山東省土木工程防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室，山東青島，266590；2. 北京市勘察設(shè)計(jì)研究院有限公司，北京，100038；3. 濟(jì)南城建集團(tuán)有限公司，山東濟(jì)南，250000)

基于隧道穿越處于復(fù)雜應(yīng)力場(chǎng)與滲流場(chǎng)環(huán)境的富水破碎帶時(shí)存在發(fā)生重大突水事故的安全隱患，通過(guò)對(duì)破碎巖體的滲流特點(diǎn)進(jìn)行研究，建立孔隙顆粒介質(zhì)流失的滲流模型；基于連續(xù)介質(zhì)力學(xué)和變質(zhì)量動(dòng)力學(xué)理論，推導(dǎo)飽和破碎巖體變質(zhì)量滲流?變形耦合理論模型；以福建漳州梁山隧道L7富水破碎帶為工程背景，分析圍巖的滲流場(chǎng)、應(yīng)力場(chǎng)與位移場(chǎng)分布特性，并總結(jié)隧道斷層破碎帶的突水塌陷機(jī)理。研究結(jié)果表明：斷層破碎帶突水實(shí)質(zhì)上是圍巖的力學(xué)平衡和地下水的滲流平衡因施工擾動(dòng)發(fā)生急劇變化，引起圍巖應(yīng)力重分布及地下水能量釋放；隧道施工揭露斷層后，巖體顆粒隨孔隙空間的流體發(fā)生遷移形成新的滲流通道，導(dǎo)致地下水在水頭壓力作用下向工程臨空面涌出，形成漏斗形的滲水區(qū)域；隨著滲流作用時(shí)間的延長(zhǎng)，地下水和巖土體逐漸流失，隧道上方的破碎巖體發(fā)生嚴(yán)重的滑移變形，形成橢圓形塌陷區(qū)域，與現(xiàn)場(chǎng)實(shí)際塌陷破壞規(guī)律基本吻合。本文提出的滲流?變形耦合模型對(duì)理解破碎巖體滲流力學(xué)機(jī)制和深埋隧道突水災(zāi)害的預(yù)防設(shè)計(jì)具有參考價(jià)值。

深埋隧道；巖體；斷層破碎帶；突水機(jī)理；流固耦合

隨著現(xiàn)代社會(huì)資源開發(fā)與基礎(chǔ)設(shè)施建設(shè)的快速發(fā)展，越來(lái)越多的隧道工程不得不穿越或修建在軟弱破碎地層中。破碎巖體長(zhǎng)期處于復(fù)雜的應(yīng)力場(chǎng)和滲流場(chǎng)環(huán)境中，巖土體的性質(zhì)較差，受地表水和地下水的滲流作用以及開挖和支護(hù)過(guò)程的擾動(dòng)影響，極易產(chǎn)生突水涌泥、巖土體塌陷等地質(zhì)災(zāi)害問(wèn)題[1?4]。據(jù)統(tǒng)計(jì)，在隧道工程施工過(guò)程中，由突涌水災(zāi)害造成的重大安全事故達(dá)到80%，其中斷層型突水事故占30%[5]。隧道工程突涌水事故不僅會(huì)給隧道施工帶來(lái)困難、造成嚴(yán)重的工期延誤和經(jīng)濟(jì)損失，甚至?xí)T發(fā)水資源枯竭，引起地表塌陷等環(huán)境地質(zhì)災(zāi)害，造成人員傷亡，威脅到社會(huì)穩(wěn)定與經(jīng)濟(jì)發(fā)展。近年來(lái)，在大秦鐵路軍都山隧道、渝懷鐵路武隆隧道、張集鐵路舊堡隧道等工程建設(shè)過(guò)程中，因?yàn)榻衣镀扑閹?dǎo)致突水涌泥事故，造成了嚴(yán)重的經(jīng)濟(jì)損失和人員傷亡[4?7]。中國(guó)隧道建設(shè)的重心長(zhǎng)期處于地形地質(zhì)條件極端復(fù)雜的山區(qū)，一大批具有巖體破碎、高地應(yīng)力、強(qiáng)滲透壓等特點(diǎn)的深長(zhǎng)隧道正在興建[5?6]。目前，人們對(duì)有關(guān)巖溶隧道突水風(fēng)險(xiǎn)與治理措施進(jìn)行了很多研究[7?10]，但對(duì)于富水軟弱破碎帶隧道的突水機(jī)理研究較少。探究破碎巖體滲流與變形相耦合的特性，從水力耦合的角度揭示破碎巖體的滲流特性及其災(zāi)變規(guī)律，對(duì)深長(zhǎng)隧道穿越富水破碎帶的重大水害發(fā)生機(jī)理與預(yù)防設(shè)計(jì)具有重要的理論和實(shí)踐意義。

1 飽和破碎巖體滲流模型

破碎巖體孔隙空間由天然孔隙(有機(jī)質(zhì)孔和無(wú)機(jī)質(zhì)孔)、天然裂隙、水力裂隙和人工裂隙等空間組成，屬于孔隙?裂隙復(fù)合介質(zhì)，但考慮到其充填物較松散，將其進(jìn)行等效處理后視為孔隙介質(zhì)。液體以流動(dòng)狀態(tài)賦存于孔隙空間中，將巖體內(nèi)的軟弱部分逐漸轉(zhuǎn)化為流體中的懸浮顆粒，因此，破碎巖體可以看作由巖體骨架、孔隙液體以及流體中的懸浮顆粒3種介質(zhì)構(gòu)成。在斷層破碎帶、巖溶管道等不良地質(zhì)揭露后，孔隙空間的懸浮顆粒逐漸從微觀滲流通道流失，引起巖體滲透特性改變，最終誘發(fā)突水災(zāi)害。一方面，微小的固體顆粒在水力梯度作用下，隨滲流作用不斷產(chǎn)生并流失，造成巖體的孔隙率和滲透率不斷增加，使巖體骨架的變形加大，巖體的破碎程度加??；另一方面，巖體的損傷與變形加劇，導(dǎo)致孔隙空間的滲流速度加快，懸浮顆粒的流動(dòng)更加通暢，流失速率加快。在巖體變形發(fā)展和顆粒滲流流失的共同作用下，破碎巖體的性質(zhì)逐漸變差，這就是流固耦合效應(yīng)。

為了研究懸浮顆粒遷移作用下破碎巖體的流固耦合問(wèn)題，建立如圖1所示的滲流模型，并進(jìn)行以下基本假定：

1) 巖體處于完全飽和狀態(tài)，孔隙空間充滿流體介質(zhì)。

2) 流體的滲流速度緩慢，且符合Darcy律。

3) 忽略流體對(duì)巖體骨架及懸浮顆粒的溶解作用，并忽略滲流過(guò)程中顆粒與骨架、顆粒與顆粒間的碰撞。

4) 滲流過(guò)程中僅考慮顆粒在水力梯度作用下的遷移效應(yīng)，忽略顆粒因濃度梯度變化而發(fā)生的擴(kuò)散效應(yīng)。

圖1 破碎巖體滲流模型

流體中懸浮顆粒的質(zhì)量變化可以采用顆粒的質(zhì)量濃度進(jìn)行描述[11]。記巖體的總體積為，孔隙率為，孔隙液體與巖體骨架的密度分別為w和s，懸浮顆粒占流體介質(zhì)的體積分?jǐn)?shù)為，則顆粒的質(zhì)量濃度sg為

則流體介質(zhì)中懸浮顆粒的質(zhì)量濃度變化率sg為

2 巖體變質(zhì)量滲流?變形耦合模型

2.1 懸浮顆粒流失的質(zhì)量守恒方程

將巖體孔隙空間內(nèi)的滲流分解為三維坐標(biāo)方向，取如圖2所示的巖石微元體，則沿著軸方向，單位時(shí)間內(nèi)流入微元體的懸浮顆粒質(zhì)量為

式中：sx為懸浮顆粒的絕對(duì)速度在軸方向的分量。

圖2 微元體滲流示意圖

考慮流體滲流的連續(xù)性，單位時(shí)間內(nèi)沿著軸方向流出微元體的顆粒質(zhì)量m2為

那么，在方向凈流入微元體的顆粒質(zhì)量m為

同理，可得和方向凈流入的顆粒質(zhì)量分別為

由顆粒質(zhì)量濃度變化引起的微元體中顆粒的質(zhì)量變化量為

根據(jù)質(zhì)量守恒定律，有

此即飽和破碎巖體懸浮顆粒流失的質(zhì)量守恒方程，將其寫為矢量形式：

2.2 懸浮顆粒的運(yùn)動(dòng)方程

流體中懸浮顆粒的遷移可以看作是滲流作用的效果，此時(shí)液體介質(zhì)與懸浮顆粒間會(huì)產(chǎn)生相對(duì)運(yùn)動(dòng)，并發(fā)生動(dòng)量的傳遞。懸浮顆粒在流體中受復(fù)雜的作用 力[12?15]，忽略部分微小作用力對(duì)懸浮顆粒運(yùn)動(dòng)的影響，將力系簡(jiǎn)化后認(rèn)為顆粒的運(yùn)動(dòng)主要受浮力、重力和曳力的作用，懸浮顆粒的運(yùn)動(dòng)方程為

式中：s和s分別為懸浮顆粒的密度與直徑；D為拽力系數(shù)；w為孔隙液體的密度；為懸浮顆粒與液體間的相對(duì)速度。

2.3 液體介質(zhì)的運(yùn)動(dòng)方程

將破碎巖體等效為連續(xù)多孔介質(zhì)，孔隙空間的滲流滿足Darcy定律，即

式中：i為液體的滲流速度；為巖體的滲透系數(shù)(m2/(Pa·s))；為水力梯度。常密度流體在均質(zhì)各向同性巖體中滲流的情況下，液體介質(zhì)的運(yùn)動(dòng)方程具有如下形式：

式中：為孔隙水壓力；為含顆粒流體的動(dòng)力黏度。流體的動(dòng)力黏度與純液體的動(dòng)力黏度0有如下關(guān) 系[16]：

2.4 巖體骨架的運(yùn)動(dòng)方程

記微元體受到應(yīng)力為，位移為s，體積力為M，則固體介質(zhì)的動(dòng)量守恒方程為

應(yīng)力、應(yīng)變與位移之間的關(guān)系為：

式中：和分別為拉梅系數(shù)與剪切模量。將式(15)和(16)代入式(14)，可得巖體骨架的運(yùn)動(dòng)方程：

2.5 耦合作用方程

為了實(shí)現(xiàn)模型中液體滲流與巖體骨架變形的動(dòng)態(tài)平衡，還需要建立顆粒流失過(guò)程中材料性質(zhì)變化的相關(guān)方程。

易知某時(shí)刻巖體骨架微元d的質(zhì)量為。假定在時(shí)間[，+d]內(nèi)，巖體顆粒的流失質(zhì)量為，則有

化簡(jiǎn)后，得

已知破碎巖體孔隙率與滲透率之間滿足[12]

其中：0為初始滲透率；為孔隙率；0為初始孔隙率。孔隙空間中液體質(zhì)點(diǎn)的速度與滲流速度s之間滿足Dupuit?Forchheimer關(guān)系式[17]：

式(8)，(9)，(11)，(17)，(19)，(20)和(21)構(gòu)成了飽和破碎巖體中懸浮顆粒流失的變質(zhì)量滲流?變形耦合模型基本方程。

3 數(shù)值模擬分析實(shí)例

3.1 工程概況

梁山隧道位于福建省漳州市，設(shè)計(jì)為客貨兩用單洞雙線隧道，是廈深鐵路全線的控制工程。隧道全長(zhǎng)9.89 km，進(jìn)口里程為DK94+000 km，出口里程為DK103+888 km。隧址區(qū)主要地層為燕山早期侵入花崗巖，隧道最大埋深約為680.0 m，開挖高度約15.1 m，跨度約16.1 m。

梁山隧道施工至洞身DK96+505~DK96+533 km段揭露1條L7富水破碎帶，經(jīng)探測(cè)發(fā)現(xiàn)[18]：斷層破碎帶寬為21~25 m，傾角為75°~85°，與隧道平面夾角約61°；L7斷層破碎帶夾層物質(zhì)主要為全風(fēng)化花崗巖以及輝綠巖和閃長(zhǎng)玢巖，軟弱夾層主要呈砂土狀，含水量極高。施工記錄顯示[18?19]：該隧道施工至破碎帶處1月內(nèi)，先后發(fā)生了4次嚴(yán)重突水涌泥事故，總涌泥量達(dá)20 000 m3，導(dǎo)致地表塌陷成1個(gè)25 m深的橢圓形陷坑。針對(duì)破碎帶突水災(zāi)害情況，施工方采用“水平旋噴超前加固圈+拱墻大管棚超前支護(hù)+掌子面水平旋噴改良預(yù)加固”的方案進(jìn)行處理。

3.2 耦合模型的程序設(shè)計(jì)

在Visual Studio 2012環(huán)境中，借助C++語(yǔ)言開發(fā)平臺(tái)，將上述理論模型編譯為DLL(動(dòng)態(tài)鏈接庫(kù))，嵌入有限差分軟件FLAC3D中進(jìn)行二次開發(fā)。由于應(yīng)力場(chǎng)的變化影響了滲流場(chǎng)，滲透系數(shù)的改變又會(huì)引起液體、顆粒與骨架的形變，可見水力作用是相互耦合的。對(duì)每加載時(shí)步，其計(jì)算流程如下。

1) 由初始條件求解單元的初始應(yīng)力狀態(tài)及水壓力。

2) 設(shè)置初始滲透率0與初始孔隙率0，計(jì)算孔隙液體、顆粒與骨架的速度等參量。

3) 在計(jì)算過(guò)程中對(duì)滲透系數(shù)參數(shù)進(jìn)行修正，并根據(jù)其修正值k+1，計(jì)算孔隙率的變化值n+1。

3.3 計(jì)算模型和模擬方案

根據(jù)工程地質(zhì)勘察資料，隧道掌子面整體位于斷層破碎帶內(nèi)部，因此，選取隧道與破碎帶相交部位作為計(jì)算區(qū)域建立數(shù)值計(jì)算模型。為了防止軟弱破碎圍巖開挖引起過(guò)大變形，導(dǎo)致計(jì)算過(guò)程不能收斂，建模時(shí)采用變更后的支護(hù)方案，僅從圍巖的滲流場(chǎng)、應(yīng)力場(chǎng)與位移場(chǎng)分布規(guī)律分析突水機(jī)理?？紤]到破碎帶對(duì)圍巖開挖的擾動(dòng)影響較大，模型寬度方向兩側(cè)各取隧道跨度的5倍，高度方向向下取隧道高度的5倍。為了簡(jiǎn)化計(jì)算，模型的上邊界取為初始地下水位(地表以下90 m)，并將水位線以上的巖層作為地面超載施加于上部邊界。

計(jì)算模型的上部邊界固定孔隙水壓力為0，左、右兩邊界和底部邊界均設(shè)為不透水邊界；模型的力學(xué)邊界采用位移邊界條件，底部邊界約束豎向位移，上部邊界為自由邊界，左、右兩端采用水平位移約束。具體計(jì)算模型如圖3所示。

圍巖的物理力學(xué)參數(shù)由工程勘察報(bào)告提供，并根據(jù)參數(shù)折減的要求進(jìn)行相應(yīng)修正。梁山隧道圍巖及襯砌的力學(xué)參數(shù)如表1所示。數(shù)值計(jì)算中圍巖與襯砌的力學(xué)模型均采用Mohr?Coulomb模型。

數(shù)據(jù)單位：m

表1 梁山隧道圍巖及襯砌物理力學(xué)參數(shù)

3.4 滲流場(chǎng)分布特性

為了解隧道開挖前后圍巖滲流場(chǎng)的變化情況，提取隧道拱頂、拱底和邊墻節(jié)點(diǎn)的孔隙水壓力變化時(shí)程曲線，如圖4所示。從圖4可以看出：在隧道開挖過(guò)程中，由于地下水體能量的瞬間釋放，各監(jiān)測(cè)點(diǎn)的孔隙水壓力均大幅度降低；此外，圍巖孔隙水壓力的變化幅度具有一定的規(guī)律，從拱底節(jié)點(diǎn)至拱頂節(jié)點(diǎn)其變化程度逐漸增大。

圖5所示為隧道開挖后圍巖的孔隙水壓力云圖。從圖5可以看出：由于開挖部分孔壓消散，隧道周邊圍巖的孔隙水壓力嚴(yán)重下降，造成滲流場(chǎng)發(fā)生改變；地下水在水頭壓力作用下，形成1個(gè)漏斗形狀的滲水區(qū)域。這與富水破碎帶巖層隧道開挖的實(shí)際情況相符。

根據(jù)模型各單元的流速矢量計(jì)算開挖斷面的涌水量，結(jié)果如圖6所示。從圖6可見：隧道開挖斷面的涌水速度急劇增加，并在短時(shí)間內(nèi)趨于穩(wěn)定，計(jì)算過(guò)程中最大涌水量為16.72 m3/h。這說(shuō)明地下水的滲流通道在施工擾動(dòng)影響下迅速變大，在隧道揭露位置形成較大的水排泄通道，導(dǎo)致地下水在隧道內(nèi)快速排泄。

1—拱頂節(jié)點(diǎn)；2—拱底節(jié)點(diǎn)；3—邊墻節(jié)點(diǎn)。

圖5 隧道開挖后孔隙水壓力云圖

圖6 隧道斷面涌水量變化曲線

3.5 應(yīng)力場(chǎng)分布特性

由于計(jì)算模型的上邊界取為地下水位線，在計(jì)算過(guò)程中，各巖層始終處于完全飽和狀態(tài)(地下水飽和度為1)，因此，在分析應(yīng)力場(chǎng)計(jì)算結(jié)果時(shí)，既可以采用有效應(yīng)力分析，也可以采用總應(yīng)力分析。圖7所示為隧道開挖后圍巖的豎向總應(yīng)力分布云圖。

圖7 隧道開挖后豎向總應(yīng)力云圖

從圖7可以看出：由于隧道開挖，巖體內(nèi)部的力學(xué)平衡狀態(tài)發(fā)生急劇變化，導(dǎo)致圍巖應(yīng)力場(chǎng)重新分布；隧道局部圍巖甚至出現(xiàn)了應(yīng)力集中現(xiàn)象，這在兩側(cè)邊墻的底角部位尤為明顯；隧道開挖對(duì)應(yīng)力場(chǎng)的擾動(dòng)影響隨距離的增加逐漸減弱，遠(yuǎn)離隧道的圍巖，應(yīng)力狀態(tài)基本上仍處于隨深度變化遞增的趨勢(shì)。

3.6 位移場(chǎng)分布特性

圖8所示為隧道開挖后拱頂、拱底和邊墻節(jié)點(diǎn)的豎向位移變化時(shí)程曲線。從圖8可以看出：由于開挖使隧道附近圍巖受到的支撐作用急劇變化，隧道拱頂發(fā)生瞬時(shí)沉降、拱底發(fā)生瞬時(shí)底臌現(xiàn)象，而隧道邊墻節(jié)點(diǎn)的瞬時(shí)變形不明顯；隨著滲流時(shí)間增長(zhǎng)，各監(jiān)測(cè)點(diǎn)相繼發(fā)生沉降，且沉降量越來(lái)越大；不同部位監(jiān)測(cè)點(diǎn)的時(shí)程曲線規(guī)律基本一致，但當(dāng)計(jì)算最終趨于平衡后，各監(jiān)測(cè)點(diǎn)的最大沉降量有所差異，其中拱頂部位的沉降量最大。

1—拱頂節(jié)點(diǎn)；2—拱底節(jié)點(diǎn)；3—邊墻節(jié)點(diǎn)。

提取隧道開挖后圍巖的豎向變形云圖，如圖9所示。由圖9可見：計(jì)算區(qū)域內(nèi)所有位移均為負(fù)值，表明隧道附近圍巖的位移方向豎直向下；隧道兩側(cè)花崗巖的位移變化較平緩，且大致處于隨深度增加變形遞減的趨勢(shì)；斷層破碎帶區(qū)域巖體破碎且受開挖擾動(dòng)影響較大，導(dǎo)致發(fā)生嚴(yán)重的巖土體位移變形；斷層破碎帶處巖體發(fā)生整體滑移變形，在隧道正上方形成橢圓形的塌陷區(qū)域，這恰與梁山隧道斷層上部地表塌陷的實(shí)際破壞規(guī)律相符。

圖9 隧道開挖后豎向變形云圖

4 深埋隧道富水?dāng)鄬油凰畽C(jī)理分析

4.1 深埋隧道富水?dāng)鄬油凰蛩?/p>

深埋隧道穿越富水破碎帶時(shí)在施工過(guò)程中易發(fā)生突水災(zāi)害，其影響因素可以歸結(jié)為2個(gè)方面：斷層破碎帶特殊的水文地質(zhì)條件和隧道開挖與支護(hù)過(guò)程的施工擾動(dòng)。

4.1.1 水文地質(zhì)條件

斷層破碎帶特殊的水文地質(zhì)條件是近斷層隧道施工中發(fā)生重大水害的根本原因。受多期構(gòu)造運(yùn)動(dòng)的作用，斷層地質(zhì)構(gòu)造非常復(fù)雜，巖體節(jié)理裂隙異常發(fā)育，自穩(wěn)能力極弱。此外，斷層破碎帶的含水構(gòu)造與導(dǎo)水構(gòu)造，形成了隧道周圍潛在的突水通道，具有很強(qiáng)的致災(zāi)能力，在很大程度上決定了突水風(fēng)險(xiǎn)。

4.1.2 工程擾動(dòng)

開挖和支護(hù)過(guò)程的施工擾動(dòng)是隧道突水事故發(fā)生的直接原因。以往的工程實(shí)踐證明[5]：在沒有受到施工擾動(dòng)的情況下，原始地質(zhì)條件的滲流場(chǎng)和地應(yīng)力場(chǎng)都處于一定的平衡狀態(tài)。也就是說(shuō)，導(dǎo)致隧道突水事故的直接原因不是構(gòu)造交匯區(qū)的控水優(yōu)勢(shì)，而是工程擾動(dòng)改變了地下水的滲流狀態(tài)，引起構(gòu)造帶地層變形活化。

4.2 深埋隧道富水?dāng)鄬油凰輽C(jī)理

結(jié)合工程實(shí)際與數(shù)值計(jì)算結(jié)果分析，認(rèn)為斷層突水實(shí)質(zhì)上是破碎圍巖的力學(xué)平衡和地下水的滲流平衡因隧道的施工擾動(dòng)而發(fā)生急劇變化，引起應(yīng)力重新分布以及地下水能量的瞬間釋放，并以流體狀態(tài)快速涌向工程臨空面的一種動(dòng)力破壞現(xiàn)象。其突水塌陷機(jī)理可以圖10進(jìn)行概括。

隧道施工揭露斷層后，附近圍巖的滲流場(chǎng)、應(yīng)力場(chǎng)及位移場(chǎng)均不斷發(fā)生變化。在地下水滲透力的拖曳下，突水口附近的巖體顆粒受到的動(dòng)水壓力超過(guò)顆粒間的內(nèi)聚力，巖體顆粒的極限平衡狀態(tài)被打破，隨孔隙空間的流體發(fā)生遷移，形成新的滲流通道。隨著滲流作用時(shí)間延長(zhǎng)，滲流通道迅速變大，使得揭露位置成為地下水的排泄通道。對(duì)于深埋隧道，上部巖土體可以有效阻止地表水對(duì)地下水的補(bǔ)給，因此，地下水的快速排泄將導(dǎo)致地下水位急劇下降，引起水力梯度顯著增加。此時(shí)，基巖面附近的巖土體產(chǎn)生瞬時(shí)負(fù)壓，巖體受到的滲透力也迅速增大，局部軟弱巖體會(huì)隨水流一起涌出，最終引起質(zhì)變，發(fā)生隧道內(nèi)的突水和涌泥事故。隨著突水涌泥的發(fā)展，越來(lái)越多的地下水和巖土體流失，上部破碎巖體受到周圍完整巖體的支撐力越來(lái)越小，在真空吸蝕力的作用下，軟弱破碎帶的巖土體就會(huì)發(fā)生塌陷。

圖10 隧道斷層突水塌陷機(jī)理

5 結(jié)論

1) 進(jìn)行深埋隧道富水軟弱帶地段的開挖分析時(shí)，應(yīng)該考慮裂隙巖體滲流場(chǎng)與應(yīng)力場(chǎng)之間的耦合效應(yīng)，不能把二者分開單獨(dú)考慮。在構(gòu)建飽和破碎巖體的流固耦合效應(yīng)模型時(shí)，對(duì)破碎圍巖的滲流過(guò)程進(jìn)行合理假設(shè)是關(guān)鍵。建立在不同假設(shè)條件下的理論推導(dǎo)，所得到的巖體骨架、液體介質(zhì)與懸浮顆粒的運(yùn)動(dòng)方式以及顆粒流失過(guò)程中液體滲流與巖體骨架變形的耦合方程是不同的。

2) 基于連續(xù)介質(zhì)力學(xué)和變質(zhì)量動(dòng)力學(xué)理論，建立的巖體變質(zhì)量滲流?變形耦合模型符合工程實(shí)際情況，推導(dǎo)出顆粒流失的質(zhì)量守恒方程、介質(zhì)運(yùn)動(dòng)方程和考慮材料性質(zhì)變化的耦合效應(yīng)方程，能夠從理論上解決模型中液體滲流與巖體骨架變形之間的動(dòng)態(tài)平衡關(guān)系。

3) 計(jì)算得到了開挖過(guò)程中圍巖的滲流場(chǎng)、應(yīng)力場(chǎng)與位移場(chǎng)分布特性。隧道開挖后圍巖應(yīng)力場(chǎng)與滲流場(chǎng)的改變使地下水在水頭壓力作用下向工程臨空面滲透，形成1個(gè)漏斗形狀的滲水區(qū)域；同時(shí)，隧道上方的破碎巖體發(fā)生嚴(yán)重的滑移變形，形成1個(gè)橢圓形的塌陷區(qū)域。對(duì)比分析數(shù)值計(jì)算結(jié)果和實(shí)際觀測(cè)結(jié)果基本相符，進(jìn)一步說(shuō)明了滲流?變形耦合模型的合理性和正確性。

4) 斷層突水實(shí)質(zhì)上是破碎圍巖的力學(xué)平衡和地下水的滲流平衡因隧道的施工擾動(dòng)而發(fā)生急劇變化，導(dǎo)致圍巖應(yīng)力重分布以及地下水能量瞬間釋放。隧道施工揭露斷層后，打破了巖土體的極限平衡狀態(tài)，巖體顆粒隨滲流作用迅速流失，形成新的地下水排泄通道，局部軟弱巖體在真空吸蝕力和滲透力的作用下隨水流一起涌出，引發(fā)隧道內(nèi)的突水和涌泥事故。

參考文獻(xiàn)：

[1] JEON S, KIM J, SEO Y, et al. Effect of a fault and weak plane on the stability of a tunnel in rock:a scaled model test and numerical analysis[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(Suppl 1): 658?663.

[2] KUN M, ONARGAN T. Influence of the fault zone in shallow tunneling: A case study of Izmir Metro Tunnel[J]. Tunnelling and Underground Space Technology, 2013, 33(1): 34?45.

[3] 李術(shù)才, 周宗青, 李利平, 等. 巖溶隧道突水風(fēng)險(xiǎn)評(píng)價(jià)理論與方法及工程應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報(bào), 2013, 32(9): 1858?1867. LI Shucai, ZHOU Zongqing, LI Liping, et al. Risk evaluation theory and method of water inrush in karst tunnels and its applications[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(9): 1858?1867.

[4] LI Shucai, ZHOU Zongqing, LI Liping, et al. Risk assessment of water inrush in karst tunnels based on attribute synthetic evaluation system[J]. Tunnelling and Underground Space Technology, 2013, 38(3): 50?58.

[5] 石少帥. 深長(zhǎng)隧道充填型致災(zāi)構(gòu)造滲透失穩(wěn)突涌水機(jī)理與風(fēng)險(xiǎn)控制及工程應(yīng)用[D]. 濟(jì)南: 山東大學(xué)土建與水利學(xué)院, 2014: 1?30. SHI Shaoshuai. Study on seepage failure mechanism and risk control of water inrush induced by filled disaster structure in deep-long tunnel and engineering applications[D]. Jinan: Shandong University. School of Civil Engineering, 2014: 1?30.

[6] 陳國(guó)順. 富水軟弱帶地段深埋隧道受力特征及施工技術(shù)分析[D]. 長(zhǎng)沙: 中南大學(xué)土木工程學(xué)院, 2013: 7?35. CHEN Guoshun. Research on mechanics characteristic and construction technology for deep tunnel in watery and weak stratum[D]. Changsha: Central South University. School of Civil Engineering, 2013: 7?35.

[7] 劉招偉, 何滿潮, 王樹仁. 圓梁山隧道巖溶突水機(jī)理及防治對(duì)策研究[J]. 巖土力學(xué), 2006, 27(2): 228?232. LIU Zhaowei, HE Manchao, WANG Shuren. Study on karst water burst mechanism and prevention countermeasures in Yuanliangshan tunnel[J]. Rock and Soil Mechanics, 2006, 27(2): 228?232.

[8] 趙延林, 曹平, 萬(wàn)文, 等. 巷道前伏承壓溶洞突水災(zāi)變流固耦合分析[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 45(5): 1598?1604. ZHAO Yanlin, CAO Ping, WAN Wen, et al. Fluid-solid coupling analysis of water bursting catastrophe from concealed confined karst cave before roadway[J]. Journal of Central South University (Science and Technology), 2014, 45(5): 1598?1604.

[9] 范威, 王川, 金曉文, 等. 吉蓮高速公路鐘家山隧道涌突水條件分析[J]. 水文地質(zhì)工程地質(zhì), 2015(2): 38?43. FAN Wei, WANG Chuan, JIN Xiaowen, et al. Water inrush condition analysis of the Znongjiashan tunnel in the Jilian Highway[J]. Hydrogeology and Engineering Geology, 2015(2): 38?43.

[10] ALIJA S, TORRIJO F J, QUINTA-FERREIRA M. Geological engineering problems associated with tunnel construction in karst rock masses: the case of Gavarres tunnel (Spain)[J]. Engineering Geology, 2013, 157(4): 103?111.

[11] 姚邦華, 茅獻(xiàn)彪, 魏建平, 等. 考慮顆粒遷移的陷落柱流固耦合動(dòng)力學(xué)模型研究[J]. 中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 2014, 43(1): 30?35. YAO Banghua, MAO Xianbiao, WEI Jianping, et al. Study on coupled fluid-solid model for collapse columns considering the effect of particle transport[J]. Journal of China University of Mining & Technology, 2014, 43(1): 30?35.

[12] 姚邦華. 破碎巖體變質(zhì)量流固耦合動(dòng)力學(xué)理論及應(yīng)用研究[D]. 徐州: 中國(guó)礦業(yè)大學(xué)深部巖土力學(xué)與地下工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 2012: 62?75. YAO Banghua. Research on variable mass fluid-solid coupling dynamic theory of broken rock mass and application[D]. Xuzhou: China University of Mining & Technology. State Key Laboratory for Geomechanics and Deep Underground Engineering, 2012: 62?75.

[13] MARK W D. Recent results of experiments with Saffman-Taylor flow[J]. Nonlinear Evolution and Chaotic Phenomena, 1988, 176: 313?317.

[14] MCCLOUD K V, MAHER J V. Experimental perturbations to Saffman-Taylor flow[J]. Physics Reports, 1995, 260(3): 139?185.

[15] 倪晉仁, 王光謙, 熊育武, 等. 泥石流的結(jié)構(gòu)兩相流模型: Ⅱ.應(yīng)用[J]. 地理學(xué)報(bào), 1998, 53(1): 77?85. NI Jinren, WANG Guangqian, XIONG Yuwu, et al. Conceptual two-phase flow model of debris flow : Ⅱ. Application[J]. Acta Geographica Sinica, 1998, 53(1): 77?85.

[16] 費(fèi)祥俊. 高濃度渾水的黏滯系數(shù)(剛度系數(shù))[J]. 水利學(xué)報(bào), 1982(3): 57?63. FEI Xiangjun. Coefficient of velocity for muddy water with high concentration[J]. Journal of Hydraulic, 1982(3): 57?63.

[17] YOUNGS E G, SMART P, HERBERTSON J G. Drainage design: soil physical principles[M]. New York: Springer Publishing Company, 1992: 25?60.

[18] 周鑫. 廈深鐵路梁山隧道L7深埋富水軟弱帶超前預(yù)加固體系支護(hù)機(jī)理分析[J]. 現(xiàn)代隧道技術(shù), 2013, 50(4): 138?145. ZOU Xin. On the support mechanism of the pre-reinforcement system for the L7 deep, weak water-rich zone passed through by the Liangshan tunnel[J]. Modern Tunnelling Technology, 2013, 50(4): 138?145.

[19] 林福地. 水平旋噴加固在廈深鐵路梁山隧道斷層突水涌泥段處理中的應(yīng)用[J]. 現(xiàn)代隧道技術(shù), 2012, 49(Suppl): 322?329. LIN Fudi. On the application of horizontal jet grouting for water inrush through fault by the Liangshan Tunnel[J]. Modern Tunnelling Technology, 2012, 49(Suppl): 322?329.

(編輯 陳燦華)

Water burst mechanism of deep buried tunnel passing through weak water-rich zone

LI Tingchun1, LYU Lianxun1, 2, DUAN Huiling1, CHEN Wei1, 3

(1. Disaster Prevention and Mitigation of Shandong Provincial Key Laboratory of Civil Engineering,Shandong University of Science and Technology, Qingdao 266590, China;2. BGI Engineering Consultants Ltd., Beijing 100038, China;3. Jinan Urban Construction Group Co. Ltd., Jinan 250000, China)

Considering that there will be security risk of inrush accident when tunnel passes through the weak water-rich zone with complex underground stress field and seepage field environment, a seepage model of porous medium was established through the research of the seepage characteristics of broken rock mass. And theoretical model of fluid-solid coupling of saturated broken rocks was deduced based on the theory of dynamics of variable mass and continuum mechanics. Taking the L7 weak water-rich zone of the Liangshan Tunnel in Zhangzhou, Fujian Province, as engineering background, the distribution features of stress field, seepage field and displacement field of the surrounding rock were analyzed, and inrush collapse mechanism of the fault fracture zone was summarized. The results show that water inrush is actually caused by the stress distribution of surrounding rock and energy release of the groundwater, which results from the dramatic change of mechanical equilibrium of rock mass and seepage balance of the groundwater for construction disturbance. When the fault is exposed, rock particles loss in the pore space to form new seepage channel, leading to welling of groundwater under the pressure of water head to form funnel-shaped water seepage area. With the increase of the seepage time, groundwater and rock mass experience gradual loss and the fractured rock mass on the top of tunnel have serious slip deformation, forming oval subsidence area. This conclusion basically coincides with practical situation. The fluid-solid interaction model proposed has high practical significance to understand the seepage mechanics and prevention design of deep buried tunnel water inrush disaster.

deep buried tunnel; rock mass; fault fracture zone; water inrush mechanism; fluid-solid coupling

10.11817/j.issn.1672-7207.2016.10.025

U457+.5

A

1672?7207(2016)10?3469?08

2015?10?10；

2015?12?21

國(guó)家自然科學(xué)基金資助項(xiàng)目(51279096, 51409154)；山東省自然科學(xué)基金資助項(xiàng)目(ZR2012EEM030)(Projects(51279096, 51409154) supported by the National Natural Science Foundation of China; Project(ZR2012EEM030) supported by the Natural Science Foundation of Shandong Province)

李廷春，博士，教授，博士生導(dǎo)師，從事裂隙巖石力學(xué)特征研究；E-mail：tchli_sd@163.com