朱彥澤，盧名燊，劉森宇，程 磊

（1.國網(wǎng)新源集團(tuán)有限公司富春江水力發(fā)電廠，浙江省杭州市 311504；2.武漢大學(xué)水利水電學(xué)院，湖北省武漢市 430072；3.武漢大學(xué)水資源工程與調(diào)度全國重點(diǎn)實(shí)驗(yàn)室，湖北省武漢市 430072）

0 引言

水庫作為一種徑流調(diào)節(jié)的工程手段，在實(shí)現(xiàn)水能資源的合理高效利用上發(fā)揮著重要作用[1]。隨著我國水電行業(yè)的飛速發(fā)展，大批電站水庫相繼建成和投入使用，各大流域逐步形成了梯級水庫群的格局[2，3]?？茖W(xué)開展梯級水庫群聯(lián)合優(yōu)化調(diào)度研究，對于提高水資源利用效率增發(fā)電量、發(fā)揮梯級水庫群的綜合效益最大化具有重要的科學(xué)意義和應(yīng)用價(jià)值[4-6]。

新安江和富春江水庫是錢塘江流域重要的兩座梯級電站，也是華東電網(wǎng)的主力調(diào)峰電廠，擔(dān)負(fù)著華東電網(wǎng)的調(diào)峰、調(diào)頻、事故備用等任務(wù)，對電網(wǎng)的安全、穩(wěn)定運(yùn)行起著重要作用。但在初設(shè)階段，各水庫根據(jù)防洪標(biāo)準(zhǔn)確定設(shè)計(jì)洪水及防洪庫容，并沒有考慮梯級水庫的聯(lián)合調(diào)度問題[7]。此外，傳統(tǒng)的調(diào)度模型已無法滿足水庫群防洪、發(fā)電、灌溉及生態(tài)綜合效益最優(yōu)等需求，以往的調(diào)度方式難以適應(yīng)當(dāng)前復(fù)雜水庫群運(yùn)行管理的實(shí)際需要[8，9]。

鑒于此，本研究以錢塘江流域中上游河段為研究對象，以“新安江—富春江”梯級水電站為調(diào)控主體，根據(jù)錢塘江梯級水電站實(shí)際需求，建立協(xié)調(diào)發(fā)電量及棄水量的梯級水庫聯(lián)合調(diào)度模型，以DPSA 算法求得不同典型年的最優(yōu)運(yùn)行策略，為錢塘江流域梯級水庫調(diào)度運(yùn)行工作提供決策參考。

1 研究區(qū)域及數(shù)據(jù)資料

1.1 梯級水電站資料

新安江水庫位于錢塘江上游，是錢塘江流域最大的控制性樞紐工程，水庫總庫容216.3 億m3，防洪庫容9.47 億m3，是以發(fā)電為主，兼有防洪等綜合利用效益的大（1）型水庫。富春江水庫位于新安江水庫下游的錢塘江干流富春江上，水庫總庫容8.76 億m3，防洪庫容僅0.78 億m3，是以發(fā)電為主，兼有航運(yùn)、灌溉等綜合利用效益的大（2）型水庫，二者擔(dān)負(fù)著華東電網(wǎng)一定的調(diào)峰、調(diào)頻和事故等任務(wù)，各水庫電站的基本參數(shù)如表1 所示。按電站設(shè)計(jì)，富春江水庫為日調(diào)節(jié)水庫，不承擔(dān)防洪、滯洪任務(wù)。富春江壩址洪峰流量的削減，主要取決于上游新安江水庫對蘭江洪水的錯峰調(diào)節(jié)。因此需要通過實(shí)施新安江、富春江水庫聯(lián)合調(diào)度，在汛期洪水來臨前新安江水庫預(yù)騰庫容，為下游富春江水庫洪水錯峰，保障下游河段防洪安全。

表1 錢塘江梯級電站基本參數(shù)Table 1 Basic parameters of qiantang river cascade reservoirs

續(xù)表

1.2 徑流資料

錢塘江梯級水庫徑流資料較齊全，新安江水庫具有1960 ～2022 年共63 年的歷史徑流資料，富春江水庫具有1969 ～2022 年共54 年的歷史徑流資料，通過對長系列年徑流量數(shù)據(jù)進(jìn)行頻率計(jì)算得到的典型年年徑流量具有足夠的代表性。本研究根據(jù)長系列徑流資料進(jìn)行頻率計(jì)算，得到了特豐水年（5%）、豐水年（25%）、平水年（50%）和枯水年（75%）4類典型年的平均流量，分別為461、392、320m3/s 和261m3/s，根據(jù)歷史資料選擇流量相近的2019、1989、1976 年和2006年等日歷年作為典型年。

1.3 計(jì)算條件

（1）調(diào)度期水位上下限。取調(diào)度期為1 年，旬為單位時段，以典型年1 月上旬起算，至12 月下旬結(jié)束，共計(jì)36 個計(jì)算時段。據(jù)新安江水庫原設(shè)計(jì)資料，水庫上游淳安縣移民高程為108.40m，但在實(shí)際運(yùn)行中，為平衡社會各方需求以及保證上游人民生命財(cái)產(chǎn)安全，新安江水庫非汛期水位保持在年消落水位102.00m 附近運(yùn)行，因此本研究中新安江水庫非汛期水位上限取實(shí)際運(yùn)行水位最大值，汛期水位上限取汛限水位，下限取死水位86.00m。富春江水庫為日調(diào)節(jié)水庫，在旬尺度中水位維持正常蓄水位23.00m。

（2）初、末水位。在新安江水庫的起調(diào)水位選擇中，需要綜合考慮發(fā)電量增加和洪水風(fēng)險(xiǎn)之間的平衡，本研究根據(jù)新安江水庫歷史運(yùn)行規(guī)律和華東水新處的建議，針對特豐水年、豐水年、平水年和枯水年不同來水情形，分別設(shè)置起調(diào)水位為98.00m、99.00m、100.00m 和101.00m，末水位則統(tǒng)一取100.00m。富春江水庫初、末水位統(tǒng)一取正常蓄水位23.00m。

2 梯級水庫優(yōu)化調(diào)度模型

2.1 目標(biāo)函數(shù)

梯級水庫優(yōu)化調(diào)度模型有兩個目標(biāo)函數(shù)，分別為梯級發(fā)電量最大和梯級棄水量最小，具體可表示為：

式中：E——梯級水庫總發(fā)電量，億kWh；

S——梯級水庫總棄水量，億m3；

T——計(jì)算期內(nèi)時段總數(shù)；

m——水庫個數(shù)；

Δt——計(jì)算時段；

Ni,t——i水庫t時段的出力；

qi,t——i水庫t時段的棄水流量。

2.2 約束條件

2.2.1 水量平衡約束

式中：Vi,t-1、Vi,t——i水庫t時段初、末的蓄水容積，m3；

Ii,t——t時段水庫天然入庫流量，m3/s；

Qi,t——t時段水庫平均出庫流量，m3/s；

Ei,t——t時段水庫蒸發(fā)、滲漏流量，m3/s。

2.2.2 水位約束

式中：Zi,t+1——i水庫在t+1 時段的水位值，m；

ZLi,t+1——該水庫在t+1 時段的水位下限值，m；

ZUi,t+1——該水庫在t+1 時段的水位上限值，m。

2.2.3 出庫流量約束

式中：QLi,t——水庫放水量下限，m3/s；

QUi,t——水庫放水量上限，m3/s。

2.2.4 水庫之間的水量平衡約束

式中：QJi,t——時段t內(nèi)上水庫i至下水庫i+1 區(qū)間流量，m3/s。

2.2.5 電站出力約束

式中：PLi,t——水電站出力下限，MW；

PUi,t——水電站出力上限，MW。

2.2.6 水庫初始邊界條件

式中：Zi,b——調(diào)度期初水庫蓄水位，m；

Zi,e——調(diào)度期末水庫蓄水位，m。

2.3 優(yōu)化方法

為求解多目標(biāo)優(yōu)化問題，通過權(quán)重法將目標(biāo)函數(shù)式（1）和式（2）合成單個目標(biāo)，即：

式中：α與β——發(fā)電量與棄水量的權(quán)重系數(shù)，考慮新安江和富春江水庫調(diào)度實(shí)際需求，設(shè)置式（8）中的發(fā)電目標(biāo)及棄水目標(biāo)權(quán)重為0.9 和0.1。

梯級水庫聯(lián)合優(yōu)化調(diào)度屬于典型的多階段決策問題，鑒于錢塘江梯級電站中，新安江水庫具有多年調(diào)節(jié)性能，富春江水庫僅有日調(diào)節(jié)功能，該梯級電站可概化為單庫二級梯級水電站發(fā)電系統(tǒng)。逐次逼近動態(tài)規(guī)劃（DPSA）算法是一種有效的降維方法[10]，在精度和效率方面均能滿足求解該梯級水庫調(diào)度問題的要求，故本研究采用DPSA 算法對錢塘江流域梯級水庫的中長期聯(lián)合優(yōu)化調(diào)度模型進(jìn)行求解。DPSA 算法流程見圖1，主要步驟如下。

圖1 DPSA 求解水庫優(yōu)化調(diào)度流程圖Figure 1 Flow chart of optimal reservoir operation based on DPSA

（1）基本資料輸入：輸入的數(shù)據(jù)包括水庫水位、庫容、發(fā)電系數(shù)、保證出力、最大最小流量、水庫的上下游關(guān)系、入流信息等。

（2）迭代計(jì)算：保持其他水庫運(yùn)行狀態(tài)不變的情況下，每次對一個水庫進(jìn)行庫容離散，運(yùn)用DP 進(jìn)行優(yōu)化。

（3）更新水庫優(yōu)化后的庫容、水位、出庫流量等運(yùn)行狀態(tài)。

（4）用DP 繼續(xù)優(yōu)化下一個水庫，直到遍歷完梯級的所有水庫，計(jì)算梯級水庫總發(fā)電量，進(jìn)行下一次迭代計(jì)算。

（5）不斷重復(fù)（2）以后的步驟，直到達(dá)到設(shè)定的迭代次數(shù)，或者結(jié)果滿足精度要求為止，輸出求解結(jié)果。

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

本研究分別針對特豐水年、豐水年、平水年和枯水年不同頻率典型年來水，采用DPSA 算法求解梯級水庫多目標(biāo)聯(lián)合優(yōu)化調(diào)度模型，得到錢塘江梯級水庫中長期優(yōu)化成果，并與常規(guī)調(diào)度方案進(jìn)行對比，如表2 和表3 所示。研究結(jié)果表明：

表2 各典型年梯級水庫優(yōu)化調(diào)度與單庫調(diào)度發(fā)電量對比Table 2 Comparison of optimal scheduling of gradient reservoirs and single reservoir scheduling of power generation in each typical year

表3 各典型年梯級水庫優(yōu)化調(diào)度與單庫調(diào)度棄水量對比Table 3 Comparison of optimal scheduling of gradient reservoirs and single reservoir scheduling of water abandonment for each typical year

（1）多目標(biāo)聯(lián)合優(yōu)化調(diào)度顯著提高了梯級水庫的發(fā)電量，并有效減少了下游水庫泄洪棄水，提高流域水資源利用效率。梯級水庫聯(lián)合優(yōu)化調(diào)度特豐水年、豐水年、平水年和枯水年與常規(guī)梯級調(diào)度相比，增發(fā)電量分別為1.04 億、0.81 億、0.59億和0.57 億kWh，增發(fā)率為2.9%、2.8%、2.7%和3.0%；減少棄水量分別為11.6 億、9.0 億、6.5 億m3和4.9 億m3，減少率為15.6%、14.2%、17.5%和62.0%。

（2）多目標(biāo)聯(lián)合優(yōu)化調(diào)度可以增強(qiáng)梯級水庫的補(bǔ)償效益，更好地發(fā)揮水電系統(tǒng)的整體潛力，獲得比單庫優(yōu)化更大的經(jīng)濟(jì)效益。在特豐水年和豐水年，上游新安江水庫通過在汛期洪水前為富春江水庫預(yù)騰庫容，大幅度減少洪水期富春江水庫因新安江泄洪疊加的水量，明顯增加富春江水庫發(fā)電效益；在平水年和枯水年，新安江水庫通過適當(dāng)抬高發(fā)電水頭，也一定程度上提高了發(fā)電效益。

對梯級水庫在不同典型年下聯(lián)合運(yùn)行和單獨(dú)運(yùn)行的運(yùn)行調(diào)度過程進(jìn)行具體分析?？紤]富春江水庫調(diào)節(jié)能力小，在旬尺度上水位維持不變，故僅重點(diǎn)分析新安江水庫的調(diào)度過程。新安江水庫各典型年調(diào)度過程如圖2 所示。

圖2 新安江水庫特豐、豐、平、枯典型年調(diào)度過程（一）Figure 2 Dispatch process of Xin'anjiang Reservoir in extreme wet，wet，normal and dry typical years（No.1）

圖2 新安江水庫特豐、豐、平、枯典型年調(diào)度過程（二）Figure 2 Dispatch process of Xin'anjiang Reservoir in extreme wet，wet，normal and dry typical years（No.2）

特豐水年來水情景中，新安江水庫全年來水集中在梅汛期，該時段錢塘江暴雨頻現(xiàn)，形成的洪水峰高量大，與蘭江洪水疊加后在富春江水庫壩址形成特大洪峰，給下游帶來巨大的防汛壓力。為減輕下游富春江水庫防汛壓力，在特豐水年來水情景時，新安江水庫的調(diào)度策略為：①汛前，以水定電，保證水庫水位處于保證出力區(qū)內(nèi)運(yùn)行，并根據(jù)主汛期預(yù)報(bào)水量進(jìn)行不同程度預(yù)泄；②汛期開始至梅汛期前，根據(jù)下游富春江斷面的洪水預(yù)報(bào)情況控制泄流，有洪水時減少泄流為富春江水庫錯峰，無洪水時大發(fā)滿發(fā)，在入梅前將水位預(yù)泄至102.00m；③梅汛期，保證不超過自身汛限水位（106.50m）的同時盡可能減少富春江水庫的泄洪棄水；④汛后，根據(jù)來年豐枯預(yù)報(bào)結(jié)果以及年度實(shí)際運(yùn)行調(diào)度情況將水位降至目標(biāo)水位。

豐水年和平水年來水情景中，新安江水庫在主汛期內(nèi)中小洪水頻發(fā)，為應(yīng)對多場次洪水，在豐水年和平水年情景時，新安江水庫的調(diào)度策略為：①汛前，按照正常發(fā)電調(diào)度方式運(yùn)行，汛前水位逐步平穩(wěn)消落至目標(biāo)水位，保證進(jìn)入汛期前有足夠防洪庫容；②主汛期，遭遇洪水時水庫下泄流量以不超過汛限水位為目標(biāo)，并同時減少下游富春江水庫棄水，無洪水時水庫下泄流量按下游富春江水庫不棄水為目標(biāo)，并盡可能為下一場洪水騰出防洪庫容；③汛后，平緩地將水位消落至目標(biāo)水位。

枯水年來水情景中，新安江水庫除春汛和梅汛期外，其余大部分時間流量較低，需要通過水庫調(diào)度以確保供水和發(fā)電的可持續(xù)性，在枯水年情景時，新安江水庫的調(diào)度策略為：①非汛期，減少泄流提高發(fā)電水頭，充分提高發(fā)電效益，但水位不超過上基本調(diào)度線；②汛期，根據(jù)降雨預(yù)報(bào)和河流流量監(jiān)測，及時調(diào)整水庫的蓄水量，在保證自身防洪安全的前提下為富春江水庫錯峰，洪末攔蓄水位，增發(fā)洪水電量。富春江水庫為大流域小水庫，調(diào)控能力有限，因此錢塘江流域梯級水庫中富春江水庫不同典型年的調(diào)度策略為：①非洪水期，以水定電，保持水庫水位處于正常蓄水位（23.00m）運(yùn)行，提高發(fā)電效率；②洪水期，與新安江水庫實(shí)行聯(lián)合調(diào)度，同時在保證防洪安全的前提下，利用可調(diào)庫容在洪前預(yù)泄騰庫，洪末攔蓄水位，增發(fā)洪水電量；在洪峰抵達(dá)前降低水位，發(fā)揮滯洪削峰作用，減輕下游防汛壓力，實(shí)現(xiàn)防洪減災(zāi)的社會效益。通過富春江水庫與新安江水庫實(shí)施聯(lián)合調(diào)度，充分發(fā)揮了上游新安江水庫的調(diào)蓄作用，保障富春江干流桐廬以下河段的防洪安全，并有效提高梯級電站的發(fā)電效益。

4 結(jié)論

本文以錢塘江流域梯級水庫電站為研究對象，對新安江—富春江梯級水庫進(jìn)行中長期優(yōu)化調(diào)度研究，建立協(xié)調(diào)發(fā)電量及棄水量的梯級水庫聯(lián)合調(diào)度模型，采用DPSA 算法進(jìn)行求解，得到以下結(jié)論：

（1）多目標(biāo)聯(lián)合優(yōu)化調(diào)度顯著提高了梯級水庫的發(fā)電效益。聯(lián)合優(yōu)化調(diào)度在不同典型年相較常規(guī)調(diào)度增發(fā)電量0.57 億～1.04 億kWh，減少棄水4.9 億～11.6 億m3，平均提高發(fā)電量2.85%，經(jīng)濟(jì)效益顯著。

（2）多目標(biāo)聯(lián)合優(yōu)化調(diào)度增強(qiáng)梯級水庫的補(bǔ)償效益。通過實(shí)施梯級水庫聯(lián)合調(diào)度，充分發(fā)揮了上游新安江水庫的調(diào)蓄作用，汛期前為下游富春江水庫預(yù)騰庫容，大幅度減少洪水期間富春江水庫因新安江泄洪疊加的水量，有效提高梯級電站的發(fā)電效益和防洪安全。

（3）多目標(biāo)聯(lián)合優(yōu)化調(diào)度為新安江水庫和富春江水庫在豐平枯典型年下提供了相應(yīng)的調(diào)度策略。調(diào)度結(jié)果建議特豐水年、豐水年和平水年新安江水庫在汛前以水定電，保持常規(guī)水位運(yùn)行，并在汛期前根據(jù)當(dāng)年來水預(yù)報(bào)提前預(yù)泄，保證汛期防洪庫容，汛后則根據(jù)來年豐枯預(yù)報(bào)結(jié)果以及實(shí)際運(yùn)行調(diào)度情況將水位降至目標(biāo)水位，枯水年則需要適當(dāng)減少泄流提高發(fā)電水頭以確保供水和發(fā)電的可持續(xù)性。