CN115691705A - 一种气体泄漏扩散监测方法及系统 - Google Patents
一种气体泄漏扩散监测方法及系统 Download PDFInfo
- Publication number
- CN115691705A CN115691705A CN202211410390.0A CN202211410390A CN115691705A CN 115691705 A CN115691705 A CN 115691705A CN 202211410390 A CN202211410390 A CN 202211410390A CN 115691705 A CN115691705 A CN 115691705A
- Authority
- CN
- China
- Prior art keywords
- gas
- mass
- leakage
- leaked
- leaking
- 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.)
- Granted
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Examining Or Testing Airtightness (AREA)
Abstract
本发明提供一种气体泄漏扩散监测方法及系统,其包括如下步骤:确定泄漏气体的流动状态;根据泄漏气体的流动状态获取泄漏气体源强;根据泄漏气体的气团类型选择气团浓度计算参数;根据获取的气团浓度计算参数计算泄漏气体气团浓度;根据泄漏气体的气团浓度获取气体泄漏影响参数,并根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型,并将泄漏气体的气团的三维仿真模型在GIS地图上进行加载。本发明基于源强估算技术对泄漏气体的初始形态进行准确预判,确保后续将其输入扩散监测模型后,可以对气体扩散进行准确的形态预测,以保证扩散监测结果的准确性。
Description
技术领域
本发明涉及气体泄漏监测技术领域,特别涉及一种气体泄漏扩散监测方法及系统。
背景技术
有毒有害气体广泛存在于化学生产过程中和人们的日常生活中,一旦使用或维护不当则极易造成泄漏,往往引发火灾、爆炸、灼伤、中毒、窒息等重大事故。
气体泄漏事故发生后,需尽快确认气体泄漏的影响参数并监测气体泄漏的形态变化、方向走向等,以便采取正确果断的应急措施,将危险控制在可控范围,将损失降低到最低。
目前,在气体扩散监测方面,主要采用的模型是高斯模型,但其无法准确估算气体泄漏的初始源强,从而导致最终获得的结果不准确。
发明内容
本发明实施例的目的在于提供一种气体泄漏扩散监测方法及系统,以解决上述问题。
为实现上述目的,本发明提供如下技术方案:
一方面,提供一种气体泄漏扩散监测方法,其包括如下步骤:
确定泄漏气体的流动状态;
根据泄漏气体的流动状态获取泄漏气体源强;
根据泄漏气体的气团类型选择气团浓度计算参数;
根据获取的气团浓度计算参数计算泄漏气体气团浓度;
根据泄漏气体的气团浓度获取气体泄漏影响参数,并根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型,并将泄漏气体的气团的三维仿真模型在GIS地图上进行加载。
优选的,根据公式(1)判断泄漏气体的流动状态;
其中,p为气体输送管道内的压力(Pa);p0为气体输送管道外的压力(Pa);k为比热比;
当公式(1)成立时,确定泄漏气体处于音速流动状态,若不成立,确定泄漏气体处于亚音速流动状态。
优选的,当泄漏气体处于音速流动状态时,根据公式(2)获取泄漏气体源强Q:
其中,Cd为气体泄漏系数;M为泄漏气体的摩尔质量;T为泄漏气体的初始热力学温度;A为气体泄漏点的面积;R为热力学常数。
优选的,当泄漏气体处于亚音速流动状态时,根据公式(3)获取泄漏气体源强Q:
其中,Cd为气体泄漏系数;M为泄漏气体的摩尔质量;T为泄漏气体的初始热力学温度;A为气体泄漏点的面积;R为热力学常数。
优选的,泄漏气体的气团类型确定方法如下:
若泄漏气体的相对密度小于或等于1,且为连续泄漏时所形成的气团为羽状气团;
若泄漏气体的相对密度比空气密度大,则确定该泄漏气体为重气,当重气泄漏后形成的重气气团与空气之间的密度差小于0.001kg/m3时,则认为重气气团转变成中性状态,且进入被动气体扩散阶段,进入被动气体扩散阶段后所形成的气团为重气体作用后气团。
优选的,如泄漏气体的气团类型为羽状气团,则气团浓度计算参数包括:
泄漏气体离开管道的气体泄漏点后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;
泄漏气体源强Q;
泄漏气体在顺风向的扩散系数σx;
泄漏气体在侧风向的扩散系数σy;
泄漏气体在垂直顺风向方向的扩散系数σz;
排口高度H。
优选的,如为羽状气团,则按公式(4)计算泄漏气体气团浓度C(x,y,z):
优选的,如泄漏气体的气团类型为为重气体作用后气团,则气团浓度计算参数包括:
任一气体泄漏点的泄漏气体浓度C;
泄漏气体离开管道的气体泄漏点后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;
泄漏气体源强Q;
泄漏气体在顺风向的扩散系数σx;
泄漏气体在侧风向的扩散系数σy;
泄漏气体在垂直顺风向方向的扩散系数σz;
排口高度H;
风速u;
时间t。
优选的,如为重气体作用后气团,则按公式(5)计算泄漏气体气团浓度C(x,y,z,t):
还提供一种气体泄漏扩散监测系统,其包括:
泄漏气体流动状态判断单元,其用于在确定产生泄漏点后,判断泄漏气体的流动状态;
音速泄漏气体源强计算单元,其用于获取处于音速流动状态的泄漏气体的泄漏气体源强;
亚音速泄漏气体源强计算单元,其用于获取处于音速流动状态的泄漏气体的泄漏气体源强;
泄漏气体气团浓度计算单元,其用于计算气团类型为羽状气团的泄漏气体的气团浓度,以及,计算气团类型为重气体作用后气团的泄漏气体的气团浓度;
气体泄漏影响参数获取单元,其用于根据泄漏气体的气团浓度获取气体泄漏影响参数;
三维建模单元,其用于根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型;
加载单元,其用于在GIS地图上加载所述泄漏气体的气团的三维仿真模型。
本发明的技术效果:本发明基于源强估算技术对泄漏气体的初始形态进行准确预判,确保后续将其输入扩散监测模型后,可以对气体扩散进行准确的形态预测,以保证扩散监测结果的准确性,同时以三维仿真+GIS的方式展示监测结果,由此可直观获取气体泄漏影响范围,以便及时做出正确的应急处置。
附图说明
图1是本发明气体泄漏扩散监测方法的步骤流程图;
图2是本发明中xyz距离的示意图;
图3是本发明中顺风向、侧风向以及垂直顺风向方向的示意图;
图4是本发明中在GIS地图上加载气团三维仿真模型的效果图;
图5是本发明气体泄漏扩散监测系统的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
实施例1:
请参照图1,本实施例提供了一种气体泄漏扩散监测方法,其具体包括如下步骤:
S1、采用人工巡视或者监测两阀室之间的管道压强变化(压强变小即说明存在泄漏点)等方式确定管道是否产生气体泄漏点;;
S2、由于泄漏气体在音速与亚音速状态下流动时的泄漏气体源强计算公式完全不同,因此需要先确定泄漏气体的流动状态,其中,所述泄漏气体包括天然气、H2、N2、O2、CO2、SO2、NH3、CO等,具体的,需根据公式(1)判断泄漏气体的流动状态;
其中,p为气体输送管道内的压力(Pa);p0为气体输送管道外的压力(Pa);k为比热比,即定压比热与定容比热的比值,其取值范围参见表1(本实施例中,k的取值为1.306):
表1比热比取值范围表
序号 | 热力学温度(T<sub>0</sub>(K)) | 比热比(k) |
1 | 0 | 1.401 |
2 | 200 | 1.390 |
3 | 400 | 1.368 |
4 | 600 | 1.349 |
5 | 800 | 1.332 |
6 | 1000 | 1.320 |
7 | 1500 | 1.306 |
当公式(1)成立时,确定泄漏气体处于音速流动状态,若不成立,确定泄漏气体处于亚音速流动状态;
S3、当泄漏气体处于音速流动状态时,根据公式(2)获取泄漏气体源强Q:
当泄漏气体处于亚音速流动状态时,根据公式(3)获取泄漏气体源强Q:
其中,Cd为气体泄漏系数,其取值范围与气体泄漏点的形状有关,两者之间的关系具体参见表2;
表2气体泄漏系数取值范围表
M为泄漏气体的摩尔质量;
T为泄漏气体的初始热力学温度(K);
A为气体泄漏点的面积(m2);
R为热力学常数,即8.314J/(mol·K)。
S4、根据泄漏气体的气团类型选择气团浓度计算参数;
具体的,由于泄漏气体的气团类型与泄漏气体的类别以及密度有关,因此,泄漏气体的气团类型确定方法如下:
若泄漏气体的相对密度小于或等于1,且为连续泄漏时所形成的气团为羽状气团;
若泄漏气体的相对密度比空气密度大,则确定该泄漏气体为重气,当重气泄漏后经过重气扩散阶段形成的重气气团与空气之间的密度差小于0.001kg/m3时,则认为重气气团转变成中性状态,且进入被动气体扩散阶段,进入被动气体扩散阶段后所形成的气团为重气体作用后气团;
进一步的,如泄漏气体的气团类型为羽状气团,则气团浓度计算参数包括:
泄漏气体离开管道的气体泄漏点100后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;且距离x、y、z单位均为米;
泄漏气体源强Q(kg/s),需要说明的是,需根据泄漏气体的流动状态来选择对应的泄漏气体源强Q;
泄漏气体在顺风向(也即管道轴向)的扩散系数σx;
泄漏气体在侧风向(也即管道左右方向)的扩散系数σy;
泄漏气体在垂直顺风向方向(也即垂直于管道上下方向)的扩散系数σz,顺风向、侧风向、垂直顺风向方向的示意图如图3所示;
需要说明的是,扩散系数仅与泄漏气体类型有关,同一种泄漏气体的扩散系数σx、σy、σz相同,且扩散系数σx、σy、σz的取值范围如表3所示:
表3泄漏气体的扩散系数表
序号 | 气体类型 | 扩散系统 |
1 | H<sub>2</sub> | 0.511 |
2 | N<sub>2</sub> | 0.132 |
3 | O<sub>2</sub> | 0.178 |
4 | CO<sub>2</sub> | 0.138 |
5 | SO<sub>2</sub> | 0.103 |
6 | NH<sub>3</sub> | 0.200 |
7 | H<sub>2</sub>O | 0.220 |
8 | HCL | 0.130 |
排口高度H(m),排口高度即泄漏气体在空中到达的最高点与泄漏点之间的距离;
如泄漏气体气团类型为重气体作用后气团,则气团浓度计算参数包括:
任一气体泄漏点100的泄漏气体浓度C(kg/m3);
泄漏气体离开管道的气体泄漏点100后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;且距离x、y、z单位均为米;
泄漏气体源强Q(kg/s);
泄漏气体在顺风向(也即管道轴向)的扩散系数σx;
泄漏气体在侧风向(也即管道左右方向)的扩散系数σy;
泄漏气体在垂直顺风向方向(也即垂直于管道上下方向)的扩散系数σz,顺风向、侧风向、垂直顺风向方向的示意图如图3所示;
排口高度H(m);
风速u(m/s);
时间t(s);
S5、根据获取的气团浓度计算参数计算泄漏气体气团浓度,具体包括:
如为羽状气团,则按公式(4)计算泄漏气体气团浓度C(x,y,z)(kg/m3):
如为重气体作用后气团,则按公式(5)计算泄漏气体气团浓度C(x,y,z,t)(kg/m3):
S6、根据泄漏气体的气团浓度获取气体泄漏影响参数,且所述气体泄漏影响参数包括:泄漏气体体积、泄漏气体气带宽度等其中的一种或几种;
将气体泄漏影响参数载入Cesium软件的三维粒子组件模块,以时间为主轴,自动建模生成泄漏气体的气团的三维仿真模型;
根据不同的业务场景,在GIS地图上加载所述泄漏气体的气团的三维仿真模型(如图4所示),以展现气体扩散监测结果。
实施例2:
本实施例提供了一种用于实现实施例1所述气体泄漏扩散监测方法的气体泄漏扩散监测系统,如图5所示,其包括:
泄漏气体流动状态判断单元1,其用于在确定产生泄漏点后,根据公式(1)判断泄漏气体的流动状态;
音速泄漏气体源强计算单元2,其用于根据公式(2)获取处于音速流动状态的泄漏气体的泄漏气体源强;
亚音速泄漏气体源强计算单元3,其用于根据公式(3)获取处于音速流动状态的泄漏气体的泄漏气体源强;
泄漏气体气团浓度计算单元4,其用于根据公式(4)计算气团类型为羽状气团的泄漏气体的气团浓度,以及,根据公式(5)计算气团类型为重气体作用后气团的泄漏气体的气团浓度;
气体泄漏影响参数获取单元5,其用于根据泄漏气体的气团浓度获取气体泄漏影响参数;
三维建模单元6,其用于根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型;
加载单元7,其用于在GIS地图上加载所述泄漏气体的气团的三维仿真模型。
综上所述,本发明基于源强估算技术对泄漏气体的初始形态进行准确预判,确保后续将其输入扩散监测模型后,可以对气体扩散进行准确的形态预测,以保证扩散监测结果的准确性,同时以三维仿真+GIS的方式展示监测结果,由此可直观获取气体泄漏影响范围,以便及时做出正确的应急处置
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种气体泄漏扩散监测方法,其特征在于,包括如下步骤:
确定泄漏气体的流动状态;
根据泄漏气体的流动状态获取泄漏气体源强;
根据泄漏气体的气团类型选择气团浓度计算参数;
根据获取的气团浓度计算参数计算泄漏气体气团浓度;
根据泄漏气体的气团浓度获取气体泄漏影响参数,并根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型,并将泄漏气体的气团的三维仿真模型在GIS地图上进行加载。
5.如权利要求1所述的气体泄漏扩散监测方法,其特征在于,泄漏气体的气团类型确定方法如下:
若泄漏气体的相对密度小于或等于1,且为连续泄漏时所形成的气团为羽状气团;
若泄漏气体的相对密度比空气密度大,则确定该泄漏气体为重气,当重气泄漏后形成的重气气团与空气之间的密度差小于0.001kg/m3时,则认为重气气团转变成中性状态,且进入被动气体扩散阶段,进入被动气体扩散阶段后所形成的气团为重气体作用后气团。
6.如权利要求5所述的气体泄漏扩散监测方法,其特征在于,如泄漏气体的气团类型为羽状气团,则气团浓度计算参数包括:
泄漏气体离开管道的气体泄漏点后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;
泄漏气体源强Q;
泄漏气体在顺风向的扩散系数σx;
泄漏气体在侧风向的扩散系数σy;
泄漏气体在垂直顺风向方向的扩散系数σz;
排口高度H。
8.如权利要求5所述的气体泄漏扩散监测方法,其特征在于,如泄漏气体的气团类型为为重气体作用后气团,则气团浓度计算参数包括:
任一气体泄漏点的泄漏气体浓度C;
泄漏气体离开管道的气体泄漏点后,沿管道轴向扩散的最大距离x、沿垂直于管道左右方向扩散的最大距离y、沿垂直于管道上下方向扩散的最大距离z;
泄漏气体源强Q;
泄漏气体在顺风向的扩散系数σx;
泄漏气体在侧风向的扩散系数σy;
泄漏气体在垂直顺风向方向的扩散系数σz;
排口高度H;
风速u;
时间t。
10.一种气体泄漏扩散监测系统,其特征在于,包括:
泄漏气体流动状态判断单元,其用于在确定产生泄漏点后,判断泄漏气体的流动状态;
音速泄漏气体源强计算单元,其用于获取处于音速流动状态的泄漏气体的泄漏气体源强;
亚音速泄漏气体源强计算单元,其用于获取处于音速流动状态的泄漏气体的泄漏气体源强;
泄漏气体气团浓度计算单元,其用于计算气团类型为羽状气团的泄漏气体的气团浓度,以及,计算气团类型为重气体作用后气团的泄漏气体的气团浓度;
气体泄漏影响参数获取单元,其用于根据泄漏气体的气团浓度获取气体泄漏影响参数;
三维建模单元,其用于根据气体泄漏影响参数生成泄漏气体的气团的三维仿真模型;
加载单元,其用于在GIS地图上加载所述泄漏气体的气团的三维仿真模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211410390.0A CN115691705B (zh) | 2022-11-11 | 2022-11-11 | 一种气体泄漏扩散监测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211410390.0A CN115691705B (zh) | 2022-11-11 | 2022-11-11 | 一种气体泄漏扩散监测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115691705A true CN115691705A (zh) | 2023-02-03 |
CN115691705B CN115691705B (zh) | 2023-07-25 |
Family
ID=85051470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211410390.0A Active CN115691705B (zh) | 2022-11-11 | 2022-11-11 | 一种气体泄漏扩散监测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115691705B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105701322A (zh) * | 2016-03-22 | 2016-06-22 | 中国特种设备检测研究院 | 一种油气管道事故仿真系统及方法 |
CN106021817A (zh) * | 2016-06-22 | 2016-10-12 | 西南石油大学 | 一种海相气田高含硫天然气集输装置泄漏快速模拟方法 |
CN108038275A (zh) * | 2017-11-28 | 2018-05-15 | 哈尔滨理工大学 | 一种天然气管道泄漏声场的数值仿真及特性分析方法 |
CN110633339A (zh) * | 2018-06-06 | 2019-12-31 | 中国石油化工股份有限公司 | 一种基于gis的石化企业气体泄漏连续动态模拟展示方法 |
CN114239193A (zh) * | 2021-10-20 | 2022-03-25 | 北京市燃气集团有限责任公司 | 燃气管道泄漏仿真模拟方法及装置 |
-
2022
- 2022-11-11 CN CN202211410390.0A patent/CN115691705B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105701322A (zh) * | 2016-03-22 | 2016-06-22 | 中国特种设备检测研究院 | 一种油气管道事故仿真系统及方法 |
CN106021817A (zh) * | 2016-06-22 | 2016-10-12 | 西南石油大学 | 一种海相气田高含硫天然气集输装置泄漏快速模拟方法 |
CN108038275A (zh) * | 2017-11-28 | 2018-05-15 | 哈尔滨理工大学 | 一种天然气管道泄漏声场的数值仿真及特性分析方法 |
CN110633339A (zh) * | 2018-06-06 | 2019-12-31 | 中国石油化工股份有限公司 | 一种基于gis的石化企业气体泄漏连续动态模拟展示方法 |
CN114239193A (zh) * | 2021-10-20 | 2022-03-25 | 北京市燃气集团有限责任公司 | 燃气管道泄漏仿真模拟方法及装置 |
Non-Patent Citations (2)
Title |
---|
李云云: "高斯烟羽模型的改进及在危化品泄漏事故模拟中的应用", pages 2 * |
高峰: "毒性物质泄漏扩散范围的研究", pages 3 * |
Also Published As
Publication number | Publication date |
---|---|
CN115691705B (zh) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110197049B (zh) | 一种基于瞬变反问题的非金属管道泄漏定位方法 | |
CN101714353A (zh) | 声音衰减系统和方法 | |
CN207231786U (zh) | 一种风阀综合性能测试装置 | |
CN108280849A (zh) | 一种综合管廊燃气泄漏浓度场预测矫正及泄漏率估计方法 | |
CN202793710U (zh) | 一种汽车整车的气密性检测装置 | |
CN103034752B (zh) | 空调管道噪声预测系统和方法 | |
CN110348112B (zh) | 一种计算天然气放空管道最大流速的自适应约束方法 | |
CN109210387B (zh) | 一种基于数学模型的气体管道泄漏检测定位的方法 | |
Le Roux et al. | Reduced-scale study of wind influence on mean airflows inside buildings equipped with ventilation systems | |
CN115691705A (zh) | 一种气体泄漏扩散监测方法及系统 | |
Song et al. | Research on a push-pull industrial trough-side exhaust hood based on CFD simulation and experiment | |
CN112485036B (zh) | 一种测试用油雾发生装置 | |
CN113587410A (zh) | 风量测量装置、地铁空调系统和控制方法 | |
CN117073933A (zh) | 一种用于燃气管网的管道压力监测方法及系统 | |
Song et al. | A study on the pollutant control effect of a new push-pull exhaust hood under different pollutant velocities | |
CN116230106A (zh) | 一种有害气体浓度预测方法和应用 | |
CN116341407A (zh) | 一种试车台空气动力学设计方法 | |
CN114239193B (zh) | 燃气管道泄漏仿真模拟方法及装置 | |
CN108984948B (zh) | 一种冷却塔出口雾羽状态评价方法 | |
CN112257209A (zh) | 基于固定监测点数据的室内污染监测系统 | |
Aboelezz | Low speed wind tunnel design and optimization using computational techniques and experimental validation | |
Pečinka et al. | Preliminary design of a low-cost mobile test cell for small gas turbine engines | |
CN112781846A (zh) | 防火门压降及气流再生噪音的测试方法、防火门检测方法 | |
CN207751705U (zh) | 一种负离子发生器标定设备 | |
JP2767088B2 (ja) | 気体状物質の漏洩地点推定方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |