CN109669355A - 基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 - Google Patents
基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 Download PDFInfo
- Publication number
- CN109669355A CN109669355A CN201811524593.6A CN201811524593A CN109669355A CN 109669355 A CN109669355 A CN 109669355A CN 201811524593 A CN201811524593 A CN 201811524593A CN 109669355 A CN109669355 A CN 109669355A
- Authority
- CN
- China
- Prior art keywords
- time
- control
- module
- fuel quantity
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001816 cooling Methods 0.000 title claims abstract description 30
- 239000000446 fuel Substances 0.000 claims abstract description 66
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 59
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 239000007789 gas Substances 0.000 claims description 40
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 39
- 239000003546 flue gas Substances 0.000 claims description 39
- 239000011159 matrix material Substances 0.000 claims description 31
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 claims description 6
- 238000009529 body temperature measurement Methods 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 description 36
- 238000004088 simulation Methods 0.000 description 28
- 230000008030 elimination Effects 0.000 description 4
- 238000003379 elimination reaction Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000002803 fossil fuel Substances 0.000 description 2
- 241000254158 Lampyridae Species 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法，基于ARX模型的预测误差计算模块根据燃料量和转速进行预测误差计算，然后接入广义预测控制模块；第二模块根据转速设定值和实际转速的偏差以及第三模块的输入，计算出最佳的燃料量；基于Hammerstein模型的预测误差计算模块根据燃料量、烟气旁路阀阀门开度和冷水温度进行预测误差计算，然后接入基于Hammerstein模型的非线性广义预测解耦控制模块；第四模块根据冷水温度设定值与实际值的偏差、第二模块计算的燃料量控制序列和第五模块输入，计算最佳的烟气旁路阀阀门开度。本发明有效克服了系统存在的复杂特性，提高了系统的灵活性和安全性。
Description
技术领域
本发明涉及热工自动控制领域，尤其是一种基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法。
背景技术
随着经济的不断发展，许多国家在减缓气候变化和减少能源供应方面面临严峻挑战。减少温室气体排放和提高能源效率的一个核心解决方案是使用分布式能源系统。分布式能源系统通常由多种模块化和小规模技术组成，位于最终用户附近，可视为传统集中式电网的重要补充。它们具有许多优点，例如低传输损耗、低排放、以及具有多种能源的灵活性，包括化石燃料和可再生能源。
作为一种关键的分布式能源系统形式，微型燃气轮机结合冷电联供系统(MGT-based combined cooling and power system,MGT-CCP)由于高效节能、低排放和运营成本节约在世界范围内得到越来越多关注。由于微型燃气轮机废弃热量能够被连续利用作为热水器或者冷却装置的热源，MGT-CCP系统的平均能源效率可以达到80％，然而传统的化石燃料电厂是30％—35％。另外相比传统的集中式发电系统，MGT-CCP系统体积小，因此可以灵活地安装在小型住宅区或商业区。
目前对MGT-CCP系统的研究主要集中在静态系统配置、运行优化和性能评价等方面。虽然近年来对MGT-CCP系统的动态特性进行了研究，但MGT-CCP系统的控制方法仍然是传统的PI/PID控制。为了提高PID控制器性能，模糊逻辑、神经网络、萤火虫算法等智能方法被使用来整定PID控制器的参数。然而，由于MGT-CCP系统的动态特性复杂，如非线性、多变量强耦合、未知干扰等，基于单输入、单输出回路分离的PID方法已不再适用于满足性能要求，即使PID的参数得到很好的整定。
发明内容
发明目的：为解决上述问题，本发明提出一种基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法，来改善MGT-CCP系统的控制性能。
技术方案：为实现上述发明目的，本发明采用以下技术方案：
基于广义预测控制的微型燃气轮机冷电联供控制系统，包括：微型燃气轮机冷电联供系统模块、广义预测控制模块、基于ARX模型的预测误差计算模块、基于Hammerstein模型的非线性广义预测解耦控制模块和基于Hammerstein模型的预测误差计算模块；基于ARX模型的预测误差计算模块根据燃料量和微型燃气轮机转速进行预测误差计算，然后接入广义预测控制模块；广义预测控制模块根据微型燃气轮机转速设定值和微型燃气轮机实际转速的偏差以及基于ARX模型的预测误差计算模块的输入，计算出最佳的燃料量；基于Hammerstein模型的预测误差计算模块根据燃料量、烟气旁路阀阀门开度和冷水温度进行预测误差计算，然后接入基于Hammerstein模型的非线性广义预测解耦控制模块；基于Hammerstein模型的非线性广义预测解耦控制模块根据冷水温度设定值与实际值的偏差、广义预测控制模块计算的燃料量控制序列和基于Hammerstein模型的预测误差计算模块输入，计算最佳的烟气旁路阀阀门开度。
可选的，微型燃气轮机冷电联供系统模块(1)包括全回热微型燃气轮机和溴化锂单效吸收式制冷机系统。
可选的，ARX模型结构为：
A(z)yt＝B(z)ut-1+δt/△；
其中，A(z)＝1+a1z-1+a2z-2+...+anz-n，n和nb分别是多项式A(z)和B(z)的阶次，z-1是后移算子，an和是模型的参数，是t时刻的转速，ut-1是t-1时刻的燃料量，δt是随机干扰，△是差分算子；
可选的，Hammerstein模型结构为：
其中，nc和nd分别是多项式和的阶次，z-1是后移算子，和是模型的参数，f1和f2是非线性函数，是t-j时刻的烟气旁路阀开度，是t时刻的冷水温度，是随机干扰。
可选的，广义预测控制模块中燃料量计算公式为：
ut＝ut-1+K(GTQG+R)-1GT(W-L-het)；
其中，ut是t时刻的燃料量，K＝[1 0 … 0]1×M，M是GPC的控制步数，G和L是GPC的参数矩阵，Q是GPC的误差加权矩阵，R是GPC的控制加权矩阵，W是GPC的参考轨迹矩阵，T是矩阵转置符号，P是GPC的预测步数；
基于ARX模型的预测误差计算模块中转速预测误差计算公式为：
其中，et是t时刻的基于ARX模型的转速预测误差，是t时刻的转速测量值；
基于Hammerstein模型的预测误差计算模块中预测误差矩阵计算公式为：
其中，是t时刻的预测误差矩阵，是t时刻的转冷水温度测量值。
可选的，基于Hammerstein模型的非线性广义预测解耦控制模块中烟气旁路阀开度计算公式为：
其中，表示多项式f2的逆，是t时刻的烟气旁路阀开度， 是HNGPDC的控制步数，和是HNGPDC的参数矩阵，是HNGPDC的误差加权矩阵，是HNGPDC的控制加权矩阵，是HNGPDC的参考轨迹矩阵， 是HNGPDC的预测步数；
冷水温度预测计算公式为：
其中，是预测的t+j时刻的冷水温度的预测值，和是HNGPDC的控制器系数，gj,0,gj,1,…,gj,j-1是GPC中的模型系数，是t+j-1时刻的烟气旁路阀开度增量，△ut+j-1是t+j-1时刻的燃料量增量。
本发明还提供了基于广义预测控制的微型燃气轮机冷电联供控制方法，包括以下步骤：
(1)设置GPC控制器参数和HNGPDC控制器参数；
(2)采集MGT-CCP系统的转速数据，利用GPC控制器计算燃料量，对转速进行控制；
(3)利用HNGPDC控制器计算烟气旁路阀开度，对冷水温度进行控制。
进一步的，步骤(1)中GPC控制器的参数为：预测时域P＝20，控制时域M＝10，误差权系数控制权系数
HNGPDC控制器的参数为：预测时域控制时域误差权系数控制权系数
进一步的，步骤(2)包括以下步骤：
(21)建立燃料量和转速之间的自回归各态历经(ARX)模型：
A(z)yt＝B(z)ut-1+δt/△；
其中，A(z)＝1+a1z-1+a2z-2+...+anz-n，n和nb分别是多项式A(z)和B(z)的阶次，z-1是后移算子，a1,a2,…,an和是模型的参数，yt是t时刻的转速，ut-1是t-1时刻的燃料量，δt是随机干扰，△是差分算子；
(22)通过下式计算转速t+j时刻的预测值
其中，参数gj,0,…,gj,j-1根据Gj的系数获得，△ut+j-1是t+j-1时刻的燃料量的增量，Ej和Fj通过下面式子求解：
Ej+1＝Ej+fj,0z-1,E1＝1.
(23)根据转速测量值和预测模型输出值，对转速t+j时刻的预测值进行修正：
其中，et是t时刻的转速预测误差， 是t时刻的转速测量值，yt是t时刻的预测模型输出值；
(24)计算最优燃料量ut，进而对转速进行控制：
ut＝ut-1+K(GTQG+R)-1GT(W-L-het)；
其中，ut是t时刻的燃料量，ut-1是t-1时刻的燃料量；K＝[1 0 … 0]1×M；W是参考轨迹矩阵，W＝[ωt+1 … ωt+P]T，ωt+P是转速参考轨迹；lj(t)＝(Gj-gj,0-gj,1z-1-…-gj,j-1z-(j-1))△ut+j-1+Fjyt。
进一步的，步骤(3)包括以下步骤：
(31)建立燃料量、烟气旁路阀开度和冷水温度之间的Hammerstein模型：
其中，nc和nd分别是多项式和的阶次，z-1是后移算子，和是模型的参数，f1和f2是非线性函数，是t-j时刻的烟气旁路阀开度，是t时刻的冷水温度，是随机干扰；
(32)计算冷水温度预测值：
其中，参数根据的系数获得，△ut+j-1是t+j-1时刻的烟气旁路阀开度的增量， 和通过下面式子求解：
(33)根据转速测量值、冷水温度测量值和转速预测模型输出值、冷水温度预测模型输出值，对冷水温度t+j时刻的预测值进行修正：
其中， 是t时刻的预测误差矩阵， 是t时刻的冷水温度测量值；
(34)计算最优烟气旁路阀门开度进而对冷水温度进行控制：
其中，是t时刻的燃料量，是t-1时刻的燃料量； W是参考轨迹矩阵， 是转速参考轨迹；
有益效果：与现有技术相比，本发明具有以下优点：
(1)在MGT-CCP控制系统中第二模块能提高转速的跟踪速度，第四模块能较好的克服MGT-CCP系统的非线性和强耦合，第三模块和第五模块能补偿未知干扰的影响；
(2)与常规的PI控制方法相比，本发明方法转速和冷水温度跟踪速度快、且无超调，提高了MGT-CCP系统的灵活性；
(3)与常规的PI控制方法相比，本发明方法受未知干扰影响明显降低，未知干扰消除时间更短，增强了MGT-CCP系统的安全性。
附图说明
图1是本发明基于广义预测控制的微型燃气轮机冷电联供控制系统结构示意图；
图2是实验1中转速仿真结果拟合曲线；
图3是实验1中冷水温度仿真结果拟合曲线；
图4是实验1中燃料量仿真结果拟合曲线；
图5是实验1中烟气旁路阀开度仿真结果拟合曲线；
图6是实验2中转速仿真结果拟合曲线；
图7是实验2中冷水温度仿真结果拟合曲线；
图8是实验2中燃料量仿真结果拟合曲线；
图9是实验2中烟气旁路阀开度仿真结果拟合曲线；
图10是实验3中燃料量输入端未知扰动；
图11是实验3中烟气旁路阀输入端未知扰动；
图12是实验3中转速仿真结果拟合曲线；
图13是实验3中冷水温度仿真结果拟合曲线；
图14是实验3中燃料量仿真结果拟合曲线；
图15是实验3中烟气旁路阀开度仿真结果拟合曲线。
具体实施方式
下面结合附图和具体实施例对本发明的技术方案进行详细说明。
如图1所示，基于广义预测控制的微型燃气轮机冷电联供控制系统，包括：微型燃气轮机冷电联供系统模块1(简称第一模块)、广义预测控制模块(Generalized predictivecontrol，GPC)2(简称第二模块)、基于ARX模型的预测误差计算模块3(简称第三模块)、基于Hammerstein模型的非线性广义预测解耦控制模块(Hammerstein model based nonlineargeneralized predictive decoupling control，HNGPDC)4(简称第四模块)和基于Hammerstein模型的预测误差计算模块5(简称第五模块)；第三模块根据燃料量和微型燃气轮机实际转速进行预测误差计算，然后接入第二模块；第二模块根据微型燃气轮机转速设定值和微型燃气轮机实际转速的偏差以及第三模块输出的预测误差，计算出最佳的燃料量；第五模块根据燃料量、烟气旁路阀阀门开度和冷水温度进行预测误差计算，然后接入第四模块；第四模块根据冷水温度设定值与实际值的偏差、第二模块计算的燃料量控制序列和第五模块输出的预测误差，计算最佳的烟气旁路阀阀门开度。
微型燃气轮机冷电联供系统模块1是由全回热微型燃气轮机和溴化锂单效吸收式制冷机系统组成。
基于广义预测控制的微型燃气轮机冷电联供控制方法，包括以下步骤：
(1)设置控制器参数；
GPC控制器(广义预测控制模块)的参数：预测时域P＝20，控制时域M＝10，误差权系数控制权系数
HNGPDC控制器(基于Hammerstein模型的非线性广义预测解耦控制模块)的参数：预测时域控制时域误差权系数控制权系数
(2)采集MGT-CCP系统的转速数据，利用GPC控制器计算燃料量，对转速进行控制，具体包括以下步骤：
(21)建立燃料量和转速之间的ARX模型：
A(z)yt＝B(z)ut-1+δt/△ (1)；
其中，A(z)＝1+a1z-1+a2z-2+...+anz-n，n和nb分别是多项式A(z)和B(z)的阶次，z-1是后移算子，a1,a2,…,an和是模型的参数，yt是t时刻的转速，ut-1是t-1时刻的燃料量，δt是随机干扰，△是差分算子。
(22)通过下式计算转速t+j时刻的预测值
其中，参数gj,0,…,gj,j-1根据Gj的系数获得，△ut+j-1是t+j-1时刻的燃料量的增量，Ej和Fj通过下面式子求解：
Ej+1＝Ej+fj,0z-1,E1＝1.
(23)根据转速测量值和预测模型输出值yt，对转速t+j时刻的预测值进行修正：
其中，et是t时刻的转速预测误差， 是t时刻的转速测量值，yt是t时刻的预测模型输出值。
(24)计算最优燃料量ut，进而对转速进行控制：
ut＝ut-1+K(GTQG+R)-1GT(W-L-het) (4)；
其中，ut是t时刻的燃料量，ut-1是t-1时刻的燃料量；K＝[1 0 … 0]1×M；W是参考轨迹矩阵，W＝[ωt+1 … ωt+P]T，ωt+P是转速参考轨迹；lj(t)＝(Gj-gj,0-gj,1z-1-…-gj,j-1z-(j-1))△ut+j-1+Fjyt。
(3)利用HNGPDC控制器计算烟气旁路阀开度，对冷水温度进行控，具体包括以下步骤：
(31)建立燃料量、烟气旁路阀开度和冷水温度之间的Hammerstein模型：
其中，nc和nd分别是多项式和的阶次，z-1是后移算子，和是模型的参数，f1和f2是非线性函数，是t-j时刻的烟气旁路阀开度，是t时刻的冷水温度，是随机干扰。
(32)计算冷水温度预测值：
其中，参数根据的系数获得，△ut+j-1是t+j-1时刻的烟气旁路阀开度的增量， 和通过下面式子求解：
(33)根据转速测量值、冷水温度测量值和转速预测模型输出值、冷水温度预测模型输出值，对冷水温度t+j时刻的预测值进行修正：
其中， 是t时刻的预测误差矩阵， 是t时刻的冷水温度测量值。
(34)计算最优烟气旁路阀门开度进而对冷水温度进行控制：
其中，是t时刻的燃料量，是t-1时刻的燃料量； W是参考轨迹矩阵， 是转速参考轨迹；
(4)下一个采样周期，重复步骤(2)和(3)。
为了对比说明本发明方法的优越性，设计了常规的PI控制器，参数选择如下：转速PI控制器参数，比例系数KP1＝8.3×10-8，积分系数KI1＝7.6×10-9；冷水温度PI控制器参数，比例系数KP2＝40，积分系数KI2＝0.13。
为了全面验证本发明的有效性，设计了下面3组实验：
实验1：冷水温度和转速跟踪实验，仿真结果如图2-图5所示。图2是实验1中转速仿真结果拟合曲线；图3是实验1中冷水温度仿真结果拟合曲线；图4是实验1中燃料量仿真结果拟合曲线；图5是实验1中烟气旁路阀开度仿真结果拟合曲线；从图2和图3可以看出，本发明方法相比常规的PI控制器，冷负荷和电负荷跟踪速度快、且无超调。从图2和图3可以看出，本发明方法相比常规的PI控制器，控制器输出波动更小，有利于执行器保护。
实验2：解耦效果验证实验，转速设定值先阶跃变化，然后以正弦形式变化，冷水温度设定值固定不变，仿真结果如图6-图9所示。图6是实验2中转速仿真结果拟合曲线；图7是实验2中冷水温度仿真结果拟合曲线；图8是实验2中燃料量仿真结果拟合曲线；图9是实验2中烟气旁路阀开度仿真结果拟合曲线。从仿真结果可以看出，本发明方法控制下冷水温度波动明显减小，表明本发明方法解耦效果相比常规的PI控制有明显改善。
实验3：未知输入扰动实验，先后在燃料量和烟气旁路阀加图10和图11所示的未知扰动，仿真结果如图12-图15所示。图12是实验3中转速仿真结果拟合曲线；图13是实验3中冷水温度仿真结果拟合曲线；图14是实验3中燃料量仿真结果拟合曲线；图15是实验3中烟气旁路阀开度仿真结果拟合曲线。从仿真结果可以看出，本发明方法相比常规的PI控制器，对于未知输入扰动消除时间明显缩短、且偏离设定值幅度明显减小，大幅提高了MGT-CCP系统的安全性。
综合实验1、实验2和实验3的结果，本发明方法相比现有技术，转速和冷水温度跟踪速度快、且无超调，提高了MGT-CCP系统的灵活性；另外本发明方法抗干扰性能比现有控制技术更好，受未知干扰影响更小，未知干扰消除时间更短，增强了MGT-CCP系统的安全性。
总之，本发明的一种基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法，控制系统包括：微型燃气轮机冷电联供系统模块、广义预测控制模块、基于ARX模型的预测误差计算模块、基于Hammerstein模型的非线性广义预测解耦控制模块和基于Hammerstein模型的预测误差计算模块；第三模块根据燃料量和微型燃气轮机转速进行预测误差计算，然后接入第二模块；第二模块根据微型燃气轮机转速设定值和微型燃气轮机实际转速的偏差以及第三模块的输入，计算出最佳的燃料量；第五模块根据燃料量、烟气旁路阀阀门开度和冷水温度进行预测误差计算，然后接入第四模块；第四模块根据冷水温度设定值与实际值的偏差、第二模块计算的燃料量控制序列和第五模块输入，计算最佳的烟气旁路阀阀门开度；本发明系统有效克服MGT-CCP系统存在的非线性、强耦合、存在未知干扰等复杂特性。本发明方法相比现有技术，转速和冷水温度跟踪速度快、且无超调，提高了MGT-CCP系统的灵活性；另外本发明方法抗干扰性能比现有控制技术更好，受未知干扰影响更小，未知干扰消除时间更短，增强了MGT-CCP系统的安全性。
Claims (10)
1.基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于，包括：微型燃气轮机冷电联供系统模块(1)、广义预测控制模块(2)、基于ARX模型的预测误差计算模块(3)、基于Hammerstein模型的非线性广义预测解耦控制模块(4)和基于Hammerstein模型的预测误差计算模块(5)；基于ARX模型的预测误差计算模块(3)根据燃料量和转速进行预测误差计算，然后接入广义预测控制模块(2)；广义预测控制模块(2)根据转速设定值和实际转速的偏差以及基于ARX模型的预测误差计算模块(3)的输入，计算出最佳的燃料量；基于Hammerstein模型的预测误差计算模块(5)根据燃料量、烟气旁路阀阀门开度和冷水温度进行预测误差计算，然后接入基于Hammerstein模型的非线性广义预测解耦控制模块(4)；基于Hammerstein模型的非线性广义预测解耦控制模块(4)根据冷水温度设定值与实际值的偏差、广义预测控制模块(2)计算的燃料量控制序列和基于Hammerstein模型的预测误差计算模块(5)输入，计算最佳的烟气旁路阀阀门开度。
2.根据权利要求1所述的基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于：微型燃气轮机冷电联供系统模块(1)包括全回热微型燃气轮机和溴化锂单效吸收式制冷机系统。
3.根据权利要求1所述的基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于，ARX模型结构为：
A(z)yt＝B(z)ut-1+δt/△；
其中，A(z)＝1+a1z-1+a2z-2+...+anz-n，n和nb分别是多项式A(z)和B(z)的阶次，z-1是后移算子，an和是模型的参数，是t时刻的转速，ut-1是t-1时刻的燃料量，δt是随机干扰，△是差分算子。
4.根据权利要求1所述的基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于，Hammerstein模型结构为：
其中，nc和nd分别是多项式和的阶次，z-1是后移算子，和是模型的参数，f1和f2是非线性函数，是t-j时刻的烟气旁路阀开度，是t时刻的冷水温度，是随机干扰。
5.根据权利要求1所述的基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于，广义预测控制模块(2)中燃料量计算公式为：
ut＝ut-1+K(GTQG+R)-1GT(W-L-het)；
其中，ut是t时刻的燃料量，K＝[1 0 … 0]1×M，M是GPC的控制步数，G和L是GPC的参数矩阵，Q是GPC的误差加权矩阵，R是GPC的控制加权矩阵，W是GPC的参考轨迹矩阵，T是矩阵转置符号，P是GPC的预测步数；
基于ARX模型的预测误差计算模块(3)中转速预测误差计算公式为：
其中，et是t时刻的基于ARX模型的转速预测误差，是t时刻的转速测量值；
基于Hammerstein模型的预测误差计算模块(5)中预测误差矩阵计算公式为：
其中，是t时刻的预测误差矩阵，是t时刻的转冷水温度测量值。
6.根据权利要求1所述的基于广义预测控制的微型燃气轮机冷电联供控制系统，其特征在于，基于Hammerstein模型的非线性广义预测解耦控制模块(4)中烟气旁路阀开度计算公式为：
其中，f2 -1表示多项式f2的逆，是t时刻的烟气旁路阀开度， 是HNGPDC的控制步数，和是HNGPDC的参数矩阵，是HNGPDC的误差加权矩阵，是HNGPDC的控制加权矩阵，是HNGPDC的参考轨迹矩阵， 是HNGPDC的预测步数；
冷水温度预测计算公式为：
其中，是预测的t+j时刻的冷水温度的预测值， 和是HNGPDC的控制器系数，gj,0,gj,1,…,gj,j-1是GPC中的模型系数，是t+j-1时刻的烟气旁路阀开度增量，△ut+j-1是t+j-1时刻的燃料量增量。
7.基于广义预测控制的微型燃气轮机冷电联供控制方法，其特征在于，包括以下步骤：
(1)设置GPC控制器参数和HNGPDC控制器参数；
(2)采集MGT-CCP系统的转速数据，利用GPC控制器计算燃料量，对转速进行控制；
(3)利用HNGPDC控制器计算烟气旁路阀开度，对冷水温度进行控制。
8.根据权利要求7所述的一种基于广义预测控制的微型燃气轮机冷电联供控制方法，其特征在于，步骤(1)中GPC控制器的参数为：预测时域P＝20，控制时域M＝10，误差权系数控制权系数
HNGPDC控制器的参数为：预测时域控制时域误差权系数控制权系数
9.根据权利要求7所述的基于广义预测控制的微型燃气轮机冷电联供控制方法，其特征在于，步骤(2)包括以下步骤：
(21)建立燃料量和转速之间的自回归各态历经(ARX)模型：
A(z)yt＝B(z)ut-1+δt/△；
其中，A(z)＝1+a1z-1+a2z-2+...+anz-n，n和nb分别是多项式A(z)和B(z)的阶次，z-1是后移算子，a1,a2,…,an和是模型的参数，yt是t时刻的转速，ut-1是t-1时刻的燃料量，δt是随机干扰，△是差分算子；
(22)通过下式计算转速t+j时刻的预测值
其中，参数gj,0,…,gj,j-1根据Gj的系数获得，△ut+j-1是t+j-1时刻的燃料量的增量，Ej和Fj通过下面式子求解：
Ej+1＝Ej+fj,0z-1,E1＝1.
(23)根据转速测量值和预测模型输出值yt，对转速t+j时刻的预测值进行修正：
其中，et是t时刻的转速预测误差， 是t时刻的转速测量值，yt是t时刻的预测模型输出值。
(24)计算最优燃料量ut，进而对转速进行控制：
ut＝ut-1+K(GTQG+R)-1GT(W-L-het)；
其中，ut是t时刻的燃料量，ut-1是t-1时刻的燃料量；K＝[1 0 … 0]1×M；W是参考轨迹矩阵，W＝[ωt+1 … ωt+P]T，ωt+P是转速参考轨迹；lj(t)＝(Gj-gj,0-gj,1z-1-…-gj,j-1z-(j-1))△ut+j-1+Fjyt。
10.根据权利要求7所述的基于广义预测控制的微型燃气轮机冷电联供控制方法，其特征在于，步骤(3)包括以下步骤：
(31)建立燃料量、烟气旁路阀开度和冷水温度之间的Hammerstein模型：
其中，nc和nd分别是多项式和的阶次，z-1是后移算子，和是模型的参数，f1和f2是非线性函数，是t-j时刻的烟气旁路阀开度，是t时刻的冷水温度，是随机干扰；
(32)计算冷水温度预测值：
其中，参数根据的系数获得，△ut+j-1是t+j-1时刻的烟气旁路阀开度的增量， 和通过下面式子求解：
(33)根据转速测量值、冷水温度测量值和转速预测模型输出值、冷水温度预测模型输出值，对冷水温度t+j时刻的预测值进行修正：
其中， 是t时刻的预测误差矩阵， 是t时刻的冷水温度测量值；
(34)计算最优烟气旁路阀门开度进而对冷水温度进行控制：
其中，是t时刻的燃料量，是t-1时刻的燃料量； W是参考轨迹矩阵， 是转速参考轨迹；
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811524593.6A CN109669355B (zh) | 2018-12-13 | 2018-12-13 | 基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811524593.6A CN109669355B (zh) | 2018-12-13 | 2018-12-13 | 基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109669355A true CN109669355A (zh) | 2019-04-23 |
| CN109669355B CN109669355B (zh) | 2021-10-19 |
Family
ID=66143772
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811524593.6A Active CN109669355B (zh) | 2018-12-13 | 2018-12-13 | 基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109669355B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488610A (zh) * | 2019-09-04 | 2019-11-22 | 东南大学 | 一种基于鲁棒模糊预测控制的微型燃气轮机热电联供系统热负荷控制方法 |
| CN111596546A (zh) * | 2020-05-19 | 2020-08-28 | 东南大学 | 一种微型燃气轮机热电联供系统抗扰预测控制方法 |
| CN111694274A (zh) * | 2020-05-19 | 2020-09-22 | 东南大学 | 一种基于扰动反馈补偿的热工过程h无穷控制系统 |
| CN111736464A (zh) * | 2020-05-19 | 2020-10-02 | 东南大学 | 一种基于广义扩增状态观测器的热工过程h无穷控制方法 |
| CN111812979A (zh) * | 2020-06-17 | 2020-10-23 | 东南大学 | 一种应用于双效溴化锂吸收式制冷系统的多变量广义预测控制方法 |
| CN115407660A (zh) * | 2022-09-01 | 2022-11-29 | 西安热工研究院有限公司 | 基于动态矩阵控制的火电厂烟气换热器控制方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101667013A (zh) * | 2009-09-04 | 2010-03-10 | 天津大学 | 微型燃气轮机冷电联供分布式供能系统优化运行控制方法 |
| CN105676647A (zh) * | 2016-03-11 | 2016-06-15 | 东南大学 | 一种基于状态观测器的冷热电三联供系统预测控制方法 |
| CN106970533A (zh) * | 2017-05-18 | 2017-07-21 | 吉林大学 | 一种基于rbf‑arx模型稳定参数估计的非线性系统建模方法 |
| CN107807524A (zh) * | 2017-10-20 | 2018-03-16 | 东南大学 | 一种微型燃气轮机冷热电三联供系统鲁棒自适应控制方法 |
| CN108345571A (zh) * | 2018-01-18 | 2018-07-31 | 东南大学 | 微型燃气轮机冷热电联供系统预测控制方法与控制系统 |
-
2018
- 2018-12-13 CN CN201811524593.6A patent/CN109669355B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101667013A (zh) * | 2009-09-04 | 2010-03-10 | 天津大学 | 微型燃气轮机冷电联供分布式供能系统优化运行控制方法 |
| CN105676647A (zh) * | 2016-03-11 | 2016-06-15 | 东南大学 | 一种基于状态观测器的冷热电三联供系统预测控制方法 |
| CN106970533A (zh) * | 2017-05-18 | 2017-07-21 | 吉林大学 | 一种基于rbf‑arx模型稳定参数估计的非线性系统建模方法 |
| CN107807524A (zh) * | 2017-10-20 | 2018-03-16 | 东南大学 | 一种微型燃气轮机冷热电三联供系统鲁棒自适应控制方法 |
| CN108345571A (zh) * | 2018-01-18 | 2018-07-31 | 东南大学 | 微型燃气轮机冷热电联供系统预测控制方法与控制系统 |
Non-Patent Citations (5)
| Title |
|---|
| CHEN CHEN 等: "Improving Simultaneous Cooling and Power Load-Following Capability for MGT-CCP Using Coordinated Predictive Controls", 《ENERGIES》 * |
| 杨得金 等: "火电厂CO2捕集系统的LSSVM-Hammerstein建模方法", 《工程热物理学报》 * |
| 潘蕾 等: "微燃机三联供机组的多变量鲁棒自适应负荷跟踪控制", 《中国电机工程学报》 * |
| 潘蕾 等: "微燃机-冷热电联供机组的Hammerstein模型", 《东南大学学报（自然科学版）》 * |
| 赵德材: "微燃机-LiBr制冷机分布式能源系统建模及非线性控制", 《CNKI》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488610A (zh) * | 2019-09-04 | 2019-11-22 | 东南大学 | 一种基于鲁棒模糊预测控制的微型燃气轮机热电联供系统热负荷控制方法 |
| CN110488610B (zh) * | 2019-09-04 | 2022-05-10 | 东南大学 | 一种基于鲁棒模糊预测控制的微型燃气轮机热电联供系统热负荷控制方法 |
| CN111596546A (zh) * | 2020-05-19 | 2020-08-28 | 东南大学 | 一种微型燃气轮机热电联供系统抗扰预测控制方法 |
| CN111694274A (zh) * | 2020-05-19 | 2020-09-22 | 东南大学 | 一种基于扰动反馈补偿的热工过程h无穷控制系统 |
| CN111736464A (zh) * | 2020-05-19 | 2020-10-02 | 东南大学 | 一种基于广义扩增状态观测器的热工过程h无穷控制方法 |
| CN111736464B (zh) * | 2020-05-19 | 2022-04-26 | 东南大学 | 一种基于广义扩增状态观测器的热工过程h无穷控制方法 |
| CN111596546B (zh) * | 2020-05-19 | 2022-06-03 | 东南大学 | 一种微型燃气轮机热电联供系统抗扰预测控制方法 |
| CN111694274B (zh) * | 2020-05-19 | 2022-06-14 | 东南大学 | 一种基于扰动反馈补偿的热工过程h无穷控制系统 |
| CN111812979A (zh) * | 2020-06-17 | 2020-10-23 | 东南大学 | 一种应用于双效溴化锂吸收式制冷系统的多变量广义预测控制方法 |
| CN115407660A (zh) * | 2022-09-01 | 2022-11-29 | 西安热工研究院有限公司 | 基于动态矩阵控制的火电厂烟气换热器控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109669355B (zh) | 2021-10-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109669355B (zh) | 基于广义预测控制的微型燃气轮机冷电联供控制系统及控制方法 | |
| CN103225799B (zh) | 一种火电厂主蒸汽温度控制方法 | |
| CN107290968B (zh) | 一种多能互补综合能源系统的协调优化控制方法 | |
| CN102621945B (zh) | 基于火力发电机组最优工况的效率动态寻优运行闭环优化控制方法 | |
| CN102401371B (zh) | 一种基于多变量预测控制的再热汽温优化控制方法 | |
| CN102494336B (zh) | 一种循环流化床锅炉燃烧过程多变量控制方法 | |
| CN106786677B (zh) | 一种互联电力系统分布式动态矩阵频率控制方法 | |
| Zeng et al. | Optimization control for the coordinated system of an ultra-supercritical unit based on stair-like predictive control algorithm | |
| CN101609312A (zh) | 采用汽耗率修正滑压运行曲线的发电机组优化控制方法 | |
| Beus et al. | Application of an adaptive model predictive control algorithm on the Pelton turbine governor control | |
| CN103439941A (zh) | 一种燃气内燃机冷热电三联供系统优化运行方法 | |
| Çetin et al. | Development of PID based control strategy in maximum exergy efficiency of a geothermal power plant | |
| Ma et al. | Dynamic control method of flue gas heat transfer system in the waste heat recovery process | |
| CN110932320A (zh) | 自动发电控制系统分布式模型预测控制器设计方法 | |
| Yang et al. | Coordination control for Integrated Solar Combined Cycle thermoelectric coupling | |
| Xu et al. | Design of type-2 fuzzy fractional-order proportional-integral-derivative controller and multi-objective parameter optimization under load reduction condition of the pumped storage unit | |
| CN107807524A (zh) | 一种微型燃气轮机冷热电三联供系统鲁棒自适应控制方法 | |
| CN108131238A (zh) | 一种抑制水锤压力波动的新型pid控制方法 | |
| CN106855691A (zh) | 用于超临界火电机组机炉系统负荷跟踪的双层控制系统 | |
| CN114371619A (zh) | 一种mgt-cchp变工况动态能效优化控制方法 | |
| Lu et al. | A two-dimensional design and synthesis method for coordinated control of flexible-operational combined cycle of gas turbine | |
| CN113555907A (zh) | 一种考虑非理想通信情况的分布式多能系统滚动优化控制方法 | |
| Albanesi et al. | Optimization of the start-up procedure of a combined cycle power plant | |
| CN212565736U (zh) | 一种二次再热主蒸汽温度串级控制系统 | |
| CN103345161B (zh) | 废塑料裂解炉余热烘干装置压力控制方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |