程 濤，何萬成，劉書寶

(1.中國電建集團北京勘測設計研究院有限公司，北京 100024；2.國網新源控股有限公司，北京 100761)

0 概 述

在天然河流上修建水電站將會改變河道原來的水沙輸移狀態(tài)，水庫泥沙淤積問題的研究對保障工程運行安全、延長工程使用壽命意義重大。在工程設計中，水庫泥沙淤積多采用一維水沙數(shù)學模型計算，但一維模型僅能給出斷面平均的水沙要素，無法準確模擬地形及水流形態(tài)復雜的近壩區(qū)域泥沙沖淤情況[1-3]。而平面二維水沙模型能較清楚的展示河床的平面變形，因此被廣泛的應用于近壩區(qū)域的河床演變模擬計算[4- 6]。

青藏高原地區(qū)某水電站為日調節(jié)水電站，混合式開發(fā)，正常蓄水位3 206 m，死水位3 203 m，調節(jié)庫容0.16億m3，電站機組包括壩后生態(tài)泄流機組和右岸引水發(fā)電機組。壩址流域上游自然條件較差，植被稀疏，是所在河流的主要產沙區(qū)，水電站庫沙比為8，屬于泥沙問題較嚴重的水庫。

河段輸沙量年內分配不均勻，主要集中在汛期(6月～9月)，占全年輸沙量的97%。其中，7、8兩月占全年的76%。為了延長水庫使用壽命，保證工程的長期正常運行，攔沙、排沙設施的設計和運行管理是本工程需要研究解決的關鍵問題[7- 8]。電站所在區(qū)域已建水庫泥沙淤積研究成果較少，有必要開展水沙模擬研究，為該工程攔排沙設施的設計和最優(yōu)泥沙調度方案的制訂提供參考依據(jù)。為了能夠較全面地模擬地形及水流形態(tài)復雜的近壩區(qū)域泥沙沖淤情況，本文建立平面二維水沙數(shù)學模型研究該區(qū)域水庫泥沙淤積狀況。

1 二維水沙運動模型

1.1 模型選擇

對于地形復雜河段的平面二維水沙運動數(shù)值模擬，多采用正交曲線變換方法，其中一般曲線變換網格生成較靈活，不受計算網格必須嚴格正交的限制。根據(jù)研究電站庫區(qū)河段特性建立一般曲線坐標系平面二維水沙模型，采用控制體積法與SIMPLE 算法對方程進行離散和求解。

二維水沙運動模型主要由水動力模塊和泥沙輸送模塊構成，各模塊基本方程的數(shù)學表達式如下：

(1)

(2)

(3)

q22HSη)]-JαW(S-S*)

(4)

-Jβ(gb-gb*)

(5)

(6)

式中，H為水深，m；u和v為x和y方向的流速，m/s；M=uh，N=vh；Z為水位，m；n為曼寧糙率系數(shù)；D為紊動粘性系數(shù)；ρ為水體密度，g/cm3；S為含沙量；S*為挾沙力；ω為泥沙顆粒沉速，m/s；γ′為泥沙干密度，kg/m3；α為懸移質恢復飽和系數(shù)；β為推移質泥沙恢復飽和系數(shù)；ε為泥沙擴散系數(shù)；ωK、SK、S*K分別對應分組泥沙顆粒沉速、含沙量和挾沙力；gb、gb*分別為單寬推移質輸沙率和有效輸沙率；gbxK、gbyK分別為x和y方向的分組推移質輸移率。

1.2 模型構建

電站壩區(qū)河段平面二維水沙數(shù)學模型上游計算區(qū)域基于壩址前河道地形條件，選為壩前長約4 km的河道，壩址為出口控制斷面，計算區(qū)域網格總數(shù)為14 080。其中，河道主流縱向網格數(shù)為176個，橫向網格數(shù)為80個。利用庫區(qū)的1∶2 000地形資料對網格進行插值，生成的二維水沙數(shù)學模型模擬的地形及樞紐布置示意見圖1。

圖1 數(shù)值模擬地形及樞紐布置示意(單位：m)

根據(jù)壩址實測水位數(shù)據(jù)和上游水文站流量觀測成果率定得到河段糙率為0.04。鑒于該水電站尚未建設，并沒有庫區(qū)淤積的實測資料，不足以根據(jù)泥沙淤積的實測數(shù)據(jù)對上述模型進行驗證。因此，以上游水文站近十年水沙系列作為邊界控制條件，利用模型進行天然河床沖淤計算，挾沙力系數(shù)及指數(shù)分別取值為0.124、1.05。計算結果顯示，河床基本維持天然不沖不淤狀態(tài)，模型適用性較好，可用于電站泥沙淤積計算。

1.3 邊界條件

電站壩址處無實測泥沙資料，以上游水文站實測輸沙率資料為基礎，分析計算水庫入庫沙量。庫區(qū)泥沙淤積計算采用的水沙系列為豐、平、枯代表年水沙過程進行循環(huán)來模擬進口水沙過程，該系列水沙資料循環(huán)排列形成20年系列作為電站庫區(qū)上游計算的水沙資料。懸移質、河床質泥沙顆粒級配采用上游水文站泥沙取樣分析成果。上游進口代表年流量和輸沙率數(shù)據(jù)見表1。

表1 上游進口徑流及輸沙率成果

水電站運行方式：汛期(6月～9月)水庫維持在汛期排沙運行水位運行，泄洪優(yōu)先開啟底孔及生態(tài)泄放孔沖沙，其他時段水庫按日調節(jié)方式運行，水位在正常蓄水位和死水位之間變動。

1.4 計算方案

該電站屬大流量、低水頭電站，能量指標對水頭變化較為敏感。鑒于所在河流輸沙特點，懸移質主要集中在汛期，沙峰與洪峰基本對應，水庫泥沙調度方式應充分利用汛期大流量時降低水位排沙，減少庫區(qū)泥沙淤積，結合電站引水防沙、保持日調節(jié)庫容需求，與上游梯級電站的同步排沙運行考慮，擬定汛期3 202、3 203 m和3 204 m 3個排沙水位方案進行比較。

2 模擬成果

2.1 壩前淤積高程

根據(jù)各方案20年淤積計算結果，電站運行15年后出庫沙量占入庫沙量的比例均達到90%，計算河段基本達到沖淤平衡。電站運行20年后，各方案電站廠房大機組進水口、生態(tài)機組進水口前泥沙淤積剖面見圖2、3，壩前泥沙淤積高程隨排沙水位升高而升高。其中，排沙水位3 204 m方案的壩前淤沙高程接近該水電站廠房進水口高程，排沙水位3 202、3 203 m方案淤沙高程基本滿足電站發(fā)電運行要求。

圖2 廠房大機組進水口前泥沙淤積剖面(垂直河道方向)

圖3 生態(tài)機組進水口前泥沙淤積剖面(平行河道方向)

各方案對電站發(fā)電量、有效庫容以及上游梯級防洪的影響情況見表2，排沙水位愈高，電站多年平均發(fā)電量愈大，有效庫容損失率愈高，水庫回水對上游梯級防洪影響越大。在水庫運行20年后，各方案有效庫容損失率均不超過3%，剩余日調節(jié)庫容均滿足運行要求。在遭遇設計洪水時，各汛期排沙水位方案對上游梯級電站廠房防洪均有不同程度的影響，其中3 202、3 203 m方案影響較小。綜合上述，結合電站引水防沙、盡量保持日調節(jié)庫容，特別是從盡量減少對上游銜接梯級水電站廠房的防洪影響以及梯級整體效益角度考慮，排沙水位不宜高于3 203 m，推薦該水電站排沙水位為3 203 m。

表2 各排沙水位方案運行影響成果

2.2 過機含沙量分析

排沙水位3 203 m方案廠房過機含沙量變化情況見圖4，水庫運行初期入庫泥沙以淤積為主，排沙比相對較小，隨著水庫運行時間的增長，過機含沙量逐漸增大，電站運行15年后庫區(qū)泥沙達到沖淤平衡狀態(tài)，排沙比大于90%，過機含沙量也基本與天然入庫懸移質含沙量一致。

圖4 廠房過機含沙量變化過程

3 結 語

(1)根據(jù)排沙水位比選計算成果，水庫排沙水位愈高，電站多年平均發(fā)電量愈大，有效庫容損失率愈高，水庫回水對上游梯級防洪影響越大。結合電站引水防沙、盡量保持日調節(jié)庫容，特別是從盡量減少對上游銜接梯級電站的防洪影響以及梯級整體效益角度考慮，排沙水位不宜高于3 203 m，推薦該水電站排沙水位為3 203 m。

(2)隨著水庫運行時間的增長，過機含沙量逐漸增大，電站運行15年后庫區(qū)泥沙達到沖淤平衡狀態(tài)，排沙比大于90%，過機含沙量也基本與天然懸移質含沙量一致。

(3)電站非汛期沖沙底孔開啟頻率較低時底孔前易淤積堵塞，建議非汛期沖沙底孔應保持一定開啟頻率或定期進行人工清淤。

總之，數(shù)值模擬結果可為電站運行提供參考；同時，電站運行單位要加強對工程附近河段的水下地形監(jiān)測，以保證電站安全運行。