CN111641226A - 一种计及自动需求响应的建筑型微电网光伏利用率提高方法 - Google Patents
一种计及自动需求响应的建筑型微电网光伏利用率提高方法 Download PDFInfo
- Publication number
- CN111641226A CN111641226A CN202010425075.XA CN202010425075A CN111641226A CN 111641226 A CN111641226 A CN 111641226A CN 202010425075 A CN202010425075 A CN 202010425075A CN 111641226 A CN111641226 A CN 111641226A
- Authority
- CN
- China
- Prior art keywords
- game
- adr
- strategy
- power
- demand
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000004044 response Effects 0.000 title claims abstract description 53
- 230000006872 improvement Effects 0.000 title claims description 10
- 238000005457 optimization Methods 0.000 claims abstract description 46
- 238000010248 power generation Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 230000007246 mechanism Effects 0.000 claims abstract description 10
- 208000001613 Gambling Diseases 0.000 claims description 18
- 230000005611 electricity Effects 0.000 claims description 16
- 238000004146 energy storage Methods 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 8
- 238000007726 management method Methods 0.000 claims description 8
- 230000008901 benefit Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000010485 coping Effects 0.000 claims description 3
- 230000002542 deteriorative effect Effects 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000013468 resource allocation Methods 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000009472 formulation Methods 0.000 claims 1
- 239000000203 mixture Substances 0.000 claims 1
- 230000002452 interceptive effect Effects 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 47
- 238000009434 installation Methods 0.000 description 15
- 238000004088 simulation Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 5
- 230000007423 decrease Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005622 photoelectricity Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000012085 test solution Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- 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/003—Load forecast, e.g. methods or systems for forecasting future load 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
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
- H02J3/322—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means the battery being on-board an electric or hybrid vehicle, e.g. vehicle to grid arrangements [V2G], power aggregation, use of the battery for network load balancing, coordinated or cooperative battery charging
-
- 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/388—Islanding, i.e. disconnection of local power supply from the network
-
- 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/22—The renewable source being solar 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
- 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
-
- 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/40—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation wherein a plurality of decentralised, dispersed or local energy generation technologies are operated simultaneously
-
- 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
- H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
- H02J2310/10—The network having a local or delimited stationary reach
- H02J2310/12—The local stationary network supplying a household or a building
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/10—Photovoltaic [PV]
-
- 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/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种计及自动需求响应的建筑型微电网光伏利用率提高方法,包括以下步骤:对一天的时间进行离散化处理,均分为多个时段;构建建筑微电网结构模型及其调度模型;基于虚拟成本和虚拟实时价格机制,为博弈参与者制定非合作ADR博弈框架;各博弈参与者独立制定策略并进行交互博弈,求解各自的成本最小化优化策略,从而达到所有参与者最优的目标状态;最后利用非合作博弈机制模拟电力交易过程,以博弈方法实现对MIP集中式策略的改进,证明了所提方法的有效性。本发明建立了一个应用于建筑型微电网系统的非合作ADR博弈框架,旨在提高光伏发电的利用率以达到自主供需平衡的状态,并同时降低MG系统的总成本。
Description
技术领域
本发明涉及一种计及自动需求响应的建筑型微电网光伏利用率提高方法。
背景技术
由于当前一次能源日益枯竭、环境问题日益严峻,而能源需求却持续增长,以光伏能源(PV)为代表的可再生能源的战略地位不断提升。但是,可再生能源来源(RES)的随机性以及微电网(MG)需求侧的能源消耗一定程度上成为其利用率的制约因素之一。
需求响应(DR)指用户对于电价变化或不同激励做出的能量使用策略的调整,被认为是最有效的解决方案,可以提高需求响应资源(DRR)的效率以平衡供需,并实现智能电网中的成本降低和负荷特性优化。随着能源互联网和电力市场的快速发展,自动需求响应(ADR)可以动态地调整负荷需求,以响应变化的电价或奖励金,这代表了需求响应的最新实现。
MG的基本职责是以优化调度方式维持供需平衡,而ADR的主要目标是减少高峰负荷和平衡供需两侧。为了实现带有RES的家庭或商业建筑MG系统的供需之间的自我平衡,一种有效的ADR方法是非常必要的,包括集中式和分散式两种方法。在集中式方法中,需要收集来自用户的所有信息,包括私人数据,以获得针对大规模用户的最佳调度计划。由于隐私暴露,交换数据流量以及解析难度,这将出现难以实际应用的问题。而在分散式方法中,解决方案之一是采用多智能体系统(MAS)在MG中建立能源市场,提供激励机制,以促进DRR和分布式发电机的管理。以分散方式研究ADR方法的一种分析工具是博弈论,在智能电网视觉中,它已应用于DRR的调度,能源管理和MG之间的能源交易,可对上述集中式方法进行改进。
发明内容
为了克服现有技术的不足,本发明面向包括建筑物负载、储能蓄电池(BatteryEnergy Storage,BES)和多个插入式电动汽车(Plug-in Electric Vehicle,PEV)的商业建筑构成的微电网系统,提出了一种计及自动需求响应的建筑型微电网光伏利用率提高方法。通过构建建筑型微电网结构模型并制定博弈参与者的非合作ADR博弈框架,以博弈方法实现对MIP集中式策略的改进,特别考虑高渗透率PV的情况,着重于定量分析DRR对储能型蓄电池系统(Battery Energy Storage System,BESS)的影响分配,旨在提高光伏发电的利用率以达到自主供需平衡的状态,并同时降低MG系统的总成本。
为了实现上述目的,本发明的技术方案包括以下步骤:
一种计及自动需求响应的建筑型微电网光伏利用率提高方法,该方法包括以下步骤:
S1:将一天离散的分为K个等长的周期,时间间隔为Δt,并用k表示任意时段,κ={1,2,…,K}表示操作周期的集合,将N+表示为PEV和BESS的集合,N+={1,2,…,N};
S2:构建建筑型微电网结构模型,包括供应侧,能量服务平台ESP和需求侧,在此基础上建立一个基于混合整数规划MIP的集中式调度模型,并确定该调度模型中的约束条件;
S3:基于虚拟成本和虚拟实时价格vRTP机制,为博弈参与者制定非合作ADR博弈框架,以实现优化目标;
S4:在博弈设定的策略空间中,各博弈参与者根据初始状态和消费需求,考虑设备运行约束和系统约束,独立进行优化决策求解各自的成本最小化优化策略;
S5:证明博弈模型的纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性,达到所有参与者全局最优的目标状态,并证明该纳什均衡解与基于MIP的操作模型的最优解相吻合,以博弈方法实现对MIP集中式策略的改进;
S6:由系统判断是否达到纳什均衡,若是,则输出最终优化策略集合作为所有参与者的优化结果,求解完成;若否,返回步骤S4根据更新后的状态信息重新进行优化。
进一步,在所述步骤S2中,建筑微电网结构模型及其调度模型的构建包括以下过程:
S2-1:将建筑微电网分为:供应侧、ESP和需求侧,供应侧:分布式光伏机组和BESS,并与电网进行能量交换;ESP:对负荷需求进行预测、汇总、更新和广播;需求侧:需求侧负荷可分为两组,即可转移负荷和不可转移负荷;
S2-2:构建建筑微电网结构模型:
1)构建供应侧PV模型
光伏发电是利用半导体“光生伏打效应”将太阳辐射能直接转换为电能的发电方式,单个光伏电池的输出功率很小,因此光伏发电单元将大量的光伏电池串并联,形成光伏阵列;光伏系统通常使用最大功率点跟踪MPPT来使光伏模块在变化的环境中以最大功率点工作,采用现有的光伏电池发电器模型,参数会随着环境的变化而变化,需根据实际的光照强度和温度确定,则k时段的光伏阵列输出功率表示为:
式中,Uk和Ik分别为单个光伏组件在k时段的电压和电流;为Ca光伏组件的串联个数;Cb为光伏组件串的并联个数;τPV为损失因子;
2)构建需求侧PEV模型
用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,设该微电网接入的车辆集合为N,则车辆规模为n=N,对于任意车辆l∈N,其相关参数为:
式中,分别表示车辆l接入微电网的时间和预期离开时间;分别表示车辆动力电池的起始SOC和离开微电网时的期望SOC,SOC表示电池剩余能量与电池容量的比值,因此有 表示PEV电池容量;分别表示额定充、放电功率;
S2-3:在此基础上建立一个基于混合整数规划MIP的集中式调度模型:
S2-4:在建立一个目标函数后,要对其进行约束以符合实际情况,在这个目标函数下,需要考虑以下约束:
1)系统功率平衡约束:
2)能量交换约束:
3)PEV约束:
4)BESS约束:
其中,表示在时刻k从电网中获得的能量,它受能量交换平衡的约束,即从电网获电时,有上限当向电网反向输电时,有上限 和 都为二进制数。Pl,k为第l个PEV的功率,和为PEV充电功率的上下限,和为PEV放电功率的上下限为包含的周期集合,长度为Nl。Sl,k表示第l个PEV的工作状态,有上限和下限η(Pl,k)表示第l个PEV的功率交换效率,表示第l个PEV的电池容量;和分别为k时刻BESS的充放电功率,ηBc和ηBd为其充放电效率;此外还考虑了电池的损耗,yB表示电池寿命,LDOD表示总次数,Lann为一年的循环次数。
在所述步骤S3中,非合作ADR博弈框架的制定包括以下内容:
S3-1:由于ADR的特殊性,ADR业务需要不同领域的主体协同配合,通过ADR系统实现需求响应业务自动化的功能,该系统包括自动需求响应服务器,需求响应运行管理系统和DR客户端;
ADR服务器基于系统内的常规负荷、可再生能源出力水平及PEV集群、储能功率信息制定虚拟实时电价,发布DR事件通知;
DR运行管理监测系统内的可再生能源出力和负荷水平,形成DR需求,并把需求发给ADR服务器;监视系统运行、监测DR实施效果;
ADR客户端向下连接DR用户的各类负荷,收集各PEV和储能的充放电迫切水平参数,并将ADR服务器制定的虚拟价格、FG信号发布给系统内的所有响应主体;
响应主体包括PEV集群和储能系统,ADR项目的针对对象,将响应主体集合设为N+;
S3-2:ADR框架的输入、输出数据:
输入数据:供电侧风、光电等可再生能源出力和储能系统的充放电功率信息采集;需求侧常规负荷、PEV集群用电信息采集;响应主体的迫切性水平参数。
输出数据:vRTP价格信号和FG补偿状态信号;
S3-3:虚拟成本表示的是日常负荷需求与光伏生产的匹配程度:
γC(Lk,k)为k时刻的虚拟成本,σ为正虚拟参数,vk是周期k的独立标准电荷,γCa为总的虚拟收益,为微电网系统中总的负载,虚拟收益与能量需求和转移负载的能力相关,第i个PEV或BESS对于总的收益的贡献为εi,其中和为消耗需求和BESS的容量,α是正相关系数;
S3-4:本发明采用拥塞价格算法定义vRTP机制;
虚拟电价vRTPk取决于微电网中的补偿需求量:
其中,vRTPk是在k时段的虚拟电价;Pi,k是在k时段i用户的供应量;a和γ均为常数;
S3-5:ADR架构旨在实现当可再生能源出力高于负荷水平时,作为响应主体的电动汽车与储能能够根据虚拟价格信号进入充电模式,及时消纳过多的可再生能源;而当可再生能源出力较低时能够根据虚拟信号的引导进行一定程度的放电,从而提高供需平衡度,同时降低系统运行成本;
ADR流程如下:
1)收集供电侧出力信息和需求侧用电信息,计算补偿需求量;
2)如果存在供需不平衡现象,则触发DR事件。收集响应主体的迫切水平参数,根据补偿供应量和补偿需求量制定虚拟实时价格信号以及FG信号,并通知微网系统中的所有响应主体;
Pi,k+1=Pi,k+φi(ULi,k-Pi,kvRTPk) (21)
式中,φi是一个影响算法收敛速度的参数;ULi,k是在k时段响应主体i的充/放电迫切水平参数;Tk表示当前k时段;
S3-6:建立非合作ADR博弈框架:
PEV和BESS的所有用户都互为博弈方,通过最大化博弈方的虚拟收益使光伏利用率达到最大,同时确定优化策略;
ρi为可行的策略集,Pi为第i个参与者的优化策略。
ρi={Pi|约束为(4)-(6),(8)-(13)}(23)
进一步,在所述步骤S4中,非合作ADR博弈过程包括以下步骤:
S4-1:以系统常规负载需求,光伏生产功率和每一个博弈方的消费需求为输入列向量,形成所有参与者的初始状态空间;
S4-2:各博弈方根据初始状态空间,同时考虑设备运行约束、系统约束以及其消费需求选择其调度策略作为初始策略,并广播给其他博弈方;
S4-3:计算并记录各博弈方初始状态下的效用函数,记为目前选用策略的效用函数,用以计算其初始策略的成本,每个博弈方都希望将自己的能源成本降到最低,并且会倾向于选择一个能够最大化自己收益的消费计划,作为应对供需不匹配的电力和其他博弈方选择的策略的最佳策略,即找到一个优化策略Pi *:
S4-4:各博弈方进行操作以期望达到目标状态,计算与记录所得k时段的效用函数,在每个时段k迭代确定其最优消费调度计划Pi,以使其虚拟利润最大化。因此,建立了各博弈方之间的ADR博弈如下:
S4-5:若所得k时段的效用函数小于目前选用策略的效用函数,则更新操作后的策略为目前选用策略,并向其他参与者广播;否则不改变目前选用策略以及不向其他参与者广播;
S4-6:若博弈方进行优化后的选用策略相较于操作前的选用策略效用有了提高,即优化后的选用策略的效用函数小于操作前的选用策略的效用函数,则更新自己的策略;否则,选择目前选用策略作为参与者的最终输出结果,即最优策略;
在每次迭代中,每一个博弈方根据(25)求解优化问题以获得其当前的最优策略,并且通过ESP交流实时信息,重复上述过程,直到博弈收敛到纳什均衡解为止;
在所述步骤S5中,博弈模型纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性证明包括以下步骤:
S5-1:证明博弈模型纳什均衡解的存在唯一性:
收益(效用函数)相对于决策变量Pi,k的二阶导数为
显然为严格凹函数,可得ADR博弈是典型的严格n人博弈,因此证明了纳什均衡解的存在唯一性;
S5-2:证明ADR博弈的唯一纳什均衡解与帕累托最优解的一致性:
和表示的是改变策略后和之前的收益,为基于MIP集中调度的最优解,因此由上式得如果即,总之,任何博弈方都不能从改变策略中获益,因此,得出结论,ADR博弈的唯一纳什均衡解与帕累托最优解是一致的,即通过ADR博弈模型达到所有参与者全局最优的目标状态;
S5-3:证明ADR博弈的最优解与基于MIP的调度优化问题的最优解相吻合:
基于MIP集中调度优化问题为如下:
因为基于MIP集中调度问题是最小化问题,那么显然
由于上述方程适用于任何一个i∈N+,遵循(27)中的条件,得出集中调度问题的最优解构成ADR博弈的纳什均衡解,因此证明了ADR博弈的最优解与基于MIP的调度优化问题的最优解相吻合,体现了以博弈方法实现对MIP集中式策略的改进。
本发明的有益效果是:
1、提出了一种基于MIP的建筑MG运行集中调度框架,在满足用户需求的同时,提高了光伏发电的利用率。
2、引入了MG系统供需不匹配的虚拟成本机制,指导PEVs和BESS的需求响应策略。
3、制定了ADR博弈模型来模拟虚拟交易过程,其中每个用户可以通过分散式的方法获得其最优调度,并给出了求解平衡解的方法和步骤,以博弈方法实现对MIP集中式策略的改进。
附图说明
图1是计及自动需求响应的建筑型微电网光伏利用率提高方法流程图。
图2是微电网最佳调度的示意图。
图3是模拟场景中一年的仿真图。
图4是3季3周的模拟结果。
图5是所提方案与现有方案的比较结果,(a):MILP方案的动态价格信号的选定曲线;(b):四种方案下的系统净负荷曲线。
图6是总成本、光伏生产和BESS安装容量之间的关系。
图7是总成本、光伏生产和BESS容量与PEV数量之间的关系。
图8是不同数量PEV下总成本和光伏生产的随机情景分析结果。
具体实施方式
下面结合附图对本发明做进一步说明。
参照图1~图8,一种计及自动需求响应的建筑型微电网光伏利用率提高方法,该方法包括以下步骤:
S1:将一天离散的分为K个等长的周期,时间间隔为Δt,并用k表示任意时段,κ={1,2,…,K}表示操作周期的集合,将N+表示为PEV和BESS的集合,N+={1,2,…,N};
S2:构建建筑型微电网结构模型,包括供应侧,能量服务平台ESP和需求侧,在此基础上建立一个基于混合整数规划MIP的集中式调度模型,并确定该调度模型中的约束条件;
S3:基于虚拟成本和虚拟实时价格vRTP机制,为博弈参与者制定非合作ADR博弈框架,以实现优化目标;
S4:在博弈设定的策略空间中,各博弈参与者根据初始状态和消费需求,考虑设备运行约束和系统约束,独立进行优化决策求解各自的成本最小化优化策略;
S5:证明博弈模型的纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性,达到所有参与者全局最优的目标状态,并证明该纳什均衡解与基于MIP的操作模型的最优解相吻合,以博弈方法实现对MIP集中式策略的改进;
S6:由系统判断是否达到纳什均衡,若是,则输出最终优化策略集合作为所有参与者的优化结果,求解完成;若否,返回步骤S4根据更新后的状态信息重新进行优化。
进一步,在所述步骤S2中,建筑微电网结构模型及其调度模型的构建包括以下过程:
S2-1:将建筑微电网分为:供应侧、ESP和需求侧,供应侧:分布式光伏机组和BESS,并与电网进行能量交换;ESP:对负荷需求进行预测、汇总、更新和广播;需求侧:需求侧负荷可分为两组,即可转移负荷和不可转移负荷;
S2-2:构建建筑微电网结构模型:
1)构建供应侧PV模型
光伏发电是利用半导体“光生伏打效应”将太阳辐射能直接转换为电能的发电方式,单个光伏电池的输出功率很小,因此光伏发电单元将大量的光伏电池串并联,形成光伏阵列;光伏系统通常使用最大功率点跟踪(MPPT)来使光伏模块在变化的环境中以最大功率点工作,采用现有的光伏电池发电器模型,参数会随着环境的变化而变化,需根据实际的光照强度和温度确定。则k时段的光伏阵列输出功率表示为:
式中,Uk和Ik分别为单个光伏组件在k时段的电压和电流;为Ca光伏组件的串联个数;Cb为光伏组件串的并联个数;τPV为损失因子;
2)构建需求侧PEV模型
用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,设该微电网接入的车辆集合为N,则车辆规模为n=N,对于任意车辆l∈N,其相关参数为:
式中,分别表示车辆l接入微电网的时间和预期离开时间;分别表示车辆动力电池的起始SOC和离开微电网时的期望SOC,SOC表示电池剩余能量与电池容量的比值,因此有 表示PEV电池容量;分别表示额定充、放电功率;
S2-3:在此基础上建立一个基于混合整数规划MIP的集中式调度模型:
S2-4:在建立一个目标函数后,要对其进行约束以符合实际情况,在这个目标函数下,需要考虑以下约束:
1)系统功率平衡约束:
2)能量交换约束:
3)PEV约束:
4)BESS约束:
其中,表示在时刻k从电网中获得的能量,它受能量交换平衡的约束,即从电网获电时,有上限当向电网反向输电时,有上限 和 都为二进制数,Pl,k为第l个PEV的功率,和为PEV充电功率的上下限,和为PEV放电功率的上下限,为包含的周期集合,长度为Nl。Sl,k表示第l个PEV的工作状态,有上限和下限η(Pl,k)表示第l个PEV的功率交换效率,表示第l个PEV的电池容量,和分别为k时刻BESS的充放电功率,ηBc和ηBd为其充放电效率;此外还考虑了电池的损耗,yB表示电池寿命,LDOD表示总次数,Lann为一年的循环次数。
在所述步骤S3中,非合作ADR博弈框架的制定包括以下过程:
S3-1:由于ADR的特殊性,ADR业务需要不同领域的主体协同配合,通过ADR系统实现需求响应业务自动化的功能,该系统包括自动需求响应服务器,需求响应运行管理系统和DR客户端等;
ADR相关的系统和设备描述如下:
ADR服务器基于系统内的常规负荷、可再生能源出力水平及PEV集群、储能功率信息制定虚拟实时电价,发布DR事件通知;
DR运行管理监测系统内的可再生能源出力和负荷水平,形成DR需求,并把需求发给ADR服务器;监视系统运行、监测DR实施效果;
ADR客户端向下连接DR用户的各类负荷,收集各PEV和储能的充放电迫切水平参数,并将ADR服务器制定的虚拟价格、FG信号发布给系统内的所有响应主体。
响应主体包括PEV集群和储能系统,ADR项目的针对对象,将响应主体集合设为N+;
S3-2:ADR框架的输入、输出数据:
输入数据:供电侧风、光电等可再生能源出力和储能系统的充放电功率信息采集;需求侧常规负荷、PEV集群用电信息采集;响应主体的迫切性水平参数。
输出数据:vRTP价格信号和FG补偿状态信号;
S3-3:虚拟成本表示的是日常负荷需求与光伏生产的匹配程度:
γC(Lk,k)为k时刻的虚拟成本,σ为正虚拟参数,vk是周期k的独立标准电荷。γCa为总的虚拟收益。为微电网系统中总的负载,虚拟收益与能量需求和转移负载的能力相关,第i个PEV或BESS对于总的收益的贡献为εi,其中和为消耗需求和BESS的容量,α是正相关系数;
S3-4:采用拥塞价格算法定义vRTP机制;
虚拟电价vRTPk取决于微电网中的补偿需求量:
其中,Pi,k是在k时段i用户的供应量;a和γ均为常数;
S3-5:ADR架构旨在实现当可再生能源出力高于负荷水平时,作为响应主体的电动汽车与储能能够根据虚拟价格信号进入充电模式,及时消纳过多的可再生能源;而当可再生能源出力较低时能够根据虚拟信号的引导进行一定程度的放电,从而提高供需平衡度,同时降低系统运行成本;
ADR流程如下:
1)收集供电侧出力信息和需求侧用电信息,计算补偿需求量;
2)如果存在供需不平衡现象,则触发DR事件。收集响应主体的迫切水平参数,根据补偿供应量和补偿需求量制定虚拟实时价格信号以及FG信号,并通知微网系统中的所有响应主体;
Pi,k+1=Pi,k+φi(ULi,k-Pi,kvRTPk) (21)
式中,φi是一个影响算法收敛速度的参数;ULi,k是在k时段响应主体i的充/放电迫切水平参数;vRTPk是在k时段的虚拟价格信号,Tk表示当前k时段;
S3-6:建立非合作ADR博弈框架:
PEV和BESS的所有用户都互为博弈方,通过最大化博弈方的虚拟收益使光伏利用率达到最大,同时确定优化策略;
ρi为可行的策略集,Pi为第i个参与者的优化策略。
ρi={Pi|约束为(4)-(6),(8)-(13)} (23)
进一步,在所述步骤S4中,非合作ADR博弈过程包括以下步骤:
S4-1:以系统常规负载需求,光伏生产功率和每一个博弈方的消费需求为输入列向量,形成所有参与者的初始状态空间;
S4-2:各博弈方根据初始状态空间,同时考虑设备运行约束、系统约束以及其消费需求选择其调度策略作为初始策略,并广播给其他博弈方;
S4-3:计算并记录各博弈方初始状态下的效用函数,记为目前选用策略的效用函数,用以计算其初始策略的成本,每个博弈方都希望将自己的能源成本降到最低,并且会倾向于选择一个能够最大化自己收益的消费计划,作为应对供需不匹配的电力和其他博弈方选择的策略的最佳策略,即找到一个优化策略Pi *:
S4-4:各博弈方进行操作以期望达到目标状态,计算与记录所得k时段的效用函数,在每个时段k迭代确定其最优消费调度计划Pi,以使其虚拟利润最大化,因此,建立了各博弈方之间的ADR博弈如下:
S4-5:若所得k时段的效用函数小于目前选用策略的效用函数,则更新操作后的策略为目前选用策略,并向其他参与者广播;否则不改变目前选用策略以及不向其他参与者广播;
S4-6:若博弈方进行优化后的选用策略相较于操作前的选用策略效用有了提高,即优化后的选用策略的效用函数小于操作前的选用策略的效用函数,则更新自己的策略;否则,选择目前选用策略作为参与者的最终输出结果,即最优策略;
在每次迭代中,每一个博弈方根据(25)求解优化问题以获得其当前的最优策略,并且通过ESP交流实时信息,重复上述过程,直到博弈收敛到纳什均衡解为止;
在所述步骤S5中,博弈模型纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性证明包括以下步骤:
S5-1:证明博弈模型纳什均衡解的存在唯一性:
收益(效用函数)相对于决策变量Pi,k的二阶导数为
显然为严格凹函数,可得ADR博弈是典型的严格n人博弈,因此证明了纳什均衡解的存在唯一性;
S5-2:证明ADR博弈的唯一纳什均衡解与帕累托最优解的一致性:
和表示的是改变策略后和之前的收益,为基于MIP集中调度的最优解,因此由上式得如果即,总之,任何博弈方都不能从改变策略中获益,因此,得出结论,ADR博弈的唯一纳什均衡解与帕累托最优解是一致的,即通过ADR博弈模型达到所有参与者全局最优的目标状态;
S5-3:证明ADR博弈的最优解与基于MIP的调度优化问题的最优解相吻合:
基于MIP集中调度优化问题为如下:
因为基于MIP集中调度问题是最小化问题,那么显然
由于上述方程适用于任何一个i∈N+,遵循(27)中的条件,得出集中调度问题的最优解构成ADR博弈的纳什均衡解,因此证明了ADR博弈的最优解与基于MIP的调度优化问题的最优解相吻合,体现了以博弈方法实现对MIP集中式策略的改进。
为使本领域技术人员更好地理解本发明,申请人以某商业建筑微电网为实施例验证本项目所提ADR优化方法的有效性。
一、测试系统和参数设置
该微电网的光伏装机容量为1500kW,我们获取了三周的光伏发电数据和不可转移的负荷需求,这些数据是从一整年的夏季,冬季和中期(例如春季或秋季)中选择的。由于这是一种商用MG测试方案,因此PEV的连接期是高峰电价时期,业主不愿为其PEV充电。因此,我们设 相关参数设置如表1所示
表1
在我们的模拟场景中,我们将ADR博弈模型的优化范围设置为一天,这是一年模拟的基础。在完成一年模拟之后,我们评估了性能ADR博弈的计算基于光伏发电利用率(PPUR,φPV)和年度总成本两个指标。
二、仿真结果与分析
我们考虑一个PEV值为n=20且BESS的安装容量为1200kWh的建筑物MG。使用本发明提出的优化模型后,为期3季的3周模拟结果如图4所示。
图4表明,在实施ADR博弈之后,需求侧的负荷曲线与供给侧的PV生产更加一致。提出的ADR博弈可以利用PEV的潜力进行辅助存储,从而可以将负载从PV功率不足的时期转移到PV输出过剩的时期,从而有效地改善MG的负载特性。
三、与现有方案的比较
为了清楚地说明所提出的博弈论ADR方案的有效性,我们重点介绍本发明与分散式智能计费方案和集中式混合整数线性规划(MILP)方案的详细比较。
1)智能充电方案:一种基于优先级的智能充电方案,其中每个PEV和BESS均根据当前系统信息独立确定其策略,以最大程度地利用PV能源。
2)MILP方案:其目标是以集中方式将每日总的用电成本降至最低。注意,MILP方案是由随时间变化的价格信号驱动的,如图5(a)所示,该价格信号可通过AMI获得。
图5(b)给出了四种方案下系统净负荷曲线的比较结果。此外,表2报告了通过提出的ADR方案以及上述两个现有方案实现的系统性能的统计分析。
表2
从图5(b)和表2可以清楚地看出,与智能充电方案和MILP方案相比,拟议的ADR方案在净负载特性,相对峰值降低(RPR)和电池利用率方面均表现出出色的性能。PV能量,有助于提高MG系统的独立运行能力。例如,拟议的ADR计划下的PPUR为95.95%,比智能充电计划和MILP计划下的PPUR分别高7.69%和3.17%。
智能充电方案是为较小规模的MG设计的,而PEV充电模式仅根据当前系统信息确定,因此很难充分考虑全局信息以最大程度地利用PV能源。结果,智能充电方案响应可用功率的能力小于所提出的方案。MILP计划着重于使整个MG的每日总成本最小化,而不是致力于最大程度地利用MG内的PV能量,因此,提出方案产生的PPUR和净负荷特性略好于MG给出的MILP计划。此外,随着该MILP方案中PEV数量的进一步增加,很容易在非高峰价格时段提供新的负荷峰值。由于这些差异,提出的ADR方案在所示的大多数性能中都优于在图5(b)和表2中所述的现有方案。
四、ADR的表现与分析
我们模拟了PEV和BESS典型配置的四种情况,并从各个方面评估了它们的性能。比较结果列于表3。
表3
1)BESS的影响
图6显示了总成本,φPV和BESS安装容量之间的关系。从图6中可以看出,随着BESS的安装容量的增加,总成本首先下降,之后上升。具体来说,当容量达到2200kWh时,MG具有良好的经济性,总成本为$25.4230×104。另一方面,随着BESS容量的增加,MG中的PPUR继续增加,呈现出单调增加的趋势。
将表1中的案例1和案例2进行比较,我们发现,尽管由于配备了BESS,BESS成本增加了$4.7776×104,但案例2的经济性仍然更好,其总成本降低了$1.6395×104。原因是,随着BESS的添加,7PV从74.67%增加到88.06%,导致购电成本降低了$3.7719×104,光伏发电补贴增加了$2.4448×104。
需要强调的是,与单一分布式光伏系统不同,在结合了光伏系统和BESS的建筑物MG中,有效减少了光伏功率的消耗,因此φPV和MG的经济性大大提高。
2)PEV的影响
将表2中的案例2与案例3进行比较,可以看出,当BESS容量相同时,PEV的最佳调度有助于提高PV能源的利用率,从而使φPV从88.06%增加到97.00%。此外,购电成本降低了$1.6148×104,光伏发电补贴提高了$1.6308×104。结果,在情况3中,总成本降低了$1.5582×104。从显然结果可以看出,基于ADR游戏的PEV调度可以有效地提高φPV并实现MG节省成本。
将表2中的案例2与案例4进行比较,可以看出,当PEV数量为20时,BESS的最佳安装容量从2200kWh降低到1200kWh,BESS的年度成本减少了$2.1599×104,因此与案例2相比,案例4的总成本降低了$1.8649×104。可以得出结论,在ADR博弈的推动下,带有V2G选件的PEV有助于降低BESS安装容量的需求,这可以对MG带来显著的经济利益。
图7显示了在不同数量的PEV的情况下,总成本,φPV和BESS安装容量之间的关系。根据图7,可以得出以下结论:
1)随着PEV数量的增加,BESS安装的最佳要求降低,总成本进一步降低。
2)在相同的PEV数下,φPV随BESS的安装容量逐渐增加,直到完全利用PV功率为止。同样,在完全利用PV功率之前,很明显,在相同的BESS容量下,具有更多PEV的MG系统具有更大的φPV。
3)尽管φPV随着BESS安装容量的增加而增加,但MG的总成本与BESS的安装容量没有正相关,但存在最佳安装容量。因此,我们可能需要根据特定要求选择合理的BESS安装容量,以实现这两个指标之间的折衷。
总之,随着V2G技术的发展,PEV能够作为移动存储单元提供有价值的辅助服务,可以代替部分固定存储。显然,PEV有助于提高光伏发电的利用率并有效降低MG的总成本。
五、敏感性分析
实际上,由于不可预测的因素,PEV的初始插件时间,计划的离开时间和初始SOC表现出很大的不确定性。因此,φPV和MG的总成本随上述不确定性而变化。为了探索不确定性对所提出的ADR博弈模型的影响,我们使用随机情景模拟方法研究了PEV数量变化较大时ADR博弈的性能。我们将BESS的安装容量设置为1200kWh。此外,为了在基于场景的随机方法的仿真精度和计算时间之间进行权衡,还生成了600个场景。
图8显示了在不同数量的PEV下总成本和φPV的随机情景分析结果,该结果表明,随着PEV数量的增加,MG的总成本不断下降并最终稳定下来。另一方面,当PEV数为60时,φPV逐渐增加并达到100%左右。此外,总成本和φPV的不确定性随PEV数而增加。但是,当PV利用率接近饱和时,两者的不确定性会降低,并且当PV电量完全耗尽时,它们的不确定性趋于零。
综上可知,本发明基于非合作博弈的自动需求响应方法促进PV能源的利用,从算例分析中,可以得出本发明所提方法能有效地促进光伏能源的利用,同时降低MG系统的总成本。并且证明了在ADR博弈下PEV的调度有助于降低BESS的安装能力,从而可以进一步为MG节省成本。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (5)
1.一种计及自动需求响应的建筑型微电网光伏利用率提高方法,其特征在于,所述方法包括以下步骤:
S1:将一天离散的分为K个等长的周期,时间间隔为Δt,并用k表示任意时段,κ={1,2,…,K}表示操作周期的集合,将N+表示为PEV和BESS的集合,N+={1,2,…,N};
S2:构建建筑型微电网结构模型,包括供应侧,能量服务平台ESP和需求侧,在此基础上建立一个基于混合整数规划MIP的集中式调度模型,并确定该调度模型中的约束条件;
S3:基于虚拟成本和虚拟实时价格vRTP机制,为博弈参与者制定非合作ADR博弈框架,以实现优化目标;
S4:在博弈设定的策略空间中,各博弈参与者根据初始状态和消费需求,考虑设备运行约束和系统约束,独立进行优化决策求解各自的成本最小化优化策略;
S5:证明博弈模型的纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性,达到所有参与者全局最优的目标状态,并证明该纳什均衡解与基于MIP的操作模型的最优解相吻合,以博弈方法实现对MIP集中式策略的改进;
S6:由系统判断是否达到纳什均衡,若是,则输出最终优化策略集合作为所有参与者的优化结果,求解完成;若否,返回步骤S4根据更新后的状态信息重新进行优化。
2.如权利要求1所述的一种计及自动需求响应的建筑型微电网光伏利用率提高方法,其特征在于,所述步骤S2中,建筑微电网结构模型及其调度模型的构建包括以下过程:
S2-1:将建筑微电网分为:供应侧、ESP和需求侧,供应侧:分布式光伏机组和BESS,并与电网进行能量交换;ESP:对负荷需求进行预测、汇总、更新和广播;需求侧:需求侧负荷可分为两组,即可转移负荷和不可转移负荷;
S2-2:构建建筑微电网结构模型:
1)构建供应侧PV模型
光伏发电是利用半导体“光生伏打效应”将太阳辐射能直接转换为电能的发电方式,单个光伏电池的输出功率很小,因此光伏发电单元将大量的光伏电池串并联,形成光伏阵列;光伏系统通常使用最大功率点跟踪MPPT来使光伏模块在变化的环境中以最大功率点工作,采用现有的光伏电池发电器模型,参数会随着环境的变化而变化,需根据实际的光照强度和温度确定,则k时段的光伏阵列输出功率表示为:
式中,Uk和Ik分别为单个光伏组件在k时段的电压和电流;为Ca光伏组件的串联个数;Cb为光伏组件串的并联个数;τPV为损失因子;
2)构建需求侧PEV模型
用一个七维行向量记录电动汽车的电池信息和客户充电需求信息,设该微电网接入的车辆集合为N,则车辆规模为n=|N|,对于任意车辆l∈N,其相关参数为:
式中,Tl in、Tl out分别表示车辆l接入微电网的时间和预期离开时间;分别表示车辆动力电池的起始SOC和离开微电网时的期望SOC,SOC表示电池剩余能量与电池容量的比值,因此有 表示PEV电池容量;Pl maxc、Pl maxd分别表示额定充、放电功率;
S2-3:在此基础上建立一个基于混合整数规划MIP的集中式调度模型:
S2-4:在建立一个目标函数后,要对其进行约束以符合实际情况,在这个目标函数下,需要考虑以下约束:
1)系统功率平衡约束:
2)能量交换约束:
3)PEV约束:
4)BESS约束:
3.如权利要求1或2所述的一种计及自动需求响应的建筑型微电网光伏利用率提高方法,其特征在于,所述步骤S3中,非合作ADR博弈框架的制定包括以下过程:
S3-1:由于ADR的特殊性,ADR业务需要不同领域的主体协同配合,通过ADR系统实现需求响应业务自动化的功能,该系统包括自动需求响应服务器,需求响应运行管理系统和DR客户端;
ADR服务器基于系统内的常规负荷、可再生能源出力水平及PEV集群、储能功率信息制定虚拟实时电价,发布DR事件通知;
DR运行管理监测系统内的可再生能源出力和负荷水平,形成DR需求,并把需求发给ADR服务器;监视系统运行、监测DR实施效果;
ADR客户端向下连接DR用户的各类负荷,收集各PEV和储能的充放电迫切水平参数,并将ADR服务器制定的虚拟价格、FG信号发布给系统内的所有响应主体;
响应主体包括PEV集群和储能系统,是ADR项目的针对对象,将响应主体集合设为N+;
S3-2:ADR框架的输入、输出数据:
输入数据:供电侧可再生能源出力和储能系统的充放电功率信息采集;需求侧常规负荷、PEV集群用电信息采集;响应主体的迫切性水平参数;
输出数据:vRTP价格信号和FG补偿状态信号;
S3-3:虚拟成本表示的是日常负荷需求与光伏生产的匹配程度:
γC(Lk,k)为k时刻的虚拟成本,σ为正虚拟参数,vk是周期k的独立标准电荷,γCa为总的虚拟收益,为微电网系统中总的负载,虚拟收益与能量需求和转移负载的能力相关,第i个PEV或BESS对于总的收益的贡献为εi,其中和为消耗需求和BESS的容量,α是正相关系数;
S3-4:采用拥塞价格算法定义vRTP机制;
虚拟电价vRTPk取决于微电网中的补偿需求量:
其中,Pi,k是在k时段i用户的供应量;a和γ均为常数;
S3-5:ADR流程如下:
1)收集供电侧出力信息和需求侧用电信息,计算补偿需求量;
2)如果存在供需不平衡现象,则触发DR事件,收集响应主体的迫切水平参数,根据补偿供应量和补偿需求量制定虚拟实时价格信号以及FG信号,并通知微网系统中的所有响应主体;
Pi,k+1=Pi,k+φi(ULi,k-Pi,kvRTPk) (21)
式中,φi是一个影响算法收敛速度的参数;ULi,k是在k时段响应主体i的充/放电迫切水平参数;Tk表示当前k时段;
S3-6:建立非合作ADR博弈框架:
PEV和BESS的所有用户都互为博弈方,通过最大化博弈方的虚拟收益使光伏利用率达到最大,同时确定优化策略;
ρi为可行的策略集,Pi为第i个参与者的优化策略;
ρi={Pi|约束为(4)-(6),(8)-(13)} (23)
4.如权利要求1或2所述的一种计及自动需求响应的建筑型微电网光伏利用率提高方法,其特征在于,所述步骤S4中,非合作ADR博弈过程包括以下步骤:
S4-1:以系统常规负载需求,光伏生产功率和每一个博弈方的消费需求为输入列向量,形成所有参与者的初始状态空间;
S4-2:各博弈方根据初始状态空间,同时考虑设备运行约束、系统约束以及其消费需求选择其调度策略作为初始策略,并广播给其他博弈方;
S4-3:计算并记录各博弈方初始状态下的效用函数,记为目前选用策略的效用函数,用以计算其初始策略的成本,每个博弈方都希望将自己的能源成本降到最低,并且会倾向于选择一个能够最大化自己收益的消费计划,作为应对供需不匹配的电力和其他博弈方选择的策略的最佳策略,即找到一个优化策略Pi *:
S4-4:各博弈方进行操作以期望达到目标状态,计算与记录所得k时段的效用函数,在每个时段k迭代确定其最优消费调度计划Pi,以使其虚拟利润最大化,因此,建立了各博弈方之间的ADR博弈如下:
S4-5:若所得k时段的效用函数小于目前选用策略的效用函数,则更新操作后的策略为目前选用策略,并向其他参与者广播;否则不改变目前选用策略以及不向其他参与者广播;
S4-6:若博弈方进行优化后的选用策略相较于操作前的选用策略效用有了提高,即优化后的选用策略的效用函数小于操作前的选用策略的效用函数,则更新自己的策略;否则,选择目前选用策略作为参与者的最终输出结果,即最优策略;
在每次迭代中,每一个博弈方根据(25)求解优化问题以获得其当前的最优策略,并且通过ESP交流实时信息,重复上述过程,直到博弈收敛到纳什均衡解为止;
5.如权利要求1或2所述的一种计及自动需求响应的建筑型微电网光伏利用率提高方法,其特征在于,所述步骤S5中,博弈模型纳什均衡解存在的唯一性以及纳什均衡与帕累托最优的一致性证明包括以下步骤:
S5-1:证明博弈模型纳什均衡解的存在唯一性:
收益相对于决策变量Pi,k的二阶导数为
显然为严格凹函数,可得ADR博弈是典型的严格n人博弈,因此证明了纳什均衡解的存在唯一性;
S5-2:证明ADR博弈的唯一纳什均衡解与帕累托最优解的一致性:
和表示的是改变策略后和之前的收益,为基于MIP集中调度的最优解,因此由上式可得如果即,总之,任何博弈方都不能从改变策略中获益,因此,我们可以得出结论,ADR博弈的唯一纳什均衡解与帕累托最优解是一致的,即通过ADR博弈模型达到所有参与者全局最优的目标状态;
S5-3:证明ADR博弈的最优解与基于MIP的调度优化问题的最优解相吻合:
基于MIP集中调度优化问题为如下:
因为基于MIP集中调度问题是最小化问题,那么显然
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010425075.XA CN111641226B (zh) | 2020-05-19 | 2020-05-19 | 一种计及自动需求响应的建筑型微电网光伏利用率提高方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010425075.XA CN111641226B (zh) | 2020-05-19 | 2020-05-19 | 一种计及自动需求响应的建筑型微电网光伏利用率提高方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111641226A true CN111641226A (zh) | 2020-09-08 |
CN111641226B CN111641226B (zh) | 2021-12-17 |
Family
ID=72331109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010425075.XA Active CN111641226B (zh) | 2020-05-19 | 2020-05-19 | 一种计及自动需求响应的建筑型微电网光伏利用率提高方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111641226B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112129995A (zh) * | 2020-09-25 | 2020-12-25 | 中国电力科学研究院有限公司 | 一种可再生能源电力生产和消纳量计量方法及系统 |
CN112202201A (zh) * | 2020-09-24 | 2021-01-08 | 上海电机学院 | 一种考虑需求响应和电动汽车的联合微电网运行策略 |
CN112636330A (zh) * | 2020-12-02 | 2021-04-09 | 上海交通大学 | 一种基于用户贡献行为的综合能源系统调峰方法 |
CN113270872A (zh) * | 2020-12-23 | 2021-08-17 | 西北工业大学深圳研究院 | 一种电网负载频率控制优化方法、系统、设备及存储介质 |
CN113283751A (zh) * | 2021-05-26 | 2021-08-20 | 西安图迹信息科技有限公司 | 基于动态博弈的商业建筑需求响应方法和系统 |
CN115438931A (zh) * | 2022-08-22 | 2022-12-06 | 成都飞机工业(集团)有限责任公司 | 一种生产线装配作业调度方法、装置、设备及介质 |
CN116565882A (zh) * | 2023-06-29 | 2023-08-08 | 北京大学 | 分布式需求响应方法、装置、系统及介质 |
CN118157110A (zh) * | 2024-03-04 | 2024-06-07 | 天津大学 | 一种微电网分布式优化方法、装置、存储介质及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104835080A (zh) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | 微电网日内调度计划混合整数规划模型的建模方法 |
CN107612041A (zh) * | 2017-08-23 | 2018-01-19 | 浙江工业大学 | 一种考虑不确定性的基于事件驱动的微电网自动需求响应方法 |
KR101811426B1 (ko) * | 2016-10-07 | 2018-01-25 | 문감사 | 태양광 발전을 이용한 통합형 전력관리 시스템 |
CN109657993A (zh) * | 2018-12-24 | 2019-04-19 | 浙江工业大学 | 一种基于非合作博弈的能源局域网储能系统自动需求响应方法 |
-
2020
- 2020-05-19 CN CN202010425075.XA patent/CN111641226B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104835080A (zh) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | 微电网日内调度计划混合整数规划模型的建模方法 |
KR101811426B1 (ko) * | 2016-10-07 | 2018-01-25 | 문감사 | 태양광 발전을 이용한 통합형 전력관리 시스템 |
CN107612041A (zh) * | 2017-08-23 | 2018-01-19 | 浙江工业大学 | 一种考虑不确定性的基于事件驱动的微电网自动需求响应方法 |
CN109657993A (zh) * | 2018-12-24 | 2019-04-19 | 浙江工业大学 | 一种基于非合作博弈的能源局域网储能系统自动需求响应方法 |
Non-Patent Citations (1)
Title |
---|
杨晓东: "需求响应参与的微电网系统调控模型与方法研究", 《中国博士学位论文全文数据库(电子期刊)》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112202201A (zh) * | 2020-09-24 | 2021-01-08 | 上海电机学院 | 一种考虑需求响应和电动汽车的联合微电网运行策略 |
CN112202201B (zh) * | 2020-09-24 | 2024-05-14 | 上海电机学院 | 一种考虑需求响应和电动汽车的联合微电网运行策略 |
CN112129995A (zh) * | 2020-09-25 | 2020-12-25 | 中国电力科学研究院有限公司 | 一种可再生能源电力生产和消纳量计量方法及系统 |
CN112129995B (zh) * | 2020-09-25 | 2023-03-21 | 中国电力科学研究院有限公司 | 一种可再生能源电力生产和消纳量计量方法及系统 |
CN112636330B (zh) * | 2020-12-02 | 2023-01-10 | 上海交通大学 | 一种基于用户贡献行为的综合能源系统调峰方法 |
CN112636330A (zh) * | 2020-12-02 | 2021-04-09 | 上海交通大学 | 一种基于用户贡献行为的综合能源系统调峰方法 |
CN113270872A (zh) * | 2020-12-23 | 2021-08-17 | 西北工业大学深圳研究院 | 一种电网负载频率控制优化方法、系统、设备及存储介质 |
CN113283751A (zh) * | 2021-05-26 | 2021-08-20 | 西安图迹信息科技有限公司 | 基于动态博弈的商业建筑需求响应方法和系统 |
CN115438931A (zh) * | 2022-08-22 | 2022-12-06 | 成都飞机工业(集团)有限责任公司 | 一种生产线装配作业调度方法、装置、设备及介质 |
CN115438931B (zh) * | 2022-08-22 | 2024-03-15 | 成都飞机工业(集团)有限责任公司 | 一种生产线装配作业调度方法、装置、设备及介质 |
CN116565882A (zh) * | 2023-06-29 | 2023-08-08 | 北京大学 | 分布式需求响应方法、装置、系统及介质 |
CN116565882B (zh) * | 2023-06-29 | 2023-09-19 | 北京大学 | 分布式需求响应方法、装置、系统及介质 |
CN118157110A (zh) * | 2024-03-04 | 2024-06-07 | 天津大学 | 一种微电网分布式优化方法、装置、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111641226B (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111641226B (zh) | 一种计及自动需求响应的建筑型微电网光伏利用率提高方法 | |
Wang et al. | Evaluation of economic benefits of virtual power plant between demand and plant sides based on cooperative game theory | |
CN109657993B (zh) | 一种基于非合作博弈的能源局域网储能系统自动需求响应方法 | |
Zheng et al. | Economic and environmental benefits of coordinating dispatch among distributed electricity storage | |
CN107423852A (zh) | 一种计及典型场景的光储联合电站优化运营方法 | |
Lim et al. | Optimal allocation of energy storage and solar photovoltaic systems with residential demand scheduling | |
Yu | System contributions of residential battery systems: New perspectives on PV self-consumption | |
CN107612041B (zh) | 一种考虑不确定性的基于事件驱动的微电网自动需求响应方法 | |
CN111682536B (zh) | 虚拟电厂参与日前双重市场的随机-鲁棒优化运行方法 | |
CN111079971A (zh) | 一种考虑车、站、网三方的充电站定价方法 | |
CN112365021A (zh) | 一种基于混合储能的区域综合能源系统规划优化方法 | |
Angenendt et al. | Evaluation of the effects of frequency restoration reserves market participation with photovoltaic battery energy storage systems and power-to-heat coupling | |
Li et al. | Two-stage community energy trading under end-edge-cloud orchestration | |
Hayati et al. | A two-stage stochastic optimization scheduling approach for integrating renewable energy sources and deferrable demand in the spinning reserve market | |
CN117639113A (zh) | 智慧微电网智能电力分配方法、装置及存储介质 | |
CN114493743B (zh) | 基于agc调频分区控制光储联合系统参与市场投标方法 | |
Fu et al. | Research on the stackelberg game method of building micro-grid with electric vehicles | |
Yang et al. | Enhancing utilization of PV energy in building microgrids via autonomous demand response | |
CN113746105B (zh) | 电力需求响应的优化控制方法、装置、设备和存储介质 | |
CN117709552A (zh) | 基于分布式能源分类聚合的虚拟电厂调度优化方法 | |
Mulleriyawage et al. | Battery system selection in DC microgrids for residential applications: An Australian case study | |
Mammoli et al. | Distributed control strategies for high-penetration commercial-building-scale thermal storage | |
CN106227986A (zh) | 一种分布式电源与智能停车场的联合部署方法及装置 | |
Hussain et al. | Energy Allocation of the Community Energy Storage System: A Contribution-Based Incentive Mechanism | |
Zhu et al. | Optimal scheduling model for virtual power plant participating in energy and regulation markets |
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 |