带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf

上传人:b*** 文档编号:4207527 上传时间:2018-09-06 格式:PDF 页数:14 大小:706.31KB
返回 下载 相关 举报
摘要
申请专利号:

CN201710152893.5

申请日:

2017.03.15

公开号:

CN106842955A

公开日:

2017.06.13

当前法律状态:

实审

有效性:

审中

法律详情:

实质审查的生效IPC(主分类):G05B 13/04申请日:20170315|||公开

IPC分类号:

G05B13/04; B01D53/14

主分类号:

G05B13/04

申请人:

东南大学

发明人:

吴啸; 沈炯; 李益国

地址:

211189 江苏省南京市江宁区东南大学路2号

优先权:

专利代理机构:

南京苏高专利商标事务所(普通合伙) 32204

代理人:

李倩

PDF下载: PDF下载
内容摘要

本发明公开了一种带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,该预测控制方法以基于化学吸附的燃烧后CO2捕集系统为被控对象,贫液流量和汽轮机低压缸抽汽量为系统控制输入量,CO2捕集率和再沸器温度为系统输出量,以火电厂烟气流量为扰动信号,基于子空间辨识方法,利用系统运行产生的数据建立预测器,在此基础上建立带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法。本发明方法在保有预测控制所有优点的前提下大幅提高系统主动克服烟气量扰动的能力,从而进一步提高CO2捕集系统对电站变负荷运行的适应性和灵活性。

权利要求书

1.带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,其特征在于:该预测控制方
法以基于化学吸附的燃烧后CO2捕集系统为被控对象,贫液流量和汽轮机低压缸抽汽量为
系统控制输入量,CO2捕集率和再沸器温度为系统输出量,以火电厂烟气流量为扰动信号,
基于子空间辨识方法,利用系统运行产生的数据建立预测器,在此基础上建立带烟气量扰
动抑制的燃烧后CO2捕集系统预测控制方法。
2.根据权利要求1所述的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,其特
征在于,具体包括如下步骤:
步骤1,火电厂变负荷运行,生成的烟气量d处于变动状态,将CO2捕集系统切换到手动状
态,以贫液流量阀门开度ua和汽轮机低压缸抽汽阀门开度信号ub为输入,对CO2捕集系统进
行激励,获取CO2捕集率ya和再沸器温度yb的开环响应数据;
步骤2,选定采样周期Ts,以为扩增输入,为输出,利用
子空间辨识方法,构建带烟气量扰动信号的CO2捕集系统子空间预测矩阵lw,lu;
步骤3,每一采样时刻,通过当下的烟气量,预估在未来一定时间内系统的CO2捕集率
和再沸器温度通过优化计算得到最优的贫液流量阀门开度ua-op和汽轮机低压缸抽汽阀
门开度信号ub-op,将其用于CO2捕集系统。
3.根据权利要求2所述的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,其特
征在于:步骤2中,构建带烟气量扰动信号的CO2捕集系统子空间预测矩阵lw,lu,具体步骤
为:
步骤A:将连续获得的从第0时刻到第2N+j-2时刻的输出数据Y和扩增输入数据分别排
列为Hankel矩阵形式:
<mrow> <mi>Y</mi> <mo>=</mo> <mo>&lsqb;</mo> <mfrac> <msup> <mi>Y</mi> <mi>p</mi> </msup> <msup> <mi>Y</mi> <mi>f</mi> </msup> </mfrac> <mo>&rsqb;</mo> <mo>=</mo> <mrow> <mo>&lsqb;</mo> <mfrac> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mn>0</mn> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mn>2</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mi>N</mi> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mi>N</mi> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mrow> <mn>2</mn> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>y</mi> <mrow> <mn>2</mn> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfrac> <mo>&rsqb;</mo> </mrow> <mo>,</mo> </mrow>
<mrow> <mover> <mi>U</mi> <mo>~</mo> </mover> <mo>=</mo> <mo>&lsqb;</mo> <mfrac> <msup> <mover> <mi>U</mi> <mo>~</mo> </mover> <mi>p</mi> </msup> <msup> <mover> <mi>U</mi> <mo>~</mo> </mover> <mi>f</mi> </msup> </mfrac> <mo>&rsqb;</mo> <mo>=</mo> <mrow> <mo>&lsqb;</mo> <mfrac> <mtable> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mn>0</mn> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mn>1</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mn>2</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mi>j</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mi>N</mi> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> <mtable> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mi>N</mi> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mi>N</mi> <mo>+</mo> <mi>j</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mn>2</mn> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mover> <mi>u</mi> <mo>~</mo> </mover> <mrow> <mn>2</mn> <mi>N</mi> <mo>+</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfrac> <mo>&rsqb;</mo> </mrow> <mo>,</mo> </mrow>
其中,N为矩阵行数,N大于CO2捕集系统阶次,j为矩阵列数,Y和分别表示输出与扩增
输入数据组成的Hankel矩阵,Yf和Yp分别表示输出数据的未来数据和过去数据,和分
别表示扩增输入数据的未来数据和过去数据,yj表示第j个输出数据,表示第j个扩增输
入数据;
步骤B:令对如下矩阵进行QR分解:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msup> <mi>W</mi> <mi>p</mi> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mover> <mi>U</mi> <mo>~</mo> </mover> <mi>f</mi> </msup> </mtd> </mtr> <mtr> <mtd> <msup> <mi>Y</mi> <mi>f</mi> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>R</mi> <mn>11</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mn>22</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>R</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mn>32</mn> </msub> </mtd> <mtd> <msub> <mi>R</mi> <mn>33</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>Q</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Q</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
获得矩阵L:

步骤C:从而获得矩阵Lw=L(:,1:N(m+l)),Lu=L(:,N(m+l)+1:end),m为输入变量维数,
l为输出变量维数,L(:,1:N(m+l))表示矩阵L的前N(m+l)列,L(:,N(m+l)+1:end)表示矩阵L
自第N(m+l)+1列之后的所有列;
步骤D:最后得到子空间预测矩阵lw=Lw(1:l,:),lu=Lu(1:l,1:m)。
4.根据权利要求2所述的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,其特
征在于:步骤3中,采用公式(1)预估在未来一段时间内系统的CO2捕集率和再沸器温度

其中,为CO2捕集系统过去N个时刻的输出和扩增输入数据组合,
为CO2捕集系统过去N个时刻
的扩增输入数据,
为未来N2个时刻的扩增输入数据,
采用如下公式计算性能指标函数J:
其中,Qf和Rf是调节输入输出控制品质的权值矩阵,
<mrow> <msub> <mi>Q</mi> <mi>f</mi> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mi>f</mi> <mi>T</mi> </msubsup> <mo>&gt;</mo> <mn>0</mn> <mo>,</mo> <msub> <mi>R</mi> <mi>f</mi> </msub> <mo>=</mo> <msubsup> <mi>R</mi> <mi>f</mi> <mi>T</mi> </msubsup> <mo>&gt;</mo> <mn>0</mn> <mo>,</mo> </mrow>
rf是未来N1时刻系统CO2捕集率和再沸器温度设定值序列,
<mrow> <msub> <mi>r</mi> <mi>f</mi> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>r</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>r</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msubsup> <mi>r</mi> <mrow> <mi>k</mi> <mo>+</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> <mi>T</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>,</mo> </mrow>
分别表示k+1时刻到k+N1时刻系统CO2捕集率ra和再沸器温度rb设定值,

是未来N1时刻系统CO2捕集率和再沸器温度预估值序列,
<mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> </mtd> <mtd> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>k</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msubsup> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>k</mi> <mo>+</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> </mrow> <mi>T</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>,</mo> </mrow>
分别表示k+1时刻到k+N1时刻系统CO2捕集率ya和再沸器温度yb预估值,

Δuf是未来N2时刻的贫液流量阀门开度信号ua和低压缸抽汽阀门开度信号ub序列
的增量,其中
CO2捕集系统贫液流量阀门和低压缸抽汽阀门开度信号u的幅值约束(umin,umax)和增量
约束(Δumin,Δumax)为:

<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>&Delta;u</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>&le;</mo> <msub> <mi>&Delta;u</mi> <mi>f</mi> </msub> <mo>&le;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mi>I</mi> </mtd> </mtr> </mtable> </mfenced> <msub> <mi>&Delta;u</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
其中,umin,umax分别表示贫液流量阀门和低压缸抽汽阀门开度信号u的最小值与最大
值,Δumin,Δumax分别表示贫液流量阀门和低压缸抽汽阀门开度信号u的最小增量与最大增
量;
每一采样时刻,将公式(1)代入公式(2),并在满足公式(3)和(4)的情况下最小化性能
指标函数J,得到最优的控制增量序列Δuf:
<mrow> <msub> <mi>&Delta;u</mi> <mi>f</mi> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&Delta;u</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>&Delta;u</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <msubsup> <mi>&Delta;u</mi> <mrow> <mi>k</mi> <mo>+</mo> <msub> <mi>N</mi> <mn>2</mn> </msub> </mrow> <mi>T</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>;</mo> </mrow>
提取最优控制增量序列Δuf中的第一块增量Δuk+1,并与当前时刻的控制作用uk相加,
得到算最优的贫液流量阀门和低压缸抽汽阀门开度信号uop=uk+Δ
uk+1。
5.权利要求1所述带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法在火电机组燃
烧后CO2捕集系统方面的应用。

说明书

带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法

技术领域

本发明涉及一种带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,属于预
测控制方法技术领域。

背景技术

随着温室效应及相关气候生态问题的日益严峻,减排CO2已成为国际社会应对气
候变化的关键举措。作为电力供应的主要设备,火电机组是CO2最稳定、最集中的排放源,世
界30%-40%,我国40%~50%的CO2排放来自于火电机组。在积极发展新能源技术、努力提
高火电机组发电效率的同时,火电机组CO2捕集被众多权威机构公认为未来30年内实现大
规模CO2减排最直接有效的技术手段。

在现有火电机组CO2捕集技术中,基于化学吸收法的燃烧后CO2捕集技术直接从电
厂燃烧后的烟气中分离CO2,具有对已有机组优秀的继承性和较好的技术适用性,是当前CO2
捕集电站采用的主流技术。由于CO2的吸附和解析涉及到一系列化学反应,燃烧后CO2捕集系
统存在很大的惯性,导致传统控制方法往往难以取得满意的控制效果。近年来预测控制算
法在CO2捕集应用中取得一定成效,然而随着火电站变负荷运行参与电网调峰日益频繁,CO2
捕集系统常规预测控制算法在应对烟气流量等扰动时,由于缺乏对烟气信号的利用,控制
效果并不理想。目前CO2捕集系统预测控制器单纯考虑解决系统约束、惯性和耦合问题,未
考虑利用可获取的烟气信号主动抗干扰。因此一种在燃烧后CO2捕集系统中加入对烟气信
号利用的预测控制算法的开发很有必要。

发明内容

发明目的:本发明所要解决的技术问题是提供一种带烟气量扰动抑制的燃烧后
CO2捕集系统预测控制方法,该预测控制方法通过利用烟气量信号,提高了CO2捕集系统的调
节品质,改善了其对火电机组变负荷运行的适应性。

本发明还要解决的技术问题是提供上述带烟气量扰动抑制的燃烧后CO2捕集系统
预测控制方法在火电机组燃烧后CO2捕集系统方面的应用。

发明内容:为解决上述技术问题,本发明所采用的技术手段为:

带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,该预测控制方法以基于
化学吸附的燃烧后CO2捕集系统为被控对象,贫液流量和汽轮机低压缸抽汽量为系统控制
输入量,CO2捕集率和再沸器温度为系统输出量,以火电厂烟气流量为扰动信号,基于子空
间辨识方法,利用系统运行产生的数据建立预测器,在此基础上建立带烟气量扰动抑制的
燃烧后CO2捕集系统预测控制方法。

上述带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,具体包括如下步骤:

步骤1,(离线准备)火电厂变负荷运行,生成的烟气量d处于变动状态,将CO2捕集
系统切换到手动状态,以贫液流量阀门开度ua和汽轮机低压缸抽汽阀门开度信号ub为输入,
对CO2捕集系统进行激励,获取CO2捕集率ya和再沸器温度yb的开环响应数据;

步骤2,(离线准备)选定采样周期Ts,以为扩增输入,
为输出,利用子空间辨识方法,构建带烟气量扰动信号的CO2捕集系统子空间
预测矩阵lw,lu;

步骤3,(在线运行)每一采样时刻,通过当下的烟气量,预估在未来一定时间内系
统的CO2捕集率和再沸器温度通过优化计算得到最优的贫液流量阀门开度ua-op和汽轮
机低压缸抽汽阀门开度信号ub-op,将其用于CO2捕集系统。

其中,步骤2中,构建带烟气量扰动信号的CO2捕集系统子空间预测矩阵lw,lu,具体
步骤为:

步骤A:将连续获得的从第0时刻到第2N+j-2时刻的输出数据Y和扩增输入数据分
别排列为Hankel矩阵形式:



其中,N为矩阵行数,N大于CO2捕集系统阶次,j为矩阵列数,Y和分别表示输出与
扩增输入数据组成的Hankel矩阵,Yf和Yp分别表示输出数据的未来数据和过去数据,和
分别表示扩增输入数据的未来数据和过去数据,yi表示第j个输出数据,表示第j个扩
增输入数据;

步骤B:令对如下矩阵进行QR分解:


获得矩阵L:


步骤C:从而获得矩阵Lw=L(:,1:N(m+l)),Lu=L(:,N(m+l)+1:end),m为输入变量
维数,l为输出变量维数,L(:,1:N(m+l))表示矩阵L的前N(m+l)列,L(:,N(m+l)+1:end)表示
矩阵L自第N(m+l)+1列之后的所有列;

步骤D:最后得到子空间预测矩阵lw=Lw(1:l,:),lu=Lu(1:l,1:m)。

其中,步骤3中,采用如下公式预估在未来一段时间内系统的CO2捕集率和再沸器
温度

其中,为CO2捕集系统过去N个时刻的输出和扩增输入数据组合,

为CO2捕集系统过去N
个时刻的扩增输入数据,

为未来N2个时刻的扩增输入数据,

采用如下公式计算性能指标函数J:

其中,Qf和Rf是调节输入输出控制品质的权值矩阵,


rf是未来N1时刻系统CO2捕集率和再沸器温度设定值序列,


分别表示k+1时刻到k+N1时刻系统CO2捕集率ra和再沸器温度rb设定值,

是未来N1时刻系统CO2捕集率和再沸器温度预估值序列,


分别表示k+1时刻到k+N1时刻系统CO2捕集率ya和再沸器温度yb预估
值,

Δuf是未来N2时刻的贫液流量阀门开度信号ua和低压缸抽汽阀门开度信号ub序列
的增量,其中

CO2捕集系统贫液流量阀门和低压缸抽汽阀门开度信号u的幅值约束(umin,umax)和
增量约束(Δumin,Δumax)为:



其中,umin,umax分别表示贫液流量阀门和低压缸抽汽阀门开度信号u的最小值与最
大值,Δumin,Δumax分别表示贫液流量阀门和低压缸抽汽阀门开度信号u的最小增量与最大
增量;

每一采样时刻,将公式(1)代入公式(2),并在满足公式(3)和(4)的情况下最小化
性能指标函数J,得到最优的控制增量序列Δuf:


提取最优控制增量序列Δuf中的第一块增量Δuk+1,并与当前时刻的控制作用uk相
加,得到算最优的贫液流量阀门和低压缸抽汽阀门开度信号

uop=uk+Δuk+1。

上述带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法在火电机组燃烧后CO2
捕集系统方面的应用。

其中,步骤2中,采样周期Ts可以用经验规则T95/Ts=5~15来选取,其中T95为过渡
过程上升到95%的调节时间;步骤3中,预测控制参数Qf、Rf、N1、N2可以根据实际控制过程中
的性能品质、计算时间等因素人为选取。

本发明带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法充分利用火电厂生
成烟气流量信号,运用预测控制思想,每一步进行一次优化求解获得系统运行最佳输入信
号:贫液流量阀门开度和低压缸抽汽阀门开度。仿真结果表明,本发明中的算法相比一般的
预测控制算法能够更加有效地抑制烟气量扰动的影响,维持CO2捕集率和再沸器温度在设
定值附近。

相比于现有技术,本发明技术方案具有的有益效果为:

本发明的预测控制方法具有良好的抗烟气量干扰能力,能够保证预测控制最优性
能的前提下,迅速消除烟气扰动对系统的影响,应用于火电站燃烧后CO2捕集系统能够有效
抑制烟气量干扰,确保CO2捕集率和再沸器温度在设定值附近;当无烟气量扰动时,本发明
比普通预测控制算法具有更好的设定值跟踪和调节性能;此外,由于本发明预测控制完全
基于数据,从而可以有效避免普通预测控制繁琐的建模过程和建模误差的影响;最后,本发
明预测控制方法总体提高了燃烧后CO2捕集系统的运行水平和对电厂负荷变动的适应性。

附图说明

图1为本发明预测控制方法的原理流程图;

图2为本发明预测控制(实线)与常规比例积分微分控制(虚线)在CO2捕集率设定
值阶跃变化下的控制效果对比图(点划线为设定值);

图3为本发明预测控制(实线)与一般预测控制(虚线)在烟气流量变化下的控制效
果对比图(点划线为设定值)。

具体实施方式

根据下述实施例,可以更好地理解本发明。然而,本领域的技术人员容易理解,实
施例所描述的内容仅用于说明本发明,而不应当也不会限制权利要求书中所详细描述的本
发明。

将本发明的预测控制方法应用于某1MW火电机组燃烧后CO2捕集系统系统仿真模
型中,控制的目标是在满足输入约束的条件下,使CO2捕集率和再沸器温度跟踪设定值。

本发明的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,基于子空间辨
识方法,利用系统输入输出和烟气量数据,预估系统在未来一段时间内的性能,在不破坏预
测控制原有设计框架的前提下主动抑制烟气量扰动作用,算法在无扰情况下,与普通预测
控制器具有更优的设定值跟踪和调节能力,本发明总体提高了CO2捕集系统的控制品质,增
强了其对火电厂负荷变动的适应性。

如图1所示,本发明的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,具
体包括如下步骤:

步骤1,在火电厂变负荷运行,生成烟气量d变动的状态下,设计30秒变化一次,持
续30000秒的贫液流量阀门开度信号ua和汽机低压缸抽汽阀门开度信号ub,对系统进行激
励,获取一系列CO2捕集率ya和再沸器温度yb的开环响应数据;

步骤2,选定采样周期Ts=30s,以为扩增输入,为
输出,利用子空间辨识方法,构建带烟气量扰动信号的CO2捕集系统子空间预测矩阵lw,lu,
具体步骤为:

A:将连续获得的1000组输出数据Y和扩增输入数据分别排列为Hankel矩阵形式
(2N+j-2=1000):




其中,N为矩阵行数,取N=10;j为矩阵列数,在硬件条件允许的情况下越大越好,Y
和分别表示输出与扩增输入数据组成的Hankel矩阵,Yf和Yp分别表示输出数据的未来数
据和过去数据,和分别表示扩增输入数据的未来数据和过去数据,yj表示第j个输出
数据,表示第j个扩增输入数据,上标f和p分别表示未来和过去,下标0,1,...,2N+j-2表
示数据的个数;

B;令对如下矩阵进行QR分解:


获得矩阵L,


C:获得矩阵Lw=L(:,1:N(m+l)),Lu=L(:,N(m+l)+1:end),m=3,m为扩增输入变
量维数,l=2,l为输入输出变量维数,L(:,1:N(m+l))表示L的前N(m+l)列,L(:,N(m+l)+1:
end)表示L自第N(m+l)+1列之后的所有列;

D:子空间矩阵lw=Lw(1:l,:),lu=Lu(1:l,1:m);

步骤3,每一采样时刻,预估在当前烟气量下,在未来一定时间内系统的CO2捕集率
和再沸器温度


其中,为CO2捕集系统过去N个时刻的输出和扩增输入数据组合,

为CO2捕集系统过去N个
时刻的扩增输入数据,为未来N2个时刻的扩增输入数据,本例中取N2=10;

步骤4,贫液流量阀门开度ua-op和汽机低压缸抽汽阀门开度信号ub-op,取式(2)性能
指标函数式:


其中,是调节输入输出控制品质的权值矩阵,
rf是未来N1时刻系统CO2捕集率和再沸器温度设定值序列,


分别表示k+1时刻到k+N1时刻系统CO2捕集率ra和再沸器温度rb设定值,

是未来N1时刻系统CO2捕集率和再沸器温度预估值序列,


分别表示k+1时刻到k+N1时刻系统CO2捕集率ya和再沸器温度yb预估
值,可由式(1)描述,取N1=10;Δuf是未来N2时刻的贫液流量阀门开度信号
和低压缸抽汽阀门开度信号序列的增量,其中

考虑CO2捕集系统阀门开度信号的幅值约束(umin=[0 0]T,umax=[0.8 0.08]T)和
增量约束(Δumin=[-0.007/s -0.006/s]T,Δumax=[0.007/s 0.006/s]T):



每一采样时刻,将(1)代入性能指标式(2),并在满足约束(3)和(4)的情况下最小
化(2),得到最优的控制增量序列输入增量提取最优
控制增量序列Δuf计中的第一块Δuk+1,并与当前时刻的控制作用uk相加,计算最优的贫液
流量阀门和低压缸抽汽阀门开度信号

uop=uk+Δuk+1 (5)

并将其施加于燃烧后CO2捕集系统;

步骤5,固定用于预估系统未来输出的子空间矩阵lw和lu,重复步骤3~4以实现连
续控制。

本实施例为了比较本发明中的带烟气量扰动抑制的燃烧后CO2捕集系统预测控制
方法、常规比例积分微分控制方法和一般预测控制方法的控制效果,做了两组仿真试验:仿
真实验1,CO2捕集系统初始捕集率稳定于0.65,在t=600s和3600s,CO2捕集率设定值从0.65
分别变化为0.85和0.75,再沸器温度设定值保持在384K不变;仿真实验2,CO2捕集系统初始
捕集率稳定于0.85,在t=50s,火电厂生成烟气流量从0.13kg/s变化为0.16kg/s。

如图2、3所示,在无烟气量变化时,CO2捕集率设定值阶跃增加或减小情况下,本发
明对燃烧后CO2捕集系统的优化控制效果曲线明显优于常规比例加积分控制器,具有满意
的设定值跟踪和调节能力。当烟气量扰动发生时,本发明的优化控制方法可以消除扰动的
影响,将CO2捕集率和再沸器温度维持在设定值上,同时相比常规扰动抑制预测控制,具有
更快速平稳的扰动抑制效果且无稳态偏差,提高了CO2捕集系统的运行品质以及对火电厂
变负荷运行的适应性。

本发明烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,利用可测得的火电
厂生成烟气流量,建立可以反映烟气量对捕集系统运行影响的高品质预估器,并在此基础
上设计预测控制器,在保有预测控制所有优点的前提下大幅提高系统主动克服烟气量扰动
的能力,从而进一步提高CO2捕集系统对电站变负荷运行的适应性和灵活性。

显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的
实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其
它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发
明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。

带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf_第1页
第1页 / 共14页
带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf_第2页
第2页 / 共14页
带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf_第3页
第3页 / 共14页
点击查看更多>>
资源描述

《带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf》由会员分享,可在线阅读,更多相关《带烟气量扰动抑制的燃烧后COSUB2/SUB捕集系统预测控制方法.pdf(14页珍藏版)》请在专利查询网上搜索。

本发明公开了一种带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法,该预测控制方法以基于化学吸附的燃烧后CO2捕集系统为被控对象,贫液流量和汽轮机低压缸抽汽量为系统控制输入量,CO2捕集率和再沸器温度为系统输出量,以火电厂烟气流量为扰动信号,基于子空间辨识方法,利用系统运行产生的数据建立预测器,在此基础上建立带烟气量扰动抑制的燃烧后CO2捕集系统预测控制方法。本发明方法在保有预测控制所有优点的前提。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 物理 > 控制;调节


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1