張子真 孫玉輝 陳昌彥 王金明
(北京市勘察設(shè)計研究院有限公司,100038,北京//第一作者,工程師)
盾構(gòu)隧道穿越粉細砂層時的地表沉降分析
張子真 孫玉輝 陳昌彥 王金明
(北京市勘察設(shè)計研究院有限公司,100038,北京//第一作者,工程師)
當盾構(gòu)拱頂遇到粉細砂層,地表沉降累積值和速率均易超標。通過繪制單線盾構(gòu)隧道地表沉降空間分布圖,發(fā)現(xiàn)沉降大的區(qū)域與粉細砂層有明顯的相關(guān)性。原因是粉細砂的粘聚力較小,自穩(wěn)性差,容易受擾動坍塌,引起較大超挖。采用有限元軟件Plaxis進行數(shù)值模擬。結(jié)果表明,相同盾構(gòu)掘進參數(shù)下,含粉細砂區(qū)域的地層地表沉降明顯偏大。
地鐵; 盾構(gòu)掘進; 粉細砂層; 地表沉降; 數(shù)值模擬
Author′s address Beijing Geotechnical Institute Engineering Consultants Co.Ltd.,100038,Beijing,China
北京地鐵16號線永豐站—永豐南站區(qū)間采用盾構(gòu)法施工。通過監(jiān)測發(fā)現(xiàn),該區(qū)間某些區(qū)域的沉降值和速率較大,連續(xù)幾天出現(xiàn)紅色預警,橙色預警更是頻繁。其中第55排監(jiān)測斷面的4個監(jiān)測點,同一天的沉降速率均超過了20 mm/d。
分析發(fā)現(xiàn),紅色預警居多的位置,盾構(gòu)拱頂上方約有1~3 m厚的粉細砂層。筆者認為,這種粉細砂層的地質(zhì)條件,決定了盾構(gòu)掘進引發(fā)的地表沉降量較大。
北京地鐵16號線永豐站—永豐南站區(qū)間,起自位于北清路路中、規(guī)劃永豐西二街西側(cè)的永豐站。該區(qū)間線路出永豐站后由西向東敷設(shè),在北清路與永豐路的交叉處以350 m的曲線半徑轉(zhuǎn)角,穿越五一渠后到達永豐路正下方;再沿永豐路由北向南敷設(shè),到達位于永豐路路中、規(guī)劃永豐南一街及永豐南三街之間的永豐南站。區(qū)間平面圖如圖1所示,地質(zhì)縱剖圖如圖2所示。拱頂土層主要為②粉土、②1粉質(zhì)黏土、③1粉質(zhì)黏土、②3粉細砂、③3粉細砂。
圖1 盾構(gòu)區(qū)間總平面圖
本文只考慮單線盾構(gòu)隧道掘進引發(fā)的地表沉降。由于右線先于左線掘進,故只考慮右線盾構(gòu)掘進所引起的地表沉降,并在收集沉降數(shù)據(jù)時盡量篩除左線盾構(gòu)的影響。從圖1和圖2可以看出,右線盾構(gòu)掘進500 m后,開始由轉(zhuǎn)彎過渡為平直段,在掘進到1 000 m之前,盾構(gòu)拱頂有3次要進入到粉細砂層中,之后盾構(gòu)將過渡到接收階段。盾構(gòu)始發(fā)、接收和轉(zhuǎn)彎時,盾構(gòu)機姿態(tài)控制難度加大,對地表沉降有加劇作用。為避開這些因素的影響,本文選擇出永豐站后的里程500~1 000 m為研究對象,分析盾構(gòu)平直段正常掘進階段,地層對地表沉降的影響。
圖2 盾構(gòu)區(qū)間地質(zhì)縱剖圖
以里程為橫軸,以右線盾構(gòu)掘進所引起的右線正上方地表最終沉降量為縱軸,繪制沉降量分布圖如圖3所示。結(jié)合地質(zhì)縱剖圖分析,可以發(fā)現(xiàn)地表沉降和地質(zhì)條件的相關(guān)性:當盾構(gòu)拱頂在粉細砂層中掘進時,相應位置的地表沉降要明顯增大。
圖3 沉降分布圖
為便于分析,圖3盾構(gòu)拱頂上方,沉降較大的三個區(qū)域,自北向南分別命名為A區(qū)、B區(qū)和C區(qū);其下對應的盾構(gòu)拱頂在粉細砂層中掘進的三個區(qū)域,分別命名為I區(qū)、II區(qū)和III區(qū);鄰近的粉質(zhì)黏土層,分別命名為甲區(qū)、乙區(qū)、丙區(qū)和丁區(qū)。
地表沉降A(chǔ)區(qū)基本上以粉細砂層I區(qū)為中心,最大沉降量為27.16 mm,盾構(gòu)掌子面上端2 m厚為③3粉細砂層,下端4 m厚為③1粉質(zhì)黏土層。根據(jù)勘察報告:③3粉細砂層中密或密實,飽和,粘聚力為0,內(nèi)摩擦角為25°,壓縮模量Es1-2=15 MPa,土層修正后圍巖分級為VI級,土石可挖性為II級;在地下水作用下易坍塌,無法形成應力拱,圍巖穩(wěn)定性差。相對于粉質(zhì)黏土,這種土用手可輕易挖掘。盾構(gòu)機刀盤邊緣線速度很快,對粉細砂層來說是一個巨大的干擾源,很容易使拱頂上方的粉細砂大量散落,形成空洞。由于其無法像粉質(zhì)黏土那樣長期維持應力拱,以致快速坍塌,造成地表沉降過大。
③1粉質(zhì)黏土層甲區(qū)和乙區(qū)上方的沉降則明顯減少了許多。如甲區(qū)粉質(zhì)黏土層上方的沉降最大為12.16 mm,約占A區(qū)地表最大沉降量的44.8%;乙區(qū)粉質(zhì)黏土層上方的沉降最大為17.35 mm,約占A區(qū)地表最大沉降量的63.9%。可見,粉質(zhì)黏土對于盾構(gòu)沉降控制是有利的。
地表沉降B區(qū)基本上以粉細砂層II區(qū)為中心,最大沉降量為30.17 mm,盾構(gòu)掌子面上端3 m厚為③3粉細砂層,大約一半在粉細砂層中,一半在粉質(zhì)黏土中。在這一區(qū)域,盾構(gòu)刀盤對粉細砂層產(chǎn)生了明顯的擾動,導致超挖嚴重,最終引起了較大的地表沉降。
③1粉質(zhì)黏土層乙區(qū)和丙區(qū)上方的沉降則明顯減少了許多。如乙區(qū)粉質(zhì)黏土層上方的沉降最大為17.35 mm,約占B區(qū)地表最大沉降量的57.5%;丙區(qū)粉質(zhì)黏土層上方的沉降最大為20.20 mm,約占B區(qū)地表最大沉降量的67.0%。
和A區(qū)相比,B區(qū)的最大沉降量為30.17 mm,大于A區(qū)最大沉降量27.16 mm,而II區(qū)穿越③3粉細砂層厚度為3 m,也大于I區(qū)穿越③3粉細砂層厚度2 m??梢?盾構(gòu)掌子面上端掘進通過粉細砂層的厚度越大,引發(fā)的地表沉降值也越大。
地表沉降C區(qū)基本上以粉細砂層III區(qū)為中心,最大沉降量為27.13 mm,盾構(gòu)掌子面上端2 m厚為③3粉細砂層,下端4 m厚為③1粉質(zhì)黏土層。在這一區(qū)域,盾構(gòu)刀盤對粉細砂層產(chǎn)生了明顯的擾動,導致超挖嚴重,最終引起了較大的地表沉降。
③1粉質(zhì)黏土層丙區(qū)和丁區(qū)上方的沉降則要明顯減少許多。如丙區(qū)粉質(zhì)黏土層上方的沉降最大為20.20 mm,約占C區(qū)地表最大沉降量的74.5%;丁區(qū)粉質(zhì)黏土層上方的沉降最大為14.04 mm,約占C區(qū)地表最大沉降量的51.8%。
和B區(qū)相比,C區(qū)的最大沉降量為27.13 mm,小于B區(qū)最大沉降量30.17 mm,與A區(qū)的最大沉降量27.16 mm相當;而III區(qū)穿越③3粉細砂層厚度為2 m,也小于II區(qū)穿越③3粉細砂層厚度3 m,與I區(qū)穿越③3粉細砂層厚度2 m相當。這說明盾構(gòu)掌子面上端掘進通過的粉細砂層厚度越大,地表沉降值也會越大。
為了更好地理解盾構(gòu)掘進時粉細砂層對地表沉降的影響,利用有限元軟件Plaxis 8.5進行數(shù)值模擬。隧道開挖對于土體來說是一個卸載的過程,如采用常規(guī)的摩爾-庫侖模型計算,由于其無法考慮卸載再壓縮模量與首次壓縮模量的差異,容易得出失真的結(jié)果。采用應變硬化模型(又稱Hardening Soil模型或HS模型),可以很好地解決這一問題。應變硬化模型既可以模擬軟土,也可以模擬硬土。
其屈服準則是摩爾-庫侖準則,并加入了一個屈服蓋帽。但與普通摩爾-庫侖模型不同的是,該模型首次壓縮和卸載再壓縮采用的是不同的模量,這和實際情況更加接近。該模型共有11個參數(shù),其中m、E50,ref、Eoed,ref、Eur,ref、vur、pref是摩爾-庫侖模型所沒有的新參數(shù)。m為本構(gòu)關(guān)系式的冪次,一般取0.4~0.7,根據(jù)軟件手冊建議,本次取0.5;E50,ref是三軸壓縮試驗在100 kPa圍壓下,偏應力與軸向應變試驗曲線中50%極限荷載所對應的割線模量;Eoed,ref為固結(jié)試驗中,荷載σ1為100 kPa所對應的切線模量;Eur,ref為三軸壓縮過程中,卸載再加載的割線模量;vur為卸載再加載泊松比,一般取0.2;pref為參考應力,一般取100 kPa。
由于目前的勘察工作只進行固結(jié)不排水三軸剪切試驗,很少進行固結(jié)排水三軸剪切試驗,且勘察報告只提供Es1-2,不提供各種土工試驗的原始數(shù)據(jù),因此無法直接從勘察報告中得到關(guān)鍵的模量參數(shù)。為此筆者查閱了相關(guān)文獻,發(fā)現(xiàn)國內(nèi)目前的研究中,文獻[1]對上海地區(qū)的四種土進行了常規(guī)三軸固結(jié)排水試驗,并計算了E50,ref、Eoed,ref、Eur,ref等與Es1-2的比例關(guān)系。本文在模擬中假設(shè)永豐站—永豐南站區(qū)間所遇到的③1粉質(zhì)黏土和③3粉細砂與文獻[1]所研究的土服從相同的規(guī)律,并根據(jù)其研究成果,推算出③1粉質(zhì)黏土和③3粉細砂可能的E50,ref、Eoed,ref、Eur,ref,輸入軟件計算,進行定性分析。
本文進行數(shù)值模擬所建立的2D模型(如圖4),寬度為400 m,厚度為75 m,盾構(gòu)直徑為6 m,管片厚度為0.35 m,中心埋深為30 m,上層土和下層土的分界面埋深30 m;根據(jù)勘察報告,取地下水位埋深為7 m。模型左、右邊界限制法向位移,底部固定。為簡化計算,本次模擬只考慮③1粉質(zhì)黏土和③3粉細砂兩種土層。土的各項參數(shù)如表1所示。表中,c為粘聚力,φ為內(nèi)摩擦角,k為滲透系數(shù),e為孔隙比,γ為濕重度,γsat為飽和重度。
表1 土體參數(shù)
幾何模型建好后,采用軟件自動程序生成網(wǎng)格。網(wǎng)格密度采用“中等”,如圖5所示。網(wǎng)格劃分后,可生成初始條件。將地下水位設(shè)定在-7 m的位置,并生成初始靜水壓力,如圖6所示。然后生成初始應力,如圖7所示。接下來進行計算。為了比較相同掘進條件、不同地層條件下地表沉降的差異情況,先后進行3次模擬。每次模擬在建立幾何模型時,按照圖4和表2對各土層進行定義。
圖4 幾何模型
圖5 網(wǎng)格劃分
圖6 初始靜水壓力
圖7 初始應力
表2 土層定義
計算步驟為:
第一步:激活“隧道管片”,挖掘管片內(nèi)的土,排干管片內(nèi)的地下水,重新計算孔隙水壓力。
第二步:雙擊“隧道中心”,設(shè)置超挖率為2%,相當于盾構(gòu)直徑減少了12 cm,以此模擬盾構(gòu)的超挖現(xiàn)象。
第三步:點擊“計算”按鈕,開始計算。
計算完成后,輸出地表沉降值,繪制沉降槽曲線,如圖8所示。
從圖8中可以看出,在相同盾構(gòu)施工條件下(超挖率),地層情況對隧道上方的沉降槽具有顯著的影響。
第1次模擬中,整個模型都是均一的③1粉質(zhì)黏土,此時沉降量較小,最大值為16.48 mm。
第2次模擬中,上層土為③3粉細砂,下層土為③1粉質(zhì)黏土,沉降明顯增大,最大值為21.01 mm。第1次模擬的最大沉降量約占第2次模擬的最大沉降量的78.4%。
圖8 沉降槽
第3次模擬中,整個模型都是均一的③3粉細砂,此時沉降最為劇烈,最大值為33.07 mm。第1次模擬的最大沉降量約占第3次模擬的最大沉降量的49.8%。
從3次模擬的結(jié)果來看,盾構(gòu)拱頂進入粉細砂層中掘進,可以明顯增大地表沉降趨勢。而盾構(gòu)掌子面進入到粉細砂層的厚度越厚,沉降增大的趨勢也越明顯。這與圖3觀察到的趨勢是一致的。
盾構(gòu)拱頂在粉細砂層中掘進,土層易受刀盤的擾動,粉細砂易散落,產(chǎn)生明顯的超挖,加大地表下沉。數(shù)值模擬結(jié)果表明,在相同超挖量的情況下,粉細砂層傾向于產(chǎn)生比粉質(zhì)黏土更大的地表沉降。在這兩個因素的共同作用下,盾構(gòu)拱頂在粉細砂層中掘進所引發(fā)的地表沉降,比在粉質(zhì)黏土中掘進時大。盾構(gòu)掘進過程中要下穿粉細砂層時,應提前做好掘進參數(shù)調(diào)整,加強同步注漿和二次注漿,控制地表沉降。同時加強現(xiàn)場監(jiān)測和巡視,使道路變形趨勢處于可控狀態(tài)。
[1] 王衛(wèi)東.基坑開挖數(shù)值分析中土體硬化模型參數(shù)的試驗研究[J].巖土力學,2012,33(5):2283-2290.
Analysis of Surface Subsidence with TBM Boreing through Fine Sand Stratum
ZHANG Zizhen, SUN Yuhui, CHEN Changyan, WANG Jinming
Boring through the fine sand stratum with TBM (tunnel boring machine) will cause great and fast surface subsidence,which may exceed the limits given by designers.By drawing the spatial distribution of surface subsidence for single shield tunnel,the obvious correlation between large area of surface subsidence and the fine sand stratum is detected,because the fine sand is vulnerable to collapse due to its weak cohesion,which results in unnecessary over-excavation.Numerical simulation with finite element software Plaxis also indicates that fine sand stratum is more susceptible to surface subsidence than clay layer when the shield driving parameters are the same.
metro; shield driving; fine sand stratum; surface subsidence; numerical simulation
TU433:U455.43
10.16037/j.1007-869x.2017.04.025
2015-05-18)