一種梯級(jí)水電站短期優(yōu)化調(diào)度方法
【專利摘要】本發(fā)明公開(kāi)了一種梯級(jí)水電站短期優(yōu)化調(diào)度方法。本發(fā)明提出一種全新的適用于梯級(jí)水電短期調(diào)度棄能消納的“平土”算法。本發(fā)明采用梯級(jí)和廠內(nèi)兩級(jí)協(xié)調(diào)控制的技術(shù)：對(duì)梯級(jí)各個(gè)水電站進(jìn)行計(jì)算得到各水庫(kù)的出庫(kù)流量和發(fā)電水頭；再以此為基礎(chǔ)，利用混合整數(shù)規(guī)劃，分別確定各個(gè)電廠的最優(yōu)開(kāi)停機(jī)計(jì)劃；進(jìn)而采用動(dòng)態(tài)規(guī)劃按照給定水量使發(fā)電量最大的目標(biāo)對(duì)機(jī)組群進(jìn)行負(fù)荷優(yōu)化分配。同時(shí)，得到各水電站的發(fā)電流量過(guò)程之后，又可以將結(jié)果反饋給梯級(jí)系統(tǒng)，梯級(jí)系統(tǒng)以此為約束條件繼續(xù)重復(fù)利用“平土”算法進(jìn)行迭代計(jì)算。本發(fā)明可以顯著緩解大規(guī)模水電系統(tǒng)的棄能矛盾，并能提高廠內(nèi)機(jī)組運(yùn)行與調(diào)度水平。
【專利說(shuō)明】一種梯級(jí)水電站短期優(yōu)化調(diào)度方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于水文學(xué)與運(yùn)籌學(xué)的交叉【技術(shù)領(lǐng)域】，更具體地，涉及一種梯級(jí)水電站短 期優(yōu)化調(diào)度方法。

【背景技術(shù)】
[0002] 水電系統(tǒng)短期優(yōu)化調(diào)度可以認(rèn)為是一個(gè)確定性優(yōu)化問(wèn)題，求解短期調(diào)度問(wèn)題有傳 統(tǒng)的數(shù)學(xué)規(guī)劃方法和智能算法。智能算法在求解該問(wèn)題時(shí)存在最優(yōu)解不穩(wěn)定的不足。目前， 尚無(wú)利用智能算法成功解決大規(guī)模水電系統(tǒng)短期優(yōu)化的調(diào)度問(wèn)題的實(shí)例。因此，當(dāng)前短期 水電調(diào)度問(wèn)題主要應(yīng)用傳統(tǒng)的線性規(guī)劃或者動(dòng)態(tài)規(guī)劃等方法求解。然而，利用傳統(tǒng)方法求 解需要對(duì)約束條件進(jìn)行線性化處理，致使與原問(wèn)題有所偏差，對(duì)于凸規(guī)劃的問(wèn)題，多采用分 段線性規(guī)劃的處理以減小誤差。然而這些方法在水庫(kù)群優(yōu)化調(diào)度尋找最優(yōu)解的過(guò)程中，容 易過(guò)早的陷入局部最優(yōu)解，導(dǎo)致不能找到理想的最優(yōu)解。
[0003] 此外，在廠內(nèi)機(jī)組經(jīng)濟(jì)運(yùn)行中，等微增率和動(dòng)態(tài)規(guī)劃被廣泛用于機(jī)組群的負(fù)荷分 配，但兩種方法都未能考慮開(kāi)、停機(jī)的時(shí)間耦合約束，應(yīng)用受到很大限制。國(guó)外有報(bào)道采用 開(kāi)、停機(jī)水量損失來(lái)限制頻繁開(kāi)停機(jī)運(yùn)行，但這種水量損失往往難于估算，更為普遍和實(shí)用 的做法是通過(guò)約束開(kāi)停機(jī)持續(xù)時(shí)間和次數(shù)來(lái)提高機(jī)組運(yùn)行的安全和穩(wěn)定性。即便如此，由 于問(wèn)題的離散、非凸和非線性特征，對(duì)于規(guī)模稍大的水電機(jī)組群的聯(lián)合開(kāi)停機(jī)和負(fù)荷優(yōu)化 問(wèn)題，求解起來(lái)也非常困難。
[0004] 從目前國(guó)內(nèi)外研究的發(fā)展趨勢(shì)看，如何考慮時(shí)間和空間約束的耦合，協(xié)調(diào)控制梯 級(jí)水電站調(diào)度和廠內(nèi)機(jī)組群負(fù)荷與開(kāi)停機(jī)優(yōu)化、并將其納入到一體化調(diào)度和控制體系中， 是進(jìn)一步亟待解決的難題。
[0005] 從報(bào)道情況看，發(fā)達(dá)國(guó)家更為關(guān)注梯級(jí)水電調(diào)度涉及生態(tài)、環(huán)境、社會(huì)的影響，對(duì) 水電棄能消納技術(shù)的研究鮮有報(bào)道。在國(guó)內(nèi)，各大電力、電網(wǎng)公司都相繼開(kāi)展了梯級(jí)水電調(diào) 度高級(jí)應(yīng)用系統(tǒng)的開(kāi)發(fā)，但由于問(wèn)題復(fù)雜、涉及面廣，實(shí)際上所開(kāi)發(fā)的系統(tǒng)在生產(chǎn)實(shí)用上還 有相當(dāng)大差距。從目前研究情況看，無(wú)論是國(guó)內(nèi)還是國(guó)外，梯級(jí)水電和廠內(nèi)經(jīng)濟(jì)運(yùn)行一般都 采用逐級(jí)控制的方式，還達(dá)到兩級(jí)一體化協(xié)調(diào)控制的水平。


【發(fā)明內(nèi)容】

[0006] 針對(duì)現(xiàn)有技術(shù)的以上缺陷或改進(jìn)需求，本發(fā)明以水電棄能消納為導(dǎo)向，模擬來(lái)回 推拉平整土地，研究出了一種全新的"平土"算法，為水庫(kù)水電站群聯(lián)合優(yōu)化調(diào)度跳出局部 最優(yōu)或者可行化修正提供一條簡(jiǎn)單、實(shí)用且高效的局部修正算法。采用梯級(jí)和廠內(nèi)兩級(jí)協(xié) 調(diào)控制技術(shù)，同時(shí)結(jié)合混合整數(shù)規(guī)劃和動(dòng)態(tài)規(guī)劃，克服大規(guī)模水電機(jī)組群開(kāi)停機(jī)和負(fù)荷優(yōu) 化的維數(shù)難題，最終獲得水電調(diào)度梯級(jí)和廠內(nèi)機(jī)組經(jīng)濟(jì)運(yùn)行一體化的發(fā)電計(jì)劃。
[0007] 本發(fā)明提供一種梯級(jí)水電站短期優(yōu)化調(diào)度方法，其總體思路是：對(duì)水電梯級(jí)和廠 內(nèi)機(jī)組群分為兩級(jí)協(xié)調(diào)控制，包括以下步驟：
[0008] 步驟1根據(jù)徑流資料以及梯級(jí)水電站的參數(shù)，建立所述梯級(jí)水電站短期優(yōu)化調(diào)度 的多目標(biāo)優(yōu)化模型；
[0009] 步驟2求解所述多目標(biāo)優(yōu)化模型的最優(yōu)解，獲得多個(gè)水庫(kù)的出庫(kù)流量和發(fā)電水頭 過(guò)程；
[0010] 步驟3將獲得的所述出庫(kù)流量和所述發(fā)電水頭過(guò)程作為廠內(nèi)水電機(jī)組群的計(jì)算 邊界條件，建立混合整數(shù)規(guī)劃模型，獲得所述廠內(nèi)水電機(jī)組群的最優(yōu)開(kāi)停機(jī)計(jì)劃；
[0011] 步驟4采用動(dòng)態(tài)規(guī)劃按照給定水量發(fā)電量最大化的目標(biāo)對(duì)所述廠內(nèi)水電機(jī)組群 進(jìn)行負(fù)荷優(yōu)化分配；
[0012] 步驟5根據(jù)獲得的所述廠內(nèi)水電機(jī)組群優(yōu)化運(yùn)行結(jié)果更新整個(gè)水電站的發(fā)電流 量上下限，反饋給上級(jí)系統(tǒng)，梯級(jí)水電調(diào)度以此為約束條件，再執(zhí)行所述步驟2,如此反復(fù)迭 代直到整個(gè)梯級(jí)水電站的棄水減小到最小和廠內(nèi)發(fā)電效益最大，獲得一體化梯級(jí)和廠內(nèi)最 優(yōu)發(fā)電計(jì)劃。
[0013] 總體而言，通過(guò)本發(fā)明所構(gòu)思的以上技術(shù)方案與現(xiàn)有技術(shù)相比，具有以下有益效 果：
[0014] 1、本發(fā)明基于成熟的數(shù)學(xué)規(guī)劃理論，不同于現(xiàn)有水電調(diào)度技術(shù)多以經(jīng)驗(yàn)調(diào)度為 主，可以應(yīng)用于實(shí)際生產(chǎn)調(diào)度中，同時(shí)克服了現(xiàn)有的諸多智能算法的調(diào)度方案生成的最優(yōu) 解不穩(wěn)定、只停留在理論層面無(wú)法應(yīng)用到實(shí)際調(diào)度生產(chǎn)中的問(wèn)題，具有實(shí)用性；
[0015] 2、本發(fā)明提出了一種全新的"平土"算法，為水電站優(yōu)化調(diào)度避免陷入局部最優(yōu)解 陷阱提供了一條全新的途徑，為當(dāng)前水電站實(shí)際調(diào)度提供理論支撐和技術(shù)支撐，提高水電 站經(jīng)濟(jì)效益；
[0016] 3、本發(fā)明采用梯級(jí)與電廠兩級(jí)調(diào)度、電廠內(nèi)開(kāi)停機(jī)和機(jī)組負(fù)荷優(yōu)化的兩階段優(yōu) 化，以及混合整數(shù)規(guī)劃和動(dòng)態(tài)規(guī)劃的應(yīng)用，克服了大規(guī)模水電調(diào)度的維數(shù)災(zāi)問(wèn)題。

【專利附圖】

【附圖說(shuō)明】
[0017] 圖1為本發(fā)明梯級(jí)水電站短期優(yōu)化調(diào)度方法的流程圖；
[0018] 圖2為本發(fā)明梯級(jí)水電棄能"平土"消納示意圖；
[0019] 圖3為本發(fā)明梯級(jí)水庫(kù)"平土"算法"前推"步驟流程圖。

【具體實(shí)施方式】
[0020] 為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合附圖及實(shí)施例，對(duì) 本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說(shuō)明。應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅僅用以解釋本發(fā)明，并 不用于限定本發(fā)明。此外，下面所描述的本發(fā)明各個(gè)實(shí)施方式中所涉及到的技術(shù)特征只要 彼此之間未構(gòu)成沖突就可以相互組合。
[0021] 圖1所示為本發(fā)明梯級(jí)水電站短期優(yōu)化調(diào)度方法的流程圖，其總體思路是：對(duì)水 電梯級(jí)和廠內(nèi)機(jī)組群分為兩級(jí)協(xié)調(diào)控制。
[0022] 步驟1根據(jù)徑流資料以及梯級(jí)水電站的參數(shù)，建立梯級(jí)電站短期優(yōu)化調(diào)度的數(shù)學(xué) 模型；
[0023] 步驟2進(jìn)行梯級(jí)水電棄能消納的計(jì)算，獲得各水庫(kù)的出庫(kù)流量過(guò)程和發(fā)電水頭過(guò) 程，在本發(fā)明實(shí)施例中，利用"平土"算法進(jìn)行梯級(jí)水電棄能消納的計(jì)算；
[0024] 步驟3將所獲得的出庫(kù)流量過(guò)程和發(fā)電水頭過(guò)程作為廠內(nèi)水電機(jī)組群的計(jì)算邊 界條件，建立混合整數(shù)規(guī)劃模型，獲得機(jī)組最優(yōu)開(kāi)停機(jī)計(jì)劃；
[0025] 步驟4采用動(dòng)態(tài)規(guī)劃按照給定水量發(fā)電量最大化的目標(biāo)對(duì)機(jī)組群進(jìn)行負(fù)荷優(yōu)化 分配；
[0026] 步驟5廠內(nèi)機(jī)組群優(yōu)化運(yùn)行結(jié)果可更新整個(gè)水電站的發(fā)電流量上下限，反饋給上 級(jí)系統(tǒng)，梯級(jí)水電調(diào)度以此為約束條件，再執(zhí)行步驟2,如此反復(fù)迭代，直到整個(gè)梯級(jí)水電站 的棄水減小到最小，從而獲得一體化梯級(jí)和廠內(nèi)最優(yōu)發(fā)電計(jì)劃，即梯級(jí)各水電廠蓄放水計(jì) 劃、機(jī)組開(kāi)停機(jī)計(jì)劃及機(jī)組出力過(guò)程，流程結(jié)束。
[0027] 圖2所示為本發(fā)明梯級(jí)水電棄能"平土"消納示意圖，圖2將結(jié)合圖1進(jìn)行描述。 本發(fā)明梯級(jí)水電棄能消納控制可分為"前推"、"后拉"、再"前推"三步進(jìn)行。前、后向修正時(shí) 需要同時(shí)優(yōu)化調(diào)節(jié)多個(gè)水庫(kù)的時(shí)段末或初蓄水，目標(biāo)是以偏離每個(gè)優(yōu)化時(shí)段的初始蓄水過(guò) 程盡可能小的方式使該時(shí)段棄水/棄能最小。當(dāng)考慮水庫(kù)間的水流滯時(shí)情況時(shí)，將上、下游 水庫(kù)的調(diào)節(jié)時(shí)段之差設(shè)定為水庫(kù)間水流滯時(shí)長(zhǎng)，消除了出庫(kù)流量對(duì)下階段調(diào)節(jié)的影響。在 本發(fā)明實(shí)施例中，水庫(kù)結(jié)構(gòu)如圖2右側(cè)所示，水庫(kù)1與水庫(kù)2為并聯(lián)的第一級(jí)水庫(kù)，水庫(kù)3、 4依次為串聯(lián)的第二、三級(jí)水庫(kù)。在本發(fā)明實(shí)施例中，以由4個(gè)串并聯(lián)水庫(kù)組成的梯級(jí)水電 站為例進(jìn)行說(shuō)明，但不以此為限。
[0028] 在上述步驟1中，將調(diào)度期內(nèi)的時(shí)間長(zhǎng)度分為T個(gè)時(shí)段，根據(jù)調(diào)度期棄水最小和迭 代量與原始解偏差最小的原則，建立梯級(jí)水庫(kù)短期調(diào)度多目標(biāo)優(yōu)化模型如下：
[0029]

【權(quán)利要求】
1. 一種梯級(jí)水電站短期優(yōu)化調(diào)度方法，其特征在于，包括以下步驟： 步驟1根據(jù)徑流資料以及梯級(jí)水電站的參數(shù)，建立所述梯級(jí)水電站短期優(yōu)化調(diào)度的多 目標(biāo)優(yōu)化模型； 步驟2求解所述多目標(biāo)優(yōu)化模型的最優(yōu)解，獲得多個(gè)水庫(kù)的出庫(kù)流量和發(fā)電水頭過(guò) 程； 步驟3將獲得的所述出庫(kù)流量和所述發(fā)電水頭過(guò)程作為廠內(nèi)水電機(jī)組群的計(jì)算邊界 條件，建立混合整數(shù)規(guī)劃模型，獲得所述廠內(nèi)水電機(jī)組群的最優(yōu)開(kāi)停機(jī)計(jì)劃； 步驟4采用動(dòng)態(tài)規(guī)劃按照給定水量發(fā)電量最大化的目標(biāo)對(duì)所述廠內(nèi)水電機(jī)組群進(jìn)行 負(fù)荷優(yōu)化分配； 步驟5根據(jù)獲得的所述廠內(nèi)水電機(jī)組群優(yōu)化運(yùn)行結(jié)果更新整個(gè)水電站的發(fā)電流量上 下限，反饋給上級(jí)系統(tǒng)，梯級(jí)水電調(diào)度以此為約束條件，再執(zhí)行所述步驟2,如此反復(fù)迭代直 到整個(gè)梯級(jí)水電站的棄水減小到最小和廠內(nèi)發(fā)電效益最大，最終獲得一體化梯級(jí)和廠內(nèi)最 優(yōu)發(fā)電計(jì)劃。
2. 如權(quán)利要求1所述的方法，其特征在于，所述步驟1將調(diào)度期內(nèi)的時(shí)間長(zhǎng)度分為T個(gè) 時(shí)段，根據(jù)調(diào)度期棄水最小和迭代量與原始解偏差最小的原則，建立所述多目標(biāo)優(yōu)化模型 如下：
其中，i，t分別表示第i個(gè)水庫(kù)和第t個(gè)時(shí)間段數(shù)；N，T分別表示所述調(diào)度期內(nèi)的水庫(kù) 數(shù)量和時(shí)段數(shù)；split表示所述第i個(gè)水庫(kù)在所述第t個(gè)時(shí)間段的棄水流量；vit表示所述第 i個(gè)水庫(kù)在所述時(shí)段t初的庫(kù)容；vf表示所述第i個(gè)水庫(kù)在所述時(shí)段t初的庫(kù)容初始解； β表示與庫(kù)容相關(guān)的系數(shù)，是水庫(kù)興利庫(kù)容與多年平均來(lái)水量之比。
3. 如權(quán)利要求2所述的方法，其特征在于，所述多目標(biāo)優(yōu)化模型需滿足以下約束條件： 庫(kù)容約束：vfSv, 出庫(kù)流量限制：(Τ4 0^(ΤΧ 發(fā)電流量限制：(4 )<I W (Λ") 水量平衡方程
其中，ν it表示時(shí)段t初水庫(kù)i的庫(kù)容，初始庫(kù)容以=if，為所述水庫(kù)i在所述調(diào) 度期起始時(shí)間觀察到的水庫(kù)庫(kù)容；f分別表示所述時(shí)段t初水庫(kù)i的庫(kù)容最大和最 小值；Qit為所述水庫(kù)i在所述時(shí)段t內(nèi)的實(shí)際下泄流量，Qit = qit+split，0廣、βΤ分別表 示所述水庫(kù)i在所述時(shí)段t內(nèi)的最小和最大下泄流量值；qit表示所述水庫(kù)i在所述時(shí)段t 的發(fā)電流量；分別表示所述水庫(kù)i在固定水頭下最小和最大發(fā)電流量， 是關(guān)于水頭的函數(shù)；hit表示所述水庫(kù)i在所述時(shí)段t的平均水頭；Ω (i)表示與所述水庫(kù)i 直接相連的水庫(kù)集合；τ k表示水庫(kù)k的水流到達(dá)與其直接相連的下一個(gè)水庫(kù)的水流時(shí)間； Iit表示所述水庫(kù)i在所述時(shí)段t的當(dāng)?shù)厝肓髁?；?t表示時(shí)段t的間隔時(shí)長(zhǎng)。
4. 如權(quán)利要求1所述的方法，其特征在于，所述步驟2中利用"平土"算法求解所述多 目標(biāo)優(yōu)化模型的最優(yōu)解，其中，所述"平土"算法具體包括以下子步驟： (2-1) "前推"過(guò)程，在T = 1，. . .，T-1+ τ _各時(shí)間段分別對(duì)梯級(jí)水庫(kù)進(jìn)行向前修正， 其中= ，u⑴是第i個(gè)水庫(kù)水流到達(dá)最后一個(gè)水庫(kù)的流達(dá)時(shí)間，時(shí)段t- u⑴ 末水庫(kù)i的庫(kù)容由初始化，然后通過(guò)迭代不斷更新t- u (i)時(shí)段末水庫(kù)i的 庫(kù)容，當(dāng)?shù)弥禎M足J時(shí)，則將得到最優(yōu)庫(kù)容其中δ為預(yù)先 設(shè)定的正值，可根據(jù)經(jīng)驗(yàn)確定迭代后與迭代初的偏差值S。求解最優(yōu)庫(kù)容<"(_的線性規(guī) 劃問(wèn)題表達(dá)為以下數(shù)學(xué)模型：
其中，時(shí)段t- u (i)的水庫(kù)初庫(kù)容表示為：;時(shí)段t滿足 u (i)彡t彡T+u (i) ; α為與偏差相關(guān)的權(quán)重系數(shù)，根據(jù)水庫(kù)調(diào)節(jié)性能取值，范圍為〇? 1. 0之間，水庫(kù)調(diào)節(jié)性能越大取值越大；I/., 和1 是水庫(kù)i在時(shí)段t初與目標(biāo)庫(kù)容 的正負(fù)偏差值； 得到t-u (i)時(shí)段末水庫(kù)的庫(kù)容值之后，時(shí)段t-u (i)的各水庫(kù)的下泄流量 ?2^ i j v{ j j 通過(guò)水量平衡方程求出，再根據(jù)公式
， 求出水庫(kù)在時(shí)段t- U Q)的水頭& ，其中，時(shí)段t- U (i)的庫(kù)容為： 9 ^/,Ι-υ?Ι)
?式中的m 在時(shí)段t-1計(jì)算得到； (2-2) "回拉"過(guò)程，在各時(shí)段Τ+τ_、Τ-1+τ_，…，2向后修正，根據(jù)已知的時(shí)段 t_ υ⑴末的庫(kù)容值，對(duì)時(shí)段t- u⑴初的庫(kù)容進(jìn)行優(yōu)化，其中u (i)彡t彡Τ+ u (i); (2-3)重復(fù)所述步驟（2-1)的"前推"過(guò)程，在各時(shí)段t = 0,1，…，T-1+τ _向前修 正，在t = T+u (i)時(shí)段需要修正所述調(diào)度期末的庫(kù)容值以減少棄水。
5. 如權(quán)利要求4所述的方法，其特征在于，所述步驟（2-1)中求解所述最優(yōu)庫(kù)容 丨1的約束條件為： 目標(biāo)庫(kù)容的偏差約束和最大最小值約束：
出庫(kù)流量限制：β=( 發(fā)電流量限制
水量平衡方程
其中，At表不時(shí)段t的間隔時(shí)長(zhǎng)。
6. 如權(quán)利要求1所述的方法，其特征在于，所述步驟3中所述廠內(nèi)水電機(jī)組群的最優(yōu)開(kāi) 停機(jī)計(jì)劃即為求以下數(shù)學(xué)模型的最優(yōu)解：
其約束條件為：
以及最小開(kāi)停機(jī)持續(xù)時(shí)間和最多開(kāi)停機(jī)次數(shù)限制， 其中，-if1為0-1變量，表示j機(jī)組在t時(shí)段是否運(yùn)行在k運(yùn)行區(qū)，運(yùn)行在k區(qū)則為1否 則為〇 ;τ為控制期離散時(shí)刻數(shù)；splt為所述t時(shí)段i水庫(kù)的棄水流量；α為權(quán)重系數(shù)，根據(jù) 水庫(kù)調(diào)節(jié)性能取值，范圍為〇?1. 〇之間，水庫(kù)調(diào)節(jié)性能越大取值越大；為〇-1變量， 表示j機(jī)組在所述t時(shí)段是否運(yùn)行在效率最高運(yùn)行區(qū)，運(yùn)行在效率最高區(qū)則為1，否則為0 ; M(i)為i水庫(kù)的機(jī)組數(shù)；q#為j機(jī)組在所述t時(shí)段的發(fā)電流量；Q#為第i個(gè)水庫(kù)的出庫(kù)流 量；K為第i個(gè)水庫(kù)水電站的機(jī)組總數(shù)；和分別表示j機(jī)組在t時(shí)段k運(yùn)行區(qū) 的發(fā)電流量上、下限。
7. 如權(quán)利要求1或6所述的方法，其特征在于，所述步驟4對(duì)所述廠內(nèi)水電機(jī)組群進(jìn)行 負(fù)荷優(yōu)化即為求以下數(shù)學(xué)模型的最優(yōu)解：
其約束條件為：
其中，為表示j機(jī)組在t時(shí)段最優(yōu)開(kāi)機(jī)計(jì)劃；p#為j機(jī)組在t時(shí)段的分配負(fù)荷。
【文檔編號(hào)】G06Q50/06GK104123589SQ201410289145
【公開(kāi)日】2014年10月29日 申請(qǐng)日期:2014年6月24日 優(yōu)先權(quán)日:2014年6月24日
【發(fā)明者】王金文, 劉雙全, 陳誠(chéng), 康傳雄 申請(qǐng)人:華中科技大學(xué), 云南電力調(diào)度控制中心