徐 敏,周建中,歐陽(yáng)文宇,覃 暉
(1.中國(guó)電建集團(tuán)貴陽(yáng)勘測(cè)設(shè)計(jì)研究院有限公司,貴州 貴陽(yáng) 550081;2.華中科技大學(xué)水電與數(shù)字化工程學(xué)院,湖北 武漢 430074)
對(duì)于以發(fā)電為主要任務(wù)的水庫(kù),在沒(méi)有水文預(yù)報(bào)或者水文預(yù)報(bào)不可靠的情況下,結(jié)合歷史徑流序列等統(tǒng)計(jì)資料,分析和制作水庫(kù)發(fā)電調(diào)度圖,能用于制訂調(diào)度計(jì)劃,并指導(dǎo)調(diào)度決策[1]。且調(diào)度圖具備表達(dá)直觀、操作簡(jiǎn)便、物理意義清晰、應(yīng)用可靠等優(yōu)點(diǎn),所以大多數(shù)具有年或年以上調(diào)節(jié)性能的水電站水庫(kù)通常采用調(diào)度圖進(jìn)行調(diào)度。然而,采用調(diào)度圖指導(dǎo)水庫(kù)調(diào)度存在著正常出力區(qū)范圍太大,難以達(dá)到全局最優(yōu)和準(zhǔn)全局最優(yōu)等問(wèn)題[2]。因此,長(zhǎng)期以來(lái),水庫(kù)發(fā)電調(diào)度圖的繪制方法一直處于逐漸改進(jìn)的過(guò)程中,不少學(xué)者均嘗試改造傳統(tǒng)調(diào)度圖,或在調(diào)度圖繪制中引入優(yōu)化思想[3-13]。但由于大多方法改變了成圖方式或調(diào)度圖本身,都不同程度地削弱了調(diào)度圖表達(dá)直觀、操作簡(jiǎn)便、物理意義清楚、應(yīng)用可靠等優(yōu)點(diǎn)。其中,對(duì)單庫(kù)調(diào)度圖的改進(jìn)或優(yōu)化沒(méi)有充分考慮梯級(jí)上下游水庫(kù)之間的徑流聯(lián)系和梯級(jí)整體效益,而梯級(jí)調(diào)度圖優(yōu)化由于缺乏實(shí)踐基礎(chǔ)而導(dǎo)致適用性不強(qiáng),因此在實(shí)際應(yīng)用中均存在一定局限性。
為克服傳統(tǒng)調(diào)度圖在制作時(shí)依賴工作人員實(shí)際調(diào)度經(jīng)驗(yàn),不夠嚴(yán)謹(jǐn)[14],且繪制過(guò)程中輸入信息以偏枯年份徑流為主,往往造成模擬運(yùn)行中水庫(kù)水量利用率偏低的缺點(diǎn),并避免已有優(yōu)化調(diào)度圖繪制方法的缺陷;本文提出了一種在繪圖過(guò)程中考慮下游電站,同時(shí)引入優(yōu)化思想控制調(diào)度圖繪制參數(shù)的年調(diào)節(jié)水電站發(fā)電調(diào)度圖繪制方法。該方法結(jié)合并改進(jìn)現(xiàn)有調(diào)度圖繪制方式,很大程度地保留了傳統(tǒng)調(diào)度圖的優(yōu)點(diǎn),能得到具有優(yōu)化屬性、實(shí)用性強(qiáng)、以及物理意義和控制作用明確的發(fā)電優(yōu)化調(diào)度圖。
對(duì)于年調(diào)節(jié)水電站水庫(kù)而言,水庫(kù)豐、枯徑流特性及自身調(diào)節(jié)能力使得水庫(kù)蓄、供水期調(diào)度目標(biāo)差異較大。因此,水電站供蓄水期分別繪制調(diào)度圖更能適應(yīng)徑流變化特點(diǎn)和水庫(kù)發(fā)電調(diào)度需求[5]。為此,首先需進(jìn)行水庫(kù)供蓄水期分期劃分。多篇文獻(xiàn)均給出了分期方法[5,15],此處不再贅述。
在分期調(diào)度圖的繪制中,基本參數(shù)包括供、蓄水期保證出力,加大、降低出力,代表年徑流,典型年徑流等,其中保證出力需優(yōu)先計(jì)算,其余參數(shù)均以其為依據(jù)產(chǎn)生。與常規(guī)繪制方法[14]不同,本文采用分期繪制調(diào)度圖方式,因此需分別針對(duì)供、蓄水期進(jìn)行計(jì)算。供水期保證出力計(jì)算方法與傳統(tǒng)調(diào)度圖一致,具體繪制方法見(jiàn)文獻(xiàn)[16]。針對(duì)蓄水期,從蓄水期初死水位開(kāi)始進(jìn)行順時(shí)序等流量調(diào)節(jié)計(jì)算,求出各蓄水期平均出力,排序后取保證率對(duì)應(yīng)出力為蓄水期保證出力。但模擬運(yùn)行中,如果按照此出力進(jìn)行計(jì)算可能導(dǎo)致一些年份由于汛期出力過(guò)大使供水期初水位達(dá)不到正常高水位,造成供水期不能充分供水而成為破壞年份,使得設(shè)計(jì)保證率無(wú)法保障[3]。為此,本文把蓄水期保證出力作為蓄水期指示出力上限,給定初始值后再依據(jù)優(yōu)化模型對(duì)指示出力進(jìn)行數(shù)值尋優(yōu)。
供蓄水期采取分別繪制并改進(jìn)加大或降低出力線的繪制方式。對(duì)選取的每一個(gè)典型徑流過(guò)程,采取等出力的方式,供水期從供水期末相應(yīng)的指示水位(死水位),逆時(shí)序計(jì)算至供水期初,回到正常蓄水位或其他要求的限制水位;蓄水期從蓄水期末相應(yīng)的指示水位(正常蓄水位),逆時(shí)序計(jì)算至蓄水期初(死水位)。典型徑流選取修正及每一個(gè)時(shí)段的計(jì)算方法詳見(jiàn)文獻(xiàn)[16]。以保證出力、指示出力為等出力值進(jìn)行計(jì)算時(shí),將所有的典型徑流過(guò)程進(jìn)行計(jì)算后,即得到一組對(duì)應(yīng)的水庫(kù)水位過(guò)程線;然后取其上、下包線,即得到供、蓄水期上、下基本調(diào)度線。
考慮當(dāng)徑流過(guò)程相同時(shí),對(duì)于逆時(shí)序計(jì)算[17]而言,通常情況下,出力越大,水位越高;出力越小,水位越低(見(jiàn)圖1)。故,當(dāng)輸入出力值較大時(shí),運(yùn)用典型徑流過(guò)程計(jì)算后,對(duì)所得水庫(kù)水位過(guò)程線取其上包線即可得位于保證出力區(qū)之上的曲線,此線與其下臨近曲線構(gòu)成一加大出力區(qū),當(dāng)水位到達(dá)此區(qū)時(shí)會(huì)加大出力使水位消落,因此可以起到控制水位的作用;同理,對(duì)于輸入出力較小的情況,取下包線可形成降低出力控制區(qū)。

圖1 加大、降低出力逆推所得調(diào)度線位置示意
因此,對(duì)各加大出力值運(yùn)用典型徑流過(guò)程進(jìn)行逆時(shí)序計(jì)算,取所得水庫(kù)水位過(guò)程線上包絡(luò)線可得各加大出力線;對(duì)各降低出力值運(yùn)用典型徑流過(guò)程進(jìn)行逆時(shí)序計(jì)算,取所得水庫(kù)水位過(guò)程線下包絡(luò)線可得各降低出力線;將各線與保證出力線整合,消除重疊部分,可得水電站發(fā)電調(diào)度圖,具體流程見(jiàn)圖2。
梯級(jí)水電站采用考慮保證率的梯級(jí)發(fā)電量最大模型,兼顧蓄末未蓄滿懲罰項(xiàng)。該模型給定調(diào)度期內(nèi)入庫(kù)徑流過(guò)程和水庫(kù)始末水位,綜合考慮各種約束條件,確定梯級(jí)各水電站水庫(kù)的發(fā)電蓄放水過(guò)程,使模擬調(diào)度期內(nèi)發(fā)電量最大。
目標(biāo)函數(shù)
(1)

圖2 改進(jìn)調(diào)度圖計(jì)算流程
式中,E為梯級(jí)水電站多年總發(fā)電量,億kW·h;N(i,m,t)為第i電站m年t時(shí)段的出力,kW;Ng(i)為梯級(jí)第i電站保證出力,kW;Nsg(i)為梯級(jí)第i電站模擬運(yùn)行統(tǒng)計(jì)保證出力,即運(yùn)行結(jié)果中年保證率對(duì)應(yīng)出力,kW;n(i)為第i水庫(kù)蓄水期末未蓄滿次數(shù);α、β為懲罰系數(shù);I為梯級(jí)電站總數(shù);M為徑流序列長(zhǎng)度,年;T為年內(nèi)單位時(shí)段總數(shù);Δt為單位時(shí)段長(zhǎng),h。
相關(guān)約束:
(1)水量平衡方程
V(i,m,t+1)=V(i,m,t)+[I(i,m,t)-Q(i,m,t)]Δt
(2)
式中,V(i,m,t)為i水庫(kù)m年t時(shí)段初始庫(kù)容,億m3;I(i,m,t)、Q(i,m,t)分別為i水庫(kù)m年t時(shí)段的平均入庫(kù)、出庫(kù)流量,m3/s。
(2)庫(kù)水位約束
Zmin(i,m,t)≤Z(i,m,t)≤Zmax(i,m,t)
(3)
式中,Zmin(i,m,t)、Zmax(i,m,t)分別為i水庫(kù)m年t時(shí)段的調(diào)度圖允許最低、最高水位,m;Z(i,m,t)為i水庫(kù)m年t時(shí)段水位,m。
(3)下泄流量約束
Qmin(i,m,t)≤Q(i,m,t)≤Qmax(i,m,t)
(4)
式中,Qmin(i,m,t)、Qmax(i,m,t)分別為i水庫(kù)m年t時(shí)段的允許最小、最大出庫(kù)流量,m3/s。
(4)電站出力約束
Nmin(i,m,t)≤N(i,m,t)≤Nmax(i,m,t)
(5)
式中,Nmin(i,m,t)、Nmax(i,m,t)分別為i水電站m年t時(shí)段的允許最小、最大出力,kW;N(i,m,t)為i水電站m年t時(shí)段的平均出力,kW。
(5)梯級(jí)上下游電站間水力聯(lián)系
I(i+1,m,t)=Q(i,m,t)+q(i+1,m,t)
(6)
式中,q(i+1,m,t)為m年t時(shí)段i水庫(kù)與i+1水庫(kù)之間的區(qū)間平均流量,m3/s。
(6)計(jì)算時(shí)段初末條件
(7)
式中,Z(i,T0)、Z(i,Tend)分別為i水庫(kù)調(diào)度期初、末水位,m;Zc(i)為指定水位,若從蓄水期初開(kāi)始模擬調(diào)度,將其設(shè)置為死水位,m。為保證模擬期間梯級(jí)水庫(kù)水量平衡,時(shí)段初、末水位應(yīng)保持一致。
建模完成后,進(jìn)行梯級(jí)水庫(kù)模擬調(diào)度,徑流式水電站按照流量控制原則運(yùn)行,在滿足各項(xiàng)約束條件下進(jìn)行計(jì)算;可調(diào)節(jié)水電站則根據(jù)調(diào)度圖進(jìn)行出力控制模擬計(jì)算,分析電站電量等指標(biāo)。結(jié)合傳統(tǒng)調(diào)度圖的使用方式[16],本文所繪調(diào)度圖運(yùn)用規(guī)則為:當(dāng)水庫(kù)時(shí)段初水位位于供水期某出力區(qū)時(shí),水電站按供水期對(duì)應(yīng)出力區(qū)出力工作;位于蓄水期某出力區(qū)時(shí),水電站按蓄水期對(duì)應(yīng)出力區(qū)出力工作。此外,模擬調(diào)度還需結(jié)合當(dāng)前水庫(kù)來(lái)水等情況,將破壞約束的調(diào)度運(yùn)行規(guī)則修正至邊界處。本文主要邊界是將調(diào)度圖最高、最低水位線設(shè)定為時(shí)段運(yùn)行水位上、下限,并控制模擬運(yùn)行中蓄水期水位不消落。因?yàn)檎{(diào)度圖的優(yōu)化實(shí)際上是通過(guò)調(diào)度結(jié)果來(lái)反饋調(diào)度圖繪制過(guò)程,因此不同的調(diào)度規(guī)則及邊界條件對(duì)最終調(diào)度圖的優(yōu)化繪制有較大的影響,采用此法繪制調(diào)度圖應(yīng)結(jié)合實(shí)際繪制目的和用途合理制定調(diào)度圖使用規(guī)則及模擬運(yùn)行條件。
考慮傳統(tǒng)調(diào)度圖繪制過(guò)程,根據(jù)水庫(kù)逆時(shí)程調(diào)節(jié)計(jì)算方法[17]可知,通常情況下,當(dāng)出力一定時(shí),入流越大,所需水庫(kù)耗水量越少,則逆序計(jì)算時(shí)所得時(shí)段初水位較低,故調(diào)度線位置越低;入流越小,則調(diào)度線位置越高。因此,若選擇豐水年徑流做代表年徑流,則繪制所得調(diào)度線相對(duì)以平、枯水年徑流輸入所繪調(diào)度線位置偏低。這使得在模擬運(yùn)行中相同水位條件下前者所繪調(diào)度圖給出的出力更大,當(dāng)來(lái)水較豐時(shí),其控制效果更優(yōu);但當(dāng)來(lái)水較枯時(shí),有蓄末蓄不滿和供水期出力破壞的風(fēng)險(xiǎn)。同理,以枯水年徑流作代表年徑流,繪制所得調(diào)度圖對(duì)于枯水年調(diào)度較優(yōu);而對(duì)于來(lái)水豐、平的年份,可能導(dǎo)致蓄水期出現(xiàn)早蓄現(xiàn)象使前期發(fā)電量過(guò)小,供水期則消落偏慢。相關(guān)分析見(jiàn)圖3。當(dāng)徑流一定,出力變化時(shí)有一致結(jié)論。

圖3 豐、枯徑流逆推所得調(diào)度圖調(diào)度結(jié)果示意示意
由以上分析可知,對(duì)于同一逆推出力值,不同輸入徑流所得調(diào)度線不同;而對(duì)于相同的徑流過(guò)程,選擇的逆推出力值不同,得到的出力區(qū)也不同;由此而得的調(diào)度圖在運(yùn)用到模擬運(yùn)行中所得調(diào)節(jié)計(jì)算結(jié)果也有差異。由于繪制及調(diào)度過(guò)程的非線性及影響因素的復(fù)雜性,該差異不易給出明確的數(shù)學(xué)解析關(guān)系式。因此,針對(duì)此類(lèi)參數(shù)優(yōu)選問(wèn)題,引入啟發(fā)式優(yōu)化搜索算法開(kāi)展相關(guān)研究,在滿足保證出力且水庫(kù)蓄末蓄滿的條件下,尋求使梯級(jí)電站多年平均發(fā)電量最大的水庫(kù)發(fā)電調(diào)度圖。通過(guò)粒子群優(yōu)化算法同時(shí)控制調(diào)度圖繪制過(guò)程中的典型年徑流個(gè)數(shù)、所選代表年平均出力或水量經(jīng)驗(yàn)頻率、調(diào)度線出力值、出力線個(gè)數(shù)等多個(gè)參數(shù),建立梯級(jí)電站聯(lián)合發(fā)電調(diào)度情形下調(diào)度圖優(yōu)化繪制模型,求解模型并運(yùn)用優(yōu)化結(jié)果確定年調(diào)節(jié)水庫(kù)調(diào)度圖各繪制參數(shù),得到水電站優(yōu)化發(fā)電調(diào)度圖。優(yōu)化調(diào)度圖繪制方法見(jiàn)圖4。

圖4 優(yōu)化調(diào)度圖計(jì)算流程
首先,設(shè)置粒子群算法各項(xiàng)參數(shù)并進(jìn)行粒子初始化。選擇種群粒子總數(shù)M;每個(gè)粒子的維數(shù)N=10,各維變量物理意義及取值范圍見(jiàn)表1。其中,n為徑流序列總數(shù);r為蓄供水期保證出力之比最大值,可取為蓄、供水期保證出力之比;Ndg為供水期保證出力;Nsg為蓄水期保證出力;Ni為裝機(jī)容量;Nmin為等流量調(diào)節(jié)計(jì)算所得最小出力。初始化每個(gè)粒子的初始位置Xi(t)和初始速度Vi(t)。即

表1 調(diào)度圖及兩發(fā)電調(diào)度圖各參數(shù)
(8)
對(duì)每一個(gè)選擇的粒子,按照各維數(shù)值計(jì)算得到相關(guān)繪制參數(shù)。其中,典型年徑流序列個(gè)數(shù)為繪制各調(diào)度線時(shí)采用的徑流序列個(gè)數(shù);供水期代表年經(jīng)驗(yàn)頻率表示第一條典型徑流過(guò)程線的頻率,再依據(jù)排頻結(jié)果向上下各選擇一半數(shù)目的徑流序列;蓄水期典型年徑流過(guò)程線選擇方法同理;保證出力乘以蓄供水期指示出力之比可得蓄水期指示出力;供水期最大加大出力值作為供水期加大出力上限;蓄水期最大加大出力值作為蓄水期加大出力上限;供水期最小降低出力值作為供水期降低出力值下限;蓄水期最小降低出力值作為蓄水期降低出力值下限。在算法運(yùn)行期間,各參數(shù)均須在可行域內(nèi)尋優(yōu)。
其次,利用梯級(jí)水庫(kù)多年歷史入庫(kù)徑流資料,建立梯級(jí)水電站聯(lián)合發(fā)電調(diào)度模型。模型以調(diào)度圖繪制目標(biāo)電站為梯級(jí)龍頭水電站,其余電站作為被補(bǔ)償電站,展開(kāi)聯(lián)合發(fā)電調(diào)度模擬運(yùn)行,從上游水庫(kù)到下游水庫(kù)依次進(jìn)行模擬運(yùn)算,求出各電站各時(shí)段水位、出庫(kù)流量、電站出力、發(fā)電量等指標(biāo)。
最后,運(yùn)行粒子群算法進(jìn)行參數(shù)尋優(yōu)。對(duì)粒子種群X(k),按上述模擬運(yùn)行方式計(jì)算每個(gè)粒子的適應(yīng)度即對(duì)應(yīng)調(diào)度圖情形下的目標(biāo)函數(shù)值;比較粒子適應(yīng)度與當(dāng)前粒子最優(yōu)值pbest下的適應(yīng)度。當(dāng)前者比后者更優(yōu),則更新個(gè)體適應(yīng)度值,同時(shí)更新pbest位置為當(dāng)前位置。否則,不更新適應(yīng)度值和pbest位置;待所有粒子均更新后,比較每個(gè)粒子pbest下的適應(yīng)度與全局粒子最優(yōu)值gbest下的適應(yīng)度,當(dāng)前者比后者更優(yōu),則更新全局適應(yīng)度值,同時(shí)更新gbest位置。否則,不更新全局應(yīng)度值和gbest位置;根據(jù)下式更新粒子速度和位置。即
(9)

ωk=0.9-(k/Kmax)×0.5
(10)
其中,Kmax表示最大迭代次數(shù)。
未達(dá)結(jié)束條件則進(jìn)行粒子群更新,評(píng)價(jià)種群X(k+1),計(jì)算每個(gè)粒子的適應(yīng)度;計(jì)算完畢,即得所設(shè)調(diào)度規(guī)則和歷史徑流輸入下的優(yōu)化發(fā)電調(diào)度圖。
針對(duì)瀾滄江流域如美-功果橋串聯(lián)梯級(jí)水電站開(kāi)展實(shí)例研究分析。該段梯級(jí)屬瀾滄江如美電站及其以下、小灣電站以上河段,此區(qū)間水能資源集中,開(kāi)發(fā)條件優(yōu)越,共布置了12個(gè)梯級(jí),是接續(xù)“西電東送”的重要能源基地之一。如美水電站為瀾滄江上游西藏段規(guī)劃“1庫(kù)7級(jí)”中的第5個(gè)梯級(jí),如果規(guī)劃建設(shè),可作為該段龍頭及控制性工程,建成后美水電站將以發(fā)電為主、兼顧其他功能。由于如美電站水庫(kù)調(diào)節(jié)庫(kù)容大,實(shí)現(xiàn)壩址徑流進(jìn)行年調(diào)節(jié),水庫(kù)正常蓄水位在可研階段設(shè)定為2 895 m,死水位為2 815 m,水庫(kù)蓄能補(bǔ)償作用顯著。調(diào)節(jié)該河段內(nèi)梯級(jí)水電的出力過(guò)程,能提高外送電能的可靠性和穩(wěn)定性。

圖5 如美水電站傳統(tǒng)發(fā)電調(diào)度示意

圖6 如美電站優(yōu)化發(fā)電調(diào)度示意
本文選用1953年~2006年共54年徑流系列資料,以月為計(jì)算時(shí)段長(zhǎng),繪制如美電站傳統(tǒng)發(fā)電調(diào)度圖[16](見(jiàn)圖5)。根據(jù)本文所建模型以及求解方法,以如美電站為梯級(jí)調(diào)控主導(dǎo)電站,按照模型隨機(jī)初始化多組初始種群進(jìn)行優(yōu)化求解,從中選擇一組較優(yōu)解作為如美電站發(fā)電優(yōu)化調(diào)度圖如圖6所示。繪制兩調(diào)度圖的各項(xiàng)參數(shù)如表1所示。圖5和圖6中縱坐標(biāo)表示如美電站運(yùn)行水位,指導(dǎo)出力為保證出力區(qū)和各加大出力、降低出力區(qū)對(duì)應(yīng)的出力值。對(duì)比圖5和圖6可看出,傳統(tǒng)調(diào)度圖控制線控制范圍較大,即意味著控制作用相對(duì)較弱,且供蓄水期出力區(qū)出力值一致,沒(méi)能考慮徑流的豐枯特性和水庫(kù)的供蓄功能特點(diǎn);而優(yōu)化調(diào)度圖不僅考慮了供蓄水期特性,還進(jìn)一步加強(qiáng)了水庫(kù)水位控制作用,而且通過(guò)優(yōu)化算法選取了各出力區(qū)數(shù)值。由表1可知,優(yōu)化發(fā)電調(diào)度圖蓄水期指示出力、降低出力值都與供水期保證出力值相同,因此在此將其統(tǒng)一設(shè)置為保證出力值。

表2 傳統(tǒng)、優(yōu)化調(diào)度圖模擬計(jì)算梯級(jí)水能年均指標(biāo)
利用長(zhǎng)系列徑流資料展開(kāi)梯級(jí)模擬調(diào)度以檢驗(yàn)兩種調(diào)度圖,檢驗(yàn)結(jié)果見(jiàn)表2。根據(jù)傳統(tǒng)調(diào)度圖對(duì)如美電站進(jìn)行模擬調(diào)度,其余各電站以徑流式電站模式運(yùn)行,如美自身年均發(fā)電量為102.96 億kW·h,梯級(jí)電站年均發(fā)電量為679.70 億kW·h,下游梯級(jí)通過(guò)未經(jīng)調(diào)節(jié)徑流進(jìn)行模擬計(jì)算得到總電量為557.25 億kW·h,因此梯級(jí)補(bǔ)償效益為19.49 億kW·h;根據(jù)優(yōu)化調(diào)度圖通過(guò)模擬運(yùn)行計(jì)算可得,如美自身年均發(fā)電量為104.73 億kW·h,梯級(jí)年均發(fā)電量為693.99 億kW·h,相對(duì)于傳統(tǒng)調(diào)度圖,分別提高了約1.72%和2.10%,補(bǔ)償效益增加了13.51 億kW·h,且兩種情況下,如美電站90%保證率對(duì)應(yīng)保證出力675 MW均能被滿足。說(shuō)明在滿足保證率要求的情形下,優(yōu)化調(diào)度圖提高了電站自身及梯級(jí)電站的發(fā)電效益。兩種方案下,如美水庫(kù)的各時(shí)段多年平均水位和平均出力過(guò)程如圖7所示。

圖7 兩種方案下如美電站多年平均運(yùn)行水位和出力過(guò)程
由表2和圖7可知:①按照傳統(tǒng)調(diào)度圖模擬運(yùn)行,前期如美電站按保證出力運(yùn)行,由于水位最高控制線較高,防棄水作用差,導(dǎo)致其水位升高過(guò)快,造成蓄水期前期發(fā)電量較小,后期雖然因高水位使得發(fā)電量增加,但是不足以彌補(bǔ)前期電量過(guò)小的情形;②按照優(yōu)化調(diào)度圖模擬運(yùn)行,如美電站蓄水期運(yùn)行水位相比傳統(tǒng)調(diào)度圖模擬運(yùn)行結(jié)果有所降低,原因是在優(yōu)化調(diào)度圖的控制作用下,為減少水庫(kù)棄水量,增加前期電量,水庫(kù)需平穩(wěn)地增長(zhǎng)水位,使水庫(kù)調(diào)豐作用得以更充分發(fā)揮,增加電站發(fā)電效益,實(shí)際上優(yōu)化調(diào)度線的最高水位線起到了防棄水線的作用;③兩種調(diào)度圖運(yùn)行模式在供水期調(diào)度結(jié)果差別不大,是因?yàn)楣┧趶搅餍蛄芯怠⒎讲钶^小,兩種調(diào)度圖均能較好地發(fā)揮水位控制作用,而優(yōu)化調(diào)度圖運(yùn)行水位略高,是為了提高電站水頭以增大發(fā)電量。由此可以看出,對(duì)于年調(diào)節(jié)水庫(kù),調(diào)度圖優(yōu)化的主要目的在于優(yōu)化蓄水期的水位控制線,在保證電站蓄末蓄滿的情況下使蓄水期各時(shí)段出力相對(duì)均勻,增大前期發(fā)電量。因此優(yōu)化調(diào)度圖的水位線具有較好的水位控制作用,可為實(shí)際調(diào)度運(yùn)行提供參考。
本文結(jié)合以往調(diào)度圖改進(jìn)和優(yōu)化繪制方式,充分考慮傳統(tǒng)調(diào)度圖的優(yōu)點(diǎn),引入智能優(yōu)化方法,建立了一種年調(diào)節(jié)水庫(kù)發(fā)電調(diào)度圖的多參數(shù)優(yōu)選繪制方法,以瀾滄江流域如美至功果橋串聯(lián)梯級(jí)水庫(kù)為研究對(duì)象,考慮梯級(jí)水電站之間的徑流聯(lián)系,建立梯級(jí)聯(lián)合調(diào)度模型,以梯級(jí)多年平均發(fā)電量最大兼顧保證出力和蓄末蓄滿為目標(biāo),進(jìn)行調(diào)度圖優(yōu)化。通過(guò)計(jì)算分析比較兩種調(diào)度圖調(diào)度方案下梯級(jí)電站發(fā)電效益及運(yùn)行特點(diǎn),得到以下結(jié)論:
(1) 依據(jù)傳統(tǒng)方法繪制水庫(kù)發(fā)電調(diào)度圖主觀性強(qiáng),繪制過(guò)程參數(shù)選取依賴人工經(jīng)驗(yàn),易得到調(diào)度結(jié)果不佳的情況,而本文提出的優(yōu)化方法可以根據(jù)徑流資料等客觀數(shù)據(jù)得到優(yōu)化發(fā)電調(diào)度圖,且經(jīng)實(shí)例驗(yàn)證,水位控制作用較好,說(shuō)明該方法具備一定潛在實(shí)用價(jià)值。
(2) 如美梯級(jí)水庫(kù)發(fā)電模擬調(diào)度表明,實(shí)測(cè)徑流資料模擬情形下,按優(yōu)化調(diào)度圖運(yùn)行可將梯級(jí)年均發(fā)電量提高約2%,優(yōu)化效果較好。表明該方法在滿足保證率的基礎(chǔ)上,能夠發(fā)揮梯級(jí)水電站的補(bǔ)償調(diào)節(jié)作用,進(jìn)一步發(fā)揮梯級(jí)水電站的發(fā)電效益。
(3) 本文尚未針對(duì)算法的收斂問(wèn)題進(jìn)行研究,而是通過(guò)多次初始化種群進(jìn)行計(jì)算并選擇了較合理的結(jié)果作為分析對(duì)象,雖然從所得解中已分析出調(diào)度圖優(yōu)化方法的有關(guān)規(guī)律,但不能完全判斷出全局最優(yōu)解,因此相關(guān)工作仍需進(jìn)一步開(kāi)展研究。