越嶺淺埋偏壓隧道支護系統(tǒng)承載特性分析

摘要：從地理、地形及地質(zhì)等方面，剖析了西南山區(qū)淺埋偏壓隧道的建造環(huán)境，針對該類隧道介質(zhì)參數(shù)、荷載的多變性及其襯砌結構力學傳遞機制的復雜性，基于襯砌結構穩(wěn)定的強度和變形控制原理，研究了應力和位移雙重等效的隧道介質(zhì)參數(shù)反分析模型及其優(yōu)化算法，并理順了其實施操作程序。然后提煉出偏壓隧道承載結構承載特性的表征參數(shù)及非均衡性刻畫指標。最后以有限元數(shù)值方法為手段，研究了某隧道淺埋偏壓段支護體系的力學響應，從中總結概化出淺埋偏壓隧道支護結構的承載特征，為類似工程設計和施工控制措施提供指導。

關鍵詞：淺埋偏壓隧道；建造環(huán)境；等效反分析；承載特性

中圖法分類號：TU458文章標識碼：A

Analysis of load-bearing characteristics of supporting system of the servants and shallow tunnel under asymmetrical pressure

SU Yong-hua，LIANG Bin，WU Wen-guo

(College of Civil Engineering ，Hunan Univ ，Changsha ，Hunan 410082 ，China )

Abstract:From the geography， topographical and geological， construction environment of shallow tunnel under asymmetrical pressure of the mountain in the southwest was analyzed. According to the denaturations of medium parameter and load， further more， the complexity of mechanical transmission mechanism of lining structure， the model of stress and displacement double equivalent back analysis for getting the parameter of tunnel medium and the optimization algorithm was studied， based on the strength and deformation control principle for the lining structure stability， and， the implementation procedures was straightened out. Then it wsa refined that the characteristics of characterization parameter and unbalanced depicting index for supporting structure of the shallow tunnel under asymmetrical pressure. In the end， the mechanical response of supporting system was studied by means of finite element numerical method， it was summaried that the bearing characteristics of support structure for the shallow tunnel under asymmetrical pressure， in order to providing guidances for similar engineering in design and construction control.

Key word: shallow tunnel under asymmetrical pressure;construction environment; equivalent back analysis; load-bearing characteristics

1 前言

在山區(qū)交通、水利水電等領域的隧道和地下結構工程中，淺埋偏壓一直是設計、支擋、施工、維護及穩(wěn)定性評價與分析理論中最具挑戰(zhàn)性的難題之一。因此吸引了許多學者從試驗、監(jiān)控、數(shù)值模擬、地層加固等不同方面進行探索，并取得了很大進展。

在試驗方面，劉鳳宣，關寶樹[1]從事了淺埋偏壓軟巖隧道圍巖的錨桿加固模型試驗研究；鐘新樵[2]通過模型試驗研究了土質(zhì)隧道形成偏壓的5大影響因素；王兵、謝錦昌[3]通過室內(nèi)模型試驗觀察了松動圍巖的力學行為；李育樞、李天斌、王棟等[4]開展了偏壓隧道的震動試驗平臺設計研究。

在地層加固方面，程新軍[5]等報告了梅河高速公路三斷嶺1號隧道鋼花管注漿預加固效果；來弘鵬[6]等報告了地表預注漿加固公路隧道淺埋偏壓破碎圍巖技術。

在開挖數(shù)值模擬方面，王祥秋[7]等基于動態(tài)監(jiān)

測與有限元數(shù)值分析，報道了某偏壓隧道圍巖位移變化規(guī)律；汪宏，蔣超[8]介紹了某淺埋偏壓隧道洞口坍方數(shù)值分析與處治結果。

此外，VIKEA[9]，SHIMOJIMAE[10]研究了淺埋偏壓連拱隧道的受力特征；楊小禮，李亮，劉寶琛[11]根據(jù)信息優(yōu)化理論建立原始數(shù)據(jù)信息與偏壓隧道結構穩(wěn)定性關系。還有眾多的技術人員，報道了大量的淺埋偏壓隧道襯砌、支護等方面的工程技術措施。

這些研究從不同側(cè)面對淺埋偏壓隧道結構穩(wěn)定機理的揭示和修建技術的進步起到了有力推動作用。但是，由于其復雜性，從根本上解決問題還有許多工作要做。因此，作者試圖從建造場地環(huán)境入手，探討西南部山區(qū)公路越嶺淺埋偏壓隧道結構問題的復雜性根源和相應的響應特征，為類似隧道結構的穩(wěn)定性分析和設計提供支撐。

2 西南山區(qū)淺埋偏壓隧道建造環(huán)境

很顯然，山區(qū)公路隧道主要功能是穿山越嶺。在隧道修筑中，隧道洞身通過山巒之間的溝壑地段或遇到地質(zhì)構造強烈影響帶時可能出現(xiàn)淺埋、偏壓或淺埋偏壓等現(xiàn)象。西南地區(qū)在地理上處于第二階梯過渡帶，在水文方面，大氣降雨量大，地表徑流復雜，淺部地層地下水量及滲流路徑變化多端。在地形上，起伏頻繁，坡度陡緩不一，溝壑縱橫。在地質(zhì)史上，屬于地質(zhì)活動多發(fā)地區(qū)，各級構造、褶皺形跡或各類殘余形跡密布，地應力場歷經(jīng)反復釋放、累積、反轉(zhuǎn)和變向等方式的變遷。在工程地質(zhì)上，地層組成復雜，根據(jù)大量的工程建設揭露，淺部地層布滿各種充填性質(zhì)和程度不一的空洞或泥、水溶洞[12]。淺部地質(zhì)成分復雜，表層和淺層通常是各種碎石、腐巖堆積體，和見水、見風后風化迅速的各類準軟巖或軟巖，在山地岸坡，更是受到自然的剝蝕，地質(zhì)體性質(zhì)變化大，不確定性程度高。

西南山區(qū)淺埋偏壓隧道的上述建造環(huán)境，從根本上造就了隧道介質(zhì)基本參數(shù)及荷載在空間的快速變化性和不確定性，以及隧道圍巖與襯砌結構力學傳遞方式及破壞機制的復雜性和綜合性。

3 隧道介質(zhì)參數(shù)的雙重等效確定方法

在隧道結構設計和穩(wěn)定性分析，首要工作是確定各類基本參數(shù)。由于巖體性質(zhì)的復雜性，其基本物理力學參數(shù)從解析上直接演繹確定是很困難的。但地質(zhì)體參數(shù)與其受荷響應兩者具有決定性的聯(lián)系。從安全理念出發(fā)，介質(zhì)受載響應表現(xiàn)為附加應力和變形。對于強度控制型，破壞根本原因在于應力過大，超過其強度極限；對于應變控制型，在于位移過大，超過其變形極限?；谧冃雾憫恚壳霸趲r土工程中已經(jīng)開發(fā)了巖土力學參數(shù)的位移等效反分析間接確定方法。

3.1 巖土參數(shù)的位移等效確定方法

位移反分析是以監(jiān)控量測位移數(shù)據(jù)為依據(jù)，以有限元數(shù)值試驗為手段，通過誤差限的控制，反復迭代調(diào)試確定隧道圍巖的有關參數(shù)。該參數(shù)可能不是地質(zhì)體的真實參數(shù)，它只與對象受荷響應產(chǎn)生的位移等效，滿足安全控制和使用要求，與現(xiàn)場監(jiān)控量測緊密聯(lián)系在一起。位移參考點通常選在拱頂、拱肩、最大跨度處等具有代表性的控制位置，如圖1所示，其基本過程與理論可參考相關文獻[13，14]。

3.2 位移-應力雙重等效反分析原理

對于西南山區(qū)公路交通淺埋偏壓隧道而言，根據(jù)上述特征可知，隧道圍巖及襯砌結構力學響應機制的復雜性和綜合性，其破壞力學行為受強度與應變的雙重控制。所以基本參數(shù)僅僅通過位移等效不能全面反映其實質(zhì)，因此提出位移-應力雙重等效確定方法，雙重等效方法基本過程如下。

監(jiān)測斷面測點布置示意圖

Fig.1 Distribution of measurement points

根據(jù)工程實際，設待求等效參數(shù)記為X=#61563;X1， X2， …， Xr#61565;T，Xk(k=1，2，…，r)為圍巖彈性模量E，泊松比μ，重度γ，內(nèi)聚力c，內(nèi)摩擦角φ等。

如圖1所示，不失一般性，對于位移等效，選取拱頂下沉(I#61475;點)、拱肩處（II#61475;，III#61475;點聯(lián)線）和最大跨度處（IV#61475;，V#61475;點聯(lián)線）收斂作為反分析的等效控制點。以各測點 的實測位移 與由計算得到的相對應測點的位移 (X)之間的誤差函數(shù)作為目標函數(shù)u(X)，即：

（1）

對于應力等效，圍巖中的應力通過軸力錨桿監(jiān)測，亦如圖1布置所示。同理，以各測點（選取拱頂Ⅰ#61475;，拱肩Ⅱ#61475;和Ⅲ#61475;，最大跨度Ⅳ#61475;和Ⅴ#61475;處）j的實測軸力 與由計算得到的相對應測點的軸力 (X)之間的誤差函數(shù)作為目標函數(shù)F(X)，即：

（2）

對于基于位移等效和應力雙重等效的反分析方法，根據(jù)數(shù)理學原理，將式(1)和式(2)無量綱化，建立如下方程：

（3）

這樣，求解圍巖等效物理力學參數(shù)即可轉(zhuǎn)化為求目標函數(shù)f(X)的極小值問題。

3.3 雙重等效實施程序及優(yōu)化算法

巖體參數(shù)是具體的工程參數(shù)，具有特定的工程條件。結合地下結構與巖體力學理論，圍巖等效參數(shù)式(3)的求解操作程序如下。

(1) 建立目標函數(shù)

根據(jù)圍巖級別(如表1所示)，確定等效待求參數(shù)xk（如圍巖重度γ，內(nèi)摩擦角φ，內(nèi)聚力c，彈性模量E及泊松比μ等）范圍的上下限xka，xkb。根據(jù)約束條件和式（3）建立優(yōu)化目標函數(shù)式（4）：

(2) 建立數(shù)值模型

綜合分析工程的幾何性質(zhì)、地層組成、地質(zhì)條件、原始應力場并根據(jù)地下結構力學理論，確定有限元數(shù)值模擬模型尺寸、邊界條件。

(3) 模擬試驗

將Xl=#61563;x1l，x2l，…，xrl#61565;輸入數(shù)值模型，啟動有限元模擬試驗，得出 (i=1，2，…，n)， (j=1，2，…，m)（l表示第l次模擬試驗）。

l=0表示初次模擬試驗，相應X0=#61563;x10，x20，…，xr0#61565;表示初次輸入模型的值。對于X0的取值，可結合表1中xk的可行域范圍，基于BQ值按線性插值確定。

(4) 一階優(yōu)化搜索求解

一階優(yōu)化方法[15]是通過對目標函數(shù)添加罰函數(shù)將約束問題轉(zhuǎn)換為非約束問題。利用目標函數(shù)和優(yōu)化變量罰函數(shù)的導數(shù)在設計空間進行搜索。在每次迭代中，通過梯度計算(用最速下降法或共軛梯度法)來確定其搜索方向，并用線性搜索法對非約束問題進行最小化。

將上述有約束問題變成無約束問題：

（5）

式中，Q（x，q）為無約束問題的目標函數(shù)；Px為優(yōu)化變量的罰因子，Pg、Ph、Pw為約束條件的罰因子，具體形式參見文獻[16]；#402;為目標函數(shù)值；#402;0為參考目標函數(shù)初始值；q為控制約束的參數(shù)。

對于第j次優(yōu)化迭代，先假定一個搜索方向矢量d(j)，則第j+1次迭代表示為：

（6）

式中，sj為線性搜索參數(shù)，采用黃金分割法和最小二乘法來確定，對sj的限制條件為：

（7）

其中，Sj*為當前線性搜索迭代步驟中最大的可能步長，Smax為線性搜索步長的最大值。

在初始迭代（j=0）時，先假定無約束目標函數(shù)的負梯度方向為搜索方向，即：

（8）

式中：q=1。

而對于其他任一步（ｊ>0），利用Polak-Ribiere的遞歸公式來確定其收斂方向，即：

（9）

（10）

當滿足收斂容差時，一階優(yōu)化迭代終止。收斂準則為： (11)

考慮到工程實際中各種誤差的存在，取ε=1#61620;10-2，即得圍巖等效力學參數(shù)。

與以往單一的位移反分析不同，模式(4)是基于變形和應力雙重等效反分析，這樣為淺埋偏壓隧道圍巖參數(shù)的有效性提供了保證。

4 支護系統(tǒng)承載特性表征參數(shù)及分析方法

在地下工程中，對于淺埋偏壓隧道通常采用曲墻拱形斷面。錨噴支護是優(yōu)先采用的支護結構形式，支護體系由圍巖和支護結構組成。由于幾何和材料的非線性復雜程度太高，對于淺埋偏壓隧道支護結構及其圍巖力學狀態(tài)的解析求解方法，目前還處于探索階段。但大量的工程實踐表明，基于有限元模擬等數(shù)值方法，可以獲得較高準確性的結果，而且代價低廉。因此本文采用有限元數(shù)值求解方法。

根據(jù)地下結構力學，對于曲墻拱形斷面，拱頂、拱肩、斷面最大跨度處及墻腳處等部位是位移、應力特征的代表點。這些部位的受力和變形(位移)狀態(tài)(圍巖位移、圍巖應力、噴射混凝土內(nèi)力和錨桿內(nèi)力等)對整個承載結構的穩(wěn)定程度具有指示性作用。對于常規(guī)應力場中的深埋隧道，受力和變形狀態(tài)在中心線兩側(cè)分布具有對稱性；而在淺埋偏壓隧道中，兩側(cè)對應位置的圍巖和支護結構的受力和變形在數(shù)值、方向和符號差異體現(xiàn)偏壓的總體特征，對應兩側(cè)的比值定義為偏壓比a：

令（12）

其中：分別表示中心線對應位置內(nèi)、外側(cè)的受力或變形大小。

同時，定義偏壓相對超差b：

令（13）

利用a和b兩個指標可對偏壓程度進行量化，其變化反映結構在偏壓條件下承載分布趨勢。

5工程實例

荊竹山隧道位于湘北某市，為單孔隧道，全長820m。采用復合式襯砌結構和新奧法施工。隧道斷面K5+250地質(zhì)剖面圖如圖2所示，左右嚴重不對稱，左低右高，自然坡度約為400，隧道覆蓋深度為11.35m~23.04m。

5.1 開挖、支護模擬和等效參數(shù)

圍巖參數(shù)雙重等效反分析和承載結構受力特征分析，本文簡化為平面問題。以K5+250斷面為對象，如圖2所示，具體反演參數(shù)有巖體變形模量E， 泊松比#61549;，重度 γ，內(nèi)聚力c及內(nèi)摩擦角#61546;。位移邊界條件為：兩側(cè)施加水平約束，下邊界施加豎向約束，上邊界為自由邊界，如圖3所示。

帶載荷和邊界條件的有限元計算模型

Fig.3 Finite element model with load and boundary onditions

對隧道開挖和支護的數(shù)值模擬分三步進行：首先對模型施加重力場和相應位移約束，得到初始地應力，然后依次模擬上臺階和下臺階開挖及支護，開挖初期支護應力釋放率取70%。

現(xiàn)場實測數(shù)據(jù)如下： = (1.76，0.263，0.916，0.324，0.602) mm，=21.284，10.158，17.345，27.335，12.956）kN。根據(jù)前述雙重等效參數(shù)及優(yōu)化算法，得到E=1.9973GPa， #61549;=0.30622， γ=20.784kN/m3，c=0.24905MPa， #61546;=25.062ordm;。

5.2 淺埋偏壓段支護結構受力特征分析

采用反演參數(shù)對K5+750的開挖和初期支護進行正演分析，主要從圍巖位移、圍巖應力、噴射混凝土內(nèi)力以及錨桿內(nèi)力等四個方面來進行。

5.2.1 圍巖位移、應力分析

各控制點的圍巖徑向位移（徑向位移指向洞身為正，背離洞身為負）如下表5所示：

圍巖大多處于受壓狀態(tài)，但在右拱肩至最大跨處圍巖則受拉。圍巖第一主應力的a值大多為負且絕對值小于1。圍巖左右兩邊的受力明顯不對稱且符號反轉(zhuǎn)。

5.2.2 噴層混凝土受力分析

對于噴層混凝土而言，一般以剪力控制為主，固不考慮其所受的軸力和彎矩。各控制點的噴層剪力如下表7所示。

噴層剪力的a值絕大多數(shù)大于1，內(nèi)側(cè)受力大于外側(cè)。內(nèi)側(cè)拱腳出現(xiàn)應力集中，為受力最為薄弱部位。

5.2.3 錨桿受力分析

各控制點的錨桿軸力（皆取錨桿的近端軸力，拉應力為正，壓應力為負）如下表8所示：

內(nèi)側(cè)拱腳處錨桿受力最大，拱頂至內(nèi)邊墻的錨桿處于受拉狀態(tài)，外拱肩、外邊墻及內(nèi)拱腳處的錨桿處于受壓狀態(tài)，而一般隧道支護結構中錨桿幾乎都是受拉且近似對稱。

隧道開挖支護完成以后，外側(cè)錨桿基本處于受壓狀態(tài)，內(nèi)側(cè)錨桿基本處于受拉狀態(tài)，可認為外側(cè)錨桿的主要作用是加固外側(cè)圍巖以抵制支護結構的擠壓，防止外側(cè)巖體失穩(wěn)；內(nèi)側(cè)錨桿的主要作用是加固拱頂及右側(cè)受拉區(qū)圍巖，防止受拉區(qū)圍巖局部失穩(wěn)。

綜上，淺埋偏壓隧道圍巖及支護結構的承載特性與一般隧道差別很大，最大的特點就是不對稱，左右兩側(cè)的受力及變形經(jīng)常相反，這就要求我們在設計及施工時慎重對待。對于偏壓程度較大的隧道（根據(jù)經(jīng)驗取a>3時），建議考慮分次分層噴射混凝土或增大噴層厚度、提高混凝土標號、合理選用鋼拱架等鋼支撐，選用強度更高的錨桿及增大錨桿間距等設計和施工措施。

6 結語

1) 本文根據(jù)位移反分析法存在的不足，提出了相對完整的適合于西南山區(qū)越嶺淺埋偏壓隧道圍巖等效力學參數(shù)確定的基于圍巖分級、圍巖位移-錨桿軸力雙重等效反分析方法及其實施程序。

2) 提出了衡量淺埋偏壓隧道偏壓程度的量化指標a和b，深入研究了越嶺淺埋偏壓隧道初期支護結構的受力及變形特征以及與一般隧道的差異。

3) 淺埋偏壓隧道支護結構的受力和變形非對稱且更為脆弱，這就要求我們在越嶺淺埋偏壓隧道設計和施工時考慮到其特殊性，根據(jù)偏壓程度對受力薄弱部位進行加強處理。

參考文獻(Reference)

[1] 劉鳳宣，關寶樹.淺埋、偏壓、軟巖隧道開挖與錨桿加固的試驗研究[J].巖石力學與工程學報，1991，10（2）:126-137.

LIU Feng-xuan，GUAN Bao-shu. Experimental study on the excavation and reinforcement of shallow-buried tunnel in soft rockmass and under uneven pressure[J]. Chinese Journal of Rock and Engineering， 1991，10（2）：126-137.

[2] 鐘新樵.土質(zhì)偏壓隧道襯砌模型試驗分析[J].西南交通大學學報，1996，31（6）：602-606.

ZHONG Xin-qiao. An analysis of model tests on earth tunnel under unsymmetrical loading[J].Journal of southwest jiaotong university，1996，31（6）：602-606.

[3] 王兵，謝錦昌.偏壓隧道模型試驗及可靠度分析[J].工程力學，1998，15（1）：85-92.

WANG Bing，XIE Jinchang. Model test and reliability analysis for an unsymmetrically loaded tunnel[J]. Engineering Mechanics， 1998，15（1）：85-92.

[4] 李育樞，李天斌，王棟等.黃草坪2#隧道洞口段減震措施的大型振動臺模型試驗研究[J].巖石力學與工程學報，2009，28（6）：1128-1136.

LI Yu-shu，LI Tian-bin，Wang Dong，et al. Large-scale shaking table test for vibration-absorption measures of portal section of Huangcaoping Tunnel No.2[J]. Chinese Journal of Rock and Engineering， 2009，28（6）：1128-1136.

[5] 程新軍，程煒.梅河高速公路三斷嶺1號隧道鋼花管注漿預加固效果的探討[J].鐵道建筑，2005，5：42-43.

CHENG Xin-jun，CHENG Wei. Study on effect of ganghua tube grouting pre-reinforcement of Sanduan Mountain Tunnel No.2 in Meihe highway[J]. Railway Engineering， 2005，5：42-43.

[6] 來弘鵬，謝永利，楊曉華.地表預注漿加固公路隧道淺埋偏壓破碎圍巖效果分析[J].巖石力學與工程學報，2008，27（11）：2309-2315.

LAI Hong-peng，XIE Yong-li，YANG Xiao-hua. Treatment effect analysis of shallow-buried crushed surrounding rocks under unsymmetrical pressure reinforced with surface pregrouting technology in highway tunnel[J]. Chinese Journal of Rock and Engineering， 2008，27（11）：2309-2315.

[7] 王祥秋，楊林德，高文華.高速公路偏壓隧道施工動態(tài)監(jiān)測與有限元仿真模擬[J].巖石力學與工程學報，2005，24（2）：284-289.

WANG Xiang-qiu，YANG Lin-de，GAO Wen-hua. Dynamic monitoring and FEM simulation analysis of an expressway tunnel with unsymmetrical loadings[J]. Chinese Journal of Rock and Engineering， 2005，24（2）：284-289.

[8] 汪宏，蔣超.淺埋偏壓隧道洞口坍方數(shù)值分析與處治[J].巖土力學，2009，30（11）：3481-3485.

WANG Hong，JIANG Chao. Numerical analysis of a collapsed portal for shallow tunnel under asymmetrical pressure and its treatment[J].Rock and Soil Mechanics， 2009，30（11）：3481-3485.

[9]VIK E A，SVERDRUP L，KELLEY A，et al.Experiences from environmental risk management of chemical grouting agents used during construction of the Romeriksporten tunnel[J]．Tunneling and Underground Space Technology，2000，15(4)：369—378．

[10] SHIMOJIMAE，TANAKAT，HOSO Y．Using sholrt-and long-term transients in seepage discharge and chemistry in a mountain tunnel to quantify fracture and matrix water fluxes[J]．Journal of Hydrology，2000，23(4)：142—161．

[11] 楊小禮，李亮，劉寶琛.偏壓隧道結構穩(wěn)定性評價的信息優(yōu)化分析[J].巖石力學與工程學報，2002，21（4）：484-488.

YANG Xiao-li，LI Liang，LIU Bao-chen. Evaluation on structure stability of unsymmetrically loaded tunnels using the theory of information optimization analysis[J]. Chinese Journal of Rock and Engineering， 2002，21（4）：484-488.

[12] 蘇永華，李翔，李志勇，等.山區(qū)公路建造環(huán)境級別的模糊粗糙評價方法[J].湖南大學學報:自然科學版，2010，37(5):1-8.

SU Yong-hua，LI Xiang，LI Zhi-yong，et al. Fuzzy rough set evaluation method for construction environment of highway engineering on mountainous area[J]. Journal of Hunan University:Natural Sciences， 2010，37(5):1-8.

[13]KAVANAGH K T，CLOUGH R W. Finite element application in the characterization of elastic solids[J]. International Journal of Solids Structures，1972，9(7)：11–23.

[14] Kirsten H A D. Determination of rock mass elastic moduli by back analysis of deformation measurement. In: Proceedings Symposium on Expliration for Rock Engineering ，Johan-nesburg，1976，1154-1160。

[15] 黎滿林，常曉林，周偉.基于有限元方法的拱壩體形優(yōu)化研究[J].水電能源科學，2003，21（2）：25-27.

LI Man-lin，CHANG Xiao-lin，ZHOU Wei.Research on configuration optimization of arch dam based on FEM[J].Water Resources and Power， 2003，21（2）：25-27.

[16] More J J，Wright S J.Optimization software guide[J]．

SIAM，Philadelphia，1993．