基于QPSO‑DMPC的反应再生系统优化控制方法与流程

本发明属于石油化工技术领域，涉及一种基于qpso-dmpc的反应再生系统优化控制方法。

背景技术：

石油化工工业在我国国民经济中占有举足轻重的地位，承担着为我国提供各种能源的重担。常规的催化裂化装置由三个部分组成，包含反应再生系统、分馏系统以及吸收稳定系统。作为催化裂化的核心部分，反应再生系统(reactionregenerationsystem，rrs)将原油经过加工，生成各种各样的轻质油产品。但现有的反应再生系统为非线性反应再生系统，存在控制精度低，且硬件负担大的问题。

技术实现要素：

鉴于上述问题，本发明的目的是提供一种基于qpso-dmpc的反应再生系统优化控制方法，以解决现有的非线性反应再生系统控制精度低且硬件负担大的问题。

本发明提供的基于qpso-dmpc的反应再生系统优化控制方法，包括：

s1：将反应再生系统的传递函数模型转化为阶跃响应模型；

s2：建立dmpc模型，dmpc模型包括开环预测模块、稳态目标计算模块和动态矩阵控制模块；

s3：利用qpso算法，在不放松约束条件的前提下，对经济优化函数进行求解；其中，约束条件包括操作变量的硬约束和软约束，被控变量的硬约束和软约束，外部目标的约束；

s4：根据qpso算法对经济优化函数求得的解获得反应再生系统的输出设定值，并与实际输出的偏差作为目标误差函数，利用qpso算法对该目标误差函数求解，获得操作变量的最佳变化量。

利用上述根据本发明提供的基于qpso-dmpc的反应再生系统优化控制方法，能够赋予粒子更高的随机性，最优粒子可以表达更大范围内的最优值，该方法不仅减小了rrs硬件负担，还能获取更优的操作变量参数，真正达到rrs的自适应最优控制。

附图说明

图1为根据本发明的基于qpso-dmpc的反应再生系统优化控制方法的流程图；

图2为根据本发明的dmpc对rrs输出的跟踪效果图；

图3为根据本发明的dmpc对rrs输入的跟踪效果；

图4为根据本发明的qpso-dmpc对rrs的输出的跟踪结果图；

图5根据本发明的qpso-dmpc对rrs的输出的跟踪结果图。

具体实施方式

在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐述了许多具体细节。然而，很明显，也可以在没有这些具体细节的情况下实现这些实施例。在其它例子中，为了便于描述一个或多个实施例，公知的结构和设备以方框图的形式示出。

名词解释

dmpc：thedouble-layerdmodelpredictivecontrol，双层模型预测控制。

qpso：quantum-behavedparticleswarmoptimizationalgorithm，量子粒子群优化算法。

图1示出了根据本发明的基于qpso-dmpc的反应再生系统优化控制方法的流程。

如图1所示，本发明提供的基于qpso-dmpc的反应再生系统优化控制方法，包括如下步骤：

s1：将rrs的传递函数模型转化为阶跃响应模型。

转换后的rrs的阶跃响应模型如下：

式(1)中，δu为操作变量的变化量，k为时间，n为模型长度，为rrs操作变量的阶跃响应系数矩阵，为rrs干扰变量的阶跃响应系数矩阵,对满足

s2：建立dmpc模型，dmpc模型包括开环预测模块、稳态目标计算模块和动态矩阵控制模块。

建立开环预测模块的过程，包括如下步骤：

s211：当δu(k+i-1)＝0、δv(k+i-1)＝0(1≤i≤p)时，设为对y(k+p|k)的预测值，其中，δν为干扰变量的变化量，p为预测时域，则有：

s212：考虑反馈校正，假设vss(k)＝vss(k-1)+δv(k)为已知，从k时刻开始，反应再生系统的操作变量不再变化时，基于式(2)得到反应再生系统的开环预测为y^ol(k+i|k)，当检测到δu(k-1)时求解得到反应再生系统的开环预测：

s213：基于式(3)与rrs的实际输出，得到误差

s214：对误差进行一阶指数平滑处理，得到：

s215：以平滑处理后的误差为基准，对rrs的输出进行反馈校正，且反馈校正在未来所有时间点都是恒定的，记为k时刻的开环动态预测值得到：

s216：结合式(4)，得到开放稳态预测：

建立稳态目标计算模块的过程，包括如下步骤：

s221：提取所有反应再生系统的操作变量和被控变量的硬约束条件与软约束条件，并合并表达为关于稳态操作变量的变化量δuss(k)的形式：

其中，为操作变量的上限，为操作变量的理想值的集合，为稳态增益矩阵，为稳态被控变量的变化量，为被控变量的理想值的集合，k为迭代次数，t为时间。

更为具体地，稳态mv的硬约束为：

在mpc控制过程中,存在mv变化速率约束0≤j≤m-1，其中，m为控制时域，则增加的稳态mv的硬约束为：

对δus(k)进行限制，则增加的稳态mv的硬约束为：

稳态cv的硬约束为：

稳态cv的软约束为：

在实际过程中，总是满足yo,h≤yo；另外，对δyss(k)进行限制，则增加的稳态cv的硬约束为

cv的新稳态值仅决定于δuss(k)的大小，而与mv动态变化路径无关，稳态预测模型为：

其中，为稳态增益矩阵；为开环稳态预测。

所有条件合并表达为关于稳态操作变量的变化量δuss(k)的形式：

s222：建立经济优化函数：

s223：放松约束条件，采用二次规划方法对式(5)进行求解，获得稳态操作变量的变化量δuss(k)。

建立动态矩阵控制模块的过程，包括如下步骤：

s231：取预测时域为p，控制时域为m，在每个时刻k，可得到：

s232：当p大于n时，y^ol(k+j|k)＝y^ol(k+n|k)，j＞n，该预测值包含预测误差的反馈校正及干扰的影响，得到：

其中，d为动态控制矩阵；

s233：在动态矩阵中，根据qpso算法对经济优化函数求得的解获得反应再生系统的输出设定值，并与实际输出的偏差作为目标误差函数，选择最小化的目标误差函数如下：

为了让预测输出尽可能地接近实际输出，以式(5)的解求得rrs的输出设定值和实际输出的误差为目标误差函数。

s234：对最小化的目标函数(6)求解，获得稳态操作变量的最佳变化量。

采用matlab7.0为仿真平台，以rrs为对象，进行dmpc算法的研究，仿真过程中，采样周期为4分钟，权重向量b＝(122211)，jmin＝-3，建模时域n＝30，作变量下限ui为0，预测控制的操作变量上限为600，被控变量下限yi为0，被控变量上限为800，稳态操作变量变化值δus(k)为100，操作变量变化值为50。

各操作变量代表名称如表1所示：

表1各操作变量代表名称

各被控变量代表名称如表2所示：

表2各被控变量代表名称

通过实验仿真，dmpc对输出的跟踪效果以及对输入的跟踪效果如图2和图3所示。

从图2和图3可以看出，在考虑各变量优先级顺序的条件下，通过放松约束条件对最佳的操作变量变化量进行求取，仿真结果表明，dmpc对rrs的输入和输出有很好的跟踪效果。然而，放松约束条件不仅对硬件设备提出了更高的要求，而且所求的最优解是通过放松约束条件后求取的最优解，并不是真正意义上的最优解。群体智能算法在不放松约束条件下，对最优化问题的求解比传统的二次规划或线性规划方法有天然的优势，因此，本发明将qpso算法引入到dmpc中。

s3：利用qpso算法，在不放松约束条件的前提下，对经济优化函数进行求解。

其中，约束条件包括操作变量的硬约束和软约束，被控变量的硬约束和软约束，外部目标的约束。

利用qpso算法，在不放松约束条件的前提下，对经济优化函数进行求解的过程如下：

s31：初始化系统参数，包括种群规模n、最大迭代次数t、随机生成n个粒子x1,x2，…,xn、粒子维数m、压缩-扩张因子α和令外部存档集q，q为空；

s32：评价每个粒子的适应度，并根据优劣对个体最优值和全局最优值进行替换；

s33：将每个粒子的当前适应度pi和个体最优适应度进行比较，如果当前适应度pi支配个体最优适应度则将当前适应度pi代替个体最优适应度否则，保留原有的个体最优适应度

s34：更新外部存档集q，将种群中所有的非支配集加入外部存档集q，并删除被支配的粒子；

s35：利用拥挤机制和禁忌算法在外部存档集q中随机选择一个粒子作为全局最优值；

s36：更新作为全局最优值的粒子的位置xij(t)，更新公式为：

eij(t)＝βpij(t)+(1-β)pgj(t)

其中，xij(t+1)表示作为全局最优值的粒子更新后的位置，pij为第i个粒子第m维的当前最优位置，pgj为全局最优位置，uij与β分别为0到1之间的随机数，α为扩张-收缩因子，会影响单个粒子的收敛性，从1至0.5随迭代次数自适应变化；

s37：判断当前全局最优解是否满足条件或者迭代次数是否达到最大迭代次数t，如果是，则输出当前全局最优解，否则，跳转至步骤s32进行重复计算，直到当前全局最优解满足条件或者迭代次数达到最大迭代次数t为止；

s38：不放松约束条件，利用qpso对建立的经济优化函数j进行求解，求得单目标下的稳态操作变量的变化量δuss(k)。

s4：根据qpso算法对经济优化函数求得的解获得反应再生系统的输出设定值，并与实际输出的偏差作为目标误差函数，利用qpso算法对该目标误差函数求解，获得操作变量的最佳变化量。

利用qpso算法对目标误差函数求解，获得操作变量的最佳变化量的过程，包括如下步骤：

s41：取预测时域为p，控制时域为m，在每个时刻k，可得到：

s42：当p大于n时，y^ol(k+j|k)＝y^ol(k+n|k)，j＞n，该预测值包含预测误差的反馈校正及干扰的影响，得到：

s43：在动态矩阵中，根据qpso算法对经济优化函数求得的解获得反应再生系统的输出设定值，并与实际输出的偏差作为目标误差函数，选择最小化的目标误差函数j(k)如下：

在动态矩阵中，根据qpso算法对经济优化函数求得的解获得反应再生系统的输出设定值的公式为：

式(5)中，yss(k)为反应再生系统的输出设定值，δuss(k)为qpso算法对经济优化函数求得的解，为稳态增益矩阵，由系统的稳态模型得到；为开环稳态预测，由辨识好的传递函数模型得到。

s44：利用qpso算法对最小化的目标函数求解，获得操作变量的最佳变化量。

步骤s44的操作过程请参照步骤s31-s37。

采用matlab7.0为仿真平台，以rrs为对象，进行各算法的研究，仿真过程中，采样周期为4分钟，建模时域n＝600，操作变量下限ui为-0.5，预测控制的操作变量上限为0.5，被控变量下限yi为-0.5，被控变量上限为0.5，稳态操作变量变化值δus(k)为0.1，操作变量变化值为0.1，b1＝[0.122]，a1＝[1020200]，j1max＝－3，j2max＝－4。各变量所代表意义如表1和表2所示，算法的参数取值如表3所示：

表3各算法参数取值表

qpso-dmpc对rrs输出的跟踪效果以及对输入的跟踪效果的如图4和图5所示。

从图4和图5中可以看出，通过设置rrs的经济优化函数，并采用qpso对rrs的经济优化问题求解，在保证经济效益的基础上，进一步对rrs的过程进行稳态控制，即采用qpso对dmpc的动态矩阵控制阶段进行求解，仿真结果表明，qpso-dmpc能对rrs的被控变量和操作变量进行跟踪，表明了qpso-dmpc算法在rrs中的有效性。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。