陳祺煒， 吳 熙， 李海峰， 羅凱明， 劉 林

(1. 東南大學電氣工程學院，江蘇 南京 210096;2. 國網(wǎng)江蘇省電力公司調(diào)度控制中心，江蘇 南京 210024)

孤網(wǎng)全黑狀態(tài)的恢復網(wǎng)架和分區(qū)優(yōu)化算法

陳祺煒1， 吳 熙1， 李海峰2， 羅凱明2， 劉 林2

(1. 東南大學電氣工程學院，江蘇 南京 210096;2. 國網(wǎng)江蘇省電力公司調(diào)度控制中心，江蘇 南京 210024)

電網(wǎng)故障后解列出的孤網(wǎng)是電力系統(tǒng)的重要組成部分，優(yōu)先恢復系統(tǒng)中的孤立電網(wǎng)具有簡化計算、加速電網(wǎng)恢復和及時帶動重要負荷的優(yōu)點。文中針對恢復的快速性和穩(wěn)定性，綜合考慮了孤網(wǎng)全黑狀態(tài)的分區(qū)策略和各分區(qū)的恢復網(wǎng)架，引入最短啟動時間算法和整數(shù)條件下的線性規(guī)劃法，提出一種全新的孤網(wǎng)分區(qū)優(yōu)化算法。該算法以機組、負荷最短啟動時間為目標函數(shù)，以分區(qū)內(nèi)的機組穩(wěn)態(tài)有功、無功功率平衡為約束，將復雜的多目標非線性規(guī)劃問題化簡為單目標多約束的線性規(guī)劃問題。文末以新英格蘭10機39節(jié)點系統(tǒng)作為算例驗證本文算法的有效性。

黑啟動；孤網(wǎng)分區(qū)；并行恢復；恢復網(wǎng)架

隨著經(jīng)濟的快速發(fā)展，交直流輸電和電網(wǎng)互聯(lián)技術(shù)不斷進步。系統(tǒng)中大量引入各種高壓電氣設(shè)備、大容量機組，導致現(xiàn)代電力系統(tǒng)日趨復雜化。分布式電源等接入電網(wǎng)供應(yīng)電源的同時，電網(wǎng)抵抗大擾動能力卻在不斷被弱化。大規(guī)模電力系統(tǒng)中，當局部電網(wǎng)發(fā)生嚴重故障時，一旦調(diào)度人員操作不當或保護裝置不正確動作，將會引發(fā)電力系統(tǒng)的連鎖反應(yīng)，擴大故障影響，導致大范圍停電甚至是全網(wǎng)的崩潰[1]。如果能在大停電事故時，采取合理快捷的黑啟動措施，及時排除故障，減小對經(jīng)濟、社會、民生的不良影響，具有非常重要的意義。

電網(wǎng)故障后繼電保護相繼動作，解列出的多個孤網(wǎng)成為電力系統(tǒng)的主要組成部分。孤網(wǎng)黑啟動就是指在全黑狀態(tài)下，不依靠大電網(wǎng)幫助，僅通過孤網(wǎng)內(nèi)部的黑啟動機組帶動待啟動機組，逐步實現(xiàn)孤網(wǎng)的恢復[2]。孤網(wǎng)黑啟動的恢復策略主要分為兩類[3]：串行恢復和并行恢復。串行恢復對孤網(wǎng)內(nèi)的待啟動機組依次恢復，屬于接力模式，過程冗雜繁瑣，恢復速度慢。并行恢復在技術(shù)允許的情況下將孤網(wǎng)的復雜機組分解為若干個分區(qū)同時黑啟動。在分區(qū)內(nèi)部，一般先恢復待啟動機組和輸電線路，形成特定的恢復網(wǎng)架，加快全網(wǎng)恢復進程。本文的黑啟動分區(qū)方案將建立在并行恢復的基礎(chǔ)上。

電力系統(tǒng)黑啟動問題一直是電力系統(tǒng)安全防御的重要課題之一，備受中外學者關(guān)注。文獻[4]提出了一種基于有序二元決策圖的分區(qū)方法，將分區(qū)模型轉(zhuǎn)換為布爾函數(shù)。文獻[5]提出了一種以節(jié)點恢復成功率為目標的分區(qū)優(yōu)化方案，并采用交叉粒子群優(yōu)化算法進行求解，提高了黑啟動成功率。文獻[6]考慮了黑啟動過程中設(shè)備投運失敗的可能性，建立了以路徑轉(zhuǎn)移系數(shù)為基礎(chǔ)的評價函數(shù)，以此對電力系統(tǒng)進行分區(qū)。文獻[7]采用圖論的傳統(tǒng)方法對電力系統(tǒng)拓撲進行優(yōu)化分區(qū)。文獻[8]則采用節(jié)點電壓相近度進行作為黑啟動機組分區(qū)依據(jù)，根據(jù)潮流計算的節(jié)點電壓、黑啟動電源的數(shù)目以及黑啟動電源所在節(jié)點的電壓獲得電壓閾值進行分區(qū)。文獻[9]以恢復電網(wǎng)中盡可能多的發(fā)電能力為目標，采用遺傳算法求解電網(wǎng)黑啟動分區(qū)。

可見，常規(guī)黑啟動方案只研究單個黑啟動機組的恢復方法，無法同時啟動多臺機組，浪費了啟動時間[10]。此外，近年來不斷引入智能算法求解電力系統(tǒng)的恢復網(wǎng)架及分區(qū)方案，存在計算量較大，結(jié)果穩(wěn)定性較差等問題，難以應(yīng)用于工程實踐。對此，文中綜合考慮孤網(wǎng)恢復網(wǎng)架和分區(qū)策略，以恢復時間最短為目標函數(shù)，使盡可能多的機組在短時間內(nèi)恢復，并使得每個分區(qū)都有且只有一個黑啟動機組，以該機組為核心構(gòu)建孤網(wǎng)各分區(qū)的恢復網(wǎng)架，為黑啟動后期的網(wǎng)架重構(gòu)奠定堅實基礎(chǔ)。

1 系統(tǒng)建模

1.1 目標函數(shù)

孤網(wǎng)內(nèi)含有黑啟動機組、待啟動機組和負荷三類節(jié)點。分區(qū)可以將這三類節(jié)點進行合理的劃分和規(guī)劃，實現(xiàn)以下3個主要目標：(1) 以最短時間恢復電力系統(tǒng)中的重要輸電線路。(2) 以最短時間恢復電網(wǎng)的發(fā)電能力。(3) 以最短時間建立恢復網(wǎng)架。因此需要將相關(guān)線路所需的恢復時間設(shè)定為線路權(quán)重(詳見1.3)，使用最短恢復時間算法對線路加權(quán)進行計算，確定黑啟動電源到任意節(jié)點的最短恢復時間。綜上所述，以任意一個待啟動機組到黑啟動機組所需的恢復時間最短為目標，待啟動機組最優(yōu)分組問題的目標函數(shù)可以描述為：

(1)

其中，m為待啟動機組數(shù)量； [s1g，s2g，…，sng]T表示編號為g的待啟動機組是否屬于黑啟動機組的布爾變量。例如，當scg=0(1

1.2 約束條件

黑啟動分區(qū)策略必需保證后續(xù)系統(tǒng)恢復過程中的安全穩(wěn)定，因此在分區(qū)過程中需要注意相關(guān)的約束條件。具體包含如下三部分：

(1)黑啟動機組約束，即黑啟動機組需要分入不同的分區(qū)中，此外，待啟動機組需要平均分入不同的分區(qū)中，保證各個分區(qū)的機組容量平衡。

(2)

(2)功率平衡約束，各分區(qū)內(nèi)的黑啟動機組容量、待啟動機組容量以及重要負荷供需平衡。

(3)

式(3)中：PGa為分區(qū)a中機組的有功容量之和；pk為負荷消耗的有功功率；α1和β1為約束系數(shù)；yak表示k號負荷節(jié)點與a號黑啟動機組的從屬關(guān)系。

(3)無功平衡約束，各分區(qū)內(nèi)線路的充電無功總和不應(yīng)大于黑啟動機組的進相無功上限，以保證系統(tǒng)的電壓平衡，此外系統(tǒng)中多余的無功需要進行補償或消耗。

(4)

式中：QGa為分區(qū)a中機組的無功容量之和；qk為負荷消耗的有功功率。

1.3 輸電線路啟動時間

結(jié)合文獻[11]所引入的時間參數(shù)，可以為電力系統(tǒng)中每一條線路設(shè)置啟動時間，以此作為線路的權(quán)重，進行1.1的目標函數(shù)計算。

線路操作人員基于經(jīng)驗，將線路的恢復時間劃分為樂觀估計時間tA，悲觀時間tB和最可能的估計時間tM，故線路恢復時間的均值和方差為：

tr=(tA+4tM+tB)/6

(5)

σr=(tB-tA)/6

(6)

假設(shè)黑啟動電源n到待啟動機組m之間由j條線路組成，則恢復時間的均值和方差為：

(7)

(8)

將啟動時間作為輸電線路的加權(quán)函數(shù)應(yīng)用于目標函數(shù)的計算中。

2 分區(qū)算法的設(shè)計

2.1 黑啟動分區(qū)個數(shù)

文獻[12]已經(jīng)提出了黑啟動分區(qū)個數(shù)的考慮原則，即在電力系統(tǒng)解列后產(chǎn)生的孤網(wǎng)中，當黑啟動機組個數(shù)n小于待啟動機組個數(shù)m時，分區(qū)數(shù)s等于黑啟動機組個數(shù)，反之等于待啟動機組個數(shù)，在此不進行證明。

s=min(n，m)

(9)

2.2 最短啟動時間算法

根據(jù)1.1的目標函數(shù)可知，孤網(wǎng)分區(qū)算法求解需要每一個機組、負荷節(jié)點相對于任意一個黑啟動機組的最短啟動時間，因此文中擬采用最短啟動時間算法進行計算，詳細步驟如下：

(1) 初始化：起始點設(shè)置為ds=0，集合ps為空，并標記起始點s，記k=s，連通圖中其他所有節(jié)點設(shè)為未標記點。

(2) 檢驗從所有已標記的點k到其直接連接的未標記的點j的距離，并設(shè)置：

dj=min[dj，dk+lkj]

(10)

其中，lkj為從點k到j(luò)的直接連接距離。

(3) 選取下一個點：從所有未標記的節(jié)點中，選取dj中最小的一個j，并將該點進行標記。

(4) 檢查是否所有的點都已經(jīng)標記，如果都標記，則退出算法，dj為所求的最短啟動時間，否則重復步驟2，3。

2.3 節(jié)點分區(qū)算法

本文1.1節(jié)的目標函數(shù)可以用線性規(guī)劃的方法進行求解，再考慮到充當變量的分區(qū)機組矩陣Sg取值非0即1，故采用0-1線性規(guī)劃法進行求解。0-1型整數(shù)線性規(guī)劃是一類特殊的整數(shù)規(guī)劃，它的變量取值僅為0或1，其模型如下：

min f=cTx

(11)

(12)

此時的決策變量被稱為為0-1變量，或布爾變量。在實際問題中，如果引進0-1變量，就可以把各種需要分別討論的線性、非線性規(guī)劃問題統(tǒng)一在一個問題中進行討論了。

2.4 分區(qū)算法的實現(xiàn)步驟

孤網(wǎng)黑啟動分區(qū)算法要求嚴格按照目標函數(shù)和約束條件對系統(tǒng)中的黑啟動機組、待啟動機組和負荷三類節(jié)點進行合理分區(qū)。本節(jié)重點旨在綜合前文的思路，建立一套完整的數(shù)學分區(qū)模型，求解各個分區(qū)的最優(yōu)恢復網(wǎng)架，使得分區(qū)內(nèi)的待啟動機組、負荷節(jié)點從黑啟動機組中獲取啟動功率的時間最短，并考慮黑啟動機組約束，使各分區(qū)的容量達到平衡。

假設(shè)孤網(wǎng)中含有n臺黑啟動機組，m臺待啟動機組，l個負荷節(jié)點。為搭建分區(qū)模型，首先構(gòu)造n行m列的分區(qū)機組矩陣S和n行l(wèi)列的分區(qū)負荷矩陣Y，矩陣S由m個n維向量Sg組成，其中Sg=[s1g，s2g，…，sng]T；矩陣Y由l個n維向量組成其中Yk=[y1k，y2k，…，ylk]T。其中g(shù)和k分別為任意待恢復機組編號和任意負荷節(jié)點編號(1≤g≤m;1≤k≤l)；向量中的元素分別表示待啟動機組、負荷節(jié)點是否屬于黑啟動機組的布爾變量，即孤網(wǎng)中所有黑啟動機組和g號待啟動機組、k號負荷節(jié)點的從屬關(guān)系?？梢娋仃嘢和Y能充分反映孤網(wǎng)內(nèi)所有待啟動機組、負荷節(jié)點和黑啟動機組之間的從屬關(guān)系。

再構(gòu)造一個n行m列的機組恢復時間矩陣B和n行l(wèi)列的負荷恢復時間矩陣C，Bg=[b1g，b2g，…，bng]T；Ck=[c1k，c2k，…，cnk]T。其中任意元素表示待啟動機組g和負荷節(jié)點k到各個黑啟動機組所需的最短啟動時間，該矩陣中數(shù)值的計算方法可以運用2.2中描述的最短啟動時間算法。

綜上所述，待啟動機組最優(yōu)分組問題可描述為：

(13)

首先求解式(13)，其中加入了黑啟動機組約束條件，即1.2中所提及的約束條件1。這樣可避免待啟動機組的分配不均勻問題，同時初步保證各分區(qū)機組容量的平衡。式(9)采用0-1線性規(guī)劃法求解，解得黑啟動機組和待啟動機組的從屬情況。

(14)

其次求解式(14)，通過負荷節(jié)點的合理分配滿足孤網(wǎng)分區(qū)的有功功率平衡和無功功率相對平衡，即1.2中所提及的約束條件2和3。文中先對待啟動機組分區(qū)后再計算每個分區(qū)的有功和無功容量，按需求對各個分區(qū)劃入負荷節(jié)點以平衡相關(guān)功率, 得到所有負荷的分區(qū)結(jié)果。該方程同樣可使用0-1線性規(guī)劃求解。詳細流程如圖1所示。

圖1 分區(qū)算法

3 算例分析

在Matlab平臺上編寫程序，算例采用ieee39節(jié)點系統(tǒng)，系統(tǒng)結(jié)構(gòu)如圖2所示。假設(shè)節(jié)點30，31和34為黑啟動機組節(jié)點，32，33，35，36，37，38，39為待啟動機組，剩下的節(jié)點即為負荷節(jié)點。

圖2 新英格蘭10機39節(jié)點系統(tǒng)分區(qū)結(jié)果

3.1 分區(qū)個數(shù)

比較黑啟動機組和待啟動機組的個數(shù)，取兩者中的最小值3為黑啟動分區(qū)個數(shù)。

3.2 機組分區(qū)

以節(jié)點30，31，34分別為分區(qū)1、分區(qū)2、分區(qū)3的源節(jié)點。針對7臺待啟動機組構(gòu)造7個恢復時間矩陣。分別以節(jié)點32，33，35，36，37，38，39為末節(jié)點進行最短啟動時間計算，求解數(shù)值放入7個3維恢復時間矩陣B1，B2，…，B7中。對于系統(tǒng)中任意一個待啟動節(jié)點，如節(jié)點37，從3個分區(qū)的黑啟動機組獲取啟動功率的時間分別為6 min，20 min，19 min，這樣節(jié)點37所對應(yīng)的恢復時間矩陣為[6，20，19]T，待啟動機組的恢復時間向量如表1所示。

表1 新英格蘭10機39節(jié)點系統(tǒng)中待恢復機組的恢復矩陣

將求解得到的恢復時間矩陣代入公式12，可以求得分區(qū)機組矩陣為：

(15)

由求解得到的分區(qū)機組矩陣可得，節(jié)點37，38，39與黑啟動節(jié)點30同屬于分區(qū)1；節(jié)點33，35，36與黑啟動節(jié)點34同屬于分區(qū)2；節(jié)點32與黑啟動節(jié)點31同屬于分區(qū)3。

3.3 負荷分區(qū)

在得到待啟動機組分區(qū)后，可以繼續(xù)進行負荷分區(qū)的計算。負荷分區(qū)計算與待啟動機組分區(qū)相類似，取系數(shù)α1，β1為0.8，α2，β2為1.2，分別以3個黑啟動節(jié)點為源點，29個負荷節(jié)點為末端進行最短啟動時間計算，構(gòu)造29個3維恢復時間矩陣，代入式(14)計算，得到的分區(qū)結(jié)果如表2所示。

表2 新英格蘭10機39節(jié)點負荷節(jié)點分區(qū)結(jié)果

對所求節(jié)點進行整理劃分，可知負荷節(jié)點1，2，3，25，26，27，28，29屬于黑啟動機組30所在的分區(qū)1；負荷節(jié)點4，5，6，7，8，9，10，11，12，13，14屬于黑啟動機組31所在的分區(qū)2；負荷節(jié)點15，16，17，18，19，20，21，22，23，24。得到如圖2所示的分區(qū)結(jié)果。

為證明本算法的有效性，將文中方法與文獻[13]所求得的分區(qū)恢復網(wǎng)架進行比較分析。文獻[13]將39號待啟動機組劃入黑啟動機組31所在的分區(qū)，啟動所需時間為15 min，而文中將39號機組劃劃入黑啟動機組30所在的分區(qū)，啟動所需時間為10 min。除此以外，文獻[13]所求得的分區(qū)1中負荷消耗有功功率超過機組發(fā)出有功17%，分區(qū)2機組發(fā)出的無功功率較負荷高出13%；文中采用了相關(guān)的約束條件進行計算，功率不平衡量可以保證在5%以內(nèi)。因此采用文中的分區(qū)方法和恢復網(wǎng)架更有利于孤網(wǎng)黑啟動。

4 結(jié)束語

將最短啟動時間算法和整數(shù)線性規(guī)劃法相結(jié)合應(yīng)用于孤網(wǎng)黑啟動分區(qū)優(yōu)化，將孤網(wǎng)中機組的恢復時間和分區(qū)的負荷容量統(tǒng)一考慮，提出一套完整的孤網(wǎng)全黑狀態(tài)下的恢復網(wǎng)架和分區(qū)優(yōu)化算法。該算法以最短時間恢復孤網(wǎng)中的重要輸電線路、發(fā)電能力、恢復網(wǎng)架為目標，求解出各機組、負荷以及輸電線路所屬的分區(qū)，構(gòu)建恢復網(wǎng)架。同時以分區(qū)內(nèi)的機組穩(wěn)態(tài)有功功率、無功功率平衡為約束，充分保證孤網(wǎng)分區(qū)的合理性和系統(tǒng)黑啟動的穩(wěn)定性。

[1] 范鈞慧，徐 楠，劉皓明. 含分布式風光電源的配電系統(tǒng)故障恢復策略[J]. 江蘇電機工程，2014，33(1)：1-4.

[2] MOREIRA C L， RESENDEFO， LOPES J A P. Using Low Voltage Microgrids for Service Restoration[J]. IEEE Transactions on Power Systems，2007，22(1)：395-403.

[3] 周云海，閔 勇. 恢復控制中的系統(tǒng)重構(gòu)優(yōu)化算法研究[J]. 中國電機工程學報，2003，23(4)：67-70.[4] 劉映尚，吳文傳，馮永青. 基于有序二元決策圖的黑啟動分區(qū)搜索策略[J]. 中國電機工程學報，2008，28(10)：26-31.

[5] 梁海平，郝 杰，顧雪平. 計及節(jié)點恢復成功率的黑啟動分區(qū)恢復方案優(yōu)化[J]. 電工技術(shù)學報，2012，27(11)：230-238.

[6] 梁海平，蘆佳碩，顧雪平. 基于恢復路徑轉(zhuǎn)移系數(shù)的電網(wǎng)黑啟動分區(qū)策略[J]. 電力系統(tǒng)及其自動化學報，2014，26(2)：1-7.

[7] 張 磊，張 璨，林振智，等. 大停電后電力系統(tǒng)黑啟動分區(qū)的兩步策略[J]. 電力自動化設(shè)備，2015，35(9)：14-21.

[8] 梁海平，顧雪平. 基于節(jié)點電壓相近度的黑啟動子系統(tǒng)劃分方法研究[J]. 電力系統(tǒng)保護與控制，2013，41(14)：81-86.

[9] 顧雪平，韓忠暉，梁海平. 電力系統(tǒng)大停電后系統(tǒng)恢復的優(yōu)化分區(qū)算法[J]. 中國電機工程學報，2009，29(10)：41-46.

[10] 門 冉，高小濤，盛昌棟. 含分布式風光電源的配電系統(tǒng)故障恢復策略[J]. 江蘇電機工程，2015，34(1)：9-12.

[11] MOTA A A，MOTA L T M，MORELATO A． Dynamic Evaluation of Reenergization Times During Power Systems Restoration[C]∥IEEE/PES Transmission & Distribution Conference & Exposition，Latin American，Sao Paulo，2004． [12] SARMADI S A N， DOBAKHSHARI A S， AZIZI S， et al. A Sectionalizing Method in Power System Restoration Based on WAMS[J]. IEEE Trans on Smart Grid，2011，02(1)：190-197.

[13] 林振智，文福拴，周 浩. 基于復雜網(wǎng)絡(luò)社團結(jié)構(gòu)的恢復子系統(tǒng)規(guī)劃算法[J]. 電力系統(tǒng)自動化，2009，33(12)：12-16.

陳祺煒

陳祺煒(1994 —)，男，江蘇南京人，碩士研究生在讀，研究方向為電力系統(tǒng)安全防御和系統(tǒng)恢復；

吳 熙(1987 —)，男，江蘇靖江人，副研究員，研究方向為電力系統(tǒng)低頻振蕩、次同步振蕩、新能源發(fā)電技術(shù)、FACTS建模與分析；

李海峰(1973 —)，男，河北懷安人，研究員級高級工程師，從事電網(wǎng)調(diào)度運行管理工作；

羅凱明(1978 —)，男，貴州惠水人，高級工程師，從事電網(wǎng)調(diào)度運行管理工作；

劉 林(1985 —)，男，江蘇泰州人，博士，從事電網(wǎng)調(diào)度運行管理工作。

Recovery Network and Partition Optimization Algorithm of Isolated Networks Under the Full Black State

CHEN Qiwei1,WU Xi1,LI Haifeng2,LUO Kaiming2,LIU Lin2

(1. School of Electrical Engineering, Southeast University, Nanjing 210096, China;2. State Grid Jiangsu Electric Power Company Dispatch and Control Center, Nanjing 210024, China)

The isolated grid which is splitting after the grid fault has become an important part of the power system. Giving priority to the isolated grid restoration has the advantage of simplifying calculation, accelerating the recovery and restoring the important load in time. It needs to be divided into several partitions in the total black state. This paper is aimed at the rapidity and stability of the restoration, comprehensively considering isolated network partitioning strategy and the restoration skeleton, introducing the shortest path algorithm and integer linear programming, proposing a new optimization algorithm for isolated network partitioning. The algorithm makes the shortest start time of unit as the objective function, constructed by the active power and reactive power balance of the unit to simplify a complex multi objective nonlinear programming problem as a single objective and multi constrained linear programming problem. New England 10-Machine 39-Bus System is used as an example to verify the effectiveness of the proposed algorithm.

black-start; isolated network partition; parallel recovery; recovery network

2016-08-26；

2016-10-12基金項目:國家電網(wǎng)公司科技項目(電網(wǎng)黑啟動過程中利用直流輸電的電網(wǎng)恢復技術(shù)研究)

TM711

A

2096-3203(2017)01-0074-05