搜索
-1-大气环境容量测算研究（20152207014环境工程李建文）摘要：中国城市的大气污染相当严重，对城市大气污染物实施总量控制将是控制城市大气环境管理污染的主要手段。人们对空气污染问题的关注促进了空气质量数值模式的发展，至今空气质量模式已经发展了三代。目前，国内外常用的模式有ISC3，ADMS，AERMOD，Models-3，CAPPS等。本文主要介绍国内外应用的一些大气环境容量模型，了解各种模型的应用效果，根据实际问题的需求选择合适的空气质量模式，为测算地区环境容量提供依据。关键词：大气，环境容量，模型Abstract:Thepeople'sconcerntotheairpollutionproblempromotedthedevelopmentofairqualitynumericalmodel.Uptothepresent,theairqualitymodelhasalreadydevelopedforthreeeras.Currently,thecommonmodelscontainISC3,ADMS,AERMOD,Models-3,CAPPSetcindomesticandinternational.Therightmodelshouldbeselectedbyactualproblem.Thecalculationmodelsfortheatmosphericcapacityareintroducedinthispaper,andtheapplicabilityofallkindsofmodelsisanalyzed.Accordingtotheneedofactualsituationchooserelevantairqualitymodel,providingrationaleforcalculatingregionalenvironmentalcapacity.Keywords:Atmosphere;Environmentalcapacity;Model.一、前言随着我国工业、交通和建筑业的蓬勃发展，以二氧化碳、氮氧化物和悬浮颗粒物为主的大气污染日趋严重，已经成为我国政府和社会共同面临的严峻问题。大气环境容量是指某一范围内达到或满足大气质量标准所允许排放的最大污染物总量。城市大气环境是一个多变量、多目标、多层次的复杂系统。大气环境规划主要以总量控制为核心，在充分了解当地大气环境背景的基础上，根据大气环境质量管理目标及污染源分布现状，结合城市经济发展和总体规划，确定不同时期大气环境容量开发和利用的最佳方案。只有设计选用较为全面地考虑各种影响因素的一套大气环境总量规划模型，才能较为准确地估计大气环境容量，从而控制大气污染物允许排放总量。目前我国常用的大气污染物总量控制方案的制定方法有A-P-2-值法、多源模拟法等，这些方法各具特色，互有所长，但都有其局限性。二、大气环境容量测算模型1、A-P值法A-P值法为国家标准《制定大气污染物排放标准的技术方法》（GB/T3840-91）提出的总量控制区排放总量限值计算公式；根据计算出的排放量限值及大气环境质量现状本底情况，确定出该区域可容许的排放量。它首先是利用基于箱模型的A值法计算出控制区的大气环境容量，然后利用P值法[1]。在区域内所有污染源排放量之和在不超过上述容量的约束条件下，确定各个点源的允许排放量。在实际过程中，控制区内的高架点源往往对控制区仅产生部分影响(一部分扩散到区外)，以及一些区外高架源对于区域内部也会产生影响。尽管可以通过考虑背景浓度将控制区周边的污染源对区域的影响考虑进来，但这种处理难免显得过于粗糙。在此基础上，王勤耕等[2]于1997年提出了一种改进的A-P值法控制法，他通过引入“影响风向”和“影响份额”的概念，合理地考虑了控制区内外高架点源对控制区的污染物总量的实际贡献。但是他在考虑区域内外高架点源影响时仅仅是以最大落地点的位置来判定污染物的影响作用，未免有点粗糙。如何更加准确地把握高架源对区域的影响是需要充分考虑的。李文慧等[3]以箱模型为基本模型，利用陕西省气象台1961—2008年的气象资料，推导出了大气环境容量宏观总量控制修正A值法，并基于修正A值法估算了西安市各区县大气环境容量，新的国家剩余容量。计算结果表明除临潼区A值大于国标A值外，其余各区县A值小于国标A值。2、ISCLT3（TheIndustrialSourceComplexLong-TermModel）长期平均模式ISCLT3属于高斯型模式，可对大气污染物浓度的月、季、年时空变化规律进行计算和分析，为环境规划、污染控制管理提供科学依据。ISCLT3模式基于高斯模型的湍流平稳、定常的基本假设。在几个月、甚至是一年的长时段内，各种风向均有可能出现，表征短时间的烟流横向扩散能力的σy已不重要，因此需要利用风向频率来计算水平浓度分布（高斯扇面平均公式）。点源排放污染物环境浓度计算：-3-ISCLT3用数值积分的方法来计算面源排放污染物的环境浓度。方法如下：ISC3模式中城市大气扩散参数选用Briggs（城市）扩散参数。烟气抬升采用Briggs计算方法。ISCLT3模式运行需要的参数主要包括：大气污染源排放参数、扩散过程的控制参数、污染气象参数、受体特征参数。根据模型计算结果和本研究提出的颗粒物容量总量控制技术方法的计算结果，我们认为天津市城区控制区的环境容量分别为：SO27.0万吨/年，NO26.2万吨/年，PM104.2万吨/年（煤烟尘），PM1012.7万吨/年（所有颗粒物源类达到二级标准时的允许排放量，相当于煤烟尘的量）。天津市城市控制区的环境容量分别为：SO219.43万吨/年，NO229.63万吨/年，PM1014.18万吨/年[4]。3、箱式模型与高斯复合模型根据污染物质在扩散场中的质量“守恒”原理，大气环境容量计算箱模式一般表述为：VdC=Qdt+VC’dt－VCdt－QDWdt－Q1dt－Q2dt－Q3dt这样，当箱体内污染物浓度达到平衡时，第i种污染物在i级大气质量目标(Cij)下的环境容量为：Qij=V(Cij－C’)+Q1+Q2+Q3。大气质量浓度场模拟，对于点源：SO2采用高斯正态烟云扩散模型；TPS采用高斯扩散—沉积模型。面源采用虚拟点源法。经计算，南通城市一区二镇三个箱体的大气环境容量以单位箱底面积(1km2)计，天生港箱体的单位箱底面积大气环境容量最大，其中：SO2为4580t/a·km2，分别是唐闸箱箱体(365t/a·km2)和城区箱体(162t/a·km2)的12.5倍和25.3倍；TSP为4879t/a·km2，分别是唐闸镇箱体(374t/a·km2)和城区箱体(158t/a·km2)的13.0倍和30.9倍。根据南通城市大气污染物排放预测情况，若一城二镇大气污染源仍按目前的排放方式，则2020年共需削减SO2为59649t/a，削减TSP为7555t/a[5]。复合模型由高斯模型与多维多箱模型组合而成。多维多箱模型建立在质量平衡的基础-4-上，综合考虑污染物的物理、化学去除机理，反映污染物在横向、纵向及高度三维方向上的扩散、污染源分布的不同和气象条件的变化及它们对周围环境的影响。考虑北京市的社会经济情况、城市布局及各种污染源分布密度等情况，在Mapinfo系统中根据北京市空间矢量地图划定多维多箱模型的各个箱体，拟定为38个箱体，依据大气混合层高度，在高度方向上将研究区域分成3个相等部分，构成1个三维多箱箱体。共组成38×3=114个箱体。对每个子箱体写出质量平衡方程，在每个主风向组下，对每个子箱体建1个方程，得到由114个方程构成的方程组。依次对4个主风向组(即N、W、S、E)建立方程组。利用Mapinfo系统的空间矢量地图，以输入污染源数据的北京市空间矢量地图图层为目标图层，用划定的多维多箱空间网格图层对其进行数值处理，在多维多箱空间网格图层属性表上生成污染源数据库，运用Matlab软件，计算出各子箱体污染物平均预测质量浓度。计算点源质量浓度时，依据已划定的各箱体位置，确定高斯模型坐标系，结合北京地区实际输入高斯模型参数，采用高斯模型计算出污染物点源平均质量浓度。研究应用试差法研究确定了北京市2002年不同达标率下的大气PM10环境容量。当达标率为80%时，北京市PM10环境容量为89228·3T/a，削减量为83670·2T/a[6]。3、ISC3模式ISC3(IndustrialSourceComplex3)模式是美国环保局开发的一个为环境管理提供支持的复合工业源空气质量扩散模式，是基于统计理论的正态烟流模式，使用的公式为目前广泛应用的稳态封闭型高斯扩散方程。ISC3模式的模拟范围小于50km，模拟物质为一次污染物，模式采用逐时的气象观测数据，来确定气象条件对烟流抬升、传输和扩散的影响[7]。ISC3扩散模型的核心是高斯烟流模型，它有如下的一些功能：（1）ISC3扩散模型可以用来模拟大气主要污染物和有毒物质及危险废弃污染物质的连续排放。-5-（2）能处理多重来源，包括点、立体、线、面和露天矿等各类源。（3）污染源的源强可按年、季、月、小时等根据需要选取设定。（4）可以计算点源排放时由于附近建筑造成的空气动力学气流下洗的影响。（5）模型包含模拟大微粒（通过干沉积）的沉积和清除的结果以及模拟气体或微粒沉积清除影响的一些算法。（6）ISC3模型使用实时气象数据来计算影响模拟地区的空气污染分布的大气条件。（7）结果输出可以是浓度、总沉积量、干沉积量或湿沉积量。4、AERMOD模式AERMOD模式是在ISC3模式基础上开发得到的，该系统以扩散统计理论为出发点，假设污染物的浓度分布在一定程度上服从高斯分布。模式系统可用于多种排放源(包括点源、面源和体源)的排放，也适用于乡村环境和城市环境、平坦地形和复杂地形、地面源和高架源等多种排放扩散情形的模拟和预测[7]。AERMOD调整模型包含的选项跟ISC3模型基本上是一样的，只有以下一些区别：（1）该模型仅计算浓度值，不能计算干、湿沉降。（2）AERMOD需要两种气象数据：一个文件包含表面标量参数，第二个文件包括垂直横断面。此两个文件由美国EPA-AERMET气象预处理程序提供。使用包括地面海拔的应用程序时，输入接受点海拔的同时用户又要输入山丘高度范围。EPA-AERMET地面预处理程序可以用来产生所有接受点位置的山丘高度和地面海拔。-6-（3）为调试的目的可以要求两种中间结果文件。第一种包括模型结果相关的信息，第二种包括气象数据的网格横断面。（4）AERMOD不支持露天矿各类源。（5）AERMOD包括另外两种输出文件。第一种按排序列出浓度；第二种输出提供弧形最大值结果和与其相关烟流属性的详细文件。研究以广州（梅州）产业转移工业园二期建设为例，通过用AERMOD模型分析工业园区建设带来的大气环境影响，并结合环境质量现状监测数据计算确定园区大气环境容量值，结果得到工业园区SO2、NO2的排放量分别为156.18t/a、145.68t/a，占环境空气容量的46.4％、69.4％。同时，研究结果表明，AERMOD模型结合了预测区域内的环境质量现状、地形地貌、污染环境气象资料等情况，计算出大气环境容量值接近工业园区环境现状情况，与A-P值计算结果对比更适合工业园区大气环境容量计算使用[8]。5、Models-3模式美国环保局研制的第3代空气质量预报和评估系统(Models-3)是由美国环保局野外研究实验室大气模式研制组研制的。该系统可用于多尺度、多污染物的空气质量的预报、评估和决策研究等。Models-3系统由中尺度气象模式、排放模式以及通用多尺度空气质量模式3大模式组成，其核心是通用多尺度空气质量模式(CMAQ)[7]。6、基于CMAQ的大气模型选取空气质量模型CMAQ计算研究区域的大气环境容量。具体方法是通过将不同污染源数据输入模型，选取1、4、7、10四个月份作为研究区域大气环境容量的代表月，根据空气质量模拟结果判断该区域各代表月某种污染物的不同达标率情况下的大气环境容量。在GIS系统的辅助