流域大型水庫群中長期水量調度研究

2022-03-01 07:34:46侯文娟

水利科技與經濟 2022年2期

侯文娟

(新疆伊犁昭蘇縣水利局，新疆伊犁 835600)

1 研究區概況

某河流全長386 km，流域面積約11 600 km2，集水面積約4 892 km2，年平均徑流量35.1×108m3。流域主要以沖擊平原和丘陵為主，流域內分布有大小30多條支流。流域內有大型水庫5座，中型水庫22座，流域水庫總庫容約22.5×108m3，控制集水面積5 000 km2。本文主要針對大型水庫的中長期水量調度問題，因此選擇調蓄能力強的大型水庫進行概化研究[1-2]，并將流域概化為5個水庫(匯水區)和3個用水單位(水資源計算分區)。5個水庫具體情況見表1，流域水庫群概化后示意圖見圖1。

表1 流域大型水庫概況

圖1 流域水庫群概化示意

2 模型構建

2.1 水量調度原則

水量調度不僅要考慮社會經濟發展，還需要兼顧生態環境的可持續發展。為了達到供水高效性、可持續性和水質安全性[3-5]，確定流域水量調度原則如下：①在滿足基本居民生活和生態用水基礎上，兼顧農業和工業用水；②各用水單位必須提交合理的用水計劃，盡可能提高水資源的利用效率；③根據水庫調水性能和用水單元自身特點制定合理的調度方式。

2.2 水量平衡模型及參數

為了將抽象的水量調度問題轉化為數學描述的框架，定義來水點、蓄水點、取水點、退水點以及節點流路，建立流域水量平衡模型，見圖2。

圖2 水量平衡模型

流域用水單元的來水節點主要包括河道和水庫，并將其用于生活用水、生態用水、農業用水和工業用水，最后產生的廢水、棄水再通過凈化后流回河道。因此，要確定調水量和調水方式，就需要確定各類型用水的耗水系數和回水系數[6-9]。本文根據已有研究，綜合確定流域用水參數情況，見表2。

表2 各類用水參數

2.3 模型構建

由于水庫群的水量調度受生活、生態、農業和工業等多方面需求的影響，因此決定了水量調度是一項復雜的綜合性目標決策問題。目標決策的最終目的是使得總的調水目標處于最優狀態，因而首先要建立的是目標函數。目標函數分為主要目標和次要目標，主要目標為供水缺率最小，次要目標為發電效益最大，兩者的計算模型為：

供水缺率：

(1)

(2)

發電效益：

(3)

式中：T為發電效益；N為水庫數量；Hi為i水庫的發電量。

構建好目標函數后，還需要對模型進行約束，約束包括兩部分：一是水庫約束；二是防洪約束。水庫約束包括水量平衡方程、水位約束、出庫流量約束、出力約束和發電水頭約束；防洪安全約束主要為水電站水位邊界。

水量平衡方程：

(4)

式中：Vi,t、Vi,t+1為第i個水庫在計算時段的初庫容和末庫容；Qi,t為平均入庫流量；Qk，t-Tki為下泄流量；Ii,t為t時段第i個水庫的入庫流量。

水位約束：

(5)

出庫流量約束：

(6)

出力約束：

(7)

發電水頭約束：

(8)

2.4 模型計算

逐步優化算法在求解多階段動態優化問題上具有獨特的優勢[10]，相比其它傳統算法能夠避免維數災的問題。因此，本文利用逐步優化算法求取目標函數的最優解，計算過程如下：

1) 確定水庫逐步優化順序，上游和干流水庫優先計算，下游和支流水庫后計算。

2) 利用DP算法得到初始解算結果。

3) 采用POA計算次數N和計算時段t。

4) 運用離散精度向量，離散各t時段的水位過程。

5) 運用動態規劃法，優化t時段相鄰時段的水位過程，t=t+1。

6) 當計算次數達到循環次數N或者達到計算精度要求后停止計算。

3 應用結果

以流域2020年調水情況為例，根據中長期流域的來水預報得到流域的最大可供水量為15.26×108m3，考慮現有供水能力，將最大供水量的90%作為地表的可用水量，其余10%作為河道基流。同時為了全面模擬2020年度水庫群的水量調度情況，還分別設計50%來水頻率(可供水量為20.56×108m3)和95%來水頻率(可供水量為12.92×108m3)兩種情況下的水量調度模擬計算。

3.1 水庫群水位變化過程

模擬得到的不同來水頻率下水庫群調水水位過程曲線(初始水位均為正常蓄水位)見圖3。

圖3 水庫群水位變化趨勢

從圖3中可以看到，位于最上游的1水庫，當處于預報來水可供水量和50%來水頻率可供水量情況下時，水庫水位走勢基本一致；當處于95%來水頻率可供水量情況下時，1-6月份和10-12月份的水庫水位較其他兩種情況時低，7-9月份的水位保持一致，且此時的水位最低，僅為309 m。2水庫位于1水庫下游，當處于預報來水可供水量和95%來水頻率可供水量情況下時，水庫水位走勢基本一致；當處于50%來水頻率可供水量情況下時，1-5月份的水庫水位較其他兩種情況時高，其余月份的水庫水位基本保持一致，且7-10月份的水位最低，為175 m。位于支流1的3水庫，當處于預報來水可供水量和50%來水頻率可供水量情況下時，1-6月份水庫水位走勢基本一致；當處于50%和95%來水頻率可供水量情況下時，9-12月份水庫水位走勢一致，7-8月份時3種情況下的水庫水位最小，均為173 m。位于支流2上游的4水庫，當處于預報來水可供水量和50%來水頻率可供水量情況下時，水庫水位走勢基本一致；當處于95%來水頻率可供水量情況下時，水庫水位略低于其它兩種情況，水庫水位最低時為7-9月份，僅為184.5 m。位于支流2下游的5水庫，當處于預報來水可供水量和95%來水頻率可供水量情況下時，水庫水位走勢基本一致；當處于50%來水頻率可供水量情況下時，水庫水位略高于其它兩種情況，水庫水位最低時為7-9月份，僅為145 m。從水位變化曲線可以看出，每年的7-9月份是調水最多的月份。

3.2 用水保證率

模擬計算分析得到不同來水頻率下各用水單位的用水保證率變化曲線見圖4。

圖4 各用水單位用水保證率

從圖4中可以看到，在預報可供水量(15.26×108m3)情況下時，位于上游的1用水單位和中游的2用水單位的用水保證率隨著月份呈動態變化，具體表現為先減小后增大的變化特征，最低用水保證率均出現在3月份，用水保證率僅為38%和37%，3用水單位的用水保證率橫跨整個2020年度均能保持在100%。在50%來水頻率可供水量(20.56×108m3)情況下，1用水單位和2用水單位的用水保證率呈逐漸升高趨勢，最低用水保證率均出現在1月份，分別為71%和88%，當4月份之后，兩個用水單位的用水保證率可達到100%，而3用水單位的用水保證率均是100%。在95%來水頻率可供水量(12.92×108m3)情況下，1用水單位和2用水單位的用水保證率呈先降低后上升的變化趨勢，最低用水保證率也是出現在3月份，分別為42%和47%，至6月份之后，兩個用水單位的用水保證率也達到100%，而3用水單位的用水保證率仍是100%。綜上分析可以發現，位于中上游的用水單位平均用水保證率要低于下游用水單位的平均用水保證率，這主要是因為3用水地區主要以生活用水為主，而在供水原則設計時，是需要優先保證生活用水的，且用水保證率最低一般出現在開春前后(3-5月份)；預報可供水量和95%來水頻率可供水量情況下的平均用水保證率要低于50%來水頻率可供水量情況下的平均用水保證率。