基于風(fēng)險(xiǎn)耦合的石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)評(píng)價(jià)

張伯瑩(中海石油技術(shù)檢測(cè)有限公司，天津 300452)

0 引言

石油儲(chǔ)罐因其中儲(chǔ)存的介質(zhì)大都具有易燃、易爆特性而易發(fā)生火災(zāi)爆炸事故，導(dǎo)致人員傷亡及財(cái)產(chǎn)損失。石油罐區(qū)風(fēng)險(xiǎn)管理的研究引起國(guó)內(nèi)外學(xué)者的重視。Moshashaei等研究并總結(jié)了浮頂油罐火災(zāi)爆炸事故原因，確定了11個(gè)主要因素及71個(gè)子因素。Halloul等采用模糊故障樹分析與專家訪談相結(jié)合的方法對(duì)原油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)進(jìn)行識(shí)別與評(píng)估。C.D.Argyropoulos提出了采用檢查表法對(duì)儲(chǔ)罐進(jìn)行風(fēng)險(xiǎn)識(shí)別。Guo Xiaofang等則預(yù)測(cè)和評(píng)價(jià)了油罐爆炸引起的熱輻射和沖擊波的影響。目前的研究尚缺少對(duì)于風(fēng)險(xiǎn)因素耦合的關(guān)注，而事故的發(fā)生多為不同風(fēng)險(xiǎn)因素間相互作用的結(jié)果。因此，本文從風(fēng)險(xiǎn)因素耦合的角度出發(fā)，評(píng)價(jià)石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)因素，并研究不同風(fēng)險(xiǎn)因素之間的耦合度，希望能夠?yàn)楣迏^(qū)安全管理提供參考。

1 建立風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系

1.1 石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)因素辨識(shí)

通過(guò)相關(guān)文獻(xiàn)的查閱，將石油罐區(qū)火災(zāi)爆事故風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)分為人的風(fēng)險(xiǎn)因素、設(shè)備風(fēng)險(xiǎn)因素、環(huán)境風(fēng)險(xiǎn)因素、管理風(fēng)險(xiǎn)因素、消防風(fēng)險(xiǎn)因素5個(gè)一級(jí)風(fēng)險(xiǎn)指標(biāo)，具體風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系如圖1所示。

圖1 石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系

1.2 基于AHP的指標(biāo)權(quán)重計(jì)算

層次分析法是一種定性與定量相結(jié)合的評(píng)價(jià)方法，在安全、環(huán)境等多個(gè)領(lǐng)域有著廣泛的應(yīng)用。

1.2.1 構(gòu)造判斷矩陣

通過(guò)行業(yè)內(nèi)若干名專家對(duì)各層指標(biāo)兩兩比對(duì)，以1～9及其倒數(shù)作為相對(duì)重要性比例標(biāo)度構(gòu)造判斷矩陣A，標(biāo)度方法見表1。

表1 1～9標(biāo)度方法

1.2.2 計(jì)算權(quán)重向量W

運(yùn)用特征向量法求得權(quán)重向量W，將權(quán)重向量右乘判斷矩陣A，有：

式中：λmax為判斷矩陣A的最大特征值；λmax存在且唯一，將求得的W進(jìn)行歸一化處理。

1.2.3 層次單排序及一致性檢驗(yàn)

計(jì)算一致性指標(biāo)CI，有：

計(jì)算一致性比例CR，有：

式中：RI取值見表2；若CR＜0.1，則判斷矩陣滿足一致性檢驗(yàn)，若CR≥0.1，則判斷矩陣不滿足一致性檢驗(yàn)，需對(duì)判斷矩陣作適當(dāng)修改。

表2 平均隨機(jī)一致性指標(biāo)

2 基于逆向云模型的打分及期望

云模型是由李德毅院士于1995年提出的一種使定性概念與定量描述之間發(fā)生轉(zhuǎn)換的模型，由期望Ex，熵En以及超熵He表述[6]。云發(fā)生器即實(shí)現(xiàn)定性概念與定量描述相互轉(zhuǎn)換發(fā)生的過(guò)程，包括正向云、逆向云、條件云三種。

逆向云發(fā)生器可將定性的概念量化，將專家打分情況轉(zhuǎn)換為數(shù)字特征值Ex、En、He，同時(shí)可通過(guò)數(shù)字特征剔除離散程度大的打分，保證打分的準(zhǔn)確性。云模型特征公式如(4)～(6)所示。

由此可實(shí)現(xiàn)對(duì)底層指標(biāo)的打分，并得出相應(yīng)期望值。

3 建立風(fēng)險(xiǎn)耦合度模型

3.1 構(gòu)造功效函數(shù)

針對(duì)石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)系統(tǒng)，其耦合度指標(biāo)的功效函數(shù)可表示為：

式中：Uij為一級(jí)底層指標(biāo)對(duì)其上級(jí)子系統(tǒng)的有序度的功效，表示各項(xiàng)指標(biāo)與目標(biāo)的一致性，Uij=∈[0-9]，其值越接近0，則一致性越差；其值越接近1，則一致性越強(qiáng)。

此時(shí)，各個(gè)子系統(tǒng)作為二級(jí)底層指標(biāo)，對(duì)其上級(jí)子系統(tǒng)的有序度的功效計(jì)算公式為：

式中：Ui為二級(jí)底層指標(biāo)對(duì)其上級(jí)子系統(tǒng)的有序度的功效；ωij為一級(jí)底層指標(biāo)的權(quán)重。由此可推出各個(gè)指標(biāo)相對(duì)上級(jí)子系統(tǒng)的功效。

3.2 構(gòu)造耦合度模型

依據(jù)物理學(xué)范疇的容量耦合概念，可以建立石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合度模型[11]。

兩個(gè)系統(tǒng)間的耦合度模型：

由于根據(jù)耦合度模型所得耦合度是與參與耦合子系統(tǒng)的數(shù)量相關(guān)的，不同數(shù)量的子系統(tǒng)間耦合度無(wú)法直接進(jìn)行對(duì)比。故為方便對(duì)比，本文均采用兩個(gè)系統(tǒng)間的耦合度模型，對(duì)于多個(gè)系統(tǒng)間的耦合，可將其拆分為兩個(gè)子系統(tǒng)，由公式(9)求兩子系統(tǒng)的組合功效值，再由兩個(gè)系統(tǒng)間的耦合度模型(10)得出耦合度，最后根據(jù)拆分方式的不同得出的耦合度求均值即為所求。

4 評(píng)價(jià)模型的應(yīng)用

文章以大連某大型石油罐區(qū)為研究對(duì)象進(jìn)行火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合評(píng)價(jià)。

4.1 構(gòu)造耦合度模型

為簡(jiǎn)化計(jì)算，運(yùn)用matlab8.3編程計(jì)算層次分析結(jié)果。求得各指標(biāo)權(quán)重見表3。

表3 風(fēng)險(xiǎn)指標(biāo)權(quán)重

4.2 專家打分及期望求值

在罐區(qū)HSE部門，技術(shù)部門，現(xiàn)場(chǎng)操作工人中共選取10人對(duì)30個(gè)底層指標(biāo)分別進(jìn)行打分，利用matlab8.3對(duì)逆向云發(fā)生器編碼并運(yùn)算，若數(shù)字特征為復(fù)數(shù)，則需要剔除離散程度大的打分并補(bǔ)入新的專家打分，打分為百分制，0～20分為輕度風(fēng)險(xiǎn)；20～40分為低度風(fēng)險(xiǎn)；40～60分為中度風(fēng)險(xiǎn)；60～80分為重度風(fēng)險(xiǎn)；80～100分為極度風(fēng)險(xiǎn)。最后所得期望結(jié)果表4。

表4 底層指標(biāo)期望

4.3 功效值計(jì)算

針對(duì)4.2節(jié)所得期望結(jié)合(7)～(9)式計(jì)算各項(xiàng)指標(biāo)的功效值。計(jì)算結(jié)果見表5。

表5 風(fēng)險(xiǎn)指標(biāo)功效

4.4 一級(jí)指標(biāo)雙因素耦合度計(jì)算

一級(jí)指標(biāo)雙因素耦合共有10種耦合方式，包括人-機(jī)耦合、人-環(huán)耦合、人-管耦合、人-消耦合、機(jī)-環(huán)耦合、機(jī)-管耦合、機(jī)-消耦合、環(huán)-管耦合、環(huán)-消耦合、管-消耦合，分別記為C21、C22、C23、C24、C25、C26、C27、C28、C29、C2X。依照公式(10)計(jì)算耦合度結(jié)果見表5。

4.5 標(biāo)多因素一級(jí)指耦合度計(jì)算

一級(jí)指標(biāo)多因素耦合包括包括人-機(jī)-環(huán)耦合、人-機(jī)-管耦合、人-機(jī)-消耦合、人-環(huán)-管耦合、人-環(huán)-消耦合、人-管-消耦合、機(jī)-環(huán)-管耦合、機(jī)-環(huán)-消耦合、機(jī)-管-消耦合、環(huán)-管-消耦合、人-機(jī)-環(huán)-管耦合、人-機(jī)-環(huán)-消耦合、人-機(jī)-管-消耦合、人-環(huán)-管-消耦合、機(jī)-環(huán)-管-消耦合、人-機(jī)-環(huán)-管-消耦合，分別記為C31、C32、C33、C34、C35、C36、C37、C38、C39、C3X、C41、C42、C43、C44、C45、C5。依 照3.2節(jié)的拆分原則及式(10)計(jì)算耦合度結(jié)果見表6。

4.6 耦合度計(jì)算結(jié)果分析

將所求得耦合度按照降序排列可得：

(1)雙因素耦合中，人-環(huán)耦合、人-管耦合的耦合度相對(duì)最高；其中有消防風(fēng)險(xiǎn)因素參與的耦合方式耦合度最小 ；不考慮消防風(fēng)險(xiǎn)因素參與的雙因素耦合方式，機(jī)、管風(fēng)險(xiǎn)因素的非全面耦合的耦合度大于其全面耦合的耦合度；人、環(huán)風(fēng)險(xiǎn)因素的非全面耦合的耦合度小于全面耦合的耦合度。

(2)三因素耦合中人-機(jī)-環(huán)耦合、人-環(huán)-管耦合、人-機(jī)-管耦合的耦合度相對(duì)最高，其中有消防風(fēng)險(xiǎn)因素參與的耦合方式耦合度最小；不考慮消防風(fēng)險(xiǎn)因素，所有耦合方式中有人的風(fēng)險(xiǎn)因素參與的三風(fēng)險(xiǎn)因素耦合方式耦合度最高。

(3)四風(fēng)險(xiǎn)因素耦合中，人-機(jī)-環(huán)-管耦合、人-機(jī)-環(huán)-消耦合的耦合度相對(duì)較高；有人的風(fēng)險(xiǎn)因素參與的耦合形式耦合度較高，有消防風(fēng)險(xiǎn)因素參與的耦合形式耦合度較低；人、機(jī)風(fēng)險(xiǎn)因素全面耦合大于其非全面耦合。

表6 耦合度計(jì)算結(jié)果

5 結(jié)語(yǔ)

(1)針對(duì)石油罐區(qū)火災(zāi)爆炸事故，在傳統(tǒng)的人-機(jī)-環(huán)-管四類風(fēng)險(xiǎn)因素基礎(chǔ)上引入消防風(fēng)險(xiǎn)因素，建立以人-機(jī)-環(huán)-管-消為一級(jí)風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)的石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)評(píng)價(jià)體系，并結(jié)合層次分析法確定各項(xiàng)指標(biāo)的權(quán)重。

(2)逆向云模型的應(yīng)用將專家對(duì)于各風(fēng)險(xiǎn)指標(biāo)的評(píng)價(jià)分?jǐn)?shù)轉(zhuǎn)化為數(shù)字特征，并對(duì)其合理性進(jìn)行檢驗(yàn)，從而降低了專家打分的離散程度，使得評(píng)價(jià)方法更科學(xué)。

(3)運(yùn)用耦合度模型對(duì)石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)因素耦合度進(jìn)行評(píng)價(jià)，定量評(píng)價(jià)了各種耦合形式的耦合度強(qiáng)弱關(guān)系，對(duì)安全管理工作具有重要意義。

(4)從耦合度模型的應(yīng)用可知，有人的風(fēng)險(xiǎn)因素參與的耦合方式耦合度較高，有消防風(fēng)險(xiǎn)因素參與的耦合方式耦合度較低。因此在該罐區(qū)安全管理工作中應(yīng)充分發(fā)揮人的主觀能動(dòng)性，提升人員的安全意識(shí)及專業(yè)素質(zhì)，避免其與其他風(fēng)險(xiǎn)因素耦合。