CN108020526A - 一种丁羟推进剂药浆组分近红外检测方法 - Google Patents
一种丁羟推进剂药浆组分近红外检测方法 Download PDFInfo
- Publication number
- CN108020526A CN108020526A CN201711326992.7A CN201711326992A CN108020526A CN 108020526 A CN108020526 A CN 108020526A CN 201711326992 A CN201711326992 A CN 201711326992A CN 108020526 A CN108020526 A CN 108020526A
- Authority
- CN
- China
- Prior art keywords
- mrow
- sample
- htpb
- mean square
- root mean
- 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.)
- Pending
Links
- 229920002121 Hydroxyl-terminated polybutadiene Polymers 0.000 title claims abstract description 87
- 239000002002 slurry Substances 0.000 title claims abstract description 53
- 239000003814 drug Substances 0.000 title claims abstract description 48
- 239000003380 propellant Substances 0.000 title claims abstract description 45
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 239000011230 binding agent Substances 0.000 claims abstract description 53
- 239000007800 oxidant agent Substances 0.000 claims abstract description 49
- 239000004014 plasticizer Substances 0.000 claims abstract description 49
- 230000001590 oxidative effect Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 36
- 239000000470 constituent Substances 0.000 claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 15
- 238000002835 absorbance Methods 0.000 claims abstract description 12
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 74
- 238000012937 correction Methods 0.000 claims description 35
- 238000010200 validation analysis Methods 0.000 claims description 34
- 239000003795 chemical substances by application Substances 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000011156 evaluation Methods 0.000 claims description 11
- 238000012795 verification Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 4
- 229910052760 oxygen Inorganic materials 0.000 claims description 4
- 239000001301 oxygen Substances 0.000 claims description 4
- 238000003825 pressing Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 9
- 125000002887 hydroxy group Chemical group [H]O* 0.000 abstract description 6
- 239000000843 powder Substances 0.000 abstract description 4
- 238000009659 non-destructive testing Methods 0.000 abstract description 3
- 239000002699 waste material Substances 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- VJHINFRRDQUWOJ-UHFFFAOYSA-N dioctyl sebacate Chemical compound CCCCC(CC)COC(=O)CCCCCCCCC(=O)OCC(CC)CCCC VJHINFRRDQUWOJ-UHFFFAOYSA-N 0.000 description 35
- 238000004458 analytical method Methods 0.000 description 10
- 239000004449 solid propellant Substances 0.000 description 10
- 239000002131 composite material Substances 0.000 description 7
- 239000007787 solid Substances 0.000 description 6
- 239000000203 mixture Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000004615 ingredient Substances 0.000 description 3
- 238000003908 quality control method Methods 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000003647 oxidation Effects 0.000 description 2
- 238000007254 oxidation reaction Methods 0.000 description 2
- 238000000985 reflectance spectrum Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- ZZIZZTHXZRDOFM-XFULWGLBSA-N tamsulosin hydrochloride Chemical group [H+].[Cl-].CCOC1=CC=CC=C1OCCN[C@H](C)CC1=CC=C(OC)C(S(N)(=O)=O)=C1 ZZIZZTHXZRDOFM-XFULWGLBSA-N 0.000 description 2
- GDDNTTHUKVNJRA-UHFFFAOYSA-N 3-bromo-3,3-difluoroprop-1-ene Chemical compound FC(F)(Br)C=C GDDNTTHUKVNJRA-UHFFFAOYSA-N 0.000 description 1
- 102000003712 Complement factor B Human genes 0.000 description 1
- 108090000056 Complement factor B Proteins 0.000 description 1
- 241000871495 Heeria argentea Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- LNDKRVOOYMEYTC-UHFFFAOYSA-N buta-1,3-dien-1-ol Chemical compound OC=CC=C LNDKRVOOYMEYTC-UHFFFAOYSA-N 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009477 glass transition Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- 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/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)
- 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
一种丁羟推进剂药浆组分近红外检测方法，采用近红外光谱测试技术，通过选择可旋转的漫反射采样技术和多点采集技术解决因丁羟药浆样品特性对测试精度的影响，通过化学计量法、分析波段选择、最佳变量数等最佳建模参数确定，建立吸光度与组分含量之间关系，通过建立准确的校正模型，形成定量组分含量方法，根据所建立的吸光度与含量之间的相关关系，通过一张被测的样品的近红外光谱即可同时给出待测样品的多组分含量定量结果，实现丁羟推进剂药浆中粘结剂HTPB、增塑剂DOS、氧化剂AP含量快速检测，具有快速、高效、无损检测的特点，能够及时发现装药生产过程中药浆组分浓度变化，实现对工艺质量稳定性和均匀性的控制，避免不必要浪费。
Description
技术领域
本发明涉及推进剂组分装药检测方法，具体是一种丁羟推进剂药浆组分近红外快速检测方法,用于对固体推进剂中的粘结剂HTPB、增塑剂DOS和氧化剂AP含量进行快速检测,为发动机装药质量提供重要保障。
背景技术
推进剂药浆组分含量是固体发动机装药工艺过程重要控制的环节，它直接影响推进剂的各项性能指标。粘结剂含量波动大小对推进剂力学性能具有决定性的影响。增塑剂不但可以降低未固化推进剂药浆的粘度，降低推进剂玻璃化温度，而且其含量波动大小对改善低温力学性能有重要作用。氧化剂对体系燃烧性能至关重要。因此将混合药浆中粘结剂(HTPB)、增塑剂(KZ)、氧化剂(AP)组分含量控制在要求的技术指标范围内是发动机装药质量的重要保证条件。
对推进剂药浆组分实施快速检测是各国推进剂制造行业共同追求的目标。复合固体推进剂药浆是含有粘合剂基体，氧化剂颗粒，金属粒子以及其它一些功能组分的复杂浆状混合物体系，药浆不透光，而且含有大量的固体颗粒。由于体系复杂，其混合过程的分析检测难度非常大。目前，国内在推进剂制造行业仅仅强调对所用原材料进行指标控制和分析，通过配方调试来确定装药参数。虽然这样做对控制发动机装药质量起到有效保障作用,但由于装药过程质量控制缺乏有效检测手段，因此当受环境等诸多因素影响时，有时的性能预示实验往往并不尽人意，浪费了大量时间、人力、物力。
目前，对于复合固体推进剂药浆组分含量的快速检测还处于简单的探索研究阶段。传统分析方法测定一个样品的多种性质或浓度数据需要多种分析仪器，成本高工作效率低，远不能适应推进剂组分快速检测的需要。因此，多年来国内外对固体推进剂装药过程质量控制快速分析虽有迫切生产需求，但在固体推进剂装药分析检测技术领域，一直均未实现对复合固体推进剂稀药浆组分含量的快速检测。
近红外分析技术具有分析快速、无损、无污、可进行多组分同时测试、易于实现在线的特点，它在农业、医药、化工等行业质量控制快速分析方面得到广泛应用，虽然从信号干扰严重、信号微弱的光谱中提取有效信息是近红外在推进剂应用的难点，但各国仍在推进剂领域积极开展近红外技术探索研究。Talanta报道，2004年，加拿大布里斯托尔宇航有限Rockwood推进剂厂Michael D.Judge利用近红外技术对固体推进剂中的四种组分均匀预混所组成的体系进行了组分含量测定。这种测定的研究体系相对简单，体系内仅有四种组分组成，与由十多种组分组成的真正固体推进剂相差甚远，且液体组分占～99％，可用透射方式进行样品光谱采集，因此其获得的光谱强度比采真正固体推进剂的光谱强度强很多。而本申请研究的却是由十多种组分组成的实际应用的固体推进剂，固体组分高达约90％的真正稀药浆。由于近红外技术属无损测试，故固体组分比例高达～90％导致的微观非均匀性、样品颜色较深使光谱强度很弱、多达十余个组分所组成体系形成了固体推进剂丁羟药浆样品特性的复杂性，形成了对模型预测精度的严重不利影响，使得固体推进剂丁羟药浆的近红外检测成为检测难点。
发明内容
为克服现有技术中存在的模型预测精度影响大的不足，本发明提出了一种丁羟推进剂药浆组分近红外检测方法。
本发明的具体过程是：
步骤1.建立样品库。
所述建立的样品库是丁羟推进剂药浆样品库，具体过程是采用实验室配制和生产现场取样两种方式获得，并将获得的样品分别用于建立校正集和建立验证集。
步骤2.采集各样品的近红外漫反射光谱。
用傅立叶变换红外光谱仪，通过旋转杯分别采集各样品的近红外漫反射光谱。
采集样品的近红外漫反射光谱时，将第一个样品分为三份，分别放入旋转杯内并平铺为1cm的厚度；用称量瓶轻压药浆以使样品紧贴杯底。分别旋转各装有样品的旋转杯，通过漫反射积分球装置对所述样品进行光谱采集，得到第一个样品的三张近红外漫反射光谱数据，取平均光谱用于数据处理。采集时，扫描次数为60～120次，分辨率为4cm-1，增益为2，近红外漫反射光谱采集范围为4000～10000cm-1。
重复所述获得第一个样品平均近红外漫反射光谱的过程，分别对其余各样品进行光谱测试；分别得到其余各样品的平均近红外漫反射光谱。
步骤3.滤波平滑预处理。
以得到的各样品的平均近红外漫反射光谱分别作为原始光谱。具体是，通过滤波对原始光谱进行平滑预处理，以消除各种干扰因素引起的光谱偏差，并对光谱图进行波长选择。波长选择时选择相关性最好的两个光谱区域。分别得到经过预处理后的各样品的平均近红外漫反射光谱。
所述滤波平滑预处理时，相关性最好的两个光谱区域分别为6800cm-1～5550cm-1光谱区域和5180cm-1～4000cm-1光谱区域。
步骤4.建立校正模型。所述的校正模型是丁羟推进剂药浆吸光度分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量之间的关系。
在建立校正模型前需确定最佳变量的数量。具有最小预测残差平方和的变量数就是模型的最佳变量数。通过公式(1)确定最佳变量数：
式中：PRESS是预测残差平方和，n是样品数；是交互验证预测值；yi是实际值。
将经过预处理后的各样品的平均近红外漫反射光谱分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量进行关联，用偏最小二乘法建立其中93个样品与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量的校正模型，形成校正集。
步骤5.校正模型的性能评价。
所述的校正模型性能包括：
光谱的吸光度分别与粘结剂HTPB、增塑剂DOS和氧化剂AP之间含量的相关系数Rk，其中R为相关系数，k＝1、2、3，分别表示粘结剂HTPB的相关系数、增塑剂DOS的相关系数和氧化剂AP相关系数；
粘结剂HTPB的校正误差均方根与验证误差均方根；
增塑剂DOS的校正误差均方根与验证误差均方根；
氧化剂AP的校正误差均方根与验证误差均方根。
在对校正模型进行性能评价时：
通过相关系数R分别评价粘结剂HTPB的预测值与实际值关系的线性程度、增塑剂DOS的预测值与实际值关系的线性程度、氧化剂AP的预测值与实际值关系的线性程度；分别得到粘结剂HTPB相关系数R1、增塑剂DOS的相关系数R2和氧化剂AP的相关系数R3。
分别计算粘结剂HTPB的校正误差均方根和验证误差均方根；
分别计算增塑剂DOS的校正误差均方根与验证误差均方根；
分别计算氧化剂AP的校正误差均方根与验证误差均方根。
当粘结剂HTPB的相关系数R1大于0.9、校正误差均方根小于0.2、验证误差均方根小于0.2时，该粘结剂HTPB校正模型满足要求；
当增塑剂DOS的相关系数R2大于0.9、校正误差均方根小于0.1、验证误差均方根小于0.1时，该增塑剂DOS校正模型满足要求；
当氧化剂AP的相关系数R3大于0.9、校正误差均方根小于0.3、验证误差均方根小于0.3时，该氧化剂AP校正模型满足要求；
在评价校正模型的性能时，通过公式(2)分别确定粘结剂HTPB1、增塑剂DOS和氧化剂AP的相关系数Rk：
该Rk值表示预测值与实际值关系的线性程度；Rk越接近1，则校正模型的预测值与实际值相关性越强。
通过公式(3)确定校正误差均方根RMSEC：
通过公式(4)确定验证误差均方根RMSEP：
上述各式中：-NIR数学模型预测值；
Ci-组分含量实际值；
-实际值的平均值；
n-建立模型所用的校正集样本数；
m-用于检验模型的验证集样本数。
步骤6.模型预测。
所述模型预测的具体过程是另采集一张样品，并得到该样品的近红外光谱。将得到的近红外光谱带入至所建立的丁羟推进剂药浆吸光度分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量之间关系的模型中，同时得到粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量。
本发明是建立一种对丁羟固体推进剂药浆粘结剂HTPB、增塑剂DOS、氧化剂AP组分的快速检测方法，为发动机装药质量提供重要保障。
本发明利用近红外测试方法，通过选择可旋转的漫反射采样技术和多点采集技术等手段解决因丁羟药浆样品特性对测试精度的不利影响，通过化学计量法、分析波段选择、最佳变量数等最佳建模参数确定，建立吸光度与组分含量之间关系，通过建立准确的校正模型，形成定量组分含量方法，根据所建立的吸光度与含量之间的相关关系，通过一张被测的样品的近红外光谱即可同时给出待测样品的多组分含量定量结果，实现丁羟推进剂药浆组分快速检测。
本发明建立了一种丁羟推进剂药浆组分近红外快速检测方法，实现了丁羟复合固体推进剂稀药浆中粘结剂HTPB、增塑剂DOS、氧化剂AP三组分含量检测,具有以下特点：
(1)测定速度极快，从得到样品到给出数据不到十分钟；
(2)分析效率高，通过被测样品的一张光谱即可进行多组分同时定量；
(3)待测样品无需预处理，不破坏样品，属无损测试，利于实现在线监测。
本发明为装药过程中丁羟推进剂药浆组分中的粘结剂、增塑剂、氧化剂组分含量提供了一种快速、高效、无损检测手段，可及时发现装药生产过程中药浆组分浓度变化，实现对工艺质量稳定性和均匀性的控制，避免不必要浪费。
附图说明
图1为实施例1中丁羟推进剂药浆HTPB含量实际值与NIR预测值之间的关系图,
图2为实施例1中丁羟推进剂药浆KZ含量实际值与NIR预测值之间的关系图。
图3为实施例1中丁羟推进剂药浆AP含量实际值与NIR预测值之间的关系图。
图4为实施例2中丁羟推进剂药浆HTPB含量实际值与NIR预测值之间的关系图。
图5为实施例2中丁羟推进剂药浆KZ含量实际值与NIR预测值之间的关系图。
图6为实施例2中丁羟推进剂药浆AP含量实际值与NIR预测值之间的关系图。
图7为本发明的流程图。
图中的1均为该实施例中校正样品在关系图中的分布示意。
具体实施方式
本发明采用近红外光谱测试技术对某配方的丁羟推进剂药浆组分粘结剂HTPB、增塑剂DOS、氧化剂AP含量快速检测。所述的粘结剂HTPB为端羟基聚丁二烯，所述的增塑剂DOS为癸二酸二异辛酯，所述的氧化剂AP为高氯酸铵。通过提取混合物光谱中与需要分析的组分含量相关的光谱信息进行关联，建立相关模型,利用该模型对未知混合物的光谱进行预测，即可得到混合物中需要分析的组分含量。
实施例1
步骤1.建立样品库。通过收集丁羟推进剂药浆样品，建立样品库。所述样品是按常规方法采用实验室配制和生产现场取样两种方式获得，并将获得的样品分别用于建立校正集和建立验证集。
本实施例中，样品集有101个样品，其中93个样品用来建立校正集，8个用来建立验证集。
步骤2.采集各样品的近红外漫反射光谱。用傅立叶变换红外光谱仪采集样品近红外漫反射光谱。
所用仪器为美国Thermo Nicolet公司的Nexus 670傅立叶变换近红外光谱仪、漫反射积分球装置、OMNIC信号采集软件及TQ analyst数据处理软件。
采集样品的近红外漫反射光谱时，将第一个样品分为三份，分别放入旋转杯内并平铺为1cm的厚度；用称量瓶轻压药浆以使样品紧贴杯底。分别旋转各装有样品的旋转杯，通过漫反射积分球装置对所述样品进行光谱采集，得到第一个样品的三张近红外漫反射光谱数据，取平均光谱用于数据处理。采集时，扫描次数为120次，分辨率为4cm-1，增益为2，近红外漫反射光谱采集范围为4000～10000cm-1。
重复所述获得第一个样品平均近红外漫反射光谱的过程，分别对其余各样品进行光谱测试；分别得到其余各样品的平均近红外漫反射光谱。共计得到101个样品的平均近红外漫反射光谱。
步骤3.滤波平滑预处理。
以得到的101个样品的平均近红外漫反射光谱分别作为原始光谱。通过Savitzky-Golay滤波对原始光谱进行平滑预处理，以消除各种干扰因素引起的光谱偏差，并对光谱图进行波长选择。波长选择时选择相关性最好的两个光谱区域，分别为6800cm-1～5550cm-1光谱区域和5180cm-1～4000cm-1光谱区域。分别得到经过预处理后的各样品的平均近红外漫反射光谱。
所选择的光谱区域中，样品光谱在7000cm-1以上时噪声较大并存在一定偏移，没有明显的吸收峰，且7400cm-1～6800cm-1和5550cm-1～5180cm-1段因受水汽吸收影响比较明显，这些波段在分析时均不予使用。
通过波长选择确定的两个光谱区域，对减少噪音信号的影响、提高运算速度和模型的稳定性是有益的。因为当采用全谱区计算时，不仅计算量大，而且在某些光谱区域，样品的光谱信息很弱，或与样品的组成、性质间缺乏相关关系，引入这些波长会造成模型精度降低甚至错误，所以在光谱分析中，需对波长点进行选择，删除一些对建模无用的干扰波长。
步骤4.建立校正模型。
所述的校正模型是丁羟推进剂药浆吸光度分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量之间的关系。
定量分析的模型质量还取决于正确选取最佳变量的数量。选择因子数太小会出现拟合不足，丢掉了样品光谱信息，但是选择最佳变量太多会导致过拟合，引入噪音，降低模型质量。因此，建立定量校正模型时，需对最佳变量数进行合理选择。
所述对最佳变量数进行合理选择的方法采用陆婉珍等在《现代近红外光谱分析技术》中提出的采用交互验证法所得的预测残差平方和的方法确定。所述交互验证的方法是每次剔除一个或多个样品，然后用剩下的样品建立模型，用该模型预测被剔除样品，重复上述过程直到每个样品都被剔除过一次。具有最小预测残差平方和的变量数就是模型的最佳变量数。最佳变量数选择为12。通过公式(1)计算PRESS值：
式中n——样品数；——交互验证预测值；yi——实际值。
本实施例中，通过交互验证的方法最终确定建模所需的最佳变量数为12。
将经过预处理后的各样品的平均近红外漫反射光谱分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量进行关联，用偏最小二乘法建立其中93个样品与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量的校正模型，形成校正集。所述93个样品标记为i，i＝1～n；本实施例中，n＝93。
剩余的8个样品作为验证集，标记为j，j＝1～m；本实施例中，m＝8。
步骤5.校正模型的性能评价。
所述的校正模型性能包括：
光谱的吸光度分别与粘结剂HTPB、增塑剂DOS和氧化剂AP之间含量的相关系数Rk，其中R为相关系数，k＝1、2、3，分别表示粘结剂HTPB的相关系数、增塑剂DOS的相关系数和氧化剂AP相关系数。
粘结剂HTPB的校正误差均方根与验证误差均方根；
增塑剂DOS的校正误差均方根与验证误差均方根；
氧化剂AP的校正误差均方根与验证误差均方根。
在对校正模型进行性能评价时：
通过相关系数R分别评价粘结剂HTPB的预测值与实际值关系的线性程度、增塑剂DOS的预测值与实际值关系的线性程度、氧化剂AP的预测值与实际值关系的线性程度。分别得到粘结剂HTPB相关系数R1、增塑剂DOS的相关系数R2和氧化剂AP的相关系数R3。
分别计算粘结剂HTPB的校正误差均方根和验证误差均方根；
分别计算增塑剂DOS的校正误差均方根与验证误差均方根；
分别计算氧化剂AP的校正误差均方根与验证误差均方根。
当粘结剂HTPB的相关系数R1大于0.9、校正误差均方根小于0.2、验证误差均方根小于0.2时，该粘结剂HTPB校正模型满足要求。
当增塑剂DOS的相关系数R2大于0.9、校正误差均方根小于0.1、验证误差均方根小于0.1时，该增塑剂DOS校正模型满足要求。
当氧化剂AP的相关系数R3大于0.9、校正误差均方根小于0.3、验证误差均方根小于0.3时，该氧化剂AP校正模型满足要求。
在评价校正模型的性能时，通过公式(2)分别确定粘结剂HTPB1、增塑剂DOS和氧化剂AP的相关系数Rk：
该Rk值表示预测值与实际值关系的线性程度。Rk越接近1，则校正模型的预测值与实际值相关性越强。
通过公式(3)确定校正误差均方根RMSEC：
通过公式(4)确定验证误差均方根RMSEP：
上述各式中：-NIR数学模型预测值；
Ci-组分含量实际值；
-实际值的平均值；
n-建立模型所用的校正集样本数；
m-用于检验模型的验证集样本数：
RMSEC和RMSEP分别代表校正集和验证集的预测值与实际值的偏差。值愈小，则模型的预测精度愈高。如果指标不好，一般需通过不断优化建模参数，以提高校正模型拟合质量。
本实施例中，校正模型性能评价结果见表1。
表1三种组分的模型参数
| 组分 | 相关系数(Rk) | 校正误差均方根(RMSEC) | 验证误差均方根(RMSEP) |
| HTPB | 0.98085 | 0.0906 | 0.0742 |
| KZ | 0.99399 | 0.0388 | 0.0370 |
| AP | 0.95765 | 0.2010 | 0.1650 |
从表中可以看出，粘结剂HTPB和增塑剂DOS的RMSEC、RMSEP均非常小，氧化剂AP稍大；三种组分的相关系数均大于0.95，KZ组分达到0.99399。
步骤6.模型预测。
另采集一张样品，并得到该样品的近红外光谱，通过所建立的校正模型，可同时得到粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量。丁羟推进剂药浆中粘结剂HTPB含量、增塑剂DOS含量和氧化剂AP含量实际值与近红外光谱预测结果分别见附图1、附图2、附图3。图中横坐标为组分含量的真实值，纵坐标为近红外光谱预测值。线性关系好，说明预测值与真实值较接近。可以看出，所建模型线性关系显著，模型校正误差和验证误差均较小，此模型预测精度较高。
实施例2
步骤1.建立样品库。通过收集丁羟推进剂药浆样品，建立样品库。所述样品是采用实验室配制和生产现场取样两种方式获得。样品集有99个样品，其中91个样品用来建立校正集，8个用来建立验证集。
步骤2.采集各样品的近红外漫反射光谱。用傅立叶变换红外光谱仪采集样品近红外漫反射光谱。
所用仪器为美国Thermo Nicolet公司的Nexus 670傅立叶变换近红外光谱仪、漫反射积分球装置、OMNIC信号采集软件及TQ analyst数据处理软件。
采集样品的近红外漫反射光谱时，将第一个样品分为三份，分别放入旋转杯内并平铺为1cm的厚度；用称量瓶轻压药浆以使样品紧贴杯底。分别旋转各装有样品的旋转杯，通过漫反射积分球装置对所述样品进行光谱采集，得到第一个样品的三张近红外漫反射光谱数据，取平均光谱用于数据处理。采集时，扫描次数为64次，分辨率为4cm-1，增益为2，近红外漫反射光谱采集范围为4000～10000cm-1。
重复所述获得第一个样品平均近红外漫反射光谱的过程，分别对其余各样品进行光谱测试；分别得到其余各样品的平均近红外漫反射光谱。共计得到99个样品的平均近红外漫反射光谱。
步骤3.滤波平滑预处理。
以得到的99个样品的平均近红外漫反射光谱分别作为原始光谱。通过Savitzky-Golay滤波对原始光谱进行平滑预处理，以消除各种干扰因素引起的光谱偏差，并对光谱图进行波长选择。波长选择时选择相关性最好的两个光谱区域，分别为6800cm-1～5550cm-1光谱区域和5180cm-1～4000cm-1光谱区域。分别得到经过预处理后的各样品的平均近红外漫反射光谱。
步骤4.建立校正模型。
通过交互验证的方法确定建模所需的最佳变量数为12。
将经过预处理后的各样品的平均近红外漫反射光谱分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量进行关联，用偏最小二乘法建立其中91个样品与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量的校正模型，形成校正集。所述91个样品标记为i，i＝1～n；本实施例中，n＝91。
剩余的8个样品作为验证集，标记为j，j＝1～m；本实施例中，m＝8。
步骤5.校正模型的性能评价。
与实施例1相同方法进行校正模型性能评价，校正模型性能评价结果见表2。
表2三种组分的模型参数
| 组分 | 相关系数(Rk) | 校正误差均方根(RMSEC) | 验证误差均方根(RMSEP) |
| HTPB | 0.98423 | 0.0805 | 0.0632 |
| KZ | 0.99370 | 0.0385 | 0.0301 |
| AP | 0.96425 | 0.1810 | 0.1900 |
步骤6.模型预测。
另采集一张待测样品的NIR光谱，调用所建立的NIR模型，经计算机处理即可同时得到HTPB、KZ、AP三组分含量。丁羟推进剂药浆中HTPB、KZ和AP含量实际值与NIR光谱预测结果分别见附图4、附图5、附图6。图中横坐标为组分含量的真实值，纵坐标为NIR光谱预测值。可以看出，所建模型线性关系显著，模型校正误差和验证误差均较小，此模型预测精度较高。
Claims (6)
1.一种丁羟推进剂药浆组分近红外检测方法，其特征在于，具体过程是：
步骤1，建立样品库；
步骤2，采集各样品的近红外漫反射光谱：
用傅立叶变换红外光谱仪，通过旋转杯分别采集各样品的近红外漫反射光谱；
步骤3.滤波平滑预处理：
以得到的各样品的平均近红外漫反射光谱分别作为原始光谱；具体是，通过滤波对原始光谱进行平滑预处理，以消除各种干扰因素引起的光谱偏差，并对光谱图进行波长选择；波长选择时选择相关性最好的两个光谱区域；分别得到经过预处理后的各样品的平均近红外漫反射光谱；
步骤4.建立校正模型：
所述的校正模型是丁羟推进剂药浆吸光度分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量之间的关系；
在建立校正模型前需确定最佳变量的数量；具有最小预测残差平方和的变量数就是模型的最佳变量数；通过公式(1)确定最佳变量数：
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>E</mi>
<mi>S</mi>
<mi>S</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<msub>
<mi>CV</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：PRESS是预测残差平方和，n是样品数；是交互验证预测值；yi是实际值；
将经过预处理后的各样品的平均近红外漫反射光谱分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量进行关联，用偏最小二乘法建立其中93个样品与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量的校正模型，形成校正集；
步骤5.校正模型的性能评价：
所述的校正模型性能包括：
光谱的吸光度分别与粘结剂HTPB、增塑剂DOS和氧化剂AP之间含量的相关系数Rk，其中R为相关系数，k＝1、2、3，分别表示粘结剂HTPB的相关系数、增塑剂DOS的相关系数和氧化剂AP相关系数；
粘结剂HTPB的校正误差均方根与验证误差均方根；
增塑剂DOS的校正误差均方根与验证误差均方根；
氧化剂AP的校正误差均方根与验证误差均方根；
在对校正模型进行性能评价时：
通过相关系数R分别评价粘结剂HTPB的预测值与实际值关系的线性程度、增塑剂DOS的预测值与实际值关系的线性程度、氧化剂AP的预测值与实际值关系的线性程度；分别得到粘结剂HTPB相关系数R1、增塑剂DOS的相关系数R2和氧化剂AP的相关系数R3；
分别计算粘结剂HTPB的校正误差均方根和验证误差均方根；
分别计算增塑剂DOS的校正误差均方根与验证误差均方根；
分别计算氧化剂AP的校正误差均方根与验证误差均方根；
当粘结剂HTPB的相关系数R1大于0.9、校正误差均方根小于0.2、验证误差均方根小于0.2时，该粘结剂HTPB校正模型满足要求；
当增塑剂DOS的相关系数R2大于0.9、校正误差均方根小于0.1、验证误差均方根小于0.1时，该增塑剂DOS校正模型满足要求；
当氧化剂AP的相关系数R3大于0.9、校正误差均方根小于0.3、验证误差均方根小于0.3时，该氧化剂AP校正模型满足要求；
该Rk值表示预测值与实际值关系的线性程度；Rk越接近1，则校正模型的预测值与实际值相关性越强；
步骤6.模型预测：
所述模型预测的具体过程是另采集一张样品，并得到该样品的近红外光谱；将得到的近红外光谱带入至所建立的丁羟推进剂药浆吸光度分别与粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量之间关系的模型中，同时得到粘结剂HTPB含量、增塑剂DOS含量、氧化剂AP含量。
2.如权利要求1所述丁羟推进剂药浆组分近红外检测方法，其特征在于，所述建立的样品库是丁羟推进剂药浆样品库，具体过程是采用实验室配制和生产现场取样两种方式获得，并将获得的样品分别用于建立校正集和建立验证集。
3.如权利要求1所述丁羟推进剂药浆组分近红外检测方法，其特征在于，
采集样品的近红外漫反射光谱时，将第一个样品分为三份，分别放入旋转杯内并平铺为1cm的厚度；用称量瓶轻压药浆以使样品紧贴杯底；分别旋转各装有样品的旋转杯，通过漫反射积分球装置对所述样品进行光谱采集，得到第一个样品的三张近红外漫反射光谱数据，取平均光谱用于数据处理；采集时，扫描次数为60～120次，分辨率为4cm-1，增益为2，近红外漫反射光谱采集范围为4000～10000cm-1；
重复所述获得第一个样品平均近红外漫反射光谱的过程，分别对其余各样品进行光谱测试；分别得到其余各样品的平均近红外漫反射光谱。
4.如权利要求1所述丁羟推进剂药浆组分近红外检测方法，其特征在于，所述滤波平滑预处理时，相关性最好的两个光谱区域分别为6800cm-1～5550cm-1光谱区域和5180cm-1～4000cm-1光谱区域。
5.如权利要求1所述丁羟推进剂药浆组分近红外检测方法，其特征在于，在评价校正模型的性能时，通过公式(2)分别确定粘结剂HTPB1、增塑剂DOS和氧化剂AP的相关系数Rk：
<mrow>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>C</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
6.如权利要求1所述丁羟推进剂药浆组分近红外检测方法，其特征在于，通过公式(3)确定校正误差均方根RMSEC：
<mrow>
<mi>R</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
<mi>C</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>i</mi>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mi>n</mi>
</mfrac>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
通过公式(4)确定验证误差均方根RMSEP：
<mrow>
<mi>R</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
<mi>P</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mi>j</mi>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>C</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mi>m</mi>
</mfrac>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
上述各式中：-NIR数学模型预测值；
Ci-组分含量实际值；
-实际值的平均值；
n-建立模型所用的校正集样本数；
m-用于检验模型的验证集样本数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711326992.7A CN108020526A (zh) | 2017-12-13 | 2017-12-13 | 一种丁羟推进剂药浆组分近红外检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711326992.7A CN108020526A (zh) | 2017-12-13 | 2017-12-13 | 一种丁羟推进剂药浆组分近红外检测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108020526A true CN108020526A (zh) | 2018-05-11 |
Family
ID=62073392
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711326992.7A Pending CN108020526A (zh) | 2017-12-13 | 2017-12-13 | 一种丁羟推进剂药浆组分近红外检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108020526A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109916848A (zh) * | 2018-12-17 | 2019-06-21 | 西安航天化学动力有限公司 | 三氟化硼三乙醇胺中硼质量分数近红外检测方法 |
| CN110702635A (zh) * | 2019-09-02 | 2020-01-17 | 内蒙合成化工研究所 | 一种近红外光谱在线检测高能粘合剂组分的方法 |
| CN111811932A (zh) * | 2020-07-16 | 2020-10-23 | 西安航天化学动力有限公司 | 复合固体推进剂单向拉伸力学性能近红外光谱分析方法 |
| CN111811991A (zh) * | 2020-07-16 | 2020-10-23 | 西安航天化学动力有限公司 | 非接触式测试复合固体推进剂密度的近红外光谱分析方法 |
| CN111929155A (zh) * | 2020-07-16 | 2020-11-13 | 西安航天化学动力有限公司 | 复合固体推进剂单向拉伸初始模量的近红外光谱分析方法 |
-
2017
- 2017-12-13 CN CN201711326992.7A patent/CN108020526A/zh active Pending
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109916848A (zh) * | 2018-12-17 | 2019-06-21 | 西安航天化学动力有限公司 | 三氟化硼三乙醇胺中硼质量分数近红外检测方法 |
| CN109916848B (zh) * | 2018-12-17 | 2021-09-07 | 西安航天化学动力有限公司 | 三氟化硼三乙醇胺中硼质量分数近红外检测方法 |
| CN110702635A (zh) * | 2019-09-02 | 2020-01-17 | 内蒙合成化工研究所 | 一种近红外光谱在线检测高能粘合剂组分的方法 |
| CN111811932A (zh) * | 2020-07-16 | 2020-10-23 | 西安航天化学动力有限公司 | 复合固体推进剂单向拉伸力学性能近红外光谱分析方法 |
| CN111811991A (zh) * | 2020-07-16 | 2020-10-23 | 西安航天化学动力有限公司 | 非接触式测试复合固体推进剂密度的近红外光谱分析方法 |
| CN111929155A (zh) * | 2020-07-16 | 2020-11-13 | 西安航天化学动力有限公司 | 复合固体推进剂单向拉伸初始模量的近红外光谱分析方法 |
| CN111811932B (zh) * | 2020-07-16 | 2024-03-26 | 西安航天化学动力有限公司 | 复合固体推进剂单向拉伸力学性能近红外光谱分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108020526A (zh) | 一种丁羟推进剂药浆组分近红外检测方法 | |
| CN1982872B (zh) | 快速检测滴丸剂水分含量的近红外漫反射光谱法 | |
| CN103018195B (zh) | 近红外光谱测定pbx炸药中pctfe含量的方法 | |
| CN101216419A (zh) | 一种快速检测黄酒品质指标的方法 | |
| CN1982874B (zh) | 快速检测滴丸剂有效成分含量的近红外漫反射光谱法 | |
| CN107589140B (zh) | 一种能量色散x射线荧光光谱检测复混肥产品中氯、磷、钾含量的方法 | |
| CN104316489B (zh) | 一种近红外光谱检测灵芝提取物掺假的方法 | |
| CN109142263B (zh) | 一种固体推进剂制备过程的组分分散均匀性在线测试方法 | |
| CN103487395A (zh) | 一种重楼药材多指标快速检测方法 | |
| CN105021564A (zh) | 基于近红外光谱分析技术测定烟叶中麦角甾醇含量的方法 | |
| Xu et al. | Rapid and simultaneous determination of the iodine value and saponification number of edible oils by FTIR spectroscopy | |
| CN103411895B (zh) | 珍珠粉掺伪的近红外光谱鉴别方法 | |
| CN108801975A (zh) | 一种微型化近红外光谱仪检测酒糟成分的光谱预处理方法 | |
| CN105372195A (zh) | 一种微量紫外分光光度计质量检测方法和检测试剂盒 | |
| CN104596979A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶纤维素的方法 | |
| CN112179871A (zh) | 一种酱类食品中己内酰胺含量无损检测的方法 | |
| CN104833652A (zh) | 一种混合材料中纤维素硝酸酯含量的快速分析方法 | |
| CN103487398B (zh) | 一种赖氨酸发酵液的分析方法 | |
| CN115144363A (zh) | 一种采用近红外光谱测定硝酸含量的方法 | |
| CN1982873B (zh) | 快速检测滴丸剂外观质量的近红外漫反射光谱法 | |
| CN105092526A (zh) | 一种基于近红外光谱技术的二元掺伪芝麻油含量的快速测定方法 | |
| CN106872398A (zh) | 一种hmx炸药水分含量快速测量方法 | |
| CN104280359A (zh) | 固体推进剂的衬层与界面反应动力学参数测定方法 | |
| CN107727602B (zh) | 中红外光谱结合向量夹角定量分析三氯蔗糖含量的方法 | |
| CN101140225B (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180511 |