CN104985003A - 一种基于数据驱动的轧制过程打滑诊断及预测的方法 - Google Patents
一种基于数据驱动的轧制过程打滑诊断及预测的方法 Download PDFInfo
- Publication number
- CN104985003A CN104985003A CN201510299592.6A CN201510299592A CN104985003A CN 104985003 A CN104985003 A CN 104985003A CN 201510299592 A CN201510299592 A CN 201510299592A CN 104985003 A CN104985003 A CN 104985003A
- Authority
- CN
- China
- Prior art keywords
- data
- rolling
- spe
- statistic
- pivot
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
一种基于数据驱动的轧制过程打滑诊断及预测的方法,利用正常轧制状态下的正常数据进行非线性分析,抽取核主元信息,建立核主元模型。将测试数据映射到核主元模型,检验测试数据的统计量T2和SPE是否超过控制限,判断轧机是否存在打滑故障。如果存在打滑故障,利用贡献图法对变量进行故障分离,准确判断打滑原因,修改相应轧制参数设定值,避免打滑现象的发生。本发明方法可在线监测轧制过程打滑故障的发生,还可对预设定的轧制规程进行诊断。
Description
技术领域
本发明涉及轧机自动控制领域,尤其是一种轧制过程打滑诊断及预防的方法。
背景技术
随着工业技术的进步,轧制设备继续向着大型化、高速化、自动化与高精度方向发展。与此同时轧件重量、轧制速度的不断增加和对产品质量要求的不断提高,轧制过程中主传动***打滑现象频繁发生。为了实现安全生产,必须采取有效的措施来避免打滑现象的发生。因此,如何避免轧机在高速轧制时发生打滑是提高轧机生产能力的一个重要课题。
到目前为止,国内外很多学者都对轧机的打滑现象做出一些研究。例如:白振华、王骏飞在“2002年薄钢板质量研讨会”上作的名为“冷连轧机打滑判断条件及其影响因素的研究”的报告,采用统计学的观点给出了打滑判断条件,并对冷连轧机打滑影响因素进行了分析,提出防治打滑的相关措施;陈俊在“2006年全国冷轧板带生产技术交流会”上作的关于“冷连轧机打滑问题的新解和思考”的报告,分析了轧制打滑是一个受张力制度、变形分配、工艺润滑冷却条件和轧制接触表面粗糙度等诸多工艺要素所控制的现象;张建成、蔡恒君等人在“2011年全国冷轧板带技术交流会”上发表的“冷连轧机打滑问题的分析”从生产实践的角度并结合相关理论阐述影响打滑的因素,提出了防治打滑的若干措施。但是,很少有通过核主元分析法对轧机在轧制过程中是否会出现打滑现象进行诊断。
发明内容
本发明目的在于提供一种容易判断、预测方法方便简单、打滑原因直观准确的基于数据驱动的轧制过程打滑诊断及预测的方法。
为实现上述目的,本发明所述方法包括以下步骤:
(1)采集数据,在轧制过程中,根据轧制理论和生产实际确定采集的正常数据,需采集的数据包括:乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度,对所采集的数据进行去噪和标准化处理;
(2)利用正常轧制状态下的正常数据进行非线性分析,计算正常状态下正常数据的核矩阵,对核矩阵进行中心化处理,建立打滑相关的核主元模型,利用自由搜索法确定核主元参数,计算正常数据的T2和SPE统计量,确定统计量T2和SPE的控制限;所述控制限对应于99%置信限;
其中,正常数据的T2统计量的表达式为:
T2=[t1,...,tp]Λ-1[t1,...,tp]T
式中,ti(i=1,...,p)是第i个非线性主元,p为主元个数,λi为特征值,Λ-1为相应特征值λi的倒数;
T2的控制限可根据F分布求得:
式中,n表示为样本数目,α表示为置信度,p为主元个数,Fp,n-p,α表示置信度为α的F分布;
SPE统计量的表达式为:
式中,φ(x)是原始样本数据x在特征空间F的映射,是特征空间的主元向量,n表示为样本数目,p为主元个数,ti(i=1,...,p)是第i个非线性主元;
正常数据的SPE的控制限为:
式中, λi为特征值,m为样本的维数,cα为标准正态分布在置信水平α下的阈值。
(3)将实际轧制过程的测试数据或轧制规程预设定值进行去噪和标准化处理,计算测试数据的核矩阵,对测试数据的核矩阵进行中心化处理,建立核主元模型,将测试数据或预设定值映射到核主元模型,计算测试数据的T2和SPE统计量;将测试数据的T2和SPE统计量与正常数据统计量T2和SPE的控制限进行比较,判断测试数据的统计量T2和SPE是否超过控制限,如果超限,则判断出现打滑故障,并输出打滑故障数据;
(4)对故障数据利用核主元贡献图法对变量进行故障分离,计算贡献量,画出贡献图,判断打滑故障原因,修改造成打滑故障的相应轧制变量,避免打滑现象发生。
进一步的,所述去噪处理采用小波去噪方法。
进一步的,利用核主元分析方法计算正常数据的T2和SPE统计量,
其中,核函数K(xi,xj)采用径向基函数,即:
式中,xi,xj∈x表示样本,c表示核主元参数;核主元参数c采用自由搜索法确定。
进一步的,实际轧制过程中的测试数据或轧制规程预设定值与历史数据的采集种类相同,均包括乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度。
应用本发明方法的工作过程大致如下:
利用正常轧制状态下的历史数据进行非线性分析,抽取核主元信息,建立核主元模型。将测试数据(预设定轧制规程的相关数据或监测的实际生产数据)映射到核主元模型,判断轧机是否存在打滑故障。如果存在打滑故障,利用核主元贡献图法对变量进行故障分离,准确判断打滑原因,修改相应轧制参数设定值,避免打滑现象的发生。
与现有技术相比,本发明方法具有如下优点:提高生产过程的安全性并使检修得到优化;在线监测轧制过程打滑故障的发生,还可对预设定的轧制规程进行诊断,如果存在问题可及时修正,防患于未然,提高生产安全性和生产效率。
附图说明
图1是本发明方法的流程图。
图2a是实施例1各采样点的T2统计量。
图2b是实施例1各采样点的SPE统计量。
图3a是实施例1的各变量T2贡献图。
图3b是实施例1的各变量SPE贡献图。
图4a是实施例1三个变量改为正常值后的T2仿真图。
图4b是实施例1三个变量改为正常值后的SPE仿真图。
具体实施方式
下面结合附图对本发明做进一步说明:
如图1所示,利用正常轧制状态下的正常数据进行非线性分析,抽取核主元信息,建立核主元模型。将测试数据映射到核主元模型,检验测试数据的统计量T2和SPE是否超过控制限,判断轧机是否存在打滑故障。如果存在打滑故障,利用核主元贡献图法对变量进行故障分离,准确判断打滑原因,修改相应轧制参数设定值,避免打滑现象的发生。
具体步骤如下:
(1)采集数据,在轧制过程中,根据轧制理论和生产实际确定采集的正常数据,需采集的数据包括:乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度,把数据分为正常数据和测试数据,对所采集的数据进行小波去噪和标准化处理;
将所有变量数据通过标准化转化为均值为0,方差为1的无量纲数值。其计算公式为:
式中,yi为标准化后的数据,xi为采集的原始数据,为平均值:n是样本数量,s为标准偏差:
(2)利用正常轧制状态下的正常数据进行非线性分析,计算正常状态下历史数据的核矩阵,对核矩阵进行中心化处理,建立打滑相关的核主元模型,利用自由搜索法确定核主元参数,计算历史数据的T2和SPE统计量,确定统计量T2和SPE的控制限;所述控制限对应于99%置信限;
正常数据的T2统计量的表达式为:
T2=[t1,...,tp]Λ-1[t1,...,tp]T
式中,ti(i=1,...,p)是第i个非线性主元,p为主元个数,λi为特征值,Λ-1为相应特征值λi的倒数;
T2的控制限可根据F分布求得:
式中,n表示为样本数目,α表示为置信度,p为主元个数,Fp,n-p,α表示置信度为α的F分布;
正常数据的SPE统计量的表达式为:
式中,φ(x)是原始样本数据x在特征空间F的映射,是特征空间的主元向量,n表示为样本数目,p为主元个数,ti(i=1,...,p)是第i个非线性主元;
SPE的控制限为:
式中, λi为特征值,m为样本的维数,cα为标准正态分布在置信水平α下的阈值;
核函数K(xi,xj)采用径向基函数:
式中,xi,xj∈x表示样本,c表示核主元参数;核主元参数c采用自由搜索法确定。
用核函数将正常轧制状态下样本变换成为核函数矩阵,将输入数据映射到高维特征空间,建立打滑相关的核主元模型,计算正常历史数据的统计量T2和SPE,确定统计量T2和SPE的控制限(对应于99%置信限)。
(3)将实际轧制过程的测试数据或轧制规程预设定值进行去噪和标准化处理,测试数据或轧制规程预设定值包括乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度;计算测试数据的核矩阵,对测试数据的核矩阵进行中心化处理,建立核主元模型,将测试数据或预设定值映射到核主元模型,计算测试数据的T2和SPE统计量;将测试数据的T2和SPE统计量与正常历史数据统计量T2和SPE的控制限进行比较,判断测试数据的统计量T2和SPE是否超过控制限,如果超限,则判断出现打滑故障,并输出打滑故障数据;
(4)利用核主元贡献图法对变量进行故障分离,计算贡献量,画出贡献图,判断打滑故障原因,对故障数据的统计量T2和SPE对每个变量求偏导数,进而得到每个变量对统计量T2和SPE的贡献率,对统计量T2和SPE影响较大的变量就是故障变量,修改造成打滑故障的相应轧制变量,避免打滑现象发生。
实施例1:以某一实际工业环境下轧制过程采集的数据进行打滑故障诊断,67采样点开始出现打滑想象,图2a和图2b是测试数据的T2和SPE统计量,从67采样点开始超限;对故障数据进行故障分离,图3a和图3b是67采样点的各变量的T2和SPE贡献图;确定有三个变量2,4,6出现故障,对这三个变量进行修改,避免故障发生。图4a和图4b是三个变量修改后的T2和SPE统计量,从图4可以看出,对三个故障变量进行修改后,T2和SPE统计量都低于控制限,可明显消除打滑故障。
以上所述的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (5)
1.一种基于数据驱动的轧制过程打滑诊断及预测的方法,其特征在于,所述方法包括以下步骤:
(1)采集数据,在轧制过程中,根据轧制理论和生产实际确定采集的正常数据,需采集的数据包括:乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度,对所采集的数据进行去噪和标准化处理;
(2)利用正常轧制状态下的正常数据进行非线性分析,计算正常状态下正常数据的核矩阵,对核矩阵进行中心化处理,建立打滑相关的核主元模型,利用自由搜索法确定核主元参数,计算正常数据的T2和SPE统计量,确定统计量T2和SPE的控制限;所述控制限对应于99%置信限;
(3)将实际轧制过程的测试数据或轧制规程预设定值进行去噪和标准化处理,计算测试数据的核矩阵,对测试数据的核矩阵进行中心化处理,建立核主元模型,将测试数据或预设定值映射到核主元模型,计算测试数据的T2和SPE统计量;将测试数据的T2和SPE统计量与正常数据统计量T2和SPE的控制限进行比较,判断测试数据的统计量T2和SPE是否超过控制限,如果超限,则判断出现打滑故障,并输出打滑故障数据;
(4)对故障数据利用核主元贡献图法对变量进行故障分离,计算贡献量,画出贡献图,判断打滑故障原因,修改造成打滑故障的相应轧制变量,避免打滑现象发生。
2.根据权利要求1所述的基于数据驱动的轧制过程打滑诊断及预测的方法,其特征在于:步骤(1)中所述去噪处理采用小波去噪方法。
3.根据权利要求1所述的基于数据驱动的轧制过程打滑诊断及预测的方法,其特征在于,步骤(2)中,正常数据的T2统计量的表达式为:
T2=[t1,...,tp]Λ-1[t1,...,tp]T
式中,ti(i=1,...,p)是第i个非线性主元,p为主元个数,λi为特征值,Λ-1为相应特征值λi的倒数;
T2的控制限可根据F分布求得:
式中,n表示为样本数目,α表示为置信度,p为主元个数,Fp,n-p,α表示置信度为α的F分布;
正常数据的SPE统计量的表达式为:
式中,φ(x)是原始样本数据x在特征空间F的映射,是特征空间的主元向量,n表示为样本数目,p为主元个数,ti(i=1,...,p)是第i个非线性主元。
SPE的控制限为:
式中, λi为特征值,m为样本的维数,cα为标准正态分布在置信水平α下的阈值。
4.根据权利要求1所述的基于数据驱动的轧制过程打滑诊断及预测的方法,其特征在于:步骤(2)中利用核主元分析方法计算正常数据的T2和SPE统计量,
其中,核函数K(xi,xj)采用径向基函数,即:
式中,xi,xj∈x表示样本,c表示核主元参数;核主元参数c采用自由搜索法确定。
5.根据权利要求1所述的基于数据驱动的轧制过程打滑诊断及预测的方法,其特征在于:步骤(3)所述实际轧制过程中的测试数据或轧制规程预设定值包括乳化液浓度、轧件温度、轧机的前后张力、轧辊速度、轧机的入口厚度和出口厚度、轧辊直径、轧件宽度、轧辊表面粗糙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510299592.6A CN104985003A (zh) | 2015-06-03 | 2015-06-03 | 一种基于数据驱动的轧制过程打滑诊断及预测的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510299592.6A CN104985003A (zh) | 2015-06-03 | 2015-06-03 | 一种基于数据驱动的轧制过程打滑诊断及预测的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104985003A true CN104985003A (zh) | 2015-10-21 |
Family
ID=54296965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510299592.6A Pending CN104985003A (zh) | 2015-06-03 | 2015-06-03 | 一种基于数据驱动的轧制过程打滑诊断及预测的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104985003A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106762452A (zh) * | 2016-12-08 | 2017-05-31 | 东北大学 | 基于数据驱动的风机主控***故障诊断和在线监测方法 |
CN109675935A (zh) * | 2019-03-06 | 2019-04-26 | 北京科技大学 | 一种变控制限的ipca轧制过程在线故障诊断方法 |
CN111069286A (zh) * | 2019-12-05 | 2020-04-28 | 唐山不锈钢有限责任公司 | 一种基于轧制过程参数分析判断的轧件打滑自动纠正方法 |
CN115591947A (zh) * | 2022-12-15 | 2023-01-13 | 太原科技大学(Cn) | 一种连轧过程板带质量分布式调控方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1655082A (zh) * | 2005-01-27 | 2005-08-17 | 上海交通大学 | 基于核主元分析的非线性故障诊断的方法 |
CN101739514A (zh) * | 2010-01-02 | 2010-06-16 | 燕山大学 | 双ucm机型的二次冷轧机组轧制工艺参数综合优化方法 |
CN101927261A (zh) * | 2010-01-02 | 2010-12-29 | 燕山大学 | 二次冷轧机组轧制模式下工艺润滑制度综合优化方法 |
CN103197663A (zh) * | 2013-03-07 | 2013-07-10 | 北京信息科技大学 | 一种故障预测方法及*** |
-
2015
- 2015-06-03 CN CN201510299592.6A patent/CN104985003A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1655082A (zh) * | 2005-01-27 | 2005-08-17 | 上海交通大学 | 基于核主元分析的非线性故障诊断的方法 |
CN101739514A (zh) * | 2010-01-02 | 2010-06-16 | 燕山大学 | 双ucm机型的二次冷轧机组轧制工艺参数综合优化方法 |
CN101927261A (zh) * | 2010-01-02 | 2010-12-29 | 燕山大学 | 二次冷轧机组轧制模式下工艺润滑制度综合优化方法 |
CN103197663A (zh) * | 2013-03-07 | 2013-07-10 | 北京信息科技大学 | 一种故障预测方法及*** |
Non-Patent Citations (2)
Title |
---|
刘美玲: "基于多元统计分析的过程***故障诊断方法研究", 《中国优秀硕士学位论文全文数据库》 * |
白振华等: "冷连轧机的打滑判断条件及打滑防治措施", 《上海金属》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106762452A (zh) * | 2016-12-08 | 2017-05-31 | 东北大学 | 基于数据驱动的风机主控***故障诊断和在线监测方法 |
CN106762452B (zh) * | 2016-12-08 | 2019-04-23 | 东北大学 | 基于数据驱动的风机主控***故障诊断和在线监测方法 |
CN109675935A (zh) * | 2019-03-06 | 2019-04-26 | 北京科技大学 | 一种变控制限的ipca轧制过程在线故障诊断方法 |
CN109675935B (zh) * | 2019-03-06 | 2020-07-31 | 北京科技大学 | 一种变控制限的ipca轧制过程在线故障诊断方法 |
CN111069286A (zh) * | 2019-12-05 | 2020-04-28 | 唐山不锈钢有限责任公司 | 一种基于轧制过程参数分析判断的轧件打滑自动纠正方法 |
CN115591947A (zh) * | 2022-12-15 | 2023-01-13 | 太原科技大学(Cn) | 一种连轧过程板带质量分布式调控方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102033523B (zh) | 基于偏最小二乘的带钢质量预测、炉况预警与故障诊断方法 | |
CN103793854B (zh) | 多重组合优化的架空输电线路运行风险信息化评估方法 | |
Ma et al. | A novel data-based quality-related fault diagnosis scheme for fault detection and root cause diagnosis with application to hot strip mill process | |
Peng et al. | Contribution rate plot for nonlinear quality-related fault diagnosis with application to the hot strip mill process | |
CN103853152B (zh) | 一种基于ar-pca的间歇过程故障监测方法 | |
CN104985003A (zh) | 一种基于数据驱动的轧制过程打滑诊断及预测的方法 | |
CN102252843B (zh) | 一种滚动轴承性能变异的评估方法 | |
CN104777830A (zh) | 一种基于kpca混合模型的多工况过程监控方法 | |
CN110262450B (zh) | 面向汽轮机的多种故障特性协同分析的故障预测方法 | |
CN106391725B (zh) | 一种适用于冷轧过程轧制压力随轧制速度变化预报方法 | |
CN106845826B (zh) | 一种基于PCA-Cpk的冷连轧生产线服役质量状态评估方法 | |
CN104731083B (zh) | 一种基于自适应特征提取的工业故障诊断方法及应用 | |
CN104050340A (zh) | 一种大型数控铣床刀具磨损程度识别方法 | |
CN107862324B (zh) | 一种基于mwspca的cbr预测模型智能化预警方法 | |
CN111209967A (zh) | 一种基于支持向量机的轧制过程板凸度预测方法 | |
CN104714537A (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
CN106483847A (zh) | 一种基于自适应ica的冷水机组故障检测方法 | |
CN107808257A (zh) | 一种对外贸易风险评估预警***及其实现方法 | |
CN111797533B (zh) | 一种核动力装置运行参数异常检测方法及*** | |
CN106528938A (zh) | 一种基于按类统计的滚动轴承可靠性可视化动态评估方法 | |
CN106569052A (zh) | 考虑实时健康状态的电力变压器可靠性评估方法 | |
CN109719138B (zh) | 一种基于数据挖掘的变形抗力唯象模型计算方法 | |
CN103487250A (zh) | 基于二维投射的煤矿设备预知维护方法 | |
CN104408525A (zh) | 作业车间调度风险的量化评估与控制方法 | |
CN104598970B (zh) | 一种爬架组的工作状态检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151021 |
|
WD01 | Invention patent application deemed withdrawn after publication |