唐 鋒（雅礱江流域水電開發(fā)有限公司，成都 610000）

基于PEEMD的水輪機(jī)尾水管壓力脈動(dòng)分析

唐 鋒
（雅礱江流域水電開發(fā)有限公司，成都 610000）

針對(duì)水輪機(jī)尾水管壓力脈動(dòng)信號(hào)表征出的非平穩(wěn)性和時(shí)變特性，本文采用一種部分集成經(jīng)驗(yàn)?zāi)B(tài)分解（Partly Ensemble Empirical Mode Decomposition，PEEMD）法對(duì)其迚行分析處理。該斱法利用排列熵對(duì)信號(hào)隨機(jī)性的敏感特性，幵成對(duì)添加幅值符號(hào)相反的白噪聲至原始信號(hào)，有效解決了集成經(jīng)驗(yàn)?zāi)B(tài)分解（Ensemble Empirical Mode Decomposition，EEMD）斱法存在的高時(shí)耗、重構(gòu)殘?jiān)胝`差問題。將該斱法應(yīng)用于水輪機(jī)尾水管壓力脈動(dòng)信號(hào)特征提取中，表明了PEEMD斱法的可行性與有效性，為實(shí)現(xiàn)水輪發(fā)電機(jī)組在線監(jiān)測(cè)與故障診斷提供了一種快速有效的信號(hào)處理與分析斱法。

水輪機(jī)；壓力脈動(dòng)信號(hào)；集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）；部分集成經(jīng)驗(yàn)?zāi)B(tài)分解（PEEMD）；信號(hào)處理

0 引言

隨著我國(guó)水利事業(yè)的興起，水輪發(fā)電機(jī)組正向著大型化、復(fù)雜化、高速化斱向發(fā)展，機(jī)組面臨的振動(dòng)安全問題日益凸顯。水輪發(fā)電機(jī)組是一個(gè)極其復(fù)雜的耦合旋轉(zhuǎn)機(jī)械系統(tǒng)，其運(yùn)行過程受到水力因素、機(jī)械因素、電磁因素三斱面的影響[1]。

水作為機(jī)組運(yùn)行過程中的直接工作介質(zhì)，水力因素導(dǎo)致的機(jī)組振動(dòng)在電站長(zhǎng)期運(yùn)行中是不可避免的問題。其中，尾水管渦帶是引起機(jī)組振動(dòng)的最主要水力原因乊一?；炝魇剿啓C(jī)在偏離最優(yōu)工況下運(yùn)行時(shí)，其轉(zhuǎn)輪出口處的旋轉(zhuǎn)分速度會(huì)在尾水管中形成低頻渦帶而產(chǎn)生壓力脈動(dòng)。這種低頻渦帶不僅會(huì)導(dǎo)致尾水管道自身的強(qiáng)烈振動(dòng)，還會(huì)引發(fā)機(jī)組軸系的周期性擺動(dòng)，甚至威脅電網(wǎng)安全[2]。及時(shí)有效地監(jiān)測(cè)機(jī)組尾水管渦帶狀況，對(duì)維護(hù)電站安全穩(wěn)定運(yùn)行、保障電網(wǎng)供電安全具有十分重要的意義。

由于機(jī)組自身結(jié)構(gòu)及運(yùn)行環(huán)境的復(fù)雜性，采集得到的尾水管壓力脈動(dòng)信號(hào)往往表征出強(qiáng)烈的非平穩(wěn)性和時(shí)變特性。傳統(tǒng)基于傅里葉變換的信號(hào)處理斱法（傅快速里葉變換、小波變換、魏格納變換）已不再適用。近年來提出的經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）在信號(hào)處理領(lǐng)域，相比小波變換表現(xiàn)出獨(dú)特的自適應(yīng)性能力，在機(jī)械故障診斷斱面得到了廣泛地應(yīng)用[3]。為了解決EMD存在的模態(tài)混疊問題，Wu和Huang于2005年創(chuàng)造性地提出了集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD），該斱法利用白噪聲時(shí)間尺度均勻性及其零均值特性有效地解決了EMD存在的模態(tài)混疊問題[4]。然而，EEMD斱法引入白噪聲的同時(shí)，也產(chǎn)生了一系列問題。包括兩關(guān)鍵參數(shù)（添加噪聲幅值、集成試驗(yàn)次數(shù)）的選擇、高時(shí)耗以及重構(gòu)信號(hào)中殘余噪聲的問題。針對(duì)參數(shù)選擇問題，Lei指出目前幵沒有合適的比較統(tǒng)一的參數(shù)選擇斱法[5]。在Wu和Huang提出EEMD斱法的同時(shí)給出的經(jīng)驗(yàn)值法，目前應(yīng)用效果不錯(cuò)[4]。針對(duì)EEMD存在的高時(shí)耗以及殘余噪聲問題，鄭提出了部分集成經(jīng)驗(yàn)?zāi)B(tài)分解（Partly Ensemble Empirical Mode Decomposition, PEEMD）法[6]。該斱法利用排列熵對(duì)信號(hào)隨機(jī)性的敏感特性，獲取原始信號(hào)包含的異常成分（噪聲成分和間歇成分），然后從原始信號(hào)中剔除異常成分，最后對(duì)剩余成分直接迚行EMD分解。該斱法避免了不必要的集成試驗(yàn)次數(shù)，大大縮短了計(jì)算時(shí)間。另外，所添加的噪聲為幅值符號(hào)相反的正負(fù)噪聲，通過多次集成平均，很好地消除了重構(gòu)信號(hào)中的殘余噪聲誤差。最終將該斱法應(yīng)用于水輪機(jī)尾水管壓力脈動(dòng)信號(hào)分析與特征提取中，實(shí)例應(yīng)用表明了該斱法的有效性，為水輪發(fā)電機(jī)組在線監(jiān)測(cè)與故障診斷提供了一種有效的信號(hào)處理與分析斱法。

1 EEMD及其存在的問題

1.1 集成經(jīng)驗(yàn)?zāi)B(tài)分解法（EEMD）

EMD斱法被認(rèn)為是自2000年來以傅里葉變換為基礎(chǔ)的線性和穩(wěn)態(tài)頻譜分析的一個(gè)重大突破[7]，它能夠依據(jù)信號(hào)自身的局部時(shí)間尺度將原始信號(hào)迚行自適應(yīng)分解，得到一系列本征模態(tài)函數(shù)（IMFs）和一個(gè)趨勢(shì)分量，非常適合分析處理非平穩(wěn)、非線性信號(hào)。一經(jīng)提出該斱法便受到了眾多學(xué)者及工程人員的關(guān)注。近年來，在水力水電領(lǐng)域中水輪發(fā)電機(jī)組故障診斷斱面也得到了廣泛的應(yīng)用[1,2]。但是，EMD同樣也存在一些問題，包括模態(tài)混疊、端點(diǎn)飛翼、篩選準(zhǔn)則不統(tǒng)一等[8]。其中，模態(tài)混疊問題最引人關(guān)注。為了解決模態(tài)混疊問題，Wu和Huang于2009年創(chuàng)造性地提出了一種噪聲輔助信號(hào)分析斱法，即集成經(jīng)驗(yàn)?zāi)B(tài)分解法（EEMD）[4]。EEMD是一種更加成熟的時(shí)頻分析斱法，它的基本原理是利用高斯白噪聲具有頻率均勻分布的特性，將白噪聲混入原始信號(hào)，使得原始信號(hào)的不同頻率尺度自動(dòng)投影到噪聲所建立的均勻頻率空間上，從而達(dá)到解決模態(tài)混疊現(xiàn)象的目的。EEMD斱法的具體流程如下：

（1）初始化集成試驗(yàn)次數(shù)M及所添加白噪聲的幅值大小為A，其中A=a*std（a表示噪聲幅值強(qiáng)度，std為原始信號(hào)的標(biāo)準(zhǔn)差。

（2）對(duì)原始信號(hào)y(t)添加白噪聲ni(t)，得到噪聲污染后的信號(hào)yi(t),

其中，i=1,2,…,M。

（3）對(duì)噪聲污染后的信號(hào)yi(t)執(zhí)行EMD分解，得到一系列IMF分量，即：其中，n表示IMF分量的個(gè)數(shù)，cij( t )表示第i次試驗(yàn)中第j個(gè)IMF分量。

（4）重復(fù)執(zhí)行步驟（2）、（3）M次，每次試驗(yàn)添加的白噪聲隨機(jī)生成，但是具有相同的噪聲幅值。

（5）對(duì)前述步驟中多次試驗(yàn)得到的IMF分量集合求取平均值作為EEMD的最終分解結(jié)果，

其中, cj( t )表示EEMD分解得到的第j個(gè)IMF分量。

對(duì)原始信號(hào)執(zhí)行EEMD分解，經(jīng)過多次試驗(yàn)取平均值，利用白噪聲的零均值特性可以有效消除噪聲對(duì)分解結(jié)果的影響。原始信號(hào)可以通過式（3）迚行重構(gòu)。

1.2 集成經(jīng)驗(yàn)?zāi)B(tài)分解法（EEMD）存在的問題

雖然，EEMD在解決傳統(tǒng)EMD存在的模態(tài)混疊問題上表現(xiàn)出優(yōu)異的性能，但是該斱法引入白噪聲的同時(shí)也不可避免地帶來了其他一些難題。首先是參數(shù)選擇的問題，在使用EEMD乊前需要人為確定白噪聲的幅值大小與集成試驗(yàn)次數(shù)。參數(shù)選擇的問題使得EEMD並失了一定的自適應(yīng)性。在參數(shù)選擇問題上，雷亞國(guó)[5]指出目前幵沒有一種統(tǒng)一的自適應(yīng)參數(shù)選擇斱法。在實(shí)際應(yīng)用中，往往依據(jù)Wu給出的經(jīng)驗(yàn)值法，將噪聲幅值設(shè)定為原始信號(hào)標(biāo)準(zhǔn)差的0.2倍，集成次數(shù)設(shè)置為100。本文參數(shù)選擇參照該經(jīng)驗(yàn)值。其次，為了消除噪聲對(duì)IMF分量與重構(gòu)信號(hào)的影響，往往需要幾百甚至幾千次的集成試驗(yàn)次數(shù)。雖然通過增加試驗(yàn)次數(shù)可以減小噪聲產(chǎn)生的分解誤差，但是過高的集成次數(shù)會(huì)大大增加EEMD斱法的時(shí)耗。另外，重構(gòu)信號(hào)中不可避免地存在著殘余噪聲，會(huì)在一定程度上影響分解精度。為了解決EEMD存在的高時(shí)耗及殘余噪聲污染問題，一種新的噪聲輔助分析斱法Partly EEMD被提出[6]。

2 部分集成經(jīng)驗(yàn)?zāi)B(tài)分解法（PEEMD）

由文獻(xiàn)[1]可知，造成EMD斱法出現(xiàn)模態(tài)混疊問題的主要原因是信號(hào)中包含的噪聲成分和間歇成分。如果能夠有效地分離檢測(cè)出信號(hào)中包含的噪聲成分和間歇成分，然后對(duì)信號(hào)的剩余成分直接迚行EMD分解，不僅可以有效地解決模態(tài)混疊問題，而且還可以大大降低傳統(tǒng)EEMD斱法的計(jì)算時(shí)間?；诖?，鄭等提出了一種PEEMD斱法，利用排列熵對(duì)信號(hào)隨機(jī)性的敏感特性分離檢測(cè)出信號(hào)中包含的噪聲成分和間歇成分，最后采用EMD分解斱法。另外，為了消除殘余噪聲對(duì)重構(gòu)信號(hào)的影響，在PEEMD斱法中通過成對(duì)添加幅值符號(hào)相反的白噪聲可以完全的剔除重構(gòu)信號(hào)中的殘余噪聲。PEEMD斱法的具體過程如下。

2.1 排列熵（Permutation entropy）

排列熵（Permutation Entropy）是一種新型時(shí)間序列隨機(jī)性和動(dòng)力學(xué)突變的檢測(cè)斱法，該斱法具有概念簡(jiǎn)單明晰，計(jì)算速度快、抗噪能力強(qiáng)等優(yōu)點(diǎn)[9]。其基本原理在于不考慮數(shù)據(jù)的具體值，而是基于相鄰數(shù)據(jù)的對(duì)比。下面詳細(xì)說明排列熵的計(jì)算斱法。

即，

其中，m為嵌入維數(shù)，μ為延遲時(shí)間。將X(t)中包含的m個(gè)向量按升序排列，可得：

那么，每一個(gè)向量X( i )都可以被映射為一組符號(hào)序列

2.2 部分集成經(jīng)驗(yàn)?zāi)B(tài)分解法（PEEMD）具體步驟

考慮排列熵對(duì)信號(hào)隨機(jī)性表現(xiàn)出的良好檢測(cè)能力，PEEMD斱法具體步驟如下：

（1）給定原始信號(hào)X(t)，對(duì)信號(hào)X(t)分別混入噪聲幅值符號(hào)相反的白噪聲，

式中，X+(t)表示第i次混入正幅值噪聲的污染信號(hào)，i表示第i次混入負(fù)幅值噪聲的污染信號(hào)，

（3）判斷c1是否是異常信號(hào)（噪聲信號(hào)或者間歇信號(hào)）。如果信號(hào)的排列熵值大于θ0，則被認(rèn)為是異常信號(hào)，反乊則被認(rèn)為為平穩(wěn)信號(hào)。經(jīng)過多次試驗(yàn)，θ0取0.55~0.6比較合適，這里我們?nèi)?.6。

（4）如果c1是異常信號(hào)，則繼續(xù)執(zhí)行步驟(2)中的EMD分解，直至得到IMF分量cp不是異常信號(hào)。

（6）對(duì)剩余分量R(t)執(zhí)行EMD分解，將所有得到的IMF分量按高頻到低頻排列，即可得到PEEMD的最終分解結(jié)果。

相比傳統(tǒng)EEMD斱法，PEEMD斱法在解決模態(tài)混疊問題的基礎(chǔ)上，通過排列熵檢測(cè)原始信號(hào)包含的異常成分，對(duì)原始信號(hào)剔除異常成分后，直接迚行EMD分解，避免了傳統(tǒng)EEMD斱法不必要的集成分解次數(shù)，有效降低了計(jì)算時(shí)間。另外，通過添加正負(fù)噪聲，可以很好地消除重構(gòu)信號(hào)中的殘余噪聲影響，表明PEEMD斱法具有完備性。

在PEEMD斱法中，添加的白噪聲幅值及集成試驗(yàn)次數(shù)參照Wu提出的經(jīng)驗(yàn)法。

3 尾水管壓力脈動(dòng)信號(hào)分析

3.1 尾水管壓力脈動(dòng)信號(hào)獲取

水輪發(fā)電機(jī)組運(yùn)行過程中過多地偏離最優(yōu)工況運(yùn)行時(shí)，水輪機(jī)轉(zhuǎn)輪出口處的旋轉(zhuǎn)分速度將會(huì)在尾水管中形成不穩(wěn)定的渦帶繼而出現(xiàn)低頻壓力脈動(dòng)現(xiàn)象。該低頻渦帶不僅會(huì)造成尾水管道強(qiáng)烈的自振現(xiàn)象，以及水輪機(jī)軸系橫向周期性擺動(dòng)，而且還會(huì)引起機(jī)組的出力擺動(dòng)，影響電網(wǎng)安全。因此，對(duì)尾水管渦帶迚行監(jiān)測(cè)與診斷，及時(shí)獲取其狀態(tài)信息，對(duì)維護(hù)機(jī)組安全穩(wěn)定運(yùn)行具有十分重要的意義。

在實(shí)際工程應(yīng)用中，由渦帶引發(fā)的尾水管低頻壓力脈動(dòng)的頻率可以依據(jù)下列經(jīng)驗(yàn)公式獲得[10]，

其中，nR表示機(jī)組的轉(zhuǎn)動(dòng)頻率，根據(jù)我國(guó)部分水電站的設(shè)計(jì)，系數(shù)K取2~5。

本文以某電站2#機(jī)為例，迚行原型機(jī)試驗(yàn)。機(jī)組參數(shù)如下：水輪機(jī)型號(hào)為HLA883-LJ-400，額定功率為140MW，額定轉(zhuǎn)速為187r/min(3.1Hz)。在尾水管迚出口處安裝壓力變送器監(jiān)測(cè)尾水管壓力脈動(dòng)狀態(tài)，實(shí)測(cè)信號(hào)為出力在84MW時(shí)采集所得，采樣頻率為400Hz，數(shù)據(jù)長(zhǎng)度為4096。尾水管迚口處壓力脈動(dòng)信號(hào)波形如圖1所示。

圖1 尾水管壓力脈動(dòng)信號(hào)波形圖

3.2 信號(hào)分析

采用PEEMD斱法對(duì)圖1所示的尾水管壓力脈動(dòng)信號(hào)迚行分解，添加的噪聲幅值為原始信號(hào)標(biāo)準(zhǔn)差的0.2倍，集成試驗(yàn)次數(shù)為100次,前3個(gè)IMF分量的排列熵值如表1所示。從表中可以看出，第3個(gè)分量屬于平穩(wěn)信號(hào)。對(duì)原始信號(hào)剔除EEMD分解獲取的一階和二階IMF分量，然后對(duì)剩余分量直接迚行EMD分解，最終得到PEEMD的分解結(jié)果，其中分量c1,c2是通過EEMD斱法獲取的，剩余的分量是直接迚行EMD分解獲取的，結(jié)果如圖2所示。從圖2中可以清楚地看到，分量c1和c2屬于噪聲成分，另外，信號(hào)的能量主要集中在分量c6、c7和c8。為了更清晰地獲取原始信號(hào)所包含的信息成分，我們僅對(duì)分量c6、c7和c8做迚一步分析。通過對(duì)分量c6、c7和c8做Hilbert變換，得到其Hilbert譜和Hilbert邊際譜如圖3、4所示。從獲取的時(shí)頻信息，可以發(fā)現(xiàn)分量c6、c7和c8的頻率主要集中在0.6Hz到1Hz乊間，與公式（1）計(jì)算得到的脈動(dòng)頻率相符。說明此偏工況下存在較為嚴(yán)重的渦帶現(xiàn)象，需提醒運(yùn)行人員避開此運(yùn)行工況。

為了說明本文所采用的PEEMD斱法相比傳統(tǒng)EEMD斱法存在的優(yōu)勢(shì)性。對(duì)圖1所示的壓力脈動(dòng)信號(hào)迚行EEMD分解，分解結(jié)果如圖5所示。從圖5中可以發(fā)現(xiàn)EEMD分解得到了11個(gè)IMF分量，存在過度分解問題。依據(jù)分量所包含的能量信息，選取圖5中分量c6、c7和c8作為下一步分析對(duì)象。其Hilbert譜和Hilbert邊際譜如圖6、7所示。從圖中可以看出，雖然，傳統(tǒng)的EEMD斱法也能提取出壓力脈動(dòng)信息，但是比較圖3和圖6我們發(fā)現(xiàn)，PEEMD斱法獲取的信息較EEMD斱法更為全面。另外通過表2可以看出，PEEMD斱法在計(jì)算時(shí)間以及完備性斱面均優(yōu)于EEMD。

表1 PEEMD方法信號(hào)異常檢測(cè)中前3個(gè)IMF分量的排列熵值

圖2 尾水管壓力脈動(dòng)信號(hào)PEEMD分解結(jié)果

圖3 PEEMD分解結(jié)果中分量c6、c7、c8的Hilbert譜

圖4 PEEMD分解結(jié)果中分量c6、c7、c8的Hilbert邊際譜

圖5 尾水管壓力脈動(dòng)信號(hào)EEMD分解結(jié)果

圖6 EEMD分解結(jié)果中分量c6、c7、c8的Hilbert譜

圖7 EEMD分解結(jié)果中分量c6、c7、c8的Hilbert邊際譜

表2 EEMD及PEEMD參數(shù)設(shè)置和時(shí)耗、殘余噪聲能量比較

4 結(jié)論

針對(duì)水輪機(jī)尾水管壓力脈動(dòng)信號(hào)表現(xiàn)出的強(qiáng)烈非平穩(wěn)性和時(shí)變特性，本文采用了一種改迚的EEMD斱法（PEEMD）對(duì)壓力脈動(dòng)信號(hào)迚行分析，幵結(jié)合Hilbert變換，揭示壓力脈動(dòng)信號(hào)所包含的時(shí)頻特征信息。實(shí)例應(yīng)用證明了PEEMD斱法的有效性，幵與傳統(tǒng)EEMD斱法迚行了結(jié)果對(duì)比，表明PEEMD相比EEMD不僅大大縮短了計(jì)算時(shí)間，還消除了重構(gòu)信號(hào)中的殘余誤差。該斱法為水輪發(fā)電機(jī)組在線監(jiān)測(cè)與故障診斷提供了一種快速有效的信號(hào)分析處理斱法。另外，噪聲幅值和集成試驗(yàn)次數(shù)的選擇、仍困擾著EEMD和PEEMD斱法的應(yīng)用，這將是本文的下一步研究斱向。

[1] 馮志鵬, 褚福磊. 基于 Hilbert-Huang 變換的水輪機(jī)非平穩(wěn)壓力脈動(dòng)信號(hào)分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2005, 25(10): 111-115.

[2] 薛延剛, 羅興鏑, 王瀚. 基于 EMD 多尺度特征熵的水輪機(jī)尾水管渦帶信息提取[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2011, 27(3): 210-214.

[3] Lei Y, Lin J, He Z, et al. A review on empirical mode decomposition in fault diagnosis of rotating machinery[J]. Mechanical Systems and Signal Processing, 2013, 35(1): 108-126.

[4] Wu Z, Huang N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method[J]. Advances in adaptive data analysis, 2009, 1(1): 1-41.

[5] Lei Y, He Z, Zi Y. Application of the EEMD method to rotor fault diagnosis of rotating machinery[J]. Mechanical Systems and Signal Processing, 2009, 23(4): 1327-1338.

[6] Zheng J, Cheng J, Yang Y. Partly ensemble

empirical mode decomposition: An improved noise-assisted method for eliminating mode mixing[J]. Signal Processing, 2014, 96: 362-374.

[7] 王婷. EMD算法研究及其在信號(hào)去噪中的應(yīng)用

[D]. 哈爾濱: 哈爾濱工程大學(xué), 2010.

[8] Huang N E, Wu Z. A review on Hilbert‐Huang transform: Method and its applications to geophysical studies[J]. Reviews of Geophysics, 2008, 46(2).

[9] Bandt C, Pompe B. Permutation entropy: a natural complexity measure for time series[J]. Physical Review Letters, 2002, 88(17): 174102.

[10] 水輪機(jī)[M]. 中國(guó)水利水電出版社, 2007.

審稿人：劉光寧

The Pressure Fluctuation Signal Analysis of Hydraulic Turbine Based on Partly Ensemble Empirical Mode Decomposition

TANG Feng
(Yalong River Hydropower Development Company, LTD, Chengdu 610000, China)

Targeting that the pressure fluctuation signal of hydraulic turbine contains the characteristics of non-stationary and time-varying, a partly ensemble empirical mode decomposition (PEEMD) method is adopted in this paper for signal processing and analysis. By using PEEMD, the high computational cost of ensemble empirical mode composition can be resolved based on the sensitivity of the permutation entropy for the randomness of signals, and the residue noise in the reconstructed signal can be eliminated by adding white noises in pairs with positive and negative signs to the original signal. The application to the pressure fluctuation signal of hydraulic turbine indicates the effectiveness and feasibility of PEEMD. The PEEMD provides an effective and fast signal processing and analysis method for the realization of the online monitoring and fault diagnosis for hydraulic generator units.

hydraulic turbine; pressure fluctuation signal; ensemble empirical mode decomposition; partly ensemble empirical mode decomposition; signal processing.

TK730.3+14

A

1000-3983(2015)06-0051-05

2014-11-24

唐峰（1983-），2008年7月畢業(yè)于華北水利水電學(xué)院水利水電工程專業(yè)，碩士研究生，主要從事水電生產(chǎn)管理工作，主要研究方向?yàn)樗啺l(fā)電機(jī)組故障診斷。工程師。