大型地下洞室群施工開挖方案及圍巖穩定分析

周樹兵，童 偉，鄧樂榮

(1.四川岷江水利電力股份有限公司，四川成都 611830;2.中國水電顧問集團成都勘測設計研究院，四川成都 610072)

1 前 言

在大型水利水電工程中，地下廠房洞室群的規模往往較大，結構復雜，大多由兩大洞室或三大洞室以及它們之間許多聯絡的通道組成。在復雜地層中開挖這類地下工程群體，洞室的開挖方案對圍巖穩定的影響至關重要。然而，地下洞室的施工開挖，從力學上講，這一過程往往是不可逆的非線性演化過程，它的最終狀態(或最終解)不是唯一的，而是與過程相關，或者說與應力路徑或應力歷史相關的［1］。因此，對已定巖體環境下，不同的地下洞室開挖方案就意味著在時空上以不同的方式對圍巖施加荷載，從而決定施工期內圍巖的洞周位移和塑性區的分布。

本文基于彈塑性圍巖理論，采用三維非線性有限元法，在確定的地應力環境條件下，研究西南某工程的主廠房、主變洞、尾調室等主要洞室分級開挖過程中圍巖的變形和破壞發育程度。同時，對3種開挖方案的異同進行分析比較，從開挖過程圍巖穩定的角度對不同的開挖方案的可行性進行了評價，并優選最佳方案。

2 工程簡介

西南某工程地下發電廠房樞紐是一大型地下洞室群體。其中主廠房洞室開挖尺寸為229.45m×27m×77.44m;主變洞開挖尺寸為179m×15.8m×19m;尾水調壓室開挖尺寸為172.5m×19.3m×99.6m。另外還有5條發電引水洞、5條母線洞、5條尾水管以及3條尾水洞，使整個地下洞室群龐大且錯綜復雜。整個地下洞室群布置區總體地形為向坡外傾斜，地勢西高東低、地面高程一般為700～820m，廠房洞室群埋深一般為260～315m，其NE側端墻距河岸約250m。地下洞室群穿過的巖層主要分別為二迭系吳家坪組(P2w)、茅口組(P1m)、棲霞組(P1q)灰巖為Ⅰ－Ⅱ類圍巖，梁山組(P1L)、志留系韓家店組(S2h)粘土巖為Ⅴ類圍巖，其它為Ⅲ－Ⅳ類圍巖。主廠房圍巖為層，除層及層圍巖中炭、泥質生屑灰巖為中硬巖外，其余均為堅硬巖。巖層走向30°～35°，傾向 NW，傾角 45°～48°，主廠房軸線與巖層走向交角40°～45°。主廠房部位分布有較多的規模較大的層間錯動帶，寬度一般為5～20cm。

3 計算模型

對地下洞室群進行了三維彈塑性有限元計算模擬，計算范圍取450m×400m×400m，計算區域包含了主廠房、主變洞、尾水調壓室、尾水管及尾水洞漸變段，同時在計算模型中模擬了穿越地下洞室群的2條層間錯動帶Fb54、Fb59和1條破碎夾層Ⅲ01。X軸與廠房縱軸線垂直，指向下游為正;Y軸鉛直向上為正;Z軸與廠房縱軸線平行，與X、Y軸成右手坐標系。本模型采用八節點六面體單元，對重點研究的部位進行網格加密。三維計算區域總共剖分177 729個單元，節點總數為187 023個。三維整體模型計算網格見圖1。地下洞室群布置見圖2。

圖1 三維整體計算網格

圖2 地下洞室群布置

計算中將計算區模型的上游邊界、左面邊界定為位移約束邊界，約束水平方向的位移;模型的底面也為約束邊界，僅約束垂直方向的位移;模型的右面邊界、下游邊界為荷載邊界;模型的頂面施加由上部巖層產生的重力。根據地勘報告并進行了三維地應力場回歸反演得出地應力場，計算時鉛直方向按自重應力場施加;水平向應力取側壓系數Kx=0.8，Kz=1.2，其中X方向為順河流方向，Z方向為垂直河流方向。

計算選用Mohr－Coulomb彈塑性材料模型。屈服準則采用Mohr－Coulomb準則。具體的計算參數如表1所示:

表1 計算采用的巖體力學及材料參數

4 初擬地下洞室群開挖方案

由于施工條件的限制，洞室群的開挖不可能全斷面一次開挖成型，而必須采取分期分塊的開挖方式，對于復雜地下洞室群的施工實際上是一個非線性加載過程，即是巖體施工力學的動態過程［2］。因此地下洞室群不同部位分期分塊開挖的組合對洞室的完建和穩定有十分重要的影響。考慮地下洞室群施工工期仿真模擬研究成果，結合施工現場具體條件以及考慮經濟、快速等因素，初步擬定3種代表性的開挖方案并進行三維彈塑性計算分析。

圖3 方案1 開挖順序示意

方案1:先開挖主廠房1步后，再同時開挖主廠房和尾水調壓室，后開挖主變室。

方案2:先開挖主廠房，再開挖主變室，后開挖尾水調壓室(擬二灘水電站地下廠房實際開挖程序)。

方案3:先開挖主廠房1步后，主廠房、主變室和尾水調壓室同時開挖。

3種開挖方案分期分塊開挖施工示意圖見圖3～5。

圖4 方案2 開挖順序示意

圖5 方案3 開挖順序示意

5 不同開挖方案的圍巖變形與破壞發育程度

5.1 圍巖的變形特征

由于地下廠房洞室群埋深較大，3種不同的開挖方案圍巖的位移場規律基本相似:由于應力釋放，圍巖產生指向洞室內部的回彈變形，頂拱及底板以豎向位移為主，邊墻以水平位移為主。同時，隨著洞室群的不斷開挖，受擾動圍巖的范圍不斷增大，洞室群最大位移量值也從低量值向高量值變化。另外，隨著開挖斷面的不斷增大和母線洞、尾水管的開挖，各個洞室的變形影響區從相對獨立發展到各洞室的變形場連接為統一的洞室群變形場，表現出等值線向洞室群中心偏轉的特征。3種方案3大洞室開挖后周邊各部位最大位移計算值見表2。

由表2中的數據可知:

(1)在3種開挖方案中，主廠房邊墻的最大位移值均發生在2號機組中心線附近下游邊墻中部，因為斷層Fb54在此部位出露，加之母線洞的開挖增加了位移變形的臨空面，所以在此處產生較大的位移。主變洞邊墻最大位移值也都出現在主變洞與母線洞的兩洞交接臨空面較大處。

表2 3大洞室開挖后周邊各部位最大位移值 cm

(2)主廠房頂拱的最大位移值和底板回彈變形值沿廠房縱軸向變化不大，變幅都在5mm之內，但是量值上都小于邊墻最大位移值。

(3)主變洞3種開挖方案中上下游邊墻最大位移值均小于頂拱的最大位移值，這是因為主變室的邊墻高度相對于主廠房和尾水調壓室較低的緣故。

(4)尾水調壓室的跨度與邊墻的比值相對于主廠房和主變室較小，因此，在3種開挖方案中尾水調壓室的頂拱和邊墻的最大位移值均小于另外兩大洞室。

(5)對比方案2和方案1、方案3可以看出，方案2的3大洞室的周邊最大位移值都較另外2種方案要大，這是由于方案2的分步開挖步數要比另外2種方案多了6步或8步，其對圍巖的擾動較大，圍巖的位移速度和應變速度越慢，從而洞室周邊位移值較大。

5.2 圍巖的破壞發育情況

圍巖的破壞發育形式主要有3種:

(1)完整巖體的拉裂破壞(可能呈現單向、雙向、三向拉裂狀態);

(2)沿優勢裂隙面方向的定向破壞;

(3)沿軟弱結構面的破壞［3］。

圍巖的穩定性在很大程度上是受結構面控制的，同時巖體及結構面參數的選取對于判斷圍巖穩定性也起決定性作用。

不同的開挖方案計算結果表明，3種方案最終的塑性區分布規律基本上相同。主廠房和主變洞的周邊均存在塑性區，尾水調壓室的高邊墻也出現塑性區。在開挖過程中，頂拱出現的塑性區分布并不隨著后期施工的影響而擴展，但邊墻的塑性區受后期開挖施工的影響較大，因為隨著整個洞室群分步開挖的進行，高邊墻和各洞室交叉部位逐漸開挖形成多個臨空面，從而使塑性區得以發展并快速擴大，并形成塑性貫通區，這說明該工程施工過程主要是高邊墻的穩定問題。各方案開挖完成后塑性區分布示意見圖6。

5.3 圍巖穩定性總體評價及方案選擇

采用三維彈塑性對3種不同的開挖方案進行分析計算表明，西南某工程地下廠房洞室群開挖后洞室周邊圍巖的位移和塑性區分布符合一般規律。頂拱、邊墻和底板等工程部位的最大位移值均在工程安全允許范圍內。但應注意以下幾個方面:

(1)開挖過程中，前期開挖斷面上的位移受后期施工的影響而變化。同時位移值的大小也受施工開挖步數的影響，開挖步數越多，對圍巖的擾動越大，圍巖的位移速度和應變速度越慢，洞室周邊位移值也較大。

圖6 各方案開挖完成后塑性區分布示意

(2)開挖頂拱時，拱頂出現的塑性區并不隨后期施工的影響而擴展，而邊墻塑性區的擴展易受其影響，甚至會出現大面積的貫穿。所以在實際施工中應注意觀測和支護。

(3)斷層對塑性區分布有一定的影響。在開挖過程中，斷層穿過的部位會首先出現塑性變形，并在此基礎上開展。從結果來看，斷層穿過的地方也是塑性區連通的部位。因此在開挖斷層出露部位時應進行預處理和加固措施。

(4)在洞室群的兩洞相交處開挖時，由于會形成較大臨空面，易產生較大的變形，特別是母線洞的開挖，極易會使主廠房和主變洞的塑性區相貫通，所以應加強觀測和重點支護。

(5)從3種開挖方案最終的塑性區面積來看，方案3的面積相對較小，所以從空間角度上方案3對圍巖穩定較為有利。

6 結 論

地下洞室群3種開挖方案計算分析成果表明，采用從上而下的分期開挖的施工方式，洞室圍巖的變形穩定性未出現異常現象。但是在洞室群的兩洞相交處，臨空面較大以及巖體節理裂隙發育的部位，圍巖穩定問題相對突出。因此，在施工中應著重監測這些部位的變位，注意及時采取支護處理。

從總體上看，3種開挖方案隨著開挖級數的增加，對圍巖的擾動越大，圍巖變位和破壞發育程度也隨之增大。因此在施工條件許可的情況下，應盡量減少開挖次數以及從空間上錯開洞室間施工的相互影響。

［1］于學馥，鄭穎人，劉懷恒，等．地下工程圍巖穩定分析［M］．北京:煤炭工業出版社，1989．

［2］鄭穎人，功效南．巖土塑性力學基礎［M］．北京:中國建筑工業出版社，1989．

［3］張學言．巖土塑性力學［M］．北京:人民交通出版社，1993．