CN114282722A - 分布式综合能源系统多目标优化方法、系统、装置及介质 - Google Patents

分布式综合能源系统多目标优化方法、系统、装置及介质 Download PDF

Info

Publication number
CN114282722A
CN114282722A CN202111595408.4A CN202111595408A CN114282722A CN 114282722 A CN114282722 A CN 114282722A CN 202111595408 A CN202111595408 A CN 202111595408A CN 114282722 A CN114282722 A CN 114282722A
Authority
CN
China
Prior art keywords
model
energy system
day
optimization
distributed
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.)
Pending
Application number
CN202111595408.4A
Other languages
English (en)
Inventor
程序
吕志鹏
王超
周珊
康琦
宋振浩
潘正阳
杨晓霞
张智慧
胡宝玉
史超
米袆仑
薛琳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Online Shanghai Energy Internet Research Institute Co ltd
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing Electric Power Co Ltd
Original Assignee
China Online Shanghai Energy Internet Research Institute Co ltd
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Online Shanghai Energy Internet Research Institute Co ltd, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Beijing Electric Power Co Ltd filed Critical China Online Shanghai Energy Internet Research Institute Co ltd
Priority to CN202111595408.4A priority Critical patent/CN114282722A/zh
Publication of CN114282722A publication Critical patent/CN114282722A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种分布式综合能源系统多目标优化方法、系统、装置及介质,方法包括:建立精细化时间尺度的分布式综合能源系统各设备元件模型,获取分布式综合能源系统的运行约束条件;分别建立分布式综合能源系统的环境效益模型、经济成本、用户舒适度精细化时间尺度模型;构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集;确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。通过在日内调度阶段引入分钟级的滚动优化机制,实现对日前调度方案的实时修正。

Description

分布式综合能源系统多目标优化方法、系统、装置及介质
技术领域
本发明属于综合能源系统优化技术领域,具体涉及一种分布式综合能源系统多目标优化方法、系统、装置及介质。
背景技术
随着电力系统中新能源的增多,其运行难度不断增大,同时为满足用户多样化的能源需求,分布式综合能源系统能够成为电网的重要补充。
现有技术的缺陷和不足:新能源的随机波动因素、用户设备多种运行工况因素对分布式综合能源系统的优化调度的影响将变得不可忽略,而目前分布式综合能源系统中,能源设备、储能设备的模型精细化程度较低,且未能实现通过精细化的模型进行日内多目标优化的实时修正。
目前分布式综合能源系统中的多目标优化调度的模型,大多针对设备的年运行成本最优形成日前的经济性优化调度模型,在对日内优化调度研究中未形成有效的解决方案,因此无法满足综合系统中多设备、复杂工况、多能耦合实时协调控制的需求。且目前分布式综合能源系统中的多目标优化模型对于用户舒适度的考虑较少,精细化程度较低,无法满足用户需求多样化、用能成本精细化控制的要求。
发明内容
本发明的目的在于提供一种分布式综合能源系统多目标优化方法、系统、装置及介质。以解决现有技术中,分布式综合能源系统中的多目标优化调度的模型,精细化程度较低,无法满足用户需求多样化、用能出成本精细化控制要求的问题。
为了实现上述目的,本发明采用如下技术方案:
本发明的第一方面,提供了一种基于日内滚动优化调度的分布式综合能源系统多目标优化方法,包括如下步骤:
建立分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;
建立小时级的日前优化调度模型和分钟级的日内优化调度模型;其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,形成分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型;所述环境效益模型以碳减排量最大为目标;所述经济性模型的优化目标中加入了弃风、弃光的惩罚成本和用户舒适度惩罚成本;所述用户舒适度模型以用户舒适度成本最低为目标;
基于所述分布式综合能源系统的约束条件、分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;
基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
可选的,所述综合能源系统各设备元件模型包括:微燃机及余热回收系统模型、分布式光伏模型、分布式风电模型、电力电子变换装置模型、储能电池模型、热泵模型、电辅热装置模型以及储热装置模型。
可选的,所述环境效益模型如下:
Figure BDA0003430359060000021
其中,EDCR分布式综合能源系统的碳减排量,
Figure BDA0003430359060000022
为时段t光伏设备的碳减排量,EWT为时段t风力发电机设备的碳减排量,
Figure BDA0003430359060000023
为时段t微燃机设备的碳减排量;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数Nt为8760个;当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
可选的,所述经济成本模型如下:
Figure BDA0003430359060000024
其中,CDOC分布式综合能源系统的运行成本,CMT,t为微燃机运行成本,Cgrid,t为电网购电成本,CES,t为储能电池运行成本,CWT,PV,t为时段t弃风、弃光惩罚成本,Ccf,t为用户舒适度惩罚成本;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数Nt为8760个;当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
可选的,所述用户舒适度模型:
Figure BDA0003430359060000031
式中,CDCC为用户舒适度成本,k为能源的种类数,Lk,t为k类能源的初始负荷功率,
Figure BDA0003430359060000032
为用户对k类能源进行需求响应后的实际用能量,yk为用户对k类能源的偏好系数;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数为8760个;当日内采用5分钟级的滚动优化,时段t为5分钟,一天中包含的时段个数为288个。
可选的,所述多目标优化模型包括:
分布式综合能源系统的日前优化目标、日内优化目标,分别以最大化年碳减排量、最小化年运行成本、最小化用户舒适度年成本,最大化日碳减排量、最小化日运行成本、最小化用户舒适度日成本为目标函数;
多目标优化模型的约束条件,包括设备出力约束、系统能量平衡约束、设备最小部分负载率约束、设备最大启动频次约束、设备爬坡约束。
可选的,基于数据包络分析法确定帕累托最优解集中各方案的优先级。
本发明的第二方面,提供了一种用于所述基于日内滚动优化调度的分布式综合能源系统多目标优化方法的系统,包括:
第一模型建立模块,用于建立分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;
第二模型建立模块,用于建立小时级的日前优化调度模型和分钟级的日内优化调度模型;其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,形成分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型;所述环境效益模型以碳减排量最大为目标;所述经济性模型的优化目标中加入了弃风、弃光的惩罚成本和用户舒适度惩罚成本;所述用户舒适度模型以用户舒适度成本最低为目标;
第三模型建立模块,用于基于所述分布式综合能源系统的约束条件、分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;
求解模块,用于基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
优先级确定模块,用于基于数据包络分析法确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
本发明的第三方面,提供了一种计算机装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
本发明的第四方面,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
本发明的有益效果如下:
本发明提供的多目标优化算法考虑污染气体排放的危害,因此增加了以碳减排量最大为目标的环境效益模型;并在经济性模型中加入了弃风、弃光的惩罚成本、用户舒适度惩罚成本;在优化目标中增加了以用户舒适度成本最低为目标的用户舒适度模型,;采用采用供能设备、电储能、热储能等精细化时间尺度的模型,能够有效支撑分布式综合能源系统多目标优化模型中日内滚动优化调度方案的实现。在线滚动优化分层优化调度模型,包括小时级的日前优化调度、分钟级的日内优化调度。通过在日内调度阶段引入分钟级的滚动优化机制,实现对日前调度方案的实时修正。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明实施例分布式综合能源系统多目标优化方法的流程图。
图2为本发明实施例中求解多目标优化模型的流程图。
图3为本发明实施例中确定帕累托最优解集中各方案的优先级的流程图。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明,旨在对本发明提供进一步的详细说明。除非另有指明,本发明所采用的所有技术术语与本申请所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式,而并非意图限制根据本发明的示例性实施方式。
本发明的第一方面,提供了一种基于日内滚动优化调度的分布式综合能源系统多目标优化方法,典型的应用场景分布式综合能源系统,包括微燃机及余热回收系统、分布式光伏、分布式风电、电力电子变换装置、储能电池、空气(地源)源热泵、电辅热装置(电锅炉)以及储热装置。
如图1所示,分布式综合能源系统多目标优化方法的具体步骤包括:
S1、建立精细化时间尺度的分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;其中,约束条件包括分布式综合能源系统的运行约束和能量平衡约束等约束条件。
S2、在线滚动分层优化模型,由小时级的日前优化调度模型、分钟级的日内优化调度模型构成。其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,因此形成分布式综合能源系统的环境效益模型、经济成本、用户舒适度精细化时间尺度模型;所述环境效益模型以碳减排量最大为目标;所述经济性模型的优化目标中加入了弃风、弃光的惩罚成本和用户舒适度惩罚成本;所述用户舒适度模型以用户舒适度成本最低为目标;
S3、基于所述分布式综合能源系统的运行约束和能量平衡约束等约束条件、分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;
S4、基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
S5、确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
具体阐述如下:
(一)方法初始化
包括分布式综合能源系统的系统架构、设备主要参数、用户负荷特征及参数、系统环境参数、决策变量以及目标函数的定义及初始设定。
(二)建立精细化时间尺度的分布式综合能源系统各设备元件模型,具体包括微燃机及余热回收系统模型、分布式光伏模型、分布式风电模型、电力电子变换装置模型、储能电池模型、热泵模型(地源热泵)、电辅热装置(电锅炉)模型以及储热装置模型。
1)精细化时间尺度的微燃机及余热回收装置模型
Figure BDA0003430359060000061
Figure BDA0003430359060000062
Figure BDA0003430359060000063
式中,PMT,t为微燃机时段t发电的出力,QMT,t为微燃机时段t排出烟气可利用的热值,IMT,t为微燃机t时刻输入燃料总能量,COPel,t为微燃机实时电效率,COPth,t为微燃机实时热效率,Δt为微燃机运行时间段,NSPMT为微燃机投入运行的个数,MT表示微燃机。fel,t为电负载率,fth,t为热负载率。Cel,t为时段t的额定电功率、Cth,t为时段t的额定热功率。c1~c4,d1~d4分别为根据微燃机参数设置的调节系数。
2)精细化时间尺度的光伏模型
光伏发电系统输出功率的数学模型可以表述为:
EPV,t=Wt×NSPPV,y×S×COPPV
式中,EPV,t为光伏发电系统时段t输出功率,Wt为对象区域辐照强度,NSPPV,y为该运行年份光伏板搭建个数,S为单个光伏板面积,COPPV为运行效率。
3)精细化时间尺度的风力发电机模型
精细化时间尺度的风力发电系统输出功率的数学模型可以表述为:
Figure BDA0003430359060000064
式中,EWT,t为风力发电系统时段t输出功率,CWT为风力发电机的切出功率,vt为环境风速,vrated为风力发电机的额定风速,vin为风力发电机的切入风速,vout为风力发电机的切出风速。a1~a4为根据风力发电机参数设置的调节系数。
4)精细化时间尺度的储能电池模型
精细化时间尺度的储能电池的模型可以表述为:
EES,t=EES,t-1+[Pst,tδst,t-Pex,tex,t]Δt
式中,EES,t为储能电池时段t的储电量,Pex,t为储能电池放电功率,Pst,t为时段t储能电池充电功率,δex,t为时段t储能电池的放电效率,δst,t为时段t储能电池的充电效率,Δt为储能电池运行的时间段。
5)精细化时间尺度的地源热泵模型
精细化时间尺度的地源热泵的数学模型可由下式描述:
Figure BDA0003430359060000071
式中,PHP,t为地源热泵时段t输出的总能量,COPHP,t为地源热泵运行效率,IHP,t为地源热泵输入总能量,CHP为地源热泵额定运行功率,T为运行时间,fHP,t为运行的负载率。NSPHP为地源热泵投入运行的个数,b1~b4为根据地源热泵参数设置的调节系数,Δt为地源热泵运行的时间段。
6)精细化时间尺度的电锅炉模型
精细化时间尺度的电锅炉的能量转换模型可由下式描述:
QEB,t=ηEB,he,t×EEB,t
式中,QEB,t为时段t输出的热功率;EEB,t为电锅炉时段t的耗电量;ηEB,he,t为电锅炉时段t的热效率。
7)精细化时间尺度的储热装置模型
精细化时间尺度的储热装置的模型可由下式描述:
HHS,t=(1-kLOSS)HHS,t-1+[Qst,tkst,t-Qex,t/kex,t]Δt
式中,HHS,t为储热装置时段t的储热量,Qex,t为时段t储热装置放电功率,Qst,t为时段t储热装置充电功率,kex,t为时段t储热装置的放热效率,kst,t为时段t储热装置的吸热效率,kLOSS为储热装置散热损失率,Δt为储热装置运行的时间段。
当日前采用1小时级的滚动优化,因此时段t为1小时,一年中包含的时段个数为8760个。当日内采用5分钟级的滚动优化,因此时段t为5分钟,一天中包含的时段个数为288个。
(三)建立环境效益、经济成本、用户舒适度模型。建立分布式综合能源系统的环境效益模型、经济成本、用户舒适度精细化时间尺度模型,其中,日前滚动优化以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标;日内滚动优化以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,。
1)精细化时间尺度的环境效益模型
分布式综合能源系统的碳减排量主要来自于使用风电、光伏等分布式能源和微燃机等天然气的使用设备,具体如下:
光伏发电设备的碳减排量:
Figure BDA0003430359060000081
式中,
Figure BDA0003430359060000082
为时段t光伏设备的碳减排量,λPV,t为时段t光伏设备的碳减排系数,PPV,t为时段t光伏设备的有功出力。
风力发电机设备的碳减排量:
Figure BDA0003430359060000083
式中,EWT为时段t风力发电机设备的碳减排量,λWT,t为风力发电机设备的碳减排系数,PWT,t为风力发电机设备的有功出力。
微燃机设备的碳减排量:
Figure BDA0003430359060000084
式中,
Figure BDA0003430359060000085
为时段t微燃机设备的碳减排量,λMT,t为微燃机设备的碳减排系数,PMT,t为微燃机设备的有功出力。
因此,精细化时间尺度的分布式综合能源系统的日碳减排量:
Figure BDA0003430359060000091
其中,EDCR分布式综合能源系统的日碳减排量,当日前采用1小时级的滚动优化,因此时段t为1小时,一年中包含的时段个数为8760个。当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
2)精细化时间尺度的经济成本模型
分布式综合能源系统的日前优化和日内滚动优化分别以年运行成本最低、日运行成本最低为优化目标,年运行成本和日运行成本均包括微燃机机组运行成本、电网购电成本、储能电池运行成本、弃风、弃光惩罚成本、用户舒适度惩罚成本,具体计算如下:
微燃机机组的成本函数:
Figure BDA0003430359060000092
式中,CMT,t为微燃机运行成本,Cf为天然气单价,PMT,t为微燃机机组时段t发出的有功功率,ηMT,t为微燃机时段t的发电效率,LNG为天然气热值,
Figure BDA0003430359060000093
为时段t微燃机的用气量。
电网购电成本:
Cgrid,t=Cgrid_buy,t·Pgrid_buy,t-Cgrid_sell,t·Pgrid_sell,t
式中,Cgrid,t为电网购电成本,Cgrid_buy,t为时段t电能购买单价,Cgrid_sell,t为时段t电能出售单价,Pgrid_buy,t为时段t从电网购买的有功功率,Pgrid_sell,t为时段t从电网出售的有功功率。
储能电池运行成本:
Figure BDA0003430359060000094
式中,CES,t为储能电池运行成本,
Figure BDA0003430359060000095
为储能电池运行成本单价,Pex,t为时段t储能电池放电功率,Pst,t为时段t储能电池充电功率。
弃风、弃光惩罚成本:
CWT,PV,t=ρWT,t(Pyc,WT,t-PWT,t)+ρPV,t(Pyc,PV,t-PPV,t)
式中,CWT,PV,t为时段t弃风、弃光惩罚成本,ρWT,t为时段t弃风惩罚成本系数,Pyc,WT,t为时段t风机预测出力,PWT,t为时段t风机实际出力,ρPV,t为时段t弃光惩罚成本系数,Pyc,PV,t为时段t光伏预测出力,PPV,t为光伏实际出力。
用户舒适度惩罚成本:Ccf,t为用户舒适度惩罚成本。
因此,精细化时间尺度的分布式综合能源系统的运行成本:
Figure BDA0003430359060000101
其中,CDOC分布式综合能源系统的日运行成本,当日前采用1小时级的滚动优化,因此时段t为1小时,一年中包含的时段个数为8760个。当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
3)精细化时间尺度的用户舒适度模型
分布式综合能源系统的日前小时级优化和日内分钟级滚动优化分别以用户舒适度年成本最低、用户舒适度日成本最低为优化目标:
Figure BDA0003430359060000102
式中,CDCC为用户舒适度成本,k为能源的种类数,Lk,t为k类能源的初始负荷功率,
Figure BDA0003430359060000103
为用户对k类能源进行需求响应后的实际用能量,yk为用户对k类能源的偏好系数,为正数,yk值越小,可中断负荷值越高。此处,K可取为5,系统中包含的能源有电、空调冷、空调热、蒸汽、生活热水。当日前采用1小时级的滚动优化,因此时段t为1小时,一年中包含的时段个数为8760个。当日内采用5分钟级的滚动优化,因此时段t为5分钟,一天中包含的时段个数为288个。
(四)基于混合整数规划方法,构建分布式综合能源系统的容量配置和运行策略的多目标优化模型,多目标优化模型包括目标函数和约束条件两部分。
1)目标函数
建立分布式综合能源系统的多目标优化模型的目标函数,分布式综合能源系统的日前优化目标、日内优化目标分别以最大化年碳减排量、最小化年运行成本、最小化用户舒适度年成本,最大化日碳减排量、最小化日运行成本、最小化用户舒适度日成本为目标函数,如下式:
max(EDCR)
min(CDOC)
min(CDCC)
2)约束条件
多目标优化模型的约束条件主要包括设备出力、系统能量平衡约束、设备最小部分负载率约束、设备最大启动频次约束、设备爬坡约束。
(1)设备运行约束
微燃机在运行过程时,每一时刻的出力都不能超过其装机容量,且满足爬坡率的上下限,具体约束如下:
Figure BDA0003430359060000111
式中,
Figure BDA0003430359060000112
为微燃机的输出功率,
Figure BDA0003430359060000113
为微燃机的额定功率,ΔPMT为微燃机输入功率变化量,
Figure BDA0003430359060000114
为微燃机爬坡率下限,
Figure BDA0003430359060000115
为微燃机爬坡率上限。
储能充放电约束,储能电池容量约束:
Figure BDA0003430359060000116
式中,EES为储能电池运行功率,PCharge为储能电池充电功率,PDis-charge为储能电池放电功率,
Figure BDA0003430359060000117
为储能电池容量,δex为储能电池的放电效率,δst为储能电池的充电效率,PCharge·PDis-charge=0表示储能电池不同时充放电,
Figure BDA0003430359060000121
为储能电池功率充放电功率下限,
Figure BDA0003430359060000122
为储能电池功率充放电功率上限。
储热充放热约束,储热容量约束:
Figure BDA0003430359060000123
式中,HHS为储热装置功率,HCharge为储热装置充热功率,HDis-charge为储热装置放热功率,
Figure BDA0003430359060000126
为储热装置容量,kex为储热装置的放热效率,kst为储热装置的吸热效率,HCharge·HDis-charge=0表示电储热装置不同时充放热,
Figure BDA0003430359060000124
为储热装置充放热功率下限,
Figure BDA0003430359060000125
为储热装置充放热功率上限。光伏、风力发电机、地源热泵、电锅炉各时刻的处理均不超过各自的装机容量。
(2)能量平衡约束
所述的能量平衡约束包括电能平衡、热能平衡以及冷能平衡,具体如下:
Pe,load(t)=PWT(t)+PPV(t)+PES(t)+PMT,e(t)+Pgrid(t)
Qhe,load(t)=QMT,he(t)+QHS,he(t)+QHP,he(t)+QEB,he(t)+Qbuy,he(t)
Qco,load(t)=QHP,co(t)
式中,Pe,load(t)、Qhe,load(t)、Qco,load(t)分别为系统中电、热、冷负荷;PWT(t)、PPV(t)、PES(t)、PMT,e(t)、Pgrid(t)分别为系统中风力发电机、光伏的电出力、储能电池与电网交互的电量、微燃机的电出力,和系统与电网交互的电量;QMT,he(t)、QHS,he(t)、QHP,he(t)、QEB,he(t)、Qbuy,he(t)分别为系统中微燃机的热出力、储热设备与热网交互的热功率、热泵、电锅炉的热出力,和系统与热网交互的热量;QHP,co(t)为热泵的冷出力。
(3)设备最大启动频次约束
由于微燃机启动时间较长,为保证系统正常运行,对微燃机启停次数加以约束,限制微燃机每次启停至多一次,具体如下:
Figure BDA0003430359060000131
Figure BDA0003430359060000132
Figure BDA0003430359060000133
Figure BDA0003430359060000134
其中,
Figure BDA0003430359060000135
为表示设备每日启停次数的二进制变量;
Figure BDA0003430359060000136
为控制微燃机开关状态的二进制变量。
(五)基于ε-约束法求解日前小时级、日内分钟级多目标在线滚动优化调度的帕累托最优解集。
ε-约束法可以将原本的多目标优化问题转化为单目标问题,通过多为偏好分析线性规划法和逼近于理想解的排序法,最终得到多目标优化的非劣解集,即帕累托前沿,如图2所示,具体步骤如下:
1)分别以经济性目标、环保性目标和用户舒适度为目标函数对系统进行单目标优化,并将求解结果作为帕累托前沿的两个端点;
2)形成日前小时级多目标优化调度,选择其中一个目标函数进行细分,在两个帕累托前沿端点的范围内根据求解需求进行分段,将该目标函数转换为若干段约束条件,在每段约束条件下,以另外一个目标函数为目标对模型进行优化求解,从而得到一系列最优解,连点成线得到帕累托前沿;
3)形成日内分钟级的滚动优化调度,选择其中一个目标函数进行细分,在两个帕累托前沿端点的范围内根据求解需求进行分段,将该目标函数转换为若干段约束条件,在每段约束条件下,以另外一个目标函数为目标对模型进行优化求解,从而得到一系列最优解,连点成线得到帕累托前沿。
(六)基于数据包络分析法确定帕累托最优解集中各方案的优先级,确定分布式综合能源系统的最优容量配置与运行策略,基于在线滚动分层优化调度,并根据日内滚动优化调度对日前调度方案进行实时修正。
作为应用于本发明的实施例,分布式综合能源系统在进行多目标优化时求解得到的一个非支配解集,解集中每一个解之间并没有优劣之分。本实施例采用多目标评价与决策方法,对比分析帕累托前沿上的所有解,进一步优中选优得到系统的综合最优解,即最优配置方案。
方法具体步骤如下:
1)建立综合能源效率最大化为评估决策单元,
Figure BDA0003430359060000141
其中,ηMT,t、ηgrid、ηpv分别为微燃机发电效率、一次能源转化为市电的效率、光伏发电效率。Eele为系统输出的电功率,EMT,t为微燃机的发电功率,Qhe系统输出的热功率,Eim为电网的购电功率,Qco系统输出的冷功率,Epv为光伏系统的发电功率。
2)最优容量配置及运行策略优化方案决策
求解各个评估目标方案的相对效率值,并根据相对效率值的大小进行排序,相对效率值越高的能源设备容量配置及运行策略优化方案相对于其他方案来说表现更好,根据日内滚动优化调度对日前调度方案进行实时修正。
本发明的第二方面,提供了一种用于所述基于日前小时级、日内分钟级在线滚动优化调度的分布式综合能源系统多目标优化方法的系统,包括:
第一模型建立模块,用于建立精细化时间尺度的分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;
第二模型建立模块,用于建立分布式综合能源系统的环境效益模型、经济成本模型、用户舒适度精细化时间尺度的在线分层优化模型。所述在线滚动分层优化模型,由小时级的日前优化调度模型、分钟级的日内优化调度模型构成。其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,因此形成分布式综合能源系统的环境效益模型、经济成本、用户舒适度精细化时间尺度模型;
第三模型建立模块,用于基于所述分布式综合能源系统的运行约束和能量平衡约束等约束条件、环境效益模型、经济成本、用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的容量配置和运行策略的多目标优化模型;
求解模块,用于基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
优先级确定模块,用于确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
本发明的第三方面,提供了一种计算机装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
本发明的第四方面,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。

Claims (10)

1.一种基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,包括如下步骤:
建立分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;
建立小时级的日前优化调度模型和分钟级的日内优化调度模型;其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,形成分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型;所述环境效益模型以碳减排量最大为目标;所述经济性模型的优化目标中加入了弃风、弃光的惩罚成本和用户舒适度惩罚成本;所述用户舒适度模型以用户舒适度成本最低为目标;
基于所述分布式综合能源系统的约束条件、分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;
基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
2.根据权利要求1所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,所述综合能源系统各设备元件模型包括:微燃机及余热回收系统模型、分布式光伏模型、分布式风电模型、电力电子变换装置模型、储能电池模型、热泵模型、电辅热装置模型以及储热装置模型。
3.根据权利要求1所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,所述环境效益模型如下:
Figure FDA0003430359050000011
其中,EDCR分布式综合能源系统的碳减排量,
Figure FDA0003430359050000012
为时段t光伏设备的碳减排量,EWT为时段t风力发电机设备的碳减排量,
Figure FDA0003430359050000013
为时段t微燃机设备的碳减排量;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数Nt为8760个;当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
4.根据权利要求2所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,所述经济成本模型如下:
Figure FDA0003430359050000021
其中,CDOC分布式综合能源系统的运行成本,CMT,t为微燃机运行成本,Cgrid,t为电网购电成本,CES,t为储能电池运行成本,CWT,PV,t为时段t弃风、弃光惩罚成本,Ccf,t为用户舒适度惩罚成本;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数Nt为8760个;当日内采用5分钟级的滚动优化时,时段t为5分钟,一天中包含的时段个数Nt为288个。
5.根据权利要求2所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,所述用户舒适度模型:
Figure FDA0003430359050000022
式中,CDCC为用户舒适度成本,k为能源的种类数,Lk,t为k类能源的初始负荷功率,
Figure FDA0003430359050000023
为用户对k类能源进行需求响应后的实际用能量,yk为用户对k类能源的偏好系数;当日前采用1小时级的滚动优化时,时段t为1小时,一年中包含的时段个数为8760个;当日内采用5分钟级的滚动优化,时段t为5分钟,一天中包含的时段个数为288个。
6.根据权利要求1所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,所述多目标优化模型包括:
分布式综合能源系统的日前优化目标、日内优化目标,分别以最大化年碳减排量、最小化年运行成本、最小化用户舒适度年成本,最大化日碳减排量、最小化日运行成本、最小化用户舒适度日成本为目标函数;
多目标优化模型的约束条件,包括设备出力约束、系统能量平衡约束、设备最小部分负载率约束、设备最大启动频次约束、设备爬坡约束。
7.根据权利要求1所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法,其特征在于,基于数据包络分析法确定帕累托最优解集中各方案的优先级。
8.一种用于所述基于日内滚动优化调度的分布式综合能源系统多目标优化方法的系统,其特征在于,包括:
第一模型建立模块,用于建立分布式综合能源系统各设备元件模型,获取分布式综合能源系统的约束条件;
第二模型建立模块,用于建立小时级的日前优化调度模型和分钟级的日内优化调度模型;其中,小时级的日前优化调度模型以分布式综合能源系统的年碳减排量最大、年运行成本最小、用户舒适度年成本最低为目标,分钟级的日内优化调度模型以分布式综合能源系统的日碳减排量最大、日运行成本最小、用户舒适度日成本最低为目标,形成分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型;所述环境效益模型以碳减排量最大为目标;所述经济性模型的优化目标中加入了弃风、弃光的惩罚成本和用户舒适度惩罚成本;所述用户舒适度模型以用户舒适度成本最低为目标;
第三模型建立模块,用于基于所述分布式综合能源系统的约束条件、分布式综合能源系统精细化时间尺度的环境效益模型、经济成本模型和用户舒适度模型,采用混合整数规划方法,构建分布式综合能源系统的最优容量配置和运行策略的多目标优化模型;
求解模块,用于基于ε-约束法求解多目标优化模型的在线滚动分层优化调度解,以得到帕累托最优解集,其中每个帕累托解对应一个方案;
优先级确定模块,用于基于数据包络分析法确定帕累托最优解集中各方案的优先级,采用优先级最高的解对应的方案作为分布式综合能源系统的最优容量配置与运行策略。
9.一种计算机装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如权利要求1至7任一项所述的基于日内滚动优化调度的分布式综合能源系统多目标优化方法。
CN202111595408.4A 2021-12-23 2021-12-23 分布式综合能源系统多目标优化方法、系统、装置及介质 Pending CN114282722A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111595408.4A CN114282722A (zh) 2021-12-23 2021-12-23 分布式综合能源系统多目标优化方法、系统、装置及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111595408.4A CN114282722A (zh) 2021-12-23 2021-12-23 分布式综合能源系统多目标优化方法、系统、装置及介质

Publications (1)

Publication Number Publication Date
CN114282722A true CN114282722A (zh) 2022-04-05

Family

ID=80998279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111595408.4A Pending CN114282722A (zh) 2021-12-23 2021-12-23 分布式综合能源系统多目标优化方法、系统、装置及介质

Country Status (1)

Country Link
CN (1) CN114282722A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116526485A (zh) * 2023-04-12 2023-08-01 浙江浙达能源科技有限公司 多能耦合分布式能源的调度优化方法、装置和存储介质
CN117689179A (zh) * 2024-01-30 2024-03-12 山东建筑大学 基于多阶段鲁棒的综合能源系统运行优化方法及系统

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116526485A (zh) * 2023-04-12 2023-08-01 浙江浙达能源科技有限公司 多能耦合分布式能源的调度优化方法、装置和存储介质
CN117689179A (zh) * 2024-01-30 2024-03-12 山东建筑大学 基于多阶段鲁棒的综合能源系统运行优化方法及系统
CN117689179B (zh) * 2024-01-30 2024-05-03 山东建筑大学 基于多阶段鲁棒的综合能源系统运行优化方法及系统

Similar Documents

Publication Publication Date Title
CN107482638B (zh) 冷热电联供型微电网多目标动态优化调度方法
CN113112087A (zh) 考虑电热负荷需求响应的综合能源系统运行成本优化方法
CN109523065B (zh) 一种基于改进量子粒子群算法的微能源网优化调度方法
CN111340274A (zh) 一种基于虚拟电厂参与的综合能源系统优化方法和系统
CN106022503A (zh) 面向具有耦合型电冷热需求的微电网容量规划方法
CN113792953B (zh) 一种虚拟电厂优化调度方法及系统
CN112464477A (zh) 计及需求响应的多能耦合综合能源运行仿真方法
CN111244939B (zh) 一种计及需求侧响应的多能互补系统两级优化设计方法
CN114282722A (zh) 分布式综合能源系统多目标优化方法、系统、装置及介质
CN113437752B (zh) 一种含混合储能的综合能源系统运行控制方法
CN112583017A (zh) 考虑储能运行约束的混合微电网能量分配方法及系统
CN110807588A (zh) 一种多能源耦合综合能源系统的优化调度方法
Dai et al. Optimization of CCHP integrated with multiple load, replenished energy, and hybrid storage in different operation modes
CN111668878A (zh) 一种可再生微能源网的优化配置方法和系统
CN110796279A (zh) 一种多能互补优化运行方法与系统
CN113435095A (zh) 一种用于综合能源系统的优化调度的方法及系统
CN107730032B (zh) 一种泛能网电冷热三联供优化调度系统及方法
CN117350419A (zh) 计及柔性负荷和碳流的园区综合能源系统优化运行方法
CN113902225A (zh) 一种综合能源系统优化方法、系统、装置及存储介质
CN113033900A (zh) 一种园区级综合能源系统容量优化配置方法和系统
CN112131712A (zh) 客户侧多能源系统多目标优化方法和系统
CN115271264A (zh) 一种工业园区能源系统调配方法及计算设备
Jintao et al. Optimized operation of multi-energy system in the industrial park based on integrated demand response strategy
CN108197412B (zh) 一种多能源耦合能量管理系统及优化方法
CN115936336B (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