考慮地層三維空間分布的盾構(gòu)掘進(jìn)模擬

張超 ，徐智文 ，劉飛香 ，陳仁朋 ?，周蘇華 ，廖金軍 ，胡偉飛

（1.湖南大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙 410082；2.建筑安全與節(jié)能教育部重點(diǎn)實(shí)驗(yàn)室（湖南大學(xué)），湖南 長(zhǎng)沙 410082；3.中國鐵建重工集團(tuán)股份有限公司，湖南 長(zhǎng)沙 410100；4.浙江大學(xué) 機(jī)械工程學(xué)院，浙江 杭州 310027）

盾構(gòu)法是城市地下空間開發(fā)的主要施工方法，被廣泛應(yīng)用于城市地下交通隧道建設(shè).盾構(gòu)掘進(jìn)本質(zhì)上是盾構(gòu)機(jī)與地層動(dòng)態(tài)耦合作用的過程.地質(zhì)環(huán)境復(fù)雜多變導(dǎo)致盾構(gòu)機(jī)服役性態(tài)難預(yù)測(cè)，是盾構(gòu)掘進(jìn)施工風(fēng)險(xiǎn)的主要誘因，極易引起地表坍塌、開挖面失穩(wěn)、盾構(gòu)機(jī)姿態(tài)偏移和刀盤結(jié)泥餅等工程事故.因此，有必要開展盾構(gòu)掘進(jìn)模擬以獲取施工過程中裝備受力狀態(tài)和地層擾動(dòng)規(guī)律［1-5］.然而，現(xiàn)有盾構(gòu)掘進(jìn)模擬通常將地質(zhì)環(huán)境簡(jiǎn)化為平行成層地層［6-9］，難以反映盾構(gòu)服役地質(zhì)環(huán)境的三維空間變異性.

復(fù)雜地質(zhì)環(huán)境的定量表征一般通過三維地質(zhì)建模實(shí)現(xiàn).三維地質(zhì)建?？商峁┑貙有畔⒌膸缀文Ｐ?，需要借助額外的算法進(jìn)一步生成有限元網(wǎng)格.然而，地質(zhì)模型的三維空間變異性給相關(guān)算法研究帶來挑戰(zhàn).目前，一些學(xué)者利用地質(zhì)構(gòu)造面、地層界面、軟弱夾層等特殊界面對(duì)地質(zhì)模型進(jìn)行分塊處理，并對(duì)獨(dú)立封閉塊體進(jìn)行網(wǎng)格劃分，成功構(gòu)建了復(fù)雜地層數(shù)值模型［10-13］.此類地層數(shù)值模型主要應(yīng)用于邊坡穩(wěn)定性分析、洞室穩(wěn)定性分析等巖土工程問題［14-16］.目前，尚無三維地質(zhì)建模技術(shù)在盾構(gòu)掘進(jìn)模擬中應(yīng)用的相關(guān)報(bào)道.

本文提出了一種考慮地層三維空間分布的盾構(gòu)掘進(jìn)模擬方法.首先利用Kriging 算法基于鉆孔數(shù)據(jù)建立三維地質(zhì)模型，提出一種基于地層分界面的條件判斷算法，克服了傳統(tǒng)方法中復(fù)雜地層曲面帶來的有限元網(wǎng)格畸變等問題，實(shí)現(xiàn)三維復(fù)雜地層的高質(zhì)量網(wǎng)格劃分；然后，編寫了相應(yīng)的復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序，實(shí)現(xiàn)復(fù)雜地層盾構(gòu)掘進(jìn)模擬的參數(shù)化建模，大幅提高建模效率；最后，以長(zhǎng)沙地鐵四號(hào)線某區(qū)間為例，討論了所提方法的適用性.

1 三維地質(zhì)建模

三維地質(zhì)模型可定量表征復(fù)雜地層的地質(zhì)信息，為構(gòu)建復(fù)雜地層有限元模型提供幾何模型.三維地質(zhì)建模包括地表及地層分界面的高程插值計(jì)算和模型可視化兩部分.

1.1 Kriging算法

在地質(zhì)統(tǒng)計(jì)學(xué)理論中，Kriging 算法是對(duì)區(qū)域化隨機(jī)變量進(jìn)行無偏最優(yōu)估計(jì)的空間插值算法.該算法在礦產(chǎn)空間分布以及地質(zhì)模型構(gòu)建方面具有廣泛應(yīng)用.下面以高程插值為例，介紹該算法的主要插值過程：將待插值點(diǎn)X0的高程定義為Z*（X0），且已知n個(gè)鉆孔樣本的數(shù)據(jù)，并將第i個(gè)鉆孔與地表或地層分界面的交點(diǎn)Xi的高程定義為Z（Xi），則Z*（X0）可表示為n個(gè)Z（Xi）的加權(quán)求和［17］，見式（1）.

式中：λi（i=1，2，…，n）為權(quán)重系數(shù).Kriging 插值過程即為通過理論變異函數(shù)γ（h）構(gòu)建的Kriging 方程組［18］求解最優(yōu)權(quán)系數(shù)的過程，方程組見式（2）.

式中：γ（hij）是由交點(diǎn)Xi、Xj的水平間距hij確定的理論變異函數(shù)值；γ（hi0）是由交點(diǎn)Xi和插值點(diǎn)X0的水平間距hi0確定的理論變異函數(shù)值；μ是拉格朗日乘子.此處，理論變異函數(shù)是兩點(diǎn)水平間距與高程變異性的函數(shù)，可通過擬合實(shí)驗(yàn)變異函數(shù)值獲取，見圖1.實(shí)驗(yàn)變異函數(shù)值則由鉆孔樣本數(shù)據(jù)計(jì)算得到.球狀變異函數(shù)和實(shí)驗(yàn)變異函數(shù)值分別見式（3）和式（4）.

圖1 球狀變異函數(shù)模型Fig.1 Spherical variogram model

式中：γ（h）為理論變異函數(shù)；C0為塊金常數(shù)；C為拱高；a為變程；h為樣本點(diǎn)間水平間距.

式中：γ*（h）為實(shí)驗(yàn)變異函數(shù)；N（h）表示水平間距為h的鉆孔樣本對(duì)數(shù)量；Xi與Xi+h分別表示上述鉆孔樣本對(duì)與地表或地層分界面的交點(diǎn).

如此，Kriging 算法首先由鉆孔樣本數(shù)據(jù)得到離散的實(shí)驗(yàn)變異函數(shù)值，并利用球狀模型等數(shù)學(xué)模型擬合離散值得到理論變異函數(shù)，再將理論變異函數(shù)代入Kriging 方程組解得插值點(diǎn)的最優(yōu)權(quán)重系數(shù)λi，最后根據(jù)公式（1）得到該點(diǎn)的高程.

1.2 建模流程

三維地質(zhì)建模流程具體包括：鉆孔樣本預(yù)處理、變異函數(shù)求解、地表及地層分界面高程計(jì)算、三維地質(zhì)模型可視化共四個(gè)步驟.

1.2.1 鉆孔樣本預(yù)處理

鉆孔樣本預(yù)處理的目的是保證所有鉆孔的地層類別和層序一致，以便采用統(tǒng)一數(shù)據(jù)格式對(duì)鉆孔樣本的地層信息進(jìn)行存儲(chǔ)和調(diào)用.根據(jù)所有鉆孔樣本提供的地層信息，對(duì)每個(gè)鉆孔揭露的地層按照沉積順序進(jìn)行編號(hào)，見圖2.當(dāng)鉆孔樣本缺失某地層時(shí)，在地層缺失處插入厚度為零的該地層，從而保證所有鉆孔的地層類別和層序相同.

圖2 鉆孔樣本預(yù)處理Fig.2 Borehole sample pretreatment

1.2.2 變異函數(shù)求解

基于預(yù)處理后的鉆孔樣本，采用公式（4）計(jì)算地表、各地層分界面的實(shí)驗(yàn)變異函數(shù)γ*（h），并利用公式（3）對(duì)離散的實(shí)驗(yàn)變異函數(shù)值進(jìn)行擬合，得到理論變異函數(shù)γ（h）.

1.2.3 地表及地層分界面高程計(jì)算

對(duì)地表或地層分界面進(jìn)行離散，并以離散的待插值點(diǎn)的水平投影作為網(wǎng)格節(jié)點(diǎn)，形成圖3 所示虛擬背景網(wǎng)格.以點(diǎn)X0高程插值為例，根據(jù)點(diǎn)X0、Xi和Xj的投影位置得到hi0和hij，代入理論變異函數(shù)γ（h）計(jì)算得γ（hi0）和γ（hij），再代入Kriging 方程組求解得到權(quán)重系數(shù)λi，求解出X0高程.按上述過程可得到地表及地層分界面上所待插值點(diǎn)的高程數(shù)據(jù).

圖3 地層分界面的虛擬背景網(wǎng)格Fig.3 Virtual background grid of stratum interface

1.2.4 三維地質(zhì)模型可視化

將地表及地層分界面上插值點(diǎn)的平面坐標(biāo)和高程數(shù)據(jù)導(dǎo)入自主搭建的可視化平臺(tái)，得到三維地質(zhì)模型，如圖4所示.

圖4 三維地質(zhì)模型示意圖Fig.4 3D geological model diagram

2 復(fù)雜地層網(wǎng)格劃分

由于復(fù)雜地層的強(qiáng)空間變異性，構(gòu)建的三維地質(zhì)模型通常具有不規(guī)則的幾何形態(tài)，難以直接生成高質(zhì)量有限元網(wǎng)格.因此，有必要針對(duì)復(fù)雜地層開發(fā)高效有限元網(wǎng)格劃分算法.

復(fù)雜地層盾構(gòu)掘進(jìn)模擬的有限元網(wǎng)格劃分難點(diǎn)主要體現(xiàn)在兩個(gè)方面.一方面，直接基于地質(zhì)模型的網(wǎng)格劃分方法在地層尖滅處易導(dǎo)致網(wǎng)格畸變，不利于數(shù)值計(jì)算收斂，見圖5（a）；另一方面，當(dāng)盾構(gòu)隧道穿越圖5（b）所示的復(fù)雜地層時(shí)，隧道邊界與地層分界面相交并形成若干不規(guī)則幾何邊界，這些不規(guī)則幾何邊界易形成畸形網(wǎng)格.

圖5 復(fù)雜地層盾構(gòu)掘進(jìn)模擬網(wǎng)格劃分難點(diǎn)Fig.5 Difficulties in grid division of complex strata

針對(duì)上述兩方面問題，本文提出一種基于地層分界面的條件判斷算法.該方法從三維地質(zhì)模型提取模型幾何邊界和地表及地層分界面高程數(shù)據(jù)，前者直接輸入網(wǎng)格劃分算法生成未賦予材料屬性的地層網(wǎng)格，后者輸入條件判斷算法以甄別地層網(wǎng)格單元的地層類別并賦予其材料屬性，最終建立復(fù)雜地層的有限元模型.

基于地層分界面的條件判斷算法流程見圖6.具體流程如下：

圖6 復(fù)雜地層網(wǎng)格劃分流程圖Fig.6 Flow chart of complex formation grid division

1）基于Gmsh編寫地層網(wǎng)格生成腳本，實(shí)現(xiàn)地層網(wǎng)格的自動(dòng)化劃分.該腳本采用六面體單元對(duì)地質(zhì)模型幾何邊界確定的建模區(qū)域進(jìn)行均勻網(wǎng)格劃分，生成初始地層網(wǎng)格，且該網(wǎng)格模型與三維地質(zhì)模型采用相同坐標(biāo)系.在初始地層網(wǎng)格中，盾構(gòu)掘進(jìn)區(qū)域設(shè)置為網(wǎng)格協(xié)調(diào)區(qū)，見圖7，得到未賦予材料屬性的地層網(wǎng)格.

圖7 網(wǎng)格協(xié)調(diào)區(qū)示意圖Fig.7 Diagram of grid coordination region

2）讀取上述地層網(wǎng)格數(shù)據(jù)，計(jì)算單元i（i=1，2，…，n；n為單元總數(shù)）的中心點(diǎn)坐標(biāo)（Xi，Yi，Zi）.

3）讀取三維地質(zhì)模型高程數(shù)據(jù)，搜索各地層界面上與單元i中心點(diǎn)水平投影重合的點(diǎn)，并將該點(diǎn)的高程命名為Zlj（j為地層界面序號(hào)；j=1，2，…，N；N為地層界面總數(shù)），其中Zlj＞Zlj+1，見圖8.

圖8 條件判斷示意圖Fig.8 Diagram of conditional judgment

4）對(duì)比各地層界面高程值Zlj與單元i中心點(diǎn)高程值Zi，進(jìn)而確定單元i的地層類別.由地表（Zlj，j=1）開始，判斷Zlj≤Zi是否成立，若成立則判斷單元i位于第j-1 層地層，繼續(xù)第5 步；否則令j=j+1，再次判斷，直到滿足Zlj≤Zi.

5）重復(fù)步驟2）～4），對(duì)地層網(wǎng)格模型中所有單元進(jìn)行逐一判斷，形成各地層的網(wǎng)格單元集.

6）對(duì)各地層單元集賦予相應(yīng)的土體物理力學(xué)參數(shù)，完成復(fù)雜地層的有限元模型構(gòu)建.

通過上述步驟可得到圖9 所示的復(fù)雜地層有限元模型，該模型包含了三維地質(zhì)模型的地質(zhì)信息.

圖9 融合三維地質(zhì)信息的復(fù)雜地層有限元模型Fig.9 Finite element model of complex strata incorporating 3D geological information

3 復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序

為實(shí)現(xiàn)盾構(gòu)掘進(jìn)快速參數(shù)化建模，基于上述方法采用C#語言編寫了復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序.該程序包括3 個(gè)模塊：復(fù)雜地層有限元建模模塊、盾構(gòu)隧道掘進(jìn)建模模塊和計(jì)算求解模塊，見圖10.復(fù)雜地層有限元建模模塊具備三維地質(zhì)建模功能和地層有限元建模功能，可利用鉆孔數(shù)據(jù)文件快速建立融合三維地質(zhì)信息的地層有限元模型.在盾構(gòu)隧道掘進(jìn)建模模塊設(shè)置隧道幾何與掘進(jìn)參數(shù)后，可自動(dòng)生成盾構(gòu)掘進(jìn)模型的inp文件.計(jì)算求解模塊通過批處理命令調(diào)用通用有限元求解器ABAQUS，并讀取盾構(gòu)掘進(jìn)模型的inp文件進(jìn)行計(jì)算，最終實(shí)現(xiàn)復(fù)雜地層條件下盾構(gòu)掘進(jìn)數(shù)值模擬.

圖10 復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序Fig.10 Pre-processing program for simulation of shield tunneling in complex strata

在復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序輸入隧道幾何和掘進(jìn)參數(shù)，即可實(shí)現(xiàn)盾構(gòu)掘進(jìn)的快速參數(shù)化建模.隧道幾何設(shè)置涵蓋隧道埋深、隧道直徑、管片尺寸等參數(shù).隧道掘進(jìn)參數(shù)設(shè)置包括掘進(jìn)步長(zhǎng)、掌子面支護(hù)系數(shù)、管片和土體材料等參數(shù).

基于上述輸入?yún)?shù)，盾構(gòu)隧道掘進(jìn)建模模塊可自動(dòng)完成以下四方面設(shè)置，并生成inp文件：

1）分段開挖設(shè)置.在每個(gè)掘進(jìn)步中，采用“生死”單元法移除掌子面前方土體單元，同時(shí)激活盾構(gòu)管片單元，實(shí)現(xiàn)分段開挖過程的模擬.

2）掌子面支護(hù)力設(shè)置.將掌子面的靜止土壓力作為支護(hù)力默認(rèn)值，沿高度方向呈梯度分布.每個(gè)掘進(jìn)步可根據(jù)實(shí)際施工資料設(shè)置掌子面支護(hù)系數(shù)，將該系數(shù)乘以默認(rèn)支護(hù)力得到實(shí)際支護(hù)力.

3）地層損失設(shè)置.盾構(gòu)掘進(jìn)引起的地層損失是導(dǎo)致地表沉降的主要因素，常以地層損失率［19］表示，并可通過整環(huán)擾動(dòng)層近似模擬［20］.通常盾體與圍巖間的間隙具備不均勻性，即拱頂間隙大于兩旁和底部間隙［21］.為近似考慮該不均勻性，本文將擾動(dòng)層分為上下兩半環(huán)分別設(shè)置材料參數(shù)，見圖11.上半環(huán)擾動(dòng)層參數(shù)取上部圍巖參數(shù)折減值；下半環(huán)擾動(dòng)層采用鄰近圍巖的材料參數(shù).

圖11 地層損失設(shè)置示意圖Fig.11 Diagram of stratum loss setting

4）材料設(shè)置.復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序包括線彈性本構(gòu)模型和摩爾庫倫本構(gòu)模型，前者用于管片和上半環(huán)擾動(dòng)層，后者用于土體和下半環(huán)擾動(dòng)層.

4 工程案例驗(yàn)證

本節(jié)以長(zhǎng)沙地鐵四號(hào)線區(qū)間隧道施工工程為例，利用前處理程序構(gòu)建盾構(gòu)地層掘進(jìn)模型，并開展掘進(jìn)模擬.將沿隧道軸線的地表沉降模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比，探討本文所提方法的有效性.

4.1 工程概況

長(zhǎng)沙地鐵四號(hào)線桐梓坡—望月湖區(qū)間隧道全長(zhǎng)1 512 m，該區(qū)間段采用土壓平衡盾構(gòu)機(jī)掘進(jìn)施工，盾構(gòu)機(jī)外徑為6.28 m.該區(qū)段為雙線隧道.本文選取先期施工的左線隧道掘進(jìn)引起的地表沉降數(shù)據(jù)進(jìn)行驗(yàn)證分析.隧道位置與鉆孔平面分布見圖12.依據(jù)鉆孔樣本數(shù)據(jù)，利用所編寫的前處理程序構(gòu)建該區(qū)域的三維地質(zhì)模型，見圖13.

圖12 隧道與鉆孔分布平面圖Fig.12 Plan view of tunnels and borehole distribution

圖13 三維地質(zhì)模型Fig.13 3D geological model

圖14 為沿隧道縱軸線方向的地質(zhì)剖面圖.該區(qū)段內(nèi)地層復(fù)雜多變，自上而下分別為素填土、硬塑黏土、強(qiáng)風(fēng)化板巖、中風(fēng)化板巖和微風(fēng)化板巖.各地層的物理力學(xué)參數(shù)見表1.

表1 地層物理力學(xué)參數(shù)表Tab.1 Physical and mechanical parameters of soils

圖14 地質(zhì)剖面圖Fig.14 Longitudinal geological profile

4.2 盾構(gòu)掘進(jìn)模擬與驗(yàn)證

為研究地層空間變異性對(duì)地表沉降的影響，采用表2 所示盾構(gòu)隧道幾何和掘進(jìn)參數(shù)，將其輸入復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序建立圖15 所示盾構(gòu)掘進(jìn)有限元模型，并開展盾構(gòu)掘進(jìn)模擬.

表2 盾構(gòu)隧道幾何和掘進(jìn)參數(shù)表Tab.2 Geometry and excavation parameter of shield tunnel

圖15 盾構(gòu)掘進(jìn)模擬有限元模型Fig.15 Shield tunneling simulation finite element model

為驗(yàn)證有限元模型的準(zhǔn)確性，選取掘進(jìn)斷面內(nèi)極軟巖占比為50%與100%的兩種工況，分別對(duì)應(yīng)711 環(huán)和730 環(huán)地表沉降監(jiān)測(cè)點(diǎn)，并對(duì)比兩種工況下地表沉降變化過程的模擬結(jié)果與實(shí)測(cè)結(jié)果，見圖16.根據(jù)盾構(gòu)施工工法，盾構(gòu)掘進(jìn)對(duì)監(jiān)控?cái)嗝娴乇沓两档挠绊懛譃閮蓚€(gè)階段，盾構(gòu)掘進(jìn)抵達(dá)監(jiān)控?cái)嗝媲胺Q為第一階段，盾構(gòu)通過監(jiān)控?cái)嗝婧蠓Q為第二階段.在地表沉降規(guī)律方面，兩種工況下，兩階段的模擬結(jié)果與實(shí)測(cè)結(jié)果趨勢(shì)均一致.在第一階段期間，地表沉降均隨掌子面與監(jiān)控?cái)嗝娴木嚯x減小而逐漸增大，而在第二階段期間，地表沉降增長(zhǎng)變緩且最終趨于穩(wěn)定.在沉降幅值方面，兩種工況下地表產(chǎn)生了超過10 mm 的差異沉降量，具體為：當(dāng)盾構(gòu)機(jī)由掘進(jìn)斷面部分極軟巖地層進(jìn)入全斷面極軟巖地層時(shí)，地表最終沉降量模擬與實(shí)測(cè)分別增加了11.67 mm 和18.02 mm，模擬與實(shí)測(cè)的誤差為35%.誤差主要來源于730環(huán)，根據(jù)施工資料，盾構(gòu)機(jī)開挖至730環(huán)停機(jī)約10 h，而有限元模型無法考慮停機(jī)時(shí)未及時(shí)注漿引起的地層損失，因此730 環(huán)最終沉降模擬值小于實(shí)測(cè)值.綜上，有限元模擬結(jié)果與實(shí)測(cè)結(jié)果趨勢(shì)一致，吻合程度好，表明有限元模型計(jì)算準(zhǔn)確.

圖16 監(jiān)測(cè)點(diǎn)地表沉降變化過程模擬與實(shí)測(cè)結(jié)果對(duì)比Fig.16 Comparison of simulation and measured results of surface subsidence change process at monitoring points

為研究地層空間變異性對(duì)地表沉降的影響，將各監(jiān)測(cè)點(diǎn)最終沉降的模擬結(jié)果和實(shí)測(cè)數(shù)據(jù)繪制于圖17 中.模擬結(jié)果顯示，地表最終沉降隨掘進(jìn)斷面內(nèi)極軟巖占比增大而顯著增加，與實(shí)測(cè)數(shù)據(jù)規(guī)律一致.綜合圖16～17 結(jié)果，掌子面接近監(jiān)測(cè)斷面前，開挖卸載作用引起掌子面前方極軟巖地層發(fā)生變形，導(dǎo)致地表發(fā)生第一階段沉降.掌子面通過監(jiān)測(cè)斷面后，極軟巖地層向隧道斷面內(nèi)發(fā)生位移，導(dǎo)致第二階段地表沉降.當(dāng)極軟巖占比增加時(shí)，上述兩階段的地表沉降均增大，進(jìn)而導(dǎo)致地表最終沉降顯著增加.

圖17 沿隧道軸線地表最終沉降模擬與實(shí)測(cè)結(jié)果對(duì)比圖Fig.17 Comparison of simulation and measured results of final settlement along tunnel axis

綜上，復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序可實(shí)現(xiàn)盾構(gòu)掘進(jìn)快速參數(shù)化建模，所建有限元模型能較準(zhǔn)確地反映復(fù)雜地層條件下盾構(gòu)掘進(jìn)過程，并能有效揭示地層空間變異性對(duì)地表沉降的影響規(guī)律.

5 結(jié) 論

本文提出了考慮地層三維空間分布的盾構(gòu)掘進(jìn)模擬方法，編寫了相應(yīng)的復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序，并以長(zhǎng)沙地鐵四號(hào)線區(qū)間隧道工程為例，驗(yàn)證了所提方法和相應(yīng)程序的有效性.主要結(jié)論如下：

1）基于地層分界面的條件判斷算法可以規(guī)避復(fù)雜地層曲面的有限元網(wǎng)格畸變問題，并且有效解決了復(fù)雜地層曲面與隧道斷面的網(wǎng)格協(xié)調(diào)問題.

2）復(fù)雜地層盾構(gòu)掘進(jìn)模擬前處理程序集成了三維地質(zhì)建模功能和復(fù)雜地層網(wǎng)格劃分算法，實(shí)現(xiàn)了復(fù)雜地層掘進(jìn)模擬的參數(shù)化建模，大幅提高了建模效率.

3）盾構(gòu)掘進(jìn)模擬結(jié)果表明，地表沉降量隨盾構(gòu)掘進(jìn)斷面內(nèi)軟巖占比增大而顯著增加，與工程實(shí)測(cè)數(shù)據(jù)規(guī)律一致，驗(yàn)證了本文所提方法的有效性.