长沙市空气自动站周边区域大气污染物排放源清单
2018-06-04张青梅罗达通罗华飞李贝睿付广义
张青梅,刘 湛*,罗达通,罗华飞,李贝睿,付广义
(1.湖南省环境保护科学研究院,湖南长沙 410004;2. 苏交科集团股份有限公司,江苏南京 210017)
在支持城市空气质量管理决策的技术体系中,构建准确、完整、更新及时的大气污染物排放源清单是识别污染来源的基础环节,也是制定污染控制策略的根本依据。目前,在城市层面开展的PM2.5来源解析、空气质量预报预警、重污染天气应急方案制定及效果评估、污染物总量减排核查核算、空气质量达标规划等工作,都需要完整的大气污染物排放源清单作为核心基础数据支撑。
大气污染物排放源清单编制在空气污染治理行动中具有重要的意义。美国自20世纪70年代开始清洁空气计划以来,逐步建立了美国国家排放清单(NEI)[1]。欧洲也推出了EMEP/CORINDAIR排放清单编制技术指南规范[2]。我国研究人员在区域大气污染物排放源清单领域也已经开展了大量工作,构建了既符合中国国情又与国际接轨的区域大气污染源排放源清单共性技术体系,将排放系数的本土化率由20%提高到了70%以上,发展了多层嵌套高分辨率区域排放源清单编制技术和方法学。清华大学建立了包括10种污染物、700多种排放源的中国多尺度大气污染物排放源清单(MEIC)[3]。已有研究中以全国尺度或城市区域尺度建立的排放源清单占多数。田贺忠[4]、王丽涛[5]等分别建立了中国NOx排放清单、中国CO人为源排放清单。Xu[6]、Wei[7]等分别建立了全国NH3、VOCs排放源清单。Zhao[8]等建立了8种污染物的全国排放源清单。Fu[9]等建立了长三角区域大气中主要污染物的排放源清单。刘湛[10]、李贝睿[11]等建立了长株潭区域生物质开放燃烧的大气污染物排放源清单和道路机动车污染物排放源清单。而在小区域尺度上,对大气污染源的排放清单的研究较少,谈佳妮[12]建立了上海市宝山区的小尺度精细化大气污染物排放源清单,提高了排放口定位数据的精确度。张骥[13]以天津市津南区空气站周边3 km为研究对象,建立了2016年天津市津南区空气站周边3 km大气污染物排放源清单。实际上较一般区域而言,掌握自动空气站周边区域的污染排放特征更为重要,因为周边区域污染源对空气自动站监测数据的影响较其他区域污染源的影响更大。因此,对空气自动站周边区域开展深入研究对于局地环境空气质量管理、区域大气重污染天气预警等具有更为直接的意义。然而,由于大气污染物排放特征的复杂性,目前已建立的大尺度排放源清单往往不能满足研究的要求。
基于长沙市环境保护管理部门对本地小尺度精细化大气污染物的排放源清单的需求,本研究根据调查搜集到的2015年的各类大气污染物排放源的活动水平信息和环境保护部发布的源排放清单编制技术指南,以长沙市环境空气自动站周边3 km范围为研究区域,建立2015年小尺度精细化大气污染物排放源清单,为各级环境管理部门和科研部门掌握空气站点附近大气污染物排放特征、解析污染成因、制定空气质量持续改善计划提供技术支持。
1 研究方法
1.1 研究区域基本情况
研究区域为长沙市10个国控空气自动站周边的3km内,以各空气站经纬度为中心,利用ArcMAP建立了3km×3km 的网格,结合实地调研和环境统计数据,获取了网格内部的所有点源信息,基本情况如表1所示。
表1 长沙自动空气站基本情况
1.2 活动水平数据获取
工业企业的活动水平信息主要来自环境统计数据,部分信息经长沙市环境保护局调查核实。活动水平数据中煤的灰分、含硫率、废气处理措施通过实地调查获得。民用源燃煤量来自《2016湖南省能源年鉴》[14]和《2016长沙市统计年鉴》[15]。
机动车保有量数据从各市统计年鉴获取,分车型、注册年代、排放控制标准车辆数据从各市环保部门机动车环保检测系统数据库获取。根据机动车燃料、车型及排放控制标准确定道路移动源所属第四级排放源分类。平均行驶里程数据参考《道路机动车大气污染物排放清单编制技术指南(试行)》中各类车型数值。农业机械燃油消耗量活动水平来自《2016湖南农村统计年鉴》[16]。工程机械保有量信息根据《2016中国工程机械年鉴》[17]由全国保有量推算得到。铁路内燃机车燃油消耗量分别根据客货周转量和货运日产量等计算得到,活动水平信息来自《2016中国交通年鉴》[18]。长沙黄花国际机场活动水平数据来自《2016民航机场生产统计公报》。内河船舶燃油消耗量通过客货周转量计算得到,活动水平信息来自《2016长沙市统计年鉴》。
农药使用活动水平信息来自《2016湖南农村统计年鉴》。汽车喷涂、其他涂层活动水平信息来自湖南省环境统计数据。染色、沥青、干洗等活动水平信息根据全国年消费量推算。去污脱脂、生活和商业溶剂使用活动水平信息来自《2016长沙市统计年鉴》中的人口数。建筑涂料采用竣工房屋面积推算,该活动水平信息来自《2016长沙市统计年鉴》。
施肥量来自《2016长沙市统计年鉴》。各区县饲养畜禽数量来自《2016湖南农村统计年鉴》。土壤排放源的活动水平为各区县耕地面积,固氮植物源的活动水平为作物种植面积,秸秆堆肥活动水平利用各区县用于堆肥的八种农作物(水稻、小麦、玉米、粗粮、棉花、豆类、花生和油菜)的作物产量、谷草比、秸秆堆肥比例三者乘积来估算,人体粪便对应的活动水平为没使用卫生厕所的成人数,采用农村常住人口数代替。其他农业源的活动水平数据来自《2016湖南农村统计年鉴》。
土壤扬尘源的活动水平为不同类型土壤的面积,数据来自湖南省国土资源规划院。道路长度通过长株潭区域道路分布的GIS底图结合遥感分类图像解译确定分区域不同等级的道路长度,道路车流量数据由长沙市环保局提供,日降水量信息来自中国气象科学数据共享服务网。施工扬尘源的活动水平数据包括施工工地的面积和施工天数,通过实地调研获得。
生物质燃料的活动水平信息来自《2016湖南农村统计年鉴》。森林火灾、草原火灾活动水平信息来自湖南林业信息网。
储存运输源的活动水平数据通过调研得到。
1.3 大气污染物排放量计算方法
长沙市空气自动站周边区域大气污染物排放量的测算包括“自下而上”的方式和“自上而下”的方式[19,20],基于长沙市大气污染源信息数据库,综合采用物料衡算法、排放系数法、实际测量法、模型估算法等方法测算各类排放源主要污染物的排放强度,建立长沙市大气污染物排放源清单,再在ArcMAP上进行网格分配,得到长沙市空气自动站周边区域内大气污染物排放源清单。排放系数主要根据长沙市的气候、经济发展特点,通过文献调研和实测相结合的方法确定。排放系数选取及排放量计算方法参考的主要技术文件包括:①城市大气污染物排放清单编制技术手册;②大气细颗粒物(PM2.5)源排放清单编制技术指南(试行);③大气可吸入颗粒物一次源排放清单编制技术指南(试行);④大气颗粒物来源解析技术指南(试行);⑤城市扬尘源排放清单编制技术指南(试行);⑥道路机动车排放清单编制技术指南(试行);⑦非道路移动污染源排放清单编制技术指南(试行);⑧生物质燃烧源大气污染物排放清单编制技术指南(试行);⑨大气挥发性有机物源排放清单编制技术指南(试行);⑩氨排放清单编制技术指南(试行)。
2 结果与讨论
2.1 空气自动站周边区域大气污染物排放源清单
本文参考贺克斌主编的《城市大气污染物排放清单编制技术手册》对污染源进行分类,根据所获取的储存运输源、废弃物处理源、工艺过程源、化石燃料固定燃烧源、农业源、生物质燃烧源、扬尘源、移动源8个源类的活动水平数据,按照前面所述的排放源清单编制技术指南进行计算,得到2015年长沙市每个环境空气自动站周边3 km范围内6类污染物NH3、NOx、PM2.5、PM10、SO2、VOCs的排放量,结果见表2。由表2可知,10个站点周边区域内的 NH3、NOx、PM2.5、PM10、SO2、VOCs排 放 量 分 别 53.65t、4 899.35t、1 846.09t、6 257.75t、5 781.67t、4 383.31t。污染物排放量最大的站点为天心区环保局自动空气站周边,NH3、NOx、PM2.5、PM10、SO2、VOCs的排放量之和为9 296.89t,远远高于其他站点周边;其次为经开区环保局自动空气站周边,NH3、NOx、PM2.5、PM10、SO2、VOCs 的排放量之和为3 405.12t;污染物排放量最小的为沙坪自动空气站周边,仅3 447.72t,这与沙坪空气自动站周边点源数量小有关。
图1 不同源类的污染物排放量分担率
2.2 不同源类污染物的分担率
不同源类对10个站点区域内污染物的分担率结果见图1。由图1可知,农业源分担了98.45%的NH3排放量,为NH3第一排放源,畜禽养殖和氮肥施用为主要的NH3排放源农业源[21],这与国内其他研究结果一致[22]。由于空气自动站均位于市区以内,人体排放氨量较小[22],尿液等进入污水处理厂中,因此人体排放活动未进行计算。移动源分担了84.24%的NOx排放量,为NOx第一排放源,李贝睿[11]的研究认为,载货汽车是NOx、PM10、PM2.5的主要贡献源,摩托车是VOCs的主要贡献源。扬尘源对PM2.5、PM10的分担率都最大,分别为60.82%和85.9%,快速的城市化和高强度的城市改造将导致施工工地的扬尘成为城市大气颗粒物的重要来源[23,24]。SO2的最大来源为化石燃料燃烧源,占比为99.54%,由表2可知,天心区环保局站点周边的SO2排放量最大,高达532.49t。这主要是因为天心区燃煤锅炉较其他区域数量要多,燃煤量大。VOCs的最大来源为工艺过程源,占比为49.88%,这是因为北汽福田汽车、三一重工、中联重科、广汽长丰汽车等大型汽车制造企业分布于经开区范围内,而湖南湘江关西涂料有限公司分布于马坡岭网格范围内,汽车制造业与涂料制造业VOCs排放量较大。
表2 各空气自动站周边区域内的大气污染物排放源清单
2.3 不同源类分担率的对比
因为大区域尺度和小区域尺度排放源类有很大差别,因此本研究与已报道的小尺度排放清单研究进行对比,结果见表3。本研究中NH3排放主要来自农业源畜禽养殖和氮肥施用,而宝山地区工业大点源的NH3排放贡献是最高的,该地区电力、热力生产业和供应业和黑色金属冶炼及压延加工业为宝山区的两大主导产业[12]。 NOx、PM2.5、PM10、SO2主要排放源与天津[13]基本保持一致,分别为移动源、扬尘源、扬尘源、化石燃料固定燃烧源。本研究中VOCs排放源主要为工艺过程源,而天津[13]为移动源。对比研究结果表明,小尺度区域内源类不同,排放的污染物不同,小尺度精细化的排放清单对于基层环境管理具有重要意义。各站点区域内NH3、NOx、PM2.5、PM10、SO2、VOCs第一排放源及其占比结果见表4。
2.4 不确定性分析
本研究排放源清单建立过程中选用的排放系数和计算方法均以原环境保护部发布的“清单编制技术指南”为准,并对部分排放系数如移动源的道路机动车污染物排放因子进行了修正,活动水平信息大部分来自政府部门发布的统计年鉴和统计公报,数据来源比较可靠。但是,任何排放源清单的计算结果都具有一定的不确定性,主要原因在于:一方面统计数据本身会存在测量偏差、随机偏差和人为偏差,如统计年鉴数据一般只对规模以上企业进行统计,这些偏差会给大气污染物排放量的计算结果带来一定的不确定性;另一方面由于排放清单的建立涉及的企业众多,各企业的燃料的含硫量、灰分、生产工艺水平、污染控制技术不同,但在排放系数的选定时一般取用参考值或者依据行业平均水平确定,因此也会给排放源清单的计算结果带来不确定性。更精确的不确定性分析需要针对每个源类进行分析,本课题组对生物质开放燃烧源、移动源中机动车尾气、人为氨排放源采用蒙特卡罗法进行了不确定性分析[10,11,19]。
3 结论
(1)2015年长沙市10个自动空气站周边3 km内大气污染源的排放 NH3、NOx、PM2.5、PM10、SO2、VOCs总 量 分 别 为 53.65t、4 899.35t、1 846.09t、6 257.75t、5 781.67t、4 383.31t,其中建筑工地和加油站点源排放的NH3、NOx、PM2.5、PM10、SO2、VOCs总量分别为 6.86t、17.74t、1 240.1t、5 531.61t、96.85t、1 791.97t。
表3 不同源类分担率的对比
表4 不同污染物第一排放源及排放比例
(2)NH3的第一排放源是农业源, NOx的第一排放源是移动源,PM2.5、PM10的第一排放源都是扬尘源,污染比较大的道路是主要干道;SO2的第一排放源是化石燃料燃烧源,VOCs的第一排放源是移动源。
(3)近年来长沙市空气污染多以PM2.5、PM10、O3为首要污染物,要改善空气质量,保障空气站空气环境中PM2.5、PM10、O3监测数据达标,必须要减少一次源排放,根据排放源清单和分担率研究结果可知,改善空气质量应从道路交通优化、燃煤控制、建筑施工场地管控三个方面入手。天心区环保局站点、经开区环保局站点、马坡岭站点还应加强北汽福田汽车、三一重工、中联重科、广汽长丰汽车等工业企业工艺过程污染控制,减少工艺过程源PM2.5、VOCs排放。
[1] JONES D L, ADAMS L, GOODENOW D, et al. Emission Inventory Improvement Program(EIIP) area source committee. Report for November-December [R]. 1994.
[2] European Environment Agency. EMEP/CORINAIR Emission Inventory Guidebook [EB/OL]. (2017-12-06). https://www.eea.europa.eu/publications/EMEPCORINAIR5/.
[3] ZHANG Q. Study on Regional Fine PM Emissions and Modeling in China[D]. Beijing: Tsinghua University, 2005.
[4] 田贺忠, 郝吉明, 卢永琪, 等. 中国氮氧化物排放清单及分布特征[J]. 中国环境科学, 2001, 21(6) 493-497.
[5] 王丽涛, 张强, 郝吉明, 等. 中国大陆CO人为源排放清单[J]. 环境科学学报, 2005, 25(12): 1580-1585.
[6] XU P, KOLOUTSOU V, SOTIRIA R, et al. Projection of NH3emissions from manure generated by livestock production in China to 2030 under six mitigation scenarios[J]. Science of the Total Environment, 2017, 607: 78-86.
[7] WEI W, WANG S X, CHATANI S, et al. Emission and speciation of non-methane volatile organic compounds from anthropogenic sources in China[J].Atmospheric Environment, 2008, 42(20): 4976-4988.
[8] ZHAO B, WANG S X, WANG J D, et al. Impact of national NOxand SO2control policies on particulate matter pollution in China[J].Atmospheric Environment., 2013, 77: 453-463.
[9] FU X, WANG S X, RAN L, et al. Estimating NH3emissions from agricultural fertilizer application in China using the bi-directional CMAQ model coupled to an agro-ecosystem model[J]. Atmospheric Chemistry and Physics, 2015, 15(12): 6637-6649.
[10] 刘湛, 李贝睿, 尤翔宇, 等. 长株潭区域生物质开放燃烧的大气污染物排放清单及不确定性分析[J]. 环境污染与防治, 2016,38(1): 23-29.
[11] 李贝睿, 刘湛, 尤翔宇, 等. 长株潭区域机动车尾气排放清单及特征分析[J]. 环境科学与技术, 2016, 39(11): 167-173.
[12] 谈佳妮, 余琦, 马蔚纯, 等. 小尺度精细化大气污染源排放清单的建立: 以上海宝山区为例[J]. 环境科学学报, 2014, 34(5):1099-1108.
[13] 张骥, 徐媛, 刘茂辉, 等. 天津津南区大气污染物小尺度精细化源清单[J]. 环境科学与技术, 2017, 40(8): 210-215.
[14] 湖南省统计局. 2016湖南省能源年鉴[M]. 北京: 中国统计出版社, 2016.
[15] 长沙市统计局. 2016长沙统计年鉴[M]. 北京: 中国统计出版社, 2016.
[16] 湖南省统计局. 2016湖南农村统计年鉴[M]. 北京: 中国统计出版社, 2016.
[17] 中国机械工业年鉴.2016中国工程机械工业年鉴[M]. 北京: 机械工业出版社, 2016.
[18] 中国交通年鉴.2016中国交通年鉴[M]. 北京:交通出版社, 2016.
[19] 王书肖, 邱雄辉, 张强, 等. 我国人为源大气污染物排放清单编制技术进展及展望[J]. 环境保护. 2017, 21(45):21-26.
[20] 杨静. 珠江三角洲2012年大气排放源清单建立与时空分配改进研究[D]. 广州: 华南理工大学, 2015.
[21] 尤翔宇, 刘湛, 张青梅, 等. 长株潭地区人为源氨排放清单及分布特征[J]. 环境科学, 2016, 37(1): 94-101.
[22] 沈兴玲, 尹沙沙, 郑君瑜, 等. 广东省人为源氨排放清单及减排潜力研究[J].环境科学学报, 2014, 34(1): 43-53.
[23] 邹本东, 徐子优, 华蕾, 等. 因子分析法解析北京市大气颗粒物PM10的来源[J]. 中国环境监测, 2007, 23(2): 79-85.
[24] 徐谦, 李令军, 赵文慧, 等. 北京市建筑施工裸地的空间分布及扬尘效应[J]. 中国环境监测, 2015, 31(5): 78-85.