CN109142263B - 一种固体推进剂制备过程的组分分散均匀性在线测试方法 - Google Patents
一种固体推进剂制备过程的组分分散均匀性在线测试方法 Download PDFInfo
- Publication number
- CN109142263B CN109142263B CN201810768568.6A CN201810768568A CN109142263B CN 109142263 B CN109142263 B CN 109142263B CN 201810768568 A CN201810768568 A CN 201810768568A CN 109142263 B CN109142263 B CN 109142263B
- Authority
- CN
- China
- Prior art keywords
- component
- sample
- content
- near infrared
- correction
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 51
- 239000004449 solid propellant Substances 0.000 title claims abstract description 45
- 238000002360 preparation method Methods 0.000 title claims abstract description 24
- 238000010998 test method Methods 0.000 title description 17
- 238000000034 method Methods 0.000 claims abstract description 114
- 238000012937 correction Methods 0.000 claims abstract description 104
- 238000012360 testing method Methods 0.000 claims abstract description 84
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 55
- 238000001228 spectrum Methods 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 238000012795 verification Methods 0.000 claims description 56
- 239000007787 solid Substances 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 17
- 238000002835 absorbance Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 12
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 8
- 238000007906 compression Methods 0.000 claims description 8
- 229910052736 halogen Inorganic materials 0.000 claims description 7
- 150000002367 halogens Chemical class 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 7
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 claims description 7
- 229910052721 tungsten Inorganic materials 0.000 claims description 7
- 239000010937 tungsten Substances 0.000 claims description 7
- 239000000020 Nitrocellulose Substances 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 229920001220 nitrocellulos Polymers 0.000 claims description 6
- 239000004575 stone Substances 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 229920000642 polymer Polymers 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 4
- 239000005062 Polybutadiene Substances 0.000 claims description 4
- GYZLOYUZLJXAJU-UHFFFAOYSA-N diglycidyl ether Chemical compound C1OC1COCC1CO1 GYZLOYUZLJXAJU-UHFFFAOYSA-N 0.000 claims description 4
- 229920002857 polybutadiene Polymers 0.000 claims description 4
- 238000003908 quality control method Methods 0.000 claims description 4
- 238000011426 transformation method Methods 0.000 claims description 4
- -1 3-azidomethyl-3-methylepoxybutane Chemical compound 0.000 claims description 3
- 238000013528 artificial neural network Methods 0.000 claims description 2
- 230000008033 biological extinction Effects 0.000 claims description 2
- 238000013098 chemical test method Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 230000002068 genetic effect Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 238000012843 least square support vector machine Methods 0.000 claims description 2
- 238000012628 principal component regression Methods 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims description 2
- 229910052724 xenon Inorganic materials 0.000 claims description 2
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 claims description 2
- 239000000463 material Substances 0.000 description 11
- 238000000862 absorption spectrum Methods 0.000 description 5
- 238000009614 chemical analysis method Methods 0.000 description 5
- 239000007800 oxidant agent Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000000945 filler Substances 0.000 description 4
- 230000001590 oxidative effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000004014 plasticizer Substances 0.000 description 3
- 229920002121 Hydroxyl-terminated polybutadiene Polymers 0.000 description 2
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 2
- 230000001070 adhesive effect Effects 0.000 description 2
- 230000002860 competitive effect Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000001856 Ethyl cellulose Substances 0.000 description 1
- ZZSNKZQZMQGXPY-UHFFFAOYSA-N Ethyl cellulose Chemical compound CCOCC1OC(OC)C(OCC)C(OCC)C1OC1C(O)C(O)C(OC)C(CO)O1 ZZSNKZQZMQGXPY-UHFFFAOYSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 150000001540 azides Chemical class 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 229920001249 ethyl cellulose Polymers 0.000 description 1
- 235000019325 ethyl cellulose Nutrition 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 239000000825 pharmaceutical preparation Substances 0.000 description 1
- 238000012545 processing 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/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (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)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种固体推进剂制备过程的组分分散均匀性在线测试方法，包括以下步骤：首先采集样品集的近红外光谱数据；之后测取样品集样品各组分的含量实测值；然后建立组分含量的近红外校正预测模型；之后取待测样品并获取待测样品各组分的含量预测值；最后判断待测样品组分分散均匀性，不断重复最后两个步骤实现对所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。本发明以近红外光谱模型结合实时测试的光谱数据实现固体推进剂组分含量及分散均匀性的在线测试。本发明方法操作简单、快速、准确，解决了传统方法无法实现在固体推进剂制备过程中对组分含量及组分分散均匀性等在线测试的难题。
Description
技术领域
本发明属于固体推进剂的检测领域，特别是一种固体推进剂制备过程的组分分散均匀性在线测试方法。
背景技术
固体推进剂是火箭发动机的动力来源，固体推进剂在成型加工过程中质量控制和检测非常重要，传统的检测方法是对加工成型后的成品进行各种理化性能测试后综合判定其质量，检测周期通常至少需要2-5天的时间，过程质量控制也主要通过工艺参数的控制来保证，在线或快速的过程质量测试方法非常少。随着火箭系统对固体推进剂综合质量的不断提高，需要开发快速或在线的质量检测方法和测试系统，保证产品生产的过程质量。
中国实用新型专利CN201621369363.3公开了一种玻璃纤维短切纱均匀性连续在线检测系统，该系统是基于图像识别技术对物料的均匀性进行检测的，其原理不同于近红外光谱测试技术。中国药事2008年22卷第7期560-563页的文献公开了微型近红外分析仪用于药物制剂混合中的在线检测，其是采用移动窗口标准偏差法直接分析物料的混合均匀性，实施判定混合终点，该方法主要用于制药混合环节的测试，对于工序较多且每个工序的物料状态都有很大差异测固体推进剂制备过程是不适合的。
近红外光谱技术是近年来快速发展的一项无损检测技术，可以用于含有氢元素的有机材料的物化性能分析，具有操作简单、样品无需前处理、易于实现快速分析的优势，已广泛应用到医药、石化、食品及其他工农业领域。但至今尚没有针对不规则的固体推进剂的组分及物料分散均匀性的在线检测系统和对应方法。
发明内容
本发明要解决的技术问题在于提供一种固体推进剂制备过程的组分分散均匀性在线测试方法。
实现本发明目的的技术解决方案为：一种固体推进剂制备过程的组分分散均匀性在线测试方法，包括以下步骤：
步骤1、准备不同配方的固体推进剂在制过程的固体样品，随机选取若干固体样品作为样品集，在一定的光谱测量条件下利用近红外光谱仪采集样品集的近红外光谱数据；
步骤2、采用常规检测方法测得样品集中每个样品各组分的含量实测值；
步骤3、建立组分含量的近红外校正预测模型；
步骤4、取待测样品，按照步骤1中的光谱测量条件采集待测样品的近红外光谱数据，并导入步骤3建立的近红外校正预测模型中，经模型运算，获得待测样品各组分的含量预测值；
步骤5、根据步骤4获得的待测样品各组分含量预测值和步骤2中每个样品各组分的含量实测值，判断待测样品组分分散均匀性；
步骤6、重复步骤4和步骤5，实现对所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
本发明与现有技术相比，其显著优点为：1)本发明的方法步骤简单，节省人力、物力成本；2)本发明方法适用于各种固体推进剂物料体系的组分和分散均匀性在线测试，具有检测速度快、分析效率高、稳定性及重复性好的优势；3)本发明可以在1～5秒完成固体推进剂的主要成分分析，降低了检测成本的同时大大提高了测试效率；4)本发明的方法在检测过程不需要使用任何化学试剂，不会对环境造成任何污染，解决了传统方法无法实现在固体推进剂制备过程中对组分含量及组分分散均匀性等在线测试的难题。
下面将结合附图说明对本发明作进一步详细描述。
附图说明
图1为本发明固体推进剂制备过程的组分分散均匀性在线测试方法的流程图。
图2为本发明实施例1中主要成份为NC、含能增塑剂和高能填料的固体推进剂物料的近红外光谱数据图。
图3为本发明实施例1固体推进剂的近红外预测模型效果图。其中，图(a)为NC组分的实测值与预测值的线性关系图，图(b)为含能增塑剂组分的实测值与预测值的线性关系图，图(c)为高能填料组分的实测值与预测值的线性关系图。
图4为本发明实施例2中主要成份为端羟基聚丁二烯、高能氧化剂和RDX的固体推进剂物料的近红外光谱数据图。
图5为本发明实施例2固体推进剂的近红外预测模型效果图。其中，图(a)为端羟基聚丁二烯组分的实测值与预测值的线性关系图，图(b)为高能氧化剂组分的实测值与预测值的线性关系图，图(c)为RDX组分的实测值与预测值的线性关系图。
具体实施方式
结合图1，本发明一种固体推进剂制备过程的组分分散均匀性在线测试方法，包括以下步骤：
步骤1、准备不同配方的固体推进剂在制过程的固体样品，随机选取若干固体样品作为样品集，在一定的光谱测量条件下利用近红外光谱仪采集样品集的近红外光谱数据。其中，固体推进剂中含有聚丁二烯、硝化纤维素、聚叠氮缩水甘油醚、3-叠氮甲基-3-甲基环氧丁烷聚合物、3,3-双环氧丁烷聚合物中的一种或两种以上的复合物；光谱测量条件为：固体样品的平铺厚度为1cm～5cm，近红外光谱仪的扫描次数设定为5～20次。在一定的光谱测量条件下利用近红外光谱仪采集样品集的近红外光谱数据，具体为：
步骤1-1、光源照射平铺的固体样品，且使得固体样品上的光斑直径为20mm～100mm；所述光源为卤钨灯、InGaAs基LED灯、氙气灯及超连续激光光源中的一种；
步骤1-2、近红外光谱仪的检测器接收固体样品漫反射的光，并对其进行检测，获得波长范围为700nm～2500nm的近红外光谱数据；每个样品重复测试5～20次得到相应的近红外光谱数据组，近红外光谱数据以吸光度形式存储；所述光谱仪的检测器为InGaAs基阵列检测器。
步骤2、采用常规检测方法测得样品集中各组分的含量实测值。其中，常规检测方法为化学测试方法、色谱法、物理分离或元素分析方法中的一种或几种。
步骤3、建立组分含量的近红外校正预测模型。具体为：
步骤3-1、对步骤1采集到的近红外光谱数据进行预处理；其中，预处理采用的是多元散射校正法、矢量归一化法、标准正态变量变换法、最大-最小归一化法、Savitzky-Golay卷积平滑法、Savitzky-Golay一阶导数法、小波变换法和NCL法中的一种或多种方法；
步骤3-2、对步骤1获得的样品集进行划分，分为校正集和验证集；其中，对样品集进行划分采用的方法为随机法、Kennad-Stone法、SPXY法中的一种；
步骤3-3、将校正集样品的近红外光谱数据与步骤2获得的每个样品各组分的含量实测值相对应，通过变量压缩方法挑选最佳建模变量即最佳建模所用的波长变量；其中，变量压缩方法采用竞争自适应重加权法、无信息变量消除法、组合区间偏最小二乘法、组合区间偏最小二乘法、遗传算法中的一种或两种；
步骤3-4、通过化学计量学建模方法建立样品各组分的定量校正模型；其中，化学计量学建模方法采用偏最小二乘法、主成分回归、最小二乘支持向量机或人工神经网络法中的一种；
步骤3-5、利用验证集样品对步骤3-4建立的各组分的定量校正模型进行外部验证，获得组分含量的近红外校正预测模型；具体为：
步骤3-5-1、将步骤1得到的验证集样品的近红外光谱数据导入步骤3-4建立的各组分的定量校正模型中，获得验证集样品各组分的预测值；
步骤3-5-2、将步骤3-5-1得到的验证集样品各组分的预测值与步骤2测得的验证集样品各组分的含量实测值相比，进行模型准确度检验，即对组分的定量校正模型进行外部验证；具体为：
求取验证集样品各组分的预测值与对应的实测值之间的组分含量相对偏差，之后对验证集所有样品的组分含量相对偏差进行求平均，获得平均组分含量相对偏差；若平均组分含量相对偏差在设定的产品质量控制所需的范围内，则所述组分的定量校正模型可用，反之则重复步骤3，优化建模条件直到所述定量校正模型可用；最终获得的可用的组分的定量校正模型即为组分含量的近红外校正预测模型。
步骤4、取待测样品，按照步骤1中的光谱测量条件采集待测样品的近红外光谱数据，并导入步骤3建立的近红外校正预测模型中，经模型运算，获得待测样品各组分的含量预测值。
步骤5、根据步骤4获得的待测样品各组分含量预测值和步骤2中每个样品各组分的含量实测值，判断待测样品组分分散均匀性。具体为：
步骤5-1、根据步骤4获得的待测样品各组分含量预测值与步骤2中每个样品各组分含量实测值，求取待测样品各组分含量与步骤2中每个样品各组分含量的组分含量相对偏差；求取组分含量相对偏差所用公式为：
步骤5-2、对步骤5-1获得的所有组分含量相对偏差进行求平均，获得平均组分含量相对偏差；平均组分含量相对偏差越小，组分的分散均匀性越好，反之组分的分散均匀性越差。
步骤6、重复步骤4和步骤5，实现对所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
进一步地，步骤3中所述建立近红外校正预测模型过程中建立各组分的定量校正模型具体为：
步骤3-1’、采用多元散射校正和Savitzky-Golay卷积平滑法对近红外光谱数据进行预处理；
步骤3-2’、按照6:4～9:1的比例范围，采用Kennard-Stone法对步骤1获得的样品集进行划分，分为校正集和验证集；
步骤3-3’、将校正集样本的近红外光谱数据与步骤2中的各组分的含量实测值相对应，通过组合区间偏最小二乘法挑选最佳建模变量；
步骤3-4’、采用偏最小二乘法建立样品中各组分含量的定量校正模型，各组分含量的定量校正模型的矩阵表达式为：
Y＝BX+E
式中，Y为n个组分在m个样品中的浓度矩阵，其对应各组分的含量实测值，称为应变量；B为消光系数矩阵；X为m个样品中在波长点1到k的吸光度矩阵，称为自变量；E为残差矩阵；
其中，
进一步地，步骤3中所述建立近红外校正预测模型过程中建立各组分的定量校正模型具体为：
步骤3-1”、采用最大-最小归一化法和Savitzky-Golay一阶导数法对近红外光谱数据进行预处理；
步骤3-2”、按照7：3～9:1的比例范围，采用随机法对步骤1获得的样品集进行划分，分为校正集和验证集；
步骤3-3”、将校正集样本的近红外光谱数据与步骤2中的各组分的含量实测值相对应，通过竞争自适应重加权法挑选最佳建模变量；
步骤3-4”、采用偏最小二乘法建立样品中各组分含量的定量校正模型，各组分含量的定量校正模型的矩阵表达式为：
Y＝BX+E。
下面结合实施例对本发明作进一步详细描述。
实施例1
选用含乙基纤维素的固体推进剂进行测试，制备不同配方由含硝化纤维素(NC)、含能增塑剂(B)、高能固体填料(C)三种主成分混制而成的不规则片状样品，样品的颗粒直径为10mm～30mm。下面是具体的测试方法和过程：
1)准备12组不同配方的固体推进剂在制过程的固体样品，随机选取8组样品作为样品集，利用近红外光谱仪采集样品集的近红外光谱数据；采用漫反射方式测试光谱数据，样品装到样品杯中，平铺厚度为1cm～2cm，漫反射光由检测器进行检测，扫描次数设定为20次，波长范围为1100nm～2300nm，每个样品重复测试20次，得到160条光谱数据以吸光度形式存储。采集近红外吸收光谱所用光谱仪的检测器为InGaAs基阵列检测器，测试所用的光源为卤钨灯，光源输出到测试样品上的光斑直径范围为100mm。图2为该组样品测试的近红外光谱数据原始吸光度曲线汇总图。
2)采用常规化学分析测试方法测得样品集中各组分的含量实测值，得到三种组分的含量如下表1所示。
表1样品集组分含量测试结果
3)将步骤1测试得到的160条吸光度曲线，按照7:1的比例采用Kennard-Stone法进行样品集划分，分为校正集和验证集，并将步骤1)中得到的校正集中的近红外光谱数据与步骤2)中的组分实际测试结果相对应，采用多元散射校正和Savitzky-Golay卷积平滑对测试的光谱数据进行预处理；变量压缩方法选用组合区间偏最小二乘法，并通过偏最小二乘法化学计量学建模方法建立样品中各组分含量的定量校正模型，用验证集样品对所建立的定量校正模型进行外部验证，最终得到组分含量的近红外校正预测模型。
对建立的近红外定量校正模型进行外部验证的步骤：将步骤1)得到的验证集各样品的近红外光谱数据导入所建立的定量校正模型中，得到验证集样品各组分的预测值，将其与步骤2)中所述验证集中各样品组分的含量实测值相比，进行模型准确度检验，预测值(％)与实测值(％)平均相对偏差的阈值设定为1.0％，实际测试结果为0.75％，则所得定量校正模型可用。三种组分含量的近红外预测模型效果图如图3所示，从图可知三种组分的近红外预测值与实测值均匀分布在对角线上，说明模型能够对组分含量进行准确预测。各样品的主成分预测结果及偏差数据如下表2所示。
表2样品组分含量及均匀性验证测试结果
4)选取组分含量在预定范围的其他待测样品，按照步骤1)中的光谱测量条件采集样品的近红外光谱数据，导入校正预测模型中，经模型运算，得到主要组分的含量预测值。
5)由步骤4)得到的组分含量预测值与步骤2)中的组分含量实测值的计算组分含量相对偏差，所述的组分含量相对偏差为含量预测值(％)与含量实测值(％)差值的绝对值与含量实测值(％)之比，采用各组分含量相对偏差的平均值大小来判断组分的分散均匀性，组分含量相对偏差为0.88％，说明组分的分散均匀性较好。
6)重复步骤4和步骤5)，实现所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
实施例2
采用含聚叠氮缩水甘油醚、高能氧化剂和RDX为主要成分的配方体系，在物料混合过程进行在线测试。下面是具体的建模测试方法和过程：
1)准备16组不同配方的固体推进剂在制过程的固体样品，随机选取8组样品作为样品集，利用近红外光谱仪采集样品集的近红外光谱数据；采用漫反射方式测试光谱数据，样品装到样品杯中，平铺厚度为1cm～2cm，漫反射光由光谱仪的检测器进行检测，扫描次数设定为20次，波长范围为1000nm～1800nm，每个样品重复测试20次，得到160条光谱数据以吸光度形式存储。采集近红外吸收光谱所用光谱仪的检测器为InGaAs基阵列检测器，测试所用的光源为卤钨灯，光源输出到测试样品上的光斑直径范围为20mm。图4为本实施例固体推进剂物料的近红外光谱数据汇总曲线图。
2)采用常规化学分析测试方法测得样品集中各组分的含量实测值。测试结果如下表3所示。
表3样品集组分含量测试结果
3)将步骤1测试得到的160条吸光度曲线，按照5:1的比例采用随机法进行样品集划分，分为校正集和验证集，并将步骤1)中得到的校正集中的近红外光谱数据与步骤2)中的组分实际测试结果相对应，采用最大-最小归一化法和Savitzky-Golay一阶导数法对测试的光谱数据进行预处理；变量压缩方法选用竞争自适应重加权法，并通过偏最小二乘法化学计量学建模方法建立样品中各组分含量的定量校正模型，用验证集样品对所建立的定量校正模型进行外部验证，最终得到组分含量的近红外校正预测模型。
对建立的近红外定量校正模型进行外部验证的步骤：将步骤1)得到的验证集各样品的近红外光谱数据导入所建立的定量校正模型中，得到验证集样品各组分的预测值，将其与步骤2)中所述验证集中各样品组分的含量实测值相比，进行模型准确度检验，预测值(％)与实测值(％)平均相对偏差的阈值设定为1.0％，实际测试结果为0.50％，则所得定量校正模型可用。三种组分含量的近红外预测模型效果图如图5所示。
4)选取组分含量在预定范围的其他待测样品，按照步骤1)中的光谱测量条件采集样品的近红外光谱数据，导入校正预测模型中，经模型运算，得到主要组分的含量预测值。
5)由步骤4)得到的组分含量预测值与步骤2)中的组分含量实测值的计算组分含量相对偏差，所述的组分含量相对偏差为含量预测值(％)与含量实测值(％)差值的绝对值与含量实测值(％)之比，采用各组分相对偏差的平均值大小来判断组分的分散均匀性，组分含量相对偏差为0.45％，说明组分的分散均匀性较好。
6)重复步骤4和步骤5)，实现所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
实施例3
采用含聚丁二烯粘合剂、高能氧化剂和RDX为主要成分的配方体系，在物料混合过程进行在线测试。下面是具体的建模测试方法和过程：
1)准备12组不同配方的固体推进剂在制过程的固体样品，随机选取8组样品作为样品集，利用近红外光谱仪采集样品集的近红外光谱数据；采用漫反射方式测试光谱数据，样品装到样品杯中，平铺厚度为2cm～3cm，漫反射光由光谱仪的检测器进行检测，扫描次数设定为20次，波长范围为900nm～1700nm，每个样品重复测试20次，得到160条光谱数据以吸光度形式存储。采集近红外吸收光谱所用光谱仪的检测器为InGaAs基阵列检测器，测试所用的光源为卤钨灯，光源输出到测试样品上的光斑直径范围为50mm。
2)采用常规化学分析测试方法测得样品集中各组分的含量实测值。
测试结果如下表4所示：
表4样品集组分含量测试结果
| 样品编号 | 聚丁二烯粘合剂 | 高能氧化剂 | RDX |
| 1 | 31.7 | 54.9 | 12.6 |
| 2 | 23.8 | 56.7 | 19.0 |
| 3 | 23.8 | 47.1 | 28.5 |
| 4 | 22.2 | 41.8 | 35.5 |
| 5 | 19.6 | 40.7 | 39.2 |
| 6 | 17.0 | 42.1 | 40.6 |
| 7 | 19.4 | 33.6 | 46.5 |
| 8 | 18.6 | 27.4 | 53.6 |
3)将步骤1测试得到的160条吸光度曲线，按照4:1的比例采用随机法进行样品集划分，分为校正集和验证集，并将步骤1)中得到的校正集中的近红外光谱数据与步骤2)中的组分实际测试结果相对应，采用测试得到的吸光度原始数据进行建模分析；通过偏最小二乘法化学计量学建模方法建立样品中各组分含量的定量校正模型，用验证集样品对所建立的定量校正模型进行外部验证，最终得到组分含量的近红外校正预测模型。
对建立的近红外定量校正模型进行外部验证的步骤：将步骤1)得到的验证集各样品的近红外光谱数据导入所建立的定量校正模型中，得到验证集样品各组分的预测值，将其与步骤2)中所述验证集中各样品组分的含量实测值相比，进行模型准确度检验，预测值(％)与实测值(％)平均相对偏差的阈值设定为1.0％，实际测试结果为0.82％，则所得定量校正模型可用。
4)采用组分含量在预定范围的其他待测样品，按照步骤1)中的光谱测量条件采集样品的近红外光谱数据，导入校正预测模型中，经模型运算，得到主要组分的含量预测值。
5)由步骤4)得到的组分含量预测值与步骤2)中的组分含量实测值的计算组分含量相对偏差，所述的组分含量相对偏差为含量预测值(％)与含量实测值(％)差值的绝对值与含量实测值(％)之比，采用各组分相对偏差的平均值大小来判断组分的分散均匀性，组分含量相对偏差为0.75％，说明组分的分散均匀性较好。
6)重复步骤4和步骤5)，实现所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
实施例4
采用与实施例1相同配比的样品进行测试，样品的颗粒直径为30mm～50mm。下面是具体的测试方法和过程：
1)准备20组不同配方的固体推进剂在制过程的固体样品，随机选取8组样品作为样品集，利用近红外光谱仪采集样品集的近红外光谱数据；采用漫反射方式测试光谱数据，样品装到样品杯中，平铺厚度为2cm～3cm，漫反射光由光谱仪的检测器进行检测，扫描次数设定为20次，波长范围为900nm～1700nm，每个样品重复测试20次，得到160条光谱数据以吸光度形式存储。采集近红外吸收光谱所用光谱仪的检测器为InGaAs基阵列检测器，测试所用的光源为卤钨灯，光源输出到测试样品上的光斑直径范围为50mm。
2)采用常规化学分析测试方法测得样品集中各组分的含量实测值。
3)将步骤1测试得到的160条吸光度曲线，按照5:1的比例采用Kennard-Stone法进行样品集划分，分为校正集和验证集，并将步骤1)中得到的校正集中的近红外光谱数据与步骤2)中的组分实际测试结果相对应，采用矢量归一化法对测试的光谱数据进行预处理；变量压缩方法选用组合区间偏最小二乘法，并通过偏最小二乘法化学计量学建模方法建立样品中各组分含量的定量校正模型，用验证集样品对所建立的定量校正模型进行外部验证，最终得到组分含量的近红外校正预测模型。
对建立的近红外定量校正模型进行外部验证的步骤：将步骤1)得到的验证集各样品的近红外光谱数据导入所建立的定量校正模型中，得到验证集样品各组分的预测值，将其与步骤2)中所述验证集中各样品组分的含量实测值相比，进行模型准确度检验，预测值(％)与实测值(％)平均相对偏差的阈值设定为1.0％，实际测试结果为0.72％，则所得定量校正模型可用。
4)采用组分含量在预定范围的其他待测样品，按照步骤1)中的光谱测量条件采集样品的近红外光谱数据，导入校正预测模型中，经模型运算，得到主要组分的含量预测值。
5)由步骤4)得到的组分含量预测值与步骤2)中的组分含量实测值的计算组分含量相对偏差，所述的组分含量相对偏差为含量预测值(％)与含量实测值(％)差值的绝对值与含量实测值(％)之比，采用各组分相对偏差的平均值大小来判断组分的分散均匀性，组分含量相对偏差为0.78％，说明组分的分散均匀性较好。
6)重复步骤4和步骤5)，实现所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
实施例5
采用含聚叠氮缩水甘油醚、3-叠氮甲基-3-甲基环氧丁烷聚合物和高能填料为主要成分的配方体系，在物料混合过程进行在线测试。下面是具体的建模测试方法和过程：
1)准备15组不同配方的固体推进剂在制过程的固体样品，随机选取10组样品作为样品集，利用近红外光谱仪采集样品集的近红外光谱数据；采用漫反射方式测试光谱数据，样品装到样品杯中，平铺厚度为1cm～2cm，漫反射光由光谱仪的检测器进行检测，扫描次数设定为20次，波长范围为900nm～1700nm，每个样品重复测试20次，得到200条光谱数据以吸光度形式存储。采集近红外吸收光谱所用光谱仪的检测器为InGaAs基阵列检测器，测试所用的光源为卤钨灯，光源输出到测试样品上的光斑直径范围为50mm。
2)采用常规化学分析测试方法测得样品集中各组分的含量实测值。
3)将步骤1测试得到的200条吸光度曲线，按照5:1的比例采用SPXY法进行样品集划分，分为校正集和验证集，并将步骤1)中得到的校正集中的近红外光谱数据与步骤2)中的组分实际测试结果相对应，采用最大-最小归一化法对测试的光谱数据进行预处理；变量压缩方法选用竞争自适应重加权法，并通过偏最小二乘法化学计量学建模方法建立样品中各组分含量的定量校正模型，用验证集样品对所建立的定量校正模型进行外部验证，最终得到组分含量的近红外校正预测模型。
对建立的近红外定量校正模型进行外部验证的步骤：将步骤1)得到的验证集各样品的近红外光谱数据导入所建立的定量校正模型中，得到验证集样品各组分的预测值，将其与步骤2)中所述验证集中各样品组分的含量实测值相比，进行模型准确度检验，预测值(％)与实测值(％)平均相对偏差的阈值设定为1.0％，实际测试结果为0.50％，则所得定量校正模型可用。
4)采用组分含量在预定范围的其他待测样品，按照步骤1)中的光谱测量条件采集样品的近红外光谱数据，导入校正预测模型中，经模型运算，得到主要组分的含量预测值。
5)由步骤4)得到的组分含量预测值与步骤2)中的组分含量实测值的计算组分含量相对偏差，所述的组分含量相对偏差为含量预测值(％)与含量实测值(％)差值的绝对值与含量实测值(％)之比，采用各组分相对偏差的平均值大小来判断组分的分散均匀性，组分含量相对偏差为0.45％，说明组分的分散均匀性较好。
6)重复步骤4和步骤5)，实现所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
通过以上实施例可以看出，本发明以近红外光谱模型结合实时测试的光谱数据实现固体推进剂生产过程中组分含量及分散均匀性在线测试，相比于现有的检测方法所存在的测试周期长、操作复杂等缺陷，具有方便、准确、环境适应性强等优点。
Claims (10)
1.一种固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，包括以下步骤：
步骤1、准备不同配方的固体推进剂在制过程的固体样品，随机选取若干固体样品作为样品集，在一定的光谱测量条件下利用近红外光谱仪采集样品集的近红外光谱数据；
步骤2、采用常规检测方法测得样品集中每个样品各组分的含量实测值；
步骤3、建立组分含量的近红外校正预测模型；
步骤4、取待测样品，按照步骤1中的光谱测量条件采集待测样品的近红外光谱数据，并导入步骤3建立的近红外校正预测模型中，经模型运算，获得待测样品各组分的含量预测值；
步骤5、根据步骤4获得的待测样品各组分含量预测值和步骤2中每个样品各组分的含量实测值，判断待测样品组分分散均匀性；
步骤6、重复步骤4和步骤5，实现对所测试的样品组分含量和组分分散均匀性的在线测试分析和实时显示。
2.根据权利要求1所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤1所述固体推进剂中含有聚丁二烯、硝化纤维素、聚叠氮缩水甘油醚、3-叠氮甲基-3-甲基环氧丁烷聚合物、3,3-双环氧丁烷聚合物中的一种或两种以上的复合物。
3.根据权利要求1所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤1中所述光谱测量条件为：固体样品的平铺厚度为1cm～5cm，近红外光谱仪的扫描次数设定为5～20次；所述在一定的光谱测量条件下利用近红外光谱仪采集样品集的近红外光谱数据，具体为：
步骤1-1、光源照射平铺的固体样品，且使得固体样品上的光斑直径为20mm～100mm；所述光源为卤钨灯、InGaAs基LED灯、氙气灯及超连续激光光源中的一种；
步骤1-2、近红外光谱仪的检测器接收固体样品漫反射的光，并对其进行检测，获得波长范围为700nm～2500nm的近红外光谱数据；每个样品重复测试5～20次得到相应的近红外光谱数据组，近红外光谱数据以吸光度形式存储；所述光谱仪的检测器为InGaAs基阵列检测器。
4.根据权利要求1所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤2中所述常规检测方法为化学测试方法或物理分离中的一种或几种。
5.根据权利要求1所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤3所述建立组分含量的近红外校正预测模型具体为：
步骤3-1、对步骤1采集到的近红外光谱数据进行预处理；
步骤3-2、对步骤1获得的样品集进行划分，分为校正集和验证集；
步骤3-3、将校正集样品的近红外光谱数据与步骤2获得的每个样品各组分的含量实测值相对应，通过变量压缩方法挑选最佳建模变量即最佳建模所用的波长变量；
步骤3-4、通过化学计量学建模方法建立样品各组分的定量校正模型；
步骤3-5、利用验证集样品对步骤3-4建立的各组分的定量校正模型进行外部验证，获得组分含量的近红外校正预测模型。
6.根据权利要求5所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，
(1)步骤3-1所述对近红外光谱数据进行预处理采用的是多元散射校正法、矢量归一化法、标准正态变量变换法、最大-最小归一化法、Savitzky-Golay卷积平滑法、Savitzky-Golay一阶导数法、小波变换法和NCL法中的一种或多种方法；
(2)步骤3-2所述对样品集进行划分采用的方法为随机法、Kennad-Stone法、SPXY法中的一种；
(3)步骤3-3所述变量压缩方法采用竞争自适应重加权法、无信息变量消除法、组合区间偏最小二乘法、组合区间偏最小二乘法、遗传算法中的一种或两种；
(4)步骤3-4所述化学计量学建模方法采用偏最小二乘法、主成分回归、最小二乘支持向量机或人工神经网络法中的一种。
7.根据权利要求5或6所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤3中所述建立近红外校正预测模型过程中建立各组分的定量校正模型具体为：
步骤3-1’、采用多元散射校正和Savitzky-Golay卷积平滑法对近红外光谱数据进行预处理；
步骤3-2’、按照6:4～9:1的比例范围，采用Kennard-Stone法对步骤1获得的样品集进行划分，分为校正集和验证集；
步骤3-3’、将校正集样本的近红外光谱数据与步骤2中的各组分的含量实测值相对应，通过组合区间偏最小二乘法挑选最佳建模变量；
步骤3-4’、采用偏最小二乘法建立样品中各组分含量的定量校正模型，各组分含量的定量校正模型的矩阵表达式为：
Y＝BX+E
式中，Y为n个组分在m个样品中的浓度矩阵，其对应各组分的含量实测值，称为应变量；B为消光系数矩阵；X为m个样品中在波长点1到k的吸光度矩阵，称为自变量；E为残差矩阵；
其中，
8.根据权利要求7所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤3中所述建立近红外校正预测模型过程中建立各组分的定量校正模型具体为：
步骤3-1”、采用最大-最小归一化法和Savitzky-Golay一阶导数法对近红外光谱数据进行预处理；
步骤3-2”、按照7：3～9:1的比例范围，采用随机法对步骤1获得的样品集进行划分，分为校正集和验证集；
步骤3-3”、将校正集样本的近红外光谱数据与步骤2中的各组分的含量实测值相对应，通过竞争自适应重加权法挑选最佳建模变量；
步骤3-4”、采用偏最小二乘法建立样品中各组分含量的定量校正模型，各组分含量的定量校正模型的矩阵表达式为：
Y＝BX+E。
9.根据权利要求5所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤3-5所述利用验证集样品对步骤3-4建立的定量校正模型进行外部验证，获得组分含量的近红外校正预测模型，具体为：
步骤3-5-1、将步骤1得到的验证集样品的近红外光谱数据导入步骤3-4建立的各组分的定量校正模型中，获得验证集样品各组分的预测值；
步骤3-5-2、将步骤3-5-1得到的验证集样品各组分的预测值与步骤2测得的验证集样品各组分的含量实测值相比，进行模型准确度检验，即对组分的定量校正模型进行外部验证；具体为：
求取验证集样品各组分的预测值与对应的实测值之间的组分含量相对偏差，之后对验证集所有样品的组分含量相对偏差进行求平均，获得平均组分含量相对偏差；若平均组分含量相对偏差在设定的产品质量控制所需的范围内，则所述组分的定量校正模型可用，反之则重复步骤3，优化建模条件直到所述定量校正模型可用；最终获得的可用的组分的定量校正模型即为组分含量的近红外校正预测模型。
10.根据权利要求1所述的固体推进剂制备过程的组分分散均匀性在线测试方法，其特征在于，步骤5具体为：
步骤5-1、根据步骤4获得的待测样品各组分含量预测值与步骤2中每个样品各组分含量实测值，求取待测样品各组分含量与步骤2中每个样品各组分含量的组分含量相对偏差；所述求取组分含量相对偏差所用公式为：
步骤5-2、对步骤5-1获得的所有组分含量相对偏差进行求平均，获得平均组分含量相对偏差；平均组分含量相对偏差越小，组分的分散均匀性越好，反之组分的分散均匀性越差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810768568.6A CN109142263B (zh) | 2018-07-13 | 2018-07-13 | 一种固体推进剂制备过程的组分分散均匀性在线测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810768568.6A CN109142263B (zh) | 2018-07-13 | 2018-07-13 | 一种固体推进剂制备过程的组分分散均匀性在线测试方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109142263A CN109142263A (zh) | 2019-01-04 |
| CN109142263B true CN109142263B (zh) | 2020-12-18 |
Family
ID=64800543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810768568.6A Active CN109142263B (zh) | 2018-07-13 | 2018-07-13 | 一种固体推进剂制备过程的组分分散均匀性在线测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109142263B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110411978B (zh) * | 2019-07-26 | 2020-05-22 | 中国农业大学 | 一种物料均匀性检测设备及检测方法 |
| CN110715902A (zh) * | 2019-08-29 | 2020-01-21 | 超威电源有限公司 | 一种铅酸蓄电池正铅膏均匀性的检测方法 |
| CN111426715B (zh) * | 2020-05-15 | 2023-03-31 | 云南中烟工业有限责任公司 | 一种基于xps一步法测量烟丝加香均匀性的方法 |
| CN111487217A (zh) * | 2020-05-25 | 2020-08-04 | 西安近代化学研究所 | 浇注pbx炸药药浆固化剂混合均匀性无损检测方法 |
| CN115096943B (zh) * | 2022-06-21 | 2024-04-26 | 郑州磨料磨具磨削研究所有限公司 | 一种流体磨具均匀性的无损检测装置及测试方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419277A (zh) * | 2011-08-16 | 2012-04-18 | 中国检验检疫科学研究院 | 防晒化妆品中纳米二氧化钛的识别检测方法 |
| UA86830U (uk) * | 2013-07-23 | 2014-01-10 | Львовский Национальный Университет Имени Ивана Франка | Спосіб визначення розподілу за розмірами частинок металу у дисперсійних системах |
| CN203941107U (zh) * | 2014-06-11 | 2014-11-12 | 华南理工大学 | 一种聚合物熔体性质的拉曼光谱在线测量装置 |
| CN203981581U (zh) * | 2014-06-11 | 2014-12-03 | 华南理工大学 | 一种聚合物熔体性质的近红外光谱在线测量装置 |
| CN104964947A (zh) * | 2015-06-12 | 2015-10-07 | 安徽华润金蟾药业股份有限公司 | 一种黄芪配方颗粒混合过程在线质量控制方法 |
| CN105699258A (zh) * | 2016-01-28 | 2016-06-22 | 华侨大学 | 一种细骨料在线检测装置及方法 |
-
2018
- 2018-07-13 CN CN201810768568.6A patent/CN109142263B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419277A (zh) * | 2011-08-16 | 2012-04-18 | 中国检验检疫科学研究院 | 防晒化妆品中纳米二氧化钛的识别检测方法 |
| UA86830U (uk) * | 2013-07-23 | 2014-01-10 | Львовский Национальный Университет Имени Ивана Франка | Спосіб визначення розподілу за розмірами частинок металу у дисперсійних системах |
| CN203941107U (zh) * | 2014-06-11 | 2014-11-12 | 华南理工大学 | 一种聚合物熔体性质的拉曼光谱在线测量装置 |
| CN203981581U (zh) * | 2014-06-11 | 2014-12-03 | 华南理工大学 | 一种聚合物熔体性质的近红外光谱在线测量装置 |
| CN104964947A (zh) * | 2015-06-12 | 2015-10-07 | 安徽华润金蟾药业股份有限公司 | 一种黄芪配方颗粒混合过程在线质量控制方法 |
| CN105699258A (zh) * | 2016-01-28 | 2016-06-22 | 华侨大学 | 一种细骨料在线检测装置及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于NIR光谱的固体推进剂代用料混合特性;程士超等;《固体火箭技术》;20161231;第39卷(第3期);第373-377页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109142263A (zh) | 2019-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109142263B (zh) | 一种固体推进剂制备过程的组分分散均匀性在线测试方法 | |
| CN103018195A (zh) | 近红外光谱测定pbx炸药中pctfe含量的方法 | |
| CN108732127B (zh) | 一种检测烟丝中各组分掺配比例的方法 | |
| CN109540837B (zh) | 近红外快速检测苎麻叶片木质纤维素含量的方法 | |
| CN101271063A (zh) | 气体红外光谱分析支持向量机校正模型方法的流程 | |
| CN108020526A (zh) | 一种丁羟推进剂药浆组分近红外检测方法 | |
| CN108444976B (zh) | 一种基于拉曼光谱的天然气热值测量方法 | |
| CN104596979A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶纤维素的方法 | |
| CN104596975A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶木质素的方法 | |
| WO2020248961A1 (zh) | 一种无参考值的光谱波数选择方法 | |
| CN111579523A (zh) | 一种利用近红外光谱法快速检测孕马尿中结合雌激素含量的方法 | |
| CN112179871A (zh) | 一种酱类食品中己内酰胺含量无损检测的方法 | |
| CN103760130B (zh) | 近红外光谱测定复方麝香注射液中吐温-80含量的方法 | |
| CN104596976A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶蛋白质的方法 | |
| CN104316492A (zh) | 近红外光谱测定马铃薯块茎中蛋白质含量的方法 | |
| CN106872398A (zh) | 一种hmx炸药水分含量快速测量方法 | |
| Yan et al. | Rapid detection of cAMP content in red jujube using near-infrared spectroscopy | |
| AU2017222946A1 (en) | A method and apparatus for the detection of the presence of mycotoxins in cereals | |
| CN106645012A (zh) | 成品汽柴油中酯类化合物快速定量分析方法 | |
| CN110865044A (zh) | 一种鉴别掺白油有机硅产品的光谱分析方法 | |
| CN115436315A (zh) | 一种基于近红外光谱的水泥外加剂检测方法 | |
| Dong et al. | Nondestructive method for analysis of the soybean quality | |
| CN110174362B (zh) | 一种检测中性糖和酸性糖含量的方法 | |
| CN112595691A (zh) | 基于近红外拉曼光谱融合的薰衣草精油特征组分定量分析模型的建立方法和定量分析方法 | |
| CN110702635A (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 |