有色金属科学与工程  2015, Vol. 6 Issue (4): 81-85, 120
文章快速检索     高级检索
引用本文
0
鄢发明, 艾光华, 吴彩斌, 李晓东, 石志中, 毛文明. 钨矿石磨矿动力学研究[J]. 有色金属科学与工程, 2015, 6(4): 81-85, 120. DOI:10.13264/j.cnki.ysjskx.2015.04.017.
YAN Faming, AI Guanghua, WU Caibin, LI Xiaodong, SHI Zhizhong, MAO Wenming. On the grinding kinetics of tungsten ores[J]. Nonferrous Metals Science and Engineering, 2015, 6(4): 81-85, 120.DOI: 10.13264/j.cnki.ysjskx.2015.04.017.
钨矿石磨矿动力学研究[PDF全文]
鄢发明1, 艾光华1, 吴彩斌1 , 李晓东2, 石志中2, 毛文明2    
1. 江西理工大学江西省矿业工程重点实验室，江西 赣州 341000;
2. 湖南柿竹园有色金属有限责任公司，湖南 郴州 423037
基金项目：江西省自然科学基金资助项目（20122BAB216020）
通信作者：吴彩斌（1972-），男，博士，教授，主要从事磨矿过程优化与节能降耗方面的研究工作，E-mail：wushirle@sina.com
作者简介： 鄢发明（1989-），男，硕士研究生，主要研究方向为精确化磨矿与磨矿节能降耗，E-mail：1286225018@qq.com
收稿日期：2014-11-17
摘要：建立和预测磨矿过程的产品粒度分布数学模型有助于指导选矿厂磨矿实践.文中在实验室条件下对柿竹园钨矿石进行不同磨矿时间的分批磨矿试验.结果表明，矿粒在0.1 mm及其以下粒级符合一阶磨矿动力学，在0.1 mm以上粒级符合n阶磨矿动力学.一阶磨矿动力学参数k与粒径d呈幂函数关系，其动力学方程为R=R0exp[-（0.014 42+9.538 65d1.761 28）t]；n阶磨矿动力学参数中，k与粒径d呈幂函数关系，n与粒径d呈对数函数关系，其动力学方程为R=R0exp[-（-0.134 1+1.151 42d0.681 69）· t0.718 18-0.386 56lnd].利用建立的磨矿动力学方程计算出的理论值与磨矿实际结果进行比较，其吻合程度较高，为选矿厂通过调节磨机给矿量和分级机返砂量来控制磨矿产品粒度奠定理论基础.
关键词： 钨矿石    磨矿动力学    粒度分布    磨矿过程    
On the grinding kinetics of tungsten ores
YAN Faming1, AI Guanghua1, WU Caibin1 , LI Xiaodong2, SHI Zhizhong2, MAO Wenming2    
1. Jiangxi Key Laboratory of Mining Engineering, Jiangxi University of Science and Technology, Ganzhou 341000, China;
2. Hunan Shizhuyuan Nonferrous Metals Limit Corp, Chenzhou 423037, China
Abstract: The mathematical model of product particle size distribution of grinding process is established and predicted, which will be to the benefit of guiding mill grinding. Batch grinding experiments on Shizhyuan tungsten ores in different grinding time in the laboratory were made to show that particles conform to first order grinding kinetics when the particle size was no more than 0.1 millimeter, however, the particles conformed to n order grinding kinetics when the particle size was more than 0.1 millimeter. There was a power function relationship between the parameter k of first order grinding kinetics and diameter d, which distributed in function of R=R0exp[-(0.014 42+9.538 65d1.761 28)t]; at n order grinding kinetics, the parameter k and diameter d were consistent with power function relationship, and the parameter n and diameter d were consistent with logarithmic function relationship, which distributed in function of R=R0exp [-(-0.134 1+1.151 42d0.681 69)t0.718 18-0.386 56lnd]. Through the comparison with theoretical values which were calculated by grinding kinetics equations with the grinding actual results, actual values had high degree of coincidence with it. Therefore the study establishes a theoretical basis for controlling the grinding product particle sizes by adjusting the feed quantity of mill and sand return of grader.
Key words: tungsten ores    grinding kinetics    size distribution    grinding process    

磨矿作业在选矿厂中占有非常重要的地位，其能耗约占整个选矿厂能耗的50 %左右[1].磨矿是一个复杂的高能耗、低效率的作业，其过程还不是完全清楚，为了更好地进行磨矿作业，研究磨矿动力学很有必要[2-4].目前，对矿石磨碎的动力学模型有许多种，主要分为2类，一类是粒度及其相关量与时间及其相关量的动力学模型，另一类是总体平衡动力学模型[5].本文选用E.W.戴维斯提出的动力学模型来分析讨论钨矿石在磨矿过程中的行为，此模型研究的是被磨物料磨碎速率与其磨矿时间的关系，而物料的磨碎速率遵循磨矿动力学方程，磨矿动力学方程的参数k和n在一定程度上可以反映物料被磨碎的快慢以及为选择性磨矿和现场生产指标评价提供依据[6-15].磨矿动力学方程式还可以用来预测不同磨矿时间磨矿产品的粒度分布[16-18].在很多情况n阶磨矿动力学较一阶磨矿动力学更符合实际，因此大多情况下，此动力学模型均是用n阶磨矿动力学来分析矿石在磨矿过程中的行为[19-21]，而本文针对湖南柿竹园钨矿石特性，将其磨矿产品分为不大于0.1 mm粒级的细粒级产品和大于0.1 mm粒级的粗粒级产品，分别探讨其磨矿动力学规律并建立数学模型，为选矿厂磨矿实践提供指导.

1 试验 1.1 矿石性质

试验所用钨矿石取自湖南柿竹园多金属选厂，主要矿物为白钨矿、黑钨矿、假象及半假象黑钨矿，其中白钨矿与黑钨矿的比例为7:3.白钨矿密度为6.1 t/m3，解离度平行中等，普氏硬度系数f为4.5~5，属中等硬度；黑钨矿密度为7.30 t/m3，解离度平行完全，普氏硬度系数f为4~5.5，属中等硬度.钨矿石不规则抗压强度为400 kg/cm2左右，即换算成普氏硬度系数f为4左右.

1.2 试验条件

在XMB三辊四筒磨机上分别进行了钨矿石不同磨矿时间（2 min、3 min、4 min、5 min、6 min、7 min）的分批磨矿试验.所用磨筒尺寸为φ135 mm×146 mm，磨筒转速为85 r/min，磨矿介质为钢球，磨矿浓度为70 %，介质充填率为45 %，磨矿量为200 g.

1.3 粒度测定

采用筛孔尺寸为0.074 mm、0.1 mm、0.15 mm、0.3 mm、0.5 mm、1.0 mm、2.0 mm的套筛测定物料0.074~2.00 mm间的各粒级分布；采用激光粒度分析仪测定物料≤0.074 mm的各粒级分布.

1.4 试验方法

对不同磨矿时间下的磨矿产品粒度分布，用一阶磨矿动力学进行验证，并采用Origin软件进行拟合，求解不同粒级下磨矿动力学参数k值，观察并比较求解不同粒级的磨矿动力学参数k值时的拟合程度R2，确定磨矿产品符合一阶磨矿动力学的粒级范围，并建立一阶磨矿动力学方程式[22-24].对于磨矿产品不在一阶磨矿动力学粒级范围内的，同样采用Origin软件进行求解不同粒级的磨矿动力学参数k和n值，并在此粒级范围内建立n阶磨矿动力学方程式[25].最后分别运用所求得的各自粒级范围内的一阶或n阶磨矿动力学方程式进行理论计算，将计算结果与实际磨矿结果进行比较，来验证本研究所求解的磨矿动力学方程式的准确性，为选矿厂磨矿提供理论依据[26-28].

2 磨矿动力学原理

磨矿动力学是指被磨物料的磨碎速率与磨矿时间的关系，其n阶磨矿动力学方程式为：

(1)

式（1）中：R为磨矿为t时刻磨矿产品中粗粒级产率；R0为原矿中该粗粒级比率；t为磨矿时间；k为动力学参数；n为动力学参数，当n=0时称为零阶动力学，n=1时称为一阶动力学[12, 29].

对式（1）进行变换得：

(2)
(3)

式（1）中，动力学参数n和k分别是粒度d的函数，表达式如式（4）和式（5）所示[30]：

(4)
(5)

式（4）、式（5）中c0，c1，a0，a1及x1，x2分别为待定参数.

3 结果与讨论 3.1 一阶磨矿动力学 3.1.1 磨矿动力学参数k的求解

在磨矿条件为磨筒转速170 r/min，矿浆浓度为70 %，介质充填率为45 %下，依次对物料进行不同时长（2 min、3 min、4 min、5 min、6 min、7 min）的磨矿，结果如表 1所示.

表1 不同磨矿时间下磨矿产品筛上累积粒度分布/% Table 1 Cumulative oversize distribution of grinding products under different grinding time /%
点击放大

由表 1可知，0.1 mm及其以下粒级筛上累积产率随着时间的延长逐渐减少，且趋势较平缓；0.1 mm以上粒级筛上累积产率在5 min之前减少的较快，而5 min之后减少的较慢.这表明粗粒级矿粒相对于较细粒级矿粒更好磨，随着磨矿时间的延长，各类粒级矿石裂缝、缺陷越少，磨细越困难，因而使得矿石的可磨性逐渐变差[8].

由表 1可计算出磨矿产品不同时间不同粒级的ln（R0/R）值，结果如表 2所示.

表2 不同磨矿时间下不同粒级的ln（R0/R）值 Table 2 ln (R0/R) of different particle sizes under different grinding time
点击放大

根据表 2绘制不同粒级的ln（R0/R）与时间t的关系曲线，如图 1所示.并根据式（1）对不同粒级的ln（R0/R）与时间t的关系曲线进行线性拟合，得到钨矿石不同粒级下的磨矿动力学参数k值，结果如表 3所示.

图 1 不同粒级下ln（R0/R）与时间t的关系曲线 Fig. 1 Relationship curve between the ln (R0/R) values of different particle sizes and grinding time t

表3 不同粒级的磨矿动力学参数k Table 3 Grinding kinetics parameters k of different particle sizes
点击放大

由表 3可知：随着粒级的减小，k值逐渐变小，矿粒被磨碎速率越小[5]，这是因为随着矿石粒径的减少，剩余矿粒表面裂缝相对较少，强度较高.矿粒在0.1 mm及其以下粒级的线性相关系数R2均在0.99以上，而在0.1 mm以上粒级的线性相关系数R2均在0.96以上，拟合程度不高.因此矿粒在0.1 mm及其以下粒级符合一阶磨矿动力学，而在0.1 mm以上粒级则符合n阶磨矿动力学.

3.1.2 一阶磨矿动力学方程的建立

在上述基础上，将表 3中不同粒级的磨矿动力学参数k在0.1 mm及其以下粒级进行曲线拟合，求出钨矿石的磨矿动力学参数k的表达式，结果如图 2所示.

图 2 磨矿动力学参数k与d的关系曲线 Fig. 2 Relationship curve between the grinding kinetics parameters k and d

由图 2可知，参数Ki呈幂函数的形式，可表示为：

(6)

利用Origin软件对式（6）进行曲线拟合，求得参数a0、a1和x，从而求得钨矿石的一阶磨矿动力学参数k的表达式为：

(7)

将式（7）代入式（1），可建立钨矿石在0.1 mm及其以下粒级的一阶磨矿动力学方程为：

(8)
3.2 n阶磨矿动力学 3.2.1 磨矿动力学参数n和k的求解

由上述分析可知，钨矿石在0.1 mm以上粒级时不符合一阶磨矿动力学.为验证其是否符合n阶磨矿动力学，由表 1中0.1 mm以上粒级数据通过式（3）进行计算，得到不同粒级的lg（lgR0/R）与t的关系，结果如表 4所示.

表4 不同磨矿时间下不同粒级lg（lgR0/R）值 Table 4 lg (lgR0/R) of different particle sizes under different grinding time
点击放大

根据表 4绘制不同粒级的lg（lgR0/R）与lgt的关系曲线，如图 3所示.根据式（2）对不同粒级的lg（lgR0/R）与lgt的关系曲线进行线性拟合，得到不同粒级的动力学参数n值，进一步利用式（3）计算动力学参数k值，结果如表 5所示.

图 3 不同粒级的lg（lgR0/R）与lgt的关系曲线图 Fig. 3 Relationship curve between the lg (lgR0/R) values of different particle sizes and lgt

表5 不同粒级的磨矿动力学参数n和k Table 5 Grinding kinetics parameters n and k of different particle sizes
点击放大

由表 5可以看出，随着粒级的减小，矿粒的动力学参数n值逐渐增大，k值逐渐减小.这是因为随着矿石粒径的减少，矿粒表面裂缝相对越少，矿粒强度越高，因而表现为随着矿石粒径逐渐减小，矿石的可磨性逐渐变差直至变化不大；随着矿石粒径的减少，矿石逐渐变脆，矿石的均匀性逐渐变差[16].

3.2.2 n阶磨矿动力学方程的建立

根据表 5中数据绘制动力学参数n和k与粒径d的关系曲线，并进行拟合，如图 4所示.

图 4 动力学参数n和k与粒径d的关系曲线图 Fig. 4 Relationship curve between grinding kinetics parameters n and k and diameter d

由图 4可知，曲线k=f（d）呈幂函数形式，而曲线n=f（d）更接近于对数形式，如式（9）和式（10）所示.

(9)
(10)

根据表 5数据，利用Origin软件分别对式（9）和式（10）进行曲线拟合，求得未知参数c0，c1，a0，a1和x，进而求得钨矿石的磨矿动力学参数k和n的表达式分别为：

(11)
(12)

将式（11）和式（12）代入式（1）得在0.1 mm以上粒级的n阶磨矿动力学方程为：

(13)
3.3 磨矿动力学方程的验证

为了验证所建立的磨矿动力学方程对柿竹园钨矿石磨矿的适用性，将式（8）和式（13）建立的0.1 mm及其以下粒级与0.1 mm以上粒级的磨矿动力学方程进行计算，得到各自粒级范围内不同磨矿时间下的磨矿产品筛上累积产率理论值，结果如表 6所示.

表6 不同磨矿时间下的磨矿产品筛上累积产率理论值/% Table 6 Theoretical values of cumulative oversize distribution of grinding products under different grinding time /%
点击放大

将表 6数据与表 1数据进行对比，得到利用磨矿动力学方程计算的不同磨矿时间下磨矿产品不同粒级的筛上累积产率理论值与实际磨矿试验值的偏差值，结果如表 7所示.

表7 筛上累积产率理论值与实际磨矿试验值的偏差值△R /% Table 7 Deviation values△R of theoretical values and actual grinding test values of cumulative oversize /%
点击放大

从表 7可以看出，所建立的磨矿动力学方程对柿竹园钨矿石实际磨碎规律的吻合程度较高，与其它磨矿动力学模型相比，该磨矿动力学方程相对结构简单.因此，利用所建立的磨矿动力学方程为选矿厂预测磨矿产品粒度分布及进一步为通过调节磨机给矿量和分级机返砂量来控制磨矿产品粒度分布情况提供了理论依据.

4 结论

1）对柿竹园钨矿石而言，矿粒在0.1 mm及其以下粒级符合一阶磨矿动力学，在0.1mm以上粒级则符合n阶磨矿动力学.

2）一阶磨矿动力学参数k与粒径d呈幂函数关系，其一阶磨矿动力学方程为R=R0exp[-（0.014 42+ 9.538 65d1.761 28）t].

3）n阶磨矿动力学参数中，k与粒径d呈幂函数关系，n与粒径d呈对数函数关系，其n阶磨矿动力学方程为R=R0exp[-（-0.134 1+1.151 42d0.68169）t0.71818-0.38656lnd].

4）所建立的磨矿动力学方程相对结构简单，而且所得到的磨矿产品筛上累积产率理论值与磨矿试验值吻合程度较高，可以用来指导选矿厂磨矿实践.

参考文献
[1] 马天雨, 桂卫华, 王雅琳, 等. 基于分批试验的工业球磨机粒级分布预测模型[J]. 北京工业大学学报, 2012, 38(9): 1281–1286.
[2] Tönshoff H K, Friemuth T, Becker J C. Process monitoring in grinding[J]. CIRP Annals-Manufacturing Technology, 2002, 51(2): 551–571. DOI: 10.1016/S0007-8506(07)61700-4.
[3] Mizonov V E, Berthiaux H, Zhukov V P, et al. Application of multi -dimensional markov chains to model kinetics of grinding with internal classification[J]. International Journal of Mineral Pro -cessing, 2004, 74(增刊1): 307–315.
[4] 刘全军, 姜美光. 碎矿与磨矿技术发展及现状[J]. 云南冶金, 2012, 41(5): 21–28.
[5] 王志辉. 建陶原料粉磨特性的试验研究[J]. 福建建材, 2004(1): 17–19.
[6] 段希祥. 磨矿动力学参数与磨矿时间的关系研究[J]. 昆明工学院学报, 1988, 13(5): 23–33.
[7] 侯英, 丁亚卓, 印万忠, 等. 磨矿动力学参数对磨矿速度的影响[J]. 东北大学学报, 2013, 34(5): 708–711.
[8] 李茂林, 崔瑞, 王非, 等. 阶段磨选磨矿粒度划分的理论分析与计算[J]. 矿冶工程, 2011, 31(1): 27–32.
[9] 刘全军, 唐荣, 文书明. 石英-磁铁矿体系选择性磨细研究[J]. 昆明工学院学报, 1994, 19(3): 85–89.
[10] 田金星. 石墨及其混合物料的磨矿动力学行为[J]. 中国有色金属学报, 1996, 6(4): 47–50.
[11] 田金星, 杨华, 尚保平. 石墨-方解石混合物料方解石组分的碎裂参数特性[J]. 金属矿山, 2000(2): 41–43.
[12] 叶贤东, 文书明, 张文彬. 超临速磨机磨矿动力学研究[J]. 有色金属, 2004, 56(1): 83–85.
[13] 吴明珠, 张成名. 磨矿动力学在磨矿操作中的应用[J]. 有色金属, 1989(5): 10–15.
[14] 侯英, 印万忠, 丁亚卓, 等. 不同破碎方式下产品磨矿特性的对比研究[J]. 有色金属, 2014(1): 5–8.
[15] 印万忠, 吴凯, 王纪镇, 等. 破碎方式对紫金山铜金矿石可磨性及浮选的影响[J]. 福州大学学报, 2014, 42(2): 321–326.
[16] 王力, 张常法, 张军, 等. 助磨剂对煤沥青磨矿动力学的影响[J]. 山东科技大学学报, 2008, 27(5): 23–26.
[17] 刘晓, 徐超, 李波. 煤沥青水浆干法制备的磨矿试验研究[J]. 广东化工, 2011, 39(23): 82–84.
[18] 王舒娅, 龙光明, 祁米, 等. 磨矿时间对天青石粒度和锶质量分数分布的影响[J]. 化工矿物与加工, 2011(7): 8–10.
[19] 王敏, 丁浩. 搅拌磨湿法超细磨矿过程动力学的研究[J]. 中国非金属矿工业导刊, 1999(2): 21–23.
[20] 陈炳辰, 周凌嘉, 仲崇波. 邦德功指数模拟计算方法的研究[J]. 金属矿山, 1990(8): 36–39.
[21] 仲崇波, 王成功, 陈炳辰. 邦德棒磨功指数模拟计算方法的研究[J]. 中国矿山, 1997, 6(2): 48–52.
[22] Ahmadi R, Hashemzadehfini M, Amiri Parian M. Rapid determination of bond rod-mill work index by modeling the grinding kinetics[J]. Advanced Powder Technology, 2013(24): 441–445.
[23] Bozkurt V, Özgür I. Dry grinding kinetics of colemanite[J]. Powder Technology, 2007, 176: 88–92. DOI: 10.1016/j.powtec.2007.02.015.
[24] Ipek H. The effects of grinding media shape on breakage rate[J]. Minerals Engineering, 2006(19): 91–93.
[25] Ian Bailon-Poujol, Jean-Paul Bailon, Gilles L'Espérance. Ballmill grinding kinetics of master alloys for steel powder metallurgy applications[J]. Powder Technology, 2011, 210: 267–272. DOI: 10.1016/j.powtec.2011.03.028.
[26] Heekyu Choi, Woong Lee, Seongsoo Kim. Application of grinding kinetics analysis of inorganic powders by a stirred ball mill[J]. Korean Journal of Chemical Engineering, 2009, 26(6): 1806–1812. DOI: 10.1007/s11814-009-0239-y.
[27] Fuerstenau D W, De A, Kapur P C. Linear and nonlinear particle breakage processes in comminution[J]. International Journal of Mineral Processing, 2004, 74: 317–327. DOI: 10.1016/j.minpro.2004.04.006.
[28] Teke E, Yekeler M, Ulusoy U, et al. Kinetics of dry grinding of industrial minerals: calcite and barite[J]. International Journal of Mineral Processing, 2002, 67(1): 29–42.
[29] 段希祥. 碎矿与磨矿[M]. 北京: 冶金工业出版社 , 2012.
[30] 陈丙辰. 磨矿原理[M]. 北京: 冶金工业出版社 , 1989.