初珠立 楊自春

海軍工程大學(xué) 船舶與動力學(xué)院，湖北 武 漢 430033

船用鍋爐主過熱蒸汽管路破損的貝葉斯網(wǎng)絡(luò)安全性定量分析

初珠立 楊自春

海軍工程大學(xué) 船舶與動力學(xué)院，湖北 武 漢 430033

針對貝葉斯網(wǎng)絡(luò)方法在處理大規(guī)模復(fù)雜系統(tǒng)上的優(yōu)點，結(jié)合當(dāng)前船用鍋爐的安全性定量分析特點，提出將貝葉斯網(wǎng)絡(luò)方法應(yīng)用于船用鍋爐的安全性分析。通過探討貝葉斯網(wǎng)絡(luò)及故障樹分析方法原理，嘗試將貝葉斯網(wǎng)絡(luò)與故障樹分析方法結(jié)合起來，主要突出貝葉斯網(wǎng)絡(luò)方法可以實現(xiàn)數(shù)據(jù)的自動更新及由后驗信息到先驗信息的反向推理能力。最后將該方法應(yīng)用于船用鍋爐的安全性定量分析中，通過更新觀測部件故障狀態(tài)實現(xiàn)了整個系統(tǒng)的自動更新，并自動查找出最大概率事件，證實該方法提高了安全性分析效率，使實際操作更加準(zhǔn)確、快捷，極大地的降低了事故帶來的風(fēng)險。

貝葉斯網(wǎng)絡(luò)；安全性定量分析；反向推理；船用鍋爐

1 引言

憑借其良好的特性，貝葉斯網(wǎng)絡(luò)（BN）已經(jīng)越來越廣泛和深入地應(yīng)用于安全性與可靠性領(lǐng)域。Sankaran等人利用BN對復(fù)雜系統(tǒng)中順序失效和相關(guān)失效進行描述［1］，提高了安全性分析能力；同時，還研究了在獲得實際數(shù)據(jù)條件下，通過BN的證據(jù)傳播特性來檢驗安全性模型的有效性［2］。Bobbioa等人研究了將故障樹（FT）映射成BN的方法，并利用BN對依靠系統(tǒng)進行可靠性分析［3］。

在船用鍋爐安全性定量分析中廣泛存在不確定性和復(fù)雜關(guān)聯(lián)性問題，而目前的分析方法諸如故障樹分析（FTA），故障模式、影響及危害性分析，事件樹分析，概率風(fēng)險分析等方法又有一定的局限性。從推理機制和系統(tǒng)狀態(tài)描述上來看，近年來發(fā)展起來的貝葉斯網(wǎng)絡(luò)技術(shù)與故障樹有很大的相似性。而且貝葉斯網(wǎng)絡(luò)還具有描述事件多態(tài)性和故障邏輯關(guān)系非確定性的能力，更加適合于對復(fù)雜系統(tǒng)的安全性進行分析。本文試圖采用貝葉斯網(wǎng)絡(luò)技術(shù)對故障樹進行分析，利用故障樹的直觀性同時結(jié)合貝葉斯網(wǎng)絡(luò)對故障事件較強的邏輯描述能力，對船用鍋爐的安全性進行定量分析。

2 貝葉斯網(wǎng)絡(luò)原理簡述

貝葉斯網(wǎng)絡(luò)是一種有向無環(huán)圖，它的具體定義如下［4］。 設(shè) U ＝ ｛ x1，……，xn｝，n≥1 是一個由變量組成的集合。設(shè)Bs是一個建立在變量集U上的網(wǎng)絡(luò)結(jié)構(gòu)，且Bs是一個有向無環(huán)圖（DAG）；設(shè)Bp是U中變量所組成的條件概率，且有Bp＝｛p（u是Bs中u的雙親集合。那么，由Bs和Bp共同組成的網(wǎng)絡(luò)就叫做貝葉斯網(wǎng)絡(luò)，這個網(wǎng)絡(luò)所表示的聯(lián)合概率概率分布為

3 故障樹分析方法

FTA［5］是一種自上而下的分析技術(shù)，它對導(dǎo)致不希望事件發(fā)生的故障進行并行和有序的綜合分析，研究每種故障的原因，以給出不希望事件在給定環(huán)境條件下可能發(fā)生的概率，利用該方法既可以進行定性分析也可以進行定量的分析［6］。在定量分析中，F(xiàn)TA可給出各產(chǎn)品層次單個零件所需的可靠性數(shù)據(jù)，并用于確定故障概率和風(fēng)險評價。FTA的定量分析是定性分析的繼續(xù)和發(fā)展。隨著工程項目研制的深入，系統(tǒng)的細節(jié)信息逐步豐富。定量分析正是在獲取故障樹中個底事件的發(fā)生概率后，在所有底事件相互獨立的條件下，求出頂事件發(fā)生概率以及確定每個底事件的重要度。但是相較于貝葉斯網(wǎng)絡(luò)而言，F(xiàn)TA方法將所有的事件定義為兩種狀態(tài)：即正常和失效（分別用0和1來表示），這就限制了其應(yīng)用范圍，且對于多態(tài)事件的處理將不再適合或者有較大的誤差。

4 基于FTA的貝葉斯網(wǎng)絡(luò)構(gòu)建

相關(guān)文獻［7-9］闡述了故障樹與貝葉斯網(wǎng)絡(luò)之間的網(wǎng)絡(luò)拓撲結(jié)構(gòu)不變的關(guān)系，介紹了故障樹與貝葉斯網(wǎng)絡(luò)之間的拓撲映射關(guān)系。圖1所示為與、或、非以及2／3表決邏輯關(guān)系向貝葉斯網(wǎng)絡(luò)的轉(zhuǎn)化示例。 另外，Helge Langseth 等［11］將其方法進行了歸納和總結(jié)，給出了故障樹的各種門向貝葉斯網(wǎng)絡(luò)轉(zhuǎn)化確定方法，具體過程如下：

步驟1：故障樹中的每個底事件對應(yīng)貝葉斯網(wǎng)絡(luò)中一個節(jié)點，重復(fù)事件只建立一個節(jié)點。

步驟2：根據(jù)故障樹中各個底事件的已知概率來確定貝葉斯網(wǎng)絡(luò)中相應(yīng)節(jié)點的概率。

步驟3：故障樹中節(jié)點間的邏輯關(guān)系，可通過相應(yīng)貝葉斯網(wǎng)絡(luò)中添加條件概率表來表示。

步驟4：按照故障樹中頂事件，中間事件和底事件之間的連接關(guān)系建立貝葉斯網(wǎng)絡(luò)中節(jié)點之間的連接。

步驟5：按照故障樹中與、或、非及表決邏輯關(guān)系確定貝葉斯網(wǎng)絡(luò)中非根節(jié)點條件概率分布。

由于貝葉斯網(wǎng)絡(luò)已經(jīng)有很成熟的算法計算節(jié)點的聯(lián)合概率分布和在各種證據(jù)下的條件概率分布［10］，因而在構(gòu)建了系統(tǒng)的貝葉斯網(wǎng)絡(luò)之后，就可以很方便地進行概率安全評估，包括計算各個后果發(fā)生的概率、各底事件的重要度以及其他信息［11］。

5 實例：船用鍋爐的安全性分析

5.1 故障樹的建立

作為船舶動力系統(tǒng)的核心部件，鍋爐起著舉足輕重的作用。在鍋爐的日常使用和維護過程中，會出現(xiàn)很多復(fù)雜的故障，比如爐膛隔離功能減弱或喪失、鍋爐煙道絕熱保溫性能下降或喪失、爐部活動底座自由膨脹受阻、鍋爐鍋部漏泄、安全閥主閥打不開、內(nèi)爐板過熱變形、主（弱）過熱蒸汽管束破損、管子變形下垂、鼓包和產(chǎn)氣量不足等，其中主過熱蒸汽管束破損發(fā)生頻率較高，并且與安全性聯(lián)系較為緊密，現(xiàn)將這一事故作為故障事件，進行安全性分析，建立如圖2所示故障樹。圖2中各事件編號及含義（Q值為各故障事件的概率值）見表1和表2。其中M2（蒸汽過熱）故障下面各事件關(guān)系為X3（高溫氧化）、X4（壓力調(diào)節(jié)器故障）與M3（工況急劇變化）事件為邏輯與關(guān)系，同時與M4（燃油霧化不良造成偏燒）構(gòu)成邏輯或關(guān)系，這一點將在貝葉斯網(wǎng)絡(luò)的邏輯關(guān)系中反應(yīng)出來。表2中故障樹底事件的定量信息數(shù)據(jù)是在部分參考文獻［12］的基礎(chǔ)上，查閱相關(guān)設(shè)備運行信息和可靠性信息等技術(shù)資料，同時結(jié)合操作人員和專家意見的基礎(chǔ)上給出的。

5.2 故障樹向貝葉斯網(wǎng)絡(luò)的轉(zhuǎn)化及整合

根據(jù)上述原理，將故障樹轉(zhuǎn)化為貝葉斯網(wǎng)絡(luò)，如圖3所示。

貝葉斯網(wǎng)絡(luò)自動分析結(jié)果如圖4所示。貝葉斯網(wǎng)絡(luò)較故障樹分析的優(yōu)點之一體現(xiàn)在部分數(shù)據(jù)更改后整個系統(tǒng)會自動更新并計算出結(jié)果。如上例中將節(jié)點M2（蒸汽過熱）更新后系統(tǒng)將自動更新，結(jié)果如圖5所示，將更新前后的相關(guān)節(jié)點值列入表3進行分析 （考慮到軟件BayesiaLab的計算精確性及各節(jié)點概率值較小這一點，將各Q值經(jīng)過處理后輸入軟件進行計算；未列出節(jié)點為非相關(guān)節(jié)點，數(shù)據(jù)無變化）。從表中可以直觀的看出當(dāng)節(jié)點 M2更新后， 節(jié)點 X1，X3，X4，X5，X6，X7，X8，X9，X10，X11，X12，X13的值都有改變。 其中 X10（回油閥磨損過大）的值改變值最為顯著，其次是X7（噴油器結(jié)焦）、X9（燃油總管破損）、X12（汽輪燃油泵排量降低），可以理解為引起M2（蒸汽過熱）事件發(fā)生的最可能因素是回油閥磨損過大，其次可能是噴油器結(jié)焦、燃油總管破損、汽輪燃油泵排量降低。在實際操作中如果出現(xiàn)蒸汽過熱現(xiàn)象，應(yīng)該優(yōu)先檢測上述事件。

表1 底事件表

表2 結(jié)果事件表

表3 節(jié)點M2更新后相關(guān)節(jié)點的變化情況

5 結(jié)束語

利用故障樹分析方法與貝葉斯網(wǎng)絡(luò)的對應(yīng)關(guān)系，建立故障樹分析圖并向貝葉斯網(wǎng)絡(luò)轉(zhuǎn)化，最后結(jié)合船用鍋爐主過熱蒸汽管束破損這一故障進行分析，綜合先驗信息及后驗信息進行評定，能夠很好地確定最大概率事件，為船用鍋爐的操作，維修帶來了較大的改進，在事故發(fā)生時能更快的確定誘發(fā)事件，為及時正確地采取安全性防范措施提供了保證。該方法驗證了貝葉斯網(wǎng)絡(luò)方法是合理的、完備的，省略了繁雜的最小割集計算，可以自動計算出未知節(jié)點的概率值，并通過數(shù)據(jù)模擬計算演示貝葉斯網(wǎng)絡(luò)在后驗概率推理方面的強大優(yōu)勢，即系統(tǒng)可以根據(jù)數(shù)據(jù)的變動進行反向推理并實現(xiàn)自動更新，相較FTA方法更為先進，其在裝備的安全性定量分析中將有越來越廣泛的應(yīng)用。

［1］SANKARAN M， ZHANG R X， NATASHA S.Bayesian networks for system reliability reassessment ［J］.Structural Safety，2001，23（6）：231－251.

［2］SANKARAN M，RAMESH R.Validation of reliability computational models using bayes networks［J］.Reliability Engineering and System Safety，2005，87（10）：223－232.

［3］BOBBIOA A，PORTINALEA L，MINICHINOB M，et al.Improving the analysis of dependable systems by mapping fault trees into Bayesian networks ［J］.Reliability Engineering and System Safety， 2001，71（12）：249－260.

［4］廉捷.貝葉斯網(wǎng)絡(luò)構(gòu)造方法及應(yīng)用研究［D］.北京：北京交通大學(xué)，2007.

［5］周忠寶，周經(jīng)倫，金光，等.基于貝葉斯網(wǎng)絡(luò)的概率安全評估方法研究 ［J］.系統(tǒng)工程學(xué)報，2006，12 （21）：636-643，667.

［6］董豆豆.基于故障樹的系統(tǒng)安全風(fēng)險實時監(jiān)測方法［J］.國防科技大學(xué)學(xué)報，2006，9（28）：112.

［7］BOBBIO A，PORTINALE L，MINICHINO M，et al.Comparing fault trees and bayesian networks for dependability analysis［C］//The 18th International Conferenceon Computer Safety，Reliability and Security，Toulouse，F(xiàn)rance：1999：310－322.

［8］王廣彥，馬志軍，胡起偉.基于貝葉斯網(wǎng)絡(luò)的故障樹分析［J］.系統(tǒng)工程理論與實踐，2004，24（6）：76－83.

［9］BOBBIO A，PORTINALE L，MINICHINO M，et al.Improving the analysis of dependable systems by mapping fault trees into bayesian networks ［J］.Reliability Engineering and System Safety，2001，71（3）：249－260.

［10］謝斌，張明珠，嚴(yán)于鮮.貝葉斯網(wǎng)絡(luò)對故障樹方法的改進［J］.綿陽師范學(xué)院學(xué)報，2004，23（2）：29－33.

［11］HELGE L，LUIGI P.Bayesian network in reliability［J］.Reliability Engineering and System Safety，2007，2 （1）：92－108.

［12］林潔.GO－FLOW 原理及其計算機輔助技術(shù)［D］.長沙：國防科學(xué)技術(shù)大學(xué)，2003.

Quantitative Safety Analysis on the Damage of Marine Boiler Steam Gas Pipelines Based on Bayesian Network

Chu Zhu－li Yang Zi－chun
College of Naval Architecture and Power，Naval University of Engineering，Wuhan 430033，China

Bayesian Network （BN） has advantages in dealing with the large scale complex system and can be used to undertake safety analysis of marine boiler system.Attempts have been made to conduct such analysis by combining the BN theory with the Fault Tree Analysis （FTA） method.In this study，particular attention was given to the automatic updating ability of the BN method as well as its capability of reverse reasoning from posteriori information to apriori information.Using this methods， the whole system was updated automatically by upgrading the FT status of detection units and therefore the events of maximum probability were determined.The results indicate that using this method can make the safety assessment more accurately and rapidly and also have the potential to improve the analysis＇s efficiency and greatly mitigate the risks brought by the accidents.

Bayesian network；quantitative safety analysis； reverse reasoning；marine boiler

TK229.7

A

1673－3185（2010）05－59－05

10.3969/j.issn.1673-3185.2010.05.012

2009－08-24

××裝備預(yù)研基金項目（9140A27050106JB11）

初珠立（1985－），男，碩士研究生。研究方向：艦船動力及熱力系統(tǒng)的科學(xué)管理。E-mail：hgchuzhuli＠163.com

楊自春（1968－），男，教授，博士生導(dǎo)師。研究方向：熱力系統(tǒng)的可靠性與故障診斷