大气环境容量测算模型简介、方法、内容、过程

1、大气环境容量测算模型简介说明：本部分内容是“重点城市大气环境容量核定工作方案”中提到的各推荐模型的简介，主要目的是为了使各城市了解各模型的功能和基本原理，同时，了解如选用该模型，都需要准备哪些输入数据，以便各城市根据本市的实际情况，提前准备。第一部分 大气扩散烟团轨迹模型1 大气扩散烟团轨迹模型简介该模型由国家环境保护总局环境规划院开发。烟团扩散模型的特点是能够对污染源排放出的“烟团”在随时间、空间变化的非均匀性流场中的运动进行模拟，同时保持了高斯模型结构简单、易于计算的特点，模型包括以下几个主要部分。1.1 三维风场的计算首先利用风场调整模型，得到各预测时刻的风场，由于烟团模型中释放烟团的时

2、间步长比观测间隔要小得多，为了给出每个时间步长的三维风场，我们采用线性插值的方法，利用前后两次的观测风场内插出其间隔时间内各个时间步长上的三维风场，内插公式如下：Vi=V(t1)+V(t2)-V(t1)×n=(t2-t1)Dtin式中： V(t1)、V(t2)分别为第1和第2个观测时刻的风场值；Dt烟团释放时间步长；n为t1、t2间隔内的时间步长数目；Vi表示t1、t2间隔内第i个时间步长上的风场值。1.2 烟团轨迹的计算位于源点的某污染源，在t0时刻释放出第1个烟团，此烟团按t0时刻源点处的风向风速运行，经一个时间步长Dt后在t1时刻到达P11，经过的距离为D11，从t1开始，第一

3、个烟团按P11处t1时刻的风向风速走一个时间步长，在t2时刻到达P12，其间经过距离D12，与此同时，在t1时刻从源点释放出第2个烟团，按源点处t1时刻的风向风速运行，在t2时刻到达P22，其经过的距离为D22，以此类推，从t0时刻经过j个Dt，到tj时刻共释放出了j个烟团，这时，这j个烟团的中心分别位于Pij，i=1，2，j，设源的坐标为（Xs，Ys，Zs(t)），Zs(t)为t时刻烟团的有效抬升高度，Pij的坐标为（Xij，Yij，Zij），u、v分别为风速在X、Y方向的分量，则有如下计算公式：t1时刻：X11=Xs+Ut0,Xs,Ys,Zs(t0)×DtY11=Ys+Vt0,X

4、s,Ys,Zs(t0)×DtZ11=Zs+Wt0,Xs,Ys,Zs(t0)×DtD1=D11=(X11-Xs)2+(Y11-Ys)2t2时刻：X12=X11+Ut1,X11,Y11,Z11×DtY12=Y11+Vt1,X11,Y11,Z11×DtZ12=Z11+Wt1,X11,Y11,Z11×DtD12=D11+D12=D11+(X12-X11)2+(Y12-Y11)2X22=Xs+Ut1,Xs,Ys,Zs(t1)×DtY22=Ys+Vt1,Xs,Ys,Zs(t1)×DtZ22=Zs+Wt1,Xs,Ys,Zs(t1)

5、5;Dt2D2=D22= (X22-Xs)2+(Y22-Ys)2以此类推，到tj时刻，共释放出j个烟团，这些烟团最后的中心位置分别在Pij，Xij，Yij，Zij，i=1，2， j，对于第i个烟团有：Xij=Xi(j-1)+Utj-1,Xi(j-1),Yi(j-1),Zi(j-1)×DtYij=Yi(j-1)+Vtj-1,Xi(j-1),Yi(j-1),Zi(j-1)×DtZij=Zi(j-1)+Wtj-1,Xi(j-1),Yi(j-1),Zi(j-1)×DtDi=j åDk=1jik=Dij-1+(Xij-Xi(j-1)2+(Yij-Yi(j-1)2D

6、ij为i个烟团从源点释放后到tj时刻所经过的距离。1.3 浓度公式由前一个小节的计算，已找到由S点(Xs，Ys)的污染源释放出来的所有烟团在第j个时刻所处的位置，这样S处的污染源在第j个时刻在地面某接受点R(X、Y、0)处造成的浓度就是所有i个烟团的浓度贡献之和。考虑中心位于Pij的烟团对R点的浓度贡献，则有：Ci=CXQs(2p)2sxsyszCX×CY×CZ×Cb×Cdé(X-Xij)2ù=EXPê-ú×22sxêúëûé(Y-Yij)2ùC

7、Y=EXPê-ú22syêúëûCb=EXP(-b×jDt)é(Vd×jDt)2ùCd=EXPê-ú22szëû式中：Qs源强，mg/s；sx、sy、sz：X方向、Y方向、Z方向的大气扩散参数，m；Cx、Cy、Cz：X、Y、Z方向扩散项，Cz在后面给出算式；Cb为污染物转化项，b为转化率，1/s；Cd为污染物沉降项，Vd为沉降速率，m/s。由于考虑到烟团对混合层的穿透作用及混合层对烟团的反射作用，垂直扩散项分以下几种情况讨论：当混合层高为零时（即无混合层时

8、）有：é(Z+Zij)2ùé(Z-Zij)2ùCz=EXPê-ú+EXPê-ú 222sz2szêúêúëûëû计算地面浓度时，Z=0，则有：Cz=EXP(-2Zij2sz2)当混合层高度Zi不为零时，垂直扩散项分以下几种情况计算。设排放源几何高度为hs，混合层高度为Zi，令Zi'=Zi-hs，设烟气抬升高为Dh（烟气抬升高度用“国标HJ/T2.2-93”推荐的模式计算），我们可定义烟气穿透率：p=1.5-Zi'Dh，按

9、不同的P值，分别计算Cz。23Zi'时，认为污染物全在混合层内，按封闭性扩散式计算，当P=0，即Dh£即污染物在混合层与地面间多次反射。Cz=n=-NåNEXP(Zij-2nZi)22sz2式中：N为反射次数，一般取为N=4即可。当P1时，即Dh>2Zi' 时，认为污染物完全穿透混合层，并在混合层以上的稳定层中扩散，由混合层的阻挡而不能到达地面，这时令Cz=0。当0<P<1，即Zi'<Dh<2Zi'时认为是部分穿透情形，这时有部分污染物抬32升到混合层以上，而（1-P）部分被封闭在混合层以内，Cz按下式计算：Cz

10、=Cz1+Cz2许多文献认为穿透到混合层以上的污染物被阻挡后不能向地面扩散，当地区大气层结处于中性偏稳定结构时，混合层对污染物的阻挡作用并不是很强，这时可设计成让这部分烟团在Zij(Zij=hs+Dh)高度上向下扩散，则有：Cz1=P×EXP(-2Zij2sz2)而（1-P）部分的烟团在Zi处按封闭扩散：NCz2=(1-P)n=-NåEXP-(Zij-2nZi)22sz21.4 大气扩散参数1.4.1 有风时扩散参数y、z的确定（0.5h取样时间）(1) 平原地区农村及城市远郊区的扩散参数选取方法如下：A、B、C级稳定度直接由表1.4-1和表1.4-2查算，D、E、F级稳定

11、度则需向不稳定方向提半级后由表1.4-1和表1.4-2查算。(2) 工业区或城区中的点源，其扩散参数选取方法如下：A、B级不提级，C级提到B级，D、E、F级向不稳定方向提一级，再按表1.4-1和表1.4-2查算。(3) 丘陵山区的农村或城市，其扩散参数选取方法同工业区。1.4.2 小风和静风（U101.5m/s) 时，0.5h取样时间的扩散参数按表1.4-3选取表1.4-3 小风和静风扩散参数的系数g01、g02(sx=sy=g01,sz=g02T)1.5 烟气抬升公式1.5.1 有风时，中性和不稳定条件的烟气抬升高度H（m）（1）当烟气热释放率Qh大于或等于是2100KJ/s，且烟气温度与环

12、境温度的差值T大于或等于35K时，H采用下式计算：hDH=noQn1Hn2U-1Qh=0.35PaQuDTTsDT=Ts-Ta式中： no烟气热状况及地表系数，见表1.5-1；n1烟气热释放率指数，见表1.5-1； n2排气筒高度指数，见表1.5-1； Qh烟气热释放率，KJ/s；H排气筒距地面几何高度，m ，超过去240m时，取H=240m； Pa大气压力，hPa； Qv实际排烟率，m3/s；T烟气出口温度与环境温度差，K； Ts烟气出口温度，K； Ta环境大气温度，K；U排气筒出口处平均风速，m/s 。表1.5-1 no、n1、n2的选取（2）当1700 kJ/sQh2100KJ/s时，D

13、H=DH1+(DH2-DH)Qh-1700400DH1=21.5VsD+0.01Qh/U-0.048(Qh-1700)/U()式中： Vs排气筒出口处烟气排出速度，m/s；D排气筒出口直径，m；H2按（1）方法计算，no、n1、n2按表1.5-1中Qh值较小的一类选取；Qh ，U与（1）中的定义相同。DH=2(1.5VsD+0.01Qh)/U （3）当Qh1700kJ/s或者T35K时，1.5.2 有风时,稳定条件按下式计算烟气抬升高度H(m)。DH=Qh1/3ædTaö-1/3ç÷+0.0098 çdZ÷Uèø1

14、/31.5.3 静风和小风时，按下式计算烟气抬升高度H(m).。DH=5.50Qh1/4ædTaö+0.0098÷çèdZø-3/8 但dTadZ取值不宜小于0.01K/m。2 模型运行所需数据数据文件1：共四行：第1行：X方向网格点的数目（MX），Y方向网格点的数目（MY），Z方向风的观测数据层数目（MZ），最大有效烟团数（NT，默认110），污染源数目（MSC），气象观测小时数目（NTimes），稳定度数目（NeleTa，默认24）；第2行：烟团的时间步长（分，默认30.0）第3行：X方向步长（m，默认1000.0），Y方向步长（

15、m，默认1000.0），大气压力(hPa，默认1013.25)，规划区类型（1农村，2城市）；第4行：化学转化率1/s，沉降速率m/s；数据文件2：网格点上的背景浓度值，单位mg/m3；(C1(I,J),I=1,MX),J=1,MY)数据文件3：网格点上的高程（地形值）， 单位m；数据文件4：共六部分。第1部分：污染源数据，共12列，MSC行；第2部分：1500m高程以下各层风的高程，共MZ个；例如：10，50，100，150，200，300，400，500，700，900，1200,1500第3部分：1500m高程以下各温度层的高程，共NEleTa个；例如：0，50，100，150，200，

16、250，300，350，400，450，500，550，600，650，700，750，800，900，1000，1100，1200，1300，1400，1500第4部分：取样时间（024），共Ntime个；例如：6 8 11 14 17 19 21（第1天）30 32 35 38 41 43 45（第2天）54 56 59 62 65 67 69（第3天）78 80 83 86 89 91 93（第4天）102 104 107 110 113 115 117134 137 139 141 （因下雨，有效时次减少） 150 152 155 158 161 163 165174 176 179

17、182 185 187 189198 200 203 206 209 211 213222 224 227 230 233 235 237246 248 251 254 257 259 261270 272 275 278 281 283 285294 296 299 302 305 307 309318 320 323 326 329 331 333342 344 346 348 350 352 354 356 358 360 362 364 366 368 370 372 374 376 378 380 （最后两天时次加密，故有10个数）第5部分：相对于各取样时间的稳定度，共Ntime个；

18、（以数字表示：A=1，B=2，BC=3，C=4，CD=5，D=6，DE=7，E=8，F=9） 例如：8 6 6 4 6 6 69 2 4 6 6 6 69 2 4 6 6 6 66 6 4 6 6 6 66 6 6 4 6 6 66 6 6 66 4 4 2 6 6 88 4 2 4 6 8 68 6 4 6 6 6 88 4 4 6 6 6 88 4 2 6 6 6 68 4 6 6 6 6 66 6 6 6 6 6 66 6 6 6 6 6 98 8 4 4 2 6 6 9 88 8 8 6 6 4 6 6 6第6部分：0-24小时平均逐时混合层高度（m），例如：317.2，285.5，2

19、53.8，222.1，190.4，158.7，127.0，227.0，327.0，455.7，584.3，713.0，733.0，753.0，773.0，695.0，617.0，539.0，507.3，475.6，443.9，412.2，380.5，348.8 数据文件5：1500m高程以下各温度层的温度值（）：共NEleTa列，NTimes行。数据文件6：1500m高程以下各层风的风速，共有NTimes个文件。每个文件格式为：第1行：月， 日， 时；第2行开始N （风的层号）U （风速在x方向上的分量，单位m/s，共MX列）V （风速在y方向上的分量，单位m/s，共MX列）U、 V各MY行N

20、=1，2，MZ第二部分 A-P值法（GB/T 3840-91） 1 A-P值法简介A-P值法为国家标准制定大气污染物排放标准的技术方法（GB/T 3840-91）提出的总量控制区排放总量限值计算公式；根据计算出的排放量限值及大气环境质量现状本底情况，确定出该区域可容许的排放量。1.1 总量控制区内大气污染物排放总量限值的计算方法1.1.1总量控制区污染物排放总量的限值由式（1）计算：Qak=åQi=1naki (1)式中：Qak-总量控制区某种污染物年允许排放总量限值，104t；Qaki-第i功能区某种污染物年允许排放总量限值，104t；n-功能区总数；i-总量控制区内各功能分区的编

21、号；a-总量下标；k-某种污染物下标。1.1.2 各功能区污染物排放总量限值由式（2）计算：Qaki=AkinSiS （2）S=åSi （3）i=1式中： Qaki-第i功能区某种污染物年允许排放总量限值，104t；S-总量控制区总面积，km2；Si-第i功能区面积，km2；Aki-第i功能区某种污染物排放总量控制系数，104t·a-1·km-1，计算方法见1.1.3。1.1.3 各类功能区内某种污染物排放总量控制系数 由式（4）计算：Aki=ACki （4）式中： Aki-第i功能区某种污染物排放总量控制系数，104t·a-1·km-1；Ck

22、i-GB 3095等国家和地方有关大气环境质量标准所规定的与第i功能区类别相应的年日平均浓度限值，mg·m-3；A-地理区域性总量控制系数，104·km2·a-1，可参照表1.1-1所列数据选取。Aki亦可按（GB/T 3840-91）附录A2方法求取。1.1.4 总量控制区内低架源（几何高度低于30m的排气筒排放或无组织排放源）大气污染物年排放总量限值由式（5）计算：Qbk=åQi=1nbki（5）式中： Qbk-总量控制区内某种污染物低架源年允许排放总量限值，104t；Qbki-第i功能区低架源某种污染物年允许排放总量限值，104t，其计算方法见1.

23、1.5；b-低架源排放总量下标。1.1.5 各功能区低架源污染物排放总量限值按式（6）计算。Qbki=aQaki (6)式中： Qbki-第i功能区低架源某种污染物年允许排放总量限值，104t；Qaki-第i功能区某种污染物年允许排放总量限值，104t； a-低架源排放分担率，见表1.1-1。1.1.6 总量控制区内点源（几何高度大于等于30m的排气筒）污染物排放率限值由式（7）计算：Qpki=Pki´He2´10-6 （7）式中： Qpki-第i功能区内某种污染物点源允许排放率限值，t·h-1；Pki-第i功能区内某种污染物点源排放控制系数，t·h-1

24、·m-2 ，计算方法见1.1.7；He-排气筒有效高度，m，计算方法见1.1.11。1.1.7 点源排放控制系数按式（8）计算：Pki=bki´bk´P´Cki (8)式中： Pki -第i功能区内某种污染物点源排放控制系数，t·h-1·m-2； ki -第i功能区某种污染物的点源调整系数，计算方法见1.1.8； k -总量控制区内某种污染物的点源调整系数，计算方法见1.1.9； Cki -见1.1.3定义，但使用日平均浓度限值，mg·m-3; P -地理区域性点源排放控制系数，见表1.1-1。 1.1.8 各功能区点源调整

25、系数按式（9）计算：bki=(Qaki-Qbki)/Qmki (9)式中： ki -见1.1.7定义，若 >1则取 =1；Qaki-见1.1.2定义; Qbki-见1.1.4定义；Qmki-第i功能区内某种污染物所有中架点源（几何高度大于或等于30m、小于100m的排气筒）年允许排放的总量,104t;1.1.9 总量控制区点源调整系数按式（10）计算：bk=(Qak-Qbk)/(Qmk+Qek) (10)式中： k -见1.1.7定义，若 >1则取 =1；Qak -见1.1.1定义；Qbk -见1.1.4定义；Qmk -总量控制区内某种污染物所有中架点源（见1.1.8定义）年允许排

26、放的总量，104t；Qek -总量控制区内某种污染物所有高架点源（几何高度大于或等于100m的排气筒）年允许排放的总量,104t。1.1.10 实际排放总量超出限值后的削减原则是尽量削减低架源总量Qbk及Qbki，使得k和ki接近或等于1，然后再按1.1.7的方法计算点源排放控制系数Pki。1.1.11 排气筒有效高度 按式（11）计算：He=H+DH （11）式中： H-排气筒距地面几何高度，m。超过240m时取H=240m； H -烟气抬升高度，m。计算公式见第一部分1.5烟气抬升公式。1.1.12 点源大气污染物排放浓度（1h平均）限值按式（20）计算：Cpki=2.78´Qp

27、ki´Qv-1´105 (20)式中： Cpki -第i功能区内允许点源烟囱出口处排放的某种大气污染物（1小时平均）浓度限值，mg·m-3；Qpki-见1.1.6定义；Qv-实际排烟率，m3/s。2 所需的输入数据1）总量控制区面积（主要指建成区，不包括大量农田和荒地，它们可作为非控制区，参考控制区执行）；2）总量控制区内的功能分区的面积（若全市空气质量标准皆为二级标准，可按行政区）；3）功能分区的空气质量控制浓度（国家空气质量SO2地面浓度标准年日平均浓度限值及日平均浓度限值）；4）环境平均风速；5）各点源的烟囱几何高度、出口内径、烟气温度、烟气出口速度、源强；

28、6）面源排放面积、平均高度、源强。第三部分 ADMS大气扩散模型软件简介 ADMS大气扩散模型软件是由英国剑桥环境研究公司开发的，分“ADMS-评价”、“ADMS-工业”、“ADMS-城市”等独立系统。其中，“ADMS-城市”版是大气扩散模型系统(ADMS)系列中的最复杂的一个系统。模拟城市区域来自工业，民用和道路交通的污染源产生的污染物在大气中的扩散，ADMS-城市模型用点源，线源，面源，体源和网格源模型来模拟这些污染源。经设计，可以考虑到的扩散问题包括最简单的(例如，一个孤立的点源或单个道路源)到最复杂的城市问题(例如，一个大型城市区域的多个工业污染源，民用和道路交通污染排放)。 它对研究

29、大气质量管理措施特别有用，例如计算先进技术的引进，低排污区对污染状况的影响，燃料的改变，限制车速的设计对空气质量的影响等。ADMS-城市可以作为一个独立的系统使用，也可以与一个地理信息系统联合使用。 ADMS-城市与MapInfo以及ESRI的ArcView可以完全有机的连接。 我们推荐将ADMS-城市与这两个地理信息系统中的任何一种一起使用。 因为这样可以使用数字地图数据，CAD制图和/或航片真实直观地设置您的污染问题。在所使用的不同类型的地图数据上，生成如等值平面图的输出和作报告用的硬拷贝图形等。ADMS-城市与其它用于城市地区的大气扩散模型的一个显著的区别是ADMS-城市 应用了现有的基

30、于MoninObukhov长度和边界层高度描述边界层结构的参数的最新物理知识。其它模型使用Pasquill稳定参数的不精确的边界层特征定义。在这个最新的方法中，边界层结构被可直接测量的物理参数定。这使得随高度的变化而变化的扩散过程可以更真实地表现出，所获取的污染物的浓度的预测结果通常是更精确，更可信。1 适用范围现将ADMS系列模型的应用范围简要描述如下1.1 ADMS - Screen (ADMS-筛选)ADMS-筛选适合用于快速计算来自单个点源的污染物地表浓度并将输出结果自动与中国和其它大气质量标准比较.用户界面简单易懂，只需对污染源和输出结果要求进行最少量的数据输入。ADMS-筛选特别适

31、合于恶劣(最坏)情况下对烟囱源的初步评价,以及对新建工厂的可行性研究进行法律规定的环境影响评价。1.2 ADMS - Industrial (ADMS-工业)ADMS-工业 可计算来自点源,线源，面源和体源的污染浓度.这套系统包括如下特点: 气象预处理模型，干湿沉降，复杂地形的影响，建筑物和海岸线的影响，烟羽可见度，放射性和化学模快; 并可计算短期(秒)内的污染高峰浓度值，如对臭味的预测. 这一系统已于地理信息系统(GIS)连接，易于分析模型结果。ADMS-工业是为计算更详细的一个或多个工业污染源的空气质量影响而 设计的。1.3 ADMS - EIA (ADMS-环评)ADMS-环评可以处理一

32、个区域或城市所有的污染源类型，这包括工业源，道路交通源，面源，体源和网格源等; 适用于区域空气环境评价，并与地理信息系统相连接。1.4 ADMS - Urban (ADMS-城市)ADMS-城市是ADMS模型系列中最复杂的一个。它用于计算来自大区 域和城市的污染浓度或空气质量.此系统可包含所用类型的污染排放源：点源，面源和来自道路的污染.除了具有ADMS-环评的所有特征外，此系统还包一个光化学模型和一个完整连接的地理信息系统(GIS).地理信息系统(GIS)可允许用户在城市地图上显示高分变率的污染浓度图。这可以阐明不同污染源对空气质量的影响，包括道路附近的污染浓度高峰值。ADMS-城市是为详细

33、评价城市区域的大气质量而设计的.这包括空气质量管理战略的发展和城市规化评价。这一套系统同时也用于空气质量预报。2 输入数据需求：源、气象数据、计算点2.1 排污数据污染源类型，污染源位置和污染排放率，流量，烟囱的排放温度，烟囱高度和烟囱直径等。模型带有的排污因子可以使从车流量，平均车速数据来计算道路交通的排污率.污染排放数据可以存储在标准的数据库中或者可有地理信息系统(GIS)系统中读取。2.2 气象数据常见的通用气象数据，如风速，风向，温度，云盖度或者一些演算出来的数值(如莫尼-奥布克夫长度，边界层高度等)。2.3 地形数据除了模拟平坦地形的情形，模拟山地时可以输入合适的地形数据，地表粗糙度

34、数据，街区窄谷的高度，道路宽度和位置。3 输出结果输出结果包括污染物在关心点或关心区的平均浓度，平均时间从10分钟，每小时或到年平均都可。比十分钟平均值还小的短期平均(如秒，分)可由模型的紊动模块处理。输出结果还包括干湿沉降和放射性行为。根据空气质量标准，可以计算从出超标的次数。也可以输出短期紊动，百分位值，滚动平均，概率分布等。气象预处理模块还将产生一个气象输出文件。ADMS生成的输出数字文件为 (.csv) 文件格式。这种文件格式用微软Excel 的电子数据表或在如Windows计事本中可以很容易地打开。ADMS的输出结果可以产生与地理信息系统中的数字地图数据相结合的污染物的等值浓度图。在

35、地表水平的单一点源和烟羽中心线变量的X-Y绘图工具也包括在模型里。 4 模式类型模型是一个三维高斯模型，以高斯分布公式为主计算污染浓度，但在非稳定条件下的垂直扩散使用了倾斜式的高斯模型(Carruthers et al.，1991). 烟羽扩散的计算使用了当地边界层的参数.化学模块中使用了远处传输的轨迹模型和箱式模型。5 污染物种类模型可处理各种基本气态污染物(SO2，NOX，NO2，CO，VOC，苯化物，芳香烃 )，臭氧，可吸入悬浮颗粒物 PM10，PM2.5，总悬浮颗粒物等等。 6 计算源的种类和个数和计算点的个数可详细地模拟3000个网格污染源，1500 个 道 路 污 染 源 和 15

36、00 个 工业污染源 ( 由点，线，面和体污染源) 能够被同时模拟。 由于可将较小的污染源集成为网格污染源,更多非常大数量的污染源在模型实际运行中都可被考虑进去。除了指定的计算点的个数一次可以有50个外，在常规网格输出时，可输出10404 (102X102) 个计算点的值。在智能型网格输出时，除了常规常规网格输出时的计算点外，还可以外加最多5000个计算点。7 烟气过程的处理ADMS使用一个朗格拉日烟羽抬升模块. 烟羽抬升模块预测抬升轨迹和因为热的气态材料的排放对污染浓度的稀释. 其机理是一个顶盖内嵌模型，包括对逆温渗透的处理. 建筑物影响模块基于两条烟羽的交汇，使用平均气流轨迹数据计算烟羽的

37、扩散(Robins et al.，1997). 复杂地形模块是基于线性化的流场模型(Carruthers et al.,1988)。8 水平风的处理对每小时而言，假设风是稳定的。根据不同的边界层参数廓线，垂直方向是不同的。当边界层或复杂地形模块使用时，水平风的描述也相应变化。 9 垂直风速除了在建筑物和复杂地形模块使用了的情况下，垂直风速为零。10 垂直扩散垂直扩散参数从当地的计算出的涡流垂直参数和浮力频率参数演算出来。 11 水平扩散水平扩散从当地的(在烟羽平均高度)计算出的的涡流水平参数演算出来。 12 化学转化化学模块可以计算大气中的一氧化氮(NO)，二氧化氮(NO2)，臭氧(O3)和挥

38、发性有机化合物之间的化学反应. 化学反应使用了GRS (Generic Reaction Set)机理. 包括的化学反应式主要有(1) ROC + hu® RP + ROC (2) RP + NO ® NO2 (3) NO2 + hu ® NO + O3 (4) NO + O3 ® NO2 (5) RP + RP ® RP (6) RP + NO2 ® SGN (7) RP + NO2 ® SNGN其中ROC = Reactive Organic Compounds RP = Radical Pool SGN = Stabl

39、e Gaseous Nitrogen products SNGN = Stable Non-Gaseous Nitrogen products另外，化学反应模块还包含了一个箱式模型。这是GRS功能的延续，它可以用来代表在很大的城区出现在大气边界层内的物理和化学过程。当模拟大型区域时， 地表的臭氧浓度将与背景值不一样。有些时候，一个更具有代表性的并且有空间分布的臭氧场可能要求输入ADMS中的GRS模块。这个臭氧场由箱式模型计算出来， 连接进了ADMS。 它允许ADMS 的输出计算网格将座落在一个更大的空，它的目的是给主要的ADMS化学反应计算提供一个更具有代表性的臭氧场。13 物理清除：重力沉降

40、，干湿沉降对颗粒物的干沉降影响使用了阻力公式，沉降速度考虑到污染物在大气表面层，穿过层流底层到达地面所受污染物阻力的总和，附加重力沉降. 对气态污染物，表面层阻力的计算使用了跟气体性质:活泼的，非活泼的，和惰性的相关的公式. 湿沉降的计算使用定义了的下洗率,和降雨量相关，14 模式比较评估研究资料ADMS模型验证包括:与标准现场，实验室和数字数据系列的结果比较参加欧州共同体的近范围扩散模型的系列技术研讨会与由HMIP(现英国环境署)资助的研究项目中存档的LIDAR数据比较在英国环境与交通部(DETR)所资助的遍布全英国各城市的大气质量评价与回顾的研究项目中，评价扩散模型的表现，并与自动监测数据

41、进行比较。在ADMS-城市中的街道窄谷模型使用了基于丹麦的交通模型(OSPM)。OSPM已被单独验证过。第四部分 ISC-AERMOD模型软件简介 ISC-AERMOD大气扩散模型是由美国Lakes环境公司开发、美国环保署（EPA）推荐的大气扩散模型软件，有界面已经汉化的版本。ISC-AERMOD大气扩散模型软件由ISCST3(Industrial Source Ccomplex Short Term Model)、AERMOD(AMS/EPA Regulatory Model)和ISC-PRIME(Industrial Source Complex - Plume Rise Model En

42、hancement)三大模型组成。1 各模型的功能1.1 ISCST3模型ISCST3扩散模型的核心是高斯烟流模型，它有如下的一些功能：（1）ISCST3扩散模型可以用来模拟大气主要污染物和有毒物质及危险废弃污染物质的连续排放。（2）能处理多重来源，包括点、立体、线、面和露天矿等各类源。（3）污染源的源强可按年、季、月、小时等根据需要选取设定。（4）可以计算点源排放时由于附近建筑造成的空气动力学气流下洗的影响。（5）模型包含模拟大微粒（通过干沉积）的沉积和清除的结果以及模拟气体或微粒沉积清除影响的一些算法。（6）ISCST3模型使用实时气象数据来计算影响模拟地区的空气污染分布的大气条件。（7）

43、结果输出可以是浓度、总沉积量、干沉积量或湿沉积量。1.2 AERMOD模型AERMOD调整模型包含的选项跟ISCST3模型基本上是一样的，只有以下一些区别：（1） 该模型仅计算浓度值，不能计算干、湿沉降。（2） AERMOD需要两种气象数据：一个文件包含表面标量参数，第二个文件包括垂直横断面。此两个文件由美国EPA-AERMET气象预处理程序提供。使用包括地面海拔的应用程序时，输入接受点海拔的同时用户又要输入山丘高度范围。EPA-AERMET地面预处理程序可以用来产生所有接受点位置的山丘高度和地面海拔（3） 为调试的目的可以要求两种中间结果文件。第一种包括模型结果相关的信息，第二种包括气象数据

44、的网格横断面。（4） AERMOD不支持露天矿各类源。（5） AERMOD包括另外两种输出文件。第一种按排序列出浓度；第二种输出提供弧形最大值结果和与其相关烟流属性的详细文件。1.3 ISC-PRIME模型ISC-PRIME模型专门针对当烟流经过两个建筑物时产生气流下洗的基本特点而设计。主要有以下两点：（1） 根据紊流经过时增加的烟流扩散系数。（2） 在建筑物下风向的下降流线和尾流增加的夹杂物的共同作用所造成烟流上升的减少。PRIME算法是与ISCST3模型结合在一起的。因此把它叫做ISC-PRIME模型。ISC-PRIME模型包含了三个新的出入数据文件，用来说明建筑物和烟囱的关系，它们是：（

45、1） 顺风向的建筑物长度。（2） 从烟囱到建筑物迎风面中心的顺风向距离。（3） 从烟囱到建筑物迎风面中心的垂直风向距离。2 所需的输入数据上述三个模块除了已经介召的各自的特点外，它们通用的输入数据主要包含以下内容：2.1 地表气象数据1-3年逐时的风向、风速、对应的风向风速的高度值、地表摩擦速度、对流速度、边界层上空500米层面上的垂直位温梯度、对流产生界面层高度、机械产生界面层高度、Monin-Obukhov长度、表面粗糙长度、温度及取得温度值时的高度。2.2 横截面气象数据横截面气象数据包含了在定点塔上每层的观测结果，或者从NWS数据得来的一层次上的观测结果。该数据文件包含每个小时的一个或

46、多个记录，主要有：年、月、日、小时、测量高度、当前水平上的风向、当前水平上的风速、当前水平上的温度、风向的标准方差、垂直风速的标准方差。2.3 污染源输入数据点源：烟囱底部中心坐标、高度、出口内径、出口烟气温度、烟气出口速度、污染物连续排放量等；面源：坐标、排放高度、排放速度、面源长度、宽度等。其它还可计算体源、露天矿、圆形面、多边形面、火炬、线源等。（略）第五部分EIAA环评助手简介“EIAA大气环评助手”是宁波环科院六五软件工作室开发的软件。是以HJ/T2.2-93 导则-大气环境、JTJ005-96公路建设项目环评规范-大气部分、中国环境影响评价培训教材等文献中推荐的模型和计算方法作为主

47、要框架，内容涵盖了导则中的全部要求，并进行了适当地拓展与加深。EIAA的2.5版采用了面向项目和面向模型两种界面，是两个独立的程序模块，以适应不同的使用习惯和计算要求。在采用“项目预测”模块时，每一个环评项目建立一个独立的“*.Prj”文件，用以保存该项目中用户输入的所有数据、所有方案组合以及所有方案的计算结果。包含了数据预处理、后处理功能，由程序内部选择计算方法。特点是功能强、自动化程度高，但计算过程不透明、操作步骤繁琐。若采用“模型预测”模块，每一个扩散模型都有相对应的程序模块，而且可以同时打开任意多个相同的窗口，便于对比查对。特点是针对性强，使用方便、灵活，便于对扩散过程进行详细研究，但

48、不含数据的预处理和后处理功能，对某些叠加还需手动进行。“模型预测”中的每一个功能模块，均有RTF格式的说明文档与之相对应。这样，对于使用者来说，既可以即时地查看该功能模块的意义、来源，又可以方便地将这些文档直接插入字处理器中。因为这些文档常常包括了大量复杂的公式，这样做就大大提高了写作报告书的速度。此外，EIAA中也提供了电子表格和图形处理功能。能够输出浓度的平面分布图或轴线变化图。1 参数输入与数据预处理可以输入地形数据、现状监测浓度、关心点位置以及当地实测的扩散参数计算系数、风速幂指数、大气温度梯度。这些参数作为本项目的预测方案中的缺省参数。可以从风向、风速和云量的气象观测数据中，按需要统

49、计出多年或某年、某季、某月的风频、各风向风速、稳定度频率、污染系数或者联合频率；也可以从多年风频、风速、稳定度频率数据，按一定的风速段估算出全年联合频率；也可直接输入一个已统计好的联合频率。可以方便绘制各种玫瑰图。可以输入任意个点源、面源、体源、线源和区域面源，同时考虑任意多个污染物。这些输入的参数或经过预处理的气象参数，均可为所有预测方案所共用。 2预测方案、参数选项和预测计算可以用已输入的环境参数、污染气象和污染源，构造各种各样的预测计算方案。程序按输入参数的情况自动选择有风、小风静风、正常排放、非正常排放、点源、面源、体源、区域面源、线源（CALINE4模式）和颗粒物沉降等模式进行计算。

50、可以计算年、季、月、日均浓度。可以同时计算任意多个计算方案，每个方案里可以有任意个污染源。对每一个计算方案，可以控制和调整模式中的每一个参数，在现实可能发生的条件下，允许用户选择计算的模式。例如，可以让一个正常排放的点源在E稳定度下风速为3 m/s时发生薰烟或海岸线薰烟，并输入相关的参数。可选择多个计算方案一起计算，在计算过程中允许调用其它程序。自动保存各方案计算结果。3预测结果分析对于预测计算结果，可以查看各接受点地面高程及其等高线图，各接受点的背景浓度及其分布图，各污染源的浓度和总的浓度及其分布图，各污染源的分担率及其分布图，各污染源或总的浓度的平均评价指数和超标面积。还可以任意改变各污染

51、源的排放率（排放强度）以观察不同排放率下的浓度变化情况。也可查看任意一个横截面或竖截面上的浓度变化图。对所有表格中的数据可以方便地进行各种运算，输出，绘图或打印。用绘图员可以方便地对给定的数据按特定的方式绘制等值线图、X-Y表图、玫瑰图，并且可以编辑修改图形中的每个元素。等值线图也接受评价底图，可将地图（评价底图）输入EIAA，让EIAA将曲线直接绘制在评价底图上。全部图形均以矢量方式打印和保存。对每一个预测计算方案，可以用实测值与预测值进行符合度分析，验证模式的准确性，分析出现偏离的原因，帮助调整有关参数。对每一个预测计算方案，可以按各关心点给定的控制浓度，以总量控制的原则对各污染源进行平权

52、削减分析，给出各污染源的平权污染削减方案。4 其它模块与项目预测不同，模型预测按工业源扩散（代表点源、面源和体源）、公路交通汽车尾气扩散（代表线源）和区域面源扩散（代表城市区域面源）分成三个独立的程序模块，再将薰烟、海岸线薰烟以及其它特殊模式分别孤立开来，由用户按需要自已选择计算模块。此外，“工具”中也提供了一系列具有专一功能的独立模块。第六部分 区域大气污染物总量控制模型软件简介区域大气污染物总量控制模型和软件是国家“九五”重点攻关课题研究成果，也是为了配合全国大气污染物总量控制制度实施而开发的大气污染物总量控制管理软件，本研究课题已通过国家鉴定，并获得国家环保总局科技进步二等奖。 1 基本

53、内容包括基础计算程序组和适用于不同的情况的五组总量模型程序组。（1）基础程序计算组。可以选择各种不同的大气扩散参数，风速扩线指数等和计算参数，确定大气污染物基础排放量和一些基础计算：其中有单源地面浓度计算、考虑混合层的单源地面浓度计算考虑地形影响的单源地面浓度计算和粒物的地面浓度，计算区间和计算参数可以任意选择，并且马上给出地面浓度和画出地面浓度曲线。（2）总量控制计算。根据不同情况采用不同的程序组计算区域大气总量。共有五组总量模型程序组。第一组 采用联合风频的平原地区计算气态污染物总量控制计算；第二组 采用联合风频的平原地区计算颗粒物总量控制计算；第三组 采用逐时气象资料的山区计算气态污染物总量控制计算；第四组 采用逐时气象资料的平原地区计算气态污染物总量控制计算； 第五组 采用逐时气象资料的平原地区计算颗粒物总量控制计算。每组程序可以根据控制区地面浓度，按照软件提示自动或主动修正削减有关大气污染源，直到面浓度达到规定