气溶胶单粒子化学成分在线聚类分析文档格式.docx

气溶胶单粒子化学成分在线聚类分析文档格式.docx

《气溶胶单粒子化学成分在线聚类分析文档格式.docx》由会员分享，可在线阅读，更多相关《气溶胶单粒子化学成分在线聚类分析文档格式.docx（8页珍藏版）》请在冰点文库上搜索。

本算法的质谱分析软件基本满足大气气溶胶单粒子在线聚类分析的要求。

【关键词】大气气溶胶;

气溶胶单粒子;

化学成分;

气溶胶激光飞行时间质谱仪;

自适应共振神经网络

　　AbstractToquicklyprocessthelaserdesorption/ionizationmassspectrawhichweregeneratedbyaerosollasertimeofflightmassspectrometer（ALTOFMS）duringitsworkingandextractthevaluableinformation,anadaptiveresonancetheorybasedneuralnetwork,ART2aalgorithm,wassuccessfullyappliedtotherealtimeclassificationofsingleparticlesofNaCl,themixtureofNaClandCaCl2aerosol,melamineandtheatmosphericaerosol.Experimentalresultsshowedthatwhenvigilancefactorwas0.1andlearningratewas0.05,ART2aalgorithmcouldsuccessfullyrevealtheaerosolparticlescategories.Besidesthesuccessfulprobabilitywasnear100%,thecentroidmassspectraforthesingleparticleclasseswereobtained,whichcouldrepresentthecharacteristicofsingleparticleclassesremarkably.ThenumberofNaClsingleparticleclassesasafunctionofvigilancefactorwasalsodiscussed.TheresultshowedthatthenumberofNaClparticleclasseswaslargerwhenvigilancefactorreached0.8;

furthermoretheclassificationwasmuchmoreprecise.ThemassspectraacquisitionandcontrolsoftwareusingART2acanbasicallymeettherequirementsofrealtimeclassificationofatmosphericaerosolsingleparticles.

　　KeywordsAtmosphericaerosol;

Aerosolsingleparticle;

Chemistrycomposition;

Aerosollasertimeofflightmassspectrometer;

Adaptiveresonanceneuralnetwork

　　1引言

　　大气气溶胶在全球气候变化、全球和区域污染过程、降低能见度等方面起到重要作用，甚至影响人类的生存环境。

获知气溶胶单粒子的粒径和化学成分将有助于解析气溶胶的来源、传输过程及其在大气中的化学反应[1～3]。

气溶胶激光飞行时间质谱仪（Aerosollasertimeofflightmassspectrometer，ALTOFMS）能够对单个气溶胶粒子的粒径和化学成分进行实时测量[4,5]。

该仪器在使用过程中会产生海量数据（每分钟达上百个质谱数据），依靠人工处理数据不仅速度慢、劳动强度大，还易出现误差，极大地削弱了仪器实时测量的优势。

因此，开发高效全面的在线质谱分析软件非常必要。

　　目前已用于质谱数据分析的方法多为离线分析，如等级聚类分析（HCA）、主成分分析（PCA）、模糊C均值（FCM）[5]、自组织特征映射神经网络（SOM）和自适应共振算法（ART2a）[6～8]。

本课题组利用FCM和ART2a算法对气溶胶质谱进行离线聚类分析[5,8]，取得了较好的效果，但无法满足实时在线快速提取信息的要求。

目前，在线聚类分析质谱数据方面的工作鲜有报道。

　　ART2a神经网络算法与其它算法相比，其最大的特点是能快速识别遇到过的质谱模式类型，即存储在数据库中的聚类中心所代表的模式类，并且能够根据环境输入，自动将不同于已知质谱模式类但又具有相似特征的输入归为新的一类，并将学习到的这个新类的聚类中心保存下来，作为一种新的模式类参与以后的分类。

简单而言，该算法具有无需监督、可灵活控制分类精细程度、快速有效识别已记忆质谱模式类以及根据环境变化自适应快速识别和学习新类的特点，非常适合气溶胶粒子化学成分的在线聚类分析。

本研究在前期离线聚类分析基础上，利用ART2a算法开发实时的化学成分分析软件，将其加载到ALTOFMS质谱采集控制系统中，并对实验室环境下和大气中的气溶胶单粒子进行在线分类，对粒子分类后的分类信息在采集控制系统中实时显示。

实验证明，此算法可对气溶胶单粒子的化学成分进行在线聚类分析，并可以较好地运用到整机系统中。

　　2实验部分

　　2.1实验装置

　　本实验是在自行研制的气溶胶激光飞行时间质谱仪上完成的。

该仪器的工作原理及仪器装置详见文献[3～5,9]。

采用空气动力学透镜加两级差分系统进样、双束连续激光测量气溶胶单粒子粒径和激光飞行时间质谱仪进行化学成分分析。

使用德维尔比斯（Devilbiss40#）产生气溶胶粒子，方法参见文献[8]。

　　2.2数据预处理

　　在线聚类分析中的ART2a算法是离线聚类分析中ART2a算法的简化版本，前者省略了后者中的重复迭代过程，缩小了后者中学习速率的作用范围。

这是为了保证存放在数据库中的聚类中心的稳定性，其余步骤相同，见文献[8]。

简化的目的是在保证数据库中聚类中心稳定的情况下，通过牺牲少量的分类精确度，实现实时分类当前采集到的气溶胶粒子。

在线仿真时单粒子的分析速度可以达到50个/s，而单粒子的采集速率仅为5个/s，因此，在线聚类分析的速度完全可以跟得上粒子的采集速率。

　　3结果与讨论

　　3.1NaCl气溶胶单粒子在线聚类分析

　　根据离线聚类分析的经验值[7,9]，本实验中警戒阈值和学习速率分别取0.1和0.05。

　　即便数据库中不存放NaCl质谱聚类中心，软件也能自行学习，将粒子质谱很好地聚成3类：

图1所示为NaCl单粒子LDI正离子质谱在线聚类中心，Class1在m/z23（Na+）和m/z81（Na235Cl+）处有明显的特征峰，占总质谱数的57.5%。

将其存入数据库中，命名为NaClⅠ类;

Class2中NaCl特征峰全部向后平移了1Da，从而被聚成一个新类，占25.7%，这种峰值漂移不是新质量数谱峰的出现，而是离子在穿越质谱仪的离子引出区和加速区极板时发生部分离子吸附，引起极板间电压起伏造成的[11]。

正是这种飞行时间质谱峰值的漂移导致了两个来自相同物质的粒子在某个或全部质量数上相差1Da;

Class3虽然完全包含了前面Class1中的所有离子种类，但它比Class1多出了m/z46（Na2+），m/z83（Na237Cl+）和m/z139（Na335Cl2+）3个离子特征峰，同样属于NaCl类，记为NaClⅡ类，占3.5%。

这与质谱数据的采集次数无关，是NaCl质谱种类本身的多样性造成的。

命中质谱113个，能够成功聚类的有98个，另外15个质谱（占总质谱数的13.3%）因Na+信号峰太强已经过饱和，没有进行聚类，如果改变软件设置的增益使信号减弱，可以全部分类。

当数据库中加载了前面获得的两类聚类中心（NaClⅠ和NaClⅡ）时，同样可以很好地聚成3类：

NaClⅠ，NaClⅡ和Newcluster1。

其中Newcluster1即为图1中的Class2。

由此可见,软件的自主学习能力是非常强大的。

　　3.2警戒阈值对NaCl粒子在线聚类数的影响

　　本实验考察了警戒阈值参数。

实验发现，加大警戒阈值，可以增加在线分类数、提高分类的精细程度。

图2为不同警戒阈值与聚类数的关系。

当警戒阈值0.7时，可以在线分为3类（见3.1节）;

警戒阈值≥0.8时，聚类数目逐渐增多;

警戒阈值=0.95时，聚类数增加到12类，且聚类中心更加精细。

由于某个质谱峰相差1Da都会被细化并归到不同的类别中，表明软件具有很好的学习能力。

　　3.3NaCl和CaCl2气溶胶粒子混合物在线聚类分析

　　用Devilbiss40#分别产生多分散的NaCl气溶胶粒子和CaCl2气溶胶粒子，将其混合在10L棕色瓶中，用ALTOFMS进行检测。

实验共获得255个有效质谱数据，分为5类，各聚类的粒子个数（所占比例）、主要质谱峰和物质类型如表1所示。

表1基于ART2a算法的聚类概述（略）

　　本实验的前提条件是所加载的数据库中只包含NaCl类而没有CaCl2类：

当实验中采集到的质谱特征峰（如23，46，81和83）只属于已有数据库中事先存放好的NaCl类时，算法能快速识别出来，自动归属到NaCl类中;

当采集到的质谱峰（m/z40，44，75和77）不属于NaCl类时，ART2a算法会自动识别和学习特征峰并形成聚类中心Newcluster1、2和3。

实验结束后，只需分析3个聚类中心，并根据经验判断Newcluster2和3可能属于CaCl2类，因为这两个聚类中心包含了CaCl2粒子的最常见的特征质谱峰（40Ca2+，44Ca2+，40Ca35Cl+和40Ca37Cl+），能较好地符合CaCl2颗粒的激光解吸附电离特性;

Newcluster1中的质谱峰主要以40和23为主，由前面的分析就可以推断出粒子可能是NaCl和CaCl2的混合物。

当NaCl和CaCl2以气溶胶形式相互混合时，会发生凝并，进入ALTOFMS的气溶胶粒子产生的质谱有一部分以Na+和Ca+形式共同存在。

由于离子的空间电荷效应导致质谱峰展宽[12]，仪器分辨率有时比较差，NaCl与CaCl2气溶胶粒子混合后的飞行时间质谱中的Na+和Ca+峰跨越了几个质量数，因此获得的聚类中心在m/z40附近有m/z35，37和38的峰，并不是杂质的影响。

因此，当采集到的粒子质谱数据中出现新的质谱模式类时，ART2a算法可以很好地适应环境输入，具有自动识别和学习新的质谱模式类型的能力。

利用此算法开发的在线质谱分析软件可以完成质谱数据的在线预处理，反应速度快、误差小，从而使得需要人工处理的数据量明显减少，这极大提高了质谱仪实时检测气溶胶粒子的效率。

　　3.4三聚氰胺气溶胶单粒子在线聚类分析

　　以三聚氰胺为例研究有机气溶胶单粒子的在线聚类分析。

图3a为三聚氰胺（C3H6N6）典型的正、负离子解吸/电离（LDI）质谱。

正离子质谱中，m/z127为分子离子（[C3H6N6+H]+），m/z85为三聚氰胺开环分子重排失去H2NCN的碎片离子峰（[C3H6N6+H-H2NCN]+），m/z68为m/z85脱去NH3的碎片离子峰（[C3H6N6+H-H2NCN-NH3]+），m/z60应为三聚氰胺开环分子重排失去HNCNCN的碎片离子峰（[C3H6N6+H-HNCNCN]+），m/z43应为m/z60脱去NH3的碎片离子峰（[C3H6N6+H-HNCNCNNH3]+）。

负离子质谱中m/z26是可以标识的，为CN-。

图3b为三聚氰胺单粒子的在线聚类中心：

负离子质谱在m/z26,41和66处有明显信号，这与其质谱完全吻合;

正离子质谱在m/z18,23,39和43处有明显信号，m/z127处没有信号，原因是大多数有机气溶胶粒子的LDI谱峰复杂多样，且常存在金属离子的干扰，分子离子峰本身质量数比较大，探测效率很低（三聚氰胺分子离子峰的探测效率仅有0.2%），软件未分析出来，这就是牺牲少量的分类精确度，换取质谱的实时分类。

随着实验方法和实验条件的完善和提高，分子离子峰的探测效率会相应优化，此软件可以很好地应用到有机气溶胶粒子的实时分类。

　　3.5大气气溶胶单粒子化学成分在线聚类分析

　　在前面研究的基础上，对大气气溶胶单粒子进行实时在线聚类分析。

图4为实验过程中截取的质谱采集控制系统软件界面。

在质谱采集控制系统[13]基础上增加了实时的化学成分分析模块，对粒子分类后的分类信息在图5左侧中间的列表控件中实时显示。

当采集到的质谱不属于已加载的质谱数据库时，在线聚类名称以newcluster1,newcluster2,newcluster3……命名。

命中的103个粒子被很好地聚为4类，被成功分类的粒子数占命中粒子总数的100%。

　　4结论

　　本实验利用ART2a算法编制软件功能模块，实现了对气溶胶激光飞行时间质谱数据的在线分类。

运用此算法对实验室环境下和大气中的气溶胶单粒子进行了在线聚类分析。

当警戒阈值为0.1、学习速率为0.05时，能够快速区分大气气溶胶粒子，被成功分类的粒子数占命中粒子总数的100%。

获得的每种粒子的聚类中心都能较好地代表该类物质特征，对于有机气溶胶单粒子而言，随着实验方法和实验条件的完善和提高，分子离子峰的探测效率会相应优化，此软件可以很好地应用到有机气溶胶粒子的实时分类。

利用ART2a算法开发的软件模块具有很好的自主学习和自动识别学习新类的能力，可以应用于气溶胶激光飞行时间质谱仪海量数据的在线聚类分析。

【参考文献】

　　1SuessDT,PratherKA.Chem.Rev.,1999,99（10）:

3007～3035

　　2XUSiQi（徐思琦）,XIEZhouQing（谢周清）,LIUWei（刘崴）,YANGHongXia（杨红霞）,LIBing（李冰）.ChineseJ.Anal.Chem.（分析化学）,2010,38

（2）:

219～224

　　3XIAZhuHong（夏柱红）,FANGLi（方黎）,ZHENGHaiYang（郑海洋）,KONGXiangHe（孔祥和）,ZHOULiuZhu（周留柱）,GUXueJun（顾学军）,ZHUYuan（朱元）,ZHANGWeiJun（张为俊）.ChineseJ.Anal.Chem.（分析化学）,2004,32（7）:

973～976

　　4ZHOULiuZhu（周留柱）,ZHUYuan（朱元）,GUOXiaoYong（郭晓勇）,ZHAOWenWu（赵文武）,ZHENGHaiYang（郑海洋）,GUXueJun（顾学军）,FANGLi（方黎）,ZHANGWeiJun（张为俊）.ChineseJ.Anal.Chem.（分析化学）,2005,33（11）:

1551～1554

　　5GUOXiaoYong（郭晓勇）,FANGLi（方黎）,ZHAOWenWu（赵文武）,GUXueJun（顾学军）,ZHENGHaiYang（郑海洋）,ZHANGWeiJun（张为俊）.SpectroscopyandSpectralAnalysis（光谱学与光谱分析）,2008,28（8）:

1713～1717

　　6SongXH,HopkePK,FergensonDP,PratherKA.Anal.Chem.,1999,71（4）:

860～865

　　7MurphyDM,MiddlebrookAM,WarshawskyM.AerosolSci.Tech.,2003,37（4）:

382～391

　　8LINYing（林莺）,GUOXiaoYong（郭晓勇）,GUXueJun（顾学军）,XIAWeiWei（夏玮玮）,ZHENGHaiYang（郑海洋）,ZHANGWeiJun（张为俊）,FANGLi（方黎）.SpectroscopyandSpectralAnalysis（光谱学与光谱分析）,2009,29（3）:

580～584

　　9GuoXY,ZhaoWW,HuangMQ,FangL,ZhangWJ.OpticaApplicata,2007,37（3）:

229～235

　　10DenisJP,KevinPR,AnthonySW,DavidBK,MurrayVJ.Anal.Chem.,2001,73（10）:

2338～2344

　　11ZHANGShuDong（张树东）,XUMingKun（许明坤）,KONGXiangHe（孔祥和）.SpectroscopyandSpectralAnalysis（光谱学与光谱分析）,2007,27（4）:

799～802

　　12SHIYong（石勇）,ZHOUXiaoGuo（周晓国）,ZHAOYingGuo（赵英国）,LIUShiLin（刘世林）,MAXingXiao（马兴孝）.ActaPhys.Chim.Sin.（物理化学学报）,2005,21（11）:

1269～1273

　　13GUOXiaoYong（郭晓勇）,ZHUYuan（朱元）,ZHOULiuZhu（周留柱）,ZHAOWenWu（赵文武）,ZHOUZhiGuang（周志广）,ZHENGHaiYang（郑海洋）,FANGLi（方黎）,ZHANGWeiJun（张为俊）.ChineseJournalofQuantumElectronic（量子电子学报）,2006,23

（2）:

217～221