中国安全科学学报

工况序号	泄漏孔径/ mm	泄漏流量/ (kg/s)	等效直径/ mm	环境风速/ (m/s)	环境风向	风机风速/ (m/s)
1	2	0.068 8	23.25	0	—	—
2	2	+X	—
3	5	+X	—
4	8	+X	—
5	0	—	5
6	0	—	10
7	5	0.430 0	58.12	0	—	—
8	10	1.720 3	116.24	0	—	—

工况序号	泄漏孔径/ mm	泄漏流量/ (kg/s)	等效直径/ mm	环境风速/ (m/s)	环境风向	风机风速/ (m/s)
1	2	0.068 8	23.25	0	—	—
2	2	+X	—
3	5	+X	—
4	8	+X	—
5	0	—	5
6	0	—	10
7	5	0.430 0	58.12	0	—	—
8	10	1.720 3	116.24	0	—	—

曲线序号	坐标
1	(9.5, 7.6, 108)	(24, 7.6, 108)
2	(9.5, 7.6, 106)	(9.5, 7.6, 106)
3	(9.5, 7.6, 104)	(24, 7.6, 104)

曲线序号	坐标
1	(9.5, 7.6, 108)	(24, 7.6, 108)
2	(9.5, 7.6, 106)	(9.5, 7.6, 106)
3	(9.5, 7.6, 104)	(24, 7.6, 104)

加氢机氢气泄漏扩散数值模拟研究

PDF下载

谭路遥 , 姚勇征 ^** , 潘傲澜 , 胡茂炜

中国安全科学学报 | 安全技术与工程 2026,36(5): 224-233

收起

中国安全科学学报 |安全技术与工程 2026 , 36 (5) : 224 -233

加氢机氢气泄漏扩散数值模拟研究

全屏

[Author(id=1278415586585723652, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, orderNo=0, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=tanluyao3915@126.com, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1278415586749301510, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415586585723652, language=EN, stringName=Luyao Tan, firstName=Luyao, middleName=null, lastName=Tan, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1278415586904490759, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415586585723652, language=CN, stringName=谭路遥, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083, bio={"img":"luOl8Py9q3Sw76YzDK23zw==","content":"

谭路遥 (2000—),女,湖南株洲人,博士研究生,主要研究方向为加氢站氢泄漏燃爆风险防控等。E-mail:tanluyao3915@126.com。

"}, bioImg=luOl8Py9q3Sw76YzDK23zw==, bioContent=

谭路遥 (2000—),女,湖南株洲人,博士研究生,主要研究方向为加氢站氢泄漏燃爆风险防控等。E-mail:tanluyao3915@126.com。

, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1278415584794755840, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, xref=null, ext=[AuthorCompanyExt(id=1278415584807338753, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China), AuthorCompanyExt(id=1278415584840893186, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083)])]), Author(id=1278415586979988233, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, orderNo=1, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=yaoyz@cumtb.edu.cn, emailSecond=null, emailThird=null, correspondingAuthor=1, authorType=1, ext={EN=AuthorExt(id=1278415587269395211, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415586979988233, language=EN, stringName=Yongzheng Yao, firstName=Yongzheng, middleName=null, lastName=Yao, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=^**, address=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1278415587328115468, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415586979988233, language=CN, stringName=姚勇征, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=^**, address=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1278415584794755840, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, xref=null, ext=[AuthorCompanyExt(id=1278415584807338753, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China), AuthorCompanyExt(id=1278415584840893186, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083)])]), Author(id=1278415587412001550, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, orderNo=2, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1278415587512664848, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415587412001550, language=EN, stringName=Aolan Pan, firstName=Aolan, middleName=null, lastName=Pan, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1278415587705602833, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415587412001550, language=CN, stringName=潘傲澜, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1278415584794755840, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, xref=null, ext=[AuthorCompanyExt(id=1278415584807338753, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China), AuthorCompanyExt(id=1278415584840893186, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083)])]), Author(id=1278415587785294611, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, orderNo=3, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=null, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1278415588120838933, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415587785294611, language=EN, stringName=Maowei Hu, firstName=Maowei, middleName=null, lastName=Hu, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1278415588208919318, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, authorId=1278415587785294611, language=CN, stringName=胡茂炜, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1278415584794755840, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, xref=null, ext=[AuthorCompanyExt(id=1278415584807338753, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China), AuthorCompanyExt(id=1278415584840893186, tenantId=1146029695717560320, journalId=1146031787341344770, articleId=1278415572169900763, companyId=1278415584794755840, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083)])])]

谭路遥, 姚勇征^**, 潘傲澜, 胡茂炜

作者信息

中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083

通讯作者:

^** 姚勇征(1991—),男,河北承德人,博士,副教授,主要从事新能源风险防控研究。E-mail: yaoyz@cumtb.edu.cn。

作者简介:

谭路遥 (2000—),女,湖南株洲人,博士研究生,主要研究方向为加氢站氢泄漏燃爆风险防控等。E-mail:tanluyao3915@126.com。

Numerical simulation of hydrogen diffusion from hydrogen dispenser leakage

Luyao Tan, Yongzheng Yao^**, Aolan Pan, Maowei Hu

Affiliations

School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China

出版时间: 2026-05-28 doi: 10.16265/j.cnki.issn1003-3033.2026.05.0886

文章导航

摘要

收起

为降低加氢机加注过程中因泄漏引发的事故风险,利用Fluent软件模拟35 MPa加氢机氢泄漏扩散过程,研究加氢区罩棚结构下氢泄漏扩散特性,分析泄漏孔径、风速及局部通风等因素对氢气摩尔分数、可燃区域分布的影响。研究结果表明:当加氢软管泄漏孔径为2 mm时,不会在罩棚下表面形成可燃区域,而当泄漏孔径达到5或10 mm时,可在罩棚下表面形成可燃区域,且罩棚下表面摩尔分数最高点在平行于射流中轴线附近。其中,泄漏孔径为10mm,垂直泄漏方向罩棚下侧摩尔分效呈现高斯分布;当环境风向垂直于泄漏方向,风速为2或8 m/s时,可有效减少泄漏加氢机附近氢气积聚;而风速为5 m/s时,在障碍物附近产生涡旋结构,诱导氢气积聚,增加燃爆风险;无环境风条件下,在加氢区加入风机进行局部通风,5或10 m/s均可在2 s内将流场内氢气云摩尔分数稀释至可燃范围以下,而10 m/s的风机风速对稀释泄漏源前方氢气摩尔分数作用更大。

关键词

加氢机 / 氢气泄漏 / 氢气扩散 / 数值模拟 / 泄漏孔径 / 环境风速 / 局部通风

Abstract

收起

In order to reduce accident risks caused by leakage during the refueling process, the hydrogen leakage and diffusion behavior of a 35 MPa hydrogen dispenser was numerically simulated using Ansys Fluent. The characteristics of hydrogen leakage and diffusion under the canopy structure in the refueling zone were investigated. The effects of leakage diameters, ambient wind velocity, and local ventilation on hydrogen concentration distribution and the evolution of flammable areas were analyzed. The results show that when the leakage diameter of filling hose is 2 mm, no flammable area is formed on the underside of canopy. However, when the leakage diameter increases to 5 mm and 10 mm, a flammable area can develop on canopy underside. The location of the highest hydrogen concentration on underside of canopy is concentrated near the axis parallel to the jet direction. Specifically, when leakage diameter is 10 mm, the hydrogen concentration on underside of canopy along vertical leakage direction exhibits a Gaussian distribution. When the ambient wind is perpendicular to leakage direction, wind velocities of 2 m/s and 8 m/s can effectively reduce hydrogen accumulation near the leakage hydrogen dispenser. In contrast, at the wind speed of 5 m/s, a vortex structure was formed near obstacles, leading to hydrogen accumulation and increasing the risk of fire and explosion. Under no ambient wind conditions, local ventilation is provided in the refueling zone. When the ventilation velocity reaches 5 m/s and 10 m/s, the hydrogen cloud concentration within the flow field can be successfully diluted to below the flammable limit within 2 s. Moreover, a ventilation velocity of 10 m/s shows a more pronounced effect in reducing the hydrogen concentrations in front of the leakage source.

Key words

hydrogen dispenser / hydrogen leakage / hydrogen diffusion / numerical simulation / leakage diameter / ambient wind velocities / local ventilation

引用本文

谭路遥, 姚勇征, 潘傲澜, 胡茂炜. 加氢机氢气泄漏扩散数值模拟研究. 中国安全科学学报, 2026 , 36 (5) : 224 -233 . DOI: 10.16265/j.cnki.issn1003-3033.2026.05.0886

Luyao Tan, Yongzheng Yao, Aolan Pan, Maowei Hu. Numerical simulation of hydrogen diffusion from hydrogen dispenser leakage[J]. China Safety Science Journal, 2026 , 36 (5) : 224 -233 . DOI: 10.16265/j.cnki.issn1003-3033.2026.05.0886

正文

收起

0 引言

收起

氢能作为一种清洁、高效的可再生能源,是全球“脱碳时代”的能源新焦点,已被30多个经济体纳入能源发展体系^[1-3]。自2020年我国提出“碳达峰、碳中和”目标以来,大力推进以氢能为代表的绿色交通体系建设。加氢站是绿色交通建设的重要基础设施,截至2024年6月,我国(含港澳台地区)建成加氢站507座,位居全球首位^[4]。加氢站中加氢机运行压力高,易受氢脆效应、疲劳损伤及人员不当操作影响,导致氢气泄漏发生。加之加氢区多为平顶罩棚,形成受限区域,造成氢气积聚,有扩大成燃爆事故的可能。因此,研究加氢区氢气泄漏扩散特性,有利于减少加注过程中的事故风险,助力加氢站安全运营。

在泄漏扩散规律研究方面,国内外学者多采用数值模拟方法研究泄漏条件、环境因素对氢气泄漏扩散的影响。戴君豪等^[5]采用Flacs研究了加氢站长管拖车储氢温度对氢气扩散的影响,发现储氢温度越低,氢射流在撞击隔爆墙后,可燃氢气云达到稳定的时间越长,扩散范围也逐渐增大。Qian Jinyuan等^[6]采用Fluent软件研究了加氢站储氢罐泄漏位置与环境风效应对氢气泄漏扩散影响,发现泄漏位置与障碍物之间越窄,可燃氢气云轮廓越不规则,与泄漏相反的风向可能导致可燃气云轮廓增大。在此基础上,Li Xiaojuan等^[7]进一步研究了可变风速与恒定风速对可燃氢气云轮廓与氢浓度分布影响。Liang Yang等^[8]针对储氢罐、长管拖车等泄漏场景,研究了泄漏方向、风向及风速等因素对事故后果的影响。大多数学者考虑到试验危险性,采用氦气代替氢气开展缩尺试验研究。Swain等^[9]通过试验验证了采用氦气代替氢气试验的准确性和可行性。田苗等^[10]搭建了缩尺试验平台,研究停车场封闭条件、泄漏流量、泄漏孔径对氢气扩散的影响。牛奕等^[11]搭建了加氢站缩尺寸试验,利用氦气代替氢气,研究泄漏量、环境风速对氢气泄漏扩散影响。为降低氢气泄漏引发燃爆事故的可能性,部分学者采用通风技术降低氢气泄漏危险。Tamura等^[12]采用鼓风机通风降低燃料电池乘用车周围氢气摩尔分数,发现鼓风机风速达到10 m/s时,可有效降低车辆周围氢气摩尔分数。Xie Hong等^[13]使用Fluent研究鼓风机对氢燃料电池汽车泄漏降低氢气摩尔分数作用,发现体积流量740 m³/min的鼓风机对于2 000 NL/min氢气泄漏率稀释作用良好。

综上所述,现有研究多集中于长管拖车、储氢罐泄漏,对加氢机氢气泄漏扩散行为研究较少,且缺乏加氢区构筑物等对气流扩散影响分析。此外,多数学者主要分析环境风速、风向对氢气云稀释作用,对于加氢站泄漏场景,应用局部通风降低氢气泄漏危险性的研究不足。因此,笔者拟以典型高压气态加氢站为研究对象,采用Fluent数值模拟软件研究加氢区罩棚结构下氢气扩散行为,分析泄漏孔径、环境风速及机械通风等因素对氢气扩散的影响,以期为加氢站建设、日常运维等提供理论依据。

1 加氢机氢气泄漏扩散模型

收起

1.1 物理模型

基于北京市某加氢站实际空间布局,构建加氢站全尺寸模型,计算域大小为L(长度)×W(宽度)×H(高度)=120 m×60 m×12 m,如图1所示。加氢站内主要设备及构筑物如图2所示。为减少模型的复杂程度及网格数量,简化加氢站内主要的涉氢设备,针对加氢软管泄漏场景,假设加氢枪连接端的加氢软管发生泄漏,泄漏处等效为圆孔,泄漏位置圆心坐标为(17.78,1.8,92.96)m。在泄漏源正前方,相隔9.4 m的位置有一加氢机,该加氢机中心点坐标为(17.58,1.45,102.66)m。另外,在泄漏源斜前方位置还有一罩棚支柱,支柱中心点坐标为(15.38,3.95,102.58)m。

1.2 数值计算模型

氢气泄漏扩散满足质量守恒、动量守恒、能量守恒,也遵循组分守恒定律^[14]。泄漏扩散过程受到障碍物影响,流动过程复杂,属于湍流流动。由于Realizable k-ε湍流模型可较好地预测圆孔射流,因此选择Realizable k-ε湍流模型^[15]。此外,采用半隐式压力耦合方程算法,对流项采用二阶迎风格式。

质量守恒方程为:

(1)$\frac{\partial \rho }{\partial t}+▽·\left(\rho v\right)=0$

式中:ρ为流体密度,kg/m³;t为时间,s;v为速度矢量,m/s。

动量守恒方程为:

(2)$ \frac{\partial(\rho \boldsymbol{v})}{\partial t}+\nabla \cdot(\rho \boldsymbol{v} \boldsymbol{v})=-\nabla p+\nabla \cdot \overline{\bar{\tau}}+\rho \boldsymbol{g}$

式中:p为静压,Pa;g为重力加速度,m/s²;$\stackrel{̿}{\tau }$为应力张量,Pa。

能量守恒方程为:

(3)$ \begin{array}{c}\nabla \cdot(\boldsymbol{v}(\rho E+p))= \\\nabla \cdot\left(k_{e f f} \nabla T-\sum_{n} h_{n} \boldsymbol{J}_{n}+\left(\overline{\bar{\tau}}_{e f f} \cdot \boldsymbol{v}\right)\right)\end{array}$

式中:E为单位质量流体的能量,J/kg;k_e为有效传热系数,W/(m·K);T为流体的热力学温标,K;h_n为组分n的单位质量焓,J/kg;$\mathrm{ }\stackrel{\text{'}}{ }$ $ \boldsymbol{J}_{n}$为组分n的质量扩散通量,kg/(m²·s);${\stackrel{̿}{\tau }}_{e}$为有效应力张量,Pa。

组分输运方程为:

(4)$ \frac{\partial\left(\rho Y_{n}\right)}{\partial t}+\nabla \cdot\left(\rho \boldsymbol{v} || Y_{n}\right)=-\nabla \cdot \boldsymbol{J}_{n}$

(5)$ \boldsymbol{J}_{n}=-\rho\left(D_{n, m}+\frac{\mu_{t}}{\rho S c}\right) \nabla Y_{n}$

式中:Y_n为组分n的质量分数,%;D_n,m 为流体扩散系数,m²/s;氢气取6.1×10^-5m²/s;Sc为湍流施密特数,取0.7。

Realizable k-ε湍流模型及k的控制方程见下式:

(6)$\begin{array}{c}\frac{\partial(\rho k)}{\partial t}+\nabla \cdot(\rho k \boldsymbol{v})=\nabla \cdot\left[\left(\mu+\frac{\mu_{t}}{\sigma_{k}}\right) \nabla k\right]+ \\G_{k}+G_{b}-\rho \varepsilon-Y_{M}\end{array}$

式中:k为湍流动能,m²/s²;ε为湍流耗散率,m²/s³;Pr为k的普朗特数,取1.0;G_k为平均速度梯度产生的湍流动能,kg/(m·s³);G_b为浮力产生的湍流动能,kg/(m·s³);Y_M为可压缩湍流中脉动扩张的贡献项,kg/(m·s³)。

ε的控制方程为:

(7)$\begin{array}{c}\frac{\partial(\rho \varepsilon)}{\partial t}+\nabla \cdot(\rho \varepsilon v)=\nabla \cdot\left[\left(\mu+\frac{\mu_{t}}{\sigma_{\varepsilon}}\right) \nabla \varepsilon\right]+\rho C_{1} S \varepsilon \\+C_{1 \varepsilon} \frac{\varepsilon}{k} C_{3 \varepsilon} G_{b}-\rho C_{2} \frac{\varepsilon^{2}}{k+\sqrt{\nu \varepsilon}}\end{array}$

式中:μ为流体分子动力黏度,Pa·s;μ_t湍流黏度,Pa·s;Pr为ε的普朗特数,取1.2;C_1ε、C₂为模型常数,分别取1.44,1.90;C₁为经验常数,取值范围为0.43~1;S为平均应变率张量的模,1/s;ν为湍流运动黏度,m²/s;C_3ε为浮力相关常数。

为验证以上数值计算模型的准确性,将数值模拟结果与Pitts等^[16]的试验数据对比,如图3所示。由对比结果可知:氢气摩尔分数相差较小,因此,所使用的数值计算模型可准确计算氢气泄漏扩散过程。

1.3 边界条件设置

将泄漏位置设置为质量流量入口,泄漏方向沿Z轴正方向。当考虑环境风速时,将加氢站左侧设置为速度入口,其他方向设置为压力出口,加氢站其他不发生泄漏部分均设置为壁面。

由于实际氢气泄漏孔尺寸与加氢站工作单元尺寸差距大,采用泄漏孔实际尺寸建模,将导致网格数量大。其次,高压氢气泄漏到环境中属于欠膨胀射流过程,在泄漏孔外形成马赫盘等复杂激波结构,在计算过程中需采用较小的松弛因子和时间步长^[17]。为简化数值模拟过程,采用Birch等^[18]提出的虚喷嘴模型计算激波结构后泄漏氢气的伪直径、质量流量等。

氢气泄漏扩散受泄漏孔径、风速等影响较大,因此,考虑3种泄漏孔径,并研究3种不同环境风速和2种不同局部通风风速的影响,共设置以下8种不同泄漏场景,见表1。

1.4 网格独立性验证

为减少网格数量对计算结果的影响,采用非结构化网格划分计算域,共划分4组网格,网格数量分别为2 448 628、2 904 858、3 260 626、4 507 210。在计算域中随机选取(17.78,1.8,95)m作为监测点,使用上述数值计算模型,得到不同网格数量下监测点的氢摩尔分数变化。网格独立性验证如图4所示,当网格数量为3 260 626与4 507 210时,所对应曲线氢气摩尔分数相差较小,考虑到节省计算资源,选择网格数量为3 260 626进行后续计算模拟,泄漏孔采用局部面加密,网格尺寸为0.002 m,泄漏源周围5 m×5 m×5 m的位置,采用影响体(Body of Influence)加密,加密尺寸为0.08 m。

2 氢气泄漏数值模拟结果分析

收起

2.1 泄漏孔径对氢气泄漏扩散影响

在无环境风条件下,选取泄漏孔径d为2、5、10 mm,研究泄漏孔径对氢气泄漏扩散的影响(图5),泄漏25 s后,泄漏孔径5与10 mm可燃区域可扩大至加氢区范围以外,而泄漏孔径为2 mm,由于泄漏量较少且氢射流初始动量较小,在撞击前方障碍物后动量迅速衰减,可燃区域主要位于前侧加氢机和罩棚支柱附近。

不同泄漏孔径可燃氢气云体积如图6所示,泄漏25 s后,泄漏孔径越大,产生可燃氢气云体积越大。d=2 mm共产生可燃氢气云体积19.39 m³,d=5 mm产生296.46 m³,d=10 mm产生1 155.23 m³。泄漏孔径10 mm可燃氢气云在5~10 s间先呈现下降趋势,10~25 s重新恢复增长趋势。可燃氢气云“先降后升”趋势是泄漏前期氢射流近场稀释作用与后期罩棚区域氢气积聚所致。

泄漏孔径10 mm下氢气摩尔分数分布(t=5~10 s)如图7a—图7c所示,近场区域由于湍流强度大,氢气浓度被迅速稀释至可燃下限以下,因此,可燃氢气云体积先呈阶段性下降趋势。泄漏孔径10 mm下氢气摩尔分数分布(t=15~25 s)如图7d—图7f所示,随着泄漏持续进行,泄漏氢气撞击前方加氢机及支柱后,开始向罩棚下扩散并积聚,形成较多的可燃氢气云,另外有部分氢气云扩散到罩棚区域外,不再受罩棚区域限制,持续卷吸空气,致使远场可燃氢气云体积持续增加。因此,t=10 s后可燃氢气云体积呈现上升趋势。

35 MPa加氢机泄漏后,泄漏氢气在动量主导下快速撞击到前侧障碍物(加氢机与罩棚支柱),射流方向发生偏转,泄漏氢气向罩棚下侧扩散。为研究氢射流撞击障碍物后,罩棚下侧氢气摩尔分数变化,在加氢区罩棚下侧壁面处设置3条监测曲线,曲线设置如图8所示,曲线坐标见表2。

泄漏孔径2 mm,罩棚下侧监测曲线浓度趋势总体呈现出“左侧高-中间次之-右侧低”趋势,并沿射流中轴线方向表现为距泄漏孔越远,氢浓度越低。2 mm泄漏孔径在扩散过程中,动量较弱,浮力主导作用强。泄漏孔径2 mm监测曲线氢气摩尔分数如图9所示,氢射流撞击前方加氢机后,动量削减并向周围扩散。左侧为罩棚支柱,进一步削弱了泄漏氢气动量,氢气在左侧滞留时间久,且积累量大。而右侧无障碍物阻挡,不易在罩棚下侧形成积聚,因此,形成了“左侧高-中间次之-右侧低”的摩尔分数趋势。

泄漏孔径5 mm,罩棚下侧Z=106、108 m摩尔分数分布呈中间高-左边次之-右边低,而Z=104 m,则表现为左侧高-中间次之-右侧最低,均表现为沿射流中轴线方向离泄漏孔越远,氢摩尔分数越高。泄漏孔径5 mm监测曲线氢气摩尔分数如图10所示,5 mm泄漏孔径形成的氢射流量较大,氢射流在撞击前方加氢机后,仍保持较大动量向前扩散,在浮力作用的共同作用下,泄漏氢气在泄漏源正前方罩棚下积聚,沿壁面向左右两侧横向扩散。因左侧由于罩棚支柱阻碍作用,氢气摩尔分数略高于右侧。

泄漏孔径10 mm,监测曲线氢气摩尔分数均呈现高斯分布,且沿射流中轴线方向大致表现为距泄漏孔越远,摩尔分数越高。泄漏孔径10 mm形成的氢射流动量远大于前2种孔径,受障碍物阻碍作用减小。泄漏孔径10 mm监测曲线摩尔分数分布如图11所示,氢射流撞击前侧加氢机后,大量氢气在中间罩棚下表面积聚,氢气摩尔分数达到10%以上,氢气摩尔分数沿罩棚左右两侧衰减。此外,10 mm泄漏孔径形成的氢射流,湍流强度大,加速了与周围空气的混合,使罩棚下表面氢气摩尔分数分布趋于均匀。与d=2 mm相比,d=5和10 mm时,罩棚下侧监测曲线沿射流中轴线方向表现为距泄漏源越远,罩棚下侧氢气摩尔分数越高。这是由于氢射流在撞击前方加氢机后,受动量主导作用大于浮力主导作用,随着氢射流向前扩散,动量逐渐衰减,浮力主导作用增强,氢气开始停滞,使得下游区域氢气摩尔分数高于近源区域。因此,5和10 mm泄漏孔径,沿射流中轴线方向大致表现为距泄漏源越远,罩棚下侧氢气摩尔分数越大。

2.2 环境风速对氢气泄漏扩散影响

以2 mm泄漏孔径为例,设置垂直泄漏方向的环境风,即沿X轴正方向,环境风速大小分别为2、5、8 m/s,分析不同环境风速对氢气扩散的影响。t=25 s,不同环境风速作用下可燃区域分布,如图12所示。v=0 m/s,氢射流撞击前方加氢机,在加氢机之间形成可燃区域,还有部分氢气受罩棚支柱阻碍,在罩棚支柱一侧形成了可燃区域。沿X轴正方向环境风速为2 m/s,削弱了氢射流沿Z轴方向动量,浮力主导作用增强,可燃区域高度抬升,并向前方加氢机与罩棚支柱间偏移。当环境风速增大至5 m/s,氢气射流在向上扩散过程中,撞击到前方加氢机,使可燃区域被前方加氢机限制在一侧,未继续扩大到前方。沿X轴正方向环境风速增大至8 m/s,可燃区域的高度上升,未被前方加氢机限制,向罩棚支柱一侧倾斜。

为研究不同环境下的风速,泄漏源附近气流运动情况,截取t=25 s,泄漏孔横截面(X=17.78 m)的不同环境风速作用下流线分布,如图13所示。当沿X轴正方向环境风速为2 m/s时,大部分可燃氢气云扩散高度超过前方加氢机,再向罩棚下侧扩散。环境风速为8 m/s时,泄漏氢气沿Z轴正方向动量减小较多,未扩散到泄漏源前侧加氢机处,而在浮力主导下向上扩散。当沿X轴正方向环境风速为5 m/s时,泄漏氢气沿Z轴正方向动量被削减后,扩散高度仍可超过前方加氢机,与环境风速2 m/s相比,继续向前扩散的氢气较少,大部分泄漏氢气都积聚在前方加氢机一侧。这是由于在泄漏源前方加氢机处形成回流区,诱导氢气积聚在泄漏源与前方加氢机之间,浓度达到可燃下限,增大了燃爆风险。前方加氢机附近回流区的形成是由于泄漏氢气扩散速度快,在泄漏源与前方加氢机之间形成负压区。而5 m/s的环境风速(沿X轴正方向)削弱了泄漏氢气沿Z轴正方向动量,在浮力的主导下向上扩散,与向负压区运动的外界气流相遇,卷吸外界气流往下运动,被前方加氢机阻挡后沿加氢机侧壁向上扩散,形成气流循环,产生回流区。因此,沿X轴正方向5 m/s的环境风速增加了泄漏源前侧加氢机附近的燃爆风险。

2.3 局部通风对氢气泄漏扩散影响

不利的环境风速致使氢气在障碍物附近积聚,为避免氢气局部积聚,在加注区罩棚下安装防爆风机,采用局部通风稀释加氢区内可燃区域。风机位于泄漏源与前侧加氢机中部,为正方形出风口(面积为1.44 m², Z=96.8~98 m),设置局部风速u为5、10 m/s,泄漏25 s后,启动加氢区罩棚下风机。不同局部通风风速作用下可燃区域分布如图14所示。当局部通风风速u为5 m/s时,风机开启1 s后,即t=26 s,可燃区域分布仅出现在前侧加氢机附近。当风机开启2 s后,即t=27 s,流体域内无可燃区域分布。局部通风风速u为10 m/s,流体域内可燃区域也在风机开启的2 s内消失。

为研究局部通风条件下,泄漏源附近气流运动情况。截取风机开启2 s时(t=27 s),泄漏孔横截面(X=17.78 m)处流线分布,如图15所示。当局部通风风速为5 m/s,排出气流向下作用于可燃区域,与氢气混合后撞击前方加氢机,改变扩散方向,向罩棚下侧运动,再沿罩棚壁面向外扩散。风机开启2 s时,仅有少量摩尔分数1%~2%氢气位于加氢机一侧,还有部分积聚在罩棚下侧,气体摩尔分数降低至1%左右。而局部通风风速增大至10 m/s,增强了出风口下方湍流强度,积聚在2加氢机之间的气体少于5 m/s的情况,氢摩尔分数在1%左右,罩棚下侧氢气摩尔分数分布情况相似。因此采用局部通风方式可增强泄漏源附近湍流强度,迅速稀释可燃区域摩尔分数。

由图14与图15可知:局部通风对降低泄漏源附近氢气摩尔分数作用较大。在相同稀释时间内,5 m/s的风机风速与10 m/s的局部通风风速作用下,都存在积聚氢气。为分析不同局部通风风速对加氢机附近氢气摩尔分数分布影响,沿泄漏源前方加氢机竖直方向布置一列监测点(Z=102 m),位于射流中轴线正前方,监测点布置高度范围为0.5~2.5 m,间隔为0.5 m,监测点位置如图16所示。t=25 s泄漏停止,并启动风机。当风速为5 m/s,监测点3—监测点5由于靠近泄漏源高度,且泄漏氢气在浮力作用下向上扩散,导致监测点3至监测点5摩尔分数高于下侧监测点1与监测点2摩尔分数。在t=25 s后,监测点摩尔分数均快速下降,在开启风机6 s后,即t=31 s,监测点3—监测点5氢气浓度低于加氢站传感器报警阈值(10%可燃下限)。当局部通风风速为10 m/s时,监测点浓度衰减速度快于5 m/s的情况,在t=28 s前,监测点2—监测点5的氢摩尔分数降至0.5%左右,t=31 s时,监测点2—监测点5摩尔分数均在加氢站传感器报警阈值以下。

3 结论

收起

1) 在无环境风情况下,5与10 mm泄漏孔径所形成的可燃区域超出加氢区范围,而2 mm泄漏孔径形成的可燃区域仅分布在泄漏源与前方加氢机之间。不同泄漏孔径下,罩棚下侧氢摩尔分数分布趋势不同:泄漏孔径为2 mm时,不会在罩棚下表面形成可燃区域;当泄漏孔径增大至5与10 mm时,垂直泄漏方向下,罩棚下侧表面摩尔分数最高点位于射流中轴线平行处附近;当泄漏孔径为10 mm时,罩棚下监测曲线摩尔分数呈现高斯分布。

2) 当泄漏孔径为2 mm时,加入垂直泄漏方向环境风速,即沿X轴正方向,环境风速为2与8 m/s可有效减少泄漏源附近氢摩尔分数;而风速为5 m/s,在前方加氢机附近产生回流区,诱导了氢气积聚。

3) 无环境风条件下,在加氢区罩棚下安装风机,可加速泄漏源附近氢气摩尔分数稀释。局部通风风速为5与10 m/s时,风机开启的2 s,流场内氢气云摩尔分数均降至可燃范围以下。在相同时间内,风机风速为10 m/s时,对流场内氢气云摩尔分数稀释效果更佳,前侧加氢机附近监测点摩尔分数衰减速度更快。

基金

收起

北京市科技计划“揭榜挂帅”项目(Z231100003823020)
北京市高层次创新创业人才支持计划科技新星计划项目(20240484576)
中央高校基本科研业务费资助(2024ZKPYAQ05)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

Huang

Ying

, Zhou

, Zhong

Ruohan

, et al. Hydrogen energy development in China: potential assessment and policy implications[J]. International Journal of Hydrogen Energy, 2024, 49: DOI: 10.1016/j.ijhydene.2023.10.176.

[2]

Yao

Yongzheng

, Tan

Luyao

, Chen

Fang

, et al. Hydrogen energy industry in China: The current status, safety problems, and pathways for future safe and healthy development[J]. Safety Science, 2025, 186: DOI: 10.1016/j.ssci.2025.106808.

[3]

彭善碧, 罗雪. 埋地掺氢天然气管道泄漏扩散数值模拟研究[J]. 中国安全科学学报, 2024, 34(3): 63-69.

Peng

Shanbi

, Luo

Xue

. Numerical simulation of leakage and difusion in buried hydrogen-blended natural gas pipeline[J]. China Safety Science Journal, 2024, 34(3): 63-69.

[4]

2024中国氢能产业发展半年报[EB/OL]. (2024-07-05). https://pc.h2cn.org.cn/3vdja5/202407/e38c086f01d3ebf90e2602ea09e95c87.htm. https://pc.h2cn.org.cn/3vdja5/202407/e38c086f01d3ebf90e2602ea09e95c87.htm

[5]

戴君豪, 杨石刚, 方秦, 等. 储氢温度对加氢站泄漏爆炸事故的影响研究[J]. 防灾减灾工程学报, 2024, 44(2): 293-300,371.

Dai

Junhao

, Yang

Shigang

, Fang

Qin

, et al. Study on the influence of hydrogen storage temperature on leakage and explosion incidents at hydrogen refueling stations[J]. Journal of Disaster Prevention and Mitigation Engineering, 2024, 44(2): 293-300.

[6]

Qian

Jinyuan

, Li

Xiaojuan

, Gao

Zhixin

, et al. A numerical study of hydrogen leakage and diffusion in a hydrogen refueling station[J]. International Journal of Hydrogen Energy, 2020, 45(28): 14 428-14 439.

[7]

Xiaojuan

, Xu

Yixiang

, Qian

Jinyuan

, et al. Effect of wind condition on unintended hydrogen release in a hydrogen refueling station[J]. International Journal of Hydrogen Energy, 2021, 46(7): 5 537-5 547.

[8]

Liang

Yang

, Pan

Xiangmin

, Zhang

Cunman

, et al. The simulation and analysis of leakage and explosion at a renewable hydrogen refuelling station[J]. International Journal of Hydrogen Energy, 2019, 44(40): 22 608-22 619.

[9]

Swain

M R

, Grilliot

E S

, Swain

M N

. The application of a hydrogen risk assessment method to vented spaces[M]. Boston: Springer, 2002: 163-173.

[10]

田苗, 王伟, 晏明海, 等. 室外立体停车场氢泄漏扩散特性及危险性分析[J]. 中国安全科学学报, 2025, 35(7): 176-183.

Tian

Miao

, Wang

Wei

, Yan

Minghai

, et al. Analysis on hydrogen leakage diffusion characteristics and risks in outdoor multi-story park[J]. China Safety Science Journal, 2025, 35(7): 176-183.

[11]

牛奕, 刘静雨, 房韬. 加氢站氢气泄漏扩散特性的小尺寸实验与模拟研究[J]. 中国安全生产科学技术, 2024, 20(5): 130-137.

Niu

, Liu

Jingyu

, Fang

Tao

. Small-scale experiment and simulation study of hydrogen leakage and diffusion characteristics in hydrogen fueling station[J]. Journal of Safety Science and Technology, 2024, 20(5): 130-137.

[12]

Tamura

, Takeuchi

, Sato

. Effectiveness of a blower in reducing the hazard of hydrogen leaking from a hydrogen-fueled vehicle[J]. International Journal of Hydrogen Energy, 2014, 39(35): 20 339-20 349.

[13]

Xie

Hong

, Li

Xuefang

, Christopher

D M

. Emergency blower ventilation to disperse hydrogen leaking from a hydrogen-fueled vehicle[J]. International Journal of Hydrogen Energy, 2015, 40(25): 8 230-8 238.

[14]

Cui

Weiyi

, Yuan

Yupeng

, Tong

Liang

, et al. Numerical simulation of hydrogen leakage diffusion in seaport hydrogen refueling station[J]. International Journal of Hydrogen Energy, 2023, 48(63): 24 521-24 535.

[15]

Pan

Xuhai

, Xue

Zhenming

, Wang

Shiqi

, et al. Numerical investigation on the effects of smoke barriers on hydrogen dispersion in a full-scale underground garage[J]. International Journal of Hydrogen Energy, 2024, 73: DOI: 10.1016/j.ijhydene.2024.06.012.

[16]

Pitts

M W

, Yang

J C

, Blais

, et al. Dispersion and burning behavior of hydrogen released in a full-scale residential garage in the presence and absence of conventional automobiles-supplemental video materials[J]. International Journal of Hydrogen Energy, 2012, 37(22): 17 457-17 469.

[17]

黄腾. 大型地下车库内氢气泄漏与扩散研究[D]. 济南: 山东大学, 2022.

Huang

Teng

. Study of hydrogen leakage and dispersion in a large underground parking garage[D]. Jinan: Shandong University, 2022.

[18]

Birch

A D

, Hughes

D J

, Swaffield

. Velocity decay of high pressure jets[J]. Combustion Science and Technology, 1987, 52: DOI: 10.1080/00102208708952575.

2026年第36卷第5期

PDF下载

175

引用本文

BibTeX

文章信息

doi: 10.16265/j.cnki.issn1003-3033.2026.05.0886

接收时间：2025-12-10
首发时间：2026-06-29
出版时间：2026-05-28

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-12-10
修回日期：2026-02-11

基金

北京市科技计划“揭榜挂帅”项目(Z231100003823020)

北京市高层次创新创业人才支持计划科技新星计划项目(20240484576)

中央高校基本科研业务费资助(2024ZKPYAQ05)

作者信息

中国矿业大学(北京) 应急管理与安全工程学院, 北京 100083

通讯作者:

^** 姚勇征(1991—),男,河北承德人,博士,副教授,主要从事新能源风险防控研究。E-mail: yaoyz@cumtb.edu.cn。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/zgaqkxxb/CN/10.16265/j.cnki.issn1003-3033.2026.05.0886

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT