孫佳佳

（江蘇省安全生產(chǎn)科學(xué)研究院，江蘇 南京 210042）

近幾年，隨著國(guó)家推進(jìn)企業(yè)入園進(jìn)區(qū)，企業(yè)分布呈現(xiàn)區(qū)域集中式，更加規(guī)模化和集中化。但也正因如此，園區(qū)企業(yè)密度增加，企業(yè)的突發(fā)事故更容易對(duì)本企業(yè)和周邊企業(yè)產(chǎn)生聯(lián)鎖事故，發(fā)生特大安全事故，從而引發(fā)特大突發(fā)性環(huán)境污染事件。2013年8月13日天津港爆炸事故正是由于硝化棉燃燒引發(fā)的硝酸銨等129種危險(xiǎn)化學(xué)物質(zhì)的聯(lián)鎖爆炸事故，對(duì)局部區(qū)域的大氣環(huán)境、水環(huán)境和土壤環(huán)境造成了不同程度的污染[1]。一些學(xué)者曾經(jīng)提出突發(fā)性事故的多米諾效應(yīng)，但在目前的環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)中，多米諾效應(yīng)很少被提及，處于被忽視的地位。本文通過(guò)石化倉(cāng)儲(chǔ)企業(yè)的安全及環(huán)境特點(diǎn)，探討環(huán)境風(fēng)險(xiǎn)多米諾效應(yīng)的實(shí)際應(yīng)用，希望對(duì)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)工作提出有益的幫助。

1 環(huán)境多米諾效應(yīng)

1.1 環(huán)境多米諾效應(yīng)概念

Cozzanid等在2005年對(duì)相關(guān)的文獻(xiàn)做過(guò)詳細(xì)的總結(jié)，提出相對(duì)完整的多米諾效應(yīng)的定義，即由初始時(shí)間發(fā)生，引發(fā)鄰近一個(gè)或多個(gè)設(shè)備發(fā)生事故，從而導(dǎo)致總事故后果比初始事故后果更加嚴(yán)重的現(xiàn)象為多米諾效應(yīng)[2]。

在化工企業(yè)或石化倉(cāng)儲(chǔ)企業(yè)中，若一個(gè)工藝或存儲(chǔ)設(shè)備、單元發(fā)生事故時(shí)，可能會(huì)伴隨其它工藝或存儲(chǔ)設(shè)備、單元的破壞，從而引發(fā)二次、三次甚至更為嚴(yán)重的事故。當(dāng)突發(fā)性事故多米諾事故對(duì)環(huán)境產(chǎn)生危害時(shí)，即為環(huán)境風(fēng)險(xiǎn)多米諾事故。

環(huán)境風(fēng)險(xiǎn)多米諾效應(yīng)是指一定區(qū)域內(nèi)存在多個(gè)環(huán)境風(fēng)險(xiǎn)源，由某些風(fēng)險(xiǎn)因素影響引發(fā)初始事故，當(dāng)初始事故釋放足夠能量或物質(zhì)后，與外界環(huán)境因素共同作用，使鄰近的環(huán)境風(fēng)險(xiǎn)源受到影響而發(fā)生事故，兩個(gè)或多個(gè)環(huán)境風(fēng)險(xiǎn)事故相互作用，連鎖發(fā)生，造成有毒有害物質(zhì)大量進(jìn)入周?chē)h(huán)境。

1.2 多米諾效應(yīng)發(fā)生條件

一般來(lái)說(shuō)，多米諾事故的發(fā)生必須具備以下條件。

（1）距離條件：初始事故的影響范圍內(nèi)需存在敏感單元；

（2）能量條件：初始事故的破壞效應(yīng)，即能量足夠大能夠?qū)е掠绊懛秶鷥?nèi)的敏感單元被破壞；

（3）概率條件：初始事故以及所引發(fā)的二次事故的發(fā)生破壞概率應(yīng)當(dāng)足夠大，即為可能發(fā)生事故。

當(dāng)上述三個(gè)條件同時(shí)滿足時(shí)，才會(huì)觸發(fā)多米諾效應(yīng)。

1.3 多米諾效應(yīng)作用機(jī)理及閾值研究

觸發(fā)多米諾效應(yīng)的事故主要有三種形式，分別為泄露、火災(zāi)及爆炸。其中泄露引發(fā)的機(jī)理主要為有毒有害物質(zhì)泄漏導(dǎo)致操作人員中毒，使操作人員神志不清，現(xiàn)場(chǎng)設(shè)備無(wú)人看管等情況，從而無(wú)法進(jìn)行正常操作，引發(fā)多米諾效應(yīng)。由于該情況主觀因素較大，無(wú)法進(jìn)行定量判定，故本文重點(diǎn)分析火災(zāi)及爆炸這兩種形式。

1.3.1 火災(zāi)機(jī)理及閾值研究

火災(zāi)一般分為池火和火球兩種事故類型。造成多米諾事故一般有兩種方式：一是火焰直接接觸，二是熱輻射造成設(shè)備失效。

池火是易燃液體泄露，在地面形成液池，液池被點(diǎn)燃形成火災(zāi)。當(dāng)形成的池火完全包圍目標(biāo)設(shè)備，與設(shè)備直接接觸時(shí)，較易引發(fā)多米諾事故。當(dāng)池火距離目標(biāo)設(shè)備存在一定距離，無(wú)法直接接觸目標(biāo)設(shè)備時(shí)，可能會(huì)通過(guò)熱輻射破壞目標(biāo)設(shè)備。

熱輻射造成設(shè)備失效需要兩個(gè)條件：熱輻射通量和時(shí)間。Valerio Cozzani等研究人員考慮目標(biāo)設(shè)備類型和火災(zāi)持續(xù)時(shí)間，對(duì)熱輻射的多米諾效應(yīng)進(jìn)行數(shù)值模擬，并與相關(guān)文獻(xiàn)進(jìn)行對(duì)比分析，總結(jié)出若火災(zāi)持續(xù)時(shí)間在10 min以上，對(duì)于常壓容器的破壞閾值為15 kW/m2，高壓容器破壞閾值為50 kW/m2[3]。

1.3.2 爆炸機(jī)理及閾值研究

當(dāng)企業(yè)發(fā)生蒸汽云爆炸或沸騰液體擴(kuò)散蒸汽爆炸（BLEVE）事故時(shí)，會(huì)產(chǎn)生沖擊波及拋射碎片等能量載體。這些能量載體與目標(biāo)設(shè)備相互作用，可能造成目標(biāo)設(shè)備的破壞。

沖擊波破壞效應(yīng)可根據(jù)沖擊波超壓破壞閾值進(jìn)行定量評(píng)估。Valerio Cozzani和Ernesto Salzano[3]兩位學(xué)者在經(jīng)過(guò)文獻(xiàn)資料分析比較的基礎(chǔ)上，對(duì)超壓引發(fā)多米諾效應(yīng)進(jìn)行研究。他們將設(shè)備分為常壓容器、壓力容器、加長(zhǎng)容器、小容器四類，并針對(duì)不同容器確定破壞閾值。具體見(jiàn)表1。

表1 超壓破壞閾值

并且，爆炸產(chǎn)生的拋射物也同樣能引發(fā)多米諾效應(yīng)，李驍驊[4]對(duì)液化石油氣儲(chǔ)罐爆炸產(chǎn)生的碎片的破壞作用進(jìn)行了研究，拋射物引發(fā)多米諾效應(yīng)的參數(shù)包括穿透厚度、沖擊能、碰撞概率，研究中表明常見(jiàn)的拋射物拋射距離閾值為300 m。

1.4 多米諾環(huán)境風(fēng)險(xiǎn)計(jì)算

1.4.1 泄露污染源計(jì)算

當(dāng)物料以液體形式泄漏到地面形成液池，在液池表面氣流運(yùn)動(dòng)作用下發(fā)生質(zhì)量蒸發(fā)現(xiàn)象，從而擴(kuò)散進(jìn)入大氣。

根據(jù)《環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)導(dǎo)則》1附錄A.2[5]，液體質(zhì)量蒸發(fā)速率可以由下式計(jì)算得出：

式中：MW—質(zhì)量蒸發(fā)速率，kg/s；

a，n—大氣穩(wěn)定度系數(shù)；

Ps—液體表面蒸汽壓，Pa；

M—物質(zhì)摩爾質(zhì)量，g/mol；

R—通用氣體系數(shù)，J/（mol·k）；

Ta—周?chē)h(huán)境溫度，K；

u—風(fēng)速，m/s；

r—液池半徑，m。

表2 大氣穩(wěn)定度系數(shù)取值

1.4.2 影響預(yù)測(cè)模型

采用多煙團(tuán)模式，計(jì)算公式：

式中：C（x，y，o）：下風(fēng)向地面（x，y）坐標(biāo)處的空氣中污染物濃度，mg/m3；

x0，y0，z0：煙團(tuán)中心坐標(biāo)；

Q：事故期間煙團(tuán)的排放量；

σx，σy，σz：為 x、y、z方向的擴(kuò)散參數(shù)，m。 常取σx=σy；

對(duì)于瞬時(shí)或短時(shí)間事故，可采用下述變天條件下多煙團(tuán)模式：

Q'： 煙團(tuán)排放量 （mg），Q'=Q·Δt為釋放率，mg/s； Δt為時(shí)段長(zhǎng)度，s；

σx,eff，σy,eff，σz,eff：煙團(tuán)在 w 時(shí)段沿 x、y 和 z方向的等效擴(kuò)散參數(shù)（m），可由下式估算：

各個(gè)煙團(tuán)對(duì)某個(gè)關(guān)心點(diǎn)t小時(shí)的濃度貢獻(xiàn)，按下式計(jì)算：

式中：n為需要跟蹤的煙團(tuán)數(shù)，可由下式確定：

式中，f為小于1的系數(shù)，可根據(jù)計(jì)算要求確定。

2 案例研究

2.1 某公司基本情況

某石化倉(cāng)儲(chǔ)企業(yè)位于東部沿海，企業(yè)建有38座地上立式鋼制儲(chǔ)罐，總罐容19.15×104m3，分為四個(gè)罐組，1號(hào)（9座儲(chǔ)罐）及2號(hào)罐組（12座儲(chǔ)罐）存儲(chǔ)易燃易爆液體，3號(hào)罐組（5座儲(chǔ)罐）存儲(chǔ)液化烴，4號(hào)罐組（12座儲(chǔ)罐）存儲(chǔ)有毒且易燃液體。

企業(yè)位于園區(qū)，周邊500 m范圍內(nèi)無(wú)居民；東側(cè)為倉(cāng)儲(chǔ)企業(yè)B，該企業(yè)儲(chǔ)罐中距評(píng)價(jià)企業(yè)儲(chǔ)罐最近為147 m，南側(cè)70 m為海域。

2.2 事故的風(fēng)險(xiǎn)源強(qiáng)

企業(yè)4個(gè)罐組的38座儲(chǔ)罐中，3號(hào)罐組主要存儲(chǔ)丁二烯、丙烯、環(huán)氧丙烷，為甲A類液化烴及閃點(diǎn)極低的甲B類易燃液體；4號(hào)罐組存儲(chǔ)苯及丙烯腈，屬于高毒化學(xué)品。相對(duì)1號(hào)及2號(hào)存儲(chǔ)甲B類易燃液體罐組，3號(hào)及4號(hào)罐組的危險(xiǎn)性較大，且3號(hào)及4號(hào)罐組相鄰，較易發(fā)生多米諾事故，因此，本文選擇3號(hào)及4號(hào)罐組的儲(chǔ)罐進(jìn)行風(fēng)險(xiǎn)源強(qiáng)分析。企業(yè)源強(qiáng)見(jiàn)表3，布置圖見(jiàn)圖1。

表3 企業(yè)源強(qiáng)表

圖1 庫(kù)區(qū)3號(hào)及4號(hào)罐組平面布置圖

2.3 風(fēng)險(xiǎn)計(jì)算

2.3.1 火災(zāi)爆炸風(fēng)險(xiǎn)多米諾影響分析

根據(jù)企業(yè)源強(qiáng)及位置進(jìn)行火災(zāi)爆炸風(fēng)險(xiǎn)計(jì)算，計(jì)算模型如下。

2.3.1.1 蒸汽云爆炸計(jì)算模型

根據(jù) 《安全評(píng)價(jià)員實(shí)用手冊(cè)》第235頁(yè)8.2.2.3易燃?xì)怏w、液體泄漏蒸氣云爆炸事故模擬分析與評(píng)價(jià)（1）①中的計(jì)算公式。

①WTNT=αWfQf/QTNT

式中：WTNT——蒸氣云的TNT當(dāng)量，kg

a——蒸氣云的TNT當(dāng)量系數(shù)

Wf——蒸氣云爆炸中燃燒的總質(zhì)量，kg

Qf——物質(zhì)的燃燒熱，kJ/kg

QTNT——TNT的爆熱，kJ/kg

②N=WTNT/M

式中：N——蒸氣云的TNT摩爾量，mol

M——TNT的分子量，kg

（2）具有可燃性的化學(xué)品燃燒后放出的熱量計(jì)算公式

Q=QfW

式中：Q—物質(zhì)燃燒后放出的熱量，kJ

Qf—物質(zhì)的燃燒熱，kJ/kg

W—物質(zhì)的質(zhì)量，kg

2.3.1.2 池火災(zāi)模擬計(jì)算

（1）確定火焰高度

計(jì)算池火焰高度的經(jīng)驗(yàn)公式如下：

式中：L 為火焰高度（m），D 為池直徑（m），為燃燒速率（kg/m2·s），ρ0為空氣密度（kg/m3），g 為引力常數(shù)。

（2）計(jì)算火焰表面熱通量

假定能量由圓柱形火焰?zhèn)让婧晚?xiàng)部向周?chē)鶆蜉椛?，用下式?jì)算火焰表面的熱通量。

式中：q0——為火焰表面的熱通量，kw/m2；

△Hc——為燃燒熱，kJ/kg；

π——為圓周率；

f——為熱輻射系數(shù)（可取為0.15）；

mf——為燃燒速率，kg/（m2·s），其它符號(hào)同前。

（3）目標(biāo)接收到的熱通量的計(jì)算

目標(biāo)接收到的熱通量q（r）的計(jì)算公式為：

式中：q （r）為目標(biāo)接收到的熱通量（kw/m2），q0為由式（3）計(jì)算的火焰表面的熱通量（kw/m2），r為目標(biāo)到油區(qū)中心的水平距離 （m），V為視角泵數(shù)。

根據(jù)企業(yè)庫(kù)區(qū)存儲(chǔ)物質(zhì)的特點(diǎn)，應(yīng)用中國(guó)安全生產(chǎn)科學(xué)研究院定量風(fēng)險(xiǎn)評(píng)估軟件對(duì)庫(kù)區(qū)的火災(zāi)爆炸進(jìn)行預(yù)測(cè)計(jì)算，給出風(fēng)險(xiǎn)源多米諾事故影響范圍，具體見(jiàn)表4。

表4 多米諾事故影響半徑

通過(guò)表4可以看出，存儲(chǔ)丁二烯、丙烯的301、302、303、304儲(chǔ)罐管道破裂或容器破裂時(shí)會(huì)引發(fā)云爆，多米諾影響半徑在55 m以上；存儲(chǔ)苯的 401、402、406、407、408、409 儲(chǔ)罐管道破裂或容器破裂時(shí)會(huì)引發(fā)池火事故，多米諾影響半徑在20 m以上。3號(hào)罐組、4號(hào)罐組內(nèi)相鄰儲(chǔ)罐均在多米諾半徑范圍內(nèi)；3號(hào)罐組發(fā)生爆炸事故時(shí)也會(huì)影響到4號(hào)罐組。

從物質(zhì)毒性分析，3號(hào)罐組內(nèi)存放丙烯、丁二烯等液化烴類物質(zhì)，泄露爆炸后產(chǎn)生的氣體為CO、CO2、非甲烷甲烴等，物質(zhì)的毒性較??；4號(hào)罐組中存放的均為高毒物質(zhì)，泄露后對(duì)環(huán)境有較大危害。企業(yè)在生產(chǎn)過(guò)程中采取火災(zāi)報(bào)警、緊急切斷、降溫噴淋等一系列的風(fēng)險(xiǎn)控制及應(yīng)急措施，儲(chǔ)罐發(fā)生二次、三次多米諾事故的概率極低，因此，本文主要分析對(duì)環(huán)境的危害性，選擇毒性較高、容量較大的402號(hào)苯儲(chǔ)罐發(fā)生池火事故產(chǎn)生的多米諾效應(yīng)對(duì)環(huán)境產(chǎn)生的影響。

2.3.2 多米諾事故環(huán)境風(fēng)險(xiǎn)影響分析

402號(hào)苯儲(chǔ)罐發(fā)生容器整體破裂或管道完全破裂事故時(shí)，會(huì)發(fā)生池火事故，多米諾影響半徑在32 m，東側(cè)的401號(hào)苯儲(chǔ)罐、西側(cè)的403號(hào)丙烯腈儲(chǔ)罐、南側(cè)的407、408、409號(hào)苯儲(chǔ)罐均在其影響范圍內(nèi)。因多米諾影響半徑內(nèi)儲(chǔ)罐是否發(fā)生事故受風(fēng)向、儲(chǔ)罐類型、儲(chǔ)罐安全裝置穩(wěn)定性等多方面因素影響，影響范圍內(nèi)的儲(chǔ)罐均發(fā)生事故的可能性極低，因此選擇影響范圍內(nèi)毒性最大的403號(hào)丙烯腈儲(chǔ)罐發(fā)生事故進(jìn)行預(yù)測(cè)。

根據(jù)泄露液體形成的液池面積、表面風(fēng)速、環(huán)境溫度、大氣穩(wěn)定度，物質(zhì)飽和蒸汽壓、摩爾質(zhì)量等參數(shù)，計(jì)算泄露的苯及丙烯腈的質(zhì)量蒸發(fā)速率，根據(jù)質(zhì)量蒸發(fā)速率計(jì)算污染物在大氣環(huán)境中的擴(kuò)散情況。

液池面積選取儲(chǔ)罐質(zhì)量蒸發(fā)速率見(jiàn)表5，事故發(fā)生后30 min時(shí)污染物下風(fēng)向濃度見(jiàn)圖2及3，超標(biāo)情況見(jiàn)表6。

表5 液體質(zhì)量蒸發(fā)計(jì)算結(jié)果

圖2 苯（30 min）下風(fēng)向濃度分布

圖3 丙烯腈（30 min）下風(fēng)向濃度分布

表6 污染物下風(fēng)向超標(biāo)情況

從圖2、圖3及表6的可以看出，污染物對(duì)下風(fēng)向的影響較大，最大落地濃度均超過(guò)半致死濃度數(shù)倍至數(shù)十倍，主要原因是本次預(yù)測(cè)選取了較為極端的最不利情況，即儲(chǔ)罐完全爆炸，泄露物料均在防火堤內(nèi)，面源高度較低。從預(yù)測(cè)結(jié)果可以看出，在E-F穩(wěn)定度條件下對(duì)環(huán)境的影響較大，苯的半致死濃度范圍為79.6 m，丙烯腈的半致死濃度范圍為983.7 m。

3 結(jié)論

目前，環(huán)境多米諾效應(yīng)在環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)中處于被忽視的地位，本文提供一個(gè)評(píng)價(jià)思路以供參考。從環(huán)境多米諾的研究現(xiàn)狀來(lái)看，缺乏一套系統(tǒng)的事故預(yù)測(cè)篩選機(jī)制，如何在可能發(fā)生的多米諾事故中選擇有代表性、實(shí)際預(yù)測(cè)意義的預(yù)測(cè)工況和預(yù)測(cè)因子，是我們目前可以研究的方向。建議在今后的研究中，制定一系列的標(biāo)準(zhǔn)規(guī)范，以篩選具有預(yù)測(cè)價(jià)值的多米諾事故。