基于網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)的臨界割集的識別①

孟祥俠， 張有東， 李祥忠， 潘廣林

(1.華北科技學(xué)院機(jī)電工程系， 北京 101601； 2.吉林石油集團(tuán)有限責(zé)任公司熱電廠， 松原 138000)

基于網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)的臨界割集的識別①

孟祥俠1， 張有東1， 李祥忠2， 潘廣林2

(1.華北科技學(xué)院機(jī)電工程系， 北京 101601； 2.吉林石油集團(tuán)有限責(zé)任公司熱電廠， 松原 138000)

準(zhǔn)確識別臨界割集能為監(jiān)測、控制系統(tǒng)在網(wǎng)絡(luò)中合理布點(diǎn)提供關(guān)鍵的依據(jù)。該文基于網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和參數(shù)，利用耦合加權(quán)的思想形成網(wǎng)絡(luò)中割集權(quán)系數(shù)，根據(jù)割集權(quán)系數(shù)的大小即可識別網(wǎng)絡(luò)中最易導(dǎo)致系統(tǒng)失穩(wěn)的“瓶頸”環(huán)節(jié)，并從穩(wěn)定的角度揭示了輸電網(wǎng)絡(luò)本身固有的分區(qū)分層的結(jié)構(gòu)特點(diǎn)對系統(tǒng)暫態(tài)穩(wěn)定的影響。研究結(jié)果表明：網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和參數(shù)對系統(tǒng)的穩(wěn)定性起支配的作用。對New England 10機(jī)系統(tǒng)進(jìn)行仿真分析，驗證了所提方法的有效性。

耦合； 加權(quán)； 割集權(quán)系數(shù)； 臨界割集； 暫態(tài)穩(wěn)定

近年來，由于電力系統(tǒng)跨區(qū)輸電、跨區(qū)聯(lián)網(wǎng)的形成，電網(wǎng)結(jié)構(gòu)愈加復(fù)雜化，動態(tài)穩(wěn)定問題日益突出[1，2]。因此，準(zhǔn)確地識別出制約系統(tǒng)暫態(tài)穩(wěn)定的“瓶頸”環(huán)節(jié)，以采取有效的監(jiān)測、控制措施來提高系統(tǒng)的穩(wěn)定性，是許多研究者關(guān)注的一個熱點(diǎn)。

從網(wǎng)絡(luò)入手研究電力系統(tǒng)的暫態(tài)穩(wěn)定性問題最早始于20世紀(jì)80年代初，Bergen和Hill在文獻(xiàn)[3]中，利用網(wǎng)絡(luò)元件的聯(lián)接強(qiáng)弱定義了割集的脆弱性指標(biāo)，此文雖然計及了網(wǎng)絡(luò)元件對暫態(tài)穩(wěn)定的影響，但對輸電網(wǎng)絡(luò)本身固有的分區(qū)分層的結(jié)構(gòu)特點(diǎn)沒有深入的認(rèn)識。文獻(xiàn)[4]進(jìn)一步從網(wǎng)絡(luò)結(jié)構(gòu)出發(fā)，揭示網(wǎng)絡(luò)結(jié)構(gòu)與系統(tǒng)穩(wěn)定性的關(guān)系。文獻(xiàn)[5]利用狀態(tài)變量的耦合強(qiáng)弱對電力系統(tǒng)進(jìn)行了動態(tài)分割，并指出網(wǎng)絡(luò)的動態(tài)分割主要決定于網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)和參數(shù)，此方法的提出有助于深入理解網(wǎng)絡(luò)結(jié)構(gòu)本身固有的分區(qū)分層的結(jié)構(gòu)特點(diǎn)，但未能從割集整體揭示網(wǎng)絡(luò)結(jié)構(gòu)脆弱性與系統(tǒng)暫態(tài)穩(wěn)定性的關(guān)系。文獻(xiàn)[6～10]利用支路勢能法，通過分析故障后支路勢能在網(wǎng)絡(luò)中的分布與變化特點(diǎn)識別網(wǎng)絡(luò)中的危險“斷面”，確定臨界割集，但此方法必須依賴全系統(tǒng)的仿真計算。

在此基礎(chǔ)上，本文進(jìn)一步利用耦合加權(quán)的思想形成網(wǎng)絡(luò)中割集權(quán)系數(shù)，根據(jù)割集權(quán)系數(shù)的大小即可識別出網(wǎng)絡(luò)中最易導(dǎo)致系統(tǒng)失穩(wěn)的“瓶頸”環(huán)節(jié)，并從穩(wěn)定的角度揭示了輸電網(wǎng)絡(luò)本身固有的分區(qū)分層的結(jié)構(gòu)特點(diǎn)對系統(tǒng)暫態(tài)穩(wěn)定的影響。

1 導(dǎo)納陣耦合的思想

1.1 耦合思想

(1)

方向為由j指向i，其中Δ是Y陣的行列式，Δji是Δ去掉j行i列的子行列式；邊權(quán)yji方向為由i指向j，其一階循環(huán)積

(2)

可用來衡量節(jié)點(diǎn)i和j的耦合程度，并對應(yīng)矩陣G的j行i列。

圖1 相互耦合的兩節(jié)點(diǎn)

一階循環(huán)積gji的耦合程度，可從行列式的角度對其進(jìn)行解釋，n階行列式Y(jié)陣按第j行展開為

Δ=yj1(-1)j+1Δj1+…+yjn(-1)j+nΔjn=

(gj1+…+gjn)Δ

(3)

導(dǎo)納陣Y及其逆陣Y-1分別為

(4)

(5)

因此可以將行矢量

i=1，2，…，n

(6)

對應(yīng)矩陣G的第i行，將G陣記為

G=Y⊙Y-1

(7)

并定義G陣為導(dǎo)納陣Y的耦合矩陣。

1.2 用導(dǎo)納陣Y形成耦合陣G的步驟

為保證劃分后的子系統(tǒng)內(nèi)的各臺發(fā)電機(jī)之間有直接的電氣聯(lián)系，或者說能保證按區(qū)域劃分的原則，首先應(yīng)將系統(tǒng)中的發(fā)電機(jī)節(jié)點(diǎn)重新編號，新編號的節(jié)點(diǎn)順序應(yīng)盡量使各臺發(fā)電機(jī)在地理位置上互相鄰接。具體步驟如下。

步驟1將系統(tǒng)中所有負(fù)荷節(jié)點(diǎn)的功率轉(zhuǎn)換為導(dǎo)納

(8)

使負(fù)荷節(jié)點(diǎn)轉(zhuǎn)變?yōu)榉亲⑷牍?jié)點(diǎn)。

Y0=AYbAT

(9)

形成節(jié)點(diǎn)導(dǎo)納矩陣Y0。

步驟3消去Y0陣的非注入節(jié)點(diǎn)，只保留發(fā)電機(jī)節(jié)點(diǎn)，所獲得的納陣記為Y。

步驟4求阻抗矩陣Z=Y-1，進(jìn)而求出

G=Y⊙Z

(10)

2 利用耦合值加權(quán)形成割集權(quán)系數(shù)

2.1 加權(quán)思想

G陣的每一個元都對應(yīng)一個數(shù)值，該數(shù)值稱為此元對應(yīng)邊的權(quán)。以下從割集的角度來分析，對于一個具有n個節(jié)點(diǎn)的系統(tǒng)，不失一般性，不妨選取這樣一個割集，此割集將系統(tǒng)分為m和n-m個節(jié)點(diǎn)組成的子系統(tǒng)a和b。

定義1子系統(tǒng)a中的任一節(jié)點(diǎn)i(i=1，2，…，m)對子系統(tǒng)b的分配權(quán)是指子系統(tǒng)b中諸節(jié)點(diǎn)和節(jié)點(diǎn)i之間的權(quán)之和，即

(11)

而子系統(tǒng)a和b所對應(yīng)割集的權(quán)系數(shù)為

(12)

圖2給出了節(jié)點(diǎn)分配權(quán)以及割集權(quán)系數(shù)的幾何解釋。其中：αib為節(jié)點(diǎn)i與子系統(tǒng)b中諸節(jié)點(diǎn)的

權(quán)之和，用來衡量節(jié)點(diǎn)i和子系統(tǒng)b的耦合程度；αab為a中諸節(jié)點(diǎn)與子系統(tǒng)b的節(jié)點(diǎn)分配權(quán)之和，可用來衡量子系統(tǒng)a和b的耦合強(qiáng)弱，αab值小則子系統(tǒng)a和b耦合弱，αab值大則a和b耦合強(qiáng)。

圖2 節(jié)點(diǎn)分配權(quán)αib及割集權(quán)系數(shù)αab

在對應(yīng)的矩陣

G=

中，αab則為矩陣G的非對角子塊Gab中所有元素的和。

2.2 利用加權(quán)形成割集權(quán)系數(shù)

由于重新編號的節(jié)點(diǎn)順序在地理位置上盡可能互相鄰接，因此可以按新編號的節(jié)點(diǎn)順序，對系統(tǒng)中各相鄰節(jié)點(diǎn)進(jìn)行優(yōu)化組合，圖3給出了重新編號的節(jié)點(diǎn)順序。節(jié)點(diǎn)1的組合為(1)，(1，2)，(1，3)，(1，4)，(1，2，3)，(1，2，4)，(1，2，5)，…，但不會出現(xiàn)類似(1，6)或(1，2，7)這種在地理位置上不相鄰的組合；節(jié)點(diǎn)2的組合為(2)，(2，3)，(2，5)，(2，6)，(2，3，6)，(2，5，6)，…，同樣也不會出現(xiàn)類似(2，7)或(2，3，8)這種組合，以此類推。

圖3 新編號的節(jié)點(diǎn)順序

(14)

修正為

(15)

圖4 割集權(quán)系數(shù)算法實現(xiàn)框圖

值得一提的是，對于節(jié)點(diǎn)相鄰但編號無法相鄰的情況，圖3中的節(jié)點(diǎn)1與3位置相鄰，但編號不相鄰。處理此種情況可以通過對耦合陣G中g(shù)ij進(jìn)行搜索，判斷其是否為零，因為G陣同Y陣一樣為稀疏矩陣，當(dāng)j>i+1且gij≠0時，說明節(jié)點(diǎn)i與j在位置上鄰接，此時可以將矩陣G中節(jié)點(diǎn)j處理成節(jié)點(diǎn)i+1。

3 利用割集權(quán)系數(shù)識別網(wǎng)絡(luò)中臨界割集

3.1 臨界割集的概念

若在系統(tǒng)中的任一位置發(fā)生故障，使得系統(tǒng)失去穩(wěn)定，系統(tǒng)失穩(wěn)可能呈現(xiàn)為兩群或多群振蕩，不管是兩群還是多群，其群間聯(lián)絡(luò)線通常都不是由一個割集組成，而是對應(yīng)多個割集，但其中只有一個割集被“撕裂”，也即此割集包含的所有支路兩端相角差均超過180°，此割集稱為“臨界割集”。圖5示意地給出了機(jī)群A與機(jī)群B之間的臨界割集。

圖5 機(jī)群A與機(jī)群B的臨界割集

3.2 利用割集權(quán)系數(shù)識別網(wǎng)絡(luò)中臨界割集

由于gij可用來衡量節(jié)點(diǎn)i與j的耦合程度，而

αab又是gij的函數(shù)，因此割集權(quán)系數(shù)αab可用來衡量子系統(tǒng)a和b之間的耦合強(qiáng)弱，利用其值的大小即可識別網(wǎng)絡(luò)中的脆弱環(huán)節(jié)。

由于從穩(wěn)定分析的觀點(diǎn)，主要對系統(tǒng)中各臺發(fā)電機(jī)的轉(zhuǎn)子角δi隨時間變化的函數(shù)感興趣，因此使用只保留發(fā)電機(jī)內(nèi)節(jié)點(diǎn)的割集權(quán)系數(shù)分析系統(tǒng)穩(wěn)定性是合理的，它雖然掩蓋了網(wǎng)絡(luò)的拓?fù)湫再|(zhì)，但從宏觀上可快速地識別出制約系統(tǒng)暫態(tài)穩(wěn)定的“瓶頸”環(huán)節(jié)。

然后在“瓶頸”環(huán)節(jié)處結(jié)合系統(tǒng)失穩(wěn)時群間聯(lián)絡(luò)線各支路相角差的變化特點(diǎn)，即可準(zhǔn)確識別出網(wǎng)絡(luò)中“臨界割集”。

4 算例與仿真分析

以New England 10機(jī)系統(tǒng)為算例，表1給出了10機(jī)系統(tǒng)部分割集權(quán)系數(shù)的大小。時域仿真采用中國電科院研制的綜合程序(PSASP)。在給定的潮流方式下，#30發(fā)電機(jī)出口發(fā)生三相瞬時性故障，切除時間為0.31 s時，系統(tǒng)呈現(xiàn)為{39}、{30,31,32,33,34,35,36,37,38}兩群失穩(wěn)，其部分支路兩端相角差變化曲線如圖6所示，其中，支路1-2、8-9兩端相角差超過180°而趨于無界，而所有其他支路兩端相角差均在有界范圍內(nèi)變化，因此，由支路1-2、8-9構(gòu)成的割集為臨界割集。為了不失一般性，表2給出了不同故障位置系統(tǒng)呈現(xiàn)出的失穩(wěn)模式，圖7為對應(yīng)的臨界割集。

表1 10機(jī)系統(tǒng)部分割集權(quán)系數(shù)

表2 10機(jī)系統(tǒng)不同故障位置的失穩(wěn)模式

由表2可以看出， 盡管大部分的故障并未發(fā)生在支路1-2、8-9或26-28、26-29構(gòu)成的臨界割集或其附近的割集上，但都導(dǎo)致系統(tǒng)在此兩個割集處“撕裂”。

之所以在網(wǎng)絡(luò)中有“系統(tǒng)在此處發(fā)生故障，而在彼處的割集被撕裂”，并導(dǎo)致系統(tǒng)失穩(wěn)的現(xiàn)象產(chǎn)生，在一定程度上是受網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)的制約。由表1割集權(quán)系數(shù)的大小可以看出，在機(jī)群{39}與{30，31，32，33，34，35，36，37，38}之間割集權(quán)系數(shù)最小，為0.022 323；機(jī)群{38}與{30,31，32，33，34，35，36，37，39}之間割集權(quán)系數(shù)次之，為0.035 485，說明{39}與其剩余機(jī)群、{38}與其剩余機(jī)群之間聯(lián)系相對較弱，此兩個部位為網(wǎng)絡(luò)上相對脆弱的環(huán)節(jié)。

而位于脆弱環(huán)節(jié)中的割集{1-2、8-9}和{26-28、26-29}在結(jié)構(gòu)上又分別處在{39}與其剩余機(jī)群、{38}與其剩余機(jī)群之間聯(lián)絡(luò)線的中部，容易成為振蕩中心所落之部位，因而在系統(tǒng)受到擾動后，最易在此兩個環(huán)節(jié)失穩(wěn)。

當(dāng)然，系統(tǒng)的暫態(tài)穩(wěn)定性不僅僅取決于網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和參數(shù)，還與故障位置、發(fā)電機(jī)的轉(zhuǎn)動慣量以及初始潮流等多種因素有關(guān)，但網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)和參數(shù)對系統(tǒng)的暫態(tài)穩(wěn)定性起支配的作用。

圖6 #39機(jī)與系統(tǒng)失穩(wěn)時部分支路兩端相角差變化

圖7 10機(jī)系統(tǒng)的臨界割集

5 結(jié)語

本文基于輸電網(wǎng)絡(luò)本身固有的分區(qū)分層的結(jié)構(gòu)特點(diǎn)，利用耦合加權(quán)的思想形成網(wǎng)絡(luò)中割集權(quán)系數(shù)，根據(jù)割集權(quán)系數(shù)的大小即可識別出網(wǎng)絡(luò)中最易導(dǎo)致系統(tǒng)失穩(wěn)的脆弱環(huán)節(jié)。此方法不依賴于數(shù)值仿真計算，方法簡單，計算量小，且可從網(wǎng)絡(luò)整體觀察問題，有助于解決大規(guī)?；ヂ?lián)電力系統(tǒng)長期制約暫態(tài)穩(wěn)定的“瓶頸”問題，為監(jiān)測、控制系統(tǒng)在網(wǎng)絡(luò)中的合理布點(diǎn)提供了一定的依據(jù)。

[1] 趙遵廉(Zhao Zunlian).中國電網(wǎng)的發(fā)展與展望(Development and prospect of power system in China)[J].中國電力(Electric Power), 2004,37(1):1-6.

[2] 薛禹盛(Xue Yusheng). 綜合防御由偶然故障演化為電力災(zāi)難——北美“8·14”大停電的警示(The way from a simple contingency to system-wide disaster-lessons from the Eastern Interconnection Blackout in 2003)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems), 2003, 27(18):1-5,37.

[3] Bergen A R, Hill D J. A structure preserving model for power system stability analysis[J].IEEE Trans on Power Apparatus and Systems,1981, 100(1):25-35.

[4] 趙金利,張群華,余貽鑫,等(Zhao Jinli, Zhang Qunhua, Yu Yixin,etal).輸電網(wǎng)網(wǎng)架結(jié)構(gòu)的譜聚類分析算法(Spectral clustering approach for structure analysis of transmission networks)[J].電力系統(tǒng)及其

自動化學(xué)報(Proceedings of the CSU-EPSA), 2009, 21(4):8-11,75.

[5] 王雨蓬, 馬昭彥(Wang Yupeng, Ma Zhaoyan).電力系統(tǒng)的動態(tài)分割(Dynamic partition of power system)[J].中國電機(jī)工程學(xué)報(Proceedings of the CSEE), 2001, 21(12): 45-49.

[6] 穆鋼,蔡國偉,胡哲,等(Mu Gang, Cai Guowei, Hu Zhe,etal).機(jī)網(wǎng)結(jié)合的暫態(tài)穩(wěn)定評價方法——關(guān)鍵割集組法(Network-machine combined transient stability assessment method——key cutsets method)[J].清華大學(xué)學(xué)報：自然科學(xué)版(Journal of Tsinghua University：Science and Technology),1997,37(7):97-101.

[7] 蔡國偉,穆鋼,程浩忠,等(Cai Guowei, Mu Gang, Cheng Haozhong,etal).基于支路暫態(tài)勢能和兩端電壓識別臨界機(jī)群的新方法(Approach to identify the critical cluster of machines based on line transient potential energy and voltage)[J].中國電力(Electric Power),2002,35(5):40-44.

[8] 蔡國偉,穆鋼, Chan K W,等(Cai Guowei, Mu Gang, Chan K W,etal).基于網(wǎng)絡(luò)信息的暫態(tài)穩(wěn)定性定量分析——支路勢能法(Branch potential energy method for power system transient stability assessment based on network dynamic variables)[J].中國電機(jī)工程學(xué)報(Proceedings of the CSEE),2004,24(5):1-6.

[9] Cai G W, Chan K W, Cheng H Z,etal.Transient instability detection using line transient potential energy[J].Electric Power Components and Systems, 2004, 32(6): 627-643.

[10]蔡國偉,孟祥俠,劉濤(Cai Guowei, Meng Xiangxia, Liu Tao).電力系統(tǒng)振蕩中心的暫態(tài)能量解析(Study on transient energy around power system oscillation center)[J].電網(wǎng)技術(shù)(Power System Technology),2005,29(8):30-34.

IdentificationofCriticalCut-setBasedonNetworkStructureandParameter

MENG Xiang-xia1， ZHANG You-dong1， LI Xiang-zhong2， PAN Guang-lin2

(1.Mechanical and Electrical Engineering Department, North China Institute of Science and Technology, Beijing 101601, China； 2.Thermal Power Plant of Jilin Petroleum Group Corporation, Songyuan 138000, China)

Accurately identifying critical cut-set can provide very important basis for monitor/control system allocation in network.Based on network topological structure and parameter,cut-set weighted coefficient in network is formed by means of coupling and weighted thinking,"bottleneck"section which is most liable to make system losing stability in network is identified according to magnitude of cut-set weighted coefficient,the law which the transmission network itself inherntly hierarchical sturcture characteristics affect the transient stability of power system is revealed in terms of stability.The results show that topological structure and parameter of network dominate stability of power system.Simulation results and analysis on New England 10 machines system testify the validation of the proposed method.

coupling； weighted； cut-set weighted coefficient； critical cut-set； transient stability

2010-07-29

2010-09-17

TM711

A

1003-8930(2011)01-0074-06

孟祥俠(1973-)，女，碩士，講師，主要從事電力系統(tǒng)穩(wěn)定與控制方面的教學(xué)和科研工作。Email:mengxiangxia@ ncist.edu.cn 張有東(1970-)，男，副教授，主要從事電力系統(tǒng)控制方面的教學(xué)和科研工作。Email:youdongzhang@ncist.edu.cn 李祥忠(1969-)，男，工程師，主要從事發(fā)電廠電氣分廠技術(shù)與管理工作。Email:lxz660297@163.com