基于模型预测控制的综合能源系统鲁棒优化运行方法与流程

本公开属于分布式综合能源系统优化运行技术领域，涉及一种基于模型预测控制的综合能源系统鲁棒优化运行方法。

背景技术：

本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。

综合能源系统是综合区域内石油、天然气、电能、热能等多种能源实现多种能源的协调规划，优化运行和互补互济，在满足系统内多元化用能需求的同时，要有效地提升能源利用效率，促进能源可持续发展的新型一体化的能源系统。由于其高效的能源利用率使得综合能源系统在近几年越来越受重视，且逐步成为研究热点。

据发明人了解，现在的综合能源系统的优化运行中的缺陷：现阶段大部分综合能源系统的优化运行均采用日前优化的方法，由于新能源以及部分负荷的快速变化性，单一的日前优化调度并不能快速应对源、荷的变化，这就需要综合能源系统滚动优化，已有的滚动优化的文献是基于点预测的确定性预测值优化系统出力，未考虑源-荷预测的误差对系统运行的影响。

技术实现要素：

本公开为了解决上述问题，提出了一种基于模型预测控制的综合能源系统鲁棒优化运行方法，本公开对园区综合能源系统进行优化调度，可减小负荷以及可再生能源出力预测误差的影响，经济可靠地满足系统用能需求。

根据一些实施例，本公开采用如下技术方案：

一种基于模型预测控制的综合能源系统鲁棒优化运行方法，包括以下步骤：

设置预测时域和滚动优化时域；

更新源荷历史数据；

考虑源-荷随机性，采用高斯过程回归方法进行源-荷预测，获得下一预测域内的预测值和误差范围；

采用min-max模型预测控制方法，建立鲁棒滚动优化模型，根据预测模型的预测数据和预测误差范围，进行未来滚动优化时域的设备出力的鲁棒滚动优化，得到未来滚动优化时域的运行指令。

作为可选择的实施方式，考虑源-荷随机性，采用高斯过程回归方法进行源-荷预测的过程包括：在预测时域内建立带有噪声的回归模型，通过历史数据得到观测值的先验分布，计算观测值与预测值的联合先验分布，进而计算预测值的后验分布。

作为可选择的实施方式，建立鲁棒滚动优化模型的具体过程包括：

表示综合能源系统在t时刻电能平衡关系、热能平衡关系和冷能平衡关系；

构建综合能源系统各设备在t+1的状态函数；

滚动优化阶段以在源-荷预测误差最大情况下系统运行成本最低为目标，构建目标函数。

作为可选择的实施方式，建立鲁棒滚动优化模型时，包括功率约束条件，即综合能源系统的发电机组、可再生能源发电、燃气锅炉、吸收式制冷机和电制冷机的功率约束。

作为可选择的实施方式，建立鲁棒滚动优化模型时，还包括储能装置的爬坡约束。

作为可选择的实施方式，建立鲁棒滚动优化模型时，还包括蓄电池和储热罐的最大容量约束。

作为可选择的实施方式，所述目标函数为：

pgas、pgrid(t)分别表示天然气价格和分时电价；gpgu(t)、gb(t)、egrid(t)分别表示发电机组，燃气锅炉用气量和电网购电量；u、w、m分别表示优化变量、源-荷预测误差和滚动优化时域。

一种基于模型预测控制的综合能源系统鲁棒优化运行系统，包括：

参数设置模块，被配置为设置预测时域和滚动优化时域；

更新模块，被配置为更新源荷历史数据；

源-荷预测模块，被配置为考虑源-荷随机性，采用高斯过程回归方法进行源-荷预测，获得下一预测域内的预测值和误差范围；

鲁棒滚动优化模块，被配置为采用min-max模型预测控制方法，建立鲁棒滚动优化模型，根据预测模型的预测数据和预测误差范围，进行未来滚动优化时域的设备出力的鲁棒滚动优化，得到未来滚动优化时域的运行指令。

一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行所述的一种基于模型预测控制的综合能源系统鲁棒优化运行方法。

一种终端设备，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行所述的一种基于模型预测控制的综合能源系统鲁棒优化运行方法。

与现有技术相比，本公开的有益效果为：

本公开在源-荷预测环节采用高斯过程回归预测，根据实时数据在线滚动预测，在给出预测值的同时给出预测误差范围，将源-荷随机性考虑在内。根据当前实时数据预测下一预测域内的预测值和误差范围，起到反馈校正的作用。

本公开在滚动优化阶段考虑源-荷预测误差采用min-max模型预测控制优化方法，在源-荷预测误差最大的情况下系统运行成本最低，提高系统的鲁棒性。

本公开的源-荷预测误差范围由高斯过程回归滚动预测给出，不同阶段的误差范围不同，在提高系统鲁棒性的同时降低系统的保守性。

附图说明

构成本公开的一部分的说明书附图用来提供对本公开的进一步理解，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。

图1是综合能源系统结构框图；

图2是鲁棒滚动优化时间窗口示意图；

图3是系统优化流程图。

具体实施方式：

下面结合附图与实施例对本公开作进一步说明。

应该指出，以下详细说明都是例示性的，旨在对本公开提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本公开的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

本公开提出了一种基于模型预测控制的综合能源系统鲁棒优化运行策略。源-荷预测采用高斯过程回归方法，在给出预测值的同时给出预测值的误差范围，将预测随机误差考虑在内，滚动优化阶段采用min-max模型预测控制方法，目的是使系统在源-荷预测误差最大的情况下成本最低。源-荷预测值的范围通过滚动预测获得，使得系统具有较强鲁棒性的同时减小系统的保守型。系统结构图如图1所示。

源-荷预测采用高斯过程回归模型在线滚动预测，在预测时域p内给出预测值的同时给出预测值的误差范围，可以将源-荷随机性考虑在内。

滚动优化环节采用min-max模型预测控制方法以预测误差最大情况下系统运行成本最低为优化目标，保证系统的鲁棒性。根据预测模型的预测数据和预测误差范围优化未来m时长的设备出力并执行第一步的优化结果。滚动优化的时间窗口如图2所示。

基于鲁棒模型预测控制的优化运行过程如图3所示。具体包括：

1、源-荷预测

(1)在预测域p内建立带有噪声的回归模型：

y＝f(x)+ξ,ξ～n(0,σ²)(1)

其中f(x)表示潜函数，ξ为服从n(0,σ²)的噪声。

(2)通过历史数据得到观测值y的先验分布：

y～n(0,k(x,x)+σ²i)(2)

其中k(x,x)为对称的正定协方差矩阵，i为n维的单位矩阵。

(3)计算观测值y与预测值y的联合先验分布：

其中k(xi,xj)是用核函数测量的xi与xj的相关性。x*表示测试集，x表示训练集。

(4)预测值y的后验分布为：

其中预测值：方差为：

2、建立鲁棒滚动优化模型

滚动优化阶段以在源-荷预测误差最大情况下系统运行成本最低为目标，目标函数如式(5)所示：

pgas、pgrid(t)分别表示天然气价格和分时电价；gpgu(t)、gb(t)、egrid(t)分别表示发电机组，燃气锅炉用气量和电网购电量；u、w、m分别表示优化变量、源-荷预测误差和滚动优化时域。

综合能源系统在t时刻电能平衡可表示为式(6)：

epgu(t)+egrid(t)+ese(t)+eren(t)+e²err(t)＝eec(t)+el(t)+e¹err(t)(6)

其中egrid为电网购电量，epgu为内燃发电机组的发电量，eren为风电的发电量，ese为蓄电池的放电量，eec为电制冷机的用电量，el为用户电负荷,为可再生能源预测误差，为电负荷预测误差。

综合能源系统供t时刻热能平衡可表示为式(7)：

hext(t)+hb(t)+hsh(t)＝hab(t)+hl(t)+herr(t)(7)

其中hext(t)、hb(t)、hsh(t)、hab(t)和hl(t)分别表示t时刻发电机组的余热回收的热量、燃气锅炉产生的热量、储热罐的供热量、吸收式制冷机输入的热量和用户的热负荷需求；herr(t)表示t时刻热负荷的预测误差。

综合能源系统供t时刻冷能平衡可表示为式(8)：

cec(t)+cab(t)＝cl(t)+cerr(t)(8)

cab(t)、cec(t)和cl(t)分别表示为：吸收式制冷机出力、电制冷机出力和用户的冷负荷需求；cerr(t)表示t时刻冷负荷的预测误差。

蓄电池(t+1)时刻的状态xe(t+1)可表示为：

xe(t+1)＝(1-τ)×xe(t)-ηba×ese(t)(9)

其中τ为蓄电池的自放电系数，ηba为储放电效率。

储热罐(t+1)时刻的状态xh(t+1)可表示为：

xh(t+1)＝(1-μ)×xh(t)-ηt×hsh(t)(10)

其中μ为储热罐的自放热系数，ηt为储放热效率。

由于储能设备的动态特性，考虑预测时域为p滚动优化时域为m时，蓄电池模型可写为：

其中

其中a＝1-τ，储热罐模型与此类似，此处不再赘述。

考虑到系统实际运行情况后，优化中除需满足能量平衡的等式约束外，还应加入如下不等式约束。

0≤epgu(t)≤epgumax(12)

0≤eren(t)≤erenmax(13)

0≤hb(t)≤hbmax(14)

0≤cab(t)≤cabmax(t)(15)

0≤cec(t)≤cecmax(t)(16)

其中epgu(t)、eren(t)、hb(t)、cab(t)和cec(t)分别表示t时刻发电机组、可再生能源发电、燃气锅炉、吸收式制冷机和电制冷机的功率，epgumax、erenmax、hbmax、cabmax(t)和cecmax(t)分别表示发电机组、可再生能源发电、燃气锅炉、吸收式制冷机和电制冷机的额定功率。表示储能装置的爬坡约束。和分别表示蓄电池和储热罐的最大容量。

通过建立以上基于源-荷高斯过程回归预测的综合能源系统鲁棒模型预测控制优化运行策略，考虑了可再生能源和负荷的随机性；min-max优化与模型预测控制相结合可以提高系统的鲁棒性。通过对min-max问题的转换，该优化问题可以由yalmip编程j进行求解实现综合能源系统的鲁棒运行。

本领域内的技术人员应明白，本公开的实施例可提供为方法、系统、或计算机程序产品。因此，本公开可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本公开可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本公开是参照根据本公开实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅为本公开的优选实施例而已，并不用于限制本公开，对于本领域的技术人员来说，本公开可以有各种更改和变化。凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。

上述虽然结合附图对本公开的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本公开的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。