高斯烟羽模型

模型假设：

1、坐标系

高斯模型的坐标系如图2.1所示，原点为排放点（若为高架源，原点为排放

点在地面的投影），x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x

轴的左侧，z轴垂直于水平面xoy,向上为正向。在此坐标系下烟流中心线或烟

流中心线在xoy面的投影与x轴重合。

图2.1高斯模式的坐标系

2、模型假设

（1）污染物的浓度在y、z轴上的分布是高斯分布（正态分布）的；

（2）污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀

分布；

（3）扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射;

（4）泄漏气体是理想气体，遵守理想气体状态方程；

（5）在水平方向，大气扩散系数呈各向同性;

(6)取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，

不随地点、时间变化而变化；

(7)地面对泄漏气体起全反射作用，不发生吸收或吸附作用；

(8)整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。

3、模型公式推导

由正态分布假设可以导出下风向任意一点X(x,y,z)处泄漏气体浓度的函

ay22

X(x9y9z)=A(x)e-e^

由概率统计理论可以写出方差的表达式为:

Lxdy

Cz2Xdz

Jo

由假设可以写出源强的积分公式:

uXdydz

ocJ-OQ

式中：%,、q为泄漏气体在y、z方向分布的标准差，单位为m;X(x,y,z)

为任一点处泄漏气体的浓度，单位为kg/n/；u为平均风速，单位为m/s;Q为

源强(即泄漏速度)，单位为kg/s;

将(1)式代入(2)式，积分可得:

1

a-

2矣

(4)

将（1）式和（4）式代入（3）式，积分可得:

A(x)=——⑸

再将（4）式和（5）式代入（1）式，可得:

z2

X(jv,y,z)=-------十五(6)

2mt(yv(y^24

上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的

存在，烟羽的扩散是有界的。根据假设可以把地面看做一镜面，对泄漏气体起全

反射作用，并采用像源法处理，原理如图2.2所示。可以把任一点p处的浓度看

做两部分的贡献之和：一部分是不存在地面时所造成的泄漏物浓度；一部分是由

于地面反射作用增加的泄漏物浓度。该处的泄漏物浓度即相当于不存在地面时由

位于（0,0,H）的实源和位于（0,0,-H）的像源在P点处所造成的泄漏物浓度之

和。

浓度，单位为kg/d；Q为源强(即源释放速率)，单位为kg/s;u为平均风速,

单位为m/s;区为水平扩散参数，单位为明％为垂直扩散参数，单位为叫t

为泄漏后是时间，单位为s;H为泄漏源有效高度，单位为明y为横向距离，单

位为m;z为垂直方向距离，单位为m。

如：式中，令z=0,即可得到地面气体浓度计算公式：

X(x,y,0,”)=—————exp-—(―)2exp(11)

,如巴巴L2%J[2a：)

令尸z=0,即可得到地面轴线气体浓度计算公式：

X(x,0,0,H)=Qexp(-^)(12)

其中，X(x,y,0)为下风向x米、横向y米处的地面扩散气体浓度，单位为

kg/Y;若令y=0,则可以得到下风向中心线上的浓度分布。

4、泄漏源有效高度(烟云抬升高度的计算)：

以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的

时候气云中心的离地高度。实际上，泄漏源有效高度就等于泄漏源几何高度加泄

漏烟云抬升高度。

影响烟云抬升高度的因素有很多，主要包括：泄漏气体的初始速度和方向、

初始温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳定

度。

有效源高：H=HS+\H

其中，“S为泄漏源几何架高，△”为烟云抬升高度。

实验表明，泄漏源抬升高度可以用下面公式近似计算：

(i3)

AH=2AVsd/V

其中，匕是气云释放速度，单位为m/s；d是泄漏出口直径，单位为m；V

为环境风速，单位为m/s；

(13)式是20世纪70年代末、80年代初，Wilson根据管道破裂泄漏实验所

得的实验公式。实验时气体的喷射方向与风速垂直并且气体喷射路径上无障碍

物。实验表明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱2以内。

计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度相

加就得到了泄漏源有效高度。

5、扩散系数的选取：

扩散系数,、%.、%的大小与大气湍流结构、离地高度、地面粗糙度、泄

漏持续时间、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。大气的湍

流结构和风速在大气稳定度中考虑。大气稳定度由1。米高度以上的风速、白天

的太阳辐射或夜间的云量等参数决定。

按照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。

表1Pasquill大气稳定度的确定

白天日照夜间条件

地面风速

阴天且云天空有量

(m/s)强中等弱

层薄，或为3/8

低空云量

为4/8

<2AA-BB

2~3A-BBCEF

3〜4BB-CCDE

4~6CC-DDDD

>6CDDDD

表2R照强度的确定

天空云层的情况日照角日照角日照角

>60°<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

的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表

3O

表3扩散系数的计算方法

大气稳定度

%巴

AO.22/a+O.OOOlx)”0.2x

B0.16x/(l+0.000lx)°50.12x

CO.!lx/(l+O.OOOlx)°$O.O8x/(l+0.0002K严