一种微型燃气轮机冷热电三联供系统的协调控制方法技术方案

本发明专利技术公开了微型燃气轮机冷热电三联供系统的协调控制方法，包括离线阶段和在线阶段；离线阶段包括：步骤1-1，简化被控系统；步骤1-2，在MGT-CCHP系统处于稳定运行状态时，依次在各输入变量中加入高斯白噪声扰动，采集输入变量和输出变量的数据作为模型辨识的数据；步骤1-3，根据扰动实验的数据，利用子空间辨识的方法辨识出该三输入三输出系统的离散状态空间模型；在线阶段包括：步骤2-1，观测当前时刻k的状态向量x(k)的值；步骤2-2，据当前时刻k的状态向量x(k)、k-1时刻的输入向量u(k-1)以及辨识出的状态空间模型，推导出未来P个时刻的输出向量的表达式；步骤2-3，优化性能指标；步骤2-4，利用二次规划的方法，求解优化性能指标最小的优化问题，得到当前时刻的控制增量。

【技术实现步骤摘要】

本专利技术涉及一种工业控制方法，特别是一种基于预测控制算法的微型燃气轮机冷热电三联供系统的协调控制方法。
技术介绍
基于微型燃气轮机的冷热电三联供系统(MicroGasTurbineBasedCombinedCoolingHeatingandPowersystem,以下简称为MGT-CCHP系统)是一种具有代表性的分布式能源的形式，由微型燃气轮机、供热系统、制冷系统组成。动力系统利用燃料燃烧产生的热量来发电，发电后产生的烟气余热被送至供热系统和制冷系统实现余热利用。传统的火力发电的能量转换效率在30-35％左右，而三联供系统对低能量品位的烟气余热再利用，其平均能量转换效率可高达80％，具有广阔的应用前景。MGT-CCHP系统的控制包括微型燃气轮机、双效吸收式制冷机和集成系统三个方面，即分别研究微型燃气轮机侧、供热制冷系统侧和整体的动态特性并建模控制。燃气轮机的控制比较成熟，燃气轮机的各重要设备厂商均为各自的设备提供了相对完善的控制系统，学者们也提出使用PID控制器以及一些先进控制器，如神经网络控制、模糊控制、预测控制等来解决燃气轮机控制问题。然而对于微型燃气轮机控制系统的研究基本上是对燃气轮机的横向延伸，忽略了微型然轮机中余热利用需求、变转速运行状况以及回热器惯性等需求。吸收式制冷机控制系统主要采用PID控制器。对于双效溴化锂吸收式制冷机控制系统的研究存在以下问题：一是吸收式制冷机大多数孤立运行，缺乏与微型燃气轮机组成联供系统时的研究；二是大部分吸收式制冷机的负荷输出单一(仅供冷或供热)，调节量少(仅热源调节)，控制系统普遍为单输入单输出结构；三是吸收式制冷机中有很多热交换设备，热惯性较大。上述研究主要针对各装置独立运行，而针对各装置组成的MGT-CCHP集成系统的控制研究较少。组合后的MGT-CCHP系统是一个有多个输出负荷需求(发电、供热、制冷)的多变量系统，有较强的耦合度，因此应当从协调控制的角度来考虑。针对强耦合的特点，有学者提出了解耦PID控制的方法，即先将该多变量系统解耦成多个单变量系统，并分别对解耦后的单变量系统采用PID控制。解耦PID可以很好地解决多变量的问题。根据分析，MGT-CCHP系统是一个多输入多输出的大惯性系统，有输入量约束。而传统的PID控制方法在大惯性以及输入约束的处理方面较弱。
技术实现思路
专利技术目的：本专利技术所要解决的技术问题是针对传统PID控制在处理微型燃气轮机冷热电三联供系统中被控量约束以及大惯性问题的不足，提供一种微型燃气轮机冷热电三联供系统的协调控制方法。为了解决上述技术问题，本专利技术公开了一种一种微型燃气轮机冷热电三联供系统的协调控制方法，包括离线阶段和在线阶段两个部分；离线阶段包括：步骤1-1，简化被控系统；步骤1-2，在MGT-CCHP系统处于稳定运行状态时，依次在各输入变量中加入高斯白噪声扰动，采集输入变量和输出变量的数据作为模型辨识的数据；步骤1-3，根据扰动实验的数据，利用子空间辨识的方法辨识出该三输入三输出系统的离散状态空间模型；在线阶段包括：步骤2-1，利用卡尔曼滤波器来观测当前时刻k的状态向量x(k)的值；步骤2-2，据当前时刻k的状态向量x(k)、k-1时刻的输入向量u(k-1)以及辨识出的状态空间模型，推导出未来P个时刻的输出向量的表达式；步骤2-3，优化性能指标；步骤2-4，利用二次规划的方法，求解优化性能指标最小的优化问题，得到当前时刻的控制增量。步骤1-1中，简化被控系统包括：保持冷媒水流量、活热水流量以及相应入口温度三者恒定；保证冷剂泵始终全开；制冷机中溶液循环泵处于自动运行模式；微型燃气轮机切换到带变速负载运行；将微型燃气轮机冷热电三联供系统方简化为一个三输入变量三输出变量的系统，三个输出变量是汽轮机转速yn、冷媒水温度yTclo以及生活热水温度yThwo，三个输入变量分别是微型燃气轮机的燃料流量控制阀门umbf、再热阀门开度ure以及双效溴化锂吸收式制冷机中的高压冷剂阀开度uhgr；输入变量服从以下约束：0≤umbf,ure,uhgr≤1；选定微型燃气轮机的燃料流量控制阀门umbf、再热阀门开度ure以及双效溴化锂吸收式制冷机中的高压冷剂阀开度uhgr这三个输入量作为控制系统的控制量。步骤1-3中，三输入三输出系统的离散线性定常状态空间模型为：x(k+1)=Ax(k)+Bu(k)y(k)=Cx(k)]]>公式(1)，式中，u(k)∈Rm，y(k)∈Rl分别是k时刻的输入列向量、输出列向量，输入向量u＝[umbf,ure,uhgr]T，输出向量y＝[yn,yTclo,yThwo]T，x(k)∈Rn是k时刻的状态向量，x(k+1)∈Rn是k+1时刻的状态向量；m、l和n分别是输入向量、输出向量和状态向量的维数，Rm，Rl和Rn分别是m、l和n维实空间，分别表示全体m、l和n维实向量构成的集；A是系统矩阵，表示了系统内部状态变量之间的联系；B是控制矩阵，表示各个输入变量如何控制状态变量；C为输出矩阵即观测矩阵，输出变量反映状态变量。采用的状态空间模型中的直接转移矩阵为零矩阵，该模型的参数A，B和C矩阵均可通过步骤1-3辨识得到，本专利技术中采用的模型辨识方法是子空间辨识方法。步骤2-2中，状态空间模型为公式(1)。步骤2-2中，推导出未来P个时刻的输出向量的表达式为：Y^(k)=Sxx(k)+Su1u(k-1)+SuΔU(k)]]>公式(2)，式中，是在k时刻对未来P个时刻(k+1时刻到k+P时刻)的输出向量值的预估值；ΔU(k)＝[Δu(k),Δu(k+1),…,Δu(k+M-1)]T是未来M个时刻(k时刻到k+M-1时刻)的输入增量值，Sx∈RlP×n,Su1∈RlP×m和Su∈RlP×mM分别是由状态空间模型推导出的预测矩阵。Sx∈RlP×n表示当前时刻(k时刻)的状态向量x(k)如何影响未来P个时刻的输出向量预估值；Su1∈RlP×m表示上一时刻(k-1时刻)的输入向量u(k-1)如何影响未来P个时刻的输出向量预估值；Su∈RlP×mM表示未来M个时刻(k时刻到k+M-1时刻)的输入增量值如何影响未来P个时刻的输出向量预估值。RlP×n，RlP×m和RlP×mM分别是lP×n,lP×m和lP×mM维实向量空间，分别表示全体lP×n,lP×m和lP×mM维实矩阵构成的集。m,l和n分别是输入向量、输出向量和状态向量的维数。P和M分别是预测时域和控制时域。Sx=CACA2...CAPT,]]>Su1=CBCB+本文档来自技高网...

【技术保护点】
一种微型燃气轮机冷热电三联供系统的协调控制方法，其特征在于，包括离线阶段和在线阶段两个部分；离线阶段包括：步骤1‑1，简化被控系统；步骤1‑2，在MGT‑CCHP系统处于稳定运行状态时，依次在各输入变量中加入高斯白噪声扰动，采集输入变量和输出变量的数据作为模型辨识的数据；步骤1‑3，根据扰动实验的数据，利用子空间辨识的方法辨识出该三输入三输出系统的离散状态空间模型；在线阶段包括：步骤2‑1，利用卡尔曼滤波器来观测当前时刻k的状态向量x(k)的值；步骤2‑2，据当前时刻k的状态向量x(k)、k‑1时刻的输入向量u(k‑1)以及辨识出的状态空间模型，推导出未来P个时刻的输出向量的表达式；步骤2‑3，优化性能指标；步骤2‑4，利用二次规划的方法，求解优化性能指标最小的优化问题，得到当前时刻的控制增量。

【技术特征摘要】
1.一种微型燃气轮机冷热电三联供系统的协调控制方法，其特征在于，包括离线阶段
和在线阶段两个部分；
离线阶段包括：
步骤1-1，简化被控系统；
步骤1-2，在MGT-CCHP系统处于稳定运行状态时，依次在各输入变量中加入高斯白噪声
扰动，采集输入变量和输出变量的数据作为模型辨识的数据；
步骤1-3，根据扰动实验的数据，利用子空间辨识的方法辨识出该三输入三输出系统的
离散状态空间模型；
在线阶段包括：
步骤2-1，利用卡尔曼滤波器来观测当前时刻k的状态向量x(k)的值；
步骤2-2，据当前时刻k的状态向量x(k)、k-1时刻的输入向量u(k-1)以及辨识出的状态
空间模型，推导出未来P个时刻的输出向量的表达式；
步骤2-3，优化性能指标；
步骤2-4，利用二次规划的方法，求解优化性能指标最小的优化问题，得到当前时刻的
控制增量。
2.根据权利要求1所述的方法，其特征在于，步骤1-1中，简化被控系统包括：
保持冷媒水流量、活热水流量以及相应入口温度三者恒定；保证冷剂泵始终全开；制冷
机中溶液循环泵处于自动运行模式；微型燃气轮机切换到带变速负载运行；
将微型燃气轮机冷热电三联供系统方简化为一个三输入变量三输出变量的系统，三个
输出变量是汽轮机转速yn、冷媒水温度yTclo以及生活热水温度yThwo，三个输入变量分别是微
型燃气轮机的燃料流量控制阀门umbf、再热阀门开度ure以及双效溴化锂吸收式制冷机中的
高压冷剂阀开度uhgr；输入变量服从以下约束：0≤umbf,ure,uhgr≤1；选定微型燃气轮机的燃
料流量控制阀门umbf、再热阀门开度ure以及双效溴化锂吸收式制冷机中的高压冷剂阀开度
uhgr这三个输入量作为控制系统的控制量。
3.根据权利要求1所述的方法，其特征在于，步骤1-3中，三输入三输出系统的离散线性
定常状态空间模型为：
x(k+1)=Ax(k)+Bu(k)y(k)=Cx(k)]]>公式(1)，
式中，u(k)∈Rm，y(k)∈Rl分别是k时刻的输入列向量、输出列向量，输入向量u＝[umbf,
ure,uhgr]T，输出向量y＝[yn,yTclo,yThwo]T，x(k)∈Rn是k时刻的状态向量，x(k+1)∈Rn是k+1时
刻的状态向量；m、l和n分别是输入向量、输出向量和状态向量的维数，Rm，Rl和Rn分别是m、l
和n维实空间，分别表示全体m、l和n维实向量构成的集；A是系统矩阵，表示了系统内部状态
变量之间的联系；B是控制矩阵，表示各个输入变量如何控制状态变量；C为输出矩阵即观测
矩阵，输出变量反映状态变量。
4.根据权利要求3所述的方法，其特征在于，步骤2-2中，状态空间模型为公式(1)。
5.根据权利要求3所述的方法，其特征在于，步骤2-2中，推导...

【专利技术属性】
技术研发人员：朱明娟，沈炯，吴啸，李益国，张俊礼，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32