郭巖巖
(國網(wǎng)河南省電力公司鄭州供電公司互聯(lián)網(wǎng)部(數(shù)據(jù)中心),河南 鄭州 450052)
在電網(wǎng)結(jié)構(gòu)優(yōu)化、故障恢復(fù)和網(wǎng)絡(luò)重構(gòu)過程中,各種開關(guān)狀態(tài)的組合將形成不同類型的電網(wǎng)拓撲,包括回路和孤島[1-3]。近些年來,電網(wǎng)設(shè)備管理逐步向農(nóng)網(wǎng)用戶與低壓用戶延伸,電網(wǎng)設(shè)備數(shù)量也逐漸上升,極大地增加了電網(wǎng)動態(tài)拓撲結(jié)構(gòu)識別與故障分析的壓力[4]。因此,為了及時識別復(fù)雜電網(wǎng)的動態(tài)拓撲故障,必須區(qū)分不同類型的電網(wǎng)拓撲,以保證電網(wǎng)運行安全。
Waarde等[5]利用測量的輸入/輸出數(shù)據(jù)識別異構(gòu)網(wǎng)絡(luò)的動態(tài)拓撲結(jié)構(gòu),根據(jù)電網(wǎng)測量數(shù)據(jù)確認重建拓撲結(jié)構(gòu)的可識別性條件;分析SISO系統(tǒng)同質(zhì)異構(gòu)網(wǎng)絡(luò)的可識別特性,利用廣義西爾維斯特方程識別拓撲網(wǎng)絡(luò)結(jié)構(gòu)。缺點是無法識別孤點故障。劉超等[6]提出一種基于高級量測體系的圖模型近鄰估計的配電網(wǎng)拓撲辨識方法,將相鄰時刻電壓幅值量測之差視作高斯隨機變量,構(gòu)建概率圖模型的精度矩陣估計模型;采用近鄰估計算法求解圖模型精度矩陣;采用生成樹算法重建出配電網(wǎng)拓撲。缺點是該方法能夠識別出的故障較少。高澤璞等[7]提出了基于知識圖譜的電網(wǎng)拓撲結(jié)構(gòu)辨識技術(shù),但隨著電網(wǎng)拓撲結(jié)構(gòu)中電源點數(shù)量增長,其識別率會逐漸降低。劉廣一等[8]提出的基于圖計算的電網(wǎng)結(jié)構(gòu)分析技術(shù),可實時識別電網(wǎng)拓撲結(jié)構(gòu),識別效率高,但識別精度稍低。
為了提高電網(wǎng)動態(tài)拓撲結(jié)構(gòu)的識別準確率,本文提出一種基于粒子群優(yōu)化算法的電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別技術(shù)。
本文根據(jù)數(shù)據(jù)源的特點,在分析電網(wǎng)動態(tài)拓撲結(jié)構(gòu)特點的基礎(chǔ)上,設(shè)計了電網(wǎng)動態(tài)拓撲圖數(shù)據(jù)庫的構(gòu)建方法,使電網(wǎng)動態(tài)拓撲圖數(shù)據(jù)庫更加實用和專業(yè)。電網(wǎng)動態(tài)拓撲結(jié)構(gòu)圖數(shù)據(jù)庫的具體構(gòu)建過程如圖1所示。
圖1 電網(wǎng)動態(tài)拓撲結(jié)構(gòu)圖數(shù)據(jù)庫構(gòu)建過程圖
由圖1可知,具體的構(gòu)建過程如下:
1)基于電網(wǎng)拓撲結(jié)構(gòu)的CIM(common information mode),創(chuàng)建各物理設(shè)備的聯(lián)系,并對物理設(shè)備的本身屬性予以考量,選取自頂向下的圖數(shù)據(jù)庫構(gòu)建方式。
2)對于采用結(jié)構(gòu)化數(shù)據(jù)庫的電網(wǎng)信息系統(tǒng),由于結(jié)構(gòu)化數(shù)據(jù)庫的知識易于提取,因此在提取數(shù)據(jù)庫數(shù)據(jù)的知識時,還需要提取數(shù)據(jù)庫中的字段。
3)通過知識融合,從多個信息系統(tǒng)結(jié)構(gòu)化數(shù)據(jù)庫中提取知識,可以提高每個信息系統(tǒng)中同一實體屬性的完整性。
4)建立電網(wǎng)動態(tài)拓撲圖數(shù)據(jù)庫中的物理設(shè)備本體模型。該模型本體屬于物理設(shè)備的抽象概念,包含了物理設(shè)備的所有屬性??梢匀藶榈囟x本體之間的關(guān)系,并盡可能減少本體知識集合。
5)在本體創(chuàng)建的基礎(chǔ)上,創(chuàng)建的本體對應(yīng)數(shù)據(jù)庫中提取的知識,兩者融合形成數(shù)據(jù)實體;通過訓(xùn)練樣本數(shù)據(jù)構(gòu)造實體之間的屬性和關(guān)系。
6)組合構(gòu)成后的實體、屬性與關(guān)系為三元組,三元組的形式包括實體—關(guān)系—屬性與實體—關(guān)系—實體兩種,將所形成的三元組儲存到相應(yīng)的圖數(shù)據(jù)結(jié)構(gòu)內(nèi),實現(xiàn)電網(wǎng)動態(tài)拓撲結(jié)構(gòu)圖數(shù)據(jù)庫的構(gòu)建。
在所構(gòu)建的電網(wǎng)動態(tài)拓撲結(jié)構(gòu)圖數(shù)據(jù)庫基礎(chǔ)上,本文所設(shè)計的電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別系統(tǒng)由圖數(shù)據(jù)庫構(gòu)建層與電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別層構(gòu)成,其中圖數(shù)據(jù)庫構(gòu)建層主要包括數(shù)據(jù)源獲取單元、數(shù)據(jù)轉(zhuǎn)換處理單元及數(shù)據(jù)存儲單元?;诹W尤簝?yōu)化算法的電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別技術(shù)框架如圖2所示。
圖2 電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別技術(shù)框架
在圖2中,圖數(shù)據(jù)庫構(gòu)建層通過數(shù)據(jù)源獲取單元,在結(jié)構(gòu)化數(shù)據(jù)庫內(nèi)獲取知識數(shù)據(jù)源,由數(shù)據(jù)轉(zhuǎn)換處理單元對所獲取的數(shù)據(jù)源實施相應(yīng)的轉(zhuǎn)換處理后,存儲到數(shù)據(jù)存儲單元中,實現(xiàn)圖數(shù)據(jù)庫的構(gòu)建;電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別層重點依據(jù)圖數(shù)據(jù)庫內(nèi)經(jīng)過圖數(shù)據(jù)化之后的電網(wǎng)動態(tài)拓撲結(jié)構(gòu)數(shù)據(jù),識別電網(wǎng)動態(tài)拓撲結(jié)構(gòu)。
運用深度優(yōu)先遍歷算法智能識別電網(wǎng)動態(tài)拓撲結(jié)構(gòu)的總體思路是將電網(wǎng)潮流首端視為電源點,由電源點開始,逐個遍歷標記,并對各個相鄰點實施逐個遍歷。如果遍歷過程中某個負荷節(jié)點被訪問兩次或搜尋到其余電源點,則可判別該電網(wǎng)動態(tài)拓撲結(jié)構(gòu)內(nèi)存在環(huán)路;如果遍歷完成后依然存在未被訪問過的負荷節(jié)點,那么可判斷該負荷節(jié)點屬于電源孤立節(jié)點;如果在完成遍歷之后沒有出現(xiàn)環(huán)路或孤立節(jié)點,則此電網(wǎng)動態(tài)拓撲結(jié)構(gòu)屬于常規(guī)運行的開式、輻射網(wǎng)狀。
為了提升深度優(yōu)先遍歷算法識別電網(wǎng)動態(tài)拓撲結(jié)構(gòu)的效率,先采用粒子群優(yōu)化算法提升圖中節(jié)點的訪問速度。
設(shè)原數(shù)組visited[0,1,…,n-1]為粒子群,設(shè)每個節(jié)點為單個粒子i,節(jié)點的路徑訪問環(huán)路判斷系數(shù)為d,則粒子群在優(yōu)化計算中的訪問速度為(vi,d),計算過程如下:
(vi,d)={ed/p(ed,i)|ed∈Ed)}
(1)
式中:(vi,d)為單個節(jié)點被優(yōu)化的訪問速度集合;Ed為路徑中所有節(jié)點的集合;ed為Ed中的個體元素;p(ed,i)為ed的概率。若c×rand≥1,則將概率取為1,否則直接保留為c×rand,以此判斷某個體元素被訪問的概率:
(2)
式中:c為粒子群算法的學(xué)習(xí)因數(shù);rand為隨機概率;A,B分別為兩個節(jié)點的最優(yōu)位置;p(ed,B)為由最優(yōu)位置節(jié)點及概率p(ed)組成的速度集合。將兩集合中最大的概率max(p(ed,1),p(ed,2))加入到新的速度集合中,生成v1,d,v2,d,分別代表速度集合中不同的速度元素,此時速度集合可轉(zhuǎn)化為:
(3)
式中:p(ed,1),p(ed,2)分別為位置節(jié)點為1和2時的被訪問概率。
比較隨機概率rand與vi,d被選取的概率,將vi,d中大于rand的概率值組成集合,得到候選集合cut(vi,d),再從中得到最優(yōu)速度節(jié)點位置元素集:
path[vi,d]={vi,d|vi,d/cut(vi,d)∈
vi,d&cut(vi,d)≥rand}
(4)
式中:path[vi,d]為擁有遍歷最優(yōu)速度的節(jié)點訪問環(huán)路函數(shù)。
在path[vi,d]上,采用深度優(yōu)先遍歷算法智能識別電網(wǎng)動態(tài)拓撲結(jié)構(gòu)的主要過程,如圖3所示。
圖3 電網(wǎng)動態(tài)拓撲結(jié)構(gòu)智能識別過程圖
以某市電力公司2016年的歷史電網(wǎng)拓撲數(shù)據(jù)為例,運用本文所提技術(shù)實施電網(wǎng)動態(tài)拓撲結(jié)構(gòu)識別,通過分析識別結(jié)果檢驗本文技術(shù)的實際應(yīng)用效果。該電力公司電網(wǎng)拓撲數(shù)據(jù)可達11億頂點與22億無向邊的全聯(lián)通圖數(shù)據(jù)。
由電力公司2016年歷史電網(wǎng)拓撲數(shù)據(jù)中隨機選出一個復(fù)雜電網(wǎng),此電網(wǎng)網(wǎng)架結(jié)構(gòu)內(nèi)包含4個電源點與46個負荷節(jié)點,該電網(wǎng)網(wǎng)架結(jié)構(gòu)如圖所示。圖4中U1~U4為4個電源點。
圖4 復(fù)雜電網(wǎng)網(wǎng)架結(jié)構(gòu)圖
在該復(fù)雜電網(wǎng)網(wǎng)架結(jié)構(gòu)的動態(tài)變化過程中,通過本文技術(shù)對其動態(tài)拓撲結(jié)構(gòu)實施智能識別,所識別出存在不同故障的動態(tài)拓撲結(jié)構(gòu)如圖5所示,其中識別出的故障負荷節(jié)點通過空心圓(?)標注。
通過圖5的識別結(jié)果可知,圖5(a)中,該復(fù)雜電網(wǎng)動態(tài)拓撲結(jié)構(gòu)內(nèi)存在一類環(huán)路,此種狀況下拓撲結(jié)構(gòu)中的負荷節(jié)點11、12、43、42、36、30、37、31、32、33、34屬于雙向電源供電,電源點分別為U2和U3;圖5(b)中,復(fù)雜電網(wǎng)動態(tài)拓撲結(jié)構(gòu)內(nèi)存在二類環(huán)路,此時拓撲結(jié)構(gòu)內(nèi)的負荷節(jié)點14、15、36、42、43、12、11均由電源點U2供電;圖5(c)中,該復(fù)雜電網(wǎng)動態(tài)拓撲結(jié)構(gòu)中存在孤點,孤點為負荷節(jié)點36,可知此時沒有電源對負荷節(jié)點36實施供電。
圖5 不同電網(wǎng)動態(tài)拓撲結(jié)構(gòu)識別結(jié)果
在固定時間內(nèi)分別運用本文技術(shù)、基于知識圖譜的電網(wǎng)拓撲結(jié)構(gòu)辨識技術(shù)(文獻[7]技術(shù))、基于圖計算的電網(wǎng)結(jié)構(gòu)分析技術(shù)(文獻[8]技術(shù))對電力公司2016年的歷史電網(wǎng)拓撲數(shù)據(jù)實施電網(wǎng)動態(tài)拓撲結(jié)構(gòu)識別,檢驗各種技術(shù)識別拓撲結(jié)構(gòu)中存在孤點、一類環(huán)路與二類環(huán)路3種故障的識別率并實施對比分析,結(jié)果見表1。
表1 各技術(shù)識別率對比結(jié)果
通過表1可得出,本文技術(shù)、文獻[7]技術(shù)與文獻[8]技術(shù)識別三種故障的平均識別率依次為98.75%、95.99%、93.75%,可見相較其他兩種技術(shù)而言,本文技術(shù)整體識別率更高。
本文構(gòu)建了電網(wǎng)動態(tài)拓撲結(jié)構(gòu)圖數(shù)據(jù)庫,用此儲存電網(wǎng)動態(tài)拓撲結(jié)構(gòu)有關(guān)數(shù)據(jù);基于圖數(shù)據(jù)庫中的圖數(shù)據(jù),利用粒子群算法提升圖中節(jié)點的訪問速度,優(yōu)化了深度優(yōu)先遍歷算法,可以智能識別電網(wǎng)動態(tài)拓撲結(jié)構(gòu),進而智能、有效地識別出電網(wǎng)動態(tài)拓撲結(jié)構(gòu)中所存在的孤點、一類環(huán)路及二類環(huán)路故障,且整體識別率高,提高了電網(wǎng)動態(tài)拓撲結(jié)構(gòu)的識別準確率,為及時消除電網(wǎng)故障、保障電網(wǎng)運行安全提供了參考,在電網(wǎng)優(yōu)化運行方面具有實際應(yīng)用意義。但由于條件限制,本文能夠識別的故障有限,未來可以針對電網(wǎng)其他故障類型進行研究。