徐 兵 伍文文

(1.西南科技大學(xué)環(huán)境與資源學(xué)院,四川 綿陽 621010; 2.武漢理工大學(xué)土木工程與建筑學(xué)院,湖北 武漢 430070)

基于隨機(jī)介質(zhì)理論的盾構(gòu)隧道地表沉降研究★

徐 兵1伍文文2

(1.西南科技大學(xué)環(huán)境與資源學(xué)院,四川 綿陽 621010; 2.武漢理工大學(xué)土木工程與建筑學(xué)院,湖北 武漢 430070)

介紹了隨機(jī)介質(zhì)理論的方法原理及反分析數(shù)值算法，并使用MATLAB編制了反分析程序，結(jié)合具體實(shí)例使用編制的程序?qū)?shù)據(jù)進(jìn)行了反分析參數(shù)統(tǒng)計(jì)，給出了該工程的監(jiān)測范圍建議值，并指出了地層損失率較高的若干斷面，并建議及時(shí)進(jìn)行二次注漿保證安全。

地表沉降，隨機(jī)介質(zhì)理論，收斂模式，地層損失率

0 引言

地下空間開挖常常引起地表的塌陷災(zāi)害，這種情況在隧道開挖時(shí)更加普遍。由于隧道的開挖常常導(dǎo)致上方地表沉降塌陷，嚴(yán)重時(shí)導(dǎo)致建筑物開裂，造成大量的財(cái)產(chǎn)損失，因此研究地下隧道開挖引起的地表沉降規(guī)律就顯得尤為重要。常用的方法有經(jīng)驗(yàn)公式法，其中有代表性的為Peck公式[1]，其主要思想認(rèn)為地下隧道開挖引起的地表沉降可近似看做一條正態(tài)分布曲線，具體見式(1)。

(1)

其中，Smax為開挖斷面最大沉降值；x為沉降點(diǎn)到斷面中心的距離；i為沉降槽寬度系數(shù)。很多學(xué)者均使用該方法對各地隧道開挖引起的地表沉降進(jìn)行了研究和統(tǒng)計(jì)[2,3]。該方法的主要參數(shù)為Smax和i，并沒有考慮到地層參數(shù)的影響，因此雖然應(yīng)用簡單但有一定的局限性。本文參考了隨機(jī)介質(zhì)理論對某地下隧道開挖引起的地表沉降進(jìn)行了研究，給出了該隧道開挖引起地表沉降的相關(guān)參數(shù)，對地表沉降的預(yù)測和控制提供了建議和參考，該方法充分考慮了地層主要影響角的影響，相比Peck公式法更為科學(xué)。

1 隨機(jī)介質(zhì)理論原理

隨機(jī)介質(zhì)理論[4]最早由波蘭籍學(xué)者Litwinszyn提出，用于研究煤礦開挖引起的地表位移問題。后來由陽軍生、劉寶琛[5,6]引入國內(nèi)并深入研究發(fā)展，被普遍用于地下隧道開挖引起的土體變位問題。另外，該方法還可以廣泛應(yīng)用于地下降水、隧道冷凍法施工等過程造成的地表位移研究。隨機(jī)介質(zhì)理論的具體原理見圖1，該理論認(rèn)為地下隧道開挖引起的地表沉降為隧道斷面收縮區(qū)域的各個(gè)微元塌陷效果的疊加，圖1中Ω為巷道擾動前的原始邊界，ω為巷道穩(wěn)定后的收斂邊界，圖1中dξdη為兩邊界之間收斂區(qū)域內(nèi)的微元體，地表沉降可以分解為兩個(gè)邊界之間收斂區(qū)域中的微元體引起沉降所疊加的結(jié)果。β為地層的主要影響角，微元沿其地層主要影響角β擴(kuò)展到地表的點(diǎn)就是該微元塌陷所造成的地表沉降邊界，對應(yīng)的影響半徑為圖1中的r，最終的整體沉降為各個(gè)微元在收斂范圍內(nèi)沉降的積分，具體見式(2)。

(2)

2 斷面收斂模式

由于隨機(jī)介質(zhì)理論需要考慮隧道斷面的收斂范圍，因此對于不同的收斂形態(tài)所反分析得到的參數(shù)也不一樣。常用的收斂模式分為兩種：均勻收斂和不均勻收斂，以圓形斷面隧道為例，其收斂模式分別見圖2。

從圖2中可以看出，均勻收斂模式表現(xiàn)為斷面均勻向中心收縮ΔR距離，而不均勻收斂則表現(xiàn)為斷面頂部向下收縮2ΔR距離，而底部則不收縮，即相當(dāng)于在均勻收斂的基礎(chǔ)上向下移動了2ΔR距離。兩種收斂模式中收斂范圍空間位置不同，反分析結(jié)果會有所差別，而一般認(rèn)為不均勻收斂更符合實(shí)際情況[7]，因此本文考慮不均勻收斂模式。

3 隨機(jī)介質(zhì)理論參數(shù)反分析

由式(2)可知根據(jù)收斂范圍以及地層主要影響角β即可求得最終的沉降曲線，但是在實(shí)際應(yīng)用中的過程恰好相反，需要應(yīng)用實(shí)際測量得到的斷面沉降數(shù)據(jù)根據(jù)式(2)過程來反分析出參數(shù)，以便得到該工程中的經(jīng)驗(yàn)參數(shù)來指導(dǎo)后續(xù)施工過程中的沉降預(yù)測和控制。由于式(2)的被積函數(shù)很復(fù)雜，只能求得數(shù)值解，可采用Gauss-Legendre積分實(shí)現(xiàn)。

(3)

4 應(yīng)用實(shí)例

表1 各斷面埋深及反分析結(jié)果

武漢地鐵四號線羅家港—園林路區(qū)間隧道采用直徑6m的單孔圓形盾構(gòu)機(jī)進(jìn)行開挖，該地區(qū)屬于長江一級階地，自上而下分別為人工填土、軟粘土、砂土和泥質(zhì)砂巖，盾構(gòu)穿越的土層主要為粉細(xì)砂土。施工過程中發(fā)現(xiàn)盾構(gòu)穿越地表沉降較為明顯，表現(xiàn)為地表明顯下沉，影響區(qū)域內(nèi)部分居民住房外墻有裂縫。為保障施工安全，擬對該地質(zhì)條件下盾構(gòu)開挖引起的地表沉降規(guī)律進(jìn)行研究，數(shù)據(jù)選擇過程中排除了數(shù)據(jù)不完整的監(jiān)測斷面，共選取了16個(gè)完整斷面的沉降監(jiān)測數(shù)據(jù)，每個(gè)斷面有7個(gè)沉降監(jiān)測點(diǎn)。隨機(jī)介質(zhì)理論的收斂模式選取不均勻收斂，所有16個(gè)斷面埋深H及參數(shù)反分析結(jié)果見表1。因篇幅限制，僅以前4個(gè)斷面為例，按照前述反分析方法得出相應(yīng)的沉降斷面反分析擬合曲線，見圖3。

從圖3可以看出隨機(jī)介質(zhì)理論沉降曲線可很好地描述隧道開挖引起的地表沉降規(guī)律，擬合效果好。根據(jù)表1的參數(shù)反分析統(tǒng)計(jì)結(jié)果，地層主要影響角β的平均值為44.75°，收斂參數(shù)ΔR的平均值為29.6 mm。而根據(jù)隨機(jī)介質(zhì)理論原理，隧道斷面收斂邊界擴(kuò)展到地表的最大影響半徑可廣義表示為：

(4)

而對于圓形隧道，部分學(xué)者則給出了考慮斷面具體形狀的更為準(zhǔn)確的圓形隧道沉降影響半徑公式[9]：

r=Rsinβ+(H+Rcosβ)cotβ

(5)

由于已經(jīng)得到該地區(qū)盾構(gòu)隧道開挖引起的地表沉降的地層主要影響角β的平均值，則在實(shí)際工程中可以參考該平均值結(jié)合隧道的埋深進(jìn)行地表沉降范圍的評估，便于沉降監(jiān)測點(diǎn)的布置以及觀測。為保守起見，統(tǒng)一選取區(qū)間較大的埋深20 m來計(jì)算沉降半徑。根據(jù)式(5)可算得該隧道開挖地表沉降半徑為16.7 m，該半徑可作為工程中監(jiān)測點(diǎn)布置范圍，可有效指導(dǎo)安全監(jiān)測的進(jìn)行。

另外，根據(jù)收斂參數(shù)ΔR可以得到該隧道各個(gè)斷面的地層損失率，其計(jì)算方法見式(6)，該損失率可以作為二次注漿的參考標(biāo)準(zhǔn)指導(dǎo)盾構(gòu)的二次注漿來控制地表沉降。對于一級階地地區(qū)，建議地層損失率應(yīng)控制在2%以下，否則會引起較大的危害。按照式(6)得到的各個(gè)斷面的地層損失率見表2。從表2中可知，斷面3，6,11～15的地層損失率均超過2%，應(yīng)及時(shí)進(jìn)行二次注漿控制地表沉降，確保施工安全。

η (6)

5 結(jié)語

通過MATLAB編制了隨機(jī)介質(zhì)理論的參數(shù)反分析程序，對武漢地鐵某區(qū)間的16個(gè)斷面沉降數(shù)據(jù)進(jìn)行了反分析和統(tǒng)計(jì)，結(jié)果表明該工程地區(qū)的地層主要影響角平均值為44.75°，收斂參數(shù)ΔR的平均值為29.6 mm。為保守起見，根據(jù)公式建議該隧道地表沉降監(jiān)測的半徑范圍為16.7 m。根據(jù)反分析得到的收斂參數(shù)ΔR計(jì)算得到各個(gè)斷面的地層損失率，并與建議的地層損失率2%的上限進(jìn)行了對比，指出了需要進(jìn)行及時(shí)二次注漿的斷面，為隧道工程的安全提供了保障和參考。

[1] PECK R B.Deep excavations and tunnelling in soft ground[C]// Proceedings of the 7th International Conference on Soil Mechanics and Foundation Engineering. Mexico：[s.n.].1969：225-290.

[2] 陳春舒，夏元友.Peck公式在長江一級階地盾構(gòu)中的適用性研究[J].武漢理工大學(xué)學(xué)報(bào),2013,35(9):85-90.

[3] 胡 斌,莫 云,胡新麗,等.Peck法在武漢地鐵隧道地表沉降預(yù)測中的適用性分析[J].工程勘察,2012,40(7):6-10.

[4] Litwiniszyn J.Fundamental principles of the mechanics of stochastic medium[C]//Proceedings of the 3rd Conference on Theoretical Applied Mechanics.Bangalore,India:[sn].1957:18-26.

[5] 陽軍生，劉寶琛.擠壓式盾構(gòu)隧道施工引起的地表移動及變形[J].巖土力學(xué),1998,19(3):10-13.

[6] 陽軍生，劉寶琛.城市隧道施工引起的地表移動及變形[M].北京：中國鐵道出版社，2002.

[7] 韓 煊,李 寧.隧道開挖不均勻收斂引起地層位移的預(yù)測模型[J].巖土工程學(xué)報(bào),2007,29(3):347-352.

[8] 鞏敦衛(wèi)，潘鳳萍.自適應(yīng)遺傳算法理論及應(yīng)用[M].北京：中國礦業(yè)大學(xué)出版社,2003.

[9] 陳春舒.地層位移模型在淺埋隧道工程中的適用性研究[D].武漢:武漢理工大學(xué),2013.

Research on surface settlement induced by shield tunnel based on stochastic medium theory★

Xu Bing1Wu Wenwen2

(1.SchoolofEnvironmentandResource,SouthwestUniversityofScienceandTechnology,Mianyang621010,China; 2.SchoolofCivilEngineeringandArchitecture,WuhanUniversityofTechnology,Wuhan430070,China)

The stochastic medium theory and it’s numerical algorithms of inverse analysis are introduced. The inverse analysis program is also developed based on MATLAB. Combined with data measured in a case, parameter inverse analysis are carried out with the developed program. The surface monitoring radius is recommended. And some sections are alarmed because the ground loss ratios are all above the recommended standard and secondary grouting are recommended to ensure safety.

surface settlement, stochastic medium theory, convergence form, ground loss ratio

1009-6825(2015)07-0168-03

2014-12-30 ★：四川省教育廳科研項(xiàng)目(項(xiàng)目編號：15ZB0124)；綿陽市科技計(jì)劃項(xiàng)目(項(xiàng)目編號：14S-02-6)

徐 兵(1976- )，男，講師； 伍文文(1990- )，女，在讀碩士

U455.43

A