一种基于特征向量和层次聚类的结晶器漏钢预报方法

摘要

本发明属于钢铁冶金连铸检测技术领域，提供一种基于特征向量和层次聚类的结晶器漏钢预报方法。该预报方法分别提取黏结漏钢、正常工况历史数据以及在线实测数据的温度特征向量，建立特征向量样本集；对样本集进行归一化处理，并进行层次聚类；此后检查和判断在线提取的特征向量是否从属于漏钢类簇，进而识别和预报结晶器漏钢。本发明回避了报警阈值等参数繁琐的调试和修改环节，克服了以往漏钢预报方法的人为依赖性，具有良好的鲁棒性和迁移性；通过温度特征提取，不仅可准确识别黏结漏钢温度模式，避免漏报并显著降低了误报次数，还能够极大压缩数据计算量和运算时间，确保在线预报的实时性。

说明书

技术领域

本发明属于钢铁冶金连铸检测技术领域，涉及一种基于特征向量和层次聚类的结晶器漏钢预报方法。

背景技术

黏结漏钢是指连铸过程中弯月面附近较薄的初生坯壳发生破裂，钢液渗出后与结晶器铜板接触发生黏结，随结晶器振动和铸坯下移，黏结反复撕裂-愈合且不断下移，当其移出结晶器出口后失去了铜板的支撑约束，钢液溢出，造成漏钢。漏钢不仅危及现场操作人员的安全，严重损坏连铸设备，同时将导致连铸生产被迫中断，设备维修和生产成本大幅上升。因此，针对漏钢的在线监控和预报始终是连铸过程控制的重中之重，对保障生产顺行具有重要意义。

目前，现有的漏钢预报方法主要通过在结晶器铜板上嵌入安装测温热电偶，根据热电偶温度的变化来监控和判别铸坯与铜板之间是否发生黏结。实践表明，黏结发生时热电偶温度随时间的变化率及其幅值均与正常工况存在明显的差异，因此，可通过提取和归纳热电偶温度及其速率、幅值的共性特征，结合逻辑判断或神经网络方法，区分和识别黏结的典型温度特征，在线预报结晶器漏钢。

从现有漏钢预报技术的实际应用情况来看，基于逻辑判断的漏钢预报方法对连铸设备、工艺参数和物性参数的依赖性较强，当工艺调整和拉速提升时，阈值变动大，导致误报率和漏报率大幅上升；神经网络方法对学习和训练样本的要求较高，样本不全或无效时都会严重影响其预报效果，模型的迁移能力较低。实践中为尽力避免漏报，在设计预报算法或设置报警阈值时都会预留出一定余赋，在工艺调整、设备维护更换等情况下，预报阈值需人为频繁调试和修正，若设置不当或调整不及时，误报也很难准确识别和剔除，甚至导致漏报，是目前漏钢预报系统面临的常见问题。

专利CN106980729A公开了一种基于混合模型的连铸漏钢预报方法，该方法采集和存储现场所有热电偶温度的实时数据，判断每个热电偶温度在时间序列上的变化是否符合黏结漏钢时的温度变化波形，保存判断结果为Y(i,j,t)，如果Y(i,j,t)在设定的阈值范围内则标记为热电偶异常，并统计当前热电偶所在行和上一行异常热电偶的数目，比较输出的异常热电偶总数与黏结警告和黏结报警的热电偶数目阈值。该方法侧重单支热电偶的温度变化和异常电偶的叠加数目，对黏结漏钢重要的“时滞性”并没有涉及，使得准确捕捉漏钢的“时滞”和“温度倒置”特征较为困难，影响了该方法的预报精度。

专利CN102554171A公开了一种连铸漏钢预报方法，将自适应算法引入BP神经网络，以实现网络结构的自动优化，并提出了一种基于逻辑判断和神经网络的摩擦力监控漏钢预报模型，有机结合温度监控的准确性与摩擦力监控的灵敏性，建立了以温度监控为主、摩擦力监控为辅的预报机制。该方法设有温度监控单偶、组偶以及摩擦力逻辑判断、时序网络四处报警点，且均有黄、橙、红三种级别，每处报警点均需设置相应的阈值，预报参数众多，使得该方法漏钢预报的普适性和迁移性较差。

专利CN108580827A公开了一种基于凝聚层次聚类预报结晶器漏钢的方法，通过建立黏结漏钢和正常工况样本库，从黏结漏钢样本集和正常工况样本集中各自随机取选取等量样本，与在线实测温度样本构成随机样本集，对该随机样本集实施层次聚类，之后，检测在线实测温度样本是否属于黏结漏钢类簇，以此识别和预报漏钢。该发明摆脱了在预报过程中对人为定义参数的依赖性，仅利用黏结漏钢和正常工况温度的特征判断在线实测温度样本是否包含属于漏钢。但该方法的局限在于：1)在所有漏钢样本集中随机选择部分样本进行聚类，无法兼顾全部样本的共性特征和单个样本的个性特征，即聚类后样本特征覆盖不全面；2)温度数据经预处理后直接用于聚类，数据量和计算量均较大，影响在线预报的实时性。针对上述问题，本发明一方面通过温度特征提取，降低数据的维度和运算量，另一方面，采用包括漏钢、正常工况的全体样本进行层次聚类，以兼顾样本的共性特征和个性特征，在确保预报准确性的前提下，提升在线运算速度和实时性。

发明内容

本发明的目的是为了克服现有预报方法在漏钢特征提取方面存在的不足，提出一种基于特征向量和层次聚类的结晶器漏钢预报方法，通过提取全部黏结漏钢样本、正常工况历史数据以及在线实测数据的温度特征向量，建立特征向量库；对特征向量库进行归一化处理和层次聚类；此后检查和判断在线提取的特征向量是否从属于漏钢类簇，进而识别和预报结晶器漏钢。

为达到上述目的，本发明的技术方案如下：

一种基于特征向量和层次聚类的结晶器漏钢预报方法，包括以下步骤：

第一步，提取黏结漏钢特征向量

(1)获取黏结漏钢历史温度数据：标记黏结位置所在电偶列第一排电偶温度最高时刻，并选取其前M秒、后N-1秒共计M+N秒的温度数据；

(2)提取和构建第一排、第二排电偶温度特征向量：分别以每列电偶为单位，提取黏结漏钢时同列电偶温度沿浇铸方向的变化特征，具体包括：

1st_Rising_Amplitude：第一排温度上升幅值；

1st_Rising_V_Max：第一排温度上升速率的最大值；

1st_Falling_V_Ave：第一排温度下降速率均值；

2nd_Rising_V_Max：第二排温度上升速率的最大值；

1st_2nd_Time_Lag：温度上升时滞，即第二排与第一排热电偶温度开始上升时刻的时间间隔；

由此构建特征向量：

s＝[1st_Rising_Amplitude,1st_Rising_V-Max,1st-Falling-V-Ave,2nd-Rising-V-Max,1st-2nd-Time-Lag]

第二步，提取正常工况特征向量

(1)获取正常工况历史温度数据：任意截取连续M+N秒的温度数据；提取和构建第一排、第二排电偶温度特征向量：分别以每列电偶为单位，提取正常工况下同列电偶温度沿浇铸方向的变化特征，具体方法同第一步步骤(2)；

第三步，提取在线实时温度特征

(1)实时采集和获取结晶器内弧宽面、外弧宽面、左侧窄面、右侧窄面铜板各行、列热电偶当前时刻及之前M+N-1秒，共计M+N秒的温度数据；

(2)提取和构建第一排、第二排电偶温度特征向量：分别以每列电偶为单位，提取在线实时温度同列电偶温度沿浇铸方向的变化特征，具体方法同第一步步骤(2)；

第四步，建立特征向量库

(1)依据第一、第二、第三步提取的黏结漏钢、正常工况与在线实测的温度特征建立特征向量样本库D；

(2)对特征向量样本库D作归一化处理，得到特征向量集S，将归一化后的在线实测温度特征记为s_new，具体归一化方法如下：

其中，x_ij为特征向量集S中第i个特征向量第j维特征的值，x_jmax、x_jmin分别为所有特征向量第j维特征的最大值和最小值，|S|表示特征向量集S中向量的总数。

第五步，特征向量层次聚类

(1)对第四步得到的特征向量集S实施层次聚类，具体过程包括：

1.1)将特征向量集S中的每个向量s看作一个初始类簇C_i＝{s_i}，建立类簇集合C＝{C₁,C₂,...,C_k}；其中s_i表示S中第i个向量，C_i表示第i个类簇，i＝1,2,...,k，k表示特征向量集S中向量的总数；

1.2)计算和确定类簇集合C中任意两个类簇C_p、C_q的距离：

d(C_p,C_q)＝min(dist(C_pi,C_qj))

其中，C_pi为类簇C_p中第i个特征向量，C_qj为类簇C_q中第j个特征向量，dist(C_pi,C_qj)表示特征向量C_pi、C_qj的欧氏距离；计算类簇C_p、C_q中任意两个向量的距离，取距离的最小值min作为类簇C_p和C_q的距离。

1.3)标记经步骤1.2)计算后距离最小的两个类簇C_m和C_n，将C_m和C_n合并成为一个新的类簇C_{m,n}并加入集合C，同时删除原有的类簇C_m和C_n，经过类簇的添加和删除后，此时集合C中类簇的总数目减1；

1.4)循环执行步骤1.2)～1.3)，当类簇集合C中仅剩两个类簇时，结束循环，完成聚类过程；

(2)检查聚类结果是否满足以下判定条件，即：

所有黏结漏钢特征向量的90％以上属于同一个类簇，且该类簇中正常工况特征向量所占的比率低于20％；

满足此条件则将该类簇记为漏钢类簇C_breakout，另一个类簇记为正常工况类簇C_normal；否则，重新执行步骤(1)、(2)，直至由漏钢、正常工况以及实测温度构成的特征向量集聚类结果满足上述判定条件为止；

第六步，漏钢识别与报警

判断新特征向量s_new是否属于类簇C_breakout，如果是，则发出漏钢警报；否则，继续执行第三、四、五、六步。

上述预报漏钢的方法适用于板坯、方坯、圆坯、异型坯等连铸坯漏钢识别。

本发明的有益效果是：所提出的一种基于特征向量和层次聚类的结晶器漏钢预报方法，回避了报警阈值等参数繁琐的调试和修改环节，克服了以往漏钢预报方法的人为依赖性，具有良好的鲁棒性和迁移性；通过温度特征提取，不仅可准确识别黏结漏钢温度模式，避免漏报并显著降低了误报次数，还极大压缩了数据计算量和运算时间，确保在线预报的实时性。

附图说明

图1是四张结晶器铜板与热电偶分布示意图；

图2是黏结漏钢时温度及其变化率的特征向量提取示意图；

图3是正常工况下温度及其变化率的特征向量提取示意图；

图4是特征向量集层次聚类及漏钢识别与报警流程图；

图5是在线实测温度1；

图6是包含在线实测温度特征向量1的层次聚类及漏钢预报结果；

图7是在线实测温度2；

图8是包含在线实测温度特征向量2的层次聚类及漏钢预报结果。

具体实施方式

下面通过具体的实施例，并结合附图对本发明做进一步的阐述。

本发明主要由六个部分构成：提取黏结漏钢特征向量、提取正常工况特征向量、提取在线实时温度特征向量、建立特征向量库、特征向量层次聚类、漏钢识别与报警。

图1是四张结晶器铜板与热电偶分布示意图。板坯连铸结晶器由四张铜板构成，包括外弧宽面铜板、左侧窄面铜板、内弧宽面铜板、右侧窄面铜板长度L为900mm，分别在四张铜板距结晶器上口L1为210mm、L2为325mm的水平横截面上布置2行测点，外弧宽面铜板和内弧宽面铜板上每行布置19列热电偶，两个电偶间距L3为150mm，两张宽面铜板各布置38支热电偶；左侧窄面铜板和右侧窄面铜板位于中心线各布置一列热电偶，两张窄面铜板各布置2支热电偶。四张铜板布置电偶总数共计80支，每支电偶至结晶器铜板热面距离相等。

步骤一、提取黏结漏钢特征向量

(1)获取黏结漏钢历史温度数据：标记黏结位置所在电偶列第一排电偶温度最高时刻，并选取其前15秒、后9秒共计25秒的温度数据；

对于黏结漏钢历史温度，选取30例温度样本。

(2)提取和构建第一排、第二排电偶温度特征向量

图2为黏结漏钢温度及其温度变化率的特征向量提取示意图。由图2可知，特征向量的提取分别以每列电偶为单位，提取黏结漏钢时同列电偶温度沿浇铸方向的变化特征，具体特征及其提取方法如下：

2.1)计算同一测点处温度数据在5秒内的变化率，即：

式中，r∈[1,2]分别表示第一排、第二排热电偶，T_(r)i表示第r排热电偶第i时刻的温度数值，v_(r)i表示第r排热电偶第i时刻温度变化率的数值。

2.2)确定温度上升时的起始温度及其对应的时刻：

首先，获取包含当前时刻及之前共计25秒内温度的最大值T_max及其对应的时刻t_max，即：

T_max＝max(T_i),i＝1,2,...,25，

然后，从T_max前向遍历，获取其之前温度的最小值T_min及时刻t_min，即：

T_min＝min(T_i),i＝1,2,...,t_max

则取T_min为温度上升时的起始温度，其对应的时刻记为t_min。

2.3)提取相应的特征构建特征向量：

1st_Rising_Amplitude：第一排温度上升幅值，即分别标记第一排温度开始上升及达到最大时的温度，计算两温度间的差值，得到温度上升幅值，单位为℃：

1st-Rising-Amplitude＝T_(1)max-T_(1)min

1st_Rising_V_Max：第一排温度上升速率的最大值，即提取步骤2.1)所得第一排热电偶温度变化率的最大值，单位为℃/s：

1st_Rising_V_Max＝max(v₍₁₎)

1st_Falling_V_Ave：第一排温度下降速率均值，即计算步骤2.1)所得第一排热电偶温度变化率小于0的温度速率的均值，单位为℃/s；

2nd_Rising_V_Max：第二排温度上升速率的最大值，即提取步骤2.1)所得第二排热电偶温度变化率的最大值，单位为℃/s；

2nd_Rising_V_Max＝max(v₍₂₎)

1st_2nd_Time_Lag：温度上升时滞，即分别标记第二排与第一排温度开始上升时对应的时刻，取二者的时间间隔作为温度上升时滞，单位为s：

1st_2nd_Time_Lag＝t_(2)min-t_(1)min

由此构建特征向量：

s＝[1st_Rising_Amplitude,1st_Rising_V_Max,1st_Falling_V_Ave,2nd_Rising_V_Max,1st_2nd-Time-Lag]

步骤二、提取正常工况特征向量

(1)获取正常工况历史温度数据：任意截取连续25秒的温度数据；

(2)提取和构建第一排、第二排电偶温度特征向量

图3为正常工况下温度及其温度变化率的特征向量提取示意图。由图3可知，特征向量的提取分别以每列电偶为单位，提取正常工况下同列电偶温度沿浇铸方向的变化特征，特征及其提取方法与步骤一(2)相同。

对于正常工况温度，选取30例温度样本。

图4是特征向量集层次聚类和漏钢识别与判定流程图。从图中可以看出，特征向量集层次聚类和漏钢识别与判定主要包含以下步骤：

步骤三、提取在线实时温度特征向量

(1)实时采集和获取结晶器内弧宽面、外弧宽面、左侧窄面、右侧窄面铜板各行、列热电偶当前时刻及之前24秒，共计25秒的温度数据；

(2)提取和构建第一排、第二排电偶温度特征向量：

特征向量的提取分别以每列电偶为单位，提取在线实测时同列电偶温度沿浇铸方向的变化特征，特征及其提取方法与步骤一(2)相同。

步骤四、建立特征向量库

(1)依据步骤一、步骤二、步骤三提取的黏结漏钢、正常工况与在线实测的温度特征建立特征向量样本库D，共计61例；

(2)对特征向量样本库D作归一化处理，得到特征向量集S，并将归一化后的在线实测温度特征向量记为s_new。具体归一化方法如下：

其中，x_ij为特征向量集S中第i个特征向量第j维特征的值，x_jmax、x_jmin分别为所有61特征向量第j维特征的最大值和最小值。

步骤五，特征向量层次聚类

(1)对步骤四得到的特征向量集S实施层次聚类；

1.1)将特征向量集S中的每个向量s看作一个初始类簇C_i＝{s_i}，建立类簇集合C＝{C₁,C₂,...,C₈₁}；其中s_i表示S中第i个向量，C_i表示第i个类簇，i＝1,2,…,61。

1.2)计算和确定类簇集合C中任意两个类簇C_p、C_q的距离：

d(C_p,C_q)＝min(dist(C_pi,C_qj))

其中，C_pi为类簇C_p中第i个特征向量，C_qj为类簇C_q中第j个特征向量，dist(C_pi,C_qj)表示特征向量C_pi、C_qj的欧氏距离；计算类簇C_p、C_q中任意两个向量的距离，取距离的最小值min作为类簇C_p和C_q的距离。

1.3)标记经步骤1.2)计算后距离最小的两个类簇C_m和C_n，将C_m和C_n合并成为一个新的类簇C_{m,n}并加入集合C，同时删除原有的类簇C_m和C_n，经过类簇的添加和删除后，此时集合C中类簇的总数目减1；

1.4)循环执行步骤1.2)～1.3)，当类簇集合C中仅剩两个类簇时，结束循环，完成聚类过程；

(2)检查聚类结果是否满足以下判定条件，即：

所有黏结漏钢特征向量的90％以上属于同一个类簇，且该类簇中正常工况特征向量所占的比率低于20％；

满足此条件则将该类簇记为漏钢类簇C_breakout，另一个类簇记为正常工况类簇C_normal；否则，重新执行步骤(1)、(2)，直至由漏钢、正常工况以及实测温度构成的特征向量集聚类结果满足上述判定条件为止；

步骤六，漏钢识别与报警

判断新特征向量s_new是否属于类簇C_breakout，如果是，则发出漏钢警报；否则，继续执行步骤三、四、五、六。

图5表示在线实测温度1的第一、二排电偶温度图。图中右侧垂直线表示在线检测时的当前时刻，连同该时刻之前的24个时刻，共计25秒的温度数据。对在线温度1进行特征提取后得到的特征向量为：

s_new＝[0.4,0.18,-1.18,1.72,0]，

图6是包含在线实测温度特征向量1，即包含s_new的特征向量集层次聚类及漏钢预报结果图。从图中可以看出，经归一化和层次聚类后，特征向量集聚类为两个类簇：左侧的类簇包含了所有标号为1的黏结漏钢特征向量和5个标号为2的正常工况样本，该类簇中黏结漏钢特征向量所占比率大于黏结漏钢样本总数的90％，正常工况特征向量所占比率小于正常工况样本总数的20％，即6个特征向量，符合类簇判定条件，特征向量集层次聚类成功，所以将该类簇记为漏钢类簇C_breakout；则另一个包含标号2较多的类簇记为正常工况类簇C_normal。从图6中亦可以看出，在线实测温度1经特征提取得到的特征向量s_new，即标号为“N”的样本，在聚类后属于正常工况类簇C_normal，并不属于黏结漏钢类簇C_breakout，因此判断为正常工况，则继续更新温度序列，执行步骤三、四、五、六。

图7表示在线实测温度2的第一、二排电偶温度图。图中右侧垂直线表示在线检测时的当前时刻，连同该时刻之前的24个时刻，共计25秒的温度数据。对在线温度2进行特征提取后得到的特征向量为：

s_new＝[7.4,1.06,-0.29,1.36,12]，

图8是包含在线实测温度特征向量2，即包含s_new的特征向量集层次聚类及漏钢预报结果图。从图中可以看出，经归一化和层次聚类后，特征向量集聚类为两个类簇：左侧的类簇包含了所有标号为1的黏结漏钢样本和5个标号为2的正常工况样本，该类簇中黏结漏钢特征向量所占比率大于黏结漏钢样本总数的90％，正常工况特征向量所占比率小于正常工况样本总数的20％，即6个特征向量，符合类簇判定条件，特征向量集层次聚类成功，所以将该类簇记为漏钢类簇C_breakout；则另一个包含标号2较多的类簇记为正常工况类簇C_normal。从图8中亦可以看出，在线实测温度2经特征提取得到的特征向量s_new，即标号为“B”的样本，在聚类后属于黏结漏钢类簇C_breakout，所以判定为漏钢，发出漏钢警报。

以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。