2016, Vol. 37
2.国土资源部 多金属矿评价与综合利用重点实验室，河南 郑州 450006;
3.中国地质科学院 郑州矿产综合利用研究所，河南 郑州 450006
2.Key Laboratory for Polymetallic Ores′ Evaluation and Utilization, Ministry of Land and Resources, Zhengzhou 450006, China;
3.Zhengzhou Institute of Multipurpose Utilization of Mineral Resources, China Geological Survey, Zhengzhou 450006, China
磨矿作业不仅能耗及材料消耗高，而且产品质量直接影响到后续选别作业的指标[1].磨矿过程是矿料粒度减小的过程，通过磨矿动力学来研究粒度减小的规律，可以计算磨机生产率，并对循环负荷的影响和实际磨矿过程的情况等做出理论上的分析和判断[2-3].
高压辊磨机是基于料层粉碎原理设计的一种高效率超细碎设备，在金属矿山逐渐得到广泛应用[4-6].前期研究发现高压辊磨机粉碎产品较常规破碎产品内部结构疏松，裂纹数量多，比表面积大，所需的磨矿时间短，可磨性好[7].本文在前期研究的基础上，对不同破碎方式下的磨矿动力学参数和磨矿过程规律进行了对比分析.
1 磨矿动力学理论分析磨矿过程可以视为一个粒度离散-时间连续的动态平衡过程，磨矿动力学数学模型表达式如式(1)所示，在实验室条件下，碎裂速率参数和碎裂分布参数可以通过单粒级批次磨矿试验来确定[8-10].
|
(1) |
式中：mi(t)为经过时间t后第i粒级残留物的质量分数，%；Ki(t)为经过时间t后第i粒级碎裂选择参数，%·min-1；bi, j为经过时间t后第i粒级碎裂分布参数，无单位.
磨矿环境的变化所产生的任何一个bi, j都会产生一个速率参数K的假值，K和bi, j之间的耦合性非常强.鉴于此，粒级累计碎裂速率函数经常被采用，该碎裂速率函数并不需要分布函数bi, j就可以独立求解[11]，其表达式见式(2)：
|
(2) |
式中：R为经过t时间后粗级别残留物的质量分数，%；t为磨矿时间，min；k为比例系数，决定于磨矿条件；负号表示粗级别减少.
对式(2)积分可得到式(3)：
|
(3) |
式(3)就是一阶磨矿动力学方程.但是通过试验验证发现，m阶动力学方程更符合实际情况，m阶动力学方程为
|
(4) |
将式(4)代入式(2)，可以得到磨矿速度的方程：
|
(5) |
对磨矿速度v求导数，可以得到磨矿加速度a的表达式：
|
(6) |
上述式中：t为磨矿时间，min；R0为被磨物料中粗级别的原始质量分数，%；k，m为比例系数；v为时间为t时的磨矿速度，%/min.
2 试验设备与样品试验选用CLM-25-10型高压辊磨机(HPGR)和实验型圆锥破碎机(CC)(动锥直径为100 mm).试验样品为鞍山式赤铁矿，有用矿物主要是赤铁矿，含有少量磁铁矿，TFe品位为27.02%.原矿d80为18.5 mm.
磨矿机型号为XMB-70A型三辊四筒球磨机，磨矿时间t分别为3, 6, 9, 12, 15 min，每个磨矿产品筛分分析后进行粒度分布拟合.
3 试验结果分析不同磨矿时间下磨矿产品的粒度分布如表 1所示.根据表 1构建不同破碎方式下的m阶动力学方程，按照式(4)对不同粒度的m和k值进行拟合，如表 2所示.
|
表 1 不同磨矿时间下磨矿产品的粒度分布 Table 1 Size distribution of grinding products in different time |
|
|
表 2 试验测试的磨矿动力学参数值 Table 2 Breakage rate parameters obtained from laboratory tests |
将表 2中的m，k值代入式(5)中，可以求出粒度瞬时磨矿速度，粒级瞬时磨矿速度(vi)也是视为粒度瞬时磨矿速度之差，以v(0.2~0.15)为例：
|
应用MATLAB软件对粒级瞬时磨矿速度进行可视化分析，如图 1~图 7所示.
|
图 1 -2.0+0.2 mm磨矿速度曲线 Fig.1 Grinding specific rate curves for -2.0+0.2 mm size |
|
图 2 -0.2+0.15 mm磨矿速度曲线 Fig.2 Grinding specific rate curves for -0.2 +0.15 mm size |
从图 1~图 7可以看出，除上限粒级(i0)(-2.0+0.2 mm)和下限粒级(in)(-0.031 mm)外，ij粒级磨矿速度曲线(图 2~图 6)理论上均存在两个极值点，两个极值时间点(t1，t2)将磨矿过程分为三个阶段.
|
图 3 -0.15+0.105 mm磨矿速度曲线 Fig.3 Grinding specific rate curves for -0.15+0.105 mm size |
|
图 4 -0.105+0.074 mm磨矿速度曲线 Fig.4 Grinding specific rate curves for -0.105+0.074 mm size |
|
图 5 -0.074+0.043 mm磨矿速度曲线 Fig.5 Grinding specific rate curves for -0.074+0.043 mm size |
|
图 6 -0.043+0.031 mm磨矿速度曲线 Fig.6 Grinding specific rate curves for -0.043+0.031 mm size |
|
图 7 -0.031 mm磨矿速度曲线 Fig.7 Grinding specific rate curves for -0.031 mm size |
当磨矿时间0 < t≤t1，aij>0，vij正向骤增.当t=t1时，加速度aij=0，vij达到正向最大值.就磨碎概率(PG)而言，i< j粗粒级的PG( < j)远大于ij粒级的PG(j).该阶段可视为粗粒级快碎磨碎阶段.
当磨矿时间t1 < t < t2时，aij < 0，vij正向减小后负向增大.该阶段i< j粒级继续被磨碎的同时ij粒级矿粒开始得到有效磨碎，ij粒级含量进入一个有进有出的动态变化过程.PG( < j)降低而PG(j)增加.当vij=0，PG( < j)=PG(j)，此时，ij粒级的含量变化到达动态平衡点，ij粒级的含量最大，设此时对应的时间点为t0.当磨矿时间t1 < t < t0时，PG( < j)>PG(j)，ij粒级的含量保持增加的状态；当磨矿时间t0 < t < t2时，PG( < j) < PG(j)，ij粒级的含量进入减少的状态.当t=t2时，加速度aij=0，vij达到负向最大值，i< j粗粒级矿粒绝大多数已被磨碎，PG(j)同时达到最大.该阶段可视为中粒级动态碎磨阶段.
当磨矿时间t>t2时，由于i< j粗粒级矿粒绝大多数被磨碎，ij粒级增加量约为0，ij粒级进入单一的磨碎过程，粒级含量逐渐减少，|vij|逐渐降低至0.该阶段可视为中粒级单一碎磨阶段.
令各粒级的磨矿速度v=0, 加速度a=0，可以求出不同破碎产品磨矿速度曲线的t1, t2和t0，如表 3所示.
|
|
表 3 不同破碎方式下磨矿速度曲线的t1, t2和t0 Table 3 t1, t2 and t0 values in different crushing modes |
从表 3可以看出，各个粒级的t1 < < t2，这说明矿粒在磨矿过程中绝大多数时间处在中粒级动态碎磨阶段.随着粒级范围的降低，t2和t0时间逐渐延长；而t1时间先延长，当粒级范围小到-0.043+0.031 mm后缩短.
此外，与圆锥破碎机产品相比，高压辊磨机产品各个粒级的t1, t2和t0显著缩短，这说明高压辊磨机更早进入中粒级动态磨矿阶段，同时也会更早出现细粒级过粉碎(t1(HPGR)=1.013 4 < < t1(CC)=7.873 8).
4 结论1)磨矿过程中“粗粒级快速碎磨-中粒级动态碎磨-中粒级单一碎磨”三个阶段与给料的破碎方式无关.粗粒级快速碎磨阶段中，i < j粗粒级矿粒被快速磨碎后导致i≥j粒级含量急剧增加，而ij粒级矿粒并未得到有效粉碎.
中粒级动态碎磨阶段中，ij粒级含量进入一个有进有出的动态变化过程.中粒级单一碎磨，ij粒级含量逐渐减少，直至磨矿速度为0.
2)试验中发现，-0.043+0.031 mm粒级要比-0.043 mm+0.105粒级更早进入碎磨状态，导致下限粒级(-0.031 mm)含量在磨矿初期就会快速增加，容易出现细粒级过粉碎.
| [1] |
段希祥.
碎矿与磨矿(第三版)[M]. 北京: 冶金工业出版社, 2013 .
( Duan Xi-xiang. Crushing and grinding (3rd ed)[M]. Beijing: Metallurgical Industry Press, 2013 . ) |
| [2] | Matijasic G, Kurajica S. Grinding kinetics of amorphous powder obtained by sol-gel process[J]. Powder Technology , 2010, 197 (3) : 165–169. DOI:10.1016/j.powtec.2009.09.010 |
| [3] | Ozkan A, Yekeler M, Calkaya M. Kinetics of wet grinding of zeolite in a steel ball mill in comparison to dry grinding[J]. International Journal of Mineral Processing , 2009, 90 (1/2/3/4) : 67–73. |
| [4] | Maxton D, Morley C, Bearman R. A quantification of the benefits of high pressure rolls crushing in an operating environment[J]. Minerals Engineering , 2003, 16 (9) : 827–838. DOI:10.1016/S0892-6875(03)00179-1 |
| [5] | Van der Meer F P, Gruendken M. Flowsheet considerations for optimal use of high pressure grinding rolls[J]. Minerals Engineering , 2010, 23 (9) : 663–669. DOI:10.1016/j.mineng.2009.09.012 |
| [6] | Ozcan O, Aydogan N A, Benzer H. Effect of operational parameters and recycling load on the high pressure grinding rolls performance[J]. International Journal of Mineral Processing , 2015, 136 (1) : 20–25. |
| [7] | Han Y X, Liu L, Yuan Z T. Comparison of low-grade hematite product characteristics in a high pressure grinding roller and jaw crusher[J]. Mineral & Metallurgical Processing , 2012, 29 (2) : 75–80. |
| [8] | Austin L G, Shoji K, Bell D. Rate equations for non-linear breakage in mills due to material effects[J]. Powder Technology , 1982, 31 (1) : 127–133. DOI:10.1016/0032-5910(82)80014-4 |
| [9] | Fuerstenau D W, Phatak P B, Kapur P C, et al. Simulation of the grinding of coarse/fine (heterogeneous) systems in ball mill[J]. International Journal of Mineral Processing , 2011, 99 (3) : 32–38. |
| [10] | Fuerstenau D W, Abouzeid A Z M, Phatak P B. Effect of particulate environment on the kinetics and energetics of dry ball milling[J]. International Journal of Mineral Processing , 2010, 97 (1/2/3/4) : 52–58. |
| [11] | Acar C, Hosten C. Grinding kinetics of steady-state feeds in locked-cycle dry ball milling[J]. Powder Technology , 2013, 249 (8) : 274–281. |