2019, Vol. 40
2. 广西大学 资源环境与材料学院, 广西 南宁 530004
2. School of Resources Environment and Materials, Guangxi University, Nanning 530004, China
在锡石多金属硫化矿选矿实践中, 存在锡石过磨和硫化矿欠磨的矛盾, 因而, 寻找一种评价或者缓解锡石多金属硫化矿磨矿矛盾的途径尤其重要, 磨矿技术效率可以用来评价所需粒级的实现效率, 以及评价和优化锡石过磨和硫化矿欠磨过程.GRNN是建立在非参数估计基础上的一种非线性回归径向基网络, 即广义回归神经网络模型[1].近年来, 针对复杂问题的智能优化研究得到很大发展, 粒子群优化算法是文献[2-5]提出的一种新型算法, 应用较为广泛.通过粒子群算法对GRNN模型参数进行优化, 这一算法保留了GRNN非线性映射能力强和容错性高的优势, 并解决了GRNN拓扑结构不确定性问题, 具有较高的预测精度和稳定性.本文首次将GRNN模型和粒子群算法用于磨矿领域, 对锡石多金属硫化矿磨矿技术效率预测模型进行了研究和优化.
1 试验和讨论 1.1 原料矿物组成和成分试验矿样取自广西车河选矿厂, 矿样经干燥、破碎、混匀、筛分后制得若干矿样.锡石多金属硫化矿中的主要硫化矿成分是磁黄铁矿和闪锌矿, 铅矿物主要是脆硫锑铅矿, 锡矿物主要是锡石, 脉石矿物主要是云母和石英.
1.2 试验方法磨矿技术效率的定义为粒度a的磨矿效率减去过磨粒度b的磨矿效率.磨矿技术效率可以评价磨矿效果, 磨矿技术效率越高, 说明磨矿合格粒级含量越大, 磨矿效果越理想.磨矿技术效率数学表达式为
|
(1) |
式中:E为磨矿技术效率, %; γ为磨矿产物中小于粒度a的产率, %; γ1为给矿中小于粒度a的产率, %; γ2为给矿中小于粒度b的产率, %; γ3为磨矿产品中小于粒度b的产率, %.将锡石多金属硫化矿分为锡石和硫化矿二元结构, 并选用-0.23+0.075 mm为锡石的所需粒级, 选用-0.15+0.038 mm为硫化矿的所需粒级.
1.3 试验结果与分析 1.3.1 磨矿时间对锡石和硫化矿磨矿技术效率的影响磨矿时间为2, 4, 6, 8, 10 min, 干矿质量分数为65 %, 充填率为32 %, 入磨粒度为-1.7 mm时, 磨矿时间对锡石和硫化矿磨矿技术效率的影响如图 1所示.可知, 在锡石和硫化矿的磨矿时间为8 min时, 磨矿技术效率最大.
|
图 1 磨矿时间对锡石和硫化矿磨矿技术效率的影响 Fig.1 Effect of grinding time on grinding technology efficiency of cassiterite and sulfide ore |
干矿质量分数为35 %, 45 %, 55 %, 65 %, 75 %, 85 %, 磨矿时间为8 min, 充填率为32 %, 入磨粒度为-1.7 mm时, 干矿质量分数对锡石和硫化矿磨矿技术效率的影响如图 2所示.可知, 锡石和硫化矿在干矿质量分数为65 %时, 磨矿技术效率最大.
|
图 2 干矿质量分数对锡石和硫化矿磨矿技术效率的影响 Fig.2 Effect of dry ore mass fraction on grinding technology efficiency of cassiterite and sulfide ore |
充填率为22 %, 27 %, 32 %, 37 %, 42 %, 磨矿时间为8 min, 干矿质量分数为65 %, 入磨粒度为-1.7 mm时, 充填率对锡石和硫化矿磨矿技术效率影响如图 3所示.可知, 锡石和硫化矿在充填率为42 %时, 磨矿技术效率最大.
|
图 3 充填率对锡石和硫化矿磨矿技术效率影响 Fig.3 Effect of filling rate on grinding technology efficiency of cassiterite and sulfide ore |
试验结果表明, 锡石和硫化矿二元结构所对应的磨矿技术效率最优参数分别为磨矿时间8 min、干矿质量分数65 %, 充填率42 %.虽然锡石和硫化矿各自磨矿技术效率不同, 但其对应的最优参数一致, 因此, 分别研究锡石和硫化矿对应的磨矿技术效率预测模型也具有互相验证稳定性的作用.
2 粒子群算法优化GRNN回归模型 2.1 GRNN回归模型对于x, y两个随机变量, 分别对应的观测值为X和Y, 假设x与y联合概率密度为f(x, y), 若f(x, y)服从正态分布, 则数学期望表达式为
|
(2) |
式中:n为学习样本数目; σ为光滑因子; X为网络输入变量;Xi为第i个神经元对应的学习样本.
2.2 粒子群算法粒子群优化算法基本原理是基于自然界中鸟群通过信息交流和互相协作而实现集体性能最优化.把优化问题P看作S维空间中只有速度和位置的粒子, 所以粒子在S维空间中都对应着受优化函数控制的适应值F, 且所有粒子可以通过控制两个极值来控制和调节自身的位置(Li)和速度(Vi), 即个体极值(Lpi)和集体极值(Lgi).通过个体极值和集体极值粒子群不断更新速度和位置, 其更新公式分别为
|
(3) |
|
(4) |
式中:i, j分别为第i个粒子的第j纬度; k表示迭代次数; c1和c2为学习因子; r1(k), r2(k)为相互独立且是[0, 1]之间随机值.集体位置是集体最优位置时为最终解, 迭代终止.本文借助粒子群法对Matlab中GRNN模型光滑因子σ进行优化, 通过优化光滑因子使得GRNN实际值与输出值均方差最小.
3 基于GRNN-粒子群算法的磨矿技术效率模型预测与验证 3.1 基于GRNN-粒子群算法的预测模型学习样本在球磨过程中, 运动机制理论上可分为泻落式、抛落式和离心式三种运动形态, 然而, 实际磨矿过程却非常复杂, 影响因素众多, 主要包括磨机参数、矿石性质和操作因素等[6-9].针对磨矿技术效率, 本文选取磨矿时间、干矿质量分数和充填率作为研究主控因子, 锡石和硫化矿预测模型学习样本如表 1所示.
|
表 1 锡石和硫化矿预测模型学习样本 Table 1 Sample prediction models of cassiterite and sulfide ore |
为验证粒子群优化GRNN模型的可靠性和适用性, 基于上述磨矿优化试验, 选取2组具有代表性试验进行锡石和硫化矿磨矿技术效率模型预测, 如表 2所示.不同粒子群参数对锡石和硫化矿磨矿技术效率模型预测结果的影响如图 4和图 5所示.图中, k表示迭代次数, 分别取k1=10, k2=30, k31=50, k4=100;学习因子c=1.0, 1.5, 2.0, 2.5.
|
|
表 2 验证数据 Table 2 Validation data |
|
图 4 锡石预测模型第一组和第二组验证数据结果 Fig.4 Validation results of the first and second set data about cassiterite prediction model (a)—第一组；(b)—第二组. |
|
图 5 硫化矿预测模型第一组和第二组验证数据结果 Fig.5 Validation results of the first and second set data about sulfide ore prediction model (a)—第一组；(b)—第二组. |
由图 4和图 5可知, 在不同学习因子和迭代次数下, 锡石和硫化矿磨矿技术效率模型预测值与试验值差保持在±0.01范围内, 预测值与试验值表现出较高的一致性, 说明参数调节对预测结果影响较小, 人为干扰因素影响极小, 模型准确性与稳定性较高.这表明, 在磨矿实践过程中, 可通过粒子群优化GRNN建立预测模型, 有效调节适宜磨矿粒度.本研究为磨矿技术效率高效智能预测提供了一种新的方法.
4 结论1) 为高效评价和优化锡石多金属硫化矿磨矿实践中锡石过磨和硫化矿欠磨过程, 提出了一种基于Matlab和粒子群算法的磨矿技术效率预测模型及其优化方法.
2) 锡石和硫化矿二元结构所对应的磨矿技术效率最优参数皆为磨矿时间8 min、干矿质量分数65 %、充填率42 %, 虽然锡石和硫化矿各自磨矿技术效率不同, 但其对应的最优参数一致.
3) 在不同学习因子和迭代次数的条件下, 磨矿技术效率模型预测均值与试验值差在±0.01 %范围内, 预测结果与试验值表现出较高的一致性, 说明参数调节对预测结果影响较小, 人为干扰因素影响极小, 模型准确性与稳定性较高.
| [1] |
Specht D F.
The general regression neural network-rediscovered[J]. Neural Networks, 1993(6): 1033–1034.
|
| [2] |
Neelakantam V V, Tapabrata R.
Optimum design of Yagi-Uda antennas using computational intelligence[J]. IEEE Transaction on Antennas and Propagation, 2004, 52(7): 1811–1818.
DOI:10.1109/TAP.2004.831338 |
| [3] |
Motiian S, Aghababaie M, Soltanian-Zadeh H.Particle swarm optimization(PSO) of power allocation in cognitive radio systems with interference constraints broadband network and multimedia technology[C]//IEEE Congress on Evolutionary Computation.New Orleans: IEEE Xplore, 2011: 558-562.
https://ieeexplore.ieee.org/document/6155997 |
| [4] |
James K, Eberhart R C.Particle swarm optimization[C]//IEEE International Conference on Neural Networks.New York, 1995: 1942-1948.
https://ieeexplore.ieee.org/document/488968 |
| [5] |
Eberhart R C, Shi Y.Particle swarm optimization: developments, applications and resources [C]//IEEE Conference on Evolutionary Computation.Seoul: IEEE Computer Society, 2002: 81-86.
https://ieeexplore.ieee.org/document/934374 |
| [6] |
Fuerstenau D W, Phatak P B, Kapur P C, et al.
Simulation of the grinding of coarse/fine(heterogeneous) systems in a ball mill[J]. International Journal of Mineral Processing, 2011, 99: 32–36.
DOI:10.1016/j.minpro.2011.02.003 |
| [7] |
Matijasic G, Kurajica S.
Grinding kinetics of amorphous powder obtained by sol-gel process[J]. Powder Technology, 2010, 197: 165–168.
DOI:10.1016/j.powtec.2009.09.010 |
| [8] |
Bazin C, Obiang P.
Should the slurry density in a grinding mill be adjusted as a function of grinding media size[J]. Minerals Engineering, 2007, 20: 810–814.
DOI:10.1016/j.mineng.2007.01.017 |
| [9] |
Schmidt J, Plata M, Troger S, et al.
Production of polymer particles below 5 μm by wet grinding[J]. Powder Technology, 2012, 228: 84–88.
DOI:10.1016/j.powtec.2012.04.064 |