留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析

丰睿 黄成臣 高寒 郑慧君 申亚梅 罗坤

丰睿, 黄成臣, 高寒, 郑慧君, 申亚梅, 罗坤. 后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析[J]. 浙江农林大学学报, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
引用本文: 丰睿, 黄成臣, 高寒, 郑慧君, 申亚梅, 罗坤. 后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析[J]. 浙江农林大学学报, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
FENG Rui, HUANG Chengchen, GAO Han, ZHENG Huijun, SHEN Yamei, LUO Kun. Decoding VOCs and PM10 for Hangzhou in the post-G20 period[J]. Journal of Zhejiang A&F University, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
Citation: FENG Rui, HUANG Chengchen, GAO Han, ZHENG Huijun, SHEN Yamei, LUO Kun. Decoding VOCs and PM10 for Hangzhou in the post-G20 period[J]. Journal of Zhejiang A&F University, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022

后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析

doi: 10.11833/j.issn.2095-0756.2019.04.022
基金项目: 

国家自然科学基金资助项目 51476144

详细信息
    作者简介: 丰睿, 工程师, 博士, 从事能源和环境工程研究。E-mail:11727081@zju.edu.cn
    通信作者: 罗坤, 教授, 博士, 从事能源及环境工程领域复杂多相流和反应的理论及数值模拟研究。E-mail:zjulk@zju.edu.cn
  • 中图分类号: X511

Decoding VOCs and PM10 for Hangzhou in the post-G20 period

  • 摘要: 后G20时期(2016年9月至2017年12月),臭氧成为杭州夏季首要大气污染物,可吸入颗粒物(PM10)成为杭州冬季首要大气污染物。利用杭州市环境检测朝晖站和环境监测下沙站空气质量数据观测值分析55种挥发性有机化合物(VOC)的臭氧形成潜力,并使用美国国家环境保护局研发的第3代大气质量预估模型(WRF-CMAQ)对杭州市中心PM10的污染数值进行模拟分析,提出了后G20时期PM10和臭氧的最佳管控方案。结果显示:在后G20时期杭州市中心和城郊,芳香烃所占臭氧形成潜力比例最高,其次为烯烃和烷烃,最后为乙炔。在市区和城郊,乙烯、间/对二甲苯、甲苯、丙烯和乙苯依次为生成臭氧最主要的5种挥发性有机化合物,因此管控这5种挥发性有机化合物能最大程度的减少臭氧质量浓度。2017年杭州市中心春夏秋冬四季的PM10中分别有50%,32%,48%和45%源于杭州市以外污染源的长距离跨地域传输。其中,春秋冬三季通过长距离传输来到杭州市区的PM10主要源于市区的北方,夏季则主要源于市区的西南方。2017年,杭州市区的交通源、工业源、生活源与农业源对本地排放PM10的贡献率分别为62.6%,27.8%,7.3%和2.3%。当杭州本地的工业源和交通源在2016年的基础上分别减排15%和5%时,杭州市中心的PM10年均质量浓度可低于70 μg·m-3,达到二类区的标准。
  • 图  1  杭州市朝晖站和下沙站月均TVOC体积分数

    Figure  1  Monthly average TVOC concentration at Zhaohui and Xiasha stations

    图  2  2017年1,4,7,10月长江三角洲PM10质量浓度的时空分布示意图

    Figure  2  Time-spacing map of PM10 over Yangtze River Delta in 2017

    图  3  2017年1,4,7,10月关停杭州市本地污染源后的PM10时空分布示意图

    Figure  3  Time-spacing map of PM10 in 2017 when Hangzhou local sources are shut down

    图  4  杭州市2017年风玫瑰图

    Figure  4  Wind-rose over Hangzhou in 2017

    表  1  WRF参数方案

    Table  1.   Parameters for WRF

    参数 选择的方案
    微观物理进程 WSM6[1]
    长波辐射 RRTM[2]
    短波辐射 Goddard[3]
    表层 MM5 similarity[4]
    陆地表面进程 Noah[5]
    边界层 YSU[6]
    积雨云参数化进程 Grell-Devenyi[7]
    下载: 导出CSV

    表  2  2017年度各季节不同种类VOC排放比例

    Table  2.   Seasonal contribution of different subgroup VOC in 2017

    季节 烷烃排放比例/% 烯烃排放比例/% 芳香烃排放比例/% 炔烃排放比例/%
    朝晖站 下沙站 朝晖站 下沙站 朝晖站 下沙站 朝晖站 下沙站
    59.7 70.7 15.4 5.5 24.0 22.6 0.9 1.2
    73.1 69.3 9.6 8.0 16.6 21.3 0.7 1.4
    62.6 58.0 7.9 10.0 28.9 30.9 0.6 1.1
    60.4 65.2 10.2 8.8 28.9 24.4 0.51.6
    下载: 导出CSV

    表  3  杭州本地污染源对市中心PM10的贡献率

    Table  3.   Hangzhou local sector contribution for PM10 at city center

    月份贡献率/%
    农业源 工业源 生活源 交通源
    1 2.6 61.1 8.4 27.9
    4 2.1 62.6 6.7 28.6
    7 3.0 60.2 7.2 29.6
    10 1.5 66.4 6.8 25.3
    平均 2.3 62.6 7.327.8
    下载: 导出CSV
  • [1] XIAO Huadong, SUN Jing, BIAN Xiaofeng, et al. GPU acceleration of the WSM6 cloud microphysics scheme in GRAPES model[J]. Comput Geosci, 2013, 59:156-162.
    [2] FOUNTOUKIS C, MARTÍN-POMARES L, PEREZ-ASTUDILLO D, et al. Simulating global horizontal irradiance in the Arabian Peninsula:sensitivity to explicit treatment of aerosols[J]. Sol Energ, 2018, 163:347-355.
    [3] MOHR C, TAO Weikuo, CHERN J, et al. The NASA goddard multi-scale modeling framework-land information system:global land/atmosphere interaction with resolved convection[J]. Environ Modell Software, 2013, 39:103-115.
    [4] DENG Tao, CHEN Yang, WAN Qilin, et al. Comparative evaluation of the impact of GRAPES and MM5 meteorology on CMAQ prediction over Pearl River Delta, China[J]. Particuology, 2018, 40:88-97.
    [5] KISHNE A, YIMAN Y, MORGAN C, et al. Evaluation and improvement of the default soil hydraulic parameters for the Noah Land Surface Model[J]. Geoderma, 2017, 285:247-259.
    [6] DONG Haotian, CAO Shuyang, TAKEMI T, et al. WRF simulation of surface wind in high latitudes[J]. J Wind Eng Ind Aerodyn, 2018, 179:287-296.
    [7] ZEPKA G S, PINTO O, SARAIVA A C V. Lightning forecasting in southeastern Brazil using the WRF model[J]. Atmos Res, 2014, 135/136:344-362.
    [8] GUENTHER A, KARL T, HARLEY P, et al. Estimates of global terrestrial isoprene emissions using MEGAN (model of emissions of gases and aerosols from nature)[J]. Atmos Chem Phys, 2006, 6(11):3181-3210.
    [9] ZHANG Q, STREETS D G, CARMICHAEL G R, et al. Asian emissions in 2006 for the NASA INTEX-B mission[J]. Atmos Chem Phys, 2009, 9(14):5131-5153.
    [10] LI Meng, ZHANG Qiang, KUROKAWA J, et al. MIX:a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP[J]. Atmos Chem Phys, 2017, 17:935-963.
    [11] HUANG Chuangfeng, CHEN Changhong, LI L, et al. The study of emission inventory of anthropogenic air pollutants and VOC species in the Yangtze River Delta region, China[J]. Atmos Chem Phys, 2011, 11(1):4105-4120.
    [12] LI Kangwei, CHEN Linghong, YING Fang, et al. Meteorological and chemical impacts on ozone formation:a case study in Hangzhou, China[J]. Atmos Res, 2017, 196(1):40-52.
    [13] CARTER W P L. Development of ozone reactivity scales for volatile organic compounds[J]. Air Waste, 1994, 44(7):881-899.
    [14] WANG Qiaoli, LI Sujing, DING Minli, et al. VOC emission characteristics and priority control analysis based on VOC emission inventories and ozone formation potentials in Zhoushan[J]. Atmos Environ, 2018, 182:234-241.
    [15] BOYLAN J, RUSSEL A. PM and light extinction model performance metrics, goals, and criteria for three-dimensional air quality models[J]. Atmos Environ, 2006, 40(26):4946-4959.
    [16] CHO S, MORRIS R, MCEACHEM P, et al. Emission sources sensitivity study for ground-level ozone and PM2.5 due to oil sands development using air quality modeling system:part Ⅱ-source apportionment modeling[J]. Atmos Environ, 2012, 55:542-556.
  • [1] 贺丹, 原江琴, 高振飞, 李朝梅, 雷雅凯, 刘艺平, 董娜琳.  郑州市13种常绿灌木叶面微结构对大气颗粒物滞留能力的影响 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20240415
    [2] 吴沁娇, 宋艳冬, 陶士杰, 王丽, 周如意, 陈武, 潘心禾, 周宇峰, 周国模.  丽水白云国家森林公园5种典型林分挥发性有机化合物释放特征及其与环境因子的关系 . 浙江农林大学学报, 2023, 40(5): 930-939. doi: 10.11833/j.issn.2095-0756.20220676
    [3] 陶士杰, 宋艳冬, 吴沁娇, 王丽, 周如意, 吴电, 周宇峰, 周国模.  毛竹林挥发性有机化合物释放特征及其与环境因子的关系 . 浙江农林大学学报, 2023, 40(5): 940-950. doi: 10.11833/j.issn.2095-0756.20230249
    [4] 胡书婧, 张汝民.  挥发性有机化合物在植物适应胁迫及生理生态中的作用 . 浙江农林大学学报, 2022, 39(6): 1378-1387. doi: 10.11833/j.issn.2095-0756.20220180
    [5] 张天然, 郑文革, 章银柯, 黄芳, 李晓璐, 袁楚阳, 于慧, 晏海, 邵锋.  杭州市临安区4种绿地内细颗粒物中重金属污染特征 . 浙江农林大学学报, 2021, 38(4): 737-745. doi: 10.11833/j.issn.2095-0756.20200558
    [6] 张恒, 孙子健, 张运林, 金森, 王玉霞.  不同距离气象数据对细小可燃物含水率预测模型精度的影响 . 浙江农林大学学报, 2018, 35(3): 529-536. doi: 10.11833/j.issn.2095-0756.2018.03.019
    [7] 蒋冬月, 李永红, 沈鑫.  芸香叶片和花瓣释放挥发性有机物成分及其变化规律 . 浙江农林大学学报, 2018, 35(3): 572-580. doi: 10.11833/j.issn.2095-0756.2018.03.025
    [8] 杨建飞, 宁莉萍, 杨了, 王杰, 钱钰滢, 陈颐萱.  黑壳楠木材构造特征及挥发性有机物成分 . 浙江农林大学学报, 2018, 35(5): 927-934. doi: 10.11833/j.issn.2095-0756.2018.05.018
    [9] 蔡宙霏, 陈雅奇, 许馨露, 王小东, 汪俊宇, 张汝民, 高岩.  4个桂花品种开花进程释放VOCs动态变化分析 . 浙江农林大学学报, 2017, 34(4): 608-619. doi: 10.11833/j.issn.2095-0756.2017.04.006
    [10] 王琦, 刘华红, 王彬, 张汝民, 高岩.  7种槭树释放挥发性有机化合物组分分析 . 浙江农林大学学报, 2016, 33(3): 524-530. doi: 10.11833/j.issn.2095-0756.2016.03.022
    [11] 刘庆, 童森淼, 马建义.  顶空固相微萃取-气相色谱-质谱联用法分析竹醋液挥发性化合物 . 浙江农林大学学报, 2014, 31(2): 308-314. doi: 10.11833/j.issn.2095-0756.2014.02.022
    [12] 王琦, 王丹, 张汝民, 高岩.  日本紫藤开花进程中挥发性有机化合物组分与含量的变化 . 浙江农林大学学报, 2014, 31(4): 647-653. doi: 10.11833/j.issn.2095-0756.2014.04.023
    [13] 刘芳, 许改平, 吴兴波, 丁倩倩, 郑洁, 张汝民, 高岩.  干旱-复水处理对迷迭香挥发性有机化合物释放规律的影响 . 浙江农林大学学报, 2014, 31(2): 264-271. doi: 10.11833/j.issn.2095-0756.2014.02.015
    [14] 张洁, 郭金星, 张汝忠, 王星星, 张晓玲, 刘琳, 侯平, 张汝民.  东魁杨梅果实储藏期挥发性有机化合物成分的变化 . 浙江农林大学学报, 2012, 29(1): 143-150. doi: 10.11833/j.issn.2095-0756.2012.01.024
    [15] 马楠, 周帅, 林富平, 高岩, 张汝民.  5种绿篱植物挥发性有机化合物成分分析 . 浙江农林大学学报, 2012, 29(1): 137-142. doi: 10.11833/j.issn.2095-0756.2012.01.023
    [16] 周帅, 马楠, 林富平, 张汝民, 高岩.  樟树花挥发性有机化合物日动态变化分析 . 浙江农林大学学报, 2011, 28(6): 986-991. doi: 10.11833/j.issn.2095-0756.2011.06.025
    [17] 左照江, 张汝民, 朱金胡, 温国胜, 侯平, 高岩.  冷蒿VOCs对4种植物种子萌发和幼苗生长的影响 . 浙江农林大学学报, 2009, 26(1): 76-82.
    [18] 黄益江, 陈杜满, 刘玉明, 刘淑荣.  二氧化硫污染对马尾松生长影响的相关研究 . 浙江农林大学学报, 1998, 15(2): 127-130.
    [19] 樊后保.  树皮对大气污染和降水酸度的指示作用 . 浙江农林大学学报, 1996, 13(2): 136-140.
    [20] 姜培坤, 蒋秋怡, 徐秋芳, 钱新标, 张钦相.  杉木檫树根际土壤有机化合物研究 . 浙江农林大学学报, 1994, 11(3): 235-240.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2019.04.022

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2019/4/810

图(4) / 表(3)
计量
  • 文章访问数:  3997
  • HTML全文浏览量:  735
  • PDF下载量:  59
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-07-31
  • 修回日期:  2018-10-24
  • 刊出日期:  2019-08-20

后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析

doi: 10.11833/j.issn.2095-0756.2019.04.022
    基金项目:

    国家自然科学基金资助项目 51476144

    作者简介:

    丰睿, 工程师, 博士, 从事能源和环境工程研究。E-mail:11727081@zju.edu.cn

    通信作者: 罗坤, 教授, 博士, 从事能源及环境工程领域复杂多相流和反应的理论及数值模拟研究。E-mail:zjulk@zju.edu.cn
  • 中图分类号: X511

摘要: 后G20时期(2016年9月至2017年12月),臭氧成为杭州夏季首要大气污染物,可吸入颗粒物(PM10)成为杭州冬季首要大气污染物。利用杭州市环境检测朝晖站和环境监测下沙站空气质量数据观测值分析55种挥发性有机化合物(VOC)的臭氧形成潜力,并使用美国国家环境保护局研发的第3代大气质量预估模型(WRF-CMAQ)对杭州市中心PM10的污染数值进行模拟分析,提出了后G20时期PM10和臭氧的最佳管控方案。结果显示:在后G20时期杭州市中心和城郊,芳香烃所占臭氧形成潜力比例最高,其次为烯烃和烷烃,最后为乙炔。在市区和城郊,乙烯、间/对二甲苯、甲苯、丙烯和乙苯依次为生成臭氧最主要的5种挥发性有机化合物,因此管控这5种挥发性有机化合物能最大程度的减少臭氧质量浓度。2017年杭州市中心春夏秋冬四季的PM10中分别有50%,32%,48%和45%源于杭州市以外污染源的长距离跨地域传输。其中,春秋冬三季通过长距离传输来到杭州市区的PM10主要源于市区的北方,夏季则主要源于市区的西南方。2017年,杭州市区的交通源、工业源、生活源与农业源对本地排放PM10的贡献率分别为62.6%,27.8%,7.3%和2.3%。当杭州本地的工业源和交通源在2016年的基础上分别减排15%和5%时,杭州市中心的PM10年均质量浓度可低于70 μg·m-3,达到二类区的标准。

English Abstract

丰睿, 黄成臣, 高寒, 郑慧君, 申亚梅, 罗坤. 后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析[J]. 浙江农林大学学报, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
引用本文: 丰睿, 黄成臣, 高寒, 郑慧君, 申亚梅, 罗坤. 后G20时期杭州市挥发性有机物和可吸入颗粒物的特征分析[J]. 浙江农林大学学报, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
FENG Rui, HUANG Chengchen, GAO Han, ZHENG Huijun, SHEN Yamei, LUO Kun. Decoding VOCs and PM10 for Hangzhou in the post-G20 period[J]. Journal of Zhejiang A&F University, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
Citation: FENG Rui, HUANG Chengchen, GAO Han, ZHENG Huijun, SHEN Yamei, LUO Kun. Decoding VOCs and PM10 for Hangzhou in the post-G20 period[J]. Journal of Zhejiang A&F University, 2019, 36(4): 810-817. doi: 10.11833/j.issn.2095-0756.2019.04.022
  • 二十国集团(G20)领导人第十一次峰会以后,杭州市如何治理和管控大气污染存在巨大挑战。2017年全年数据显示,杭州市中心可吸入颗粒物(PM10)质量浓度年均值为76.1 μg·m-3,超过GB 3095-2012《环境空气质量标准》国家环境二类标准8.7%,其中2017年1和2月PM10质量浓度均值达107.2 μg·m-3,12月均值高达124.9 μg·m-3,成为年度冬季首要污染物;夏季8 h最大值臭氧的平均质量浓度高达150.2 μg·m-3,成为年度夏季首要污染物。为了在2022年亚运会前空气质量达到国家环境二级标准,本研究分析了杭州市2016年9月至2017年12月夏季首要大气污染物臭氧最主要前体物挥发性有机化合物(VOC)的臭氧生成潜力(OPF),并运用WRF-CMAQ空气质量模型拟合分析了冬季首要污染物PM10的时空分布、长距离跨地域传输率和本地不同污染源贡献率,旨在提出杭州市切实可行的大气污染防治方案。

    • 数据来源于杭州市环境监测站中的朝晖站和下沙站空气质量观测值。朝晖站(30.288°N, 120.159°E)位于杭州市中心的朝晖五区朝晖实验小学旁,测量高度约18 m,距离杭州市中心地标性建筑环球中心和西湖文化广场约0.8 km,东侧0.3 km为杭州最主要的南北向主干道上塘高架,西侧和北侧0.1 km为步行街。周围污染源主要是生活源和交通源,周围5.0 km内无大型工业污染源,反映了杭州市中心的污染状况。下沙站(30.306°N, 120.348°E)位于杭州下沙科教园区的浙江育英职业技术学院旁,距离老杭州市区约11.6 km,测量高度约20 m,东侧距离杭州绕城高速约1.3 km,北侧距离德胜快速路约1.8 km,南侧距离下沙路约0.1 km,主要污染源为生活源和交通源,反映了杭州市城郊结合部的污染状况。观测时段为2016年7月至2017年12月,有效数据率为98.7%。下沙站位于朝晖站正东方向约18.3 km。由于2017年全年杭州市的东风频率约12%,西风仅为7%,所以两站之间的污染物区域传输相对较小。

      环境监测数据包括每小时的臭氧(O3)、二氧化氮(NO2)、二氧化硫(SO2)、细颗粒物(PM2.5)和可吸入颗粒物(PM10)质量浓度以及日均共计55种VOC的体积分数。VOC的测量仪器为荷兰SYNSPEC公司生产的Syntech Spectras GC955在线气相色谱仪。其他污染物数据由美国Thermo Fisher Scientific公司生产的49i型臭氧分析仪、42i型氮氧化物分析仪、43i型二氧化硫分析仪以及PDM-3700型PM检测仪测量。为确保数据的准确性,隔半个月对所有检测仪器进行1次维护。气象数据包括每小时的气温、相对湿度、风速、风向、降水量、气压和露点温度及每日的日照时数等,数据来源于杭州市气象观测站。

    • 应用美国国家环境保护局研发的第3代气象模型系统的中尺度气象及区域空气质量耦合模型——WRF-CMAQ,研究污染物的长距离跨区域传输及杭州市区内的不同污染源贡献率。模拟过程包括3个步骤:(1)WRF模型提供模拟气象数据,包括气温、相对湿度和风速风向等;(2)污染物清单提供自然源和人为源的排放信息;(3)将步骤1和2中的信息导入CMAQ系统中,计算出某一时间段内污染物的时空分布。WRF-CMAQ的数值模拟计算在浙江大学清洁能源利用国家重点实验室具有2 688个中央处理器(CPU)和5 376个GB的随机存取存储器(RAM)的IBM集群计算机上进行。

    • 气象模拟应用WRFv3.4版本:诸如海拔高度、气压、气温、露点温度和风向风速等地理或气象条件源于美国环境预测中心(NCEP)发布的WRF输入边界层数据库(https://rda.ucar.edu/datasets/ds083.2/),全球高空大气观测数据库(https://rda.ucar.edu/datasets/ds351.0/)以及全球地表气象观测数据库(https://rda.ucar.edu/datasets/ds461.0/)。所有数据库都具有6 h时间分辨率和1°×1°空间分辨率。为确保边界层对污染物模拟的影响最小化,设置WRF模拟的空间比MEGAN和CMAQ略大。表 1列出了气象模拟的WRF参数方案。

      表 1  WRF参数方案

      Table 1.  Parameters for WRF

      参数 选择的方案
      微观物理进程 WSM6[1]
      长波辐射 RRTM[2]
      短波辐射 Goddard[3]
      表层 MM5 similarity[4]
      陆地表面进程 Noah[5]
      边界层 YSU[6]
      积雨云参数化进程 Grell-Devenyi[7]
    • 由陆地生态系统到大气中气体和气溶胶的自然排放清单用自然源气体及气溶胶排放模型MEGAN[8]估算,主要包含自然源产生的氨气,PM10,VOC和可吸入颗粒物等。MEGAN的全球土地覆盖分布数据的空间分辨率为1 km,区域分辨率可达300 m。MEGAN的代码可以免费获取,存储格式为ARC-GIS或net-CDF。运用v2.04版本的MEGAN系统结合土地覆盖数据和WRF模拟出的气象条件,进行自然界生物排放量的估算。土地覆盖数据由3个部分组成:(1)叶面积指数。基于MODIS数据库的土地面积除以叶面积的总面积估算,其空间分辨率为30 s。(2)植物功能类型。来自MODIS3数据库以及全国第2次森林普查;植物分为阔叶、针叶、草本(包括农作物和草地)以及灌木4类。(3)排放因子。不同季节自然源排放污染物总量的差异。

      对于人为污染源清单,在长江三角洲以外地区,运用清华大学在2016年制定的东亚地区的清单[9-10],包括SO2,CO,PM2.5和PM10等26种主要污染物的污染源,所有污染源都被格栅化纳入0.5°×0.5°分辨率的数据库。长江三角洲地区的人为排放源污染物清单,运用上海环境科学院于2016年制定的清单[11]。运用全国范围的清单可详细计算杭州PM10的长距离传输率。

    • CMAQ有5个模块:(1)初始条件处理器,为模拟提供初始浓度场;(2)边界条件处理器,为模拟提供边界浓度场;(3)光解率处理器,提供区域光解速率;(4)气象化学接口处理器,将WRF输出的气象场转换为下一步模拟所需的数据格式;(5)CMAQ化学传输模型,CMAQ的核心,运用大气物理化学模型计算区域传输率以及计算最终的污染物浓度。

      本研究采用的CMAQ版本为5.0.2,结合朗伯投影坐标系进行模拟,运用3层嵌套网格技术:第1层覆盖整个中国,坐标原点为34°N,110°E,网格分辨率为36 km,网格数为173×136;第2层覆盖中国东部,坐标原点为32°N,116°E,网格分辨率为12 km,网格数为135×228;第3层覆盖长江三角洲区域,坐标原点设置在杭州市所在的30°N,120°E,网格分辨率为4 km,网格数为150×174。垂直方向上模拟的高度为100 hPa,网格为14个,越是靠近地面网格越细。分别选择CB05和AERO6作为气相化学反应机理和气溶胶机理的模型参数化方案。

    • 由于杭州市臭氧质量浓度主要受其前体污染物VOC的控制[12],因此通过观测值分析杭州市VOC体积分数的分布情况(图 1)。2016年7月底至9月底,杭州市政府关停了市区绝大多数排放VOC的企业及餐饮店,所以2016年8和9月朝晖站的TVOC体积分数分别仅为7.94和9.63 nL·L-1。同时间,下沙站的TVOC分别为19.53和22.58 nL·L-1。可见杭州市中心的TVOC绝大多数是由工业源及生活源产生。图 1显示:后G20时期,杭州市的TVOC体积分数呈现9-11月偏高,其他月份较为稳定的特点。2017年,朝晖站和下沙站的TVOC年均体积分数分别为27.62和29.53 nL·L-1,市中心和城郊结合部的TVOC质量分数差异不明显。2017年1-8月,两站点TVOC体积分数均保持稳定;而8-9月,两站点的TVOC体积分数均呈现大幅增长态势;8-12月波动起伏较大。根据2017年度杭州市气象数据,夏季TVOC较高是由于7-8月低降雨量减少了湿沉降以及城市热岛效应改变了大气环流并降低了平均风速,产生了不利于污染物扩散的条件;9-12月的波动较大则是由于秋冬季的市政建设产生了大量的TVOC。

      图  1  杭州市朝晖站和下沙站月均TVOC体积分数

      Figure 1.  Monthly average TVOC concentration at Zhaohui and Xiasha stations

    • 将55种VOC分成四大类,其中烷烃29种,烯烃10种,芳香烃15种以及乙炔。2017年全年,朝晖站的烷烃、烯烃、芳香烃和乙炔分别占TVOC的63.1%,11.2%,25.1%和0.6%;下沙站的烷烃、烯烃、芳香烃和乙炔分别占TVOC的65.6%,8.3%,24.8%和1.3%。可见2017年杭州市区和城郊结合部的不同种类VOC在TVOC的占比大致相似。采用天文季节的四季等长划分方法,按照春(3-5月),夏(6-8月),秋(9-11月),冬(12月-翌年2月)划分四季,各季节不同种类VOC的占比如表 2所示。由表 2可见:除了春季的烷烃和烯烃之外,市区和城郊的不同属VOC排放占比并无太大差距。杭州市区春季烷烃质量浓度较少而烯烃偏高,城郊则正好相反。不同种类VOC的大气反应活性并不相同。单位体积分数的VOC生成臭氧的能力称为臭氧形成潜力(F)。

      $$ F = \sum\limits_{i = 1}^{\rm{N}} {{V_i}} \times R。 $$ (1)

      表 2  2017年度各季节不同种类VOC排放比例

      Table 2.  Seasonal contribution of different subgroup VOC in 2017

      季节 烷烃排放比例/% 烯烃排放比例/% 芳香烃排放比例/% 炔烃排放比例/%
      朝晖站 下沙站 朝晖站 下沙站 朝晖站 下沙站 朝晖站 下沙站
      59.7 70.7 15.4 5.5 24.0 22.6 0.9 1.2
      73.1 69.3 9.6 8.0 16.6 21.3 0.7 1.4
      62.6 58.0 7.9 10.0 28.9 30.9 0.6 1.1
      60.4 65.2 10.2 8.8 28.9 24.4 0.51.6

      式(1)中:Vi为第i种VOC的体积分数占比;R为该种VOC的最大增量反应活性,由CARTER[13]测量并估算。某一大类的VOC的总臭氧形成潜力等于所有种类VOC的臭氧形成潜力之和[14]。结果表明:后G20时期朝晖站的总臭氧形成潜力比下沙站高出约7.4%,即杭州市中心和城郊的TVOC对臭氧生成的贡献大致相等。朝晖站的烷烃、烯烃、芳香烃和乙炔分别占该站总臭氧形成潜力的24.1%,41.5%,34.3%和0.1%;下沙站的烷烃、烯烃、芳香烃和乙炔分别占该站总臭氧形成潜力的28.1%,34.4%,37.3%和0.2%。由此可见,在市区烯烃是生成臭氧最主要的前体物,其次为芳香烃;在城郊芳香烃是生成臭氧最主要的前体物,其次为烯烃。虽然烷烃在2个站点的质量浓度都是最高的,但是对于生成臭氧的贡献却比烯烃和芳香烃少。此外,后G20时期,在杭州市区和城郊间/对二甲苯、甲苯、乙烯、乙苯和丙烯依次为臭氧形成潜力值最大的单个种类VOC,所以控制这5种VOC对于降低杭州臭氧质量浓度至关重要。

    • 运用WRF-CMAQ大气模型,模拟了分别代表冬、春、夏、秋的2017年1,4,7和10月杭州市朝晖站PM10的质量浓度,以评估杭州市中心区域四季不同污染物跨区域传输比例以及不同本地排放源对本地排放PM10的贡献率。

    • 用分数偏差均值(B)和分数误差均值(E)评估WRF-CMAQ模拟的优劣:

      $$ B=\frac{1}{N} \sum\limits_{i=1}^{N} \frac{S-O}{O+S / 2} \times 100 \%; $$ (2)
      $$ E=\frac{1}{N} \sum\limits_{i=1}^{N} \frac{|S-O|}{O+S / 2} \times 100 \%。 $$ (3)

      式(2)~(3)中:n为模拟的次数;S为模拟值;O为观测值。BOYLAN等[15]定义:当模拟值B≤±30%和E≤+50%时,该大气模型可以被认为是精确的。对于朝晖站PM10的模拟结果显示:4个月的模拟值均为B≤±20%和E≤+35%,因此模拟的精度达到了国际上定量的要求。

      图 2为2017年1,4,7,10月的平均PM10质量浓度的时空分布,覆盖了包括杭州市在内的长江三角洲地区。由图 2可见:杭州市的PM10污染冬季最为严重,市区范围内出现了2个PM10质量浓度大于130 μg·m-3的红色极值区域;春季PM10污染开始缓解,夏季污染最轻,秋季PM10质量浓度又开始回升。红色区域显示杭州市区PM10质量浓度极大值区域,为市中心与位于城东的大江东工业园区和建筑工地,与实际情况相符。从长江三角洲范围来看,PM10质量浓度的高值区域位于工业城市上海和苏州及其附近,而位于杭州市南方的浙南山区普遍PM10较低。

      图  2  2017年1,4,7,10月长江三角洲PM10质量浓度的时空分布示意图

      Figure 2.  Time-spacing map of PM10 over Yangtze River Delta in 2017

    • 采用置零法[16],通过关闭杭州市本地总污染源清单以及分别关闭本地不同污染源清单的方式,来评估2017年1,4,7,10月长距离传输和本地源对PM10的贡献率。长距离传输的计算公式为:

      $$ {P_i} = \frac{{{C_{{\rm{zero}}}}}}{{{C_{\rm{b}}}}} \times 100\%。 $$ (4)

      式(4)中:Pi为第i个月长距离传输所占的比例;Cb为打开全国所有污染源的模拟值;Czero为关闭杭州本地总污染源清单后的模拟值。结果显示:在杭州市中心区域,2017年1,4,7和10月的PM10中分别有45%,50%,32%和48%源于杭州市之外的污染源长距离跨区域传输。跨区域传输的PM10除7月外,其他3个月都接近或达50%。因此要防控PM10,仅仅控制杭州本地的污染源是不够的。图 3显示:当关停杭州市本地PM10污染源之后,杭州市区的PM10污染程度有所下降,但仍较为严重。图 4显示:2017年杭州市秋冬季节的主导风向为正北风,分别占总风向频率的40%和28%;夏季40%为西南风,且平均风速全年最低;春季风向较为多变,其中正北方风向占比20%。由于杭州市PM10主要排放源皆位于市区的北方,因此春秋冬三季PM10通过气团水平位移的长距离区域传输比例(45%~50%)比夏季(32%)要高,与模拟值相符:春秋冬季跨区域传输的PM10主要源于杭州市区的北方,夏季长距离传输的PM10主要源于市区的西南方。图 2中PM10各月的扩散方向符合图 4各季的风向。

      图  3  2017年1,4,7,10月关停杭州市本地污染源后的PM10时空分布示意图

      Figure 3.  Time-spacing map of PM10 in 2017 when Hangzhou local sources are shut down

      图  4  杭州市2017年风玫瑰图

      Figure 4.  Wind-rose over Hangzhou in 2017

      WRF-CMAQ清单中,把污染源分成农业源(包括秸秆燃烧)、工业源、生活源和交通源。通过置零法逐个关停本地4类污染源,分析杭州市本地污染源对PM10的贡献。计算公式为:

      $$ {P_j} = \frac{{{C_{\rm{b}}} - {C_j}}}{{{C_{\rm{b}}} - {C_{{\rm{zero}}}}}} \times 100\%。 $$ (5)

      式(5)中:Pj为第j个本地污染源对PM10的贡献率;Cj为关停第j个本地污染源后的PM10质量浓度。如表 3所示:2017年全年杭州市中心区域的4类污染源贡献率基本稳定。对PM10贡献最大的为工业源,占62.6%;其次是交通源,占27.8%。因此,要控制杭州市区的PM10污染,应着重控制交通源和工业源。计算显示,当杭州本地的工业源和交通源在2016年的基础上分别减排15%和5%时,杭州市中心PM10的年均质量浓度可低于70 μg·m-3,达到GB 3095-2012二类区标准。

      表 3  杭州本地污染源对市中心PM10的贡献率

      Table 3.  Hangzhou local sector contribution for PM10 at city center

      月份贡献率/%
      农业源 工业源 生活源 交通源
      1 2.6 61.1 8.4 27.9
      4 2.1 62.6 6.7 28.6
      7 3.0 60.2 7.2 29.6
      10 1.5 66.4 6.8 25.3
      平均 2.3 62.6 7.327.8
    • 本研究通过观测值分析,并结合WRF-CMAQ的数值模拟,提出了杭州市后G20前亚运会时期,管控市区最主要大气污染物臭氧和PM10的方法。结论如下:在杭州市中心,烯烃是生成臭氧的最主要的VOC大类,其次为芳香烃和烷烃,最后为炔烃;在城郊结合部,芳香烃是生成臭氧的最主要的VOC大类,其次为烯烃和烷烃,最后为炔烃。其中,乙烯、间/对二甲苯、丙烯和乙苯依次为生成臭氧最主要的5种VOC。由于杭州的臭氧主要受到VOC的影响,管控这5种VOC是降低杭州市夏季臭氧质量浓度最有效的方法。2017年杭州市中心春夏秋冬四季的PM10中分别有50%,32%,48%和45%源于杭州市之外污染源的长距离传输。其中,春秋冬三季通过长距离传输来到杭州市区的PM10主要源于市区的北方,而夏季主要源于市区的西南方。2017年杭州市区的工业源、交通源、生活源与农业源对杭州本地排放的PM10的贡献率分别为62.6%,27.8%,7.3%和2.3%。当杭州本地的工业源和交通源在2016年的基础上分别减排15%和5%时,杭州市中心的PM10年均质量浓度可低于70 μg·m-3,达GB 3095-2012对二类区要求。

      由于TVOC约1 200种,在空气污染物中常见的约100~110种。因此,未来的工作重点是,除了对含碳氢元素的烃类VOC进行测量,也应对酮类、醛类、醚类等含氧VOC,以及卤代烃类的VOC进行测量,从而可以提出更完备的臭氧防控方案。此外,用于此次模拟的PM10清单的空间分辨率为4 km,虽然模拟的精度达到了国际标准,但是对杭州城市内小范围的模拟还显粗糙。因此,如何通过制作分辨率更高的清单以及实时更新城市下垫面数据,计算杭州主城区内不同街道之间的区域传输率是下一步工作的难点和重心。

参考文献 (16)

目录

    /

    返回文章
    返回