王 輝，王 鵬，梁明純

(1.北京科技大學(xué) 土木工程系，北京 100083；2.北京科技大學(xué) 城市地下空間工程北京市重點(diǎn)實(shí)驗(yàn)室，北京 100083； 3.泛華建設(shè)集團(tuán)有限公司，北京 100070)

地鐵隧道施工會(huì)破壞周圍地質(zhì)體原有的平衡，誘發(fā)應(yīng)力重分布，從而引起地層位移與地表不均勻沉降，這種現(xiàn)象在大斷面淺埋暗挖隧道施工中表現(xiàn)得尤為突出[1-4]。為保證隧道的安全施工，需對(duì)軟弱圍巖、自穩(wěn)能力差的隧道洞口或洞身段進(jìn)行預(yù)支護(hù)。主要的預(yù)支護(hù)方法有：超前錨桿法[5]、冷凍法[6-7]、管棚法[8]、水平旋噴注漿法[9]、超前小導(dǎo)管注漿法[10]，其中超前小導(dǎo)管注漿法具有錨桿和棚架雙重作用，其經(jīng)濟(jì)性、靈活性和良好的加固效果使其在隧道穿越淺埋、軟弱地層等不良地段的施工過(guò)程中發(fā)揮了重要作用[11-12]。因此，確定合理的超前小導(dǎo)管預(yù)注漿參數(shù)，以增強(qiáng)超前小導(dǎo)管的支護(hù)效果，改善掌子面的受力情況，使隧道開(kāi)挖過(guò)程中不同部位的地層應(yīng)力、位移等指標(biāo)都減少至最小，對(duì)工程順利開(kāi)展具有重要意義。

本文基于哈爾濱地鐵三號(hào)線湘江路站—會(huì)展中心站區(qū)間(湘—會(huì)區(qū)間)大斷面地鐵隧道暗挖段工程，采用FLAC 3D軟件獲得不同超前小導(dǎo)管預(yù)注漿參數(shù)下地鐵隧道開(kāi)挖引起的地表沉降縱向和橫向分布規(guī)律。然后設(shè)計(jì)正交試驗(yàn)，采用FLAC 3D對(duì)9組試驗(yàn)方案進(jìn)行模擬，以掌子面后方14 m處截面最大地表沉降和掌子面處地表沉降釋放率為評(píng)價(jià)指標(biāo)，進(jìn)一步分析不同超前小導(dǎo)管參數(shù)對(duì)地表沉降的影響程度，為后續(xù)類似地鐵隧道工程的設(shè)計(jì)和優(yōu)化提供參考。

1 工程概況

哈爾濱地鐵三號(hào)線湘—會(huì)區(qū)間自湘江路站起，由南向北沿紅旗大街穿行至?xí)怪行恼局梗L(zhǎng)為721.633 m，采用淺埋暗挖法施工。隧道橫斷面位置關(guān)系如圖1所示。大斷面隧道段為長(zhǎng)84 m、寬14.15 m、高10.62 m的馬蹄形隧道，隧道頂板埋深為10.90 m，采用雙側(cè)壁導(dǎo)坑法施工。大斷面隧道段土層的物理力學(xué)參數(shù)見(jiàn)表1，可見(jiàn)該段隧道范圍內(nèi)土體強(qiáng)度低，隧道在開(kāi)挖過(guò)程中出現(xiàn)塑性破壞的可能性較大，因此需要使用超前小導(dǎo)管對(duì)隧道進(jìn)行超前支護(hù)。

圖1 隧道橫斷面位置關(guān)系(單位：m)

表1 土層參數(shù)

2 隧道與小導(dǎo)管模型的建立

2.1 模型和邊界條件

根據(jù)工程經(jīng)驗(yàn)與圣維南原理，本模型沿隧道開(kāi)挖方向取84 m，水平方向取40 m，豎直方向取51 m，上邊緣取至地表。初期襯砌采用Liner結(jié)構(gòu)單元模擬，小導(dǎo)管采用Pile單元模擬，同一坐標(biāo)的不同結(jié)構(gòu)單元節(jié)點(diǎn)采用Link單元?jiǎng)傂赃B接。

模型的左右邊界、前后邊界以及底面施加法向位移約束，上部取自由邊界。由于隧道本身和地質(zhì)條件的對(duì)稱性，隧道的受力基本上也是對(duì)稱的，為了簡(jiǎn)化計(jì)算，利用FLAC 3D建立單側(cè)隧道模型，如圖2所示。

圖2 單側(cè)隧道模型

選取Mohr-Coulumb彈塑性模型進(jìn)行數(shù)值分析，采用應(yīng)力釋放法模擬支護(hù)施加前產(chǎn)生的變形。根據(jù)設(shè)計(jì)，小導(dǎo)管管間距離按2倍導(dǎo)管直徑選取，縱向隨鋼支撐逐榀打設(shè)，每榀鋼支撐間距為1 m。小導(dǎo)管布置模型如圖3所示。模型開(kāi)挖步長(zhǎng)為1 m，開(kāi)挖順序按照左右導(dǎo)洞對(duì)稱開(kāi)挖考慮，具體為：①側(cè)面上導(dǎo)洞開(kāi)挖；②側(cè)面上導(dǎo)洞開(kāi)挖至7 m后側(cè)面下導(dǎo)洞開(kāi)挖；③側(cè)導(dǎo)洞完成開(kāi)挖后中上洞開(kāi)挖；④中上洞開(kāi)挖至7 m后中下洞開(kāi)挖。

圖3 小導(dǎo)管布置模型

2.2 小導(dǎo)管與初期支護(hù)材料參數(shù)設(shè)置

超前小導(dǎo)管采用等效剛度的方法考慮，用慣性矩作為計(jì)算權(quán)重，初期支護(hù)與小導(dǎo)管材料的參數(shù)見(jiàn)表2。

表2 初期支護(hù)與小導(dǎo)管材料物理力學(xué)參數(shù)

3 小導(dǎo)管預(yù)注漿參數(shù)對(duì)地表沉降的影響

3.1 小導(dǎo)管預(yù)注漿參數(shù)設(shè)置

為了研究小導(dǎo)管沿開(kāi)挖方向投影長(zhǎng)度(簡(jiǎn)稱投影長(zhǎng)度)、注漿半徑及徑向加固范圍這3個(gè)參數(shù)對(duì)大斷面隧道開(kāi)挖地表沉降的影響，假設(shè)以下3種工況進(jìn)行數(shù)值模擬分析。

工況1：超前小導(dǎo)管注漿半徑和徑向加固范圍均為0.8 m，小導(dǎo)管投影長(zhǎng)度分別為3，4，5 m。

工況2：超前小導(dǎo)管投影長(zhǎng)度為4 m，徑向加固范圍為1.2 m，小導(dǎo)管注漿半徑分別為0.8，1.2，1.5 m。

工況3：超前小導(dǎo)管投影長(zhǎng)度為4 m，注漿半徑為1.8 m，小導(dǎo)管徑向加固范圍分別為0.8，1.2，1.5 m。

超前小導(dǎo)管徑向加固范圍即為超前小導(dǎo)管沿徑向投影長(zhǎng)度，其值應(yīng)與注漿半徑相等，但本工程采用超前小導(dǎo)管注漿和深孔注漿2種加固形式，因此，注漿半徑應(yīng)該大于徑向加固范圍。但如果注漿效果較差，則會(huì)出現(xiàn)注漿半徑小于徑向加固范圍的情況。

3.2 不同小導(dǎo)管預(yù)注漿參數(shù)下地表縱向沉降規(guī)律

根據(jù)以上3種工況模擬得到不同小導(dǎo)管預(yù)注漿參數(shù)下大斷面隧道側(cè)導(dǎo)洞與中洞開(kāi)挖完成后的地表沉降縱向分布曲線，見(jiàn)圖4。

圖4 不同小導(dǎo)管參數(shù)下地表沉降縱向分布曲線

由圖4可知：

1)地鐵隧道開(kāi)挖引起的地表沉降沿縱向呈倒S形分布，先期沉降集中發(fā)生于掌子面前方20 m內(nèi)，在掌子面前方10 m至掌子面區(qū)段迅速增大；掌子面后方約15 m后沉降速率開(kāi)始變緩，到30 m后呈現(xiàn)緩慢變化趨勢(shì)。掌子面處地表沉降量約占總沉降量的1/2。

2)側(cè)導(dǎo)洞開(kāi)挖完成時(shí)地表最大沉降位于掌子面后方35～40 m，中洞開(kāi)挖完成時(shí)最大地表沉降位于掌子面后方8～10 m，之后表現(xiàn)為回彈趨勢(shì)。這是因?yàn)橹卸撮_(kāi)挖后下部土體重力荷載降低，隨后初期支護(hù)閉合成環(huán)，土體回彈帶動(dòng)整條隧道上?。涣硗?，數(shù)值模擬中Mohr-Coulumb模型的彈性模量設(shè)為定值也會(huì)增大回彈現(xiàn)象。此外，超前加固條件越好，回彈表現(xiàn)越明顯。

3)隨著小導(dǎo)管投影長(zhǎng)度、注漿半徑及徑向加固范圍的增大，掌子面前方土體的沉降變化速率逐漸變緩，掌子面處地表沉降變小，占總沉降比例也逐漸減小。相比較而言，小導(dǎo)管投影長(zhǎng)度影響最大。

3.3 不同小導(dǎo)管預(yù)注漿參數(shù)下地表橫向沉降規(guī)律

圖5 不同小導(dǎo)管參數(shù)下地表沉降橫向分布曲線

如圖4所示，中洞開(kāi)挖完成時(shí)表現(xiàn)出的回彈趨勢(shì)與實(shí)際略有不符，所以，截取掌子面后方14 m處截面(此處縱向分布曲線尚未表現(xiàn)回彈趨勢(shì)或表現(xiàn)不明顯)得到不同小導(dǎo)管預(yù)注漿參數(shù)下大斷面隧道中洞開(kāi)挖完成時(shí)的地表沉降橫向分布曲線，見(jiàn)圖5。

1)由圖5(a)可知，小導(dǎo)管投影長(zhǎng)度分別為3，4，5 m 時(shí)對(duì)應(yīng)的掌子面后方14 m處截面的最大地表沉降分別為33.9，27.8，21.3 mm，沉降槽寬度分別為20，19，18 m。顯然增加小導(dǎo)管投影長(zhǎng)度對(duì)控制最大地表沉降有明顯改善作用，并使沉降槽寬度減小。

2)隨著注漿半徑及徑向加固范圍的增大，最大地表沉降逐漸減小，同時(shí)沉降槽寬度略有減小，但相比增加小導(dǎo)管投影長(zhǎng)度，增大注漿半徑及徑向加固范圍對(duì)控制地表沉降的改善效果不明顯。

4 基于正交試驗(yàn)的多因素敏感性分析

4.1 影響因素及正交試驗(yàn)設(shè)計(jì)

為了進(jìn)一步確定地鐵隧道開(kāi)挖時(shí)不同超前小導(dǎo)管預(yù)注漿參數(shù)對(duì)地表沉降的影響程度，選取小導(dǎo)管投影長(zhǎng)度、注漿半徑和徑向加固范圍作為主要敏感因素，每個(gè)因素選取3個(gè)水平進(jìn)行比較，見(jiàn)表3。

表3 敏感因素水平

在不考慮交互作用時(shí)，選用正交表L9(34)安排試驗(yàn)組合，設(shè)計(jì)9組試驗(yàn)方案，見(jiàn)表4。其中，方案6，8，9中，小導(dǎo)管徑向加固范圍小于注漿半徑，這與現(xiàn)場(chǎng)超前加固方法相對(duì)應(yīng)；而方案4，5，7中，小導(dǎo)管徑向加固范圍大于注漿半徑，可以對(duì)應(yīng)注漿效果較差的情況。

表4 試驗(yàn)方案及結(jié)果

圖6 不同方案地表沉降橫向分布曲線

圖7 不同方案地表沉降縱向分布曲線

采用FLAC 3D對(duì)各組試驗(yàn)方案進(jìn)行模擬。正交試驗(yàn)選取中洞開(kāi)挖完成時(shí)掌子面后方14 m處截面地表沉降橫向分布曲線(見(jiàn)圖6)，獲得該截面地表最大沉降值(參見(jiàn)表4)。而為了規(guī)避中洞開(kāi)挖模擬時(shí)明顯的回彈趨勢(shì)，選用導(dǎo)洞開(kāi)挖完成時(shí)的地表沉降縱向分布曲線(見(jiàn)圖7)，計(jì)算掌子面處地表沉降釋放率(參見(jiàn)表4)。

4.2 試驗(yàn)結(jié)果分析

4.2.1 極差分析

分別對(duì)掌子面后方14 m處截面最大地表沉降和掌子面處地表沉降釋放率結(jié)果進(jìn)行極差分析，結(jié)果見(jiàn)表5。可以看出，小導(dǎo)管投影長(zhǎng)度對(duì)掌子面后方14 m處截面最大地表沉降和掌子面處沉降釋放率影響最大，其次是注漿半徑，而徑向加固范圍的影響最小。

表5 極差分析結(jié)果

注：Rj為第j個(gè)因素各水平的綜合平均值的極差，j=(A,B,C)。

4.2.2 方差分析

為了確定試驗(yàn)結(jié)果差異來(lái)自于不同敏感因素的相應(yīng)水平還是來(lái)自試驗(yàn)誤差，取顯著水平α=0.1，對(duì)掌子面后方14 m處截面最大地表沉降和掌子面處地表沉降釋放率進(jìn)行方差分析，結(jié)果見(jiàn)表6。

表6 方差分析結(jié)果

由表6可知，小導(dǎo)管投影長(zhǎng)度和注漿半徑對(duì)掌子面后方14 m處截面最大地表沉降影響顯著；小導(dǎo)管徑向加固范圍方差大于誤差，但F值小于臨界值，因此對(duì)最大地表沉降有影響但不顯著。小導(dǎo)管投影長(zhǎng)度對(duì)掌子面處地表沉降釋放率影響顯著；注漿半徑方差大于誤差，但F值小于臨界值，因此注漿半徑對(duì)掌子面處地表沉降釋放率有影響但不顯著；小導(dǎo)管徑向加固范圍方差等于誤差，可視為誤差。

5 結(jié)論

1)以哈爾濱地鐵三號(hào)線湘—會(huì)區(qū)間為實(shí)例，采用FLAC 3D模擬得出不同的超前小導(dǎo)管預(yù)注漿參數(shù)下地表沉降縱向和橫向分布規(guī)律，發(fā)現(xiàn)改善小導(dǎo)管投影長(zhǎng)度、注漿半徑及徑向加固范圍可以有效地減小地表最大沉降量。

2)通過(guò)設(shè)計(jì)正交試驗(yàn)，以掌子面后方14 m處截面最大地表沉降和掌子面處地表沉降釋放率為評(píng)價(jià)指標(biāo)，通過(guò)極差分析和方差分析得出：小導(dǎo)管沿開(kāi)挖方向投影長(zhǎng)度對(duì)地表沉降的影響最為顯著，其次是注漿半徑，而徑向加固范圍的影響最小。