数学建模高斯扩散模型

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
高斯扩散模式的一般适用条件是:①地面开阔平坦,性质均匀,下垫面以上大气湍 流稳定;②扩散处于同一大气温度层结中,扩散范围小于 10km;③扩散物质随空气一起 运动,在扩散输送过程中不产生化学反应,地面也不吸收污染物而全反射;④平均风向 和风速平直稳定,且 u>1~2m/s。
高斯扩散模式适应大气湍流的性质,物理概念明确,估算污染浓度的结果基本 上能与实验资料相吻合,且只需利用常规气象资料即可进行简单的数学运算,因此使用 最为普遍。
由于利用常规气象资料便能确定帕斯奎尔大气稳定度,因此 P-G 扩散曲线简便实 用。但是,P-G 扩散曲线是利用观测资料统计结合理论分析得到的,其应用具有一定的 经验性和局限性。σy 是利用风向脉动资料和有限的扩散观测资料作出的推测估计,σz 是在近距离应用了地面源在中性层结时的竖直扩散理论结果,也参照一些扩散试验资料 后的推算,而稳定和强不稳定两种情况的数据纯系推测结果。一般,P-G 扩散曲线较适 用于近地源的小尺度扩散和开阔平坦的地形。实践表明,σy 的近似估计与实际状况比 较符合,但要对地面粗糙度和取样时间进行修正;σz 的估计值与温度层结的关系很大, 适用于近地源的 lkm 以内的扩散。因此,大气扩散参数的准确定量描述仍是深入研究的 课题。
令 则当
地面最大浓度值 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 可以有效地防止污染物在地面某一局部区域的聚积。
当污染物到达地面后被全部反射时,可以按照全反射原理,用“像源法”来求 解空间某点 k 的浓度。图 5-10 中 k 点的浓度显然比大空间点源扩散公式(5-19)计
算值大,它是位于(0,0,H)的实源在 k 点扩散的浓度和反射回来的浓度的叠加。反射 浓度可视为由一与实源对称的位于(0,0,-H)的像源(假想源)扩散到 k 点的浓度。 由图可见,k 点在以实源为原点的坐标系中的垂直坐标为(z-H),则实源在 k 点扩散的浓 度为式(5-19)的坐标沿 z 轴向下平移距离 H:
式(5-25)和式(5-26)是在估算大气污染时经常选用的计算公式。由于它们是

σy/σz=常数的假定下得到的,应用于小尺度湍流扩散更合适。除了极稳
定或极不稳定的大气条件,通常可设σy/σz=2 估算最大地面浓度,其估算值与孤立高
架点源(如电厂烟囱)附近的环境监测数据比较一致。通过理论或经验的方法可得σz=f
如果排放源相对较高,而且高度相差较大,也可假定 z 方向上有一虚拟点源,由源
的最初垂直分布的标准差确定 ,再由 求出 ,由 出σy,最后代入式(5-33)求出地面浓度。
求出σz,由(x0+x) 求
第二种扩散模式假定污染物浓度均匀分布在面源的 y 方向,且扩散后的污染物 全都均匀分布在长为π(x0+x)/8 的弧上,如图 5-12 所示。因此,利用式(5-32)求 σy 后,由稳定度级别应用 P-G 曲线图查出 x0,再由(x0+x) 查出σz,则面源下风向任 一点的地面浓度由下式确定:
第一种扩散模式假定污染物排放量集中在各面源单元的形心上。由假设①可得:
32)
(5-
由确定的大气稳定度级别和上式求出的 ,应用 P-G 曲线图(见下节)可查取 xo。再由(x0+x)分布查出σy 和σz,则面源下风向任一处的地面浓度由下式确定:
wenku.baidu.com
(5-33)
上式即为点源扩散的高斯模式(5-24),式中 H 取面源的平均高度,m。
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)。
式中 C—空间点(x,y,z)的污染物的浓度,mg/m3; A(x)—待定函数;
(5-16)
σy、σz—分别为水平、垂直方向的标准差,即 y、x 方向的扩散参数,m。 由守恒和连续假设条件③和④,在任一垂直于 x 轴的烟流截面上有:
(5-17) 式中 q—源强,即单位时间内排放的污染物,μg/s;
二、高斯扩散模式
(一)连续点源的扩散 连续点源一般指排放大量污染物的烟囱、放散管、通风口等。排放口安置在地面的 称为地面点源,处于高空位置的称为高架点源。 1. 大空间点源扩散
高斯扩散公式的建立有如下假设:①风的平均流场稳定,风速均匀,风向平直; ②污染物的浓度在 y、z 轴方向符合正态分布;③污染物在输送扩散中质量守恒;④污 染源的源强均匀、连续。
(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,表明污染物以在大气
(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,则在点源 下风向任一点的浓度分布函数为:
中得以完全扩散。
2.高架点源扩散
在点源的实际扩散中,污染物可能受到地面障碍物的阻挡,因此应 当考虑地面对扩散的影响。处理的方法是,或者假定污染物在扩散过程中 的质量不变,到达地面时不发生沉降或化学反应而全部反射;或者污染物 在没有反射而被全部吸收,实际情况应在这两者之间。
(1)高架点源扩散模式。点源在地面上的投影点 o 作为坐标原点,有 效源位于 z 轴上某点, z=H。高架有效源的高度由两部分组成,即 H=h+Δh,其中 h 为排放口的有效高度,Δh 是热烟流的浮升力和烟气以一定速度竖直离开排放口的冲力 使烟流抬升的一个附加高度,如图 5-10 所示。
(二)连续线源的扩散
当污染物沿一水平方向连续排放时,可将其视为一线源,如汽车行驶在平坦开 阔的公路上。线源在横风向排放的污染物浓度相等,这样,可将点源扩散的高斯模式对 变量 y 积分,即可获得线源的高斯扩散模式。但由于线源排放路径相对固定,具有方向 性,若取平均风向为 x 轴,则线源与平均风向未必同向。所以线源的情况较复杂,应当 考虑线源与风向夹角以及线源的长度等问题。
(5 -27)
其浓度是大空间连续点源扩散式(5-19)或地面无反射高架点源扩散式(5-20)在 H =0 时的两倍,说明烟流的下半部分完全对称反射到上部分,使得浓度加倍。若取 y 与 z 等于零,则可得到沿 x 轴线上的浓度分布:
(5-28)
如果污染物到达地面后被完全吸收,其浓度即为地面无反射高架点源扩散式(5- 20)在 H=0 时的浓度,也即大空间连续点源扩散式(5-19)。
如果风向和线源的夹角β>45,无限长连续线源下风向地面浓度分布为:
29)
(5-
当β<45时,以上模式不能应用。如果风向和线源的夹角垂直,即β=90,可得:
(5 -30)
对于有限长的线源,线源末端引起的“边缘效应”将对污染物的浓度分布有很 大影响。随着污染物接受点距线源的距离增加,“边源效应”将在横风向距离的更远处 起作用。因此在估算有限长污染源形成的浓度分布时,“边源效应”不能忽视。对于横 风向的有限长线源,应以污染物接受点的平均风向为 x 轴。若线源的范围是从 y1 到 y2, 且 y1<y2,则有限长线源地面浓度分布为:
(x)的具体表达式,代入(5-25)可求出最大浓度点离源的距离 xmax,具体可查阅我
国 GB3840—91《制定地方大气污染物排放标准的技术方法》。
3. 地面点源扩散
对于地面点源,则有效源高度 H=0。当污染物到达地面后被全部反射时,可令式(5 -22)中 H=0,即得出地面连续点源的高斯扩散公式:
§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;②面源单元 向下风向扩散的浓度可用虚拟点源在下风向造成的同样的浓度所代替。根据污染物在面 源范围内的分布状况,可分为以下两种虚拟点源扩散模式:
相关文档
最新文档