天然气泄漏扩散状态模拟及后果分析

王秀清 李书营 李建龙

(郑州铁路职业技术学院 河南 郑州 450052)

近年来，随着环保要求的提高及生活条件的改善，气体燃料，如人工煤气、液化气、天然气等逐渐成为城市居民生活燃料的主体，同时在工业领域的应用范围也在不断扩大。气体燃料的特点是易泄漏和易燃易爆，当这些可燃气体在室内有限空间内发生泄漏后，就会增加火灾和爆炸发生的危险性。同时，由于燃料的比重不同于空气，以及门窗等有限通风条件的存在，导致室内燃料浓度分布是一个非均匀的变化过程，而气体燃料的着火与爆炸的发生是有一定的浓度界线的。因此，系统地研究室内燃气泄漏的扩散过程与发展，爆炸发生的危险性等，对防止室内火灾的发生及控制具有重要的意义。

1 数值计算模型

1．1 几何模型

图1 调压室几何模型

本文的天然气泄漏点在坐标值为(10，8．5，0．5)处，调压室的窗户上配有轴流风机。

1．2 数学模型

空气与天然气在调压室进行混合，相应的基本控制方程如下［1］:

1．3 网格与计算求解

数值计算采用非结构化网格，采用ANSYSICEM CFD11．0划分网格，数值计算采用数值软件ANSYS CFX11．0。

燃气泄漏源给定泄漏速度、泄漏压力和泄漏量，燃气管壁和空气进口均处理为绝热边界条件。具体的边界条件及初始条件如下:泄漏源的温度

图2 计算网格全局图

T=273k，环境温度 T=288k，轴流风机速度:V=0．67m/s，标准大气压力 P0=1．013125MPa。其中泄漏量与泄漏速度的计算如下［2］:

式中:w为泄漏速度，ρ0为标准状态下燃气的密度，kg/m3，a 取0．64。

当在临界状态或超临界状态时，即［3］

按式(5)计算泄漏速度，当v处于超临界状态时，计算出的速度为超音速。由于孔口泄漏，无扩压段，最大泄漏速度只能达到临界状态的音速，所以当v达到超临界状态时，泄漏速度应按临界状态的速度计算［4］:

式中:ρ1=1．32kg/m3

通过公式(5)－(7)，可以算出当泄漏源压力为中压 P=0．4MPa，燃气泄漏速度 v=142．67m/s，泄漏质量流速Qm=0．00306kg/s;当泄漏源压力为高压P=1．6MPa，燃气泄漏速度 v=285．34m/s，泄漏质量流速Qm=0．00613kg/s。对于判断数值计算收敛残差的设置，取均方根残差RMS=104。

2 计算结果与比较

2．1 当屋顶开启量为60%

2．2 当屋顶开启量为80%

2．3 当屋顶开启量为100%

2．4 总结

从上述的计算结果图可以得出在不同屋盖开启量以及不同压力下天然气泄漏稳定后的最大浓度、最大温度及最大流速，下面的表格是对它们的具体比较:

水平方向 垂直方向极大质量分数(%)极大速度(m/s)极大质量分数(%)极大速度(m/s)100% 1．93 5．72 2．24 6．75 60% 1．92 2．86 2．24 3．37高压 60% 1．91 5．73 2．24 6．74中压 80% 1．91 2．86 2．24 3．37高压 80% 1．91 5．73 2．24 6．74中压 100% 1．92 2．86 2．24 3．37高压中压

3 结论

本文用CFX软件对地下调压间天然气泄漏进行了稳态模拟，得到了在不同压力下，不同屋盖开启量的情况下天然气的浓度场，流场分布，并对结果进行了详细地分析，得出以下几点结论:

1．天然气泄漏后，刚开始属于长瘦形，然后长度缩短，但高度却增加，变得饱满起来，而且变得上下不对称，这也是由于天然气比空气轻，泄漏后在浮力作用下，向上空扩散造成的。

2．在天然气沿水平方向泄漏和空气混合后的极大质量分数均为1．92%左右，沿垂直方向均为2．24%;而混合后沿水平方向的极大速度在中压时均为2．86m/s，高压时为5．73m/s;沿垂直方向的极大速度在中压时为 3．37m/s，高压为 6．74m/s。

3．泄漏以后调压间的温度非常均匀，最大温度均为288K，和环境温度一样基本没有变化。

4．随着屋盖开启量的增大，天然气扩散加快。

5．天然气的爆炸浓度范围在2．84% ～8．92%(质量百分比)，而通过模拟可得在以上各种情况下天然气泄漏达到稳定以后的最大浓度都低于爆炸浓度范围。由于轴流风机的抽气量为2．144m3/s，而泄漏源的最大泄漏量为0．0081m3/s，则轴流风机的抽气量为泄漏量的265倍，并且部分气体从楼梯窗户和屋顶开启部分扩散出去，因此泄漏的天然气很快就扩散出去，基本不会聚集在室内，使之达不到爆炸的浓度。

6．泄漏后不仅调压室内不会达到爆炸，而且备件间和监控室也不会达不到爆炸极限。

［1］王英．室内燃气泄漏扩散状态模拟及后果分析［D］．重庆大学硕士学位论文，2007．

［2］田贯三．管道燃气泄漏过程动态模拟的研究［J］．山东建筑工程学院院报，1999，14(4):57 －59．

［3］潘旭海，蒋军成．事故泄漏源模型研究与分析［J］．南京工业大学学报，2004，24(1):107 －109．

［4］于暢，田贯三．城市燃气泄漏强度计算模型的探讨［J］．山东建筑大学学报，2007，22(6):543 －545．