首页 / 专利库 / 燃料种类 / 辛烷值 / 一种判断辛烷值模型维护更新的方法

一种判断辛烷值模型维护更新的方法

阅读:487发布:2020-05-14

专利汇可以提供一种判断辛烷值模型维护更新的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种判断 辛烷值 模型维护更新的方法,其特征在于采用 光谱 距离 阈值 作为分界点来判断待测样本是否为模型界外样本,是则进行报警以提醒工作人员需对现有模型进行维护更新。收到报警 信号 后,工作人员对当前待测组分进行标准化验,并添加至模型中,以提高模型的预测 精度 。本发明提出的按需维护更新的方法相对于传统定时维护更新的方法,在保证模型预测精度的前提下,大大节约了时间人 力 资源。,下面是一种判断辛烷值模型维护更新的方法专利的具体信息内容。

1.一种判断辛烷值模型维护更新的方法,其特征在于采用光谱距离阈值作为分界点来判断待测样本是否为模型界外样本,是则进行报警以提醒工作人员需对现有模型进行维护更新;所述判断待测样本是否为模型界外样本是指判断待测样本是否脱离模型的覆盖范围,具体步骤如下:
(1)对待测样本和模型的光谱库中的光谱进行常规预处理;
(2)计算光谱库的光谱重心Score,即光谱库光谱的平均光谱;
(3)采用主成分分析,计算预处理后待测样本光谱、光谱库光谱和光谱库的光谱重心Score的得分向量;
(4)基于得分向量的第一和第二主成分,计算光谱库光谱到光谱库的光谱重心Score的光谱距离,并确定最大光谱距离Dmax;
(5)基于得分向量的第一和第二主成分,利用光谱库的光谱重心Score、最大光谱距离Dmax,以及待测样本确定待测样本的邻近区域;
(6)在邻近区域内查找待测样本的邻近光谱库样本,组成邻近光谱库;
(7)计算邻近区域的光谱重心Snewcore;
(8)基于得分向量的第一和第二主成分,计算邻近光谱库光谱、光谱库的光谱重心Score分别与邻近区域的光谱重心Snewcore的光谱距离,并找出其中最大光谱重心距离Dnewmax,确定光谱距离阈值DST,确定方法如下式:
DST=μ×Dnewmax
其中μ为常数,1.0≤μ≤1.1;
(9)基于得分向量的第一和第二主成分,计算待测样本与邻近区域的光谱重心Snewcore的光谱样本距离DPA;
(10)判断光谱样本距离DPA是否大于光谱距离阈值DST,若大于,则待测样本为模型界外样本,需给出报警,否则不需给出报警。
2.根据权利要求1所述的一种判断辛烷值模型维护更新的方法,其特征在于步骤(5)中所述的邻近区域,确定方法如下:
(501)基于得分向量的第一和第二主成分建立二维坐标系,其中第一主成分为横坐标,第二主成分为纵坐标,以光谱库的光谱重心Score所在的点为端点,连接待测样本所在的点,得到射线;
(502)设定初始度α,以光谱库的光谱重心Score为轴中心,以射线为对称轴,对称轴长为最大光谱距离Dmax,旋转对称轴,旋转角度为±α,得到扇形区域即为所述的邻近区域。
3.根据权利要求2所述的一种判断辛烷值模型维护更新的方法,其特征在于所述初始角度α小于5°。
4.根据权利要求2所述的一种判断辛烷值模型维护更新的方法,其特征在于步骤(6)中在邻近区域内查找待测样本的邻近光谱库样本,具体步骤如下:
(601)扇形区域的两个端点和一个顶点在二维坐标系中的坐标分别为O1(x1,y1)、O2(x2,y2)和O3(x3,y3),光谱库中第i样本的坐标Oi(xi,yi);
(602)以O1、O2和O3三点确定三角形,求出三角形三边的长,然后求三角形面积A;
(603)采用步骤(602)的方法分别计算Oi与O1、O2、O3任意两点所围成的三角形面积Ai12、Ai13、Ai23;
(604)判断A是否等于Ai12、Ai13、Ai23之和,若是,则此光谱库中第i样本为待测样本的邻近样本Sneighbor,否则不是。
5.根据权利要求1所述的一种判断辛烷值模型维护更新的方法,其特征在于步骤(7)中所述的邻近区域的光谱重心Snewcore,采用下式计算:
式中,n为邻近光谱库样本的总数,Snewcore为邻近区域的光谱重心,Score为光谱库的光谱重心,Sneighbor为待测样本的邻近样本。
6.根据权利要求1所述的一种判断辛烷值模型维护更新的方法,其特征在于所述的最大光谱距离Dmax、最大光谱重心距离Dnewmax、光谱样本距离DPA,采用的是欧式距离,计算公式如下:
其中Di,j是指第i个光谱与第j个光谱间的欧式距离;ti为第i个样本光谱的得分向量,tj为第j个样本光谱的得分向量。

说明书全文

一种判断辛烷值模型维护更新的方法

技术领域

[0001] 本发明涉及炼油企业汽油性质检测,尤其是汽油性质预测模型的维护,具体是一种判断辛烷值模型维护更新的方法。

背景技术

[0002] 近红外光谱定量分析技术,又称“黑匣子”分析技术,是一种间接的测量方法,即通过对样品光谱和其质量参数进行关联建立起预测模型,然后通过预测模型和未知质量参数的样品光谱来预测样品的组成和性质。
[0003] 在汽油标准中,辛烷值(RON)列于车用汽油规格的首项,是汽油关键质量指标之一。为了能够及时掌握调合汽油及调合组分油的辛烷值变化,为炼厂调合及生产调控提供有效依据,目前一般采用近红外光谱定量分析技术对其进行预测。但是,一般情况下建立的预测模型涵盖的调合汽油或调合组分油的辛烷值范围有限,不能囊括全部情况。当生产工艺发生变化时,汽油或调合油的辛烷值可能超出预测模型的范围,导致预测不准,不能给生产人员提供正确的指导。

发明内容

[0004] 为解决现有技术存在的问题,本发明提出了一种判断辛烷值模型维护更新的方法,该方法是指当待测样本因工艺改变导致辛烷值发生突变,采用光谱距离阈值作为分界点来判断待测样本是否为模型界外样本,是则进行报警以提醒工作人员需对现有模型进行维护更新;所述判断待测样本是否为模型界外样本是指判断待测样本是否脱离模型的覆盖范围,具体步骤如下:
[0005] (1)对待测样本和模型的光谱库中的光谱进行常规预处理;
[0006] (2)计算光谱库的光谱重心Score,即光谱库光谱的平均光谱;
[0007] (3)采用主成分分析,计算预处理后待测样本光谱、光谱库光谱和光谱库的光谱重心Score的得分向量;
[0008] (4)基于得分向量的第一和第二主成分,计算光谱库光谱到光谱库的光谱重心Score的光谱距离,并确定最大光谱距离Dmax;
[0009] (5)基于得分向量的第一和第二主成分,利用光谱库的光谱重心Score、最大光谱距离Dmax,以及待测样本确定待测样本的邻近区域;
[0010] (6)在邻近区域内查找待测样本的邻近光谱库样本,组成邻近光谱库;
[0011] (7)计算邻近区域的光谱重心Snewcore;
[0012] (8)基于得分向量的第一和第二主成分,计算邻近光谱库光谱、光谱库的光谱重心Score分别与邻近区域的光谱重心Snewcore的光谱距离,并找出其中最大光谱重心距离Dnewmax,确定光谱距离阈值DST,确定方法如下式:
[0013] DST=μ×Dnewmax
[0014] 其中μ为常数,1.0≤μ≤1.1;
[0015] (9)基于得分向量的第一和第二主成分,计算待测样本与邻近区域的光谱重心Snewcore的光谱样本距离DPA;
[0016] (10)判断光谱样本距离DPA是否大于光谱距离阈值DST,若大于,则待测样本为模型界外样本,需给出报警,否则不需给出报警。
[0017] 优选的,步骤(5)中所述的邻近区域,确定方法如下:
[0018] (501)基于得分向量的第一和第二主成分建立二维坐标系,其中第一主成分为横坐标,第二主成分为纵坐标,以光谱库的光谱重心Score所在的点为端点,连接待测样本所在的点,得到射线;
[0019] (502)设定初始度α,以光谱库的光谱重心Score为轴中心,以射线为对称轴,对称轴长为最大光谱距离Dmax,旋转对称轴,旋转角度为±α,得到扇形区域即为所述的邻近区域。
[0020] 更优的,所述初始角度α小于5°。
[0021] 优选的,步骤(6)中在邻近区域内查找待测样本的邻近光谱库样本,具体步骤如下:
[0022] (601)扇形区域的两个端点和一个顶点(即光谱库的光谱重心Score)在二维坐标系中的坐标分别为O1(x1,y1)、O2(x2,y2)和O3(x3,y3),光谱库中第i样本的坐标Oi(xi,yi);
[0023] (602)以O1、O2和O3三点确定三角形,求出三角形三边的长,然后求三角形面积A;
[0024] (603)采用步骤(602)的方法分别计算Oi与O1、O2、O3任意两点所围成的三角形面积Ai12、Ai13、Ai23;
[0025] (604)判断A是否等于Ai12、Ai13、Ai23之和,若是,则此光谱库中第i样本为待测样本的邻近样本Sneighbor,否则不是。
[0026] 优选的,步骤(7)中所述的邻近区域的光谱重心Snewcore,采用下式计算:
[0027]
[0028] 式中,n为邻近光谱库样本的总数,Snewcore为邻近区域的光谱重心,Score为光谱库的光谱重心,Sneighbor为待测样本的邻近样本。
[0029] 优选的,所述的最大光谱距离Dmax、最大光谱重心距离Dnewmax、光谱样本距离DPA,采用的是欧式距离,计算公式如下:
[0030]
[0031] 其中Di,j是指第i个光谱与第j个光谱间的欧式距离;ti为第i个样本光谱的得分向量,tj为第j个样本光谱的得分向量。
[0032] 有益效果:
[0033] 本发明提出了一种判断辛烷值模型维护更新的方法,当待测样本为模型外样本时,能及时提醒工作人员组分油或调合油由于工艺的改变,其辛烷值发生了突变,采用当前模型进行预测可能结果不准。收到报警信号后,工作人员对当前待测组分进行标准化验,并添加至模型中,以提高模型的预测精度。本发明提出的按需维护更新的方法相对于传统定时维护更新的方法,在保证模型预测精度的前提下,大大节约了时间人资源。附图说明
[0034] 图1是判断辛烷值模型维护更新的步骤。
[0035] 图2是S_Zorb-72光谱、光谱库光谱和光谱库重心Score的二维分布图。
[0036] 图3是S_Zorb-72待测样本的邻近区域。
[0037] 图4是S_Zorb-73待测样本的邻近区域。

具体实施方式

[0038] 下面结合附图和实施案例对本发明作进一步的说明。
[0039] 本发明以某炼化企业汽油调合业务的实际实施情况并结合一具体算例,给出详细的计算过程和具体的操作流程。该企业目前参与调合的组分油有S_Zrob汽油、己烷轻石、# # # #重整抽余油、MTBE等,调合油有国V92、95和国IV93 、97,为给调合生产提供有效调控依据,该企业目前采用近红外光谱定量分析技术对组分油及调合汽油的性质进行预测。但由于生产装置不稳定,导致S_Zrob汽油、己烷轻石、重整抽余油、MTBE等组分油的辛烷值会产生较大波动,导致待测样本性质预测不准。
[0040] 本专利以S_Zrob汽油为例,介绍当待测样本超出辛烷值模型预测范围时给出警报的过程。S_Zorb辛烷值预测模型建模样本共71个,预处理方法为基线校正和矢量归一-1 -1 -1化,采用有效波数段范围为7800-9000cm ,5100-6030cm ,6400-6480cm ,模型光谱库中样本编号及辛烷值如表1所示,即S_Zorb-1~S_Zorb-71,其中编号为S_Zorb-72和S_Zorb-73为两个待测样本。
[0041] 表1 S_Zorb样本编号及对应辛烷值
[0042]
[0043] 首先判断待测样本S_Zorb-72是否为模型界外样本,将S_Zorb-72和模型光谱库中的光谱S_Zorb-1~S_Zorb-71进行基线校正和矢量归一化,然后计算光谱库的光谱重心Score,即光谱库光谱的平均光谱。采用主成分分析,计算预处理后待测样本光谱、光谱库光谱和光谱库重心Score的第一和第二主成分,如表2所示:
[0044] 表2待测样本S_Zorb-72光谱、光谱库光谱和光谱库重心Score的第一、第二主成分[0045]
[0046]
[0047] 优选的实施例中,基于表1数据,以得分向量的第一主成分为横坐标,第二主成分为纵坐标,将待测样本光谱、光谱库光谱和光谱库重心Score绘在二维图上展示,如图2所示。
[0048] 优选的实施例中,基于得分向量的第一、第二主成分,采用欧式距离公式计算光谱库的光谱重心Score到所有光谱库光谱的距离,查找最大光谱距离为Dmax=0.1481。以光谱库的光谱重心Score所在的点为端点,连接待测样本所在的点,得到射线。以光谱库的光谱重心Score为顶点,以射线为对称轴,对称轴长为最大光谱距离Dmax,旋转对称轴,优选的实施例中,旋转角度为±5°,得到扇形区域即S_Zorb-72的邻近区域如图3所示。扇形的两个端点坐标分别为[-0.1127-0.0954]和[-0.1276-0.0743],扇形顶点即光谱库的光谱重心Score。
[0049] 在扇形区域内查找待测样本S_Zorb-72的邻近光谱库样本为S_Zorb-3、S_Zorb-17、S_Zorb-24、S_Zorb-58、S_Zorb-59共计5个,采用下式
[0050]
[0051] 计算扇形区域即邻近光谱库的重心Snewcore,具体数据如表3所示:
[0052] 表3 S_Zorb-72待测样本的邻近光谱库样本
[0053]编号 第一主成分 第二主成分
S_Zorb-3 -0.0406 -0.0267
S_Zorb-17 -0.0234 -0.0139
S_Zorb-24 -0.0418 -0.0284
S_Zorb-58 -0.0446 -0.0299
S_Zorb-59 -0.0426 -0.0261
S_Zorb-72 -0.0296 -0.0210
Score 0.0004 0.0003
Snewcore -0.0191 -0.0124
[0054] 计算邻近区域的光谱重心(邻近光谱库重心)Snewcore到所有邻近光谱库光谱及光谱库的光谱重心Score的距离,查找最大光谱距离为Dnewmax=0.0310,将μ设为1.1,采用下式计算光谱距离阈值DST:
[0055] DST=1.1*0.0310=0.0341
[0056] 优选的实施例中,采用欧氏距离公式计算待测样本S_Zorb-72到邻近区域的光谱重心(邻近光谱库重心)Snewcore的距离DPA=0.0136,很明显DPA小于DST,因此S_Zorb-72,是模型界内样本,不需要给出警报。
[0057] 采用上述方法,判断待测样本S_Zorb-73是否为模型界外样本,表4为经过基线校正和矢量归一化后待测样本S_Zorb-73光谱、光谱库光谱和光谱库的光谱重心Score得分向量的第一、第二主成分,基于第一、第二主成分计算光谱库光谱到光谱库的光谱重心Score的最大光谱距离Dmax为0.1481。
[0058] 表4待测样本S_Zorb-73光谱、光谱库光谱和光谱库重心Score的第一、第二主成分[0059]
[0060] 图4为S_Zorb-73待测样本的邻近区域,在邻近区域内S_Zorb-73的邻近光谱只有一个,即S_Zorb-16,坐标为(0.0299-0.0260)。通过S_Zorb-16和光谱库的光谱重心Score计算得到邻近区域的光谱重心(邻近区域光谱库重心)Snewcore,其坐标为(0.0147-0.0131)。计算邻近区域内的最大光谱重心距离Dnewmax,即S_Zorb-16到Snewcore的光谱距离为0.0200,优选的实施例中,将μ设为1.1则光谱距离阈值DST为0.0220。待测样本S_Zorb-73到Snewcore的光谱距离DPA=0.1104,明显大于DST,因此S_Zorb-73为模型界外样本,需要给出报警。
[0061] 尽管本发明已经参照附图和优选实例进行了说明,但是,对于本领域的技术人员来说,本发明可以有各种更改和变化。本发明的各种更改、变化,和等同物有所附的权利要求书的内容涵盖。
[0062] 本发明未涉及技术均与现有技术相同,或可采用现有技术实现。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈