基于小波能量譜分析的變壓器振動特性

李少飛，馬興明，毛新宇

(國網(wǎng)黑龍江省電力有限公司大慶供電公司，黑龍江 大慶 163000)

0 引 言

電力變壓器作為電力系統(tǒng)的關(guān)鍵設(shè)備，也是最為昂貴、復(fù)雜的設(shè)備之一。由于電力變壓器處于電網(wǎng)的中心位置，運(yùn)行環(huán)境復(fù)雜，且時常遭受各種不良運(yùn)行工況的沖擊，一旦發(fā)生故障，極有可能造成突發(fā)大面積停電，甚至引起爆炸、火災(zāi)等事故，帶來較大的直接和間接經(jīng)濟(jì)損失[1-3]。因此對變壓器的運(yùn)行狀況進(jìn)行長期監(jiān)測，及時評估變壓器狀態(tài)，能夠有效減少突發(fā)故障帶來的損失。

變壓器振動信號能夠有效反映變壓器的運(yùn)行狀況及內(nèi)部狀態(tài)[4-6]，因此對變壓器振動問題的研究及對振動信號的分析是目前的研究熱點。文獻(xiàn)[5]利用變壓器振動測試系統(tǒng)，對在變壓器油箱表面測得的鐵芯及繞組的振動信號進(jìn)行研究，在理論研究影響鐵芯及繞組振動各因素的基礎(chǔ)上，實驗驗證了在空載狀態(tài)下油箱表面測得的鐵芯振動信號的基頻幅值與施加電壓的平方以及繞組振動信號的基頻幅值與流過繞組的負(fù)載電流平方都呈較好的線性關(guān)系。文獻(xiàn)[6]基于變壓器鐵芯及繞組噪聲及振動機(jī)理，建立振動噪聲測量系統(tǒng)，測量變壓器實際振動及噪聲信息，根據(jù)指紋成像噪聲源定位和振動信號特征值測試結(jié)果，對變壓器內(nèi)部故障進(jìn)行定位。文獻(xiàn)[7]針對變壓器繞組多種故障并發(fā)的工況，在分析變壓器繞組振動機(jī)理的基礎(chǔ)上，提出一種基于集合經(jīng)驗?zāi)Ｊ椒纸獾恼駝有盘柼崛》椒?，對變壓器繞組振動信號進(jìn)行分解得到各階本征模函數(shù)并構(gòu)造特征矢量，將該特征矢量作為變壓器繞組狀態(tài)識別的判據(jù)。以上研究對變壓器振動問題進(jìn)行了相關(guān)的建模與理論分析，但是對于變壓器長期帶電運(yùn)行下的振動監(jiān)測未進(jìn)行相關(guān)的研究。

本文針對變壓器穩(wěn)態(tài)運(yùn)行振動問題，建立變壓器振動分析原理及模型，測量變壓器帶電運(yùn)行狀態(tài)下箱體不同位置的振動信號，基于小波分解原理對振動信號進(jìn)行分解重構(gòu)及能量譜提取，實現(xiàn)對變壓器振動的辨識及內(nèi)部狀態(tài)的診斷。

1 變壓器振動機(jī)理

變壓器振動主要由變壓器本體(鐵芯、繞組)的振動和冷卻裝置的振動通過機(jī)械連接及絕緣油傳遞到變壓器表面。在變壓器箱體表面可采集到的各種振動信號中，油泵、風(fēng)扇等冷卻裝置的振動頻率與鐵芯、繞組振動有明顯區(qū)別，易于分離[8-9]。因此，分析變壓器箱體表面振動信號可得，振動信號主要是由變壓器本體產(chǎn)生的。

1.1 鐵芯振動機(jī)理

鐵芯振動主要原因是硅鋼片的磁致伸縮[10-11]。根據(jù)電磁感應(yīng)原理及鐵磁材料磁致伸縮特性，可得鐵芯磁致伸縮引起的振動加速度a1為

式中：L為鐵磁材料原始尺寸；εs為硅鋼片的飽和磁致伸縮率；Us為電源電壓幅值；N為一次側(cè)繞組匝數(shù)；S為鐵芯橫截面積；Bs為鐵芯飽和磁感應(yīng)強(qiáng)度；ω激勵頻率。

振動加速度在材料及運(yùn)行條件不變的情況下與電壓平方呈正比關(guān)系，振動頻率為激勵頻率的2倍。變壓器正常運(yùn)行時，激勵為工頻激勵，因此鐵芯振動的基頻為100 Hz，但是受到鐵磁材料的非線性及鐵芯內(nèi)框和外框的磁通路徑長短不同的影響，使得鐵芯內(nèi)出現(xiàn)以100 Hz為基頻的高次諧波。

1.2 繞組振動機(jī)理

處于磁場中的載流導(dǎo)體要承受機(jī)械力的作用，當(dāng)變壓器繞組中通過電流時，由于電流與漏磁場的作用，在繞組內(nèi)產(chǎn)生電磁機(jī)械力，其大小由通過繞組的電流與漏磁場的磁通密度大小共同決定[12-13]。作用在繞組線圈上的電動力F為

式中：p為電動力系數(shù)；Im為負(fù)載電流幅值；φ0為負(fù)載電流初相位。

根據(jù)變壓器繞組的機(jī)械結(jié)構(gòu)，常將變壓器繞組等效為彈簧質(zhì)量系統(tǒng)，根據(jù)動力學(xué)定理可得繞組振動加速度由一個穩(wěn)態(tài)分量及逐漸衰減的分量組成，在理想狀態(tài)下穩(wěn)定運(yùn)行的變壓器繞組的振動加速度信號是電源頻率的2倍，對于工頻網(wǎng)絡(luò)即為100 Hz。振動加速度大小正比于繞組電流的平方。而由于加速度受絕緣材料彈性系數(shù)影響，絕緣材料彈性系數(shù)不是常數(shù)，具有較強(qiáng)的非線性特性，導(dǎo)致繞組振動在負(fù)載較大時呈明顯的非線性特征且出現(xiàn)較多的高次諧波。

2 小波變換原理

目前針對變壓器振動信號的分析多采用快速傅里葉變換(Fast Fourier Transform, FFT)，相比快速傅立葉變換而言，小波變換(Wavelet Transform, WT)本質(zhì)上是一種窗口大小固定不變而其時頻窗可以改變的時頻局部化分析方法，在高頻部分具有較高的時間分辨率和較低的頻率分辨率，在低頻部分具有較高的頻率分辨率和較低的時間分辨率[14-15]，小波變換原理圖如圖1所示。

圖1 小波變換原理圖

正是因為WT具有可變精度特性，從而使其既有足夠的時間分辨率對信號中的短時高頻成分進(jìn)行分析，又有良好的頻率分辨率對信號中的低頻信號進(jìn)行分析。對于非平穩(wěn)信號的分析，WT有著FFT難以替代的優(yōu)勢。WT的基本思想同F(xiàn)FT相似，均為采用一族函數(shù)來表示一個信號或函數(shù)。不同的是FFT變換時將一連續(xù)的、平穩(wěn)的時域信號表示為不同頻率的正弦函數(shù)的線性疊加，而WT則是通過對基波函數(shù)的伸縮和平移來構(gòu)成的。

對于能量有限信號，滿足

其小波變換定義為以函數(shù)簇ψa,b(t)為積分核的積分變換，如下式：

式中：a為伸縮因子；b為平移因子；a-0.5為歸一化常數(shù)，用來保證變換的能量守恒。

函數(shù)簇ψa,b(t)為小基波函數(shù)ψ(t)經(jīng)過伸縮和平移產(chǎn)生的小波函數(shù)簇：

變換前后，信號能量守恒，即

由上式可以看出，WT是一種線性變換，物理含義就是用一簇頻率不同的震蕩函數(shù)ψa,b(t)作為窗口函數(shù)對信號x(t)進(jìn)行掃描和平移，此時WT在某種意義上類似于FFT變換；但不同的是，WT的時域和頻域分辨率與頻率有關(guān)，在高頻段，小波變換能達(dá)到高時域分辨率，而頻域分辨率交叉，在低頻段時相反，這是快速傅里葉變換不具備的。

3 振動信號故障分析

變壓器的振動信號頻率成分十分豐富,適宜采用WT分析振動信號。因為FFT僅得到頻域上一些離散頻點的信息,會導(dǎo)致離散頻點外的信息丟失,而小波變換的結(jié)果是表征各個頻段的時域信息,且時頻窗口大小隨尺度變化,可實現(xiàn)對低頻信號采用大時窗、對高頻信號采用小時窗的分析要求。利用小波理論對變壓器的振動信號進(jìn)行多分辨率分析,并對各個頻段能量進(jìn)行量化處理,可提取出振動信號的頻段-能量對應(yīng)關(guān)系,作為變壓器故障診斷的依據(jù)。變壓器正常運(yùn)行時,其能量大部分集中在100～300 Hz內(nèi)[16-17],當(dāng)變壓器發(fā)生故障或存在潛在故障時,高頻成分增加,導(dǎo)致高頻能量比例增大、主頻段能量減小。將正常運(yùn)行狀態(tài)下的各頻段能量作為參考指紋向量,用實時監(jiān)測到的數(shù)據(jù)同此參考指紋量作比較,當(dāng)主頻段能量小于某一特定值或高頻(>600 Hz)分量能量高于某一特定值時,診斷為變壓器故障。當(dāng)然,此特征值的具體數(shù)值,因變壓器的類型、容量而異,需要通過試驗并結(jié)合現(xiàn)場運(yùn)行經(jīng)驗來確定。

采用JF2001-T壓電式加速度傳感器以及NI9234數(shù)據(jù)采集儀對1臺型號為S11-M-10/10的變壓器樣機(jī)進(jìn)行振動測試，采樣頻率為2 kHz。測試接線原理圖如圖2所示。

圖2 測試接線示意圖

對變壓器箱體不同位置的振動信號進(jìn)行采樣，將采樣數(shù)據(jù)經(jīng)由振動信號采樣設(shè)備進(jìn)行模數(shù)轉(zhuǎn)換，對轉(zhuǎn)換后的數(shù)字信號進(jìn)行WT變換。振動原始信號如圖3所示。

圖3 原始振動信號

由圖3可知，變壓器振動信號具有明顯的周期性，但由于振動信號為較多周期信號的疊加，因此僅依靠原始振動信號幅值變化很難對振動信號進(jìn)行有效識別。對原始振動信號進(jìn)行FFT變化，得到振動信號頻域結(jié)果，如圖4所示。

圖4 FFT變換結(jié)果圖

不難看出，變壓器正常運(yùn)行時，繞組振動信號集中于100 Hz，并含有少量以100 Hz為基頻的的高次諧波，但由于FFT變換無法區(qū)分振動信號的時域差別，會導(dǎo)致對變壓器運(yùn)行的錯誤判斷，因此對Z2采樣信號采用db9小波進(jìn)行4層分解重構(gòu)，并提取重構(gòu)信號能量譜，得到如圖5所示結(jié)果。

由圖5(b)可知，第4層的高頻分量d4與原始信號的總趨勢相吻合，無需再進(jìn)行下一層的小波分解，第4層對應(yīng)于0～250 Hz頻段;其低頻分量a4(對應(yīng)于0～100 Hz)幅度較小，可見振動信號的主頻段為100～200 Hz。

由圖5(a)可知，每進(jìn)行1次小波分解，對應(yīng)的低、高通濾波器實際頻帶減半，即實現(xiàn)2分頻，經(jīng)過4層分解后，其相應(yīng)的低、高頻能量占比分別為EL、EH，1-4層能量占比分別為17.32%、7.21%、11.26%、64.21%。由此可見第4層(0～250 Hz)分解重構(gòu)信號能量占比較大，第4層EH約為61.69，說明變壓器振動信號主要集中于100～200 Hz，而第3層(250～500 Hz)分解重構(gòu)信號能量占比約為11.26%，表明變壓器中頻振動信號占比仍較大，對于第1、2層分解重構(gòu)信號，其能量占比受到外界環(huán)境以及測量誤差的影響，其能量譜結(jié)果對基于振動信號監(jiān)測的變壓器運(yùn)行狀態(tài)來說意義較小，因此可以采用第4層分解重構(gòu)結(jié)果作為變壓器振動信號指紋量，以監(jiān)測變壓器運(yùn)行狀態(tài)并對變壓器運(yùn)行狀態(tài)進(jìn)行快速有效的診斷。進(jìn)一步研究變壓器不同故障下能量譜變化，得到如表1所示結(jié)果。

圖5 小波變換結(jié)果

表1 不同運(yùn)行狀態(tài)能量譜特征

由表1可以看出，當(dāng)變壓器發(fā)生故障時，變壓器振動小波分解重構(gòu)信號能量譜變化較為明顯，以4層小波分解第3、4層分解重構(gòu)信號能量占比為指紋量進(jìn)行分析可以看出：

1)短路故障下低頻信號能量占比增加24.91%，中頻信號能量占比下降16.78%。分析其原因，當(dāng)變壓器短路時，工頻電流激增，使得100 Hz分量下振動加速度明顯增加，使第4層分解重構(gòu)信號能量占比增加。繞組松動故障下低頻信號能量占比下降6.29%，中頻信號能量占比下降36.67%。分析其原因，當(dāng)變壓器繞組松動故障時，振動加速度高次諧波增加，使得低頻信號能量占比下降。三相不平衡運(yùn)行故障下繞組振動加速度能量占比變化情況與繞組松動故障基本一致，其原因也是受到不平衡運(yùn)行諧波的影響，使得低頻信號能量占比下降。

2)通過對振動信號的分解重構(gòu)及能量占比特征值的提取，能夠較好地判斷變壓器的運(yùn)行故障，并及時制定解決措施。

4 結(jié) 語

1)變壓器正常運(yùn)行時，100～200 Hz振動信號分量能量占比較高，且包含以100 Hz為基頻的中高頻振動信號，當(dāng)變壓器故障時，不同頻率信號分量能量占比變化。

2)基于小波分解重構(gòu)的振動信號能量譜能夠有效識別變壓器運(yùn)行故障及內(nèi)部結(jié)構(gòu)問題，為基于振動信號分析的變壓器故障診斷問題提供了輔助決策。