CN114184576A - 一种基于分子光谱流程工业在线测量方法及系统 - Google Patents
一种基于分子光谱流程工业在线测量方法及系统 Download PDFInfo
- Publication number
- CN114184576A CN114184576A CN202111224106.6A CN202111224106A CN114184576A CN 114184576 A CN114184576 A CN 114184576A CN 202111224106 A CN202111224106 A CN 202111224106A CN 114184576 A CN114184576 A CN 114184576A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- current stage
- component
- spectral data
- target sample
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 99
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims abstract description 57
- 239000003550 marker Substances 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 29
- 238000012216 screening Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000000491 multivariate analysis Methods 0.000 claims abstract description 14
- 238000010905 molecular spectroscopy Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 41
- 238000007781 pre-processing Methods 0.000 claims description 15
- 230000008030 elimination Effects 0.000 claims description 11
- 238000003379 elimination reaction Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 10
- 238000000513 principal component analysis Methods 0.000 claims description 7
- 230000002441 reversible effect Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 230000002860 competitive effect Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000000556 factor analysis Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000000750 progressive effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 18
- 238000002329 infrared spectrum Methods 0.000 abstract description 13
- 238000009776 industrial production Methods 0.000 abstract description 6
- 238000005259 measurement Methods 0.000 abstract description 5
- 230000000875 corresponding effect Effects 0.000 description 26
- 239000000463 material Substances 0.000 description 20
- 238000004519 manufacturing process Methods 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 14
- 238000012937 correction Methods 0.000 description 12
- 239000000203 mixture Substances 0.000 description 5
- 238000010986 on-line near-infrared spectroscopy Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000036961 partial effect Effects 0.000 description 3
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000001237 Raman spectrum Methods 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000002211 ultraviolet spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/33—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N2021/3125—Measuring the absorption by excited molecules
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于分子光谱流程工业在线测量方法及系统，涉及流程工业在线测量技术领域，该方法包括获取目标样本当前阶段对应的光谱数据；基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。本发明无需建模，极大的缩短了在线近红外光谱技术的调试建模时间，对于流程工业生产具有重要意义。
Description
技术领域
本发明涉及流程工业在线测量技术领域，特别是涉及一种基于分子光谱流程工业在线测量方法及系统。
背景技术
流程工业智能制造是过程信息化与自动化的融合。其中，信息化占主导地位，而过程信息化的核心则是过程信息感知技术。与成熟的诸多过程信息(如温度、流量等)感知技术相比，过程物料组成及关键性质的过程分析技术则是流程工业智能制造亟待发展的信息感知技术。近红外光谱可从分子水平上反映物质组成与结构信息，是理想的过程分析信号。随着光纤、光学材料和仪器制造技术的快速发展，可以方便、快速地获取流程工业过程物料的近红外光谱信号，与化学计量学方法结合，可以实现工艺过程中物料组成和多种物化性质的实时检测(即过程物料信息感知技术)，对于提高流程工业过程优化与控制具有重要作用。基于近红外光谱的过程分析技术(ProcessAnalytical Technology，PAT)，已广泛用于石化、制药和半导体制造等流程工业领域，对其工业生产质量提升具有明显作用，有着光明的应用前景。
目前近红外光谱的过程分析技术主要是用反向区间PLS建立定量校正模型，然后将未知样本的近红外光谱数据带入定量校正模型中以预测在线生产物料的组分含量或者物理性质。建立一个稳健的、准确的且能够用于在线监测的定量校正模型大概需要收集1000个样本，甚至更多，且这1000个样本的性质要有一定的分散性，能够涵盖生产中的所有物料性质的变化。一般建立一个定量校正模型需要一年以上的时间。因此建立定量校正模型的过程将耗费大量的人力物力，耗时也较长。实际上，在工业生产中，在线测量的目的是为先进控制提供有效参数，实现生产的平稳运行，而在控制实施过程中往往只关心组分性质参数的变化趋势。因此一种无须建模即能够满足先进控制的近红外监测技术对于近红外光谱技术的应用以及成本的节约具有重要意义。
发明内容
本发明的目的是提供一种基于分子光谱流程工业在线测量方法及系统，极大的缩短了在线近红外光谱技术的调试建模时间，对于流程工业生产具有重要意义。
为实现上述目的，本发明提供了如下方案：
一种基于分子光谱流程工业在线测量方法，包括：
获取目标样本当前阶段对应的光谱数据；
基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
可选的，所述标记组分光谱的确定过程为：
获取历史光谱数据；所述历史光谱数据为在一段时间内采集的标定样本对应的光谱数据；所述目标样本和所述标定样本为在同一管线上的采集的物质；所述标定样本的个数为5～10个；
计算每个所述标定样本中标记组分对应的浓度值；
将所述历史光谱数据与每个所述标定样本中标记组分对应的浓度值进行对应组合，进而形成光谱矩阵；
采用波长变量筛选算法，对预处理后的光谱矩阵中的特征波长进行筛选，并记录所述特征波长对应的索引位置；
将预处理后的历史光谱数据与所述索引位置对应的数值进行组合，形成新矩阵；
采用多元分析手段对所述新矩阵进行处理，得到主成成分光谱矩阵；
采用渐进因子分析算法或者简化自模式算法对所述主成成分光谱矩阵进行处理，得到所述主成成分光谱矩阵的初始值；
基于交替最小二乘算法和所述主成成分光谱矩阵的初始值进行迭代优化，得到标记组分光谱。
可选的，所述基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图，具体包括：
对所述当前阶段对应的光谱数据进行预处理；
根据预处理后的当前阶段对应的光谱数据和所述索引位置对应的数值，得到目标波长；
基于所述目标波长和标记组分光谱，采用广义逆矩阵计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
可选的，所述基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图，具体包括：
对所述当前阶段对应的光谱数据进行预处理；
根据预处理后的当前阶段对应的光谱数据和标记组分光谱，采用正交投影和斜投影计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
可选的，所述预处理的方法包括导数、平滑、SNV、MSC、消除常数偏移、减去一条直线、离散小波变换、连续小波变换的一种或者几种的组合。
可选的，所述多元分析手段为主成分分析算法或者多元曲线分辨算法。
可选的，所述波长变量筛选算法为反向区间PLS法、无信息变量消除法、蒙特卡洛-无信息变量消除法或竞争自适用重加权采样法。
一种基于分子光谱流程工业在线测量系统，包括：
光谱数据获取模块，用于获取目标样本当前阶段对应的光谱数据；
浓度变化曲线绘制模块，用于基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种基于分子光谱流程工业在线测量方法及系统，根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到标记组分光谱，然后基于目标样本当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图，无需建模，极大的缩短了在线近红外光谱技术的调试建模时间，对于流程工业生产具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种基于分子光谱流程工业在线测量方法的流程示意图；
图2为本发明标记组分光谱确定过程示意图；
图3为本发明标记组分浓度预测过程示意图；
图4为本发明生产过程近红外光谱数据示意图；
图5为本发明CARS选择变量频率分布和所选择特征波长位置分布图；
图6为本发明近红外趋势预测结果示意图；
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在线近红外光谱技术在使用过程中需要建立校正模型，该校正模型建立周期长，耗费大量人力物力，且校正模型需要周期性维护，成本高。本发明的目的是为了解决在线近红外光谱技术使用过程需要建立校正模型的问题，使用非建模的方式预测物料的变化趋势，为先进过程控制提供可靠的性质参数，从而使生产平稳运行。本发明的使用极大的缩短了在线近红外光谱技术的调试建模时间，对于流程工业生产具有重要意义。
本发明的主要步骤如下：
1)使用合适的光谱处理手段(包括导数、消除常数偏移、去趋势校正、多元散射校正、标准正态变量变换等)对光谱数据进行预处理，消除部分无关信息的影响。
2)使用合适的变量筛选方法(反向区间PLS、无信息变量消除法、蒙特卡洛-无信息变量消除法、竞争自适用重加权采样(CARS)等)对变量进行筛选，选择最相关的变量进行后续运算。
3)使用主成分分析(PCA)或者多元曲线分辨(MCR)的方法从历史光谱数据中提取“纯组分信号”。
4)未知待分析信号向“纯组分信号”做正交投影计算、斜投影或者广义逆计算，计算物料的变化趋势。
下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
目前在线近红外光谱技术需要将近红外光谱与实验室分析数据通过化学计量学软件进行关联，建立校正模型。在应用时将校正模型导入在线分析软件，在线分析软件采集到未知物料性质的光谱时即可预测该物料性质。建立稳定可靠的模型一般需要至少1000个具有代表性样本，由于在线生产工艺相对稳定，样本性质变化较小，因此建模周期长，严重影响在线近红外光谱技术在流程工业中的应用。
请参见图1，本实施例提供的一种基于分子光谱流程工业在线测量方法，包括：
步骤101：获取目标样本当前阶段对应的光谱数据。
步骤102：基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图。
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
本实施例使用变量筛选算法结合PCA或者MCR等多元分析手段分析历史光谱数据。本实施例变量筛选过程与传统的建模相比，本实施例使用5～10个代表性样本进行变量筛选，选择最具代表性的特征波长用于趋势分析。
1)在生产线上安装在线近红外光谱分析系统，每隔一定时间(1min、3min、10min等)收集一条光谱数据，将一段时间内收集的光谱数据组成光谱矩阵D1，即历史光谱数据
2)稳定生产时在管线上取一个样本进行实验室分析，得到该测点物料最关心的一个或者几个性质值，本实施例所述性质值为浓度值；在工艺波动时在管线上取多个样本进行实验室分析，得到相同的性质值，并记录每个样本的详的取样时间，总共记录5～10个样本的性质值。
3)根据取样时间将实验室分析值，即浓度值与历史光谱数据相对应，这一系列操作后组成光谱矩阵D2，组分浓度为c2。
4)对光谱矩阵D2进行预处理，其目的是消除光谱漂移、噪音、基线等干扰因素的影响，预处理方法为导数、平滑、SNV、MSC、消除常数偏移、减去一条直线、离散小波变换、连续小波变换的一种或者几种的组合。
5)使用波长变量筛选算法(反向区间PLS法、无信息变量消除法、蒙特卡洛-无信息变量消除法、竞争自适用重加权采样(CARS)法)对预处理后的光谱矩阵D2进行处理，选择最合适的波长变量，即特征波长，并记录该特征波长的索引位置，其目的是找到与待测组分浓度最相关的波长变量。
6)对光谱矩阵D1进行与步骤4)相同的预处理；选择步骤5)特征波长的索引位置的数值组成新矩阵D3。
7)采用多元分析手段对新矩阵D3进行处理，得到主成成分光谱矩阵。
8)对主成成分光谱矩阵使用渐进因子分析(EFA)或者使用简化自模式(SIMPLISMA)算法估算初始值。
9)假设主成成分光谱矩阵为各“纯组分”光谱与其浓度乘积的线性加和后组成的光谱矩阵。
在封闭性、非负性等一种或者多种的约束条件下使用交替最小二乘(als)算法进行如公式(2)和公式(3)的迭代计算。
D＝CST (1)；
ST＝C+D (2)；
C＝D(ST)+ (3)；
当迭代计算前后两次光谱的残差小于给定阈值或者迭代次数达到100次停止迭代。
其中，D为原始光谱矩阵，C为浓度矩阵，S为组成体系的纯组分光谱，即标记组分光谱，C+和(ST)+分别是C和ST的广义逆矩阵。
10)当光谱仪收集到新的未知浓度信息的光谱数据时，即获取到目标样本当前阶段对应的光谱数据，根据步骤4)的预处理方法进行预处理，并根据步骤5)记录的索引位置数值得到新的向量d，将公式(3)中的D改为d，即可计算该光谱的组分浓度信息。
11)随着时间的变化，不断收集未知物料信息的光谱，求得组分浓度信息，一系列组分浓度信息即可组成浓度变化趋势，即根据预处理后的当前阶段对应的光谱数据和所述索引位置对应的数值，得到目标波长；基于所述目标波长和标记组分光谱，采用广义逆矩阵计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
或者根据公式(4)和公式(5)进行斜投影和正交投影计算，得到组分浓度信息。
ES|H＝S(STPH ⊥S)-1STPH ⊥ (4)。
PH ⊥＝I-PH＝I-H(HTH)-1HT (5)。
式中，S为被测组分的纯光谱信号，H为背景光谱信号，I为与PH维数相同的单位矩阵，T为矩阵的转置，PH ⊥为正交投影算子，ES|H为斜投影算子。
11)随着时间的变化，不断收集未知物料信息的光谱，求得组分浓度信息，一系列组分浓度信息即可组成浓度变化趋势，即根据预处理后的当前阶段对应的光谱数据和标记组分光谱，采用正交投影和斜投影计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
其中，步骤(1)～步骤(9)为历史数据建立模型步骤，其流程图如图2所示。步骤(10)～步骤(11)为预测未来数据步骤，具体计算流程如图3所示。
实施例二
下面通过一条具体生产线来说明本发明技术方案
1)在某流程工业生产线上安装傅里叶型在线近红外光谱(波长范围10000～4000cm-1，分辨率8cm-1，扫描次数32次)分析系统，每1分钟收集一条光谱数据，将一段时间内收集的光谱数据组成光谱矩阵D1。
2)稳定生产时，在管线上取一个样进行实验室分析，得到该测点物料最关心的一个性质值(本测点为M2含量)；在工艺波动时，取样进行M2含量的实验室分析，取样分析时记录详细的取样时间，总共记录10个样本的性质值。
3)根据取样时间将实验室分析值(即性质值)与光谱相对应，这一系列光谱组成光谱矩阵D2，组分性质为c2。
4)对光谱矩阵D2进行一阶导数(Savitzky-Golay法，拟合次数2次，窗口宽度为17)和SNV预处理。
一阶导数计算公式为：
其中，窗口宽度为2m+1，A为归一化常数，x′i为光谱数据xi的导数，wj为相应的导数系数，在窗口宽度确定后，wj就确定了。
SNV的计算公式为：
其中，
m为变量个数，i＝1,2,3,…n，n为样本数目。
5)使用CARS选择最合适的波长变量，并记录该变量的索引位置。
6)对D1进行与步骤4)相同的预处理并对所有数据取绝对值处理；选择步骤5)变量索引位置的数值组成新的矩阵D3。
7)对D3做主成分分析计算主成分贡献率，根据主成分贡献率确定保留的组分数。
8)对D3使用简化自模式(SIMPLISMA)算法估算初始值。
9)按照式(2)和(3)进行多元曲线分辨计算，对S设定非负性约束，对C设定封闭性和非负性约束。当前后两次迭代S的平均标准偏差小于0.001时或者迭代次数达到100次时停止迭代。迭代计算后的S即用于下一步计算。
10)当光谱仪收集到新的未知浓度信息的光谱，根据步骤4)预处理方法进行预处理并求绝对值，根据5)记录的索引位置数值得到新的向量d，根据公式(4)进行斜投影计算得到各个分辨光谱的组分浓度信息。
11)再在线取样分析用于验证分辨光谱的准确性，取10个样本按照步骤10进行计算，计算各个分辨光谱计组分浓度信息和取样分析组分浓度的相关系数，相关系数最大的浓度趋势用作该物料的主要浓度变化趋势。
12)随着时间的变化，不断收集未知物料信息的光谱，求得组分浓度信息，一系列组分浓度信息即可组成浓度变化趋势。
对光谱进行导数和SNV预处理可以消除光谱由加性效应和乘性效应引起的的噪音，Savitzky-Golay法不仅可以对光谱求导数，还可以对光谱进行平滑处理，从而消除部分仪器和环境变化引起的高频噪音，SNV可以消除光程微量变化、杂散光等影响。因此Savitzky-Golay导数结合SNV方法对在线连续收集的近红外光谱预处理，可以最大程度的消除噪音信息。
CARS具有“优胜劣汰，从中取优”的特点,首先以指数函数快速递减的方式快速剔除大部分无用变量，再将剩余变量随机组合选出最优变量组合，所选择的变量能够最大限度的代表组分信息，从而排除其他无关变量的干扰。因此该发明选择使用CARS选择少数几个变量用于迭代计算。
多元曲线分辨(MCR)的目的是将各个纯组分的光谱及其在体系中的含量从复杂的混合物光谱中分辨出来，与PCA不同，MCR-ALS作为MCR的重要手段在实施过程中通常施加一定的约束条件，如单峰性、封闭性、非负性，通过迭代计算分辨所得的组分光谱具有真实的物理和化学意义。因此本发明选择MCR-ALS用于求取所关心组分的“真实光谱”。
斜投影算法是通过数据向量空间描述化学复杂体系，其数据向量空间M被分为两部分：被测变量的向量子空间S和被测组分之外变量，即向量相邻子空间H。子空间S指体系中被测组分的纯信号，即投影方向。相邻子空间H指体系中除被测组分外的其余组分信号，即投影空间，有M＝S⊕H。斜投影算子ES|H是将待测组分纯信号从混合信号中分离出来。
如图4所示，收集了一系列在线生产的近红外光谱数据。把10条具有实验室化验数据的光谱进行预处理，并选择特征波长，其预处理光谱和位置如图5所示，图5(a)中是CARS计算的选中变量频率分布直方图，图5(b)是预处理光谱和选中变量的位置。使用选中变量进行MCR-ALS计算，对未来一天的数据进行预测，其预测趋势图如图6所示，可以发现两者预测趋势图几乎完全一致，趋势预测结果噪音更小，更有利于控制过程。
实施例三
本实施例提供了一种基于分子光谱流程工业在线测量系统，包括：
光谱数据获取模块，用于获取目标样本当前阶段对应的光谱数据；
浓度变化曲线绘制模块，用于基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
此外，预处理方法中SNV可以换成MSC、基线校正等；一阶导数可以换成二阶导数；导数方法可以换成Norris求导法、差分求导法或者连续小波变换法。
变量筛选方法中CARS可以换成反向区间偏最小二乘、无信息变量消除(UVE)、蒙特卡洛-无信息变量消除(MC-UVE)、遗传算法、相关系数法等。
本发明是在线光谱分析的数据处理方法，因此也适用于拉曼光谱、紫外光谱和中红外光谱等分子光谱。
本发明交替最小二乘(ALS)法计算时，给定的初始值为初始光谱矩阵，也可以给定初始浓度矩阵
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于分子光谱流程工业在线测量方法，其特征在于，包括：
获取目标样本当前阶段对应的光谱数据；
基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
2.根据权利要求1所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述标记组分光谱的确定过程为：
获取历史光谱数据；所述历史光谱数据为在一段时间内采集的标定样本对应的光谱数据；所述目标样本和所述标定样本为在同一管线上的采集的物质；所述标定样本的个数为5～10个；
计算每个所述标定样本中标记组分对应的浓度值；
将所述历史光谱数据与每个所述标定样本中标记组分对应的浓度值进行对应组合，进而形成光谱矩阵；
采用波长变量筛选算法，对预处理后的光谱矩阵中的特征波长进行筛选，并记录所述特征波长对应的索引位置；
将预处理后的历史光谱数据与所述索引位置对应的数值进行组合，形成新矩阵；
采用多元分析手段对所述新矩阵进行处理，得到主成成分光谱矩阵；
采用渐进因子分析算法或者简化自模式算法对所述主成成分光谱矩阵进行处理，得到所述主成成分光谱矩阵的初始值；
基于交替最小二乘算法和所述主成成分光谱矩阵的初始值进行迭代优化，得到标记组分光谱。
3.根据权利要求2所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图，具体包括：
对所述当前阶段对应的光谱数据进行预处理；
根据预处理后的当前阶段对应的光谱数据和所述索引位置对应的数值，得到目标波长；
基于所述目标波长和标记组分光谱，采用广义逆矩阵计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
4.根据权利要求1所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图，具体包括：
对所述当前阶段对应的光谱数据进行预处理；
根据预处理后的当前阶段对应的光谱数据和标记组分光谱，采用正交投影和斜投影计算当前阶段目标样本中标记组分的各个浓度值，并根据所有所述浓度值绘制当前阶段目标样本中标记组分的浓度变化曲线图。
5.根据权利要求2、3或者4所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述预处理的方法包括导数、平滑、SNV、MSC、消除常数偏移、减去一条直线、离散小波变换、连续小波变换的一种或者几种的组合。
6.根据权利要求1或者2所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述多元分析手段为主成分分析算法或者多元曲线分辨算法。
7.根据权利要求2所述的一种基于分子光谱流程工业在线测量方法，其特征在于，所述波长变量筛选算法为反向区间PLS法、无信息变量消除法、蒙特卡洛-无信息变量消除法或竞争自适用重加权采样法。
8.一种基于分子光谱流程工业在线测量系统，其特征在于，包括：
光谱数据获取模块，用于获取目标样本当前阶段对应的光谱数据；
浓度变化曲线绘制模块，用于基于所述当前阶段对应的光谱数据和标记组分光谱，绘制当前阶段目标样本中标记组分的浓度变化曲线图；
其中，所述标记组分光谱是根据波长变量筛选算法和多元分析手段对历史光谱数据进行处理得到的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111224106.6A CN114184576B (zh) | 2021-10-19 | 2021-10-19 | 一种基于分子光谱流程工业在线测量方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111224106.6A CN114184576B (zh) | 2021-10-19 | 2021-10-19 | 一种基于分子光谱流程工业在线测量方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114184576A true CN114184576A (zh) | 2022-03-15 |
| CN114184576B CN114184576B (zh) | 2024-08-13 |
Family
ID=80539810
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111224106.6A Active CN114184576B (zh) | 2021-10-19 | 2021-10-19 | 一种基于分子光谱流程工业在线测量方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114184576B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5121337A (en) * | 1990-10-15 | 1992-06-09 | Exxon Research And Engineering Company | Method for correcting spectral data for data due to the spectral measurement process itself and estimating unknown property and/or composition data of a sample using such method |
| US5532487A (en) * | 1994-11-23 | 1996-07-02 | E. I. Du Pont De Nemours And Company | Near-infrared measurement and control of polyamide processes |
| US20160000378A1 (en) * | 2014-05-02 | 2016-01-07 | David R. Hall | Human Health Property Monitoring System |
| WO2017045296A1 (zh) * | 2015-09-14 | 2017-03-23 | 上海创和亿电子科技发展有限公司 | 在线近红外样本量确定方法 |
| WO2017084119A1 (zh) * | 2015-11-19 | 2017-05-26 | 江南大学 | 一种具有无测点温度补偿功能的近红外物性参数测量方法 |
| CN107064054A (zh) * | 2017-02-28 | 2017-08-18 | 浙江大学 | 一种基于cc‑pls‑rbfnn优化模型的近红外光谱分析方法 |
| US20190049297A1 (en) * | 2015-11-19 | 2019-02-14 | Jiangnan University | Method for Correcting Measuring-Point-Free Temperature Compensation Model During Online Application of Near Infrared Spectrum Analyzer |
| CN110320165A (zh) * | 2019-08-08 | 2019-10-11 | 华南农业大学 | 香蕉可溶性固形物含量的可见/近红外光谱无损检测方法 |
| US20200049623A1 (en) * | 2018-08-13 | 2020-02-13 | Oil Crops Research Institute, Chinese Academy Of Agricultural Sciences | Nir spectroscopy method for fatty acid content of oilseeds |
| CN111999258A (zh) * | 2020-07-03 | 2020-11-27 | 桂林理工大学 | 一种面向光谱基线校正的加权建模局部优化方法 |
| CN113049507A (zh) * | 2021-03-09 | 2021-06-29 | 三峡大学 | 多模型融合的光谱波长选择方法 |
-
2021
- 2021-10-19 CN CN202111224106.6A patent/CN114184576B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5121337A (en) * | 1990-10-15 | 1992-06-09 | Exxon Research And Engineering Company | Method for correcting spectral data for data due to the spectral measurement process itself and estimating unknown property and/or composition data of a sample using such method |
| US5532487A (en) * | 1994-11-23 | 1996-07-02 | E. I. Du Pont De Nemours And Company | Near-infrared measurement and control of polyamide processes |
| US20160000378A1 (en) * | 2014-05-02 | 2016-01-07 | David R. Hall | Human Health Property Monitoring System |
| WO2017045296A1 (zh) * | 2015-09-14 | 2017-03-23 | 上海创和亿电子科技发展有限公司 | 在线近红外样本量确定方法 |
| WO2017084119A1 (zh) * | 2015-11-19 | 2017-05-26 | 江南大学 | 一种具有无测点温度补偿功能的近红外物性参数测量方法 |
| US20190049297A1 (en) * | 2015-11-19 | 2019-02-14 | Jiangnan University | Method for Correcting Measuring-Point-Free Temperature Compensation Model During Online Application of Near Infrared Spectrum Analyzer |
| CN107064054A (zh) * | 2017-02-28 | 2017-08-18 | 浙江大学 | 一种基于cc‑pls‑rbfnn优化模型的近红外光谱分析方法 |
| US20200049623A1 (en) * | 2018-08-13 | 2020-02-13 | Oil Crops Research Institute, Chinese Academy Of Agricultural Sciences | Nir spectroscopy method for fatty acid content of oilseeds |
| CN110320165A (zh) * | 2019-08-08 | 2019-10-11 | 华南农业大学 | 香蕉可溶性固形物含量的可见/近红外光谱无损检测方法 |
| CN111999258A (zh) * | 2020-07-03 | 2020-11-27 | 桂林理工大学 | 一种面向光谱基线校正的加权建模局部优化方法 |
| CN113049507A (zh) * | 2021-03-09 | 2021-06-29 | 三峡大学 | 多模型融合的光谱波长选择方法 |
Non-Patent Citations (3)
| Title |
|---|
| 屈东胜 等: "基于波长调制光谱的多参数测量方法研究", 红外与毫米波学报, vol. 35, no. 04, 31 December 2016 (2016-12-31), pages 470 - 476 * |
| 李斌 等: "基于可见-近红外光谱的茄子叶绿素荧光参数F_v/F_m预测方法", 光谱学与光谱分析, no. 09, 7 September 2020 (2020-09-07), pages 2834 - 2839 * |
| 胡军 等: "基于CARS波段筛选的甲醇汽油中甲醇含量中红外光谱检测", 激光与光电子学进展, no. 23, 31 December 2019 (2019-12-31), pages 281 - 287 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114184576B (zh) | 2024-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105300923B (zh) | 一种近红外光谱分析仪在线应用时无测点温度补偿模型修正方法 | |
| CN107817223A (zh) | 快速无损实时预测原油性质模型的构建方法及其应用 | |
| CN110161013B (zh) | 基于机器学习的激光诱导击穿光谱数据处理方法和系统 | |
| CN107703097B (zh) | 利用近红外光谱仪构建快速预测原油性质的模型的方法 | |
| WO2017084119A1 (zh) | 一种具有无测点温度补偿功能的近红外物性参数测量方法 | |
| Bro et al. | Theory of net analyte signal vectors in inverse regression | |
| CN101915744A (zh) | 物质成分含量的近红外光谱无损检测方法及装置 | |
| CN105319198A (zh) | 基于拉曼光谱解析技术的汽油苯含量预测方法 | |
| CN104596975A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶木质素的方法 | |
| CN103134763B (zh) | 由红外光谱预测原油密度的方法 | |
| CN104596979A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶纤维素的方法 | |
| CN109283153B (zh) | 一种酱油定量分析模型的建立方法 | |
| Li et al. | A feasibility study on quantitative analysis of low concentration methanol by FT-NIR spectroscopy and aquaphotomics | |
| CN114112978A (zh) | 原位红外光谱检测结晶过程溶液浓度的小波函数标定建模方法 | |
| Andrade‐Garda et al. | Partial Least‐Squares Regression | |
| CN105466885A (zh) | 基于无测点温度补偿机制的近红外在线测量方法 | |
| CN107356556A (zh) | 一种近红外光谱定量分析的双集成建模方法 | |
| CN116662751A (zh) | 一种基于主成分分析与杠杆值法剔除异常样本的烟叶含水率检测方法 | |
| CN114184576B (zh) | 一种基于分子光谱流程工业在线测量方法及系统 | |
| CN104596982A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶果胶的方法 | |
| CN114414524A (zh) | 快速检测航空煤油性质的方法 | |
| CN110398488B (zh) | 应用于激光诱导击穿光谱的非线性定量分析方法 | |
| CN103134762B (zh) | 由红外光谱预测原油氮含量的方法 | |
| CN112861412A (zh) | 基于近红外光谱主成分和神经网络的生物质挥发分含量测量与建模方法 | |
| CN105259135B (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 |