关于《风速对综合管廊天然气管舱泄漏扩散影响的数值模拟》研究的学术报告
本报告旨在向学术界同仁介绍周宁、任福平、陈兵、李雪、袁雄军、赵会军等研究者于《天然气化工—C1化学与化工》2020年第45卷上发表的一项数值模拟研究。该研究主要来自常州大学石油工程学院（周宁、任福平、李雪、袁雄军、赵会军）以及中国安全生产科学研究院（陈兵）。这项工作的核心目标是探究通风风速对地下综合管廊中天然气管舱内发生气体泄漏时扩散行为的影响，其研究成果对于保障城市地下生命线工程的安全具有重要的理论与应用价值。
从学术背景来看，该研究属于公共安全工程、流体力学及计算传热学（Computational Fluid Dynamics, CFD）的交叉领域。随着中国城镇化进程的加速，地下综合管廊作为集约化敷设市政管线的工程形式日益普及。然而，管廊内部，尤其是容纳高压天然气的独立舱室（即“管舱”），一旦发生管道泄漏，极易形成爆炸性气体环境，构成重大安全风险。尽管已有研究关注开放空间或埋地管道的燃气泄漏，但针对地下综合管廊这一特定受限空间，其内部的泄漏扩散规律，特别是主动通风（换气）条件的影响，尚不明确。本研究正是在此背景下展开，旨在通过数值模拟方法，系统揭示不同通风风速下管舱内天然气（主要成分为甲烷）泄漏后的扩散特性、浓度分布演化规律及可燃气体监测报警响应，从而为管廊通风系统的优化设计、泄漏应急响应策略的制定提供科学依据。
该研究的详细工作流程严谨而系统，主要分为以下几个关键步骤： 首先，建立物理与数学模型。研究者以实际工程规范（GB 50838-2015）为依据，构建了一个长200米、断面尺寸为2.2米×2.6米的典型天然气管舱三维几何模型。泄漏口设置在管道中心位置，舱体一端设进风口，另一端设排风口。研究选择了两种典型的泄漏场景：小孔泄漏（孔径15mm）和大孔泄漏（孔径45mm）。在数学模型方面，研究采用了基于雷诺平均纳维-斯托克斯方程（RANS）的Realizable k-ε湍流模型来模拟舱内复杂的湍流流动，同时耦合组分输运模型（Species Transport Model）来追踪甲烷气体在空气中的扩散与混合过程。控制方程包括连续性方程、动量方程、能量方程及组分输运方程，通过求解这些方程来获得流场与浓度场的时空分布。
其次，设定模拟工况与边界条件。研究核心变量为通风风速。研究者依据管廊通风换气次数的常规设计值（最大12次/小时），通过换算公式确定了对应的入口风速范围，并设置了从0 m/s（无风）到3.81 m/s的一系列风速工况进行对比分析。泄漏源被设定为持续泄漏，模拟了从泄漏初始到相对稳定状态的动态过程。
第三，执行数值计算与数据提取。利用商业CFD软件（如Fluent）对上述模型和工况进行求解计算。模拟重点关注了几个关键物理过程：1）甲烷云团的时空演化形态；2）舱内速度场、湍流场的结构，特别是涡旋的产生与发展；3）沿管舱长度方向布置的多个可燃气体探测器（模拟中为监测点）处甲烷浓度随时间的变化，以此确定达到爆炸下限20%（即报警浓度）所需的时间（报警时间）。
第四，进行系统的结果分析与机理探讨。研究对海量的模拟输出数据进行了深入分析，不仅描述了现象，更致力于揭示现象背后的流体力学机理。例如，通过分析流线图和涡量场，解释了无风状态下对称涡对的生成、舱顶壁面对射流的反射作用，以及有风状态下下风向区域大涡破碎、湍流增强等物理机制。同时，将浓度场分布与报警时间数据相关联，量化了风速对泄漏危险区域范围及监测系统响应速度的影响。
本研究取得了若干重要而细致的发现，结果部分逻辑严密，环环相扣： 在无风（风速为0 m/s）工况下，模拟清晰地显示，泄漏初期甲烷气体呈“蘑菇云”状对称扩散。由于缺乏横向气流的驱散作用，气体主要受初始射流动量、浮力及与舱顶壁面相互作用的影响。小孔泄漏由于泄漏动量相对较小，气体在浮力作用下向上运移，在舱顶受阻后向两侧蔓延，导致甲烷浓度在竖直方向上出现明显的分层现象。相比之下，大孔泄漏射流动量巨大，气体能迅速冲击舱顶并沿顶壁向更远距离扩散，分层现象不如小孔泄漏明显。流场分析表明，泄漏口附近因速度梯度产生初始涡对，随后在舱顶反射区形成回流区，这些涡结构的演变显著影响着气体的卷吸与混合过程。从安全监测角度看，报警时间沿管舱长度呈“V”型分布，即泄漏点附近报警最快，向两侧逐渐延迟。大孔泄漏的初始报警时间（0.7秒）远快于小孔泄漏（1.1秒），凸显了泄漏孔径对早期预警的关键影响。
在引入通风（有风）工况后，气体扩散模式发生了根本性转变。天然气云团不再对称分布，而是整体向下风向“漂移”，其形态被描绘为“蜗牛”状或“水母”状。上风向区域的甲烷被新鲜空气迅速稀释，浓度显著降低；而下风向区域则因气流的输送和湍流混合的增强，成为气体聚集的主要区域。流场显示，通风改变了涡旋结构，下风向区域的大尺度涡团失稳并破碎成更多小尺度涡旋，这种湍流强度的增大加剧了气体与空气的混合（卷吸作用），但也可能导致局部浓度起伏。最重要的安全相关发现是，通风风速对报警时间产生了规律性影响。随着风速增大，上风向探测器报警大幅延迟甚至不再报警（因为浓度始终低于阈值）；而下风向探测器的报警时间，则与探测器到泄漏点的距离呈现出近似线性的增长关系。研究指出了一个关键的安全阈值：当通风风速达到或超过3.81米/秒时，无论是小孔还是大孔泄漏，整个管舱内的甲烷浓度在任何位置都无法累积到爆炸下限的20%，这意味着在此风速下，常规的可燃气体探测器将不会触发报警。这一发现对通风系统的安全设计具有直接指导意义。
基于以上结果，本研究得出了明确结论：通风风速是控制管廊天然气管舱泄漏后果的关键工程参数。它不仅能改变泄漏气体的空间分布形态（从对称分布变为向下风向漂移），还能通过改变流场结构（涡旋演变、湍流强度）影响混合效率。从安全工程角度，研究定量地揭示了风速与监测报警响应之间的关系，并找到了一个可能使探测器“失效”的临界风速值。这警示我们，在管廊通风系统设计与运行中，并非风速越高越好，需要综合考虑稀释效果与监测可靠性，寻求最优通风策略。
本研究的亮点突出体现在以下几个方面：首先，研究对象的特殊性：聚焦于新兴且安全风险高的地下综合管廊天然气管舱这一特定受限空间，填补了该领域细致化研究的空白。其次，机理探究的深度：不仅呈现了浓度分布的宏观结果，更深入结合流场（涡旋、湍动能）分析，从流体力学本质上解释了不同风速下扩散模式差异的成因，体现了从现象描述到机理揭示的深化。最后，明确的工程指导价值：研究结论直接关联工程实践，特别是关于报警时间分布和临界通风风速的发现，为管廊通风标准的细化、泄漏监测探头的最佳布置以及应急预案的制定提供了关键的数据支持和理论依据。该工作将CFD数值模拟技术与公共安全需求紧密结合，是一项具有重要理论意义和实用价值的优秀研究。