氨是大气中唯一的碱性气体[1],是参与大气氮循环的重要组成成分. 氨与SO2、 NOx等通过大气化学反应生成的硫酸铵、 硝酸铵等二次颗粒物是大气细颗粒物(PM2.5)的重要组成部分[2]. 此外,氨的排放还会破坏甲烷氧化加剧温室效应[3],NH4+的湿沉降还会导致湖泊富营养化[4]. 因此,随着大气灰霾污染日益严重,国内外学者对氨排放特征的研究也日益深入.
国外在氨排放领域的研究开始较早,对不同排放源和排放特征都有大量深入的研究. 美国环保署在AP-42大气污染物排放因子库的基础上,制定了排放源清单改进计划,编制了包括NH3等多种污染物在内的国家排放源清单,建立了美国国家排放源清单数据库[5]. 欧洲环保署设计了一套针对欧洲国家的排放量计算方法和框架体系,编制了包含NH3等污染物的1980-2005年欧洲国家排放源清单,并实现年度动态更新[6]. Streets等[7]根据美国航空航天局的TRACE-P项目建立了2000年亚洲地区1°×1°的排放源清单,采用排放因子法估算了氨的排放量. 我国在氨排放研究领域起步较晚,已有的研究主要集中在源排放清单及排放特征分析上,王文兴等[1]计算了1991年全国氨的排放情况,研究了氨排放强度地理分布; 孙庆瑞等[8]不仅估算了我国氨排放量,还测定了北京、 湖南等地的氨浓度; Huang等[9]建立了2006年全国氨排放清单,分析了其时空分布规律; 不同学者也分别对长三角[10]、 珠三角[11]、 四川[12]等地区的人为源氨排放特征进行了研究,建立了区域尺度的氨排放清单,部分学者还估算了不同控制措施的减排潜力.
长株潭城市群包括长沙、 株洲、 湘潭三市,是长江中游城市群的重要组成部分,是国家“资源节约型和环境友好型”社会建设综合配套改革试验区,也是我国典型的冶金、 化工产业密集区[13]. 长株潭三市沿湘江谷地呈“品”字型分布,两两相距不足40 km. 近年来随着长株潭地区工业化和城市化的快速推进,能源消费日益增多,大气污染相互叠加,区域性复合型的大气污染特征逐渐显现,以细颗粒物为代表的大气灰霾过程时有发生. 早期针对全国范围的大尺度氨排放清单已经不能满足区域空气质量预报预警和大气雾霾等相关工作的要求. 本研究综合国内外近年来的研究成果,采用“自下而上”的编制方法,通过较为详实的活动水平和本地化的排放因子数据,建立了长株潭地区2013年人为源氨排放清单,并进行了高分辨率的空间分配,以期为区域大气复合污染控制和空气质量改善提供基础数据.
1 材料与方法 1.1 研究区域及研究对象本研究以2013年为基准年,研究区域为长株潭城市群,包括长沙、 株洲、 湘潭这3个地级市所属的23个县区. 估算的人为氨排放源包括农田生态系统、 畜禽养殖业两大主要排放源和生物质燃烧、 人体排泄、 化工生产、 废物处理和移动源等其它排放源.
1.2 活动水平数据本研究农田生态系统和畜禽养殖业排放源活动水平数据来自文献[14],其中,集约化畜禽养殖源活动水平数据来自实际调查. 非农业源活动水平数据来自相应统计年鉴,其中,生物质燃烧排放源活动水平数据来自文献[14, 15],人体排泄等排放源活动水平数据来自文献[16, 17, 18, 19],化工生产源、 烟气脱硝和移动源活动水平数据来自文献[20],污水处理和固废填埋排放源活动水平数据来自文献[21].
1.3 计算方法及排放系数本研究采用排放系数法对长株潭地区人为源氨排放量进行估算,计算公式为:
排放系数选取综合参考国内外研究成果,主要采用文献[22]中推荐的排放系数本地化计算方法,根据长株潭地区的具体情况计算得到本地化的排放系数. 在具体估算过程中,不同排放源的估算公式、 排放系数确定方法、 活动水平数据处理方法等有所不同,下面分别介绍.
1.3.1 农田生态系统农田生态系统的氨排放源包括氮肥施用、 土壤本底、 固氮植物和秸秆堆肥.
氮肥施用后,氮通过物理化学反应以氨的形式挥发到大气中,是大气中氨的一个重要来源[23]. 氮肥施用后氨的挥发与氮肥的施用种类、 施用方式、 施用量、 土壤的温度和酸碱度等都有关系[1]. 氮肥施用中氨的排放系数计算公式为:
长株潭地区不同氮肥施用份额及基准排放系数见表 1所示,其中以全国不同类型氮肥的年生产量作为长株潭地区不同氮肥的施用比例[10],氮肥施用的平均基准排放系数是根据不同氮肥的施用比例和基准排放系数确定.
对于农田生态系统中的其他排放源,其排放系数如表 2所示.
秸秆堆肥活动水平为可用于堆肥的8种农作物,包括水稻、 小麦、 玉米、 粗粮、 棉花、 豆类、 花生和油菜的田间堆肥的秸秆量. 秸秆量按下式计算:
畜禽养殖业中的氨排放主要由动物排泄物释放,包括圈舍饲养、 户外放牧、 粪便存储处理和后续施肥这4个阶段[26]. 根据湖南省农业统计资料[14],本研究将畜禽养殖源细分为牛、 羊、 肉猪、 母猪、 蛋鸡、 肉鸡、 鸭、 鹅和马驴骡9类,其中牛、 羊、 母猪、 蛋鸡和马、 驴、 骡等饲养周期大于1年的,养殖数量按统计资料中的年底存栏数计算,对于肉猪、 肉鸡、 肉鸭、 肉鹅等肉用家禽,养殖数量按统计资料中的年内出栏数计算. 其中畜禽养殖分为集约化畜禽养殖和散养,养殖阶段分为户外、 圈舍、 存储和施肥这4个阶段,排泄物种类分为液态和固态分别计算. 不同养殖阶段和不同排泄物种类铵态氮(TAN)的计算公式为:
A户外=TAN户外
A圈舍=TAN室内×X
A存储=A圈舍-A圈舍×EF圈舍
A施肥= [A存储-A存储×EF存储-ENN损失]×(1-R饲料)
其它行业主要包括生物质燃烧、 化工生产、 人体排放、 废物处理和道路移动源. 其排放系数见表 4所示.
(1)生物质燃烧源
主要包括森林大火、 草原大火、 秸秆燃烧、 薪柴燃烧.
森林火灾消耗生物量计算公式为:
长株潭地区草地面积较小,且年内没有草原火灾发生,因此不进行计算.
秸秆露天焚烧量根据秸秆产量的比例计算,秸秆露天焚烧比例一般取为20%,燃烧率一般取为90%.
(2)化工生产
主要包括合成氨和氮肥生产.
(3)人体排放
人体排放氨主要通过人体呼吸、 汗液和粪尿排出[11],由于人体呼吸和汗液排放量较小,因此本研究中仅计算人体粪尿排泄物的氨排放,城镇家庭大多有卫生设施,相应的氨排放计算在污水处理中,人体排放活动水平仅计算卫生设施不完善的农村人口.
(4)废物处理
包括污水处理、 固废填埋和烟气脱硝过程. 烟气脱硝主要包括发电供热企业和水泥生产企业等.
(5)道路移动源
随着城市机动车保有量的持续增长,尾气排放的氨量逐渐增加. 本研究根据车辆类型分为轻型汽油、 轻型柴油、 重型汽油、 重型柴油和摩托车这5类.
长株潭地区是长江中游城市群的重要组成部分,是湖南省经济增长的核心引擎. 2013年该地区常住人口为1.40×107人,其中农村常住人口4.95×106人,耕地总面积为6.34×103 km2,全年主要农作物总产量为5.98×106 t,其中水稻的产量为5.56×106 t,按式(3)计算得到的长株潭地区全年各类农作物秸秆产生量为7.97×106 t. 2013年长株潭地区汽车保有量为2.27×106辆,其中摩托车1.14×106辆.
2013年长株潭地区全年氮肥施用量(实物量)6.57×105 t. 氮肥施用量较大的县区是湘潭县、 宁乡县、 湘乡市和攸县,分别占长株潭地区氮肥施用总量的23.17%、 13.29%、 11.87%和10.77%.
2013年长株潭地区各类畜禽集约化养殖比例约为43%. 牛、 羊、 母猪、 蛋鸡的年末存栏数分别为3.89×105头、 1.19×106头、 1.09×106头、 2.27×107头,生猪、 肉鸡、 肉鸭(鹅)的年内出栏数分别为1.81×107头、 7.44×107头、 1.53×107头. 宁乡县、 浏阳市和醴陵市是畜禽养殖量最大的县区,其各类畜禽养殖总量占长株潭地区总量的31.70%、 13.93%和10.81%.
2.2 长株潭地区人为氨源排放清单长株潭地区大气氨源排放清单如表 5所示,从中可以看出,长株潭地区2013年大气氨源NH3排放总量为72.62 kt,畜禽养殖业是长株潭地区最大的氨排放源,排放量为42.60 kt,占氨排放总量的58.66%. 农田生态系统是第二大排放源,排放量为21.54 kt,占氨排放总量的29.66%,两者排放量之和占氨排放总量的88.32%,说明农牧业是长株潭地区大气氨排放的主要贡献源,主要原因在于长株潭地处亚热带,是我国重要的畜禽养殖和农作物种植基地; 另外碳酸氢铵在我国氮肥施用中占比较大,而其又较容易以氨的形式挥发,也造成了氨排放量较大. 另一方面,非农业源氨的排放也不可忽视,人体排放、 废物处理和化工生产的氨排放量分别占总量的5.36%、 2.13%和1.90%.
畜禽养殖是长株潭地区最主要的氨排放源,各类牲畜氨排放量占畜禽养殖业氨排放总量的百分比如图 1所示. 从中可知,肉牛、 蛋鸡和肉猪是排放量最大的3种畜禽,氨排放量分别为11.18、 9.12和7.85 kt,分别占畜禽养殖业氨排放总量的26.26%、 21.40%和18.43%,而肉鸡、 肉鸭、 肉鹅和马、 驴、 骡等的贡献率相对较低. 从养殖数量上来看,鸡、 鸭、 鹅和猪是长株潭地区2013年养殖量最大的牲畜种类,而肉鸡、 肉鸭、 肉鹅的分担率较低是因为其年均氮排泄量远低于生猪、 牛等牲畜.
农田生态系统是长株潭地区第二大氨排放源,其中氮肥施用的贡献率最大,占农田生态系统氨排放量的92.44%. 主要是由于我国氮肥施用以碳铵和尿素为主,两者使用比例约90%[27],而这两种氮肥的氨挥发率较高,欧美等发达国家多施用氨挥发率较低的硝铵和硫铵,加之我国大多采用表面撒施为主的施肥方式,更加加重了氮肥施用过程中氨的挥发.
2.3 人为氨源排放的地域分布根据长株潭地区各区县不同排放源的活动水平和NH3的排放因子,计算得到各区县人为源氨排放清单,如图 2所示. 由于人类活动水平的差异,各区县大气氨的排放呈现出一定的地域分布特征. 从图 2中可以看出,宁乡县、 湘潭县和浏阳市是长株潭城市群氨排放量最高的城市,排放量分别为1.27×104、 9.32×103和8.74×103 t,分别占长株潭地区氨排放总量的17.49%、 12.82%和12.02%. 宁乡县和浏阳市主要是由于其为长株潭地区重要的畜禽养殖基地,畜禽养殖量大造成的,其畜禽养殖业的氨排放量分别占各自氨排放总量的70.09%和76.97%,湘潭县是由于其化肥施用量大,每亩耕地施肥以氨折纯量计高达38 kg,导致农田生态系统排放的氨总量偏大.
从排放强度来看,石峰区和岳塘区的排放强度最高,分别达到9.14 t ·km-2和5.01 t ·km-2,远高于其他区县,主要是由于其为长株潭城市群重要的工业生产基地,辖区内集中分布了多家大型化工、 钢铁、 有色和火电企业,工业生产和烟气脱硝过程中排放的氨量较大.
2.4 大气氨源排放清单比较通过文献收集,总结不同学者对国内不同地区大气氨源排放的研究,与本研究结果进行对比,如表 6所示. 表明畜禽养殖和氮肥施用是最主要的人为氨排放源,其次是人体排放.
从排放强度上来看,关于湖南省的氨排放强度,本研究得到的长株潭地区平均氨排放强度为2.59 t ·km-2,高于文献[1, 8]的结果,低于董文煊等[28]的研究结果. 与国内外其他研究结果对比来看,本研究结果略高于日本[29](1.10 t ·km-2)和欧洲[30](1.27 t ·km-2),接近四川省[12]排放强度(2.12 t ·km-2),低于长三角[10](4.20 t ·km-2)和广东省[26](3.69 t ·km-2).
从不同源强的贡献率上来看,本研究中长株潭地区的畜禽养殖源氨排放贡献率达到了58.60%,与珠三角[11](62.1%)、 四川省[12](62.31%)和湖南省[1](64%)的研究结果相近,以氮肥施用为主的农田生态系统氨排放贡献率为29.73%,也与珠三角[11](21.7%)、 四川省[12](23.14%)和湖南省[1](18%)的研究结果相近. 第三大排放源人体粪便氨的排放贡献率为5.36%,与广东省[26](6.0%)和四川省[12](6.39%)的结果接近.
不同研究得到的氨排放量不同,一方面是因为在不同年份和地区,社会经济发展、 人口数量、 人们物质需求和农牧业结构等客观条件不同,导致活动数据有所差异. 另一方面,不同学者估算的排放源类别和采用的排放因子不同,也造成了排放清单结果的较大差异. 早期的排放清单选取的排放系数多引用国外数据,然而,随着我国学者对氨排放研究的不断深入,逐渐采用实测建立本地化排放系数,以减小氨排放清单不确定性. 本研究中以氮肥施用为主的农田生态系统的贡献率偏低,就是因为长株潭地区的土壤以红壤为主,pH值一般在4.5-6.0之间,年平均地面温度为18.9℃,根据文献[22]计算得到的酸性土壤中氮肥的基准排放系数与碱性土壤相比,相差5.84倍,导致氮肥施用的贡献率下降.
本清单在前人研究的基础上,尽可能选择本地化排放因子,扩大排放源覆盖面,与不同学者的研究结果对比,本研究估算的排放量较为合理,排放强度与其他地区相比基本一致. 因此,研究结果能较好地反映长株潭城市群人为源氨的实际排放情况.
2.5 氨排放清单的不确定性分析排放源清单估算分析过程中不确定性来源[33]主要包括活动水平数据的缺失和排放系数的本地化程度不足等. 本研究中,畜禽养殖业活动水平数据主要来自相关统计年鉴,集约化畜禽养殖场信息来自现状调研,排放因子按照不同粪便管理阶段分别计算,其不确定性相对较小. 排放系数选取文献[22]中推荐的计算方法依据本地化的参数计算得到,不确定性较小. 氮肥施用是本研究中主要的不确定来源,氮肥施用量采用自上而下进行分配,没有根据不同氮肥划分,排放因子是折算的平均值,故其不确定性较大. 人体排放、 污水处理和固废填埋、 道路移动源均采用自下而上的方法进行估算,活动水平数据相对较为可靠,不确定性主要来自于排放系数. 生物质燃烧中家用薪材和秸秆的燃烧量是通过估算得出,然后采用自上而下的方法进行分配,存在着一定的误差. 此外,秸秆露天燃烧、 森林火灾和草原火灾的燃烧量计算中的相关参数和排放系数的选取,由于缺少本地实测数据,只能参考国外或其他地区研究结果,森林火灾活动水平是按照各县市区的林地面积对湖南省的火灾森林面积分配得到的,这也增大了估算结果的不确定性.
本研究采用Monte Carlo方法对不同人为源及总排放清单的不确定性进行分析,结果如表 7所示.
总体而言,估算过程中不确定性的主要来源为活动水平数据的选取和排放因子的选取. 从Monte Carlo模拟结果来看,农田生态系统、 人体排放、 废物处理和生物质燃烧不确定度较高,畜禽养殖和化工生产的不确定度较低. 农田生态系统主要是由于基础统计数据中不同氮肥施用比例、 含氮百分比和秸秆堆肥比例过于笼统,不够明确造成的,人体排放是由于农村人口未使用厕所的人数不明确造成的,废物是由于排放因子过于老化和单一造成,生物质燃烧是由于湖南地区秸秆露天焚烧和家庭燃烧基础数据不够全面及部分排放因子单一造成的.
2.6 氨排放清单的空间分配 2.6.1 空间分配方法空间分配是在空间特征表征数据的基础上采用基于GIS的计算方法[34]将一类污染源的排放量分配到对应的网格中. 本研究基于长株潭地区污染源空间分布规律和空间特征表征数据,建立了3 km×3 km的长株潭地区网格化人为源氨排放清单. 点源的空间分配直接根据排放源的经纬度将排放量分配到点源所在的网格中,面源和线源根据空间表征数据建立与模型区域相一致的空间分配因子进行分配. 常见的空间特征表征数据包括人口密度分布数据、 土地利用分类数据、 道路网数据等.
2.6.2 空间分布特征基于上述空间分配方法,得到了长株潭地区2013年人为源氨排放空间分布,如图 3所示. 从中可以看出,长株潭地区的氨排放主要分布在西部,主要原因是西部地区人口密度大,距离中心城区近,畜禽饲养量较大. 而东北部和东南部地区由于相对人口密度低,氨排放量也较小. 排放量较大的点源主要是一些大型的电厂、 化肥厂等工业点源和大型的畜禽养殖场. 由于农田生态系统和畜禽养殖业是大气氨排放的主要来源,而城区周边集中了大量的集约化畜禽养殖场,因此城区周边区域的氨排放量相对较高.
(1)本研究通过收集汇总长株潭城市群各区县的农业、 工业生产基础信息以及生活能源消耗数据,建立了2013年长株潭城市群人为源氨排放清单. 长株潭地区2013年人为源氨排放总量为72.62 kt,排放强度为2.59t ·km-2.
(2)畜禽养殖和农田生态系统是长株潭地区氨排放的主要来源,排放量分别为42.60 kt和21.54 kt,分别占氨排放总量的58.66%和29.66%. 畜禽源中,猪和牛是最大的排放源,分别占畜禽养殖氨排放总量的29.52%和26.70%.
(3)宁乡县、 湘潭县和浏阳市为排放量最大的3个区县,石峰区和岳塘区为排放强度最大的两个县区. 无论从排放量和排放强度上看,宁乡县和湘潭县都居前列.
(4)从空间分布特征来看,主城区氨排放量较小,城区周边排放量较大,排放量较大的网格主要是大型点源.
[1] | 王文兴, 卢筱凤, 庞燕波, 等. 中国氨的排放强度地理分布[J]. 环境科学学报, 1997, 17 (1): 2-7. |
[2] | Goebes M D, Strader R, Davidson C. An ammonia emission inventory for fertilizer application in the United States[J]. Atmospheric Environment, 2003, 37 (18): 2539-2550. |
[3] | 杨志鹏. 基于物质流方法的中国畜牧业氨排放估算及区域比较研究[D]. 北京: 北京大学, 2008. |
[4] | 叶雪梅, 郝吉明, 段雷, 等. 中国主要湖泊营养氮沉降临界负荷的研究[J]. 环境污染与防治, 2002, 24 (1): 54-58. |
[5] | 郑君瑜, 王水胜, 黄志炯, 等. 区域高分辨率大气排放源清单建立的技术方法与应用[M]. 北京: 科学出版社, 2014. |
[6] | European Environment Agency. EMEP/EEA air pollutant emission inventory guidebook-2009[EB/OL]. http://www.eea.europa.eu/publications/emep-eea-emission-inventory-guidebook-2009, 2009-06-19. |
[7] | Streets D G, Bond T C, Carmichael G R, et al. An inventory of gaseous and primary aerosol emissions in Asia in the year 2000[J]. Journal of Geophysical Research, 2003, 108 (D21): 8809. |
[8] | 孙庆瑞, 王美蓉. 我国氨的排放量和时空分布[J]. 大气科学, 1997, 21 (5): 590-598. |
[9] | Huang X, Song Y, Li M M, et al. A high-resolution ammonia emission inventory in China[J]. Global Biogeochemical Cycles, 2012, 26 (1): GB1030. |
[10] | 董艳强, 陈长虹, 黄成, 等. 长江三角洲地区人为源氨排放清单及分布特征[J]. 环境科学学报, 2009, 29 (8): 1611-1617. |
[11] | 尹沙沙, 郑君瑜, 张礼俊, 等. 珠江三角洲人为氨源排放清单及特征[J]. 环境科学, 2010, 31 (5): 1146-1151. |
[12] | 冯小琼, 王幸锐, 何敏, 等. 四川省2012年人为源氨排放清单及分布特征[J]. 环境科学学报, 2015, 35 (2): 394-401. |
[13] | 杨晴, 张凯, 柴发合, 等. 长株潭城市群秋季大气颗粒物及其重金属元素污染特征[J]. 环境科学研究, 2013, 26 (6): 590-597. |
[14] | 湖南农村经济调查队. 2014湖南农村统计年鉴[M]. 北京: 中国统计出版社, 2014. 91-141. |
[15] | 湖南省林业厅. 2014年全省林业灾害控制情况[EB/OL]. http://www.hnforestry.gov.cn/MoreList_Infoshow.aspx?ID=425580, 2015-04-13. |
[16] | 湖南省统计局. 2014湖南统计年鉴[M]. 北京: 中国统计出版社, 2014. 625-676. |
[17] | 长沙市统计局. 2014长沙统计年鉴[M]. 北京: 中国统计出版社, 2014. 187-398. |
[18] | 株洲市统计局. 2014株洲统计年鉴[M]. 北京: 中国统计出版社, 2014. 403-513. |
[19] | 湘潭市统计局. 2014湘潭统计年鉴[M]. 北京: 中国统计出版社, 2014. 199-319. |
[20] | 湖南省环境保护厅. 2014湖南环境统计年鉴[M]. 北京: 中国统计出版社, 2014. |
[21] | 中华人民共和国住房和城乡建设部. 中国城市建设统计年鉴2013[M]. 北京:中国统计出版社, 2014. 494-545. |
[22] | 环境保护部. 大气氨源排放清单编制技术指南[EB/OL]. http://www.zhb.gov.cn/gkml/hbb/bgg/201408/W020140828351293771578.pdf, 2014-08-19. |
[23] | Zhang Y, Dore A J, Ma L, et al. Agricultural ammonia emissions inventory and spatial distribution in the North China Plain[J]. Environmental Pollution, 2010, 158 (2): 490-501. |
[24] | 湖南省农业厅. 湖南土壤[M]. 北京: 中国农业出版社, 1989. 247-266. |
[25] | 国家发展改革委, 农业部, 财政部. 关于印发"十二五"农作物秸秆综合利用实施方案的通知[EB/OL]. http://www.mof.gov.cn/zhengwuxinxi/zhengcefabu/201112/t20111221_617842.htm, 2011-11-29. |
[26] | 沈兴玲, 尹沙沙, 郑君瑜, 等. 广东省人为源氨排放清单及减排潜力研究[J]. 环境科学学报, 2014, 34 (1): 43-53. |
[27] | 苏芳, 丁新泉, 高志岭, 等. 华北平原冬小麦-夏玉米轮作体系氮肥的氨挥发[J]. 中国环境科学, 2007, 27 (3): 409-413. |
[28] | 董文煊, 邢佳, 王书肖. 1994-2006年中国人为源大气氨排放时空分布[J]. 环境科学, 2010, 31 (7): 1457-1463. |
[29] | Kannari A, Tonooka Y, Baba T, et al. Development of multiple-species 1km×1km resolution hourly basis emissions inventory for Japan[J]. Atmospheric Environment, 2007, 41 (16): 3428-3439. |
[30] | Asman W A H. Ammonia emission in Europe: updated emission and emission variations[R]. The Netherlands: National Institute of Public Health and Environmental Protection, 1992. |
[31] | 广东省环境监测中心站, 香港特别行政区政府环境保护署. 珠江三角洲地区空气污染物排放清单编制手册[R]. 广州: 广东省环境保护局, 2005. |
[32] | 陈晓秋, 陈进生, 吴水平, 等. 海峡西岸城市群大气复合污染成因与调控机制[M]. 北京: 科学出版社, 2014. 246-253. |
[33] | 钟流举, 郑君瑜, 雷国强, 等. 大气污染物排放源清单不确定性定量分析方法及案例研究[J]. 环境科学研究, 2007, 20 (4): 15-20. |
[34] | Sutton M A, Place C J, Eager M, et al. Assessment of the magnitude of ammonia emission in the United Kingdom[J]. Atmospheric Environment, 1995, 29 (12): 1393-1411. |