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

高斯扩散模式在瞬间排放空气污染物模拟中的应用 摘要：在文章中提出应用高斯模式模拟和预测在瞬间排放状况下空气污染等级，用FORTRAN语言编写的高斯模式程序还可应用于区域污染影响评价中，模式不仅可以从GIS中输入数据而且还可以应用GIS格式输出结果。 关键词：高斯模式 空气污染 地理信息系统 瞬时污染源 浓度场 瞬间排放是指工业企业或电厂的事故性污染物排放，如贮油罐或输油管道发生事故等。排放的污染物污染了空气、土壤、地面及地下水，影响植被和影响环境。 模拟瞬间空气污染要求得到污染区域面积、污染浓度和等级、污染预测等。 本文提出用高斯模式的解析解来模拟和预测瞬间排放空气污染状况。基于烟羽扩散上的解析公式求解的高斯模式非常广泛的应用于评价区域污染状况。高斯数学模式作为一个污染物扩散的基础模式被国际原子能机构广泛推广。 从瞬间点源污染源排放的污染物，其转换和扩散可以用以下的扩散方程来表示： +div(CV)=?(K?C)+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) 式中：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的夹角为α。这样我们可以将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方向的扩散参数 H表示混合层高度 x0, y0, z0表示坐标系原点 h有效烟囱高度（如，烟气抬升后的高度） 在高斯模式中，σy和σz是重要的参数，可以用帕斯奎尔-吉福德扩散参数曲线来确定，也可将P-G扩散曲线用近似幂函数式表示，这样方便用电子计算机计算大气污染浓度分布。该幂函数的确定与大气稳定度的划分有关。 在高斯模式中，城市和开阔平原的布里吉斯扩散参数以烟云的宽度和高度来表示横风向和垂直向的扩散，σy 和σz的计算方法见表1。 大气稳定度的划分方法见表2。 表1 布里吉斯扩散参数 稳定度 σx, σy (m) σz (m) 开阔平原 A 0.22x(1 + 0.0001x)?1/2 0.2x B 0.16x(1 + 0.0001x)?1/2 0.12x C 0.11x(1 + 0.0001x)?1/2 0.08x(1 + 0.0002x)?1/2 D 0.08x(1 + 0.0001x)?1/