高斯烟羽模型 的参数选取和地形影响

模型假设：1、 坐标系高斯模型的坐标系如图2.1所示，原点为排放点(若为高架源，原点为排放点在地面的投影)，x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x轴的左侧，z轴垂直于水平面xoy,向上为正向。在此坐标系下烟流中心线或烟流中心线在xoy面的投影与x轴重合。图2.1高斯模式的坐标系2、模型假设污染物的浓度在y、z轴上的分布是高斯分布(正态分布)的；污染源的源强是连续且均匀的，初始时刻重金属云团内部的浓度、温度呈均匀分布；扩散过程中不考虑重金属云团内部温度的变化，忽略热传递、热对流及热辐射；泄漏重金属颗粒气体是理想气体，遵守理想气体状态方程；在水平方向，大气扩散系数呈各向同性；取x轴为平均扩散方向，整个扩散过程中地势一致且不发生较大的变化地面对重金属颗粒发生吸收或吸附作用；整个过程中，重金属颗粒云团内部不发生任何化学反应等。3、模型公式推导由正态分布假设可以导出下风向任意一点X(x,y,z)处重金属气体浓度的函数为：X(x,y,z)=A(x)e一ay2e—bz2由概率统计理论可以写出方差的表达式为：卜y2XdyfXdy0卜Xdz由假设可以写出源强的积分公式3)Q=卜卜uXdydz3)—g—g式中：b、b为泄漏气体在y、z方向分布的标准差，单位为m；X(x,y,z)yz为任一点处重金属颗粒气体的浓度，单位为kg/m3；u为平均风速，单位为m/sQ为源强(即重金属排放速度)，单位为kg/s；将(1)式代入(2)式，积分可得：1a— 2b2y1a— 2b2yb二丄2b2z将(1)式和4)4)Q式代入(3)式，积分可得Qt PA(x)=2兀ubb 2nutbb 2nSbbyz yz yzP为污染源在传播过程中的源强总量；S为传播的距离再将(4)式和(5)式代入(1)式，可得：(\y2.z25)PX(x,y,z)= exp2兀SbbyzC+ (6)I2b2 2b2I'y z丿-上式为无界空间连续点源扩散的高斯模型公式改进公式。其中：X(x,y,z)为下风向x米、横向y米、地面上方z米处的扩散的气体浓度,单位为kg/m3；Q为源强(即源释放速率)，单位为kg/s；u为平均风速，单位为m/s；t为泄漏后是时间，单位为s；H为泄漏源有效高度，单位为m；y为横向距离，单位为m；z为垂直方向距离，单位为m。如：式中，令z=0，即可得到地面气体浓度计算公式：PX(xPX(x,y,0,H)= exp兀Sddyz-丄(丄)2dyexp无界空间的烟羽模型的解(11)其中，X(x,y,O)为下风向x米、横向y米处的地面扩散重金属颗粒浓度,单位为kg/m3。5、扩散系数的选取：扩散系数d、d、d的大小与大气湍流结构、离地高度、地面粗糙度、泄xyz漏持续时间、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。大气的湍流结构和风速在大气稳定度中考虑。大气稳定度由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。表1Pasquill大气稳定度的确定地面风速(m/s)白天日照夜间条件强中等弱阴天且云层薄，或低空云量为4/8天空云量为3/8<2AA-BB2~3A-BBCEF3~4BB-CCDE4~6CC-DDDD>6CDDDD表2日照强度的确定天空云层的情况日照角>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的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表3。表3扩散系数的计算方法大气稳定度A0.22x/(]+^QQQ\x^0.2xB0」6x/(l+O.OOOU)°S0A2xCOJU/(l+0.0001x)fli0.08x/(l+0.0002^)ft5DO.D8jr/(l+0,OOOU)0i0.06x/(l+0.00l5x)°5E0.06x/(l+0.000lx)°S0.03x/d+0.0003a)F0.04a/(]+0.000[x)°JO.Qf6V(1+0.0003^)4、地势的作用在该模型的影响：以上式中的重金属颗粒云团有效高度是指泄漏气体形成的气云基本上变成水平状的时候气云中心的离地高度。在上式中假定为污染源的地势海拔。实际上，泄漏源有效高度就等于泄漏源几何高度加泄漏烟云抬升高度。影响烟云抬升高度的因素有很多，主要包括：泄漏气体的初始速度和方向、初始温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳定度。有效源高：H有效源高：H二H+AHS其中，Hs为泄漏源几何架高，AH为烟云抬升高度实验表明，泄漏源抬升高度可以用下面公式近似计算：AH=2.4Vd/VS

其中，V是气云释放速度，单位为m/s；d是泄漏出口直径，单位为m；VS为环境风速，单位为m/s；当地形高于烟囱有效高度；(x,0,h,H)=texp(x,0,h,H)=texpyz(H-h》 T——2c2Z.(12)h为地形高度，亦计算点地面高于烟囱底的高度，当h>H时，烟流中心直TT接落于地面计算点