刁 塑，劉俊勇，劉友波，劉挺堅，丁理杰，BAZARGAN Masoud

（1.四川大學 電氣信息學院，四川 成都 610065；2.四川電力科學研究院，四川 成都 610072；3.阿爾斯通電網研究與技術中心，英國 斯塔福德 ST17 4LX）

0 引言

近年來，電網規(guī)模日益擴大，監(jiān)視控制手段不斷加強，但在世界范圍內仍出現(xiàn)了十余起大規(guī)模停電事故，給所在國家造成了重大損失，人們在事后也對其進行了反思與總結[1-3]。相關研究[4-12]通過結合復雜網絡、自組織臨界性等新興系統(tǒng)科學理論，對大規(guī)模電網進行電氣拓撲分析與連鎖故障模擬，探究使電網發(fā)生大停電事故或趨于自組織臨界態(tài)的重要因素。在連鎖故障擴散與傳播時，電網結構與薄弱環(huán)節(jié)的狀態(tài)往往起到關鍵作用，揭示當前電網結構、元件狀態(tài)與連鎖故障之間關系的研究有一定實際意義。

目前國內外對于關鍵或脆弱元件辨識與電網連鎖故障模型的研究均較為豐富，文獻[4-6]基于電網加權拓撲模型對高介數的元件進行攻擊，獲取其對電網脆弱性的影響，對后續(xù)復雜電網研究具有啟發(fā)式的指導；文獻[7-9]提出了物理意義更加明確的關鍵線路評價指標，克服了復雜網絡最短路徑與潮流特性不符的問題；文獻[10]提出了綜合考慮概率與結構重要度的線路故障的預測指標；熵是衡量系統(tǒng)混亂程度的重要物理量；文獻[11-14]等結合電力系統(tǒng)自身特點，提出潮流熵、脆性風險熵、能量熵等廣義角度上的信息熵，從不同方面研究元件脆弱性評估與系統(tǒng)自組織臨界態(tài)辨識，但同時考慮電網結構量和狀態(tài)量隨連鎖故障過程變化的較少。在故障模型方面，很多學者利用復雜系統(tǒng)理論來解釋連鎖故障引發(fā)大停電問題[15]，提出的方法有OPA模型[16]、CASCADE模型[17]、隱性故障模型[18]等。在驗證上述脆弱元件辨識的方法時，往往使用較為簡單的故障模型，如隨機攻擊、蓄意攻擊等，且大多強調其指標合理性，應用性驗證不足，通常所提方法具有較好辨識效果，但未利用辨識結果對電網調度運行、規(guī)劃建設的啟示進行更深入的研究。連鎖故障模型一般通過調整系統(tǒng)參數如線路隨機故障概率、隱性故障概率與過載容忍因子等觀察其對電網自組織臨界性的影響，而未計及關鍵元件對連鎖故障的影響。

在建筑學方面有針對抗震結構的研究[19]，本文借鑒建筑學關于結構和承載的思想，為研究電網固有結構、元件狀態(tài)及其轉移能力與系統(tǒng)連鎖故障抵御能力之間的關系，結合復雜網絡理論建立了區(qū)別于傳統(tǒng)潮流與輸電容量的線路承載量與承載能力的數學模型，運用熵理論獲取全局承載系數分布，通過IEEE 118節(jié)點系統(tǒng)連鎖故障仿真驗證承載分布對電網自組織臨界性的影響，并針對其特點和影響在運行與建設2個技術層面提出抑制電網大停電事故的策略與方法，通過對比加入防御策略前后的實際電網連鎖故障仿真并結合某省網實際情況，驗證了本文理論的工程適用性。

1 電氣結構承載模型

根據電網網架結構、線路特性、節(jié)點負載等，提出線路承載量與承載能力的概念。與傳統(tǒng)潮流和輸電能力的區(qū)別在于，承載模型在電網被拓撲化簡的前提下，仍能反映連鎖故障過程中電網結構與線路承載狀態(tài)的變化情況。為反映電網承載均勻性，本文采用熵理論來評價電網線路承載系數分布有序程度。

1.1 網絡結構拓撲建模

應用復雜網絡理論研究電網特性前，應獲取其網絡拓撲結構，本文對電網拓撲獲取原則與假設如下。

a.節(jié)點包含2類：發(fā)電機節(jié)點集合Ng，包含電網中所有的發(fā)電機組；廣義負荷節(jié)點集合Nl，包含電網中所有變電站及負荷節(jié)點。

b.網絡中的邊集合包含電網中的所有線路和變壓器支路，權重為線路的電抗或變壓器漏抗。

c.合并操作，如合并雙回路輸電線和站內同電壓等級母線節(jié)點等；不計并聯(lián)電容支路。

1.2 電氣結構承載能力

文獻[7-8]基于電路方程提出電氣介數Be的概念，克服了傳統(tǒng)加權介數模型假設潮流只沿網絡最短路徑流動的不足。由于發(fā)電機出力和負荷值對于線路自身的傳輸能力并無直接影響，為單純反映線路承載能力大小與結構及自身參數的關系，定義電氣結構介數Bs如下：

其中，Ii（m，n）為在節(jié)點i與參考節(jié)點直接注入單位電流時，線路（m，n）上流過的電流。

電氣結構介數表示采用滿足基爾霍夫定律的全網單位電流分布值疊加來衡量線路在結構上的關鍵程度，其物理意義為僅考慮加權拓撲結構下，由發(fā)電節(jié)點往負荷節(jié)點傳輸功率時，線路所承擔的輸電任務，其值越大說明在“發(fā)電-負荷”節(jié)點對間傳輸單位功率時所攤功率越多。

輸電線路由于電壓等級和材質等原因具有不同的輸電容量，同桿多回并架線路的可靠性也未計入網絡拓撲結構中，為同時考慮上述因素，本文定義線路的承載能力 Wc（m，n）如下：

其中，F(xiàn)max（m，n）為線路（m，n）的輸電容量；N為實際線路回數。

承載能力可理解為單位電氣結構介數傳輸功率，其物理意義在于：線路輸電能力與其所承擔的輸電任務應有所對應，即越關鍵的線路所需傳輸容量越大，若不滿足這一原則，承載能力值將較小，說明該線路承擔輸電任務的能力較弱。式（2）中出現(xiàn)線路回數是由于重要輸電通道通常采用雙回及以上的多回線路進行電能輸送，而在網絡拓撲化簡時，未計及多回輸電線路容量的提升。

1.3 電氣結構承載量

當某條線路斷開后，原來傳輸的功率通常會發(fā)生轉移，受影響最大的則是該條線路兩端節(jié)點所包含的線路集合，本文基于復雜網絡理論中節(jié)點與線路的負載關系[20-21]，定義線路的承載量由其兩端節(jié)點功率度數所決定。

首先定義節(jié)點功率度數km如下：

其中，dm為節(jié)點m的傳統(tǒng)度數；Fm為與節(jié)點m相連的線路上的有功功率。

為簡單說明其物理含義，以圖1所示IEEE 9節(jié)點電網為例，說明傳統(tǒng)度數與功率度數的區(qū)別及物理意義。計算結果如表1所示（采用標幺值計算，僅列出聯(lián)絡節(jié)點）。

圖1 IEEE 9節(jié)點網絡拓撲示意圖Fig.1 Topology of IEEE 9-bus system

表1 IEEE 9節(jié)點系統(tǒng)傳統(tǒng)度數與功率度數計算結果Table 1 Calculative traditional degree and power degree for different nodes of IEEE 9-bus system

從表1可得，傳統(tǒng)度數僅能反映節(jié)點網絡拓撲信息，節(jié)點 L4、L7、L9與 L5、L6、L8這 2 組節(jié)點內部度數無差異，僅反映節(jié)點拓撲連接程度一樣。節(jié)點功率度數能夠表征節(jié)點有功吞吐密度，節(jié)點L6的功率度數大于除L9外的所有連接點，原因是其聯(lián)絡線少且負荷重，即在結構和狀態(tài)兩方面均為重載單元。功率度數較大的節(jié)點故障并退出運行，對周圍線路的影響將大于其他節(jié)點。

在復雜網絡理論中，邊的負載一定程度上受其兩端節(jié)點負載及自身權重影響[20]，本文定義線路的承載量取決于兩端節(jié)點的功率度數的乘積：

其中，θ為常數，本文中取θ=1。

文獻[20]證明了在加權網絡中，邊的負載與其兩端節(jié)點度數的乘積成線性相關，即θ取1時，是最優(yōu)負載分配原則。用功率度數km替代傳統(tǒng)節(jié)點度，可彌補純拓撲結構下度數不能反映節(jié)點電氣特性的不足。承載量是描述線路兩端節(jié)點功率吞吐密度的量，若輸電線路兩端節(jié)點的功率度數較大，意味著該線路自身負載較重或位于大功率匯集區(qū)域，是承擔電氣任務較重的元件，一旦故障則會對周圍線路造成較大影響和沖擊。

1.4 電氣結構承載熵

由上節(jié)定義可知，當系統(tǒng)發(fā)生連鎖故障并導致線路開斷和電網結構變化時，承載量和承載能力都將計及這些變化。受文獻[11-14]的啟發(fā)，本文也引入熵作為全局觀測量，通過上述能夠同時考慮結構量與狀態(tài)量變化的線路承載模型，提出衡量全網線路承載狀態(tài)分布混亂程度的電氣結構承載熵。

設線路i承載能力為Wi，實際承載量為Wci，線路集合為L，則線路i的承載系數μi為：

與傳統(tǒng)負載率相比，承載系數能反映計及電網結構與狀態(tài)參數下線路承載情況，而并非單純的潮流負載情況。承載能力并不表示最大承載量，對μi進行歸一化處理，即：

給定常數序列 U=｛U1，U2，…，Uk，…，Uv｝，v 表示數列總項數，本文取 U=｛0，0.01，0.02，…，1.00｝。 用 nk表示承載系數 μ′i?（Uk，Uk+1]的線路條數，不同承載系數區(qū)間內的線路條數概率為：

其中，P（k）為承載系數 μ′i?（Uk，Uk+1]的線路數占總線路數的比例；v-1表示所定義常數區(qū)間的間隔總數。

定義全網電氣結構承載熵為：

其中，C一般為常數，本文中C取ln10。

全網電氣結構承載熵是根據線路承載能力與信息熵理論提出的評價全網所有線路承載情況分布有序程度的指標。采用信息熵理論對原有線路承載系數進行處理后，可以得知：承載熵值越小時，全網的承載系數處于同一區(qū)間的線路條數較多，這說明整個網絡大多數線路所承擔的輸電任務跟自身輸電能力及所處位置是相符合的；承載熵值越大，那么線路承載系數則較為分散地分布在各個數值區(qū)間，說明某些線路容量很大，但其輸送功率較小且位于并不關鍵的位置，亦可能是某些線路容量較小，但輸送功率較大且處于結構關鍵位置。電網在結構上可能存在重載區(qū)域，具有潛在風險。

2 考慮隱性故障的連鎖故障仿真

為研究承載結構特性對于電網連鎖故障的影響，本文采用文獻[18]提出的基于隱性故障模型，研究系統(tǒng)處于不同初始承載熵水平下連鎖故障停電規(guī)模與其累積概率分布的關系，并統(tǒng)計負荷損失數據，揭示在連鎖故障發(fā)生過程中負荷損失情況與承載熵的關系。

2.1 隱性故障模型

隱性故障模型[18]仿真流程如下：

a.設定隱性故障模型中系統(tǒng)的各項參數，如隱性故障概率h、發(fā)電機旋轉備用比例g、線路i的最大容量系數ci、控制措施參數Tc；

b.隨機選擇一條線路發(fā)生故障，使其跳開；

c.網絡孤島檢測與處理，計算直流潮流；

d.判斷是否有線路i的潮流Fi≥Fimax，若有潮流越限的線路，則使其跳開，返回步驟c；

e.每當有線路跳開（過載或保護誤動作），由隱性故障機理確定下一次可能發(fā)生隱性故障的線路，如有線路發(fā)生隱性故障，則使其因保護誤動作而跳開，返回步驟c；

f.當跳開線路條數Ttrip≥Tc時，模擬電網進行調度并進行優(yōu)化控制，并將Ttrip置0，但仍需判斷有無隱性故障發(fā)生；

g.如無線路繼續(xù)跳開，則判定連鎖故障結束，統(tǒng)計本次仿真中負荷損失情況。

2.2 隱性故障仿真與結果分析

應用上述隱性故障模型對IEEE 118節(jié)點系統(tǒng)進行仿真，其中發(fā)電機旋轉備用取10%，線路隱性故障概率取0.02，當跳開線路達到3次時進行優(yōu)化控制，仿真次數為5000次。

2.2.1 不同初始承載熵下仿真結果與分析

為研究承載熵對電網連鎖故障的影響，本文通過調整線路最大輸電容量來獲取同一平均負載率下不同初始承載熵的系統(tǒng)狀態(tài)[12]，并研究其對電網自組織臨界性的影響。雖然在電網結構不變的情況下，通過調節(jié)負載率分布僅能使承載熵在1.5～2.1間變化，但已可看出初始承載熵的不同對電網自組織臨界性具有一定影響。

圖2（a）、（b）分別為平均負載率為 0.4和 0.6時，不同初始承載熵下停電規(guī)模累積概率分布曲線。本文僅討論在平均負載率不變的初始條件下，承載熵與停電規(guī)模累積概率分布曲線的關系。從圖中可以看出，隨著承載熵值增加，停電規(guī)模累積概率分布曲線逐漸呈現(xiàn)冪率特征。承載熵不能完全表征當前系統(tǒng)是否已處于自組織臨界態(tài)，但隨著初始承載熵的不斷增大，大停電事故累積概率分布曲線末端右移，逐漸出現(xiàn)“長尾”的現(xiàn)象，意味著系統(tǒng)發(fā)生大停電事故的概率有所增大。平均負載率為0.4時冪率特征變化現(xiàn)象較為明顯，而當平均負載率不斷提高時，承載熵所反映的趨勢有所減緩，當平均負載率達到0.8以上時，系統(tǒng)已處于極度重載的情況，無論承載量如何分布，只要有故障發(fā)生，都極易引起連鎖故障。

2.2.2 故障過程中負荷損失變化分析

取停電規(guī)模和故障線路條數均超過40%的停電事故進行統(tǒng)計，結果如圖3、圖4所示。圖3（a）、（b）分別表示平均負載率為0.4與0.6時，連鎖故障過程中潮流熵和承載熵的變化情況。圖4（a）、（b）分別表示平均負載率為0.4與0.6時，滿足篩選條件的事故每切除一條線路時的單次負荷損失量的平均結果，反映了連鎖故障過程中負荷損失的變化情況。

圖2 不同平均負載率和承載熵下停電累積概率分布Fig.2 Outage cumulative probabilistic distribution for different average power flow rates and different carrying entropies

圖3 連鎖故障過程中潮流熵和承載熵的變化趨勢Fig.3 Trend of power flow entropy and carrying entropy during cascading failure

圖4 連鎖故障過程中負荷損失變化趨勢Fig.4 Trend of load loss during cascading failure

對比連鎖故障仿真各時期承載熵、潮流熵以及負荷損失，可得以下現(xiàn)象：

①在連鎖故障初期，觀察圖4可見電網并未發(fā)生大的負荷損失，但圖3中潮流熵與承載熵都迅速增大，表明系統(tǒng)安全性正迅速下降；

②在連鎖故障中期，圖3中潮流熵值與承載熵值均處于較高水平，最大承載熵值分別為5.07和5.08，線路潮流和承載量分布較為混亂，同時，圖4中出現(xiàn)較大負荷損失，系統(tǒng)處于連鎖故障傳播擴散時期，系統(tǒng)安全性較為脆弱；

③在連鎖故障末期，電網結構較初始狀態(tài)已有較大變化，承載熵最終降低為3.51與4.20，觀察圖4可見負荷損失也逐漸減小，表明系統(tǒng)趨于穩(wěn)定，說明承載熵能夠反映結構和狀態(tài)的變化，而僅反映線路潮流分布均勻性的潮流熵則不具有這一特性，仍然處于較高水平。

2.2.3 基于隱性故障模型的仿真結果啟示

通過上述仿真，本文認為承載熵能夠反映當前系統(tǒng)的自組織臨界性質，通過電網建設或改變運行方式能夠使承載熵在一定程度內變化，在允許范圍內調整線路承載量合理分布能夠對電網發(fā)生大停電事故有一定抵制作用；另一方面，較傳統(tǒng)潮流熵而言，承載熵計及了連鎖故障過程中電網結構不斷變化的特性，較好地反映了系統(tǒng)當前對連鎖故障的抑制能力，系統(tǒng)在發(fā)生連鎖故障時，若在優(yōu)化調度中能考慮線路承載分布合理性，將有助于抑制連鎖故障事故繼續(xù)傳播。

3 考慮承載結構的運行與規(guī)劃防御策略

由仿真所得啟示，為降低電網發(fā)生連鎖故障事故概率，對大停電事故進行防御，本文考慮承載結構，提出在電網運行和建設兩方面的連鎖故障防御策略。

3.1 故障過程中的電網運行防御策略

在基于直流潮流的連鎖故障模型中，一般有計及調度行為的最優(yōu)控制策略，通常其目標函數為調整時使發(fā)電變化和負荷損失之和最小[16，22]：

其中，pg與Pgk分別為調整后與第k天時發(fā)電機節(jié)點功率注入量；pl與Plk分別為調整后與第k天時負荷節(jié)點功率注入量；S為切負荷成本折算系數，通常取較大數值以盡量保證經濟性。

本文提出在連鎖故障過程中考慮承載系數分布的電網運行最優(yōu)控制策略，使承載量更加合理地分配到線路上，在允許范圍內適度調整承載熵，讓系統(tǒng)對連鎖故障具有一定抵制作用。

由于熵是度量整個系統(tǒng)有序程度的量，計算過程中涉及對變量進行區(qū)間劃分后的概率統(tǒng)計，將其直接寫入目標函數并使其最小，可能引起大量功率轉移。借鑒熵權系數法思想，根據每條線路承載系數區(qū)間分布統(tǒng)計結果，用式（7）中常數數列區(qū)間（Uk，Uk+1]中包含線路數量nk作為該區(qū)間內線路的熵權系數，意義在于：處于同一區(qū)間的線路條數越多，說明系統(tǒng)承載量分布越有序，承載熵值越?。蝗舫休d熵值越大，則說明包含線路條數較少的區(qū)間越多，線路承載系數分布越混亂；為在一定程度上控制承載熵值，理應更傾向于調整熵權系數較低的區(qū)間內線路的潮流，即著力調整承載“受力”過重或過輕的線路。在優(yōu)化調度調整潮流時，也應滿足盡量降低潮流轉移量來減小其對系統(tǒng)沖擊的原則。

本文基于當前潮流分布，將線路熵權系數作為潮流調整量的權重，可確定在不出現(xiàn)大規(guī)模潮流轉移前提下，通過調整少量線路潮流來控制承載熵的優(yōu)化原則。將考慮熵權系數的潮流變化量計入目標函數，并作歸一化處理后，最優(yōu)控制策略修改如下：

系統(tǒng)總功率平衡：

線路潮流限值：

基爾霍夫定律：

其中，Δpgmax為節(jié)點g功率注入變化最大量；njk為線路j的熵權系數；fj為線路j調整后流過的潮流值；Fjk為當前線路潮流值；Δfjmax為潮流變化最大量；Ak由第k天的電網導納矩陣計算求得；Fmjkax為第k天線路j的最大傳輸容量；pmaxi為發(fā)電機節(jié)點i的最大功率輸出。

與傳統(tǒng)的目標函數對比，式（10）增加了對于線路潮流變化量的限制，目的是控制承載熵的變化幅度，這將削弱系統(tǒng)控制發(fā)電機出力變化和負荷損失的能力，可能造成在優(yōu)化計算時增加小量負荷損失，但卻有助于在發(fā)生故障后對于電網的承載狀態(tài)進行有效控制，降低線路承載分布混亂程度，這也是符合連鎖故障防御原則的。

3.2 電網建設中的防御策略

諸多研究表明網絡結構初始狀態(tài)與自組織臨界性也有密切聯(lián)系，因此在電網建設規(guī)劃階段也可加入對初始承載熵的控制，以減少電網發(fā)生大停電事故的概率。本文不考慮新節(jié)點和新線路的生長模型，僅根據當前負荷水平及發(fā)電機容量增長來制定考慮承載量分布的線路升級擴容方案。

考慮電網建設的連鎖故障模型通常采取固定模式或比例進行線路改造，如OPA模型在慢動態(tài)中僅針對故障斷開的線路進行擴容，改進OPA模型考慮了線路負載率的情況，在擴容前進行潮流計算，對負載率超過ε的線路進行固定比例擴容。本文認為在電網規(guī)劃建設過程中，可考慮承載量分布情況，制定更合理的擴容方案，原因在于負載率超過ε的線路承載情況亦可能有較大差異，如處于大功率匯集區(qū)的線路比普通線路更易引起連鎖故障，擴容需求更大。

為與傳統(tǒng)擴容方案進行對比，同時消除不同擴容策略所確定擴容總量不同帶來的影響，本文基于OPA模型中慢動態(tài)建設過程提出考慮承載熵的與傳統(tǒng)方案等量擴容的主動防御策略。

根據擴容需求負載率ε確定需擴容線路集合E，按容量增長因子λ確定與傳統(tǒng)擴容方案相等的擴容總量Q。

負荷變化限值：

機組出力限值：

其中，F(xiàn)jkmax為第k天線路j的最大傳輸容量。

由線路承載系數與熵權系數確定集合E內各條線路的擴容分配值Gj，其意義在于承載較重且熵權系數小的線路擴容需求大且對承載熵影響也大，按擴容需求進行線路擴容，可在一定程度上控制承載系數分布，緩解線路承載壓力。

各條線路擴容增量由擴容總量按各條線路擴容分配值占總分配值的比例分攤確定。

其中，ΔFmaxj（k+1）為第 k+1天線路 j的最大傳輸容量的增量。

4 計及防御策略后的連鎖故障仿真

4.1 加入防御策略的連鎖故障模型

上節(jié)提出了在電網運行和建設2個方面的防御策略，為將其應用于連鎖故障仿真，需將第2節(jié)中的連鎖故障模型擴展為帶有電網建設過程的模型。本文在隱性故障模型和改進OPA模型[22]的基礎上，結合所提防御策略，提出在連鎖故障和電網建設過程中加入對承載熵的控制策略的連鎖故障仿真模型。式（10）—（15）組成運行防御策略并計入內層快動態(tài)過程，式（16）—（18）組成建設防御策略并計入外層慢動態(tài)過程，具體仿真流程如圖5所示。

圖5 加入防御策略的連鎖故障仿真流程圖Fig.5 Flowchart of cascading failure simulation with defense strategy

a.設定連鎖故障模型中的各項參數：線路隨機故障概率τ，負荷增長率β，線路容量增長因子δ，保護裝置誤動概率σ，保護裝置拒動概率ξ，擴容需求負載率ε，優(yōu)化控制成功概率η，負荷損失代價因子L，慢動態(tài)天數d。

b.執(zhí)行快動態(tài)過程，模擬由于隨機因素導致線路以概率τ開斷，孤島檢測處理后，計算直流潮流，根據線路越限情況與優(yōu)化控制成功概率η按式（10）—（15）進行優(yōu)化調度。

c.根據保護動作概率確定是否切除線路。

d.若有線路被切除，則進入步驟b，否則統(tǒng)計快動態(tài)過程數據，進入步驟e。

e.執(zhí)行慢動態(tài)過程，若當前天數k小于設定天數d，則根據負荷增長因子β確定第k+1天負荷水平及發(fā)電機最大出力，否則仿真結束。

f.進行直流潮流計算與優(yōu)化計算，并根據線路功率越限與切負荷情況，按備用容量比例確定發(fā)電機出力分攤增量。

g.再次進行潮流計算，按式（16）—（18）對負載率超過ε的線路進行升級擴容，模擬電網規(guī)劃與建設，慢動態(tài)過程結束，轉入步驟b。

4.2 實際電網連鎖故障模擬

根據本文所提加入防御策略后的連鎖故障模型，對西南某省實際電網按第4.1節(jié)所述流程進行仿真。觀察加入計及承載結構的防御策略改善承載分布后對于電網大停電事故累積概率的影響，并針對該省實際情況，對線路承載進行排序與分析。

該省電網含680個節(jié)點（含發(fā)電機節(jié)點與廣義負荷節(jié)點）、977條支路（含輸電線路與等效變壓器支路），2012年峰值負荷約為19907MW，2013年峰值負荷約為23324MW，可計算全年負荷增長率為0.1716，折算可得快動態(tài)中日負荷增長率約為1.00047，模型中參數見表2。

表2 連鎖故障仿真模型參數Table 2 Parameters of cascading failure simulation model

圖6為加入主動防御策略前后停電規(guī)模累積概率分布曲線，從圖中看出，加入防御策略后，小規(guī)模停電事故概率雖有略微上升，但大規(guī)模停電事故概率有所改善，在一定程度上說明計及承載熵的防御策略利用少量負荷損失作為代價，對承載分布進行有效控制，增加系統(tǒng)對于連鎖故障的抑制能力，在系統(tǒng)承受范圍內降低了大規(guī)模停電事故發(fā)生的概率。

另一方面，對電網線路承載情況進行統(tǒng)計，列出承載較重的線路如表3所示，并繪制如圖7所示的全網承載系數分布圖，用以觀察電網線路承載分布情況。通過表3可知，從線路負載率來看，500 kV與220 kV輸電線路負載率并不高，甚至某些線路的負載率很低，但其承載系數較大，說明其在電氣結構上處于“承載受力”較重的區(qū)域或自身“承載受力”較重，而傳統(tǒng)方法從負載率角度出發(fā)則無法反映這一現(xiàn)象。觀察圖7可得，當前全網承載分布有很不均勻的區(qū)域，正好對應表3中大部分線路所在區(qū)域，該區(qū)域是該省網的重要負荷區(qū)域，屬于存在大量負荷以及大功率匯集的地區(qū)，與本文第1節(jié)對于線路承載量的定義相符合。

表3 承載系數排名前10的線路Table 3 Carrying coefficient of top ten lines

圖7 某省電網承載系數分布圖Fig.7 Carrying coefficient distribution of a provincial power grid

圖8所示為該省網500 kV網架結構圖，從圖8中可知，譚家、龍王、桃鄉(xiāng)、丹景、蜀州等節(jié)點及其之間的線路構成了該省網南北輸電斷面，且負荷較重，是電氣結構上的重載區(qū)域。除尖山與桃鄉(xiāng)的聯(lián)絡線外，其余線路承載系數均排名靠前。從圖中可以看出，尖山—桃鄉(xiāng)線路雖處于關鍵結構位置，但因其有3條并架回路，保證其傳輸能力與其關鍵程度相匹配，故不是“受力”較重線路，與本文第2節(jié)對于線路承載系數的定義相符。在慢動態(tài)過程中，對線路進行升級擴容時，本文策略考慮了以上因素，同時計及負載率與承載系數再進行擴容，使電網更加堅強，也在一定程度上提升了系統(tǒng)對于連鎖故障的抵御能力，具有一定合理性。

圖8 某省500 kV網架結構Fig.8 Topology of a provincial 500 kV power grid

5 結論

本文在分析電網復雜網絡結構和電氣狀態(tài)的基礎上，結合熵理論提出揭示全網承載分布特性的承載熵，并分析了承載分布對于電網連鎖故障的影響，仿真表明：系統(tǒng)初始承載熵值能夠在一定程度上反映當前電網狀態(tài)，初始電網線路承載系數分布與大停電發(fā)生概率具有一定關聯(lián)性；在電網運行過程中，承載熵值變化與負荷損失變化情況一致，說明計及電網結構特性的承載熵在電網架構發(fā)生較大變化時仍能夠較好地反映當前系統(tǒng)安全性。

基于以上2點啟示，本文提出了同時在電網運行和規(guī)劃建設環(huán)節(jié)中針對連鎖故障的防御策略，對電網發(fā)生大規(guī)模停電事故進行預防控制。針對實際電網仿真表明：采用文中所提的2個層面的防御策略后，大停電事故概率有所下降，考慮承載系數的擴容策略可使全網承載分布更合理，將阻止少量關鍵或脆弱線路成為連鎖故障觸發(fā)源，驗證了本文指標體系與防御策略具有一定的工程適用性。