朱連臣 王渭明

(1.中鐵二十五局集團(tuán)有限公司，510600，廣州；2.山東科技大學(xué)土木工程與建筑學(xué)院，266590，青島;3.中鐵二十五局集團(tuán)第五工程有限公司，266101,青島//第一作者，高級(jí)工程師)

由于地下工程施工受地下水影響較大，因此關(guān)于滲流場(chǎng)的研究尤為重要。文獻(xiàn)[1]結(jié)合Trasee數(shù)值方法，對(duì)防水壩中滲流場(chǎng)作用下巖體的受力狀態(tài)進(jìn)行了模擬研究。文獻(xiàn)[2-3]通過(guò)室內(nèi)實(shí)驗(yàn)、理論分析、數(shù)值計(jì)算和現(xiàn)場(chǎng)試驗(yàn)等手段，探討了復(fù)合地層盾構(gòu)施工在流固耦合作用下開挖面失穩(wěn)機(jī)理、破壞演化特性、地層內(nèi)力場(chǎng)的分布及地表沉降槽曲線的變化等問(wèn)題。文獻(xiàn)[4]充分考慮了土壓盾構(gòu)動(dòng)態(tài)施工全過(guò)程，建立了離散元模型，研究了隧道掘進(jìn)對(duì)上軟下硬地層的擾動(dòng)特征。文獻(xiàn)[5]通過(guò)數(shù)值模型和室內(nèi)試驗(yàn)，對(duì)隧道開挖面穩(wěn)定性上限進(jìn)行了分析研究，并提出相應(yīng)的注漿加固措施。

綜上所述，近年來(lái)國(guó)內(nèi)外學(xué)者對(duì)不同施工方法、不同地層，以及不同水土壓力采用不同算法等情況下的隧道圍巖、開挖面及地表建筑物的沉降、應(yīng)力場(chǎng)、孔隙水壓力場(chǎng)進(jìn)行了總結(jié)。但絕大多數(shù)工程未考慮富水流塑軟硬復(fù)合等復(fù)雜地層流固耦合效應(yīng)下的盾構(gòu)施工穩(wěn)定性及控制方法，計(jì)算體系尚不完整。本文以隧道洞周圍巖的孔隙水壓力場(chǎng)、位移場(chǎng)，以及盾構(gòu)掘進(jìn)引起的地表和河堤的沉降作為隧道穩(wěn)定性判定依據(jù)，利用有限差分軟件對(duì)富水礫石軟硬復(fù)合地層盾構(gòu)隧道施工進(jìn)行了詳細(xì)模擬分析。

1 工程概況

長(zhǎng)沙地鐵3號(hào)線烈士公園東站—絲茅沖站區(qū)間隧道(以下簡(jiǎn)為“烈絲區(qū)間隧道”)距離烈士公園東站約1.13 km，區(qū)間隧道右線邊線的東側(cè)距離瀏陽(yáng)河大堤約36 m。區(qū)間隧道線間距約為12 m，區(qū)間側(cè)穿瀏陽(yáng)河大堤段長(zhǎng)度約為72 m。具體位置關(guān)系如圖1所示。

圖1 長(zhǎng)沙地鐵3號(hào)線烈絲區(qū)間隧道與瀏陽(yáng)河位置關(guān)系圖

烈絲區(qū)間隧道側(cè)穿瀏陽(yáng)河大堤段采用盾構(gòu)法施工，為雙向隧道。盾構(gòu)機(jī)外徑為6.0 m，內(nèi)徑為5.4 m，每環(huán)襯砌管片厚度為0.3 m，寬為1.5 m，管片采用C50高強(qiáng)度混凝土。

2 數(shù)值模擬分析

2.1 模型的建立

本次數(shù)值模擬以長(zhǎng)沙地鐵3號(hào)線烈絲區(qū)間隧道側(cè)穿瀏陽(yáng)河大堤段為研究對(duì)象，模型橫向范圍自兩隧道盾構(gòu)中心軸線向左右兩側(cè)各取36 m；豎向范圍自地表至隧道拱頂平均埋深取14.2 m，自隧底向下取23.8 m；縱向范圍為側(cè)穿瀏陽(yáng)河大堤段,取72 m。因此模型的長(zhǎng)度、寬度和高度分別為72 m、96 m、44 m。

采用有限差分軟件FLAC 3D對(duì)隧道施工進(jìn)行三維建模分析。三維計(jì)算模型如圖2所示。

圖2 烈絲區(qū)間隧道三維數(shù)值模型

2.2 計(jì)算參數(shù)

結(jié)合勘察單位和設(shè)計(jì)單位提供的詳細(xì)地質(zhì)勘探報(bào)告及設(shè)計(jì)圖紙，確定了模型中各地層的物理力學(xué)參數(shù)和流體計(jì)算參數(shù)，詳見表1和表2。

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

表2 模型中流體計(jì)算參數(shù)

2.3 計(jì)算工況

為突出流固耦合效應(yīng)對(duì)隧道穩(wěn)定性的影響，在計(jì)算中分3種工況進(jìn)行模擬計(jì)算。其中，工況1為管片剛拼裝完成；工況2和3均為管片拼裝完成后，經(jīng)過(guò)一段時(shí)間的滲流作用，地下水已經(jīng)形成滲流場(chǎng)作用于管片上，這兩種工況應(yīng)考慮流固耦合效應(yīng)。

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

2.4.1 圍巖孔隙水壓力場(chǎng)分析

圖3為隧道開挖36 m后所處橫斷面的圍巖孔隙水壓力圖。由圖3a)和圖3b)可以清楚地看出，無(wú)論左線隧道還是右線隧道先行施工，所得到的圍巖孔隙水壓力分布形狀均較為相似。由圖3c)可知，當(dāng)左線隧道和右線隧道同時(shí)開挖后，所得到的圍巖孔隙水壓力形成以兩個(gè)隧道為中心的漏斗狀，且隧道中線處的孔隙水壓力呈中間高、兩側(cè)低的形式分布。

2.4.2 隧道洞周位移分析

對(duì)左、右線隧道洞周在兩種不同工況作用下的洞周不同位置的位移進(jìn)行數(shù)據(jù)提取整理發(fā)現(xiàn)，在工況1下，左線隧道拱頂、拱肩和拱腰處的位移均大于工況2下同位置的位移。由地下水上浮作用產(chǎn)生的隧道洞周圍巖的豎向位移值占總位移值的20%～30%，因此在實(shí)際工程中地下水對(duì)隧道洞周圍巖位移的影響不可忽略。

圖3 橫斷面圍巖孔隙水壓力場(chǎng)

2.4.3 地表位移分析

選取整體沉降量較大的右線隧道中心地表沉降值做豎向位移變化曲線，如圖4所示。

圖4將整個(gè)掘進(jìn)過(guò)程分為盾構(gòu)臨近期、盾構(gòu)通過(guò)期、盾尾脫出期和固結(jié)變形期等4個(gè)時(shí)期，通過(guò)對(duì)比分析發(fā)現(xiàn)，工況2時(shí)地表隆起值較工況1時(shí)更大，分析認(rèn)為產(chǎn)生這種變形的主要原因是由盾構(gòu)面的支護(hù)壓力對(duì)監(jiān)測(cè)面附近的土體產(chǎn)生了應(yīng)力擾動(dòng)導(dǎo)致的。盾構(gòu)通過(guò)期位移值由正變負(fù)，說(shuō)明地表在該時(shí)期處于下沉趨勢(shì)，產(chǎn)生這種變形的主要原因是由盾殼對(duì)監(jiān)測(cè)面附近土體的剪切摩擦作用導(dǎo)致的。在盾尾脫出期，兩種工況下地表的沉降值迅速增大，其中工況2的沉降值增速更大。固結(jié)變形期主要是監(jiān)測(cè)面的土體以及地下水位在受到擾動(dòng)變化后達(dá)到平衡狀態(tài)的恢復(fù)期，這部分變形緩慢進(jìn)行，產(chǎn)生變形的主要原因是由地層超孔隙水壓力消散產(chǎn)生的固結(jié)沉降和由土體蠕動(dòng)產(chǎn)生的蠕變沉降導(dǎo)致的。

3 隧道洞周圍巖加固方案及模擬分析

3.1 隨道洞周圍巖加固措施

左、右線隧道側(cè)穿瀏陽(yáng)河大堤區(qū)間的地層多為富水流塑性地層，提出采用φ800 mm三重管旋噴樁密排的加固措施。加固范圍為：水平加固至隧道結(jié)構(gòu)輪廓線外1.0 m，豎直加固深度為隧道頂以上3.0 m、隧道底以下2 m，局部淤泥質(zhì)地層較厚的區(qū)域應(yīng)相應(yīng)增加加固深度1.0～2.0 m。

3.2 加固方案的模擬工況

工況4為隧道洞周圍巖采取加固措施，不考慮流固耦合作用；工況5為隧道洞周圍巖采取加固措施，考慮流固耦合作用。

3.3 加固前后模擬結(jié)果對(duì)比

3.3.1 地表位移分析

圖5為工況4和工況5下地表沉降對(duì)比圖。由圖5可知，兩種工況下整個(gè)隧道斷面的沉降趨勢(shì)基本一致，呈對(duì)稱分布，其中隧道中線處地表沉降最大，沿隧道兩側(cè)沉降值逐漸減小。

圖5 兩種工況下地表沉降值對(duì)比圖

3.3.2 隧道洞周圍巖主應(yīng)力分析

不同工況下隧道拱頂、拱腰、拱腳、拱底等關(guān)鍵部位圍巖的最大主應(yīng)力和最小主應(yīng)力，如圖6所示。

由圖6可知，各工況下隧道洞周圍巖的最大和最小主應(yīng)力曲線趨勢(shì)較為相似，均以隧道中線為對(duì)稱軸呈左右對(duì)稱分布。其中，工況2主應(yīng)力值較工況4、工況5整體變化更加陡峭，這說(shuō)明該工況下隧道洞周圍巖應(yīng)力場(chǎng)中受到孔隙水壓力的沖擊比較大。

圖6 隧道洞周圍巖主應(yīng)力分布曲線

4 地表沉降實(shí)測(cè)數(shù)據(jù)與模擬數(shù)據(jù)對(duì)比分析

選取右側(cè)隧道地表沉降現(xiàn)場(chǎng)實(shí)測(cè)值與工況4、工況5的地表沉降模擬值進(jìn)行對(duì)比分析，如圖7所示。

圖7 地表沉降實(shí)測(cè)值與模擬值對(duì)比圖

由圖7可知，各工況下，地表沉降模擬值曲線的整體走勢(shì)和實(shí)測(cè)值基本相符，均沿隧道中線左右對(duì)稱，且整個(gè)沉降槽呈正態(tài)分布，從而驗(yàn)證了數(shù)值模擬結(jié)果符合工程實(shí)際規(guī)律。

5 結(jié)論

本文采用有限差分軟件FLAC 3D，對(duì)長(zhǎng)沙地鐵3號(hào)線烈絲區(qū)間側(cè)穿瀏陽(yáng)河大堤段盾構(gòu)隧道進(jìn)行了三維數(shù)值模擬，得到以下結(jié)論：

1)考慮流固耦合效應(yīng)時(shí)，盾構(gòu)開挖面孔隙水壓力場(chǎng)呈漏斗狀分布，孔隙水的流動(dòng)導(dǎo)致圍巖應(yīng)力場(chǎng)和位移場(chǎng)受到影響，與工況2數(shù)據(jù)對(duì)比，具有明顯的區(qū)別。

2)盾構(gòu)開挖、襯砌管片安裝和壁后注漿完成后，圍巖洞周在拱頂出現(xiàn)沉降，在拱底處出現(xiàn)隆起。對(duì)比工況1和工況2，由于地下水的上浮作用，隧道拱頂、拱肩和拱腰處的位移值在流固耦合作用下減小，隧道拱底和拱腳處的位移值在流固耦合作用下增大。

3)在盾構(gòu)臨近期，開挖面土體受到支護(hù)應(yīng)力的作用，產(chǎn)生應(yīng)力擾動(dòng)和地下水位下降，從而導(dǎo)致地表隆起。瀏陽(yáng)河河堤路基的最大沉降值為2.39 mm，地表沉降不符合規(guī)范值。而施作三重管旋噴樁加固地層后，沉降值符合規(guī)范要求，說(shuō)明該加固措施效果明顯。

4)盾構(gòu)隧道左、右兩線整體沉降模擬值和實(shí)測(cè)值趨勢(shì)基本一致，整體沉降曲線以隧道中線為中心呈對(duì)稱分布，但模擬值普遍偏小于實(shí)測(cè)值。