CN113138248B - 特征图谱提取方法及料液配制质量稳定性的检测方法 - Google Patents

特征图谱提取方法及料液配制质量稳定性的检测方法 Download PDF

Info

Publication number
CN113138248B
CN113138248B CN202110446263.5A CN202110446263A CN113138248B CN 113138248 B CN113138248 B CN 113138248B CN 202110446263 A CN202110446263 A CN 202110446263A CN 113138248 B CN113138248 B CN 113138248B
Authority
CN
China
Prior art keywords
characteristic
spectrum
ion
feed liquid
data matrix
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
Application number
CN202110446263.5A
Other languages
English (en)
Other versions
CN113138248A (zh
Inventor
胡芸
金晶
李永青
杨敬国
阮艺斌
彭黔荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Tobacco Guizhou Industrial Co Ltd
Original Assignee
China Tobacco Guizhou Industrial Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Tobacco Guizhou Industrial Co Ltd filed Critical China Tobacco Guizhou Industrial Co Ltd
Priority to CN202110446263.5A priority Critical patent/CN113138248B/zh
Publication of CN113138248A publication Critical patent/CN113138248A/zh
Application granted granted Critical
Publication of CN113138248B publication Critical patent/CN113138248B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/8686—Fingerprinting, e.g. without prior knowledge of the sample components
    • 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/62—Detectors specially adapted therefor
    • G01N30/72—Mass spectrometers
    • 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
    • G01N2030/022—Column chromatography characterised by the kind of separation mechanism
    • G01N2030/025—Gas chromatography

Landscapes

  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Library & Information Science (AREA)
  • Engineering & Computer Science (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种特征图谱提取方法，包括：获取气相离子迁移谱图的数据矩阵；在数据矩阵中选取与特征区域相关的特征数据矩阵，特征数据矩阵中位于同一行的数据表征在同一气相色谱时间但不同离子漂移时间的各离子强度，特征数据矩阵中位于同一列的数据表征在同一离子漂移时间但不同气相色谱时间的各离子强度；根据特征数据矩阵中各行的数据，构建色谱特征图谱；根据特征数据矩阵中各列的数据，构建迁移谱特征图谱；基于色谱特征图谱和迁移谱特征图谱构建特征图谱。该方法能够快速、简便地从气相离子迁移谱图中提取出特征图谱，简化了数据分析的运算量。本发明还提供了一种料液配制质量稳定性的检测方法。

Description

特征图谱提取方法及料液配制质量稳定性的检测方法
技术领域
本发明涉及气相离子迁移技术应用领域，特别涉及特征图谱提取方法及料液配制质量稳定性的检测方法。
背景技术
烟用料液可以调节烟气的酸碱度，能够改善烟支的吸味，提高抽吸舒适性，适当增加卷烟香气以及掩盖和去除烟叶本身带来的不良杂气和刺激性等，因此料液的种类和料液质量的稳定性会直接影响卷烟产品的口感以及质量。
现有技术中，为区分物质种类或检测物质的质量是否稳定一般采用对待检测物质的气相离子迁移谱图(也称GC-IMS谱图)中的特征区域进行分析来判断。例如专利申请号为201710948020.5的关于卷烟质量稳定性评判的一篇中国专利，即是通过获取不同批次卷烟样品的GC-IMS谱图，然后从这些谱图中中选取出被分析牌号卷烟的特征分析区域，而后对不同批次卷烟的特征分析区域进行分析，鉴别不同批次卷烟挥发性成分稳定性，从而判断卷烟的质量稳定性的。再比如专利申请号为CN201911403118.8的关于烟叶种类识别的一篇中国专利，也是利用对待测烤烟烟叶样品进行顶空分析和气相色谱-离子迁移谱分析，获得GC-IMS谱图，通过比对谱图中的特征分析区域来鉴别烤烟烟叶种类。在上述专利中，均是通过在各个GC-IMS谱图中查找对比各特征峰的变化以分析物质种类或稳定性的，然而GC-IMS谱图中的特征峰的信息为三维信息，这导致特征峰的查找比对过程繁琐，计算量大，且耗时较长。
发明内容
本发明的主要目的在于解决现有技术中关于在各GC-IMS谱图中查找比对特征峰的过程繁琐，计算量大，且耗时较长的问题。
为实现上述目的，本发明的实施例提供了一种特征图谱提取方法，能够快速、简便地从气相离子迁移谱图中提取出特征图谱，简化了比对和数据分析的运算量。具体地，该方法用于从气相离子迁移谱图中提取特征图谱，包括：
获取气相离子迁移谱图的数据矩阵；
在数据矩阵中选取与特征区域相关的特征数据矩阵，特征数据矩阵中位于同一行的数据表征在同一气相色谱时间但不同离子漂移时间的各离子强度，特征数据矩阵中位于同一列的数据表征在同一离子漂移时间但不同气相色谱时间的各离子强度；
根据特征数据矩阵中各行的数据，构建色谱特征图谱；
根据特征数据矩阵中各列的数据，构建迁移谱特征图谱；
基于色谱特征图谱和迁移谱特征图谱构建特征图谱。
采用本方案，可以在一定程度上简化数据分析的工作量，为利用GC-IMS谱图辨识物质种类以及物质稳定性提供了便利。并且该方法首先对待检测物质的气相色谱进行分离，然后再获得离子迁移谱，最大程度地表达样品的指纹特征，确保样品数据的全面性、真实性及溯源性。
作为本发明的一个具体实施方式，根据特征数据矩阵中各行的数据，构建色谱特征图谱包括：
对特征数据矩阵的每一行中的各个数据进行累加求和，获得各行相对应的行离子强度；
根据各行离子强度，构建色谱特征图谱；
根据特征数据矩阵中各列的数据，构建迁移谱特征图谱包括：
对特征数据矩阵的每一列中的各个数据进行累加求和，获得各列相对应的列离子强度；
根据各列离子强度，构建迁移谱特征图谱。
作为本发明的一个具体实施方式，色谱特征图谱和迁移谱特征图谱均为向量，其中，色谱特征图谱中的元素为按照气相色谱时间依次排列的行离子强度，迁移谱特征图谱中的元素为按照离子漂移时间依次排列的列离子强度。
作为本发明的一个具体实施方式，特征图谱为向量，特征图谱中的元素由色谱特征图谱中的各元素和迁移谱特征图谱中的各元素组成。
作为本发明的一个具体实施方式，特征区域为离子漂移时间为1.06ms～2.00ms，气相色谱时间为0～200s所对应的区域。
作为本发明的一个具体实施方式，上述特征图谱提取方法用于料液种类识别和/或料液配制质量稳定性检测。
相应地，本发明还提供了一种料液配制质量稳定性的检测方法，包括：
利用顶空-气相色谱-离子迁移谱技术获取不同批次料液样品的气相离子迁移谱图；
利用如权利要求1至5任一项的特征图谱提取方法对各待检测料液样品的气相离子迁移谱图进行分析，获得各料液样品的特征图谱；
计算各料液样品的特征图谱的模，并基于各特征图谱的模计算相对标准偏差；
根据相对标准偏差判断不同批次料液配制质量是否稳定。
作为本发明的一个具体实施方式，利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
顶空进样条件设定为：孵化温度50℃、孵化时间5min、进样体积500μl、进样针温度80℃。
作为本发明的一个具体实施方式，利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
气相色谱-离子迁移谱仪器条件设定为：色谱柱为多毛细管分离柱，固定相OV-5、色谱柱温40℃、漂移管温度45℃、进样口温度80℃、进样器-色谱柱管路温度65℃、载气为高纯氮气。
作为本发明的一个具体实施方式，利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
载气流速采用程序升速模式，其中，0～1min载气速度为2ml/min，1～3min载气速度由2ml/min逐渐上升到50ml/min，3～5min载气速度由50ml/min逐渐上升到150ml/min，5～20min载气速度维持在150ml/min；迁移谱漂移气体流量150ml/min。
附图说明
图1是本发明提供的特征图谱提取方法的流程图；
图2是本发明提供的不同牌号烟用料液的气相色谱-离子迁移谱对比图；
图3是本发明提供的基于GC-IMS谱图获得的气相色谱方向的色谱特征图谱；
图4是本发明提供的基于GC-IMS谱图获得的离子漂移方向的迁移谱特征图谱；
图5是本发明提供的不同牌号料液的特征图谱的主成分分析图；
图6是本发明提供的不同牌号不同批次料液的特征图谱的模的得数示意图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式，本领域技术人员可由本说明书所揭示的内容轻易地了解本发明的其他优点及功效。虽然本发明的描述将结合较佳实施例一起介绍，但这并不代表此发明的特征仅限于该实施方式。恰恰相反，结合实施方式作发明介绍的目的是为了覆盖基于本发明的权利要求而有可能延伸出的其它选择或改造。为了提供对本发明的深度了解，以下描述中将包含许多具体的细节。本发明也可以不使用这些细节实施。此外，为了避免混乱或模糊本发明的重点，有些具体细节将在描述中被省略。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
应注意的是，在本说明书中，相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明的实施方式作进一步地详细描述。
如图1所示，本发明提供了一种特征图谱提取方法，可以用于从气相离子迁移谱图中提取特征图谱，具体地，该方法包括如下步骤：
S1：获取气相离子迁移谱图的数据矩阵。
具体地，可以利用气相色谱-离子迁移谱仪器(比如德国GAS仪器公司的型号为FlavourSpec的顶空-气相色谱-离子迁移谱仪器)获取得到待检测物质(诸如料液、卷烟、植物油等等)的气相离子迁移谱图，然后利用软件导出气相离子迁移谱图的数据矩阵(也可称之为气相色谱-离子迁移谱数据矩阵)。需要说明的是，由气相离子迁移谱图获得的数据矩阵包括了离子漂移时间、气相色谱时间和离子强度等信息，为后续分析提供了更加丰富的化学信息。可选地，该数据矩阵可以为CSV格式。具体地，该数据矩阵的行方向表征了在某个色谱时间的迁移谱信息，该矩阵的列方向表示了在某个漂移时间点的色谱信息，或者也可以理解为该数据矩阵中的每个数据元素表征的是在某个色谱时间点和某一漂移时间点的离子强度，其中，同一行的数据表征的是在同一气相色谱时间但不同离子漂移时间的各离子强度，同一列的数据表征的是在同一离子漂移时间但不同气相色谱时间的各离子强度。具体地，数据矩阵的表现形式可以为：
Figure BDA0003037019070000041
其中，X代表气相离子迁移谱的数据矩阵，xij代表该矩阵中位于第i行第j列的数据元素，该数据元素表征了在某一特定气相色谱时间和某一特定离子漂移时间的离子强度。其中，第一行至第m行是按照气相色谱时间排列的，第一列至第n列是按照离子漂移时间排列的。
S2：在数据矩阵中选取与特征区域相关的特征数据矩阵，特征数据矩阵中位于同一行的数据表征在同一气相色谱时间但不同离子漂移时间的各离子强度，特征数据矩阵中位于同一列的数据表征在同一离子漂移时间但不同气相色谱时间的各离子强度。
具体地，如上面所提到的，数据矩阵中的各个数据元素是与气相色谱时间和离子漂移时间相关的离子强度。图2为不同牌号烟用料液的挥发性和半挥发性成分的气相离子迁移谱图，由图2可看出，气相离子迁移谱早期离子漂移时间1.0ms处的第一条竖线是水反应离子峰(即RIP峰)，当挥发性化合物存在于离子迁移的电离区域时，水反应离子峰的强度就会降低或消失。RIP峰的两侧的每一个点代表一种挥发性有机物，颜色浅表示浓度较小，颜色深表示浓度越大。并且特征峰数量的多少能表现挥发性成分的多少，挥发性组分的含量可以用峰强度表示，例如在图2中，牌号A料液的特征峰的个数较多，因此牌号A的挥发性成分较多，而牌号C的挥发性成分较少。从图2中还可以看出，不同牌号的料液的挥发性成分有一定差异，体现在图2中的特征峰的位置及强度不同，因此可以利用特征峰强度和其位置判断料液种类是否相同。其中特征峰的位置即特征峰所处的离子漂移时间和气相色谱时间。但图2中是三维信息，由图2中的谱图不能便捷得到特征峰的位置和强度信息，因此，申请人尝试将含有特征峰的特征区域内的数据提取出，构建成了一个新的特征数据矩阵，可选地，特征区域为离子漂移时间为1.06ms～2.00ms，气相色谱时间为0～200s所对应的区域，即特征数据矩阵中的各个数据元素为离子漂移时间为1.06ms～2.00ms，气相色谱时间为0～200s这一范围内所对应的离子强度。即从数据矩阵中提取出离子漂移时间介于1.06ms～2.00ms且气相色谱时间介于0～200s之间的各个离子强度，构建新的特征数据矩阵。具体地，该特征数据矩阵的形式可以为
Figure BDA0003037019070000051
其中，Y代表特征数据矩阵，yij代表处于第i行第j列的数据元素，同样的，第一行至第a行是按照气相色谱时间排列的，第一列至第b列是按照离子漂移时间排列的，且a小于或等于m，b小于或等于n。
S3：根据特征数据矩阵中各行的数据，构建色谱特征图谱。
示例性的，该色谱特征图谱的构建方法包括：
S31：对特征数据矩阵的每一行中的各个数据进行累加求和，获得各行相对应的行离子强度；
S32：根据各个行离子强度，构建气相色谱特征图谱。
以上面提到的特征数据矩阵Y为例，各行对应的行离子强度的计算公式为：
Figure BDA0003037019070000061
其中，qc代表第c行所对应的行离子强度。此时，色谱特征图谱可以由各个行所对应的行离子强度构造而成，其可以表示为：
q＝(q1,q2,q3,…qa)
其中，q代表色谱特征图谱，其为向量形式，该色谱特征图谱中的各个元素是按照气相色谱时间依次排列的各个行离子强度。
在具体实施时，可以根据计算得到的色谱特征图谱中的各个元素绘制为曲线，如图3所示，图中四条曲线分别为根据牌号A、牌号B、牌号C、牌号D四种料液的色谱特征图谱绘制的曲线，分别表示各牌号料液的行离子强度随气相色谱时间的变化曲线。
S4：根据特征数据矩阵中各列的数据，构建迁移谱特征图谱。
示例性的，迁移谱特征图谱的构建方法包括：
S41：对特征数据矩阵的每一列中的各个数据进行累加求和，获得各列相对应的列离子强度；
S42：根据各列离子强度，构建迁移谱特征图谱。
仍以上面提到的特征数据矩阵Y为例，各列对应的列离子强度的计算公式为：
Figure BDA0003037019070000062
其中，ld代表第d列所对应的列离子强度。此时，迁移谱特征图谱可以由各个列所对应的列离子强度构造而成，其可以表示为：
l＝(l1,l2,l3,…lb)
其中，l代表迁移谱特征图谱，其为向量形式，该迁移谱特征图谱中的各个元素ld是按照离子漂移时间依次排列的各个列离子强度。
在具体实施时，可以根据计算得到的迁移谱特征图谱中的各个元素绘制曲线，如图4所示，图中四条曲线分别为根据牌号A、牌号B、牌号C、牌号D四种料液的迁移谱特征图谱绘制的曲线，分别表示各牌号料液的列离子强度随离子漂移时间的变化曲线。
S5：基于色谱特征图谱和迁移谱特征图谱构建特征图谱。
即，为了充分利用色谱特征图谱和迁移谱特征图谱，可以采用数据融合的方式将上述计算所得的色谱特征图谱和迁移谱特征图谱合并为一体进行分析。
示例性的，该特征图谱为向量，特征图谱中的元素由色谱特征图谱中的各元素和迁移谱特征图谱中的各元素组成。
以色谱特征图谱q，迁移谱特征图谱l为例，特征图谱即是把色谱特征图谱q和迁移谱特征图谱l这两个向量组合起来，具体地，特征图谱t可以表示为：
t＝(q1,q2,q3,…qa,l1,l2,l3,…lb)
当然，特征图谱也可以表示为：
t＝(l1,l2,l3，…lb,q1,q2,q3,…qa)
即在，特征图谱的构建过程中对色谱特征图谱和迁移谱特征图谱的排列顺序不做限制。
可选地，本发明提供的特征图谱提取方法还可以应用于料液种类的识别和/或料液配制质量稳定性的检测。
这样通过上述方法，可以将利用矩阵表征的GC-IMS谱图降维为利用向量表征的特征图谱，在一定程度上简化了数据分析的工作量，为利用GC-IMS谱图辨识物质种类以及物质稳定性提供了便利。并且该方法首先对待检测物质的气相色谱进行分离，然后再获得离子迁移谱，通过顶空-气相色谱-离子迁移谱对料液挥发性成分进行无损分析，最大程度地表达样品的指纹特征，确保样品数据的全面性、真实性及溯源性。
相应的，本发明还提供了一种料液配制质量稳定性的检测方法，包括如下步骤：
A1：利用顶空-气相色谱-离子迁移谱技术获取不同批次料液样品的气相离子迁移谱图；
A2：利用上述提到的特征图谱提取方法对各待检测料液样品的气相离子迁移谱图进行分析，获得各料液样品的特征图谱；
A3：计算各料液样品的特征图谱的模，并基于各特征图谱的模计算相对标准偏差；
A4：根据相对标准偏差判断不同批次料液配制质量是否稳定。
因为计算得到的样品的特征图谱实际上是一个向量，因此这里的所说的特征图谱的模其实可以理解为求解向量的大小，可以利用公式M＝||t||2计算特征图谱的模，其中，M表示特征图谱的模，t表示特征图谱。
具体地，特征图谱的模可以在某种程度上表征离子的强度值，因此可以利用各个特征图谱的模计算他们之间的相对标准偏差，若该相对标准偏差小于阈值相对标准偏差(比如可以设置为5％)，则认为该料液的配制质量稳定，否则，该料液的配制质量不稳定。
示例性的，利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
顶空进样条件设定为：孵化温度50℃、孵化时间5min、进样体积500μl、进样针温度80；
气相色谱-离子迁移谱仪器条件设定为：色谱柱为多毛细管分离柱，固定相OV-5、色谱柱温40℃、漂移管温度45℃、进样口温度80℃、进样器-色谱柱管路温度65℃、载气为高纯氮气；
载气流速采用程序升速模式，其中，0～1min载气速度为2ml/min，1～3min载气速度由2ml/min逐渐上升到50ml/min，3～5min载气速度由50ml/min逐渐上升到150ml/min，5～20min载气速度维持在150ml/min；迁移谱漂移气体流量150ml/min。
在该条件下，料液的挥发性物质成分分离度好，检测到的挥发性成分多，能快速采集到料液中的挥发性成分。
相应的，本发明还提供了一种料液种类的识别方法，包括以下步骤：
获取各待检测料液的气相离子迁移谱图；
利用上述提到的特征图谱提取方法对获取的各气相离子迁移谱图进行分析，获得各待检测料液的特征图谱；
采用主成分分析方法对特征图谱进行可视化分析，对各待检测料液进行分类识别。
即，经过主成分分析后，得到主成分的得分矢量，选取前两个主成分，根据第一主成分的得分对第二主成分的得分，画出前述检测料液样品在二维平面上的分布情况，由图形可以直观地看出个样品的聚集情况。这样通过观察样品的离散情况，能够用特征图谱很好地区分不同牌号的料液。如果是同种类样品，样品聚集在一起，反之，样品则离得较远。
具体地，获取各待检测料液的气相离子迁移谱图的实验条件可以为：
顶空进样条件设定为：孵化温度50℃、孵化时间5min、进样体积500μl、进样针温度80；
气相色谱-离子迁移谱仪器条件设定为：色谱柱为多毛细管分离柱，固定相OV-5、色谱柱温40℃、漂移管温度45℃、进样口温度80℃、进样器-色谱柱管路温度65℃、载气为高纯氮气；
载气流速采用程序升速模式，其中，0～1min载气速度为2ml/min，1～3min载气速度由2ml/min逐渐上升到50ml/min，3～5min载气速度由50ml/min逐渐上升到150ml/min，5～20min载气速度维持在150ml/min；迁移谱漂移气体流量150ml/min。
在该条件下，料液的挥发性物质成分分离度好，检测到的挥发性成分多，能快速采集到料液中的挥发性成分。
实施例
(1)料液样品采集制备
收集配制后的烟用料液样品：从香糖料厨房收集4个不同牌号不同批次的料液成品合计51个，每个样品稀释100倍后，各取500ul，分别置于20mL顶空进样瓶中，压盖密封待测。
(2)获取各个料液样品的GC-IMS谱图
A、分析仪器：采用德国GAS仪器公司的顶空-气相色谱-离子迁移谱仪器(型号为FlavourSpec)；
B、实验条件：
顶空进样条件设定为：孵化温度50℃、孵化时间5min、进样体积500μl、进样针温度80；
气相色谱-离子迁移谱仪器条件设定为：色谱柱为多毛细管分离柱(MCC)，固定相OV-5、色谱柱温40℃、漂移管温度45℃、进样口温度80℃、进样器-色谱柱管路温度65℃、载气为高纯氮气；
载气流速采用程序升速模式，其中，0～1min载气速度为2ml/min，1～3min载气速度由2ml/min逐渐上升到50ml/min，3～5min载气速度由50ml/min逐渐上升到150ml/min，5～20min载气速度维持在150ml/min；迁移谱漂移气体流量150ml/min。
在上述实验条件下，利用上述仪器获得了各个料液样品的GC-IMS谱图，如图2所示。
(4)特征图谱提取，图3和图4分别表示的是根据色谱特征图谱绘制的曲线和根据迁移谱特征图谱绘制的曲线。
(5)对色谱特征图谱和迁移谱特征图谱融合后的特征图谱做主成分分析，识别烟用料液种类。图5为基于色谱特征图谱和迁移谱特征图谱融合后的不同牌号料液的特征图谱的主成分分析图。由该图可以看出，不同牌号的料液样品的主成分得分明显不同，其相同牌号样品的主成分得分聚集在一起，不同牌号样品的则相对远离，由图5中样品主成分得分的聚集程度，可以知道，这次检测的不同牌号不同批次的料液成品共有四种牌号，与实际情况相符，因此本发明提供的特征图谱提取方法可以成功应用于料液种类的识别，且数据分析过程被极大简化，提高了识别效率。
(6)计算色谱特征图谱和迁移谱特征图谱融合后的特征图谱的模，根据各特征图谱的模计算相对标准偏差。图6为基于色谱和迁移谱特征图谱融合后的不同牌号料液的特征图谱的模，通过计算得到牌号A、B、C、D不同批次的料液的相对标准偏差为3.13％，1.90％，1.15％和2.32％。以阈值误差为5％为例，由上述计算结果可知，牌号A、B、C、D的料液配制质量均稳定，但是牌号A不同批次间的料液的质量变化相对较大，而牌号C不同批次间的料液的质量变化相对较小。这与实际情况也相符，因此本发明提供的特征图谱提取方法也可以成功应用于料液配制质量稳定性的检测评价，且利用该方法后，料液稳定性评价过程中的数据分析的运算量也被极大简化，提高了评价效率。
虽然通过参照本发明的实施方式，已经对本发明进行了图示和描述，但本领域的普通技术人员应该明白，以上内容是结合具体的实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。本领域技术人员可以在形式上和细节上对其作各种改变，包括在不偏离本发明的精神和范围的情况下做出若干简单推演或替换。

Claims (8)

1.一种特征图谱提取方法，用于从气相离子迁移谱图中提取特征图谱，其特征在于，包括：
获取所述气相离子迁移谱图的数据矩阵；
在所述数据矩阵中选取与特征区域相关的特征数据矩阵，所述特征数据矩阵中位于同一行的数据表征在同一气相色谱时间但不同离子漂移时间的各离子强度，所述特征数据矩阵中位于同一列的数据表征在同一离子漂移时间但不同气相色谱时间的各离子强度；
根据所述特征数据矩阵中各行的数据，构建色谱特征图谱；
根据所述特征数据矩阵中各列的数据，构建迁移谱特征图谱；
基于所述色谱特征图谱和所述迁移谱特征图谱构建所述特征图谱；
其中，所述色谱特征图谱和所述迁移谱特征图谱均为向量，所述色谱特征图谱中的元素为按照所述气相色谱时间依次排列的所述行离子强度，所述迁移谱特征图谱中的元素为按照所述离子漂移时间依次排列的所述列离子强度；
所述根据所述特征数据矩阵中各行的数据，构建色谱特征图谱包括：
对所述特征数据矩阵的每一行中的各个数据进行累加求和，获得各行相对应的行离子强度；
根据各所述行离子强度，构建色谱特征图谱；
所述根据所述特征数据矩阵中各列的数据，构建迁移谱特征图谱包括：
对所述特征数据矩阵的每一列中的各个数据进行累加求和，获得各列相对应的列离子强度；
根据各所述列离子强度，构建迁移谱特征图谱。
2.如权利要求1所述的特征图谱提取方法，其特征在于，所述特征图谱为向量，所述特征图谱中的元素由所述色谱特征图谱中的各元素和所述迁移谱特征图谱中的各元素组成。
3.如权利要求1所述的特征图谱提取方法，其特征在于，所述特征区域为所述离子漂移时间为1.06ms～2.00ms，所述气相色谱时间为0～200s所对应的区域。
4.如权利要求1至3任一项所述的特征图谱提取方法，其特征在于，用于料液种类识别和/或料液配制质量稳定性检测。
5.一种料液配制质量稳定性的检测方法，其特征在于，包括：
利用顶空-气相色谱-离子迁移谱技术获取不同批次料液样品的气相离子迁移谱图；
利用如权利要求1至3任一项所述的特征图谱提取方法对各待检测料液样品的所述气相离子迁移谱图进行分析，获得各所述料液样品的特征图谱；
计算各所述料液样品的特征图谱的模，并基于各所述特征图谱的模计算相对标准偏差；
根据所述相对标准偏差判断不同批次料液配制质量是否稳定。
6.如权利要求5所述的方法，其特征在于，所述利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
顶空进样条件设定为：孵化温度50℃、孵化时间5min、进样体积500μl、进样针温度80℃。
7.如权利要求5所述的方法，其特征在于，所述利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
气相色谱-离子迁移谱仪器条件设定为：色谱柱为多毛细管分离柱，固定相OV-5、色谱柱温40℃、漂移管温度45℃、进样口温度80℃、进样器-色谱柱管路温度65℃、载气为高纯氮气。
8.如权利要求7所述的方法，其特征在于，所述利用顶空-气相色谱-离子迁移谱技术获取各待检测料液样品的气相离子迁移谱图包括：
载气流速采用程序升速模式，其中，0～1min载气速度为2ml/min，1～3min载气速度由2ml/min逐渐上升到50ml/min，3～5min载气速度由50ml/min逐渐上升到150ml/min，5～20min载气速度维持在150ml/min；迁移谱漂移气体流量150ml/min。
CN202110446263.5A 2021-04-25 2021-04-25 特征图谱提取方法及料液配制质量稳定性的检测方法 Active CN113138248B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110446263.5A CN113138248B (zh) 2021-04-25 2021-04-25 特征图谱提取方法及料液配制质量稳定性的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110446263.5A CN113138248B (zh) 2021-04-25 2021-04-25 特征图谱提取方法及料液配制质量稳定性的检测方法

Publications (2)

Publication Number Publication Date
CN113138248A CN113138248A (zh) 2021-07-20
CN113138248B true CN113138248B (zh) 2022-08-02

Family

ID=76813530

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110446263.5A Active CN113138248B (zh) 2021-04-25 2021-04-25 特征图谱提取方法及料液配制质量稳定性的检测方法

Country Status (1)

Country Link
CN (1) CN113138248B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640445B (zh) * 2021-08-11 2024-06-11 贵州中烟工业有限责任公司 基于图像处理的特征峰识别方法及计算设备、存储介质
CN113533581B (zh) * 2021-08-11 2022-08-02 贵州中烟工业有限责任公司 基于信息熵分析物质质量的方法及计算机可读存储介质
CN113610817B (zh) * 2021-08-11 2024-03-26 贵州中烟工业有限责任公司 一种特征峰识别方法及计算设备、存储介质
CN113960152B (zh) * 2021-10-18 2023-06-06 北方工业大学 一种基于离子迁移谱的痕量爆炸物识别方法及系统
CN114113391B (zh) * 2021-11-24 2024-03-26 贵州中烟工业有限责任公司 一种膨胀介质含量检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0299652A1 (en) * 1987-07-15 1989-01-18 Hewlett-Packard Company Baseline correction for chromatography
CN105181852A (zh) * 2015-10-12 2015-12-23 西安石油大学 一种鉴别延胡索生熟药材的方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101534933B (zh) * 2006-05-26 2013-03-27 沃特世科技公司 关于n维数据的离子检测和参数估计
US20130080073A1 (en) * 2010-06-11 2013-03-28 Waters Technologies Corporation Techniques for mass spectrometry peak list computation using parallel processing
CN103198331B (zh) * 2013-03-25 2017-02-08 江苏易谱恒科技有限公司 基于pca_lda分析的多谱图特征融合识别方法
CN107402192A (zh) * 2017-03-03 2017-11-28 广西中烟工业有限责任公司 一种快速分析香精香料质量稳定性的方法
CN109655566A (zh) * 2017-10-12 2019-04-19 贵州中烟工业有限责任公司 一种鉴别卷烟挥发性成分稳定性的方法
CN109655530A (zh) * 2017-10-12 2019-04-19 贵州中烟工业有限责任公司 一种鉴别香精香料质量差异的方法
CN109655529A (zh) * 2017-10-12 2019-04-19 贵州中烟工业有限责任公司 一种鉴别卷烟真伪的方法
CN109655532A (zh) * 2017-10-12 2019-04-19 贵州中烟工业有限责任公司 一种对卷烟分类鉴别的方法
CN109668992A (zh) * 2017-10-13 2019-04-23 贵州中烟工业有限责任公司 一种霉变烟叶识别方法
US11398373B2 (en) * 2018-07-25 2022-07-26 Micromass Uk Limited User interface for ion mobility separation device
CN110097127A (zh) * 2019-05-07 2019-08-06 杭州麦迪特检测技术服务有限公司 基于核主成分分析的黄芪谱图特征提取与产地鉴别方法
CN111141860B (zh) * 2019-12-19 2023-01-10 中国农业科学院茶叶研究所 基于离子迁移谱的绿茶香气动态指纹分析及香型判别方法
CN111044638A (zh) * 2019-12-31 2020-04-21 贵州中烟工业有限责任公司 一种对不同品种烤烟烟叶分类鉴别的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0299652A1 (en) * 1987-07-15 1989-01-18 Hewlett-Packard Company Baseline correction for chromatography
CN105181852A (zh) * 2015-10-12 2015-12-23 西安石油大学 一种鉴别延胡索生熟药材的方法

Also Published As

Publication number Publication date
CN113138248A (zh) 2021-07-20

Similar Documents

Publication Publication Date Title
CN113138248B (zh) 特征图谱提取方法及料液配制质量稳定性的检测方法
CN109655532A (zh) 一种对卷烟分类鉴别的方法
CN110441423A (zh) 一种测定粮食香气成分的方法及其系统
CN111060642A (zh) 一种对同品种不同产地烟叶分类鉴别的方法
CN103063815A (zh) 一种辛味中药材产地的快速鉴别方法
Zhu et al. Volatile-based prediction of sauvignon blanc quality gradings with static headspace–gas chromatography–ion mobility spectrometry (SHS–GC–IMS) and interpretable machine learning techniques
CN111044638A (zh) 一种对不同品种烤烟烟叶分类鉴别的方法
CN103134850A (zh) 一种基于特征香气的茶叶品质快速检测装置及检测方法
CN109738569B (zh) 一种基于spme-gc-ms的高仿卷烟鉴别方法
CN112578039B (zh) 一种燕窝中掺假物的分析鉴定方法
CN104316635A (zh) 水果风味品质的快速鉴别方法
CN101487826B (zh) 一种识别白酒香型和/或等级的方法
Vest Nielsen et al. Full second-order chromatographic/spectrometric data matrices for automated sample identification and component analysis by non-data-reducing image analysis
Pérez-Jiménez et al. Application of untargeted volatile profiling and data driven approaches in wine flavoromics research
CN109655566A (zh) 一种鉴别卷烟挥发性成分稳定性的方法
CN104034820A (zh) 一种快速鉴别樱桃酒品牌的方法
CN109655533A (zh) 一种鉴别香精香料类型的方法
CN111505174A (zh) 一种鉴别真假烟的方法
CN113640445B (zh) 基于图像处理的特征峰识别方法及计算设备、存储介质
CN113295813B (zh) 一种料液施加均匀性的检测方法
CN113075316B (zh) 一种靖西大果山楂酒窖藏时间的鉴定方法
CN108445134B (zh) 酒类产品鉴别的方法
CN103293257A (zh) 一种卷烟烟气色谱数据保留时间漂移mwfft校准方法
CN110887921B (zh) 一种高效快速分析杜仲叶及其发酵产物特征挥发性成分的方法
CN106908527A (zh) 一种鉴别荔枝蜜产地的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant