高斯扩散模式在瞬间排放空气污染物模拟中的应用

1、高斯扩散模式在瞬间排放空气污染物模拟中的应用摘要：在文章中提出应用高斯模式模拟和预测在瞬间排放状况下空气污染等级，用FORTRAN语言编写的高斯模式程序还可应用于区域污染影响评价中，模式不仅可以从GIS中输入数据而且还可以应用GIS格式输出结果。关键词：高斯模式 空气污染 地理信息系统 瞬时污染源 浓度场瞬间排放是指工业企业或电厂的事故性污染物排放，如贮油罐或输油管道发生事故等。排放的污染物污染了空气、土壤、地面及地下水，影响植被和影响环境。模拟瞬间空气污染要求得到污染区域面积、污染浓度和等级、污染预测等。本文提出用高斯模式的解析解来模拟和预测瞬间排放空气污染状况。基于烟羽扩散上的解析公式求解

2、的高斯模式非常广泛的应用于评价区域污染状况。高斯数学模式作为一个污染物扩散的基础模式被国际原子能机构广泛推广。从瞬间点源污染源排放的污染物，其转换和扩散可以用以下的扩散方程来表示：+div(CV)=(KC)+Ri +Q(t t0)(x x0)(y y0)(z z0) (1)式中：C(x, y, z, t)为污染物浓度 V为风速 K为扩散系数 R为污染物光化学转化率 Q为污染物排放量 x0, y0, z0为污染源相对坐标在一定的风速和扩散系数条件下，公式（1）有其高斯扩散模式的解析解。因此，污染物浓度值C由点源污染源的高度H决定。H在高斯扩散模式中由下述公式计算：C（x,y,z,t）= (2)式

3、中：t为时间 Q为排放量 u，v，w为风速分别在x，y，z方向的分量 x, y, z分别在x，y，z方向的扩散系数 h为点源高度 H为混合层高度高斯模式中，假设X轴与风向方向一致，Z轴铅直向上，V=W=0。公式（2）可以转化为以下形式。C（x,y,z,t）= (3)从公式(3)我们可以看出，每一个烟团需要用不同的坐标系进行计算，当我们计算多源的污染浓度时，我们需要用到几个坐标系，这样计算起来很复杂。因此，公式（3）必须做相应的转化到同一个坐标系中。我们建立一个相对的坐标系，I表示原点，坐标轴为和（见1）。并以I为原点建立第二个坐标系，LX表示X轴，其方向与风向一致，LX与OX的夹角为。这样我们

4、可以将XIY坐标系转换为XOY坐标系，从而得到图中M点的坐标。 M y Y X O x 图1：坐标轴的转换X = x cos + y sin = (x ) cos + (y ) sin ,Y = x sin + y cos = (x ) sin + (y ) cos 我们仅对Y坐标轴进行考虑，等于360 ，我们得到最终的公式：Y = (x ) sin + (y ) cos (4)我们将（4）式中的sin和cos用风速u和v替换后，得到点源在特定坐标系中的污染物浓度计算公式：C（x,y,z,t）= (5)式中：t表示时间 Q表示源强 u, v表示风速 x, y, z表示表示X，Y，Z方向的扩散参

5、数 H表示混合层高度 x0, y0, z0表示坐标系原点 h有效烟囱高度（如，烟气抬升后的高度）在高斯模式中，y和z是重要的参数，可以用帕斯奎尔-吉福德扩散参数曲线来确定，也可将P-G扩散曲线用近似幂函数式表示，这样方便用电子计算机计算大气污染浓度分布。该幂函数的确定与大气稳定度的划分有关。在高斯模式中，城市和开阔平原的布里吉斯扩散参数以烟云的宽度和高度来表示横风向和垂直向的扩散，y 和z的计算方法见表1。大气稳定度的划分方法见表2。表1 布里吉斯扩散参数稳定度x, y (m)z (m)开阔平原A0.22x(1 + 0.0001x)1/20.2xB0.16x(1 + 0.0001x)1/20.

6、12xC0.11x(1 + 0.0001x)1/20.08x(1 + 0.0002x)1/2D0.08x(1 + 0.0001x)1/20.06x(1 + 0.0015x)1/2E0.06x(1 + 0.0001x)1/20.03x(1 + 0.0003x)1F0.04x(1 + 0.0001x)1/20.016x(1 + 0.0003x)1城市A-B0.32x(1 + 0.0004x)1/20.24x(1 + 0.001x)C0.22x(1 + 0.0004x)1/20.2xD0.16x(1 + 0.0004x)1/20.14x(1 + 0.0003x)1/2E-F0.11x(1 + 0.0

7、004x)1/20.08x(1 + 0.0015x)1/2表中X为计算点在X方向轴上的距离。表2 帕斯奎尔大气稳定度划分10米高度上的风速日间，太阳辐射强度夜间，云量强中等弱有云无云6CDDDD烟气抬升高度的计算理论已经非常广泛，通过直观的观测和理论研究，可以得到。烟气的抬升首先是由烟气的初始动量和浮力及环境空气的湍流引起，在中性和不稳定的稳定度下它能引起烟气有很高的抬升高度，烟气的抬升相当于加高了烟囱的高度，这个高度也叫做有效烟囱高度。有三个参数决定气体从污染源的排放状况，它们分别是：烟囱的物理结构（如出口直径）、烟气排放状况以及污染物的化学性质。现在研究出很多计算烟气有效高度的公式。在稳定

8、空气环境中计算烟气抬升高度现在使用得较广泛的是布里吉斯公式烟气抬升公式。 （6）式中：h表示烟气抬升高度 h表示烟囱高度 表示烟囱出口处风速 F表示浮力通量 S表示混合层参数F按下式计算： （7）式中：g表示重力加速度 Ts表示排放温度Ta表示环境温度Vs表示烟气速度ds表示烟囱直径模型ASC文件GRIDGIS浓度分布图浓度贡献图2 GIS与模型之间的相互作用混合层参数是由大气温度层结a所决定：在中性或不稳定状况下，烟气抬升高度可以近似用以下公式进行计算。 （8）式中F与u与公式（6）的含义相同。由于布里吉斯使用了有关排放源的可获得的数据资料，因此布里吉斯公式相对于其它公式而言更为适用。用FO

9、RTRAN编写的扩散模型程序已达到很高的水平了，程序包括以下方面：SIGY和SIGZ-输入标准值WIND-输入污染物排放时的风向及风速RBD1-输入污染物的特性及排放特征，即：污染源位置、源强、和烟囱直径等。G2-浓度计算程序有关的模型和GIS有很多种，在本文中，我们使用的是应用较多和最简单的一种。在此，模型和GIS系统是二个独立的系统，用计算机文件将它们联系起来。因此，模型从GIS系统的文件中读取输入的数据（如污染源数据），并且也以GIS能够接受的格式作为输出结论。（见图2）。GIS系统有绘图和注释的功能，其绘图的功能是用各种数据库来体现图层。也就是一个整体的图层分为各种形式的专题图层，例如街道、高速公路、绿化、水面等等。每一个图层都由一个或多个包含种各种信息的专题图层所组成：如铁路、道路、河流、湖库、公园、污染源等等。注释数据表是由*.VAT，*.PAT，*.AAT的文件确定的，这些表包括了图层的不同方面（如点、线、面），其