本發(fā)明涉及電力網(wǎng)絡(luò)分塊層面，具體涉及一種基于快速分裂算法的節(jié)點(diǎn)遷移網(wǎng)絡(luò)分塊優(yōu)化方法。

背景技術(shù)：

電力系統(tǒng)并行計(jì)算的基本思想是將電力網(wǎng)絡(luò)按照地理或網(wǎng)絡(luò)拓?fù)溥B接結(jié)構(gòu)，分割成若干個(gè)子網(wǎng)絡(luò)和一個(gè)協(xié)調(diào)網(wǎng)絡(luò)，而子網(wǎng)絡(luò)彼此之間通過一組與協(xié)調(diào)網(wǎng)絡(luò)對應(yīng)的邊界量進(jìn)行聯(lián)系。此時(shí)，子網(wǎng)絡(luò)內(nèi)部相關(guān)運(yùn)算可實(shí)行并行化計(jì)算。其中，作為分割方法之一的節(jié)點(diǎn)分裂法在具有將電力系統(tǒng)網(wǎng)絡(luò)劃分，以實(shí)現(xiàn)并行計(jì)算優(yōu)勢的同時(shí)，也會(huì)因無法客觀合理地確定分裂節(jié)點(diǎn)集，而出現(xiàn)無法擴(kuò)展應(yīng)用的瓶頸。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明所要解決的技術(shù)問題是提供一種基于快速分裂算法的節(jié)點(diǎn)遷移網(wǎng)絡(luò)分塊優(yōu)化方法，其應(yīng)用于電力系統(tǒng)潮流的并行計(jì)算時(shí)，將變得更加完整、系統(tǒng)，即為節(jié)點(diǎn)分裂算法的割裂節(jié)點(diǎn)成員的確定提供一種客觀依據(jù)。

為解決上述技術(shù)問題，本發(fā)明采用的技術(shù)方案是：

一種基于快速分裂算法的節(jié)點(diǎn)遷移網(wǎng)絡(luò)分塊優(yōu)化方法，包括以下步驟：

步驟1：基于快速分裂算法擇取薄弱線路，即：計(jì)算電力網(wǎng)絡(luò)中所有線路的邊聚類系數(shù)；按邊聚類系數(shù)，對網(wǎng)絡(luò)當(dāng)前包含所有線路進(jìn)行升序排序，獲取當(dāng)前邊聚類系數(shù)最小的線路，記作curLine后并添加至薄弱線路集合中；從當(dāng)前網(wǎng)絡(luò)中移除curLine，重復(fù)排序添加步驟，直至電力網(wǎng)絡(luò)所有節(jié)點(diǎn)均退化為一個(gè)獨(dú)立節(jié)點(diǎn)社團(tuán)；

步驟2：將切割之路所連接節(jié)點(diǎn)，按照一定原則，擇取當(dāng)前切割支路將其從所對應(yīng)的子網(wǎng)絡(luò)內(nèi)部遷移到協(xié)調(diào)網(wǎng)絡(luò)內(nèi)部，即撕裂節(jié)點(diǎn)集合；

所述原則是指：在逐步將當(dāng)前薄弱線路所對應(yīng)兩端節(jié)點(diǎn)分別作為撕裂節(jié)點(diǎn)時(shí)：1)優(yōu)先遷移鄰域連接度大的節(jié)點(diǎn)；2)優(yōu)先遷移不增加以孤島形式存在節(jié)點(diǎn)的撕裂節(jié)點(diǎn)；3)若兩端節(jié)點(diǎn)鄰域連接度相同，且均不增加孤島節(jié)點(diǎn)，隨機(jī)選取節(jié)點(diǎn)作為撕裂節(jié)點(diǎn)；4)此前已作為撕裂節(jié)點(diǎn)的節(jié)點(diǎn)，均執(zhí)行鎖定操作，且后續(xù)新增撕裂節(jié)點(diǎn)若被鎖定在撕裂節(jié)點(diǎn)集合中，則轉(zhuǎn)至下一薄弱線路的撕裂節(jié)點(diǎn)的擇取；5)每獲取一個(gè)撕裂節(jié)點(diǎn)時(shí)，需采用深度優(yōu)先搜索算法判斷當(dāng)前網(wǎng)絡(luò)分割格局，存儲(chǔ)當(dāng)前各子網(wǎng)絡(luò)和協(xié)調(diào)網(wǎng)絡(luò)所包含節(jié)點(diǎn)編號(hào)和線路信息；

步驟3：定義分割目標(biāo)條件函數(shù)K，表征分割的合理性：

式中，符號(hào)“|”表征條件前提，該符號(hào)右邊為前提條件，左邊對應(yīng)最終目標(biāo)條件函數(shù)值；V_conet為協(xié)調(diào)網(wǎng)絡(luò)所包含的撕裂節(jié)點(diǎn)數(shù)目；V_{net_i}、V_{net_j}為劃分后的第i、j個(gè)子網(wǎng)絡(luò)內(nèi)部所包含節(jié)點(diǎn)的總數(shù)目。

進(jìn)一步的，在撕裂過程中，設(shè)定對應(yīng)子網(wǎng)絡(luò)集合中各子網(wǎng)絡(luò)規(guī)模占據(jù)完整網(wǎng)絡(luò)規(guī)模的閾值thres，若低于此閾值的子網(wǎng)絡(luò)數(shù)目達(dá)到設(shè)定的比例系數(shù)prop，則終止撕裂過程。

與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果是：其應(yīng)用于電力系統(tǒng)潮流、短路電流等問題的并行計(jì)算應(yīng)用時(shí)，將變得更加完整、系統(tǒng)，即為節(jié)點(diǎn)分裂算法的割裂節(jié)點(diǎn)成員的確定提供一種客觀依據(jù)。

附圖說明

圖1為支路切割-節(jié)點(diǎn)撕裂轉(zhuǎn)換示意圖。

圖2為節(jié)點(diǎn)遷移網(wǎng)絡(luò)分塊優(yōu)化流程圖。

圖3系統(tǒng)分區(qū)格局示意圖。

具體實(shí)施方式

下面結(jié)合附圖和具體實(shí)施方式對本發(fā)明作進(jìn)一步詳細(xì)的說明。本發(fā)明方法首先以快速分裂算法確定電力網(wǎng)絡(luò)拓?fù)滏溄雨P(guān)系薄弱的線路集合；其次，引入“遷移”的概念，將該集合中薄弱線路轉(zhuǎn)換為分裂節(jié)點(diǎn)，以確定分裂節(jié)點(diǎn)集合和分割格局集合；最后，依據(jù)分割目標(biāo)條件函數(shù)，擇取電力網(wǎng)絡(luò)最優(yōu)分割格局。詳述如下：

一、基于快速分裂算法的薄弱線路擇取

復(fù)雜網(wǎng)絡(luò)理論在電力系統(tǒng)中的應(yīng)用已日趨成熟，如電網(wǎng)結(jié)構(gòu)辨識(shí)、無功分區(qū)、網(wǎng)絡(luò)動(dòng)力學(xué)分析和500/220kV電磁環(huán)網(wǎng)解環(huán)等。復(fù)雜網(wǎng)絡(luò)理論中的快速分裂算法作為社團(tuán)發(fā)現(xiàn)方面的開創(chuàng)性算法，是一種典型的分裂算法。其基本思想采用一種新的標(biāo)準(zhǔn)——邊聚類系數(shù)取代GN分裂算法中的邊介數(shù)來衡量連接兩個(gè)不同社團(tuán)的邊。仿照節(jié)點(diǎn)聚類系數(shù)的定義，邊聚類系數(shù)定義為實(shí)際包含該邊的三角形的數(shù)目與所有可能包含該邊的三角形的數(shù)目之比。具體的，對于連接節(jié)點(diǎn)i和節(jié)點(diǎn)j的邊，它的聚類系數(shù)定義為

C_ij.3＝(Z_ij.3+1)/(min[k_i-1,k_j-1]) (1)

式中：C_ij.3表示包含該邊的三角形的數(shù)目；k_i、k_j則對應(yīng)節(jié)點(diǎn)i、j的節(jié)點(diǎn)度；其中，min[k_i-1,k_j-1]則表示兩者的較小值。

由式(1)可知，社團(tuán)內(nèi)部邊的聯(lián)系比較緊密，因此存在很多三角形，即它的邊聚類系數(shù)值C_ij.3較大，而連接兩個(gè)社團(tuán)間的邊包含在極少的三角形內(nèi)或者不被任何三角形所包含，因此C_ij.3值較小。由此可見，C_ij.3大小即可作為一種衡量邊的社團(tuán)間連接性的標(biāo)準(zhǔn)。該值所對應(yīng)邊即為所述薄弱線路。

因此，本發(fā)明結(jié)合快速分裂算法以獲取表征電力系統(tǒng)網(wǎng)絡(luò)拓?fù)渎?lián)系脆弱的薄弱線路集合，包括步驟：1)計(jì)算電力網(wǎng)絡(luò)中所有線路的邊聚類系數(shù)；2)按邊聚類系數(shù)，對網(wǎng)絡(luò)當(dāng)前包含所有線路進(jìn)行降序排序，獲取當(dāng)前邊聚類系數(shù)最小的線路，記作curLine后并添加至薄弱線路集合中；3)從當(dāng)前網(wǎng)絡(luò)中移除curLine，重復(fù)步驟2)，直至電力網(wǎng)絡(luò)所有節(jié)點(diǎn)均退化為一個(gè)獨(dú)立節(jié)點(diǎn)社團(tuán)。

二、節(jié)點(diǎn)撕裂分割轉(zhuǎn)換

依據(jù)快速分裂算法所獲取的薄弱線路集合后，將其逐步視作電力網(wǎng)絡(luò)初始支路切割對象。為得到撕裂節(jié)點(diǎn)形式下的分割，則需要將切割支路轉(zhuǎn)換為撕裂節(jié)點(diǎn)。轉(zhuǎn)換基本思想為將切割之路所連接節(jié)點(diǎn)，按照一定的原則，擇取其一從所對應(yīng)的子網(wǎng)絡(luò)內(nèi)部遷移到協(xié)調(diào)網(wǎng)絡(luò)內(nèi)部，即撕裂節(jié)點(diǎn)集合。

圖1表示一個(gè)簡單實(shí)例。圖1(a)是兩條支路切割分割的網(wǎng)絡(luò)，確定排序在前兩位的切割支路為(v2,v4)和(v5,v4)，邊聚類系數(shù)分別為0.67和0.50。要將其轉(zhuǎn)換為節(jié)點(diǎn)撕裂分割需要的切割對象，有多種轉(zhuǎn)換方法，圖1顯示了可能的2種轉(zhuǎn)換，圖1(b)將子網(wǎng)絡(luò)1中的節(jié)點(diǎn)v2、v5遷移至撕裂節(jié)點(diǎn)集合中；類似地，圖1(c)將子網(wǎng)絡(luò)2中的節(jié)點(diǎn)v4遷移至撕裂節(jié)點(diǎn)集合中。兩者同為節(jié)點(diǎn)撕裂分割，均將完整網(wǎng)絡(luò)劃分成子網(wǎng)絡(luò)1、2和協(xié)調(diào)網(wǎng)絡(luò)形式。其中，撕裂節(jié)點(diǎn)集合分別為{v4}、{v2、v5}。

不同的轉(zhuǎn)換方法所形成的協(xié)調(diào)網(wǎng)絡(luò)和子網(wǎng)絡(luò)，所包含的節(jié)點(diǎn)成員各不相同。其中，按{v4}撕裂后的子網(wǎng)絡(luò)1、2和協(xié)調(diào)網(wǎng)絡(luò)所包含節(jié)點(diǎn)成員分別為{v1、v2、v3、v5}、{v6、v7、v8}和{v4}；按{v2、v5}撕裂后的子網(wǎng)絡(luò)1、2和協(xié)調(diào)網(wǎng)絡(luò)所包含節(jié)點(diǎn)成員分別為{v1、v3}、{v4、v6、v7、v8}和{v2、v5}。兩者撕裂的區(qū)別在于前者僅需遷移v4節(jié)點(diǎn)，而后者需遷移兩個(gè)節(jié)點(diǎn)，才實(shí)現(xiàn)網(wǎng)絡(luò)分塊的目的，顯然增加了操作的復(fù)雜性。

為了得到最佳的遷移，必須在遷移過程中計(jì)及一定的原則，以此優(yōu)化遷移組合，才能使分割目標(biāo)函數(shù)最優(yōu)。進(jìn)一步考慮五個(gè)原則，具體指在逐步將當(dāng)前薄弱線路所對應(yīng)兩端節(jié)點(diǎn)分別作為撕裂節(jié)點(diǎn)時(shí)：1)優(yōu)先遷移鄰域連接度大的節(jié)點(diǎn)。鄰域連接度判斷的前提為薄弱線路k、k+1條的逐步遍歷；2)優(yōu)先遷移不增加以孤島形式存在節(jié)點(diǎn)的撕裂節(jié)點(diǎn)；3)若兩端節(jié)點(diǎn)鄰域連接度相同，且均不增加孤島節(jié)點(diǎn)，隨機(jī)選取節(jié)點(diǎn)作為撕裂節(jié)點(diǎn)；4)此前已作為撕裂節(jié)點(diǎn)的節(jié)點(diǎn)，均執(zhí)行鎖定操作，且后續(xù)新增撕裂節(jié)點(diǎn)若被鎖定在撕裂節(jié)點(diǎn)集合中，則轉(zhuǎn)至下一薄弱線路的撕裂節(jié)點(diǎn)的擇??；5)每獲取一個(gè)撕裂節(jié)點(diǎn)時(shí)，需采用深度優(yōu)先搜索算法判斷當(dāng)前網(wǎng)絡(luò)分割格局，存儲(chǔ)當(dāng)前各子網(wǎng)絡(luò)和協(xié)調(diào)網(wǎng)絡(luò)所包含節(jié)點(diǎn)編號(hào)和線路信息。

在撕裂過程中，需設(shè)定對應(yīng)子網(wǎng)絡(luò)集合中各子網(wǎng)絡(luò)規(guī)模占據(jù)完整網(wǎng)絡(luò)規(guī)模的閾值thres，若低于此閾值的子網(wǎng)絡(luò)數(shù)目達(dá)到設(shè)定的比例系數(shù)prop，則終止撕裂過程。閾值設(shè)定依據(jù)網(wǎng)絡(luò)規(guī)模而定，且此閾值的設(shè)定對網(wǎng)絡(luò)劃分影響較大，即劃分后可能會(huì)出現(xiàn)較多小規(guī)模的子網(wǎng)絡(luò)或獨(dú)立節(jié)點(diǎn)。而增加過多過小規(guī)模網(wǎng)絡(luò)的數(shù)目，將因各子網(wǎng)絡(luò)數(shù)據(jù)的處理而顯得繁瑣，且浪費(fèi)進(jìn)程資源。此時(shí)，可采用合并的方式適當(dāng)聚集各個(gè)小規(guī)模的子網(wǎng)絡(luò)，以形成中等規(guī)模的網(wǎng)絡(luò)，合并后新增子網(wǎng)絡(luò)規(guī)模需達(dá)至此設(shè)定閾值，方可作為最終子網(wǎng)絡(luò)。

對于原則2)，系統(tǒng)網(wǎng)絡(luò)劃分目的為實(shí)現(xiàn)潮流并行計(jì)算，以提高運(yùn)算效率。新增孤島節(jié)點(diǎn)作為一個(gè)子網(wǎng)絡(luò)進(jìn)行計(jì)算時(shí)，如短路電流計(jì)算時(shí)該網(wǎng)絡(luò)各節(jié)點(diǎn)相對于短路節(jié)點(diǎn)的自阻抗或互阻抗元素的計(jì)算量顯然很小，且相比孤島節(jié)點(diǎn)，其它子網(wǎng)絡(luò)的規(guī)模仍會(huì)很大，計(jì)算速度仍取決于大規(guī)模子網(wǎng)絡(luò)的自阻抗、互阻抗相關(guān)變量的運(yùn)算。因此，孤立節(jié)點(diǎn)過多既浪費(fèi)進(jìn)程資源，也無法提高并行運(yùn)算的效率。

三、分割優(yōu)化目標(biāo)條件函數(shù)

分割優(yōu)化目標(biāo)是使與潮流計(jì)算對應(yīng)的各個(gè)子進(jìn)程負(fù)載盡量均衡，即各子網(wǎng)絡(luò)內(nèi)部并行計(jì)算量大致平衡，如此各個(gè)子任務(wù)執(zhí)行時(shí)間將不依賴于規(guī)模過大子網(wǎng)絡(luò)的自阻抗、互阻抗相關(guān)變量的計(jì)算；且達(dá)到網(wǎng)絡(luò)客觀劃分的前提下，協(xié)調(diào)網(wǎng)絡(luò)所包含的撕裂節(jié)點(diǎn)盡量少。因此，可定義如式(1)所示的分割目標(biāo)條件函數(shù)K，以表征分割的合理性。式(1)中，符號(hào)“|”表征條件前提，該符號(hào)右邊為前提條件，左邊對應(yīng)最終目標(biāo)條件函數(shù)值。

式中：V_conet為協(xié)調(diào)網(wǎng)絡(luò)所包含的撕裂節(jié)點(diǎn)數(shù)目；V_{net_i}、V_{net_j}為劃分后的第i、j個(gè)子網(wǎng)絡(luò)內(nèi)部所包含節(jié)點(diǎn)的總數(shù)目。

為驗(yàn)證所提節(jié)點(diǎn)遷移網(wǎng)絡(luò)分塊方法的有效性，以IEEE 30節(jié)點(diǎn)系統(tǒng)作為測試算例。采用快速分裂算法獲取該網(wǎng)絡(luò)的薄弱線路集合，如表1所示。

表1 IEEE 30節(jié)點(diǎn)系統(tǒng)部分線路邊聚類系數(shù)

鑒于線路較多，表1僅列寫出對排序在前8條線路薄弱線路的邊聚類系數(shù)，表中排序順序自上而下，自左向右。在此基礎(chǔ)上，在設(shè)定thres、prop分別為0.3、0.95后，將該系統(tǒng)逐步轉(zhuǎn)換成子網(wǎng)絡(luò)集合和協(xié)調(diào)網(wǎng)絡(luò)并存的形式，進(jìn)一步結(jié)合分割條件目標(biāo)函數(shù)K，確定最優(yōu)分割格局，即8個(gè)子網(wǎng)絡(luò)和一個(gè)協(xié)調(diào)網(wǎng)絡(luò)，對應(yīng)K值最大，為0.0909。最后，將小規(guī)模網(wǎng)絡(luò)進(jìn)行適當(dāng)合并，最終系統(tǒng)分區(qū)格局由如圖3中表征，即子網(wǎng)絡(luò)1、2，分別涵蓋節(jié)點(diǎn)成員為{7、10、17、20、21、16、19、18、14}、{23、24、25、26、28、8、29、30、11、13}；協(xié)調(diào)網(wǎng)絡(luò)包含節(jié)點(diǎn)為{3、4、6、9、12、15、22、27、2、5、1}。