田 偉,王文星,李 浩
(浙江省交通規劃設計研究院,浙江杭州310006)
在展線困難的山區公路中,連體隧道是最常見的短隧道形式。在偏壓地形中,連體隧道通常采用半明半暗工法進洞,但是左右線洞口里程往往難以一致,于是存在先開挖內側埋深較大側隧道還是先開挖外側埋深較淺側隧道的兩難選擇。由于施工工序不合理引發的滑坡、塌方等工程事故時有發生。因此,洞口施工階段是偏壓連拱隧道施工的主要風險點。
由于地質體具有結構復雜、空間單元復雜等特征,傳統基于極限平衡理論的分析方法往往勉為其難。地質建模及有限元數值分析方法考慮了邊坡巖體的非均質和不連續性,可以給出巖體的應力、應變大小與分布,避免了極限平衡分析法中將滑體視為剛體而過于簡化的缺點,能近似地從應力應變去分析邊坡的變形破壞機制。針對小凈距隧道進洞施工的穩定性研究較多,如陳思陽采用FLAC3D和強度折減法模擬分析短臺階、CD、CRD三種方法對大斷面黃土偏壓隧道邊坡穩定性的影響,并推薦采用CRD工法[1];石熊等進一步對CRD工法的各分部施工工序進行了分析,并建議山體外側先開挖[2],這一結論與王薇等的研究結論吻合[3];中巖墻(中夾巖)對小凈距隧道的穩定性至關重要,因此王帥帥、吳德興、侯瑞彬等學者開展了針對性的研究,包括中巖墻的全過程受力特征、合理厚度等[4-6]。
但偏壓連拱隧道的施工工序和受力更加復雜,更有必要建立精確的三維地質體結構模型,進行地質體的全過程施工力學分析。為此,以杭新景高速公路半塢隧道為例,建立邊仰坡體和隧道結構的三維地質結構模型,實現對滑坡體結構的三維可視化功能,通過對各個施工階段的位移、剪應力、剪應變、邊坡穩定性系數等指標的對比分析,提出合理的施工工序。
半塢隧道位于浙西侵蝕剝蝕中低山區,地層主要為第四系殘坡積和震旦系下統志棠組第三段粉砂巖、細砂巖。第四系殘坡積含角礫粉質粘土和粘性土角礫。含角礫粉質粘土,灰黃色,軟可塑狀,角礫含量20%~30%,局部達到40%~50%,粒徑以0.5~3cm為主,棱角狀,成份以粉砂巖為主,偶含有少量植物根系;粘性土角礫,灰黃色,稍密狀,碎石含量60%~70%,粒徑2~5cm,棱角狀,余為角礫及粘性土,膠結性較好,局部粘性土含量較高,達50%。震旦系下統志棠組第三段(Z1z3)粉細砂巖,灰色,中厚層—中薄層狀結構,巖質較硬。
根據杭新景高速公路半塢隧道的地質勘探資料及結構設計資料,建立了邊仰坡滑體的三維地質結構幾何模型(圖1),實現了對滑坡體結構的三維可視化功能,其中隧道線路軸線方向為X軸、水平垂直于線路軸線的方向為Y軸,豎向為Z軸。導入到ABAQUS軟件后對其進行三維網格劃分,優先采用四節點四面體單元和八節點六面體單元,劃分后的單元總數為254712個,節點總數為78459個。計算分析所采用的模型邊界條件:側向邊界采用水平法向約束;模型底部邊界采用三向固定約束。計算初始應力條件:因考慮到滑坡體主要為地表淺部巖土體,可以忽略構造應力影響,故在計算過程中僅考慮自重應力產生的初始應力條件。材料計算參數如表1所示,并采用摩爾—庫倫理想彈塑性模型。

圖1 半塢隧道及地層幾何模型

表1 隧道結構及地層物理力學參數表

表2 2個方案的施工工序表
由于隧道線路軸線與地形等高線夾角較小、地形偏壓顯著,左右洞暗洞進洞位置不一致。需要分析不同的進洞順序下,邊坡體的應力與位移場及相應的安全系數。暫且定義2種施工工序方案——方案A:左洞為先行洞,然后開挖右洞;方案B:右洞為先行洞,然后開挖左洞。2個方案的總體施工工序如表2所示。
以右洞先行為例,施工步可進一步細化為:開挖中導洞→澆筑中隔墻→右洞外側邊坡臺階開挖及支護→右洞外側擋墻澆筑→右洞立鋼拱架、并澆筑護拱→右面及修整。洞拱部暗挖并施做初期支護→右洞暗洞落底、仰拱開挖及澆筑回填→澆筑右洞二次襯砌、拱背土石回填→左洞拱部暗洞開挖、并施做初期支護→左洞暗洞落底、仰拱開挖及澆筑回填→澆筑左洞二次襯砌→溝槽、路
為研究不同的施工工序對邊坡體和隧道變形與穩定的影響,首先采用強度折減法計算得到了。然后分析比較2種工序方案施工過程中的位移場、剪應力場、最大剪應變增量、邊仰坡穩定性系數等的施工響應。
Y方向位移反映了邊坡的變形,當先行洞開挖及支護完成時,Y方向位移主要集中在左洞擋墻及洞口處,采用A、B工序方案時分別為5.8mm和1.4mm;X方向位移反映了仰坡的變形,采用A、B工序方案時分別為3.5mm和1.1mm。
施工完成后,Y方向位移主要集中在右洞洞頂及左洞入口處,采用A、B工序方案時分別為10.6mm和9.4mm;對應的X方向位移分別為5.8mm和4.7mm。可見,B工序方案對邊仰坡體的擾動程度要小于A方案。B工序方案施工完成后擋墻與襯砌結構合位移如圖2所示。

圖2 B工序方案施工完成后擋墻與襯砌支護結構合位移
土體剪應力值是摩爾庫倫模型中衡量土體是否將發生破壞的主要指標。當先行洞開挖及支護完成時,A工序方案的最大剪應力和剪應變主要發生在左洞擋墻墻趾處,分別約為0.6MPa和0.009,B工序方案的最大剪應力和剪應變主要發生在右洞洞頂處,分別為0.5MPa和0.007。兩洞施工完畢后,最大剪應力和剪應變分布范圍和量值基本一致,分別為0.8MPa和0.012。
2種工序方案對隧道結構剪應力的影響較小,當先行洞開挖及支護完成時,最大剪應力主要發生在明洞和先行洞半明半暗洞的洞腰部位,量值在655~665kPa之間。施工完成后,擋墻與隧道最大剪應力分布區域與范圍基本一致,最大值為1.0MPa。
2種工序方案對隧道位移有一定的影響,當先行洞開挖及支護完成時A、B工序方案的最大合位移值分別為1.21cm(左洞暗洞入口拱頂及擋墻位置處)和0.87cm(右洞洞頂處)。施工完畢后擋墻與隧洞支護結構合位移分布范圍與大小基本一值,A、B工序方案的最大合位移值分別為1.87cm和1.71cm。
邊仰坡穩定性分析采用強度折減法,其定義的安全系數為巖土體的實際抗剪強度與臨界破壞時的折減后抗剪強度的比值,其主要原理是對強度指標凝聚力c和摩擦角φ進行折減[(式(1)、式(2)],然后對邊坡進行數值計算分析,不斷地增加折減系數,反復計算,直至其達到臨界破壞,此時得到的折減系數即為安全系數FS。

式中:Ftrial——折減系數;
CF——折減之后的凝聚力;
φF——折減之后的摩擦角。
于是,可以計算初始狀態和各個施工步對應的邊仰坡穩定性系數(圖3),其中,天然狀態下邊仰坡體的初始穩定性系數為1.14;隨著1~5號抗滑樁、左幅路基回填及地基處理、挖除邊坡外土體并進行坡面防護等加固措施的開展,穩定性系數一度上升到1.39;而后,隨著外側邊坡削坡施工削弱了土體原有的穩定狀態,穩定性系數下降到1.23,仍大于初始狀態;施工左右明洞并回填洞頂后,穩定性系數又提高到1.42;隨著隧道開挖,穩定性系數再次下降,當先行洞開挖及支護完成時,A、B工序方案對應的穩定性系數分別為1.32和1.39;施工完畢后,A、B工序方案對應的穩定性系數分別為1.25和1.32。因此,2種工序方案下,邊坡體整體安全系數均較大,邊坡體均處于穩定狀態。然而不同的工況下,邊坡體的變形值與范圍存在一定差異,從控制變形的角度來考慮,建議以右洞為先行洞,然后開挖左洞。

圖3 2個方案的邊仰坡穩定性系數對比圖
通過建立半塢隧道邊仰坡三維地質結構模型,并對各個施工階段的位移、剪應力、剪應變、邊坡穩定性系數等進行全過程分析,對比分析了施工工序的影響。
(1)抗滑樁對于提高邊坡穩定的效果明顯,從數值上看提高了22%,對于偏壓嚴重的洞口是提高邊坡穩定性、減少刷坡高度和土方量的有效手段。
(2)施工過程中最危險的階段為澆筑偏壓當前的削坡開挖階段,施工過程中應重點關注。
(3)從控制變形的角度來考慮,建議以右洞為先行洞,然后開挖左洞。
半塢隧道在計算分析的指導下,按照信息化施工的理念精心施工,加強施工過程中隧道、坡體表面和深部變形監測,及時動態反饋監測信息,評估坡體穩定,確保了隧道安全、如期建成。
[1]陳思陽,朱彥鵬,李忠.大斷面黃土偏壓隧道開挖側向邊坡穩定性影響分析[J].現代隧道技術,2014(1):82-89.
[2]石熊,張家生,劉寶琛.大斷面淺埋偏壓隧道CRD法施工工序研究[J].現代隧道技術,2015(3):193-199.
[3]王薇,鄒江海,潘文碩,等.不同施工順序對陡坡偏壓小凈距隧道圍巖穩定性的影響研究[J].中國安全生產科學技術 ,2016(8):28-33.
[4]王帥帥,高波,全曉娟,等.不同地質條件下淺埋偏壓小凈距隧道圍巖穩定性研究[J].水文地質工程地質,2014(3).
[5]吳德興,李偉平,謝寶超,等.陡坡偏壓小凈距隧道施工擾動空間效應研究[J].鐵道科學與工程學報,2016(5):906-913.
[6]侯瑞彬,申玉生,陳明奎.基于強度折減法的淺埋偏壓小凈距隧道圍巖穩定性分析[J].鐵道標準設計,2014(4):55-60.