空气污染预报实习论文.doc

平流扩散及随机游动方法预报
污染物浓度的比较分析
学生姓名:叶丹大气物理 20061320005 指导教师:王咏薇
摘要本文应用平流扩散及随机游动这两种方法分别预报北京方庄小区污染物浓度分布,并用Grads绘制各自预报所得的污染物浓度分布图,以此对这两种方法做简单的对比分析。
关键词平流扩散随即游动污染物浓度预报
引言
如今,大型工厂越来越多,环境污染也越发严重,当工厂、机车等排放的污染物进入大气层后,通过扩散稀释、转化、迁移和清除在大气层内散布开来,如何了解污染物的扩散规律对我们控制和预报污染物的变化是非常有用的。
目前采用的预报方法分为统计预报和数值预报两种。而平流扩散是一种数值预报方式,是以梯度输送理论为基础求解湍流扩散方程的方法,在大气扩散理论处理和应用中起着重要作用,当前,中小尺度数值模式大多采用高阶闭合,而大尺度区域的扩散输送问题仍然采用K闭合。
而随机游动是一种新的扩散模拟方法,又称蒙特卡洛模拟或马尔科夫链模拟。随机游动扩散模拟使统计理论的发展不必受平稳和均匀假设的限制,随着高速、大容量电子计算机的发展,这种基于大量粒子轨迹统计的数值处理的扩散模拟,有可能较为经济地实施运行,因而开始得到广泛的应用和发展。
二、方法介绍
1、平流扩散
基本原理:将雷诺平均方程中的湍流粘性应力项采用湍流半经验理论进行闭合,按照欧拉方法处理扩散问题。
主要方法:采用梯度输送理论,即由湍流运动引起的局地质量通量与该地被扩散的平均浓度梯度成正比,方向相反。
,, ①,,为x,y,z三个方向的湍流扩散系数。
再由连续方程,,将速度和浓度写成平均值和脉动值之和,取雷诺平均后将①式代入,得普遍形式的湍流扩散方程
再写成三维显式差分格式:
【1】
再由初始条件得到预报时间的污染物浓度状况。
(1)
(2) 边界条件无梯度边界条件
(3) 扩散系数的计算 Kx=Ky=Kz=×tke2/ep/
2、随机游动
基本原理:湍流具有高度的随机性,将随机函数和随机场理论引入湍流扩散研究,发展形成统计理论处理途径。
主要方法:对随机的大气扩散行为,用大量粒子的随机游动方式来模拟,即用大量标记粒子的施放来表征污染物的连续排放,让他们在流场中按平均风输送,同时又用一系列随机位移来模拟大气扩散,这样就表达了平流输送和湍流扩散两种作用。这些质点在空间和时间上的总体分布,构成空气污染物的散布图。
标记粒子的输运速度由平均速度和脉动速度组成,即
⑴式中. 为粒子的平均速度; 为粒子的脉动速度。
脉动速度分量为相关分量和随机分量(又称蒙特卡洛分量) 之和,即
⑵式中为时间步长; 为拉格朗日自相关系数; 是随机分量,它可以表达为
⑶式中为速度脉动标准差; 为具有零平均值和单位标准差的随机数。
对于拉格朗日自相关系数,一般取为指数形式,即
⑷式中为拉格朗日积分时间尺度。
式(2) ~ (4) 的物理意义是: 把脉动扩散看成是一个连续过程,它遵从马尔可夫假设,即粒子与周围连续交换动量。在确定马尔可夫扩散过程以后,即可通过计算施放粒子的轨迹来求得某区域的污染物浓度。
在模拟计算时,可将模拟区域划分为网格,则网格内污染物的瞬时浓度分布为
(5) 式中 i , j , k