電力擾動信號的小波形態(tài)學(xué)去噪方法

陳曉娟，李松寒，隋吉生，王 麗

(1.東北電力大學(xué)信息工程學(xué)院，吉林 吉林132012;2.吉林省電力有限公司，吉林 長春130021)

0 前言

隨著人們生活水平的不斷提高，對電能質(zhì)量也提出了越來越高的要求。由于大容量電力設(shè)備和用電整流設(shè)備的投入以及其他非線性負荷的影響，不可避免地會發(fā)生各種電能質(zhì)量擾動，其主要表現(xiàn)形式為電壓凹陷、降落和瞬時擾動等。任何擾動都會對電能質(zhì)量造成破壞。因此，如何找出擾動現(xiàn)象的規(guī)律，并及時地采取相應(yīng)措施以避免擾動及其帶來的損失日益成為提高電能質(zhì)量中的重要課題。為了更好地監(jiān)測電力系統(tǒng)的運行狀態(tài)，大量的電力監(jiān)測裝置廣泛應(yīng)用在電力系統(tǒng)中。但采集到的信號通常含有噪聲，為了精確地反映被采集信號量的特征，必須對其進行去噪處理，通常采用小波變換的方法［1］。但有些擾動信號的小波變換系數(shù)值可能與噪聲的小波變換系數(shù)值相近，且低于閾值被當作噪聲去掉，導(dǎo)致信號部分失真。而基于形態(tài)學(xué)的方法［2］，由于其計算速度快、實時性強的優(yōu)點，應(yīng)用于信號消噪也取得了較好的效果，但并未應(yīng)用在電力擾動信號的處理。本文就如何取得更好的去噪效果進行了討論，并對以上幾種方法進行了分析和仿真，并創(chuàng)新地提出了小波形態(tài)學(xué)相結(jié)合的去噪方法，該算法的實質(zhì)是采用形態(tài)學(xué)的方法消除大部分的噪聲，然后再利用特殊處理后的閾值經(jīng)過小波變換對其做平滑處理，從而使信號的失真度大大降低，而且能夠有效消除噪聲。

1 小波去噪和形態(tài)學(xué)原理

1.1 小波去噪原理

在實際電力系統(tǒng)發(fā)生擾動時，通過監(jiān)測裝置采集到的擾動信號都是含噪的非平穩(wěn)信號，信號在任一時刻附近表現(xiàn)出的頻率特性都很重要，僅從時域或頻域上來分析是不夠的。由于傳統(tǒng)去噪方法的種種限制，小波理論在信號去噪領(lǐng)域逐漸受到了許多學(xué)者的重視。目前，小波消噪的方法主要分為小波分解與重構(gòu)法、非線性小波變換閾值法、平移不變量法和小波變換模極大值法。其中，非線性小波變換閾值法應(yīng)用范圍比較廣泛，通常適用于處理的信號中含有隨機白噪聲的情況，盡可能多的保留原始信號的重要信息，也能對原始信號進行很好的估計，并且計算速度也比較快。因此，本文以此方法為例來說明小波去噪的基本原理。其具體的處理步驟如下［1］:

(Ⅰ)選擇合適的小波基函數(shù)，運用小波變換的方法對含噪的信號進行N層分解(小波分解的層數(shù)需根據(jù)噪聲的相關(guān)類型進行確定)。

(Ⅱ)小波分解高頻系數(shù)的閾值量化。根據(jù)相關(guān)的閾值選取規(guī)則，對以上進行小波分解尺度下的高頻系數(shù)進行閾值量化處理。閾值量化處理大體分為軟閾值量化處理方法和硬閾值量化處理方法兩種。

(Ⅲ)對經(jīng)過上述處理過后的底層低頻系數(shù)和各層的高頻系數(shù)進行一維小波重構(gòu)。

在以上的這3個步驟中，最關(guān)鍵的是小波基的選取以及閾值的選取和量化，這在某種程度上關(guān)系到信號去噪質(zhì)量的好壞。

1.2 數(shù)學(xué)形態(tài)學(xué)的基本原理

數(shù)學(xué)形態(tài)學(xué)是一種非常有效的非線性濾波技術(shù)。運用形態(tài)學(xué)在進行信號處理時，需不斷移動結(jié)構(gòu)元素以便提取信號中的有用信息。與傳統(tǒng)濾波技術(shù)相比，該方法易于實現(xiàn)且計算時間少。目前，數(shù)學(xué)形態(tài)學(xué)主要應(yīng)用在電力系統(tǒng)中的繼電保護［3］和信號去噪［4］等方面。由于電力系統(tǒng)信號分析和處理中所涉及到的信號通常是一維信號，本文就一維離散情況下的多值形態(tài)變換(包括膨脹運算、腐蝕運算以及此兩種運算級聯(lián)組合所構(gòu)成的運算)進行分析。

設(shè)輸入信號序列和結(jié)構(gòu)元素分別由f(n)和g(n)代表，則f(n)關(guān)于g(n)的膨脹和腐蝕具體計算公式分別如下:

式中，⊕、Θ分別為膨脹、腐蝕的運算符號。

形態(tài)開運算和閉運算是基于膨脹和腐蝕運算的基礎(chǔ)上進行展開的，其具體計算公式分別如下:

式中，°、●分別為形態(tài)開運算和形態(tài)閉運算的運算符號。

在大多數(shù)形態(tài)濾波器的構(gòu)造中，開運算通常被用于正脈沖噪聲的去除，而閉運算則是常被用于負脈沖噪聲消除。考慮到開運算和閉運算去除噪聲的類型，Maragos則采用基于開閉運算的級聯(lián)組合方式，即形態(tài)開－閉(OC)和閉－開(CO)濾波器。其定義式如下:

由于開運算的反擴展性和閉運算的擴展性會造成經(jīng)過開－閉濾波器處理后的信號幅度偏小，而經(jīng)過閉－開濾波器處理后的信號幅度偏大，所以單獨使用它們并不能獲得良好的濾波效果。為了獲得更好的去噪效果，通常采用開－閉濾波器和閉－開濾波器的平均組合方式，其定義式如下:

經(jīng)過形態(tài)濾波器處理效果的好壞不僅取決于上述形態(tài)運算的方式，而且還取決于結(jié)構(gòu)元素的形狀和長度［4］。經(jīng)過大量的現(xiàn)場實驗證明:平滑的半圓形結(jié)構(gòu)元素通常被用來抑制電力信號中常見的白噪聲;三角形的結(jié)構(gòu)元素則常被用來進行含有脈沖噪聲的信號處理［3］。另外，所選取的結(jié)構(gòu)元素長度過大或者過小都不能達到很好去噪效果，同時也會破壞原信號的結(jié)構(gòu)。文獻［5］指出采用長度為信號周期長度的1/50～1/20的結(jié)構(gòu)元素是比較適合的。

2 小波形態(tài)學(xué)融合算法

小波變換在非平穩(wěn)信號的處理上效果顯著，形態(tài)學(xué)則在時域的處理上具有很好的效果且較為簡單。本文將此兩種算法融合進行消噪，算法的具體步驟如下:

(Ⅰ)選擇合適的小波基函數(shù)和適當?shù)某叨龋瑢⒑性肼暤臄_動信號進行N層小波分解(本文是進行3層小波分解)，提取各層的小波系數(shù)。

(Ⅱ)根據(jù)相應(yīng)結(jié)構(gòu)元素選擇規(guī)則［6］對各尺度下的小波系數(shù)進行選取結(jié)構(gòu)元素，然后采用交替混合型形態(tài)濾波器進行濾波，進而得到新的小波系數(shù)。

(Ⅲ)根據(jù)相關(guān)閾值求解方法［5，7］進行閾值的求解，再將此閾值與系數(shù)k(0＜k＜1)進行相乘。然后，再利用軟閾值方法對新的每一層的小波系數(shù)進行閾值量化處理。其中系數(shù)k的值是采用自適應(yīng)方法求得，其具體計算方法如下:

3 仿真及結(jié)果分析

如圖1所示的擾動信號為例，然后對其疊加白噪聲和隨機背景噪聲。從圖1中可以看出:模擬出的擾動信號已經(jīng)被噪聲完全給淹沒了，不能確定其中是否含有擾動信號，更不能直觀地看出哪些是有用的信息。

3.1 形態(tài)學(xué)的信號去噪

本文運用形態(tài)學(xué)的信號去噪處理，所選用的結(jié)構(gòu)元素為扁平結(jié)構(gòu)元素、長度為13交替濾波器(這樣可以解決信號處理后的單向偏移問題)。從其處理結(jié)果(見圖2a)所示，運用數(shù)學(xué)形態(tài)學(xué)處理過后的信號存在很多毛刺，而且丟失了很重要的原信號信息。

3.2 小波方法的信號去噪

考慮到中國配電網(wǎng)線路上的實際噪聲情況，本文應(yīng)用非線性小波變換閾值法［1］對擾動信號進行處理。根據(jù)電力系統(tǒng)暫態(tài)信號分析中小波基的選擇相關(guān)原則［8］，選取的小波基是與所模擬擾動信號波形相似的db5小波對含噪的擾動信號進行4層分解。去噪閾值選取按照‘heursure’和‘minimaxi’兩種閾值選擇規(guī)則進行設(shè)定，其仿真結(jié)果如圖2b、圖2c所示。此兩種方法在去除噪聲的同時還去除掉了原信號的一些特征，不利于對擾動信號進行后期的處理。

圖1 原始擾動信號及含噪擾動信號

3.3 小波形態(tài)學(xué)融合的信號去噪

本文選擇與擾動信號相似的小波基函數(shù)db5并對含噪的信號進行3層分解，然后對各尺度下的小波系數(shù)選取長度為10、13、11的扁平結(jié)構(gòu)元素，采用交替混合型形態(tài)學(xué)濾波器進行濾波，得到新的小波系數(shù)。利用軟閾值對新的每一層的小波系數(shù)進行閾值處理，然后運用閾值處理過后的小波系數(shù)來重構(gòu)消噪后的信號，其結(jié)果如圖2d所示。從圖2中可以看出:該方法不僅有效地去除了噪聲，而且還保留信號在突變位置的相關(guān)信息，失真度很小。

以上只是通過主觀的方法比較，可看出基于小波－形態(tài)學(xué)融合的去噪方法明顯好于其他的去噪方法。但是為了更好地從原理上說明，采用客觀的誤差評價方法，其計算公式如下:

其中，Y(i)與X(i)為原始信號和消噪后的信號。

利用公式(9)對這3種方法的誤差進行估計，計算結(jié)果如表1所示。

圖2 運用各種方法消噪后信號波形

表1 各種方法的評價指標

從表1可以看出:基于小波－形態(tài)學(xué)融合的去噪方法優(yōu)于其他兩種方法。

4 結(jié)論

在電力系統(tǒng)暫態(tài)信號分析應(yīng)用中，必須根據(jù)所要求的精度、速度等相關(guān)因素選取合適的信號處理方法以達到滿意的處理效果。單獨使用小波或形態(tài)學(xué)進行去噪，都會造成原信號的部分失真。本文將此兩種方法進行融合以便達到很好的去噪效果。仿真結(jié)果證明:融合的去噪方法能很好地保留原信號的重要信息，且信號的失真度非常小，有很強的實用性。

［1］何慧，顧潔.電能質(zhì)量信號的小波閾值去噪［J］.電力系統(tǒng)及其自動化學(xué)報，2010，22(2):103－108.

［2］尹文琴，劉前進.數(shù)學(xué)形態(tài)學(xué)在電力系統(tǒng)中的應(yīng)用綜述［J］.繼電器，2007，35(19):76－83.

［3］李庚銀，羅艷，周明，等.基于數(shù)學(xué)形態(tài)學(xué)和網(wǎng)格分形的電能質(zhì)量擾動檢測及定位［J］.中國電機工程學(xué)報，2006，26(3):25－30.

［4］張文斌，楊辰龍，周曉軍.形態(tài)濾波方法在振動信號降噪中的應(yīng)用［J］.浙江大學(xué)學(xué)報，2009，43(11):2096－2099.

［5］歐陽森，宋政湘，陳德桂，等.小波軟閾值去噪技術(shù)在電能質(zhì)量檢測中的應(yīng)用［J］.電力系統(tǒng)自動化，2002，26(19):56－60.

［6］趙昭，劉利林，張承學(xué)，等.形態(tài)學(xué)濾波器結(jié)構(gòu)元素選取原則研究與分析［J］.電力系統(tǒng)保護與控制，2009，37(14):21－25.

［7］張磊，潘泉，張洪才，等.小波域濾波閾值參數(shù) c的選?。跩］.電子學(xué)報，2001，29(3):400－402.

［8］何正友，錢清泉.電力系統(tǒng)暫態(tài)信號分析中小波基的選擇原則［J］.電力系統(tǒng)自動化設(shè)備，2003，25(10):45－48.