高斯烟羽模型.doc

模型假设:1、坐标系高斯模型的坐标系如图2.1所示，原点为排放点(若为高架源，原点为排放点在地面的投影)，x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x轴的左侧，z轴垂直于水平面xoy，向上为正向在此坐标系下烟流中心线或烟流中心线在xoy面的投影与x轴重合图2」為斯模式的坐标系2、模型假设(1) 污染物的浓度在y、z轴上的分布是高斯分布(正态分布)的；(2) 污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布；(3) 扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射;(4) 泄漏气体是理想气体，遵守理想气体状态方程；(5) 在水平方向，大气扩散系数呈各向同性;(6) 取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变,不随地点、时间变化而变化；(7) 地面对泄漏气体起全反射作用，不发生吸收或吸附作用；(8) 整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等3、模型公式推导由正态分布假设可以导出下风向任意一点X(x,y,z)处泄漏气体浓度的函数为:2b2X(x,y,z)A(x)eayebz(1)由概率统计理论可以写出方差的表达式为:y2Xdyo(2)(3)oXdyz2XdzoXdzo由假设可以写出源强的积分公式:uXdydz为任一点处泄漏气体的浓度，单位为kg/m3；u为平均风速，单位为m/s；Q为源强(即泄漏速度)，单位为kg/s；式中：y、z为泄漏气体在y、z方向分布的标准差，单位为m；X(x,y,z)将(1)式代入(2)式，积分可得:将（1）式和（4）式代入（3）式，积分可得:(5)A(x)再将（4）式和（5）式代入（1）式，可得:X(x,y,z)z2—expz(6)上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的存在，烟羽的扩散是有界的。

根据假设可以把地面看做一镜面，对泄漏气体起全反射作用，并采用像源法处理，原理如图2.2所示可以把任一点p处的浓度看做两部分的贡献之和：一部分是不存在地面时所造成的泄漏物浓度；一部分是由于地面反射作用增加的泄漏物浓度该处的泄漏物浓度即相当于不存在地面时由位于（0,0,H）的实源和位于（0,0,-H）的像源在P点处所造成的泄漏物浓度之俱資主输＜实遒»俱資主输＜实遒»他厨）斥84二F匕-z)^2.2像源法原理示意图其中，实源的贡献为:QXi")百-exp(z2-iy、/■)exp(-y1(zH)2)2*)(7)其中，像源的贡献为:-exp(z半)exp(2y(8)则该处的实际浓度为:X(x,y,z)X/x,y,z)X(x,y,z)(9)由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为:Q1y2X(x,y,z,t,H)exp(2)(10)(10)2uyz22[exp(27)exp(舟j)]222其中：X（x,y,z）为下风向x米、横向y米、地面上方z米处的扩散的气体浓度，单位为kg/m3；Q为源强（即源释放速率），单位为kg/s;u为平均风速，单位为m/s；y为水平扩散参数，单位为mz为垂直扩散参数，单位为m；t为泄漏后是时间，单位为s；h为泄漏源有效高度，单位为my为横向距离，单位为mz为垂直方向距离，单位为m如：式中，令z=0，即可得到地面气体浓度计算公式:X(x,y,0,H)—expzexpH2(11)令y=z=0，即可得到地面轴线气体浓度计算公式:X(x,0,0,H)-exp(zH2(12)其中，X（x,y,0）为下风向x米、横向y米处的地面扩散气体浓度，单位为kg/m3；若令y=0,则可以得到下风向中心线上的浓度分布。

4、泄漏源有效高度（烟云抬升高度的计算）：以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的时候气云中心的离地高度实际上，泄漏源有效高度就等于泄漏源几何高度加泄漏烟云抬升高度影响烟云抬升高度的因素有很多，主要包括：泄漏气体的初始速度和方向、初始温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳定度有效源高：HHsH其中，Hs为泄漏源几何架高，H为烟云抬升高度实验表明，泄漏源抬升高度可以用下面公式近似计算：H2.4VSd/V（13）其中，Vs是气云释放速度，单位为m/s;d是泄漏出口直径，单位为m；V为环境风速，单位为m/s；（13）式是20世纪70年代末、80年代初，Wilson根据管道破裂泄漏实验所得的实验公式实验时气体的喷射方向与风速垂直并且气体喷射路径上无障碍物实验表明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱2以内计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度相加就得到了泄漏源有效高度5、扩散系数的选取：扩散系数x、y、z的大小与大气湍流结构、离地高度、地面粗糙度、泄漏持续时间、抽样时间间隔、风速以及离开泄漏源的距离等因素有关。

大气的湍流结构和风速在大气稳定度中考虑大气稳定度由10米高度以上的风速、白天的太阳辐射或夜间的云量等参数决定按照Pasquill的分类方法，随着气象条件稳定性的增加，大气稳定度可以分为A、B、C、DE、F六类其中A、BC三类表示气象条件不稳定，E、F两类表示气象条件稳定，D类表示中性气象条件，也就是说气象条件的稳定性在稳定和不稳定之间AB、C三种类型的稳定度中，A类表示气象条件极其不稳定，B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定E和F两种类型的稳定度中，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定大气稳定度具体分类方法见下表1、表2表1Pasquill大气稳定度的确定地面风速（m/s）白天日照夜间条件强中等弱阴天且云层薄，或低空云量为4/8天空云量为3/82AA-BB2~3A-BBCEF3~4BB-CCDE4~6CC-DDDD6CDDDD表2日照强度的确定角召"O日八6日>O53角>照且日O60VO53角v照且日O5>44云4薄为协量空云高空或天强等中弱、咼7为层871里云I云，13空兀2天肓为等中弱弱高为层m1里云34云，21空兀匸弱弱弱表1和表2中的云量是指当地天空层覆盖率。

例如，云量为3/8是指当地3/8的天空有云层覆盖日照角是指当地太阳光线与地平线之间的夹度例如,一般来说，随着大气稳定度的增加，扩散系数减小根据Hanna和Drivas的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表3表3扩散系数的计算方法大气稳定度AB0.16jc/(UO,(XJ01x)^amCOJIi/ci+OOflOIrf10.08x/(l+0.0002D1o,oaV(i+氏two"尸606#口“,001$工严E0.06a/(1+Q.ooouro.oix/(i+afloojur)FOr04x/(l+G.000lx)fl<0.016jr/(l+fl0003i)。