基于EEMD分解與相關(guān)系數(shù)的變壓器振動(dòng)信號(hào)特征量提取

徐舒蓉，劉君，曾華榮，呂黔蘇，李堃，楊旗，趙莉華

(1.貴州電網(wǎng)有限責(zé)任公司電力科學(xué)研究院，貴州 貴陽(yáng)，550002；2.南方電網(wǎng)有限責(zé)任公司防冰減災(zāi)重點(diǎn)實(shí)驗(yàn)室，貴州 貴陽(yáng)，550002；3.四川大學(xué)電氣工程學(xué)院，四川 成都，610065)

變壓器振動(dòng)信號(hào)的提取可實(shí)現(xiàn)與被采集設(shè)備無(wú)電氣連接，并且具有較強(qiáng)抗干擾能力、靈敏度也比較高，因此基于振動(dòng)信號(hào)的變壓器機(jī)械故障診斷法是近年來(lái)變壓器在線檢測(cè)技術(shù)的研究熱點(diǎn)[1-3]。

基于變壓器振動(dòng)信號(hào)的故障診斷技術(shù)研究主要包括有限元仿真建模與試驗(yàn)研究、信號(hào)處理兩大方向。Berler等學(xué)者通過(guò)試驗(yàn)采集變壓器空載、負(fù)載工況下的振動(dòng)信號(hào)，通過(guò)分析其時(shí)域特征來(lái)評(píng)價(jià)變壓器的線圈壓緊情況[4]；Belen Garcia，Juan Carlos Burgos 等學(xué)者定量研究了電網(wǎng)頻率、電壓與電流大小、溫度等因素對(duì)變壓器振動(dòng)信號(hào)波形的影響規(guī)律，并分析了振動(dòng)的傳播路徑變化情況[5-7]；汲勝昌團(tuán)隊(duì)也研究了變壓器振動(dòng)與空載電壓、負(fù)載電流的關(guān)系，并且通過(guò)仿真搭建變壓器箱體振動(dòng)信號(hào)傳播模型，運(yùn)用該模型初步成功診斷出繞組變形故障，驗(yàn)證了模型的有效性[8-14]。信號(hào)處理方面，傅里葉變換(FFT)是分析信號(hào)頻域特征最基本的算法，反映了信號(hào)全局頻率的分布[15-18]，有研究者提出基于變壓器振動(dòng)信號(hào)基頻占比以及高次諧波占比等作為變壓器狀態(tài)診斷的特征量[19-20]。隨著信號(hào)處理技術(shù)的進(jìn)一步發(fā)展，學(xué)者提出基于經(jīng)驗(yàn)?zāi)B(tài)分解的希爾伯特黃變換和基于小波變換的時(shí)頻分析法[21-22]。

需要指出的是，目前各研究者得出的結(jié)論因試驗(yàn)方法、測(cè)試條件、樣品差異等因素導(dǎo)致其研究結(jié)果沒(méi)有統(tǒng)一性，并且由于缺乏大量數(shù)據(jù)支撐而無(wú)法下定論，并且已有的研究基本只能實(shí)現(xiàn)正常、異常兩種判斷，對(duì)于利用振動(dòng)信號(hào)判斷變壓器故障類(lèi)型還亟須研究。

本文在相關(guān)學(xué)者的研究基礎(chǔ)上，提出一種基于相關(guān)系數(shù)法與改進(jìn)集合經(jīng)驗(yàn)?zāi)Ｊ椒纸?Ense-mble Empirical Mode Decomposition，EEMD)的頻率、峭度、幅值、均值等特征量提取方法。旨在為建立相對(duì)統(tǒng)一的變壓器振動(dòng)信號(hào)故障診斷方案提供理論參考。首先，根據(jù)實(shí)際經(jīng)驗(yàn)設(shè)計(jì)濾波器對(duì)信號(hào)進(jìn)行濾波，對(duì)濾波之后的振動(dòng)信號(hào)進(jìn)行EEMD分解得到一系列IMF(Intrinsic Mode Function，IMF)分量，運(yùn)用相關(guān)系數(shù)法求取敏感IMF分量并對(duì)敏感IMF分量進(jìn)行重構(gòu)的信號(hào)處理方法，求重構(gòu)信號(hào)的特征頻率、峭度、幅值與均值特征量以構(gòu)成特征向量。最后，運(yùn)用該方法進(jìn)行了試驗(yàn)變壓器及不同故障狀態(tài)下在運(yùn)變壓器的振動(dòng)信號(hào)特征向量提取，研究其變化規(guī)律。

1 變壓器振動(dòng)信號(hào)試驗(yàn)數(shù)據(jù)采集

運(yùn)用PCB 356A16壓電式加速度傳感器與NI9234數(shù)據(jù)收集儀采集型號(hào)為S11-M-10/10的試驗(yàn)變壓器振動(dòng)信號(hào)，設(shè)置采樣頻率25.6kHz。振動(dòng)信號(hào)采集試驗(yàn)接線方法、信號(hào)采集點(diǎn)的布置分別如圖1、圖2所示。

圖1 變壓器振動(dòng)信號(hào)采集接線圖Fig. 1 Vibration signal acquisition system

如圖1所示，通過(guò)調(diào)壓器給試驗(yàn)變壓器一次側(cè)供電。圖2中，將傳感器分別布置在箱體頂部和正面三相繞組所對(duì)應(yīng)位置，測(cè)試該試驗(yàn)變壓器空載、負(fù)載工況下的振動(dòng)數(shù)據(jù)。

對(duì)于在運(yùn)電力變壓器的振動(dòng)信號(hào)采集，其油箱頂部位置是高壓危險(xiǎn)區(qū)，通常情況下將傳感器安裝于變壓器正面。對(duì)于三相變壓器，傳感器分別布置在A、B、C三相中部位置；對(duì)于單相變壓器，將傳感器分別布置在變壓器正面箱體的上、中、下三個(gè)位置。

圖2 變壓器振動(dòng)信號(hào)采集點(diǎn)的布置圖Fig.2 Test point distribution map

2 EEMD分解算法、相關(guān)系數(shù)法

2.1 EEMD分解算法

2.1.1 EEMD分解過(guò)程

集合經(jīng)驗(yàn)?zāi)Ｊ椒纸?EEMD)是一種將白噪聲作為輔助來(lái)處理數(shù)據(jù)的方法[23-24]。其主要過(guò)程：在被分析信號(hào)中添加同樣的高斯白噪聲，利用高斯白噪聲頻域譜均勻分布的特點(diǎn)使待分析信號(hào)自動(dòng)分布到適當(dāng)?shù)膮⒖汲叨戎小Ｈ缓?，基于白噪聲信?hào)的零均值特點(diǎn)，經(jīng)過(guò)多次平均使加入的白噪聲消除，則整體平均值為最終結(jié)果。其中，待分析信號(hào)不會(huì)因此而發(fā)生改變[25]。

EEMD是EMD的改進(jìn)算法，EEMD算法沒(méi)有EMD經(jīng)驗(yàn)?zāi)B(tài)混疊缺陷，但是卻具有EMD自適應(yīng)性、正交性、完備性等優(yōu)點(diǎn)。EEMD算法分為如下幾個(gè)步驟[26-27]：

(1)參數(shù)歸一化：由待分析信號(hào)的標(biāo)準(zhǔn)偏差確定集成次NE、所需添加白噪聲的幅值a。

(2)同樣根據(jù)EMD算法對(duì)添加了白噪聲后的信號(hào)xm進(jìn)行分解：

xm=x+nm

(1)

式(1)中：nm為添加的白噪聲，x為待分析信號(hào)。分解后的結(jié)果稱(chēng)為本征模函數(shù)(IMFs)，記為ci,m，分解后余下的非零信號(hào)rm稱(chēng)為殘余函數(shù)。

(3)然后，設(shè)置m=m+1，循環(huán)上述步驟(2)，直到m等于初始時(shí)設(shè)置的集成次數(shù)NE。

(4)求取上述分解得到的IMFs總體平均值，該值即為EEMD分解結(jié)果。

2.1.2 EEMD參數(shù)設(shè)定

EEMD中最關(guān)鍵的兩個(gè)參數(shù)是所添加白噪聲的幅值a與集成的次數(shù)NE，EEMD分解結(jié)果的性能好與否與它們直接關(guān)聯(lián)，因此需對(duì)其進(jìn)行合理選擇。集成的次數(shù)NE、白噪聲的幅值a與標(biāo)準(zhǔn)差σ之間有關(guān)系[28]：

(2)

EEMD分解時(shí)，所添加的白噪聲幅值越小，則分解得到的IMF分量與原始信號(hào)的誤差越?。坏?，幅值小，信噪比高，添加的噪聲對(duì)極點(diǎn)的選取影響小，失去補(bǔ)充尺度的作用。應(yīng)合理選擇其幅值，一般a取0.2。在確定所添加白噪聲幅值的情況下，集成的次數(shù)越多，則因噪聲產(chǎn)生的誤差越小。但是次數(shù)多計(jì)算量大，因此一般集成次數(shù)設(shè)置為100。

2.2 敏感IMF分量提取

信號(hào)經(jīng)過(guò)EEMD分解后，會(huì)產(chǎn)生一系列IMF分量，各IMF分量與原信號(hào)的相關(guān)性各不相同，習(xí)慣把相關(guān)性強(qiáng)的IMF視為“敏感IMF分量”。在工程運(yùn)用中，需挑選與原信號(hào)相關(guān)性強(qiáng)，即承載主要信息的敏感IMF分量進(jìn)行分析。本文根據(jù)IMF分量與原信號(hào)相關(guān)系數(shù)的大小得到敏感IMF，相關(guān)系數(shù)>0.5則認(rèn)為是敏感分量[29-30]。

2.3 振動(dòng)信號(hào)特征量的選取

本文選擇特征頻率、峭度、幅值和均值作為振動(dòng)信號(hào)的特征量。

(1)特征頻率(fc)

理論研究表明，變壓器鐵心、繞組線圈振動(dòng)信號(hào)基波頻率為100Hz，由于非線性影響因素導(dǎo)致200Hz、300Hz等倍頻分量的產(chǎn)生。當(dāng)變壓器鐵心和繞組故障時(shí)，主頻率發(fā)生變化，因此本文將振動(dòng)信號(hào)頻譜中幅值較大的頻率定義為特征頻率。

(2)峭度(k)

峭度實(shí)質(zhì)為一個(gè)數(shù)學(xué)統(tǒng)計(jì)量，可以體現(xiàn)波形的尖峰度，為無(wú)量綱參數(shù)。峭度一般用字母k表示，其定義式如公式(3)所示[16]：

(3)

其中，μ為信號(hào)x的平均值，σ為信號(hào)x的標(biāo)準(zhǔn)差，E(t)為變量t的期望值。當(dāng)變壓器狀態(tài)發(fā)生改變時(shí)，其振動(dòng)信號(hào)時(shí)域波形的尖峰度k也會(huì)出現(xiàn)相應(yīng)的變化。

(3)幅值(m)

變壓器振動(dòng)信號(hào)是周期信號(hào)，其幅值m也是一個(gè)重要參數(shù)(單位是g，即重力加速度)，變壓器運(yùn)行狀態(tài)異常情況下，幅值也會(huì)發(fā)生變化。

(4)均值(q)

對(duì)于變壓器振動(dòng)信號(hào)，均值可體現(xiàn)波形的對(duì)稱(chēng)程度(單位為g)，正弦信號(hào)均值應(yīng)該為0，信號(hào)均值越趨近0則說(shuō)明該信號(hào)畸變較小，即正、負(fù)半周對(duì)稱(chēng)性較好。

2.4 特征量提取流程

基于快速譜峭度圖法、EEMD分解、相關(guān)系數(shù)法的變壓器振動(dòng)信號(hào)特征量提取流程如圖3示，其流程為：

(1)根據(jù)實(shí)際經(jīng)驗(yàn)理論設(shè)計(jì)濾波器對(duì)信號(hào)進(jìn)行濾波。

(2)對(duì)(1)中濾波后的信號(hào)進(jìn)行EEMD分解，得到一系列IMF分量。

(3)運(yùn)用相關(guān)系數(shù)計(jì)算法獲取各IMF分量與原信號(hào)的相關(guān)系數(shù)大小，取系數(shù)>0.5的分量為敏感IMF分量。

(4)重構(gòu)(3)中的敏感IMF分量，求重構(gòu)信號(hào)的特征頻率、峭度、幅值、均值特征量，構(gòu)成變壓器振動(dòng)信號(hào)特征向量。

圖3 特征量提取流程圖Fig.3 The flow chart of feature extraction

3 試驗(yàn)變壓器的振動(dòng)信號(hào)特征向量提取

3.1 正常振動(dòng)信號(hào)特征量提取

提取試驗(yàn)變壓器空載振動(dòng)信號(hào)特征量，空載時(shí)二次繞組無(wú)電流，一次繞組電流很小，此時(shí)主要為鐵心振動(dòng)。由圖2可知，位于正面中部的第 B5測(cè)點(diǎn)以最短路徑獲取高、低壓側(cè)的振動(dòng)信息，為最佳監(jiān)測(cè)點(diǎn)。這里對(duì)B5測(cè)點(diǎn)的信號(hào)進(jìn)行分析，圖4展示了振動(dòng)信號(hào)的原始波形、頻譜。

由圖4可知，鐵心振動(dòng)信號(hào)波形畸變較嚴(yán)重，存在<50Hz的低頻干擾分量及少量基頻的倍頻分量，這是由于環(huán)境中干擾信號(hào)疊加以及鐵心材料非線性導(dǎo)致。

(a)原始波形(a) Original waveform

(b)原始波形頻譜(b) Spectrum of oringinal waveform圖4 B5測(cè)點(diǎn)振動(dòng)信號(hào)原始波形與頻譜圖Fig.4 The original signal waveform and spectrum of vibration a B5

以B5測(cè)點(diǎn)振動(dòng)信號(hào)的4000個(gè)數(shù)據(jù)為一個(gè)數(shù)據(jù)段，根據(jù)實(shí)際經(jīng)驗(yàn)理論設(shè)計(jì)濾波器對(duì)原信號(hào)進(jìn)行濾波，圖5為濾波后的信號(hào)波形，對(duì)該波形進(jìn)行EEMD分解，得到圖6所示的一系列IMF分量。表1為IMF分量與原信號(hào)的相關(guān)系數(shù)值。

圖5 濾波后的振動(dòng)信號(hào)波形Fig.5 Signal waveform after filtering

(a)IMF1～I(xiàn)MF6 (b)IMF7～I(xiàn)MF12圖6 重構(gòu)信號(hào)的EEMD分解后得到的IMFFig.6 The decomposition IMFs using EEMD

據(jù)表1可選取相關(guān)系數(shù)>0.5的IMF6、IMF7、IMF8這三個(gè)分量為敏感分量。根據(jù)三個(gè)分量對(duì)信號(hào)進(jìn)行重構(gòu)并求包絡(luò)譜，圖7為重構(gòu)后的波形。

(a)重構(gòu)后的信號(hào)波形(a)Reconstructed signal waveform

(b)重構(gòu)信號(hào)頻譜(b)Reconstructed signal spectrum圖7 重構(gòu)信號(hào)波形及重構(gòu)信號(hào)頻譜圖Fig.7 Reconstructed signal waveform and reconstructed signal spectrum

據(jù)圖7可知其時(shí)域波形比圖5更平滑、清晰；該變壓器空載正常運(yùn)行時(shí)的特征頻率fc=100Hz，消除了原始信號(hào)頻譜圖中存在的低頻干擾，其能量聚集性更高。最后求該重構(gòu)信號(hào)的峭度k、幅值m、均值q，計(jì)算結(jié)果分別為：k=2.64，m=0.0069g，q=4.9×10-5g。綜合四個(gè)參數(shù)構(gòu)成試驗(yàn)變壓器正?？蛰d運(yùn)行時(shí)振動(dòng)信號(hào)的特征向量：

Tkz=[100 2.64 0.0069 4.9×10-5]

3.2 變壓器鐵心松動(dòng)后空載振動(dòng)信號(hào)

模擬變壓器鐵心松動(dòng)故障，進(jìn)行空載狀態(tài)振動(dòng)信號(hào)采集，同取B5測(cè)點(diǎn)進(jìn)行分析。其原始時(shí)域波形及頻譜如圖8所示。由圖8可知其時(shí)域波形與正常振動(dòng)信號(hào)波形相似；但是其主頻增大為200Hz。不僅存在低頻干擾信號(hào)，并且包含高頻干擾分量。

圖9是濾波后的信號(hào)波形，同理，為了進(jìn)一步消除低頻干擾，對(duì)該信號(hào)進(jìn)行EEMD分解，提取敏感分量進(jìn)行信號(hào)重構(gòu)。求得IMFs與原信號(hào)的相關(guān)系數(shù)結(jié)果如表1。

(a)原始波形(a)Original waveform

(b)原始波形頻譜(b)Spectrum of oringinal waveform圖8 故障狀態(tài)原始信號(hào)波形與頻譜Fig.8 The original signal waveform and spectrum in fault status

圖9 空載故障振動(dòng)信號(hào)濾波后信號(hào)波形Fig.9 The signal waveform after filtering

表1 IMFs與原信號(hào)的相關(guān)系數(shù)Tab.1 The correlation coefficient between IMFs and original signal

據(jù)表1可選取相關(guān)系數(shù)>0.5的IMF7、IMF8、IMF9為敏感分量對(duì)信號(hào)進(jìn)行重構(gòu)，并求其包絡(luò)譜，結(jié)果如圖10，由圖10可知重構(gòu)后的信號(hào)時(shí)域波形更加平滑，鐵心松動(dòng)變壓器空載時(shí)的特征頻率f=200Hz，存在400Hz等倍頻分量。與圖8中頻譜圖對(duì)比，消除了高頻干擾成分，能量聚集性更好、倍頻分辨率更高。計(jì)算得峭度、幅值、均值分別為：k=3.60，m=0.0083g，q=0.0004g。綜上得試驗(yàn)變壓器空載鐵心松動(dòng)故障振動(dòng)信號(hào)的特征向量：

Tkzgz=[200 3.60 0.0083 0.0004]

與特征向量Tkz相比較，Tkzgz特征頻率發(fā)生變化，峭度、幅值與均值呈增大趨勢(shì)。

(a)重構(gòu)后的信號(hào)波形(a)Reconstructed signal waveform

(b)重構(gòu)信號(hào)頻譜(b)Reconstructed signal spectrum圖10 重構(gòu)信號(hào)波形、重構(gòu)信號(hào)頻譜圖Fig.10 The reconstruct signal waveform and its spectrum

3.3 試驗(yàn)變壓器空載正常、鐵心松動(dòng)故障下特征量

綜上分析，可將試驗(yàn)變壓器正常、鐵心松動(dòng)故障狀態(tài)下的空載振動(dòng)信號(hào)的特征頻率、峭度、幅值、均值構(gòu)成的特征向量總結(jié)如表2所示。正常狀態(tài)下，振動(dòng)信號(hào)的特征頻率均為100Hz，與理論相符；鐵心松動(dòng)故障時(shí)，特征頻率增大為200Hz，100Hz幅值減小幾乎減小為零，并且出現(xiàn)400Hz分量；其峭度值增大、振動(dòng)幅值、均值都增大?？梢?jiàn)通過(guò)用本文提出的方法求振動(dòng)信號(hào)特征量，可判斷該變壓器出現(xiàn)異常。

表2 試驗(yàn)變壓器正常與故障狀態(tài)下特征向量Tab.2 The eigenvectors of test transformer under normal and fault state

4 實(shí)際工程應(yīng)用

4.1 鐵心磁路故障

將本文提出的提取特征量方法運(yùn)用到現(xiàn)場(chǎng)運(yùn)行變壓器中。對(duì)某變電站型號(hào)為ODFPS-250000/500 的500kV單相變壓器正常狀態(tài)與鐵心磁路故障狀態(tài)進(jìn)行振動(dòng)信號(hào)測(cè)試與分析。

4.1.1 正常運(yùn)行狀態(tài)下的特征向量提取

該變壓器正常運(yùn)行時(shí)原始時(shí)域波形與頻譜如圖11所示，其中100Hz的整數(shù)倍頻信號(hào)較豐富，這與變壓器負(fù)載率、溫度、電流畸變率等影響因素有關(guān)。

濾波后的波形如圖12所示。通過(guò)相關(guān)系數(shù)法選取EEMD分解后的敏感IMF分量，相關(guān)系數(shù)計(jì)算結(jié)果如表3。

(a)原始波形(a)Original waveform

(b)原始信號(hào)頻譜(b)Spectrum of oringinal signal圖11 實(shí)際運(yùn)行變壓器正常振動(dòng)信號(hào)波形及頻譜Fig.11 The normal vibration signal waveform and frequency spectrum of an actual operation transformer

圖12 濾波后時(shí)域信號(hào)波形Fig12 Signal waveform after filtering

表3 IMFs與原信號(hào)的相關(guān)系數(shù)Tab.3 The correlation coefficient between IMFs and original signal

由表3可選取相關(guān)系數(shù)>0.5的IMF5、IMF6、IMF7、IMF8這四個(gè)分量進(jìn)行重構(gòu)后并求其包絡(luò)譜，結(jié)果如圖13所示。與圖11中原始頻譜相比，經(jīng)過(guò)濾波和重構(gòu)后的頻譜能量聚集性更好，倍頻分辨率更高，消除了大部分低頻與高頻干擾成分。該變壓器正常運(yùn)行時(shí)特征頻率范圍是100Hz～500 Hz，主頻為200Hz，400Hz與500Hz幅值也較大。在MATLAB中求得其峭度、幅值、均值分別為：k=2.73，m=0.169g，q=0.0001g。綜上可得該變壓器正常情況下的特征向量Tzc：

Tzc=[200 2.73 0.169 0.0001]

(a)重構(gòu)后的信號(hào)波形(a)Reconstructed signal waveform

(b)重構(gòu)信號(hào)頻譜(b)Reconstructed signal spectrum圖13 重構(gòu)振動(dòng)信號(hào)波形及重構(gòu)信號(hào)頻譜Fig.13 Reconstructed signal waveform and reconstructed signal spectrum

4.1.2 磁路故障運(yùn)行狀態(tài)下的特征向量提取

該變壓器發(fā)生磁路故障下的振動(dòng)原始波形與頻譜如圖14所示，發(fā)生磁路故障后的振動(dòng)信號(hào)振動(dòng)幅度變大，每個(gè)周期內(nèi)的波形出現(xiàn)明顯的尖峰；振動(dòng)頻譜成分變復(fù)雜，出現(xiàn)大量100Hz的倍頻分量及50Hz的奇次倍頻分量。

(a)原始波形(a)Original waveform

(b)原始信號(hào)頻譜(b)Spectrum of oringinal signal圖14 實(shí)際運(yùn)行變壓器磁路故障振動(dòng)信號(hào)及頻譜Fig.14 The waveform and spectrum of actual operation transformer vibration signal under magnetic circuit fault

圖15 濾波后的振動(dòng)信號(hào)波形Fig.15 Signal waveform after filtering

表4 IMFs與原信號(hào)的相關(guān)系數(shù)Tab.4 The correlation coefficient between IMFs and original signal

對(duì)原信號(hào)濾波之后的波形如圖15，對(duì)圖15中信號(hào)波形進(jìn)行EEMD分解并求相關(guān)系數(shù)，結(jié)果如表4所示。

據(jù)表4可取IMF5、IMF6、IMF7為敏感分量進(jìn)行重構(gòu)并求其包絡(luò)譜，結(jié)果如圖16，可知重構(gòu)信號(hào)去掉了部分相關(guān)性小的分量，該變壓器磁路故障時(shí)，主頻仍然為200Hz，但出現(xiàn)大量50Hz的奇次倍頻分量，可判斷這是該變壓器磁路故障下出現(xiàn)的最具特征的頻率成分，這里以(200Hz，250Hz)作為該故障狀態(tài)下的頻率特征。同理求得該重構(gòu)信號(hào)的峭度、幅值、均值為：k=3.72，m=0.266g，q=0.0003g。綜上所述得該變壓器磁路故障下的特征向量Tclgz為：

Tclgz=[(200, 250) 3.72 0.266 0.0003]

(a)重構(gòu)后的信號(hào)波形(a)Reconstructed signal waveform

(b)重構(gòu)信號(hào)頻譜(b)Reconstructed signal spectrum圖16 濾波信號(hào)、重構(gòu)信號(hào)波形及重構(gòu)信號(hào)頻譜圖Fig.16 Reconstructed signal waveform and reconstructed signal spectrum

4.1.3 實(shí)測(cè)變壓器正常與磁路故障下特征向量對(duì)比

綜上所述，得該變壓器正常與磁路故障下的特征向量如表5所示。與正常狀態(tài)下Tzc相比，Tclgz出現(xiàn)了250Hz頻率特征，峭度增大1.01，振動(dòng)幅值增大0.56倍，均值增大2倍?？蓪⑻卣飨蛄縏clgz可作為該變壓器磁路故障的判據(jù)。

表5 實(shí)際運(yùn)行變壓器正常與故障狀態(tài)下特征向量Tab.5 The feature vectors of actual operation transformer in normal and fault state

4.2 鐵心多點(diǎn)接地故障

對(duì)某220kV變電站型號(hào)為SFPSZ10-1800 00/220的三相變壓器正常狀態(tài)與鐵心多點(diǎn)接地故障狀態(tài)進(jìn)行振動(dòng)信號(hào)測(cè)試與分析。

4.2.1 正常狀態(tài)

運(yùn)用本文處理方法對(duì)該變壓器正常運(yùn)行時(shí)振動(dòng)信號(hào)進(jìn)行濾波、提取相關(guān)系數(shù)并重構(gòu)后的時(shí)域波形與頻譜如圖17。

(a)重構(gòu)信號(hào)時(shí)域波形(a)Time domain waveform of reconstructed signal

(b)重構(gòu)信號(hào)頻域波形(b)Frequency domain waveform of reconstructed signal圖17 重構(gòu)信號(hào)波形及重構(gòu)信號(hào)頻譜圖Fig.17 Reconstructed signal waveform and reconstructed signal spectrum

據(jù)圖17可知，經(jīng)處理后該變壓器正常狀態(tài)下該變壓器時(shí)域波形清晰干凈。主頻為300Hz，低頻幅值較大，大于400Hz的分量較少。計(jì)算其峭度、幅值、均值得其特征向量為：

Tzc=[300Hz 2.15 0.10 0.00075]

4.2.2 鐵心多點(diǎn)接地故障

同理得到該變壓器發(fā)生鐵心多點(diǎn)接地故障時(shí)重構(gòu)后的振動(dòng)信號(hào)時(shí)域及頻譜圖18，可知故障后的信號(hào)波形存在較多毛刺、頻率范圍整體擴(kuò)大，基頻的倍頻分量幅值明顯增大并且主頻為增大為500Hz，700Hz處的幅值也較大。計(jì)算其峭度、幅值、均值得鐵心多點(diǎn)接地故障下的特征向量為：

Tgz=[500Hz 2.83 0.20 0.00115]

該故障情況下不僅主頻增大，峭度、幅值、均值都增大。但峭度值變化較小。

(a)重構(gòu)信號(hào)時(shí)域波形(a)Time domain waveform of reconstructed signal

(b)重構(gòu)信號(hào)頻域波形(b)Frequency domain waveform of reconstructed signal圖18 鐵心多點(diǎn)接地故障重構(gòu)信號(hào)波形及重構(gòu)信號(hào)頻譜Fig.18 Reconstructed signal waveform and reconstructed signal spectrum under Iron core multi-point grounding fault

5 結(jié)論

綜上所述，本文運(yùn)用一種基于相關(guān)系數(shù)法與EEMD分解的頻率、峭度、幅值、均值等特征量提取方法，對(duì)試驗(yàn)變壓器正常、鐵心振動(dòng)信號(hào)進(jìn)行特征向量提取及分析，并將其運(yùn)用到實(shí)際運(yùn)行的變壓器中，研究磁路故障、鐵心多點(diǎn)接地故障下的特征量。得結(jié)論如下：

(1)EEMD克服了EMD模態(tài)混疊的缺陷，有效地將振動(dòng)信號(hào)分解為不同頻率的IMF分量；相關(guān)系數(shù)法可正確提取敏感IMF。經(jīng)過(guò)濾波與重構(gòu)后的變壓器振動(dòng)信號(hào)保留了重要分量，濾除了大部分的低頻、高頻干擾成分，使信號(hào)頻譜能量聚集性更好，倍頻分辨率更高。

(2)運(yùn)用本文提出的特征提取流程對(duì)試驗(yàn)變壓器空載正常、空載鐵心松動(dòng)狀態(tài)下的振動(dòng)信號(hào)進(jìn)行分析，結(jié)果表明在該變壓器空載時(shí)特征頻率為100Hz，峭度值集中在2.0～2.6之間；空載鐵心松動(dòng)狀態(tài)下特征頻率為200Hz，峭度值、幅值與均值都增大；本文提出的特征量能有效區(qū)分變壓器的故障狀態(tài)。

(3)經(jīng)實(shí)例研究分析可知在不同故障下，本文提取的特征向量中各參數(shù)均呈現(xiàn)增大趨勢(shì)，但是其增大的幅度因故障類(lèi)型而異。磁路故障時(shí)主頻大小的變化不明顯，但是峭度值變化較大；鐵心多點(diǎn)接地故障下峭度值變化較小而主頻則增大為500Hz；均值和幅值的增大倍數(shù)在不同故障狀態(tài)下的差異不太明顯。

綜上，本文研究的變壓器振動(dòng)信號(hào)特征量提取方法不僅可判斷變壓器是否異常，并且有望區(qū)分不同故障類(lèi)型的信號(hào)，進(jìn)而實(shí)現(xiàn)變壓器故障類(lèi)型診斷，對(duì)實(shí)際工程中變壓器的故障診斷具有一定的參考價(jià)值。