[0046] 以裂解反应炉温度为实际对象,以燃料阀门的开度为输入,以裂解反应炉的温度为输出,来建立裂解反应炉温度的模型。
[0047] 本发明方法的步骤包括:
[0048] 步骤1、采集过程的实时运行数据,建立过程对象改进的神经网络模型,具体步骤如下:
[0049] 1-1.根据建立模型,综合主成分分析方法,得到给定子集p的最优解如下形式:
[0050]
[0051] 其中,S为数据的协方差矩阵,S2协方差矩阵的平方,SP中选定S中p个变量的p×p子阵。 为求矩阵的逆矩阵,tr()求括号中矩阵的迹。
[0052] 将其转换为最小化问题表示为J1,形式如下:
[0053] J1=1/f1
[0054] 1-2.将输入输出数据与预测输出 之间的映射关系,即模型结构表示为如下形式:
[0055]
[0056] 其中,x(k)=[y(k-1),…,y(k-n),u(k-1),…,u(k-m)]是阀门开度的观测输入,y(k-1),…,y(k-n)分别是k-1,…,k-n时刻裂解炉温度输出,u(k-1),…,u(k-m)分别是k-1,…,k-m时刻选取的干扰。n和m分别是输出和输入的最大阶次,f表示模型的映射关系。
[0057] 1-3.引入隐藏神经元,模型预测输出 表示为如下形式:
[0058]
[0059] 其 中 , ] 是 连 接 隐 藏 神 经 元 到 输 出 层 的 权 值 ,nh是隐藏节点的个数。φi(||x(k)||)是隐层中第i个神经元输出,形式如下:
[0060]
[0061] ||x(k)-ci||是x(k)与ci之间的欧几里德距离,||x(k)||为x(k)的模, 与σi∈R分别是高斯函数的中心向量和宽度,exp是指数函数。
[0062] 1-4.给出N1个训练数据样本,Y1=[y1(1),…,y1(N1)],U=[u(1),…,u(N1)],可以通过递归方法计算权重系数:
[0063]
[0064] 其中,Y1是训练数据,U是主成分分析方法所选取的干扰,0<μ<1是遗忘因子,P(k)、P(k-1)是k,k-1时刻的正定的协方差矩阵,P(0)=α2I,I是一个(n+m)×(n+m)单位矩5
阵,α是一个足够大的实数。设置为10,ω(0)=ε,ε是一个充分小的n+m的实矩阵,设置为10-3,K(k)是一个权重矩阵,fi(k)是第i个神经元k时刻训练输出,ωi(k-1)是第i个神经元k-1时刻训练权重,T是转置符号。
[0065] 1-5.建立的改进模型被训练,其建模精度通过使用训练和测试数据根据数学误差处理进行评估:
[0066]
[0067] 其中,J2是经数学处理误差目标,,k=1,…,N1, 是其预测输出。y2(k)是测试数据,k=1,…,N2, 是其预测输出。
[0068] 步骤2、基于改进的遗传算法的变量选择和神经网络建模,改进的遗传算法同时优化了两个目标J1、J2,然后设计了编码方法和变量选择,选择模型的结构和参数优化的各种算法以解决问题,其步骤是:
[0069] 2-1.初始化种群大小Np、最大进化代数G、算子概率Pc、Pm,首先对神经网络模型参数进行编码,为了简单起见,输入层中的n被设置为2,而根据现有知识将一个输入变量的m设置为1。隐藏层中的神经元数nh及其高斯函数参数ci,σi,i=1,…,nh被优化,1≤nh≤H,H是隐层的最大隐藏层数。设计了不同变量选择的编码和改进的神经网络,并导出第i条染色体Ci,形式如下所示:
[0070]
[0071] 其中1≤i≤Np,Np是群体大小。
[0072] 2-2.根据矩阵知识导出[1,nh]行中的元素,如下所示:
[0073]
[0074] σj=rwmax 1≤j≤nh
[0075] 其中r在[0.01,1]中随机产生的系数,umax和umin是最大和最小输入,ymax和ymin是最大和最小输出,wmax是高斯基函数的最大宽度,设置为umax、ymax的较大值。
[0076] 最后一行Ci代表3-8列被选择,他们是由特定的编码模式表示,其有效位为[3-8]。例如编码cH+1:
[0077] cH+1=[0 0 0 0 1 1 0 1 1 0]
[0078] 这意味着u3,u4,u6被选择,列c5,c6,c8是高斯函数的有效中心。一旦Ci获得,确定了改进的神经网络的结构和参数,然后由步骤1-4中方法使用训练数据获得权值ω。
[0079] 2-3.改进的遗传算法并选择算子
[0080] 在改进的遗传算法中使用排序算法,得到排序和拥挤距离。排第一的个体被视为精英,被选为父母。为了保持群体的多样性,具有相同值的J1和J2的个体被认为是一个个体。排第1的个体选择到父母群体中,直到超过群体规模。通过下降排序比较当前前面的拥挤距离,并且将具有较大拥挤距离的个体选择到父母群体中。如果大小仍然小于设定的群体规模,在改进算法的基础上,来选择J1和J2剩余种群数的一半。在选定的群体中通过优化的遗传方法,产生后代。
[0081] 2-4类染色体联会过程,产生遗传算法的算子,在Ci行和C′i行之间执行概率为pc的算子,见图1,其中,在Ci和C′i之间在[1,9]之间随机生成交叉位置。基函数的参数发生变化,后代中所选择的变量也发生了变化。但是隐藏节点的数量不能改变。
[0082] 2-5.步骤2-1中的元素以概率Pm突变。当实现变异算子时,根据步骤2-2生成元素,其中的元素执行逻辑非运算,即1到0和0到1得到新的结构,然后可以获得改进的神经网络模型和不同的变量。
[0083] 2-6.延长剪切算子
[0084] 如果隐藏神经元的数目小于2,则随机添加神经元新元素[1,H-2]之间的随机数,并根据步骤2-2计算新神经元的元素。在Ci中只有一个不为零的神经元将被剪切,并且隐藏的神经元的数量减少。
[0085] 2-7.循环重复优化搜素,依照步骤1-1到2-6,直到达到允许的最大进化代数结束优化搜索计算,得到改进后的遗传算法优化后的染色体,经解码后得到优化后的神经网络模型的参数。
[0086] 综上,本发明针对当前对于裂解反应炉温度的动态特性,使用神经网络具有良好的逼近速度,同时可以提高压力预测模型的精度,又可以简化模型结构。基于实际过程提出一种新型神经网络来改善模型精度和简化其结构。通过改进的主成分分析方法,建立在自然选择和自然遗传学基础上的迭代自适应随机全局优化搜索算法,能够解决许多传统优化方法不能解决的难题。通过选取合适的遗传算子,将改进的遗传算法用于优化神经网络模型,对裂解反应炉温度应用改进的主成分分析方法和神经网络建模精度,使用优化的进化算法来解决问题,通过和神经网络模型相结合,既能迅速逼近裂解反应炉实际温度,又保证了模型响应能力强结构简单的特点。