任宇超,徐永能,喬 僑
(南京理工大學, 南京 210094)
隨著全國人口不斷增加,隨之而來的交通負擔也越來越重,緩解交通負擔刻不容緩。 隨著時代的發展,各個國家漸漸發現城市軌道交通在緩解交通擁堵問題中的積極作用,城市軌道交通作為一種必不可少的交通方式, 在解決交通擁擠、優化空間結構、改善自然環境、優化土地資源等方面具有顯著優勢,和它特有的便捷、可靠、安全、準時的功能特征, 在國內外發達城市中迅速發展起來。軌道交通信號設備是保證列車運行安全,實現行車指揮和列車運行現代化,提高運輸效率的關鍵系統設備[1]。地鐵車門系統安裝在列車上,是列車平穩運營、避免事故發生的保證,一旦其出現故障,輕則導致晚點、清客下線或抽線,重則造成危害乘客人身安全的重大事故。因此,地鐵車門系統故障預測至關重要。
故障樹分析通常應用于具有故障機制確定和故障邏輯關系的系統,因為它對系統的故障狀態做出了許多假設。故障樹中的事件只有兩種狀態:故障和正常[2]。 但是許多事件都是多態的,例如電機,除了故障和正常狀態。在額定速度以下還存在其他故障條件,并且故障樹難以描述具有多態性的事件。故障樹中的邏輯門描述了確定性邏輯關系,這要求上級事件和下一級事件之間存在明確的因果關系[3]。然而,對于許多復雜的系統,它們失敗的可能性很多,并且上層事件和下一個事件之間沒有明確的因果關系。在這種情況下,使用概率方法更合適,但邏輯門不具有描述概率的能力[4]。
相比故障樹分析,貝葉斯網絡對于故障的預測更具有準確性,貝葉斯網絡可以有效的解決故障樹刪除的內容:裝備理論與裝備技術二態性的問題。貝葉斯網絡是把故障的原因與故障表現作為表現內容,在故障征兆出現后可以依照節點之間的相關關系與概率值,計算得出各類設備產生故障的可能性,最后綜合處設備診斷工作的結果[5]。
故障樹分析(FTA)是一種分析方法,它將系統故障的原因從整體劃分為部分致密化。
FTA 是具有系統最不理想事件的分析(頂級事件)的目標。 確定系統中可能出現的組件之間的邏輯連接以及故障,環境變化,人為錯誤和其他因素(各種中間事件,基本事件等)和系統故障。
FTA 一般是基于如下假設的:
1) FT 中事件是二態的,即正常或故障;
2) 各個事件是彼此獨立的;
3) 故障和導致該故障的因果的關系是通過與門(AND)或者或門(OR)表達的,也有用其他邏輯門進行表達的,如非門(NOT)、異或門(XOR)、禁門等[6]。
從以上假設看,FT 存在一定的固有的局限性。
貝葉斯網絡是一種用圖表示知識的方法,并且是可以計算的概率模型[7]。貝葉斯網絡表示了變量之間的概率依賴關系,它主要由兩部分組成:有向無環圖;條件概率表[8]。
貝葉斯網絡用節點變量表達各個信息要素; 用連接節點之間的有向邊表達各個信息要素之間的關聯關系; 用條件概率表達各個信息要素之間的影響程度[9]。
推理是指根據已經發生或提供證據的事件做出一系列判斷的過程。貝葉斯網絡具有堅實的概率論基礎,因此有許多算法支持推理的發展。 每種算法在測量復雜性,應用普遍性和測量結果的準確性方面是不同的。貝葉斯網絡的推理一般分為兩個方向:精確推理和近似推理[10]。
聯合樹算法(Junction Tree,JT)相比其他的維耶斯推理算法,更加準確快速。
消息的傳遞過程可以分為兩個方向:信念收集(Collect Belief)方向與信念散播(Distribute Belief)方向。由于消息傳遞方案的不同,聯合樹算法可分為Shenoy-Shafer算法和Hugin算法[11]。其中,基于Hugin聯合樹推理算法分析故障貝葉斯網絡的具體過程如下:
1) 有向無環圖到聯合樹JT 的結構轉變。將有向圖轉換為無向圖:使用公共子節點連接每個父節點,并將所有有向邊更改為無向邊; 將無向圖三角化;將三角化的圖轉換為樹;找到該圖的根并查找最大的生成樹以獲得最終結果。
2) 聯合樹JT 的信念初始化。為了使聯合樹JT 滿足相應的約束,JT 中的所有簇(或“簇”)節點和相應的分區節點都被賦予信念潛力的初始值。 首先,將每個節點的置信勢φx的初始值設置為1,并且對于任何變量V,可以將對應的父節點記錄為Pa(V)。在包含V和Pa(V)的團(或“簇”)中,相應的條件概率為P(V|Pa(V)),則有φx←φxP(V|Pa(V))。該想法是通過使用故障貝葉斯網絡的條件概率來初始化聯合樹JT 的潛在功能。此后,在沒有給出證據時計算邊緣概率時,需要全局轉移和邊緣化操作;在計算給定證據的條件概率時,有必要引入觀察證據和歸一化操作。
如果JT滿足約束,則變量X的以下關系成立:
其中,Ci表示簇節點的信念勢,Sj表示分隔節點的信念勢。
3) JT的信念傳遞與吸收。為了使JT 一致,有必要在初始化之后執行與JT相反方向的置信傳輸。 選擇JT 中的根節點,實現信念在JT 中傳播和吸收的信念,從而使JT 的信念一致。
當信念從群集節點傳遞到分離節點時,存在
當信念從分離節點傳遞到群集節點時,存在
分析地鐵車門系統的故障信息,形成了由31個基礎故障事件,18個中間故障事件和一個頂部故障事件組成的轉向架故障系統[12]。
對各故障事件進行分類編號;同時,結合概率論和數理統計知識,首先在故障信息樣本的基礎上,計算每個基本故障事件的發生概率, 然后根據每個級別的中間事件的條件概率表達式確定。 具體見表1~表3。

表1 基礎故障事件信息

表2 中間故障事件信息

表3 故障頂事件信息
根據地鐵車門系統及其下級子系統劃分,分別建立對應的轉向架系統故障樹及其子系統故障樹,以驅動電機故障樹為例,圖1為驅動電機故障樹。

圖1 驅動電機故障樹示意圖
以地鐵門系統電機系統的故障為例,建立故障樹FT并將其轉換為故障貝葉斯網絡FBN。采用貝葉斯網絡的聯合樹推理方法對故障診斷進行分析。
在圖1 所示的故障樹中,M15視為頂事件,M1~M9為中間事件,X1、X6~X17、X31為底事件,邏輯門均為邏輯或門。其中:
X1表示EDCU故障X17表示電源切除開關誤動作
X6表示電機自身故障X31表示ATC信號設備故障
X7表示EDCU插頭松動M1表示信號產生失敗
X8表示EDCU程序版本低M2表示信號傳輸失敗
X9表示ATC信號設備故障M3表示無開門信號
X10表示開門按鈕接觸不良M4表示無開門使能信號
X11表示開門繼電器故障M5表示有關門信號
X12表示使能繼電器故障M6表示無零速信號
X13表示關門按鈕接觸故障M7表示信號錯誤
X14表示關門繼電器故障M8表示EDCU電源斷開
X15表示速度傳感器故障M9表示EDCU未起作用
X16表示速度傳感器故障M15表示驅動電機不動作
根據FT到FBN的轉換算法,轉換后轉換后的FBN如圖2所示。

圖2 驅動電機故障貝葉斯網絡
在確定貝葉斯網絡之后,可以基于觀察到的數據做出一系列假設。 這些假設改變了網絡中的先驗數據并逐漸影響整個網絡。 聯合樹算法將有向圖轉換為樹,這降低了計算難度。
將圖2貝葉斯網絡有向圖轉化為無向圖,并進行無向圖三角化。由于該無向圖不存在超過3個點的環。
確定簇節點,最終完成聯合樹的構建。在驅動電機聯合樹上進行信念傳遞,首先對其進行初始化。選取M18M9X6作為根節點,信念從子節點M1X9X10、M2X7X11、M4X12X31、M5X13X14、M6X15X16、M8X1X17傳遞到根節點M15M9X6,傳遞方向如圖3中φi所示,對應信念為φ~φ8;信念從根節點M15M9X6傳遞到各子節點M1X9X10、M2X7X11、M4X12X31、M5X13X14、M6X15X16、M8X1X17處,傳遞方向如圖3中Ψi所示,對應信念為Ψ1~Ψ8。

圖3 驅動電機聯合樹信念傳播示意圖
信念吸收的算式為
信念傳遞的算式為
通過在JT上進行信念初始化和信念傳遞,JT既滿足了約束性也達到了全局一致性。
以地鐵車門系統中VCU故障M16為例。可知,VCU故障中的緊急解鎖行程開關故障X2、鎖閉行程開關故障X18及關閉行程開關故障X19是輪對系統的3個主要故障模式,而且VCU 故障M16與3個節點故障是由邏輯或門相連的,即表示X2、X18、X19中任何一個故障的發生均有可能導致M16的發生,因此,可通過3個故障底事件的先驗概率分析和后驗概率計算來驗證對故障模式重要度的判斷,X2、X18、X19的先驗概率如表4所示,M16的條件概率如表5所示。

表4 基礎故障節點先驗概率

表5 VCU 故障節點的條件概率
首先,分析了3個故障底部事件的先驗概率。它可以在一定程度上反映每個故障底部事件對公共故障模式的影響程度,但受樣本采樣的限制。 并不能夠直接說明其重要程度。 可以看出,在3個底部事件中,緊急解鎖開關故障X2具有高達0.15的先驗概率。鎖定行程開關故障X18和閉合行程開關故障X19的先驗概率分別為0.10和0.06。兩者的先驗概率值相對接近,因此重要性等級不能僅由先驗概率完全確定。
測算結果如表6所示。

表6 基礎故障節點后驗概率
通過分析先驗概率和條件概率,可以看出,在計算故障后驗概率后,可以看出后驗概率的概率值較高。根據后驗概率表中的數據,相應事件對系統故障發生的影響越大,可以看出,輪對系統故障的3個底部事件的重要性順序為:鎖止行程開關故障X18,緊急解鎖行程開關故障X2,關閉行程開關故障X19。
將故障樹分析與貝葉斯網絡相結合,形成故障貝葉斯網絡,詳細說明了故障樹向故障貝葉斯網絡進行轉化的過程。以地鐵車門系統中基礎制動系統為例,進行基于聯接樹算法的貝葉斯網絡精確推理理論的實例分析。實現了故障貝葉斯網絡模型的Matlab仿真,進行了故障因果關系推理和故障重要度判斷分析。