蔡海兵,黃以春,李 陽

( 安徽理工大學(xué) 土木建筑學(xué)院,安徽 淮南 232001)

基于隨機(jī)介質(zhì)理論的土體融沉預(yù)測及其參數(shù)敏感性分析

蔡海兵,黃以春,李 陽

( 安徽理工大學(xué) 土木建筑學(xué)院,安徽 淮南 232001)

為預(yù)測人工凍結(jié)法施工引起的地表融沉,以廣州某地鐵為研究背景,考慮土體壓力的變化,運(yùn)用隨機(jī)介質(zhì)理論建立馬蹄形水平凍結(jié)融沉計(jì)算模型,分析對(duì)地表融沉有影響的各參數(shù)的敏感性。由計(jì)算結(jié)果得出,水平凍結(jié)引起的地表豎向融沉呈正態(tài)分布,在隧道中心處達(dá)到最大值；參數(shù)敏感性分析得出各參數(shù)敏感度從大到小依次為凍結(jié)壁厚度、土體主要影響角正切值、隧道埋深、融化與壓密系數(shù)；最大豎向融沉隨著凍結(jié)壁厚度、土體主要影響角正切值、融化與壓縮系數(shù)的增大而增大,隨著隧道埋深的增大而減小。

地鐵；人工凍結(jié)法；土壓力；融沉；隨機(jī)介質(zhì)理論；參數(shù)分析

1 概述

人工凍結(jié)輔助施工法具有安全、可靠、無污染等特點(diǎn)，被廣泛運(yùn)用到礦井建設(shè)、隧道修筑、深基坑開挖等領(lǐng)域[1-2]，起到土體的臨時(shí)加固及隔絕地下水等作用。

該工法具有優(yōu)越性的同時(shí)也具有缺陷：在凍結(jié)時(shí)期水的體積會(huì)膨脹9%[3]，而“分凝作用”又會(huì)導(dǎo)致水分的遷移，從而極大地加劇這種膨脹；在融化時(shí)期冰的體積會(huì)減小，且土體在上部荷載的作用下發(fā)生壓縮沉降，這種效應(yīng)稱為“凍脹”和“融沉”。一般情況下土體的融沉要比凍脹值大許多，這會(huì)對(duì)周圍的建筑物、地下管線等構(gòu)造物的安全造成危害[4-5]。因此，對(duì)人工凍結(jié)引起的土體融沉進(jìn)行準(zhǔn)確的預(yù)測就顯得尤為重要。我國的專家學(xué)者在這方面研究取得了許多的成果，如程樺[6]、李大勇[7]等采用有限元對(duì)人工凍結(jié)工程凍融效應(yīng)進(jìn)行了數(shù)值模擬分析。但由于凍土是土、冰、水、氣四相體，物理力學(xué)性質(zhì)復(fù)雜，采用數(shù)值模擬時(shí)涉及到的參數(shù)眾多且不容易得到確定，而采用隨機(jī)介質(zhì)理論進(jìn)行分析計(jì)算能得到合理的結(jié)果且實(shí)施起來簡單，更適合在實(shí)際工程中應(yīng)用。如蔡海兵[8]等基于該理論建立了隧道水平凍結(jié)施工期的融沉歷時(shí)預(yù)測模型，研究了隧道水平凍結(jié)引起的地表融沉量大小及其隨時(shí)間的變化規(guī)律。劉波[9]等研究了土體壓力對(duì)地鐵解凍后地表沉降的影響，得出了融化沉降與壓縮沉降的計(jì)算方法及土體壓力對(duì)壓縮沉降的影響結(jié)果。

上述研究都得出了一些有意義的結(jié)論與成果，但不足之處是均將土體壓力視為恒值，未考慮隧道周圍的土體壓力在豎直方向上的變化，在建立模型時(shí)將隧道形狀簡化為圓形或矩形，這與實(shí)際情況不符，特別是對(duì)于大斷面不規(guī)則的隧道將很大程度地影響計(jì)算精度。本文考慮土體壓力隨深度的變化及隧道的實(shí)際形狀，采用隨機(jī)介質(zhì)理論和一維沉降公式建立了馬蹄形隧道水平凍結(jié)的地表沉降預(yù)測模型；分析了融化參數(shù)對(duì)地表沉降量的影響情況，得出各參數(shù)與地表沉降量之間的關(guān)系。

2 計(jì)算模型

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

20世紀(jì)50年代，波蘭學(xué)者Litwiniszyn[10]提出隨機(jī)介質(zhì)理論。該理論將土體視為一種隨機(jī)介質(zhì)，假設(shè)隧道、煤巷等在開挖時(shí)土體不產(chǎn)生固結(jié)、壓縮，則造成地層沉降的原因是由于開挖空間在周圍水土壓力作用下“收斂”造成的。而整個(gè)“收斂區(qū)域”可以看成由無數(shù)個(gè)“收斂單元”組成，所有“收斂單元”造成的沉降之和等于整個(gè)“收斂區(qū)域”的沉降。假設(shè)收斂單元為dεdηdξ，由此單元造成的沉降稱為“單元沉降”，對(duì)其在收斂區(qū)域內(nèi)積分即可得到整個(gè)“收斂區(qū)域”造成的沉降。如圖1所示，由隨機(jī)介質(zhì)理論得出因該單元“收斂”造成對(duì)應(yīng)的地表豎向沉降

(1)

對(duì)于融沉問題，地表的沉降形狀和結(jié)果與土體開挖相似，即公式(1)亦適用于地表融沉的計(jì)算。一般將水平凍結(jié)工程解凍時(shí)引起的地表沉降視為平面二維問題。將上述公式在y軸上從-∞到+∞積分得到

(2)

式中，β為地層的主要影響角，與土體的物理力學(xué)性質(zhì)有關(guān)。

圖1 單元開挖示意

2.2 一維沉降理論

前蘇聯(lián)學(xué)者崔托維奇[12]在多年的凍土試驗(yàn)基礎(chǔ)上提出了著名的凍土融化一維沉降理論。他認(rèn)為凍土融化時(shí)會(huì)產(chǎn)生兩個(gè)變化：一是由于凍土的膠凝材料“冰”在融化時(shí)體積會(huì)減小10%，二是融土本身蓬松多孔，在上部土體壓力的作用下發(fā)生壓縮固結(jié)。因此凍土融化時(shí)產(chǎn)生的沉降由兩部分組成，在一維條件下土體的沉降量公式為

(3)

2.3 土體融沉模型

當(dāng)凍結(jié)壁發(fā)展到預(yù)定厚度時(shí)凍結(jié)工程進(jìn)入維護(hù)凍結(jié)期，便可在凍結(jié)壁的保護(hù)下施工。工程結(jié)束后便停止凍結(jié)，將凍結(jié)壁解凍。實(shí)際工程中一般分為自然解凍和強(qiáng)制解凍，這兩種方法只是溫度場分布與融化速度不同，最終的融化沉降量一致。此時(shí)隧道的襯砌已經(jīng)完成，在四周水土壓力的作用下已融土體被壓縮，此時(shí)會(huì)產(chǎn)生“土體損失”，反映在地層上即土體的沉降。這種沉降是隨機(jī)的，和土體開挖產(chǎn)生沉降的形狀與結(jié)果相同，因此可以采用隨機(jī)介質(zhì)理論的方法進(jìn)行計(jì)算。由于隧道襯砌的彈性模量遠(yuǎn)遠(yuǎn)大于土體，可將其視為不變形體，四周壓力使融土向隧道中心方向收縮，且深度不同土體收縮量亦不相同。即凍結(jié)壁融化收縮的區(qū)域可認(rèn)為是從凍結(jié)壁外邊緣R到融后邊界Rs，如圖2所示，在凍結(jié)壁上取一微元有

(4)

(5)

(6)

式中，μ為土體側(cè)向壓力系數(shù)；γ為土體平均重度；h為隧道軸線埋深。

由一維沉降公式得出考慮土體壓力隨深度變化情況下的凍土融沉量為

(7)

由隨機(jī)介質(zhì)理論得凍結(jié)壁融化引起的地表位移為

(8)

其中，Rs=R-ΔS

圖2 凍結(jié)壁受力分析

3 實(shí)例計(jì)算與分析

圖3 隧道凍結(jié)設(shè)計(jì)斷面(單位：mm)

在建立計(jì)算模型時(shí)考慮到凍結(jié)管為馬蹄形布置，若將其簡單等效為圓形會(huì)造成與實(shí)際情況的偏差，影響計(jì)算精度。故將凍結(jié)壁分為4個(gè)不同半徑的異心圓弧分段計(jì)算，保證了模型的準(zhǔn)確性。各段凍結(jié)壁的參數(shù)見表1。

表1 凍結(jié)壁參數(shù)

由一維融化沉降公式得到每部分凍結(jié)壁在不同深度土壓力作用下的沉降量為

(9)

將表1中各參數(shù)代入公式得

(10)

土體收斂后的邊界為

(11)

所以各部分凍結(jié)壁的積分區(qū)域?yàn)?/p>

采用線性疊加原理得到總體融沉量為

(12)

由于公式中存在不可積函數(shù)，所以采用Maple數(shù)學(xué)計(jì)算軟件編寫凍脹程序進(jìn)行計(jì)算，結(jié)果見圖4。

圖4 地表豎向融沉分布規(guī)律

從計(jì)算結(jié)果可以看出，水平凍結(jié)工程解凍后引起的地表豎向融沉以隧道軸線為對(duì)稱軸呈正態(tài)分布。其中在隧道中心處沉降最大，隨著距離軸線距離增加融沉減小，超過一定范圍后融沉很小可不計(jì)。最大融沉值達(dá)到110.86 mm。此隧道埋深淺、斷面大，地表對(duì)沉降非常敏感，實(shí)際施工時(shí)應(yīng)在解凍時(shí)采取跟蹤注漿等措施減小地表沉降。雖然此計(jì)算是在不考慮注漿的情況下得出的，但是用于沉降預(yù)測及對(duì)相應(yīng)的預(yù)防措施進(jìn)行論證是十分有用的。此預(yù)測模型也有待改進(jìn)的地方：實(shí)際工程地質(zhì)條件復(fù)雜，地層往往是由多層不同的土層組成，各土層的物理力學(xué)性質(zhì)差異較大，此模型沒有對(duì)土體進(jìn)行分層考慮而是將其視為整體；該計(jì)算結(jié)果表示最終沉降量大小，實(shí)際沉降則是時(shí)間的函數(shù)，對(duì)于砂性土在解凍后短時(shí)間內(nèi)便可達(dá)到最終沉降量。黏性土的沉降具有滯后性，沉降量一直在變化，需要幾個(gè)月甚至幾年才能逐漸穩(wěn)定。此模型更適合用于砂性土體的融化沉降預(yù)測。

4 參數(shù)敏感性分析

參數(shù)敏感性分析是通過改變相關(guān)參數(shù)的數(shù)值來分析系統(tǒng)模型受這些參數(shù)變動(dòng)影響大小的一種方法。假設(shè)系統(tǒng)模型為y=f(x1,x2,…,xn)，x為系統(tǒng)模型的參數(shù)。當(dāng)x偏離基準(zhǔn)參數(shù)x*時(shí)產(chǎn)生偏移值Δx，同理y的偏移值為Δy，則參數(shù)的敏感性系數(shù)表示為k=δy/δx，δx和δy分別為參數(shù)x與系統(tǒng)模型輸出值y的變化率，δx=Δx/x*，δy=Δy/y*。各參數(shù)敏感性系數(shù)的大小代表了該參數(shù)對(duì)系統(tǒng)模型的重要程度。

表2 基準(zhǔn)參數(shù)

各基準(zhǔn)參數(shù)合理的上下浮動(dòng)范圍為[-30%，-20%，-10%，0，10%，20%，30%]，為了計(jì)算方便，假設(shè)凍結(jié)壁為圓形且土體壓力為定值，取凍結(jié)壁外半徑為5 m。將上述各組參數(shù)帶入地表豎向融沉計(jì)算公式中，得到參數(shù)變化下地表豎向融沉，如圖5～圖9所示。

圖5 主要影響角變化下地表融沉分布規(guī)律

圖6 隧道埋深變化下地表融沉分布規(guī)律

圖7 凍結(jié)壁厚度變化下地表融沉分布規(guī)律

圖8 壓縮系數(shù)變化下地表融沉分布規(guī)律

圖9 融化系數(shù)變化下地表融沉分布規(guī)律

從圖5～圖9中得知，隨著凍結(jié)壁厚度、融化系數(shù)、壓縮系數(shù)的增大，地表各點(diǎn)的豎向融沉都在增大，且影響范圍也隨之增大；當(dāng)土體主要影響角正切值增大時(shí)，在隧道軸線附近處(-8 m

從圖10可知，參數(shù)變化率與隧道軸線處地表融沉變化率基本呈線性關(guān)系，直線斜率的絕對(duì)值大小代表了該參數(shù)的敏感性系數(shù)。參數(shù)的敏感性系數(shù)從大到小依次為凍結(jié)壁厚度>主要影響角正切值>隧道埋深>融化系數(shù)=壓縮系數(shù)。通過參數(shù)變化率與地表豎向融沉關(guān)系可知，進(jìn)行人工凍結(jié)設(shè)計(jì)時(shí)在保證安全與施工要求的前提下盡可能減小凍結(jié)壁厚度，合理地增大隧道埋深可以減少融沉效應(yīng)；在施工過程中通過跟蹤注漿甚至改良土體的方法亦可減小融沉。

圖10 參數(shù)變化率與地表融沉關(guān)系曲線

5 結(jié)論

(1)廣州某地鐵水平凍結(jié)法施工引起的地表豎向融沉曲線近似正態(tài)分布，在隧道中心處的地表融沉值最大，達(dá)到110.86 mm。隨著距中心距離的增大融沉逐漸減小，當(dāng)超出一定影響范圍后值接近于零，可忽略不計(jì)。

(2)各參數(shù)對(duì)地表融沉的敏感度從大到小依次為凍結(jié)壁厚度、土體主要影響角正切值、隧道埋深、土體融化與壓縮系數(shù)。其中，當(dāng)凍結(jié)壁厚度、主要影響角正切值、土體融化與壓縮系數(shù)增大時(shí)，隧道最大豎向融沉也隨之增大，當(dāng)隧道埋深增大時(shí)最大豎向融沉減小。

(3)實(shí)際工程中需要控制地表的最大融沉量，通過減小凍結(jié)壁厚度，適當(dāng)增加隧道埋深，進(jìn)行跟蹤注漿等方法可以有效地減小土體的融沉。

[1]武冰冰，蔡海兵，程樺.地鐵隧道盾構(gòu)洞門的垂直凍結(jié)加固技術(shù)[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì)，2011(9):80-84.

[2]周曉敏，土夢恕，張頂立，崔海濤.地層凍結(jié)技術(shù)在北京地鐵施工中的應(yīng)用分析[J].巖土工程界，2001，5(3):61-64.

[3]馬芹永.人工凍結(jié)法的理論與施工技術(shù)[M].北京:人民交通出版社，2007.

[4]李方正.人工地層凍結(jié)的環(huán)境效應(yīng)及其工程對(duì)策[J].公路交通科技，2004，24(1):67-70.

[5]陶敬德，王明年，劉大剛.凍結(jié)法隧道施工引起的地表移動(dòng)及變形預(yù)測[J].現(xiàn)代隧道技術(shù)，2006，43(6):45-50.

[6]程樺，臧華.人工地層水平凍結(jié)凍脹效應(yīng)準(zhǔn)耦合數(shù)值分析[J].巖土工程學(xué)報(bào)，2003，25(1):87-90.

[7]李大勇，趙少飛，胡向東.越江隧道泵房凍結(jié)施工三維數(shù)值模擬[J].巖土力學(xué)，2004(9):67-70.

[8]蔡海兵，彭立敏，鄭騰龍.隧道水平凍結(jié)施工期地表融沉的歷時(shí)預(yù)測模型[J].巖土力學(xué)，2014，35(2):504-509.

[9]劉波，李東陽，廖建軍.土體壓力對(duì)地鐵旁通道解凍后地表沉降的影響[J].煤炭學(xué)報(bào)，2011，36(4):551-555.

[10]Litwinisyn J. Fundamental principles of the mechanics of stochastic medium[C]∥Proceedings of the 3rd Conference on Theoretical Applied Mechanics， Bangalore， India， 1957:18-26.

[11]劉寶琛，張家生.近地表開挖引起的地表沉降的隨機(jī)介質(zhì)方法[J].巖石力學(xué)與工程學(xué)報(bào)，1995，14(4):289-295.

[12]崔托維奇H A.凍土力學(xué)[M].北京:科學(xué)出版社，1985.

Prediction of Surface Thawing Settlement Based on Stochastic Medium Theory and Parameter Sensitivity Analysis

CAI Hai-bing, HUANG Yi-chun, LI Yang

(School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan Anhui 232001, China)

To predict surface thawing settlement caused by artificial freezing, the horseshoe horizontal freezing calculation model on a subway in Guangzhou is established. Stochastic medium theory is employed to calculate the thawing settlement in view of the change of earth pressure and the analysis of the sensitivity of parameters that have influence on the surface frost heave is conducted. The calculation results show that the surface vertical thawing settlement caused by horizontal freezing tends to be in normal distribution, and the largest values are in the center of the tunnel. The sensitivity of the parameters from high to low are the frozen wall thickness, the tangent value of main influencing angle, the buried depth of tunnel, the coefficient of thaw, the coefficient of thaw compression. The maximum vertical displacement increases with the increasing of frozen wall thickness, tangent value of the main influencing angle, coefficient of thaw and coefficient of thaw compression, but decreases with the increasing of the depth of tunnel.

Subway; Artificial freezing; Soil pressure; Surface thawing settlement; Stochastic medium theory; Parameter sensitivity

2014-09-13；

2014-09-27

國家自然科學(xué)基金項(xiàng)目(51208004)

蔡海兵(1980—)，男，副教授，博士研究生，E-mail：haibingcai@163.com。

1004-2954(2015)08-0107-05

U455.49

A

10.13238/j.issn.1004-2954.2015.08.023