-
二十国集团(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参数方案。
-
由陆地生态系统到大气中气体和气溶胶的自然排放清单用自然源气体及气溶胶排放模型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。
-
将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.5 1.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较低。
-
采用置零法[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
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.3 27.8
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,达到二类区的标准。Abstract: In the post-G20 period (September 2016 to December 2017), the dominant atmospheric pollutant in winter was PM10 and in summer was tropospheric ozone. For the purpose of better regulating PM10 and ground-level ozone, chemical reactivity of 55 kinds of volatile organic compounds (VOCs) was analyzed at Zhaohui and Xiasha Environmental Monitoring Stations using ozone formation potential (OFP). The third generation atmospheric model WRF-CMAQ developed by the United States Environmental Protection Agency (USEPA) was also used to simulate time-space distributions of PM10 in the city center. Results showed that alkenes were the primary precursor in the downtown area, followed by aromatics, alkanes, and alkynes; whereas, aromatics were the primary precursor in the suburbs, followed by alkenes, alkanes, and alkynes. Ethylene, m/p-xylene, toluene, propylene, and ethylbenzene were ranked as the top five VOCs for OFP in both downtown and suburban areas of Hangzhou. For PM10 in the city center, regional transport of pollutants accounted for about 50% in spring, 32% in summer, 48% in autumn, and 45% in winter. The regional transported PM10 mainly came from the north in spring, autumn, and winter but from the southwest in summer. Local Hangzhou PM10 levels in 2017 included local industry (62.6%), traffic (27.8%), residence (7.3%), and agriculture (2.3%). Consequently, controlling the top five VOC species was the best way to alleviate ground-level ozone, and based on 2016 statistics, shutting down 15% of the local industry and reducing 5% of the local traffic would reduce the annual PM10 level to 70 μg·m-3 meeting the second national environmental air quality standard.
-
表 1 WRF参数方案
Table 1. Parameters for WRF
表 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.5 1.6 表 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.3 27.8 -
[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. -
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2019.04.022