徐 剛,夏 甜

(1.三峽大學(xué)水利與環(huán)境學(xué)院,湖北 宜昌 443002; 2.四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610065)

我國在長(zhǎng)江流域、黃河流域和各大中小型流域相繼形成了梯級(jí)開發(fā)的模式。流域梯級(jí)電站以及跨流域梯級(jí)電站集中向電網(wǎng)供電是未來流域梯級(jí)電站聯(lián)合調(diào)度的發(fā)展趨勢(shì)。國內(nèi)水庫調(diào)度圖常采用常規(guī)制作方法,即以典型年徑流作為入流,進(jìn)行等流量或等出力模擬調(diào)度[1],取歷年蓄水指示線包絡(luò)線作為調(diào)度線。王平[2]根據(jù)前人的研究經(jīng)驗(yàn),指出了水庫發(fā)電調(diào)度圖常規(guī)制作方法存在的問題,并提出了改進(jìn)方法,采用長(zhǎng)系列時(shí)歷法計(jì)算降低出力線,優(yōu)先滿足水電站降低出力的要求,然后從降低出力線至正常蓄水位之間的區(qū)域中劃分出保證出力區(qū)和加大出力區(qū),協(xié)調(diào)了水庫調(diào)度圖的各個(gè)調(diào)度區(qū)。Xie等[3]基于求解得出的大量確定性優(yōu)化調(diào)度結(jié)果分析得出發(fā)電調(diào)度規(guī)則,模擬調(diào)度結(jié)果滿足經(jīng)濟(jì)效益高的要求。隨著動(dòng)態(tài)規(guī)劃法、非線性規(guī)劃法以及智能算法在水庫優(yōu)化調(diào)度中的應(yīng)用,采用等出力或等流量徑流調(diào)節(jié)計(jì)算方法進(jìn)行水庫徑流調(diào)節(jié)顯得過于保守,浪費(fèi)了一些本可以利用的水資源。另外由于代表年的代表性不足,采用常規(guī)制作方法編制的水庫調(diào)度圖,不能保證準(zhǔn)確滿足設(shè)計(jì)保證率以及枯水期降低出力要求,需要采用長(zhǎng)系列資料進(jìn)行反復(fù)檢驗(yàn)和修正,而這種修正也是經(jīng)驗(yàn)性的,并沒有準(zhǔn)確、快捷的求解方法。

為了改善水庫調(diào)度圖常規(guī)制作方法,更好地運(yùn)用水庫調(diào)度圖指導(dǎo)梯級(jí)電站水庫聯(lián)合調(diào)度,本文基于確定性優(yōu)化調(diào)度模型的水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法制作水庫調(diào)度圖,在此基礎(chǔ)上綜合運(yùn)用水庫常規(guī)調(diào)度圖、水庫調(diào)度函數(shù)等水庫調(diào)度規(guī)則進(jìn)行梯級(jí)電站水庫聯(lián)合調(diào)度,并以岷江雜谷腦流域梯級(jí)電站聯(lián)合調(diào)度為實(shí)例對(duì)本文所提出的方法進(jìn)行驗(yàn)證。

1 梯級(jí)電站水庫聯(lián)合調(diào)度

1.1 確定性優(yōu)化調(diào)度模型

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

假定水庫以發(fā)電功能為主,此時(shí)優(yōu)化目標(biāo)為最大化發(fā)電量,即

(1)

式中:E為發(fā)電量;T為計(jì)算期的總時(shí)數(shù);A為出力系數(shù);Qt為第t時(shí)段電站引用流量;Ht為第t時(shí)段平均水頭;Δt為計(jì)算時(shí)段長(zhǎng)。

1.1.2 約束條件

水庫的水量平衡約束、水庫的蓄水量約束、電站引用流量約束、電站出力約束分別為

Vt+1=Vt+(It-Qt-St)Δt

(2)

Vt,min≤Vt≤Vt,max

(3)

Qt,min≤Qt≤Qt,max

(4)

Nmin≤AQtHt≤Nmax

(5)

式中:Vt、Vt+1分別為水庫第t時(shí)段初、末的蓄水量;It為第t時(shí)段水庫入庫流量;St為第t時(shí)段棄水流量;Vt,min為水庫發(fā)電允許的最少蓄水量,一般對(duì)應(yīng)死庫容;Vt,max為水庫發(fā)電允許的最大蓄水量,一般為汛限水位或正常蓄水位對(duì)應(yīng)的蓄水量;Qt,min為發(fā)電機(jī)組過機(jī)流量下限值;Qt,max為發(fā)電機(jī)組過機(jī)流量上限值;Nmin為水電站出力下限值;Nmax為水電站出力上限值,可綜合考慮機(jī)組額定出力及調(diào)峰要求等確定。另外,以上所有變量都必須滿足非負(fù)約束。

1.1.3 優(yōu)化算法

水庫優(yōu)化調(diào)度中的優(yōu)化算法有很多,主要以動(dòng)態(tài)規(guī)劃法[4]、POA逐次逼近算法[5-6]、GA遺傳算法[7-8]和模擬群集智能搜索算法[9-10]為代表。下面以模糊群集智能搜索算法中的蟻群算法[11]為例,介紹模擬群集智能搜索算法在水庫優(yōu)化調(diào)度中的應(yīng)用。

a. 設(shè)置蟻群算法參數(shù):種群大小M、最大進(jìn)化代數(shù)T、信息素?fù)]發(fā)參數(shù)ρ(0<ρ<1)、信息啟發(fā)式因子α、期望值啟發(fā)式因子、信息素痕跡初值τ0、閾值參數(shù)q0(0≤q0≤1)。

b. 初始蟻群:隨機(jī)生成m只人工螞蟻及路徑,并將路徑添加到路徑集合S中,初始化τi=τ0(i=1,2,…,m),τi代表路徑上的累積信息素。

c. 構(gòu)建路徑:先對(duì)每只螞蟻根據(jù)狀態(tài)轉(zhuǎn)移規(guī)則選擇新路徑,然后在[0,1]區(qū)間中產(chǎn)生隨機(jī)數(shù)q,再將q與q0進(jìn)行比較,當(dāng)q≤q0時(shí),對(duì)人工螞蟻所在的路徑Sold隨機(jī)地進(jìn)行變異,即對(duì)水庫水位離散序列Z={z1,z2,…,zi,…,zT}中的zi在搜索空間限定范圍內(nèi)隨機(jī)地進(jìn)行變異操作;當(dāng)q>q0時(shí),依據(jù)高概率原則在路徑集合S中選擇路徑s,其中概率P按式(6)計(jì)算:

(6)

其中ηis=Ei-Es(i,s∈S)

式中:τi、τs分別為路徑i、s相關(guān)的信息素痕跡;Ei、Es分別為路徑i、s的目標(biāo)函數(shù)值(發(fā)電量);ηis為用于評(píng)價(jià)螞蟻從路徑i向路徑s轉(zhuǎn)移的局部啟發(fā)函數(shù)。

d. 更新信息素:計(jì)算各螞蟻當(dāng)前路徑i的目標(biāo)函數(shù)值Ei(i=1,2,…,m)。對(duì)當(dāng)前路徑集合S按式(7)更新信息素痕跡:

τi=(1-ρ)τi+ρΔτi(i∈S)

(7)

e. 判斷是否滿足迭代終止條件和達(dá)到最大迭代次數(shù)。若是,則結(jié)束迭代,輸出結(jié)果;否則返回構(gòu)建路徑步驟,重新進(jìn)行迭代計(jì)算。

1.2 水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法

1.2.1 水庫調(diào)度圖改進(jìn)與優(yōu)化原則和思路

針對(duì)水庫調(diào)度圖常規(guī)制作方法的不足,對(duì)常規(guī)制作方法進(jìn)行的改進(jìn)與優(yōu)化遵循以下原則和思路:

a. 采用長(zhǎng)系列時(shí)歷法計(jì)算。由于常規(guī)方法計(jì)算的調(diào)度圖不能保證完全滿足設(shè)計(jì)保證率的要求,現(xiàn)有規(guī)范均要求再用長(zhǎng)系列資料進(jìn)行檢驗(yàn)和修正,所以改進(jìn)方法應(yīng)直接采用長(zhǎng)系列時(shí)歷資料計(jì)算調(diào)度圖,以協(xié)調(diào)徑流調(diào)節(jié)和水庫調(diào)度兩階段計(jì)算成果。

b. 采用上文描述的確定性優(yōu)化調(diào)度模型進(jìn)行徑流調(diào)節(jié)計(jì)算。隨著優(yōu)化算法的大規(guī)模運(yùn)用,如仍然采用常規(guī)方法求解歷年的水庫水位蓄水指示線,不符合水電站經(jīng)濟(jì)運(yùn)行要求。

c. 協(xié)調(diào)水庫不同調(diào)度區(qū)的關(guān)系。首先確定調(diào)度圖降低出力和降低出力限制線,優(yōu)先滿足降低出力時(shí)不會(huì)發(fā)生破壞,在降低出力線以上確定保證出力區(qū)和出力防破壞線,最后在出力防破壞線以上確定不同的加大出力線和加大出力區(qū)。

d. 不同調(diào)節(jié)性能水庫采用統(tǒng)一的計(jì)算方法。針對(duì)具有不完全年調(diào)節(jié)能力以上調(diào)節(jié)能力的水庫,采用統(tǒng)一的計(jì)算方法。

e. 水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法以計(jì)算機(jī)為工具。改進(jìn)計(jì)算方法將采用計(jì)算機(jī)實(shí)現(xiàn)自動(dòng)計(jì)算,方便地試算出合適的保證出力值、出力保證率、降低出力值和加大出力值。

1.2.2 水庫調(diào)度圖改進(jìn)與優(yōu)化制作步驟

步驟1 求解降低出力和降低出力限制線。選取長(zhǎng)系列中的枯水年,通過試算確定枯水期的降低出力,當(dāng)降低出力小到一定程度時(shí),枯水期各時(shí)段出力均不被破壞,取此時(shí)水庫蓄水指示線上包絡(luò)線作為枯水期降低出力限制線。

步驟2 求解保證出力和保證出力區(qū)。選取長(zhǎng)系列中的平水年,通過試算求解設(shè)定的保證率對(duì)應(yīng)的出力作為保證出力,取歷年正常出力時(shí)段的水庫蓄水指示線下包絡(luò)線作為正常出力時(shí)段的降低出力限制線,正常出力時(shí)段的降低出力限制線至水庫正常蓄水位之間的區(qū)域作為可以發(fā)出不小于保證出力的區(qū)域。以正常出力時(shí)段降低出力限制線至第1根加大出力線之間的區(qū)域作為保證出力區(qū),第1根加大出力線也是出力防破壞線。

步驟3 求解加大出力和加大出力區(qū)。當(dāng)水庫有加大出力需求時(shí),需從可以發(fā)出不小于保證出力的區(qū)域中分解出不同的加大出力區(qū)。選取長(zhǎng)系列中的豐水年作為入流資料,設(shè)置加大出力值,取歷年同期水庫蓄水指示線下包絡(luò)線作為加大出力線,加大出力線以上區(qū)域?yàn)榧哟蟪隽^(qū)。預(yù)想出力線為預(yù)想出力情況下歷年水庫蓄水指示線下包絡(luò)線,預(yù)想出力線以上、正常蓄水位以下的區(qū)域?yàn)轭A(yù)想出力區(qū)。

步驟4 調(diào)度線修正。為避免調(diào)度線交叉,應(yīng)確保加大出力線在降低出力線之上。

1.3 梯級(jí)電站水庫調(diào)度規(guī)則的綜合運(yùn)用

通常在流域梯級(jí)電站中,梯級(jí)電站龍頭水庫的調(diào)節(jié)能力較強(qiáng),整個(gè)梯級(jí)電站的總出力與龍頭水庫的調(diào)蓄狀態(tài)相關(guān)性較強(qiáng),而下游電站調(diào)節(jié)能力一般較弱,有些甚至只具有日調(diào)節(jié)能力,其水位在長(zhǎng)時(shí)間內(nèi)變化不大,繪制水位-出力調(diào)度圖有所不便,同時(shí)下游電站出力和上游梯級(jí)電站出力有關(guān), 可以采用調(diào)度函數(shù)[12-13]作為下游電站決策方案。因此,將流域中龍頭電站的出力決策作為流域梯級(jí)電站聯(lián)合調(diào)度的核心決策,下游電站出力決策作為核心決策的分解決策。流域梯級(jí)電站聯(lián)合調(diào)度決策方案包括龍頭電站調(diào)度圖和下游電站調(diào)度函數(shù)，不同決策方案的計(jì)算流程如圖1所示。

圖1 流域梯級(jí)電站聯(lián)合調(diào)度不同決策方案計(jì)算流程

2 雜谷腦流域梯級(jí)電站聯(lián)合調(diào)度實(shí)例

雜谷腦流域獅子坪、薛城和古城3個(gè)梯級(jí)電站采取聯(lián)合調(diào)度方式進(jìn)行集中控制調(diào)度。獅子坪電站為雜谷腦河干流的龍頭水庫電站,電站裝機(jī)容量195 MW,水庫正常蓄水位2 540.0 m,死水位2 460.0 m,調(diào)節(jié)庫容1.189億m3,具有年調(diào)節(jié)能力。薛城電站裝機(jī)容量138 MW,水庫正常蓄水位1 709.5 m,相應(yīng)庫容114.8萬m3,汛期排沙運(yùn)行水位1 704.0 m,水庫死水位1 704.0 m,調(diào)節(jié)庫容62.3萬m3,具有日調(diào)節(jié)能力。古城電站裝機(jī)容量168 MW,水庫正常蓄水位1 554.5 m,相應(yīng)庫容93.9萬m3,汛期排沙運(yùn)行水位1 550.0 m,水庫死水位1 550.0 m,調(diào)節(jié)庫容58萬m3,具有日調(diào)節(jié)能力。

2.1 獅子坪電站水庫常規(guī)調(diào)度圖的改進(jìn)與優(yōu)化

遵循水庫調(diào)度圖改進(jìn)與優(yōu)化原則和思路,利用水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法繪制獅子坪電站水庫調(diào)度圖。該方法采用C#計(jì)算機(jī)程序?qū)崿F(xiàn),基于Visual Studio 2005平臺(tái),采用語言編程環(huán)境開發(fā)出基于蟻群算法的水庫優(yōu)化調(diào)度模型作為程序徑流調(diào)節(jié)計(jì)算模塊,優(yōu)化調(diào)度模塊的輸入為長(zhǎng)系列徑流資料,輸出為水庫歷年蓄水過程線,再取相應(yīng)的蓄水指示線包絡(luò)線作為調(diào)度線。獅子坪電站水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法的具體步驟如下：

2.1.1 求解降低出力和降低出力限制線

選取長(zhǎng)系列中經(jīng)驗(yàn)頻率大于70%的枯水年作為入流,初定1—4月最小出力下限值為60 MW,11—12月為80 MW,汛期不設(shè)最小值,通過不斷調(diào)整枯水期時(shí)段的最小出力下限值,試算出枯水期時(shí)段的降低出力,這一降低出力應(yīng)能使長(zhǎng)系列的 各枯水期時(shí)段都能滿足最小出力而不被破壞。當(dāng)1—4月降低出力為35 MW、11—12月降低出力為50 MW時(shí),長(zhǎng)系列資料中所有年份枯水期時(shí)段出力全部不被破壞,試算結(jié)束。取整個(gè)長(zhǎng)系列徑流調(diào)節(jié)過程的枯水期時(shí)段蓄水指示線的上包絡(luò)線作為對(duì)應(yīng)的降低出力限制線,在該指示線以上可以發(fā)大于降低出力的出力,指示線以下必須發(fā)小于降低出力的出力?？菟跁r(shí)段降低出力限制線和死水位之間的區(qū)域?yàn)榻档统隽^(qū)。

圖2 獅子坪電站水庫調(diào)度圖

2.1.2 求解保證出力和保證出力區(qū)

選取長(zhǎng)系列中經(jīng)驗(yàn)頻率在30%～70%之間的平水年作為入流,根據(jù)已經(jīng)確定的降低出力,把相應(yīng)破壞時(shí)段的最小出力下限值分別設(shè)定為各自的降低出力,其他正常出力時(shí)段的最小出力下限值設(shè)置為195 MW,也就是滿發(fā),利用上文介紹的優(yōu)化調(diào)度模型進(jìn)行長(zhǎng)系列徑流調(diào)節(jié)計(jì)算,通過不斷調(diào)整最小出力下限值,保證率也會(huì)相應(yīng)增大,當(dāng)出力下限值設(shè)置為120 MW時(shí),保證率可以達(dá)到82.8%。取試算的長(zhǎng)系列所有年份蓄水過程線對(duì)應(yīng)的水庫蓄水指示線下包絡(luò)線作為正常出力時(shí)期降低出力限制線。 以正常出力時(shí)期降低出力限制線至第1根加大出力線之間的區(qū)域作為保證出力區(qū)。

2.1.3 求解加大出力線和加大出力區(qū)

選取長(zhǎng)系列中經(jīng)驗(yàn)頻率小于30%的豐水年作為入流,獅子坪電站最有可能加大出力的時(shí)段為5—8月,利用優(yōu)化調(diào)度模型進(jìn)行徑流調(diào)節(jié)計(jì)算,設(shè)定豐水期時(shí)段的出力下限值為已經(jīng)確定的降低出力,正常時(shí)段5—8月出力下限值設(shè)置為195 MW,經(jīng)過試算正常工作時(shí)段電站出力為195 MW時(shí)不會(huì)造成破壞,取該徑流調(diào)節(jié)過程的所有蓄水指示線下包絡(luò)線作為加大出力線。分別設(shè)置加大出力值為195 MW、170 MW和145 MW,繪制出相應(yīng)的加大出力線。 不同的加大出力線之間的區(qū)域分別對(duì)應(yīng)不同的加大出力區(qū),其中 195 MW加大出力區(qū)為預(yù)想出力區(qū)。

2.1.4 調(diào)度線的修正

若調(diào)度線上個(gè)別點(diǎn)造成調(diào)度線交叉,利用計(jì)算機(jī)對(duì)該點(diǎn)進(jìn)行修正,確保調(diào)度線不交叉。利用水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法繪制的獅子坪電站水庫調(diào)度圖如圖2所示。

2.2 薛城和古城電站水庫運(yùn)行策略的擬定

根據(jù)水庫調(diào)度圖改進(jìn)與優(yōu)化制作方法,計(jì)算出獅子坪電站水庫調(diào)度圖,形成雜谷腦梯級(jí)電站聯(lián)合調(diào)度決策方案中的核心決策。將歷年來水作為輸入,以獅子坪電站水庫調(diào)度圖作為調(diào)度規(guī)則進(jìn)行模擬調(diào)度,求解獅子坪電站1977—2005年的調(diào)度過程。通過分析獅子坪電站29年的模擬調(diào)度數(shù)據(jù),采用多元回歸分析分別確定薛城電站和古城電站與獅子坪電站出力和天然來水之間逐月的相關(guān)關(guān)系,形成了薛城電站和古城電站的水庫調(diào)度函數(shù),即雜谷腦梯級(jí)電站中的分解決策一和分解決策二。薛城電站和古城電站調(diào)度函數(shù)分別為

Nxc=aQt+bNszp+c

(8)

Ngc=aQt+bNszp+c

(9)

式中:a、b、c為參數(shù),其值見表1;Nxc為薛城電站出力;Ngc為古城電站出力;Nszp為獅子坪電站出力;Qt為天然來水流量。

表1 薛城和古城電站調(diào)度函數(shù)參數(shù)的取值及復(fù)相關(guān)系數(shù)

2.3 雜谷腦流域梯級(jí)電站聯(lián)合調(diào)度方案的檢驗(yàn)

2.3.1 獅子坪電站調(diào)度圖檢驗(yàn)

以全部歷史來水系列作為輸入,以調(diào)度圖作為調(diào)度規(guī)則,逐年計(jì)算出獅子坪電站水庫的蓄泄過程、出力過程等,計(jì)算出水庫按照調(diào)度圖運(yùn)行所具有的保證率。獅子坪電站除了最大限度地滿足較枯水期破壞時(shí)段的出力外,正常工作時(shí)段平均出力均大于120 MW。因此獅子坪電站按照改進(jìn)與優(yōu)化制作方法繪制水庫調(diào)度圖是合理可行的。

2.3.2 薛城和古城電站調(diào)度函數(shù)的檢驗(yàn)

在對(duì)變量進(jìn)行線性回歸分析時(shí),采用最小二乘法進(jìn)行參數(shù)估計(jì),復(fù)相關(guān)系數(shù)R2為回歸平方和與總離差平方和的比值,R2越大,模型越精確,回歸效果越顯著。擬合了薛城和古城電站水庫調(diào)度函數(shù),通過計(jì)算各函數(shù)的R2來檢驗(yàn)函數(shù)的擬合效率,R2計(jì)算結(jié)果如表1所示,可見大部分月份的R2值都在0.9以上,由于12月處于豐枯銜接時(shí)段,來水偶然性較大,導(dǎo)致線性相關(guān)性不如其他月份顯著。對(duì)R2進(jìn)行F檢驗(yàn),計(jì)算公式如下:

(10)

式中:n為樣本數(shù)目;k為自變量維度。

計(jì)算得出F=6.45,通過查詢F分布表[14]得知F0.01=5.53,F>F0.01表明復(fù)相關(guān)關(guān)系在置信度水平α=0.01上顯著,因而本例中調(diào)度函數(shù)的擬和效果較好。

2.3.3 雜谷腦梯級(jí)電站水庫聯(lián)合調(diào)度方案及其檢驗(yàn)

通過將雜谷腦流域梯級(jí)電站聯(lián)合調(diào)度決策方案分解為獅子坪電站水庫調(diào)度圖、薛城和古城電站調(diào)度函數(shù),可以實(shí)現(xiàn)梯級(jí)電站聯(lián)合調(diào)度。根據(jù)獅子坪電站水庫調(diào)度圖當(dāng)前時(shí)段初的出力,求解獅子坪下泄流量,利用薛城、古城電站調(diào)度函數(shù)計(jì)算薛城電站和古城電站的發(fā)電出力,計(jì)算獅子坪、薛城和古城3個(gè)電站出力總和。

分別計(jì)算了優(yōu)化調(diào)度、改進(jìn)與優(yōu)化調(diào)度圖(采用改進(jìn)與優(yōu)化制作方法繪制的水庫調(diào)度圖)結(jié)合調(diào)度函數(shù)以及常規(guī)調(diào)度圖結(jié)合調(diào)度函數(shù)3種梯級(jí)電站聯(lián)合調(diào)度方式的梯級(jí)電站多年平均發(fā)電量,計(jì)算結(jié)果分別為21.124億kW·h、20.451億kW·h和17.632億kW·h。結(jié)果表明: 改進(jìn)與優(yōu)化調(diào)度圖結(jié)合調(diào)度函數(shù)方式的梯級(jí)電站多年平均發(fā)電量只比優(yōu)化調(diào)度方式的發(fā)電量降低了3.3%,比常規(guī)調(diào)度圖結(jié)合調(diào)度函數(shù)方式提高了16.0%; 改進(jìn)與優(yōu)化方法繪制的水庫調(diào)度圖結(jié)合調(diào)度函數(shù)方式在保證獅子坪電站水庫調(diào)度圖設(shè)計(jì)保證率較高的同時(shí),還顯著增加了梯級(jí)電站的多年平均發(fā)電量。

3 結(jié) 語

針對(duì)水庫調(diào)度圖常規(guī)制作方法的不足,用長(zhǎng)系列法替代典型年法, 按照先滿足破壞時(shí)段降低出力的要求繪制降低出力線,后繪制出力防破壞線和加大出力線的方式,協(xié)調(diào)了水庫調(diào)度圖的不同調(diào)度區(qū),滿足了降低出力的要求,提高了水庫調(diào)度圖的保證率和實(shí)用性;基于水庫確定性優(yōu)化調(diào)度模型的徑流調(diào)節(jié),優(yōu)化了水庫調(diào)度方式。采用改進(jìn)與優(yōu)化調(diào)度圖結(jié)合調(diào)度函數(shù)的梯級(jí)電站聯(lián)合調(diào)度方式,在提高調(diào)度圖保證率的同時(shí),也增加了梯級(jí)電站聯(lián)合調(diào)度的多年平均發(fā)電量。

參考文獻(xiàn)：

[1] 張秀菊.水電站保證出力的計(jì)算方法[J].水力發(fā)電,2006,32(10):21-23.(ZHANG Xiuju. Calculation method of firm output of hydropower station[J].Water Power,2006,32(10):21-23.(in Chinese)

[2] 王平.水庫發(fā)電調(diào)度圖常規(guī)計(jì)算方法的問題和改進(jìn)建議[J].水利水電工程設(shè)計(jì),2007,26(3):26-29.(WANG Ping. The problems and the improving suggestions for conventional calculation method of power generation dispatching drawings of reservoir[J].Design of Water Resources & Hydroelectric Engineering,2007,26(3):26-29.(in Chinese)

[3] XIE Wei, JI Changming, YANG Zijun, et al. Short-term power generation scheduling rules for cascade hydropower stations based on hybrid algorithm[J].Water Science and Engineering, 2012, 5(1):46-58.

[4] 路志宏,魏守平.變狀態(tài)空間動(dòng)態(tài)規(guī)劃法在水庫優(yōu)化調(diào)度中的應(yīng)用[J].水電能源科學(xué),2003,21(2):55-57.(LU Zhihong, WEI Shouping. Status variable space dynamic programming and its application for optimal operation of cascaded hydroelectric systems[J].Water Resources and Power,2003,21(2):55-57.(in Chinese)

[5] 宗航,李承軍,周建中,等.POA算法在梯級(jí)水電站短期優(yōu)化調(diào)度中的應(yīng)用[J].水電能源科學(xué),2003,21(1):46-48. (ZONG Hang, LI Chengjun, ZHOU Jiangzhong, et al. Research and application for short-term cascaded hydroelectric scheduling based on progressive optimality algorithm[J].Water Resources and Power,2003,21(1):46-48.(in Chinese)

[6] 周佳,馬光文,張志剛.基于改進(jìn)POA算法的雅礱江梯級(jí)水電站群中長(zhǎng)期優(yōu)化調(diào)度研究[J],水力發(fā)電學(xué)報(bào),2010,29(3):18-22.(ZHOU Jia, MA Guangwen, ZHANG Zhigang. Study on the mid-long term optimal dispatching of cascaded hydropower stations on Yalong River based on POA modified adaptive algorithm[J].Journal of Hydroelectric Engineering,2010,29(3):18-22.(in Chinese)

[7] 劉攀,郭生練,李瑋,等.遺傳算法在水庫調(diào)度中的應(yīng)用綜述[J].水利水電科技進(jìn)展,2006,29(4):78-83.(LIU Pan, GUO Shenglian, LI Wei, et al. A review of application of genetic algorithm to reservoir operation[J].Advances in Science and Technology of Water Resources,2006,29(4):78-83.(in Chinese)

[8] 何向陽,周建中,張勇傳,等.基于改進(jìn)NSGA-Ⅱ的梯級(jí)水電站多目標(biāo)發(fā)電優(yōu)化調(diào)度[J]. 武漢大學(xué)學(xué)報(bào):工學(xué)版,2011,44(6):715-719.(HE Xiangyang, ZHOU Jiangzhong, ZHANG Yongchuan, et al. Multiobjective optimal dispatching of cascade hydropower stations using improved NSGA-Ⅱ[J].Engineering Journal of Wuhan University,2011,44(6):715-719.(in Chinese)

[9] 朱記偉,張洪波,辛琛,等.水庫調(diào)度PSO優(yōu)化模型及求解方法[J].系統(tǒng)工程,2010,28(1):105-112.(ZHU Jiwei, ZHANG Hongbo, XIN Sheng, et al. PSO Model and calculation method to reservoir optimal operation[J].Systems Engineering,2010,28(1):105-112.(in Chinese)

[10] 向波,紀(jì)昌明,羅慶松.免疫粒子群算法及其在水庫優(yōu)化調(diào)度中的應(yīng)用[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2008,36(2):198-201.(XIANG Bo, JI Changming, LUO Qinsong. Immune particle swarm optimization algorithm and its application in reservoir operation optimization[J].Journal of Hohai University:Natural Sciences,2008,36(2):198-201.(in Chinese)

[11] 徐剛,馬光文,梁武湖,等.蟻群算法在水庫優(yōu)化調(diào)度中的應(yīng)用[J].水科學(xué)進(jìn)展,2005,16(3):397-400. (XU Gang, MA Guangwen, LIANG Wuhu, et al. Application of ant colony algorithm to reservoir optimal operation[J].Advances in Water Science,2005,16(3):397-400.(in Chinese)

[12] 權(quán)先璋,李承軍,張士軍,等.水電站優(yōu)化線性調(diào)度規(guī)則研究[J]. 華中理工大學(xué)學(xué)報(bào),1999,27(12):36-38.(QUAN Xianzhang, LI Chengjun, ZHANG Shijun, et al.A study on the optimal linear dispatching rules for hydropower station[J].Journal of Huazhong University of Science and Technology,1999,27(12):36-38.(in Chinese)

[13] 李承軍,陳畢勝,張高峰.水電站雙線性調(diào)度規(guī)則研究[J].水力發(fā)電學(xué)報(bào),2005,24(1):11-15.(LI Chengjun, CHEN Bisheng, ZHANG Gaofeng. Study on bilinear dispatching rule of hydropower station[J].Journal of Hydroelectric Engineering,2005,24(1):11-15.(in Chinese)

[14] 盛驟,謝式千,潘承毅.概率論與數(shù)理統(tǒng)計(jì)[M].3版.北京:高等教育出版社,2001.