耿　萍，丁　梯，何　悅，晏啟祥

(西南交通大學(xué)交通隧道工程教育部重點實驗室，成都　610031)

?

節(jié)理巖體中隧道開挖后滲流場重分布影響因素研究

耿萍，丁梯，何悅，晏啟祥

(西南交通大學(xué)交通隧道工程教育部重點實驗室，成都610031)

摘要：在富水地區(qū)的節(jié)理巖體中開挖隧道會打破原有的滲流平衡，形成新的滲流場，造成突水等災(zāi)害?？紤]流固耦合作用，運用離散元理論，建立垂直交叉節(jié)理巖體模型，以洞周測點的流量為指標(biāo)，從地應(yīng)力、水位高度和節(jié)理法向剛度等方面進(jìn)行隧道開挖后節(jié)理巖體滲流場重分布的影響因素研究，并從開挖前后節(jié)理隙寬和水壓力變化分析影響產(chǎn)生的原因，結(jié)論可為進(jìn)一步研究圍巖穩(wěn)定性和突水災(zāi)害等提供參考。

關(guān)鍵詞：富水隧道；節(jié)理巖體；滲流；離散元；流固耦合

1概述

隧道在富水地區(qū)施工時常伴有高壓涌水等問題，富水地區(qū)隧道開挖時的滲流問題也是學(xué)者們關(guān)注的重點[1-3]。

在長期的地質(zhì)作用過程中，巖體常伴有不同發(fā)育程度的節(jié)理產(chǎn)生，由于其滲透系數(shù)遠(yuǎn)大于巖體本身，因而成為了地下水的主要流動通道。在節(jié)理巖體中進(jìn)行隧道施工時，洞周節(jié)理處的滲流問題是考慮的重點。此時，一方面滲流產(chǎn)生滲透壓力，減少了節(jié)理面上的有效法向應(yīng)力，從而影響了巖體的變形和穩(wěn)定；另一方面，滲流特性又依賴于節(jié)理的應(yīng)力和變形，由此一來節(jié)理巖體變形和滲流便形成了復(fù)雜的耦合作用[4-7]。

隧道開挖后節(jié)理巖體的變形和滲流特性會受到諸多因素的影響，并進(jìn)一步影響二者的耦合作用。文獻(xiàn)[8-10]對富水地區(qū)隧道開挖后的流固耦合作用進(jìn)行了研究；文獻(xiàn)[11-13]分析了節(jié)理走向、傾角、構(gòu)造應(yīng)力等對圍巖穩(wěn)定性的影響；文獻(xiàn)[14-15]分析了圍巖加固注漿對隧道開挖后含裂隙巖體滲流場重分布的影響；文獻(xiàn)[16]對節(jié)理在不同接觸狀態(tài)下的滲流特性進(jìn)行了模型試驗研究。

本文運用離散元理論，建立垂直交叉節(jié)理巖體模型進(jìn)行數(shù)值分析，從地應(yīng)力、水位高度和節(jié)理法向剛度等因素對隧道開挖后節(jié)理巖體滲流場重分布的影響進(jìn)行研究，可為進(jìn)一步研究圍巖穩(wěn)定性和突水災(zāi)害等提供參考。

2基于離散元的流固耦合算法

離散元理論的基礎(chǔ)是假設(shè)介質(zhì)為離散塊體的集合，因而某一塊體不會受到鄰接塊體變形協(xié)調(diào)的約束，但該塊體自身需滿足平衡方程。液體流動的數(shù)值計算是利用區(qū)域結(jié)構(gòu)來實現(xiàn)的，即對于一個封閉系統(tǒng)來說，存在一個區(qū)域網(wǎng)絡(luò)，每一個區(qū)域假設(shè)充滿具有統(tǒng)一壓力的液體，并且通過接觸與相鄰塊體相連。在圖1中，區(qū)域①～⑤被接觸點A～F分割，區(qū)域⑤是一個空間，區(qū)域①、③和④代表節(jié)理，區(qū)域②處于兩條節(jié)理的交界處，接觸點是塊體相互作用力傳遞的點。

注:①～⑤表示區(qū)域，A-F表示接觸點。圖1　節(jié)理間的流體模型

區(qū)域中液體壓力假設(shè)根據(jù)靜水梯度線性變化，相鄰區(qū)域間的流動由壓力差計算決定。根據(jù)接觸類型的不同，流動速度有兩種計算方法。

注：秒表示計算時間步，不與實際時間等價圖2　各影響因素下洞周各測點的流量值

(1)對于點-點接觸或點-邊接觸，

液體由壓力為p1的區(qū)域流到壓力為p2的區(qū)域的流動速度q由式(1)可得。

(1)

其中，kc為單點接觸滲透系數(shù)；

(2)

式中，ρw為流體密度；g為重力加速度；y2、y1為區(qū)域中心的重力方向坐標(biāo)值。

(2)對于邊-邊接觸，可定義接觸長度，為左邊區(qū)域的一半加上右邊區(qū)域的一半。流動速度q由式(3)可得

(3)

其中，kj為節(jié)理滲透系數(shù)；a為接觸隙寬；l為區(qū)域間接觸所分配的長度。

3模型建立及計算工況

依托某工程實例，建立垂直交叉節(jié)理模型(傾角45°/135°)，節(jié)理間距D=1.5 m，隧道直徑7.4 m。模型上部采用應(yīng)力邊界條件，左、右和下部采用位移邊界條件，同時在四周施加水力邊界條件，巖體及節(jié)理相關(guān)參數(shù)見表1。計算工況按影響因素分3類，共16組，見表2。

表1　巖體相關(guān)參數(shù)

表2　計算工況

4滲流影響因素分析

記錄節(jié)理與隧道相交的20個測點的流量時程，得到各測點穩(wěn)定后的流量值并求和得洞周總流量值，如圖2所示。各工況下的流量極值點都出現(xiàn)在拱腳及拱肩處即隧道橫截面共軛45°方向，說明所研究的3種影響因素均不影響洞周流量大小的分布規(guī)律。

4.1　地應(yīng)力的影響

圖2(a)顯示：隨著豎向地應(yīng)力增加，洞周各測點的流量都逐漸減小。當(dāng)豎向地應(yīng)力為2 MPa時，洞周最大流量和總流量分別為0.095、1.54 m3/s；當(dāng)豎向地應(yīng)力增長到12 MPa時，洞周最大流量和總流量降低到0.020 8、0.26 m3/s。

通過對比各測點的隙寬和水壓力隨地應(yīng)力的變化趨勢，發(fā)現(xiàn)主要有兩種類型，見圖3和圖4。

圖3　隙寬和水壓力隨地應(yīng)力變化趨勢(類型1)

圖4　隙寬和水壓力隨地應(yīng)力變化趨勢(類型2)

圖3和圖4均顯示：初始狀態(tài)，隙寬隨地應(yīng)力增大逐漸減小但最終恒定，水壓力始終保持不變；開挖穩(wěn)定后的水壓力隨地應(yīng)力增大有微弱減小。

隙寬變化的規(guī)律不唯一，這是由于地應(yīng)力會影響節(jié)理法向力，開挖后節(jié)理法向力、水壓力和隙寬又相互影響。當(dāng)?shù)貞?yīng)力較小時，節(jié)理法向力較小，若不足以彌補開挖后水壓力的釋放，則使隙寬減小，反之便會使隙寬增大；同理隨著地應(yīng)力增大，開挖后節(jié)理巖體法向應(yīng)力釋放足以彌補水壓力減小造成的隙寬減小，從而使隙寬增大，且增量隨地應(yīng)力增大而增大。但無論哪種情況，開挖后的隙寬總是隨地應(yīng)力的增大而減小，洞周流量也因此減小。

4.2　水位高度的影響

圖2(b)顯示：隨著水位高度的增加，洞周各測點的流量均都有所增加。當(dāng)水位高度為75 m時，洞周最大流量和總流量分別為0.025 7、0.45 m3/s；當(dāng)水位高度為175 m時，洞周最大流量和總流量升高到0.091 9、1.14 m3/s。

各測點的隙寬和水壓力隨地應(yīng)力的變化趨勢一致，如圖5所示。

圖5　隙寬和水壓力隨水位高度變化趨勢

圖5顯示：初始狀態(tài)時隙寬和水壓力均隨水位高度增長而增大，但開挖穩(wěn)定后，二者隨水位高度增長均僅有微弱增長。由此表明開挖穩(wěn)定后的節(jié)理隙寬，即圍巖變形不受初始水壓力的控制，因為此時水壓力瞬間釋放，已經(jīng)遠(yuǎn)小于圍巖壓力，不足以影響隙寬。開挖穩(wěn)定后隙寬和水壓力基本保持不變，使流量隨水位高度增長而增大。

4.3　節(jié)理法向剛度的影響

圖2(c)顯示：隨著節(jié)理法向剛度增加，洞周各測點的流量及洞周總流量都逐漸增加。當(dāng)節(jié)理法向剛度為7.5 GN/m時，洞周最大流量和總流量分別為0.019 3、0.195 m3/s；當(dāng)節(jié)理法向剛度為22.5 GN/m時，洞周最大流量和總流量升高到0.067 2、1.09 m3/s。

各測點的隙寬和水壓力隨地應(yīng)力的變化趨勢也一致，如圖6所示。

圖6　隙寬和水壓力隨節(jié)理法向剛度變化趨勢

圖6顯示：初始狀態(tài)時，地應(yīng)力不變，節(jié)理隙寬隨法向剛度增大而增大，且增大趨勢逐漸減緩，而初始水壓力則不隨之變化；開挖穩(wěn)定后，不同法向剛度下節(jié)理隙寬均增大，且增量隨法向剛度增大而減小。由于不同法向剛度下的節(jié)理隙寬差別較大，而水壓力釋放量基本保持不變，因此使得洞周測點的流量受法向剛度的影響較大，且可假設(shè)當(dāng)法向剛度達(dá)到無窮大時節(jié)理隙寬最大且將保持不變，此時洞周流量為最大。

5結(jié)論

(1)隧道開挖后，水壓力和節(jié)理巖體圍巖壓力得到釋放，洞周滲流場重新分布，地應(yīng)力、水位高度和節(jié)理法向剛度對于洞周滲流均有影響，但均不影響洞周流量大小的分布規(guī)律，拱腳和拱肩即隧道橫截面共軛45°方向為流量極值點。

(2)地應(yīng)力的增大使開挖前后的節(jié)理隙寬均減小，從而使開挖后洞周流量隨之減小；水位高度影響初始水壓力，不影響開挖后的洞周節(jié)理隙寬和水壓力，因此洞周流量隨水位高度增長而增大；節(jié)理法向剛度直接影響開挖前后的節(jié)理隙寬，法向剛度越大，開挖前后的隙寬越大，因此洞周流量也越大。

(3)節(jié)理巖體中隧道開挖后洞周流量大小主要取決于水壓力釋放量和開挖后節(jié)理隙寬。水位高度影響前者，地應(yīng)力和節(jié)理法向剛度影響后者。本文結(jié)論可為進(jìn)一步研究節(jié)理巖體下隧道施工圍巖穩(wěn)定性和圍巖突水災(zāi)害等問題提供參考。

參考文獻(xiàn)：

[1]楊會軍，王夢恕.深埋長大隧道滲流數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報，2006，25(3):511-519.

[2]王建秀，胡力繩，張金，等.高水壓隧道圍巖滲流-應(yīng)力耦合作用模式研究[J].巖土力學(xué)，2008(S)：237-240.

[3]趙建宇.中天山隧道高壓富水段涌水量探測與連通試驗研究[J].鐵道標(biāo)準(zhǔn)設(shè)計，2013(5):83-86.

[4]吉小明，王宇會.隧道開挖問題的水力耦合計算分析[C]∥第八屆全國青年巖石力學(xué)與工程學(xué)術(shù)大會暨第八屆全國巖石破碎與粉碎學(xué)術(shù)大會論文集.2005:848-852.

[5]肖維民，夏才初，鄧榮貴.巖石節(jié)理應(yīng)力-滲流耦合試驗系統(tǒng)研究進(jìn)展[J].巖石力學(xué)與工程學(xué)報，2014(S2).

[6]陳勝宏，王鴻儒，熊文林.節(jié)理巖體滲流和變形耦合性質(zhì)的研究[J]. 武漢水利電力學(xué)院學(xué)報，1990(1).

[7]高新強，仇文革.隧道襯砌外水壓力計算方法研究現(xiàn)狀與進(jìn)展[J].鐵道標(biāo)準(zhǔn)設(shè)計，2004(4):84-87.

[8]陳衛(wèi)忠，楊建平，楊家?guī)X，等.裂隙巖體應(yīng)力滲流耦合模型在壓力隧道工程中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報，2006，25(12):2384-2391.

[9]李地元，李夕兵，張偉，等.基于流固耦合理論的連拱隧道圍巖穩(wěn)定性分析[J].巖石力學(xué)與工程學(xué)報，2007，26(5):1056-1064.

[10]高健，張義同，喬金麗.滲透力對隧道開挖面穩(wěn)定性影響分析[J].巖土工程學(xué)報，2009(10):1547-1553.

[11]Exploration for rock engineering: proceedings of the Symposium on Exploration for Rock Engineering， Johannesburg， 1-5 November 1976[M]. AA Balkema， 1976.

[12]蔣爵光，李雋蓬，錢惠國.在不同構(gòu)造應(yīng)力作用下節(jié)理巖體隧道的穩(wěn)定性[J].西南交通大學(xué)學(xué)報，1982(1):20-23.

[13]任德惠，張平.不同傾角結(jié)構(gòu)面對巷道穩(wěn)定性的影響[J].煤炭學(xué)報，1988(3):11-17.

[14]高新強，仇文革，孔超.高水壓隧道修建過程中滲流場變化規(guī)律試驗研究[J].中國鐵道科學(xué)，2013(1):50-58.

[15]張成平，張頂立，王夢恕，等.高水壓富水區(qū)隧道限排襯砌注漿圈合理參數(shù)研究[J].巖石力學(xué)與工程學(xué)報，2007，26(11):2270-2276.

[16]夏才初，王偉，曹詩定.節(jié)理在不同接觸狀態(tài)下的滲流特性[J].巖石力學(xué)與工程學(xué)報，2010，29(7):1297-1306.

Research on Seepage Field Redistribution of Tunnel in Fractured Rock Mass

GENG Ping， DING Ti， HE Yue， YAN Qi-xiang

(Key laboratory of Transportation Tunnel Engineering， Ministry of Education， School of Civil Engineering，

Southwest Jiaotong University， Chengdu 610031， China)

Abstract：Tunneling in water-affluent area is likely to break the balance of seepage and form new seepage field， which may result in water outburst. Considering fluid-solid coupled interaction， this paper studies the effect of different crustal stress， water level and joint normal stiffness on seepage of tunnel in fractured rock. Causes of seepage changes resulted from different impact factors are analyzed by analyzing alterations of the width of fissure and water pressure before and after excavation. The conclusions may provide references for further researches on rock stability and water irruption.

Key words：Water-affluent tunnel; Jointed rock mass; Seepage; Discrete element; Fluid-solid coupled interaction

作者簡介：耿萍(1964—)，女，教授，博士，E-mail：gengp@scjm.gov.cn。

基金項目：國家自然科學(xué)基金(51378433，51578457)

收稿日期：2015-06-08； 修回日期：2015-06-17

中圖分類號：U451

文獻(xiàn)標(biāo)識碼：ADOI:10.13238/j.issn.1004-2954.2016.01.024

文章編號：1004-2954(2016)01-0111-03