蘇俊 方超 孫猛
摘 ?要:傳統(tǒng)風(fēng)電機(jī)組故障診斷預(yù)警算法準(zhǔn)確度較低,預(yù)警速度較慢,為此,提出基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷預(yù)警算法。首先進(jìn)行故障診斷數(shù)據(jù)預(yù)處理,基于工況辨識(shí),選擇合適的風(fēng)電機(jī)組運(yùn)行參數(shù),采用溫度預(yù)測(cè)的方法,利用觀測(cè)溫度與預(yù)期溫度的殘差來(lái)反映設(shè)備的運(yùn)行狀態(tài),使用MSET建立風(fēng)機(jī)齒輪箱溫度模型,確定預(yù)警指標(biāo),設(shè)計(jì)基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷預(yù)警算法。實(shí)驗(yàn)結(jié)果表明,設(shè)計(jì)算法能有效減少診斷預(yù)警時(shí)間,且其準(zhǔn)確度更高。
關(guān)鍵詞:工況辨識(shí);風(fēng)電機(jī)組;故障診斷;預(yù)警算法
中圖分類(lèi)號(hào):TP277 ? 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):2096-4706(2021)12-0033-04
Abstract: The traditional fault diagnosis and early warning algorithm of wind turbine has low accuracy and slow early warning speed. Therefore, a fault diagnosis and early warning algorithm of wind turbine based on condition identification is proposed. Firstly, preprocess the fault diagnosis data, select the appropriate operating parameters of the wind turbine based on the condition identification, use the temperature prediction method to reflect the operating state of the equipment through the residual between the observed temperature and the expected temperature, use MSET to establish the temperature model of the wind turbine gearbox, determine the early warning index, and design the fault diagnosis and early warning algorithm of the wind turbine based on the condition identification. Experimental results show that the designed algorithm can effectively reduce the diagnosis and early warning time, and its accuracy is higher.
Keywords: condition identification; wind turbine; fault diagnosis; early warning algorithm
0 ?引 ?言
國(guó)家日益重視環(huán)境保護(hù),大力支持新能源的發(fā)展,風(fēng)力發(fā)電就是一種基于新能源的發(fā)電技術(shù),其可以利用環(huán)境的自然特性產(chǎn)生并供給電力,可以極大地減小能源消耗,減輕環(huán)境污染。由于其這些優(yōu)點(diǎn),目前已經(jīng)成為各國(guó)的主流發(fā)電方式[1]。但由于風(fēng)力發(fā)電的特殊性,某些地區(qū)雖具有豐富的風(fēng)力資源,但由于其自然環(huán)境惡劣導(dǎo)致機(jī)組設(shè)備經(jīng)常受損,風(fēng)力發(fā)電機(jī)故障頻發(fā),直接影響風(fēng)力發(fā)電的效率,也對(duì)風(fēng)力發(fā)電的安全性造成威脅,因此需要提前對(duì)風(fēng)電機(jī)組進(jìn)行故障診斷預(yù)警,從而快速診斷風(fēng)電機(jī)組故障,提升發(fā)電效率。
因此,在這種情況下,相關(guān)研究人員研究風(fēng)電機(jī)組故障產(chǎn)生的原因[2],并根據(jù)故障原因利用傳統(tǒng)的檢測(cè)技術(shù)將故障信號(hào)進(jìn)行轉(zhuǎn)換,然后提取故障信息特征,并處理故障信息,從而實(shí)現(xiàn)風(fēng)電機(jī)組的故障預(yù)警。傳統(tǒng)的故障預(yù)警方法主要可以預(yù)警中小規(guī)模的風(fēng)電機(jī)組,當(dāng)對(duì)大規(guī)模風(fēng)電機(jī)組進(jìn)行預(yù)警時(shí),由于傳統(tǒng)方法使用的算法過(guò)程較為簡(jiǎn)單,不能滿(mǎn)足故障診斷和故障預(yù)警的要求。因此需要在確保故障預(yù)測(cè)準(zhǔn)確性的同時(shí)發(fā)出預(yù)警,進(jìn)行實(shí)時(shí)性檢測(cè),基于此,本文提出一種基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷預(yù)警模型,并在現(xiàn)有風(fēng)電機(jī)組故障診斷預(yù)警算法的基礎(chǔ)上,設(shè)計(jì)一種新的故障診斷預(yù)警算法。
1 ?基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷預(yù)警算法設(shè)計(jì)
1.1 ?進(jìn)行故障診斷數(shù)據(jù)預(yù)處理
故障處理預(yù)警平臺(tái)記錄的風(fēng)電機(jī)組運(yùn)行數(shù)據(jù)量非常多,包含風(fēng)電機(jī)組在各種異常運(yùn)行情況下的數(shù)據(jù)點(diǎn),這些異常數(shù)據(jù)點(diǎn)影響研究結(jié)果的可靠性。MSET模型必須輸入風(fēng)電機(jī)組正常運(yùn)行期間的健康數(shù)據(jù),以便刪除這些異常數(shù)據(jù)點(diǎn)[3]。本文在特定時(shí)間選擇后選擇適當(dāng)數(shù)量的數(shù)據(jù),檢查所選數(shù)據(jù)是否滿(mǎn)足風(fēng)電機(jī)組的全部運(yùn)行條件,并去除設(shè)備輸出為零或負(fù)值或風(fēng)速為零的樣本點(diǎn)。如果風(fēng)速低于切入風(fēng)速并大于切出風(fēng)速,則這些數(shù)據(jù)滿(mǎn)足異常運(yùn)行條件下的數(shù)據(jù)點(diǎn)要求。
同時(shí),不同運(yùn)行參數(shù)的數(shù)據(jù)維度不同,最大值和最小值相差很大[4]。例如,風(fēng)速范圍通常為3 m/s~25 m/s,發(fā)電機(jī)轉(zhuǎn)速范圍通常為1 000 rpm~1 800 rpm,絕對(duì)值差超過(guò)10倍,忽略這種情況會(huì)導(dǎo)致建模時(shí)出現(xiàn)不匹配現(xiàn)象,進(jìn)而影響模型的計(jì)算精度。為此,減少數(shù)據(jù)之間不同維度的影響,需要將各個(gè)運(yùn)行參數(shù)的值控制在統(tǒng)一范圍內(nèi),歸一化公式如式(1)所示:
1.2 ?選取風(fēng)電機(jī)組故障診斷預(yù)警參數(shù)
在風(fēng)電機(jī)組記錄的數(shù)據(jù)中,齒輪箱溫度、滾動(dòng)軸承振動(dòng)對(duì)風(fēng)電機(jī)組的相關(guān)參數(shù)產(chǎn)生一定影響,其數(shù)據(jù)范圍與風(fēng)電機(jī)組的故障狀況以及運(yùn)行時(shí)間相關(guān)[5]。因此,在確定風(fēng)電機(jī)組的運(yùn)行狀態(tài)和合理劃分風(fēng)電機(jī)組的運(yùn)行工況之前,需要選擇合適的風(fēng)電機(jī)組運(yùn)行參數(shù)。所選參數(shù)直接或間接地影響風(fēng)力發(fā)電機(jī)的狀況,也間接證明風(fēng)電機(jī)組的故障范圍。從設(shè)備振動(dòng)頻率、齒輪箱溫度、設(shè)備SCADA、振動(dòng)加速度、風(fēng)速、功率、發(fā)電機(jī)轉(zhuǎn)矩之間的關(guān)系可知,風(fēng)速、功率和發(fā)電機(jī)轉(zhuǎn)矩對(duì)齒輪箱溫度和檢測(cè)設(shè)備有明顯的影響。
利用風(fēng)電機(jī)組發(fā)電時(shí),其預(yù)警參數(shù)也與風(fēng)力發(fā)電機(jī)的類(lèi)型有關(guān),變速器的風(fēng)力發(fā)電機(jī)組主要受到其轉(zhuǎn)動(dòng)速度影響,同時(shí),其與轉(zhuǎn)動(dòng)加速度關(guān)系密切,一般呈正比例關(guān)系。即轉(zhuǎn)動(dòng)加速度越高,預(yù)警參數(shù)范圍越大。變速箱轉(zhuǎn)速越高其溫度越高,當(dāng)風(fēng)電機(jī)組轉(zhuǎn)速增加時(shí),齒輪箱承受的載荷也隨之增加,發(fā)電機(jī)滾動(dòng)軸承轉(zhuǎn)速也在增加,振動(dòng)越來(lái)越強(qiáng),齒輪箱溫度也隨之升高。
1.3 ?建立MSET故障診斷預(yù)警模型
目前,風(fēng)電發(fā)電機(jī)組故障預(yù)警數(shù)據(jù)監(jiān)測(cè)與分析,利用溫度預(yù)測(cè)方法,其主要是利用觀測(cè)溫度與預(yù)期溫度之間的剩余量來(lái)反映裝置的運(yùn)行狀態(tài)[6]。如果剩余誤差超過(guò)設(shè)定值,觸發(fā)報(bào)警裝置,使用系統(tǒng)誤差溫度預(yù)測(cè)對(duì)實(shí)際溫度進(jìn)行準(zhǔn)確估計(jì)。通過(guò)MSET建模方法確定風(fēng)機(jī)齒輪箱溫度模型,并使用SCADA數(shù)據(jù)之間的復(fù)雜相關(guān)性,描述設(shè)備運(yùn)行期間齒輪箱溫度和其他參數(shù)之間的隱含關(guān)系,以獲得準(zhǔn)確的預(yù)警值。根據(jù)溫度觀測(cè)值與估計(jì)值之間殘差的統(tǒng)計(jì)性質(zhì),需要確定合適的殘差誤差值,以確定齒輪箱誤差的預(yù)警范圍。基于工況辨識(shí)故障預(yù)警流程如圖1所示。
如圖1所示。左側(cè)部分使用MSET進(jìn)行工況拆分建立正確的訓(xùn)練模型,設(shè)置每個(gè)子任務(wù)的殘差閾值指標(biāo)并發(fā)布預(yù)警信息異常率指標(biāo),右側(cè)部分確定觀測(cè)數(shù)據(jù)的工作狀態(tài)并計(jì)算相應(yīng)的工作狀態(tài)指數(shù)用于確定異常值,并使用滑動(dòng)窗口方法計(jì)算相應(yīng)時(shí)間內(nèi)異常值的比率。
故障預(yù)警模型使用MSET技術(shù),因其主要在預(yù)估領(lǐng)域有巨大作用,因此將其應(yīng)用在故障預(yù)警模型的建立中,可以識(shí)別風(fēng)電機(jī)組各個(gè)設(shè)備的故障信息,并將其合并到模型中進(jìn)行統(tǒng)一反饋。除此之外,其基于數(shù)據(jù)驅(qū)動(dòng)的有線性模型和非線性函數(shù)。線性建模方法在風(fēng)電機(jī)組故障傳感器的標(biāo)定、系統(tǒng)狀態(tài)監(jiān)測(cè)、故障報(bào)警和診斷等方面得到了廣泛而有效的應(yīng)用。使用模型估算和分析靜態(tài)運(yùn)行歷史數(shù)據(jù)中運(yùn)行參數(shù)之間的相關(guān)性[7],并通過(guò)在包含相關(guān)參數(shù)的變量之間創(chuàng)建非線性模型來(lái)估算設(shè)備的運(yùn)行狀態(tài)。
建立MSET故障診斷預(yù)警模型的前提是創(chuàng)建適當(dāng)?shù)倪^(guò)程矩陣D,即利用故障預(yù)警檢測(cè)參數(shù),建立恰當(dāng)?shù)倪^(guò)程矩陣,以覆蓋風(fēng)電機(jī)組的各個(gè)設(shè)備的運(yùn)行范圍,運(yùn)行中的歷史數(shù)據(jù)能夠真實(shí)地反映風(fēng)電機(jī)組的運(yùn)行狀態(tài),其組成的變量集之間是相互關(guān)聯(lián)并耦合的。由相互關(guān)聯(lián)的變量組成的觀察向量可以用來(lái)描述設(shè)備的運(yùn)行狀態(tài)。設(shè)定風(fēng)電機(jī)組運(yùn)行過(guò)程中有m個(gè)狀態(tài),則過(guò)程矩陣D可定義為:
D=[X(1) X(2)…X(m)] ? ? ? ? ? ? ? ? (2)
其中,建立的過(guò)程矩陣可以反映風(fēng)電機(jī)組在各個(gè)時(shí)期的運(yùn)行變化,除此之外,風(fēng)電機(jī)組在運(yùn)行過(guò)程中還需要將轉(zhuǎn)換變量提取出來(lái),并計(jì)算轉(zhuǎn)換參數(shù),將矩陣結(jié)構(gòu)重組,帶入選取的參數(shù),進(jìn)行二次轉(zhuǎn)換。
MSET建模的主要步驟是確定操作參數(shù)。其算法的核心在于歷史SCADA數(shù)據(jù)組成的MSET狀態(tài)矩陣,建模選擇的數(shù)據(jù)具備的正常運(yùn)行條件是覆蓋風(fēng)力發(fā)電機(jī)組,以及所選擇的運(yùn)行參數(shù)是否與溫度變化有很強(qiáng)的相關(guān)性是選擇的核心。變速箱型號(hào)也對(duì)該模型有一定的影響,可以選擇對(duì)輸出影響較大的不同類(lèi)型的變量作為模型的輸入變量。當(dāng)故障預(yù)警模型的變量因素較多時(shí),所有變量均可作為模型的輸入?yún)?shù),這明顯增加了模型的復(fù)雜度,導(dǎo)致模型的計(jì)算量顯著增加,增加計(jì)算時(shí)間,降低精度。模型的主觀篩選具有局限性。采用灰色關(guān)聯(lián)度(GRA),選擇合理的運(yùn)行參數(shù)作為一組模型變量,應(yīng)用灰色關(guān)聯(lián)度計(jì)算各種輸入變量,可以實(shí)現(xiàn)模型輸入變量的有效、全面、綜合選擇,無(wú)須人工主觀篩選,提高模型的適應(yīng)性,實(shí)現(xiàn)定性驗(yàn)證分析基于此建立的故障診斷預(yù)警模型W如式(3)所示。
1.4 ?制定風(fēng)電機(jī)組故障預(yù)警指標(biāo)
使用MSET計(jì)算預(yù)警指標(biāo),確定兩者之間的殘差[8],當(dāng)殘差在閾值范圍內(nèi),則判斷該殘差為一個(gè)穩(wěn)定點(diǎn),如果殘差過(guò)大,則該殘差不能判斷為一個(gè)點(diǎn)。閾值指標(biāo)是在MSET模型離線訓(xùn)練期間確定的。其余部分用作訓(xùn)練數(shù)據(jù),輸入模型用于計(jì)算殘值。為了確定預(yù)警過(guò)程中的異常值,評(píng)價(jià)預(yù)警指標(biāo)。由于非參數(shù)模型在預(yù)測(cè)模型輸入變量的過(guò)程中,存在不確定性。如果閾值較低,說(shuō)明異常檢測(cè)算法對(duì)正常運(yùn)行狀態(tài)信息過(guò)于敏感,易造成錯(cuò)誤判斷。在確定殘差的平均值時(shí),需要對(duì)閾值賦予一定的系數(shù),確定預(yù)警指標(biāo)E如式(4)所示。
依據(jù)風(fēng)機(jī)運(yùn)行的隨機(jī)性能和擾動(dòng)性能可知,若運(yùn)行過(guò)程中出現(xiàn)異常波動(dòng)均屬于正?,F(xiàn)象。如果呈現(xiàn)下降趨勢(shì),則可以判斷當(dāng)前風(fēng)力發(fā)電機(jī)運(yùn)行狀況不正常,運(yùn)行過(guò)程存在隱患。因此,當(dāng)判斷某一特定點(diǎn)或少數(shù)點(diǎn)異常時(shí),立即發(fā)出預(yù)警信息,觸發(fā)誤報(bào),增加誤報(bào)率。因此,為降低誤報(bào)率,補(bǔ)充異常率作為預(yù)警指標(biāo),發(fā)布預(yù)警信息。假設(shè)在一段時(shí)間內(nèi)連續(xù)計(jì)算N個(gè)變速箱溫度的觀測(cè)和估計(jì)殘差,可以將其中的Na判斷為異常值。
1.5 ?故障診斷預(yù)警算法
在風(fēng)力發(fā)電機(jī)組的實(shí)際運(yùn)行中,依據(jù)實(shí)際運(yùn)行條件確定殘差閾值,結(jié)合溫度模型試驗(yàn)結(jié)果和以往經(jīng)驗(yàn)建立殘差閾值。利用任務(wù)條件識(shí)別中處理的有效數(shù)據(jù)構(gòu)建MSET模型歷史矩陣,將通過(guò)任務(wù)條件識(shí)別得到的4個(gè)任務(wù)條件所包含的數(shù)據(jù)作為驗(yàn)證樣本。使用設(shè)計(jì)的方法,將所有未分離工況的數(shù)據(jù)和四種工況的數(shù)據(jù)作為驗(yàn)證樣本,設(shè)置溫度估計(jì)模型,計(jì)算殘差平均值,確定相應(yīng)的預(yù)警閾值。根據(jù)該閾值,結(jié)合式(1)~(4),設(shè)計(jì)的基于工況辨識(shí)的故障診斷預(yù)警算法U如式(5)所示。
式(5)中,B代表元素差值,C代表控制范圍,使用設(shè)計(jì)的算法將風(fēng)電機(jī)組實(shí)際運(yùn)行工況進(jìn)行劃分,將空間數(shù)據(jù)單獨(dú)存儲(chǔ),作為后續(xù)工作的故障預(yù)警數(shù)據(jù)源,為基于工況辨識(shí)下的風(fēng)電機(jī)組故障預(yù)警提供基礎(chǔ)。
2 ?實(shí)驗(yàn)
為檢測(cè)本文設(shè)計(jì)的風(fēng)電機(jī)組故障診斷預(yù)警算法的準(zhǔn)確度,使用傳統(tǒng)的預(yù)警算法和本文設(shè)計(jì)的預(yù)警算法進(jìn)行對(duì)比實(shí)驗(yàn),檢測(cè)其與標(biāo)準(zhǔn)的故障預(yù)警齒輪軸溫的區(qū)別,風(fēng)電機(jī)組故障預(yù)警系統(tǒng)程序運(yùn)行如圖2所示。
2.1 ?實(shí)驗(yàn)準(zhǔn)備
各節(jié)點(diǎn)的硬件配置為:CPU其內(nèi)存為128 GB,硬盤(pán)360 GB,故障截取的并行比用式(6)計(jì)算。
式(6)中,F(xiàn)代表并行比,H代表故障節(jié)點(diǎn)數(shù),T代表預(yù)警反應(yīng)時(shí)間,實(shí)驗(yàn)用前置機(jī)及服務(wù)器獲取的數(shù)據(jù),整合為XML格式,根據(jù)并行比公式,列出故障樣本的截取情況如表1所示。
根據(jù)表1截取的故障樣本,將多次實(shí)驗(yàn)的平均值進(jìn)行對(duì)比。
2.2 ?實(shí)驗(yàn)結(jié)果與討論
記錄本文設(shè)計(jì)算法與傳統(tǒng)算法故障診斷預(yù)警的齒輪軸溫,與標(biāo)準(zhǔn)的故障診斷預(yù)警齒輪軸溫對(duì)比,繪制曲線擬合圖,如圖3所示。
由圖3可知,本文設(shè)計(jì)的基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷預(yù)警算法診斷的故障齒輪軸溫與標(biāo)準(zhǔn)值擬合度高,能有效減少預(yù)警時(shí)間,優(yōu)于傳統(tǒng)的故障診斷預(yù)警方法。
3 ?結(jié) ?論
綜上所述,風(fēng)力能源的使用對(duì)維持我國(guó)生態(tài)環(huán)境可持續(xù)發(fā)展有重要意義,因此維持風(fēng)電機(jī)組的正常工作狀態(tài),使其可以穩(wěn)定地工作,對(duì)能源保護(hù)有重要意義,因此本文設(shè)計(jì)了基于工況辨識(shí)的風(fēng)電機(jī)組故障診斷方法,實(shí)驗(yàn)證明其與標(biāo)準(zhǔn)數(shù)值擬合度高,優(yōu)于傳統(tǒng)的故障診斷預(yù)警方法,具有一定的應(yīng)用價(jià)值,為后續(xù)的故障診斷預(yù)警提供參考。
參考文獻(xiàn):
[1] 靳志杰,霍志紅,許昌,等.基于特征選擇和XGBoost的風(fēng)電機(jī)組故障診斷 [J].可再生能源,2021,39(3):353-358.
[2] 鄧子豪,李錄平,劉瑞,等.基于SCADA數(shù)據(jù)特征提取的風(fēng)電機(jī)組偏航齒輪箱故障診斷方法研究 [J].動(dòng)力工程學(xué)報(bào),2021,41(1):43-50.
[3] 王紹平,王冰,曹智杰,等.海上風(fēng)電機(jī)組分布式控制系統(tǒng)故障診斷與定位研究 [J].可再生能源,2020,38(10):1343-1348.
[4] 孫群麗,周瑛,劉長(zhǎng)良.基于LARS特征選擇的風(fēng)電機(jī)組故障診斷的研究 [J].可再生能源,2020,38(10):1349-1354.
[5] 常興邦,詹俊,盛利.基于粒子濾波與BP神經(jīng)網(wǎng)絡(luò)的風(fēng)電機(jī)組主軸承溫度故障診斷研究 [J].技術(shù)與市場(chǎng),2020,27(6):11-15.
[6] 高峰,鄧星星,劉強(qiáng),等.大型風(fēng)電機(jī)組電動(dòng)變槳系統(tǒng)變槳角度故障診斷 [J].太陽(yáng)能學(xué)報(bào),2020,41(5):98-106.
[7] 鄧子豪,李錄平,劉瑞,等.大型風(fēng)電機(jī)組電機(jī)驅(qū)動(dòng)型主動(dòng)偏航系統(tǒng)故障診斷技術(shù)概述 [J].太陽(yáng)能,2020(4):34-41.
[8] 孫鶴旭,孫澤賢,張靖軒.基于并行模糊C-均值聚類(lèi)的風(fēng)電機(jī)組發(fā)電機(jī)故障診斷研究 [J].太陽(yáng)能學(xué)報(bào),2020,41(3):8-14.
作者簡(jiǎn)介:蘇俊(1983.01—),男,漢族,江蘇鹽城人,中級(jí)工程師,碩士研究生,研究方向:嵌入式開(kāi)發(fā)、風(fēng)機(jī)健康管理。