CN109086469B - 基于递归神经网络与偏好信息的铝电解建模与优化方法 - Google Patents

基于递归神经网络与偏好信息的铝电解建模与优化方法 Download PDF

Info

Publication number
CN109086469B
CN109086469B CN201810193122.5A CN201810193122A CN109086469B CN 109086469 B CN109086469 B CN 109086469B CN 201810193122 A CN201810193122 A CN 201810193122A CN 109086469 B CN109086469 B CN 109086469B
Authority
CN
China
Prior art keywords
individual
aluminum
aluminum electrolysis
function
production process
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810193122.5A
Other languages
English (en)
Other versions
CN109086469A (zh
Inventor
易军
白竣仁
吴凌
陈雪梅
周伟
陈实
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201810193122.5A priority Critical patent/CN109086469B/zh
Publication of CN109086469A publication Critical patent/CN109086469A/zh
Application granted granted Critical
Publication of CN109086469B publication Critical patent/CN109086469B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks

Abstract

本发明公开了一种基于递归神经网络与偏好信息的铝电解建模与优化方法，首先利用递归神经网络对铝电解生产过程进行建模，然后决策者设定期望目标值，再利用偏好多目标量子个体群算法对生产过程模型进行优化，得到各决策变量的一组最满足决策者期望的最优解以及该最优解对应的电流效率、槽电压、全氟化物排放量和吨铝能耗。利用差分进化算法中变异、交叉和选择操作，对决策变量进行偏好寻优，以此确定铝电解生产过程中工艺参数的最优值，可有效提高电流效率，降低槽电压，减少温室气体排放量和吨铝能耗，满足决策者偏好的同时，达到节能减排的目的。

Description

基于递归神经网络与偏好信息的铝电解建模与优化方法
技术领域
本发明属于最优控制领域，具体涉及一种基于递归神经网络与偏好信息的铝电解建模与优化方法。
背景技术
环保型铝电解生产过程长期以来都倍受重视，但非常具有挑战性。在电解铝工业中，最终目标是在电解槽平稳运行的基础上，提高电流效率、降低槽电压以及减少全氟化物、减少吨铝能耗的排放量。然而，铝电解槽参数较多，并且参数间呈现出非线性、强耦合性，给铝电解生产过程建模带来了较大难度，而递归神经网络具有很强的非线性映射能力，适用于解决非线性系统建模问题，为铝电解生产过程建模提供了新的思路。对于四个目标，同时实现则非常困难，因为目标相互之间存在冲突的现象，因此可引入决策者的偏好信息，设定期望目标，灵活调整不同目标之间的权重，利用偏好R-PMDE优化算法进行变量优化。R-PMDE是在DE算法的基础上，引入偏好R支配方法。DE是一种经典的进化算法，该算法简单、运算速度快、进化过程可直接用方程描述，因而被广泛应用于多个领域。
发明内容
本申请通过提出一种基于递归神经网络与偏好信息的铝电解建模与优化方法，以解决现有技术中铝电解生产过程中因无法获得最优工艺参数而导致的耗能巨大、效率低且严重污染环境的技术问题。
本发明的目的是这样实现的：
一种基于递归神经网络与偏好信息的铝电解建模与优化方法，包括如下步骤：
S1：选择对电流效率、槽电压以及全氟化物排放量有影响的控制参数构成决策变量X＝[x1,x2,···,xM]，M为所选控制参数的个数；
S2：选定铝电解工业现场，采集N组决策变量X1,X2,···,XN及其对应的电流效率y1,y2,···,yN，槽电压z1,z2,···,zN，以及全氟化物排放量作s1,s2,···,sN和吨铝能耗c1,c2,···,cN为数据样本，以每一组决策变量Xi作为输入，以对应的电流效率yi、槽电压zi以及全氟化物si和吨铝能耗ci作为输出，利用递归神经网络对样本进行训练、检验，建立四个铝电解槽生产过程模型；
S3：利用基于R支配的偏好多目标差分进化算法，根据决策者预先设定的期望值作为参考点，建立基于R支配的严格偏序关系，对步骤S2所得的四个生产过程模型进行优化，得到一组最满足决策者期望的决策变量Xbest及其对应的电流效率ybest、槽电压zbest以及全氟化物sbest和吨铝能耗cbest；
S4：按照步骤S3所得的最优决策变量Xbest中的控制参数来控制步骤S2中所选定的铝电解工业现场，使其达到节能减排降耗的目的。
优选地，步骤S1中，所述控制参数包括系列电流、下料次数、分子比、出铝量、铝水平、电解质水平、槽温。
优选地，步骤S2中，以电流效率作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
优选地，步骤S2中，以槽电压作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
优选地，步骤S2中，以全氟化物排放量作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
优选地，步骤S2中，以吨铝能耗作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
优选地，步骤S3包括以下步骤：
S31：根据R支配的偏好关系，评价每个个体的适应度，并根据优劣对个体最优值和全局最优值进行替换：
S32：更新种群内个体的基因信息，包括变异操作、交叉操作和选择操作。
优选地，步骤S31包括以下步骤：
S311：初始化系统参数，包括种群规模R，最大迭代次数T，随机生成n个个体x1,x2,···,xn，令外部存档集Q为空；
S312：决策者设定偏好目标参考点r(yp,zp,sp,cp)，即电流效率、槽电压、全氟化物排放量和吨铝能耗四个目标的期望值；
S313：对于每一个个体x，计算其适应度及其与参考点的距离：
Figure BDA0001592291060000031
Figure BDA0001592291060000041
其中，fj(x)是个体x在第j目标上的适应度值，ωj是第j个目标的权重，
Figure BDA0001592291060000042
是第j个目标值的上界，
Figure BDA0001592291060000043
是第j个目标值的下界；
S314：判断任意两个个体xi与xk之间的优劣关系，若xiPareto支配xk，即
Figure BDA0001592291060000044
则认为xi优于xk，反之亦成立；若两者之间无Pareto支配关系，则计算其偏好比较指标D(xi,xk,r)：
Figure BDA0001592291060000045
若D(xi,xk,r)＜-δ，则认为xi R支配xk，即
Figure BDA0001592291060000046
xi优于xk；若0≥D(xi,xk,r)≥-δ，则认为xi与xk相互之间非R支配，即，两者可视作等价，无优劣之分，其中δ∈[0,1]，是预先设定的阈值，
S315：确定个体的最优基因pbesti，在系统初始化时，个体最优基因设为该个体的初始基因xi；在下一次迭代后，基于S314提出的R支配关系，对个体的新基因xi与pbesti进行优劣比较，优秀者保存为pbesti，
S316：更新外部存档集Q，对种群中相互之间非R支配的个体加入存档集Q，删除被支配的个体；
S317：利用拥挤机制和禁忌算法在外部存档集Q中随机选择一个个体作为全局最优基因。
优选地，步骤S32包括以下步骤：
S321：对种群进行变异操作，对于每一个体xi，在种群中随机选择其他三个不同的个体xr1,xr2,xr3，将其中任意两个体形成的差分矢量通过比例因子F缩放后加到第三个个体上，以此产生变异个体，公式如下：
vi＝xr1+F·(xr2-xr3)，i≠r1≠r2≠r3
其中，r1,r2,r3为从集合{1,2,…,n}中随机选择的互不相同的整数，并且每进行一次变异，这些整数都会重新随机选取；
S322：目标个体xi与其变异个体vi进行交叉操作，生成试验个体ui；
S323：试验个体ui与目标个体xi进行选择操作，以确定哪一个个体进入下一代；
S324：判断当前全局最优解是否满足条件或者迭代次数是否达到最大迭代次数T，如果是，则输出当前全局最优解，否则，跳转至步骤S321进行重复计算，直到当前全局最优解满足条件或者迭代次数达到最大迭代次数T。
由于采用了上述技术方案，本发明首先利用递归神经网络对铝电解生产过程进行建模，然后决策者设定期望目标值，再利用偏好多目标量子个体群算法对生产过程模型进行优化，得到各决策变量的一组最满足决策者期望的最优解以及该最优解对应的电流效率、槽电压、全氟化物排放量和吨铝能耗。利用差分进化算法(DE)中变异、交叉和选择操作，对决策变量进行偏好寻优，以此确定铝电解生产过程中工艺参数的最优值，可有效提高电流效率，降低槽电压，减少温室气体排放量和吨铝能耗，满足决策者偏好的同时，达到节能减排的目的。
附图说明
图1为本发明的方法流程图；
图2为CF4排放量预测结果图；
图3为CF4排放量预测误差图
图4为电流效率预测结果图；
图5为电流效率预测误差图；
图6为槽电压排放量预测结果图；
图7为槽电压排放量预测误差图；
图8为吨铝能耗预测结果图；
图9为吨铝能耗预测误差图。
具体实施方式
如图1所示，一种基于递归神经网络与偏好信息的铝电解建模与优化方法，包括如下步骤：
S1：选择对电流效率、槽电压以及全氟化物排放量有影响的控制参数构成决策变量X＝[x1,x2,···,xM]，M为所选控制参数的个数。
本实施例是通过统计铝电解生产过程中对电流效率、槽电压以及全氟化物排放量和吨铝能耗有影响的原始变量，并从中确定对电流效率、槽电压以及全氟化物排放量和吨铝能耗影响大的参数作为决策变量X。
通过对实际工业生产过程中测量参数进行统计，得到对电流效率、槽电压以及全氟化物排放量和吨铝能耗最大的变量为：流x1、下料次数x2、分子比x3、出铝量x4、铝水平x5、电解质水平x6、槽温x7共7个变量。
S2：选定铝电解工业现场，采集N组决策变量X1,X2,···,XN及其对应的电流效率y1,y2,···,yN，槽电压z1,z2,···,zN，以及全氟化物排放量作s1,s2,···,sN和吨铝能耗c1,c2,···,cN为数据样本，以每一组决策变量Xi作为输入，以对应的电流效率yi、槽电压zi以及全氟化物si和吨铝能耗ci作为输出，利用递归神经网络对样本进行训练、检验，建立四个铝电解槽生产过程模型。为了满足建模需求，步骤S2中的递归神经网络包括输入层、隐藏层和输出层。
所述四个铝电解槽生产过程模型包括：
针对电流效率所构建的生产过程模型而言，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000；
针对槽电压所构建的生产过程模型而言，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000；
针对全氟化物所构建的生产过程模型而言，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
针对吨铝能耗所构建的生产过程模型而言，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
在本实施例中，采集重庆天泰铝业有限公司系列电解槽中的223#槽电解槽2013年全年生产数据以及2014年前40天数据，共计405组数据，其中，2013年全年生产数据作为建模训练样本，2014年的40组数据作为测试样本。数据样本如下表1所示。
表1数据样本
Figure BDA0001592291060000071
Figure BDA0001592291060000081
在递归神经网络设计中，由于存在递归信号，网路状态随时间的变化而变化，因此除了隐层节点数，学习速率也同样影响着神经网络模型的稳定性与准确性，是神经网络设计中的重难点。
隐层的节点数的设定由试凑法获得：
Figure BDA0001592291060000082
式中，p为隐层神经元节点数，n为输入层神经元数，m为输出层神经元数，k为1-10之间的常数。
最佳学习速率取值为：
Figure BDA0001592291060000083
Figure BDA0001592291060000084
本例中递归神经网络的设置参数如下表2所示。
表2递归神经网络设置参数
Figure BDA0001592291060000085
Figure BDA0001592291060000091
神经网络的训练过程中主要按照以下步骤进行：
设置Xk＝[xk1,xk2,···,xkM](k＝1,2,···,N)为输入矢量，N为训练样本个数，
Figure BDA0001592291060000092
为第g次迭代时输入层M与隐层I之间的权值矢量，WJP(g)为第g次迭代时隐层J与输出层P之间的权值矢量为Yk(g)＝[yk1(g),yk2(g),···,ykP(g)],(k＝1,2,…,N)为第g次迭代时网络的实际输出，dk＝[dk1,dk2,···,dkP],(k＝1,2,…,N)为期望输出；
步骤S2中建立铝电解生产过程模型具体包括如下步骤：
S21：初始化，迭代次数g初值设为0，WMI(0)、WJP(0)均为为(0,1)区间的随机值；S22：随机输入样本Xk；
S23：对输入样本Xk，前向计算递归神经网络每层神经元的输入信号和输出信号
S24：根据期望输出dk和实际输出Yk(g)，计算误差E(g)；
S25：判断误差E(g)是否满足要求，如不满足，则进入步骤S26，如满足，则进入步骤S29；
S26：判断迭代次数g+1是否大于最大迭代次数，如大于，则进入步骤S29，否则，进入步骤S27；
S27：对输入样本Xk反向计算每层神经元的局部梯度；
网络输出层节点误差为：e(k)＝d(k)-y(k),e(k)为网络期望输出，y(k)为网络实际输出。
通过计算输出层节点误差对各层的权值变化率为：
Figure BDA0001592291060000101
Figure BDA0001592291060000102
其中βij(0)＝0；i＝1,2,···,n1；j＝1,2,···,n0。
Figure BDA0001592291060000103
δi(0)＝0；i＝1,2,···,n1。
其中
Figure BDA0001592291060000104
分别表示隐含层第i个节点的输入及输出；n0、n1分别为输出层和隐含层节点数；
Figure BDA0001592291060000105
分别表示关联层、输出层、隐含层权值。
S28：网络权值修正计算公式为：
Figure BDA0001592291060000106
其中w(k)可为
Figure BDA0001592291060000107
式中w(k)可代表输出层、隐含层或输入层的权值，η为学习速率，令g＝g+1，跳转至步骤S23；
S29：判断是否完成所有的训练样本，如果是，则完成建模，否则，继续跳转至步骤S22。
通过上述循环过程，可得到递归神经网络预测效果如图2、3、4、5、6、7、8、9所示。优化模型的建立是铝电解生产过程优化的基础，模型精度直接影响优化结果。通过对图2、3、4、5、6、7、8、9分析可知，经递归神经网络训练，电流效率的最大预测误差为0.41％，槽电压的最大预测误差为0.08％，四氟化碳CF4排放量预测误差-1.20％，吨铝能耗预测误差为0.81％，模型预测精度高，满足建模要求。
S3：利用基于R支配的偏好多目标差分进化算法，即R-PMDE算法，根据决策者预先设定的期望值(参考点)，建立基于R支配的严格偏序关系，对步骤S2所得的四个生产过程模型进行优化，得到一组最满足决策者期望的决策变量Xbest及其对应的电流效率ybest、槽电压zbest以及全氟化物sbest和吨铝能耗cbest；
在铝电解生产过程模型基础上，利用R-PMDE算法在各决策变量范围内对其进行优化，各变量具体变化范围如表3所示。
表3各变量取值范围
Figure BDA0001592291060000111
优选地，步骤S3中的R-PMDE算法包括以下步骤：
S31：根据R支配的偏好关系，评价每个个体的适应度，并根据优劣对个体最优值和全局最优值进行替换。
进一步地，步骤S31包括以下步骤：
S311：初始化系统参数，包括种群规模R，最大迭代次数T，随机生成n个个体x1,x2,···,xn，令外部存档集Q为空；
S312：决策者设定偏好目标参考点r(yp,zp,sp,cp)，即电流效率、槽电压、全氟化物排放量和吨铝能耗四个目标的期望值；
S313：对于每一个个体x，计算其适应度及其与参考点的距离：
Figure BDA0001592291060000121
Figure BDA0001592291060000122
其中，fj(x)是个体x在第j目标上的适应度值，ωj是第j个目标的权重，
Figure BDA0001592291060000123
是第j个目标值的上界，
Figure BDA0001592291060000124
是第j个目标值的下界；
S314：判断任意两个个体xi与xk之间的优劣关系，若xiPareto支配xk，即
Figure BDA0001592291060000125
则认为xi优于xk，反之亦成立；若两者之间无Pareto支配关系，则计算其偏好比较指标D(xi,xk,r)：
Figure BDA0001592291060000126
若D(xi,xk,r)＜-δ，则认为xi R支配xk，即
Figure BDA0001592291060000127
xi优于xk；若0≥D(xi,xk,r)≥-δ，则认为xi与xk相互之间非R支配，其中δ∈[0,1]，是预先设定的阈值，
S315：确定个体的最优基因pbesti，在系统初始化时，个体最优基因设为该个体的初始基因xi；在下一次迭代后，基于S314提出的R支配关系，对个体的新基因xi与pbesti进行优劣比较，优秀者保存为pbesti，
S316：更新外部存档集Q，对种群中相互之间非R支配的个体加入存档集Q，删除被支配的个体；
S317：利用拥挤机制和禁忌算法在外部存档集Q中随机选择一个个体作为全局最优基因。
S32：更新种群内个体的基因信息，包括变异操作、交叉操作和选择操作。
进一步地，步骤S32包括以下步骤：
S321：对种群进行变异操作，对于每一个体xi，在种群中随机选择其他三个不同的个体xr1,xr2,xr3，将其中任意两个体形成的差分矢量通过比例因子F缩放后加到第三个个体上，以此产生变异个体，公式如下：
vi＝xr1+F·(xr2-xr3)，i≠r1≠r2≠r3
其中，r1,r2,r3为从集合{1,2,…,n}中随机选择的互不相同的整数，并且每进行一次变异，这些整数都会重新随机选取；
S322：目标个体xi与其变异个体vi进行交叉操作，生成试验个体ui；以二项式交叉方式为例，首先设置交叉概率常数CR，对于M维变量中的每一维变量j，若生成的[0,1]的随机数小于或等于CR，则进行交叉操作。这种交叉策略可概括为：
Figure BDA0001592291060000131
其中，randij为[0,1]均匀分布的随机数，决定了第i个试验个体的第j个元素由变异个体还是目标个体贡献。
S323：试验个体ui将与目标个体xi进行竞争，以确定哪一个个体进入下一代。采用贪婪的选择策略，以最小化优化为例，其选择操作如下：
Figure BDA0001592291060000132
如果试验个体的目标函数值小于或等于相应目标个体的目标函数值，则试验个体代替目标个体进入下一代。
S324：判断当前全局最优解是否满足条件或者迭代次数是否达到最大迭代次数T，如果是，则输出当前全局最优解，否则，跳转至步骤S321进行重复计算，直到当前全局最优解满足条件或者迭代次数达到最大迭代次数T。
通过上述步骤对铝电解生产过程进行优化可得100组最优的决策变量与对应的输出值，选取其中最合理的3组列于下表4中。
表4最佳生产参数
y<sub>1</sub> y<sub>2</sub> y<sub>3</sub> y<sub>4</sub> x<sub>1</sub> x<sub>2</sub> x<sub>3</sub> x<sub>4</sub> x<sub>5</sub> x<sub>6</sub> x<sub>7</sub>
99.24 3635 3.65 10835.15 1649 628 2.55 1210 16.5 14.5 942
98.13 3682 3.58 11527.21 1653 627 2.38 1200 17 15 924
95.37 3605 3.68 10478.52 1670 617 2.47 1090 17.5 15.5 935
利用最佳运行参数与2013年全年记录的平均值进行对比可知，电流效率提高了3.99％、槽电压降低了158mv，CF4排放量降低了0.38kg，吨铝能耗降低了1219.27KWh/t-Al。
S4：按照步骤S3所得的最优决策变量Xbest中的控制参数来控制步骤S2中所选定的铝电解工业现场，使其达到节能减排降耗的目的。
本申请的上述实施例中，通过提供一种基于递归神经网络与偏好信息的铝电解建模与优化方法，首先利用递归神经网络对铝电解生产过程进行建模，然后决策者设定期望目标值，再利用偏好多目标量子个体群算法对生产过程模型进行优化，得到各决策变量的一组最满足决策者期望的最优解以及该最优解对应的电流效率、槽电压、全氟化物排放量和吨铝能耗。利用差分进化算法DE中变异、交叉和选择操作，对决策变量进行偏好寻优，以此确定铝电解生产过程中工艺参数的最优值，可有效提高电流效率，降低槽电压，减少温室气体排放量和吨铝能耗，满足决策者偏好的同时，达到节能减排的目的。
最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。

Claims (6)

1.一种基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，包括如下步骤：
S1：选择对电流效率、槽电压以及全氟化物排放量有影响的控制参数构成决策变量X＝[x1,x2,···,xM]，M为所选控制参数的个数；
S2：选定铝电解工业现场，采集N组决策变量X1,X2,···,XN及其对应的电流效率y1,y2,···,yN，槽电压z1,z2,···,zN，以及全氟化物排放量s1,s2,···,sN和吨铝能耗c1,c2,···,cN为数据样本，以每一组决策变量Xi作为输入，分别以对应的电流效率yi、槽电压zi以及全氟化物si和吨铝能耗ci作为输出，利用递归神经网络对样本进行训练、检验，建立四个铝电解槽生产过程模型；
S3：利用基于R支配的偏好多目标差分进化算法，根据决策者预先设定的期望值作为参考点，建立基于R支配的严格偏序关系，对步骤S2所得的四个生产过程模型进行优化，得到一组最满足决策者期望的决策变量Xbest及其对应的电流效率ybest、槽电压zbest以及全氟化物sbest和吨铝能耗cbest；
步骤S3包括以下步骤：
S31：根据R支配的偏好关系，评价每个个体的适应度，并根据优劣对个体最优值和全局最优值进行替换；
S32：更新种群内个体的基因信息，包括变异操作、交叉操作和选择操作；
步骤S31包括以下步骤：
S311：初始化系统参数，包括种群规模R，最大迭代次数T，随机生成n个个体x1,x2,···,xn，令外部存档集Q为空；
S312：决策者设定偏好目标参考点r(yp,zp,sp,cp)，所述偏好目标参考点包括电流效率、槽电压、全氟化物排放量和吨铝能耗四个目标的期望值；
S313：对于每一个个体x，计算其适应度及其与参考点的距离：
Figure FDA0003877115380000021
ωj∈[0,1],
Figure FDA0003877115380000022
其中，fj(x)是个体x在第j目标上的适应度值，第ωj是第j个目标的权重，
Figure FDA0003877115380000023
是第j个目标值的上界，
Figure FDA0003877115380000024
是第j个目标值的下界；
S314：判断任意两个个体xi与xk之间的优劣关系，若xiPareto支配xk，即
Figure FDA0003877115380000026
则认为xi优于xk，反之亦成立；若两者之间无Pareto支配关系，则计算其偏好比较指标D(xi,xk,r)：
Figure FDA0003877115380000025
若D(xi,xk,r)＜-δ，则认为xi R支配xk，即
Figure FDA0003877115380000027
xi优于xk；若0≥D(xi,xk,r)≥-δ，则认为xi与xk相互之间非R支配，即，两者可视作等价，无优劣之分，其中δ∈[0,1]，是预先设定的阈值；
S315：确定个体的最优基因pbesti，在系统初始化时，个体最优基因设为该个体的初始基因xi；在下一次迭代后，基于S314提出的R支配关系，对个体的新基因xi与pbesti进行优劣比较，优秀者保存为pbesti；
S316：更新外部存档集Q，对种群中相互之间非R支配的个体加入存档集Q，删除被支配的个体；
S317：利用拥挤机制和禁忌算法在外部存档集Q中随机选择一个个体作为全局最优基因；
步骤S32包括以下步骤：
S321：对种群进行变异操作，对于每一个体xi，在种群中随机选择其他三个不同的个体xr1,xr2,xr3，将其中任意两个体形成的差分矢量通过比例因子F缩放后加到第三个个体上，以此产生变异个体，公式如下：
vi＝xr1+F·(xr2-xr3)，i≠r1≠r2≠r3
其中，r1,r2,r3为从集合{1,2,…,n}中随机选择的互不相同的整数，并且每进行一次变异，这些整数都会重新随机选取；
S322：目标个体xi与其变异个体vi进行交叉操作，生成试验个体ui；
S323：试验个体ui与目标个体xi进行选择操作，以确定哪一个个体进入下一代；
S324：判断当前全局最优解是否满足条件或者迭代次数是否达到最大迭代次数T，如果是，则输出当前全局最优解，否则，跳转至步骤S321进行重复计算，直到当前全局最优解满足条件或者迭代次数达到最大迭代次数T；
S4：按照步骤S3所得的最优决策变量Xbest中的控制参数来控制步骤S2中所选定的铝电解工业现场，使其达到节能减排降耗的目的。
2.根据权利要求1所述的基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，步骤S1中，所述控制参数包括系列电流、下料次数、分子比、出铝量、铝水平、电解质水平、槽温。
3.根据权利要求1所述的基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，步骤S2中，以电流效率作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
4.根据权利要求1所述的基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，步骤S2中，以槽电压作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
5.根据权利要求1所述的基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，步骤S2中，以全氟化物排放量作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Tansig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
6.根据权利要求1所述的基于递归神经网络与偏好信息的铝电解建模与优化方法，其特征在于，步骤S2中，以吨铝能耗作为输出，建立铝电解槽生产过程模型，其输入层采用10个神经元节点，隐藏层采用15个神经元节点，输出层采用1个神经元节点，输入层到隐藏层之间传递函数为Logsig函数，隐藏层到输出层之间的函数为Purelin函数，样本训练时的迭代次数为1000。
CN201810193122.5A 2018-03-09 2018-03-09 基于递归神经网络与偏好信息的铝电解建模与优化方法 Active CN109086469B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810193122.5A CN109086469B (zh) 2018-03-09 2018-03-09 基于递归神经网络与偏好信息的铝电解建模与优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810193122.5A CN109086469B (zh) 2018-03-09 2018-03-09 基于递归神经网络与偏好信息的铝电解建模与优化方法

Publications (2)

Publication Number Publication Date
CN109086469A CN109086469A (zh) 2018-12-25
CN109086469B true CN109086469B (zh) 2022-11-11

Family

ID=64839576

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810193122.5A Active CN109086469B (zh) 2018-03-09 2018-03-09 基于递归神经网络与偏好信息的铝电解建模与优化方法

Country Status (1)

Country Link
CN (1) CN109086469B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109935304A (zh) * 2019-03-15 2019-06-25 重庆科技学院 一种基于用户体验测试系统的自闭症拥抱机智能设计建模与决策参数优化方法
CN113379068B (zh) * 2021-06-29 2023-08-08 哈尔滨工业大学 基于结构化数据的深度学习架构搜索方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007036653A1 (de) * 2007-07-25 2009-02-05 Varta Microbattery Gmbh Elektroden und Lithium-Ionen-Zellen mit neuartigem Elektrodenbinder
CN101984441A (zh) * 2010-10-27 2011-03-09 哈尔滨工业大学 基于eda技术的电子系统多目标可靠性容差设计方法
CN103345559A (zh) * 2013-07-10 2013-10-09 重庆科技学院 铝电解过程电解槽工艺能耗的动态演化建模方法
CN105302973A (zh) * 2015-11-06 2016-02-03 重庆科技学院 基于moea/d算法的铝电解生产优化方法
CN105302976A (zh) * 2015-11-06 2016-02-03 重庆科技学院 基于spea2算法的铝电解生产优化方法
CN105404926A (zh) * 2015-11-06 2016-03-16 重庆科技学院 基于bp神经网络与mbfo算法的铝电解生产工艺优化方法
CN105426960A (zh) * 2015-11-06 2016-03-23 重庆科技学院 基于bp神经网络与mbfo算法的铝电解节能减排控制方法
CN105809297A (zh) * 2016-05-18 2016-07-27 西南石油大学 一种基于多目标差分进化算法的火电厂环境经济调度方法
CN106403950A (zh) * 2016-06-24 2017-02-15 电子科技大学 一种基于参考信息的码垛机器人轨迹规划方法
CN107480813A (zh) * 2017-07-27 2017-12-15 河海大学 基于多目标混沌遗传算法的流域水资源优化配置方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7634360B2 (en) * 2003-09-23 2009-12-15 Prediction Sciences, LL Cellular fibronectin as a diagnostic marker in stroke and methods of use thereof
US20080104001A1 (en) * 2006-10-27 2008-05-01 Kipp James E Algorithm for estimation of binding equlibria in inclusion complexation, host compounds identified thereby and compositions of host compound and pharmaceutical

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007036653A1 (de) * 2007-07-25 2009-02-05 Varta Microbattery Gmbh Elektroden und Lithium-Ionen-Zellen mit neuartigem Elektrodenbinder
CN101984441A (zh) * 2010-10-27 2011-03-09 哈尔滨工业大学 基于eda技术的电子系统多目标可靠性容差设计方法
CN103345559A (zh) * 2013-07-10 2013-10-09 重庆科技学院 铝电解过程电解槽工艺能耗的动态演化建模方法
CN105302973A (zh) * 2015-11-06 2016-02-03 重庆科技学院 基于moea/d算法的铝电解生产优化方法
CN105302976A (zh) * 2015-11-06 2016-02-03 重庆科技学院 基于spea2算法的铝电解生产优化方法
CN105404926A (zh) * 2015-11-06 2016-03-16 重庆科技学院 基于bp神经网络与mbfo算法的铝电解生产工艺优化方法
CN105426960A (zh) * 2015-11-06 2016-03-23 重庆科技学院 基于bp神经网络与mbfo算法的铝电解节能减排控制方法
CN105809297A (zh) * 2016-05-18 2016-07-27 西南石油大学 一种基于多目标差分进化算法的火电厂环境经济调度方法
CN106403950A (zh) * 2016-06-24 2017-02-15 电子科技大学 一种基于参考信息的码垛机器人轨迹规划方法
CN107480813A (zh) * 2017-07-27 2017-12-15 河海大学 基于多目标混沌遗传算法的流域水资源优化配置方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
"Analysis of asymmetric neural networks with time delay and its applications in associative memories";Jin Liu;《2017 International Workshop on Complex Systems and Networks (IWCSN)》;20171208;全文 *
"Anode effect prediction of aluminum electrolysis using GRNN";Kaibo Zhou;《2015 Chinese Automation Congress (CAC)》;20151127;全文 *
"Operating Parameters Optimization for the Aluminum Electrolysis Process Using an Improved Quantum-Behaved Particle Swarm Algorithm";Jun Yi;《IEEE Transactions on Industrial Informatics》;20171207;全文 *
"基于决策偏好的进化多目标水库防洪调度算法研究";魏艳夕;《中国优秀硕士学位论文全文数据库》;20170216;全文 *
"基于拥挤距离排序的铝电解多目标优化";易军;《仪器仪表学报》;20151115;第36卷(第11期);全文 *
"基于量子优化的铝电解多目标控制系统理论研究";李界家;《轻金属》;20160620;全文 *

Also Published As

Publication number Publication date
CN109086469A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
CN109085752B (zh) 基于角度支配关系的铝电解偏好多目标优化算法
CN108445756B (zh) 基于ar支配关系的铝电解节能减排智能控制方法
CN105302973A (zh) 基于moea/d算法的铝电解生产优化方法
CN105447567B (zh) 基于bp神经网络与mpso算法的铝电解节能减排控制方法
CN105404926B (zh) 基于bp神经网络与mbfo算法的铝电解生产工艺优化方法
CN113282122B (zh) 一种商用建筑能耗预测优化方法及系统
CN108846526A (zh) 一种二氧化碳排放量预测方法
CN107992645B (zh) 基于混沌-烟花混合算法的污水处理过程软测量建模方法
CN109086469B (zh) 基于递归神经网络与偏好信息的铝电解建模与优化方法
CN105321000A (zh) 基于bp神经网络与mobfoa算法的铝电解工艺参数优化方法
CN105404142B (zh) 基于bp神经网络与mbfo算法的铝电解多参数控制方法
CN111832817A (zh) 基于mcp罚函数的小世界回声状态网络时间序列预测方法
CN110110447B (zh) 一种混合蛙跳反馈极限学习机带钢厚度预测方法
CN101893852B (zh) 一种pH中和过程的多目标建模方法
CN105426959B (zh) 基于bp神经网络与自适应mbfo算法的铝电解节能减排方法
Yüzgeç Performance comparison of differential evolution techniques on optimization of feeding profile for an industrial scale baker’s yeast fermentation process
CN108363303B (zh) 基于ar偏好信息的差分进化铝电解多目标优化方法
CN105302976A (zh) 基于spea2算法的铝电解生产优化方法
CN105334824A (zh) 基于nsga-ⅱ算法的铝电解生产优化方法
CN109100995B (zh) 基于决策者偏好信息的铝电解节能减排优化方法
CN115963731B (zh) 一种基于改进遗传算法的指挥控制体系网络结构优化方法
CN108984813B (zh) 基于递归神经网络与角度偏好的铝电解建模与优化方法
CN113486952B (zh) 一种基因调控网络的多因素模型优化方法
CN105426960B (zh) 基于bp神经网络与mbfo算法的铝电解节能减排控制方法
Bose et al. A novel interpolation based missing value estimation method to predict missing values in microarray gene expression data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant