任 巖，黃 今，胡雷鳴，李孝凱，陳紅平

(1.華北水利水電大學(xué)電力學(xué)院，河南 鄭州 450046；2.可再生能源電力技術(shù)湖南省重點(diǎn)實(shí)驗(yàn)室(長沙理工大學(xué))，湖南 長沙 410114；3.江西洪屏抽水蓄能有限公司，江西 宜春 330600；4.國網(wǎng)湖南省電力有限公司水電公司，湖南 長沙 410004)

0 引 言

水輪發(fā)電機(jī)組的檢測(cè)和維修[1-3]是保障機(jī)組運(yùn)行穩(wěn)定的關(guān)鍵措施，在機(jī)組運(yùn)行過程中，背景噪聲強(qiáng)大，且可能含有間歇信號(hào)和脈沖信號(hào)的干擾，為了在檢測(cè)過程中提取有用信號(hào)，對(duì)機(jī)組信號(hào)去噪[4-5]顯得至關(guān)重要。

為了分析處理非平穩(wěn)信號(hào)，一系列信號(hào)分析方法相繼被提出。1998年Huang等人提出了一種全新的信號(hào)時(shí)頻分析方法EMD，EMD[6-7]基于信號(hào)本身極值點(diǎn)分布，將信號(hào)分解為若干個(gè)固有模態(tài)函數(shù)(IMF)和一個(gè)殘余分量，從而對(duì)信號(hào)進(jìn)行全面分析，但由于其缺乏理論結(jié)構(gòu)，因此當(dāng)信號(hào)中存在噪聲，脈沖等異?，F(xiàn)象時(shí)會(huì)引起模態(tài)混疊效應(yīng)等。針對(duì)上述問題，EEMD[8-10]、CEEMD[11-13]等改進(jìn)方法依次被提出，但都依然存在一定程度的模態(tài)混疊問題。而FDM[14-15]作為一種新的時(shí)頻分析方法，以傅里葉分解作為理論基礎(chǔ)，在分解過程中對(duì)極值點(diǎn)沒有求解需求，因此不會(huì)出現(xiàn)經(jīng)驗(yàn)?zāi)B(tài)分解過程中產(chǎn)生的模態(tài)混疊現(xiàn)象，同時(shí)由于并未引進(jìn)白噪聲，因此也不存在白噪聲的殘留。文獻(xiàn)[16]提出了采用傅里葉分解(FDM)和1.5維Teager能量譜的方法來對(duì)強(qiáng)背景噪聲下對(duì)滾動(dòng)軸承微弱故障特征進(jìn)行提取。

排列熵(permutation entropy，PE)[17-19]作為一種檢測(cè)隨機(jī)性時(shí)間序列和動(dòng)力學(xué)突變行為的方法，針對(duì)包含噪聲等干擾信號(hào)的復(fù)雜信號(hào)，有著計(jì)算簡(jiǎn)單、抗噪能力強(qiáng)等優(yōu)點(diǎn)，被廣泛應(yīng)用在故障信號(hào)特征的提取。

因此本文提出一種結(jié)合FDM和排列熵的去噪方法。利用FDM在整個(gè)傅里葉域自適應(yīng)地搜尋解析傅里葉固有頻帶函數(shù)(AFIBFs)，從而將非平穩(wěn)、非線性信號(hào)自適應(yīng)分解為若干個(gè)傅里葉固有頻帶函數(shù)和一個(gè)殘余分量；再通過計(jì)算各頻帶函數(shù)的排列熵，對(duì)比設(shè)定的值，選取含噪少的頻帶函數(shù)，對(duì)其進(jìn)行重構(gòu)。通過仿真和實(shí)例的驗(yàn)證，證明該方法去噪后獲得的分析信號(hào)更趨近于真實(shí)不含噪的信號(hào)，去噪效果較好。

1 基本原理

基于傅里葉分解[14-16](FDM)理論和排列熵[17-19]算法，提出基于FDM和排列熵相結(jié)合的水輪機(jī)信號(hào)去噪方法，具體步驟如下：

(1)利用傅里葉分解(FDM)對(duì)水輪機(jī)含噪信號(hào)進(jìn)行自適應(yīng)分解，分解為若干個(gè)不同頻段的FIBF和一個(gè)殘余分量之和。

(2)確定排列熵時(shí)間延遲t和嵌入維度m，并且根據(jù)公式計(jì)算各個(gè)頻段FIBF的熵值。

(3)通過實(shí)驗(yàn)研究，設(shè)定閾值為0.5作為判斷標(biāo)準(zhǔn)。若FIBF的排列熵量大于0.5，則為包含較多噪聲的信號(hào)，將其舍去；若排列熵值小于0.5，則為符合條件的有用信號(hào)，將其留下。

(4)最后將符合條件的有用信號(hào)進(jìn)行重構(gòu)，即可得到去除背景噪聲后的水輪機(jī)信號(hào)。

2 仿真信號(hào)分析

構(gòu)造水輪機(jī)仿真信號(hào)如式(1)所示，原始加噪信號(hào)和不加噪信號(hào)如圖1、2。

圖1 原始加噪信號(hào)

圖2 原始不加噪信號(hào)

(1)

將原始加噪信號(hào)進(jìn)行FDM分解，按照頻率由高到低依次分解為30個(gè)模態(tài)函數(shù)分量(如圖3所示)。因篇幅限制，現(xiàn)挑選10個(gè)典型的模態(tài)分量波形圖。觀察圖3可發(fā)現(xiàn)信號(hào)被有效地分解，u2、u28、u29特征頻率清晰，無模態(tài)混疊現(xiàn)象，分別對(duì)應(yīng)仿真信號(hào)的x(1)、x(2)、x(3)。之后通過已設(shè)定的延遲時(shí)間和時(shí)間序列，將分解得到的頻帶函數(shù)計(jì)算獲得排列熵，部分分量排列熵如表1所示。

表1 FDM各分量排列熵

圖3 FDM分解

通過多次改變?cè)肼暣笮?duì)仿真信號(hào)進(jìn)行分析，發(fā)現(xiàn)包含較多有用特征的信號(hào)排列熵值普遍偏小，因此設(shè)定閾值R=0.5。對(duì)上述固有頻帶函數(shù)進(jìn)行篩選，發(fā)現(xiàn)u3～u27的排列熵值大于0.5，將其舍去。這不僅符合前面對(duì)FDM分解圖的分析，同時(shí)通過排列熵的計(jì)算進(jìn)行更精確的篩選，從而獲得需要重構(gòu)的固有頻帶函數(shù)。

為進(jìn)行對(duì)比分析，將仿真信號(hào)進(jìn)行EMD分解(如圖4所示)。由圖4可知，原始加噪信號(hào)經(jīng)過分解獲得頻率由高到低的u1～u11模態(tài)分量，其中u1～u3為明顯的噪聲信號(hào)，u4～u6出現(xiàn)一定程度的模態(tài)混疊。之后將各模態(tài)分量同樣進(jìn)行排列熵計(jì)算，部分結(jié)果如表2所示。

圖4 EMD分解

表2 EMD各分量排列熵

由表2可知，u1～u4排列熵大于設(shè)定值0.5，應(yīng)將其舍去。對(duì)比圖4中的分析，舍去了可能為有用信號(hào)的u4，并且選取的重構(gòu)信號(hào)中u5、u6存在一定的模態(tài)混疊現(xiàn)象，可能會(huì)影響對(duì)信號(hào)的去噪效果。

接著根據(jù)FDM-排列熵、EMD-排列熵的處理結(jié)果，選取含噪少的頻帶函數(shù)，對(duì)仿真信號(hào)進(jìn)行重構(gòu)，獲得去噪后的信號(hào)圖，如圖5、6所示。對(duì)比圖2與圖5，發(fā)現(xiàn)不管是在波形的走向還是特征頻率的體現(xiàn)都和原波形圖有著極高的相似度，說明失真度很??；同時(shí)波形非常光滑，表明去噪效果很好。對(duì)比圖2和圖6，一方面在波形上可以看出明顯的失真現(xiàn)象，重構(gòu)后的波形圖遺失了部分有用信息，同時(shí)在端點(diǎn)處還出現(xiàn)了一定程度的模態(tài)混疊；另一方面波形雖較圖1光滑度顯著提高，但是相比較FDM-排列熵處理后的波形(圖5)還是可以看到一定程度的粗糙。

圖5 FDM-排列熵去噪信號(hào)

圖6 EMD-排列熵去噪信號(hào)

為了能夠更科學(xué)地體現(xiàn)2種方法去噪效果的好壞，本文引進(jìn)了相關(guān)系數(shù)和信噪比。通過計(jì)算，F(xiàn)DM-排列熵相關(guān)系數(shù)為0.997 1；EMD-排列熵相關(guān)系數(shù)為0.972 1，低于FDM-排列熵，與原信號(hào)相似程度低，重構(gòu)后的波形缺失部分有用信號(hào)。FDM-排列熵信噪比為13.67；EMD-排列熵信噪比為11.89，信噪比也低于FDM-排列熵，說明EMD-排列熵相對(duì)來說去噪效果較差，符合之前對(duì)波形圖的分析。

3 實(shí)例分析

本文以國內(nèi)某水電站為例展開機(jī)組試驗(yàn)。機(jī)組參數(shù)如下：水輪機(jī)型號(hào)為HL-702-LJ-410，水輪機(jī)設(shè)計(jì)出力為63 MW，額定轉(zhuǎn)速為136 r/min。在水輪機(jī)上導(dǎo)X相所采擺度數(shù)據(jù)，采樣頻率為500 Hz，數(shù)據(jù)長度為100 000。截取其中一段進(jìn)行分析，信號(hào)波形如圖7所示。

圖7 實(shí)測(cè)水輪機(jī)信號(hào)

采用FDM對(duì)水輪機(jī)采集信號(hào)進(jìn)行處理，其中表3

表3 水輪機(jī)采集信號(hào)各模態(tài)函數(shù)排列熵

為水輪機(jī)采集信號(hào)部分頻帶函數(shù)排列熵，選擇u1、u26、u27、u28、u29進(jìn)行重構(gòu)。圖8為FDM重構(gòu)后的波形。由圖8可知重構(gòu)后的波形光滑，無明顯的毛刺，并且各個(gè)特征信號(hào)均無缺失，因此可說明達(dá)到了較好的去噪效果。

圖8 FDM-排列熵去噪后信號(hào)

4 結(jié) 論

(1)通過對(duì)EMD-排列熵和FDM-排列熵分別進(jìn)行仿真信號(hào)和水輪機(jī)實(shí)例分析，發(fā)現(xiàn)EMD-排列熵對(duì)振動(dòng)信號(hào)去噪過程中會(huì)出現(xiàn)模態(tài)混疊現(xiàn)象，F(xiàn)DM-排列熵并未出現(xiàn)。

(2)將EMD-排列熵和FDM-排列熵重構(gòu)后信號(hào)以相關(guān)系數(shù)和信噪比作為評(píng)價(jià)指標(biāo)進(jìn)行對(duì)比，發(fā)現(xiàn)FDM-排列熵較EMD-排列熵相關(guān)系數(shù)、信噪比更高，波形更光滑，說明該方法去噪效果較好。