宋玉婷， 巢萬(wàn)里

(湖南省交通科學(xué)研究院， 湖南 長(zhǎng)沙　410015)

?

連拱隧道滲流防護(hù)措施的數(shù)值模擬研究

宋玉婷， 巢萬(wàn)里

(湖南省交通科學(xué)研究院， 湖南 長(zhǎng)沙410015)

[摘要]孔隙水滲流作用對(duì)隧道圍巖穩(wěn)定性影響顯著，尤其當(dāng)隧址所處圍巖孔隙水滲流作用較為強(qiáng)烈時(shí)，影響更甚。因此，必要的排水措施如泄水洞的布設(shè)就具有十分重要的意義。以某隧道為工程背景，借助有限差分軟件FLAC3D對(duì)地層由多期洪積扇堆積的某隧道進(jìn)行初始地應(yīng)力狀態(tài)、無(wú)泄水洞布設(shè)以及泄水洞布設(shè)3種工況進(jìn)行數(shù)值模擬研究。通過(guò)對(duì)圍巖滲流狀態(tài)及泄水洞布設(shè)前后的圍巖應(yīng)力及穩(wěn)定狀態(tài)的數(shù)值模擬結(jié)果進(jìn)行的比較分析，確定泄水洞的功效以驗(yàn)證泄水洞的必要性以及泄水洞位置的合理性。通過(guò)對(duì)滲流、應(yīng)力場(chǎng)、位移場(chǎng)及塑性區(qū)的對(duì)比分析，認(rèn)為泄水洞的布設(shè)能有效降低洞周孔隙水壓力、圍巖應(yīng)力和洞周位移變形量，以及泄水洞后巖土體位移、應(yīng)力及工后沉降，泄水洞的合理布設(shè)將顯著改善降水效果。

[關(guān)鍵詞]連拱隧道； 滲流； 數(shù)值模擬； 泄水洞

0引言

孔隙水滲流作用是影響圍巖自穩(wěn)能力和隧道穩(wěn)定性的重要因素[1-3]。掌子面開(kāi)挖引起圍巖孔隙水滲流的重分布，進(jìn)而影響圍巖應(yīng)力場(chǎng)分布，而應(yīng)力場(chǎng)的重分布又將對(duì)滲流場(chǎng)的分布造成影響，并最終形成多場(chǎng)的耦合作用，對(duì)圍巖和隧道支護(hù)結(jié)構(gòu)造成不良影響[4，5]。在孔隙水充裕的富水圍巖地區(qū)，隧道穩(wěn)定性受其影響更甚，防排水及圍巖加固措施的不當(dāng)與不及時(shí)均將將使圍巖應(yīng)力增大，穩(wěn)定性降低，導(dǎo)致隧道突水、突泥甚至掌子面坍塌的災(zāi)害發(fā)生[6，7]。

隧道滲流場(chǎng)的研究手段主要包括理論分析、模型試驗(yàn)和現(xiàn)場(chǎng)監(jiān)測(cè)，其中模型試驗(yàn)?zāi)茌^方便地模擬現(xiàn)場(chǎng)，對(duì)工程進(jìn)行較好的實(shí)際還原，因而在巖土及隧道工程中被普遍認(rèn)可并被廣泛使用。然而，由于滲流場(chǎng)模擬等方面困難，目前并沒(méi)有一定數(shù)量可供參考的文獻(xiàn)。以跨瓊州海峽海底隧道為背景，譚忠盛等對(duì)隧道襯砌在海水滲流環(huán)境下的防排水方式及不同防排水方式下的滲流場(chǎng)分布規(guī)律進(jìn)行了模型試驗(yàn)研究，但不足的是試驗(yàn)不夠全面，具有一定局限性，無(wú)法反應(yīng)隧道襯砌的全受力狀態(tài)[8，9]。

黃甫明等對(duì)不同襯砌注漿方式及量級(jí)條件下的圍巖滲流場(chǎng)分布情況進(jìn)行了研究，得到了孔隙水壓力與滲流場(chǎng)隨襯砌布設(shè)及注漿量級(jí)的相互關(guān)系及規(guī)律[10-12]。

李術(shù)才等以膠州灣海底隧道為背景，采用自行研發(fā)的設(shè)備和傳感器，進(jìn)行了水下隧道的模型試驗(yàn)研究[13，14]。上述研究局限在襯砌外側(cè)或者注漿圈內(nèi)外范圍，將圍巖、注漿圈、襯砌作為一個(gè)整體考慮的水壓分布規(guī)律研究，目前未見(jiàn)報(bào)導(dǎo)。

泄水洞作為一種隧道排水措施，常用于富水隧道施工中，其在排出圍巖孔隙水的同時(shí)將改善隧道滲流場(chǎng)和應(yīng)力場(chǎng)的分布，有助于圍巖和襯砌穩(wěn)定性[15-17]。深入了解掌子面開(kāi)挖過(guò)程中泄水洞的布設(shè)數(shù)量和方式對(duì)隧道滲流場(chǎng)和應(yīng)力場(chǎng)的影響規(guī)律，對(duì)于富水及水下隧道的設(shè)計(jì)和施工具有重要意義。

為研究泄水洞的布設(shè)對(duì)隧道滲流作用的防護(hù)及對(duì)圍巖穩(wěn)定性提高的效果，以某隧道為工程背景，運(yùn)用FLAC3D有限差分法數(shù)值模擬軟件，對(duì)隧道布設(shè)與不布設(shè)泄水洞兩種工況進(jìn)行分析，比較兩種工況下圍巖滲流、應(yīng)力及沉降的不同，分析泄水洞布設(shè)的有效性。

1工程概況

某山嶺鐵路隧道起訖里程為CK615+065～CK633+470，全長(zhǎng)18.405 km，屬特長(zhǎng)隧道。隧道場(chǎng)地地質(zhì)條件復(fù)雜，地層由多期洪積扇堆積而成，與粘土層形成多層承壓水。隧址范圍內(nèi)的洪積扇主要賦存兩層承壓水，下層承壓水水位高于上層承壓水水位，隧道穿過(guò)下層含水層，水位距離為隧道底板為120～160 m；隧道穿過(guò)上層含水層，水位距離隧道底板為50～80 m。

根據(jù)現(xiàn)場(chǎng)地表出露第三系地層形成的邊坡，在無(wú)水狀態(tài)或滲水狀態(tài)下整體可以自穩(wěn)，長(zhǎng)期暴露局部出現(xiàn)剝落或坍塌。在隧道開(kāi)挖過(guò)程中，在巖土結(jié)構(gòu)未破壞且在無(wú)水狀態(tài)下，掌子面應(yīng)可以整體自穩(wěn)，局部膠結(jié)較差的可能出現(xiàn)剝落或坍塌；由于隧道通過(guò)該段地下水位較高，在掌子面在滲水條件是否自穩(wěn)，首先考慮施工擾動(dòng)是否破壞了巖土體結(jié)構(gòu)，同時(shí)由于排水可能導(dǎo)致土體內(nèi)部細(xì)顆粒成分的流失，造成巖土體結(jié)構(gòu)破壞；其次評(píng)估地下水的流動(dòng)產(chǎn)生滲透力對(duì)掌子面穩(wěn)定性的影響。有必要在隧道開(kāi)挖支護(hù)前首先對(duì)高承壓水進(jìn)行排水泄壓，為此對(duì)泄水洞排水泄壓效果進(jìn)行探討。主要研究?jī)?nèi)容：恰當(dāng)選擇本構(gòu)模型與各項(xiàng)參數(shù)，對(duì)泄水洞設(shè)置與否三種工況(即場(chǎng)地初始狀態(tài)、有隧道無(wú)泄水洞、有隧道有泄水洞)進(jìn)行數(shù)值模擬，并比較分析兩種不同工況下的滲流、應(yīng)力及沉降情況。確定泄水洞的功效以驗(yàn)證泄水洞的必要性以及泄水洞位置的合理性。

2模型建立及參數(shù)選擇

2.1FLAC3D在滲流計(jì)算中的數(shù)值表述

FLAC3D提供了巖土體的流體運(yùn)動(dòng)法則和無(wú)滲流場(chǎng)的流動(dòng)零模型，可以對(duì)巖土體中孔隙水的瞬時(shí)流動(dòng)作用進(jìn)行模擬。不同的巖土體材料可以有不同的流動(dòng)模型和法則，模擬過(guò)程可以獨(dú)立進(jìn)行亦可以根據(jù)巖土體本構(gòu)關(guān)系對(duì)孔隙水滲流效應(yīng)的約束結(jié)合進(jìn)行。允許將任何力學(xué)與熱學(xué)模型進(jìn)行耦合計(jì)算的同時(shí)允許熱膨脹的發(fā)生?；谝陨嫌悬c(diǎn)，本次數(shù)值模擬研究采用FLAC3D有限差分法軟件進(jìn)行，將取得較好的計(jì)算結(jié)果。

2.2工況定義及模型建立

根據(jù)隧道地質(zhì)勘探資料，CK628+800～K631+800段為粘質(zhì)黃土、粗圓礫土與粉質(zhì)粘土，基巖為灰?guī)r，圍巖級(jí)別V級(jí)。分別穿越上下層高承壓水，可能存在發(fā)生突水風(fēng)險(xiǎn)，為最不利斷面。隧道開(kāi)挖方式采用傳統(tǒng)的礦山法施工，其開(kāi)挖方法為臺(tái)階法。工況見(jiàn)表1。

表1 工況定義Table1 Casedefinition分類(lèi)工況狀態(tài)工況一不設(shè)置泄水洞工況二設(shè)置泄水洞

原隧道設(shè)計(jì)典型斷面如圖1所示，泄水洞位置在仿真計(jì)算中暫擬為該中部，為達(dá)到較高精度，模型

圖1　隧道典型斷面Figure 1　Typical section of the tunnel

尺寸寬×高為255×271 m，劃分為9 146個(gè)網(wǎng)格單元，建立FLAC3D數(shù)值模型如圖2所示。

圖2　隧道計(jì)算模型Figure 2　Numerical model for multiple-arch tunnel

2.3模型參數(shù)選擇

參考隧道有關(guān)地質(zhì)物探資料及相關(guān)現(xiàn)行規(guī)范的基礎(chǔ)上，取得各材料物理力學(xué)計(jì)算參數(shù)如表2～表4所示。其中，巖土體材料皆為干密度，注漿加固A為工況掌子面注漿孔數(shù)28個(gè)時(shí)加固體物理力學(xué)參數(shù)，注漿加固B為工況掌子面注漿孔數(shù)96個(gè)時(shí)。

表2 巖土層的物理力學(xué)參數(shù)Table2 Physical-mechanicalparametersofrocklayers材料密度/(kg·m-3)彈性模量E/GPa泊松比v上層含水層23404.50.26下層含水層24805.00.25粘聚力/kPa內(nèi)摩擦角/(°)滲透系數(shù)/(m·s-1)45355.90×10-650364.86×10-6

表3 邊墻砂漿錨桿Cable結(jié)構(gòu)單元力學(xué)參數(shù)Table3 Cablestructuralunitmechanicsparametersofmortarbolt密度/(kg·m-3)彈性模量/Pa水泥漿粘結(jié)力/(N·m-1)水泥漿摩擦角/(°)水泥漿剛度/Pa水泥漿外圈周長(zhǎng)/m截面積/m2抗壓強(qiáng)度/N抗拉強(qiáng)度/N78504.50×10101.15×107353.70×1080.770.03141.00×1062.50×106

表4 鋼拱架和網(wǎng)噴初支shell結(jié)構(gòu)單元力學(xué)參數(shù)表Table4 Shellstructuralunitmechanicsparametersofsteelarchandinitialreinforcement密度/(kg·m-3)彈性模量/Pa泊松比厚度/m250010.50.250.2

2.4數(shù)值模擬實(shí)現(xiàn)過(guò)程

為對(duì)比分析泄水洞在滲流隧道中的作用，分別進(jìn)行布設(shè)與不布設(shè)泄水洞的數(shù)值模擬，步驟如下：

① 不設(shè)置泄水洞。

初始平衡狀態(tài)→開(kāi)挖隧道斷面巖芯→開(kāi)挖面噴射混凝土進(jìn)行初期支護(hù)(含臨時(shí)仰拱)。

② 設(shè)置泄水洞。

初始平衡狀態(tài)→開(kāi)挖泄水洞→泄水洞開(kāi)挖面支護(hù)→開(kāi)挖隧道斷面→隧道開(kāi)挖面噴射混凝土進(jìn)行初期支護(hù)(含臨時(shí)仰拱)。

數(shù)值模擬計(jì)算中僅考慮初支，二襯主要作為強(qiáng)度安全儲(chǔ)備在計(jì)算中不予反映。

2.5邊界條件及初始應(yīng)力狀態(tài)

本模型中模型邊界條件分為位移邊界條件、滲流邊界條件和透水邊界條件。根據(jù)工程概況和模型定義，邊界條件分別設(shè)置如下：

位移邊界條件：模型上底面為自由面，下底面約束縱向沉降，左右兩個(gè)方向約束X方向位移、前后兩個(gè)方向約束Y方向位移。

滲流邊界條件：FLAC3D提供模型的不透水邊界，即模型各邊界是不透水的，邊界空壓自由變化。

透水邊界條件：透水條件和飽和度由邊界空壓決定，且空壓的絕對(duì)值必須大于流體抗拉強(qiáng)度。

按上述原則建立某連拱隧道滲流數(shù)值模擬模型，畜平衡后初始應(yīng)力條件和孔隙水壓力場(chǎng)見(jiàn)圖3。

(a) 應(yīng)力場(chǎng)(b) 孔隙水壓力場(chǎng)

3數(shù)值計(jì)算結(jié)果分析

3.1滲流分析

本模型單個(gè)臺(tái)階開(kāi)挖與支護(hù)所花時(shí)間為10 h，定義一個(gè)工作循環(huán)為60 h，即滲流場(chǎng)為一個(gè)工作循環(huán)即60 h后的計(jì)算結(jié)果，如圖4所示。圖4中(a)可以看出開(kāi)挖后空隙水壓力場(chǎng)變化明顯，圍繞開(kāi)挖部分形成了漏斗狀得低孔隙水壓力區(qū)域，開(kāi)挖面附近水壓力等勢(shì)面密集。在中透水的泥灰?guī)r中，開(kāi)挖后僅60 h內(nèi)滲流場(chǎng)即產(chǎn)生了較大變化，可見(jiàn)如不采取止水措施，將會(huì)引起隧道突水涌水，開(kāi)挖施工將無(wú)法進(jìn)行。比較兩種工況條件下，隧道開(kāi)挖后，洞周均一致出現(xiàn)漏斗狀低孔隙水壓力區(qū)域，工況一與工況二的洞周最大孔隙水壓力分別為1.6 MPa，1.5 MPa，工況二的洞周最大孔隙水壓力小于工況一，且設(shè)置泄水洞時(shí)隧道的開(kāi)挖面附近水壓力等勢(shì)面較未設(shè)置泄水洞有較大程度的改善，可見(jiàn)泄水洞降壓效果初步顯現(xiàn)。

(a) 工況一斷面圖(未設(shè)置泄水洞)(b) 工況一大樣圖(未設(shè)置泄水洞)

(c) 工況二斷面圖(設(shè)置泄水洞)(d) 工況二大樣圖(設(shè)置泄水洞)

圖4各工況中泄水洞降壓效果已顯現(xiàn)，但具體分析滲流趨勢(shì)和降壓效果還需考慮圍巖滲流矢量場(chǎng)。圖5各工況下隧道圍巖滲流矢量圖中顯示，兩種工況隧道圍巖滲流矢量場(chǎng)分布基本一致，具體滲流速度值有所不同。工況二(設(shè)置泄水洞)的滲流速度大小相對(duì)未設(shè)置泄水洞情況弱微減少。

(a) 工況一滲流矢量圖(未設(shè)置泄水洞)

(b) 工況一滲流云圖(未設(shè)置泄水洞)

(c) 工況二滲流矢量圖(設(shè)置泄水洞)

(d) 工況二滲流云圖 (設(shè)置泄水洞)

圖6為各工況下隧道斷面涌水量歷史記錄，2種工況隧道斷面涌水量變化趨勢(shì)基本一致，最大涌水量值有所不同。工況一的每延米斷面最大涌水量為2.576×10-3m3/s，工況二為2.214×10-3m3/s，工況二(設(shè)置泄水洞)的最大涌水量相對(duì)未設(shè)置泄水洞情況減少近14%。由此可見(jiàn)，泄水洞的降壓效果較明顯。

(a) 涌水量歷史記錄(未設(shè)置泄水洞)(b) 涌水量歷史記錄(設(shè)置泄水洞)

綜合上述分析可知：隧道穿越承壓水層時(shí)預(yù)先設(shè)置泄水洞后，降壓效果較明顯。

3.2應(yīng)力場(chǎng)分析

圖7可以看出，隧道開(kāi)挖后兩種工況的洞周壓應(yīng)力分布大致相同，均分布在隧道側(cè)壁，拱頂上方出現(xiàn)低應(yīng)力卸荷泡區(qū)。但工況一中側(cè)壁最大壓應(yīng)力為15.6 MPa，工況二為18.5 MPa?？梢?jiàn)設(shè)置泄水洞存在二次圍巖擾動(dòng)效應(yīng)，使得削弱圍巖壓應(yīng)力集中方面的效果不明顯。

(a) 工況一(未設(shè)置泄水洞)(b) 工況二(設(shè)置泄水洞)

圖8中開(kāi)挖后2種工況的洞周拉應(yīng)力分布規(guī)律基本一致，均集中分布在拱肩和墻角，并且與洞周壓應(yīng)力場(chǎng)一樣，在拱頂上方出現(xiàn)低應(yīng)力卸荷泡區(qū)。但工況二較工況一應(yīng)力場(chǎng)分布更趨于均勻，且工況二的最小主應(yīng)力最大值較工況一弱減小?？梢?jiàn)泄水洞的降壓效果對(duì)削弱圍巖拉應(yīng)力集中有一定作用，但效果不明顯。

(a) 工況一(未設(shè)置泄水洞)(b) 工況二(設(shè)置泄水洞)

綜合圖7與圖8的分析結(jié)果：設(shè)置泄水洞能一定程度改善圍巖應(yīng)力分布狀態(tài)，主拉應(yīng)力的拱頂與拱腳部位和主壓應(yīng)力的側(cè)壁是高壓水地層隧道施工應(yīng)重點(diǎn)關(guān)注的部位，必須加強(qiáng)隧道關(guān)鍵部位的監(jiān)控，掌握圍巖變化情況，保證施工安全。

3.3位移場(chǎng)分析

隧道掘進(jìn)時(shí)，對(duì)圍巖的開(kāi)挖擾動(dòng)，致使洞周?chē)鷰r松動(dòng)而出現(xiàn)洞周位移，其中又以拱頂沉降和仰拱隆起的豎向位移，為工程施工中所重點(diǎn)關(guān)注。圖9中各種工況條件下，皆出現(xiàn)較大的豎向位移，具體數(shù)值可歸納如表5所示。

(a) 工況一(未設(shè)置泄水洞)(b) 工況二(設(shè)置泄水洞)

表5 各工況條件下斷面豎向位移值Table5 Settlementsofthe2differentconditions位置沉降數(shù)值/mm誤差/%工況一工況二工況二拱頂55.4121.19-62仰拱22.3111.47-49

根據(jù)表5中2種工況的比較情況，可以看出：泄水洞降壓能有效減少隧道開(kāi)挖斷面的變形量。

圖10中各種工況條件側(cè)壁相對(duì)位移具體數(shù)值同樣可歸納如表6所示。

(a) 工況一(未設(shè)置泄水洞)(b) 工況二(設(shè)置泄水洞)

Figure 10Relative horizontal displacements of the 2 different conditions

從表6可以看出：同樣驗(yàn)證了泄水洞降壓能顯著減少洞周位移變形量。

綜上分析，設(shè)置泄水洞與否的位移場(chǎng)變化有明顯差別。未設(shè)置泄水洞時(shí)洞周產(chǎn)生位移變形較大，設(shè)置泄水洞能顯著減少洞周位移變形量。

表6 各工況條件下斷面水平相對(duì)位移(洞周收斂)Table6 Relativehorizontaldisplacementsofthe2differentconditions位置沉降數(shù)值/mm誤差/%工況一工況二工況二左側(cè)壁95.3320.27-79右側(cè)壁95.5819.50-80

3.4塑性區(qū)場(chǎng)分析

圖11各工況塑性區(qū)分布圖顯示，工況一與工況二的洞周大范圍圍巖均出現(xiàn)過(guò)歷史剪切屈服，拱頂、拱肩、拱腳與仰拱還出現(xiàn)過(guò)歷史張拉屈服，且隧道間的巖柱出現(xiàn)剪切屈服破壞。工況二應(yīng)力屈服狀態(tài)的圍巖范圍相對(duì)有所擴(kuò)展。可見(jiàn)，設(shè)置泄水洞后對(duì)圍巖的二次擾動(dòng)影響圍巖的穩(wěn)定性。2種工況條件下不采取有效的堵水加固措施，隧道洞周?chē)鷰r容易失穩(wěn)。

(a) 工況一(未設(shè)置泄水洞)(b) 工況二(設(shè)置泄水洞)

4結(jié)論

運(yùn)用FLAC3D對(duì)隧道穿越高承壓水地段的開(kāi)挖施工過(guò)程進(jìn)行數(shù)值模擬，針對(duì)未設(shè)置泄水洞、設(shè)置泄水洞兩種工況條件下，對(duì)比分析各工況的開(kāi)挖擾動(dòng)力學(xué)響應(yīng)與流固耦合效應(yīng)，得到如下結(jié)論：

① 圍巖滲流場(chǎng)分布狀態(tài)復(fù)雜，具有間斷、各向異性的特點(diǎn)，尤其是在隧道開(kāi)挖圍巖原平衡被打破之后，滲流場(chǎng)需要經(jīng)歷重分布并最終平衡的過(guò)程。計(jì)算中認(rèn)為未經(jīng)開(kāi)挖的圍巖孔隙水壓力呈水平層狀分布，并沿深度增加。

② 隧道開(kāi)挖后，區(qū)內(nèi)地下水滲流狀態(tài)發(fā)生改變，地下水向隧道內(nèi)排泄。滲流場(chǎng)數(shù)值模擬結(jié)果顯示，隨著松動(dòng)圈巖體滲透系數(shù)的減小，滲入隧道的水量減小，區(qū)內(nèi)地下水位降幅和孔隙水壓力降幅變緩。

③ 掌子面開(kāi)挖后，原水壓力平衡被打破，圍巖滲流狀態(tài)被改變，孔隙水流向隧道洞室，數(shù)值模擬結(jié)果顯示，滲流狀態(tài)與入洞水量的大小與滲流系數(shù)有關(guān)，即襯砌支護(hù)與防排水措施對(duì)滲流破壞的防治具有重要意義。

④ 隧道開(kāi)挖引起洞周較大范圍圍巖屈服破壞，洞周?chē)鷰r發(fā)生掌子面坍塌和涌水突泥等風(fēng)險(xiǎn)的可能性較大，需進(jìn)行堵水加固處理。

⑤ 設(shè)置泄水洞后巖土體位移、應(yīng)力及工后沉降結(jié)果較明顯改善，說(shuō)明將泄水洞設(shè)置位置是合適的，降水效果明顯；設(shè)置泄水洞能有效降低洞周孔隙水壓力、圍巖應(yīng)力和洞周位移變形量，結(jié)合鋼拱架與圍巖注漿加固處理，可以滿足工程要求。

[參考文獻(xiàn)]

[1]王夢(mèng)恕.隧道與地下工程技術(shù)及其發(fā)展[M].北京：北京交通大學(xué)出版社，2004.

[2]Louis C.Rock Hydraulics in Rock Mechanics [M].New York，Wiley，1974.

[3]劉招偉.圓梁山隧道巖溶突水機(jī)理及其防治對(duì)策[D].北京：中國(guó)地質(zhì)大學(xué)，2004.

[4]王遇國(guó).巖溶隧道突水災(zāi)害與防治研究[D].北京：中國(guó)鐵道科學(xué)研究院，2010.

[5]孫少銳，吳繼敏，魏繼紅，等.大跨度雙連拱隧道開(kāi)挖滲流場(chǎng)特征研究[J].勘察科學(xué)技術(shù)，2005(6)：14-19.

[6]楊會(huì)軍，王夢(mèng)恕.深埋長(zhǎng)大隧道滲流數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報(bào)，2006，25(3)：512-519.

[7]Tsang C. R. Coupled behavior of rock joints. Proc International Symposium on Rock Joints.1990(6)：505-518.

[8]石新棟，皇甫明，譚忠盛，等.跨瓊州海峽隧道方案的探討[J].隧道建設(shè)，2010(6)：625-628.

[9]譚忠盛，洪開(kāi)榮，萬(wàn)姜林，等.軟硬不均地層復(fù)合盾構(gòu)的研究及掘進(jìn)技術(shù)[J].巖石力學(xué)與工程學(xué)報(bào)，2006，25(z2)：3945-3952.

[10]皇甫明.暗挖海底隧道圍巖穩(wěn)定性及支護(hù)結(jié)構(gòu)受力特征的研究[D].北京：北京交通大學(xué)，2005.

[11]皇甫明，譚忠盛，王夢(mèng)恕，等.暗挖海底隧道滲流量的解析解及其應(yīng)用[J].中國(guó)工程科學(xué)，2009，07(7)：66-70.

[12]石新棟，皇甫明，譚忠盛，等.跨瓊州海峽隧道方案的探討[J].隧道建設(shè)，2010(6)：625-628.

[13]李術(shù)才，宋曙光，李利平，等.海底隧道流固耦合模型試驗(yàn)系統(tǒng)的研制及應(yīng)用 [J].巖石力學(xué)與工程學(xué)報(bào)，2013(05)：883-890.

[14]李術(shù)才，王凱，李利平，等.海底隧道新型可拓展突水模型試驗(yàn)系統(tǒng)的研制及應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2014，33(12)：2409-2418.

[15]Pinder G F，Gray W G.Finite element simulation in surface and subsurface hydrology[M].New York：Academic Press，1977.

[16]Heuer R E.Estimating rock tunnel water inflow[J].Proceedings-Rapid Excavation and Tunnel conference，1995，7：41-60.

[17]Gerber R E，Boyce J I，Howard K W F.Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard[J].Hydrogeology Journal，2001，9(1)：60-78.

Numerical Simulation on Seepage Safeguard Procedures for Multiple-arch Tunnel

SONG Yuting， CHAO Wanli

(Hunan Communications Research Intitute， Changsha， Hunan 410015， China)

[Abstract]Pore water seepage effect has a significant influence on tunnel rock stability，especially in heavy pore water conditions.It is of great importance for measures such as drainage holes in such conditions.Based on Zhong Tiaoshan multiple-arch tunnel and numeral simulation software FLAC3Dof finite difference method，this paper conducts a numeral simulation research of different working conditions both with and without the drainage holes on tunnel stability.Necessity and effects of drainage holes are certified by rock stress and stability status analysis.According to analysis of seepage field，stress field and displacement field and plastic zones，result shows that drainage hole can reduce pore water pressure，rock stresses and settlements significantly，namely precipitation effect can be improved by drainage holes.

[Key words]multiple-arch tunnel； seepage； numerical simulation； drainage holes

[收稿日期]2016-04-15

[基金項(xiàng)目]湖南省科學(xué)技術(shù)廳科技計(jì)劃項(xiàng)目(2013SK5080)；湖南省交通科技項(xiàng)目(201430)

[作者簡(jiǎn)介]宋玉婷(1983-)，女，湖南長(zhǎng)沙人，工程師，主要從事路橋?qū)I(yè)設(shè)計(jì)研究。

[中圖分類(lèi)號(hào)]U 453.6+1

[文獻(xiàn)標(biāo)識(shí)碼]A

[文章編號(hào)]1674-0610(2016)03-0111-06