邵晨穎，劉友波，邵安海，邱 高，高紅均，劉俊勇

（1. 四川大學(xué)電氣工程學(xué)院，四川省成都市 610065；2. 國網(wǎng)福建超高壓公司，福建省廈門市 361004）

0 引言

精確的拓撲信息是配電網(wǎng)各類分析計算的基礎(chǔ)，對配電網(wǎng)的安全運行、靈活控制及經(jīng)濟規(guī)劃有重要意義。近年來，電動汽車等新型設(shè)備在配電網(wǎng)中的滲透率逐漸增加，分布式能源（distributed energy resource，DER）大規(guī)模接入，為實現(xiàn)靈活可靠供電，配電網(wǎng)拓撲變化愈發(fā)頻繁［1-3］。在實際中，配電網(wǎng)測量數(shù)據(jù)經(jīng)常發(fā)生部分缺失，通信延遲和故障、數(shù)據(jù)丟包等問題會導(dǎo)致數(shù)據(jù)采集與監(jiān)控（supervisory control and data acquisition，SCADA）系統(tǒng)和高級量測體系采集的數(shù)據(jù)缺乏精確性與同步性［4］，這也給其拓撲的實時魯棒辨識帶來挑戰(zhàn)。同步相量測量裝置（phasor measurement unit，PMU）能夠精確測量同步電壓電流相量，但成本較高。如何在有限投資和局部可觀條件下實現(xiàn)抗噪、抗數(shù)據(jù)缺失的配電網(wǎng)拓撲魯棒辨識顯得愈發(fā)重要。

不少學(xué)者基于數(shù)據(jù)建立物理判據(jù)［5-7］或模型［8-10］驅(qū)動拓撲辨識，如文獻［5-7］以節(jié)點電壓的關(guān)聯(lián)性指導(dǎo)拓撲辨識；文獻［8-10］將拓撲辨識簡化為混合整數(shù) 線 性 規(guī) 劃（mixed-integer linear programming，MILP）問題。與大多數(shù)利用電壓數(shù)據(jù)的方法不同，文獻［10］使用饋線電流互感器獲取電流數(shù)據(jù)進行拓撲辨識。饋線電流互感器能實時獲得精確電流相量，成本低且無須停電安裝［11-13］，在設(shè)備投資成本方面極具優(yōu)勢。然而，基于數(shù)據(jù)統(tǒng)計關(guān)系的方法容易受DER 接入容量和類型的影響，難以適用于有源配電網(wǎng)；基于物理模型的算法則計算量大，時效性差。

為提高拓撲辨識效率，基于人工智能的無模型辨識方法成為新的研究方向。文獻［14-18］利用監(jiān)督學(xué)習(xí)挖掘量測信息和拓撲間的關(guān)聯(lián)關(guān)系。文獻［19-21］篩選出最有效的電壓量測，以盡量少的數(shù)據(jù)實現(xiàn)高精度的辨識。但是，傳統(tǒng)監(jiān)督學(xué)習(xí)對未知拓撲的識別效果較差，且不少方法需要使用全量測信息或裝設(shè)一定數(shù)量的PMU 才能得到準(zhǔn)確的辨識結(jié)果，經(jīng)濟性較低。

因此，基于少量低成本饋線電流互感器的有源配電網(wǎng)無模型拓撲辨識仍有待深入研究。條件生成對 抗 網(wǎng) 絡(luò)（conditional generative adversarial network，CGAN）已被應(yīng)用于電力系統(tǒng)負荷預(yù)測、量測數(shù)據(jù)重建等電力領(lǐng)域中［22］，而利用其進行拓撲辨識的研究尚少。原始CGAN 存在梯度消失和模式坍塌問題，訓(xùn)練難度較大，為此，本文提出一種基于梯度懲罰優(yōu)化條件生成對抗網(wǎng)絡(luò)（wasserstein CGAN with gradient penalty，WCGAN-GP）的配電網(wǎng)拓撲辨識算法，其特點在于：1）CGAN 利用數(shù)據(jù)參數(shù)化模型在高維超平面中仿射物理模型，類似于專家思維但魯棒性更優(yōu)，對已知拓撲和未知拓撲都能有較好的辨識效果，經(jīng)梯度優(yōu)化改進的CGAN 有效提高了訓(xùn)練的穩(wěn)定性；2）利用觀測窗中的時序電流數(shù)據(jù)，進一步提升拓撲識別準(zhǔn)確度，增強了算法的抗噪和抗數(shù)據(jù)缺失魯棒性；3）基于有限電流觀測量的辨識方式大幅降低了電力系統(tǒng)的設(shè)備投資成本。

1 傳統(tǒng)配電網(wǎng)拓撲辨識模型

1.1 物理模型驅(qū)動的配電網(wǎng)拓撲辨識模型

物理模型驅(qū)動的配電網(wǎng)拓撲辨識模型為：

式 中：i∈Nnode，Nnode為 節(jié) 點 集 合；k∈K，j∈Ni，node，K和Ni，node分別為裝設(shè)量測裝置的節(jié)點集合和與節(jié)點i相連的節(jié)點集合；Xk，m和Xk分別為第k個量測裝置測得的特征值和其對應(yīng)潮流計算值；Ii和Iij分別為節(jié)點i的注入電流和線路ij的電流；Vi和yij分別為節(jié)點i的電壓和線路ij的導(dǎo)納值；I*i為節(jié)點i注入電流的共軛；Si為節(jié)點i所帶負荷；sij和bi分別為線路ij上開關(guān)的狀態(tài)和節(jié)點i與電網(wǎng)其余部分的連接狀態(tài)。

該模型本質(zhì)上是將采集到的量測值與潮流計算值相匹配，解出所有線路的二進制變量sij以推出網(wǎng)絡(luò)拓撲。該優(yōu)化問題是一個混合整數(shù)非線性規(guī)劃問題，求解的時間和空間復(fù)雜度都較高。因此，一般將其轉(zhuǎn)化為更易處理的MILP 形式后再求解，但轉(zhuǎn)化后其時效性問題仍未得以解決。

1.2 監(jiān)督學(xué)習(xí)驅(qū)動的配電網(wǎng)拓撲辨識模型

基于監(jiān)督學(xué)習(xí)的配電網(wǎng)拓撲辨識方法的內(nèi)在機理在于通過離線訓(xùn)練挖掘量測信息與拓撲之間的映射關(guān)系，在線辨識時以實時量測為輸入，輸出預(yù)測拓撲，相關(guān)模型為：

式中：Te和T分別為預(yù)測拓撲和真實拓撲；Nnn(·)為監(jiān)督學(xué)習(xí)網(wǎng)絡(luò)（如深度神經(jīng)網(wǎng)絡(luò)、深度置信網(wǎng)絡(luò)、卷積神經(jīng)網(wǎng)絡(luò)等）的輸出函數(shù)；α為損失項系數(shù)；wn為第n層網(wǎng)絡(luò)權(quán)重；M為網(wǎng)絡(luò)層數(shù)；Lloss為損失函數(shù)。

這種通過直接學(xué)習(xí)輸入與輸出映射的方式容易導(dǎo)致算法的泛化能力不足，針對訓(xùn)練集外的未知拓撲形式識別效果較差。

2 基于生成對抗網(wǎng)絡(luò)的配電網(wǎng)拓撲辨識

2.1 生成對抗網(wǎng)絡(luò)驅(qū)動的配電網(wǎng)拓撲辨識

為了克服上述傳統(tǒng)物理模型的時效性問題和監(jiān)督學(xué)習(xí)的泛化性問題，考慮采用生成對抗網(wǎng)絡(luò)進行配電網(wǎng)拓撲辨識。生成對抗網(wǎng)絡(luò)驅(qū)動的配電網(wǎng)拓撲辨識與監(jiān)督學(xué)習(xí)類似，依托于機器學(xué)習(xí)技術(shù)，其辨識速度要優(yōu)于傳統(tǒng)物理模型驅(qū)動的配電網(wǎng)辨識模型。不同的是，生成對抗網(wǎng)絡(luò)不直接學(xué)習(xí)量測信息與拓撲間的映射關(guān)系，而是訓(xùn)練2 個優(yōu)化目標(biāo)相互對立的網(wǎng)絡(luò)，間接挖掘到拓撲的分布情況，從而增強算法的泛化能力。

生成對抗網(wǎng)絡(luò)由生成器G和判別器D組成。生成器旨在學(xué)習(xí)輸入隨機噪聲與真實樣本分布之間的關(guān)聯(lián)關(guān)系，由隨機噪聲產(chǎn)生符合給定樣本分布的數(shù)據(jù)，以此欺騙判別器。判別器的任務(wù)在于盡可能分辨生成器產(chǎn)生的樣本和真實樣本。通過交替訓(xùn)練，G和D在對抗中不斷優(yōu)化。然而，生成對抗網(wǎng)絡(luò)只能隨機生成拓撲而無法生成某種特定條件下對應(yīng)的拓撲，不能直接使用生成對抗網(wǎng)絡(luò)進行拓撲辨識。

2.2 CGAN 辨識拓撲基本原理

CGAN 在生成器和判別器的輸入中額外加入條件信息。判別器不僅需要區(qū)分生成器生成數(shù)據(jù)和真實數(shù)據(jù)，還需判斷數(shù)據(jù)是否與條件信息相對應(yīng)。CGAN 的目標(biāo)函數(shù)為：

式中：D(·)和G(·)分別為判別器和生成器的輸出；x和z分 別 為 真 實 樣 本 和 噪 聲；Ex，Pdata(·)和Ez，Pz(·)分 別為x在真實樣本分布中采樣和z在噪聲分布中采樣的期望值；c為條件信息。

考慮到配電網(wǎng)不同拓撲在不同負荷場景下的線路電流觀測量也有所不同，故配電網(wǎng)拓撲與電流、負荷信息間存在一定的對應(yīng)關(guān)系。如果將電流相量和其他量測數(shù)據(jù)作為條件信息，相應(yīng)的拓撲信息作為樣本輸入，則CGAN 的生成器就能夠?qū)W習(xí)到與量測數(shù)據(jù)相關(guān)的拓撲樣本分布情況，從而生成預(yù)測拓撲。在配電網(wǎng)拓撲辨識場景下，CGAN 的作用機理如圖1 所示。

圖1 CGAN 原理圖Fig.1 Schematic diagram of CGAN

3 基于改進生成對抗網(wǎng)絡(luò)的配電網(wǎng)拓撲辨識

3.1 梯度懲罰優(yōu)化的CGAN

一方面，CGAN 使用Jensen-Shannon（JS）散度來衡量真實數(shù)據(jù)分布和生成器生成數(shù)據(jù)分布之間的距離，而任意2 個無重合的分布之間的JS 散度恒為log 2，當(dāng)判別器訓(xùn)練到一定程度，真實拓撲分布和生成器估計拓撲分布無重合時，生成器的損失函數(shù)將保持恒定，造成梯度消失而無法繼續(xù)訓(xùn)練，此時模型將失去辨識能力；另一方面，在CGAN 訓(xùn)練過程中，判別器與生成器很難同時收斂，容易出現(xiàn)模式崩塌［23］，導(dǎo)致生成器只能生成少數(shù)幾種拓撲，無法實現(xiàn)電流測量信息與拓撲的正確對應(yīng)。

不 同 于JS 散 度，Wasserstein 距 離 對2 個 無 重 疊的分布仍能精準(zhǔn)反映它們的遠近，從而為生成器擬合拓撲數(shù)據(jù)分布提供精確的訓(xùn)練方向［24］。因此，采用Wasserstein 距離定義生成器的損失函數(shù)可以有效解決上述問題。同時，使用梯度懲罰實現(xiàn)Wasserstein 距離中的1-Lipschitz 條件限制，從而得到WCGAN-GP 的優(yōu)化函數(shù)為：

式 中：λ為懲罰 項系數(shù)；xdata和xg分別為真實拓撲數(shù)據(jù)和生成器生成拓撲數(shù)據(jù)；ε為一個隨機數(shù)，服從［0，1］上 的 均 勻 分 布；Ex?(·) 為x? 的 期 望 函 數(shù)；||?x?D(x?)||2為將x?輸入判別器后其梯度的二范數(shù)。

3.2 改進生成對抗網(wǎng)絡(luò)驅(qū)動的配電網(wǎng)拓撲辨識模型

利用物理模型的拓撲辨識方法計算效率低、模型更新困難，基于監(jiān)督學(xué)習(xí)的數(shù)據(jù)驅(qū)動方法雖然識別速度快，但無法準(zhǔn)確辨識訓(xùn)練集外的未知拓撲。為解決上述問題，提出了一種基于WCGAN-GP 的配電網(wǎng)拓撲辨識算法。WCGAN-GP 通過學(xué)習(xí)真實拓撲樣本的分布，把握拓撲數(shù)據(jù)的分布特征，將電流和負荷量測作為條件變量引入，挖掘配電網(wǎng)有限觀測量與拓撲之間復(fù)雜的非線性關(guān)系，對訓(xùn)練集中存在和不存在的拓撲都能精確辨識。 將原始WCGAN-GP 模型應(yīng)用到配電網(wǎng)拓撲辨識場景中，得到該模型的數(shù)學(xué)形式，如式（7）所示，其原理圖如圖2 所示。

圖2 適用于配電網(wǎng)拓撲辨識的WCGAN-GP 模型Fig.2 WCGAN-GP model of topology identification for distribution network

式中：t∈H，H為裝設(shè)饋線電流互感器的線路集合；F和Fe分別為判別器對真實拓撲和生成器生成拓撲輸出的判別結(jié)果；It，m和θt，m分別為第t條裝設(shè)饋線電流互感器的線路電流幅值和相角測量值；PL，m和QL，m分別為節(jié)點總有功和無功偽測量值；EGP為式（5）中的懲罰項。

1）數(shù)據(jù)集構(gòu)建：在不同光-荷場景下對真實拓撲進行潮流計算，得到拓撲對應(yīng)的部分線路電流幅值、相角和節(jié)點總負荷構(gòu)成的配電網(wǎng)量測數(shù)據(jù)集。

2）離線訓(xùn)練：將量測數(shù)據(jù)作為條件信息與服從高斯分布的隨機噪聲輸入生成器，生成器生成估計拓撲；再將真實拓撲、估計拓撲和量測數(shù)據(jù)輸入判別器，其判別結(jié)果反饋給生成器和判別器，二者根據(jù)反饋結(jié)果分別更新自身參數(shù)進行優(yōu)化，提高各自的生成能力和判別能力，最終得到基于WCGAN-GP 的拓撲辨識器。

3）在線辨識：將饋線電流互感器實時量測的電流數(shù)據(jù)和節(jié)點總負荷偽測量數(shù)據(jù)輸入訓(xùn)練好的模型中，生成器將會生成預(yù)測拓撲。

配電網(wǎng)拓撲變化（如計劃停運、故障隔離與復(fù)電、重構(gòu)轉(zhuǎn)供等）頻率往往遠小于量測設(shè)備采集頻率，因此，在實際中可以在2 個相鄰的拓撲狀態(tài)間采集到連續(xù)時間斷面的量測數(shù)據(jù)。對于數(shù)據(jù)驅(qū)動的辨識方法而言，量測誤差的存在容易影響其辨識精度。為提升模型對量測數(shù)據(jù)誤差的魯棒性，在量測誤差較大時也能保證辨識的準(zhǔn)確率，采用連續(xù)觀測窗內(nèi)的多個斷面的數(shù)據(jù)進行拓撲辨識。雖然開關(guān)狀態(tài)從一個測量斷面到另一個測量斷面是固定的，但是由于每個節(jié)點所帶負載和光伏出力的隨機性，電流測量值也是隨機變化的。在這種情況下，每個測量斷面都是一個新的隨機場景，為估計拓撲提供了冗余，有利于更好地緩解測量誤差的影響。

3.3 模型結(jié)構(gòu)設(shè)計

在配電網(wǎng)拓撲估計中，多斷面的量測數(shù)據(jù)是一個時序面板數(shù)據(jù)，為了更好地挖掘其中隱含的信息，提高WCGAN-GP 的穩(wěn)定性、收斂速度和生成器生成數(shù)據(jù)的質(zhì)量［25］，引入一維卷積神經(jīng)網(wǎng)絡(luò)（onedimensional convolutional neural network，1D-CNN）構(gòu)造生成器模型，判別器采用全連接網(wǎng)絡(luò)。WCGAN-GP 的具體結(jié)構(gòu)參數(shù)如附錄A 表A1 所示。由于生成器需要產(chǎn)生0-1 變量表示的估計拓撲，故使用Sigmoid 作為最后一層全連接層的激活函數(shù)。

3.4 數(shù)據(jù)的生成與處理

3.4.1 數(shù)據(jù)生成

WCGAN-GP 的訓(xùn)練集和測試集數(shù)據(jù)生成步驟如下。

步驟1：基于文獻［26］中提出的配電網(wǎng)輻射狀約束，對有效開關(guān)狀態(tài)組合進行輪詢，生成一定數(shù)量的輻射狀配電網(wǎng)拓撲；

步驟2：利用拉丁超立方采樣（Latin hypercube sampling，LHS），生成負荷和光伏出力在基準(zhǔn)值的80%～120%范圍內(nèi)波動的光-荷場景；

步驟3：將拓撲與光-荷場景組合，分別對各個樣本進行潮流計算，得到線路電流幅值、相角和節(jié)點負荷結(jié)果；

步驟4：在樣本潮流計算結(jié)果上加上高斯噪聲干擾，以增強WCGAN-GP 的抗噪性能。

3.4.2 數(shù)據(jù)預(yù)處理

1）數(shù)據(jù)標(biāo)準(zhǔn)化

數(shù)據(jù)在輸入WCGAN-GP 前，由于尺度、量級、單位有所不同，需要對其進行標(biāo)準(zhǔn)化處理，將各維數(shù)據(jù)變換為服從均值為0，方差為1的正態(tài)分布，表達式為：

式中：bnorm和b分別為標(biāo)準(zhǔn)化后和標(biāo)準(zhǔn)化前的數(shù)據(jù)；μ為均值；σ為方差。

2）缺失值處理

實時數(shù)據(jù)采集時，可能因通信問題或設(shè)備故障而造成部分數(shù)據(jù)缺失。由于模型的輸入數(shù)據(jù)中包含了一個連續(xù)觀測窗內(nèi)的多個斷面數(shù)據(jù)，這些斷面數(shù)據(jù)都對應(yīng)于同一個拓撲。受益于饋線電流互感器較高的采集頻率，收集足夠斷面數(shù)據(jù)所需時間較短，節(jié)點負荷在這期間一般不會出現(xiàn)較大波動，斷面數(shù)據(jù)間的差異不會太大。因此，采用其余時間斷面數(shù)據(jù)的平均值補全缺失數(shù)據(jù)，有利于降低數(shù)據(jù)缺失對拓撲辨識帶來的影響。

3.5 模型訓(xùn)練

WCGAN-GP 模型的具體訓(xùn)練步驟如下。

步驟1：根據(jù)3.4 節(jié)生成大量數(shù)據(jù)并對其進行預(yù)處理，劃分出訓(xùn)練集、驗證集、已知拓撲測試集和未知拓撲測試集。

步驟2：根據(jù)3.3 節(jié)搭建G和D并隨機初始化。

步驟3：利用訓(xùn)練集數(shù)據(jù)，先訓(xùn)練ndiscr次D（固定G參數(shù)，更新D參數(shù)），然后訓(xùn)練1 次G（固定D參數(shù)，更新G參數(shù)），生成器和判別器的損失函數(shù)分別如式（9）和式（10）所示，二者根據(jù)其損失函數(shù)進行隨機梯度下降更新。

式中：LG和LD分別為生成器和判別器的損失函數(shù)。

步驟4：重復(fù)步驟3，交替訓(xùn)練G和D，直至達到納什平衡，即LD不再下降時結(jié)束訓(xùn)練。

3.6 模型評價指標(biāo)

為了評估模型的辨識效果，模型的評價指標(biāo)為：

式中：η為拓撲辨識的準(zhǔn)確率；nCT為辨識正確的拓撲數(shù)；nET為辨識錯誤的拓撲數(shù)。

4 算例分析

仿真中使用OpenDSS 實現(xiàn)對配電網(wǎng)的潮流計算，利用Python 的Keras+Tensorflow 架構(gòu)搭建并訓(xùn)練WCGAN-GP，所使用的處理器型號為Intel Core i5-11400，主頻為2.60 GHz，內(nèi)存為16 GB。在改造的IEEE 33 節(jié)點和IEEE 69 節(jié)點系統(tǒng)上分別進行測試，系統(tǒng)結(jié)構(gòu)圖如附錄A 圖A1 所示。

為了考慮DER 的影響，在原始結(jié)構(gòu)中的部分節(jié)點接入一定容量的光伏，配置如附錄A 表A2 所示。使用含恒阻抗、恒電流、恒功率的靜態(tài)負荷（ZIP）模型模擬系統(tǒng)中各類負荷的存在，且滿足占比PZ∶PI∶PP=QZ∶QI∶QP=0.2∶0.3∶0.5。由文獻［10］可知，為了實現(xiàn)配電網(wǎng)拓撲準(zhǔn)確辨識，每個環(huán)路上都至少應(yīng)配置一個饋線電流互感器。因此，在IEEE 33 節(jié)點系統(tǒng)放置了5 個饋線電流互感器。IEEE 69 節(jié)點系統(tǒng)結(jié)構(gòu)更加復(fù)雜，為了加速WCGAN-GP 的訓(xùn)練和收斂，在每個環(huán)路上都放置了2 個饋線電流互感器。

根據(jù)3.4.1 節(jié)，樣本生成情況如附錄A 表A3 所示。其中，訓(xùn)練集和驗證集中不包含未知拓撲測試集中的拓撲，只包含已知拓撲測試集中的拓撲。對訓(xùn)練集樣本隨機加入5%、10%、15%、20%這4 種高斯噪聲，以增強模型的泛化性能。WCGAN-GP訓(xùn)練參數(shù)取值如表A4 所示。

為了驗證所提算法的優(yōu)越性，將本文方法（WCGAN-GP）與貝葉斯網(wǎng)絡(luò)［14］（Bayasian network，BN）、融合注意力機制的卷積神經(jīng)網(wǎng)絡(luò)［17］（attention mechanism and convolutional neural network，ACNN）、深 度 神 經(jīng) 網(wǎng) 絡(luò)［19］（deep neural network，DNN）3 種監(jiān)督學(xué)習(xí)算法和基于物理模型的MILP［10］一起進行多項測試并對比測試結(jié)果。

4.1 準(zhǔn)確性測試

線路電流相量可由饋線電流互感器測量，根據(jù)可用互感器的不同類型，電流幅值的測量誤差可達1%至3%，電流相角的誤差可達1°至5°［10］。節(jié)點總負荷偽測量值可由SCADA 系統(tǒng)預(yù)測得到。因此，考慮實際量測中可能存在的誤差，在測試集中加入不同高斯噪聲組合，以e1表示饋線電流互感器量測的電流幅值和相角誤差，e2表示節(jié)點負荷偽測量誤差，模擬在不同電流測量誤差和偽測量誤差組合下模型的識別性能，已知拓撲測試集上的測試結(jié)果如圖3 所示。

圖3 不同誤差組合下識別效果Fig.3 Identification effects of different error combinations

由圖3 可知，5 種算法在測量誤差較小的情況下都能有較好的辨識效果，WCGAN-GP 的辨識性能最好，是唯一能實現(xiàn)100%辨識準(zhǔn)確率的算法。此外，基于物理模型的MILP 和基于傳統(tǒng)監(jiān)督學(xué)習(xí)的ACNN、DNN、BN 在測量誤差增大后，辨識準(zhǔn)確率都有明顯下降，而WCGAN-GP 即使在e1=10%和e2=50%的情況下，仍能保證93.6%以上的準(zhǔn)確率，可見測量誤差對WCGAN-GP 辨識拓撲的影響較小。這主要是因為，一方面，所提模型使用了連續(xù)觀測窗內(nèi)的多斷面量測數(shù)據(jù)，有效減少測量誤差可能帶來的估計錯誤；另一方面，生成對抗網(wǎng)絡(luò)中的判別器相當(dāng)于提供了一種根據(jù)任務(wù)、數(shù)據(jù)集不同而自適應(yīng)的損失，與損失形式相對固定的傳統(tǒng)監(jiān)督學(xué)習(xí)方式相比，這種自適應(yīng)的生成對抗損失訓(xùn)練得到的模型精度更高。

4.2 魯棒性測試

為了驗證所提WCGAN-GP 模型的魯棒性，在e1=1%和e2=5%的情況下，分別測試了DER 滲透率、數(shù)據(jù)缺失和負荷模型對拓撲辨識準(zhǔn)確率的影響。

4.2.1 DER 滲透率的影響

考慮實際系統(tǒng)中接入的DER，測試了DER 滲透率在10%～30%下各種算法的辨識效果，測試結(jié)果如附錄A 圖A2 所示。各方法對DER 滲透率的魯棒性都較強，但整體來看WCGAN-GP 的識別效果更好，DER 滲透率不超過25%時，準(zhǔn)確率都在99%以上。隨著DER 滲透率增至30%，系統(tǒng)潮流出現(xiàn)反向，靠近接入DER 節(jié)點的線路電流由末端流向首端，造成WCGAN-GP 模型的識別效果略有下降。

4.2.2 數(shù)據(jù)缺失的影響

實際采集量測數(shù)據(jù)時可能會出現(xiàn)數(shù)據(jù)丟失的情況?？紤]到拓撲變化過于頻繁時可能無法獲得足夠的斷面數(shù)據(jù)，對WCGAN-GP 在數(shù)據(jù)缺失和斷面缺失的情況下進行測試。缺失率定義為：

式中：ρ為數(shù)據(jù)缺失率或斷面缺失率；Nl為缺失的數(shù)據(jù)量或斷面數(shù)；N為總數(shù)據(jù)量或斷面數(shù)。

測試時，隨機選取一定數(shù)量的數(shù)據(jù)置零模擬實際中存在的缺失數(shù)據(jù)。不同缺失率下的辨識準(zhǔn)確率如圖4 所示。

圖4 不同數(shù)據(jù)缺失率下的識別效果Fig.4 Identification effects of different data missing rates

顯然，數(shù)據(jù)缺失率小于70% 時，WCGAN-GP的斷面缺失辨識準(zhǔn)確率在99%以上，數(shù)據(jù)缺失辨識準(zhǔn)確率在86%以上。這是因為數(shù)據(jù)缺失會使其冗余度降低，采用一個連續(xù)觀測窗內(nèi)其余斷面數(shù)據(jù)的平均值補全缺失數(shù)據(jù)可以盡可能降低數(shù)據(jù)缺失的影響，保證冗余度，使其對數(shù)據(jù)缺失的魯棒性較強。MILP 基于物理模型，對數(shù)據(jù)依賴性強，因此數(shù)據(jù)稍有缺失都會對辨識結(jié)果產(chǎn)生較大的影響；BN 雖然隨著缺失率的升高識別準(zhǔn)確率整體變化不大，但其整體識別效果不如WCGAN-GP；DNN 和ACNN 由于沒有數(shù)據(jù)補全預(yù)處理，故只要出現(xiàn)數(shù)據(jù)缺失就完全無法正確辨識拓撲。

4.2.3 負荷模型的影響

配電網(wǎng)中存在各種類型的負荷，為了研究不同負荷占比對辨識效果的影響，對系統(tǒng)接入不同比例的ZIP 負荷，測試結(jié)果如附錄A 表A5 和表A6 所示。可見，WCGAN-GP 模型的識別效果在各種ZIP 負荷場景下只有輕微變化，且準(zhǔn)確率均高于其余算法。

4.3 對未知拓撲的適應(yīng)性測試

為了證明WCGAN-GP 模型對訓(xùn)練集中不存在的拓撲類型的辨識適應(yīng)性，用未知拓撲測試集對各種方法進行測試，結(jié)果如圖5 所示。

圖5 各算法對未知拓撲識別效果Fig.5 Identification effects of each algorithm for unknown topologies

圖5 中，BN 選擇拓撲庫中可能性最高的拓撲作為預(yù)測拓撲，由于拓撲庫中不存在未知拓撲，其辨識準(zhǔn)確率始終為0；而DNN 和ACNN 作為傳統(tǒng)監(jiān)督學(xué)習(xí)算法，直接學(xué)習(xí)量測信息與拓撲間的映射關(guān)系，對未知拓撲的適應(yīng)性較差。WCGAN-GP 則通過學(xué)習(xí)真實拓撲的分布特征，并將其與量測數(shù)據(jù)關(guān)聯(lián)起來，生成器可以輸出和訓(xùn)練樣本中的真實拓撲在某些特征（由判別器提取得到）上保持一致的預(yù)測拓撲，在未知拓撲測試集上具有更高的準(zhǔn)確率，解決了傳統(tǒng)監(jiān)督學(xué)習(xí)算法在未知拓撲辨識上的缺陷與不足，這一點在量測誤差增大后愈發(fā)明顯。

4.4 時效性分析

配電網(wǎng)拓撲辨識對于實時性具有一定的要求，其辨識時間必須小于拓撲變化周期。使用饋線電流互感器采集電流數(shù)據(jù)時，其采集周期可達毫秒級［13］。一般情況下，WCGAN-GP 在一個拓撲變化周期內(nèi)完全可以收集到足夠多的斷面數(shù)據(jù)。在這個基礎(chǔ)上，對比WCGAN-GP 和其他算法的在線識別時間，結(jié)果如附錄A 表A7 所示。WCGAN-GP 的辨識速度雖不如傳統(tǒng)監(jiān)督學(xué)習(xí)，但與MILP 相比仍有較大的優(yōu)勢。WCGAN-GP 在系統(tǒng)規(guī)模變大時，識別時間變化不明顯，MILP 的識別時間則與系統(tǒng)節(jié)點數(shù)有很大關(guān)系。

4.5 經(jīng)濟性分析

目前，配電網(wǎng)分布式同步相量測量裝置（distribution PMU，DPMU）的配置成本遠高于饋線電流互感器［10，27］。DPMU 需要考慮停電安裝、調(diào)試造成的經(jīng)濟損失，饋線電流互感器不僅安裝時無須停電，不存在停電損失，而且安裝簡便，只需要將其鉗裝在中壓架空線或開關(guān)柜電纜相應(yīng)位置即可。不同拓撲辨識方法所須裝設(shè)的量測設(shè)備類型及數(shù)目如附錄A 表A8 所示。可見，從經(jīng)濟性出發(fā)，MILP 在實際應(yīng)用中所需的設(shè)備投資成本最低，WCGANGP 次之，但綜合辨識準(zhǔn)確性、魯棒性、時效性和經(jīng)濟性來看，WCGAN-GP 的性能最好，有助于實現(xiàn)高性價比的拓撲辨識。

5 結(jié)語

本文將CGAN 應(yīng)用于有源配電網(wǎng)拓撲辨識場景中，提出了基于局部連續(xù)電流相量和WCGANGP 的拓撲辨識方法。利用拓撲信息在不同光-荷場景下潮流計算得到的部分線路電流幅值、相角和系統(tǒng)總有功、無功負荷值離線訓(xùn)練模型，并運用饋線電流互感器采集實時數(shù)據(jù)實現(xiàn)在線辨識拓撲。針對IEEE 33 節(jié)點、IEEE 69 節(jié)點算例系統(tǒng)，通過分析WCGAN-GP 的拓撲識別性能并與其他算法進行對比測試，得到結(jié)論如下：

1）傳統(tǒng)監(jiān)督學(xué)習(xí)能通過挖掘量測數(shù)據(jù)與拓撲之間的映射關(guān)系來準(zhǔn)確預(yù)測一定泛化域內(nèi)的已知拓撲，但在泛化域外的未知拓撲測試集上辨識效果欠佳，而WCGAN-GP 能通過歷史數(shù)據(jù)學(xué)習(xí)到與電流和負荷信息相關(guān)的拓撲數(shù)據(jù)分布特征，對未知拓撲識別的準(zhǔn)確率更高，適應(yīng)性更強；

2）算法使用了連續(xù)觀測窗內(nèi)的多斷面數(shù)據(jù)，對應(yīng)于同一拓撲的多個場景數(shù)據(jù)為拓撲辨識提供了冗余，同時多斷面的量測為缺失值的補全提供了數(shù)據(jù)基礎(chǔ)。與其他算法相比，WCGAN-GP 有效提升了算法的識別性能，并大幅降低了實時數(shù)據(jù)測量誤差或缺失的影響，保證在輸入數(shù)據(jù)誤差或缺失率較大時仍能實現(xiàn)高精度的拓撲辨識，表明了本文方法在配電網(wǎng)拓撲辨識時有著較好的魯棒性；

3）基于有限電流觀測數(shù)據(jù)的WCGAN-GP 模型只需裝設(shè)少量低成本的非接觸式饋線電流互感器即可實現(xiàn)拓撲辨識，在實際應(yīng)用中經(jīng)濟性較高。

后續(xù)研究工作中，將考慮饋線電流互感器的布點優(yōu)化問題，旨在以最少的測量點數(shù)量得到最佳的辨識效果，以支撐覆蓋一定供電區(qū)域的網(wǎng)格化有源配電網(wǎng)全局拓撲高準(zhǔn)確度快速辨識。

附錄見本刊網(wǎng)絡(luò)版（http：//www.aeps-info.com/aeps/ch/index.aspx），掃英文摘要后二維碼可以閱讀網(wǎng)絡(luò)全文。