基于MIC特征提取與BO-CatBoost的航空發(fā)動(dòng)機(jī)RUL預(yù)測(cè)

李東君, 李 亞, 李東文, 朱貴富

(1.昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院,昆明,650504;2.昆明理工大學(xué)信息化建設(shè)管理中心,昆明,650504)

航空發(fā)動(dòng)機(jī)作為飛機(jī)的核心部件之一,內(nèi)部結(jié)構(gòu)復(fù)雜,在運(yùn)行期間,其健康狀況易受環(huán)境及內(nèi)部構(gòu)造等多種因素的影響。相關(guān)研究表明,在飛機(jī)因機(jī)械故障引起的飛行事故中,發(fā)動(dòng)機(jī)故障導(dǎo)致的飛行事故數(shù)量占比最高[1]。飛機(jī)發(fā)動(dòng)機(jī)一旦發(fā)生故障或失效,將會(huì)嚴(yán)重威脅人們的生命健康及財(cái)產(chǎn)安全。為提升發(fā)動(dòng)機(jī)運(yùn)行的安全性、可靠性及穩(wěn)定性,并有效保障人們的生命及財(cái)產(chǎn)安全,對(duì)發(fā)動(dòng)機(jī)開展視情維修與健康監(jiān)測(cè)具有重要意義。

故障預(yù)測(cè)和健康管理(prognostics and health management,PHM)技術(shù)主要利用機(jī)械設(shè)備在運(yùn)行期間產(chǎn)生的各類數(shù)據(jù),運(yùn)用數(shù)據(jù)處理與數(shù)據(jù)分析等手段,實(shí)現(xiàn)對(duì)復(fù)雜設(shè)備的健康狀態(tài)檢測(cè)、預(yù)測(cè)管理及維修策略支持,有效降低故障發(fā)生率和維修成本[2]。航空發(fā)動(dòng)機(jī)剩余使用壽命(remaining useful life,RUL)維護(hù)預(yù)測(cè)是PHM的關(guān)鍵任務(wù)之一。RUL預(yù)測(cè)的主要目標(biāo)是根據(jù)機(jī)械設(shè)備健康狀態(tài)的有效信息,對(duì)設(shè)備到達(dá)安全操作限制的剩余壽命時(shí)間進(jìn)行估計(jì)[3]。由于航空發(fā)動(dòng)機(jī)長(zhǎng)期工作在條件惡劣的環(huán)境中,為保證發(fā)動(dòng)機(jī)運(yùn)行的安全性,需要定期對(duì)發(fā)動(dòng)機(jī)進(jìn)行維護(hù),如何有效利用發(fā)動(dòng)機(jī)的歷史退化信息對(duì)發(fā)動(dòng)機(jī)的健康狀態(tài)進(jìn)行RUL預(yù)測(cè)是當(dāng)下研究的重點(diǎn)。

目前,因預(yù)測(cè)原理的差異,國(guó)內(nèi)外關(guān)于航空發(fā)動(dòng)機(jī)的RUL預(yù)測(cè)方法主要分為2類:基于物理失效模型的RUL預(yù)測(cè)方法和基于數(shù)據(jù)驅(qū)動(dòng)的RUL預(yù)測(cè)方法[4]。雖然基于物理失效模型的RUL預(yù)測(cè)方法具有良好的預(yù)測(cè)精度,但由于航空發(fā)動(dòng)機(jī)內(nèi)部系統(tǒng)構(gòu)造極為復(fù)雜,在對(duì)部件進(jìn)行退化過(guò)程建模時(shí)需要掌握部件原理的專業(yè)知識(shí),使得預(yù)測(cè)難度增大且普適性較差?；跀?shù)據(jù)驅(qū)動(dòng)的RUL預(yù)測(cè)方法無(wú)需建立復(fù)雜的數(shù)學(xué)模型,它借助機(jī)器學(xué)習(xí)手段挖掘和分析退化數(shù)據(jù)隱藏的內(nèi)部信息來(lái)獲得數(shù)據(jù)間內(nèi)在的關(guān)聯(lián)特征,且預(yù)測(cè)的性能和準(zhǔn)確率也較高[5-6]。

隨著人工智能技術(shù)及智能算法的飛速發(fā)展,基于數(shù)據(jù)驅(qū)動(dòng)的RUL預(yù)測(cè)方法逐漸成為國(guó)內(nèi)外學(xué)者研究的熱點(diǎn)。文獻(xiàn)[7]對(duì)歷史監(jiān)測(cè)數(shù)據(jù)進(jìn)行融合,將貝葉斯理論(bayesian)應(yīng)用到航空發(fā)動(dòng)機(jī)RUL預(yù)測(cè)中。文獻(xiàn)[8]利用最大信息系數(shù)(maximal information coefficient,MIC)對(duì)短期負(fù)荷數(shù)據(jù)進(jìn)行特征選擇,將篩選的特征輸入構(gòu)建的模型中進(jìn)行預(yù)測(cè)。文獻(xiàn)[9]提出了一種利用皮爾遜相關(guān)系數(shù)方法對(duì)發(fā)動(dòng)機(jī)數(shù)據(jù)子集進(jìn)行協(xié)變量篩選獲得最優(yōu)協(xié)變量表達(dá)式。文獻(xiàn)[10]采用類別特征梯度提升算法(categorical boosting,CatBoost)構(gòu)建孔隙壓力預(yù)測(cè)模型,有效提升了預(yù)測(cè)精度。文獻(xiàn)[11]提出了一種在聯(lián)邦學(xué)習(xí)框架中使用全局健康退化表示(a global health degradation representation ,GHDR)的新的RUL預(yù)測(cè)方法。

對(duì)于復(fù)雜機(jī)械設(shè)備的RUL預(yù)測(cè)問(wèn)題,雖然現(xiàn)有文獻(xiàn)所提方法在航空發(fā)動(dòng)機(jī)剩余壽命方面已經(jīng)取的一定的成績(jī),但綜合考慮發(fā)動(dòng)機(jī)自身復(fù)雜的物理結(jié)構(gòu)及運(yùn)行環(huán)境的惡劣條件等因素,對(duì)于發(fā)動(dòng)機(jī)的RUL預(yù)測(cè)還存在一些待解決的:發(fā)動(dòng)機(jī)傳感器由于受外部環(huán)境因素的影響,使監(jiān)測(cè)的數(shù)據(jù)受到噪聲干擾,無(wú)法有效去除環(huán)境噪聲;針對(duì)航空發(fā)動(dòng)機(jī)監(jiān)測(cè)參數(shù)非線性特點(diǎn),現(xiàn)有方法無(wú)法充分提取對(duì)發(fā)動(dòng)機(jī)運(yùn)行周期影響較大的關(guān)鍵特征,在RUL預(yù)測(cè)準(zhǔn)確度上仍有提升空間。

針對(duì)以上問(wèn)題,本文提出了一種基于MIC特征提取與貝葉斯優(yōu)化類別特征梯度提升(bayesian optimization categorical boosting,BO CatBoost)相結(jié)合的發(fā)動(dòng)機(jī)RUL預(yù)測(cè)方法?？紤]到不同監(jiān)測(cè)特征對(duì)發(fā)動(dòng)機(jī)壽命影響的差異,引入MIC算法計(jì)算各個(gè)特征對(duì)發(fā)動(dòng)機(jī)RUL的相關(guān)性,篩選出相關(guān)性較大的特征作為CatBoost預(yù)測(cè)模型的輸入,引入貝葉斯優(yōu)化(bayesian optimization,BO)算法對(duì)CatBoost預(yù)測(cè)模型中的超參數(shù)進(jìn)行訓(xùn)練和調(diào)優(yōu),得到RUL預(yù)測(cè)結(jié)果的最優(yōu)值。利用評(píng)價(jià)指標(biāo)均方根誤差(root mean square error,RMSE)、判定系數(shù)(coefficient of determination,R2)和平均絕對(duì)誤差(mean absolute error,MAE)對(duì)預(yù)測(cè)模型進(jìn)行性能分析和評(píng)估,有效驗(yàn)證了本文所提方法的可行性。

1 RUL預(yù)測(cè)方法

1.1 最大信息系數(shù)

最大信息系數(shù)MIC以互信息(mutual information,MI)為基礎(chǔ),是一種用于衡量2個(gè)變量之間線性或非線性相關(guān)性強(qiáng)弱的算法[12]。設(shè)X={x1,x2,…,xn}與Y={y1,y2,…,yn}分別為數(shù)據(jù)集中的隨機(jī)變量,n為樣本數(shù)量,則X與Y之間的MI為:

(1)

式中:p(x,y)為X與Y之間的聯(lián)合概率密度;p(x)與p(y)分別為X與Y之間的邊緣概率密度。

MIC克服了MI在計(jì)算連續(xù)變量的聯(lián)合概率密度函數(shù)困難的缺陷,能最大程度地找到兩變量之間的相關(guān)性[13]。MIC的計(jì)算公式為:

式中:B為樣本數(shù)量;N為樣本變量;I(x;y)為x與y之間的MI。2個(gè)變量間的MIC值越接近1,則其相關(guān)性越強(qiáng),MIC∈[0,1]。

1.2 CatBoost算法

類別特征梯度提升算法(CatBoost)是在梯度提升決策樹(gradient boosting decision tree,GBDT)算法的基礎(chǔ)上改進(jìn)的算法[14]。CatBoost算法是由類別特征和梯度提升組成的一種高準(zhǔn)確性梯度提升框架。該算法以對(duì)稱決策樹作為基學(xué)習(xí)器,解決了梯度偏差及預(yù)測(cè)偏移問(wèn)題,有效防止模型的過(guò)度擬合。該算法將數(shù)據(jù)的類別特征進(jìn)行編碼,每層分裂時(shí),設(shè)置分裂閾值,并將所有的類別特征與指定的特征進(jìn)行組合,參與下一層分裂。分裂結(jié)束后,使用梯度無(wú)偏來(lái)估計(jì)預(yù)測(cè)偏移,尋找最優(yōu)目標(biāo)。

建立分類特征樹的過(guò)程中,因需要考慮監(jiān)測(cè)數(shù)據(jù)之間的相關(guān)性,本文引入MIC算法來(lái)計(jì)算監(jiān)測(cè)參數(shù)之間的非線性關(guān)聯(lián)性。MIC能有效避免互信息在計(jì)算連續(xù)變量的聯(lián)合概率密度函數(shù)困難的問(wèn)題,且最大程度地挖掘2個(gè)變量之間的非線性相關(guān)性。

本文使用CatBoost回歸算法來(lái)解決發(fā)動(dòng)機(jī)RUL的預(yù)測(cè)問(wèn)題。設(shè)原始數(shù)據(jù)集為|D|={(x1,y1),(x2,y2),…,(xn,yn)},則σ=(σ1,σ2,…,σn)為|D|經(jīng)過(guò)重新排序后的序列狀態(tài)。

(3)

式中:xσj,k為數(shù)據(jù)集σj的第k個(gè)特征;p為先驗(yàn)概率,用來(lái)減少噪聲數(shù)據(jù)的干擾;?為大于0時(shí)的權(quán)重系數(shù)值,用于調(diào)節(jié)p的影響程度。

1.3 貝葉斯優(yōu)化算法

貝葉斯優(yōu)化(BO)算法是一種基于搜索函數(shù)和高斯過(guò)程(gaussian processes,GP)的參數(shù)更新優(yōu)化算法,根據(jù)定義的目標(biāo)函數(shù)迭代評(píng)估新的參數(shù)[15]。利用高斯過(guò)程進(jìn)行調(diào)參,設(shè)置最優(yōu)參數(shù)后,不斷迭代更新先驗(yàn)值和模型參數(shù),直到找到最優(yōu)的超參數(shù)組合。本文使用貝葉斯算法中的GP及采集函數(shù)(acquisition function,AF)優(yōu)化CatBoost中的超參數(shù)。GP回歸表達(dá)式:

f(x)～GP(m(x),K(x,x))

(4)

式中:K(x,x)為協(xié)方差矩陣;m(x)為均值向量函數(shù)。

采集函數(shù)對(duì)樣本的候選值進(jìn)行評(píng)估后得到最優(yōu)解。其計(jì)算公式為:

式中:σ(x)為GP的方差;u(x)為樣本均值。

φ(*)為標(biāo)準(zhǔn)正態(tài)分布的累積分布函數(shù),α為超參數(shù)。貝葉斯優(yōu)化CatBoost算法超參數(shù)流程如圖1所示。

圖1 貝葉斯優(yōu)化CatBoost算法超參數(shù)流程圖

1.4 評(píng)價(jià)指標(biāo)

為更好地驗(yàn)證本文所提RUL預(yù)測(cè)模型的準(zhǔn)確性,選取實(shí)驗(yàn)中2個(gè)常用的模型性能評(píng)價(jià)指標(biāo),即判定系數(shù)R2[16]、均方根誤差RMSE[17-19]與平均絕對(duì)誤差MAE[20]來(lái)衡量模型的預(yù)測(cè)性能。

RMSE通常用于度量模型的預(yù)測(cè)值和真實(shí)值之間的總體偏差,其值越小,則預(yù)測(cè)性能精度越高。RMSE的計(jì)算式為:

(6)

R2用于反映由樣本回歸線做出解釋的離差平方和中的比重,其值越接近1,說(shuō)明模型擬合度越好,其表達(dá)式如下:

(7)

MAE表示絕對(duì)誤差的平均值,其計(jì)算式為:

(8)

2 基于MIC與BO-CatBoost的預(yù)測(cè)模型框架

本文所提出的航空發(fā)動(dòng)機(jī)RUL預(yù)測(cè)模型框架總體介紹:利用MIC算法篩選出能夠表征發(fā)動(dòng)機(jī)退化性能的傳感器監(jiān)測(cè)參數(shù),為降低監(jiān)測(cè)參數(shù)量綱對(duì)預(yù)測(cè)模型的影響,將監(jiān)測(cè)參數(shù)歸一化在[0,1]之間。構(gòu)建基于貝葉斯超參數(shù)尋優(yōu)的CatBoost預(yù)測(cè)模型進(jìn)行訓(xùn)練,將測(cè)試集輸入訓(xùn)練結(jié)束的預(yù)測(cè)模型中預(yù)測(cè)發(fā)動(dòng)機(jī)的RUL。并通過(guò)評(píng)價(jià)指標(biāo)RMSE、R2、MAE來(lái)評(píng)價(jià)預(yù)測(cè)模型的性能。

1)數(shù)據(jù)預(yù)處理

在數(shù)據(jù)預(yù)處理階段,利用MIC算法分析各個(gè)監(jiān)測(cè)參數(shù)與發(fā)動(dòng)機(jī)運(yùn)行壽命間的相關(guān)性強(qiáng)弱,篩選出相關(guān)性較強(qiáng)的監(jiān)測(cè)參數(shù)用于實(shí)驗(yàn)驗(yàn)證;為降低監(jiān)測(cè)參數(shù)量綱對(duì)預(yù)測(cè)模型的影響,對(duì)監(jiān)測(cè)參數(shù)進(jìn)行歸一化處理。

2)CatBoost模型訓(xùn)練

構(gòu)建基于CatBoost算法的發(fā)動(dòng)機(jī)RUL預(yù)測(cè)訓(xùn)練模型,設(shè)置超參數(shù)及其尋優(yōu)范圍;利用BO算法對(duì)CatBoost訓(xùn)練模型中的超參數(shù)進(jìn)行優(yōu)化,尋優(yōu)過(guò)程中返回RMSE最小值及尋優(yōu)參數(shù)的取值。利用篩選出的最優(yōu)超參數(shù)作為預(yù)測(cè)模型的最終超參數(shù)組合。

3)剩余壽命預(yù)測(cè)

將經(jīng)過(guò)數(shù)據(jù)預(yù)處理后的監(jiān)測(cè)參數(shù)輸入到構(gòu)建的BO-CatBoost預(yù)測(cè)模型中進(jìn)行訓(xùn)練,返回預(yù)測(cè)評(píng)價(jià)結(jié)果RMSE,R2及MAE的值,實(shí)現(xiàn)航空發(fā)動(dòng)機(jī)的RUL預(yù)測(cè)。

3 實(shí)驗(yàn)驗(yàn)證

3.1 數(shù)據(jù)集介紹

本文用于實(shí)驗(yàn)的數(shù)據(jù)是來(lái)自美國(guó)國(guó)家航空航天局(national aeronautics and space administration,NASA)發(fā)布的渦扇發(fā)動(dòng)機(jī)商用模塊化航空推進(jìn)系統(tǒng)仿真數(shù)據(jù)集(commercial modular aero-propulsion system simulation,C-MAPSS)[21]。該數(shù)據(jù)集中包含4個(gè)子數(shù)據(jù)集FD001～FD004,每個(gè)子數(shù)據(jù)集都由訓(xùn)練集與測(cè)試集組成。訓(xùn)練集中包含了發(fā)動(dòng)機(jī)從初始運(yùn)行狀態(tài)到磨損失效后的全壽命周期數(shù)據(jù),測(cè)試集僅包含發(fā)動(dòng)機(jī)在故障發(fā)生前的部分運(yùn)行周期數(shù)據(jù)[22]。每個(gè)子數(shù)據(jù)集中包含飛行高度、馬赫數(shù)與油門桿解算器角度3種操作條件,運(yùn)行周期及21種航空發(fā)動(dòng)機(jī)傳感器監(jiān)測(cè)參數(shù)。本文選取FD001～FD004共4個(gè)子數(shù)據(jù)集進(jìn)行發(fā)動(dòng)機(jī)RUL實(shí)驗(yàn)。C-MAPSS數(shù)據(jù)集見表1。

表1 C-MAPSS數(shù)據(jù)集

圖2展示了FD001數(shù)據(jù)集中100臺(tái)發(fā)動(dòng)機(jī)的最大運(yùn)行周期分布情況。從圖2可知,在FD001數(shù)據(jù)集100臺(tái)發(fā)動(dòng)機(jī)中最小的運(yùn)行周期為128,最大的運(yùn)行周期為362,其他發(fā)動(dòng)機(jī)的運(yùn)行周期大部分分布在[145,250]范圍內(nèi)。

圖2 FD001數(shù)據(jù)集100臺(tái)發(fā)動(dòng)機(jī)的最大運(yùn)行周期分布情況

3.2 數(shù)據(jù)預(yù)處理

3.2.1 特征提取

不同的監(jiān)測(cè)參數(shù)對(duì)構(gòu)建的訓(xùn)練模型具有不同程度的影響,為有效提升發(fā)動(dòng)機(jī)RUL預(yù)測(cè)的精度,采用基于MIC的特征選擇方法對(duì)影響航空發(fā)動(dòng)機(jī)壽命周期的監(jiān)測(cè)參數(shù)進(jìn)行篩選。表2為FD001不同監(jiān)測(cè)參數(shù)與發(fā)動(dòng)機(jī)RUL的MIC計(jì)算結(jié)果。由表2可知,編號(hào)為1、5、10、16、18與19的監(jiān)測(cè)參數(shù)MIC值為0,且編號(hào)6的值接近于0,說(shuō)明這7個(gè)監(jiān)測(cè)參數(shù)對(duì)航空發(fā)動(dòng)機(jī)的壽命運(yùn)行周期相關(guān)性極弱,將這7種監(jiān)測(cè)參數(shù)從21種監(jiān)測(cè)參數(shù)中剔除,利用剩余的14種監(jiān)測(cè)進(jìn)行航空發(fā)動(dòng)機(jī)的RUL預(yù)測(cè)實(shí)驗(yàn)。

3.2.2 數(shù)據(jù)預(yù)處理

因航空發(fā)動(dòng)機(jī)不同的監(jiān)測(cè)參數(shù)具有不同的量綱,為了縮小監(jiān)測(cè)參數(shù)數(shù)值之間的差異,提高預(yù)測(cè)的效率及準(zhǔn)確率,本文選取最小最大歸一化公式對(duì)發(fā)動(dòng)機(jī)監(jiān)測(cè)參數(shù)進(jìn)行歸一化處理。其計(jì)算公式如下[23]:

3.3 基于貝葉斯優(yōu)化的CatBoost模型

因CatBoost訓(xùn)練模型中的超參數(shù)取值范圍不同,會(huì)對(duì)發(fā)動(dòng)機(jī)的RUL預(yù)測(cè)值帶來(lái)不同程度的影響。本文使用BO算法對(duì)CatBoost回歸預(yù)測(cè)模型中的超參數(shù)進(jìn)行尋優(yōu)。搭建BO-CatBoost訓(xùn)練模型并設(shè)置超參數(shù)取值范圍。表3為BO-CatBoost超參數(shù)尋優(yōu)結(jié)果。

表3 BO-CatBoost超參數(shù)尋優(yōu)結(jié)果

為獲得最佳的超參數(shù)值,本文分別訓(xùn)練預(yù)測(cè)模型中的超參數(shù)iterations與learning_rate在不同取值組合下的評(píng)估指標(biāo)RMSE值。圖3為預(yù)測(cè)模型不同超參數(shù)值的RMSE。由圖3可知,當(dāng)iterations為999,learning_rate為0.07時(shí),模型訓(xùn)練效果最好,RMSE值最小。

圖3 預(yù)測(cè)模型不同超參數(shù)值的RMSE

3.4 實(shí)驗(yàn)預(yù)測(cè)結(jié)果

圖4展示了FD001～FD004數(shù)據(jù)集下的航空發(fā)動(dòng)機(jī)RUL預(yù)測(cè)結(jié)果。圖4中,紅色曲線為發(fā)動(dòng)機(jī)的RUL真實(shí)值,藍(lán)色曲線為RUL預(yù)測(cè)值。由于FD001與FD003為在單一操作條件和故障模式下采集的數(shù)據(jù),發(fā)動(dòng)機(jī)數(shù)量少,預(yù)測(cè)結(jié)果較為稀疏。FD002與FD004為在多工況環(huán)境中采集的數(shù)據(jù),發(fā)動(dòng)機(jī)數(shù)量較多,與其他2個(gè)數(shù)據(jù)集的RUL預(yù)測(cè)相比,預(yù)測(cè)難度具有一定挑戰(zhàn)性。由對(duì)比結(jié)果可知,發(fā)動(dòng)機(jī)的RUL真實(shí)值與RUL預(yù)測(cè)值比較貼合,預(yù)測(cè)誤差較小,說(shuō)明本文所提模型的預(yù)測(cè)效果較好。

(a)FD001測(cè)試集

(b)FD002測(cè)試集

(c)FD003測(cè)試集

(d)FD004測(cè)試集

3.5 比較分析

為全面評(píng)估本文所提算法的性能,實(shí)驗(yàn)過(guò)程中選擇RMSE、MAE及R23個(gè)性能評(píng)價(jià)指標(biāo)來(lái)衡量各個(gè)預(yù)測(cè)模型的性能。表4展示了不同預(yù)測(cè)模型在FD001～FD004數(shù)據(jù)集中的RUL預(yù)測(cè)結(jié)果。由表4可知,與其他預(yù)測(cè)方法相比,采用MIC-BO-CatBoost方法進(jìn)行預(yù)測(cè)時(shí),評(píng)價(jià)指標(biāo)RMSE與MAE的值最小,R2也最接近于1。

圖5展示了MIC-BO-CatBoost預(yù)測(cè)模型在FD001～FD004數(shù)據(jù)集中的預(yù)測(cè)誤差箱線圖。箱子中間的實(shí)線表示RUL預(yù)測(cè)的期望,箱子的規(guī)模越小,表示模型在預(yù)測(cè)過(guò)程中的RUL預(yù)測(cè)結(jié)果的不確定性越低,準(zhǔn)確率越高。由圖5可知,本文所提預(yù)測(cè)模型MIC-BO-CatBoost在FD001與FD003數(shù)據(jù)集上的箱線圖規(guī)模較小,預(yù)測(cè)精確度較高。

表4 不同預(yù)測(cè)模型在FD001～FD004數(shù)據(jù)集中的RUL預(yù)測(cè)結(jié)果

圖5 MIC-BO-CatBoost預(yù)測(cè)模型在FD001～FD004數(shù)據(jù)集中的預(yù)測(cè)誤差箱線圖

為了對(duì)比不同模型在FD001～FD004數(shù)據(jù)集中的預(yù)測(cè)效果,本文構(gòu)建嶺(Ridge)回歸[24]、K近鄰回歸(KNN)算法[25]、極端梯度提升(XGBoost)算法[26]與本文所提方法進(jìn)行比較分析。分別選取FD001～FD004數(shù)據(jù)集中的第32號(hào)、3號(hào)、76號(hào)及56號(hào)發(fā)動(dòng)機(jī)進(jìn)行RUL預(yù)測(cè)。圖6展示了FD001～FD004在不同預(yù)測(cè)模型下的RUL預(yù)測(cè)結(jié)果?；疑珔^(qū)間代表RUL預(yù)測(cè)分布的95%置信區(qū)間。由圖6可知,其他預(yù)測(cè)模型的RUL預(yù)測(cè)結(jié)果大部分位于RUL真實(shí)值之上或之下,說(shuō)明預(yù)測(cè)結(jié)果表現(xiàn)出滯后或超前預(yù)測(cè),而本文所提方法的RUL預(yù)測(cè)值緊密圍繞RUL真實(shí)值波動(dòng),擬合程度較高。與其他3種預(yù)測(cè)模型的對(duì)比結(jié)果可知,本文所提方法的RUL預(yù)測(cè)值與真實(shí)值的偏差最小,且預(yù)測(cè)結(jié)果基本都被95%置信區(qū)間覆蓋,有效說(shuō)明本文所提方法能夠更好地用于航空發(fā)動(dòng)機(jī)的剩余使用壽命預(yù)測(cè)。

(a)測(cè)試集FD001第32號(hào)發(fā)動(dòng)機(jī)

(b)測(cè)試集FD002第3號(hào)發(fā)動(dòng)機(jī)

(c)測(cè)試集FD003第76號(hào)發(fā)動(dòng)機(jī)

(d)測(cè)試集FD004第56號(hào)發(fā)動(dòng)機(jī)

4 結(jié)論

本文針對(duì)航空發(fā)動(dòng)機(jī)非線性特征提取困難及預(yù)測(cè)精度不高等問(wèn)題,提出了一種基于MIC-BO-CatBoost的航空發(fā)動(dòng)機(jī)預(yù)測(cè)模型,使用C-MPASS發(fā)動(dòng)機(jī)退化數(shù)據(jù)集進(jìn)行驗(yàn)證和分析,得到以下結(jié)論:

1)通過(guò)MIC方法分析各個(gè)監(jiān)測(cè)參數(shù)對(duì)發(fā)動(dòng)機(jī)壽命運(yùn)行周期的相關(guān)性,提取能表征發(fā)動(dòng)機(jī)退化過(guò)程的非線性特征。

2)利用CatBoost算法對(duì)RUL預(yù)測(cè)模型進(jìn)行訓(xùn)練,有效解決了預(yù)測(cè)中的梯度偏差及預(yù)測(cè)偏移問(wèn)題。

3)使用BO-CatBoost算法建立預(yù)測(cè)模型,通過(guò)在FD001～FD004數(shù)據(jù)集上進(jìn)行實(shí)驗(yàn)評(píng)估后,結(jié)果表明與其他預(yù)測(cè)模型相比,本文所提方法的預(yù)測(cè)精度更加貼近發(fā)動(dòng)機(jī)的RUL真實(shí)值,證明了MIC-BO-CatBoost預(yù)測(cè)模型的有效性。