不同风速下不同高度(1、3、5、10m)平面上甲烷辨识体积分数分布云图分别见图5～7。由图5可知，当风速为1m／s时，在z＝1m平面上，近壁面甲烷辨识体积分数聚积宽度仅为2m，但随着高度的增大甲烷辨识体积分数分布范围向两侧扩大。在z＝5m平面上，聚积宽度达到l2m。在z＝10m处仍有较大聚积宽度，甚至能越过障碍物扩散到其背面。

 

 

 

由图6可知，当风速增大至5m／s时，在z＝1m平面上，障碍物近壁面甲烷辨识体积分数聚积宽度达到l0m，危险体积分数范围分布高度比风速为1m／s时有所降低，约在5m以下。虽然在z＝10m处仍可观测到甲烷辨识体积分数，但不足以构成威胁。

由图7可知，当风速为10m／s时，在z＝1m平面上，障碍物近壁面甲烷辨识体积分数聚积宽度达到18m，但危险体积分数范围分布高度明显下降，在z＝5m处已观测不到危险体积分数范围，分布高度在z＝3m以下。

①在低风速(1 m／s)时，天然气泄漏对2～3层住户威胁最大；在中风速(5m／s)时，天然气泄漏对1～2层住户威胁最大；在高风速(10m／s)时，天然气泄漏对1层住户威胁最大，且易在近地面处迅速扩散。

②当风速较小时，天然气不易扩散，影响范围比较大。当风速增大时，天然气易被稀释到爆炸下限以下，虽然扩散更广，但造成的危害反而减小。因此，在低风速时，更要严密关注，以免造成更大的损失。

③在天然气发生泄漏后，受到威胁的楼层居民应尽快疏散，不要接打电话、开关电器，以免产生电火花，更不能吸烟、使用打火机。

④本文的模拟是建立在简化物理模型上的，在实际中，住宅小区地形更为复杂，建筑外形及大气参数更为多变。我们需进一步借助其他工具以及根据实际情况按比例设计物理模型描绘复杂地形，提高模拟的准确性。

 

[1]孙洁．室外管道燃气泄漏扩散模拟与可视化研究(硕士学位论文)[D]．重庆：重庆交通大学，2012：1-5．

[2]刘欣．可燃气泄漏扩散的CFD研究[J]．计算机与应用化学，2012，29(2)：185-190．

[3]李百战，张甫仁．室内天然气泄漏三维浓度场的CFD模拟[J]．煤气与热力，2012，32(1)：B09-Bl2．

[4]程浩力，刘德俊，龙世华，等．城市燃气管道泄漏的CFD模拟[J]．油气储运，2011(8)：647-651．

[5]吴晋湘，牛坤，闫运忠．燃气连续性泄漏扩散规律的研究[J]．河北工业大学学报，2007，36(3)：1-6．

[6]肖淑衡．厂区燃气等危险性气体泄漏扩散的模拟研究(硕士学位论文)[D]．广州：广州大学，2006：7-8．

[7]LIU X P，NIU J L，KWOK K C S．Analysis of coneentration fluctuations in gas dispersion around high-rise building for different incident wind directions[J]．Journal of Hazardous Materials，2011(3)：1623-1632．

[8]YANG Rui，ZHANG Jing，SHEN Shifei，et al．Numerical investigation of the impact of different configurations and aspect rations on dense gas dispersion in urban street canyons[J]．Tsinghua Science and Technology，2007(3)：345-351．

[9]陶文铨．数值传热学[M]．3版．西安：西安交通大学出版社，2001：332-350．