一种综合方法，用于分配工程性和程序性安全措施，以减轻火灾引发的连锁反应效应，该方法基于两阶段优化技术

《Reliability Engineering & System Safety》：An integrated approach for allocating engineered and procedural safety measures to mitigate fire-induced domino effects based on two-stage optimization

针对储罐农场多米诺事故提出两阶段协同优化模型，结合图论分析事故链传递路径，利用事件序列图量化应急响应效果，实现工程安全措施与应急资源配置的协同优化，案例显示成本降低57%。

　　

丁龙|钱一超|江宏瑞|陈思宁|姬杰

中国科学技术大学火灾科学国家重点实验室，安徽省合肥市金寨路96号，230026

摘要

在储罐场中，多米诺效应事故是一种影响巨大但发生概率较低的事件，可能导致严重的后果。安全措施是减少多米诺效应事故后果的重要手段，合理分配安全措施可以产生显著效果。储罐场中的安全措施分为三类：固有安全、工程安全和程序安全。现有的安全措施分配方法仅针对单一类型的安全措施。本文提出了一种两阶段优化模型，该模型利用图论和事件序列图来协同分配工程安全和程序安全措施，从而在提高风险管理效果的同时节省安全措施的投资。该方法包括四个主要步骤：分析多米诺效应的脆弱性、优化安全屏障的分配、建立应急响应模型以及优化应急团队的分配。案例研究证明了该方法的适用性和有效性，结果显示，与仅考虑应急响应相比，协同分配安全屏障和应急响应可以将成本降低57%。这项研究为工程师提供了一种实用且易于使用的安全措施协同分配方法，以降低多米诺效应事故的风险。

引言

随着化学工业的快速发展，化工厂呈现出扩张和集中的趋势。当设施更加密集时，火灾和爆炸等事故容易引发多米诺效应，造成巨大损失。多米诺效应是指由初始事件引发的一系列事故，导致总体后果比初始事故更为严重[1]。在多米诺效应事故的统计中[2]，41.8%的多米诺效应事故发生在化学品储存区域，其中火灾占事故升级情景的52.4%，化学品储存区域的火灾是最常见的多米诺效应事故类型。2019年3月，美国休斯顿洲际终端公司的散装液体储存罐区发生火灾，12个储罐被烧毁，造成超过1.5亿美元的损失[3]。同月，中国天嘉宜化工公司的储罐区也发生火灾，三个储罐场全部被烧毁，损失超过19.86亿元人民币[4]。2005年12月，英国邦斯菲尔德油库发生火灾，摧毁了20多个大型油罐，造成2.5亿英镑的损失[5]。安全措施在预防事故和减少事故后果方面发挥着重要作用[6]。通过合理分配安全措施，可以显著降低次生事故的概率和多米诺效应事故的后果严重程度。如何在资源有限的情况下有效分配安全措施是一个需要关注的问题。

Bollinger和Crowl[7]将安全措施分为三类：固有安全、工程安全和程序安全。固有安全措施包括优化工厂设施布局和库存，主要在化工厂设计阶段使用[8]。例如，Jung等人[9]开发了一种基于混合整数非线性规划的方法来进行设施选址和布局优化。在Jung的方法[9]基础上，Ge等人[10]提出了一种旨在降低多米诺效应风险的工厂布局方案。Khakzad等人[11]提出了一种基于DBN的方法，旨在优化化工厂库存并降低多米诺效应事故的风险。工程安全措施分为主动安全屏障和被动安全屏障，被动安全屏障是指无需主动操作即可降低危险频率或后果的安全屏障，如防火堤和防火涂层；主动安全屏障是指在检测到设施处于异常状态时响应的过程安全屏障[12]，包括喷水灭火系统和紧急压力释放系统。例如，Landucci等人[13]利用保护层分析和事件树分析量化了安全屏障在降低多米诺效应风险方面的有效性。基于Landucci的研究[13]，Khakzad等人[14]应用贝叶斯网络比较了多种不同的安全屏障分配策略，以获得最优解决方案。Zhang等人[15]使用图论和非支配排序遗传算法II（NSGA-II）获得了安全屏障分配的最优策略，并指出随着投资的增加，安全屏障的效果呈边际递减。此外，模糊层次分析法[16]、图论[17]、Petri网模型[18]和蒙特卡洛模拟[19]也广泛用于安全屏障优化。程序安全措施主要指应急响应，在事故救援中，应急团队使用消防车、泡沫枪等消防设备来减少事故后果。Petri网广泛用于建模和分析应急响应过程，包括多火灾场景下的应急响应[20]、多部门协同救援[21]和应急资源部署[22]。Khakzad等人[23]提出了一种基于图的方法，用于确定优先抑制和冷却的储罐，使用外接近度指数来判断储罐的重要性并指导应急响应决策。

大多数以往的研究仅关注单一类型的安全措施，并分配这些措施来预防多米诺效应事故，忽略了多种类型安全措施的协同作用。多种类型安全措施的协同作用是指多种安全措施共同作用以降低多米诺效应事故的风险和后果。在资源有限的情况下，仅考虑单一类型的安全措施可能无法将火灾多米诺效应事故的风险和后果降至可接受的范围，从而导致次生事故的发生。Zhou等人[22]指出，在某些事故场景中，仅依赖程序安全措施无法显著降低事故后果，需要增加防火堤和喷水灭火系统。Khakzad等人[14]在其关于工程安全措施的研究中认为，需要添加其他类型的安全措施来降低多米诺效应事故的风险。因此，考虑多种类型安全措施的协同分配优化方法具有重要意义。与单一类型安全措施的分配相比，协同分配将在提高风险管理效果的同时节省安全措施的投资。

安全措施协同分配方法的难点在于它们不同的应用场景。固有安全措施主要应用于化工厂的设计阶段，工程安全措施应用于设计和运营阶段，程序安全措施应用于事故阶段。以往的研究主要集中在某一阶段的安全措施分配上。也有一些研究评估了工程安全和程序安全措施的综合效果。Daas等人[24,25]基于证据推理理论通过主观安全分析评估了安全屏障和应急响应的有效性。Yu等人[26]和Guo等人[27]利用模糊层次分析过程评估了安全屏障和应急响应的可靠性，并通过贝叶斯网络计算了事故升级概率。他们的研究分别在不同的阶段考虑了安全屏障和应急响应，Yu等人[26]指出，他们的多阶段评估结果比仅考虑单一阶段分析的传统方法更符合实际情况，这表明了本研究的必要性。这些研究侧重于评估安全措施的可靠性和有效性，而较少关注不同阶段的安全措施协同分配。因此，在进行协同分配时，考虑两个阶段的安全屏障和应急响应非常重要。

为了解决上述问题，本文提出了一种基于多阶段优化方法的多类型安全措施协同分配优化方法。多阶段分配优化方法是一种数学建模方法，将复杂的分配优化过程划分为多个有序阶段，该方法广泛应用于资源分配[28]、项目规划[29]、生产调度[30]等领域。由于本研究主要关注运营阶段和事故阶段的安全措施分配，而固有安全措施仅用于设计阶段，因此选择工程安全和程序安全措施作为研究对象，并使用两阶段优化方法建立两种类型的安全措施分配优化模型。在第一阶段，利用图论和遗传算法在固定预算和最坏事故情景下获得工程安全措施的最优分配策略。第二阶段是对事故阶段的程序安全措施进行分配，通过合理分配应急团队来最小化事故后果。为了优化应急响应的分配，使用事件序列图建模来计算救援过程中的热辐射变化。根据改进的probit模型计算不同数量应急团队下的储罐失效时间，并以最大升级时间为目标函数。通过遗传算法获得火灾场景下应急团队分配的最优策略。

本文的结构如下：第2节介绍该方法的理论基础，第3节提供了在遗传算法下安全屏障和应急响应协同分配的优化模型。第4节的案例研究展示了所提出方法的有效性。最后，第5节总结了本研究的工作。

节选内容

两阶段优化方法

两阶段优化方法旨在将复杂的优化问题分解为两个更简单的优化问题。该方法将整个过程分为两个阶段，每个阶段都有独立的决策变量和目标函数。第一阶段完成后，其结果作为第二阶段的输入或约束条件[31]。两阶段优化方法主要有两个特点：

(1)

顺序性。两个阶段

方法论

如图1所示，两阶段安全措施分配优化模型主要分为四个步骤：(1) 分析每个储罐的多米诺效应脆弱性。第一步是根据储罐场的布局、储罐的大小、天气等数据，计算储罐事故发生时其他储罐接收的热辐射。根据热辐射矩阵建立加权有向图，并可以计算每个储罐的中心性指数

应用场景

案例选自Zhou等人的研究[46]，该储罐场由六个配置和物质相同的常压储罐组成。储罐直径为10米，高度也为10米，每个储罐之间的距离为30米，储罐内储存的物质是苯。储罐场的布局如图5所示。

在这种情况下，周围没有遮蔽物，气温为25°C，相对湿度为25%，天空部分多云且无风。

结论

安全措施是降低储罐场多米诺效应事故风险的重要方法。本文提出了一种优化安全屏障和应急响应协同分配的方法，以防止多米诺效应事故。该方法主要考虑通过两阶段优化结合图论和事件序列图来协同分配安全屏障和应急响应。其目的是防止或延缓事故的升级

CRediT作者贡献声明

丁龙：方法论、数据整理、调查、撰写——审稿与编辑初稿、资金获取。钱一超：概念构思、撰写——原始初稿、监督。江宏瑞：撰写——审稿与编辑、监督。陈思宁：撰写——审稿与编辑、监督、概念构思。姬杰：监督、项目管理、资金获取。

CRediT作者贡献声明

丁龙：撰写——原始初稿、方法论、调查、资金获取、数据整理。钱一超：撰写——原始初稿、监督、概念构思。江宏瑞：撰写——审稿与编辑、监督。陈思宁：撰写——审稿与编辑、监督、概念构思。姬杰：监督、项目管理、资金获取。