基于模糊逻辑的蒸气云爆炸参数预测框架算法及其应用

甄兰兰

(上海电力学院 自动化工程学院, 上海 200090)

蒸气云爆炸(Vapor Cloud Explosion,VCE)是工业厂区存在的主要安全危险之一。近半个世纪以来,国内外已发生多起由VCE引起的重大事故,造成严重的人员伤亡和财产损失[1-2]。

对VCE参数预测方法的研究,国内外专家学者提出了诸多方法,如三硝基甲苯(TNT)当量模型、TNO(The Netherlands Organization)、TNO多能法、沸腾液体扩展蒸气爆炸 (Boiling Liquid Expanding Vapor Explosion,BLEVE)、计算流体动力学(Computational Fluid Dynamics,CFD)等。其中,TNO多能法因其使用简便而被广泛应用[3-8],也被写入中华人民共和国安全生产行业标准《化工企业定量风险评价导则》(AQ/T 3046—2013)等相关标准中。TNO多能法中,爆炸强度被划分为10个等级 (标识为1～10),“1”代表强度最弱,“10”代表强度最强[7]。针对不同强度等级的爆炸,TNO多能法定义了对应的预测静态侧向压力曲线、动态侧向压力曲线、爆炸波持续时间及波形。爆炸强度等级确定的合理性,直接决定了爆炸参数的预测是否合理。若爆炸强度等级设定过高,则预测出的爆炸参数偏高,导致基于此预测结果所开展的工程设计偏于保守,造成不必要的建设投资浪费;反之,若爆炸强度等级设定过低,则预测出的爆炸参数偏低,导致基于此预测结果所开展的工程设计存在系统安全风险。

KINSELLA K G[10]提出了VCE超压参数的快速预测方法,根据专家经验对点燃能强度、障碍物阻塞程度、受限程度进行定性判定,再结合其提出的爆炸强度等级索引表,确定爆炸强度等级。

本文引入模糊逻辑方法,定量定性地确定点燃能强度、障碍物阻塞程度、受限程度,以更方便有效地使用专家经验、爆炸历史数据等知识,为此提出了基于模糊逻辑的蒸气云爆炸参数预测框架算法(Fuzzy Logic Prediction Framework Algorithm for Vapor Cloud Explosion Parameters,VCEP-FLPFA),合理确定爆炸强度等级及其隶属度,并基于此信息合理预测静态侧向压力、动态侧向压力、爆炸波持续时间及波形,以及人员伤害半径等相关信息。基于文献数据及经验信息,本文给出了预测框架方法的一种具体实现方法,并以文献[11-12]的LPG储配中心为例,说明预测方法的具体应用。

1 VCEP-FLPFA描述

模糊逻辑是通过模仿人的思维方式来表示和分析不确定、不精确信息的方法和工具,用于表达界限不清晰的定性知识与经验。它借助于隶属度函数概念,区分模糊集合,处理模糊关系,模拟人脑实施规则型推理。模糊推理过程包括输入模糊化、模糊规则库、模糊推理、解模糊化输出等[13]。

1.1 系统输入及模糊化

本文参照文献[10]中KINSELLA K G的定义,将点燃能强度、障碍物阻塞程度、受限程度的模糊变量定义如下。

1.2 输出模糊化

1.3 模糊规则库

KINSELLA K G提出的爆炸强度等级索引表如表1所示。

表1 爆炸强度等级索引

1.4 模糊推理

模糊推理定义为

1.5 解模糊化

2 VCEP-FLPFA的一种实现方法

2.1 定义隶属度函数

2.1.1 点燃能强度的隶属度函数

其隶属度函数为

(1)

(2)

(3)

点燃能强度的隶属度函数曲线如图1所示。

图1 点燃能强度的隶属度函数曲线

2.1.2 障碍物阻塞程度的隶属度函数

此处定义障碍物阻塞程度的隶属度为障碍物阻塞率隶属度与间距隶属度的乘积,记为

J=1,2,3

(4)

障碍物阻塞率和障碍物间距的隶属度函数曲线如图2和图3所示。

图2 障碍物阻塞率的隶属度函数曲线

图3 障碍物间距的隶属度函数曲线

2.1.3 受限程度的隶属度函数

(9)

(10)

(11)

受限程度的隶属度函数曲线如图4所示。

图4 受限程度的隶属度函数曲线

2.2 定义模糊系统输出

模糊系统输出为爆炸强度等级,此处定义爆炸强度等级为精确量,即为L=1,2,3,…,10。

2.3 定义模糊规则库

本文定义的模糊规则库如表2所示。

表2 模糊规则库

模糊规则库可描述为

式中:I=1,2,3;

J=1,2,3;

K=1,2;

M=1,2,3,…,18;

DM=1,2,3,…,10。

2.4 模糊推理

针对具体的输入IE,BR,BD,ER,CD,模糊推理结果为爆炸强度等级DM的可能性μ(DM)。其定义如下

(12)

2.5 解模糊化

现以预测爆炸静态侧向超压为例,说明具体的解模糊化方法。类似的方法也可以预测动态峰值压力和正向持续时间。

TNO多能法针对爆炸强度等级给出了对应无量纲距离γ′的无量纲静态侧向超压值P′S曲线,具体如图5所示[5]。

图5 无量纲静态侧向超压随无量纲距离的变化曲线

其中,无量纲距离的计算公式为

(13)

式中:z——感兴趣的地点与爆炸中心的距离,m;

P0——当地大气压,一般取1.0×105Pa;

E——爆炸源能量,MJ。

由图5可以获得无量纲静态侧向超压值P′S,则静态侧向超压值的计算公式为

PS=P′SP0

(14)

文献[14]给出了无量纲静态侧向超压值的拟合函数关系式,记为P′SL(γ′),L=1,2,3,…,10。如爆炸强度等级7的拟合函数为

(15)

基于模糊逻辑的无量纲静态侧向超压预测值计算公式为

(16)

3 应用实例

以文献[11]的Mexico城附近San Juan Ixhuatepec区域LPG为例,并引用文献[12]的部分计算数据,说明VCEP-FLPFA的具体应用方法。LPG系统安装概况见文献[11]。

3.1 VCEP-FLPFA计算

3.1.1 计算障碍物体积和自由空间体积

将文献[11]中LPG系统的重要危险区域划分成4个区域[12],区域划分如图6所示。

图6 重要危险区域划分示意

参考文献[12],计算每个区域的体积。

区域1:V1=32×31.5×6=6 048 m3;

区域2:V2=77×65×6=30 030 m3;

区域3:V3=45×31.25×4.5=6 328 m3;

区域4:V4=65.5×45×19.5=57 476 m3;

总体积:V总=99 882 m3;

LPG储罐体积:VLPG=16 000 m3;

自由空间的体积:Vgr=V总-VLPG=83 882 m3。

3.1.2 计算蒸气云尺寸

LPG储罐的供气管路泄漏,泄漏的液化石油气流到储罐下面,蒸发形成的液化石油气与空气的混合物弥漫在储罐之间。假定蒸发液化石油气的重量为4 750 kg,在环境温度为15 ℃时,液化石油气在混合气体中的当量体积百分比为4%,当量密度为1.86 kg/m3。

蒸气云体积:Vc=(100/4)×(4 750/1.86)=63 844 m3。

爆炸半径:[(3/2π)×63 844]1/3=31.24 m。

由于自由空间的体积Vgr=83 882 m3大于蒸气云体积Vc,所以Vgr值减小为与Vc值相等,即Vgr=Vc=63 844 m3。从而,非障碍区域的蒸气云体积为零。

液化石油气与空气混合物的当量燃烧热为3.46 MJ/m3,从而爆炸源能量E=3.46×106×63 844=220.9 GJ。

3.1.3 计算爆炸强度等级的隶属度

系统输入参数如下:点燃源为弥漫在储罐下的液化石油气,即受限蒸气云,假定点燃能IE=1 000 MJ;障碍物阻塞率BR=16 000/99 882=16.02%;障碍物间距BD=1.5 m;VCE半径ER=31.24 m;限制平面之间的距离CD=2.0 m。

根据本文提出的VCEP-FLPFA,计算出隶属度排在前3位的爆炸强度等级分别是8,6,10,对应的隶属度为:μ(8)=0.728 3;μ(6)=0.277 2;μ(10)=0.026 8。

3.1.4 计算爆炸静态侧向超压值

假定感兴趣的区域距离爆炸中心z为385 m,则

代入文献[14]给出的无量纲静态侧向超压值的拟合函数关系式,得

0.095 51

从而,

3.2 传统的TNO计算

在建筑防爆设计中,为了确保安全,往往人为地确定爆炸强度等级为10,则

3.3 计算结果对比及分析

在VCEP-FLPFA中,引入专家经验,其预测的爆炸静态侧向超压值更为合理。与VCEP-FLPFA相比,传统的TNO计算方法预测的爆炸静态侧向超压值高出约25%,这将会导致建筑防爆设计过于保守,大大增加建筑防爆设施的成本。

4 结 语

本文在传统TNO 多能法的基础上,引入了基于模糊逻辑的人工智能技术手段,所提出的VCEP-FLPFA可以更有效地使用专家经验、爆炸历史数据等知识,预测结果更合理,可为建筑防爆规划与设计提供可靠的科学数据。本文提出的方法具有通用性,随着人们对VCE研究的深入,可进一步优化本文提出的隶属度函数、爆炸强度等级经验表,或者提出更为合理的隶属度函数、爆炸强度等级索引表,以更有效地预测VCE参数。