![点击分享此内容可以赚币 分享](/master/images/share_but.png)
基于I-GWO-BP神经网络的紫金山铜矿浮选回收率预测研究.pdf
《基于I-GWO-BP神经网络的紫金山铜矿浮选回收率预测研究.pdf》由会员分享,可在线阅读,更多相关《基于I-GWO-BP神经网络的紫金山铜矿浮选回收率预测研究.pdf(9页珍藏版)》请在咨信网上搜索。
1、基于 IGWOBP 神经网络的紫金山铜矿浮选回收率预测研究夏永涛1,马英强1,2,印万忠1,衷水平1,2,鲁军3,张德文3,詹殷权11.福州大学 紫金地质与矿业学院,福建 福州 350108;2.福建省新能源金属绿色提取与高值利用重点实验室,福建 福州 350108;3.紫金矿业集团股份有限公司,福建 上杭 364200中图分类号:TD952.1文献标识码:A文章编号:10010076(2023)03005109DOI:10.13779/ki.issn1001-0076.2023.03.005摘要为克服传统测量浮选回收率方式存在的低效率、滞后性等问题,结合紫金山硫化铜矿浮选厂生产情况,采用基于
2、MI(Mutual Information)互信息法对选厂原矿品位、丁铵黑药用量等浮选条件因子进行特征选择,在此基础上,建立了基于BP(Back Propagation)、GWOBP(Grey Wolf OptimizerBack Propagation)、IGWOBP(ImprovedGrey Wolf OptimizerBackPropagation)的三种浮选回收率预测模型,并选取紫金山硫化铜矿浮选车间生产数据进行神经网络训练与验证试验,分析了浮选回收率预测模型的准确性。结果表明:相较于基于 BP、GWOBP 的浮选回收率预测模型而言,基于 IGWOBP 的浮选回收率预测模型具有更大的相
3、关系数和更小的均方误差根,说明该模型泛化拟合能力更强,对浮选回收率的预测值在很大程度上逼近于真实值,预测精度更高。本研究结果可为实现浮选回收率高效、准确、自动的在线预测技术开发提供支持。关键词浮选回收率;互信息法;BP 神经网络;改进的灰狼算法 浮选回收率是衡量选矿生产质量的重要指标之一。通过人工化验、仪器测量、计算获得的实时回收率数据可帮助工人及时调整浮选条件进而改善浮选指标,但这种方式效率低、成本高、误差大。同时,通过传统的物理化学反应和内部因素来建立浮选回收率预测模型有很大难度。因此,一些学者通过软测量建模的研究方法,以期实现对浮选回收率的快速、自动、准确预测。尤其是近年来,人工神经网络
4、在浮选回收率预测方面应用逐渐变多。刘青等1在大量现场生产数据的基础上,分别采用多元线性回归和 BP 神经网络建模精准实现了对金矿选矿厂浮选回收率的预测。刘利敏等2利用 BP 神经网络对矿浆浓度、pH 值、充气量大小、药剂用量、泡沫层厚度与浮选回收率之间的数学关系进行了拟合,其输出值能以较高精度逼近实际样本数据。王衡嵩等3建立了磨矿细度、戊基黄药用量、碳酸钠用量、硫酸铜用量和硫化钠用量与浮选回收率的反向传播神经网络(BPNN)模型,通过验证实验表明:浮选回收率较响应曲面分析(RSM)模型有较好的改善。张勇等人4基于混沌蚁群算法预优化 BP 神经网络实现对作业回收率的精准预测。冉宇等人5对泡沫图像
5、进行多种方式处理来确定泡沫大小、统计泡沫大小分布规律,采用 PDF 泡沫概率统计方法进行浮选面积的统计,最后应用 BP 神经网络建立了回收率预测模型。目前来看,在人工神经网络中,尤其是 BP 神经网络实现对浮选回收率的预测更加普遍、高效。但是 BP 神经网络的初始权值、阈值和数据集本身的适配性对于网络预测有很大影响。因此,在数据集确定的情况下进一步加强 BP 神经网络自身泛化能力是实现浮选回收率精准预测的关键。为进一步加强浮选回收率预测模型的准确度,本文提出对原矿品位、磨矿细度、丁铵黑药用量、石灰用量和松醇油用量等浮选条件因子进行特征选择,引入基于 Tent 混沌映射初始种群,应用自定义非线性
6、收敛因子更新策略的改进灰狼算法对 BP 神经网络进行了优化设计,以紫金山硫化铜矿浮选厂生产数据为依据,构建了浮选回收率预测模型,并进行了基于 BP、GWOBP、IGWOBP 神经网络的浮选回收率验证实验对比分析。收稿日期:2023 02 11基金项目:国家自然科学基金项目(51804081);福建省自然科学基金项目(2019J01253);福州大学本科生科研训练计划项目(28212)通信作者:马英强(1983),男,博士,副教授,硕士生导师,主要从事矿物浮选技术与理论方面的研究工作,Email:。第 3 期矿产保护与利用No.32023 年 6 月Conservation and Utiliz
7、ation of Mineral ResourcesJun.20231基于浮选条件因子的特征选择矿山浮选车间工艺流程会受到多因素干扰,因此在大批量浮选工艺数据处理问题上,会出现数据集维度冗余的问题。选择关键的浮选条件因子可以在保证不降低浮选回收率预测准确性的前提下,简化数据集,降低模型学习难度,加快收敛速度。本文在紫金山硫化铜矿浮选车间记录 5 个浮选条件因子和浮选回收率,共 760 组浮选车间数据,利用互信息法对浮选条件因子的重要性进行了量化处理,并在Matlab2021a软件中,将计算结果可视化,结果如图 1 所示。0.00.10.20.30.4各浮选条件因子与浮选回收率之间的互信息值原矿
8、品位精矿品位磨矿细度石灰用量浮选条件因子松醇油用量丁铵黑药用量图 1各浮选条件因子与浮选回收率之间的互信息值Fig.1 The value of mutual information between each flotation condition factor and flotation recovery 由图 1 分析可得,在六个浮选条件因子中,原矿品位与浮选回收率之间的互信息值最高,磨矿细度、精矿品位与浮选回收率之间的互信息值最低,石灰用量、丁铵黑药用量、松醇油用量与浮选回收率之间的互信息值适中,而精矿品位、磨矿细度与浮选回收率之间的互信息值最低。据此可知,原矿品位、石灰用量、丁铵黑药用
9、量和松醇油用量四个浮选条件因子对浮选回收率的影响较大,因此选择这四个浮选条件因子作为浮选回收率预测模型的输入变量。2BP 神经网络模型浮选工业车间的数据体量大、偏差多,要实现各浮选条件因子与浮选回收率之间的数学建模,BP 神经网络的非线性映射能力和多特征融合能力具备一定的优势,且对样本中的异常值敏感度不高,相比于支持向量机(Support Vector Machine,SVM)等其他机器学习方法更具有适用性。2.1基于 BP 神经网络的浮选回收率预测模型2.1.1输入、输出层的设计随着紫金山硫化铜矿浮选厂的日常运转,浮选车间生产过程中积累了庞大的工业生产数据集,这为建立浮选回收率预测模型提供了
10、基础数据。为了构建与训练模型,本文选取了包含原矿品位、磨矿细度、石灰用量、丁铵黑药用量、松醇油用量和浮选回收率等浮选工艺参数在内的部分数据集,其中 600 组用于网络训练,50 组用于网络泛化能力验证,通过对网络不断地学习、训练,最终可建成较为可靠的浮选回收率预测模型。根据前述对于浮选条件因子的特征选择结果可知,在神经网络建模过程中,需要将原矿品位、石灰用量、丁铵黑药用量和松醇油用量这四个对浮选回收率影响较大的变量作为四个神经网络输入层的神经元,网络输出是浮选回收率,将浮选回收率这个因变量作为网络输出层的单一神经元。2.1.2隐层神经元的设计隐层神经元数量的确定对神经网络的性能影响很大,若隐层
11、结点数过少,网络的学习能力达不到要求;若过多,会增加网络结构的复杂性,且更易陷入局部最值。根据隐层神经元确定经验公式6(1):Z=X+Y+a(1)XYa式中:Z、分别表示输入、输出神经元数目,表示 110 的整数。XY由上节可知,为 4,为 1,从上述推断隐层神经元数目位于(3,12)。将此经验公式与 MATLAB 仿真结果结合起来,依次设置隐层神经元为 312 的整数,进行多次网络训练,记录每次训练网络所得到的均方误差,其结果如表 1 所示。由表 1 可知,当隐含层神经元个数为 10 时,均方误差 MSE 为 0.029 027,网络的训练效果最佳,据此建立的 BP 神经网络结构如图 2 所
12、示。2.2基于 GWOBP 神经网络的浮选回收率预测模型在基于 BP 神经网络的浮选回收率模型预测中,52 矿产保护与利用2023 年网络的初始值是随机的,因此在网络训练过程中,可能会得到局部最优的浮选回收率值。为了最大程度地使浮选回收率预测值逼近真实值,本文提出应用灰狼算法来优化基于 BP 神经网络的浮选回收率预测模型。灰狼优化算法是由 Mirjalili 等6在 2014 年提出的。根据该算法,本文设计将每个灰狼个体的位置信息比作基于 BP 神经网络的浮选回收率模型的初始权值、阈值,每迭代一次,将原矿品位、石灰用量、丁铵黑药用量和松醇油用量输入模型完成浮选回收率的预测,并将得到网络误差值标
13、记为灰狼个体的适应度值。按照该适应度值的排名划分 狼、狼、狼,剩下的灰狼个体为 狼,此时的 狼、狼、狼、狼的位置信息即为若干基于 BP 神经网络的浮选回收率模型的初始权值、阈值,此后 狼根据、和 狼的指示不断更新位置信息,完成抓捕、追踪、围猎和攻击猎物的过程,这也是获得基于 BP 神经网络的浮选回收率模型最优初始权值、阈值的过程,并将算法终止后的 狼的位置信息作为基于 GWOBP 神经网络浮选回收率预测模型的初始权值、阈值,将四个浮选条件因子输入进行训练。利用灰狼算法寻优基于 BP神经网络的浮选回收率模型初始权值、阈值过程如下:图 2BP 神经网络结构Fig.2 Structure diagr
14、am of BP neural network 假设灰狼种群大小为 N,捕猎范围为 d 维,其中第 i 只灰狼的位置信息可表示为:Xi=(Xi 1,Xi 2,Xi 3,Xi d),整个狩猎过程如下所示:包围猎物数学模型为:D=?CXp(t)X(t)?(2)X(t+1)=Xp(t)AD(3)Xp(t)tX(t)tC式中:为第 物的位置;为第 代灰狼个体的位置;为摆动因子。C=2r1(4)式中:r1为0,1之间的随机数;A 为收敛因子。A=2ar2a(5)a=2(1tTmax)(6)式中:Tmax为最大迭代数目;r2为0,1之间的随机数;a 为控制参数,其值随着迭代次数的增加而线性减小,amax=
15、2,amin=0。狩猎过程数学模型为:D=?C1X(t)X(t)?(7)D=?C2X(t)X(t)?(8)D=?C3X(t)X(t)?(9)表 1 不同隐含层节点数对应的网络均方误差Table 1 Network mean squared error corresponding to thenumber of nodes in different hidden layers隐含层节点数相应的均方误差30.037 56640.035 89950.036 76560.031 69970.029 84580.031 2590.035 295100.029 027110.031 674120.029
16、306第 3 期夏永涛,等:基于 IGWOBP 神经网络的紫金山铜矿浮选回收率预测研究 53 由式(7)(9)可计算出其余 灰狼个体与灰狼、的距离。X1(t)=X(t)A1D(10)X2(t)=X(t)A2D(11)X3(t)=X(t)A3D(12)由式(10)(12)计算出其余 灰狼位置更新的方向。X(t+1)=X1(t)+X2(t)+X3(t)3(13)由式(13)更新灰狼个体的位置。当灰狼算法寻优基于 BP 神经网络的浮选回收率模型初始权值、阈值的过程结束,将此时 狼的位置信息作为基于 GWOBP 神经网络浮选回收率预测模型的最优初始权值、阈值,并将原矿品位、石灰用量、丁铵黑药用量和松醇
17、油用量四个浮选条件因子输入网络,进行模型训练,进一步提升浮选回收率预测精度,使浮选回收率预测值更逼近真实值,也为实现浮选回收率高效、准确、自动的在线预测技术的开发提供支持。2.3基于 IGWOBP 神经网络的浮选回收率预测模型在基于 GWOBP 神经网络的浮选回收率模型预测中,最佳网络初始权值、阈值对于寻找全局最优的浮选回收率值至关重要,但是传统灰狼算法存在后期全局搜索能力较差、易陷入局部最优以及初始种群多样性不足等问题。因此针对这些局限性,本文提出基于 Tent 混沌映射初始化种群和非线性收敛因子的两种更新策略来帮助寻找最佳网络初始值,最大程度地使浮选回收率预测值逼近真实值。2.3.1基于
18、Tent 混沌映射初始化种群更新策略文献中多采用 Logistic 映射78产生混沌序列,但是 Logistic 映射产生的混沌序列均匀性较差,且多在(0,0.1)和(0.9,1)两个区间取值,因此计划使用运算速度更快、均匀性更好的 Tent 混沌映射910,取迭代次数为 100,具体在 Matlab 2021a 中仿真实现如图 3 所示。T(t+1)=T(t)/a,0 T(t)a1T(t)/a,a T(t)1 时狼群进行全局搜索,|A|1 时进行局部搜索,因为 A 始终在(a,a)之间变化,在式(7)中的参数 a 呈线性递减,导致算法后期容易出现全局搜索能力较差、结果易早熟、陷入局部最优等问
19、题11。为了实现在后期最优解邻域精确搜索,需保证参数 a 的变化呈现先平缓后急促的幅度,提出如式(15)的非线性收敛因子更新策略,其函数图像基于 Matlab 的实现如图 4 所示。01002003004005000.00.51.01.52.0迭代次数收敛因子更新值 I-GWO GWO图 4收敛因子变化趋势Fig.4 Change trend of convergence factor a=a(ori)a(ori)1(e1)2(etTmax1)2(15)a(ori)其中:为收敛因子的初始值,t 为当前迭代次数,Tmax为最大迭代次数。从图 4 可得,改进后的非线性收敛因子曲线呈现先平缓后急促的
20、趋势,便于后期算法陷入局部最优时跳出,提高寻找全局最优值的能力。在后期急促阶段,通过小幅改变灰狼个体位置信息,可在基于 IGWOBP神经网络的浮选回收率模型预测中,将原矿品位、石灰用量、丁铵黑药用量和松醇油用量作为输入变量,可以获得变幅较小的网络初始权值、阈值并进行训练,0.00.20.40.60.81.002468101214161820频数混沌映射值图 3Tent 混沌映射序列Fig.3 Sequence of Tent chaotic map 54 矿产保护与利用2023 年实现后期最优解邻域的精确搜索。结合上述两种改进策略组成的 IGWO 算法,可以优化基于 BP 神经网络的浮选回收率
21、预测模型,该模型预测流程如图 5 所示。3数据预处理3.1数据去噪处理紫金山铜矿浮选车间工艺流程受到多因素的影响,会导致大批量浮选数据存在异常值,这对于网络的训练会产生很大干扰。因此,本文引入了 Hampel滤波器对数据进行去噪处理,可以实现异常值的检测和删除,且不会过度平滑数据,保证原始数据的真实性。由前述图 1 分析可知,浮选回收率预测模型的输入变量为原矿品位、石灰用量、丁铵黑药用量和松醇油用量。因此,将各浮选条件因子数据去噪处理,对去噪前后的部分数据进行对比分析,结果如图 6 所示。由图 6 分析可得,将去噪前后的原矿品位、石灰用量、丁铵黑药用量和松醇油用量的样本点特征值进行对比分析发现
22、,去噪前的各浮选条件因子的样本点特征值离散化程度较高,而去噪后的各浮选条件因子的样本点特征值线性化程度更高,且离群值、异常值明显减少,特征值之间的量纲差异减小,这有利于后续对浮选回收率模型的训练,加快网络收敛,提高对浮选回收率预测的准确性。3.2数据归一化处理当浮选工艺参数特征值的范围差异较大时,若浮选回收率预测模型发生网络梯度更新时会处于来回震荡的状态,很大程度上影响着浮选回收率预测模型的训练效率、预测精度等,基于此要对浮选各工艺参数数据集进行归一化处理,消除各浮选工艺参数间的量纲差异,帮助浮选回收率预测模型实现快速收敛。用公式(16)可以实现对数据集的线性变换,使各浮选工艺参数的特征值映射
23、到 0,1 之间。x=xxminxmaxxmin(16)式中:x指的是归一化后的各浮选工艺参数的特征值,x 指的是各浮选工艺参数的特征初始值,xmin指的是各浮选工艺参数的最小特征初始值,xmax指的是各浮选工艺参数的最大特征初始值,实现归一化后的各浮选工艺参数数据如图 7 所示。由图 7 可知,归一化处理后的原矿品位、石灰用量、丁铵黑药用量、松醇油用量和浮选回收率等浮选工艺参数的特征值范围都稳定在 0,1 区间内,这样可以使浮选样本数据集中的奇异样本减少,并且可以有效改善网络训练效果、加快网络收敛速度,更大程度地使得浮选回收率网络预测值逼近于真实值。4网络模型的训练与验证实验为了对浮选回收率
24、预测模型进行训练与测试,本文预先采集了 2022 年 19 月的紫金山硫化铜矿浮选厂的 760 组生产数据,经特征参数选择和数据预处理之后,保留了 650 组包含原矿品位、石灰用量、丁铵黑药用量、松醇油用量和浮选回收率等浮选工艺参数的 图 5基于 IGWOBP 神经网络的浮选回收率预测模型流程Fig.5 Flowchart of flotation recovery prediction model based onIGWOBP neural network第 3 期夏永涛,等:基于 IGWOBP 神经网络的紫金山铜矿浮选回收率预测研究 55 归一化数据,随机选取其中 600 组数据集作为训练
25、集,50 组数据集作为验证实验测试集。4.1网络模型的训练本文对前述三种浮选回收率神经网络预测模型进行训练。在训练基于 BP 的浮选回收率预测模型之前,网络按照如表 2 所示设置基本参数,初次迭代时,随机生成网络权值、节点阈值,在样本点输出端形成输出值与期望值之间的误差信号,其沿着反向原路径传播,并按照学习函数对网络权值和节点阈值进行更新,不断地重复此过程,直到满足最大迭代次数或者学习精度后,停止训练,得到最终的网络权值、节点阈值,此时基于 BP 的浮选回收率预测模型进入验证实验环节。基于 GWOBP 的浮选回收率预测模型和基于 IGWOBP 的浮选回收率预测模型通过灰狼算法得到若干初始网络权
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 GWO BP 神经网络 紫金山 铜矿 浮选 回收率 预测 研究
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。