鄧慧瓊，羅杰,3，王曉銘，劉昊，李培強(qiáng)，李晨晨

(1.福建工程學(xué)院電子電氣與物理學(xué)院，福州市 350108；2.智能電網(wǎng)仿真分析與綜合控制福建省高校工程研究中心(福建工程學(xué)院)，福州市 350108；3.國網(wǎng)永安市供電公司，福建省永安市 366000)

0 引 言

現(xiàn)代電網(wǎng)已向著高度智能化、復(fù)雜化、大容量的互聯(lián)網(wǎng)趨勢(shì)發(fā)展，促進(jìn)社會(huì)生活水平和經(jīng)濟(jì)快速發(fā)展的同時(shí)，也伴隨著一系列故障問題。某個(gè)小范圍內(nèi)的故障就可能迅速涉及整個(gè)電網(wǎng)，進(jìn)而發(fā)展為大停電事故。因此，吸取事故的經(jīng)驗(yàn)教訓(xùn)，研究連鎖故障發(fā)展過程并制定完善的連鎖故障防御決策，對(duì)提升電網(wǎng)運(yùn)行的安全性和可靠性具有重要應(yīng)用價(jià)值。

現(xiàn)階段針對(duì)連鎖故障建模方面的研究，最引人注目的是基于自組織邊界理論的有關(guān)模型[1-6]，如基于電力系統(tǒng)薄弱環(huán)節(jié)的OPA模型[1-2]、基于元件級(jí)聯(lián)失效的CASCADE模型[3-4]、基于交直流潮流的Manchester模型[5-6]等。這些模型對(duì)連鎖故障引發(fā)大停電事故的內(nèi)在機(jī)理進(jìn)行探討，揭示大停電事故具有自組織邊界特性，對(duì)預(yù)防連鎖故障具有極大參考價(jià)值。

眾多研究表明在連鎖故障的早期階段, 繼電保護(hù)的動(dòng)作以及潮流重新分配是使電網(wǎng)運(yùn)行狀態(tài)進(jìn)一步惡化的主因。因此，及時(shí)采取預(yù)防控制措施顯得尤為重要。現(xiàn)階段針對(duì)連鎖故障控制方面的研究，文獻(xiàn)[7]提出一種基于多智能體系統(tǒng)理論的廣域協(xié)同預(yù)控制系統(tǒng)，防止區(qū)域內(nèi)故障引發(fā)連鎖故障。文獻(xiàn)[8]提出一種基于不完全信息多階段對(duì)策的連鎖故障預(yù)防控制模型,為及時(shí)阻斷故障提供了依據(jù)。這些基于確定性的控制策略從諸多角度考慮了安全性，但對(duì)經(jīng)濟(jì)性考慮不足。文獻(xiàn)[9-10]引入多階段動(dòng)態(tài)博弈論分析復(fù)雜系統(tǒng)連鎖性故障的過程，給出連鎖故障事故鏈的風(fēng)險(xiǎn)評(píng)估方法，從而確定最佳的預(yù)防控制方案。這些基于不確定性的控制策略不便于確定具體的故障位置和事故的發(fā)展態(tài)勢(shì)，通常只能用于安全評(píng)估。文獻(xiàn)[11-12]建立了考慮安全性和經(jīng)濟(jì)性的連鎖故障協(xié)調(diào)控制模型，通過求解得到發(fā)電出力調(diào)整和切負(fù)荷策略。這些文獻(xiàn)雖給出了兼顧安全性和經(jīng)濟(jì)性的預(yù)防控制策略，但都是基于不確定性來評(píng)價(jià)安全，且未對(duì)初始故障進(jìn)行篩選，這樣做可能會(huì)引起計(jì)算量過大的問題。同時(shí)，以上的各種控制策略中，對(duì)于連鎖故障的模擬，雖然考慮了繼電保護(hù)的動(dòng)作，但通常都未與保護(hù)的動(dòng)作行為相結(jié)合，這不利于反映真實(shí)的連鎖故障觸發(fā)。

針對(duì)現(xiàn)有研究的不足，本文針對(duì)連鎖故障的觸發(fā)階段，將故障中的繼電保護(hù)動(dòng)作特性加以考慮，從而使連鎖故障的描述更貼近實(shí)際。在此基礎(chǔ)上，以電網(wǎng)實(shí)際運(yùn)行狀態(tài)與可觸發(fā)連鎖故障邊界狀態(tài)之間的最短距離來衡量安全性，以發(fā)電運(yùn)行成本來衡量經(jīng)濟(jì)性，通過支路脆弱性評(píng)估方法選取初始故障，進(jìn)而構(gòu)建針對(duì)不同初始故障場(chǎng)景下的連鎖故障預(yù)防控制模型。

1 連鎖故障的數(shù)學(xué)表現(xiàn)形式

本文針對(duì)一種常見的連鎖故障模式展開研究：初始故障切斷后，伴隨著電網(wǎng)發(fā)生潮流轉(zhuǎn)移，剩余支路因相應(yīng)電氣量進(jìn)入后備保護(hù)的動(dòng)作區(qū)而導(dǎo)致連鎖故障的發(fā)生。

本文以各支路都配置距離型后備保護(hù)為例，假設(shè)在某一電網(wǎng)中，某時(shí)刻其支路La發(fā)生初始故障被切斷后，剩余任一支路Lb是否會(huì)發(fā)生連鎖故障，可根據(jù)其配置的距離保護(hù)檢測(cè)到的電氣量是否進(jìn)入保護(hù)動(dòng)作區(qū)來判別。設(shè)支路Lb介于節(jié)點(diǎn)i和節(jié)點(diǎn)j之間，其i側(cè)配置有距離型后備保護(hù)，根據(jù)保護(hù)整定值和測(cè)量阻抗，可定義如式(1)所示的方程：

(1)

設(shè)支路Lb在電網(wǎng)中的編號(hào)為b，除初始故障外的剩余支路共有n條，將所有剩余支路都考慮進(jìn)來，進(jìn)一步可給出表達(dá)式如式(2)所示:

(2)

由式(2)可知，從電網(wǎng)層面上來看，當(dāng)Z<0時(shí)，此時(shí)電網(wǎng)處于連鎖故障狀態(tài)；當(dāng)Z>0時(shí)，此時(shí)電網(wǎng)處于安全狀態(tài)；而當(dāng)Z=0時(shí)，此時(shí)電網(wǎng)剛好處于連鎖故障的邊界狀態(tài)。顯然，式(2)給出了判別電網(wǎng)連鎖故障的數(shù)學(xué)表現(xiàn)形式。

2 初始故障的選取

連鎖故障往往由一個(gè)簡(jiǎn)單的初始故障開始，繼而引發(fā)后續(xù)故障，因此，選取出系統(tǒng)中較關(guān)鍵的初始故障并進(jìn)行預(yù)防控制顯得尤為重要。本節(jié)基于電網(wǎng)拓?fù)浣Y(jié)構(gòu)和運(yùn)行狀態(tài)，給出脆弱指標(biāo)的定義和評(píng)估方法，以此確定初始故障。

1)嚴(yán)重度指標(biāo)。

根據(jù)前文分析，在不考慮后備保護(hù)動(dòng)作的不確定性和其他閉鎖條件的情況下，可利用剩余支路受故障支路切斷后的擾動(dòng)影響，評(píng)估故障支路的嚴(yán)重度。設(shè)故障支路La在電網(wǎng)中的編號(hào)為a，將式(1)進(jìn)一步處理后可得到嚴(yán)重度指標(biāo)為：

(3)

式中：wi、wj為權(quán)重，本文認(rèn)為支路i側(cè)和j側(cè)對(duì)嚴(yán)重度指標(biāo)K(a)同等重要，所以權(quán)重均取0.5。K(a)綜合了兩側(cè)的距離保護(hù)，將剩余支路兩側(cè)的最大擾動(dòng)程度定義為故障支路的嚴(yán)重度指標(biāo)。

2)加權(quán)潮流熵指標(biāo)。

本文從支路抗沖擊能力和系統(tǒng)潮流轉(zhuǎn)移的層面出發(fā)，利用支路負(fù)載率對(duì)潮流熵指標(biāo)進(jìn)行加權(quán)，加權(quán)后的潮流熵指標(biāo)考慮了支路負(fù)載率和支路總數(shù)對(duì)系統(tǒng)運(yùn)行狀態(tài)的影響，克服了傳統(tǒng)模型只計(jì)及轉(zhuǎn)移潮流分布均衡性的問題。定義故障支路La切斷后的加權(quán)潮流熵為[13]：

(4)

式中：λab為故障支路La切斷后對(duì)支路Lb的潮流沖擊率；γab為故障支路La切斷后對(duì)支路Lb的負(fù)載率。為了體現(xiàn)過載危險(xiǎn)性，引入過載系數(shù)μ，當(dāng)γab<1時(shí)，μ取為1；當(dāng)γab≥1時(shí)，μ取為10。支路Lb受電網(wǎng)其他支路潮流沖擊影響越大且自身抵抗潮流沖擊的能力越弱，支路故障的可能性就越大。

3)綜合脆弱指標(biāo)。

為了更全面、合理地評(píng)價(jià)支路的脆弱性，本文進(jìn)一步考慮網(wǎng)絡(luò)效能指標(biāo)S(a)[14]、支路介數(shù)B(a)[15]和電壓偏移指標(biāo)U(a)[16]，將各分項(xiàng)指標(biāo)歸一化處理后，加權(quán)得到綜合脆弱指標(biāo)F(a)為：

F(a)=w1K′(a)+w2V′(a)+w3S′(a)+w4B′(a)+w5U′(a)

(5)

4)脆弱性評(píng)估方法。

基于前述分析，將各分項(xiàng)指標(biāo)歸一值和綜合脆弱指標(biāo)一并作為輸入量，通過聚類算法將支路按其指標(biāo)的大小及組合進(jìn)行劃分，從而找出脆弱性較強(qiáng)的支路，并將其視為初始故障。其中，本文擬采用聚類算法中的模糊C均值算法(fuzzy C-means algorithm, FCM)，F(xiàn)CM是通過目標(biāo)函數(shù)優(yōu)化迭代找到一組中心矢量，使各樣本到其距離的平方和最小，以此確定最佳分類數(shù)和聚類中心[18]。

3 連鎖故障預(yù)防控制模型

針對(duì)前述選取出的初始故障，本節(jié)將介紹如何建立考慮安全性和經(jīng)濟(jì)性的預(yù)防控制模型。

3.1 系統(tǒng)安全性

一般來說，在連鎖故障的觸發(fā)階段，初始故障切斷前后的電網(wǎng)節(jié)點(diǎn)注入功率基本不變。由此可知，對(duì)于給定的電網(wǎng)，決定電網(wǎng)是否發(fā)生連鎖故障的主要因素是初始故障切斷前的節(jié)點(diǎn)注入功率。進(jìn)一步地，依據(jù)節(jié)點(diǎn)注入功率可給出一種分析電網(wǎng)對(duì)于連鎖故障安全裕度指標(biāo)[19]。

為敘述方便，以一個(gè)含兩臺(tái)發(fā)電機(jī)的系統(tǒng)為例，假設(shè)系統(tǒng)的無功功率已取得平衡，可以不予考慮，且負(fù)荷功率不變，此時(shí)系統(tǒng)的狀態(tài)主要決定于兩發(fā)電機(jī)節(jié)點(diǎn)上的有功功率P1和P2(分別對(duì)應(yīng)于第1臺(tái)發(fā)電機(jī)節(jié)點(diǎn)和第2臺(tái)發(fā)電機(jī)節(jié)點(diǎn))?，F(xiàn)以P1為橫坐標(biāo)，以P2為縱坐標(biāo)，將此系統(tǒng)任一狀態(tài)下(P1，P2)對(duì)應(yīng)的點(diǎn)畫在坐標(biāo)圖中，便可得到如圖1所示的情形。

圖1 電網(wǎng)對(duì)于連鎖故障的安全裕度Fig.1 Security margin of the power grid for cascading failure

(6)

若推廣到多于兩臺(tái)發(fā)電機(jī)的系統(tǒng)，并同時(shí)考慮有功功率和無功功率，則需將式(6)中的有功功率向量換成復(fù)數(shù)形式，用電網(wǎng)的任一節(jié)點(diǎn)注入功率向量S來表示，可寫成式(7)的形式：

(7)

進(jìn)一步地，面向選取出的多個(gè)初始故障，對(duì)應(yīng)的安全裕度也應(yīng)是多個(gè)，故將式(7)處理后可得到針對(duì)多個(gè)初始故障作用下的安全性指標(biāo)：

R′=min(Re)

(8)

式中：e為初始故障的編號(hào)；Re為第e個(gè)初始故障對(duì)應(yīng)的安全裕度；R′為各初始故障對(duì)應(yīng)的安全裕度最小值。經(jīng)過處理后，所得到的安全性指標(biāo)可表征在多個(gè)初始故障作用下的電網(wǎng)安全裕度。若電網(wǎng)的運(yùn)行狀態(tài)處于以s為圓心，以R′為半徑的安全范圍內(nèi)，則電網(wǎng)是安全的。

3.2 系統(tǒng)經(jīng)濟(jì)性

本文主要研究如何通過預(yù)防控制方法提高電網(wǎng)對(duì)于連鎖故障的安全裕度，避免電網(wǎng)觸發(fā)連鎖故障，對(duì)于預(yù)防控制不能奏效而需要采取切負(fù)荷等緊急控制措施的情形，限于篇幅，本文暫不考慮。因此，定義發(fā)電運(yùn)行成本作為系統(tǒng)的經(jīng)濟(jì)性指標(biāo)，發(fā)電運(yùn)行成本C可表示為：

(9)

式中：G為參與調(diào)整的發(fā)電機(jī)集合；Pg為發(fā)電機(jī)g的有功出力；ag、bg、cg為其運(yùn)行成本系數(shù)。

3.3 數(shù)學(xué)模型

本文以發(fā)電機(jī)的節(jié)點(diǎn)注入功率作為控制變量，建立了考慮安全性和經(jīng)濟(jì)性的連鎖故障預(yù)防控制模型，該模型可分為兩個(gè)層次(分別稱為內(nèi)層和外層)。其中，外層模型通過調(diào)整當(dāng)前運(yùn)行點(diǎn)至新的位置，使其對(duì)應(yīng)的電網(wǎng)安全裕度盡可能大，而內(nèi)層模型根據(jù)外層模型傳遞進(jìn)來的運(yùn)行點(diǎn)數(shù)據(jù)，計(jì)算出該運(yùn)行點(diǎn)對(duì)應(yīng)的電網(wǎng)安全裕度，外層模型同時(shí)還適當(dāng)兼顧了經(jīng)濟(jì)性的考量。

3.3.1 外層優(yōu)化模型

在外層模型中，先確定一個(gè)電網(wǎng)當(dāng)前運(yùn)行點(diǎn)，通過調(diào)整當(dāng)前運(yùn)行點(diǎn)至新的位置，使其對(duì)應(yīng)的電網(wǎng)安全裕度盡可能大，發(fā)電運(yùn)行成本盡可能小，其目標(biāo)函數(shù)F1和F2可表示為：

(10)

式中：C為發(fā)電運(yùn)行成本。

約束條件為：

1)電網(wǎng)的安全約束。

由式(2)可知，在S的作用下，當(dāng)初始故障切斷后，電網(wǎng)處于不發(fā)生連鎖故障的安全狀態(tài)條件：

(11)

2)潮流約束。

電網(wǎng)在安全狀態(tài)下，還需要保證初始故障前后滿足潮流約束條件。其中，等式約束為系統(tǒng)節(jié)點(diǎn)潮流約束；而不等式約束包括：發(fā)電機(jī)有功出力和無功出力約束、節(jié)點(diǎn)電壓約束、支路傳輸功率約束。為方便敘述，將其縮寫成式(12)的形式：

(12)

3.3.2 內(nèi)層優(yōu)化模型

在內(nèi)層模型中，根據(jù)外層模型傳遞進(jìn)來的當(dāng)前運(yùn)行點(diǎn)數(shù)據(jù)，計(jì)算出該運(yùn)行點(diǎn)對(duì)應(yīng)的電網(wǎng)安全裕度R′，其目標(biāo)函數(shù)F3可表示為：

F3=R′

(13)

約束條件為：

1)電網(wǎng)的邊界約束。

內(nèi)層模型是以電網(wǎng)實(shí)際運(yùn)行狀態(tài)與可觸發(fā)的連鎖故障邊界狀態(tài)之間的最短距離為目標(biāo)，所以在S的作用下，當(dāng)初始故障切斷后，電網(wǎng)需滿足剛好處于連鎖故障邊界狀態(tài)的條件：

(14)

2)潮流約束。

電網(wǎng)在連鎖故障邊界狀態(tài)下，同樣需要滿足一定的潮流約束，給出類似式(12)的約束條件為：

(15)

4 模型的求解思路

本文利用內(nèi)外層嵌套的優(yōu)化算法求解預(yù)防控制模型，內(nèi)層模型采用粒子群算法(particle swarm optimization, PSO)[20]，外層模型采用基于自適應(yīng)網(wǎng)格的多目標(biāo)粒子群(multi-objective PSO based on adaptive grid，AG-MOPSO)算法[21]。在求解過程中，內(nèi)外層之間的關(guān)鍵變量及結(jié)果將根據(jù)需要進(jìn)行交互傳遞，交互體現(xiàn)為：外層優(yōu)化每確定一個(gè)當(dāng)前運(yùn)行點(diǎn)，都將該運(yùn)行點(diǎn)數(shù)據(jù)傳遞至內(nèi)層優(yōu)化中，由內(nèi)層優(yōu)化計(jì)算出該運(yùn)行點(diǎn)對(duì)應(yīng)的電網(wǎng)安全裕度，并將其傳遞到外層優(yōu)化中。這一交互操作具有時(shí)序關(guān)系，在確定每一個(gè)當(dāng)前運(yùn)行點(diǎn)時(shí)，先從外層操作開始，然后進(jìn)入內(nèi)層操作，最后再返回到外層操作，如此循環(huán)往復(fù)地進(jìn)行交替操作，但在每一步的操作時(shí)間上，并沒有嚴(yán)格的限制。

4.1 粒子更新準(zhǔn)則

PSO和AG-MOPSO算法中，每個(gè)粒子的位置和速率都是按照式(16)進(jìn)行更新[20]：

(16)

式中：yk(t+1)和vk(t+1)分別為粒子k在第t+1次的位置和速率；yp.k(t)為粒子的個(gè)體最佳位置；ybset為粒子群體的全局最佳位置；w為慣性權(quán)重，其取值從0.9至0.1按線性規(guī)律遞減；β1和β2為加速系數(shù)，其取值均為2；r1和r2均為0至1區(qū)間內(nèi)的隨機(jī)數(shù)。

在使用算法時(shí)，每一個(gè)粒子對(duì)應(yīng)于一個(gè)節(jié)點(diǎn)注入功率組合，而初始粒子群以電網(wǎng)當(dāng)前運(yùn)行狀態(tài)下的節(jié)點(diǎn)注入功率S為基礎(chǔ)，通過在S中的各個(gè)元素上隨機(jī)增加一個(gè)微小值來獲得。在迭代過程中，對(duì)于初始故障前的潮流方程，擬采用電網(wǎng)當(dāng)前運(yùn)行狀態(tài)下的潮流解去求解，以加快計(jì)算效率；而對(duì)于初始故障后的潮流方程，本文擬采用牛頓-拉夫遜法求解，以取得較精確的計(jì)算結(jié)果。

4.2 適應(yīng)度函數(shù)

本文引入懲罰函數(shù)對(duì)模型中的約束問題做進(jìn)一步處理[22], 該方法的具體過程為：在外層模型中，將式(11)中的Z>0寫為p(x)>0的形式，將式(12)中的所有潮流等式約束和不等式約束分別寫成d(x)=0和h(x)≥0的形式；在內(nèi)層模型中，將式(14)中的Z=0寫成z(x)=0，將式(15)中的所有潮流等式約束和不等式約束分別寫成m(x)=0和g(x)≥0的形式，加入懲罰因子，利用目標(biāo)函數(shù)和約束函數(shù)共同構(gòu)造如下所示的適應(yīng)度函數(shù)：

(17)

式中：M1、M2、M3均為懲罰因子；f1、f2為外層模型的適應(yīng)度函數(shù)；F1、F2為外層模型的目標(biāo)函數(shù)；f3為內(nèi)層模型的適應(yīng)度函數(shù)；F3為內(nèi)層模型的目標(biāo)函數(shù)。

4.3 折衷最優(yōu)解的選取

(18)

各Pareto解的標(biāo)準(zhǔn)化隸屬度計(jì)算公式為：

(19)

式中:l為Pareto解的個(gè)數(shù)；λr為第r個(gè)Pareto解的綜合隸屬度，選擇λr最大的解作為折衷最優(yōu)解。

4.4 算法流程

本文基于PSO算法和AG-MOPSO算法對(duì)雙層優(yōu)化的預(yù)防控制模型進(jìn)行求解。對(duì)于PSO算法的優(yōu)化過程，其具體步驟詳見文獻(xiàn)[20]，限于篇幅，僅給出如圖2所示的總體算法流程。

圖2 算法流程Fig.2 Flowchart of the algorithm

5 算例分析

本文以IEEE 39節(jié)點(diǎn)系統(tǒng)和IEEE 118節(jié)點(diǎn)系統(tǒng)為算例，在Matlab軟件上進(jìn)行仿真驗(yàn)證。IEEE 39節(jié)點(diǎn)系統(tǒng)包含9臺(tái)發(fā)電機(jī)、1臺(tái)平衡機(jī)和46條支路，IEEE 118節(jié)點(diǎn)系統(tǒng)包含53臺(tái)發(fā)電機(jī)、1臺(tái)平衡機(jī)和186條支路，各元件參數(shù)和節(jié)點(diǎn)數(shù)據(jù)取自BPA軟件。為便于計(jì)算，一般采用標(biāo)幺值pu表示。PSO和AG-MOPSO算法的參數(shù)設(shè)置：算法最大迭代次數(shù)均取為200，種群數(shù)均取為40。

5.1 IEEE 39節(jié)點(diǎn)系統(tǒng)仿真結(jié)果分析

對(duì)于式(5)中的各分項(xiàng)指標(biāo)權(quán)重，根據(jù)熵權(quán)法計(jì)算得出結(jié)果為：w1=0.238 4，w2=0.351 6，w3=0.107 9，w4=0.136 8，w5=0.165 2。依次切斷非發(fā)電機(jī)出口支路，將各分項(xiàng)指標(biāo)的歸一值和綜合脆弱指標(biāo)作為輸入量進(jìn)行聚類劃分，最終初始故障被分為3類，各類聚類中心如表1所示。

表1 IEEE 39節(jié)點(diǎn)的聚類中心Table 1 Cluster center of IEEE 39 bus

由表1可知，A1類各項(xiàng)指標(biāo)的聚類中心都較大，因此，將A1類支路視為初始故障。A1類中共有4條支路，分別為：L2-3(介于節(jié)點(diǎn)2和節(jié)點(diǎn)3之間的支路，以此類推)、L15-16、L21-22、L26-27，其各項(xiàng)指標(biāo)排序結(jié)果如表2所示。

表2 IEEE 39節(jié)點(diǎn)的脆弱支路Table 2 Vulnerable branches of IEEE 39 bus system

由表2可知，各分項(xiàng)指標(biāo)單獨(dú)排序與綜合脆弱指標(biāo)F(a)的排序結(jié)果差異較大，綜合脆弱指標(biāo)排序靠前的支路往往同時(shí)包含脆弱性強(qiáng)和脆弱性弱的指標(biāo)，說明單一指標(biāo)評(píng)價(jià)支路脆弱性考慮不夠全面，可能忽略某些薄弱環(huán)節(jié)，而綜合脆弱指標(biāo)同時(shí)考慮網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和支路切斷后的運(yùn)行狀態(tài)，具有一定優(yōu)越性。

由于實(shí)際電網(wǎng)規(guī)模大，支路數(shù)量多，相應(yīng)約束條件也多，將電網(wǎng)連鎖故障預(yù)防控制分析的范圍縮小至某些脆弱的初始故障，顯著提升了計(jì)算效率，降低了進(jìn)行預(yù)防控制的難度。因此，在選取出初始故障的基礎(chǔ)上，按照預(yù)防控制模型求解出在多個(gè)初始故障作用下的Pareto最優(yōu)前沿(Pareto optimal front, POF)如圖3所示。

圖3 IEEE 39節(jié)點(diǎn)的POFFig.3 POF from IEEE 39 bus system

圖3最終得到一條由Pareto解構(gòu)成的POF曲線，通過隸屬度函數(shù)計(jì)算各Pareto解的隸屬度，根據(jù)綜合隸屬度的值找出折衷最優(yōu)解，其最優(yōu)解取值：電網(wǎng)安全裕度為8.19 pu，發(fā)電運(yùn)行成本為1 045.86 pu。其最優(yōu)解對(duì)應(yīng)的發(fā)電機(jī)有功出力結(jié)果如表3所示,表中的功率數(shù)據(jù)若為正數(shù)，則表示該功率實(shí)際是注入相應(yīng)節(jié)點(diǎn)，若為負(fù)數(shù)，則相反。

表3 IEEE 39節(jié)點(diǎn)的預(yù)防控制優(yōu)化結(jié)果Table 3 Optimization results of preventive control in IEEE 39 bus system MW

進(jìn)一步地，為了有效衡量預(yù)防控制的效果，表4給出了預(yù)防控制前后的對(duì)比情況。

表4 預(yù)防控制前后對(duì)比結(jié)果Table 4 Comparison results before and after preventive control pu

由對(duì)比結(jié)果可知，在不同初始故障的沖擊下，控制后的電網(wǎng)安全性指標(biāo)和經(jīng)濟(jì)性指標(biāo)均優(yōu)于控制前。由此說明本文提出的預(yù)防控制策略能在滿足系統(tǒng)各種約束條件的前提下，以更低的經(jīng)濟(jì)成本獲得更高的安全裕度，保證了電網(wǎng)能同時(shí)兼顧安全性和經(jīng)濟(jì)性。

5.2 IEEE 118節(jié)點(diǎn)系統(tǒng)仿真結(jié)果分析

為了進(jìn)一步驗(yàn)證本文方法的可行性，本節(jié)以系統(tǒng)規(guī)模較大的IEEE 118節(jié)點(diǎn)系統(tǒng)做算例分析。IEEE 118節(jié)點(diǎn)系統(tǒng)上的算例處理方法與前述IEEE 39節(jié)點(diǎn)系統(tǒng)相同，根據(jù)熵權(quán)法計(jì)算權(quán)重結(jié)果為：w1=0.201 7，w2=0.063 5，w3=0.324 7，w4=0.011 7，w5=0.398 4。利用FCM算法將初始故障分為2類，各類聚類中心如表5所示。

由表5可知，B1類各項(xiàng)指標(biāo)的聚類中心都較大，因此，將B1類支路視為初始故障。B1類中共有14條支路，其各項(xiàng)指標(biāo)按排序結(jié)果如表6所示。

表5 IEEE 118節(jié)點(diǎn)的聚類中心Table 5 Cluster center of IEEE 118 bus system

表6 IEEE 118節(jié)點(diǎn)的脆弱支路Table 6 Vulnerable branches of IEEE 118 bus system

進(jìn)一步地，針對(duì)上述選出的初始故障，按預(yù)防控制模型求解出的POF如圖4所示。

圖4 IEEE 118節(jié)點(diǎn)的POFFig.4 POF from IEEE 118 bus system

由圖3和圖4的仿真結(jié)果可知，本文所建的預(yù)防控制模型適用于各節(jié)點(diǎn)系統(tǒng)，算法得到的Pareto解構(gòu)成一條POF曲線，POF曲線的均勻延展性較好，每個(gè)Pareto解具有不同的發(fā)電運(yùn)行成本和電網(wǎng)安全裕度，POF曲線由左至右，所對(duì)應(yīng)Pareto解的電網(wǎng)安全裕度和發(fā)電運(yùn)行成本值都是由小到大。

IEEE 118節(jié)點(diǎn)系統(tǒng)計(jì)算出的折衷最優(yōu)解取值：電網(wǎng)安全裕度為14.11 pu，發(fā)電運(yùn)行成本為1 985.91 pu。限于篇幅，表7給出了IEEE 118節(jié)點(diǎn)系統(tǒng)中最優(yōu)解對(duì)應(yīng)的部分發(fā)電機(jī)有功出力結(jié)果。表8給出了預(yù)防控制前后的對(duì)比結(jié)果。

表7 IEEE 118節(jié)點(diǎn)的預(yù)防控制優(yōu)化結(jié)果Table 7 Optimization results of preventive control in IEEE 118 bus system MW

表8 預(yù)防控制前后對(duì)比結(jié)果Table 8 Comparison results before and after preventive control pu

由上述算例結(jié)果可知，面向所有的關(guān)鍵初始故障，經(jīng)過預(yù)防控制后，電網(wǎng)的安全性和經(jīng)濟(jì)性都得到一定的提升。此時(shí)電網(wǎng)進(jìn)入一個(gè)更安全的運(yùn)行狀態(tài)，當(dāng)電網(wǎng)節(jié)點(diǎn)注入功率在以s為圓心，以安全裕度值為半徑的安全范圍內(nèi)任意變動(dòng)時(shí)，針對(duì)不同的初始故障沖擊，電網(wǎng)對(duì)于連鎖故障的觸發(fā)都是安全的，即不會(huì)觸發(fā)連鎖故障?？傊疚姆椒杀阌谟^察電網(wǎng)的安全情況以及電網(wǎng)節(jié)點(diǎn)注入功率的合理變化范圍，為運(yùn)行人員提供了平衡安全性和經(jīng)濟(jì)性的最優(yōu)發(fā)電機(jī)出力策略。

6 結(jié) 論

本文在考慮安全性和經(jīng)濟(jì)性的基礎(chǔ)上，給出了預(yù)防電網(wǎng)連鎖故障的數(shù)學(xué)模型和具體求解思路，并進(jìn)行了算例驗(yàn)證，主要的結(jié)論如下：

1)本文所構(gòu)建的預(yù)防控制模型實(shí)質(zhì)上是動(dòng)態(tài)交互、相互影響的雙層優(yōu)化模型。通過求解雙層模型得到最優(yōu)發(fā)電機(jī)出力策略，使電網(wǎng)能同時(shí)兼顧安全性和經(jīng)濟(jì)性的要求。

2)對(duì)于安全性的考慮，主要是以電網(wǎng)當(dāng)前運(yùn)行狀態(tài)與可觸發(fā)連鎖故障的邊界狀態(tài)之間的最短距離來衡量安全性，其具體形式為兩種狀態(tài)下的電網(wǎng)節(jié)點(diǎn)注入功率向量之間的距離，這種形式可以直觀地量化和比較控制前后電網(wǎng)對(duì)于連鎖故障觸發(fā)的安全裕度。

3)對(duì)于初始故障的選取，定義綜合考慮網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和運(yùn)行狀態(tài)的脆弱性指標(biāo)，利用FCM算法進(jìn)行分類，保證了選出重要的初始故障。

總之，本文方法能準(zhǔn)確有效地選取出初始故障，并按照提高電網(wǎng)安全裕度和保證經(jīng)濟(jì)成本的思路進(jìn)行預(yù)防控制，為進(jìn)一步的研究提供借鑒。