(東北石油大學(xué) 機(jī)械科學(xué)與工程學(xué)院，黑龍江大慶 163318)

0 引言

由于石油儲(chǔ)罐中儲(chǔ)存的石油以及其他物質(zhì)具有易燃、易爆炸、易揮發(fā)、易產(chǎn)生靜電、易擴(kuò)散等特性，因此石油罐區(qū)易發(fā)生火災(zāi)爆炸事故，造成巨大的人員傷亡及財(cái)產(chǎn)損失。

國(guó)內(nèi)外對(duì)于石油罐區(qū)安全領(lǐng)域已有一定的研究。Moshashaei等[1]通過研究浮頂油罐火災(zāi)爆炸事故，總結(jié)火災(zāi)爆炸事故原因，并運(yùn)用分析網(wǎng)絡(luò)及Super Decisions軟件進(jìn)行加權(quán)和優(yōu)先級(jí)排序，確定了浮頂罐火災(zāi)爆炸事故的11個(gè)主要標(biāo)準(zhǔn)及71個(gè)子標(biāo)準(zhǔn)。Halloul等[2]采用模糊故障樹分析與專家訪談相結(jié)合的方法對(duì)原油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)進(jìn)行識(shí)別與評(píng)估。Argyropoulos等[3]通過與希臘的工業(yè)安全專家進(jìn)行討論，提出在塞韋索的指令的框架下、采用檢查表法對(duì)儲(chǔ)罐進(jìn)行風(fēng)險(xiǎn)識(shí)別，并根據(jù)導(dǎo)致罐區(qū)發(fā)生事故后果的嚴(yán)重性進(jìn)行排序。Guo等[4]則針對(duì)現(xiàn)有的儲(chǔ)油罐區(qū)環(huán)境風(fēng)險(xiǎn)，預(yù)測(cè)和評(píng)價(jià)了油罐爆炸引起的熱輻射和沖擊波的影響。包其富等[5]采用層次分析-模糊綜合評(píng)價(jià)方法實(shí)現(xiàn)單個(gè)儲(chǔ)罐的安全狀況分級(jí)，并結(jié)合區(qū)域定量風(fēng)險(xiǎn)評(píng)價(jià)方法實(shí)現(xiàn)石油罐區(qū)整體風(fēng)險(xiǎn)評(píng)價(jià)。任常興等[6]引入火災(zāi)風(fēng)險(xiǎn)嚴(yán)重度指數(shù)評(píng)價(jià)石油罐區(qū)火災(zāi)風(fēng)險(xiǎn)相對(duì)影響水平及等級(jí)。牛淑婷[7]提出主要的引發(fā)石油罐區(qū)火災(zāi)爆炸事故的風(fēng)險(xiǎn)因素,并制定相應(yīng)的針對(duì)性措施,加強(qiáng)安全管理。

分析國(guó)內(nèi)外學(xué)者的研究可以看出，當(dāng)前對(duì)于石油罐區(qū)的風(fēng)險(xiǎn)評(píng)估以及風(fēng)險(xiǎn)管理的研究大多從單因素角度出發(fā)，而沒有從風(fēng)險(xiǎn)因素耦合的角度進(jìn)行研究，且石油罐區(qū)火災(zāi)爆炸事故多為不同風(fēng)險(xiǎn)因素相互作用、相互耦合而發(fā)生的。因此，為有效避免石油罐區(qū)發(fā)生火災(zāi)爆炸事故，應(yīng)從風(fēng)險(xiǎn)耦合的角度出發(fā)，研究致使罐區(qū)發(fā)生火災(zāi)爆炸事故的多種風(fēng)險(xiǎn)因素的耦合特性，建立罐區(qū)風(fēng)險(xiǎn)因素耦合規(guī)律，為石油罐區(qū)實(shí)現(xiàn)在風(fēng)險(xiǎn)閾值內(nèi)的動(dòng)態(tài)平衡提供理論依據(jù)。

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

通過查閱相關(guān)文獻(xiàn)[1-10]，將石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)分為人員風(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)因素，具體風(fēng)險(xiǎn)因素辨識(shí)見表1。

表1 風(fēng)險(xiǎn)因素辨識(shí)

2 石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合類型及機(jī)理

風(fēng)險(xiǎn)耦合是指系統(tǒng)內(nèi)部可導(dǎo)致系統(tǒng)穩(wěn)定性失衡的因素之間相互影響、相互促進(jìn)，致使系統(tǒng)風(fēng)險(xiǎn)值上升或下降的現(xiàn)象[11]。如果數(shù)種風(fēng)險(xiǎn)因素發(fā)生耦合后導(dǎo)致系統(tǒng)風(fēng)險(xiǎn)值上升，此現(xiàn)象屬于風(fēng)險(xiǎn)正耦合；如果數(shù)種風(fēng)險(xiǎn)因素發(fā)生耦合后導(dǎo)致系統(tǒng)風(fēng)險(xiǎn)值降低，則此現(xiàn)象稱之為風(fēng)險(xiǎn)負(fù)耦合。本文主要研究風(fēng)險(xiǎn)因素進(jìn)行耦合過程后致使系統(tǒng)風(fēng)險(xiǎn)值上升的現(xiàn)象，即風(fēng)險(xiǎn)正耦合現(xiàn)象。

2.1 石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合類型

將上文提到的石油罐區(qū)系統(tǒng)火災(zāi)爆炸事故風(fēng)險(xiǎn)(分為人、機(jī)、環(huán)、管、消5個(gè)一級(jí)風(fēng)險(xiǎn)因素)分別記為a,b,c,d,e。依照此分類，石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合可分為單因素耦合、雙因素耦合及多因素耦合。

單因素耦合是指同一個(gè)一級(jí)風(fēng)險(xiǎn)因素分類下的二級(jí)風(fēng)險(xiǎn)因子之間發(fā)生耦合。這種耦合形式通常屬于最常見的耦合形式，但耦合后影響程度較小，一般不會(huì)超出系統(tǒng)整體風(fēng)險(xiǎn)閾值。

雙因素耦合是指2個(gè)不同一級(jí)風(fēng)險(xiǎn)因素分類下的二級(jí)風(fēng)險(xiǎn)因子之間發(fā)生耦合。包括人-機(jī)耦合風(fēng)險(xiǎn)、人-環(huán)耦合風(fēng)險(xiǎn)、人-管耦合風(fēng)險(xiǎn)、人-消耦合風(fēng)險(xiǎn)、機(jī)-環(huán)耦合風(fēng)險(xiǎn)、機(jī)-管耦合風(fēng)險(xiǎn)、機(jī)-消耦合風(fēng)險(xiǎn)、環(huán)-管耦合風(fēng)險(xiǎn)、環(huán)-消耦合風(fēng)險(xiǎn)、管-消耦合風(fēng)險(xiǎn)，分別記為T21(a,b)，T22(a,c)，T23(a,d)，T24(a,e)，T25(b,c)，T26(b,d)，T27(b,e)，T28(c,d)，T29(c,e)，T2Ⅹ(d,e)。

多因素耦合是指3個(gè)及3個(gè)以上不同一級(jí)風(fēng)險(xiǎn)因素分類下的二級(jí)風(fēng)險(xiǎn)因子之間發(fā)生耦合。包括人-機(jī)-環(huán)耦合風(fēng)險(xiǎn)、人-機(jī)-管耦合風(fēng)險(xiǎn)、人-機(jī)-消耦合風(fēng)險(xiǎn)、人-環(huán)-管耦合風(fēng)險(xiǎn)、人-環(huán)-消耦合風(fēng)險(xiǎn)、人-管-消耦合風(fēng)險(xiǎn)、機(jī)-環(huán)-管耦合風(fēng)險(xiǎn)、機(jī)-環(huán)-消耦合風(fēng)險(xiǎn)、機(jī)-管-消耦合風(fēng)險(xiǎn)、環(huán)-管-消耦合風(fēng)險(xiǎn)、人-機(jī)-環(huán)-管耦合風(fēng)險(xiǎn)、人-機(jī)-環(huán)-消耦合風(fēng)險(xiǎn)、人-機(jī)-管-消耦合風(fēng)險(xiǎn)、人-環(huán)-管-消耦合風(fēng)險(xiǎn)、機(jī)-環(huán)-管-消耦合風(fēng)險(xiǎn)、人-機(jī)-環(huán)-管-消耦合風(fēng)險(xiǎn),分別記為T31(a,b,c),T32(a,b,d),T33(a,b,e),T34(a,c,d),T35(a，c，e),T36(a,d,e),T37(b,c,d),T38(b,c,e),T39(b,d,e),T3Ⅹ(c,d,e),T41(a,b,c,d),T42(a,b,c,e),T43(a,b,d,e),T44(a,c,d,e),T45(b,c,d,e),T51(a,b,c,d,e)。

2.2 石油罐區(qū)風(fēng)險(xiǎn)耦合機(jī)理

石油罐區(qū)系統(tǒng)中存在諸多風(fēng)險(xiǎn)因素，其中單個(gè)風(fēng)險(xiǎn)因素往往并不會(huì)直接致使石油罐區(qū)發(fā)生火災(zāi)爆炸事故，而是各種風(fēng)險(xiǎn)因素之間發(fā)生耦合導(dǎo)致原油罐區(qū)整體風(fēng)險(xiǎn)值上升至超出系統(tǒng)自身安全水平閾值，從而導(dǎo)致石油罐區(qū)發(fā)生火災(zāi)爆炸事故。國(guó)內(nèi)學(xué)者大部分只對(duì)原油罐區(qū)火災(zāi)爆炸事故基于多方面風(fēng)險(xiǎn)因素進(jìn)行致因分析，研究單個(gè)風(fēng)險(xiǎn)因素與事故之間的聯(lián)系及影響，卻忽略了風(fēng)險(xiǎn)因素之間的耦合作用。

圖1 石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合機(jī)理示意

如圖1所示，在石油罐區(qū)系統(tǒng)中存在人、機(jī)、環(huán)、管、消5個(gè)子系統(tǒng)，各子系統(tǒng)分類下的二級(jí)風(fēng)險(xiǎn)因子之間相互促進(jìn)、疊加，發(fā)生耦合，擴(kuò)大原有風(fēng)險(xiǎn)值，使其突破子系統(tǒng)安全防御閾值并與其他子系統(tǒng)的風(fēng)險(xiǎn)因子再次進(jìn)行耦合的過程中產(chǎn)生聚類突變，經(jīng)過耦合過程，打破原系統(tǒng)的風(fēng)險(xiǎn)平衡狀態(tài)，或超過石油罐區(qū)系統(tǒng)的風(fēng)險(xiǎn)閥值整體水平，最終使風(fēng)險(xiǎn)因素耦合后的整體風(fēng)險(xiǎn)水平高于風(fēng)險(xiǎn)耦合前的風(fēng)險(xiǎn)水平，從而發(fā)生火災(zāi)爆炸事故。

3 風(fēng)險(xiǎn)耦合模型的構(gòu)建

N-K模型起源于信息理論，是一種解決復(fù)雜系統(tǒng)計(jì)算的通用模型。其中包含參數(shù)N和K，N是指構(gòu)成系統(tǒng)要素的個(gè)數(shù)，K是指網(wǎng)絡(luò)中相互依賴關(guān)系的個(gè)數(shù)，其中0≤K≤N-1[12]。

3.1 確定石油罐區(qū)風(fēng)險(xiǎn)耦合信息交互公式

通過對(duì)石油罐區(qū)系統(tǒng)運(yùn)行過程中的5類風(fēng)險(xiǎn)因素交互信息進(jìn)行計(jì)算，對(duì)由于風(fēng)險(xiǎn)耦合作用所形成的新的風(fēng)險(xiǎn)狀態(tài)進(jìn)行評(píng)價(jià)，通過各類風(fēng)險(xiǎn)耦合方式發(fā)生次數(shù)判斷各類耦合方式的耦合概率。即若某種耦合形式發(fā)生的次數(shù)越多，則該耦合方式的耦合發(fā)生概率越大；通過耦合值的大小反映耦合風(fēng)險(xiǎn)值的大小，即某種耦合形式所得耦合值越大，表明該耦合方式耦合風(fēng)險(xiǎn)值越高，發(fā)生火災(zāi)爆炸事故的概率越高。交互信息計(jì)算公式如下：

·log2[Phijkl/(Ph....·P.i...

·P..j..·P...k.·P....l)]

(1)

其中,h=1,2,…,H;i=1,2,…,I;j=1,2,…,J;k=1,2,…,K;l=1,2,…,L。

式中Phijkl——石油罐區(qū)系統(tǒng)中人的風(fēng)險(xiǎn)處于h狀態(tài),設(shè)備風(fēng)險(xiǎn)處于i狀態(tài)，環(huán)境風(fēng)險(xiǎn)處于j狀態(tài)，管理風(fēng)險(xiǎn)處于k狀態(tài)，消防風(fēng)險(xiǎn)處于l狀態(tài)時(shí)，5種風(fēng)險(xiǎn)發(fā)生耦合的概率。

T值是對(duì)于石油罐區(qū)系統(tǒng)火災(zāi)爆炸風(fēng)險(xiǎn)耦合的數(shù)量化評(píng)估。T值越大，表明在該狀態(tài)下石油罐區(qū)系統(tǒng)安全水平越低，發(fā)生火災(zāi)爆炸事故的概率就越大。

3.2 構(gòu)建石油罐區(qū)火災(zāi)爆炸事故雙因素風(fēng)險(xiǎn)耦合模型

石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合多以局部耦合形式發(fā)生，局部耦合可以由兩種不同類型的風(fēng)險(xiǎn)因素進(jìn)行耦合，例如人-機(jī)耦合形式，即人的風(fēng)險(xiǎn)因素與設(shè)備風(fēng)險(xiǎn)因素進(jìn)行耦合致使石油罐區(qū)系統(tǒng)安全平衡狀態(tài)失穩(wěn)，從而發(fā)生火災(zāi)爆炸事故，其他類型風(fēng)險(xiǎn)因素在此次耦合過程中不具有決定性影響，這種耦合形式稱為局部耦合。因此，可以通過計(jì)算發(fā)生兩種風(fēng)險(xiǎn)類型之間的信息交互來(lái)研究局部耦合情況。以人-機(jī)耦合為例構(gòu)建風(fēng)險(xiǎn)耦合模型如式(2)，同理可構(gòu)建全部雙因素風(fēng)險(xiǎn)耦合模型。

(2)

3.3 構(gòu)建石油罐區(qū)火災(zāi)爆炸事故多因素風(fēng)險(xiǎn)耦合模型

同理，局部耦合也可以由3種或4種不同類型風(fēng)險(xiǎn)因素進(jìn)行耦合，例如人-機(jī)-環(huán)耦合，此種耦合方式中管理風(fēng)險(xiǎn)與消防風(fēng)險(xiǎn)無(wú)決定性影響；人-機(jī)-環(huán)-管耦合，此種耦合方式中消防風(fēng)險(xiǎn)無(wú)決定性影響。以人-機(jī)-環(huán)耦合為例構(gòu)建三因素風(fēng)險(xiǎn)耦合模型如式(3)，以人-機(jī)-環(huán)-管耦合為例構(gòu)建四因素風(fēng)險(xiǎn)耦合模型如式(4)，同理可構(gòu)建全部多因素風(fēng)險(xiǎn)耦合模型。

·P.i...·P..j..)]

(3)

·P.i...·P..j..·P...k .)]

(4)

4 石油罐區(qū)火災(zāi)爆炸風(fēng)險(xiǎn)耦合實(shí)證分析

4.1 風(fēng)險(xiǎn)耦合計(jì)算

根據(jù)相關(guān)資料[8]及實(shí)地調(diào)研統(tǒng)計(jì)出1972～2018年國(guó)內(nèi)外石油罐區(qū)火災(zāi)爆炸事故案例82起，通過對(duì)事故統(tǒng)計(jì)數(shù)據(jù)分析，得出風(fēng)險(xiǎn)耦合次數(shù)及風(fēng)險(xiǎn)耦合頻率，風(fēng)險(xiǎn)未發(fā)生用“0”表示，風(fēng)險(xiǎn)發(fā)生用“1”表示，則5類風(fēng)險(xiǎn)分類共有32種可能的風(fēng)險(xiǎn)耦合形式結(jié)果見表2。

表2 1972～2018年石油罐區(qū)火災(zāi)爆炸風(fēng)險(xiǎn)耦合次數(shù)及頻率

為計(jì)算T值，首先需要計(jì)算以下不同方式風(fēng)險(xiǎn)耦合的概率。

(1)單因素耦合不同情況下發(fā)生的概率。

人的風(fēng)險(xiǎn)因素不參與風(fēng)險(xiǎn)耦合時(shí)石油罐區(qū)發(fā)生火災(zāi)爆炸事故的概率為P0....=P00000+P01000+P00100+P00010+P00001+P01100+P01010+P01001+P00110+P00101+P00011+P01110+P01101+P01011+P00111+P01111=0.0000+0.0000+0.0976+0.0000+0.0000+0.0488+0.0122+0.0000+0.0244+0.0366+0.0000+0.0122+0.0366+0.0366+0.0122+0.0244=0.3416,同理可求得其余結(jié)果。計(jì)算結(jié)果見表3。

表3 單因素風(fēng)險(xiǎn)變化概率

(2)雙因素耦合不同情況下發(fā)生的概率。

人的風(fēng)險(xiǎn)因素與設(shè)備風(fēng)險(xiǎn)因素均不參與風(fēng)險(xiǎn)耦合時(shí)石油罐區(qū)發(fā)生火災(zāi)爆炸事故的概率為P00...=P00000+P00100+P00010+P00001+P00110+P00101+P00011+P00111=0.0000+0.0976+0.0000+0.0000+0.0244+0.0366+0.0000+0.0122=0.1708,同理可求得其余結(jié)果,計(jì)算結(jié)果如表4所示。

表4 雙因素風(fēng)險(xiǎn)變化概率

(3)多因素耦合不同情況下發(fā)生的概率。

人的風(fēng)險(xiǎn)因素與設(shè)備風(fēng)險(xiǎn)因素及消防風(fēng)險(xiǎn)因素均不參與風(fēng)險(xiǎn)耦合時(shí)石油罐區(qū)發(fā)生火災(zāi)爆炸事故的概率為P00..0=P00000+P00100+P00010+P00110=0.0000+0.0976+0.0000+0.0244=0.1220。同理可以求得其余結(jié)果，計(jì)算結(jié)果見表5。

表5 多因素風(fēng)險(xiǎn)變化概率

(續(xù)表5)

通過已構(gòu)建風(fēng)險(xiǎn)耦合模型,可計(jì)算出石油罐區(qū)各項(xiàng)風(fēng)險(xiǎn)耦合的T值，計(jì)算結(jié)果如下：

T21=0.023467427,T22=0.156375293,

T23=0.100106004,T24=0.039581477,

T25=0.016828638,T26=0.000016532,

T27=0.003747548,T28=0.217325925,

T29=0.014844345,T2Ⅹ=0.002121561,

T31=0.272014739,T32=0.204785229,

T33=0.120081086,T34=0.447377565,

T35=0.252879588,T36=0.220010281,

T37=0.331434599,T38=0.095687131,

T39=0.079904587,T3Ⅹ=0.260460098,

T41=0.621386225,T42=0.319736817,

T43=0.276255606,T44=0.515568031,

T45=0.365189486,T5=0.721525291。

對(duì)所得T值按降序排列可得：

T5>T41>T44>T34>T45>T37>T42>T43>T31>T3Ⅹ>T35>T36>T28>T32>T22>T33>T23>T38>T39>T24>T21>T29>T25>T27>T2Ⅹ>T26

按耦合因素?cái)?shù)量升序作T值分布圖，如圖2所示。

圖2 不同類型風(fēng)險(xiǎn)耦合T值分布圖

求雙因素耦合T值、三因素耦合T值、四因素耦合T值平均值得：

4.2 計(jì)算結(jié)果分析

(1)隨著參與耦合風(fēng)險(xiǎn)數(shù)量的增加，T值總體呈上升趨勢(shì)；且隨著耦合風(fēng)險(xiǎn)種類的增加，T值均值遞增，即參與風(fēng)險(xiǎn)耦合的風(fēng)險(xiǎn)因素越多，石油罐區(qū)發(fā)生火災(zāi)爆炸事故風(fēng)險(xiǎn)越大。管理者應(yīng)找出關(guān)鍵風(fēng)險(xiǎn)因素，消弱風(fēng)險(xiǎn)因素風(fēng)險(xiǎn)水平，減少參與耦合的風(fēng)險(xiǎn)種類。

(2)在雙因素風(fēng)險(xiǎn)耦合中，環(huán)-管耦合、人-環(huán)耦合、人-管耦合的T值相對(duì)較大，具有較強(qiáng)的耦合性；在三因素風(fēng)險(xiǎn)耦合中，人-環(huán)-管耦合、機(jī)-環(huán)-管耦合、人-機(jī)-環(huán)耦合具有較強(qiáng)的耦合性；在四因素風(fēng)險(xiǎn)耦合中，人-機(jī)-環(huán)-管耦合與人-環(huán)-管-消耦合具有較強(qiáng)的耦合性。以上耦合形式應(yīng)作為控制石油罐區(qū)火災(zāi)爆炸風(fēng)險(xiǎn)優(yōu)先考慮的耦合形式。

(3)由單因素變化概率值及風(fēng)險(xiǎn)耦合T值可以得出，在石油罐區(qū)火災(zāi)爆炸5個(gè)風(fēng)險(xiǎn)因素中，有人的風(fēng)險(xiǎn)因素參與的風(fēng)險(xiǎn)耦合次數(shù)最多，且風(fēng)險(xiǎn)耦合的T值較大。

5 結(jié)語(yǔ)

(1)將消防風(fēng)險(xiǎn)因素加入石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)因素分類，并將其引入風(fēng)險(xiǎn)耦合T值計(jì)算，提高了石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合分析的全面性。

(2)對(duì)石油罐區(qū)火災(zāi)爆炸事故風(fēng)險(xiǎn)耦合作用進(jìn)行分析，建立了耦合視角下火災(zāi)爆炸事故致因模型，指出石油罐區(qū)發(fā)生火災(zāi)爆炸事故是由致使風(fēng)險(xiǎn)擴(kuò)大導(dǎo)致的，并利用N-K模型對(duì)各類耦合形式進(jìn)行度量。

(3)在進(jìn)行石油罐區(qū)的風(fēng)險(xiǎn)管理過程中，應(yīng)避免多個(gè)因素發(fā)生耦合作用，當(dāng)某一因素出現(xiàn)危險(xiǎn)源時(shí)，其他因素存在風(fēng)險(xiǎn)的可能也很大，在及時(shí)消除該風(fēng)險(xiǎn)因素的同時(shí)，應(yīng)避免與其他潛在風(fēng)險(xiǎn)因素發(fā)生耦合效應(yīng)，尤其注意避免具有較強(qiáng)耦合性的耦合方式發(fā)生。

(4)在石油罐區(qū)火災(zāi)爆炸風(fēng)險(xiǎn)耦合中，人是關(guān)鍵耦合因素，應(yīng)充分發(fā)揮人的主觀能動(dòng)性，提高管理強(qiáng)度，并加強(qiáng)對(duì)設(shè)備因素，環(huán)境因素以及消防因素的管理，從而降低罐區(qū)發(fā)生火災(zāi)爆炸的頻率。