CN103940934B - 一种分析混合物成分的方法 - Google Patents
一种分析混合物成分的方法 Download PDFInfo
- Publication number
- CN103940934B CN103940934B CN201410113885.6A CN201410113885A CN103940934B CN 103940934 B CN103940934 B CN 103940934B CN 201410113885 A CN201410113885 A CN 201410113885A CN 103940934 B CN103940934 B CN 103940934B
- Authority
- CN
- China
- Prior art keywords
- spectrogram
- analysis
- spectrum
- pure
- potpourri
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000203 mixture Substances 0.000 title claims abstract description 145
- 238000004458 analytical method Methods 0.000 title claims abstract description 55
- 238000001228 spectrum Methods 0.000 claims abstract description 146
- 238000000034 method Methods 0.000 claims abstract description 134
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 51
- 238000004587 chromatography analysis Methods 0.000 claims abstract description 19
- 230000014759 maintenance of location Effects 0.000 claims abstract description 3
- 239000000126 substance Substances 0.000 claims description 27
- 238000001819 mass spectrum Methods 0.000 claims description 17
- 238000005457 optimization Methods 0.000 claims description 15
- 239000000463 material Substances 0.000 claims description 10
- 238000005516 engineering process Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 6
- 238000005481 NMR spectroscopy Methods 0.000 claims description 2
- 238000004811 liquid chromatography Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 abstract description 18
- 238000000926 separation method Methods 0.000 abstract description 6
- 238000005070 sampling Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 14
- 239000003350 kerosene Substances 0.000 description 9
- UFWIBTONFRDIAS-UHFFFAOYSA-N Naphthalene Chemical compound C1=CC=CC2=CC=CC=C21 UFWIBTONFRDIAS-UHFFFAOYSA-N 0.000 description 8
- 150000001875 compounds Chemical class 0.000 description 8
- 229940094933 n-dodecane Drugs 0.000 description 8
- 239000013598 vector Substances 0.000 description 8
- SNRUBQQJIBEYMU-UHFFFAOYSA-N dodecane Chemical compound CCCCCCCCCCCC SNRUBQQJIBEYMU-UHFFFAOYSA-N 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 238000004817 gas chromatography Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- -1 n-undecane hydrocarbon Chemical class 0.000 description 5
- 239000000470 constituent Substances 0.000 description 4
- 239000003814 drug Substances 0.000 description 4
- 238000002451 electron ionisation mass spectrometry Methods 0.000 description 4
- 230000008676 import Effects 0.000 description 4
- 238000002203 pretreatment Methods 0.000 description 4
- 150000001335 aliphatic alkanes Chemical class 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 description 3
- 238000004949 mass spectrometry Methods 0.000 description 3
- 238000012067 mathematical method Methods 0.000 description 3
- RSJKGSCJYJTIGS-UHFFFAOYSA-N undecane Chemical compound CCCCCCCCCCC RSJKGSCJYJTIGS-UHFFFAOYSA-N 0.000 description 3
- 240000002234 Allium sativum Species 0.000 description 2
- 101100001678 Emericella variicolor andM gene Proteins 0.000 description 2
- 235000004611 garlic Nutrition 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229930014626 natural product Natural products 0.000 description 2
- 125000002524 organometallic group Chemical group 0.000 description 2
- 239000000376 reactant Substances 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- CVOFKRWYWCSDMA-UHFFFAOYSA-N 2-chloro-n-(2,6-diethylphenyl)-n-(methoxymethyl)acetamide;2,6-dinitro-n,n-dipropyl-4-(trifluoromethyl)aniline Chemical compound CCC1=CC=CC(CC)=C1N(COC)C(=O)CCl.CCCN(CCC)C1=C([N+]([O-])=O)C=C(C(F)(F)F)C=C1[N+]([O-])=O CVOFKRWYWCSDMA-UHFFFAOYSA-N 0.000 description 1
- 241000234282 Allium Species 0.000 description 1
- 235000002732 Allium cepa var. cepa Nutrition 0.000 description 1
- 206010038743 Restlessness Diseases 0.000 description 1
- 150000004945 aromatic hydrocarbons Chemical class 0.000 description 1
- 125000003118 aryl group Chemical group 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 238000009614 chemical analysis method Methods 0.000 description 1
- 238000013375 chromatographic separation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000008267 milk Substances 0.000 description 1
- 210000004080 milk Anatomy 0.000 description 1
- 235000013336 milk Nutrition 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
- G01N30/8679—Target compound analysis, i.e. whereby a limited number of peaks is analysed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8631—Peaks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Library & Information Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及一种分析混合物成分的方法,包括:(1)将混合物样品,利用色谱技术进行分离,得到初步的色谱图;(2)对步骤(1)得到的初步的色谱图中的任意一个区间,根据需要进行采样,得到一系列在不同保留时间的谱图,称之为混合物谱图;(3)利用熵最小算法的系列方法,从步骤(2)得到的混合物谱图中,计算得到各重建纯谱及其对应的伪浓度。本发明提供的方法,将色谱技术和熵最小算法结合,主要克服了现有分析方法对成分分离的依赖,使得分离不再重要;不需要彻底地分离混合物中各个成分,同时也达到了熵最小算法使用的先决条件;并且方法通用、速度快、高效、成本低,对人员要求不高,具有巨大的应用前景。
Description
技术领域
本发明涉及化学分析领域,具体地,涉及一种分析混合物成分的方法。
背景技术
现有分析方法中,分析混合物中的各成分主要是利用各种色谱分析方法,即利用混合物中各成分在各种色谱中流动相和固定相上的化学、物理性质的不同,使得各成分在色谱中移动速度产生差别,从而使各成分在不同的时间经过检测器(如紫外、红外、质谱等),获得各种谱图,从而分离、分辨混合物中的各成分。
混合物全成分的分析现状:
对复杂混合物中各成分进行分析,常用分析方法是在多次不同条件(例如改变温度、改变温度变化方式、改变流速、使用不同的固定相(如更换色谱柱)和不同的流动相等)的色谱分析中,每次将一部分的混合物分离,然后综合各个不同条件下的结果,得到较为完全的分析结果。现有方法经过多次试验,仍然不能完全知道是否全部分离了成分。
现有分析方法对复杂混合物的全成分分离非常耗时、耗力和耗费,而且不通用。通常对一种复杂混合物的分析,例如大蒜中成分的分析,需要多年的研究,同时需要几百万甚至几千万的投入,并需要很多专业人员。而且,由于不同混合物中的成分不同(matrix),所以现有分析方法也不是通用的。例如大蒜成分的分析方法,对于洋葱就不适合了,需要重复分析,还是需要大量资金、时间和人员的投入。
现有分析方法对复杂混合物全成分分析的最大困难为:分离,就是如何将混合物中各个成分一个一个分离开。对于稍微复杂的混合物,将混合物在几次分析中全部分离开是不可能的,出现重叠谱图(混合物谱图)是不可避免的,更不要说有着成千上万成分的中药和天然产物了。如何尽快地分离开混合物中的各个成分,从而对其中各个成分进行分析,这是个世界性的难题。
标的物的分析现状:
对于含有标的物的分析,例如食品、药物或环境安全等的分析,现有分析方法一般需要将标的样品尽量提纯,或者将样品中的标的物做物化处理,例如氘代反应等;这些前处理一般都比较繁琐、费事。然后,处理好的样品需要经过色谱分离,将标的物从背景化合物中(martrix)分离开,然后才能定性和定量。
由于各种标的物不同,导致检测、分析各种标的物的过程中,所需要的仪器和条件等都不一样,所以一般的分析中心只能对于一些特定的标的物做分析检测,而且分析成本都会比较高。
标的物的分析,其根本的问题,和混合物全成分分析的根本问题一样,都是分离的问题。所以,解决了分离的问题,就是解决分析混合物成分的核心问题。
熵最小算法(EntropyMinimizationAlgorithms)
在分析仪器和各种分析方法中,大量应用了各种数学方法,帮助解决分析中的各种问题,例如基线的问题等。这些数学方法统称为化学计量法。
有一类化学计量法主要是用于从混合物谱图中提取出纯成分的信息(例如SIMPLSMA,IPCA,OPA-ALS)。但是这些方法需要一些先决条件和信息,例如重叠谱图中每个纯成分的峰型是对称的(例如高斯对称),每个纯成分峰之间需要有一定的间隔,重叠程度也不能太大等;另外还需猜测混合物中有几个成分等。由于先决条件太多,这些数学方法对混合物中的纯成分分离帮助有限,只能起辅助作用。尤其是对于复杂混合物,例如中药和天然产物的样品,里面的成分信息是猜不出的,分析产生的重叠谱图也不能达到这些方法的应用要求。
熵最小算法(EM)是基于ShannonEntropy(香农熵)发展而来的。在没有提供任何成分信息的情况下,只需要混合谱数据,就能从混合谱中提取出每个纯谱。ShannonEntropy原理最早发表于1948年[C.E.Shannon,TheBellSystemTechnicalJournal,27(1948)379-423.],它是一个信息领域的学术用语,用于衡量随机参数的不确定性。
最早将ShannonEntropy应用于化学分析的是MarcGarland[Y.Z.Zeng,M.Garland,AnalyticaChimicaActa,359(1998)303-310.],他在2002年发表了BTEM(Band-TargetEntropyMinimization)[W.Chew,E.Widjaja,M.Garland,Organometallics,21(2002)1982-1990.],即锁峰熵最小方法。该方法是用于分离简单的红外混合谱。BTEM需要人为指导操作,主要应用于研究领域。
2003年张华俊等发表了tBTEM方法(WeightedTwo-BandTargetEntropyMinimization[H.J.Zhang,M.Garland,Y.Z.Zeng,P.Wu,JAmSocMassSpectrom,14(2003)1295-1305.]),即加权多重锁峰熵最小方法,该方法开始主要应用在质谱上,同时一种自动分析、去除重合谱的方法也一同发表了。
2006年张华俊等发表了MREM[H.J.Zhang,W.Chew,M.Garland,AppliedSpectroscopy,61(2007)1366-1372.](Multi-ReconstructionEntropyMinimization)方法,即多重建熵最小方法。该方法用局部优化方法代替了全局优化方法,而且不需人为指定搜索范围,真正实现了自动搜索纯谱的功能。
2009年熵最小方法顺利应用于紫外谱图解析[F.Gao,H.J.Zhang,L.F.Guo,M.Garland,ChemometricsandIntelligentLaboratorySystems,95(2009)94-100.]。由于紫外谱图非常宽,峰形不对称,各成分重叠非常厉害,即使2~3个成分的紫外混合物谱图也是非常难以解析的。
将熵最小系列方法(BTEM、tBTEM和MREM等)应用于混合谱中,能在不知道任何成分的情况下,仅利用混合物的谱图信息,就能从混合谱中重建出纯谱。它们已经成功应用于红外、紫外、质谱和核磁共振图谱等谱图中。
EM系列方法应用的前提是:1)混合谱谱图的数目应大于混合谱中成分的数目;例如,一个混合物中含有10个成分,那么经计算机采样得到的混合谱谱图的数目就应该大于10;2)各成分在每个混合谱谱图中的浓度比例,应该是不同的;如果混合物谱图中的任意两个或者多个成分的比例在不同的采样时间都保持一致,那么这几个成分对于熵最小算法来说,就是不可分辨的。现有的文献中,各种EM系列方法的应用,是利用检测器检测一个动态反应的不同过程,在不同的反应时间采样,从而得到很多不同的混合谱,从而使得各个混合谱中各成分的比例不一致。例如:在一个不平衡的化学反应过程中,于20个不同的时间采样,得到20个混合谱;由于反应随着时间推进,所以反应物和生成物的浓度随着时间发生变化,从而导致反应物和生成物的在各个采样点的浓度比例不一致。
然而实际的化学分析中,一个混合物样品,例如一个中药样品或者一个牛奶样品,里面各个成分的浓度比例是恒定的,所以就无法满足上述EM系列方法应用的两个条件。如何从一个恒定浓度比例的样品中,得到很多不同比例的混合物谱图,从而利用EM方法进行分析,这又是一个亟待解决的问题。
发明内容
为了解决上述技术问题,本发明的目的是提供一种分析混合物成分的方法。该方法将色谱分离技术和熵最小系列算法结合起来,可以快速分析混合物的各成分。
本发明提供的分析混合物成分的方法,包括:
(1)将混合物样品,利用色谱技术进行分离,得到初步的色谱图;
(2)对步骤(1)得到的初步的色谱图中的任意一个区间,根据需要进行采样,得到一系列在不同保留时间的谱图,称之为混合物谱图(即混合谱);
(3)利用熵最小算法的系列方法(即EM系列方法),从步骤(2)得到的混合物谱图中,计算得到各重建纯谱及其对应的伪浓度。
其中,步骤(1)中,得到的初步的色谱图中,包括纯成分(已经分离出的成分)的峰和混合物的重叠峰。
进一步地,对于纯成分的峰,利用现有常规方法已经能够进行分析与鉴定,一般不需要利用本发明方法后续的步骤进行分析和处理,当然使用者也可以利用本发明方法后续的步骤进行分析;对于混合物的重叠峰,需要利用本发明上述方法后续的步骤(上述步骤(2)和步骤(3)进行分析。
其中,步骤(1)得到的初步的色谱图中的任意一个区间,是指初步的色谱图的一段,或者全部。
进一步地,步骤(2)中的数据采样,指将需要的区间数据从仪器产生的数据中提取。采样是否均匀由仪器本身决定。
其中,步骤(1)中,色谱技术为:气相色谱技术、液相色谱技术或者两者的组合。
其中,步骤(1)中,利用色谱技术进行分离,经检测器得到混合物的谱图。
进一步地,所述检测器选自:质谱、红外、紫外、荧光、分光光度仪和核磁共振检测器中的一种或几种的组合。
其中,步骤(3)中,熵最小算法的系列方法是个总称,其包括基于ShannonEntropy(香农熵)的各种方法,包括但不限制于BTEM(Band-TargetEntropyminimization)、tBTEM(WeightedTwo-BandTargetEntropyMinimization)或MREM(Multi-ReconstructionEntropyMinimization))等方法中的一种或几种的组合。
进一步地,利用BTEM或者tBTEM方法,每次运算,可以得到一个重建纯谱;利用MREM方法,每次运算,可以得到多个重建纯谱。
本发明中的“重建纯谱”指的是:利用熵最小算法的系列方法(即各种熵最小算法),从混合物谱图中,经过计算得到的纯谱。
本发明中的“各重建纯谱对应的伪浓度”指的是:利用熵最小算法的系列方法,在得到重建纯谱的时候,同时会计算得到一个和该重建纯谱相关的浓度(即叫做“伪浓度”)。
其中,步骤(3)中,对于含有未知物的混合物样品的分析,需要尽可能知道步骤(2)得到的混合物谱图中的成分信息的时候,利用熵最小算法的系列方法,在不同的参数下,反复多次计算后,得到尽可能多的重建纯谱。
上述方法中:
不同的参数:指锁定不同的峰,使用不同的目标函数和罚函数或使用不同的优化方法等。如:在一个m/z从100-200的混合物谱图中(间隔是1m/z),该数据共有101个m/z通道;可以将每个通道锁定,然后计算,或者任意两个通道进行计算;或者只对通道值大于30%最大峰的通道进行计算等。
计算多次:通过不同的参数设定,进行计算,每次都得到一个重建纯谱(BTEM或者tBTEM方法)或者多个重建纯谱(MREM方法)。由于不同的参数设定,导致多次的计算。
尽可能多的重建纯谱:由于多次计算后,会有很多重建纯谱,有些纯谱是相同的,所以可以利用算法去除,最后每个不同的纯谱都只有一个。
其中,步骤(3)中,对于含有标的物/已知物的混合物样品的分析,基于标的物/已知物的标准谱图,利用熵最小算法的系列方法计算,构建使用不同参数的重建方法,得到该标的物/已知物的重建纯谱;在多次的重建中,如果重建纯谱和标的物/已知物的标准谱图相比较是一致的,而且该重建纯谱的伪浓度是有意义的,则该混合物中含有该标的物/已知物,反之则不含有。
进一步地,多次的重建,例如:根据标的物纯谱,对于某个数据通道的值>30%最高值的,就进行锁住,然后做计算;或者可以每个通道都锁住一次,然后进行计算。
进一步地,重建纯谱的伪浓度是有意义的,是根据经验,在一般的色谱柱分析中,一般峰型是较对称的单峰,而不是如同随机值一样的。
本发明所述的方法,在步骤(3)后,可以将得到的重建纯谱和标准数据库中的谱图进行比对,确认它们的信息;若比对后重建纯谱和标准数据库中的某谱图一致,则说明重建纯谱代表的是标准数据库中的某谱图所代表的成分。若比对后重建纯谱和标准数据库中的谱图都不一致,则表明发现了标准数据库中没有记录的化学成分;后续可以再寻找其他数据库进行比对,或根据需要利用后续的各种分析方法鉴定该成分。
本发明提供的方法,有益效果如下:
1)主要克服了现有分析方法对成分分离的依赖,使得分离不再重要,克服了现有分析方法的技术瓶颈;
2)将色谱技术和熵最小算法结合,克服了彼此的缺点;
3)将色谱技术和熵最小算法结合连用后,不需要彻底地分离混合物中各个成分,同时也达到了熵最小算法使用的先决条件;
4)扩大了熵最小算法的应用范围,可以分析各成分比例恒定的混合物样品;
5)可以分析混合物样品中的未知成分,也可以分析含有已知成分/标的物的混合物,进而可以对各成分进行定性与定量;
6)是一个通用的分析方法,不同的混合物样品均适用;
7)速度快,因为不需要花非常多的时间在分离上;成本低,不需要十分优良的分离设备,对人员要求也低。
总之,本发明提供的分析混合物成分的方法,将色谱技术和熵最小算法结合,主要克服了现有分析方法对成分分离的依赖,使得分离不再重要;不需要彻底地分离混合物中各个成分,同时也达到了熵最小算法使用的先决条件;并且方法通用、速度快、高效、成本低,对人员要求不高,具有巨大的应用前景。
附图说明
图1为EM系列方法应用于混合物样品(全部都是未知物)的分析流程图。
图2为EM系列方法应用于含有标的物/已知物的混合物样品的分析流程图。
图3为实施例1中,航空煤油样品的TIC谱图,出峰时间为14~15分钟。
图4为实施例1中,航空煤油样品,从出峰14~15分钟的TIC谱图,总共采样了176次数据的混合谱,集中显示的3-D质谱图。
图5-1为实施例1中,线型正十一烷烃的EI-MS纯谱谱图(从日本AIST数据库获得)。
图5-2为实施例1中,线型正十二烷烃的EI-MS纯谱谱图(从日本AIST数据库获得)。
图6-1为实施例1中,通过锁定峰m/z=57和170,运用EM方法计算得出的纯谱。
图6-2为实施例1中,通过锁定峰m/z=57和170,运用EM方法计算得出的纯谱对应的伪浓度。
图7-1为实施例1中,通过锁定m/z=128,运用EM方法计算得出的纯谱。
图7-2为实施例1中,通过锁定m/z=128,运用EM方法计算得出的纯谱对应的伪浓度。
图8为实施例2中,挥发物样品的TIC谱图,m/z=800-850秒。
图9为实施例2中,挥发性物质样品的3-D质谱图,TIC=826.22-831.707秒。
图10-1为实施例2中,当锁住m/z=91的峰,用EM计算出的重建纯谱。
图10-2为实施例2中,当锁住m/z=91的峰,用EM计算出的重建纯谱对应的伪浓度。
图11为实施例2中,锁住m/z=91的峰后,重建纯谱和NIST数据库中的标准谱图相匹配的结果;其中丰度正值部分为重建纯谱,丰度负值部分为NIST数据库中标准谱图。
图12-1为实施例2中,当锁住m/z=71的峰,用EM计算出的重建纯谱。
图12-2为实施例2中,当锁住m/z=71的峰,用EM计算出的重建纯谱对应的伪浓度。
图13为实施例2中,锁住m/z=71的峰后,重建纯谱和NIST数据库中的标准谱图相匹配的结果;其中丰度正值部分为重建纯谱,丰度负值部分为NIST数据库中标准谱图。图14-1为实施例2中,当锁住m/z=105的峰,用EM计算出的重建纯谱。
图14-2为实施例2中,当锁住m/z=105的峰,用EM计算出的重建纯谱对应的伪浓度。
图15为实施例2中,锁住m/z=105的峰后,重建纯谱和NIST数据库中的标准谱图相匹配的结果;其中丰度正值部分为重建纯谱,丰度负值部分为NIST数据库中标准谱图。
图16-1为实施例2中,当锁住m/z=57和85的峰,用EM计算出的重建纯谱。
图16-2为实施例2中,当锁住m/z=57和85的峰,用EM计算出的重建纯谱对应的伪浓度。
图17为实施例2中,锁住m/z=57和85的峰后,重建纯谱和NIST数据库中的标准谱图相匹配的结果;其中丰度正值部分为重建纯谱,丰度负值部分为NIST数据库中标准谱图。
图18为实施例2中,经过矫正的4个重建纯谱的浓度,和重建总TIC浓度,以及、实际TIC浓度的汇总比较图。
图19为值空间和优化方法示例:在一次优化中,BTEM/tBTEM熵最小算法只会找到最下面的一个纯谱;而MREM算法就会找到全部的多个谱图。
具体实施方式
以下实施例用于说明本发明,但不用来限制本发明的范围。
本发明所用的物料,均为可以从市场上购买得到的常规物料。本发明中未提及的其他操作步骤,均为本领域的常规操作。
以气相色谱-质谱连用方法(GC-MS)为例(流程图如图1所示,含有未知物的混合物样品的分析):
将一个混合物样品注入气相色谱(GC)仪,混合物中含有A、B、C……和Z等成分,由于各成分自身物理化学性质(极性、分子大小等)的不同,它们在色谱柱中移动的速度不一样,从而导致各个成分到达检测器(质谱)的时间先后不同,而出现不同的峰。图1中,峰1、2和3是纯成分(已经分离出的成分)的峰;而峰4是一个混合物的重叠峰,里面的成分没有被分离开。现有的分析方法,需将一个成分和其他成分分离开,然后才能定性(例如峰1、2和3)。对于没有分离开的成分(例如重叠峰4),现有常规方法就需要将同样的样品重新在不同条件下做分离,经过反复非常多的次数,才能将各个成分在不同的分离过程中分离开。对于不同的样品,由于成分不同,所以分析方法不能通用。现有的分析方法,大部分的时间和精力都花在如何分离上,而且需要复杂的样品前处理过程,以减少样品中的成分数,从而降低色谱分离难度。
本发明的方法中,先用色谱技术(GC)分离混合物样品中的成分,经检测器(质谱)得到混合物的谱图(即初步的色谱图)。对于纯成分的峰(已经分离出的成分的峰,例如峰1、2和3),不再做任何处理。而对于重叠峰4这样的混合物峰,本发明就利用数据采样,即在不同的采样点进行数据采样(采样不一定是均匀的),从而得到一系列的混合谱(即混合物谱图);然后利用EM方法对混合谱进行多次操作,从而得到不同的重建纯谱及其对应的伪浓度。
将一个成分浓度比例固定的样品,利用色谱的分离能力,经过色谱技术分离后,在检测器端,就能得到一系列的含有不同成分、不同浓度比例的纯成分或者混合物的谱图。对于未分离好的混合成分的重叠峰,经过计算机多次采样(在不同时间点进行数据采样),得到的混合物谱图数目就会大于重叠峰所代表的混合物中含有的成分数。不同的化合物在色谱柱中的移动速度,不论色谱柱的分离能力高低,总有差异,所以在任何一个采样点,各个成分之间的比例是不一样的,这样就满足了EM系列方法应用的前提。
本发明的方法中,分离混合物的各成分主要是依赖EM系列方法重建出各成分的纯谱,所以对色谱技术中色谱柱的分离能力要求不高,但是好的色谱分离是欢迎的,因为这样能够减少混合物重叠峰,从而减少需要EM系列方法数学处理的量。对于不同的混合物样品,它们中的成分千差万别,同一种色谱分离方法将导致不同的分离结果,这是现有分离技术的问题,也是现有分离技术对不同样品不通用的原因;但是本发明的技术由于只是对混合物重叠峰操作,所以不同混合物样品的不同分离结果,对本发明的方法没有进行区分;也就是说,本发明的方法只要对不同的混合物重叠峰处理就行,而不在乎它们出现在什么地方。所以本发明所提供的方法,是一个通用的分析方法,对样品的成分组成情况不敏感;同时本发明的方法也可以大大降低样品前处理的过程;因为未经过前处理的样品,经色谱分离,会多一些混合物重叠峰,只是增加计算量,本发明方法对此没有任何敏感性。
利用本发明的技术,从混合谱(混合物谱图)中能够重建得到纯谱,本发明的方法并不能知道这些纯谱代表的是什么物质。在有可能的情况下,例如在运用气相色谱-电子电离质谱方法中(GC-EI-MS),由于EI-MS有庞大的成分数据库(例如NIST质谱数据库),所以将本发明方法计算得到的纯谱谱图和该成分数据库以及其他的成分数据库进行比对,就能容易地知道的重建出的纯谱代表的物质是什么了。
利用本发明的技术,对于混合物样品中的标的物/已知物分析(即含有标的物/已知物的混合物样品的分析)也非常有用,只需要利用标的物/已知物的标准谱图,设计EM系列方法的应用方法(例如锁住标的物/已知物的特征峰或特征谱图区间等),在没有分离好的谱图上,就可以计算出是否含有标的物/已知物和该标的物/已知物的伪浓度;在多次的重建中,如果在某一次结果中,重建纯谱和标的物/已知物的标准谱图是一致的,而且该重建纯谱的伪浓度是有意义的,则可以判断该混合物中含有该标的物,反之则不含有。这样可以大大降低对分离的要求,加快速度;该应用对于质检、检测和监测行业都非常有用。以气相色谱-质谱连用方法(GC-MS)为例,流程图如图2所示的为含有标的物/已知物的混合物样品的分析。
本发明的方法,将EM系列方法应用在标的物/已知物分析中,大大降低了对分离的要求,由于不需要将标的物/已知物和其他成分分开,所以对于不同的混合物样品,也就变得不敏感了;即不再需要设定各种不同的分析程序,所以本发明的方法是一个快速、通用的方法。
一、本发明中所使用的熵最小算法(BTEM,tBTEM和MREM)解释如下:
当一个样品,经过GC-MS后,经计算机采样得到一套数据Ak×υ,其中k是计算机采样次数,例如从出峰时间10-12分钟之间,每隔0.1秒采样,共采样21次。υ是数据通道数,例如被采样质谱的质量电荷比(m/z)是10-100,间隔是1个m/z,共91个数据通道。
下面,将算法的标准程序按照先后次序列出:
1)将数据集Ak×υ按照方程(1)进行奇异值分解(SVD),然后将右奇异矩阵VT υ×υ没有物理意义的部分除去,将对角矩阵Sk×υ中的全零部分去除,得到VT k×υ和Sk×k,矩阵U没有用到。
2)对矩阵VT k×υ中的每个行向量进行观察,那些看起来像背景噪声的行向量可以去除,例如从第j行向量以后看起来都是噪声,我们就把(j+1)-k行向量都去掉,得到VT j×υ.
3)从那些保留下来的VT j×υ中,通过检视这些j个行向量,使用者可能对某一个或几个峰,或者是某一区间有兴趣;例如峰m/z=91,或者m/z=90-100感兴趣;使用者就可以将这些峰或区间锁定(target),并开始运算。
4)首先,一个随机向量T1×j由计算机生成,该T向量将由所使用的优化算法自动更新,直到找到最优值。通过T,一个和纯谱有关的向量aest通过方程(2)得到。aest将根据T的变化而变化,直到优化结束方程;最后得到的aest我们认为就是纯谱。方程(2)中的S矩阵(见步骤1)可用,也可不用,根据实际需要而定。
5)将aest用锁定的峰或区间中的最高峰,按照方程(3)单位化。单位化后的向量,我们叫。这个过程就是我们的“锁峰”过程。其中,a′和a′′是锁定的峰和区间。
6)利用的得到的ahat构建优化方法的目标函数和罚函数,方程(4)是总的目标和法函数,具体目标函数和罚函数,根据不同类型的谱图,有不同的表达方式,详细请参考如下相关文献(文献①、②和③)。其中P函数是罚函数。通过方程(5),计算得到伪浓度信息。
7)检测目标函数的值,看是否已经达到优化终点;如果达到了,那么就输出纯谱aest,和伪浓度CkX1。如果没有达到优化终点,那么就将由优化方法,重新生成一个T1Xj,返回到步骤4)到7)重新计算,直到优化结束。
重复步骤3)到7),就可以锁定不同的峰或区间,运行算法,得到不同的纯谱和相关浓度,然后就可得到混合谱中的其他,或者是全部纯谱信息。
上述文献为:
①Chew,W.,E.Widjaja,andM.Garland,Band-targetentropyminimization(BTEM):Anadvancedmethodforrecoveringunknownpurecomponentspectra.applicationtotheFTIRspectraofunstableorganometallicmixtures.Organometallics,2002.21(9):p.1982-1990.
②Zhang,H.J.,etal.,Weightedtwo-BandTargetEntropyMinimizationforthereconstructionofpurecomponentmassspectra:Simulationstudiesandtheapplicationtorealsystems.JournaloftheAmericanSocietyforMassSpectrometry,2003.14(11):p.1295-1305.
③Zhang,H.,W.Chew,andM.Garland,Themulti-reconstructionentropyminimizationmethod:Unsupervisedspectralreconstructionofpurecomponentsfrommixturespectra,withouttheuseofaPrioriinformation.AppliedSpectroscopy,2007.61(12):p.1366-1372.
二、优化方法的使用:
熵最小系列算法(BTEM,tBTEM和MREM),都需要用到优化方法(得到T)。
首先,熵最小方法利用SVD后的VT数据,利用优化目标方程(4),构建一个多维值空间;然后利用优化方法,在这个值空间里搜索纯谱。对于锁住的不同峰或区间,其构建的多维值空间不同。
如图19所示,在搜索方式上,BTEM/tBTEM认为全局最小点是纯谱的位置,所以用的是全局优化方法,每次优化,只能得到一个纯谱。后来张华俊发现(上述文献③),在值空间中的各个局部最小点,都对应一个纯谱,所以研发了MREM;MREM使用局部优化方法,每次优化,在每个局部最小点,都有一个纯谱;所以一次优化,能够得到多个纯谱,所以速度大大加快,能力大大加强。将BTEM,tBTEM,MREM方法混合使用,效果更强;能够找到100%被其他谱图淹没的谱图,而且能够应用于紫外这样的高难度体系。
实施例1:在气质联用分析仪上应用本发明的方法分析航空煤油
一、实验条件
一个航空煤油样品(市场上购得),在安捷伦的气相色谱质谱联用分析仪上分析,色谱柱型号是HP5-MS,温度控制是先30℃恒温5分钟,然后从30℃升温到300℃,升温速度是10℃/分钟。由于航空煤油是个复杂的混合物,里面主要是很多结构类似的烷烃和芳香烃成分,还有很多成分属于同分异构体,所以混合程度非常复杂;即该航空煤油体系中各个成分用气相色谱质谱联用分析仪无法完全分离开。
该混合物经过上述仪器分析,得到TIC(TotalIonCurrent)谱图,其中从出峰14~15分钟的TIC谱图见图3。对这个时间段的谱图进行数据采样,总共采样了176次数据的混合谱,得到的一系列混合物谱图见图4(显示在一张3-D质谱图中)。其中质谱数据点是从50~200m/z,间隔为1m/z(每次采样时间是指一个特定的TIC时间,但是每次采样后的质谱,数据点就是50~200m/z,间隔为1m/z)。
二、快速发现航空煤油混合物中的已知成分
航空煤油样品中,含有很多烷烃。从日本的AIST(NationalInstituteofAdvancedIndustrialScienceandTechnology,Japan.)有机物数据库中找到了正十一烷烃的EI-MS标准谱图(见图5-1)。其中正十一烷的最高峰在m/z=57处,分子量峰在156m/z处。另外一个烷烃,正十二烷的最高峰也是在m/z=57处,而分子量峰在170m/z处。
对于采样好的混合物数据,应用EM算法,希望了解在TIC=14~15min的区间是否有正十一烷。EM算法将峰m/z=57单独锁定(具体用BTEM方法),或者将峰m/z=57和156一起锁定(具体用tBTEM方法),每种锁定方法都经过5次反复运算。结果发现,应用两种算法,每次计算出的重建纯谱,与该成分的标准谱图进行比对,没有找到相一致的图谱。所以认为,在TIC=14~15分钟之间,没有正十一烷烃。这个结果和实际实验结果吻合,因为用纯的正十一烷,在相同实验条件下做标定,该成分实际出峰时间在12.75分钟左右。
对于TIC=14~15min的区间,是否含有正十二烷,也利用上述EM方法试图找出。将峰m/z=57单独锁定,或者将峰m/z=57和170一起锁定,每种锁定方法都经过5次反复运算;每次计算出的重建纯谱和正十二烷的标准谱图(图5-2)相符。
其中锁定峰m/z=57和170的计算出的纯谱见图6-1,计算出的伪浓度峰显示(见图6-2),正十二烷烃出峰时间在14.37分钟左右。从上述图3实验结果中可以看出,TIC=14.33~14.43分钟之间,有一个很明显的混合物的重叠峰,这个重叠峰至少包含有两个成分,而且该峰重合度太大了,已经分辨不出其中的纯谱。而根据上述EM的计算结果,运用本发明方法则顺利的将正十二烷的纯谱重建出来,而且同时得到该正十二烷的伪浓度。用纯的正十二烷在相同的实验条件下做标定,其出峰时间就在14.38分钟左右。和EM计算的结果吻合。
通过该实施例,若想知道混合物样品中是否含有已知成分或标的物,对其进行分析,只需要根据已知物/标的物的标准谱图,在没有分离好的混合物谱图上,通过上述方法就能很快得出是否有该已知物/标的物;如果有的话,可以计算得出它们的浓度(需要根据浓度标准曲线来标定)和出峰时间是多少。
因此,本发明可用于快速分析,对于质检和质监方面非常有用;样品不需要分离好,不同的样品,分析方法类似;也不需要复杂的样品前处理,所以本发明效率非常高。
三、快速发现混合物中的未知成分
经过查看实验得出的TIC峰(图3),发现在TIC=14.25min左右有个峰,但是该峰和左右的几个峰重叠了,无法得到该峰所代表的纯成分的峰。在该TIC=14.25min的地方,采样得到的质谱混合谱上(采样方法同上),观察到有个m/z=128的峰。利用EM方法,锁定峰m/z=128,重建纯谱和其伪浓度如图7-1和图7-2所示。
但是,如果没有和数据库中纯成分的标准谱图进行比对,无法知道该重建纯谱代表什么化合物。当将重建纯谱和NIST标准数据库中的图谱匹配,确定了该重建纯谱和芳香烃萘的标准质谱图谱一样,从而确定在TIC=14.25min的峰是化合物萘的峰。在同样的条件下,用纯萘做标定后,萘的TIC出峰时间和本实施例的伪浓度峰出峰时间相吻合。
本实施例中上述快速发现混合物中的未知成分的实验说明,对于混合物样品中的未知物,并不需要很好的色谱分离,就能在未知是什么成分的情况下,仅仅根据混合谱中给出的信息(例如质谱峰),就可以快速找出未知成分;进一步若通过和标准谱图数据库进行比对,就能知道各成分的结构信息了。
实施例2:在气质联用分析仪上应用本发明的方法分析挥发性物质
一、实验条件:
一个海藻油经加氢催化得到的航空煤油样品(通过常规方法实验得到),在安捷伦的气相色谱质谱联用分析仪上分析,色谱柱型号是HP5-MS,温度控制是先30℃恒温5分钟,然后从30℃升温到300℃,升温速度是10℃/分钟。
该混合物经过上述仪器分析,得到TIC(TotalIonCurrent)谱图,其中从出峰800~850秒的TIC谱图见图8。图8中可见,在827~832秒的区间有峰,但是峰型不对称,不知这个峰是否是纯成分的峰。于是对827~832秒的谱图进行数据采样,总共采样了17次数据,得到的一系列混合物谱图见图9(显示在一张3-D质谱图中)。其中质谱数据点是从50~150m/z,间隔为0.5m/z。
二、熵最小算法应用:发现未知成分并定性
锁峰m/z=91:
利用熵最小算法(具体为BTEM算法),将峰m/z=91锁住后进行计算,得到的重建纯谱和其对应的伪浓度如图10-1和图10-2所示。可以看出,重建纯谱的伪浓度,有比较好的对称性。
重建纯谱的数据,用MREM软件导出成*.msp格式后,导入Agilent(安捷伦)的气相质谱,和仪器自带的NIST数据库中的标准谱图比对后,结果如图11所示。匹配结果显示重建纯谱代表的是化合物为C10H12,结构式为:
锁峰m/z=71:
利用熵最小算法(具体为BTEM算法),将峰m/z=71锁住后进行计算,得到的重建纯谱和其对应的伪浓度如图12-1和图12-2所示。可以看出,重建纯谱的伪浓度,有不错的对称性。
重建纯谱的数据,用MREM软件导出成*.msp格式后,导入Agilent的气相质谱,和仪器自带的NIST数据中的标准谱图比对后,结果如图13所示。匹配结果显示重建纯谱代表的化合物是C5H11Br,结构式为:
锁峰m/z=105:
利用熵最小算法(具体为BTEM方法),将峰m/z=105锁住后进行计算,得到的重建纯谱和其对应的伪浓度如图14-1和图14-2所示。可以看出,重建纯谱的伪浓度分布有点奇怪(不太对称),但是也是有可能的。
重建纯谱的数据,用MREM软件导出成*.msp格式后,导入Agilent的气相质谱,和仪器自带的NIST数据中的标准谱图比对后,结果如图15所示。匹配结果显示重建纯谱代表的化合物是C11H16,结构式为:
锁峰m/z=57和85:
利用熵最小算法(具体为tBTEM算法),将峰m/z=57和85锁住后进行计算,得到的重建纯谱和对应的伪浓度如图16-1和图16-2所示。
重建纯谱的数据,用MREM软件导出成*.msp格式后,导入Agilent的气相质谱,和仪器自带的NIST数据中的标准谱图比对后,结果如图17所示。匹配结果显示重建纯谱代表的化合物是C11H26,结构式为:
三、定量分析:
将上述四个重建纯谱的伪浓度经过矫正后,然后将它们的重建TIC浓度和实验得出的实际的TIC浓度比较(如图18所示)。结果显示,重建总TIC浓度和实际的TIC浓度所代表的峰的形状比较一致。它们的面积数值如下表:
在图18中,各曲线的含义分别表示为:粗实线代表实际(或原始)TIC浓度;粗虚线表示4个重建浓度的简单加和,用于和实际TIC浓度做比较;线1代表锁峰m/z=57和85后经过矫正的重建纯谱的TIC浓度;线2代表锁峰m/z=91后经过矫正的重建纯谱的浓度;线3代表锁峰m/z=71后经过矫正的重建纯谱的浓度;线4代表锁峰m/z=105后经过矫正的重建纯谱的浓度。
四、结果讨论:
利用化学分析方法,发现混合物中的未知成分,一直是比较困难而繁琐的事情。本实施例说明,只利用混合物谱图中的质谱峰(例如m/z=105),然后利用熵最小算法,就能方便地发现未知成分,得到未知成分的重建纯谱;然后利用数据库中的标准谱图进行比对,就能知道这些未知成分的纯谱所代表的物质。本实施例也显示,利用熵最小方法,能够很好的重建各成分的TIC浓度;因此就能够用常用的气相质谱方法定量各个成分的浓度。
实施例3:MREM算法在气质联用分析仪上应用本发明的方法分析挥发性物质
在和实施例2中相同的数据上,应用MREM方法,得到和实施例2中相同的结果。
MREM方法应用参数:
锁峰区间:整个区间,m/z=50~150(也就是没有指定是否用哪个特定的峰,不必指定特定的峰,让MREM方法自己找)。锁峰个数:每次只锁定一个。
优化方法:用于局部优化的模拟退火方法(上述文献③)。
优化次数:30次。
当用MREM算法,在实施例2的数据中,用如上参数后,能够得到30个重建纯谱,将其中的相同的重建纯谱去除后,得到和实施例2一样的结果,也就是4个纯谱,得到的TIC结果也和图18一样。实施例4:MREM和tBTEM算法连用,在气质联用分析仪上应用本发明的方法分析挥发性物质
在和实施例2中相同的数据上,连用MREM和tBTEM方法其中MREM算法的用处是:用其锁住区间的功能,和局部优化的功能;tBTEM算法的用处是:用其可锁住多重区间的功能。得到和实施例2中相同的结果。
MREM方法应用参数:
锁峰区间:两个区间,m/z=50~100,和m/z=101~150(也就是利用了tBTEM的锁住多重峰的功能,利用了MREM方法的锁住区间的功能(不必指定特定的峰))。
锁峰个数:每个区间一个,共两个峰(这个就是tBTEM的特点)。
优化方法:用于局部优化的模拟退火方法(上述文献③)。
优化次数:30次。
当连用MREM和tBTEM算法,在实施例2的数据中,用如上参数后,能够得到30个重建纯谱(重建纯谱包括纯的和复合谱),将其中相同的重建纯谱去除,并去除了复合谱后(详细算法参考上述文献②),得到和实施例2一样的结果,也就是4个纯谱,得到的TIC结果也和图18一样。
虽然,上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (11)
1.一种分析混合物成分的方法,包括:
(1)将恒定浓度比例的混合物样品,利用色谱技术进行分离,得到初步的色谱图;
(2)对步骤(1)得到的初步的色谱图中的任意一个区间,根据需要进行采样,得到一系列在不同保留时间的谱图,称之为混合物谱图;
(3)利用熵最小算法的系列方法,从步骤(2)得到的混合物谱图中,计算得到各重建纯谱及其对应的伪浓度;
所述伪浓度为利用熵最小算法的系列方法,在得到重建纯谱的时候,同时会计算得到一个和该重建纯谱相关的浓度;
对于含有未知物的混合物样品的分析,需要尽可能知道步骤(2)得到的混合物谱图中的成分信息的时候,利用熵最小算法的系列方法,在不同的参数下,反复多次计算后,得到尽可能多的重建纯谱;
所述不同的参数采用锁定不同的峰,使用不同的目标函数和罚函数或使用不同的优化方法。
2.根据权利要求1所述的分析混合物成分的方法,其特征在于,步骤(1)中,得到的初步的色谱图中,包括纯成分的峰和混合物的重叠峰。
3.根据权利要求2所述的分析混合物成分的方法,其特征在于,对于混合物的重叠峰,利用步骤(2)和步骤(3)进行分析。
4.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,步骤(1)得到的初步的色谱图中的任意一个区间,是指初步的色谱图的一段,或者全部。
5.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,步骤(1)中,色谱技术为:气相色谱技术、液相色谱技术或者两者的组合。
6.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,步骤(1)中,利用色谱技术进行分离,经检测器得到混合物初步的色谱图。
7.根据权利要求6所述的分析混合物成分的方法,其特征在于,所述检测器选自:质谱、红外、紫外、荧光、分光光度仪和核磁共振检测器中的一种或几种的组合。
8.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,步骤(3)中,熵最小算法的系列方法包括BTEM、tBTEM或MREM方法中的一种或几种的组合。
9.根据权利要求8所述的分析混合物成分的方法,其特征在于,利用BTEM或者tBTEM方法,每次运算,可以得到一个重建纯谱;利用MREM方法,每次运算,可以得到多个重建纯谱。
10.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,步骤(3)中,对于含有标的物/已知物的混合物样品的分析,基于标的物/已知物的标准谱图,利用熵最小算法的系列方法计算,构建使用不同参数的重建方法,得到该标的物/已知物的重建纯谱;在多次的重建中,如果重建纯谱和标的物/已知物的标准谱图相比较是一致的,而且该重建纯谱的伪浓度是有意义的,则该混合物中含有该标的物/已知物,反之则不含有。
11.根据权利要求1~3任意一项所述的分析混合物成分的方法,其特征在于,在步骤(3)后将得到的重建纯谱和标准数据库中的谱图进行比对,确认它们的信息;若比对后重建纯谱和标准数据库中的某谱图一致,则说明重建纯谱代表的是标准数据库中的某谱图所代表的成分。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410113885.6A CN103940934B (zh) | 2014-03-25 | 2014-03-25 | 一种分析混合物成分的方法 |
US15/120,974 US20170059537A1 (en) | 2014-03-25 | 2015-02-12 | Method for Analyzing Mixture Components |
SG11201607851PA SG11201607851PA (en) | 2014-03-25 | 2015-02-12 | Method for analyzing mixture components |
PCT/CN2015/072818 WO2015143963A1 (zh) | 2014-03-25 | 2015-02-12 | 一种分析混合物成分的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410113885.6A CN103940934B (zh) | 2014-03-25 | 2014-03-25 | 一种分析混合物成分的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940934A CN103940934A (zh) | 2014-07-23 |
CN103940934B true CN103940934B (zh) | 2016-04-27 |
Family
ID=51188688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410113885.6A Active CN103940934B (zh) | 2014-03-25 | 2014-03-25 | 一种分析混合物成分的方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20170059537A1 (zh) |
CN (1) | CN103940934B (zh) |
SG (1) | SG11201607851PA (zh) |
WO (1) | WO2015143963A1 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104267108A (zh) * | 2014-09-23 | 2015-01-07 | 广西壮族自治区药用植物园 | 一种基于mrem化学计量法与色谱联用的中药成分分析方法 |
CN105842330B (zh) * | 2015-09-09 | 2018-06-05 | 张华俊 | 一种质谱检测及解析方法 |
CN107389779B (zh) * | 2016-11-01 | 2020-05-12 | 苏州芷宁信息科技有限公司 | 一种化学结构的快速分析方法 |
CN106525890B (zh) * | 2016-11-29 | 2017-10-10 | 中国地质调查局油气资源调查中心 | 一种合成与分解混合气体水合物的实验装置 |
CN107727602B (zh) * | 2017-10-09 | 2020-04-07 | 广西科技大学 | 中红外光谱结合向量夹角定量分析三氯蔗糖含量的方法 |
CN108375639B (zh) * | 2018-02-27 | 2021-01-29 | 苏州芷宁信息科技有限公司 | 一种快速建立样品中组分质谱数据库的方法 |
CN108663437B (zh) * | 2018-06-19 | 2020-11-27 | 苏州芷宁信息科技有限公司 | 一种快速筛查标的物的方法 |
SG10202002459UA (en) * | 2020-03-17 | 2021-10-28 | Shimadzu Corp | Peak analyzing method and waveform processing device |
CN112947332B (zh) * | 2021-02-04 | 2022-12-16 | 威高国科质谱医疗科技(天津)有限公司 | 一种基于模拟退火的三重四极质谱仪参数优化方法 |
CN113945530B (zh) * | 2021-10-19 | 2024-06-14 | 中国计量科学研究院 | 气体浓度检测方法和质谱仪 |
CN115248271A (zh) * | 2022-06-07 | 2022-10-28 | 浙江理工大学绍兴生物医药研究院有限公司 | 鼠尾草挥发性成分定性分析方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7972770B2 (en) * | 2005-08-16 | 2011-07-05 | Memorial Sloan-Kettering Cancer Center | Methods of detection of cancer using peptide profiles |
WO2007112597A1 (en) * | 2006-04-05 | 2007-10-11 | National Research Council Of Canada | Blind extraction of pure component mass spectra from overlapping mass spectrometric peaks |
AU2009231804B2 (en) * | 2008-03-31 | 2015-06-04 | Boston Medical Center Corporation | Predictive marker for topoisomerase I inhibitors |
EP2350926A2 (en) * | 2008-11-24 | 2011-08-03 | Institut Ruder Boskovic | Method of and system for blind extraction of more than two pure components out of spectroscopic or spectrometric measurements of only two mixtures by means of sparse component analysis |
US20120057743A1 (en) * | 2009-01-09 | 2012-03-08 | Chemimage Corporation | System and method for fluorescence guided ingredient specific particle sizing |
US20120041689A1 (en) * | 2009-01-09 | 2012-02-16 | Chemimage Corporation | System and method for particle detection in spectral domain |
WO2012091680A1 (en) * | 2010-12-30 | 2012-07-05 | Nanyang Technological University | A device for controlled release of a bioactive agent |
WO2012112119A1 (en) * | 2011-02-16 | 2012-08-23 | Agency For Science, Technology And Research (A*Star) | Method and apparatus for analyzing a liquid sample by raman spectroscopy |
CN104267108A (zh) * | 2014-09-23 | 2015-01-07 | 广西壮族自治区药用植物园 | 一种基于mrem化学计量法与色谱联用的中药成分分析方法 |
-
2014
- 2014-03-25 CN CN201410113885.6A patent/CN103940934B/zh active Active
-
2015
- 2015-02-12 SG SG11201607851PA patent/SG11201607851PA/en unknown
- 2015-02-12 US US15/120,974 patent/US20170059537A1/en not_active Abandoned
- 2015-02-12 WO PCT/CN2015/072818 patent/WO2015143963A1/zh active Application Filing
Non-Patent Citations (3)
Title |
---|
Application of the BTEM family of algorithms to reconstruct individual UV-vis spectra from multi-component mixtures;Feng Gao等;《Chemometrics and Intelligent Laboratory Systems》;20090115;第95卷(第1期);第94~100页 * |
The multi-reconstruction entropy minimization method: unsupervised spectral reconstruction of pure components from mixture spectra, without the use of a Priori information;Huajun Zhang等;《Applied Spectroscopy》;20071231;第61卷(第12期);第1366~1372页 * |
Weighted two-band target entropy minimization for the reconstruction of pure component mass spectra: simulation studies and the application to real systems;Huajun Zhang等;《Journal of the American Society for Mass Spectrometry》;20031130;第14卷(第11期);第1295~1305页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2015143963A1 (zh) | 2015-10-01 |
CN103940934A (zh) | 2014-07-23 |
SG11201607851PA (en) | 2016-11-29 |
US20170059537A1 (en) | 2017-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103940934B (zh) | 一种分析混合物成分的方法 | |
US10302607B2 (en) | Method for detailed and bulk classification analysis of complex samples using vacuum ultra-violet spectroscopy and gas chromatography | |
CN103913534B (zh) | 一种天然气中系列烃类化合物碳同位素分析方法 | |
Martin-Drumel et al. | Automated microwave double resonance spectroscopy: A tool to identify and characterize chemical compounds | |
Pierce et al. | Predicting percent composition of blends of biodiesel and conventional diesel using gas chromatography–mass spectrometry, comprehensive two-dimensional gas chromatography–mass spectrometry, and partial least squares analysis | |
JP6557737B2 (ja) | 質量スペクトルの検出及び解析方法 | |
CN103235057A (zh) | 一种利用气相色谱-质谱不解析化合物鉴别白酒原产地的方法 | |
Rüger et al. | Comprehensive chemical comparison of fuel composition and aerosol particles emitted from a ship diesel engine by gas chromatography atmospheric pressure chemical ionisation ultra-high resolution mass spectrometry with improved data processing routines | |
Boegelsack et al. | Development of retention time indices for comprehensive multidimensional gas chromatography and application to ignitable liquid residue mapping in wildfire investigations | |
Krishnamurthy | Complete reduction to amplitude frequency table (CRAFT)—a perspective | |
Fernando et al. | Visible/infrared dissociation of NO3: Roaming in the dark or roaming on the ground? | |
CN104267108A (zh) | 一种基于mrem化学计量法与色谱联用的中药成分分析方法 | |
Wei et al. | Authentication of chaste honey adulterated with high fructose corn syrup by HS-SPME-GC-MS coupled with chemometrics | |
CN115004019A (zh) | 一种基于晶体结构数据库的材料解析方法、***及应用 | |
CN102798683B (zh) | 一种通用的气相色谱-质谱全组分定量分析方法 | |
Wilkins | Multidimensional GC for qualitative IR and MS of mixtures | |
CN113933442B (zh) | 全二维气相色谱-质谱联用数据分析方法、***及应用 | |
JP2024526079A (ja) | 質量スペクトル中の分子種を同定するための方法及び装置 | |
Laws et al. | Velocity map imaging spectroscopy of C2H− and C2D−: A benchmark study of vibronic coupling interactions | |
CN114577966B (zh) | 一种mscc结合调制峰归类的gc×gc指纹快速比较方法 | |
CN107219321B (zh) | 一种混合质谱筛除方法 | |
CN109239247A (zh) | 一种液相色谱无对照品定性分析方法 | |
CN104849392B (zh) | 相对书写时间鉴定方法 | |
US20230288384A1 (en) | Method for determining small molecule components of a complex mixture, and associated apparatus and computer program product | |
CN112986307B (zh) | 核磁共振浓度排序谱的构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200814 Address after: Room 408, 4 / F, building B9, No. 11, Kaiyuan Avenue, Huangpu District, Guangzhou City, Guangdong Province 510000 (R & D office only) Patentee after: Guangzhou molecular information technology Co., Ltd Address before: 215600 hi tech Innovation Service Center, 1 Cathay Pacific Road, Zhangjiagang, Jiangsu, Suzhou A-326 Patentee before: Zhang Huajun |
|
TR01 | Transfer of patent right |