一种微型燃气轮机热电联供系统抗扰预测控制方法
技术领域
本发明涉及微型燃气轮机热电联供系统控制技术，特别是涉及一种微型燃气轮机热电联供系统抗扰预测控制方法。
背景技术
分布式能源系统相较于传统集中供能系统，能在用户端直接向用户提供各种形式能量，实现能量的梯级利用。微型燃气轮机热电联供(MGT-CHP)系统是一种重要的分布式能源形式，具有清洁、高效和结构灵活等特点。然而，由于微型燃气轮机热电联供系统的动态具有大热惯性、多变量强耦合、输入约束、非线性和未知干扰等复杂特性，PID方法已经无法满足性能要求。现有的微型燃气轮机热电联供系统的预测控制系统不能保证闭环系统的稳定性；且现有的微型燃气轮机热电联供系统的控制方法无法在控制器设计阶段考虑控制量约束，且不能同时调节热负荷和电负荷。
发明内容
发明目的：本发明的目的是提供一种微型燃气轮机热电联供系统抗扰预测控制方法，来提高系统抗扰性能，进而提高微型燃气轮机热电联供系统控制品质。
技术方案：本发明的微型燃气轮机热电联供系统抗扰预测控制方法，包括以下步骤：
(1)采集微型燃气轮机热电联供系统的运行数据，包括燃料量、一次水旁路阀开度、转速和二次供水温度；
(2)分别辨识转速系统和温度系统的状态空间模型，然后建立热电联供系统的整体控制模型；
(3)基于热电联供系统整体控制模型建立广义扩增状态观测器；
(4)基于广义扩增状态观测器建立稳定预测控制的预测模型；
(5)建立稳定预测控制器，调节燃料量和一次水旁路阀开度，控制转速和二次供水温度。
进一步的，步骤(2)包括以下步骤：
(21)基于燃料量和转速的数据建立转速系统的过程模型，其形式如下：
其中，x1(k)是在离散k时刻的转速系统的状态量，y1(k)是在离散k时刻的转速，u1(k)是在离散k时刻的燃料量，(A1,B1,C1,a1,b1)表示转速系统的参数；
(22)基于燃料量、一次旁路阀开度和二次供水温度系统建立温度系统的过程模型，其形式如下：
其中，x2(k)是在离散k时刻的温度系统的状态量，y2(k)是在离散k时刻的二次供水温度，u(k)在离散k时刻的微型燃气轮机热电联供系统控制量，u(k)＝[u1(k) u2(k)]T，T表示矩阵转置符号，u2(k)是在离散k时刻的一次水旁路阀开度，(A2,B2,C2,a2,b2)表示温度系统的参数；
(23)建立热电联供系统的整体控制模型，其形式如下：
其中，02×1是2行1列的零矩阵，(A,B,C,a,b)表示热电联供系统控制模型的参数。
进一步的，步骤(3)包括以下步骤：
(31)建立如下的带有扰动项的热电联供系统模型：
其中，x(k)在离散k时刻的微型燃气轮机热电联供系统状态量，y(k)在离散k时刻的微型燃气轮机热电联供系统输出量，y(k)＝[y1(k) y2(k)]T，d(k)是状态集总扰动，v(k)是输出集总扰动，Bd和Cv分别是状态集总扰动矩阵和输出集总扰动矩阵；
(32)建立广义扩增状态观测器，其形式如下：
其中，符号“^”表示估计值，0和I是零矩阵和单位矩阵，L是观测器增益，△表示增量，△s表示s次增量；
(33)求解观测器增益L＝M-1N，矩阵M和N通过求解下面可行性问题得到：
其中，M、N和S是待求的矩阵，S是正定对称矩阵，
进一步的，步骤(4)包括以下步骤：
(41)求解补偿扰动后的超临界火电机组稳态输入ut和状态xt：
其中yr是转速和温度的设定值；
(42)建立稳定预测控制的预测模型，其形式如下：
其中，
进一步的，步骤(5)包括以下步骤：
(51)建立稳定预测控制的无穷时域性能指标，其形式如下：
其中，Q和R分别是状态和控制输入调节矩阵参数；
(52)将步骤(51)的无穷时域性能指标转化为如下的优化问题：
受限于：
其中，
(53)计算离散k时刻的超临界火电机组输入
有益效果：与现有技术相比，本发明具有以下有益效果：
(1)提出的抗扰预测控制方法，将系统的非线性、模型失配和扰动集总一个扰动，利用广义扩增状态观测器对估计集总扰动和状态量，通过前馈补偿得到预测模型，然后设计了稳定预测控制器，在保证系统稳定性的同时，有效地提高了系统抗扰性能；
(2)仿真实验表明，提出的抗扰预测控制方法有效地解决了微型燃气轮机热电联供系统控制存在的控制难点，控制效果优于常规的预测控制和PID控制，具有跟踪速度快、超调量小和抗模型失配和干扰能力强的优点。
附图说明
图1为本发明方法步骤流程图；
图2为本发明实施例中微型燃气轮机热电联供系统示意图；
图3(a)为本发明方法在实验1中微型燃气轮机热电联供系统转速控制结果；
图3(b)为本发明方法在实验1中微型燃气轮机热电联供系统温度控制结果；
图3(c)为本发明方法在实验1中燃料量变化曲线；
图3(d)为本发明方法在实验1中阀门变化曲线；
图4(a)为本发明方法在实验2中微型燃气轮机热电联供系统转速控制结果；
图4(b)为本发明方法在实验2中微型燃气轮机热电联供系统温度控制结果；
图5为本发明方法在实验3中微型燃气轮机热电联供系统转速和温度控制结果；
图6(a)为本发明方法在实验4中微型燃气轮机热电联供系统转速控制结果；
图6(b)为本发明方法在实验4中微型燃气轮机热电联供系统温度控制结果。
具体实施方式
下面结合附图和具体实施例对本发明技术方案做进一步的说明。
微型燃气轮机热电联供系统是一个二输入二输出系统，它的输入量包括燃料量u1和一次水旁路阀开度u2，输出量包括微型燃气轮机转速y1和二次供水温度y2。微型燃气轮机转速只受燃料量影响，二次供水温度同时受到燃料量和一次水旁路阀开度影响，因此根据它的特定物理结构，可拆分成转速系统和温度系统。
如图1所示，本发明一种微型燃气轮机热电联供系统抗扰预测控制方法，包括以下步骤：
(1)采集微型燃气轮机热电联供系统的运行数据，包括燃料量、一次水旁路阀开度、转速和二次供水温度；
(2)分别辨识转速系统和温度系统的状态空间模型，然后建立微型燃气轮机热电联供系统的整体控制模型；
(21)基于燃料量和转速的数据建立转速系统的过程模型，其形式如下：
其中，x1(k)是在离散k时刻的转速系统的状态量，y1(k)是在离散k时刻的转速，u1(k)是在离散k时刻的燃料量，(A1,B1,C1,a1,b1)表示转速系统的参数；
(22)基于燃料量、一次旁路阀开度和二次供水温度系统建立温度系统的过程模型，其形式如下：
其中，x2(k)是在离散k时刻的温度系统的状态量，y2(k)是在离散k时刻的二次供水温度，u(k)在离散k时刻的微型燃气轮机热电联供系统控制量，u(k)＝[u1(k) u2(k)]T，T表示矩阵转置符号，u2(k)是在离散k时刻的一次水旁路阀开度，(A2,B2,C2,a2,b2)表示温度系统的参数；
(23)建立微型燃气轮机热电联供系统的整体控制模型，其形式如下：
其中，02×1是2行1列的零矩阵，(A,B,C,a,b)表示热电联供系统控制模型的参数。
(3)基于微型燃气轮机热电联供系统整体控制模型建立广义扩增状态观测器；
(31)建立如下的带有扰动项的热电联供系统模型：
其中，x(k)在离散k时刻的微型燃气轮机热电联供系统状态量，y(k)在离散k时刻的微型燃气轮机热电联供系统输出量，y(k)＝[y1(k) y2(k)]T，d(k)是状态集总扰动，v(k)是输出集总扰动，Bd和Cv分别是状态集总扰动矩阵和输出集总扰动矩阵；
(32)建立广义扩增状态观测器，其形式如下：
其中，符号“^”表示估计值，0和I是零矩阵和单位矩阵，L是观测器增益，△表示增量，△s表示s次增量；
(33)求解观测器增益L＝M-1N，矩阵M和N通过求解下面可行性问题得到：
其中，M、N和S是待求的矩阵，S是正定对称矩阵，
(4)基于广义扩增状态观测器建立稳定预测控制的预测模型；
(41)求解补偿扰动后的超临界火电机组稳态输入ut和状态xt：
其中yr是转速和温度的设定值；
(42)建立稳定预测控制的预测模型，其形式如下：
其中，
(5)基于步骤(4)得到的预测模型建立稳定预测控制器，调节燃料量和一次水旁路阀开度，控制转速和二次供水温度。
(51)建立稳定预测控制的无穷时域性能指标，其形式如下：
其中，Q和R分别是状态和控制输入调节矩阵参数；
(52)将步骤(51)的无穷时域性能指标转化为如下的优化问题：
受限于：
其中，
(5-3)计算离散k时刻的超临界火电机组输入
实施例。
为了使本发明要解决的技术问题、技术方案和有益效果更加清楚，下面结合附图及具体的实施例进行详尽的描述。
在此以某80kw微型燃气轮机和115kw供热系统组成的微型燃气轮机热电联供系统为例进行仿真实验，其结构图如图2所示，并具体说明本发明的设计方法及实施方案。图2所示的微型燃气轮机热电联供系统由微型燃气轮机和供热系统组成。微型燃气轮机的主要部件包括压气机、燃烧室、透平、回热器、发电机等。微型燃气轮机的工作过程如下：压气机从外界吸入空气，压缩后经过回热器再送入燃烧室，同时天然气燃料也被喷入燃烧室与高温的压缩空气混合，在定压下燃烧；生成的高温高压烟气进入透平做功，推动叶轮带动压气机叶轮一同旋转，同时推动发电机转子旋转产生电能；做功后的乏气通过回热器，再经风机排出。供热系统的主要部件为管壳式换热器以及水-水板式热交换器。微型燃气轮机排出的气体经过烟气换热器与来自冷却塔的一次水换热之后排空；高温的一次水通过一次网管道进入水-水板式热交换器，加热二次水产生热量，自身降温后重新回流进冷却塔完成循环。系统的操作变量是输入燃料量，输出变量是热水温度。
采集微型燃气轮机热电联供系统的运行数据，包括燃料量、一次水旁路阀开度、转速和二次供水温度，利用n4sid算法分别辨识得到转速系统和温度系统的状态空间模型如下：
其中，u＝[uF uV]T，
提出的抗扰预测控制器的设计参数为：采样时间，扰动矩阵Bd＝diag{1,1,1,1,1}，状态加权系数Q＝diag{200000,100,5,1,1}，输入加权系数R＝diag{2,1000}，自由控制步数N＝2，w＝[0.01；0.01；0.01；0.01；0.01]，输入约束umin＝[0；0]，umax＝[0.009；100]，△u1,max＝0.0003，△u1,min＝-0.0003。
为了验证提出方法的优越性，同时将模型预测控制(MPC)和PID控制实施到图2所示的微型燃气轮机热电联供系统。
实验1：大范围设定值跟踪控制，转速和温度的设定值变化设置如下：
实验2：抗阶跃干扰实验。假设系统运行在(51400r/min,82.3℃)操作点，如下的未知输入扰动被加到系统的输入端：
实验3：抗斜波干扰实验。假设系统运行在(63652r/min,82℃)操作点，如下的未知输出扰动被加到系统的输出端：
实验4：抗模型失配实验。这边假设辨识的系统发生如下三种类型的失配。
失配1：
失配2：
失配3：
实验1的仿真结果在图3(a)～3(d)给出。实验2的仿真结果在图4(a)～4(b)给出。实验3的仿真结果在图5给出。实验4的仿真结果在图6(a)～6(b)给出。提出的抗扰预测控制方法有效地解决了微型燃气轮机热电联供系统存在的控制难点，控制效果优于模型预测控制和PID控制，具有跟踪速度快，超调量小，抗模型失配和干扰能力强优点。
【一种微型燃气轮机热电联供系统抗扰预测控制方法】相关文章：
