环境科学  2023, Vol. 44 Issue (6): 3130-3141   PDF    
引用本文
王艳妮, 杨敬婷, 黄贤峰, 程燕, 陆标, 顾兆林. 贵州省生物质燃烧源大气污染物排放清单[J]. 环境科学, 2023, 44(6): 3130-3141.
WANG Yan-ni, YANG Jing-ting, HUANG Xian-feng, CHENG Yan, LU Biao, GU Zhao-lin. Emission Inventory of Airborne Pollutants from Biomass Combustion in Guizhou Province[J]. Environmental Science, 2023, 44(6): 3130-3141.
贵州省生物质燃烧源大气污染物排放清单
王艳妮1,2, 杨敬婷2, 黄贤峰2, 程燕1, 陆标3, 顾兆林1     
1. 西安交通大学人居环境与建筑工程学院, 西安 710054;
2. 贵州省环境科学研究设计院, 贵阳 550081;
3. 贵州大学资源与环境工程学院, 贵阳 550025
收稿日期: 2022-07-04; 修订日期: 2022-08-19
基金项目: 中央引导地方科技发展资金项目(黔科中引地[2022]4056); 贵州省基础研究计划项目(黔科合基础-ZK[2022]一般239)
作者简介: 王艳妮(1982~), 女, 高级工程师, 主要研究方向为大气污染源排放清单, E-mail: 215521519@qq.com
通信作者: 顾兆林, E-mail: guzhaoln@mail.xjtu.edu.cn
摘要: 为准确掌握贵州省生物质燃烧源大气污染物的排放状况, 基于收集资料和实地调查结合的方式获取活动水平, 引用文献和本地实测数据结合的方式选取排放系数, 采用排放系数法结合GIS技术, 建立了贵州省2019年3 km×3 km生物质燃烧源9种大气污染物排放清单.结果表明: ①全省生物质燃烧源CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC的排放量分别为: 293505.53、14781.19、4146.11、8501.07、45025.70、39463.58、41879.31、6832.33和15134.74 t.户用生物质炉具CO、SO2、NH3和BC的排放量贡献率最大, 秸秆露天焚烧NOx、VOCs、PM2.5、PM10和OC的排放量贡献率最大. ②各市(州)生物质燃烧源排放的大气污染物分布不均衡, 主要集中于黔东南苗族侗族自治州.各二级源污染物排放量变化特征不同, 月排放集中在2、3、4和12月, 小时排放在每日14:00~15:00最大. ③排放清单仍具有一定不确定性, 后续研究中应进一步提升活动水平数据获取的准确性, 继续开展更细分的排放系数本地化实测, 完善贵州省生物质燃烧源排放清单, 为大气污染协同治理提供科学依据.
关键词: 生物质燃烧源      大气污染物      排放清单      贵州省      不确定性分析     
Emission Inventory of Airborne Pollutants from Biomass Combustion in Guizhou Province
WANG Yan-ni1,2 , YANG Jing-ting2 , HUANG Xian-feng2 , CHENG Yan1 , LU Biao3 , GU Zhao-lin1     
1. School of Human Settlement and Civil Engineering, Xi'an Jiaotong University, Xi'an 710054, China;
2. Guizhou Institute of Environmental Sciences Research and Design, Guiyang 550081, China;
3. College of Resources and Environment Engineering, Guizhou University, Guiyang 550025, China
Abstract: With the objective to ascertain the emissions of biomass combustion in Guizhou, the activity levels were measured through data collection and field surveys, and the emission factors were acquired using actual monitoring data and data cited from previous literature. A 3 km×3 km-gridded emission inventory of nine air pollutants from biomass combustion sources in Guizhou Province in 2019 was developed in combination with GIS technology. The results showed that the total emissions of CO, NOx, SO2, NH3, VOCs, PM2.5, PM10, BC, and OC in Guizhou were estimated to be 293 505.53, 14 781.19, 4 146.11, 8 501.07, 45 025.70, 39 463.58, 41 879.31, 6 832.33, and 15 134.74 t, respectively. The distribution of atmospheric pollutants emitted by biomass combustion sources in different cities was noticeably uneven, being mainly concentrated in Qiandongnan Miao and Dong Autonomous Prefecture. The analysis of variation characteristics of emissions indicated that the monthly emissions were concentrated in February, March, April, and December, and the hourly emissions peaked daily from 14:00 to 15:00. Some uncertainty remained in the emission inventory. It is necessary to perform in-depth analyses of the accuracy of obtaining activity-level data, localize the emission factors through more combustion experiments in subsequent research for improving the emission inventory of air pollutants from biomass combustion in Guizhou Province, and provide a basis for the cooperative governance of the atmospheric environment.
Key words: biomass combustion      air pollutants      emission inventory      Guizhou Province      uncertainly analysis     

在发展中国家, 燃烧生物质是能源利用的一个重要途径[1], 但同时会产生大量污染物, 是主要大气污染来源之一[2~5], 这些污染物对气候变化、区域空气质量和人类健康的影响日益受到关注[6~13], 中国也不例外, 如杨光义[14]对东北地区生物质露天燃烧源排放进行大气环境影响评估, 2015年秋季和春季PM2.5浓度的贡献率分别为52.0%和24.6%, 是东北地区两次重污染事件的一个重要原因; 徐敬等[15]模拟结果显示秸秆燃烧导致华北地区PM2.5浓度月平均值上升幅度普遍超过3 μg·m-3; 河南、山东和天津等地的秸秆燃烧在合适风场的作用下会严重影响北京, 可导致丰台及通州等地PM2.5小时浓度上升超过17 μg·m-3.

与煤炭、石油和天然气等化石能源燃烧以集中燃烧排放为主不同, 生物质燃烧源具有分散性、隐蔽性和开放性等特点, 调查难度远高于化石能源, 制定全面和详细的反映区域生物质燃烧现状的排放清单对于大气污染防治具有重要的价值.近些年针对中国省域或城市大气污染源排放清单研究表明, 生物质燃烧源排放大气污染物占人为源的比例不容忽视, 鄂州市[16]、嘉兴市[17]、天水市[18]、鹤壁市[19]、济南市[20]、兰州市[21]、资阳市[22]、山东省[23]和四川省[24]的CO、VOCs、PM2.5、PM10、BC和OC排放量约占5.07%~78.88%.对于我国生物质燃烧源的排放特征, Shi等[25]估算了东北地区2001~2017年生物质开放燃烧CO2排放量, 秸秆露天焚烧贡献最大(68%), 其次是森林火灾和草原火灾; Zhao等[26]基于卫星观测核算了东北地区生物质燃烧排放的温室气体和气溶胶, 秸秆露天焚烧是其主要来源, 每年的农作物收获后和种植前排放达到峰值; Wu等[27, 28]编制了大陆地区生物质燃烧排放左旋葡聚糖(LG)排放清单, 2014年排放量贡献率为39.20%, 并实测了玉米芯本地化排放情况; Cui等[29]的研究可知大陆地区秸秆露天焚烧亚硝酸(HONO)排放主要集中在华北和东北地区, 在3、4、6和10月达到峰值.目前没有针对贵州省的生物质燃烧排放核算, 田贺忠等[3]研究了中国生物质燃烧大气污染物排放清单, 其中贵州省的活动水平数据来源于统计年鉴、排放因子来源于已有研究成果的平均值, 而“已有研究结果”不包含贵州本地化研究结果.其他学者核算中国生物质燃烧排放时包含贵州省也未开展大样本、实地详细调查[2, 4, 5, 30~34], 无本地化实测数据, 对于贵州省生物质排放源的核算较为粗略, 未见在不确定性方面的敏感度分析, 从时间上基准年最新仅为2017年.

优良生态环境是贵州在新时代西部大开发上闯新路最大的发展优势和竞争优势, 截至2022年3月, 全省9个中心城市环境空气质量平均优良天数比率为96.8%[35].然而, 贵州是一个喀斯特地貌典型发育的山区省份, 非煤炭产区的广大农村地区尤其是喀斯特地区的能源消费以传统生物质能为主[36~42], 秸秆供能占农村生活用能的30.00%~35.00%[40], 另外还有相当一部分就地焚烧, 生物质燃烧源排放的大气污染物不容小觑.贵州省2015年大气污染源排放清单[43]表明生物质燃烧源排放污染物在人为源中占比较大, CO、VOCs、PM2.5、BC和OC高达26.56%、30.10%、23.16%、36.47%和47.45%.目前尚未见针对贵州省的生物质燃烧源排放研究的报道, 因此, 本研究统计调查了贵州省2019年各市(州)生物质燃烧的活动水平数据, 计算建立了9种大气污染物CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC的排放清单, 利用ArcGIS和蒙特卡洛方法进行空间分配和不确定性分析.

1 材料与方法

本文以2019年为基准年, 以贵州省9个市州(贵阳市、安顺市、黔东南苗族侗族自治州、黔南布依族苗族自治州、铜仁市、黔西南布依族苗族自治州、毕节市、六盘水市和遵义市)为研究区域.根据文献[44, 45]中推荐的生物质燃烧源计算方法, 估算得到了生物质锅炉、户用生物质炉具(秸秆、薪柴)和生物质开放燃烧(森林火灾、草原火灾、秸秆露天焚烧)的9项大气污染物CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC排放量.

采用公式(1)计算生物质燃烧产生某一种大气污染物的排放量(Ei).

(1)

式中, A为排放源活动水平, t; EF为排放系数, g·kg-1; i为某一种大气污染物; j为贵州省各市(州); k为生物质燃烧类型(生物质锅炉、户用生物质炉具、森林火灾、草原火灾和秸秆露天焚烧); m为燃料、植被带、草地或秸秆类型.

对于生物质锅炉, 在安装除尘器等污染控制设施情况下, 排放系数EF由公式(2)计算得到：

(2)

式中, EF0为污染物产生系数, g·kg-1; η为污染控制设施的去除效率, %.

生物质锅炉活动水平包括生物质锅炉的地理位置、所用生物质燃料类型、生物质燃料年消费量和污染控制设施的去除效率.户用生物质炉具活动水平主要包括常住人口户数、生物质燃烧用户比例、秸秆及薪柴用户比例、不同类型秸秆比例和生物质燃料户均消费量等.生物质开放燃烧活动水平包括森林火灾受灾面积、草原火灾受灾面积、植被带森林类型、草地类型和农作物产量等.本次源清单编制过程中, 活动水平来源于社区或村庄填报、研究人员抽样入户调查、MODIS火点数据、统计年鉴、环境统计和其他统计数据等, 其具体调查内容和获取途径见表 1.由于六盘水市、毕节市、黔西南布依族苗族自治州和遵义市的煤炭资源储量约占贵州省已发现总资源储量的95.51%[46], 这些市(州)为全省主要产煤区, 居民生活能源以煤炭为主, 因此选择非主要产煤市(州)贵阳市、安顺市、黔东南苗族侗族自治州、黔南布依族苗族自治州和铜仁市作为社区或村庄填报范围来调查生物质使用情况, 填报之前开展培训, 并建立工作群通过网络答疑.选择生物质燃烧用户数占比前两位的黔东南苗族侗族自治州、黔南布依族苗族自治州和省会贵阳市作为研究人员抽样入户调查范围, 采用分层随机抽样, 即样本散布在市(州), 以区(市、县)作为调查样本单元, 以户为单位进行调查.活动水平数据质量控制包括正确性、一致性和完整性检验, 经过多级数据审核、利用SPSS软件分析抽样效果、与宏观数据对比等多种方式, 确保数据真实有效.

表 1 生物质燃烧源活动水平调查内容和数据获取途径 Table 1 Survey of activity level for biomass combustion and data acquiring method

本次清单编制中, 除户用生物质炉具薪柴、辣椒、大豆和玉米秸秆的BC和OC排放系数来自笔者实验室的本地化实测结果(未发表的数据), 各二级源类各污染物排放系数的选取主要来自已有研究, 见表 2.

表 2 生物质燃烧源排放系数/g·kg-1 Table 2 Emission factors for biomass combustion/g·kg-1

2 结果与讨论 2.1 2019年贵州省生物质燃烧大气污染物排放清单分析

根据排放系数法核算, 建立2019年贵州省生物质燃烧大气污染物排放清单(表 3).2019年贵州省生物质燃烧源CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC污染物的总排放量分别为293 505.53、14 781.19、4 146.11、8 501.07、45 025.70、39 463.58、41 879.31、6 832.33和15 134.74 t.目前没有针对贵州省生物质燃烧源大气污染物的排放核算, 有学者核算中国生物质燃烧排放时包含贵州省[2~5, 30~34], 本研究与其他研究结果在污染物排放量上有一定差异(表 4), 这与采用的研究方法、活动水平的数据来源和排放系数的选取途径等有关, 但总体上具有一定可比性.笔者核算贵州省2015年排放清单时包含生物质燃烧源[43], 本研究结果与该生物质燃烧源CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC排放量(分别为508 567.98、17 232.09、6 993.25、5 772.5、64 619.6、65 538.92、67 988.72、10 088.40和26 143.51 t)的比值分别是0.58、0.86、0.59、1.47、0.70、0.60、0.62、0.68和0.58, 主要原因是2019年排放清单开展了生物质用户比例和户均生物质燃烧量等实地调查, 降低了参与计算的总生物质燃烧量, 户用生物质炉具秸秆和薪柴的使用情况更为精准, 以及森林和草原火灾受灾面积比2015年大幅度减少.

表 3 2019年贵州省生物质燃烧源排放清单/t Table 3 Emission inventory of Guizhou Province in 2019/t

表 4 贵州省生物质燃烧排放结果对比/t Table 4 Comparison of results of biomass combustion emission in Guizhou Province/t

2019年贵州省生物质燃烧源各二级源类的排放量贡献率见图 1, 户用生物质炉具CO、SO2、NH3和BC的排放量贡献率最高, 分别为53.27%、52.00%、77.78%和59.69%; 秸秆露天焚烧NOx、VOCs、PM2.5、PM10和OC的排放量贡献率最大, 分别为53.75%、61.02%、56.22%、55.66%和77.16%.这与其他城市的研究结果一致, 如卢轩等[68]的研究发现, 郑州市秸秆家用和秸秆露天焚烧分别占生物质燃烧源VOCs的43.79%和41.24%; 张立斌等[69]的研究发现监利县生物质开放燃烧VOCs、NH3、PM10、PM2.5和OC的排放量贡献率最高, 生物质炉灶CO和BC的排放量贡献率最高; 高玉宗等[70]的研究建立了2018年西宁市生物质燃烧源大气污染物排放清单, 结果表明秸秆露天焚烧CO、NOx、VOCs、PM2.5、PM10、BC和OC的排放量贡献率最高, 户用生物质NH3和SO2的排放量贡献率最大; 刘春蕾等[71]在2015年南京市生物质燃烧源大气污染物排放清单中指出户用生物质炉具PM10、PM2.5、BC、OC、CO、VOCs和NH3的排放量贡献率最高.

图 1 生物质燃烧源大气污染物各二级源排放量贡献率 Fig. 1 Contribution rate of emissions from secondary sources of air pollutants from biomass combustion sources

2.2 空间分布特征

表 5和图 2分别为2019年贵州省生物质燃烧源CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC这9项污染物不同市(州)的排放量和贡献率.对于生物质锅炉, 根据排污设施经纬度将大气污染物排放量直接定位到相应网格; 对于户用生物质炉具, 由农村人口矢量数据计算出各市(州)农村人口数, 进而计算出各市(州)单位农村人口数对应的大气污染物排放量, 最后与网格化后3 km×3 km分辨率的行政区划数据连接计算, 得到相应网格的户用生物质炉具排放; 同理, 由森林、草原受灾矢量数据和旱地矢量数据计算, 分别得到相应网格的森林、草原火灾排放和秸秆露天焚烧排放; 将各网格中各二级源大气污染物排放量叠加, 利用地理信息系统ArcGIS工具[72], 得到贵州省生物质燃烧源大气污染物排放空间分布如图 3所示.

表 5 生物质燃烧源大气污染物不同市(州)排放量/t Table 5 Emission inventory of biomass combustion in different cities or prefectures/t

图 2 生物质燃烧源大气污染物不同市(州)排放量贡献率 Fig. 2 Emission contributions of air pollutants from biomass combustion sources in different cities

图 3 2019年贵州省生物质燃烧源大气污染物空间分布 Fig. 3 Space distribution characteristics of biomass combustion in Guizhou Province in 2019

表 5、图 2和图 3表明, 2019年贵州省生物质燃烧源排放区域主要集中在黔东南苗族侗族自治州, CO、NOx、SO2、NH3、VOCs、PM2.5、PM10和BC的排放量贡献率最大.这是因为, 黔东南苗族侗族自治州的户用生物质炉具的薪柴消费量(143.97×104 t), 远大于其他市(州)的薪柴消费量(13.48×104 t~91.71×104 t).而对于OC的排放量贡献率是遵义市最大, 原因是其户用生物质炉具的高粱秸秆消费量、秸秆露天焚烧的高粱、油菜和辣椒秸秆消费量比其他市(州)大, 分别为0.12×104、24.35×104、27.06×104和152.47×104 t.各污染物排放的空间贡献率市(州)排序略有区别, 具体如下.①CO：黔东南苗族侗族自治州(19.02%)＞遵义市(15.39%)＞黔南布依族苗族自治州(14.84%)＞毕节市(13.60%)＞铜仁市(13.40%)＞黔西南布依族苗族自治州(7.47%)＞贵阳市(5.54%)＞安顺市(5.54%)＞六盘水市(5.21%); ②NOx：黔东南苗族侗族自治州(15.22%)＞遵义市(14.83%)＞毕节市(14.61%)＞黔南布依族苗族自治州(13.18%)＞贵阳市(12.94%)＞铜仁市(10.45%)＞黔西南布依族苗族自治州(7.55%)＞安顺市(6.08%)＞六盘水市(5.14%); ③SO2：黔东南苗族侗族自治州(19.19%)＞黔南布依族苗族自治州(14.53%)＞遵义市(13.40%)＞铜仁市(12.88%)＞毕节市(11.80%)＞贵阳市(11.42%)＞黔西南布依族苗族自治州(6.48%)＞安顺市(5.77%)＞六盘水市(4.53%); ④NH3：黔东南苗族侗族自治州(24.54%)＞铜仁市(16.02%)＞黔南布依族苗族自治州(15.17%)＞遵义市(12.59%)＞毕节市(11.85%)＞黔西南布依族苗族自治州(5.53%)＞安顺市(5.01%)＞贵阳市(4.66%)＞六盘水市(4.64%); ⑤VOCs：黔东南苗族侗族自治州(17.37%)＞遵义市(16.78%)＞黔南布依族苗族自治州(14.55%)＞毕节市(13.76%)＞铜仁市(12.31%)＞黔西南布依族苗族自治州(7.97%)＞安顺市(6.11%)＞贵阳市(6.04%)＞六盘水市(5.10%); ⑥PM2.5：黔东南苗族侗族自治州(17.79%)＞遵义市(16.18%)＞毕节市(14.75%)＞黔南布依族苗族自治州(14.52%)＞铜仁市(12.60%)＞黔西南布依族苗族自治州(8.05%)＞安顺市(5.75%)＞六盘水市(5.49%)＞贵阳市(4.87%); ⑦PM10：黔东南苗族侗族自治州(17.84%)＞遵义市(16.17%)＞毕节市(14.81%)＞黔南布依族苗族自治州(14.50%)＞铜仁市(12.68%)＞黔西南布依族苗族自治州(7.95%)＞安顺市(5.78%)＞六盘水市(5.45%)＞贵阳市(4.82%); ⑧BC：黔东南苗族侗族自治州(20.55%)＞黔南布依族苗族自治州(15.06%)＞遵义市(14.51%)＞铜仁市(14.24%)＞毕节市(13.72%)＞黔西南布依族苗族自治州(7.05%)＞安顺市(5.40%)＞六盘水市(5.22%)＞贵阳市(4.26%); ⑨OC：遵义市(17.91%)＞毕节市(16.35%)＞黔南布依族苗族自治州(14.24%)＞黔东南苗族侗族自治州(13.71%)＞铜仁市(10.51%)＞黔西南布依族苗族自治州(9.46%)＞安顺市(6.27%)＞六盘水市(5.92%)＞贵阳市(5.62%).

2.3 时间分布特征

贵州省生物质燃烧源大气污染物排放时间分配方法见表 6, 其分配结果具体为：①从生物质燃烧源排放月分配系数可见(图 4), 生物质锅炉月排放量较平均, 可以看出企业基本处于常年稳定生产状态; 户用生物质炉具(秸秆和薪柴)大气污染物排放具有明显的季节性变化特征, 采暖季12月至次年2月的各污染物排放量高于非采暖季4~10月的排放, 主要是因为采暖季供暖加大了对户用生物质燃烧的需求量, 从而增加了污染物的排放量, 3月和11月属于采暖季和非采暖季的交替季; 森林和草原火灾多发生在3月, 其次是2、4、8和12月, 其他月份的分配系数之和仅为0.2; 秸秆露天焚烧月排放量变化明显, 多集中在2~4月, 其分配系数合计为0.6, 其次是8月和12月, 其他月份的分配系数之和仅为0.2, 这与农作生产有关, 具有较强的季节性[67]. ②从生物质燃烧源小时分配系数可见(图 5), 生物质锅炉在每日06:00~07:00和17:00~20:00的排放量逐渐达到最大; 户用生物质炉具(秸秆和薪柴)在07:00~09:00和18:00~20:00排放量较大, 与当地居民炊事习惯一致; 森林和草原火灾主要集中在14:00~15:00, 对应的分配系数合计为0.5; 秸秆露天焚烧多集中在11:00~12:00和14:00~15:00, 其分配系数合计为0.7, 这与农民劳作时段有关.

表 6 生物质燃烧源大气污染物排放时间分配方法 Table 6 Time distribution method of atmospheric pollutant emissions from biomass combustion sources

图 4 生物质燃烧源排放月分配系数 Fig. 4 Monthly emission coefficient of biomass combustion

图 5 生物质燃烧源排放小时分配系数 Fig. 5 Hourly emission coefficient of biomass combustion

2.4 不确定性分析

在排放清单的定量估算中, 不确定性来源于活动水平和排放系数的选取.不确定性分析可用于重要污染源信息的甄别, 评估排放清单的可靠性[69].本研究采用蒙特卡洛方法评估排放总量的置信区间, 量化分析2019年贵州省生物质燃烧源排放清单的不确定性.使用Excel软件中的Oracle Crystal Ball插件进行蒙特卡洛随机模拟抽样, 在输入量定义假设中, 定义户用生物质炉具(秸秆和薪柴)的活动水平为正态分布, 其他燃烧源的活动水平为均匀分布; 主要来源是本地实测的户用生物质炉具薪柴、辣椒、大豆和玉米秸秆的BC和OC排放系数为正态分布, 来源是指南或文献的其他燃烧源各污染物排放系数为均匀分布.污染物排放量设为输出量, 建立输出量和输入量之间的变量关系, 通过10 000次模拟计算, 获得输出量概率分布在95%置信区间的不确定性范围(表 7).

表 7 排放清单的不确定性 Table 7 Uncertainty of emission inventory

由表 7可知, 2019年贵州省生物质燃烧源大气污染物排放清单的不确定性总体范围为-60.18%~144.67%, 主要由活动水平和排放系数的不确定性共同导致, 其中, 户用生物质炉具薪柴的VOCs排放的不确定性范围最大, 为-48.32%~144.24%; 户用生物质炉具薪柴的NOx排放的不确定性范围最小, 为-34.85%~79.51%.对于生物锅炉排放清单, 由于活动水平来自环境统计调查[47], 可信度较高, 活动水平对不确定性的敏感度比排放系数低, 其中涉及去除效率的污染物(如PM2.5、PM10、BC和OC)活动水平敏感度比其他污染物(如CO、NOx、SO2、NH3和VOCs)高, 这是因为活动水平的变量(即去除效率)增加, 导致不确定性范围变大.对于户用生物质炉具排放清单, 薪柴、辣椒、大豆和玉米秸秆的BC和OC排放系数来源为本地实测数据, 有效降低了其不确定性, 对应的排放系数敏感度也显著降低.对于森林和草原火灾, 活动水平主要来自政府部门提供的数据[58, 59], 活动水平对不确定性的敏感度比排放系数低, 而对于秸秆露天焚烧, 由于排放系数细化到不同类型的秸秆, 降低了排放系数对不确定性的敏感度, 且低于活动水平的敏感度.

3 结论

(1) 2019年贵州省生物质燃烧源CO、NOx、SO2、NH3、VOCs、PM2.5、PM10、BC和OC的排放量分别为: 293 505.53、14 781.19、4 146.11、8 501.07、45 025.70、39 463.58、41 879.31、6 832.33和15 134.74 t.从排放量贡献率来看, 户用生物质炉具的CO、SO2、NH3和BC最大, 分别为: 53.27%、52.00%、77.78%和59.59%; 秸秆露天焚烧的NOx、VOCs、PM2.5、PM10和OC最大, 分别为: 53.75%、61.02%、56.22%、55.66%和77.16%.

(2) 各市(州)排放的污染物分布不均衡, 主要集中于黔东南苗族侗族自治州, CO、NOx、SO2、NH3、VOCs、PM2.5、PM10和BC的排放量贡献率最大, 分别为: 19.02%、15.22%、19.19%、24.54%、17.37%、17.79%、17.84%和20.55%.OC的排放量贡献率则是遵义市最大(17.91%).

(3) 各二级源污染物排放量变化特征不同.生物质锅炉月排放量较平均; 户用生物质炉具呈现出明显的季节性变化, 采暖季污染物排放量远高于非采暖季; 森林和草原火灾多发生在3月, 其次是2、4、8和12月; 秸秆露天焚烧集中在2~4月.生物质锅炉在06:00~07:00和17:00~20:00的排放量较大; 户用生物质炉具在07:00~09:00和18:00~20:00的排放量较大; 森林和草原火灾集中在14:00~15:00; 秸秆露天焚烧集中在11:00~12:00和14:00~15:00.

(4) 本研究建立的排放清单能够较好地代表贵州省生物质燃烧源大气污染物的排放现状, 但仍然存在一定不确定性(-60.18%~144.67%), 主要由活动水平和排放系数的不确定性共同导致.后续研究中应进一步提升活动水平数据获取的准确性, 继续开展排放更细分的系数本地化实测, 从而降低不确定性, 完善贵州省生物质燃烧源排放清单, 为贵州省大气污染协同治理提供科学依据.

参考文献
[1] Demirbas A. Potential applications of renewable energy sources, biomass combustion problems in boiler power systems and combustion related environmental issues[J]. Progress in Energy and Combustion Science, 2005, 31(2): 171-192. DOI:10.1016/j.pecs.2005.02.002
[2] 田贺忠, 郝吉明, 陆永琪, 等. 中国生物质燃烧排放SO2, NOx量的估算[J]. 环境科学学报, 2002, 22(2): 204-208.
Tian H Z, Hao J M, Lu Y Q, et al. Evaluation of SO2 and NOx emissions resulted from biomass fuels utilization in China[J]. Acta Scientiae Circumstantiae, 2002, 22(2): 204-208.
[3] 田贺忠, 赵丹, 王艳. 中国生物质燃烧大气污染物排放清单[J]. 环境科学学报, 2011, 31(2): 349-357.
Tian H Z, Zhao D, Wang Y. Emission inventories of atmospheric pollutants discharged from biomass burning in China[J]. Acta Scientiae Circumstantiae, 2011, 31(2): 349-357. DOI:10.13671/j.hjkxxb.2011.02.018
[4] 陆炳, 孔少飞, 韩斌, 等. 2007年中国大陆地区生物质燃烧排放污染物清单[J]. 中国环境科学, 2011, 31(2): 186-194.
Lu B, Kong S F, Han B, et al. Inventory of atmospheric pollutants discharged from biomass burning in China continent in 2007[J]. China Environmental Science, 2011, 31(2): 186-194.
[5] Zhou Y, Xing X F, Lang J L, et al. A comprehensive biomass burning emission inventory with high spatial and temporal resolution in China[J]. Atmospheric Chemistry and Physics, 2017, 17(4): 2839-2864. DOI:10.5194/acp-17-2839-2017
[6] Takami K, Shimadera H, Uranishi K, et al. Impacts of biomass burning emission inventories and atmospheric reanalyses on simulated PM10 over Indochina[J]. Atmosphere, 2020, 11(2). DOI:10.3390/atmos11020160
[7] Puliafito S E, Bolaño-Ortiz T, Berná L, et al. High resolution inventory of atmospheric emissions from livestock production, agriculture, and biomass burning sectors of Argentina[J]. Atmospheric Environment, 2020, 223. DOI:10.1016/j.atmosenv.2019.117248
[8] Mataveli G A V, de Oliveira G, Seixas H T, et al. Relationship between biomass burning emissions and deforestation in Amazonia over the last two decades[J]. Forests, 2021, 12(9). DOI:10.3390/F12091217
[9] Sahu S K, Mangaraj P, Beig G, et al. Quantifying the high resolution seasonal emission of air pollutants from crop residue burning in India[J]. Environmental Pollution, 2021, 286. DOI:10.1016/j.envpol.2021.117165
[10] Le H A, Phuong D M, Linh L T. Emission inventories of rice straw open burning in the Red River Delta of Vietnam: evaluation of the potential of satellite data[J]. Environmental Pollution, 2020, 260. DOI:10.1016/j.envpol.2020.113972
[11] Le H A, Khoi N Q, Mallick J. Integrated emission inventory and modelling to assess the distribution of particulate matters from rice straw open burning in Hanoi, Vietnam[J]. Atmospheric Pollution Research, 2022, 13(5). DOI:10.1016/j.apr.2022.101416
[12] Das B, Bhave P V, Puppala S P, et al. A model-ready emission inventory for crop residue open burning in the context of Nepal[J]. Environmental Pollution, 2020, 266. DOI:10.1016/j.envpol.2020.115069
[13] Pino-Cortés E, Carrasco S, Díaz-Robles L A, et al. Emission inventory processing of biomass burning from a global dataset for air quality modeling[J]. Air Quality, Atmosphere & Health, 2022, 15(4): 721-729.
[14] 杨光义. 东北地区生物质露天燃烧源排放及其对大气环境的影响与评估[D]. 长春: 中国科学院大学(中国科学院东北地理与农业生态研究所), 2020.
Yang G Y. Emissions of open biomass burning and its impacts on atmospheric environment in the Northeastern China[D]. Changchun: Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, 2020.
[15] 徐敬, 寇星霞, 李梓铭. 秸秆燃烧排放对北京及其周边地区PM2.5浓度影响的数值模拟[J]. 气候与环境研究, 2018, 23(5): 587-595.
Xu J, Kou X X, Li Z M. Modeling the impact of emissions from crop residue burning on PM2.5 concentration in Beijing and its surrounding areas[J]. Climatic and Environmental Research, 2018, 23(5): 587-985.
[16] 倪紫琳, 周亚端, 张银菊, 等. 鄂州市高时空分辨率大气污染源排放清单的建立[J]. 环境科学与技术, 2021, 44(2): 90-96.
Ni Z L, Zhou Y D, Zhang Y J, et al. Establishment of high resolution atmospheric emission inventory in Ezhou city[J]. Environmental Science & Technology, 2021, 44(2): 90-96.
[17] 马强, 宫密秘, 杨晓霞, 等. 嘉兴市人为源大气污染物排放清单研究[J]. 环境污染与防治, 2020, 42(6): 660-665.
Ma Q, Gong M M, Yang X X, et al. Air pollutants emission inventory of anthropogenic sources in Jiaxing[J]. Environmental Pollution & Control, 2020, 42(6): 660-665.
[18] 庞可, 张芊, 马彩云, 等. 天水市人为源挥发性有机物排放清单及特征研究[J]. 环境科学与管理, 2021, 46(10): 101-106.
Pang K, Zhang Q, Ma C Y, et al. Emission inventory and characteristic analysis of anthropogenic volatile organic compounds in Tianshui City[J]. Environmental Science and Management, 2021, 46(10): 101-106.
[19] 范西彩, 张新民, 张晓红, 等. 鹤壁市大气挥发性有机物源排放清单研究[J]. 中国环境科学, 2021, 41(2): 558-565.
Fan X C, Zhang X M, Zhang X H, et al. Research on the emission inventory of volatile organic compounds in Hebi City, Henan Province[J]. China Environmental Science, 2021, 41(2): 558-565.
[20] 高素莲, 闫学军, 张荣荣, 等. 济南市人为源挥发性有机物排放清单及特征[J]. 环境科学与技术, 2021, 44(S1): 191-198.
Gao S L, Yan X J, Zhang R R, et al. Emission inventory and characteristics of anthropogenic volatile organic compounds in Jinan[J]. Environmental Science & Technology, 2021, 44(S1): 191-198.
[21] 郭文凯, 李光耀, 陈冰, 等. 兰州市高分辨率人为源排放清单建立及在WRF-Chem中应用评估[J]. 环境科学, 2021, 42(2): 634-642.
Guo W K, Li G Y, Chen B, et al. Establishment of a high-resolution anthropogenic emission inventory and its evaluation using the WRF-Chem Model for Lanzhou[J]. Environmental Science, 2021, 42(2): 634-642.
[22] 陈东, 陈军辉, 钱骏, 等. 资阳市大气污染源排放清单研究[J]. 四川环境, 2021, 40(1): 68-75.
Chen D, Chen J H, Qian J, et al. Research on air pollutant emission inventory in Zi Yang City[J]. Sichuan Environment, 2021, 40(1): 68-75.
[23] 江沛玉, 陈晓玲, 李玲玉. 山东省高分辨率大气颗粒物排放清单[J]. 环境污染与防治, 2021, 43(9): 1124-1132.
Jiang P Y, Chen X L, Li Y L. High-resolution emission inventory of atmospheric particulate matter in Shandong[J]. Environmental Pollution & Control, 2021, 43(9): 1124-1132.
[24] 周子航, 邓也, 谭钦文, 等. 四川省人为源大气污染物排放清单及特征[J]. 环境科学, 2018, 39(12): 5344-5358.
Zhou Z H, Deng Y, Tan Q W, et al. Emission inventory and characteristics of anthropogenic air pollutant sources in the Sichuan Province[J]. Environmental Science, 2018, 39(12): 5344-5358.
[25] Shi Y S, Gong S Y, Zang S Y, et al. High-resolution and multi-year estimation of emissions from open biomass burning in Northeast China during 2001-2017[J]. Journal of Cleaner Production, 2021, 310. DOI:10.1016/j.jclepro.2021.127496
[26] Zhao H M, Yang G Y, Tong D Q, et al. Interannual and seasonal variability of greenhouse gases and aerosol emissions from biomass burning in Northeastern China constrained by satellite observations[J]. Remote Sensing, 2021, 13(5). DOI:10.3390/rs13051005
[27] Wu J, Kong S F, Zeng X, et al. First high-resolution emission inventory of levoglucosan for biomass burning and non-biomass burning sources in China[J]. Environmental Science & Technology, 2021, 55(3): 1497-1507.
[28] Wu J, Kong S F, Yan Y Y, et al. Neglected biomass burning emissions of air pollutants in China-views from the corncob burning test, emission estimation, and simulations[J]. Atmospheric Environment, 2022, 278. DOI:10.1016/j.atmosenv.2022.119082
[29] Cui L L, Li R, Fu H B, et al. Nitrous acid emission from open burning of major crop residues in mainland China[J]. Atmospheric Environment, 2021, 244. DOI:10.1016/j.atmosenv.2020.117950
[30] 曹国良, 张小曳, 王丹, 等. 中国大陆生物质燃烧排放的污染物清单[J]. 中国环境科学, 2005, 25(4): 389-393.
Cao G L, Zhang X Y, Wang D, et al. Inventory of atmospheric pollutants discharged from biomass burning in China continent[J]. China Environmental Science, 2005, 25(4): 389-393.
[31] Zhang X H, Lu Y, Wang Q G, et al. A high-resolution inventory of air pollutant emissions from crop residue burning in China[J]. Atmospheric Environment, 2019, 213: 207-214.
[32] Yin L F, Du P, Zhang M S, et al. Estimation of emissions from biomass burning in China (2003-2017) based on MODIS fire radiative energy data[J]. Biogeosciences, 2019, 16(7): 1629-1640.
[33] 张晓荟. 中国秸秆焚烧大气污染物高分辨率排放特征研究[D]. 南京: 南京大学, 2019.
Zhang X H. High-resolution characteristics of air pollutant emissions from crop residue burning in China[D]. Nanjing: Nanjing University, 2019.
[34] Li R M, Chen W W, Zhao H M, et al. Inventory of atmospheric pollutant emissions from burning of crop residues in China based on satellite-retrieved farmland data[J]. Chinese Geographical Science, 2020, 30(2): 266-278.
[35] 贵州省生态环境厅. 贵州省生态环境质量月报(2022年3月)[EB/OL]. https://sthj.guizhou.gov.cn/hjsj/hjzlsjzx_5802731/hjzlyb_5802733/202204/P020220428544739326382.pdf, 2022-04-28.
[36] 蔡景行, 陆海, 贾田, 等. 贵阳市主要农作物秸秆资源空间分布与利用[J]. 农技服务, 2020, 37(12): 104-106.
[37] 王璐, 严如江. 贵州省兴仁市农作物秸秆利用情况调查分析[J]. 贵州农业科学, 2019, 47(12): 134-137.
Wang L, Yan R J. Survey and analysis of crop straw utilization in Xingren City, Guizhou[J]. Guizhou Agricultural Science, 2019, 47(12): 134-137.
[38] 郭应军, 熊康宁, 安裕伦, 等. 中国西南石漠化地区农村能源消费结构研究[J]. 农业工程学报, 2019, 35(3): 226-234.
Guo Y J, Xiong K N, An Y L, et al. Study on rural household energy consumption structure in rocky desertification region of Southwest China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(3): 226-234.
[39] 王琦, 熊康宁, 盈斌, 等. 喀斯特地区农村能源消费特征及影响因素[J]. 中国农业资源与区划, 2019, 40(3): 146-157.
Wang Q, Xiong K N, Ying B, et al. Consumption characteristics of rural energies in Karst Area and their influencing factors[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2019, 40(3): 146-157.
[40] 张雅蓉, 李渝, 蒋太明, 等. 贵州主要农作物秸秆资源分布及综合利用现状[J]. 贵州农业科学, 2015, 43(8): 262-267.
Zhang Y R, Li Y, Jiang T M, et al. Distribution and comprehensive utilization of straw resources of main crops in Guizhou[J]. Guizhou Agricultural Sciences, 2015, 43(8): 262-267.
[41] 陈志祯, 黄瑞, 黄光谋, 等. 罗甸县农作物秸秆资源调查与应用展望[J]. 当代畜牧, 2015(33): 23-24.
[42] 唐亮, 刘肇军, 熊康宁, 等. 喀斯特地区农村家庭能源消费结构及效应分析——以贵州毕节撒拉溪示范区为例[J]. 贵州师范大学学报(自然科学版), 2013, 31(2): 96-102.
Tang L, Liu Z J, Xiong K N, et al. Study on structure and effects of rural household energy in rocky desertification area-A case in salaxi of Guizhou[J]. Journal of Guizhou Normal University (Natural Sciences), 2013, 31(2): 96-102.
[43] 贵州省环境科学研究设计院. 贵州省大气污染源排放清单成果汇总报告[R]. 贵阳: 贵州省环境科学研究设计院, 2017. 273-290.
[44] 贺克斌. 城市大气污染源排放清单编制技术手册[EB/OL]. https://max.book118.com/html/2018/1116/8107025055001133.shtm, 2018-11-16.
[45] 环境保护部. 关于发布《大气可吸入颗粒物一次源排放清单编制技术指南(试行)》等5项技术指南的公告[EB/OL]. https://www.mee.gov.cn/gkml/hbb/bgg/201501/t20150107_293955.htm, 2014-12-31.
[46] 贵州省环境监测中心站调查研究项目组. 贵州省典型地区村镇散煤燃烧调查与空气质量影响研究[R]. 贵阳: 贵州省环境监测中心站, 2019. 13-14.
[47] 贵州省生态环境厅. 2019贵州环境统计[R]. 贵阳: 贵州省生态环境厅, 2019.
[48] 贵阳市统计局, 国家统计局贵阳调查队. 贵阳统计年鉴2020[M]. 北京: 中国统计出版社, 2020.
Bureau of Statistics of Guiyang, NBS Survey Office in Guiyang. Guiyang statistical yearbook 2020[M]. Beijing: China Statistics Press, 2020.
[49] 六盘水市统计局. 六盘水2020年统计年鉴[EB/OL]. http://tjj.gzlps.gov.cn/bmxxgk/zfxxgk/fdzdgknr/zdly/tjnj_5712553/202106/t20210615_68512543.html, 2021-06-15.
[50] 遵义市统计局, 国家统计局遵义调查队. 2020年遵义统计年鉴[EB/OL]. http://www.zunyi.gov.cn/sjfb/tjxx/tjnj_5860125/202011/t20201106_68833238.html, 2020-11-06.
[51] 安顺市统计局, 国家统计局安顺调查队. 2020年安顺统计年鉴[EB/OL]. http://www.anshun.gov.cn/zfsj/tjnj/202112/t20211207_71948788.html, 2021-12-07.
[52] 毕节市统计局, 国家统计局毕节调查队. 毕节统计年鉴2020[M]. 北京: 中国统计出版社, 2021.
[53] 铜仁市统计局. 2020年铜仁统计年鉴[EB/OL]. http://www.trs.gov.cn/zfsj/tjnj/202106/t20210602_68376344.html, 2021-06-02.
[54] 黔西南州统计局, 国家统计局黔西南调查队. 2019年黔西南州统计年鉴[EB/OL]. http://www.qxn.gov.cn/zwgk/zfjg/ztjj_5135341/bmxxgkml_5135344/tjgb_5135353/202108/t20210825_69768473.html, 2021-08-25.
[55] 黔东南州统计局. 2019年黔东南统计年鉴[EB/OL]. http://tjj.qdn.gov.cn/tjsj/sjcx_57102/qdntjnj/202203/t20220316_73006053.html, 2022-03-16.
[56] 黔南州统计局. 黔南统计年鉴2020[M]. 北京: 中国统计出版社, 2020.
Qiannan Autonomous Prefectural Bureau of Statistics. Qiannan statistical yearbook 2020[M]. Beijing: China Statistics Press, 2020.
[57] 贵州省统计局, 国家统计局贵州调查总队. 贵州统计年鉴2020[M]. 北京: 中国统计出版社, 2020.
Guizhou Province Bureau of Statistics, NBS Survey Office in Guizhou. Guizhou statistical yearbook 2020[M]. Beijing: China Statistics Press, 2020.
[58] 贵州省林业局. 2019年贵州林业年鉴[EB/OL]. https://www.guizhou.gov.cn/dcgz/stgz/lynj/202109/t20210915_70433528.html, 2021-09-15.
[59] 贵州省应急厅. 贵州省森林火灾过火面积统计表[R]. 2022.
[60] GB/T 17297-1998, 中国气候区划名称与代码气候带和气候大区[S].
[61] 国家林业和草原局. 中国林业和草原统计年鉴2020[M]. 北京: 中国林业出版社, 2021.
[62] 贵州省第三次全国国土调查领导小组办公室, 贵州省自然资源厅, 贵州省统计局. 贵州省第三次全国国土调查主要数据公报[EB/OL]. https://zrzy.guizhou.gov.cn/wzgb/zwgk/zdlyxxgk/zrzy/dcjc/202201/t20220121_72378280.html, 2021-12-28.
[63] NASA. Fire information for resource management system[EB/OL]. https://firms.modaps.eosdis.nasa.gov/download/, 2019-12-31.
[64] 唐喜斌, 黄成, 楼晟荣, 等. 长三角地区秸秆燃烧排放因子与颗粒物成分谱研究[J]. 环境科学, 2014, 35(5): 1623-1632.
Tang X B, Huang C, Lou S R, et al. Emission factors and PM chemical composition study of biomass burning in the Yangtze River delta region[J]. Environmental Science, 2014, 35(5): 1623-1632.
[65] 何敏, 王幸锐, 韩丽, 等. 四川省秸秆露天焚烧污染物排放清单及时空分布特征[J]. 环境科学, 2015, 36(4): 1208-1216.
He M, Wang X R, Han L, et al. Emission inventory of crop residues field burning and its temporal and spatial distribution in Sichuan Province[J]. Environment Science, 2015, 36(4): 1208-1216.
[66] 王文杰, 田杰, 张勇, 等. 我国典型秸秆露天燃烧黑碳的高时间分辨率排放特征及其排放清单[J]. 地球环境学报, 2019, 10(2): 190-200.
Wang W J, Tian J, Zhang Y, et al. High time resolution emission characteristics and emission inventory of black carbon in typical crop residues burning in China[J]. Journal of Earth Environment, 2019, 10(2): 190-200.
[67] 杨夏捷, 马远帆, 鞠园华, 等. 华南农产品主产区2005—2014年秸秆露天燃烧污染物排放估算及时空分布[J]. 农业环境科学学报, 2018, 37(2): 358-368.
Yang X J, Ma Y F, Ju Y H, et al. Temporal and spatial distribution of air pollutants emitted from field burning of straw crops in Southern China during 2005-2014[J]. Journal of Agro-Envirnomental Science, 2018, 37(2): 358-368.
[68] 卢轩, 张瑞芹, 韩跞锎. 郑州市VOCs组分排放清单及其臭氧生成潜势[J]. 环境科学, 2020, 41(10): 4426-4435.
Lu X, Zhang R Q, Han L K. Emission inventory of VOCs components in Zhengzhou and their ozone formation potential[J]. Environmental Science, 2020, 41(10): 4426-4435.
[69] 张立斌, 黄凡, 张银菊, 等. 监利县大气污染源排放清单及特征研究[J]. 环境科学与技术, 2020, 43(9): 182-189.
Zhang L B, Huang F, Zhang Y J, et al. Emission inventory and characteristics of atmospheric pollution sources in Jianli County[J]. Environmental Science & Technology, 2020, 43(9): 182-189.
[70] 高玉宗, 姬亚芹, 林孜, 等. 西宁市生物质燃烧源大气污染物排放清单[J]. 环境科学, 2021, 42(12): 5585-5593.
Gao Y Z, Ji Y Q, Lin Z, et al. Pollutant emission inventory of biomass combustion sources in Xining City[J]. Environmental Science, 2021, 42(12): 5585-5593.
[71] 刘春蕾, 杨峰. 南京市生物质燃烧源大气污染物排放清单及其特征[J]. 安徽农学通报, 2017, 23(16): 95-97.
Liu C L, Yang F. Emission inventories of atmospheric pollutants discharged from biomass burning in Nanjing City[J]. Anhui Agricultural Science Bulletin, 2017, 23(16): 95-97.
[72] 王艳妮, 秦俊虎, 杨敬婷, 等. 基于入户调查的贵阳市生活燃煤排放清单[J]. 环境科学学报, 2021, 41(5): 1809-1817.
Wang Y N, Qin J H, Yang J T, et al. Household survey based study on emission inventory of residential coal burning in city of Guiyang[J]. Acta Scientiae Circumstantiae, 2021, 41(5): 1809-1817.