[0043] 如图1所示，本发明的具体步骤是：
[0044] 步骤1:获取基础数据，其中包括危险源种类、状态、存量、分布情况以及储罐类型。
[0045] 步骤2：化工园区网格化。
[0046] 当评估区域确定时，采用相同大小的网格将整个化工园区划分为如图 2所示的以笛卡尔坐标系(n×m)中的二维空间。其中园区内危化品储罐已标注。在图1中，储罐位置可以分别用坐标(2,5)和(5,9)表示。如果每个网格的长度为100米，则图1中任意两个网格之间的距离可以用它们的坐标表示如下：
[0047]
[0048] 式中，dij为网格i与网格j之间的距离(米)，(xi,yi)和(xi,yj)分别为网格i与网格j的坐标；
[0049] 步骤3：设备故障频率分析。
[0050] 化工园区引发初始事故的原因主要有外部原因、设备失效和人为原因。外部原因可能是自然灾害或多米诺效应。本发明采用了设备故障引起的初始事故和多米诺效应引起的二次事故进行风险评估。这是化工园区中最常见的事故模式。
[0051] 常见的设备失效总结如下：
[0052] 1)基本过程控制系统(BPCS)失效；
[0053] 2)振动引起的失效(例如，旋转设备)；
[0054] 3)保养/维修不足，包括替代不当建造材料所造成的失效；
[0055] 4)控制支持系统失效(例如，电力、仪表气源)。
[0056] 5)磨损，疲劳，或腐蚀所造成的容器或管道失效；
[0057] 6)设计，规范，或制造缺陷所造成的容器或管道失效；
[0058] 7)超压(例如，热膨胀，清理/吹洗管道)或负压(真空失稳)；
[0059] 8)高温(如火灾暴露，冷却损失)或低温造成脆性断裂(如自动制冷设备，低温环境)引起的失效。
[0060] 目前，设备失效频率分析最常用的方法是基于数据库。事故发生的频率主要由事故的数量和设备暴露的密度来计算，计算方法如下：
[0061]
[0062] 式中，fd为事故频率(年)；Na是事故数量(年)；Nd为装置数量；Td为设备的暴露时间(年)。
[0063] 主要失效数据库如下：
[0064] 1)TNO陆上定量风险评估导则中泄漏频率统计；
[0065] 2)海上设备可靠性数据库；
[0066] 3)美国化学工程师学会过程安全中心(CCPS)过程设备可靠性数据库；
[0067] 4)美国知名风险评估公司Exida数据库；
[0068] 5)英国HSE海上事故数据库；
[0069] 步骤4：基于多米诺效应的事故后果分析。
[0070] ③计算个人死亡概率
[0071] 事故后果可分为人员伤亡，经济损失和环境污染等。考虑经济损失及环境污染在物料特性、设备布置、资产分布等诸多方面存在很多不确定性，本发明仅从人员伤亡角度来量化事故后果。
[0072] 不同损伤因素下的人体损伤概率单位计算公式如下：
[0073]
[0074] 式中，Yh为人体损伤概率单位；I为热辐射的强度(千瓦每平方米)；C为有毒气体浓度(ppm)；Ps为静态超压峰值(千帕)；te为人体暴露时间(分钟)。
[0075] 个体死亡概率计算公式如下：
[0076]
[0077] 式中，V为个体死亡概率，V∈[0,1]；u为积分变量。
[0078] ④多米诺效应分析
[0079] 多米诺效应定义为：始于一个单元，通过热辐射、冲击波超压、抛射碎片等影响临近单元设备，从而加重了后果的事件。引起多米诺效应的事故类型可分为火灾和爆炸两类。其中，对装置的危害主要有强烈的热辐射、冲击波超压。多米诺效应的拓展方式如图3所示。
但是，并不是所有的火灾和爆炸都能引起多米诺效应，火灾和爆炸的能量必须达到一定的阈值。因此，基于多米诺效应的事故后果分析更关注第二次事故发生的概率，即多米诺效应概率。
[0080] 不同损伤因素下的设备损伤概率单位计算公式如下：
[0081]
[0082] ln1(tff)＝‑1.128ln(I)‑2.667×10‑5V+9.87
[0083] ln2(tff)＝‑0.97ln(I)‑8.835V0.032
[0084] 式中，Yd为设备损伤概率单位；ttf为设备失效所需时间(秒)；V为二次设备容积(立方米)；Ps为静态超压峰值(千帕)；I为热辐射的强度(千瓦每平方米)。
[0085] 多米诺效应概率计算如下：
[0086]
[0087] 式中，P为多米诺效应概率；x为积分变量；
[0088] 步骤5：计算个人风险，个人风险是指长期生活在无防护危险区的人受到伤害的频率。它是可接受的社会安全风险的最小单元。任何一个网格中心的个人风险的计算如下所示：
[0089]
[0090] 式中，IRi为网格i的中心处的个人风险值；fs为设备s的失效频率；Vsi为设备s发生事故对网格i的中心处人员的个体死亡概率；Pss′为设备s发生事故引发设备s′失效的多米诺效应概率；
[0091] 步骤6：计算社会风险，社会风险是指事故发生概率与事故造成的伤亡人数之间的关系。任何一个网格中心的社会风险的计算如下所示：
[0092]
[0093] 式中，SRi为网格i的中心处的社会风险值；ρi为网格i内的人口密度；
[0094] 步骤7：通过离散法和插值法得到网格内其余区域的个人风险值以及社会风险值。
