航空發(fā)動機性能參數(shù)的混沌識別與預(yù)測

邸亞洲，高 峰，王小飛，曲建嶺

（海軍航空工程學(xué)院青島校區(qū) 山東 青島266041）

航空發(fā)動機性能參數(shù)的混沌識別與預(yù)測

邸亞洲，高 峰，王小飛，曲建嶺

（海軍航空工程學(xué)院青島校區(qū) 山東 青島266041）

性能參數(shù)監(jiān)控是航空發(fā)動機監(jiān)控的重要手段之一，而對性能參數(shù)進行預(yù)測可以提前掌握航空發(fā)動機在未來時刻的性能狀況，從而預(yù)防和排除故障。本文首先介紹了改進的加權(quán)一階局域混沌預(yù)測算法，然后對航空發(fā)動機性能參數(shù)（轉(zhuǎn)差率S）進行了混沌識別，最后采用改進的加權(quán)一階局域預(yù)測算法對航空發(fā)動機性能參數(shù)進行了混沌預(yù)測。實驗結(jié)果表明，改進的加權(quán)一階局域預(yù)測算法具有很好的學(xué)習(xí)能力和較高的預(yù)測精度，適用于航空發(fā)動機性能參數(shù)監(jiān)控。

性能參數(shù)；發(fā)動機監(jiān)控；混沌識別；混沌預(yù)測

由于內(nèi)在品質(zhì)因素和外部環(huán)境因素共同作用的影響，航空發(fā)動機的性能會隨飛行時間的累積逐步衰退，包括轉(zhuǎn)子氣流通道變形、慣性時間減小、渦輪落壓比下降和排氣溫度裕度降低等[1-3]。盡管零部件的性能衰退不一定會導(dǎo)致故障發(fā)生，但是其變化過程可以由特性參數(shù)表現(xiàn)出來，而且在不同發(fā)動機有所區(qū)別[4]。因此，通過分析和預(yù)測發(fā)動機特征參數(shù)來監(jiān)控其性能演化情況，對于預(yù)防和排除故障，進而保證飛行安全，具有重要意義[5-7]。

文中將混沌預(yù)測理論引入航空發(fā)動機性能參數(shù)監(jiān)控，主要介紹了改進的加權(quán)一階局域混沌預(yù)測算法，并提取了轉(zhuǎn)差率S這一航空發(fā)動機性能特征參數(shù)，對其進行了混沌識別和預(yù)測。

1 改進的加權(quán)一階局域預(yù)測法

局域預(yù)測法相對于其它混沌預(yù)測算法具有柔韌性好、擬合速度快等優(yōu)點，在工程應(yīng)用中非常廣泛[8]，其中加權(quán)一階局域預(yù)測是最常用的經(jīng)典方法，它以相點之間的歐氏距離來確定鄰近點對預(yù)測參考點的影響權(quán)值系數(shù)[9]，算法原理如下：

設(shè){xi，i=1，2，…，N}為給定的時間序列，選擇適合的延遲時間τ和嵌入維m，然后按照如下方式重構(gòu)相空間X（t）=（x（t），x（t+τ），…，x（t+（m-1）τ））∈Rm，其中有t=1，2，…，M，M=N-（m-1）τ為相點數(shù)量。

設(shè)Xk為參考相點，它與前i（i＜k-1）步相點之間的歐氏距離為di=‖Xk-Xi‖。尋找di的最小值dmin，并選擇距離最近的n（n一般取m+1）個鄰近點Xki，i= 1，2，…n。則Xki相對Xk的權(quán)重系數(shù)為：

其中，a為調(diào)整參數(shù)。為了預(yù)測未來相Xk+1，使用線性模型來擬合鄰近點和它們未來相

按下式用最小二乘法計算待定參數(shù)a和b：

然后可得f的近似映射f′（Xki）=a+bXki，而Xk+1通常由下式估算：

文獻[10]指出吸引子局部前r步相點之間的時間相關(guān)性和空間相關(guān)性在局域預(yù)測中具有重要作用，因而分別采用了向量夾角余弦度量和歐氏距離度量，來計算吸引子軌跡前r步相點在時間軸上的關(guān)聯(lián)度和空間關(guān)聯(lián)度，并將其融合起來計算兩個相點之間的綜合相關(guān)性程度，綜合相關(guān)度定義如下：

設(shè)Xk為參考點，Xn為其鄰近點且滿足|k-n|＞T，T為混沌時間序列的平均周期，其坐標(biāo)分別為：

其中 Xk-1、Xn-1分別為 Xk、Xn的前一個時刻相點，而Xk-1、Xn-1的一步演化向量被定義為：

由于一步演化向量的方向代表了單個相點在時間軸上的演化方向，故將在時間軸上Xk-1與Xn-1之間的相關(guān)性φ（k-1，n-1）定義為：

而Xk-1與 Xn-1之間的空間相關(guān)性被定義為 φ（k-1，n-1），它仍然由歐氏距離來度量：

Xk與Xn之間綜合相關(guān)性κ（k，n）通過結(jié)合其前r步相點之間的時間相關(guān)性和空間相關(guān)性來計算。設(shè)Xk-i、Xn-i（i=0，1，2，…，r）分別為Xk與Xn的前 r步相點，則它們之間的時間及空間相關(guān)性分別為：

當(dāng)i=0時，由于不存在一步演化向量，取φ（k，n）=1。顯然φ（k-i，n-i）越小，Xk-i與Xn-i的一步演化向量的夾角余弦ωi越大，它們在時間軸上的相關(guān)性越?。欢眨╧-i，n-i）越小，Xk-i與Xn-i的空間相關(guān)性越大。Xk與Xn之間綜合相關(guān)性κ（k，n）定義為：

然后以κkn取代Xk與Xn之間的歐氏距離‖Xk-Xn‖來確定Xn的權(quán)重系數(shù)Wn。

2 發(fā)動機性能特征參數(shù)

在雙軸渦輪噴氣式發(fā)動機中，高壓轉(zhuǎn)子轉(zhuǎn)速n2與低壓轉(zhuǎn)子轉(zhuǎn)速n1之比，稱為轉(zhuǎn)差率[11]，用符號S表示，顯然S=n2/n1。當(dāng)發(fā)動機進口氣流總溫恒定時，若緩收油門，則高/低壓轉(zhuǎn)速均下降，但由于高壓渦輪的功率減小程度要小于低壓渦輪功，故高壓轉(zhuǎn)速的減小程度必然要比低壓轉(zhuǎn)速小，結(jié)果導(dǎo)致轉(zhuǎn)差率增大[12]。如圖1中實線所示，此時發(fā)動機處于穩(wěn)態(tài)。若發(fā)動機處于過渡態(tài)：加速時，低壓轉(zhuǎn)子比高壓轉(zhuǎn)子加速慢，轉(zhuǎn)差率比平衡態(tài)的高，如圖中虛線2所示；減速時，高壓轉(zhuǎn)子轉(zhuǎn)速下降得比低壓轉(zhuǎn)子的快，轉(zhuǎn)差率低于平衡態(tài)值，如圖中虛線3所示[12]。

伴隨著飛行時間的增加，發(fā)動機渦輪效率和壓氣機增壓比下降等會不斷下降，進而導(dǎo)致相同飛行條件下的轉(zhuǎn)差率出現(xiàn)變化。因而結(jié)合歷史數(shù)據(jù)對同一臺發(fā)動機的轉(zhuǎn)差率進行監(jiān)控和預(yù)測研究，不僅能夠準(zhǔn)確掌握發(fā)動機轉(zhuǎn)子的性能趨勢變化情況，而且及時發(fā)現(xiàn)可能出現(xiàn)的故障，從而保證飛行安全。

在航空發(fā)動機的性能趨勢監(jiān)控研究中，通常提取穩(wěn)定巡航態(tài)下的性能參數(shù)進行對比研究[2]。因此，提取某飛行架次中處于穩(wěn)定巡航態(tài)下的轉(zhuǎn)差率S參數(shù)如圖2所示（為方便下文表述，記圖2中的序列為TS）。

圖1 轉(zhuǎn)差率與高壓轉(zhuǎn)子轉(zhuǎn)速相互關(guān)系

圖2 發(fā)動機穩(wěn)定巡航態(tài)下的轉(zhuǎn)差率S曲線

從上圖可以看出，穩(wěn)定巡航態(tài)下的轉(zhuǎn)差率S變化不大，總體變化幅度在0.1左右，說明發(fā)動機進口氣壓比較穩(wěn)定，渦輪風(fēng)扇工作狀況良好。如果轉(zhuǎn)差率S的變化幅度比較大，說明發(fā)動機轉(zhuǎn)子系統(tǒng)的流量和功率匹配出現(xiàn)問題，需要對渦輪風(fēng)扇進行檢修。

3 發(fā)動機性能參數(shù)的混沌識別與預(yù)測

航空發(fā)動機是一個極其復(fù)雜的機械系統(tǒng)，采用混沌理論研究其動力學(xué)行為是當(dāng)下的研究熱點，且發(fā)動機狀態(tài)參數(shù)具有混沌特征已經(jīng)得到了有效證實[13]。在此基礎(chǔ)上，本節(jié)采用Lyapunov指數(shù)法來判別發(fā)動機性能參數(shù)的混沌性質(zhì)，并采用改進的加權(quán)一階局域預(yù)測法對性能參數(shù)進行了預(yù)測。

3.1 轉(zhuǎn)差率S的混沌識別

Lyapunov指數(shù)法是混沌識別的常用方法之一。Lyapunov指數(shù)與重構(gòu)相空間中混沌吸引子軌跡收縮或擴展的性質(zhì)相關(guān)聯(lián)，用于表示平均迭代的指數(shù)分離中的指數(shù)。在Lyapunov指數(shù)譜中，最大的Lyapunov指數(shù)決定了吸引子軌道發(fā)散即覆蓋整個吸引子的快慢，也就是定量表示“蝴蝶效應(yīng)”強弱的量[14]。在正的Lyapunov指數(shù)方向上，吸引子軌道迅速分離，系統(tǒng)的長期演化行為對初始條件非常敏感，表明運動是混沌的[15]。因此最大Lyapunov指數(shù)為正，被視為判斷混沌性質(zhì)的重要條件。

首先采用互信息法分別計算提取的TS的重構(gòu)參數(shù)延遲時間τ，計算結(jié)果如圖3所示?；バ畔⒘侩S延遲時間變化的第一個局部極小值時刻即為最佳延遲時間，由圖3可以看出τ=4。

圖3 轉(zhuǎn)差率S延遲時間

在的基礎(chǔ)上，采用Grassberger-Procaccia法計算原始數(shù)據(jù)的嵌入維數(shù)m，臨界距離依次為r=e0.2i-5（i= 1，2，…20）m=2，3，…8，計算降噪前后對應(yīng)嵌入維數(shù)的關(guān)聯(lián)積分曲線，結(jié)果如圖4（a）所示。將圖4（a）標(biāo)度區(qū)域內(nèi)（即-3到-1）的關(guān)聯(lián)積分點擬合成直線，得到關(guān)聯(lián)維數(shù)隨嵌入維增加變化的曲線，如圖4（b）所示?？梢钥闯?，關(guān)聯(lián)維數(shù)在m=5時達到穩(wěn)定，此時TS的飽和關(guān)聯(lián)維數(shù)為2.312。

圖4 轉(zhuǎn)差率S嵌入維

最后在m=5，τ=4的基礎(chǔ)上，采用小數(shù)據(jù)量法計算TS的最大Lyapunov指數(shù)。圖5為重構(gòu)的混沌吸引子的所有相點及其最鄰近點在r步演化后的平均距離y（r）變化情況，圖5中擬合直線的斜率即為最大Lyapunov指數(shù)，經(jīng)計算得λ1=0.074，因此可以判定TS為混沌序列。

圖5 轉(zhuǎn)差率S最大Lyapunov指數(shù)

3.2 轉(zhuǎn)差率S的混沌預(yù)測

以TS的前501個樣本為訓(xùn)練樣本，以后100個樣本為測試樣本，然后采用第2節(jié)介紹的方法重構(gòu)混沌吸引子，并采用加權(quán)一階局域預(yù)測法和改進方法分別進行預(yù)測，預(yù)測結(jié)果如圖6所示。

圖6 轉(zhuǎn)差率S預(yù)測結(jié)果

由上圖可以看出，改進方法相對于加權(quán)一階局域預(yù)測法具有更高的預(yù)測精度 （計算得均方根誤差為0.108）和可靠預(yù)測步長，能夠有效跟蹤和預(yù)測轉(zhuǎn)差率S的演化趨勢，對于發(fā)動機性能趨勢長期監(jiān)控具有重要意義。

4 結(jié) 論

1）改進的加權(quán)一階局域混沌預(yù)測算法綜合考慮了相點的鄰近點演化形態(tài)對其未來趨勢預(yù)測的影響，相對于加權(quán)一階局域預(yù)測法具有更高的預(yù)測精度。

2）特征參數(shù)轉(zhuǎn)差率S與航空發(fā)動機轉(zhuǎn)子系統(tǒng)的性能衰退狀況有關(guān)，且該參數(shù)的演化過程具有混沌特性，可以采用混沌理論來對其進行監(jiān)測研究。

3）改進的加權(quán)一階局域預(yù)測算法在預(yù)測轉(zhuǎn)差率S時具有較高的精度，可以用來進行發(fā)動機性能參數(shù)的長期預(yù)測研究。

[1]Armstrong J B，Simon D L.Implementation of an integrated on-board aircraft diagnostic architecture [R].Reston，VA:AIAA，2011-5859，2011.

[2]Sanja G.Aircraft turbine engine control research at NASA glenn research center[R].NASA/TM-2013-217821，NASA，2013.

[3]趙連春，杰克，馬丁利.飛機發(fā)動機控制—設(shè)計、系統(tǒng)分析和健康監(jiān)視[M].北京:航空工業(yè)出版社，2011.

[4]Simo D L，Litt J S.A data filter for identifying steady-state operating points in engine flight data for condition monitoring applications[J].Journal of Engineering for Gas Turbines and Power，2011，133:1-8.

[5]Linda C.NASA aviation safety reporting system [R].NASA/ARC-E-DAA-TN2844，NASA，2014.

[6]Volphni A.Enhanced self turning on-board realtime model for aircraft engine performance health tracking[R].NASA/CR-2008-215272，Washington: NASA，2008.

[7]譚巍，李冬，樊照遠，等.基于模型信息熵的航空發(fā)動機性能評估和可靠性分析 [J].航空發(fā)動機，2011，37（5）:45-48.

[8]王桓.電力時間序列的混沌識別與短期預(yù)測 [D].長沙:湖南大學(xué)，2009.

[9]Sun L Y.Prediction of multivariable chaotic time series with local polynomial fitting [J].Computers and Mathematics with Application， 2010， 59: 737-744.

[10]QU Jian-ling，WANG Xiao-fei，QIAO Yu-chuan. An improved localweighted linearprediction model for chaotic times series[J].Chinese Physics. Letters，2014，31（1）:1-9.

[11]周勝明，曲建嶺，高峰.非線性支持向量回歸在航空發(fā)動機性能參數(shù)預(yù)測中的應(yīng)用[J].測控技術(shù)，2012，2:35-41.

[12]何立明.飛機推進系統(tǒng)原理[M].北京:國防工業(yè)出版社，2006.

[13]張科星.基于混沌時間序列理論的民航發(fā)動機性能參數(shù)預(yù)測模型研究[D].天津:中國民航大學(xué)，2009.

[14]Vaccari D A，Wang H K.Multivariate polynomial regression for identification of chaotic time series [J].Mathematical and Computer Modelling of Dynamical Systems，2007，13（4）:395-412.

[15]Samanta B.Prediction of chaotic time series using computational intelligence[J].Expert Systems with Application，2011（38）:11406-11411.

Chaotic recognition and forecasting of aeroengine performance parameter

DI Ya-zhou，GAO Feng，WANG Xiao-fei，QU Jian-ling
（Naval Aeronautical Engineering Institute Qingdao Branch，Qingdao 266041，China）

Performance parameter monitoring is a key means of aeroengine monitoring，and performance parameter forecasting can be used to obtain the future performance condition of aeroengine，thus preventing and eliminating faults.This article firstly introduces the improved local weighted linear chaotic forecast model briefly，then the aeroengine character parameter（i.e.rotor speed ratio S）is recognized as chaotic，at last the aeroengine character parameter is forecasted by using the improved local weighted linear forecast model.Experimental results show that the improved local weighted linear forecast model has good learning capability and high forecasting accuracy， which is suitable to aeroengine character parameters monitoring.

performance parameter；aeroengine monitoring；chaotic recognition；forecasting；chaotic forecasting

TN98

：A

：1674－6236（2017）03-0141-04

2016-05-29稿件編號：201605278

國家自然科學(xué)基金（51505491）

邸亞洲（1968—），男，河北石家莊人，碩士，副教授。研究方向：飛參數(shù)據(jù)研究及應(yīng)用。