【精品】空气污染预报实习论文.doc

【精品】空气污染预报实习论文.doc平流扩散及随机游动方法预报污染物浓度的比较分析学生姓名:叶丹大气物理20061320005 指导教师:王咏薇摘要 本文应用平流扩散及随机游动这两种方法分别预报北京方庄小区污染物浓度分布,并用Grads绘制各自预报所得的污染物浓度分布图, 平流扩散 随即游动污染物浓度预报—、引言如今,大型工厂越来越多,环境污染也越发严重,当工厂、机车等排放的污染物进入大气层示,通过扩散稀释、转化、迁移和清除在大气层内散布开來,如何了解污染物的扩散规律对我们控制和预报污染物的变化是非常有用的。目前采用的预报方法分为统计预报和数值预报两种。而平流扩散是一种数值预报方式,是以梯度输送理论为基础求解湍流扩散方程的方法,在大气扩散理论处理和应用屮起着重要作用,当前,屮小尺度数值模式大多采用高阶闭合,而大尺度区域的扩散输送问题仍然采用K闭合。而随机游动是一种新的扩散模拟方法,又称蒙特卡洛模拟或马尔科夫链模拟。随机游动扩散模拟使统计理论的发展不必受平稳和均匀假设的限制,随着高速、人容量电子计算机的发展,这种基于大量粒了轨迹统计的数值处理的扩散模拟,有可能较为经济地实施运行,因而开始得到广泛的应用和发展。二、方法介绍1、平流扩散基木原理:将雷诺平均方程屮的湍流粘性应力项采用湍流半经验理论进行闭合,按照欧拉方法处理扩散问题。主要方法:采川梯度输送理论,即由湍流运动引起的局地质量通量与该地被扩散的平均浓度梯度成正比,方向相反。xyz三个方向的湍流扩散系数。① K、,,等dydz,将速度和浓度写成平均值和脉动值Z和,取雷诺平均示将①式代入,得普遍形式的湍流扩散方程再写成三维显式差分格式:qS+N心”+12Ax几厂2q:+q二u^-u^十R必-2q;+心2ArAv2再由初始条件得到预报时间的污染物浓度状况。=Ky=、随机游动基木原理:湍流具有高度的随机性,将随机函数和随机场理论引入湍流扩散研究,发展形成统计理论处理途径。主要方法:对随机的大气扩散行为,用大量粒了的随机游动方式来模拟,即用大量标记粒了的施放来表征污染物的连续排放,让他们在流场屮按平均风输送,同时乂用一系列随机位移来模拟大气扩散,这样就表达了平流输送和湍流扩散两种作用。这些质点在空间和时间上的总体分布,构成空气污染物的散布图。标记粒了的输运速度由平均速度和脉动速度纟fl成,即U=U+U'(1);为粒子的脉动速度。脉动速度分帚为相关分最和随机分最(又称蒙特卡洛分最)Z和,即U\t+M)=U\t)RL^t)+p(2)式屮△/为时间步长;Rg)为拉格朗口自相关系数;。是随机分量,它可以表达为p=1-/?;(人»巧(3)式屮b为速度脉动标准差;了为具有零平均值和单位标准差的随机数。对于拉格朗口白相关系数,一般取为指数形式,即/?z(Ar)=exp(—)⑷式屮7;=PR(t)dt为拉格朗口积分时间尺度。Tl式(2)〜(4)的物理意义是:把脉动扩散看成是一个连续过程,它遵从马尔可夫假设,即粒了与周I韦I连续交换动量。在确定马尔可夫扩散过稈以后,即可通过计算施放粒了的轨迹来求得某区域的污染物浓度。在模拟计算时,可将模拟区域划分为网格,则网格内污染物的瞬时浓度分布为Cijk=<筹 (5)式屮i,j,k表示网格的编号,•Q为污染物排放总量,Bq;N为释放的标记粒子总数,•刀为Z时刻统计每个网格单元的粒了数;△孑为网格的体积。在瞬时浓度的基础上还可以计算时间积分浓度。此时不但要考虑网格中的粒子数量,同时也要考虑粒了在网格屮的停留时间。计算式如下N(6)Ar式屮Tp,(i,j,k)为第门个粒了在网格",j,p=i""/VAVk)屮的停留时间。利用随机游动模型计算污染物在人气屮的扩散传输过稈是非常清晰的。其屮最核心的是通过式(2)〜(4)求得不同时刻的粒子脉动速度,而在脉动速度的求取屮,拉格朗LIH相关系数乩(△力、速度脉动标准差。以及随机数$是:3大要素,在下一-部分的模型应用过稈屮会有具体描述。在获得脉动速度示就可以跟踪粒了的运动轨迹,记录它们的随机位移,最后通过式(5)和(6)统计这些粒了在时间和空间上的总体分布,从而得出污染物的扩散规律。【2】三、研究区域以北京城区东南部的位于南二环附近的方庄小区(其屮含道路、楼房、绿化带及其它公共设施)为主要研究对象,预报x方向约1340m,y方向约1300m,z方向1550m的区域内一小时后的污染状况。以下为2002年1刀1口8点时(初始时刻)方庄小区的流场和污染情况。图1:方庄小区建筑物高度分布(单位:m)图2: