趙 悅，方彥軍，董政呈

?

基于狀態(tài)識(shí)別的經(jīng)驗(yàn)?zāi)B(tài)分解法火電廠運(yùn)行數(shù)據(jù)預(yù)處理

趙 悅，方彥軍，董政呈

（武漢大學(xué)動(dòng)力與機(jī)械學(xué)院，湖北 武漢 430072）

針對(duì)火電機(jī)組運(yùn)行監(jiān)測(cè)數(shù)據(jù)量大且復(fù)雜情況下的去噪問題，提出一種基于狀態(tài)識(shí)別的改進(jìn)經(jīng)驗(yàn)?zāi)B(tài)分解去噪（SREMD）算法。該算法以經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）方法為基礎(chǔ)，首先運(yùn)用基于拉伊達(dá)（Pauta）準(zhǔn)則的濾波方法去除顯著異常數(shù)據(jù)，然后根據(jù)數(shù)據(jù)連續(xù)變化率動(dòng)態(tài)識(shí)別機(jī)組運(yùn)行狀態(tài)，最后根據(jù)機(jī)組運(yùn)行狀態(tài)的穩(wěn)態(tài)過程和過渡過程分別進(jìn)行針對(duì)性EMD去噪，以適應(yīng)火電機(jī)組的狀態(tài)切換特性。將該算法用于實(shí)際機(jī)組運(yùn)行數(shù)據(jù)，結(jié)果表明，本文算法有效完成了火電機(jī)組監(jiān)測(cè)數(shù)據(jù)的去噪預(yù)處理，在保持信號(hào)整體趨勢(shì)的基礎(chǔ)上能達(dá)到更好的去噪效果。

火電機(jī)組；數(shù)據(jù)預(yù)處理；去噪；濾波；Pauta準(zhǔn)則；EMD；狀態(tài)識(shí)別

火電機(jī)組運(yùn)行數(shù)據(jù)是由現(xiàn)場(chǎng)成千上萬個(gè)傳感器獲得，這些傳感器實(shí)際工作環(huán)境復(fù)雜，很多傳感器在高溫、高壓等極端條件下工作，其精度會(huì)受到影響，出現(xiàn)故障的概率也較大。此時(shí)傳感器采集到的數(shù)據(jù)就會(huì)存在各種異常值，若對(duì)這些數(shù)據(jù)不做相應(yīng)的預(yù)處理，后續(xù)數(shù)據(jù)的使用就會(huì)受很大的影響。因此，采集數(shù)據(jù)需進(jìn)行異常值檢測(cè)、噪聲剔除等數(shù)據(jù)預(yù)處理操作。

火電機(jī)組運(yùn)行數(shù)據(jù)具有如下特點(diǎn)：1）測(cè)量點(diǎn)多，數(shù)據(jù)量龐大；2）測(cè)量環(huán)境復(fù)雜多樣，數(shù)據(jù)夾雜大量噪聲；3）數(shù)據(jù)類型多樣，具有較強(qiáng)的動(dòng)態(tài)性[1-3]。傳統(tǒng)數(shù)據(jù)預(yù)處理方法直接應(yīng)用于電廠運(yùn)行數(shù)據(jù)預(yù)處理時(shí)會(huì)出現(xiàn)很多問題，達(dá)不到預(yù)期效果。如常用的量程檢測(cè)法，只可識(shí)別出少量的異常數(shù)據(jù)；簡(jiǎn)單濾波算法只可剔除少量誤差較大的噪聲數(shù)據(jù)。而大部分的大數(shù)據(jù)預(yù)處理方法不能直接用于電廠大數(shù)據(jù)的預(yù)處理，當(dāng)前關(guān)于電廠運(yùn)行數(shù)據(jù)的大數(shù)據(jù)分析側(cè)重于參數(shù)研究、運(yùn)行優(yōu)化等方面[4-11]，對(duì)電廠數(shù)據(jù)的濾波去噪等針對(duì)性預(yù)處理的研究則較少。

一些新的數(shù)據(jù)預(yù)處理方法被提出：文獻(xiàn)[12]采用SSD（steady state detection）算法和EWMA（exponentially weighted moving average）濾波進(jìn)行穩(wěn)態(tài)檢測(cè)并結(jié)合改進(jìn)的方法實(shí)現(xiàn)了穩(wěn)態(tài)檢測(cè)數(shù)據(jù)的預(yù)處理；文獻(xiàn)[13]采用聚類算法實(shí)現(xiàn)了電力數(shù)據(jù)的異常、毛刺以及極大極小值的處理與優(yōu)化；文 獻(xiàn)[14]針對(duì)電廠測(cè)量數(shù)據(jù)一般出自冗余測(cè)量系統(tǒng)的特點(diǎn)，考慮不同變量的測(cè)量精度，根據(jù)基于最大似然原理的數(shù)據(jù)校正模型進(jìn)行顯著性測(cè)量誤差檢測(cè)，并通過校正得到合理數(shù)據(jù)。但這些方法大多只針對(duì)電廠靜態(tài)數(shù)據(jù)，對(duì)動(dòng)態(tài)性較強(qiáng)的數(shù)據(jù)處理效果不佳；且僅能修正少數(shù)顯著異常數(shù)據(jù)，數(shù)據(jù)中各種復(fù)雜因素引起的噪聲未能得到很好地處理。

異常值檢測(cè)是數(shù)據(jù)預(yù)處理的基礎(chǔ)，其中發(fā)展最成熟的是基于統(tǒng)計(jì)的異常值檢測(cè)技術(shù)[15-16]。Pauta準(zhǔn)則適用于數(shù)量大的數(shù)據(jù)異常值處理，該方法較簡(jiǎn)便?；诖说臑V波方法是以隨機(jī)噪聲假設(shè)為基礎(chǔ)，通過計(jì)算均方誤差設(shè)置閾值，有效剔除異常數(shù)據(jù)，并根據(jù)數(shù)據(jù)在時(shí)間上的連續(xù)性，對(duì)數(shù)據(jù)進(jìn)行修復(fù)，從而達(dá)到濾波效果。但該方法僅對(duì)明顯的異常數(shù)據(jù)進(jìn)行過濾，對(duì)剩余大部分噪聲過濾效果不佳。對(duì)具有非平穩(wěn)特性數(shù)據(jù)的降噪而言，一般的中值濾波、滑動(dòng)平均濾波等方法效果不明顯，取而代之的是經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）[17-20]、小波分解等方法。傳統(tǒng)的EMD方法存在模態(tài)混疊的問題，即一個(gè)單獨(dú)的IMF信號(hào)中會(huì)含有全異的尺度，且相同的尺度會(huì)存在于不同的IMF信號(hào)中，對(duì)于存在狀態(tài)切換的火電機(jī)組數(shù)據(jù)，傳統(tǒng)的EMD方法在去噪的同時(shí)會(huì)在一定程度上導(dǎo)致信號(hào)的損失，使信號(hào)只保留基本趨勢(shì)而損失很多可能很重要的細(xì)節(jié)信號(hào)，存在噪聲剔除不凈且不能很好跟蹤數(shù)據(jù)趨勢(shì)的缺點(diǎn)。另外，實(shí)驗(yàn)表明傳統(tǒng)EMD方法容易受信號(hào)中的異常值影響，將異常值引入固有模態(tài)函數(shù)，從而使EMD重構(gòu)后的信號(hào)明顯偏離原信號(hào)。

因此，本文以Pauta準(zhǔn)則和EMD方法為基礎(chǔ)，充分結(jié)合火電機(jī)組運(yùn)行監(jiān)測(cè)數(shù)據(jù)的特點(diǎn)提出一種改進(jìn)的狀態(tài)識(shí)別濾波算法，即SREMD算法。該算法在拉伊達(dá)準(zhǔn)則濾波的基礎(chǔ)上，通過識(shí)別機(jī)組的運(yùn)行狀態(tài)，將數(shù)據(jù)分為穩(wěn)態(tài)過程數(shù)據(jù)和過渡過程數(shù)據(jù)，并運(yùn)用EMD法進(jìn)行針對(duì)性處理。在充分結(jié)合電廠數(shù)據(jù)特點(diǎn)的同時(shí)發(fā)揮EMD法處理非平穩(wěn)數(shù)據(jù)的優(yōu)勢(shì)，對(duì)火電機(jī)組數(shù)據(jù)進(jìn)行去噪等處理。

1 改進(jìn)狀態(tài)識(shí)別濾波算法

1.1 基本思想

在進(jìn)行EMD分解前先對(duì)原始數(shù)據(jù)進(jìn)行異常值剔除來消除顯著異常值對(duì)EMD重構(gòu)的影響，即需先進(jìn)行Pauta準(zhǔn)則濾波。圖1為火電廠給水溫度與發(fā)電功率原始數(shù)據(jù)。由圖1可知，火電機(jī)組數(shù)據(jù)存在動(dòng)態(tài)性較強(qiáng)的特點(diǎn)，且火電機(jī)組為一個(gè)狀態(tài)切換系統(tǒng)，即機(jī)組運(yùn)行狀態(tài)可以根據(jù)整體變化率分為穩(wěn)態(tài)運(yùn)行過程和狀態(tài)切換過渡過程2種狀態(tài)。

圖1 給水溫度與發(fā)電機(jī)功率原始數(shù)據(jù)

因此，本文提出一種穩(wěn)定狀態(tài)識(shí)別算法，其基本思想是通過調(diào)節(jié)窗口寬度，根據(jù)數(shù)據(jù)的連續(xù)變化率識(shí)別出過渡過程和穩(wěn)態(tài)過程。對(duì)此計(jì)算如下：

式中：Δx為相鄰兩數(shù)據(jù)的變化率；Δmax為最大變化率；Δmax為最小變化率；為設(shè)定窗口寬度，其最優(yōu)值通過最小化復(fù)合評(píng)價(jià)指標(biāo)來確定；12為閾值參數(shù)。

長(zhǎng)度的窗口內(nèi)，相鄰兩數(shù)據(jù)變化率的均值若大于窗口內(nèi)最大變化率與閾值參數(shù)1的乘積，則可判斷此段為上升過渡過程；同樣相鄰兩數(shù)據(jù)變化率的均值小于窗口內(nèi)最小變化率與閾值參數(shù)2的乘積，則此時(shí)處于數(shù)值下降過渡過程。否則，判定為平穩(wěn)階段，即為穩(wěn)態(tài)過程。

狀態(tài)識(shí)別后，針對(duì)穩(wěn)態(tài)過程中信號(hào)整體趨于平穩(wěn)，但存在高頻噪聲的特點(diǎn)，進(jìn)行局部EMD分解，并進(jìn)一步設(shè)計(jì)重構(gòu)尺度對(duì)穩(wěn)態(tài)過程數(shù)據(jù)進(jìn)行重構(gòu)。同時(shí)對(duì)原數(shù)據(jù)直接進(jìn)行EMD濾波。之后，取直接EMD濾波的過渡過程的數(shù)據(jù)與重構(gòu)后的穩(wěn)態(tài)過程段數(shù)據(jù)進(jìn)行平滑對(duì)接，從而實(shí)現(xiàn)2種過程識(shí)別并進(jìn)行針對(duì)性處理。

1.2 算法流程

采用SREMD算法對(duì)火電機(jī)組運(yùn)行數(shù)據(jù)進(jìn)行濾波處理時(shí)，首先導(dǎo)入機(jī)組運(yùn)行數(shù)據(jù)，采用Pauta準(zhǔn)則濾波，清除顯著異常值數(shù)據(jù)；然后根據(jù)式(1)識(shí)別出穩(wěn)態(tài)數(shù)據(jù)和過渡過程數(shù)據(jù)；對(duì)穩(wěn)態(tài)數(shù)據(jù)使用EMD方法進(jìn)行針對(duì)性分解和重構(gòu)，清除噪聲數(shù)據(jù)；同時(shí)對(duì)過渡數(shù)據(jù)進(jìn)行直接EMD濾波；最后將穩(wěn)態(tài)數(shù)據(jù)與過渡過程數(shù)據(jù)重新組合并平滑化，完成數(shù)據(jù)的整體濾波。其算法流程如圖2所示。

圖2 SREMD算法流程

1.3 性能評(píng)價(jià)指標(biāo)

通過不同的性能指標(biāo)[21]來評(píng)價(jià)SREMD算法的去噪性能。常用的性能指標(biāo)包括：1）均方根誤差RMSE，可以體現(xiàn)信號(hào)的整體偏差信息；2）信噪比SNR，體現(xiàn)噪聲信息對(duì)整體信號(hào)的影響；3）平滑度，體現(xiàn)信號(hào)的局部變異信息，即是否有較多的局部突變。各性能指標(biāo)計(jì)算公式如下。

式中signal和noise分別表示信號(hào)和噪聲的能量。

一般而言，當(dāng)數(shù)據(jù)真值已知時(shí)，均方根誤差越小越好，信噪比越大越好，平滑度越小越好。由于火電機(jī)組在線運(yùn)行的真實(shí)數(shù)據(jù)未知，測(cè)量數(shù)據(jù)本身帶有大量噪聲，故單獨(dú)的指標(biāo)無法直接通過與真實(shí)信號(hào)的比較來反映去噪效果。

經(jīng)過分析發(fā)現(xiàn)，均方根誤差和信噪比都是反映信號(hào)細(xì)節(jié)信息的指標(biāo)，而平滑度則反映信號(hào)的整體趨勢(shì)特征，且均方根誤差與平滑度和信噪比之間呈負(fù)相關(guān)。因此，本文從數(shù)據(jù)整體的角度考慮，通過構(gòu)造復(fù)合性能評(píng)價(jià)指標(biāo)對(duì)算法進(jìn)行綜合性的定量評(píng)價(jià)，具體計(jì)算為

根據(jù)復(fù)合指標(biāo)可以計(jì)算SREMD算法的最優(yōu)窗口寬度為

為了更加合理地確定復(fù)合指標(biāo)中各指標(biāo)的權(quán)值，本文通過指標(biāo)的變異系數(shù)定權(quán)方法[22]來確定各指標(biāo)權(quán)值。變異系數(shù)也稱為標(biāo)準(zhǔn)差率，是標(biāo)準(zhǔn)差與均值的比值，能反映指標(biāo)數(shù)值的變異程度，變異系數(shù)越大的指標(biāo)越難實(shí)現(xiàn)，因此權(quán)值越大。

式中，CVRMSE表示均方誤差的變異系數(shù)，RMSE、RMSE分別表示指標(biāo)的標(biāo)準(zhǔn)差和均值。

則

2 實(shí)例分析

以某火電機(jī)組發(fā)電機(jī)功率24 h的實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)進(jìn)行實(shí)驗(yàn)仿真，并從直觀效果和量化性能指標(biāo)2個(gè)方面對(duì)算法效果進(jìn)行比較分析。對(duì)Pauta準(zhǔn)則濾波、EMD法與本文SREMD算法對(duì)數(shù)據(jù)的預(yù)處理效果進(jìn)行驗(yàn)證分析。

2.1 3種算法去噪聲效果定性比較

對(duì)于火電廠大數(shù)據(jù)而言，Pauta準(zhǔn)則濾波能較好地剔除的異常值，EMD法對(duì)高頻數(shù)據(jù)有較好地消除噪聲的效果，但在處理火電機(jī)組數(shù)據(jù)時(shí)存在明顯的局限性，因此本文提出的改進(jìn)的火電機(jī)組數(shù)據(jù)SREMD算法是很有必要的。

圖3為基于Pauta準(zhǔn)則的數(shù)據(jù)濾波。由圖3可見，原始數(shù)據(jù)中存在明顯的毛刺，如圖3中的3處圓圈標(biāo)記處，這些數(shù)據(jù)不滿足Pauta準(zhǔn)則的條件，能夠予以剔除。而對(duì)于幅值較小的其他隨機(jī)噪聲，濾波效果欠佳。

圖3 基于Pauta準(zhǔn)則的數(shù)據(jù)濾波

圖4為不同重構(gòu)尺度的EMD濾波。由圖4可見，不同重構(gòu)尺度下EMD濾波效果不同。重構(gòu)系數(shù)過小時(shí)濾波效果不佳，過大又會(huì)導(dǎo)致重構(gòu)后的數(shù)據(jù)波形在很多細(xì)節(jié)處嚴(yán)重失真。

對(duì)比圖4中不同重構(gòu)系數(shù)發(fā)現(xiàn)：=4時(shí)噪聲去除效果不佳；=6時(shí)會(huì)導(dǎo)致原始數(shù)據(jù)信息損失，不能較好地反映原始數(shù)據(jù)信號(hào)的趨勢(shì)；5時(shí)濾波效果最佳。因此，選取=5時(shí)的EMD濾波。此時(shí)對(duì)火電機(jī)組發(fā)電機(jī)功率數(shù)據(jù)的噪聲濾除效果較好，但依然不能很好地跟蹤數(shù)據(jù)的變化趨勢(shì)。所以本文結(jié)合2種方法，在應(yīng)用Pauta準(zhǔn)則去除顯著異常值后再在狀態(tài)識(shí)別的基礎(chǔ)上針對(duì)性應(yīng)用EMD濾波。

圖4 不同重構(gòu)尺度的EMD濾波

根據(jù)式(2)—式(5)，計(jì)算不同窗口寬度下SREMD算法去噪聲復(fù)合性能指標(biāo)，通過尋找復(fù)合指標(biāo)的最小值來確定最優(yōu)的窗口寬度。由表1可知，當(dāng)SREMD算法的窗口寬度=19時(shí)，復(fù)合指標(biāo)值最小，即算法的綜合去噪聲性能達(dá)到最優(yōu)。因此，選取窗口寬度=19作為SREMD算法的參數(shù)。

表1 不同窗口寬度的復(fù)合指標(biāo)

Tab.1 The composite indexes with different window widths

圖5為發(fā)電機(jī)功率的相對(duì)穩(wěn)態(tài)過程識(shí)別。圖6為基于狀態(tài)識(shí)別的SREMD濾波。SREMD算法是先對(duì)原始數(shù)據(jù)進(jìn)行Pauta濾波，剔除明顯噪聲，然后結(jié)合整體EMD濾波與穩(wěn)態(tài)過程局部EMD濾波，既很好地剔除了高頻噪聲，又很好地反映數(shù)據(jù)整體趨勢(shì)，效果突出，特別是對(duì)于穩(wěn)態(tài)區(qū)域，在保證數(shù)據(jù)整體趨勢(shì)的情況下，對(duì)于高頻噪聲的處理效果顯著。且與傳統(tǒng)EMD法相比而言，SREMD法在分解重構(gòu)前進(jìn)行了Pauta濾波，消除了異常值對(duì)分解重構(gòu)的影響，避免了直接使用EMD法時(shí)出現(xiàn)的固有模態(tài)函數(shù)受異常值影響從而導(dǎo)致數(shù)據(jù)重構(gòu)后偏離原始信號(hào)的問題。

圖5 穩(wěn)態(tài)過程識(shí)別與濾波

圖6 基于狀態(tài)識(shí)別的SREMD濾波

2.2 性能指標(biāo)定量分析

表2為3種算法的歸一化均方根誤差、信噪比、平滑度與復(fù)合性能指標(biāo)對(duì)比結(jié)果。根據(jù)表2的數(shù)據(jù)，就單項(xiàng)指標(biāo)而言，Pauta濾波的均方誤差最小而信噪比最高，這是因?yàn)閮H去除了顯著異常的噪聲，與原始數(shù)據(jù)差異很小。EMD濾波的平滑度指標(biāo)最好，是因?yàn)槠渲苯尤コ烁哳l部分的所有本征母函數(shù)，剩余部分就是頻率較低的成分，但是這樣處理會(huì)導(dǎo)致很多信息丟失。SREMD濾波的3項(xiàng)單性能指標(biāo)均較好，其中基于狀態(tài)識(shí)別的SREMD濾波算法復(fù)合指標(biāo)最佳，即其綜合去噪聲能力最強(qiáng)。

表2 3種算法性能指標(biāo)比較

Tab.2 The performance index comparison between three algorithms

3 結(jié) 論

本文將EMD法與Pauta準(zhǔn)則相結(jié)合，針對(duì)火電機(jī)組運(yùn)行監(jiān)測(cè)數(shù)據(jù)的特性，提出了一種基于狀態(tài)識(shí)別的SREMD濾波算法。該算法首先采用Pauta準(zhǔn)則修復(fù)數(shù)據(jù)中的顯著異常值，去除異常數(shù)據(jù)對(duì)本征模函數(shù)的影響，通過對(duì)數(shù)據(jù)進(jìn)行機(jī)組運(yùn)行狀態(tài)識(shí)別，并運(yùn)用EMD方法對(duì)平穩(wěn)狀態(tài)數(shù)據(jù)進(jìn)行針對(duì)性分解重構(gòu)去除高頻噪聲，將平穩(wěn)狀態(tài)濾波數(shù)據(jù)與直接EMD濾波得到的過渡過程數(shù)據(jù)平滑對(duì)接，從而在保證信號(hào)趨勢(shì)的前提下較好地去除噪聲。

本文最后設(shè)計(jì)了基于變異指數(shù)的加權(quán)復(fù)合評(píng)價(jià)指標(biāo)，根據(jù)指標(biāo)確定SREMD算法的最優(yōu)窗口寬度，并采用實(shí)際數(shù)據(jù)進(jìn)行驗(yàn)證SREMD算法的有效性。結(jié)果表明，該算法能夠更有效對(duì)帶有狀態(tài)切換的火電機(jī)組數(shù)據(jù)進(jìn)行去噪，并能夠很好的保持機(jī)組數(shù)據(jù)整體趨勢(shì)。

［1］ 楊世銘, 陶文銓. 傳熱學(xué)[J]. 4版. 北京: 高等教育出版社, 2006: 37-45.

YANG Shiming, TAO Wenquan. Heat transfer[M]. 4th ed. Beijing: Higher Education Press, 2006: 37-45.

［2］ 陳雪. 熱工測(cè)量數(shù)據(jù)預(yù)處理方法研究[D]. 北京: 華北電力大學(xué), 2014: 6-9.

CHEN Xue. Research on data preprocessing method for thermal parameters[D]. Beijing: North China Electric Power University, 2014: 6-9.

［3］ FAMILI A, SHEN W M, RICHARD W, et al. Data preprocessing and intelligent data analysis[J]. Intelligent Data Analysis, 1997, 1(1): 3-23.

［4］ WANG C L, ZHANG C L, ZHANG P T. Denoising algorithm based on wavelet adaptive threshold[J]. Physics Procedia, 2012, 24(1): 678-685.

［5］ 朱龍霞, 肖明美, 栗濤, 等. 基于時(shí)間序列分析的火電機(jī)組運(yùn)行優(yōu)化研究[J]. 河北工業(yè)科技, 2017(2): 127-131.

ZHU Longxia, XIAO Mingmei, LI Tao, et al. Research of operation optimization of thermal power unit based on time series analysis[J]. Hebei Industrial Technology, 2017(2): 127-131.

［6］ 劉麗平, 葉春, 忻建華. 電廠在線性能分析及故障診斷系統(tǒng)[J]. 熱力發(fā)電, 2003, 32(11): 76-77.

LIU Liping, YE Chun, XIN Jianhua. Power plant on-line performance analysis and fault diagnosis system[J]. Thermal Power Generation, 2003, 32(11): 76-77.

［7］ 馬海林, 孫秋紅, 徐欣航. 基于數(shù)據(jù)挖掘技術(shù)的電廠運(yùn)行參數(shù)目標(biāo)值優(yōu)化[J]. 河北工業(yè)科技, 2012(5): 303-307.

MA Hailin, SUN Qiuhong, XU Xinhang. Study on optimization of power plant running parameters based on data minning[J]. Hebei Industrial Technology, 2012(5): 303-307.

［8］ 李輝, 陳教超, 司風(fēng)琪, 等. 數(shù)據(jù)預(yù)處理技術(shù)在電廠運(yùn)行優(yōu)化系統(tǒng)中的應(yīng)用研究[J].華東電力, 2007(11): 110-113.

LI Hui, CHEN Jiaochao, SI Fengqi, et al. Application of data preprocessing techonlogy to power plant perfor- mance optimization systems[J]. East China Electric Power, 2007(11): 110-113.

［9］ 王一蓉, 王瑞杰, 陳文剛, 等. 基于遺傳優(yōu)化的調(diào)控系統(tǒng)缺失數(shù)據(jù)填補(bǔ)算法[J].電力系統(tǒng)保護(hù)與控制, 2016, 44(21): 182-186.

WANG Yirong, WANG Ruijie, CHEN Wengang, et al. A missing data filling algorithm for dispatching and control system based on genetic optimization[J]. Power System Protection and Control, 2016, 44(21): 182-186.

［10］ 史林軍, 周佳佳, 溫榮超, 等. 基于經(jīng)驗(yàn)?zāi)B(tài)分解與小波分析相結(jié)合的風(fēng)電功率平滑控制[J]. 電力系統(tǒng)保護(hù)與控制, 2016, 44(24): 9-16.

SHI Linjun, ZHOU Jiajia, WEN Rongchao, et al. Power smoothing control of wind power based on combination of empirical mode decomposition and wavelet analysis[J]. Power System Protection and Control, 2016, 44(24): 9-16.

［11］ 楊大勇, 葛琪, 董永超, 等. 基于K均值聚類的光伏電站運(yùn)行狀態(tài)模式識(shí)別研究[J]. 電力系統(tǒng)保護(hù)與控制, 2016, 44(14): 25-30.

YANG Dayong, GE Qi, DONG Yongchao, et al. Research on operation state pattern recognition of PV station based on the principle of K-means clustering[J]. Power System Protection and Control, 2016, 44(14): 25-30.

［12］ 楊雯宇. 基于穩(wěn)態(tài)檢測(cè)的電廠數(shù)據(jù)預(yù)處理研究[D]. 北京: 華北電力大學(xué), 2017: 6-26.

YANG Wenyu. Analysis on power plant data preprocessing based on steady state detection[D]. Beijing: North China Electric Power University, 2017: 6-26.

［13］ 蘇舟, 李燦, 姚李孝, 等. 電力負(fù)荷數(shù)據(jù)預(yù)處理研究及應(yīng)用[J]. 電網(wǎng)與清潔能源, 2017, 33(5): 40-43.

SU Zhou, LI Can, YAO Lixiao, et al. Research and application of pretreatment of electrical load data[J]. Power System and Clean Energy, 2017, 33(5): 40-43.

［14］ 劉福國, 王學(xué)同, 蘇相河, 等. 基于系統(tǒng)測(cè)量冗余的電廠異常運(yùn)行數(shù)據(jù)檢測(cè)與校正[J]. 中國電機(jī)工程學(xué)報(bào), 2003(7): 204-207.

LIU Fuguo, WANG Xuetong, SU Xianghe, et al. Detection and correction on the abnormal operation data based on measurement in a power plant[J]. Proceedings of the CSEE, 2003(7): 204-207.

［15］ GUMBELL E J. Statistics of extremes[J]. Applied Ocean Research, 1958, 43(1): 1-30.

［16］ ZHAO J, ZHANG G, LA SCALA M, et al. Multistage phasor-aided bad data detection and identification[C]. IEEE Power & Energy Society General Meeting, 2015 (GHTCE).

［17］ HUANG N E.Review of empirical mode decomposition analysis[J]. Proceedings of SPIE, 2001(4391): 71-79.

［18］ REHMAN N, MANDIC D P. Multivariate empirical mode decomposition[J]. Proceedings Mathematical Physical & Engineering Sciences, 2010, 466(2117): 1291-1302.

［19］ 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[J]. Proceedings Mathematical Physical & Engineering Sciences, 1998, 454(1971): 903-995.

［20］ 劉霞, 黃陽, 黃敬, 等. 基于經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)的小波熵閾值地震信號(hào)去噪[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2016, 46(1): 262-269.

LIU Xia, HUANG Yang, HUANG Jing, et al. Wavelet entropy threshold seismic signal denoising based on empirical mode decomposition (EMD)[J]. Journal of Jilin University (Earth Science Edition), 2016, 46(1): 262-269.

［21］ 陶珂, 朱建軍. 多指標(biāo)融合的小波去噪最佳分解尺度選擇方法[J]. 測(cè)繪學(xué)報(bào), 2012, 41(5): 749-755.

TAO Ke, ZHU Jianjun. A hybrid indicator for deter- mining the best decomposition scale of wavelet denoising[J]. Chinese Journal of Chemical Engineering, 2012, 41(5): 749-755.

［22］ VOLKER D. Rhythm and speech rate: a variation coefficient for Δe[J]. Language and Language-Processing, 2006(1): 231-241.

Operating data preprocessing using EMD method with state recognition for thermal power plants

ZHAO Yue, FANG Yanjun, DONG Zhengcheng

(School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China)

The data from the operating thermal power units are very large and contains complex noise, which is hard to eliminate the noise thoroughly. To solve this problem, an improved empirical mode decomposition filtering algorithm based on the state recognition (SREMD) was proposed. On the basis of empirical mode decomposition (EMD), this SREMD algorithm firstly eliminates the significantly abnormal data by using a filtering method which is based on the Pauta criterion. Then, it recognizes the running states of the unit according to continuous change rate of the data. Finally, it performs corresponding EMD denoising according to the steady state and transition process of unit operation, to adapt to the state switching characteristics. Moreover, this algorithm was applied for actual unit operation data. The results show that, the proposed SREMD algorithm effectively completed the pre-processing of the monitoring data of thermal power units and achieved a better denoising effect on the basis of maintaining the overall trend of the signal.

thermal power unit, data preprocessing, denoising, filtering, Pauta criterion, EMD, state recognition

National Natural Science Foundation of China (51707153); China Postdoctoral Science Foundation (2017M612499)

趙悅(1994—)，女，碩士研究生，主要研究方向?yàn)橹悄芸刂评碚撗芯考皯?yīng)用，2017202080058@whu.edu.cn。

TM932

B

10.19666/j.rlfd.201803045

趙悅, 方彥軍, 董政呈. 基于狀態(tài)識(shí)別的經(jīng)驗(yàn)?zāi)B(tài)分解法火電廠運(yùn)行數(shù)據(jù)預(yù)處理[J]. 熱力發(fā)電, 2019, 48(1): 49-54. ZHAO Yue, FANG Yanjun, DONG Zhengcheng. Operating data preprocessing using EMD method with state recognition for thermal power plants[J]. Thermal Power Generation, 2019, 48(1): 49-54.

2018-03-16

國家自然科學(xué)基金項(xiàng)目(51707153)；中國博士后科學(xué)基金(2017M612499)

方彥軍(1957—)，男，博士，教授，主要研究方向?yàn)榛痣姀S熱工控制技術(shù)，yjfang@whu.edu.cn。

（責(zé)任編輯 杜亞勤）