劉桂然

（國電聯(lián)合動力技術(shù)有限公司，風電設(shè)備及控制國家重點實驗室，北京100039）

1 概述

主軸軸承性能直接影響到風電機組可靠性，設(shè)計、材料、熱處理、潤滑等任何一個環(huán)節(jié)的偏差都會導致軸承失效，而如何在主軸軸承發(fā)生早期失效時進行故障診斷和預(yù)測，是風電行業(yè)需要亟待解決的一個重要課題。

人工神經(jīng)網(wǎng)絡(luò)是適應(yīng)工程需要發(fā)展起來的一門交叉學科，近年來人工神經(jīng)網(wǎng)絡(luò)已開始用于風電領(lǐng)域，進行軸承、齒輪箱等部件的故障，本文介紹了一種基于BP- 神經(jīng)網(wǎng)絡(luò)的風電主軸軸承溫度在線預(yù)測和診斷的方法。

2 BP- 神經(jīng)網(wǎng)絡(luò)算法

2.1 變量的定義

在三層BP- 神經(jīng)網(wǎng)絡(luò)中，假定輸入神經(jīng)元個數(shù)為M，隱含層神經(jīng)元個數(shù)為I，輸出層神經(jīng)元個數(shù)為J，輸入層第m 個神經(jīng)元為xm，隱含層第i 個神經(jīng)元為ki，輸出層第j 個神經(jīng)元為yj，從xm到ki的連接權(quán)值為ωmi，從ki到y(tǒng)i的連接權(quán)值記ωij，隱含層傳遞函數(shù)為sigmoid 函數(shù)，輸出層傳遞函數(shù)采用線性函數(shù)，網(wǎng)絡(luò)實際輸出為Y（n）=[諄1J，諄2J…，諄JJ]，期望輸出為d（n）=[d1，d2，…，dJ]，n 為迭代次數(shù)，第n 次迭代誤差為ej（n）=dj（n）-Yj（n）。

2.2 工作信號正向傳播

2.3 網(wǎng)絡(luò)權(quán)值調(diào)整

權(quán)值調(diào)整采用LM算法，LM算法是利用標準的數(shù)值優(yōu)化技術(shù)的快速算法，權(quán)值修正公式：

J 為雅克比矩陣，I 為單位矩陣，η 為學習率參數(shù)，介于0-1 之間。

3 故障模型建立

BP- 神經(jīng)網(wǎng)絡(luò)是前向神經(jīng)網(wǎng)絡(luò)的核心，廣泛應(yīng)用于分類識別、逼近和回歸，BP- 神經(jīng)網(wǎng)絡(luò)由1 個輸入層、若干隱含層和1 個輸出層組成，每層均由若干神經(jīng)元組成，各相鄰神經(jīng)元之間多為全連接。

3.1 故障模型

主軸軸承包括浮動軸承和止推軸承，浮動軸承只承受徑向力，止推軸承需要同時承載軸向力和徑向徑向力，止推軸承，主軸軸承失效主要指止推軸承的失效。溫度是反映軸承運行狀態(tài)的最直接的指標，當軸承存在異常磨損、潤滑不足或裝配問題時，會首先表現(xiàn)為軸承溫度異常，主軸軸承故障主要以止推軸承為主，本模型以止推軸承溫度Tz作為預(yù)測輸出變量。

影響止推軸承溫度的因素包括發(fā)電功率、運行時間和環(huán)境溫度等，發(fā)電功率P 直接影響軸承溫度的變化，作為模型的一個輸入變量。浮動軸承溫度Tf作為模型第二個輸入變量，原因在于浮動軸承溫度已經(jīng)包含了運行時間和環(huán)境溫度的影響。

變量定義：

基于上述討論，建立止推軸承溫度預(yù)測BP- 網(wǎng)絡(luò)模型如圖1所示，其中：

X1=（P1，P2，P3，…Pn），表示發(fā)電機功率P

X2=（Tf1，Tf2，Tf3…Tfn），表示浮動軸承溫度Tf

Y=（Tz1，Tz2，Tz3…Tzn），表示待預(yù)測的止推軸承溫度Tz

3.2 樣本處理

模型以某2MW機組為研究對象，采集過去1 年的歷史數(shù)據(jù)作為模型訓練樣本。為避免功率和運行時間離散因素對預(yù)測結(jié)果造成干擾，樣本只考慮滿發(fā)功率時的軸承狀態(tài)，對低于額定功率的狀態(tài)數(shù)據(jù)進行剔除，獲得有效樣本數(shù)據(jù)共22763 組。

3.3 模型訓練

按功能將上述樣本隨機分為三類，其中訓練樣本占70%，驗證樣本和測試樣本各15%。訓練完成后進行模型精度分析，包括均方誤差（MSE）和R 值分析。

圖1 網(wǎng)絡(luò)預(yù)測模型

R 用于衡量目標數(shù)據(jù)（期望輸出）與實際輸出之間的相關(guān)性，R 為1 說明兩者完全相符，R 為0 說明數(shù)據(jù)完全隨機，從圖2 可以看到三類數(shù)據(jù)的R 值均為0.93 左右，說明實際輸出和期望輸出相關(guān)性較好。MSE 值用于衡量目標數(shù)據(jù)（期望輸出）與實際輸出之間的偏差，MSE 越小，偏差越小，MSE 為0，說明沒有誤差。從圖2 可知，模型MSE 結(jié)果為5.8 左右，但從圖3 的誤差直方圖來看，誤差在5℃以內(nèi)的數(shù)據(jù)占比95.3%，從風電主軸軸承溫度預(yù)測角度來說，5℃在可接受誤差范圍之內(nèi)。

圖2 MSE 和R

圖3 誤差直方圖

4 故障預(yù)測

4.1 溫度預(yù)測

模型訓練完成后，從某風場機組中隨機抽選做為預(yù)測對象，采集該機組1 個月內(nèi)功率滿發(fā)數(shù)據(jù)，并進行止推軸承溫度預(yù)測，預(yù)測結(jié)果如圖4 所示。A6MBGB為實際采集的止推軸承溫度數(shù)據(jù)的真實值，yA6 為采用所建立的BP- 神經(jīng)網(wǎng)絡(luò)預(yù)測的到的結(jié)果，預(yù)測值與真實值趨勢一致，雖然存在一定偏差，但曲線的基本走勢和趨勢基本一致。

4.2 偏差分析

對于不同機組來說，除軸承制造參數(shù)差異外，裝配、潤滑等多種因素都會影響軸承的預(yù)測結(jié)果，造成結(jié)果偏差，整體來看預(yù)測值更接近真實值。

圖4 同風場機組溫度曲線預(yù)測

對預(yù)測修正曲線進行回歸分析，見圖5。可以看到R 值為0.90881，預(yù)測輸出與期望輸出具有很好的一致性，預(yù)測誤差在5℃以上的數(shù)據(jù)有30 個左右，占預(yù)測樣本總數(shù)的3%。97%的數(shù)據(jù)在期望范圍之內(nèi)，預(yù)測精度可以滿足要求。

圖5 回歸分析

5 結(jié)論

本文以BP 神經(jīng)網(wǎng)絡(luò)為基礎(chǔ)，提出了一種風電主軸軸承故障預(yù)測和診斷方法，基于機組歷史數(shù)據(jù)，構(gòu)建了風電主軸軸承的溫度預(yù)測模型，并成功實現(xiàn)了對其它機組的溫度預(yù)測。

本模型的實際價值在于可實現(xiàn)對風電軸承早期故障的預(yù)測，對軸承潛在問題的進行預(yù)警，提示對軸承和機組進行檢查，并采取相應(yīng)的解決方案，如采取軸承清洗、換油，調(diào)整軸承受力狀態(tài)等措施。進而解決軸承高溫問題，降低軸承故障帶來軸承更換損失，及停機造成的發(fā)電量損失。

總體來看，該模型的預(yù)測精度滿足實際工程要求，是一種理想的風電主軸軸承溫度預(yù)測和診斷方法，可用于風場中主軸軸承的故障診斷和應(yīng)用。