环境科学  2022, Vol. 43 Issue (10): 4380-4391   PDF    
浙江省2008~2018年人为源氨排放清单及分布特征
方利江1, 杨一群2, 叶观琼2     
1. 舟山市生态环境保护技术中心, 舟山 316000;
2. 浙江大学海洋学院, 舟山 316000
摘要: 以浙江省为研究区域,通过收集11个地级市各类氨排放的活动水平数据,采用排放因子法建立了2008~2018年浙江省人为源氨排放清单,并利用ArcGIS进行1 km×1 km空间网格分配.结果表明,2008~2018年浙江省人为源氨排放量总体呈现下降趋势,年均下降率约3.97%;2018年浙江省氨排放量为108.52 kt,其中农业源为90.02 kt,非农业源为18.50 kt,排放强度为1.03 t ·km-2;杭州市、嘉兴市和温州市的氨排放量高于其他城市,分别占全省总量的14.72%、11.86%和11.80%;空间分布特征显示,氨排放主要分布在浙江省北部区域,呈现出"北高南低"的排放趋势;不确定性分析表明,氨排放量模拟平均值为108.37 kt,在95%置信区间的不确定度范围为-5.40%~5.60%.
关键词: 浙江省      人为源      氨排放清单      年际变化      空间分布特征      不确定性分析     
Inventory and Distribution Characteristics of Anthropogenic Ammonia Emissions in Zhejiang Province from 2008 to 2018
FANG Li-jiang1 , YANG Yi-qun2 , YE Guan-qiong2     
1. Zhoushan Ecology and Environment Protection Technical Center, Zhoushan 316000, China;
2. Ocean College, Zhejiang University, Zhoushan 316000, China
Abstract: As ammonia is the main component of PM2.5, long time series of ammonia emission characteristics are an important basis for studying the historical causes of PM2.5 pollution. In this study, the activity data of various anthropogenic ammonia emissions from 11 cities were collected in Zhejiang. The anthropogenic ammonia emissions inventory in Zhejiang was established using emission factors, and then a 1 km×1 km spatial grid distribution was carried out using ArcGIS software. The results showed that from 2008 to 2018, the ammonia emissions from anthropogenic sources in Zhejiang exhibited a downward trend, with an average annual decline rate of approximately 3.97%. The ammonia emissions were 108.52 kt in 2018, and the emission intensity was 1.03 t·km-2, in which there was 90.02 kt from agricultural sources and 18.50 kt from non-agricultural sources. The ammonia emissions of Hangzhou, Jiaxing, and Wenzhou were higher than those of the other cities, accounting for 14.72%, 11.86%, and 11.80% of the total ammonia emissions, respectively. The spatial distribution characteristics showed that ammonia emissions were mainly distributed in the northern part of Zhejiang, showing an emission trend of "high in the north and low in the south." Uncertainty analysis showed that the simulated average value of ammonia emissions was 108.37 kt, and the uncertainty range in the 95% confidence interval was -5.40%-5.60%.
Key words: Zhejiang Province      anthropogenic sources      ammonia emissions inventory      interannual variation      spatial distribution characteristics      uncertainty analysis     

氨(NH3)是大气中唯一的天然碱性气体和第三高含量的含氮化合物[1], 气态氨浓度的增加会严重影响区域空气质量、人类健康和生态系统的氮循环[2].排放到大气中的NH3可以与二氧化硫(SO2)和氮氧化物(NOx)等酸性物质经化学反应转化为硝酸铵(NH4NO3)、硫酸铵[(NH4)2SO4]、硫酸氢铵(NH4HSO4)和氯化铵(NH4Cl)等铵盐, 成为PM2.5的重要组成部分, 是造成大气PM2.5污染的重要排放来源[3].形成的铵盐可以远距离传送到较远地方, 通过干湿沉降返回土壤或水体中, 导致土壤酸化和水体富营养化甚至生物多样性下降, 从而引起区域性大气环境污染问题[4, 5].有研究表明, 1990~2013年间, 全球氨排放对PM2.5贡献率从25%增加到32%[6].在美国约有47%的PM2.5由铵盐贡献[7], 在中国约占29.8%[8], 在欧洲为30%[9].如果NH3排放减少50%, 将导致铵盐减少22%, 硝酸盐减少39%, PM2.5总量则减少15%[10].

近年来, 随着对气溶胶污染形成和作用机制研究的不断深入, 国内学者对氨排放越发关注, 建立了国家级、区域级和省市级等不同尺度的氨排放清单.Fu等[11]估算了1980~2016年中国的氨排放量, 从4.7 Tg·a-1增长到11 Tg·a-1. Yu[12]建立了长三角地区2006~2014年的氨排放清单, 排放量范围在981.65~1 014.30 kt.王琛等[13]认为河南省2013年的氨排放量为1 035.3 kt, 排放强度为6.4 t·km-2.全面掌握氨排放清单和空间分布规律, 有助于了解氨排放对当地气溶胶污染的影响, 是阐明区域性复合型空气污染问题的关键.

浙江省位于中国东部沿海, 是国家确定的蓝天保卫战重点区域之一, 尽管在“十三五”期间, 通过实施“打赢蓝天保卫战三年行动计划”等系列措施, 大气污染防治攻坚战取得了阶段性胜利, 2020年设区市ρ(PM2.5)平均值为25 μg·m-3, 比2015年降低43.2%[14], 但以PM2.5为代表的大气复合型污染天气仍时有发生, “十四五”时期进一步改善环境空气质量, 实现全省“全面消除重污染天气, 基本消除中度污染天气”的目标依旧面临较大压力.目前, 浙江省的氨排放清单已有相关研究成果[15, 16], 但缺乏长时间序列和较全面排放源的评估计算.本文采用“自下而上”的评估方法, 建立了2008~2018年浙江省人为源氨排放清单, 并首次计算了以港口船舶为典型代表的非道路移动源氨排放量, 研究结果对解析浙江省大气污染的历史成因、溯源和作用机制等具有积极意义, 也可为浙江省大气污染治理和空气质量管理决策提供有力的数据支持和科学参考.

1 材料与方法 1.1 研究区域与对象

研究区域为浙江省, 包括杭州市、宁波市、温州市、嘉兴市、湖州市、绍兴市、金华市、衢州市、舟山市、台州市和丽水市等11个地级市.估算的人为源包括农业源和非农业源共计11种排放源.其中, 农业源包括畜禽养殖、氮肥施用、土壤本底、固氮植物和生物质燃烧; 非农业源包括人体排放、工业生产、废物处理、燃料燃烧、道路移动源和非道路移动源.

1.2 估算方法和数据来源

以地级市为基础单元, 对浙江省各排放源的氨排放量进行估算[17], 计算公式如下:

式中, E为排放量; ij分别为地区和排放源; A为活动水平; EF为排放因子; γ为氮-大气氨转换系数, 畜禽养殖取值1.214, 其他行业取值1.0.

文中所涉及的活动水平数据来自政府官方统计资料, 其中畜禽养殖、固氮植物、人体排放、工业生产、道路移动源和非道路移动源的活动水平取自各地市统计年鉴, 氮肥施用、土壤本底、生物质燃烧、废物处理的活动水平取自文献[18], 燃料燃烧的活动水平取自文献[19], 氨排放因子除个别单独注释引用外, 其余均取自文献[17].

1.2.1 畜禽养殖

畜禽养殖排放的氨主要由动物排泄物所释放, 是大气NH3排放的重要来源之一, 约占总排放量的47.5%~58%[20, 21].粪便管理包括户外、圈舍内、粪便存储处理和后续施肥共4个阶段, 各阶段的总铵态氮量依据文献[17]进行计算.根据浙江省养殖情况, 养殖方式主要分为散养和集约化, 畜禽种类主要以常见的饲养品种为计算口径, 分为: 母猪、生猪、奶牛、肉牛、绵羊、肉羊、兔、肉兔、蛋鸡、蛋鸭、肉鸡、肉鸭和肉鹅这13类.对于饲养周期大于1 a的畜禽, 以“年末存栏数”作为活动水平; 饲养周期小于1 a的畜禽, 以“年内出栏数”作为活动水平.散养时畜禽排泄物在室内户外各占50%, 集约化养殖条件下畜禽排泄物在室内户外分别占100%和0[17].猪、牛、家禽和羊的非集约化比例分别为10%、5%、5%和50%[22], 考虑兔的养殖环境和养殖方式与家禽相似, 其非集约化比例以5%计.因文献[17]未明确兔类在不同排泄物管理阶段和不同粪便形态下的排放系数, 因此兔和肉兔的氨排放因子参考GAINS模型[23], 通过物质流循环的畜禽排放计算方法获得其氨排放因子分别为0.58 kg·(d·只)-1和0.24 kg·(d·只)-1.浙江省家禽饲养类别中, 鸡、鸭和鹅分别约占58%、32%和4%, 蛋鸭和蛋鸡存栏比为6∶4[24].单位畜禽排泄量、含氮量和铵态氮比例见表 1.

表 1 畜禽粪便排泄物铵态氮量的估算参数1) Table 1 Parameters used in estimating annual TAN excretion per animal for each livestock class

1.2.2 氮肥施用

氮肥施用后会发生一系列物理化学反应, 以气态氨的形式从土壤挥发到大气造成很大的氮亏损, 成为大气氨的第二个重要来源[26].浙江省在农用化肥方面施用的种类为氮肥、磷肥、钾肥和复合肥, 涉及氨排放的主要为氮肥和复合肥.为进一步提高清单精度, 实现排放因子本地化, 从文献[27]中“各地区农作物化肥投入情况”获悉, 浙江省施用氮肥和复合肥的农作物主要为稻谷(包括早籼稻、粳稻、晚籼稻)和油菜籽, 施用的氮肥种类主要包括尿素、碳铵和其他氮肥; 复合肥主要包括二铵和三元素复合肥.

根据各地市氮肥和复合肥施用量, 并结合文献[27]中农作物不同氮肥的每亩折纯用量及其年种植面积, 计算出研究期历年氮肥和复合肥各成分的施用比例, 结果如表 2所示.浙江省土壤类型以酸性为主, 因此采用校正后的本地化排放因子, 尿素、碳铵、其他氮肥和二铵、三元素复合肥的排放因子分别为38.8、27.5、1.8和3.2、7.4 kg·t-1 [15], 其中二铵和三元素复合肥的氮挥发率与硝铵和硫铵相差较小[28], 排放因子分别引用文献[15]中的硝铵和硫铵.

表 2 浙江省历年不同化肥施用份额/% Table 2 Different fertilizer application shares in Zhejiang province over the years/%

1.2.3 土壤本底

土壤中的生物参与有机物分解排泄出氨或易于水解成氨的氮类化合物, 产生一定量的氨挥发[29].选用各地市耕地面积作为活动水平, 排放因子为180 kg·km-2.

1.2.4 固氮植物

浙江省广泛种植的固氮作物主要包括大豆、花生和绿肥这3类, 选用播种面积作为活动水平.排放因子分别为105、120和135 kg·km-2.

1.2.5 人体排放

人体主要通过排泄物排放氨, 排放量与人口数量和处理条件密切相关[30].考虑到农村地区卫生处理设施存在不完善, 以各地市统计年鉴中的农村人口数量和农村卫生厕所普及率作为活动水平, 排放因子取值0.787 kg·人-1.

1.2.6 生物质燃烧

生物质燃烧包含秸秆灶膛燃烧、薪柴燃烧、秸秆露天焚烧(田间堆肥)、草原火灾和森林火灾.考虑浙江省的地形地貌及植被类型, 氨排放不涉及草原火灾.自2014年以来, 浙江省大力推进秸秆“五化”利用, 杜绝露天焚烧, 至2016年浙江省秸秆年综合利用率达92%[31], 由此将秸秆处置利用划分为两个阶段, 分别为露天焚烧和田间堆肥.

(1) 2008~2015年, 主要估算秸秆露天焚烧量, 计算公式如下:

式中, E为秸秆露天焚烧的氨排放量; ij分别为地区和农作物种类; P为农作物产量; N为谷草比; B为秸秆干燥比; F为秸秆露天焚烧的燃烧率; R为农作物露天焚烧比例; EF为秸秆露天焚烧的氨排放因子.相关计算参数见表 3.

表 3 秸秆露天焚烧氨排放的估算参数[28, 31, 35] Table 3 Parameters used for emission estimation from straw burning

(2) 2016~2018年, 主要估算秸秆田间堆肥量, 计算公式如下:

式中, E为秸秆田间堆肥的氨排放量; P为农作物产量; N为谷草比; R为秸秆堆肥比例, 取值6.40%[32]; EF为秸秆堆肥的氨排放因子, 取值为0.32 kg·t-1.

(3) 森林火灾的生物质燃烧量, 计算公式如下:

式中, E为森林火灾的氨排放量; A为森林火灾受害面积; B为森林的平均地上生物量, 取值为221.94 t·hm-2[33]; η为森林的燃烧效率, 取值为0.33[34]; EF为森林火灾的氨排放因子, 取值为2.9 g·kg-1[33].

1.2.7 工业生产

氮肥生产和合成氨是主要的工业氨排放行业.根据全国排污许可证管理信息平台(http://permit.mee.gov.cn/)中的氮肥生产和合成氨企业情况, 明确全省工业过程氨排放源分布.合成氨排放因子为0.01 kg·t-1, 氮肥生产为5.0 kg·t-1.

1.2.8 废物处理

废物处理主要分为污水处理、垃圾填埋、垃圾焚烧、垃圾堆肥和烟气脱硝等.由火力发电、钢铁、水泥等工业企业生产过程中燃煤消耗引起的烟气脱硝设施的氨排放量已包含在“燃料燃烧”的工业源中, 因此未予重复计算.氨排放因子分别为污水(0.003 g·m-3)、填埋(0.56 kg·t-1)、堆肥(1.275 kg·t-1)和焚烧(0.21 kg·t-1).

1.2.9 燃料燃烧

燃料燃烧主要为工业和生活消耗的煤炭、油品和天然气等燃烧而产生的氨排放.根据文献[19]中“浙江能源平衡表”获取上述3种燃料的工业和生活消耗量, 其中工业燃料消耗量包括终端消费量和用于加工转换其他能源的消费量, 生活燃料消耗量只涉及终端消费量.由于各地市统计年鉴对煤炭、油品和天然气等能源消耗数据较为有限, 只统计了规模以上的企业, 无法反映全行业全领域能源消耗水平, 因此各地市工业燃料消耗量通过全省活动水平数据, 结合文献[36]中“各市国民经济主要指标”的“各市工业生产总值”按占全省比例进行折算.各地市生活燃料消耗量根据其年末常住人口数占全省比例进行折算.燃料燃烧的氨排放因子如表 4.

表 4 燃料燃烧源氨排放因子[37] Table 4 Ammonia emission factors from fuel combustion sources

1.2.10 道路移动源

机动车排放的氨由尾气催化净化装置将氮氧化物转化为N2和O2过程中的二次产物产生[38].相比于人为源的氨排放总量, 机动车产生的氨排放虽然相对较小, 但随着经济发展和机动车保有量的逐年攀升, 机动车尾气的氨排放较其他排放源呈现更为明显的增长趋势.机动车流量与NH3浓度的日变化特征相似, 对NH3早高峰的贡献率可达40%[39].本文将机动车分为大型客车、小型客车、大型货车、小型货车和摩托车这5大类, 历年保有量由各地市统计年鉴获得.排放因子和年均行驶里程见表 5.

表 5 各车型年均行驶里程和排放因子[40] Table 5 Annual average mileages for various types of vehicles and emission factors

1.2.11 非道路移动源

非道路移动源包括工程机械、农业机械、小型通用机械、柴油发电机组、船舶、铁路内燃机车和飞机等.考虑到浙江省作为海洋大省, 海洋和内陆港口较为发达, 拥有宁波舟山港、嘉兴港、温州港和台州港, 以及体现与海洋港口联动发展的义乌国际陆港、湖州港、金华兰溪港、衢州港和丽水青田港等内河(陆)港, 因此本文主要将沿海和内河船舶的氨排放量作为非道路移动源的典型, 其他形式的非道路移动源由于机械保有量较小, 日常使用频率较低, 其氨排放量较小, 暂不予以估算.

对于内河和沿海船舶氨排放量, 基于燃油消耗量计算, 计算公式如下:

式中, E为船舶的氨排放量; Y为船舶燃油消耗量; EF为船舶的氨排放因子, 货船取值为0.8 mg·kg-1[41], 客船取值为4.6 mg·kg-1[42]; Z为船舶客运周转量; Z为船舶货物周转量; YX为船舶油耗系数, 推荐取50[43].

2 结果与讨论 2.1 浙江省氨排放年际变化趋势

根据上述研究方法和活动水平数据, 计算获得2008~2018年浙江省人为源氨排放量的年际变化趋势和各排放源贡献率, 如表 6图 1所示.可以看出, 氨排放量总体呈下降趋势, 从2008年的162.74 kt下降到2018年的108.52 kt, 年均下降率约3.97%.其中, 2008~2012年间氨排放量比较稳定, 年均值为165.97 kt, 介于162.74~166.58 kt; 2013~2018年间则逐年明显下降, 年均下降率6.89%, 最大年下降率9.52%, 出现在2014年.农业源是氨排放的主导因素, 约占历年氨排放总量的82.96%~91.22%, 从2008年的148.20 kt下降到2018年的90.02 kt, 其中畜禽养殖由117.72 kt下降至67.43 kt, 年均下降率5.42%, 占氨排放总量62.14%~73.74%, 氮肥施用由19.68 kt下降至16.84 kt, 占氨排放总量11.85%~15.52%.非农业源的氨排放量波动较小, 介于14.54~18.50 kt.废物处理、道路移动源因污水和垃圾无害化处理能力提升及民用车辆保有量提高, 氨排放量逐年上升, 分别从2.96 kt增长至5.11 kt和2.43 kt增长至9.86 kt.人体排放由4.91 kt减少至0.10 kt, 这与农村卫生厕所普及率的大幅度提高有关, 从2008年的83.75%提高至2018年的99.65%.非道路移动源的贡献最小, 介于1.31~3.55 t, 可见船舶氨排放量在各类排放源中影响很小. 2014年是人为源氨排放量明显下降的转折年, 较2013年下降14.96 kt, 其中畜禽养殖下降了14.15 kt, 分析原因主要是2013年底, 浙江省作出了“五水共治”的决策部署, 打响畜禽养殖污染治理攻坚战, 在2014年, 全省有畜牧业的87个县(市、区)全部重新调整划定禁、限养区, 关停搬迁禁、限养区养殖场7.46万个, 对年存栏50头以上规模畜禽养殖场开展污染整治, 因此导致畜禽养殖氨排放量明显减少.

表 6 2008~2018年浙江省各排放源氨排放量清单/t Table 6 Ammonia emission inventory of various sources in Zhejiang province from 2008 to 2018/t

图 1 2008~2018年浙江省氨排放量变化趋势及各排放源氨排放贡献率 Fig. 1 Trend of ammonia emissions and contribution for different sources in Zhejiang province from 2008 to 2018

2.2 浙江省各地市氨排放年际变化趋势

图 2描述了2008~2018年各地市氨排放量变化趋势.结果表明, 各地市氨排放量呈现先稳定后下降的趋势, 2014年受政府治水治污政策影响, 氨排放量明显下降.杭州市和嘉兴市是浙江省历年氨排放量最大的地市, 年平均氨排放量分别为20.18 kt和19.45 kt.嘉兴市从2008年的25.08 kt下降到2018年的12.87 kt, 在2012~2016年呈现急剧下降态势, 年均下降率达15.29%, 主要是由于在此期间大量“散乱污”畜禽养殖场关闭, 畜禽养殖规模快速减小, 年末生猪存栏数从2012年的273.10万头下降至2016年的18.34万头, 畜禽养殖氨排放量由20.09 kt下降到8.17 kt.杭州市氨排放量表现为先升后降趋势, 从2008年的20.82 kt上升至2013年的23.07 kt, 随后又逐渐下降至2018年的15.98 kt, 引起这一波动主要为畜禽养殖排放源的变化.丽水市和舟山市氨排放量较小且较为稳定, 分别介于5.75~9.29 kt和0.70~1.40 kt.宁波市、温州市、湖州市、绍兴市、金华市、衢州市和台州市氨排放量较为接近, 年平均值在11.72~15.52 kt.

图 2 2008~2018年浙江省不同地市氨排放量变化趋势 Fig. 2 Trends of ammonia emissions for different cities in Zhejiang province from 2008 to 2018

2.3 2018年浙江省及各地市氨排放特征

2018年浙江省氨排放清单和排放强度如表 7图 3所示.全年氨排放量为108.52 kt, 其中农业源氨排放量为90.02 kt, 非农业源为18.50 kt, 分别占氨排放总量的82.96%和17.04%.农业源中的畜禽养殖、氮肥施用为最大贡献源, 氨排放量为67.43 kt和16.84 kt, 分别占氨排放总量的62.14%和15.52%, 两者共占农业源93.61%.非农业源中的道路移动源是第三大贡献源, 氨排放量为9.86 kt, 占比9.09%, 其中客车氨排放量为9.35 kt、货车为0.36 kt、摩托车为0.15 kt.由于现有民用车辆以燃油为主, 机动车尾气排放较大, 因此大力发展新能源汽车替代燃油大小型客车, 对氨减排具有重要意义.废物处理是非农业源氨排放的另一个重要来源, 占比4.71%.非道路移动源的氨排放量相较其他排放源最小, 仅为3.55 t, 集中在舟山市和宁波市, 分别为1.33 t和1.25 t.

表 7 2018年浙江省人为源氨排放清单/t Table 7 Anthropogenic ammonia emission inventory in Zhejiang province in 2018/t

图 3 2018年浙江省各地市各排放源氨排放量及氨排放强度 Fig. 3 Different ammonia emissions and emission intensities for each city in Zhejiang province in 2018

从各地市氨排放情况看, 杭州市、嘉兴市和温州市高于其他城市, 分别占全省氨排放总量14.72%、11.86%和11.80%.舟山市受土地面积制约且多为海岛地形, 各种产业规模较小, 其氨排放量为全省最小(0.70 kt), 仅占0.64%.农业源是主要贡献源, 占71.31%~93.88%.畜禽养殖氨排放最大的是杭州市和衢州市, 分别为9.69 kt和9.34 kt, 主要是由于生猪年内出栏数较大所致, 分别达219.23万头和176.58万头.杭州市的废物处理氨排放量居全省第一, 达1.41 kt, 为宁波市的近2倍, 原因是在污水处理量接近情况下(分别为72 987×104m3和63 241×104 m3), 受不同垃圾填埋量所引起, 杭州市的垃圾填埋量是宁波市的近10倍, 分别为177.53×104 t和17.50×104 t.工业生产氨排放主要分布在杭州市、衢州市、绍兴市和金华市, 氨排放量分别达481.66、189.70、95.37和66.00 t.

从氨排放强度看, 嘉兴市最高, 为3.05 t·km-2, 其次为湖州市、衢州市、绍兴市和温州市, 分别为1.53、1.33、1.07和1.06 t·km-2, 高于全省平均值(1.03 t·km-2), 丽水市最低, 为0.33 t·km-2.各地市氨排放强度由氨排放总量和土地面积构成, 嘉兴市氨排放量居全省首位, 但其土地面积仅占全省4.00%, 因此氨排放强度较大.杭州市虽然具有较大的氨排放量, 但土地面积占全省15.95%, 氨排放强度总体不大, 仅为0.95 t·km-2.丽水市氨排放量较小, 而土地面积最大, 占全省16.35%, 氨排放强度最小(0.33 t·km-2).舟山市是全省土地面积最小的地市, 仅1 459 km2, 氨排放强度较小(0.48 t·km-2).

2.4 2018年主要贡献源氨排放特征

畜禽养殖是浙江省氨排放最大贡献源, 在各地市的氨排放总量占比达到40.57%~79.45%, 随着全省粪污综合利用及处置率和规模化畜禽养殖场粪污处理设施配套比例要求的不断提高, 畜禽养殖的氨减排存在较大空间, 故本文针对不同畜禽种类的氨排放贡献作进一步分析, 结果如图 4所示.可见, 禽类养殖是最大的畜禽养殖排放贡献源, 氨排放量达29.08 kt, 占畜禽养殖氨排放总量43.12%, 其中蛋禽占比33.79%, 肉禽占比9.33%.浙江省蛋禽养殖数量较肉禽少, 但由于蛋禽养殖周期以365 d计, 长于肉禽的50~70 d, 因此蛋禽氨排放量较肉禽高.从养殖种类看, 蛋鸭氨排放贡献率最大, 占畜禽养殖氨排放总量19.61%, 其次为生猪(19.38%)、蛋鸡(14.18%)、肉牛(11.16%)、绵羊(9.99%)和肉羊(8.47%), 共占畜禽养殖氨排放总量82.79%, 母猪、奶牛和兔类氨排放贡献率较低, 均在3%以下.猪肉、牛肉和禽蛋是浙江省居民日常生活食用的必需品, 具有较大的饲养规模, 氨排放贡献也较高.相比肉羊, 肉牛年内出栏数为8.15万头, 远少于肉羊的135.06万头, 但由于肉牛粪尿等排泄物的铵态氮量[99.6 g·(d·头)-1]远高于肉羊[12.10 g·(d·头)-1], 因此肉牛仍具有较高的氨排放量.在实现畜禽养殖氨减排途径方面, 宜推广应用低蛋白日粮及粪污密闭或覆盖贮存、种养结合与粪便深施等技术[44], 同时合理控制猪、牛和羊等大牲畜的养殖规模, 实行宜养区以地定畜措施, 将对氨排放的有效减少发挥较大作用.

图 4 2018年不同畜禽养殖种类氨排放贡献率 Fig. 4 Contribution of ammonia emissions from various types of livestock in Zhejiang province in 2018

氮肥施用是浙江省氨排放第二大贡献源, 在各地市的氨排放总量占比达到10.08%~26.26%, 图 5为各地市不同氮肥类型的氨排放情况.尿素是各地市氮肥施用的主要贡献源, 占氮肥氨排放总量的82.59%~94.33%, 施用量达38.15×104 t.碳铵贡献率最小, 占3.05%~3.50%, 主要是因为浙江省年均气温较高, 碳铵的分解挥发损失量比尿素大, 同时施肥过程中需深施等多种原因, 致使其施用范围较小, 施用量较低(1.99×104 t).嘉兴市和绍兴市因尿素产生的氨排放量最高, 分别为2 683.23 t和2 101.10 t, 由三元素复合肥产生的氨排放最大的为宁波市和金华市, 分别为238.05 t和236.45 t.舟山市农业种植面积有限, 氮肥施用产生的氨排放量在全省中最小(95.08 t).在实现氮肥施用氨减排途径方面, 建议推广测土配方施肥和新型肥料应用等氮肥减量技术与模式.

图 5 2018年浙江省各地市不同氮肥类型氨排放量及贡献率 Fig. 5 Ammonia emissions and contribution from different fertilizer types for each city in Zhejiang province in 2018

2.5 2018年浙江省氨排放空间分布特征

本文基于ArcGIS空间分析技术, 综合考虑人口空间分布、农村居民点分布、土地利用类型、GDP空间分布和道路网等空间特征标准数据, 利用多因子权重分配法, 将各地市为基本统计单元的不同活动水平氨排放量展布到栅格单元, 建立了1 km×1 km网格化人为源氨排放清单, 相关栅格数据来源于中国科学院地理科学与资源研究所资源环境科学数据中心(https://www.resdc.cn/), 分配结果如图 6所示.从中可知, 浙江省人为源氨排放空间分布总体较为分散, 主要在浙江省北部区域, 呈现出“北高南低”的排放趋势, 高值区集中于环杭州湾带, 分布在嘉兴市、湖州市东部、杭州市东北部、绍兴市北部和宁波市中北部, 单位网格氨排放量大于10.00 t.台州市和温州市的沿海区域、金华市和衢州市的中部地区也存在部分高值区, 舟山市和丽水市高值区较少, 仅为零星分布.高值区的分布与浙江省杭嘉湖、宁绍等平原区域一致, 这是由于平原地区主要是粮食作物的主产区与畜禽养殖区, 而城市中心、丘陵和山地等区域农业活动水平较少, 因此氨排放量较小.

图 6 2018年浙江省氨排放空间分布(1 km×1 km) Fig. 6 Spatial distribution of ammonia emissions in Zhejiang province in 2018

2.6 不同地区氨排放清单对比

本文对比了与其他学者在江苏[40]、山东[45]和安徽[46]等地区相近基准年的氨排放清单(表 8).结果表明, 畜禽养殖和氮肥施用是各地区氨排放主要贡献源, 基本维持在氨排放总量的78.31%~85.90%.本文估算的氨排放量相较这些地区均处于较低水平, 分别是其氨排放量的20.45%、15.17%和17.11%.差异较大的氨排放源主要为氮肥施用和畜禽养殖, 其中氮肥施用除了与当地年施用量有关外, 还与所选用的各类型氮肥的排放因子有关, 本文基于浙江省酸性土壤条件选择的氨排放因子远小于其他中性或碱性土壤地区, 例如文献[40]中江苏省尿素和碳铵的氨排放因子分别为21.3%和17.4%, 是本文的5.49倍和6.33倍.畜禽养殖则主要与各地区不同的养殖品种及规模大小直接相关.总体看来, 氨排放量的高低体现了不同地区间经济发展水平、人口规模、物质需求、农牧业结构和自然条件等的差异.

表 8 与其他地区氨排放研究结果比较/kt Table 8 Comparison of ammonia emissions with other regions/kt

对比相同基准年的浙江省氨排放结果, 本文估算的2017年浙江省氨排放量为115.04 kt, 与赵睿东等[16]的研究结果相近, 相对偏小6.96 kt, 差异较大的氨排放源主要体现在氮肥施用、畜禽养殖、人体排放和废物处理等这4种.其中, 氮肥施用相比偏小13.40 kt, 主要是由于本文将各类型氮肥的施用份额实现本地化, 且所选用的校正排放因子也存在差异.畜禽养殖相比高出36.55 kt, 是因为本文在估算氨排放量较大的畜禽种类时, 考虑到出栏的肉牛和肉羊饲养周期大于1 a[47], 因此将其饲养周期以365 d计, 从而导致畜禽养殖氨排放量较大.人体排放相比偏小13.84 kt, 则是由于本文在同样基于农村人口基础上, 额外考虑了农村卫生厕所普及率, 使其更加符合实际, 浙江省2017年农村卫生厕所普及率已达98.64%[18].废物处理相比偏小16.48 kt, 是因为本文对烟气脱硝相关工业企业燃煤产生的氨排放已在“燃料燃烧源”中予以计算, 因此未对烟气脱硝过程再次估算, 且燃料燃烧源中工业燃煤消耗的排放因子(0.014 kg·t-1)较烟气脱硝(0.155 kg·t-1)相差一个数量级, 致使同等煤炭消耗量所计算的烟气脱硝氨排放量明显高于燃料燃烧源.

2.7 不确定性分析

本文假设氨排放估算过程中的活动水平和排放因子符合对数正态分布, 采用蒙特卡罗模拟方法, 以各类排放源活动水平和排放因子为参数, 选用Oracle Crystal Ball软件随机抽样运算10 000次, 定量分析了在95%置信区间下的2018年浙江省氨排放不确定度, 如表 9所示.结果表明, 模拟的氨排放总量平均值为108.37 kt, 略低于本文估算的氨排放量0.15 kt, 不确定度范围为-5.40%~5.60%, 各排放源的不确定度在-18.25%~21.13%之间.氮肥施用平均氨排放量为16.84 kt, 与估算值相同, 不确定度为-16.33%~18.48%.畜禽养殖和燃料燃烧的不确定度较小, 分别为-6.82%~7.23%和-8.37%~8.97%, 主要是由于本文估算的畜禽养殖氨排放考虑了浙江省比较全面的共13类畜禽饲养种类和室内散养、室外散养和集约化等3种饲养方式以及户外、圈舍、粪便存储处理、施肥这4个阶段的畜禽粪便管理; 燃料燃烧氨排放则对煤炭、油品和天然气等3种燃料区分了工业和生活两个方面的消费情况, 总体看来本文较为全面地细分了活动水平.对比长三角[12] (-55%~60%)、京津冀[47] (-17%~27%)和青海东部城市群[48](-9.8%~17.1%)等其他研究区域, 本文模拟结果较好, 所建立的清单不确定度明显下降.

表 9 95%置信区间下的2018年浙江省氨排放不确定度/% Table 9 Uncertainty of ammonia emissions under 95% confidence interval in Zhejiang province of 2018/%

3 结论

(1) 2008~2018年浙江省人为源氨排放量总体呈现下降趋势, 从2008年的162.74 kt下降到2018年的108.52 kt, 年均下降率为3.97%.农业源是浙江省氨排放的主要贡献源, 占氨排放总量的82.96%~91.22%; 非农业源比较稳定, 氨排放总量介于14.54~18.50 kt.杭州市和嘉兴市是历年氨排放量最大的地市, 年平均分别为20.18 kt和19.45 kt.

(2) 2018年浙江省农业源氨排放量为90.02 kt, 非农业源18.50 kt.杭州市、嘉兴市和温州市的氨排放量高于其他城市, 分别占全省氨排放总量的14.72%、11.86%和11.80%.浙江省氨排放强度为1.03 t·km-2, 最高的是嘉兴市, 为3.05 t·km-2.禽类养殖氨排放量达到29.08 kt, 占畜禽养殖氨排放总量的43.12%.尿素氨排放量约14.80 kt, 占各地市氮肥氨排放总量的82.59%~94.33%.

(3)从空间分布特征看, 浙江省人为源氨排放主要分布在北部区域, 呈现“北高南低”的排放趋势, 高值区集中在嘉兴市、湖州市东部、杭州市东北部、绍兴市北部和宁波市中北部, 舟山市和丽水市高值区较少.

(4)不确定性结果显示, 模拟的2018年浙江省氨排放总量平均值为108.37 kt, 在95%置信区间的不确定度范围为-5.40%~5.60%, 各排放源的不确定度在-18.25%~21.13%之间, 畜禽养殖和燃料燃烧的不确定度较小.

参考文献
[1] Ayers G P, Gras J L. Ammonia gas concentrations over the Southern Ocean[J]. Nature, 1980, 284(5756): 539-540. DOI:10.1038/284539a0
[2] Chen S H, Cheng M M, Guo Z, et al. Enhanced atmospheric ammonia (NH3) pollution in China from 2008 to 2016: evidence from a combination of observations and emissions[J]. Environmental Pollution, 2020, 263. DOI:10.1016/j.envpol.2020.114421
[3] Behera S N, Sharma M. Investigating the potential role of ammonia in ion chemistry of fine particulate matter formation for an urban environment[J]. Science of the Total Environment, 2010, 408(17): 3569-3575. DOI:10.1016/j.scitotenv.2010.04.017
[4] Krupa S V. Effects of atmospheric ammonia (NH3) on terrestrial vegetation: a review[J]. Environmental Pollution, 2003, 124(2): 179-221. DOI:10.1016/S0269-7491(02)00434-7
[5] 王文兴, 卢筱凤, 庞燕波, 等. 中国氨的排放强度地理分布[J]. 环境科学学报, 1997, 17(1): 2-7.
Wang W X, Lu X F, Pang Y B, et al. Geographical distribution of NH3 emission intensities in China[J]. Acta Scientiae Circumstantiae, 1997, 17(1): 2-7. DOI:10.3321/j.issn:0253-2468.1997.01.001
[6] Gu B J, Zhang L, Van Dingenen R, et al. Abating ammonia is more cost-effective than nitrogen oxides for mitigating PM2.5 air pollution[J]. Science, 2021, 374(6568): 758-762. DOI:10.1126/science.abf8623
[7] Anderson N, Strader R, Davidson C. Airborne reduced nitrogen: ammonia emissions from agriculture and other sources[J]. Environment International, 2003, 29(2-3): 277-286. DOI:10.1016/S0160-4120(02)00186-1
[8] 薛文博, 许艳玲, 唐晓龙, 等. 中国氨排放对PM2.5污染的影响[J]. 中国环境科学, 2016, 36(12): 3531-3539.
Xue W B, Xu Y L, Tang X L, et al. Impacts of ammonia emission on PM2.5 pollution in China[J]. China Environmental Science, 2016, 36(12): 3531-3539. DOI:10.3969/j.issn.1000-6923.2016.12.002
[9] Hamaoui-Lague L, Meleux F, Beekmann M, et al. Improving ammonia emissions in air quality modelling for France[J]. Atmospheric Environment, 2014, 92: 584-595. DOI:10.1016/j.atmosenv.2012.08.002
[10] Megaritis A G, Fountoukis C, Charalampidis P E, et al. Response of fine particulate matter concentrations to changes of emissions and temperature in Europe[J]. Atmospheric Chemistry and Physics, 2013, 13(6): 3423-3443. DOI:10.5194/acp-13-3423-2013
[11] Fu H, Luo Z B, Hu S Y. A temporal-spatial analysis and future trends of ammonia emissions in China[J]. Science of the Total Environment, 2020, 731. DOI:10.1016/j.scitotenv.2020.138897
[12] Yu X N, Shen L, Hou X H, et al. High-resolution anthropogenic ammonia emission inventory for the Yangtze River Delta, China[J]. Chemosphere, 2020, 251. DOI:10.1016/j.chemosphere.2020.126342
[13] 王琛, 尹沙沙, 于世杰, 等. 河南省2013年大气氨排放清单建立及分布特征[J]. 环境科学, 2018, 39(3): 1023-1030.
Wang C, Yin S S, Yu S J, et al. A 2013-based atmospheric ammonia emission inventory and its characteristic of spatial distribution in Henan Province[J]. Environmental Science, 2018, 39(3): 1023-1030.
[14] 浙江省发展和改革委员会. 浙江省生态环境保护"十四五"规划[EB/OL]. http://fzggw.zj.gov.cn/art/2021/7/12/art_1229123366_2310857.html, 2021-07-12.
[15] 余飞翔, 晁娜, 吴建, 等. 浙江省2013年农业源氨排放清单研究[J]. 环境污染与防治, 2016, 38(10): 41-46.
Yu F X, Chao N, Wu J, et al. Research on agricultural ammonia emission inventory of Zhejiang Province in 2013[J]. Environmental Pollution & Control, 2016, 38(10): 41-46.
[16] 赵睿东, 于兴娜, 侯新红, 等. 浙江省人为源氨排放清单建立及分布特征[J]. 环境科学, 2020, 41(9): 3976-3984.
Zhao R D, Yu X N, Hou X H, et al. Establishment and spatial distribution of anthropogenic ammonia emission inventory in Zhejiang province[J]. Environmental Science, 2020, 41(9): 3976-3984.
[17] 环境保护部. 大气氨源排放清单编制技术指南(试行)[EB/OL]. https://www.mee.gov.cn/gkml/hbb/bgg/201408/W020140828351293771578.pdf, 2014-08-20.
[18] 浙江省统计局. 浙江自然资源与环境统计年鉴-2009-2019[M]. 北京: 中国统计出版社, 2009-2019.
[19] 国家统计局能源统计司. 中国能源统计年鉴-2009-2019[M]. 北京: 中国统计出版社, 2009-2019.
[20] Xu P, Zhang Y S, Gong W W, et al. An inventory of the emission of ammonia from agricultural fertilizer application in China for 2010 and its high-resolution spatial distribution[J]. Atmospheric Environment, 2015, 115: 141-148. DOI:10.1016/j.atmosenv.2015.05.020
[21] Ye Z L, Guo X R, Cheng L, et al. Reducing PM2.5 and secondary inorganic aerosols by agricultural ammonia emission mitigation within the Beijing-Tianjin-Hebei region, China[J]. Atmospheric Environment, 2019, 219. DOI:10.1016/j.atmosenv.2019.116989
[22] 刘波, 童仪, 李安, 等. 长江三角洲地区畜禽养殖业氨排放清单研究[J]. 生态与农村环境学报, 2018, 34(11): 1042-1049.
Liu B, Tong Y, Li A, et al. Study on ammonia emission inventory of livestock over Yangtze River Delta[J]. Journal of Ecology and Rural Environment, 2018, 34(11): 1042-1049. DOI:10.11934/j.issn.1673-4831.2018.11.012
[23] Klimont Z, Winiwarter W. Integrated ammonia abatement-modelling of emission control potentials and coats in GAINS[R]. IR-11-027. Laxenburg, Austria: International Institute for Applied Systems Analysis (ⅡASA), 2011.
[24] 俞国乔, 何世山, 宋美娥. 浙江省家禽业发展现状及对策建议[J]. 浙江农业科学, 2014(12): 1922-1924. DOI:10.3969/j.issn.0528-9017.2014.12.052
[25] 杨志鹏. 基于物质流方法的中国畜牧业氨排放估算及区域比较研究[D]. 北京: 北京大学, 2008.
Yang Z P. Estimation of ammonia emission from livestock in China based on mass-flow method and regional comparison[D]. Beijing: Peking University, 2008.
[26] Harrison R, Webb J. A review of the effect of N fertilizer type on gaseous emissions[J]. Advances in Agronomy, 2001, 73: 65-108.
[27] 国家发展和改革委员会价格司. 全国农产品成本收益资料汇编-2009-2019[M]. 北京: 中国统计出版社, 2009-2019.
[28] 高榕. 山东省农业源氨排放清单及其减排潜力研究[D]. 济南: 济南大学, 2018.
Gao R. Inventories of the ammonia emission from agriculture and the potential mitigation for future emission[D]. Ji'nan: University of Jinan, 2018.
[29] 覃思, 孔少飞, 吴剑, 等. 1996-2016年湖北省氨排放时空差异及影响因素[J]. 中国环境科学, 2020, 40(4): 1403-1413.
Qin S, Kong S F, Wu J, et al. Spatial-temporal diversities of ammonia emissions and impacting factors in Hubei Province from 1996 to 2016[J]. China Environmental Science, 2020, 40(4): 1403-1413. DOI:10.3969/j.issn.1000-6923.2020.04.002
[30] Sutton M A, Dragosits U, Tang Y S, et al. Ammonia emissions from non-agricultural sources in the UK[J]. Atmospheric Environment, 2000, 34(6): 855-869. DOI:10.1016/S1352-2310(99)00362-3
[31] 裘一冰, 叶波, 陆若辉. 浙江省农作物秸秆饲料化利用现状的分析与探讨[J]. 浙江农业科学, 2018, 59(2): 331-332.
[32] 李飞, 邵建均. 浙江省农作物秸秆资源台账制度建设现状及对策建议[J]. 浙江农业科学, 2020, 61(8): 1640-1641, 1646.
Li F, Shao J J. Status of agricultural resources accounting construction and its countermeasures[J]. Journal of Zhejiang Agricultural Sciences, 2020, 61(8): 1640-1641, 1646.
[33] 刘元隆. 福建省人为源氨排放清单研究[D]. 厦门: 厦门大学, 2016.
Liu Y L. Research on ammonia emission inventory in Fujian Province[D]. Xiamen: Xiamen University, 2016.
[34] Ito A, Penner J E. Global estimates of biomass burning emissions based on satellite imagery for the year 2000[J]. Journal of Geophysical Research: Atmospheres, 2004, 109(D14): D14S05.
[35] 彭立群, 张强, 贺克斌. 基于调查的中国秸秆露天焚烧污染物排放清单[J]. 环境科学研究, 2016, 29(8): 1109-1118.
Peng L Q, Zhang Q, He K B. Emissions inventory of atmospheric pollutants from open burning of crop residues in China based on a national questionnaire[J]. Research of Environmental Sciences, 2016, 29(8): 1109-1118.
[36] 浙江省统计局, 国家统计局浙江调查总队. 浙江统计年鉴-2009-2019[M]. 北京: 中国统计出版社, 2009-2019.
[37] 沈兴玲, 尹沙沙, 郑君瑜, 等. 广东省人为源氨排放清单及减排潜力研究[J]. 环境科学学报, 2014, 34(1): 43-53.
Shen X L, Yin S S, Zheng J Y, et al. Anthropogenic ammonia emission inventory and its mitigation potential in Guangdong Province[J]. Acta Scientiae Circumstantiae, 2014, 34(1): 43-53.
[38] Heeb N V, Forss A M, Brühlmann S, et al. Three-way catalyst-induced formation of ammonia-velocity-and acceleration-dependent emission factors[J]. Atmospheric Environment, 2006, 40(31): 5986-5997. DOI:10.1016/j.atmosenv.2005.12.035
[39] Gu M N, Pan Y P, Walters W W, et al. Vehicular emissions enhanced ammonia concentrations in winter mornings: insights from diurnal nitrogen isotopic signatures[J]. Environmental Science & Technology, 2022, 56(3): 1578-1585.
[40] 侯新红, 于兴娜, 沈丽, 等. 2013-2017年江苏省人为源氨排放清单的建立及特征[J]. 环境科学, 2019, 40(11): 4862-4869.
Hou X H, Yu X N, Shen L, et al. Establishment and characteristics of an artificial ammonia emissions inventory in Jiangsu Province from 2013 to 2017[J]. Environmental Science, 2019, 40(11): 4862-4869.
[41] 黄学良, 张洲, 杨威强, 等. 广州港船舶停泊工况排放因子实测及排放量初步估算[J]. 环境科学, 2017, 38(8): 3162-3168.
Huang X L, Zhang Z, Yang W Q, et al. Emission factors and preliminary emission estimates of air pollutants from ships at berth in the Guangzhou Port[J]. Environmental Science, 2017, 38(8): 3162-3168.
[42] Cooper D A. Exhaust emissions from high speed passenger ferries[J]. Atmospheric Environment, 2001, 35(24): 4189-4200. DOI:10.1016/S1352-2310(01)00192-3
[43] 环境保护部. 非道路移动源大气污染物排放清单编制技术指南(试行)[EB/OL]. https://www.mee.gov.cn/gkml/hbb/bgg/201501/W020150107594587960717.pdf, 2014-12-31.
[44] 刘学军, 沙志鹏, 宋宇, 等. 我国大气氨的排放特征、减排技术与政策建议[J]. 环境科学研究, 2021, 34(1): 149-157.
Liu X J, Sha Z P, Song Y, et al. China's atmospheric ammonia emission characteristics, mitigation options and policy recommendations[J]. Research of Environmental Sciences, 2021, 34(1): 149-157.
[45] 赵国梁, 常景云, 成杰民, 等. 2008-2018年山东省人为源氨排放清单研究[J]. 环境污染与防治, 2021, 43(2): 259-265.
Zhao G L, Chang J Y, Cheng J M, et al. Anthropogenic ammonia emission inventory in Shandong Province from 2008 to 2018[J]. Environmental Pollution & Control, 2021, 43(2): 259-265.
[46] 王康宏, 沈丽, 赵睿东. 2006-2017年安徽省人为源氨排放清单及特征[J]. 大气科学学报, 2020, 43(3): 547-556.
Wang K H, Shen L, Zhao R D. Anthropogenic ammonia emission inventory and its spatial distribution in Anhui Province from 2006 to 2017[J]. Transactions of Atmospheric Sciences, 2020, 43(3): 547-556.
[47] 叶枝兰. 京津冀地区氨排放历史特征及减排潜力预测研究[D]. 北京: 北京工业大学, 2020.
Ye Z L. Historical characteristics and prediction of mitigation potential of ammonia emissions within the Beijing-Tianjin-Hebei Region, China[D]. Beijing: Beijing University of Technology, 2020.
[48] 张津建, 李广英, 张元勋, 等. 青海东部城市群大气氨排放清单研究[J]. 环境科学研究, 2021, 34(9): 2113-2121.
Zhang J J, Li G Y, Zhang Y X, et al. Ammonia emission inventory in eastern urban agglomeration of Qinghai Province[J]. Research of Environmental Sciences, 2021, 34(9): 2113-2121.