顧衛(wèi)祥,王先發(fā),徐靖楠,趙冰冰,童建強
(1.南京工程學院,江蘇 南京 211167;2.國網(wǎng)安徽省電力有限公司東至縣供電公司,安徽 池州 247200)
隨著經(jīng)濟的發(fā)展,世界各國對能源的需求與日俱增,大力推廣可再生能源發(fā)電已成為共識,我國風能資源豐富,是未來可再生能源的重要組成部分。目前,我國風力發(fā)電初具規(guī)模,裝機容量與日俱增。風力發(fā)電在迅速發(fā)展的同時,也伴隨著巨大的挑戰(zhàn)。風電機組內部結構相對比較復雜,工作環(huán)境比較惡劣,因此風電機組故障頻發(fā),發(fā)電可靠性能得不到保證,運維成本高。為了提高電力系統(tǒng)的穩(wěn)定性,提高風電機組發(fā)電的可靠性,減少運維成本,目前很多學者致力于風電機組故障的研究[1-2]?,F(xiàn)階段的研究大都是針對齒輪箱故障判斷尋求最優(yōu)算法,例如改進灰色算法和時間序列法[3],相關向量機和遺傳算法結合[4],改進的 Hilbert 變換包絡解調法[5],階次 分析法[6]以及基于錦標賽選擇策略的雙種群人工蜂群算法[7]等,但是都存在一定的局限性。文獻[8]提出了一種基于馬爾科夫鏈并利用歷史檢查維護數(shù)據(jù)的齒輪箱軸承狀態(tài)評估和剩余壽命預測方法,該方法降低了對大量監(jiān)測數(shù)據(jù)的依賴性且不需要安裝在線監(jiān)測設備,但是準確率有待提高。文獻[9]用人工免疫算法對BP神經(jīng)網(wǎng)絡進行改進,用于風機齒輪箱故障趨勢的預測,以單位時間內的維護成本率最小為優(yōu)化目標,同時在維護模型中考慮了維護和置換所用時間以及回復改善因子,進行優(yōu)化求解,但是其優(yōu)化目標較為單一具有局限性。文獻[10]提出了一種基于隨機子空間識別方法的齒輪箱故障預測算法,所提方法雖然能夠預測出故障的發(fā)生,但在識別出具體故障方面還需要完善。
針對風機齒輪箱軸承故障診斷,提出了一種將離散粒子群優(yōu)化算法 (Discrete Particle Swarm Optimization,DPSO)與基于多尺度小波核函數(shù)的核極限學習機 (Multi-scale Kernel Extreme Learning Machine,MKELM)相結合的新型DPSO-MKELM算法,并進行了實驗分析,結果表明該算法具有更高的分類精度和較快的收斂速度。
粒子群優(yōu)化算法(Particle Swarm Optimization,PSO)起源于對鳥類捕食行為的研究,于1995年被Kennedy 和 Eberhart兩位博士首次提出[11]。
PSO通過模擬鳥類捕食行為解決優(yōu)化模型,通過食物位置確定優(yōu)化模型的最優(yōu)解,通過鳥群和食物的距離確定尋優(yōu)的能力。PSO將鳥群的個體都看作是粒子,即潛在的解。每個粒子在多維空間飛行并且粒子的好壞由適應度函數(shù)決定。粒子在飛行過程中根據(jù)自身的飛行經(jīng)驗以及周圍優(yōu)質粒子的飛行特征,不斷調整飛行狀態(tài),從而達到群體最優(yōu)狀態(tài),最終找到的全局最優(yōu)解[12]。該算法易于實現(xiàn)、收斂速度快、優(yōu)化能力強,被廣泛用于各種領域。
1.2.1 改進權重因子
權重因子對PSO算法性能的影響很大。如果較大時,可以增強粒子全局搜索能力,有利于增加種群的多樣性[13];如果較小時,可以增強粒子局部搜索能力,有利于改善算法的收斂性。在傳統(tǒng)PSO算法中,取在(0,1)之間[14],不利于增加種群的多樣性和改善算法的收斂性。因此,對權重因子的改進如下:
1.2.2 改進學習因子
學習因子c1和c2分別代表粒子的自我學習和社會學習,分別控制算法的全局搜索和局部搜索。使用動態(tài)學習因子 c1、c2,其中 c1和 c2之和不大于 4[15]。在算法搜索前期,設置的c1值大于c2值,粒子主要依靠自身的歷史經(jīng)驗對全局進行搜索,避免陷入局部最優(yōu)狀態(tài);在搜索后期,設置的值小于值,粒子主要依靠種群的優(yōu)勢信息進行搜索,最終粒子聚攏到全局最優(yōu)區(qū)域。因此,對學習因子c1、c2的改進如下:
1.2.3 適應度函數(shù)
適應度函數(shù)是PSO算法的關鍵[16],關系到算法的收斂速度,甚至對最優(yōu)解的選取也有影響。使用的適應度函數(shù)如式(3)所示。
式中:K 是迭代次數(shù);R(λ)為本次迭代值;R(0)為期望值。
極限學習機(Extreme Learning Machine,ELM)是由Huang等人在2006年提出[17],是一種新型的單隱層前饋神經(jīng)網(wǎng)絡。極限學習機相比于神經(jīng)網(wǎng)絡有泛化能力強、分類精度高、運算速度快等優(yōu)點,由于運行過程無需迭代并且隱含層節(jié)點隨機產(chǎn)生,能夠避免結果容易陷入局部最優(yōu)的缺陷[18],被廣泛應用于人工智能以及模式識別等領域。
軸承故障診斷需要多層映射,即需要很多隱藏層節(jié)點,而傳統(tǒng)的極限學習機只有單層映射,只有一個隱藏層,不能滿足需求。因此,將核函數(shù)理論用于極限學習機構造核極限學習機KELM。每個樣本的隱藏層輸出可以看作是的線性映射,這種線性映射可以呈現(xiàn)ax+b或徑向基函數(shù) (Radial Basis Function,RBF)的形式。
因此,
根據(jù)核函數(shù)的理論,可以構造隱式映射來代替式(4)中的根據(jù)核函數(shù)的理論,可以構造隱式映射來代替式(4)中的內積,即構造核函數(shù)來替代HHT,具體步驟如下所示:
令 HHT(i,j)=K(xi,xj),則有
式中:c是正則化系數(shù)。
ELM解決方案的公式可以寫成
相較于傳統(tǒng)的極限學習機,改進后的核極限學習機具有更強大的函數(shù)逼近功能和次線性分類能力[19-20]。
與支持向量機(Support Vector Machine,SVM)類似,凡是滿足Mercer條件的函數(shù),都可以作為核極限學習機的核函數(shù),例如Poly核函數(shù)、RBF核函數(shù)、Morlet核函數(shù)等。為了簡化算法,將SVM內核應用于ELM。令母小波函數(shù)為h(x),其可伸縮性和平移因子分別為a和b,則小波基函數(shù)可表示為
根據(jù)張量積理論,多維小波函數(shù)可以寫成多個一維小波函數(shù)的張量積:
根據(jù)式(10),可以將核函數(shù)轉化為
小波核函數(shù)具有非線性映射強的特點,可以高精度逼近任意函數(shù)。但在實際復雜的情況下,單一的小波核函數(shù)構造的核極限學習機不能滿足應用需求,如數(shù)據(jù)異變或不規(guī)則,樣本量較大,樣品不均勻分布等。因此,可將多個不同尺度的小波核函數(shù)相疊加構成組合核函數(shù)。疊加形式如下所示。
式(13)中尺度因子an主要作用是將不通尺度的小波核函數(shù)疊加,以獲得擁有強大分類能力的組合核函數(shù)。多尺度小波核函數(shù)構造如下
式(14)可以縮寫為
多尺度小波核極限學習機比單一的小波核極限學習機在參數(shù)選擇方面范圍更大,削弱了參數(shù)的影響,提高了識別效果,增強了逼近能力,使用更方便。
采用基于離散粒子群優(yōu)化算法(DPSO)與核極限學習機(MKELM)的一種新型的DPSO-MKELM算法解決風機齒輪箱故障診斷,其流程如圖1所示。
圖1 PSO-MKELM模型流程
針對某實際風電場故障部件軸承正常、軸承故障、保持架損壞等3種狀態(tài),采集振動信號,每種狀態(tài)各采集60組數(shù)據(jù)。其中對每種狀態(tài)各取50組數(shù)據(jù)用于訓練PSO-MKELM模型,剩下30組用于測試誤差,并與文獻[21]中的高斯核函數(shù)的RPF-KELM以及傳統(tǒng)的ELM算法進行了對比,結果如圖2~4所示。圖中“○”代表真實故障,“*”代表預測故障,只有兩者重合時代表預測正確,不重合代表預測錯誤;故障類別數(shù)字1代表軸承正常,數(shù)字2代表軸承故障,數(shù)字3代表保持架損壞。
從圖2~4可以看出,所提出的PSO-MKELM在3個故障識模型表現(xiàn)最突出,正確率達到96.667%。而文獻[21]所提出的RBF-KELM模型識別正確率為93.333%,傳統(tǒng)的ELM算法識別正確率只有90%,因此PSO-MKELM模型能夠很好地完成軸承故障識別的任務。與此同時,多尺度核函數(shù)參數(shù)選擇比單核函數(shù)參數(shù)選擇更復雜,因此,多尺度小波核函數(shù)具有更好的泛化能力和更高的分類精度。
圖2 RPF-KELM模型預測結果
圖3 ELM模型預測結果
圖4 DPSO-MKELM模型預測結果
針對PSO算法的迭代過程中容易陷入局部最優(yōu)狀態(tài),提出了一種改進型的DPSO算法。將滿足Mercer條件的多尺度小波核函數(shù)作為極限學習機的核函數(shù),并將其與DPSO結合,構建了DPSO-MKELM模型。在風輪機齒輪軸承故障診斷中的測試表明,多尺度小波核極限學習機在故障分類方面具有較高的精度,性能優(yōu)于RPF-KELM和ELM,具有更高的應用價值。