铁矿氧化球团回转窑三维温度场仿真模型

范晓慧，李俊，陈许玲，王祎，甘敏

(中南大学 资源加工与生物工程学院，湖南 长沙，410083)

我国的钢铁工业发展迅速，2011年我国钢产量已达到6.83×108t。目前，高碱度烧结矿配加酸性球团矿和(或)块矿是国内外公认的合理炉料结构。与烧结相比，球团生产过程能耗低，产品的强度大、品位高，其产量从2005年的0.59×108t提高到2010年的1.4×108t，在高炉炉料中的比例从10%提高到15.45%。链篦机-回转窑氧化球团生产工艺具有原料适应性强、生产规模大、生产的产品质量均匀且强度高、可用煤作燃料等众多优点，在我国得到迅速发展，2010年其产能占球团生产总量的56.51%，单机规模最高达5×106t/a。链篦机-回转窑氧化球团生产中，球团矿质量与回转窑的温度有直接关系，然而，由于回转窑处于高温及旋转状态，在线直接测量温度很困难，因此，窑内温度分布的研究一直是回转窑研究的重点。国内外的学者已对此进行一系列的研究，如邱夏陶等[1-5]分别建立了回转窑温度场模型。但目前所取得的研究成果大多是一维模型，不能全面地反映窑内的温度分布。为此，本文作者建立回转窑三维温度场模型，预测回转窑窑内球团和气体温度分布，以便为指导回转窑生产提供理论依据。

1 回转窑三维温度场模型

1.1 传热和传质分析

回转窑内煤粉燃烧产生大量的热烟气，在烟气和球团逆向运动的过程中，烟气以对流和辐射的方式对球团和窑壁进行加热，窑内壁以辐射和传导的方式对球团进行加热，窑外壁则以对流和辐射的方式向环境的散热。同时，回转窑内球团也在不断的运动中，其在横截面内的运动分为表面活动层和固定层 2个区域[6-9]，如图1所示。固定层中球团仅作圆周运动，球团之间并无相对运动，固定层球团的运动可表示为u=ωr[10]；表面活动层内由于球团的随机滚落运动，球团混合现象较强烈，需采用质量守恒和动量守恒定律进行运动分析。

由于烟气的扩散速度快，温度在横截面内相对均匀，因此，可以用轴向传热模型来描述烟气在各个横截面的温度分布。而对于球团，由于表面活动层和固定层的存在，球团存在径向和轴向的温度梯度，可以用二维传热模型描述横截面的温度分布，再结合用轴向传热模型即可描述球团在各个横截面的温度分布，进而构成回转窑三维温度场模型。

图1 横截面球团流动示意图Fig.1 Flow of pellets in cross-section

1.2 三维温度场模型

1.2.1 模型的建立

(1) 轴向传热模型。轴向传热模型是在假设过程为稳态，忽略各个横截面间的轴向传热，烟气、窑壁温度及黑度各处均匀，球团黑度各处均匀的基础上建立的。

将回转窑沿轴向分成n个区域，任一区域内存在烟气能量平衡、球团能量平衡、窑壁能量平衡，其平衡方程式分别如下：

(2) 横截面球团的二维传热模型。在回转窑工艺中，由于球团产生沿轴向向前的运动主要发生在表面活动层，因此窑内轴向的能量传递也主要是由于表面活动层内球团的运动进行的。则表面活动层内球团的能量方程可表示为[11]：

固定层球团的能量方程可表示为：

1.2.2 模型系数的计算

(2) 球团、气体和窑壁间的辐射传热系数为：

(3) 球团和窑壁间的导热系数为：

(4) 窑外壁与环境之间的传热系数。窑外壁向周围环境的散热包括对流传热和辐射传热，总传热系数为2者之和，即：

其中：D为窑的内径，m；De为当量直径，m；Dsh为窑外壁直径，m；η为填充率；ug烟气速度，m·s-1；νg为烟气的运动黏度，m2·s-1；ω为窑转速，rad·s-1；σ为波尔兹曼常数，W/(m2·K4)；εg，εb，εw和εsh分别为烟气、球团、窑内壁、窑外壁的黑度；εgw，εgb和εwb为系统有效黑度，计算方法见文献[3]；λb和λg分别为球团、烟气的导热系数，W·m-1·K-1；ρb为球团的密度，kg·m-3；Cb为球团的比热容，J·kg-1·K-1；τ为球团与窑壁接触的时间，s；Rp为球团直径，m；δ为气膜厚度，m；γ为计算气膜导热引入的调节系数[12]；Gr为格拉晓夫数；Pr为普朗特数；Δh的取值可由红外热像仪所获得的实际球团温度来确定，作为对传热系数的修正值。

(5) 煤粉燃烧放热量为[13-14]：

其中：Li为区域i处火焰距喷嘴距离，m；Lf为火焰长度，m，其取值可参考文献[15]。

2 系统软件开发

采用Visual Studio.Net，Matlab以及OpenCV作为开发工具，结合 ADO数据库技术，根据某球团厂实际生产情况开发铁矿氧化球团回转窑三维温度场仿真系统。

系统具有以下功能:

(1) 根据现场实际工况，在线计算轴向烟气、球团的温度分布，并在系统中进行动态显示；可以直观地了解窑内热工状况。

(2) 系统显示窑内主要横截面(窑头、窑中和窑尾)的温度分布，使操作人员能更直观地了解球团料层内部的温度变化。

本文利用2011年8月28日—9月3日在某氧化球团厂回转窑检测的窑头烟气的数据进行系统验证。运行结果表明：系统仿真计算的窑头烟气温度误差在7%以内时，命中率在90%以上(如表1所示)，因此，本系统能很好地反映窑内的温度分布。

表1 窑头烟气温度命中率统计Table 1 Hit rate statistics of flue gas temperature at lower end of kiln

3 结论

(1) 回转窑窑内烟气和球团成逆向运动；在任一横截面内，烟气扩散快，混合充分；而对于球团，固定层中球团仅作圆周运动，球团之间并无相对运动，表面活动层内球团做随机滚落运动，球团混合比较强烈。

(2) 根据烟气和球团在回转窑内不同的运动特点，用轴向传热模型来描述窑内烟气温度分布，用横截面二维传热模型和轴向传热模型描述球团的温度分布，建立回转窑三维温度场模型。

(3) 采用Visual Studio.Net，Matlab以及OpenCV混合编程的方式开发回转窑三维温度场仿真系统，应用现场实际生产数据验证，模型计算烟气温度命中率在90%以上，说明系统能很好地反映窑内的温度分布。

[1]邱夏陶, 韩小良. 回转窑传热数学模型及其优化[J]. 钢铁,1994, 26(6): 66-70.QIU Xia-tao, HAN Xiao-liang. Mathematical model of heat transfer and optimum for rotary kiln[J]. Iron and Steel, 1994,26(6): 66-70.

[2]Martins M A, Oliveira L S. Modeling and simulation of petroleum coke calcination in rotary kilns[J]. Fuel, 2001, 80(11):1611-1622.

[3]马爱纯, 周孑民, 孙志强, 等. 氧化铝熟料窑内一维传热模型[J]. 冶金能源, 2004, 23(1): 23-26.MA Ai-chun, ZHOU Jie-min, SUN Zhi-qiang, et al. A one-dimensional thermal model for the alumina clinker kiln[J].Energy for Metallurgical Industry, 2004, 23(1): 23-26.

[4]易正明. 氧化铝回转窑热工分析与控制应用研究[D]. 长沙:中南大学能源科学与工程学院, 2007: 60-93.YI Zheng-ming. Application and study on thermal analysis and control of alumina rotary kiln[D]. Changsha: Central South University. School of Energy Science and Engineering, 2007:60-93.

[5]王春华, 陈文仲, 贾冯睿, 等. 回转窑内传热及燃烧过程的数值模拟[J]. 化工学报, 2010, 61(6): 1379-1384.WANG Chun-hua, CHEN Wen-zhong, JIA Feng-rui, et al.Numerical simulation of heat transfer and combustion in carbon rotary kiln[J]. Journal of Chemical Industry and Engineering(China), 2010, 61(6): 1379-1384.

[6]Dhanjal S K, Barr P V, Watkinson A P. The rotary kiln: An investigation of bed heat transfer in the transverse plane[J].Metallurgical and Materials Transactions B, 2004, 35(6):1059-1070.

[7]刘刚, 池涌, 蒋旭光, 等. 颗粒物料在回转窑内的运动特性模型[J]. 浙江大学学报(工学版), 2007, 41 (7): 1195-1200.LIU Gang, CHI Yong, JIANG Xu-guang, et al. Model of axial transport of particles in rotary kiln[J]. Journal of Zhejiang University(Engineering Science), 2007, 41(7): 1195-1200.

[8]王春华, 陈文仲, 梁增英, 等. 炭素回转窑内物料运动形式的确定[J]. 工业炉, 2008, 30(4): 34-36.WANG Chun-hua, CHEN Wen-zhong, LIANG Zeng-ying, et al.Determination of solid particles motion types in rotary kilns[J].Industrial Furnace, 2008, 30(4): 34-36.

[9]李少华, 张轩, 张立栋, 等. 回转窑内颗粒轴向混合运动的数值模拟[J]. 中国粉体技术, 2011, 17(2): 23-26.LI Shao-hua, ZHANG Xuan, ZHANG Li-dong, et al. Numerical simulation of particles axial mixing in rotary kiln[J]. China Power Science and Technology, 2011, 17(2): 23-26.

[10]Boateng A A, Barr P V. A thermal model for the rotary kiln including heat transfer within the bed[J]. International Journal of Heat and Mass Transfer, 1996, 39(10): 2131-2147.

[11]张志霄. 废轮胎回转窑热解特性及应用研究[D]. 杭州: 浙江大学热能工程研究所, 2004: 116-117.ZHANG Zhi-xiao. Investigation of waste tire pyrolysis in rotary kiln and the application project[D]. Hangzhou: Zhejiang University. Institute for Thermal Power Engineering, 2004:116-117.

[12]Martin H. Theories of heat transfer coefficient heating wall and particle[J]. Drying Tech, 1982,15(2): 179-183.

[13]刘方, 宋正昶, 杨丽, 等. 水泥回转窑温度分布研究[J]. 冶金能源, 2009, 28(5): 21-25.LIU Fang, SONG Zheng-chang, YANG Li, et al. Study of the temperature distribution in cement rotary kilns[J]. Energy for Metallurgical Industry, 2009, 28(5): 21-25.

[14]雷先明, 肖友刚, 蒋寿生. 耦合工况下回转窑热工过程综合数学模型[J]. 化学工程, 2009, 37(10): 18-21.LEI Xian-ming, XIAO You-gang, JIANG Shou-sheng.Integrated mathematical model of thermal processes of rotary kiln under coupling condition[J]. Chemical Engineering(China),2009, 37(10): 18-21.

[15]Beer J M, Chigier N A. Combustion Aerodynamics[M]. London:Applied Science Publishers, 1972: 39-48.