CN114239193B - 燃气管道泄漏仿真模拟方法及装置 - Google Patents

燃气管道泄漏仿真模拟方法及装置 Download PDF

Info

Publication number
CN114239193B
CN114239193B CN202111218747.0A CN202111218747A CN114239193B CN 114239193 B CN114239193 B CN 114239193B CN 202111218747 A CN202111218747 A CN 202111218747A CN 114239193 B CN114239193 B CN 114239193B
Authority
CN
China
Prior art keywords
gas pipeline
gas
leakage
evaluation unit
point
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
Application number
CN202111218747.0A
Other languages
English (en)
Other versions
CN114239193A (zh
Inventor
刘慧�
车明
马旭卿
郭保玲
张玉星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Gas Group Co Ltd
Original Assignee
Beijing Gas Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Gas Group Co Ltd filed Critical Beijing Gas Group Co Ltd
Priority to CN202111218747.0A priority Critical patent/CN114239193B/zh
Publication of CN114239193A publication Critical patent/CN114239193A/zh
Application granted granted Critical
Publication of CN114239193B publication Critical patent/CN114239193B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205Re-meshing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Graphics (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明公开了一种燃气管道泄漏仿真模拟方法及装置,该方法包括:对燃气管道进行评价单元网格的划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点;构建燃气管道泄漏扩散状态模型;设置燃气管道的泄漏参数;基于所述燃气管道的泄漏参数及所述燃气管道泄漏扩散状态模型,计算各观察点的泄漏气体浓度,生成显示信息,以根据所述显示信息可视化各观察点的泄漏气体浓度变化。本发明建立燃气管道三维模型并划分网格,选择合适的观察点,分析燃气管道被破坏后泄漏扩散情况,模拟设定燃气管道参数条件下的燃气管道被第三方破坏后的泄漏过程,分析通风条件对天然气扩散的影响情况。

Description

燃气管道泄漏仿真模拟方法及装置
技术领域
本发明涉及燃气管道泄漏技术领域,特别地涉及一种燃气管道泄漏仿真模拟方法及装置。
背景技术
燃气管道作为重要的城镇基础设施,已经遍布城镇的大街小巷。目前燃气管道的铺设方式主要为埋地型,管道长期与土壤接触,受到其中水分和化学物质的氧化腐蚀,以及泥石流和土壤沉降的因素的影响,造成管壁减薄破裂导致气体泄漏。
当燃气管道受到的外部载荷发生破坏时会导致燃气泄漏,由于城镇燃气管道一般位于人员密集的区域,燃气泄漏很容易引起火灾甚至爆炸,将造成严重的人员伤亡和财产损失,因此国内外有大量针对燃气泄漏及后果分析的研究。
现有的管道燃气泄漏分析方法有:(1)假设了天然气管道的泄漏几何区域,基于有限容积原理和多孔介质理论提出了天然气的泄漏控制方程,并且使用软件模拟了埋地天然气管道的泄漏过程;(2)分析了埋地天然气管道的泄漏规律以及泄漏气体在土壤中扩散的主要影响因素,研究了泄漏导致的爆炸和中毒事故的特点;(3)使用液体泄漏及扩散实验平台研究了气体的泄漏和扩散过程,得到了天然气管道泄漏时的变化规律,并且根据实验数据修正了理论计算公式。
上述的管道燃气泄漏分析方法主要侧重于燃气泄漏和扩散过程的模拟,缺乏结合实际工程的泄漏后果的应用分析,没有给出较全面的城市燃气管道第三方破坏引发泄漏流程模拟。
发明内容
有鉴于此,本发明提出一种燃气管道泄漏仿真模拟方法及装置,针对燃气管道不同的泄漏位置、风向等条件下的泄漏浓度,模拟泄漏过程。
本发明第一方面提供一种燃气管道泄漏仿真模拟方法,该方法包括:对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点;构建燃气管道泄漏扩散状态模型;设置燃气管道的泄漏参数;基于所述燃气管道的泄漏参数及所述燃气管道泄漏扩散状态模型,计算各观察点的泄漏气体浓度,生成显示信息,以根据所述显示信息可视化各观察点的泄漏气体浓度变化。
进一步的,当燃气管道内介质的压力与环境压力的比值大于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure BDA0003311737390000021
式中,P0为环境压力;P为燃气管道内介质的压力;k为气体的绝热指数;Cdg为气体泄漏系数;M为气体摩尔质量;R为气体常数;T为气体温度;A为泄漏点面积;r为泄漏点半径。
进一步的,当燃气管道内介质的压力与环境压力的比值小于等于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure BDA0003311737390000022
式中,P0为环境压力;P为燃气管道内介质的压力;k为气体的绝热指数;Cdg为气体泄漏系数;M为气体摩尔质量;R为气体常数;T为气体温度;A为泄漏点面积;r为泄漏点半径。
进一步的,所述对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点的步骤,包括:设置燃气管道上泄漏点位置;基于燃气管道上泄漏点位置,对燃气管道的各个不同区域进行不同的评价单元网格划分;构建各评价单元内燃气管道三维模型;根据燃气管道上泄漏点位置,在各评价单位内燃气管道三维模型的泄漏点上方设置多个不同观察点;在各评价单位内燃气管道三维模型的两个通风口中分别设置一观察点。。
进一步的,所述燃气管道的泄漏参数包括泄漏时长、通风条件及风速中一个或多个。
进一步的,所述计算所述燃气管道上各观察点的泄漏气体浓度,生成显示信息的步骤,包括:获取燃气管道的压力、泄漏点直径及气体温度;基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径、气体温度及泄漏参数,计算燃气管道泄漏气体的质量流量;根据各观察点的位置,计算各观察点体积;基于所述燃气管道泄漏气体的质量流量及各观察点体积,计算各观察点的泄漏气体浓度;根据各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
本发明第二方面提供一种燃气管道泄漏仿真模拟装置,该装置包括:管道三维模型建立模型,用于对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点;泄漏扩散状态模型构建模块,用于构建燃气管道泄漏扩散状态模型;燃气泄漏参数设置模块,用于设置燃气管道的泄漏参数;泄漏气体浓度模拟模块,用于基于所述燃气管道的泄漏参数及所述燃气管道泄漏扩散状态模型,计算各观察点的泄漏气体浓度,生成显示信息,以根据所述显示信息可视化各观察点的泄漏气体浓度变化。
进一步的,所述燃气管道三维模型对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点的步骤,包括:设置燃气管道上泄漏点位置;基于燃气管道上泄漏点位置,对燃气管道的各个不同区域进行不同的评价单元网格划分;构建各评价单元内燃气管道三维模型;根据燃气管道上泄漏点位置,在各评价单位内燃气管道三维模型的泄漏点上方设置多个不同观察点;在各评价单位内燃气管道三维模型的两个通风口中分别设置一观察点。。
进一步的,所述泄漏气体浓度模拟模块计算所述燃气管道上各观察点的泄漏气体浓度,生成显示信息的步骤,包括:获取燃气管道的压力、泄漏点直径及气体温度;基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径、气体温度及泄漏参数,计算燃气管道泄漏气体的质量流量;根据各观察点的位置,计算各观察点体积;基于所述燃气管道泄漏气体的质量流量及各观察点体积,计算各观察点的泄漏气体浓度。
进一步的,所述泄漏气体浓度模拟模块还用于:根据各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
上述的燃气管道泄漏仿真模拟方法及装置,首先,建立实际尺寸的燃气管道三维模型并划分网格,选择合适的观察点进行模拟,分析燃气管道被破坏后泄漏扩散情况;模拟设定燃气管道参数条件下的燃气管道被第三方破坏后的泄漏过程,分析通风条件对天然气扩散的影响情况。
附图说明
为了说明而非限制的目的,现在将根据本发明的优选实施例、特别是参考附图来描述本发明,其中:
图1是本申请一实施例提供的燃气管道泄漏仿真模拟方法的流程图;
图2是地下燃气管道三维模型的剖面图;
图3(a)是管道泄漏5s的气体浓度变化图;
图3(b)是管道泄漏20s的气体浓度变化图;
图3(c)是管道泄漏60s的气体浓度变化图;
图3(d)是管道泄漏180s的气体浓度变化图;
图3(e)是可爆炸气体浓度变化图;
图4(a)和图4(b)是管道泄漏180s的泄漏气体浓度随通风时间变化图;
图5是本申请另一实施例提供的燃气管道泄漏仿真模拟装置的结构图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施例对本发明进行详细描述。需要说明的是,在不冲突的情况下,本发明的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。
图1是本申请一实施例提供的燃气管道泄漏仿真模拟方法的流程图。该管道燃气泄漏仿真模拟方法针对燃气管道不同的泄漏位置、风向等条件下的泄漏浓度,模拟泄漏过程。
请参阅附图1,该燃气管道泄漏仿真模拟方法包括:
S100,对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点。
具体地,步骤S100的具体实现方式如下:
S101,设置燃气管道上泄漏点位置,基于燃气管道上泄漏点位置对燃气管道的各个不同区域进行不同的评价单元网格划分。
网格是构成运算的基本单元,会影响模拟结果的准确程度和运算时间,因此在燃气管道三维模型建立后,对燃气管道三维模型中各个区域进行不同的网格划分,并反复调整。网格的设定既要符合数值模型的要求,也要确保三维模型的科学性。在三维建模和划分网格时要相互协调,保证运算过程的顺利进行。
本实施例仅考虑不同通风情景对气体泄漏的影响,将燃气管道泄漏点位置不作为变量,本实施例将泄漏点位置设置为(1.2,22.3,1.5),将泄漏点附近区域设为核心区域,核心区域范围为min(1,22,1.3),max(1.5,22.5,1.8),将该区域内的网格细分为0.05m,在距泄漏点较远处将网格沿着min(0,4,0),max(3.5,57,10)方向拉伸。
S102,利用FLACS构建各评价单元内燃气管道三维模型。
以地下燃气管道为例,建立的地下燃气管道三维模型的剖面图如图2所述,综合管廊主体为2.5×2.2×60m3的矩形狭长舱室,每间隔50m设有人员进出口及通风口,管廊内部设有DN500mm的天然气管道,每隔10m有一个支撑柱。
S103,根据燃气管道上泄漏点位置,在各评价单元内燃气管道三维模型的泄漏点上方设置多个不同观察点,并在各评价单元内燃气管道三维模型的两个通风口中分别设置一观察点。
根据燃气管道泄漏点位置(1.2,22.3,1.5),在泄漏点上方设置了1个观察点P1,沿管道Y方向5、10、20、25、40、50、55m处设置观察点P1、P2、P3、P4、P5、P6、P7和P8,观察气体在管道内的扩散情况。分别在两个通风口中设置一个观察点P9,观察气体是否会扩散至通风口。
S200,构建燃气管道泄漏扩散状态模型。
基于燃气管道中天然气气体质量流量与质量流动状态的关系构建燃气管道泄漏扩散状态模型,其中,燃气管道泄漏扩散状态模型的表达式如下:
当燃气管道内介质的压力与环境压力的比值
Figure BDA0003311737390000061
时,气体被定义为亚音速流动,泄漏的质量流量计算公式为:
Figure BDA0003311737390000062
A=πr2 (9)
式中,P0为环境压力,Pa;P为燃气管道内介质的压力,Pa;k为气体的绝热指数,对于多原子气体天然气,取k=1.3;Cdg为气体泄漏系数,与泄漏点形状有关,本实施例取1.0;M为气体摩尔质量,取0.016kg/mol;R为气体常数,取8.3144J/(mol.K);T为气体温度,取293K;A为泄漏点面积,m2;r为泄漏点半径,m。
Figure BDA0003311737390000071
时,气体被定义为音速流动,泄漏的质量流量计算公式为:
Figure BDA0003311737390000072
S300,设置燃气管道的泄漏参数。
其中,燃气管道的泄漏参数包括泄漏时长、通风条件及风速等。
结合实际工况,根据城市综合管廊工程技术规范中通风系统的要求,设置5种场景进行模拟,5种场景下的泄漏参数见表1。
表1泄漏场景参数
场景 泄漏时长/s 通风 风速
1 180
2 180 1m/s
3 180 3m/s
4 180 泄漏120s后开启 1m/s
5 180 泄漏120s后开启 3m/s
S400,基于燃气管道的泄漏参数,计算燃气管道上各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度变化。
具体地,步骤S400的具体实现方式如下:
S401,获取燃气管道压力、泄漏点直径、泄漏点位置及气体温度,基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径及泄漏参数,计算燃气管道泄漏气体的质量流量。
以地下燃气管道为例,燃气管道压力为P=0.25MPa,泄漏点位置定为(1.2,22.3,1.5),泄漏点直径为20mm,计算泄漏孔的面积是0.000314m2
基于燃气管道泄漏扩散状态模型,利用燃气管道的压力、泄漏点面积、泄漏时间180s及温度20℃,计算燃气管道泄漏气体的质量流量0.1897kg/s。
S402,根据燃气管道上各观察点的位置,计算各观察点体积V。
S403,基于燃气管道泄漏气体的质量流量及各观察点体积,计算燃气管道各观察点的泄漏气体浓度。
以观察点Pi为例,观察点Pi的泄漏气体浓度Ci等于燃气管道泄漏气体的质量流量除以观察点Pi的体积。
S404,根据燃气管道各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
在本实施例中,显示信息可为曲线图,根据燃气管道各观察点的泄漏气体浓度,生成各观察点的泄漏气体浓度随时间变化曲线图。
以地下燃气管道,泄漏场景1(泄漏时长/s、无通风、无风速)为例,管道泄漏180s的气体浓度变化如图3(a)-图3(e)所示。泄漏开始时,由于管道压力作用使天然气从泄漏点喷出,喷射气体触碰到管廊顶板,发生“顶棚射流”现象,然后向两边扩散。中部位气体呈对称分布(图3(a)),随后,气体迅速向四周扩散与空气混合,并不断与壁面发生碰撞,使得气体流动方向发生改变,开始向管廊底部扩散。受管道阻碍,管道处聚集了部分浓度的泄漏气体。当气体运动到管廊底面到达管廊底部壁面后,气体流动的动能受到损失,速度减小,同时由于天然气密度低,又重新沿着两侧壁面流动到管廊上部,因此管廊上下部天然气浓度大于中下部天然气浓度,如图图3(b)所示。泄漏发生60s时,有少部分气体沿着通风口的缝隙扩散至人员进出口,使泄漏点左侧更容易扩散。随着泄漏时间的延长,到180s时管廊内充满了泄漏气体,人员进出口中泄漏气体的浓度也有所提高,此时管廊内部分泄漏气体浓度达到15%以上。已知天然气爆炸极限为5~15%,图3(e)为天然气浓度在5%至15%浓度范围内的气体云团分布,可以发现管廊内和人员进出口内都存在可爆炸气体,若遇火源将发生爆炸事故。此时管廊内的泄漏气体体积为20×2.5×2.2m3,人员进出口内泄漏气体体积为2×2.5×6.5。图4(a)和图4(b)为观察点泄漏气体浓度随时间的变化曲线,可以看出,泄漏180s后,P1、P2、P3、P4和P6观测泄漏点气体浓度在0.05~0.19之间,泄漏气体浓度在爆炸极限之内,人员进出口P9泄漏气体浓度达到0.08,高于爆炸下限。
上述的燃气管道泄漏仿真模拟方法,利用FLACS建立实际尺寸的燃气管道三维模型并划分网格,选择合适的观察点进行模拟,分析埋地管道被破坏后泄漏扩散情况;模拟设定燃气管道参数条件下的燃气管道被第三方破坏后的泄漏过程,分析通风条件对天然气扩散的影响情况。
图5是本申请另一实施例提供的燃气管道泄漏仿真模拟装置的流程图。该管道燃气泄漏仿真模拟装置针对燃气管道不同的泄漏位置、风向等条件下的泄漏浓度,模拟泄漏过程。
请参阅附图5,该燃气管道泄漏仿真模拟装置500包括:
燃气管道三维模型建立模型501,用于对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点。
具体地,燃气管道三维模型建立模型501对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点的具体实现方式如下:
燃气管道三维模型建立模型501设置燃气管道泄漏点位置,基于燃气管道泄漏点位置对燃气管道的各个不同区域进行不同的评价单元网格划分。
网格是构成运算的基本单元,会影响模拟结果的准确程度和运算时间,因此在燃气管道三维模型建立后,对燃气管道三维模型中各个区域进行不同的网格划分,并反复调整。网格的设定既要符合数值模型的要求,也要确保三维模型的科学性。在三维建模和划分网格时要相互协调,保证运算过程的顺利进行。
本实施例仅考虑不同通风情景对气体泄漏的影响,将燃气管道泄漏点位置不作为变量,本实施例将泄漏点位置设置为(1.2,22.3,1.5),将泄漏点附近区域设为核心区域,核心区域范围为min(1,22,1.3),max(1.5,22.5,1.8),将该区域内的网格细分为0.05m,在距泄漏点较远处将网格沿着min(0,4,0),max(3.5,57,10)方向拉伸。
燃气管道三维模型建立模型501利用FLACS构建各评价单元内燃气管道三维模型。
以地下燃气管道为例,建立的地下燃气管道三维模型的剖面图如图2所述,综合管廊主体为2.5×2.2×60m3的矩形狭长舱室,每间隔50m设有人员进出口及通风口,管廊内部设有DN500mm的天然气管道,每隔10m有一个支撑柱。
燃气管道三维模型建立模型501根据燃气管道泄漏点位置,在泄漏点上方设置多个不同观察点。
以燃气管道泄漏点位置(1.2,22.3,1.5)为例,在泄漏点上方设置了1个观察点P1,沿管道Y方向5、10、20、25、40、50、55m处设置观察点P1、P2、P3、P4、P5、P6、P7和P8,观察气体在管道内的扩散情况。分别在两个通风口中设置一个观察点P9,观察气体是否会扩散至通风口。
泄漏扩散状态模型构建模块502,用于构建燃气管道泄漏扩散状态模型。
泄漏扩散状态模型构建模块502基于燃气管道中天然气气体质量流量与质量流动状态的关系构建燃气管道泄漏扩散状态模型,其中,燃气管道泄漏扩散状态模型的表达式如下:
Figure BDA0003311737390000101
时,气体被定义为亚音速流动,泄漏的质量流量计算公式为:
Figure BDA0003311737390000102
式中,P0为环境压力,Pa;P为燃气管道内介质的压力,Pa;k为气体的绝热指数,对于多原子气体天然气,取k=1.3;Cdg为气体泄漏系数,与泄漏点形状有关,本实施例取1.0;M为气体摩尔质量,取0.016kg/mol;R为气体常数,取8.3144J/(mol.K);T为气体温度,取293K;A为泄漏点面积,m2;r为泄漏点半径,m。
Figure BDA0003311737390000103
时,气体被定义为音速流动,泄漏的质量流量计算公式为:
Figure BDA0003311737390000111
泄漏参数设置模块503,用于设置燃气管道的泄漏参数。
其中,燃气管道的泄漏参数包括泄漏时长、通风条件及风速等。
泄漏气体浓度模拟模块504,用于基于燃气管道的泄漏参数,计算燃气管道上各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度变化。
具体地,泄漏气体浓度模拟模块504基于燃气管道的泄漏参数,计算燃气管道上各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度变化的具体实现方式如下:
泄漏气体浓度模拟模块504获取燃气管道压力、泄漏点直径及泄漏点位置,基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径及泄漏参数,计算燃气管道泄漏气体的质量流量。
以地下燃气管道为例,燃气管道压力为P=0.25MPa,泄漏点位置定为(1.2,22.3,1.5),泄漏点直径为20mm,计算泄漏孔的面积是0.000314m2
基于燃气管道泄漏扩散状态模型,利用燃气管道的压力、泄漏点面积、泄漏时间180s及温度20℃,计算燃气管道泄漏气体的质量流量0.1897kg/s。
泄漏气体浓度模拟模块504根据燃气管道上各观察点的位置,计算各观察点体积V。
泄漏气体浓度模拟模块504根据燃气管道泄漏气体的质量流量及各观察点体积,计算燃气管道各观察点的泄漏气体浓度。
以观察点Pi为例,观察点Pi的泄漏气体浓度Ci等于燃气管道泄漏气体的质量流量除以观察点Pi的体积。
泄漏气体浓度模拟模块504根据燃气管道各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
上述的燃气管道泄漏仿真模拟装置,通过利用FLACS建立实际尺寸的燃气管道三维模型并划分网格,选择合适的观察点进行模拟,分析埋地管道被破坏后泄漏扩散情况;模拟设定燃气管道参数条件下的燃气管道被第三方破坏后的泄漏过程,分析通风条件对天然气扩散的影响情况。
上述具体实施方式,并不构成对本发明保护范围的限制。本领域技术人员应该明白的是,取决于设计要求和其他因素,可以发生各种各样的修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明保护范围之内。

Claims (6)

1.一种燃气管道泄漏仿真模拟方法,其特征在于,包括:
对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点;
构建燃气管道泄漏扩散状态模型,当燃气管道内介质的压力与环境压力的比值大于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure FDA0003884610980000011
当燃气管道内介质的压力与环境压力的比值小于等于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure FDA0003884610980000012
式中,P0为环境压力;P为燃气管道内介质的压力;k为气体的绝热指数;Cdg为气体泄漏系数;M为气体摩尔质量;R为气体常数;T为气体温度;A为泄漏点面积;r为泄漏点半径;
设置燃气管道的泄漏参数;
基于所述燃气管道的泄漏参数及所述燃气管道泄漏扩散状态模型,计算各观察点的泄漏气体浓度,生成显示信息,以根据所述显示信息可视化各观察点的泄漏气体浓度变化;
其中,所述对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点的步骤,包括:
设置燃气管道上泄漏点位置;基于燃气管道上泄漏点位置,对燃气管道的各个不同区域进行不同的评价单元网格划分;构建各评价单元内燃气管道三维模型;根据燃气管道上泄漏点位置,在各评价单位内燃气管道三维模型的泄漏点上方设置多个不同观察点;在各评价单位内燃气管道三维模型的两个通风口中分别设置一观察点。
2.根据权利要求1所述的燃气管道泄漏仿真模拟方法,其特征在于,所述燃气管道的泄漏参数包括泄漏时长、通风条件及风速中一个或多个。
3.根据权利要求1所述的燃气管道泄漏仿真模拟方法,其特征在于,所述计算所述燃气管道上各观察点的泄漏气体浓度,生成显示信息的步骤,包括:
获取燃气管道的压力、泄漏点直径及气体温度;
基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径、气体温度及泄漏参数,计算燃气管道泄漏气体的质量流量;
根据各观察点的位置,计算各观察点体积;
基于所述燃气管道泄漏气体的质量流量及各观察点体积,计算各观察点的泄漏气体浓度;
根据各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
4.一种燃气管道泄漏仿真模拟装置,其特征在于,包括:
燃气管道三维模型建立模块,用于对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点;
泄漏扩散状态模型构建模块,用于构建燃气管道泄漏扩散状态模型,当燃气管道内介质的压力与环境压力的比值大于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure FDA0003884610980000031
当燃气管道内介质的压力与环境压力的比值小于等于设置阈值时,所述燃气管道泄漏扩散状态模型为:
Figure FDA0003884610980000032
式中,P0为环境压力;P为燃气管道内介质的压力;k为气体的绝热指数;Cdg为气体泄漏系数;M为气体摩尔质量;R为气体常数;T为气体温度;A为泄漏点面积;r为泄漏点半径;
泄漏参数设置模块,用于设置燃气管道的泄漏参数;
泄漏气体浓度模拟模块,用于基于所述燃气管道的泄漏参数及所述燃气管道泄漏扩散状态模型,计算各观察点的泄漏气体浓度,生成显示信息,以根据所述显示信息可视化各观察点的泄漏气体浓度变化;
其中,所述燃气管道三维模型建立模块对燃气管道进行评价单元网格划分,建立各评价单元内燃气管道三维模型,并对各评价单元内燃气管道三维模型设置观察点的步骤,包括:
设置燃气管道上泄漏点位置;基于燃气管道上泄漏点位置,对燃气管道的各个不同区域进行不同的评价单元网格划分;构建各评价单元内燃气管道三维模型;根据燃气管道上泄漏点位置,在各评价单位内燃气管道三维模型的泄漏点上方设置多个不同观察点;在各评价单位内燃气管道三维模型的两个通风口中分别设置一观察点。
5.根据权利要求4 所述的燃气管道泄漏仿真模拟装置,其特征在于,所述泄漏气体浓度模拟模块计算所述燃气管道上各观察点的泄漏气体浓度,生成显示信息的步骤,包括:
获取燃气管道的压力、泄漏点直径及气体温度;
基于燃气管道泄漏扩散状态模型及燃气管道的压力、泄漏点直径、气体温度及泄漏参数,计算燃气管道泄漏气体的质量流量;
根据各观察点的位置,计算各观察点体积;
基于所述燃气管道泄漏气体的质量流量及各观察点体积,计算各观察点的泄漏气体浓度。
6.根据权利要求5所述的燃气管道泄漏仿真模拟装置,其特征在于,所述泄漏气体浓度模拟模块还用于:
根据各观察点的泄漏气体浓度,生成显示信息,以根据显示信息可视化各观察点的泄漏气体浓度随时间变化趋势。
CN202111218747.0A 2021-10-20 2021-10-20 燃气管道泄漏仿真模拟方法及装置 Active CN114239193B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111218747.0A CN114239193B (zh) 2021-10-20 2021-10-20 燃气管道泄漏仿真模拟方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111218747.0A CN114239193B (zh) 2021-10-20 2021-10-20 燃气管道泄漏仿真模拟方法及装置

Publications (2)

Publication Number Publication Date
CN114239193A CN114239193A (zh) 2022-03-25
CN114239193B true CN114239193B (zh) 2023-02-03

Family

ID=80743142

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111218747.0A Active CN114239193B (zh) 2021-10-20 2021-10-20 燃气管道泄漏仿真模拟方法及装置

Country Status (1)

Country Link
CN (1) CN114239193B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115691705B (zh) * 2022-11-11 2023-07-25 武汉数信科技有限公司 一种气体泄漏扩散监测方法及系统
CN118517643A (zh) * 2024-07-25 2024-08-20 海纳云物联科技有限公司 用于燃气泄漏定位的方法及装置、计算机可读存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590336A (zh) * 2017-09-13 2018-01-16 哈尔滨理工大学 燃气管道泄漏对内部流场影响的数值模拟方法
CN108038275A (zh) * 2017-11-28 2018-05-15 哈尔滨理工大学 一种天然气管道泄漏声场的数值仿真及特性分析方法
CN108954016A (zh) * 2017-11-09 2018-12-07 北京市燃气集团有限责任公司 基于增强现实技术的燃气管道泄漏处置系统
CN110630909A (zh) * 2019-09-04 2019-12-31 常州大学 一种燃气管道射流火焰作用范围预测方法
CN112070273A (zh) * 2020-08-11 2020-12-11 深圳市燃气集团股份有限公司 一种燃气扩散预测方法、智能终端及存储介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111611723B (zh) * 2020-06-05 2023-05-23 中国石化销售股份有限公司华南分公司 城市建筑群输送管道泄漏扩散数据生成方法、装置及介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590336A (zh) * 2017-09-13 2018-01-16 哈尔滨理工大学 燃气管道泄漏对内部流场影响的数值模拟方法
CN108954016A (zh) * 2017-11-09 2018-12-07 北京市燃气集团有限责任公司 基于增强现实技术的燃气管道泄漏处置系统
CN108038275A (zh) * 2017-11-28 2018-05-15 哈尔滨理工大学 一种天然气管道泄漏声场的数值仿真及特性分析方法
CN110630909A (zh) * 2019-09-04 2019-12-31 常州大学 一种燃气管道射流火焰作用范围预测方法
CN112070273A (zh) * 2020-08-11 2020-12-11 深圳市燃气集团股份有限公司 一种燃气扩散预测方法、智能终端及存储介质

Also Published As

Publication number Publication date
CN114239193A (zh) 2022-03-25

Similar Documents

Publication Publication Date Title
CN114239193B (zh) 燃气管道泄漏仿真模拟方法及装置
CN107590336B (zh) 燃气管道泄漏对内部流场影响的数值模拟方法
Lu et al. Study on leakage and ventilation scheme of gas pipeline in tunnel
CN108038275A (zh) 一种天然气管道泄漏声场的数值仿真及特性分析方法
Walton et al. Large-eddy simulation of pollution dispersion in an urban street canyon—Part I: comparison with field data
CN111611723B (zh) 城市建筑群输送管道泄漏扩散数据生成方法、装置及介质
Toja-Silva et al. On the urban geometry generalization for CFD simulation of gas dispersion from chimneys: Comparison with Gaussian plume model
CN203533216U (zh) 燃气管网泄漏源实时定位分析系统
EP4195091A1 (en) Target available model-based environment prediction method and apparatus, program, and electronic device
CN103196038A (zh) 燃气管网泄漏源实时定位分析方法及系统
CN107505267B (zh) 一种气体探测器布点分析方法及装置
CN102902860A (zh) 基于cfd技术的工作场所职业暴露模拟分析方法
Dai et al. Evaluation of fast fluid dynamics with different turbulence models for predicting outdoor airflow and pollutant dispersion
CN104361157B (zh) 一种建筑物间风环境的评价方法
CN114117732A (zh) 燃气管道泄漏后爆炸模拟方法及装置
CN109614638A (zh) 一种非直接建模的城市风环境cfd模拟方法
Meroney et al. Virtual reality in wind engineering: the windy world within the computer. J
KR101907248B1 (ko) 확률적 가스폭발 시나리오 산출 시스템 및 이를 사용한 확률적 가스폭발 시나리오 산출 방법
CN110879919B (zh) 一种爆炸作用下毒物扩散的分段式模拟方法
Ejaz et al. Numerical evaluation of separation efficiency in the diverging T-junction for slug flow
Salamonowicz et al. Numerical modelling of dispersion of ammonia and chlorine in urban areas during emergency accident
CN118313312A (zh) 燃气管道泄漏模拟方法、装置及电子设备
CN114330050A (zh) 抽水蓄能电站地下厂房通风廊道的温湿度数值模拟方法
Li et al. Investigation of urban natural gas pipeline leak and resulting dispersion in a semi-closed space: A case of accident in Shiyan, China
CN116341231B (zh) 人工增雨催化剂扩散轨迹的模拟方法、装置、服务器和存储介质

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