CN104531924B - 一种高炉炉料分布实时预报系统及方法 - Google Patents
一种高炉炉料分布实时预报系统及方法 Download PDFInfo
- Publication number
- CN104531924B CN104531924B CN201410797862.1A CN201410797862A CN104531924B CN 104531924 B CN104531924 B CN 104531924B CN 201410797862 A CN201410797862 A CN 201410797862A CN 104531924 B CN104531924 B CN 104531924B
- Authority
- CN
- China
- Prior art keywords
- furnace
- chute
- charge level
- charge
- 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
Links
Landscapes
- Blast Furnaces (AREA)
Abstract
本发明提供一种高炉炉料分布实时预报系统及方法,包括采集数据器、数据处理器、初始料面设定器、同档位下新料面检测器、多档位下新料面检测器、料面下降模块、料面分布信息生成器和布料控制器;方法包括获取当前高炉生产过程工艺数据;处理高炉生产过程工艺数据中的炉料重量及化验数据;设定初始料面;构造相同布料档位下的新料面;形成多档位下新料面;构造下降料面;计算料层分布信息并绘制料层分布图像;根据高炉炉料分布实时预报的料层分布信息指导生产操作。本发明采用了机理模型与数据驱动相结合的方法,可以根据矿焦分布炉况的判断和布料模型的输出结果,对高炉矿焦分布进行合理调整,抑制炉况波动,降低炉况异常引起的休风率。
Description
技术领域
本发明属于高炉炉料分布数字模拟技术领域,特别涉及一种高炉炉料分布实时预报系统及方法。
背景技术
高炉冶炼是一个连续的生产过程,全过程是炉料自上而下,煤气自下而上的相互接触过程中完成的。炉料从炉顶装入炉内,从风口鼓入由热风炉加热到1000~1300℃的热风,炉料中焦炭在风口前与鼓风中的氧发生燃烧反应,产生高温还原气体,在炉内上升过程中加热缓慢下降的炉料,并还原铁矿石中的铁氧化物为金属铁。矿石升到一定温度后软化,熔融滴落,矿石中未被还原的物质形成熔渣,实现渣铁分离。已熔化的渣铁聚集在炉缸内,发生诸多反应,最后调整铁的成分和温度达到终点,定期从炉内排放炉渣和生铁。
高炉生产过程特点有:
1)过程的复杂性。高炉是目前为止最复杂的一种冶金反应器,发生在高炉内部,特别是高炉下部的各种物理化学现象,至今还没有得到完全和充分的认识。
2)检测信息不完全。具体表现在检测点少,且多局限于过程的边界(炉顶、渣铁、风口、炉身静压力等),高炉整个冶炼过程的炉内反应和变化无法直接观察,检测手段和方法受到一定的限制。
3)生产过程大滞后,被控过程的响应速度缓慢。
4)可控范围狭窄。这有两方面的含义,一是因为铁水质量必须满足炼钢等用户的严格要求,使得高炉的各种控制参数的可调范围相对比较小;二是炉况必须早调、小调,才不至于发生过大的波动,否则炉况将会急剧恶化导致失控。
高炉的这些特点,使得高炉的过程控制,特别是实时在线控制成为一项特别复杂的工作。虽然生产人员已经在自动化过程控制方面积累了相当多的经验,也取得了一定的成效,但是由于炼铁工艺的要求越来越精细,使得如何进行科学布料成为生产过程中重要的关键技术问题。
研究高炉布料分布规律,并建立相应的高炉布料数学模型,对于企业可以增加产量、节能、降低燃料消耗、提高煤气利用率、提高生铁质量、降低生铁成本、延长一代炉龄寿命等;对于社会可以减少CO2的排放量从而减少对环境的影响;必将产生巨大的经济效益和社会效益。
目前高炉配有专家系统,具有炉料控制、出铁管理等功能,但还没有炉料的料层实时分布状况进行预报模型,无法得到准确的矿焦分布,无法判断炉料下降、煤气气流分布及料层塌落等因素对料层分布的影响,因此,开发炉料分布模型,对维护高炉炉型、提高高炉合理布料、降低燃料比,都是非常必要的。关于炉料分布模型已有许多,模拟手段多采用几何方法,大都只针对炉喉区域半边料层进行模拟,对于整个炉体区域料层实时分布状况研究的很少,不能满足现场复杂炉料分布情况的数值模拟要求。北京首钢自动化提出的一种基于智能算法的高炉布料数值模拟方法(申请号:201210055516.9)提给出了炉喉区料层分布的数值模拟问题,该方法在计算中采用的三段法料面形状,该方法在远离炉喉区的位置与实际的偏差较大,而且对于某些炉况特殊时期的加料方法形成的料层形状的计算存在不足;北京科技大学的程树森提出了基于光学栅格的三维料面显示(无钟高炉布料测试新技术及料面三维图像重建)无法给出由炉喉至炉体上部的其他多层料面的形状。
发明内容
针对现有技术存在的不足,本发明提供一种高炉炉料分布实时预报系统及方法。
本发明的技术方案是:
一种高炉炉料分布实时预报系统,包括采集数据器、数据处理器、初始料面设定器、同档位下新料面检测器、多档位下新料面检测器、料面下降模块、料面分布信息生成器和布料控制器;
采集数据器:获取当前高炉生产过程工艺数据,包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
数据处理器:处理高炉生产过程工艺数据中的炉料重量及化验数据,包括缺失数据补齐和奇异数据校正;
初始料面设定器:采用N条线段组合而成的分段函数设定初始料面;
同档位下新料面检测器:根据获取的当前高炉生产过程工艺数据,构造相同布料档位下的新料面;
多档位下新料面检测器:形成多档位下新料面,即在不同布料档位下构造相同布料档位下的新料面,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
料面下降模块:构造下降料面,即当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
料面分布信息生成器:计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
布料控制器:根据高炉炉料分布实时预报的料层分布信息指导生产操作。
一种高炉炉料分布实时预报方法,包括以下步骤:
步骤1:获取当前高炉生产过程工艺数据,包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
步骤2:处理高炉生产过程工艺数据中的炉料重量及化验数据,包括缺失数据补齐和奇异数据校正;
缺失数据补齐采用均值的方法:计算前M次布料时该数据的平均值,并作为当前该数据值。
奇异数据校正采用上下界限定的方法:根据每项数据的物理特性和工艺要求,给定一个上界和下界,如果当前数据超出该项数据的上界与下界之间的范围,则采用前次布料时该数据的值进行校正。
步骤3:采用N条线段组合而成的分段函数设定初始料面;
步骤4:构造相同布料档位下的新料面;
步骤4.1:根据料流阀开度信息和原料粒度,构造炉料出料罐运动速度检测模型;
步骤4.2:根据炉料出料罐运动速度、Y形管与水平方向的夹角、Y形管斜面长度构造炉料出Y形管运动速度检测模型;
步骤4.3:根据炉料出Y形管运动速度、Y形管与水平方向的夹角、溜槽与水平方向的倾角、中心喉管长度构造炉料落入溜槽时的运动速度检测模型;
步骤4.4:将料流作为质点,根据炉料落入溜槽时的运动速度、溜槽旋转速度、溜槽与水平方向的倾角、溜槽有效长度、溜槽倾动距、构造炉料出溜槽运动速度检测模型;
步骤4.5:构造分料流速度检测模型:在相同布料档位下将料流划分为m股,各股料流出溜槽时水平速度和垂直速度不相同,各股料流与水平方向倾角相同且与溜槽与水平方向的倾角相等,根据溜槽有效长度、溜槽倾动距、溜槽与水平方向的倾角、溜槽固定点位置标高、炉喉平台位置标高计算溜槽末端的坐标,根据每股料流间距、溜槽与水平方向的倾角以及溜槽末端的坐标计算第i股料流出溜槽的坐标,根据炉料出溜槽运动速度计算各股料流运动速度;
步骤4.6:根据炉料空区运动时间和各股料流运动速度构造各股料流的空区运动轨迹模型;
步骤4.7:构造新料面:将各股料流的空区运动轨迹与上次布料料面的交点作为当前料流的初始落点,根据每股料流的初始落点、内堆角、外堆角及上次布料料面各个线段的斜率形成不同的堆,直到形成的堆所围成的体积等于每股炉料实际撒落的体积,即得到新形成料面,完成该档位下布料;
步骤5:形成多档位下新料面:在不同布料档位下重复步骤4,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
步骤6:构造下降料面:当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
在相同时间内,在保证炉喉及炉身上半部每层炉料下降过程中体积不变的前提下,探尺对应的料面点的位移与速度之比等于半径方向上该料面上其余各点的位移与速度之比;
步骤7:计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
步骤7.1:生成料面特征,包括料面漏斗宽度、料面平台宽度和料面漏斗深度;
步骤7.2:生成径向矿焦比:
步骤8:根据高炉炉料分布实时预报的料层分布信息指导生产操作;
步骤8.1:根据生成的径向矿焦比设定径向矿焦比的上限、下限;
步骤8.2:判断径向矿焦比是否超出设定的径向矿焦比的上限、下限:若边缘径向矿焦比低于设定的径向矿焦比的下限,则压边,即增加靠近高炉边缘档位的布料圈数;边缘径向矿焦比高于设定的径向矿焦比的上限,则松边,即减少靠近高炉边缘档位的布料圈数;漏斗径向矿焦比低于设定的径向矿焦比的下限,则抑制中心,即增加靠近高炉中心档位的布料圈数;漏斗径向矿焦比高于设定的径向矿焦比的上限,则疏导中心,即减少靠近高炉中心档位的布料圈数。
有益效果:
基于高炉实际生产流程的特点和要求,本发明提供了一种高炉炉料实时分布预报系统及方法,适用于高炉炉体区域复杂工艺过程中对炉料的料层分布实时状况进行预报。采用了机理模型与数据驱动相结合的方法,可以根据矿焦分布炉况的判断和布料模型的输出结果,对高炉矿焦分布进行合理调整,有效指导高炉布料操作,抑制炉况波动,降低炉况异常引起的休风率,提高煤气利用率和冶炼强度,减少焦比和综合焦比,改善高炉冶炼状况,提高料面透气性,改善煤气利用率,实现高炉生产稳定运行,从而达到促进炉况顺行、降低焦比的目的。使煤气与炉料进行充分接触,以最大可能地利用煤气的热能和化学能,为高炉能够优质、低耗、高产、长寿、稳定地生产提供有力的手段。
附图说明
图1为本发明具体实施方式的高炉炉料分布实时预报系统框图;
图2为本发明具体实施方式的高炉炉料分布实时预报方法流程图;
图3为本发明具体实施方式的构造相同布料档位下的新料面流程图;
图4为本发明具体实施方式的高炉炉料分布实时预报工艺图;
图5为本发明具体实施方式的初始料面设定图;
图6为本发明具体实施方式的炉料流宽度划分示意图;
图7为本发明料面信息结果图,(a)为按十字测温点位置计算矿焦比曲线,(b)为按等面积计算矿焦比曲线,(c)为按环位计算矿焦比曲线,(d)为漏斗的宽度和深度,(e)为平台宽度。
具体实施方式
下面结合附图对本发明的具体实施方式做详细说明。
本实施方式是针对一个大型钢铁企业的炼铁厂进行高炉炉料分布实时预报。
一种高炉炉料分布实时预报系统,如图1所示,包括采集数据器、数据处理器、初始料面设定器、同档位下新料面检测器、多档位下新料面检测器、料面下降模块、料面分布信息生成器和布料控制器;
采集数据器:获取当前高炉生产过程工艺数据,包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
数据处理器:处理高炉生产过程工艺数据中的炉料重量及化验数据,包括缺失数据补齐和奇异数据校正;
初始料面设定器:采用N条线段组合而成的分段函数设定初始料面;
同档位下新料面检测器:根据获取的当前高炉生产过程工艺数据,构造相同布料档位下的新料面;
多档位下新料面检测器:形成多档位下新料面,即在不同布料档位下构造相同布料档位下的新料面,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
料面下降模块:构造下降料面,即当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
料面分布信息生成器:计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
布料控制器:根据高炉炉料分布实时预报的料层分布信息指导生产操作。
如图4所示,DAU001-DAU003分别为料流阀开度信息、溜槽倾角及圈数信息、炉料重量及化验信息的数据采集器;DAU004为探尺料线深度信息的数据采集器;DPU001为数据处理器;PPD001为新料面检测器;FCM001为料面下降修正处理器;MPM001为料面信息分析器,用来计算料面特征信息及布料控制策略。高炉炉料分布实时预报工艺如图4所示。
一种高炉炉料分布实时预报方法,如图2所示,包括以下步骤:
步骤1:获取当前高炉生产过程工艺数据,包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
高炉布料制度数据如表1所示,炉料成分及比例数据如表2所示,高炉设备数据如表3所示;
表1 高炉布料制度(O/C:矿石/焦炭)
表2 炉料成分及比例
矿石成分 | 比例 | 焦炭成分 | 比例 |
四烧 | 75.69% | 中块 | 38.35% |
峨球 | 22.48% | 新焦炉 | 40.50% |
焦丁 | 1.84% | 离石 | 21.15% |
表3 高炉设备数据表
步骤2:处理高炉生产过程工艺数据中的炉料重量及化验数据,包括缺失数据补齐和奇异数据校正;
缺失数据补齐采用均值的方法:计算前M次布料时该数据的平均值,并作为当前该数据值。
奇异数据校正采用上下界限定的方法:根据每项数据的物理特性和工艺要求,给定一个上界和下界,如果当前数据超出该项数据的上界与下界之间的范围,则采用前次布料时该数据的值进行校正。
步骤3:采用N条线段组合而成的分段函数设定初始料面;
初始料面作为料面形状的起始,在经过一段时间布料过程之后,趋于最终稳定状态。模型应具有收敛性。
采用N条线段组合而成的分段函数来描述初始料面,如图5所示,通过设定A、B、C、D点的横纵坐标,线段BC与水平面夹角alfa,线段DE与水平面夹角beta,来设定多种初始料面。
料面可以表达为N条线段组合而成的分段函数,则当前要布料的料面函数y(x)表示如下形式:
其中,{(x1,y1),…,(xN,yN)}是分段函数的端点坐标集合,x1=0表示当前要布料的料面函数起始端点在炉中心线上,xN=r表示当前要布料的料面函数末端点在炉墙上。
步骤4:构造相同布料档位下的新料面,如图3所示;
步骤4.1:根据料流阀开度信息和原料粒度,构造炉料出料罐运动速度检测模型;
炉料出料罐的初速度v0:
其中:λ:原料系数;
D:料流阀开度(m);
d:原料的粒度(m);
设定原料为均匀混料,n种原料配比为p1,p2,……,pn;原料粒度分别为d1,d2,……,dn,则d=p1×d1+p2×d2+……pn×dn;求出不同料流阀开度下,炉料出料罐的初速度v0。
步骤4.2:根据炉料出料罐运动速度、Y形管与水平方向的夹角、Y形管斜面长度构造炉
料出Y形管运动速度检测模型;
炉料出Y形管运动速度v1:
其中:v0:炉料出料罐时的初速度;
αY:Y形管与水平方向的夹角;
lY:Y形管斜面长度;
μ:Y形管的摩擦系数;
步骤4.3:根据炉料出Y形管运动速度、Y形管与水平方向的夹角、溜槽与水平方向的
倾角、中心喉管长度构造炉料落入溜槽时的运动速度检测模型;
炉料以v1从Y形管流出经过中心喉管落到溜槽上。在此过程里,炉料可能会与中心喉管发生碰撞;在接触溜槽时也会与溜槽发生碰撞,从而有速度衰减和变化,假设炉料碰撞溜槽后,只存在沿溜槽方向的速度。炉料在此过程受到重力作用,由炉料运动方程可以求出炉料落入溜槽时速度v2。
其中:v1:炉料出Y形管时速度;
αY:Y形管与水平方向的夹角;
α:溜槽与水平方向的倾角;
H0:中心喉管长度;
b:溜槽倾动距;
k:速度衰减系数;
步骤4.4:将料流作为质点,根据炉料落入溜槽时的运动速度、溜槽旋转速度、溜槽与水平方向的倾角、溜槽有效长度、溜槽倾动距、构造炉料出溜槽运动速度检测模型;
炉料离开溜槽运动速度为v3:
其中:v2:炉料出中心喉管运动速度;
ω:溜槽旋转速度;
α:溜槽与水平方向的倾角;
lc:溜槽有效长度;
b:溜槽倾动距;
μ:溜槽的摩擦系数;
k:速度衰减系数;
炉料离开溜槽末端时,在各方向的分速度分别为:
步骤4.5:构造分料流速度检测模型:在相同布料档位下将料流划分为m股,如图6所示,各股料流出溜槽时水平速度和垂直速度不相同,各股料流与水平方向倾角相同且与溜槽与水平方向的倾角相等,根据溜槽有效长度、溜槽倾动距、溜槽与水平方向的倾角、溜槽固定点位置标高、炉喉平台位置标高计算溜槽末端的坐标,根据每股料流间距、溜槽与水平方向的倾角以及溜槽末端的坐标计算第i股料流出溜槽的坐标,根据炉料出溜槽运动速度计算各股料流运动速度;
模型考虑料流宽度对布料的影响,在相同档位下,将料流划分为m股。各料流出溜槽时档位相同,水平速度和垂直速度不相同。靠近炉喉中心处料流的速度慢,靠近炉墙料流速度快。第i股料流速度公式:
x0=lc·cosα-b·sinα
y0=(gc-gh)-lc·sinα-b·cosα
xi0=x0+i·dis·sinα
yi0=y0
其中:
gc:溜槽悬挂点标高;
gh:炉喉平台标高;
α:溜槽与水平方向的倾角;
lc:溜槽有效长度;
b:溜槽倾动距;
ai:实际数据驱动参数;
vsi:第i股料流的水平速度;
vhi:第i股料流的垂直速度;
dis:每股料流间距;
x0:溜槽末端的水平坐标;
y0:溜槽末端的垂直坐标;
xi0:第i股料流出溜槽的水平坐标;
yi0:第i股料流出溜槽的水平坐标;
步骤4.6:根据炉料空区运动时间和各股料流运动速度构造各股料流的空区运动轨迹模型;
炉料离开溜槽末端后的运动,在竖直方向上只受到重力作用。炉料粒子在只受重力的情况下走行轨迹为抛物线。以炉料空区运动时间为自变量的空区运动轨迹可表示为:
x=xi0+vhi·t
其中:
xi0:第i股料流出溜槽的水平坐标;
yi0:第i股料流出溜槽的水平坐标;
t:炉料空区运动时间;
vsi:第i股炉料在空区垂直方向的运动速度;
vhi:第i股炉料在空区水平方向的运动速度;
步骤4.7:构造新料面:将各股料流的空区运动轨迹与上次布料料面的交点作为当前料流的初始落点,根据每股料流的初始落点、内堆角、外堆角及上次布料料面各个线段的斜率形成不同的堆,直到形成的堆所围成的体积等于每股炉料实际撒落的体积,即得到新形成料面,完成该档位下布料;
每股料流在空区运动后落在旧料面上,落点会随着炉料的不断撒落不断改变,直到该股炉料全部撒落。初始落点定义为步骤4.6中求得的炉料轨迹最先接触旧料面时的交点,根据每股炉料初始落点、内堆角、外堆角及旧料面各个线段的斜率可形成不同的堆,随着时间的推进,形成的堆不断变化,直到形成的堆所围成的体积等于每股炉料实际撒落的体积。此时就可得到新形成料面,并完成该档位下布料。
步骤5:形成多档位下新料面:在不同布料档位下重复步骤4,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
步骤6:构造下降料面:当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
布料过程中料面实际上是不断连续下降的。并且由于受炉身角的影响,料面下降速度在炉喉半径方向上不一致:靠近炉中心线下降速度较慢;靠近炉墙下降速度较快。料面形状将在下降过程中趋于平坦。同时,料层厚度也由于炉身直径的增大而变薄。
根据体积守恒原则,料层厚度与炉身直径有关,料面点坐标与速度分布有关。当炉料在炉喉区域时,垂直向下运动,运动距离等于平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量。而平移量的计算根据体积守恒,由上至下逐层计算获得。其中,首层平移量根据探尺测量位置可知。
其中:s:料面向下平移量;
θ:炉身角;
在相同时间内,在保证炉喉及炉身上半部每层炉料下降过程中体积不变的前提下,探尺对应的料面点的位移与速度之比等于半径方向上该料面上其余各点的位移与速度之比;
根据料面在炉喉内的下降速度分布值,可计算炉料在下降过程中的平移量。下降速度可根据线性函数计算,也可根据历史数据拟合成二次函数计算。
其中:st:探尺向下平移量;
vt:探尺平均下降速度;
vr:当前位置点r的下降速度;
vr=a1′+a2′·r′-a3′·r′2 (11)
其中:r′:半径方向上离高炉中心的距离;
R0:炉喉半径;
R1:探尺与炉中心线的距离;
vp:炉料的平均下降速度;
ai:实际数据驱动参数;
步骤7:计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
步骤7.1:生成料面特征,包括料面漏斗宽度、料面平台宽度和料面漏斗深度;
料面特征是描述线性料面的特征。根据各个料面点的坐标值计算相关信息。料面漏斗宽度是料面内侧峰值点与炉中心线之间的距离。料面平台宽度是料面内侧峰值点与料面外侧峰值点之间的距离。料面漏斗深度是料面内侧峰值点与漏斗中心之间的高度差。
rwall ji=last(pji)
rcl ji=first(pji)
rterraceR ji=front(rwall ji) (12)
rterraceL ji=behind(rcl ji)
twidthji=rterraceR ji-rterraceL ji
fwidthji=rterraceL ji-rcl ji (13)
fheightji=zterraceL ji-zcl ji
其中:pji:料面点坐标集合;
rwall ji:料面炉墙坐标点水平坐标;
rcl ji:料面漏斗中心坐标点水平坐标;
rterraceR ji:料面内侧峰值点水平坐标;
rterraceL ji:料面外侧峰值点水平坐标;
twidthji:料面平台宽度;
fwidthji:料面中心漏斗宽度;
fheightji:料面中心漏斗深度;
步骤7.2:生成径向矿焦比:
模型中选用径向矿焦比衡量透气性。矿焦比定义如下:当最后几批炉料装入炉内后,在炉喉半径方向上各点处矿石层厚度与焦炭层厚度的比值与矿焦密度比的乘积,它反映了径向上的矿焦质量比例。
根据以档位划分的料面函数计算矿焦比值,可以对应十字测温点给出相应值。设最近n′批炉料,则可以求出矿焦比FO/C(r′):
其中:fo(r′)i:第i层矿石料面函数;
fc(r′)i:第i层焦炭料面函数;
ρo:矿石堆比重;
ρc:焦炭堆比重;
n′:定义炉料批数;
r′:半径方向上离高炉中心的距离;
计算如下三类径向矿焦比值:边缘、平台和漏斗区域的径向矿焦比值(均值);等面积区域(边缘、中间、中心)径向矿焦比值(均值);十字测温位置对应点的径向矿焦比值。
BorderO/C=Average(FO/C(r)) rterraceR≤r≤rwall
TerraceO/C=Average(FO/C(r)) rterraceL≤r≤rterraceR
FunnelO/C=Average(FO/C(r)) rcl≤r≤rterraceL
Temperature(i)O/C=FO/C(rTi) rTi∈{rT1,rT2,…,rT5}
本实施方式生成的料面特征如图7所示,其中(a)为按十字测温点位置计算矿焦比曲线,(b)为按等面积计算矿焦比曲线,(c)为按环位计算矿焦比曲线,(d)为漏斗的宽度和深度,(e)为平台宽度。
步骤8:根据高炉炉料分布实时预报的料层分布信息指导生产操作;
步骤8.1:根据生成的径向矿焦比设定径向矿焦比的上限、下限;
步骤8.2:判断径向矿焦比是否超出设定的径向矿焦比的上限、下限:若边缘径向矿焦比低于设定的径向矿焦比的下限,则压边,即增加靠近高炉边缘档位的布料圈数;边缘径向矿焦比高于设定的径向矿焦比的上限,则松边,即减少靠近高炉边缘档位的布料圈数;漏斗径向矿焦比低于设定的径向矿焦比的下限,则抑制中心,即增加靠近高炉中心档位的布料圈数;漏斗径向矿焦比高于设定的径向矿焦比的上限,则疏导中心,即减少靠近高炉中心档位的布料圈数。
Claims (5)
1.一种高炉炉料分布实时预报系统,其特征在于,包括采集数据器、数据处理器、初始料面设定器、同档位下新料面检测器、多档位下新料面检测器、料面下降模块、料面分布信息生成器和布料控制器;
采集数据器:获取当前高炉生产过程工艺数据后传至数据处理器,高炉生产过程工艺数据包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
数据处理器:对高炉生产过程工艺数据中的炉料重量及化验数据进行处理,包括缺失数据补齐和奇异数据校正;
初始料面设定器:根据数据处理器传来的数据,采用N条线段组合而成的分段函数设定初始料面;
同档位下新料面检测器:根据获取的当前高炉生产过程工艺数据,在设定的初始料面基础上构造相同布料档位下的新料面;
多档位下新料面检测器:在相同布料档位下的新料面基础上形成多档位下新料面,即在不同布料档位下构造相同布料档位下的新料面,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
料面下降模块:在构造多档位下新料面后构造下降料面,即当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
料面分布信息生成器:根据多档位下新料面检测器构造的多档位新料面和料面下降模块构造的下降料面,计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
布料控制器:根据料面分布信息生成器实时预报的料层分布信息指导生产操作。
2.采用权利要求1所述的高炉炉料分布实时预报系统的高炉炉料分布实时预报方法,其特征在于:包括以下步骤:
步骤1:获取当前高炉生产过程工艺数据,包括炉料重量及化验数据、料线信息、料流阀开度信息、高炉本体数据和布料矩阵;
炉料重量及化验数据包括:原料种类、原料体积、原料粒度、原料自然堆角和原料堆比重;
料线信息包括:探尺下降速度、矿石料线位置和焦炭料线位置;
料流阀开度信息包括:矿石料流阀开度和焦炭料流阀开度;
高炉本体数据包括:溜槽固定点位置标高、中心喉管长度、炉喉平台位置标高、炉喉半径、炉身角、Y形管与水平方向的夹角、Y形管斜面长度、探尺零点位置标高、探尺距炉中心距离、溜槽长度、溜槽有效长度、溜槽倾动距、溜槽旋转速度、溜槽摩擦系数和十字测温点水平坐标;
布料矩阵:布料档位和布料圈数,布料档位即溜槽与水平方向的倾角;
步骤2:处理高炉生产过程工艺数据中的炉料重量及化验数据,包括缺失数据补齐和奇异数据校正;
步骤3:设定初始料面;
步骤4:构造相同布料档位下的新料面;
步骤4.1:根据料流阀开度信息和原料粒度,构造炉料出料罐运动速度检测模型;
步骤4.2:根据炉料出料罐运动速度、Y形管与水平方向的夹角、Y形管斜面长度构造炉料出Y形管运动速度检测模型;
步骤4.3:根据炉料出Y形管运动速度、Y形管与水平方向的夹角、溜槽与水平方向的倾角、中心喉管长度构造炉料落入溜槽时的运动速度检测模型;
步骤 4.4:将料流作为质点,根据炉料落入溜槽时的运动速度、溜槽旋转速度、溜槽与水平方向的倾角、溜槽有效长度、溜槽倾动距、构造炉料出溜槽运动速度检测模型;
步骤4.5:构造分料流速度检测模型:在相同布料档位下将料流划分为m股,各股料流出溜槽时水平速度和垂直速度不相同,各股料流与水平方向倾角相同且与溜槽与水平方向的倾角相等,根据溜槽有效长度、溜槽倾动距、溜槽与水平方向的倾角、溜槽固定点位置标高、炉喉平台位置标高计算溜槽末端的坐标,根据每股料流间距、溜槽与水平方向的倾角以及溜槽末端的坐标计算第i股料流出溜槽的坐标,根据炉料出溜槽运动速度计算各股料流运动速度;
步骤4.6:根据炉料空区运动时间和各股料流运动速度构造各股料流的空区运动轨迹模型;
步骤4.7:构造新料面:将各股料流的空区运动轨迹与上次布料料面的交点作为当前料流的初始落点,根据每股料流的初始落点、内堆角、外堆角及上次布料料面各个线段的斜率形成不同的堆,直到形成的堆所围成的体积等于每股炉料实际撒落的体积,即得到新形成料面,完成该档位下布料;
步骤 5:形成多档位下新料面:在不同布料档位下重复步骤4,形成多档位下新料面,其中在进行档位间切换时溜槽倾角每增加1°,两档位之间的过渡圈内包含Q圈炉料体积,该过渡圈内的体积占用下一档位体积总值;
步骤 6:构造下降料面:当炉料在炉喉区域时,垂直向下运动,运动距离等于炉料平移量;当炉料在炉身区域时,炉料沿炉身角方向运动,水平运动距离等于平移量水平分量,垂直运动距离等于平移量垂直分量;
在相同时间内,在保证炉喉及炉身上半部每层炉料下降过程中体积不变的前提下,探尺对应的料面点的位移与速度之比等于半径方向上该料面上其余各点的位移与速度之比;
步骤7:计算料层分布信息并绘制料层分布图像,从而进行高炉炉料分布实时预报,料层分布信息包括径向矿焦比和料面特征信息;
步骤7.1:生成料面特征,包括料面漏斗宽度、料面平台宽度和料面漏斗深度;
步骤7.2:生成径向矿焦比:
步骤 8:根据高炉炉料分布实时预报的料层分布信息指导生产操作;
步骤8.1:根据生成的径向矿焦比设定径向矿焦比的上限、下限;
步骤8.2:判断径向矿焦比是否超出设定的径向矿焦比的上限、下限:若边缘径向矿焦比低于设定的径向矿焦比的下限,则压边,即增加靠近高炉边缘档位的布料圈数;边缘径向矿焦比高于设定的径向矿焦比的上限,则松边,即减少靠近高炉边缘档位的布料圈数;漏斗径向矿焦比低于设定的径向矿焦比的下限,则抑制中心,即增加靠近高炉中心档位的布料圈数;漏斗径向矿焦比高于设定的径向矿焦比的上限,则疏导中心,即减少靠近高炉中心档位的布料圈数。
3.根据权利要求2所述的高炉炉料分布实时预报方法,其特征在于:步骤3中采用N条线段组合而成的分段函数设定初始料面。
4.根据权利要求2所述的高炉炉料分布实时预报方法,其特征在于:缺失数据补齐采用均值的方法:计算前M次布料时该数据的平均值,并作为当前该数据值。
5.根据权利要求2所述的高炉炉料分布实时预报方法,其特征在于:奇异数据校正采用上下界限定的方法:根据每项数据的物理特性和工艺要求,给定一个上界和下界,如果当前数据超出该项数据的上界与下界之间的范围,则采用前次布料时该数据的值进行校正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410797862.1A CN104531924B (zh) | 2014-12-18 | 2014-12-18 | 一种高炉炉料分布实时预报系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410797862.1A CN104531924B (zh) | 2014-12-18 | 2014-12-18 | 一种高炉炉料分布实时预报系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104531924A CN104531924A (zh) | 2015-04-22 |
CN104531924B true CN104531924B (zh) | 2016-10-19 |
Family
ID=52847543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410797862.1A Active CN104531924B (zh) | 2014-12-18 | 2014-12-18 | 一种高炉炉料分布实时预报系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104531924B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105734191A (zh) * | 2014-12-10 | 2016-07-06 | 上海宝钢工业技术服务有限公司 | 万向节传动机构布料器的布料踪迹显示装置及方法 |
CN105243190B (zh) * | 2015-09-14 | 2018-08-07 | 内蒙古科技大学 | 一种高炉布料过程料面输出形状的建模方法 |
CN105950807B (zh) * | 2016-06-02 | 2018-10-16 | 燕山大学 | 一种多信息融合的高炉布料过程料面形状建模方法 |
CN106011350B (zh) * | 2016-07-05 | 2018-05-04 | 东北大学 | 一种高炉布料过程料面形状估计方法及系统 |
KR101858871B1 (ko) * | 2016-12-23 | 2018-06-27 | 주식회사 포스코 | 장입물 프로파일링 장치 |
TWI642789B (zh) * | 2017-04-25 | 2018-12-01 | 中國鋼鐵股份有限公司 | Method for obtaining distribution of coke ratio of material layer after blast furnace feeding |
CN107590333A (zh) * | 2017-09-07 | 2018-01-16 | 北京金恒博远科技股份有限公司 | 高炉内部料层的仿真方法及装置 |
CN107609261A (zh) * | 2017-09-07 | 2018-01-19 | 北京金恒博远科技股份有限公司 | 高炉布料过程的仿真方法、装置及系统 |
JP7176560B2 (ja) * | 2018-03-28 | 2022-11-22 | Jfeスチール株式会社 | 高炉操業方法 |
CN108470173B (zh) * | 2018-04-17 | 2021-06-15 | 北京矿冶科技集团有限公司 | 一种矿石颗粒的分割方法 |
CN113139275B (zh) * | 2021-03-22 | 2022-08-19 | 浙江大学 | 一种基于多层矿焦比分布模型的高炉炉喉温度估计方法 |
CN113549723B (zh) * | 2021-07-13 | 2022-09-20 | 张家港宏昌钢板有限公司 | 基于并罐高炉的偏斜检测及校正方法、装置及存储介质 |
CN113656966A (zh) * | 2021-08-17 | 2021-11-16 | 北京神网创新科技有限公司 | 一种高炉无钟炉顶在线布料模型仿真方法 |
TWI796041B (zh) * | 2021-12-10 | 2023-03-11 | 中國鋼鐵股份有限公司 | 高爐爐頂之料倉之秤重系統的動態壓力補償方法 |
CN115044719B (zh) * | 2022-06-13 | 2023-09-22 | 武汉钢铁有限公司 | 一种降料面休风判断料面位置的方法 |
CN115584370B (zh) * | 2022-09-01 | 2023-10-13 | 武安市裕华钢铁有限公司 | 高炉大渣量冶炼中改善软熔带透气性的方法 |
CN116254419A (zh) * | 2023-03-10 | 2023-06-13 | 安徽华铂再生资源科技有限公司 | 一种铅蓄电池回收用侧吹炉的出渣结构 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629286A (zh) * | 2012-02-24 | 2012-08-08 | 北京首钢自动化信息技术有限公司 | 一种基于智能算法的高炉布料数值模拟方法 |
-
2014
- 2014-12-18 CN CN201410797862.1A patent/CN104531924B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629286A (zh) * | 2012-02-24 | 2012-08-08 | 北京首钢自动化信息技术有限公司 | 一种基于智能算法的高炉布料数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
首钢炼铁高炉专家系统的开发与应用;周检平 等;《冶金自动化》;20081231(第S2期);530-534 * |
高炉无料钟炉顶布料料流轨迹检测与三维图像重建;杜鹏宇 等;《冶金自动化》;20091130;第33卷(第6期);1-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN104531924A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104531924B (zh) | 一种高炉炉料分布实时预报系统及方法 | |
CN104133945B (zh) | 一种高炉布料过程径向矿焦比的控制方法 | |
Li et al. | Numerical investigation of burden distribution in ironmaking blast furnace | |
CN106011350B (zh) | 一种高炉布料过程料面形状估计方法及系统 | |
CN102559963A (zh) | 一种高炉矿焦比分布调整的模型试验方法 | |
CN106957935B (zh) | 一种高炉内部软熔带形状的软测量方法 | |
CN107723399A (zh) | 一种高炉布料智能监测系统和调整方法 | |
CN108396136A (zh) | 一种烧结混合料制备方法及其布料方法与装置 | |
CN103981308A (zh) | 一种钒钛矿冶炼高炉的装料新方法 | |
KR101668584B1 (ko) | 고로로의 원료 장입 방법 | |
CN107245540B (zh) | 一种高炉布料过程径向料层厚度分布的控制策略 | |
CN103320151B (zh) | 一种试验用焦炉捣固机及其均密度控制方法 | |
CN114216349B (zh) | 一种基于编码解码网络的烧结终点预报方法 | |
Mitra | Modeling of burden distribution in the blast furnace | |
CN104302788B (zh) | 向高炉装入原料的方法 | |
Wu et al. | Optimization of coke ratio for the second proportioning phase in a sintering process base on a model of temperature field of material layer | |
CN209857663U (zh) | 一种烧结矿立式冷却机布料装置 | |
CN105975647A (zh) | 一种高耐磨铸造磨球成分设计与热处理制度的可视化方法 | |
CN102095526B (zh) | 一种基于烧结热量损失计算的环冷机烟气温度预测方法 | |
CN104215064B (zh) | 用于铁矿粉直接还原的直形隧道窑 | |
CN113656966A (zh) | 一种高炉无钟炉顶在线布料模型仿真方法 | |
Zhao et al. | The model of prediction of blast furnace gas output | |
CN206069960U (zh) | 一种烧结用混料装置 | |
US20240167111A1 (en) | Method for producing pig iron | |
Zhou et al. | Development of Physical Simulation Experiment Technology for Blast Furnace Coal Injection Ironmaking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |