云团扩散模型之欧阳地创编

1、欧阳地创编1 云团扩散模型根据物质泄漏后所形成的气云的物理性质的不同，可以将描述气云扩散的模型分为非重气云模型和重气云模型两种1.1 非重气云模型高斯模型是一种常用的非重气扩散模型，高斯烟羽(plume model)模型又称高架点连续点源扩散模型，适用于连续源的扩散，即连续源或泄放时间大于或等于扩散时间的扩散。高斯烟团(puff model)模型适用于短时间泄漏的扩散，即泄放时间相对于扩散时间比较短的情形，如突发性泄放等。若假设气体云内空间上的分布为高斯分布，则地面地处风向的烟团浓度分布算式为5-13。式中，c(x，y，h)点(x，y，h)处浓度值，mg/m ；q源强，即单位时问的排放量，mg

2、/s；u环境平均风速，m/s；x，y，z扩散参数；h源高(烟团高度)，m；x下方向到泄漏原点的距离，m；y，z侧风方向、垂直向上方向离泄漏原点的距离，m。高斯模式的实际应用效果很大程度上依赖于如何给定模式中的一些参数，尤其要注意源强、扩散参数等的确定。源强与污染物的物理化学属性、扩散方式、释放点的地理环境等有关。扩散参数表征大气边界层3欧阳地创编欧阳地创编内湍流扩散的强弱，是高斯模式的一项重要数据。高斯扩散模式所描述的扩散过程(实质上也包含了在实际应用中对高斯模式的一些限制)主要有：1)下垫面平坦、开阔、性质均匀，平均流场稳定，不考虑风场的切变。2)扩散过程中，污染物本身是被动、保守的，即污染

3、物和空气无相对运动，且扩散过程中污染物无损失、无转化，污染物在地面被反射。3)扩散在同一温度层结中发生，平均风速大于1.0 m/s。4)适用范围一般小于 1020 km。1.2 重气云模型由于重气本身的特殊性，在重气扩散领域也有大量基于不同理论的模型。鉴于重气扩散与中性或浮性气体扩散有着明显的区别，目前国内外已开发大量的不同复杂程度的重气扩散模型，如箱模型、相似模型、ltahgdm 模型、cfd 模型等。1.2.1 箱(box)模型箱模型是指假定浓度、温度和其他场，在任何下风横截面处为矩形分布等简单形状，这里的矩形分布是指在某些空间范围内场是均匀的，而在其他地方为零。该类模型预报气云的总体特征

4、，如平均半径、平均高度和平均气云温度，而不考虑其在空间上的细节特征。重气效应消失后其行为表现为被动气体扩散，所以该类模型还包括被动扩散的高斯模型及对它的修正。1.2.2 层流及湍流大气环境中的重气扩散(ltahgdm)模型ltahgdm 模型(heavy gas dispersion model in lsaminar and turbulent atmosphere层流及湍流大气环境中的重气扩散模型)以箱模型为基础，结合虚点源模型，能描述重气泄漏扩散整个过程。模型同三维有限元模型相比，具有形式简单、原始输入数据运算速度快等优点。ltahgdm 模型的建立基于以下几点假设：1)危险性气体初时泄

5、漏时，其外形呈正圆柱形(h=2r)。2)初始时刻泄漏源即此核电站内部的浓度、温度呈均匀分布。欧阳地创编欧阳地创编3)扩散过程不考虑泄漏源即此核电站内部温度的变化，忽略热传递、热对流及热辐射。4)泄漏气体认为是理想气体，遵守理想气体状态方程。5)在水平方向上，大气扩散系数呈各向同性。6)整个扩散过程中风速的大小、方向保持不变。7)地面对泄漏气体不吸收。8)整个过程中不发生任何化学反应等。放射性气体的扩散受其自身重力沉降引起的湍流及周围大气的湍流的双重影响。随着扩散的进行，放射性气体的浓度被稀释，重气效应逐渐消失，大气湍流逐渐成为控制此放射性气体扩散的主要因素。假设此放射性气体排出时的半径为r，高

6、度为 h。认为放射性气体排出的静压头等于空气的动力拖拽，则其径向尺寸变化率为在式(2)等温流动或式(3)扩散气体与空气具有相同的摩尔比热及地面加热可以忽略的非等温流动情况下，b 可认为是一常数，其值等于b0。欧阳地创编欧阳地创编重气云团的顶部空气卷吸和侧面空气卷吸对于云团的稀释是非常重要的。在垂直方向，由于云团顶部的空气卷吸和重力沉降的作用，使云团在垂直方向上的浓度分布呈现出从顶部到底部逐渐变大的高斯分布；在水平方向，由于侧面空气卷吸的作用，云团边缘也会形成高斯状浓度分布区，但由于重力沉降的原因，云团半径逐渐变大，侧面空气卷吸作用不会很快影响到云团内部，因此，可以假设在半径为rc 的区域内，浓

7、度均匀分布。云团内部的浓度可表示为随着重气效应的消失，大气湍流逐渐控制云团的扩散，rc 逐渐变小，最终为零。此时整个云团内部成高斯分布，可按照高斯烟团模型进行相关计算。判断重气云团向非重气云团转变的可以利用尉准则，当ri小于临界 richardsion 数时，重气云团已经转变为非重气云团。文中ric 取为 0.1。欧阳地创编欧阳地创编由于是在大气湍流环境下的扩散，因此，扩散系数来自于重气沉降引起的湍流扩散和环境湍流扩散两方面：一般认为云团高度就是箱模型中所假设的圆柱形的高度，即：关于 ra的计算，c.s.matthias 通过理论及实验分析，给出了如下的计算公式：式中，rc云团核心半径，m；h

8、云团高度，m；v云团体积，m ；0云团初始密度，kg/m ；33欧阳地创编欧阳地创编v0云团初始体积，m ；h0云团初始高度，m；t云团扩散时间，s；l云团特征尺寸，m；rg重力沉降引起的径向扩散系数，m；zg重力沉降引起的垂直扩散系数，m；a1云团重力沉降系数；ric临界 richardsion 数；ririchardsion 数；r云团半径，m；g重力加速度，m/s ；a空气密度，kg/m ；云团内部密度，kg/m ；r0云团初始半径，m；d0云团初始直径，m；0云团与周围空气初始密度差；云团扩散特征时间，s；r，z预测点圆柱坐标，m；3323欧阳地创编欧阳地创编ra大气湍流引起的径向扩散

9、系数，m；za大气湍流引起的垂直扩散系数，m；a2，c1，c2经验常数；ril特征 richardsion 数；u环境风速，m/s。2 系统设计及功能总体设计的任务是根据目标系统的物理模型确定一个合理的软件系统的体系结构。该易燃易爆毒性气体扩散模拟系统分为高斯模型模块、box 模型模块、ltahgdm 模型模块，其中：1)高斯模型模块由扩散浓度随距离变化的模拟、带有最小安全距离和扩散浓度值的模拟、固定距离浓度值计算模块组成。2)box 模型模块、ltahgdm 模型扩散半径随时间变化的模拟、扩散浓度随时间变化模拟、扩散浓度随距离变化模拟、固定距离浓度值计算模块组成。具体系统的功能结构如图1 所

10、示。欧阳地创编欧阳地创编图 1 扩散模拟系统功能结构系统主功能界面及高斯模型、ltahgdm 模型模拟界面如图 2、图 3、图 4 所示。图 2 系统主界面欧阳地创编欧阳地创编图 3 高斯模型模拟界面欧阳地创编欧阳地创编图 4 ltahgdm 模型模拟界面3 软件应用3.1 高斯模型的应用3.1.1 初始条件以氯气为例，假设某化工厂室外有一储罐，罐内压力为0.9 mpa，温度为 15，分子量为 0.03545kg/mol，绝热指数 1.310，假设由于罐体破裂发生连续型泄漏，泄漏口面积为0.02 m ，在一个阴天的夜晚储罐发生泄漏，有效泄漏高度为6 m，根据当地气象条件，风速为2.1 m/s。

11、3.1.2 假设条件设风速方向为 x 轴方向，泄漏源中心地面投影为坐标点，假定流场稳定，则扩散符合烟羽模型。2欧阳地创编欧阳地创编假设该大气稳定度为 d，泄漏源强为 5.341 kg/s。3.1.3 模拟计算1)在图 3 中的相应的文本框中输入对应的参数，如物质选择为氯气，泄漏源强为5.341 kg/s，平均风速为 2.1 m/s，有效泄漏高度为 6 m，选择大气稳定度为 d，点击不同的按钮，就可得到相应的模拟结果，如在固定高度输入 1.5 m，点击“下风向固定高度不同距离扩散浓度值”按钮，其结果如图5 所示，曲线表示下风向 1.5 m 高处不同距离的扩散浓度。图 5 下风向 1.5 m 高处

12、的扩散浓度曲线2)保持以上参数，点击“查看最高允许浓度并显示最小安全距离”按钮，即可显示所评价物质的最高允许浓度，如氯气的最大允许浓度1 mg/m ，并根据此浓度模拟出安全疏散所需要的最小安全距离，如图 6 所示。3欧阳地创编欧阳地创编图 6 人群疏散的最小安全距离3)保持以上参数，输入相应的下风向距离，即可计算固定高度在该距离下的具体浓度。如输入下风向距离 125 m，点击“确定”按钮，即可得出该距离下的浓度值0.0021 kg/m ，如图 7 所示。3欧阳地创编欧阳地创编图 7 125m 处的浓度值3.2 ltahgdm 模型的应用3.2.1 初始条件以大连市某韩资企业内的液化气瓶组站发生

13、泄漏为例，该瓶组站内共有50 kg 液化天然气钢瓶 8台，选取其最危险状态即液化气钢瓶破裂导致瓶组站内的所有液化气全部瞬时泄漏，相关气象资料根据该公司提供的资料查得。由于lng 主要成分甲烷的质量分数在90%以上，天然气泄漏后很难计算混合物的相关状态，因此，将 lng 看作甲烷计算。3.2.2 模拟计算1)在图 4 的相应的文本框中输入对应的参数，如初始半径为4 m，初始高度为 8 m，云团初始浓度为 100 mg/m ，气云密度为 3 kg/m ，空气密度为 1.29 kg/m ，云团重力沉降系数为0.7，点击“查看扩散半径随时间变化图”按钮，即得出云团扩散半径随时间变化的模拟曲线，如图8

14、所示。333图 8 扩散半径随时间变化的模拟欧阳地创编欧阳地创编2)保持以上参数，在下风向距离文本框中输入数值，如15 m，点击“查看扩散浓度时间变化图”按钮，即可得出相同距离15 m 下，不同扩散时间上的浓度扩散模拟图，如图9 所示。图 9 不同时间上的浓度扩散模拟图3)保持以上参数，在云团扩散时间文本框中输入数值，如2 s，点击“查看扩散浓度随时间变化图”按钮，即可得出在相同扩散时间2 s 下，不同下风向的浓度扩散模拟图，如图10 所示。4)保持以上参数，在下风向距离文本框中输入数值，如25 m，点击“计算”按钮即可得出25 m 处的浓度值为 1.0479 mg/m ，如图 11 所示。3

15、欧阳地创编欧阳地创编图 10 不同距离的浓度散模拟图图 11 25m 处对应的浓度值4 结论1)运用非重气扩散模型中的高斯模型和重气扩散模型中的box 模型、ltahgdm 模型进行软件开欧阳地创编欧阳地创编发，得出了相关物理量变化曲线。2)系统实现了非重气云团扩散模型中的高斯模型模块，该模块以氯气工厂储灌泄漏扩散为例，完成了固定高度下不同距离上的扩散浓度值曲线的模拟和不同物质最小安全距离曲线的模拟，同时求出了下风向 125 m 处的浓度 0.0021 kg/m 。3)系统实现了重气云团扩散模型中的ltahgdm 模型模块，该模块以大连市某韩资企业内的液化气瓶组站发生泄漏为例，完成了扩散半径随

16、时问变化曲线的模拟、固定时间扩散浓度随距离变化曲线的模拟以及对固定距离扩散浓度随时问变化曲线的模拟，同时计算出固定扩散时间下风向25 m 处的浓度 1.0479mg/m 。4)使用的语言工具是 java 语言，传统的软件往往与具体的实现环境(操作系统，如 linux，win-dows)有关，一旦环境有所变化就需要对软件进行一番改动，耗时耗力，而java 编写的软件能在执行码上兼容。这样，只要计算机提供了java 解释器，java 编写的软件就能在其上运行。参考文献1刘诗飞，詹予忠. 重大危险源辨识及危害后果分析m. 北京：化学工业出版社，20042吴宗之，高进东. 危险评价方法与应用m. 北京

17、：冶金工业出版社，20013柳静献，常德强，林秀丽等. 作业场所职业危害监管信息管理系统开发j. 中国安全科学学报，2008，18(3)：1181224孟亦飞，蒋军成. 化工装置火灾、爆炸、毒物扩散危险快速辨识方法j. 中国安全科学学报，2007，17(6)：1081135潘旭海，蒋军成. 重气云团瞬时泄漏扩散的数值模拟研究j. 化学工程，2003，31(1)：35396张启平，麻德贤. 危险物泄漏扩散过程的重气效应j. 北京化工大学学报，1998，25(3)：86907蒋军成，潘旭海. 一种描述重气泄漏扩散过程的新型模型j. 南京工业大学学报，2002，24(1)：414633欧阳地创编欧阳地创编8潘旭海，蒋军成. 事故泄漏源模型研究与分析j. 南京工业大学学报，2002，24(1)：1051109王文娟，刘剑锋. 危险性气体泄漏扩散数学模拟研究j. 工业安全与环保，2006，32(11)：232410黄斌，陈建中，殷志伟. 液氯泄漏事故扩散范围预测j. 湖南有色金属，2007，23(5)：454811胡世明，吴宗之，张政. 危险物质意外泄漏的重气扩散数学模拟j. 劳动保护科学技术，2000，20(2)：313812h.w.m.witlox. the hegadas model for ground level heavy