董興輝， 張 光， 程友星， 王 帥

(1．華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，北京102206；2．河南理工大學(xué) 電氣工程與自動(dòng)化學(xué)院，河南焦作 454000)

風(fēng)電機(jī)組運(yùn)行環(huán)境惡劣，受氣象等多種不確定因素的影響，容易出現(xiàn)性能與狀態(tài)劣化，頻繁出現(xiàn)故障。風(fēng)電機(jī)組的關(guān)鍵部件一旦失效，檢修時(shí)間較長(zhǎng)，不僅影響發(fā)電量，也增加了風(fēng)電場(chǎng)的運(yùn)維成本[1-3]。風(fēng)機(jī)軸承作為風(fēng)電機(jī)組關(guān)鍵部件，在機(jī)組機(jī)械故障中占很高比例，發(fā)電機(jī)、齒輪箱的機(jī)械故障中約80%是由軸承失效而導(dǎo)致的[4]。較早掌握風(fēng)電機(jī)組軸承的劣化程度和劣化趨勢(shì)，實(shí)時(shí)了解風(fēng)電機(jī)組狀態(tài)，在可預(yù)知情況下合理安排檢修，有助于提高整機(jī)的可靠性和利用率。

近年來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)風(fēng)電機(jī)組軸承健康研究集中在劣化程度與狀態(tài)評(píng)估2方面。安學(xué)利等[5-7]基于振動(dòng)參數(shù)量計(jì)算風(fēng)電機(jī)組主軸承、齒輪箱軸承和發(fā)電機(jī)軸承等的劣化程度；李輝等[1,8]基于溫度參數(shù)量研究風(fēng)電軸承劣化程度，并在此基礎(chǔ)上進(jìn)一步對(duì)軸承進(jìn)行狀態(tài)評(píng)估，但沒(méi)有對(duì)機(jī)組軸承劣化發(fā)展趨勢(shì)開(kāi)展預(yù)測(cè)研究。

筆者基于數(shù)據(jù)采集與監(jiān)視控制(SCADA)系統(tǒng)監(jiān)測(cè)數(shù)據(jù)，應(yīng)用改進(jìn)的集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)和時(shí)間序列神經(jīng)網(wǎng)絡(luò)，通過(guò)建立預(yù)測(cè)模型，研究風(fēng)電機(jī)組軸承健康狀態(tài)溫度劣化趨勢(shì)。首先，利用SCADA監(jiān)測(cè)的有關(guān)參數(shù)(如風(fēng)速、功率、發(fā)電機(jī)轉(zhuǎn)速和環(huán)境溫度等)與軸承溫度的相關(guān)系數(shù)，歸一化得到影響權(quán)重，基于溫度特征量構(gòu)建軸承健康劣化度模型；然后利用改進(jìn)后的EEMD對(duì)劣化趨勢(shì)進(jìn)行分解，得到一系列相對(duì)平穩(wěn)的本征模態(tài)函數(shù)(IMF)分量和一個(gè)剩余分量；最后，應(yīng)用時(shí)間序列神經(jīng)網(wǎng)絡(luò)對(duì)各類分量建立預(yù)測(cè)模型，疊加所有預(yù)測(cè)分量，得到最終預(yù)測(cè)結(jié)果。

1 軸承健康劣化度模型

1.1 軸承劣化度

選用軸承運(yùn)行溫度參數(shù)作為風(fēng)電機(jī)組健康裝填評(píng)價(jià)指標(biāo)，采用相對(duì)劣化度d′(t)來(lái)表征風(fēng)電機(jī)組軸承當(dāng)前的相對(duì)劣化程度，表達(dá)式如下：

(1)

1.2 軸承健康狀態(tài)模型

軸承健康狀態(tài)模型的建立如下。

(1) 樣本數(shù)據(jù)。

SCADA數(shù)據(jù)中，與風(fēng)電機(jī)組軸承溫度相關(guān)的主要參數(shù)有風(fēng)速、功率、發(fā)電機(jī)轉(zhuǎn)速和環(huán)境溫度[8]。從SCADA系統(tǒng)中選取n組健康的、同一時(shí)刻的軸承溫度、風(fēng)速、有功功率、發(fā)電機(jī)轉(zhuǎn)速和環(huán)境溫度歷史數(shù)據(jù)，構(gòu)建樣本集{(Vi,Pi,Ni,Ti，環(huán),Ti)}，i=1,2,…，n，其中Vi表示風(fēng)速，Pi表示有功功率，Ni表示轉(zhuǎn)速，Ti，環(huán)表示環(huán)境溫度，Ti表示軸承溫度。

(2) 相關(guān)系數(shù)和權(quán)重。

相關(guān)系數(shù)C反映了變量X(V,P,N,T環(huán))對(duì)Y(T)的相關(guān)程度，記作C(X,Y)。

(2)

將相關(guān)系數(shù)歸一化處理后，得到變量X(V,P,N,T環(huán))對(duì)Y(T)的影響權(quán)重A(a1,a2,a3,a4)，其中：

(3)

式中：aj(Xj,Y)為第j類變量對(duì)Y(T)的影響權(quán)重；k為變量的總數(shù)，即k=4。

(3) 工況劃分。

(4)

式中：σ為不同工況下溫度偏離方差。

(4) 劣化度。

(5)

2 軸承健康劣化趨勢(shì)預(yù)測(cè)

風(fēng)電機(jī)組軸承健康劣化趨勢(shì)具有明顯的非平穩(wěn)性特性，這直接影響到其預(yù)測(cè)精度，需要提前處理將劣化趨勢(shì)分解成相對(duì)平穩(wěn)的分量。

經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)可以將復(fù)雜的非平穩(wěn)信號(hào)分解為一系列相對(duì)平穩(wěn)的本征模態(tài)函數(shù)分量和一個(gè)剩余分量，一定程度上避免了人為因素對(duì)分解結(jié)果的干預(yù)。但在某些情況下，EMD會(huì)產(chǎn)生模態(tài)混疊，出現(xiàn)分解后有局部不平穩(wěn)和不連貫現(xiàn)象，因此往往得不到理想的分解結(jié)果。EEMD加入高斯白噪聲，利用高斯白噪聲頻率均勻分布的統(tǒng)計(jì)特性，使信號(hào)在不同尺度上具有連續(xù)性，避免了模態(tài)混疊現(xiàn)象[10-11]。

2.1 EEMD模態(tài)分解

將多組具有有限振幅的高斯白噪聲序列ni(t)添加到原始劣化趨勢(shì)上，共進(jìn)行m次，即

di(t)=d(t)+ni(t),i=1,2,…,m

(6)

式中：di(t)為t時(shí)段第i次加入高斯白噪聲后的劣化度。

對(duì)di(t)進(jìn)行EMD分解，得到n個(gè)本征模態(tài)函數(shù)cij(t)和一個(gè)ri(t)。其中，cij(t)表示第i次加入高斯白噪聲后，分解獲得的第j個(gè)IMF，j=1，2，…，n。

將獲得的IMF進(jìn)行總體平均運(yùn)算，得到基于EEMD的IMF：

(7)

2.2 改進(jìn)的EEMD方法

由于風(fēng)電機(jī)組軸承健康劣化趨勢(shì)受各種噪聲的影響，主要是隨機(jī)噪聲和脈沖噪聲。而EEMD分解方法僅能去除信號(hào)中的隨機(jī)噪聲，脈沖噪聲還有待處理。利用小波變換對(duì)信號(hào)進(jìn)行去噪，能濾除高頻噪聲部分，重構(gòu)低頻部分，保留有價(jià)值的信息[12]，因此在對(duì)信號(hào)進(jìn)行EEMD分解之前需進(jìn)行小波去噪。具體步驟如下：(1)選擇合適的小波基，并確定分解層數(shù)n，然后對(duì)信號(hào)進(jìn)行n層小波分解；(2)將高頻系數(shù)全部置零；(3)對(duì)第n層的低頻系數(shù)和第1～第n層經(jīng)過(guò)閾值處理后的高頻系數(shù)進(jìn)行重構(gòu)，得到去噪后的信號(hào)。

2.3 組合預(yù)測(cè)

分別疊加高頻分量和低頻分量得到組合高頻信號(hào)和組合低頻信號(hào)。應(yīng)用時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型對(duì)組合高頻信號(hào)、組合低頻信號(hào)和殘余量進(jìn)行預(yù)測(cè)，進(jìn)一步疊加各預(yù)測(cè)結(jié)果得到最終預(yù)測(cè)結(jié)果(見(jiàn)圖1)。

3 算法驗(yàn)證

3.1 數(shù)據(jù)準(zhǔn)備

以華北某風(fēng)電場(chǎng)56號(hào)風(fēng)電機(jī)組的發(fā)電機(jī)后軸承為研究對(duì)象，該風(fēng)機(jī)切入風(fēng)速為3 m/s，切出風(fēng)速為24 m/s。隨機(jī)選擇2015-06-01—2016-03-30采樣頻率為10 min的發(fā)電機(jī)后軸承溫度、風(fēng)速和有功功率的SCADA數(shù)據(jù)，篩選后共有32 595組健康的樣本數(shù)據(jù)，如表1所示。

圖1 改進(jìn)的EEMD-時(shí)間序列神經(jīng)網(wǎng)絡(luò)組合預(yù)測(cè)流程圖

Fig.1 Flow chart of improved EEMD-time series neural network combined forecasting

表1 健康數(shù)據(jù)集

3.2 相關(guān)性計(jì)算

應(yīng)用式(2)計(jì)算各參數(shù)與軸承溫度的相關(guān)系數(shù)，通過(guò)式(3)歸一化處理后的各相關(guān)系數(shù)作為各變量影響軸承劣化度的權(quán)重，結(jié)果見(jiàn)表2。

3.3 健康狀態(tài)劣化趨勢(shì)

從32 595組樣本數(shù)據(jù)中，選取32 395組數(shù)據(jù)作為建模樣本，其余的200組數(shù)據(jù)作為測(cè)試樣本。

表2各參數(shù)與軸承溫度的相關(guān)系數(shù)

Tab.2Correlationcoefficientbetweeneachparameterandbearingtemperature

參數(shù)風(fēng)速功率發(fā)電機(jī)轉(zhuǎn)速環(huán)境溫度C0.5570.5530.4530.136權(quán)重0.3280.3250.2670.080

圖2 健康狀態(tài)劣化趨勢(shì)

圖3 加權(quán)后健康狀態(tài)劣化趨勢(shì)

3.4 模態(tài)分解

從圖3可以看出，風(fēng)機(jī)軸承健康狀態(tài)劣化趨勢(shì)具有復(fù)雜的非平穩(wěn)性，為了提高狀態(tài)劣化趨勢(shì)預(yù)測(cè)的準(zhǔn)確性，利用改進(jìn)的EEMD將其分解成若干個(gè)相對(duì)平穩(wěn)的分量，分解結(jié)果見(jiàn)圖4。圖中c1～c6為不同尺度的較平穩(wěn)的IMF，c7為趨勢(shì)分量。

通過(guò)計(jì)算得到各個(gè)分量的均值(見(jiàn)表3)。從表3可以看出，各個(gè)分量的均值在c5處出現(xiàn)了較大的波動(dòng)，有明顯偏離0的趨勢(shì)，因此c1～c4為高頻分量，c5～c6為低頻分量。

3.5 預(yù)測(cè)結(jié)果

時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型能夠較好地適應(yīng)信號(hào)的非線性和非平穩(wěn)性，具有較高的預(yù)測(cè)精度。因此，最終選擇時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型來(lái)獲得各類組合分量的預(yù)測(cè)結(jié)果。

圖4 基于改進(jìn)的EEMD軸承健康狀態(tài)劣化趨勢(shì)分解結(jié)果

Fig.4 Decomposed results of degradation trend for wind turbine based on modified EEMD

表3 各個(gè)分量的均值

選用前160組數(shù)據(jù)進(jìn)行訓(xùn)練和測(cè)試，后40組數(shù)據(jù)進(jìn)行預(yù)測(cè)。各類組合分量分別進(jìn)行時(shí)間延遲L=3的時(shí)間序列神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)，組合高頻信號(hào)、組合低頻信號(hào)和殘余量的預(yù)測(cè)結(jié)果如圖5所示。從圖5可以看出，隨著軸承健康狀態(tài)劣化度疊加分量的平穩(wěn)性越好和規(guī)律性越強(qiáng)，預(yù)測(cè)精度越高。

將各分量的預(yù)測(cè)結(jié)果進(jìn)行疊加，得到風(fēng)機(jī)軸承健康狀態(tài)劣化趨勢(shì)的最終預(yù)測(cè)結(jié)果。并與其他的預(yù)測(cè)方法進(jìn)行對(duì)比，結(jié)果如圖6所示。由圖6可知，所提出的預(yù)測(cè)模型可以有效跟蹤風(fēng)機(jī)軸承劣化趨勢(shì)。

4 模型誤差分析

采用平均絕對(duì)誤差(MAE)、平均絕對(duì)百分比誤差(MAPE)、均方根誤差(RMSE)和均方根百分比誤差(MSPE)等誤差指標(biāo)來(lái)度量模型精度。同時(shí)，也為了驗(yàn)證本文模型的準(zhǔn)確性，分別采用人工神經(jīng)網(wǎng)絡(luò)(BPNN)模型、徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(RBFNN)模型、支持向量機(jī)(SVM)模型和時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型對(duì)不同機(jī)組的發(fā)電機(jī)后軸承健康狀態(tài)劣化趨勢(shì)進(jìn)行預(yù)測(cè)，未改進(jìn)的EEMD預(yù)測(cè)誤差如表4所示。

由表4可知，單看MAE和RMSE 2個(gè)誤差指標(biāo)，對(duì)于未使用EEMD方法分解的預(yù)測(cè)模型，即BPNN模型、RBFNN模型、SVM模型和時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型，通過(guò)對(duì)3臺(tái)機(jī)組的發(fā)電機(jī)后軸承狀態(tài)劣化趨勢(shì)結(jié)果來(lái)看，時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型誤差都是最小的，說(shuō)明對(duì)于較復(fù)雜的非線性和非平穩(wěn)性信號(hào)，時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型的預(yù)測(cè)精度相對(duì)其他傳統(tǒng)預(yù)測(cè)模型更高。

(a) 組合高頻信號(hào)預(yù)測(cè)結(jié)果

(b) 組合低頻信號(hào)預(yù)測(cè)結(jié)果

(c) 殘余量預(yù)測(cè)結(jié)果

圖6 最終預(yù)測(cè)結(jié)果

而在使用EEMD方法分解之后，即EEMD+時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型，3臺(tái)機(jī)組的各項(xiàng)誤差指標(biāo)都有明顯改善。從表4中的均值誤差指標(biāo)分析，EE-MD+時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型的幾個(gè)誤差值相比時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型降低了50%多，表明EEMD結(jié)合傳統(tǒng)參數(shù)預(yù)測(cè)模型的組合預(yù)測(cè)能夠有效提高預(yù)測(cè)精度。

表4 5種模型的誤差指標(biāo)

充分考慮小波去噪對(duì)信號(hào)的影響，分別對(duì)EEMD分解前的信號(hào)和EEMD分解后的各分量進(jìn)行去噪，預(yù)測(cè)誤差如表5所示。由表5可知，在對(duì)風(fēng)機(jī)軸承健康狀態(tài)劣化趨勢(shì)進(jìn)行EEMD分解之前進(jìn)行小波去噪，能提高預(yù)測(cè)精度，而對(duì)EEMD的各分量進(jìn)行小波去噪，并不能提高預(yù)測(cè)精度，反而在一定程度上降低了預(yù)測(cè)精度。

5 結(jié) 論

(1) 利用SCADA監(jiān)測(cè)風(fēng)速、功率、轉(zhuǎn)速和環(huán)境溫度參數(shù)，計(jì)算這些參數(shù)與軸承溫度(劣化度)的相關(guān)系數(shù)，歸一化處理記作該參數(shù)對(duì)軸承劣化度的影響權(quán)重。在此基礎(chǔ)上，基于溫度特征量構(gòu)建的軸承健康劣化度模型，可有效描述軸承的劣化程度。

(2) 應(yīng)用改進(jìn)后的EEMD結(jié)合時(shí)間序列神經(jīng)網(wǎng)絡(luò)模型分析風(fēng)機(jī)軸承健康狀態(tài)劣化趨勢(shì)，將具有非平穩(wěn)性特性的劣化趨勢(shì)分解為相對(duì)平穩(wěn)和規(guī)律的分量。相對(duì)于其他預(yù)測(cè)方法多進(jìn)行了2次預(yù)測(cè)，但能很好地跟蹤劣化趨勢(shì)，得到更高的預(yù)測(cè)精度。

表5改進(jìn)的EEMD組合預(yù)測(cè)模型的誤差指標(biāo)

Tab.5ErrorindexesofmodifiedEEMDcombinedforecastingmodel

機(jī)組誤差指標(biāo)EEMD+小波小波+EEMD38號(hào)MAE0.03330.0171RMSE0.04590.0327MAPE0.11220.0614MSPE0.16550.126256號(hào)MAE0.02910.0091RMSE0.03610.0109MAPE0.06720.0215MSPE0.08640.026687號(hào)MAE0.02110.0066RMSE0.03840.0129MAPE0.22320.0661MSPE1.36560.1161均值MAE0.02780.0109RMSE0.04010.0188MAPE0.13420.0497MSPE0.53920.0896

(3) 改進(jìn)的EEMD結(jié)合時(shí)間序列神經(jīng)網(wǎng)絡(luò)的組合預(yù)測(cè)模型，可更準(zhǔn)確地預(yù)測(cè)風(fēng)電機(jī)組軸承健康狀態(tài)劣化趨勢(shì)，有利于風(fēng)電場(chǎng)制定檢修策略和開(kāi)展可靠性評(píng)估。

：

[1] 李輝, 胡姚剛, 唐顯虎, 等. 并網(wǎng)風(fēng)電機(jī)組在線運(yùn)行狀態(tài)評(píng)估方法[J].中國(guó)電機(jī)工程學(xué)報(bào), 2010, 30(33): 103-109.

LI Hui, HU Yaogang, TANG Xianhu, et al. Method for on-line operating conditions assessment for a grid-connected wind turbine generator system[J].ProceedingsoftheCSEE, 2010, 30(33): 103-109.

[2] TCHAKOUA P, WAMKEUE R, OUHROUCHE M, et al. Wind turbine condition monitoring: state-of-the-art review, new trends, and future challenges[J].Energies, 2014, 7(4): 2595-2630.

[3] YANG Wenxian, TAVNER P J, CRABTREE C J, et al. Cost-effective condition monitoring for wind turbines[J].IEEETransactionsonIndustrialElectronics, 2010, 57(1):263-271.

[4] 蘇連成, 李興林, 李小俚, 等. 風(fēng)電機(jī)組軸承的狀態(tài)監(jiān)測(cè)和故障診斷與運(yùn)行維護(hù)[J].軸承, 2012(1): 47-53.

SU Liancheng, LI Xinglin, LI Xiaoli, et al. Condition monitoring and fault diagnosis and operating maintenance systems for wind turbine bearings[J].Bearing, 2012(1): 47-53.

[5] 安學(xué)利, 唐擁軍, 王允. 基于健康樣本的風(fēng)電機(jī)組滾動(dòng)軸承狀態(tài)評(píng)估[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào), 2015, 13(1):48-53.

AN Xueli, TANG Yongjun, WANG Yun. State evaluation of wind turbine roller bearing based on normal samples[J].JournalofChinaInstituteofWaterResourcesandHydropowerResearch, 2015, 13(1):48-53.

[6] 安學(xué)利, 唐擁軍, 吳光軍, 等. 基于健康樣本的風(fēng)電機(jī)組狀態(tài)評(píng)估與異常檢測(cè)系統(tǒng)[J].大電機(jī)技術(shù), 2015(2):14-17.

AN Xueli, TANG Yongjun, WU Guangjun, et al. Condition assessment and anomaly detection system of wind turbine based on healthy samples[J].LargeElectricMachineandHydraulicTurbine, 2015(2):14-17.

[7] 滕偉, 姜銳, 張陽(yáng)陽(yáng), 等. 基于多尺度包絡(luò)譜圖的直驅(qū)式風(fēng)電機(jī)組軸承故障特征提取[J].動(dòng)力工程學(xué)報(bào), 2017, 37(5): 373-378.

TENG Wei, JIANG Rui, ZHANG Yangyang, et al.Bearing fault feature extraction for a direct drive wind turbine using multi-scale enveloping spectrogram[J].JournalofChineseSocietyofPowerEngineering, 2017,37(5): 373-378.

[8] 李輝, 胡姚剛, 李洋, 等. 基于溫度特征量的風(fēng)電機(jī)組關(guān)鍵部件劣化漸變概率分析[J].電力自動(dòng)化設(shè)備, 2015, 35(11): 1-7, 19.

LI Hui, HU Yaogang, LI Yang, et al. Gradual deterioration probability analysis based on temperature characteristic parameters for critical components of wind turbine generator system[J].ElectricPowerAutomationEquipment, 2015, 35(11): 1-7, 19.

[9] 劉長(zhǎng)良, 閆蕭. 基于工況辨識(shí)和變分模態(tài)分解的風(fēng)電機(jī)組滾動(dòng)軸承故障診斷[J].動(dòng)力工程學(xué)報(bào), 2017, 37(4): 273-278, 334.

LIU Changliang, YAN Xiao. Fault diagnosis of a wind turbine rolling bearing based on variational mode decomposition and condition identification[J].JournalofChineseSocietyofPowerEngineering, 2017, 37(4): 273-278, 334.

[10] 王賀, 胡志堅(jiān), 張翌暉, 等. 基于聚類經(jīng)驗(yàn)?zāi)B(tài)分解和最小二乘支持向量機(jī)的短期風(fēng)速組合預(yù)測(cè)[J].電工技術(shù)學(xué)報(bào), 2014, 29(4): 237-245.

WANG He, HU Zhijian, ZHANG Yihui, et al.A hybrid model for short-term wind speed forecasting based on ensemble empirical mode decomposition and least squares support vector machines[J].TransactionsofChinaElectrotechnicalSociety, 2014, 29(4): 237-245.

[11] AN Xueli, JIANG Dongxiang, LI Shaohua, et al. Application of the ensemble empirical mode decomposition and Hilbert transform to pedestal looseness study of direct-drive wind turbine[J].Energy, 2011, 36(9): 5508-5520.

[12] 崔玉敏, 和衛(wèi)星, 高彩虹. 基于小波變換及HHT的軸承故障診斷應(yīng)用研究[J].計(jì)算機(jī)仿真, 2010, 27(10): 294-296, 300.

CUI Yumin, HE Weixing, GAO Caihong. Application of bearing fault diagnosis based on wavelet and HHT[J].ComputerSimulation, 2010, 27(10): 294-296, 300.