楊 會 剛

(大唐四川集控中心，四川 成都 610091)

1 概 述

長河壩、黃金坪梯級電站的送出線路均是通過500 kV線路送至康定變電站，同屬四川電網(wǎng)“康甘斷面”。目前，“康甘斷面”的水電站總裝機(jī)容量已近8 000 MW，其中長河壩、黃金坪梯級電站裝機(jī)容量3 450 MW，占“康甘斷面”總裝機(jī)容量的43%。而“康甘斷面”送出能力為4 800 MW，汛期兩級電站將面臨閑置50%裝機(jī)容量，產(chǎn)生大量棄水。

長河壩、黃金坪電站作為梯級電站，上下游水力聯(lián)系緊密，但兩級電站的運行方式仍是獨立運行，沒有形成有效的聯(lián)動機(jī)制，負(fù)荷不匹配給兩級電站的水庫調(diào)度工作帶來了極大的困難。以兩級電站2018年1～3月的運行情況為例，負(fù)荷批復(fù)不匹配情況頻發(fā)。兩級電站日發(fā)電水量差在500萬m3(黃金坪可調(diào)庫容的30%)以上的情況多達(dá)16 d。其中，按照1月19日、3月6日的批復(fù)負(fù)荷，黃金坪在死水位附近日發(fā)電水量比長河壩多1 100萬m3(黃金坪可調(diào)庫容的70%)，水量入不敷出，水位持續(xù)走低，耗水率偏高。

大渡河干流規(guī)劃雙江口水電站為上游控制性水庫，電站裝機(jī)容量2 000 MW，水庫正常蓄水位2 500 m，庫容為28.97億m3，調(diào)節(jié)庫容19.17億m3，具有年調(diào)節(jié)能力。設(shè)計的水庫調(diào)度原則是12月～次年4月為供水期，庫水位于4月底消落至死水位。目前，雙江口電站處于建設(shè)期，計劃于2024年投產(chǎn)。

長河壩、黃金坪梯級電站設(shè)計與雙江口電站聯(lián)合運行，長河壩電站水位在1 680～1690 m之間運行，發(fā)電耗水率為1.8～1.9 m3/kWh。目前，長河壩電站枯水期來水偏少，發(fā)電耗水率維持在1.8～2.3 m3/kWh，其經(jīng)濟(jì)性有待提高。

2 聯(lián)動機(jī)制和水位控制

2.1 出力聯(lián)動機(jī)制

基于長河壩、黃金坪梯級電站工程和機(jī)組特性，結(jié)合NHQ曲線和多年運行數(shù)據(jù)，通過水量平衡計算，建立兩級電站在不同水位下的出力匹配關(guān)系，見表1。

2.2 水位控制研究

表1 長河壩、黃金坪梯級電站不同水位出力比例表

2.2.1 水位控制方案

長河壩電站設(shè)計的運行方式為：全年在死水位和正常蓄水位之間做周調(diào)節(jié)運行，即水位在1 680～1 690 m間運行；當(dāng)遭遇特枯水年時，可降至極限死水位1 650 m，做季調(diào)節(jié)運行。

目前，長河壩電站的運行條件較設(shè)計工況有兩個較大的差異，具體如下：

(1)“康甘斷面”最大送出能力為3 200 MW，長河壩電站按照斷面裝機(jī)容量的比例計算，平均最大出力為1 100 MW。

(2)由于雙江口水庫尚未投運，長河壩電站枯水期12月～次年4月平均入庫流量較雙江口投運后偏少151 m3/s左右。

由于運行條件發(fā)生重大變化，長河壩電站的運行方式需結(jié)合實際情況進(jìn)行研究。本文對長河壩電站1952～2017年的長系列徑流資料進(jìn)行逐月排頻，形成25%、50%、75%來水頻率和多年平均共計4個典型徑流序列，按照“康甘斷面”受限和實際發(fā)電能力兩種模式，分別計算汛前將長河壩水位均勻消落至1 650 m和1 680 m兩種方案的年發(fā)電量情況，結(jié)果見表2。

“康甘斷面”受限時，最大出力1 100 MW，庫水位按1 650 m控制較經(jīng)濟(jì)。其中，枯水年(75%)較1 680 m方案增加發(fā)電量0.46億kWh，增幅0.63%；中水年(50%)增加0.74億kWh，增幅0.98%；豐水年(25%)增加0.88億kW.h，增幅1.12%；多年平均來水增加發(fā)電量0.79億kWh，增幅1.03%。來水越豐，年發(fā)電量增加越明顯。

表2 長河壩電站不同水位控制方案的發(fā)電量計算表 /億kWh

“康甘斷面”不受限時，最大出力2 600 MW，庫水位按1 680 m控制較經(jīng)濟(jì)。其中，枯水年(75%)較1 650 m方案增加發(fā)電量3.4億kWh，增幅3.02%；中水年(50%)增加2.65億kWh，增幅2.11%；豐水年(25%)增加1.01億kWh，增幅0.74%；多年平均來水增加發(fā)電量1.72億kWh，增幅1.32%。來水越枯，年發(fā)電量增加越明顯。

2.2.2 水位控制結(jié)論

在“康甘斷面”送出受限沒有解決前，長河壩電站水位按1 650 m控制較為經(jīng)濟(jì)。其中，長河壩電站4月份的平均入庫流量較為關(guān)鍵，當(dāng)4月平均來水小于518 m3/s時，1 650 m方案較1 680 m方案發(fā)電量增加顯著。

最優(yōu)的1 650 m運行方式為：汛期6～9月在1 680 m附近運行，10月初回蓄至1 690 m，10月～次年3月維持在1 690 m運行，4月均勻消落至1 650 m，5月來水大于發(fā)電時回蓄至1 680 m。

3 優(yōu)化調(diào)度模型研究

目前，求解梯級水電站優(yōu)化運行的方法主要有線性規(guī)法、逐步優(yōu)化算法(POA)、動態(tài)規(guī)劃(DP)和遺傳算法(GA)等。由于長河壩、黃金坪梯級電站的運行與上游雙江口、猴子巖電站的運行方式密不可分，優(yōu)化模型需向上游進(jìn)行延伸計算，采用動態(tài)規(guī)劃模型易陷入“維數(shù)災(zāi)”，遺傳算法收斂速度慢[1-2]。本文采用逐步優(yōu)化算法(POA)建立長河壩、黃金坪梯級電站優(yōu)化調(diào)度模型，通過賦予合理的初始決策條件，加快模型迭代收斂速度，實現(xiàn)最優(yōu)化[3-4]。

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

長河壩、黃金坪梯級電站以年發(fā)電量最大為優(yōu)化調(diào)度目標(biāo)，目標(biāo)函數(shù)如下：

(1)

式中E為梯級電站最大化的年發(fā)電量；Mt為第t時段的小時數(shù)；n為電站個數(shù)；T為年時段數(shù)；η為出力系數(shù)；Qi,t為第t時段第i個水電站的發(fā)電流量；Hi,t為第t時段第i水電站的水頭。

3.2 約束條件

(1)水量平衡約束。

Vi,t+1=Vi,t+(qi,t-Qi,t-Si,t)Δt?t∈T

(2)

式中Vi,t+1為第t時段末第i個水庫的蓄水量；Vi,t為第t時段初第i個水庫的蓄水量；qi,t為第t時段第i個水庫的入庫流量；Qi,t為第t時段第i個水電站的發(fā)電流量；Si,t為第t時段第i個水電站的棄水流量；Δt為計算時段長度。

(2)運行水位約束。

ZiminZi,tZimax?t∈T

(3)

式中Zimin為第i個水庫最低運行水位；Zi,t為第t時段第i水庫的水位；Zimax為第i個水庫最高運行水位。

(3)電站出力約束。

Piminη·Qi,t·Hi,tPimax?t∈T

(4)

式中Pimin、Pimax分別為第i個電站的第t時段的最小出力和最大出力。

(4)出庫流量約束。

QiminQi,tQi,tmax?t∈T

(5)

式中Qimin為第i電站最小下泄流量要求；Qi,tmax為第i電站在第t時段最大的下泄能力，包括最大發(fā)電流量和最大泄洪流量。

(5)非負(fù)約束。各變量均為非負(fù)值。

3.3 模型求解

逐步優(yōu)化算法將多階段的問題分解為多個兩階段問題，解決兩階段問題的方法是：固定當(dāng)前其他階段的變量，對所選兩階段的決策變量進(jìn)行搜索尋優(yōu)，在解決當(dāng)前兩階段問題后再計算下一個兩階段，將上次的結(jié)果作為下次優(yōu)化的初始條件，進(jìn)行尋優(yōu)，反復(fù)循環(huán)，直到收斂為止[5-6]。假設(shè)梯級電站的調(diào)度期為1年，初始時刻為5月初，終止時刻為次年4月末，將一年離散為T時段，梯級電站數(shù)為n，電站序號為i(0≤i≤n-1)，則POA算法的計算步驟如下：

(1)確定初始條件和約束條件。

(2)按照電站順序，依次對第i個電站尋優(yōu)。固定第0時刻和第2時刻的水位Zi,0和Zi,2不變，調(diào)整第1時刻的水位Zi,1，使第0和1兩時段的發(fā)電量最大。

(3)按照電站順序，依次對第i個電站下一時刻進(jìn)行尋優(yōu)。固定第1時刻和第3時刻的水位Zi,1和Zi,3不變，調(diào)整第2時刻的水位Zi,2，使第1和2兩時段的發(fā)電量最大。

(4)重復(fù)步驟(3)，直到終止時刻為止。從而得到初始條件和約束條件下的梯級各水庫水位過程、發(fā)電流量過程和梯級總電量。

(5)以前次求得的各水庫過程線為初始軌跡，重新回到步驟(2)，直到相鄰兩次迭代求得的發(fā)電量增量達(dá)到精度要求為止。

3.4 模型計算

為充分發(fā)揮水庫群的綜合利用效益，實現(xiàn)梯級水庫的分?jǐn)偡篮椤㈠e峰蓄水的聯(lián)合調(diào)度功能，在保證水庫群汛期防洪的前提下，按照后汛期從上到下依次蓄水，枯水期按照從上到下依次消落的原則進(jìn)行模型計算[7]如下：

(1)水位控制原則。雙江口：5～9月控制在2 480 m附近運行，9月中旬蓄至2 500 m，10～11月維持2 500 m運行，12月率先開始均勻消落，4月中旬消落至2 420 m。猴子巖：5～9月控制在1 835 m附近運行，9月下旬蓄至1 842 m，10月～次年3月維持1 842 m運行，4月底消落至1 835 m。長河壩：5～9月控制在1 680 m附近運行，9月底蓄至1 690 m，10月～次年4月維持1 690 m運行，4月底消落至1 680 m。黃金坪：配合長河壩的運行方式在1 472～1 476 m之間運行。

(2)最小下泄流量控制。為保證下游用水需求，雙江口按不低于121 m3/s下泄，猴子巖按不低于160 m3/s下泄，長河壩按不低于166.5 m3/s下泄，黃金坪按不低于168 m3/s下泄。

將各電站的水位庫容關(guān)系曲線、水位流量關(guān)系曲線、泄洪能力曲線作為模型的基本計算參數(shù)，并對黃金坪電站1952～2017年的來水資料進(jìn)行還原計算，消除猴子巖和長河壩電站對天然徑流的調(diào)節(jié)影響，從而推算出雙江口電站的逐年來水過程，作為初始徑流序列，輸入模型進(jìn)行計算。各電站歷年的計算過程均滿足水位、流量、出力等各項約束，計算結(jié)果合理?！叭龓焖募墶彪娬径嗄昶骄鶃硭膬?yōu)化調(diào)度計算成果見表3。

3.5 模型驗證

表3 “三庫四級”電站多年平均來水優(yōu)化調(diào)度計算成果表

為檢驗計算成果的合理性，本文以神經(jīng)網(wǎng)絡(luò)模型對逐步優(yōu)化模型的計算成果進(jìn)行校核。人工神經(jīng)網(wǎng)絡(luò)是人腦生理機(jī)制的模擬，善于處理不完善的聯(lián)想記憶、特征缺損模式識別、規(guī)則自動學(xué)習(xí)等問題，在水庫群優(yōu)化調(diào)度中有著廣泛的應(yīng)用[8]。

以雙江口入庫流量、雙江口時段初水位、雙江口-猴子巖區(qū)間流量、猴子巖時段初水位、猴子巖-長河壩區(qū)間流量和長河壩時段初水位為模型輸入；以長河壩時段末水位為模型輸出，建立三層BP網(wǎng)絡(luò)模型，見圖1。

即：f(x)=1/[1+E(x)]。

輸出層各節(jié)點的輸出為：

(6)

圖1 三層BP網(wǎng)絡(luò)結(jié)構(gòu)圖

模型輸出層僅一個節(jié)點，節(jié)點輸出為：

(7)

輸出層輸出值與期望輸出的誤差為：

(8)

將歷年徑流資料、各電站工程曲線以及約束條件輸入人工神經(jīng)網(wǎng)絡(luò)模型進(jìn)行迭代計算，當(dāng)輸出層輸出值與期望輸出之間的誤差最小時，模型達(dá)到最優(yōu)。將人工神經(jīng)網(wǎng)絡(luò)模型對長河壩、黃金坪梯級電站1952～2017年的發(fā)電量計算結(jié)果與逐步優(yōu)化模型計算結(jié)果進(jìn)行對比(表4)，逐步優(yōu)化模型的最優(yōu)性和收斂性結(jié)果如下：

表4 長河壩、黃金坪梯級電站逐步優(yōu)化模型與人工神經(jīng)網(wǎng)絡(luò)模型年發(fā)電量計算成果對比表 /億kWh

(1)最優(yōu)性。兩級電站的逐步優(yōu)化模型1952～2017年平均發(fā)電量為154.15億kWh，較神經(jīng)網(wǎng)絡(luò)模型多0.61億kWh，偏多0.4%；較兩級電站設(shè)計發(fā)電量149.11 kWh多5.04億kWh，偏多3.4%，逐步優(yōu)化模型的計算結(jié)果最優(yōu)性表現(xiàn)突出。

(2)收斂性。在1952～2017年共計66年的系列優(yōu)化計算中，逐步優(yōu)化模型在94%的年份發(fā)電量計算收斂性最優(yōu)。僅有1959、1973、2002、2006、2007年五個特枯年(來水頻率90%以上)發(fā)電量計算結(jié)果，較人工神經(jīng)網(wǎng)絡(luò)模型偏少0.2～0.9%，整體收斂性較好。

總體來說，從本次對長河壩、黃金坪梯級電站1952～2017年的發(fā)電量計算結(jié)果來看，逐步優(yōu)化模型計算較穩(wěn)定，收斂性較好。

4 結(jié) 語

基于現(xiàn)實需要，開展長河壩、黃金坪梯級水庫優(yōu)化調(diào)度研究。筆者從梯級電站的工程特性入手，通過探索出力聯(lián)動機(jī)制，基于逐步優(yōu)化算法(POA)建立優(yōu)化調(diào)度模型，以1952～2017年的徑流資料為基礎(chǔ)，對梯級電站的年度發(fā)電量進(jìn)行優(yōu)化計算，并通過BP人工神經(jīng)網(wǎng)絡(luò)模型對計算結(jié)果進(jìn)行檢驗。為兩級電站的科學(xué)調(diào)度了提供了參考，為流域水庫群聯(lián)合調(diào)度提供了思路。針對優(yōu)化調(diào)度存在的問題，建議如下：

優(yōu)化電網(wǎng)運行方式，加快外送通道建設(shè)。一方面制定合理的特高壓直流外送方式，減少對康

甘、攀西、九石雅斷面的輸電能力限制；另一方面盡快啟動康定至蜀都改接和串補工程，加快省內(nèi)特高壓交流環(huán)網(wǎng)和跨省直流通道建設(shè)，全力做好跨省跨區(qū)水電外送工作。

加快龍頭水庫建設(shè)，統(tǒng)籌流域聯(lián)合調(diào)度。加快推進(jìn)雙江口水庫建設(shè)，研究防洪、發(fā)電效益評價指標(biāo)，統(tǒng)籌水庫群和單體水庫的利益關(guān)系，建立公平合理的聯(lián)合調(diào)度補償機(jī)制。

建立流域水情信息、水文預(yù)報、氣象預(yù)報共享機(jī)制，提高來水預(yù)報精度，為流域水庫群聯(lián)合調(diào)度提供技術(shù)支持。

結(jié)合最新的水情自動測報以及氣象預(yù)報技術(shù)，研究水庫的中小洪水調(diào)度模式以及汛期水位動態(tài)控制方案。