CN109632691B - 一种原油精细物性的近红外快速分析方法 - Google Patents

一种原油精细物性的近红外快速分析方法 Download PDF

Info

Publication number
CN109632691B
CN109632691B CN201811629590.9A CN201811629590A CN109632691B CN 109632691 B CN109632691 B CN 109632691B CN 201811629590 A CN201811629590 A CN 201811629590A CN 109632691 B CN109632691 B CN 109632691B
Authority
CN
China
Prior art keywords
crude oil
sample
detected
samples
physical properties
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811629590.9A
Other languages
English (en)
Other versions
CN109632691A (zh
Inventor
陈夕松
蒋立沫
苏曼
周文娟
江广兰
梅彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NANJING RICHISLAND INFORMATION ENGINEERING CO LTD
Original Assignee
NANJING RICHISLAND INFORMATION ENGINEERING CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NANJING RICHISLAND INFORMATION ENGINEERING CO LTD filed Critical NANJING RICHISLAND INFORMATION ENGINEERING CO LTD
Priority to CN201811629590.9A priority Critical patent/CN109632691B/zh
Publication of CN109632691A publication Critical patent/CN109632691A/zh
Application granted granted Critical
Publication of CN109632691B publication Critical patent/CN109632691B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3577Investigating 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating 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

一种原油精细物性的近红外快速分析方法
技术领域
本发明涉及石油化工领域的原油物性分析方法,具体为一种采用近红外技术快速分析原油中烃族组成这类精细物性的方法。
背景技术
近红外光谱分析技术可实现原油性质的快速检测,但现有原油性质检测中,关注的还是原油整体物理性质,如石脑油收率、硫含量、酸值等等。然而,近年来随着炼油技术的发展,研究发现即使相同的整体物性,由于详细的化合物分子类型和碳分布等不同,其加工过程中反应和转化规律差别很大,进而最佳加工路线迥异。以原油组分中价值最高的石脑油为例,其所含的正构烷烃宜作为化工乙烯原料,而非正构烷烃(主要是异构烷烃、环烷烃和芳烃)宜作为炼油重整原料进行加工。因此,为保证乙烯装置和重整装置都能在额定负荷下生产,并最大程度发挥石脑油的利用效率,在即使分析清楚原油中石脑油收率的情况下,还应进一步弄清其正构烷烃及非正构烷烃的准确含量。
上述含量虽然可以在实验室内采用气相色谱等技术化验检测,但周期长,费用高。近红外光谱分析技术结合常规的化学计量学方法,如偏最小二乘PLS,虽然可以进行快速分析,但PLS这类常规方法对性质分布的均匀性、建模样本数量等都有严格要求,否则建模预测精度较差。而原油是自然形成的,难以做到样本性质均匀分布,部分原油与样本库中其它原油性质甚至有很大差别,表现为与其相似的原油样本数量很少,难以采用PLS建模。因此,需要研发一种快速分析技术,实现对原油烃族组成这类精细物性的快速、准确分析检测。
发明内容
本发明提出一种原油精细物性的近红外快速分析方法,针对原油烃族组成这类精细物性,采用加权距离方法建立近红外预测模型,对待测样本进行快速分析,具体包括以下步骤:
1)获取待测原油样本的近红外光谱数据;
2)对待测原油样本和原油库中样本的光谱数据进行基线校正、波数段截取和矢量归一化;
3)采用主成分分析方法计算吸光度矩阵的得分矩阵,提取累积贡献率达85%对应的得分向量组成新的得分矩阵;
4)基于新的得分矩阵,以待测样本为中心建立边长比例固定的框,进行相似样本选取;
5)统计框内相似样本数量num,设定相似样本阈值M,对待测样本的不确定度进行判断,依据规则选取相似样本建立校正集;
6)基于选取的校正集样本,采用加权距离方法建立原油精细物性预测模型,并对待测样本进行快速分析。
优选的,原油的精细物性指原油中所含链烷烃、烯烃、环烷烃和芳烃这类烃族组成的质量分数。
优选的,相似样本边界框的建立规则,根据主成分分析的结果选取前n个主成分,建立比例为n:(n-1):……:1的n维框,其中n为步骤3)中满足累积贡献率的主成分数量。
具体的,步骤5)中校正集样本的建立规则:
若相似样本数量num≥M,则该待测样本的不确定度为1,依据与待测样本的距离由近及远选取M个相似样本组成校正集;
若相似样本数量0.5M≤num<M,则该待测样本的不确定度为2,选取num个相似样本组成校正集;
若相似样本数量num<0.5M,则该待测样本的不确定度为3,不进行建模预测。
具体的,步骤6)中,采用加权距离方法计算原油精细物性:
Figure GDA0002818663120000021
式中,p为待测原油的精细物性数据;dj为第j个校正集样本与待测原油间的光谱距离,n为校正集样本的个数;pj为第j个校正集样本的性质值。
具体的所述的光谱距离dj的计算方法,采用相似样本和待测样本的吸光度距离来代表样本间的光谱距离,计算方法如下:
Figure GDA0002818663120000022
式中,dj指待测样本与第j个相似样本间的光谱距离,j=1,2,...,n,n为校正集样本的个数;Xi为待测样本的第i个波数下的吸光度;1≤i≤k,k为光谱预处理后截取的波数点总数;Xji为第j个相似样本在第i个波数下的吸光度;
Figure GDA0002818663120000031
为光谱距离系数,
Figure GDA0002818663120000032
取常数20000。
优选的相似样本阈值M=10。
有益效果:
本发明提出一种原油精细物性的近红外快速分析方法,基于近红外技术,采用加权距离方法建立原油精细物性的预测模型,可对原油中链烷烃、烯烃、环烷烃和芳烃等烃族组成这类精细物性快速分析,有助于生产中对石油烃类的定向转化,实现对石油资源的最佳高效利用。
附图说明
图1原油精细物性的近红外快速分析的流程图
具体实施案例
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本案例以原油PONA值中的链烷烃(P)为例,对某炼化企业原油库中的原油样本,首先采用近红外光谱仪扫描获取样本光谱数据,然后运用本发明提出的加权距离局部建模技术,对待测原油链烷烃的质量分数进行预测。下面以马林原油(MAL)和沙重原油(SAH)为待测原油样本,给出待测原油样本链烷烃质量分数预测的详细过程。图1为基于近红外光谱的原油精细物性预测的流程图。具体的实施步骤如下:
(1)获取原油样本的PONA值和光谱数据
选取原油库中的典型原油样本41个,样本的烃族组成通过气相色谱法测得,样本的链烷烃、环烷烃和芳烃的质量分数如表1所示,烯烃因含量很少,表1中略去。
表1原油库中样本的PONA值
Figure GDA0002818663120000033
Figure GDA0002818663120000041
(2)光谱数据预处理
根据扫描得到的光谱,选取6400cm-1和9200cm-1两个波数点进行基线校正,并截取波数段4000~4740cm-1内的光谱数据进行归一化处理。选取其中186个波数对应的吸光度,在此波数段范围内建立模型。
(3)采用主成分分析计算得分矩阵
根据186个光谱吸光度组成的矩阵,采用MATLAB自带的函数princomp函数计算其得分矩阵,由于前两个主成分的累积贡献率达88.6%,已经超过85%,所以选取前两个主成分对应的得分向量组成新的得分矩阵,如表2所示。
表2原油样本光谱得分向量
Figure GDA0002818663120000042
待测样本为MAL原油和SAH原油,经过主成分分析后,得到待测样本光谱的前两个主成分,如表3所示:
表3待测原油样本光谱得分向量
Figure GDA0002818663120000043
Figure GDA0002818663120000051
(4)进行相似样本选取
由于选取了前两个主成分,所以建立长宽比例为2:1的二维矩形框,本案例的矩形框的长取为0.4,宽为0.2。基于主成分分析方法得到的光谱得分向量以及设定的矩形框,以待测原油样本的得分向量坐标为矩形框的中心,进行相似样本选择。
(5)建立校正样本集
本案例中相似样本阈值M=10,当相似样本数量M≥10时,该待测样本的不确定度为1;当相似样本数量5<M≤10时,不确定度为2;若相似样本数量M<4时,不确定度为3,不进行建模预测。基于以上规则,两个待测原油样本统计数据如表4所示:
表4待测原油样本不确定度
原油名称 相似样本个数 是否大于阈值M 是否小于阈值0.5M 不确定度
MAL 5 2
SAH 8 2
由表4的统计结论可知,两个待测原油样本的不确定度均为2,可以进行建模预测,将全部选取相似样本组成各自的校正样本集,并根据式(1)计算待测样本与校正集样本的光谱距离,分别如表5和表6所示。
表5MAL原油校正集样本链烷烃性质和光谱距离
序号 原油名称 链烷烃质量百分比(wt%) 与待测样本的光谱距离
1 SRS 4.2109 36.7434
2 DAL 4.1041 22.5460
3 NAP 7.7796 48.2229
4 RCD 11.2666 18.6518
5 STN 12.3450 19.4925
表6SAH原油校正集样本链烷烃性质和光谱距离
Figure GDA0002818663120000052
Figure GDA0002818663120000061
(6)建立链烷烃性质预测模型
基于表5和表6的数据,采用式(2)的加权距离局部建模方法分别计算MAL原油和SAH原油的链烷烃质量百分比,并与实际值比较,结果如表7所示。
表7待测样本预测结果
待测原油名称 实际值(wt%) 预测值(wt%) 偏差(wt%)
MAL 8.2738 8.5960 0.3222
SAH 8.1596 7.1531 1.0064
为与常规的化学计量学方法进行比较,采用偏最小二乘方法PLS建立链烷烃性质预测模型,与基于加权距离的局部建模方法所建立的模型预测效果进行性能比较,结果如表8所示。
表8基于两种方法建模预测结果比较
Figure GDA0002818663120000062
由表8可知,与基于PLS方法相比较,基于加权距离的局部建模方法预测的链烷烃质量百分比误差较小。特别是在样本数较少的情况下,结果更为明显,可用于原油链烷烃质量百分比性质值的预测。可见,基于该方法可以快速分析原油的烃族组成,有助于炼化企业选择合适的加工路线,从而实现对石油烃类的定向转化。

Claims (5)

1.一种原油精细物性的近红外快速分析方法,其特征在于针对原油烃族组成这类精细物性,采用加权距离方法建立近红外预测模型,对待测样本进行快速分析,具体包括以下步骤:
1)获取待测原油样本的近红外光谱数据;
2)对待测原油样本和原油库中样本的光谱数据进行基线校正、波数段截取和矢量归一化;
3)采用主成分分析方法计算吸光度矩阵的得分矩阵,提取累积贡献率达85%对应的得分向量组成新的得分矩阵;
4)基于新的得分矩阵,以待测样本为中心建立边长比例固定的框,进行相似样本选取;
5)统计框内相似样本数量num,设定相似样本阈值M,对待测样本的不确定度进行判断,依据规则选取相似样本建立校正集;
6)基于选取的校正集样本,采用加权距离方法建立原油精细物性预测模型,并对待测样本进行快速分析;
采用加权距离方法计算原油精细物性:
Figure FDA0002818663110000011
式中,p为待测原油的精细物性数据;dj为第j个校正集样本与待测原油间的光谱距离,n为校正集样本的个数;pj为第j个校正集样本的精细物性值;
光谱距离dj的计算方法,采用相似样本和待测样本的吸光度距离来代表样本间的光谱距离,计算方法如下:
Figure FDA0002818663110000012
式中,dj指待测样本与第j个相似样本间的光谱距离,j=1,2,...,n,n为校正集样本的个数;Xi为待测样本的第i个波数下的吸光度;1≤i≤k,k为光谱预处理后截取的波数点总数;Xji为第j个相似样本在第i个波数下的吸光度;
Figure FDA0002818663110000014
为光谱距离系数,
Figure FDA0002818663110000013
取常数20000。
2.根据权利要求1所述的一种原油精细物性的近红外快速分析方法,其特征在于原油的精细物性指原油中所含链烷烃、烯烃、环烷烃和芳烃这类烃族组成的质量分数。
3.根据权利要求1所述的一种原油精细物性的近红外快速分析方法,其特征在于相似样本边界框的建立规则,根据主成分分析的结果选取前n个主成分,建立比例为n:(n-1):……:1的n维框,其中n为步骤3)中满足累积贡献率的主成分数量。
4.根据权利要求1所述的一种原油精细物性的近红外快速分析方法,其特征在于步骤5)中校正集样本的建立规则:
若相似样本数量num≥M,则该待测样本的不确定度为1,依据与待测样本的距离由近及远选取M个相似样本组成校正集;
若相似样本数量0.5M≤num<M,则该待测样本的不确定度为2,选取num个相似样本组成校正集;
若相似样本数量num<0.5M,则该待测样本的不确定度为3,不进行建模预测。
5.根据权利要求4所述的一种原油精细物性的近红外快速分析方法,其特征在于相似样本阈值M=10。
CN201811629590.9A 2018-12-28 2018-12-28 一种原油精细物性的近红外快速分析方法 Active CN109632691B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811629590.9A CN109632691B (zh) 2018-12-28 2018-12-28 一种原油精细物性的近红外快速分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811629590.9A CN109632691B (zh) 2018-12-28 2018-12-28 一种原油精细物性的近红外快速分析方法

Publications (2)

Publication Number Publication Date
CN109632691A CN109632691A (zh) 2019-04-16
CN109632691B true CN109632691B (zh) 2021-03-30

Family

ID=66079113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811629590.9A Active CN109632691B (zh) 2018-12-28 2018-12-28 一种原油精细物性的近红外快速分析方法

Country Status (1)

Country Link
CN (1) CN109632691B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111829976A (zh) * 2019-04-18 2020-10-27 中国石油化工股份有限公司 一种由原油近红外光谱预测其汽油馏分烃族组成的方法
CN112710625B (zh) * 2020-12-16 2022-07-12 南京富岛信息工程有限公司 一种粒料和粉料状固体聚合物性质快速检测方法
CN113433088B (zh) * 2021-06-25 2024-08-02 南京富岛信息工程有限公司 一种原油长输管线混油段精细化监测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388123A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种由近红外光谱预测原油性质的方法
CN105466884A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 一种通过近红外光谱识别原油种类及其性质的方法
CN108009559A (zh) * 2016-11-02 2018-05-08 哈尔滨工业大学 一种基于空谱联合信息的高光谱数据分类方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160195507A1 (en) * 2015-01-05 2016-07-07 Saudi Arabian Oil Company Characterization of crude oil and its fractions by fourier transform infrared spectroscopy (ftir) analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388123A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种由近红外光谱预测原油性质的方法
CN105466884A (zh) * 2014-09-04 2016-04-06 中国石油化工股份有限公司 一种通过近红外光谱识别原油种类及其性质的方法
CN108009559A (zh) * 2016-11-02 2018-05-08 哈尔滨工业大学 一种基于空谱联合信息的高光谱数据分类方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于粒子群算法实现原油加工配方优化;蒋立沫等;《炼油与化工》;20180131;第29卷(第1期);1-5 *

Also Published As

Publication number Publication date
CN109632691A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN109632691B (zh) 一种原油精细物性的近红外快速分析方法
US9846147B2 (en) Prediction of refining characteristics of oil
CA2740845C (en) Estimating detailed compositional information from limited analytical data
US10082807B2 (en) Method of preparing a certified target product from a mixture of components by spectral analysis
CN107976419B (zh) 一种由油品近红外光谱预测其性质的方法
US11905468B2 (en) Assemblies and methods for enhancing control of fluid catalytic cracking (FCC) processes using spectroscopic analyzers
CN113724800A (zh) 一种基于分子级机理模型与大数据技术的催化裂化装置模拟预测方法
CN109324014B (zh) 一种自适应的原油性质近红外快速预测方法
US20150041634A1 (en) Characterization and prediction of jet fuel quality
CN105987886B (zh) 近红外光谱测定加氢裂化尾油性质的方法
CN109507352B (zh) 一种预测石化生产中任意流股分子组成的方法
Cui et al. Computer-aided gasoline compositional model development based on GC-FID analysis
CN111899812A (zh) 一种石油加工装置的产物模拟方法、系统、设备和存储介质
CN109580918B (zh) 一种用于预测石脑油的分子组成的方法
CN110763651B (zh) 一种适合生产润滑油基础油的原油调配方法
Doicin et al. Octane number estimation using neural networks
WO2020139515A1 (en) Method for determining the composition and properties of hydrocarbon fractions by spectroscopy or spectrometry
US9963647B2 (en) Method to optimize crude slate for optimum hydrodesulfurization performance
CN113793652A (zh) 一种基于分段智能优选的谱图化学计量学解析方法
CN111103257B (zh) 一种快速预测蒸汽裂解液相油品中单环芳烃收率的方法
CN113782112A (zh) 一种石油馏分组成模型确定方法及装置
CN111103258A (zh) 近红外光谱预测蒸汽裂解汽油中单环芳烃收率模的方法
Baldrich Ferrer et al. Detailed chemical characterization of petroleum middle fractions by chemometrics analysis of their ultraviolet spectrum
CN109580917B (zh) 一种用于预测石脑油的分子组成的方法
US20220299951A1 (en) Compositional property estimation models relating to processes and related methods

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