发明名称 内部热耦合精馏塔的非线性模型预测控制系统及方法
摘要 一种内部热耦合精馏塔的非线性模型预测控制系统,包括与内部热耦合精馏塔直接连接的现场智能仪表、用于存放历史数据的数据库及上位机,所述现场智能仪表、数据库及上位机依次相连,所述的上位机包括用以运用非线性动态优化在线运算得到当前时刻控制变量的输出值的非线性模型预测控制器,所述非线性模型预测控制器包括组分推断模块、非线性模型模型预测控制模块。本发明较传统PID等控制系统,不仅实现了对两端组分的平稳操作,并且具有更好的动态控制效果。
申请公布号 CN101887262A 申请公布日期 2010.11.17
申请号 CN201010214267.2 申请日期 2010.06.30
申请人 浙江大学 发明人 刘兴高;王成裕
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 1.一种内部热耦合精馏塔的非线性模型预测控制系统,包括与内部热耦合精馏塔直接连接的现场智能仪表、用于存放历史数据的数据库及上位机,所述现场智能仪表、数据库及上位机依次相连,其特征在于:所述的上位机包括非线性模型预测控制器,所述的非线性模型预测控制器包括:组分推断模块:用于依据检测得到的温度与压强数据推断组分,其算式为(1)、(2):<maths num="0001"><![CDATA[<math><mrow><msub><mi>Y</mi><mn>1</mn></msub><mo>=</mo><mfrac><mi>&alpha;</mi><mrow><mi>&alpha;</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><msup><mn>10</mn><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mfrac><mi>b</mi><mrow><msub><mi>T</mi><mn>1</mn></msub><mo>+</mo><mi>c</mi></mrow></mfrac><mo>)</mo></mrow></msup><mrow><mrow><mo>(</mo><mi>&alpha;</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>P</mi><mi>r</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>Xn</mi><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mi>s</mi></msub><mi>&alpha;</mi></mrow><mrow><msup><mn>10</mn><mrow><mo>(</mo><mi>a</mi><mo>-</mo><mfrac><mrow><msub><mi>T</mi><mi>n</mi></msub><mo>+</mo><mi>c</mi></mrow><mi>b</mi></mfrac><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>&alpha;</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mrow><mi>&alpha;</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,Y<sub>1</sub>为塔顶产品中的轻产品组分,X<sub>n</sub>为塔底产品中的轻产品组分,P<sub>r</sub>为精馏段压强,P<sub>s</sub>为提馏段压强,T<sub>1</sub>、T<sub>n</sub>分别为塔顶、塔底温度,α为相对挥发度,a、b、c为安托尼常数;非线性模型预测控制模块,用于编写并储存用微分方程描述的内部热耦合精馏塔机理模型,并支持对关键过程参数进行在线调整,求解如下非线性动态优化命题:目标函数:<maths num="0003"><![CDATA[<math><mrow><munder><mi>min</mi><mrow><mi>u</mi><mrow><mo>(</mo><mo>&CenterDot;</mo><mo>)</mo></mrow></mrow></munder><mi>J</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mi>t</mi><mrow><mi>t</mi><mo>+</mo><msub><mi>T</mi><mi>p</mi></msub></mrow></msubsup><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>1</mn></msub><msub><mi>J</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;</mi><mn>2</mn></msub><msub><mi>J</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>z</mi><mo>,</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;</mi><mn>3</mn></msub><msub><mi>J</mi><mn>3</mn></msub><mrow><mo>(</mo><msub><mi>Vu</mi><mi>MV</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mi>d&tau;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>约束条件:<img file="FDA0000022890370000014.GIF" wi="543" he="81" />(3)0=g(x(t),z(t),u(t),p)    (4)u<sup>LB</sup>≤u(τ)≤u<sup>UB</sup>,<maths num="0004"><![CDATA[<math><mrow><mo>&ForAll;</mo><mi>&tau;</mi><mo>&Element;</mo><mo>[</mo><mi>t</mi><mo>,</mo><mi>t</mi><mo>+</mo><msub><mi>T</mi><mi>c</mi></msub><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>u(τ)=u(t+T<sub>c</sub>),<maths num="0005"><![CDATA[<math><mrow><mo>&ForAll;</mo><mi>&tau;</mi><mo>&Element;</mo><mo>[</mo><mi>t</mi><mo>,</mo><msub><mi>T</mi><mi>c</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><msub><mi>T</mi><mi>p</mi></msub><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>x<sup>LB</sup>≤x(τ)≤x<sup>UB</sup>,<maths num="0006"><![CDATA[<math><mrow><mo>&ForAll;</mo><mi>&tau;</mi><mo>&Element;</mo><mo>[</mo><mi>t</mi><mo>,</mo><mi>t</mi><mo>+</mo><msub><mi>T</mi><mi>p</mi></msub><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,T<sub>p</sub>和T<sub>c</sub>分别是预测时域和控制时域,目标值J由三部分组成,第一部分J<sub>1</sub>表征将来预测输出与目标轨迹的偏差值,e为将来预测输出与目标轨迹的偏差值,J<sub>2</sub>表征于过程对象输入输出相关的经济指标,x为微分变量,z为代数变量,u为控制变量,J<sub>3</sub>表征控制变量的变化量Vu<sub>MV</sub>,μ<sub>1</sub>、μ<sub>2</sub>、μ<sub>3</sub>分别为权重系数,t为时间,τ为暂态时间变量,<img file="FDA0000022890370000018.GIF" wi="84" he="73" />为微分变量x的一阶导数,p为过程参数,f表示微分方程组,g表示代数方程组;公式(3)是非线性过程模型,公式(4)是过程的初始状态,公式6)是控制变量的上下限约束,公式(7)假设控制变量的值在大于等于控制时域,小于等于预测时域的区域内保持不变,公式(8)是过程微分变量的上下限,在每个采样时刻,求解上述非线性动态优化命题,得到最优的控制变量的值,实现系统的在线滚动优化;在线滚动优化,运用非线性动态优化在线运算得到当前时刻的控制变量,进料热状况q(k),精馏塔压强Pr(k)的值。
地址 310027 浙江省杭州市西湖区浙大路38号