基于机理和数据驱动的转炉输入—输出混合模型 | [PDF全文] |
钢铁行业是我国的国民经济支柱产业,同时也是能源消耗大户,能源消耗量占全国总能耗的15 %[1].钢铁行业节能一直是研究的重点,从20世纪80年代的单体性技术节能到目前以“能源转换功能”为主要指导时期的系统节能[2]. 2010年殷瑞钰院士指出应注重与铁素流相关的碳素能量流的输入/输出特点和能量流网络的构建[3].目前对于钢铁企业中能量流的研究较多,蔡九菊教授[4]针对物质和能量流动过程,构建了冶金工序的物质流图和能量流图;苍大强[5]对冶金工业节能新技术进行了深入分析;赵业清[6]利用混合petri网建模方法,构建了钢铁企业能量流网络系统混合Petri网模型;王懿等[7]对钢铁企业能量流网络的特点进行详细分析并指出未来发展趋势;郑忠等[8]对钢铁企业信息系统进行分析,探讨了基于信息流的物质流和能量流的协同优化机制;税烺等[9]利用火用分析和传统热平衡分析对冶金生产的余能回收进行研究;孙彦广等[10]借助能量流网络仿真对钢铁工业中的能源介质进行优化调配;张利娜等[11]对转炉石灰石替代石灰造渣进行了研究;常立忠等[12]借助BP神经网络建立了转炉静态模型,对转炉终点C含量进行预报研究;严良涛等[13]以基于遗传算法的核偏最小二乘回归方法,建立转炉终点C含量预测模型;Xu L F等[14]采用支持向量机建立了转炉终点的预测模型;其他学者也分别采用不同的方法对转炉的终点控制进行了研究[15-18].随着节能研究的深入,钢铁企业节能空间进一步缩减,因此需要对能量流网络进行精细化控制研究以达到优化节能的目的,其中冶金工序的输入-输出模型作为能量流网络精细化控制的底层支持.虽然目前对于能量流网络的研究较多,但是缺乏对冶金工序输入-输出模型的研究,因此文中以转炉工序为例,以机理和数据驱动相结合的方法构建转炉输入-输出模型,为能量流网络的精细化控制提供研究基础.
1 模型构建转炉在钢铁流程中作为一个中转装置,用于承载来自上一工序的能量流,并在该装置中产生一个输入-输出行为,在承载上一工序的能量流时,同时接受外加物质流.按照输入物质/能量流,外加物质/能流量、输出物质/能量流的流程模式对转炉工序进行解析,转炉工序的输入为:
1)下游工序返回流:废钢
2)上游工序顺流:来自KR预处理工序的铁水
3)外加消耗流:辅料(造渣剂、冷却剂)、合金、能源介质(氧气)
输出为:
1)顺流下游工序:钢水
2)返流上游工序:无
3)回收处理流:炉渣、炉气
4)直接排放流:烟尘
在工序解析的基础上,将已知量视为模型的输入,未知量作为模型的输出,则转炉工序的输入-输出模型如图 1所示,由于转炉工序的合金料是在出钢过程中加入,于是该模型不考虑合金项.
1.1 模型思路
该模型将机理方法与数据方法相结合,对模型输出项进行计算.先采用机理方法对模型的输出项进行机理计算,再采用以实际生产数据训练好的神经网络对出钢温度进行预测,以预测温度和机理计算时采用的目标温度进行对比,若差值的绝对值≤10 ℃,则满足要求,输出各计算项,若差值的绝对值大于10 ℃,则不满足要求,返回机理模型对热损失参数进行修正,直至差值达到要求,热损失参数的修正按照式(1)进行.
$ {\omega _{i + 1}} = {\omega _i} + \frac{{c{m_{{\rm{steel - LD}}}}\Delta T}}{{{Q_{{\rm{heat, input}}}}}} $ | (1) |
式(1)中:ωi+1为修正后的热损失参数;ωi为修正前的热损失参数;c为钢水比热容,kJ/(kg·℃);msteel为钢水量,kg;△T为数据模型的预测温度与目标钢水温度之差,△T=T1-T0,℃;Qheat,input为机理模型计算的热收入,kJ.
模型的具体流程如图 2所示.
1.2 机理模型控制方程
转炉冶炼造渣剂主要有石灰、轻烧白云石、氧化球团等.石灰是转炉冶炼的主要造渣料,它的主要成分为CaO,具有很强的脱磷和脱硫能力,石灰的加入量主要依据铁水中[Si]含量和终渣碱度R来确定:
$ {m_{{\rm{lime - LD}}}} = \frac{{R{\rm{ \times }}\sum {{\rm{Si}}{{\rm{O}}_2}} - \sum {{\rm{CaO}}} }}{{{\omega _{{\rm{CaO, lime}}}} - R \times {\omega _{{\rm{Si}}{{\rm{O}}_{\rm{2}}}{\rm{, lime}}}}}} $ | (2) |
式(2)中:mlime-LD为石灰加入量,kg;R为炉渣碱度;∑SiO2为除石灰外的其他物料带入的及氧化生成SiO2量总和,kg;∑CaO为除石灰外的其他物料带入的CaO量之和,kg;ωCaO,lime为石灰中CaO的成分含量,%;ωSiO2,lime为石灰中SiO2的成分含量,%.
轻烧白云石的主要成分为MgO和CaO,可以替代一部分石灰,加入轻烧白云石的主要目的是保证渣中有一定含量的MgO,减少前期酸性渣对转炉炉衬的侵蚀.轻烧白云石的加入量根据炉渣中所要求的MgO含量来确定,渣中MgO含量一般控制在7 %~10 %.
$ {m_{{\rm{dolomite - LD}}}} = \frac{{{m_{{\rm{slag - LD}}}} \times {\omega _{{\rm{MgO, slag}}}} - \sum {{\rm{MgO}}} }}{{{\omega _{{\rm{MgO, dolomite}}}}}} $ | (3) |
式(3)中:mdolomite-LD为轻烧白云石加入量,kg;mslag-LD为转炉炉渣量,kg;ωMgO,slag为炉渣MgO成分含量,%;∑MgO为除轻烧白云石外其他物料带入的MgO含量,kg;ωMgO,dolomite为轻烧白云石中MgO的成分含量,%.
氧化球团作为一种优良的冷却剂,不仅可以吸收富余热量,还能减少铁水的氧化,球团的加入量根据热平衡来计算:
$ {m_{{\rm{pllet - LD}}}} = \frac{{{Q_{{\rm{heat, input}}}} - \sum {{Q_{{\rm{heat, input}}}}} }}{{{q_{{\rm{pellt}}}}}} $ | (4) |
式(4)中:mpellet-LD为球团加入量,kg;∑Qheat,input为除球团吸热外的其他热支出之和,kJ;qpellt为单位质量球团吸热量,kJ/kg.
转炉氧气的消耗主要用于铁水中元素的氧化,元素氧化所需的氧大部分来源于氧枪吹入,小部分来源于球团供氧,氧气的吹氧量为:
$ {V_{{{\rm{O}}_{\rm{2}}}{\rm{-LD}}}} = \frac{{{V_{{{\rm{O}}_{\rm{2}}}{\rm{, 铁水中元素氧化}}}}\;\;\;\;\;\;-{V_{{{\rm{O}}_{\rm{2}}}{\rm{, pellet}}}}}}{{{\eta _{{O_2}}}}} $ | (5) |
式(5)中:VO2-LD为转炉吹氧量,m3;VO2,铁水中元素氧化为铁水中元素氧化耗氧量,m3;VO2,pellet为冷固球团供氧量,m3;ηO2为氧气利用系数.
转炉冶炼吹损主要有非Fe元素氧化损失、渣中铁损失、烟尘铁损失、喷溅铁损,假设球团中铁氧化物全被还原为Fe,则出钢量为:
$ \begin{array}{l} {m_{{\rm{steel - LD}}}} = {m_{{\rm{iron - LD}}}} + {m_{{\rm{scrap - LD}}}} - {m_{{\rm{Fe, slag}}}} - \sum {{m_{\Delta \left[i \right]}}} - \\ \;\;\;\;\;\;\;\;\;\;\;\;\;{m_{{\rm{Fe, dust}}}} + {m_{{\rm{Fe, pellet}}}} \end{array} $ | (6) |
式(6)中:miron-LD为入炉铁水量,t;mscrap-LD为入炉废钢量,t;mFe,slag为渣中全铁量,t;∑m△[i]为非Fe元素氧化量总和,t,其中i为[C]、[Si]、[Mn]、[P];mFe,dust为烟尘中全铁量,t;mFe,pellet为氧化球团中全铁量,t.
转炉炉渣主要来源于造渣剂和球团中非Fe成分固态氧化物、铁水中各元素氧化物,渣量的计算公式为:
$ {m_{{\rm{slag - LD}}}} = \frac{{{m_{{\rm{oxid, lime}}}} + {m_{{\rm{oxid, dolomite}}}} + {m_{{\rm{oxid, pellet}}}} + {m_{{\rm{oxid, iron}}}}}}{{1 - {\omega _{{\rm{FeO, slag}}}} - {\omega _{{\rm{F}}{{\rm{e}}_{\rm{2}}}{{\rm{O}}_{\rm{3}}}{\rm{, slag}}}}}} $ | (7) |
式(7)中:mslag-LD为炉渣量,kg;moxid,lime为石灰中固态氧化物总量,kg;moxid,dolomite为轻烧白云石中固态氧化物总量,kg;moxid,pellet为球团中非Fe固态氧化物总量,kg;moxid,iron为铁水中非Fe元素氧化生成物总量,kg;ωFeO,slag为渣中FeO成分含量,%;ωFe2O3,slag为渣中Fe2O3成分含量,%.
转炉炉气的主要成分为CO和CO2、N2,还有少许SO2、H2O和自由O2,其中CO和CO2来源于铁水中C的氧化,N2、自由O2来源于吹氧带入,SO2、H2O来源于铁水中S元素氧化和外加物料,冶炼一炉钢产生的炉气量为:
$ {V_{{\rm{gas}}}} = {V_{{\rm{CO}}}} + {V_{{\rm{C}}{{\rm{O}}_{\rm{2}}}}} + {V_{{{\rm{N}}_{\rm{2}}}}} + {V_{{\rm{S}}{{\rm{O}}_{\rm{2}}}}} + {V_{{{\rm{H}}_{\rm{2}}}{\rm{O}}}} + {V_{{{\rm{O}}_{\rm{2}}}{\rm{, free}}}} + {V_{{{\rm{O}}_{\rm{2}}}{\rm{, steel}}}} $ | (8) |
式(8)中:Vgas为炉气量,m3;VCO为C氧化生成的CO量,m3;VCO2为C氧化生成的CO2量,m3;VN2为氧气带入的N2量,m3;VSO2为铁水氧化生产和辅料中元素生成SO2量,m3;VH2O为辅料带入的水汽量,m3;VO2,free为自由氧量,m3;VO2,steel为钢中溶氧量,m3.
在传统的计算模型中,一般是先假定炉渣量,计算出轻烧白云石,然后计算出石灰,球团等,由于轻烧白云石能替代部分石灰,所以在计算中需要重复迭代计算.此模型计算对传统计算模型进行改进,不假定渣量,将上述式(2)、式(3)、式(4)转换为以辅料成分和已知参数为系数的方程组见式(9),进而转换为式(10)所示的矩阵形式进行求解,因为渣量、钢水量、炉气量贯穿在式(4)中,所以式(5)~式(8)在辅料计算后再单独计算.
$ \left\{ \begin{array}{l} {a_{11}}{m_{{\rm{lime - LD}}}} + {a_{12}}{m_{{\rm{dolomite - LD}}}} + {a_{13}}{m_{{\rm{pellet - LD}}}} = {b_1}\\ {a_{21}}{m_{{\rm{lime - LD}}}} + {a_{22}}{m_{{\rm{dolomite - LD}}}} + {a_{23}}{m_{{\rm{pellet - LD}}}} = {b_2}\\ {a_{31}}{m_{{\rm{lime - LD}}}} + {a_{32}}{m_{{\rm{dolomite - LD}}}} + {a_{33}}{m_{{\rm{pellet - LD}}}} = {b_3} \end{array} \right. $ | (9) |
$ \left[\begin{array}{l} \;\;{m_{{\rm{lime-LD}}}}\\ {m_{{\rm{dolomite-LD}}}}\\ \;\;{m_{{\rm{pellet-LD}}}} \end{array} \right] = {\left[\begin{array}{l} {a_{11}}\;{a_{12}}\;{a_{13}}\\ {a_{21}}\;{a_{22}}\;{a_{23}}\\ {a_{31}}\;{a_{32}}\;{a_{33}} \end{array} \right]^{ - 1}}\left[\begin{array}{l} {b_1}\\ {b_2}\\ {b_3} \end{array} \right] $ | (10) |
以S钢厂的4号转炉为研究对象,对该转炉一年的实际生产数据进行分析,为了避免不同钢种带来的差异影响,此模型仅以该厂生产最多的SPHC-W1钢种为例,采集到原始数据1 311组,删除数据缺失项后,再采用拉依达准则(准则)对数据进行处理,删除掉异常数据后剩余813组.在该转炉实际生产中,造渣剂的加入种类有石灰、轻烧白云石、白云石、氧化球团、含碳菱镁球、灰石渣6种,813组数据中,加入各种辅料的炉次占比情况分布如图 3所示.
从图(3)中可以看出,在统计的813炉数据中,加入石灰和氧化球炉次占比100 %,加入轻烧白云石的炉次占比92.4 %,加入含碳菱镁球、灰石渣、白云石的炉次占比均低于13 %,因此可知,生产过程中加入的物料主要有石灰、轻烧白云石、氧化球团3种,加入其它辅料的炉次视为物料水平不稳定,于是保留物料水平稳定的炉次699组,对该699炉分别计算它们的物料平衡和热平衡,计算时相关参数设定情况为:①C生成CO和CO2的比例为9:1;②炉气温度:1 520 ℃;③炉衬侵蚀量:0.3 %;④烟尘损失量1.6 %,其中FeO占比77 %,Fe2O3占比20 %;⑤氧气纯度:99.5 %;⑥炉气中自由氧含量:0.05 %;⑦根据炉渣化验情况,确定渣中FeO占比17 %,Fe2O3占比7 %;⑧对钢水终点O和C数据拟合,得到终点O(×106)与终点C(%)的关系:
$ \left[{\rm{O}} \right] = 52.743 \times {\left[{\rm{C}} \right]^{ - 0.832}} $ | (11) |
通过物料平衡和热平衡的计算,冶炼相关参数结果统计见表 1.
点击放大 |
其中:变异系数=标准差/平均值,用以表征某项目的波动情况;
热效率=(钢水物理热+废钢物理热+氧化球团物理热)/热收入.
从表 1中可以看出氧气利用率、渣中MgO和热效率3个参数比较稳定,模型中直接取其平均值作为模型中冶炼参数.碱度由于铁水中[Si]含量的波动而产生较大的波动,另外碱度还与铁水中的[P]含量和铁水温度有关,所以根据实际渣料加入量计算出的碱度与铁水中[Si]、[P]以及铁水温度T进行多元回归,得到式(12)所示关系,于是模型中碱度参数按照式(12)进行赋值.
$ R = 14.67 \times \left[{\rm{P}} \right] - 0.003 \times T - 9.67 \times \left[{{\rm{Si}}} \right] + 9.93 $ | (12) |
在采用神经网络对出钢温度进行预测时,首先需要确定输入参数,由机理分析可知,热收入主要有铁水物理热和元素氧化热,除去钢水物理热外,热支出有冷却剂(废钢、氧化球团)吸热、炉渣热、热损失,铁水物理热由铁水温度和铁水量决定,元素氧化热由铁水成分含量、目标钢水成分含量决定,炉渣热包含造渣剂的吸热量,热损失则与转炉设备有关,因此可知出钢温度与铁水量、铁水成分、目标钢水成分、冷却剂量、造渣剂加入量以及炉龄都有关系,为了减少输入变量个数,对铁水成分和目标钢水成分进行计算,以元素百分氧化量替代相应输入变量,最终确定以铁水量、C氧化量、Si氧化量、Mn氧化量、P氧化量、铁水温度、石灰加入量、轻烧白云石加入量、废钢、氧化球团加入量、吹氧量、炉龄作为输入变量,以650炉生产数据作为训练数据,49炉作为测试数据.利用MATLAB软件建立一个多输入单输出的3层BP神经网络,输入层节点数和输出层节点数分别为输入的变量个数和输出的变量个数,即输入层节点数为12,输出层节点数为1,以logsig函数和tansig函数作为传递函数,traingdx作为训练函数,通过试凑法确定较优隐含层节点数为9,测试结果如图 4所示.
测试结果精度汇总于表 2.
点击放大 |
用机理方法对上述49炉终点钢水温度进行计算,误差在±10 ℃范围内的炉次只有5炉,占比为10.20 %,对比表 2知,神经网络预测精度在±10 ℃误差范围内为75.51 %,明显优于机理计算,于是将训练好的神经网络模型保存,用于后续模型计算.
3 模型验证将相关参数取值代入模型中,并采用C#编程语言,将模型的输出项控制方程程序化,在计算出钢温度时,直接调用训练好的神经网络,以机理计算出的输出参数作为数据驱动模型的输入变量,对出钢进行预测计算,并将结果与目标温度进行对比,若误差未满足要求则返回调整热损失参数继续计算,直至误差在要求范围内.简易的程序界面如图 5所示.
为了验证文中模型的可行性,以数据驱动模型中预测精度在±10 ℃误差范围内的36炉数据进行验证计算.对该36炉数据分别采用单一的机理模型和本文的混合模型进行计算,结果见表 3和表 4.
点击放大 |
点击放大 |
生产数据中未采集单炉炉气量和炉渣量数据,因此无法对比;出钢量由于钢包容量限制的原因,钢水产量与出钢量并不一定相等,于是模型结果不做比较.由表 3可知,在相同误差范围内,该文混合模型的命中率相较于机理模型都有所提高,石灰加入量命中率提高11.1 %,轻烧白云石加入量命中率提高8.3 %,氧化球团加入量命中率提高8.3 %;转炉吹氧量主要用于铁水中碳氧化,与渣料加入量的多少关系不大,因此2种模型计算出的结果差距不大.
4 结论1)构建了基于机理和数据驱动的转炉输入-输出混合模型,以数据驱动模型对机理模型进行校正,提高了模型精度.
2)以C#语言将模型程序化,以简易界面可视化显示,并用该程序对36炉数据进行输入-输出计算,在±1 t误差范围内,石灰加入量命中率为86.1 %,轻烧白云石加入量命中率为88.8 %,球团加入量命中率为77.7 %,相较于纯机理模型分别提高了11.1 %、8.3 %、8.3 %;吹氧量命中率在±7 %误差范围内为91.7 %,比纯机理模型提高2.9 %.
[1] |
王俊岭, 张新社. 中国钢铁工业经济增长、能源消耗与碳排放脱钩分析[J].
河北经贸大学学报, 2017, 38(4): 77–82.
|
[2] |
殷瑞钰. 从开放系统、耗散结构到钢厂的能量流网络化集成[J].
中国冶金, 2010, 20(8): 1–14.
|
[3] |
殷瑞钰. 钢铁制造流程的能量流行为和能量流网络问题[J].
工程研究-跨学科视野中的工程, 2010, 2(1): 1–4.
|
[4] |
蔡九菊, 王建军, 陆钟武, 等. 钢铁企业物质流与能量流及其相互关系[J].
东北大学学报(自然科学版), 2006, 27(9): 979–982.
|
[5] |
苍大强. 国内外冶金工业源头节能减排的新方法、新技术[J].
有色金属科学与工程, 2015, 6(6): 1–6.
|
[6] |
赵业清. 基于混合Petri网的钢铁企业能量流网络模型研究[J].
冶金自动化, 2014, 38(5): 27–32.
|
[7] |
王懿, 芦永明. 钢铁企业能量流网络研究现状与发展趋势[J].
冶金能源, 2016, 35(3): 3–8.
|
[8] |
郑忠, 黄世鹏, 李曼琛, 等. 钢铁制造流程的物质流和能量流协同优化[J].
钢铁研究学报, 2016, 28(4): 1–7.
|
[9] |
税烺, 贺东风, 艾立翔, 等. 冶金生产余能回收的一种新的能量分析法[J].
有色金属科学与工程, 2012, 3(1): 43–48.
|
[10] |
孙彦广, 梁青艳, 李文兵, 等. 基于能量流网络仿真的钢铁工业多能源介质优化调配[J].
自动化学报, 2017, 43(6): 1065–1079.
|
[11] |
张利娜, 潘贻芳, 袁章福, 等. 转炉复吹与石灰石造渣行为控制技术的研究[J].
有色金属科学与工程, 2015, 6(3): 16–21.
|
[12] |
常立忠, 李正邦. 基于BP神经网络的转炉静态模型[J].
炼钢, 2006, 22(6): 41–44.
|
[13] |
严良涛, 李鸣, 杨大勇. 基于GA-KPLSR的转炉终点碳含量的预测研究[J].
控制工程, 2017, 24(5): 923–926.
|
[14] |
XU L F, LI W, ZHANG M, et al. A model of basic oxygen furnace (BOF) end-point prediction based on spectrum information of the furnace flame with support vector machine (SVM)[J].
Optik -International Journal for Light and Electron Optics, 2011, 122(7): 594–598. DOI: 10.1016/j.ijleo.2010.04.018. |
[15] |
胡燕. 基于粗糙集的转炉炼钢知识发现及终点控制模型研究[D]. 重庆: 重庆大学, 2013.
http://cdmd.cnki.com.cn/Article/CDMD-10611-1013044413.htm |
[16] |
王正, 孙彦辉, 赵志超, 等. 200 t转炉终点控制对脱磷的影响[C]//2015年全国特钢年会论文集. 北京: 冶金工业出版社, 2015: 44-52.
|
[17] |
ZHANG B, XUE Z L, LIU K, et al. Development and application of prediction model for end-point manganese content in converter based on data from sub-lance[J].
Advanced Materials Research, 2013, 683: 497–503. DOI: 10.4028/www.scientific.net/AMR.683. |
[18] |
DONG X L, DONG S. The converter steelmaking end point prediction model based on RBF neural network[J].
Applied Mechanics & Materials, 2014, 577: 98–101. |