《高斯烟羽模型》由会员分享，可在线阅读，更多相关《高斯烟羽模型（7页珍藏版）》请在金锄头文库上搜索。
1、高斯烟羽模型Prepared on 22 November 2020LYOl模型假设：1、坐标系 高斯模型的坐标系如图所示，原点为排放点(若为高架源，原点为排放点在地面 的投影)，X轴正向为风速方向，y轴在水平面上垂直于X轴，正向在x轴的左侧，z 轴垂直于水平面xoy,向上为正向。在此坐标系下烟流中心线或烟流中心线在xoy面 的投影与x轴重合。2、模型假设(1) 污染物的浓度在y、z轴上的分布是高斯分布(正态分布)的；(2) 污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分 布；(3) 扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射；(4) 泄漏气体是理想气体，
2、遵守理想气体状态方程；(5) 在水平方向，大气扩散系数呈各向同性；(6) 取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，不随 地点、时间变化而变化；(7) 地面对泄漏气体起全反射作用，不发生吸收或吸附作用；(8) 整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。3、模型公式推导由正态分布假设可以导出下风向任意一点X (x,y,z)处泄漏气体浓度的函数为：X (x, y, z) = A( x)e - ay2 e - bz21)由概率统计理论可以写出方差的表达式为：y 2 Xdy t Xdy02)Q2由假设可以写出源强的积分公式：Q J 卜 uXdydz-X -X式中：
3、Q、Q为泄漏气体在y、Z方向分布的标准差，单位为m；x (x,y,z )为 yz任一点处泄漏气体的浓度，单位为 kg/m3； u 为平均风速，单位为 m/s； Q 为源强即泄漏速度)，单位为 kg/s；将(1)式代入(2)式，积分可得：12q 2y4)12q 2 z将(1)式和(4)式代入(3)式，积分可得：(x) = Q-2兀UQ Qyz再将(4)式和(5)式代入(1)式，可得：X (x, y, z)=Q2兀uo oexpyz(y22o 2y、 z2 2o 2 Iz上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的存 在，烟羽的扩散是有界的。根据假设可以把地面看做一镜面，对泄
4、漏气体起全反射作 用，并采用像源法处理，原理如图所示。可以把任一点 p 处的浓度看做两部分的贡献 之和：一部分是不存在地面时所造成的泄漏物浓度；一部分是由于地面反射作用增加 的泄漏物浓度。该处的泄漏物浓度即相当于不存在地面时由位于（0,0,H）的实源和 位于（0,0,-H）的像源在P点处所造成的泄漏物浓度之和。其中，实源的贡献为：Q /I y 21 (z - H )27)X (x, y, z) =exp(-) exp(-)i2兀uo o2 o 2 F 2 o 2y zyz其中，像源的贡献为：X (x, y,z) = -一Qexp(-2兰)exp(-2(Z + H)22兀uo o2 o 22 o
5、 2y zyz则该处的实际浓度为：X (x, y, z) = X (x, y, z) + X (x, y, z)12由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为：Q1 y 2X (x, y, z, t, H )=exp( -) x2兀 u o o2 o 210)y zy( 1 (z - H)2( 1 (z + H)2o2zexp( -) + exp(-)2 o 22 o 2z其中：X （x,y,z）为下风向x米、横向y米、地面上方z米处的扩散的气体浓 度，单位为kg/m3 ； Q为源强（即源释放速率），单位为kg/s ； u为平均风速，单位为m/s ； b为水平扩散参数，单位为m
6、； b为垂直扩散参数，单位为m ； t为泄漏后 yz是时间，单位为 s； H 为泄漏源有效高度，单位为 m； y 为横向距离，单位为 m； z 为 垂直方向距离，单位为 m。如：式中，令z=0,即可得到地面气体浓度计算公式:X (x,y,O, H) = 一Q兀ub b yzexp - )2yexp12b 2 丿z(11)令y=z=O,即可得到地面轴线气体浓度计算公式:QH 2X (x,0,0, H) =Qexp(-)(12)兀ub b2b 2y zz其中，X (x,y,O)为下风向x米、横向y米处的地面扩散气体浓度，单位为kg/m3 ；若令y=0，则可以得到下风向中心线上的浓度分布4、泄漏源有
7、效高度(烟云抬升高度的计算)：以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的时候气云中心的离地高度。实际上，泄漏源有效高度就等于泄漏源几何高度加泄漏烟云抬升 高度。高度。影响烟云抬升高度的因素有很多，主要包括：泄漏气体的初始速度和方向、初始 温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳定度。有效源高：H = HS +AH其中，H s为泄漏源几何架高，AH为烟云抬升高度。实验表明，泄漏源抬升高度可以用下面公式近似计算：AH = 2.4V d /V(13)S( 13 )其中，V是气云释放速度，单位为m/s ； d是泄漏出口直径，单位为m；V为环S境风速，单位为
8、 m/s；(13)式是20世纪70年代末、80年代初，Wilson根据管道破裂泄漏实验所得的 实验公式。实验时气体的喷射方向与风速垂直并且气体喷射路径上无障碍物。实验表 明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱 2以内。计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度相加就 得到了泄漏源有效高度。5、扩散系数的选取：扩散系数&、b、b的大小与大气湍流结构、离地高度、地面粗糙度、泄漏持 xyz续时间、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。大气的湍流结构和 风速在大气稳定度中考虑。大气稳定度由 10米高度以上的风速、白天的太阳辐射或 夜间的云量等参数决定。
9、按照 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 大气稳定度的确定地面风速(m/s)白天日照夜间条件强中等弱阴天且云 层薄，或 低空云量为
10、4/8天空云量为3/8AA-BBEFA-BBCDEBB-CCDDCC-DDDDCDD表 2 日照强度的确定天空云层的情况日照角 60日照角 35日照角 15 且 35天仝云量为4/8, 或咼空有薄云强中等弱天空云量为5/8 7/8,云层高度为 2134-4877m中等弱弱天空云量为5/87/8,云层高度2134m弱弱弱表1 和表2 中的云量是指当地天空层覆盖率。例如，云量为 3/8 是指当地 3/8 的天空有云层覆盖。日照角是指当地太阳光线与地平线之间的夹度。例如，阳光垂直照 射地卖弄时的日照角为 90。一般来说，随着大气稳定度的增加，扩散系数减小。根据 Hanna 和 Drivas 的建 议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表 3。表 3 扩散系数的计算方法
