劉國(guó)楠，張遠(yuǎn)榮，肖文海

（中國(guó)鐵道科學(xué)研究院深圳研究設(shè)計(jì)院，廣東 深圳 518034）

0 前言

在沿海、沿江地區(qū)常遇到地下水量豐富且強(qiáng)透水性的砂層，盾構(gòu)開(kāi)挖會(huì)引起砂層孔隙水快速滲透和砂土屈服變形導(dǎo)致地面沉降。當(dāng)?shù)孛嬗薪ㄖ飼r(shí)，建筑物基礎(chǔ)下方的砂層受擾動(dòng)屈服后，荷載作用促使土體孔壓升高，進(jìn)一步加快孔隙水的滲流，對(duì)建筑物地基的影響問(wèn)題更為突出。研究富水砂層中的盾構(gòu)開(kāi)挖與地面建筑物沉降的關(guān)系，有助于在工程中合理選取盾構(gòu)施工參數(shù)，指導(dǎo)盾構(gòu)機(jī)施工操作，確定合理的地層預(yù)加固方案。

本文結(jié)合深圳地鐵2號(hào)線某工點(diǎn)工程實(shí)際，采用流固耦合水土模型和FLAC3D程序，用數(shù)值模擬的方法分析研究富水砂層盾構(gòu)施工對(duì)地面建筑物的影響問(wèn)題。取得的結(jié)果可指導(dǎo)預(yù)加固處理設(shè)計(jì)和盾構(gòu)施工。

1 工程實(shí)例

深圳地鐵2號(hào)線東角頭站至海上世界站區(qū)間，雙線隧道長(zhǎng)約1600m，場(chǎng)地原屬于海漫灘地貌，后填海形成陸地，現(xiàn)為城市鬧市區(qū)。該區(qū)間下有厚度8～10m的富水粗砂、礫砂層，地下水位在地面下1 m，水量豐富。地面建筑物密集，以天然地基的多層磚混住宅為主。本文選該區(qū)間的一段進(jìn)行盾構(gòu)施工數(shù)值模擬。

隧道為單圓盾構(gòu)襯砌結(jié)構(gòu)，管片外徑6.0m，厚度0.3m，環(huán)寬1.5m，襯砌分6塊，采用錯(cuò)縫拼裝。

本文研究段為地面建筑物重點(diǎn)保護(hù)段，布置地面和建筑物沉降監(jiān)測(cè)點(diǎn)如圖1所示。

圖1 模擬工點(diǎn)平面圖Fig.1 Layout plane of simulated project.

2 數(shù)值模擬方法

2.1 建立數(shù)值模型

2.1.1 幾何模型

選自東向西里程為ZK4＋340m～＋430段作為本次數(shù)值模擬對(duì)象。依據(jù)勘察資料，用FLAC3D軟件建立場(chǎng)地模型［1］，計(jì)算模型尺寸和地層構(gòu)造如圖2所示。定義X方向?yàn)闄M向，Y方向?yàn)樗淼揽v向，Z方向?yàn)樨Q向。

2.1.2 本構(gòu)模型和結(jié)構(gòu)單元

圖2 計(jì)算模型Fig.2 Computational model.

巖土體實(shí)體單元采用摩爾－庫(kù)倫模型；盾構(gòu)開(kāi)挖的土體單元采用空模型；盾構(gòu)鋼殼采用三維殼體單元；預(yù)制混凝土襯砌管片采用襯砌單元；間隙單元采用軟單元來(lái)模擬，其彈性模量定義為10kPa；漿液?jiǎn)卧獜椥阅Ａ侩S凝固時(shí)間增加取值從1MPa增至40MPa。

分析采用有效應(yīng)力法，并考慮水土的流－固相互作用。假定孔隙水滲流符合達(dá)西定律，可以得出孔隙水滲流的連續(xù)性表達(dá)式［2－3］：

式中：k為土的滲透系數(shù)；t為時(shí)間；γw為水的容重。結(jié)合土體運(yùn)動(dòng)平衡關(guān)系、變形幾何關(guān)系和有效應(yīng)力下的本構(gòu)關(guān)系可推導(dǎo)出位移和孔壓表示的滲流固結(jié)協(xié)調(diào)方程式：

式中：G為砂土的剪切模量；v為泊松比；p為孔隙水壓力；ρ為與飽和度相關(guān)的砂土密度；u、v和w 為單元位移函數(shù)和為單元位移速率函數(shù)。

2.1.3 邊界條件和初始條件

模型底面為不透水邊界且所有節(jié)點(diǎn)固定；模型左側(cè)和右側(cè)約束水平方向運(yùn)動(dòng)，正面和背面約束縱向運(yùn)動(dòng)；所有側(cè)面設(shè)定為孔壓不變的透水邊界。場(chǎng)地初始孔隙水壓地下1m為零，以下為靜水壓力。在模型正上方設(shè)定有一座7層樓高的建筑物，基礎(chǔ)面積6m×21m。基礎(chǔ)附近應(yīng)力7.5×104Pa。

2.1.4 材料參數(shù)

盾構(gòu)開(kāi)挖數(shù)值模擬部分材料采用的力學(xué)參數(shù)如 表1所示。

表1 模擬計(jì)算參數(shù)

2.2 場(chǎng)地初始模型

在上述模型的基礎(chǔ)上，求解出場(chǎng)地初始模型，此狀態(tài)代表盾構(gòu)施工前的孔壓和場(chǎng)地應(yīng)力狀態(tài)。場(chǎng)地地層的初始孔壓為靜水壓力，從地下1m處0Pa線性增大至3.1×105Pa。圖3是場(chǎng)地初始模型中心橫剖面有效應(yīng)力分布，在建筑基礎(chǔ)下方存在附加應(yīng)力擴(kuò)散形式。

圖3 中心橫剖面有效應(yīng)力Fig.3 Effective stress distribution on the cross profile of model center.

2.3 隧道開(kāi)挖模擬

對(duì)盾構(gòu)開(kāi)挖施工過(guò)程采取分節(jié)的方式模擬計(jì)算，每節(jié)開(kāi)挖并進(jìn)行流－固耦合計(jì)算后，通過(guò)相關(guān)單元的變更替換來(lái)模擬盾構(gòu)機(jī)推進(jìn)過(guò)程。參考盾構(gòu)開(kāi)挖相關(guān)技術(shù)見(jiàn)文獻(xiàn)［4］，本文模擬盾構(gòu)施工時(shí)對(duì)不同位置單元節(jié)點(diǎn)施加的主要荷載如圖4、圖5所示。

圖4 盾構(gòu)機(jī)荷載示意圖Fig.4 Sketch of load on the shield.

2.4 模擬結(jié)果對(duì)比

圖5 盾構(gòu)機(jī)施加荷載示意圖Fig.5 Sketch of shield applied load.

圖6標(biāo)識(shí)了地表沉降計(jì)算點(diǎn)、實(shí)測(cè)點(diǎn)與孔壓計(jì)算點(diǎn)位置，其中D1、D2、D3為地表建筑物計(jì)算沉降點(diǎn)；H1、H2、H3為地面沉降計(jì)算點(diǎn)；106、107為地表建筑觀測(cè)點(diǎn)；Y77、Y78為地面沉降實(shí)測(cè)點(diǎn)；A、B、C、E、F、G為孔壓計(jì)算點(diǎn)。

圖6 實(shí)測(cè)點(diǎn)與計(jì)算點(diǎn)位置Fig.6 Position of measuring pionts and calculation points.

D2點(diǎn)在盾構(gòu)施工過(guò)程中的沉降發(fā)展過(guò)程如圖7中的虛線所示。盾構(gòu)推進(jìn)時(shí)，該點(diǎn)先隆起，隆起量約10mm；盾構(gòu)開(kāi)挖通過(guò)時(shí)地表開(kāi)始產(chǎn)生沉降，沉降量約20～30mm；盾構(gòu)離開(kāi)后，地表沉降趨于穩(wěn)定。該變化規(guī)律與地面實(shí)測(cè)的沉降基本一致（圖7中實(shí)線）。上述對(duì)比說(shuō)明，本文采用的盾構(gòu)施工對(duì)地面沉降影響的數(shù)值模擬方法合理可信。

圖7 隧道計(jì)算與實(shí)測(cè)沉降Fig.7 Calculated and observed subsidences on the axis and the across section of the tunnel.

縱剖面上深度6m的B點(diǎn)和深度為10m的F點(diǎn)計(jì)算孔壓變化曲線如圖8所示。盾構(gòu)推進(jìn)時(shí)，掌子面孔隙水壓力先上升，盾構(gòu)完全通過(guò)時(shí)，孔隙水壓力消散，隧道管片安裝后，孔隙水壓力逐漸恢復(fù)且趨于穩(wěn)定?？克淼拦绊斘恢迷浇?，孔壓變化越大。

圖8 隧道縱剖面B和F點(diǎn)孔隙水壓力消散變化曲線Fig.8 Calculated pore pressures curves at point B and Fon the axis section of the tunnel.

3 對(duì)地表建筑物的影響分析

按圖9所示的天然地基建筑物位于隧道縱向中心線上方、距隧道一倍埋深、距隧道二倍埋深位置三種工況，分別進(jìn)行盾構(gòu)施工影響數(shù)值模擬分析。

取模型中間1／3作為計(jì)算域，在計(jì)算時(shí)選取部分觀測(cè)點(diǎn)如圖10所示。

3.1 工況一

盾構(gòu)開(kāi)挖時(shí)隧道中心縱剖面孔壓分布如圖11所示。盾構(gòu)頂推作用使得掌子面前方土體孔壓升高，掌子面排水導(dǎo)致盾構(gòu)上方拱頂土體孔壓下降。

圖9 三種計(jì)算工況建筑物位置示意圖Fig.9 Position Sketch of three buildings.

圖10 計(jì)算模型觀測(cè)點(diǎn)示意圖Fig.10 Sketch of measuring points of model.

圖11 盾構(gòu)開(kāi)挖時(shí)隧道中心縱剖面的孔壓分布圖Fig.11 Pore pressure distribution on the axis profile of the tunnel.

隧道中心縱剖面地表各點(diǎn)沉降變化曲線如圖12所示，從觀測(cè)點(diǎn)A6、A11和A16沉降變化曲線可以看出開(kāi)挖面到達(dá)Y＝45m處時(shí)，建筑物基礎(chǔ)先微傾斜；繼續(xù)開(kāi)挖到過(guò)72小時(shí)后，即基礎(chǔ)正下方隧道挖掘完成后，建筑物基礎(chǔ)沉降速率加大；經(jīng)過(guò)120個(gè)小時(shí)，即隧道挖掘結(jié)束后，兩端沉降速率變緩，沉降值相近。

隧道開(kāi)挖后地層產(chǎn)生變形的原因主要有以下幾方面：一是隧道開(kāi)挖后隧道圍巖應(yīng)力得到一定的釋放，產(chǎn)生徑向收縮，這也是隧道開(kāi)挖地層損失的重要原因；二是在富水砂層中由于隧道開(kāi)挖產(chǎn)生開(kāi)挖失水，砂土中孔隙水壓力減小，有效應(yīng)力增大，導(dǎo)致土骨架被壓縮，發(fā)生整體沉降變形；三是隧道開(kāi)挖對(duì)隧道周邊的土體產(chǎn)生擾動(dòng)，導(dǎo)致土體承載力下降，土體發(fā)生膨脹、擠密、甚至流動(dòng)等變形。

圖12 隧道中心縱剖面地表各點(diǎn)沉降曲線Fig.12 Ground settlement curves of some calculation points on the ground surface of axis profile of the tunnel.

如圖13所示，受建筑物基礎(chǔ)附加應(yīng)力的影響，由于基礎(chǔ)下方土體受到擾動(dòng)和外荷載共同作用導(dǎo)致基礎(chǔ)沉降較大，其沉降量是周邊地面沉降量的7～10倍?；A(chǔ)下沉和地表土體向中間靠攏。

圖13 模型中心橫剖面豎向位移等值線圖Fig.13 Z－disp.contours in cross section of the tunnel.

3.2 工況二

工況二建筑物基礎(chǔ)位置在隧道右側(cè)斜上方45°方向上，該方向恰是隧道開(kāi)挖后圍巖剪應(yīng)力主方向，另外基礎(chǔ)荷載產(chǎn)生的剪應(yīng)力主方向在左下方45°方向上。兩個(gè)荷載的剪應(yīng)力主方向相向使基礎(chǔ)與隧道之間的土體容易受剪切破壞，塑性破壞區(qū)域甚至貫通。圖14所示該部位的土體橫向和豎向變形較大，由于該部位土體向基礎(chǔ)左下方移動(dòng)，勢(shì)必?cái)D壓基礎(chǔ)左側(cè)土體和隧道正上方土體的向上隆起，這些部位土體的隆起值比右邊土體的隆起值大。因基礎(chǔ)偏離隧道，最大沉降量減小。

3.3 工況三

圖14 模型中心橫剖面橫向位移平面等值線圖Fig.14 X－disp.contours in the cross section of the tunnel.

由于工況三的建筑物基礎(chǔ)偏離隧道較遠(yuǎn)，相比工況一與工況二而言，基礎(chǔ)下方土體受隧道開(kāi)挖影響較小，基礎(chǔ)附近孔壓變化較小?；A(chǔ)下方土體附近的橫向位移值比工況一大，比工況二小。如圖15所示計(jì)算結(jié)果表明基礎(chǔ)豎向最大沉降值約5mm。

圖15 模型中心橫剖面豎向位移平面等值線圖Fig.15 Z－disp.contours in the cross section of the tunnel.

3.4 三種工況下影響比較分析

三種工況下，豎向位移分布規(guī)律的共同點(diǎn)是基礎(chǔ)下方土體沉降隨深度增加沉降逐漸減小，基礎(chǔ)中心的沉降比基礎(chǔ)邊緣的沉降大。

如圖16所示，工況一、工況二及工況三的基礎(chǔ)中心最大沉降分別為22.9mm、7.01mm 和5.3 mm。工況一基礎(chǔ)沉降大，基礎(chǔ)兩側(cè)附近土體隆起值也較大，工況二的基礎(chǔ)左側(cè)附近土體的隆起值比右側(cè)附近土體的隆起值大2～3mm。由于工況三的基礎(chǔ)中心離隧道中心線較遠(yuǎn)，盾構(gòu)開(kāi)挖對(duì)該工況基礎(chǔ)兩側(cè)土體影響較小，基礎(chǔ)兩側(cè)土體隆起值相近，量值較小。

圖16 隧道橫剖面上地表沉降曲線Fig.16 Ground settlements in the cross section of the tunnel.

如圖17所示，工況一基礎(chǔ)中心在隧道中心線的位置上，場(chǎng)地條件對(duì)稱，橫向位移量值很??；工況二的基礎(chǔ)在隧道右上斜45°方向上，基礎(chǔ)橫向位移較大，方向指向隧道中心縱剖面；工況三的基礎(chǔ)中心線遠(yuǎn)離隧道中心線，橫向位移介于工況一和工況二之間。

圖17 基礎(chǔ)中心下方地層橫向位移沿深度分布圖Fig.17 Lateral displacement curves with depth under the foundation.

如圖18所示，工況一基礎(chǔ)下方豎向沉降近23 mm，從地下6m開(kāi)始減小。從地下10m到12m，由于土體受隧道開(kāi)挖脫空的影響較大，土體豎向位移值逐漸增大。工況二和工況三的基礎(chǔ)下方豎向位移從基礎(chǔ)底部開(kāi)始逐漸減小，到地下5m處豎向沉降基本穩(wěn)定在2mm以內(nèi)。

圖18 基礎(chǔ)中心下方地層豎向位移沿深度分布圖Fig.18 Vertical displacement curves with depth under the foundation.

4 結(jié)語(yǔ)

根據(jù)深圳地鐵2號(hào)線盾構(gòu)過(guò)富水砂層施工實(shí)例，借鑒相關(guān)文獻(xiàn)［5］盾構(gòu)數(shù)值模擬經(jīng)驗(yàn)，用FLAC3D建立與工程相符的數(shù)值模型。在相對(duì)隧道的三個(gè)典型位置上分別設(shè)定同樣的建筑物，比對(duì)盾構(gòu)過(guò)富水砂層后建筑物受到的影響得到以下結(jié)論：

（1）在富水砂層上的建筑物下穿盾構(gòu)開(kāi)挖隧道時(shí)，由于砂層的強(qiáng)透水性和易液化的特性，建筑物工后沉降明顯，盾構(gòu)開(kāi)挖過(guò)程重點(diǎn)監(jiān)測(cè)建筑物下方砂層降水和孔壓變化，沿開(kāi)挖方向在建筑物上密集布置差異沉降監(jiān)測(cè)點(diǎn)。

（2）對(duì)于橫向偏離隧道中心一倍隧道埋深的建筑物，工后建筑物水平移動(dòng)最大，易受剪切破壞，盾構(gòu)開(kāi)挖下穿時(shí)，應(yīng)加強(qiáng)建筑物在垂直隧道開(kāi)挖方向上的差異沉降和地表的差異沉降的監(jiān)測(cè)。

（3）建筑物基礎(chǔ)中心離隧道中心橫向距離超過(guò)兩倍隧道埋深時(shí)，盾構(gòu)開(kāi)挖對(duì)建筑影響較小，重點(diǎn)監(jiān)測(cè)建筑物下方砂層降水即可。

［1］LAC3DFast Lagranggian Analysis of Continua User＇s Guide［M］.Intasca Consoulting Group Iinc，2008.

［2］陳仲頤 周景星，王洪瑾.土力學(xué)［M］.北京：清華大學(xué)出版社，1992.

［3］FLAC3DFluid－Mechanical Interacting Manual［M］.Intasca Consulting Group Linc，2008.

［4］陳饋，洪開(kāi)榮，吳學(xué)松.盾構(gòu)施工技術(shù)［M］.北京：人民交通出版社，2009.

［5］張海波，劉國(guó)楠，高俊合.盾構(gòu)近距離掘進(jìn)對(duì)橋梁樁基的影響分析［J］.鐵道建筑，2007，08：37－40.