基于時間順序脈沖神經(jīng)膜系統(tǒng)的電網(wǎng)故障診斷方法

林德垠

（西華大學(xué)電氣與電子信息學(xué)院，四川 成都 610039）

0 引言

隨著電網(wǎng)規(guī)模不斷地擴大，故障的復(fù)雜性也不斷增加，利用故障數(shù)據(jù)進行電網(wǎng)故障診斷可以快速、準(zhǔn)確地找出故障區(qū)域或故障設(shè)備，對迅速恢復(fù)供電具有重要意義。近年來，在電網(wǎng)智能故障診斷的研究課題中，許多國內(nèi)外專家與學(xué)者相繼提出多種故障診斷方法，例如專家系統(tǒng)、人工神經(jīng)網(wǎng)絡(luò)、Petri網(wǎng)、解析模型和脈沖神經(jīng)膜系統(tǒng)等。目前，基于脈沖神經(jīng)膜系統(tǒng)（Spiking Neural P System，SNPS）的故障診斷方法憑借其邏輯清晰嚴(yán)密、推理過程簡單明了的優(yōu)勢受到廣泛關(guān)注。

警報信息中的時間屬性是電網(wǎng)發(fā)生故障時的一個重要特征，其包括豐富的故障相關(guān)信息，如何利用警報信息中的時間屬性對研究電網(wǎng)故障診斷具有重要意義。文獻[6]將溯因推理應(yīng)用在電網(wǎng)故障診斷中，解釋了發(fā)生故障的原因，并且還原了發(fā)生故障的過程。文獻[7]提出了一種基于時間序列相似性匹配的電力系故障診斷方法，該方法利用警報假說時間序列與實際警報時間序列之間的距離來識別故障元件與故障類型。文獻[8]利用保護裝置的時間關(guān)聯(lián)特性將時間模糊Petri網(wǎng)運用到電網(wǎng)故障診斷中。文獻[9]為了充分利用警報信息的時序特性，提出了一種時態(tài)模糊推理實數(shù)脈沖神經(jīng)網(wǎng)絡(luò)系統(tǒng)。文獻[10]將警報信息中保護裝置動作的因果關(guān)系與脈沖神經(jīng)膜系統(tǒng)結(jié)合起來，提高了故障診斷過程中處理警報信息不確定性的能力。

脈沖神經(jīng)膜系統(tǒng)是膜計算領(lǐng)域的一個重要分支，它是受生物神經(jīng)元通過脈沖形式存儲、傳遞和交換信息啟發(fā)而提出的一類特殊神經(jīng)膜系統(tǒng)。SNPS的圖形化診斷模型不僅可以直觀、形象地表達、推理故障知識，而且其強大的分布式并行計算能力能夠有效地處理故障警報信息的不確定性。然而，SNPS在電網(wǎng)故障診斷過程中處理時序信息的能力還較為薄弱，沒有充分利用警報信息的時間屬性以及信息之間的時間約束特性。因此，該文在脈沖神經(jīng)膜系統(tǒng)的基礎(chǔ)上，利用警報信息中斷路器、繼電器動作的時間屬性以及保護裝置之間的動作順序和時序約束關(guān)系建立基于時間順序?qū)崝?shù)脈沖神經(jīng)膜系統(tǒng)（Time Sequence Spiking Neural P System with Real Numbers，rTSSNPS）的故障診斷模型。同時，為了簡化對警報信息進行時序處理的過程，該文結(jié)合故障診斷模型重新設(shè)置推理算法，使算法具有處理警報信息時間屬性的能力，從而提升脈沖神經(jīng)膜系統(tǒng)處理保護裝置拒動、誤動、時標(biāo)錯誤和信息丟失等不確定性的能力，進而減少運算步驟，提高故障診斷效率。

1 rTSSNPS

1.1 rTSSNPS定義

可以將1個rTSSNPS定義為多元組，如公式（1）所示。

σ為命題神經(jīng)元集合，σ=（θ,T,λ,r,w）（1）（為命題神經(jīng)元的個數(shù)）；σ為規(guī)則神經(jīng)元集合，σ=（θ,c,T,λ,r,w）（1）（為規(guī)則神經(jīng)元的個數(shù)）。其中，θ為第個命題神經(jīng)元內(nèi)部脈沖的脈沖值，θ為[0，1]的實數(shù)；θ為第個規(guī)則神經(jīng)元內(nèi)部脈沖的脈沖值，θ為[0，1]的實數(shù)；T為命題神經(jīng)元的時間信息脈沖，每個命題神經(jīng)元都會產(chǎn)生1個時間信息；T為各個命題神經(jīng)元之間的時間距離，表示上個命題神經(jīng)元點火到下個命題神經(jīng)元激發(fā)的時間距離約束；r為命題神經(jīng)元σ的點火規(guī)則，其形式為/α→α（為點火條件；為神經(jīng)元；為傳遞的脈沖），和為[0，1]的實數(shù)；r為規(guī)則神經(jīng)元σ的點火規(guī)則，其形式為/α→α（為產(chǎn)生的新脈沖），和為[0，1]的實數(shù)；λ和λ分別為命題神經(jīng)元和規(guī)則神經(jīng)元的點火閥值，它們都是[0，1]的實數(shù)；w和w分別為命題神經(jīng)元σ的輸出權(quán)重向量和規(guī)則神經(jīng)元σ的輸出權(quán)重向量；c為規(guī)則神經(jīng)元的真值，c為[0，1]的實數(shù)，與σ對應(yīng)的模糊產(chǎn)生式規(guī)則的確定性因子相等。該脈沖神經(jīng)膜系統(tǒng)涉及4種類型的神經(jīng)元，即一類命題神經(jīng)元（用符號“”表示）和三類規(guī)則神經(jīng)元，即“常規(guī)”“與”和“或”規(guī)則神經(jīng)元，分別用字母“general”“and”和“or”來表示。

1.2 模糊產(chǎn)生式規(guī)則

類型一（常規(guī)規(guī)則）為（-c）：IFp（θ，T，T，w）THENp（θ，T），相應(yīng)的rTSSNPS模型如圖1（a）所示。其中，p與p分別為規(guī)則R的前件命題與后件命題；c為規(guī)則R的確定因子，c為[0，1]的實數(shù)；θ和θ分別為命題p和p產(chǎn)生的真值，θ和θ為[0，1]的實數(shù)；T、T分別為完成命題p和命題p所消耗的時間，也可以表示命題神經(jīng)元σ和命題神經(jīng)元σ點火的時間點；T為命題p到命題p所消耗的時間，也可以表示命題神經(jīng)元σ到命題神經(jīng)元σ的時間距離；w為命題p的權(quán)重值。

命題p的真值和時間如公式（2）~公式（3）所示。

其運算符號的解釋如公式（4）所示。

式中：為（）維度的邏輯矩陣；、為的時間矩陣；和分別為時間矩陣和中的時間元素；A為矩陣中的邏輯元素；、為該矩陣的行和列。

類型二 （復(fù)合“與”規(guī)則）為（=c）：IF（，，T，）and…andp（θ，T，T，w）THENp（θ，T），相應(yīng)的rTSSNPS如圖1（b）所示。其中，，…，p為規(guī)則R的前件命題；p為規(guī)則R的后件命題；c為規(guī)則R的確定性因子，c為[0，1]的實數(shù)；，…，θ分別為，…，p的真值，，…，θ為[0，1]的實數(shù)；θ為命題p的真值，θ為[0，1]的實數(shù)；，…，T分別為完成命題，…，p所消耗的時間，也可以分別表示命題神經(jīng)元，…，σ點火的時間點；T=[T，…，T]分別為命題，…，p到命題p所消耗的時間，也可以分別表示命題神經(jīng)元，…，σ到命題神經(jīng)元σ的時間距離；，…，w分別為命題，…，p的權(quán)重值。

命題p的真值和時間如公式（5）~公式（6）所示。

類型三（復(fù)合“或”）為（-c）：IF（，，T，）or …orp（θ，T，T，w）THENp（θ，T），相應(yīng)的rTSSNPS如圖1（c）所示。其中，，…，p為規(guī)則R的前件命題；p為規(guī)則R的后件命題；c為規(guī)則R的確定因子，c為[0，1]的實數(shù)；，…，θ分別為命題，…，p的真值，，…，θ為[0，1]的實數(shù)；θ為命題p的真值，θ為[0，1]的實數(shù)；，…，T分別為完成命題，…，p所消耗的時間，也可以分別表示命題神經(jīng)元，…，σ點火的時間點；T=[T，…，T]分別為命題，…，p到命題p所消耗的時間，也可以分別表示命題神經(jīng)元，…，σ到命題神經(jīng)元σ的時間距離；，…，w分別為命題，…，p的權(quán)重值。

圖1 rTSSNPS 3種模糊產(chǎn)生式規(guī)則

命題p的真值和時間如公式（7）~公式（8）所示。

其運算符號的解釋如公式（9）所示。

式中：、為×的實數(shù)矩陣；、為×的時間矩陣；為×的邏輯矩陣；max為最大值，、為該矩陣的行和列。

1.3 時間脈沖矩陣推理算法

為了能使算法準(zhǔn)確、高效地表達和推力知識，減少運算的復(fù)雜程度，更加契合故障診斷過程中保護裝置動作的邏輯關(guān)系，該文根據(jù)線路和母線的故障診斷模型將推理算法分為算法一和算法二。其中，算法一完成線路送端與受端雙向并行矩陣推理運算以及母線各個出線方向多向并行矩陣推理運算任務(wù)；算法二綜合線路送端與受端、母線各個出線方向輸出結(jié)果對故障元件進行識別。時序脈沖矩陣推理算法如下。

算法一如下。

輸入：W，W，T，，λ，λ，，，，

輸出：神經(jīng)元的時間與真值

令推理步驟=0

根據(jù)輸入的初始條件計算根據(jù)、、、計算

則推理步數(shù)=+1

While （<3）do

for 每個命題神經(jīng)元 do

if 點火規(guī)則={α，θ≥λ，1}滿足，then該命題神經(jīng)元點火并根據(jù)公式

計算和

if 命題神經(jīng)元存在突觸后規(guī)則神經(jīng)元 then 該命題神經(jīng)元向其突觸后規(guī)則神經(jīng)元發(fā)送1個脈沖和時間

else 該命題神經(jīng)元僅積累脈沖和時間

for 每個規(guī)則神經(jīng)元do

if 點火規(guī)則={α，δ≥λ，1}滿足，

then 該神經(jīng)元點火并根據(jù)公式

計算θ與，并向其突觸后命題神經(jīng)元發(fā)送1個脈沖與時間

其中，θ=[,…,θ]為命題神經(jīng)元的實數(shù)脈沖值矩陣θ（=1,…,） 為[0，1]的實數(shù)，表示第個命題神經(jīng)元的脈沖值；δ=[,…,δ]為規(guī)則神經(jīng)元的實數(shù)脈沖值矩陣，δ（=1,…,）為[0，1]的實數(shù)，表示第個規(guī)則神經(jīng)元的脈沖值；=（,…,c）為規(guī)則神經(jīng)元的對角矩陣，c（=1,…,） 為[0，1]的實數(shù)，表示第個規(guī)則神經(jīng)元的真值；T=[,…,t]為當(dāng)前所在推理步驟命題神經(jīng)元的時間矩陣，該矩陣包括個命題神經(jīng)元的動作時間，t（=1,…,）為時間信息，表示第個命題神經(jīng)元點火的時間；T=[t,…,t]為規(guī)則神經(jīng)元的時間距離矩陣，t（=1,…,）時間常量，表示規(guī)則神經(jīng)元前件命題神經(jīng)元的點火到規(guī)則神經(jīng)元后件命題神經(jīng)元激發(fā)時消耗的時間；W為規(guī)則神經(jīng)元到命題神經(jīng)元的突觸權(quán)重矩陣；W、W、W分別為命題神經(jīng)元到“常規(guī)”規(guī)則神經(jīng)元、“與”規(guī)則神經(jīng)元以及“或”規(guī)則神經(jīng)元的突觸權(quán)重矩陣；λ=[λ,…,λ]為命題神經(jīng)元的點火閥值矩陣，λ（=1,…,） 為[0，1]的實數(shù)，表示第個命題神經(jīng)元的點火閥值）；λ=[λ,…,λ]為規(guī)則神經(jīng)元的點火閥值矩陣，λ（=1,…,）為[0，1]的實數(shù)，表示第個命題神經(jīng)元的點火閥值。

算法二如下。

輸入：W,W,T,,,

輸出：神經(jīng)元時間與真值

綜合算法一的輸出結(jié)果令推理步驟=3，計算

輸出：θ=[θ+（T）]/2和T

如果給出rTSSNPS中輸入神經(jīng)元的初始脈沖值、時間點約束矩陣以及各個命題神經(jīng)元之間的時間距離約束矩陣，那么該系統(tǒng)通過執(zhí)行算法一和算法二就能夠獲得其余命題神經(jīng)元的真值與時間點約束區(qū)間。θ=[，…，θ]為命題神經(jīng)元的實數(shù)脈沖值矩陣，θ（=1，…，）為[0，1]的實數(shù)；δ=[，…，δ]為規(guī)則神經(jīng)元的實數(shù)脈沖值矩陣；δ（=1，…，）為第個規(guī)則神經(jīng)元的脈沖值，δ（=1，…，）為[0，1]的實數(shù)；=diag（，…，c）為規(guī)則神經(jīng)元的對角矩陣；c（=1，…，）為第個規(guī)則神經(jīng)元的真值，c（=1，…，）為[0，1]的實數(shù)；T=[，…，t]為當(dāng)前所在推理步驟命題神經(jīng)元的時間矩陣，該矩陣包括個命題神經(jīng)元的動作時間；t（=1，…，）為時間信息（第個命題神經(jīng)元點火的時間）；T=[t，…，t]T為規(guī)則神經(jīng)元的時間距離矩陣；t（=1，…，）為時間常量（規(guī)則神經(jīng)元前件命題神經(jīng)元點火到規(guī)則神經(jīng)元后件命題神經(jīng)元激發(fā)時消耗的時間）；W為規(guī)則神經(jīng)元到命題神經(jīng)元的突觸權(quán)重矩陣（如果規(guī)則神經(jīng)元σ到命題神經(jīng)元σ存在突觸，則w∈（0，1]，其等于突觸（，）的輸出權(quán)重，否則w=0；W，W，W分別為命題神經(jīng)元到“常規(guī)”規(guī)則神經(jīng)元、“與”規(guī)則神經(jīng)元以及“或”規(guī)則神經(jīng)元的突觸權(quán)重矩陣（如果命題神經(jīng)元σ到規(guī)則神經(jīng)元σ存在突觸，則w∈（0，1]，其等于突觸（，）的輸出權(quán)重，否則w=0；λ=（λ，…，λ）為命題神經(jīng)元的點火閥值矩陣；λ（=1，…，）為第個命題神經(jīng)元的點火閥值，λ（=1，…，）為[0，1]的實數(shù)；λ=（λ，…，λ）為規(guī)則神經(jīng)元的點火閥值矩陣；λ（=1，…，）為第個命題神經(jīng)元的點火閥值，λ（=1，…，）為[0，1]的實數(shù)。

在該算法中，各個乘法算子定義如公式（10）所示。

式中：為×的實數(shù)矩陣；、為×的時間矩陣；為一種運算推理過程。

對故障診斷模型來說，已知原因事件實際動作時間（t）及通過時間推理得到的原因事件的動作時間（t），利用公式（10）判斷推理過程是否準(zhǔn)確，如公式（11）所示。

式中：為×的實數(shù)矩陣；為×的時間矩陣；為一種運算推理過程。

對故障診斷模型來說，利用公式（11）檢查事件經(jīng)過時序信息處理是否因故障而動作，如公式（12）~公式（13）所示。

式中：、為×的實數(shù)矩陣；為×的邏輯矩陣；avg為數(shù)值取平均數(shù)。

式中：為×的實數(shù)矩陣；為×的邏輯矩陣；plus為數(shù)值相加。

2 基于rTSSNPS的電網(wǎng)故障診斷模型

該文簡化利用時序信息的過程，將保護裝置動作的時間屬性融入故障診斷模型中，使神經(jīng)元不僅能夠傳遞脈沖值，而且能夠傳遞時間信號；同時，通過建立診斷模型增強拓?fù)浣Y(jié)構(gòu)的適應(yīng)能力，降低計算的復(fù)雜程度，提高了故障診斷的效率。

2.1 利用時序信息

當(dāng)電網(wǎng)發(fā)生故障時，保護裝置動作有嚴(yán)密的邏輯關(guān)系，以線路或母線故障為起點，保護裝置動作的先后順序如下：主保護繼電器及其對應(yīng)的斷路器動作，近后備保護繼電器及其對應(yīng)斷路器動作，遠后備保護繼電器及其對應(yīng)斷路器動作。主保護、近后備保護和遠后備保護動作的時間約束分別為10 ms~40 ms、310 ms~340 ms和450 ms~510 ms，斷路器到各相應(yīng)保護繼電器的時間約束為20 ms~60 ms。

結(jié)合rTSSNPS故障診斷模型，為了簡化處理警報信息中時間信息的過程，根據(jù)文獻[11]并利用t及其對應(yīng)時間距離約束（t，t）得到該事件的前驅(qū)事件以及前驅(qū)事件的時間約束，如公式（15）所示。式中：（t）為該事件的前驅(qū)事件的時間約束；（t）為該事件的時間約束。

2.2 故障診斷模型

以IEEE39節(jié)點系統(tǒng)為例，建立線路L診斷模型（圖2）。該文根據(jù)保護裝置動作時間先后的邏輯關(guān)系建立故障診斷模型（當(dāng)診斷電網(wǎng)故障時，根據(jù)從元件故障到繼電器動作再到斷路器跳閘的動作順序和時間約束建立診斷模型），每個命題神經(jīng)元都會被賦予1個脈沖值和1個時間點約束。其中，脈沖值表征該命題神經(jīng)元的置信度，時間點約束表征該命題神經(jīng)元點火的時間約束。此外，命題神經(jīng)元之間都會有1個時間區(qū)間約束，表征從上一個命題神經(jīng)元點火到下一個命題神經(jīng)元激發(fā)所消耗的時間。

圖2 線路L0414的rTSSNPS的故障診斷模型

為了更直接、更清晰地將警報信息展示在建立的診斷模型中，用~表示警報信息中保護裝置的動作時間，用（t，t）、（t，t）等表示保護裝置之間的時間約束。與文獻[10]不同，該文提出的故障診斷模型直接將警報信息的時間點約束矩陣、保護裝置之間的時間距離約束矩陣和保護裝置的初始置信度作為輸入，將元件故障概率和故障約束時間作為輸出，簡化處理警報信息中時間信息的過程，從而提高故障診斷的效率。

同時，該文為了使診斷模型適應(yīng)電網(wǎng)拓?fù)浣Y(jié)構(gòu)的變化，根據(jù)動作時間的先后邏輯建立故障元件的診斷模型。當(dāng)線路發(fā)生故障時，在線路的送端與受端分別建立主保護子網(wǎng)、近后備保護子網(wǎng)以及遠后備保護子網(wǎng)；當(dāng)母線發(fā)生故障時，在母線的出線方向分別建立主保護子網(wǎng)、遠后備保護子網(wǎng)。當(dāng)拓?fù)浣Y(jié)構(gòu)發(fā)生變化時，線路故障診斷模型只需要增加或減少遠后備保護及其對應(yīng)的斷路器，母線故障診斷模型只需要增加或減少出線方向的診斷子網(wǎng)絡(luò)， 保證故障診斷模型拓?fù)浣Y(jié)構(gòu)、突觸連接矩陣、輸入矩陣及輸出矩陣的變動較小。

2.3 診斷方法流程圖

基于rTSSNPS的故障診斷方法的具體操作流程，如圖3所示。

圖3 故障診斷流程圖

步驟一：確定故障區(qū)域中的可疑元件集={，，…x}。在SCADA系統(tǒng)中讀取有關(guān)保護裝置的動作信息，利用結(jié)線分析法判斷停電區(qū)域中的可疑故障元件集。

步驟二：建立故障診斷模型。利用保護裝置邏輯動作規(guī)則和警報信息的時間信息分別針對各個可疑故障元件x建立基于rTSSNPS的故障診斷模型。

步驟三：獲取輸入神經(jīng)元脈沖值和時間矩陣。根據(jù)表1 和表2的保護裝置動作的可信度與保護裝置為未動作的可信度，并結(jié)合實際動作信息賦予輸入神經(jīng)元脈沖值，利用警報信息的時間信息賦予并確定輸入神經(jīng)元的時間矩陣，再根據(jù)保護動作邏輯規(guī)則確定各命題神經(jīng)元之間的時間距離矩陣。

表1 保護裝置的動作可信度值

表2 保護裝置的未動作可信度值

步驟四：判斷故障元件。對故障區(qū)域中的各個可疑元件x執(zhí)行時間脈沖矩陣推理算法，并計算輸出神經(jīng)元的時間點約束區(qū)間及其脈沖值。其中，設(shè)置輸出神經(jīng)元的故障閥值為0.5，判定數(shù)值大于0.5的元件為故障元件，否則無故障元件。

步驟五：輸出診斷結(jié)果。將步驟四推理得到故障元件及相關(guān)時間區(qū)間與警報信息進行對比，判斷保護裝置是否存在拒動、誤動以及漏報等情況。

3 算例分析

基于IEEE39節(jié)點系統(tǒng)模擬不同場景下的電網(wǎng)故障診斷過程，以驗證該文提出的模型的有效性和容錯性。這些場景包括狀態(tài)信息不完整和不確定的復(fù)雜故障，警報信息見表3。將該文提出的方法與文獻[10]、文獻[12]中的2種經(jīng)典故障診斷方法進行對比，結(jié)果見表4。同時，將該方法與近年來提出的故障診斷方法進行對比分析，結(jié)果見表5。

表3 保護裝置故障信息

表4 算例識別結(jié)果與比較

表5 該文所采用的方法與現(xiàn)有3種方法的對比

表4中算例一~算例二為單一故障，算例三~算例六為多重故障。對算例一、算例三、算例五和算例六來說，3種方法都能準(zhǔn)確識別故障元件，然而文獻[12]未能明確不同故障元件的時間點約束與時間區(qū)間約束的關(guān)聯(lián)屬性，其在算例二中錯誤診斷母線B發(fā)生故障，在算例四中錯誤診斷L發(fā)生故障。從表5中可以看出，文獻[3]不能識別時序錯誤的動作信息，而文獻[10]與文獻[12]雖然能夠正確地對保護裝置進行評價，但是由于2種方法都需要先對警報信息進行時序處理，因此診斷速度較慢。

對表4與表5進行分析可知，該文提出的故障診斷方法，充分利用警報信息中保護裝置的時序信息以及保護裝置之間的時序約束關(guān)系，并將其在建立的診斷模型上直觀地表現(xiàn)出來，實現(xiàn)了時序信息與故障置信度的并行推理，能夠快速處理繼電器/斷路器拒動、誤動以及警報信息缺失等復(fù)雜情況，進而準(zhǔn)確識別故障元件，提高故障診斷的效率，具有較高的容錯性（以算例二和算例三為例進行詳細(xì)說明）。

算例二：單重故障。在發(fā)生故障后，調(diào)度中心收到的警告信息為B（31 ms）、L（27 ms）、L（32 ms）、L（472 ms）、L（466 ms）、CB（73 ms）、CB（504 ms）以及CB（516 ms）。采用基于rTSSNPS的故障識別方法對該案例進行故障診斷的流程如下。

步驟一：確定故障區(qū)域可疑元件集。根據(jù)警報信息利用結(jié)線分析法搜索故障區(qū)域得到可疑故障元件集={B，L，B}。

步驟二：建立故障診斷模型。根據(jù)可疑故障元件集={B，L，B}分別針對線路L、母線B和母線B建立基于rTSSNPS的故障診斷模型（由于篇幅有限，算例二以線路L為例進行詳細(xì)說明，其故障診斷模型如圖2所示）。

步驟三：獲取輸入神經(jīng)元脈沖值和時間矩陣。該文以線路L送端診斷模型為例進行推理說明，輸入神經(jīng)元脈沖值和時間矩陣分別如公式（16）~公式（20）所示。

式中：為零矩陣（不同下角標(biāo)對應(yīng)不同的運算規(guī)則）；?為警報信息未包括該保護裝置動作的時間信息。

相關(guān)突出權(quán)重矩陣如公式（21）所示。

步驟四：判斷故障元件。線路L送端診斷模型執(zhí)行時間脈沖矩陣推理算法，當(dāng)=0時，推理得到的結(jié)果如公式（22）所示。

當(dāng)=1時，推理得到的結(jié)果如公式（23）所示。

當(dāng)=2時，推理得到的結(jié)果如公式（24）所示。

同理，可以得到其受端診斷模型的輸出脈沖值與時間點約束分別為0.694 7和-38 ms~16 ms，則執(zhí)行時序脈沖矩陣推理算法二得到輸出結(jié)果為0.897 2和-13 ms~16 ms。因此，判斷線路L為故障元件，其故障概率為0.8972，故障時刻為-13 ms~16 ms。類似的利用相關(guān)報警信息構(gòu)建母線B和B的rTSSNPS故障診斷模型，2條母線故障概率分別為0.267 4和0.100 4。

步驟五：輸出診斷結(jié)果。根據(jù)步驟四對故障元件的判斷并結(jié)合SCADA系統(tǒng)中獲得的警報信息可以推出線路L在-13 ms~16 ms時刻發(fā)生故障。其中，在線路L發(fā)生故障后，線路兩端主保護正確動作，送端斷路器CB跳閘，但是線路受端CB拒動，由線路受端遠后備保護L、L動作跳開兩端斷路器CB、CB，而B誤動，B無故障。

算例三：多重故障。故障發(fā)生后調(diào)度中心收到的警告信息為B（25 ms）、L（27 ms）、L（32 ms）、L（475 ms）、CB（52 ms）、CB（65 ms）、CB（65 ms）、CB（65 ms）以及CB（521 ms）。

算例三的診斷過程與算例二相同，首先，確定故障區(qū)域可疑元件集={L，B，L，L}。其次，分別建立可疑元件的故障診斷模型，并根據(jù)警報信息獲取輸入神經(jīng)元脈沖值與時間矩陣，根據(jù)矩陣推理算法得到可疑故障元件的故障置信度和時間點約束分別為（L，0.9557，（-8 ms~17 ms））、（B，0.9028，（-11 ms~15 ms））、（L，0.1200，?）以及（L，0.1343，?）。最后，根據(jù)故障元件的判斷以及相關(guān)警報信息可知，線路L在-8 ms~17 ms時間段發(fā)生故障，線路兩端主保護正確動作并作用于相應(yīng)斷路器跳閘；母線B在-11 ms~15 ms時間段發(fā)生故障，母線B的L和L出線方向主保護正確動作跳開相應(yīng)的斷路器，但是母線B的L出線方向斷路器CB拒動，其遠后備保護L動作跳開斷路器CB，線路L、L無故障。

4 結(jié)論

該文提出了一種改進的時間順序?qū)崝?shù)脈沖神經(jīng)膜系統(tǒng)電網(wǎng)故障診斷模型和算法，根據(jù)警報信息的時間屬性與保護裝置之間的時序約束邏輯關(guān)系建立故障診斷模型， 減少其拓?fù)浣Y(jié)構(gòu)隨著網(wǎng)絡(luò)變化的調(diào)整幅度；將繼電器和斷路器的初始置信度和其動作時間作為時序脈沖矩陣推理算法的輸入，將元件故障概率和元件故障約束時間區(qū)間作為輸出，省去保護裝置時序約束一致性判斷步驟；將時間脈沖矩陣推理算法與rTTSNPS故障診斷模型相結(jié)合，實現(xiàn)并行推理運算，減少了計算的復(fù)雜程度。該方法充分運用保護裝置動作的時間邏輯關(guān)系，簡化時序推理過程，拓?fù)浣Y(jié)構(gòu)適應(yīng)性強， 可以提高故障診斷效率，還可以處理警報信息缺失、時標(biāo)錯誤，保護裝置拒動誤動等問題。