旋轉(zhuǎn)機(jī)械運(yùn)行穩(wěn)定性劣化的高階累積量特征提取方法

蔣章雷　徐小力

北京信息科技大學(xué)現(xiàn)代測控技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，北京，100192

?

旋轉(zhuǎn)機(jī)械運(yùn)行穩(wěn)定性劣化的高階累積量特征提取方法

蔣章雷徐小力

北京信息科技大學(xué)現(xiàn)代測控技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，北京，100192

摘要：研究了基于高階累積量理論的1.5維譜、四階累積量對角切片譜方法在提取運(yùn)行穩(wěn)定性劣化特征方面的性能;從敏感性、趨勢性、差異性、一致性角度，討論了兩種特征提取方法的趨勢預(yù)測適用性?；跇?gòu)建的轉(zhuǎn)子實(shí)驗(yàn)臺進(jìn)行轉(zhuǎn)子系統(tǒng)多種劣化類型下不同劣化程度狀態(tài)的特征提取實(shí)驗(yàn),檢驗(yàn)兩種方法作為特征提取手段的性能,解決了風(fēng)電機(jī)組傳動系統(tǒng)運(yùn)行穩(wěn)定性劣化的狀態(tài)診斷、劣化趨勢預(yù)測中特征提取方法的選擇缺少理論依據(jù)的問題。

關(guān)鍵詞：風(fēng)電機(jī)組；穩(wěn)定性劣化；1.5維譜；四階累積量對角切片譜；特征提取

0引言

我國一直鼓勵、支持風(fēng)電等戰(zhàn)略性新興產(chǎn)業(yè)的發(fā)展，《能源發(fā)展“十二五”規(guī)劃》中提出大力發(fā)展風(fēng)能等可再生能源，加快風(fēng)能資源的分散開發(fā)利用[1]?！讹L(fēng)電發(fā)展“十二五”規(guī)劃》中提出:加快風(fēng)電產(chǎn)業(yè)技術(shù)升級，提高風(fēng)電的技術(shù)性能和產(chǎn)品質(zhì)量，使風(fēng)電成為具有較強(qiáng)國際競爭力的重要戰(zhàn)略性新興產(chǎn)業(yè)[2]。

隨著風(fēng)電機(jī)組的建設(shè)和運(yùn)行，風(fēng)電機(jī)組故障數(shù)量逐年增加，由此造成的經(jīng)濟(jì)損失越來越嚴(yán)重。風(fēng)電機(jī)組傳動系統(tǒng)包括傳動軸、齒輪箱等，是風(fēng)電機(jī)組故障多發(fā)部件[3]，從一系列公開發(fā)表的資料和文獻(xiàn)來看，其狀態(tài)診斷及趨勢預(yù)測研究日益成為國內(nèi)外學(xué)者研究的熱點(diǎn)，而特征提取方法是研究的難點(diǎn)[4-6]。高階累積量是解決非線性、非平穩(wěn)、非高斯、非最小相位、非因果信號的主要手段，利用高階累積量理論提取不同于傳統(tǒng)時頻分析方法的信號特征已成為機(jī)電設(shè)備狀態(tài)診斷研究中的一項(xiàng)重要內(nèi)容[7-8]。

本文將高階累積量理論應(yīng)用于運(yùn)行穩(wěn)定性劣化特征提取，考察高階累積量方法作為特征提取手段的性能，解決風(fēng)電機(jī)組傳動系統(tǒng)運(yùn)行穩(wěn)定性劣化的狀態(tài)診斷、劣化趨勢預(yù)測中特征提取方法的選擇缺少理論依據(jù)的問題。

1高階累積量的基本理論

1.1高階累積量的定義

高階累計(jì)量是指階數(shù)大于2的累計(jì)量,包括高階累積量、高階累積量譜等。

令單個隨機(jī)變量x具有概率密度f(x),定義其第一特征函數(shù)如下：

(1)

定義單個隨機(jī)變量x的累積量生成函數(shù)(第二特征函數(shù))如下：

ψ(ω)=lnΦ(ω)

(2)

定義單個隨機(jī)變量x的k階累積量為第二特征函數(shù)的k階導(dǎo)數(shù)在原點(diǎn)的值：

(3)

將單個隨機(jī)變量的高階累積量定義進(jìn)行推廣,得到隨機(jī)向量的高階矩與高階累積量的定義[9]:

令x=(x1,x2,…,xk)T是一隨機(jī)向量,且ω=(ω1,ω2,…,ωk),則隨機(jī)向量的第一特征函數(shù)為Φ(ω1,ω2,…,ωk)=E{ej(ω1x1+ω2x2+…+ωkωk)},對其求r=v1+v2+…+vk次偏導(dǎo)數(shù),則有

(4)

若令ω1=ω2…=ωk=0,則由式(4)得到隨機(jī)向量x的r階矩的定義：

(5)

參照隨機(jī)變量的第二特征函數(shù)定義,得到隨機(jī)向量的第二特征函數(shù)：

ψ(ω1,ω2,…,ωk)=lnΦ(ω1,ω2,…,ωk)

(6)

則隨機(jī)向量的r階累積量定義為

(7)

其中,v1+v2+…+vk=r,特別地，取v1=v2=…=vk=1,得到k階矩和k階累積量，記作:mk=m1,1…,1=mom(x1,x2,…,xk),ck=c1,1,…,1=cum(x1,x2,…,xk)。

設(shè)數(shù)列{x(n)}為零均值k階平穩(wěn)隨機(jī)過程，則其k階矩定義為

mkx(τ1,τ2,…,τk-1)=mom(x(n),x(n+τ1),

…,x(n+τn-1))

(8)

式中,τ1,τ2,…,τk-1為延時。

其k階累積量定義為

mkx(τ1,τ2,…,τk-1)=cum(x(n),

x(n+τ1),…,x(n+τn-1))

(9)

1.2高階累積量譜的定義

理論上,高斯過程大于二階的累積量譜為零,而非零累積量譜提供了一種非高斯性度量，因此工程上通常以高階累積量譜作為分析非高斯隨機(jī)過程的主要數(shù)學(xué)工具[9-10]。

若式(9)中高階累積量ckx(τ1,τ2,…,τk-1)是絕對可和的，即

(10)

則k階累積量譜定義為k階累積量的k-1維Fourier變換，即

(11)

根據(jù)式(9),三階累積量定義為

c3x(τ1,τ2)=cum(x(n),x(n+τ1)，x(n+τ2))

(12)

根據(jù)式(11)，三階累積量譜定義為

(13)

習(xí)慣上，高階累積量譜簡稱高階譜或多譜，我們稱三階譜為雙譜，稱四階譜為三譜。特別地,S2x(ω)稱為數(shù)列{x(n)}的功率譜。

2運(yùn)行穩(wěn)定性劣化的高階累積量特征提取方法

高階累積量及其高階累積量譜是多維函數(shù)，維數(shù)的增加將導(dǎo)致計(jì)算復(fù)雜性大大提高，在對實(shí)際信號作高階累積量估計(jì)時，多維的累積量計(jì)算起來非常復(fù)雜；同時對于三維以上的數(shù)據(jù)體難以用圖形進(jìn)行表示，因此必須降低高階累積量的維數(shù)。Mendel[11]提出利用高階累積量一維片段及其一維Fourier變換作為從高階累積量中提取有用信息的方法,即對高階累積量作切片處理，以降低高階累積量及其高階累積量譜的維數(shù),這主要包括1.5維譜以及四階累積量對角切片譜方法。

2.11.5維譜的定義及性質(zhì)分析

考慮一個非高斯平穩(wěn)隨機(jī)過程x(n),它的三階累積量為c3x(τ1,τ2)=cum(x(n),x(n+τ1),x(n+τ2)),定義三階累積量一維切片為

(14)

式(14)稱為三階累積量一維對角切片[12]。

定義三階累積量一維對角切片譜為

(15)

式(15)為三階累積量一維對角切片的一維Fourier變換,定義為三階累積量對角切片譜，或稱為1.5維譜。

1.5維譜方法僅使用了三階累積量的對角線切片，而且只進(jìn)行了一維Fourier變換，因此可以大大減小計(jì)算量。

2.2四階累積量對角切片譜的定義及性質(zhì)分析

一個非高斯平穩(wěn)過程x(n)的四階累積量的一維切片定義如下：

(16)

x(n)的四階累積量的對角切片C4x(τ,τ,τ)的Fourier變換,定義為x(n)的四階累積量對角切片譜：

(17)

3高階累積量特征提取方法的劣化特征提取性能實(shí)驗(yàn)

風(fēng)電機(jī)組傳動系統(tǒng)故障大部分是傳動系統(tǒng)的趨勢性故障，都會經(jīng)歷發(fā)生、發(fā)展直至惡化的階段，其狀態(tài)劣化過程對應(yīng)著各個劣化狀態(tài)之間的遷移過程。為了準(zhǔn)確描述劣化狀態(tài)之間的遷移過程，找到狀態(tài)劣化發(fā)展趨勢與劣化特征之間的對應(yīng)關(guān)系，需要模擬傳動系統(tǒng)的多種劣化類型，每一種劣化類型下需要包含不同劣化程度狀態(tài)。為進(jìn)行風(fēng)場風(fēng)電機(jī)組難以實(shí)施的各種劣化實(shí)驗(yàn)，基于典型傳動形式構(gòu)建轉(zhuǎn)子實(shí)驗(yàn)臺。

設(shè)置轉(zhuǎn)子系統(tǒng)多種劣化類型下不同劣化程度狀態(tài)，并采集振動數(shù)據(jù)。劣化類型包括：轉(zhuǎn)子不平衡狀態(tài)、基礎(chǔ)松動狀態(tài)、動靜碰摩狀態(tài)、轉(zhuǎn)子不平衡與動靜碰摩的耦合故障狀態(tài)、基礎(chǔ)松動與動靜碰摩的耦合故障狀態(tài)、轉(zhuǎn)子不平衡與基礎(chǔ)松動的耦合故障狀態(tài)等，每種劣化類型下包含不同的劣化程度。

3.1實(shí)驗(yàn)設(shè)備及實(shí)驗(yàn)數(shù)據(jù)的采集

實(shí)驗(yàn)數(shù)據(jù)采集現(xiàn)場如圖1所示，轉(zhuǎn)子轉(zhuǎn)速設(shè)定為1700 r/min,以近似模擬額定風(fēng)速下某型風(fēng)電機(jī)組齒輪箱輸出軸的轉(zhuǎn)速；數(shù)據(jù)采集系統(tǒng)選用丹麥B&K公司的PULSE系統(tǒng),采集儀型號為3050-B-060,為六通道數(shù)據(jù)采集儀;通過加速度振動傳感器采集振動信號,傳感器型號為INV9832A,為三向加速度振動傳感器;采樣頻率設(shè)定為10 240 Hz,每組振動數(shù)據(jù)包含2048個采樣點(diǎn)。

3.1.1轉(zhuǎn)子正常、轉(zhuǎn)子不平衡運(yùn)行狀態(tài)

轉(zhuǎn)子不平衡是風(fēng)電機(jī)組中較常見的一種故障狀態(tài)[13],故設(shè)置轉(zhuǎn)子多個不平衡運(yùn)行狀態(tài)進(jìn)行模擬實(shí)驗(yàn)。轉(zhuǎn)子上固定有旋轉(zhuǎn)圓盤，圓盤沿圓周均布著16個圓孔，通過在圓孔中擰入不同數(shù)量的配重螺釘實(shí)現(xiàn)轉(zhuǎn)子不平衡狀態(tài)下不同劣化程度的模擬：分別將兩個、四個、六個配重螺釘擰入旋轉(zhuǎn)圓盤，設(shè)定此時轉(zhuǎn)子狀態(tài)分別為輕度、中度、重度不平衡。給出每種狀態(tài)下各自旋入配重螺釘?shù)目傎|(zhì)量與旋轉(zhuǎn)圓盤質(zhì)量的百分比，以此作為不平衡劣化程度的度量標(biāo)準(zhǔn)，見表1。通過實(shí)驗(yàn)得到轉(zhuǎn)子正常運(yùn)行狀態(tài)、不平衡劣化類型下不同劣化程度狀態(tài)的振動數(shù)據(jù)。

3.1.2基礎(chǔ)松動狀態(tài)

松動故障會使機(jī)械系統(tǒng)出現(xiàn)偏心而產(chǎn)生不平衡力,當(dāng)不平衡力超過重力時,機(jī)器就會被不平衡力周期性地抬起,導(dǎo)致機(jī)械系統(tǒng)出現(xiàn)更嚴(yán)重的故障狀態(tài)[14]。松動轉(zhuǎn)子實(shí)驗(yàn)臺上軸承座支撐架與底座的緊固螺釘，以減小支撐結(jié)構(gòu)彈性剛度，模擬基礎(chǔ)松動狀態(tài)。將軸承座支撐架的一個緊固螺釘旋松一圈代表輕度基礎(chǔ)松動狀態(tài)；將緊固螺釘旋松兩圈代表中度基礎(chǔ)松動狀態(tài)。采集得到基礎(chǔ)松動狀態(tài)下不同劣化程度狀態(tài)的振動數(shù)據(jù)。

3.1.3動靜碰摩狀態(tài)

動靜碰摩除了引起轉(zhuǎn)子、定子之間的間隙增大,軸承支承磨損之外,更重要的是會引發(fā)轉(zhuǎn)軸不穩(wěn)定的反進(jìn)動,導(dǎo)致接觸面上產(chǎn)生很大的法向力和摩擦力,轉(zhuǎn)子的反向渦動不收斂，反向渦動轉(zhuǎn)速迅速增大,最終引發(fā)旋轉(zhuǎn)機(jī)械的嚴(yán)重?fù)p壞[14]。在轉(zhuǎn)子實(shí)驗(yàn)臺中利用一個頂桿來進(jìn)行動靜碰摩實(shí)驗(yàn)[15],將動靜件碰摩定性地分成四個階段：剛開始碰摩、早期尖銳型觸碰、中期半尖銳型碰摩和晚期平鈍型碰摩,動靜碰摩原理圖如圖2a所示。選擇在轉(zhuǎn)子實(shí)驗(yàn)臺上固定一個可以旋入碰摩塊的支撐架，碰摩塊與轉(zhuǎn)軸之間的間隙可通過螺紋調(diào)整，如圖2b所示，其中，最左側(cè)支撐架為可以旋入碰摩塊的支撐架。

設(shè)置動靜碰摩不同劣化程度的兩個狀態(tài),將剛開始碰摩到早期尖銳型觸碰作為輕度動靜碰摩狀態(tài),將早期尖銳型觸碰到中期半尖銳型碰摩作為中度動靜碰摩狀態(tài),采集轉(zhuǎn)子系統(tǒng)振動信號。

3.1.4耦合故障

在轉(zhuǎn)子系統(tǒng)實(shí)際運(yùn)行中,兩種或者兩種以上形式的故障有可能同時出現(xiàn)，即出現(xiàn)耦合故障，耦合故障的存在,會增加轉(zhuǎn)子系統(tǒng)故障診斷、趨勢預(yù)測的復(fù)雜性。

設(shè)置轉(zhuǎn)子輕度不平衡與輕度動靜碰摩狀態(tài)為耦合故障輕度狀態(tài)；中度不平衡與中度動靜碰摩狀態(tài)為耦合故障中度狀態(tài)。同理，設(shè)置轉(zhuǎn)子不平衡與基礎(chǔ)松動的耦合故障的輕度與中度狀態(tài)、動靜碰摩與基礎(chǔ)松動的耦合故障的輕度與中度狀態(tài)。采集每種狀態(tài)下的振動信號。

3.2特征提取方法的劣化特征提取性能分析

應(yīng)用1.5維譜、四階累積量對角切片譜方法對傳動系統(tǒng)多種劣化類型不同劣化程度數(shù)據(jù)進(jìn)行特征提取處理，顯示特征提取結(jié)果，得到其劣化特征譜圖。通過特征譜圖分析不同特征提取方法的劣化特征提取性能。

作為劣化特征提取性能的分析依據(jù)，本文提出如下判斷特征提取方法對于狀態(tài)劣化是否具有差異性、一致性、敏感性、趨勢性的標(biāo)準(zhǔn)：

(1)依據(jù)系統(tǒng)在正常運(yùn)行狀態(tài)下與不同劣化類型下的特征譜圖是否具有明顯區(qū)別，作為判斷特征提取方法是否具有差異性的標(biāo)準(zhǔn)。

(2)依據(jù)同一種劣化類型下不同劣化程度的特征譜圖中譜峰分布形態(tài)是否相似，作為判斷特征提取方法是否具有一致性的標(biāo)準(zhǔn)。

(3)依據(jù)特征譜圖中最大峰值在系統(tǒng)由正常運(yùn)行狀態(tài)遷移到輕度劣化狀態(tài)時變化是否顯著，作為判斷特征提取方法對狀態(tài)劣化是否具有敏感性的標(biāo)準(zhǔn)。

(4)依據(jù)特征譜圖中最大峰值是否隨著劣化程度的增大而不斷增大，能否表現(xiàn)出趨勢性信息,作為判斷特征提取方法對系統(tǒng)由正常運(yùn)行狀態(tài)到逐漸劣化的全過程是否具有趨勢性的標(biāo)準(zhǔn)。

顯然,若特征提取方法對于狀態(tài)劣化具有差異性、一致性，則這種方法適合用于劣化類型的診斷研究；若特征提取方法對于狀態(tài)劣化具有敏感性、趨勢性，則這種方法適合用于運(yùn)行穩(wěn)定性劣化趨勢的預(yù)測研究。

3.2.11.5維譜劣化特征提取性能分析

單一劣化類型的劣化特征提取。對轉(zhuǎn)子正常運(yùn)行狀態(tài)、不平衡運(yùn)行狀態(tài)的振動數(shù)據(jù)進(jìn)行基于1.5維譜的特征提取操作，得到轉(zhuǎn)子相關(guān)運(yùn)行狀態(tài)的特征譜圖,如圖3～圖10所示。

由圖3～圖10的特征譜圖可知：

(1)圖4所示轉(zhuǎn)子正常狀態(tài)特征譜圖中最大峰值位于224 Hz處，幅值為1942；圖6所示輕度不平衡狀態(tài)特征譜圖中最大峰值位于224 Hz處,幅值急劇增大至4519。最大峰值在系統(tǒng)由正常狀態(tài)劣化到輕度不平衡狀態(tài)時變化十分顯著,表明1.5維譜特征提取方法對狀態(tài)劣化具有敏感性,這有利于揭示系統(tǒng)早期劣化進(jìn)程。

(2)在轉(zhuǎn)子不平衡狀態(tài)特征譜圖的最大峰值明顯變化的同時，譜峰的分布形態(tài)并未出現(xiàn)明顯變化，表明同一種劣化類型下的劣化特征具有一致性，這為風(fēng)電機(jī)組傳動系統(tǒng)的狀態(tài)診斷研究提供了良好基礎(chǔ)。

(3)隨著轉(zhuǎn)子系統(tǒng)不平衡狀態(tài)的劣化,最大峰值經(jīng)歷了4519→5071→5505的變化,最大峰值均出現(xiàn)在224 Hz附近。轉(zhuǎn)子系統(tǒng)同一劣化類型下不同劣化程度的最大峰值存在明顯的趨勢性。

(4)轉(zhuǎn)子系統(tǒng)正常狀態(tài)下頻率成分比較單一,劣化到不平衡狀態(tài)后,譜圖中低頻部分與高頻部分的峰值數(shù)量明顯增多,出現(xiàn)了多譜峰現(xiàn)象。正常運(yùn)行狀態(tài)下的特征譜圖與不平衡運(yùn)行狀態(tài)下的特征譜圖具有明顯的差異性。

對基礎(chǔ)松動狀態(tài)數(shù)據(jù)、動靜碰摩狀態(tài)數(shù)據(jù)進(jìn)行特征提取,結(jié)果表明,系統(tǒng)出現(xiàn)輕度劣化后，特征譜圖中峰值數(shù)值變化明顯，特征提取方法對系統(tǒng)由正常運(yùn)行狀態(tài)遷移到輕度劣化狀態(tài)具有敏感性；同種劣化類型下不同劣化程度的特征譜圖具有一致性；系統(tǒng)由正常運(yùn)行狀態(tài)逐漸劣化的過程能夠表現(xiàn)出趨勢性。同時，轉(zhuǎn)子正常運(yùn)行狀態(tài)特征譜圖中譜峰分布形態(tài)與不同劣化類型下特征譜圖中譜峰分布形態(tài)具有比較明顯的差異性。

耦合劣化類型的劣化特征提取。利用1.5維譜特征提取方法對轉(zhuǎn)子輕度不平衡耦合輕度基礎(chǔ)松動狀態(tài)、轉(zhuǎn)子中度不平衡耦合中度基礎(chǔ)松動狀態(tài)下的數(shù)據(jù)進(jìn)行特征提取操作，得到圖11、圖12。

由轉(zhuǎn)子系統(tǒng)正常運(yùn)行狀態(tài)下的特征譜圖及圖11、圖12可看出：

(1)隨著轉(zhuǎn)子系統(tǒng)由正常運(yùn)行狀態(tài)劣化到耦合故障輕度狀態(tài),特征譜圖中最大峰值由1942(位于224 Hz處)減小到1762(位于208 Hz處),這表明該特征提取方法對于系統(tǒng)耦合劣化狀態(tài)的敏感性低于單一劣化狀態(tài)的敏感性。

(2)特征譜圖中譜峰的分布形態(tài)不盡相同：轉(zhuǎn)子系統(tǒng)正常狀態(tài)的特征譜圖只有一個明顯波峰;耦合故障中度劣化譜圖出現(xiàn)了多個波峰，表明對于耦合故障類型，該方法提取的劣化特征的一致性低于單一劣化類型的一致性。

(3)隨著轉(zhuǎn)子系統(tǒng)逐漸劣化，特征譜圖中最大峰值隨著劣化程度的增加呈現(xiàn)不斷減小的趨勢，這表明，耦合故障狀態(tài)下不同劣化程度的最大峰值的趨勢性與劣化程度的趨勢性相反。

分析不平衡耦合動靜碰摩狀態(tài)、動靜碰摩耦合基礎(chǔ)松動狀態(tài)的特征譜圖，結(jié)果表明：當(dāng)耦合故障狀態(tài)中某個單一劣化類型占主要成分時,該特征提取方法對于狀態(tài)劣化的敏感性、一致性有所提高。

3.2.2四階累積量對角切片譜劣化特征提取性能分析

三階累積量是三階中心矩,描述了概率分布的非對稱性；四階累積量是任意隨機(jī)變量的矩與正態(tài)隨機(jī)變量的同階矩的差，描述了概率分布相比正態(tài)分布平坦或者尖銳的程度[16]?？紤]到四階累積量能夠揭示不同于三階累積量的信息，因此繼續(xù)考察四階累積量對角切片譜特征提取方法的性能。

(1)單一劣化類型的劣化特征提取。對轉(zhuǎn)子正常、不平衡狀態(tài)振動數(shù)據(jù)進(jìn)行四階累積量對角切片譜操作，得到各狀態(tài)的譜圖，如圖13～圖16所示。

由圖13～圖16可知：譜圖中最大峰值在系統(tǒng)由正常運(yùn)行狀態(tài)劣化到輕度不平衡狀態(tài)時變化顯著，表明四階累積量對角切片譜特征提取方法對于系統(tǒng)狀態(tài)劣化具有敏感性；不同劣化程度的譜圖中譜線的分布形態(tài)相似，且譜圖表現(xiàn)出的頻率成份較為單一，表明同一種劣化類型下的劣化特征具有一致性；譜圖中最大峰值變化為：435→939→1575.3→705，可見隨著不平衡狀態(tài)的逐漸劣化，最大峰值沒有表現(xiàn)出確定的趨勢性信息。

對基礎(chǔ)松動狀態(tài)、動靜碰摩狀態(tài)的特征提取的研究結(jié)果表明，該特征提取方法對系統(tǒng)正常運(yùn)行狀態(tài)、不同劣化類型狀態(tài)能夠表現(xiàn)出一定的差異性；對系統(tǒng)由正常運(yùn)行狀態(tài)遷移到輕度劣化狀態(tài)具有敏感性；同種劣化類型下不同劣化程度的譜圖具有一致性；對系統(tǒng)由正常運(yùn)行狀態(tài)逐漸劣化的全過程沒有表現(xiàn)出趨勢性。

(2)耦合劣化類型的劣化特征提取。對輕度不平衡耦合輕度基礎(chǔ)松動狀態(tài)、中度不平衡耦合中度基礎(chǔ)松動狀態(tài)的振動數(shù)據(jù)進(jìn)行四階累積量對角切片譜操作，得到兩種狀態(tài)的四階累積量對角切片譜圖，如圖17、圖18所示。

由轉(zhuǎn)子正常運(yùn)行狀態(tài)的譜圖及圖17、圖18可知：四階累積量對角切片譜方法對于系統(tǒng)狀態(tài)劣化具有一定的敏感性；同種劣化類型下的特征譜圖分布形態(tài)具有一致性；不具有表征狀態(tài)劣化的趨勢性；具有比較明顯的差異性。

3.3特征提取方法的趨勢預(yù)測適用性討論

針對多種劣化類型下不同劣化程度狀態(tài)數(shù)據(jù)，基于1.5維譜和四階累積量對角切片譜方法進(jìn)行了特征提取分析，得到了兩種特征提取方法對狀態(tài)劣化的差異性、一致性、敏感性、趨勢性信息，見表2，其中方法1、2分別代表1.5維譜、四階累積量對角切片譜。

從兩種特征提取方法對于狀態(tài)劣化的敏感性、趨勢性、差異性、一致性的角度來討論這兩種特征提取方法作為穩(wěn)定性劣化趨勢預(yù)測的適用性[17-18]。

(1)對于單一劣化類型，1.5維譜特征提取方法與四階累積對角切片譜特征提取方法都具有良好的差異性、一致性,因此可以考慮融合兩種特征提取方法進(jìn)行單一劣化類型狀態(tài)的診斷研究;1.5維譜方法在敏感性、趨勢性方面的性能要優(yōu)于四階累積對角切片譜方法，因此可以利用1.5維譜方法進(jìn)行穩(wěn)定性劣化趨勢預(yù)測。

(2)對于耦合劣化類型,四階累積對角切片譜特征提取方法在差異性、一致性方面的性能要優(yōu)于1.5維譜特征提取方法,因此可以利用四階累積對角切片譜方法進(jìn)行耦合劣化類型狀態(tài)的診斷研究;1.5維譜與四階累積對角切片譜方法都不具有敏感性、趨勢性，因此這兩種方法都不適合用于耦合劣化類型穩(wěn)定性劣化趨勢預(yù)測。

4結(jié)語

本文針對風(fēng)電機(jī)組變工況、非平穩(wěn)運(yùn)行特點(diǎn)，研究了1.5維譜以及四階累積量對角切片譜方法用于特征提取的性能，基于敏感性、趨勢性、差異性、一致性衡量了特征提取方法的趨勢預(yù)測適用性，為風(fēng)電機(jī)組傳動系統(tǒng)運(yùn)行穩(wěn)定性劣化的狀態(tài)診斷、劣化趨勢預(yù)測中特征提取方法的選擇提供了理論依據(jù)。

參考文獻(xiàn):

[1]能源發(fā)展“十二五”規(guī)劃[R].北京: 國家能源局, 2013.

[2]風(fēng)電發(fā)展“十二五”規(guī)劃[R].北京:國家能源局,2012.

[3]Johan R,Margareta B L.Survey of Failures in Wind Power Systems with Focus on Swedish Wind Power Plants during 1997-2005[J].IEEE Transactions on Energy Conversion,2007,22(1):167-173.

[4]王國彪,何正嘉,陳雪峰，等.機(jī)械故障診斷基礎(chǔ)研究“何去何從”[J].機(jī)械工程學(xué)報,2013,49(1): 63-72.

Wang Guobiao, He Zhengjia, Chen Xuefeng,et al. Basic Research on Machinery Fault Diagnosis-What Is the Prescription [J]. Journal of Mechanical Engineering, 2013,49(1): 63-72.

[5]于傳洋,許寶杰,吳國新,等.一種風(fēng)電機(jī)組的故障預(yù)測算法[J].北京信息科技大學(xué)學(xué)報,2013,28(4)：35-38.

Yu Chuanyang, Xu Baojie, Wu Guoxin,et al. A Fault Prediction Algorithm of Wind Turbine Generator[J].Journal of Beijing Information Science and Technology University, 2013,28(4): 35-38.

[6]王紅軍,左云波.基于局部線性降維拓?fù)淇臻g的主軸故障診斷方法[J].北京信息科技大學(xué)學(xué)報, 2014,29(2):35-38.

Wang Hongjun, ZuoYunbo. Spindle Fault Diagnosis Method Based on Logal Linear Dimension Reduction Topological Space[J]. Journal of Beijing Information Science and Technology University, 2013,28(4):35-38.

[7]Montero F E H,Medina O C.The Application of Bispectrum on Diagnosis of Rolling Element Bearings:A Theoretical Approach[J].Mechanical Systems and Signal Processing,2008,22(3):588-596.

[8]Zhou Y,Chen J,Dong G M, et al. Application of the Horizontal Slice of Cyclicbispectrum in Rolling Element Bearings Diagnosis[J].Mechanical Systems and Signal Processing,2012,26(1):229-243.

[9]張賢達(dá).時間序列分析——高階統(tǒng)計(jì)量方法[M].北京:清華大學(xué)出版社,1996.

[10]邱天爽,張旭秀,李小兵,等.統(tǒng)計(jì)信號處理——非高斯信號處理及其應(yīng)用[M].北京:電子工業(yè)出版社,2004.

[11]Mendel J M.Tutorial on Higher-order Statistics(Spectra) in Signal Processing and System Theory Theoretical Results and Some Application[J].Proceedings of the IEEE,1991，79(3):278-305.

[12]Nikias C L,Mendel J M.Signal Processing with Higher-order Spectra[J].IEEE Signal Processing Magazine,1993,10(3):10-37.

[13]Villa L F,Renones A,Peran J R,et al.Statistical Fault Diagnosis Based on Vibration Analysis for Gear Test-bench under Non-stationary Conditions of Speed and Load[J].Mechanical Systems and Signal Processing,2012,29(5):436-446.

[14]聞邦椿,武新華,丁千,等.故障旋轉(zhuǎn)機(jī)械非線性動力學(xué)的理論與試驗(yàn)[M].北京:科學(xué)出版社,2004.

[15]胡蔦慶,張雨,劉耀宗,等.轉(zhuǎn)子系統(tǒng)動靜件間尖銳碰摩時的振動特征試驗(yàn)研究[J].中國機(jī)械工程,2002,13(9): 777-780.

Hu Yuanqing,Zhang Yu,Liu Yaozong,et al.Experiment Research on Vibration Characteristics of Sharp Rup-impact between Rotor and Stator of a Rotor System[J]. China Mechanical Engineering,2002,13(9): 777-780.

[16]王書明,朱培民,李宏偉,等.地球物理學(xué)中的高階統(tǒng)計(jì)量方法[M].北京:科學(xué)出版社,2006.

[17]左云波,王西彬,徐小力.DFAR模型在旋轉(zhuǎn)機(jī)械故障預(yù)測中的應(yīng)用[J].中國機(jī)械工程,2009,20(12):1460-1463.

Zuo Yunbo,Wang Xibin,Xu Xiaoli.Application of DFAR Model to the Trend Prediction of Fault in Rotary Machines[J].China Mechanical Engineering,2009,20(12):1460-1463.

[18]蔣章雷,徐小力.旋轉(zhuǎn)機(jī)械運(yùn)行穩(wěn)定性劣化的1.5維譜特征提取方法[J].中國機(jī)械工程,2015,26(23):3208-3213.

Jiang Zhanglei,Xu Xiaoli.Feature Extraction Method of 1.5-dimensional Faced to Running Stability Deterioration of Rotating Machinery[J].China Mechanical Engineering,2015,26(23):3208-3213.

(編輯陳勇)

Feature Extraction Method Based on High Order Cumulant Faced to Running Stability Deterioration of Rotating Machinery

Jiang ZhangleiXu Xiaoli

Key Laboratory of Modern Measurement & Control Technology,Ministry of Education,Beijing Information Science and Technology University,Beijing,100192

Key words:wind turbine;stability deterioration;1.5-dimension spectrum;fourth-order cumulant diagonal slice spectrum;feature extraction

Abstract:This paper focused on feature extraction method based on high order cumulant like 1.5-dimension spectrum and fourth-order cumulant diagonal slice spectrum,and also its performance on running stability deterioration.Experimental data of varying degrees of deterioration under various types of deterioration were collected,and sensitivity property,trend property,difference property,and consistency property were proposed to evaluate different methods.It can provide theoretical basis for selecting feature extraction method for fault diagnosis and trend prediction of running stability deterioration.

收稿日期：2015-04-09

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(51275052);北京市自然科學(xué)基金資助重點(diǎn)項(xiàng)目(3131002);北京市教委科研計(jì)劃資助重點(diǎn)項(xiàng)目(KZ201311232036);現(xiàn)代測控技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室開放課題資助項(xiàng)目(KF20141123202)

作者簡介：蔣章雷,男,1983年生。北京信息科技大學(xué)現(xiàn)代測控技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室助理研究員、博士。主要研究方向?yàn)闄C(jī)電裝備運(yùn)行狀態(tài)監(jiān)控檢測技術(shù)。發(fā)表論文10余篇。徐小力，男,1951年生。北京信息科技大學(xué)現(xiàn)代測控技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室教授、博士研究生導(dǎo)師。

中圖分類號：TH17

DOI：10.3969/j.issn.1004-132X.2016.02.004