CN109492335B - 一种退火炉炉温预测方法及系统 - Google Patents

一种退火炉炉温预测方法及系统 Download PDF

Info

Publication number
CN109492335B
CN109492335B CN201811521096.0A CN201811521096A CN109492335B CN 109492335 B CN109492335 B CN 109492335B CN 201811521096 A CN201811521096 A CN 201811521096A CN 109492335 B CN109492335 B CN 109492335B
Authority
CN
China
Prior art keywords
furnace
gas flow
temperature
annealing furnace
data
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.)
Active
Application number
CN201811521096.0A
Other languages
English (en)
Other versions
CN109492335A (zh
Inventor
曹卫华
宋文硕
张永月
吴敏
袁艳
金亚利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences
Original Assignee
China University of Geosciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Geosciences filed Critical China University of Geosciences
Priority to CN201811521096.0A priority Critical patent/CN109492335B/zh
Publication of CN109492335A publication Critical patent/CN109492335A/zh
Application granted granted Critical
Publication of CN109492335B publication Critical patent/CN109492335B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D1/00General methods or devices for heat treatment, e.g. annealing, hardening, quenching or tempering
    • C21D1/26Methods of annealing
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D11/00Process control or regulation for heat treatments
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
    • C21D9/00Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor
    • C21D9/52Heat treatment, e.g. annealing, hardening, quenching or tempering, adapted for particular articles; Furnaces therefor for wires; for strips ; for rods of unlimited length

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Organic Chemistry (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Mechanical Engineering (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Thermal Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Heat Treatment Processes (AREA)

Abstract

本发明提供了一种退火炉炉温预测方法及系统,其方法包括:首先建立基于SVR的退火炉炉温预测模型,然后采用神经网络的十折交叉验证方法和采集的原始生产数据对建立的模型进行训练,最后设定实际的生产参数并结合实时采集的现场生产数据,利用训练好的炉温预测模型对退火炉炉温进行预测。本发明的有益效果是:本发明提出的技术方案可对退火炉实际生产起到很好的预测指导作用,对退火炉板温控制提供有效的控制依据,降低退火炉生产过程中钢卷超温报警比例,改善生产环境,保障生产过程安全稳定。

Description

一种退火炉炉温预测方法及系统
技术领域
本发明涉及冶金领域,尤其涉及一种退火炉炉温预测方法及系统。
背景技术
目前涟钢冷轧热镀锌退火炉采用一级PLC交叉限幅PID控制结合新日铁的二级模型机控制,自投产以来由于环境的多变性、生产工艺过程的复杂性、以及冷轧钢卷品种规格的增加,导致温度控制速度、精度、适应性无法满足实际生产需求。随着智能控制技术与计算机技术的迅速发展,退火炉计算机控制系统的应用日趋广泛,实现了生产过程自动化,退火炉温度控制水平得到显著提高。因此,退火炉中板温控制技术和控制系统的研究受到广泛的关注。
目前,对退火炉的温度控制研究方法从传统的方法到当今流行的机器学习控制方法都有一定涉及:传统PID控制作为经典控制方法,其算法简单,易于实现,具有清晰的控制结构,易于调节。通过设计PID控制器进行退火炉板温控制,该方法简单,但是不能解决频繁的规格变化带来的温度波动问题。在动态性能、抗干扰能力以及稳态精度上,PID控制已经满足不了现代工艺要求,不适用于大惯性和大滞后系统;模糊智能控制系统采用了PC总线式工业控制机,通过原有的以仪表为主要控制装置的串级比值控制系统进行改造,实现了对炉温、燃油、空气流量等的有效调节,但模糊控制精度不够高;专家控制方法是专家系统应用到退火炉温度控制系统中,利用专家系统能够模拟专家思维过程对复杂问题进行推理的特性,推理得出退火炉控制参数,然而专家系统的推理机制和知识获取难度大,推理出的专家知识的准确性不高,调整时间较长;神经网络建立连续退火炉动态数学模型,基于模型实现退火炉温度智能控制,满足了工业过程在线应用的特点,但是神经网络结构中各层神经元数目确定困难,训练开销大,参数调节困难。
由此可见,基于智能控制方法的退火炉温度控制系统还存在一定的缺陷,控制方法复杂,有待进一步的研究。而基于预测控制方法的退火炉温度控制系统能及时有效地抑制干扰对板温的影响,易于在实际生产中应用,且能够有效提高退火炉板温控制精度。因此基于预测控制方法来实现退火炉炉温控制,建立炉温预测模型,从而提高炉温预测精度。
发明内容
为了解决针对冷轧连退生产线中工况切换时多个过程变量与期望板温目标的改变,从而导致分区炉温的波动,导致影响钢卷的性能质量甚至造成长时间超出安全阈值而引发事故的问题,本发明提供了一种退火炉炉温预测方法及系统,一种退火炉炉温预测方法,主要包括以下步骤:
S101:获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;其中,n=1,2,3,…,m,m为退火炉RTF段的分区数量;
S102:对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
S103:建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到训练后的炉温预测模型;
S104:设定退火炉的出炉板温N6和总下发的煤气流量N1,同时获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5
S105:根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温Tn
进一步地,步骤S102中,预处理的步骤包括:
S201:数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
S202:采用范围标准化的方法对数据进行标准化,如公式(1)所示:
Figure GDA0002439768270000031
上式中Xij为原始生产数据;
Figure GDA0002439768270000032
为标准化后的第i个样本的第j个变量;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j的值均大于0。
进一步地,步骤S103中,训练后的炉温预测模型如公式(2)所示:
Figure GDA0002439768270000033
上式中,K(x(i),x(j))为核函数,
Figure GDA0002439768270000034
Figure GDA0002439768270000035
为非线性映射函数,用于将输入样本数据非线性地映射到更高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;约束条件:
Figure GDA0002439768270000036
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;i=1,2,3,…,n;n为样本个数;j>0,为样本i的变量,且为预设值,一个样本内有若干个变量不同的数据点;x()为第i个样本x()中的第j个变量,j=1,2,…,l;l为每个样本中变量的数量;b为偏置,通过KKT条件求得。
进一步地,步骤S103中,采用神经网络中的十折交叉验证方法对退火炉RTF段分区炉温预测模型进行训练。
进一步地,步骤S104中,所述分区煤气流量下发策略具体如下:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值。
进一步地,一种退火炉炉温预测系统,其特征在于:包括以下模块:
原始生产数据获取模块,用于获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;
预处理模块,用于对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
预测模型建立模块,用于建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到训练后的炉温预测模型;
实时生产数据获取模块,用于设定退火炉的出炉板温N6和总下发的煤气流量N1,同时获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5
炉温预测模块,用于根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温Tn
进一步地,预处理模块中,预处理的步骤包括如下单元:
数据筛选单元,用于数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
数据标准化单元,用于采用范围标准化的方法对数据进行标准化,如公式(3)所示:
Figure GDA0002439768270000041
上式中Xij为原始生产数据;
Figure GDA0002439768270000042
为标准化后的第i个样本的第j个生产数据;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j均大于0。
进一步地,预测模型建立模块中,训练后的炉温预测模型如公式(4)所示:
Figure GDA0002439768270000043
上式中,K(x(i),x(j))为核函数,
Figure GDA0002439768270000044
Figure GDA0002439768270000045
为非线性映射函数,用于将输入样本数据非线性地映射到更高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;约束条件:
Figure GDA0002439768270000051
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;i=1,2,3,…,n;n为样本个数;j>0,为样本i的变量,且为预设值,一个样本内有若干个变量不同的数据点;x(j)为第i个样本x(i)中的第j个变量,j=1,2,…,l;l为每个样本中变量的数量;b为偏置,通过KKT条件求得。
进一步地,预测模型建立模块中,采用神经网络中的十折交叉验证方法对退火炉RTF段分区炉温预测模型进行训练。
进一步地,实时生产数据获取模块中,所述分区煤气流量下发策略具体如下:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值。
本发明提供的技术方案带来的有益效果是:本发明提出的技术方案可对退火炉实际生产起到很好的预测指导作用,对退火炉板温控制提供有效的控制依据,降低退火炉生产过程中钢卷超温报警比例,改善生产环境,保障生产过程安全稳定。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例中一种退火炉炉温预测方法的流程图;
图2是本发明实施例中一种退火炉炉温预测系统的模块组成示意图;
图3是本发明实施例中SVR模型和BP模型预测结果对比示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
本发明的实施例提供了一种退火炉炉温预测方法及系统。
请参考图1,图1是本发明实施例中一种退火炉炉温预测方法的流程图,具体包括如下步骤:
S101:获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;其中,n=1,2,3,…,m,m为退火炉RTF段的分区数量;
S102:对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
预处理的步骤包括:
S201:数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
S202:采用范围标准化的方法对数据进行标准化,如公式(1)所示:
Figure GDA0002439768270000061
上式中Xij为原始生产数据;
Figure GDA0002439768270000062
为标准化后的第i个样本的第j个生产数据;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j均大于0;
S103:建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到如公式(2)所示的训练后的炉温预测模型:
Figure GDA0002439768270000063
上式中,K(x(i),x(j))为核函数,
Figure GDA0002439768270000064
Figure GDA0002439768270000065
为非线性映射函数,将输入样本数据非线性地映射到高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;s.t.
Figure GDA0002439768270000066
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;i=1,2,3,…,n;n为样本个数;j>0,为样本i的变量,一个样本内有若干个变量不同的数据点;l为每个样本中的变量总数量;b为偏置,通过KKT(Karush-Kuhn-Tucker)条件求得;
在本发明实施例中,针对涟钢的某退火炉进行实际炉温预测模型的具体建立及训练过程如下:
应用SVR算法构建由主要影响参数到最终分区炉温间的关系模型,即Tn=Fn(dn、N1、N2、N3、N4、N5、N6),其中n=1,…,7。以RTF段1区炉温预测模型为例,T1=F1(d1、N1、N2、N3、N4、N5、N6),根据支持向量机的基本原理,可以将该模型关系描述为,给定要辨识系统的7维空间的样本数据集合,如公式(3)所示:
Figure GDA0002439768270000071
非线性函数映射
Figure GDA0002439768270000072
R7→R将7维空间的输入样本数据非线性地映射到更高维特征空间。那么就在高维特征空间里构造了最优函数
Figure GDA0002439768270000073
Figure GDA0002439768270000074
这样将输入空间里的非线性函数估计问题转化到高维特征空间中的线性函数估计问题。为了获得在特征空间里的线性函数需要求解下式优化问题,该优化问题是为了获得一个结构风险最小化,如公式(4)所示:
Figure GDA0002439768270000075
上式中,C为平衡因子,控制对超出误差的样本的惩罚程度。支持向量回归机的损失函数是引入一个ε不敏感损失函数,如公式(5)所示:
Figure GDA0002439768270000076
上式中,ε为不敏感损失函数参数,它对于误差小于ε的输入量不进行惩罚,ε取值大小影响支持向量的数目。当ε取值较大时支持向量的数目下降,模型精度较低;当ε取值较小时支持向量的数目会增多,同时模型精度会提高。支持向量机通过采用ε不敏感损失函数在高维特征空间完成线性回归,同时通过最小化||ω||2来减小模型的复杂度,在考虑到允许拟合误差的情况,也引入了非负的松弛变量ξ,ξ*≥0。通过引入为Lagrange函数及Lagrange乘子α(i),α*(i),η(i),η*(i),可以将式(4)在约束条件下的函数求极值问题转换为其对应的对偶优化问题,如公式(6)所示:
Figure GDA0002439768270000081
上式中
Figure GDA0002439768270000082
径向基函数K(x(i),x(j))为核函数。
在求解式(6)之后,可以获得如公式(7)所示的回归函数:
Figure GDA0002439768270000083
上式中,b为偏置,通过KKT(Karush-Kuhn-Tucker)条件得偏置b为:
Figure GDA0002439768270000084
式(7)即为退火炉RTF段1区(整个RTF段共有7个区)的炉温预测模型,同理,也可以得到其他6个区的炉温预测模型,其他6个区的炉温预测模型与式(7)一致。
S104:设定退火炉的出炉板温N6和总下发的煤气流量N1,同时获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5;其中,板厚N2、板宽N3、和生产速度N4,通过工业PLC获取;实际进炉板温N5通过安装于退火炉入口处的红外测温仪测量获得;
所述分区煤气流量下发策略,具体为:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值。
本实施例针对涟钢某退火炉具体生产工艺进行如下举例说明(以4个区为例):
增加煤气流量时:按照4区→3区→2区→1区的顺序将煤气流量增加到本区不超过炉温上限值的最高煤气流量。即不超过4区炉温上限1280℃的情况下可以将煤气流量最高加到850Nm3/h,如果炉温即将超限或850Nm3/h已加满,则考虑将煤气流量加到3区,若3区即将超过炉温上限1280℃或650Nm3/h已满,则进行2区操作,最后进行1区操作;
减少煤气流量时,按照1区→2区→3区→4区的顺序减少至本区下限。首先从1区减少至煤气流量值200Nm3/h,再依次减少2区,3区和4区。
数据统计中存在1200-1300Nm3/h总煤气流量时,同时打开4个区,123区的煤气流量下限值为150Nm3/h,但是再细化总煤气流量区间为四个时会造成操作工人的操作困难,所以都统一到最低限为200Nm3/h,当煤气流量低于1300Nm3/h时开始关区。
S105:根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温Tn
请参阅图2,图2是本发明实施例中一种退火炉炉温预测系统的模块组成示意图,其特征在于,该系统包括顺次连接的原始生产数据获取模块11、预处理模块12、预测模型建立模块13、实时生产数据获取模块14和炉温预测模块15;
原始生产数据获取模块11,用于获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;
预处理模块12,用于对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
预测模型建立模块13,用于建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到如公式(8)所示的训练后的炉温预测模型:
Figure GDA0002439768270000091
上式中,K(x(i),x(j))为核函数,
Figure GDA0002439768270000101
Figure GDA0002439768270000102
为非线性映射函数,将输入样本数据非线性地映射到更高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;约束条件:
Figure GDA0002439768270000103
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;n=1,…,7,代表退火炉RTF段七个分区炉温的编号;i=1,2,3,…,n;n为样本个数;j>0,为变量,且为预设值,一个样本内有若干个变量不同的数据点;l为每个样本中的变量总数量;b为偏置,通过KKT(Karush-Kuhn-Tucker)条件求得;
实时生产数据获取模块14,用于设定退火炉的出炉板温N6和总下发的煤气流量N1,获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5;其中,板厚N2、板宽N3、和生产速度N4,通过工业PLC获取;实际进炉板温N5通过安装于退火炉入口处的红外测温仪测量获得;
炉温预测模块15,用于根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温。
在本实施例中,预处理的步骤包括如下单元:
数据筛选单元,用于数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
数据标准化单元,用于采用范围标准化的方法对数据进行标准化,如公式(9)所示:
Figure GDA0002439768270000111
上式中Xij为原始生产数据;
Figure GDA0002439768270000112
为标准化后的第i个样本的第j个生产数据;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j均大于0;
在本实施例中,采用神经网络中的十折交叉验证方法对退火炉RTF段分区炉温预测模型进行训练。
在本实施例中,实时生产数据获取模块中,所述分区煤气流量下发策略具体如下:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值。
在本实施例中,基于机理分析的方法,选择出影响退火炉炉温的关键因素,根据关键因素获取对应的原始生产数据,并根据关键因素,确定SVR炉温预测模型的输入。以下是机理分析的具体过程:
煤气燃烧过程必须遵循能量守恒定律,即煤气燃烧释放的热量Qg与板温变化产生的热量和炉体温度变化产生的热量Qf的总量守恒,如式(10)所示:
Qg=Qf+Qc (10)
能够直接反映炉体和钢板热量变化的参数即为由热电偶测的炉温和红外测温仪测得的板温。由于炉体为复杂的工业设备,对应吸热能力没有准确的计算方法,所以Qf的值很难直接得到,但是根据煤气燃烧热量计算公式和固体热量变化原理,Qg和Qc的计算如公式(11):
Figure GDA0002439768270000113
式中q为煤气热值且一般为常数;N1为煤气燃烧的体积;c为钢铁的比热容;m为炉中加热钢板的质量;ρ为钢铁的密度;N2和N3分别为钢板的厚度和宽度;N4t表示钢板在N4的生产速度下t时间内通过的长度;ΔTc则表示为钢板进入RTF段前后的板温变化量。
因为各区炉温Tn与炉体总吸热Qf及分区炉体吸热Qf_n直接相关,其中n=1,…,7。结合式(10)和(11),可分析得知:影响生产状况下Tn的变量包括N1、N2、N3、N4和ΔTc,而ΔTc则需要使用进炉板温N5和设定出炉板温N6计算,同时各分区所分配的煤气流量dn也会对Tn产生重要影响,所以和单分区炉温相关的变量因素主要有7个,可以据此构建RTF段分区炉温预测模型。
根据模型的输入输出结构,从现场采集了钢板在稳定工况生产时对应母钢卷的信息和实时生产的数据,再经过数据预处理后得到有效数据共413组,训练集数据包含300组,测试集数据为113组。
在使用SVR方法进行建模的同时,也利用工业现场建模广泛使用的BP神经网络方法进行对比验证,两种算法使用相同的数据集进行训练和测试,以其构建的RTF段7区炉温预测模型为例,模型在测试集的验证结果如图3所示。
退火炉RTF段炉温稳定范围较小,所以对预测精度要求较高,在现场一般要求预测误差在±5℃之内的数据占比达到80%以上;同时最大偏差不能超过20℃,否则就认为模型无效。据此标准可以对SVR模型和BP神经网络模型的预测结果进行评估,其统计结果如表1所示。
表1 SVR模型和BP模型统计结果对比
Figure GDA0002439768270000121
由表1可知,SVR模型在预测精度和最大偏差均优于BP神经网络模型且满足现场要求,所以可以投入现场使用对分区煤气流量进行多目标优化。
本发明的有益效果是:本发明提出的技术方案可对退火炉实际生产起到很好的预测指导作用,对退火炉板温控制提供有效的控制依据,降低退火炉生产过程中钢卷超温报警比例,改善生产环境,保障生产过程安全稳定。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种退火炉炉温预测方法,其特征在于:包括以下步骤:
S101:获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;其中,n=1,2,3,…,m,m为退火炉RTF段的分区数量;
S102:对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
S103:建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到训练后的炉温预测模型;
S104:设定退火炉的出炉板温N6和总下发的煤气流量N1,同时获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5
S105:根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温Tn;其中,所述的总下发的煤气流量等于所述的煤气燃烧的体积;
步骤S104中,所述分区煤气流量下发策略具体如下:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值;
其中,在减少煤气流量时,各区的最低煤气流量值统一为200Nm3/h,当煤气流量低于1300Nm3/h时开始关区。
2.如权利要求1所述的一种退火炉炉温预测方法,其特征在于:步骤S102中,预处理的步骤包括:
S201:数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
S202:采用范围标准化的方法对数据进行标准化,如公式(1)所示:
Figure FDA0002615920120000021
上式中Xij为原始生产数据;
Figure FDA0002615920120000022
为标准化后的第i个样本的第j个变量;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j的值均大于0。
3.如权利要求1所述的一种退火炉炉温预测方法,其特征在于:步骤S103中,训练后的炉温预测模型如公式(2)所示:
Figure FDA0002615920120000023
上式中,K(x(i),x(j))为核函数,
Figure FDA0002615920120000024
Figure FDA0002615920120000025
为非线性映射函数,用于将输入样本数据非线性地映射到更高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;约束条件:
Figure FDA0002615920120000026
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;i=1,2,3,…,n;n为样本个数;j>0,为样本i的变量,且为预设值,一个样本内有若干个变量不同的数据点;x(j)为第i个样本x(i)中的第j个变量,j=1,2,…,l;l为每个样本中变量的数量;b为偏置,通过KKT条件求得。
4.如权利要求1所述的一种退火炉炉温预测方法,其特征在于:步骤S103中,采用神经网络中的十折交叉验证方法对退火炉RTF段分区炉温预测模型进行训练。
5.一种退火炉炉温预测系统,其特征在于:包括以下模块:
原始生产数据获取模块,用于获取退火炉原始生产数据;所述原始生产数据包括:煤气燃烧的体积N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、设定的出炉板温N6、分配的煤气流量dn和实际炉温;
预处理模块,用于对获取的退火炉原始生产数据进行预处理,得到预处理后符合SVR模型的生产数据;
预测模型建立模块,用于建立基于SVR的退火炉RTF段分区炉温预测模型,并利用预处理后的符合SVR模型的生产数据对模型进行训练,得到训练后的炉温预测模型;
实时生产数据获取模块,用于设定退火炉的出炉板温N6和总下发的煤气流量N1,同时获取退火炉生产过程中的实时生产数据;并根据总下发的煤气流量N1和分区煤气流量下发策略,计算得到各分区分配的煤气流量dn;所述实时生产数据包括:板厚N2、板宽N3、生产速度N4、实际进炉板温N5
炉温预测模块,用于根据总下发的煤气流量N1、板厚N2、板宽N3、生产速度N4、实际进炉板温N5、出炉板温N6和各分区分配的煤气流量dn,采用训练后的炉温预测模型对退火炉的分区炉温进行预测,得到退火炉RTF段的各分区预测炉温Tn
实时生产数据获取模块中,所述分区煤气流量下发策略具体如下:
增加煤气流量时:按照退火炉各分区编号从高到低的顺序将煤气流量依次增加到各区不超过对应区炉温上限值的最高煤气流量值;
减少煤气流量时:按照退火炉各分区编号从低到高的顺序将煤气流量依次减少到各区不低于对应区炉温下限值的最低煤气流量值;
其中,在减少煤气流量时,各区的最低煤气流量值统一为200Nm3/h,当煤气流量低于1300Nm3/h时开始关区。
6.如权利要求5所述的一种退火炉炉温预测系统,其特征在于:预处理模块中,预处理的步骤包括如下单元:
数据筛选单元,用于数据筛选:采用限幅滤波的方法,剔除不符合工艺条件和异常工况的数据;
数据标准化单元,用于采用范围标准化的方法对数据进行标准化,如公式(3)所示:
Figure FDA0002615920120000041
上式中Xij为原始生产数据;
Figure FDA0002615920120000042
为标准化后的第i个样本的第j个生产数据;Xjmin和Xjmax分别表示样本集中变量j的最小和最大值;i和j均大于0。
7.如权利要求5所述的一种退火炉炉温预测系统,其特征在于:预测模型建立模块中,训练后的炉温预测模型如公式(4)所示:
Figure FDA0002615920120000043
上式中,K(x(i),x(j))为核函数,
Figure FDA0002615920120000044
Figure FDA0002615920120000045
为非线性映射函数,用于将输入样本数据非线性地映射到更高维特征空间;α(j)和α*(j)分别为Lagrange函数及Lagrange乘子;约束条件:
Figure FDA0002615920120000046
C为平衡因子,控制对超出误差的样本的惩罚程度,为预设值;x(i)=[dn(i),N1(i),N2(i),N3(i),N4(i),N5(i),N6(i)]T;i=1,2,3,…,n;n为样本个数;j>0,为样本i的变量,且为预设值,一个样本内有若干个变量不同的数据点;x(j)为第i个样本x(i)中的第j个变量,j=1,2,…,l;l为每个样本中变量的数量;b为偏置,通过KKT条件求得。
8.如权利要求5所述的一种退火炉炉温预测系统,其特征在于:预测模型建立模块中,采用神经网络中的十折交叉验证方法对退火炉RTF段分区炉温预测模型进行训练。
CN201811521096.0A 2018-12-12 2018-12-12 一种退火炉炉温预测方法及系统 Active CN109492335B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811521096.0A CN109492335B (zh) 2018-12-12 2018-12-12 一种退火炉炉温预测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811521096.0A CN109492335B (zh) 2018-12-12 2018-12-12 一种退火炉炉温预测方法及系统

Publications (2)

Publication Number Publication Date
CN109492335A CN109492335A (zh) 2019-03-19
CN109492335B true CN109492335B (zh) 2020-12-08

Family

ID=65710007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811521096.0A Active CN109492335B (zh) 2018-12-12 2018-12-12 一种退火炉炉温预测方法及系统

Country Status (1)

Country Link
CN (1) CN109492335B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110874511B (zh) * 2019-11-19 2022-09-16 武汉晨曦芸峰科技有限公司 一种智能炉料配比计算方法
CN112430727B (zh) * 2020-10-15 2022-05-20 宝钢日铁汽车板有限公司 一种连续退火炉炉温预警方法及系统
CN114058833B (zh) * 2021-11-17 2023-12-22 上海途开者机械工程技术有限公司 连续退火炉自动控制方法和系统、计算机可读存储介质以及计算机程序产品
CN115930594B (zh) * 2023-02-20 2023-05-30 科大智能物联技术股份有限公司 一种竖炉生产控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1149082A (zh) * 1996-08-27 1997-05-07 宝山钢铁(集团)公司 连续退火炉在线控制方法
CN101139652A (zh) * 2007-10-26 2008-03-12 武汉钢铁(集团)公司 罩式炉钢卷退火工艺中的离线预测方法
CN101456038A (zh) * 2009-01-08 2009-06-17 上海交通大学 热轧带钢层流冷却过程板带温度监测方法
JP4869670B2 (ja) * 2005-10-04 2012-02-08 新日本製鐵株式会社 連続焼鈍炉における板温制御方法及び装置
CN105886751A (zh) * 2016-04-13 2016-08-24 中国地质大学(武汉) 一种冷轧热镀锌退火炉板温协调控制系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AT519995B1 (de) * 2017-05-29 2021-04-15 Andritz Ag Maschf Verfahren zur Regelung der Aufwickeltemperatur eines Metallbandes

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1149082A (zh) * 1996-08-27 1997-05-07 宝山钢铁(集团)公司 连续退火炉在线控制方法
JP4869670B2 (ja) * 2005-10-04 2012-02-08 新日本製鐵株式会社 連続焼鈍炉における板温制御方法及び装置
CN101139652A (zh) * 2007-10-26 2008-03-12 武汉钢铁(集团)公司 罩式炉钢卷退火工艺中的离线预测方法
CN101456038A (zh) * 2009-01-08 2009-06-17 上海交通大学 热轧带钢层流冷却过程板带温度监测方法
CN105886751A (zh) * 2016-04-13 2016-08-24 中国地质大学(武汉) 一种冷轧热镀锌退火炉板温协调控制系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A Dynamic Data-driven Model for Predicting Strip Temperature in Continuous Annealing Line Heating Process;曹卫华等;《第37届中国控制会议论文集》;20180725;第1887-1891页 *
基于支持向量机的退火炉煤气流量预测模型;曹卫华 等;《全国冶金自动化信息网2017年会论文集》;20171025;第38-41页 *

Also Published As

Publication number Publication date
CN109492335A (zh) 2019-03-19

Similar Documents

Publication Publication Date Title
CN109492335B (zh) 一种退火炉炉温预测方法及系统
CN109253494B (zh) 一种基于热负荷预测的电蓄热装置的控制方法
CN111237989B (zh) 一种基于负荷预测的楼宇通风空调控制方法及装置
Li et al. Modeling for indoor temperature prediction based on time-delay and Elman neural network in air conditioning system
CN109062053A (zh) 一种基于多变量校正的脱硝喷氨控制方法
CN101509812B (zh) 一种冶金加热炉炉内钢坯温度分布检测的软测量方法
CN102540879A (zh) 基于群决策检索策略的多目标评价优化方法
Li et al. Coke oven intelligent integrated control system
CN113849020B (zh) 一种基于人工智能算法的钢坯升温曲线设计方法及装置
CN112330012B (zh) 一种基于迁移学习的建筑能耗预测方法及设备
Wu et al. Integrated soft sensing of coke-oven temperature
CN106200379A (zh) 一种非自衡对象的分布式动态矩阵控制方法
CN111931436A (zh) 一种基于数值模拟与神经网络的燃烧器喷口风量预测方法
CN101126917A (zh) 槽式反应器基于多核支持向量机的非线性预测控制方法
CN105955014A (zh) 分布式动态矩阵控制优化的焦炭炉炉膛压力控制方法
CN115608793A (zh) 一种机理融合数据的精轧温度调控方法
CN103529699A (zh) 一种煤气化炉系统的炉温自学习控制方法
Zhao et al. A model predictive control regulation model for radiant air conditioning system based on delay time
CN105400951B (zh) 一种用于硅钢退火的加湿混合控制系统及其控制方法
CN113887833A (zh) 一种分布式能源用户侧逐时负荷预测方法及系统
Ji et al. Application of the improved the ELM algorithm for prediction of blast furnace gas utilization rate
CN109598643B (zh) 一种退火炉分区煤气流量分配方法、设备及存储设备
CN111336828A (zh) 一种基于fcm模糊时间序列的加热炉炉温控制器
Tian Predictive control of coke oven flue temperature based on orthogonal neural network
CN112036091A (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
GR01 Patent grant