混氢天然气管道放空自燃过程数值模拟分析.pdf
《混氢天然气管道放空自燃过程数值模拟分析.pdf》由会员分享,可在线阅读,更多相关《混氢天然气管道放空自燃过程数值模拟分析.pdf(13页珍藏版)》请在咨信网上搜索。
1、第 43 卷第 9 期2023 年 9 月 149 天然气工业Natural Gas Industry引文:朱红钧,李佳男,陈俊文,等.混氢天然气管道放空自燃过程数值模拟分析J.天然气工业,2023,43(9):149-161.ZHU Hongjun,LI Jianan,CHEN Junwen,et al.Numerical simulation analysis of the spontaneous combustion process during the venting of hydrogen-mixed natural gas line pipesJ.Natural Gas Indus
2、try,2023,43(9):149-161.混氢天然气管道放空自燃过程数值模拟分析朱红钧1李佳男1陈俊文2粟华忠1唐 堂11.西南石油大学石油与天然气工程学院2.中国石油工程建设有限公司西南分公司摘要:天然气掺氢输送是实现大规模、远距离氢能转运的主要手段之一,但混氢天然气在放空过程中存在爆轰或爆燃的风险,对放空管壁施加的超压过大还会造成管壁破裂。因此,亟需明晰混氢天然气放空自燃与流场演化过程,进而量化放空管壁的一次超压。为此,利用计算流体力学方法数值模拟了混氢天然气管道放空自燃过程,对比分析了不同掺氢浓度条件下对自燃及压力波传播的影响规律,揭示了混氢天然气在阀门通道和放空管中的爆燃机制。研究
3、结果表明:高压气体在阀门通道内以压力波形式传播并不断地碰撞反射与叠加,形成马赫环结构,加热气体使温度升至自燃点,并触发自燃;在阀门通道和放空管内均出现了爆燃现象,但压力波在放空管内能量迅速衰减,一段距离后温度大幅降低,气体不再燃烧;掺氢比越大,压力波传递速度越快,自燃触发的时间越短,对放空管壁产生的一次超压也越大。结论认为,在确定的泄放压力与阀门开度工况下,降低掺氢浓度可有效减轻爆燃风险和减小壁面超压,实际工程中需结合输送经济成本和安全风险控制掺氢量在合理的浓度范围。混氢天然气管道放空自燃过程数值模拟分析取得的新认识有助于指导掺氢天然气管道的安全运行,将助力于绿色氢能的大规模混合输运。关键词:
4、混氢天然气;放空;爆燃;超压;自燃;数值模拟DOI:10.3787/j.issn.1000-0976.2023.09.015Numerical simulation analysis of the spontaneous combustion process during the venting of hydrogen-mixed natural gas line pipesZHU Hongjun1,LI Jianan1,CHEN Junwen2,SU Huazhong1,TANG Tang1(1.Petroleum Engineering School,Southwest Petroleum
5、 University,Chengdu,Sichuan 610500,China;2.Southwest Company,China Petroleum Engineering&Construction Corporation,Chengdu,Sichuan 610041,China)Natural Gas Industry,vol.43,No.9,p.149-161,9/25/2023.(ISSN 1000-0976;In Chinese)Abstract:Hydrogen-mixed natural gas transportation is one of the important wa
6、ys for large-scale and long-distance hydrogen energy transportation.In the venting process of hydrogen-mixed natural gas,however,there are risks of detonation or deflagration,which exert excessive overpressure on the pipe wall,leading to its fracture.Therefore,it is in an urgent need to clarify the
7、spontaneous combustion and flow field evolution process during the venting of hydrogen-mixed natural gas,so as to quantify the primary overpressure on the wall of the venting pipe.In this paper,the spontaneous combustion process during the venting of hydrogen-mixed natural gas pipeline is numericall
8、y simulated by using the computational fluid dynamics(CFD)method,the influence laws of different hydrogen concentrations on spontaneous combustion and pressure wave propagation are comparatively analyzed,and the deflagration mechanisms of hydrogen-mixed natural gas in valve channels and venting pipe
9、s are revealed.And the following research results are obtained.First,high-pressure gas is propagated in the form of pressure wave in valve channels.During the propagation,it undergoes continuous collision,reflection and superimposition to form Mach ring structure,which heats the gas to the spontaneo
10、us combustion temperature,resulting in spontaneous combustion.Second,the phenomenon of deflagration appears in valve channels and venting pipes.In venting pipes,however,the energy of pressure wave decays rapidly,resulting in a significant reduction of temperature and hence the disappearance of gas c
11、ombustion.Third,the larger the hydrogen concentration is,the faster the pressure wave propagates,the shorter time it takes to trigger the spontaneous combustion,and the greater the primary overpressure on the venting pipe wall.In conclusion,with given discharge pressure and valve opening,reducing th
12、e hydrogen concentration can effectively reduce the deflagration risks and the overpressure on the pipe wall.In actual engineering,the hydrogen mixing volume shall be controlled within a reasonable concentration range based on transportation economic costs and safety risks.These new understandings p
13、rovide useful guidance for the safe venting of hydrogen-mixed natural gas pipelines and the large-scale mixed transportation of green hydrogen energy.Keywords:Hydrogen-mixed natural gas;Venting;Deflagration;Overpressure;Spontaneous combustion;Numerical simulation基金项目:四川省杰出青年科学基金项目“混氢输送管道非线性耦合动力学行为与流
14、固耦合机制研究”(编号:2023NSFSC1953)、中国石油集团工程股份有限公司支撑课题“含氢输气管道系统关键技术研究”(编号:2019ZYGC-01-09)、中国石油集团工程股份有限公司支撑课题“氢/氨储运及储存设备关键技术研究”(编号:2021ZYGC-03-01)。作者简介:朱红钧,1983 年生,教授,博士研究生导师;主要从事多场多相耦合力学方面的教学与研究工作。地址:(610500)四川省成都市新都区新都大道 8 号。ORCID:0000-0001-5977-5028。E-mail:2023 年第 43 卷 150 天 然 气 工 业0引言氢能具有较高的单位质量能量密度、零碳排放量
15、、转换效率优异等显著优点,大力开发氢能是推动中国能源绿色低碳转型和实现“双碳”目标的重要战略途径之一1-3。氢能产业的大规模发展与其储运技术的突破密不可分,利用已有天然气管网进行掺氢输送是实现大规模、远距离、低成本氢能转运的重要手段,英国、荷兰、德国、美国等已开展了天然气掺氢输送的中试实验4-5,但还局限于低压配气管线或加氢站站前管线。中国天然气掺氢输送方面的研究刚起步,国家电力投资集团有限公司“朝阳可再生能源掺氢示范项目”是国内天然气掺氢的首次尝试。现有天然气掺氢输送的研究主要关注氢气掺混的均匀度、管材的选择及泄漏风险等方面,重点评价局部管段可能存在的氢致失效(氢脆)风险6-7。然而,管道放
16、空系统是输气站场安全设施的重要组成部分,放空气体膨胀波引起温度和压力的急剧变化,过大的超压易造成管道破裂,引发严重后果8-11。掺氢天然气在节流放空的过程中存在正 JouleThomson 效应和逆 JouleThomson 效应的竞争,使得混合气体的放空温度、压力变化更为剧烈,并存在爆轰与爆燃的风险,给安全生产带来巨大隐患。现有标准对爆轰与爆燃的定义仍存争议,多数标准(如氢气放空系统标准)根据燃烧波是否超过音速来区分爆燃与爆轰,袁生学等12和归明月等13认为爆轰与爆燃的区分不仅要关注燃烧波速度,还需考虑火焰的结构。因此,明确掺氢天然气管道的放空泄压过程及管内压力、温度、速度等瞬变规律,是制订
17、安全预防和应急抢险措施的依据。高压氢气在放空过程中,可以在没有任何点火源的情况下出现自燃和爆炸现象14。Dryer 等15实验观测了高压氢气的放空过程,发现在常用氢气储存和管道压力范围内,氢气经射流释放会引发自燃,与激波和边界层的相互作用有关。Bragin 等16利用Fluent 软件模拟了高压储存的氢气泄放过程,发现自燃由壁面边界层开始。Kim 等17利用纹影系统和高速摄像观测了透明管道内气体的自燃过程,实验结果表明自燃首先发生在混合区前缘,并随着冲击波向下游移动,这与 Bragin 的数值结果一致。Xu 等18数值模拟了不同放空管结构对加压氢气释放自燃的影响,发现不同类型管道的自燃过程主要
18、与入射激波和反射激波有关。Lee 等19对高压氢气管道放空过程进行了数值模拟,结果表明如果掺氢气体放空过程中产生强烈的膨胀波且空气被充分加热,管内气体自燃则主要取决于是否形成了氢气和空气的混合区。与纯氢气管道放空方面的研究相比,关于掺氢天然气管道放空的报道较少。Rudy 等20实验研究了在高压氢气中添加 5%和 10%甲烷的射流自燃过程,结果表明,与纯氢气自燃压力相比,氢气甲烷混合物自燃所需压力更大。Golovastov 等21实验研究了不同压力、不同甲烷浓度对氢气射流自燃的影响,指出气体自燃所需的时间随压力的升高而减少,随掺入甲烷浓度的增加而增加。Zeng 等22使用透明矩形管道实测了高压氢
19、气添加 2.5%甲烷对自燃和火焰传播的影响,结果发现添加甲烷增加了初始点火时间,减弱了火焰强度,降低了火焰传播速度。上述研究大多关注于射流燃烧过程,对于掺氢天然气管道在放空过程中是否出现爆轰或爆燃过程还不明确,掺氢浓度对自燃及压力波传播的影响还不清楚。因此,有针对性地开展不同掺氢浓度混氢管道的放空数值模拟,明晰是否会发生爆轰或爆燃以及它们的流场演变过程,量化放空管壁的一次超压,有助于指导掺氢天然气管道的安全放空,助力绿色氢能的大规模混合输运。1物理模型选用如图 1 所示的简化管道模型,设置阀门通道与下游放空管初始时刻被空气填充(模拟极端未置换工况,以评价爆轰与爆燃的最大风险),放空管内空气的初
20、始温度为 300 K,压力为 0.1 MPa。上游管内气体为氢气与甲烷的混合物,压强为 8.0 MPa。改变上游管内氢气浓度(体积分数为 0、50%、100%)以分析其对混合气体放空过程的影响。上游管道入口定义为压强入口(8.0 MPa),下游放空管出口设定为压强出口(0.1 MPa)。图1管道简化模型图注:L1表示上游管道长度,mm;D1表示上游管道管径,mm;L2表示阀门通道长度,mm;D2表示阀门通道管径,mm;L3表示下游放空管道长度,mm;D3表示下游放空管道管径,mm;p 表示压强,MPa。第 9 期 151 朱红钧等:混氢天然气管道放空自燃过程数值模拟分析2数值方法2.1控制方程
21、混氢气体的放空过程为多组分、非稳态的可压缩气体湍流流动,用非定常雷诺时均纳维斯托克斯(RANS)方程进行描述:(1)(2)式中下标 i、j 表示坐标轴分量(i、j=1、2);表示混合气体密度,kg/m3;t 表示流动时间,s;表示混合气体的动力黏度,Pas;表示混合气体的时均流速,m/s;ij表示雷诺应力,Pa;表示气体压强,Pa;gi表示重力加速度在坐标轴 i 上的分量,m/s2。利用可实现 k(Realizable k)双方程湍流模型对 RANS 方程组进行封闭23:(3)(4)其中,式中 k 表示湍动能,m2/s2;表示湍动能耗散率,m2/s3;t表示涡黏系数;YM表示可压缩湍流中的膨胀
22、对总耗散速率的贡献;Gk、Gb表示平均速度和浮力产生的湍流动能,m2/s2;C1、C2、C3表示常数,分别取 1.44、1.90 和 0.0924;k、表示湍流普朗特数,分别取 1.0 和 1.2。考虑到可压缩气体的温度变化,建立能量守恒方程:(5)式中 T 表示气体温度,K;表示气体传热系数,W/(m2K);cp表示比热容,J/(kgK);ST表示化学反应产生的热量以及摩擦生成的热量之和,J。此外,各气体组分在流动过程中还满足组分质量守恒方程:(6)式中 cm表示组分 m 的质量分数;Dm表示组分 m 的扩散系数,m2/s;Sm表示组分 m 的化学反应速率,kg/(m3s)。选用单步反应的有
23、限速率模型25计算反应热,以保守估算爆轰、爆燃反应带来的风险。由涡耗散模型求解流场的燃烧反应速率,反应过程中物质 j 的生成速率(Rj)为:(7)(8)式中 vj表示反应物 j 的化学计量系数;Mw,j表示反应物 j 的分子量;A、B 表示经验常数,分别取 4.0、0.5;YR表示某一反应物组分的质量分数;vR表示反应物组分的化学计量系数;Mw,R表示反应物组分的分子量;YP表示某一生成物组分的质量分数;vP表示生成物组分的化学计量系数;Mw,P表示生成物组分的分子量。计算采用基于密度的求解器和显式算法,利用Coupled 算法耦合求解压力与速度,控制方程采用二阶迎风差分格式26离散,混合气体
24、的物性参数则根据气体混合定律计算。2.2网格划分及无关性验证为保证计算效率与精度,采用结构化网格划分计算域,并在阀门通道和管轴进行网格加密,图 2 显示了以管轴为对称轴的一半计算域网格。阀门通道壁面的第一层网格高度为 0.1 mm,第一层网格的无量纲高度 y+约为 30。网格的增长率通过网格无关性测试确定,如表 1 所示,G1 G5 为 5 套不同分辨率的网格,因纯氢气放空的速度、温度、压力变化都更大,故选用纯氢气进行网格测试。计算结果如表 1 所示,表中的相对误差体现了相邻两套网格 Gi图2模型网格图2023 年第 43 卷 152 天 然 气 工 业与 G(i+1)计算结果的差异。可见,G
25、4 和 G5 两套网格计算的结果相差小于 1.1%,综合考虑计算成本与结果的准确性,选用 G4 网格的分辨率进行后续计算。在该网格分辨率下,时间步长设为 210 8 s,以保证库朗数小于 0.5。表1网格无关性验证结果表网格 网格单元数 最大温度/K相对 误差最大马赫数(Mamax)相对 误差G1 57 7332 48020.85G2 70 2112 6707.67%21.452.88%G3 98 8762 7081.42%21.952.33%G4112 3522 7180.37%22.201.14%G5150 0002 7230.18%22.441.08%2.3模型验证为了验证数值模型的准确
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 天然气 管道 放空 自燃 过程 数值 模拟 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。