CN103882176A - 一种基于数据驱动的转炉炼钢过程在线动态最优控制方法 - Google Patents
一种基于数据驱动的转炉炼钢过程在线动态最优控制方法 Download PDFInfo
- Publication number
- CN103882176A CN103882176A CN201410114943.7A CN201410114943A CN103882176A CN 103882176 A CN103882176 A CN 103882176A CN 201410114943 A CN201410114943 A CN 201410114943A CN 103882176 A CN103882176 A CN 103882176A
- Authority
- CN
- China
- Prior art keywords
- carbon element
- element content
- temperature
- liquid steel
- steel temperature
- 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.)
- Granted
Links
Images
Landscapes
- Carbon Steel Or Casting Steel Manufacturing (AREA)
Abstract
一种基于数据驱动的转炉炼钢过程在线动态最优控制方法,包括建立离线预测模型数据库;利用数据驱动的方法建立钢水温度预测模型和碳元素含量预测模型,并得到操作变量与转炉炼钢过程中钢水温度与钢水碳含量之间的对应关系;选择与当前转炉炼钢生产过程信息匹配的炉次数据集合,确定钢水温度与碳元素含量的控制参考曲线;建立转炉炼钢过程实时动态优化模型,确定各操作变量优化设定值集合;从中选择操作变量设定值进行控制操作。本发明能实现对转炉炼钢过程实时在线控制,方便操作人员根据实际工况设定选择条件,提高炼钢厂生产效率。
Description
技术领域
本发明属于冶金行业过程控制技术领域,涉及一种基于数据驱动的转炉炼钢过程在线动态最优控制方法。
背景技术
转炉炼钢的产量占目前炼钢总产量的80%以上,转炉炼钢过程钢的控制是保证冶炼钢水质量的一个重要条件。
转炉炼钢生产工艺和主要控制流程如图1所示。一定量的固体废钢和经过脱硫处理的液态铁水装入转炉炉体中。然后,垂直炉底的氧枪从上方进入炉口,氧气F011从枪口吹入,同时添加石灰石、白云石等副原料F001-F007,并在炉底吹入氮气和惰性气体F012进行搅拌,热金属与氧气接触的表面形成熔池,需要去除的碳、硅、锰等元素与氧发生反应,生成物进入炉渣或直接从烟气中排出,同时产生的能量使钢水温度升高。在吹炼到达吹炼后期时下取样得到当前钢水的温度与成分信息,根据这些信息与目标终点值的偏差计算得出还需要加入的副原料量、顶吹氧气量、底吹气体量及氧***式,进行动态阶段吹炼。加入的副原料量、顶吹氧气量、底吹气体量及氧***式即转炉炼钢操作变量。当动态阶段吹炼结束时再次取样,若温度与成分满足要求则进行合金化处理或者出钢,若不满足则继续补吹直至达到要求为止,如图8所示。
虽然操作人员已经在生产作业的控制方面积累了相当多的经验,也取得了很大的成绩,但是由于转炉生产中存在的众多生产工艺约束,而且转炉冶炼过程中高温、多种物理化学变化反应速度快、冶炼周期短等特点,生产控制非常复杂。当前的生产控制方法主要存在以下问题:
转炉炼钢生产中冶炼过程环境复杂导致很多数据无法精确测量。炼钢过程是一个非常复杂的多元多相高温状态下进行的物理化学过程,存在很多难以定量的非线性因素,并且不易获得准确实时的检测信息,所以模型的建立与过程的控制都是极为困难的。自动副枪技术是利用副枪的移动携带测量探头深入钢水中直接取样,该方法虽然测量精确,但是成本较高测量时破坏炉内冶炼环境,无法多次连续测量;声纳仪的应用是利用发射超声与接收炉腔的回声来判断炉内的冶炼状态,光谱分析仪是通过提取炉口火焰光强与图像信息来判断炉内冶炼状态,由于炉内环境复杂,化学反应剧烈,浓烟与小型喷溅的干扰直接影响炉内和炉口信息的准确性;烟气分析方法作为转炉炼钢控制的主要方法之一,是通过连续检测转炉吹炼过程中产生的烟气成分和烟气流量,来推算出炉内冶炼状态,由于检测设备位于炉顶烟道内,取样时间大滞后于反应发生的时间。
转炉炼钢生产中已有的模型与控制方法大多未考虑确定冶炼过程中的温度及成分含量的最优轨迹曲线并利用这一结果进行控制。现有的文献大多是研究转炉终点控制,在冶炼过程中的利用温度与碳等元素含量的反应走势进行控制的研究方面报道的较少。中国专利CN202401090U是利用烟气反馈的数据根据模型计算脱碳氧效率值来推定当前钢水中的碳浓度与温度来判断终点时间。当烟气数据产生异常波动的情况时模型结果会产生较大偏差。新日铁与西门子炼钢***中的温度与成分中各元素的含量走势也是以预报为目的,没有明确考虑在当前冶炼条件下温度与成分含量的走势并利用这一结果进行冶炼控制。
发明内容
针对现有技术存在的问题,本发明提供一种基于数据驱动的转炉炼钢过程在线动态最优控制方法。
本发明的技术方案是这样实现的:
一种基于数据驱动的转炉炼钢过程在线动态最优控制方法,包括以下步骤:
步骤1:建立基于数据驱动的离线预测模型数据库;
步骤1.1:建立动态优化精英炉次集合;
步骤1.1.1:在炼钢过程的吹炼后期至出钢结束,每间隔一个时间周期投掷一次碳温氧探头,测量炼钢转炉中的钢水温度和碳元素含量信息;
步骤1.1.2:将碳温氧探头测量的钢水温度与光谱分析仪测得的转炉炉口火焰温度进行比对,将碳温氧探头测量的钢水中碳元素含量与烟气分析仪得到的碳元素含量进行比对,得到炼钢转炉内的钢水实际温度与实际碳元素含量;
步骤1.1.3:选出钢水实际温度、碳元素含量与目标值相符的炉况稳定的炉次建立动态优化精英炉次数据集合,该集合中每一炉次的信息包括:测量范围内的各个时刻的钢水温度、钢水碳元素含量、顶吹氧气模式、不同时间范围内加入的氧气总量、底吹气体总量、各副原料加入的重量、烟气量、烟气中一氧化碳量、烟气中二氧化碳量、入炉铁水重量和入炉废钢重量;
步骤1.2:对动态优化精英炉次数据集合分别以钢水温度和钢水碳元素含量为预测目标按照各自属性进行分类;
步骤1.2.1:将动态优化精英炉次数据集合中的炉次按照钢种进行预划分;
步骤1.2.2:将动态优化精英炉次数据集合中同一钢种的炉次,分别以钢水温度与钢水碳元素含量为预测目标进行聚类,其中,以钢水温度为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水温度、主原料加入量四种属性划分的多个炉次数据集合,以钢水碳元素含量为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水碳元素含量、主原料加入量四种属性划分的多个炉次数据集合;
步骤2:采用最小二乘支持向量机方法,建立钢水温度预测模型和碳元素含量预测模型,钢水温度预测模型的输入和碳元素含量预测模型的输入均为当前时刻钢水温度、当前时刻钢水碳元素含量、某一时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量、各副原料加入的重量、烟气量、烟气中一氧化碳量、烟气中二氧化碳量、入炉铁水重量和入炉废钢重量,钢水温度预测模型的输出为下一时刻预测的钢水温度,碳元素含量预测模型的输出为下一时刻预测的钢水碳元素含量;
步骤3:在吹炼后期下副枪取样得到当前时刻钢水温度与钢水中各元素成分信息,得到当前时刻的钢水温度和钢水碳元素含量:
步骤4:根据钢水温度预测模型和碳元素含量预测模型确定转炉炼钢过程中的钢水温度和钢水碳元素含量;
步骤4.1:采集当前生产过程中的钢水温度与钢水碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量信息;
步骤4.2:选择与步骤4.1采集的当前生产过程中的信息匹配的步骤1.2中划分的炉次数据集合,并利用步骤2建立的钢水温度预测模型和碳元素含量预测模型对钢水温度与钢水碳元素含量进行预测;
步骤5;确定生产过程中各时刻钢水温度与碳元素含量,作为参考曲线;
步骤5.1:对动态优化精英炉次数据集合中同一钢种的炉次,以钢水温度和钢水碳元素含量为预测目标,分别按照当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量5个属性进行分类;
步骤5.2:选择与生产过程中当前时刻的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量相匹配的步骤5.1中划分的炉次数据集合作为参考炉次集合,并从中选取最相似的h炉作为参考炉次;
步骤5.3:计算各参考炉次与当前炉次的钢水温度、碳元素含量和入炉主原料加入量的不同产生的偏差数;
步骤5.4:确定各参考炉次的权重系数;
步骤5.5:确定当前炉次的未来各时刻的钢水温度与碳元素含量,作为参考曲线;
步骤6:根据烟气在线分析仪采集的转炉炼钢过程各个时刻的烟气流量、CO含量及CO2含量值与对应时刻的参考曲线中的烟气流量、CO含量及CO2含量值进行比较得到偏差值,实时校正钢水温度预测模型和碳元素含量预测模型;
步骤7:建立转炉炼钢过程实时动态优化模型,该模型以最小化当前预测的钢水温度和碳元素含量值与参考曲线的偏差为目标,该模型操作变量是当前时刻到下一时刻的时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量和各副原料加入的重量;
步骤8:确定转炉炼钢过程实时动态优化模型的各操作变量的优化设定值集合;
步骤9:根据转炉炼钢生产实际,从操作变量的优化设定值集合中选择操作变量设定值,根据操作变量设定值对转炉生产进行操作;
步骤10:若预测得到的钢水温度值与参考曲线中的参考值误差大于1.8%,或者钢水碳元素含量值与参考曲线中的参考值误差大于2.3%,则返回步骤5.1;否则,返回步骤5.2,进行转炉炼钢过程的实时动态优化与控制,直到出钢结束;
步骤11:定期对离线预测模型数据库进行更新。
所述的基于数据驱动的转炉炼钢过程在线动态最优控制方法采用的转炉炼钢过程在线动态控制***,包括炼钢转炉、氧枪及副枪,还包括光谱分析仪、烟气分析仪和投掷探头;光谱分析仪安装在炼钢转炉炉口,烟气分析仪安装在转炉烟道处,投掷探头采用碳温氧探头。
有益效果:
本发明提出了一种基于数据驱动的转炉的实时在线动态最优控制方法,在该方法中给出了转炉炼钢生产动态阶段各个时刻的操作变量的设定值,按照该设定值进行控制能够实现对转炉炼钢过程的实时在线控制,直接提高钢水温度和碳元素含量的命中率;本发明提出的最优控制方法在生产过程中实时控制钢水温度与碳元素含量,使钢水温度平均终点命中偏差小于3.1度,钢水碳元素含量的平均命中精度到达0.97%,该方法能够提供给现场操作人员多种操作设定与控制方案,方便操作人员根据实际工况需求灵活设定选择条件,进而帮助提高炼钢厂的整体生产效率。
附图说明
图1为本发明所涉及的具体转炉生产工艺和控制流程示意图;
图2为本发明具体实施方式的转炉炼钢过程在线动态最优控制方法原理图;
图3(a)为本发明具体实施方式的温度的预测值与实际值比较的结果;
图3(b)为本发明具体实施方式的钢水碳元素含量预测与实际比较误差值;
图4(a)为本发明具体实施方式的碳元素含量的实际出钢结果与目标值比较误差值;
图4(b)为本发明具体实施方式钢水温度的实际出钢结果与目标值比较结果图;
图5(a)为本发明具体实施方式计算得到的当前实际顶吹氧气流量的最优设定值与控制器输出值的比较结果图;
图5(b)为本发明具体实施方式计算得到的当前底吹气体流量的最优设定值与控制器输出值的比较结果图;
图5(c)为本发明具体实施方式计算得到的氧***式的最优设定值与控制器输出值的比较结果图;
图6(a)为本发明具体实施方式计算得到的碳元素含量值与控制参考值的比较结果图;
图6(b)为本发明具体实施方式计算得到的实际钢水温度与控制参考值的比较结果图;
图7(a)为本发明具体实施方式所需的当前实际烟气流量检测值与参考值的比较结果图;
图7(b)为本发明具体实施方式所需的当前实际CO、CO2在烟气中的含量检测值与参
考值的比较结果图;
图8为本发明所设定的转炉炼钢生产过程阶段示意图;
图9为本发明具体实施方式的基于数据驱动的转炉炼钢过程在线动态最优控制方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式做详细说明。
本实施方式的基于数据驱动的转炉炼钢过程在线动态最优控制方法采用的转炉炼钢过程在线动态控制***,包括炼钢转炉、氧枪及副枪、安装在炼钢转炉炉口的光谱分析仪和烟气分析仪,副枪采用碳温氧探头。
转炉生产工艺和控制流程如图1所示。图中FIC001-FIC007分别为加入块状石灰石、轻烧白工石、菱镁球等副原料的重量控制器;FIC011和FIC012分别为顶吹氧气的气流量控制器和底吹气体的气流量控制器;AI023为转炉副枪取样的在线分析仪,用来获取取样时刻钢水温度和碳元素含量的实际测量值;DIC013为氧***式高度控制器;FI021和AI022分别为烟气流量检测器和烟气成分在线分析仪,用来获取各个时刻烟气的流量值和烟气中CO和CO2的含量值;OIC010为转炉炼钢过程实时动态优化模型控制器。在实施过程中从动态阶段开始的t时刻,烟气在线分析仪AI023进行采样,OIC010根据烟气流量检测器FI021和烟气在线分析仪AI022的分析结果对预测模型进行校正,并得到相应的参考曲线,然后基于转炉炼钢过程实时动态优化模型,获得t时刻各控制变量的最优设定值,该设定值由底层的控制***进行实施。这些设定值将保持不变,直到下一个采样时刻t+1;然后将预测结果代替烟气在线分析仪AI023得到的钢水温度和碳元素含量的实际测量值进行下一步的操作。
本实施方式采用数据解析技术的基于数据驱动的转炉炼钢过程在线动态最优控制方法,该方法以副原料石灰石、白云石的加入时间与重量、各时刻氧气吹入量和底吹气体吹入量、氧枪控制模式为控制变量,原理如图2所示:首先,采集大量的实际炼钢生产过程中的钢水温度与钢水碳元素含量、入炉主原料加入量、目标出钢温度、目标出钢碳元素含量信息,并建立动态阶段参考炉次信息库,用数据解析的方法建立炉温与钢水碳元素含量的预报模型。然后,根据动态阶段得到的当前冶炼温度与钢水碳元素含量信息,结合目标出钢温度与目标碳含量,利用工艺规则和聚类参考炉次信息等方法来确定出最佳的吹炼温度(钢水温度)和钢水碳元素含量的走势曲线,该曲线作为本炉次控制操作的参考曲线。根据烟气在线分析仪采集的转炉炼钢过程各个时刻的烟气流量、CO含量及CO2含量值与对应时刻的参考曲线中的烟气流量、CO含量及CO2含量值进行比较得到偏差值,实时校正钢水温度预测模型和碳元素含量预测模型;建立转炉炼钢过程实时动态优化模型,确定转炉炼钢过程实时动态优化模型的各操作变量的优化设定值集合;根据转炉炼钢生产实际,从操作变量的优化设定值集合中选择操作变量设定值,根据操作变量设定值对转炉生产进行操作。本发明的最优控制方法可以直观的全面地掌握转炉炼钢现场生产情况。
本实施方式的基于数据驱动的转炉炼钢过程在线动态最优控制方法,如图9所示,包括以下步骤:
步骤1:建立基于数据驱动的离线预测数据库;
步骤1.1:建立动态优化精英炉次数据集合;
步骤1.1.1:在炼钢过程的吹炼后期至出钢结束,每间隔一个时间周期投掷一次碳温氧探头,测量炼钢转炉中的钢水温度和碳元素含量信息;
步骤1.1.2:将碳温氧探头测量的钢水温度与光谱分析仪测得的转炉炉口火焰温度进行比对,将碳温氧探头测量的钢水中碳元素含量与烟气分析仪得到的碳元素含量进行比对,得到炼钢转炉内的钢水实际温度与实际碳元素含量;
步骤1.1.3:选出钢水实际温度、碳元素含量与目标值相符的炉况稳定的炉次建立动态优化精英炉次数据集合,该集合中每一炉次的信息包括:测量范围内的各个时刻的钢水温度、钢水碳元素含量、顶吹氧气模式、不同时间范围内加入的氧气总量、底吹气体总量、各副原料加入的重量、烟气量、烟气中一氧化碳量、烟气中二氧化碳量、入炉铁水重量和入炉废钢重量;
步骤1.2:对动态优化精英炉次数据集合分别以钢水温度和钢水碳元素含量为预测目标按照各自属性进行分类;
步骤1.2.1:将动态优化精英炉次数据集合中的炉次按照钢种进行预划分;
步骤1.2.2:将动态优化精英炉次数据集合中同一钢种的炉次,分别以钢水温度与钢水碳元素含量为预测目标进行聚类,其中,以钢水温度为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水温度、主原料加入量四种属性划分的M1个炉次数据集合,以钢水碳元素含量为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水碳元素含量、主原料加入量四种属性划分的M2个炉次数据集合;
设某钢种全部样本的集合为N,所含样本的数目为n,聚类的具体步骤如下:
Step1.初始化原始数据矩阵
xim表示第i个样本第m个属性量的观测值;
Step2.对数据进行归一化处理,x′ik为x′ik归一化后的结果;
得到归一化后数据矩阵
Step3.计算数据点X′i.的密度函数Di,
其中,Xi.为矩阵第i行数据;
令Dc1=max{D1,D2,…,Dn},Xc1作为第1个聚类中心,ra为系数。
令Dck+1=max{D1,D2,…,Dn},Xck+1作为第k+1个聚类中心,η为阈值。
Step6.计算其余各点到Dck(k=1,2,…,m)的距离,根据距离大小判别所属类。
在以下过程中记预测钢水温度的集合中第i个分类为Q1i,预测钢水碳元素含量的集合中第i个分类为Q2i,对每一分类,使用插值法将离散数据拟合成连续数据,将数据完善。
本实施方式中,在动态优化精英炉次数据集合中按照钢种不同进行区分,划分出14个大类,钢种分别是GL4G、J24F、JT5q等共14种;然后,在同一钢种的炉次信息中,分别针对钢水温度和碳元素含量,按照不同属性采用聚类方法将动态优化精英炉次数据集合划分为不同性能的分类。
如钢种为GL4G的全部样本的集合为N,所含样本的数目为n=3740,聚类算法中设定阈值η=1.5,针对钢水温度的分类和针对钢水碳元素含量的分类的属性数目均为4,取m=4;在本实施例中经过计算得到M1=16,M2=20;在以下过程中记预测钢水温度的集合中第i个分类为Q1i,预测钢水碳元素含量的集合中第i个分类为Q2i。对每一分类,使用插值法将离散数据拟合成连续数据,将数据完善。
步骤2:采用最小二乘支持向量机方法(Least Squares Support Vector Machine,LSSVM),建立钢水温度预测模型Y′T(t+1)=LSSVM1{T(t),Sc(t),HI(t),N(t),O(t),GF(t),CO(t),CO2(t),V1(t),Hin1,Hin2}和碳元素含量预测模型Y′C(t+1)=LSSVM2{T(t),Sc(t),HI(t),N(t),O(t),GF(t),CO(t),CO2(t),V1(t),Hin1,Hin2};并得到操作变量与转炉炼钢过程中钢水温度与钢水碳含量之间的对应关系;钢水温度预测模型的输入和碳元素含量预测模型的输入均为t时刻的钢水温度T(t)、钢水碳元素含量SC(t)、t时刻到t+1时刻的时间范围内加入的氧气总量N(t)、顶吹氧气供气模式V1(t)、底吹气体总量O(t)、各副原料加入的重量HI(t)、烟气量GF(t)、烟气中一氧化碳量CO(t)、烟气中二氧化碳量CO2(t)、入炉铁水重量Hin1和入炉废钢重量Hin2,钢水温度预测模型的输出为下一时刻t+1时刻预测的钢水温度,碳元素含量预测模型的输出为下一时刻t+1时刻预测的钢水碳元素含量;
最小二乘支持向量机方法的步骤如下:
Step1.对于给定的训练数据集,利用一个非线性映射将样本从原空间映射到特征空间,在高维特征空间中构造最优决策函数,将非线性估计函数转化为高维特征空间中线性估计函数;
Step2.利用结构风险最小化原理,将回归问题表示成约束优化问题,用拉格朗日法求解这个优化问题;
Step3.定义核函数,将优化问题转化为求解线性方程组;
Step4.用最小二乘法求出回归系数和偏差得到非线性预测模型.
步骤3:在吹炼后期下副枪取样得到当前时刻钢水温度与碳元素含量信息,作为当前时刻的钢水温度和钢水碳元素含量。
炉次号=1368;钢种=GL4G;当前时刻的钢水温度T(t)=1608℃;当前时刻的钢水碳元素含量S(t)=0.05412%;
步骤4:根据钢水温度预测模型和碳元素含量预测模型确定转炉炼钢过程中的钢水温度和钢水碳元素含量;
步骤4.1:采集当前生产过程中的钢水温度与钢水碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量信息;
步骤4.2:选择与步骤4.1采集的当前生产过程中的信息匹配的步骤1.2中划分的炉次数据集合,并利用步骤2建立的钢水温度预测模型和碳元素含量预测模型对钢水温度与钢水碳元素含量进行预测;
计算步骤4.1采集的信息与步骤1.2中划分的各炉次数据集合之间的距离,该距离定义为当前生产状态与分类Q1i和Q2i中各样本之间差异的最小值。即针对钢水温度预报,考虑当前钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度四种属性,其中表示当前生产状态属性向量,表示分类Q1i中的一个样本的生产状态属性向量,为状态差异,其中:P为属性的总数量,P=4,x1iτ表示当前生产状态第τ属性向量值,y1iτ表示分类Q1i中的一个样本的生产状态第τ属性向量值。再次,从中选取最小的距离所对应的炉次数据集合作为当前信息状态匹配的集合。钢水碳元素含量预报与钢水温度预报相同,考虑当前钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水碳元素含量四种属性,得到最小的距离所对应的炉次数据集合m2。在计算时,如果当前预测以为钢水温度为主体,则使用分类Q1i,如果当前预测以钢水中碳元素含量为主体,则使用分类Q2i。最后,利用步骤2建立的钢水温度预测模型和碳元素含量预测模型对钢水温度与钢水碳元素含量进行预测;
当前的钢水温度与碳元素含量取步骤3中取得的数据。入炉主原料加入量=175897吨;目标出钢钢水温度=1674℃;目标出钢钢水碳元素含量=0.025782%;计算当前生产情况与步骤1.2中所获得的各分类之间的距离,该距离定义为当前生产状态与分类Q1i和Q2i中各样本之间差异的最小值。分别得到分类m1建立钢水温度的实时预测模型和分类m2建立钢水碳元素含量的实时预测模型。图3(a)表示温度的预测值与实际值比较的结果,图3(b)表示钢水碳元素含量预测与实际比较误差值,经多个炉次的比较结果,最大误差为0.49%,平均误差为0.26%。
步骤5:确定生产过程中各时刻钢水温度与碳元素含量,作为参考曲线;
步骤5.1:对动态优化精英炉次数据集合中同一钢种的炉次,以钢水温度和钢水碳元素含量为预测目标,分别按照当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量5个属性进行分类;
划分类别的方法采用步骤1.2中的方法,首先,在步骤1.1建立的动态优化精英炉次数据集合中按照钢种不同进行预化分;然后,在同一钢种的炉次信息中,以当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量5个属性在动态优化精英炉次数据集合中进行分类,在以下过程中记划分的集合中第i个分类为Ri。
本实施方式中,在动态优化精英炉次数据集合中按照钢种不同化分为14类;然后,在同一钢种的炉次信息中,以当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量5个属性在动态优化精英炉次数据集合中进行类别划分,共划分为35个类别。
步骤5.2:选择与生产过程中当前时刻的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量相匹配的步骤5.1中划分的炉次数据集合作为参考炉次集合,并从中选取最相似的h炉作为参考炉次;
按照当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度和目标出钢钢水碳元素含量这5个属性,依据步骤4的方法选择步骤5.1中对应的分类作为参考炉次集合,本实施方式选择分类R28作为参考炉次集。然后从参考炉次集合中选取h=5炉最相近的炉次,根据所选炉次各时刻钢水温度与碳元素含量的实际值。
例如最相近参考炉次
炉次号=1271;钢种=GL4G;当前时刻的钢水温度T(t)=1610℃;当前时刻的钢水碳元素含量S(t)=0.05116%;入炉主原料加入量=175488吨;目标出钢钢水温度=1673℃;目标出钢钢水碳元素含量=0.025782%;
步骤5.3:计算各参考炉次与当前炉次的钢水温度、碳元素含量和入炉主原料加入量的不同产生的偏差数;
其中:λ1,μ1,ν1分别为钢水温度偏差系数、碳元素含量偏差系数和主原料加入量偏差系数,
λ1=0.022,μ1=0.081,ν1=0.016;k为参考炉次。
步骤5.4:确定各参考炉次的权重系数;
计算参考炉次的权重系数
其中,k′为指定参考炉次,i=1,2,...,k′,...,h;(Ψ1(1),Ψ1(2),Ψ1(3),Ψ1(4),Ψ1(5))为(0.58,0.17,0.12,0.08,0.05)。
步骤5.5:确定当前炉次的未来各时刻最佳钢水温度与碳元素含量,作为参考曲线;
根据参考炉次计算值和权系数,用加权平均法计算出当前炉次的未来各时刻钢水温度T(t+1)与碳元素含量Sc(t+1)。其中Tj(t+1),Scj(t+1)是选择的参考炉次的钢水温度与碳元素含量。
动态控制阶段钢水温度和碳元素含量的计算结果如图6所示;表1为计算得到的动态控制阶段各副原料加入量及加入时间。
表1计算得到的动态控制阶段各副原料加入量及加入时间
时间 | 副原料1 | 副原料2 | 副原料3 | 副原料4 | 副原料5 | 副原料6 | 副原料7 |
0.01 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
0.03 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
… | … | … | … | … | … | … | … |
0.47 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
0.5 | 303 | 0 | 0 | 0 | 0 | 0 | 0 |
0.53 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
… | … | … | … | … | … | … | … |
1.53 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
各副原料参考值(单位:千克)
步骤6:根据烟气在线分析仪采集的转炉炼钢过程各个时刻的烟气流量、CO含量及CO2含量值与对应时刻的参考曲线中的烟气流量、CO含量及CO2含量值进行比较得到偏差值,实时校正钢水温度预测模型和碳元素含量预测模型;
首先,采集当前转炉炼钢过程的反馈信息,将由烟气在线分析仪得到的当前t时刻的烟气流量、CO含量及CO2含量值与对应t时刻的参考曲线中的烟气流量、CO含量及CO2含量值进行比较,如果相对误差e(t)均小于一个给定的阈值则不进行误差校正,Y″T(t+1)=Y′T(t+1),Y″C(t+1)=Y′C(t+1);否则,对下一时刻t+1的钢水温度与碳元素含量的预测值按照校正值进行校正,即钢水温度的校正值Y″T(t+1)=Y′T(t+1)×[1+r11×e(t)],钢水碳元素含量的校正值Y″C(t+1)=Y′C(t+1)×[1+0.5×r12×e(t)],式中Y′T(t+1)和Y′C(t+1)是根据预测模型得到的钢水温度和碳元素含量的预测值,r11和r12是在[-1.0,1.0]之间服从均匀分布的随机数,该随机数产生后,在下面的步骤7中将保持不变。比较结果如图7所示,当前实际烟气流量检测值与参考值的比较结果如图7(a)所示;当前实际CO、CO2在烟气中的含量检测值与参考值的比较结果如图7(b)所示。
步骤7:建立转炉炼钢过程实时动态优化模型,该模型以最小化当前预测的钢水温度和碳元素含量值与参考曲线的偏差为目标,该模型操作变量是当前时刻到下一时刻的时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量和各副原料加入的重量;
以最小化当前预测的钢水温度和碳元素含量值与参考曲线的偏差为目标函数,以各操作变量的可控制范围为约束,建立转炉炼钢过程实时动态优化模型,该模型操作变量是当前时刻到下一时刻的时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量和各副原料加入的重量;
优化目标为:
min(YT″(t+1)-YT(t+1))2 (11)
minξ(YC″(t+1)-YC(t+1))2 (12)
其中,ξ当YC″(t+1)≤YC(t+1)时取0,否则取eYC″(t+1)-YC(t+1),还表示了钢水碳元素含量尽量不超过规定要求值。
约束条件为:
HI(t)min≤HI(t)≤HI(t)max (13)
N(t)min≤N(t)≤N(t)max (14)
O(t)min≤O(t)≤O(t)max (15)
V1(t)min≤V1(t)≤V1(t)max (16)
Y′T(t+1)=LSSVM1{T(t),SC(t),HI(t),N(t),O(t),GF(t),CO(t),CO2(t),V1(t),Hin1,Hin2} (17)
Y′C(t+1)=LSSVM2{T(t),SC(t),HI(t),N(t),O(t),GF(t),CO(t),CO2(t),V1(t),Hin1,Hin2} (18)
YT(t+1)=T(t+1) (19)
YC(t+1)=SC(t+1) (20)
Y″T(t+1)=Y′T(t+1)×[1+r11×e(t)] (21)
Y″C(t+1)=Y′C(t+1)×[1+0.5×r12×e(t)] (22)
其中,HI(t)min和HI(t)max分别为各副原料进料流量的最小允许值和各副原料进料流量的最大允许值;N(t)min和N(t)max分别为底吹气体流量的最小允许值和最大允许值,O(t)min和O(t)max分别为氧枪供氧的最小允许流量值和最大允许流量值;V1(t)min和V1(t)max分别为氧***式中允许的最小高度和最大高度值;YT(t+1)与YC(t+1)为步骤4.5中得到的各时刻钢水温度与碳元素含量。
步骤8:确定转炉炼钢过程实时动态优化模型的各操作变量的优化设定值集合;
使用改进的粒子群算法(Particle Swarm Optimization,PSO)求解步骤7中所建立的转炉炼钢过程实时动态优化模型,该模型的解为各控制变量组成的向量X=[HI(t),N(t),O(t),V1(t)],获得各控制变量最俸取值的集合。如果一个解Xi的钢水温度和碳元素含量与参考值的差值均小于等于另一个解Xj的差值,并且至少有一个差值要取小于号,则称解Xi支配解Xj。解集合的任意两个不同的解Xi和Xj都是非劣解,即Xi不能支配解Xj,同时Xj也不能支配解Xi,同时,在算法所搜索到的可行域内不存在其它的解能够支配该集合中的任意一个解。经计算得到最好的20个非劣解集合组成精英解集G,即操作变量的优化设定值集合,全局最优解Xgj=(HI(t)gj,N(t)gj,O(t)gj,V1(t)gj)。
具体步骤如下:
Step1.产生由n个解组成的初始种群
在初始化范围内,对种群中每个粒子的位置xij=(HI(t)ij,N(t)ij,O(t)ij,V1(t)ij)和速度进行随机初始化,即解中的控制变量的值在约束(3)-(6)所要求的范围内容随机产生;i=1,2,...,n;第i个解在j维空间的值记为Xij=(HI(t)ij,N(t)ij,O(t)ij,V1(t)ij),其中HI(t)ij,N(t)ij,O(t)ij,V1(t)ij分别表示第i个解的副原料I加入量、底吹气体吹入量、氧气吹入量、氧枪供气模式,j是解的维数;并随机选取初始种群中的20个非劣解存入到精英解集合G中,精英解集合的作用是存储算法在迭代过程中所得到的非劣解;种群中的非劣解是指当前种群中不存在其它的解能够支配该解。
Step2.初始化相关参数
种群规模N=100,最大迭代次数m_nMaxIteration=100,粒子最大允许未改进代数m_nMaxNoImprove=10,单个粒子允许最大连续未改进代数m_nParticleMaxNoUpdate=5,每次迭代获得的精英解规模m_nElitePoolSize=20,惯性权重初始值ω=0.98,个体认知系数c1=2,社会认知系数c2=2,精英解学习系数c3=2,粒子更新速度上限vmax=5,下限vmin=2。
Step3.计算每个粒子的适应值。
以中心点钢水温度为优化目标,选择适值函数为中心点温度预测值与真实值之间的均方根误差最小。以中心点钢水碳元素含量为优化目标,选择适值函数为中心点碳元素含量预测值与真实值之间的均方根误差最小。
Step4.更新每个粒子的个体最优位置。
对于每个粒子,将其每次计算的适应值与自已所经历的历史最优位置的适应值相比较,若更好,则将其作为粒子的个体最好适应值,并用当前位置更新个体最优位置,pej=(HI(t)ej,N(t)ej,O(t)ej,V1(t)ej)。
Step5.更新所有粒子的全局的最优位置。
对所有粒子,将其适应值和群体内粒子经所历过的全局最优位置的适应值进行比较,若更好,则将其作为全局最好适应值,并用粒子的当前位置更新全局最优解粒子位置,pgj=(HI(t)gj,N(t)gj,O(t)gj,V1(t)gj)。
Step6.更新精英解集G
在PSO算法中引入的精英解集合G,用于存储PSO算法所得到的前20个质量最好的解,plj=(HI(t)lj,N(t)lj,O(t)lj,V1(t)lj)。用粒子的个体最俸解更新精英池中最差的精英解,其更新过程按照以下情况进行:(1)如果Xij的温度和碳元素含量与参考值的差值均大于或者等于精英解集G中的某一个解,则丢弃Xij;(2)如果Xij的温度和碳元素含量与参考值的差值均小于或者等于精英解集G中的某一个或者多个解,即Xij支配这些解,则将这些解从G中删除,并将解Xij***到G中;(3)如果Xij与G中的解互为非劣解,则将Xi***到G中。
Step7.更新粒子的位置和速度
根据式(1)和(2)速度更新公式和位置更新公式,对粒子进行更新。
速度更新公式:
vij(t+1)=ω(t+1)vij(t)+c1r1(t)(pej(t)-xij(t))+c2r2(t)(pgj(t)-xij(t))+c3r3(t+1)(plj(t)-xij(t)) (35)位置更新公式:
xij(t+1)=xij(t)+vij(t+1) (36)其中,i=1,2,...,n,j=1,2,...,w,i表示粒子编号,j表示单个粒子中维数的编号,惯性权重初始值ω=0.98,个体认知系数c1=2,社会认知系数c2=2,精英解学习系数c3=2,t表示迭代次数,r1,r2和r3是两个相对独立的[0,1]之间均匀分布的随机数。
Step8.终止准则
判断终止条件,当迭代次数达到最大迭代次数m_nMaxIteration或粒子最大允许未改进代数m_nMaxNoImprove时,满足算法的终止条件,停止,输出精英解集G,及全局最优解,Xgj=(HI(t)gj,N(t)gj,O(t)gj,V1(t)gj)。否则,转Step2.,继续进行迭代。
步骤9:根据转炉炼钢生产实际,从操作变量的优化设定值集合中选择操作变量设定值,根据操作变量设定值对转炉生产进行操作;
将步骤8所得到的精英解集G和全局最优解及与之对应的钢水温度和碳元素含量与参考值的差值,提供给现场操作人员,现场操作人员可以根据现场实际工况选择一个解。在本实施方式中,现场操作人员可提前对这些选项进行设定,可以直接选择最俸解,也可也根据实际工况变更选项,由图1中的动态优化实时控制器自动选择一个合适的解,并将当前时刻t对应的各操作变量的设定值下发到控制器,对副原料的加入、氧气流量、底吹气体流量、顶吹氧气供气模式等控制变量执行相应的调整操作。
计算得到的当前炉次动态控制阶段的最优控制参数设定值和控制器实际输出值如图5所示,计算得到的当前实际顶吹氧气流量的最优设定值与控制器输出值的比较结果如图5(a)所示;计算得到的际前底吹气体流量的最优设定值与控制器输出值的比较结果如图5(b)所示;计算得到的氧***式的最优设定值与控制器输出值的比较结果如图5(c)所示。
计算得到的实际碳元素含量和实际钢水温度分别如图6(a)和(b)所示。表2为当前炉次的动态控制阶段各副原料实际加入量及加入时间,时间单位为分钟。
表2当前炉次的动态控制阶段各副原料实际加入量及加入时间
时间 | 副原料1 | 副原料2 | 副原料3 | 副原料4 | 副原料5 | 副原料6 | 副原料7 |
0.01 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
0.03 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
… | … | … | … | … | … | … | … |
0.47 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
0.5 | 300 | 0 | 0 | 0 | 0 | 0 | 0 |
0.53 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
… | … | … | … | … | … | … | … |
1.53 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
各副原料实际加入值(单位:千克)
步骤10:若预测得到的钢水温度值与参考曲线中的参考值误差大于1.8%,或者钢水碳元素含量值与参考曲线中的参考值误差大于2.3%,则返回步骤5.1;否则,返回步骤5.2,进行转炉炼钢过程的实时动态优化与控制,直到出钢结束;
步骤11:定期对离线预测模型数据库进行更新。
定期对动态优化精英炉次数据集合进行更新,可使集合中炉次信息接近当前实际炉况,从而增强模型的自适应能力。
基于某炼钢厂的实际生产数据,从本实施例所得到的操作优化模型的解指导实际生产得到本炉次实际出钢钢水温度(图6)为1673℃,命中误差小于0.6%,实际出钢钢水碳元素含量为0.026%,满足工艺要求;多次重复本实施例的过程,所得到的多个炉次实际出钢钢水温度和碳元素含量值与目标值的比较结果图4(a)所示,误差如图4(b)所示,钢水碳元素含量的命中精度最小为0.9527%,平均命中精度为0.97%,温度最大偏差为7.52度,平均偏差为3.082度。同时该方法能够提供给现场操作人员多种操作设定与控制方案,方便操作人员根据实际工况需求灵活设定选择条件,进而帮助提高炼钢厂的整体生产效率,实验结果表明生产效率能够提高0.3-0.5%。
Claims (2)
1.一种基于数据驱动的转炉炼钢过程在线动态最优控制方法,其特征在于:包括以下步骤:
步骤1:建立基于数据驱动的离线预测模型数据库;
步骤1.1:建立动态优化精英炉次集合;
步骤1.1.1:在炼钢过程的吹炼后期至出钢结束,每间隔一个时间周期投掷一次碳温氧探头,测量炼钢转炉中的钢水温度和碳元素含量信息;
步骤1.1.2:将碳温氧探头测量的钢水温度与光谱分析仪测得的转炉炉口火焰温度进行比对,将碳温氧探头测量的钢水中碳元素含量与烟气分析仪得到的碳元素含量进行比对,得到炼钢转炉内的钢水实际温度与实际碳元素含量;
步骤1.1.3:选出钢水实际温度、碳元素含量与目标值相符的炉况稳定的炉次建立动态优化精英炉次数据集合,该集合中每一炉次的信息包括:测量范围内的各个时刻的钢水温度、钢水碳元素含量、顶吹氧气模式、不同时间范围内加入的氧气总量、底吹气体总量、各副原料加入的重量、烟气量、烟气中一氧化碳量、烟气中二氧化碳量、入炉铁水重量和入炉废钢重量;
步骤1.2:对动态优化精英炉次数据集合分别以钢水温度和钢水碳元素含量为预测目标按照各自属性进行分类;
步骤1.2.1:将动态优化精英炉次数据集合中的炉次按照钢种进行预划分;
步骤1.2.2:将动态优化精英炉次数据集合中同一钢种的炉次,分别以钢水温度与钢水碳元素含量为预测目标进行聚类,其中,以钢水温度为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水温度、主原料加入量四种属性划分的多个炉次数据集合,以钢水碳元素含量为预测目标进行聚类得到以炼钢转炉内的钢水实际温度、实际碳元素含量、实际出钢钢水碳元素含量、主原料加入量四种属性划分的多个炉次数据集合;
步骤2:采用最小二乘支持向量机方法,建立钢水温度预测模型和碳元素含量预测模型,钢水温度预测模型的输入和碳元素含量预测模型的输入均为当前时刻钢水温度、当前时刻钢水碳元素含量、某一时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量、各副原料加入的重量、烟气量、烟气中一氧化碳量、烟气中二氧化碳量、入炉铁水重量和入炉废钢重量,钢水温度预测模型的输出为下一时刻预测的钢水温度,碳元素含量预测模型的输出为下一时刻预测的钢水碳元素含量;
步骤3:在吹炼后期下副枪取样得到当前时刻钢水温度与钢水中各元素成分信息,得到当前时刻的钢水温度和钢水碳元素含量;
步骤4:根据钢水温度预测模型和碳元素含量预测模型确定转炉炼钢过程中的钢水温度和钢水碳元素含量;
步骤4.1:采集当前生产过程中的钢水温度与钢水碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量信息;
步骤4.2:选择与步骤4.1采集的当前生产过程中的信息匹配的步骤1.2中划分的炉次数据集合,并利用步骤2建立的钢水温度预测模型和碳元素含量预测模型对钢水温度与钢水碳元素含量进行预测;
步骤5:确定生产过程中各时刻钢水温度与碳元素含量,作为参考曲线;
步骤5.1:对动态优化精英炉次数据集合中同一钢种的炉次,以钢水温度和钢水碳元素含量为预测目标,分别按照当前的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量5个属性进行分类;
步骤5.2:选择与生产过程中当前时刻的钢水温度与碳元素含量、入炉主原料加入量、目标出钢钢水温度、目标出钢钢水碳元素含量相匹配的步骤5.1中划分的炉次数据集合作为参考炉次集合,并从中选取最相似的 炉作为参考炉次;
步骤5.3:计算各参考炉次与当前炉次的钢水温度、碳元素含量和入炉主原料加入量的不同产生的偏差数;
步骤5.4:确定各参考炉次的权重系数;
步骤5.5:确定当前炉次的未来各时刻的钢水温度与碳元素含量,作为参考曲线;
步骤6:根据烟气在线分析仪采集的转炉炼钢过程各个时刻的烟气流量、CO含量及CO 2 含量值与对应时刻的参考曲线中的烟气流量、CO含量及CO 2 含量值进行比较得到偏差值,实时校正钢水温度预测模型和碳元素含量预测模型;
步骤7:建立转炉炼钢过程实时动态优化模型,该模型以最小化当前预测的钢水温度和碳元素含量值与参考曲线的偏差为目标,该模型操作变量是当前时刻到下一时刻的时间范围内加入的氧气总量、顶吹氧气供气模式、底吹气体总量和各副原料加入的重量;
步骤8:确定转炉炼钢过程实时动态优化模型的各操作变量的优化设定值集合;
步骤9:根据转炉炼钢生产实际,从操作变量的优化设定值集合中选择操作变量设定值,根据操作变量设定值对转炉生产进行操作;
步骤10:若预测得到的钢水温度值与参考曲线中的参考值误差大于1.8%,或者钢水碳元素含量值与参考曲线中的参考值误差大于2.3%,则返回步骤5.1;否则,返回步骤5.2,进行转炉炼钢过程的实时动态优化与控制,直到出钢结束;
步骤11:定期对离线预测模型数据库进行更新。
2.权利要求1所述的基于数据驱动的转炉炼钢过程在线动态最优控制方法采用的转炉炼钢过程在线动态控制***,包括炼钢转炉、氧枪及副枪,其特征在于:还包括光谱分析仪、烟气分析仪和投掷探头;光谱分析仪安装在炼钢转炉炉口,烟气分析仪安装在转炉烟道处,投掷探头采用碳温氧探头。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410114943.7A CN103882176B (zh) | 2014-03-25 | 2014-03-25 | 一种基于数据驱动的转炉炼钢过程在线动态控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410114943.7A CN103882176B (zh) | 2014-03-25 | 2014-03-25 | 一种基于数据驱动的转炉炼钢过程在线动态控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103882176A true CN103882176A (zh) | 2014-06-25 |
CN103882176B CN103882176B (zh) | 2015-09-30 |
Family
ID=50951293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410114943.7A Active CN103882176B (zh) | 2014-03-25 | 2014-03-25 | 一种基于数据驱动的转炉炼钢过程在线动态控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103882176B (zh) |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104714519A (zh) * | 2014-12-30 | 2015-06-17 | 东北大学 | 一种连续退火机组生产过程工艺参数设定与在线优化方法 |
CN105242545A (zh) * | 2015-10-31 | 2016-01-13 | 南京南瑞继保工程技术有限公司 | 一种基于工艺批次模型的转炉用氧流量的预测方法 |
CN105714014A (zh) * | 2016-03-30 | 2016-06-29 | 本钢板材股份有限公司 | 转炉氧枪/加料/温度制度综合简化模型***及操作方法 |
CN105974896A (zh) * | 2016-06-07 | 2016-09-28 | 东北大学 | 一种基于信息物理融合的转炉炼钢优化控制***及方法 |
CN106843171A (zh) * | 2016-12-28 | 2017-06-13 | 沈阳化工大学 | 一种基于数据驱动方式的运行优化控制方法 |
CN107356200A (zh) * | 2017-07-03 | 2017-11-17 | 浙江大学 | 基于渣块轨迹的煤粉锅炉炉内落渣测量方法和*** |
CN107937658A (zh) * | 2017-02-08 | 2018-04-20 | 辽宁鸿盛冶金科技有限公司 | 一种转炉冶炼自动智能控制*** |
CN108008702A (zh) * | 2016-11-01 | 2018-05-08 | 北京明诚技术开发有限公司 | 智能化自动炼钢监控方法及*** |
CN108359766A (zh) * | 2018-03-02 | 2018-08-03 | 北京科技大学 | 一种基于钢水连续测温的底吹氧气转炉冶炼过程控制方法 |
CN108364095A (zh) * | 2018-02-06 | 2018-08-03 | 北京科技大学 | 基于数据挖掘的炼钢生产过程钢水质量诊断方法 |
CN108760091A (zh) * | 2018-05-29 | 2018-11-06 | 中北大学 | 基于改进检定炉的热电偶传感器动态补偿***构建方法 |
CN108921232A (zh) * | 2018-07-31 | 2018-11-30 | 东北大学 | 一种热轧带钢冷却历史数据聚类及相似性度量方法 |
CN109064064A (zh) * | 2018-09-12 | 2018-12-21 | 中南大学 | 一种基于图像分析炉门流渣的实时预报炉内渣成分方法 |
TWI660047B (zh) * | 2017-06-30 | 2019-05-21 | 日商Jfe鋼鐵股份有限公司 | 轉爐之操作監視方法及轉爐之操作方法 |
CN109857067A (zh) * | 2018-12-23 | 2019-06-07 | 北京首钢自动化信息技术有限公司 | 一种大数据环境下的炼钢多工序温度协调控制***及方法 |
CN110358888A (zh) * | 2018-03-26 | 2019-10-22 | 上海梅山钢铁股份有限公司 | 基于参考炉次分组的转炉冶炼控制装置及控制方法 |
CN110373510A (zh) * | 2019-07-16 | 2019-10-25 | 中南大学 | 一种基于双模态切换的钢水质量窄窗口控制方法及*** |
JP2020015960A (ja) * | 2018-07-26 | 2020-01-30 | 日本製鉄株式会社 | 操業条件演算装置及び操業条件演算方法 |
CN111172346A (zh) * | 2020-01-07 | 2020-05-19 | 邯郸钢铁集团有限责任公司 | 基于co体积百分含量曲线动态判断转炉炉内状态的方法 |
CN111334637A (zh) * | 2020-03-25 | 2020-06-26 | 宝钢湛江钢铁有限公司 | 转炉氧枪高度和流量方式自识别方法及*** |
CN111334636A (zh) * | 2020-03-27 | 2020-06-26 | 宝钢湛江钢铁有限公司 | 基于副枪检测的转炉终点动态控制方法及*** |
CN111635976A (zh) * | 2020-06-30 | 2020-09-08 | 南京宝希智能技术有限公司 | 一种非接触式转炉钢水温度连续测量***及检测方法 |
CN111766179A (zh) * | 2020-07-08 | 2020-10-13 | 大唐环境产业集团股份有限公司 | 基于lssvm的石灰石浆液密度测量方法、***及设备 |
CN111881993A (zh) * | 2020-08-03 | 2020-11-03 | 长沙有色冶金设计研究院有限公司 | 面向铜锍吹炼过程的操作模式多层分级匹配优化方法 |
CN111893237A (zh) * | 2020-07-08 | 2020-11-06 | 北京科技大学 | 一种转炉炼钢的熔池碳含量及温度全程实时预测方法 |
CN113454413A (zh) * | 2019-02-19 | 2021-09-28 | 杰富意钢铁株式会社 | 操作结果预测方法、学习模型的学习方法、操作结果预测装置及学习模型的学习装置 |
CN114525375A (zh) * | 2022-03-25 | 2022-05-24 | 山东钢铁股份有限公司 | 一种转炉异常炉况控制方法及*** |
CN114854929A (zh) * | 2022-05-07 | 2022-08-05 | 北京科技大学 | 一种转炉喷吹co2钢液成分与温度实时动态预测的方法 |
CN115145211A (zh) * | 2022-09-06 | 2022-10-04 | 上海桁壹高科技有限公司 | 一种非开关量设备的控制方法、***、介质和电子设备 |
TWI781626B (zh) * | 2021-05-18 | 2022-10-21 | 中國鋼鐵股份有限公司 | 高爐鐵水溫度預測方法與系統 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101845531A (zh) * | 2010-05-07 | 2010-09-29 | 北京首钢自动化信息技术有限公司 | 一种转炉冶炼终点钢水碳、温度控制***及其方法 |
CN102690923A (zh) * | 2012-06-13 | 2012-09-26 | 鞍钢股份有限公司 | 一种转炉副枪过程碳含量预报方法 |
CN102965462A (zh) * | 2011-09-02 | 2013-03-13 | 鞍钢股份有限公司 | 一种转炉冶炼终点自动控制装置及控制方法 |
CN103103309A (zh) * | 2013-01-24 | 2013-05-15 | 北京科技大学 | 一种辅助预报转炉炼钢终点的方法 |
CN202989198U (zh) * | 2012-10-30 | 2013-06-12 | 湖南镭目科技有限公司 | 一种转炉内碳含量和温度检测*** |
CN103160640A (zh) * | 2013-02-26 | 2013-06-19 | 河北钢铁股份有限公司邯郸分公司 | 一种转炉炼钢过程中动态检测炉渣中锰、磷、硫含量的方法 |
CN103451353A (zh) * | 2013-08-30 | 2013-12-18 | 武汉钢铁(集团)公司 | 基于模拟副枪静态检测的启动转炉动态控制模型的方法 |
-
2014
- 2014-03-25 CN CN201410114943.7A patent/CN103882176B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101845531A (zh) * | 2010-05-07 | 2010-09-29 | 北京首钢自动化信息技术有限公司 | 一种转炉冶炼终点钢水碳、温度控制***及其方法 |
CN102965462A (zh) * | 2011-09-02 | 2013-03-13 | 鞍钢股份有限公司 | 一种转炉冶炼终点自动控制装置及控制方法 |
CN102690923A (zh) * | 2012-06-13 | 2012-09-26 | 鞍钢股份有限公司 | 一种转炉副枪过程碳含量预报方法 |
CN202989198U (zh) * | 2012-10-30 | 2013-06-12 | 湖南镭目科技有限公司 | 一种转炉内碳含量和温度检测*** |
CN103103309A (zh) * | 2013-01-24 | 2013-05-15 | 北京科技大学 | 一种辅助预报转炉炼钢终点的方法 |
CN103160640A (zh) * | 2013-02-26 | 2013-06-19 | 河北钢铁股份有限公司邯郸分公司 | 一种转炉炼钢过程中动态检测炉渣中锰、磷、硫含量的方法 |
CN103451353A (zh) * | 2013-08-30 | 2013-12-18 | 武汉钢铁(集团)公司 | 基于模拟副枪静态检测的启动转炉动态控制模型的方法 |
Cited By (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104714519B (zh) * | 2014-12-30 | 2017-07-28 | 东北大学 | 一种连续退火机组生产过程工艺参数设定与在线优化方法 |
CN104714519A (zh) * | 2014-12-30 | 2015-06-17 | 东北大学 | 一种连续退火机组生产过程工艺参数设定与在线优化方法 |
CN105242545A (zh) * | 2015-10-31 | 2016-01-13 | 南京南瑞继保工程技术有限公司 | 一种基于工艺批次模型的转炉用氧流量的预测方法 |
CN105242545B (zh) * | 2015-10-31 | 2017-12-05 | 南京南瑞继保工程技术有限公司 | 一种基于工艺批次模型的转炉用氧流量的预测方法 |
CN105714014A (zh) * | 2016-03-30 | 2016-06-29 | 本钢板材股份有限公司 | 转炉氧枪/加料/温度制度综合简化模型***及操作方法 |
CN105974896B (zh) * | 2016-06-07 | 2018-08-28 | 东北大学 | 一种基于信息物理融合的转炉炼钢优化控制***及方法 |
CN105974896A (zh) * | 2016-06-07 | 2016-09-28 | 东北大学 | 一种基于信息物理融合的转炉炼钢优化控制***及方法 |
CN108008702A (zh) * | 2016-11-01 | 2018-05-08 | 北京明诚技术开发有限公司 | 智能化自动炼钢监控方法及*** |
CN106843171A (zh) * | 2016-12-28 | 2017-06-13 | 沈阳化工大学 | 一种基于数据驱动方式的运行优化控制方法 |
CN106843171B (zh) * | 2016-12-28 | 2019-10-01 | 沈阳化工大学 | 一种基于数据驱动方式的运行优化控制方法 |
CN107937658A (zh) * | 2017-02-08 | 2018-04-20 | 辽宁鸿盛冶金科技有限公司 | 一种转炉冶炼自动智能控制*** |
TWI660047B (zh) * | 2017-06-30 | 2019-05-21 | 日商Jfe鋼鐵股份有限公司 | 轉爐之操作監視方法及轉爐之操作方法 |
CN107356200A (zh) * | 2017-07-03 | 2017-11-17 | 浙江大学 | 基于渣块轨迹的煤粉锅炉炉内落渣测量方法和*** |
CN107356200B (zh) * | 2017-07-03 | 2021-06-29 | 浙江大学 | 基于渣块轨迹的煤粉锅炉炉内落渣测量方法和*** |
CN108364095A (zh) * | 2018-02-06 | 2018-08-03 | 北京科技大学 | 基于数据挖掘的炼钢生产过程钢水质量诊断方法 |
CN108364095B (zh) * | 2018-02-06 | 2022-05-13 | 北京科技大学 | 基于数据挖掘的炼钢生产过程钢水质量诊断方法 |
CN108359766A (zh) * | 2018-03-02 | 2018-08-03 | 北京科技大学 | 一种基于钢水连续测温的底吹氧气转炉冶炼过程控制方法 |
CN110358888B (zh) * | 2018-03-26 | 2021-06-08 | 上海梅山钢铁股份有限公司 | 基于参考炉次分组的转炉冶炼控制装置及控制方法 |
CN110358888A (zh) * | 2018-03-26 | 2019-10-22 | 上海梅山钢铁股份有限公司 | 基于参考炉次分组的转炉冶炼控制装置及控制方法 |
CN108760091A (zh) * | 2018-05-29 | 2018-11-06 | 中北大学 | 基于改进检定炉的热电偶传感器动态补偿***构建方法 |
JP2020015960A (ja) * | 2018-07-26 | 2020-01-30 | 日本製鉄株式会社 | 操業条件演算装置及び操業条件演算方法 |
JP7159677B2 (ja) | 2018-07-26 | 2022-10-25 | 日本製鉄株式会社 | 操業条件演算装置及び操業条件演算方法 |
CN108921232B (zh) * | 2018-07-31 | 2021-05-04 | 东北大学 | 一种热轧带钢冷却历史数据聚类及相似性度量方法 |
CN108921232A (zh) * | 2018-07-31 | 2018-11-30 | 东北大学 | 一种热轧带钢冷却历史数据聚类及相似性度量方法 |
CN109064064A (zh) * | 2018-09-12 | 2018-12-21 | 中南大学 | 一种基于图像分析炉门流渣的实时预报炉内渣成分方法 |
CN109064064B (zh) * | 2018-09-12 | 2021-11-23 | 中南大学 | 一种基于图像分析炉门流渣的实时预报炉内渣成分方法 |
CN109857067A (zh) * | 2018-12-23 | 2019-06-07 | 北京首钢自动化信息技术有限公司 | 一种大数据环境下的炼钢多工序温度协调控制***及方法 |
CN109857067B (zh) * | 2018-12-23 | 2021-01-05 | 北京首钢自动化信息技术有限公司 | 一种大数据环境下的炼钢多工序温度协调控制***及方法 |
CN113454413B (zh) * | 2019-02-19 | 2023-06-27 | 杰富意钢铁株式会社 | 操作结果预测方法、学习模型的学习方法、操作结果预测装置及学习模型的学习装置 |
CN113454413A (zh) * | 2019-02-19 | 2021-09-28 | 杰富意钢铁株式会社 | 操作结果预测方法、学习模型的学习方法、操作结果预测装置及学习模型的学习装置 |
CN110373510A (zh) * | 2019-07-16 | 2019-10-25 | 中南大学 | 一种基于双模态切换的钢水质量窄窗口控制方法及*** |
CN110373510B (zh) * | 2019-07-16 | 2020-04-24 | 中南大学 | 一种基于双模态切换的钢水质量窄窗口控制方法及*** |
CN111172346A (zh) * | 2020-01-07 | 2020-05-19 | 邯郸钢铁集团有限责任公司 | 基于co体积百分含量曲线动态判断转炉炉内状态的方法 |
CN111334637A (zh) * | 2020-03-25 | 2020-06-26 | 宝钢湛江钢铁有限公司 | 转炉氧枪高度和流量方式自识别方法及*** |
CN111334637B (zh) * | 2020-03-25 | 2021-07-13 | 宝钢湛江钢铁有限公司 | 转炉氧枪高度和流量方式自识别方法及*** |
CN111334636B (zh) * | 2020-03-27 | 2022-04-08 | 宝钢湛江钢铁有限公司 | 基于副枪检测的转炉终点动态控制方法及*** |
CN111334636A (zh) * | 2020-03-27 | 2020-06-26 | 宝钢湛江钢铁有限公司 | 基于副枪检测的转炉终点动态控制方法及*** |
CN111635976A (zh) * | 2020-06-30 | 2020-09-08 | 南京宝希智能技术有限公司 | 一种非接触式转炉钢水温度连续测量***及检测方法 |
CN111893237A (zh) * | 2020-07-08 | 2020-11-06 | 北京科技大学 | 一种转炉炼钢的熔池碳含量及温度全程实时预测方法 |
CN111893237B (zh) * | 2020-07-08 | 2021-11-09 | 北京科技大学 | 一种转炉炼钢的熔池碳含量及温度全程实时预测方法 |
CN111766179A (zh) * | 2020-07-08 | 2020-10-13 | 大唐环境产业集团股份有限公司 | 基于lssvm的石灰石浆液密度测量方法、***及设备 |
CN111881993A (zh) * | 2020-08-03 | 2020-11-03 | 长沙有色冶金设计研究院有限公司 | 面向铜锍吹炼过程的操作模式多层分级匹配优化方法 |
CN111881993B (zh) * | 2020-08-03 | 2024-04-12 | 长沙有色冶金设计研究院有限公司 | 面向铜锍吹炼过程的操作模式多层分级匹配优化方法 |
TWI781626B (zh) * | 2021-05-18 | 2022-10-21 | 中國鋼鐵股份有限公司 | 高爐鐵水溫度預測方法與系統 |
CN114525375A (zh) * | 2022-03-25 | 2022-05-24 | 山东钢铁股份有限公司 | 一种转炉异常炉况控制方法及*** |
CN114854929A (zh) * | 2022-05-07 | 2022-08-05 | 北京科技大学 | 一种转炉喷吹co2钢液成分与温度实时动态预测的方法 |
CN115145211A (zh) * | 2022-09-06 | 2022-10-04 | 上海桁壹高科技有限公司 | 一种非开关量设备的控制方法、***、介质和电子设备 |
CN115145211B (zh) * | 2022-09-06 | 2022-12-09 | 上海桁壹高科技有限公司 | 一种非开关量设备的控制方法、***、介质和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103882176B (zh) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103882176B (zh) | 一种基于数据驱动的转炉炼钢过程在线动态控制方法 | |
CN104630410B (zh) | 一种基于数据解析的转炉炼钢质量实时动态预测方法 | |
CN109447346B (zh) | 基于灰色预测与神经网络组合模型转炉氧耗量预测方法 | |
CN113239482B (zh) | 一种转炉后吹碳含量动态预测方法及装置 | |
CN108469180B (zh) | 基于大数据和机器学习的烧结终点预测***的建立方法 | |
CN106096788B (zh) | 基于pso_elm神经网络的转炉炼钢工艺成本控制方法及*** | |
CN108676955A (zh) | 一种转炉炼钢终点碳含量和温度控制方法 | |
CN110066895A (zh) | 一种基于Stacking的高炉铁水质量区间预测方法 | |
CN109359723A (zh) | 基于改进正则化极限学习机的转炉终点锰含量预测方法 | |
CN104894335A (zh) | 一种aod炉喷溅预报分析仪信息融合方法 | |
CN107368125B (zh) | 一种基于cbr与rbr并行混合推理的高炉炉温控制***及方法 | |
An et al. | A multi-time-scale fusion prediction model for the gas utilization rate in a blast furnace | |
CN101881981A (zh) | 一种rh钢水温度、成分闭环控制*** | |
CA2894813C (en) | Method and device for predicting, controlling and/or regulating steelworks processes | |
CN106096724B (zh) | 基于elm神经网络的转炉炼钢工艺成本控制方法及*** | |
CN106054836A (zh) | 基于grnn的转炉炼钢工艺成本控制方法及*** | |
Li et al. | Burden surface decision using MODE with TOPSIS in blast furnace ironmkaing | |
CN113761787A (zh) | 基于深度迁移网络的高炉铁水硅含量在线预测方法及*** | |
CN110427715A (zh) | 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 | |
CN108004368A (zh) | 智能化自动炼钢方法及装置 | |
Cardoso et al. | Prediction of silicon content in the hot metal using Bayesian networks and probabilistic reasoning | |
CN102373310A (zh) | 一种指导转炉补吹过程操作的方法 | |
CN102968644A (zh) | 一种氩氧精炼铁合金冶炼终点预测方法 | |
Wang et al. | Predicting the endpoint phosphorus content of molten steel in BOF by two-stage hybrid method | |
CN101947645A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |