2023年高斯烟羽扩散模型再研究.docx

高斯烟羽扩散模型再研究 篇一：高斯烟羽模型 模型假设： 1、 坐标系 高斯模型的坐标系如图2.1所示，原点为排放点（假设为高架源，原点为排放点在地面的投影），x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x轴的左侧，z轴垂直于水平面xoy，向上为正向。在此坐标系下烟流中心线或烟流中心线在xoy面的投影与x轴重合。 2、模型假设 (1)污染物的浓度在y、z轴上的分布是高斯分布（正态分布）的； (2)污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布； (3)扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射； (4)泄漏气体是理想气体，恪守理想气体状态方程； (5)在水平方向，大气扩散系数呈各向同性； (6)取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，不随地点、时间变化而变化； (7)地面对泄漏气体起全反射作用，不发生吸收或吸附作用； (8)整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反响等。 3、模型公式推导 由正态分布假设能够导出下风向任意一点X（x,y,z）处泄漏气体浓度的函数为： X(x,y,z)A(x)eay2ebz2 （1） 由概率统计理论能够写出方差的表达式为： 2y 2z Q0y2Xdy00Xdyz2Xdz （2） 0Xdz由假设能够写出源强的积分公式： uXdydz（3） 式中：y、z为泄漏气体在y、z方向分布的标准差，单位为 m；X（x,y,z）为任一点处泄漏气体的浓度，单位为 kg/m3；u为平均风速，单位为 m/s；Q为源强（即泄漏速度），单位为 kg/s； 将（1）式代入（2）式，积分可得： 1a22yb1 22z A（x）（4） 将（1）式和（4）式代入（3）式，积分可得： Q （5） 2uyz （6） 再将（4）式和（5）式代入（1）式，可得： 2y2zX（x,y,z）exp2222uyz2yzQ 上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的存在，烟羽的扩散是有界的。依照假设能够把地面看做一镜面，对泄漏气体起全反射作用，并采纳像源法处理，原理如图2.2所示。能够把任一点p处的浓度看做两局部的奉献之和：一局部是不存在地面时所造成的泄漏物浓度；一局部是由于地面反射作用增加的泄漏物浓度。该处的泄漏物浓度即相当于不存在地面时由位于（0,0,H）的实源和位于（0,0,-H）的像源在P点处所造成的泄漏物浓度之和。 其中，实源的奉献为： Q1y21(zH)2 X1(x,y,z)exp()exp()（7） 222uyz2y2z 其中，像源的奉献为： Q1y21(zH)2 X（exp()exp()（8） 2x,y,z）222uyz2y2z 那么该处的实际浓度为： X（x,y,z）X（X（ （9） 1x,y,z）2x,y,z） 由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为： Q1y2 X（x,y,z,t,H）exp()22uyz2y 1(zH)21(zH)2 [exp()exp()]222z2z （10） 其中：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为垂直方向间隔，单位为m。 如：式中，令z=0，即可得到地面气体浓度计算公式： 1y2H2 X（x,y,0,H）exp）exp22uyz2yzQ(11) 令y=z=0，即可得到地面轴线气体浓度计算公式： H2 X(x,0,0,H)exp(2) （12）uyz2zQ 其中，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、D、E、F六类。其中A、B、C三类表示气象条件不稳定，E、F两类表示气象条件稳定，D类表示中性气象条件，也确实是说气象条件的稳定性在稳定和不稳定之间。A、B、C三品种型的稳定度中，A类表示气象条件极其不稳定，B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定。E和F两品种型的稳定度中，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定。大气稳定度详细分类方法见下表1、表2。 表1和表2中的云量是指当地天空层覆盖率。例如，云量为3/8是指当地3/8的天空有云层覆盖。日照角是指当地太阳光线与地平线之间的夹度。例如，阳光垂直照射地卖弄时的日照角为90°。 一般来说，随着大气稳定度的增加，扩散系数减小。依照Hanna和Drivas的建议，化学危险品事故泄漏扩散系数与大气稳定度类型和下风向的关系如下表3。 表3 扩散系数的计算方法 篇二：高斯烟羽模型 模型假设： 1、 坐标系 高斯模型的坐标系如图2.1所示，原点为排放点（假设为高架源，原点为排放点在地面的投影），x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x轴的左侧，z轴垂直于水平面xoy，向上为正向。在此坐标系下烟流中心线或烟流中心线在xoy面的投影与x轴重合。 2、模型假设 (1)污染物的浓度在y、z轴上的分布是高斯分布（正态分布）的； (2)污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布； (3)扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射； (4)泄漏气体是理想气体，恪守理想气体状态方程； (5)在水平方向，大气扩散系数呈各向同性； (6)取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，不随地点、时间变化而变化； (7)地面对泄漏气体起全反射作用，不发生吸收或吸附作用； (8)整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反响等。 3、模型公式推导 由正态分布假设能够导出下风向任意一点X（x,y,z）处泄漏气体浓度的函数为： X(x,y,z)A(x)e ay 2 e bz 2 （1） 由概率统计理论能够写出方差的表达式为： y22 zQ yXdy 2 Xdy 2 zXdz （2） Xdz 由假设能够写出源强的积分公式： uXdydz（3） 式中：y、z为泄漏气体在y、z方向分布的标准差，单位为 m；X（x,y,z）为任一点处泄漏气体的浓度，单位为 kg/m3；u为平均风速，单位为 m/s；Q为源强（即泄漏速度），单位为 kg/s； 将（1）式代入（2）式，积分可得： a122y b1 2 2z A（x） （4） 将（1）式和（4）式代入（3）式，积分可得： Q （5） 2uyz （6） 再将（4）式和（5）式代入（1）式，可得： 2y2 z X（x,y,z）exp2222uyz2yz Q 上式为无界空间连续点源扩散的高斯模型公式，然而在实际中，由于地面的存在，烟羽的扩散是有界的。依照假设能够把地面看做一镜面，对泄漏气体起全反射作用，并采纳像源法处理，原理如图2.2所示。能够把任一点p处的浓度看做两局部的奉献之和：一局部是不存在地面时所造成的泄漏物浓度；一局部是由于地面反射作用增加的泄漏物浓度。该处的泄漏物浓度即相当于不存在地面时由位于（0,0,H）的实源和位于（0,0,-H）的像源在P点处所造成的泄漏物浓度之和。 其中，实源的奉献为： Q1y21(zH)2 exp()exp()（7） X1(x,y,z)22 2uyz2y2z 其中，像源的奉献为： Q1y21(zH)2 X（）exp()exp()（8） 2x,y,z22 2uyz2y2z 那么该处的实际浓度为： X（x,y,z）X（x,y,z）X（x,y,z） （9） 12 由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为： X（x,y,z,t,H） Q2uyz exp( 1y22 2 y )[exp( 1(zH)22 2z )exp( 1(zH)22 2z )] （10） 其中：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为垂直方向间隔，单位为m。 如：式中，令z=0，即可得到地面气体浓度计算公式： 1y2H2 X（x,y,0,H）exp）exp22 uyz2yz Q (11) 令y=z=0，即可得到地面轴线气体浓度计算公式： H2 X(x,0,0,H)exp(2) （12） uyz2z Q 其中，X（x,y,0）为下风向x米、横向y米处的地面扩散气体浓度，单位为kg/m3；假设令y=0，那么能够得到下风向中心线上的浓度分布。 4、泄漏源有效高度（烟云抬升高度的计算）： 以上式中的泄漏源有效高度是指泄漏气体构成的气云根本上变成水平状的时候气云中心的离地高度。实际上，泄漏源有效高度就等于泄漏源几何高度加泄漏烟云抬升高度。 阻碍烟云抬升高度的要素有非常多，主要包括：泄漏气体的初始速度和方向、初始温度、泄漏口直径、环境风速及风速岁高度的变化率、环境温度及大气稳定度。 有效源高： H HSH 其中，Hs为泄漏源几何架高，H为烟云抬升高度。 实验说明，泄漏源抬升高度能够用下面公式近似计算： H2.4Vd/V （13） S 其中，VS是气云释放速度，单位为m/s；d是泄漏出口直径，单位为 m；V为环境风速，单位为m/s； （13）式是20世纪70年代末、80年代初，Wilson依照管道破裂泄漏实验所得的实验公式。实验时气体的喷射方向与风速垂直同时气体喷射途径上无障碍物。实验说明，当气体喷射方向垂直向上时，预测值与实际值之比在洗漱2以内。 计算出泄漏烟云抬升高度以后，将泄漏源抬升高度与泄漏源实际几何高度相加就得到了泄漏源有效高度。 5、扩散系数的选取： 扩散系数x、y、z的大小与大气湍流构造、离地高度、地面粗糙度、泄漏持续时间、抽样