焦慧明 ，劉天琪 ，蘇學(xué)能 ，王 彪 ，白懿鵬

（1.四川大學(xué) 電氣信息學(xué)院，四川 成都 610065；2.國(guó)網(wǎng)四川省電力公司電力科學(xué)研究院，四川 成都 610072）

0 引言

電力系統(tǒng)安全穩(wěn)定導(dǎo)則［1］提出：隨著高一級(jí)電壓電網(wǎng)的建設(shè)，下級(jí)電壓電網(wǎng)應(yīng)逐步實(shí)現(xiàn)分區(qū)運(yùn)行，相鄰分區(qū)之間保持互為備用，以避免和消除影響電網(wǎng)安全穩(wěn)定的不同電壓等級(jí)的電磁環(huán)網(wǎng)。目前已有很多關(guān)于電磁環(huán)網(wǎng)分區(qū)運(yùn)行的研究，但大多依據(jù)工程經(jīng)驗(yàn)來(lái)制定分區(qū)方案［2-5］，相對(duì)主觀，且方案制定周期較長(zhǎng)。

復(fù)雜網(wǎng)絡(luò)理論在電力系統(tǒng)中的應(yīng)用日趨成熟［6］，已有學(xué)者從復(fù)雜網(wǎng)絡(luò)的視角研究電網(wǎng)分區(qū)方案的制定［7］。復(fù)雜網(wǎng)絡(luò)中的 GN（Girvan-Newman）算法［8］作為社區(qū)發(fā)現(xiàn)方面的開(kāi)創(chuàng)性算法，是一種典型的分裂算法。文獻(xiàn)［9］在結(jié)合GN算法的基礎(chǔ)上，提出了一種基于社團(tuán)結(jié)構(gòu)特性量化分析的電磁環(huán)網(wǎng)分區(qū)方法。該方法依據(jù)定量表征線(xiàn)路薄弱程度的綜合邊介數(shù)指標(biāo)，對(duì)電磁環(huán)網(wǎng)予以解環(huán)。然而，該方法在解環(huán)過(guò)程中忽略了電磁環(huán)網(wǎng)分區(qū)需計(jì)及的相關(guān)原則，如具有2臺(tái)主變的樞紐變電站需“手拉手”運(yùn)行和解環(huán)后各分區(qū)不能存在3座及以上站點(diǎn)的串供結(jié)構(gòu)等；此外，由于GN算法是不斷移除當(dāng)前綜合邊介數(shù)最大的邊，直接運(yùn)用會(huì)導(dǎo)致指定分區(qū)數(shù)目時(shí)只有一種分區(qū)方案，而此時(shí)該方案還會(huì)受到分區(qū)原則的制約，從而有可能導(dǎo)致無(wú)法獲得可行分區(qū)方案。針對(duì)分區(qū)方案的評(píng)估方法已有模糊評(píng)價(jià)模型［10］、多層次模糊評(píng)價(jià)模型［11-12］以及層次分析法［13］等。此外，主成分分析PCA（Principal Component Analysis）法［14］利用降維的思想，降低指標(biāo)集維數(shù)，同時(shí)保留指標(biāo)集中對(duì)方差貢獻(xiàn)較大的指標(biāo)，可有效確定各指標(biāo)對(duì)方案的貢獻(xiàn)，客觀評(píng)價(jià)方案的優(yōu)劣。

鑒于上述分析，在文獻(xiàn)［9］的基礎(chǔ)上，提出一種基于回溯思想改進(jìn)GN算法的電磁環(huán)網(wǎng)分區(qū)方法。首先，該方法分別采用深度優(yōu)先搜索DFS（Depth-First Search）算法［15］和圖論節(jié)點(diǎn)度檢測(cè)方式，實(shí)現(xiàn)樞紐變電站運(yùn)行方式辨識(shí)和3臺(tái)及以上串供結(jié)構(gòu)的檢測(cè)及修正；進(jìn)一步，為充分發(fā)揮GN算法的優(yōu)勢(shì)以擴(kuò)展分區(qū)方案的多樣性，采用回溯思想改進(jìn)GN算法的分區(qū)過(guò)程。其次，為提高分區(qū)方案的校驗(yàn)效率，提出結(jié)合進(jìn)程連接PSD-BPA的方式以實(shí)現(xiàn)分區(qū)方案的靜態(tài)安全分析和短路電流自動(dòng)校驗(yàn)，從而獲取校驗(yàn)成功方案的校驗(yàn)指標(biāo)；并將此指標(biāo)與衡量分區(qū)劃分結(jié)構(gòu)合理性的模塊度相結(jié)合，以作為評(píng)估指標(biāo)集；進(jìn)而，用PCA法以該指標(biāo)集衡量分區(qū)方案的優(yōu)劣。最后，將所提方法應(yīng)用于某實(shí)際地區(qū)500 kV/220 kV電磁環(huán)網(wǎng)的解環(huán)，結(jié)果表明，該方法可有效實(shí)現(xiàn)電磁環(huán)網(wǎng)的自動(dòng)分區(qū)，并確保所得方案滿(mǎn)足實(shí)際電網(wǎng)分區(qū)原則。

1 GN算法用于解環(huán)分區(qū)的缺陷及解決措施

1.1 GN算法用于解環(huán)分區(qū)的缺陷

GN算法屬于分裂層次聚類(lèi)算法，通過(guò)移除綜合邊介數(shù)最大的邊來(lái)劃分隱藏在網(wǎng)絡(luò)中的不同社團(tuán)，并將模塊度函數(shù)［16］作為網(wǎng)絡(luò)劃分的衡量標(biāo)準(zhǔn)。但在電磁環(huán)網(wǎng)分區(qū)的應(yīng)用中，沒(méi)有考慮電磁環(huán)網(wǎng)分區(qū)時(shí)應(yīng)滿(mǎn)足的原則。

以500 kV/220 kV電磁環(huán)網(wǎng)為例，220 kV側(cè)環(huán)網(wǎng)分區(qū)的方案應(yīng)綜合考慮電網(wǎng)結(jié)構(gòu)合理、調(diào)度運(yùn)行靈活和投資規(guī)模適當(dāng)?shù)纫蛩兀猸h(huán)后的供電區(qū)需滿(mǎn)足2個(gè)原則：原則1，具有2臺(tái)主變的樞紐變電站“手拉手”運(yùn)行，其中樞紐變電站指500 kV變電站；原則2，為保障分區(qū)的供電可靠性，規(guī)定不超過(guò)3座及以上的220 kV站點(diǎn)串供運(yùn)行。

因此，為使得到的分區(qū)方案具有實(shí)際意義，需對(duì)GN算法進(jìn)行改造完善。

1.2 解決措施

1.2.1 采用DFS辨識(shí)樞紐變電站的運(yùn)行方式

針對(duì)上述原則1，提出采用DFS辨識(shí)樞紐變電站的運(yùn)行方式。通過(guò)DFS來(lái)獲取當(dāng)前方案的分區(qū)情況，進(jìn)而辨識(shí)500 kV變電站的運(yùn)行方式。主要步驟如下：

a.獲取表征電網(wǎng)220 kV站點(diǎn)的鄰接關(guān)系矩陣和斷線(xiàn)割集；

b.采用DFS判斷該方案是否形成新的分區(qū)，若是，則獲取其分區(qū)數(shù)目及其包含的站點(diǎn)；

c.遍歷方案下的所有分區(qū)，獲取各分區(qū)中與樞紐變電站直接相連的站點(diǎn)；

d.由步驟c獲取各分區(qū)下的樞紐變電站數(shù)目，若數(shù)目小于2，則進(jìn)一步判斷是否為2臺(tái)主變的樞紐變電站，若是，則不滿(mǎn)足要求。

1.2.2 串供結(jié)構(gòu)的檢測(cè)

針對(duì)上述原則2，提出基于圖論節(jié)點(diǎn)度概念檢測(cè)串供結(jié)構(gòu)的方法。

圖論中某節(jié)點(diǎn)的節(jié)點(diǎn)度定義為與該節(jié)點(diǎn)相連的節(jié)點(diǎn)數(shù)目。如圖1所示，引入圖論中節(jié)點(diǎn)度的概念，串供結(jié)構(gòu)末端站點(diǎn)的節(jié)點(diǎn)度為1，記為A站，與A站相連的記為B站。假定存在3臺(tái)所組成的串供結(jié)構(gòu)，則B站的節(jié)點(diǎn)度為2，B站除了與A站相連，還應(yīng)與另一站點(diǎn)相連，記為C站。串供結(jié)構(gòu)辨識(shí)方法主要判斷由節(jié)點(diǎn)度為1的站點(diǎn)開(kāi)始，是否連續(xù)存在2個(gè)及以上相鄰站的節(jié)點(diǎn)度為2，若是，則存在串供結(jié)構(gòu)。

圖1 串供結(jié)構(gòu)的辨識(shí)Fig.1 Identification of serial-supplying structure

1.2.3 串供結(jié)構(gòu)的修正

若存在3座及以上220 kV站點(diǎn)的串供結(jié)構(gòu)，則需對(duì)其進(jìn)行修正。此時(shí)方案中的解環(huán)斷線(xiàn)需進(jìn)行調(diào)整，以消除3座及以上220 kV站點(diǎn)的串供結(jié)構(gòu)，如圖2所示。

圖2 串供結(jié)構(gòu)的修正Fig.2 Correction of serial-supplying structure

串供結(jié)構(gòu)的修正主要步驟如下：

a.獲取串供結(jié)構(gòu)列表中節(jié)點(diǎn)度為1的A站；

b.依據(jù)分區(qū)方案斷線(xiàn)信息，得到A站對(duì)應(yīng)的斷線(xiàn)編號(hào)n；

c.遍歷串供結(jié)構(gòu)列表中的站點(diǎn)，從尾部往前推移得到線(xiàn)路編號(hào)n+1，替代斷線(xiàn)信息中的線(xiàn)路n；

d.結(jié)合串供結(jié)構(gòu)檢測(cè)方法，判斷修正后的方案是否還存在3座及以上220 kV站點(diǎn)的串供結(jié)構(gòu)，若存在，則需繼續(xù)修正。

1.3 計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法流程

通過(guò)采用DFS辨識(shí)樞紐變電站的運(yùn)行方式和3臺(tái)及以上串供結(jié)構(gòu)的檢測(cè)及修正方法，完善GN算法的執(zhí)行過(guò)程，提出計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法，以便獲得具有實(shí)際意義的分區(qū)方案。計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法流程如圖3所示。

圖3 計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法流程圖Fig.3 Flowchart of GN algorithm considering partitioning principles of electromagnetic loop network

2 基于回溯思想改進(jìn)GN算法的自動(dòng)分區(qū)方法及方案評(píng)估

2.1 基于回溯思想改進(jìn)的GN算法

通過(guò)計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法雖然保證了所得方案滿(mǎn)足分區(qū)原則，但由于GN算法的執(zhí)行過(guò)程是不斷地移除當(dāng)前綜合邊介數(shù)最大的邊，導(dǎo)致指定分區(qū)數(shù)目時(shí)只有一種分區(qū)方案，而此時(shí)該方案還會(huì)受到分區(qū)原則的制約，可能導(dǎo)致無(wú)法得到可行分區(qū)方案。

為充分發(fā)揮GN算法利用綜合邊介數(shù)信息實(shí)現(xiàn)分區(qū)的優(yōu)勢(shì)，進(jìn)一步擴(kuò)展分區(qū)方案的多樣性，提出基于回溯思想改進(jìn)GN算法的分區(qū)方法。其流程如圖4所示。該方法首先依據(jù)GN算法獲取全網(wǎng)邊介數(shù)信息；然后由工程經(jīng)驗(yàn)指定最大可斷線(xiàn)數(shù)目得到斷線(xiàn)割集；最后對(duì)斷線(xiàn)割集并行從頭部、尾部進(jìn)行回溯。其中將頭部閉合線(xiàn)路過(guò)程（圖4中①部分）和尾部閉合線(xiàn)路過(guò)程（圖4中②部分）分別放在2個(gè)不同的進(jìn)程中并行執(zhí)行，以提高運(yùn)算效率。

圖4 基于回溯思想的改進(jìn)GN算法流程圖Fig.4 Flowchart of GN algorithm improved with backtracking concept

2.2 方案評(píng)估模型

2.2.1 方案校驗(yàn)及指標(biāo)

500 kV/220 kV電磁環(huán)網(wǎng)分區(qū)的目的是為了降低因500 kV電網(wǎng)結(jié)構(gòu)的加強(qiáng)而導(dǎo)致的220 kV電網(wǎng)的單相、三相短路電流水平，以保證電網(wǎng)運(yùn)行的穩(wěn)定性和供電可靠性。改進(jìn)GN算法可得到更多的分區(qū)方案，傳統(tǒng)的方案校驗(yàn)難以自動(dòng)批量校驗(yàn)方案的各項(xiàng)指標(biāo)。為提高校驗(yàn)方案的效率，進(jìn)一步驗(yàn)證方案的合理性，從而自動(dòng)甄選出滿(mǎn)足供電可靠性和短路電流裕度的方案，提出結(jié)合進(jìn)程連接PSD-BPA的方式，以實(shí)現(xiàn)分區(qū)方案的靜態(tài)安全分析和短路電流校驗(yàn)，并從中獲取校驗(yàn)成功方案的校驗(yàn)指標(biāo)。其中，校驗(yàn)指標(biāo)分為2類(lèi)。

a.供電負(fù)荷能力。具體指通過(guò)母線(xiàn)“N-1”、主變“N-1”、線(xiàn)路“N-1”、同塔“N-2”（同塔率大于 70%）等靜態(tài)安全分析，獲取各分區(qū)的制約故障信息和實(shí)際供電負(fù)荷能力的信息，并以 K1、K2、K3、K4分別表示對(duì)應(yīng)指標(biāo)。以某分區(qū)方案的母線(xiàn)“N-1”指標(biāo)K1為例，解環(huán)后形成m個(gè)分區(qū)，則該方案的K1表示為：

其中，Pi.base（i=1，2，…，m）為各分區(qū)的基準(zhǔn)負(fù)荷；Pi1（i=1，2，…，m）為各分區(qū)母線(xiàn)“N-1”的校驗(yàn)指標(biāo)。

b.短路電流校驗(yàn)。具體指與樞紐變電站直接連接的220 kV母線(xiàn)和部分開(kāi)關(guān)遮斷容量為40 kA的220 kV變電站220 kV母線(xiàn)的單相、三相短路電流的校驗(yàn)。需說(shuō)明的是，當(dāng)樞紐變電站母線(xiàn)的短路電流不超標(biāo)時(shí)，電網(wǎng)其他母線(xiàn)的電流一般也不會(huì)超標(biāo)。以K5和K6分別表示與樞紐變電站直接相連的220 kV母線(xiàn)的單相短路電流裕度指標(biāo)及其三相短路電流裕度指標(biāo)；以K7和K8分別表示開(kāi)關(guān)遮斷容量為40 kA的220 kV變電站220 kV母線(xiàn)的單相短路電流裕度指標(biāo)及其三相短路電流裕度指標(biāo)。以樞紐變電站直接相連220 kV母線(xiàn)的單相短路電流裕度指標(biāo)為例，假定有n座樞紐變電站，則該方案K5表示為：

其中，Ii.base（i=1，2，…，n）為各站斷路器的遮斷容量閾值；Iis（i=1，2，…，n）為采取解環(huán)措施后各站的實(shí)際單相短路電流值。

2.2.2 方案評(píng)估

參照文獻(xiàn)［9］定義模塊度指標(biāo) Q［16］，以客觀合理地表征分區(qū)劃分結(jié)構(gòu)的合理性，考慮到電力網(wǎng)絡(luò)的線(xiàn)路長(zhǎng)度及參數(shù)特性，選取支路導(dǎo)納模值作為邊權(quán)重，即考慮邊權(quán)重模塊度指標(biāo)越大，說(shuō)明分區(qū)內(nèi)部的連接緊密程度越高，則電力網(wǎng)絡(luò)的社團(tuán)結(jié)構(gòu)劃分得越好。假定分區(qū)數(shù)目為l，則：

其中，e 為一個(gè) l×l階的對(duì)稱(chēng)矩陣；eii（1≤i≤l）為連接分區(qū)i內(nèi)節(jié)點(diǎn)之間的邊與初始網(wǎng)絡(luò)所有邊的比值；Tr（e）為連接同一分區(qū)內(nèi)部節(jié)點(diǎn)的邊與網(wǎng)絡(luò)所有邊之間的比值；ai為連接分區(qū)i的節(jié)點(diǎn)的邊與初始網(wǎng)絡(luò)所有邊的比值。

由8個(gè)校驗(yàn)指標(biāo)和模塊度指標(biāo)Q共同構(gòu)成評(píng)估指標(biāo)集，并在此基礎(chǔ)上采用PCA法［14］處理該指標(biāo)集建立評(píng)估模型。

首先，獲取待評(píng)估方案每一指標(biāo)的均值和方差，其中，短路電流指標(biāo)以其裕度值表示，使評(píng)估指標(biāo)集內(nèi)所有指標(biāo)均為正向指標(biāo)，并使用z-score標(biāo)準(zhǔn)化方法［17］做均一化處理，以消除量綱的影響。經(jīng)過(guò)處理的數(shù)據(jù)符合標(biāo)準(zhǔn)正態(tài)分布，其轉(zhuǎn)換函數(shù)為：

其中，μ為所有樣本數(shù)據(jù)的均值；σ為所有樣本數(shù)據(jù)的標(biāo)準(zhǔn)差。

然后，獲取其方差矩陣，包括自方差和協(xié)方差，通過(guò)對(duì)協(xié)方差矩陣進(jìn)行特征分解，以得出指標(biāo)集的主成分（即特征向量）與對(duì)應(yīng)的權(quán)值（即特征值），進(jìn)而降低指標(biāo)集的維數(shù)，同時(shí)保持指標(biāo)集中對(duì)方差貢獻(xiàn)最大的指標(biāo)。

最后，形成主成分矩陣及綜合評(píng)價(jià)指標(biāo)，根據(jù)綜合評(píng)價(jià)指標(biāo)得到最優(yōu)的分區(qū)方案。

2.3 基于回溯思想改進(jìn)GN算法的電磁環(huán)網(wǎng)自動(dòng)分區(qū)方法流程

結(jié)合1.3節(jié)、2.1節(jié)、2.2節(jié)，提出基于回溯思想改進(jìn)GN算法的電磁環(huán)網(wǎng)自動(dòng)分區(qū)方法，其流程如圖5所示。

圖5 基于回溯思想改進(jìn)GN算法的電磁環(huán)網(wǎng)自動(dòng)分區(qū)方法流程圖Fig.5 Flowchart of automatic electromagnetic loop network partitioning algorithm based on backtracking-improved GN algorithm

3 算例結(jié)果與分析

3.1 算例1

以某地區(qū)的500 kV/220 kV電磁環(huán)網(wǎng)為例驗(yàn)證本文方法的可行性和有效性。該地區(qū)電網(wǎng)包含60座220 kV變電站，其中有5座與500 kV變電站直接聯(lián)系，稱(chēng)作樞紐變電站，分別記為A—E；其余220 kV變電站分別編號(hào)為1—55。值得注意的是，由于該地區(qū)500 kV變電站的主變數(shù)目均為2臺(tái)，分區(qū)時(shí)需要與其他500 kV變電站“手拉手”運(yùn)行。該地區(qū)解環(huán)前的電網(wǎng)拓?fù)淙鐖D6所示。

圖6 電網(wǎng)拓?fù)鋱DFig.6 Topology of power grid

表1為分區(qū)前部分變電站220 kV母線(xiàn)處的短路電流水平。其中，與樞紐變電站直接聯(lián)系的220 kV母線(xiàn)的斷路器開(kāi)關(guān)遮斷容量為50 kA，其余220 kV變電站的220 kV母線(xiàn)的斷路器開(kāi)關(guān)遮斷容量為40 kA。由表1可以看出，樞紐變電站220 kV母線(xiàn)的單相、三相短路電流普遍超標(biāo)。此外，21站220 kV母線(xiàn)的三相短路電流也超過(guò)了斷路器的遮斷容量。因此，有必要對(duì)該地區(qū)500 kV/220 kV電磁環(huán)網(wǎng)的220 kV環(huán)網(wǎng)進(jìn)行解環(huán)，從而達(dá)到抑制短路電流水平的目的。

表1 分區(qū)前部分站點(diǎn)220 kV母線(xiàn)處的短路電流Table 1 220 kV bus short circuit currents of some substations before partition

3.1.1 GN算法獲得的分區(qū)方案

使用改進(jìn)前的GN算法對(duì)其進(jìn)行分區(qū)的過(guò)程中，依據(jù)工程經(jīng)驗(yàn)設(shè)定最大可斷線(xiàn)路數(shù)目為10，可由GN算法得出前10條斷線(xiàn)的綜合邊介數(shù)信息，如表2所示。在斷線(xiàn)過(guò)程執(zhí)行到第6條線(xiàn)路，即24－26時(shí)，會(huì)形成3個(gè)分區(qū)，其中分區(qū)1包含A站與B站，分區(qū)2僅包含C站，分區(qū)3包含D站與E站，如圖7所示。首先該分區(qū)中具有2臺(tái)主變的500 kV變電站沒(méi)有與其他500 kV變電站“手拉手”運(yùn)行，另外在分區(qū)2中形成了C站、34站、35站、26站的串供結(jié)構(gòu)。由此說(shuō)明改進(jìn)前的GN算法沒(méi)能獲得可行的分區(qū)方案。

表2 綜合邊介數(shù)Table 2 Comprehensive edge betweenness

3.1.2 改進(jìn)GN算法獲得的分區(qū)方案

用基于回溯思想改進(jìn)GN算法進(jìn)行電磁環(huán)網(wǎng)自動(dòng)分區(qū)，不僅通過(guò)1.3節(jié)的方法保證了分區(qū)方案滿(mǎn)足電磁環(huán)網(wǎng)的分區(qū)原則，并且通過(guò)2.1節(jié)提出的回溯思想擴(kuò)展了解的多樣性，依據(jù)工程經(jīng)驗(yàn)設(shè)定最大可斷線(xiàn)路數(shù)目為10，得出的分區(qū)方案如表3所示。采用2.2.1節(jié)所提方法對(duì)方案進(jìn)行校驗(yàn)，得出方案6的母線(xiàn)“N-1”和線(xiàn)路“N-1”指標(biāo)校驗(yàn)不通過(guò)，該方案無(wú)效。

圖7 基于GN算法獲取實(shí)際電網(wǎng)的分區(qū)方案Fig.7 A practical network partitioning scheme by GN algorithm

表3 220 kV環(huán)網(wǎng)分區(qū)方案Table 3 Partitioning schemes of 220 kV loop network

在上述分區(qū)方案中，可以看出改進(jìn)GN算法相比GN算法能夠獲得更多的分區(qū)方案，并且大多數(shù)方案均可通過(guò)靜態(tài)安全分析和短路電流校驗(yàn)。

為得到最優(yōu)分區(qū)方案，由2.2.1節(jié)所提方法處理后的靜態(tài)安全分析指標(biāo)、短路電流指標(biāo)以及模塊度指標(biāo)Q，共計(jì)9個(gè)指標(biāo)，共同構(gòu)成分區(qū)方案優(yōu)劣評(píng)估指標(biāo)集。其中方案1—5的模塊度指標(biāo)Q分別為0.4247、0.4388、0.4050、0.4422 和 0.4300。

采用PCA法并設(shè)置主成分分量為85%，可得9個(gè)指標(biāo)的分量比率和累積率，如圖8所示。進(jìn)而形成反映方案優(yōu)劣的綜合評(píng)價(jià)指標(biāo)值，如表4所示。結(jié)果表明，方案2的綜合評(píng)價(jià)指標(biāo)優(yōu)于其他方案，為該算例的最優(yōu)分區(qū)方案。其余分區(qū)方案評(píng)估排名依次為方案3、方案1、方案5和方案4。

圖8 各指標(biāo)分量占全分量的比率和累積率Fig.8 Proportion of each index component and accumulative proportion

表4 綜合評(píng)價(jià)指標(biāo)值Table 4 Comprehensive evaluation index

3.2 算例2

為了進(jìn)一步驗(yàn)證方法的有效性，本文還以該地區(qū)未來(lái)新投入2座500 kV變電站后的電網(wǎng)為例進(jìn)行類(lèi)似分析并給出相應(yīng)分區(qū)方案。該電網(wǎng)拓?fù)鋱D如圖9所示，共67座220 kV變電站，其中7座樞紐變電站分別編號(hào)為A—F，其他變電站分別編號(hào)為1—60。由于A站對(duì)應(yīng)的500 kV樞紐變電站包含3臺(tái)500 kV主變，因此可獨(dú)立供電一片區(qū)域。

圖9 電網(wǎng)拓?fù)鋱DFig.9 Topology of power grid

基于回溯思想改進(jìn)的GN算法可得經(jīng)校驗(yàn)合格的有效分區(qū)方案共4個(gè)，使用2.2節(jié)所提評(píng)估方法并設(shè)置主成分分量為85%，可得方案1的綜合評(píng)價(jià)指標(biāo)值為0.8388，為最優(yōu)分區(qū)方案，具體方案信息如表5所示。

表5 220 kV環(huán)網(wǎng)分區(qū)方案Table 5 Partitioning schemes of 220 kV loop network

4 結(jié)語(yǔ)

本文提出了一種基于回溯思想改進(jìn)GN算法的電磁環(huán)網(wǎng)自動(dòng)分區(qū)方法。該方法首先針對(duì)GN算法的應(yīng)用缺陷，提出計(jì)及電磁環(huán)網(wǎng)分區(qū)原則的GN算法；同時(shí)，為充分發(fā)揮GN算法的優(yōu)勢(shì)，提出基于回溯思想的改進(jìn)GN算法，以擴(kuò)展可行分區(qū)方案的多樣性；而且，為解決人工校驗(yàn)所存在的弊端，結(jié)合進(jìn)程連接PSD-BPA的方式自動(dòng)對(duì)所得方案進(jìn)行靜態(tài)安全分析和短路電流校驗(yàn)；并將得到的指標(biāo)與模塊度指標(biāo)相結(jié)合，構(gòu)成衡量方案優(yōu)劣的評(píng)估指標(biāo)集；進(jìn)而采用PCA法對(duì)該指標(biāo)集進(jìn)行處理，以評(píng)估出最優(yōu)分區(qū)方案。最后，將本文所提方法應(yīng)用于某地區(qū)500 kV/220 kV電磁環(huán)網(wǎng)，結(jié)果表明，該方法可有效實(shí)現(xiàn)電磁環(huán)網(wǎng)自動(dòng)分區(qū)，并確保所得方案滿(mǎn)足實(shí)際電網(wǎng)分區(qū)原則。

由于GN算法每計(jì)算一次最大綜合邊介數(shù)信息都需遍歷網(wǎng)絡(luò)中所有的站點(diǎn)，在網(wǎng)絡(luò)規(guī)模大到一定程度時(shí)，其時(shí)間成本將顯著提高，因此該算法在時(shí)間效率上尚有一定的局限性，有待進(jìn)一步研究。

參考文獻(xiàn)：

［1］電力行業(yè)電網(wǎng)運(yùn)行與控制標(biāo)準(zhǔn)化技術(shù)委員會(huì).電力系統(tǒng)安全穩(wěn)定導(dǎo)則：DL 755—2001［S］.北京：中華人民共和國(guó)國(guó)家經(jīng)濟(jì)貿(mào)易委員會(huì)：2001.

［2］侯春青，鄭惠萍.2005年山西中南部500 kV及220 kV電網(wǎng)的電磁環(huán)網(wǎng)運(yùn)行方式研究［J］.電網(wǎng)技術(shù)，2005，29（10）：80-84.HOU Chunqing，ZHENG Huiping.Research on parallel operation of 500 kV and 220 kV power networks in south and central part of Shanxi Power Grid in 2005［J］.Power System Technology，2005，29（10）：80-84.

［3］張祖平，范明天，周莉梅.城市電網(wǎng)電磁環(huán)網(wǎng)的解環(huán)問(wèn)題研究［J］.電網(wǎng)技術(shù)，2008，32（19）：42-44.ZHANG Zuping，F(xiàn)AN Mingtian，ZHOU Limei.Research of breaking up electromagnetic ring in urban power network［J］.Power System Technology，2008，32（19）：42-44.

［4］孔濤，王洪濤，劉玉田，等.500 kV-220 kV電磁環(huán)網(wǎng)開(kāi)環(huán)研究［J］.電力自動(dòng)化設(shè)備，2003，23（12）：13-16.KONG Tao，WANG Hongtao，LIU Yutian，et al.Research on untying 500 /220 kV electro-magnetic loop network［J］.Electric Power Automation Equipment，2003，23（12）：13-16.

［5］潘煒，劉文穎，楊以涵，等.電磁環(huán)網(wǎng)條件下西北750 kV電網(wǎng)運(yùn)行方式的研究［J］.電網(wǎng)技術(shù)，2007，31（15）：33-38.PAN Wei，LIU Wenying，YANG Yihan，etal.Research on operating modes of 750 kV Northwest China Power Grid electromagnetically coupled with 330 kV power grid［J］.Power System Technology，2007，31（15）：33-38.

［6］魏震波.復(fù)雜網(wǎng)絡(luò)社區(qū)結(jié)構(gòu)及其在電網(wǎng)分析中的應(yīng)用研究綜述［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2015，35（7）：1567-1577.WEI Zhenbo.Overview of complex networks community structure and its applications in electric powernetwork analysis［J］.Proceedings of the CSEE，2015，35（7）：1567-1577.

［7］倪向萍，阮前途，梅生偉，等.基于復(fù)雜網(wǎng)絡(luò)理論的無(wú)功分區(qū)算法及其在上海電網(wǎng)中的應(yīng)用［J］.電網(wǎng)技術(shù)，2007，31（9）：6-12.NI Xiangping，RUAN Qiantu，MEI Shengwei，et al.A new network partitioning algorithm based on complex network theory and its application in Shanghai Power Grid［J］.Power System Technology，2007，31（9）：6-12.

［8］GIRVAN M，NEWMAN M E J.Community structure in social and biological networks［J］.Proceedings of the National Academy of Sciences，2002，99（12）：7821-7826.

［9］楊冬，劉玉田.基于網(wǎng)絡(luò)社團(tuán)結(jié)構(gòu)特性量化分析的電磁環(huán)網(wǎng)分區(qū)方法［J］.電力系統(tǒng)自動(dòng)化，2013，37（21）：66-71.YANG Dong，LIU Yutian.A new method for electromagnetic loop partitioning based on quantitative analysis ofnetwork community structure characteristics［J］.Automation of Electric Power Systems，2013，37（21）：66-71.

［10］葉華，劉玉田，牛新生.500 kV-220 kV電磁環(huán)網(wǎng)開(kāi)環(huán)方案模糊綜合評(píng)價(jià)［J］.電力自動(dòng)化設(shè)備，2006，26（7）：1-5.YE Hua，LIU Yutian，NIU Xinsheng.Fuzzy comprehensive evaluation of opening schemes for 500 kV-220 kV electromagnetic loop［J］.Electric Power Automation Equipment，2006，26（7）：1-5.

［11］于會(huì)泉.電網(wǎng)短路電流限制措施的優(yōu)化配置研究［D］.北京：華北電力大學(xué)，2011.YU Huiquan.Research on the optimization strategy of shortcircuit current limiting mehods for power networks［D］.Beijing：North China Power Electric University，2011.

［12］王佳明，劉文穎，李群炬，等.基于多層次模糊評(píng)價(jià)模型的短路電流限制措施優(yōu)化配置［J］.電網(wǎng)技術(shù)，2011，35（11）：125-129.WANG Jiaming，LIU Wenying，LI Qunju，et al.Optimal configuration of short-circuit current limitation measures based on multilevel fuzzy comprehensive evaluation model［J］.Power System Technology，2011，35（11）：125-129.

［13］張東妮.電磁環(huán)網(wǎng)解環(huán)方案綜合評(píng)價(jià)研究［D］.重慶：重慶大學(xué)，2011.ZHANG Dongni.Study of comprehensive evaluation in unlocking the electromagnetic loop network［D］.Chongqing：Chongqing University，2011.

［14］聶宏展，聶聳，喬怡，等.基于主成分分析法的輸電網(wǎng)規(guī)劃方案綜合決策［J］.電網(wǎng)技術(shù)，2010，34（6）：134-138.NIE Hongzhan，NIE Song，QIAO Yi，et al.Comprehensive decision-making of alternative transmission network planning based on principal component analysis［J］.Power System Technology，2010，34（6）：134-138.

［15］THOMAS H C，CHARLES E L，RONALD L R，et al.Introduction to algorithms［M］.3rd ed.London，UK：the MIT Press，2009：594-611.

［16］NEWMAN M E J，GIRVAN M.Finding and evaluating community structure in networks［J］.Physical Review E，2004，69（2）：026113.

［17］徐澤水.不確定多屬性決策方法及應(yīng)用［M］.北京：清華大學(xué)出版社，2004：26-27.