王顯

（滄州水利勘測規(guī)劃設(shè)計(jì)院有限公司，河北 滄州 061000）

1 研究背景

水利工程建設(shè)對優(yōu)化水資源配置，充分發(fā)揮有限水資源的價(jià)值具有十分重要的作用和意義?；诖?，我國各級政府不斷加大在水利工程建設(shè)領(lǐng)域的投入力度，建設(shè)了一大批復(fù)雜地質(zhì)環(huán)境條件下的輸水工程。在長距離輸水工程中，輸水隧洞往往是不可或缺的重要組成部分。在長距離輸水隧洞建設(shè)中，經(jīng)常需要穿越河谷地帶，從而面臨淺埋富水砂軟弱地層，給工程建設(shè)的順利進(jìn)行帶來諸多困擾。因此，在穿越上述地層的輸水隧洞施工過程中，需要采用一定的輔助措施對圍巖進(jìn)行預(yù)支護(hù)。水平旋噴技術(shù)是地下洞室工程中常用的預(yù)支護(hù)手段，通常沿地下洞室的洞口周邊以外插的方式分布和咬合，從而形成拱棚預(yù)支結(jié)構(gòu)，該方式十分適合富水砂性地層地下洞室工程的預(yù)支護(hù)，可以有效緩解涌水、涌沙地質(zhì)災(zāi)害[1]。

雖然超前水平旋噴預(yù)支護(hù)技術(shù)在富水砂性地層的地下工程建設(shè)中得到廣泛應(yīng)用，顯示出相對于超前小導(dǎo)管注漿等傳統(tǒng)預(yù)支護(hù)方式的優(yōu)越性[2]，但是，其成樁工程效果受到施工工藝、地質(zhì)環(huán)境等諸多因素的影響，進(jìn)而影響其支護(hù)效果。目前，業(yè)界關(guān)于水平旋噴拱棚的受力特征與變形規(guī)律研究尚不完善，導(dǎo)致工程實(shí)踐應(yīng)用缺乏必要的理論依據(jù)，從而限制了該加固方式的推廣和應(yīng)用[3]?；诖?，此次研究以某輸水隧洞穿越河谷段工程為依托，利用數(shù)值模擬的方式探討輸水隧洞水平旋噴拱棚撓度的變化規(guī)律，以便為相關(guān)研究和工程實(shí)踐提供必要的支持和借鑒。該段輸水隧洞為馬蹄形斷面設(shè)計(jì)，其寬度和高度分別為4.5 m 和4.8 m，水平旋噴拱棚樁的樁徑為0.5 m，相鄰樁中心距為3.5 m，共需要設(shè)置25 根水平旋噴樁。

2 模型計(jì)算

2.1 有限元模型的構(gòu)建

此次研究利用FLAC3D 有限元軟件進(jìn)行數(shù)值計(jì)算模型的構(gòu)建[4]，利用上、下臺階法模擬隧洞的開挖過程[5]。為了有效消除模型邊界效應(yīng)可能對計(jì)算結(jié)果產(chǎn)生的影響，模型的橫向尺寸為隧洞寬度的10 倍，下部深度為拱底向下5 倍的隧洞高度，模型上部模擬至地表，最終確定模型尺寸為45.0 m×40.0 m×50.0 m。

為保證計(jì)算精度，模型的網(wǎng)格劃分系數(shù)取0.5，對隧洞開挖巖土體、水平旋噴拱棚及襯砌結(jié)構(gòu)等關(guān)鍵部位進(jìn)行加密處理[6]?？紤]到對水平旋噴拱棚和周圍巖土體采用共節(jié)點(diǎn)的計(jì)算方式會導(dǎo)致拱棚和土層之間的協(xié)同變形，從而產(chǎn)生較大的計(jì)算誤差，因此，計(jì)算中應(yīng)考慮旋噴樁和巖土體的接觸效應(yīng)，利用界面單元解決上述問題。界面單元可以有效提升有限元網(wǎng)格的柔性，防止非物理應(yīng)力結(jié)果的出現(xiàn)，保證計(jì)算結(jié)果的精確性[7]。此次研究模型劃分為21 246 個(gè)單元，22 357 個(gè)節(jié)點(diǎn)。

2.2 邊界條件和計(jì)算參數(shù)

考慮到河谷內(nèi)地表水和地下水均十分豐富，因此，在模型計(jì)算過程中將各種介質(zhì)均設(shè)置為飽和狀態(tài)，且保持水位高度不變[8]。模型采用軟件中提供的完全流固耦合方法進(jìn)行管棚軸力和撓度的計(jì)算。隧洞圍巖巖土體采用實(shí)體單元模擬，其中，圍巖采用小應(yīng)變硬化本構(gòu)模型，水平旋噴拱棚和初襯采用彈性本構(gòu)模型，并將其設(shè)置為不透水材料。隧道的開挖斷面設(shè)置為透水邊界，水壓力初始值設(shè)定為0。

模型材料參數(shù)是影響計(jì)算結(jié)果的重要因素，研究中以背景工程的地質(zhì)勘查資料和相關(guān)施工規(guī)范為依據(jù)確定模型材料物理力學(xué)參數(shù)，如表1所示。

表1 模型材料物理力學(xué)參數(shù)

2.3 計(jì)算方案

結(jié)合背景工程實(shí)際和相關(guān)研究成果，本文選擇地層滲透系數(shù)、地下水位和隧洞埋深等因素進(jìn)行計(jì)算分析。地層滲透系數(shù)分別設(shè)置為50.000，5.000，0.500，0.050，0.005 m/d；地下水位設(shè)置為0（拱頂部位），1.0，2.0，3.0，4.0 m；隧道埋深設(shè)置為7.0，9.0，11.0，13.0，15.0 m。研究中固定兩個(gè)參數(shù)值不變，以撓度為表征量，對第3 個(gè)參數(shù)不同取值條件下水平旋噴拱棚的力學(xué)特征進(jìn)行模擬計(jì)算。

3 計(jì)算結(jié)果與分析

3.1 地層滲透系數(shù)

研究中保持地下水水位2.0 m、隧道埋深11.0 m不變，利用構(gòu)建的有限元模型對不同地層滲透系數(shù)條件下水平旋噴拱棚的撓度進(jìn)行模擬計(jì)算，從計(jì)算結(jié)果中提取出上臺階開挖進(jìn)尺為10.0 m 時(shí)拱棚拱頂部位的樁體撓度，并繪制出變化曲線，見圖1。由圖1 可知，隨著地層滲透系數(shù)的減小，拱頂部位的撓度最大值呈現(xiàn)出不斷增大的變化特點(diǎn)，且撓度的最大值出現(xiàn)在施工開挖掌子面的后方附近；撓度并非隨著地層滲透系數(shù)的減小呈現(xiàn)線性變化，當(dāng)滲透系數(shù)減小到一定值時(shí)，撓度增長較為有限，究其原因，主要是水平旋噴樁拱棚可以發(fā)揮顯著的預(yù)支護(hù)作用，使拱棚上部的地下水難以形成豎向滲流，只能通過水平方向進(jìn)行滲流，而地層滲透系數(shù)的減小會抑制滲流作用的發(fā)生，導(dǎo)致水平旋噴拱棚的撓度增大。

圖1 不同地層滲透系數(shù)撓度變化曲線

3.2 地下水位

研究中保持地層滲透系數(shù)0.500 m/d、隧道埋深11.0 m 不變，利用構(gòu)建的有限元模型對不同地下水位條件下水平旋噴拱棚的撓度進(jìn)行模擬計(jì)算，從計(jì)算結(jié)果中提取出上臺階開挖進(jìn)尺為10.0 m 時(shí)拱棚拱頂部位的樁體撓度，并繪制出變化曲線，見圖2。由圖2 可以看出，在不同地下水位工況下，水平旋噴拱棚拱頂部位樁體的撓度變化規(guī)律基本一致，均呈現(xiàn)出先增大后減小并逐漸趨于穩(wěn)定的變化趨勢，撓度的最大值出現(xiàn)在施工開挖掌子面后方附近；隨著地下水水位的升高，水平旋噴拱棚拱頂部位樁體的撓度呈現(xiàn)出不斷增大的變化趨勢，且地下水位越高，其對撓度的影響越大，由此可見，地下水水位對拱頂部位水平旋噴樁的撓度影響較大。為保證隧洞施工安全進(jìn)行，對于地下水水位較高的洞段，建議采取增加旋噴樁樁徑、樁間咬合厚度及注漿漿液強(qiáng)度的工程措施，同時(shí)，還應(yīng)該加強(qiáng)施工掌子面的防滲處理。

圖2 不同地下水位撓度變化曲線

3.3 隧洞埋深

研究中保持地層滲透系數(shù)0.500 m/d、地下水位2.0 m 不變，利用構(gòu)建的有限元模型對不同隧洞埋深條件下水平旋噴拱棚的撓度進(jìn)行模擬計(jì)算，從計(jì)算結(jié)果中提取出上臺階開挖進(jìn)尺為10.0 m 時(shí)拱棚拱頂部位的樁體撓度，并繪制出變化曲線，見圖3。由圖3 可以看出，在不同隧洞埋深條件下，水平旋噴拱棚拱頂部位樁體的撓度變化規(guī)律基本一致，均呈現(xiàn)出先增大后減小并逐漸趨于穩(wěn)定的變化趨勢；撓度的最大值出現(xiàn)在施工開挖掌子面后方附近；隨著隧洞埋深的增加，水平旋噴拱棚拱頂部位樁體的撓度呈現(xiàn)出不斷增大的變化趨勢，且增大的趨勢表現(xiàn)出比較顯著的線性特征，究其原因，主要是地下水位不變的情況下，隧洞埋深的增大會導(dǎo)致拱棚上部巖土體壓力和地下水作用的穩(wěn)步增加，旋噴樁撓度也呈現(xiàn)出不斷增大的變化特點(diǎn)。因此，對于埋深較大的洞段也需要采取上文提出的工程措施，保證施工安全。

圖3 不同隧洞埋深撓度變化曲線

4 結(jié)語

此次研究以具體工程為依托，利用數(shù)值模擬的方式探討和分析了輸水隧洞水平旋噴拱棚撓度影響因素及相關(guān)規(guī)律，并提出了具體的工程建議，研究結(jié)論本身對相關(guān)理論研究和工程設(shè)計(jì)具有重要的支持和借鑒作用。當(dāng)然，此次研究在水平旋噴拱棚成樁效果良好的前提下展開，沒有探究樁體缺陷的影響，也沒有考慮水平旋噴拱棚出現(xiàn)漏水情況對自身受力特征的影響。顯然，上述情況在工程施工中往往是不可避免的。因此，在今后的研究中需要針對上述情況進(jìn)行進(jìn)一步的分析研究，以獲得更為全面和準(zhǔn)確的研究成果。