发明名称 废塑料裂解炉余热烘干装置压力控制方法
摘要 本发明公开了一种废塑料裂解炉余热烘干装置压力控制方法。现在的烘干装置压力控制基本上采用传统的开环控制方式,甚至是手动的操作手段,控制参数的调节完全依赖于技术人员的经验。另外,现有基于状态空间的预测控制方法,通常采用无限时间域内渐近收敛的观测器进行系统状态估计,估计误差难以在有限时间域内收敛,因此估计精度和收敛速度有待提高。本发明采用预测控制方法进行生活垃圾废塑料裂解炉余热烘干装置的压力控制,首次设计了有限时间域内收敛的观测器,对系统状态向量进行估计用于压力输出值的预测,并通过基于递归神经网络的优化方法进行优化求解,从而提高余热烘干装置压力控制的精确性和快速性。
申请公布号 CN103345161A 申请公布日期 2013.10.09
申请号 CN201310281209.5 申请日期 2013.07.05
申请人 杭州电子科技大学 发明人 陈云;谭立华;周绍生
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1. 废塑料裂解炉余热烘干装置压力控制方法,其特征在于该方法的具体步骤是:步骤1. 建立系统的预测模型,具体是:首先,以生活垃圾废塑料裂解炉余热烘干装置的高温烟气阀阀门开度值、冷气阀阀门开度值、烘干机转速和废塑料进料量为输入控制量,以高精度数字式压力测量仪采集到的余热烘干装置的压力值为输出量,通过系统辨识方法,建立余热烘干装置压力控制系统的离散时间传递函数模型<img file="DEST_PATH_IMAGE002.GIF" wi="161" he="25" />其中<img file="DEST_PATH_IMAGE004.GIF" wi="34" he="22" />为<i>k</i>时刻由高精度数字式压力测量仪采集所得的烘干气体压力值;<img file="DEST_PATH_IMAGE006.GIF" wi="226" he="26" />为<i>k</i>时刻的控制输入变量,其中<img file="DEST_PATH_IMAGE008.GIF" wi="38" he="25" />表示<i>k</i>时刻高温烟气阀阀门开度值,<img file="DEST_PATH_IMAGE010.GIF" wi="40" he="25" />表示<i>k</i>时刻冷气阀阀门开度值,<img file="DEST_PATH_IMAGE012.GIF" wi="40" he="25" />表示<i>k</i>时刻烘干机转速,<img file="DEST_PATH_IMAGE014.GIF" wi="40" he="25" />表示<i>k</i>时刻废塑料进料量;<img file="DEST_PATH_IMAGE016.GIF" wi="46" he="25" />和<img file="DEST_PATH_IMAGE018.GIF" wi="46" he="25" />表示通过辨识得到的多项式矩阵,其形式为<img file="DEST_PATH_IMAGE020.GIF" wi="134" he="48" />,<img file="DEST_PATH_IMAGE022.GIF" wi="112" he="48" />其中<img file="DEST_PATH_IMAGE024.GIF" wi="18" he="25" />、<img file="DEST_PATH_IMAGE026.GIF" wi="18" he="25" />表示需要辨识的模型参数,<i>I</i>表示具有合适维数的单位矩阵,<img file="DEST_PATH_IMAGE028.GIF" wi="18" he="25" />、<img file="DEST_PATH_IMAGE030.GIF" wi="18" he="25" />表示采样个数;然后,将前述离散时间传递函数模型通过状态空间实现,建立下述基于状态空间描述的烘干装置压力控制系统模型<img file="DEST_PATH_IMAGE032.GIF" wi="190" he="51" />其中<img file="DEST_PATH_IMAGE034.GIF" wi="411" he="33" />为系统<i>k</i>时刻的状态向量,<img file="DEST_PATH_IMAGE036.GIF" wi="149" he="22" />为<i>k</i>时刻控制输入的增量,系数矩阵<i>A</i>、<i>B</i>和<i>C</i>分别为<img file="DEST_PATH_IMAGE038.GIF" wi="461" he="291" />步骤2. 预测输出值,具体是:由步骤1建立的基于状态空间描述的烘干装置压力控制系统模型,通过迭代计算,可得输出值<img file="DEST_PATH_IMAGE040.GIF" wi="304" he="51" />其中<img file="DEST_PATH_IMAGE042.GIF" wi="34" he="22" />为状态向量<img file="DEST_PATH_IMAGE044.GIF" wi="34" he="22" />的估计值,<i>y</i>(<i>k</i>+<i>j</i>) 表示<i>k</i>+<i>j</i>时刻的输出值,<i>j</i>为不小于1的正整数;令<i>N</i>和<i>N</i><sub><i>u</i></sub>分别表示预测时域和控制时域,则可将预测输出的状态空间模型表示为如下形式<img file="DEST_PATH_IMAGE046.GIF" wi="165" he="25" />其中<img file="DEST_PATH_IMAGE048.GIF" wi="296" he="102" /><img file="DEST_PATH_IMAGE050.GIF" wi="350" he="108" />上式中<img file="DEST_PATH_IMAGE052.GIF" wi="69" he="28" />为<i>k</i>+<i>s</i>时刻输出量的估计值,且<img file="DEST_PATH_IMAGE054.GIF" wi="84" he="22" />;由上式知,对输出量的预测依赖于系统状态向量的估计值<img file="298079DEST_PATH_IMAGE042.GIF" wi="34" he="22" />,状态估计的准确性及收敛速度将直接影响输出值<i>y</i>的预测;步骤3.有限时间域状态估计,具体是:对于可观的矩阵对<img file="DEST_PATH_IMAGE056.GIF" wi="54" he="26" />,设计两个观测器<img file="DEST_PATH_IMAGE058.GIF" wi="350" he="28" />其中<i>z</i><sub>1</sub>(<i>k</i>)和<i>z</i><sub>2</sub>(<i>k</i>)分别为<i>k</i>时刻两个观测器的状态向量,<i>L</i><sub>1</sub><i>、L</i><sub>2</sub>为待设计的观测器增益矩阵;引入滞后时间<i>d</i>,<i>d</i>为正整数,并记<img file="DEST_PATH_IMAGE060.GIF" wi="90" he="25" />,则<img file="DEST_PATH_IMAGE062.GIF" wi="240" he="54" />其中<img file="DEST_PATH_IMAGE064.GIF" wi="340" he="102" />则有<img file="DEST_PATH_IMAGE066.GIF" wi="449" he="111" />若初始状态为<img file="DEST_PATH_IMAGE068.GIF" wi="70" he="28" />,及<img file="DEST_PATH_IMAGE070.GIF" wi="159" he="28" />,可得<img file="DEST_PATH_IMAGE072.GIF" wi="310" he="87" />由于<img file="DEST_PATH_IMAGE074.GIF" wi="58" he="20" />和<img file="DEST_PATH_IMAGE076.GIF" wi="74" he="22" />,当<i>k</i>=<i>d</i>时,<img file="DEST_PATH_IMAGE078.GIF" wi="237" he="29" />对于废塑料余热烘干装置压力控制系统,选择绝对值大于零且小于1的实数标量<img file="DEST_PATH_IMAGE080.GIF" wi="108" he="26" />;利用标准的极点配置方法,设计矩阵<i>L</i><sub>1</sub>,将<img file="DEST_PATH_IMAGE082.GIF" wi="26" he="26" />的极点配置为<img file="431733DEST_PATH_IMAGE080.GIF" wi="108" he="26" />,此即<i>M</i><sub>1</sub>的模小于1;时矩阵<i>M</i><sub>1</sub>稳定,则存在非奇异矩阵<img file="DEST_PATH_IMAGE084.GIF" wi="16" he="17" />,将<img file="398814DEST_PATH_IMAGE082.GIF" wi="26" he="26" />做相似变换,得<img file="DEST_PATH_IMAGE086.GIF" wi="294" he="28" />类似的,再设计矩阵<i>L</i><sub>2</sub>,使<img file="DEST_PATH_IMAGE088.GIF" wi="142" he="25" />,其中标量<i>α</i>满足0&lt;<i>α</i>&lt;1,则可将<img file="DEST_PATH_IMAGE090.GIF" wi="26" he="26" />的极点配置为<img file="DEST_PATH_IMAGE092.GIF" wi="136" he="26" />,且有<img file="DEST_PATH_IMAGE094.GIF" wi="316" he="57" />上式表明,对任意有限大小的正整数<i>d</i>,<img file="DEST_PATH_IMAGE096.GIF" wi="238" he="30" />,即矩阵<img file="DEST_PATH_IMAGE098.GIF" wi="78" he="30" />可逆,系统状态向量<img file="676998DEST_PATH_IMAGE044.GIF" wi="34" he="22" />的估计误差将在任意的有限时间步长<i>k</i>=<i>d</i>时收敛至零;步骤4.建立性能指标,具体是:建立二次型优化目标函数<img file="DEST_PATH_IMAGE100.GIF" wi="235" he="28" />其中<i>Q</i>和<i>R</i>为压力预测输出和控制量的加权矩阵;根据实际生活垃圾废塑料裂解炉余热烘干装置压力控制系统特性,建立约束条件<img file="DEST_PATH_IMAGE102.GIF" wi="138" he="54" />其中<img file="DEST_PATH_IMAGE104.GIF" wi="124" he="25" />分别表示实际余热烘干装置控制输入量和压力值的下限和上限;考虑到步骤2建立的系统方程<img file="290382DEST_PATH_IMAGE046.GIF" wi="165" he="25" />,优化问题可表示为<img file="DEST_PATH_IMAGE106.GIF" wi="276" he="75" />其中<img file="DEST_PATH_IMAGE108.GIF" wi="295" he="144" />且<img file="DEST_PATH_IMAGE110.GIF" wi="175" he="30" />,<img file="DEST_PATH_IMAGE112.GIF" wi="154" he="28" />,<img file="DEST_PATH_IMAGE114.GIF" wi="183" he="30" />,<img file="DEST_PATH_IMAGE116.GIF" wi="160" he="28" />;利用基于递归神经网络的优化方法求解上述二次规划问题,提高在线优化求解的速度;步骤5.动态优化求解,具体是:将步骤4建立的二次优化问题利用递归神经网络进行迭代优化求解;首先,建立动态方程<img file="DEST_PATH_IMAGE118.GIF" wi="373" he="63" />其中<i>ξ</i>(<i>t</i>)是递归神经网络的状态向量,<i>γ</i>&gt;0是常数,矩阵<i>H</i>由步骤4给出,<i>P</i><sub><i>ξ</i></sub>的定义如下 <img file="DEST_PATH_IMAGE120.GIF" wi="260" he="108" />然后,利用下述优化方法进行求解:1)系统初始化;令<i>k</i>=1,设定求解时间上限<i>K</i><sub><i>m</i></sub>、预测时域<i>N</i>、控制时域<i>N</i><sub><i>u</i></sub>、采样时间<i>t</i>、加权矩阵<i>Q</i>和<i>R</i>及每一次采样时刻的递归神经网络优化计算时间<i>△t</i>,获得系统模型的状态变量<i>x</i>(<i>k</i>);2)计算递归神经网络参数;由被控系统模型和前述神经网络动态模型计算<i>F</i>、<i>G</i>、<i>H</i>、<i>a</i>和<i>H</i><sup>-1</sup>,计算<i>H</i><sup>-1</sup><i>a</i>和控制变量的下限<img file="DEST_PATH_IMAGE122.GIF" wi="32" he="25" />和上限;3)递归神经网络动态优化;根据上一步计算所得的参数值,递归神经网络开始运行求解,△<i>t</i>时间后停止,获得稳态值<i>U</i>*,即为最优的控制输入;4)预测输出值;将上一步求解的<i>U</i>*的第一个分量,即<img file="DEST_PATH_IMAGE124.GIF" wi="38" he="22" />对应的最优值的第一个分量<img file="DEST_PATH_IMAGE126.GIF" wi="50" he="25" />,作为系统的控制变量,并计算下一步的状态变量<i>x</i>(<i>k</i>+1);5)若<i>k</i>&lt;<i>K</i><sub><i>m</i></sub>,则令<i>k</i>=<i>k</i>+1,返回第2)步循环计算,否则优化过程运行结束。
地址 310018 浙江省杭州市下沙高教园区2号大街