海上掺氢天然气管道中弯头的混合机理与几何效应

《Journal of Marine Science and Engineering》：Mixing Mechanism and Geometric Effects of Elbows in Offshore Hydrogen-Blended Natural Gas Pipelines

　　

利用现有海底天然气管道进行掺氢，为海上低碳能源运输提供了一条实用途径。然而，传统T型管需要较长距离才能实现均匀混合条件，增加了分层和氢脆风险。本研究创新性地将常规管件——弯头用作掺氢天然气管道的被动混合元件，期望利用曲率诱导的二次流和涡流重组来加速均质化，无需额外的流动混合器。研究人员进行了数值模拟，以研究氢支管与弯头入口之间的距离L₁以及弯头曲率半径R_c的影响。分析了氢分布、变异系数（coefficient of variation, COV）、均匀混合路径和流动涡量，以比较弯头构型与传统T型管。结果表明，较小的L₁缩短均匀混合路径S_h，而较小的弯头曲率半径R_c产生更强的二次流，但也加剧了由离心力引起的非对称氢富集。当掺氢比（hydrogen blending ratio, HBR）≥20%时，弯头构型的混合效果更为显著，并随HBR降低而逐渐减弱。本研究中，在L₁=0和R_c/D=4的最优结构下获得了最短混合路径S_h=11.73 m，相对于传统T型管减少了67.3%。这证明了弯头结构能显著提高掺氢天然气管道的混合效率，为海上管道掺氢提供了新的解决方案。

氢能作为连接传统化石能源与可再生能源的重要桥梁，其清洁零碳特性正推动其在海洋工程等领域的应用。随着基于海上风电的“绿氢”等可再生氢生产技术发展，利用日趋成熟的海底天然气管道基础设施注入氢气并形成混合气共同输送，展现出显著的应用潜力。然而，氢气与天然气之间较大的密度差异会导致海底管道内混合不均甚至分层，这不仅可能造成局部氢浓度过高而诱发钢管氢脆，加速材料老化开裂，还会影响末端供气稳定性与系统安全。传统T型管需要较长距离才能达到均匀混合，增加了上述风险。因此，针对海底管道特点开发能在短距离内实现快速均匀混合的高效结构，具有重要的工程意义和紧迫的现实需求。

本研究创新性地将常规管件——弯头用作掺氢天然气管道的被动混合元件，期望利用曲率诱导的二次流与涡流重组加速均质化，而无需附加流动混合器。已有研究多聚焦于T型注入段，而弯头结构在氢-天然气浓度均质化方面的作用尚未被充分探索。研究人员采用数值模拟方法，考察了氢支管与弯头入口距离L₁和弯头曲率半径R_c的影响，以传统T型管为基准，比较了氢分布、变异系数（COV）、均匀混合路径和流动涡量。研究表明，较小的L₁可缩短均匀混合路径；较小的R_c能产生更强的二次流，但也会加剧离心力引起的非对称氢富集。弯头结构的混合增强效果在掺氢比（HBR）≥20%时更为显著，并随HBR降低而减弱。最优结构（L₁=0、R_c/D=4）获得最短混合路径为11.73 m，相对传统T型管缩短67.3%，证明弯头结构能显著提升掺氢天然气管道混合效率，为海上管道掺氢提供了新方案。该成果发表于《Journal of Marine Science and Engineering》。

在技术方法上，研究人员使用ANSYS FLUENT 17开展三维数值模拟。物理模型尺寸取自中国电力投资集团朝阳天然气掺氢示范项目管道网络：主管直径D=0.8 m，氢支管直径0.32 m，上游总长度20 m，出口压力1.8 MPa，甲烷入口速度5 m/s，氢气掺混比25%（体积分数），温度300 K。采用标准k-ε湍流模型求解连续性、动量、能量及组分输运方程；氢气在甲烷中的分子扩散系数设为2.88×10^-5 m²/s，湍流施密特数Sct=0.7。采用SIMPLE算法与二阶迎风离散。网格独立性验证显示Mesh 3与Mesh 4的COV平均相对差异为1.37%，选用Mesh 3。模型经Fernandes等人压降结果和Gardel关联式验证。额外的RSM模拟比较表明，湍流模型选择对相对混合增强影响有限。

关于支管与弯头入口距离的影响：在固定R_c/D=3时，改变L₁=0、5、10、15 m。流场显示，L₁=0 m时直管段流线严重缠结，弯头区形成主涡对及额外涡对，下游出口处出现瞬态四涡结构，约3 m内衰减；而L₁=15 m时下游保持较稳定的双涡结构。COV衰减曲线表明，L₁=0 m初始下降最快。截面平均涡量在L₁=0 m时高达60 s^-1，随L₁增大而减小，且沿弯头转角差异逐渐收敛。均匀混合路径S_h分别为15.53 m、26.53 m、37.23 m和26.83 m，呈非单调变化。结论是，短L₁使注入氢在流动充分发展前进入弯头，强化与二次流的相互作用，促进横向输运。

关于曲率半径的影响：固定L₁=0 m，设置R_c/D=2、3、4、5。结果显示，R_c/D=2时氢羽流强烈偏转，产生显著非对称富集；随R_c/D增大，非对称性减弱，R_c/D=5时高速区移向内弧面，浮力驱动的向上迁移贡献增大。下游流场演化表明，R_c/D=2和3的外涡约4 m后消失，R_c/D=4在3 m后形成稳定双涡，R_c/D=5则经历三涡、双涡和单涡转变。弯头内平均COV降低率随R_c/D增大而减小，截面平均涡量在R_c/D=2最高、R_c/D=5最低。均匀混合路径随R_c/D增大先减后增，分别为21.41 m、15.53 m、11.73 m、13.38 m，在R_c/D=4时最短。该结果表明，更强的二次流并不必然带来更短的混合路径，需平衡横向输运与浓度均匀性。

关于与传统T型管的比较：在相同中心线路径长度约1.6D处，弯头构型（L₁=0、R_c/D=3和4）比传统T型管呈现更宽的氢分布、更低的局部峰值浓度，并产生额外涡对形成四涡结构，而T型管主要为双涡。COV演化显示弯头进入弯曲段后下降更快。相对传统T型管，均匀混合路径分别缩短56.7%和67.3%。在不同HBR下，15%时缩短12.55%，20%时44.75%，25%时67.3%，表明弯头的相对混合优势随HBR增加而增强，尤其HBR≥20%时更明显。

讨论与结论部分指出，弯头作为现有管道组件，无需附加扰动装置即可强化氢-天然气混合。研究结论可概括为：第一，L₁对均匀混合路径的影响非单调，但将氢注入支管靠近下游弯头在研究的条件下有利于混合，管道布局设计应重视注入点与弯头的相对位置。第二，较小的R_c/D产生更强二次流，但小曲率带来的强离心效应加剧非对称富氢，因此存在最优曲率半径（本研究中R_c/D=4）以平衡输运与均匀性。第三，弯头结构在15%～25% HBR范围内均保持较短的均匀混合路径，且HBR越高优势越显著。需要指出，这些结论适用于所研究的运行条件，R_c/D=4并非普遍最优；未来需考虑更多掺氢比、流动条件和管道构型以建立通用设计准则。对于运行条件相当的海上管段，将氢注入支管靠近既有弯头提供了一种无需附加混合装置的被动方案，可有效缩短均匀混合路径。