CN107506851A - 一种多区域虚拟电厂综合能源协调调度优化模型 - Google Patents
一种多区域虚拟电厂综合能源协调调度优化模型 Download PDFInfo
- Publication number
- CN107506851A CN107506851A CN201710616739.9A CN201710616739A CN107506851A CN 107506851 A CN107506851 A CN 107506851A CN 201710616739 A CN201710616739 A CN 201710616739A CN 107506851 A CN107506851 A CN 107506851A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- regions
- power
- 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
Links
- 239000007789 gas Substances 0.000 claims abstract description 85
- 238000004146 energy storage Methods 0.000 claims abstract description 59
- 230000005611 electricity Effects 0.000 claims abstract description 49
- 238000010521 absorption reaction Methods 0.000 claims abstract description 19
- 239000002918 waste heat Substances 0.000 claims abstract description 17
- 238000011084 recovery Methods 0.000 claims abstract description 16
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 238000005057 refrigeration Methods 0.000 claims abstract description 9
- 238000009987 spinning Methods 0.000 claims description 23
- 238000010438 heat treatment Methods 0.000 claims description 17
- 238000001816 cooling Methods 0.000 claims description 15
- 230000007613 environmental effect Effects 0.000 claims description 9
- 238000005338 heat storage Methods 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000009194 climbing Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 239000003345 natural gas Substances 0.000 claims description 3
- 241000894007 species Species 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 238000006116 polymerization reaction Methods 0.000 abstract description 5
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000004220 aggregation Methods 0.000 abstract 1
- 230000002776 aggregation Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 239000000567 combustion gas Substances 0.000 description 2
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000005619 thermoelectricity Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Coloring Foods And Improving Nutritive Qualities (AREA)
Abstract
本发明提供了一种多区域虚拟电厂综合能源协调调度优化模型,用以解决聚合位于不同区域的多种分布式能源和负荷的虚拟电厂冷热电协调调度优化建模问题。该模型聚合单元包括燃气轮机、锅炉、风机、光伏、电储能、热储能、电负荷、热负荷、冷负荷以及相应的冷热电系统配套装置,如电制冷机吸收式制冷机、余热回收装置和热交换器。该模型考虑虚拟电厂内不同区域间的冷热电交互,以及区域内的电、热、冷互补问题。本模型能实现虚拟电厂内不同区域和不同类型聚合单元的冷热电协调优化调度,有效降低虚拟电厂成本。
Description
技术领域
本发明属于电力系统电源调度领域,特别涉及一种多区域虚拟电厂综合能源协调调度优化模型。
背景技术
能源关乎国家根本,近年来,随着综合能源系统、能源互联网等理念的推广,能源行业朝着高效、安全和可持续的能源利用模式不断发展。基于分布式能源的冷热电联供系统通过能量的梯级利用,实现了更高能源利用率、更低能源成本以及更好环保性能等多功能目标,成为区域综合能源系统发展的重要方向和形式。然而,包括冷热电联供系统在内的区域综合能源系统由于地理位置分散,往往都是独立运行,彼此间缺乏协调,难免存在资源配置不合理、系统整体安全性不强等问题。而虚拟电厂技术可以很好地解决上述问题,通过先进的通信、计量、控制等技术聚合不同区域的冷热电联供系统,并通过更高层面的软件构架,虚拟电厂能实现不同区域冷热电联供系统的协调优化控制,从而获得更好经济和环保效益。
与微网相比,虚拟电厂更侧重于吸引并聚合不同区域、不同类型的分布式能源参与电网调度和电力市场交易[7-8]。除传输电能外,虚拟电厂还可提供旋转备用服务,因此,虚拟电厂可参与的电力市场包括能量市场和旋转备用市场。同时参与能量和旋转备用市场可提高决策灵活性,进而获取更大收益。就目前而言,还未有虚拟电厂冷热电协调调度优化的研究,此外,现有研究考虑的虚拟电厂大多位于同一区域,极少涉及多区域虚拟电厂和虚拟电厂区域间能量流动问题。
发明内容
发明目的:提供了一种多区域虚拟电厂综合能源协调调度优化模型,解决同时参与能量市场和旋转备用市场下,聚合位于不同区域的多种分布式能源和负荷的虚拟电厂冷热电协调调度优化建模问题。该虚拟电厂模型的聚合单元包括燃气轮机、锅炉、风机、光伏、电储能、热储能、电负荷、热负荷、冷负荷等。
技术方案:一种多区域虚拟电厂综合能源协调调度优化模型,包括以下步骤:
步骤1:建立虚拟电厂内实现冷热电联供的各聚合单元模型;
步骤2:建立多区域多区域虚拟电厂综合能源协调调度优化模型的目标函数;
步骤3:建立多区域多区域虚拟电厂综合能源协调调度优化模型的约束条件。
进一步,步骤1建立虚拟电厂内实现冷热电联供的各聚合单元模型包括:
1)燃气轮机和锅炉模型。
式中:分别为t时段i区域燃气轮机电功率和余热功率;为t时段i区域锅炉热功率;分别为t时段i区域燃气轮机和锅炉天然气耗量;分别为i区域燃气轮机效率、燃气轮机散热损耗率和锅炉效率。
2)余热回收装置模型。
式中:为t时段i区域余热回收装置供热功率;为i区域余热回收装置效率。
3)电制冷机和吸收式制冷机模型。
式中:分别为t时段i区域电制冷机输入电功率和输出冷功率;分别为t时段i区域吸收式制冷机输入热功率和输出冷功率;分别为i区域电制冷机和吸收式制冷机制冷系数。
4)风机模型。
式中:vi,t分别为t时段i区域风机输出功率和实际风速;分别为i区域风机额定功率、额定风速、切入风速和切出风速。
5)光伏模型。
式中:Ii,t分别为t时段i区域光伏输出功率和太阳辐射强度;分别为i区域光伏板能量转换效率和面积。
6)储能模型。
式中:分别为t时段i区域电储能蓄电量和热储能储热量;分别为t-1时段i区域电储能蓄电量和热储能储热量;分别为i区域电储能和热储能自身能量损耗率;分别为t时段i区域电储能充、放电功率和热储能储、放热功率;分别为i区域电储能充、放电效率和热储能储、放热效率。
7)热、冷管道模型。
式中:分别为t时段i区域向j区域输送能量时i区域输出的热功率和j区域输入的热功率;分别为t时段i区域向j区域输送能量时i区域输出的冷功率和j区域输入的冷功率;σtf、σcf分别为单位距离热能和冷能损耗率;为i区域和j区域间的距离。
进一步,步骤2建立多区域虚拟电厂综合能源协调调度优化模型的目标函数包括:
该模型以虚拟电厂净利润最大化(净成本最小化)为优化目标,决策变量包括各区域在能量市场的交易量和向旋转备用市场提供的备用容量、燃气轮机电功率和锅炉热功率、燃气轮机和可中断负荷备用容量、电制冷机和吸收式制冷机功率、电储能充放电功率、热储储放热功率、中断电负荷量、冷热管道输送功率,其目标函数为:
目标函数包含4个部分,分别为i区域电力市场收益燃气轮机和锅炉成本中断负荷成本和环境成本每部分具体表达式如下。
包括能量市场和旋转备用市场的收益:
式中:分别为t时段能量市场和旋转备用市场电价;分别为t时段i区域在能量市场的交易量(正为售电,负为购电)和向旋转备用市场提供的备用容量。
包括燃气轮机和锅炉的燃料成本、燃气轮机的启停成本:
式中:为t时段天然气价格;分别为i区域燃气轮机的启动和停止成本;布尔变量分别表示t时段i区域燃气轮机是否启动、停止,是则置1,否则置0。
表示为当虚拟电厂中断用户负荷时,需支付一定的补偿费用。考虑到不同中断程度对用户的影响不同,将中断补偿价格与负荷中断等级挂钩,中断等级越高,补偿价格越高:
式中:nm为中断等级数;为第m级中断负荷补偿价格;为t时段i区域第m级中断负荷。
为从电网购电以及燃气轮机和锅炉运行时产生污染气体(包括CO2、SO2、NOx、CO)对应的环境损失和受到的罚款。
式中:ne为污染气体的种类;分别为从电网购电、燃气轮机和锅炉运行时产生第e种污染气体的排放强度;Ve、Ye分别为第e种污染气体的环境价值和罚款数量级。
进一步,步骤3建立多区域虚拟电厂综合能源协调调度优化模型的约束条件包括:
1)各区域电、热、冷功率平衡约束。
式中:分别为t时段i区域电、热、冷负荷;为i区域热交换器效率。
2)热、冷管道约束。
式中:为i区域和j区域间热能传输功率上限;布尔变量表示t时段i区域是否向j区域输送热能,是则置1,否则置0,该布尔变量保证每时段i区域和j区域间功率流动方向唯一。
由于冷管道与热管道约束条件形式一致,不再赘述。
3)燃气轮机约束。
式中:Pi mt,max、Pi mt,min分别为i区域燃气轮机电功率上、下限;为t时段i区域燃气轮机备用容量;布尔变量表示t时段i区域燃气轮机是否运行,是则置1,否则置0;分别为t-1时段i区域燃气轮机电功率;ri u、ri d分别为i区域燃气轮机向上、向下爬坡率;tr为备用服务时间;布尔变量表示t-1时段i区域燃气轮机是否运行,是则置1,否则置0。
4)储能约束。
式中:Pi esc,max、Pi esd,max分别为i区域电储能最大充、放电功率;布尔变量分别表示t时段i区域电储能是否充、放电,是则置1,否则置0;分别为i区域电储能蓄电量上、下限和始、末值;分别为初始时段和末尾时段i区域电储能蓄电量。
热储能与电储能约束的形式一致,不再赘述。
5)中断负荷约束。
式中:为i区域第m级负荷中断系数;为t时段i区域中断负荷;为t时段i区域负荷备用容量。
6)备用容量约束。
7)锅炉、余热回收装置、电制冷机和吸收式制冷机约束。
式中:分别为i区域锅炉、余热回收装置、电制冷机和吸收式制冷机功率上限。
有益效果:本发明相对于现有技术:考虑了虚拟电厂内不同区域间的冷热电交互,以及区域内的电、热、冷互补,解决了聚合位于不同区域的多种分布式能源和负荷的虚拟电厂冷热电协调调度优化建模问题。该模型能实现虚拟电厂内不同区域和不同类型聚合单元的冷热电协调优化调度,有效降低虚拟电厂成本。
附图说明
图1为本发明模型建立的流程图;
图2为多区域虚拟电厂结构示意图;
图3为虚拟电厂能源集线器结构图示意图;
图4为夏季和冬季典型日电力市场电价示意图;
图5为夏季和冬季典型日风机和光伏输出功率示意图。
具体实施方式
下面结合附图介绍本发明的一个具体实施例。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
一种多区域虚拟电厂综合能源协调调度优化模型,包括以下步骤:
步骤1:建立虚拟电厂内实现冷热电联供的各聚合单元模型;
步骤1建立虚拟电厂内实现冷热电联供的各聚合单元模型包括:
1)燃气轮机和锅炉模型。
式中:分别为t时段i区域燃气轮机电功率和余热功率;为t时段i区域锅炉热功率;分别为t时段i区域燃气轮机和锅炉天然气耗量;分别为i区域燃气轮机效率、燃气轮机散热损耗率和锅炉效率。
2)余热回收装置模型。
式中:为t时段i区域余热回收装置供热功率;为i区域余热回收装置效率。
3)电制冷机和吸收式制冷机模型。
式中:分别为t时段i区域电制冷机输入电功率和输出冷功率;分别为t时段i区域吸收式制冷机输入热功率和输出冷功率;分别为i区域电制冷机和吸收式制冷机制冷系数。
4)风机模型。
式中:vi,t分别为t时段i区域风机输出功率和实际风速;分别为i区域风机额定功率、额定风速、切入风速和切出风速。
5)光伏模型。
式中:Ii,t分别为t时段i区域光伏输出功率和太阳辐射强度;分别为i区域光伏板能量转换效率和面积。
6)储能模型。
式中:分别为t时段i区域电储能蓄电量和热储能储热量;分别为t-1时段i区域电储能蓄电量和热储能储热量;分别为i区域电储能和热储能自身能量损耗率;分别为t时段i区域电储能充、放电功率和热储能储、放热功率;分别为i区域电储能充、放电效率和热储能储、放热效率。
7)热、冷管道模型。
式中:分别为t时段i区域向j区域输送能量时i区域输出的热功率和j区域输入的热功率;分别为t时段i区域向j区域输送能量时i区域输出的冷功率和j区域输入的冷功率;σtf、σcf分别为单位距离热能和冷能损耗率;为i区域和j区域间的距离。
步骤2:建立多区域虚拟电厂综合能源协调调度优化模型的目标函数;
步骤2建立多区域虚拟电厂综合能源协调调度优化模型的目标函数包括:
该模型以虚拟电厂净利润最大化(净成本最小化)为优化目标,决策变量包括各区域在能量市场的交易量和向旋转备用市场提供的备用容量、燃气轮机电功率和锅炉热功率、燃气轮机和可中断负荷备用容量、电制冷机和吸收式制冷机功率、电储能充放电功率、热储储放热功率、中断电负荷量、冷热管道输送功率,其目标函数为:
目标函数包含4个部分,分别为i区域电力市场收益燃气轮机和锅炉成本中断负荷成本和环境成本每部分具体表达式如下。
包括能量市场和旋转备用市场的收益:
式中:分别为t时段能量市场和旋转备用市场电价;分别为t时段i区域在能量市场的交易量(正为售电,负为购电)和向旋转备用市场提供的备用容量。
包括燃气轮机和锅炉的燃料成本、燃气轮机的启停成本:
式中:为t时段天然气价格;分别为i区域燃气轮机的启动和停止成本;布尔变量分别表示t时段i区域燃气轮机是否启动、停止,是则置1,否则置0。
表示为当虚拟电厂中断用户负荷时,需支付一定的补偿费用。考虑到不同中断程度对用户的影响不同,将中断补偿价格与负荷中断等级挂钩,中断等级越高,补偿价格越高:
式中:nm为中断等级数;为第m级中断负荷补偿价格;为t时段i区域第m级中断负荷。
为从电网购电以及燃气轮机和锅炉运行时产生污染气体(包括CO2、SO2、NOx、CO)对应的环境损失和受到的罚款。
式中:ne为污染气体的种类;分别为从电网购电、燃气轮机和锅炉运行时产生第e种污染气体的排放强度;Ve、Ye分别为第e种污染气体的环境价值和罚款数量级。
步骤3:建立多区域虚拟电厂综合能源协调调度优化模型的约束条件。
步骤3建立多区域虚拟电厂综合能源协调调度优化模型的约束条件包括:
1)各区域电、热、冷功率平衡约束。
式中:分别为t时段i区域电、热、冷负荷;为i区域热交换器效率。
2)热、冷管道约束。
式中:为i区域和j区域间热能传输功率上限;布尔变量表示t时段i区域是否向j区域输送热能,是则置1,否则置0,该布尔变量保证每时段i区域和j区域间功率流动方向唯一。
由于冷管道与热管道约束条件形式一致,不再赘述。
3)燃气轮机约束。
式中:Pi mt,max、Pi mt,min分别为i区域燃气轮机电功率上、下限;为t时段i区域燃气轮机备用容量;布尔变量表示t时段i区域燃气轮机是否运行,是则置1,否则置0;分别为t-1时段i区域燃气轮机电功率;ri u、ri d分别为i区域燃气轮机向上、向下爬坡率;tr为备用服务时间;布尔变量表示t-1时段i区域燃气轮机是否运行,是则置1,否则置0。
4)储能约束。
式中:Pi esc,max、Pi esd,max分别为i区域电储能最大充、放电功率;布尔变量分别表示t时段i区域电储能是否充、放电,是则置1,否则置0;分别为i区域电储能蓄电量上、下限和始、末值;分别为初始时段和末尾时段i区域电储能蓄电量。
热储能与电储能约束的形式一致,不再赘述。
5)中断负荷约束。
式中:为i区域第m级负荷中断系数;为t时段i区域中断负荷;为t时段i区域负荷备用容量。
6)备用容量约束。
7)锅炉、余热回收装置、电制冷机和吸收式制冷机约束。
式中:分别为i区域锅炉、余热回收装置、电制冷机和吸收式制冷机功率上限。
下面以一个虚拟电厂为例介绍本发明:
为验证上述模型,以虚拟电厂模式控制长沙某一多区域冷热电联供型综合能源系统。多区域虚拟电厂结构如图1所示。该虚拟电厂系统分为工业区、商业区和居民区,每个区域均装设有一能源集线器,能源集线器结构如图2所示,其负责供应各区域的电、热、冷负荷。能源集线器内各单元参数如表1所示。
表1能源集线器内各单元参数
夏季和冬季典型日能量市场和旋转备用市场电价预测值见图3,根据当地风速和太阳辐射强度,预测风机和光伏输出功率如图4所示。旋转备用市场备用服务时间为10min。
为衡量虚拟电厂实现冷热电联供、虚拟电厂区域间通过热、冷管道相连以及虚拟电厂同时参与能量市场和旋转备用市场对虚拟电厂成本的影响,设置4种方案如表1所示。
表1 4种不同的虚拟电厂构建方案
采用上述4种方案构建夏季和冬季典型日虚拟电厂调度优化模型,所得电力市场收益($)、各部分成本($)和虚拟电厂净成本($)如表2所示。
表2 4种方案结果对比
从表2夏季和冬季典型日方案1、2、4对比可以看出,相比于方案1,方案2中虚拟电厂通过冷热电联供,有效利用燃气轮机余热,提高了燃气轮机的能源利用率,因此方案2中燃气轮机输出功率明显提高,表现为燃气轮机成本增大。同时,燃气轮机供热的增加和更多热能通过吸收式制冷机供冷也大幅度减少了锅炉输出和能量市场购电量,表现为锅炉成本降低,能量市场收益提高,最终使得虚拟电厂净成本下降。进一步,方案4中虚拟电厂区域间通过热、冷管道相连,轻负荷区域可通过管道向重负荷区域供能,从而提高了轻负荷区域燃气轮机的输出功率,因此虚拟电厂净成本更低。方案3,虚拟电厂只参与能量市场,但未参与旋转备用市场,虽然该方案在高电价时段增加燃气轮机输出功率和中断负荷量获得更大的能量市场收益,但所增净收益少于虚拟电厂在旋转备用市场的获益,其净成本高于方案4。综上所述,虚拟电厂实现冷热电联供、虚拟电厂不同区域间通过热、冷管道相连以及虚拟电厂同时参与能量市场和旋转备用市场均可有效降低虚拟电厂的净成本。
以上仿真结果验证了本发明所构模型有效性和实用性,说明本文中虚拟电厂通过同时参与能量市场和旋转备用市场、不同区域间的热冷功率流动以及各区域内不同聚合单元的电热冷功率互补,极大地提高了调度的灵活性,从而大幅度降低虚拟电厂的净成本。
Claims (4)
1.一种多区域虚拟电厂综合能源协调调度优化模型,其特征在于:包括以下步骤:
步骤1:建立虚拟电厂内实现冷热电联供的各聚合单元模型;
步骤2:建立多区域虚拟电厂综合能源协调调度优化模型的目标函数;
步骤3:建立多区域虚拟电厂综合能源协调调度优化模型的约束条件。
2.根据权利要求1所述的多区域虚拟电厂综合能源协调调度优化模型,其特征在于:所述步骤1建立虚拟电厂内实现冷热电联供的各聚合单元模型包括以下步骤:
1)燃气轮机和锅炉模型:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
</mrow>
式中:分别为t时段i区域燃气轮机电功率和余热功率;为t时段i区域锅炉热功率;分别为t时段i区域燃气轮机和锅炉天然气耗量;分别为i区域燃气轮机效率、燃气轮机散热损耗率和锅炉效率;
2)余热回收装置模型:
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>h</mi>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>h</mi>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msubsup>
</mrow>
式中:为t时段i区域余热回收装置供热功率;为i区域余热回收装置效率;
3)电制冷机和吸收式制冷机模型:
<mrow>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msubsup>
<msubsup>
<mi>K</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msubsup>
<msubsup>
<mi>K</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>p</mi>
<mi>a</mi>
</mrow>
</msubsup>
</mrow>
式中:分别为t时段i区域电制冷机输入电功率和输出冷功率;分别为t时段i区域吸收式制冷机输入热功率和输出冷功率;分别为i区域电制冷机和吸收式制冷机制冷系数;
4)风机模型:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo><</mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msubsup>
<mo>,</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>></mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msubsup>
<mo><</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo><</mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mi>r</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>w</mi>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mi>r</mi>
</msubsup>
<mo><</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo><</mo>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中:vi,t分别为t时段i区域风机输出功率和实际风速;Pi wtr、分别为i区域风机额定功率、额定风速、切入风速和切出风速;
5)光伏模型:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msubsup>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msubsup>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
式中:Ii,t分别为t时段i区域光伏输出功率和太阳辐射强度;分别为i区域光伏板能量转换效率和面积;
6)储能模型:
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
</mfrac>
</mrow>
式中:分别为t时段i区域电储能蓄电量和热储能储热量;分别为t-1时段i区域电储能蓄电量和热储能储热量;分别为i区域电储能和热储能自身能量损耗率;分别为t时段i区域电储能充、放电功率和热储能储、放热功率;分别为i区域电储能充、放电效率和热储能储、放热效率;
7)热、冷管道模型:
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&sigma;</mi>
<mrow>
<mi>t</mi>
<mi>f</mi>
</mrow>
</msup>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>W</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>f</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>f</mi>
</mrow>
</msup>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:分别为t时段i区域向j区域输送能量时i区域输出的热功率和j区域输入的热功率;分别为t时段i区域向j区域输送能量时i区域输出的冷功率和j区域输入的冷功率;σtf、σcf分别为单位距离热能和冷能损耗率;为i区域和j区域间的距离。
3.根据权利要求1所述的多区域虚拟电厂综合能源协调调度优化模型,其特征在于:所述步骤2建立多区域虚拟电厂综合能源协调调度优化模型的目标函数包括以下步骤:
该模型以虚拟电厂净利润最大化为优化目标,决策变量包括各区域在能量市场的交易量和向旋转备用市场提供的备用容量、燃气轮机电功率和锅炉热功率、燃气轮机和可中断负荷备用容量、电制冷机和吸收式制冷机功率、电储能充放电功率、热储储放热功率、中断电负荷量、冷热管道输送功率,其目标函数为:
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>i</mi>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>m</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>n</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
目标函数包含4个部分,分别为i区域电力市场收益燃气轮机和锅炉成本中断负荷成本和环境成本每部分具体表达式如下;
包括能量市场和旋转备用市场的收益:
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>m</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>t</mi>
<mrow>
<mi>e</mi>
<mi>m</mi>
</mrow>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>t</mi>
<mrow>
<mi>s</mi>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
</mrow>
式中:分别为t时段能量市场和旋转备用市场电价;为t时段i区域在能量市场的交易量,其中正为售电,负为购电,为t时段i区域向旋转备用市场提供的备用容量;
包括燃气轮机和锅炉的燃料成本、燃气轮机的启停成本:
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>n</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>t</mi>
<mrow>
<mi>n</mi>
<mi>g</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>u</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
</mrow>
式中:为t时段天然气价格;分别为i区域燃气轮机的启动和停止成本;布尔变量分别表示t时段i区域燃气轮机是否启动、停止,是则置1,否则置0;
表示为当虚拟电厂中断用户负荷时,需支付一定的补偿费用;考虑到不同中断程度对用户的影响不同,将中断补偿价格与负荷中断等级挂钩,中断等级越高,补偿价格越高:
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>m</mi>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&lambda;</mi>
<mi>m</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
式中:nm为中断等级数;为第m级中断负荷补偿价格;为t时段i区域第m级中断负荷;
为从电网购电以及燃气轮机和锅炉运行时产生污染气体对应的环境损失和受到的罚款,污染气体包括CO2、SO2、NOx、CO;
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>e</mi>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>e</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>e</mi>
</msub>
</munderover>
<msubsup>
<mi>D</mi>
<mi>e</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>e</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>e</mi>
</msub>
</munderover>
<msubsup>
<mi>D</mi>
<mi>e</mi>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>e</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>D</mi>
<mi>e</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中:ne为污染气体的种类;分别为从电网购电、燃气轮机和锅炉运行时产生第e种污染气体的排放强度;Ve、Ye分别为第e种污染气体的环境价值和罚款数量级。
4.根据权利要求1所述的多区域虚拟电厂综合能源协调调度优化模型,其特征在于:所述步骤3建立多区域虚拟电厂综合能源协调调度优化模型的约束条件包括以下步骤:
1)各区域电、热、冷功率平衡约束:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>h</mi>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>/</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>h</mi>
<mi>e</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
<mo>:</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
<mo>:</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>+</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
<mo>:</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>f</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
<mo>:</mo>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</munder>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msubsup>
</mrow>
式中:分别为t时段i区域电、热、冷负荷;为i区域热交换器效率;
2)热、冷管道约束:
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>o</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
式中:为i区域和j区域间热能传输功率上限;布尔变量表示t时段i区域是否向j区域输送热能,是则置1,否则置0,该布尔变量保证每时段i区域和j区域间功率流动方向唯一;
冷管道与热管道约束条件形式一致;
3)燃气轮机约束:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>t</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>o</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>t</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>o</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mi>u</mi>
</msubsup>
<msup>
<mi>t</mi>
<mi>r</mi>
</msup>
</mrow>
<mrow>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mi>d</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mi>u</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>o</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>o</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>u</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>o</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>o</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
</mrow>
式中:Pi mt,max、Pi mt,min分别为i区域燃气轮机电功率上、下限;为t时段i区域燃气轮机备用容量;布尔变量表示t时段i区域燃气轮机是否运行,是则置1,否则置0;分别为t-1时段i区域燃气轮机电功率;ri u、ri d分别为i区域燃气轮机向上、向下爬坡率;tr为备用服务时间;布尔变量表示t-1时段i区域燃气轮机是否运行,是则置1,否则置0;
4)储能约束:
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>d</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mo>,</mo>
<mi>min</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mn>24</mn>
</mrow>
<mrow>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
式中:Pi esc,max、Pi esd,max分别为i区域电储能最大充、放电功率;布尔变量分别表示t时段i区域电储能是否充、放电,是则置1,否则置0;分别为i区域电储能蓄电量上、下限和始、末值;分别为初始时段和末尾时段i区域电储能蓄电量;
热储能与电储能约束的形式一致;
5)中断负荷约束:
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>m</mi>
</msub>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>n</mi>
<mi>m</mi>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>k</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msubsup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
式中:为i区域第m级负荷中断系数;为t时段i区域中断负荷;为t时段i区域负荷备用容量;
6)备用容量约束:
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>s</mi>
<mi>r</mi>
<mi>m</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>t</mi>
</mrow>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
7)锅炉、余热回收装置、电制冷机和吸收式制冷机约束:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>h</mi>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mi>i</mi>
<mrow>
<mi>h</mi>
<mi>r</mi>
<mi>s</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>e</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mi>c</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>Q</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mo>,</mo>
<mi>max</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:分别为i区域锅炉、余热回收装置、电制冷机和吸收式制冷机功率上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710616739.9A CN107506851A (zh) | 2017-07-26 | 2017-07-26 | 一种多区域虚拟电厂综合能源协调调度优化模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710616739.9A CN107506851A (zh) | 2017-07-26 | 2017-07-26 | 一种多区域虚拟电厂综合能源协调调度优化模型 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107506851A true CN107506851A (zh) | 2017-12-22 |
Family
ID=60690027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710616739.9A Pending CN107506851A (zh) | 2017-07-26 | 2017-07-26 | 一种多区域虚拟电厂综合能源协调调度优化模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107506851A (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063925A (zh) * | 2018-08-16 | 2018-12-21 | 合肥工业大学 | 一种计及负荷聚合商的区域综合能源系统优化运行方法 |
CN109146320A (zh) * | 2018-09-12 | 2019-01-04 | 河海大学 | 一种考虑配电网安全性的虚拟电厂优化调度方法 |
CN109255487A (zh) * | 2018-09-21 | 2019-01-22 | 国网浙江省电力有限公司 | 一种基于标准化矩阵模型的综合能源系统优化方法 |
CN109617142A (zh) * | 2018-12-13 | 2019-04-12 | 燕山大学 | 一种cchp型微电网多时间尺度优化调度方法及系统 |
CN109713723A (zh) * | 2019-02-13 | 2019-05-03 | 广东工业大学 | 一种区域综合能源系统运行优化方法 |
CN109799708A (zh) * | 2019-01-30 | 2019-05-24 | 广东电网有限责任公司电力调度控制中心 | 一种基于最大内嵌正方体的虚拟电厂灵活性聚合方法 |
CN109859071A (zh) * | 2019-01-28 | 2019-06-07 | 上海电力学院 | 虚拟电厂多能互补的源荷储优化配置方法 |
CN110135613A (zh) * | 2018-10-23 | 2019-08-16 | 上海交通大学 | 一种基于纳什谈判的多虚拟电厂协同优化运行方案 |
CN110137957A (zh) * | 2019-05-22 | 2019-08-16 | 贵州大学 | 一种基于直流配电中心的柔性互联配电网削峰填谷方法 |
CN110263387A (zh) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法 |
CN110288152A (zh) * | 2019-06-25 | 2019-09-27 | 河海大学 | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 |
CN110286590A (zh) * | 2019-06-21 | 2019-09-27 | 河海大学 | 考虑多类型能源存储的冷热电联供系统日前优化调度方法 |
CN110619515A (zh) * | 2019-01-08 | 2019-12-27 | 魏丽 | 电子支付管理方法及存储介质 |
CN110808615A (zh) * | 2019-12-07 | 2020-02-18 | 国家电网有限公司 | 考虑不确定性的气电虚拟电厂调度优化方法 |
CN110929938A (zh) * | 2019-11-25 | 2020-03-27 | 新奥数能科技有限公司 | 能源系统的优化方法及装置、存储介质、电子装置 |
CN111179110A (zh) * | 2019-12-06 | 2020-05-19 | 清华-伯克利深圳学院筹备办公室 | 虚拟电厂可变阶聚合等效鲁棒动态模型建模方法及装置 |
CN111210119A (zh) * | 2019-12-24 | 2020-05-29 | 国网山东省电力公司枣庄供电公司 | 多种市场下vpp电、热、气优化调度模型的建立方法 |
CN111382939A (zh) * | 2020-03-06 | 2020-07-07 | 国网冀北电力有限公司 | 虚拟电厂资源优化配置方法、装置及设备 |
CN111625961A (zh) * | 2020-05-26 | 2020-09-04 | 中国科学院工程热物理研究所 | 综合能源系统协同优化运行调控方法 |
CN111953018A (zh) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | 一种分布式多能互补供能系统及功率分配方法 |
CN112084641A (zh) * | 2020-08-28 | 2020-12-15 | 国网综合能源服务集团有限公司 | 一种基于itsdp模型的园区虚拟电厂优化运行方法 |
CN116937551A (zh) * | 2023-07-19 | 2023-10-24 | 太原理工大学 | 一种电气互联虚拟电厂优化调度方法及终端 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127389A (zh) * | 2016-06-24 | 2016-11-16 | 河海大学 | 一种虚拟电厂热电联合调度鲁棒优化模型 |
-
2017
- 2017-07-26 CN CN201710616739.9A patent/CN107506851A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127389A (zh) * | 2016-06-24 | 2016-11-16 | 河海大学 | 一种虚拟电厂热电联合调度鲁棒优化模型 |
Non-Patent Citations (2)
Title |
---|
周亦洲等: "多区域虚拟电厂综合能源协调调度优化模型", 《HTTPS://KNS.CNKI.NET/KCMS/DETAIL/11.2107.TM.20170704.1537.016.HTML》 * |
王妍: "基于风光储建模的微电网及其经济优化运行研究", 《中国优秀硕士论文全文数据库》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063925A (zh) * | 2018-08-16 | 2018-12-21 | 合肥工业大学 | 一种计及负荷聚合商的区域综合能源系统优化运行方法 |
CN109063925B (zh) * | 2018-08-16 | 2021-08-17 | 合肥工业大学 | 一种计及负荷聚合商的区域综合能源系统优化运行方法 |
CN109146320A (zh) * | 2018-09-12 | 2019-01-04 | 河海大学 | 一种考虑配电网安全性的虚拟电厂优化调度方法 |
CN109146320B (zh) * | 2018-09-12 | 2021-06-04 | 河海大学 | 一种考虑配电网安全性的虚拟电厂优化调度方法 |
CN109255487A (zh) * | 2018-09-21 | 2019-01-22 | 国网浙江省电力有限公司 | 一种基于标准化矩阵模型的综合能源系统优化方法 |
CN110135613B (zh) * | 2018-10-23 | 2023-01-03 | 上海交通大学 | 一种基于纳什谈判的多虚拟电厂协同优化运行方案 |
CN110135613A (zh) * | 2018-10-23 | 2019-08-16 | 上海交通大学 | 一种基于纳什谈判的多虚拟电厂协同优化运行方案 |
CN109617142A (zh) * | 2018-12-13 | 2019-04-12 | 燕山大学 | 一种cchp型微电网多时间尺度优化调度方法及系统 |
CN110619515A (zh) * | 2019-01-08 | 2019-12-27 | 魏丽 | 电子支付管理方法及存储介质 |
CN109859071A (zh) * | 2019-01-28 | 2019-06-07 | 上海电力学院 | 虚拟电厂多能互补的源荷储优化配置方法 |
CN109859071B (zh) * | 2019-01-28 | 2023-04-25 | 上海电力学院 | 虚拟电厂多能互补的源荷储优化配置方法 |
CN109799708B (zh) * | 2019-01-30 | 2021-08-24 | 广东电网有限责任公司电力调度控制中心 | 一种基于最大内嵌正方体的虚拟电厂灵活性聚合方法 |
CN109799708A (zh) * | 2019-01-30 | 2019-05-24 | 广东电网有限责任公司电力调度控制中心 | 一种基于最大内嵌正方体的虚拟电厂灵活性聚合方法 |
CN109713723A (zh) * | 2019-02-13 | 2019-05-03 | 广东工业大学 | 一种区域综合能源系统运行优化方法 |
CN110137957A (zh) * | 2019-05-22 | 2019-08-16 | 贵州大学 | 一种基于直流配电中心的柔性互联配电网削峰填谷方法 |
CN110263387B (zh) * | 2019-05-29 | 2023-04-07 | 国网福建省电力有限公司 | 基于配电网、天然气网的能源系统运行装置 |
CN110263387A (zh) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法 |
CN110286590A (zh) * | 2019-06-21 | 2019-09-27 | 河海大学 | 考虑多类型能源存储的冷热电联供系统日前优化调度方法 |
CN110288152B (zh) * | 2019-06-25 | 2022-08-26 | 河海大学 | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 |
CN110288152A (zh) * | 2019-06-25 | 2019-09-27 | 河海大学 | 考虑电/热柔性负荷的区域综合能源系统储能配置方法 |
CN110929938B (zh) * | 2019-11-25 | 2022-09-27 | 新奥数能科技有限公司 | 能源系统的优化方法及装置、存储介质、电子装置 |
CN110929938A (zh) * | 2019-11-25 | 2020-03-27 | 新奥数能科技有限公司 | 能源系统的优化方法及装置、存储介质、电子装置 |
CN111179110A (zh) * | 2019-12-06 | 2020-05-19 | 清华-伯克利深圳学院筹备办公室 | 虚拟电厂可变阶聚合等效鲁棒动态模型建模方法及装置 |
CN111179110B (zh) * | 2019-12-06 | 2023-09-05 | 清华-伯克利深圳学院筹备办公室 | 虚拟电厂可变阶聚合等效鲁棒动态模型建模方法及装置 |
CN110808615A (zh) * | 2019-12-07 | 2020-02-18 | 国家电网有限公司 | 考虑不确定性的气电虚拟电厂调度优化方法 |
CN111210119A (zh) * | 2019-12-24 | 2020-05-29 | 国网山东省电力公司枣庄供电公司 | 多种市场下vpp电、热、气优化调度模型的建立方法 |
CN111382939A (zh) * | 2020-03-06 | 2020-07-07 | 国网冀北电力有限公司 | 虚拟电厂资源优化配置方法、装置及设备 |
CN111625961A (zh) * | 2020-05-26 | 2020-09-04 | 中国科学院工程热物理研究所 | 综合能源系统协同优化运行调控方法 |
CN111953018A (zh) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | 一种分布式多能互补供能系统及功率分配方法 |
CN112084641A (zh) * | 2020-08-28 | 2020-12-15 | 国网综合能源服务集团有限公司 | 一种基于itsdp模型的园区虚拟电厂优化运行方法 |
CN116937551A (zh) * | 2023-07-19 | 2023-10-24 | 太原理工大学 | 一种电气互联虚拟电厂优化调度方法及终端 |
CN116937551B (zh) * | 2023-07-19 | 2024-04-05 | 太原理工大学 | 一种电-气互联虚拟电厂优化调度方法及终端 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107506851A (zh) | 一种多区域虚拟电厂综合能源协调调度优化模型 | |
Li et al. | Two-stage optimal operation of integrated energy system considering multiple uncertainties and integrated demand response | |
Zhang et al. | A two-stage operation optimization method of integrated energy systems with demand response and energy storage | |
Luo et al. | Multi-objective capacity optimization of a distributed energy system considering economy, environment and energy | |
Song et al. | Multi-objective optimization of a solar hybrid CCHP system based on different operation modes | |
Pan et al. | Multi-objective and two-stage optimization study of integrated energy systems considering P2G and integrated demand responses | |
Shen et al. | Multi-objective capacity configuration optimization of an integrated energy system considering economy and environment with harvest heat | |
Ju et al. | A two-stage optimal coordinated scheduling strategy for micro energy grid integrating intermittent renewable energy sources considering multi-energy flexible conversion | |
CN107807523A (zh) | 考虑分时电价的区域能源互联网多源协调优化运行策略 | |
Wang et al. | Analysis of inlet air throttling operation method for gas turbine in performance of CCHP system under different operation strategies | |
CN109146182A (zh) | 计及多种储能的分布式冷热电联供系统的经济调度方法 | |
CN103400042A (zh) | 一种冷热电联供型微电网优化配置方法 | |
Li et al. | Optimization configuration of regional integrated energy system based on standard module | |
CN114202176A (zh) | 一种综合能源系统优化调度方法、系统、设备及存储介质 | |
CN108665188A (zh) | 一种基于优化模型的园区多能源主体综合配比方法 | |
CN108008629A (zh) | 一种多种能源互补利用系统的优化运行方法 | |
CN108052722A (zh) | 面向综合能效最优的分布式冷热电混合能源系统设计方法 | |
CN113435095A (zh) | 一种用于综合能源系统的优化调度的方法及系统 | |
Yu et al. | Cross‐regional integrated energy system scheduling optimization model considering conditional value at risk | |
Liu et al. | Capacity allocation for regional integrated energy system considering typical day economic operation | |
Liu | Research on optimal placement of low‐carbon equipment capacity in integrated energy system considering carbon emission and carbon trading | |
Wang et al. | Research on capacity and strategy optimization of combined cooling, heating and power systems with solar photovoltaic and multiple energy storage | |
Cao et al. | Optimal operation of cold–heat–electricity multi-energy collaborative system based on price demand response | |
CN115689166A (zh) | 一种区域分布式能源资源聚合利用方法及系统 | |
Zhao et al. | Multi-objective bi-level planning of the integrated energy system considering uncertain user loads and carbon emission during the equipment manufacturing process |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171222 |