基于蝴蝶結(jié)貝葉斯網(wǎng)絡(luò)的高校實(shí)驗(yàn)室火災(zāi)爆炸事故的風(fēng)險(xiǎn)評(píng)估

王 蕾，代養(yǎng)勇

（1.山東省食品藥品檢驗(yàn)研究院山東省食品藥品安全檢測(cè)工程技術(shù)研究中心，濟(jì)南 250101；2.山東農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院，山東泰安 271018）

0 引言

近些年，實(shí)驗(yàn)室事故時(shí)有出現(xiàn)，引發(fā)了社會(huì)的普遍關(guān)注。據(jù)相關(guān)媒體報(bào)道：2018 年北京交通大學(xué)東校區(qū)化學(xué)實(shí)驗(yàn)室爆炸，3 名學(xué)生死亡；2016 年上海東華大學(xué)環(huán)境與生物工程學(xué)院實(shí)驗(yàn)室發(fā)生爆炸，2 名學(xué)生重傷；2015 清華大學(xué)化學(xué)系實(shí)驗(yàn)室發(fā)生一起爆炸事故，1 名博士生當(dāng)場(chǎng)死亡。高?？蒲袑?shí)驗(yàn)室的安全事故頻發(fā)，尤其是化學(xué)實(shí)驗(yàn)室，已經(jīng)成為實(shí)驗(yàn)安全的“重災(zāi)區(qū)”。傳統(tǒng)的評(píng)價(jià)方法，如故障樹、失效模式和影響分析（FMEA）、危險(xiǎn)與可操作性分析（HAZOP）、保護(hù)層分析（LOPA），均為靜態(tài)評(píng)價(jià)方法只適用于評(píng)價(jià)某個(gè)時(shí)間點(diǎn)的風(fēng)險(xiǎn)，而無(wú)法得知整個(gè)實(shí)驗(yàn)室的風(fēng)險(xiǎn)動(dòng)態(tài)變化［1］。

目前對(duì)實(shí)驗(yàn)室安全風(fēng)險(xiǎn)評(píng)價(jià)的研究多集中在以職業(yè)健康管理體系要求進(jìn)行危險(xiǎn)源識(shí)別，然后根據(jù)歷年的數(shù)據(jù)和發(fā)生的概率人為附加分值進(jìn)行判斷。劉婷婷等［2］《重大危險(xiǎn)源辨識(shí)》（GB18218—2014）和《職業(yè)性接觸毒物危害程度分析》（GB50844—1985），重點(diǎn)進(jìn)行危險(xiǎn)源的識(shí)別；魏永前等［3］構(gòu)建實(shí)驗(yàn)室安全“雙體系”預(yù)防機(jī)制從制度建設(shè)和專項(xiàng)整治方面進(jìn)行控制，真正運(yùn)用風(fēng)險(xiǎn)評(píng)估模型進(jìn)行風(fēng)險(xiǎn)評(píng)估的鳳毛麟角。因此本文引入蝴蝶結(jié)模型（Bow-Tie model，下文簡(jiǎn)稱為BT）和貝葉斯網(wǎng)絡(luò)（Bayesian Network，下文簡(jiǎn)稱為BN）模型，通過條件概率對(duì)各事件的依賴程度進(jìn)行分析評(píng)價(jià)，利用貝葉斯網(wǎng)絡(luò)軟件進(jìn)行反向推理，對(duì)影響實(shí)驗(yàn)室火災(zāi)爆炸事故的因素進(jìn)行分析評(píng)價(jià)，實(shí)現(xiàn)對(duì)高校實(shí)驗(yàn)室火災(zāi)爆炸事故的風(fēng)險(xiǎn)預(yù)測(cè)和故障診斷的雙重功能。

1 蝴蝶結(jié)模型和貝葉斯網(wǎng)絡(luò)模型構(gòu)建

1.1 蝴蝶結(jié)（BT）模型

BT模型是表述一個(gè)事故場(chǎng)景最好的圖形化方法之一，它可以清楚地描述事故的原因和事故的結(jié)果，是將故障樹與事件樹分析方法融合為一體的事故因果分析方法。由圖1 可知，BT 模型左側(cè)是故障樹，分析有哪些危害，危害一旦釋放會(huì)造成哪些威脅，其安全屏障的作用是預(yù)防事故；右側(cè)是事件樹，表示一旦屏障和補(bǔ)救措施失效后會(huì)產(chǎn)生什么樣嚴(yán)重的后果［3］。

1.2 貝葉斯網(wǎng)絡(luò)模型（BN模型）

貝葉斯理論是處理不確定性信息的重要工具，是一種基于概率的不確定性推理方法，是目前處理非確定性問題的一種主流方法。BN 網(wǎng)絡(luò)是有向無(wú)環(huán)圖，由代表變量節(jié)點(diǎn)及連接這些節(jié)點(diǎn)有向邊構(gòu)成，用條件概率進(jìn)行表達(dá)關(guān)系強(qiáng)度，節(jié)點(diǎn)變量可以是任何問題的抽象。BN 的理論基礎(chǔ)即貝葉斯條件概率計(jì)算公式［4-7］，即

式中，P（Vi|A）為后驗(yàn)概率，i=1，2，…，n；P（A|Vi）為事件Vi發(fā)生條件下事件A 發(fā)生的概率，即條件概率；P（Vi）為先驗(yàn)概率；P（A）為事件A發(fā)生的概率。

根據(jù)貝葉斯公式結(jié)合蝴蝶結(jié)模型中分析出的要素轉(zhuǎn)換成的各個(gè)節(jié)點(diǎn)以及節(jié)點(diǎn)之間的依賴關(guān)系，其各節(jié)點(diǎn)的聯(lián)合概率為P（V），見式（2）。

式中，Vpa（vi）為節(jié)點(diǎn)Vi的所有父節(jié)點(diǎn)。

1.3 風(fēng)險(xiǎn)評(píng)估模型的建立

以近10 年來(lái)高校實(shí)驗(yàn)室發(fā)生的112 起典型事故為數(shù)據(jù)來(lái)源，對(duì)引起高校實(shí)驗(yàn)室火災(zāi)爆炸事故的因素進(jìn)行整理匯總分析，根據(jù)BT 模型分析識(shí)別出關(guān)鍵風(fēng)險(xiǎn)源，并對(duì)關(guān)鍵風(fēng)險(xiǎn)源進(jìn)行同類合并轉(zhuǎn)換為BN 模型中的節(jié)點(diǎn)，再根據(jù)BN模型簡(jiǎn)化相同事件，尋找出造成頂端事件的主要節(jié)點(diǎn)以及各種節(jié)點(diǎn)之間的關(guān)系。根據(jù)歷史數(shù)據(jù)以及專家經(jīng)驗(yàn)，確定實(shí)驗(yàn)室安全事故的先驗(yàn)概率，運(yùn)用公式進(jìn)行動(dòng)態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。

賀蕾等［8］分析出高校實(shí)驗(yàn)室事故主要為火災(zāi)性、爆炸性事故，占全部安全事故的81%。造成火災(zāi)、爆炸事故的主要危險(xiǎn)源有：人為因素、設(shè)備因素、實(shí)驗(yàn)設(shè)施以及管理措施不健全等［9］。①人為因素，主要是指違規(guī)操作、不能正確使用試劑試藥等易揮發(fā)易爆炸的化學(xué)藥品。②設(shè)備因素，尤其是壓力容器，如氫氣瓶、氧氣瓶、乙炔氣瓶等。壓力容器爆炸，不僅會(huì)造成人員傷亡事故，內(nèi)部介質(zhì)擴(kuò)散還會(huì)引發(fā)燃燒等連鎖反應(yīng)。③設(shè)施因素，部分實(shí)驗(yàn)室缺乏必備的消防設(shè)施。有些高校的化學(xué)實(shí)驗(yàn)室，尤其是年代久遠(yuǎn)的教學(xué)樓，缺乏必備火災(zāi)報(bào)警、煙霧噴淋系統(tǒng)等消防設(shè)施。④管理因素，缺乏對(duì)實(shí)驗(yàn)室人員的安全教育和技能培訓(xùn)，相關(guān)的管理制度及應(yīng)急措施不健全等。

上述危險(xiǎn)源依據(jù)BT模型［10］可以更好地識(shí)別高校實(shí)驗(yàn)室發(fā)生火災(zāi)爆炸的主要原因、采取措施的有效性以及可能產(chǎn)生的后果等因素，具體模型見圖2 所示，其中C1～C7是代表所有控制措施失效后產(chǎn)生的后果，本文重點(diǎn)針對(duì)造成火災(zāi)爆炸事故的原因進(jìn)行分析，因此不再對(duì)產(chǎn)生的后果進(jìn)行過多闡述。根據(jù)BT 模型的分析，可以很直觀的得出影響高校實(shí)驗(yàn)室發(fā)生火災(zāi)爆炸事故的主要變量，如表1 所示。

圖2 實(shí)驗(yàn)室火災(zāi)、爆炸BT模型

表1 實(shí)驗(yàn)室火災(zāi)、爆炸事故主要變量

BT模型向BN模型轉(zhuǎn)化主要是將BT模型中的事故樹和事件樹分別轉(zhuǎn)化為貝葉斯網(wǎng)絡(luò)模型中的節(jié)點(diǎn)，并使用有向弧把各個(gè)節(jié)點(diǎn)連接起來(lái)建立映射關(guān)系［11］。具體轉(zhuǎn)化見圖3。

圖3 BT模型向BN模型轉(zhuǎn)化示意圖

考慮到BN推理案例中沒有高校實(shí)驗(yàn)室火災(zāi)爆炸事故的案例，本文選擇基于專家知識(shí)以及部分文獻(xiàn)數(shù)據(jù)構(gòu)建BN模型。根據(jù)眾多文獻(xiàn)以及媒體報(bào)道引起高校實(shí)驗(yàn)室火災(zāi)爆炸的原因［12］，對(duì)BN 變量進(jìn)行優(yōu)化，對(duì)安全制度、應(yīng)急預(yù)案、違規(guī)操作、壓力容器反應(yīng)釜陳舊、設(shè)備長(zhǎng)時(shí)間運(yùn)轉(zhuǎn)、電路老化、消防設(shè)施陳舊等建立模型，其他變量引起事故的比例較少暫不建模［13］，優(yōu)化后建立的模型見圖4。

圖4 BN網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)圖

根據(jù)近10 年消防年鑒中相關(guān)數(shù)據(jù)以及專家知識(shí)可以得出各父節(jié)點(diǎn)的參考權(quán)重，其中管理制度的權(quán)重為2.7%，應(yīng)急預(yù)案的權(quán)重為1.8%，違規(guī)操作的權(quán)重為16.1%，易燃易爆化學(xué)品使用的權(quán)重為23.2%，壓力容器反應(yīng)釜陳舊的權(quán)重為8.9%，設(shè)備長(zhǎng)時(shí)間運(yùn)行的權(quán)重為16.9%，電路老化的權(quán)重為21.4%，消防設(shè)施陳舊的權(quán)重為8.9%。

BN網(wǎng)絡(luò)實(shí)質(zhì)上體現(xiàn)了各節(jié)點(diǎn)變量間的概率依賴關(guān)系，網(wǎng)絡(luò)中每一個(gè)節(jié)點(diǎn)都對(duì)應(yīng)一個(gè)條件概率表（Conditional Probability Table，CPT）［14］，表示該節(jié)點(diǎn)與其父節(jié)點(diǎn)間的概率分布關(guān)系。鑒于變量節(jié)點(diǎn)較多，選取節(jié)點(diǎn)H為例進(jìn)行說明，由此可得，具體如表2 所示。其中，“NO”代表不安全的概率，“YES”代表安全的概率。

表2 人為因素的條件概率（CPT）

根據(jù)圖3 模型中各節(jié)點(diǎn)的依賴關(guān)系，運(yùn)用Netica軟件構(gòu)建BN模型，并進(jìn)行計(jì)算，具體如圖5 所示。根據(jù)已經(jīng)發(fā)生的高校實(shí)驗(yàn)室安全事故中各個(gè)節(jié)點(diǎn)發(fā)生的概率，利用貝葉斯公式綜合計(jì)算，得到全國(guó)高校實(shí)驗(yàn)室火災(zāi)爆炸安全事故的發(fā)生概率為0.175，較全國(guó)整體火災(zāi)、爆炸發(fā)生比例高很多。由此可以看出，高校實(shí)驗(yàn)室存在很大的安全隱患，BN 模型的構(gòu)建可以為預(yù)警高校實(shí)驗(yàn)室火災(zāi)、爆炸事故提供一個(gè)動(dòng)態(tài)的分析模型，學(xué)校安全部門可以將日常安全檢查的數(shù)據(jù)輸入模型，進(jìn)行風(fēng)險(xiǎn)預(yù)警，便于采取相應(yīng)的措施。

圖5 BN模型以及各節(jié)點(diǎn)發(fā)生概率

2 結(jié)果驗(yàn)證

2.1 模型的驗(yàn)證

以2001～2016 年高校實(shí)驗(yàn)室火災(zāi)爆炸發(fā)生情況為例，2001～2016 年我國(guó)高校實(shí)驗(yàn)室共發(fā)生火災(zāi)爆炸事故93 起，從事故原因分析中可以得出易燃易爆化學(xué)品導(dǎo)致的火災(zāi)爆炸事故37.7%，電路老化5.7%，違規(guī)操作20.8%，壓力容器反應(yīng)釜陳舊22.6%，設(shè)備長(zhǎng)時(shí)間運(yùn)行8.4%，消防設(shè)施陳舊3.8%，管理制度0.8%，應(yīng)急預(yù)案0.2%［8］。按照易燃易爆化學(xué)品導(dǎo)致發(fā)生事故的比率降為5%進(jìn)行計(jì)算，將相關(guān)數(shù)據(jù)輸入模型，通過模型計(jì)算得到火災(zāi)爆炸事故的發(fā)生概率為0.159，比之前的0.175 有所下降，證明該模型有效。

2.2 節(jié)點(diǎn)后驗(yàn)概率分布

后驗(yàn)概率分析是BN 網(wǎng)絡(luò)推理中涉及最多的問題，一方面可在已知原因前提下推斷出結(jié)果發(fā)生的概率，即預(yù)測(cè)推理，同時(shí)也可以在已知結(jié)果的狀態(tài)下，找出導(dǎo)致結(jié)果發(fā)生的最可能原因，也就是所說的故障診斷［15］。本文運(yùn)用貝葉斯網(wǎng)絡(luò)的故障診斷功能，找出導(dǎo)致目標(biāo)事件發(fā)生的最可能因素，為更準(zhǔn)確的評(píng)估事件態(tài)勢(shì)做出鋪墊和提供依據(jù)。

在貝葉斯網(wǎng)絡(luò)模型中假設(shè)節(jié)點(diǎn)E 發(fā)生的概率為100%（也就是說實(shí)驗(yàn)室火災(zāi)爆炸發(fā)生的概率為100%），得出節(jié)點(diǎn)的后驗(yàn)證概率分布，具體如圖6所示。

圖6 后節(jié)點(diǎn)概率示意圖（注：YES 代表節(jié)點(diǎn)安全，NO 代表節(jié)點(diǎn)不安全）

從圖6 可以看出，當(dāng)實(shí)驗(yàn)室發(fā)生火災(zāi)爆炸事故的情況下，易燃易爆化學(xué)品使用（H4）、電路老化（B2）發(fā)生的可能性最大，概率超過40%；其次為設(shè)備長(zhǎng)時(shí)間運(yùn)行（D2）概率為30.5%、違規(guī)操作（H1）概率為27.9%；最不可能的原因?yàn)閼?yīng)急預(yù)案和管理制度失效。

2.3 模型敏感因素分析

在BN網(wǎng)絡(luò)中，子節(jié)點(diǎn)的后驗(yàn)概率會(huì)隨父節(jié)點(diǎn)的先驗(yàn)概率的變化而改變，敏感性分析將由父節(jié)點(diǎn)發(fā)生變化而引起子節(jié)點(diǎn)的參數(shù)變化程度進(jìn)行了量化，從而識(shí)別出BN模型中的關(guān)鍵因素［16］。在GeNie軟件中將根節(jié)點(diǎn)E設(shè)為目標(biāo)節(jié)點(diǎn)，執(zhí)行敏感分析，得出節(jié)點(diǎn)敏感度分析，如下圖所示，其中，節(jié)點(diǎn)的顏色深度與敏感度成正比。

由圖7 中可知，整個(gè)模型中節(jié)點(diǎn)的敏感程度可分為6 個(gè)等級(jí)，其中，高敏感為節(jié)點(diǎn)H（人員不安全行為）；較高敏感度為H4（易燃易爆化學(xué)品的使用）、B2（電路老化）、B（設(shè)施陳舊）。通過敏感性分析可以發(fā)現(xiàn)：人員不安全行為、易燃易爆化學(xué)品的使用、電路老化、設(shè)施陳舊是導(dǎo)致實(shí)驗(yàn)室發(fā)生火災(zāi)爆炸事故的主要影響因素。

圖7 BN模型變量敏感性分析

3 結(jié)論

（1）BT-BN模型兼顧了事故樹、事件樹以及貝葉斯網(wǎng)絡(luò)定量分析風(fēng)險(xiǎn)的優(yōu)點(diǎn)，可以對(duì)高?；馂?zāi)、爆炸事故的風(fēng)險(xiǎn)進(jìn)行動(dòng)態(tài)分析，克服其他體系或者檢查過程中大量零散數(shù)據(jù)的缺點(diǎn)。將日常的安全巡查、人員培訓(xùn)等相關(guān)數(shù)據(jù)信息匯總后輸入模型，可得到安全提示，以便及時(shí)采取控制措施防止發(fā)生嚴(yán)重后果。

（2）從構(gòu)建的模型中正向計(jì)算，可知易燃易爆化學(xué)品的使用是誘發(fā)高校實(shí)驗(yàn)室發(fā)生火災(zāi)爆炸事故的主要誘因。通過模型故障診斷可知，最可能導(dǎo)致事故發(fā)生的因素是易燃易爆化學(xué)品使用、線路老化。

（3）要有效避免高校實(shí)驗(yàn)室火災(zāi)、爆炸事故，就必須加強(qiáng)實(shí)驗(yàn)人員培訓(xùn)提高安全意識(shí)，完善易燃易爆化學(xué)品存儲(chǔ)條件以及化學(xué)性質(zhì)的標(biāo)志標(biāo)識(shí)，同時(shí)還要定期對(duì)實(shí)驗(yàn)室線路進(jìn)行檢查和維護(hù)。