CN113959628A - 高炉上部压力分布与炉况判定方法 - Google Patents

高炉上部压力分布与炉况判定方法 Download PDF

Info

Publication number
CN113959628A
CN113959628A CN202111226896.1A CN202111226896A CN113959628A CN 113959628 A CN113959628 A CN 113959628A CN 202111226896 A CN202111226896 A CN 202111226896A CN 113959628 A CN113959628 A CN 113959628A
Authority
CN
China
Prior art keywords
furnace
blast furnace
pressure
condition
blast
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
Application number
CN202111226896.1A
Other languages
English (en)
Other versions
CN113959628B (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.)
Pangang Group Panzhihua Steel and Vanadium Co Ltd
Pangang Group Panzhihua Iron and Steel Research Institute Co Ltd
Pangang Group Xichang Steel and Vanadium Co Ltd
Original Assignee
Pangang Group Research Institute Co Ltd
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 Pangang Group Research Institute Co Ltd filed Critical Pangang Group Research Institute Co Ltd
Priority to CN202111226896.1A priority Critical patent/CN113959628B/zh
Publication of CN113959628A publication Critical patent/CN113959628A/zh
Application granted granted Critical
Publication of CN113959628B publication Critical patent/CN113959628B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L11/00Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21BMANUFACTURE OF IRON OR STEEL
    • C21B5/00Making pig-iron in the blast furnace
    • C21B5/006Automatically controlling the process
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L13/00Devices or apparatus for measuring differences of two or more fluid pressure values
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/0028Force sensors associated with force applying means
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Abstract

本发明公开了一种高炉上部压力分布与炉况判定方法,属于冶金生产工艺技术领域。提供一种简单、高效的判定高炉炉况的高炉上部压力分布与炉况判定方法,所述的判定方法先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况。

Description

高炉上部压力分布与炉况判定方法
技术领域
本发明涉及一种高炉上部压力分布与炉况判定方法,属于冶金生产工艺技术领域。
背景技术
高炉是当前对铁氧化物还原最高效的工艺设备、炼钢工艺原材料的主要生产者,如在生产过程中出现问题,会延误和影响整个生产工艺的正常进行。但是由于在实际的生产过程中高炉内部是高温度、高压力、密闭复杂并且恶劣的环境,冶金过程很难进行实时的监测。因此,尽早的发现高炉炉况的异常波动并采取及时措施进行调整对高炉的稳定顺行至关重要。为此,国内外科研工作者开发了许多高炉专家系统综合判定高炉炉况,但是这些系统往往复杂,实用性并不高。同时高钛型高炉由于TiC、TiN、Ti(C,N)的存在,渣铁分离比普通高炉更为困难,受异常炉况的影响更大,而目前针对高钛型高炉压力分布与炉况判断的方法较少,因此研究开发了一种适用于对高钛型高炉的上部压力分布与炉况进行判定的方法便成为了本领域技术人员急需解决的技术问题。
发明内容
本发明所要解决的技术问题是:提供一种简单、高效的判定高炉炉况的高炉上部压力分布与炉况判定方法。
为解决上述技术问题所采用的技术方案是:一种高炉上部压力分布与炉况判定方法,所述的判定方法先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况,
若0.95Pc<Pcd<1.05Pc,则判定为高炉上部中心压力分布正常,炉况正常;并在此基础上通过比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势,
若Pcd≤0.95Pc或Pcd≥1.05Pc,则判定为高炉上部中心压力分布异常,高炉发生波动炉况异常,
其中,
Figure BDA0003314465740000011
式中Pc为临界压力,单位为kPa;
Figure BDA0003314465740000012
Figure BDA0003314465740000013
分别为同一时间下高炉风口鼓风压力与炉顶压力数据,单位为kPa,其中i为常数,取1、2、3、……;n为常数,表示所采集的数据总数量;
Pcd=Pw-Ptop,式中Pcd为高炉上部中心压差,单位为kPa;Pw和Ptop分别为高炉风口鼓风压力与炉顶压力数据,单位为kPa;
εmax=max(ε123,…,εm),式中
Figure BDA0003314465740000021
εm为第m层压力分布的均匀性,m取1、2、3、……,
Figure BDA0003314465740000022
为不同高度静压力均值,ΔPm为第m层静压力极差。
进一步的是,在计算高炉上部中心压差Pcd时,需要对采集的实时数据进行合理性判断,剔除不正常的数据。
上述方案的优选方式是,在计算不同高度静压力均值
Figure BDA0003314465740000023
时,需先对炉墙静压力数据进行分类编号,以位于同一高度平面的静压力数据记为同一层数据,沿风口平面向炉顶方向,以距离风口平面上方的第一层静压力数据记为编号“1”,第二层数据记为编号“2”,……,依次类推,靠近炉顶的最后一层数据记为编号“m”,然后计算不同高度层的炉墙静压力数据的均值,自风口向上到炉顶依次为
Figure BDA0003314465740000024
进一步的是,在计算第m层静压力极差ΔPm时,依据采集到的数据按ΔPm=Pmax-m-Pmin-m进行计算,式中Pmax-m和Pmin-m为第m层静压力数据的最大值和最小值,单位为kPa。
上述方案的优选方式是,临界波动值εc为高炉顺行期间一段时间的压力数据,按时间点分组依次代入公式ΔPm=Pmax-m-Pmin-m和εmax=max(ε123,…,εm),计算得到不同时间点下的压力分布均匀性指标最大值,记为εmax-1、εmax2、εmax-3、……、εmax-t,随后取上述数据的平均值并四舍五入取整后获得的。
进一步的是,计算出不同高度层的炉墙静压力数据的均值
Figure BDA0003314465740000025
后,将不同高度静压力均值
Figure BDA0003314465740000026
与对应的位置高度数据绘制成曲线并做平滑处理,寻找曲线变化的拐点位置S,并依据S的位置进行分类,
若拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况;
若拐点位置S位于炉腰,则高炉上部边缘压力获得发展,记为[B]类情况;
若拐点位置S位于炉身,则高炉上部边缘压力得到强化,记为[C]类情况,
其中,拐点位置S即为曲线上切线斜率最大的点所对应的位置。
上述方案的优选方式是,比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势时,
1)在[A]类情况下:
当εmax<1.5εc时,则可认为高炉上部压力分布正常,炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布出现异常,但由于边缘受到抑制,上部气流分布仍然表现为正常,且炉况仍可维持稳定;
当2.0εc≤εmax<2.5εc时,认为高炉上部边缘压力分布不均匀性加重,边缘气流分布出现波动,且炉况开始出现异常;
当εmax≥2.5εc时,认为高炉上部边缘压力分布非常紊乱,边缘气流分布大幅波动且炉况非常不稳定,急需调整;
2)在[B]或[C]类情况下:
当εmax<1.5εc时,认为高炉上部压力分布正常,中心正常,边缘发展或者强化,气流分布均匀,且炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布不均匀,边缘气流不稳定,且炉况开始出现异常;
Figure BDA0003314465740000031
时,认为高炉上部边缘压力分布非常紊乱,边缘气流失衡,且炉况非常不稳定,急需调整。
本发明的有益效果是:本申请提供的判断方法通过先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况,具体来说就是,若0.95Pc<Pcd<1.05Pc,则判定为高炉上部中心压力分布正常,炉况正常;并在此基础上通过比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势;若Pcd≤0.95Pc或Pcd≥1.05Pc,则判定为高炉上部中心压力分布异常,高炉发生波动炉况异常需要进行调整,达到简单、高效的判定高炉炉况的目的。
附图说明
图1为本发明高炉上部压力分布与炉况判定方法的判断逻辑图;
图2为本发明的判定方法涉及到的高炉上部不同高度压力平均值与对应高度的曲线图。
具体实施方式
如图1、图2所示是本发明提供的一种简单、高效的判定高炉炉况的高炉上部压力分布与炉况判定方法。所述的判定方法先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况,
若0.95Pc<Pcd<1.05Pc,则判定为高炉上部中心压力分布正常,炉况正常;并在此基础上通过比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势,
若Pcd≤0.95Pc或Pcd≥1.05Pc,则判定为高炉上部中心压力分布异常,高炉发生波动炉况异常,
其中,
Figure BDA0003314465740000041
式中Pc为临界压力,单位为kPa;
Figure BDA0003314465740000042
Figure BDA0003314465740000043
分别为同一时间下高炉风口鼓风压力与炉顶压力数据,单位为kPa,其中i为常数,取1、2、3、……;n为常数,表示所采集的数据总数量;
Pcd=Pw-Ptop,式中Pcd为高炉上部中心压差,单位为kPa;Pw和Ptop分别为高炉风口鼓风压力与炉顶压力数据,单位为kPa;
εmax=max(ε123,…,εm),式中
Figure BDA0003314465740000044
εm为第m层压力分布的均匀性,m取1、2、3、……,
Figure BDA0003314465740000045
为不同高度静压力均值,ΔPm为第m层静压力极差。本申请提供的判断方法通过先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况,具体来说就是,若0.95Pc<Pcd<1.05Pc,则判定为高炉上部中心压力分布正常,炉况正常;并在此基础上通过比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势;若Pcd≤0.95Pc或Pcd≥1.05Pc,则判定为高炉上部中心压力分布异常,高炉发生波动炉况异常需要进行调整,达到简单、高效的判定高炉炉况的目的。
上述实施方式中,为了提高判断准确度,并尽可能的避免出现判定的异常情况,本申请在计算高炉上部中心压差Pcd时,需要对采集的实时数据进行合理性判断,剔除不正常的数据。而且在计算不同高度静压力均值
Figure BDA0003314465740000046
时,需先对炉墙静压力数据进行分类编号,以位于同一高度平面的静压力数据记为同一层数据,沿风口平面向炉顶方向,以距离风口平面上方的第一层静压力数据记为编号“1”,第二层数据记为编号“2”,……,依次类推,靠近炉顶的最后一层数据记为编号“m”,然后计算不同高度层的炉墙静压力数据的均值,自风口向上到炉顶依次为
Figure BDA0003314465740000047
相应的,在计算第m层静压力极差ΔPm时,依据采集到的数据按ΔPm=Pmax-m-Pmin-m进行计算,式中Pmax-m和Pmin-m为第m层静压力数据的最大值和最小值,单位为kPa。
与此同时,本申请所述的临界波动值εc为高炉顺行期间一段时间的压力数据,按时间点分组依次代入公式ΔPm=Pmax-m-Pmin-m和εmax=max(ε123,…,εm),计算得到不同时间点下的压力分布均匀性指标最大值,记为εmax-1、εmax2、εmax-3、……、εmax-t,随后取上述数据的平均值并四舍五入取整后获得的。并且在计算出不同高度层的炉墙静压力数据的均值
Figure BDA0003314465740000051
后,将不同高度静压力均值
Figure BDA0003314465740000052
与对应的位置高度数据绘制成曲线并做平滑处理,寻找曲线变化的拐点位置S,并依据S的位置进行分类,
若拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况;
若拐点位置S位于炉腰,则高炉上部边缘压力获得发展,记为[B]类情况;
若拐点位置S位于炉身,则高炉上部边缘压力得到强化,记为[C]类情况,
其中,拐点位置S即为曲线上切线斜率最大的点所对应的位置。
比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势时,
1)在[A]类情况下:
当εmax<1.5εc时,则可认为高炉上部压力分布正常,炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布出现异常,但由于边缘受到抑制,上部气流分布仍然表现为正常,且炉况仍可维持稳定;
当2.0εc≤εmax<2.5εc时,认为高炉上部边缘压力分布不均匀性加重,边缘气流分布出现波动,且炉况开始出现异常;
当εmax≥2.5εc时,认为高炉上部边缘压力分布非常紊乱,边缘气流分布大幅波动且炉况非常不稳定,急需调整;
2)在[B]或[C]类情况下:
当εmax<1.5εc时,认为高炉上部压力分布正常,中心正常,边缘发展或者强化,气流分布均匀,且炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布不均匀,边缘气流不稳定,且炉况开始出现异常;
Figure BDA0003314465740000053
时,认为高炉上部边缘压力分布非常紊乱,边缘气流失衡,且炉况非常不稳定,急需调整。
综上所述,本申请提供的技术方案,仅利用高炉上部压力数据经过简单计算就能较为准确地判断高炉上部压力分布情况及炉况,简单、高效、实用,且适用于高钛型高炉。
具体实施例
本发明的目的在于,提供一种用于高炉上部压力分布与炉况判定方法,该判定方法基于高炉上部(风口以上,包括炉腹、炉腰、炉身)的侧壁静压力数据及顶压、风压数据即可快速判定普通矿冶炼高炉(或钒钛磁铁矿冶炼高炉)的上部压力分布情况及高炉炉况,达到简单、高效判定高炉炉况的目的。
技术方案:
为实现以上目的,本发明通过以下步骤判定高炉上部压力分布与炉况:
步骤一:计算判定所需临界压差
利用所判定高炉在炉况稳定顺行条件下一段时间(至少1个月以上)的高炉压力数据,包括炉顶压力、鼓风压力与高炉上部炉墙侧壁静压力数据。首先计算临界压差,计算公式如下:
Figure BDA0003314465740000061
式中:Pc为临界压力,kPa;
Figure BDA0003314465740000062
Figure BDA0003314465740000063
分别为同一时间下高炉风口鼓风压力与炉顶压力数据,kPa,其中i为常数,取1、2、3、……;n为常数,表示所采集的数据总数量。
步骤二:计算高炉上部中心压差
采集计算所需的实时数据(高炉上部压力数据)并判断其合理性,剔除不正常数据。首先,计算高炉上部中心压差,计算公式如下:
Pcd=Pw-Ptop (2)
式中:Pcd为高炉上部中心压差,kPa;Pw和Ptop分别为高炉风口鼓风压力与炉顶压力数据,kPa。
计算高炉上部边缘压力分布
首先,对炉墙静压力数据进行分类编号,以位于同一高度平面的静压力数据记为同一层数据,沿风口平面向炉顶方向,以距离风口平面上方的第一层静压力数据记为编号“1”,第二层数据记为编号“2”,……,依次类推,靠近炉顶的最后一层数据记为编号“m”。
计算不同高度层的炉墙静压力数据的均值,自风口向上到炉顶依次为
Figure BDA0003314465740000064
Figure BDA0003314465740000065
然后将不同高度静压力均值
Figure BDA0003314465740000066
与对应的位置高度数据,绘制成曲线并做平滑处理,寻找曲线变化的拐点位置S,即曲线上切线斜率最大的点所对应的位置(高度及层数编号)。
若拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况;
若拐点位置S位于炉腰,则高炉上部边缘压力获得发展,记为[B]类情况;
若拐点位置S位于炉身,则高炉上部边缘压力得到强化,记为[C]类情况。
步骤三:计算压力分布均匀性
根据表一数据,计算每层炉墙静压力数据的极差值ΔPm,计算公式如下:
ΔPm=Pmax-m-Pmin-m (3)
式中:ΔPm为第m层静压力极差,kPa;Pmax-m和Pmin-m为第m层静压力数据的最大值和最小值,kPa。
然后,计算每层压力分布的均匀性指标εm,计算公式如下:
Figure BDA0003314465740000071
式中:εm为第m层压力分布的均匀性,其中m取1、2、3、……。
计算压力分布均匀性指标的最大值εmax,即取ε1、ε2、ε3、……、εm的最大值,并记录最大值对应的层数编号mmax,以便查找问题所在区域,计算公式如下:
εmax=max(ε123,…,εm) (5)
根据步骤一的高炉顺行期间一段时间的压力数据,按时间点分组依次代入公式(3)、(4)和(5),计算得到不同时间点下的压力分布均匀性指标最大值,记为εmax-1、εmax2、εmax-3、……、εmax-t,随后取上述数据的平均值并四舍五入取整后设为临界波动值εc
步骤四:判断高炉炉况
1)若0.95Pc<Pcd<1.05Pc,则可认为高炉上部中心压力分布正常,进一步判断炉况,如下:
(1)在[A]类情况下:
当εmax<1.5εc时,则可认为高炉上部压力分布正常,炉况稳定。
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布出现异常,但由于边缘受到抑制,上部气流分布仍然表现为正常,且炉况仍可维持稳定;
当2.0εc≤εmax<2.5εc时,认为高炉上部边缘压力分布不均匀性加重,边缘气流分布出现波动,且炉况开始出现异常;
当εmax≥2.5εc时,认为高炉上部边缘压力分布非常紊乱,边缘气流分布大幅波动且炉况非常不稳定,急需调整。
(2)在[B]或[C]类情况下:
当εmax<1.5εc时,认为高炉上部压力分布正常,中心正常,边缘发展或者强化,气流分布均匀,且炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布不均匀,边缘气流不稳定,且炉况开始出现异常;
当εmax≥2.0εc时,认为高炉上部边缘压力分布非常紊乱,边缘气流失衡,且炉况非常不稳定,急需调整。
2)若Pcd≤0.95Pc或Pcd≥1.05Pc,则可认为高炉上部中心压力分布异常,高炉发生波动,炉况异常。
实施例一
为让本发明的上述特征和优点能更加明显易懂,结合附图1判定逻辑图,以某钢铁公司高炉上部压力采集数据为例,对本发明作进一步详细说明:
图1为该发明的判定逻辑图;
图2为某高炉的上部不同高度静压力数据与对应高度的曲线图。
步骤一:计算判定所需临界压差
利用所判定高炉在炉况稳定顺行条件下一段时间(至少1个月以上)的高炉压力数据,包括炉顶压力、鼓风压力与高炉上部炉墙侧壁静压力数据。首先计算临界压差,计算公式如下:
Figure BDA0003314465740000081
式中:Pc为临界压力,kPa;
Figure BDA0003314465740000082
Figure BDA0003314465740000083
分别为同一时间下高炉风口鼓风压力与炉顶压力数据,kPa,其中i为常数,取1、2、3、……;n为常数,表示用于判定的数据总数量。基于所判定高炉近一个月的压力数据,代入式(1)计算得到:Pc=155.7kPa。
步骤二:计算高炉上部中心压差
采集计算所需的实时数据(高炉上部压力数据)并判断其合理性,剔除不正常数据,如表一。首先,计算高炉上部中心压差,计算公式如下:
Pcd=Pw-Ptop (2)
式中:Pcd为高炉上部中心压差,kPa;Pw和Ptop分别为高炉风口鼓风压力与炉顶压力数据,kPa。代入表一中数据,计算得到Pcd=160.5kPa。
表一某时间点高炉上部压力实时采集数据
Figure BDA0003314465740000084
Figure BDA0003314465740000091
步骤二:计算高炉上部边缘压力分布
首先,对炉墙静压力数据进行分类编号,以位于同一高度平面的静压力数据记为同一层数据,沿风口平面向炉顶方向,以距离风口平面上方的第一层静压力数据记为编号“1”,第二层数据记为编号“2”,……,依次类推,靠近炉顶的最后一层数据记为编号“m”。此高炉m=8,即共8层。
根据表一数据,计算不同高度层的炉墙静压力数据的均值,自风口向上到炉顶依次为
Figure BDA0003314465740000092
Figure BDA0003314465740000093
然后将不同高度静压力均值
Figure BDA0003314465740000094
与对应的位置高度数据,绘制成曲线并做平滑处理,寻找曲线变化的拐点位置S,即曲线上切线斜率最大的点所对应的位置(高度),如图2所示,计算得到S≈18.6m。
由于该高炉炉腰的高度区间为15.91~18.85m,则拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况。
步骤三:计算压力分布均匀性
根据表一数据,计算每层炉墙静压力数据的极差值ΔPm,计算公式如下:
ΔPm=Pmax-m-Pmin-m (3)
式中:ΔPm为第m层静压力极差,kPa;Pmax-m和Pmin-m为第m层静压力数据的最大值和最小值,kPa。
然后,计算每层压力分布的均匀性指标εm,计算公式如下:
Figure BDA0003314465740000095
式中:εm为第m层压力分布的均匀性,其中m取1、2、3、……。
计算压力分布均匀性指标的最大值εmax,即取ε1、ε2、ε3、……、εm中的最大值,代入数据到式(5)得到εmax=15.3%,并记录最大值对应的层数编号mmax=3,属于高炉炉腰区域,计算公式如下:
εmax=max(ε123,…,εm) (5)
根据步骤一的高炉顺行期间一段时间的压力数据,按时间点分组依次代入公式(3)、(4)和(5),计算得到不同时间点下的压力分布均匀性指标最大值,记为εmax-1,εmax2,εmax-3,……,εmax-t,随后取上述数据的平均值并四舍五入取整后设为临界波动值εc,经计算εc=10.0%。
步骤四:判断高炉炉况
1)由于Pcd=160.5kPa,Pc=155.7kPa,则0.95Pc<Pcd<1.05Pc,则可认为高炉上部中心压力分布正常,进一步判断炉况,如下:
根据步骤二判断可知拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况。在此条件下,由于εmax=15.3%,εc=10.0%,则1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布出现异常,但由于边缘受到抑制,上部气流分布仍然表现为正常,且炉况仍可维持稳定,但需要现场操作者持续关注炉况,避免恶化。

Claims (7)

1.一种高炉上部压力分布与炉况判定方法,其特征在于:所述的判定方法先采集炉况稳定顺行条件下的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得炉况稳定顺行条件下的临界压差Pc和临界波动值εc,然后采集需要判定炉况高炉实际运行过程中的炉顶压力值、鼓风压力值以及高炉上部炉墙侧壁静压力值,并计算获得该高炉运行时的高炉上部中心压差Pcd和压力分布均匀性指标的最大值εmax,最后将Pcd与Pc进行比较来判断高炉运行时的炉内情况,
若0.95Pc<Pcd<1.05Pc,则判定为高炉上部中心压力分布正常,炉况正常;并在此基础上通过比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势,
若Pcd≤0.95Pc或Pcd≥1.05Pc,则判定为高炉上部中心压力分布异常,高炉发生波动炉况异常,
其中,
Figure FDA0003314465730000011
式中Pc为临界压力,单位为kPa;
Figure FDA0003314465730000012
Figure FDA0003314465730000013
分别为同一时间下高炉风口鼓风压力与炉顶压力数据,单位为kPa,其中i为常数,取1、2、3、……;n为常数,表示所采集的数据总数量;
Pcd=Pw-Ptop,式中Pcd为高炉上部中心压差,单位为kPa;Pw和Ptop分别为高炉风口鼓风压力与炉顶压力数据,单位为kPa;
εmax=max(ε123,…,εm),式中
Figure FDA0003314465730000014
εm为第m层压力分布的均匀性,m取
Figure FDA0003314465730000015
为不同高度静压力均值,ΔPm为第m层静压力极差。
2.根据权利要求1所述的高炉上部压力分布与炉况判定方法,其特征在于:在计算高炉上部中心压差Pcd时,需要对采集的实时数据进行合理性判断,剔除不正常的数据。
3.根据权利要求1所述的高炉上部压力分布与炉况判定方法,其特征在于:在计算不同高度静压力均值
Figure FDA0003314465730000016
时,需先对炉墙静压力数据进行分类编号,以位于同一高度平面的静压力数据记为同一层数据,沿风口平面向炉顶方向,以距离风口平面上方的第一层静压力数据记为编号“1”,第二层数据记为编号“2”,……,依次类推,靠近炉顶的最后一层数据记为编号“m”,然后计算不同高度层的炉墙静压力数据的均值,自风口向上到炉顶依次为
Figure FDA0003314465730000017
Figure FDA0003314465730000018
4.根据权利要求3所述的高炉上部压力分布与炉况判定方法,其特征在于:在计算第m层静压力极差ΔPm时,依据采集到的数据按ΔPm=Pmax-m—Pmin-m进行计算,式中Pmax-m和Pmin-m为第m层静压力数据的最大值和最小值,单位为kPa。
5.根据权利要求4所述的高炉上部压力分布与炉况判定方法,其特征在于:临界波动值εc为高炉顺行期间一段时间的压力数据,按时间点分组依次代入公式ΔPm=Pmax-m—Pmin-m和εmax=max(ε123,…,εm),计算得到不同时间点下的压力分布均匀性指标最大值,记为εmax-1、εmax2、εmax-3、……、εmax-t,随后取上述数据的平均值并四舍五入取整后获得的。
6.根据权利要求3、4、5所述的高炉上部压力分布与炉况判定方法,其特征在于:计算出不同高度层的炉墙静压力数据的均值
Figure FDA0003314465730000021
后,将不同高度静压力均值
Figure FDA0003314465730000022
与对应的位置高度数据绘制成曲线并做平滑处理,寻找曲线变化的拐点位置S,并依据S的位置进行分类,
若拐点位置S位于炉腹,则高炉上部边缘压力分布受到抑制,记为[A]类情况;
若拐点位置S位于炉腰,则高炉上部边缘压力获得发展,记为[B]类情况;
若拐点位置S位于炉身,则高炉上部边缘压力得到强化,记为[C]类情况,
其中,拐点位置S即为曲线上切线斜率最大的点所对应的位置。
7.根据权利要求6所述的高炉上部压力分布与炉况判定方法,其特征在于:比较εc与εmax的关系判高炉炉况在不同类型情况下的发展趋势时,
1)在[A]类情况下:
当εmax<1.5εc时,则可认为高炉上部压力分布正常,炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布出现异常,但由于边缘受到抑制,上部气流分布仍然表现为正常,且炉况仍可维持稳定;
当2.0εc≤εmax<2.5εc时,认为高炉上部边缘压力分布不均匀性加重,边缘气流分布出现波动,且炉况开始出现异常;
当εmax≥2.5εc时,认为高炉上部边缘压力分布非常紊乱,边缘气流分布大幅波动且炉况非常不稳定,急需调整;
2)在[B]或[C]类情况下:
当εmax<1.5εc时,认为高炉上部压力分布正常,中心正常,边缘发展或者强化,气流分布均匀,且炉况稳定;
当1.5εc≤εmax<2.0εc时,认为高炉上部边缘压力分布不均匀,边缘气流不稳定,且炉况开始出现异常;
Figure FDA0003314465730000023
时,认为高炉上部边缘压力分布非常紊乱,边缘气流失衡,且炉况非常不稳定,急需调整。
CN202111226896.1A 2021-10-21 2021-10-21 高炉上部压力分布与炉况判定方法 Active CN113959628B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111226896.1A CN113959628B (zh) 2021-10-21 2021-10-21 高炉上部压力分布与炉况判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111226896.1A CN113959628B (zh) 2021-10-21 2021-10-21 高炉上部压力分布与炉况判定方法

Publications (2)

Publication Number Publication Date
CN113959628A true CN113959628A (zh) 2022-01-21
CN113959628B CN113959628B (zh) 2023-06-30

Family

ID=79465916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111226896.1A Active CN113959628B (zh) 2021-10-21 2021-10-21 高炉上部压力分布与炉况判定方法

Country Status (1)

Country Link
CN (1) CN113959628B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10219317A (ja) * 1997-01-31 1998-08-18 Kawasaki Steel Corp 高炉操業方法
CN1575351A (zh) * 2001-08-24 2005-02-02 阿维扎技术公司 含有内压控制系统的常压晶片加工反应器及方法
CN102559964A (zh) * 2012-02-10 2012-07-11 中冶南方工程技术有限公司 高炉炉顶压力前馈控制方法
JP2015086461A (ja) * 2013-11-01 2015-05-07 新日鐵住金株式会社 高炉操業方法
CN111349793A (zh) * 2020-04-16 2020-06-30 沈阳东大山汇环境科技有限公司 一种高炉或竖炉的锌回收方法及其装置
CN112342327A (zh) * 2020-10-09 2021-02-09 攀钢集团研究院有限公司 一种基于理论燃烧温度控制的钒钛磁铁矿高炉冶炼方法
CN112668148A (zh) * 2020-12-04 2021-04-16 攀钢集团研究院有限公司 一种高钛型高炉上部气流分布情况与炉况判定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10219317A (ja) * 1997-01-31 1998-08-18 Kawasaki Steel Corp 高炉操業方法
CN1575351A (zh) * 2001-08-24 2005-02-02 阿维扎技术公司 含有内压控制系统的常压晶片加工反应器及方法
CN102559964A (zh) * 2012-02-10 2012-07-11 中冶南方工程技术有限公司 高炉炉顶压力前馈控制方法
JP2015086461A (ja) * 2013-11-01 2015-05-07 新日鐵住金株式会社 高炉操業方法
CN111349793A (zh) * 2020-04-16 2020-06-30 沈阳东大山汇环境科技有限公司 一种高炉或竖炉的锌回收方法及其装置
CN112342327A (zh) * 2020-10-09 2021-02-09 攀钢集团研究院有限公司 一种基于理论燃烧温度控制的钒钛磁铁矿高炉冶炼方法
CN112668148A (zh) * 2020-12-04 2021-04-16 攀钢集团研究院有限公司 一种高钛型高炉上部气流分布情况与炉况判定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"钢铁研究学报 2019年 第31卷 总目次", 钢铁研究学报, no. 12 *
吴敏;聂秋平;许永华;曹卫华;: "高炉煤气流分布性能可拓评价方法研究及其应用", 中南大学学报(自然科学版), no. 03 *
谢洪恩;秦兴国;郑魁;蒋胜;: "高钛型高炉渣熔化性温度影响因素分析", 中国冶金, no. 09 *

Also Published As

Publication number Publication date
CN113959628B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
CN104630410B (zh) 一种基于数据解析的转炉炼钢质量实时动态预测方法
CN111593155B (zh) 高炉诊断系统及方法
CN111607673B (zh) 铁前诊断系统及方法
CN103439999B (zh) 根据冷却壁温度变化进行高炉异常炉温控制的方法
CN112668148B (zh) 一种高钛型高炉上部气流分布情况与炉况判定方法
CN106521059B (zh) 用相控阵雷达测量高炉料面矿焦比来控制高炉气流分布的方法
CN111639801B (zh) 高炉炉况的评分方法及系统
CN103544273A (zh) 利用模式识别技术对炉况整体状态进行评估的方法
CN108647407A (zh) 一种转炉炼钢烟气分析定碳方法
CN102925602B (zh) 高炉操作炉型维护方法
CN110699503B (zh) 高炉炉缸状态调整策略确定方法、装置及电子设备
CN102758032B (zh) 一种实时预测高炉发生管道概率的方法
CN113959628A (zh) 高炉上部压力分布与炉况判定方法
CN110835661A (zh) 一种高炉操作炉型的判断方法
CN109022650B (zh) 一种高炉上部煤气流分布的识别方法
CN113637814B (zh) 高炉中心气流波动的调控方法
KR101225270B1 (ko) 연화융착대 근부위치 추정 방법
CN109491421B (zh) 基于冶金二次资源循环利用的智能温度控制系统的温度控制方法
CN106048117A (zh) 炼铁高炉判断炉温的简易方法
CN107907435B (zh) 一种轧辊辊身硬度分析方法
CN111286570A (zh) 利用扫描雷达调控异常操作炉型的方法
CN113265498A (zh) 一种高炉炉型管控方法
CN111026207A (zh) 一种高炉循环冷却水流量调节控制方法
CN115655999A (zh) 一种基于炉身静压力判断炉缸死料柱状态的方法及系统
CN113325811B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231222

Address after: 617000 Taoyuan street, East District, Panzhihua, Sichuan Province, No. 90

Patentee after: PANGANG GROUP PANZHIHUA IRON & STEEL RESEARCH INSTITUTE Co.,Ltd.

Patentee after: PANGANG GROUP PANZHIHUA STEEL & VANADIUM Co.,Ltd.

Patentee after: XICHANG STEEL VANADIUM CO.,LTD., PANGANG Group

Address before: 610000 Panzhihua Iron and Steel Group Research Institute Co.,Ltd.

Patentee before: PANGANG GROUP RESEARCH INSTITUTE Co.,Ltd.