馬迪迪,趙靜,林亞龍,王婧雯
(合肥市軌道交通集團(tuán)有限公司,合肥 230001)
地鐵、高速列車等有軌設(shè)備的行駛速度逐漸提升、軌道網(wǎng)密度逐漸增大,車次數(shù)量越來(lái)越多,在此背景下,如果全自動(dòng)運(yùn)行的軌旁設(shè)備出現(xiàn)故障,會(huì)打亂鐵路網(wǎng)絡(luò)原始的調(diào)度計(jì)劃,產(chǎn)生一定的影響[1,2]。為了確保地鐵、高鐵等設(shè)備運(yùn)行的安全性,需要度量軌旁設(shè)備的可靠性。相關(guān)方法一直是研究的重點(diǎn):
葉遠(yuǎn)波[3]等首先分析了導(dǎo)致設(shè)備出現(xiàn)故障和失效的因素,將嚴(yán)重度、發(fā)生度和可檢度作為可靠性度量指標(biāo),對(duì)設(shè)備的可靠性展開(kāi)評(píng)估,該方法的可靠性度量結(jié)果與實(shí)際結(jié)果不符,存在度量精度低的問(wèn)題。楊墨緣[4]等分析了設(shè)備的工作原理以及拓?fù)浣Y(jié)構(gòu),根據(jù)分析結(jié)果劃分設(shè)備的任務(wù)剖面,在馬爾科夫過(guò)程的基礎(chǔ)上建立可靠性解析模型,以此度量設(shè)備的可靠性,該方法度量設(shè)備可靠性所需的時(shí)間較長(zhǎng),存在度量效率低的問(wèn)題。
為了解決上述方法中存在的問(wèn)題,提出多尺度特征融合下的軌旁設(shè)備可靠性度量方法。
全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備的多維源信號(hào)可以包括以下內(nèi)容:
1)聲音信號(hào):例如列車或其他交通工具的噪音或警示聲等。
2)視覺(jué)信號(hào):例如行駛方向指示燈、停車信號(hào)燈、交通信號(hào)燈等可視化信號(hào)。
3)電子信號(hào):例如傳感器數(shù)據(jù)、控制指令、通信信號(hào)等。
4)環(huán)境信號(hào):例如溫度、濕度、氣壓等環(huán)境參數(shù)。
5)電力信號(hào):例如電能消耗、供電質(zhì)量等。
不同的信號(hào)源可以提供關(guān)于設(shè)備狀態(tài)、環(huán)境條件和操作指令等多方面的信息,通過(guò)多維源信號(hào)的采集和處理,設(shè)備可以從不同的角度獲取反饋信息,使得設(shè)備可以更全面地了解和響應(yīng)實(shí)際情況。但是增加了全自動(dòng)運(yùn)行系統(tǒng)的復(fù)雜度和處理難度,并可能帶來(lái)更高的數(shù)據(jù)傳輸和存儲(chǔ)需求。對(duì)此,進(jìn)行多維源信號(hào)分解與重構(gòu)。
通過(guò)將多維源信號(hào)分解為單維源信號(hào),更靈活地對(duì)每個(gè)信號(hào)進(jìn)行處理和控制。
設(shè)x(n) 表示待分解的軌旁設(shè)備信號(hào),通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解獲得設(shè)備信號(hào)的第k個(gè)IMF 分量E k(x),采用互補(bǔ)完全經(jīng)驗(yàn)?zāi)B(tài)分解方法分解軌旁設(shè)備信號(hào)的具體過(guò)程如下:
1)用E1(x)=x-Q(x)表示軌旁設(shè)備的第1 個(gè)IMF 分量,其中Q(x) 代表的是包絡(luò)計(jì)算[5,6],針對(duì)信號(hào)x i(n) ,可根據(jù)下式計(jì)算第一個(gè)殘差 1r:
式中:
<·>—取均值。
2)根據(jù)步驟(1)計(jì)算得到的殘差 1r獲得軌旁設(shè)備信號(hào)的第一個(gè)模態(tài)IMF1=x-r1;
3)同理計(jì)算信號(hào)的第二個(gè)殘差,以此獲得軌旁設(shè)備信號(hào)的第二個(gè)模態(tài)IMF2=r1-r2;
4)按照相同的方法計(jì)算軌旁設(shè)備信號(hào)的第j個(gè)殘差rj,獲得對(duì)應(yīng)的第j個(gè)模態(tài),重復(fù)上述過(guò)程,直至獲得軌旁設(shè)備信號(hào)的全部IMFs。
采用獨(dú)立成分分析方法對(duì)上述分解得到的IMF 展開(kāi)獨(dú)立分析與重構(gòu),以此獲得軌旁設(shè)備的多維源信號(hào)。
此時(shí),即使其中某個(gè)信號(hào)出現(xiàn)故障或干擾,其他信號(hào)仍然能夠正常工作。
不同的應(yīng)用場(chǎng)景和設(shè)備可能對(duì)信號(hào)有不同的處理需求,因此進(jìn)行多維信號(hào)重構(gòu)可以根據(jù)具體的應(yīng)用場(chǎng)景,提取或抑制某些特定的信號(hào)源,提高信號(hào)的可理解性、滿足特定的處理需求。
通過(guò)下述過(guò)程計(jì)算軌旁設(shè)備多維源信號(hào)對(duì)應(yīng)的模糊熵:
1)給定軌旁設(shè)備的N維復(fù)合信號(hào)Y=[y1,y2, …,yN];
式中:
m≤N- 2—相空間的維數(shù)。
3)設(shè)代表的是相空間Yi與Yj之間存在的最大絕對(duì)距離,其計(jì)算公式如下:
4)構(gòu)建模糊隸屬函數(shù)A(Y):
式中:
r代表的是相似容限度。利用函數(shù)A(Y) 計(jì)算矢量之間的相似度:
5)設(shè)置函數(shù)φm(r),其表達(dá)式如下:
將復(fù)雜的多維源信號(hào)分解成更簡(jiǎn)單的單維源信號(hào)后,更好地理解信號(hào)的特征和含義,利用函數(shù)φm(r)計(jì)算軌旁設(shè)備各維源信號(hào)的模糊熵F(m,r,N):
設(shè)定下述模糊熵閾值判據(jù):
式中:
1<k≤N/2,β(k)—升序排序公式(7)計(jì)算得到的模糊熵后,其中的第k個(gè)模糊熵值。根據(jù)公式(8)對(duì)設(shè)備多維源信號(hào)展開(kāi)判斷,排在前k個(gè)熵值對(duì)應(yīng)的各維源信號(hào)屬于有用信號(hào),剩余信號(hào)屬于高頻噪聲,將其剔除,完成軌旁設(shè)備信號(hào)的去噪處理。
去噪后的多維源信號(hào)可以為全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量提供多方面的信息和智能化的決策支持。
去噪后的多維源信號(hào)可以在不同環(huán)境條件下保持穩(wěn)定,但是對(duì)于某些運(yùn)行狀態(tài)要求較高的軌旁設(shè)備,如列車自動(dòng)監(jiān)控系統(tǒng),軌旁設(shè)備可靠性度量信號(hào)解讀的不確定性擴(kuò)大,對(duì)此,需要通過(guò)多尺度特征融合,計(jì)算軌旁設(shè)備特征序列,消除量綱對(duì)可靠性度量結(jié)果產(chǎn)生的影響,即設(shè)計(jì)完成基于灰色關(guān)聯(lián)度的可靠性度量方法。
多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法將去噪后的軌旁設(shè)備信號(hào)輸入卷積神經(jīng)網(wǎng)絡(luò)中,以此完成信號(hào)特征的提取與融合。
將軌旁設(shè)備信號(hào)輸入網(wǎng)絡(luò)的輸入層中,其次利用卷積層局部卷積處理軌旁設(shè)備信號(hào)[9,10],輸出信號(hào)特征:
式中:
y—卷積層的輸出,即學(xué)習(xí)到的軌旁設(shè)備信號(hào)特征;
ω—權(quán)重矩陣;
x—輸入的軌旁設(shè)備信號(hào);
b—卷積核對(duì)應(yīng)的偏置;
*—卷積操作;
g(·) —激活函數(shù)。
在池化層中對(duì)軌旁設(shè)備信號(hào)特征y展開(kāi)全局平均池化處理,以此降低特征維度。將降維后的特征輸入網(wǎng)絡(luò)的批標(biāo)準(zhǔn)化層中,對(duì)特征展開(kāi)歸一化處理:
式中:
yn—特征的歸一化結(jié)果;
η、β—尺度位移參數(shù);
將歸一化處理后的軌旁設(shè)備信號(hào)特征yn輸入全連接層,利用Softmax 函數(shù)對(duì)特征展開(kāi)融合處理:
式中:
g(yi)—軌旁信號(hào)的多尺度特征融合結(jié)果;
Nc—分類總數(shù)[11,12]。
多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法采用灰色關(guān)聯(lián)分析方法完成設(shè)備的可靠性度量。
根據(jù)2.1 節(jié)的特征融合結(jié)果建立軌旁設(shè)備信號(hào)的特征序列0Y,結(jié)合軌旁設(shè)備的歷史運(yùn)行數(shù)據(jù)和狀態(tài)建立參考模式特征集Y={Ys|s=1,2, …,S},其中S代表的是軌旁設(shè)備總運(yùn)行狀態(tài)[13,14]。
灰色關(guān)聯(lián)分析方法通過(guò)計(jì)算軌旁設(shè)備特征序列0Y與Y={Yx|s=1,2, …,S}中參考狀態(tài)之間的灰色關(guān)聯(lián)度κ(Y0,Ys),根據(jù)計(jì)算結(jié)果判斷軌旁設(shè)備屬于哪類運(yùn)行狀態(tài),以此實(shí)現(xiàn)設(shè)備的可靠性度量。
用[y(1),y(2), …,y(q), …,y(Q)表示軌旁設(shè)備的特征序列,其中y(q) 代表的是軌旁設(shè)備特征序列中存在的第q維特征,Q代表的是特征序列對(duì)應(yīng)的維度。
為了消除量綱對(duì)可靠性度量結(jié)果產(chǎn)生的影響,采用極差最大化變換對(duì)y(q) 展開(kāi)量化處理:
式中:
經(jīng)量化處理后獲得的軌旁設(shè)備信號(hào)特征序列表示為Y0=[y0(1),y0(2), …,y0(q), …,y0(Q)],設(shè)置關(guān)聯(lián)系數(shù)κ0,s(q),其計(jì)算公式如下:
式中:
|y0(q)-y s(q)|—y0(q)與y s(q)的絕對(duì)差值;
ι—灰色關(guān)聯(lián)度的分辨系數(shù);
根據(jù)公式(14)的計(jì)算結(jié)果,判斷軌旁設(shè)備的運(yùn)行狀態(tài),完成設(shè)備的可靠性度量。
可靠性度量各項(xiàng)參數(shù)和指標(biāo)后,可以了解全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備是否正常運(yùn)行、存在故障或處于異常狀態(tài),便于及時(shí)采取故障恢復(fù)和糾正措施。
為了驗(yàn)證多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法的整體有效性,需要對(duì)其展開(kāi)相關(guān)測(cè)試。
實(shí)驗(yàn)所用全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備及運(yùn)行信號(hào),如圖1 所示。
分析圖1 可知,全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備及實(shí)驗(yàn)環(huán)境包括:
1)軌道交通正線線路,包括正線、停車線、折返線、聯(lián)絡(luò)線、出入段轉(zhuǎn)換軌、保護(hù)區(qū)段等;
2)各種線路的長(zhǎng)度、坡度、平曲線半徑、豎曲線半徑、線路限速等;
3)軌道交通車輛段、停車場(chǎng)線路,包括機(jī)車車庫(kù)、洗車庫(kù)、列檢庫(kù)、維修庫(kù)、牽出線等;
4)車輛,包括各級(jí)牽引、制動(dòng)曲線、車輛編組、構(gòu)造速度等;
5)運(yùn)營(yíng)時(shí)刻表和排班計(jì)劃;
6)各類信號(hào)設(shè)備和信號(hào)系統(tǒng);
7)牽引供電系統(tǒng)。
原始軌旁設(shè)備信號(hào)中存在大量的噪聲,不利于信號(hào)特征提取,現(xiàn)采用多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法對(duì)其展開(kāi)去噪處理,結(jié)果如圖2 所示。
圖2 信號(hào)去噪結(jié)果
由圖2 可知,采用所提方法展開(kāi)去噪可有效消除信號(hào)中存在的高頻分量,體現(xiàn)出信號(hào)的細(xì)節(jié)特征,表明所提方法具有良好的去噪性能。
引入波形相似系數(shù)NCC對(duì)所提方法的去噪效果展開(kāi)評(píng)價(jià),NCC在區(qū)間[0,1]內(nèi)取值,NCC越大,表明去噪效果越好:
式中:
s(i) —原始軌旁設(shè)備信號(hào);
y(i) —去噪后的設(shè)備信號(hào),
N—信號(hào)數(shù)量。
通過(guò)計(jì)算,所提方法的NCC高達(dá)0.92,表明經(jīng)過(guò)所提方法去噪后,信號(hào)中的細(xì)節(jié)信息得以保留,驗(yàn)證了所提方法具有良好的去噪效果。
采用多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法展開(kāi)設(shè)備可靠性度量測(cè)試,利用F1 指標(biāo)評(píng)價(jià)上述方法的度量精度:
式中:
P—精確率;
R—召回率;
TP—度量正確的正樣本數(shù);
FP—度量錯(cuò)誤的負(fù)樣本數(shù);
FN—度量錯(cuò)誤的正樣本數(shù)。
分析圖3 可知,在設(shè)備可靠性度量測(cè)試中,所提方法的F1 指標(biāo)值遠(yuǎn)高于文獻(xiàn)[3]方法和文獻(xiàn)[4]方法,表明所提方法在度量過(guò)程中具有較高的精確率和召回率,因?yàn)樗岱椒ㄔ谠O(shè)備可靠性度量之前,分解了設(shè)備信號(hào),消除信號(hào)中存在的高頻噪聲,確保度量結(jié)果不受到噪聲的干擾,提高了所提方法的度量精度。
圖3 不同方法的F1 指標(biāo)
選取10 個(gè)軌旁設(shè)備,采用所提方法、文獻(xiàn)[3]方法和文獻(xiàn)[4]方法對(duì)其展開(kāi)可靠性度量,對(duì)比三種方法所用的時(shí)間,結(jié)果如表1 所示。
表1 度量時(shí)間
根據(jù)表1 中的數(shù)據(jù)可知,所提方法可在短時(shí)間內(nèi)完成相同數(shù)量設(shè)備的可靠性度量,表明所提方法具有較高的度量效率。
針對(duì)目前設(shè)備可靠性度量方法存在的度量精度低和效率低的問(wèn)題,提出多尺度特征融合下的全自動(dòng)運(yùn)行信號(hào)軌旁設(shè)備可靠性度量方法,采用互補(bǔ)完全經(jīng)驗(yàn)?zāi)B(tài)分解方法將復(fù)雜的多維源信號(hào)分解成更簡(jiǎn)單的單維源信號(hào),重構(gòu)軌旁設(shè)備信號(hào)的相空間,去噪處理軌旁設(shè)備信號(hào),通過(guò)多尺度特征融合,計(jì)算軌旁設(shè)備特征序列,計(jì)及參考狀態(tài)之間的灰色關(guān)聯(lián)度,改進(jìn)多尺度特征融合過(guò)程,可靠性度量各項(xiàng)參數(shù)和指標(biāo),判斷軌旁設(shè)備的運(yùn)行狀態(tài),在信號(hào)去噪、度量精度和度量效率等方面表現(xiàn)出良好的性能。