云团扩散模型

1、1云团扩散模型根据物质泄漏后所形成的气云的物理性质的不同，可以将描述气云扩散的模型分为非重气云模型和重气云模型两种5-13。1.1非重气云模型高斯模型是一种常用的非重气扩散模型，高斯烟羽(Plume model)模型又称高架点连续点源扩散模型，适用于连续源的扩散，即连续源或泄放时间大于或等于扩散时间的扩散。高斯烟团(Puff model)模型适用于短时间泄漏的扩散，即泄放时间相对于扩散时间比较短的情形， 如突发性泄放等。若假设气体云内空间上的分布为高斯分布，则地面地处风向的烟团浓度分布算式为£(*.¥、田=冬X exp 。(/一?-)+4 + 三)(0(2tt产%2%/式中

2、，c(x , y, H)点(x , y, H)处浓度值,mg/mQ-一源强，即单位时问的排放量，mg/s；u环境平均风速， m/s;bx, by, z z 扩散参数；H源高(烟团高度)，x 下方向到泄漏原点的距离，m;y, z 侧风方向、垂直向上方向离泄漏原点的距离，ni高斯模式的实际应用效果很大程度上依赖于如何给定模式中的一些参数，尤其要注意源强、扩散参 数等的确定。源强与污染物的物理化学属性、扩散方式、释放点的地理环境等有关。扩散参数表征大气边界层内湍流扩散的强弱，是高斯模式的一项重要数据。高斯扩散模式所描述的扩散过程（实质上也包含了在实际应用中对高斯模式的一些限制 ）主要有：1）下垫面平

3、坦、开阔、性质均匀，平均流场稳定，不考虑风场的切变。2）扩散过程中，污染物本身是被动、保守的，即污染物和空气无相对运动，且扩散过程中污染物无 损失、无转化，污染物在地面被反射。3）扩散在同一温度层结中发生，平均风速大于1.0 m/s 。4）适用范围一般小于 1020 km。1.2重气云模型由于重气本身的特殊性，在重气扩散领域也有大量基于不同理论的模型。鉴于重气扩散与中性或浮 性气体扩散有着明显的区别，目前国内外已开发大量的不同复杂程度的重气扩散模型，如箱模型、相似模型、LTA HGDMI 型、CFD模型等。1.2.1 箱（BOX）模型箱模型是指假定浓度、温度和其他场，在任何下风横截面处为矩形分

4、布等简单形状，这里的矩形分 布是指在某些空间范围内场是均匀的，而在其他地方为零。该类模型预报气云的总体特征，如平均半径、平 均高度和平均气云温度，而不考虑其在空间上的细节特征。重气效应消失后其行为表现为被动气体扩散，所 以该类模型还包括被动扩散的高斯模型及对它的修正。1.2.2 层流及湍流大气环境中的重气扩散（LTA HGDM）1型LTA HGDIMI 型（Heavy Gas Dispersion Model in Lsaminar and Turbulent Atmosphere层流及湍流大气环境中的重气扩散模型）以箱模型为基础，结合虚点源模型，能描述重气泄漏扩散整个过程。模型同三维有限元模

5、型相比，具有形式简单、原始输入数据运算速度快等优点。LTA HGDMI型的建立基于以下几点假设：1）危险性气体初时泄漏时，其外形呈正圆柱形（H=2R）。2）初始时刻泄漏源即此核电站内部的浓度、温度呈均匀分布。3）扩散过程不考虑泄漏源即此核电站内部温度的变化，忽略热传递、热对流及热辐射。4）泄漏气体认为是理想气体，遵守理想气体状态方程。5）在水平方向上，大气扩散系数呈各向同性。6）整个扩散过程中风速的大小、方向保持不变。7）地面对泄漏气体不吸收。8）整个过程中不发生任何化学反应等。放射性气体的扩散受其自身重力沉降引起的湍流及周围大气的湍流的双重影响。随着扩散的进行， 放射性气体的浓度被稀释，重气

6、效应逐渐消失，大气湍流逐渐成为控制此放射性气体扩散的主要因素。假设此放射性气体排出时的半径为R,高度为H。认为放射性气体排出的静压头等于空气的动力拖拽，则其径向尺寸变化率为¥=%产作（2）b-gV/r（3）在式（2）等温流动或式（3）扩散气体与空气具有相同的摩尔比热及地面加热可以忽略的非等温流动情况下，b可认为是一常数，其值等于 boo(4)b0 =gVoAo/TF其中，& =(PLP,)P。对式(2)积分,得到：/=g+2组”力(5)云团体积为ym邛+。命其中，公“了重气云团的顶部空气卷吸和侧面空气卷吸对于云团的稀释是非常重要的。在垂直方向，由于云团顶部的空气卷吸和重力沉降

7、的作用，使云团在垂直方向上的浓度分布呈现出 从顶部到底部逐渐变大的高斯分布；在水平方向，由于侧面空气卷吸的作用，云团边缘也会形成高斯状浓度分布区，但由于重力沉降的原因，云团半径逐渐变大，侧面空气卷吸作用不会很快影响到云团内部，因此，可以假设在半径为Rc的区域内，浓度均匀分布。云团内部的浓度可表示为C = C© 叼(7)C = :©同"户的(8)式中，gs代可由牛顿迭代法求解：R2 =+ (2ir)l/2Hca jt + 2a,(9)Rc逐渐变小，最终为零。此时整个云团内部当R小于临界 Richardsion 数时，重气云团已随着重气效应的消失，大气湍流逐渐控制云团

8、的扩散, 成高斯分布，可按照高斯烟团模型进行相关计算。判断重气云团向非重气云团转变的可以利用尉准则，经转变为非重气云团。文中Ric取为0.1。由于是在大气湍流环境下的扩散，因此，扩散系数来自于重气沉降引起的湍流扩散和环境湍流扩散两方面:(10)般认为云团高度就是箱模型中所假设的圆柱形的高度，即:(J, =小 nJ(H)九 二% (犬-勺)J关于 ha的计算，C.S.Matthias 通过理论及实验分析，给出了如下的计算公式:其中,其中.<7父="刖1 +6)Ui/q 附"茕恒喂吊=心(12)式中，风可利用如下公式求得:Ri =p/o- L 2y心；1.2)式中，RC云

9、团核心半径，mH云团高度，m；v云团体积，m5；p 0云团初始密度，kg/m3；V)云团初始体积，m3;H0云团初始高度，mt 云团扩散时间，s;l云团特征尺寸，mb rg 重力沉降引起的径向扩散系数，成b zg重力沉降引起的垂直扩散系数，成ai云团重力沉降系数；Ric临界 Richardsion 数；RiRichardsion 数；R云团半径，m;g重力加速度，m/s2;p a空气密度，kg/m3； 云团内部密度，kg/m3；R云团初始半径，m；D)云团初始直径，m； 0云团与周围空气初始密度差； 云团扩散特征时间，s;r, z预测点圆柱坐标，m;dra 大气湍流引起的径向扩散系数，m；za

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

11、模拟系统功能结构系统主功能界面及高斯模型、LTA HGDhgl型模拟界面如图 2、图3、图4所示。LMT/II模型图2系统主界面高斯IK型x菱商双工w*大弋rt员容立家小空3HL .图3高斯模型模拟界面ITI嘏巾里LT4-IICW*更¥微槽疝计算-iMKAibKffr«RH(M*aM -f，立口）事用*为5停餐MTMKA图4 LTA HGDMII型模拟界面3软件应用3.1 高斯模型的应用3.1.1 初始条件以氯气为例，假设某化工厂室外有一储罐，罐内压力为 0.9 MPa,温度为15C,分子量为0.03545 kg/mol ,绝热指数1.310 ,假设由于罐体破裂发生连续型泄

12、漏，泄漏口面积为0.02 m2,在一个阴天的夜晚储罐发生泄漏，有效泄漏高度为6 m,根据当地气象条件，风速为 2.1 m/s 。3.1.2 假设条件设风速方向为x轴方向，泄漏源中心地面投影为坐标点，假定流场稳定，则扩散符合烟羽模型。假 设该大气稳定度为 D,泄漏源强为5.341 kg/s 。3.1.3 模拟计算1）在图3中的相应的文本框中输入对应的参数，如物质选择为氯气，泄漏源强为5.341 kg/s ,平均风速为2.1 m/s ,有效泄漏高度为 6 m,选择大气稳定度为D,点击不同的按钮，就可得到相应的模拟结果，如在固定高度输入1.5 m,点击“下风向固定高度不同距离扩散浓度值”按钮，其结果

13、如图5所示，曲线表示下风向1.5 m高处不同距离的扩散浓度。图5下风向1.5 m高处的扩散浓度曲线2）保持以上参数，点击“查看最高允许浓度并显示最小安全距离”按钮，即可显示所评价物质的最 高允许浓度，如氯气的最大允许浓度1 mg/R3,并根据此浓度模拟出安全疏散所需要的最小安全距离，如图 6所示。图6人群疏散的最小安全距离3）保持以上参数，输入相应的下风向距离，即可计算固定高度在该距离下的具体浓度。如输入下风 向距离125 m,点击“确定”按钮，即可得出该距离下的浓度值0.0021 kg/m 3,如图7所示。图7 125m处的浓度值3.2 LTA HGD膜型的应用3.2.1 初始条件以大连市某

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

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

16、。3 0 5 0 5. iM z 11 1 Ou d旦要谈S因浓度随卜风向距离变化曲线3.00 5 1015 2025 加 40 45 50 5560 65 7075 80 85 9095 100时间图10不同距离的浓度散模拟图结果板撤-速距离下浓度为：L0479 mMmmm图11 25m处对应的浓度值4结论Box模型、LTA HGD膜型进行软件开发,i）运用非重气扩散模型中的高斯模型和重气扩散模型中的得出了相关物理量变化曲线。2)系统实现了非重气云团扩散模型中的高斯模型模块，该模块以氯气工厂储灌泄漏扩散为例，完成 了固定高度下不同距离上的扩散浓度值曲线的模拟和不同物质最小安全距离曲线的模拟，

17、同时求出了下风向125 m 处的浓度 0.0021 kg/m 3。3)系统实现了重气云团扩散模型中的LTA- HGD喉型模块，该模块以大连市某韩资企业内的液化气瓶组站发生泄漏为例，完成了扩散半径随时问变化曲线的模拟、固定时间扩散浓度随距离变化曲线的模拟以 及对固定距离扩散浓度随时问变化曲线的模拟，同时计算出固定扩散时间下风向25 m处的浓度1.0479 mg/nt4)使用的语言工具是 JAVA语言，传统的软件往往与具体的实现环境(操作系统，如Linux ,win-dows)有关，一旦环境有所变化就需要对软件进行一番改动，耗时耗力，而JAVA编写的软件能在执行码上兼容。这样，只要计算机提供了JA

18、VA解释器，JAVA编写的软件就能在其上运行。参考文献1刘诗飞，詹予忠.重大危险源辨识及危害后果分析M.北京：化学工业出版社，2004 2吴宗之，高进东.危险评价方法与应用M.北京：冶金工业出版社，20013柳静献，常德强，林秀丽等.作业场所职业危害监管信息管理系统开发J.中国安全科学学报,2008, 18(3) : 118 1224孟亦飞，蒋军成.化工装置火灾、爆炸、毒物扩散危险快速辨识方法J.中国安全科学学报，2007, 17(6) : 108 113J. 化学工程，2003, 31(1) : 35-5潘旭海，蒋军成.重气云团瞬时泄漏扩散的数值模拟研究39-906张启平，麻德贤危险物泄漏扩散过程的重气效应J.北京化工大学学报，