发明名称 一种废塑料裂解反应釜的压力平衡控制方法
摘要 本发明公开了一种废塑料裂解反应釜的压力平衡控制方法。目前废塑料裂解反应釜的压力平衡控制基本上采用标准PID控制等较简单的方式,控制效果不够理想。本发明采用预测控制和动态前馈补偿相结合的复合控制方法进行废塑料裂解反应釜的压力平衡控制,首先在系统运行平稳且没有可测扰动输入的情况下设计预测控制器,然后在系统具有可测扰动的情形下通过串联低通滤波器的方法设计严格正则的动态前馈控制器,最后将上述两类控制器相结合,获得前馈-反馈复合控制律,从而提高废塑料裂解反应压力控制的控制精度和抗干扰性能,有助于裂解反应的高效和安全运行。
申请公布号 CN103336539A 申请公布日期 2013.10.02
申请号 CN201310281206.1 申请日期 2013.07.05
申请人 杭州电子科技大学 发明人 王俊宏;陈云;张乐;郭云飞
分类号 G05D16/20(2006.01)I 主分类号 G05D16/20(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1. 一种废塑料裂解反应釜的压力平衡控制方法,其特征在于该方法的具体步骤是:步骤一:平稳运行、不存在扰动输入情况下的预测控制设计a. 建立系统稳态模型;首先,以生活垃圾废塑料裂解反应釜转速、过剩空气系数、烟气排出速度、进料速度、出料速度和催化剂供应量为控制输入量,以压力测量仪采集到的裂解反应釜内的压力值为输出量,在被控系统平稳运行、不存在扰动输入的情况下,通过系统辨识方法,建立反应釜压力平衡控制的离散时间传递函数模型<img file="2013102812061100001DEST_PATH_IMAGE002.GIF" wi="161" he="25" />其中<img file="2013102812061100001DEST_PATH_IMAGE004.GIF" wi="34" he="22" />为<i>k</i>时刻由压力传感器测量所得的裂解反应釜内的压力值;<img file="2013102812061100001DEST_PATH_IMAGE006.GIF" wi="294" he="26" />为<i>k</i>时刻的控制输入向量,6个控制输入分量依次表示<i>k</i>时刻的反应釜转速、过剩空气系数、烟气排出速度、进料速度、出料速度和催化剂供应量;<img file="2013102812061100001DEST_PATH_IMAGE008.GIF" wi="46" he="25" />和<img file="2013102812061100001DEST_PATH_IMAGE010.GIF" wi="46" he="25" />表示通过辨识得到系统参数<img file="2013102812061100001DEST_PATH_IMAGE012.GIF" wi="134" he="48" />,<img file="2013102812061100001DEST_PATH_IMAGE014.GIF" wi="112" he="48" />其中<img file="2013102812061100001DEST_PATH_IMAGE016.GIF" wi="18" he="25" />、<img file="2013102812061100001DEST_PATH_IMAGE018.GIF" wi="18" he="25" />表示需要辨识的模型参数,<i>I</i>表示具有合适维数的单位矩阵,<img file="2013102812061100001DEST_PATH_IMAGE020.GIF" wi="18" he="25" />、<img file="2013102812061100001DEST_PATH_IMAGE022.GIF" wi="18" he="25" />表示采样个数;然后,可将前述反应釜压力控制系统模型表示为下述状态空间模型<img file="2013102812061100001DEST_PATH_IMAGE024.GIF" wi="190" he="51" />其中<img file="2013102812061100001DEST_PATH_IMAGE026.GIF" wi="375" he="30" />为状态向量,<img file="DEST_PATH_IMAGE028.GIF" wi="149" he="22" />为<i>k</i>时刻控制输入的增量,系数矩阵<i>A</i>、<i>B</i>和<i>C</i>分别为<img file="DEST_PATH_IMAGE030.GIF" wi="461" he="291" />b. 预测输出值;由a步建立的裂解反应釜压力控制系统模型,通过迭代计算,可得<i>k</i>+<i>j</i>时刻的输出值为<img file="DEST_PATH_IMAGE032.GIF" wi="304" he="51" />其中<img file="DEST_PATH_IMAGE034.GIF" wi="34" he="22" />为向量<img file="DEST_PATH_IMAGE036.GIF" wi="34" he="22" />的估计值,<i>j</i>为正整数;进一步,可将预测输出表示为下述形式<img file="DEST_PATH_IMAGE038.GIF" wi="167" he="25" />其中<img file="DEST_PATH_IMAGE040.GIF" wi="353" he="220" />其中<i>N</i>和<i>N</i><sub><i>u</i></sub>分别表示预测时域和控制时域,<img file="DEST_PATH_IMAGE042.GIF" wi="69" he="28" />为<i>k</i>+<i>s</i>时刻输出量的预测值,且<img file="DEST_PATH_IMAGE044.GIF" wi="84" he="22" />;由上式可知,输出量的预测依赖于状态向量<img file="356518DEST_PATH_IMAGE036.GIF" wi="34" he="22" />的估计值<img file="108573DEST_PATH_IMAGE034.GIF" wi="34" he="22" />;可通过设计一个渐近估计器,用于状态向量的估计;利用标准极点配置方法,设计一个维数适当的观测器增益矩阵<i>L</i>,使得矩阵(<i>A</i>-<i>LC</i>)的所有特征值都小于1,则该观测器是渐近收敛的;此时,<img file="931036DEST_PATH_IMAGE036.GIF" wi="34" he="22" />的估计值<img file="617976DEST_PATH_IMAGE034.GIF" wi="34" he="22" />可按下式计算<img file="DEST_PATH_IMAGE046.GIF" wi="301" he="25" />c. 优化控制律求解;建立二次型优化目标函数<img file="DEST_PATH_IMAGE048.GIF" wi="340" he="54" />其中<img file="DEST_PATH_IMAGE050.GIF" wi="36" he="22" />和<img file="DEST_PATH_IMAGE052.GIF" wi="37" he="22" />为标量系数;当取<img file="DEST_PATH_IMAGE054.GIF" wi="57" he="22" />和<img file="DEST_PATH_IMAGE056.GIF" wi="64" he="22" />时(<img file="DEST_PATH_IMAGE058.GIF" wi="17" he="22" />为定常的标量),优化目标函数可写为<img file="DEST_PATH_IMAGE060.GIF" wi="396" he="34" />利用二次规划(QP)问题的求解方法可得最优控制输入为<img file="DEST_PATH_IMAGE062.GIF" wi="238" he="28" />其中<i>I</i>表示维数合适的单位矩阵;在控制实施过程中,只有控制序列<img file="DEST_PATH_IMAGE064.GIF" wi="38" he="22" />的第一个分量,即<img file="DEST_PATH_IMAGE066.GIF" wi="44" he="22" />,被送到实际废塑料裂解反应釜,使其产生相应的控制效果,并在下一采样时刻重复上述优化求解的计算;步骤二:存在可测输入扰动情形下的动态前馈补偿设计在裂解反应过程中,尤其是启停机、运行模式变化、废塑料品质和成分变化、天气变化等诸多原因,都使得系统存在过剩空气系数、供热烟气温度、烟气排出速度、进料速度、出料速度等可测的输入量扰动;在这种情况下,通过设计动态前馈控制器,可以对相应的可测扰动进行动态补偿,提高系统的抗干扰性能;首先,在系统存在可扰动输入的情况下,通过系统辨识方法,建立下述反应釜压力平衡控制的传递函数模型<img file="DEST_PATH_IMAGE068.GIF" wi="294" he="28" />其中<img file="957821DEST_PATH_IMAGE004.GIF" wi="34" he="22" />、<img file="DEST_PATH_IMAGE070.GIF" wi="34" he="22" />分别<i>k</i>时刻的压力测量值和控制输入向量,具体形式前面已给出;<img file="DEST_PATH_IMAGE072.GIF" wi="42" he="25" />为<i>k</i>时刻的可测噪声向量;<img file="990149DEST_PATH_IMAGE008.GIF" wi="46" he="25" />、<img file="983513DEST_PATH_IMAGE010.GIF" wi="46" he="25" />、<img file="DEST_PATH_IMAGE074.GIF" wi="48" he="25" />表示通过辨识得到的系统参数<img file="941105DEST_PATH_IMAGE012.GIF" wi="134" he="48" />,<img file="865067DEST_PATH_IMAGE014.GIF" wi="112" he="48" />,<img file="DEST_PATH_IMAGE076.GIF" wi="134" he="48" />其中<img file="326136DEST_PATH_IMAGE016.GIF" wi="18" he="25" />、<img file="444395DEST_PATH_IMAGE018.GIF" wi="18" he="25" />、<img file="DEST_PATH_IMAGE078.GIF" wi="20" he="25" />表示需要辨识的模型参数,<i>I</i>表示具有合适维数的单位矩阵,<img file="358125DEST_PATH_IMAGE020.GIF" wi="18" he="25" />、<img file="85778DEST_PATH_IMAGE022.GIF" wi="18" he="25" />、<img file="DEST_PATH_IMAGE080.GIF" wi="18" he="25" />表示采样个数;然后,引入被控对象没有扰动时的传递函数<img file="DEST_PATH_IMAGE082.GIF" wi="45" he="25" />和被控对象扰动通道的传递函数<img file="DEST_PATH_IMAGE084.GIF" wi="61" he="27" />,可得<img file="DEST_PATH_IMAGE086.GIF" wi="253" he="28" />其中<img file="DEST_PATH_IMAGE088.GIF" wi="178" he="30" />,<img file="DEST_PATH_IMAGE090.GIF" wi="186" he="30" />;将需要设计的前馈控制增益表示为<img file="DEST_PATH_IMAGE092.GIF" wi="62" he="27" />,则加入前馈补偿环节后,系统的输出<img file="DEST_PATH_IMAGE094.GIF" wi="33" he="24" />可表示为<img file="DEST_PATH_IMAGE096.GIF" wi="394" he="54" />显然,当<img file="DEST_PATH_IMAGE098.GIF" wi="56" he="24" />时,有<img file="DEST_PATH_IMAGE100.GIF" wi="289" he="28" />;若再令<img file="DEST_PATH_IMAGE102.GIF" wi="209" he="27" />,则可设计相应的前馈控制器增益,具体形式为<img file="DEST_PATH_IMAGE104.GIF" wi="206" he="30" />因此,动态前馈控制器的输出为<img file="DEST_PATH_IMAGE106.GIF" wi="344" he="30" />最后,进行正则的、可实现的前馈控制器设计;对于许多实际的控制对象,由于上述被控对象没有扰动时的传递函数<img file="DEST_PATH_IMAGE108.GIF" wi="54" he="28" />是严格正则的,导致其逆<img file="DEST_PATH_IMAGE110.GIF" wi="64" he="30" />是非正则的和不可实现的,而且<img file="227784DEST_PATH_IMAGE110.GIF" wi="64" he="30" />的非正则性将导致控制输入为高频时前馈控制器<img file="DEST_PATH_IMAGE112.GIF" wi="62" he="28" />的幅值太大;为了设计严格正则的动态前馈控制器,将通过串联下述低通滤波器的方法消除<img file="500634DEST_PATH_IMAGE112.GIF" wi="62" he="28" />的高频分量<img file="DEST_PATH_IMAGE114.GIF" wi="408" he="44" />其中<i>T</i>为采样周期,<img file="DEST_PATH_IMAGE116.GIF" wi="98" he="24" />,<img file="DEST_PATH_IMAGE118.GIF" wi="77" he="24" />可依次取为100、120、140、160…等值;步骤三:预测控制-动态前馈补偿的复合控制设计综合步骤一和步骤二设计的预测控制器和动态前馈控制器,采用前馈-反馈复合控制方法对系统进行校正,实现废塑料裂解反应釜压力平衡的有效控制;根据步骤一和步骤二,最终可得系统的预测控制-动态前馈复合控制律为<img file="DEST_PATH_IMAGE120.GIF" wi="553" he="40" />其中<img file="DEST_PATH_IMAGE122.GIF" wi="125" he="25" />。
地址 310018 浙江省杭州市下沙高教园区2号大街