徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng)研發(fā)

肖　楊,鄺錄章

(五凌電力有限公司,湖南長沙410004)

?

徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng)研發(fā)

肖楊,鄺錄章

(五凌電力有限公司,湖南長沙410004)

摘要：為滿足水電廠實(shí)時調(diào)度需要,構(gòu)建了以發(fā)電量最大為目標(biāo)的徑流式水電廠實(shí)時優(yōu)化調(diào)度模型,采用離散動態(tài)規(guī)劃法逐步逼近尋優(yōu)計(jì)算,開發(fā)出徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng)。該系統(tǒng)可預(yù)報(bào)計(jì)算未來24 h的入庫流量過程,并可實(shí)時滾動計(jì)算滿足約束條件的未來24 h內(nèi)各時段最優(yōu)發(fā)電運(yùn)行過程,從而量化指導(dǎo)水電廠實(shí)時優(yōu)化調(diào)度。通過對6座徑流式水電廠的應(yīng)用,入庫流量計(jì)算準(zhǔn)確率達(dá)98%以上,平均每年可優(yōu)化增加發(fā)電量4 500萬kW·h。

關(guān)鍵詞：徑流式水電廠；實(shí)時優(yōu)化調(diào)度；離散動態(tài)規(guī)劃法

1工程概況

由五凌電力有限公司(以下簡稱“五凌電力”)開發(fā)建設(shè)、運(yùn)行管理的沅水、資江、湘江上12座水電站已全部接入遠(yuǎn)程集中控制運(yùn)行,如何提高水能利用效率,創(chuàng)新水庫調(diào)度工作,創(chuàng)造最大的發(fā)電效益,成為水電廠運(yùn)行管理及水庫調(diào)度研究的重點(diǎn)。目前,盡量抬升運(yùn)行庫水位、降低發(fā)電耗水率的水庫運(yùn)行思維已成為運(yùn)行調(diào)度人員的共識,但如何具體量化調(diào)度期內(nèi)各時段的庫水位及出力值,使得發(fā)電效益最大,亟需在實(shí)際生產(chǎn)運(yùn)行中深入研究。

水庫優(yōu)化調(diào)度理論研究已較為成熟,但須將理論成果轉(zhuǎn)化為生產(chǎn)實(shí)踐[1]。本文結(jié)合以往理論研究[1-7],根據(jù)流域水庫生產(chǎn)運(yùn)行實(shí)際,以調(diào)度期內(nèi)發(fā)電量最大為目標(biāo),在滿足約束條件下,采用離散動態(tài)規(guī)劃法,開發(fā)徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng),以期計(jì)算來水過程,實(shí)時得出水庫最優(yōu)運(yùn)行過程,為五凌電力集中控制運(yùn)行的各徑流式水電廠提供量化的指導(dǎo)依據(jù)。五凌電力所屬徑流式水電廠設(shè)計(jì)運(yùn)行參數(shù)見表1。

2模型構(gòu)建

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

以日為調(diào)度周期,小時(或30min、15min)為計(jì)算時段,在已知水電站水庫入流過程的條件下,求水電站出力過程和相應(yīng)的水庫蓄泄?fàn)顟B(tài)變化過程,使調(diào)度周期內(nèi)的發(fā)電量最大(現(xiàn)已取消豐枯、峰谷電價)。其目標(biāo)函數(shù)為

(1)

表1各徑流式水電廠設(shè)計(jì)運(yùn)行參數(shù)

水電廠正常蓄水位/m死水位/m總庫容/億m3調(diào)節(jié)庫容/億m3裝機(jī)容量/萬kW裝機(jī)臺數(shù)/臺單機(jī)容量/萬kW設(shè)計(jì)水頭/m掛治3223200.58590.0706153520.70凌津?yàn)?149.16.340.4627938.50東坪96.5920.1950.0927.241.86.80株溪口87.586.50.50.0527.441.856.10馬跡塘55.753.71.030.10594.231.46.55近尾洲6665.14.60.2436.31832.1066.80

式中 ,E為水電站調(diào)度期內(nèi)的總發(fā)電量,kW·h；Ni為i時段電站的出力,kW；Δt為時段長度,s；T為調(diào)度期的總時段數(shù)；qi為第i時段的發(fā)電流量,m3/s；ρi為第i時段的平均耗水率,m3/kW·h。其中,ρi=f(Hi),Hi=Z上i-Z下i,Z下i= f(qi),Z上i=f(Vi),Vi=Vi-1+(Qi-qi)Δt。式中 ,Hi為第i時段的水頭,m；ρi=f(Hi)為由水頭查耗水率；Z上i為第i時段庫水位,m；Vi為第i時段庫容,m3；Z上i=f(Vi)為由庫容查庫水位；Z下i為第i時段尾水位,m；Z下i=f(qi)為由出庫流量查尾水位；Qi為第i時段入庫流量,m3/s。

2.2約束條件

(1)水量平衡約束：Vi=Vi-1+(Qi-qi)Δt

(2)

(2)機(jī)組出力約束：Nmin≤Ni≤Nmax

(3)

(3)發(fā)電流量約束：qmin≤qi≤qmax

(4)

(4)水位約束：Zmin≤Z上i≤Zmax

(5)

式中,Nmin為機(jī)組限制最小出力(各單機(jī)最小限制出力之和),kW；Nmax為機(jī)組限制最大出力(各單機(jī)最大限制出力之和),kW；qmin為最小發(fā)電流量,m3/s；qmax為最大發(fā)電流量,一般為各單機(jī)最大過流量之和,m3/s；Zmin為限制最低庫水位,一般取死水位,m；Zmax為限制最高庫水位,一般取正常蓄水位,m。

3計(jì)算方法

3.1目標(biāo)函數(shù)的求解方法

目標(biāo)函數(shù)的優(yōu)化求解采用離散動態(tài)規(guī)劃法[7- 8]逐步逼近尋優(yōu)計(jì)算,其計(jì)算分為兩個過程：

(1)初始調(diào)度線選取。首先初選步長ΔZ1,通過動態(tài)規(guī)劃法求得初始優(yōu)化調(diào)度線{Z0(0),Z1(0),…,ZT(0)}。

(2)優(yōu)化調(diào)度計(jì)算。在初始優(yōu)化調(diào)度線上各取若干個增量ΔZ2,形成一個庫水位廊道,除最初、最末時段外,其余各時段形成3條調(diào)度線。用動態(tài)規(guī)劃法求得該廊道內(nèi)的優(yōu)化調(diào)度線{Z0(1),Z1(1),…,ZT(1)}。如果該優(yōu)化調(diào)度線各點(diǎn)均落在廊道內(nèi)的非邊界點(diǎn)上,或者雖在邊界點(diǎn)上但卻是該時段末(初)所允許的最大(小)水位,則該調(diào)度線即為所求最優(yōu)解。否則,在{Z0(1),Z1(1),…,ZT(1)}上下各取若干個增量ΔZ2,形成新的廊道,再用動態(tài)規(guī)劃法在該廊道內(nèi)尋優(yōu),直至優(yōu)化調(diào)度線全部為非廊道邊界點(diǎn)和最大(小)庫水位點(diǎn)所組成為止[8-9]。

3.2單庫動態(tài)規(guī)劃計(jì)算過程

(6)

利用式(6)進(jìn)行水庫調(diào)度需對狀態(tài)變量進(jìn)行離散化。對調(diào)度期內(nèi),假設(shè)庫水位在正常蓄水位和死水位之間連續(xù)變化,則可將水位劃分為m-1個網(wǎng)格,共m個點(diǎn),步長為ΔZ,如圖1所示。

圖1　動態(tài)規(guī)劃網(wǎng)格

運(yùn)用動態(tài)規(guī)劃法計(jì)算：已知日調(diào)節(jié)水庫的死水位、正常蓄水位、水位-庫容關(guān)系、出庫流量-尾水位關(guān)系、初始水位Z0,以及預(yù)報(bào)的各時段入庫流量Qi,其求解計(jì)算流程如圖2所示。

初選步長ΔZ1=0.1m,依據(jù)上述動態(tài)規(guī)劃法計(jì)算流程求得初始優(yōu)化調(diào)度線(水位精度為0.1m)；而后選取若干個增量ΔZ2=0.01m,再用上述動態(tài)規(guī)劃法在廊道內(nèi)選優(yōu)。由于ΔZ2=0.01m,優(yōu)化計(jì)算的各時段庫水位離理論最優(yōu)過程不超過0.01m,而目前五凌電力水情測報(bào)的水位精度為0.01m,故求得的優(yōu)化過程為目前測報(bào)精度下的最優(yōu)過程。

采用VB6.0編程工具,開發(fā)實(shí)時優(yōu)化調(diào)度系統(tǒng)。依托五凌電力辦公網(wǎng)絡(luò)環(huán)境,讀取五凌電力梯級水庫調(diào)度自動化系統(tǒng)oracle數(shù)據(jù)庫中日96點(diǎn)發(fā)電計(jì)劃、實(shí)時庫水位、尾水位等相關(guān)數(shù)據(jù),存入本地?cái)?shù)據(jù)庫,從而進(jìn)行實(shí)時優(yōu)化計(jì)算。

4參數(shù)及原始數(shù)據(jù)選取

系統(tǒng)優(yōu)化計(jì)算所需參數(shù)主要包括來流量、限制出力及水位、水頭-耗水率關(guān)系、水頭-單機(jī)限制出力關(guān)系、水位-庫容關(guān)系、出庫流量-尾水位關(guān)系等。其中,水頭-耗水率、出庫流量-尾水位采用最新的水庫實(shí)際運(yùn)行數(shù)據(jù)進(jìn)行擬合,圖3、4分別為凌津?yàn)┧姀S水頭-耗水率和出庫流量-尾水位的擬合關(guān)系；機(jī)組檢修、限制出力、限制水位等約束均根據(jù)運(yùn)行中實(shí)際要求進(jìn)行修改；期初庫水位、尾水位從五凌電力梯級水庫調(diào)度自動化系統(tǒng)oracle數(shù)據(jù)庫中實(shí)時讀取,期末庫水位推薦為能蓄至的最高水位(上限為正常蓄水位)。

圖3　凌津?yàn)┧姀S水頭-耗水率關(guān)系擬合

圖4　凌津?yàn)┧姀S出庫流量-尾水位關(guān)系擬合

對于模型優(yōu)化計(jì)算,需準(zhǔn)確可靠的入庫流量過程。根據(jù)電網(wǎng)下達(dá)的上游水電廠日96點(diǎn)計(jì)劃(從五凌電力梯級水庫調(diào)度自動化系統(tǒng)oracle數(shù)據(jù)庫中實(shí)時自動獲取),采用上游電廠“水頭-耗水率曲線”計(jì)算出發(fā)電流量過程,若上游電廠有棄水,讀取棄水流量,再加上區(qū)間流量,即得到24h小時的實(shí)時入庫流量過程[10]。

5實(shí)例計(jì)算及分析

5.1系統(tǒng)評價分析

選用6座徑流式水電廠2012年～2014年實(shí)際水情資料,各水電廠按預(yù)期控制水位運(yùn)行,可算得各電廠預(yù)期增加效益(見表2),與實(shí)際運(yùn)行相比,平均每年可增加發(fā)電量4 500萬kW·h。

表2徑流式水電廠預(yù)期年增加效益分析

電站實(shí)際年均發(fā)電量/億kW·h優(yōu)化年均發(fā)電量/億kW·h預(yù)期增加電量/億kW·h預(yù)期效益增加率/%掛治3.6413.7650.1243.4凌津?yàn)?1.77212.0140.2422.1東坪2.6722.6950.0230.9株溪口2.8192.8540.0351.3馬跡塘1.9751.9960.0211.1近尾洲3.3083.3150.0070.2

為進(jìn)一步對比說明,選取2015年9月28日凌津?yàn)┧姀S實(shí)際運(yùn)行情況與系統(tǒng)優(yōu)化計(jì)算進(jìn)行比較(期末水位相同)：系統(tǒng)優(yōu)化電量363.02萬kW·h,實(shí)際發(fā)電量358.18萬kW·h,增發(fā)電量4.84萬kW·h,增發(fā)1.35%；優(yōu)化調(diào)度系統(tǒng)計(jì)算的日均入庫流量為1 784m3/s,而實(shí)際測報(bào)計(jì)算的日均入庫流量為1 774m3/s,僅相差0.56%(準(zhǔn)確率達(dá)99.44%)。圖5為系統(tǒng)優(yōu)化計(jì)算結(jié)果及過程與實(shí)際運(yùn)行對比。

圖5　2015年9月28日凌津?yàn)┧姀S系統(tǒng)優(yōu)化結(jié)果與實(shí)際運(yùn)行過程比較

由圖5可知,在來流量增加前的入庫流量較小階段(3∶00～5∶00),應(yīng)適當(dāng)維持一定出力,一方面可在水頭較高情況下適當(dāng)發(fā)電；另一方面,可以適當(dāng)“騰庫”,在流量增加階段出力不必增加過快,避免明顯的抬高尾水位而降低發(fā)電水頭。在來流量持續(xù)增加階段(6∶00～10∶00),應(yīng)根據(jù)入庫流量情況適當(dāng)加大出力,由于之前的“騰庫”,所以該階段不必過大、過快加大出力,以便維持較高運(yùn)行水頭；否則,該階段出力快速、大幅增加,抬高了尾水位,從而明顯降低了運(yùn)行水頭,使得該部分較大出力時段的耗水率偏大,整體發(fā)電效益達(dá)不到最優(yōu)。在來流量較大且處于穩(wěn)定階段(18∶00～20∶00),盡量高水位運(yùn)行,合理安排出力維持出入庫平衡。在來流量減小階段(20∶00～24∶00),應(yīng)根據(jù)入庫情況適當(dāng)逐步降低出力以維持高水位。

所以,水庫運(yùn)行不能只關(guān)注高水位運(yùn)行,應(yīng)根據(jù)來水情況優(yōu)化出力安排,降低整體耗水率,提高發(fā)電效益。通過該系統(tǒng),可快速計(jì)算得出各種入庫流量變化情況下整體發(fā)電效益最大的各時段出力安排與運(yùn)行庫水位。

5.2實(shí)時優(yōu)化調(diào)度及分析

徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng)以目前時刻為起點(diǎn),進(jìn)行未來24h的實(shí)時優(yōu)化計(jì)算,得出可視化的水庫最優(yōu)運(yùn)行過程,指導(dǎo)水庫實(shí)時運(yùn)行。以凌津?yàn)┧姀S為例進(jìn)行計(jì)算說明。

現(xiàn)以當(dāng)前時刻(11月17日08∶15)凌津?yàn)┧姀S進(jìn)行實(shí)時優(yōu)化調(diào)度計(jì)算：點(diǎn)擊“凌津?yàn)┧姀S實(shí)時優(yōu)化調(diào)度系統(tǒng)”,顯示各計(jì)算參數(shù),如五強(qiáng)溪日平均入庫流量、五凌區(qū)間流量、期末推薦庫水位、限制最高、最低庫水位等,各參數(shù)均可根據(jù)實(shí)際發(fā)生的情況進(jìn)行修改；然后點(diǎn)擊“讀取參數(shù)”,各時段入庫流量、檢修臺數(shù)等亦可進(jìn)行修改；最后點(diǎn)擊“開始”按鈕進(jìn)行優(yōu)化計(jì)算,得到優(yōu)化計(jì)算結(jié)果界面如圖6所示。由于設(shè)置了“時間間隔”,故系統(tǒng)可在設(shè)定的時間間隔,自動重新以當(dāng)前時刻的參數(shù)、條件進(jìn)行未來24h的優(yōu)化計(jì)算,從而實(shí)現(xiàn)實(shí)時滾動的優(yōu)化計(jì)算。

圖6　2015年11月17日08∶15凌津?yàn)┧姀S實(shí)時優(yōu)化調(diào)度計(jì)算結(jié)果

圖6中計(jì)算結(jié)果顯示了庫水位、出力、入庫流量、出庫流量、水頭等水庫運(yùn)行關(guān)鍵過程,以及可發(fā)電量、實(shí)發(fā)電量、蓄能電量、耗水率等優(yōu)化計(jì)算值,其中：實(shí)發(fā)電量為由當(dāng)前時刻到期末推薦水位的調(diào)度期內(nèi)優(yōu)化計(jì)算發(fā)電量,蓄能電量為期末、期初庫容之差的水量按當(dāng)日平均發(fā)電耗水率計(jì)算的發(fā)電量,可發(fā)電量為實(shí)發(fā)電量與蓄能電量之和。

同時,水庫在運(yùn)行過程中,若需指定時段的庫水位,則可在“讀取參數(shù)”步驟中,選擇所需的計(jì)算時段和設(shè)置其控制水位,然后開始計(jì)算便可得到該時間內(nèi)的最優(yōu)過程。即該系統(tǒng)可根據(jù)運(yùn)行需要,選擇24h內(nèi)任意結(jié)束的計(jì)算時段、滿足庫水位、機(jī)組運(yùn)行臺數(shù)等需要的最優(yōu)運(yùn)行過程。

因?yàn)楸鞠到y(tǒng)采用離散動態(tài)規(guī)劃法逐步逼近尋優(yōu)求解步長為0.01m,優(yōu)化計(jì)算結(jié)果能達(dá)到目前水情測報(bào)精度(0.01m)下的最優(yōu)過程,且經(jīng)計(jì)算比較,與使用常規(guī)動態(tài)規(guī)劃法直接求解的結(jié)果相同；同時,當(dāng)計(jì)算時段長選為1h時,本系統(tǒng)計(jì)算時間只需0.250s(圖6),而使用動態(tài)規(guī)劃直接計(jì)算則需2.520s,若計(jì)算時段長選為30min或15min時,使用動態(tài)規(guī)劃法直接計(jì)算耗時會更長。所以本系統(tǒng)節(jié)省計(jì)算時間顯著且能達(dá)到最優(yōu)值。

6結(jié)語

開發(fā)的徑流式水電廠實(shí)時優(yōu)化調(diào)度系統(tǒng),可根據(jù)上游電廠發(fā)電計(jì)劃和水庫運(yùn)行情況以及區(qū)間來水情況,較為準(zhǔn)確計(jì)算入庫流量過程,結(jié)合實(shí)際運(yùn)行的條件和要求設(shè)定各參數(shù)值,能快速計(jì)算出最優(yōu)的水庫運(yùn)行過程,并能按設(shè)定的時間間隔實(shí)時滾動向前更新計(jì)算,能很好地實(shí)時量化指導(dǎo)水電廠的負(fù)荷安排及申請,從而獲得較好的經(jīng)濟(jì)效益,使水庫優(yōu)化調(diào)度理論較好地應(yīng)用于實(shí)際生產(chǎn)運(yùn)行。

參考文獻(xiàn)：

[1]鄧?yán)? 張璇, 楊永生, 等. 流域水資源調(diào)度研究綜述[J]. 水利經(jīng)濟(jì), 2011, 29(6)： 23- 27.

[2]唐明, 馬光文, 陶春華, 等. 水電站短期優(yōu)化調(diào)度模型的探討[J]. 水力發(fā)電, 2007, 35(5)： 88- 90.

[3]王興菊, 趙然杭. 水庫多目標(biāo)優(yōu)化調(diào)度理論及其應(yīng)用研究[J]. 水利學(xué)報(bào), 2003(3)： 104- 109.

[4]費(fèi)如君, 盧毓偉, 樊紹華, 等. 梯級水電站短期優(yōu)化調(diào)度算法研究及軟件開發(fā)[J]. 水力發(fā)電, 2013, 39(7)： 86- 88．

[5]黃益芬. 水電站水庫優(yōu)化調(diào)度研究[J]. 水力發(fā)電, 2002(5)： 64- 66．

[6]倪建軍, 徐立中, 李臣明, 等. 水庫調(diào)度決策研究綜述[J]. 水利水電科技進(jìn)展, 2004, 24(6)： 63- 66．

[7]暢建霞, 黃強(qiáng), 王義民. 水電站水庫優(yōu)化調(diào)度幾種方法的探討[J]. 水電能源科學(xué), 2000, 18(3)： 19- 22.

[8]鄒進(jìn), 劉可真. 水資源系統(tǒng)運(yùn)行與優(yōu)化調(diào)度[M]. 北京： 冶金工業(yè)出版社, 2006．

[9]楊侃. 水庫優(yōu)化調(diào)度中增量動態(tài)規(guī)劃收斂性研究[J]. 水科學(xué)進(jìn)展, 1995, 6(S1)： 23- 29.

[10]李益民, 段佳美. 水庫調(diào)度[M]. 2014年版. 北京： 中國電力出版社, 2014．

(責(zé)任編輯高瑜)

收稿日期：2015- 11- 19

作者簡介：肖楊(1986—),男,河南固始人,助理工程師,碩士,主要從事流域水資源優(yōu)化調(diào)度工作．

中圖分類號：TV697.1

文獻(xiàn)標(biāo)識碼：A

文章編號：0559- 9342(2016)04- 0097- 04

DevelopmentofReal-timeOptimizationDispatchingSystemforRun-offHydropowerStation

XIAOYang,KUANGLuzhang

(WulingPowerCorporation,Changsha410004,Hunan,China)

Abstract:In order to meet the requirements of hydropower station real-time scheduling, the real-time optimal scheduling model of run-off hydropower station is built with the goal of maximizing power generation. The discrete dynamic programming method is adopted in model calculation. The real-time optimal scheduling system can forecast the flow process of next 24 hours, and calculate the power generation process of next 24 hours which satisfies the requirements of constraint conditions in the ways of real-time and rolling. The system has been used in six run-off hydropower stations, and the results show that the accuracy rate of inflow calculation is more than 98% and the average annual power generation increases 45 million kW·h．

Key Words:run-off hydropower station; real-time optimization dispatching; discrete dynamic programming