劉 晨，曾冠銘，王 松，余來斌，楊 濤，田荊濤，梁超恒，盧亞婷

(1.中國建筑第八工程局有限公司，上海 200000； 2.中建八局有限公司南方公司， 廣西 南寧 530022；3.西京學院 土木工程學院，陜西 西安 710123)

突水涌泥是富水軟弱地層(如花崗巖風化槽)修建隧道時較為常見的地質災害之一[1-4]。京九線歧嶺隧道[5]修建過程中就遭遇大范圍的花崗巖風化槽地層，多次出現突水涌泥災害，施工洞口多次被廢棄，由此也被稱為“天字一號工程”。廣西大浦高速公路平田隧道最大涌水量約1 800 m3/h，最大涌泥量達2 500 m3；包茂高速均昌隧道[6]最大涌水量達1 200 m3/h。正在修建的廣西巴馬-憑祥公路大新至憑祥段也遇到近十處巖溶、風化巖等各種不良地質，給施工帶來極大挑戰(zhàn)，嚴重影響工程安全施工與進度。

此類地層極易誘發(fā)突水突泥災害主要源于其巖體顆粒松散，遇水易崩解。特別在高水壓作用下，崩解松散顆粒物極易隨水分遷移流失，進而引起地層孔隙、滲透特性增強，加劇涌水通道發(fā)展。該滲流侵蝕特性與裂隙巖體突水有明顯不同[7-11]?；诖耍瑒⒔鹑冉Y合理論分析、室內試驗和現場分析角度從變質量遷移特性及非線性流態(tài)角度分析其災害機制[12-13]，提供了花崗巖風化槽地層隧道滲流侵蝕涌水涌泥模型，為災害模擬和治理提供了理論依據。

實踐表明，預注漿技術是改進地層水力特性常用方法[14-16]。但注漿后，由于高水力梯度的滲透侵蝕作用，隧道掌子面和未注漿區(qū)域的破碎巖層仍有可能漸進損失破裂。由此也常常需要預留一定的安全厚度防止高水壓擊穿巖體。過去，參考帷幕注漿止?jié){墻設計厚度，工程界常采用3 m的厚度作為預留的安全厚度?？紤]到現場巖石以及注漿巖體裂隙孔隙的存在，與澆筑混凝土有比較大的差異。因此，預留的安全厚度也需要重新探討。

本文依托廣西某小間距隧道，采用滲流侵蝕耦合突水突泥理論模型，通過Comsol多場耦合系統(tǒng)開展隧道安全厚度研究。在獲得隧道孔隙、滲透、涌水涌泥量基礎上，評估合適的安全厚度，為后續(xù)工程提供參考。

1 變質量突水突泥理論模型

為了描述全強風化花崗巖涌水過程出現的非線性滲流和變質量特性，基于溶質運移和滲流力學導出了滲流侵蝕耦合理論模型。模型的主要假設、定義及控制方程如下：

1.1 基本假設

(1) 假定全強風化花崗巖土體由水相(f)，骨架顆粒相(s)、可動顆粒相(fs)三相組成，且可動顆粒相假定隨水流在骨架空隙間移動。

(2) 假定可動顆粒相和水在任意時刻的速度均相同，不考慮水流攜沙過程中的能量損失。

(3) 假定全強風化花崗巖土體孔隙為有效孔隙，且孔隙被水完全填充。

(4) 假定骨架顆粒相為剛性，同時水相不可壓縮。

1.2 參量定義

(1) 各相體積分數定義為：na=dVa/dV, 其中a對應的三相(包括水相(f)，骨架顆粒相(s)、可動顆粒相(fs))。dVa是a相體積，dV是三相的總體積。

(2) 各相局部密度定義為：ρa=dma/dV, 其中dma是a相的質量。各相實密度則定義為：ρa′=dma/dVa。

孔隙率定義為：φ=dVV/dV=nf+nfs，細顆粒濃度定義為：c=dVfs/dVV=nfs/(nf+nfs),其中dVV是孔隙的體積。

1.3 控制方程

(1) 質量平衡方程。根據多孔介質流動動力學理論，多相流體的質量守恒方程為：

(1)

式中:va′和ma′是a相流速和質量速率。上述方程第一項表述a相的質量隨時間變化率，第二項是累計量，右邊項表述質量遷移速率。

(2) 孔隙率演化方程。根據石油開采過程中的出砂侵蝕現象及大壩管涌侵蝕現象，引進如下描述固體顆粒在水流侵蝕作用下的典型孔隙率演化方程[12]，用以描述花崗巖崩解細顆粒在高壓水的遷移流失特征：

(2)

式中:λ是方程系數;φ為孔隙率;φm為孔隙率演化穩(wěn)定值[12]; |q| 是體積通量。

(3) 滲流方程。滲流方程采用最常用的Darcy滲流[13]。

(3)

(4)

式中:k0為介質初始滲透率。

為此，由上述質量方程、孔隙率演化方程及滲流方程組成的耦合滲流侵蝕模型，可用于描述突水突泥的實時演化過程。

2 工程背景

2.1 工程概況

依托的隧道工程位于廣西東南部的構造侵蝕型中低山地貌區(qū)，地形起伏較大，溝谷發(fā)育，隧道埋深最大達到450 m，最低也超過90 m；設計隧型為分離式小凈距，左右兩洞車道中心線間距30 m，凈距17 m，隧道凈空(孔-寬×高)為1-10.75×5 m。隧道左洞起止樁號為DK6+455—DK10+725，長4 270 m，右洞起止樁號為CK6+477—CK10+765，長4 288 m。采用新奧法施工，初期支護為錨桿噴射混凝土支護、二襯采用現澆混凝土支護。穿越地層主要為不同風化程度的加里東期花崗巖，開挖涌水量揭露較大的區(qū)域集中在K7+500—K8+160附近。該區(qū)域隧道埋深約100 m，巖體主要為強度極低的全強風化花崗。在此地層施工時面臨多次突水突泥地質災害，涌水量最大達到1 280 m3/h，最大涌泥量超過8 000 m3。

2.2 水文地質條件

根據地質勘探和地球物理勘探，對隧道所處地質環(huán)境進行了評估，如圖1所示。隧道沿線主要穿越不同風化等級的花崗巖地層。在滲涌水較大區(qū)域(K7+620—K8+200)，其巖性主要是中風化、強風化和全風化花崗巖。對于強風化花崗巖，其節(jié)理非常發(fā)育，巖體質量很差。巖層的滲透性和儲水能力都很高。而全風化花崗巖經過劇烈風化后，形成以大量的石英石顆粒和黏土礦物成分特征。其巖石質量很差，巖體結構在擾動下極易崩解。同時，由于及其發(fā)育的節(jié)理、孔隙存在，此類地層儲水能力及導水能力較強?，F場鉆孔及物探揭示此區(qū)域處于強富水區(qū)域，單孔涌水量揭露達到400 m3/h，水壓也高達0.7 MPa～1.5 MPa，如圖2所示。根據室內試驗和現場壓水試驗，得到關鍵的物理力學參數如強度及滲透特性，如表1所示。

圖1 隧道地質剖面圖

表1 風化花崗巖物理力學參數

圖2 鉆孔涌水量

當隧道揭露此類地層時，受高水壓水作用，全強風化花崗巖顆粒極易受水流侵蝕，從而形成滲漏通道。此外，由于高風險地區(qū)的所處的“V”型谷地特征，水資源補給充足且迅速，一旦發(fā)生容易造成大規(guī)模的地下水流入和地面塌陷。

3 隧道臨界防突厚度數值研究

對于此類富水軟弱地層，現場采用帷幕注漿進行預加固。帷幕注漿加固圈厚度為5 m。注漿完成后，根據檢查孔、P-Q-t曲線等方法檢測達到開挖質量后，開始開挖工作。安全厚度是既保證隧道安全開挖，又是達到較好的經濟效益(縮短工期)的關鍵指標。下文采用仿真手段對掌子面前方所需的安全厚度進行研究。

3.1 模型建立

根據上述理論模型，采用Comsol多場耦合軟件進行方程求解，并開展全強風化花崗巖隧道開挖臨界防突厚度研究，為實際工程提供參考。根據實際工程條件，建立如圖3簡化模型。模型地層從上往下簡化成全強風化層、中風化。并在掌子面前方拱頂注漿圈與未注漿圈交界處布置監(jiān)測點，用以觀察臨界防突厚度不同時交界處孔隙率、涌水量等變化規(guī)律，從而反應隧道突水演化特性。

圖3 臨界防突厚度計算模型

模型的邊界條件為：滲流邊界為頂部位置為地下水位處(y=0)，水壓力沿深度線性分布，揭露的掌子面為大氣壓力?？紫堵始翱蓜宇w粒濃度邊界為：模型四周為狄氏邊界，隧道出口為開放邊界，供顆粒自由遷移。

根據室內及現場試驗，得到仿真模型的初始孔隙率、滲透率等參數，見表2。其中，隧道全強風化花崗巖原巖的滲透率、孔隙率分別為Kr、φr，注漿后滲透率為Kg、φg。中風化花崗巖的滲透率、孔隙率分別為Kz、φz。注漿后孔隙率按照Kozeny-Carmen方程進行反算。通過模擬開挖過程，研究防突厚度僅剩8 m～2 m的滲流侵蝕發(fā)展過程。計算時間為12 h。

表2 臨界防突厚度計算模型參數

3.2 模型結果與分析

3.2.1 孔隙率變化規(guī)律

圖4為隧道掘進后剩余不同防突厚度h的孔隙率演化云圖。從圖4可以看出，當防突厚度h>6 m時，經歷長時間的滲流侵蝕作用，圍巖孔隙率發(fā)展很小。h=5 m時，掌子面前方有部分顆粒開始遷移，有微通道逐步發(fā)育。h=4 m時，孔隙率顯著發(fā)展，并在掌子面前方逐步形成滲流通道。防突厚度進一步減少至4 m以下時，掌子面前方土體幾乎全部流失，已形成顯著的突水通道，通道貫穿至未注漿區(qū)域，表明極可能誘發(fā)突水突泥災害。

圖4 不同臨界防突厚度下孔隙率演化云圖

孔隙率增長源于地層顆粒遷移。同時，而孔隙率變化速率亦與滲流速度直接相關。防突厚度較大時，水力梯度較小，結合達西定律可知滲流速度也較小。由此引起的孔隙率變化則相對更小。相反，防突厚度較低時，高水力梯度誘發(fā)的高滲流速度加快孔隙率變化。

3.2.2 滲透率變化規(guī)律

圖5顯示滲透性隨時間呈現三階段變化，即初始的緩慢增長，隨后的快速增長及最后的穩(wěn)定階段。初始時，此時巖土體顆粒流失較少，引起的孔隙變化較小，致使巖土體滲透特性變化較為緩慢。隨著顆粒逐漸遷移，引起孔隙率顯著增長(增幅可達到100倍)，根據Kozeny-Carmen方程，可知滲透率將大幅度增長。此時滲流也極可能進入非線性滲流，表明突水突泥災害風險大幅增大。進入第三階段后，由于大量的可動顆粒在第二階段已經流失，此時不再有顯著的顆粒遷移，滲透特性也不再顯著增強。

從圖5還可得到，防突厚度大于4 m時，滲透率增長幅度較??；而低于4 m時，突水演化開始進入快速增長階段，對應的滲透率增長也極為顯著，表明發(fā)生災害風險大幅增加。具體來說，防突厚度在8 m、6 m、5 m、4 m、3 m、2 m時的滲透率最終分別為7.4×10-14/m2、1.5×10-13/m2、2.0×10-13/m2、2.3×10-13/m2、9.9×10-12/m2、1.3×10-10/m2，較地層初始滲透率分別增長了1.5倍、3倍、4倍、4.6倍、198倍、810倍。

圖5 監(jiān)測點滲透率時變曲線

3.2.3 涌水量變化規(guī)律

圖6為掌子面每延米涌水量時變演化曲線?？梢钥闯?，防突厚度大于4 m時，最終涌水量均低于1.5 m3/(m·h)，但降至3 m時，涌水量呈指數型增長，經歷8 h即開始迅速增長，最終達到14.4 m3/(m·h)。繼續(xù)降至2 m，涌水量發(fā)展速率繼續(xù)加快，6 h后即進入快速增長階段，最終涌水量近20 m3/(m·h)。12 h后，防突厚度為8 m、6 m、5 m、4 m、3 m、2 m的隧道每延米涌水量分別為0.2 m3/(m·h)、0.3 m3/(m·h)、0.3 m3/(m·h)、1.5 m3/(m·h)、14.4 m3/(m·h)、17.8 m3/(m·h)。顯然，防突厚度從4 m降至3 m時，涌水量將大幅增大，直接導致突水災害的發(fā)生。而防突厚度較大(>5 m)，對涌水量的降低效果不明顯。因此，從涌水量角度分析，預留4 m～5 m防突厚度較為可行。

圖6 不同防突厚度下掌子面每延米涌水量時變曲線

3.2.4 涌泥量變化規(guī)律

圖7為不同防突厚度下的每延米涌泥量曲線，可以看出，涌泥量及其發(fā)展速率隨防突厚度降低顯著增大。具體來說，防突厚度較大時(>5 m)，經歷10 h后涌泥量方有所增加，而防突厚度較小時(3 m～4 m)，涌泥量經歷6 h后開始增長，而防突厚度僅為2 m時，經歷2 h后隧道涌泥量即顯著增長。

圖7 不同防突厚度下掌子面每延米涌泥量時變曲線

12 h后，防突厚度為8 m、6 m、5 m、4 m、3 m、2 m的每延米涌泥量分別為2.4 m3/m、3.6 m3/m、6.7 m3/m、11.4 m3/m、19.2 m3/m、69.8 m3/m?？梢钥闯觯磕嗔侩S防突厚度減小一直增大，并在3 m降至2 m時涌泥量增長最為顯著，涌泥量增加了近4倍；其次為4 m降至3 m，涌泥量增長近2倍。而防突厚度較大時(>5 m)，涌泥量小得多，不足10 m3/m。從涌泥量角度同樣得到，防突厚度在4 m～5 m對于涌泥量的控制是較為經濟有效的。

通過上述關鍵參數時變規(guī)律分析，可以得到安全厚度在低于4 m時，突水突泥演化極易進入快速演化階段，引起孔隙、滲透等關鍵參數大幅增長，誘發(fā)突水突泥災害。大于5 m時，孔隙率、滲透率等關鍵參數能得到較為顯著控制，不再顯著發(fā)展，沒能形成涌水通道，滲流也大多處在第一階段。從經濟安全角度考慮，對此類巖體隧道掌子面前面開挖，建議預留4 m～5 m的安全厚度。

4 結 論

基于溶質運移理論及多孔介質滲流，給出了可考慮質量遷移和流態(tài)特征的滲流侵蝕耦合理論模型。并基于此理論模型，以均昌隧道為背景，開展隧道開挖的最小安全厚度研究。主要結論如下：

(1) 計算結果表明：在防突厚度在小于4 m時，短時間內掌子面前方顆粒即發(fā)生顯著流失，極易形成突水通道，涌水涌泥量大大超過開挖設計的標準。

(2) 當防突厚度大于4 m時，尤其是大于5 m時，地層圍巖孔隙率、滲流通道得到有效遏制，掌子面的涌水涌泥均得到有效控制。

(3) 從安全、經濟、技術角度分析，防突厚度取4 m～5 m是較為合適的。