高新強(qiáng)，艾旭峰，孔　超

(1.石家莊鐵道大學(xué) 土木工程學(xué)院，河北 石家莊　050043；2.石家莊鐵道大學(xué) 道路與鐵道工程安全保障省部共建教育部重點(diǎn)實(shí)驗(yàn)室，河北 石家莊　050043；3.西南交通大學(xué) 土木工程學(xué)院， 四川 成都　610031；4.西南交通大學(xué) 交通隧道工程教育部重點(diǎn)實(shí)驗(yàn)室，四川 成都　610031)

在高水壓富水區(qū)修建隧道是隧道工程界的技術(shù)難題。大量的工程實(shí)踐證明：在高水壓富水區(qū)，圍巖的穩(wěn)定性、地下涌(突)水和襯砌水壓力是關(guān)系到隧道能否修建成功的主要因素。同時(shí)，在當(dāng)今環(huán)保要求較高的條件下，還必須重視隧道修建對地下水環(huán)境的影響[1-2]。

許多學(xué)者通過試驗(yàn)[3-4]或數(shù)值模擬[5-12]的方法，對隧道滲流問題進(jìn)行了研究。目前，結(jié)合裂隙巖體隧道工程的開挖、注漿堵水、結(jié)構(gòu)防排水等特點(diǎn)對裂隙巖體隧道滲流場所進(jìn)行的數(shù)值分析多采用有限元連續(xù)介質(zhì)方法，而裂隙巖體隧道滲流模型中應(yīng)包含離散介質(zhì)和連續(xù)介質(zhì)，因此僅采用連續(xù)介質(zhì)進(jìn)行數(shù)值模擬存在諸多問題。針對這些問題，本文采用離散元法，基于三維離散元軟件3DEC，僅考慮滲流場，忽略應(yīng)力場與滲流場的耦合作用，建立裂隙巖體隧道滲流分析數(shù)值模型，分析隧道主要施工過程(開挖、噴混凝土、襯砌)中隧道滲流場的變化規(guī)律和特征，并分析圍巖裂隙寬度、注漿圈厚度、注漿圈滲透系數(shù)及噴混凝土滲透系數(shù)等主要滲透參數(shù)對隧道滲流場的影響規(guī)律，為裂隙巖體隧道水荷載的確定提供理論支持。

1　離散元分析模型

1.1　模型尺寸

計(jì)算模型在x，y，z軸上的長度分別為200，190，4 m。隧道中心距模型的頂面100 m，距模型的底面90 m，距模型的左、右側(cè)面均為100 m，如圖1所示。計(jì)算模型中建有垂直于隧道軸向(z軸方向)的貫通裂隙，圖1(b)中模型側(cè)面不同顏色的交界處即為裂隙，假定裂隙為光滑平行，裂隙寬1 mm。

1.2　水壓力邊界條件

計(jì)算模型頂面邊界外水壓力為0，前后端面為不透水邊界，其余邊界按靜水壓力邊界施加。水壓力邊界條件如圖2所示。

圖1　滲流數(shù)值模型

圖2　水壓力邊界條件(單位：MPa)

2　隧道施工過程中滲流場的變化規(guī)律

2.1　初始滲流場

隧道未開挖情況下的初始滲流場如圖3所示。由圖3(a)可知：初始滲流場的水壓力僅與高度位置相關(guān)，即處于靜水壓力狀態(tài)。由圖3(b)可知：滲流速度為0，表明在外邊界不透水的情況下初始滲流場為靜水場。

圖3　初始滲流場

2.2　毛洞滲流場

未注漿情況下開挖隧道，模型外部水壓力邊界條件不變，內(nèi)部在開挖輪廓面處設(shè)置1個(gè)固定孔隙水壓力為0的滲流邊界條件，滲流平衡后毛洞的滲流場如圖4所示。由圖4可知：隧道開挖后，裂隙中的地下水向隧道中心流動，在開挖暴露面上流速最大，在隧道周圍形成降水漏斗，導(dǎo)致地下水大量流失(16 600 m3·(m·d)-1)，隧道上部地下水水位明顯下降。

圖4　毛洞滲流場

計(jì)算過程中對隧道各部位徑向圍巖中水壓力的分布進(jìn)行監(jiān)測，監(jiān)測線位置如圖5所示。

圖5　水壓力監(jiān)測線位置

滲流平衡后，裂隙中各監(jiān)測線上不同部位的水壓力如圖6所示，圖中橫坐標(biāo)表示監(jiān)測線上的點(diǎn)到開挖輪廓面的徑向距離。由圖6可知：隧道開挖后對其上部和隧道周邊0～30 m范圍內(nèi)的滲流場影響最為顯著；到開挖輪廓面的徑向距離越遠(yuǎn)開挖對水壓力的影響越小。

2.3　施作5 m注漿圈并開挖后隧道滲流場

注漿后，注漿區(qū)域裂隙的滲透性將減小，因此模型中通過減小注漿區(qū)域裂隙的寬度來模擬注漿效果。設(shè)注漿區(qū)域裂隙寬度為0.05 mm，為未注漿區(qū)域裂隙寬度的1/20，注漿圈區(qū)域裂隙的滲透系數(shù)為未注漿區(qū)域裂隙滲透系數(shù)的1/400。注漿圈厚度為5 m，內(nèi)外部水壓力邊界條件不變。

滲流平衡后隧道滲流場如圖7所示。由圖7(a)可知：裂隙中的水壓力受影響的范圍主要集中在注漿圈區(qū)域，注漿圈以外區(qū)域基本恢復(fù)靜水壓力狀態(tài)。由圖7(b)可知：水壓力在開挖輪廓面處為0，然后沿徑向遞增，在注漿圈范圍以外逐漸趨近靜水壓力，說明注漿圈起到了堵水的目的。由圖7(c)可知：裂隙中的地下水在注漿圈區(qū)域內(nèi)較注漿圈以外區(qū)域流速大，注漿圈區(qū)域流速在0.02～0.06 m·s-1范圍內(nèi)，此時(shí)，由裂隙滲入隧道的水量約為9.5 m3·(m·d)-1，約為注漿前滲水量(16 600 m3·(m·d)-1)的1/1 660。

圖6　毛洞滲流場各監(jiān)測線上各點(diǎn)的水壓力

圖7　施作5 m注漿圈并開挖后隧道滲流場

2.4　施作噴混凝土后隧道滲流場

3DEC 4.10版本中，流體僅能在裂隙中滲透，塊體是不透水的。因此，在模擬噴混凝土等連續(xù)介質(zhì)的滲透性時(shí)，需要用裂隙介質(zhì)的滲透性來等效。模型中噴混凝土塊體的滲透性采用對應(yīng)裂隙的滲透性來等效，等效過程中主要考慮流體徑向的滲流，用多條徑向裂隙來模擬，如圖8所示。

圖8　噴混凝土層數(shù)值模型

噴混凝土的厚度為0.28 m，滲透系數(shù)為1×10-8m·s-1；模型外部滲流邊界條件不變；在噴混凝土內(nèi)表面設(shè)置1個(gè)固定孔隙水壓力為0的滲流邊界條件，模擬施作噴混凝土后的滲流狀態(tài)。

施作噴混凝土后隧道滲流場如圖9所示。由圖9(a)可知：注漿圈與噴混凝土接觸面上水壓力顯著增加，增加至0.40～0.45 MPa，而注漿圈和圍巖分界面與施作噴混凝土前相比未發(fā)生顯著變化。由圖9(b)可知：施作噴混凝土后，對滲流場的主要影響范圍在注漿圈區(qū)域和噴混凝土區(qū)域，在噴混凝土區(qū)域(0～0.28 m)內(nèi)，水壓力變化迅速，其變化范圍約為0.45 MPa；注漿圈區(qū)域(0.28～5 m)內(nèi)，水壓力變化值相比噴混凝土施作前有所減小，由施作前的約1 MPa減小為施作后的約0.55 MPa。由圖9(c)可知：施作噴混凝土后，由于噴混凝土滲透系數(shù)小，起到了限流作用，注漿圈區(qū)域內(nèi)裂隙中地下水的流速有所減小，流速在0.015～0.04 m·s-1范圍內(nèi)，流動方向基本未發(fā)生變化。

施作噴混凝土后由裂隙滲入隧道的總水量約為5.6 m3·(m·d)-1，約為噴混凝土施作前滲水量(9.5 m3·(m·d)-1)的3/5。

圖9　施作噴混凝土后隧道滲流場

2.5　施作二次襯砌后隧道滲流場

模型中在左、右墻腳位置各設(shè)置φ10 cm縱向盲管，縱向盲管與垂直裂隙貫通。模型外部滲流邊界條件不變，內(nèi)部在縱向盲管出口位置設(shè)置固定孔隙水壓力為0的排水邊界條件，模擬全排的情況。

滲流平衡后隧道滲流場如圖10所示。由圖10(a)可知：施作二次襯砌后，噴混凝土區(qū)域內(nèi)，水壓力的變化范圍與施作二次襯砌前相比有所增大，約為0.6 MPa，注漿圈區(qū)域水壓力變化范圍相應(yīng)減小，約為0.4 MPa；由圖10(b)可知：對滲流場的主要影響范圍為注漿圈和噴混凝土區(qū)域；由圖10(c)可知：施作二次襯砌后，裂隙中的地下水經(jīng)過圍巖、注漿圈、噴混凝土，在二次襯砌背后向縱向盲管處匯集、排出。此時(shí)，由裂隙滲入隧道的總水量約為3.5 m3·(m·d)-1，約為施作二次襯砌前滲水量(5.6 m3·(m·d)-1)的7/10。

圖10　施作二次襯砌后隧道滲流場

3　滲透參數(shù)對隧道滲流場的影響

圍巖裂隙寬度、注漿圈厚度、注漿圈滲透系數(shù)及噴混凝土滲透系數(shù)等滲透參數(shù)對隧道施工過程中隧道滲流場會產(chǎn)生影響。

模型采用施作噴混凝土后的滲流離散元數(shù)值模型，其各參數(shù)的初始值見表1。

3.1　巖體裂隙寬度對隧道滲流場的影響

圍巖裂隙寬度范圍取0.2～1.8 mm，每次變化0.2 mm，其他參數(shù)不變，不同圍巖裂隙寬度時(shí)的水壓力和滲水量如圖11和圖12所示，圖中“1點(diǎn)”指監(jiān)測線1與噴混凝土和圍巖接觸面的交點(diǎn)處，其余類同。由圖11可知：隨著裂隙寬度的增大，噴混凝土與注漿后圍巖接觸面處的水壓力也逐漸增大，水壓力變化范圍約為0.2 MPa。由圖12可知：滲入隧道的水量隨著圍巖裂隙寬度的增加而增大，其變化趨勢與水壓力變化趨勢基本一致，滲水量變化值約為0.3 m3·(m·d)-1。

表1　模型參數(shù)的初始值

圖11　水壓力與圍巖裂隙寬度的關(guān)系

圖12　隧道滲水量與圍巖裂隙寬度的關(guān)系

3.2　注漿圈厚度對隧道滲流場的影響

注漿圈厚度的取值范圍為0～20 m，每次增加2 m，其他參數(shù)保持不變。水壓力與注漿圈厚度的關(guān)系如圖13所示。由圖13可知：隨著注漿圈厚度的增加，水壓力逐漸減??；對裂隙巖體進(jìn)行注漿后，注漿圈的厚度在0～8 m之間變化時(shí)對噴混凝土背后的水壓力的影響最為顯著，水壓力變化值約為0.7 MPa；隨著注漿圈厚度的進(jìn)一步增大，影響減弱?？紤]到施工水平和經(jīng)濟(jì)效益，堵水時(shí)建議注漿圈厚度采用3～8 m。

圖13　水壓力與注漿圈厚度的關(guān)系

隧道滲水量與注漿圈厚度的關(guān)系如圖14所示。由圖14可知：隨著注漿圈厚度的增加，隧道滲水量逐漸減小，其變化趨勢與水壓力的變化趨勢基本一致；注漿圈厚度在0～8 m之間變化時(shí)對隧道滲水量影響較為顯著，滲水量變化值約為8 m3·(m·d)-1；超過8 m以后，影響減弱，滲水量變化值僅為2 m3·(m·d)-1左右。

圖14　隧道滲水量與注漿圈厚度的關(guān)系

3.3　注漿圈滲透系數(shù)對隧道滲流場的影響

分析中注漿圈滲透系數(shù)采用等效滲透系數(shù)，逐漸減小圍巖裂隙寬度，以實(shí)現(xiàn)注漿圈滲透系數(shù)由大到小的變化。

水壓力與注漿圈等效滲透系數(shù)的關(guān)系如圖15所示。由圖15可知：各監(jiān)測點(diǎn)的水壓力隨著注漿圈等效滲透系數(shù)的增加而增加，注漿圈等效滲透系數(shù)小于2.0×10-7m·s-1時(shí)對噴混凝土背后水壓力影響顯著，注漿圈內(nèi)水壓力變化值約為0.9 MPa，當(dāng)注漿圈等效滲透系數(shù)大于2.0×10-7m·s-1之后，水壓力變化較為平緩，變化值約為0.1 MPa；當(dāng)注漿圈滲透系數(shù)大于1.0×10-7m·s-1時(shí)，噴混凝土背后的水壓力較大。因此采用注漿堵水時(shí)，建議注漿圈等效滲透系數(shù)小于1.0×10-7m·s-1。

圖15　水壓力與注漿圈等效滲透系數(shù)的關(guān)系

隧道滲水量與注漿圈等效滲透系數(shù)的關(guān)系如圖16所示。由圖16可知：隧道滲水量隨著注漿圈等效滲透系數(shù)的增大而增加，其變化趨勢與噴混凝土背后水壓力變化趨勢基本一致；注漿圈等效滲透系數(shù)小于2.0×10-7m·s-1時(shí)，對隧道滲水量影響顯著，滲水量的變化值約為11 m3·(m·d)-1；注漿圈等效滲透系數(shù)大于2.0×10-7m·s-1后，影響大大減弱，在2.0×10-7～1.0×10-6m·s-1范圍內(nèi)，滲水量的變化值不足2 m3·(m·d)-1。由此可知，注漿圈等效滲透系數(shù)小于2.0×10-7m·s-1時(shí)堵水效果最為顯著。

圖16　隧道滲水量與注漿圈等效滲透系數(shù)的關(guān)系

3.4　噴混凝土滲透系數(shù)對隧道滲流場的影響

噴混凝土厚度多在20～30 cm之間，其變化范圍較小，對滲流場的影響不顯著，故僅分析噴混凝土滲透系數(shù)對隧道滲流場的影響。

水壓力與噴混凝土滲透系數(shù)的關(guān)系如圖17所示。由圖17可知：噴混凝土滲透系數(shù)在小于4×10-8m·s-1時(shí)，對其背后水壓力影響顯著，水壓力變化值最大可達(dá)0.9 MPa；當(dāng)噴混凝土滲透系數(shù)大于4×10-8m·s-1(此時(shí)，噴混凝土滲透系數(shù)與注漿圈滲透系數(shù)的比值約為1.54，即兩者滲透系數(shù)基本相當(dāng))后，對水壓力的影響逐漸減弱，水壓力變化值很小，僅為0.1 MPa左右。

圖17　水壓力與噴混凝土滲透系數(shù)的關(guān)系

隧道滲水量與噴混凝土滲透系數(shù)的關(guān)系如圖18所示。由圖18可知：在圍巖裂隙寬度為1 mm，注漿圈等效滲透系數(shù)為2.6×10-8m·s-1的情況下，噴混凝土的滲透系數(shù)在0～4.0×10-8m·s-1范圍內(nèi)變化時(shí)，對隧道滲水量的影響最為顯著，其變化值約8 m3·(m·d)-1；當(dāng)噴混凝土滲透系數(shù)大于4.0×10-8m·s-1后，影響大大減弱，滲水量變化值不足2 m3·(m·d)-1。隨著噴混凝土滲透系數(shù)繼續(xù)增大，隧道滲水量趨于此種情況下的上限值，約為10 m3·(m·d)-1。

圖18　隧道滲水量與噴混凝土滲透系數(shù)的關(guān)系

分析中，裂隙巖體、注漿圈、噴混凝土三者滲透性的初始比是恒定的，當(dāng)三者的滲透系數(shù)之間初始比變化時(shí)，噴混凝土背后水壓力或隧道滲水量的變化規(guī)律基本一致，影響顯著的范圍有所不同，即曲線的陡緩不同。文獻(xiàn)[11—14]做過類似的分析計(jì)算。

4　結(jié)　論

(1)高水壓富水區(qū)裂隙巖體隧道施工過程中滲流場特征：①隧道開挖前，初始滲流場為靜水場；②在未注漿情況下開挖隧道，對隧道滲流場影響顯著，地下水大量流失(16 600 m3·(m·d)-1)，隧道上部地下水的水位面明顯下降；③注漿后開挖隧道，對注漿圈(厚5 m)范圍內(nèi)的滲流場影響顯著，對注漿圈外裂隙圍巖滲流場影響不顯著，隧道滲水量約為9.5 m3·(m·d)-1，約為注漿前的1/1 660；④施作噴混凝土后，對噴混凝土和注漿圈范圍內(nèi)滲流場影響顯著，對注漿圈外裂隙圍巖滲流場影響不顯著，隧道滲水量降低為5.6 m3·(m·d)-1，約為施作前的3/5；⑤施作二次襯砌后，噴混凝土內(nèi)水壓力變化范圍有所增大，注漿圈區(qū)域內(nèi)水壓力變化范圍相應(yīng)減小，隧道滲水量降低為3.8 m3·(m·d)-1，約為施作前的7/10。

(2)圍巖裂隙寬度、注漿圈厚度、注漿圈等效滲透系數(shù)及噴混凝土滲透系數(shù)等對裂隙巖體隧道滲流場的影響規(guī)律：①圍巖裂隙寬度在0.2～1.8 mm之間時(shí)，噴混凝土背后的水壓力及隧道滲流量均隨著裂隙寬度的增大而增大，但變化范圍較??；②隨著注漿圈厚度的增加，噴混凝土背后水壓力和隧道滲水量均減小，注漿圈厚度在0～8 m之間時(shí)，對噴混凝土背后水壓力和隧道滲水量影響均顯著，而超過8 m后，影響大為減弱；③隨著注漿圈等效滲透系數(shù)的增加，噴混凝土背后水壓力和隧道滲水量均增加，注漿圈等效滲透系數(shù)小于2.0×10-7m·s-1時(shí)，對噴混凝土背后水壓力和隧道滲水量影響顯著，而大于2.0×10-7m·s-1后，影響大為減弱；④噴混凝土背后的水壓力隨著噴混凝土滲透系數(shù)的增加而減小，而隧道滲水量隨著噴混凝土滲透系數(shù)的增加而增大，噴混凝土滲透系數(shù)小于4.0×10-8m·s-1時(shí)，對噴混凝土背后水壓力和隧道滲水量影響均顯著，而大于4.0×10-8m·s-1后，影響不再明顯。

[1]高新強(qiáng). 高水壓山嶺隧道襯砌水壓力分布規(guī)律研究[D]. 成都：西南交通大學(xué)，2005.

(GAO Xinqiang. Study on the Distribution Rule of Water Press Upon Lining in Subject to High Hydraulic Pressure Mountain Tunnel[D]. Chengdu：Southwest Jiaotong University, 2005. in Chinese)

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

(GAO Xinqiang, QIU Wenge. Current Situation and Development of Computation for Water Pressure outside Tunnel Lining[J]. Railway Standard Design, 2004(12):84-87. in Chinese)

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

(GAO Xinqiang, QIU Wenge, KONG Chao. Test Study on the Variation Law of Seepage Field during the Construction Process of High Water Pressure Tunnel[J]. China Railway Science, 2013，34(1):50-58. in Chinese)

[4]蔣中明，馮樹榮，傅勝，等. 某水工隧洞裂隙巖體高水頭作用下的滲透性試驗(yàn)研究[J]. 巖土力學(xué)，2010，31(3):673-676.

(JIANG Zhongming，F(xiàn)ENG Shurong，F(xiàn)U Sheng， et al. Test Study of Osmotic Behavior of Fractured Rock Mass of Water Tunnel under High Water Pressure[J]. Rock and Soil Mechanics, 2010，31(3):673-676. in Chinese)

[5]孫玉杰，鄔愛清，張宜虎，等. 基于離散單元法的裂隙巖體滲流與應(yīng)力耦合作用機(jī)制研究[J]. 長江科學(xué)院院報(bào)，2009，26(10):62-66，70.

(SUN Yujie, WU Aiqing, ZHANG Yihu，et al. Study on Coupling Mechanism of Stress and Fluid Flow in Fractured Rock Mass Based on Discrete Element Method[J]. Journal of Yangtze River Scientific Research Institute, 2009，26(10):62-66，70. in Chinese)

[6]房倩，張頂立，黃明琦，等. 基于連續(xù)介質(zhì)模型的海底隧道滲流問題分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2007，26(增2):3776-3784.

(FANG Qian, ZHANG Dingli, HUANG Mingqi, et al. Analysis of Seepage Problem Induced by Subsea Tunnel Excavation Based on Continuum Medium Model[J]. Chinese Journal of Rock Mechanics and Engineering, 2007，26(Supplement 2):3776-3784. in Chinese)

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

(CHEN Weizhong, YANG Jianping, YANG Jialing，et al. Hydromechanical Coupled Model of Jointed Rock Mass and Its Application to Pressure Tunnels[J]. Chinese Journal of Rock Mechanics and Engineering, 2006，25(12):2384-2391.in Chinese)

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

(YANG Huijun, WANG Mengshu. Numerical Simulation of Seepage in Deep-Buried Long and Big Tunnels [J]. Chinese Journal of Rock Mechanics and Engineering, 2006，25(3):511-519. in Chinese)

[9]吉小明. 隧道開挖的雙重介質(zhì)水力耦合模型有限元分析[J]. 巖土力學(xué)，2006，27(增刊):305-310.

(JI Xiaoming. Dual-Porosity Media Hydromechanical Analyses of Tunnel Excavation [J]. Rock and Soil Mechanics, 2006，27(Supplement):305-310. in Chinese)

[10]耿萍，于本昌，何悅，等. 隧洞開挖后滲流量影響因素分析 [J]. 鐵道建筑，2014(4):54-57.

(GENG Ping, YU Benchang,HE Yue， et al. Influence Factors Analysis of Seepage Flow after Tunnel Excavation [J]. Railway Engineering, 2014(4):54-57. in Chinese)

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

(ZHANG Chengping, ZHANG Dingli, WANG Mengshu, et al. Study on Appropriate Parameters of Grouting Circle for Tunnels with Limiting Discharge Lining in High Water Pressure and Water-Enriched Region[J]. Chinese Journal of Rock Mechanics and Engineering, 2007，26(11):2270-2276. in Chinese)

[12]李鵬飛，張頂立，趙勇，等. 海底隧道復(fù)合襯砌水壓力分布規(guī)律及合理注漿加固圈參數(shù)研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2012，31(2):280-288.

(LI Pengfei, ZHANG Dingli, ZHAO Yong，et al. Study of Distribution Law of Water Pressure Acting on Composite Lining and Reasonable Parameters of Grouting Circle for Subsea Tunnel [J]. Chinese Journal of Rock Mechanics and Engineering, 2012，31(2):280-288. in Chinese)