高斯烟羽模型

1、高斯烟羽模型模型假设：1、坐标系高斯模型的坐标系如图2.1所示，原点为排放点(若为高架源，原点为排 放点在地面的投影)，x轴正向为风速方向，y轴在水平面上垂直于x轴，正向 在x轴的左侧，z轴垂直于水平面xoy，向上为正向。在此坐标系下烟流中心线 或烟流中心线在xoy面的投影与x轴重合。9 2,!再斯模式的坐标系2、模型假设(1) 污染物的浓度在y、z轴上的分布是高斯分布(正态分布)的；(2) 污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀 分布；(3) 扩散过程中不考虑云团内部温度的变化， 忽略热传递、热对流及热辐射;(4) 泄漏气体是理想气体，遵守理想气体状态方程；(5) 在

2、水平方向，大气扩散系数呈各向同性；(6) 取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变， 不随地点、时间变化而变化；(7) 地面对泄漏气体起全反射作用，不发生吸收或吸附作用；(8) 整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。3、模型公式推导由正态分布假设可以导出下风向任意一点X(x,y,z )处泄漏气体浓度的函数为：X(x, y, z) A(x)e ay e bZ由概率统计理论可以写出方差的表达式为:0 y2XdyXdy0z2Xdz2 0zXdz0由假设可以写出源强的积分公式:(3)uXdydz式中：y、 z为泄漏气体在y、Z方向分布的标准差，单位为m； X（

3、x,y,z为任一点处泄漏气体的浓度，单位为 kg/m 3 ； u为平均风速，单位为m/s ； Q为 源强（即泄漏速度），单位为kg/s ；将（1）式代入（2）式，积分可得:将（1）式和A( x)12z（4）式代入（3）式，积分可得:Q(4)再将（4）式和U y z(5)式代入(1)式，可得:(5)X( x, y, z)Qexp2 U y z2 2y z2 22 y2 z(6)上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面 的存在，烟羽的扩散是有界的。根据假设可以把地面看做一镜面，对泄漏气体 起全反射作用，并采用像源法处理，原理如图2.2所示。可以把任一点p处的浓度看做两部分的

4、贡献之和：一部分是不存在地面时所造成的泄漏物浓度；一 部分是由于地面反射作用增加的泄漏物浓度。该处的泄漏物浓度即相当于不存 在地面时由位于（0,0,H ）的实源和位于（0,0,-H ）的像源在P点处所造成的泄 漏物浓度之和俱資主输 实遒他。厨）斥8 4二F匕-z)附2.2像源法原理示意图其中，实源的贡献为：Xl(x,y,z) 血旳(21 y；)exp(y22 -)z1(z H)2其中，像源的贡献为：QX2（x, y, z）exp（2 U y z劣)exp(2 y(8)则该处的实际浓度为：X（x, y,z）X/x, y,z）X2（x, y,z）由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公

5、式为:2X （x,y,z,t,H ） Qexp（ 比）2 U y z2222 -）z地面上方(9)exp（寸讣）exp（2 z 2其中：X（x,y,z ）为下风向x米、横向y米、体浓度，单位为kg/m3 ； Q为源强（即源释放速率）,单位为(10)z米处的扩散的气kg/s ; u为平均风t为泄漏后是时间，单位为s； h为泄漏源有效高度，单位为 m y为横向距离, 单位为m z为垂直方向距离，单位为 m。女口：式中，令z=0,即可得到地面气体浓度计算公式：X(x, y,0, H) expz2y2exp（11）令y=z=0，即可得到地面轴线气体浓度计算公式:X(x,0,0,H)exp(H2（12）

6、其中，X（x,y,0 ）为下风向x米、横向y米处的地面扩散气体浓度，单位 为kg/m3;若令y=0,则可以得到下风向中心线上的浓度分布。4、泄漏源有效高度（烟云抬升高度的计算）：以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的 时候气云中心的离地高度。实际上，泄漏源有效高度就等于泄漏源几何高度加 泄漏烟云抬升高度。影响烟云抬升高度的因素有很多，主要包括：泄漏气体的初始速度和方向、 初始温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳 定度。有效源高：H Hs H其中，H $为泄漏源几何架高，H为烟云抬升高度。实验表明，泄漏源抬升高度可以用下面公式近似计算：(13

7、)H2.4VSd /V其中，Vs是气云释放速度，单位为 m/s； d是泄漏出口直径，单位为 m； V 为环境风速，单位为 m/s；（13）式是20世纪70年代末、80年代初，Wils on根据管道破裂泄漏实验 所得的实验公式。实验时气体的喷射方向与风速垂直并且气体喷射路径上无障 碍物。实验表明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱 以内。计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度5、扩散系数的选取:扩散系数x、 y、 z的大小与大气湍流结构、离地高度、地面粗糙度、泄漏持续时间、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。大气的湍流结构和风速在大气稳定度

8、中考虑。大气稳定度由10米高度以上的风速、白天的太阳辐射或夜间的云量等参数决定。按照Pasquill的分类方法，随着气象条件稳定性的增加，大气稳定度可以 分为A、B、C、D E、F六类。其中A B C三类表示气象条件不稳定，E、F两 类表示气象条件稳定，D类表示中性气象条件，也就是说气象条件的稳定性在稳 定和不稳定之间。A B、C三种类型的稳定度中，A类表示气象条件极其不稳定， B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定。E和F两种类 型的稳定度中，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定。大 气稳定度具体分类方法见下表1、表2。表1 Pasquill大气稳定度的确定

9、地面风速（m/s）白天日照夜间条件强中等弱阴天且云 层薄，或 低空云量为4/8天空云量为3/82AA-BB2 3A-BBCEF34BB-CCDE4 6CC-DDDD6CDDDD表2日照强度的确定角 召“ O 日八6 日O53 角 照且 日O60VO53 角v 照且 日O544云 4薄 为协 量空 云高 空或 天强等 中弱、咼77 为层87 1里云4-4 云，3 空偲21 天肓为等 中弱弱高 为层m 1里云34 云，21 空偲t 天肓B弱弱弱表1和表2中的云量是指当地天空层覆盖率。例如，云量为3/8是指当地3/8的天空有云层覆盖。日照角是指当地太阳光线与地平线之间的夹度。例如， 阳光垂直照射地卖弄时的日照角为 90。一般来说，随着大气稳定度的增加，扩散系数减小。根据Hanna和Drivas的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下 表3。表3扩散系数的计算方法大气稳定度JA022(“000010.2xB0.16x/(l + 0.0001x)5O.12jC0.11x