环境科学学报  2016, Vol. 36 Issue (2): 630-637
2014年1月北京市大气重污染过程单颗粒物特征分析    [PDF全文]
刘浪1, 王燕丽1, 杜世勇2, 侯鲁健2, 杨文1, 张文杰1, 韩斌1, 白志鹏1     
1. 中国环境科学研究院, 环境基准与风险评估国家重点实验室, 北京 100012;
2. 济南市环境保护科学研究院, 济南 250014
摘要: 利用在线单颗粒物气溶胶质谱仪(SPAMS)对2014年1月北京市典型大气重污染过程进行了连续监测,分析了具有正负离子质谱信息的颗粒物共2248225个.同时,利用ART-2a神经网络分类方法并结合Matlab统计分析,将具有质谱信息的颗粒物归为10类,分别为:矿尘类颗粒物(Dust)、元素碳颗粒物(EC)、有机碳颗粒物(OC)、元素碳和有机碳混合颗粒物(ECOC)、钠钾颗粒物(NaK)、富钾颗粒物(K)、含氮有机物(KCN)、高分子有机物(MOC,Macromolecular OC)、多环芳烃类颗粒物(PAHs)和重金属类颗粒物(Metal).结合PM2.5质量浓度数据和HYSPLIT 4.0后向轨迹模型结果,将观测时间段划分为3个典型污染过程和1个清洁过程.结果显示,重污染期间OC、MOC和PAHs为最主要的颗粒物类型.最后,本文还比对分析了污染过程和清洁期间颗粒物的混合状态,结果表明,污染过程中硫酸盐和硝酸盐较清洁期间更容易与碳质颗粒物结合.
关键词: 北京市    单颗粒物质谱    重污染过程    混合状态    
Single particle analysis during heavy air pollution episodes in January 2014 in Beijing
LIU Lang1, WANG Yanli1, DU Shiyong2, HOU Lujian2, YANG Wen1, ZHANG Wenjie1, HAN Bin1, BAI Zhipeng1     
1. State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012;
2. Jinan Academy of Environmental Sciences, Jinan 250014
Abstract: To investigate the composition and possible sources of particles in Beijing, a single particle mass spectrometer (single particle aerosol mass spectrometry, SPAMS) was deployed to measure changes of single particle species and sizes in Beijing in January 2014. A total of 2,248,225 particles with both positive and negative spectrums were collected and characterized in combination with the ART-2a neural network algorithm. Ten types of particles were classified as dust particles (Dust), elemental carbon(EC), organic carbon(OC), EC and OC combined particles(ECOC), K-Na containing particles(NaK), K-containing particles (K), organic nitrogen and potassium containing particles(KCN), macromolecular organic carbon(MOC), polycyclic aromatic hydrocarbons(PAHs), metal-containing particles(Metal). Three haze pollution episodes and one clean period were observed based on PM2.5 mass concentration and back trajectory results from Hysplit-4 models. OC, MOC and PAHs were the major components of single particles during the three haze pollution periods, which showed obviously higher abundance compared with clean periods. Results from the mixing state of secondary species with different types of aged particles showed that sulfate and nitrate were more easily mixed with carbon-containing particles during haze pollution events than in clean periods.
Key words: Beijing    single particle aerosol mass spectrometry    haze pollution    mixing state    
1 引言(Introduction)

大气颗粒物来源广泛(Prather,2009),其中包含危害人体健康的有毒化学成分(Pope et al., 2006).北京近年来大气重污染过程频发,传统的离线滤膜采样方法时间分辨率低,因此,实时在线分析单个气溶胶组分信息对于了解大气颗粒物污染特征具有一定的优势.

自20世纪90年代起,国外开始利用单颗粒气溶胶质谱分析大气颗粒物.例如,Moffet等(2008)利用单颗粒气溶胶飞行时间质谱仪(ATOFMS)研究了墨西哥北部大气颗粒物的类型和混合状态;Gard等(1998)利用ATOFMS直接观测大气颗粒物,研究沿海城市灰霾形成原因;Silva等(1999)利用ATOFMS得到了南加州8种生物质燃烧排放颗粒物的粒径和化学特征,并发现左旋葡聚糖能够作为生物质燃烧颗粒物的标识物;Healy等(2012)利用ATOFMS研究了不同粒径下元素碳颗粒物的源和混合状态,将元素碳颗粒物分成4类:生物质燃烧排放、交通排放、内部混合OC和SO42-及内部混合OC和NO3-.国内利用单颗粒质谱研究大气颗粒物起步较晚,2007年以后,国内开始了单颗粒物气溶胶质谱的研究工作.Yang等(2009)利用ATOFMS研究了上海夏季草酸的颗粒物混合状态和形成机制,发现生物质燃烧是颗粒物中草酸的主要一次源;Bi等(2011)利用单颗粒气溶胶质谱仪(SPAMS)分析了2010年4—5月间珠三角区域城市亚微米气溶胶的粒径和化学组成,认为生物质燃烧为亚微米颗粒的一个显著来源;张国华等(2013)利用SPAMS对比分析了广州地区春秋季亚微米含碳颗粒的混合状态,比较了不同季节硝酸盐的形成机制,发现秋季颗粒物中硝酸盐更容易与生物质和OC结合,春季硝酸盐的分布比较均匀;Li等(2014)利用SPAMS研究了北京市2011年春季大气颗粒物特征,将观测期间颗粒物分为富钾颗粒(39.79%)、含碳颗粒(32.7%)和重金属颗粒(19.2%),并研究了沙尘过程中大气颗粒物的粒径分布和化学组分特征.基于SPAMS仪器的特点,本文应用该仪器分析北京市2014年1月典型大气重污染过程中主要的大气单颗粒物类型、组分变化及混合特性,以期为了解北京市大气重污染特征提供基础数据支持.

2 实验和方法(Experiments and methods) 2.1 采样点信息

采样地点位于北京市朝阳区北五环外中国环境科学研究院内(东经116°24′55″,北纬40°2′29″),采样头高度距离地面约5 m.采样时间为2014年1月16—29日.同期,PM2.5小时平均质量浓度采用奥体中心监测站公布的数据(http://www.pm25.in),其他气象数据包括温度、相对湿度、风速、风向和能见度来自于网站(http://www.wundergr and .com/).

2.2 采样方法及数据分析

本文利用在线单颗粒物气溶胶质谱仪(SPAMS, 0515,广州禾信分析仪器有限公司)对大气颗粒物进行连续采样,该仪器利用空气动力学透镜将分散的大气颗粒物聚焦为一准直线,颗粒物逐个进入仪器,通过双激光测径系统和飞行时间质量分析器对每个颗粒物的粒径和化学成分进行测定,SPAMS采样流量为100 mL · min-1,可测粒径范围为200~2000 nm,其原理介绍见文献(Li et al.,2011).

正式采样前采用已知粒径(0.23、0.32、0.51、0.73、0.96、1.4和2.0 μm)的聚苯乙烯标准小球(PSL)对SPAMS粒径进行校正.流程如下:①将标准小球溶于超纯水中,混合均匀后,使用气溶胶发生器(Model 9032,TSI,USA)将产生的喷雾(含有标准小球)经硅胶扩散干燥管干燥后进入SPAMS;②测定标准小球在SPAMS两束测径激光之间的飞行时间,根据粒径与飞行时间之间的关系拟合得到粒径校正方程(图 1),粒径校正系数R2>0.99.实际采样中,单个颗粒物的飞行时间可由仪器直接测得,然后根据粒径校正方程得到颗粒物的粒径.

图1 SPAMS粒径校正曲线 Fig.1 Size calibration of SPAMS
3 结果与讨论(Results and discussion) 3.1 重污染过程的识别

图 2给出了采样期间各气象要素、PM2.5质量浓度小时均值及SPAMS数据时间序列.根据PM2.5质量浓度小时均值数据(图 2a),在观测期间出现了5次PM2.5质量浓度小时均值大于100 μg · m-3的污染过程,小时均值最高达到375 μg · m-3,初步将这5个过程定为污染过程,分别表示为P1、P2、P3、P4、P5,而相对的清洁过程对应的PM2.5质量浓度小时均值维持在50 μg · m-3以下的低浓度水平.图 2b为SPAMS数据,包括测径颗粒物个数(Size)和打击到的颗粒物数(Hit).观测期间各气象参数见图 2c~e.

图2 各气象要素、PM2.5质量浓度及SPAMS数据时间序列(a.PM2.5质量浓度;b.SPAMS测径颗粒数(Size)、打击颗粒数(Hit);c.温度与相对湿度;d.风向与风速;e.能见度) Fig.2 Time series of meteorological parameters and SPAMS data(a. PM2.5 mass concentration; b. Particle numbers of SPAMS Sized & Hit particle; c. Temperature and Relative humidity; d.Wind direction and Wind speed; e.Visibility)

为研究不同污染过程中颗粒物的来源情况,本文利用NOAA HYSPLIT 4.0模式对观测期间的气团来源进行分析,选取距离地面500 m高度处的气团,模拟气团36 h的后向轨迹(图 3),每天计算4个时次(02:00、08:00、14:00、20:00,北京时),对P1~P5进行进一步归类.

图3 采样点500 m(AGL)处气团36 h后项轨迹聚类结果 Fig.3 Daily backward(-36 h)air mass trajectories in Beijing 500 m above sea level

气团后向轨迹显示,观测期间影响北京地区的气团主要有3类,即低空东南向局地传输(35%)、中低空偏北向区域传输(40%)及高空西北向跨区域传输(25%).观测期间模拟出60条气团后向轨迹线,分析发现,前3次污染过程(P1、P2、P3)中53%的气团属于低空东南向局地传输,47%的气团属于中低空偏北向区域传输;第4次污染过程(P4)中54%的气团属于低空东南向局地传输,37%的气团属于中低空偏北向区域传输,9%的气团属于高空西北向跨区域传输;第5次污染过程(P5)中85%的气团属于低空东南向局地传输,15%的气团属于高空西北向跨区域传输.根据污染气团传输方向,结合PM2.5质量浓度结果,将观测期间5次污染过程归为3类,即①Pa:16日14时至19日16时(受第1类和第2类气团影响,P1、P2和P3);②Pb:22日16时至24日12时(受3类气团共同影响,P4);③Pc:26日20时至27日14时(受第1类和第3类气团影响,P5).为了对比分析污染过程与非污染过程颗粒物化学组分的变化情况,将20日02时至21日14时定为典型清洁天(Pd),该时间段气团来自高空西北向跨区域传输.

3.2 颗粒物分类结果

本文利用自适应共振神经网络算法(ART-2a),根据颗粒物正负离子质谱图的相似程度,将具有相似谱图特征的颗粒物归为一类.采样期间共测得颗粒物2777460个,其中,同时具有正负离子质谱信息的颗粒物共2248225个,分类结果为10类:矿尘类颗粒物(Dust)、元素碳颗粒物(EC)、有机碳颗粒物(OC)、元素碳和有机碳混合颗粒物(ECOC)、钠钾颗粒物(NaK)、富钾颗粒物(K)、含氮有机物(KCN)、高分子有机物(MOC,Macromolecular OC)、多环芳烃类颗粒物(PAHs)和重金属类颗粒物(Metal).各类颗粒物平均质谱图见图 4.

图4 各类颗粒物正负离子平均质谱图 Fig.4 Average positive and negative mass spectra for main particle classes

图 4所示,矿尘类颗粒物(Dust)质谱图中含有56Fe+40Ca+56CaO+79PO3-76SiO3-等地壳元素特征峰(图 4a).元素碳颗粒物(EC)正负离子质谱图(图 4b)中主要以碳簇离子[Cn](m/z=12、24、36、48、60、72、84等,且24C2-36C3-48C4-有较高的信号强度)为主,EC排入大气后表面会不断附着 18NH4+46NO2-62NO3-80SO3-96SO4-97HSO4-125H(NO3)2- 等二次成分,在大气中老化的过程使其逐渐变性为二次气溶胶颗粒(牟莹莹,2013),因此,该类颗粒物质谱图中有一定的46NO2-62NO3-97HSO4-离子峰.有机碳颗粒物(OC)阳离子质谱峰主要为大量的有机碳离子碎片(图 4c),如37C3H+51C4H3+63C5H3+74C6H+285C7H+等,阴离子中97HSO4-离子峰最强,反映出OC颗粒易于和大气中的硫酸盐成分结合.此外,其他二次成分如硝酸盐(m/z=-46、-62)和26CN-也有很强的峰值.元素碳和有机碳混合颗粒物(ECOC)同时具有EC和OC的特征峰,且阴离子中硝酸盐的峰值强度高于硫酸盐(图 4d),二次离子组分(26CN-46NO2-62NO3-97HSO4-)的存在表明ECOC在大气中经过了老化过程.图 4e为钠钾类颗粒物(NaK)质谱图,这类颗粒物阳离子中23Na+39K+的离子信号很明显,其他阳离子没有特别明显的峰值出现;阴离子中46NO2-62NO3-为特征的离子峰信号较强,96SO4-97HSO4-26CN-24C2-等阴离子也存在一定的信号峰,但峰值较弱.图 4f为富钾类颗粒物(K)质谱图,这类颗粒物阳离子质谱图中39K+的峰值非常强,而其他阳离子的信号很弱甚至没有,阴离子中46NO2-62NO3-的信号较强;39K+是生物质燃烧排放颗粒的最好标志物,此外,左旋葡萄糖离子碎片(45CHO2-59C2H3O2-71C3H3O-)是生物质燃烧的标志性物质,但由于其在大气中很容易被氧化而消失(Pratt et al.,2011),因而只有很少的颗粒物能够检测出包含左旋葡萄糖离子碎片.图 4g为含氮有机物(KCN)质谱图,其标志性离子峰为26CN-58[C2H5NHCH2]+59[N(CH3)3]+(张国华,2013;Moffet et al., 2008).图 4h为重金属类颗粒物(Metal),该类颗粒物有很强的206,207,208Pb+信号,存在较弱的51V+67VO+63Cu+等重金属元素信号;阴离子中有较强的硝酸盐离子信号.高分子有机物(MOC)正负离子质谱图见图 4i,MOC正谱图中存在大量质荷比大于100的大分子离子碎片峰,且这些大分子离子碎片峰质荷比相差14(CH2),阴离子中二次离子(46NO2-62NO3-97HSO4-)的信号较强.多环芳烃类颗粒物(PAHs)正谱图中质荷比主要集中在189、202和228,并伴有大量有机物碎片峰(图 4j).

3.3 各类颗粒物数量随粒径和时间变化特征

各类颗粒物数量及其随粒径和时间的变化特征见图 5. SPAMS打击到的颗粒物个数结果列于表 1.观测期间,打击到的颗粒物粒径主要分布在0.4~0.7 μm之间,且在0.5~0.6 μm处打击到的颗粒物总数最多(1243912个).整个粒径段OC、PAHs、MOC、ECOC及K类颗粒物个数最多,分别为611437、388622、479334、244281和217765个.对于不同类型的颗粒物,SPAMS对其检测的主要粒径段也不同.例如,OC颗粒粒径主要集中在0.45~0.65 μm(407413个),PAHs颗粒物粒径主要分布在0.45~0.7 μm段(268643个),MOC颗粒物粒径主要分布在0.45~0.65 μm(355722个).

图5 各类颗粒物数量随粒径和时间变化(a.各类颗粒物的数量随粒径变化;b.各类颗粒物的数量随时间变化) Fig.5 Temporal variation of number of size-resolved particles(a. variation of sizes; b. temporal variation)

表1 各类颗粒物粒径分布统计结果 Table 1 Statistics of different types of size-resolved particle

图 5b为各类颗粒物数量随时间的变化趋势,在重污染期间,OC、PAHs和MOC颗粒个数有明显的增加,即该3类颗粒物与重污染过程有较好的时间一致性,表明OC、PAHs和MOC颗粒物是影响北京市重污染过程最主要的颗粒物.

3.3 不同污染过程颗粒物组分变化特征

图 6是3个典型污染过程(Pa、Pb、Pc)和典型清洁期间(Pd)不同类型颗粒物数量变化,以及各类颗粒物所占总颗粒物数的百分比变化.

图6 3个典型污染过程及清洁过程各类颗粒物数量 Fig.6 Number variation of different particle types during 3 typical air pollution events and clean period

经统计得出各污染过程中各类颗粒物个数及数量百分比,具体见表 2. 3个污染过程中,颗粒物总数均比清洁过程有非常明显的增加,分别增加379410、208699和406701个.不同污染过程中,颗粒物主要类型也不尽相同.Pa期间,OC、PAHs、ECOC和K类颗粒占较大比例,其数量百分比分别为22.5%、18%、15%和14%;Pb期间,OC(27.2%)和PAHs(23.4%)数量百分比较大;Pc期间则以OC(31.3%)和MOC(28.1%)颗粒为主,且MOC颗粒物在Pc期间有明显增加.结合污染过程的气象条件发现,观测期间,低风速静稳天气条件多伴随高污染过程,一方面是由于静稳天气不利于污染物的清除和输送,继而污染物累积最终导致重污染事件的发生;另一方面,静稳天气可促进大气中复杂的物理化学过程从而生成二次有机气溶胶(SOA)(Guo et al.,2014),气-粒转化作用加强.

表2 不同污染过程各类颗粒物数量和所占数量百分比统计结果 Table 2 Statistical number and percentage results of different particle types during different pollution events

本文观测结果显示,PAHs颗粒物数量百分比较大,尤其是在Pb中,占23.4%,这与其他研究结果有所差异,其他研究普遍认为大气颗粒物中的PAHs含量较低(金银龙等,2011);于国光等(2008)研究了北京市大气气溶胶中的PAHs,结果显示,PAHs质量浓度为0.10~36.71 ng · m-3.目前,对于大气PAHs的分析手段主要为滤膜采样分析,样品在保存和运输的过程中很容易引起PAHs的损失.本文所采用的仪器为在线分析仪器,这可能是导致本文PAHs占比结果偏大的原因之一.

3.4 颗粒物混合状态

混合状态反映了不同颗粒物上各化学组分的含量情况.颗粒物中二次颗粒物组分的存在能够在一定程度上反映出颗粒物在大气中经历的老化过程,因此,本文选择将二次颗粒物组分作为衡量颗粒物混合状态的依据.在单颗粒物研究中,二次颗粒物组分通常选择铵盐、硫酸盐、硝酸盐和含氮有机物产生的质谱峰(m/z=±18、-97、-46/-62和-26),以及二次有机气溶胶43[C2H3O]+产生的质谱特征峰(张国华,2013).图 7为污染过程和清洁期间颗粒物的混合状态.

图7 污染过程(a、b、c)和清洁过程(d)各类颗粒物与二次气溶胶的标志质谱峰的混合状态(不同颜色表示二次气溶胶的标志质谱峰在各颗粒物类型中的比重) Fig.7 Mixing state of secondary markers on the various particle classes(Color represents number fraction of particle classes(y-axis)that contain secondary makers)

对比分析清洁过程和污染过程各类颗粒物中二次成分的混合情况,对于PAHs和Dust颗粒物,二次成分的含量百分比在清洁过程和污染过程并没有明显的变化,但其他类型颗粒物中二次成分的含量均表现出一定的差异.例如,对于OC、NaK和K颗粒物,污染过程中硫酸盐的含量百分比相比清洁过程的高,而铵根的含量百分比较清洁过程却有下降的趋势,表明污染期间OC、NaK和K颗粒与硫酸盐的结合能力较强,在清洁过程中OC、NaK和K颗粒与铵盐的结合能力虽会增强,但总的混合状态还是以硫酸盐和硝酸盐为主.对于Metal颗粒物,污染过程其与硫酸盐和二次有机碳结合能力较清洁过程增强,而与铵盐和硝酸盐的结合能力并未有明显变化,观测期间,硝酸盐在重金属颗粒中均占有最大的比例.对于MOC颗粒物,硝酸盐、硫酸盐和二次有机碳均占很大的比例,污染过程中,硫酸盐与MOC颗粒物的结合能力有所增加,反映出重污染期间MOC更易于硫酸盐混合.对于EC和ECOC颗粒,污染过程中硫酸盐和硝酸盐比重明显增加.对于KCN颗粒物,污染过程中硫酸盐和硝酸盐的混合程度较清洁过程有明显增加,而铵盐和二次有机碳变化不大,表明污染过程中KCN更易与硫酸盐和硝酸盐结合.

观测期间,硝酸盐和硫酸盐在各类颗粒物中均占有很大的比重,其次为二次有机碳成分,而铵盐的比重相对较小.污染过程颗粒物中的硫酸盐和硝酸盐比重较清洁过程会更大,尤其是在碳质颗粒物中(EC、ECOC、MOC和OC),表明在污染过程中硫酸盐和硝酸盐更容易与含碳颗粒物结合.

4 结论(Conclusions)

1)本文利用SPAMS对北京市2014年1月大气颗粒物进行连续观测,并利用ART-2a分类方法将颗粒物分为10类,发现颗粒物粒径主要分布在0.4~0.7 μm之间.

2)研究选取了观测时段的3个典型污染过程和一个典型清洁过程,分析了不同污染过程和清洁过程颗粒物组分的变化特征,结果显示,重污染期间颗粒物类型主要为OC、MOC和PAHs.

3)本文比对分析了典型污染过程和清洁过程颗粒物的混合状态,结果显示,观测期间硝酸盐和硫酸盐在各类颗粒物中均占有很大的比重,尤其在碳质颗粒物中,而且污染期间含碳颗粒物中的硫酸盐和硝酸盐较清洁期间比重大,表明污染过程硫酸盐和硝酸盐更容易与含碳颗粒物结合.

参考文献
[1] Bi X,Zhang G,Li L,et al.2011.Mixing state of biomass burning particles by single particle aerosol mass spectrometer in the urban area of PRD,China[J].Atmospheric Environment,45(20):3447-3453
[2] Gard E E,Kleeman M J,Gross D S,et al.1998.Direct observation of heterogeneous chemistry in the atmosphere[J].Science,279(5354):1184-1187
[3] Guo S,Hu M,Zamora M L,et al.2014.Elucidating severe urban haze formation in China[J].Proceedings of the National Academy of Sciences:111(49):17373-17378
[4] Healy R M,Sciare J,Poulain L,et al.2012.Sources and mixing state of size-resolved elemental carbon particles in a European megacity:Paris[J].Atmospheric Chemistry and Physics,12(4):1681-1700
[5] 金银龙,李永红,常君瑞,等.2011.我国五城市大气多环芳烃污染水平及健康风险评价[J].环境与健康杂志,(9):758-761
[6] Li L,Huang Z,Dong J,et al.2011.Real time bipolar time-of-flight mass spectrometer for analyzing single aerosol particles[J].International Journal of Mass Spectrometry,303(2):118-124
[7] Li L,Li M,Huang Z,et al.2014.Ambient particle characterization by single particle aerosol mass spectrometry in an urban area of Beijing[J].Atmospheric Environment,94:323-331
[8] Moffet R C,Foy B D,Molina L A,et al.2008.Measurement of ambient aerosols in northern Mexico City by single particle mass spectrometry[J].Atmospheric Chemistry and Physics,8(16):4499-4516
[9] 牟莹莹.2013.利用单颗粒气溶胶飞行时间质谱研究气溶胶的老化与混合状态及定量计算[D].上海:华东理工大学
[10] Pope Iii C A,Dockery D W.2006.Health effects of fine particulate air pollution:lines that connect[J].Journal of the Air & Waste Management Association,56(6):709-742
[11] Prather K A.2009.Our current understanding of the impact of aerosols on climate change[J].Chem Sus Chem,2(5):377-379
[12] Pratt K A,Murphy S M,Subramanian R,et al.2011.Flight-based chemical characterization of biomass burning aerosols within two prescribed burn smoke plumes[J].Atmospheric Chemistry and Physics,11(24):12549-12565
[13] Silva P J,Liu D,Noble C A,et al.1999.Size and chemical characterization of individual particles resulting from biomass burning of local Southern California species[J].Environmental Science & Technology,33(18):3068-3076
[14] Yang F,Chen H,Wang X,et al.2009.Single particle mass spectrometry of oxalic acid in ambient aerosols in Shanghai:Mixing state and formation mechanism[J].Atmospheric Environment,43(25):3876-3882
[15] 于国光,王铁冠,吴大鹏,等.2008.北京市大气气溶胶中多环芳烃的研究[J].中国矿业大学学报,(1):72-78
[16] 张国华.2013.城市大气气溶胶混合状态的初步研究[D].广州:中国科学院广州地球化学研究所
[17] Zhang G,Bi X,Li L,et al.2013.Mixing state of individual submicron carbon-containing particles during spring and fall seasons in urban Guangzhou,China:a case study[J].Atmospheric Chemistry and Physics,13(9):4723-4735