盧 薇 ，王 娟 ，劉曉慶

（1.中國水電顧問集團(tuán)成都勘測(cè)設(shè)計(jì)研究院，四川 成都 610072；2.上海勘測(cè)設(shè)計(jì)研究院，上海 200434；3.西華大學(xué)能源與環(huán)境學(xué)院，四川 成都 610039）

1 工程概況

某水電站位于四川省阿壩藏族羌族自治州馬爾康縣、金川縣境內(nèi)大渡河上游東源（主源）足木足河,與西源（次源）綽斯甲河匯合口以下1～6 km河段，壩址上距馬爾康縣城約45 km，下距金川縣城約43 km，經(jīng)馬爾康至成都的公路距離約400 km，有國道G317線、省道S211線通達(dá)。

該水電站是大渡河流域水電梯級(jí)開發(fā)的上游控制性水庫，是大渡河流域水電梯級(jí)開發(fā)的關(guān)鍵性工程之一。壩址處控制流域面積39330 km2，年徑流量166億m3，多年平均流量527 m3/s。壩址區(qū)河谷屬高山深切曲流河谷，兩岸山體雄厚，河谷深切，谷坡陡峻；出露地層巖性主要為∶燕山早期木足渡似斑狀黑云鉀長花崗巖和晚期可爾因二云二長花崗巖；河床覆蓋層深厚（最厚約67.8 m）；壩址區(qū)無區(qū)域性斷裂切割。電站水庫正常蓄水位2500 m，對(duì)應(yīng)庫容約27.32億m3，具有年調(diào)節(jié)能力，電站裝機(jī)容量2000 MW，年發(fā)電量81.28億kW·h，設(shè)計(jì)枯水年枯期平均出力48.6萬kW。土質(zhì)心墻堆石壩最大壩高約312 m，樞紐最大泄洪流量約8000 m3/s。

攔河大壩采用礫石土心墻堆石壩，礫石土心墻堆石壩，壩頂高程2508.00 m，壩頂長度約642 m。河床部位心墻底高程2202.00 m，最大壩高312 m（含6 m厚基座），壩頂寬度16.00 m。上游壩坡為1∶2.0，下游壩坡1∶1.90。心墻頂高程2508.00 m，寬4.00 m，與防浪墻進(jìn)行可靠連接，并高于最高靜水位（2504.51 m）；心墻上、下游坡初選均為1∶0.2，心墻底部與新鮮基巖相連，最大斷面心墻底高程2202.00 m，順河向?qū)挾葹?28.00 m。

2 有限元模型的建立

計(jì)算模型所選剖面為大壩的最大橫剖面，如圖1所示，計(jì)算模型包括上下游堆石體、上下游過濾層、上下游反慮層、心墻、上下游覆蓋層9個(gè)材料區(qū)，計(jì)算模型坐標(biāo)分別為X軸順河指向下游為正，Y軸沿壩高鉛直向上為正，計(jì)算選取范圍分別為：水平方向自距上游壩坡300 m處至距下游壩坡300 m處，鉛直方向自壩頂至基巖面；邊界條件：上下游坡面為水頭邊界，AB，CD為默認(rèn)不透水邊界，BC為基巖不透水邊界；其剖分為8334個(gè)節(jié)點(diǎn)，9648個(gè)單元。

2.1 計(jì)算參數(shù)

水庫特征水位：正常蓄水位為304 m，死水位為224 m，下游水位55.21 m，壩體材料參數(shù)見表1。

2.2 計(jì)算工況的選取

在正常蓄水位下，固定上下游水位，采用飽和與非飽和的方法，考慮材料滲透系數(shù)的各向異性來進(jìn)行滲流分析，工況組合如下：

工況 1：Kx=Ky

工況 2：Kx=2Ky

工況 3：Kx=3Ky

工況 4：Kx=4Ky

工況 5：Kx=5Ky

表1 壩體材料參數(shù)表

3 滲流計(jì)算結(jié)果分析

以Kx=Ky（工況 1）的為代表，如圖1、圖 2表示了飽和與非飽和滲流計(jì)算浸潤線、流速矢量圖。

圖1 （Kx=Ky）非飽和滲流計(jì)算的浸潤線、流速矢量圖

圖2 （Kx=Ky）飽和滲流計(jì)算的浸潤線、流速矢量圖

由圖1和圖2的對(duì)比，在考慮非飽和區(qū)的穩(wěn)定滲流場(chǎng)中，浸潤線（零壓線）附近的飽和區(qū)的滲流水，會(huì)穿過浸潤線而進(jìn)入非飽和區(qū)，而在只考慮飽和區(qū)的穩(wěn)定滲流分析中，通常將浸潤線當(dāng)作滲流場(chǎng)最上部的一根流線，沿著此線沒有流量的流入和流出。

在考慮飽和與非飽和兩種計(jì)算方法的情況下，對(duì)所得的等壓線分布圖（圖略）進(jìn)行比較可見，由于心墻和上游壩殼的滲透系數(shù)相差較大，水頭相當(dāng)于直接作用在心墻上，心墻上游的等壓線幾乎呈水平線分布。在非飽和穩(wěn)定滲流中，由于非飽和區(qū)的存在，心墻下游壩殼也存在水壓，等壓線分布于整個(gè)壩體，增加水平方向的滲透系數(shù)，更多的水流向下游壩殼，負(fù)壓區(qū)逐漸減小,負(fù)等壓線向下游推移。而飽和穩(wěn)定滲流中，由于無水流入非飽和區(qū)，心墻下游壩殼部分無等壓線分布，隨著水平滲透系數(shù)的增加，整個(gè)壩體的等壓線分布無明顯變化。

此外，表2顯示了不同計(jì)算工況時(shí)單寬滲流量、逸出點(diǎn)Y坐標(biāo)、逸出坡降。

從表2可以看出，不考慮非飽和區(qū)影響，Kx=Ky時(shí)，單寬流量、心墻下游溢出點(diǎn)坐標(biāo)、逸出點(diǎn)坡降分別為 3.1912 m3/d、56.214 m、1.428，隨著水平滲透系數(shù)的增加，當(dāng)Kx=5Ky時(shí)，滲流量為17.613 m3/d，增大較明顯，心墻下游邊坡逸出點(diǎn)略有增大，為57.179 m，浸潤線有所抬升；逸出坡降略有減小，為1.271。當(dāng)考慮非飽和區(qū)的影響時(shí)，在各個(gè)對(duì)應(yīng)工況下，其單寬流量、逸出點(diǎn)高程大于不考慮非飽和區(qū)影響時(shí)的值，滲透坡降小于不考慮非飽和區(qū)影響時(shí)的值，說明材料的各向異性在考慮非飽和區(qū)影響的條件下對(duì)單寬流量影響較大。

表2 不同工況時(shí)單寬滲流量、逸出點(diǎn)Y坐標(biāo)、逸出坡降對(duì)照表

4 結(jié)語

1）考慮非飽和區(qū)的影響和不考慮非飽和區(qū)影響，采用兩種方法計(jì)算所得等壓線分布分布存在較大差異，由于非飽和區(qū)的存在，流水會(huì)穿過浸潤線而進(jìn)入非飽和區(qū)，下游壩殼區(qū)域也存在滲流。前者計(jì)算出的單寬流量、心墻下游坡逸出高程比后者大，而坡降小于后者，而實(shí)際工程中土體都存在非飽和區(qū)，所以用非飽和理論來進(jìn)行滲流分析計(jì)算更為貼近工程實(shí)際。

2）考慮滲透參數(shù)的各向異性進(jìn)行滲流分析，隨著水平滲透系數(shù)的增加，滲透坡降變化不大，均在允許坡降范圍內(nèi)。單寬流量增大較明顯，對(duì)其影響較大，這需要根據(jù)興利的要求來選取排水設(shè)施，同時(shí)也應(yīng)做好防滲措施。實(shí)際工程中，存在著土體不均勻性、材料各向異性的特點(diǎn)，所以對(duì)材料滲透參數(shù)的各向異性來進(jìn)行滲流分析較為合理。