何海丹，賈燕冰，劉睿瓊

(太原理工大學(xué)電氣與動力工程學(xué)院，太原市030024)

0 引言

近年來世界范圍內(nèi)的大停電事故頻繁發(fā)生，經(jīng)分析得出少數(shù)關(guān)鍵線路故障是誘導(dǎo)連鎖故障發(fā)生的主要因素［1-2］，因此，有效辨識電力系統(tǒng)關(guān)鍵線路能預(yù)防災(zāi)難性事故的發(fā)生，繼而提高電力系統(tǒng)穩(wěn)定性。

辨識系統(tǒng)關(guān)鍵線路的指標(biāo)有潮流介數(shù)［3-6］、潮流轉(zhuǎn)移度［7］、電氣介數(shù)［8］、最大流判據(jù)［9］等。文獻(xiàn)［3-4］根據(jù)線路被母線間最短路徑所經(jīng)過的次數(shù)反映該線路被利用情況作為判斷該線路是否為關(guān)鍵線路的指標(biāo);文獻(xiàn)［5-6］基于潮流分布理論定義潮流介數(shù)指標(biāo)，衡量輸電線路在發(fā)電和負(fù)荷間功率傳輸?shù)倪^程中被利用的程度，繼而識別系統(tǒng)關(guān)鍵線路。針對在系統(tǒng)故障的同時常伴隨潮流大范圍轉(zhuǎn)移的特征，文獻(xiàn)［7］提出潮流轉(zhuǎn)移度指標(biāo);文獻(xiàn)［8-9］通過疊加系統(tǒng)支路電流絕對值分量來評估電網(wǎng)結(jié)構(gòu)的關(guān)鍵環(huán)節(jié)。

文獻(xiàn)［3-4］的模型和方法都是假設(shè)母線間功率只按最短路徑流動(包括電氣距離)，忽略了電網(wǎng)節(jié)點(diǎn)間功率按Kirchoff定律傳輸?shù)氖聦?shí)，難以反映故障的嚴(yán)重性。文獻(xiàn)［5-9］雖克服了上述不足，但指標(biāo)的計算需要系統(tǒng)潮流分布數(shù)據(jù)，有一定的局限性，且上述文獻(xiàn)都是單方面在系統(tǒng)結(jié)構(gòu)上辨識關(guān)鍵線路，未考慮支路風(fēng)險性，不夠全面。

鑒于上述問題，本文基于功率分布理論定義發(fā)電機(jī)－線路功率分布因子和負(fù)荷－線路功率分布因子，綜合系統(tǒng)內(nèi)“發(fā)電－負(fù)荷對”定義支路功率分布因子介數(shù)指標(biāo)，從系統(tǒng)結(jié)構(gòu)上識別線路重要度;綜合功率分布因子介數(shù)和支路概率重要度，提出支路綜合重要度指標(biāo)，其從結(jié)構(gòu)和概率2個方面識別系統(tǒng)支路的重要度。

上述文獻(xiàn)關(guān)鍵線路識別指標(biāo)的計算需要以系統(tǒng)線路潮流分布和系統(tǒng)當(dāng)前運(yùn)行方式為基礎(chǔ)，本文提出的指標(biāo)計算只需要系統(tǒng)基本參數(shù)，無須考慮潮流計算速度及收斂性等問題，且可提前識別系統(tǒng)關(guān)鍵線路，更簡單快捷。

1 結(jié)構(gòu)重要度

基于電流分布理論，得出發(fā)電機(jī)與負(fù)荷間直接的功率傳輸關(guān)系［10］。

設(shè)電力網(wǎng)絡(luò)具有n個節(jié)點(diǎn)，q個電源，l條支路。在潮流斷面確定的情況下，電源以電流源描述記為根據(jù)功率分布理論，由電網(wǎng)節(jié)點(diǎn)電壓方程，可得:

式中:Zik為阻抗矩陣 Z 中的元素;為第i個節(jié)點(diǎn)的電壓;為第m個電源的電流;為電流源的共軛。

為了區(qū)分指定電源和非指定電源引入變量m，其代表的意義和k相同。

由電源k提供的功率S·iji←ks為

1.1 發(fā)電－支路功率分布因子

定義發(fā)電－支路功率分布因子

正常運(yùn)行時系統(tǒng)節(jié)點(diǎn)電壓標(biāo)幺值約等于1，近似取Ui=Uk=1，得到:

1.2 負(fù)荷－支路功率分布因子

根據(jù)發(fā)電－支路功率分布因子可知任意支路中由某臺發(fā)電機(jī)所提供的功率。任意一臺發(fā)電機(jī)在某條支路上的分布功率計算式為

支路的總傳輸功率等于各個發(fā)電機(jī)對其提供功率的疊加，即

類似式(2)、(3)，可得出負(fù)荷m消耗的功率和由i端流入支路ij中的傳輸功率向負(fù)荷m提供的功率，同式(4)、(5)推導(dǎo)過程，可得負(fù)荷 － 支路功率分布因子為

式中:i為支路的注入節(jié)點(diǎn);m為負(fù)荷節(jié)點(diǎn)。

1.3 功率分布因子介數(shù)

定義支路功率分布因子介數(shù):已知電網(wǎng)系統(tǒng)中有

個電源，f個負(fù)荷，則系統(tǒng)有q×f個發(fā)電－負(fù)荷對，定義發(fā)電負(fù)荷對(Gk，Lm)在支路ij上的發(fā)電－負(fù)荷對功率分布因子Fij(k→m)，即支路ij的功率中由電源Gk提供的且最終流入負(fù)荷Lm的功率占支路ij總傳輸功率的比值:

對于和發(fā)電機(jī)和負(fù)荷直連節(jié)點(diǎn)的線路權(quán)重定義為

式中:Sks是單臺發(fā)電機(jī)的發(fā)電量為系統(tǒng)總發(fā)電量是負(fù)荷m對支路ij汲取的功率;G集是和發(fā)電機(jī)直連線路集合;L集是和負(fù)荷直連線路集合。

電網(wǎng)中每條支路承擔(dān)的不同發(fā)電－負(fù)荷對之間功率傳輸任務(wù)越大，可認(rèn)為其在電網(wǎng)結(jié)構(gòu)上的重要程度越高。

綜合各個發(fā)電負(fù)荷對定義支路功率分布因子介數(shù)為

支路結(jié)構(gòu)重要度由電網(wǎng)中所有發(fā)電機(jī)發(fā)出功率、所有負(fù)荷消耗的功率、支路傳輸功率、經(jīng)過支路的發(fā)電－負(fù)荷對功率分布因子絕對值之和，以及線路傳輸容量裕度共同決定。

以支路發(fā)電負(fù)荷對功率分布因子絕對值之和描述其在功率傳輸所承擔(dān)的作用，以支路傳輸容量裕度描述系統(tǒng)的受干擾能力，兩者的比值描述支路的結(jié)構(gòu)重要度。支路結(jié)構(gòu)重要度因子越大，說明該支路對于電網(wǎng)結(jié)構(gòu)穩(wěn)定性越重要，其在電網(wǎng)中的地位越重要。

2 概率重要度

經(jīng)驗(yàn)表明，線路負(fù)載率越高，傳輸功率波動越大，越容易發(fā)生故障［11］。本文根據(jù)文獻(xiàn)［11］提出的支路概率重要度指標(biāo)，選取線路故障的影響因素為:線路的負(fù)載率、線路隱性故障概率、線路潮流變化量、線路歷史故障率。

線路概率重要度定義為

式中:mij是支路ij負(fù)載率;χij是支路ij的潮流變化率;hij是支路ij隱性故障率;αij是支路ij歷史故障率;o有2種狀態(tài)0、1，0表示初始狀態(tài)，1表示非初始狀態(tài);w1、w2、w3、w4是各項(xiàng)指標(biāo)的權(quán)重。

3 綜合重要度

定義支路綜合重要度Yij為

支路結(jié)構(gòu)重要度或支路概率重要度越大，說明支路對于提高系統(tǒng)穩(wěn)定性地位越重要。本文將支路的結(jié)構(gòu)重要度及概率重要度結(jié)合確定支路綜合重要度因子，該因子越大，說明該支路在電網(wǎng)中發(fā)生故障的概率越高，故障后影響范圍越大。

4 算例

在Matlab平臺上搭建模型，以IEEE 39節(jié)點(diǎn)為例進(jìn)行仿真，該系統(tǒng)共有10臺發(fā)電機(jī)、19個負(fù)荷點(diǎn)和46條線路，其拓?fù)浣Y(jié)構(gòu)如圖1所示。

圖1 IEEE 39節(jié)點(diǎn)系統(tǒng)接線圖Fig.1 IEEE 39-bus system wiring diagram

4.1 功率分布介數(shù)

根據(jù)公式(11)將46條輸電線路的線路介數(shù)降序排列后繪制于圖2。圖中橫坐標(biāo)表示輸電線路編號，縱坐標(biāo)為對應(yīng)輸電線路的線路功率分布介數(shù)。

如圖所示，按本文方法選取線路功率分布介數(shù)指標(biāo)排序在前10位的線路，并將結(jié)果與已有模型［10］進(jìn)行對比，如表1所示。

文獻(xiàn)［10］基于電路方程，克服了加權(quán)介數(shù)模型假設(shè)母線間潮流只沿最短路徑流動的不足，能有效反映所有“發(fā)電－負(fù)荷節(jié)點(diǎn)對”對線路的真實(shí)利用情況，從排序在前10位的線路中得出線路16－19(27)、16－17(26)、15－16(25)、2－3(3)和 26－27(42)共同屬于關(guān)鍵線路。由于本文還考慮了和發(fā)電機(jī)、負(fù)荷直連的線路的重要性，線路46、17都屬于和發(fā)電機(jī)直聯(lián)的線路。例如線路29－38(46)不加發(fā)電機(jī)權(quán)重時，潮流量是8.166，介數(shù)排序第15位，不屬于關(guān)鍵線路，鑒于該線路和發(fā)電機(jī)直連，考慮權(quán)重后，線路介數(shù)為9.242，介數(shù)排序第9，提高了其關(guān)鍵性。線路3－39(17)同理分析。線路32、6排序靠前的原因是對應(yīng)節(jié)點(diǎn)上消耗的負(fù)荷量大(分別是680 MW、822 MW)，且是系統(tǒng)中最大兩處負(fù)荷。例如線路19－20(32)，將20節(jié)點(diǎn)負(fù)荷等效成節(jié)點(diǎn)對地導(dǎo)納前，線路的線路介數(shù)是1.963，介數(shù)排序第41，考慮權(quán)重后，線路介數(shù)是 14.291，介數(shù)排序第 2。線路3－4(6)同理分析可得。

圖2 IEEE 39節(jié)點(diǎn)系統(tǒng)線路功率分布介數(shù)Fig.2 Line power distribution betweenness of IEEE 39-bus system

表1 支路功率分布介數(shù)排序Table 1 Branch power distribution betweenness sorting

切負(fù)荷量是反映系統(tǒng)可靠性的指標(biāo)，本文選取最優(yōu)切負(fù)荷量［12］作為驗(yàn)證電力系統(tǒng)受破壞程度，為了避免因線路潮流上限值設(shè)定的不同而影響切負(fù)荷量，這里將所有線路的潮流上限都設(shè)定為650 MW。

分別斷開2種方法得出的排序前7位結(jié)構(gòu)重要度高的線路，計算線路開斷后的系統(tǒng)最優(yōu)切負(fù)荷量如圖3所示。

圖3 2種方法切負(fù)荷量對比Fig.3 Load-shedding comparison in two methods

由圖3可見:按本文方法得出的切負(fù)荷折線均高于文獻(xiàn)［10］的切負(fù)荷折線，說明基于功率分布介數(shù)的切負(fù)荷量比基于文獻(xiàn)［10］的電氣介數(shù)的切負(fù)荷量高，則本文方法能夠更好地評估系統(tǒng)支路重要度，且本文指標(biāo)無需潮流計算，更簡單快速。

4.2 支路綜合重要度指標(biāo)

將本文功率分布介數(shù)指標(biāo)和文獻(xiàn)［11］提出的支路概率重要度指標(biāo)相結(jié)合，計算初始系統(tǒng)的綜合重要度指標(biāo)，如表2所示。

表2 初始狀態(tài)系統(tǒng)支路綜合重要度排序Table 2 Branch integrated importance sorting in the initial state

從表2可知，支路16－19是綜合重要度最高的線路。在模擬一階故障時，選取支路16－19開斷，綜合重要度較高的支路是16－17和16－21，都是和故障線路有相同節(jié)點(diǎn)的線路，16－。19斷線后，切斷了發(fā)電機(jī)33、34向大系統(tǒng)供電的路徑，影響發(fā)電機(jī)32周圍支路潮流的變換，因此相比較支路10－13和13－14在新運(yùn)行方式下的功率分布介數(shù)要高，使得綜合重要度值也高。表3為一階故障下支路綜合重要度排序。

表3 一階故障下支路綜合重要度排序Table 3 Branch integrated importance sorting under first-order fault

5 結(jié)論

本文以支路重要度來衡量各個支路的關(guān)鍵性，基于發(fā)電負(fù)荷對功率分布因子提出支路功率介數(shù)指標(biāo)，從結(jié)構(gòu)上識別關(guān)鍵線路。為了更全面地評估系統(tǒng)支路重要度，從結(jié)構(gòu)和概率兩方面綜合考慮，提出支路綜合重要度指標(biāo)。通過IEEE 39節(jié)點(diǎn)算例表明，采用本文所提方法，可以快速準(zhǔn)確地識別電力系統(tǒng)中的關(guān)鍵線路，尤其對于規(guī)模龐大的互聯(lián)大電網(wǎng)來說，潮流計算收斂難，所需時間長，本文所提方法可避免進(jìn)行潮流計算，具有很大的應(yīng)用前景。

［1］張文亮，周孝信，印永華，等．華北–華中–華東特高壓同步電網(wǎng)構(gòu)建和安全性分析［J］．中國電機(jī)工程學(xué)報，2010，30(16):1-5．Zhang Wenliang，Zhou Xiaoxin，Yin Yonghua，et al．Composition and security analysis of"North China-Central China-East China"UHV synchronous power grid［J］．Proceedings of the CSEE，2010，30(16):1-5．

［2］梁志峰，葛睿，董昱，等．印度“7.30”、“7.31”大停電事故分析及對我國電網(wǎng)調(diào)度運(yùn)行工作的啟示［J］．電網(wǎng)技術(shù)，2013，37(7):1841-1848．Liang Zhifeng，Ge Rui，Dong Yu，et al．Analysis of large-scale blackout occurred on July 30 and July 31，2012 in India and its lessons to China's power grid dispatch and operation［J］．Power System Technology，2013，37(7):1841-1848．

［3］王亮，劉艷，顧雪平，等．綜合考慮節(jié)點(diǎn)重要度和線路介數(shù)的網(wǎng)絡(luò)重構(gòu)［J］．電力系統(tǒng)自動化，2010，34(12):29-33．Wang Liang，Liu Yan，Gu Xueping，et al．Skeleton-network reconfiguration based on node importance and line betweenness［J］．Automation of Electric Power Systems，2010，34(12):29-33．

［4］ Chen G，Dong Z Y，Hill D J，et al．An improved model for structural vulnerability of analysis of power networks［J］．Physica A:Statistical Mechanics and Its Applications，2009，388(19):4259-4266．

［5］梁才，劉文穎，但揚(yáng)清，等．輸電線路的潮流介數(shù)及其在關(guān)鍵線路識別中的應(yīng)用［J］．電力系統(tǒng)自動化，2014，38(8):35-40．Liang Cai，Liu Wenying，Dan Yangqing，et al．Flow betweenness of transmission lines and its application in critical line identification［J］．Automation of Electric Power Systems，2014，38(8):35-40．

［6］劉文穎，梁才，徐鵬，等．基于潮流介數(shù)的電力系統(tǒng)關(guān)鍵線路辨識［J］．電機(jī)工程學(xué)報，2013，33(31):90-98．Liu Wenying，Liang Cai，Xu Peng，et al．Identification of critical line in power systems based on flow betweenness［J］．Proceedings of the CSEE，2013，33(31):90-98．

［7］曾文凱，文勁宇，程時杰，等．復(fù)雜電網(wǎng)連鎖故障下的關(guān)鍵線路辨識［J］．電機(jī)工程學(xué)報，2014，34(7):1103-1112．Zeng Wenkai， Wen Jinyu， Cheng Shijie， et al． Critical line identification of complex power system in cascading failure［J］．Proceedings of the CSEE，2014，34(7):1103-1112．

［8］鞠文云，李銀紅．基于最大流傳輸貢獻(xiàn)度的電力網(wǎng)關(guān)鍵線路和節(jié)點(diǎn)辨識［J］．電力系統(tǒng)自動化，2012，9(1):6-12．Ju Wenyun，Li Yinhong．Identification of critical lines and nodes in power grid based on maximum flow transmission contribution degree［J］．Automation of Electric Power Systems，2012，9(1):6-12．

［9］徐林，王秀麗，王錫凡．電氣介數(shù)及其在電力系統(tǒng)關(guān)鍵線路識別中的應(yīng)用［J］．中國電機(jī)工程學(xué)報，2010，30(1):33-39．Xu Lin，Wang Xiuli，Wang Xifan．Electric betweenness and its application in vulnerable line identification in power system［J］．Proceedings of the CSEE，2010，30(1):33-39．

［10］顏麗，鮑海．基于電流分布的電網(wǎng)功率分布因子的計算［J］．中國電機(jī)工程學(xué)報，2010，31(1):80-85．Yan Li，Bao Hai．Algorithm of power distribution factor based on current distribution［J］．Proceedings of the CSEE，2010，31(1):80-85．

［11］王濤，王興武，顧雪平，等．基于概率及結(jié)構(gòu)重要度的電力系統(tǒng)事故鏈模型與仿真［J］．電力自動化設(shè)備，2013，33(7):51-57．Wang Tao，Wang Xingwu，Gu Xueping，et al．Critical lines identification of power grid based on the fault chain model［J］．High Voltage Engineering，2014，40(8):2306-2323．

［12］李文沅．電力系統(tǒng)風(fēng)險評估模型、方法和應(yīng)用［M］．北京:科學(xué)出版社，2006．