高斯烟羽模型

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

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

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

4、；一部分是由于地面反射作用增加的泄漏物浓度。该处的泄漏物浓度即相当于不存在地面时由 位于（0,0,H ）的实源和位于（0,0,-H ）的像源在P点处所造成的泄漏物浓度之 和。其中，实源的贡献为：Xjx, y,z)exp(爲)expe(z-)2、 2(7)其中，像源的贡献为：2Q1 y2exp（-;-y）exp（- 2X2(x, y, z)二1 (z H)23 y；z则该处的实际浓度为：X（ x,y,z） X （x,y,z） X 2（x,y,z）(8)(9)由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为：X(x, y,z,t, H )=Q1 y22廿，（三尹）(10)y zy1 （Z

5、H ）2 丄1 （z+ H ）2exp（2） exp（2）2 2 crzz其中：X（ x,y,z ）为下风向x米、横向y米、地面上方z米处的扩散的气体 浓度，单位为kg/m3 ； Q为源强（即源释放速率），单位为kg/s ； u为平均风速,单位为m/s；二y为水平扩散参数，单位为 m二z为垂直扩散参数，单位为 m； t为泄漏后是时间，单位为s； h为泄漏源有效高度，单位为 m y为横向距离，单 位为m z为垂直方向距离，单位为 mb女口：式中，令z=0，即可得到地面气体浓度计算公式：X (x, y,0,H)二Q二 U；y；z(11)令y=z=0，即可得到地面轴线气体浓度计算公式:QHX(x,0

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

7、明，泄漏源抬升高度可以用下面公式近似计算：H = 2.4Vsd/V( 13)其中，Vs是气云释放速度，单位为 m/s； d是泄漏出口直径，单位为 m； V 为环境风速，单位为m/s；(13)式是20世纪70年代末、80年代初，Wilson根据管道破裂泄漏实验 所得的实验公式。实验时气体的喷射方向与风速垂直并且气体喷射路径上无障碍 物。实验表明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱2以内。计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度相 加就得到了泄漏源有效高度。5、扩散系数的选取：扩散系数匚X、二y、二z的大小与大气湍流结构、离地高度、地面粗糙度、泄漏持续时间、

8、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。大气的湍流结构和风速在大气稳定度中考虑。大气稳定度由10米高度以上的风速、白天的太阳辐射或夜间的云量等参数决定。按照Pasquill的分类方法，随着气象条件稳定性的增加，大气稳定度可以分为A、B、C、DE、F六类。其中A、BC三类表示气象条件不稳定，E、F两类表示气象条件稳定，D类表示中性气象条件，也就是说气象条件的稳定性在稳 定和不稳定之间。A B、C三种类型的稳定度中，A类表示气象条件极其不稳定， B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定。E和F两种类 型的稳定度中，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定。大 气稳定度具体分类方法见下表 1表2。表1 Pasquill大气稳定度的确定地面风速（m/s）白天日照夜间条件强中等弱阴天且云 层薄，或 低空云量为4/8天空云量为3/8AA-BBEFA-BBCDEBB-CCDDCC-DDDDCDD表2日照强度的确定天空云层的情况日照角 60 日照角 35 日照角 15。且 35 天空云量为4/8， 或咼空有薄云强中等弱天空云量为5/8-7/8，云层咼度为 2134-4877m中等弱弱天空云量为5/8-7/8，云层咼 度 2134m弱厶At *日.冃丄、【才rrL十1=弱弱表1和表2中的云量是指当地天空层覆盖率。例如，云量为3/8是指当