金氰化浸出过程混合建模及仿真算法分析研究.pdf

第4 2卷第4期(总第1 9 0期)2 0 2 3年8月湿法冶金H y d r o m e t a l l u r g yo fC h i n aV o l.4 2N o.4(S u m.1 9 0)A u g.2 0 2 3金氰化浸出过程混合建模及仿真算法分析研究莫文水(广西现代职业技术学院,广西 河池 5 4 7 0 0 0)摘要:针对某金氰化浸出过程,首先建立了多级动态机制模型,并在此基础上采用基于免疫原理的R B F神经网络数据模型的学习算法,计算实际生产过程中难以测量的动力学反应速度,从而估算动态机制模型中的未知参数,再与物料守恒方程串联,建立了预测金浸出率的串联型混合模型。其次与基于免疫原理的R B F神经网络建立的纯数据模型对比仿真来验证模型的有效性;最后根据实际生产过程的实际值与预测值之间的偏差进行模型更新,并通过仿真分析验证了更新机制的准确性。结果表明:串联混合模型大大提高了浸出过程浸出率的预测精度,模型的更新机制提高了模型的精准性与泛化能力。关键词:金;浸出;建模;仿真;算法中图分类号:T F 8 0 3.2 1 文献标识码:A 文章编号:1 0 0 9-2 6 1 7(2 0 2 3)0 4-0 4 2 9-0 7 D O I:1 0.1 3 3 5 5/j.c n k i.s f y j.2 0 2 3.0 4.0 1 7收稿日期:2 0 2 3-0 1-1 0基金项目:2 0 2 1年度广西高校中青年教师科研基础能力提升项目(2 0 2 1 KY 1 4 1 8)。第一作者简介:莫文水(1 9 7 6),男,壮族,本科,副教授,主要研究方向为信息可视化、U I设计研究。引用格式:莫文水.金氰化浸出过程混合建模及仿真算法分析研究J.湿法冶金,2 0 2 3,4 2(4):4 2 9-4 3 5.目前,湿法冶金浸出过程反应模型通常依据物料守恒、能量守恒及动力学方程建立1-5,但大部分研究依据理论仿真得出结论,与实际生产情况有偏差,试验参数无法固定,很难准确获取数据,易导致研究结论与实际生产情况不相符6-1 2。针对某金氰化浸出过程金浸出率无法精准预测的问题,首先建立了单级浸出动态模型,再基于单级模型建立多级浸出动态模型,并仿真验证多级浸出动态模型浸出率的变化趋势;其次用基于免疫原理的R B F神经网络数据模型的学习算法模拟动力学反应速度方程,估算动态机制模型中的未知参数,与物料守恒方程串联,建立了浸出率预测的串联混合模型,并根据模型的预测误差进行模型更新;最后通过在实际工业生产中的应用,验证串联混合模型的准确性。1 动态模型的建立1.1 单级和多级浸出动态模型的建立某金浸出工艺选择气力连续搅拌槽反应器(C S T R)1 3为浸出槽,以N a C N溶液作浸出剂,采用多级串联方式进行氰化浸出。提前假设搅拌槽内温度分布均匀、矿浆搅拌均匀、反应过程中不放热、忽略反应器中的物料隔离且浸出槽中矿浆p H恒定8,1 3-1 4。根据金氰化浸出过程的物料守恒和动力学反应,建立单级动态模型1 2。物料守恒方程如下:固相金守恒方程,dws,idt=Qs,iMs,i(ws,i1-ws,i)-r(A u)i;(1)液相金守恒方程,dwl,idt=Ql,iMl,iwl,i1-wl,i -Ms,iMl,ir(A u)i;(2)液相氰守恒方程,d(C N-)idt=Ql,iMl,i(C N-)i-1-(C N-)i)-Q(C N-)iMl,i-r(C N-)i。(3)式中:i浸出槽级数;ws,i固相金品位,m g/k g;wl,i液相金品位,m g/k g;(C N-)i液相氰根离子 质 量 浓 度,m g/L;Qs,i矿 浆 固 相 流 量,k g/h;Ql,i矿浆液相流量,k g/h;Q(C N-)i浸出槽中添加的氰化钠流量,m g/h;Ms,i浸出槽中固相滞留质量,k g;Ml,i浸出槽中液相滞留质 湿法冶金 2 0 2 3年8月量,k g;r(A u)i金 溶 解 速 度,m g/(k gh);r(C N-)i氰离子的消耗速度,m g/(k gh)。将动力学反应方程1 4与实际样本数据相结合,采用最小二乘法计算模型1 5参数:r(A u)i=(1.1 31 0-3-4.3 71 0-1 1(d-)2.9 3)(ws,i-ws,(d-)2.1 3(C N-)i0.9 6 1(O)0.2 2 8i;(4)r(C N-)i=1.6 91 08d-0.5 4 7-6.4 0(C N-)3.7 1i;(5)ws,=0.3 5 7(1-1.4 9 e-1.7 6 1 0-2d-)。(6)式中:d-动力学反应方程常数;(O)i溶解氧质量浓度,m g/L;ws,浸出渣中最小固相金品位,m g/k g。假设浸出过程稳定,基于Qs,i和Ql,i质量守恒可知:Qs,i=Qs,i-1;(7)Ql,i=Ql,i-1。(8)假设浸出槽中所有反应物充分混合,Qs,i和Ql,i之间的转化关系式为:Ql,i=Qs,i1w,i-1 。(9)式中:w,i矿浆浓度,k g/L。忽略物料隔离,假设Qs,i和Ql,i与矿浆有相同的i:i=ViQs,is+Ql,il10 0 0。(1 0)式中:s矿浆中固相密度,g/c m3;l矿浆中液相密度,g/c m3。由式(9)、(1 0)分别得固相、液相滞留量为:Ms,i=Qs,ii;(1 1)Ml,i=Ql,ii。(1 2)综上,式(1)(1 2)组成了单级浸出动态模型。使用O D E 4 5计算该动态模型得到ws,i,再由式(1 3)计算金浸出率:i=ws,i-1-ws,iws,i-11 0 0%。(1 3)基于单级浸出动态模型,将串联的每一级浸出槽的输出变量作为下一级浸出槽的输入变量,建立多级浸出动态机制模型,从而得到每一级浸出槽中矿浆的固相金品位,之后由式(1 4)计算整个浸出过程金总浸出率:N=ws,0ws,Nws,01 0 0%。(1 4)式中:ws,0矿浆初始固相金品位,m g/k g;ws,NN级浸出后矿浆的固相金品位,m g/k g。1.2 多级浸出动态模型的仿真分析采用六级浸出槽的浸出率试验数据对多级浸出动态模型进行模拟,并仿真分析不同输入变量对浸出率变化趋势的影响1 4-1 6。模型中相关变量取值见表1。表1 多级浸出动态模型中相关变量取值变量取值Qs,l/(k gh-1)25 0 0Q(C N-)l/(m gh-1)51 06ws,0/(m gk g-1)5wl,0/(m gk g-1)0(C N-)0/(m gk g-1)2 5 0d/m8 0w/(gc m-3)3 9V/m36 8s/(gc m-3)2.8l/(gc m-3)1(O)/(m gL-1)7 分别以矿石流量、矿浆质量浓度、氰化钠添加流量、溶解氧质量浓度及矿石平均粒径为输入变量,其余变量按照表1取值,仿真分析5个变量对金浸出率的影响,试验结果如图15所示。图1 矿石流量对金浸出率的影响图2 矿浆质量浓度对金浸出率的影响034第4 2卷第4期莫文水:金氰化浸出过程混合建模及仿真算法分析研究图3 氰化钠添加流量对金浸出率的影响图4 溶解氧质量浓度对金浸出率的影响图5 矿石平均粒径对金浸出率的影响 由图1看出:金浸出率随矿石流量增大而降低。因为矿石流量增大易导致矿浆在浸出槽中滞留时间变短,使矿浆无法与浸出剂充分反应,造成金浸出率降低。由图2看出:金浸出率随矿浆质量浓度增大而提高,矿浆质量浓度增至2 0k g/L后,金浸出率趋于稳定。这是由于矿浆质量浓度过大,黏度增大,导致金浸出不完全。由图3看出:金浸出率随氰化钠添加流量增加先快速升高;氰化钠添加流量增至5 0m g/h后,金浸出率升高幅度趋缓慢。由图4看出:金浸出率随溶解氧质量浓度增加而提高,溶解氧质量浓度增至5m g/L以后,金浸出率的增长幅度趋缓。由图5看出:浸出率随矿石平均粒径增大而降低。矿石平均粒径增大会导致化学反应不充分,使得浸出率减小。通过以上对动态机制模型中各输入变量对浸出率影响的仿真分析,帮助我们更好的分析了浸出率的影响因素情况及动态特性,为实际生产过程提供正确的理论指导。2 基于R B F神经网络的串联混合模型实际浸出过程中,动力学反应速度无法直接测量,因此,以浸出过程动态模型为基础,采用径 向 基 函 数 神 经 网 络(R B F NN)数 据 建 模法1 7-2 1建立浸出过程的动力学反应速度模型,估算浸出过程动态模型中的未知参数,再与物料守恒方程串联,从而建立预测浸出率的串联混合模型。2.1 径向基函数神经网络R B F神经网络是一种具有单隐层的三层前馈网络,包括输入层、隐含层和输出层。本研究主要讨论神经网络的训练算法,包括基础设置与公式。设输入层、隐含层、输出层的节点数分别为n、m、s,有N组输入输出样本对应(X(i),Y(i),(i=1,2,N),其中X(i)是n维向量,Y(i)是s维向量,输入层连接隐含层的权值设为1,n维径向基函数则作为隐含层节点的激活函数:(X)=(X(i)cj)。(1 5)式中:j=1,2,m;cj基函数的中心;欧几里得范数。径向基函数通常选取G a u s s i o n函数:(X)=(X(i)-cj)=e x p-122jX(i)-cj2 ,(j=1,2,m),(1 6)式中,j第j个以cj为中心的高斯函数宽的向量。一般输出可表示为:y(i)k=mj=1j k(i)j+k,(k=1,s)。(1 7)式中:j k隐含层节点j连接输出层节点k的权值;(i)j隐含层节点的输出值;k输出层节点k的阈值。式(1 7)的矩阵形式为:Y=W,RP(m+1),WR(m+1)s。(1 8)134 湿法冶金 2 0 2 3年8月式中:Y=y(1)1y(1)2y(1)sy(2)1y(2)2y(2)sy(P)1y(P)2y(P)s ,=(1)1(1)2(1)m1(2)1(2)22m1(P)1(P)2(P)m1 ,W=1 11 21s2 12 22sm1m2m s12s 。R B F神经网络算法是目前常用的混合学习算法之一,主要包含两个阶段。第一阶段是非监督学习阶段,需要得出径向基函数中的数据中心和宽度,通常采用K-m e a n s聚类算法。首先定义参数:P个样本X1,X2,XP,M个隐节点数,c(k)是第k次迭代的中心。其次 随 机 确 定M个 互 不 相 同 的 向 量c1(0),c2(0),cM(0)作为初始聚类中心;随后计算输入样本点与选择的聚类中心点的欧式距离并比较得出距离最小值的节点:dj(k)=Xp-cj(k),(p=1,2,P;j=1,2,M)。(1 9)令r表示该隐含层节点下标,根据式(2 0)与聚类中心的最小欧几里得距离对每一个输入样本Xp进行归类:r(Xp)=m i nXp-cj(k),(p=1,2,P)。(2 0)式中:r(Xp)Xp被分为r类。同理将 所 有 样 本 分 为M个 子 集U1(k),U2(k),UM(k),代表M个聚类域。采用竞争学习的规则对各个聚类中心进行调整:cj(k+1)=cj(k)+Xp-cj(k),j=rcj(k),jr(2 1)式中:学习率,取值范围为0,1。其次,令k=k+1,判断c的改变量是否小于预设阈值,是则利用聚类中心之间的距离确定径向基函数的扩展常数;否则通过式(1 9)重新计算欧氏距离并寻找距离最小值的隐含层节点,并按照以上步骤继续计算直至c的改变量小于预设阈值。第二阶段是通过最小二乘法计算连接权值并确定隐含层节点数目。设第j个隐含层节点的输出为p j=(Xp-cj),所以隐含层的输出矩阵为=(p j)PM,教师信号为d;设R B F神经网络的待定权值为W=(w1,w2,wM),隐含层输出矩阵为F(X)=W=d,利用的伪逆+得到W=+d,其中+=(T)-1T。混合学习算法虽有一定优点,但也存在R B F神经网络的隐含层数需要提前根据试验确定,其数值的不同会很大程度影响聚类的性能等缺点。2.2 基于免疫原理的R B F神经网络模型学习算法针对混合学习算法存在的不足,研究人员提出了免疫网络模型1 2,改进后可用于选择R B F神经网络隐含层数据中心(隐含层数),无需提前试验确定。假设一组向量作为输入数据X=x1,x2,xN,其中向量xi=xi1,xi2,xi pTRp(i=1,2,N)。确定R B F神经网络的数据中心需得到一个新的数据集作为聚类的输入数据,设该数据集为C=c1,c2,ch,cj=c1j,c2j,cpjTRp,j=1,2,h,其中h远小于N。首先定义抗原与抗体间的相互作用强度为亲和力,设输入数据xi和数据中心cj之间的亲和力为ai j;定义抗体与抗体间的相互作用强度为相似度,设数据中心ci与数据中心cj之间的相似度为si j;其中ai j与si j的表达式分别为:ai j=11+xi-cj,(i=1,2,N;j=1,2,h);(2 2)si j=11+ci-cj,(i,j=1,2,h)。(2 3)式中:欧几里得范数。当xi=cj时,ai j=1为最大;如果xi-cj0,则0ai j1。同理,当ci=cj时,si j=1为最大;如果ci-cj0,则0si js,保留该数据中心,反之则清除。3)将全部数据中心记忆集m 合并生成总数据集M。计算M中所有数据中心间的相似度si j,同理利用si js清除不同数据集中相似的数据中心,完成免疫系统中的网络抑制过程,生成最终数据集M。4)最后,随机选取一定数量的新数据中心替换C中对应数量亲和力低的数据中心,重复以上步骤。达到预设的最大递推次数后的M 即为所求的R B F数据中心集。综上,即可确定R B F神经网络的隐含层数。随后,采用递推最小二乘法求解线性方程组Y=W中的W,即可估计输出层权值。计算P(k):P(k)=1P(k-1)P(k-1)Z(k)ZT(k)P(k-1)+ZT(k)P(k-1)Z(k)。(2 5)式中:P(k)逆相关矩阵;Z(k)隐含层节点的输出值;遗忘因子,一般取01。更新网络权值:W(k)=W(k-1)+g(k)y(k)-ZT(k)W(k-1)。(2 6)式中:W(k)权值矢量;y(k)神经网络的期望输出。给出终止条件:J(k)=J(k-1)+12y(k)-ZT(k)W(k-1)。(2 7)式中:J(k)误差;误差满足的迭代停止值。若|J(k)-J(k-1)|,则此时的W为最终值;反之,转到计算P(k)值继续迭代。由此,确定了输出层的权值。2.3 仿真分析基于某金精炼厂浸出工艺采样数据,对串联混合模型进行仿真验证。为了方便对比分析,建立一个基于免疫原理R B F神经网络的纯数据模型(R B F NN数据模型)。将输入与输出数据之间拟合成非线性函数关系,模型的方程式如下:ak=f(Ql,Ql,0,Qs,Qs,0,Q(C N)0,d-,C0,Q(C N-)i),(i=1,2,3,4)。(2 8)对串联混合模型仿真验证的变量取值参见表1,免疫学习算法中R B F NN数据模型参数取值见表2。表2 免疫学习算法中R B F N N数据模型参数取值参数取值亲和力阈值d1相似度阈值s0.2 5百分比1 0%学习率0.9 5递推终止次数5R B F神经网络的中心数1 5R B F宽度5 两种模型预测效果的评价指标如下:1)最大绝对误差(MA E)MA E=m a xNi=1y-iyi;(2 9)2)均方根误差(RM S E)RM S E=1nni=1(yi-yi)2。(3 0)选择2 0 0组数据作为样本集,仿真分析浸出过程金浸出率的动态变化过程;再选取2 5组浸出达到稳态后的数据作为样本进行验证,结果如图6、7所示。图6 金浸出率动态变化的预测结果334 湿法冶金 2 0 2 3年8月图7 稳态金浸出率的预测结果 串联混合模型和R B F NN数据模型的误差对比结果见表3。可以看出:基于串联混合模型的MA E和RM S E都小于R B F NN数据模型的对应数据。说明串联混合模型预测值精度和泛化能力更高,在实际工业中应用效果会更好。表3 串联混合模型和R B F N N数据模型的误差对比结果模型动态MA E动态RM S E稳态MA E稳态RM S E串联混合模型0.0 1 430.0 0 860.0 0 960.0 0 41R B F NN数据模型0.0 2 670.0 2 310.0 1 730.0 1 283 串联混合模型的更新及仿真分析由于浸出过程具有时变性,模型精度会随时间延长而下降,若要精准预测,所建立模型需要具有更新能力。因此,需要通过对比预测值yp(t)和实际值y(t)之间误差判断是否需要进行模型更新。具体操作步骤为:1)初始化相关参数并设定标准误差值E=0.0 1;2)结合初始数据通过串联混合模型得到金浸出率预测值,将输入输出数据添加到历史数据中;3)再次计算得到金浸出率预测值;4)计算误差e(t)=|y(t)-yp(t)|,若e(t)E,则模型不变;反之,重复以上步骤。5)返回3)循环操作。选取2 0 0组数据作为样本集,仿真分析浸出过程金浸出率的动态变化过程;选取2 5组数据进行测 试,验 证 混 合 模 型 更 新 策 略 的 有 效 性。预测结果如图8、9所示。串联混合模型和改进混合模型误差对比结果见表4。由图8、9和表4看出:模型更新后的预测值更接近实际值,说明对串联 混 合 模 型 进 行 更 新,可 提 高 模 型 的 准确性。图8 金浸出率动态变化的预测结果图9 稳态金浸出率的预测结果表4 串联混合模型和更新串联混合模型的误差对比模型动态MA E动态RM S E稳态MA E稳态RM S E串联混合模型0.0 1 360.0 0 790.0 0 820.0 0 43更新串联混合模型0.0 0 850.0 0 610.0 0 530.0 0 264 结论针对金氰化浸出过程建立了一种具有高预测精度的串联型混合模型,并通过计算实际值与预测值之间的误差对模型进行更新,提高了模型的准确性。但由于实际生产过程具有一定的复杂性,所建立串联混合模型在实际应用过程中仍有一定的局限性,还有待进一步研究改进。参考文献:1 刘元清.湿法冶金浸出过程建模与优化控制研究D.沈阳:东北大学,2 0 1 8.2 李勇,彭伟,刘洪波,等.从贵州某金矿石中氰化浸出金试验研究J.湿法冶金,2 0 1 9,3 8(2):8 4-8 7.3 黄彪林,黄李金鸿,曾祥荣,等.用硫氰酸铵从金精矿生物氧化渣中浸出金J.湿法冶金,2 0 2 1,4 0(6):4 5 1-4 5 6.4 牟兴建.湿法冶金金氰化浸出过程建模与优化D.沈阳:东北大学,2 0 1 4.434第4 2卷第4期莫文水:金氰化浸出过程混合建模及仿真算法分析研究5 顾群.湿法冶金下游流程的分层优化及全流程优化仿真算法分析研究J.湿法冶金,2 0 2 3,4 2(2):2 0 5-2 1 4.6 胡广浩,毛志忠,何大阔.浸出过程浸出率预测与在线优化J.系统工程学报,2 0 1 1,2 6(4):5 2 4-5 3 0.7 胡广浩,毛志忠,杨菲.基于选择性集成算法的浸出率混合预测模型J.仪器仪表学报,2 0 1 3,3 4(5):1 0 4 9-1 0 5 3.8 白璐.湿法冶金浸出过程非线性预测控制方法的研究D.沈阳:东北大学,2 0 1 7.9 张俊,毛志忠,贾润达.金氰化浸出过程动力学反应速度估计J.东北大学学报(自然科学版),2 0 1 2,3 3(1 1):1 5 3 8-1 5 4 1.1 0 邱洋.氯酸钾与漂白粉协同浸出贵州难浸金矿的机制研究D.贵阳:贵州大学,2 0 2 2.1 1 周天驰.湿法冶金浸出过程机制模型J.有色冶金设计与研究,2 0 1 4,3 5(4):4 0-4 2.1 2 王帆.湿法冶金浸出过程混合建模与优化控制研究D.沈阳:东北大学,2 0 1 9.1 3 张俊,毛志忠,贾润达,等.湿法冶金浸出过程金浓度的混合预报模型J.仪器仪表学报,2 0 1 3,3 4(3):6 2 9-6 3 4.1 4 张俊.金氰化浸出过程建模与自适应实时优化方法的研究D.沈阳:东北大学,2 0 1 5.1 5 陈溥.基于核偏最小二乘法的湿法冶金萃取过程建模J.湿法冶金,2 0 2 1,4 0(3):2 6 0-2 6 6.1 6 刘茜,程靖,梁建勋.混合碰撞建模方法及其试验验证J.机械工程学报,2 0 2 2,5 8(1):1 1 6-1 2 3.1 7 陈兆鑫.城市污水处理过程混合建模与控制D.沈阳:沈阳工业大学,2 0 2 2.1 8 杨思雨,王正盛,李伟.求解大型离散不适定问题的预处理分数阶T i k h o n o v正则化方法J.南京航空航天大学学报:英文版,2 0 2 2,3 9(增刊1):1 0 6-1 1 2.1 9 陈玉红.R B F网络在时间序列预测中的应用研究D.哈尔滨:哈尔滨工程大学,2 0 0 9.2 0 陆晓琴,王喜.新型自适应R B F神经网络应用于P P P模式风险控制研究J.江苏商论,2 0 1 8(1):9 9-1 0 2.2 1 D A S GU P T A D,YU S H,N I N O F.R e c e n ta d v a n c e si na r t i f i c i a li mm u n es y s t e m s:m o d e l sa n da p p l i c a t i o n sJ.A p p l i e dS o f tC o m p u t i n gJ o u r n a l,2 0 1 0(2):1 5 7 4-1 5 8 7.H y b r i dM o d e l i n ga n dS i m u l a t i o nA l g o r i t h mo fG o l dC y a n i d eL e a c h i n gP r o c e s sMO W e n s h u i(G u a n g x iM o d e r nP o l y t e c h n i cC o l l e g e,H e c h i 5 4 7 0 0 0,C h i n a)A b s t r a c t:F i r s t l y,am u l t i-s t a g ed y n a m i cm e c h a n i s m m o d e lw a se s t a b l i s h e df o rag o l dc y a n i d el e a c h i n gp r o c e s s.O nt h i sb a s i s,t h el e a r n i n ga l g o r i t h mo fR B Fn e u r a ln e t w o r kd a t am o d e lb a s e do ni mm u n ep r i n c i p l ew a su s e dt oc a l c u l a t e t h ed y n a m i cr e a c t i o nr a t ew h i c hw a sd i f f i c u l t t om e a s u r e i nt h ea c t u a lp r o d u c t i o np r o c e s s,s oa st oe s t i m a t et h eu n k n o w np a r a m e t e r s i nt h ed y n a m i cm e c h a n i s m m o d e l,a n dt h e nc o n n e c t e dw i t ht h em a t e r i a l c o n s e r v a t i o ne q u a t i o n i ns e r i e s.T h e r e f o r e,as e r i e sm i x e dm o d e lw a se s t a b l i s h e dt op r e d i c t t h el e a c h i n gr a t e.S e c o n d l y,t h em o d e lw a sc o m p a r e dw i t ht h ep u r ed a t am o d e lb a s e do n t h e i mm u n ep r i n c i p l eo fR B Fn e u r a l n e t w o r k t ov e r i f y t h e e f f e c t i v e n e s so f t h em o d e l.F i n a l l y,t h em o d e lw a su p d a t e da c c o r d i n gt ot h ed e v i a t i o nb e t w e e nt h ea c t u a lv a l u ea n dt h ep r e d i c t e dv a l u e i nt h ea c t u a l p r o d u c t i o np r o c e s s,a n dt h ea c c u r a c yo f t h eu p d a t i n gm e c h a n i s m w a sv e r i f i e db ys i m u l a t i o na n a l y s i s.T h er e s u l t s s h o wt h a t t h e s e r i e sm i x i n gm o d e l g r e a t l y i m p r o v e s t h ep r e d i c t i o na c c u r a c yo f t h el e a c h i n gr a t e i nt h e l e a c h i n gp r o c e s s,a n dt h eu p d a t i n gm e c h a n i s mo f t h em o d e l i m p r o v e s t h ea c c u r a c ya n dg e n e r a l i z a t i o na b i l i t yo f t h em o d e l.K e yw o r d s:g o l d;l e a c h i n g;m o d e l i n g;s i m u l a t i o n;a l g o r i t h m534