环境科学学报  2017, Vol. 37 Issue (1): 95-103
大气气溶胶碳质组分在线分析仪的研制和应用    [PDF全文]
薛瑞, 曾立民 , 吴宇声, 张远航    
北京大学环境科学与工程学院, 环境模拟与污染控制国家重点实验室, 北京 100871
摘要: 大气气溶胶碳质组分主要包括有机碳(Organic Carbon,OC)和元素碳(Elemental Carbon,EC),是大气细粒子的重要组成部分,对大气环境及人体健康具有非常重要的意义.针对目前商品化碳质气溶胶在线分析仪存在的一些问题,研制出大气气溶胶碳质组分在线分析仪(PKU-OCEC),实现了对OC和EC的小时浓度测量.仪器同时采用热分解-光学透射(Thermo Optical Transmission,TOT)法和热分解-光学反射(Thermo Optical Reflection,TOR)法分离OC和EC,检测限为0.29 μg·m-3(以C计)(采样时间40 min,采样流量8 L·min-1),测量范围最大值为211.50 μg·m-3(以C计),可以满足不同环境的大气观测要求.于2016年1月23-31日采用该仪器对北京怀柔大气进行监测,并与商品化在线碳质气溶胶分析仪(RT-4 OCEC,Sunset Lab Inc.)和多角度吸收光度计(5012 MAAP,Thermo Scientific)进行比对.结果显示,PKU-OCEC测得的TC(Total Carbon,TC=OC+EC)和OC浓度均与Sunset-OCEC对应结果显著相关(TC:r=0.986,OC:r=0.984;p<0.01),说明两台仪器测量结果具有很好的一致性;Sunset-OCEC测得的EC数据有39.39%低于0.013 μg·m-3(以C计),剔除该部分异常值后,其测得的EC数据与PKU-OCEC数据的线性拟合方程由CEC_Sunset=0.9278×CEC_PKU-0.0704(TOT法)和CEC_Sunset=0.9476×CEC_PKU-0.2665(TOR法)变为CEC_Sunset=0.6810×CEC_PKU+1.2842(TOT法)和 CEC_Sunset=0.7035×CEC_PKU+1.1179(TOR法),新的方程产生较大正截距,说明Sunset-OCEC的检测限偏高;PKU-OCEC测得的EC数据与MAAP测得的BC数据显著相关,其线性拟合度(R2=0.8868(TOT法)、R2=0.8614(TOR法))大于Sunset-OCEC与MAAP的拟合度(R2=0.8047);PKU-OCEC分析仪TOT法和TOR法测得的OC浓度相当(COC_TOR=0.9963×COC_TOTR2=0.9959).综合结果表明,PKU-OCEC分析仪测量数据可靠,运行稳定且操作简便,可以用于各种环境条件下的大气监测.
关键词: 气溶胶碳质组分在线分析仪     有机碳元素碳分析仪     多角度吸收光度计(MAAP)     监测技术    
Development and application of an on-line carbonaceous aerosol analyzer
XUE Rui, ZENG Limin , WU Yusheng, ZHANG Yuanhang    
State Joint Key Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing 100871
Received 15 Mar 2016; received in revised from 28 Apr 2016; accepted 28 Apr 2016
the Demonstration Project on Development and Application of on-line Monitoring Equipment for Chemical Composition of Atmospheric Fine Particulate Matter of Ministry of Environment Protection of China (No.2013YQ060569),the Strategic Priority Research Program of Chinese Academy of Sciences(No.12700-452-10161-029) and the National High Technology Research and Development Program(No.SS2014AA063007)
Biography: XUE Rui(1990—),female,E-mail:sherry13@yeah.net
*Corresponding author: E-mail:lmzeng@pku.edu.cn
Abstract: Carbonaceous aerosols, i.e. organic carbon (OC) and elemental carbon (EC), are major composition of atmospheric fine particles, and have great impact on atmospheric environment and human health. In order to improve the commercial on-line carbonaceous aerosol analyzer, an on-line carbonaceous aerosol analyzer (PKU-OCEC) was developed, which is capable in measuring the hourly mass concentration of OC and EC in the atmospheric particles. The instrument combined both Thermo Optical Transmission (TOT) method and Thermo Optical Reflection (TOR) method to separate OC and EC, with the detection limit of 0.29 μg C·m-3(with a sample flow of 8 L·min-1 and a sample time of 40 min). The maximum measurable concentration is 211.50 μg C·m-3, which satisfies the requirement of the observations in various atmospheric environments. The instrument was used to measure the ambient particles at a rural site at Huairou in the northeast of Beijing during 23-31 January 2016. Inter-comparison was conducted between newly developed instrument, the commercial on-line carbonaceous aerosol analyzer (RT-4 OCEC, Sunset Lab Inc.) and Multi-Angle Absorption Photometer (5012 MAAP, Thermo Scientific). The concentration of Total Carbon(TC,TC=OC+EC) and OC were well correlated between PKU-OCEC and Sunset-OCEC (with the correlation coefficients of TC:r=0.986, OC:r=0.984; p<0.01), exhibiting good consistency. 39.39% of the EC data measured by Sunset-OCEC were below 0.013 μgC·m-3. When removing the outliers, the liner regression equations of the two sets of the data from PKU-OCEC and Sunset-OCEC changed from CEC_Sunset=0.9278×CEC_PKU-0.0704(TOT method) and CEC_Sunset=0.9476×CEC_PKU-0.2665(TOR method) to CEC_Sunset=0.6810×CEC_PKU+1.2842(TOT method) and CEC_Sunset=0.7035×CEC_PKU+1.1179(TOR method). The larger positive intercepts of the corrected liner regression equations indicated higher detection limit of Sunset-OCEC. The EC concentrations measured by PKU-OCEC correlated better with those by MAAP, with the correlation coefficients of R2=0.8868(TOT method) and R2=0.8614(TOR method), than that between Sunset-OCEC and MAAP(R2=0.8047). OC concentrations measured by TOT and TOR method with PKU-OCEC were consistent (COC_TOR=0.9963×COC_TOT,R2=0.9959). Our results indicate that PKU-OCEC performs well, and is capable in producing reliable measurements. It can therefore be used to measure ambient OC and EC concentrations under various atmospheric environments.
Key words: on-line carbonaceous aerosol analyzer     OCEC analyzer     MAAP     monitoring techniques    
1 引言(Introduction)

大气气溶胶中的碳质组分通常占大气细粒子质量浓度的10%~70%,是大气细粒子的重要组成成分(Turpin et al.,1991Steinfeld,1998Polidori et al.,2006).其可分为3大类:有机碳(Organic Carbon,OC)、元素碳(Elemental Carbon,EC)和碳酸碳(Carbonate Carbon,CC)(Chow et al.,2002Cao et al.,2003).通常来说,OC是指一种含有上百种有机化合物(如脂肪族、芳香族、酸类等)的混合体,主要来源于一次燃烧过程的排放、生物排放,以及气态有机污染物参与光化学反应和气粒转化过程的排放(Offenberg et al.,2002);EC通常是指化石燃料或生物质等含碳物质经过不完全燃烧后直接排放的一种结晶度很低的无定形碳(Tanner et al.,1982),其具有很强的吸光性,能吸收可见光、近红外和近紫外光波段(Chen et al.,2004);CC主要存在于土壤和煤矿飞尘的粗粒子中,其质量浓度远小于OC和EC,故一般被忽略(Gray et al.,1986).

大气气溶胶中的碳质组分能够对全球气候、大气能见度、人体健康等方面产生影响:碳质气溶胶可以作为凝结核改变大气中云滴的浓度和寿命,影响地球辐射平衡(Matsumoto et al.,1997);EC可以吸收从红外到紫外全波段的光,从而增加正辐射强迫,显著降低区域大气的能见度(Bond et al.,2006);碳质气溶胶还容易通过呼吸作用进入人体肺部,破坏其结构和功能,引发慢性呼吸道疾病甚至癌变(Nemmar et al.,2002).

在大气气溶胶碳质组分的研究中,准确测定OC、EC的浓度,对研究大气化学反应和污染物的源解析具有重要意义,成为当今环境监测领域的热点.然而,OC、EC并没有绝对的定义,其划分在很大程度上依赖于相应的检测方法(Chow et al.,2001;2004;2005),目前主要有热分解法、光学法和热分解-光学分析法.其中,热分解法主要依据OC和EC热学性质的不同对两者进行区分,但无法解决OC高温碳化导致分割点不准确的问题(Huntzicker et al.,1982).光学法主要利用光学原理测定吸光性碳质气溶胶即EC(通常把光学法测量的EC称作BC)的含量(郇宁等,2005),如Thermo公司生产的5012型多角度吸收光度计(Multi Angel Absorption Photometer,MAAP),通过测量采集到玻璃纤维膜上的BC对光的吸收和散射,经校正系数折算后得到气溶胶中BC的含量(Hitzenberger et al.,1999);然而,光学法无法得到OC的结果.

热学法中OC的碳化产物具有吸光性,在热学法的基础上使用光学方法对这部分裂解碳进行校正,这种热学分解+光学校正的方法即为目前最为成熟的热分解-光学分析法.它基于以下两个假设(吴琳等,2010):第一,由OC碳化生成的EC(Pyrolytically generated EC,PEC)先于颗粒物中原有的EC(Native EC,NEC)氧化分解;第二,PEC在分析过程中的吸光系数始终等同于NEC.然而研究表明,PEC的吸光系数在分析过程中会发生一定的变化,且在大多数情况下和NEC的吸光系数并不相等;有氧加热时,PEC和NEC的氧化分解过程并非各自独立,而是交错进行的.因此,热分解-光学校正法仍具有一定误差.尽管如此,其对OC碳化做了较好的修正,目前应用最为广泛.热分解-光学分析法中最为成熟是NIOSH/TOT(The National Institute for Occupational Safety and Health/Thermo Optical Transmission)法和IMPROVE/TOR(Interagency Monitoring of Protected Visual Environments/Thermo Optical Reflection)法,典型的使用这两种方法的仪器分别是美国Sunset实验室的离线和在线OCEC分析仪(Birch et al.,1996)及沙漠所的离线OCEC分析仪(Chow et al.,1993).在线分析仪克服了传统离线分析仪时间分辨率低、样品不易储存、易受人为干扰等缺陷,可以研究时间尺度较短的大气化学过程,具有更大的应用市场.

在实际应用中发现,目前商品化的在线OCEC分析仪价格昂贵、维护成本高,在结构和操作方面也存在一些不合理的地方(刘巧铃,2010):仪器操作软件和数据处理软件独立,不能实现对以往数据的快速查看;通信串口、升温程序和载气流量等参数的更改只能通过调用“.par”类型的文件完成,不仅操作复杂,还容易因字符位置输入错误导致软件无法正常工作;仪器仅使用TOT法分离OC和EC,而很多文献采用的是TOR法,单一的分离方法不利于实验室之间的结果比对;甲烷内标气在仪器分析过程中全程通入,造成了极大浪费.基于以上问题,本实验室研制出一套大气气溶胶碳质组分在线分析仪(PKU-OCEC),同时采用TOT和TOR法分离OC和EC,以期为进一步了解我国气溶胶的污染现状,治理大气污染,改善空气质量提供支持.

2 仪器构造与原理(Structure and principle of the instrument) 2.1 系统构成与配件

气溶胶碳质组分在线分析仪由自动控制系统、气路系统、加热系统、激光系统等组成,其中,自动控制系统包括控制电路、单片机和工控电脑(PPC-070C,灵江工控计算机技术有限公司);气路系统包括大流量质量流量控制器(Mass Flow Controller,MFC;MFC-20LAIR,天津吉思特仪器仪表有限公司)、小流量MFC(MFC-100HE,天津吉思特仪器仪表有限公司)、大流量二通常闭电磁阀(UF06,Uni-D公司)、小流量二通常闭电磁阀、二位六通电磁阀、球阀(KLD20S,KLD公司)、非色散红外(Non-Dispersive Infra-Red Spectroscopy,NDIR)CO2气体检测器(Model DX6220,RMT公司)、压力传感器(Honeywell公司)、PM2.5旋风切割头(URG公司)、挥发性有机物(Volatile Organic Compounds,VOCs)去除管、采样泵(HP-40V,AIRTEC公司)等;激光系统包括激光发射器(3.3V DC,200 mW,中心发射波长660 nm)、窄带滤光片(波长660 nm,带宽8 nm)、硅光电池激光接收器(BPX65,OSRAM公司)等;加热系统包括数字温控仪(E5CC,欧姆龙公司)、K型热电偶、加热炉丝(1200 W,钢花牌)等.此外,仪器还包含其他部件如风机(MQ12032,闽泉电机)、散热风扇、固态继电器、开关电源等;仪器使用氦气(He体积分数99.999%以上)、氦氧混合气(He/O2混合气,O2体积分数为10%)和氦甲烷混合气(He/CH4混合气,CH4体积分数为5%)高压钢瓶气,减压阀调节钢瓶出口气体压力至0.2 MPa.

2.2 仪器工作原理

仪器同时采用TOT法和TOR法分离OC和EC,首先,在8 L·min-1的流量下采集大气样品,颗粒物样品被收集到石英纤维膜上,接着通入He气,并在无氧环境下逐级加热石英膜,样品中的OC挥发出来(一部分OC高温碳化生成EC),进入氧化炉管并被MnO2催化氧化为CO2;再向载气中通入He/O2混合气,并继续对石英膜逐级升温,样品中的EC被氧化分解并逸出,同样进入氧化炉管最终转化为CO2;OC和EC转化生成的CO2均进入NDIR检测器被检测.整个过程中,激光发射器始终发射一束波长660 nm的激光,经石英膜透射和反射后到达透射激光接收器和反射激光接收器;无氧加热石英膜,一部分OC发生碳化,使得透射和反射激光减弱;随着He切换成He/O2,EC被氧化分解,透射和反射激光又逐渐增强,当恢复到无氧加热开始前的光强时即认为到达OC、EC的分割点,即:该时刻之前热分解出的碳质组分是OC,之后为EC.升温程序结束后,向定量环中通入并充满He/CH4混合气,混合气被MnO2催化氧化为CO2并被NDIR检测器定量.通过比较OC、EC和He/CH4混合气的CO2响应积分面积,得到样品中OC、EC的含量,再根据采样时间和采样流量即可算出大气中OC、EC的浓度.

2.3 仪器结构组成

图 1所示,MFC1~MFC4分别控制He/CH4混合气、He/O2混合气、He主气和He吹扫气的流量;二通电磁阀1和2共同控制He/O2气的通断,两者始终处于相反的工作状态;二通电磁阀3控制气体进入NDIR检测器;球阀、二通电磁阀4、MFC5、限流阀和采样泵一起控制大气样品的采集,其中,限流阀维持采样流量稳定,减小采样时气流变化对石英纤维膜位置的影响;二位六通电磁阀装有定量环,可以储存固定体积的He/CH4气,切换六通电磁阀可以将定量环中的He/CH4气注入系统;压力传感器监测整个分析过程中腔体的压力,以判断系统是否出现漏气,保证仪器稳定运行.解析-氧化炉管结构如图 2所示,包括解析炉主管、解析炉副管、氧化炉管和支管;主管和副管连接处内部烧有一石英环,用来放置石英纤维膜;主管内套有一进样管,进样管和石英环一起将石英纤维膜固定;解析炉副管内部烧有一密封石英片;主管和副管连接处两侧分别有一氧化炉管和一支管;氧化炉管的中间部分内部填充分析纯级MnO2催化剂;石英环及氧化炉管中间位置对应的炉管外部均缠绕有加热炉丝;K型热电偶1伸入支管内部测量石英纤维膜温度;K型热电偶2贴放在氧化炉管外部,测量催化剂温度,两K型热电偶分别与两台温控仪相连;为了保证石英膜实际温度最大程度接近升温程序,仪器采用比例-微分-积分结合脉冲宽度调制方法对加热炉丝进行控制.

图 1 仪器结构图 Fig. 1 Structure of the instrument

图 2 解析-氧化炉管结构图 Fig. 2 Structure of the thermal desorption and oxidation furnace tube

激光校正系统如图 2所示,包括反射激光校正系统和透射激光校正系统,接头1内装有一激光发射器,其发出的激光穿过石英纤维膜、接头2内的滤光片后进入透射激光接收器;与此同时,另一部分激光被石英纤维膜反射,穿过接头1内的滤光片后进入反射激光接收器.

仪器通过工控电脑、单片机和控制电路三者结合来实现自动化运行,控制电路为硬件提供工作电压和DO、AO控制信号,同时采集和放大相应的检测信号供单片机进行A/D 转换,并为工控电脑读取和识别.利用计算机语言编写仪器的自动化控制软件,软件嵌套有状态显示窗口、参数设定窗口、标准曲线窗口、数据记录窗口、手动积分窗口等,界面友好,操作简单,不仅能够实时显示仪器的运行状态,自主设定采样时间、流量、升温程序等分析条件,还能输入蔗糖多点标定结果并生成标准曲线,自动计算并保存样品OC和EC的浓度信息.

2.4 仪器工作流程

仪器工作流程包括采样过程、吹扫过程、OC分析过程、EC分析过程、甲烷定量过程、甲烷标定过程、冷却过程和待机过程,具体如下.

采样过程:球阀、二通电磁阀4、MFC5和采样泵开启,空气样品通过PM2.5旋风切割头、VOCs去除管、球阀进入解析-氧化炉管,颗粒物累积在石英膜上,剩余气体通过支管被采样泵排出;解析炉维持在20 ℃,并开启风机对其降温;氧化炉维持在500 ℃.

吹扫过程:MFC2、MFC3、MFC4、二通电磁阀2和二通电磁阀3开启,其余阀体和风机关闭;He气一路通过MFC4、支管进入解析-氧化炉管,另一路通过MFC2和MFC3,经六通阀进入解析-氧化炉管;两路He气一起将残留在管内的杂质气体吹出,保持系统处于无氧环境;氧化炉开始升温并维持在850 ℃.

OC分析过程:各阀体状态维持不变;解析炉开始程序升温,OC在无氧环境中被逐步热解析出来,进入氧化炉管氧化为CO2,再进入NDIR检测器被定量;一部分OC在高温下碳化为EC.

EC分析过程:二通电磁阀开启,二通电磁阀2关闭,He/O2混合气通过MFC2进入解析-氧化炉管;解析炉继续程序升温,EC氧化并逸出,最终转化成CO2并被NDIR检测器定量.

甲烷定量过程:MFC1打开,He/CH4混合气进入并充满定量环,该过程维持45 s,多余He/CH4混合气从排空端排出.

甲烷标定过程:MFC1关闭,He/CH4混合气停止通入定量环;六通电磁阀切换到下一位,定量环中的He/CH4气被吹入解析-氧化炉管并被氧化为CO2,再进入NDIR检测器被定量;解析炉停止加热,风机开启并对解析炉降温.

冷却过程:六通电磁阀恢复原始位置,风机继续对解析炉降温,当温度降至70 ℃时,分析阶段结束,软件计算并保存数据;仪器进入待机过程,等待下一次采样分析.

3 实验与方法(Experiments and methods) 3.1 条件实验

为了保证仪器的监测效果,设计一系列条件实验,对仪器各方面性能进行测试,确定其各个实验参数,保证仪器工作的稳定性.

3.1.1 MFC稳定性测试

仪器测得的OC、EC浓度由NDIR检测器测得的CO2浓度响应-时间积分面积得出,因此,气路系统流量的稳定性对测定结果至关重要.仪器使用MFC对气体流量进行控制,为了测试MFC的准确性,分别用安捷伦气体流量计(Electronic Flowmeter 5067-0023,Agilent公司)和皂膜流量计(Gilibrator-2,SENSIDYNE公司)对仪器使用的小流量MFC和大流量MFC进行标定,结果表明:MFC和流量计测量结果显著相关(r>0.999,p<0.01);随机选取3个单点流量对上述各MFC的标定曲线进行验证,单点偏差均在0.5%以内,证明仪器使用的MFC具有稳定准确的流量控制能力.

3.1.2 甲烷内标气稳定性测试

在分析的最后阶段,系统会向定量环中注入一定量的He/CH4内标气,这种引入内标物的方式可以减少外界干扰或者仪器运行不稳定导致的系统偏差,还可以通过内标气体的积分面积反映仪器运行的稳定性.本实验统计了184个样品He/CH4内标气对应的CO2浓度响应积分面积,如图 3所示,积分面积的平均值为31033.96,相对标准偏差为4.20%(小于5%),满足实验要求,说明仪器内标气及氧化系统稳定.

图 3 甲烷内标气稳定性实验 Fig. 3 Stability test of CO2 integral area corresponding to the internal standard gas of methane
3.1.3 仪器的标定曲线

分析开始前,需使用已知浓度的蔗糖溶液対仪器进行标定,得到蔗糖含碳量和CO2响应系数(响应系数为蔗糖氧化产物CO2的浓度积分面积与甲烷内标气氧化产物CO2的浓度积分面积的比值)的工作曲线.用微量注射器分别移取0.6、2、4、6、8、10 μL的蔗糖溶液(浓度为4.230μg·μL-1,以C计),对应大气碳质气溶胶浓度为7.93、26.44、52.88、79.31、105.75、132.19 μg·m-3(以C计)(采样流量8 L·min-1,采样时间40 min),该浓度范围覆盖正常环境中碳质气溶胶浓度范围(郇宁等,2005);将溶液滴加至已高温灼烧过(去除OC和EC)的石英膜上并进行分析,每个标准浓度点至少进行3次平行实验,结果如图 4所示,标定曲线显示出很好的线性拟合度(R2=0.9990),进一步计算得甲烷定量环含量为24.10 μg(以C计).

图 4 仪器标定曲线 Fig. 4 Calibration curve of the instrument
3.1.4 仪器的检测限和测量范围

用0.6 μL高纯水(Milli-Q,MILOPORE公司)连续测定6次,得到仪器最低检测限(3δ,其中,δ为空白样的标准偏差)为0.093 μg(以C计),相当于环境大气碳质气溶胶浓度0.29 μg·m-3(以C计)(采样流量8 L·min-1,采样时间40 min);增加3.1.3 节中蔗糖溶液的用量,得到满足标定曲线线性拟合度R2≥0.9900的仪器最大测量值67.68 μg(以C计),相当于环境大气碳质气溶胶浓度211.50 μg·m-3(以C计),该浓度满足大气环境样品的分析要求.

3.2 仪器比对实验

为了检验仪器在外场观测中的稳定性和数据可靠性,将PKU-OCEC与商品化仪器进行比对,比对地点设在中国科学院大学雁栖湖校区的测试楼楼顶(116.6870°E,40.4154°N).观测平台濒临京加路,该路段为101国道的一部分,平台向东1.6 km有怀柔北火车站,除此之外附近无明显污染源.本实验使用PKU-OCEC分析仪与Sunset-OCEC分析仪、MAAP进行比对;PKU-OCEC与Sunset-OCEC并排放置,MAAP与两者相距不超过50 m;PKU-OCEC和Sunset-OCEC均使用8 L·min-1采样流量的PM2.5切割头,切割头距地面约20 m,距楼顶约2 m.

比对时间为2016年1月23—31日,其中,27日0:00—14:00 Sunset-OCEC出现仪器故障导致数据缺失,27日10:00—14:00 PKU-OCEC因电路跳闸导致数据缺失,其余时间3台仪器均处于连续工作的稳定状态.PKU-OCEC与Sunset-OCEC均采用整点周期采样,采样时间40 min,采样膜选用双层Pall石英纤维膜,并选择Quartz升温程序(表 1);使用前用蔗糖标准溶液对两台仪器进行多点标定,标定结果与蔗糖实际浓度偏差均小于5%;MAAP时间分辨率为5 min,采样膜选用Whatman玻璃纤维膜.

表 1 quartz升温程序 Table 1 Heating program of quartz
4 结果与讨论(Results and discussion)

观测期间包含一个较为完整的大气污染过程,具有较好的代表性,1月28日出现重污染,该日1:00 TEOM测得观测期间最大PM2.5小时平均浓度107.0 μg·m-3,随后污染物浓度下降.图 5为PKU-OCEC与Sunset-OCEC测得的TC浓度的时间序列,可以看出,两者变化趋势非常一致,28日1:00两仪器均测得观测期间TC最大浓度.进一步统计Sunset-OCEC测得的EC结果,发现其37.10%的浓度数据为0 μg·m-3(以C计,下同),2.29%的浓度数据大于0 μg·m-3且≤0.013 μg·m-3,60.71%的浓度数据≥1.21 μg·m-3,没有数据落在0.013~1.21 μg·m-3范围内;由于0.013 μg·m-3的浓度很低,可近似为0,且在0.013~1.21 μg·m-3范围内没有数据,可以认为Sunset-OCEC测得的浓度低于0.013 μg·m-3的EC数据为异常数据.将这部分异常数据剔除后,对两台仪器测得的TC浓度进行线性拟合,得到方程CTC_Sunset=0.8851×CTC_PKU(R2=0.9536);保留直线截距,得到方程CTC_Sunset=0.94×CTC_PKU-1.326(R2=0.9577,r=0.986,p<0.01),可以看出,两台仪器的TC测量结果显著相关且线性拟合度较高,其测量精度基本一致,满足比对测量的要求;线性拟合直线截距为-1.326,即当PKU-OCEC测得的TC浓度低于1.41 μg·m-3时,Sunset-OCEC对应的线性拟合结果为0(负值取0),说明其检测限可能偏高.

图 5 观测期间TC比对时间序列 Fig. 5 ime series of TC concentration during the observation period

对PKU-OCEC分析仪和Sunset-OCEC分析仪测得的OC浓度进行线性拟合,得到方程COC_Sunset=0.9480×COC_PKU-1.4685(R2=0.9691,r=0.984,p<0.01)(TOT法)和COC_Sunset=0.9271×COC_PKU -1.0679(R2=0.9687,r=0.984,p<0.01)(TOR法),可以看出,PKU-OCEC两种方法测得的OC结果均与Sunset-OCEC显著相关且有较好的线性拟合度;进一步对PKU-OCEC两种方法测得的OC浓度进行线性拟合,得到方程COC_TOR=0.9963×COC_TOT(R2=0.9959,r=0.998,p<0.01),说明两种方法测得的结果显著相关且无明显差异,即分析过程中样品膜表面和厚度方向的黑度变化较为一致.

对PKU-OCEC分析仪和Sunset-OCEC分析仪测得的EC浓度进行线性拟合(图 6),得到方程CEC_Sunset= 0.9278×CEC_PKU-0.0704(R2=0.8263)(TOT法)和CEC_Sunset= 0.9476×CEC_PKU-0.2665(R2=0.8237)(TOR法),两台仪器EC测量结果的线性拟合度略低于TC和OC;如图 6所示,大量数据点落在X轴附近(已用圆圈标注),即Sunset-OCEC测得的EC浓度低于0.013 μg·m-3(以C计)的异常数据;将该部分异常数据剔除后,两台仪器的线性拟合方程为CEC_Sunset= 0.6810×CEC_PKU+1.284(TOT法)和CEC_Sunset= 0.7035×CEC_PKU+1.1179(TOR法)(图 7),产生的正截距进一步说明Sunset-OCEC检测限偏高;查阅仪器说明书可知,Sunset-OCEC分析仪检测限为0.5 μg·m-3(以C计),可以推定本次观测使用的Sunset-OCEC分析仪的检测限高于说明书参数.

图 6 PKU-OCEC分析仪和Sunset-OCEC分析仪测得的EC浓度比对(a.PKU-OCEC TOT法;b.PKU-OCEC TOR法) Fig. 6 Comparison of EC concentration between PKU-OCEC and Sunset-OCEC(a.TOT method with PKU-OCEC,b.TOR method with PKU-OCEC)

图 7 PKU-OCEC分析仪和Sunset-OCEC分析仪测得的EC浓度(剔除异常值后)的比对(a.PKU-OCEC TOT法;b.PKU-OCEC TOR法) Fig. 7 Comparison of EC concentration between PKU-OCEC and Sunset-OCEC(without the outliers)(a.TOT method with PKU-OCEC,b.TOR method with PKU-OCEC)

将PKU-OCEC分析仪和Sunset-OCEC分析仪测得的EC浓度分别与MAAP测得的BC浓度进行拟合,结果如图 8所示,图 8c中圆圈内数据点即为Sunset-OCEC测得的EC异常数据点,剔除后如图 8d所示.可以看出,两台仪器测得的EC与MAAP测得的BC均显著相关,PKU-OCEC测得的EC与MAAP测得的BC的线性拟合度R2分别为0.8868(TOT法)和0.8614(TOR法),高于Sunset-OCEC与MAAP的拟合结果;两台仪器测得的EC与MAAP测得的BC的拟合直线斜率均为1.2左右,即BC浓度普遍高于EC浓度,这和Sheridan等(2005)的研究结果一致,其认为光学法测量的结果实际上是EC和吸光性OC的总和;但在EC浓度较低时,PKU-OCEC和Sunset-OCEC均有EC结果高于BC的情况出现,推测其原因,一方面可能是仪器在低浓度时测量结果受噪音响应影响较大;另一方面,由于污染源并不稳定,不同浓度EC对应的气溶胶混合状态可能发生变化,Watson等(2005)指出,当EC从气溶胶球体的核心迁移到表面后,其消光系数可能减小一个数量级;因此,若观测期间低浓度EC对应的气溶胶消光系数较低,可能导致MAAP按理论消光系数计算得到的BC浓度低于实际值.

图 8 PKU-OCEC分析仪和Sunset-OCEC分析仪测得的EC浓度与MAAP测得的BC浓度比对(a.PKU-OCEC TOT法,b.PKU-OCEC TOR法,c.Sunset-OCEC原始数据,d. Sunset-OCEC剔除异常点后数据) Fig. 8 Inter-comparison of EC(BC)concentrations measured with PKU-OCEC,Sunset-OCEC and MAAP(a. TOT method with PKU-OCEC,b. TOR method with PKU-OCEC,c. Sunset-OCEC original data,d. Sunset-OCEC data without the outliers)
5 结论(Conclusions)

1) 基于目前商品化在线OCEC分析仪在结构和操作上存在的一些问题,研制了一套大气气溶胶碳质组分在线分析仪,仪器同时采用TOT法和TOR法分离OC和EC,方便实验室之间的数据比对;自动化控制软件可以快速实现参数更改及数据查询功能;气路设计保证了分析阶段系统内载气流量的恒定,且将定量环通气时间缩短为45 s,节约了95%以上的He/CH4气,极大地延长了He/CH4钢瓶气的使用寿命;条件实验表明,仪器流量控制系统及内标气体响应稳定,蔗糖标定曲线线性拟合度高(R2=0.9990),仪器检测限为0.29 μg·m-3(以C计)(采样流量8 L·min-1,采样时间40 min),测量范围最大值为211.5 μg·m-3(以C计),可以满足不同环境的大气观测要求.

2) 利用该仪器于2016年1月23—31日对北京怀柔大气进行监测,并与Sunset-OCEC分析仪和MAAP进行数据比对,结果显示:PKU-OCEC测得的TC、OC浓度和Sunset-OCEC的对应结果显著相关(TC:r=0.986,OC:r=0.984,p<0.01),说明两台仪器测量结果具有很好的一致性;Sunset-OCEC测得的观测期间EC数据有39.39%浓度近似为0 μg·m-3,剔除该部分异常值后,PKU-OCEC与Sunset-OCEC相应数据的线性拟合方程由CEC_Sunset=0.9278×CEC_PKU-0.0704(TOT法)和CEC_Sunset=0.9476×CEC_PKU-0.2665(TOR法)变为CEC_Sunset=0.6810×CEC_PKU+1.2842(TOT法)和CEC_Sunset=0.7035×CEC_PKU+1.1179(TOR法),新的方程产生出较大正截距,说明Sunset-OCEC的检测限偏高;在与MAAP的比对中发现,PKU-OCEC测得的EC与MAAP测得的BC显著相关,其线性拟合度(R2=0.8868(TOT法),R2=0.8614(TOR法))大于Sunset-OCEC和MAAP的拟合度(R2=0.8047);PKU-OCEC分析仪两种方法测得的OC浓度相当且显著相关(COC_TOR= 0.9963×COC_TOTR2=0.9959,r=0.998,p<0.01),说明分析过程中所测样品在膜表面和厚度方向的黑度变化较为一致.

3) 综合结果表明,PKU-OCEC分析仪在实际外场观测中和商品化仪器有较好的比对效果,测量数据可靠,检测限较低,运行稳定且操作简便,可以广泛地用于各种条件下的环境监测.

参考文献
[${referVo.labelOrder}] Birch M E, Cary R A. 1996. Elemental carbon-based method for monitoring occupational exposures to particulate diesel exhaust[J]. Aerosol Science &Technology, 25(3) : 221–241.
[${referVo.labelOrder}] Bond T C, Bergstrom R W. 2006. Light absorption by carbonaceous particles:An investigative review[J]. Aerosol Science &Technology, 40(40) : 27–67.
[${referVo.labelOrder}] Cao J J, Lee S C, Ho K F, et al. 2003. Characteristics of carbonaceous aerosol in Pearl River Delta Region,China during 2001 winter period[J]. Atmospheric Environment, 37(11) : 1451–1460. DOI:10.1016/S1352-2310(02)01002-6
[${referVo.labelOrder}] Chen L W A, Chow J C, Watson J G, et al. 2004. Modeling reflectance and transmittance of quartz-fiber filter samples containing elemental carbon particles:Implications for thermal/optical analysis[J]. Journal of Aerosol Science, 35(6) : 765–780. DOI:10.1016/j.jaerosci.2003.12.005
[${referVo.labelOrder}] Chow J C, Watson J G. 2002. PM2.5 carbonate concentrations at regionally representative interagency monitoring of protected visual environment sites[J]. Journal of Geophysical Research-Atmospheres, 107(D21) : 8346–8349.
[${referVo.labelOrder}] Chow J C, Watson J G, Chen L W A, et al. 2004. Equivalence of elemental carbon by thermal/optical reflectance and transmittance with different temperature protocols[J]. Environmental Science &Technology, 38(16) : 4414–4422.
[${referVo.labelOrder}] Chow J C, Watson J G, Crow D, et al. 2001. Comparison of IMPROVE and NIOSH carbon measurements[J]. Aerosol Science and Technology, 34(1) : 23–34. DOI:10.1080/02786820119073
[${referVo.labelOrder}] Chow J C, Watson J G, Louie P K K, et al. 2005. Comparison of PM2.5 carbon measurement methods in Hong Kong,China[J]. Environmental Pollution, 137(2) : 334–344. DOI:10.1016/j.envpol.2005.01.006
[${referVo.labelOrder}] Chow J C, Watson J G, Pritchett L C, et al. 1993. The dri thermal/optical reflectance carbon analysis system:description,evaluation and applications in U. U.S.Air quality studies[J]. Atmospheric Environment, 27 : 1185–1201. DOI:10.1016/0960-1686(93)90245-T
[${referVo.labelOrder}] Gray H A, Cass G R, Huntzicker J J, et al. 1986. Characteristics of atmospheric organic and elemental carbon particle concentrations in Los Angeles[J]. Environmental Science &Technology, 20(6) : 580–589.
[${referVo.labelOrder}] HitzenberGer R, Larson S, Dillner A, et al. 1999. Intercomparison of measurement methods for black carbon aerosols[J]. Atmospheric Environment, 33(17) : 2823–2833. DOI:10.1016/S1352-2310(98)00360-4
[${referVo.labelOrder}] Huntzicker J J,Johnson R L,Shah J J,et al.1982.Analysis of Organic and Elemental Carbon in Ambient Aerosols by a Thermal-Optical Method//Wolff G T,Klimisch R L.Particulate Carbon:Atmospheric Life Cycle[M].Boston,MA:Springer.79-88
[${referVo.labelOrder}] 刘巧铃. 2010.在线气溶胶碳质组分(OCEC)分析仪的研制[D].北京:北京大学.19-22
[${referVo.labelOrder}] Matsumoto K, Tanaka H, Nagao I, et al. 1997. Contribution of particulate sulfate and organic carbon to cloud condensation nuclei in the marine atmosphere[J]. Geophysical Research Letters, 24(6) : 655–658. DOI:10.1029/97GL00541
[${referVo.labelOrder}] Nemmar A, Hoet P H M, Vanquickenborne B, et al. 2002. Passage of inhaled particles into the blood circulation in humans[J]. Circulation, 105(4) : 411–414. DOI:10.1161/hc0402.104118
[${referVo.labelOrder}] Offenberg J H, Baker J E. 2002. The influence of aerosol size and organic carbon content on gas/particle partitioning of polycyclic aromatic hydrocarbons (PAHs)[J]. Atmospheric Environment, 36(7) : 1205–1220. DOI:10.1016/S1352-2310(01)00427-7
[${referVo.labelOrder}] Polidori A, Turpin B J, Ho-Jin L, et al. 2006. Local and regional secondary organic aerosol:Insights from a year of semi-continuous carbon measurements at Pittsburgh[J]. Aerosol Science &Technology, 40(10) : 861–872.
[${referVo.labelOrder}] Steinfeld J I. 1998. Atmospheric chemistry and physics:From air pollution to climate change[J]. Environment:Science and Policy for Sustainable Development, 40(7) : 26–26.
[${referVo.labelOrder}] Sheridan, P J, Arnott W P, et al. 2005. The reno aerosol optics study:An evaluation of aerosol absorption measurement methods[J]. Aerosol Sci Tech, 39(1) : 1–16. DOI:10.1080/027868290901891
[${referVo.labelOrder}] Tanner R L, Gaffney J S, Phillips M F. 1982. Determination of organic and elemental carbon in atmospheric aerosol samples by thermal evolution[J]. Anal Chem, 9(9) : 1627–1630.
[${referVo.labelOrder}] Turpin B J, Huntzicker J J. 1991. Secondary formation of organic aerosol in the Los Angeles basin:A descriptive analysis of organic and elemental carbon concentrations[J]. Atmospheric Environment Part A General Topics, 25(2) : 207–215. DOI:10.1016/0960-1686(91)90291-E
[${referVo.labelOrder}] Watson J G, Chow J C, Chen L W A. 2005. Summary of organic and elemental carbon/black carbon analysis methods and intercomparisons[J]. Aerosol &Air Quality Research, 5(1) : 65–102.
[${referVo.labelOrder}] 吴琳, 冯银厂, 叶文媛, 等. 2010. 大气颗粒物中碳组分测定结果比较:元素分析和热光反射方法[J]. 环境科学研究, 2010, 23(12) : 1481–1487.
[${referVo.labelOrder}] 郇宁, 曾立民, 邵敏. 2005. 气溶胶中有机碳及元素碳分析方法进展[J]. 北京大学学报 (自然科学版), 2005, 41(6) : 957–964.