CN111668836B - 一种源-网-荷-储海量设备的聚合与分层调控方法 - Google Patents

一种源-网-荷-储海量设备的聚合与分层调控方法 Download PDF

Info

Publication number
CN111668836B
CN111668836B CN202010582962.8A CN202010582962A CN111668836B CN 111668836 B CN111668836 B CN 111668836B CN 202010582962 A CN202010582962 A CN 202010582962A CN 111668836 B CN111668836 B CN 111668836B
Authority
CN
China
Prior art keywords
flexgas
aggregator
gas
nth
time
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.)
Active
Application number
CN202010582962.8A
Other languages
English (en)
Other versions
CN111668836A (zh
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.)
Huazhong University of Science and Technology
Guangdong Power Grid Co Ltd
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202010582962.8A priority Critical patent/CN111668836B/zh
Publication of CN111668836A publication Critical patent/CN111668836A/zh
Application granted granted Critical
Publication of CN111668836B publication Critical patent/CN111668836B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • 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
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28—Arrangements for balancing of the load in a network by storage of energy
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects
    • 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
    • Y02E70/00—Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Power Engineering (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种源‑网‑荷‑储海量设备的聚合与分层调控方法，针对海量FlexGas用户广泛在接入气电联合系统的调控需求，建立了FlexGas负荷的气‑电能量转换单体模型，并考虑到海量FlexGas负荷的异质性，进一步进行聚合，得到海量FlexGas负荷的聚合模型，然后分别基于所得单体模型和聚合模型，建立了上层的考虑风电不确定性的上层调度模型和下层控制模型，共同构成了基于海量FlexGas负荷的源网荷储协调调度模型，从而实现对源网荷储进行协同调控，不仅能够消纳风电，提高气电联合系统的经济性，还可以实现对海量FlexGas负荷的有效控制。

Description

一种源-网-荷-储海量设备的聚合与分层调控方法
技术领域
本发明属于电气工程领域，更具体地，涉及一种源-网-荷-储海量设备的聚合与分层调控方法。
背景技术
风力发电技术的快速发展，对电力系统的运行提出了一系列新的挑战。电制气(P2G)技术由于可以利用电能生产天然气并将其储能在天然气网络中，为消纳风电提供了有效的解决方案，并进一步加强了电网与天然气系统的耦合。与此同时，气电联合系统调度中的需求响应也引起了越来越多的兴趣。FlexGas是一种横跨电力，天然气和供暖系统的新兴技术，它可以实现用户在使用能量间的转换，并且在需求响应方面具有很大的潜力。但是由于偏好和生活方式的广泛差异，海量的FlexGas用户之间的用能行为可能会存在很大的差异，故研究一种源-网-荷-储海量设备的聚合与分层调控方法存在重要的意义。
现有的源-网-荷-储海量设备的聚合与分层调控方法主要采用“采集+集中”的模式，然而调度中心直接控制大量异构柔性负载是不现实的，并且该方法未考虑FlexGas负荷为气电联合系统提供的需求响应，弃风成本较高，无法满足海量用户广泛接入气电联合系统的调控需求，故亟需同时兼容现有调度体系和创新市场机制的海量终端设备协同调控技术体系。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供一种源-网-荷-储海量设备的聚合与分层调控方法，用以解决现有技术由于未考虑FlexGas负荷而导致的弃风成本较高的技术问题。
为了实现上述目的，第一方面，本发明提供了一种源-网-荷-储海量设备的聚合与分层调控方法，包括以下步骤：
S1、根据FlexGas用户的热需求与其燃气锅炉的供热功率和电锅炉的供热功率之间的供需平衡关系，分别为各FlexGas聚合商所控制的各FlexGas用户的供热设备建立FlexGas负荷的气-电能量转换模型；
S2、分别将各FlexGas聚合商所控制的各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型进行聚合，为各FlexGas聚合商建立能量连续的海量FlexGas负荷的聚合模型；
S3、通过最小化气电联合系统运行成本，且将风电因素引入到约束条件，结合各FlexGas聚合商的海量FlexGas负荷的聚合模型，建立上层调度模型；通过最小化调度中心向FlexGas聚合商下发的电量调度需求与FlexGas聚合商实际的电量聚合结果之差，结合各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型，为各FlexGas聚合商建立下层控制模型；综合上层调度模型和各FlexGas聚合商的下层控制模型得到源网荷储的协调调控模型；
S4、采集FlexGas用户燃气锅炉和电锅炉的转换效率；并分别计算各FlexGas聚合商所控制的各FlexGas用户的热需求以及各FlexGas聚合商总的热需求；将FlexGas用户燃气锅炉和电锅炉的转换效率、各FlexGas聚合商所控制的各FlexGas用户的热需求、各FlexGas聚合商总的热需求以及气电联合系统的参数，输入到源网荷储的协调调控模型中进行求解，得到各FlexGas聚合商所控制的每一个FlexGas用户设备的调控结果。
进一步优选地，第n个FlexGas聚合商所控制的第m个FlexGas用户的FlexGas负荷的气-电能量转换模型为：
Figure BDA0002553101720000031
Figure BDA0002553101720000032
Figure BDA0002553101720000033
Figure BDA0002553101720000034
其中，
Figure BDA0002553101720000035
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求，
Figure BDA0002553101720000036
和
Figure BDA0002553101720000037
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的供热功率和电锅炉的供热功率，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure BDA0002553101720000038
和
Figure BDA0002553101720000039
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，
Figure BDA00025531017200000310
和
Figure BDA00025531017200000311
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉供热功率和电锅炉供热功率的最大值，
Figure BDA00025531017200000312
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用燃气锅炉供热，
Figure BDA00025531017200000313
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用电锅炉供热。
进一步优选地，第n个FlexGas聚合商的海量FlexGas负荷的聚合模型为：
Figure BDA00025531017200000314
Figure BDA00025531017200000315
Figure BDA0002553101720000041
其中，
Figure BDA0002553101720000042
为t时刻第n个FlexGas聚合商总的热需求，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure BDA0002553101720000043
和
Figure BDA0002553101720000044
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure BDA0002553101720000045
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，cd为电能的折扣价格，
Figure BDA0002553101720000046
和
Figure BDA0002553101720000047
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure BDA0002553101720000048
为t时刻第n个FlexGas聚合商燃气锅炉总能耗值的最大值，
Figure BDA0002553101720000049
为t时刻第n个FlexGas聚合商电锅炉总能耗值的最大值。
进一步优选地，上层调度模型的目标函数为：
Figure BDA00025531017200000410
其中，NT为调度的总时段数，
Figure BDA00025531017200000411
为t时刻火电机组运行成本，
Figure BDA00025531017200000412
为t时刻天然气的消耗成本，
Figure BDA00025531017200000413
为t时刻火电机组、P2G机组以及GFU机组的向上备用容量总成本，
Figure BDA00025531017200000414
为t时刻火电机组、P2G机组以及GFU机组的向下备用容量总成本，Nagg为FlexGas聚合商的个数，
Figure BDA00025531017200000415
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，ps为风电场景s的概值，Ns为风电预测场景个数，Nb为电力系统中交流电网的节点个数，cl和cw分别为切负荷和弃风的惩罚价格，
Figure BDA00025531017200000416
为t时刻在第s个风电预测场景下第i个交流电网节点处的切负荷量，
Figure BDA00025531017200000417
为t时刻在第s个风电预测场景下第i个交流电网节点处的弃风量，
Figure BDA0002553101720000051
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向上备用部署总成本，
Figure BDA0002553101720000052
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向下备用部署总成本。
进一步优选地，上层调度模型的约束包括：
Figure BDA0002553101720000053
Figure BDA0002553101720000054
Figure BDA0002553101720000055
其中，cd为电能的折扣价格，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure BDA0002553101720000056
和
Figure BDA0002553101720000057
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure BDA0002553101720000058
为t时刻第n个FlexGas聚合商燃气锅炉总能耗值的最大值，
Figure BDA0002553101720000059
为t时刻第n个FlexGas聚合商电锅炉总能耗值的最大值。
进一步地，第n个FlexGas聚合商的下层控制模型为：
Figure BDA00025531017200000510
Figure BDA00025531017200000511
Figure BDA00025531017200000512
Figure BDA00025531017200000513
Figure BDA0002553101720000061
Figure BDA0002553101720000062
其中，Nk为第n个FlexGas聚合商所控制的FlaxGas用户个数，NT为调度的总时段数，
Figure BDA0002553101720000063
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的耗电量，
Figure BDA0002553101720000064
为t时刻第n个FlexGas聚合商电锅炉总能耗值，
Figure BDA0002553101720000065
和
Figure BDA0002553101720000066
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的供热功率和电锅炉的供热功率，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure BDA0002553101720000067
和
Figure BDA0002553101720000068
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，
Figure BDA0002553101720000069
和
Figure BDA00025531017200000610
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉供热功率和电锅炉供热功率的最大值，
Figure BDA00025531017200000611
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用燃气锅炉供热，
Figure BDA00025531017200000612
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用电锅炉供热，
Figure BDA00025531017200000613
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求。
第二方面，本发明还提供了一种存储介质，当计算机读取所述指令时，使所述计算机执行本发明第一方面所提供的源-网-荷-储海量设备的聚合与分层调控方法。
总体而言，通过本总体而言，通过本发明所构思的以上技术方案，能够取得以下有益效果：
1、本发明提出了一种源-网-荷-储海量设备的聚合与分层调控方法，针对海量FlexGas用户广泛在接入气电联合系统的调控需求，建立了FlexGas负荷的气-电能量转换单体模型，并考虑到海量FlexGas负荷的异质性，进一步进行聚合，得到海量FlexGas负荷的聚合模型，然后分别基于所得单体模型和聚合模型，建立了上层的考虑风电不确定性的上层调度模型和下层控制模型，共同构成了基于海量FlexGas负荷的源网荷储协调调度模型，该调度模型通过对源网荷储进行协同调控，能够消纳风电，降低电联合系统的弃风成本，大大提高了气电联合系统的经济性。
2、本发明所提出的一种源-网-荷-储海量设备的聚合与分层调控方法，通过将海量的FlexGas负荷的气-电能量转换模型进行聚合，并基于上层调度模型确定气电联合系统运行成本最小时各FlexGas聚合商电锅炉的总能耗值，然后通过下层控制模型，确定各FlexGas聚合商所控制的各调度时段各FlexGas聚合商所控制的每一个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，实现对海量FlexGas负荷的有效控制。
附图说明
图1为本发明实施例1提供的源-网-荷-储海量设备的聚合与分层调控方法流程图；
图2为本发明实施例1提供的IEEE39节点电力系统和27节燃气网结构图；
图3为本发明实施例1提供的海量FlexGas负荷聚合后的调控结果图；
图4为本发明实施例1提供的分别采用不考虑FlexGas负荷的调控方法和本发明所提出的方法进行调控后的弃风成本结果图；
图5为本发明实施例1提供的采用本发明所提出的调控方法的聚合误差示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实施例1、
一种源-网-荷-储海量设备的聚合与分层调控方法，如图1所示，包括以下步骤：
S1、根据FlexGas用户的热需求与其燃气锅炉的供热功率和电锅炉的供热功率之间的供需平衡关系，分别为各FlexGas聚合商所控制的各FlexGas用户的供热设备建立FlexGas负荷的气-电能量转换模型；
FlexGas技术通过将电锅炉安装在燃气锅炉的热水回流管道上实现供暖用户在用电或者用气间的切换。当电锅炉打开时，回流的热水达到设置的温度，并向燃气锅炉传递“停止运行”的信号，然后燃气锅炉停止供暖。在配置有FlexGas装置的房屋中，暖气可由燃气锅炉消耗燃气或电锅炉消耗电产生；具体的，第n个FlexGas聚合商所控制的第m个FlexGas用户的FlexGas负荷的气-电能量转换模型为：
Figure BDA0002553101720000081
其中，
Figure BDA0002553101720000082
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求，t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的供热功率
Figure BDA0002553101720000083
和电锅炉的供热功率
Figure BDA0002553101720000084
满足：
Figure BDA0002553101720000085
其中，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure BDA0002553101720000086
和
Figure BDA0002553101720000087
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，
Figure BDA0002553101720000091
和
Figure BDA0002553101720000092
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉供热功率和电锅炉供热功率的最大值，
Figure BDA0002553101720000093
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用燃气锅炉供热，
Figure BDA0002553101720000094
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用电锅炉供热。
S2、分别将各FlexGas聚合商所控制的各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型进行聚合，为各FlexGas聚合商建立能量连续的海量FlexGas负荷的聚合模型；
具体的，对于装设有FlexGas装置的房屋，燃气锅炉和电锅炉的能耗值等于零或供暖需求，但是，对于包含了海量FlexGas用户的FlexGas聚合商，总用电量或者用气量可以实现从零到总供暖需求值间连的续变换，因此，基于FlexGas负荷的气-电能量转换模型及其大量性和异质性进行聚合，得到第n个FlexGas聚合商的海量FlexGas负荷的聚合模型为：
Figure BDA0002553101720000095
Figure BDA0002553101720000096
其中，
Figure BDA0002553101720000097
为t时刻第n个FlexGas聚合商总的热需求，
Figure BDA0002553101720000098
和
Figure BDA0002553101720000099
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure BDA00025531017200000910
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，cd为电能的折扣价格。进一步地，
Figure BDA00025531017200000911
Nk为每个FlexGas聚合商控制的FlaxGas用户个数。
具体的，t时刻第n个FlexGas聚合商燃气锅炉总能耗值
Figure BDA00025531017200000912
和电锅炉总能耗值
Figure BDA0002553101720000101
满足：
Figure BDA0002553101720000102
其中，
Figure BDA0002553101720000103
为t时刻第n个FlexGas聚合商燃气锅炉总能耗值的最大值，
Figure BDA0002553101720000104
为t时刻第n个FlexGas聚合商电锅炉总能耗值的最大值。与考虑海量FlexGas负荷的气-电能量转换模型参与气电联合系统的集中调度相比，聚合模型中消除了数百个二进制变量，可以大大减少调度中心的计算负担。
S3、通过最小化气电联合系统运行成本，且将风电因素引入到约束条件，结合各FlexGas聚合商的海量FlexGas负荷的聚合模型，建立上层调度模型；通过最小化调度中心向FlexGas聚合商下发的电量调度需求与FlexGas聚合商实际的电量聚合结果之差，结合各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型，为各FlexGas聚合商建立下层控制模型；综合上层调度模型和各FlexGas聚合商的下层控制模型得到源网荷储的协调调控模型；
具体的，为了最小化气电联合系统的运行成本，以气电联合系统运行成本最小为上层目标。由于电力系统和燃气网络通过P2G机组，GFU机组和FlexGas聚合商进行相互连接和耦合，为了应对风电的不确定性，随机规划考虑了电力系统和天然气网络的安全约束，以共同优化气电联合系统的发电量以及火电机组、P2G机组和GFU机组的备用容量，得到上层调度模型的目标函数为：
Figure BDA0002553101720000105
其中，NT为调度的总时段数，
Figure BDA0002553101720000106
为t时刻火电机组运行成本，
Figure BDA0002553101720000107
为t时刻天然气的消耗成本，
Figure BDA0002553101720000111
为t时刻火电机组、P2G机组以及GFU机组的向上备用容量总成本，
Figure BDA0002553101720000112
为t时刻火电机组、P2G机组以及GFU机组的向下备用容量总成本，Nagg为FlexGas聚合商的个数，
Figure BDA0002553101720000113
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，ps为风电场景s的概值，Ns为风电预测场景个数，Nb为电力系统中交流电网的节点个数，cl和cw分别为切负荷和弃风的惩罚价格，
Figure BDA0002553101720000114
为t时刻在第s个风电预测场景下第i个交流电网节点处的切负荷量，
Figure BDA0002553101720000115
为t时刻在第s个风电预测场景下第i个交流电网节点处的弃风量，
Figure BDA0002553101720000116
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向上备用部署总成本，
Figure BDA0002553101720000117
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向下备用部署总成本；
具体的，
Figure BDA0002553101720000118
的计算方式如下：
Figure BDA0002553101720000119
Figure BDA00025531017200001110
Figure BDA00025531017200001111
Figure BDA00025531017200001112
Figure BDA00025531017200001113
Figure BDA00025531017200001114
其中，Nb为电力系统中交流电网的节点个数，
Figure BDA00025531017200001115
为t时刻第i个交流电网节点处火电机组的开启价格，ui,t表示火电机组是否处于正在开启状态，
Figure BDA0002553101720000121
为t时刻第i个交流电网节点处火电机组的关闭价格，vi,t用于表示火电机组是否处于正在关闭的状态，
Figure BDA0002553101720000122
为t时刻第i个交流电网节点处火电机组的煤耗成本，Ngd为电力系统中气网的负荷节点个数，cg为燃气价格，
Figure BDA0002553101720000123
为t时刻电力系统中第k个气网负荷节点处的能耗值，
Figure BDA0002553101720000124
为t时刻P2G机组在第k个气网负荷节点处产生的气体质量流量，
Figure BDA0002553101720000125
为t时刻FlexGas聚合商在第k个气网负荷节点处消耗的气体质量流量，
Figure BDA0002553101720000126
分别为t时刻第i个交流电网节点处火电机组、GFU机组、P2G机组的上备用价格，
Figure BDA0002553101720000127
分别为t时刻第i个交流电网节点处火电机组、GFU机组、P2G机组的上备用容量，
Figure BDA0002553101720000128
分别为t时刻第i个交流电网节点处火电机组、GFU机组、P2G机组的下备用价格，
Figure BDA0002553101720000129
分别为t时刻第i个交流电网节点处火电机组、GFU机组、P2G机组的下备用容量，
Figure BDA00025531017200001210
分别为t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、GFU机组、P2G机组的上备用价格，
Figure BDA00025531017200001211
Figure BDA00025531017200001212
分别为t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、GFU机组、P2G机组的上备用容量，
Figure BDA00025531017200001213
分别为t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、GFU机组、P2G机组的下备用价格，
Figure BDA00025531017200001214
分别为t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、GFU机组、P2G机组的下备用容量。
进一步地，上层调度模型的目标为两个阶段的运行成本之和最小，第一阶段的运行成本包括火电机组运行成本、天然气消耗成本、向上和向下的备用容量成本，以及FlexGas聚合商的折扣成本；第二阶段的运行成本包括在各风电场景下的弃风和切负荷成本，以及在各风电场景里下火电机组、P2G机组、GFU机组向上备用和向下备用的部署成本。调度中心将风电预测场景作为输入，并通过FlexGas聚合商提供的运行参数，然后根据两阶段随机优化制定气电联合系统的最优调度方案。
进一步地，第一阶段表示日前调度，其约束条件包括：
(1.1)电力系统的功率平衡约束：
Figure BDA0002553101720000131
其中，Nb为电力系统中交流电网的节点个数，Pi,t表示t时刻第i个交流电网节点处火电机组的出力，
Figure BDA0002553101720000132
表示t时刻第i个交流电网节点处GFU机组的出力，Wi,t为t时刻第i个交流电网节点处的风电出力，
Figure BDA0002553101720000133
为t时刻第i个交流电网节点处GFU机组的出力，
Figure BDA0002553101720000134
为t时刻第i个交流电网节点处FlaxGas消耗的电量，Li,t为t时刻第i个交流电网节点处的负荷值。
(1.2)火电机组最小开/关机时间约束：
Figure BDA0002553101720000135
其中，Ii,t为t时刻第i个交流电网节点处火电机组的启停状态，
Figure BDA0002553101720000136
和
Figure BDA0002553101720000137
分别表示t时刻第i个交流电网节点处火电机组的最小开机和关机时间。
(1.3)火电机组燃料成本的松弛约束：
Figure BDA0002553101720000138
其中，
Figure BDA0002553101720000139
为t时刻第i个交流电网节点处火电机组的燃料成本，火电机组的二次煤耗成本可进行分段线性化，ai,r、bi,r分别第i个交流电网节点处煤耗成本函数分段线性化后的第r段的斜率和截距。
(1.4)火电机组、GFU机组和P2G机组的启停状态约束：
Figure BDA0002553101720000141
其中，
Figure BDA0002553101720000142
用于表示机组是否处于正在开启状态，
Figure BDA0002553101720000143
用于表示机组是否处于正在关闭的状态，
Figure BDA0002553101720000144
用于表示t时刻第i个交流电网节点处机组的启停状态；其中，“.”为通配符，用于指代火电机组、GFU机组或P2G机组，这里的机组为火电机组、GFU机组或P2G机组，每种机组均需要满足上述约束。
(1.5)火电机组、GFU机组和P2G机组的出力范围约束：
Figure BDA0002553101720000145
其中，
Figure BDA0002553101720000146
用于表示t时刻第i个交流电网节点处机组的出力，
Figure BDA0002553101720000147
用于表示第i个交流电网节点处机组出力的下限，
Figure BDA0002553101720000148
用于表示第i个交流电网节点处机组出力的上限；其中，“.”为通配符，用于指代火电机组、GFU机组或P2G机组，这里的机组为火电机组、GFU机组或P2G机组，每种机组均需要满足上述约束。
(1.6)火电机组、GFU机组和P2G机组的爬坡约束：
Figure BDA0002553101720000149
其中，
Figure BDA00025531017200001410
用于表示第i个交流电网节点处机组向上爬坡的上限，
Figure BDA00025531017200001411
用于表示第i个交流电网节点处机组向下爬坡的上限；其中，“.”为通配符，用于指代火电机组、GFU机组或P2G机组，这里的机组为火电机组、GFU机组或P2G机组，每种机组均需要满足上述约束。
(1.7)火电机组和GFU机组预留的上备用和下备用容量约束：
Figure BDA00025531017200001412
其中，
Figure BDA0002553101720000151
用于表示t时刻第i个交流电网节点处机组提供的上备用容量，
Figure BDA0002553101720000152
用于表示t时刻第i个交流电网节点处机组提供的下备用容量，
Figure BDA0002553101720000153
用于表示第i个交流电网节点处机组向上爬坡的上限，
Figure BDA0002553101720000154
用于表示第i个交流电网节点处机组向下爬坡的上限，
Figure BDA0002553101720000155
用于表示t时刻第i个交流电网节点处机组的启停状态，
Figure BDA0002553101720000156
用于表示第i个交流电网节点处机组出力的上限，
Figure BDA0002553101720000157
用于表示t时刻第i个交流电网节点处机组的出力，
Figure BDA0002553101720000158
用于表示第i个交流电网节点处机组出力的下限；其中，“.”为通配符，用于指代火电机组或GFU机组，这里的机组为火电机组或GFU机组，每种机组均需要满足上述约束。
(1.8)P2G机组预留的上备用和下备用容量约束：
Figure BDA0002553101720000159
其中，
Figure BDA00025531017200001510
为t时刻第i个交流电网节点处P2G机组的上备用容量，
Figure BDA00025531017200001511
为t时刻第i个交流电网节点处P2G机组的下备用容量，
Figure BDA00025531017200001512
为第i个交流电网节点处P2G机组向上爬坡的上限，
Figure BDA00025531017200001513
为第i个交流电网节点处P2G机组向下爬坡的上限，
Figure BDA00025531017200001514
为t时刻第i个交流电网节点处P2G机组的启停状态，
Figure BDA00025531017200001515
为t时刻第i个交流电网节点处P2G机组的出力，
Figure BDA00025531017200001516
用于表示第i个交流电网节点处P2G机组出力的下限，
Figure BDA00025531017200001517
用于表示第i个交流电网节点处P2G机组出力的上限。
(1.9)传输线的功率上下限约束：
Figure BDA00025531017200001518
其中，sfl,i为线路l相对第i个交流电网的灵敏度矩阵，Pi,t表示t时刻第i个交流电网节点处火电机组的出力，
Figure BDA0002553101720000161
表示t时刻第i个交流电网节点处GFU机组的出力，Wi,t为t时刻第i个交流电网节点处的风电出力，
Figure BDA0002553101720000162
为t时刻第i个交流电网节点处P2G机组的出力，
Figure BDA0002553101720000163
为t时刻第i个交流电网节点处FlaxGas消耗的电量，Li,t为t时刻第i个交流电网节点处的负荷值，
Figure BDA0002553101720000164
为电力系统交流电网线路长期运行允许的潮流上限。
(1.10)燃气网络的气节点上的气源或气负荷的质量流量约束：
Figure BDA0002553101720000165
其中，Ngd为电力系统中气网的负荷节点个数，
Figure BDA0002553101720000166
表示t时刻第k个气网负荷节点处的负荷，
Figure BDA0002553101720000167
为t时刻FlexGas聚合商在第k个气网负荷节点处消耗的气体质量流量，
Figure BDA0002553101720000168
为t时刻第k个气网负荷节点处P2G机组的能量，
Figure BDA0002553101720000169
为t时刻第k个气网负荷节点处GFU机组的能量，其中，
Figure BDA00025531017200001610
为第i个交流电网节点处的P2G机组的能量转换效率，
Figure BDA00025531017200001611
为第i个交流电网节点处的GFU机组的能量转换效率。需要说明的是，气电联合系统包含气网和电网，火电机组和电力负荷属于电网，气负荷属于气网，GFU，P2G和FlexGas聚合商均连接在电网和气网，这里，符号i表示连接的电网节点，符号k表示连接的气网节点。
进一步地，第二阶段为对应在不同风电场景下的再调度，其约束条件包括：
(2.1)针对风电场景的约束，具体包括：
(2.1.1)各风电场景下的功率再平衡约束：
Figure BDA00025531017200001612
其中，Nb为电力系统中交流电网的节点个数，Pi,t,s、
Figure BDA00025531017200001613
分别为t时刻在第s个风电预测场景下第i个交流电网节点处的火电机组、GFU机组、P2G机组的出力，Wi,t,s为t时刻在第s个风电预测场景下第i个交流电网节点处的风电出力，Li,t,s为t时刻在第s个风电预测场景下第i个交流电网节点处的负荷值，
Figure BDA0002553101720000171
为t时刻在第s个风电预测场景下第i个交流电网节点处的弃风量，
Figure BDA0002553101720000172
为t时刻在第s个风电预测场景下第i个交流电网节点处的切负荷量，
Figure BDA0002553101720000173
为t时刻第i个交流电网节点处FlexGas聚合商消耗的电量。
(2.1.2)火电机组、P2G机组和GFU机组在每个风电场景下的部署备用约束：
Figure BDA0002553101720000174
其中，
Figure BDA0002553101720000175
用于表示t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、P2G机组或GFU机组提供的上备用容量，
Figure BDA0002553101720000176
用于表示t时刻在第s个风电预测场景下第i个交流电网节点处火电机组、P2G机组或GFU机组提供的下备用容量，
Figure BDA0002553101720000177
用于表示t时刻第i个交流电网节点处火电机组、P2G机组或GFU机组提供的上备用容量，
Figure BDA0002553101720000178
用于表示t时刻第i个交流电网节点处火电机组、P2G机组或GFU机组提供的下备用容量，其中“.”为通配符，用于指代火电机组、GFU机组或P2G机组，每种机组均需要满足上述约束。
(2.1.3)在第s个风电预测场景下，火电机组和GFU机组的出力需满足爬坡约束：
Figure BDA0002553101720000179
其中，
Figure BDA0002553101720000181
用于表示t时刻在第s个风电预测场景下第i个交流电网节点处机组的出力，
Figure BDA0002553101720000182
表示t时刻第i个交流电网节点处机组的出力，
Figure BDA0002553101720000183
为t时刻在第s个风电预测场景下第i个交流电网节点处机组的上备用容量，
Figure BDA0002553101720000184
为t时刻在第s个风电预测场景下第i个交流电网节点处机组的下备用容量，
Figure BDA0002553101720000185
用于表示第i个交流电网节点处机组向上爬坡的上限，
Figure BDA0002553101720000186
用于表示第i个交流电网节点处机组向下爬坡的上限；其中，“.”为通配符，用于指代火电机组或GFU机组，这里的机组为火电机组或GFU机组，每种机组均需要满足上述约束。
(2.1.4)在第s个风电预测场景下，P2G机组的出力需满足爬坡约束：
Figure BDA0002553101720000187
其中，
Figure BDA0002553101720000188
用于表示t时刻在第s个风电预测场景下第i个交流电网节点处P2G机组的出力，
Figure BDA0002553101720000189
表示t时刻第i个交流电网节点处P2G机组的出力，
Figure BDA00025531017200001810
为t时刻在第s个风电预测场景下第i个交流电网节点处P2G机组的上备用容量，
Figure BDA00025531017200001811
为t时刻在第s个风电预测场景下第i个交流电网节点处P2G机组的下备用容量。
(2.1.5)电力系统输电线路的传输功率约束：
Figure BDA00025531017200001812
(2.1.6)弃风约束：
Figure BDA00025531017200001813
(2.1.7)切负荷约束：
Figure BDA00025531017200001814
(2.2)针对燃气网络的约束，具体包括：
(2.2.1)燃气的动量传输过程约束：
Figure BDA0002553101720000191
其中，MFab,t+1和MEab,t+1分别表示t+1时刻燃气管道首端Fab和末端Eab的气体质量流量，Aab、Lab、dab分别表示燃气管道ab的横截面积、长度和直径，pa,t+1和pb,t+1分别表示t+1时刻燃气管道两端节点a、b的气压，Δt为调度时间间隔，λ为摩擦系数，
Figure BDA0002553101720000192
为燃气管道ab中的气体平均流速。
(2.2.2)燃气管道中的物质平衡过程：
Figure BDA0002553101720000193
其中，c为音速。
(2.2.3)进出燃气管道的质量流量之间的平衡约束：
Figure BDA0002553101720000194
其中，a∈(.)c表示以节点c为末端的燃气管道上的节点，b∈c(.)表示以节点c为首端的燃气管道上的节点，MEac,t和MFcb,t分别表示t时刻燃气管道首端Fac和末端Ecb的气体质量流量，Mc,t为t时刻节点c处的质量流量。需要说明的是，气网有包括道和节点，管道的两端分别是两个节点，对于节点c来说，连接节点c的所有管道流向c的质量流量要等于节点c出产生或者消耗的质量流量。
(2.2.4)燃气管道源节点处的气压约束：
psa,t＝psa,0
其中，psa,t为t时刻燃气管道源节点sa的气压，psa,0表示初始时刻燃气管道源节点sa的气压，源节点的气压在所有时刻保持不变。
(2.2.5)燃气管道节点气压范围约束：
Figure BDA0002553101720000201
其中，
Figure BDA0002553101720000202
为燃气管道节点a的气压下限，
Figure BDA0002553101720000203
为燃气管道节点a的气压上限。
(2.2.6)燃气管道质量和流量的范围约束：
Figure BDA0002553101720000204
其中，
Figure BDA0002553101720000205
为燃气管道首端气体质量流量下限，
Figure BDA0002553101720000206
为燃气管道首端气体质量流量上限，
Figure BDA0002553101720000207
为燃气管道末端气体质量流量下限，
Figure BDA0002553101720000208
为燃气管道末端气体质量流量上限。
(2.2.7)风电场景下再调度的气源或气负荷的实际质量流量约束：
Figure BDA0002553101720000209
其中，
Figure BDA00025531017200002010
为t时刻在第s个风电预测场景下第k个气网负荷节点处P2G机组的能量，
Figure BDA00025531017200002011
为t时刻在第s个风电预测场景下第k个气网负荷节点处GFU机组的能量，
Figure BDA00025531017200002012
为t时刻在第s个风电预测场景下第k个气网负荷节点处P2G机组的能量，
Figure BDA00025531017200002013
表示t时刻在第s个风电预测场景下第i个交流电网节点处GFU机组的出力，
Figure BDA00025531017200002014
为t时刻在第s个风电预测场景下第i个交流电网节点处P2G机组的出力。
除此之外，结合FlexGas聚合商的海量FlexGas负荷的聚合模型，上层调度模型的约束条件还包括：
Figure BDA00025531017200002015
Figure BDA0002553101720000211
Figure BDA0002553101720000212
进一步地，气电联合系统包含气网和电网，火电机组和电力负荷属于电网，气负荷属于气网，GFU，P2G和FlexGas聚合商均连接在电网和气网，每一个FlexGas聚合商均会与一个气网负荷节点以及一个交流电网节点相连，记第n个FlexGas聚合商与第k个气网负荷节点相连，则t时刻第n个FlexGas聚合商燃气锅炉总能耗值与其连接的第k个气网负荷节点处消耗的气体质量流量相等，即
Figure BDA0002553101720000213
同样的，记第n个FlexGas聚合商与第i个交流电网节点相连，则第n个FlexGas聚合商电锅炉总能耗值与其连接的第i个交流电网节点处FlaxGas聚合商消耗的电量相等，即
Figure BDA0002553101720000214
进一步地，下层目标与海量FlexGas用户的聚合和控制相关，假设每个FlexGas机组均安装了智能控制器并与FlexGas聚合商进行双向通信，则FlexGas聚合商不仅可以获取各个FlexGas用户的操作参数，还可以向其释放控制信号。以调度中心向FlexGas聚合商下发的电量调度需求(即最小化上层目标所得的结果——第n个FlexGas聚合商电锅炉总能耗值
Figure BDA0002553101720000215
)与FlexGas聚合商实际电量聚合结果之差最小为下层目标，即
Figure BDA0002553101720000216
进一步地，为了使
Figure BDA0002553101720000217
最小，即是使所有FlexGas用户的总耗电量最大(即
Figure BDA0002553101720000218
)。结合FlexGas负荷的气-电能量转换模型，得到第n个FlexGas聚合商的下层控制模型为：
Figure BDA0002553101720000219
Figure BDA0002553101720000221
Figure BDA0002553101720000222
Figure BDA0002553101720000223
Figure BDA0002553101720000224
Figure BDA0002553101720000225
其中，Nk为第n个FlexGas聚合商所控制的FlaxGas用户个数，NT为调度的总时段数，
Figure BDA0002553101720000226
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的耗电量，
Figure BDA0002553101720000227
为t时刻第n个FlexGas聚合商电锅炉总能耗值。通过求解上述模型可以得到第n个FlexGas聚合商所控制的各FlexGas用户的燃气锅炉的耗气量和电锅炉的耗电量。采用同样的方法，即可求解得到各FlexGas聚合商所控制的各FlexGas用户的燃气锅炉的耗气量和电锅炉的耗电量。
S4、采集FlexGas用户燃气锅炉和电锅炉的转换效率；并分别计算各FlexGas聚合商所控制的各FlexGas用户的热需求以及各FlexGas聚合商总的热需求；将FlexGas用户燃气锅炉和电锅炉的转换效率、各FlexGas聚合商所控制的各FlexGas用户的热需求、各FlexGas聚合商总的热需求以及气电联合系统的参数，输入到源网荷储的协调调控模型中进行求解，得到各FlexGas聚合商所控制的每一个FlexGas用户设备的调控结果。
具体的，根据热力学公式，可以得到t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求为：
Figure BDA0002553101720000231
其中，
Figure BDA0002553101720000232
和
Figure BDA0002553101720000233
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户家的室内温度和室外温度，Rn,m为第n个FlexGas聚合商所控制的第m个FlexGas用户房屋的等效热阻，Cn,m为第n个FlexGas聚合商所控制的第m个FlexGas用户房屋的等效热容。
Figure BDA0002553101720000234
用来表示通过功率调节第m个FlexGas用户家中t时刻室内温度的动态过程，由于室内温度最终将达到温度设定点，并且在一定的热功率下保持恒定，故在预测FlexGas用户的热量需求时该动态过程可以忽略不计。因此，第m个FlexGas用户t时刻的热需求可以由设定点温度和室外温度之间的差异确定，具体为：
Figure BDA0002553101720000235
进一步根据
Figure BDA0002553101720000236
即可计算得到第n个FlexGas聚合商总的热需求。
将各FlexGas聚合商总的热需求发送给调度中心，调度中心将FlexGas用户燃气锅炉和电锅炉的转换效率、所得各FlexGas聚合商总的热需求以及气电联合系统的参数，输入到源网荷储的协调调控模型的上层调度模型进行求解，得到各FlexGas聚合商电锅炉总能耗值，并发送给各对应的FlexGas聚合商，FlexGas聚合商将所得FlexGas聚合商电锅炉总能耗值以及气电联合系统的参数输入到源网荷储的协调调控模型的下层控制模型中，求解得到各FlexGas聚合商所控制的每一个FlexGas用户设备的调控结果，即各调度时段各FlexGas聚合商所控制的每一个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量。其中，气电联合系统的参数包括交流电网、燃气网络、火电机组、P2G机组、GFU机组、风电机组、参与需求响应的FlexGas负荷商聚合的参数。具体的，交流电网参数包括：交流电网节点个数Nb，线路长期运行允许潮流上限
Figure BDA0002553101720000241
燃气网络包括：负荷节点个数Ngd，各节点的气压上、下限
Figure BDA0002553101720000242
各管道的物质流量上、下限
Figure BDA0002553101720000243
火电机组、P2G机组、GFU机组参数包括：出力上、下限
Figure BDA0002553101720000244
最大上、下爬坡率
Figure BDA0002553101720000245
火电机组最小启、停时间
Figure BDA0002553101720000246
线性化煤耗成本的斜率ai,r和截距bi,r，P2G机组和GFU机组的能量转换效率
Figure BDA0002553101720000247
和
Figure BDA0002553101720000248
风电机参数包括：风电机组出力预测值Wi,t和场景值Wi,t,s。参与需求响应的FlexGas负荷商聚合的参数包括：FlexGas聚合商个数Nagg，每个FlexGas聚合商所控制的FlaxGas用户个数Nk。
为了详细说明本发明所提出的方法，下面结合实施例进行详述，如图2所示为所选取的IEEE39节点电力系统和27节燃气网结构图。在电力系统的母线17、22、32和33处分别安装了2个P2G机组和2个GFU机组，并分别对应于燃气网络中的节点27、8、13和22。每个P2G机组和GFU机组的容量分别为60MW和200MW。此外，该电力系统包括6个火电机组，2个风电场，以及21条负荷母线和29条输电线路。天然气系统包含2个源节点，6个负荷节点和27条管道。六个FlexGas聚合商分别连接在燃气网络的6个负荷节点和电力系统中的母线7、12、18、24、27和31上。每个FlexGas聚合商可控制的FlexGas用户数为500。切负荷和弃风的惩罚价格分别为200$/MWh和1500$/MWh。按照上述步骤S1-S4进行建模求解，具体的，采用MATLAB通过商业软件GUROBI对源网荷储的协调调控模型中进行求解的进行求解，得到海量FlexGas负荷聚合后的调控结果如图3所示。为了进一步说明本发明所提出的方法，分别采用不考虑FlexGas负荷的调控方法和本发明所提出的方法进行调控，得到图4所示的结果，图4表现了FlexGas负荷参与需求响应在气电联合系统中的风电消纳效果，从图中可以看出，在各个调度时段，相比于现有的不考虑FlexGas负荷的调控方法，本发明所提出方法，通过考虑海量FlexGas负荷，为气电联合系统提供需求响应，弃风成本更低。图5验证了分层调控下聚合模型的有效性，从图中可以看出，在各个调度时段，采用本发明所提出的源-网-荷-储海量设备的聚合与分层调控方法的聚合误差均为0，聚合效果良好。
实施例2、
一种存储介质，当计算机读取所述指令时，使所述计算机执行本发明实施例1所提供的源-网-荷-储海量设备的聚合与分层调控方法。具体过程同上，这里不做赘述。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，包括以下步骤：
S1、根据FlexGas用户的热需求与其燃气锅炉的供热功率和电锅炉的供热功率之间的供需平衡关系，分别为各FlexGas聚合商所控制的各FlexGas用户的供热设备建立FlexGas负荷的气-电能量转换模型；
S2、分别将各FlexGas聚合商所控制的各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型进行聚合，为各FlexGas聚合商建立能量连续的海量FlexGas负荷的聚合模型；
S3、通过最小化气电联合系统运行成本，且将风电因素引入到约束条件，结合各FlexGas聚合商的海量FlexGas负荷的聚合模型，建立上层调度模型；通过最小化调度中心向FlexGas聚合商下发的电量调度需求与FlexGas聚合商实际的电量聚合结果之差，结合各FlexGas用户供热设备的FlexGas负荷的气-电能量转换模型，为各FlexGas聚合商建立下层控制模型；综合上层调度模型和各FlexGas聚合商的下层控制模型得到源网荷储的协调调控模型；
S4、采集FlexGas用户燃气锅炉和电锅炉的转换效率；并分别计算各FlexGas聚合商所控制的各FlexGas用户的热需求以及各FlexGas聚合商总的热需求；将FlexGas用户燃气锅炉和电锅炉的转换效率、各FlexGas聚合商所控制的各FlexGas用户的热需求、各FlexGas聚合商总的热需求以及气电联合系统的参数，输入到源网荷储的协调调控模型中进行求解，得到各FlexGas聚合商所控制的每一个FlexGas用户设备的调控结果。
2.根据权利要求1所述的源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，第n个FlexGas聚合商所控制的第m个FlexGas用户的FlexGas负荷的气-电能量转换模型为：
Figure FDA0003147513110000021
Figure FDA0003147513110000022
Figure FDA0003147513110000023
Figure FDA0003147513110000024
其中，
Figure FDA0003147513110000025
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求，
Figure FDA0003147513110000026
和
Figure FDA0003147513110000027
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的供热功率和电锅炉的供热功率，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure FDA0003147513110000028
和
Figure FDA0003147513110000029
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，
Figure FDA00031475131100000210
和
Figure FDA00031475131100000211
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉供热功率和电锅炉供热功率的最大值，
Figure FDA00031475131100000212
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用燃气锅炉供热，
Figure FDA00031475131100000213
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用电锅炉供热。
3.根据权利要求1所述的源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，第n个FlexGas聚合商的海量FlexGas负荷的聚合模型为：
Figure FDA00031475131100000214
Figure FDA00031475131100000215
Figure FDA0003147513110000031
其中，
Figure FDA0003147513110000032
为t时刻第n个FlexGas聚合商总的热需求，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure FDA0003147513110000033
和
Figure FDA0003147513110000034
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure FDA0003147513110000035
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，cd为电能的折扣价格，
Figure FDA0003147513110000036
和
Figure FDA0003147513110000037
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure FDA0003147513110000038
为t时刻第n个FlexGas聚合商燃气锅炉总能耗值的最大值，
Figure FDA0003147513110000039
为t时刻第n个FlexGas聚合商电锅炉总能耗值的最大值。
4.根据权利要求1所述的源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，所述上层调度模型的目标函数为：
Figure FDA00031475131100000310
其中，NT为调度的总时段数，Ft th为t时刻火电机组运行成本，Ft gas为t时刻天然气的消耗成本，Ft ur为t时刻火电机组、P2G机组以及GFU机组的向上备用容量总成本，Ft dr为t时刻火电机组、P2G机组以及GFU机组的向下备用容量总成本，Nagg为FlexGas聚合商的个数，
Figure FDA00031475131100000311
为t时刻第n个FlexGas聚合商使用电能时的折扣成本，ps为风电场景s的概值，Ns为风电预测场景个数，Nb为电力系统中交流电网的节点个数，cl和cw分别为切负荷和弃风的惩罚价格，
Figure FDA00031475131100000312
为t时刻在第s个风电预测场景下第i个交流电网节点处的切负荷量，
Figure FDA0003147513110000041
为t时刻在第s个风电预测场景下第i个交流电网节点处的弃风量，
Figure FDA0003147513110000042
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向上备用部署总成本，
Figure FDA0003147513110000043
为t时刻在第s个风电预测场景下火电机组、P2G机组以及GFU机组的向下备用部署总成本。
5.根据权利要求4所述的源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，所述上层调度模型的约束包括：
Figure FDA0003147513110000044
Figure FDA0003147513110000045
Figure FDA0003147513110000046
其中，cd为电能的折扣价格，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure FDA0003147513110000047
和
Figure FDA0003147513110000048
分别为t时刻第n个FlexGas聚合商燃气锅炉总能耗值和电锅炉总能耗值，
Figure FDA0003147513110000049
为t时刻第n个FlexGas聚合商燃气锅炉总能耗值的最大值，
Figure FDA00031475131100000410
为t时刻第n个FlexGas聚合商电锅炉总能耗值的最大值。
6.根据权利要求1所述的源-网-荷-储海量设备的聚合与分层调控方法，其特征在于，第n个FlexGas聚合商的下层控制模型为：
Figure FDA00031475131100000411
Figure FDA00031475131100000412
Figure FDA0003147513110000051
Figure FDA00031475131100000515
Figure FDA0003147513110000052
Figure FDA0003147513110000053
其中，Nk为第n个FlexGas聚合商所控制的FlaxGas用户个数，NT为调度的总时段数，
Figure FDA0003147513110000054
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的耗电量，
Figure FDA0003147513110000055
为t时刻第n个FlexGas聚合商电锅炉总能耗值，
Figure FDA0003147513110000056
和
Figure FDA0003147513110000057
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的供热功率和电锅炉的供热功率，ηg-h和ηe-h分别为FlexGas用户的燃气锅炉和电锅炉的转换效率，
Figure FDA0003147513110000058
和
Figure FDA0003147513110000059
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉的耗气量和电锅炉的耗电量，
Figure FDA00031475131100000510
和
Figure FDA00031475131100000511
分别为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户燃气锅炉供热功率和电锅炉供热功率的最大值，
Figure FDA00031475131100000512
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用燃气锅炉供热，
Figure FDA00031475131100000513
表示t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户是否用电锅炉供热，
Figure FDA00031475131100000514
为t时刻第n个FlexGas聚合商所控制的第m个FlexGas用户的热需求。
7.一种存储介质，其特征在于，当计算机读取指令时，使所述计算机执行权利要求1-6任意一项所述的源-网-荷-储海量设备的聚合与分层调控方法。
CN202010582962.8A 2020-06-23 2020-06-23 一种源-网-荷-储海量设备的聚合与分层调控方法 Active CN111668836B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010582962.8A CN111668836B (zh) 2020-06-23 2020-06-23 一种源-网-荷-储海量设备的聚合与分层调控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010582962.8A CN111668836B (zh) 2020-06-23 2020-06-23 一种源-网-荷-储海量设备的聚合与分层调控方法

Publications (2)

Publication Number Publication Date
CN111668836A CN111668836A (zh) 2020-09-15
CN111668836B true CN111668836B (zh) 2021-08-31

Family

ID=72389507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010582962.8A Active CN111668836B (zh) 2020-06-23 2020-06-23 一种源-网-荷-储海量设备的聚合与分层调控方法

Country Status (1)

Country Link
CN (1) CN111668836B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549256A (zh) * 2018-04-08 2018-09-18 中国电力科学研究院有限公司 一种需求响应数字物理混合仿真方法和系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107784382A (zh) * 2016-08-31 2018-03-09 北京南瑞电研华源电力技术有限公司 基于能源路由器的用户侧能源互联网规划方法
CN109063925B (zh) * 2018-08-16 2021-08-17 合肥工业大学 一种计及负荷聚合商的区域综合能源系统优化运行方法
CN110619425B (zh) * 2019-08-06 2022-03-04 国网山东省电力公司经济技术研究院 一种考虑源网荷储差异特性的多功能区综合能源系统协同规划方法
CN110968827B (zh) * 2019-11-13 2023-06-27 国家电网有限公司 一种多区域综合能源系统优化配置方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549256A (zh) * 2018-04-08 2018-09-18 中国电力科学研究院有限公司 一种需求响应数字物理混合仿真方法和系统

Also Published As

Publication number Publication date
CN111668836A (zh) 2020-09-15

Similar Documents

Publication Publication Date Title
Lv et al. Model predictive control based robust scheduling of community integrated energy system with operational flexibility
Wu et al. Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers
CN104734168B (zh) 一种基于电热联合调度的微电网运行优化系统及方法
CN109919399B (zh) 一种综合能源系统日前经济调度方法及系统
CN111737884B (zh) 一种含多种清洁能源微能源网多目标随机规划方法
CN113725915A (zh) 一种考虑可再生能源不确定性和热惰性的乡村电热综合能源系统运行优化方法
Huang et al. A multi-rate dynamic energy flow analysis method for integrated electricity-gas-heat system with different time-scale
CN113256045A (zh) 考虑风光不确定性的园区综合能源系统日前经济调度方法
CN113792969A (zh) 考虑气网动态特性和电转气综合能源系统优化调度方法
CN116341881B (zh) 计及热网灵活性的电-热系统鲁棒超前调度方法及系统
CN112465236B (zh) 一种考虑综合满意度的社区综合能源系统调度方法
CN111985165B (zh) 一种计及智能楼宇热储特性的电热综合能源系统调度方法
Cheng et al. A day-ahead scheduling of large-scale thermostatically controlled loads model considering second-order equivalent thermal parameters model
Gu et al. Bi-level decentralized optimal economic dispatch for urban regional integrated energy system under carbon emission constraints
CN114139837A (zh) 一种考虑双层碳排放优化分配模型的区域多系统双层分散优化调度方法
Cui et al. Optimal operation of CCHP microgrids with multiple shiftable loads in different auxiliary heating source systems
CN112862158B (zh) 一种综合能源系统日前优化调度方法及系统
Chen et al. Flexible dispatching method for park‐level integrated energy systems considering energy cascade utilization
CN111047069B (zh) 一种考虑热网储能特性的热电联合调度灵活性评估方法
CN113885367A (zh) 一种楼宇能效管理控制系统及方法
CN110992206B (zh) 一种多源电场的优化调度方法及系统
CN111668836B (zh) 一种源-网-荷-储海量设备的聚合与分层调控方法
Yang et al. Dual-layer flexibility dispatching of distributed integrated energy systems incorporating resilient heating schemes based on the standardized thermal resistance method
Wang et al. Lyapunov optimization based online energy flow control for multi-energy community microgrids
CN113077173B (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
TR01 Transfer of patent right

Effective date of registration: 20230426

Address after: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037

Patentee after: HUAZHONG University OF SCIENCE AND TECHNOLOGY

Patentee after: GUANGDONG POWER GRID Co.,Ltd.

Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037

Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY

TR01 Transfer of patent right