数学建模高斯扩散模型

二、高斯扩散模式
（一）连续点源的扩散 连续点源一般指排放大量污染物的烟囱、放散管、通风口等。排放口安置在地面的 称为地面点源，处于高空位置的称为高架点源。 1. 大空间点源扩散
高斯扩散公式的建立有如下假设：①风的平均流场稳定，风速均匀，风向平直； ②污染物的浓度在 y、z 轴方向符合正态分布；③污染物在输送扩散中质量守恒；④污 染源的源强均匀、连续。
20）
（5－
k 点在以像源为原点的坐标系中的垂直坐标为(z＋H)，则像源在 k 点扩散的浓度为 式（5－19）的坐标沿 z 轴向上平移距离 H：
21） 由此，实源 Cs 与像源 Cx 之和即为 k 点的实际污染物浓度：
（5－
（5－22）
若污染物到达地面后被完全吸收，则 Cx＝0，污染物浓度 C（x，y，z，H）＝Cs， 即式（5－20）。
§4-2 高斯扩散模式
ū — 平均风速； Q—源强是指污染物排放速率。与空气中污染物质的浓度成正比，它是研究空气 污染问题的基础数据。通常：
（ⅰ）瞬时点源的源强以一次释放的总量表示； （ⅱ）连续点源以单位时间的释放量表示； （ⅲ）连续线源以单位时间单位长度的排放量表示； （ⅳ）连续面源以单位时间单位面积的排放量表示。 δy—侧向扩散参数，污染物在 y 方向分布的标准偏差，是距离 y 的函数，m； δz—竖向扩散参数，污染物在 z 方向分布的标准偏差，是距离 z 的函数，m； 未知量—浓度 c、待定函数 A(x)、待定系数 a、b； 式①、②、③、④组成一方程组，四个方程式有四个未知数，故方程式可解。
扩散参数σy、σz 是表示扩散范围及速率大小的特征量，也即正态分布函数的标准 差。为了能较符合实际地确定这些扩散参数，许多研究工作致力于把浓度场和气象条件 结合起来，提出了各种符合实验条件的扩散参数估计方法。其中应用较多的由是帕斯奎 尔(Pasquill) 和吉福特(Gifford)提出的扩散参数估算方法，也称为 P－G 扩散曲线， 如图 5－13 和图 5－14 所示。由图可见，只要利用当地常规气象观测资料，由表 5－1 查取帕斯奎尔大气稳定度等级，即可确定扩散参数。扩散参数σ具有如下规律：①σ随 着离源距离增加而增大；②不稳定大气状态时的σ值大于稳定大气状态，因此大气湍流 运动愈强，σ值愈大；③以上两种条件相同时，粗糙地面上的σ值大于平坦地面。
式中，s1＝y1/σy，s2＝y2/σy，积分值可从正态概率表中查出。 （三）连续面源的扩散
（5－31）
当众多的污染源在一地区内排放时，如城市中家庭炉灶的排放， 可将它们作为面源来处理。因为这些污染源排放量很小但数量很大，若 依点源来处理，将是非常繁杂的计算工作。
常用的面源扩散模式为虚拟点源法，即将城市按污染源的分布和高低不同划分为若 干个正方形，每一正方形视为一个面源单元，边长一般在 0.5～10km 之间选取。这种方 法假设：①有一距离为 x0 的虚拟点源位于面源单元形心的上风处，如图 5-12 所示，它 在面源单元中心线处产生的烟流宽度为 2y0＝4.3σy0，等于面源单元宽度 B；②面源单元 向下风向扩散的浓度可用虚拟点源在下风向造成的同样的浓度所代替。根据污染物在面 源范围内的分布状况，可分为以下两种虚拟点源扩散模式：
如果排放源相对较高，而且高度相差较大，也可假定 z 方向上有一虚拟点源，由源
的最初垂直分布的标准差确定 ，再由 求出 ，由 出σy，最后代入式（5－33）求出地面浓度。
求出σz，由(x0＋x) 求
第二种扩散模式假定污染物浓度均匀分布在面源的 y 方向，且扩散后的污染物 全都均匀分布在长为π(x0＋x)／8 的弧上，如图 5-12 所示。因此，利用式（5－32）求 σy 后，由稳定度级别应用 P－G 曲线图查出 x0，再由(x0＋x) 查出σz，则面源下风向任 一点的地面浓度由下式确定：
（5 －27）
其浓度是大空间连续点源扩散式（5－19）或地面无反射高架点源扩散式（5－20）在 H ＝0 时的两倍，说明烟流的下半部分完全对称反射到上部分，使得浓度加倍。若取 y 与 z 等于零，则可得到沿 x 轴线上的浓度分布：
（5－28）
如果污染物到达地面后被完全吸收，其浓度即为地面无反射高架点源扩散式（5－ 20）在 H＝0 时的浓度，也即大空间连续点源扩散式（5－19）。
中得以完全扩散。
2．高架点源扩散
在点源的实际扩散中，污染物可能受到地面障碍物的阻挡，因此应 当考虑地面对扩散的影响。处理的方法是，或者假定污染物在扩散过程中 的质量不变，到达地面时不发生沉降或化学反应而全部反射；或者污染物 在没有反射而被全部吸收，实际情况应在这两者之间。
（1）高架点源扩散模式。点源在地面上的投影点 o 作为坐标原点，有 效源位于 z 轴上某点， z＝H。高架有效源的高度由两部分组成，即 H＝h＋Δh，其中 h 为排放口的有效高度，Δh 是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力 使烟流抬升的一个附加高度，如图 5－10 所示。
（二）连续线源的扩散
当污染物沿一水平方向连续排放时，可将其视为一线源，如汽车行驶在平坦开 阔的公路上。线源在横风向排放的污染物浓度相等，这样，可将点源扩散的高斯模式对 变量 y 积分，即可获得线源的高斯扩散模式。但由于线源排放路径相对固定，具有方向 性，若取平均风向为 x 轴，则线源与平均风向未必同向。所以线源的情况较复杂，应当 考虑线源与风向夹角以及线源的长度等问题。
令 则当
地面最大浓度值 Cmax 及其离源的距离 xmax 可以由式（5－24）求导并取极值得到。 ，由于σy、σz 均为 x 的未知函数，最简单的情况可假定σy/σz＝常数，
25） 时，得地面浓度最大值
（5－
26）
（5－
由式（5－25）可以看出，有效源 H 越高， xmax 处的σz 值越大，而σz∝xmax，则 Cmax 出现的位置离污染源的距离越远。式（5－26）表明，地面上最大浓度 Cmax 与有效源高度 的平方及平均风速成反比，增加 H 可以有效地防止污染物在地面某一局部区域的聚积。
高斯扩散模式的一般适用条件是：①地面开阔平坦，性质均匀，下垫面以上大气湍 流稳定；②扩散处于同一大气温度层结中，扩散范围小于 10km；③扩散物质随空气一起 运动，在扩散输送过程中不产生化学反应，地面也不吸收污染物而全反射；④平均风向 和风速平直稳定，且 u＞1～2m/s。
高斯扩散模式适应大气湍流的性质，物理概念明确，估算污染浓度的结果基本 上能与实验资料相吻合，且只需利用常规气象资料即可进行简单的数学运算，因此使用 最为普遍。
（x）的具体表达式，代入（5－25）可求出最大浓度点离源的距离 xmax，具体可查阅我
国 GB3840—91《制定地方大气污染物排放标准的技术方法》。
3. 地面点源扩散
对于地面点源，则有效源高度 H＝0来自当污染物到达地面后被全部反射时，可令式（5 －22）中 H＝0，即得出地面连续点源的高斯扩散公式:
如果风向和线源的夹角β＞45，无限长连续线源下风向地面浓度分布为：
29）
（5－
当β＜45时，以上模式不能应用。如果风向和线源的夹角垂直，即β＝90，可得：
（5 －30）
对于有限长的线源，线源末端引起的“边缘效应”将对污染物的浓度分布有很 大影响。随着污染物接受点距线源的距离增加，“边源效应”将在横风向距离的更远处 起作用。因此在估算有限长污染源形成的浓度分布时，“边源效应”不能忽视。对于横 风向的有限长线源，应以污染物接受点的平均风向为 x 轴。若线源的范围是从 y1 到 y2， 且 y1＜y2，则有限长线源地面浓度分布为：
（2）地面全部反射时的地面浓度。实际中，高架点源扩散问题中最关心的是 地面浓度的分布状况，尤其是地面最大浓度值和它离源头的距离。在式（5－22）中， 令 z＝0，可得高架点源的地面浓度公式：
（5－23） 上式中进一步令 y＝0 则可得到沿 x 轴线上的浓度分布：
（5－24） 地面浓度分布如图图 5－11 所示。y 方向的浓度以 x 轴为对称轴按正态 分布；沿 x 轴线上，在污染物排放源附近地面浓度接近于零，然后顺风向不断 增大，在离源一定距离时的某处，地面轴线上的浓度达到最大值，以后又逐渐减小。
图 5－9 所示为点源的高斯扩散模式示意图。有效源位于坐标原点 o 处，平均风向 与 x 轴平行，并与 x 轴正向同向。假设点源在没有任何障碍物的自由空间扩散，不考虑 下垫面的存在。大气中的扩散是具有 y 与 z 两个坐标方向的二维正态分布，当两坐标方 向的随机变量独立时，分布密度为每个坐标方向的一维正态分布密度函数的乘积。由正 态分布的假设条件②，参照正态分布函数的基本形式式（5－15），取μ＝0，则在点源 下风向任一点的浓度分布函数为：
第一种扩散模式假定污染物排放量集中在各面源单元的形心上。由假设①可得：
32）
（5－
由确定的大气稳定度级别和上式求出的 ，应用 P－G 曲线图（见下节）可查取 xo。再由(x0＋x)分布查出σy 和σz，则面源下风向任一处的地面浓度由下式确定：
（5－33）
上式即为点源扩散的高斯模式（5－24），式中 H 取面源的平均高度，m。
（5 －34）
三、扩散参数及烟流抬升高度的确定 高斯扩散公式的应用效果依赖于公式中的各个参数的准确程度，尤其是扩散参数σ y、σz 及烟流抬升高度Δh 的估算。其中，平均风速 u 取多年观测的常规气象数据；源 强 q 可以计算或测定，而σy、σz 及Δh 与气象条件和地面状况密切相关。
1. 扩散参数σy、σz 的估算
u—平均风速，m/s。 将式（5－16）代入式（5－17）, 由风速稳定假设条件①，A 与 y、z 无关，考
虑到
③和④，积分可得待定函数 A（x）：
（5－18） 将式（5－18）代入式（5－16），得大空间连续点源的高斯扩散模式
19）
（5－
式中，扩散系数σy、σz 与大气稳定度和水平距离 x 有关，并随 x 的增大而增加。当 y ＝0，z＝0 时，A（x）＝C（x，0,0），即 A（x）为 x 轴上的浓度，也是垂直于 x 轴截 面上污染物的最大浓度点 Cmax。当 x→∞，σy 及σz→∞，则 C→0，表明污染物以在大气