ＥＭＤ方法在局部放電超聲信號提取中的應(yīng)用

盧　麗

摘 要：超聲波定位法是一種對變壓器內(nèi)部不同部位放電點進行檢測的行之有效的方法。由于變壓器所處現(xiàn)場存在大量電磁干擾，采集到的信號包含大量噪聲，因此必須做相應(yīng)的預處理。變壓器局部放電超聲波定位中，首要任務(wù)是對超聲波信號的提取，其準確性將直接影響到定位的實現(xiàn)。EMD方法可以將非線性、非平穩(wěn)信號分解為不同頻率成分的內(nèi)稟模態(tài)函數(shù)，從而提取超聲波信號。以陜西省高壓電氣設(shè)備局部放電定位項目為背景，詳細分析了EMD方法的原理、實現(xiàn)過程，并用實例測試仿真，驗證了該方法的有效性與準確性。

關(guān)鍵詞：EMD；內(nèi)稟模態(tài)函數(shù)；局部放電；超聲檢測

中圖分類號：TP23文獻標識碼：B

文章編號：1004-373X(2009)03-137-03

Application of EMD Method in Extracting Ultrasonic Signals of Transformer Partial Discharge

LU Li

(Shaanxi Regional Electric Power Group,Xi′an,710061,China)

Abstract：The ultrasonic wave locating is an effective metheod on the examination to the transformer interior electric discharge spot of different parts.The transformer at the scene exsits massive electromagnetic interference and the gathered signal contains massive noise,so the signal must make corresponding pretreatment.In the ultrasonic location of transformer partial discharge,the most important is to extract the ultrasonic signals.Its precision directly affects the accuracy of the location.The Empirical Mode Decomposition(EMD) method can decompose the nonlinear and non-stationary signals to Intrinsic Mode Functions(IMF) with different frequency,then the ultrasonic signals are extracted.Based on the Shaanxi Project of High Vol-tage Instruments,the principle and realization process of EMD method are analysed,the correctness and stability are verified.

Keywords：EMD；IMF；partial discharge;ultrasonic examination

0 引 言

變壓器內(nèi)部局部放電在線監(jiān)測和定位對變壓器的故障檢修十分重要。超聲波定位法是一種對變壓器內(nèi)部不同部位放電點進行檢測的行之有效的方法。由于變壓器所處現(xiàn)場存在大量電磁干擾，采集到的信號包含大量噪聲，因此必須做相應(yīng)的預處理。

基于經(jīng)驗的模態(tài)分解(Empirical Mode Decomposition,EMD)方法可以有效地對非平穩(wěn)信號中的各頻率成分進行分離，提取出超聲波信號。其主要過程是采用EMD方法將信號分解為若干個內(nèi)稟模態(tài)函數(shù)(Intrinsic Mode Function，IMF)分量之和，然后對每個IMF分量進行時頻率分析。該方法實現(xiàn)速度快，提取波形精度高，對非平穩(wěn)、非線性信號具有良好的時頻聚集性。本文對EMD方法的原理、實現(xiàn)過程進行了詳細分析，并通過仿真數(shù)據(jù)驗證了算法的有效性與準確性。

1 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) 在任意時刻，有局部極大值點形成的上包絡(luò)線和由局部極小值點形成的下包絡(luò)線的平均值為零，即上、下包絡(luò)線相對于時間軸局部對稱。

內(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ù)量的極值點和過零點，在相鄰的兩個過零點之間只有一個極值點，而且上、下包絡(luò)線關(guān)于時間軸局部對稱，任何兩個模態(tài)之間是相互獨立的；任何時候，一個信號都可以包含許多內(nèi)稟模態(tài)函數(shù)，如果模態(tài)函數(shù)相互重疊，便形成復雜信號。在此假設(shè)的基礎(chǔ)上，可以采用EMD方法通過下面的步驟對任何信號x(t)進行分解：

(1) 確定信號的局部極值點，用三次樣條插值［4］將所有的局部極大值點連接形成上包絡(luò)線。

(2) 用三次樣條插值將所有的局部極小值點連接形成下包絡(luò)線，上、下包絡(luò)線應(yīng)該包絡(luò)所有的數(shù)據(jù)點。

(3) 上、下包絡(luò)線的平均值記為m1，求出:

x(t)－m1=h1

(1)

理想地，如果h1是一個IMF，那么h1就是x(t)的第1個IMF分量。

(4) 如果h1不滿足IMF的條件，把h1作為原始數(shù)據(jù)，重復步驟(1)~(3)，得到上、下包絡(luò)線的平均值m11，再判斷h11=h1－m11是否滿足IMF的條件，如不滿足，則重循環(huán)k次，得到h1(k－1)－m1k=h1k，使得h1k滿足IMF的條件。記c1=h1k，則c1為信號x(t)的第1個滿足IMF條件的分量。

(5) 將c1從x(t)中分離出來，得到:

r1=x(t)－c1

(2)

將r1作為原始數(shù)據(jù)重復步驟(1)~(4)，得到x(t)的第2個滿足IMF條件的分量c2，重復循環(huán)n次，得到信號x(t)的n個滿足IMF條件的分量。這樣就有:

r1－c2=r2髍n－1－cn=rn

(3)

當rn成為一個單調(diào)函數(shù)不能再從中提取滿足IMF條件的分量時，循環(huán)結(jié)束。這樣由式(2)和式(3)得到:

x(t)=∑ni=1ci+rn

(4)

式中，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分量ci，取c1～c6分別如圖2～圖7所示，圖8中r6為殘余分量，其中，橫坐標為采樣點數(shù)，縱坐標為信號幅值。

圖1 初始信號x(t)

圖2 第1個IMF分量c1

圖3 第2個IMF分量c2

圖4 第3個IMF分量c3

圖5 第4個IMF分量c4

圖6 第5個IMF分量c5

圖7 第6個IMF分量c6

圖8 殘余分量r6

從圖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è)備局部放電定位項目為背景，通過理論分析和實驗仿真，充分說明了在對采集到的初始信號進行超聲波信號的提取應(yīng)用中，EMD方法具有比較好的有效性和準確性,可以取得較高質(zhì)量的超聲波信號，為后期的定位計算做足準備。

參考文獻

［1］Huang N E,Shen Z,Long S R,et al.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Non-stationary Time Series Analysis.Proc.R.Soc.Lond.A,1998,454:903-995.

［2］Huang N E,Shen Z,Long S R.A New View of Nonlinear Water Waves:The Hilbert Spectrum.Annu.Rev.Fluid Mech.,1999,31:417-457.

［3］于德介,程軍圣,楊宇.機械故障診斷的Hilbert-Huang變換方法.北京：科學出版社,2006.

［4］David Kincaid,Ward Cheney.Numerical Analysis:Mathematics of Scientific Computing.3版.北京:機械工業(yè)出版社,2005.

［5］李海濤,鄧櫻.Matlab程序設(shè)計教程.北京：高等教育出版社，2002.

［6］何文林,趙壽生,劉浩軍.大容量電力變壓器局部放電定位試驗研究［J］.電力設(shè)備,2006,5(10):35-37.

［7］Judd M D，Cleary G P，Bennoch C J，et al.Power Transformer Monitoring Using UHF Sensors：Site Trials［A］.Proceedings of the 2006 IEEE International Symposium on Electrical Insulation［C］.Boston,2006:145-149.

［8］Raja K，Devaux F，Lelaidier S.Recognition of Discharge Sources Using UHF PD Signatures［J］.IEEE Electrical Insulation Magazine，2007,18(5)：8-14.

［9］葛為民.變壓器局部放電綜合判斷分析［J］.高壓電器,2006,42(2)：144-148.

［10］龔慶武，劉美觀，左克鋒，等.GPS同步采樣裝置中防止干擾GPS秒脈沖信號的措施［J］.電力系統(tǒng)自動化，2006,24(1):45-47.

作者簡介

盧 麗 女，1975年出生，陜西米脂人，碩士研究生。主要研究方向為電力電氣工程。

注：本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文。