基于非线性退化的轧辊剩余寿命预测实现方法技术

一种基于非线性退化的轧辊剩余寿命预测实现方法，其特征在于，通过构建基于核函数的非线性Wiener过程退化模型并初始化参数后，采用轧辊健康状态历史序列对模型进行训练；在在线阶段采用基于非线性Wiener过程退化模型的完整表达式，对待测轧辊退化数据进行计算，得到预测输出及其概率分布，并根据预测输出及其概率分布推导得轧辊剩余寿命的概率密度函数，实现剩余寿命预测。本发明专利技术通过对轧辊剩余寿命进行预测与分析，能够实现对轧钢生产线从当前的定期维护过渡到基于条件的维护，提高轧钢产品质量并降低生产成本。品质量并降低生产成本。品质量并降低生产成本。

Realization method of roll residual life prediction based on nonlinear degradation

【技术实现步骤摘要】
基于非线性退化的轧辊剩余寿命预测实现方法


[0001]本专利技术涉及的是一种轧钢生产领域的技术，具体是一种基于核函数
‑
Wiener过程的非线性退化的轧辊剩余寿命预测实现方法。

技术介绍

[0002]轧钢生产线是典型的连续工业生产线，即除了启停和异常情况之外，生产物料以连续状态通过整个生产流程，期间生产设备均处于运转状态。轧辊是轧机的关键部件和大型耗材，其更换冷却及维护操作会打断钢铁的连续生产进程，进而影响轧钢生产线的生产效率和经济效益。针对轧辊失效过程建立退化模型并实现剩余寿命预测，对于钢铁企业降低轧辊磨削损耗、提升轧线生产效率及产品质量具有重要意义。

技术实现思路

[0003]本专利技术针对现有技术需要利用大量的数据对相关模型进行训练才能进行操作，且提供的参考信息有限的缺陷，考虑轧辊退化路径中存在的随机性和非线性特点，提出一种基于非线性退化的轧辊剩余寿命预测实现方法，通过对轧辊剩余寿命进行预测与分析，能够实现对轧钢生产线从当前的定期维护过渡到基于条件的维护，提高轧钢产品质量并降低生产成本。
[0004]本专利技术是通过以下技术方案实现的：
[0005]本专利技术涉及一种基于非线性退化的轧辊剩余寿命预测实现方法，通过构建基于核函数的非线性Wiener过程退化模型并初始化参数后，采用轧辊健康状态历史序列对模型进行训练；在在线阶段采用基于非线性Wiener过程退化模型的完整表达式，对待测轧辊退化数据进行计算，得到预测输出及其概率分布，并根据预测输出及其概率分布推导得轧辊剩余寿命的概率密度函数，实现剩余寿命预测。
[0006]所述方法具体包括：
[0007]第一步，从更换表单/维修表单收集轧辊直径和轧制公里数等历史数据，并去除无效数据，获得轧辊健康状态历史序列将作为模型的输入数据；
[0008]第二步，构建非线性Wiener过程退化模型；
[0009]第三步，构建基于核函数的非线性Wiener过程退化模型的一般形式；
[0010]第四步，根据轧辊实测数据选取效果最好的核函数嵌入模型，得到非线性Wiener过程退化模型的具体形式；
[0011]第五步，基于轧辊退化历史数据和贝叶斯公式，推导模型未知参数w和σ的后验分布，获得边际似然函数，并通过最大化边缘似然函数，最终得到参数w和σ的最大似然估计值。
[0012]第六步，通过迭代估计算法求得非线性Wiener过程退化模型的近似解，作为模型参数的估计值，生成非线性Wiener过程退化模型的完整表达式；
[0013]所述的迭代求解，根据核函数性质删除迭代过程中部分权重和核函数，通过此过
程实现模型的稀疏性，提升了工程分析效率以适用轧辊性能状态的在线监测。
[0014]第七步：基于非线性Wiener过程退化模型的完整表达式，对待测轧辊退化数据进行计算，得到预测输出及其概率分布；
[0015]第八步：基于Wiener过程的首达时间服从逆高斯分布这一特性，根据预测输出及其概率分布推导得轧辊剩余寿命的概率密度函数，实现剩余寿命预测。
[0016]本专利技术涉及一种实现上述方法的系统，包括：数据预处理单元、模型构建单元、退化建模单元以及剩余寿命预测单元，其中：数据预处理单元根据生产线原始历史数据，进行数据清洗处理，得到模型的输入数据；模型构建单元根据处理完的轧辊实测数据，基于核函数的Wiener过程模型，进行贝叶斯参数估计，得到非线性Wiener过程退化模型；退化建模单元依据非线性Wiener过程退化模型计算得到未来时刻的预测输出及其概率分布，剩余寿命预测单元根据非线性Wiener过程退化模型和预测输出及其概率分布，得到轧辊当前剩余寿命的概率密度函数，从而提高设备运转效率与安全性并实现视情维护。技术效果
[0017]本专利技术考虑随机过程的同时结合核函数能够捕捉轧辊的复杂非线性退化路径的特点，构建基于核函数的Wiener过程退化模型，推导得到的失效分布蕴含首达时的概念，包含了可能失效的概率分布范围，为制定轧辊维修和备件计划提供更具体的决策依据，相较非线性Wiener过程和相关向量机算法，本专利技术在拟合优度上分别提高109％和14％，剩余寿命预测误差小于5％，具有较优秀的建模能力与预测效果。管理人员可以根据建模结果和剩余寿命制定维护计划，并可根据产线压力在预测区间内进行动态调节。
附图说明
[0018]图1为本专利技术流程图；
[0019]图2和图3为实施例效果示意图。
具体实施方式
[0020]如图1所示，本实施例结合某钢铁公司1580热轧生产线现场数据，进行基于非线性退化的轧辊剩余寿命预测，具体包括：
[0021]第一步，从更换表单/维修表达收集有效的轧辊直径和轧制公里数历史数据，并去除无效数据，结合轧辊健康指标获得轧辊的健康状态时间序列。
[0022]第二步，构建经典非线性Wiener过程退化模型X(t)＝ηΛ(t)+σB(t)，其中：σB(t)是扩散函数，Λ(t)为时间尺度转换函数，常用的为指数函数t
λ
，不适用轧辊的复杂非线性过程，需引入形式更灵活的漂移函数。
[0023]第三步，引入核函数解决轧辊退化过程中的非线性表征问题，替代经典非线性Wiener过程中的漂移函数Λ(t)，构建基于核函数的Wiener过程退化模型其中：核函数的一般形式为k(t，t
i
)＝Λ(t)
·
Λ(t
i
)，w
i
是权重系数，w0为截距，t
i
为i时刻健康状态时间序列的对应值；
[0024]第四步，根据轧辊实测数据的预测效果选取效果最好的核函数嵌入模型，具体为：使用均方根误差(RMSE)预测值与实际值间的均方根误差、用平均绝对误差(MAE)表示绝对
误差，其值越小则表明算法预测精度越高、用可决系数(CD)拟合优度，即相关系数的平方；通过对三者(CD的倒数)叠加评价模型的性能，其中：通过对三者(CD的倒数)叠加评价模型的性能，其中：其中：x(t)表示实际值，为模型预测值；
[0025]表1常用核函数的表达式
[0026]第五步，基于轧辊退化历史数据和贝叶斯公式，基于核函数
‑
Wiener过程退化模型的独立增量为x＝(Δx1，
…
，Δx
n
)
T
，其联合概率密度函数服从多元高斯分布，即：，其联合概率密度函数服从多元高斯分布，即：，其联合概率密度函数服从多元高斯分布，即：其中：其中其中∑是元素的协方差矩阵，∑(i，j)＝σ
2 min(t
i
，t
j
)；
[0027]对未知参数w和σ的后验分布进行推导，获得边际似然函数，并通过最大化边缘似然函数，最终得到参数w和σ的最大似然估计值，具体包括：
[0028]5.1)为避免模型出现过拟合现象，定义参数w的先验分布为均值为0的高斯分布，即：其中α＝(本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于非线性退化的轧辊剩余寿命预测实现方法，其特征在于，通过构建基于核函数的非线性Wiener过程退化模型并初始化参数后，采用轧辊健康状态历史序列对模型进行训练；在在线阶段采用基于非线性Wiener过程退化模型的完整表达式，对待测轧辊退化数据进行计算，得到预测输出及其概率分布，并根据预测输出及其概率分布推导得轧辊剩余寿命的概率密度函数，实现剩余寿命预测。2.根据权利要求1所述的基于非线性退化的轧辊剩余寿命预测实现方法，其特征是，所述的非线性Wiener过程退化模型，具体为：其中：σB(t)为扩散函数；核函数的一般形式为k(t，t
i
)＝Λ(t)
·
Λ(t
i
)，Λ(t)为时间尺度转换函数，w
i
是权重系数，w0为截距，w＝(w0，w1，
…
，w
n
)
T
，t
i
为i时刻健康状态时间序列的对应值。3.根据权利要求1所述的基于非线性退化的轧辊剩余寿命预测实现方法，其特征是，所述的轧辊健康状态历史序列，通过更换表单/维修表单收集轧辊直径和轧制公里数等历史数据，并去除无效数据，获得轧辊健康状态历史序列。4.根据权利要求1所述的基于非线性退化的轧辊剩余寿命预测实现方法，其特征是，所述的初始化参数是指：基于轧辊退化历史数据和贝叶斯公式，推导模型未知参数w和σ的后验分布，获得边际似然函数，并通过最大化边缘似然函数，最终得到参数w和σ的最大似然估计值。5.根据权利要求1所述的基于非线性退化的轧辊剩余寿命预测实现方法，其特征是，所述的训练是指：通过迭代估计算法，即根据核函数性质删除迭代过程中部分权重和核函数求得非线性Wiener过程退化模型的近似解，作为模型参数的估计值，生成非线性Wiener过程退化模型的完整表达式，从而实现模型的稀疏性，提升了工程分析效率以适用轧辊性能状态的在线监测。6.根据权利要求4所述的基于非线性退化的轧辊剩余寿命预测实现方法，其特征是，所述的初始化参数，具体包括：5.1)为避免模型出现过拟合现象，定义参数w的先验分布为均值为0的高斯分布，即：其中α＝(α0，α1，
…
，α
n
)是w的超参数向量；5.2)定义α和σ2的超先验分布服从Gamma分布，即p(β)＝Gamma(β|c，d)，a＝b＝c＝d＝10
‑4；5.3)基于贝叶斯公式，推导未知参数w、α和σ的后验分布为p(w，α，σ2|x)＝p(w|x，α，σ2)p(α，σ2|x)；5.4)最大化边际似然函数p(α，σ2|x)，使p(x|α，σ2)最大，得：)最大，得：其中其中7.根据权利要求5所述的基于非线性退化的轧辊剩余寿...

【专利技术属性】
技术研发人员：王汉禹，陈震，周笛，陈兆祥，潘尔顺，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：