数学建模高斯扩散模型.docVIP

§4-2高斯扩散模式

ū—平均风速；????????????

Q—源强是指污染物排放速率。与空气中污染物质的浓度成正比，它是研究空气污染问题的基础数据。通常：

?（ⅰ）瞬时点源的源强以一次释放的总量表示；

?（ⅱ）连续点源以单位时间的释放量表示；

?（ⅲ）连续线源以单位时间单位长度的排放量表示；

?（ⅳ）连续面源以单位时间单位面积的排放量表示。

?δy—侧向扩散参数，污染物在y方向分布的标准偏差，是距离y的函数，m；

?δz—竖向扩散参数，污染物在z方向分布的标准偏差，是距离z的函数，m；

未知量—浓度c、待定函数A(x)、待定系数a、b；

式①、②、③、④组成一方程组，四个方程式有四个未知数，故方程式可解。

二、高斯扩散模式

（一）连续点源的扩散

连续点源一般指排放大量污染物的烟囱、放散管、通风口等。排放口安置在地面的称为地面点源，处于高空位置的称为高架点源。

1.大空间点源扩散

???高斯扩散公式的建立有如下假设：①风的平均流场稳定，风速均匀，风向平直；②污染物的浓度在y、z轴方向符合正态分布；③污染物在输送扩散中质量守恒；④污染源的源强均匀、连续。

图5－9所示为点源的高斯扩散模式示意图。有效源位于坐标原点o处，平均风向与x轴平行，并与x轴正向同向。假设点源在没有任何障碍物的自由空间扩散，不考虑下垫面的存在。大气中的扩散是具有y与z两个坐标方向的二维正态分布，当两坐标方向的随机变量独立时，分布密度为每个坐标方向的一维正态分布密度函数的乘积。由正态分布的假设条件②，参照正态分布函数的基本形式式（5－15），取μ＝0，则在点源下风向任一点的浓度分布函数为：

?????????（5－16）

式中?C—空间点（x，y，z）的污染物的浓度，mg/m3；

?????A（x）—待定函数；

?????σy、σz—分别为水平、垂直方向的标准差，即y、x方向的扩散参数，m。

???由守恒和连续假设条件③和④，在任一垂直于x轴的烟流截面上有：

???????????????????????（5－17）

式中?q—源强，即单位时间内排放的污染物，μg/s；

?????u—平均风速，m/s。

???将式（5－16）代入式（5－17）,由风速稳定假设条件①，A与y、z无关，考虑到③和④，积分可得待定函数A（x）：

?????????????????????????????????????????（5－18）

将式（5－18）代入式（5－16），得大空间连续点源的高斯扩散模式

?

???????????（5－19）

式中，扩散系数σy、σz与大气稳定度和水平距离x有关，并随x的增大而增加。当y＝0，z＝0时，A（x）＝C（x，0,0），即A（x）为x轴上的浓度，也是垂直于x轴截面上污染物的最大浓度点Cmax。当x→∞，σy及σz→∞，则C→0，表明污染物以在大气中得以完全扩散。

???2．高架点源扩散

???在点源的实际扩散中，污染物可能受到地面障碍物的阻挡，因此应当考虑地面对扩散的影响。处理的方法是，或者假定污染物在扩散过程中的质量不变，到达地面时不发生沉降或化学反应而全部反射；或者污染物在没有反射而被全部吸收，实际情况应在这两者之间。

（1）高架点源扩散模式。点源在地面上的投影点o作为坐标原点，有效源位于z轴上某点，z＝H。高架有效源的高度由两部分组成，即H＝h＋Δh，其中h为排放口的有效高度，Δh是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力使烟流抬升的一个附加高度，如图5－10所示。

???当污染物到达地面后被全部反射时，可以按照全反射原理，用“像源法”来求解空间某点k的浓度。图5－10中k点的浓度显然比大空间点源扩散公式（5－19）计

高斯扩散模式的一般适用条件是：①地面开阔平坦，性质均匀，下垫面以上大气湍流稳定；②扩散处于同一大气温度层结中，扩散范围小于10km；③扩散物质随空气一起运动，在扩散输送过程中不产生化学反应，地面也不吸收污染物而全反射；④平均风向和风速平直稳定，且u＞1～2m/s。

???高斯扩散模式适应大气湍流的性质，物理概念明确，估算污染浓度的结果基本上能与实验资料相吻合，且只需利用常规气象资料即可进行简单的数学运算，因此使用最为普遍。

（二）连续线源的扩散

???当污染物沿一水平方向连续排放时，可将其视为一线源，如汽车行驶在平坦开阔的公路上。线源在横风向排放的污染物浓度相等，这样，可将点源扩散的高斯模式对变量y积分，即可获得线源的高斯扩散模式。但由于线源排放路径相对固定，具有方向性，若取平均风向为x轴，则线源与平均风向未必同向。所以线源的情况较复杂，应当考虑线源与风向夹角以及线源的长度等问题。

???如果风向和线源的夹角β＞45?，无限长连续线源下风向地面浓度分布