CN107976419A - 一种由油品近红外光谱预测其性质的方法 - Google Patents
一种由油品近红外光谱预测其性质的方法 Download PDFInfo
- Publication number
- CN107976419A CN107976419A CN201610919947.1A CN201610919947A CN107976419A CN 107976419 A CN107976419 A CN 107976419A CN 201610919947 A CN201610919947 A CN 201610919947A CN 107976419 A CN107976419 A CN 107976419A
- Authority
- CN
- China
- Prior art keywords
- sample
- spectrum
- property
- formula
- tested
- 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
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000001228 spectrum Methods 0.000 claims abstract description 195
- 238000010561 standard procedure Methods 0.000 claims abstract description 3
- 239000003921 oil Substances 0.000 claims description 57
- 238000002835 absorbance Methods 0.000 claims description 34
- 230000003595 spectral effect Effects 0.000 claims description 30
- 230000000875 corresponding effect Effects 0.000 claims description 28
- 239000010779 crude oil Substances 0.000 claims description 24
- 239000003502 gasoline Substances 0.000 claims description 21
- 239000002283 diesel fuel Substances 0.000 claims description 15
- 230000000996 additive effect Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 claims description 11
- TVMXDCGIABBOFY-UHFFFAOYSA-N octane Chemical compound CCCCCCCC TVMXDCGIABBOFY-UHFFFAOYSA-N 0.000 claims description 10
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical group [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 9
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 7
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 239000010687 lubricating oil Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 239000003209 petroleum derivative Substances 0.000 claims description 3
- 239000005864 Sulphur Substances 0.000 claims description 2
- 150000001336 alkenes Chemical class 0.000 claims description 2
- 150000004945 aromatic hydrocarbons Chemical class 0.000 claims description 2
- 239000010426 asphalt Substances 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
- 229910052760 oxygen Inorganic materials 0.000 claims description 2
- 239000001301 oxygen Substances 0.000 claims description 2
- 230000010076 replication Effects 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 239000001993 wax Substances 0.000 claims description 2
- 238000012067 mathematical method Methods 0.000 claims 1
- 239000002893 slag Substances 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 9
- 239000000047 product Substances 0.000 description 37
- 239000002253 acid Substances 0.000 description 7
- 229910052717 sulfur Inorganic materials 0.000 description 7
- 239000011593 sulfur Substances 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 230000000052 comparative effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008033 biological extinction Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 150000001335 aliphatic alkanes Chemical class 0.000 description 1
- 238000004523 catalytic cracking Methods 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 238000010183 spectrum analysis 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/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
- 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
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
技术领域
本发明为一种由光谱预测油品性质的方法,具体地说,是一种由油品近红外光谱预测其性质的方法。
背景技术
CN102374975A提出了一种新的油品性质预测方法—库光谱拟合方法(LibrarySpectra Fitting Method),这种方法基于油品的近红外光谱库和光谱拟合技术,其基本原理是:光谱相似的样品的性质也相似,通过光谱库中的一个或多个光谱对未知待测样品的光谱进行拟合,然后根据参与拟合光谱的油品的性质计算出待测样品的性质。
上述库光谱拟合方法的化学实质在于未知样品可由一组库样品按一定比例混合而成,因此,未知样品的待测性质可通过库样品的性质按照混合比例计算得出。对于线性加和的性质如组分含量可通过简单的线性组合计算,对于非线性的性质如凝点和粘度等则需要用到调合规则。
上述方法存在的一个显著问题是光谱拟合采用的是光谱吸光度的绝对值,计算量大,也没有针对性地比较待测样品光谱与库样品光谱之间的差异,预测结果的准确性有待提高。
发明内容
本发明的目的是提供一种由油品近红外光谱预测其性质的方法,该法通过测定待测油品的近红外光谱,利用差谱分析,可快速、准确地预测待测油品的性质。
本发明提供的由油品近红外光谱预测其性质的方法,包括如下步骤:
(1)将油品按性质分类,每类收集一组样品,每组样品数量至少为200个,测定每个样品的近红外光谱,再用常规方法按油品类别测定每个样品需要预测的性质数据,
(2)对每个样品的近红外光谱进行二阶微分,将样品的吸光度与其测定的性质数据对应,按油品类别分别建立近红外光谱数据库,
(3)测定待测样品的近红外光谱,并进行二阶微分,取其吸光度为待测样品光谱向量x,
(4)根据待测油品样品的类型,选择相应类别的近红外光谱数据库,
(5)按式①计算待测样品与所选近红外光谱数据库中每个样品吸光度的欧式距离,选出欧式距离最小的1个库光谱p,为与待测样品最相近的库样品,其对应的性质数据为h,
式①中,dj为待测样品光谱向量x与数据库中第j个样品光谱向量之间的距离,xi为待测样品光谱第i个波长点处的吸光度,vj,i为数据库中第j个样品第i个波长点处的吸光度,m为组成样品光谱的吸光度的采样点数,
(6)按式②将待测样品光谱x减去所选数据库中与其最相近样品的光谱p,得到待测样品的差谱△x,
Δx=x-p ②
(7)对所选数据库中每个样品,按式③计算其与光谱p的差谱△vi,
△vi=vi–p i=1,2,3,…,k, ③
式③中,k为所选近红外光谱数据库的样品数,vi为所选数据库中第i个样品的吸光度,i为所选数据库中第i个样品,
(8)对所选数据库中的每个样品,按式④计算其性质与最相近样品的性质h的数据差△qi,
△qi=qi–h i=1,2,3,…,k, ④
式④中,qi为所选数据库中第i个样品的性质数据,
(9)按下述方法用所选库光谱的差谱对待测样品的差谱△x进行拟合,
a)按式⑤表示待测样品的差谱△x,求得参与拟合的库样品的拟合系数:
式中,ai为第i个数据库光谱对应的拟合系数,
拟合系数ai采用非负约束最小二乘法求出,即满足如下目标函数:
b)将上述求得的拟合系数ai中所有的非零拟合系数提取出来,按式⑥进行归一化处理,得到归一化的拟合系数bi,
式⑥中,g为非零拟合系数的库光谱个数,
c)按式⑦计算待测样品的拟合差光谱:
d)按式⑧计算待测样品差光谱的拟合度s:
式⑧中,△xj为待测样品差谱的第j个波长点的吸光度,为拟合差谱第j个波长点的吸光度,m为组成样品差谱的吸光度的采样点数,
(10)若s大于设定的阈值,按式⑩由参与拟合的库样品对应的性质数据计算待测样品的性质,
式⑩中,为预测的待测样品的性质数据,bi为参与拟合的库样品光谱的归一化拟合系数,△qi为参与拟合的库样品光谱对应的性质数据差,h为与待测样品最相近的库样品对应的性质数据。
本发明方法通过已知油品与待测样品光谱的差谱及性质数据之间的差,对待测油品的差谱进行拟合,由参与拟合的已知油品的性质预测待测样品,提高了预测速度和预测准确性。
具体实施方式
本发明方法利用不同油品光谱之间差异(称为差谱)包含的信息,体现油品之间的性质或组成差异,由已知油品光谱的差谱对待测油品进行拟合运算,更加突出了油品光谱之间差异,具有更为准确的预测结果。
本发明方法首先利用模式识别方法从已知油品光谱库中选出与待测油品最相似的1个库光谱,分别将油品库中所有的光谱都减去这个库光谱,得到库光谱差谱,再将待测油品的光谱减去这个库光谱,得到待测油品光谱的差谱,然后用这些差谱对待测油品光谱的差谱进行拟合运算,由参与拟合的样品的性质预测待测油品的性质。
本发明方法(1)步是将油品分类,每类油品均收集一组样品,每类油品的样品数量在200个以上。再按标准方法测其性质数据,以建立红外光谱数据库。所述油品按性质可分成汽油、柴油、润滑油、渣油或原油,因此,建立的红外光谱数据库也分成上述种类。测定的性质数据对汽油样品,预测的性质为汽油密度、辛烷值、芳烃、烯烃、苯和氧含量中的一种或多种;对柴油样品,预测的性质为十六烷值、密度、冷滤点和多环芳烃含量中的一种或多种;对原油样品,预测的性质为原油密度、残炭、硫、蜡、胶质、沥青质含量和实沸点蒸馏收率中的一种或多种。
(2)步是测定收集的油品样品的近红外光谱进行二阶微分,将吸光度与样品对应的性质建立不同油品类别的近红外光谱数据库。选取油品光谱吸光度建库时,对不同的油品选取的光谱范围有所不同。当所述油品为汽油或柴油等轻质油品时,测定近红外光谱选取的光谱范围优选为10000~4000cm-1。当所述油品为润滑油、渣油或原油等重质油品时,测定近红外光谱选取的光谱范围优选为7000~4000cm-1。
(3)、(4)步为测定待测样品的近红外光谱,其测定条件应与(2)步建立近红外光谱数据库测定光谱的条件相同。然后根据待测油品的类型,选择相应的近红外光谱数据库。
(5)步为采用欧式距离方法,从近红外光谱数据库中选出与待测样品距离最小的光谱,为与待测样品最相近的1个库样品的光谱p。欧式距离采用式①计算,其中所述的采样点数m为在光谱区间按波长点采集吸光度的数目。
(6)步为计算待测样品光谱与数据库中与其最相近的光谱p的差谱△x,(7)步为分别计算数据库中每个库样品光谱与待测样品最接近的光谱p的差谱△vi,
(8)步计算数据库中的每个样品的性质与待测样品最接近的样品光谱p对应的性质的数据差值。
(9)步为用库光谱的差谱对待测样品的差谱进行拟合,采用非负约束最小二乘法求出非零的拟合系数,由具有非零的拟合系数的样品差谱对待测样品的差谱进行拟合,算出拟合度进行判断。
本发明方法中,拟合度s是判断库光谱对待测样品拟合完全与否的指标,该值越大说明拟合度越高,由拟合样品计算得到的预测性质越准确。如果s小于设定的阈值,说明拟合不完全,即待测样品不能完全由库样品光谱拟合表示,无法准确对其性质进行预测。
本发明所述拟合度阈值sv通过光谱测量重复性确定,具体方法为:选取一个样品重复测定三次近红外光谱,对每次测定的光谱进行二阶微分处理后取其吸光度,然后计算每次测定的光谱与平均光谱的差谱,按式⑨计算任意两次测定的差谱之间的伪拟合度(sr)值,取最大的sr值,乘以系数0.75即为阈值,
式⑨中,△x’j和△x”j为两次重复测定的光谱与平均光谱之间差谱的第j个波长点的吸光度,m为组成样品差谱的吸光度的采样点数。
若拟合计算得到的s大于等于设定的阈值,则通过式⑩计算待测样品的性质数据。对于满足线性加合性的性质,可直接用式⑩预测其性质。
对于性质数据不服从线性加和性的性质,则需要将其按照一定的规则进行数学变换,变为可线性加和性的数据后,再按式⑩得到预测数据,再将其按转换为线性加合性的逆运算得到呈非线性加合性的性质预测值。如原油的50℃粘度即呈非线性加合性,不能用式⑩直接得到,需将粘度v用式转变成有线性加合性的粘度系数C,
再按式⑩由参与拟合的原油样品的粘度系数C计算得到待测原油的粘度系数再将其代入式得到待测原油的粘度。
下面通过实例进一步详细说明本发明,但本发明并不限于此。
实例1
预测原油的密度、酸值、硫含量、残炭和粘度值。
(1)建立原油样品的近红外光谱数据库
收集有代表性的原油样品450个,原油品种基本覆盖了世界主要原油产区。测定原油样品的近红外光谱,选取7000~4000cm-1光谱范围的吸光度,对其进行二阶微分处理。
分别用SH/T 0604、GB/T7304、GB/T 17040、GB/T17144、GB/T11137测定每个样品的密度、酸值、硫含量、残炭和50℃粘度。对于50℃粘度,由于其不具有线性加和性,需要将50℃粘度值转换成50℃粘度系数值C,C=1000lglg(ν+0.8),将处理后的近红外光谱与其对应的五种性质数据建立原油近红外光谱数据库。
(2)计算拟合度的阈值sv
取1个原油样品,重复测定三次近红外光谱,选取7000~4000cm-1光谱范围的吸光度,进行二阶微分处理后,计算三次测定光谱的平均光谱,然后计算每次测定的光谱与平均光谱的差谱,按式⑨计算任意两次测定差谱之间的伪拟合度(sr)值,取最大的sr值,乘以系数0.75即为sv,本实例计算的sv=2.78。
(3)从原油近红外光谱数据库中选取与待测样品最相近的库样品
取待测样品A,按与(1)步相同的方法测定其近红外光谱,取7000~4000cm-1光谱范围的吸光度,进行二阶微分。按照式①计算其与原油近红外光谱数据库中各样品光谱的欧式距离,欧式距离最小的为第28号样品,为与待测样品最相近的库样品,其光谱用v28表示,性质数据用q28表示。
(4)构建原油近红外光谱的差谱库和性质数据差库
按式③将近红外光谱数据库中每个库样品的光谱减去光谱v28,得到每个样品的差谱;再按式④将每个样品的性质数据减去性质数据q28,得到每个样品的性质数据差。用所有样品的差谱与对应的性质数据差建立原油近红外光谱差谱库和性质数据差库。
(5)对待测样品的光谱差进行拟合
按式②,将待测样品A的光谱x减去光谱v28,得到待测样品A的差谱△x。由建立的原油近红外光谱差谱库和性质数据差库,按照本发明(9)步a)的方法,按式⑤对差谱△x进行拟合,将计算出的不为零的拟合系数及对应库光谱编号列于表1,将这5个库光谱的拟合系数代入式⑥可计算归一化的拟合系数bi,再按式⑦加和得到拟合的待测样品的光谱,用式⑧计算出拟合差谱的拟合度s,s为4.25,大于拟合度的阈值sv,说明待测样品A得到了较好的拟合,可以用拟合系数计算待测样品的性质数据。
表1
(6)对待测样品的性质进行预测分析
由归一化的拟合系数和参与拟合的库光谱样品对应的性质数据,由式⑩计算待测样品的下列性质数据,密度、酸值、硫含量和残炭。其中,为待测样品密度、酸值、硫含量和残炭的预测值,bi为参与拟合的库样品归一化的拟合系数,△qi为参与拟合的库光谱样品对应的密度、酸值、硫含量和残炭值的差值,h为第28号库样品的性质数据q28。
对于50℃粘度,先按式计算参与拟合的库样品的粘度系数,将其代入式⑩,计算出为待测样品的50℃粘度系数再由公式得到最终的50℃粘度预测值。本发明方法对待测原油样品A的密度、酸值、硫含量、残炭和50℃粘度的预测结果见表2。
对比例1
按CN102374975A的方法预测原油样品A的密度、酸值、硫含量、残炭和50℃粘度,结果见表2。
表2
实例2
预测重整汽油辛烷值。
(1)建立重整汽油样品的近红外光谱数据库
收集有代表性的重整汽油样品1687个,重整汽油样品基本覆盖了多种工艺和催化剂的产物,其研究法辛烷值的分布范围为93.2~104.6。测定重整汽油样品的近红外光谱,选取10000~4000cm-1光谱范围的吸光度,对其进行二阶微分处理。用GB/T5487方法测定每个样品的研究法辛烷值(RON)。将处理后得到的重整汽油样品的近红外光谱与其对应的RON建立汽油近红外光谱数据库。
(2)计算拟合度的阈值sv
取1个重整汽油样品,重复测定三次近红外光谱,选取10000~4000cm-1光谱范围的吸光度,进行二阶微分处理后,计算三次测定光谱的平均光谱,然后计算每次测定的光谱与平均光谱的差谱,按式⑨计算任意两次测定差谱之间的伪拟合度(sr)值,取最大的sr值,乘以系数0.75即为sv,本实例计算的sv=4.86。
(3)从汽油近红外光谱性质数据库中选取与待测样品最相近的库样品
取1个待测样品B,按与(1)步相同的方法测定其近红外光谱,取10000~4000cm-1光谱范围的吸光度,进行二阶微分。按照式①计算其与原油近红外光谱数据库中各样品光谱的欧式距离,欧式距离最小的为第596号样品,为与待测样品最相近的库样品,其光谱用v596表示,辛烷值性质数据用q596表示。
(4)构建汽油近红外光谱差谱库和性质数据差库
按式③将近红外光谱数据库中每个库样品的光谱减去光谱v596,得到每个样品的差谱;再按式④将每个样品的性质数据减去性质数据q596,得到每个样品的性质数据差。用所有样品的差谱与对应的性质数据差建立汽油近红外光谱差谱库和性质数据差库。
(5)对待测样品的光谱差进行拟合
按式②,将待测样品B的光谱x减去光谱v596,得到待测样品B的差谱△x。由建立的汽油近红外光谱差谱库和性质数据差库,按照本发明(9)步a)的方法,按式⑤对差谱△x进行拟合,将计算出的不为零的拟合系数及对应库光谱编号列于表3,将这7个库光谱的拟合系数代入式⑥可计算归一化的拟合系数bi,再按式⑦加和得到拟合的待测样品的光谱,用式⑧计算出拟合差谱的拟合度s,s为9.13,大于拟合度的阈值sv,说明待测样品B得到了较好的拟合,可以用拟合系数计算待测样品的性质数据。
表3
(6)对待测样品的性质进行预测分析
由归一化的拟合系数和参与拟合的库光谱样品对应的性质数据,由式⑩计算待测样品的研究法辛烷值。其中,为待测样品辛烷值的预测值,bi为参与拟合的库样品归一化的拟合系数,△qi为参与拟合的库光谱样品对应的辛烷值的差值,h为第596号库样品的辛烷值数据q596。预测结果见表4。
对比例2
按CN102374975A的方法预测原油样品B的研究法辛烷值,结果见表4。
表4
实例3
预测柴油多环芳烃含量。
(1)建立柴油样品的近红外光谱数据库
收集有代表性的柴油样品482个,柴油样品的类型包括直馏柴油、催化裂化柴油、加氢柴油。测定柴油样品的近红外光谱,选取10000~4000cm-1光谱范围的吸光度,对其进行二阶微分处理。用SH/T0606方法测定每个样品的多环芳烃含量。将处理后得到的柴油近红外光谱与其对应的多环芳烃含量建立柴油近红外光谱数据库。
(2)计算拟合度的阈值sv
取1个柴油样品,重复测定三次其近红外光谱,选取10000~4000cm-1光谱范围的吸光度,进行二阶微分处理后,计算三次测定光谱的平均光谱,然后计算每次测定的光谱与平均光谱的差谱,按式⑨计算任意两次测定差谱之间的伪拟合度(sr)值,取最大的sr值,乘以系数0.75即为sv,本实例计算的sv=3.17。
(3)从柴油近红外光谱性质数据库中选取与待测样品最相近的库样品
取1个待测样品C,按与(1)步相同的方法测定其近红外光谱,取10000~4000cm-1光谱范围的吸光度,进行二阶微分。按照式①计算其与柴油近红外光谱数据库中各样品光谱的欧式距离,欧式距离最小的为第315号样品,为与待测样品最相近的库样品,其光谱用v315表示,多环芳烃含量性质数据用q315表示。
(4)构建柴油近红外光谱差谱库和性质数据差库
按式③将近红外光谱数据库中每个库样品的光谱减去光谱v315,得到每个样品的差谱;再按式④将每个样品的性质数据减去性质数据q315,得到每个样品的性质数据差。用所有样品的差谱与对应的性质数据差建立汽油近红外光谱差谱库和性质数据差库。
(5)对待测样品的光谱差进行拟合
按式②,将待测样品C的光谱x减去光谱v315,得到待测样品C的差谱△x。由建立的柴油近红外光谱差谱库和性质数据差库,按照本发明(9)步a)的方法,按式⑤对差谱△x进行拟合,将计算出的不为零的拟合系数及对应库光谱编号列于表5,将这6个库光谱的拟合系数代入式⑥可计算归一化的拟合系数bi,再按式⑦加和得到拟合的待测样品的光谱,用式⑧计算出拟合差谱的拟合度s,s为4.22,大于拟合度的阈值sv,说明待测样品C得到了较好的拟合,可以用拟合系数计算待测样品的性质数据。
表5
(6)对待测样品的性质进行预测分析
由归一化的拟合系数和参与拟合的库光谱样品对应的性质数据,由式⑩计算待测样品的多环芳烃含量。其中,为待测样品多环芳烃含量的预测值,bi为参与拟合的库样品归一化的拟合系数,△qi为参与拟合的库光谱样品对应的多环芳烃含量的差值,h为第315号库样品的多环芳烃含量性质数据q315。预测结果见表6。
对比例3
按CN102374975A的方法预测原油样品C的多环芳烃含量,结果见表6。
表6
由表2、4、6可知,本发明方法较之对比例,预测准确性更高。
Claims (7)
1.一种由油品近红外光谱预测其性质的方法,包括如下步骤:
(1)将油品按性质分类,每类收集一组样品,每组样品数量至少为200个,测定每个样品的近红外光谱,再用标准方法按油品类别测定每个样品需要预测的性质数据,
(2)对每个样品的近红外光谱进行二阶微分,将样品的吸光度与其测定的性质数据对应,按油品类别分别建立近红外光谱数据库,
(3)测定待测样品的近红外光谱,并进行二阶微分,取其吸光度为待测样品光谱向量x,
(4)根据待测油品样品的类型,选择相应类别的近红外光谱数据库,
(5)按式①计算待测样品与所选近红外光谱数据库中每个样品吸光度的欧式距离,选出欧式距离最小的1个库光谱p,为与待测样品最相近的库样品,其对应的性质数据为h,
式①中,dj为待测样品光谱向量x与数据库中第j个样品光谱向量之间的距离,xi为待测样品光谱第i个波长点处的吸光度,vj,i为数据库中第j个样品第i个波长点处的吸光度,m为组成样品光谱的吸光度的采样点数,
(6)按式②将待测样品光谱x减去所选数据库中与其最相近样品的光谱p,得到待测样品的差谱△x,
Δx=x-p ②
(7)对所选数据库中每个样品,按式③计算其与光谱p的差谱△vi,
△vi=vi–p i=1,2,3,…,k, ③
式③中,k为所选近红外光谱数据库的样品数,vi为所选数据库中第i个样品的吸光度,i为所选数据库中第i个样品,
(8)对所选数据库中的每个样品,按式④计算其性质与最相近样品的性质h的数据差△qi,
△qi=qi–h i=1,2,3,…,k, ④
式④中,qi为所选数据库中第i个样品的性质数据,
(9)按下述方法用所选库光谱的差谱对待测样品的差谱△x进行拟合,
a)按式⑤表示待测样品的差谱△x,求得参与拟合的库样品的拟合系数:
式中,ai为第i个数据库光谱对应的拟合系数,
拟合系数ai采用非负约束最小二乘法求出,即满足如下目标函数:
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>|</mo>
<mo>|</mo>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&Delta;v</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>,</mo>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
b)将上述求得的拟合系数ai中所有的非零拟合系数提取出来,按式⑥进行归一化处理,得到归一化的拟合系数bi,
式⑥中,g为非零拟合系数的库光谱个数,
c)按式⑦计算待测样品的拟合差谱:
d)按式⑧计算待测样品差谱的拟合度s:
式⑧中,△xj为待测样品差谱的第j个波长点的吸光度,为拟合差谱第j个波长点的吸光度,m为样品差谱的吸光度的采样点数,
(10)若s大于设定的阈值,按式⑩由参与拟合的库样品对应的性质数据计算待测样品的性质,
式⑩中,为预测的待测样品的性质数据,bi为参与拟合的库样品光谱的归一化拟合系数,△qi为参与拟合的库样品光谱对应的性质数据差,h为与待测样品最相近的库样品对应的性质数据。
2.按照权利要求1所述的方法,其特征在于所述油品按性质分成汽油、柴油、润滑油、渣油或原油。
3.按照权利要求1所述的方法,其特征在于当所述油品为汽油或柴油时,测定近红外光谱选取的光谱范围为10000~4000cm-1。
4.按照权利要求1所述的方法,其特征在于当所述油品为润滑油、渣油或原油时,测定近红外光谱选取的光谱范围为7000~4000cm-1。
5.按照权利要求1所述的方法,其特征在于(10)步所述阈值的确定方法为:选取一个样品重复测定三次近红外光谱,对每次测定的光谱进行二阶微分处理后取其吸光度,然后计算每次测定的光谱与平均光谱的差谱,按式⑨计算任意两次测定的差谱之间的伪拟合度(sr)值,取最大的sr值,乘以系数0.75即为阈值,
式⑨中,△x′j和△x″j为两次重复测定的光谱与平均光谱之间差谱的第j个波长点的吸光度,m为组成样品差谱的吸光度的采样点数。
6.按照权利要求1所述的方法,其特征在于计算待测样品非线性加和性的性质数据的方法是:将参与拟合的光谱对应的性质数据通过数学方法转换成有线性加合性的数据,将其代入式⑩得到预测数据,再将其按转换为线性加合性的逆运算得到呈非线性加合性的性质预测值。
7.按照权利要求1所述的方法,其特征在于所述的油品为汽油时,预测的性质为汽油密度、辛烷值、芳烃、烯烃、苯和氧含量中的一种或多种;所述的油品为柴油时,预测的性质为柴油十六烷值、密度、冷滤点和多环芳烃含量中的一种或多种;所述的油品为原油时,预测的性质为原油密度、残炭、硫、蜡、胶质、沥青质含量和实沸点蒸馏收率中的一种或多种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610919947.1A CN107976419B (zh) | 2016-10-21 | 2016-10-21 | 一种由油品近红外光谱预测其性质的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610919947.1A CN107976419B (zh) | 2016-10-21 | 2016-10-21 | 一种由油品近红外光谱预测其性质的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107976419A true CN107976419A (zh) | 2018-05-01 |
CN107976419B CN107976419B (zh) | 2020-06-16 |
Family
ID=62004526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610919947.1A Active CN107976419B (zh) | 2016-10-21 | 2016-10-21 | 一种由油品近红外光谱预测其性质的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107976419B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109283153A (zh) * | 2018-11-26 | 2019-01-29 | 佛山市海天(高明)调味食品有限公司 | 一种酱油定量分析模型的建立方法 |
CN109406436A (zh) * | 2018-10-30 | 2019-03-01 | 甘肃省交通规划勘察设计院股份有限公司 | 基于红外光谱大数据的道路沥青快速检测预警方法 |
CN110376157A (zh) * | 2019-07-19 | 2019-10-25 | 应急管理部天津消防研究所 | 一种判定泡沫塑料种类及其阻燃性能的检测分析方法 |
CN110763649A (zh) * | 2018-07-26 | 2020-02-07 | 中国石油化工股份有限公司 | 一种由近红外光谱及性质选择目标原油调合配方的方法 |
CN110987858A (zh) * | 2019-11-30 | 2020-04-10 | 江苏万标检测有限公司 | 一种利用神经网络数据模型进行油品快速检测的方法 |
CN111044482A (zh) * | 2018-10-11 | 2020-04-21 | 中国石油化工股份有限公司 | 原油调配方法 |
CN111257277A (zh) * | 2018-11-30 | 2020-06-09 | 湖南中烟工业有限责任公司 | 一种基于近红外光谱技术的烟叶相似度判定方法 |
CN111829976A (zh) * | 2019-04-18 | 2020-10-27 | 中国石油化工股份有限公司 | 一种由原油近红外光谱预测其汽油馏分烃族组成的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102374975A (zh) * | 2010-08-19 | 2012-03-14 | 中国石油化工股份有限公司 | 一种利用近红外光谱预测油品物性数据的方法 |
-
2016
- 2016-10-21 CN CN201610919947.1A patent/CN107976419B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102374975A (zh) * | 2010-08-19 | 2012-03-14 | 中国石油化工股份有限公司 | 一种利用近红外光谱预测油品物性数据的方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110763649A (zh) * | 2018-07-26 | 2020-02-07 | 中国石油化工股份有限公司 | 一种由近红外光谱及性质选择目标原油调合配方的方法 |
CN110763649B (zh) * | 2018-07-26 | 2022-01-04 | 中国石油化工股份有限公司 | 一种由近红外光谱及性质选择目标原油调合配方的方法 |
CN111044482A (zh) * | 2018-10-11 | 2020-04-21 | 中国石油化工股份有限公司 | 原油调配方法 |
CN111044482B (zh) * | 2018-10-11 | 2022-07-12 | 中国石油化工股份有限公司 | 原油调配方法 |
CN109406436A (zh) * | 2018-10-30 | 2019-03-01 | 甘肃省交通规划勘察设计院股份有限公司 | 基于红外光谱大数据的道路沥青快速检测预警方法 |
CN109283153A (zh) * | 2018-11-26 | 2019-01-29 | 佛山市海天(高明)调味食品有限公司 | 一种酱油定量分析模型的建立方法 |
CN109283153B (zh) * | 2018-11-26 | 2021-02-26 | 佛山市海天(高明)调味食品有限公司 | 一种酱油定量分析模型的建立方法 |
CN111257277A (zh) * | 2018-11-30 | 2020-06-09 | 湖南中烟工业有限责任公司 | 一种基于近红外光谱技术的烟叶相似度判定方法 |
CN111257277B (zh) * | 2018-11-30 | 2023-02-17 | 湖南中烟工业有限责任公司 | 一种基于近红外光谱技术的烟叶相似度判定方法 |
CN111829976A (zh) * | 2019-04-18 | 2020-10-27 | 中国石油化工股份有限公司 | 一种由原油近红外光谱预测其汽油馏分烃族组成的方法 |
CN110376157A (zh) * | 2019-07-19 | 2019-10-25 | 应急管理部天津消防研究所 | 一种判定泡沫塑料种类及其阻燃性能的检测分析方法 |
CN110987858A (zh) * | 2019-11-30 | 2020-04-10 | 江苏万标检测有限公司 | 一种利用神经网络数据模型进行油品快速检测的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107976419B (zh) | 2020-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107976419A (zh) | 一种由油品近红外光谱预测其性质的方法 | |
CN105388123B (zh) | 一种由近红外光谱预测原油性质的方法 | |
CN102374975B (zh) | 一种利用近红外光谱预测油品物性数据的方法 | |
US10228708B2 (en) | Method for optimising a mixture of components by means of spectral analysis | |
RU2502984C2 (ru) | Способ прогнозирования устойчивости технологического потока углеводородов с использованием ближних инфракрасных спектров | |
KR100902806B1 (ko) | 정제 공급원료 또는 정제 공정 제품의 휴대용 분석 장치 | |
CN105466884B (zh) | 一种通过近红外光谱识别原油种类及其性质的方法 | |
WO2017007845A1 (en) | Method for correlating physical and chemical measurement data sets to predict physical and chemical properties | |
AU2006206585B2 (en) | Method for modification of a synthetically generated assay using measured whole crude properties | |
Pinto et al. | Use of 1H NMR and chemometrics to detect additives present in the Brazilian commercial gasoline | |
US10241040B2 (en) | Method for characterizing a product by means of topological spectral analysis | |
CN107976416B (zh) | 一种由红外光谱预测原油性质的方法 | |
US20160061720A1 (en) | Method for Characterising a Product by Means of Topological Spectral Analysis | |
JP2012246490A5 (zh) | ||
US11415568B2 (en) | Systems and methods for implicit chemical resolution of vacuum gas oils and fit quality determination | |
CN109724938B (zh) | 由近红外光谱预测润滑油基础油性质的方法 | |
CN107976420B (zh) | 一种由近红外光谱预测混合原油组成的方法 | |
CN110763651B (zh) | 一种适合生产润滑油基础油的原油调配方法 | |
Morales-Medina et al. | Prediction of density and viscosity of colombian crude oils from chromatographic data | |
US9915606B2 (en) | Method for transferring between spectrometers | |
CN109724939A (zh) | 由近红外光谱预测加氢尾油性质的方法 | |
Tertyshna et al. | Forming of crude oil mixtures with increased yield of target fractions | |
WO2024011687A1 (zh) | 一种油品物性快评模型建立方法及装置 | |
CN117476127A (zh) | 利用光谱信息预测油品物性的方法及装置 | |
CN109724937A (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 |