高斯烟羽模型

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

2、过程中风速的大小、方向保持不变， 不随地点、时间变化而变化；地面对泄漏气体起全反射作用，不发生吸收或吸附作用；整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。3、模型公式推导由正态分布假设可以导出下风向任意一点Xx,y,z处泄漏气体浓度的函 数为：1X(x, y, z) = A(x)e ay 2 e bz21由概率统计理论可以写出方差的表达式为:js y 2 Xdyj Xdyz 2 Xdzjs Xdz由假设可以写出源强的积分公式:Q = j-s式中：b y、b z为泄漏气体在y、z方向分布的标准差，单位为m；Xx,y,z为任一点处泄漏气体的浓度，单位为kg/m3； u为平均风速，单

3、位为m/s； Q为源强即泄漏速度，单位为kg/s；将1式代入2式，1a =将1式代入2式，1a =2b 2yb =上 2b 2z将1式和4式代入3式，积分可得:人（x）= Q2兀ub b,5式代入1式，可得：再将4式和积分可得:4X(X(x, y, z)=exp2兀ub b上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的 存在，烟羽的扩散是有界的。根据假设可以把地面看做一镜面，对泄漏气体起全 反射作用，并采用像源法处理，原理如图2.2所示。可以把任一点p处的浓度看 做两部分的贡献之和：一部分是不存在地面时所造成的泄漏物浓度；一部分是由 于地面反射作用增加的泄漏物浓度。该处的泄

4、漏物浓度即相当于不存在地面时由 位于0,0,H的实源和位于0,0,-H的像源在P点处所造成的泄漏物浓度之 和。图2.2像源法原理示意图其中，实源的贡献为：Q / 1 y 21 (z - H )2Xy，z) =exp(-2旅)exp(-2-)y zyz其中，像源的贡献为：/、 Q / 1 y21(z + h)2X*y，z)= 2Sexp(-2ge硕-2-) y zyz则该处的实际浓度为：X(x, y, z)= X(x, y, z)+ X x, y, z)由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为：Q z 1 y2exp( - _ 二)x TOC o 1-5 h z 2兀 u b b

5、2 b 2( 1 (z - H)2: 1 (z + H)2：1exp( -) + exp(-) HYPERLINK l bookmark5 o Current Document 2 b 22 b 2其中：Xx,y,z为下风向x米、横向y米、地面上方z米处的扩散的气体 浓度，单位为kg/m3； Q为源强即源释放速率，单位为kg/s； u为平均风速， 单位为峪；。为水平扩散参数，单位加；。为垂直扩散参数，单位为m；七为泄漏后是时间，单位为s； H为泄漏源有效高度，单位为m； y为横向距离，单 位为m； z为垂直方向距离，单位为m。如：式中，令z=0，即可得到地面气体浓度计算公式：X(x, X(x,

6、 y,0, H)= Q 兀ub blz y、exp -()22 b1-y(H 2exp 2b 2z(11)令y=z=0，即可得到地面轴线气体浓度计算公式:QH 2X (x,0,0, H) = -_2- exp(- -)12y zz其中，Xx,y,0为下风向x米、横向y米处的地面扩散气体浓度，单位为 kg/m3；若令y=0，则可以得到下风向中心线上的浓度分布。4、泄漏源有效高度烟云抬升高度的计算：以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的 时候气云中心的离地高度。实际上，泄漏源有效高度就等于泄漏源几何高度加泄 漏烟云抬升高度。影响烟云抬升高度的因素有很多，主要包括：泄漏气体

7、的初始速度和方向、 初始温度、泄漏口直径、环境风速与风速岁高度的变化率、环境温度与大气稳定 度。有效源高：H = HS +AH其中，H s为泄漏源几何架高，AH为烟云抬升高度。 实验表明，泄漏源抬升高度可以用下面公式近似计算：AH = 2.4V d / V13其中，,是气云释放速度，单位为m/s； d是泄漏出口直径，单位为m； V 为环境风速，单位为m/s；13式是20世纪70年代末、80年代初，Wilson根据管道破裂泄漏实验 所得的实验公式。实验时气体的喷射方向与风速垂直并且气体喷射路径上无障碍 物。实验表明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱2以内。计算出泄漏烟云抬升高度

8、以后，将泄漏源抬升高度与泄漏源实际几何高度相 加就得到了泄漏源有效高度。5、扩散系数的选取：扩散系数b x、b y、b z扩散系数b x、按照Pasquill的分类方法，随着气象条件稳定性的增加，大气稳定度可以分为A、B、C、D、E、F六类。其中A、B、C三类表示气象条件不稳定，E、F两 类表示气象条件稳定，D类表示中性气象条件，也就是说气象条件的稳定性在稳 定和不稳定之间。A、B、C三种类型的稳定度中，A类表示气象条件极其不稳定， B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定。E和F两种类 型的稳定度中，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定。大 气稳定度具体分类方

9、法见下表1、表2。表1 Pasquill大气稳定度的确定地面风速m/s白天日照夜间条件强中等弱阴天且云 层薄，或 低空云量 为4/8天空云量 为3/8 6CDDDD表2日照强度的确定天空云层的情况日照角 60日照角35日照角15。且35天空云量为4/8，或高空有薄云强中等弱天空云量为 5/8-7/8，云层高 度为 2134-4877m中等弱弱天空云量为 5/8-7/8，云层高 度2134m弱弱弱表1和表2中的云量是指当地天空层覆盖率。例如，云量为3/8是指当地 3/8的天空有云层覆盖。日照角是指当地太阳光线与地平线之间的夹度。例如， 阳光垂直照射地卖弄时的日照角为90。一般来说，随着大气稳定度的增加，扩散系数减小。根据Hanna和Drivas 的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表 3。表3扩散系数的计算方法太气稳定度A0.22a/(I 十叫OOLr)妲OJj：B0.1女/（1+。颂吹）”D.2aCOHj