O 引 言
變壓器內(nèi)部局部放電在線監(jiān)測和定位對變壓器的故障檢修十分重要。超聲波定位法是一種對變壓器內(nèi)部不同部位放電點進行檢測的行之有效的方法。由于變壓器所處現(xiàn)場存在大量電磁干擾，采集到的信號包含大量噪聲，因此必須做相應的預處理。
基于經(jīng)驗的模態(tài)分解(Empirical Mode Decomposi—tion，EMD)方法可以有效地對非平穩(wěn)信號中的各頻率成分進行分離，提取出超聲波信號。其主要過程是采用EMD方法將信號分解為若干個內(nèi)稟模態(tài)函數(shù)(Intrin—sic Mode Function，IMF)分量之和，然后對每個IMF分量進行時頻率分析。該方法實現(xiàn)速度快，提取波形精度高，對非平穩(wěn)、非線性信號具有良好的時頻聚集性。本文對EMD方法的原理、實現(xiàn)過程進行了詳細分析，并通過仿真數(shù)據(jù)驗證了算法的有效性與準確性。
l EMD方法
從物理學的角度來看，信號可分為單分量和多分量信號兩大類。單分量信號在任意時刻都只有一個頻率，稱為信號的瞬時頻率。多分量信號則在某些時刻具有各自的瞬時頻率。瞬時頻率可以很好地表示信號頻率隨時間變化的情況。EMD方法通過將信號分解為若干個IMF分量之和，從而分析各分量的時頻特性。
1．1 內(nèi)稟模態(tài)函數(shù)IMF
在EMD變換中，為了計算瞬時頻率，定義了內(nèi)稟模態(tài)函數(shù)。一個內(nèi)稟模態(tài)函數(shù)必須滿足下面兩個條件：
(1)在整個數(shù)據(jù)段內(nèi)，極值點的個數(shù)和過零點的個數(shù)必須相等或相差最多不能超過一個；
(2)在任意時刻，有局部極大值點形成的上包絡線和由局部極小值點形成的下包絡線的平均值為零，即上、下包絡線相對于時間軸局部對稱。
內(nèi)稟模態(tài)函數(shù)反映信號內(nèi)部固有的波動性，在它的每一個周期上，僅僅包含一個波動模態(tài)，不存在多個波動模態(tài)混疊的現(xiàn)象。
1．2 EMD方法——“篩分”過程
對于內(nèi)稟模態(tài)函數(shù)，可以用Hilbert變換構(gòu)造解析信號，然后求出瞬時頻率。而對不滿足內(nèi)稟模態(tài)函數(shù)條件的普通信號，先要采用EMD方法將其分解。這個分解過程基于一個基本的假設(shè)：任何復雜的信號都是由一些不同的內(nèi)稟模態(tài)函數(shù)組成，每一內(nèi)稟模態(tài)函數(shù)不論是線性或是非線性、非平穩(wěn)的，都具有相同數(shù)量的極值點和過零點，在相鄰的兩個過零點之間只有一個極值點，而且上、下包絡線關(guān)于時間軸局部對稱，任何兩個模態(tài)之間是相互獨立的；任何時候，一個信號都可以包含許多內(nèi)稟模態(tài)函數(shù)，如果模態(tài)函數(shù)相互重疊，便形成復雜信號。在此假設(shè)的基礎(chǔ)上，可以采用EMD方法通過下面的步驟對任何信號x(t)進行分解：
(1)確定信號的局部極值點，用三次樣條插值將所有的局部極大值點連接形成上包絡線。
(2)用三次樣條插值將所有的局部極小值點連接形成下包絡線，上、下包絡線應該包絡所有的數(shù)據(jù)點。
(3)上、下包絡線的平均值記為ml，求出：
理想地，如果h1是一個IMF，那么h1就是x(t)的第1個IMF分量。
(4)如果h1不滿足IMF的條件，把h1作為原始數(shù)據(jù)，重復步驟(1)～(3)，得到上、下包絡線的平均值m11，再判斷h11=h1一m11是否滿足IMF的條件，如不滿足，則重循環(huán)k次，得到hi(k-1)-m1k=h1k，使得h1k滿足IMF的條件。記c1=h1k，則c1為信號x(t)的第1個滿足IMF條件的分量。
(5)將c1從x(t)中分離出來，得到：
將r1作為原始數(shù)據(jù)重復步驟(1)～(4)，得到x(t)的第2個滿足IMF條件的分量c2，重復循環(huán)n次，得到信號x(t)的n個滿足IMF條件的分量。這樣就有：
當rn成為一個單調(diào)函數(shù)不能再從中提取滿足IMF條件的分量時，循環(huán)結(jié)束。這樣由式(2)和式(3)得到：
式中，rn稱為殘余函數(shù)，代表信號的平均趨勢。
EMD的分解過程其實是一個“篩分”過程，在“篩分”的過程中，不僅消除了模態(tài)波形的疊加，而且使波形輪廓更加對稱。EMD方法從特征時間尺度出發(fā)，首先把信號中特征時間尺度最小的模態(tài)分離出來，然后分離特征時間尺度較大的模態(tài)函數(shù)，最后分離特征時間尺度最大的分量，因此可以把EMD方法看成是一組高通濾波器。
2．實驗仿真
在變壓器局部放電在線檢測項目中，聲源信號頻率為40~300 kHz，采樣頻率為5 MHz，采樣點n=2 048，噪聲為高斯白噪聲，在頻率為50 Hz的工頻信號上，采集到的初始信號x(t)用Matlab模擬，如圖1所示，對其用EMD方法進行分解，得出的各個IMF分量c1，取c1～c6分別如圖2～圖7所示，圖8中r6為殘余分量，其中，橫坐標為采樣點數(shù)，縱坐標為信號幅值。
從圖1～圖8可以看出，對初始信號(圖1)很難確定出其有效成分，更不能直接用來計算時延。經(jīng)過EMD分解后，總共得到了10個IMF分量，這里取了前6個，其余的均作為殘余函數(shù)處理。c1～c3基本形狀相似，對照實驗數(shù)據(jù)，可以解釋為白噪聲；c4可以解釋為工頻；而c5則可以清楚地表現(xiàn)出超聲波信號；對c6意義不大，而殘余分量r6范圍也可以接受。這樣一來，對所采集到的信號經(jīng)過預處理，可以直接提取出最有價值的超聲波信號部分，為后面的定位計算做了很好的鋪墊。
3 結(jié) 語
本文以陜西省高壓電氣設(shè)備局部放電定位項目為背景，通過理論分析和實驗仿真，充分說明了在對采集到的初始信號進行超聲波信號的提取應用中，EMD方法具有比較好的有效性和準確性，可以取得較高質(zhì)量的超聲波信號，為后期的定位計算做足準備。
下載文檔
