王 莎，朱俊松，王萬千，王蔚楠，陳建康

立洲水電站調(diào)壓井不同施工方案三維非線性有限元分析

王 莎1，朱俊松2，王萬千1，王蔚楠1，陳建康1

（1.四川大學(xué)水利水電學(xué)院，成都 610065；2.中國市政工程西南設(shè)計研究總院，重慶 610081）

就調(diào)壓井而言，其施工期的圍巖穩(wěn)定、結(jié)構(gòu)受力特性是設(shè)計中須深入研究的關(guān)鍵技術(shù)問題之一。以立洲水電站調(diào)壓井工程為研究對象，應(yīng)用三維非線性有限元法對其不同方案的開挖支護(hù)效應(yīng)、應(yīng)力變形規(guī)律、圍巖穩(wěn)定及襯砌結(jié)構(gòu)內(nèi)力分布等開展了較為深入系統(tǒng)的研究，揭示該工程調(diào)壓井不同施工方案井壁圍巖、邊坡的應(yīng)力變形、塑性區(qū)分布規(guī)律及其穩(wěn)定性，提出了推薦施工方案、建議優(yōu)化措施等。上述研究成果為該工程優(yōu)化設(shè)計提供了重要的參考依據(jù)。關(guān) 鍵 詞：調(diào)壓井；三維有限元；應(yīng)力變形；圍巖穩(wěn)定

在“十二五”期間，中國西部將有一批水電工程投入建設(shè)和運(yùn)營，其顯著特點：工程區(qū)水文與地質(zhì)環(huán)境復(fù)雜，高地應(yīng)力，高地震烈度，壩高庫大，地下廠房、調(diào)壓井等規(guī)模宏大。因此調(diào)壓井施工期的圍巖力學(xué)特性、受力條件及與之對應(yīng)的應(yīng)力變形特征和失穩(wěn)破壞規(guī)律、襯砌結(jié)構(gòu)受力特性也成為調(diào)壓井設(shè)計中的重大難點。因此，運(yùn)用有限元數(shù)值模擬對軟弱地質(zhì)條件下調(diào)壓井圍巖穩(wěn)定性及結(jié)構(gòu)的受力特性進(jìn)行研究具有工程指導(dǎo)意義。

1 工程概況

立洲水電站引水系統(tǒng)采用“一洞一室三管”的布置形式。調(diào)壓井及壓力管道位于廠房西側(cè)山脊上，南側(cè)為龔家溝，北東及東側(cè)為木里河。該工程調(diào)壓井為露頂圓筒阻抗式，最大開挖直徑為27.4 m，井筒高度為147.25 m；其周圍圍巖屬穩(wěn)定性很差的Ⅳ和Ⅴ類圍巖，即井筒自上而下97 m范圍內(nèi)為Ⅴ類圍巖，易發(fā)生傾倒破壞、楔形體剪切滑移破壞。針對立洲調(diào)壓井開挖斷面較大且圍巖自穩(wěn)能力差，結(jié)構(gòu)受力復(fù)雜等條件，運(yùn)用三維非線性有限元數(shù)值模擬，擬定3種方案進(jìn)行調(diào)壓室施工期圍巖穩(wěn)定性研究，通過分析，提出最優(yōu)施工方案，為該工程的設(shè)計和施工提供有價值的借鑒和參考。

2 調(diào)壓井圍巖開挖過程及支護(hù)措施模擬

2.1 調(diào)壓井圍巖分級開挖模擬方法

立洲調(diào)壓井開挖采用“開挖釋放荷載”方式，即先加載再開挖來模擬開挖作用效應(yīng)。首先由自重作用下的初始地應(yīng)力場確定第1級開挖邊界上的初始應(yīng)力｛ˉσ0｝，然后去掉第1級開挖巖體（一般將模型單元的剛度設(shè)置為0），將｛ˉσ0｝中開挖巖體作用于開挖面的法向正應(yīng)力及切向剪應(yīng)力反向施加于開挖邊界，具體計算第1級開挖釋放荷載｛q1｝見下式［1］，即

式中BT為巖體節(jié)點的應(yīng)力轉(zhuǎn)置矩陣。

可得到第1級開挖釋放荷載作用形成的洞室?guī)r體第1級附加位移場與附加應(yīng)力場｛▽u1｝，｛▽σ1｝，因此，第1級開挖完成時邊坡的一次應(yīng)力場和位移場為［2］：

由此可計算第2～N級開挖釋放荷載｛q2｝～｛qn｝，以及相應(yīng)分級開挖產(chǎn)生的附加位移場｛▽u2｝～｛▽un｝與附加應(yīng)力場｛▽σ2｝～｛▽σn｝，由此可推出第2～N級開挖完成時的各級應(yīng)力場為

每級開挖單獨(dú)作用形成的位移場為

第i級開挖井筒圍巖累積位移場為

調(diào)壓室圍巖在上述開挖過程中，只要確定了井筒巖體開挖高度和開挖面邊界，運(yùn)用有限元計算軟件［3］，即可求得開挖邊界上的等效釋放荷載、附加位移場和應(yīng)力場。

2.2 調(diào)壓井分級開挖支護(hù)措施的模擬

調(diào)壓井作為地下洞室開挖的一種，錨桿和錨索作為一種臨時和永久支護(hù)已得到廣泛應(yīng)用。加錨巖體的變形和強(qiáng)度參數(shù)會有較大提高，可提高錨桿加固區(qū)域圍巖的強(qiáng)度參數(shù)來模擬錨桿的加固作用。通常將錨桿加固區(qū)的單元的強(qiáng)度參數(shù)凝聚力C值提高，本文提高15%～30%。

在立洲調(diào)壓井分級開挖過程中，邊坡圍巖屬于覆蓋層，巖性差，需施加預(yù)應(yīng)力錨索來提高邊坡圍巖強(qiáng)度；對于預(yù)應(yīng)力錨索的模擬，采用施加一對鎖固力。

在進(jìn)行立洲調(diào)壓井的開挖過程模擬時，根據(jù)調(diào)壓井開挖完成后井壁破壞部位及發(fā)育程度，確定系統(tǒng)錨桿的布置方案。并在每級開挖完成后及時進(jìn)行系統(tǒng)錨桿支護(hù)，系統(tǒng)錨桿支護(hù)參數(shù)根據(jù)各調(diào)壓井周邊圍巖和邊坡的圍巖類型各不相同［4］，立洲調(diào)壓井主要設(shè)置Φ25，＠1.5 m×1.5 m，L＝6.0 m；Φ28，＠1.5 m× 1.5 m，L＝9.0 m的2種錨桿。

在調(diào)壓井分級開挖過程中，混凝土噴層的模擬是將它及時同級噴護(hù)，并與圍巖協(xié)調(diào)變形，共同承擔(dān)圍巖變形產(chǎn)生的應(yīng)力。

混凝土襯砌結(jié)構(gòu)的模擬是在各級開挖過程中，對立洲調(diào)壓井位于井筒上部的鎖口段滯后一級將混凝土襯砌形成，并與圍巖及混凝土噴層協(xié)調(diào)變形，共同承擔(dān)圍巖變形產(chǎn)生的應(yīng)力；而井筒襯砌結(jié)構(gòu)的模擬，采用“先開挖及時錨噴支護(hù)，后襯砌施工”的方式進(jìn)行模擬。

3 調(diào)壓井不同施工方案分級開挖和支護(hù)程序

針對立洲調(diào)壓室開挖內(nèi)徑大，所屬圍巖為Ⅴ或Ⅳ類，巖體質(zhì)量差等特點，擬定調(diào)壓室采用3種方案進(jìn)行分級開挖，見表1。

方案一，無噴錨支護(hù)條件下分級開挖；方案二、三，分級開挖的同時，及時跟進(jìn)噴錨支護(hù)，鎖口段混凝土襯砌滯后一級施工程序。

表1 立洲水電站調(diào)壓室計算方案Table 1 Computational schemes for the surge chamber of Lizhou Hydropower Station

其中調(diào)壓室開挖過程方案一分為27級，平均每8 m為一開挖層；方案二、三為28級；根據(jù)擬定的分級開挖順序和進(jìn)度安排進(jìn)行有限元模擬計算，具體開挖支護(hù)程序見圖1。分別計算模擬各方案調(diào)壓井自上而下的分級開挖過程，分析開挖施工各步驟圍巖的位移、應(yīng)力特征及可能出現(xiàn)的破壞部位、發(fā)育深度及分布范圍。

圖1 立洲調(diào)壓井分級開挖流程Fig.1 Stepped excavation for the surge tank of Lizhou Hydropower Station

4 有限元模擬范圍及材料參數(shù)

4.1 模擬范圍

立洲調(diào)壓井三維有限元計算選取調(diào)壓井結(jié)構(gòu)和較大范圍圍巖體作為整體研究對象［5］，其中有限元計算模型鉛直向底部取至1 810.00 m高程，約1.35倍井高，上部延伸至地表，以調(diào)壓井中心軸線為界，上游取210 m，下游側(cè)取250 m；前后側(cè)各取140.0 m（約5倍筒徑），井臺高程前后側(cè)計算邊界圍巖厚度約120.0 m。有限元計算坐標(biāo)系選定為：

x軸沿調(diào)壓井引水隧洞方向（上下游側(cè)），軸向方位N56.2°W；

y軸垂直調(diào)壓井引水隧洞水流方向（前后側(cè)），軸向方位N3°E；

z軸與x和y垂直，且沿井筒高程方向鉛直向上。

根據(jù)橫、縱剖面巖層分界線，地形等高線及調(diào)壓井輪廓，并考慮分級開挖程序進(jìn)行有限元三維建模；離散中鎖口混凝土、噴層、襯砌結(jié)構(gòu)和圍巖采用空間8節(jié)點等參實體單元，系統(tǒng)錨桿采用只計入軸向剛度的錨桿單元模擬。整個計算模型共剖分為100 327個節(jié)點和99 680個單元。壓井三維有限元計算網(wǎng)格見圖2所示。

4.2 材料參數(shù)

根據(jù)調(diào)壓室地勘資料和地質(zhì)剖面圖揭示，立洲調(diào)壓井2 058.47 m高程以上97.0 m范圍為Ⅴ類圍巖，（單位彈抗系數(shù)k0＝500 MPa／m，變形模量E0＝0.5～1 GPa），2 058.47 m高程以下至調(diào)壓井阻抗板底部為Ⅳ類圍巖，圍巖條件仍然很差。有限元計算模擬了調(diào)壓室的所處的三維地形、覆蓋層和幾個風(fēng)化層巖帶等。調(diào)壓室井筒襯砌混凝土、鎖口混凝土按C25混凝土考慮，貼坡和噴層混凝土為C20。其力學(xué)參數(shù)匯總于表2、表3。

圖2 立洲水電站調(diào)壓井有限元計算網(wǎng)格Fig.2 Finite element computing girds for the surge tank of Lizhou Hydropower Station

表2 調(diào)壓井圍巖與覆蓋層物理力學(xué)參數(shù)地質(zhì)建議值Table 2 Suggested values of physical and mechanical parameters for the surrounding rock and overburden layer of the surge tank

表3 錨桿、混凝土、鋼筋等材料物理力學(xué)指標(biāo)Tab le 3 Physical and mechanical indexes for the anchor rod，concrete，and steel

根據(jù)設(shè)計提供資料，可以確定立洲調(diào)壓室內(nèi)底板高程為2 011.25 m，最高涌浪水位高程為2 145.72 m，井外最高水位為2 083.45 m。

5 成果分析

通過對立洲水電站調(diào)壓井工程三維非線性有限元模擬分析，研究成果已應(yīng)用于實際工程設(shè)計，主要成果總結(jié)如下。

5.1 調(diào)壓井不同開挖方案圍巖變形分析

調(diào)壓室四周圍巖均向井內(nèi)變形，同一高程變形表現(xiàn)明顯不對稱性，上游側(cè)大于下游側(cè)，前、后側(cè)基本相等，呈“橢圓型”變形。變形受開挖面釋放荷載、井周圍巖類別及各側(cè)埋深控制。

調(diào)壓井不同方案下，縱向變形分布規(guī)律如表4所示。

表4 開挖高2 088.00 m三方案y＝0縱剖面洞壁最大變形對比Table 4 Com parison ofmaxim um longitudinal deformation of the lining in three schemes（h＝2 088.00 m and y＝0）

3個方案開挖變形規(guī)律大致相同，井壁變形上游側(cè)大于下游側(cè)，量值隨開挖高程降低先增后減，至高程2 088 m附近達(dá)到最大，主要是該高程附近區(qū)域處于風(fēng)化層與微新基巖交界帶，圍巖質(zhì)量變化較大所致。井筒開挖支護(hù)約束作用明顯，井壁水平向變形相對無支護(hù)減小，開挖面鉛直向變形變化甚小。鎖口段與井筒支護(hù)段，方案二、三變形規(guī)律與方案一基本一致，但都較方案一量值有所減小。

5.2 不同開挖方案應(yīng)力分析

3個方案大主應(yīng)力σ1整體上隨高程降低而遞增。方案一開挖面應(yīng)力最大值12.55 MPa（第27級開挖高程2 008.75 m與隧洞交匯處）；方案二噴層承受圍巖壓力，應(yīng)力達(dá)34.14 MPa；方案三應(yīng)力值相比方案二減小明顯，最大值18.81 MPa（高程2 008.75 m與隧洞交匯處）。

從小主應(yīng)力σ3分布來看，方案一井筒拉應(yīng)力分布范圍一般在開挖層中部，量值較小，不超過－0.5 MPa；方案二拉應(yīng)力普遍大于方案一，主要分布在當(dāng)級開挖層中部圍巖及鎖口段與圍巖交匯處，最大值－10.41 MPa發(fā)生在第15級開挖高程2 112m井壁上游側(cè)；方案三與方案二應(yīng)力分布規(guī)律相同，最大值－1.42 MPa發(fā)生在高程2 155 m上游側(cè)鎖口處。

3個方案應(yīng)力分布規(guī)律基本相同，同一高程井周應(yīng)力分布存在明顯不對稱性，即上游側(cè)大于下游側(cè)，左右側(cè)基本相同，噴層及鎖口砼內(nèi)壁壓應(yīng)力水平一般大于外壁。方案二在結(jié)構(gòu)突變處（鎖口砼厚度有變化）拉應(yīng)力很大，疑有應(yīng)力集中現(xiàn)象，方案三拉應(yīng)力值變化比較穩(wěn)定。

5.3 3個方案塑性區(qū)分布規(guī)律

3個方案圍巖塑性發(fā)育、分布規(guī)律大體相同，即井周圍巖塑性區(qū)分布不對稱，上游側(cè)圍巖埋深大，開挖卸荷效應(yīng)明顯，塑性區(qū)比較發(fā)育；下游側(cè)圍巖埋深相對較淺，無明顯塑性區(qū)；井壁左右兩側(cè)圍巖埋深接近，塑性區(qū)發(fā)育較為對稱；塑性區(qū)發(fā)育深度方案一最大，方案二、三基本相同，較方案一有大幅減小。從點安全系數(shù)等值線分布圖（見圖3至圖5）可知，開挖邊坡上游存在大量小于1.2的區(qū)域，支護(hù)、錨固效果明顯，提高了調(diào)壓室圍巖的穩(wěn)定性和安全性。從圖3至圖5可見，邊坡存在大量塑性區(qū)，建議增加抗滑樁。

圖3 方案一第27級（最后一級）開挖y＝0橫剖面塑性區(qū)分布、等值線圖Fig.3 Plastic zone distribution and contour map of the horizontal profile at the last step（27th）in scheme 1（y＝0）

圖4 方案二第28級（最后一級）開挖y＝0縱剖面塑性區(qū)分布、等值線圖Fig.4 Plastic zone distribution and contour map of the longitudinal profile at the last step（28th）in scheme 2（y＝0）

圖5 方案三第28級（最后一級）開挖y＝0縱剖面塑性區(qū)分布、等值線圖Fig.5 Plastic zone distribution and contour map of the longitudinal profile at the last step（28th）in scheme 3（y＝0）

6 結(jié) 論

本文結(jié)合立洲水電站調(diào)壓井工程三維非線性數(shù)值模擬研究課題，擬定3種施工方案，并對不同方案下的井筒結(jié)構(gòu)應(yīng)力、變形，以及圍巖穩(wěn)定性進(jìn)行了較為深入的分析和探討。分析結(jié)果表明：圍巖變形，塑性區(qū)分布，方案二、三均與方案一有明顯相似的分布規(guī)律，但在量值上均較方案一有所減小。在應(yīng)力分布方面，方案二的小主應(yīng)力普遍大于方案一，方案三的應(yīng)力分布雖與方案二規(guī)律相同，但比方案二?。淮笾鲬?yīng)力方案三比方案二明顯減小。方案三施工過程襯砌厚度與錨桿長度均小于方案一、二，節(jié)約了施工成本。綜合考慮以上因素，選取方案三為最優(yōu)施工方案。上述研究成果，為該工程的優(yōu)化設(shè)計提供了重要的參考依據(jù)。

［1］ 周榮祥，王秋生.用ANSYS模擬巖土開挖的兩種方法［J］.甘肅科技，2005，21（10）：135－138.（ZHOU Rong-xiang，WANG Qiu-sheng.Two ANSYS-Based Methods for Excavation Simulation［J］.Gansu Science and Technology，2005，21（10）：135－138.（in Chinese））

［2］ 毛擁政.四川省尼日河開建橋電站工程調(diào)壓井圍巖及支護(hù)結(jié)構(gòu)穩(wěn)定性分析及優(yōu)化研究［D］.西安：西安理工大學(xué)，2006.（MAO Yong-zheng.Stability Analysis and Optimization of the Surge Shaft Surrounding Rock and Supporting Structure for Kaijianqiao Hydropower Station on Niri River in Sichuan Province［D］.Xi’an：Xi’an University of Technology，2006.（in Chinese））

［3］ 楊興國，符文熹，曹 竹.復(fù)雜地質(zhì)條件下大型調(diào)壓井關(guān)鍵施工技術(shù)數(shù)值仿真優(yōu)化分析［J］.四川大學(xué)學(xué)報（工程科學(xué)版），2006，38（4）：5－9.（YANG Xing-guo，F(xiàn)UWen-xi，CAO Zhu.Numerical Modeling on Key Construction Schemes of Large-scale Surge Tank under Complex Geological Conditions［J］.Journal of Sichuan University（Engineering Science Edition），2006，38（4）：5－9.（in Chinese））

［4］ 李啟升.高邊墻地下調(diào)壓室襯砌結(jié)構(gòu)工作性態(tài)研究［D］.南京：河海大學(xué)，2006.（LIQi-sheng.Work Behavior of Lining Structure of Underground Surge-chamber with High SideWall［D］.Nanjing：Hohai University，2006.（in Chinese））

［5］ LEE PL，ZHOUW，CAMERON IT.Constrained Generic Model Control of a Surge Tank［J］.Computers and Chemical Engineering，1991，15（3）：191－195.

（編輯：曾小漢）

3-D Nonlinear Finite Element Analysis of Surge Tank of Lizhou Hydropower Station in Different Construction Schemes

WANG Sha1，ZHU Jun-song2，WANGWan-qian1，WANGWei-nan1，CHEN Jian-kang1
（1.College ofWater Resources and Hydropower，Sichuan University，Chengdu 610065，China；2.Southwest China Municipal Engineering Design Research Institute，Chongqing 610081，China）

The surrounding rock stability and structural load bearing features of surge tank in construction period are key technical issues in the design.The surge tank of Lizhou hydropower station in different construction schemes is taken as an example to research the excavation supporting effect，stress deformation law，surrounding rock stability，lining structure stress distribution and reinforcement.The research is carried out through 3-D nonlinear finite elementmethod.The stress deformation and plastic zone distribution and stability of the lining surrounding rock，and slope stress in different schemes are revealed.Construction scheme and optimizationmeasures are recommended.The research result provides reference for the optimization design of the project.

surge tank；3-D finite element；stress deformation；surrounding rock stability

TV732.51

A

1001－5485（2012）12－0083－05

10.3969／j.issn.1001－5485.2012.12.017 2012，29（12）：83－87

2011－10－19；

2011－11－22

王 莎（1987－），女，遼寧鞍山人，碩士研究生，主要從事山地災(zāi)害與防治研究，（電話）15982296885（電子信箱）694193545＠qq.com。