CN107860743A - 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 - Google Patents
利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 Download PDFInfo
- Publication number
- CN107860743A CN107860743A CN201710985799.8A CN201710985799A CN107860743A CN 107860743 A CN107860743 A CN 107860743A CN 201710985799 A CN201710985799 A CN 201710985799A CN 107860743 A CN107860743 A CN 107860743A
- Authority
- CN
- China
- Prior art keywords
- sample
- crude oil
- near infrared
- property
- spectrum
- 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
- 239000000523 sample Substances 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000010779 crude oil Substances 0.000 claims abstract description 114
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 50
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 239000003921 oil Substances 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 23
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 238000004611 spectroscopical analysis Methods 0.000 claims abstract description 4
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 23
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 22
- 229910052717 sulfur Inorganic materials 0.000 claims description 20
- 239000011593 sulfur Substances 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 17
- 239000002253 acid Substances 0.000 claims description 14
- 238000012549 training Methods 0.000 claims description 14
- 229910052594 sapphire Inorganic materials 0.000 claims description 12
- 239000010980 sapphire Substances 0.000 claims description 12
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical group [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 11
- 229910052757 nitrogen Inorganic materials 0.000 claims description 11
- 239000013307 optical fiber Substances 0.000 claims description 11
- 230000003287 optical effect Effects 0.000 claims description 10
- 239000010426 asphalt Substances 0.000 claims description 9
- 239000002184 metal Substances 0.000 claims description 9
- 229910052751 metal Inorganic materials 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 238000002835 absorbance Methods 0.000 claims description 7
- 238000009835 boiling Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000004821 distillation Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 229910052799 carbon Inorganic materials 0.000 claims 1
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 230000005856 abnormality Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 239000005864 Sulphur Substances 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000010561 standard procedure Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- UFWIBTONFRDIAS-UHFFFAOYSA-N Naphthalene Chemical compound C1=CC=CC2=CC=CC=C21 UFWIBTONFRDIAS-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012628 principal component regression Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000005292 vacuum distillation 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
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
技术领域
本发明涉及利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用。
背景技术
原油作为炼化企业的最主要原料,一方面,原油的需求量剧增、进口量扩大、价格居高不下且波动频繁;另一方面,原油产品存在着性质劣质化、种类丰富、同名油前后性质有差异、装置进料要求高、混兑油性质难掌握等诸多特点。这些给炼化企业带来了巨大的压力。及时获得当前原油的性质评价数据——即原油快速评价,将为原油贸易、原油输送、原油调合、原油加工、全厂生产计划、生产调度等生产过程优化提供支撑。原油评价包含指标众多,如密度、残炭、酸值、硫含量、氮含量、蜡含量、沥青质含量和实沸点蒸馏曲线等。采用传统的评价方法,存在分析时间长、处理繁琐、仪器要求高、劳动强度大等现象,已不能满足实际应用的需求。
近红外分析技术是目前最有前景且应用最广泛的快速分析方法之一。近几年来光纤在近红外光谱技术领域的应用使近红外光谱技术从实验室走向现场,光纤化学和热稳定性、对电磁干扰不敏感、传输信号能量集中、灵敏度高、价格低廉等优点,使得近红外光谱仪可以在恶劣、危险的环境中进行远距离快速在线分析。近红外光谱仪的光纤测样附件,按照测样的形式可分为透射器件和漫反射器件;从样品的引入方式来看,有探头型和流通池型。基于光纤结构的近红外探头结构非常灵活,可按硬件要求设计。从现场应用效果看,光纤结构的近红外探头越来越受到急需远程和在线分析的各行业的青睐。
原油组分复杂,属于黏稠深色液体。原油的待测性质多,并且其近红外光谱吸收带较宽且重叠严重。在实际测量中,光谱分析系统探头的结构非常关键。光纤探头的性质和结构不同对测量信噪比有很大的影响。光纤探头的窗片或透镜表面被污染则会影响光通量使灵敏度降低,测试过程中有外来光的干扰则会使检测的信噪比和灵敏度下降。常规透射式光纤探头在测量深色原油时携带样品信息不够,并且在实际应用中流通池易被黏稠原油粘附,进而造成样品谱图失真,模型预测精度低、现场仪器维护工作量大等问题,影响实际投用效果。
因设计原理的不同,反射式光纤探头结合近红外光谱仪,在原油的近红外光谱检测上将一定程度上改善现有的分析不足。根据通过光纤技术远程采集信号,建立原油的近红外反射光谱数据库,利用光谱预处理技术以及近红外建模技术,可以快速获取原油性质,就有可能成为原油等深色重质油品物化性质的在线快速测定的一种良好手段。同时,由于近红外分析仪是二次测量仪表,即近红外分析仪并不是直接测量物质性质,必须先建立待测物质的属性与近红外光谱之间的数学模型然后根据模型来测量物质属性。因此,可以预想到,一种兼顾实用性、实时性、稳定性和良好预测精度的原油快速评价方法的发明,将倍受青睐。
发明内容
鉴于上述问题,本发明提出了一种利用反射式近红外光纤探头构建快速预测原油性质的模型以及利用该模型快速预测原油性质的方法。该方法通过选取合适型号的反射式近红外光纤探头和离线/在线近红外分析仪,采用将光纤探头直接插入原油样品的测量方式,快速获取原油近红外光谱图。通过一阶微分对原油近红外光谱图进行预处理,采用主成分分析方法(PCA)对预处理后的光谱样本进行筛选,利用偏最小二乘法(PLS)建立的原油定量分析模型,基于此模型可实现未知原油性质的快速预测分析。该方法与其他的原油性质测量方法相比,无需复杂的样品预处理,具有操作简单、探头维护量小、近红外光谱信噪比高、定量分析模型精度高等特点,可快速预测原油的性质,在工业在线应用时具有较好的前景。
本发明提供的构建基于近红外反射式光纤探头检测原油性质的模型的方法包括以下步骤:
步骤一:构建原油校正集,测定校正集中原油的性质;
步骤二:利用反射式光纤探头采集原油近红外光谱;
步骤三:利用一阶微分对步骤二获得的原油近红外光谱进行预处理;
步骤四:利用主成分分析方法对预处理后的光谱数据进行训练集样本的选择;
步骤五:根据原油待测性质和原油光谱数据集,选择特定的光谱频段;和
步骤六:利用偏最小二乘法建立原油性质与近红外光谱数据之间的数学关联模型。
在一个或多个实施方案中,步骤一中,用于构建校正集的原油20℃的密度在0.7-1.1g/cm3的范围内,硫含量在0.03%-5.50%的范围内,酸值在0.01-12.00mgKOH/g的范围内。
在一个或多个实施方案中,所述原油性质包括密度、残炭、酸值、硫含量、氮含量、蜡含量、胶质含量、沥青质含量和实沸点数据中的一个或多个。
在一个或多个实施方案中,步骤二包括,将校正集样品放置于30℃温度下(维持恒定),待原油样品温度达到稳定状态后,测定该原油样品的近红外光谱数据。
在一个或多个实施方案中,步骤二中,采用的反射式近红外光纤探头含有14根光纤,7根输出,7根输入,探头外部加装有一个底部由蓝宝石材质制成光窗材料的金属套筒,套筒可拆卸。
在一个或多个实施方案中,步骤二中,利用配置反射式近红外光纤探头的离线/在线近红外分析仪,将该反射式近红外光纤探头直接插入原油,探头外部套筒底的蓝宝石光窗材料被原油浸没即可的方式,测得原油近红外光谱数据。
在一个或多个实施方案中,步骤二中,扫描范围为4000-12500cm-1,扫描次数为10-400次。
在一个或多个实施方案中,步骤三包括,利用一阶微分对步骤二获得的波数范围为12500~4000cm-1区域的原油样本近红外光谱图进行预处理,消除基线和背景干扰,建立初始训练集。
在一个或多个实施方案中,步骤三所述的预处理为S-G一阶导数,用以消除背景干扰。
在一个或多个实施方案中,步骤四包括,对样本光谱进行主成分(PCA)分析后,利用主成分得分作为特征变量,再利用样本间的欧式距离来选择样本,确定最终的训练集。
在一个或多个实施方案中,步骤四包括,采用剔除的方式对校正样本进行选择,其过程为:首先对样本光谱进行主成分(PCA)分析,然后利用主成分得分作为特征计算每个样本与邻近样本之间的欧式距离,并根据样本分布的密度确定阈值;对于每一个样本,剔除与其距离小于阈值的样本,从而剔除冗余样本,剩余的样本作为校正样本;部分样本之间的欧式距离过小造成样本冗余,通过设定阈值并且剔除与其距离小于阈值的样本的方法来消除冗余样本。
在一个或多个实施方案中,步骤四所述的所有样本两两之间的欧式距离dij按下式计算:
式中,dij:第i个样本和第j个样本之间的欧式距离;
m:为主成分数;
xi:为第i个样本的主成分得分;
xj:为第j个样本的主成分得分。
在一个或多个实施方案中,步骤五中,对于密度选择波数范围为4247-6107cm-1、对于残炭选择波数范围4599-5456cm-1和6099-7500cm-1、对于酸值选择波数范围4599-7500cm-1、对于硫含量选择波数范围4599-7500cm-1、对于氮含量选择波数范围4500-6600cm-1、对于蜡含量选择波数范围4500-6600cm-1、对于胶质含量选择波数范围4500-6600cm-1、对于沥青质含量选择波数范围4500-6600cm-1和对于实沸点蒸馏数据选择波数范围为4599-7500cm-1。
在一个或多个实施方案中,步骤六中,所述原油性质包括密度、残炭、酸值、硫含量、氮含量、蜡含量、胶质含量、沥青质含量和实沸点数据中的一个或多个。
在一个或多个实施方案中,步骤六中,所述数学关联模型如下式所示:
y=a0+a1x1+a2x2+…+anxn
其中,y为预测的性质,ai为模型参数,xi为光谱第i个波数点的吸光度。
在一个或多个实施方案中,所述数学关联模型如下建立:
(1)对光谱矩阵X和浓度矩阵Y作标准化变换,变换后的矩阵分别记为V和U;
(2)计算V矩阵的权重ω′=u′V/u′u;
(3)对权重向量进行归一化ω′new=ωold′/(ωold′ωold)1/2;
(4)估计V矩阵得分向量t=Vω′;
(5)计算U矩阵的载荷q′=t′U/t′t;
(6)产生U矩阵的得分向量u=Uq/q′q;
比较u(old)与u(new),如果||u(old)-u(new)||<阈值,表明以收敛,迭代停止,否则转到第一步继续进行迭代;
(7)计算标量b用以内部关联b=u′t/t′t;
(8)计算V矩阵的载荷p′=t′v/t′t;
(9)计算V和U矩阵的残差E=V-tp′,F=U-uq′。
(10)计算预测标准差SEP,如果SEP大于预期精度,则表明最佳维数已得到,否则对下一维进行计算,则可得到最终的系数矩阵
B=W(P'W)-1Q'。
本发明还包括一种利用反射式近红外光纤探头快速预测原油性质的方法,所述方法包括:
步骤A:利用反射式光纤探头采集原油样品的近红外光谱;
步骤B:利用一阶微分对步骤A获得的原油近红外光谱进行预处理;
步骤C:根据原油待测性质和原油光谱数据集,选择特定的光谱频段;和
步骤D:利用以下数学关联模型预测该原油样品的相关性质:
y=a0+a1x1+a2x2+…+anxn
其中,y为预测的性质,ai为模型参数,xi为光谱第i个波数点的吸光度。
在一个或多个实施方案中,所述原油性质包括密度、残炭、酸值、硫含量、氮含量、蜡含量、胶质含量、沥青质含量和实沸点数据中的一个或多个。
在一个或多个实施方案中,步骤A包括,将原油样品放置于30℃温度下,并维持温度恒定,待原油样品温度达到稳定状态后,测定该原油样品的近红外光谱数据。
在一个或多个实施方案中,步骤A中,采用的反射式近红外光纤探头含有14根光纤,7根输出,7根输入,探头外部加装有一个底部由蓝宝石材质制成光窗材料的金属套筒,套筒可拆卸。
在一个或多个实施方案中,步骤A中,利用配置反射式近红外光纤探头的离线/在线近红外分析仪,将该反射式近红外光纤探头直接插入原油,探头外部套筒底的蓝宝石光窗材料被原油浸没即可的方式,测得原油近红外光谱数据。
在一个或多个实施方案中,步骤A中,扫描范围为4000-12500cm-1,扫描次数为10-400次。
在一个或多个实施方案中,步骤B包括,利用一阶微分对步骤A获得的波数范围为12500~4000cm-1区域的原油样本近红外光谱图进行预处理,消除基线和背景干扰,建立初始训练集。
在一个或多个实施方案中,步骤B所述的预处理为S-G一阶导数,用以消除背景干扰。
在一个或多个实施方案中,步骤C中,对于密度选择波数范围为4247-6107cm-1、对于残炭选择波数范围4599-5456cm-1和6099-7500cm-1、对于酸值选择波数范围4599-7500cm-1、对于硫含量选择波数范围4599-7500cm-1、对于氮含量选择波数范围4500-6600cm-1、对于蜡含量选择波数范围4500-6600cm-1、对于胶质含量选择波数范围4500-6600cm-1、对于沥青质含量选择波数范围4500-6600cm-1和对于实沸点蒸馏数据选择波数范围为4599-7500cm-1。本发明的有益效果如下:
本发明测试方式简单、快速、实用,利用近红外分析仪并配置含有14根光纤的反射式近红外探头,将配套可拆卸的底部带有蓝宝石光窗材料的金属套筒的反射式近红外光纤探头直接插入待测点的原油样品中,利用近红外光谱仪快速测定原油性质。与传统的测量方法相比,本发明方法大大缩短了检测时间,减少了人力、物力。测试过程中无需使用任何试剂对原油样品处理,不损坏样品;与其他近红外测量方式相比,不需要把样品从待测点取出送到近红外分析仪进而避免了繁琐的原油预处理系统,仅通过光纤的加长,即可实现原油样品原位、实时分析。同时,带蓝宝石光窗材料金属套筒的反射式探头的使用,可以有效减少原油黏稠而粘附在光纤探头上影响样品实时性的现象。在此基础上,本发明采用的综合建模方法,即,基于采用一阶微分对采集到的原油样本近红外光谱图进行预处理,通过主成分分析技术筛选样本,并选择原油光谱图合适的波数范围,利用偏最小二乘法建立原油属性值与其近红外光谱数据之间的数学模型,可实现未知原油属性值的快速预测分析。利用该方法,预处理后的近红外光谱信噪比高,建立的模型精度高,可以检测原油密度、残炭、酸值、硫含量、氮含量、蜡含量、胶质含量、沥青质含量和实沸点蒸馏数据。该发明可为原油性质监控、储运、调合、常减压蒸馏装置操作等各项与原油有关业务的优化提供支撑。
附图说明
图1:基于反射式近红外光纤探头检测原油样品近红外光谱数据。(a)实验过程;(b)近红外光的路线示意图。
图2:基于反射光纤探头快速预测原油性质方法的总流程图。
图3:原始的原油近红外光谱图。
图4:预处理后的原油近红外谱图。
图5:PCA分析主成分。
图6:剔除异常点的样本点图。
图7:近红外原油API回归模型。
图8:近红外原油硫含量回归模型。
具体实施方式
图1显示本发明基于反射式近红外光纤探头检测样品近红外光谱数据实验过程。图2为本发明预测原油性质的总流程图,具体包括以下步骤:
(1)构建原油校正集,利用标准分析方法测定样品相关性质;
(2)利用近红外分析仪和反射式探头,将配套可拆卸的底部带有蓝宝石光窗材料的金属套筒的反射式近红外光纤探头,直接插入待测点的原油样品中测得近红外光谱数据;
(3)利用一阶微分对光谱进行预处理;
(4)利用主成分分析技术对预处理后的光谱图进行特征变量筛选,并剔除异常样本;
(5)根据建模项目,确定近红外波数范围;和
(6)利用偏最小二乘法建立性质校正模型。
图2的流程图中还显示了利用所建的模型,结合步骤(2)、(3)、(5),测试未知原油性质的步骤。
下文将对这些步骤进行详细描述。应理解,在本发明范围内中,本发明的上述各技术特征和在下文(如实施例)中具体描述的各技术特征之间都可以互相组合,从而构成优选的技术方案。
一、构建原油校正集,测定校正集中原油的性质
可收集不同种类的原油样本,通常覆盖石蜡基原油、中间基原油和环烷基原油等。通常,所收集的原油样本数量不少于50个。较好的是,所收集的原油样本的密度(20℃)、硫含量和酸值指标分别控制在0.7~1.1g/cm3、0.03%~5.50%和0.01~12.00mgKOH/g的范围之内。然后利用传统的标准方法测得所收集原油的多个性质属性,如密度、残炭、氮含量、硫含量、酸值、盐含量、蜡含量、胶质含量、沥青质含量和实沸点蒸馏数据等,并记录数据。
二、采集原油近红外光谱
可选取合适型号的离线或在线近红外光谱仪,配套可拆卸的底部带有蓝宝石光窗材料的金属套筒的反射式近红外光纤探头进行近红外光谱扫描,采用将反射式近红外光纤探头直接插入原油样品(原油样品温度维持在30℃以下的某个恒定温度)的测量方式,获得每份样本的近红外光谱图。例如,可将原油样品放置于30℃温度下,并维持温度恒定,待原油样品温度达到稳定状态后,测定该原油样品的近红外光谱数据。
本发明采用的反射式近红外光纤探头可含有14根光纤,7根输出,7根输入,探头外部可加装有一个底部由蓝宝石材质制成光窗材料的金属套筒,套筒可拆卸。
每张光谱图扫描时间为100-400次,取平均值。光谱扫描范围为4000-12500cm-1,分辨率16-32cm-1。示例性的原油预处理光谱见图3。
三、利用一阶微分对步骤二获得的原油近红外光谱进行预处理
该预处理包括对校正集每份样品的12500-4000cm-1的谱区进行一阶微分处理,消除基线和背景干扰,提高分辨率和灵敏度。预处理后,可建立初始训练集。示例性的预处理后的原油近红外谱图见图4。
四、利用主成分分析技术对预处理后的光谱图进行特征变量筛选,并剔除异常样本
本步骤是对步骤三中预处理过后的光谱进行校正样本的选择,从样本库中选择代表性强的样本建立模型,这样不仅可以提高模型的建立速度、减少模型库的储存空间,更为重要的是它便于模型的更新和维护。
在建立训练集之前,应首先剔除异常样本,这些异常样本可能含有异常化学组分或组分浓度较为极端,与其他样本存在显著差异。若这些样本点参与模型的建立,会影响分析模型的准确性和稳健性。
在本步骤中采用剔除的方式对训练样本进行选择,其基本过程为首先对样本光谱进行主成分(PCA)分析后,利用主成分得分(示例性的PCA分析主成分得分可见图5)作为特征计算每个样本与邻近样本之间的欧式距离,并根据样本分布的密度确定阈值。对于每一个样本,剔除与其距离小于阈值的样本,从而剔除冗余样本,剩余的样本作为训练样本。部分样本之间的欧式距离过小造成样本冗余,通过设定阈值并且剔除与其距离小于阈值的样本的方法来消除冗余样本。示例性的剔除异常样本点后训练集可见图6。
上述中所有样本两两之间的欧式距离dij按下式计算:
式中,dij:第i个样本和第j个样本之间的欧式距离;
m:为主成分数;
xi:为第i个样本的主成分得分;
xj:为第j个样本的主成分得分。
五、根据建模项目,确定近红外波数范围
本步骤对训练集中的光谱样本进行波数选择。随着对偏最小二乘等方法的深入研究,发现通过筛选特征波数或区间有可能得到更好地定量模型。通过波数选择可以简化模型,并且通过波数选择可以剔除不相关的变量,得到预测能力更强,稳健性更好的模型。
通常,波数选择范围为4000-12500cm-1之内的任何有限波数范围,可为多个波数段的组合。在本发明的某些实施方案中,对于密度选择波数范围为4247-6107cm-1、对于残炭选择波数范围4599-5456cm-1和6099-7500cm-1、对于酸值选择波数范围4599-7500cm-1、对于硫含量选择波数范围4599-7500cm-1、对于氮含量选择波数范围4500-6600cm-1、对于蜡含量选择波数范围4500-6600cm-1、对于胶质含量选择波数范围4500-6600cm-1、对于沥青质含量选择波数范围4500-6600cm-1和对于实沸点蒸馏选择波数范围为4599-7500cm-1。
六、利用偏最小二乘法建立原油性质与近红外光谱数据之间的数学关联模型
偏最小二乘回归与主成分回归相比,不仅仅考虑了光谱矩阵,同时也考虑了浓度矩阵的影响。此步骤利用经过预处理和波数选择的训练集中的原油样本的近红外光谱图与性质数值建立模型。该方法的数学关联模型如下式所示:
y=a0+a1x1+a2x2+…+anxn
其中,y为预测的性质,ai为模型参数,xi为光谱第i个波数点的吸光度。
在一个或多个实施方案中,所述数学关联模型如下建立:
(1)对光谱矩阵X和浓度矩阵Y作标准化变换,变换后的矩阵分别记为V和U;
(2)计算V矩阵的权重ω′=u′V/u′u;
(3)对权重向量进行归一化ω′new=ωold′/(ωold′ωold)1/2;
(4)估计V矩阵得分向量t=Vω′;
(5)计算U矩阵的载荷q′=t′U/t′t;
(6)产生U矩阵的得分向量u=Uq/q′q;
比较u(old)与u(new),如果||u(old)-u(new)||<阈值,表明以收敛,迭代停止,否则转到第一步继续进行迭代;
(7)计算标量b用以内部关联b=u′t/t′t;
(8)计算V矩阵的载荷p′=t′v/t′t;
(9)计算V和U矩阵的残差E=V-tp′,F=U-uq′。
(10)计算预测标准差SEP,如果SEP大于预期精度,则表明最佳维数已得到,否则对下一维进行计算,则可得到最终的系数矩阵
B=W(P′W)-1Q′。
本发明在预测待测原油样本的性质时,首先采用本发明步骤二中所述的方法测定待测原油样本的近红外光谱图,然后采用步骤三所述的方法对待测原油样本的近红外光谱图进行预处理,之后根据步骤五选择的波数范围选择变量,并利用步骤六中所建立的定量分析模型预测待测原油的相关性质。
下面通过实施例对本发明进行具体描述。有必要在此指出的是,以下实施例只用于对本发明作进一步说明,不能理解为对本发明保护范围的限制,该领域的专业技术人员根据本发明的内容作出的一些非本质的改进和调整,仍属于本发明的保护范围。
实施例1
以下以硫含量预测的实施例来说明本发明具体步骤包括:
步骤一:采集不同种类的原油样本156个,122个作为校正集,34个作为验证集。
步骤二:样品温度控制在30℃,选用BRUKER布鲁克近红外光谱仪,以及配套可拆卸的底部带有蓝宝石光窗材料金属套筒的BRUKER布鲁克反射式近红外光纤探头,进行试验测定。通过将近红外光纤探头直接插入各原油样本的方式,测定原油样本的近红外光谱,光谱范围扫描范围为4000-12500cm-1,分辨率32cm-1,累计扫描次数400次。并且按照传统标准方法测量原油样本的硫含量和API。图3为原始的原油近红外光谱图。可以看到,原始光谱的基线漂移严重,谱峰重叠严重。
步骤三:选取8000-4000cm-1谱区范围的吸光度,对其进行一阶微分预处理,建立原油样本近红外光谱矩阵。图4为预处理之后的谱图。
步骤四:对预处理后的原油样本采用剔除的方式进行训练样本的选择,首先对预处理后的原油样本光谱进行主成分分析后,利用主成分得分(图5)作为特征计算每个样本与邻近样本之间的欧式距离,并对于每一个样本,比较其与邻近样本之间的欧氏距离与阈值1.8,剔除与其距离小于阈值的样本,本例中剔除样本14,27,35,86,91,101,111,117,122从而剔除冗余样本,剩余的样本作为训练样本。最终,选中103个训练样本构成原油光谱训练样本集(图6)。
步骤五:根据待测项目为硫含量,优选的波数区间为4599-7500cm-1;根据待测项目为API,优选的波数区间为4247-6107cm-1。
步骤六:运用偏最小二乘法建立硫含量值与近红外光谱的回归模型,预测属性与近红外光谱之间的关系如下式
y=a0+alxl+a2x2+…+anxn
其中:y为预测属性,ai为模型参数,xi为光谱第i个波数点的吸光度。
所建的硫含量回归模型如图7,所建的API回归模型如图8。利用验证集对模型进行验证:对待测原油样本首先根据步骤二测得近红外光谱,根据步骤三中的方法,对8000-4000cm-1谱区内对其进行一阶微分预处理,之后选择4599-7500cm-1(硫含量)或4247-6107cm-1(API)的波数区间,最后利用步骤六中建立的模型对其进行预测。硫含量模型的决定系数达0.9228,交互验证均方误差为0.373。API模型的决定系数达0.9674,交互验证均方误差为1.36预测值与实际值的比较结果如下表1所示,预测过程快速、简单,预测结果准确。
表1:原油API预测值与实际值结果对比
表2:原油硫含量预测值与实际值结果对比
Claims (10)
1.一种利用反射式近红外光纤探头构建快速预测原油性质的模型方法,其特征在于,所述方法包括以下步骤:
步骤一:构建原油校正集,测定校正集中原油的性质;
步骤二:利用反射式光纤探头采集原油近红外光谱;
步骤三:利用一阶微分对步骤二获得的原油近红外光谱进行预处理;
步骤四:利用主成分分析方法对预处理后的光谱数据进行训练集样本的选择;
步骤五:根据原油待测性质和原油光谱数据集,选择特定的光谱频段;
步骤六:利用偏最小二乘法建立原油性质与近红外光谱数据之间的数学关联模型。
2.如权利要求1所述的方法,其特征在于,
所述步骤一中用于构建校正集的原油20℃的密度在0.7-1.1g/cm3的范围内,硫含量在0.03%-5.50%的范围内,酸值在0.01-12.00mgKOH/g的范围内;和/或
所述原油性质包括密度、残炭、酸值、硫含量、氮含量、蜡含量、胶质含量、沥青质含量和实沸点数据中的一个或多个。
3.如权利要求1-2中任一项所述的方法,其特征在于,所述步骤二包括,将校正集样品放置于30℃温度下,维持温度恒定,待原油样品温度达到稳定状态后,测定该原油样品的近红外光谱数据。
4.如权利要求1-3中任一项所述的方法,其特征在于,所述步骤二中,利用配置反射式近红外光纤探头的离线/在线近红外分析仪,将该反射式近红外光纤探头直接插入原油,探头外部套筒底的蓝宝石光窗材料被原油浸没即可的方式,测得原油近红外光谱数据;
优选地,所述反射式近红外光纤探头含有14根光纤,7根输出,7根输入,探头外部加装有一个底部由蓝宝石材质制成光窗材料的金属套筒,套筒可拆卸;
优选地,检测时的扫描范围为4000-12500cm-1,扫描次数为10-400次。
5.如权利要求1-4中任一项所述的方法,其特征在于,所述步骤三包括,利用一阶微分对所述步骤二获得的波数范围为12500~4000cm-1区域的原油样本近红外光谱图进行预处理,消除基线和背景干扰,建立初始训练集;
优选地,所述预处理为S-G一阶导数。
6.如权利要求1-5中任一项所述的方法,其特征在于,所述步骤四包括,对样本光谱进行主成分(PCA)分析后,利用主成分得分作为特征变量,再利用样本间的欧式距离来选择样本,确定最终的训练集;
优选地,步骤四包括,采用剔除的方式对校正样本进行选择,其过程为:首先对样本光谱进行主成分(PCA)分析,然后利用主成分得分作为特征计算每个样本与邻近样本之间的欧式距离,并根据样本分布的密度确定阈值;对于每一个样本,剔除与其距离小于阈值的样本,从而剔除冗余样本,剩余的样本作为校正样本;部分样本之间的欧式距离过小造成样本冗余,通过设定阈值并且剔除与其距离小于阈值的样本的方法来消除冗余样本;
更优选地,所述所有样本两两之间的欧式距离dij按下式计算:
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
式中,dij:第i个样本和第j个样本之间的欧式距离;
m:为主成分数;
xi:为第i个样本的主成分得分;
xj:为第j个样本的主成分得分。
7.如权利要求1-6中任一项所述的方法,其特征在于,所述步骤五中,对于密度选择波数范围为4247-6107cm-1、对于残炭选择波数范围4599-5456cm-1和6099-7500cm-1、对于酸值选择波数范围4599-7500cm-1、对于硫含量选择波数范围4599-7500cm-1、对于氮含量选择波数范围4500-6600cm-1、对于蜡含量选择波数范围4500-6600cm-1、对于胶质含量选择波数范围4500-6600cm-1、对于沥青质含量选择波数范围4500-6600cm-1、对于实沸点蒸馏数据选择波数范围为4599-7500cm-1。
8.如权利要求1-7中任一项所述的方法,其特征在于,所述步骤六的数学关联模型如下式所示:
y=a0+a1x1+a2x2+…+anxn
其中,y为预测的性质,ai为模型参数,xi为光谱第i个波数点的吸光度。
9.如权利要求8所述的方法,其特征在于,所述数学关联模型如下建立:
(1)对光谱矩阵X和浓度矩阵Y作标准化变换,变换后的矩阵分别记为V和U;
(2)计算V矩阵的权重ω′=u′V/u′u;
(3)对权重向量进行归一化ω′new=ωold′/(ωold′ωold)1/2;
(4)估计V矩阵得分向量t=Vω′;
(5)计算U矩阵的载荷q′=t′U/t′t;
(6)产生U矩阵的得分向量u=Uq/q′q;
比较u(old)与u(new),如果||u(old)-u(new)||<阈值,表明以收敛,迭代停止,否则转到第一步继续进行迭代;
(7)计算标量b用以内部关联b=u′t/t′t;
(8)计算V矩阵的载荷p′=t′v/t′t;
(9)计算V和U矩阵的残差E=V-tp′,F=U-uq′。
(10)计算预测标准差SEP,如果SEP大于预期精度,则表明最佳维数已得到,否则对下一维进行计算,则可得到最终的系数矩阵
B=W(P′W)-1Q′。
10.一种利用反射式近红外光纤探头快速预测原油性质的方法,其特征在于,所述方法包括:
步骤A:利用反射式光纤探头采集原油样品的近红外光谱;
步骤B:利用一阶微分对步骤A获得的原油近红外光谱进行预处理;
步骤C:根据原油待测性质和原油光谱数据集,选择特定的光谱频段;和
步骤D:利用以下数学关联模型预测该原油样品的相关性质:
y=a0+a1x1+a2x2+…+anxn
其中,y为预测的性质,ai为模型参数,xi为光谱第i个波数点的吸光度;
优选地,所述步骤A如权利要求3或4所述;所述步骤B如权利要求5所述;所述步骤C如权利要求7所述。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710985799.8A CN107860743A (zh) | 2017-10-20 | 2017-10-20 | 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710985799.8A CN107860743A (zh) | 2017-10-20 | 2017-10-20 | 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107860743A true CN107860743A (zh) | 2018-03-30 |
Family
ID=61697620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710985799.8A Pending CN107860743A (zh) | 2017-10-20 | 2017-10-20 | 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107860743A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958267A (zh) * | 2017-11-21 | 2018-04-24 | 东南大学 | 一种基于光谱线性表示的油品性质预测方法 |
CN111398208A (zh) * | 2020-04-22 | 2020-07-10 | 谭贵平 | 一种利用近红外技术对传统工艺翡翠进行快速鉴定的方法 |
CN112149858A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 一种利用原油谱学信息预测渣油馏分sara组成的方法 |
CN116343933A (zh) * | 2022-11-22 | 2023-06-27 | 中石油云南石化有限公司 | 对原油油段性质进行跟踪、控制并重定义的方法及用途 |
CN116413236A (zh) * | 2023-02-27 | 2023-07-11 | 西南石油大学 | 一种用于钻井返出液总烃含量的检测装置及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893561A (zh) * | 2010-07-13 | 2010-11-24 | 中国人民解放军总后勤部油料研究所 | 一种润滑油新油质量近红外光谱快速测定方法 |
CN104897607A (zh) * | 2015-06-18 | 2015-09-09 | 北京工商大学 | 便携式近红外光谱食品建模与快速检测一体化方法和系统 |
CN106066215A (zh) * | 2016-07-29 | 2016-11-02 | 武汉理工大学 | 一种蓝宝石高温传感器 |
-
2017
- 2017-10-20 CN CN201710985799.8A patent/CN107860743A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893561A (zh) * | 2010-07-13 | 2010-11-24 | 中国人民解放军总后勤部油料研究所 | 一种润滑油新油质量近红外光谱快速测定方法 |
CN104897607A (zh) * | 2015-06-18 | 2015-09-09 | 北京工商大学 | 便携式近红外光谱食品建模与快速检测一体化方法和系统 |
CN106066215A (zh) * | 2016-07-29 | 2016-11-02 | 武汉理工大学 | 一种蓝宝石高温传感器 |
Non-Patent Citations (2)
Title |
---|
王玉: "《药品检验》", 31 October 2011, 中国医药科技出版社 * |
陈兰珍 等: "《蜂蜜近红外光谱技术》", 31 January 2012 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107958267A (zh) * | 2017-11-21 | 2018-04-24 | 东南大学 | 一种基于光谱线性表示的油品性质预测方法 |
CN107958267B (zh) * | 2017-11-21 | 2021-04-27 | 东南大学 | 一种基于光谱线性表示的油品性质预测方法 |
CN112149858A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 一种利用原油谱学信息预测渣油馏分sara组成的方法 |
CN111398208A (zh) * | 2020-04-22 | 2020-07-10 | 谭贵平 | 一种利用近红外技术对传统工艺翡翠进行快速鉴定的方法 |
CN111398208B (zh) * | 2020-04-22 | 2023-10-31 | 谭贵平 | 一种利用近红外技术对传统工艺翡翠进行快速鉴定的方法 |
CN116343933A (zh) * | 2022-11-22 | 2023-06-27 | 中石油云南石化有限公司 | 对原油油段性质进行跟踪、控制并重定义的方法及用途 |
CN116343933B (zh) * | 2022-11-22 | 2023-09-01 | 中石油云南石化有限公司 | 对原油油段性质进行跟踪、控制并重定义的方法及用途 |
CN116413236A (zh) * | 2023-02-27 | 2023-07-11 | 西南石油大学 | 一种用于钻井返出液总烃含量的检测装置及方法 |
CN116413236B (zh) * | 2023-02-27 | 2024-03-22 | 西南石油大学 | 一种用于钻井返出液总烃含量的检测装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107860743A (zh) | 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 | |
CN107817223A (zh) | 快速无损实时预测原油性质模型的构建方法及其应用 | |
CN105300923B (zh) | 一种近红外光谱分析仪在线应用时无测点温度补偿模型修正方法 | |
CN107703097A (zh) | 利用衰减全反射探头与近红外光谱仪构建快速预测原油性质的模型的方法及其应用 | |
CN101706307B (zh) | 基于透射光谱的输油管道油品界面的检测方法及其装置 | |
CN105486658B (zh) | 一种具有无测点温度补偿功能的近红外物性参数测量方法 | |
CN109668858A (zh) | 基于近红外光谱检测发酵过程生物量和组分浓度的方法 | |
CN101413885A (zh) | 一种快速定量蜂蜜品质的近红外光谱方法 | |
CN104897607A (zh) | 便携式近红外光谱食品建模与快速检测一体化方法和系统 | |
CN101769866B (zh) | 一种牛奶成分的检测装置和方法 | |
CN102590129B (zh) | 近红外检测花生中氨基酸含量的方法 | |
CN106932378A (zh) | 一种基于拉曼光谱的炼厂酸性气成分的在线检测系统和方法 | |
CN109324013A (zh) | 一种利用高斯过程回归模型构建原油性质近红外快速分析的方法 | |
CN102004097A (zh) | 一种基于主导因素结合偏最小二乘法的煤质在线检测方法 | |
CN109324014A (zh) | 一种自适应的原油性质近红外快速预测方法 | |
CN102306236B (zh) | 基于角度度量的多变量分析方法 | |
CN109738380A (zh) | 一种土壤盐渍化程度的高光谱遥感判断方法 | |
CN103134770B (zh) | 消除水分对近红外光谱检测土壤全氮含量影响的方法 | |
CN106442474A (zh) | 一种基于偏最小二乘法的水泥生料三率值测量方法 | |
CN101349638A (zh) | 果蔬维生素c含量的光谱快速无损检测方法 | |
CN108693139A (zh) | 电子烟油理化指标的近红外光谱预测模型建立方法及应用 | |
CN104849234A (zh) | 基于近红外光谱分析吡虫啉原药主成分含量的测定方法 | |
CN103063599B (zh) | 由透射红外光谱预测原油密度的方法 | |
CN103134764B (zh) | 由透射红外光谱预测原油实沸点蒸馏曲线的方法 | |
CN103105369A (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 |
Application publication date: 20180330 |
|
RJ01 | Rejection of invention patent application after publication |