多源信息融合的故障元件診斷*

李郁俠　張凱華　譚雅嵐,　李　石　何小軍　陳　波

(1.西安理工大學(xué)　西安　710048)(2.國(guó)家電網(wǎng)陜西省漢中局　漢中　723099)

?

多源信息融合的故障元件診斷*

李郁俠1張凱華1譚雅嵐1,2李石2何小軍2陳波2

(1.西安理工大學(xué)西安710048)(2.國(guó)家電網(wǎng)陜西省漢中局漢中723099)

摘要分析了故障診斷信息的主要來(lái)源,對(duì)比了各個(gè)信息的相同點(diǎn)和不同點(diǎn),建立了面向元件的故障診斷模型,該模型簡(jiǎn)化了以往的故障診斷計(jì)算的復(fù)雜度,使得故障診斷的效率可以顯著提高,提出了新的故障診斷方法。該方法首先通過(guò)計(jì)算各個(gè)元件的貝葉斯疑似度獲得可疑故障元件,并結(jié)合了基于綜合阻抗的故障線路識(shí)別方法,能夠?qū)收暇€路具體相別進(jìn)行識(shí)別,然后提出保護(hù)動(dòng)作的時(shí)序約束,對(duì)保護(hù)動(dòng)作信息進(jìn)行識(shí)別,從而提高故障診斷的準(zhǔn)備性。通過(guò)SIMULINK驗(yàn)證了模型的有效性及正確性。

關(guān)鍵詞故障診斷; 綜合阻抗; 時(shí)序約束; 貝葉斯疑似度

Class NumberTM41

1引言

作為經(jīng)濟(jì)持續(xù)發(fā)展的支柱型行業(yè),電力工業(yè)對(duì)于國(guó)計(jì)民生的重要性日漸凸顯。電力系統(tǒng)的安全性與穩(wěn)定性在我國(guó)的發(fā)展歷程中將會(huì)一直占據(jù)極高的戰(zhàn)略地位。因此電網(wǎng)在發(fā)生事故后,快速診斷并隔離故障,及時(shí)恢復(fù)對(duì)用戶的可靠安全供電,成為智能電網(wǎng)的一個(gè)研究方向,電網(wǎng)故障診斷方法的研究一直是國(guó)內(nèi)外研究的重點(diǎn)課題。

多類(lèi)型故障數(shù)據(jù)的應(yīng)用促進(jìn)了信息融合技術(shù)在電網(wǎng)故障診斷問(wèn)題中的研究。在當(dāng)前研究中,根據(jù)所用基本原理的不同主要采用WAMS采集電氣量信息結(jié)合小波變換的多源信息融合電網(wǎng)故障診斷方法[1～5]。然而在系統(tǒng)運(yùn)行總存在很多不確定因素導(dǎo)致元件發(fā)生故障的異常情況,按照基于結(jié)合繼電保護(hù)信息優(yōu)化模型[6～7]。對(duì)于線路部分存在很大問(wèn)題,不能精確的診斷故障元件。為了解決這一問(wèn)題,本文在結(jié)合了基于綜合阻抗的故障線路識(shí)別方法下,提出了改進(jìn)的優(yōu)化模型,利用有向二分圖法降低了故障診斷過(guò)程運(yùn)算的復(fù)雜度,提高了診斷速度。

2基于面向的有向二分圖壓

有向二分圖法是對(duì)貝葉斯網(wǎng)絡(luò)的簡(jiǎn)化,可以簡(jiǎn)單明了地表達(dá)出故障與征兆之間的關(guān)系[8]。有向二分圖法具有一定的建模能力,降低了運(yùn)算的復(fù)雜度,能夠彌補(bǔ)貝葉斯網(wǎng)絡(luò)在故障診斷過(guò)程中的一些不足。

有向二分圖由三部分組成： 1) 故障節(jié)點(diǎn)和征兆節(jié)點(diǎn)組成的集合V,V=F∪S,其中S為征兆集,F為故障集； 2) 有向邊集合E,E=F×S,由故障指向征兆； 3) 權(quán)重p(s|f)組成的集合PF×S,PF×S={p(s|f)|f∈F,s∈S},其中p(s|f)指在故障f發(fā)生的條件下,征兆s發(fā)生的概率值,0≤p(s|f)≤1。

圖1　有向二分圖

圖1是由三個(gè)故障節(jié)點(diǎn)、四個(gè)征兆節(jié)點(diǎn)組成的概率加權(quán)二分圖,其中故障集F={f1,f2,f3},征兆集S={s1,s2,s3,s4}。征兆s1可能由故障f1引起,記為Domain(s1)={f1},類(lèi)似的有Domain(s2)={f1,f2},同理有Domain(s3)={f1,f3},Domain(s4)={f2,f3}。故障f1可能引起征兆s1,s2,s3,記為Symptom(f1)={s1,s2,s3},同理有Symptom(f2)={s2,s4},Symptom(f3)={s3,s4}。圖中Pf=(0.004,0.007,0.003)為各故障發(fā)生的先驗(yàn)概率,而故障f發(fā)生的條件下,征兆s發(fā)生的概率值

3貝葉斯疑似度

在上述模型中定義貝葉斯疑似度[9],記作Bsd(f,SN),其代表的意義為觀察到的征兆SN導(dǎo)致某故障f發(fā)生的概率與信息完備條件下該故障發(fā)生的絕對(duì)概率的比值,其值越大,該故障發(fā)生的可能性越大。可由下式計(jì)算[10]：

(1)

其中,

(2)

在診斷過(guò)程中應(yīng)對(duì)Domain(SN)中的每個(gè)故障進(jìn)行貝葉斯疑似度計(jì)算,其中SN對(duì)于為診斷模型中的征兆。以故障元件為對(duì)象,按照從大到小的順序?qū)ζ湄惾~斯疑似度排序,直到Symptom(f)覆蓋了所有觀察到的征兆SN。則可知在該觀察征兆下,按照故障發(fā)生的可能性大小,即可得到最優(yōu)故障假設(shè)集F。

4診斷流程

4.1故障診斷流程

由于故障診斷系統(tǒng)是基于調(diào)度端開(kāi)發(fā)的,因此輸入故障診斷模型的調(diào)度信息均來(lái)自于調(diào)度系統(tǒng)。故障發(fā)生后,上傳至調(diào)度系統(tǒng)的信息具有以下特點(diǎn)： 1) 冗余信息太多,有效信息獲取困難； 2) 關(guān)鍵信息的不確定性即存在的漏報(bào),對(duì)于故障診斷的結(jié)果影響非常大。因此,為了得到快速、可信的診斷結(jié)果,必須對(duì)冗余信息及丟失信息進(jìn)行處理,故障診斷系統(tǒng)才能輸出正確的診斷結(jié)果。本文采用如圖2所示的征兆信息識(shí)別過(guò)程對(duì)來(lái)自調(diào)度端的數(shù)據(jù)進(jìn)行處理。

圖2　信息識(shí)別流程

在圖2的故障信息識(shí)別流程中,首先從電網(wǎng)調(diào)度端獲取調(diào)度數(shù)據(jù),通過(guò)對(duì)電氣量進(jìn)行小波分析,從小波熵的變化判斷得到故障發(fā)生時(shí)刻,從該時(shí)刻開(kāi)始在固定的時(shí)間窗內(nèi)獲取調(diào)度信息。其中通過(guò)開(kāi)關(guān)量獲取停電區(qū)域,通過(guò)電氣量獲取綜合阻抗判斷故障發(fā)生在那個(gè)區(qū)域,對(duì)保護(hù)動(dòng)作信息進(jìn)行判斷,剔除掉保護(hù)信息中的冗余信息。然后進(jìn)行征兆信息識(shí)別,現(xiàn)對(duì)征兆信息的識(shí)別步驟進(jìn)行說(shuō)明。

步驟1固定時(shí)間窗內(nèi)接收警報(bào)信息,包括開(kāi)關(guān)量和保護(hù)量

步驟2動(dòng)作信息識(shí)別

輸入步驟1中接收到的警報(bào)信息,按照保護(hù)動(dòng)作的時(shí)序約束識(shí)別警報(bào)信息。

步驟3診斷模型識(shí)別故障

將步驟2中識(shí)別的警報(bào)信息輸入診斷模型,得到診斷結(jié)果

步驟4故障選相

輸入診斷結(jié)果,若其中含有線路,則通過(guò)基于綜合阻抗原理的線路選相識(shí)別具體故障相別,若不含有線路,則直接輸出步驟3中的診斷結(jié)果。

現(xiàn)以圖所示系統(tǒng)發(fā)生故障后征兆信息識(shí)別過(guò)程為例,說(shuō)明上述征兆信息的識(shí)別過(guò)程。以在固定“時(shí)間窗口”中電力系統(tǒng)中SOE(Sequence Of Event)記錄的具有時(shí)序邏輯的信息作為故障診斷系統(tǒng)的數(shù)據(jù)處理依據(jù),可提高故障診斷的準(zhǔn)確性。“時(shí)間窗口”設(shè)置太短易漏掉關(guān)聯(lián)信息,設(shè)置太長(zhǎng)會(huì)影響診斷效率,一般根據(jù)運(yùn)行人員的經(jīng)驗(yàn)值整定為3s～10s[11],在實(shí)際的診斷系統(tǒng)中該時(shí)間窗口值可以靈活設(shè)定。在本文中時(shí)間窗口設(shè)定為5s。

4.2診斷模型初始化

1) 故障發(fā)生的概率p(f)

通過(guò)國(guó)網(wǎng)公司對(duì)繼電保護(hù)與安全自動(dòng)裝置運(yùn)行情況分析,可知各個(gè)元件故障發(fā)生的統(tǒng)計(jì)概率,即統(tǒng)計(jì)周期內(nèi)發(fā)生故障的次數(shù)Nf與運(yùn)行的設(shè)備數(shù)量N之間的比值即為故障發(fā)生的概率,可作為故障發(fā)生的先驗(yàn)概率。不同電壓等級(jí)的變壓器發(fā)生故障的概率一般情況下不相等,如對(duì)于某地區(qū)電網(wǎng)統(tǒng)計(jì)數(shù)據(jù),如表1所示。

表1　變壓器運(yùn)行情況分析

根據(jù)診斷對(duì)象的不同,選取具體的先驗(yàn)概率初始化。

2) 故障f發(fā)生條件下,征兆s出現(xiàn)的概率p(s|f)

以變壓器保護(hù)動(dòng)作的條件概率為例,根據(jù)統(tǒng)計(jì)概率可得表2[12]。

表2　變壓器保護(hù)動(dòng)作情況統(tǒng)計(jì)

由表中數(shù)據(jù)可知,變壓器保護(hù)故障時(shí)正確動(dòng)作的概率p(s|f)=79.05%。

4.3實(shí)例驗(yàn)證

以圖3所示的三機(jī)九節(jié)點(diǎn)電力系統(tǒng)為例,說(shuō)明本文提出的故障診斷系統(tǒng)的診斷過(guò)程。

圖3　三機(jī)九節(jié)點(diǎn)系統(tǒng)圖

步驟1固定時(shí)間窗內(nèi)進(jìn)行接收警報(bào)信息。

如圖3所示系統(tǒng)故障后,在固定時(shí)間窗口內(nèi)遙信信息如表3所示。

表3　遙信信息

由表2即可得到征兆SN={L5lm(t1)=1,CB15(t2)=1,B7S(t3)=1,L5lb(t4)=1},現(xiàn)根據(jù)本文中所述方法進(jìn)行故障診斷。

步驟2保護(hù)動(dòng)作信息識(shí)別

由于L5的主保護(hù)和后備保護(hù)均指向L5,因此合并為一個(gè)保護(hù),可得S0={L5lm(t1)=1,CB15(t2)=1,B7S(t3)=1}

1) 全局完備化

根據(jù)網(wǎng)絡(luò)拓?fù)潢P(guān)系和繼電保護(hù)關(guān)聯(lián)結(jié)點(diǎn)對(duì)信息進(jìn)行完備化處理,若某一結(jié)點(diǎn)的關(guān)聯(lián)結(jié)點(diǎn)出現(xiàn)在征兆信息中,則對(duì)其狀態(tài)進(jìn)行修正。

在S0中存在兩個(gè)保護(hù),通過(guò)網(wǎng)絡(luò)拓?fù)洳樵兊玫脚c這兩個(gè)保護(hù)有關(guān)聯(lián)的斷路器分別為：L5lm(CB14,CB15),B7S=(CB7,CB14,CB13),但其狀態(tài)并沒(méi)有上傳,故對(duì)斷路器CB14,CB7,CB13的狀態(tài)進(jìn)行再次讀取,讀取結(jié)果為動(dòng)作,則將其修正為CB13=1,CB7=1。則有S1={L5lm(t1)=1,CB15(t2)=1,B7S(t3)=1,CB13=1,CB7=1}。

如與B7相關(guān)的斷路器失靈保護(hù)CB14S相關(guān)聯(lián)的結(jié)點(diǎn)斷路器CB14未動(dòng)作,CB13動(dòng)作,CB5動(dòng)作,但是未上傳CB14的狀態(tài),則應(yīng)選擇將其作為征兆信息之一,但狀態(tài)未知。故將S1修正為S2,S2={L5lm(t1)=1,CB15(t2)=1,B7S(t3)=1,CB13=1,CB7=1,CB14=[0]}。其中狀態(tài)[0]表示該信息缺失,即斷路器CB14的動(dòng)作狀態(tài)未知。

2) 局部完備化

信息上傳過(guò)程中的缺失信息有可能是故障診斷的關(guān)鍵信息,必須對(duì)這類(lèi)信息進(jìn)行修正,即征兆信息的局部完備化。貝葉斯網(wǎng)絡(luò)則是一種概率意義上的估計(jì),能夠根據(jù)已有的先驗(yàn)知識(shí)和信息,計(jì)算出缺失信息結(jié)點(diǎn)在不同狀態(tài)下發(fā)生的概率,推斷缺失信息的狀態(tài)。

如對(duì)于征兆集為S2,斷路器CB14的動(dòng)作信號(hào)缺失,則假設(shè)與其相關(guān)聯(lián)結(jié)點(diǎn)在已觀察到狀態(tài)下,CB14不同狀態(tài)的概率,將其概率最大時(shí)的狀態(tài)作為修正的狀態(tài)。如p(CB14=1|CB14s=1,L5lm=1)=0.2055,p(CB14=0|L5lm=1,CB14s=1)=0.7955,故將CB14修正為0,則最終得到的征兆集為S3={L5lm(t1)=1,CB15(t2)=1,B7S(t3)=1,CB13=1,CB7=1,CB14=0}。

步驟3貝葉斯疑似度的故障診斷模型

1) 征兆信息完備化

根據(jù)已知的先驗(yàn)概率初始化故障診斷模型如圖4所示。

2) 貝葉斯疑似度計(jì)算

由式(1)得：Bsd(fL5,SN)=0.9764,Bsd(fCB14,SN)=1。

3) 得到最優(yōu)故障假設(shè)集

由于Symptom(fT1),Symptom(fCB2)覆蓋征兆集S2,因此故障集F={fCB14,fL5}。

圖4　故障診斷模型

4) 輸出故障診斷結(jié)果

由于最優(yōu)故障假設(shè)集中,斷路器拒動(dòng)與線路故障形成助判的關(guān)系,因此,診斷結(jié)果為變壓器故障,斷路器拒動(dòng)。該診斷結(jié)果與實(shí)際運(yùn)行相符。

步驟4判斷故障元件中是否含有線路

根據(jù)步驟3的結(jié)果,其中含有L5,因此采用基于綜合阻抗的故障線路選相技術(shù)對(duì)線路的故障相別進(jìn)行識(shí)別[13]。

三機(jī)九節(jié)點(diǎn)系統(tǒng),在線路Line5,Line6兩端都配備IED裝置,可以對(duì)電壓、電流進(jìn)行采集。對(duì)于線路Line5,Line6作以下區(qū)域劃分。

圖5　區(qū)域劃分

圖6　L5接地短路故障的綜合阻抗

區(qū)域ⅠABCIm(Zcd)-310.8546.99-310.85Ph(Zcd)-69.1876.05-69.18區(qū)域ⅡABCIm(Zcd)-636.6123.46-636.61Ph(Zcd)-89.9872.26-89.98區(qū)域ⅢABCIm(Zcd)-636.61-637.67-636.61Ph(Zcd)-90.00-90.83-90.00

由圖6的仿真結(jié)果可知,通過(guò)區(qū)域劃分后,利用Matlab仿真得出Ⅰ區(qū)B相阻抗角為76.05度,阻抗虛部為46.99,Ⅱ區(qū)B相阻抗角為72.26,阻抗虛部為23.46,最后判斷為L(zhǎng)5發(fā)生B相接地短路。

5結(jié)語(yǔ)

本文對(duì)故障元件診斷模型進(jìn)行了實(shí)例驗(yàn)證,結(jié)果正確可信,相較于以往的故障診斷方法具體優(yōu)點(diǎn)有以下突出優(yōu)點(diǎn)：

1) 采用有向二分圖法可有效降低故障診斷過(guò)程運(yùn)算的復(fù)雜程度,提高診斷速度；

2) 根據(jù)保護(hù)動(dòng)作時(shí)序及貝葉斯概率對(duì)征兆信息進(jìn)行完備化處理,解決了征兆信息中存在的不確定性問(wèn)題；

3) 模型初始化概率來(lái)源于歷史統(tǒng)計(jì)資料,數(shù)據(jù)真實(shí)可靠；

4) 對(duì)于故障線路能夠進(jìn)行故障相別診斷。

隨著智能電網(wǎng)建設(shè)的不斷展開(kāi)以及新能源在電網(wǎng)中的接入,電力系統(tǒng)運(yùn)行的不穩(wěn)定因素也在逐步上升,使得越來(lái)越多新問(wèn)題不斷涌現(xiàn)并成為故障診斷領(lǐng)域亟待解決的新問(wèn)題。本文提出的故障診斷方法是對(duì)提升電網(wǎng)故障診斷能力的有益嘗試,雖然對(duì)故障診斷方法進(jìn)行了改進(jìn),仍需在工程應(yīng)用中不斷改善。

參 考 文 獻(xiàn)

[1] 文清豐.基于多信息融合的電網(wǎng)診斷技術(shù)研究[D].北京:華北大力大學(xué),2014.

WEN Qingfeng. Grid diagnosis based on information fusion technology research[D]. Beijing: University of North China Vigorously,2014.

[2] 趙冬梅.基于多信息源的電網(wǎng)故障診斷方法研究[D].北京:華北電力大學(xué),2005.

ZHAO Dongmei. Research on power grid fault diagnosis method based on multi information source[J]. North China Electric Power University,2005.

[3] 蘇廣寧,張沛超,胡炎,等.基于多源信息的電網(wǎng)故障診斷新方法[J].電力系統(tǒng)自動(dòng)化,2012,36(1):61-65.

SU Guangning, ZHANG Peichao, HU Yan, et al. A new method of power network fault diagnosis based on multi-source information[J]. The Automation of Power System,2012,36(1):61-65.

[4] 朱傳柏,郭創(chuàng)新,曹一家.基于調(diào)度綜合數(shù)據(jù)平臺(tái)的大規(guī)模電網(wǎng)分層故障診斷[J].電力系統(tǒng)自動(dòng)化,2009,33(1):51-55.

ZHU Chuanbai, GUO Chuangxin, CAO Yijia. Hierarchical fault diagnosis of large-scale power grid based on scheduling integrated data platform[J]. Automation of Electric Power Systems,2009,33(1):51-55.

[5] 楊勝春,姚建國(guó),高宗和,等.基于調(diào)度大二次系統(tǒng)的智能化電網(wǎng)調(diào)度輔助決策的研究[J].電網(wǎng)技術(shù),2006,30(增刊):176-180.

YANG Shengchun, YAO Jianguo, GAO Zonghe, et al. Based on power dispatching automation system of intelligent grid scheduling decision[J]. Power System Technology,2006,30(Suppl):176-180.

[6] Ramsay F R F, Nethercot W. Engine ignition faults, and a new method of diagnosis and location running[J]. Power Engineering, Journal of the Institution of Electrical Engineers,1946,93(36):629-635.

[7] SUN J, QIN S Y, SONG Y H. Fault diagnosis of electric power systems based on fuzzy Petri nets[J]. IEEE Transactions on Power Systems,2004,19(4):2053-2059.

[8] 周曙,王曉茹,錢(qián)清泉.基于貝葉斯網(wǎng)絡(luò)的分布式電網(wǎng)故障診斷方法[J].電網(wǎng)技術(shù),2010,34(9):76-81.ZHOU Shu, WANG Xiaoru, QIAN Qingquan. Based on Bayesian network distributed power grid fault diagnosis method[J]. Power System Technology,2010,34(9):76-81.

[9] 徐大可,汲勝昌,李彥明.變壓器繞組變形在線監(jiān)測(cè)的理論研究[J].高電壓技術(shù),2000,26(3):38-43.XU Dake, JI Shengchang, LI Yanming. Theoretical study on the online monitoring of transformer winding deformation[J]. High Voltage Technology,2000,26(3):38-43.

[10] 索南加樂(lè),劉凱,粟小華,等.輸電線路綜合阻抗縱聯(lián)保護(hù)新原理[J].電力系統(tǒng)自動(dòng)化,2008,3:36-41.

SUONAN Jiale, LIU Kai, SU Xiaohua, et al. Transmission line impedance longitudinal link protection principle[J]. Automation of Electric Power Systems,2008,3:36-41.

[11] 王汝言,吳晴,熊余,等.基于貝葉斯征兆解釋度的鏈路故障定位算法[J].計(jì)算機(jī)應(yīng)用研究,2013,30(3):712-719.

WANG Ruyan, WU Qing, XIONG Yu, et al. A link fault location algorithm based on Bayesian sign interpretation degree[J]. Computer Application Research,2013,30(3):712-719.

收稿日期:2016年1月7日,修回日期:2016年2月13日

作者簡(jiǎn)介:李郁俠,男,博士,研究方向:電力系統(tǒng)狀態(tài)監(jiān)測(cè)與故障診斷。張凱華,男,碩士,研究方向:電力系統(tǒng)分析與優(yōu)化運(yùn)行。

中圖分類(lèi)號(hào)TM41

DOI:10.3969/j.issn.1672-9722.2016.07.013

Fault Diagnosis of Multi-source Information Fusion

LI Yuxia1ZHANG Kaihua1TAN Yalan1,2LI Shi2HE Xiaojun2CHEN Bo2

(1. Xi’an University of Technology, Xi’an710048)(2. Gird State of Han Zhong, Hanzhong723099)

AbstractThe main source of fault diagnosis information is analyzed, and the fault diagnosis model based on different information is established. The fault diagnosis model is simplified. The efficiency of fault diagnosis can be significantly improved. In this method, the suspicious fault components are obtained by calculating the Bayesian suspected degree of various components and the fault line identification method based on the integrated impedance can be identified, and the timing constraint of the protection action is put forward. The validity and correctness of the model are verified by SIMULINK.

Key Wordsfault diagnosis, synthetic impedance, timing constraint