一种稀土萃取分离过程组分含量区间控制方法技术

一种稀土萃取分离过程组分含量区间控制方法，所述方法针对稀土萃取分离过程各流量/各组分含量过程控制特点，建立稀土萃取过程回声状态网络模型；提出广义预测控制的稀土萃取分离多组分含量的区间控制方法，实现稀土萃取分离多组分含量的区间控制。传统方法采用萃取过程平衡状态下的软测量模型即静态模型，难以实现萃取过程组分含量在线预测以及难以建立精确的控制模型，从而影响稀土组分含量跟踪控制的效果。本发明专利技术控制方法，根据区间控制策略进行调整，优化计算，得到稀土萃取过程的准确控制量，使稀土萃取过程组分含量满足区间控制要求，保证了两端出口产品的质量。本发明专利技术适用于稀土萃取过程建模和优化控制。

【技术实现步骤摘要】

本专利技术涉及多组分稀土萃取分离过程组分含量区间控制，属稀土湿法冶金过程控制领域。
技术介绍
稀土是不可再生资源，是发展高新技术不可缺少的重要元素，它对经济和社会发展起重要的作用。稀土萃取过程中各级组分含量分布及两端出口产品纯度的检测与控制，一直是稀土湿法冶金过程研究领域的核心和难题。实现稀土萃取过程各级组分含量分布的优化控制，对稳定生产工艺、提升产品质量，降低原材料能耗，提高企业经济效益具有重要的现实意义。研究稀土萃取过程组分含量分布控制，首先要建立稀土萃取分离过程模型。针对稀土萃取过程的建模，可分为静态建模和动态建模两种。静态模型，如串级萃取平衡计算模型，不能反应萃取分离过程动态变化特性，难以作为过程控制模型。动态建模，如一种简化的双线性动态模型，该模型在局部范围预测精度较高，但当工况发生较大变化时，仍存在预测误差偏大问题。目前稀土分离主要采用溶剂萃取法，而稀土萃取分离一般具有多组分、产品纯度要求高等特点，萃取过程机理复杂、是一个多变量、强耦合、大滞后、强非线性过程。为了能实现稀土萃取分离过程组分含量分布控制，必须建立能体现出上述复杂特征的动态模型。在稀土萃取过程组分含量控制方法发展过程中，PID控制以及带参数自适应的改进PID算法，都在一定程度上取得令人满意的效果。但稀土萃取过程是一个复杂的多变量、大滞后、非线性动态系统，PID算法难以获得理想的效果。基于模糊控制、专家系统和鲁棒自适应控制等稀土萃取过程智能控制算法，取得了比PID更好的控制效果，但是模糊控制方法和专家系统控制方法主要是模拟专家实际经验，当萃取工业运行环境发生变化或受到随机干扰时，难以在线调整；鲁棒自适应控制器是从全局稳定性的角度而设计的且忽略了系统动态性能。此外，上述控制方法主要以萃取段或洗涤段监测点的组分含量设定值作为调控目标，但未涉及稀土萃取分离过程区间控制。
技术实现思路
本专利技术的目的是，对稀土萃取分离复杂的非线性运行过程建立有效回声状态神经网络模型描述，通过模型预测输出和期望输出组分含量之间的偏差设计性能指标函数，设计多变量预测控制器以及区间控制策略，实现稀土萃取分离多组分含量的区间控制。本专利技术的技术方案是：所述方法针对稀土萃取分离过程各流量/各组分含量过程控制特点，建立稀土萃取过程回声状态网络模型；提出广义预测控制的稀土萃取分离多组分含量的区间控制方法，实现稀土萃取分离多组分含量的区间控制。本专利技术基于串级萃取理论的稀土萃取过程回声状态网络模型：由于稀土元素间分离系数较小，稀土工业生产中普遍把若干萃取槽串联起来，使被萃物质与水相、有机相多次接触，从而实现稀土元素间的有效分离，同时得到两个或两个以上高纯度、高收率的产品。根据工业现场经验，料液流量大小通常由产品产量决定，在萃取过程中基本保持恒定。由此可将稀土萃取过程描述为如下非线性函数关系式，即：y1(k)=f1[y1(k-1),u1(k),u2(k)]+ζ1y2(k)=f2[y2(k-1),u1(k),u2(k)]+ζ2---(1)]]>式中，y1为萃取段监测级组分含量，y2为洗涤段监测级组分含量，u1为萃取剂流量，u2为洗涤剂流量。萃取段和洗涤段监测级组分含量满足：y1min≤y1(t)≤y1maxy2min≤y2(t)≤y2max---(2)]]>其中，y1min、y1max、y2min和y2max分别为两端监测点组分含量上下限。然后，利用神经网络思想，将稀土萃取过程建模，回声状态网络(EchoStateNetworks,简称ESNs)是一种特殊类型的神经网络，其基本思想是使用大规模随机连接的递归网络，取代经典神经网络中的中间层，从而简化网络的训练过程。式(1)等效成式(3)的模型，即稀土萃取过程回声状态网络模型为：x(k+1)=f(Wxx(k)+Winu(k+1)+Wbacky(k))y(k+1)=fout(WoutTx(k+1))---(3)]]>式中，Wx，Win，Wback，分别为表示内部状态变量、输入和输出到内部状态的连接权值矩阵；维数分别为N×K，N×N，N×L，Wout表示状态储备池输出对于输出层的连接权值矩阵，维数为L×N，f＝(f1,f2,…,fN)为内部神经元激活函数，这里为双曲正切激励函数。fout是输出函数，为线性函数。ESNs将状态储备池设计成包含数量几百个(常取100～1000之间)神经单元，储备池的神经元为随机的稀疏连接(通常保持1％～5％连接)，其矩阵的谱半径小于1。在网络训练过程中，连接到储备池的连接权矩阵Wx，Win，Wback随机产生，一经产生就固定不变。而连接到输出的连接权矩阵Wout需要通过训练得到，因为状态变量、输入和输出与输出之间是线性关系，所以通常这些连接权只需通过求解线性回归问题得到。本专利技术稀土萃取分离多组分含量的区间控制方法如下：本专利技术提出基于广义预测控制的稀土萃取多组分含量区间的控制方法来实现稀土萃取过程组分含量区间控制，其框图如图2所示。采用预测模型架构，参考轨迹yr(k)为监测级组分含量的设定区间值，通过建立稀土萃取过程回声状态网络模型，模型输出ym(k)与实际工业现场的输出y(k)得到误差值em(k)，未来时刻的预测输出ym(k+i)和当前时刻误差值进行反馈校正神经网络预测模型的各个参数，然后再利用反馈校正后的网络进行下一步预测，得到反馈校正后的预测值yp(k+i)，判断预测值是否满足给定区间范围，然后根据区间控制策略进行调整，优化计算，得到过程的控制量，进行控制稀土萃取过程。在预测控制中采用如下的二次型性能指标，在预测时域内求解优化问题，以模型输出和期望输出的组分含量偏差设计目标函数：minJ(k)=Σi=1PΣj=1nqj(i)[yrj(k+i)-yj(k+i)]2+Σi=1MΣj=1mrj(i)[Δuj(k+i-1)]2---(4)]]>式(4)中，yj(k+i)为组分含量预测模型第j个输出变量超前i步的最优预测值；△uj(k+j-1)为控制增量序列；yrj(k+i)为第j个输出变量未来参考轨迹；P，M分别表示预测长度、控制长度；qj(i)，rj(i)分别为误差加权系数和控制加权系数。m、n为输出，输入变量个数，对象为2输入2输出系统，m、n都为2。将式(4)用矩阵形式表示为：minJ=||Yr-Y||Q2+||ΔU||R2---(5)]]>式(5)中，Y为组分含量预测输出，Yr为系统实际运行轨迹，△U为系统控制增量，P为模型预测长度，M为模型控制长度，其中：Y＝Y0+G△U，控制矩阵G∈R(2*N)×(2*N)。初始值Y0＝[y01(k+1),y01(k+2),…,y01(k+P),y02(k+1),y02(k+2),…,y02(k+P)]T，Yr＝[yr1(k+1),yr1(k+2),…,yr1(k+P),yr2(k+1)，yr2(k+2),…,yr2(k+P)]T，Q为误差权系数矩阵，Q＝diag(Q1,Q2,…,Qm)，R为控制权系数矩阵，R＝diag(R1,R2,…,Rn)，一般为单位阵，其本文档来自技高网...

【技术保护点】
一种稀土萃取过程组分含量区间控制方法，其特征是，所述方法针对稀土萃取分离过程各流量/各组分含量过程控制特点，建立稀土萃取过程回声状态网络模型；提出广义预测控制的稀土萃取分离多组分含量的区间控制方法，实现稀土萃取分离多组分含量的区间控制。

【技术特征摘要】
1.一种稀土萃取过程组分含量区间控制方法，其特征是，所述方法针对稀土萃取分离过程各流量/各组分含量过程控制特点，建立稀土萃取过程回声状态网络模型；提出广义预测控制的稀土萃取分离多组分含量的区间控制方法，实现稀土萃取分离多组分含量的区间控制。2.根据权利要求1所述一种稀土萃取过程组分含量区间控制方法，其特征是，所述回声状态网络模型为：x(k+1)=f(Wxx(k)+Winu(k+1)+Wbacky(k))y(k+1)=fout(WoutTx(k+1))---(1)]]>其中，式中，Wx，Win，Wback，分别为表示内部状态变量、输入和输出到内部状态的连接权值矩阵；维数分别为N×K，N×N，N×L，Wout表示状态储备池输出对于输出层的连接权值矩阵，维数为L×N，f＝(f1,f2,…,fN)为内部神经元激活函数，这里为双曲正切激励函数；fout是输出函数，为线性函数；ESNs将状态储备池设计成包含数量100～1000个神经单元，储备池的神经元为随机的稀疏连接，通常保持1％～5％连接，其矩阵的谱半径小于1；在网络训练过程中，连接到储备池的连接权矩阵Wx，Win，Wback随机产生，一经产生就固定不变；而连接到输出的连接权矩阵Wout需要通过训练得到，因为状态变量、输入和输出与输出之间是线性关系，所以通常这些连接权只需通过求解线性回归问题得到。3.根据权利要求1所述一种稀土萃取过程组分含量区间控制方法，其特征是，所述稀土萃取分离多组分含量的区间控制方法，以模型输出和期望输出的组分含量偏差设计目标函数：minJ(k)=Σi=1PΣj=1nqj(i)[yrj(k+i)-yj(k+i)]2+Σi=1MΣj=1mrj(i)[Δuj(k+i-1)]2---(2)]]>式(2)中，yj(k+i)为组分含量预测模型第j个输出变量超前i步的最优预测值；Δuj(k+j-1)为控制增量序列；yrj(k+i)为第j个输出变量未来参考轨迹；P，M分别表示预测长度、控制长度；qj(i)，rj(i)分别为误差加权系数和控制加权系数。m、n为输出，输入变量个数，对象为2输入2输出系统，m、n都为2。将式(2)用矩阵形式表示为：minJ=||Yr-Y||Q2+||ΔU||R2---(3)]]>式(3)中，Y为组分含量预测输出，Yr为系统实际运行轨迹，ΔU为系...

【专利技术属性】
技术研发人员：杨辉，何丽娟，陆荣秀，朱建勇，
申请(专利权)人：华东交通大学，
类型：发明
国别省市：江西;36