CN113283079A - 一种高炉风口回旋区边界的计算及实时监测方法 - Google Patents
一种高炉风口回旋区边界的计算及实时监测方法 Download PDFInfo
- Publication number
- CN113283079A CN113283079A CN202110566801.4A CN202110566801A CN113283079A CN 113283079 A CN113283079 A CN 113283079A CN 202110566801 A CN202110566801 A CN 202110566801A CN 113283079 A CN113283079 A CN 113283079A
- Authority
- CN
- China
- Prior art keywords
- boundary
- convolution
- area
- blast
- depth
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000012544 monitoring process Methods 0.000 title claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 34
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 4
- 239000000571 coke Substances 0.000 claims description 51
- 238000011160 research Methods 0.000 claims description 15
- 238000007664 blowing Methods 0.000 claims description 14
- 239000002245 particle Substances 0.000 claims description 14
- 238000013178 mathematical model Methods 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 7
- 230000033001 locomotion Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 9
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 13
- 230000014509 gene expression Effects 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000003723 Smelting Methods 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Manufacture Of Iron (AREA)
- Blast Furnaces (AREA)
Abstract
本发明提供一种高炉风口回旋区边界的计算及实时监测方法,涉及高炉炼铁工艺技术领域。该方法首先根据高炉风口回旋区的形成原理,建立回旋区的深度计算模型,进而得到回旋区深度的计算公式,获得回旋区深度的变化规律;再通过高炉风口回旋区的深度模型建立高炉风口回旋区的边界模型,确定回旋区边界的计算公式;然后获得建模参数,分析建模参数对回旋区边界模型的影响,确定影响回旋区边界的主要参数;最后利用回旋区边界计算公式求出回旋区的高度;当回旋区高度或回旋区深度超出设定范围时,通过调节鼓风风压和鼓风风量使回旋区高度或深度恢复至正常范围内。该方法能够实时监测回旋区深度和回旋区边界的变化情况,为高炉的实际生产提供安全指导。
Description
技术领域
本发明涉及高炉炼铁工艺技术领域,尤其涉及一种高炉风口回旋区边界的计算及实时监测方法。
背景技术
在高炉炼铁生产中,通过高炉风口高温高速的空气被鼓入高炉内,由于鼓风的作用,在风口前沿附近形成一个焦炭在其内作循环运动的区域即为高炉风口回旋区。其中风口回旋区位于竖炉炉体下部风口前端,是由焦炭和喷入炉内的辅助燃料与鼓风中的氧进行剧烈的燃烧反应形成的。鼓风和煤气的混合气流在此区域内循环运动,同时伴随着细小焦炭颗粒和未燃尽煤粉的高速旋转,以及碎焦在回旋过程中的烧尽过程。回旋区大小与鼓风参数及原燃料条件等因素有直接关系,是热能和气体还原剂的发源地,为整个高炉生产提供热量和能量的补给,风口回旋区的深度和内部复杂的物理、化学反应决定了高炉中煤气流的一次分布及上部炉料的下降状态,反映了焦炭的燃烧状态,是炉况顺行的基础,在冶炼过程中起着至关重要的作用。
高炉风口回旋区的特征的研究主要可以分为两大方面:回旋区特征的直接法研究和回旋区特征的间接法研究。一是高炉风口回旋区特征的直接研究法是通过对表示高炉回旋区放的相关参数的直接检测来进行研究,主要集中在对回旋区的大小、形状及温度等参数的直接测量,但是存在仪表设备易受炉内实际环境影响而导致测量结果波动较大,同时仪器成本较高,而且无法达到实时监测的目的,无法在中小企业中完全普及。直接研究法又分为经验观察法研究和实测法研究;二是高炉回旋区特征的间接研究法即模型研究法,包括以下两个方面:其一是通过建立高炉风口回旋区的物理参数实验模型,针对高炉风口回旋区特征,在模型上进行实验检测,但是由于回旋区内部存在但由于回旋区内部的反应复杂多变,冷态模型不能很好的反映实际的回旋区内部状态;较为常用的方法是依据回旋区运动过程中动量、质量和热量的传输建立欧拉数学模型求解,但是采用现有的欧拉模型建模过程复杂,需要的参数较多,计算困难,花费时间较长,难以实现实时监测的目的的反应,实验模型不能很好的反映实际的回旋区内部状态,较为常用的方法是依据回旋区运动过程中动量、质量和热量的传输建立欧拉数学模型求解,但是采用现有的欧拉模型建模过程复杂,需要的参数较多,计算困难,花费时间较长,难以实现实时监测的目的;其二是利用建立了的高炉风口回旋区的二维或者三维数学模型,对回旋区内的化学反应过程进行数值模拟,从而达到研究回旋区的特征变化规律的目的。
发明内容
本发明要解决的技术问题是针对上述直接测量法和实验模型法技术存在的不足并在对机理数学模型上进行改进,提供一种高炉风口回旋区边界的计算及实时监测方法,在回旋区深度模型的基础上通过在回旋区边界处各点建立二力的平衡方程,能够高效实时的求解回旋区的边界的变化情况,获得回旋区边界的变化规律,研究回旋区内部参数对回旋区深度以及高度的影响,以及通过控制鼓风参数调节回旋区的深度以及高度的变化,为高炉的稳定运行提供可靠保障。
为解决上述技术问题,本发明所采取的技术方案是:一种高炉风口回旋区边界的计算及实时监测方法,包括以下步骤:
步骤1、根据高炉风口回旋区的形成原理,建立回旋区的深度计算模型,进而得到回旋区深度的计算公式,获得回旋区深度的变化规律;
当回旋区空穴处于稳定运动状态时,取回旋区内部最深处一微元区域A为研究对象,此时A处在鼓风气体冲力与焦炭层阻力共同作用下达到平衡,依据二力平衡建立回旋区的深度计算模型,用以求解回旋区深度的变化情况,如下公式所示:
其中,FA表示区域A处鼓风气体冲力,FB表示区域A处焦炭层阻力,ρg0表示标准态下的鼓风密度,Vg表示区域A处鼓风风量,ST表示风口面积,DR表示回旋区深度,DT表示风口直径,α为常数,用于表示回旋区深度与宽度之间的关系,P表示鼓风风压,Tm表示回旋区温度,SP表示区域A处焦炭粒子总截面积,ρP表示焦炭密度,VP表示区域A处焦炭粒子体积,g表示重力加速度,DPR表示回旋区边界前焦炭直径,DPR=0.6DPc,DPc表示入炉前焦炭直径,Tw表示鼓风温度,P表示鼓风风压;
由上述回旋区深度的计算模型得到回旋区深度的计算公式为:
其中,K、β均为待定系数,ρg表示实际风温风压下的鼓风密度;
步骤2、通过高炉风口回旋区的深度计算模型建立高炉风口回旋区的边界模型,确定回旋区边界的计算公式;
当回旋区内部运动处于稳定状态时,任意取回旋区边界点B为研究对象,此时边界点B处在鼓风气体冲力与焦炭层阻力共同作用下达到平衡,依据二力平衡建立回旋区边界处任意一点的数学模型,用以求解回旋区边界的变化情况:
其中,
MB∝ρg·Vg·(1-σL)
SB∝(1-σ((DR-x)a+hb))·SA
其中,FD表示边界点B处鼓风气体冲力,Fb表示边界点B处焦炭层对鼓风气流的阻力,MB表示B处鼓风质量流率,UB代表B处风速,SB为B处的截面积,SA表示区域A处的截面积,ρg表示实际风温风压下的鼓风密度,P0表示1标准大气压,L表示从回旋区最深处到回旋区边界上任意一点的曲线距离,σ表示鼓风在回旋区边界上的损失率,H表示高炉的总高度,h表示B点到风口的垂直距离,a、b为回旋区边界的形状参数,I表示边界点B处的截面包括的焦炭颗粒总数,为回旋区边界点B处第i个焦炭颗粒的直径;
因为回旋区最深处到回旋区上的任意一点的边界上是曲线形状,所以设定L=(DR-x)a+hb,其中,回旋区边界点B相对于风口的坐标为(x,h);
由上述回旋区边界点处的数学模型进一步得到回旋区边界的计算公式为:
K1·x2-2K1·DR·x+K1·h2-K2·h=K3
其中,
其中,K3和K4为待定系数;
步骤3、获得建模参数,分析建模参数对回旋区边界模型的影响,确定影响回旋区边界的主要参数;
根据回旋区边界模型以及建模涉及到的相关参数,由回旋区边界变化规律可知鼓风参数中的鼓风风压P和鼓风风量Vg以及回旋区温度Tm以及焦炭的鼓风损失率σ对于回旋区边界的变化有着直接的影响作用;由回旋区边界的计算公式可知,降低鼓风风压P、提高鼓风风量Vg、增大回旋区温度Tm以及降低回旋区焦炭处的鼓风损失率都有利于回旋区向中心发展,增加回旋区的深度;其中,鼓风风压P、鼓风风量Vg可调节,回旋区温度Tm的变化受多种因素的影响,无法直接调节使其变化;
步骤4、利用回旋区边界计算公式求出回旋区的高度;
对回旋区边界计算公式的变量h进行求导,并将求导的结果等于零就求出回旋区在竖直方向上的最大值,即回旋区的高度,如下公式所示:
其中,G是回旋区的高度,K5为待定系数;
步骤5、当回旋区高度或回旋区深度超出设定范围时,通过调节鼓风风压P和鼓风风量Vg使回旋区高度或深度恢复至正常范围内。
采用上述技术方案所产生的有益效果在于:本发明提供的一种高炉风口回旋区边界各点的计算及实时监测方法,首先确定回旋区深度模型,并在该模型的基础上,利用物理原理以及化学原理推导出回旋区边界上任意一点的数学模型,通过这个数学模型就可以直观的反应高炉风口回旋区的口袋形状,并成功引入回旋区温度这个变量,获取影响回旋区深度和边界的两个鼓风参数为鼓风风压和鼓风温度,研究这两个参数对回旋区深度和回旋区边界变化的影响,以及当回旋区深度和边界发生不合理变化时,通过及时调节这两个参数,使回旋区深度和边界恢复至正常范围内;本发明方法能够实时监测回旋区深度和回旋区边界的变化情况,为高炉的实际生产提供安全指导,以保障企业的生命财产安全。
附图说明
图1为本发明实施例提供的高炉风口回旋区形成示意图;
图2为本发明实施例提供的回旋区的深度计算模型建模示意图,其中,(a)为回旋区深度机理模型建模整体示意图,(b)为风口的正向截面图;
图3为本发明实施例提供的回旋区边界模型建模示意图;
图4为本发明实施例提供的回旋区深度变化规律示意图;
图5为本发明实施例提供的回旋区高度变化规律示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
高炉风口回旋区形成如图1所示,高炉炼铁过程中,鼓风离开风口时具有很强的动能,它吹动风口前的焦炭并与之发生燃烧反应,在风口前缘形成一疏松且近似呈椭圆形的气相空穴。另外风口前的煤气流以回旋区为放射中心,分别沿长径向炉缸中心发展,沿短径向两侧发展,与此同时自空穴上部和两侧不断有新的焦炭补充近来使得焦碳在空腔内作回旋运动,这个区域为高炉风口回旋区。本实施例以4000m3的高炉为例,采用本发明的高炉风口回旋区边界的计算及实时监测方法实现对该高炉回旋区边界的计算及实时监测。
本发明实施例中,在回旋区深度计算模型的基础上研究对象为回旋区边界处上任意一点B,通过对B处进行受力分析并与最深点进行比较,建立回旋区边界模型,并基于回旋区边界模型进行回旋区边界的监测。
一种高炉风口回旋区边界各点的计算及实时监测方法,包括以下步骤:
步骤1、根据高炉风口回旋区的形成原理,建立回旋区的深度计算模型,如图2所示,进而得到回旋区深度的计算公式,获得回旋区深度的变化规律;
当回旋区内部运动处于稳定状态时,取回旋区内部最深处一区域A为研究对象,此时A处在鼓风气体冲力与焦炭层阻力共同作用下达到平衡,依据二力平衡建立回旋区的深度计算模型,用以求解回旋区深度的变化情况,如下公式所示:
其中,FA表示区域A处鼓风气体冲力,FB表示区域A处焦炭层阻力,ρg0表示标准态下的鼓风密度,Vg表示区域A处鼓风风量,ST表示风口面积,DR表示回旋区深度,DT表示风口直径,α为常数,用于表示回旋区深度与宽度之间的关系,P表示鼓风风压,Tm表示由CCD测温仪所测回旋区温度,SP表示区域A处焦炭粒子总截面积,ρP表示焦炭密度,VP表示区域A处焦炭粒子体积,g表示重力加速度,DPR表示回旋区边界前焦炭直径,DPR=0.6DPc,DPc表示入炉前焦炭直径,ρg表示实际风温风压下的鼓风密度,Tw表示鼓风温度,P表示鼓风风压;
由上述回旋区的深度计算模型得到回旋区深度的计算公式为:
其中,K、β均为待定系数;
本实施例在经过大量试验验证后,最终求得回旋区深度的计算公式为:
步骤2、通过高炉风口回旋区的深度计算模型建立高炉风口回旋区的边界模型,确定回旋区边界的计算公式;
当高炉回旋区的风口鼓风从风口吹到最深处时,鼓风会向上、向下继续行进知道直到风口回旋区的最高处,所以可以把从高炉风口回旋区最深处作为向上鼓风的风源继续推导模型,鼓风从风口回旋区最深处开始向上吹出风口回旋区的空腔直到回旋区的最高处,在高炉风口回旋区深度计算模型的建立公式的基础上,由于鼓风在回旋区边界的损失再利用力平衡原理可以推导出边高炉风口回旋区界各点的边界模型。
高温鼓风从风口吹入后,其行程近似为一个圆锥形管道到达回旋区最深处A点,之后鼓风转而向上运动,在距离回旋区风口水平距离x米、距离回旋区风口垂直距离h米的位置形成一块微元区域B,其截面积为SB,在B处鼓风气体冲力FD与焦炭层对其的阻力Fb共同作用下达到平衡状态,此时认为回旋区处于稳定运行状态,满足表达式
FD=Fb (2)其中,鼓风气体冲力FD的大小主要与鼓风质量和鼓风速度成正比,即:
式中,MB表示处鼓风质量流率,具体含义为单位时间内通过该截面的流体质量,即:
MB∝ρg·Vg·(1-σL) (4)
其中,ρg表示鼓风密度,Vg表示鼓风风量,L表示从回旋区最深处到回旋区边界上任意一点的曲线距离,σ表示鼓风在回旋区边界上的损失率。
对L进行进一步的讨论,因为回旋区最深处到回旋区上的任意一点的边界上是曲线形状,所以假设L=(DR-x)a+hb,其中,(x,h)为B点相对于风口的坐标,a、b为曲线参数。
UB代表B处风速,B处的风速与鼓风风压成正比,与回旋区内部的实际温度的也成正比,即:
式中,MB表示B处鼓风质量流率,ρg表示鼓风密度,SB为B处截面积,P为鼓风风压,P0为1标准大气压,Tm为测温仪实测回旋区温度,T0=298K。
SB为B处截面积,其表达式为:
式中,WR为回旋区宽度;
SP为B处焦炭粒子总截面积,其表达式为:
式中,DPR为回旋区边界前焦炭直径,DPc表示入炉前焦炭直径,DPR=0.6DPc。
联立(3)、(4)式可以进一步推出B处鼓风气体冲力FD的表达式为:
实验证明回旋区深度DR和回旋区宽度WR之间的关系为:
式中,α为一常数,α=0.2~0.6;
设定ST表示风口面积,其表达式为:
联立(5)、(8)、(9)式可以推出:
最后,将表达式(10)代入表达式(7)可以得到B处鼓风气体冲力FD的最终表达式为:
Fb表示B处焦炭层对鼓风气流的阻力,它主要由焦炭的压实程度决定,近似认为它与焦炭颗粒的重力和鼓风风压成正比,即:
式中,ρP表示焦炭密度,VP表示A处焦炭粒子体积,g表示重力加速度,DPR表示回旋区边界前焦炭直径,P表示鼓风风压,H表示高炉的总高度,h表示B点到风口的垂直距离;鼓风密度ρg又可以进一步表示为:
式中,μg0表示标准态下的鼓风密度,Tw表示鼓风温度,P表示鼓风风压,P0表示1标准大气压。
最后将表达式(11)、(12)、(13)代入表达式(1),可以求得回旋区边界的计算公式为:
K1·x2-2K1·DR·x+K1·h2-K2·h=K3 (15)
其中,
其中,K3、K4均为待定系数;
利用上述回旋区边界的计算公式,可以进一步建立回旋区的数学模型,研究回旋区边界的变化趋势。
步骤3、获得建模参数,分析建模参数对回旋区边界模型的影响,确定影响回旋区边界的主要参数;
根据如图3所示的回旋区边界模型以及建模涉及到的相关参数,由回旋区边界变化规律可知鼓风参数中的鼓风风压P和鼓风风量Vg以及回旋区温度Tm对于回旋区边界的变化有着直接的影响作用;由回旋区边界的计算公式可知,降低鼓风风压P、提高鼓风风量Vg、增大回旋区温度Tm有利于回旋区向中心发展,增加回旋区的深度,使得回旋区边界向外扩张,这满足实际高炉冶炼过程中回旋区边界的变化规律,同时初步验证了上述风口回旋区边界计算公式的合理性;其中,鼓风风压P、鼓风风量Vg可调节,回旋区温度Tm的变化受多种因素的影响,无法直接调节使其变化;
本实施例中,回旋区边界模型建模过程中用到的相关参数如表1所示:
表1回旋区边界模型建模过程中用到的相关参数
步骤4、利用回旋区边界计算公式求出回旋区的高度;
对回旋区边界计算公式的变量h进行求导,并将求导的结果等于零就求出回旋区在竖直方向上的最大值,即回旋区的高度,如下公式所示:
其中,G是回旋区的高度,K5是待定系数;
本实施例在经过大量试验验证后,最终求得K5=2.12,且回旋区高度的计算公式为:
步骤5、当回旋区高度或回旋区深度超出设定范围时,通过调节鼓风风压P和鼓风风量Vg使回旋区高度或深度恢复至正常范围内。
本实施例中,高炉平稳运行下回旋区深度变化和高度变化如图4所示,其中建模所用参数和数据全部来自于4000m3的高炉,在实际的高炉冶炼中4000m3的高炉正常情况下的回旋区深度在1.6m至1.8m之间,回旋区高度在1.1m到1.7m之间,由图4和图5可以看出,此时的高炉生产处于正常状态,回旋区深度在1.65m至1.75m内波动,回旋区高度在1.1m至1.68m内波动,均处在正常范围内。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (5)
1.一种高炉风口回旋区边界的计算及实时监测方法,其特征在于:包括以下步骤:
步骤1、根据高炉风口回旋区的形成原理,建立回旋区的深度计算模型,进而得到回旋区深度的计算公式,获得回旋区深度的变化规律;
步骤2、通过高炉风口回旋区的深度计算模型建立高炉风口回旋区的边界模型,确定回旋区边界的计算公式;
步骤3、获得建模参数,分析建模参数对回旋区边界模型的影响,确定影响回旋区边界的主要参数;
步骤4、利用回旋区边界计算公式求出回旋区的高度;
步骤5、当回旋区高度或回旋区深度超出设定范围时,通过调节鼓风风压和鼓风风量使回旋区高度或深度恢复至正常范围内。
2.根据权利要求1所述的一种高炉风口回旋区边界的计算及实时监测方法,其特征在于:所述步骤1的具体方法为:
当回旋区空穴处于稳定运动状态时,取回旋区内部最深处一微元区域A为研究对象,此时A处在鼓风气体冲力与焦炭层阻力共同作用下达到平衡,依据二力平衡建立回旋区的深度计算模型,用以求解回旋区深度的变化情况,如下公式所示:
其中,FA表示区域A处鼓风气体冲力,FB表示区域A处焦炭层阻力,ρg0表示标准态下的鼓风密度,Vg表示区域A处鼓风风量,ST表示风口面积,DR表示回旋区深度,DT表示风口直径,α为常数,用于表示回旋区深度与宽度之间的关系,P表示鼓风风压,Tm表示回旋区温度,SP表示区域A处焦炭粒子总截面积,ρP表示焦炭密度,VP表示区域A处焦炭粒子体积,g表示重力加速度,DPR表示回旋区边界前焦炭直径,DPR=0.6DPc,DPc表示入炉前焦炭直径,Tw表示鼓风温度,P表示鼓风风压;
由上述回旋区深度的计算模型得到回旋区深度的计算公式为:
其中,K、β均为待定系数,ρg表示实际风温风压下的鼓风密度。
3.根据权利要求2所述的一种高炉风口回旋区边界的计算及实时监测方法,其特征在于:所述步骤2的具体方法为:
当回旋区内部运动处于稳定状态时,任意取回旋区边界点B为研究对象,此时边界点B处在鼓风气体冲力与焦炭层阻力共同作用下达到平衡,依据二力平衡建立回旋区边界处任意一点的数学模型,用以求解回旋区边界的变化情况:
其中,
MB∝ρg·Vg·(1-σL)
SB∝(1-σ((DR-x)a+hb))·SA
其中,FD表示边界点B处鼓风气体冲力,Fb表示边界点B处焦炭层对鼓风气流的阻力,MB表示B处鼓风质量流率,UB代表B处风速,SB为B处的截面积,SA表示区域A处的截面积,ρg表示实际风温风压下的鼓风密度,P0表示1标准大气压,L表示从回旋区最深处到回旋区边界上任意一点的曲线距离,σ表示鼓风在回旋区边界上的损失率,H表示高炉的总高度,h表示B点到风口的垂直距离,a、b为回旋区边界的形状参数,I表示边界点B处的截面包括的焦炭颗粒总数,为回旋区边界点B处第i个焦炭颗粒的直径;
因为回旋区最深处到回旋区上的任意一点的边界上是曲线形状,所以设定L=(DR-x)a+hb,其中,回旋区边界点B相对于风口的坐标为(x,h);
由上述回旋区边界点处的数学模型进一步得到回旋区边界的计算公式为:
K1·x2-2K1·DR·x+K1·h2-K2·h=K3
其中,
其中,K3和K4为待定系数。
4.根据权利要求3所述的一种高炉风口回旋区边界的计算及实时监测方法,其特征在于:所述步骤3的具体方法为:
根据回旋区边界模型以及建模涉及到的相关参数,由回旋区边界变化规律可知鼓风参数中的鼓风风压P和鼓风风量Vg以及回旋区温度Tm以及焦炭的鼓风损失率σ对于回旋区边界的变化有着直接的影响作用;由回旋区边界的计算公式可知,降低鼓风风压P、提高鼓风风量Vg、增大回旋区温度Tm以及降低回旋区焦炭处的鼓风损失率都有利于回旋区向中心发展,增加回旋区的深度;其中,鼓风风压P、鼓风风量Vg可调节,回旋区温度Tm的变化受多种因素的影响,无法直接调节使其变化。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110566801.4A CN113283079B (zh) | 2021-05-24 | 2021-05-24 | 一种高炉风口回旋区边界的计算及实时监测方法 |
US18/555,827 US20240211662A1 (en) | 2021-05-24 | 2021-12-03 | Calculating and real-time monitoring method for boundary of blast furnace tuyere raceway |
PCT/CN2021/135430 WO2022247209A1 (zh) | 2021-05-24 | 2021-12-03 | 一种高炉风口回旋区边界的计算及实时监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110566801.4A CN113283079B (zh) | 2021-05-24 | 2021-05-24 | 一种高炉风口回旋区边界的计算及实时监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113283079A true CN113283079A (zh) | 2021-08-20 |
CN113283079B CN113283079B (zh) | 2024-03-19 |
Family
ID=77281131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110566801.4A Active CN113283079B (zh) | 2021-05-24 | 2021-05-24 | 一种高炉风口回旋区边界的计算及实时监测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240211662A1 (zh) |
CN (1) | CN113283079B (zh) |
WO (1) | WO2022247209A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113283078A (zh) * | 2021-05-24 | 2021-08-20 | 东北大学 | 一种高炉风口回旋区深度的计算及实时监测方法 |
CN114606354A (zh) * | 2022-03-31 | 2022-06-10 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区高度的方法 |
CN114752717A (zh) * | 2022-03-31 | 2022-07-15 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区宽度的方法 |
CN114908202A (zh) * | 2022-04-27 | 2022-08-16 | 鞍钢股份有限公司 | 利用炉缸取样手段来计算整体炉缸活跃区域大小的方法 |
WO2022247209A1 (zh) * | 2021-05-24 | 2022-12-01 | 东北大学 | 一种高炉风口回旋区边界的计算及实时监测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102758039A (zh) * | 2011-04-28 | 2012-10-31 | 宝山钢铁股份有限公司 | 一种定量评价风口回旋区活跃性的方法 |
KR20140016635A (ko) * | 2012-07-30 | 2014-02-10 | 현대제철 주식회사 | 고로의 장입물 강하 판단 방법 |
CN106815405A (zh) * | 2016-12-20 | 2017-06-09 | 武汉钢铁股份有限公司 | 一种高炉风口回旋区深度计算方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113283079B (zh) * | 2021-05-24 | 2024-03-19 | 东北大学 | 一种高炉风口回旋区边界的计算及实时监测方法 |
-
2021
- 2021-05-24 CN CN202110566801.4A patent/CN113283079B/zh active Active
- 2021-12-03 WO PCT/CN2021/135430 patent/WO2022247209A1/zh active Application Filing
- 2021-12-03 US US18/555,827 patent/US20240211662A1/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102758039A (zh) * | 2011-04-28 | 2012-10-31 | 宝山钢铁股份有限公司 | 一种定量评价风口回旋区活跃性的方法 |
KR20140016635A (ko) * | 2012-07-30 | 2014-02-10 | 현대제철 주식회사 | 고로의 장입물 강하 판단 방법 |
CN106815405A (zh) * | 2016-12-20 | 2017-06-09 | 武汉钢铁股份有限公司 | 一种高炉风口回旋区深度计算方法及*** |
Non-Patent Citations (2)
Title |
---|
侴宏晔 等: "竖炉风口回旋区简化模型的研究", 上海金属, vol. 26, no. 02, pages 21 - 24 * |
侴宏晔 等: "高炉风口回旋区简化模型的探讨", 上海大学学报(自然科学版), vol. 10, no. 01, pages 9 - 12 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113283078A (zh) * | 2021-05-24 | 2021-08-20 | 东北大学 | 一种高炉风口回旋区深度的计算及实时监测方法 |
WO2022247209A1 (zh) * | 2021-05-24 | 2022-12-01 | 东北大学 | 一种高炉风口回旋区边界的计算及实时监测方法 |
CN113283078B (zh) * | 2021-05-24 | 2024-04-26 | 东北大学 | 一种高炉风口回旋区深度的计算及实时监测方法 |
CN114606354A (zh) * | 2022-03-31 | 2022-06-10 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区高度的方法 |
CN114752717A (zh) * | 2022-03-31 | 2022-07-15 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区宽度的方法 |
CN114606354B (zh) * | 2022-03-31 | 2023-03-17 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区高度的方法 |
CN114752717B (zh) * | 2022-03-31 | 2023-03-17 | 鞍钢股份有限公司 | 一种利用炉缸取样手段分析风口回旋区宽度的方法 |
CN114908202A (zh) * | 2022-04-27 | 2022-08-16 | 鞍钢股份有限公司 | 利用炉缸取样手段来计算整体炉缸活跃区域大小的方法 |
CN114908202B (zh) * | 2022-04-27 | 2023-03-17 | 鞍钢股份有限公司 | 利用炉缸取样手段来计算整体炉缸活跃区域大小的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113283079B (zh) | 2024-03-19 |
WO2022247209A1 (zh) | 2022-12-01 |
US20240211662A1 (en) | 2024-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113283079A (zh) | 一种高炉风口回旋区边界的计算及实时监测方法 | |
CN104212924B (zh) | 一种高炉气流分布的检测方法 | |
CN104531924B (zh) | 一种高炉炉料分布实时预报***及方法 | |
CN111174569A (zh) | 一种在线预测回转窑内煅烧段烟气温度的方法及*** | |
CN111575421B (zh) | 利用扫描雷达调整软融带形状的方法 | |
Zhang et al. | Gas-particle flow and combustion characteristics of pulverized coal injection in blast furnace raceway | |
Dong et al. | Gas-solid flow and shaft injected gas penetration in an oxygen blast furnace analyzed using a three-dimensional DEM-CFD coupling mathematical model | |
CN107515999B (zh) | 一种用于确定高炉炉顶煤气流运动与粉尘分布的方法 | |
CN105045994B (zh) | 一种优化带式焙烧机热工制度的方法 | |
CN105950806B (zh) | 一种用于确定高炉炉顶内的低粉尘区域的方法以及装置 | |
CN113283078B (zh) | 一种高炉风口回旋区深度的计算及实时监测方法 | |
CN116227387A (zh) | 一种煤气化炉管道冲蚀预测及调控方法 | |
CN110427693B (zh) | 带有喷枪的工业反应器内流体仿真计算方法 | |
Zhao et al. | Three-dimensional numerical simulation of flow and splash behavior in an oxygen coal combustion melting and separating furnace | |
Zhang et al. | Numerical investigation of gas-solid heat transfer process and parameter optimization in shaft kiln for high-purity magnesia | |
CN116497166B (zh) | 一种高炉炉顶布料工艺流程的智能控制方法 | |
Tang et al. | Mathematical modeling of the argon oxygen decarburization converter exhaust gas system at the reduction stage | |
JP2003129146A (ja) | 精鉱バーナの制御方法 | |
Lyalyuk et al. | Determining the gas trajectory in blast-furnace injection of pulverized coal | |
CN117690343A (zh) | 确定炉顶煤气循环-氧气高炉煤气流分布状态的冷态实验装置及方法 | |
CN111599415A (zh) | 一种基于计算机仿真的高炉数字化***实现方法 | |
Gao et al. | CFD-DEM simulation of raceway size and mechanical characteristics of industrial scale blast furnace | |
CN116497166A (zh) | 一种高炉炉顶布料工艺流程的智能控制方法 | |
Odenthal et al. | Recent SMS siemag developments in BOF steelmaking | |
WO2023070692A1 (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 |