CN109086539A - 液氨储罐泄露伤害半径获取方法和设备 - Google Patents
液氨储罐泄露伤害半径获取方法和设备 Download PDFInfo
- Publication number
- CN109086539A CN109086539A CN201810917779.1A CN201810917779A CN109086539A CN 109086539 A CN109086539 A CN 109086539A CN 201810917779 A CN201810917779 A CN 201810917779A CN 109086539 A CN109086539 A CN 109086539A
- Authority
- CN
- China
- Prior art keywords
- radius
- storage tank
- liquid ammonia
- leakage
- ammonia storage
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
本发明实施例提供一种液氨储罐泄露伤害半径获取方法和设备,包括:基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。本发明实施例提供的液氨储罐泄露伤害半径获取方法和设备,通过设置判断液氨储罐泄露形成的云团是否为重气云团,并基于判断结果,获取对应的伤害半径,能够将重气云团扩散和非重气云团的扩散有机结合起来,使得危险物质的扩散过程成为一个连续的变化过程,进而准确获取液氨储罐泄漏浓度扩散的伤害半径。
Description
技术领域
本发明实施例涉及液氮储罐泄露事故处理领域,更具体地,涉及一种液氨储罐泄露伤害半径获取方法和设备。
背景技术
现在,很多工业企业在生产过程中会使用液氨进行制冷,液氮储罐用途广泛,但是,一旦发生泄漏事故,泄漏出的氨气可能会对周边环境和居民造成严重的危害。当事故发生时的人员疏散是应对这种突发事故的有效方案之一,在人员疏散方案选择时,如何建立科学有效的泄漏浓度扩散分析方法,即如何获取液氨储罐泄露的伤害半径对于事故应急响应和辅助决策具有重要意义。
目前,液氨泄漏普遍采用湍流的高斯/重气扩散模型进行计算,获取在当时的事故泄露条件下泄漏气体的影响范围(包括死亡区,重伤区和影响区)。但是,这种不分阶段的伤害半径计算获取方法不够准确。
发明内容
针对现有技术存在的问题,本发明实施例提供一种液氨储罐泄露伤害半径获取方法和设备。
本发明实施例提供一种液氨储罐泄露伤害半径获取方法,包括:基于液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
本发明实施例提供一种液氨储罐泄露伤害半径获取设备,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中:所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行上述获取方法。
本发明实施例提供的液氨储罐泄露伤害半径获取方法和设备,通过设置判断液氨储罐泄露形成的云团是否为重气云团,并基于判断结果,获取对应的伤害半径,能够将重气云团扩散和非重气云团的扩散有机结合起来,使得危险物质的扩散过程成为一个连续的变化过程,进而准确获取液氨储罐泄漏浓度扩散的伤害半径。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明液氨储罐泄露伤害半径获取方法实施例的流程图;
图2为本发明另一液氨储罐泄露伤害半径获取方法实施例的流程图;
图3为本发明实施例中的液氨储罐泄露伤害半径获取设备的框架示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明液氨储罐泄露伤害半径获取方法实施例的流程图,如图1所示,包括:S101、基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;S102、基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
具体地,步骤S101中,所述液氨储罐为装有液氨的储罐。所述重气云团为由于重气效应在底面形成的薄而宽的云团。
进一步地,步骤S102中,所述基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径,具体包括:
若所述液氨储罐泄露形成的云团为重气云团,获取对应于重气云团的伤害半径;若所述液氨储罐泄露形成的云团为非重气云团,获取对应于非重气云团的伤害半径。
进一步地,获取了所述伤害半径后,本发明实施例中还包括基于所述伤害半径获取伤害范围,具体为基于泄露原点和所述伤害半径,获取对应的伤害范围。
本发明实施例提供的液氨储罐泄露伤害半径获取方法,通过设置判断液氨储罐泄露形成的云团是否为重气云团,并基于判断结果,获取对应的伤害半径,能够将重气云团扩散和非重气云团的扩散有机结合起来,使得危险物质的扩散过程成为一个连续的变化过程,进而准确获取液氨储罐泄漏浓度扩散的伤害半径。
需要说明的是,在进行本发明实施例所述的液氨储罐泄露伤害半径获取方法之前,还包括:获取液氨储罐泄露的工况参数,用以进行后续的计算。
所述工况参数包括:物质信息、气象信息、周边地貌、泄漏事故、事故场景、液位高度信息。其中,物质信息包括:储存温度、储存压力、储罐容积、附属管道内径、充装百分比;气象信息包括:风向、风速、大气稳定度;周边地貌:郊区、市区;可能泄漏设备包括:管道、阀门、罐体本身;事故场景:小孔泄漏、中孔泄漏、大孔泄漏、完全破裂、自定义事故场景。其中,小孔对应孔径小于等于5mm,中孔对应孔径为5-25mm,大孔对应孔径为25-100mm,完全破裂对应孔径大于100mm。
基于上述实施例,所述基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径,具体包括:若所述液氨储罐泄露形成的云团为重气云团,则基于液氨储罐泄露初始时刻重气云团内部浓度、预设伤害扩散浓度和液氨储罐瞬时泄漏形成的重气云团的初始体积,获取对应的待确定轻伤半径、待确定重伤半径和待确定死亡半径,其中,所述预设伤害扩散浓度包括预设轻伤扩散浓度、预设重伤扩散浓度和预设死亡扩散浓度;通过ε准则,获取所述重气云团让位于环境湍流的第一扩散距离;若所述第一扩散距离大于所述待确定轻伤半径,并且所述第一扩散距离大于所述待确定重伤半径,并且所述第一扩散距离大于所述待确定死亡半径,则将所述待确定轻伤半径确定为轻伤半径,将所述待确定重伤半径确定为重伤半径,将所述待确定死亡半径确定为死亡半径。
具体地,基于液氨储罐泄露初始时刻重气云团内部浓度、预设伤害扩散浓度和液氨储罐瞬时泄漏形成的重气云团的初始体积,获取对应的待确定轻伤半径、待确定重伤半径和待确定死亡半径,具体包括:
通过下式获取对应的待确定轻伤半径、待确定重伤半径和待确定死亡半径:
其中,x为对应的待确定轻伤半径、待确定重伤半径或待确定死亡半径;C为对应的预设轻伤扩散浓度、预设重伤扩散浓度或预设死亡扩散浓度;C0为液氨储罐泄露初始时刻重气云团内部浓度;Vo为液氨储罐瞬时泄漏形成的重气云团的初始体积。
进一步地,随着空气的不断卷入,重气云团的密度将不断减小,重气坍塌引起的重气云团的扩散将逐步让位于环境湍流引起的扩散,通过ε准则,通过下式获取所述重气云团让位于环境湍流的第一扩散距离:
其中,xf为第一扩散距离;g为重力加速度;Vo为液氨储罐瞬时泄漏形成的重气云团的初始体积;ρ0为云团的初始密度;ρa为环境空气密度;εcr=0.01。
进一步地,将所述待确定轻伤半径确定为轻伤半径,将所述待确定重伤半径确定为重伤半径,将所述待确定死亡半径确定为死亡半径之后还包括:基于泄露原点,获取对应的轻伤范围、重伤范围和死亡范围。
所述轻伤范围为当液氨储罐泄露时,可能受到轻伤的范围,所述重伤范围为当液氨储罐泄露时,可能受到重伤的范围,所述死亡范围为当液氨储罐泄露时,可能受到非常严重的伤害甚至失去生命的范围。
基于上述实施例,若所述第一扩散距离小于等于所述待确定轻伤半径,或者所述第一扩散距离小于等于所述待确定重伤半径,或者所述第一扩散距离小于等于所述待确定死亡半径,则通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径。
需要说明的是,所述第一扩散距离小于等于所述待确定轻伤半径,或者所述第一扩散距离小于等于所述待确定重伤半径,或者所述第一扩散距离小于等于所述待确定死亡半径,解释为第一扩散距离小于等于待确定轻伤半径、第一扩散距离小于等于待确定重伤半径和第一扩散距离小于等于待确定死亡半径三个条件满足其一即可进行后续的步骤。
具体地,随着空气的不断卷入,液氨储罐泄漏扩散的相关计算应分为二个阶段:初始以气雾状态时的扩散;然后随着空气的大量卷入,密度不断变小,逐渐变为湍流的中性气体。即随着空气的不断卷入,重气云团会逐渐扩散为非重气云团,扩散成的非重气云团基于环境湍流引发会产生第二扩散距离。
需要说明的是,形成重气云团,直至扩散成非重气云团是一个持续的过程。
本发明实施例提供的液氨储罐泄露伤害半径获取方法,不仅可以模拟出稳定状态下的液氨持续泄漏扩散态势,还可以动态的模拟出液氨持续泄漏扩散随时间变化的过程。
基于上述实施例,所述基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径,具体包括:若所述液氨储罐泄露形成的云团为非重气云团,则通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径。
基于上述实施例,所述基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团,之前还包括:获取所述液氨储罐泄露的闪蒸量。
在上述实施例的基础上,所述获取所述液氨储罐泄露的闪蒸量,具体包括:通过伯努利方程,获取所述液氨储罐的泄漏量;基于所述泄漏量,通过闪蒸计算,获取所述液氨储罐泄露的闪蒸量。
通过伯努利方程,通过下式获取所述液氨储罐的泄漏量:
其中,Q0为液氨储罐的泄漏量;Cd为液体泄漏系数;A为裂口面积;ρ为泄漏液体密度;P为初始罐压;P0为环境压力;g为重力加速度;h为裂口之上液位高度。
基于所述泄漏量,通过闪蒸计算,通过下式获取所述液氨储罐泄露的闪蒸量:
其中,Mv为闪蒸量;Cp为泄漏液体的定压比热;Q0为液氨储罐的泄漏量;t为泄漏时间;T为泄漏前液体储存温度;T0为液体在常压下的沸点;H为液体的气化热。
基于上述实施例,所述基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团,具体包括:基于所述液氨储罐泄露的闪蒸量,获取所述液氨储罐泄漏形成的云团的初始体积;基于所述液氨储罐泄漏形成的云团的初始体积、环境风速、重力加速度、云团的初始密度和环境空气密度,判断所述液氨储罐泄露形成的云团是否为重气云团。
具体地,所述若所述液氨储罐泄露形成的云团为重气云团,通过下述公式判断,具体为:
若下式成立,则所述液氨储罐泄露形成的云团为非重气云团,对应的,若下式不成立,则所述液氨储罐泄露形成的云团为重气云团,式为:
其中,Vo为液氨储罐泄漏形成的云团的初始体积;V为环境风速;g为重力加速度;ρ0为云团的初始密度;ρa为环境空气密度。
其中,Vo为液氨储罐泄漏形成的云团的初始体积;Mv为闪蒸量;M为氨的分子量;T0为液体常压下沸点。
基于上述实施例,所述通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径,具体包括:通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离;根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度;获取和所述预设轻伤扩散浓度的差值的绝对值最小的泄露最大浓度的第一位置,并且,获取和所述预设重伤扩散浓度的差值的绝对值最小的泄露最大浓度的第二位置,并且,获取和所述预设死亡扩散浓度的差值的绝对值最小的泄露最大浓度的第三位置;将所述第一位置和所述泄露原点之间的距离作为轻伤半径,并且,将所述第二位置和所述泄露原点之间的距离作为重伤半径,并且,将所述第三位置和所述泄露原点之间的距离作为死亡半径。
具体地,根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度,具体包括:
泄露最大浓度分为两种情况,第一种情况是液氨储罐完全破裂造成的瞬时泄漏,第二种情况是管道、储罐、阀门的大孔、中孔、小孔泄漏以及管道的完全破裂造成的连续泄漏。
处于第一种情况时,根据所述第二扩散距离,通过下式获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度:
其中,Q为排放的物料质量;σx、σy和σz为x、y、z方向扩散参数;Vt为瞬时排放时的体积;t为瞬时排放的时刻;x为下风向距离;y为横风向距离;z为离地面的距离。
处于第二种情况时,根据所述第二扩散距离,通过下式获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度:
其中,Q为排放的物料质量;σx、σy和σz为x、y、z方向扩散参数;x为下风向距离;y为横风向距离;z为离地面的距离。
进一步地,所述第一位置、所述第二位置和所述第三位置均为由x、y表示的位置坐标。
进一步地,预设轻伤扩散浓度优选设置为360mg/m3,预设重伤扩散浓度优选设置为1390mg/m3,预设死亡扩散浓度优选设置为3500mg/m3。
具体地,所述第一位置、所述第二位置和所述第三位置具体可通过下述步骤获取:
首先,设定初始的步长step。
其次,从所述本地坐标系的坐标(0,b0)开始,沿着x轴的正方向,按照所述步长step递增和递减取x值,根据高斯模型计算出不同泄漏时刻下氨气的预设轻伤扩散浓度(360mg/m3)、预设重伤扩散浓度(1390mg/m3)、预设死亡扩散浓度(3500mg/m3)对应的y值。
最后,将氨气浓度下计算出来的(x,y)值连起来就形成了对应的轻伤半径、重伤半径、死亡半径。
基于上述实施例,所述通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径,具体包括:通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离;根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度;获取和预设轻伤扩散浓度的差值的绝对值最小的泄露最大浓度的第一位置,并且,获取和预设重伤扩散浓度的差值的绝对值最小的泄露最大浓度的第二位置,并且,获取和预设死亡扩散浓度的差值的绝对值最小的泄露最大浓度的第三位置;将所述第一位置和所述泄露原点之间的距离作为轻伤半径,并且,将所述第二位置和所述泄露原点之间的距离作为重伤半径,并且,将所述第三位置和所述泄露原点之间的距离作为死亡半径。
具体地,根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度,具体包括:
泄露最大浓度分为两种情况,第一种情况是液氨储罐完全破裂造成的瞬时泄漏,第二种情况是管道、储罐、阀门的大孔、中孔、小孔泄漏以及管道的完全破裂造成的连续泄漏。
处于第一种情况时,根据所述第二扩散距离,通过下式获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度:
其中,Q为排放的物料质量;σx、σy和σz为x、y、z方向扩散参数;Vt为瞬时排放时的体积;t为瞬时排放的时刻;x为下风向距离;y为横风向距离;z为离地面的距离。
处于第二种情况时,根据所述第二扩散距离,通过下式获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度:
其中,Q为排放的物料质量;σx、σy和σz为x、y、z方向扩散参数;x为下风向距离;y为横风向距离;z为离地面的距离。
进一步地,所述第一位置、所述第二位置和所述第三位置均为由x、y表示的位置坐标。
进一步地,预设轻伤扩散浓度优选设置为360mg/m3,预设重伤扩散浓度优选设置为1390mg/m3,预设死亡扩散浓度优选设置为3500mg/m3。
具体地,所述第一位置、所述第二位置和所述第三位置具体可通过下述步骤获取:
首先,设定初始的步长step。
其次,从所述本地坐标系的坐标(0,b0)开始,沿着x轴的正方向,按照所述步长step递增和递减取x值,根据高斯模型计算出不同泄漏时刻下氨气的预设轻伤扩散浓度(360mg/m3)、预设重伤扩散浓度(1390mg/m3)、预设死亡扩散浓度(3500mg/m3)对应的y值。
最后,将氨气浓度下计算出来的(x,y)值连起来就形成了对应的轻伤半径、重伤半径、死亡半径。
作为一个优选的实施例,图2为本发明另一液氨储罐泄露伤害半径获取方法实施例的流程图,请参考图2,本发明实施例中所述的液氨储罐泄露伤害半径获取方法包括:
首先,获取液氨储罐泄露的工况参数。
其次,通过伯努利方程,获取所述液氨储罐的泄漏量;基于所述泄漏量,通过闪蒸计算,获取所述液氨储罐泄露的闪蒸量。
然后,基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团。
若所述液氨储罐泄露形成的云团为重气云团,则基于液氨储罐泄露初始时刻重气云团内部浓度、预设伤害扩散浓度和液氨储罐瞬时泄漏形成的重气云团的初始体积,获取对应的待确定轻伤半径、待确定重伤半径和待确定死亡半径,其中,所述预设伤害扩散浓度包括预设轻伤扩散浓度、预设重伤扩散浓度和预设死亡扩散浓度;
通过ε准则,获取所述重气云团让位于环境湍流的第一扩散距离;
判断所述第一扩散距离是否大于所述待确定轻伤半径、并且判断所述第一扩散距离是否大于所述待确定重伤半径,并且判断所述第一扩散距离是否大于所述待确定死亡半径,若所述第一扩散距离大于所述待确定轻伤半径,并且所述第一扩散距离大于所述待确定重伤半径,并且所述第一扩散距离大于所述待确定死亡半径,则将所述待确定轻伤半径确定为轻伤半径,将所述待确定重伤半径确定为重伤半径,将所述待确定死亡半径确定为死亡半径;
若所述第一扩散距离小于等于所述待确定轻伤半径,或者所述第一扩散距离小于等于所述待确定重伤半径,或者所述第一扩散距离小于等于所述待确定死亡半径,则通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径;
若所述液氨储罐泄露形成的云团为非重气云团,则通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径。
基于上述实施例,图3为本发明实施例中的液氨储罐泄露伤害半径获取设备的框架示意图。请参考图3,本发明实施例提供一种液氨储罐泄露伤害半径获取设备,包括:处理器(processor)310、通信接口(Communications Interface)320、存储器(memory)330和总线340,其中,处理器310,通信接口320,存储器330通过总线340完成相互间的通信。处理器310可以调用存储器330中的逻辑指令,以执行如下方法,包括:基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
本发明实施例公开一种计算机程序产品,所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,所述计算机程序包括程序指令,当所述程序指令被计算机执行时,计算机能够执行上述各方法实施例所提供的获取方法,例如包括:基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
基于上述实施例,本发明实施例提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行上述各方法实施例所提供的获取方法,例如包括:基于所述液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
本领域普通技术人员可以理解:实现上述设备实施例或方法实施例仅仅是示意性的,其中所述处理器和所述存储器可以是物理上分离的部件也可以不是物理上分离的,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如U盘、移动硬盘、ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
本发明实施例提供的液氨储罐泄露伤害半径获取方法和设备,通过设置判断液氨储罐泄露形成的云团是否为重气云团,并基于判断结果,获取对应的伤害半径,能够将重气云团扩散和非重气云团的扩散有机结合起来,使得危险物质的扩散过程成为一个连续的变化过程,进而准确获取液氨储罐泄漏浓度扩散的伤害半径。本发明实施例提供的液氨储罐泄露伤害半径获取方法不仅可以模拟出稳定状态下的液氨持续泄漏扩散态势,还可以动态的模拟出液氨持续泄漏扩散随时间变化的过程。本发明实施例对典型的重气或高斯扩散模型进行改进,将重气云团扩散和非重气云团的高斯扩散有机结合起来,使得危险物质(氨)的扩散过程成为一个连续的变化过程,进而准确确定液氨储罐泄漏浓度扩散的扩散影响半径,也就可根据扩散影响半径指导采取相应的控制措施。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种液氨储罐泄露伤害半径获取方法,其特征在于,包括:
基于液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团;
基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径。
2.根据权利要求1所述的获取方法,其特征在于,所述基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径,具体包括:
若所述液氨储罐泄露形成的云团为重气云团,则基于液氨储罐泄露初始时刻重气云团内部浓度、预设伤害扩散浓度和液氨储罐瞬时泄漏形成的重气云团的初始体积,获取对应的待确定轻伤半径、待确定重伤半径和待确定死亡半径,其中,所述预设伤害扩散浓度包括预设轻伤扩散浓度、预设重伤扩散浓度和预设死亡扩散浓度;
通过ε准则,获取所述重气云团让位于环境湍流的第一扩散距离;
若所述第一扩散距离大于所述待确定轻伤半径,并且所述第一扩散距离大于所述待确定重伤半径,并且所述第一扩散距离大于所述待确定死亡半径,则将所述待确定轻伤半径确定为轻伤半径,将所述待确定重伤半径确定为重伤半径,将所述待确定死亡半径确定为死亡半径。
3.根据权利要求2所述的获取方法,其特征在于,若所述第一扩散距离小于等于所述待确定轻伤半径,或者所述第一扩散距离小于等于所述待确定重伤半径,或者所述第一扩散距离小于等于所述待确定死亡半径,则通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径。
4.根据权利要求1所述的获取方法,其特征在于,所述基于所述液氨储罐泄露形成的云团是否为重气云团的判断结果,获取对应的伤害半径,具体包括:
若所述液氨储罐泄露形成的云团为非重气云团,则通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径。
5.根据权利要求1所述的获取方法,其特征在于,所述基于液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团,之前还包括:
获取所述液氨储罐泄露的闪蒸量。
6.根据权利要求1所述的获取方法,其特征在于,所述基于液氨储罐泄露的闪蒸量,判断所述液氨储罐泄露形成的云团是否为重气云团,具体包括:
基于所述液氨储罐泄露的闪蒸量,获取所述液氨储罐泄漏形成的云团的初始体积;
基于所述液氨储罐泄漏形成的云团的初始体积、环境风速、重力加速度、云团的初始密度和环境空气密度,判断所述液氨储罐泄露形成的云团是否为重气云团。
7.根据权利要求3所述的获取方法,其特征在于,所述通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径,具体包括:
通过高斯模型,获取所述重气云团扩散为非重气云团并基于环境湍流引发的第二扩散距离;
根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度;
获取和所述预设轻伤扩散浓度的差值的绝对值最小的泄露最大浓度的第一位置,并且,获取和所述预设重伤扩散浓度的差值的绝对值最小的泄露最大浓度的第二位置,并且,获取和所述预设死亡扩散浓度的差值的绝对值最小的泄露最大浓度的第三位置;
将所述第一位置和所述泄露原点之间的距离作为轻伤半径,并且,将所述第二位置和所述泄露原点之间的距离作为重伤半径,并且,将所述第三位置和所述泄露原点之间的距离作为死亡半径。
8.根据权利要求4所述的获取方法,其特征在于,所述通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离,并根据所述第二扩散距离获取对应的轻伤半径、重伤半径和死亡半径,具体包括:
通过高斯模型,获取所述非重气云团基于环境湍流引发的第二扩散距离;
根据所述第二扩散距离,获取所述液氨储罐泄露的泄露原点周围任一位置的泄露最大浓度,其中,所述泄露最大浓度为瞬时泄漏最大浓度或者连续泄漏最大浓度;
获取和预设轻伤扩散浓度的差值的绝对值最小的泄露最大浓度的第一位置,并且,获取和预设重伤扩散浓度的差值的绝对值最小的泄露最大浓度的第二位置,并且,获取和预设死亡扩散浓度的差值的绝对值最小的泄露最大浓度的第三位置;
将所述第一位置和所述泄露原点之间的距离作为轻伤半径,并且,将所述第二位置和所述泄露原点之间的距离作为重伤半径,并且,将所述第三位置和所述泄露原点之间的距离作为死亡半径。
9.根据权利要求5所述的获取方法,其特征在于,所述获取所述液氨储罐泄露的闪蒸量,具体包括:
通过伯努利方程,获取所述液氨储罐的泄漏量;
基于所述泄漏量,通过闪蒸计算,获取所述液氨储罐泄露的闪蒸量。
10.一种液氨储罐泄露伤害半径获取设备,其特征在于,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至9任一所述的获取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810917779.1A CN109086539A (zh) | 2018-08-13 | 2018-08-13 | 液氨储罐泄露伤害半径获取方法和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810917779.1A CN109086539A (zh) | 2018-08-13 | 2018-08-13 | 液氨储罐泄露伤害半径获取方法和设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109086539A true CN109086539A (zh) | 2018-12-25 |
Family
ID=64834350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810917779.1A Pending CN109086539A (zh) | 2018-08-13 | 2018-08-13 | 液氨储罐泄露伤害半径获取方法和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109086539A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111125919A (zh) * | 2019-12-26 | 2020-05-08 | 清华大学 | 疏散区域计算方法、装置、电子设备和存储介质 |
CN111314430A (zh) * | 2020-01-19 | 2020-06-19 | 北京戴纳实验科技有限公司 | 一种气罐信息采集系统 |
CN113128755A (zh) * | 2021-03-31 | 2021-07-16 | 东南大学 | 液氢泄露扩散范围的预测方法及系统 |
CN114441513A (zh) * | 2020-11-03 | 2022-05-06 | 中国石油化工股份有限公司 | 一种含硫化氢设备泄漏扩散距离的估算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011192883A (ja) * | 2010-03-16 | 2011-09-29 | Toshiba Corp | プロセスシミュレーションをコンピュータに実行させるプログラム |
CN205244839U (zh) * | 2015-12-22 | 2016-05-18 | 湖北瑞辰祥生态肥业有限公司 | 一种液氨泄露报警及处理系统 |
CN105868426A (zh) * | 2015-01-20 | 2016-08-17 | 中国科学院上海高等研究院 | 一种基于烟团模型的重气持续泄露模拟方法 |
-
2018
- 2018-08-13 CN CN201810917779.1A patent/CN109086539A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011192883A (ja) * | 2010-03-16 | 2011-09-29 | Toshiba Corp | プロセスシミュレーションをコンピュータに実行させるプログラム |
CN105868426A (zh) * | 2015-01-20 | 2016-08-17 | 中国科学院上海高等研究院 | 一种基于烟团模型的重气持续泄露模拟方法 |
CN205244839U (zh) * | 2015-12-22 | 2016-05-18 | 湖北瑞辰祥生态肥业有限公司 | 一种液氨泄露报警及处理系统 |
Non-Patent Citations (3)
Title |
---|
宇德明: "《易燃、易爆、有毒危险品储运过程定量风险评价》", 30 June 2000 * |
彭林: "液氨储罐事故泄漏环境风险评价探讨", 《广东化工》 * |
黄德寅等: "化工厂制冷装置液氨泄漏急性职业中毒定量风险评估", 《OCCUP HEALTH & EMERG RESCUE》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111125919A (zh) * | 2019-12-26 | 2020-05-08 | 清华大学 | 疏散区域计算方法、装置、电子设备和存储介质 |
CN111314430A (zh) * | 2020-01-19 | 2020-06-19 | 北京戴纳实验科技有限公司 | 一种气罐信息采集系统 |
CN111314430B (zh) * | 2020-01-19 | 2022-08-05 | 北京戴纳实验科技有限公司 | 一种气罐信息采集系统 |
CN114441513A (zh) * | 2020-11-03 | 2022-05-06 | 中国石油化工股份有限公司 | 一种含硫化氢设备泄漏扩散距离的估算方法 |
CN114441513B (zh) * | 2020-11-03 | 2024-02-09 | 中国石油化工股份有限公司 | 一种含硫化氢设备泄漏扩散距离的估算方法 |
CN113128755A (zh) * | 2021-03-31 | 2021-07-16 | 东南大学 | 液氢泄露扩散范围的预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109086539A (zh) | 液氨储罐泄露伤害半径获取方法和设备 | |
CN108932394B (zh) | 有毒气体泄漏危害范围的确定方法 | |
CN108122051B (zh) | 一种基于无人机探测的危险介质泄漏过程实时动态预测方法 | |
CN103605822B (zh) | 石油化工码头液态化学品管道泄漏危险性动态分级方法 | |
CN108171373A (zh) | 一种化工园区毒气泄露逃生路径规划方法 | |
CN104182588A (zh) | 一种重气持续泄漏扩散的动态模拟方法 | |
CN108108525A (zh) | 基于地理信息系统的气体泄漏事故模拟推演方法及装置 | |
CN112989717B (zh) | 一种加氢站火灾事故热辐射评估方法及系统 | |
CN102096768A (zh) | 危险化学品泄漏污染评估方法 | |
CN106503395A (zh) | 一种lng储罐泄露的定量模拟分析实现方法 | |
CN105868426B (zh) | 一种基于烟团模型的重气持续泄露模拟方法 | |
CN107871025B (zh) | 基于改进人工蜂群算法的气体传感器优化部署方法及系统 | |
Khan et al. | TORAP—a new tool for conducting rapid risk-assessments in petroleum refineries and petrochemical industries | |
Niazi et al. | Effect of wind directions on human injury and fatality risk modeling due to vapor cloud explosion in offshore platforms | |
Yan et al. | Analysis on typical parameters of bodies of revolution affecting the sonic boom | |
CN106523908A (zh) | 一种lng储罐泄漏浓度扩散分析实现方法 | |
Xin et al. | Dynamic simulation and quantitative risk assessment of indoor heavy gas diffusion | |
CN116341890A (zh) | 一种海洋油气平台安全风险量化评估分级技术方法 | |
CN106769496A (zh) | 一种lng储罐泄漏爆炸超压强度分析实现方法 | |
CN111624303B (zh) | 一种有害气体泄漏检测方法、系统及计算机存储介质 | |
Kang | An examination on the dispersion characteristics of boil-off gas in vent mast exit of membrane type LNG carriers | |
Yuzhong et al. | LNG accident consequence QRA considering domino effect | |
CN112016738B (zh) | 一种吸入性有毒物质的危害预测方法 | |
Abbassi et al. | Application of integrated risk-based approach to safety design of high speed LNG powered vessels | |
Godbole et al. | Atmospheric Dispersion of CO2 Following the Full-Scale Burst Test of a High-Pressure CO2 Pipeline |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181225 |