基于EEMD樣本熵的直流牽引網(wǎng)故障電流識(shí)別

冷 月，楊洪耕，王智琦

（四川大學(xué)電氣信息學(xué)院，四川 成都 610065）

基于EEMD樣本熵的直流牽引網(wǎng)故障電流識(shí)別

冷 月，楊洪耕，王智琦

（四川大學(xué)電氣信息學(xué)院，四川 成都 610065）

針對(duì)地鐵直流牽引網(wǎng)的振蕩電流容易引起繼電保護(hù)系統(tǒng)頻繁誤動(dòng)的問(wèn)題，提出一種結(jié)合總體平均經(jīng)驗(yàn)?zāi)B(tài)分解(ensemble empirical mode decomposition，EEMD)和樣本熵的直流牽引網(wǎng)振蕩電流與短路故障電流識(shí)別方法。利用EEMD方法對(duì)直流牽引網(wǎng)的饋線電流信號(hào)進(jìn)行分解，求取各固有模態(tài)函數(shù)（intrinsic mode function，IMF）分量的樣本熵值，并將計(jì)算結(jié)果求和，進(jìn)而獲得反映直流牽引網(wǎng)運(yùn)行狀態(tài)信息的特征量。通過(guò)對(duì)典型饋線電流信號(hào)進(jìn)行分析計(jì)算，可知EEMD和樣本熵相結(jié)合的特征提取方法可以有效地區(qū)分直流牽引網(wǎng)振蕩電流與短路故障電流。算例分析驗(yàn)證該方法的有效性。

短路故障電流；振蕩電流；樣本熵；總體平均經(jīng)驗(yàn)?zāi)B(tài)分解；直流牽引網(wǎng)

0 引 言

近年來(lái)，隨著我國(guó)各大城市交通擁擠問(wèn)題越來(lái)越嚴(yán)重，城市軌道交通為多數(shù)城市所青睞。然而，由于城市軌道交通運(yùn)載能力提升、行車密度增加及多車啟動(dòng)等原因，致使直流牽引網(wǎng)出現(xiàn)振蕩電流。

目前，城市軌道交通直流牽引網(wǎng)的饋線主保護(hù)采用di/dt-ΔI（電流變化率-電流增量）保護(hù)，它能很好地區(qū)分短路電流和啟動(dòng)電流[1]。而牽引網(wǎng)振蕩電流的增量和上升率特征均接近于短路電流[2]，致使di/dt-ΔI保護(hù)不能有效地識(shí)別這兩種電流，容易造成饋線保護(hù)頻繁誤動(dòng)作。因此，要改善直流牽引網(wǎng)饋線保護(hù)頻繁誤動(dòng)問(wèn)題，需在研究牽引網(wǎng)振蕩電流和短路故障電流波形特點(diǎn)的基礎(chǔ)上，尋求新的特征提取方法。文獻(xiàn)[3]利用小波包變換計(jì)算牽引網(wǎng)電流的能量譜實(shí)現(xiàn)故障診斷；文獻(xiàn)[4]用經(jīng)驗(yàn)?zāi)B(tài)分解（empirical mode decomposition，EMD），對(duì)牽引網(wǎng)饋線電流進(jìn)行處理，將分解后的余量斜率作為特征量對(duì)牽引網(wǎng)運(yùn)行狀態(tài)進(jìn)行判斷；文獻(xiàn)[5]采用EMD算法對(duì)直流牽引網(wǎng)電流信號(hào)進(jìn)行分解，構(gòu)建EMD多尺度能量熵，并將其作為特征量對(duì)牽引網(wǎng)振蕩電流和故障電流進(jìn)行識(shí)別。文獻(xiàn)[6]對(duì)直流牽引網(wǎng)電流信號(hào)進(jìn)行相空間重構(gòu)和關(guān)聯(lián)維數(shù)計(jì)算，并將關(guān)聯(lián)維數(shù)作為特征量進(jìn)行故障識(shí)別。

本文以城市軌道交通直流牽引網(wǎng)饋線電流為研究對(duì)象，通過(guò)分析直流牽引網(wǎng)短路故障電流和振蕩電流的波形特征，結(jié)合振蕩電流和短路故障電流非線性的特點(diǎn)，針對(duì)di/dt-ΔI保護(hù)在識(shí)別直流牽引網(wǎng)振蕩電流方面存在的不足，提出一種基于EEMD樣本熵的直流牽引網(wǎng)振蕩電流與故障電流識(shí)別方法。將饋線電流信號(hào)通過(guò)EEMD方法分解為有限個(gè)IMF分量，分別計(jì)算其樣本熵值，將計(jì)算結(jié)果求和得到特征量，并用它來(lái)定量刻畫直流牽引網(wǎng)的運(yùn)行狀態(tài)特征。通過(guò)算例分析，驗(yàn)證了該方法能敏感地從波形信號(hào)中提取反映直流牽引網(wǎng)運(yùn)行狀態(tài)信息的特征量，準(zhǔn)確識(shí)別振蕩電流和短路故障電流。

1 直流牽引網(wǎng)短路故障電流和振蕩電流

由文獻(xiàn)[2，5-6]可知，牽引網(wǎng)短路故障電流的波形特征如下：

1）故障發(fā)生后，電流按指數(shù)規(guī)律快速上升至峰值；

2）電流上升前不存在振蕩；

3）電流波形比較粗糙，且無(wú)明顯的規(guī)則性。

牽引網(wǎng)振蕩電流的波形特征如下：

1）屬于典型的非線性非平穩(wěn)信號(hào)；

2）電流快速上升前存在小幅度振蕩，且振蕩周期約為100ms；

3）所包含的頻率分量少于短路故障電流；

4）隨著機(jī)車移動(dòng)或擾動(dòng)結(jié)束，振蕩消失，因此不能將其判為故障電流。

直流牽引網(wǎng)振蕩電流在幅值增量和變化率上與短路故障電流差異不大，但是振蕩電流波形更具規(guī)則性[7]。

2 EEMD及樣本熵理論

2.1 EEMD理論

EMD是一種處理非平穩(wěn)信號(hào)的新方法[8-10]，它通過(guò)“篩分”的方式將復(fù)雜信號(hào)按不同時(shí)間尺度自適應(yīng)地分解成一系列IMF分量和一個(gè)殘余分量之和。雖然EMD方法已被成功應(yīng)用于諸多領(lǐng)域，但仍存在模態(tài)混疊等缺點(diǎn)[8，11-13]。為此，Huang和Wu提出了EEMD方法，其本質(zhì)是加高斯白噪聲的多次EMD[14]。

EEMD的具體分解步驟[14]如下：

1）在原始信號(hào)x（t）中加入多組隨機(jī)高斯白噪聲信號(hào)ni（t），得到加噪后的總體信號(hào)為

式中：xi（t）——第i次加噪后的總體信號(hào)；

k——加入白噪聲的幅值系數(shù)。

2）將信號(hào)xi（t）進(jìn)行EMD分解[9-10]得

式中cij（t）為第i次加噪后分解得到的第j個(gè)IMF分量。

3）得到最終的IMF分量cj（t）和殘余分量raν（t）：

式中M為加入的高斯白噪聲信號(hào)的組數(shù)，取100。

2.2 樣本熵理論

近似熵[15]是一種統(tǒng)計(jì)量化信號(hào)復(fù)雜程度的無(wú)量綱指標(biāo)，只需要較少的數(shù)據(jù)就可以得到合理穩(wěn)健的估計(jì)值。但由于近似熵存在固有的對(duì)自身數(shù)據(jù)段的比較過(guò)程，所以計(jì)算一定會(huì)產(chǎn)生偏差。針對(duì)近似熵的缺陷，Richman提出了樣本熵[16-17]算法。

設(shè)長(zhǎng)度為N的數(shù)據(jù)序列為{x（i），i=1，2，…，N}，則可通過(guò)以下步驟[16-17]計(jì)算樣本熵：1）用原數(shù)據(jù)序列{x（i）}組成一組m維向量X（i）：

式中i=1，2，…，N-m+1；m為模式維數(shù)。

2）計(jì)算X（i）與X（j）間的距離d（i，j）：

3）按照設(shè)定的相似容限r(nóng)，對(duì)每個(gè)i值進(jìn)行統(tǒng)計(jì)，將滿足條件d（i，j）＜r的數(shù)目記為n。計(jì)算n與向量總數(shù)N-m+1的比值

5）m增加1，重復(fù)上述步驟得Cm+1（r）。

6）該數(shù)據(jù)序列樣本熵的理論值為

而實(shí)際應(yīng)用中

一般情況下，取m=2，r=（0.1～0.25）σ（σ為原始數(shù)據(jù)序列的標(biāo)準(zhǔn)差）。

下面采用4組分析信號(hào)來(lái)驗(yàn)證樣本熵算法在特征提取中所具有的能力。定義如圖1所示的4組理想信號(hào)：x1（t）=sin（2πf1t）；x2（t）=0.5sin（2πf2t）；x3（t）= x1（t）+x2（t）；x4（t）=x1（t）+x2（t）+N。其中：f1=50 Hz，f2= 150Hz，N為均勻白噪聲。本文取m=2，r=0.2σ。

圖1 4組分析信號(hào)波形

表1為計(jì)算所得信號(hào)x1（t）、x2（t）、x3（t）和x4（t）的樣本熵值。

通過(guò)表1可以看出，x1（t）為規(guī)則的基波正弦，所以其樣本熵值較小；x2（t）也為基波正弦信號(hào)，故其樣本熵值與x1（t）相近；而x3（t）由x1（t）和x2（t）疊加而成，信號(hào)變復(fù)雜，故其樣本熵值也相應(yīng)增大；x4（t）在x3（t）的基礎(chǔ)上又疊加了白噪聲，故其近似熵值繼續(xù)增大。由此可見(jiàn)，樣本熵值大小與信號(hào)復(fù)雜程度有關(guān)，與理論分析一致。

表1 4組分析信號(hào)的樣本熵

3 基于EEMD樣本熵的直流牽引網(wǎng)振蕩電流與故障電流識(shí)別方法

基于EEMD樣本熵的直流牽引網(wǎng)振蕩電流與故障電流識(shí)別步驟為：1）采樣直流牽引網(wǎng)饋線電流信號(hào)；2）對(duì)采樣信號(hào)進(jìn)行EEMD分解；3）求取各IMF分量的樣本熵值，并求和得到特征量E；4）根據(jù)E的大小識(shí)別振蕩電流與短路故障電流。

3.1 短路故障電流

利用圖像數(shù)字化軟件得到文獻(xiàn)[5]中圖2（a）保護(hù)啟動(dòng)前200 ms的短路故障電流數(shù)據(jù)作為分析對(duì)象，并利用Matlab繪制曲線，如圖2所示。

圖2 保護(hù)啟動(dòng)前的短路故障電流波形

圖2所示信號(hào)的EEMD分解結(jié)果，如圖3所示。添加100組幅值為0.05的高斯白噪聲信號(hào)到原始信號(hào)中，從圖3可以看出，各個(gè)IMF分量使得短路故障電流信號(hào)的波形特征在不同的時(shí)間尺度下顯現(xiàn)出來(lái)。

圖3 短路故障電流的EEMD分解

計(jì)算經(jīng)EEMD分解得到的各IMF分量相對(duì)應(yīng)的樣本熵值，并將各樣本熵值求和，可以得到短路故障電流的特征量E為1.942。

3.2 振蕩電流

利用圖像數(shù)字化軟件得到文獻(xiàn)[5]中圖3保護(hù)啟動(dòng)前200ms的振蕩電流數(shù)據(jù)作為分析對(duì)象，并利用Matlab繪制曲線，如圖4所示。

同樣，圖4對(duì)所示信號(hào)加入100組白噪聲信號(hào)，每組白噪聲信號(hào)的幅值設(shè)為0.05，進(jìn)行EEMD分解。圖5中，第1行至第5行分別為IMF分量c1~c5，反映了從高頻到低頻不同時(shí)間尺度的振蕩電流信號(hào)的波形特征；第6行為殘余分量r。

圖4 保護(hù)啟動(dòng)前的振蕩電流波形

圖5 振蕩電流的EEMD分解

同理，分別計(jì)算IMF分量c1~c5各自的樣本熵值，并求和得到振蕩電流的特征量E為0.169。

4 算例分析

為了進(jìn)一步驗(yàn)證本文方法的有效性，本文參考文獻(xiàn)[18-21]，利用Matlab/Simulink進(jìn)行實(shí)驗(yàn)，得到多例電流信號(hào)。用兩種方法（方法1為EMD多尺度能量熵，方法2為本文方法），計(jì)算電流信號(hào)的特征量。表2所示為5例短路故障電流的特征量；表3所示為5例振蕩電流的特征量。

表2 短路故障電流特征量

表3 振蕩電流特征量

由表2和表3可以看出，對(duì)于多例電流信號(hào)，分別用兩種方法進(jìn)行計(jì)算，短路故障電流信號(hào)的特征量都比1大，而振蕩電流信號(hào)的特征量均比1小。分析表明，信號(hào)的樣本熵值與其規(guī)則程度緊密相關(guān)。而將EEMD分解后各IMF分量樣本熵之和作為特征量，相當(dāng)于在多個(gè)時(shí)間尺度對(duì)電流信號(hào)進(jìn)行分析，通過(guò)多個(gè)尺度，可以得到饋線電流信號(hào)中更深層次的信息，更好地對(duì)電流信號(hào)進(jìn)行分類。因此，用EEMD樣本熵值作為特征量對(duì)牽引網(wǎng)短路故障電流和振蕩電流進(jìn)行識(shí)別是可行的。

利用折線圖來(lái)更直觀地表示兩種電流信號(hào)特征量之間的關(guān)系，如圖6所示。

圖6 電流信號(hào)的特征量

由圖可知，雖然通過(guò)方法1和本文方法（方法2）得到的短路電流的特征量均大于1，振蕩電流的特征量都小于1，但本文方法得到的結(jié)果偏離1的程度更大，即得到的兩種電流信號(hào)的特征量差異更大，更有利于信號(hào)的識(shí)別。此外，利用方法1對(duì)多組電流信號(hào)進(jìn)行分析，得到的結(jié)果波動(dòng)較大，而本文方法得到的結(jié)果波動(dòng)較小。顯然，本文方法能更有效地識(shí)別直流牽引網(wǎng)短路故障電流和振蕩電流。

采用本文方法，在保護(hù)啟動(dòng)前200 ms內(nèi)，直流牽引網(wǎng)短路故障電流信號(hào)的EEMD樣本熵值E∈[1.871，1.932]，直流牽引網(wǎng)振蕩電流信號(hào)的EEMD樣本熵值E∈[0.151，0.173]。因此，在200ms的數(shù)據(jù)窗長(zhǎng)內(nèi)，定義EEMD樣本熵值作為特征量且整定閾值設(shè)為1，將計(jì)算得到的特征量與1進(jìn)行比較可判斷電流信號(hào)的類型，比1大則為短路故障電流信號(hào)。為增加保護(hù)算法的可靠性，避免保護(hù)誤動(dòng)，可以將整定閾值設(shè)為1.4，特征量＞1.4，則為短路故障電流信號(hào)。

5 結(jié)束語(yǔ)

地鐵直流牽引網(wǎng)中出現(xiàn)的振蕩電流容易造成牽引網(wǎng)饋線主保護(hù)頻繁誤動(dòng)。通過(guò)分析牽引網(wǎng)短路故障電流和振蕩電流的波形形態(tài)差異，本文提出了一種EEMD與樣本熵相結(jié)合的直流牽引網(wǎng)短路故障電流與振蕩電流識(shí)別方法，為直流牽引網(wǎng)故障診斷提供了一種新的思路。通過(guò)對(duì)直流牽引網(wǎng)饋線電流信號(hào)的分析計(jì)算可以看出，在直流牽引網(wǎng)振蕩電流與短路故障電流識(shí)別方面，本文方法更具有效性和準(zhǔn)確性。將本文方法作為饋線主保護(hù)的后備保護(hù)，可有效防止由振蕩電流造成的直流牽引網(wǎng)di/dt-ΔI保護(hù)誤動(dòng)。

[1]HUANG S R，CHEN B N．Harmonic study of the Le Blanc transformerforTaiwan railway’selectrification system[J]．IEEE Transactions on Power Delivery，2002，17（2）：495-499．

[2]祝冰心．直流牽引網(wǎng)故障電流識(shí)別算法研究[D]．北京：北京化工大學(xué)，2014．

[3]李夏青，左麗．直流牽引網(wǎng)故障電流能量譜及特征矢量[J]．電工技術(shù)學(xué)報(bào)，2010，25（11）：164-168．

[4]倪濤，李夏青．基于EMD分解的直流牽引供電系統(tǒng)短路故障識(shí)別[J]．北京石油化工學(xué)院學(xué)報(bào)，2012，20（2）：13-17．

[5]田行軍，李夏青，李運(yùn)華．直流牽引網(wǎng)振蕩電流與故障電流波形識(shí)別[J]．電工技術(shù)學(xué)報(bào)，2013，28（11）：247-253．

[6]祝冰心，李夏青，田行軍，等．基于關(guān)聯(lián)維數(shù)的直流牽引網(wǎng)故障識(shí)別[J]．北京石油化工學(xué)院學(xué)報(bào)，2014，22（1）：37-43．

[7]LI Y H，TIAN X J，LI X Q．Nonlinear second-order model and equilibrium point characteristic analysis of DC traction powersupplysystem[J]．ISRN Applied Mathematics，2013（11）：1-8．

[8]HUANG N E，SHEN Z，LONG S R．The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]．Proceedingsofthe RoyalSociety ofLondon：SeriesA，1998，454（1971）：903-995．

[9]安相璧，郭正，陳成法，等．基于EMD與振動(dòng)綜合度的再制造發(fā)動(dòng)機(jī)磨合質(zhì)量評(píng)價(jià)研究[J]．中國(guó)測(cè)試，2015，41（12）：129-133．

[10]張淑清，翟欣沛，董璇，等．EMD及Duffing振子在小電流系統(tǒng)故障選線方法中的應(yīng)用[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33（10）：161-167．

[11]文成，周傳德.基于改進(jìn)EMD和形態(tài)濾波的滾動(dòng)軸承故障診斷[J]．中國(guó)測(cè)試，2016，42（1）：135-139．

[12]崔心瀚，馬立元，魏忠林，等．EMD改進(jìn)方法研究及其在燃?xì)廨啓C(jī)工頻特征提取中的應(yīng)用[J]．中國(guó)測(cè)試，2016，42（1）：121-127．

[13]孫文軍，芮國(guó)勝，張馳，等.機(jī)械故障的級(jí)聯(lián)共振和EMD聯(lián)合參數(shù)辨識(shí)方法[J].電子測(cè)量與儀器學(xué)報(bào)，2016，30（3）: 352-360.

[14]WU Z，HUANG N E．Ensemble empirical mode decomposition：anoise-assisted dataanalysismethod[J]．Advances in adaptive data analysis，2009，1（1）：1-41．

[15]董青迅，李興源，穆子龍，等．基于近似熵的次同步諧振檢測(cè)方法[J]．電網(wǎng)技術(shù)，2011，35（11）：103-106．

[16]張學(xué)清，梁軍，張熙，等．基于樣本熵和極端學(xué)習(xí)機(jī)的超短期風(fēng)電功率組合預(yù)測(cè)模型[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33（25）：33-40．

[17]鄭近德，程軍圣，楊宇.基于多尺度熵的滾動(dòng)軸承故障診斷方法[J].湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2015（5）：38-41.

[18]龔?fù)⒅荆绷鳡恳╇娤到y(tǒng)數(shù)學(xué)模型與短路計(jì)算研究[D]．北京：北京交通大學(xué)，2009．

[19]喻樂(lè).城市軌道交通供電系統(tǒng)建模與直流饋線保護(hù)的研究[D]．北京：北京交通大學(xué)，2012．

[20]王奎鵑．直流牽引供電系統(tǒng)電磁暫態(tài)仿真[D]．北京：北京化工大學(xué)，2013．

[21]王建民．直流牽引供電系統(tǒng)振蕩機(jī)理分析及抑制策略[D]．北京：北京化工大學(xué)，2015．

（編輯：李妮）

Fault current identification of DC traction network based on EEMD sample entropy

LENG Yue，YANG Honggeng，WANG Zhiqi
（School of Electrical Engineering and Information，Sichuan University，Chengdu 610065，China）

Aiming at the frequent malfunctions of relay protection system caused by the oscillation current of Metro DC traction network，an identification approach of the oscillation current and the short-circuit fault current in DC traction network combined ensemble empirical mode decomposition（EEMD）with sample entropy is presented.In this method，the feeder current is decomposed by using EEMD method，and the sum of the sample entropy of each intrinsic mode function（IMF）is figured out.Thus，the feature which involves information on the operation state of DC traction network is acquired.From the analytical calculation of the typical feeder currents，it can be found outthatthefeatureextraction method combined EEMD and sampleentropycan effectively distinguish the oscillation current and the short-circuit fault current of DC traction network. Examples verify the effectiveness of the method.

short-circuit fault current；oscillation current；sample entropy；ensemble empirical mode decomposition；DC traction network

A

：1674-5124（2016）12-0095-05

10.11857/j.issn.1674-5124.2016.12.020

2016-05-08；

：2016-06-12

冷 月（1992-），女，重慶市人，碩士研究生，專業(yè)方向?yàn)檐壍澜煌恳╇娤到y(tǒng)保護(hù)與控制。