1、气体扩散浓度计算模型介绍,华东理工大学 沈艳涛 2006.8.31,第一部分 扩散过程与模型分类介绍,相关背景污染性泄露,大气污染性泄露的形式： 自然方面：火山喷发的有害气体，某些物质自燃或在一定条件下产生的有毒气体，环境微生物产生的某些气体 日常生活方面：生活用煤产生的含氮硫氧气体 石化燃料动力的交通车辆产生的尾气将在一定气候下生成光化学雾 工业用气体的泄漏，特别是化学工业用到的大量的有毒有害，易燃易爆的气体 其他方面产生的一些气体及烟尘,可能造成的伤害,1、SOx,NOX,光气：口腔，呼吸道与肺部病变，皮肤病 2、液NH3,液Cl2等在管道破口喷射引起冻伤及化学毒性与环境危害 3、爆炸性气
2、体的爆炸性危害 4、,常见的泄露形式: 管道破损后的连续喷射烟羽,常见的泄露源: 爆炸形成瞬时泄露烟团,扩散过程研究,不同性质气体在不同条件下表现出不同的特征 观察者对过程特征的选取,重气扩散过程,四个阶段 初始阶段：物质从容器泄漏出，形成气云后在本身的惯性力和外界风速的作用下，上升变形； 重力沉降阶段和空气卷吸阶段：当气云初始动量消失后，重力占主导地位。由于云团与周围空气间的密度差，导致重气塌陷，沿地表面拓展，引起云团厚度的降低和径向尺寸的增大，而在大气湍流的作用下外界空气进入云团，即空气卷吸，云团被稀释，同时由于初始泄漏云团与周围环境的温度差异而进行热量交换； 非重气扩散转变：随着云团的稀
3、释冲淡，重气效应逐渐消失，重气扩散转变为非重气扩散； 大气湍流扩散阶段(被动扩散)：即大气湍流对云团的扩散起支配作用。,轻气扩散过程,物质从容器泄漏出，形成气云后在本身的惯性力和外界风速的作用下，上升变形。,中性气扩散过程,两个阶段 初始阶段：物质从容器泄漏出，形成气云后在本身的惯性力和外界风速的相互作用； 大气湍流扩散阶段：即大气湍流对云团的扩散起支配作用。,过程中变异性问题,源与边界的差异性及弱化 温度差异 密度差异,气体泄漏扩散研究方法,试验法 风洞实验法 模型法,试验法,试验数据,问题特点,模型,验证,特征提取与模化,风洞实验,比例,气体扩散浓度计算模型分类,重气泄漏扩散的数值模拟方法
4、依据各自的建模原理以及复杂程度可分为五类 : 第一类：唯像模型 第二类：箱及相似模型 第三类 浅层模式模型 第四类：三维模式模型 第五类 随机游走模式模型,第一类：唯像模型,唯像模型是由一系列图表或简单关系式来描述扩散行为的。 Britter and McQuaid在重气扩散手册中推荐了一套简单而实用的方程式和列线图，称之为B&M模型，他们是收集了许多重气扩散的实验室和现场实验的研究结果，以无因次的形式将数据连线并绘制成与数据匹配的曲线或列线图。,模型特点与适用,该模型比较简单，属于经验模型，外延性较差，可以用于确定工厂警戒线处产生主要影响的基本物理因素。 德国的VDI模型也采用了与BM模型类
5、似的处理方法。,B&M模型表达式,Cm，C0分别为气云横截面上的平均浓度、初始浓度，kgm-3； Vc0为连续烟流释放的初始气云体积流量，m3s-1； Vi0为瞬时烟团释放的初始气云体积，m3； u为10m高处的风速，ms-1； g0为初始的折算重力项，g0g(0-a)/a，0，a分别为初始气云密度和外界空气密度； fc,fi普遍化无因次函数。,第二类：箱及相似模型,箱模型和相似模型都是假定浓度、温度等在任何下风向横截面均满足一个简单分布，箱模型假定浓度、稳定等在箱内是均匀分布的，其它区域为0；而相似模型则假定模型内符合相似分布（如高斯分布）等简单形状。,模型常使用的相关假定,(1)危险性气体
6、初始泄漏时，其外形呈正圆柱形 ( H2R)或在某规则区域正态分布； (2)初始时刻云团内部的浓度、温度呈均匀分布； (3)扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射； (4)泄漏气体是理想气体，遵守理想气体状态方程； (5)在水平方向，大气扩散系数呈各向同性； (6)整个扩散过程中风速的大小、方向保持不变； (7)地面对泄漏气体不吸收； (8)整个过程中不发生任何化学反应等。,箱模型实例（by Van Ulden，1970）,对于重气瞬时泄漏形成的云团，一般把箱模型看作为一个圆柱形，如Van Ulden（1970年）提出将重气烟团当作一个初始体积为V0，初始高度为H0，初始
7、半径为R0的圆柱形，高度和半径随时间变化，与被动扩散的高斯模型相比，主要改进是考虑到云团的重力沉降现象，即在重力作用下，云团下沉，半径R增加，同时高度H减小。,箱模型实例（by Van Ulden，1970）,r，a为气云的“参考”密度和空气密度，kgm-3； K为常数。,对于重气连续泄漏形成的烟羽，一般把箱模型看作一个矩形，如Jagger在Fryer & Kaiser提出的烟团模型DENZ的基础上，开发了相应的烟流模型CRUNCH，用来模拟稳态连续泄漏。 模型假定高为H、宽为2L的矩形截面，原先半径和高度随时间变化的微分方程变成半宽和高度随下风距离变化的方程，原先径向重力扩散速度变成了侧向重
8、力扩散速度。,箱模型对重气研究基本假定,箱模型：重性气云早期扩展,扩散的过程中还考虑到周边空气的卷吸。早期的研究者在Van Ulden的重气云实验基础上，提出空气由模型的顶部卷吸进来是占主导作用的，随着更深入的认识，很多研究者，如Hanna & Drivas和Mcquaid都一致认为空气是从模型的顶部和侧面同时卷吸进来的，卷吸的速度受Richardson、纵向湍流速度、大气稳定度、风速、摩擦风速等影响。由于不同的箱模型采用了不同的空气卷吸参数，从而导致了不同的模式计算结果的差别是很大的。,随着云团的稀释冲淡过程，重气效应逐步地消失，当重气扩散转变为非重气扩散时，大气湍流对云团的扩散起支配作用，
9、云团的高度、半径及运行状态完全取决于大气湍流特性，实际上气体的浓度分布开始接近为高斯形状，仍然假定为均匀就不再合理。因此箱模型通常都有从均匀气云向高斯分布的转折点，即重气扩散向非重气扩散的转折点，采用理查逊数、沉降速度和速度尺度的关系，或者运用云团密度与周围空气的密度差来判断。,箱模型：重性气向非重气的转折,箱模型其他研究情况介绍,Manju（1995）在总结以前学者的研究基础上，开发了T Heavy Gas Models模型，可以用于模拟重气瞬时泄漏扩散和连续泄漏扩散。 扩散模型包括了重力沉降、空气卷吸、云团受热和向非重气云团过渡。 对扩散过程中的重力沉降系数、顶部卷吸系数和侧面卷吸系数的取
10、值进行了分析比较，提出了建议取值，并提出利用云团密度与周围空气密度差小于0.001Kg/m3来判断云团是否过渡为非重气云团。,模型验证情况,T Heavy Gas Models瞬时泄漏扩散模型对Thorney Island Tests系列试验下风向不同距离的泄漏物质最大浓度进行了模拟验证，T Heavy Gas Models连续泄漏扩散模型对Maplin Sands Tests系列试验下风向不同距离的泄漏物质最大浓度进行了模拟验证，两个试验的模拟结果都是较好的，基本上反映了重气的扩散情形。,验证试验,相似模型介绍,相似模型主要是针对HEGADAS以及以HEGADAS为基础开发的模型。相似模型是
11、对箱模型概念的扩展，考虑了气云内部浓度和速度的分布，并采取了湍流扩散系数而非空气卷吸速度的方法。 壳牌公司HEGADAS模型是HGSYSTEM系统软件包的重要组成部分，HEGADAS模型即有处理稳定连续释放的定常态版本，也有预报来自液化气液池蒸发在中等或高风下扩散的瞬间版本。 DEGADIS模型是在HEGADAS模型基础上作的改进，是美国海岸警卫队和气体研究院开发的。,相似模型特点与适用,箱及相似模型具有概念清晰、计算量较小等优点，可为危险评价、应急救援、制定控制措施等提供指导。但其自身也存在着局限性，如假定速度和浓度的相似分布，模拟的精度较差，重气云团向非重气云团过渡也存在着很大的不确定性。
12、,已开发的相似模型简介,SAFER、TRACE模型 （在Kaiser和Walker提出模型的基础上开发的） CONSEQ、PHAST、WHAZAN、SAFETI模型 （在Cox和Carpenter提出模型的基础上开发的） DENZ、DRIFT、CIGALE 2、SLOPEFMI模型 （在Fryer和Kaiser提出模型的基础上开发的） HEGADAS、HEGABOX、HGSYSTEM模型（Colenbrander） T HEAVY GAS MODELS模型（MANJU MOHAN , T. S. PANWAR 和M. P. SINGH） CHARM、ELOE模型（Eidsvik） 等等,第三类
13、 浅层模式模型,浅层模式模型，是对重气扩散的控制方程加以简化来描述其物理过程，是对于三维模式模型和简单箱模型的折衷。它是基于浅层理论（浅水近似）推广得到的，模型采用了厚度平均变量来描述流场特征，有利于考虑复杂地形的重气扩散情况。,浅层模式模型思想,该模型需要计算气云的宽度和高度，是拟三维的。侧风浓度分布应用相似分布确定，气云与环境大气的混合运用卷吸概念处理。不少专家对浅层模型进行了进一步开发，Wheatley & Webber对卷吸和热量传递的浅层模型进行了推导。Errnak等将浅层模型发展为SLAB模型，包括求解质量、组分、下风动量、侧风动量和能量的侧风平均守恒方程，以及气云宽度方程和理想气
14、体状况方程。,浅层模式模型适用性,常用于非互溶的流体中,此类模型有：,SHALLOW（Webber等） TWODEE（Hankin和Britter） DISPLAY1, DISPLAY2（wrtz等） 等等,第四类：三维模式模型,该模型采用计算流体力学（CFD）方法模拟重气得到三维非定常态湍流流动过程。这种数值方法是通过建立各种条件下的基本守恒方程（包括质量、动量、能量及组分等），结合一些初始条件和边界条件，运用数值计算理论和方法，求解NavierStokes方程，实现预报真实过程各种场的分布。,FEM3介绍,比较有代表的是FEM3（3-D Finite Element Model）模型，该模
15、型在1979年提出，最初是为了模拟LNG的突发性泄放，获得比较好的结果。 该模型所用有限元解法是由Galerkin法改进而来，主要可解不定常的连续性方程、动量方程、热量方程、扩散方程以及理想方程气体状态方程，用K理论（梯度输运理论）来处理湍流问题。,已经开发的模型,在此基础上一些科研人员已经开发了一些模型，有：SIGMET（英国）、TRANSLOC（Schnatz和Flothmann）、FEM3（Chan等）、ZEPHRY（Hertel和Teuscher）、MARIAH（Taft）、HEAVYGAS（Deaves）、MERCURE（Riou）、ADREA-HF（Bartzis）、MERADIS
16、（Vergison等）等等26。,三维模式模型优点,三维模式模型克服了箱及相似模型中辨别和模拟重力下沉、空气卷吸、气云受热等物理效应时所遇到的许多问题，能更好地描述重气在大气湍流运动中的物理现象，具有广泛的通用性，并且能够描述障碍物或明显的地形变化的复杂扩散过程。 该方法成本低廉，物理场景数学模型搭建灵活性强，因而日益受到研究者的重视。当前的研究不仅在于建立深入刻画过程的数学模型如N-S方程组及其封闭性边界构成的模型（如k-模式），计算过程如网格划分，收敛性讨论等也得到了深入的研究。与此同时，适用于工程目的的专业软件如fluent等已经受到了业界的认可和推崇。,三维模式模型缺点,此类模型明显的
17、缺点是模拟计算比较复杂，且有关数据计算要花费大量的计算机时。但是，随着计算机工业和计算理论的迅速发展和对过程认识的加深，这些问题在逐步得到解决。,三维模式模型表达式,建立对气体流动和扩散过程进行描述，包括如下方程 连续性方程 动量守恒方程 能量方程 密度方程(或气体状态方程) 湍流流动模型 初始条件和边界条件地确定对于三维重气流动和扩散过程的控制微分方程组有如下特点： 非线性 多变量 强耦合 定解条件复杂 所以不太可能用解析法求解，更实际的是将微分方程组离散化以后，用数值迭代法求解,三维模式模型解法,计算示例,计算结果,第五类 随机游走模式模型,长期以来，随机游走模型主要应用于中性气体扩散特征
18、模拟，随着计算能力的提高，随机游走模型已大量应用于点源扩散的数值模拟，逐渐发展到了浮力气体扩散，以及重气扩散模拟。,随机游走模式模型优缺点,随机游走模式在描述痕量物质在湍流中的传输具有许多优点：方程原理相对简单，易于计算机编程的实现；湍流的统计特性易于结合进方程中，易于实现湍流场特征模拟。 由于气象模式提供了粒子随机游走所依赖的平均风场及湍流特征资料，因此，精确的风场资料对于随机游走模式的正确模拟是至关重要的。而由于风场精度的限制，通常只能得到一些定性的结果 而难以得到定量的浓度特征。,HLY模型,化工部劳动保护研究所于1996年开发 泄漏源模型、5个泄漏模式、6个扩散模式和伤害分区模型, 其
19、中扩散的模型是在箱模型的基础上开发,EPA的相关研究,国内研究单位,中国经贸委安全工程中心（研究院） 北京环科院 南京工业大学 华东理工大学 中国科技大学 大连理工大学 其他,气体扩散模型选用考虑因素,扩散特性的影响因素： 密度差异 温度差异 源强度与类型 大气流动状况 地形建筑物等影响 其他特征,模型的选用原则,根据：问题特征的描述 所有的相关条件,相似准则,相似准则的推导 相似变换法 相似定数法 积分类比法 相似准则的物理意义 相似准则的性质 相似准则都是无量纲量 由微分方程直接导出的相似准则都具有明确的物理意义 相似系统对应点（稳态）和现象转换点（转换瞬间，非稳态）相似准则相等 任何两个
20、相似准则相乘或相除还是相似准则。 相似准则数值的大小反映了过程中质的变化。,模型的精度估计,解决问题的全面性 计算方法,第二部分 常见实用泄漏扩散模型介绍,几种典型源强估算,估算参数 1、泄漏速度 2、泄漏量 3、泄漏持续时间 4、泄漏物质的形态,常见的泄漏源： 工艺单元中液体经小孔泄漏 储罐中液体经小孔泄漏 液体经管道泄漏 气体或蒸气经小孔泄漏 闪蒸液体的泄漏 易挥发液体蒸发 其他,液体经小孔泄漏的源模式,机械能守恒方程：,不可压,a1,无机械功,D,一般孔流系数C01；; 安全起见C01.0,储罐中液体经小孔泄漏,z0 A,物质守恒方程,A0,液体经管道泄漏,查手册选取阻力系数，计算一般要
21、试算,气体或蒸气经小孔泄漏,理想气体状态方程绝热方程机械能守恒方程孔流系数定义,高压气体或蒸气经小孔泄漏时,临界压力点：流速音速,实际流速与压力关系,点源 烟流模式 烟团模式 特殊气象条件下的扩散模式 面源 等效点源 箱模式 手算模式 线源,短期平均浓度的计算 长期平均浓度的计算 颗粒物扩散计算 复杂下垫面地区的扩散计算,常见的扩散计算类型,烟流模式,有风（风速大于1m/s），风向稳定，连续点源地面浓度（z0），公式为下式,公式中的各参数可以参考相关文献。比如大气环境影响评价实用技术,地面轴线最大浓度求法,烟团模式,风速在静风时（风速小于1m/s）,特殊气象条件下的扩散模式,漫烟型扩散,面源,等效点源法,等效点源坐标：,箱模式,线源,平行风模式,谢谢大家！,
