牛國(guó)臣，呂波漾

(中國(guó)民航大學(xué) 機(jī)器人研究所,天津 300300)

0 引言

近年來(lái)我國(guó)民航運(yùn)輸機(jī)場(chǎng)的旅客吞吐量持續(xù)增長(zhǎng)，行李托運(yùn)業(yè)務(wù)的辦理效率影響旅客的出行體驗(yàn)，為增強(qiáng)機(jī)場(chǎng)運(yùn)輸系統(tǒng)的工作能力、響應(yīng)智慧機(jī)場(chǎng)建設(shè)號(hào)召[1-2]，自助行李托運(yùn)設(shè)備開始在國(guó)內(nèi)機(jī)場(chǎng)廣泛應(yīng)用。自助行李托運(yùn)設(shè)備能為旅客辦理自助值機(jī)和行李托運(yùn)業(yè)務(wù)，為保障設(shè)備的健康運(yùn)行，采集設(shè)備監(jiān)測(cè)數(shù)據(jù)，及時(shí)了解設(shè)備部件系統(tǒng)的運(yùn)行狀態(tài)，指導(dǎo)維護(hù)人員采取相應(yīng)措施具有重要的意義[3]。

自助行李托運(yùn)設(shè)備是復(fù)雜機(jī)電設(shè)備，現(xiàn)有的設(shè)備狀態(tài)評(píng)估方法可分為三類：基于傳感器數(shù)據(jù)的退化模型法、基于運(yùn)行日志的數(shù)據(jù)驅(qū)動(dòng)法以及混合模型法[4-6]。退化模型法是指采集部件工作時(shí)的振動(dòng)、溫度等信息，利用設(shè)備的物理特性構(gòu)建退化模型[7-9]，常用于軸承、電機(jī)轉(zhuǎn)子等機(jī)械結(jié)構(gòu)較為簡(jiǎn)單的機(jī)械設(shè)備，任子強(qiáng)等人[10]針對(duì)單一傳感器監(jiān)測(cè)方法存在效率低、精度低的缺點(diǎn)，提出了一種多傳感器融合的航空發(fā)動(dòng)機(jī)壽命預(yù)測(cè)方法，但仍然存在監(jiān)測(cè)數(shù)據(jù)利用不充分、模型適應(yīng)性弱的不足。數(shù)據(jù)驅(qū)動(dòng)法則是利用統(tǒng)計(jì)學(xué)或機(jī)器學(xué)習(xí)等理論[11]，根據(jù)運(yùn)行日志數(shù)據(jù)建立設(shè)備退化映射關(guān)系[12-14]，文獻(xiàn)[15]提出一種基于長(zhǎng)短期時(shí)間特性的時(shí)間卷積特征神經(jīng)網(wǎng)絡(luò)用于預(yù)測(cè)設(shè)備剩余壽命，王衛(wèi)華等人[16]提出一種基于日志聚類的多故障預(yù)測(cè)方法，使用層次聚類算法挖掘與故障時(shí)間相關(guān)的時(shí)間序列，進(jìn)行故障預(yù)測(cè)的同時(shí)平衡準(zhǔn)確率與召回率，但此方法可能會(huì)生成相同的故障預(yù)測(cè)規(guī)則，工程應(yīng)用上有一定困難?；诨旌夏Ｐ偷脑O(shè)備狀態(tài)評(píng)估方法尚處于研究階段，文獻(xiàn)[17]通過(guò)對(duì)電梯事件型數(shù)據(jù)進(jìn)行分析，提出一種基于比例風(fēng)險(xiǎn)模型與機(jī)器學(xué)習(xí)混合的電梯剩余壽命預(yù)測(cè)方法。文獻(xiàn)[18]利用設(shè)備健康狀態(tài)信息預(yù)測(cè)剩余使用壽命，提出基于剩余壽命預(yù)測(cè)的維修與備件訂購(gòu)聯(lián)合策略，以降低設(shè)備檢修成本和備件成本。

由于自助行李托運(yùn)設(shè)備實(shí)際投入時(shí)間不長(zhǎng)，產(chǎn)生的狀態(tài)型數(shù)據(jù)類型不足夠豐富，而傳統(tǒng)設(shè)備狀態(tài)評(píng)估方法未能充分結(jié)合狀態(tài)型和事件型兩類數(shù)據(jù)進(jìn)行評(píng)估，因此本文提出一種融合Cox模型與維納過(guò)程的設(shè)備狀態(tài)評(píng)估方法，首先構(gòu)建基于Cox模型的狀態(tài)突變模型，獲得在風(fēng)險(xiǎn)因素影響下設(shè)備發(fā)生突變失效的概率；然后構(gòu)建基于維納過(guò)程的狀態(tài)漸變模型，獲得表征退化的復(fù)合退化指標(biāo)，預(yù)測(cè)設(shè)備的綜合健康狀態(tài)值；最后結(jié)合設(shè)備維護(hù)經(jīng)驗(yàn)得到設(shè)備整體的定性狀態(tài)，充分利用兩類數(shù)據(jù)對(duì)設(shè)備健康狀態(tài)進(jìn)行評(píng)估，提高了數(shù)據(jù)利用率，彌補(bǔ)因僅用狀態(tài)型數(shù)據(jù)評(píng)估造成的精度不足。

1 狀態(tài)評(píng)估體系設(shè)計(jì)

1.1 設(shè)備監(jiān)測(cè)數(shù)據(jù)分析

為確保自助行李托運(yùn)設(shè)備的正常運(yùn)行，工作人員會(huì)定期巡檢設(shè)備，記錄包含維護(hù)時(shí)間、設(shè)備編號(hào)、異常情況和采取措施等事件型巡檢信息。同時(shí)，自助行李托運(yùn)設(shè)備有自檢功能，設(shè)備自檢系統(tǒng)可以監(jiān)測(cè)子設(shè)備的各類反饋數(shù)據(jù)，包含反映子設(shè)備工作是否正常的事件型數(shù)據(jù)，以及反映子設(shè)備性能良莠的狀態(tài)型數(shù)據(jù)。為便于分析，選取影響自助行李托運(yùn)設(shè)備正常工作的關(guān)鍵子設(shè)備，歸為值機(jī)交互子系統(tǒng)、通道擺閘門子系統(tǒng)和行李運(yùn)輸子系統(tǒng)，分別編號(hào)1，2，3。其中，值機(jī)交互子系統(tǒng)監(jiān)測(cè)打印機(jī)、觸摸屏等部件的反饋數(shù)據(jù)；通道擺閘門子系統(tǒng)監(jiān)測(cè)驅(qū)動(dòng)電機(jī)、擺閘門等部件的反饋數(shù)據(jù)；行李運(yùn)輸子系統(tǒng)監(jiān)測(cè)行李運(yùn)輸機(jī)、傳輸皮帶等部件的反饋數(shù)據(jù)。自助行李托運(yùn)設(shè)備各子系統(tǒng)的詳細(xì)監(jiān)測(cè)指標(biāo)如表1所示。

表1 自助行李托運(yùn)設(shè)備監(jiān)測(cè)指標(biāo)

1.2 設(shè)備狀態(tài)評(píng)估流程

基于自助行李托運(yùn)設(shè)備的兩類監(jiān)測(cè)數(shù)據(jù)，構(gòu)建突變狀態(tài)模型和漸變狀態(tài)模型，設(shè)計(jì)設(shè)備狀態(tài)評(píng)估方法，如圖1所示。突變狀態(tài)模型將子系統(tǒng)事件型監(jiān)測(cè)數(shù)據(jù)與對(duì)應(yīng)設(shè)備運(yùn)行時(shí)間作為輸入，獲得受風(fēng)險(xiǎn)事件協(xié)變量影響下設(shè)備正常運(yùn)行的生存函數(shù)與當(dāng)前時(shí)刻失效概率。結(jié)合各機(jī)場(chǎng)投入使用的同類設(shè)備失效概率分布給出每個(gè)子系統(tǒng)的失效概率閾值，判斷其在當(dāng)前時(shí)刻條件下是否發(fā)生失效。如果子系統(tǒng)沒(méi)有發(fā)生功能失效,則轉(zhuǎn)至對(duì)應(yīng)漸變狀態(tài)模型。漸變狀態(tài)模型將子系統(tǒng)多維狀態(tài)型監(jiān)測(cè)數(shù)據(jù)作為輸入，定義復(fù)合退化指標(biāo)表征子系統(tǒng)的退化量，建立基于維納過(guò)程的子系統(tǒng)性能退化模型預(yù)測(cè)其健康狀態(tài)值。根據(jù)設(shè)備運(yùn)行狀態(tài)對(duì)照區(qū)間得到子系統(tǒng)運(yùn)行狀態(tài)。最后綜合各子系統(tǒng)運(yùn)行狀態(tài)得到自助行李托運(yùn)設(shè)備整體運(yùn)行狀態(tài)，分別為正常、注意、警告、失效4種狀態(tài)。

圖1 自助行李托運(yùn)設(shè)備狀態(tài)評(píng)估體系

2 突變狀態(tài)模型構(gòu)建

2.1 Cox回歸模型

Cox回歸是由英國(guó)統(tǒng)計(jì)學(xué)家D.R.Cox提出的一種回歸模型[19]，又稱比例風(fēng)險(xiǎn)回歸模型，可同時(shí)分析眾多因素對(duì)項(xiàng)目生存周期的影響。模型基本形式為：

h(t,X)=h0(t)eβ1X1+β2X2+…+βmXm

(1)

式中，X=[X1,X2,…,Xm]為風(fēng)險(xiǎn)事件協(xié)變量；β=[β1,β2,…,βm]為回歸系數(shù)。h0(t)為基準(zhǔn)風(fēng)險(xiǎn)函數(shù)，表達(dá)式為：

(2)

式中，t為時(shí)間，n為觀測(cè)量的總數(shù)，λ為t時(shí)刻存在的風(fēng)險(xiǎn)值。

對(duì)風(fēng)險(xiǎn)函數(shù)積分可得到累積風(fēng)險(xiǎn)函數(shù)：

(3)

累積風(fēng)險(xiǎn)函數(shù)的值代表設(shè)備運(yùn)行至當(dāng)前時(shí)刻，且受協(xié)變量影響時(shí)，子系統(tǒng)即將發(fā)生失效異常的概率，使用其作為評(píng)價(jià)子系統(tǒng)是否發(fā)生突變的指標(biāo)。

2.2 模型輸入分析與參數(shù)求解

自助行李托運(yùn)設(shè)備事件型監(jiān)測(cè)數(shù)據(jù)中存在的缺少完整隨訪時(shí)間和生存結(jié)果的數(shù)據(jù)，稱為刪失數(shù)據(jù)。根據(jù)刪失狀態(tài)E可以將事件型數(shù)據(jù)分為兩類：記錄完整異常的數(shù)據(jù)，刪失狀態(tài)E=0；未記錄到異常發(fā)生卻截?cái)嗟臄?shù)據(jù)，刪失狀態(tài)E=1。

從設(shè)備事件型監(jiān)測(cè)數(shù)據(jù)中提取每個(gè)子系統(tǒng)的風(fēng)險(xiǎn)事件協(xié)變量，定義子系統(tǒng)j的風(fēng)險(xiǎn)事件協(xié)變量為Xj=[Xj1,Xj2,Xj3,Xj4,Xj5,Xj6]。它由兩部分組成，第一部分與外部環(huán)境有關(guān)，對(duì)應(yīng)的回歸系數(shù)βj1-βj2代表了不同設(shè)備運(yùn)行環(huán)境對(duì)子系統(tǒng)失效發(fā)生概率的影響；第二部分與設(shè)備業(yè)務(wù)狀態(tài)有關(guān)，對(duì)應(yīng)回歸系數(shù)βj3-βj6代表了不同部件運(yùn)行狀態(tài)對(duì)子系統(tǒng)失效發(fā)生概率的影響。各子系統(tǒng)風(fēng)險(xiǎn)事件協(xié)變量說(shuō)明如表2～4所示。

表2 值機(jī)交互子系統(tǒng)協(xié)變量

表3 通道擺閘門子系統(tǒng)協(xié)變量

表4 行李運(yùn)輸系統(tǒng)協(xié)變量

使用上述數(shù)據(jù)訓(xùn)練比例風(fēng)險(xiǎn)模型，結(jié)合刪失數(shù)據(jù)的模型似然函數(shù)表達(dá)式為：

(4)

式中，I為刪失狀態(tài)E=0的數(shù)據(jù)集合，J為刪失狀態(tài)E=1的數(shù)據(jù)集合。使用極大似然估計(jì)可得到模型參數(shù)β,γ,η估算值。

3 漸變狀態(tài)模型構(gòu)建

3.1 子系統(tǒng)復(fù)合退化指標(biāo)定義

令xi,j(t)表示第i(i=1,2,...,M)臺(tái)自助行李托運(yùn)設(shè)備中第j(j=1,2,3)子系統(tǒng)在t時(shí)刻融合多維狀態(tài)型監(jiān)測(cè)數(shù)據(jù)的復(fù)合退化指標(biāo)，表征該時(shí)刻設(shè)備子系統(tǒng)的性能退化量，M為被監(jiān)測(cè)自助行李托運(yùn)設(shè)備的數(shù)量。因此子系統(tǒng)復(fù)合退化指標(biāo)可表示為：

xi,j(t)=Zi,j(t)·Wj

(5)

式中，Zi,j(t)=[zi,j,1(t),…,zi,j,k(t),…,zi,j,S(t)]為第i臺(tái)設(shè)備中子系統(tǒng)j狀態(tài)型監(jiān)測(cè)數(shù)據(jù)向量，每個(gè)子系統(tǒng)有S個(gè)狀態(tài)型監(jiān)測(cè)數(shù)據(jù)，k=1,2,...,S，zi,j,k(t)表示第i臺(tái)設(shè)備中子系統(tǒng)j的第k個(gè)狀態(tài)型指標(biāo)在t時(shí)刻的監(jiān)測(cè)值；Wj=[ωj,1,…,ωj,k,…,ωj,s]T是融合系數(shù)向量，衡量各監(jiān)測(cè)值在融合過(guò)程中權(quán)重。

由于設(shè)備性能退化過(guò)程存在隨機(jī)性，每臺(tái)自助行李托運(yùn)設(shè)備的失效閾值存在差異。為減少設(shè)備漸變狀態(tài)監(jiān)測(cè)的不確定性，令所求子系統(tǒng)失效閾值與其他設(shè)備子系統(tǒng)真實(shí)失效時(shí)刻復(fù)合退化指標(biāo)的方差最小，得到設(shè)備失效時(shí)刻復(fù)合退化指標(biāo)的平均值為最優(yōu)解。定義P為設(shè)備子系統(tǒng)的失效閾值，表達(dá)式如下：

(6)

3.2 基于維納過(guò)程的子系統(tǒng)性能退化建模

采用維納過(guò)程對(duì)自助行李托運(yùn)設(shè)備子系統(tǒng)性能的退化過(guò)程建模，模型表達(dá)式為：

(7)

(8)

(9)

3.3 復(fù)合退化指標(biāo)融合系數(shù)求解

對(duì)基于維納過(guò)程的設(shè)備子系統(tǒng)隨機(jī)退化過(guò)程模型，定義隨機(jī)變量Ti,j表示設(shè)備i子系統(tǒng)j首次達(dá)到失效閾值P的時(shí)間為：

(10)

可以得到Ti,j的數(shù)學(xué)期望表達(dá)式為:

(11)

定義目標(biāo)函數(shù)：

(12)

式中，Γi,j為設(shè)備i子系統(tǒng)j的真實(shí)失效時(shí)間。將式(8)和式(11)代入式(12)，可以得：

(13)

采用非線性規(guī)劃方法尋找最優(yōu)解，得到設(shè)備子系統(tǒng)j復(fù)合退化指標(biāo)的融合系數(shù)Wj。

3.4 設(shè)備漸變狀態(tài)評(píng)估

定義Li,j,n為設(shè)備i子系統(tǒng)j在tn時(shí)刻的健康狀態(tài)值，其表達(dá)式為：

(14)

根據(jù)定義可知Li,j,n的數(shù)值越高，設(shè)備i子系統(tǒng)j在tn時(shí)刻的性能越好。由tn+ln=Ti,j可以求解得到Li,j,n。選取tn時(shí)刻值機(jī)交互子系統(tǒng)、通道擺閘門子系統(tǒng)和行李運(yùn)輸子系統(tǒng)中的最小健康狀態(tài)值表征當(dāng)前時(shí)刻自助行李托運(yùn)設(shè)備整體的健康度：

Hi=minLi,j

(15)

式中,Hi為最終得到的設(shè)備綜合健康狀態(tài)值。

4 設(shè)備綜合狀態(tài)評(píng)估

構(gòu)建突變狀態(tài)模型與漸變狀態(tài)模型后，可實(shí)現(xiàn)設(shè)備在t時(shí)刻的實(shí)時(shí)狀態(tài)評(píng)估，其步驟為：

1)首先進(jìn)入突變狀態(tài)模型，提取風(fēng)險(xiǎn)事件協(xié)變量，結(jié)合設(shè)備運(yùn)行時(shí)間作為輸入，分別輸出3個(gè)關(guān)鍵子系統(tǒng)的失效風(fēng)險(xiǎn)概率。結(jié)合設(shè)備制造商指導(dǎo)經(jīng)驗(yàn)設(shè)定子系統(tǒng)狀態(tài)突變閾值分別為80%、83%和83%。若設(shè)備3個(gè)子系統(tǒng)失效風(fēng)險(xiǎn)概率都低于對(duì)應(yīng)狀態(tài)突變風(fēng)險(xiǎn)閾值，即都未發(fā)生失效狀態(tài)突變，進(jìn)入漸變狀態(tài)模型評(píng)估設(shè)備整體健康狀態(tài)，否則認(rèn)為設(shè)備發(fā)生失效。

2)使用狀態(tài)型監(jiān)測(cè)數(shù)據(jù)作為漸變狀態(tài)模型的輸入，得到當(dāng)前時(shí)刻的設(shè)備綜合健康狀態(tài)值,其值越大設(shè)備性能越佳，定義S={S1,S2,S3,S4}4種設(shè)備運(yùn)行狀態(tài)，分別為“失效，警告，注意，正常”。定性標(biāo)準(zhǔn)如表5所示。

表5 運(yùn)行狀態(tài)定性標(biāo)準(zhǔn)

定性描述設(shè)備的運(yùn)行狀態(tài)，及時(shí)調(diào)整維護(hù)策略，提高設(shè)備的使用效率。

5 實(shí)驗(yàn)分析

為驗(yàn)證本文提出的自助行李托運(yùn)設(shè)備狀態(tài)評(píng)估方法，使用商用模塊化航空推進(jìn)系統(tǒng)仿真數(shù)據(jù)集(C-MAPSS，commercial modular aero-propulsion system simulation)和自助行李托運(yùn)設(shè)備運(yùn)行數(shù)據(jù)集作為設(shè)備狀態(tài)評(píng)估方法的輸入來(lái)驗(yàn)證模型的有效性。

5.1 C-MAPSS數(shù)據(jù)集驗(yàn)證

C-MAPSS數(shù)據(jù)集包含四組子數(shù)據(jù)集，對(duì)應(yīng)不同工況下飛機(jī)渦扇發(fā)動(dòng)機(jī)機(jī)電系統(tǒng)模擬產(chǎn)生的退化數(shù)據(jù)[20]。數(shù)據(jù)包括引擎單元號(hào)、時(shí)間戳、3種配置變量以及21個(gè)傳感器監(jiān)測(cè)數(shù)據(jù)。在本實(shí)驗(yàn)中，首先使用訓(xùn)練集數(shù)據(jù)建立突變和漸變狀態(tài)模型，然后使用測(cè)試集數(shù)據(jù)驗(yàn)證模型精確性。

5.1.1 訓(xùn)練集數(shù)據(jù)建模

C-MAPSS數(shù)據(jù)集中不區(qū)分子系統(tǒng)，故該實(shí)驗(yàn)中子系統(tǒng)下標(biāo)j省略。提取發(fā)動(dòng)機(jī)設(shè)置協(xié)變量X=[X1,X2,X3]作為Cox回歸模型的輸入，協(xié)變量說(shuō)明如表6所示。

表6 發(fā)動(dòng)機(jī)設(shè)置協(xié)變量說(shuō)明

使用數(shù)據(jù)集中每臺(tái)發(fā)動(dòng)機(jī)運(yùn)行至失效的時(shí)間周期值作為目標(biāo)隨訪時(shí)間輸入，建模得到發(fā)動(dòng)機(jī)累積風(fēng)險(xiǎn)概率曲線，單臺(tái)發(fā)動(dòng)機(jī)累積風(fēng)險(xiǎn)概率隨運(yùn)行時(shí)間增長(zhǎng)的變化曲線如圖2所示。

圖2 不同配置條件下的累積風(fēng)險(xiǎn)函數(shù)

圖2中橫軸表示發(fā)動(dòng)機(jī)運(yùn)行周期，縱軸表示發(fā)動(dòng)機(jī)將發(fā)生失效的風(fēng)險(xiǎn)概率。顯然，受設(shè)置協(xié)變量的影響，不同運(yùn)行配置下發(fā)動(dòng)機(jī)的失效概率不同。

原始狀態(tài)型數(shù)據(jù)經(jīng)過(guò)高斯濾波與標(biāo)準(zhǔn)化處理后作為本文方法漸變狀態(tài)模型的輸入Zi(t)。設(shè)定融合系數(shù)初值為Winit=[1,1,...,1,1]T，采用非線性規(guī)劃方法得到最優(yōu)融合系數(shù)：

W1=[0.0084 0.1574 -0.0944 1.5073

-0.0016 0.0052 -0.5198 0.4752 1.2362

-0.0010 1.6497 -2.1277 0.2016 -0.8603

1.9436 -0.0080 1.2963 0.0001 0.0043

-0.8646 -1.2576]T

圖3 發(fā)動(dòng)機(jī)復(fù)合退化指標(biāo)退化效果

5.1.2 模型精度分析

提取測(cè)試集數(shù)據(jù)的設(shè)置協(xié)變量X，令每臺(tái)發(fā)動(dòng)機(jī)的當(dāng)前運(yùn)行周期與對(duì)應(yīng)剩余壽命的加和ti=t_testi,N+RULi作為目標(biāo)隨訪時(shí)間，輸入由訓(xùn)練集構(gòu)建的突變狀態(tài)模型獲得測(cè)試集中發(fā)動(dòng)機(jī)失效時(shí)刻的風(fēng)險(xiǎn)概率，經(jīng)統(tǒng)計(jì)設(shè)實(shí)驗(yàn)中發(fā)動(dòng)機(jī)突變概率閾值為85%，若高于該閾值，則認(rèn)為判斷正確，結(jié)果如表7所示。

表7 突變模型測(cè)試集驗(yàn)證

若發(fā)動(dòng)機(jī)當(dāng)前周期的失效風(fēng)險(xiǎn)概率未超過(guò)設(shè)定閾值，則進(jìn)入漸變狀態(tài)模型評(píng)估發(fā)動(dòng)機(jī)健康狀態(tài)。將測(cè)試集多維傳感器數(shù)據(jù)作為漸變狀態(tài)模型的輸入，得到每臺(tái)發(fā)動(dòng)機(jī)的健康狀態(tài)值，與每臺(tái)發(fā)動(dòng)機(jī)真實(shí)剩余壽命進(jìn)行對(duì)比，其中測(cè)試集一的100臺(tái)發(fā)動(dòng)機(jī)的評(píng)估效果如圖4所示。

圖4 健康狀態(tài)評(píng)估值與真實(shí)值對(duì)比

為定量分析與評(píng)價(jià)預(yù)測(cè)的精確度，本文采用均方根誤差(RMSE，root mean squared error)作為評(píng)價(jià)指標(biāo)來(lái)描述誤差值大小，并與基于單一傳感器11號(hào)傳感器監(jiān)測(cè)值的設(shè)備狀態(tài)評(píng)估方法進(jìn)行對(duì)比，4個(gè)測(cè)試集得到的RMSE值如表8所示。

表8 漸變模型預(yù)測(cè)精度評(píng)價(jià)

可以發(fā)現(xiàn)，多維監(jiān)測(cè)傳感器數(shù)據(jù)融合的設(shè)備漸變狀態(tài)模型比起單傳感器評(píng)估模型有較高的精度，可以細(xì)化設(shè)備的運(yùn)行健康狀況。

5.2 自助行李托運(yùn)設(shè)備運(yùn)行數(shù)據(jù)集驗(yàn)證

為進(jìn)一步驗(yàn)證方法的可行性，使用大興機(jī)場(chǎng)南航值機(jī)區(qū)域的自助行李托運(yùn)設(shè)備在2019年6月到2020年6月期間的監(jiān)測(cè)數(shù)據(jù)進(jìn)行實(shí)際效果驗(yàn)證。

5.2.1 數(shù)據(jù)處理與分析

對(duì)照1.1節(jié)自助行李托運(yùn)設(shè)備監(jiān)測(cè)指標(biāo)表對(duì)人工巡檢數(shù)據(jù)信息排序并與設(shè)備自檢系統(tǒng)獲得的數(shù)據(jù)合并，得到70組共計(jì)245 280條監(jiān)測(cè)數(shù)據(jù)。按照2.2節(jié)中子系統(tǒng)協(xié)變量提取操作對(duì)設(shè)備事件型監(jiān)測(cè)數(shù)據(jù)處理，得到各子系統(tǒng)的風(fēng)險(xiǎn)事件協(xié)變量。以位于靠近進(jìn)站口的一臺(tái)運(yùn)行時(shí)長(zhǎng)為364小時(shí)的自助行李托運(yùn)設(shè)備為例，處理后的事件型數(shù)據(jù)與風(fēng)險(xiǎn)事件協(xié)變量如表9和表10所示。

表9 子系統(tǒng)事件型數(shù)據(jù)

表10 子系統(tǒng)協(xié)變量

自助行李托運(yùn)設(shè)備狀態(tài)型監(jiān)測(cè)數(shù)據(jù)進(jìn)行濾波和標(biāo)準(zhǔn)化處理，部分原始數(shù)據(jù)如表11～13所示。

表11 值機(jī)交互子系統(tǒng)狀態(tài)型數(shù)據(jù)

表12 通道擺閘門子系統(tǒng)狀態(tài)型數(shù)據(jù)

表13 行李運(yùn)輸子系統(tǒng)狀態(tài)型數(shù)據(jù)

5.2.2 方法有效性驗(yàn)證

構(gòu)建基于Cox回歸的設(shè)備突變狀態(tài)模型，分別得到3個(gè)子系統(tǒng)的回歸系數(shù)βj1-βj6，如表14所示。

表14 各子系統(tǒng)回歸系數(shù)

回歸系數(shù)為負(fù)值時(shí)，表示該協(xié)變量取值相比于參照類或基準(zhǔn)風(fēng)險(xiǎn)會(huì)使風(fēng)險(xiǎn)概率升高，系數(shù)越小代表協(xié)變量產(chǎn)生風(fēng)險(xiǎn)的占比越高，可知行李運(yùn)輸子系統(tǒng)出現(xiàn)卡包事件會(huì)使設(shè)備產(chǎn)生失效風(fēng)險(xiǎn)的概率增加。再基于各子系統(tǒng)狀態(tài)型數(shù)據(jù)構(gòu)建子系統(tǒng)復(fù)合退化指標(biāo)和基于維納過(guò)程的漸變狀態(tài)模型，得到模型參數(shù)α和σ。構(gòu)建模型后，將用于測(cè)試的自助行李托運(yùn)設(shè)備監(jiān)測(cè)信息作為輸入量，輸入設(shè)備狀態(tài)評(píng)估模型得到設(shè)備綜合健康狀態(tài)值與定性狀態(tài)。

某臺(tái)設(shè)備在2020年6月15日10時(shí)失效，工作人員對(duì)其進(jìn)行了維護(hù)。使用該設(shè)備失效前10天的監(jiān)測(cè)數(shù)據(jù)作為本文狀態(tài)評(píng)估方法的輸入，得到運(yùn)行時(shí)間為[490,740]時(shí)，設(shè)備失效風(fēng)險(xiǎn)概率未超過(guò)突變閾值。得到設(shè)備綜合健康狀態(tài)值與對(duì)應(yīng)狀態(tài)，對(duì)比使用基于單一事件型數(shù)據(jù)中打印機(jī)數(shù)據(jù)設(shè)備狀態(tài)評(píng)估的結(jié)果，效果如圖5所示。

圖5 設(shè)備狀態(tài)評(píng)估結(jié)果

圖5中左坐標(biāo)軸對(duì)應(yīng)本文方法獲得的設(shè)備綜合健康狀態(tài)值曲線，右坐標(biāo)軸對(duì)應(yīng)單一指標(biāo)評(píng)估的二值判斷曲線?？芍趩我恢笜?biāo)的設(shè)備狀態(tài)評(píng)估方法因無(wú)法細(xì)化設(shè)備運(yùn)行狀態(tài)，評(píng)估結(jié)果容易造成虛警，本文提出的融合Cox回歸與維納過(guò)程的設(shè)備狀態(tài)評(píng)估方法可有效利用事件型數(shù)據(jù)和狀態(tài)型數(shù)據(jù)，并對(duì)自助行李托運(yùn)設(shè)備狀態(tài)進(jìn)行評(píng)估，為及時(shí)調(diào)整設(shè)備維護(hù)策略提供了決策依據(jù)。

6 結(jié)束語(yǔ)

針對(duì)自助行李托運(yùn)設(shè)備運(yùn)行狀態(tài)的綜合評(píng)估，本文以設(shè)備事件型數(shù)據(jù)和狀態(tài)型數(shù)據(jù)為研究對(duì)象，建立監(jiān)測(cè)指標(biāo)體系。提出一種設(shè)備健康狀態(tài)評(píng)估方法。構(gòu)建基于Cox回歸的狀態(tài)突變模型，獲得風(fēng)險(xiǎn)事件協(xié)變量影響下的設(shè)備失效概率；定義子系統(tǒng)復(fù)合退化指標(biāo)，建立維納退化模型模擬設(shè)備狀態(tài)漸變過(guò)程，預(yù)測(cè)設(shè)備綜合健康狀態(tài)值；最終得到整個(gè)設(shè)備的定性狀態(tài)和相應(yīng)維護(hù)策略。經(jīng)實(shí)驗(yàn)驗(yàn)證，本文方法與單一指標(biāo)方法相比，提高設(shè)備數(shù)據(jù)利用率，不僅減少了虛警現(xiàn)象，而且提高了自助行李托運(yùn)設(shè)備的狀態(tài)評(píng)估精度。