胡楠

【摘 要】本文結(jié)合某隧道穿越山體沖溝的松散黃土層出現(xiàn)的塌方事故，在現(xiàn)場(chǎng)地質(zhì)調(diào)研和監(jiān)控測(cè)量的基礎(chǔ)上，應(yīng)用有限差分軟件FLAC3D，根據(jù)對(duì)現(xiàn)場(chǎng)地形和施工工序步驟的數(shù)值模擬，計(jì)算分析隧道周邊圍巖和土層隨著隧道開(kāi)挖過(guò)程受力狀態(tài)的變化。結(jié)果顯示：通過(guò)正臺(tái)階法方式開(kāi)挖進(jìn)入淺埋區(qū)域時(shí)，圍巖變形發(fā)展迅速、拱頂沉降和地表沉降值過(guò)大，遠(yuǎn)超過(guò)正常變形范圍，且存在潛在的剪切帶。因此可認(rèn)為，隧道在進(jìn)入淺埋區(qū)域后的開(kāi)挖會(huì)導(dǎo)致圍巖的破壞和地表土體的沉陷。分析結(jié)果與現(xiàn)場(chǎng)測(cè)量數(shù)據(jù)基本吻合，為后期加固施工提供了理論依據(jù)。

【關(guān)鍵詞】松散黃土層；數(shù)值模擬；受力狀態(tài)；圍巖變形；地表土沉降

0 引言

隨著隧道建設(shè)的快速發(fā)展，越來(lái)越多的隧道線路需要穿越淺埋黃土區(qū)。受埋深淺、黃土承載力有限的限制，淺埋黃土隧道施工時(shí)容易出現(xiàn)支護(hù)變形過(guò)大、支護(hù)開(kāi)裂甚至冒頂塌方的情況[1-2]。20世紀(jì)70年代以來(lái)，有限元法已經(jīng)廣泛地應(yīng)用于土體的動(dòng)力分析中。近年來(lái)，拉格朗日元由于能解決大變形問(wèn)題而倍受青睞，美國(guó)Itasca公司推出的FLAC3D能夠很好地進(jìn)行動(dòng)力分析[3]。

FLAC3D是連續(xù)介質(zhì)快速拉格朗日差分分析方法（Fast Lagrangian Analysis of Continuum）的英文縮寫(xiě)。美國(guó)Itasca Consulting Group Inc. 將此方法用于巖土體的工程力學(xué)計(jì)算中，并于1986年開(kāi)發(fā)出Flac[4]。目前，該軟件從二維平面分析拓展到三維空間分析，已成為處理功能強(qiáng)大的新一代軟件—Flac3D。

Flac3D程序已成為目前巖土力學(xué)計(jì)算中的重要數(shù)值方法之一，它采用ANSI C++語(yǔ)言編寫(xiě)，能夠進(jìn)行土質(zhì)、巖石和其它材料的三維結(jié)構(gòu)受力特性模擬和塑性流動(dòng)分析，并可通過(guò)調(diào)整三維網(wǎng)格中的多面體單元來(lái)擬合實(shí)際的結(jié)構(gòu)[5-6]。Flac3D采用的顯示拉格朗日算法和混合—離散分區(qū)技術(shù)能夠非常準(zhǔn)確地發(fā)現(xiàn)模擬材料中的塑性破壞和流動(dòng)。由于無(wú)需形成剛度矩陣，因此，采用較小的計(jì)算資源就能夠求解大范圍的三維巖土工程問(wèn)題[7]。

本文以高速公路隧道淺埋段地表塌陷為實(shí)例，利用數(shù)值模擬軟件Flac3D建立計(jì)算模型[8]，對(duì)地表塌陷原因、圍巖破壞位置、加固前后圍巖受力狀態(tài)進(jìn)行了詳細(xì)研究和對(duì)比[9]。并結(jié)合現(xiàn)場(chǎng)監(jiān)控量測(cè)和數(shù)據(jù)分析，驗(yàn)證了數(shù)值模擬結(jié)果，預(yù)測(cè)淺埋黃土段隧道圍巖長(zhǎng)期變形，確保注漿后隧道的穩(wěn)定和長(zhǎng)期運(yùn)營(yíng)安全[10-11]。

1 工程背景

隧道出口段為V級(jí)圍巖，圍巖破碎，埋深，呈層塊構(gòu)造，節(jié)理裂隙發(fā)育。隧道發(fā)生坍塌，初期支護(hù)全部垮塌并形成冒頂，地表下沉開(kāi)裂形成凹坑。

通過(guò)現(xiàn)場(chǎng)勘察和塌方原因分析，認(rèn)為導(dǎo)致塌方發(fā)生的原因有：（1） 隧道右洞掌子面為紅褐色粘土夾砂礫，夾石，拱部土質(zhì)松散；（2）坍方范圍埋深較淺，其中最淺處埋深僅，且地表存在偏壓；（3）本季雨水較多，雨水下滲導(dǎo)致土體重量增加，導(dǎo)致塌方冒頂。為了及時(shí)恢復(fù)施工秩序，決定對(duì)塌方區(qū)域進(jìn)行地表注漿加固處理，以提高圍巖強(qiáng)度，確保隧道安全、順利通過(guò)淺埋黃土段。

2 模型的建立

采用有限差分軟件Flac3D，根據(jù)現(xiàn)場(chǎng)地形及地層建立數(shù)值計(jì)算模型，模型長(zhǎng)140 m，寬60 m，高100 m；其中隧道寬14 m，高11.3 m。圍巖從上至下分別為：粉土層、砂土層、碎石土層及V級(jí)圍巖。隧道垂直穿過(guò)淺埋區(qū)域，淺埋段隧道拱頂埋深15m，淺埋區(qū)域沿掘進(jìn)方向長(zhǎng)40 m。數(shù)值模型如圖1所示。

模型邊界條件為：底部固定約束、四周約束法向位移、頂部為自由邊界。圍巖材料特性遵循Mohr-Coulomb屈服準(zhǔn)則，初支通過(guò)實(shí)體單元模擬，按均質(zhì)彈性體考慮。選取六面體單元對(duì)模型進(jìn)行網(wǎng)格劃分，假設(shè)模型在初始應(yīng)力條件下進(jìn)行開(kāi)挖，開(kāi)挖前進(jìn)行計(jì)算使模型不平衡率低于1e-5，得初始應(yīng)力場(chǎng)。計(jì)算荷載僅考慮圍巖自重及圍巖開(kāi)挖釋放荷載，根據(jù)現(xiàn)場(chǎng)工程勘探報(bào)告選取模型計(jì)算參數(shù)見(jiàn)表1。

3 開(kāi)挖模擬分析

當(dāng)隧道通過(guò)淺埋區(qū)域時(shí)，分別采用正臺(tái)階法開(kāi)挖，分析開(kāi)挖過(guò)程中圍巖的變形發(fā)展情況。開(kāi)挖按2m一循環(huán)開(kāi)挖，40m長(zhǎng)淺埋區(qū)域一個(gè)分為20個(gè)開(kāi)挖循環(huán)，每個(gè)開(kāi)挖循環(huán)計(jì)算1000時(shí)步。

圖2為圍巖剪應(yīng)變?cè)隽吭茍D。如圖所示，在采用正臺(tái)階法通過(guò)淺埋區(qū)域時(shí)，剪應(yīng)變?cè)隽糠植寂c采用全斷面法無(wú)異，同樣潛在剪切帶已經(jīng)從地表至開(kāi)挖面連通，存在剪切破壞可能。

臺(tái)階法開(kāi)挖工序多于全斷面法，開(kāi)挖時(shí)間亦更長(zhǎng)，在Flac3D計(jì)算程序中，圍巖的變形隨著時(shí)步的增加會(huì)不斷增加（變形速率不為0）。故采用正臺(tái)階法開(kāi)挖時(shí)，通過(guò)數(shù)值模擬得到的圍巖變形反而更大。這也說(shuō)明在采用正臺(tái)階法開(kāi)挖時(shí)，圍巖變形發(fā)展迅速而不能保持穩(wěn)定。開(kāi)挖過(guò)程中1#～6#監(jiān)測(cè)斷面拱頂監(jiān)測(cè)點(diǎn)數(shù)據(jù)如圖3所示。隨著隧道的開(kāi)挖，各監(jiān)測(cè)斷面拱頂沉降快速發(fā)展，其中3#監(jiān)測(cè)斷面（RK182+275）處拱頂下沉值達(dá)7.26m，已經(jīng)嚴(yán)重超出了圍巖變形破壞的控制范圍，圍巖出現(xiàn)大規(guī)模變形并導(dǎo)致破壞。圖4為開(kāi)挖過(guò)程中A～E監(jiān)測(cè)點(diǎn)的地表沉降監(jiān)測(cè)數(shù)據(jù)圖，其中B監(jiān)測(cè)點(diǎn)（RK182+280）的沉降值最大，達(dá)6.62m。

綜合上述分析，隧道采用正臺(tái)階法開(kāi)挖通過(guò)淺埋區(qū)域時(shí)，圍巖的變形趨勢(shì)、塑性狀態(tài)與剪應(yīng)變?cè)隽糠植寂c全斷面法開(kāi)挖時(shí)差異不大。圍巖變形發(fā)展依然迅速，出現(xiàn)塑性屈服并導(dǎo)致破壞。同時(shí)，由于圍巖變形速率較大并保持恒定，數(shù)值計(jì)算中拱頂沉降值以及地表沉降值隨著計(jì)算時(shí)步的增加而增加。

4 結(jié)論

基于隧道工程實(shí)際和Flac3D數(shù)值計(jì)算程序，對(duì)雙側(cè)壓淺埋段隧道開(kāi)挖過(guò)程進(jìn)行了數(shù)值模擬分析。

分析隧道在淺埋段趨于開(kāi)挖過(guò)程中圍巖的變化。選取全斷面法、正臺(tái)階法、正臺(tái)階留核心土法三種不同的開(kāi)挖方式進(jìn)行掘進(jìn)，分析淺埋段圍巖在開(kāi)挖過(guò)程中的變形規(guī)律及應(yīng)力狀態(tài)。結(jié)果顯示：通過(guò)三種方式開(kāi)挖進(jìn)入淺埋區(qū)域時(shí)，圍巖變形發(fā)展迅速、拱頂沉降和地表沉降值過(guò)大，遠(yuǎn)超過(guò)正常變形范圍，且存在潛在的剪切帶。因此可認(rèn)為，隧道在進(jìn)入淺埋區(qū)域后的開(kāi)挖會(huì)導(dǎo)致圍巖的破壞和地表土體的沉陷。

【參考文獻(xiàn)】

[1]劉耀儒，劉元高，周維垣，等.應(yīng)用 FLAC 方法進(jìn)行動(dòng)力分析[J].巖石力學(xué)與工程學(xué)報(bào)，2001，20（2）：1518-1522.

[2]張翾，張頂立，李鵬飛.黃土隧道地表塌陷原因分析與施工對(duì)策研究[J].隧道建設(shè)，2010，30（4）：430-433.

[3]喬明燦，樂(lè)金朝，牛向飛.隧道地表坍陷深層注漿理論研究與工程實(shí)踐[J].鐵道建筑，2011（9）：80-83.

[4]姜洪濤，譚兵，臧萬(wàn)軍，等.高速公路隧道淺埋段地表深層注漿施工技術(shù)[J].現(xiàn)代隧道技術(shù)，2008，45（3）：65-68.

[5]陳金山，王智學(xué)，戴前鋒.靈隱寺隧道軟弱圍巖淺埋偏壓段地表注漿加固技術(shù)[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì)，2003（1）：11-12.

[6]劉凡亮，王勇.地表注漿加固技術(shù)在淺埋松散段隧道中的應(yīng)用[J].四川建筑， 2012，32（2）：193-195.

[7]胡玉林.大南山隧道 淺埋段地表注漿加固技術(shù)[J].中國(guó)科技信息，2010（23）：64-66.

[8]譚毓浚.南嶺隧道下連溪地表坍陷的整治[J].鐵道建筑，1983（06）.

[9]王曉川.建筑基礎(chǔ)的注漿加固[J].探礦工程，1985（06）：61-62.

[10]李曉紅，姜德義，劉春，等.公路隧道穿越采空區(qū)治理技術(shù)研究[J].巖土力學(xué)， 2005，26（6）：910-914.

[11]James S. Mellett. Ground penetrating radar applications in engineering， environmental management， and geology[J]. Journal of Applied Geophyscis， 1995（33）： 157-166.

[責(zé)任編輯：楊玉潔]