CN115860211A - 一种基于局部在线建模的铸坯质量预测方法 - Google Patents
一种基于局部在线建模的铸坯质量预测方法 Download PDFInfo
- Publication number
- CN115860211A CN115860211A CN202211513791.9A CN202211513791A CN115860211A CN 115860211 A CN115860211 A CN 115860211A CN 202211513791 A CN202211513791 A CN 202211513791A CN 115860211 A CN115860211 A CN 115860211A
- Authority
- CN
- China
- Prior art keywords
- data set
- data
- casting blank
- local
- sample
- 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
- 238000005266 casting Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000009749 continuous casting Methods 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000004519 manufacturing process Methods 0.000 claims description 38
- 238000012545 processing Methods 0.000 claims description 20
- 238000010606 normalization Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 11
- 238000007637 random forest analysis Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 7
- 238000002790 cross-validation Methods 0.000 claims description 5
- 229910000831 Steel Inorganic materials 0.000 abstract description 14
- 239000010959 steel Substances 0.000 abstract description 14
- 238000003723 Smelting Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 10
- 238000001816 cooling Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 238000005096 rolling process Methods 0.000 description 5
- 229910000677 High-carbon steel Inorganic materials 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000011112 process operation Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005021 gait Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 108010074864 Factor XI Proteins 0.000 description 1
- 101100533306 Mus musculus Setx gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003749 cleanliness Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005204 segregation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于局部在线建模的铸坯质量预测方法,属于钢铁冶炼技术领域,能够解决传统离线全局模型预测铸坯质量时自适应能力差、预测精度不稳定的问题;该方法以支持向量回归模型为局部模型进行局部在线建模,实现连铸过程中对铸坯质量的预测;局部在线建模的样本数据集为:采用相似度函数选择出的与待测数据最相关的历史数据集;在相似度计算时,在传统欧氏距离的基础上引入工艺参数特征权重和样本时间特征权重来强化样本数据与待测数据的相关性,并根据样本相似度从历史数据集中选出最相关数据构建成样本数据集用于局部建模。
Description
技术领域
本发明涉及钢铁冶炼技术领域,尤其涉及一种基于局部在线建模的铸坯质量预测方法。
背景技术
近年来,我国钢铁行业正处于产业结构调整和优化的关键时期。在日趋激烈的市场竞争中,不断提升产品质量是钢铁企业成功发展的必经之路。连铸工序是钢铁生产流程中承上启下的重要环节,而铸坯质量直接影响最终钢铁产品质量。在连铸生产过程中,钢水洁净度较差、工艺操作不稳定、设备状况不佳等因素往往造成不同程度的铸坯质量缺陷。常见的铸坯质量问题主要包括表面裂纹、夹渣、气孔等表面缺陷和中心裂纹、偏析、夹杂、疏松、缩孔等内部缺陷。若能基于连铸实时状况对铸坯质量进行在线预测,就能有效减少或避免严重缺陷铸坯的发生率,保障后续轧制工序的顺利进行,进而达到最终产品质量要求。
早期的铸坯质量预测主要基于连铸生产机理来构建数学模型。不过,连铸生产过程的钢水成分、工艺操作、设备状况等影响因素与铸坯质量缺陷之间具有强耦合、非线性及不确定性关系,不少机理研究并未通透,导致数学机理模型应用范围十分有限。随着数字化、信息化技术的蓬勃发展,以数据挖掘和神经网络为代表的智能化技术应运而生,逐步走向工业应用阶段,为解决铸坯质量在线预测提供新的思路。
现以应用与钢铁领域连铸坯质量预测中的方法有BP神经网络、支持向量机、极限学习机等,此类方法的共同特点都是先在离线状态下对历史数据进行训练得到全局模型,再基于离线模型结合待测样本输出预测结果。然而,连铸生产是具有时变性、非线性、强耦合特征的复杂过程,时间推移、设备更新、工艺改进等客观因素往往导致模型特征与实际工况不相匹配。而上述的铸坯质量预测模型属于离线全局模型,模型建立后不会随着生产过程的改变而自适应调整,无法跟踪生产状态的变化,极易导致模型性能下降甚至失效。因此,亟须开发一种自适应能力更好、预测精度更稳定的铸坯质量预测方法。
即时学习算法是一种具有较高自适应能力的人工智能技术。与传统全局学习算法不同,即时学习算法采用“边建模、边预测”的运行方式,通过相似性计算选出与待测样本最相关的历史样本在线构建局部模型,实现长时间的高精度稳定预测。相似样本的选择或者样本权重的计算是即时学习算法的关键步骤,在很大程度上影响模型的预测精度。传统的样本相似性度量函数:欧氏距离、马氏距离、曼哈顿距离等,而基于上述度量函数选择的相似样本仅考虑了查询样本与历史样本输入变量的相似度,并未考虑变量之间的相关性以及生产过程的时变性,可能导致模型的预测精度不稳定。
因此,需要研究一种基于局部在线建模的铸坯质量预测方法来应对现有技术的不足,以解决或减轻上述一个或多个问题。
发明内容
有鉴于此,本发明提供了一种基于局部在线建模的铸坯质量预测方法,能够解决传统离线全局模型预测铸坯质量时自适应能力差、预测精度不稳定的问题。
本发明提供一种基于局部在线建模的铸坯质量预测方法,所述预测方法以支持向量回归模型为局部模型进行局部在线建模,实现在连铸过程中对铸坯质量的预测;局部在线建模的样本数据集为:采用相似度函数选择出的与待测数据最相关的历史数据集;
在相似度计算时,在传统欧氏距离的基础上引入工艺参数特征权重、样本时间特征权重来强化样本数据与待测数据的相关性,并根据样本相似度从历史数据集中选出最相关数据构建成样本数据集用于局部建模。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述预测方法的具体步骤包括:
S1、获取连铸生产过程中与铸坯质量相关的工况数据以及对应的质量数据,作为初始历史数据集,
S2、对所述初始历史数据集进行预处理得到标准样本数据集;
S3、根据所述标准样本数据集建立离线随机森林模型,并通过交叉验证的方式计算得到各工艺参数的特征权重;
S4、对待测数据进行预处理,根据S3得到的各工艺参数的特征权重计算待测数据与所述标准样本数据集中所有样本的加权欧氏距离;
S5、在所述标准样本数据集中引入时间特征权重,根据所述时间特征权重和所述加权欧氏距离计算得到待测数据与所述标准样本数据集中所有样本的相似度权值;
S6、根据得到的相似度权值从所述标准样本数据集中选择出满足要求的样本数据,构建成相似样本数据集;
S7、以得到的所述相似样本数据集为训练数据构建局部模型;
S8、将待测数据输入构建好的局部模型中进行预测,得到预测结果。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S2中的预处理包括:空值处理、异常值处理和归一化处理。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S1中质量数据的获得方式为通过低倍检测采集得到。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S4中的预处理包括归一化处理。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S6的具体内容包括:对相似度权值进行排序,选择相似度权值最大的前Nk个样本数据,构建成相似样本数据集;Nk为大于1的正整数。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S8得到预测结果后,丢弃当前的局部模型。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述预测方法的步骤还包括:
S9、得到待测数据对应的真实样品,并通过低倍检测方式得到真实质量数据,将该待测数据和其真实质量数据添加到初始历史数据集中。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,针对下一待测数据,重复进行步骤S2-S8。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,步骤S1中通过连铸生产过程跟踪系统实时采集工况数据;
所述连铸生产过程跟踪系统能够实现铸坯生产过程数据与铸坯时间方向的精确匹配。
与现有技术相比,上述技术方案中的一个技术方案具有如下优点或有益效果:本发明在样本相似性度量策略下,通过在线局部建模方法,实时查询历史数据库中的相似样本,选取最优的相似样本进行建模,能够自适应地根据连铸生产过程非线性、时变性的工况为系统在线建立当前工作点的局部预测模型,避免了离线全局模型的局限性,有效解决了预测模型中模型参数的在线更新问题;
上述技术方案中的另一个技术方案具有如下优点或有益效果:本发明根据连铸生产过程跟踪系统实时采的生产过程的工况数据构建历史数据库,并采取滚动训练方式基于历史数据库进行建模,这样既能够重复利用有用的历史样本数据,又能不断根据最新的工况数据即时选取最优的相似样本;
上述技术方案中的另一个技术方案具有如下优点或有益效果:本发明在每一个待测样本预测时刻建立新的局部模型,不保留旧模型的参数,从而过去时刻出现的干扰不会影响当前时刻局部预测模型的预测准确性,提高预测模型的鲁棒性、自适应能力和稳定性。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有技术效果。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明一个实施例提供的基于局部在线建模的铸坯质量预测方法的架构示意图;
图2是本发明一个实施例提供的原始数据集划分示意图;
图3是本发明一个实施例提供的在随机森林算法3次交叉验证下各输入变量的特征权重图;
图4是本发明一个实施例提供的在线局部模型和离线支持向量回归模型对65号高碳钢铸坯三角区裂纹等级预测的结果对比图;
图5是本发明一个实施例提供的折算后在线局部模型和离线支持向量回归模型对65号高碳钢铸坯三角区裂纹等级预测的结果对比图。
具体实施方式
为了更好的理解本发明的技术方案,下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明针对传统铸坯质量预测大多基于全局建模方法,存在模型结构和参数相对固定、自适应能力较差、难以在场景复杂、工况多变的连铸过程得到长期稳定应用的问题,提供一种基于局部在线建模的铸坯质量预测方法,以支持向量回归为局部模型,利用相似度函数选择与待测样本最相关的历史数据作为局部建模的子样本集进行局部在线建模,实现在连铸过程中铸坯质量的预测。在相似度计算中基于传统欧氏距离引入工艺参数特征权重、样本时间特征权重来强化样本数据与待测数据的相关性,并根据样本相似度构建相似样本集进行局部建模,在各待测样本预测过程中,采取滚动训练的方式来模拟连铸真实生产过程。此建模方法能够表征当前系统状态,也能很好地适应连铸生产过程数据的时变性、非线性性特点,行之有效地保证了在连铸过程中对铸坯质量的实时精准预测。该预测方法详细的步骤包括:
步骤1:通过连铸生产过程跟踪系统实时采集与储存连铸生产过程中与铸坯质量相关的工况数据集,表示为:X=[x1,x2,...,xn]T,xi∈Rm,i=1,2,…,n,n为样本数量,m为样本维度,R为实数集;通过实验室低倍检测采集的数据,得到各样本对应的真实质量等级值Y=[y1,y2,...,yn]T∈Rn。
步骤2:将采集到的数据作为初始历史数据集Dorg=[X;YT]∈R(m+1)×n,对初始历史数据集Dorg进行数据预处理,得到预处理的数据集DL;
预处理方法包括:空值处理、异常值处理和归一化处理。
步骤3:根据标准化后的初始历史数据集DL建立离线随机森林模型,并通过交叉验证的方式计算得到各输入变量的权重Wvar=[w1,w2,...,wm]T∈Rm。
步骤4:对于新采集的查询样本xnew∈Rm,按照公式(3)进行归一化处理。
步骤5:根据输入变量权重矩阵Wvar计算查询样本xnew与数据集XL中所有样本的加权欧氏距离si。
步骤6:针对连铸生产时变性的特点,在数据集的样本数据中引入时间特征权重wt,并将样本时间特征权重wt和步骤5中的查询样本与数据集XL中所有样本的加权欧氏距离进行融合,计算查询样本与数据集XL中所有样本的相似度权值
步骤7:根据相似度计算结果,在初始历史数据集DL中选择相似度权值最大的前Nk个样本构建相似样本数据集DR。
步骤8:根据相似样本数据集DR建立局部模型;局部建模结构采用支持向量回归模型;
步骤9:查询样本xnew输入局部模型,模型输出预测值f(x)new,并丢弃此局部模型;
步骤10:当查询样本xnew所对应的真实输出值ynew由低倍检测得到时,采取滚动训练的方式将样本[xnew,ynew]添加到初始历史数据集DL中,并更新数据库;否则,维持初始历史数据集DL中所包含的样本不变。
步骤11:下一个待测样本xnew+1进入预测模型重复步骤2-10。
实施例1:
该实施例提供的基于局部在线建模的铸坯质量预测方法,其架构如图1所示。预测方法的具体步骤包括:
步骤1:通过连铸生产过程跟踪系统实时采集与储存连铸生产过程中与铸坯质量相关的工况数据集X=[x1,x2,...,xn]T,xi∈Rm,i=1,2,…,n,n为样本数量,m为样本维度;通过实验室低倍检测采集的数据,得到各样本对应的真实质量等级值Y=[y1,y2,...,yn]T∈Rn。
工况数据获取首先是铸坯生产过程数据与铸坯时间方向的精确匹配。为此,基于铸坯生产过程跟踪系统,计算铸坯任意位置在各个装置(结晶器、二冷各区等)的关键时间节点。根据铸坯头、尾在各个装置的时刻,截取这段时间的生产过程数据,与铸坯时间方向进行匹配。然后再通过时空变换,将铸坯时间方向的数据变换到铸坯长度方向,实现铸坯生产过程数据与铸坯长度方向的精确匹配。
铸坯在结晶器的工况数据与时间方向进行匹配:
铸坯在二冷各区的工况数据与时间方向进行匹配:
式中:lJ-C为结晶器到切割机的距离;lCn-C为二冷区各个冷却区到切割机的距离;tC为铸坯切割时刻;tJ为铸坯在结晶器浇注时刻;tCn为铸坯在二冷区各个冷却区的冷却时刻;v为拉速。
步骤2:将采集到的数据作为初始历史数据集Dorg=[X;YT]∈R(m+1)×n,对初始历史数据集Dorg进行数据预处理,预处理方法包括:空值处理、异常值处理,归一化处理。
空值处理:空值即缺失值,采用均值替换缺失值;
异常值处理:对采集的工艺参数值设置边界范围,当所采集的数据值超出边界范围时,视为异常信号,将其直接删除,不予采用;
Min-Max标准化公式为:
式中,函数min{Dorg}、max{Dorg}分别表示计算矩阵各行的最小值和最大值。
步骤3:根据标准化后的初始历史数据集DL建立离线随机森林模型,并通过交叉验证的方式计算得到各输入变量的权重Wvar=[w1,w2,...,wm]T∈Rm。
随机森林是以K个决策树{h(X,Θk),k=1,2,...,K}为基本分类单元,进行集成学习后得到的一个组合分类器。由于构建每个决策树时,随机抽取训练样本和属性子集的过程是独立的,且总体是一样的,因此参数集{Θk,k=1,2,...,K}是一个独立同分布的随机向量。其中,随机森林求变量权重评分的步骤:
初始历史数据集DL有n个查询的样本,表示各个步态参数的变量依次为F1,F2,…,Fm。应用自助法重采技术有放回地抽取K个新的自助样本集,在此过程中得到k个分类回归树,每次未被抽中的样本则组成了k个袋外数据(Out-of-bag,OOB),该部分数据样本作为测试样本用于评估各个步态参数在分类中的重要性,具体如下:
①用自助样本可以得到若干个树形分类器,同时对应的OOB进行分类,得到k个OOB样本中的每一个样本的投票数,记作v1,v2,…,vk。
②将变量Fi的数值在k个OOB中的顺序随机打乱,形成新的OOB测试集,根据样本判别的正确数,所得到的结果可表示为:
③v1,v2,…,vk与矩阵(3)对应的第i行向量相减,求和平均后得到各变量Fi的权重评分,即:
步骤4:对于新采集的查询样本xnew∈Rm,按照公式(3)进行归一化处理。
步骤5:根据输入变量权重矩阵Wvar计算查询样本xnew与数据集XL中所有样本的加权欧氏距离:
式中,si为查询样本与数据集XL中样本的加权欧氏距离,wi为权重因子,H为查询样本的特征向量,hi为历史样本数据集中样本i的特征向量。
步骤6:针对连铸生产时变性的特点,在样本数据中引入时间特征权重wt,确保分配最相关、最合适的样本进行在线局部建模。并将样本时间特征权重wt和步骤5中的查询样本与数据集XL中所有样本的加权欧氏距离进行融合,计算查询样本与数据集XL中所有样本的相似度权值
定义1时间权向量的熵I。时间权向量的熵是一个度量值,熵值越大,则说明数据集中各样本时间权重之间的差异越小;反之亦然。
定义2时间度λ。时间度λ∈[0,1],λ的大小体现对时序数据的重视程度(见表1),即当λ越小时,表明越注重离当前工况点时间接近的数据;反之亦然。
确定数据集XL中各样本时间权重的准则:在给定时间度λ的值下,以尽可能地挖掘样本的信息和兼顾被评价对象在时许上的差异信息为标准来寻找适合该样本集的时间权向量。此方法可归纳为求解非线性规划问题:
表1时间度λ的标度参考表
步骤7:根据相似度计算结果,在初始历史数据集DL中选择相似度权值最大的前Nk个样本构建相似样本数据集DR。
步骤8:根据相似样本数据集DR建立局部模型,局部建模结构采用支持向量回归模型。
支持向量回归基于支持向量机发展而来,广泛应用于非线性问题的建模和预测。该方法通过引入最大误差值ε,规定预测结果误差的绝对值大于ε的样本参与损失值计算。支持向量回归问题的目标函数为:
式中,样本为(x,y),f(x)为模型输出,y为期望输出,η、b为待定模型参数,ε为偏差最大允许值,C为正则化常数;p为相似样本数据集DR样本数量;lε(f(x)-y)为ε-不敏感损失函数;
进而,通过引入拉格朗日乘子,由拉格朗日乘子法得到拉格朗日函数L,在满足KKT最优化条件下,将回归模型带入拉格朗日函数L,并令其主要变量偏导数的值为零,则支持向量回归可表示为:
步骤9:查询样本xnew输入局部模型,模型输出预测值f(x)new,并丢弃此局部模型。
步骤10:当查询样本xnew所对应的真实输出值ynew由低倍检测得到时,采取滚动训练的方式将样本[xnew,ynew]添加到初始历史数据集DL中,并更新数据库;否则,维持初始历史数据集DL中所包含的样本不变。
步骤11:下一个待测样本xnew+1进入预测模型重复步骤2-10。
实施例2:
以下结合国内某钢厂65号高碳钢铸坯三角区裂纹为例,验证基于即时学习算法的铸坯质量预测方法的有效性。在钢铁生产流程中,铸坯质量直接影响最终钢铁产品质量,三角区裂纹是铸坯典型的内部质量缺陷,主要分布在板坯的三角区部位,裂纹的方向垂直于板坯的窄面,其形成原因与连铸过程的钢水成分、工艺参数、工艺操作、设备状况等因素密切相关。
结合现场调研和前人研究,选取影响铸坯三角区裂纹形成的7个关键因素作为输入变量,具体见表2。同时,以三角区裂纹等级为预测模型输出值,根据裂纹的严重程度可分为0、0.5、1、1.5、2共5种等级。
表2铸坯质量缺陷影响因素
采集2022年1月-6月国内某钢厂65号高碳钢铸坯三角区裂纹的90组实际生产数据作为原始数据集。为贴近实际生产状况,原始数据集划分如图2所示,根据生产时间序列将原始数据集划分为两部分,前70组时间较早的数据作为历史数据集,后20组时间较近的数据用于模拟待测样本。在线局部建模时,最初前70组数据作为初始历史数据集,当各待测样本完成预测后,依次添加到初始历史数据集中,直至20个待测样本全部完成预测。在近邻样本选择时,用于在线局部建模的样本容量Nk设置为30。
以铸坯三角区裂纹为预测对象,在随机森林3次交叉验证下定量分析表1中7个关键输入变量的权重,分析结果如图3所示。根据输入变量权重矩阵Wvar计算查询样本xnew与数据集XL中所有样本的加权欧氏距离;同时根据样本时间特征权重计算方式,通过多次仿真试验对比确定相关参数:时间度λ=0.3,可调参数r=0.4,求解数据集XL中各样本的时间权重wt,并将样本时间特征权重wt和查询样本与数据集XL中所有样本的加权欧氏距离进行融合,计算查询样本与数据集XL中所有样本的相似度权值在初始历史数据集中选取相似度权值最大的前Nk=30个样本进行局部在线建模,局部模型结构采用支持向量回归模型。
表3为使用不同的建模方法对铸坯三角区裂纹等级的预测结果比较,即离线支持向量回归模型与本发明的在线局部建模方法的比较。本发明使用平均绝对误差MAE、均方误差MSE、决定系数R2、准确率A等多种误差度量方式作为综合评价指标。
表3两种建模方式预测的评价结果对比
在实际铸坯质量检测中,铸坯三角区裂纹的评级方法采用0、0.5、1、1.5、2级阶梯型评价,而回归模型得到的结果为连续的数值,结果如图4所示。为与实际铸坯质量检测结果相匹配,回归模型的预测结果采用表4所示的对应关系进行折算。在此基础上,图5所示为预测值按实际评价方法取值后的两种模型对比情况。不难发现,在线局部建模的预测值与真实值重合的个数明显多于离线支持向量回归模型。具体来说,局部模型的预测准确率为0.9,而离线支持向量回归模型的预测准确率仅为0.65,这主要得益于本发明的在线局部建模方式充分考虑了连铸生产过程非线性、时变性的特点,即将变量之间的相关性以及生产过程的时变性融入样本相似度的计算,强化局部建模的样本数据与待测样本的相关性。
表4预测结果与实际评级之间对应关系
以上对本申请实施例所提供的一种基于局部在线建模的铸坯质量预测方法,进行了详细介绍。以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的商品或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种商品或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的商品或者系统中还存在另外的相同要素。“大致”是指在可接收的误差范围内,本领域技术人员能够在一定误差范围内解决所述技术问题,基本达到所述技术效果。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。
Claims (10)
1.一种基于局部在线建模的铸坯质量预测方法,其特征在于,所述预测方法以支持向量回归模型为局部模型进行局部在线建模,实现在连铸过程中对铸坯质量的预测;局部在线建模的样本数据集为:采用相似度函数选择出的与待测数据最相关的历史数据集;
在相似度计算时,在传统欧氏距离的基础上引入工艺参数特征权重和样本时间特征权重来强化样本数据与待测数据的相关性,并根据样本相似度从历史数据集中选出最相关数据构建成样本数据集用于局部建模。
2.根据权利要求1所述的基于局部在线建模的铸坯质量预测方法,其特征在于,所述预测方法的具体步骤包括:
S1、获取连铸生产过程中与铸坯质量相关的工况数据以及对应的质量数据,作为初始历史数据集;
S2、对所述初始历史数据集进行预处理得到标准样本数据集;
S3、根据所述标准样本数据集建立离线随机森林模型,并通过交叉验证的方式计算得到各工艺参数的特征权重;
S4、对待测数据进行预处理,根据S3得到的各工艺参数的特征权重计算待测数据与所述标准样本数据集中所有样本的加权欧氏距离;
S5、在所述标准样本数据集中引入时间特征权重,根据所述时间特征权重和所述加权欧氏距离计算得到待测数据与所述标准样本数据集中所有样本的相似度权值;
S6、根据得到的相似度权值从所述标准样本数据集中选择出满足要求的样本数据,构建成相似样本数据集;
S7、以得到的所述相似样本数据集为训练数据构建局部模型;
S8、将待测数据输入构建好的局部模型中进行预测,得到预测结果。
3.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S2中的预处理包括:空值处理、异常值处理和归一化处理。
4.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S1中质量数据的获得方式为通过低倍检测采集得到。
5.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S4中的预处理包括归一化处理。
6.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S6的具体内容包括:对相似度权值进行排序,选择相似度权值最大的前Nk个样本数据,构建成相似样本数据集;Nk为大于1的正整数。
7.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S8得到预测结果后,丢弃当前的局部模型。
8.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,还包括:
S9、得到待测数据对应的真实样品,并通过低倍检测方式得到真实质量数据,将该待测数据和其真实质量数据添加到初始历史数据集中。
9.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,针对下一待测数据,重复进行步骤S2-S8。
10.根据权利要求2所述的基于局部在线建模的铸坯质量预测方法,其特征在于,步骤S1中通过连铸生产过程跟踪系统实时采集工况数据;
所述连铸生产过程跟踪系统能够实现铸坯生产过程数据与铸坯时间方向的精确匹配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211513791.9A CN115860211B (zh) | 2022-11-29 | 一种基于局部在线建模的铸坯质量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211513791.9A CN115860211B (zh) | 2022-11-29 | 一种基于局部在线建模的铸坯质量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115860211A true CN115860211A (zh) | 2023-03-28 |
CN115860211B CN115860211B (zh) | 2024-10-25 |
Family
ID=
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116258087A (zh) * | 2023-05-15 | 2023-06-13 | 矿冶科技集团有限公司 | 冰铜品位软测量方法、装置、电子设备及存储介质 |
CN116307938A (zh) * | 2023-05-17 | 2023-06-23 | 成都瑞雪丰泰精密电子股份有限公司 | 一种加工中心进给系统健康状态评估方法 |
CN116449788A (zh) * | 2023-05-30 | 2023-07-18 | 湖南科技大学 | 一种钢厂铸坯生产过程故障检测与产能优化的方法 |
CN116823054A (zh) * | 2023-06-30 | 2023-09-29 | 浙大城市学院 | 一种基于数据挖掘的铝合金铸件制造质量实时预测方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107290965A (zh) * | 2017-08-01 | 2017-10-24 | 浙江大学 | 基于局部加权贝叶斯网络的自适应软测量预测方法 |
CN108827905A (zh) * | 2018-04-08 | 2018-11-16 | 江南大学 | 一种基于局部加权Lasso的近红外模型在线更新方法 |
CN109840362A (zh) * | 2019-01-16 | 2019-06-04 | 昆明理工大学 | 一种基于多目标优化的集成即时学习工业过程软测量建模方法 |
JP2019087152A (ja) * | 2017-11-09 | 2019-06-06 | 新日鐵住金株式会社 | 製造プロセスの状態予測装置、方法及びプログラム、並びに製造プロセスの制御システム |
US10678196B1 (en) * | 2020-01-27 | 2020-06-09 | King Abdulaziz University | Soft sensing of a nonlinear and multimode processes based on semi-supervised weighted Gaussian regression |
CN111985723A (zh) * | 2020-08-28 | 2020-11-24 | 浙江工业大学 | 基于即时最小二乘支持向量回归离心泵外特性的预测方法 |
CN112418538A (zh) * | 2020-11-30 | 2021-02-26 | 武汉科技大学 | 一种基于随机森林分类的连铸坯夹杂预测方法 |
CN113012766A (zh) * | 2021-04-27 | 2021-06-22 | 昆明理工大学 | 一种基于在线选择性集成的自适应软测量建模方法 |
CN114117919A (zh) * | 2021-11-29 | 2022-03-01 | 中国石油大学(华东) | 基于样本协同表示的即时学习的软测量建模方法 |
CN114239397A (zh) * | 2021-12-16 | 2022-03-25 | 浙江大学 | 基于动态特征提取与局部加权深度学习的软测量建模方法 |
CN114528764A (zh) * | 2022-02-18 | 2022-05-24 | 清华大学 | 基于整体优化的即时学习的软测量建模方法及装置 |
CN114841073A (zh) * | 2022-05-17 | 2022-08-02 | 中国石油大学(华东) | 基于局部标签传播的即时学习半监督软测量建模方法 |
CN115114968A (zh) * | 2020-10-09 | 2022-09-27 | 武汉科技大学 | 基于随机森林的铸坯裂纹智能预测方法 |
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107290965A (zh) * | 2017-08-01 | 2017-10-24 | 浙江大学 | 基于局部加权贝叶斯网络的自适应软测量预测方法 |
JP2019087152A (ja) * | 2017-11-09 | 2019-06-06 | 新日鐵住金株式会社 | 製造プロセスの状態予測装置、方法及びプログラム、並びに製造プロセスの制御システム |
CN108827905A (zh) * | 2018-04-08 | 2018-11-16 | 江南大学 | 一种基于局部加权Lasso的近红外模型在线更新方法 |
CN109840362A (zh) * | 2019-01-16 | 2019-06-04 | 昆明理工大学 | 一种基于多目标优化的集成即时学习工业过程软测量建模方法 |
US10678196B1 (en) * | 2020-01-27 | 2020-06-09 | King Abdulaziz University | Soft sensing of a nonlinear and multimode processes based on semi-supervised weighted Gaussian regression |
CN111985723A (zh) * | 2020-08-28 | 2020-11-24 | 浙江工业大学 | 基于即时最小二乘支持向量回归离心泵外特性的预测方法 |
CN115114968A (zh) * | 2020-10-09 | 2022-09-27 | 武汉科技大学 | 基于随机森林的铸坯裂纹智能预测方法 |
CN112418538A (zh) * | 2020-11-30 | 2021-02-26 | 武汉科技大学 | 一种基于随机森林分类的连铸坯夹杂预测方法 |
CN113012766A (zh) * | 2021-04-27 | 2021-06-22 | 昆明理工大学 | 一种基于在线选择性集成的自适应软测量建模方法 |
CN114117919A (zh) * | 2021-11-29 | 2022-03-01 | 中国石油大学(华东) | 基于样本协同表示的即时学习的软测量建模方法 |
CN114239397A (zh) * | 2021-12-16 | 2022-03-25 | 浙江大学 | 基于动态特征提取与局部加权深度学习的软测量建模方法 |
CN114528764A (zh) * | 2022-02-18 | 2022-05-24 | 清华大学 | 基于整体优化的即时学习的软测量建模方法及装置 |
CN114841073A (zh) * | 2022-05-17 | 2022-08-02 | 中国石油大学(华东) | 基于局部标签传播的即时学习半监督软测量建模方法 |
Non-Patent Citations (3)
Title |
---|
任明仑等: ""灰铸铁抗拉强度预测的局部加权线性回归建模"", 《电子测量与仪器学报》, vol. 33, no. 3, 15 March 2019 (2019-03-15), pages 65 - 71 * |
彭开香等: ""一种基于即时学习的局部模型集成转炉炼钢终点预测方法"", 《冶金自动化》, vol. 46, no. 2, 31 March 2022 (2022-03-31), pages 110 - 117 * |
陈恒志等: ""基于极限学习机(ELM)的连铸坯质量预测"", 《工程科学学报》, vol. 40, no. 7, 15 July 2018 (2018-07-15), pages 815 - 821 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116258087A (zh) * | 2023-05-15 | 2023-06-13 | 矿冶科技集团有限公司 | 冰铜品位软测量方法、装置、电子设备及存储介质 |
CN116258087B (zh) * | 2023-05-15 | 2023-08-01 | 矿冶科技集团有限公司 | 冰铜品位软测量方法、装置、电子设备及存储介质 |
CN116307938A (zh) * | 2023-05-17 | 2023-06-23 | 成都瑞雪丰泰精密电子股份有限公司 | 一种加工中心进给系统健康状态评估方法 |
CN116449788A (zh) * | 2023-05-30 | 2023-07-18 | 湖南科技大学 | 一种钢厂铸坯生产过程故障检测与产能优化的方法 |
CN116823054A (zh) * | 2023-06-30 | 2023-09-29 | 浙大城市学院 | 一种基于数据挖掘的铝合金铸件制造质量实时预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108764517B (zh) | 一种高炉铁水硅含量变化趋势预测方法、设备和存储介质 | |
CN112100745B (zh) | 基于lda理论的汽车大梁钢力学性能预测方法 | |
CN114897227B (zh) | 基于改进随机森林算法的多钢种力学性能预报方法 | |
Li et al. | Ensemble machine learning systems for the estimation of steel quality control | |
CN111444942B (zh) | 一种高炉铁水硅含量智能预报方法及系统 | |
CN114678080A (zh) | 转炉终点磷含量预测模型及构建方法、磷含量预测方法 | |
CN113011491A (zh) | 一种主成分分析协同随机森林的热连轧带钢宽度预测方法 | |
CN113780420B (zh) | 基于gru-gcn的变压器油中溶解气体浓度预测方法 | |
CN113177358A (zh) | 一种基于模糊细粒度特征提取的水泥质量软测量方法 | |
CN115169453A (zh) | 基于密度聚类协同深度残差网络的热连轧宽度预测方法 | |
CN118197463A (zh) | 基于堆叠集成策略对冶金过程终点质量进行预测的方法 | |
CN112418522B (zh) | 一种基于三支集成预测模型的工业加热炉钢温预测方法 | |
CN114117919B (zh) | 基于样本协同表示的即时学习的软测量建模方法 | |
CN115952888A (zh) | 基于多变量灰色模型的能源碳排放预测方法及系统 | |
Li et al. | A recognition method of plate shape defect based on RBF-BP neural network optimized by genetic algorithm | |
CN117592856A (zh) | 基于石化企业产出量预测的算法流程 | |
US20240002964A1 (en) | Method and system for determining converter tapping quantity | |
CN115860211B (zh) | 一种基于局部在线建模的铸坯质量预测方法 | |
CN115860211A (zh) | 一种基于局部在线建模的铸坯质量预测方法 | |
CN114943328A (zh) | 基于bp神经网络非线性组合的sarima-gru时序预测模型 | |
CN114967625A (zh) | 一种基于时空卷积潜变量优化近邻成分分析的工业过程质量预测方法 | |
CN111626465A (zh) | 一种新能源功率短期区间预测方法及系统 | |
CN114798763A (zh) | 一种加热炉出钢板坯粗轧出口温度的预测方法及系统 | |
CN111027611A (zh) | 基于动态贝叶斯网络的模糊pls建模方法 | |
CN117732886B (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 |