一种气体泄漏扩散监测方法及系统与流程

1.本发明涉及气体泄漏监测技术领域，特别涉及一种气体泄漏扩散监测方法及系统。


背景技术：

2.有毒有害气体广泛存在于化学生产过程中和人们的日常生活中，一旦使用或维护不当则极易造成泄漏，往往引发火灾、爆炸、灼伤、中毒、窒息等重大事故。
3.气体泄漏事故发生后，需尽快确认气体泄漏的影响参数并监测气体泄漏的形态变化、方向走向等，以便采取正确果断的应急措施，将危险控制在可控范围，将损失降低到最低。
4.目前，在气体扩散监测方面，主要采用的模型是高斯模型，但其无法准确估算气体泄漏的初始源强，从而导致最终获得的结果不准确。


技术实现要素：

5.本发明实施例的目的在于提供一种气体泄漏扩散监测方法及系统，以解决上述问题。
6.为实现上述目的，本发明提供如下技术方案：
7.一方面，提供一种气体泄漏扩散监测方法，其包括如下步骤：
8.确定泄漏气体的流动状态；
9.根据泄漏气体的流动状态获取泄漏气体源强；
10.根据泄漏气体的气团类型选择气团浓度计算参数；
11.根据获取的气团浓度计算参数计算泄漏气体气团浓度；
12.根据泄漏气体的气团浓度获取气体泄漏影响参数，并根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型，并将泄漏气体的气团的三维仿真模型在gis地图上进行加载。
13.优选的，根据公式(1)判断泄漏气体的流动状态；
[0014][0015]
其中，p为气体输送管道内的压力(pa)；p0为气体输送管道外的压力(pa)；k为比热比；
[0016]
当公式(1)成立时，确定泄漏气体处于音速流动状态，若不成立，确定泄漏气体处于亚音速流动状态。
[0017]
优选的，当泄漏气体处于音速流动状态时，根据公式(2)获取泄漏气体源强q：
[0018][0019]
其中，cd为气体泄漏系数；m为泄漏气体的摩尔质量；t为泄漏气体的初始热力学温度；a为气体泄漏点的面积；r为热力学常数。
[0020]
优选的，当泄漏气体处于亚音速流动状态时，根据公式(3)获取泄漏气体源强q：
[0021][0022]
其中，cd为气体泄漏系数；m为泄漏气体的摩尔质量；t为泄漏气体的初始热力学温度；a为气体泄漏点的面积；r为热力学常数。
[0023]
优选的，泄漏气体的气团类型确定方法如下：
[0024]
若泄漏气体的相对密度小于或等于1，且为连续泄漏时所形成的气团为羽状气团；
[0025]
若泄漏气体的相对密度比空气密度大，则确定该泄漏气体为重气，当重气泄漏后形成的重气气团与空气之间的密度差小于0.001kg/m3时，则认为重气气团转变成中性状态，且进入被动气体扩散阶段，进入被动气体扩散阶段后所形成的气团为重气体作用后气团。
[0026]
优选的，如泄漏气体的气团类型为羽状气团，则气团浓度计算参数包括：
[0027]
泄漏气体离开管道的气体泄漏点后，沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z；
[0028]
泄漏气体源强q；
[0029]
泄漏气体在顺风向的扩散系数σ
x
；
[0030]
泄漏气体在侧风向的扩散系数σy；
[0031]
泄漏气体在垂直顺风向方向的扩散系数σz；
[0032]
排口高度h。
[0033]
优选的，如为羽状气团，则按公式(4)计算泄漏气体气团浓度c(x，y，z)：
[0034][0035]
优选的，如泄漏气体的气团类型为为重气体作用后气团，则气团浓度计算参数包括：
[0036]
任一气体泄漏点的泄漏气体浓度c；
[0037]
泄漏气体离开管道的气体泄漏点后，沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z；
[0038]
泄漏气体源强q；
[0039]
泄漏气体在顺风向的扩散系数σ
x
；
[0040]
泄漏气体在侧风向的扩散系数σy；
[0041]
泄漏气体在垂直顺风向方向的扩散系数σz；
[0042]
排口高度h；
[0043]
风速u；
[0044]
时间t。
[0045]
优选的，如为重气体作用后气团，则按公式(5)计算泄漏气体气团浓度c(x，y，z，t)：
[0046][0047]
还提供一种气体泄漏扩散监测系统，其包括：
[0048]
泄漏气体流动状态判断单元，其用于在确定产生泄漏点后，判断泄漏气体的流动状态；
[0049]
音速泄漏气体源强计算单元，其用于获取处于音速流动状态的泄漏气体的泄漏气体源强；
[0050]
亚音速泄漏气体源强计算单元，其用于获取处于音速流动状态的泄漏气体的泄漏气体源强；
[0051]
泄漏气体气团浓度计算单元，其用于计算气团类型为羽状气团的泄漏气体的气团浓度，以及，计算气团类型为重气体作用后气团的泄漏气体的气团浓度；
[0052]
气体泄漏影响参数获取单元，其用于根据泄漏气体的气团浓度获取气体泄漏影响参数；
[0053]
三维建模单元，其用于根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型；
[0054]
加载单元，其用于在gis地图上加载所述泄漏气体的气团的三维仿真模型。
[0055]
本发明的技术效果：本发明基于源强估算技术对泄漏气体的初始形态进行准确预判，确保后续将其输入扩散监测模型后，可以对气体扩散进行准确的形态预测，以保证扩散监测结果的准确性，同时以三维仿真 gis的方式展示监测结果，由此可直观获取气体泄漏影响范围，以便及时做出正确的应急处置。
附图说明
[0056]
图1是本发明气体泄漏扩散监测方法的步骤流程图；
[0057]
图2是本发明中xyz距离的示意图；
[0058]
图3是本发明中顺风向、侧风向以及垂直顺风向方向的示意图；
[0059]
图4是本发明中在gis地图上加载气团三维仿真模型的效果图；
[0060]
图5是本发明气体泄漏扩散监测系统的结构示意图。
具体实施方式
[0061]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。
[0062]
实施例1：
[0063]
请参照图1，本实施例提供了一种气体泄漏扩散监测方法，其具体包括如下步骤：
[0064]
s1、采用人工巡视或者监测两阀室之间的管道压强变化(压强变小即说明存在泄漏点)等方式确定管道是否产生气体泄漏点；；
[0065]
s2、由于泄漏气体在音速与亚音速状态下流动时的泄漏气体源强计算公式完全不同，因此需要先确定泄漏气体的流动状态，其中，所述泄漏气体包括天然气、h2、n2、o2、co2、so2、nh3、co等，具体的，需根据公式(1)判断泄漏气体的流动状态；
[0066][0067]
其中，p为气体输送管道内的压力(pa)；p0为气体输送管道外的压力(pa)；k为比热比，即定压比热与定容比热的比值，其取值范围参见表1(本实施例中，k的取值为1.306)：
[0068]
表1比热比取值范围表
[0069]
序号热力学温度(t0(k))比热比(k)101.40122001.39034001.36846001.34958001.332610001.320715001.306
[0070]
当公式(1)成立时，确定泄漏气体处于音速流动状态，若不成立，确定泄漏气体处于亚音速流动状态；
[0071]
s3、当泄漏气体处于音速流动状态时，根据公式(2)获取泄漏气体源强q：
[0072][0073]
当泄漏气体处于亚音速流动状态时，根据公式(3)获取泄漏气体源强q：
[0074][0075]
其中，cd为气体泄漏系数，其取值范围与气体泄漏点的形状有关，两者之间的关系具体参见表2；
[0076]
表2气体泄漏系数取值范围表
[0077]
[0078][0079]
m为泄漏气体的摩尔质量；
[0080]
t为泄漏气体的初始热力学温度(k)；
[0081]
a为气体泄漏点的面积(m2)；
[0082]
r为热力学常数，即8.314j/(mol
·
k)。
[0083]
s4、根据泄漏气体的气团类型选择气团浓度计算参数；
[0084]
具体的，由于泄漏气体的气团类型与泄漏气体的类别以及密度有关，因此，泄漏气体的气团类型确定方法如下：
[0085]
若泄漏气体的相对密度小于或等于1，且为连续泄漏时所形成的气团为羽状气团；
[0086]
若泄漏气体的相对密度比空气密度大，则确定该泄漏气体为重气，当重气泄漏后经过重气扩散阶段形成的重气气团与空气之间的密度差小于0.001kg/m3时，则认为重气气团转变成中性状态，且进入被动气体扩散阶段，进入被动气体扩散阶段后所形成的气团为重气体作用后气团；
[0087]
进一步的，如泄漏气体的气团类型为羽状气团，则气团浓度计算参数包括：
[0088]
泄漏气体离开管道的气体泄漏点100后，沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z；且距离x、y、z单位均为米；
[0089]
泄漏气体源强q(kg/s)，需要说明的是，需根据泄漏气体的流动状态来选择对应的泄漏气体源强q；
[0090]
泄漏气体在顺风向(也即管道轴向)的扩散系数σ
x
；
[0091]
泄漏气体在侧风向(也即管道左右方向)的扩散系数σy；
[0092]
泄漏气体在垂直顺风向方向(也即垂直于管道上下方向)的扩散系数σz，顺风向、侧风向、垂直顺风向方向的示意图如图3所示；
[0093]
需要说明的是，扩散系数仅与泄漏气体类型有关，同一种泄漏气体的扩散系数σ
x
、σy、σz相同，且扩散系数σ
x
、σy、σz的取值范围如表3所示：
[0094]
表3泄漏气体的扩散系数表
[0095]
序号气体类型扩散系统1h20.5112n20.1323o20.1784co20.1385so20.1036nh30.2007h2o0.2208hcl0.130
[0096]
排口高度h(m)，排口高度即泄漏气体在空中到达的最高点与泄漏点之间的距离；
[0097]
如泄漏气体气团类型为重气体作用后气团，则气团浓度计算参数包括：
[0098]
任一气体泄漏点100的泄漏气体浓度c(kg/m3)；
[0099]
泄漏气体离开管道的气体泄漏点100后，沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z；且距离x、y、z单位均为米；
[0100]
泄漏气体源强q(kg/s)；
[0101]
泄漏气体在顺风向(也即管道轴向)的扩散系数σ
x
；
[0102]
泄漏气体在侧风向(也即管道左右方向)的扩散系数σy；
[0103]
泄漏气体在垂直顺风向方向(也即垂直于管道上下方向)的扩散系数σz，顺风向、侧风向、垂直顺风向方向的示意图如图3所示；
[0104]
排口高度h(m)；
[0105]
风速u(m/s)；
[0106]
时间t(s)；
[0107]
s5、根据获取的气团浓度计算参数计算泄漏气体气团浓度，具体包括：
[0108]
如为羽状气团，则按公式(4)计算泄漏气体气团浓度c(x，y，z)(kg/m3)：
[0109][0110]
如为重气体作用后气团，则按公式(5)计算泄漏气体气团浓度c(x，y，z，t)(kg/m3)：
[0111][0112]
s6、根据泄漏气体的气团浓度获取气体泄漏影响参数，且所述气体泄漏影响参数包括：泄漏气体体积、泄漏气体气带宽度等其中的一种或几种；
[0113]
将气体泄漏影响参数载入cesium软件的三维粒子组件模块，以时间为主轴，自动建模生成泄漏气体的气团的三维仿真模型；
[0114]
根据不同的业务场景，在gis地图上加载所述泄漏气体的气团的三维仿真模型(如图4所示)，以展现气体扩散监测结果。
[0115]
实施例2：
[0116]
本实施例提供了一种用于实现实施例1所述气体泄漏扩散监测方法的气体泄漏扩散监测系统，如图5所示，其包括：
[0117]
泄漏气体流动状态判断单元1，其用于在确定产生泄漏点后，根据公式(1)判断泄漏气体的流动状态；
[0118]
音速泄漏气体源强计算单元2，其用于根据公式(2)获取处于音速流动状态的泄漏气体的泄漏气体源强；
[0119]
亚音速泄漏气体源强计算单元3，其用于根据公式(3)获取处于音速流动状态的泄
漏气体的泄漏气体源强；
[0120]
泄漏气体气团浓度计算单元4，其用于根据公式(4)计算气团类型为羽状气团的泄漏气体的气团浓度，以及，根据公式(5)计算气团类型为重气体作用后气团的泄漏气体的气团浓度；
[0121]
气体泄漏影响参数获取单元5，其用于根据泄漏气体的气团浓度获取气体泄漏影响参数；
[0122]
三维建模单元6，其用于根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型；
[0123]
加载单元7，其用于在gis地图上加载所述泄漏气体的气团的三维仿真模型。
[0124]
综上所述，本发明基于源强估算技术对泄漏气体的初始形态进行准确预判，确保后续将其输入扩散监测模型后，可以对气体扩散进行准确的形态预测，以保证扩散监测结果的准确性，同时以三维仿真 gis的方式展示监测结果，由此可直观获取气体泄漏影响范围，以便及时做出正确的应急处置
[0125]
在不冲突的情况下，本文中上述实施例及实施例中的特征可以相互结合。
[0126]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。