摘　要：采用计算流体力学(CFD)方法，在天然气持续泄漏条件下，数值模拟风速对障碍物迎风面天然气体积分数分布的影响，分析庭院燃气管道泄漏对迎风面住户的危害。

关键词：天然气泄漏  计算流体力学  庭院管道  迎风面住户

Abstract：The influence of wind speed on natural gas volume fraction distribution at windward side of obstacles is numerically simulated by computational fluid dynamics(CFD)method under the condition of natural gas continuous leakage．The hazard of courtyard gas pipeline leakage to households at windward side is analyzed．

Keywords：natural gas leakage；computational fluid dynamics；courtyard pipeline；households at windward side

 

进入20世纪后期，我国能源结构开始发生变化，天然气需求量迅速增加。天然气具有热值高、来源广、环境友好等优点，但天然气泄漏后易扩散，易导致人员中毒甚至发生爆炸，对人民的生命财产造成危害，因此研究燃气泄漏后的扩散规律具有重要意义^[1]。

现代计算机技术的迅速发展使得采用计算流体动力学(Computational Fluid Dynamics，CFD)方法研究各种流体的流动、传热、传质、燃烧等过程成为可能^[2]。与传统的实验和理论分析方法相比，CFD方法具有成本低且条件设置简便等优点。近年来，国内外学者开始采用CFD方法对市政管道、储罐、室内管道燃气泄漏问题进行研究^[3-4]。

目前，也有很多关于气体泄漏和扩散动态模拟的研究成果发表：吴晋湘等人^[5]利用CFD方法对燃气持续泄漏后的扩散现象进行了数值模拟。刘欣^[2]利用CFD方法对海上油气平台可燃气体泄漏进行了模拟计算。肖淑衡^[6]利用泄漏危险性气体在大气风流中运动的数值计算模型，对天然气等轻质气体泄漏扩散运动进行了数值模拟。X．P．Liu等人^[7]对高层建筑周围燃气泄漏受到不同风速的影响进行了分析。

本文采用CFD方法，研究庭院天然气管道在持续泄漏条件下，不同风速对障碍物迎风面天然气体积分数分布的影响。

研究区域的物理模拟见图1，模拟庭院燃气管道持续泄漏对迎风面住户的危害。图中红色矩形体为障碍物，模拟1～3层住户的迎风面，蓝色圆点为天然气泄漏口，图中数值单位为m。坐标系原点O位于底面中心点。障碍物底面中心点位于(0，20，0)，长度(x轴方向)为20m，厚度(y轴方向)为1m，高度(z轴方向)为10m。泄漏口圆心位于(0，25，0)，直径为20mm，燃气流向向上，空气来流方向指向y轴负方向。模型底面为地面，地面设置为固体边界，模型其他表面为自由出口压力边界，设为常压。燃气泄漏压力稳定在2000Pa，燃气质量流量为60.58kg／h，采用稳态数值模拟，以甲烷代替天然气，风速分别选取l、5、10m／s。

 

①泄漏控制方程

天然气泄漏扩散可认为是多组分气体相互作用的湍流，泄漏的天然气与空气混合形成爆炸性混合气体，混合气体各组分之间不发生化学反应，泄漏控制方程满足连续性方程、动量方程、组分方程^[8]。

②流动方程

选择湍流计算中应用范围最广的k-e湍流方程，通过增加单位质量流体湍流脉动动能后和单位质量流体脉动动能耗散率占的控制方程使湍流时均方程组封闭，标准方程形式见文献[9]。

由于甲烷爆炸极限(体积分数)为5％～l5％，因此本文给出的障碍物迎风面甲烷体积分数分布云图可辨识的甲烷体积分数(称为甲烷辨识体积分数)范围设定为0～15％，并将甲烷爆炸极限定义为危险体积分数范围。

①yOz平面的体积分数分布

不同风速下yOz平面上甲烷辨识体积分数分布云图(视线方向为x轴负方向)见图2～4。由图2可知，当风速为1m／s时，障碍物迎风面危险体积分数范围分布在7～10m高度范围内。当风速为5m／s时，分布在0～5m高度范围内。当风速为l0m／s时，分布在0～3m高度范围内。

 

 

②不同高度上的体积分数分布