发明名称 一种炼油工业过程的部分解耦非最小化模型预测控制方法
摘要 本发明涉及一种炼油工业过程的部分解耦非最小化模型预测控制方法。传统的控制手段精度不高,导致后续生产控制参数不稳定、产品合格率和装置效率较低。本发明方法首先基于炼油工业过程模型建立部分解耦状态空间模型,挖掘出基本的过程特性;然后基于该部分解耦状态空间模型建立预测控制回路;最后通过计算预测控制器的参数,将过程对象整体实施预测控制。本发明方法弥补了传统控制的不足,方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。本发明提出的控制技术可以有效减少工艺参数与实际工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的工艺参数达到严格控制。
申请公布号 CN102520617A 申请公布日期 2012.06.27
申请号 CN201110454456.1 申请日期 2011.12.30
申请人 杭州电子科技大学 发明人 张日东;薛安克;陈云;杨成忠;彭冬亮
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.一种炼油工业过程的部分解耦非最小化模型预测控制方法,其特征在于该方法包括以下步骤:(1)利用炼油工业过程模型建立部分解耦状态空间模型,具体方法是:首先采集炼油工业过程的输入输出数据,利用该数据建立输入输出模型如下:<img file="2011104544561100001DEST_PATH_IMAGE002.GIF" wi="57" he="21" />其中<img file="2011104544561100001DEST_PATH_IMAGE004.GIF" wi="18" he="21" />、<img file="2011104544561100001DEST_PATH_IMAGE006.GIF" wi="16" he="18" />、<img file="2011104544561100001DEST_PATH_IMAGE008.GIF" wi="20" he="21" />为三个变量,分别是:<img file="2011104544561100001DEST_PATH_IMAGE010.GIF" wi="381" he="206" />,<img file="2011104544561100001DEST_PATH_IMAGE012.GIF" wi="53" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE014.GIF" wi="54" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE016.GIF" wi="20" he="12" />,<img file="2011104544561100001DEST_PATH_IMAGE018.GIF" wi="60" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE020.GIF" wi="58" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE022.GIF" wi="58" he="26" />,<img file="62968DEST_PATH_IMAGE016.GIF" wi="20" he="12" />,<img file="2011104544561100001DEST_PATH_IMAGE024.GIF" wi="64" he="26" />表示过程的多项式方程,<img file="2011104544561100001DEST_PATH_IMAGE026.GIF" wi="173" he="25" />分别为输入、输出数据,所述的输入输出数据为数据采集器中存储的数据;进一步将上述方程通过克莱姆方程处理为<img file="2011104544561100001DEST_PATH_IMAGE028.GIF" wi="70" he="42" />其中,<img file="2011104544561100001DEST_PATH_IMAGE030.GIF" wi="18" he="18" />是<img file="446545DEST_PATH_IMAGE004.GIF" wi="18" he="21" />的行列式数值,<img file="2011104544561100001DEST_PATH_IMAGE032.GIF" wi="21" he="25" />是将<img file="892219DEST_PATH_IMAGE030.GIF" wi="18" he="18" />的第<img file="2011104544561100001DEST_PATH_IMAGE034.GIF" wi="10" he="18" />列替换成<img file="629231DEST_PATH_IMAGE008.GIF" wi="20" he="21" />获得的行列式数值;将上述过程模型展开得到:<img file="2011104544561100001DEST_PATH_IMAGE036.GIF" wi="166" he="25" />其中,<img file="2011104544561100001DEST_PATH_IMAGE038.GIF" wi="14" he="16" />是得到的模型阶次,<img file="2011104544561100001DEST_PATH_IMAGE040.GIF" wi="109" he="25" />和<img file="2011104544561100001DEST_PATH_IMAGE042.GIF" wi="14" he="18" />为对角矩阵,<img file="2011104544561100001DEST_PATH_IMAGE044.GIF" wi="205" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE046.GIF" wi="201" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE048.GIF" wi="241" he="52" />将过程模型通过后移算子<img file="2011104544561100001DEST_PATH_IMAGE050.GIF" wi="16" he="17" />处理成过程的状态空间表示方式:<img file="2011104544561100001DEST_PATH_IMAGE052.GIF" wi="218" he="49" />其中,<img file="2011104544561100001DEST_PATH_IMAGE054.GIF" wi="70" he="25" />、<img file="2011104544561100001DEST_PATH_IMAGE056.GIF" wi="62" he="22" />分别是第<img file="2011104544561100001DEST_PATH_IMAGE058.GIF" wi="34" he="20" />时刻的变量值,<img file="2011104544561100001DEST_PATH_IMAGE060.GIF" wi="529" he="26" />,<img file="2011104544561100001DEST_PATH_IMAGE062.GIF" wi="17" he="18" />为取转置符号;<img file="2011104544561100001DEST_PATH_IMAGE064.GIF" wi="376" he="220" /><img file="2011104544561100001DEST_PATH_IMAGE066.GIF" wi="194" he="26" /><img file="2011104544561100001DEST_PATH_IMAGE068.GIF" wi="185" he="25" /><img file="2011104544561100001DEST_PATH_IMAGE070.GIF" wi="21" he="25" />为一单位矩阵;定义一过程期望输出为<img file="2011104544561100001DEST_PATH_IMAGE072.GIF" wi="33" he="22" />,并且输出误差<img file="2011104544561100001DEST_PATH_IMAGE074.GIF" wi="33" he="22" />为:<img file="2011104544561100001DEST_PATH_IMAGE076.GIF" wi="118" he="22" />进一步得到第<img file="775785DEST_PATH_IMAGE058.GIF" wi="34" he="20" />时刻的输出误差<img file="2011104544561100001DEST_PATH_IMAGE078.GIF" wi="52" he="22" />为:<img file="2011104544561100001DEST_PATH_IMAGE080.GIF" wi="346" he="25" />其中,<img file="2011104544561100001DEST_PATH_IMAGE082.GIF" wi="52" he="22" />为第<img file="831728DEST_PATH_IMAGE058.GIF" wi="34" he="20" />时刻的过程期望输出;     最后定义一个新的复合状态变量:<img file="2011104544561100001DEST_PATH_IMAGE084.GIF" wi="109" he="49" />  将上述处理过程综合为一个部分解耦的过程模型:<img file="2011104544561100001DEST_PATH_IMAGE086.GIF" wi="250" he="22" />其中,<img file="DEST_PATH_IMAGE088.GIF" wi="53" he="22" />为第<img file="745458DEST_PATH_IMAGE058.GIF" wi="34" he="20" />时刻的复合状态变量,并且<img file="DEST_PATH_IMAGE090.GIF" wi="290" he="52" />(2)基于该部分解耦状态空间模型设计预测控制器,具体方法是:a.定义该预测函数控制器的目标函数为:<img file="DEST_PATH_IMAGE092.GIF" wi="354" he="48" />其中,<img file="DEST_PATH_IMAGE094.GIF" wi="17" he="18" />是预测步长,<img file="DEST_PATH_IMAGE096.GIF" wi="22" he="18" />是预测步长,<img file="DEST_PATH_IMAGE098.GIF" wi="21" he="26" />、<img file="DEST_PATH_IMAGE100.GIF" wi="20" he="26" />是加权矩阵,<img file="DEST_PATH_IMAGE102.GIF" wi="56" he="22" />,<img file="DEST_PATH_IMAGE104.GIF" wi="57" he="22" />分别为第<img file="DEST_PATH_IMAGE106.GIF" wi="37" he="22" />时刻的复合变量和输入变量;b.定义控制变量的作用范围为<img file="DEST_PATH_IMAGE108.GIF" wi="138" he="22" />c.计算控制器的参数,具体是:首先定义<img file="DEST_PATH_IMAGE110.GIF" wi="371" he="107" /><img file="DEST_PATH_IMAGE112.GIF" wi="249" he="112" />然后依据下式计算控制向量<img file="DEST_PATH_IMAGE114.GIF" wi="29" he="20" />:<img file="DEST_PATH_IMAGE116.GIF" wi="268" he="25" />其中,<img file="DEST_PATH_IMAGE118.GIF" wi="17" he="22" />,<img file="DEST_PATH_IMAGE120.GIF" wi="16" he="18" />是依据控制要求设定的两个矩阵,<img file="DEST_PATH_IMAGE122.GIF" wi="26" he="18" />是依据过程期望输出设定的输出向量。
地址 310018 浙江省杭州市下沙高教园区2号大街