基于MapInfo多米诺事故热辐射叠加分析与设计

常州大学信息科学与工程学院 刘 俊

常州大学石油工程学院 李俊华 嵇 炜

基于MapInfo多米诺事故热辐射叠加分析与设计

常州大学信息科学与工程学院 刘 俊

常州大学石油工程学院 李俊华 嵇 炜

引入了多米诺的定义，给出事故多米诺效应定量评价的具体程序，在触发多米诺效应发生条件和发生模式的基础上，重点研究了多米诺事故热辐射叠加算法分析，利用Mapinfo设计了多米诺事故热辐射叠加分析软件以及相应的设计步骤，解决较为复杂的重大事故多米诺效应概率叠加计算问题。研究结果可以根据发生概率更加形象地动态地预测事故进一步发展的趋势，为救援指挥提供决策，对相关园区的安全管理和救援决策具有重要指导意义。

多米诺效应；热辐射叠加；MapInfo；概率

1.引言

在我国经济的迅猛发展过程中，各地化工企业也不断增多，各地方政府为了对化工企业更好地进行管理，通常都是将大量的化工企业集中到化工园区。化工园区的建设虽然促进当地经济和化工产业发展，但是同时也带来了严峻的安全问题。由于园区中多数企业生产、储存、使用、运输着大量高温高压、易燃易爆、有毒有害的危险化学品，稍有不慎就有可能发生事故，破坏性很大。园区内某重大危险源一旦发生火灾、爆炸，极易造成事故蔓延，影响邻近装置甚至会波及相邻企业，引发灾难性的多米诺效应事故。

国内外多米诺效应事故比普通事故造成了根据严重的危害和损失，如：1984年墨西哥城发生液化石油气槽车爆炸碎片引发连锁爆炸事故，烧毁建筑面积27公顷，造成544人死亡，1800多人受伤，120万人疏散[1]。2008年8月26日清晨6时45分，广西河池市宜州郊区最大的有机化工企业广西维尼纶集团有限责任公司，有机车间突然发生爆炸，爆炸引发的火灾导致各仓库内装有甲醇、乙炔、醋酸乙烯、液氯等易燃易爆物品的储存罐连环爆炸，造成20人死亡，60人受伤，疏散工厂生活区的居民及附近15个自然屯的群众共计1.1万多人[2]。

由血淋淋的教训可见，多米诺效应的研究十分有必要。在欧盟颁布的塞韦索法令II[3-4](Directive 96/82/EC)中规定，为了完善风险控制、风险评价和预防措施程序，要求有重大危险源的相关企业提交的安全报告中必须要有厂内和厂外多米诺效应分析的内容，以及可能产生多米诺效应的企业之间要共享危险源信息等，这一规定显示出国外对多米诺事故的重视。目前在国内，关于重大危险源的多米诺效应方面研究主要在理论研究方面较多，但是因为化工园区里的企业在空间地域中分布随机性很大，一旦发生事故，单从理论模型有时很难为园区的企业进行实时定量的分析，通过基于Mapinfo的地理信息系统可以将模型与空间地域相结合，能够有效地对多米诺效应事故进行实时定量地分析，在救援决策支持中起到重要的作用。

2.多米诺效应概述

2.1 多米诺效应原理

由初始事故引起的2次或2次以上的事故从而造成严重后果的现象称为事故多米诺效应[5]。

多米诺事故发生的必要条件：

(1)事故所形成破坏效应的范围之内存在敏感装置；

(2)破坏效应强度足够大以致装置失效；

表1 火灾事故的多米诺效应的临界值表

(3)破坏的发生存在极大的可能性[6]。

多米诺效应模式：初始单元发生火灾或爆炸，所释放的能量作用于二级或三级单元，在满足上述条件的情况下导致2次事故，以此类推导致3次事故的发生，这种现象将继续传播直到下一个区域内没有危险源为止，多米诺效应模式见图1。

能够触发多米诺效应发生的事故类型主要有：①火灾；②爆炸；③毒物泄漏；④以上3种情况同时或交叉发生[7]。

2.2 定量评价多米诺事故风险的流程

V.Cozzani等人在前期研究的基础上，根据QRA的思想提出了定量评价多米诺事故风险的流程，如图2所示。从图2可看出，定量评价多米诺事故风险的关键技术为：临界值标准和扩展概率的确定及可能二次事故组合的发生概率和后果计算。

2.3 多米诺效应临界值的确定

目前多米诺效应研究中主要考虑了热辐射、超压和抛射破片这三种扩展向量对二次目标设备的损坏。但是笔者目前主要针对热辐射进行基于MapInfo的定量分析，其余两种基于MapInfo的定量分析将在后续的研究中讨论分析。

从图2看出，多米诺效应一般利用临界值标准来辨识初始事故是否可能引起周围设备的损坏。如果初始事故在给定位置处的物理量值大于设定的临界值，则初始事故就有可能对该处的设备造成一定程度的损坏，有必要进一步对其产生的多米诺效应进行深入分析，否则认为初始事故不会对该处的设备造成影响，没有必要再分析，因此临界值标准在多米诺效应的定量评价中占据着非常重要的地位。表1总结了火灾事故的多米诺效应的参考临界值[8]。

2.4 多米诺效应热辐射概率分析

当二级单元吸收源于一级单元事故释放出的热辐射能量后，会引起单元内的设备温度升高、压力增大和设备材质的许用应力等物理特性的变化。一旦设备压力超过它的许用压力或者器壁所受应力超过其许用应力就有发生爆炸或破裂的可能性。

热辐射多米诺效应的概率模型，即“失效(ttf)”及ttf VS．热辐射简化模型的概率模型[9]：

常压圆柱形立式容器：

高压圆柱形水平容器：

式中：

tff：设备失效时间，s；

I：作用于目标设备的热辐射强度，kW/M2；

V：设备体积或容积，M3；

由上可知，热辐射多米诺效应的概率计算是以热辐射强度和设备体积为基础的。

3.多米诺事故热辐射叠加算法分析

3.1 热辐射强度的叠加

(1)如图3所示，假设在A处发生一级池火事故为例(不失一般性，其他火灾或爆炸事故公式不再一一列举)。可以通过以下公式算出池火的热辐射通量：

设液池为一半径为r的圆形池，则液池燃烧时放出的总热通量Q为：

式中，

r—液池半径，m；

h—火焰高度，m；

η—效率因子，可取0.13-0.35；

HC—液体燃烧热，J/kg；

假设全部辐射热都是从液池中心点的一个微小的球面发出的则在距液池中心某一距离的入射热辐射强度I为：

式中，Q—总热辐射通量，W；

tc—空气导热系数；

x—对象点到液池中心距离。

(2)通过公式(2-1)和(2-2)可以算出池火事故热辐射强度，并根据危险程度划分若干个等级的危险区域，如图3所示A1区、A2区、A3区，各区的热辐射强度均值分别为IA1、IA2、IA3。假设B点有一个危险源在A3区的影响范围之内，可以通过公式(1-1)或(1-2)算出该点的多米诺效应的概率，其他各点也可以一一算出其概率。

(3)若B点因热辐射影响而发生多米诺效应产生二级事故，通过同上原理，可以再为B点再划分若干个等级的危险区域，如图所示B1区、B2区、B3区，各区的热辐射强度均值分别为IB1、IB2、IB3。这样必然会产生A点和B点热辐射强度的叠加区域1区、2区、3区、4区，叠加效果如图4所示，叠加区域的辐射强度如下：

叠加1区辐射强度：I1=IA2+IB3

叠加2区辐射强度：I2=IA3+IB3

叠加3区辐射强度：I3=IA3+IB2

叠加4区辐射强度：I4=IA3+IB1

如若C处有危险源，C处位于叠加2区范围内，分别受到A3区和B3区的辐射强度的叠加为I2，通过公式(1-1)或(1-2)算出该点的多米诺效应的概率。其他处危险源概率类似。

3.2 热辐射区域图元的叠加算法

(1)以图4为例分析，B处事故为A处事故发生多米诺效应的下级事故，需判断新发生B事故各等级危险区域与当前已发生的A事故各等级危险区域叠加，如下：

①判断B3环形图元是否与A1图元相交，不相交；

②判断B3环形图元是否与A2环形图元相交，结果为相交，交集为1号图元；

③执行B3环形图元减去1号交集图元操作，剩余部分图元判断是否与A3环形图元相交，结果为相交，交集为2号图元；

④判断B2环形图元是否与A1图元相交，不相交；

⑤判断B2环形图元是否与A2环形图元相交，不相交；

⑥判断B2环形图元是否与A3环形图元相交，结果为相交，交集为3号图元；

⑦判断B3图元是否与A1图元相交，不相交；

⑧判断B3图元是否与A2环形图元相交，不相交；

⑨判断B3图元是否与A3环形图元相交，结果为相交，交集为4号图元；

(2)由以上分析可以归纳出热辐射区域图元的叠加算法如图5所示，其中：

Feature(i)为多米诺效应中最近发生事故热辐射的某个等级的区域；

Feature(j)为已发生事故的热辐射的某个区域或叠加区域；

FeatureTemp为临时图元；

FeatureTemp∩Feature(j)为两个图元的交集；

4.实例研究

4.1 基于Mapinfo的多米诺事故热辐射分析软件设计流程

由于重大事故多米诺效应的分析计算较为复杂，涉及到众多的影响因素和相关理论，为便于问题的解决，实际研究应在理论分析的基础上加以合理简化。利用Delphi7.0+Mapinfo开发了热辐射事故多米诺效应分析软件，主要程序流程如图6所示。

4.2 图层算法设计

4.2.1 地图图层设计

在Mapinfo地图的图层中增加一个temp图层如图7所示，在将整个热辐射区算法中的所用到的图元至于其中。在该图层中建立如图8所示表结构，ID为图元的编号，ColorCol为图元颜色的编号，RValues为该图元热辐射值。

4.2.2 叠加图元分割设计算法实现

本程序中最核心内容的为叠加图元的分割算法，利用Delphi7.0+Mapinfo算法实现如下

5.结论

目前，国内外对多米诺效应事故热辐射叠加尤其是定量方面的研究尚处于初步阶段，对于破片和超压叠加的研究仍需进一步深入进行。笔者在国内外最新研究的基础上，对基于MapInfo的多米诺事故热辐射叠加进行了分析与研究，并对多米诺事故热辐射叠加算法进行了模拟分析，开发了多米诺事故热辐射叠加风险评估软件。

由于我国在重大危险源的监管中很少考虑多米诺效应，又鉴于多米诺效应难以动态预测和控制，基于MapInfo的多米诺事故热辐射叠加分析软件可以更加直观形象地反映事故的状态，并根据发生概率动态预测事故进一步发展的趋势，为救援指挥提供决策．笔者所给分析方法或可为相关安全评价提供一定的参考，对相关园区的安全管理和救援决策具有重要意义。

[1]C.M.Pietersen.Analysis 0f the LPG Disaster in Mexico City.Journal of Hazardous Materials.1988,20:85-107

[2]罗文胜.宜州化工厂突发连环爆炸[J].财经,2008,18:30-31.

[3]V.Cozzani,E.Salzano.The quantitative assessment of domino effects caused by overpressure Part I.Probit models[J].Journal of hazardous Materials.2004,A107:67-80.

[4]V.Cozzani,E.Salzano.The quantitative assessment of domino effects caused by overpressure Part II.Case studies[J].Journal of hazardous Materials.2004,A107:81-94.

[5]Lees F P.Loss prevention in process industries.Butterworths,l-3,London,1996.

[6]Latha P,Gautem G,Raghavan KV.Strategies For quantif i cation of thermally initiated cascade effects.Loss Prevention Industries 1992,5(1):15-21.

[7]Khan F I,Abbasi S A.Estamation of Probabiliyies and likely,consequences of a chain of accidents (domino effect )in manali industrial complex.Journal of Cleaner Production,2001(9):493-508．

[8]V.Cozzani,G.Antonioni,E.Salzano. Quantitative assessment of domino scenarios by a GIS-based software tool[J].Journal 0f Loss Prevention in the Process industries.2006,19:468-477.

[9]V.Cozzani,G.Gubinelli,E.Salzano,etal.The assessment of risk caused by domino effect in quantitative area risk analysis.Journal of Hazardous Materials.2005,A127:14-30.

刘俊（1979—），男，江苏江都人，大学本科，常州大学信息科学与工程学院讲师，研究方向：计算机技术。