最优BP多模丙烯聚合生产过程最优软测量仪表和方法

摘要

本发明公开了一种最优BP多模丙烯聚合生产过程最优软测量仪表，包括丙烯聚合生产过程、现场智能仪表、控制站、存放数据的DCS数据库、最优BP多模优化最优软测量仪以及熔融指数软测量值显示仪。现场智能仪表及控制站与丙烯聚合生产过程相连，与DCS数据库相连；最优软测量仪与DCS数据库及软测量值显示仪相连。所述的最优BP多模优化最优软测量仪包括模型更新模块、数据预处理模块、PCA主成分分析模块、神经网络模型模块、最优BP多模优化模块、模型更新模块。以及提供了一种用软测量仪表实现的软测量方法。本发明实现在线测量、在线参数优化、软测量速度快、模型自动更新、抗干扰能力强、精度高。

主权项

一种最优BP多模丙烯聚合生产过程最优软测量仪表，包括丙烯聚合生产过程、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的DCS数据库、最优BP多模优化最优软测量仪以及熔融指数软测量显示仪，所述现场智能仪表、控制站与丙烯聚合生产过程连接，所述现场智能仪表、控制站与DCS数据库连接，所述DCS数据库与最优BP多模优化最优软测量仪的输入端连接，所述最优BP多模优化最优软测量仪的输出端与熔融指数软测量显示仪连接，其特征在于：所述最优BP多模优化最优软测量仪包括：(1)、数据预处理模块，用于将从DCS数据库输入的模型输入变量进行预处理，对输入变量中心化，即减去变量的平均值；再进行归一化处理，即除以变量值的变化区间；(2)、PCA主成分分析模块，用于将输入变量预白化处理及变量去相关，通过对输入变量施加一个线性变换实现，即主成分由C＝xU得到，其中x为输入变量，C为主成分得分矩阵，U为载荷矩阵；若对原始数据进行重构，由x＝CU^T计算，其中上标T表示矩阵的转置；当选取的主成分数目小于输入变量的变量个数时，x＝CU^T+E，其中E为残差矩阵；(3)、神经网络模型模块，用于采用BP神经网络、通过误差函数最小化来完成输入到输出的一种高度非线性映射，映射中保持拓扑不变性；需要建立若干子神经网络，第一个子BP网络的训练目标是预报结果与实际结果差距J₁最小；$J_1=\frac{1}{N}\underset{j=1}{\overset{N}{\Sigma }}{\left(F_1\right(x_j)-d(x_j\left)\right)}^2---\left(1\right)$N为样本数目，x为输入变量，j为样本点序号，F₁(·)为子网络预报结果，d(·)为实际结果；从第二个子网络开始，训练目标变为使得网络的预报误差尽可能小，同时网络的预报结果与之前的网络预报结果又尽可能大的差异，目标函数如下：$J_i=\frac{1}{N}\underset{j=1}{\overset{N}{\Sigma }}{\left(F_i\right(x_j)-d(x_j\left)\right)}^2-\frac{\lambda }{N}\underset{j=1}{\overset{N}{\Sigma }}{\left(F_i\right(x_j)-F(x_j\left)\right)}^2---\left(2\right)$J_i为前i个子网络的训练目标，F_i(·)为第i个网络的预报结果；d(·)为实际结果；F(·)为前i‑1个子网络的综合结果；λ为调节参数，N为样本数目；训练的终止条件为将得到的新的子网络加入多模神经网络后，网络群的预报误差不再减小；(4)、最优BP多模优化模块，采用一种智能连续空间蚁群算法对每个BP网络进行训练和最优化，具体步骤为：(a)算法初始化，根据待优化的BP网络构造出初始的解集S＝(s₁,s₂,···,s_n)，n为初始解的个数，s_n为第n个初始解，确定蚁群的大小M，设置蚁群寻优算法迭代次数的阈值MaxGen并初始化蚁群寻优的迭代次数序号gen＝0；(b)计算出解集S对应的适应度值G_i(i＝1,2,···,n)，适应度值越大代表解越好；再根据下式确定解集中每个解被取到作为蚂蚁寻优的初始解的概率P_i(i＝1,2,···,n)$P_a\left(k\right)=\frac{G_a}{{\Sigma }_{a=1}^n{G}_i},(a=1,2,...,n)---\left(3\right)$n为初始解的个数，s_n为第n个初始解，k为迭代次数；初始化执行寻优算法的蚂蚁编号a＝0；(c)蚂蚁a选取S中的一个解作为寻优的初始解，选取规则是根据P来做轮盘选；(d)蚂蚁a在选取的初始解的基础上进行寻优，找到更好的解s_a'；(e)如果a<M，则a＝a+1，返回步骤c；否则继续向下执行步骤f；(f)如果gen<MaxGen，则gen＝gen+1，使用步骤d中所有蚂蚁得到的更好的解取代S中的对应解，返回步骤b；否则向下执行步骤g；(g)计算出解集S对应的适应度值G_a(a＝1,2,···,n)，选取适应度值最大的解作为算法的最优解，结束算法并返回；每一只蚂蚁在它选定的初始解的基础上寻优时会循环固定的次数，如果本次循环得到了更好的解，则在下次循环中会基于该解并保持搜索方向不变；否则在下次循环中仍基于原来的解但会调整搜索方向；同时随着整个蚁群寻优代数的增加，蚂蚁搜索的步长会智能地减小，以适合整个蚁群寻优的收敛：del_k＝Random·k^r  (4)式中，del_k为蚂蚁第k代迭代的初始步长，k为迭代代数，Random为随机向量，r为负常实数；对于解集S中长期不被蚂蚁选作寻优初始解的那些解，会采用遗传算法中的变异和交叉策略进行处理，从而提高算法的全局寻优性能；然后，对步骤(3)中的每个子网络按如下公式，进行构建$O\left(x\right)=\frac{1}{I}\underset{j=1}{\overset{I}{\Sigma }}{F}_j\left(x\right)---\left(5\right)$I为总的子网络数目，x为输入变量，O(·)为模型输出，F_j(·)为第j个子网络的输出；即最终多模BP神经网络的预报结果为各个子网络预报结果的平均值；所述最优BP多模优化最优软测量仪还包括：模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新神经网络模型；所述的最优软测量仪表，训练子BP神经网络，然后将其的构建起来形成神经网络群；由于子网络的选取标准是预报误差小、与其他的子网络差异大，所以这些预报效果好、又各不相同的子神经网络的综合预报效果能够具有更好的预报精度和稳定性；同时，PCA主成分分析模块中，PCA方法实现输入变量的预白化处理，能够简化神经网络模型的输入变量，进而提高模型的性能。