第 32 卷 第 1 期
2018. 03
沈 阳 化 工 大 学 学 报
JOURNAL OF SHENYANG UNIVERSITY OF CHEMICAL TECHNOLOGY
Vol. 32 No. 1
Mar. 2018
文章编号:
2095 - 2198( 2018) 01 - 0001 - 08
化工过程实时优化策略研究综述
袁德成1,2 , 潘多涛1,2 , 曾 静1,2
( 1. 辽宁省化工控制技术重点实验室,辽宁 沈阳 110142;
2. 沈阳化工大学 信息工程学院,辽宁 沈阳 110142)
摘 要: 综述化工过程实时优化方法的研究现状和发展动态,内容包括工业自动化系统的一般
结构、实时优化基本原理、实时优化的设计方法及需解决的一些问题等. 实时优化功能介于工厂自
动化系统的计划 / 调度与基础控制功能中间,扮演实现预期经济目标的角色,即在满足安全、质量、
环保和设备约束的条件下,降低消耗或提高利润. 实时优化是把实际测量信息反馈到基于模型的
最优规划问题中,借助串联结构,实现包括多个相关基础反馈回路在内的最优反馈控制功能. 将测
量信息用于提高模型准确度、修正最优规划问题、直接构造控制输入等,不断发展形成了实时优化
的不同设计方法. 算法的收敛性、在线性能评估、集成化应用支撑技术开发等代表了实时优化发展
的方向.
关键词: 实时优化; 在线优化; 模型适应; 非线性规划
doi: 10. 3969 / j. issn. 2095 - 2198. 2018. 01. 001
中图分类号:
TP273
文献标识码: A
现代化学工业过程自动化系统的主要任务
包括建模、控制、优化和监视等. 其中,控制任务
力发电[5]、电网能源管理[6]等新领域也显示出
应用潜力. RTO 持续获得重视的原因在于这样
是调整关键的工艺变量维持其设定值不变或跟
的技术能带来可观的经济收益,尤其对大规模企
踪其变化; 优化任务则是追求局部或全局工艺在
严格约束条件内实现经济收益类指标的最佳. 同
时配置控制和优化任务的架构有非线性模型预
测控制[1]、经济性模型预测控制等模式[2],但这
样的安排在建模、计算和性能监视等方面所需开
业,即使在产量上 1 % 的改善也可能在年收益上
获益巨大[7]. RTO 的设计目标是在满足安全、质
量、环保和设备约束的条件下,降低消耗或提高
利润,在当今市场竞争激烈的条件下,这样的应
用技术正是众多投资者关注的重点.
销大. 更实用的方法是把控制和优化任务分层处
本文从 工 业 自 动 化 系 统 的 一 般 结 构、RTO
置,即由位于上层的实时优化功能( real-time op-
基本原理、模型更新机制、系统构造等方面,综述
timization,缩写为 RTO) 依据某类性能指标和约
RTO 的各种设计方法及优缺点,最后将需要解
束条件,通过求解最优规划给出关键工艺变量的
决的一些问题和结论进行说明.
设定值,再由位于下层的控制( 回路) 功能实现
在扰动或不确定性环境中跟踪这些设定值的最
优轨迹. 这样的配置在目前工业应用广泛的 DCS
系统上也相对容易实现. RTO 在流程工业包括
炼油、化工、食品、生物、造纸、矿冶等领域的应用
一直受到重视[3],近年来在楼宇节能[4]、高空风
收稿日期:
2017 - 12 - 08
基金项目: 国家自然科学基金项目( 61503257)
1 工业自动化系统的架构设计
一个复杂技术系统的结构设计至关重要. 工
业自动化系统的结构设计除考虑自身的硬软件
体系设计之外,被控过程的工艺和装备配置、控
制系统组成、工艺与控制并行交互设计等都是成
作者简介: 袁德成( 1960 - ) ,男,内蒙古阿拉善人,教授,博士,主要从事化工过程建模、控制与优化的研究.
2
沈 阳 化 工 大 学 学 报
2018 年
功实现预期目标必须统筹考虑的重要事项,已经
引起相关领域学者的研究兴趣[8 - 9]. 基础控制层
的基本任务是选择被控变量、操作变量和确定它
集成了来自工艺装置的生产实时数据和供应链
管理的数据,某一层的操作则是利用从下一层获
取的信息和接受上一层的指令,通过某种优化策
们之间如何配对,有系统的经验法则可遵循. 当
略计算给出下一层的设定值,即基于串级-反馈
上层实施 RTO 功能时,两层的设计需要统筹考
控制的思想. 目前支撑分层型过程自动化系统架
虑,最终决策也并不唯一. 在 RTO 系统中包含的
构功能实现的有效方法是基于稳态模型的手段,
约束条件区分为等式约束和不等式约束两类,优
但如果各功能层所依赖的模型不一致,就有可能
化操作是在安全、质量和产量要求( 硬性) 确保
无法实现既定目标. 这也是有关实时优化技术中
实现后再追求的设计目标,因此,首先是从操作
被研究者关注最多的问题.
变量集合中选择出一些控制变量确保满足等式
约束条件,其次是用剩余的操作变量实现 RTO.
不等式约束的处理相对复杂,主要是因为在最优
操作点上可能出现等式约束和不等式约束两种
情形,分 别 被 称 为 活 动 型 ( active ) 和 消 极 型
( inactive) [10]. 如果活动型约束在整个控制进程
中与工艺条件和扰动无关,则可以把它等同为前
述的等式约束,用于控制硬性指标; 若活动型约
束集合在控制进程中有变化,则需要变更控制结
构,即一些回路伺服跟踪,一些回路持续调节,从
图 1 分层型过程自动化系统
原有的控制结构快速过渡到新的最优操作. 在自
动化系统结构设计时,充分考虑约束条件的影响
并加以有效利用,能得到事半功倍的效果. 在设
Fig. 1 Hierarchical plant automation systems
2 RTO 基本原理
计阶段如何从约束条件集合中检出活动型约束,
RTO 问题的原型属于基于模型设计优化控
在控制进程中如何在线监视约束集合活动性的
制策略,其基本组成如图 2 所示. RTO 与一般最
变迁等问题,都是值得深入研究的课题.
分层型过程自动化系统架构,其一般组成如
图 1 所示[11],其中实时优化位于基础控制层和
计划调度层之间,协调计划调度任务和生产装置
的操作. 位于最上层的计划功能,生成经济预报
和生产目标,而调度功能则决定如何执行所选择
的计划,这层功能实现的关键是确定可行解域;
位于最下层的基础控制层则利用各类闭环控制
算法,例如 PID、MPC 等,实现安全、产量和质量
等基本生产目标,这些控制回路的设定点一般是
具体的工艺参数,如温度、压力、流量、液位和物
质组分等,与一个企业追求的经济指标之间并没
有直接显性的关系. 实时优化层任务的重点在于
考虑装置运行的经济效益,即控制目标可能是能
耗最小或收益最大,而要实现这样一个要求往往
需要基础层若干个控制回路协同动作. 时间尺度
从上层的周或天到底层的秒或毫秒级不等,空间
尺度可能跨越单个装置、全流程甚至更广,信息
优规划问题的不同之处在于优化算法中在线融
合现场测量信号,包括性能指标及约束条件计
算,即用实际测量改善不确定性影响.
图 2 RTO 系统的基本功能块组成
Fig. 2 Basic function blocks of a RTO system
3
( 1)
目标函数
min
u
Φp( u) = ( u,yp( u) )
约束条件
G p,i( u) = gi( u,yp( u) ) ≤0, i = 1,…,ng
其中下标 p 表示实际的被控对象,决策变量
( 设定点) u ∈ Rnu,输出变量 yp ∈ Rny,性 能 指 标
Φp: Rnu × Rny→R,输入输出上的约束条件G p: Rnu
× Rny→Rng. 在实际问题中,过程真实的输入输
出映射关系并不能精确得到,往往只有近似模型
( 2)
可用:
y = f( u,θ)
( 3)
其中模型输出 y∈Rny,模型参数 θ∈Rnθ,输
入输出映射 f: Rnu × Rnθ →Rny. 借助这个模型,可
以构造基于模型的优化问题( B) 求取近似解:
目标函数
min
u
Φ( u,θ) = ( u,f( u,θ) )
( 4)
约束条件
G i ( u,θ) = gi ( u,f( u,θ) ) ≤0, i = 1,…,ng
( 5)
在上述目标函数和约束条件中,均包含了未
知参数 θ,一般假设: 函数 和 f 连续,结构已知
且结构不随参数 θ 的估计而变化.
对于模型问题 B,最优解存在的一阶必要条
2. 1 RTO 算法中信息预处理
RTO 一般包括模型更新和基于模型的优化
两个模块,但原始数据的处理和优化计算结果施
加前的深度分析都是 RTO 成功应用必须完成的
基础性工作. RTO 是基于系统稳态模型的优化
操作[12],即通过测量数据分析、判断系统到达了
稳态才开始启动 RTO. 稳态检测一般基于统计
过程控 制 方 法,例 如 用 Shewart 图、Hotelling T2
检验等监视过程变量的波动是否在其控制线内,
如果相邻时刻的样品均值在统计上相等,则视为
过程处于稳态. 在工业应用场合,来自过程的测
量数据常受到高频、零均值随机噪音的污染,需
要用低通滤波器,例如滑动平均滤波、指数加权
滑动平均滤波等方法,降低噪音的影响,为后续
模型更新提供高质量的有效数据. RTO 一般优
化的是经济指标,往往需要物质成分的测量分析
数据,人为操作因素引起的过失误差或某一时刻
数据丢失的事件都可能发生,必须通过统计测量
调整检出过失误差数据并加以剔除. 数据调整则
是通过求解二次优化问题,把测量数据调整到满
足物质和能量守恒定律( 稳态模型) [13].
在测量数据验真之后进行的模型更新,即调
整稳态模型中某些参数以适应当前采集到的测
量数据,原理简单,但具体实现要处理的棘手问
题很多,例如: ( 1) 在稳态模型中选择哪些参数
进行适应调整; ( 2) 依靠当前的测量数据这些参
数是否可以辨识( 激励信号是否充分) ; ( 3) 这些
参数的变化对后续进行的经济优化问题影响是
否敏感; ( 4) 测量这些参数的传感器位置和精度
是否适当等问题. 把数据调整和参数估计结合在
一起,同时求解也是可行的方案,但伴随的计算
病态问题需要在计算进程中加以诊断. 如果同一
过程积累了多组数据,也可以用于提高参数估计
第 1 期
袁德成,等: 化工过程实时优化策略研究综述
件( NCO) ,也称为 KKT 条件:
G i ( u,θ) = gi ( u,f( u,θ) ) ≤0,
i = 1,…,ng
μi≥0, i = 1,…,ng
Φ( u,θ)
u
ng
+ ∑
i = 1
μi
G i ( u,θ)
u
= 0
μi G i ( u,θ) = 0, i = 1,…,n
对 于 实 际 问 题 A,也 可 导 出 形 式 相 似 的
KKT 条件[14]. 如果进一步假设问题 A 只具有一
个全局最优解,一个至关重要的问题是: 通过迭
g
( 6)
的精度.
代求解问题 B 得到的最优解,最终能否逼近问
一旦得到更新后的模型参数,则可以通过求
题 A 的最优解. 对于所有基于模型的控制和优
解非线 性 规 划 问 题,计 算 出 控 制 器 的 最 优 设
化问题,这是其解决方案必须回答的共性问题.
定值.
2. 2 RTO 基本问题描述
要匹配 KKT 条件,除约束条件外,还涉及到目标
函数梯度和约束条件梯度的匹配问题. 通过输出
变量的测量值准确估计函数梯度值是近几年研
针对一个实际控制对象,实施的优化问题
究的一个热点.
( A) 表述如下:
4
沈 阳 化 工 大 学 学 报
2018 年
2. 3 RTO 适用条件与限制
解可能不可行. 而且,哪些模型参数能用参数估
计的方法不断适应,哪些模型参数要保持常数等
RTO 适用的场合包括: ( 1) 除了优先保证
问题,并没有明确直接的确定方法. 关键的模型
工艺安全、产率和质量指标实现的控制变量之
参数应该能从可用的测量数据加以辨识,但非线
外,还存在许多工艺变量可用于优化操作,即可
性模型形式、数据质量不佳、参数搜索边界过大
操控的变量个数要远大于需要控制的目标个数;
等因素将导致参数估计复杂化. 因此,在设计和
( 2) 用于优化操作的变量对利润的影响要显著,
研究 RTO 方法时一个感兴趣的问题是在增强可
如果找不到这样的变量,则实时优化无从做起;
( 3) 扰动出现的间隔与进行 RTO 决策需要的时
行性和最优性时减轻对模型精度以及模型更新
的需求[15].
间开销应该在同一时间尺度上; ( 4) 确定最优操
二是模型更新与优化的相互作用问题. 如果
作点的程序不能太复杂.
模型和实际对象在结构上完全一致,仅在参数上
使用 RTO 的限制包括: ( 1) 批过程没有运
存在失配但可辨识,则经过模型更新和优化可以
行稳态,用于连续操作的 RTO 方法不适合这样
收敛到实际对象的最优点. 如果在结构上存在失
一大类问题; ( 2) 在 RTO 方法实施中需要过程
配,即使两步法计算最终收敛,但是不是达到了
处于稳态操作,可能需要花费足够长的时间等待
实际对象的最优点无法保证. 注意在问题 B 中,
下一个稳态的到来; ( 3) RTO 方法优化的是稳
未知参数 θ 同时出现在目标函数和约束条件,故
态经济指标,在过程操作变迁中的损失或利益则
通过更新模型一般无法同时影响目标函数和约
没有考虑. 有扰动才需要在线的优化操作,扰动
束条件. 因此,对于规划问题 B,重新定义模型准
出现则使得稳态模型与真实过程特性之间必然
失配,RTO 方法的复杂性和工业应用的限定条
确度是必要的. 除模型和对象匹配外,优化问题
中相关函数的一阶、二阶导数也需要匹配[16].
件基本源于此.
3 基于过程模型适应的 RTO 设计
在 2. 2 节中论述的 RTO 经典两步递推设计
法中,更新模型需要进行参数估计,有了模型之
后再进行优化计算得到新的操作点,即 RTO 依
赖过程模型完成最优操作条件的计算.
在第 k 次递推时,通过实际输出测量值 yp
进行参数估计:
θk = arg min
θ
yp( uk) - f( uk,θ)
( 7)
一旦得到参数 θk,则可求解最优问题 B 获
得当前的最优操作指令. 原理简明,但要成功实
施除数据验真和最优指令等模块需要仔细设计
外,以下问题也必须考虑:
一是来自真实过程的简化模型表示以及不
确定性. 用工厂操作数据更新过程模型,对过程
输出预报会更加准确. RTO 系统性能取决于模
型表达过程行为的准确程度. 问题是开发精确模
型往往很困难. 在有模型失配和不可测扰动时,
RTO 提供的最优解可能是实际过程的次优解.
更糟糕的是当约束不能被模型准确预报时,最优
三是由模型 / 对象对应的 KKT 条件匹配问
题而引出的模型结构设计问题. 如果参数更新能
导致模型 / 对象对应的 KKT 条件完全匹配,则
RTO 系统能保证收敛到实际问题的最优点,即
在参数更新时也要考虑如何能比较出两组 KKT
条件是否匹配的问题. 待估计参数的个数、函数
梯度估计、模型结构形式选择等问题,都需要在
模型更新时加以考虑,以留下足够的柔性方便后
续 KKT 条件是否匹配的判定[17].
4 基于优化模型适应的 RTO 设计
4. 1 集成系统优化与参数估计( ISOPE) 方法
先考虑仅在控制输入上有上下限约束的情
况: uL ≤u≤uU . 在第 k 次 RTO 递推时,用已知的
控制输入 uk 和对象输出 yp ( uk ) ,在模型匹配条
件 y( uk,θk) = yp ( uk ) 下,求解参数估计问题,得
到更新后的参数 θk. 进一步假设对象输出梯度
yp( uk)
存 在,在 目 标 函 数 梯 度 计 算 中 引 入
u
ISOPE 修正子 λk∈Rnu:
第 1 期
袁德成,等: 化工过程实时优化策略研究综述
5
( λk) T =
y
( uk,y( uk,θk) )
[
yp( uk)
u
-
]
y( uk,θk)
u
( 8)
i = 1,…,ng. 这些修正项采用了测量值和预测值
之差的形式,文献[14]证明了在一定条件下如
果问题 E 的解收敛,则其解也是实际问题 A 的
基于参数 θk 和修正子 λk,可以构造修正的
一个 KKT 点. 在实际对象 / 模型失配的环境中,
修正项适应的方法与 3 节的模型适应方法相比
具有显著的优点,不需要直接更新模型,付出的
代价是需要估计目标函数和约束条件的梯度. 梯
度估计的方法目前主要有稳态摄动法和动态摄
动法两类[19],稳态摄动法至少需要( nu + 1) 稳态
操作点才能估计出有关梯度,而动态摄动法则是
利用相邻稳态过渡期间的变化信息估计梯度. 如
果两个 RTO 递推的等待时间较长,用动态摄动
法估计梯度比较合适,但缺点是要求稳态过渡期
间的数据包含足够的信息量.
近几年关于修正子适应方法的研究比较活
跃,在优化问题( E) 基础上提出了输出变量修正
子适应方法[20]、方向导数修正子适应方法[21]、
二阶修正子适应方法[22],以及应用于批过程的
动态实时优化策略等[23].
5 基于控制适应的 RTO 设计
前述的模型适应方法着眼点放在参数精确
估计方面,而修正项适应方法则把模型失配问题
统筹放在规划问题架构下求解. 不论采用何种方
法,目的都是要计算出决策变量轨迹,而本节论
述的控制适应方法正是直接针对决策变量设计
以实现 RTO 目标的策略[14].
5. 1 自寻优控制[24]
如果在有扰动时不改变受控变量的设定值,
仍然能实现最优性能的损失在可接受的范围内,
则这样的控制结构被称为自寻优控制( self-opti-
最优问题( C) :
u*
k + 1 = arg min
u
( u,y( u,θk) ) + λT
k u
( 9)
约束条件: uL ≤u≤uU
在具体实施中,为防止控制的过激作用,常
用一阶指数滤波器:
uk + 1 = uk + K( u*
如果存在过程约束 gi ( u,y( u,θk ) ) ≤0,i =
k + 1 - uk)
( 10)
1,…,ng,同样可构造修正子:
( λk) T = (
y
( uk,y( uk,θk) ) +
g
y
μT
k
[
( uk,y( uk,θk) ) ) ·
yp( uk)
u
-
y( uk,θk)
u
]
类似的可以构造最优问题( D) :
u*
k + 1 = arg min
( u,y( u,θk) ) + λT
k u
u
( 11)
( 12)
约束条件: gi ( u,y( u,θk) ) ≤0,i = 1,…,ng
4. 2 修正子适应方法
在不确定性和扰动存在的情况下,如何使静
态模型更准确一直是 RTO 应用的主攻方向. 正
如 3 节所述,参数的准确估计受到诸多因素限
制,近年来出现了新的研究思路,即在最优问题
B 的目标函数和约束条件中,分别加上适当的线
性修正项,重构优化问题,不再进行直接的模型
更新计算[18]. 加上修正项后的优化问题( E) 定
义为:
u*
k + 1 = arg min
u
Φm,k( u) = Φ( u) +
( λΦ) T ( u - uk)
约束条件:
G m,i,k( u) = G i ( u) + εG i
( λG i
k +
k ) T ( u - uk) ≤0, i = 1,…,ng
其中,在目标函数中引入梯度修正项 ( λΦ
[
Φp( uk)
u
-
Φ( uk)
u
]
束修正项 εG i
k = G p,i ( uk) - G i ( uk) ,i = 1,…,ng 和
约束梯度修正项( λG i
k ) T =
G p,i ( uk)
u
-
G ( uk)
u
,
( 13)
mization control) . 在分层控制结构中,这样的控
制结构设计非常有意义,因为在基础控制层有众
多的基本控制回路,其设定值一般可以作为上一
层的自由度,即优化层的控制变量. 从经济角度
看这些回路的设定值需要给出最佳设定,但问题
是: 全部的设定值都用于优化操作,还是只选择
( 14)
k ) T =
,而在 约 束 条 件 中 引 入 约
其中一部分; 选择应遵守什么样的原则等并不明
晰. 在基础控制层,按照分散化控制架构,许多学
者给出了控制结构设计的程序: ( 1) 选择受控变
量; ( 2) 选择控制变量; ( 3) 选择测量变量; ( 4)
6
沈 阳 化 工 大 学 学 报
2018 年
设计控制配置,即受控变量和控制变量的配对;
的效益大.
( 5) 控制算法或规律( 例如 PID、LQG 等) . 这样
和其他的全流程控制结构设计方法一样,自
的设计准则也可以用于由优化层和基础控制层
寻优控制设计给出的是设计方法,针对具体案
组成的串级型控制系统,但问题要复杂得多. 对
例,需要针对所有可能的扰动,评价不同控制配
于控制或调节而言,要实现的目标是某个具体的
置对性能指标损失的影响,显然要面对的优化问
物理量,容易确定; 对于优化而言,要实现的目标
题规模很大. 结合具体过程控制的机理和操作特
则需要以目标函数的形式给出. 在控制中,受控
征,控制配对有一些经验规则可以遵循,面向单
变量与其设定值之间的增益一般不能改变符号,
元操作的大型仿真工具 Aspen 也提供了相关支
而在优化中则与具体的场景有关,不一定要求增
撑平台. 另外,应用自寻优控制的一个公开问题
益符号不变.
是如何找到一组输出变量或其组合,以组成反馈
优化问题( F) 表述如下:
目标函数 min
u
Φ( u,d)
( 15)
约束条件 g( u,d) ≤0
其中 u 是 nu 维控制向量,需要决策,d 是扰
动. 在稳态情况下,决策变量是扰动的函数,目标
( 16)
函数可以进一步表达为:
min
u
Φ( u,d) = Φ( uopt ( d) ,d) = Φopt ( d)
( 17)
对于假设已知的扰动d0 ,则可以应用各种非
线性优化算法求解,得到最优设定值:
c s = uopt ( d0 )
实际操作变量 u 和最优操作变量 uopt ( d) 作
用到受控对象上时,产生的性能指标一般不同,
( 18)
用性能指标差或损失更能比较两种策略产生的
控制系统.
5. 2 极值搜索控制[25 - 26]
极值 搜 索 控 制 ( extremum-seeking control,
ESC) 是非线性自适应控制和优化算法的融合,
以简单的结构把先进控制技术和 RTO 技术合并
为单一技术,在驱动过程朝着估计利润最大化方
向发展的同时,也实现了参数的连续估计. 从外
部施加的正弦信号,持续激励 ESC 系统往极值
点移动,而高通滤波器又抑制了低频信号在回路
中传播,对最终收敛到最优点起到了保证作用.
ESC 已在小规模系统仿真中显示出非常好的性
能,需要做的工作是在工业现场加以应用和检
验.
效果:
5. 3 最优必要条件( NCO) 跟踪控制[27 - 29]
L( u,d) = Φ( u,d) - Φ( uopt ( d) ,d)
反馈控制对性能指标 Φ 的影响可以进一步
( 19)
表达为:
ΔΦ =[Φ( unom ,d = 0) - Φ( unom ,di) ]+
[Φ( unom ,di ) - Φ( uopt,di) ]+
[Φ( uopt,di ) - Φ( ucon,di ) ]
其中第一项表示出现扰动时操作变量不变
( 20)
产生的性能指标损失; 第二项表示用最优操作策
略主动适应扰动时产生的性能指标差; 第三项表
示用最优操作和反馈控制补偿扰动时产生的性
能指标差. 如果上式中第一项比后两项之和还要
大得多,则表明实施 RTO 几乎没有经济收益; 如
果这三项均很小,则表明实施最优操作和反馈控
制对于抑制扰动、保持性能指标满意的作用不明
显; 如果第三项比较显著,这是应用 RTO 最好的
情形,即自适应调整设定值比固定设定值能获取
自寻优控制把优化问题转换为反馈控制系
统问题,即在有扰动进入或模型失配时不是变更
设定值,而是让反馈回路发挥作用. 对于许多非
线性过程,挑选某些输出变量或其线性组合以实
现自寻优控制往往会遇到困难. 最优必要条件跟
踪控制是选择活性约束函数以及其梯度函数的
重构为自寻优控制的输出变量,令其设定值为
零,通过反馈控制的方式直接强迫一阶最优必要
条件实现. 这类控制方法需要确定活动性约束集
合、对活动性约束进行跟踪以及推动梯度函数趋
向零等步骤.
基于控制适应的 RTO 实现方法避免了在线
求解优化数值解的时间开销,但这类方法需要活
动性约束集合在优化进程中保持不变. 如果活动
性约束集合是变化的,则监督和调整将使基于控
制适应的方法变得更加复杂.
第 1 期
袁德成,等: 化工过程实时优化策略研究综述
7
6 发展方向
( 1) RTO 系统的收敛性分析. 迄今为止,关
于 RTO 系统收敛性的研究进展缓慢. 由于 RTO
系统是通过不断递推逼近最优点,在有外部扰动
进入或实际对象 / 模型存在失配的环境中,逼近
最优点的进程可能受到干扰,最终能否到达实际
优化问题的最优点目前还缺乏一般结论,而这又
是 RTO 实现在线闭环的关键. 另外,在最优性质
方面也有许多问题需要探索,例如在模型适应两
步法中,是通过提高模型精度逐步逼近实际问题
的最优点,但目前关于模型精度和收敛到全局最
优点之间的关系还没有建立. 在几乎所有 RTO
系统的实现方法中,当实际对象 / 模型存在失配
时,匹配实际对象 / 模型的 KKT 条件是非常复杂
的,而要判断 RTO 系统是否能达到全局最优,这
又是必须完成的工作.
( 2) RTO 系统的性能评估技术. 在 RTO 系
统实施进程中,当过程内部特性变迁或有外部扰
动进入时,RTO 系统的性能可能变坏,此时需要
重新设计. 像普通控制器性能评估技术[30]一样,
RTO 系统也需要性能评估技术,这是保证 RTO
长期运行的关键技术,否则在投运初期,RTO 系
统可能运行良好,一旦维护不到位,不但品质下
降,也将影响工业界认可 RTO 技术的程度. 这方
面相关研究文献较少[31],有值得探索的空间.
( 3) RTO 设计、运行维护的可视化平台支
撑技术. 实践已表明: 只有工业实际应用有迫切
需求且能被现场工程师理解掌握的先进控制技
临前所未有的机遇和挑战. 过程强化技术正在变
革传统工艺路径和装备结构的设计理论和方法,
过程控制应用建模、仿真、控制和优化技术既保
证了工艺与装备设计基本目标的实现,也为提升
物质和能源效率发挥了越来越重要的作用. 虽然
RTO 提出至今已近 40 年,可从研发、设计、实施
到维护,均需要专业人员才能完成,要设计出像
PID 控制那样工程师皆会用的“傻瓜型”RTO 控
制器,还需付出更多努力. 可以预期,市场的巨大
需求、相关学科的持续发展,都会促使 RTO 理
论、方法和技术不断完善,与 PID 和 M PC 一起
成为工业自动化的主流技术.
参考文献:
[1] COPP D A,HESPANHA J P. Simultaneous Nonlin-
ear M odel Predictive Control and State Estimation
[J]. Automation,2017,77( 3) : 143 - 154.
[2] ELLIS M ,DURAND H,CHRISTOFIDES P D. A
Tutorial Review of Economic M odel Predictive Con-
trol M ethods[J]. Journal of Process Control,2014,
24( 8) : 1156 - 1178.
[3] M ARCHETTI A G,FRANCOIS G,FAULWASSER
T,et al. M odifier Adaptation for Real-Optimiza-
tion—M ethods and Applications [J]. Processes,
2016,4( 55) : 1 - 35.
[4] ZAVALA V M . Real-Time Optimization Strategies
for Building Systems[J]. Ind Eng Chem Res,2013,
52( 9) : 3137 - 3150.
[5] COSTELLO S,FRANCOIS G,BONVIN D. Real-
Time Optimization for Kites[J]. IFAC Proceedings
Volumes,2013,46( 12) : 64 - 69.
术,才能在工业生产过程中得到长期稳定的使
[6] ELSIED M ,OUKAOUR A,YOUSSEF T,et al. An
用. 目前 RTO 技术主要还是掌握在学术界手中,
Advanced Real Time Energy M anagement System
不论是定制化设计,还是接受更多工业实践的检
for M icrogrids[J]. Energy,2016,114: 742 - 752.
验和评价,都需要集成参数估计、信息处理、非线
性规划、仿真、性能评价、非线性代数方程连续性
分析等 功 能 于 一 体 的 可 视 化 平 台,有 效 支 撑
RTO 系统的研发和应用.
7 结 论
[7] YOUNG R E. Petroleum Refining Process Control
and Real-Time Optimization [J]. IEEE Control
Syst. M ag. ,2006,26( 6) : 73 - 83.
[8] MARLIN T E,HRYMAK A N. Real-Time Operations
Optimization of Continuous Processes[J]. In AIChE
Symposium Series-CPC-V,1997,93: 156 - 164.
[9] BRDYS M A,TATJEWSKI P. Iterative Algorithms
综述了 RTO 系统的构造原理、研究现状以
for M ultilayer Optimizing Control[M]. London: Im-
及潜在的发展方向. 正在推进的节能降耗、低碳
经济、可持续发展等政策导向,使得流程工业面
perial College Press,2005: 37 - 66.
[10] BAZARAA M S,SHERALI H D,SHETTY C M .
8
沈 阳 化 工 大 学 学 报
2018 年
Nonlinear Programming: Theory and Algorithms
[M]. 3 rd ed. New Jersey: John Wiley and Sons,
Time Optimization[J]. Journal of Process Control,
2016,39: 64 - 76.
2006: 163 - 193.
[22] FAULWASSER T,BONVIN D. On the Use of Sec-
[11] FORBES J F,M ARLIN T E,YIP W S. Real-Time
ond-Order M odifiers for Real-Time Optimization
Optimization: Status, Issues, and Opportunities
[J]. IFAC Proceedings Volumes,2014,47( 3) : 7622
[C]/ / Sunggyu Lee. Encyclopedia of Chemical Pro-
- 7628.
cessing,[S. l. ]: CRC Press,2005: 2585 - 2598.
[23] COSTELLO S,FRAN C,OIS G,SRINIVASAN B,et
[12] ENGELL S. Feedback Control for Optimal Process
al. M odifier Adaptation for Run-to-Run Optimization
Operation[J]. Journal of Process Control,2007,17
of Transient Processes[J]. IFAC Proceedings Vol-
( 3) : 203 - 219.
umes,2011,44( 1) : 11471 - 11476.
[13] NARASIM HAN S,JORDACHE C. Data Reconcilia-
[24] SKOGESTAD S. PLANTWIDE CONTROL:
the
tion & Gross Error Detection: an Intelligent Use of
Search for
the Self-Optimizing Control Structure
Process Data[M]. Houston: Gulf Publishing Com-
[J]. Journal of Process Control,2000,10( 5) : 487 -
pany,2000: 75 - 133.
507.
[14] CHACHUAT B,SRINIVASAN B,BONVIN D. Ad-
[25] ARIYUR K B,KRSTIC M . Real-Time Optimization
aptation Strategies for Real-time Optimization[J].
by Extremum-Seeking Control[J]. Lecture Notes in
Computer & Chemical Engineering,2009,33 ( 10) :
Control & Information Sciences,2003,5 ( 2 ) : xii,
1557 - 1567.
236.
[15] FORBES J F,M ARLIN T E. M odel Accuracy for E-
[26] GUAY M ,ZHANG T. Adaptive Extremum Seeking
conomic Optimizing Controllers:
the Bias Update
Control of Nonlinear Dynamic Systems w ith Para-
Case[J]. Ind Eng Chem Res,1994,33( 8) : 1919 -
metric Uncertainties[J]. Automatica,2003,39 ( 7) :
1929.
1283 - 1294.
[16] FORBES J F,M ARLIN T E. Design Cost: A Sys-
[27] FRANCOIS G,SRINIVASAN B,BONVIN D. Use
tematic Approach to Technology Selection for M od-
el-Based Real-Time Optimization Systems [J].
Computer and Chemical Engineering,1996,20 ( 6 /
7) : 717 - 734.
of M easurements for Enforcing the Necessary Con-
ditions of Optimality in the Presence of Constraints
and Uncertainty [J]. Journal of Process Control,
2005,15( 6) : 701 - 712.
[17] FORBES J F,M ARLIN T E,M ACGREGOR J F.
[28] CHACHUAT B,M ARCHETTI A,BONVIN D.
M odel Adequacy Requirements for Optimizing Plant
Operations[J]. Computers and Chemical Engineer-
Process Optimization Via Constraints Adaptation
[J]. Journal of Process Control,2008,18( 3 /4) : 244
ing,1994,18( 6) : 497 - 510.
- 257.
[18] M ARCHETTI A,CHACHUAT B,BONVIN D.
[29] SRINIVASAN B,BIEGLER L T,BONVIN D.
M odifier-Adaptation M ethodology
for Real-Time
Optimization[J]. Ind Eng Chem Res,2009,48
( 13) : 6022 - 6033.
[19] ZHANG Y,FORBES J F. Performance Analysis of
Perturbation-Based M ethods for Real-Time Optimi-
zation[J]. The Canadian Journal of Chemical Engi-
neering,2006,84: 209 - 218.
[20] GROS S,SRINIVASAN B,BONVIN D. Optimizing
Control Based on Output Feedback[J]. Comput.
Chem. Eng. ,2009,33( 1) : 191 - 198.
[21] COSTELLO S,FRANCOIS G,BONVIN D. A Di-
rectional M odifier-Adaptation Algorithm for Real-
Tracking the Necessary Conditions of Optimality
w ith Changing Set of Active Constraints Using a
Barrier-Penalty Function[J]. Computers & Chemical
Engineering,2008,32( 3) : 572 - 579.
[30] HUANG B,KADALI R. Dynamic M odeling,Predic-
tive Control and Performance M onitoring[M]. Lon-
don: Springer,2008: 145 - 175.
[31] DARBY M L,NIKOLAOU M ,JONES J,et al.
RTO: An Overview and Assessment of Current Prac-
tice[J]. Journal of Process Control,2011,21 ( 6 ) :
874 - 884.
( 下转第 90 页)