2. 西安市环境监测站, 西安 710119;
3. 陕西省环境监测中心, 西安 710054
2. Xi'an Environmental Monitoring Station, Xi'an 710119, China;
3. Shaanxi Environmental Monitoring Center, Xi'an 710054, China
碳组分[主要指元素碳(EC)和有机碳(OC)]是大气颗粒物的重要组成部分, 影响光的散射与吸收, 为大气化学过程提供载体和催化活性点位, 对环境空气质量、大气能见度、全球气候变化、身体健康[1~3]都有着重要的影响.目前我国已开始对部分城市的碳组分展开研究, 主要集中在量浓度、时空分布特征、二次有机碳估算、来源解析以及对人体健康影响等方面[4~7], 但对于人口密集的关中区域却研究较少, 因此, 有必要对关中地区碳组分污染特征展开研究.
本研究通过对2017年秋冬季关中地区细颗粒物(PM2.5)中OC、EC进行观测, 分析颗粒物中碳组分特征, 采用PMF受体模型进行源解析, 对认识关中地区碳组分的区域污染特征及形成机制有着重要的意义, 以期为关中地区颗粒物污染控制提供科学依据.
1 材料与方法 1.1 样品采集本研究选取关中地区的西安、渭南、铜川、宝鸡、咸阳为典型城市代表, 并在西安城市运动公园、渭南师范学院、铜川新区环境保护局、宝鸡市环境保护局、陕西中医药大学设置采样点见表 1.秋冬季采样时间分别为2017年9月4日至11月7日、2017年11月13日至2018年1月19日, 采样周期14 d(当日10:30至次日08:30).
![]() |
表 1 采样点位情况 Table 1 Details for the sampling sites |
采样仪器为TH-150F型智能中流量采样器(武汉天虹仪表有限责任公司, 中国), 采样流量100 L·min-1.采样滤膜使用石英滤膜(Φ90 mm, Whatman QMA, 英国).执行《环境空气颗粒物(PM2.5)手工监测方法(重量法)技术规范》(HJ 656-2013)标准.气象条件采用当地气象局同期数据.
1.2 样品分析 1.2.1 样品称重分析采样前后滤膜的称重使用精度为十万分之一(10 μg)的电子天平(Mettler Toledo XAl05, 瑞士)进行称量, 称量条件为恒温(25℃±1℃)、恒湿(40%±5%), 平衡时间为24 h, 每张滤膜多次称量取平均值, 保证误差值小于±40 μg.颗粒物浓度为颗粒物质量除以标况体积(0℃, 101.3 kPa).
1.2.2 碳组分分析碳组分分析使用Model 5L-NDIR型OC/EC分析仪(Sunset Laboratory Inc), 采用热光透射法(Thermal-Optical Transmission, TOT), 遵循美国EPA推荐的NIOSH870分析方法[8].
2 结果与讨论 2.1 细颗粒物浓度分布特征关中地区观测期间各点位细颗粒物及OC、EC浓度如图 1所示, PM2.5的平均浓度为(87.39±46.54) μg·m-3, 与我国《环境空气质量标准》(GB 3095-2012)规定的二级日均标准(PM2.5为75 μg·m-3)相比, 整体来说受到细颗粒物的污染较为严重, XY、XA和WN的浓度分别是标准浓度的1.39、1.38和1.23倍, 而TCH和BJ点位的颗粒物浓度均较低, 环境空气质量较好.各点位PM2.5浓度的空间分布为XY>XA>WN>BJ>TCH, XA与XY污染物浓度水平相当, 这与监测点位的地理位置相近、存在着相同的扩散条件和区域性污染结构密切相关.
![]() |
图 1 关中地区细颗粒物及其OC、EC浓度分布特征 Fig. 1 Mass concentrations of PM2.5, OC, and EC in the Guanzhong area |
观测期间关中地区大气PM2.5颗粒物中OC和EC的平均浓度分别为(14.48±7.86) μg·m-3和(2.27±0.95) μg·m-3, OC和EC浓度的空间分布与细颗粒物基本相同. XY点位碳组分浓度最高, 细颗粒中OC浓度达到(17.8±9.66) μg·m-3, EC浓度达到(2.5±1.18) μg·m-3, 该点位东临沣河, 植被丰富, 北部有交通要道, 且咸阳工业较多, 明显存在多种污染源贡献.
2.2.2 OC、EC浓度国内外对比元素碳和有机碳目前没有标准的浓度限值, 可以相比于其它城市判断关中地区的污染状况, 如表 2所示, 关中地区PM2.5中OC和EC的浓度低于较北的城市济南、北京和天津, 北方城市冬季需要供暖, 污染排放量增加, 碳组分的浓度普遍较高.而天津是重要的工业城市, 化工行业、制造业等第二产业发达, 产业链接多为汽车集疏运输形式实现[9], 北京和济南第二产业多, 有着密集的公路和航空线路, 燃料燃烧污染和交通污染可能是这些城市OC和EC浓度高的主要原因.关中地区高于南方及沿海城市广州、香港及国外城市莱切(意大利)、洛杉矶(美国), 关中地区火电和煤化工行业集聚, 主要以煤炭消耗为主的单一能源结构, 煤炭消费量占终端能源消费量的比重为72.5%[10], 加之盆地地形结构, 不易于扩散, 碳组分污染严重.本研究监测浓度低于田鹏山等[11]于2012年冬季监测的数据, 低于王帆等[12]2013年西安市秋冬季监测数据, 虽然碳组分测定方法及升温协议上存在差别, 但其总碳组分(TC)的含量基本上相同[13], 说明近些年来碳组分的污染程度明显改善, 这与关中地区区域联控政策和一系列大气污染防治工作如燃煤锅炉淘汰改造、清洁取暖及燃煤替代、散乱污企业整治、机动车限行和建筑施工扬尘管控等密切相关.
![]() |
表 2 关中地区OC和EC浓度与国内外城市对比1) Table 2 Comparison of OC and EC concentrations with other cities |
2.2.3 OC和EC在PM2.5中占比及季节分布
OC和EC在细颗粒物中的占比如表 3所示, 比例范围分别为6.54%~31.54%和1.08%~7.10%, 平均占比分别为18.04%和2.99%. TCH点位碳组分浓度及占比均低于其它点位, 这与铜川的海拔高度与地理位置有关, 关中地区为南北两山夹一谷的地形, 主要风向为东北风, 秋冬季地面风速低, 弱东风可将东部部分污染传输至西部, 即WN、XY、XA和BJ点位方向的污染物积累, 又由于关中地区的口袋形地形, 没有西北强风的影响, 污染物将在盆地内徘徊, 不易扩散, 呈明显的中部浓度高, 周边相对较低特征, TCH点位所受影响较小.
![]() |
表 3 关中地区OC和EC在PM2.5中平均占比1)/% Table 3 Percentage of OC and EC in PM2.5 for the Guanzhong area/% |
关中地区各点位秋冬季节的颗粒物污染特征如图 2和图 3所示, 各点位OC和EC浓度季节分布明显, 呈秋季低冬季高的特点.冬季混合层降低, 环境容量减少, 污染物排放量增加, 不易扩散导致污染加重.且冬季低挥发性前体物低温情况下可在气相中凝聚并凝结于颗粒物表面, 同时可溶性气体在颗粒物中的凝结也可增加OC的浓度. PM2.5碳组分占比中, XA和WN点位秋冬季变化较小, 而XY冬季细颗粒碳组分占比明显下降, 可能原因是观测期间降雪天气居多, 加之EC具有良好的表面吸附活性, 比表面积较大的雪花在降落过程中更易黏附或者凝结碳气溶胶, 沉降于地面从而减少碳气溶胶的含量, 同时降雪天气常常伴随着冷锋过境的影响[21].
![]() |
图 2 OC和EC在PM2.5中占比的季节分布 Fig. 2 Seasonal variation of OC and EC ratios in PM2.5 |
![]() |
图 3 PM2.5及其OC和EC浓度的季节分布 Fig. 3 Seasonal variation of OC, EC, and PM2.5concentrations |
不同污染源排放的颗粒物中OC与EC比值不同, 因此在一定程度上能够定性判断污染的来源[22].有研究结果表明, 烟煤/沥青煤燃烧的OC/EC的值为0.32~13.6, 无烟煤燃烧的比值为1.49~14.8, 平均比值为6.3[23, 24], 机动车尾气污染为4.1[25], 由于生物质燃烧的种类较多, 国内外对生物质燃烧污染的OC/EC值存在差异, Cachier等[26]的研究结果为9, Chen等[24]得到国内部分地区生物质燃烧中OC/EC的比值为13~15, 关中地区细颗粒物中OC/EC比值范围为3.14~10.76, 说明关中地区燃煤污染和机动车污染源是主要的污染贡献源.如图 4所示, 各点位颗粒物中OC/EC比值有着明显的差异, 这也与当地不同的污染源结构有关.
![]() |
图 4 PM2.5中OC/EC值分布 Fig. 4 Distribution of OC/EC ratios in PM2.5 |
空间分布上, XY的比值最高, 其次是BJ、XA和WN, TCH的比值最低, 季节分布来看, 除TCH外其它点位冬季明显地高于秋季, TCH点位秋冬季污染结构未发生较为明显的变化, 也可能监测期间有良好的扩散条件.对关中地区细颗粒物中秋冬季OC与EC进行相关性分析可知, 秋季相关系数R2(P < 0.05)为0.799, 冬季R2 (P < 0.05)为0.792, 均表现为显著性相关, 说明关中地区秋冬季细颗粒物中OC与EC有着较为相同的来源.
2.3 碳组分污染来源解析 2.3.1 碳组分特征分析本研究采用NIOSH870协议, 通过热光透射法在不同的温度梯度下能够得到10种碳组分, 在一定程度上能够表现源谱特征, 实现颗粒物中碳组分的来源解析.各点位碳组分浓度分布如图 5所示, 分布顺序为OC1>EC2>EC3>OC4>EC1>OC2>OC3>EC4>EC6>EC5, OC1在碳组分中比例最高为23.44%, EC2和EC3的比例均约为15%, EC5和EC6的比例均小于1%.各点位中组分浓度分布情况较为一致, 比例存在明显的差异, 其中EC3和OC1在TCH点位的比例高, OC2在XA点位的比例最高, EC2在XY点位的比例较高.季节分布中, OC1、OC2和OC3的比例秋冬季基本保持不变, OC4、EC4、EC5和EC6的比例低, 季节特征明显, 秋季高于冬季. EC1、EC2和EC3则表现为秋季低冬季高.
![]() |
图 5 关中地区细颗粒中各类碳组分浓度 Fig. 5 Concentrations of carbon components in PM2.5 |
正矩阵因子分解法(positive matrix factorization, PMF)目前已经广泛地应用于气溶胶及其组分污染来源的定量解析, 本研究采用美国国家环境保护署(EPA)发布的PMF 5.0模型, 输入的组分种类为OC1、OC2、OC3、OC4、EC1、EC2、EC3、EC4、EC5、EC6和TC, 不确定性设置为相应浓度的10%, 运行前需对异常值数据进行排除, 以尽量减少数据误差.
本研究通过对不同因子数量的对比, 得出4类(factor=4)因子运行的结果中Qrobust/Qexp的减少量最小为1.84, 因此选择4类污染因子进行关中地区颗粒物污染来源解析[27]. PMF模型模拟计算的数据与实际监测数据的相关性如表 4所示, 每组物种除EC5外其它R2值均大于0.85, TC的相关性达到0.99, 数据值较为准确.为确定本次运行结果的稳定性, 需对PMF模型结果进行误差评估分析(error etimation, EE), 结果显示, 运行4种污染因子, DISP方法误差分析中, Q值的变化率小于0.1%, 因子的交换个数为0, BS方法误差分析中, 因子映射率(mapping)均大于80%, 说明模型运行结果较为稳定[27].
![]() |
表 4 监测数据与PMF模型模拟计算结果相关性分析 Table 4 Linear fitting between monitoring data and simulation results of the PMF model |
根据不同污染因子的特征谱图可以判断污染因子的类型[8, 28], 观测期间关中地区PM2.5中碳组分各类污染源特征谱图如图 6, 因子1中OC1(52.15%)、OC2(68.17%)、OC3(58.61%)和EC1(50.82%)的贡献较大, 主要是较易挥发的有机物, 将其识别为生物质燃烧与燃煤污染源; 因子2中EC1(44.04%)、EC2(43.25%)和OC1(19.37%)的贡献比例大, 将其识别为汽油车尾气污染源; 因子3中EC3(24.22%)和EC4(59.59%)的贡献比例较大, 将其识别为柴油车污染; 因子4中OC4(15.5%)、EC5(32.84%)和EC6(83.51%)的贡献比例较大, 还有少量的EC4和EC3存在, 将其识别为道路扬尘类污染源, 该类源可能含有碳酸盐类的其它扬尘污染, 也可能为石英膜在高温条件下自身的物质被分解出来造成的, 有待进一步研究.
![]() |
图 6 细颗粒物碳组分不同污染来源成分谱图 Fig. 6 Contributions of different carbonaceous components to factors 1-4 |
监测期间大气颗粒物细颗粒物中碳组分的源解析结果如图 7所示, 饼图为该城市点位观测期间各污染源的平均贡献分布.碳组分主要的污染贡献来源是生物质燃烧与燃煤源, 所占比例为48.63%, 其次是汽油车尾气源(23.07%), 柴油车尾气源(18.82%), 道路扬尘类污染(9.47%).从季节分布来看, 生物质燃烧与燃煤污染在冬季明显增加, 由于冬季供暖开始, 散煤燃烧以及城市供热系统等污染排放量增加.汽油车尾气污染贡献有所增加, 可能是冬季低温条件下机动车冷启动引起的排放量增加, 柴油车尾气排放污染与道路扬尘类污染均有所减少, 可能受到了冬季“禁土”与柴油车限行政策的影响.
![]() |
图 7 细颗粒物中碳组分污染来源PMF解析结果 Fig. 7 Results of the PMF modeling of carbon components and pollution sources |
从点位分布来看, XA、WN、BJ和XY点位生物质燃烧与燃煤是主要贡献源, 其中XY影响最大, 而TCH点位主要是机动车尾气源污染, 其受到的道路扬尘污染贡献比例最大为13.28%, 因扩散条件良好, TCH点位污染物混合及化学反应条件不足, 直接排放和尘类污染影响相对较大. XY点位秋冬季生物质与燃煤燃烧源贡献变化相对较小, 污染程度均高. BJ点位冬季生物质与燃煤燃烧贡献最大, 可能受到污染传输累积于此.不同城市点位的污染结构有所不同, 这与城市的地理位置、气象扩散条件、发展水平、能源结构等有很大的关系.
3 结论(1) 观测期间关中地区5个点位细颗粒物浓度的空间分布差异明显, XY、WN和XA观测点位颗粒物浓度均超标, 而BJ和TCH空气质量相对较好.细颗粒物中OC和EC的浓度分别为(14.48±7.86) μg·m-3和(2.27±0.95) μg·m-3, 占比分别为18.04%和2.99%, 与其他城市相比污染较为严重.碳组分含量的空间分布为XY>WN>XA>BJ>TCH, 且季节差异明显冬季含量高于秋季.
(2) 关中地区细颗粒物中OC/EC比值范围为3.14~10.76, 说明燃煤污染和机动车污染源是主要的污染贡献源, OC与EC的相关性显著, 有着较为相同的污染来源.
(3) 关中地区细颗粒物中各点位各类碳组分的平均浓度分布序列为OC1>EC2>EC3>OC4>EC1>OC2>OC3>EC4>EC6>EC5, OC1组分比例最高为23.44%, EC2和EC3的比例约为15%. PMF模型源解析结果表明, 该地区碳组分污染来源主要有生物质燃烧与燃煤污染源、汽油车尾气污染、柴油车尾气污染、道路扬尘污染源4类, 其贡献率分别为48.63%、23.07%、18.82%和9.47%, 其中铜川受到的机动车尾气污染源的影响最大, 其它点位生物质燃烧与燃煤污染贡献最大.燃煤类污染源贡献季节变化最为明显, 这与关中地区供暖有直接关系.
[1] | 竹涛, 徐东耀, 于妍. 大气颗粒物控制[M]. 北京: 化学工业出版社, 2013. |
[2] |
张承中, 丁超, 周变红, 等. 西安市冬、夏两季PM2.5中碳气溶胶的污染特征分析[J]. 环境工程学报, 2013, 7(4): 1477-1481. Zhang C Z, Ding C, Zhou B H, et al. Pollution characteristics of carbon aerosol in PM2.5 of Xi'an during winter and summer[J]. Chinese Journal of Environmental Engineering, 2013, 7(4): 1477-1481. |
[3] | Sharma S, Brook J R, Cachier H, et al. Light absorption and thermal measurements of black carbon in different regions of Canada[J]. Journal of Geophysical Research:Atmospheres, 2002, 107(D24): 4771. DOI:10.1029/2002JD002496 |
[4] | 陶李.西安市大气中PM2.5及其碳气溶胶污染特征研究[D].西安: 西安建筑科技大学, 2011. |
[5] | Cao J J, Lee S C, Ho K F, et al. Spatial and seasonal variations of atmospheric organic carbon and elemental carbon in Pearl River Delta Region, China[J]. Atmospheric Environment, 2004, 38(27): 4447-4456. DOI:10.1016/j.atmosenv.2004.05.016 |
[6] |
唐小玲, 毕新慧, 陈颖军, 等. 不同粒径大气颗粒物中有机碳(OC)和元素碳(EC)的分布[J]. 环境科学研究, 2006, 19(1): 104-108. Tang X L, Bi X H, Chen Y J, et al. Study on the size distribution of organic carbon (OC) and elemental carbon (EC) in the aerosol[J]. Research of Environmental Sciences, 2006, 19(1): 104-108. DOI:10.3321/j.issn:1001-6929.2006.01.026 |
[7] |
杜博涵, 黄晓锋, 何凌燕, 等. 宁波市PM2.5中碳组分的时空分布特征和二次有机碳估算[J]. 环境科学, 2015, 36(9): 3128-3134. Du B H, Huang X F, He L Y, et al. Seasonal and spatial variations of carbon fractions in PM2.5 in Ningbo and the estimation of secondary organic carbon[J]. Environmental Science, 2015, 36(9): 3128-3134. |
[8] | Bian Q J, Alharbi B, Shareef M M, et al. Sources of PM2.5 carbonaceous aerosol in Riyadh, Saudi Arabia[J]. Atmospheric Chemistry and Physics, 2018, 18(6): 3969-3985. DOI:10.5194/acp-18-3969-2018 |
[9] | 宋天月.浅谈天津市主导产业[DB/OL]. https://wenku.baidu.com/view/70561f780029bd64793e2c5b.html, 2016-03-07. |
[10] | 陕西省统计局.新常态下陕西能源生产平稳发展能源消费结构改善[DB/OL]. http://www.shaanxitj.gov.cn/site/1/html/126/131/138/17854.htm, 2018-05-23. |
[11] |
田鹏山, 曹军骥, 韩永明, 等. 关中地区冬季PM2.5中碳气溶胶的污染特征及来源解析[J]. 环境科学, 2016, 37(2): 427-433. Tian P S, Cao J J, Han Y M, et al. Pollution characteristics and sources of carbonaceous aerosol in PM2.5 during winter in Guanzhong area[J]. Environmental Science, 2016, 37(2): 427-433. |
[12] |
王帆, 韩婧, 张俊, 等. 西安城区大气PM2.5中有机碳与元素碳的污染特征[J]. 环境保护科学, 2015, 41(2): 80-85. Wang F, Han J, Zhang J, et al. Simple analysis of pollution characteristics of organic carbon and elemental carbon in PM2.5 in the atmosphere of Xi'an city[J]. Environmental Protection Science, 2015, 41(2): 80-85. DOI:10.3969/j.issn.1004-6216.2015.02.017 |
[13] | Chow J C, Watson J G, Crow D, et al. Comparison of IMPROVE and NIOSH carbon measurements[J]. Aerosol Science and Technology, 2001, 34(1): 23-34. DOI:10.1080/02786820119073 |
[14] | Ji D S, Zhang J K, He J, et al. Characteristics of atmospheric organic and elemental carbon aerosols in urban Beijing, China[J]. Atmospheric Environment, 2016, 125: 293-306. DOI:10.1016/j.atmosenv.2015.11.020 |
[15] | Zhou S Z, Yang L X, Gao R, et al. A comparison study of carbonaceous aerosols in a typical North China Plain urban atmosphere:seasonal variability, sources and implications to haze formation[J]. Atmospheric Environment, 2017, 149: 95-103. DOI:10.1016/j.atmosenv.2016.11.009 |
[16] | Lai S C, Zhao Y, Ding A J, et al. Characterization of PM2.5 and the major chemical components during a 1-year campaign in rural Guangzhou, Southern China[J]. Atmospheric Research, 2016, 167: 208-215. DOI:10.1016/j.atmosres.2015.08.007 |
[17] | Zhou S Z, Wang T, Wang Z, et al. Photochemical evolution of organic aerosols observed in urban plumes from Hong Kong and the Pearl River Delta of China[J]. Atmospheric Environment, 2014, 88: 219-229. DOI:10.1016/j.atmosenv.2014.01.032 |
[18] |
元洁, 刘保双, 程渊, 等. 2017年1月天津市区PM2.5化学组分特征及高时间分辨率来源解析研究[J]. 环境科学学报, 2018, 38(3): 1090-1101. Yuan J, Liu B S, Cheng Y, et al. Study on characteristics of PM2.5 and chemical components and source apportionment of high temporal resolution in January 2017 in Tianjin urban area[J]. Acta Scientiae Circumstantiae, 2018, 38(3): 1090-1101. |
[19] | Cesari D, Merico E, Dinoi A, et al. Seasonal variability of carbonaceous aerosols in an urban background area in Southern Italy[J]. Atmospheric Research, 2018, 200: 97-108. DOI:10.1016/j.atmosres.2017.10.004 |
[20] | Shirmohammadi F, Hasheminassab S, Saffari A, et al. Fine and ultrafine particulate organic carbon in the Los Angeles basin:trends in sources and composition[J]. Science of the Total Environment, 2016, 541: 1083-1096. DOI:10.1016/j.scitotenv.2015.09.133 |
[21] | 热比古丽·达木拉.乌鲁木齐大气多环芳烃、有机碳/元素碳粒径及季节性分布[D].乌鲁木齐: 新疆大学, 2013. |
[22] | Turpin B J, Huntzicker J J. Identification of secondary organic aerosol episodes and quantitation of primary and secondary organic aerosol concentrations during SCAQS[J]. Atmospheric Environment, 1995, 29(23): 3527-3544. DOI:10.1016/1352-2310(94)00276-Q |
[23] | Tian J, Ni H Y, Cao J J, et al. Characteristics of carbonaceous particles from residential coal combustion and agricultural biomass burning in China[J]. Atmospheric Pollution Research, 2017, 8(3): 521-527. DOI:10.1016/j.apr.2016.12.006 |
[24] | Chen Y J, Zhi G R, Feng Y L, et al. Measurements of emission factors for primary carbonaceous particles from residential raw-coal combustion in China[J]. Geophysical Research Letters, 2006, 33(20). |
[25] | Cao J J, Wu F, Chow J C, et al. Characterization and source apportionment of atmospheric organic and elemental carbon during fall and winter of 2003 in Xi'an, China[J]. Atmospheric Chemistry and Physics, 2005, 5(11): 3127-3137. DOI:10.5194/acp-5-3127-2005 |
[26] | Cachier H, Brémond M P, Buat-Ménard P. Carbonaceous aerosols from different tropical biomass burning sources[J]. Nature, 1989, 340(6232): 371-373. DOI:10.1038/340371a0 |
[27] | Brown S G, Eberly S, Paatero P, et al. Methods for estimating uncertainty in PMF solutions:examples with ambient air and water quality data and guidance on reporting PMF results[J]. Science of the Total Environment, 2015, 518. |
[28] | Watson J G, Chow J C, Lowenthal D H, et al. Differences in the carbon composition of source profiles for diesel- and gasoline-powered vehicles[J]. Atmospheric Environment, 1994, 28(15): 2493-2505. DOI:10.1016/1352-2310(94)90400-6 |