高斯烟羽模型

模型假设 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 处泄漏气体浓度的 函数为 22 bzay eexAzyxX 1 由概率统计理论可以写出方差的表达式为 2 0 0 2 2 0 0 2 2 Xdz Xdzz Xdy Xdyy z y 由假设可以写出源强的积分公式 3 uXdydzQ 式中 为泄漏气体在 y z 方向分布的标准差 单位为 y z m X x y z 为任一点处泄漏气体的浓度 单位为 kg m u 为平均风速 单 3 位为 m s Q 为源强 即泄漏速度 单位为 kg s 将 1 式代入 2 式 积分可得 4 2 2 2 1 2 1 z y b a 将 1 式和 4 式代入 3 式 积分可得 5 zy u Q xA 2 再将 4 式和 5 式代入 1 式 可得 6 2 2 2 2 22 exp 2 zyzy zy u Q zyxX 上式为无界空间连续点源扩散的高斯模型公式 然而在实际中 由于地面 的存在 烟羽的扩散是有界的 根据假设可以把地面看做一镜面 对泄漏气体 起全反射作用 并采用像源法处理 原理如图 2 2 所示 可以把任一点 p 处的 浓度看做两部分的贡献之和 一部分是不存在地面时所造成的泄漏物浓度 一 部分是由于地面反射作用增加的泄漏物浓度 该处的泄漏物浓度即相当于不存 在地面时由位于 0 0 H 的实源和位于 0 0 H 的像源在 P 点处所造成的泄 漏物浓度之和 其中 实源的贡献为 2 1 exp 2 1 exp 2 2 2 2 2 1 zyzy Hzy u Q zyxX 7 其中 像源的贡献为 8 2 1 exp 2 1 exp 2 2 2 2 2 2 zyzy Hzy u Q zyxX 则该处的实际浓度为 9 zyxXzyxXzyxX 21 由以上条件公式可得到高架连续点源扩散的高斯烟羽模型公式为 10 2 1 exp 2 1 exp 2 1 exp 2 2 2 2 2 2 2 zz yzy HzHz y u Q HtzyxX 其中 X x y z 为下风向 x 米 横向 y 米 地面上方 z 米处的扩散的气 体浓度 单位为 kg m Q 为源强 即源释放速率 单位为 kg s u 为平均风 3 速 单位为 m s 为水平扩散参数 单位为 m 为垂直扩散参数 单位为 y z m t 为泄漏后是时间 单位为 s H 为泄漏源有效高度 单位为 m y 为横向距 离 单位为 m z 为垂直方向距离 单位为 m 如 式中 令 z 0 即可得到地面气体浓度计算公式 11 2 2 2 2 exp 2 1 exp 0 zyzy Hy u Q HyxX 令 y z 0 即可得到地面轴线气体浓度计算公式 2 exp 0 0 2 2 zzy H u Q HxX 12 其中 X x y 0 为下风向 x 米 横向 y 米处的地面扩散气体浓度 单位 为 kg m 若令 y 0 则可以得到下风向中心线上的浓度分布 3 4 泄漏源有效高度 烟云抬升高度的计算 以上式中的泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状的 时候气云中心的离地高度 实际上 泄漏源有效高度就等于泄漏源几何高度加 泄漏烟云抬升高度 影响烟云抬升高度的因素有很多 主要包括 泄漏气体的初始速度和方向 初始温度 泄漏口直径 环境风速及风速岁高度的变化率 环境温度及大气稳 定度 有效源高 HHH S 其中 为泄漏源几何架高 为烟云抬升高度 s H H 实验表明 泄漏源抬升高度可以用下面公式近似计算 13 VdVH S 42 其中 是气云释放速度 单位为 m s 是泄漏出口直径 单位为 S Vd 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 Pasquill 大气稳定度的确定 白天日照夜间条件 地面风速 m s 强中等弱 阴天且云 层薄 或 低空云量 为 4 8 天空云量 为 3 8 2 3 2 4 3 6 4 6 A A B B C C A B B B C C D D B C C D D E D D D F E D D 表 2 日照强度的确定 天空云层的情况日照角 60 日照角 35 日照角 15 且 35 天空云量为 4 8 或高空有薄云 强中等弱 天空云量为 5 8 7 8 云层高度为 2134 4877m 中等弱弱 天空云量为 5 8 7 8 云层高度 2134m 弱弱弱 表 1 和表 2 中的云量是指当地天空层覆盖率 例如 云量为 3 8 是指当地 3