CN103914622B - 一种化学品泄漏快速预测预警应急响应决策方法 - Google Patents
一种化学品泄漏快速预测预警应急响应决策方法 Download PDFInfo
- Publication number
- CN103914622B CN103914622B CN201410136961.5A CN201410136961A CN103914622B CN 103914622 B CN103914622 B CN 103914622B CN 201410136961 A CN201410136961 A CN 201410136961A CN 103914622 B CN103914622 B CN 103914622B
- Authority
- CN
- China
- Prior art keywords
- leakage
- wind direction
- sensor
- scene
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000000126 substance Substances 0.000 title claims abstract description 17
- 238000012216 screening Methods 0.000 claims abstract description 32
- 238000012549 training Methods 0.000 claims abstract description 31
- 238000009792 diffusion process Methods 0.000 claims abstract description 28
- 238000004088 simulation Methods 0.000 claims abstract description 24
- 230000007935 neutral effect Effects 0.000 claims abstract description 23
- 230000001537 neural effect Effects 0.000 claims abstract description 14
- 230000000505 pernicious effect Effects 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000013528 artificial neural network Methods 0.000 claims description 19
- 238000005457 optimization Methods 0.000 claims description 13
- 238000007689 inspection Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000012502 risk assessment Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000000052 comparative effect Effects 0.000 claims description 2
- 231100000614 poison Toxicity 0.000 abstract description 4
- 230000007096 poisonous effect Effects 0.000 abstract description 4
- 238000003062 neural network model Methods 0.000 abstract 2
- 238000013524 data verification Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 38
- 230000009931 harmful effect Effects 0.000 description 16
- 231100000331 toxic Toxicity 0.000 description 13
- 230000002588 toxic effect Effects 0.000 description 13
- ZAMOUSCENKQFHK-UHFFFAOYSA-N Chlorine atom Chemical compound [Cl] ZAMOUSCENKQFHK-UHFFFAOYSA-N 0.000 description 5
- 229910052801 chlorine Inorganic materials 0.000 description 5
- 239000000460 chlorine Substances 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 230000001235 sensitizing effect Effects 0.000 description 4
- 239000002341 toxic gas Substances 0.000 description 4
- KZBUYRJDOAKODT-UHFFFAOYSA-N Chlorine Chemical compound ClCl KZBUYRJDOAKODT-UHFFFAOYSA-N 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000001739 density measurement Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000012852 risk material Substances 0.000 description 2
- YGYAWVDWMABLBF-UHFFFAOYSA-N Phosgene Chemical compound ClC(Cl)=O YGYAWVDWMABLBF-UHFFFAOYSA-N 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Alarm Systems (AREA)
- Emergency Alarm Devices (AREA)
Abstract
本发明一种化学品泄漏快速预测预警应急响应决策方法,将扩散模型模拟与神经网络和气体传感器***相结合,应用于工业园区有害气体泄漏的快速预警及辅助决策,包括:园区风险因素识别,用于识别可能发生的各类泄露事故;数值模拟,对所有可能发生的事故进行模拟,得到有害气体的影响范围;数据筛选,根据实际的传感器布局提取和重组数值模拟结果中有效部分;神经网络训练,利用筛选后的数据对特定神经网络模型进行训练,以获得针对特定工业园区及周边条件的模型参数,并使用冗余数据参数验证;传感器***与神经网络模型整合,将模型与传感器DCS结合起来。
Description
技术领域
本发明属于工业生产中的安全预警技术领域,特别涉及一种化学品泄漏快速预测预警应急响应决策方法。
背景技术
很多流程工业在生产过程中会使用或者产生一些对人体有害的有毒、有害气体(如氯气,光气等),这些工业园区一旦发生有毒有害气体泄漏事故,泄漏出的有毒有害气体可能会对周边一定范围内的人类造成严重的危害。在有毒有害气体泄漏事故发生时,泄漏的物质和大致泄漏位置能够比较容易地确定,但有害气体的泄漏量或者泄漏速率则很难在现场获得。在有限的时间利用有限的信息预测有毒气体的扩散趋势和影响范围是事故应急响应过程的重要环节。目前,传统的数值扩散预测模型,包括高斯烟羽模型,计算流体力学模型,以及一些较为成熟的扩散模拟软件都需要使用者提供详细的泄漏源泄露速率,并且需要一定的时间计算才能得到结果。因此,泄漏源的泄露速率或者泄漏量以及较长的计算时间限制了这些传统气体扩散模型在事故应急响应以及辅助决策过程中的应用,它们更多地被用于事故发生之后的调查分析。由于历史原因,我国相当一部分可能发生有毒有害气体泄露事故的工业园区周边3~5km范围内即有居民区存在,当发生严重有毒有害气体泄露事故时,这些有毒有害气体很可能扩散至工业园区界区之外,对居民区产生威胁。在没有快速有效地预测有毒有害气体扩散范围的方法时,决策者往往以最坏的情况考虑有害气体的覆盖范围,而最坏情况的预测结果往往意味着政府部门需要疏散几万甚至十几万人,这是非常不现实的。因此一种快速有效地预测有害气体扩散范围的方法对于事故应急响应和辅助决策具有重要意义。
处理有毒有害气体的工厂都会有针对性地在有毒有害气体储罐附近设置一个或多个有害气体传感器,用于监控是否有气体泄漏。这些气体传感器均与控制中心相连,提供有毒气体泄漏报警信息。但目前很多工业园区均没有完全发挥气体传感器的作用,事故发生时,应急响应人员只能通过传感器获得泄漏气体的种类,瞬时浓度,而不能通过这些信息快速地预测泄漏气体的分布情况,进而采取有效的控制措施,制定合理的疏散计划。
目前,世界上普遍使用的事故后果分析方法主要包括使用不同种类的数值模型(高斯烟羽模型,计算流体力学模型(CFD)以及商业模型PHAST,FLACS等对已经发生的事故进行重现,展示在当时的事故泄露条件下泄漏气体的影响范围(包括死亡区,重伤区和影响区),并研究大气扩散条件对泄漏气体扩散范围和浓度分布的影响。但使用数值模型分析的缺点在于必须知道泄漏源的源强(泄露速率)以及泄露形式(***泄露/孔径泄露等),结合气象参数和传质扩散方程进行模拟,模型复杂,计算耗时长,不能用于事故状态下的实时或者快速预测。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种化学品泄漏快速预测预警应急响应决策方法,对工业园区进行风险分析,泄露情景模拟,适量补充并优化园区现有的有毒有害气体传感器***,将园区的有毒气体报警***与气体扩散预测分析相结合,为事故应急响应决策提供技术支持。
为了实现上述目的,本发明采用的技术方案是:
一种化学品泄漏快速预测预警应急响应决策方法,包括风险因素识别、情景数值模拟、模拟结果筛选、神经网络训练以及训练结果与传感器***集成共五个阶段,其中:
风险因素识别包括:识别化工园区的风险因素,量化各个风险要素并根据识别的风险要素及其取值范围,合理取值并组合各种可能发生的泄露情景;
情景数值模拟阶段包括:将园区风险因素识别步骤中组成的所有可能发生的泄露情景进行模拟,以获得不同泄露情景下泄露气体的影响范围;
模拟结果筛选包括:
对一种包括4个气体传感器的布局方案中的关键参数进行优化,利用模拟得到的泄露气体浓度分布进行简单优化,得到最适传感器布局;
根据确定的传感器布局方案,提取在传感器布局可测量风向范围内的虚拟探测数据以及与之匹配的环境敏感点泄漏气体扩散数据,按照神经网络的训练逻辑准备训练数据和校验数据;
神经网络训练步骤包括:
建立用于函数拟合的前向神经网络,将准备好的训练数据作为神经网络训练的输入输出,计算网络参数;
将检验数据输入部分输入神经网络,结果部分与神经网络预测结果对比,评估预测精度;
神经网络与传感器和园区控制***集成步骤:
将各个风险源的模拟分析以及训练结果与风险源和环境敏感点的地理位置信息以数据库和调用程序的形式结合,并提供与传感器***的数据接口,实现从事故发生—传感器报警—模型快速预测—辅助决策的工作流程。
所述园区风险识别部分由以下步骤组成:
步骤一、利用风险评估识别风险因子及事故情景因素并判定可能发生的各类泄漏情景的必备要素,包括风险要素和事故情景要素。这些要素包括并不限于风险物质种类(有毒/易燃的气体或易挥发液体)、储量/用量、储存位置(空间坐标)、储存形式(压力,温度,设备)、风速、湿度、温度、大气稳定度、周边环境敏感点位置(空间坐标)等。
步骤二、分别确定各个风险因子以及事故情境因素的取值范围(针对连续变量),在各个因素的取值范围内分别以一定步长划分数值序列(步长不得超过该因素取值范围的10%),最后将不同的风险因子/事故情景因素的各个取值组合成大量的可能发生的泄露情景。如共有n种风险因素,各个风险因素有Ni种取值(Ni≥1,整数),则共有种泄露情景。
所述情景数值模拟阶段使用的数值模型包括:高斯及类高斯扩散模型(GaussianPlume Model)、计算流体力学(CFD)模型以及PHAST、FLACS、SLAB、ALOHA、HGSYSTEM等整合模型。最后得到J种泄露情景模拟结果。
所述模拟结果筛选步骤需要对一组4个传感器进行布局优化,且4个传感器采用正多边形结构,顺序编号为1~4,相邻编号传感器的间距相等,优化步骤如下:
步骤一、根据泄漏源的位置以及园区布局,传感器的检测限,确定基本布局参数L1,d和α的取值范围。其中L1为1号传感器距离泄漏源的距离,d为传感器间距,α为1-2和1-4传感器夹角的1/2.
步骤二、根据风频信息确定主导风向,传感器对称布置在主导风向下风向上,实际风向与主导风向夹角为β。并取一个较大值β0做为优化区间,将其离散为实际风向序列(β1,β2,…,βM),序列长度为M。(一般取β0=50°,将区间[-50°,50°]离散为间隔1°的风向序列)
步骤三、对布局参数L1,d,α在取值范围内离散成参数序列,长度分别为P1,P2,P3,组成种布局方案。
步骤四、对第k种布局方案(k=1,2,…,K)使用数据筛选算法,针对第j种情景模拟结果(j=1,2,…,J),筛选出该布局方案在风向为βm(m=1,2,…,M)时的有效布局数目
步骤五、对第k种布局方案(k=1,2,…,K)计算不同风向βm(m=1,2,…,M)时的情景适用率记录满足m=1,2,,…,M的最大风向范围Wk,取满足{max(Wk),k=1,2,…,K}的布局为最优化布局。
所述模拟结果筛选步骤中判断有效布局的标准为:对于任意传感器布局,对特定风向下的每种泄漏情景,4个传感器中有3个或3个以上被烟羽覆盖则称该布局为在该风向和泄露情景条件下的有效布局。
所属模拟结果筛选步骤中,情景适用率为(m=1,2,…,Mk=1,2,…,K)。
所述模拟结果数据筛选部分包含以下筛选算法:
步骤一、检查传感器布局参数L1,d,α,实际风向β,下风向敏感源位置(L2,Y2)以及泄露情景序号J等是否符合要求
步骤二、对第j个泄露情景(j=1,2,…,J)模拟结果的气体扩散距离,浓度分布,扩散时间和烟羽宽度进行插值,确定筛选起始点
步骤三、计算参数γ,θ2,θ4,W1,和W2,其中γ为2,4号传感器以泄漏源为顶点和主导风向夹角,θ2,θ4为2、4号传感器以泄漏源为顶点和实际风向的夹角,W1为环境敏感点距离烟羽中心线的垂直距离,W2为环境敏感点处的烟羽半宽
步骤四、检查传感器布局在该泄露情景和实际风向下是否为有效布局,若是,继续;若不是,跳过并记录此泄露情景,继续步骤二
步骤五、提取3个有效传感器位置的污染源浓度和扩散时间。
步骤六、检查泄露气体是否能够扩散至下风向环境敏感点,即计算环境敏感点在第j个泄露情景中的位置参数Qj。如果能,继续步骤七,如果不能,记录此序号并继续步骤二。
步骤七、提取下风向环境敏感点处有害气体的浓度和扩散时间。
步骤八、将所有提取到的传感器位置浓度和时间值以及下风向敏感点浓度和时间值整理并保存。
所述的判断泄漏气体是否能扩散至下风向环境敏感点的方法为:计算环境敏感点在气体泄露烟羽中的位置参数Qj=W1/W2,若Qj≤1,则环境敏感点受到泄漏气体扩散影响,反之则不受影响。
所述神经网络训练部分实现过程如下:
步骤一、确定网络输入,以现场可得的参数作为输入,包括储存压力(p),风速(v),风向(β),温度(T),湿度(h),大气稳定度(S),风险源下风向3台气体传感器的有效浓度测量值(c1,c2,c3)和报警时间(t1,t2,t3),下风向敏感点的位置(L2,Y2)
步骤二、确定网络输出,以下风向环境敏感点处有害气体到达的浓度(c)和时间(t)作为输出。
步骤三、网络内部结构确定,使用前向神经网络(BP网络),包含一个非线性隐含层(sigmoid传递函数)和一个线性输出层(linear传递函数),网络不包含反馈回路。
步骤四、按照输入输出要求,使用数据筛选算法准备输入和输出数据矩阵以及冗余的验证数据。(矩阵的列为输入输出各要素,行为不同情景的输入输出数据)
步骤五、使用MATLAB神经网络工具箱,按照单次统一训练方法对神经网络使用输入输出数据训练,得到网络参数。
步骤六、使用冗余输入数据作为神经网络输入,比较网络输出和冗余输出的差异,评估预测精度。
所述神经网络与传感器***集成部分包括如下步骤:
步骤一、建立单一风险源与主导风向下风向的环境敏感点的区域地图。
步骤二、对某单一风险源进行风险分析,识别可能的泄露情景并对泄露情景进行模拟。
步骤三、执行数据筛选步骤的传感器布局优化,获得该传感器布局在主导风向下风向的最大适用角度。
步骤四、对该风险源主导风向下风向最大使用角度范围内的所有环境敏感点进行数据筛选,得到一组神经网络训练的输入输出矩阵。
步骤五、使用数据筛选得到的神经网络输入输出矩阵进行训练和预测精度评估。
步骤六、将训练得到的神经网络参数制作成能够自动调用的应用程序,使传感器的DCS数据能够通过数据接口采集并使用应用程序计算。
步骤七、真实事故发生时,传感器报警—启动预制好的应用程序—预测泄露扩散对环境敏感点的影响(包括是否影响环境敏感点或者泄漏气体到达环境敏感点的时间和浓度)—决策人员判断是否需要疏散。
与现有技术相比,本发明提出了一种处理有毒有害气体工厂的传感器布局的简单优化方法,并引进使用前期风险识别,数值模拟以及神经网络训练的快速预测方法,以便在事故发生时能够及时有效地预测有毒有害气体的扩散范围,为园区周围居民区人员疏散和保护提供决策支持。在园区现有有毒有害气体传感器的基础上,针对园区生产工艺以及周边环境,完成对园区生产过程中风险识别,泄露情景模拟,完善传感器布局并进行神经网络训练,克服了有毒气体泄漏事故发生时报警信息不足以支持事故应急响应的不足,保证了方法的适用性和准确性。
附图说明
图1为本发明一个实施例方法的流程图。
图2为本发明一个实施例方法应用情景示意图。
图3为本发明一个实施例方法得到的传感器布局优化曲线。
图4为本发明一个实施例方法用于数据筛选的基本算法
图5为本发明一个实施例方法应用的神经网络结构图。
图6为本发明一个实施例方法应用的神经网络冗余数据验证。
图7为本发明一个实施例方法应用的神经网络对敏感区域是否受影响的判断。
具体实施方式
下面结合附图和实施例详细说明本发明的实施方式。
图1是根据本发明的一个实施例的方法的流程图,根据本发明的一个实施例的方法包括:
—风险因素识别步骤
利用风险评估识别风险因子及事故情景因素并判定可能发生的各类泄漏情景的必备要素,包括风险要素和事故情景要素。这些要素包括并不限于风险物质种类(有毒/易燃的气体或易挥发液体)、储量/用量、储存位置(空间坐标)、储存形式(压力,温度,设备)、风速、湿度、温度、大气稳定度、周边环境敏感点位置(空间坐标)等。然后分别确定各个风险因子以及事故情境因素的取值范围(针对连续变量),在各个因素的取值范围内分别以一定步长划分数值序列(步长不得超过该因素取值范围的10%),最后将不同的风险因子/事故情景因素的各个取值组合成大量的可能发生的泄露情景。如共有n种风险因素,各个风险因素有Ni种取值(Ni≥1,整数),则共有种泄露情景。
—泄露情景模拟步骤
该步骤使用传统气体扩散模型,对前一步骤中识别出的各类泄露情景进行数值模拟,取得不同情况下的有毒有害气体扩散数据及分布情况。可用的模型包括高斯及类高斯扩散模型(Gaussian Plume Model)、计算流体力学(CFD)模型以及PHAST、FLACS、SLAB、ALOHA、HGSYSTEM等。
—有效数据筛选步骤
图2所示为有效数据筛选步骤的风险源—烟羽—传感器—下风向敏感点的位置关系图。模拟结果筛选步骤需要对一组4个传感器进行布局优化,且4个传感器采用正多边形结构,顺序编号为1~4,相邻编号传感器的间距相等。优化步骤如下:
步骤一、根据泄漏源的位置以及园区布局,传感器的检测限,确定基本布局参数L1,d和α的取值范围。其中L1为1号传感器距离泄漏源的距离,d为传感器间距,α为1-2和1-4传感器夹角的1/2.
步骤二、根据风频信息确定主导风向,传感器对称布置在主导风向下风向上,实际风向与主导风向夹角为β。并取一个较大值β0做为优化区间,将其离散为实际风向序列(β1,β2,…,βM),序列长度为M。(一般取β0=50°,将区间[-50°,50°]离散为间隔1°的风向序列)
步骤三、对布局参数L1,d,α在取值范围内离散成参数序列,长度分别为P1,P2,P3,组成种布局方案。
步骤四、对第k种布局方案(k=1,2,…,K)使用数据筛选算法,针对第j种情景模拟结果(j=1,2,…,J),筛选出该布局方案在风向为βm(m=1,2,…,M)时的有效布局数目其中有效布局定义为:对于任意传感器布局,对特定风向下的每种泄漏情景,4个传感器中有3个或3个以上被烟羽覆盖(即能探测到烟羽浓度)的布局方案。
步骤五、对第k种布局方案(k=1,2,…,K)计算不同风向βm(m=1,2,…,M)时的情景适用率(m=1,2,…,M k=1,2,…,K),记录满足m=1,2,,…,M的最大风向范围Wk,取满足{max(Wk),k=1,2,…,K}的布局为最优化布局。
图3为由L1=[40,60,80]m,d=[20,30,40]m和α=[30°,45°,60°]的取值组成的27中布局方案的布局优化结果(图3d)。在情景使用率为95%的条件下,传感器布局方案L1=40m,d=20m和α=60°为最优化布局,该布局的最大风向适用范围为﹣20°~20°。确定传感器最优化布局之后,使用图4所示的筛选算法进行模拟结果筛选:
步骤一、检查传感器布局参数L1,d,α,实际风向β,下风向敏感源位置(L2,Y2)以及泄露情景序号j等是否符合要求
步骤二、对第j个泄露情景(j=1,2,…,J)模拟结果的气体扩散距离,浓度分布,扩散时间和烟羽宽度进行插值,确定筛选起始点
步骤三、计算参数γ,θ2,θ4,W1,和W2,其中γ为2,4号传感器以泄漏源为顶点和主导风向夹角,θ2,θ4为2、4号传感器以泄漏源为顶点和实际风向的夹角,W1为环境敏感点距离烟羽中心线的垂直距离,W2为环境敏感点处的烟羽半宽
步骤四、检查传感器布局在该泄露情景和实际风向下是否为有效布局(有效传感器数目大于3),若是,继续;若不是,跳过并记录此泄露情景,继续步骤二
步骤五、提取3个有效传感器位置的污染源浓度和扩散时间。
步骤六、检查泄露气体是否能够扩散至下风向环境敏感点,即计算环境敏感点在第j个泄露情景中的位置参数Qj=W1/W2。如果能,即Qj≤1,继续步骤七,如果不能,即Qj>1,记录此序号并继续步骤二。
步骤七、提取下风向环境敏感点处有害气体的浓度和扩散时间。
步骤八、将所有提取到的传感器位置浓度和时间值以及下风向敏感点浓度和时间值整理并保存。
—神经网络训练步骤
神经网络训练部分实现过程如下:
步骤一、确定网络输入,以现场可得的参数作为输入,包括储存压力(p),风速(v),风向(β),温度(T),湿度(h),大气稳定度(S),风险源下风向3台气体传感器的有效浓度测量值(c1,c2,c3)和报警时间(t1,t2,t3),下风向敏感点的位置(L2,Y2)
步骤二、确定网络输出,以下风向环境敏感点处有害气体到达的浓度(c)和时间(t)作为输出。
步骤三、网络内部结构确定,如图5所示,使用前向神经网络(BP网络),包含一个非线性隐含层(sigmoid传递函数)和一个线性输出层(linear传递函数),网络不包含反馈回路。
步骤四、按照输入输出要求,使用数据筛选算法准备输入和输出数据矩阵以及冗余的验证数据。(矩阵的列为输入输出各要素,行为不同情景的输入输出数据)
步骤五、使用MATLAB神经网络工具箱,按照单次统一训练方法对神经网络使用输入输出数据训练,得到网络参数。
步骤六、使用冗余输入数据作为神经网络输入,比较网络输出和冗余输出的差异,评估预测精度。如图6显示了一个氯气泄露实施例的神经网络预测情况。图中所有子图的纵坐标均为特定泄露情境下,根据3个传感器报警信息、天气条件以及下风向敏感区域的位置,使用神经网络对下风向敏感区域的
(a)氯气的浓度(ppm)
(b)氯气到达敏感位置的时间(s)
(c)敏感位置烟羽的半宽W2(m)
(d)敏感位置距离实际风向的垂直距离W1(m)
而图6a~6d的横坐标则是Phast模拟的结果,理论上,图中散点的拟合结果应为直线y=x,实际结果显示,只有氯气浓度的神经网络预测结果相比Phast模拟结果偏小(斜率约为0.93),氯气到达时间、W1和W2的预测结果都非常精确。
利用图6c和6d得到的W1和W2的值可以预测下风向敏感区域是否受氯气泄漏的威胁,如图7所示,横坐标为用到的泄漏情景的数量:
(a)神经网络浓度预测值与Phast模拟值的相对误差
(b)根据图6中神经网络预测的W1和W2计算出的W1/W2比值
(c)根据Phast模拟以及烟羽中选点得到的W1/W2比值
(d)将预测W1/W2和理论W1/W2值取整后相减得到的逻辑结果
根据(a)可以发现神经网络预测值与理论值相比的最大误差不超过30%,图7b和7c非常接近,其差别在7d中显示:在2780例泄露情景中,有2例泄露情景下风向敏感点理论上受到氯气泄露威胁,但神经网络预测的结果显示其没有受到威胁,称为“漏报”,漏报率为0.072%;另一方面,只有14例泄露情景中谢分享敏感点理论上没有受到威胁,但神经网络预测显示其受到了威胁,称为“误报”,误报率为0.504%。
步骤七、将训练得到的神经网络参数制作成能够自动调用的应用程序,使传感器的DCS数据能够通过数据接口采集并使用应用程序计算。
步骤八、真实事故发生时,传感器报警—启动预制好的应用程序—预测泄露扩散对环境敏感点的影响(包括是否影响环境敏感点或者泄漏气体到达环境敏感点的时间和浓度)—决策人员判断是否需要疏散。
Claims (9)
1.一种化学品泄漏快速预测预警应急响应决策方法,包括风险因素识别、情景数值模拟、模拟结果筛选、神经网络训练以及训练结果与传感器***集成共五个阶段:
风险因素识别包括:识别化工园区的风险因素,量化各个风险要素并根据识别的风险要素及其取值范围,合理取值并组合各种可能发生的泄露情景;
情景数值模拟阶段包括:将园区风险因素识别步骤中组成的所有可能发生的泄露情景进行模拟,以获得不同泄露情景下泄露气体的影响范围;
模拟结果筛选包括:
对一种包括4个气体传感器的布局方案中的关键参数进行优化,利用模拟得到的泄露气体浓度分布进行简单优化,得到最适传感器布局;
根据确定的传感器布局方案,提取在传感器布局可测量风向范围内的虚拟探测数据以及与之匹配的环境敏感点泄漏气体扩散数据,按照神经网络的训练逻辑准备训练数据和校验数据;
神经网络训练步骤包括:
建立用于函数拟合的前向神经网络,将准备好的训练数据作为神经网络训练的输入输出,计算网络参数;
将检验数据输入部分输入神经网络,结果部分与神经网络预测结果对比,评估预测精度;
神经网络与传感器和园区控制***集成步骤:
将各个风险源的模拟分析以及训练结果与风险源和环境敏感点的地理位置信息以数据库和调用程序的形式结合,并提供与传感器***的数据接口,实现从事故发生—传感器报警—模型快速预测—辅助决策的工作流程;
其特征在于,所述神经网络与传感器***集成部分包括如下步骤:
步骤一、建立单一风险源与主导风向下风向的环境敏感点的区域地图;
步骤二、对某单一风险源进行风险分析,识别可能的泄露情景并对泄露情景进行模拟;
步骤三、执行数据筛选步骤的传感器布局优化,获得该传感器布局在主导风向下风向的最大适用角度;
步骤四、对该风险源主导风向下风向最大使用角度范围内的所有环境敏感点进行数据筛选,得到一组神经网络训练的输入输出矩阵;
步骤五、使用数据筛选得到的神经网络输入输出矩阵进行训练和预测精度评估;
步骤六、将训练得到的神经网络参数制作成能够自动调用的应用程序,使传感器的DCS数据能够通过数据接口采集并使用应用程序计算;
步骤七、真实事故发生时,传感器报警—启动预制好的应用程序—预测泄露扩散对环境敏感点的影响—决策人员判断是否需要疏散。
2.根据权利要求1所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述园区风险因素识别部分由以下步骤组成:
步骤一、利用风险评估识别风险因子及事故情景因素并判定可能发生的各类泄漏情景的必备要素,包括风险要素和事故情景要素;
步骤二、分别确定各个风险因子以及事故情境因素的取值范围,在各个因素的取值范围内分别以一定步长划分数值序列,其中步长不得超过该因素取值范围的10%,最后将不同的风险因子/事故情景因素的各个取值组合成大量的可能发生的泄露情景,当共有n种风险因素时,各个风险因素有Ni种取值,其中Ni≥1,且Ni为整数,则共有种泄露情景。
3.根据权利要求2所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述情景数值模拟阶段使用的数值模型包括:高斯及类高斯扩散模型(Gaussian PlumeModel)、计算流体力学(CFD)模型以及PHAST、FLACS、SLAB、ALOHA以及HGSYSTEM整合模型,最后得到J种泄露情景模拟结果。
4.根据权利要求2所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述模拟结果筛选步骤对一组4个传感器进行布局优化,且4个传感器采用正多边形结构,顺序编号为1~4,相邻编号传感器的间距相等,优化步骤如下:
步骤一、根据泄漏源的位置以及园区布局,传感器的检测限,确定基本布局参数L1,d和α的取值范围,其中L1为1号传感器距离泄漏源的距离,d为传感器间距,α为1-2和1-4传感器夹角的1/2;
步骤二、根据风频信息确定主导风向,传感器对称布置在主导风向的下风向上,实际风向与主导风向夹角为β,并取一个较大值β0做为优化区间,将其离散为实际风向序列(β1,β2,…,βM),序列长度为M;
步骤三、对布局参数L1,d,α在取值范围内离散成参数序列,长度分别为P1,P2,P3,组成种布局方案;
步骤四、对第k种布局方案使用数据筛选算法,针对第j种情景模拟结果,筛选出该布局方案在风向为βm时的有效布局数目其中k=1,2,…,K;j=1,2,…,J;m=1,2,…,M;
步骤五、对第k种布局方案计算不同风向βm时的情景适用率记录满足的最大风向范围Wk,取满足{max(Wk),k=1,2,…,K}的布局为最优化布局。
5.根据权利要求4所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述模拟结果筛选步骤中判断有效布局的标准为:对于任意传感器布局,对特定风向下的每种泄漏情景,4个传感器中有3个以上被烟羽覆盖则称该布局为在该风向和泄露情景条件下的有效布局。
6.根据权利要求4或5所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所属模拟结果筛选步骤中,情景适用率为其中m=1,2,…,M;k=1,2,…,K。
7.根据权利要求6所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述模拟结果数据筛选部分包含以下筛选算法:
步骤一、检查传感器布局参数L1,d,α,实际风向β,下风向敏感源位置(L2,Y2)以及泄露情景序号J是否符合要求;
步骤二、对第j个泄露情景模拟结果的气体扩散距离,浓度分布,扩散时间和烟羽宽度进行插值,确定筛选起始点,其中j=1,2,…,J;
步骤三、计算参数γ,θ2,θ4,W1,和W2,其中γ为2号或4号传感器以泄漏源为顶点与主导风向的夹角,θ2为2号传感器以泄漏源为顶点与实际风向的夹角,θ4为4号传感器以泄漏源为顶点与实际风向的夹角,W1为环境敏感点距离烟羽中心线的垂直距离,W2为环境敏感点处的烟羽半宽;
步骤四、检查传感器布局在该泄露情景和实际风向下是否为有效布局,若是,继续;若不是,跳过并记录此泄露情景,继续步骤二;
步骤五、提取3个有效传感器位置的污染源浓度和扩散时间;
步骤六、检查泄露气体是否能够扩散至下风向环境敏感点,即计算环境敏感点在第j个泄露情景中的位置参数Qj,如果能,继续步骤七,如果不能,记录此序号并继续步骤二;
步骤七、提取下风向环境敏感点处有害气体的浓度和扩散时间;
步骤八、将所有提取到的传感器位置浓度和时间值以及下风向敏感点浓度和时间值整理并保存。
8.根据权利要求7所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述的判断泄漏气体是否能扩散至下风向环境敏感点的方法为:计算环境敏感点在气体泄露烟羽中的位置参数Qj=W1/W2,若Qj≤1,则环境敏感点受到泄漏气体扩散影响,反之则不受影响。
9.根据权利要求1所述的化学品泄漏快速预测预警应急响应决策方法,其特征在于,所述神经网络训练部分实现过程如下:
步骤一、确定网络输入,以现场可得的参数作为输入,包括储存压力(p),风速(v),风向(β),温度(T),湿度(h),大气稳定度(S),风险源下风向3台气体传感器的有效浓度测量值(c1,c2,c3)和报警时间(t1,t2,t3),下风向敏感点的位置(L2,Y2);
步骤二、确定网络输出,以下风向环境敏感点处有害气体到达的浓度(c)和时间(t)作为输出;
步骤三、网络内部结构确定,使用前向神经网络即BP网络,包含一个非线性隐含层和一个线性输出层,非线性隐含层取sigmoid传递函数,线性输出层取linear传递函数,网络不包含反馈回路;
步骤四、按照输入输出要求,使用数据筛选算法准备输入和输出数据矩阵以及冗余的验证数据;
步骤五、使用MATLAB神经网络工具箱,按照单次统一训练方法对神经网络使用输入输出数据训练,得到网络参数;
步骤六、使用冗余输入数据作为神经网络输入,比较网络输出和冗余输出的差异,评估预测精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410136961.5A CN103914622B (zh) | 2014-04-04 | 2014-04-04 | 一种化学品泄漏快速预测预警应急响应决策方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410136961.5A CN103914622B (zh) | 2014-04-04 | 2014-04-04 | 一种化学品泄漏快速预测预警应急响应决策方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103914622A CN103914622A (zh) | 2014-07-09 |
CN103914622B true CN103914622B (zh) | 2017-07-07 |
Family
ID=51040298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410136961.5A Active CN103914622B (zh) | 2014-04-04 | 2014-04-04 | 一种化学品泄漏快速预测预警应急响应决策方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103914622B (zh) |
Families Citing this family (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376361A (zh) * | 2014-10-15 | 2015-02-25 | 南京航空航天大学 | 一种基于bp神经网络算法的核事故源项反演方法 |
CN104346657A (zh) * | 2014-11-10 | 2015-02-11 | 张江华 | 一种用于危险化学品泄漏的优化疏散方法及*** |
CN104680025B (zh) * | 2015-03-12 | 2017-08-25 | 重庆科技学院 | 基于遗传算法极限学习机的抽油机参数优化方法 |
CN104680023B (zh) * | 2015-03-12 | 2017-08-29 | 重庆科技学院 | 基于多目标决策的抽油机参数优化方法 |
CN104792935B (zh) * | 2015-03-20 | 2017-03-15 | 北京华泰诺安探测技术有限公司 | 一种核生化威胁云分析*** |
CN104792468B (zh) * | 2015-04-07 | 2016-04-06 | 中国石油大学(华东) | 一种炼油装置气体检测报警仪优化布置方法 |
CN104880494A (zh) * | 2015-06-02 | 2015-09-02 | 大连市环境监测中心 | 环境监测终端及环境应急监测快速部署*** |
CN104950708B (zh) * | 2015-06-02 | 2018-04-13 | 大连市环境监测中心 | 一种空中环境监测终端及环境应急监测部署*** |
CN106267620B (zh) * | 2015-06-10 | 2019-06-28 | 张国基 | 化学泄漏处理车 |
WO2017086981A1 (en) * | 2015-11-19 | 2017-05-26 | Halliburton Energy Services, Inc. | Method and system for monitoring and predicting gas leak |
CN106770930B (zh) * | 2016-11-24 | 2019-04-16 | 北京蓝色星语科技有限公司 | 一种使用无人机进行危险品检测的方法及*** |
CN106644284B (zh) * | 2016-12-16 | 2019-07-12 | 北京伟瑞迪科技有限公司 | 一种用于化工厂装置的自动泄漏监测与修复*** |
US10663367B2 (en) | 2017-07-28 | 2020-05-26 | International Business Machines Corporation | Adaptive sensing for a gas leak detection |
CN107607671A (zh) * | 2017-08-24 | 2018-01-19 | 慧感(上海)物联网科技有限公司 | 一种气体泄漏源定位方法 |
CN108091111B (zh) * | 2017-11-20 | 2018-11-13 | 淮阴工学院 | 一种油罐车油气泄漏智能预警*** |
CN108596364B (zh) * | 2018-03-29 | 2021-09-21 | 杭州电子科技大学 | 一种化工园区重大危险源动态预警方法 |
CN108733881B (zh) * | 2018-04-13 | 2021-10-12 | 同济大学 | 应对公共建筑室内突发空气污染的处理方法 |
CN108733893A (zh) * | 2018-04-25 | 2018-11-02 | 同济大学 | 耦合深度学习方法的公共建筑突发污染的溯源 |
CN108932379A (zh) * | 2018-06-26 | 2018-12-04 | 山东睿益环境技术有限公司 | 一种人员疏散应急处置方法及其应用 |
CN110874807A (zh) * | 2018-09-03 | 2020-03-10 | 中国石油化工股份有限公司 | 吸入性毒物泄漏应急响应处置卡及其使用方法和用途 |
CN109324155B (zh) * | 2018-10-31 | 2022-03-01 | 中国石油天然气股份有限公司 | 炼化装置尾气污染物在线监测数据预警方法及装置 |
CN109737309B (zh) * | 2019-01-09 | 2021-05-04 | 慧感(上海)物联网科技有限公司 | 一种基于风险识别的风险源定漏溯源方法及其*** |
CN110046834B (zh) * | 2019-05-13 | 2022-06-24 | 苏州大学 | 工作场所有害气体健康风险定量评估***及其评估方法 |
CN110531030B (zh) * | 2019-08-16 | 2021-11-12 | 北京慧辰资道资讯股份有限公司 | 一种基于物联网大数据大气污染物分析的方法及装置 |
CN110472359B (zh) * | 2019-08-21 | 2023-06-23 | 东北大学 | 基于数值模拟技术的高温熔融铝液与水接触***预防预警方法 |
CN110674983A (zh) * | 2019-09-05 | 2020-01-10 | 辽宁工程技术大学 | 一种基于copula函数尾部关联分析的工作面瓦斯预警方法 |
EP4031894B1 (de) * | 2019-09-18 | 2023-06-07 | Bayer Aktiengesellschaft | Vorhersage von mrt-aufnahmen durch ein mittels überwachten lernens trainiertes vorhersagemodell |
CN111105086A (zh) * | 2019-12-16 | 2020-05-05 | 茂盟(上海)工程技术股份有限公司 | 一种冷热电负荷预测可视化布局***及方法 |
CN111144055B (zh) * | 2019-12-27 | 2024-02-02 | 苏州大学 | 城市环境中有毒重气泄漏浓度分布确定方法、装置及介质 |
CN111157680B (zh) * | 2019-12-31 | 2021-10-26 | 北京辰安科技股份有限公司 | 室内挥发性物质的泄漏溯源方法及装置 |
CN111127279B (zh) * | 2019-12-31 | 2023-07-21 | 中国科学院合肥物质科学研究院 | 一种基于监测实验数据库的核应急决策***、方法及存储介质 |
CN111583067A (zh) * | 2020-03-31 | 2020-08-25 | 西南交通大学 | 一种城市地下大空间施工安全预警与应急决策方法及*** |
CN113821989A (zh) * | 2020-06-18 | 2021-12-21 | 中国石油化工股份有限公司 | 应急处置仿真风险设备建模方法、装置及存储介质 |
CN112330082B (zh) * | 2020-08-13 | 2023-08-15 | 青岛科技大学 | 一种基于动态机理模型的智能定量风险评估方法 |
EP4195091A4 (en) * | 2020-09-11 | 2024-04-03 | Siemens Ltd. China | METHOD AND APPARATUS FOR PREDICTING A MODEL-BASED ENVIRONMENT AVAILABLE ON A TARGET, PROGRAM AND ELECTRONIC DEVICE |
CN112270502B (zh) * | 2020-11-17 | 2021-06-01 | 北京三维天地科技股份有限公司 | 基于人工智能技术的环境应急任务协同处置平台 |
WO2022126544A1 (zh) * | 2020-12-17 | 2022-06-23 | 西门子(中国)有限公司 | 确定污染源的方法、装置和计算机可读存储介质 |
CN113095552B (zh) * | 2021-03-29 | 2023-12-26 | 中国海洋石油集团有限公司 | 一种海洋平台泄漏可燃气云团体积预测方法及*** |
CN113159337A (zh) * | 2021-04-08 | 2021-07-23 | 柳州东风容泰化工股份有限公司 | 一种邻氯苯酚泄露应急处理方法和装置 |
CN113468815A (zh) * | 2021-07-12 | 2021-10-01 | 南京航空航天大学 | 一种基于高斯扩散模型和深度神经网络的液态危化品挥发浓度估算方法 |
CN113624921B (zh) * | 2021-08-04 | 2024-03-29 | 杭州谱育科技发展有限公司 | 多模式综合污染溯源方法 |
CN114023043A (zh) * | 2021-10-13 | 2022-02-08 | 厦门理工学院 | 一种泄露气体安全预警方法 |
CN114001867A (zh) * | 2021-10-19 | 2022-02-01 | 北京伟瑞迪科技有限公司 | 面向园区气体泄漏的监测及应对方法 |
CN114263855B (zh) * | 2021-11-19 | 2024-04-26 | 合肥工业大学 | 一种天然气运输管道泄漏的预测方法及其应用 |
CN114241723A (zh) * | 2021-12-30 | 2022-03-25 | 武汉市三藏科技有限责任公司 | 一种基于高斯大气扩散用的环境监测报警装置及方法 |
CN115019482B (zh) * | 2022-08-04 | 2022-11-11 | 武汉鼎业环保工程技术有限公司 | 基于大数据的危险气体巡回预警方法及装置 |
CN115600643B (zh) * | 2022-10-17 | 2023-06-09 | 中国科学技术大学 | 一种毒害气体快速预测方法及*** |
CN115994431A (zh) * | 2022-12-14 | 2023-04-21 | 中国石油大学(华东) | 基于风险评价的氢气探测器空间布置方案优化方法及*** |
CN116403165B (zh) * | 2023-06-05 | 2023-08-04 | 交通运输部水运科学研究所 | 一种危险化学品泄露应急处理方法、装置及可读存储介质 |
CN117290639B (zh) * | 2023-11-24 | 2024-02-27 | 深圳大学 | 相变储能罐中材料泄漏定位检测的方法及*** |
CN117789422A (zh) * | 2024-02-26 | 2024-03-29 | 江西依爱弘泰消防安全技术有限公司 | 一种可燃气体报警控制***及方法 |
CN117959902B (zh) * | 2024-04-01 | 2024-06-04 | 大连华邦化学有限公司 | 一种基于数据反馈的气体纯化功能预警***及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2317311A2 (en) * | 2009-09-11 | 2011-05-04 | SELEX Galileo Ltd | Sensing network and method |
CN102313793A (zh) * | 2011-09-07 | 2012-01-11 | 中国环境科学研究院 | 一种化工园区空气中多组分气体污染物的监控*** |
-
2014
- 2014-04-04 CN CN201410136961.5A patent/CN103914622B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2317311A2 (en) * | 2009-09-11 | 2011-05-04 | SELEX Galileo Ltd | Sensing network and method |
CN102313793A (zh) * | 2011-09-07 | 2012-01-11 | 中国环境科学研究院 | 一种化工园区空气中多组分气体污染物的监控*** |
Non-Patent Citations (2)
Title |
---|
工业园区事故风险评价研究;孙爱军;《http://wenku.***.com/view/e6e7bcf1551810a6f4248600.html?from=search》;20130724;第36-53页 * |
翟卓韬等.基于MapX的泄露事故救援决策***.《计算机与应用化学》.2007,第24卷(第4期),第498-502页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103914622A (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103914622B (zh) | 一种化学品泄漏快速预测预警应急响应决策方法 | |
Wu et al. | Smart detection of fire source in tunnel based on the numerical database and artificial intelligence | |
KR101885350B1 (ko) | 대기오염 통합 관리 시스템 | |
RU2680652C2 (ru) | Способ комплексного мониторинга природной среды | |
CN102967689B (zh) | 一种污染源监测布点方法 | |
KR101956188B1 (ko) | 유해화학물질의 다층적 확산분석정보 제공방법 | |
CN104750949B (zh) | 一种气体泄漏中毒个人风险定量评估及表征方法 | |
CN112633553B (zh) | 燃气管线-危化企业耦合隐患辨识与风险评估方法及*** | |
Mousavi et al. | A sensitivity analysis of parameters affecting the hydrogen release and dispersion using ANOVA method | |
KR20220058673A (ko) | 환경시설의 대기오염가스 확산을 모니터링하기 위한 가상물리-기반 3차원 확산 예측지도 작성 시스템 및 그 방법 | |
Vázquez-Román et al. | A CFD-based approach for gas detectors allocation | |
Efthimiou et al. | CFD-RANS prediction of the dispersion of a hazardous airborne material released during a real accident in an industrial environment | |
CN108280849A (zh) | 一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法 | |
Jeon et al. | Development of surrogate model using CFD and deep neural networks to optimize gas detector layout | |
Rad et al. | A risk-based methodology for optimum placement of flammable gas detectors within open process plants | |
Ma et al. | Location of contaminant emission source in atmosphere based on optimal correlated matching of concentration distribution | |
KR20180117026A (ko) | 산업단지 내 위험확산 방지 방법 | |
CN115906706A (zh) | 一种危化品罐区多米诺综合风险高精度分析方法 | |
Lee et al. | A risk-based sensor placement methodology | |
CN105510535B (zh) | 基于现场实验的化工园区气体传感器扇形优化部署方法 | |
Qi et al. | Flexible alarming mechanism of a general GDS deployment for explosive accidents caused by gas leakage | |
US20200117840A1 (en) | Injection of simulated sources in a system of networked sensors | |
Dong et al. | Optimization of sensor deployment sequences for hazardous gas leakage monitoring and source term estimation | |
Ngae et al. | Optimization of an urban monitoring network for emergency response applications: An approach for characterizing the source of hazardous releases | |
CN116384255B (zh) | 一种基于多源数据融合的园区危险态势感知方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |