一种基于数据驱动的轧制过程打滑诊断及预测的方法
【技术领域】
[0001] 本发明涉及轧机自动控制领域，尤其是一种轧制过程打滑诊断及预防的方法。
【背景技术】
[0002] 随着工业技术的进步，乳制设备继续向着大型化、高速化、自动化与高精度方向发 展。与此同时轧件重量、乳制速度的不断增加和对产品质量要求的不断提高，乳制过程中主 传动系统打滑现象频繁发生。为了实现安全生产，必须采取有效的措施来避免打滑现象的 发生。因此，如何避免轧机在高速轧制时发生打滑是提高轧机生产能力的一个重要课题。
[0003] 到目前为止，国内外很多学者都对轧机的打滑现象做出一些研宄。例如：白振华、 王骏飞在"2002年薄钢板质量研讨会"上作的名为"冷连轧机打滑判断条件及其影响因素 的研宄"的报告，采用统计学的观点给出了打滑判断条件，并对冷连轧机打滑影响因素进行 了分析，提出防治打滑的相关措施；陈俊在"2006年全国冷轧板带生产技术交流会"上作的 关于"冷连轧机打滑问题的新解和思考"的报告，分析了轧制打滑是一个受张力制度、变形 分配、工艺润滑冷却条件和轧制接触表面粗糙度等诸多工艺要素所控制的现象；张建成、蔡 恒君等人在"2011年全国冷轧板带技术交流会"上发表的"冷连轧机打滑问题的分析"从生 产实践的角度并结合相关理论阐述影响打滑的因素，提出了防治打滑的若干措施。但是，很 少有通过核主元分析法对轧机在轧制过程中是否会出现打滑现象进行诊断。
【发明内容】
[0004] 本发明目的在于提供一种容易判断、预测方法方便简单、打滑原因直观准确的基 于数据驱动的轧制过程打滑诊断及预测的方法。
[0005] 为实现上述目的，本发明所述方法包括以下步骤：
[0006] (1)采集数据，在轧制过程中，根据轧制理论和生产实际确定采集的正常数据，需 采集的数据包括：乳化液浓度、乳件温度、乳机的前后张力、乳辊速度、乳机的入口厚度和出 口厚度、乳辊直径、乳件宽度、乳辊表面粗糙度，对所采集的数据进行去噪和标准化处理；
[0007] (2)利用正常轧制状态下的正常数据进行非线性分析，计算正常状态下正常数据 的核矩阵，对核矩阵进行中心化处理，建立打滑相关的核主元模型，利用自由搜索法确定核 主元参数，计算正常数据的T 2和SPE统计量，确定统计量T2和SPE的控制限；所述控制限对 应于99 %置信限；
[0008] 其中，正常数据的T2统计量的表达式为：
[0009] T2= A^ 1Lt1, ...,tp]T
[0010] 式中，ti(i = I, . . . , p)是第i个非线性主元，P为主元个数，A i为特征值，A H为 相应特征值Ai的倒数；
[0011] T2的控制限可根据F分布求得：
[0013] 式中，n表示为样本数目，a表示为置信度，P为主元个数，Fp，n_p， a表示置信度为 a的F分布；
[0014] SPE统计量的表达式为：
[0016] 式中，巾（X)是原始样本数据X在特征空间F的映射，^p(X)是特征空间的主元向 量，n表示为样本数目，p为主元个数，h (i = 1，. . .，p)是第i个非线性主元；
[0017] 正常数据的SPE的控制限为：
Ca为标准正态分布在置信水平a下的阈值。
[0020] (3)将实际轧制过程的测试数据或轧制规程预设定值进行去噪和标准化处理，计 算测试数据的核矩阵，对测试数据的核矩阵进行中心化处理，建立核主元模型，将测试数据 或预设定值映射到核主元模型，计算测试数据的T 2和SPE统计量；将测试数据的T 2和SPE 统计量与正常数据统计量T2和SPE的控制限进行比较，判断测试数据的统计量T 2和SPE是 否超过控制限，如果超限，则判断出现打滑故障，并输出打滑故障数据；
[0021] (4)对故障数据利用核主元贡献图法对变量进行故障分离，计算贡献量，画出贡献 图，判断打滑故障原因，修改造成打滑故障的相应轧制变量，避免打滑现象发生。
[0022] 进一步的，所述去噪处理采用小波去噪方法。
[0023] 进一步的，利用核主元分析方法计算正常数据的T2和SPE统计量，
[0024] 其中，核函数K(Xi，Xj)采用径向基函数，即：
[0026] 式中，Xi, XjG X表示样本，c表示核主元参数；核主元参数c采用自由搜索法确定。
[0027] 进一步的，实际轧制过程中的测试数据或轧制规程预设定值与历史数据的采集种 类相同，均包括乳化液浓度、乳件温度、乳机的前后张力、乳辊速度、乳机的入口厚度和出口 厚度、乳辊直径、乳件宽度、乳辊表面粗糙度。
[0028] 应用本发明方法的工作过程大致如下：
[0029] 利用正常轧制状态下的历史数据进行非线性分析，抽取核主元信息，建立核主元 模型。将测试数据（预设定轧制规程的相关数据或监测的实际生产数据）映射到核主元模 型，判断轧机是否存在打滑故障。如果存在打滑故障，利用核主元贡献图法对变量进行故障 分离，准确判断打滑原因，修改相应轧制参数设定值，避免打滑现象的发生。
[0030] 与现有技术相比，本发明方法具有如下优点：提高生产过程的安全性并使检修得 到优化；在线监测轧制过程打滑故障的发生，还可对预设定的轧制规程进行诊断，如果存在 问题可及时修正，防患于未然，提高生产安全性和生产效率。
【附图说明】
[0031] 图1是本发明方法的流程图。
[0032] 图2a是实施例1各采样点的T2统计量。
[0033] 图2b是实施例1各采样点的SPE统计量。
[0034] 图3a是实施例1的各变量T2贡献图。
[0035] 图3b是实施例1的各变量SPE贡献图。
[0036] 图4a是实施例1三个变量改为正常值后的T2仿真图。
[0037] 图4b是实施例1三个变量改为正常值后的SPE仿真图。
【具体实施方式】
[0038] 下面结合附图对本发明做进一步说明：
[0039] 如图1所示，利用正常轧制状态下的正常数据进行非线性分析，抽取核主元信息， 建立核主元模型。将测试数据映射到核主元模型，检验测试数据的统计量T 2和SPE是否超 过控制限，判断轧机是否存在打滑故障。如果存在打滑故障，利用核主元贡献图法对变量进 行故障分离，准确判断打滑原因，修改相应轧制参数设定值，避免打滑现象的发生。
[0040] 具体步骤如下：
[0041] (1)采集数据，在轧制过程中，根据轧制理论和生产实际确定采集的正常数据，需 采集的数据包括：乳化液浓度、乳件温度、乳机的前后张力、乳辊速度、乳机的入口厚度和出 口厚度、乳辊直径、乳件宽度、乳辊表面粗糙度，把数据分为正常数据和测试数据，对所采集 的数据进行小波去噪和标准化处理；
[0042] 将所有变量数据通过标准化转化为均值为0,方差为1的无量纲数值。其计算公式 为：
[0045] (2)利用正常轧制状态下的正常数据进行非线性分析，计算正常状态下历史数据 的核矩阵，对核矩阵进行中心化处理，建立打滑相关的核主元模型，利用自由搜索法确定核 主元参数，