CN107942682B - 非高斯系统的动态经济性能优化与控制一体化设计方法 - Google Patents
非高斯系统的动态经济性能优化与控制一体化设计方法 Download PDFInfo
- Publication number
- CN107942682B CN107942682B CN201711400399.2A CN201711400399A CN107942682B CN 107942682 B CN107942682 B CN 107942682B CN 201711400399 A CN201711400399 A CN 201711400399A CN 107942682 B CN107942682 B CN 107942682B
- Authority
- CN
- China
- Prior art keywords
- economic
- vector
- control
- random
- layer
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Feedback Control In General (AREA)
Abstract
非高斯系统的动态经济性能优化与控制一体化设计方法,属于随机系统优化控制领域,包括:步骤一、非高斯扰动下的非线性随机系统模型;步骤二、建立两层动态优化控制框架结构;步骤三、基于统计信息的性能指标选取;步骤四、动态经济性能优化与模型预测控制器设计。本发明与传统的两层RTO‑MPC在上层采用稳态模型相比,本发明充分考虑了过程的动态特性,强调所得最优设定轨迹的可达性,更具有实际工程意义。
Description
技术领域
本发明具体涉及了一种在非高斯随机噪声影响下系统的经济性能优化与控制一体化设计方法,属于随机系统优化控制领域。
背景技术
在实际系统中,随机干扰广泛存在,尤其在工业过程当中,例如化工过程中反应堆的温度浓度的随机因素。因此,随机系统的优化控制一直是系统控制领域研究的热点问题。基于系统变量高斯性的假设,通过优化和控制均值、方差取得了大量的研究成果,并已广泛应用到工程实际中。然而,实际工业过程变量往往是非高斯的,这使得对于非高斯分布控制系统的研究显得更为有必要。受造纸过程中经典案例的启发,王宏教授首次提出了随机分布控制理论,为非高斯随机系统的控制问题研究开辟了一个全新的思路。
先进控制方法的实施需要增加企业投资成本,能否在项目实施前给出准确的投运成本和经济收益分析,是企业关心的一个重要问题。在先进控制的运行中,随着过程操作条件和设备等变化,先进控制系统的性能会随时间的推移而退化,无法达到预期的经济效益。因此,保持先进控制系统投运以后的高效、优质运行,及时对先进控制系统进行性能评估和维护是非常有必要的。先进控制策略对于经济性的提高,是通过减小过程随机不确定性,优化操作点来实现的。其对于企业经济效益的提高主要表现在:例如化工过程反应物转化为产品的过程中增加产品产量、增加高附加值产品生产配额、降低能耗、减少再加工消耗。因此,为了保证企业的经济效益,应该使所设计的先进控制系统始终在最优操作点运行。在这样的基础之下,研究非高斯系统的动态经济性能优化与模型预测控制一体化设计方法就显得很有理论意义和实际应用价值。
目前,针对非高斯随机系统控制和经济性能优化研究已经取得了比较多的成果。但就当下比较热门的研究成果来说,仍然存在着以下的一些问题:
(1)在非高斯随机控制方面,在已有的科研成果中,设定值通常取做常值或预先设定的参考轨迹。但事实上,设定值/参考轨迹需要根据实际经济性能指标进行设定,是随时改变的一条运动轨迹,因而这一处理方案有待进一步考虑的;
(2)同样的,还是在非高斯随机控制方面,已有文献主要是构建基于统计信息的瞬时性能指标,得到单步控制律,不能预测未来的控制作动趋势,具有一定的保守性;
(3)在传统的实时优化和模型预测控制两层经济性过程优化中,实时优化没有考虑过程的动态特性,因此不能保证所得到设定值的可达性。
发明内容
本发明解决的技术问题是:针对典型的化工反应过程在遭受非高斯干扰(反应堆温度及反应物浓度)影响的情况下,设计一个两层结构的动态经济性能优化与模型预测控制策略。上层充分考虑系统的动态特性,根据经济效益指标设定优化目标函数(反应过程中的损耗函数),确定最优操作点即最优的温度和浓度值,作为下层模型预测控制动态优化的设定值或参考轨迹。下层模型预测控制采用基于统计信息的非高斯随机控制策略跟踪上层优化得到的设定值(反应堆浓度和温度)或温度和浓度的变化轨迹。
为了解决上述技术问题,本发明通过以下步骤实现:
步骤一、非高斯扰动下的非线性随机系统模型描述。
非高斯非线性随机系统可描述为如下的状态空间形式:
xk+1=f(xk,uk,ωk)
其中,是状态向量(例如化工反应过程中的反应堆温度、反应物浓度),是操作输入向量(热量、体积流量),是随机干扰向量(反应堆温度和浓度产生的偏差),可以是高斯分布也可以是非高斯分布。输入满足如下的约束条件定义系统的实际状态轨迹xk和参考轨迹之间的误差为
步骤二、建立两层动态优化控制框架结构,如图1所示。
两层动态优化控制策略的实施过程如下:在k时刻,上层的经济模型(化工过程反应产生的损耗)预测控制(economic model predictive control,EMPC)接收过程的状态反馈并通过求解经济性能优化问题得到经济性最优状态轨迹即温度和浓度的最优设定。下层的模型预测控制(model predictive control,MPC)通过优化基于统计信息的性能指标函数得到最优控制律,使得过程状态(温度和浓度)跟踪上层得到的最优经济性时变状态轨迹。记上层EMPC和下层MPC的预测时域分别为NE和N。因为下层MPC需要在整个运行周期和预测时域上均有可跟踪的最优轨迹,因此NE应当足够大,使得k+NE>k+N。
上述动态经济优化和控制过程可以归纳总结为以四步:
1)在k时刻,EMPC接收来自非高斯随机过程的系统状态xk;
3)下层MPC计算得到最优控制律,使得系统状态跟踪最优状态轨迹;
4)未来时刻重复上述步骤1)到3)。
步骤三、基于统计信息的性能指标选取,具体如下:
由于过程系统受非高斯随机干扰的影响,因此需要采用除均值、方差之外更一般的统计信息进行随机性度量。下面将对上层EMPC和下层MPC的性能指标分别进行介绍。
(1)上层EMPC性能指标的选取:
经济性能指标需要根据化工反应过程进行选取,如考虑经济效益,反应物转化损耗。这里将经济性度量用一般形式表示为L(xk,uk)。由于xk是非高斯随机变量,因此L(xk,uk)也是非高斯随机变量,需要考虑其随机特性。因此,EMPC性能指标可选取为
其中,E[·]和H[·]分别为随机变量的均值和二阶Renyi熵,可分别通过下面的非参数估计方法进行计算得到。
(2)MPC性能指标的选取:
在下层的跟踪控制部分,需要设计控制律(例如反应物的体积流量和热量的有效控制),使得过程状态(如实际的反应温度和浓度)跟踪上层优化得到的设定轨迹,也就是让ek尽可能接近零。考虑到其非高斯特性,这里引入生存信息势(Survival InformationPotential,SIP)构建性能指标函数,即
步骤四、动态经济性能优化与模型预测控制器设计,具体如下:
(1)上层EMPC优化问题表述如下
其中,U为控制输入的可行域,αi为充分小的正实数。其中,是与时间有关的经济性度量函数。第一个约束条件是用来得到过程状态的系统模型;第二个概率约束表示状态可以在一定的概率范围内超出其最大最小约束边界;第三个是对控制能量的约束。
(2)下层MPC可归结为如下的优化问题
采用可行的优化算法如随机梯度法,即可完成最优控制律的设计。
本发明与现有的技术相比,具有如下的优点:
(1)考虑到先进控制策略实施的经济性要求,本发明提出了两层的动态经济优化与控制(如典型化工过程的温度和浓度的跟踪控制)(EMPC-MPC)的一体化设计方案。与传统的两层RTO-MPC在上层采用稳态模型相比,本发明充分考虑了过程的动态特性,强调所得最优设定轨迹的可达性,更具有实际工程意义。
(2)由于实际工业过程不可避免受随机噪声影响,且一般为非高斯随机噪声,本发明在经济性能优化和控制器设计时均考虑其带来的影响,比已有方法中忽略噪声影响或假设噪声服从高斯分布更具有一般性和实际意义(如充分考虑到化工反应过程的经济效益以降低实际的损失)。
(3)本发明采用熵、SIP等统计信息充分刻画非高斯随机向量的随机特性,建立经济性能指标和跟踪控制指标。
(4)与已有的非高斯随机控制方法相比,本发明充分考虑了控制输入的随机特性,设计了多步预测控制策略。
附图说明
图1是本发明的动态经济性能优化与控制方法的结构框图。
具体实施方式
步骤一、一个非等温的连续均匀搅拌反应器系统(continuous stirred tankreactor,CSTR),一个基础的反应从A到B,反应堆的进料为纯的A,A的体积流速为F,温度为T0+ΔT0,摩尔浓度为CA0+ΔCA0。反应堆的夹套被用来为反应堆移除热量。根据能量守恒定律上述系统的动力学方程:
其中,CA是进料A的浓度,T是反应堆温度,Q是流体热量,VR是反应堆体积,ΔH是反应热量,k0和E是指数因子和反应活化能。Cp和q分别指比热容和反应器内流体密度。y=[T,CA],y1=T,y2=CA是系统输出,u1=q和u2=F操作输入,ΔCA0和ΔT0为进料浓度和温度的非高斯随机扰动。
定义系统的实际输出y和参考轨迹yE之间的误差为e=y-yE。
步骤二、对上述系统进行离散化,根据下面的步骤建立两层动态优化控制框架结构:
1)在k时刻,上层EMPC接收来自非高斯随机过程的系统输出yk;
3)下层MPC计算得到最优控制律,使得系统输出跟踪最优设定轨迹;
4)未来时刻重复上述步骤1)到3)。
步骤三、基于统计信息的性能指标选取,具体如下:
(1)上层EMPC性能指标的选取:
考虑到进料投入、产物和反应温度因素,给出经济成本函数:
L(yk,uk,k)=A1u1 2(k)-A2u2(CA0-CA)+A3(y1(k)-365)2
其中,A1,A2,A3为权值参数,考虑到系统输出的随机特性,要使上述经济成本函数最小,选定经济性能指标为:
其中,和分别是输出yk+i的二阶Renyi熵和期望值,是Lk+i的概率密度函数。这里,采用非参数估计法计算得到性能指标的值。设是经济成本函数Lk+i的200个可能的采样值,则通过Parzen窗方法可得到二阶Renyi熵的计算公式:
经济成本函数Lk+i的期望可由下式计算得到
(2)下层MPC性能指标的选取:
下层MPC的目的是设计控制律,使得系统输出跟踪上层优化得到的设定轨迹,因此性能指标选择为
步骤四、动态经济性能优化与模型预测控制器设计,具体如下:
上层EMPC可通过求解下列优化问题实现:
下层MPC可归结为如下的优化问题
采用随机梯度法,即可完成非等温的连续均匀搅拌反应器系统温度和浓度的最优控制。
Claims (1)
1.非高斯系统的动态经济性能优化与控制一体化设计方法,其特征是包括:
步骤一、建立非高斯扰动下的非线性随机过程系统模型;
xk+1=f(xk,uk,ωk)
其中,f(·)是表示系统动态特性的非线性向量函数,是系统的状态向量,即化工反应过程中的反应堆温度、反应物浓度,表示实nx维列向量空间,是操作输入向量,即热量、体积流量,表示实nu维列向量空间,是随机干扰向量,即反应堆温度和浓度产生的偏差,表示实nω维列向量空间;操作输入向量满足如下的约束条件ui是uk的第i个元素,是ui的边界值;系统的状态向量xk与经济性最优状态轨迹之间的误差为
步骤二、建立两层动态优化控制框架结构:
1)在k时刻,上层的经济模型预测控制接收来自非高斯随机过程系统的状态向量xk,k是即时时间;
3)下层的模型预测控制计算得到最优控制律,使得系统状态跟踪最优状态轨迹;
4)未来时刻重复上述内容1)到3);
步骤三、基于统计信息的性能指标选取,具体如下:
(1)上层的经济模型预测控制性能指标的选取:
其中,L(xk+i,uk+i)是k+i时刻上层经济模型预测控制的经济性能指标函数,E[·]和H[·]分别为随机变量的均值和二阶Renyi熵,分别通过下面的非参数估计方法进行计算得到;
(2)下层的模型预测控制性能指标的选取:
采用生存信息势构建性能指标函数,即
步骤四、动态经济性能优化与模型预测控制器设计,具体如下:
(1)上层经济模型预测控制优化问题表述如下
其中,是上层经济模型预测控制中第k+i时刻的操作输入向量,和分别是的下界值和上界值,αi为充分小的正实数;一个约束条件是用来得到过程状态的系统模型;另一个约束条件表示状态在一定的概率范围内超出其最大最小约束边界;
(2)下层模型预测控制归结为如下的优化问题
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711400399.2A CN107942682B (zh) | 2017-12-22 | 2017-12-22 | 非高斯系统的动态经济性能优化与控制一体化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711400399.2A CN107942682B (zh) | 2017-12-22 | 2017-12-22 | 非高斯系统的动态经济性能优化与控制一体化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107942682A CN107942682A (zh) | 2018-04-20 |
CN107942682B true CN107942682B (zh) | 2020-10-30 |
Family
ID=61942162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711400399.2A Active CN107942682B (zh) | 2017-12-22 | 2017-12-22 | 非高斯系统的动态经济性能优化与控制一体化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107942682B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108446529B (zh) * | 2018-06-22 | 2020-03-31 | 太原理工大学 | 基于广义互熵—dpca算法的有机朗肯循环系统故障检测方法 |
CN110703599B (zh) * | 2019-09-17 | 2022-06-07 | 太原理工大学 | 基于动态数据校正的有机朗肯循环系统控制性能优化方法 |
CN110609476B (zh) * | 2019-09-19 | 2022-04-12 | 太原理工大学 | 一种基于高斯过程模型的多变量非线性动态系统模型预测控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794117A (zh) * | 2010-02-11 | 2010-08-04 | 哈尔滨工业大学 | 一种基于机理模型的优化控制膜生物反应器运行的方法 |
CN103176408A (zh) * | 2013-02-01 | 2013-06-26 | 浙江大学 | 基于双层递阶结构的聚丙烯生产过程牌号切换控制方法 |
CN103472723A (zh) * | 2013-08-19 | 2013-12-25 | 上海交通大学 | 基于多模型广义预测控制器的预测控制方法及系统 |
CN105515029A (zh) * | 2015-12-03 | 2016-04-20 | 华北电力大学 | 飞轮储能系统的控制方法及装置 |
CN106288654A (zh) * | 2016-07-27 | 2017-01-04 | 杭州杭氧股份有限公司 | 一种氧氮精馏外压缩空分设备的快速变负荷优化控制方法 |
US9568901B2 (en) * | 2012-08-27 | 2017-02-14 | Nec Corporation | Multi-objective energy management methods for micro-grids |
CN107180279A (zh) * | 2017-06-14 | 2017-09-19 | 重庆科技学院 | 基于qpso‑dmpc的反应再生系统优化控制方法 |
-
2017
- 2017-12-22 CN CN201711400399.2A patent/CN107942682B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794117A (zh) * | 2010-02-11 | 2010-08-04 | 哈尔滨工业大学 | 一种基于机理模型的优化控制膜生物反应器运行的方法 |
US9568901B2 (en) * | 2012-08-27 | 2017-02-14 | Nec Corporation | Multi-objective energy management methods for micro-grids |
CN103176408A (zh) * | 2013-02-01 | 2013-06-26 | 浙江大学 | 基于双层递阶结构的聚丙烯生产过程牌号切换控制方法 |
CN103472723A (zh) * | 2013-08-19 | 2013-12-25 | 上海交通大学 | 基于多模型广义预测控制器的预测控制方法及系统 |
CN105515029A (zh) * | 2015-12-03 | 2016-04-20 | 华北电力大学 | 飞轮储能系统的控制方法及装置 |
CN106288654A (zh) * | 2016-07-27 | 2017-01-04 | 杭州杭氧股份有限公司 | 一种氧氮精馏外压缩空分设备的快速变负荷优化控制方法 |
CN107180279A (zh) * | 2017-06-14 | 2017-09-19 | 重庆科技学院 | 基于qpso‑dmpc的反应再生系统优化控制方法 |
Non-Patent Citations (3)
Title |
---|
Single Neuron Stochastic Predictive PID Control Algorithm for Nonlinear and Non-Gaussian Systems Using the Survival Information Potential Criterion;Ren, Mifeng等;《ENTROPY》;20160630;第18卷(第6期);1-16 * |
Statistical Information Based Single Neuron Adaptive Control for Non-Gaussian Stochastic Systems;Ren, Mifeng等;《ENTROPY》;20120731;第14卷(第7期);1154-1164 * |
非高斯系统的控制及滤波方法研究;任密蜂;《中国博士学位论文全文数据库(电子期刊)信息科技辑》;20141215(第12期);I140-72 * |
Also Published As
Publication number | Publication date |
---|---|
CN107942682A (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104375475B (zh) | 间歇反应器中间歇反应过程的优化控制方法 | |
CN107942682B (zh) | 非高斯系统的动态经济性能优化与控制一体化设计方法 | |
Mujtaba et al. | Neural network based modelling and control in batch reactor | |
CN102339040B (zh) | 精馏塔优化控制方法 | |
Zou et al. | Design of fractional order predictive functional control for fractional industrial processes | |
Kano et al. | Inferential control of distillation compositions: selection of model and control configuration | |
CN114995155B (zh) | 一种高纯度精馏过程鲁棒解耦控制系统及控制方法 | |
Tian et al. | Hybrid modeling of molten steel temperature prediction in LF | |
CN110221540B (zh) | 基于Hammerstein模型的连续搅拌反应器系统控制方法 | |
Mallaiah et al. | A simulated annealing optimization algorithm based nonlinear model predictive control strategy with application | |
Tajjudin et al. | Comparison between optimally-tuned PID with self-tuning PID for steam temperature regulation | |
CN110750049B (zh) | 带有时滞和扰动的间歇过程2d预测容错控制方法 | |
Konakom et al. | Neural network‐based controller design of a batch reactive distillation column under uncertainty | |
Shettigar J et al. | Development and validation of advanced nonlinear predictive control algorithms for trajectory tracking in batch polymerization | |
CN111221252A (zh) | 一种针对带分数阶滞后过程工业系统的预测控制器参数解析整定方法 | |
Venkatesan et al. | Controller design based on model predictive control for a nonlinear process | |
Zhou et al. | Hybrid intelligent control scheme of a polymerization kettle for ACR production | |
Zhengshun et al. | ARX-NNPLS model based optimization strategy and its application in polymer grade transition process | |
Karacan et al. | Multivariable system identification and generic model control of a laboratory scale packed distillation column | |
Rohman et al. | Performance evaluation of adaptive based model predictive control for ethylene glycol production from dimethyl oxide hydrogenation | |
CN113568309B (zh) | 用于温度场的在线时空控制方法 | |
Habobi et al. | Dynamics, Simulation, and Control of a Batch Distillation Column using Labview | |
CN203217277U (zh) | 槽式反应器基于在线支持向量机的非线性模型预测控制器 | |
Zhou et al. | A just-in-time learning based integrated IMC-ILC control strategy for batch processes | |
CN104122878A (zh) | 工业节能减排控制装置及控制方法 |
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 |