□馬翠麗□張建□萬占偉

（1河南省水利勘測設(shè)計研究有限公司；2黃河勘測規(guī)劃設(shè)計有限公司）

水庫泥沙淤積及其對工程運用的影響

□馬翠麗1□張建2□萬占偉2

（1河南省水利勘測設(shè)計研究有限公司；2黃河勘測規(guī)劃設(shè)計有限公司）

為分析泥沙淤積對白水水庫庫容及死水位選擇的影響，對水庫來水來沙和淤積形態(tài)進(jìn)行計算分析，設(shè)計水庫淤積形態(tài)，從而確定壩前淤積高程、死水位和水庫淤積后有效庫容。為了工程更安全的興利運用，采用錐體淤積形態(tài)設(shè)計。根據(jù)設(shè)計成果，白水水庫運用50a后，庫區(qū)淤積泥沙187.78萬m3，壩前淤積高程為1 218 m，正常蓄水位下有效庫容約1 094萬m3。綜合考慮泥沙淤積、灌溉效益、庫區(qū)淹沒影響及工程投資等主要因素，確定水庫死水位為1 226 m，為工程規(guī)模論證提供依據(jù)。

水庫泥沙淤積；水庫庫容；死水位

0 引言

由于近年來氣候變化異常，我國西南地區(qū)遭遇歷史罕見的特大旱災(zāi)，城鄉(xiāng)居民生活、農(nóng)業(yè)灌溉和工業(yè)用水均無法保障，急需推進(jìn)水利基礎(chǔ)設(shè)施建設(shè)，新建一批大中小型水庫工程。而水庫要保證一定的使用年限、更好地發(fā)揮效益，工程設(shè)計時需要考慮泥沙淤積引起的一系列問題如庫容損失、淤堵取水口、水質(zhì)惡化、加重維護(hù)成本等。田海濤等對中國內(nèi)地115座水庫淤積情況進(jìn)行統(tǒng)計分析，截止到2003年水庫庫容年均淤積率為0.76%，相當(dāng)于每年損失一座庫容42.30億m3的超大型水庫，其中，中小型水庫比大型水庫淤積嚴(yán)重［3］。龔朝蜀和劉家應(yīng)先后對烏江渡水電站實測泥沙淤積特點與成因進(jìn)行研究，并與原設(shè)計進(jìn)行對比，發(fā)現(xiàn)水庫運行初期壩前出現(xiàn)嚴(yán)重泥沙淤積的原因主要為上游來沙增多、泥沙較細(xì)且河道比降較大易產(chǎn)生異重流、壩前施工圍堰攔淤、水庫初期低水位運行、排沙洞底高程較高實際排沙率很小等，水庫實際為帶狀淤積形態(tài)，提出原設(shè)計的失誤和不足，為水庫工程泥沙設(shè)計總結(jié)了經(jīng)驗和教訓(xùn)。水庫泥沙淤積量和淤積形態(tài)對水庫有效庫容、特征水位、水庫回水長度及取水口位置影響很大。因此，工程規(guī)模論證時必須考慮入庫泥沙及泥沙淤積問題，以保障水庫使用年限，更好發(fā)揮工程效益。

1 白水水庫淤積設(shè)計

1.1 白水水庫基本情況

白水水庫位于北盤江右岸一級支流西泌河上游支流白水河上，壩址位于普安縣江西坡鎮(zhèn)黃泥田村附近，壩址以上集雨面積42.66 km2。工程區(qū)距離普安縣城約11 km，距離江西坡鎮(zhèn)約5 km。

水庫的建設(shè)任務(wù)主要是為解決普安縣城鄉(xiāng)生活、工業(yè)用水以及農(nóng)田和茶園的灌溉用水。設(shè)計最大壩高50.00 m，總庫容1 343萬m3，正常蓄水位1 246.40 m，對應(yīng)庫容1 276萬m3。白水水庫屬Ⅲ等、中型水庫。

1.2 水庫泥沙設(shè)計條件

1.2.1 庫區(qū)河道特征

白水水庫壩址以上集水面積42.66 km2，蓄水區(qū)河道平均比降為10.50‰。正常蓄水位1 246.40 m下，庫區(qū)平均河寬約300 m，壩址上游約1.30 km處河道束窄，河寬不足100 m，從束窄處向上游左右岸均有支溝匯入，河道又展寬，最大河寬達(dá)500 m，平均河寬約300 m。

1.2.2 入庫水沙特征

白水河流域?qū)俚湫偷纳絽^(qū)雨源型河流，徑流由降水產(chǎn)生，徑流特性與降水特性基本一致，徑流年內(nèi)分配不均，洪、枯變化較大。根據(jù)鄰近烏都河流域草坪頭水文站資料分析，5-10月徑流量占全年的85%左右，其中6-9月占全年的70%左右，11月至次年4月占全年的15%左右。年最小流量多出現(xiàn)在每年的3月和4月。

白水河流域?qū)贌o資料地區(qū)，無水文站、雨量站和氣象站。根據(jù)集水面積、下墊面因素及地理位置與設(shè)計流域的關(guān)系，本次設(shè)計選用鄰近流域拖長江草坪頭水文站作為水文計算參證站。水庫壩址徑流采用水文比擬法計算。壩址以上流域多年平均降雨量1 500 mm，多年平均徑流系數(shù)0.51。草坪頭水文站集水面積1 085 km2，控制流域多年平均降水量1 304.4 mm，多年平均徑流量6.94億m3，多年平均徑流系數(shù)0.49。經(jīng)計算，白水水庫多年平均徑流量為3 266萬m3，多年平均流量為1.04 m3/s。

白水河流域無實測泥沙資料。流域內(nèi)開墾度一般，植被較好，但由于流域內(nèi)地勢陡峭，水土流失較為嚴(yán)重。泥沙主要來源于暴雨對坡面的侵蝕以及洪水對河床的沖刷。根據(jù)貴州省懸移質(zhì)輸沙模數(shù)分區(qū)圖，白水水庫所在流域懸移質(zhì)輸沙模數(shù)為500～1 000 t/km2區(qū)域邊緣，輸沙模數(shù)比較大。根據(jù)鄰近流域烏都河草坪頭水文站實測泥沙資料統(tǒng)計，草坪頭水文站多年平均懸移質(zhì)輸沙模數(shù)為945 t/km2，結(jié)合白水河流域地形地貌、下墊面條件，確定白水河流域懸移質(zhì)輸沙模數(shù)為1 000 t/km2，年均懸移質(zhì)輸沙量為4.27萬t；為安全考慮，推移質(zhì)輸沙模數(shù)按懸移質(zhì)輸沙模數(shù)的20%計算，則年均推移質(zhì)輸沙量為0.85萬t。因此，水庫壩址處多年平均年輸沙總量為5.12萬t。

1.3 水庫淤積形態(tài)設(shè)計

1.3.1 泥沙淤積量

白水水庫為蓄水運用，以滿足附近城鎮(zhèn)生活用水和農(nóng)業(yè)灌溉用水等，未設(shè)置排沙洞，難以排泄入庫泥沙，根據(jù)Brune的水庫多年平均攔沙率經(jīng)驗曲線，計算水庫懸移質(zhì)攔沙率為93%，按最大不利設(shè)計，入庫泥沙全部攔淤庫區(qū)。水庫淤積計算年限分別按20a和50a進(jìn)行水庫淤積設(shè)計。懸移質(zhì)泥沙和推移質(zhì)泥沙干容重分別按1.30 t/m3、1.80 t/m3，水庫運行20 a、50 a累計攔沙體積分別為75.11萬m3、187.78萬m3。

1.3.2 水庫淤積形態(tài)及壩前淤積高程

1.3.2.1 縱向淤積形態(tài)

白水水庫無防洪任務(wù)，水庫汛期限制水位為正常蓄水位。水庫泥沙淤積形態(tài)主要受上游徑流變化和水庫運用方式影響。在汛期5-10月，當(dāng)水庫來水大于生活、灌溉和河道生態(tài)用水要求時，水庫蓄水運用，控制水庫水位不超過正常蓄水位；當(dāng)水庫蓄水達(dá)到正常蓄水位且入庫流量大于用水需求時，水庫按來水量下泄，發(fā)生大洪水時水庫敞泄運用；當(dāng)水庫來水小于生活、灌溉和河道生態(tài)用水要求時，水庫按綜合要求下泄水量。在非汛期11-4月份，水庫蓄水興利運用，利用正常蓄水位和死水位之間的調(diào)節(jié)庫容對入庫徑流過程進(jìn)行調(diào)節(jié)，滿足供水目標(biāo)要求。根據(jù)水庫蓄水運用方式可知，白水水庫高水位維持在正常蓄水位運用，除了上游來水較少為了維持下游灌溉和河道生態(tài)等的用水需求才會降低水位下泄水量。

水庫泥沙淤積形態(tài)分別采用清華大學(xué)和陳文彪、謝葆玲經(jīng)驗公式進(jìn)行計算判斷。

清華大學(xué)水庫泥沙淤積形態(tài)經(jīng)驗公式：

式中：a——水庫泥沙淤積形態(tài)判別指數(shù)；V——總庫容，m3；WS——年均懸移質(zhì)入庫沙量，m3；J0——河道比降，?。經(jīng)計算a=3.8，大于2.2，白水水庫為三角洲淤積形態(tài)。

另陳文彪、謝葆玲公式：

當(dāng)φ＞0.04，為三角洲淤積；當(dāng)φ＜0.04，為帶狀淤積。式中為庫水位變化程度，Δh為水庫歷年平均壩前水位變幅，m；h為水庫歷年平均壩前水深，m；根據(jù)水庫正常運行期運用方式表征入庫水沙條件，WS是多年平均年入庫懸移質(zhì)沙量，m3；W為多年平均年入庫水量，m3。經(jīng)計算φ＞0.05，判斷白水水庫為三角洲淤積形態(tài)。

但是，據(jù)調(diào)研貴州省已建和規(guī)劃設(shè)計的水庫縱剖面多為帶狀沿程均勻淤積，最終形成錐體淤積形態(tài)。因此，下面分別按三角洲和錐體淤積設(shè)計，并對設(shè)計成果進(jìn)行比較。

根據(jù)入庫淤積量和庫區(qū)地形、橫斷面測量數(shù)據(jù)，計算庫區(qū)泥沙淤積體各段的淤積厚度和淤積面比降，進(jìn)行鋪沙試算。隨著水庫運行年限增加，錐體形態(tài)與三角洲形態(tài)設(shè)計的壩前淤積高程變化幅度增加，錐體形態(tài)壩前淤積高程較大，對水庫興利運用來說是一種不利的淤積形態(tài)，而且白水河徑流年內(nèi)分配不均，洪枯季變化較大，會影響水庫初期的運行水位。為保證水庫能按規(guī)劃的興利任務(wù)運用，提供更多的興利庫容，并汲取烏江渡水庫淤積設(shè)計教訓(xùn)，白水水庫淤積形態(tài)設(shè)計成果采用錐體淤積形態(tài)。

1.3.2.2 橫向及支流淤積形態(tài)

庫區(qū)泥沙橫向淤積主要集中在天然河道主槽內(nèi)，為水平淤積形態(tài)。庫區(qū)最大支流位于庫區(qū)左岸，距壩1.41 km，河長約2.76 km，比降大于干流，支流自身來沙量很小，泥沙淤積主要來源于干流倒灌，按水平淤積形態(tài)設(shè)計，淤積高程與庫區(qū)入?yún)R口處一致。白水水庫淤積形態(tài)設(shè)計成果如圖1所示。

圖1 水庫運行后泥沙淤積形態(tài)設(shè)計示意圖

1.3.2.3 壩前淤積高程

白水水庫運用20 a后，淤積泥沙75.11萬m3，設(shè)計壩前淤積高程為1 212 m；運用50 a后，淤積泥沙187.78萬m3，設(shè)計壩前淤積高程為1 218 m。

2 泥沙淤積對工程興利運用的影響

2.1 水庫有效庫容

根據(jù)庫區(qū)1：5 000地形圖，量算得到水庫原始庫容曲線。水庫運用后隨著泥沙的累計淤積，有效庫容在不斷減少。根據(jù)水庫運用20 a、50 a后設(shè)計淤積量及設(shè)計淤積形態(tài)進(jìn)行斷面鋪沙，算得水庫運用20 a、50 a后有效庫容。水庫運用20 a和50 a后正常蓄水位對應(yīng)下的水庫有效庫容分別為1 205萬m3和1 094萬m3，庫容損失主要發(fā)生在壩前區(qū)域。

2.2 死水位選擇

白水水庫為供水水庫，不承擔(dān)其他任務(wù)。因此，死水位的選擇應(yīng)考慮水庫淤積、水工建筑物布置，取水要求、淹沒范圍、工程投資等因素綜合確定。在滿足各種要求的前提下，盡量選取較低的水位，以減少淹沒和工程規(guī)模、節(jié)省投資。

根據(jù)水庫淤積形態(tài)分析，水庫運用50 a后壩前淤積面高程為1 218.00 m，根據(jù)建筑物布置，引水鋼管的進(jìn)口底高程為1 222.50 m，應(yīng)低于死水位3.50 m，高于淤積面1.50 m。故初步確定死水位為1 226.00 m。

水庫設(shè)計時死水位的選擇還需考慮其他影響因素：如死水位對自流灌溉面積的影響、庫區(qū)泵站及庫尾橋梁、庫區(qū)淹沒范圍和工程總投資等，確定水庫死水位為1 226.00 m，根據(jù)原始庫容曲線確定對應(yīng)的死庫容為214.40萬m3。

2.3 合理性分析

為分析淤積設(shè)計成果和死水位的合理性及工程安全運行，水庫泥沙按最大不利情況即全部水平淤積到壩前設(shè)計，根據(jù)水庫原始庫容曲線和泥沙淤積量，推算水庫運行50 a壩前水平淤積高程為1 224.00 m。因此，工程選取的死水位是安全合理的。

3 結(jié)論

泥沙淤積是水庫工程設(shè)計時必須考慮的問題，特別是沒有排沙設(shè)施的水庫工程。水庫泥沙淤積形態(tài)與天然河道比降、上游來沙量、泥沙顆粒級配、水庫運用方式及泄洪排沙設(shè)施進(jìn)水口高程等諸多因素密切相關(guān)，難以做出精確計算。白水河屬于少沙河流，本次設(shè)計采用經(jīng)驗公式進(jìn)行的鋪沙試算進(jìn)行淤積設(shè)計。對于多沙河流，應(yīng)采用物理和數(shù)學(xué)模型模擬庫區(qū)泥沙輸移和沉積，結(jié)合經(jīng)驗?zāi)Ｐ秃同F(xiàn)狀工程實際調(diào)查成果，為工程設(shè)計提供更可靠的研究成果。雖然白水水庫所在流域不屬于多沙區(qū)，但水庫屬于中小型水庫，庫容有限。為了保障水庫的使用壽命，發(fā)揮最大的社會效益，建議水庫上游及庫區(qū)周圍加強(qiáng)水土保持工作，減少水土流失。此外，水庫設(shè)計時，經(jīng)濟(jì)投資允許時，增加排沙設(shè)施、攔沙堰等工程措施減少水庫泥沙淤積，從而延長水庫庫容使用年限。此外，我國西南地區(qū)現(xiàn)有水庫測淤資料較少，建議多開展現(xiàn)有水庫泥沙淤積測驗工作，為新建水庫工程提供更可靠的參考數(shù)據(jù)。

［1］劉孝盈，吳保生，于琪洋.水庫淤積影響及對策研究［J］.泥沙研究，2011(6)：37-40.

［2］郜國明，談廣鳴，李濤.多沙河流水庫泥沙研究展望［J］.浙江水利科技，2014(5)：42-46.

［3］謝金明，吳保生，劉孝盈.水庫泥沙淤積管理綜述［J］.泥沙研究，2013(6)：71-80.

TV145

B

1673-8853（2017）04-0077-03

2017-01-05

編輯：劉長垠 邢博輝

馬翠麗，（1979-），女，中級工程師，碩士，從事水文泥沙研究與河流規(guī)劃治理。