发明名称 一种误差容忍限机制下的焦化加热炉氧含量控制方法
摘要 本发明公开了一种误差容忍限机制下的焦化加热炉氧含量控制方法。本发明方法首先基于焦化加热炉氧含量实时过程数据建立过程模型,挖掘出基本的过程特性;然后基于该过程模型建立预测函数控制回路,实施预测函数控制。本发明方法弥补了传统控制的不足,并有效地方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。本发明提出的控制技术可以有效减少理想氧含量工艺参数与实际氧含量工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的氧含量工艺参数达到严格控制。
申请公布号 CN102436178A 申请公布日期 2012.05.02
申请号 CN201110373142.9 申请日期 2011.11.22
申请人 浙江大学 发明人 张建明;张日东
分类号 G05B13/04(2006.01)I;C10G9/00(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 张法高
主权项 1.一种误差容忍限机制下的焦化加热炉氧含量控制方法,其特征在于它的步骤如下:(1)采用响应曲线法设计氧含量过程的数学模型,具体步骤是:a.将过程的预测函数控制器停留在手动操作状态,操作拨盘使其输出有个阶跃变化,由记录仪表记录实际过程的输出值,将实际过程输出值y<sub>L</sub>(k)的响应曲线转换成无量纲形式<img file="FDA0000110947520000011.GIF" wi="145" he="57" />具体是:<img file="FDA0000110947520000012.GIF" wi="448" he="57" />其中,y<sub>L</sub>(∞)是预测函数控制器的输出有阶跃变化时的实际过程输出y<sub>L</sub>(k)的稳态值;b.选取满足<img file="FDA0000110947520000013.GIF" wi="286" he="57" /><img file="FDA0000110947520000014.GIF" wi="277" he="57" />的两个时间计算点k<sub>1</sub>和k<sub>2</sub>,依据下式计算预测函数控制器所需要的增益K、时间常数T、时滞参数τ:K=y<sub>L</sub>(∞)/q<sub>1</sub>T=2(k<sub>1</sub>-k<sub>2</sub>)τ=2k<sub>1</sub>-k<sub>2</sub>其中,q<sub>1</sub>为过程的预测函数控制器输出的阶跃变化幅度;(2)设计预测函数控制器,具体步骤是:c.将步骤b得到的参数转化为拉普拉斯形式的局部受控传递函数模型:<maths num="0001"><![CDATA[<math><mrow><mfrac><mrow><msub><mi>y</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>q</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><msub><mi>K</mi><mi>m</mi></msub><mrow><msub><mi>T</mi><mi>m</mi></msub><mi>s</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msup><mi>e</mi><mrow><msub><mrow><mo>-</mo><mi>&tau;</mi></mrow><mi>m</mi></msub><mi>s</mi></mrow></msup></mrow></math>]]></maths>其中,s为拉普拉斯变换算子,K<sub>m</sub>为局部受控传递函数模型的比例常数,T<sub>m</sub>为局部受控传递函数模型的时间常数,τ<sub>m</sub>为局部受控传递函数模型的时滞,y<sub>m</sub>(s)表示当前时刻过程模型的输出值的拉普拉斯变换,q<sub>1</sub>(s)表示过程模型的预测函数控制器输出的拉普拉斯变换,K<sub>m</sub>=KT<sub>m</sub>=T;τ<sub>m</sub>=τd.依据步骤c计算出的模型参数设计预测函数控制器,具体步骤是:①对该对象在采样时间T<sub>s</sub>下加一个零阶保持器离散化,得到离散模型为y<sub>m</sub>(k+1)=a<sub>m</sub>y<sub>m</sub>(k)+K<sub>m</sub>(1-a<sub>m</sub>)u(k-L)其中y<sub>m</sub>(k+1),y<sub>m</sub>(k)分别是k+1,k时刻离散模型的输出,u(k-L)是k-L时刻离散模型的输入,a<sub>m</sub>为相应的离散模型的参数,<img file="FDA0000110947520000016.GIF" wi="284" he="62" />L为相应的离散模型的时滞,L=τ<sub>m</sub>/T<sub>s</sub>;②选取预测函数控制器的参考轨迹y<sub>r</sub>(k+H),由下式来表示y<sub>r</sub>(k+H)=c(k+H)-λ<sup>H</sup>[c(k)-y(k)]c(k+H),c(k)分别为k+H,k时刻的参考轨迹,λ<sup>H</sup>为参考轨迹的参数,H为预测函数控制器的预测时间参数;③设定误差容许限β,求取控制误差E并依据步骤②得到的预测函数控制器参考轨迹求取控制量u(k)<maths num="0002"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><mfrac><mn>1</mn><mi>H</mi></mfrac><mo>|</mo><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>H</mi><mo>)</mo></mrow><mo>|</mo></mrow></math>]]></maths><img file="FDA0000110947520000022.GIF" wi="1593" he="219" />其中,H是离散时刻,y<sub>Pav</sub>(k)是k时刻离散模型的补偿输出,<img file="FDA0000110947520000023.GIF" wi="56" he="60" />是控制量计算参数,u(k-1)是k-1时刻离散模型的输入。
地址 310027 浙江省杭州市西湖区浙大路38号