架空天然气管道泄漏扩散过程特性分析.docxVIP

第49卷第4期 2012年8月化 工 设 备 与 管 道 Vol.49 第49卷第4期 2012年8月 PROCESS EQUIPMENT PIPING Aug.2012 架空天然气管道泄漏扩散过程特性分析 刘庆宝，王娇阳 (大庆油田有限责任公司天然气分公司，黑龙江大庆 163500) 摘 要：在天然气管道集输过程中，由于管道老化、腐蚀等因素引起的泄漏事故时有发生，会引起穿孔泄漏后的 天然气扩散，可能会引发火灾、中毒或爆炸。因此，进行架空天然气管道泄漏扩散过程的影响因素模拟研究及分 析，对天然气管道输送安全运营和保障人生财产安全意义重大。文章建立了架空天然气管道在迎风坡和背风坡处 发生泄漏扩散过程的物理模型，并利用CFD软件对架空天然气管道泄漏后的甲烷气体扩散过程进行了数值模拟。 结果表明：泄漏的甲烷气体扩散范围和速度非常大，在泄漏时间300 s 时泄漏距离已经达到40m, 并在40 m处形成 高浓度的甲烷气团；泄漏点在背风坡时，山坡附近形成高浓度甲烷气团。 关键词：天然气管道；泄漏；扩散 中图分类号： TQ058.81;TH 183 文献标识码：A 文章编号：1009-3281(2012)04-0061-03 管道运输具有安全密封，受地形影响小，能长 期安全稳定运行等优点，已成为石油工业中最重要的 运输工具[1]。天然气主要成分为甲烷，并含有腐蚀 性气体，其对管道内壁的化学腐蚀以及长期的外腐 蚀，使管壁失效，在管内气体压力的作用下，出现穿 孔泄漏。管道穿孔泄漏是天然气管道泄漏最常见的情 况之一 [2]。当空气中甲烷达到5%~10%时遇火会引 发爆炸。因此天然气管道发生泄漏时，有爆炸和中毒 的危险，从而对人民群众的生命和财产安全构成严重 危害。天然气泄漏规律的研究对于安全输送和环境保 护具有重要的现实意义。 目前针对天然气管道泄漏扩散过程，主要利用 半经验半理论计算公式对危险气体进行研究[3-10]。 这些公式具有计算快、编程简单等优点，其不足之处 是计算的准确性对气体扩散条件与模型假设的接近 程度过于依赖。随着计算流体力学、以及计算机技术 的发展，数值模拟研究方法开始应用于泄漏过程的研 究。利用CFD 软件对架空天然气管道在迎风坡和背 风坡处发生泄漏时的气体扩散规律进行研究，得出了 天然气泄漏气体扩散规律，为天然气的安全输送及环 境保护提供了理论依据。 1 泄漏扩散控制方程 由于管内气体压力高，在泄漏口处传输速度较 大，为湍流流动。 k- ε是在湍流方程的基础上引入一 个关于离散率的方程，是目前广泛应用的湍流模型， 计算时采用标准的k-ε 模型。 物质传输模型 的湍流脉动动能方程(K 方程): (1) G+G-pe-X+S 湍流动能耗散方程(ε方程): (2) 湍流粘度： (3) 其中： G 为平均速度梯度引起的湍动能k 的 产生项， G, 为浮力引起的湍动能k 的产生项， Ym 为可压缩湍流中脉动扩张贡献。模型中Ci=1.44, C2e=1.2,C=1,C=0.9 为经验常数。 2 数值模拟及结果分析 2.1 模拟工况 模拟区域(见图1)100m×40m, 土坡高度4m, 收稿日期：2011-12-30 作者简介：刘庆宝(1983一),男，山东省泗水人，工程师。从事 油气储运管理工作。 第49卷第4期·62 ·化 工 设 备 与 管 第49卷第4期 ·62 · 坡度45°,泄漏点(1)处直径为10 cm。 环境风速 为10m/s, 天然气泄漏速度为10 m/s。 泄漏气体为 甲烷。模拟泄漏点在迎风坡处和背风坡处发生泄漏 时的情况。 图1 模拟结构示意 2.2 边界条件 BC 边为墙壁， AD 边为自由出流边界，1点处为 速度入口。 工况1: AB 边为速度入口， CD 边为自由出流 边界。 工况2: CD 边为速度入口， AB 边为自由出流 边界。 2.3 模拟结果 2.3.1 工况1时的模拟结果 当泄漏点在迎风坡时，通过观察大气中的甲烷 分布图(见图2~4),泄漏的甲烷气体在泄漏100s 时 已经主要分布在背风坡处；在200s 时，泄漏的甲烷 气体开始聚集成气团，并以大气团扩散；在300s 时， 甲烷气团已经扩散到40m 处。 图2 100 s时大气中甲烷分布 图3 200 s时大气中甲烷分布 图4 300 s时大气中甲烷分布 2.3.2 工况2时的模拟结果 当泄漏点在背风坡时，通过观察大气中的甲烷 分布图(见图5~7),泄漏的甲烷气体在泄漏100 s 时已经分布