基于改進(jìn)自組織映射聚類算法的連鎖故障預(yù)測模型

宋玉琴， 趙 攀， 周琪瑋

(西安工程大學(xué)電子信息學(xué)院， 西安 710600)

近年來，電力系統(tǒng)逐漸向復(fù)雜化、集成化的方向發(fā)展，使得連鎖故障發(fā)生的概率大大增加，給社會經(jīng)濟(jì)發(fā)展、人民生產(chǎn)生活等許多方面帶來了嚴(yán)重的后果[1-3]。如2015年土耳其“3·31”大停電造成除部分電力供應(yīng)由伊朗提供的凡城(Van)和哈卡里(Hakkari)地區(qū)外，土耳其的電力供應(yīng)幾乎全部中斷，影響人數(shù)約7 000 萬，占全國人口總數(shù)90%，共損失負(fù)荷約32 950 MW，經(jīng)濟(jì)損失約7 億美元[4]。研究表明：由簡單故障引起失去功角穩(wěn)定或電壓穩(wěn)定，從而導(dǎo)致大面積停電事故的發(fā)生，這一原因已經(jīng)越來越少[5]。絕大多數(shù)大停電事故是由于系統(tǒng)中某條線路因故停運，引起潮流轉(zhuǎn)移與重新分配，造成其他相關(guān)線路相繼過載跳閘，形成連鎖故障[6]。因此，加強(qiáng)連鎖故障機(jī)理分析，快速、準(zhǔn)確地提前預(yù)測連鎖故障可能的發(fā)展路徑，對及時采取有效的阻斷措施從而預(yù)防大面積停電事故的發(fā)生具有重要意義。

目前，針對這一問題的研究，已經(jīng)提出了一些方法[7-11]。羅毅等[12]將簡單事故鏈理論應(yīng)用于連鎖故障的預(yù)測中，建立連鎖故障事故鏈預(yù)測模型。該方法雖然能得到連鎖故障的傳播路徑，但由于事故鏈初始故障指標(biāo)和上下級支路間關(guān)聯(lián)性指標(biāo)不完善；具體表現(xiàn)為指標(biāo)單一，不能全面反映故障的引起因素，導(dǎo)致預(yù)測準(zhǔn)確度不高。何曄等[13]為實現(xiàn)連鎖故障預(yù)測的計算量和預(yù)測事故鏈數(shù)量的平衡，首次將C均值聚類算法引入預(yù)測模型中。與以前的方法相比，該方法能在工作量與預(yù)測結(jié)果完備性之間做出較好的平衡[14]。但是C均值聚類的效果差，不能對關(guān)聯(lián)性指標(biāo)進(jìn)行準(zhǔn)確分類。丁明等[15]提出將靜態(tài)故障樹改為動態(tài)故障樹，該提議更符合大停電事故的發(fā)展規(guī)律，但使得算法復(fù)雜度大大增加，不利于推廣。為了減少算法的復(fù)雜度，馬至遠(yuǎn)等[14]在分析大停電事故原理的基礎(chǔ)上，提出只需關(guān)注關(guān)鍵線路發(fā)生故障的可能性，該提議在很大程度上提高了模型的預(yù)測速度，也符合大多數(shù)情況下大停電事故的發(fā)展規(guī)律，但這往往會忽略某些非關(guān)鍵線路突發(fā)故障而引起的大停電事故。

為解決上述問題，現(xiàn)首先設(shè)計初始故障指標(biāo)評價體系、線路間關(guān)聯(lián)性指標(biāo)評價體系，有效地解決指標(biāo)混亂、單一等問題。其次，引入改進(jìn)SOM聚類算法，對中間環(huán)節(jié)指標(biāo)進(jìn)行聚類分類；在保證分類準(zhǔn)確的基礎(chǔ)上不斷優(yōu)化算法，確保預(yù)測的故障路徑全面而不冗余。最終，以IEEE36節(jié)點模型為算例進(jìn)行仿真分析，驗證本文方法的準(zhǔn)確性和全面性。

1 故障預(yù)測模型的指標(biāo)設(shè)計

1.1 初始故障集指標(biāo)的設(shè)計

1.1.1 節(jié)點擾動對支路的潮流沖擊性指標(biāo)

電力系統(tǒng)在正常工作時處于平衡狀態(tài)，支路l的潮流為Pl0，當(dāng)某節(jié)點(設(shè)為a節(jié)點)負(fù)荷突變，引起系統(tǒng)擾動，支路l的潮流改變?yōu)镻la，節(jié)點a受到擾動后，支路l的潮流增加量(節(jié)點a對支路l的潮流沖擊)為

ΔEla=Pla-Pl0(1)

則節(jié)點a對系統(tǒng)的潮流沖擊為

式(2)中：S為系統(tǒng)支路總數(shù)。

用支路l的潮流沖擊率ηla來表示支路l在節(jié)點a對系統(tǒng)的潮流沖擊中所占的比例，即

結(jié)合式(1)～式(3)，定義節(jié)點a在系統(tǒng)的潮流分布熵為

HD(a)的大小反映了節(jié)點在發(fā)生擾動后系統(tǒng)所受到的潮流沖擊的分布特性。結(jié)合式(4)，可以得出線路l在系統(tǒng)節(jié)點受到擾動后，自身所承擔(dān)的潮流沖擊為

由式(5)可知，支路附近節(jié)點受到的擾動對系統(tǒng)潮流沖擊越大，節(jié)點的潮流分布熵越小，系統(tǒng)的潮流沖擊就越容易集中在節(jié)點附近的極少數(shù)線路上，這幾條線路就更容易發(fā)生故障。

1.1.2 輸電線路重要性指標(biāo)

一般情況下，在發(fā)電機(jī)Gm向負(fù)荷Ln傳送功率的過程中，線路P(m,n)承擔(dān)的部分越多，就表明該支路在電力系統(tǒng)中占有比較重要的位置，表明其在電力系統(tǒng)中的重要性比較高。定義一對發(fā)電負(fù)荷對(m,n)的線路潮流介數(shù)為

式(6)中：Pij(m,n)為總傳輸功率P(m,n)在線路i、j上的分量；Pij(m,n)/P(m,n)為線路i、j對發(fā)電負(fù)荷對(m,n)的輸電貢獻(xiàn)比例因子；min(Sm,Sn)為該發(fā)電負(fù)荷對的權(quán)重，該權(quán)重取實際發(fā)電量和實際負(fù)荷消耗量中的較小值，表征了輸電線路所傳輸?shù)淖畲笥杏霉β省?/p>

電網(wǎng)中每條線路都承擔(dān)了不同發(fā)電負(fù)荷對之間的功率傳送任務(wù)，從而一條線路可能被多個發(fā)電負(fù)荷對的功率傳輸路徑經(jīng)過，其被經(jīng)過的次數(shù)也反映了該線路在電網(wǎng)當(dāng)前運行方式下的重要程度[16]。將經(jīng)過該支路的所有發(fā)電負(fù)荷對所傳輸?shù)墓β蔬M(jìn)行累加，定義線路的重要性指標(biāo)為

式(7)中：G為給線路提供功率的所有發(fā)電機(jī)節(jié)點集合；L為消耗線路功率的所有負(fù)荷節(jié)點集合。

1.1.3 支路自身故障概率指標(biāo)

查閱IEEE36節(jié)點系統(tǒng)歷史運行數(shù)據(jù)，以及多次對系統(tǒng)進(jìn)行仿真實驗，得到支路自身故障概率λl。該故障主要包括線路硬件失效、線路保護(hù)或斷路器的拒動、誤動等一些小概率故障。

1.1.4 初始故障集的生成

綜合上述指標(biāo)，在每個指標(biāo)中都加入合適的閾值，定義支路l的初始故障綜合指標(biāo)為

El=η1Clη2Filη3λl(8)

式(8)中：η1、η2、η3分別為3個初始指標(biāo)的閾值。本文的3個指標(biāo)從不同的角度分析了可能引起連鎖故障的原因。為保證預(yù)測結(jié)果符合實際情況而加入閾值，在閾值的確定過程中，將每次模型的預(yù)測結(jié)果與電力系統(tǒng)分析綜合程序(power system analysis synthesis program，PSASP)中線路真實故障情況相對比，若有不同，分析原因，修改初始故障集的閾值。閾值的確定運用理論和實驗相結(jié)合的方法。由理論分析可得，指標(biāo)2的影響程度高于指標(biāo)1和指標(biāo)3。因為當(dāng)重要線路發(fā)生故障時，連接的線路眾多，對后續(xù)影響更嚴(yán)重，而指標(biāo)1和指標(biāo)3不一定會觸及到重要線路。因此，確定閾值為η1=0.3，η2=0.5，η3=0.2；代入預(yù)測模型中，預(yù)測結(jié)果與線路實際故障的結(jié)果有很大的不同(邊緣線路過多而重要線路過少，導(dǎo)致預(yù)測的故障線路條數(shù)偏少，不能反映真實的情況)。因此，提高指標(biāo)2的權(quán)重的實驗結(jié)果表明：接負(fù)載的線路過多，而重要線路依舊不全面。繼續(xù)提高指標(biāo)2的權(quán)重，實驗結(jié)果表明：幾乎全部為重要線路，而非重要線路發(fā)生故障的概率很小。經(jīng)過多次試驗，最終設(shè)定為η1=0.2，η2=0.7，η3=0.1；此時模型預(yù)測結(jié)果和實際的故障結(jié)果接近。

1.2 上下級支路間關(guān)聯(lián)性指標(biāo)的建立

1.2.1 上下級支路間耦合性指標(biāo)

定義在故障傳播中，在上級故障支路l-1發(fā)生故障前后，下級支路l的傳輸功率變化與上級支路l-1在故障前傳輸功率的比值為上級支路故障切除后對下級支路的影響度，即

1.2.2 支路在電網(wǎng)運行中的負(fù)載率指標(biāo)

支路在電網(wǎng)運行中的負(fù)載率指標(biāo)定義為支路l在支路l-1發(fā)生故障且將要傳播到下一級時，其傳輸功率與其所承擔(dān)的最大傳輸功率的比值。

式(10)中：Flmax為支路l的最大傳輸功率。

1.2.3 支路潮流轉(zhuǎn)移熵沖擊性指標(biāo)

該指標(biāo)為了衡量上級支路故障開斷造成的潮流轉(zhuǎn)移對下級支路的影響。具體定義如下：

設(shè)支路l-1開斷對支路l的傳輸裕度影響值ηl,l-1為

因此考慮支路l-1開斷對支路l傳輸裕度影響的潮流轉(zhuǎn)移熵Hl,l-1為

Hl,l-1=-ηl,l-1lnηl,l-1(12)

1.2.4 支路間關(guān)聯(lián)性綜合評價指標(biāo)

綜合上述評價指標(biāo)，定義上下級支路間關(guān)聯(lián)性指標(biāo)為

Ml=Sl×Dl×Bl(14)

1.3 事故鏈搜索結(jié)束判據(jù)

1.3.1 電力系統(tǒng)解列

電力系統(tǒng)由于故障退出運行的線路過多，系統(tǒng)無法繼續(xù)同步運行，會自動解列成各自同步運行的幾個部分，對于大電網(wǎng)的損害極大。

1.3.2 潮流計算不收斂

潮流計算在數(shù)學(xué)上可以歸結(jié)為求解非線性方程組，而求解方法離不開迭代，如果在有限次內(nèi)迭代不能得到結(jié)果，則證明系統(tǒng)潮流計算不收斂。

1.3.3 達(dá)到最大搜索深度

由北美電力可靠協(xié)會(North American electric reliability council，NERC)提供的北美電網(wǎng)連鎖故障統(tǒng)計數(shù)據(jù)可以得出停運支路數(shù)目與發(fā)生概率之間的關(guān)系為：隨著連鎖故障停運支路的增長，對應(yīng)的發(fā)生概率迅速減小，再繼續(xù)下去故障預(yù)測將失去工程意義[17]。若以概率小于10-3為準(zhǔn)，則可設(shè)定故障等級達(dá)到6級別為最深搜索深度。因此，在預(yù)測過程中，當(dāng)搜索到第6級故障時，搜索結(jié)束。

2 聚類算法的分析研究

連鎖故障的前后關(guān)聯(lián)性是有別于其他故障的重要依據(jù)。因此在每一次關(guān)聯(lián)性指標(biāo)計算完成后，依據(jù)關(guān)聯(lián)性指標(biāo)的大小進(jìn)行聚類分類，將關(guān)聯(lián)性指標(biāo)最高的一類線路作為下級故障線路集[14]。

聚類算法的選擇直接影響著預(yù)測結(jié)果。文獻(xiàn)[13]用C均值聚類對上下級支路間關(guān)聯(lián)性指標(biāo)進(jìn)行聚類分類，但該算法存在明顯的缺陷：①聚類結(jié)果對于初始聚類中心的依賴度過高，而初始聚類中心的選擇又有不確定性或盲目性，導(dǎo)致聚類結(jié)果不準(zhǔn)確。②聚類數(shù)目必須是提前已知的，這樣才能保證聚類的準(zhǔn)確性，但是在實際應(yīng)用中，線路故障數(shù)目無法確定，因此無法提前得到聚類數(shù)目，所以該方法不能保證準(zhǔn)確率。文獻(xiàn)[18]的預(yù)測模型中引入迭代自組織數(shù)據(jù)分析聚類算法(iterative self-organizing data analysis algorithm，ISODATA)，能取得較高的預(yù)測精度，但該算法計算參數(shù)多，計算過程復(fù)雜，大量參數(shù)的調(diào)整優(yōu)化使得訓(xùn)練模型時間變長，預(yù)測速度慢。本文方法將分別把這兩種聚類算法引入模型中，與改進(jìn)SOM聚類算法進(jìn)行對比分析，來確定適合該模型的聚類算法。

SOM聚類算法是一種無監(jiān)督學(xué)習(xí)的聚類算法[19]，SOM算法的核心是通過競爭學(xué)習(xí)來實現(xiàn)其自組織功能[20]。在SOM聚類算法的網(wǎng)絡(luò)訓(xùn)練參數(shù)設(shè)定中，根據(jù)樣本數(shù)，設(shè)定輸出層神經(jīng)元為6。根據(jù)故障輸出特點，設(shè)定輸出層節(jié)點為2，兩個輸出節(jié)點分別代表不易發(fā)生故障和易發(fā)生故障。隨機(jī)初始化時，從訓(xùn)練集中隨機(jī)抽取f個輸入樣本作為初始值，本文方法抽取1 048 576個樣本作為初始值。拓?fù)漕I(lǐng)域的設(shè)計如下，其中SOM的鄰域函數(shù)為

N(dj)=e-dj/[2r2(t)](15)

式(15)中：N(dj)為鄰域函數(shù)；dj為優(yōu)勝鄰域內(nèi)節(jié)點與獲勝節(jié)點的歐氏距離；r(t)為t次迭代時鄰域范圍的半徑。

權(quán)向量調(diào)整為

Wij(t+1)=Wij(t)+ΔWij=Wij(t)+η(t)N(dj)(X-Wij) (16)

式(16)中：X為輸入向量；Wij為權(quán)向量；ΔWij為權(quán)向量變化量；η(t)為學(xué)習(xí)率。

設(shè)計學(xué)習(xí)率：訓(xùn)練時采用“競爭學(xué)習(xí)”的方式，每個輸入樣本在輸出層找到一個和它最匹配的節(jié)點，稱為激活節(jié)點；然后用隨機(jī)梯度下降法更新激活節(jié)點的參數(shù)；同時，臨近激活節(jié)點的點也根據(jù)它們距離激活節(jié)點的遠(yuǎn)近而適當(dāng)?shù)馗聟?shù)。而學(xué)習(xí)率的確定也直接影響聚類結(jié)果中分類的準(zhǔn)確度，從而影響預(yù)測結(jié)果。因此，在參考其他論文學(xué)習(xí)率確定的基礎(chǔ)上，考慮到隨著樣本的變化，靜態(tài)學(xué)習(xí)率會導(dǎo)致其準(zhǔn)確度降低的問題，因此使用動態(tài)學(xué)習(xí)率；通過實驗，采用自適應(yīng)的方式，確定學(xué)習(xí)率的大小，學(xué)習(xí)率與迭代次數(shù)成反比，隨著實驗的進(jìn)行，學(xué)習(xí)率不斷更新，能更好地適應(yīng)實驗。其定義為

η(t+1)=η(0)e-t/T(17)

當(dāng)前樣本數(shù)據(jù)學(xué)習(xí)過程結(jié)束后，輸入層傳入下一組樣本數(shù)據(jù)，迭代上述學(xué)習(xí)過程，直至t=T；或者定義最小的學(xué)習(xí)率minη，當(dāng)η(t)

改進(jìn)SOM聚類算法應(yīng)用于上下級支路間的關(guān)聯(lián)性指標(biāo)聚類分類，將上下級支路關(guān)聯(lián)性指標(biāo)作為樣本輸入，通過該方法，將關(guān)聯(lián)性最高的一類指標(biāo)找出來，并對應(yīng)到相應(yīng)的下級支路中，得到潛在發(fā)生故障的下級支路。

3 連鎖故障預(yù)測模型算法流程

電力系統(tǒng)連鎖故障預(yù)測流程圖，如圖1所示。連鎖故障路徑樹的生成算法具體步驟如下。

圖1 電力系統(tǒng)連鎖故障預(yù)測流程圖Fig.1 Flow chart of power system chain fault prediction

(1)將初始故障線路作為故障樹的根節(jié)點，斷開其中一條線路，判斷系統(tǒng)所處狀態(tài)：①系統(tǒng)穩(wěn)定，通過上下級支路間關(guān)聯(lián)性指標(biāo)計算該支路與其他支路間的關(guān)聯(lián)性，并通過改進(jìn)SOM聚類算法對指標(biāo)進(jìn)行聚類分類，將聚類結(jié)果較高的一類支路作為下級故障支路集，記錄下級支路的個數(shù)為m，并計算每條支路發(fā)生故障的概率。②系統(tǒng)解列，該條事故連鎖搜索完畢，記錄其路徑及發(fā)生故障的概率。

(3)以此類推，直到獲得完整地故障樹，并且根據(jù)每級各節(jié)點發(fā)生故障的概率，計算出該故障樹中各條故障路徑的發(fā)生概率。

(4)再斷開初始故障集合中其他故障線路，重復(fù)步驟(1)～步驟(3)，直到預(yù)測完所有初始故障線路造成的潛在故障為止。

4 算例分析

將IEEE-36節(jié)點系統(tǒng)作為算例驗證本文方法的準(zhǔn)確性、快速性和全面性。其系統(tǒng)接線圖如圖2所示。根據(jù)國家標(biāo)準(zhǔn)，中國的輸電等級分類為：高壓10～220 kV，超高壓330 ～ 750 kV。仿真中，所有發(fā)電機(jī)容量為10 MW，輸出電壓為6 kV，經(jīng)過升壓變壓器升為10 kV送出。系統(tǒng)中，10 kV線路有25條，長度為70～90 km，400 kV線路有9條，長度為175～200 km，750 kV線路有2條，長度為 2 000 km。

bus(power system busbar)表示電力系統(tǒng)母線

由1.1節(jié)實驗結(jié)果可知，初始故障評價指標(biāo)的閾值為：η1=0.2，η2=0.7，η3=0.1。由實驗可得各條母線初始故障概率如圖3所示。

圖3 初始故障集合概率分布Fig.3 Probability distribution of initial failure set

由圖3可知bus9、bus16、bus19、bus21、bus22都容易發(fā)生初始故障。在驗證本文所提方法時，選擇了初始故障概率較高的母線之間的線路作為初始故障來預(yù)測下級故障的發(fā)生，其中線路用“L母線標(biāo)號”表示。在初始故障支路的具體選擇中，不僅考慮母線故障概率，而且參考圖2的系統(tǒng)接線圖，考慮各支路在電力系統(tǒng)中是否處于重要位置，選擇的部分支路為：L19-21，L19-16，L22-9，L16-21。由于文章篇幅有限，只展示bus16和bus21之間的線路L16-21作為初始故障線路所引起的后續(xù)故障傳播情況，在該線路因故停運后，首先驗證系統(tǒng)穩(wěn)定性，然后計算該條線路與剩余線路的關(guān)聯(lián)性，通過優(yōu)化SOM聚類算法得到故障指標(biāo)最高的一類線路并將其作為下級故障線路集。

在得到第二級故障線路后，分析系統(tǒng)所處狀態(tài)：①系統(tǒng)解列，則該條故障路徑搜索結(jié)束。②系統(tǒng)穩(wěn)定，則繼續(xù)重復(fù)上面步驟，在后續(xù)搜索中遵守故障搜索停止的3個條件；此外，遇到已經(jīng)故障的線路，不搜索；在三相變壓器中，本著搜索出最大故障樹的原則，兩條二次側(cè)所接線路都搜索。得到的預(yù)測故障樹如圖4所示。

圖4 連鎖故障路徑樹Fig.4 Cascading failure path tree

以L16-21線路故障作為根節(jié)點的連鎖故障具體預(yù)測路徑及各條路徑發(fā)生故障的概率如表1所示。在確定下級故障數(shù)目時，經(jīng)過多次實驗，確定上下級支路關(guān)聯(lián)性指標(biāo)的閾值為78%，只有兩級關(guān)聯(lián)度達(dá)到78%，才能確認(rèn)為下級故障支路。這可以避免出現(xiàn)文獻(xiàn)[13]所提方法的缺陷，即只選擇一條關(guān)聯(lián)支路作為下級故障時，不符合連鎖故障發(fā)展路徑的不確定性、復(fù)雜性。也可以避免文獻(xiàn)[21]所提方法的缺陷：每條線路必須選擇4條下級故障支路而造成事故鏈規(guī)模龐大，事故鏈模型計算復(fù)雜度提高，預(yù)測時間增加等。

表1 線路L16-21斷開后的預(yù)測路徑表

如圖5所示為不同方法預(yù)測中故障等級增加對負(fù)荷損失率的影響。系統(tǒng)負(fù)荷損失的計算方法為：故障預(yù)測所涉及的線路負(fù)荷與系統(tǒng)總負(fù)荷的比值，該值越大，表示故障預(yù)測所涉及的線路越多，故障等級越深入，故障預(yù)測越趨近于結(jié)束。由圖5可知，

圖5 不同方法預(yù)測中故障等級增加對負(fù)荷損失率的影響Fig.5 The influence of fault grade increase to load loss rate in different forecasting methods

若下級故障只選擇一條線路，最終6級故障后，負(fù)荷損失率只有34%，以L21-16為初始故障，只得到一條傳播路徑，為(L21-16)-(L21-22)-(L22-9)-(L9-51/11)-(L11-25)-(L25-26)；從負(fù)荷損失、故障路徑等方面，都不符合連鎖故障的特點。若下級故障必須選擇4條線路，在4級故障時就達(dá)到了100%的負(fù)荷損失，但由于模型越來越龐大，計算時間長(具體時間如圖6所示)，不利于在線故障預(yù)測。而本文方法在6級故障時達(dá)到97%的負(fù)荷損失；如前文所述，在故障達(dá)到6級后，隨著故障等級的增加，發(fā)生概率越來越小。因此，本文方法在也符合該理論。

如圖6所示為不同方法預(yù)測中時間增加對負(fù)荷損失率的影響可知，本文方法在3.245 s達(dá)到97%的負(fù)荷損失，表示此時預(yù)測已經(jīng)結(jié)束，雖然選擇一條支路作為下級故障的方法也可以在4.167 s完成預(yù)測，但是只能達(dá)到34%的負(fù)荷損失，嚴(yán)重不符合連鎖故障發(fā)展的特點。而選擇4條支路作為下級故障的方法，雖然可以達(dá)到100%的負(fù)荷損失預(yù)測，但是需要19.564 s，時間過長，而電力系統(tǒng)故障的發(fā)生時間極短，不符合實際情況。

圖6 不同方法預(yù)測中時間增加對負(fù)荷損失率的影響Fig.6 The influence of time increase to load loss rate in different forecasting methods

將C均值聚類算法、ISODATA聚類算法引入預(yù)測模型中，對3種聚類算法的預(yù)測結(jié)果進(jìn)行對比，結(jié)果如圖7所示。由圖7可知，隨著故障等級的深入，改進(jìn)SOM聚類算法可以得到37條潛在故障線路，與ISODATA聚類算法得到的33條以及與C聚類算法得到的31條相比，可以得到更全面的潛在故障線路，擴(kuò)大了故障預(yù)測的范圍。

圖7 不同聚類算法選擇中故障等級增加對潛在故障線路數(shù)目的影響Fig.7 The influence of fault level increase to the number of potential fault lines in different clustering algorithm selection

如圖8所示為不同聚類算法選擇中時間增加對負(fù)荷損失率的影響，由圖8可知，隨著仿真時間的增加，改進(jìn)SOM聚類算法可以快速地預(yù)測出故障線路，負(fù)荷損失在3.245 s達(dá)到97%，表明該算法在此時已經(jīng)完成了預(yù)測。而ISODATA聚類算法達(dá)到97%的負(fù)荷損失需要9.987 s，證明其算法參數(shù)多，運算復(fù)雜。雖然C聚類算法可以在4.879 s完成預(yù)測，但負(fù)荷損失率為92%，不能得到全面預(yù)測結(jié)果。

圖8 不同聚類算法選擇中時間增加對負(fù)荷損失率的影響Fig.8 The influence of time increase to load loss rate in different clustering algorithm selection

通過對比實驗表明，將改進(jìn)SOM聚類算法引入預(yù)測模型中，不僅能快速得到預(yù)測結(jié)果，而且兼顧了結(jié)果的全面性。能在更大范圍上預(yù)測潛在故障的發(fā)生，為相關(guān)管理人員及時處理故障提供了有力的依據(jù)。

5 結(jié)論

以電力系統(tǒng)連鎖故障發(fā)生路徑預(yù)測為研究對象，為解決傳統(tǒng)故障樹預(yù)測方法中指標(biāo)不明確、預(yù)測時間和預(yù)測全面性難以兼顧的問題，采用IEEE36節(jié)點仿真模型的數(shù)據(jù)進(jìn)行驗證。

(1)考慮到傳統(tǒng)故障樹預(yù)測模型中，線路初始指標(biāo)、上下級支路關(guān)聯(lián)性指標(biāo)不清晰的問題，綜合考慮電力系統(tǒng)動、靜態(tài)特性；暫、穩(wěn)態(tài)特性，設(shè)計評價指標(biāo)，且賦予3個初始指標(biāo)各自合適的閾值，從初始故障選擇角度提高了預(yù)測模型的準(zhǔn)確性。

(2)首次將改進(jìn)SOM聚類算法引入電力系統(tǒng)連鎖故障預(yù)測模型中，采用自適應(yīng)的方式來調(diào)整聚類過程中的學(xué)習(xí)率，提高了聚類結(jié)果的準(zhǔn)確性。對于上下級指標(biāo)的關(guān)聯(lián)性設(shè)定合適的閾值，在保證全面搜索下級關(guān)聯(lián)支路的基礎(chǔ)上，有效減少了聚類分類環(huán)節(jié)的時間，從中間環(huán)節(jié)提高了模型的預(yù)測速度。

(3)由于實驗室硬件設(shè)備有限，因此選擇IEEE36系統(tǒng)進(jìn)行仿真實驗，實驗結(jié)論適用于該類型大小的系統(tǒng)。對于實際大系統(tǒng)而言，連鎖故障發(fā)生的機(jī)理以及傳播特點與小系統(tǒng)相同，只是系統(tǒng)支路和故障等級會成倍數(shù)增加，程序運算將更加復(fù)雜，為保證快速預(yù)測更換性能更好的中央處理器(CPU)或使用圖形處理器(GPU)，將完全滿足大系統(tǒng)連鎖故障的預(yù)測要求。

(4)研究工作主要針對電力系統(tǒng)連鎖故障的預(yù)測開展，而如何對預(yù)測結(jié)果進(jìn)行風(fēng)險評估，建立可視化界面，更加清晰明確地反映電力系統(tǒng)連鎖故障的發(fā)展情況需要進(jìn)一步探究。因此，若進(jìn)行下一步研究，需要尋找合適的路徑風(fēng)險評估模型，開發(fā)建立可視化界面，從而能有效及時地預(yù)測電力系統(tǒng)連鎖故障的發(fā)生。