《高炉数学模型的进展_储满生.pdf》由会员分享，可在线阅读，更多相关《高炉数学模型的进展_储满生.pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第 17 卷第 4 期2007 年 4 月 中国冶金 China Metallurgy Vol.17,No.4 Apr.2007作者简介:储满生(1973-),男,博士后;E-mail:chumansheng ;修订日期:2006-11-22高炉数学模型的进展储满生1,郭宪臻2,沈峰满1,八木顺一郎3(1.东北大学材料与冶金学院,辽宁 沈阳 110004;2.安阳钢铁集团公司炼铁厂,河南 安阳 455004;3.東北大学多元物質科學研究所,日本 仙台 980-8577)摘 要:高炉是一个复杂的气固向流反应器。为了理解、控制和改进高炉炼铁过程,更多的努力被用于开发数学模型。全高炉动力学模型经历了
2、从一维到多维、稳态到非稳态、简单到复杂的发展过程。多流体高炉模型是全高炉动力学数学模型研究的最新成果之一。该模型基于多流体理论,可对炉内主要现象进行多维数学模拟,并能精确地预测高炉在指定条件下的操作指标。未来的高炉数学模型应朝复杂化和实用化方向发展,并拓展解析范围、提高模拟精度。关键词:高炉;数学模型;多流体理论中图分类号:TF531 文献标识码:A 文章编号:1006-9356(2007)04-0010-05Development of Blast Furnace Mathematical ModelCHU Man-sheng1,GUO Xian-zhen2,SHEN Feng-man1,Y
3、AGI Jun-ichiro3(1.School of Materials and M etallurgy,Northeastern University,Shenyang 110004,Liaoning,China;2.Iron-Making Plant,Anyang Iron and Steel Co Ltd,Anyang 455004,Henan,China;3.Institute ofMultidisciplinary Research for Advanced M aterials,Tohoku University,Sendai,980-8577,Japan)Abstract:Bl
4、ast furnace is a complex metallurgical reactor with gas-solid counter flow.In order to understand,con-trol and improve the blast furnace process,more efforts are made to develop mathematical models.In the past dec-ades,the reaction-kinetics models of total blast furnace have generally developed from
5、 one-dimensional,stead stateand simple models to mult-i dimensional,transient and complex ones.Mult-i fluid blast furnace model is one of thelatest reaction-kinetics models,which is based on mult-i fluid theory and can make mult-i dimensional treatment ofin-furnace phenomena and accurate prediction
6、of operational parameters under the given conditions.T he futuremodels of blast furnace should be more complicated and applicable after enlarging simulation range and improvingcomputation exactness.Key words:blast furnace;mathematical model;mult-i fluid theory 高炉生产是由炉顶加入炉料,从炉缸渣铁口排放渣铁;而从风口鼓入热风和喷吹煤粉,煤
7、气从炉顶逸出。所以,可将高炉看作是存在炉料下降和煤气上升 2 个逆向运动的反应器,炉内所有传输现象和反应均发生于炉料与煤气的向流运动中。高炉内多种多相物质共存且相互作用,且诸多物理化学现象同时发生,故高炉在化工领域被认为是最复杂的冶金反应器之一。随现代测控技术的发展和对实际运行高炉的炉体解剖,人们获得了大量有用数据,对高炉有了更深理解,但光凭这些仍很难详细掌握炉内的现象。因此,为更好地理解、控制和改进高炉炼铁过程,更多的努力用于开发高炉数学模型 1。根据对高炉内传输现象描述方法的不同,数学模型分为平衡理论模型、热化学模型和反应动力学模型;按所考虑的空间坐标维数,分为一维、二维和三维模型;考虑时
8、间变量与否,有稳态模型和非稳态模型 1 2。高炉炼铁技术在过去几十年中得到长足的进步。为提高产量、降低能耗和减少环境负荷,喷吹煤粉和有效利用含碳含铁废弃物等大量新技术逐渐被采用,高炉的功能不断扩大,炉内现象更趋复杂化。而建立在反应动力学和传输现象理论基础上的反应动力学数学模型是一有用的工具,可用其详细分析炉内状态和精确预测高炉操作性能。1 高炉数学模型的发展历程 高炉反应动力学模型的发展历程见图 1。最早获得发展的是高炉一维模型,且先有稳态模型,随后逐渐发展为非稳态。鞭严等人在 20 世纪 60 年代末开发的高炉稳态一维模型最具代表性 2。在鞭严的模型中,考虑了炉内主要化学反应和传热过程,模拟
9、结果给出了主要工艺变量沿高炉高度方向上的分布。后来,许多研究者仿效鞭严的建模思想建立了一系列用于解决不同问题的高炉数学模型。这些早期的高炉模型很好地把握了对局部和全高炉的能量平衡和物质平衡这一基本规律,因而在模拟高炉现象、分析操作参数对炉况和冶炼指标的影响、指导开停炉等方面获得了相当的成功。如一些模型应用于分析鼓风压力波动对高炉操作的影响、预测最低燃料比以及模拟顶压操作等实践 2。但是,对于一维高炉模型来说,过程参数被假设为径向均布,而高炉的解剖和取样分析证实了气体温度和炉内物质成分等在径向上都是不均匀的。另外,在其建模过程中,炉内物质和能量的传输过程只能通过常微分方程来描述。再加上边界值设定
10、不合理等缺陷,这些早期一维模型的预测精度和应用范围都很有限。到 80 年代,计算机技术的发展允许模型处理更大的矩阵,新建立的模型可采用偏微分方程作为其控制方程。在此期间,大量二维高炉模型被开发,其中较为知名的有 Hatano 和 Kurita 的模型,Yagi、Takeda 和 Omori 的模型,Sugiyama 和 Sugata 的BRIGHT 模型 1,2。二维模型主要被用于描述炉内更为复杂的现象,即评估操作条件对高炉操作性能和炉况的影响,分析软熔带的变化和影响,模拟和开发高炉炼铁新技术等。总体来说,这些模型对指导实际高炉操作和促进炼铁技术的进步做出了一定的贡献。到 90 年代初,一个基
11、本概念/多流体理论0被提出,即用多相流和相间双向相互作用来描述发生在炉下部的现象,且炉内物质相应通过流动机制来加以区分。因此,除了最基本的气-固-液态物质外,被炉内气流挟带的未燃煤粉被处理为一独立的粉相 3。在随后基于这个理论而发展的高炉数学模型中,根据物性的不同,液相又进一步划分为渣相和铁水相,而粉相分为静态滞留粉相和动态滞留粉相。这些模型都总称为/多流体高炉数学模型0 4 10,属于较为复杂全面的全高炉反应动力学模型。目前,多流体模型能合理处理二维和三维的问题。2 多流体高炉数学模型 多流体高炉数学模型基于多流体理论、冶金传输原理、反应动力学以及计算流体力学等理论模拟高炉炼铁过程 4,5,
12、7。模型主要考虑了从渣面到料面的整个填充区域,且认为高炉是轴对称的。多流体模型中将气相、固相(焦炭和含铁炉料)、铁水、熔渣和粉相分别视为具有各自流动机制的独立相,每相由 1 个或多个组元组成,且每个组元具有独立的成分和物性。模型中考虑的物质和相见表 1。图 1 全高炉反应动力学模型的发展历程Fig.1 Development of blast furnace mathematical models表 1 多流体模型中考虑的相和物质Table 1 Phases and chemical species considered in mult-i fluid model相物质气相CO,CO2,O2,
13、H2,H2O,N2,SiO,CH4固相块矿Fe2O3,Fe3O4,FeO,Fe,脉石(CaO,Al2O3,MgO),SiO2,H2O烧结矿Fe2O3,Fe3O4,FeO,Fe,脉石(CaO,Al2O3,MgO),SiO2,H2O球 团Fe2O3,Fe3O4,FeO,Fe,脉石(CaO,Al2O3,MgO),SiO2,H2O焦炭C,SiC,脉石(CaO,Al2O3,MgO),SiO2,H2O铁水Fe,C,Si熔渣FeO,脉石(CaO,Al2O3,M gO),SiO2粉相煤粉(废塑料)C,脉石(CaO,Al2O3,MgO),SiO2,挥发分11第 4 期 储满生等:高炉数学模型的进展 多流体模型中
14、各相间的相互作用见图 2。图中实线代表动量、质量和能量的完全相互作用,点线代表质量的传输。模型认为气-固相与其它相之间均有完全的质量、动量和能量交换,而不连续相(液相和粉相)间不进行动量交换,但通过化学反应和相变进行质量和能量的交换。图 2 多流体模型中相间相互作用示意Fig.2 Diagram of interphase interactions consideredin mult-i fluid model 多流体模型采用一系列偏微分方程来描述炉内各相的行为。由于各相之间同时且双向相互作用,故这些偏微分方程是强烈耦合的,须同时求解。所有方程可用一个统一化的形式来加以描述,即:99x(EiQ
15、iuiW)+1r99r(rEiQiviW)=99xEiSW9W9x+1r99rrEiSW9W9r+SW(1)式中 i)要考虑的气相、固相、铁水、熔渣和粉相;W)求解的变量,通过改变 W,分别代表质量、动量、能量方程以及连续性方程;Ei)相 i 的体积分数;SW)有效扩散系数,根据所求解的独立变量具有不同的意义。如对于动量方程,SW表示动力粘度;而对于能量方程,SW为导热系数;SW)源相,主要是由于化学反应、相间相互作用、外力以及相变等因素而产生。由于化学反应和相变,质量方程则具有质量源。模型中主要考虑的反应是:铁氧化物间接还原;熔渣中氧化铁(FeO)直接还原;回旋区燃烧反应;焦炭气化溶损反应;
16、水煤气反应和水煤气转变反应;硅迁移反应;渗碳反应;熔化及相变反应;水分蒸发或凝结反应等。并利用已有的研究成果对上述反应的速度进行了动力学公式化描述 4 7。不同相间因速度差的存在而产生动量交换,相间的动量交换量计算为相间速度差和相间动量传输系数的乘积。模型中针对不同体系采用不同的动量传输系数计算方法 3 7,包括:气-固:Ergun 方程;气-液:单液滴拽力方程;气-粉:经 Richardson和 Zaki 修正的单颗粒拽力方程;固-液:Kozeny-Carman 方程;固-粉:Fanning 方程;液-粉:未考虑。能量方程的源相包括反应热和相间对流换热。相间对流换热量通过相间对流换热系数、相
17、间接触面积以及相间温度差的乘积来定量估算。针对不同体系而采用的相间对流换热系数估算方法 3 7有:气-固:经 Akiyama 修正的 Ranz 方程;气-液:Mackey-Warner 关系式;气-粉:Ranz-Marshall 方程;固-液:Eckert-Drake 方程;固-粉:Emulsion模型;液-粉:未考虑。模型计算时同时求解所有方程。首先,在计算区域内利用 BFC 法 11进行网格化;随后,利用控制单元体法 12离散化所有方程;最后,采用 SIMPLE法 12和迭代矩阵法求解所有离散的方程。根据计算结果可给出:炉内温度、压力、速度、反应速度、体积分量以及相物质组成等变量的多维分布
18、;产量、利用系数、渣量、焦比、煤比、还原剂消耗总量;以及炉顶煤气利用率等操作指标的变化。利用某实际运行高炉的生产数据,对多流体数学模型的有效性和计算的精确度进行了充分验证。图 3 给出了模型预测的炉内二维温度分布与实测值(受炉内测温元件的限制,1200 e 以上温度区域未能给出)。模型基本预测了炉内温度分布走向,模拟值与测量值显示了良好的一致性。表 2 给出了模型预测值和实际高炉操作参数之间的对比情况,主要项的预测误差低于 5%。应该说多流体模型较准确预测了高炉操作指标。3 展望 高炉数学模型未来的发展方向见图 4。高炉数学模型的核心是提高数学模型的实用性和精确度。首先是扩展和丰富目前的建模体
19、系、思想和方法。为此,一些尚未解决的现象应在新的模型中加以考虑,包括:粉相和液相静态滞留和动态滞留的区分,以及 2 种滞留物质之间传输机制和速率;炉内固体物料粉化而导致的粉相产生机制及速率;脉12 中国冶金 第 17 卷(a)实测值;(b)计算值图 3 实测和模型计算的炉内温度分布(e)Fig.3 Calculated and measured temperaturedistribution in furnace表 2 高炉模型预测值和实际操作参数的对比Table 2 Comparison between predicted values andpractical operation para
20、meters 项目模型计算值实际测量值误差/%生铁产量/(kg#s-1)103.5 104.3 0.7焦比/(kg#t-1)313.6 313.0 0.2渣量/(kg#t-1)252.8 261.0 3.1煤气利用效率/%53.3 50.6 5.3炉顶煤气温度/e248.0 240.0 3.3石混合物在炉内的熔融现象和速率;铁水中微量元素的扩散机制和基于速度论的描述;液相和粉相之间的相互双向作用机制;生铁渗碳的机理和速率;/死料柱0更新的机理和速率。第 2 个应获得发展的新方向是运用粒子追踪法和随机过程理论来模拟炉内某些现象。这主要是由于实际操作中出现的如管道、悬料、滑料及风压波动等炉况失常现
21、象,都认定是不连续性和随机现象所致。最后,应针对特定对象和现场应用目标来简化图 4 高炉数学模型的未来发展方向Fig.4 Perspectives on development of blast furnace models模型。除改进数学模型本身外,更为重要的是改善人-机操作界面,促进和扩大数学模型在高炉实际生产中的应用。4 结论 (1)高炉动力学数学模型建立在反应动力学和传输现象理论基础之上,是详细分析炉内状态和精确预测高炉操作性能的有效工具。其发展经历了由一维到多维,稳态到非稳态,考虑的现象由简单到复杂的过程。(2)基于多流体理论的多流体高炉模型是一复杂全面的高炉过程模拟仿真系统,可详尽
22、模拟高炉内的主要现象,并能精确预测高炉操作指标,模型预测值与实测值显示了良好的一致性。(3)高炉数学模型的未来发展应扩展和丰富目前的建模体系、思想和方法,采用新理论,针对特定13第 4 期 储满生等:高炉数学模型的进展 对象和实际应用目标而简化模型,从而提高数学模型的实用性和精确度。参考文献:1 毕学工.高炉过程数学模型及计算机控制 M.北京:冶金工业出版社,1996.2 Omori K.Blast Furnace Phenomena and Modelling M .Lon-don:Elsevier Applied Science,1987.3 Yagi J.Mathematical M o
23、deling of the Flow of Four Fluids in aPacked Bed J.ISIJ International,1993,33(6):619-639.4 Austin P R,Nogami H,Yagi J.A Mathematical M odel forBlast Furnace Reaction Analysis Based on the Four Fluid Mod-el J.ISIJ International,1997,37(8):748-755.5 Castro J A,Nogami H,Yagi J.T ransient M athematical
24、Mod-el of Blast Furnace Base on Multi-fluid Concept,with Appl-ication to High PCI Operation J.ISIJ International,2000,40(7):637-646.6 Sungging P,Nogami H,Yagi J.Numerical Analy sis of StaticHoldup of Fine Particles in Blast Furnace J.ISIJ Internation-al,2004,42(2):304-309.7 Chu M.Study on Super High
25、 Efficiency Operations of BlastFurnace Based on Mult-i fluid M odel D.Sendai:T ohoku Un-iversity,Japan,2004.8 Chu M,Nogami H,Yagi J.Numerical Analysis on Injectionof Hydrogen Bearing Materials into Blast Furnace J.ISIJ In-ternational,2004,44(3):801-808.9 Chu M,Nogami H,Yagi J.Numerical Analysis on C
26、hargingCarbon Composite Agglomerates into Blast Furnace J.ISIJInternational,2004,44(5):510-517.10 Chu M,Nogami H,Yagi J.Numerical Analysis on BlastFurnace Performance Under Operation with T op Gas Recy-cling and Carbon Composite Agglomerates Charging J.ISIJInternational,2004,44(12):2159-2167.11 Thom
27、pson J F,Warsi Z U,Mastin C W.Numerical GridGeneration M.New York:North Holland,1985.12 Patankar S.Numerical Heat T ransfer and Fluid Flow M .New York:M cGRAW-Hill,1980.第八届全国连铸年会在海口召开 由中国金属学会连续铸钢分会举办的第八届全国连铸年会于 2007 年 3 月 22 27 日在海口召开。这次年会的主题是:求索、总结、创新、发展,/十一五0再建连铸辉煌。会议讨论的内容包括:连铸新技术新工艺的发展和应用;炉外精炼-连铸
28、-连轧,协调优化生产高质量、高附加值产品;提高铸机拉坯速度、调整结构,优化生产;节能、环保的有关问题;新连铸设备新技术发展动态;现有连铸机进一步改造问题;连铸理论研究在实际生产中的应用;有关钢水精练、中间包冶金、连铸二级软件包的研发和应用、结晶器专家系统、轻压下技术、动态自动配水;合金钢大方坯、圆坯连铸以及异型坯连铸的生产技术及设备特点等与连铸相关的技术。大会由连铸分会秘书长许宏伟致开幕词,中国金属学会副秘书长苏天森讲话。中国钢铁协会副秘书长、中国金属学会副秘书长李世骏,北京科技大学教授蔡开科等分别作了精彩的专题报告。大会后就连铸工艺技术、质量控制、连铸电磁搅拌技术及连铸机设计和改造进行了分组讨论。大会有来自大学院所,工矿企业 170 多个单位的代表共 267人,录用论文 133 篇,出版了第八届全国连铸学术会议论文集。连铸分会每四年举办一次全国连铸学术年会,从 1981 年至今已举办了七届年会,年会对我国连铸技术发展和学术研究及生产进步起到了积极的推动作用。(本刊讯)14 中国冶金 第 17 卷