发明名称 空分节能过程的比例积分控制系统及方法
摘要 一种空分节能过程的比例积分控制系统,包括与空分塔直接连接的现场智能仪表和DCS系统,所述DCS系统包括存储装置、控制站和上位机,所述现场智能仪表与存储装置、控制站和上位机相连,所述的上位机包括用以求解控制律输出操作变量值的比例积分控制器,所述比例积分控制器包括组分推断模块和比例积分控制律求解模块。本发明也提供了一种空分节能过程的比例积分控制方法。本发明能够有效快速的跟踪空分节能过程的设定值变化,针对空分节能过程的控制模型容易产生模型失配的情况,建立结构简单,控制高效的比例积分控制方案。
申请公布号 CN101901005B 申请公布日期 2012.01.25
申请号 CN201010213752.8 申请日期 2010.06.30
申请人 浙江大学 发明人 刘兴高;周叶翔
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 一种空分节能过程的比例积分控制系统,包括与空分塔直接连接的现场智能仪表和DCS系统,所述DCS系统包括存储装置、控制站和上位机,所述现场智能仪表与存储装置、控制站和上位机相连,其特征在于:所述上位机包括用以求解比例积分控制律并输出控制变量值的比例积分控制器,所述的比例积分控制器包括:组分推断模块,用以获取现场智能仪表检测到的温度,压强数据计算空分塔上塔的各塔板处的组分浓度,采用式(1)(2)得到: <mrow> <msub> <mi>X</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>N</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&times;</mo> <msub> <mi>&alpha;</mi> <mi>N</mi> </msub> <mo>&times;</mo> <msup> <mn>10</mn> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mi>N</mi> </msub> </mrow> <msub> <mi>b</mi> <mi>N</mi> </msub> </mfrac> <mo>-</mo> <msub> <mi>a</mi> <mi>N</mi> </msub> <mo>)</mo> </mrow> </msup> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>&alpha;</mi> <mi>N</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>X</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>O</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&times;</mo> <msub> <mi>&alpha;</mi> <mi>O</mi> </msub> <mo>&times;</mo> <msup> <mn>10</mn> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mi>O</mi> </msub> </mrow> <msub> <mi>b</mi> <mi>O</mi> </msub> </mfrac> <mo>-</mo> <msub> <mi>a</mi> <mi>O</mi> </msub> <mo>)</mo> </mrow> </msup> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>&alpha;</mi> <mi>O</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>其中,k为当前采样时刻,Xi,N(k)为k采样时刻空分塔上塔第i块塔板处氮的液相组分浓度,Xi,O(k)为空分塔上塔第i块塔板处氧的液相组分浓度,Pr(k)为k采样时刻上塔压强,Ti(k)为k采样时刻上塔第i块塔板处的温度,αN、αO分别为氮和氧相对于氩的相对挥发度,aN、bN、cN、aO、bO、cO为安东尼常数;比例积分控制律求解模块,用以根据当前氮和氧的液相组分浓度数据和当前时刻操作变量值求取当前的控制变量的理想改变值,控制律采用式(3)(4)得到: <mrow> <msub> <mi>P</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mo>[</mo> <msub> <mi>K</mi> <mn>11</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>X</mi> <mn>1,1</mn> </msub> <mo>*</mo> </msup> <mo>-</mo> <msub> <mi>X</mi> <mn>1,1</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mn>12</mn> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mrow> <mo>(</mo> <msup> <msub> <mi>X</mi> <mn>1,1</mn> </msub> <mo>*</mo> </msup> <mo>-</mo> <msub> <mi>X</mi> <mn>1,1</mn> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>t</mi> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mo>[</mo> <msub> <mi>K</mi> <mn>21</mn> </msub> <mrow> <mo>(</mo> <msup> <msub> <mi>X</mi> <mrow> <mn>3</mn> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mo>*</mo> </msup> <mo>-</mo> <msub> <mi>X</mi> <mrow> <mn>3</mn> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mn>22</mn> </msub> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mrow> <mo>(</mo> <msup> <msub> <mi>X</mi> <mrow> <mn>3</mn> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mo>*</mo> </msup> <mo>-</mo> <msub> <mi>X</mi> <mrow> <mn>3</mn> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>t</mi> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>其中,Pr(k‑1)为k‑1采样时刻上塔压强,k为当前采样时刻,X1,1*为塔顶液氮浓度设定值,X3,n*为塔底液氧浓度设定值,X1,1(k‑1)、X1,1(j)分别为k‑1采样时刻、j时刻的塔顶液氮的浓度测量值,X3,n(k‑1)、X3,n(j)分别为为k‑1采样时刻、j时刻的塔底液氧的浓度测量值,t为采样周期,K11,K12,K21,K22为控制律参数,q(k)为k采样时刻下塔进料热状况,q(k‑1)为k‑1采样时刻下塔进料热状况。
地址 310027 浙江省杭州市西湖区浙大路38号