CN117470800A - 一种由油品近红外光谱预测其物性的方法及装置 - Google Patents
一种由油品近红外光谱预测其物性的方法及装置 Download PDFInfo
- Publication number
- CN117470800A CN117470800A CN202210833786.XA CN202210833786A CN117470800A CN 117470800 A CN117470800 A CN 117470800A CN 202210833786 A CN202210833786 A CN 202210833786A CN 117470800 A CN117470800 A CN 117470800A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- sample
- fitting
- near infrared
- measured
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000000704 physical effect Effects 0.000 title claims abstract description 75
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 54
- 238000001228 spectrum Methods 0.000 claims abstract description 131
- 239000013598 vector Substances 0.000 claims abstract description 111
- 239000011159 matrix material Substances 0.000 claims abstract description 68
- 239000002283 diesel fuel Substances 0.000 claims abstract description 32
- 238000002835 absorbance Methods 0.000 claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 230000004069 differentiation Effects 0.000 claims abstract description 8
- 239000003921 oil Substances 0.000 claims description 63
- 238000005070 sampling Methods 0.000 claims description 13
- 230000003595 spectral effect Effects 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 10
- 150000004945 aromatic hydrocarbons Chemical class 0.000 claims description 8
- 238000004497 NIR spectroscopy Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 claims description 6
- 239000003502 gasoline Substances 0.000 claims description 6
- 239000010779 crude oil Substances 0.000 claims description 4
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical group [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 2
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 2
- 150000001336 alkenes Chemical class 0.000 claims description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 2
- 238000009835 boiling Methods 0.000 claims description 2
- 239000000084 colloidal system Substances 0.000 claims description 2
- 238000004821 distillation Methods 0.000 claims description 2
- 239000010687 lubricating oil Substances 0.000 claims description 2
- TVMXDCGIABBOFY-UHFFFAOYSA-N octane Chemical compound CCCCCCCC TVMXDCGIABBOFY-UHFFFAOYSA-N 0.000 claims description 2
- JRZJOMJEPLMPRA-UHFFFAOYSA-N olefin Natural products CCCCCCCC=C JRZJOMJEPLMPRA-UHFFFAOYSA-N 0.000 claims description 2
- 229910052760 oxygen Inorganic materials 0.000 claims description 2
- 239000001301 oxygen Substances 0.000 claims description 2
- 229910052717 sulfur Inorganic materials 0.000 claims description 2
- 239000011593 sulfur Substances 0.000 claims description 2
- 239000001993 wax Substances 0.000 claims description 2
- 239000000203 mixture Substances 0.000 abstract description 25
- 238000010561 standard procedure Methods 0.000 abstract description 5
- 230000000875 corresponding effect Effects 0.000 description 26
- 150000002430 hydrocarbons Chemical class 0.000 description 11
- 239000004215 Carbon black (E152) Substances 0.000 description 10
- UFWIBTONFRDIAS-UHFFFAOYSA-N Naphthalene Chemical compound C1=CC=CC2=CC=CC=C21 UFWIBTONFRDIAS-UHFFFAOYSA-N 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 10
- 229930195733 hydrocarbon Natural products 0.000 description 10
- -1 monocyclic aromatic hydrocarbon Chemical class 0.000 description 9
- PQNFLJBBNBOBRQ-UHFFFAOYSA-N indane Chemical compound C1=CC=C2CCCC2=C1 PQNFLJBBNBOBRQ-UHFFFAOYSA-N 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 150000004996 alkyl benzenes Chemical class 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000000513 principal component analysis Methods 0.000 description 5
- CXWXQJXEFPUFDZ-UHFFFAOYSA-N tetralin Chemical compound C1=CC=C2CCCCC2=C1 CXWXQJXEFPUFDZ-UHFFFAOYSA-N 0.000 description 5
- 150000001335 aliphatic alkanes Chemical class 0.000 description 4
- 150000001924 cycloalkanes Chemical class 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 125000002619 bicyclic group Chemical group 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 150000002468 indanes Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 125000002950 monocyclic group Chemical group 0.000 description 1
- 150000002790 naphthalenes Chemical class 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 238000010986 on-line near-infrared spectroscopy Methods 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 239000012188 paraffin wax Substances 0.000 description 1
- 238000012567 pattern recognition method Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000004808 supercritical fluid chromatography Methods 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
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
本发明实施例提供一种由油品近红外光谱预测其物性的方法及装置，该方法包括：包括收集一组具有代表性的样品，测定每个样品的近红外光谱，再用标准方法测定每个样品的性质数据，将测定样品的近红外光谱经二阶微分后与其测定的性质数据建立数据库。选取6600～8900cm‑1谱区的吸光度建立矩阵X与柴油样品对应的物性或组成数据Y建立柴油样品近红外光谱数据库，将待测柴油样品在上述谱区的吸光度构成光谱矢量x，通过PLS算法求解X＝T×P得到其对应于每一个组成或物性的得分矩阵T，将光谱矩阵X与得分矩阵T进行拼接，用数据库样品光谱的光谱‑得分差谱对待测样品的光谱‑得分差谱进行拟合，对于光谱拟合完全的待测样品，由参与拟合的库光谱的性质数据预测待测柴油的性质。该法可快速、准确地预测待测柴油性质。
Description
技术领域
本发明涉及化工领域，具体地涉及一种由油品近红外光谱预测其物性的方法及装置。
背景技术
柴油馏分馏程约150～350℃，族组成主要分为链烷烃、环烷烃、芳烃三大类组分，关于柴油烃族组成，特别是其中芳烃组成的测定目前主要的分析方法有质谱法和超临界流体色谱法。近红外光谱(NIR)分析技术是目前最有前景且应用最广泛的快速分析方法之一，在线近红外技术也是发展最迅速的过程分析技术之一。NIR主要反映含氢基团X-H(X＝C、N、O)合频和倍频的振动，具有丰富的组成和结构信息，非常适合于分析含大量有机化合物的石油及产品；且由于仪器构造简单，性能稳定，无预处理，分析速度快而被大量用于过程分析。依托标准方法建立数据库，NIR可将分析时间缩短至几分钟，在提高分析效率的同时还能减低测试费用。
CN102374975A提出了一种新的油品性质预测方法—库光谱拟合方法(LibrarySpectra Fitting Method)，这种方法基于油品的近红外光谱库和光谱拟合技术，其基本原理是：光谱相似的样品的性质也相似，通过光谱库中的一个或多个光谱对未知待测样品的光谱进行拟合，然后根据参与拟合光谱的柴油的性质计算出待测样品的性质。上述库光谱拟合方法的化学实质在于未知样品可由一组库样品按一定比例混合而成，因此，未知样品的待测性质可通过库样品的性质按照混合比例计算得出。对于线性加和的性质如组分含量可通过简单的线性组合计算，对于非线性的性质如凝点和粘度等则需要用到调合规则。
上述方法存在的一个显著问题是光谱拟合采用的是光谱吸光度的绝对值，计算耗时较长，且预测结果的准确性也有待提高。
发明内容
本发明实施例的目的是提供一种由油品近红外光谱预测其物性的方法及装置，其通过测定待测油品的近红外光谱，利用差谱分析，可快速、准确地预测待测油品的性质。
为了实现上述目的，本发明实施例提供一种由油品近红外光谱预测其物性的方法，该方法包括：测定待测样品的近红外光谱并进行二阶微分处理，选取特征谱区的吸光度构成光谱矢量；针对待测物性，通过PLS算法求解X＝T×P得到所述光谱矢量的得分向量；将所述光谱矢量与所述得分向量进行拼接，以得到待测样品的光谱-得分组合向量k；计算所述待测样品与油品近红外光谱-得分数据库内的每一样品的欧式距离，并选取欧氏距离最小的光谱-得分组合向量，作为与所述待测样品最相近的库样本的光谱-得分组合向量p，该最相近的库样本的物性为h；计算待测样品的光谱-得分组合向量k与所述光谱-得分组合向量p的差谱△k；计算所述数据库中每个库样品的光谱-得分组合向量与所述光谱-得分组合向量p的差谱△vi；利用所述差谱△vi对待测样品的差谱△k进行拟合，得到拟合系数；对于所述数据库的每一样品，计算其物性与所述最相近的库样本的物性h之间的数据差△qi；根据参与拟合的库样本对应的所述数据差△qi及所述拟合系数，对所述待测样品的所述待测物性进行预测，
其中所述油品近红外光谱-得分数据库通过以下操作而被获得：收集一组具有代表性的油品样品，测定每个样品的近红外光谱，并测定每个样品所需要预测的物性，形成矩阵Y；对每个样品的近红外光谱进行二阶微分，将样品的吸光度与其测定的物性对应，选取特征谱区建立近红外光谱数据库；以及利用PLS方法对所述近红外光谱数据库内的光谱矩阵X分别对应Y矩阵的每个物性进行分解，得到得分矩阵T，并按样品顺序将光谱矩阵X与得分矩阵T进行拼接，得到光谱-得分组合矩阵K。
其中，所述欧氏距离通过以下公式而被计算：
式①中，dj为待测样品的光谱-得分组合向量k与所述油品近红外光谱-得分数据库中第j个样品的光谱-得分组合向量之间的距离，xi为待测样品光谱-得分组合向量第i个点处的数值，vj,i为所述油品近红外光谱-得分数据库中第j个样品第i个点处的数值，m为组成样品光谱-得分组合向量的采样点数。
其中，利用所述差谱△vi对待测样品的差谱△k进行拟合，得到拟合系数，包括：
按式⑤表示待测样品的差谱△k，求得参与拟合的库样品的拟合系数：
式中，ai为第i个库样本的光谱-得分组合向量对应的拟合系数，n为参与拟合的库样本的光谱-得分组合向量的个数，
拟合系数ai采用非负约束最小二乘法求出，即满足如下目标函数：
将上述求得的拟合系数ai中所有的非零拟合系数提取出来，按式⑥进行归一化处理，得到归一化的拟合系数bi，
式⑥中，g为非零拟合系数所对应的库样本的光谱-得分组合向量的个数。
其中，该方法还包括：
按式⑦计算待测样品的拟合差谱：
按式⑧计算待测样品差谱的拟合度s：
式⑧中，△kj为待测样品差谱的第j个采样点的数值，为拟合差谱第j个采样点的数值，m为样品差谱的采样点数，
对所述待测样品的所述待测物性进行预测在所述拟合度大于设定的阈值的情况下被执行。
其中，对所述待测样品的所述待测物性进行预测包括：
按式⑩由参与拟合的库样品对应的物性，计算待测样品的物性，
式⑩中，为预测的待测样品的物性，bi为参与拟合的库样品的光谱-得分组合向量的归一化拟合系数，△qi为参与拟合的库样品光谱-得分组合向量对应的物性数据差，h为与待测样品最相近的库样品对应的物性。
其中，所述阈值的确定方法为：选取一个样品重复测定三次近红外光谱，对每次测定的光谱进行二阶微分处理后取其吸光度，然后计算每次测定的光谱与平均光谱的差谱，按式⑨计算任意两次测定的差谱之间的伪拟合度(sr)值，取最大的sr值，乘以系数0.75即为阈值，
式⑨中，△x′j和△x″j为两次重复测定的光谱与平均光谱之间差谱的第j个波长点的吸光度，m为组成样品差谱的吸光度的样点数。
其中，所述油品包括柴油、汽油、润滑油、以及渣油。
其中，所述特征谱区为6600～8900cm-1。
其中，测定近红外光谱的光谱范围为10000～4000cm-1。
其中，所述的油品为汽油时，预测的物性为汽油密度、辛烷值、芳烃、烯烃、苯和氧含量中的一种或多种；所述油品为柴油时，预测的物性为柴油十六烷值、密度、冷滤点和多环芳烃含量中的一种或多种；所述油品为原油时，预测的物性为原油密度、残炭、硫、蜡、胶质、沥青质含量和实沸点蒸馏收率中的一种或多种。
另一方面，本发明提供一种由油品近红外光谱预测其物性的装置，该装置包括：存储器；以及处理器，该处理器被配置为利用上述由油品近红外光谱预测其物性的方法。
另一方面，本发明提供一种机器可读存储介质，该机器可读存储介质上存储有指令，其特征在于，该指令在被处理器执行时使得所述处理器被配置成执行上述由油品近红外光谱预测其物性的方法。
主成分分析(PCA)算法是化学计量学领域里的一个经典算法，基本模型是X＝T×P，将光谱矩阵X分解为T矩阵(术语：得分矩阵)和P矩阵(术语：载荷矩阵)的乘积。一般来说，主成分分析(PCA)只对X矩阵(光谱矩阵)进行分解，分解为T矩阵(得分矩阵)与P矩阵(载荷矩阵)的乘积，得分矩阵的特点是前几个变量基本涵盖了原矩阵的大部分信息，则所需计算的数据量被大大压缩，因此本案中采用得分替代原光谱矩阵进行计算。但这里的得分与样品的物性、组成(y矩阵)并不能保证有最大的相关性，因此本案采用PLS方法提取主成分，因为PLS分解X矩阵的过程中要考虑与样品的物性、组成保持最大的相关性。那么得到的得分变量就与样品的物性、组成有较大相关性。那么得到的得分变量就与样品的物性、组成有较大相关性，最后得到的拟合结果即对未知样品物性、组成的预测也会更准确。
本发明方法通过已知油品与待测样品光谱-得分组合向量的差谱及物性之间的差，对待测油品的差谱进行拟合，由参与拟合的已知油品的物性预测待测样品，更加突出了油品光谱与物性之间的关系，具有更为准确的预测结果。另外，主成分分析的特点就是可以对数据进行极大的压缩，因此本案还提高了预测速度。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
图1为本发明提供的由油品近红外光谱预测其物性的方法的流程图。
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。
图1是本发明提供的由油品近红外光谱预测其物性的方法的流程图。如图1所示，本发明实施例提供一种由油品近红外光谱预测其物性的方法，包括如下步骤：
S110：测定待测样品的近红外光谱并进行二阶微分处理，选取特征谱区的吸光度构成光谱矢量；
S120：针对待测物性，通过PLS算法求解X＝T×P得到所述光谱矢量的得分向量；
这里使用PLS方法计算X＝T×P，将光谱矩阵X分解为T矩阵(术语：得分矩阵)和P矩阵(术语：载荷矩阵)，分解过程考虑Y是为了保证与Y(物性或组成)有较大相关，这里的P矩阵可以理解为新坐标轴，将待测样品的光谱矢量x直接向得到新坐标轴(P)投影即得到该待测样品的得分向量。
S130：将所述光谱矢量与所述得分向量进行拼接，以得到待测样品的光谱-得分组合向量k；
S140：计算所述待测样品与油品近红外光谱-得分数据库内的每一样品的欧式距离，并选取欧氏距离最小的光谱-得分组合向量，作为与所述待测样品最相近的库样本的光谱-得分组合向量p，该最相近的库样本的物性为h；
S150：计算待测样品的光谱-得分组合向量k与所述光谱-得分组合向量p的差谱△k；
S160：计算所述数据库中每个库样品的光谱-得分组合向量与所述光谱-得分组合向量p的差谱△vi；
S170：利用所述差谱△vi对待测样品的差谱△k进行拟合，得到拟合系数；
S180：对于所述数据库的每一样品，计算其物性与所述最相近的库样本的物性h之间的数据差△qi；
S190：根据参与拟合的库样本对应的所述数据差△qi及所述拟合系数，对所述待测样品的所述待测物性进行预测，
其中所述油品近红外光谱-得分数据库通过以下操作而被获得：
收集一组具有代表性的油品样品，测定每个样品的近红外光谱，并测定每个样品所需要预测的物性，形成矩阵Y；
对每个样品的近红外光谱进行二阶微分，将样品的吸光度与其测定的物性对应，选取特征谱区建立近红外光谱数据库；以及
利用PLS方法对所述近红外光谱数据库内的光谱矩阵X分别对应Y矩阵的每个物性进行分解，得到得分矩阵T，并按样品顺序将光谱矩阵X与得分矩阵T进行拼接，得到光谱-得分组合矩阵K。
需要说明的是，上述油品近红外光谱-得分数据库可在上述方法被执行之前预先被建立。
下面结合具体的实施情况来对本发明的技术方案进行描述，当然，本发明并不限于这些描述，本领域技术人员在阅读本案申请文件之后可预见的其他实施形式也是可行的。另外，在本申请文件中，“物性”、“性质”、“性质数据”与“组成”可互换使用，“样本”与“样品”可互换使用。
本发明方法利用不同油品光谱之间差异(称为差谱)包含的信息，体现油品之间的性质或组成差异，油品近红外光谱经偏最小二乘(PLS)计算后的得分信息，由已知油品光谱的光谱-得分向量的差谱对待测油品的光谱-得分向量的差谱进行拟合运算，体现了油品光谱与物性之间的关系，也更加突出了油品光谱之间差异，具有更为准确的预测结果。
本发明方法首先利用模式识别方法从已知油品光谱库中选出与待测油品最相似的1个库光谱-得分组合向量，分别将油品库中所有的光谱-得分组合向量都减去这个库光谱-得分组合向量，得到库光谱-得分组合向量差谱，再将待测油品的光谱-得分组合向量减去这个库光谱-得分组合向量，得到待测油品光谱-得分组合向量的差谱，然后用这些差谱对待测油品光谱-得分组合向量的差谱进行拟合运算，由参与拟合的样品的性质预测待测油品的性质。
本发明方法包括以下步骤：
(1)收集一组具有代表性的油品样品，样品数量至少为400个，测定每个样品的近红外光谱，再用标准方法测定每个样品油品需要预测的性质数据，如详细烃组成、十六烷值、密度、冷滤点和烃族组成等，形成矩阵Y，
(2)对每个样品的近红外光谱进行二阶微分，将样品的吸光度与其测定的性质数据对应，选取6600～8900cm-1谱区建立近红外光谱库X，
(3)测定待测样品的近红外光谱，并进行二阶微分，取其6600～8900cm-1谱区为待测样品光谱向量，
(4)利用PLS方法对油品近红外库光谱矩阵分别对应Y矩阵每个组成或物性进行分解，分别获取前30-100个主成分对应的得分向量，将得分矩阵T按样本顺序将X与T矩阵进行拼接，得到光谱-得分组合矩阵K，矩阵K的个数与组成或物性个数相同；以及同样方法得到待测样本的光谱-得分组合向量k，
(5)按式①计算待测样品与所选近红外光谱-得分数据库中每个样品的欧式距离，选出欧式距离最小的1个库样本p，为与待测样品最相近的库样品，其对应的性质数据为h，
式①中，dj为待测样品向量k与数据库中第j个样品光谱-得分组合向量之间的距离，xi为待测样品光谱-得分组合向量第i个点处的数值，vj,i为数据库中第j个样品第i个点处的数值，m为组成样品光谱-得分组合向量的采样点数，
(6)按式②将待测样品光谱-得分组合向量k减去所选数据库中与其最相近样品的光谱-得分组合向量p，得到待测样品的差谱△k，
Δk＝k-p ②
(7)对所选数据库中每个样品，按式③计算其与光谱-得分组合向量p的差谱△vi，
△vi＝vi–p i＝1，2，3，…，n， ③
式③中，n为所选近红外光谱-得分组合向量数据库的样品数，vi为所选数据库中第i个样品的数值，i为所选数据库中第i个样品，
(8)对所选数据库中的每个样品，按式④计算其性质与最相近样品的性质h的数据差△qi，
△qi＝qi–h i＝1，2，3，…，n， ④
式④中，qi为所选数据库中第i个样品的性质数据，
(9)按下述方法用所选库光谱-得分组合向量的差谱对待测样品的差谱△k进行拟合，
a)按式⑤表示待测样品的差谱△k，求得参与拟合的库样品的拟合系数：
式中，ai为第i个数据库光谱-得分组合向量对应的拟合系数，
拟合系数ai采用非负约束最小二乘法求出，即满足如下目标函数：
非负约束最小二乘法的具体算法参见文献：C.L.Lawson and R.J.Hanson,Solving Least Squares Problems,Prentice-Hall,Englewood Cliffs,NJ(1974)；160～165
b)将上述求得的拟合系数ai中所有的非零拟合系数提取出来，按式⑥进行归一化处理，得到归一化的拟合系数bi，
式⑥中，g为非零拟合系数的库光谱-得分组合向量个数，
c)按式⑦计算待测样品的拟合差谱：
d)按式⑧计算待测样品差谱的拟合度s：
式⑧中，△kj为待测样品差谱的第j个采样点的数值，为拟合差谱第j个采样点的数值，m为样品差谱的采样点数，
本发明方法中，拟合度s是判断库光谱-得分组合向量对待测样品拟合完全与否的指标，该值越大说明拟合度越高，由拟合样品计算得到的预测性质越准确。如果s小于设定的阈值，说明拟合不完全，即待测样品不能完全由库样品光谱-得分组合向量拟合表示，无法准确对其性质进行预测。
(10)若s大于设定的阈值，按式⑩由参与拟合的库样品对应的性质数据计算待测样品的性质，
式⑩中，为预测的待测样品的性质数据，bi为参与拟合的库样品光谱-得分组合向量的归一化拟合系数，△qi为参与拟合的库样品光谱-得分组合向量对应的性质数据差，h为与待测样品最相近的库样品对应的性质数据。
本发明所述拟合度阈值sv通过光谱-得分组合向量测量重复性确定，具体方法为：选取一个样品重复测定三次近红外光谱-得分组合向量，对每次测定的光谱-得分组合向量进行二阶微分处理后取其吸光度，然后计算每次测定的光谱-得分组合向量与平均光谱-得分组合向量的差谱，按式⑨计算任意两次测定的差谱之间的伪拟合度(sr)值，取最大的sr值，乘以系数0.75即为阈值，
式⑨中，△x′j和△x″j为两次重复测定的光谱-得分组合向量与平均光谱-得分组合向量之间差谱的第j个波长点的吸光度，m为组成样品差谱的吸光度的采样点数。
若拟合计算得到的s大于等于设定的阈值，则通过式⑩计算待测样品的性质数据。对于满足线性加合性的性质，可直接用式⑩预测其性质。
下面通过实例进一步详细说明本发明，但本发明并不限于此。
实例1
预测柴油烃族组成。
(1)建立柴油样品的近红外光谱数据库
收集有代表性的柴油样品450个，测定柴油样品的近红外光谱，选取6600～8900cm-1光谱范围的吸光度，对其进行二阶微分处理，形成柴油近红外光谱矩阵X600×450。
用标准方法SH/T 0606测定每个样品的烃族组成。将处理后的近红外光谱与其对应的性质数据建立柴油近红外光谱数据库。
(2)将柴油近红外光谱矩阵X600×450分别与烃组成矩阵Y8×450的每一个组成(链烷烃、环烷烃、烷基苯、茚满或四氢萘、萘、总单环芳烃、总双环芳烃、总三环芳烃)通过偏最小二乘算法进行主成分分析，各取库中每个柴油光谱分别对应链烷烃、环烷烃、烷基苯、茚满或四氢萘、萘、总单环芳烃、总双环芳烃、总三环芳烃的前40个主成分得分与其对应样本的二阶微分光谱拼接构建为柴油拟合数据库矩阵K1 640×450、K2 640×450、K3 640×450、K4 640×450、K5 640×450、K6 640×450、K7 640×450、K8 640×450。
(3)从柴油近红外光谱-得分组合向量数据库中选取与待测样品最相近的库样品
取待测样品A，按与(1)步相同的方法测定其近红外光谱，取6600～8900cm-1光谱范围的吸光度，进行二阶微分，并采用类似于上述的方法得到待测样品A的光谱-得分组合向量x。按照式①计算其与柴油近红外光谱-得分组合向量数据库中各样品光谱-得分组合向量的欧式距离，欧式距离最小的样本为与待测样品最相近的库样品，其光谱-得分组合向量用v表示，性质数据用q表示。
(4)构建柴油近红外光谱-得分组合向量的差谱库和性质数据差库
按式③将近红外光谱-得分组合向量数据库中每个库样品的光谱-得分组合向量减去光谱-得分组合向量v，得到每个样品的差谱；再按式④将每个样品的性质数据减去性质数据q，得到每个样品的性质数据差。用所有样品的差谱与对应的性质数据差建立柴油近红外光谱-得分组合向量差谱库和性质数据差库。对柴油拟合数据库矩阵K1 640×450、K2 640×450、K3 640×450、K4 640×450、K5 640×450、K6 640×450、K7 640×450、K8 640×450都进行相同的计算。
(5)对待测样品的光谱-得分组合向量差进行拟合
按式②，将待测样品A的光谱-得分组合向量x减去光谱-得分组合向量v，得到待测样品A的差谱△x。由建立的柴油近红外光谱-得分组合向量差谱库和性质数据差库，按照本发明(9)步a)的方法，按式⑤对差谱△x进行拟合，将计算出的不为零的拟合系数，可以用拟合系数计算待测样品的性质数据，对柴油拟合数据库矩阵K1 640×450、K2 640×450、K3 640×450、K4 640×450、K5 640×450、K6 640×450、K7 640×450、K8 640×450都进行相同的计算，得到每个性质或组成的拟合系数。其中对于单环芳烃的拟合系数对应库光谱-得分组合向量编号列于表1，
表1
(6)对待测样品的性质进行预测分析
由对应每一个组成的归一化的拟合系数和参与拟合的库光谱-得分组合向量样品对应的性质数据，由式⑩计算待测样品的下列组成数据：链烷烃、环烷烃、烷基苯、茚满或四氢萘、萘、总单环芳烃、总双环芳烃、总三环芳烃。
表2
对比例1
使用传统偏最小二乘(PLS)方法预测柴油样品A的烃组成，结果见表2。
(1)建立柴油样品的近红外光谱数据库
收集有代表性的柴油样品450个，测定柴油样品的近红外光谱，选取6600～8900cm-1光谱范围的吸光度，对其进行二阶微分处理，形成柴油近红外光谱矩阵X600×450。用标准方法SH/T 0606测定每个样品的烃族组成。
对于柴油光谱矩阵X600×450，其中450为柴油样品个数，600为近红外光谱的数据点数。与其对应的烃组成建立柴油性质矩阵Y8×450，其中450为柴油样品个数，8为柴油烃组成个数(链烷烃、环烷烃、烷基苯、茚满或四氢萘、萘、总单环芳烃、总双环芳烃、总三环芳烃)。柴油近红外光谱矩阵X600×450与柴油性质矩阵Y8×450构成柴油样品的近红外光谱数据库。
(2)将矩阵X与矩阵Y用偏最小二乘法(PLS)分别建立烃组成(链烷烃、环烷烃、烷基苯、茚满或四氢萘、萘、总单环芳烃、总双环芳烃、总三环芳烃)校正模型，对未知柴油样本A进行预测，结果见表2。
结果表明本专利方法测试结果优于传统PLS方法。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (12)
1.一种由油品近红外光谱预测其物性的方法，包括如下步骤：
测定待测样品的近红外光谱并进行二阶微分处理，选取特征谱区的吸光度构成光谱矢量；
针对待测物性，通过PLS算法求解X＝T×P得到所述光谱矢量的得分向量；
将所述光谱矢量与所述得分向量进行拼接，以得到待测样品的光谱-得分组合向量k；
计算所述待测样品与油品近红外光谱-得分数据库内的每一样品的欧式距离，并选取欧氏距离最小的光谱-得分组合向量，作为与所述待测样品最相近的库样本的光谱-得分组合向量p，该最相近的库样本的物性为h；
计算待测样品的光谱-得分组合向量k与所述光谱-得分组合向量p的差谱△k；
计算所述数据库中每个库样品的光谱-得分组合向量与所述光谱-得分组合向量p的差谱△vi；
利用所述差谱△vi对待测样品的差谱△k进行拟合，得到拟合系数；
对于所述数据库的每一样品，计算其物性与所述最相近的库样本的物性h之间的数据差△qi；
根据参与拟合的库样本对应的所述数据差△qi及所述拟合系数，对所述待测样品的所述待测物性进行预测，
其中所述油品近红外光谱-得分数据库通过以下操作而被获得：
收集一组具有代表性的油品样品，测定每个样品的近红外光谱，并测定每个样品所需要预测的物性，形成矩阵Y；
对每个样品的近红外光谱进行二阶微分，将样品的吸光度与其测定的物性对应，选取特征谱区建立近红外光谱数据库；以及
利用PLS方法对所述近红外光谱数据库内的光谱矩阵X分别对应Y矩阵的每个物性进行分解，得到得分矩阵T，并按样品顺序将光谱矩阵X与得分矩阵T进行拼接，得到光谱-得分组合矩阵K。
2.根据权利要求1所述的方法，其中，所述欧氏距离通过以下公式而被计算：
式①中，dj为待测样品的光谱-得分组合向量k与所述油品近红外光谱-得分数据库中第j个样品的光谱-得分组合向量之间的距离，xi为待测样品光谱-得分组合向量第i个点处的数值，vj,i为所述油品近红外光谱-得分数据库中第j个样品第i个点处的数值，m为组成样品光谱-得分组合向量的采样点数。
3.根据权利要求1所述的方法，其中，利用所述差谱△vi对待测样品的差谱△k进行拟合，得到拟合系数，包括：
按式⑤表示待测样品的差谱△k，求得参与拟合的库样品的拟合系数：
式中，ai为第i个库样本的光谱-得分组合向量对应的拟合系数，n为参与拟合的库样本的光谱-得分组合向量的个数，
拟合系数ai采用非负约束最小二乘法求出，即满足如下目标函数：
将上述求得的拟合系数ai中所有的非零拟合系数提取出来，按式⑥进行归一化处理，得到归一化的拟合系数bi，
式⑥中，g为非零拟合系数所对应的库样本的光谱-得分组合向量的个数。
4.根据权利要求3所述的方法，其中，该方法还包括：
按式⑦计算待测样品的拟合差谱：
按式⑧计算待测样品差谱的拟合度s：
式⑧中，△kj为待测样品差谱的第j个采样点的数值，为拟合差谱第j个采样点的数值，m为样品差谱的采样点数，
对所述待测样品的所述待测物性进行预测在所述拟合度大于设定的阈值的情况下被执行。
5.根据权利要求4所述的方法，其中，对所述待测样品的所述待测物性进行预测包括：
按式⑩由参与拟合的库样品对应的物性，计算待测样品的物性，
式⑩中，为预测的待测样品的物性，bi为参与拟合的库样品的光谱-得分组合向量的归一化拟合系数，△qi为参与拟合的库样品光谱-得分组合向量对应的物性数据差，h为与待测样品最相近的库样品对应的物性。
6.根据权利要求4所述的方法，其中所述阈值的确定方法为：选取一个样品重复测定三次近红外光谱，对每次测定的光谱进行二阶微分处理后取其吸光度，然后计算每次测定的光谱与平均光谱的差谱，按式⑨计算任意两次测定的差谱之间的伪拟合度(sr)值，取最大的sr值，乘以系数0.75即为阈值，
式⑨中，△x′j和△x″j为两次重复测定的光谱与平均光谱之间差谱的第j个波长点的吸光度，m为组成样品差谱的吸光度的样点数。
7.根据权利要求1所述的方法，其中所述油品包括柴油、汽油、润滑油、以及渣油。
8.根据权利要求1所述的方法，其中，所述特征谱区为6600～8900cm-1。
9.根据权利要求1所述的方法，其中，测定近红外光谱的光谱范围为10000～4000cm-1。
10.根据权利要求1所述的方法，其中，
所述的油品为汽油时，预测的物性为汽油密度、辛烷值、芳烃、烯烃、苯和氧含量中的一种或多种；
所述油品为柴油时，预测的物性为柴油十六烷值、密度、冷滤点和多环芳烃含量中的一种或多种；
所述油品为原油时，预测的物性为原油密度、残炭、硫、蜡、胶质、沥青质含量和实沸点蒸馏收率中的一种或多种。
11.一种由油品近红外光谱预测其物性的装置，所述装置包括：
存储器；以及
处理器，该处理器被配置为利用根据权利要求1-10中任一项权利要求所述的由油品近红外光谱预测其物性的方法。
12.一种机器可读存储介质，该机器可读存储介质上存储有指令，其特征在于，该指令在被处理器执行时使得所述处理器被配置成执行根据权利要求1-10中任一项权利要求所述的由油品近红外光谱预测其物性的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210833786.XA CN117470800A (zh) | 2022-07-14 | 2022-07-14 | 一种由油品近红外光谱预测其物性的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210833786.XA CN117470800A (zh) | 2022-07-14 | 2022-07-14 | 一种由油品近红外光谱预测其物性的方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117470800A true CN117470800A (zh) | 2024-01-30 |
Family
ID=89635218
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210833786.XA Pending CN117470800A (zh) | 2022-07-14 | 2022-07-14 | 一种由油品近红外光谱预测其物性的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117470800A (zh) |
-
2022
- 2022-07-14 CN CN202210833786.XA patent/CN117470800A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6897071B2 (en) | Topological near infrared analysis modeling of petroleum refinery products | |
| CN107976419B (zh) | 一种由油品近红外光谱预测其性质的方法 | |
| WO2003048759A1 (en) | Method for analyzing an unknown material as a blend of known materials calculated so as to match certain analytical data and predicting properties of the unknown based on the calculated blend | |
| CN105466884B (zh) | 一种通过近红外光谱识别原油种类及其性质的方法 | |
| WO2006078688A2 (en) | Method for modification of a synthetically generated assay using measured whole crude properties | |
| Khanmohammadi et al. | Feature selection strategies for quality screening of diesel samples by infrared spectrometry and linear discriminant analysis | |
| WO2006126978A1 (en) | Topological near infrared analysis modeling of petroleum refinery products | |
| CN107966499B (zh) | 一种由近红外光谱预测原油碳数分布的方法 | |
| Flumignan et al. | Multivariate calibrations in gas chromatographic profiles for prediction of several physicochemical parameters of Brazilian commercial gasoline | |
| JP2012246490A (ja) | 潤滑剤の製造を制御する方法 | |
| US11415568B2 (en) | Systems and methods for implicit chemical resolution of vacuum gas oils and fit quality determination | |
| JP2012246490A5 (zh) | ||
| CN117470800A (zh) | 一种由油品近红外光谱预测其物性的方法及装置 | |
| US6295485B1 (en) | Control of lubricant production by a method to predict a base stock's ultimate lubricant performance | |
| Swarin et al. | Prediction of gasoline properties with near-infrared spectroscopy and chemometrics | |
| Flumignan et al. | Multivariate calibrations on 1H NMR profiles for prediction of physicochemical parameters of Brazilian commercial gasoline | |
| Issa | Parametric and Nonparametric Approaches of Reid Vapor Pressure Prediction for Gasoline Containing Oxygenates: A Comparative Analysis Using Partial Least Squares, Nonlinear, and LOWESS Regression Modelling Strategies with Physical Properties | |
| CN109724938B (zh) | 由近红外光谱预测润滑油基础油性质的方法 | |
| CN109580918B (zh) | 一种用于预测石脑油的分子组成的方法 | |
| CN117476127A (zh) | 利用光谱信息预测油品物性的方法及装置 | |
| US11402323B2 (en) | Systems and processes for performance property determination using optical spectral data | |
| CN116008216A (zh) | 用于检测油液掺混的方法及设备 | |
| CN113324945A (zh) | 近红外光谱预测工业白油中芳烃含量的方法 | |
| Fodor | Analysis of petroleum fuels by midband infrared spectroscopy | |
| CN111829976A (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 |