发明名称 一种误差容忍限机制下的焦化加热炉炉膛压力控制方法
摘要 本发明公开了一种误差容忍限机制下的焦化加热炉炉膛压力控制方法。本发明方法首先基于焦化加热炉炉膛压力实时过程数据建立过程模型,挖掘出基本的过程特性;然后基于该过程模型建立比例积分控制回路;最后通过计算预测函数控制器的参数,将比例积分控制与焦化加热炉炉膛压力对象整体实施预测函数控制。本发明方法弥补了传统控制的不足,并有效地方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。本发明提出的控制技术可以有效减少理想炉膛压力工艺参数与实际炉膛压力工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的炉膛压力工艺参数达到严格控制。
申请公布号 CN102419551B 申请公布日期 2014.05.21
申请号 CN201110373209.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="FDA00003474711000012.GIF" wi="56" he="74" />(k),具体是:<img file="FDA00003474711000013.GIF" wi="465" he="74" />其中,y<sub>L</sub>(∞)是比例积分控制器的输出有阶跃变化时的实际过程输出y<sub>L</sub>(k)的稳态值;b.选取满足<img file="FDA00003474711000014.GIF" wi="580" he="73" />的两个时间计算点k<sub>1</sub>和k<sub>2</sub>,依据下式计算比例积分控制器所需要的增益K、时间常数T、时滞参数τ:K=y<sub>L</sub>(∞)/qT=2(k<sub>1</sub>-k<sub>2</sub>)τ=2k<sub>1</sub>-k<sub>2</sub>其中,q为过程的比例积分控制器输出的阶跃变化幅度;c.计算过程的比例积分控制器的参数,具体是:K<sub>c</sub>=1.2T/KτT<sub>i</sub>=2τ其中K<sub>c</sub>为比例积分控制器的比例参数,T<sub>i</sub>为比例积分控制器的积分参数;(2)设计预测函数比例积分控制器,具体步骤是:d.将过程的比例积分控制器停留在自动操作状态,操作拨盘使其输入有阶跃变化,由记录仪表记录实时过程的输出,将过程输出值y(k)的响应曲线转换成无量纲形式y<sup>*</sup>(k),具体是:y<sup>*</sup>(k)=y(k)/y(∞)其中,y(∞)是过程的比例积分控制器的输入有阶跃变化时的过程输出y(k)的稳态值;e.选取满足y(k<sub>3</sub>)=0.39,y(k<sub>4</sub>)=0.63的另两个时间计算点k<sub>3</sub>和k<sub>4</sub>,依据下式计算预测函数比例积分控制器所需要的增益K<sub>1</sub>、时间常数T<sub>1</sub>、时滞参数τ<sub>1</sub>:K<sub>1</sub>=y(∞)/q<sub>1</sub>T<sub>1</sub>=2(k<sub>3</sub>-k<sub>4</sub>)τ<sub>1</sub>=2k<sub>3</sub>-k<sub>4</sub>其中,q<sub>1</sub>为过程的比例积分控制器输入的阶跃变化幅度;f.将步骤e得到的参数转化为拉普拉斯形式的局部受控传递函数模型:<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><mn>1</mn><mrow><msub><mi>T</mi><mi>m</mi></msub><mi>s</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&tau;</mi><mi>m</mi></msub><mi>s</mi></mrow></msup></mrow></math>]]></maths>其中,s为拉普拉斯变换算子,y<sub>m</sub>(s)表示当前时刻过程模型的输出值的拉普拉斯变换,q<sub>1</sub>(s)表示过程模型的比例积分控制器输入的拉普拉斯变换;T<sub>m</sub>为局部受控传递函数模型的时间常数、τ<sub>m</sub>为局部受控传递函数模型的时滞参数:T<sub>m</sub>=T<sub>1</sub>τ<sub>m</sub>=τ<sub>1</sub>g.依据步骤f计算出的模型参数设计预测函数比例积分控制器,具体步骤是:①对所述的局部受控传递函数模型在采样时间T<sub>s</sub>下加一个零阶保持器离散化,得到离散模型为y<sub>m</sub>(k+1)=a<sub>m</sub>y<sub>m</sub>(k)+(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="FDA00003474711000024.GIF" wi="306" he="79" />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="FDA00003474711000025.GIF" wi="1615" he="236" />其中,y<sub>Pav</sub>(k)是k时刻离散模型的补偿输出,<img file="FDA00003474711000023.GIF" wi="82" he="81" />是控制量计算参数,u(k-1)是k-1时刻离散模型的输入。
地址 310027 浙江省杭州市西湖区浙大路38号