某水電站圍堰三維滲流有限元分析

劉娜，閆兆杰，溫泳

(長(zhǎng)春工程學(xué)院土木工程學(xué)院 ，長(zhǎng)春 130012)

主要研究巖土數(shù)值計(jì)算。

摘要：利用三維有限元軟件對(duì)該圍堰進(jìn)行了滲流分析，通過(guò)分析得到：圍堰在下游溢出口處水力坡降較大，容易發(fā)生水力破壞；而滲流量主要集中在防滲墻下方的覆蓋層中，這主要是懸掛式防滲墻引起的。在實(shí)際工程中一定要注意這兩點(diǎn)，以確保該圍堰的穩(wěn)定性。

關(guān)鍵詞：滲流；有限元；水力破壞；滲流量

doi:10.3969/j.issn.1009-8984.2015.03.013

收稿日期：2015-07-24

作者簡(jiǎn)介：劉娜(1982-)，女(漢)，內(nèi)蒙古，講師

中圖分類號(hào)：TV223.4獻(xiàn)標(biāo)志碼：A

0工程背景

該水電站為大渡河干流水電梯級(jí)開(kāi)發(fā)的第10級(jí)電站，壩址位于甘孜藏族自治州康定縣境內(nèi)，沿大渡河右岸省道S211線向上游約85 km至丹巴縣城，向下游約50 km至瀘定縣城。工程采用斷流圍堰、隧洞泄流、大壩基坑全年施工的導(dǎo)流方式。上游圍堰為斜墻圍堰，堰頂高程為1 530.50 m，最大堰高56.50 m，堰頂寬10.00 m，迎水面邊坡1∶2.5，背水面邊坡1∶1.75；堰基防滲采用懸掛式防滲墻型式，防滲墻深度為65 m，防滲墻厚0.8 m；堰體防滲采用斜墻型式，最大擋水水頭約54.0 m。

1堰基與堰肩地質(zhì)特征

兩岸堰肩巖體主要為淺灰色—灰白色黑云母花崗巖，穿插有細(xì)晶花崗巖脈及石英閃長(zhǎng)巖團(tuán)塊，巖體堅(jiān)硬較完整，無(wú)規(guī)模較大的控制性軟弱結(jié)構(gòu)面分布，邊坡整體穩(wěn)定。風(fēng)化卸荷總體較弱，上、下游圍堰兩岸均為弱風(fēng)化，右岸強(qiáng)卸荷水平深度0～6 m；左岸強(qiáng)卸荷水平深度15～20 m。該工程截流后即開(kāi)始圍堰防滲墻的施工及堰體的堆筑，并要求在一個(gè)枯水期內(nèi)完成圍堰施工并達(dá)到擋水要求。由于堰體堆筑量大，基礎(chǔ)防滲深度大，且要求在一個(gè)枯水期建成擋水，涉及提前發(fā)電和度汛安全等重大問(wèn)題，工期緊、難度大，故對(duì)圍堰滲流量、堰體及基礎(chǔ)滲透穩(wěn)定做了三維有限元研究。通過(guò)對(duì)堰體及防滲結(jié)構(gòu)在設(shè)計(jì)工況下滲流和滲透穩(wěn)定的分析，充分論證圍堰防滲結(jié)構(gòu)的可行性和合理性。

2計(jì)算方案

采用的坐標(biāo)：上游圍堰計(jì)算中均采用笛卡兒直角坐標(biāo)系，以順河向?yàn)閤軸，下游為正方向；以垂直向?yàn)閥軸，向上為正向；以橫河向?yàn)閦軸，向右為正向。計(jì)算坐標(biāo)原點(diǎn)選取在防滲墻軸線、河道中線、堰底平面的交匯處。

圍堰為防滲墻圍堰，堰頂高程為1 530.50 m，最大堰高為56.50 m，堰頂寬10.00 m，最大擋水頭為51.0 m。

計(jì)算區(qū)域的選擇：計(jì)算區(qū)域順河向取為450 m(堰體順河向長(zhǎng)為246 m)，垂直向自壩底以下取為120 m，其中基巖部分為50 m，橫河向取350 m寬(堰體軸線橫河向?qū)挒?74 m)。此計(jì)算中迎水面水位51.0 m，背水面水位0.00 m。防滲墻深度為65 m進(jìn)行計(jì)算。

圖1　圍堰模型

3計(jì)算結(jié)構(gòu)分析

3.1　水頭分布

從圖2可以看出在正常水位時(shí)，等水頭線向防滲墻處靠攏，過(guò)了防滲墻后逐漸散開(kāi)，說(shuō)明防滲墻對(duì)整個(gè)滲流場(chǎng)起到了良好的防滲效果。鑒于三維云圖的立體觀測(cè)效果欠佳，防滲墻兩側(cè)的水頭變化又是關(guān)心的重點(diǎn)部位，所以給出了垂直于堰軸線的順河向的直線上不同深度處的各點(diǎn)的水頭值。這條直線過(guò)點(diǎn)(-200，y，0)和(300，y，0)。如圖3所示，分別為坐標(biāo)y為-5 m，-20 m，-45 m，-50 m時(shí)，即堰基的不同高程處水頭值，從圖上可知水頭值自上游至下游逐漸降低，且在防滲墻處有突降；堰基中隨著高程線的降低，水頭線自上游至下游變化越來(lái)越平緩，在防滲墻處的突降幅度越來(lái)越小。當(dāng)堰基中高程線降低到防滲墻底部時(shí)，水頭線自上游至下游變成一條平滑曲線。

圖2　圍堰等水頭面

(a)y=-5 m處

(b)y=-20 m處

(c)y=-45 m處

(d)y=-50 m處 圖3　堰底處順河向水平線上各點(diǎn)的水頭

3.2　水力坡降

水力坡降與土體的滲透破壞關(guān)系密切，通過(guò)計(jì)算水力坡降可以較容易地判斷出容易發(fā)生水力破壞的區(qū)域。同理，取垂直于堰軸線的順河向的直線上不同深度處的各點(diǎn)水力坡降的值進(jìn)行分析。由圖4可得出以防滲墻為中心，水力坡降呈現(xiàn)中心處高，兩邊低的尖峰分布形態(tài)，在防滲墻處呈現(xiàn)頂部水力坡降高，底部水力坡降低的形態(tài)，而防滲墻的垂直下方的水力坡降又低于防滲墻中的分布形態(tài)。從上圖4可以看出，當(dāng)y=-5 m時(shí)，在下游200 m處水力坡降較大，可以達(dá)到0.9，隨著深度的增加，此處水力坡降會(huì)有所降低，但仍可達(dá)到0.25左右，所以，在下游200 m左右處需要特別小心，這位置水力坡降較大，在較淺的位置很容易發(fā)生管涌，在下游滲流溢出口處應(yīng)設(shè)置排水棱體，增加排水。通過(guò)計(jì)算說(shuō)明本工程采用的是懸掛式防滲墻，雖然沒(méi)有封閉式防滲墻的截滲效果好，但起到了延長(zhǎng)滲透路徑的作用，在一定程度上達(dá)到了降低水力坡降的效果，保證了整個(gè)堰體和堰基的滲流穩(wěn)定。

(a)-5 m處

(c)-45 m處

(b)-20 m處

(d)-50 m處

3.3　滲流量

通過(guò)分析得出上游圍堰防滲墻處截面滲流量和圍堰壩軸線處截面滲流量表。表1中左、右堰肩、堰底覆蓋層、基巖、堰體所表示的區(qū)域如圖1所示。由表1可以看出，整個(gè)圍堰區(qū)域的滲流量中，防滲墻下方的覆蓋層中的滲流量最大，占整個(gè)滲流量的絕對(duì)份額，左堰肩和右堰肩中的滲流量所占份額相對(duì)不大，這種滲流量分布規(guī)律是懸掛式防滲墻圍堰的特征。當(dāng)圍堰采用懸掛式防滲墻，尤其是防滲墻較短時(shí)，左堰肩和右堰肩中的滲流量相對(duì)很小，主要滲流量在防滲墻垂直下方的覆蓋層中，對(duì)滲流量的控制重點(diǎn)應(yīng)放在這一區(qū)域。

表1　滲流量統(tǒng)計(jì)表　　m 3/s

4結(jié)語(yǔ)

該圍堰滲流計(jì)算結(jié)果表明：此圍堰采用斜墻—懸掛式防滲墻防滲體系，懸掛式防滲墻可以起到延長(zhǎng)滲透路徑、降低水力坡降的作用。但由于懸掛式防滲墻沒(méi)有深入基巖，圍堰在下游接近出口處-5 m深的范圍內(nèi)水力坡降較大，容易發(fā)生水力破壞，所以在實(shí)際工程中特別需要注意此處，應(yīng)該在下游壩坡滲流出口處做相應(yīng)的處理，可以鋪設(shè)反濾層并設(shè)置排水棱體，以防止發(fā)生管涌破壞。由于懸掛式防滲墻不能深入基巖，主要滲流量在防滲墻垂直下方的覆蓋層中，左堰肩和右堰肩中的滲流量相對(duì)很小，所以對(duì)滲流量的控制重點(diǎn)應(yīng)放在覆蓋層中，同時(shí)應(yīng)注意滲流量過(guò)大影響施工的問(wèn)題。

參考文獻(xiàn)

[1] 王保田.懸掛式防滲墻防滲效果的模擬試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2008,27(增1):2766-2771.

[2] 楊秀竹， 陳福全，雷金山，等.懸掛式帷幕防滲作用的有限元模擬[J].巖土力學(xué)，2005，26(1)：105-107.

[3] 劉桃溪，辛全才，解曉峰，等.張溝均質(zhì)土石壩滲流和壩坡穩(wěn)定分析 [J].人民黃河，2012,34(3)：120-122.

[4] 陳賢，張立勇，王利容，等，雙河水電站攔河閘壩三維滲流有限元分析[J].人民黃河，2014,36(1)：104-106.

[5] 許季軍，張家發(fā)，程展林.堤防半封閉式防滲墻防滲機(jī)理及設(shè)計(jì)參數(shù)研究[J].人民長(zhǎng)江，2001，32(5)：42-48.

[6] 李少明.防滲墻質(zhì)量缺陷對(duì)土石壩滲流控制的影響[J].南水北調(diào)與水利科技，2012，10(5)：174-177.

The finite element analysis to 3-FEM seepage of a hydroeletric cofferdam

LIU Na,et al.

(SchoolofCivilEngineering,ChangchunInstituteofTechnology,Changchun130012,China)

Abstract：The cofferdam seepage was analyzed by three-dimensional finite element software, and the result shows that the outfall seepage gradient in down-stream of cofferdam is big，where hydraulic failure can happen easily; however the seepage quantity was concentrated in the covering layer under the diaphgragm wall. That was mainly caused by the suspended diaphgragm wall. These two points in the practical engineering must be paid attention to, with the purpose to ensure the stability of the cofferdam.

Key words：seepage; finite element; hydraulic failure; seepage quantity