CN112001752A - 基于有限理性的多虚拟电厂动态博弈交易行为分析方法 - Google Patents
基于有限理性的多虚拟电厂动态博弈交易行为分析方法 Download PDFInfo
- Publication number
- CN112001752A CN112001752A CN202010832963.3A CN202010832963A CN112001752A CN 112001752 A CN112001752 A CN 112001752A CN 202010832963 A CN202010832963 A CN 202010832963A CN 112001752 A CN112001752 A CN 112001752A
- Authority
- CN
- China
- Prior art keywords
- vpp
- power
- dso
- strategy
- bidding
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- 230000006399 behavior Effects 0.000 claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000011161 development Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims abstract description 15
- 238000011160 research Methods 0.000 claims abstract description 12
- 230000005611 electricity Effects 0.000 claims description 85
- 238000012725 vapour phase polymerization Methods 0.000 claims description 62
- 230000006870 function Effects 0.000 claims description 31
- 230000008901 benefit Effects 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 20
- 238000010248 power generation Methods 0.000 claims description 15
- 230000002860 competitive effect Effects 0.000 claims description 11
- 238000004146 energy storage Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 238000007599 discharging Methods 0.000 claims description 9
- 230000009194 climbing Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 101000805729 Homo sapiens V-type proton ATPase 116 kDa subunit a 1 Proteins 0.000 description 3
- 102100037979 V-type proton ATPase 116 kDa subunit a 1 Human genes 0.000 description 3
- 101000854873 Homo sapiens V-type proton ATPase 116 kDa subunit a 4 Proteins 0.000 description 2
- 101000806601 Homo sapiens V-type proton ATPase catalytic subunit A Proteins 0.000 description 2
- 101000850434 Homo sapiens V-type proton ATPase subunit B, brain isoform Proteins 0.000 description 2
- 101000670953 Homo sapiens V-type proton ATPase subunit B, kidney isoform Proteins 0.000 description 2
- 102100037466 V-type proton ATPase catalytic subunit A Human genes 0.000 description 2
- 102100033476 V-type proton ATPase subunit B, brain isoform Human genes 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000005612 types of electricity Effects 0.000 description 2
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0206—Price or cost determination based on market factors
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Accounting & Taxation (AREA)
- Development Economics (AREA)
- Finance (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Tourism & Hospitality (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其包括:首先,各虚拟电厂竞标主体充分考虑自身发展阶段的目标需求和研究上层运营商动态定价行为;其次,采用不同交易目标建模对多虚拟电厂进行交易行为分析,采用粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而改进自身目标使其逐步变优;最后,研究了多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解;所提动态博弈模型的求解具有良好的收敛性,为虚拟电厂制定不同交易目标参与市场交易提供新的思路和参考。
Description
技术领域
本发明属于多虚拟电厂交易行为分析技术领域,具体涉及基于有限理性的多虚拟电厂动态博弈交易行为分析方法。
背景技术
在可再生能源和电能替代的电力发展战略下,分布式电源、储能、可响应负荷在配电网层面大力推广应用,使得传统的规模化电力生产和远距离传输逐渐向分布式能源发展。但其大规模入网对电网的安全运行造成巨大的挑战。目前,国内外学者从单利益主体出发,针对虚拟电厂的发电调度、竞价模式等方面已经做了不少研究。但是随着近年来国家电网发展战略的改变,使得更多元化的社会资本涌入电力市场,各VPP将可能属于不同利益主体运行。在各利益主体自主理性的驱使下,各主体追求自身利益最大化,传统的单一主体的优化调度方法将很难适用。然而已有研究多数集中在给定交易电价的前提下,研究VPP之间的博弈关系,但仍存在如下问题:一是未考虑虚拟电厂对市场交易电价的影响,仅作为价格接受者参与市场运行,本质上仍为VPP内部优化;二是决策环境过于理想化,多是完全理性条件下的博弈研究,不仅要求主体自身有完美的理性,还要求各主体互相信任对方的理性,这在实际情况中很难达到。同时,已有文献研究的多是传统的以收益最大化或是运行成本最低为目标制定最优竞标方案,但在未来多元化市场中,各利益主体在不同的发展阶段有着不同的发展需求,仅以收益最大化或运行成本最低制定目标的交易行为不能够在完全竞争市场上的占据竞标优势。
发明内容
本发明目的在于提供基于有限理性的多虚拟电厂动态博弈交易行为分析方法,用于解决上述现有技术中存在的技术问题之一,如:现有技术中,仍存在如下问题:一是未考虑虚拟电厂对市场交易电价的影响,仅作为价格接受者参与市场运行,本质上仍为VPP内部优化;二是决策环境过于理想化,多是完全理性条件下的博弈研究,不仅要求主体自身有完美的理性,还要求各主体互相信任对方的理性,这在实际情况中很难达到。同时,已有文献研究的多是传统的以收益最大化或是运行成本最低为目标制定最优竞标方案,但在未来多元化市场中,各利益主体在不同的发展阶段有着不同的发展需求,仅以收益最大化或运行成本最低制定目标的交易行为不能够在完全竞争市场上的占据竞标优势。
为实现上述目的,本发明的技术方案是:
基于有限理性的多虚拟电厂动态博弈交易行为分析方法,包括以下步骤:
S1:各虚拟电厂竞标主体根据自身发展阶段的目标需求和研究上层运营商动态定价行为;
S2:采用不同交易目标建模对多虚拟电厂进行交易行为分析;
S3:基于粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而使自身目标逐步变优;
S4:研究多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解。
进一步的,根据VPP发展的不同阶段和竞争对手的决策行为对自身优化决策的影响,建立多VPP有不同目标的决策方案;考虑包括VPP内的购售电电量和分布式电源出力为策略,设计多种不同目标的VPP竞标决策方案,一是考虑到不同区域VPP的用户用电差异和VPP联盟目标追求不同,VPP可选择不同的交易目标和决策方案;二是多VPP在进行博弈过程中,各主体会综合其余竞争者的市场交易行为后可改变自身的行为规则,进行目标的转换;在进行博弈的过程中,通过进化学习会不断获取竞争者的信息;随着自身获取信息的不断增多所作出的策略也会逐步变优;最终直至任何主体都不能从单方面从改变竞标策略中获益,既达到稳定策略。
进一步的,多利益主体在有限理性条件下的稳定均衡策略可定义为:在N个利益主体进行动态博弈中,存在稳定策略组合X=(X1,X2,...,XN),对于任何其他策略组合Y=(Y1,Y2,...,YN)|Y≠X,存在有主体i和ε0∈(0,1)使得:
式中:I(Yi,S-i)和I(Xi,S-i)分别为主体i采用Yi和Xi策略时的目标效用情况;S-i为除去主体i以外系统其他主体的策略组合。
进一步的,VPP的博弈策略考虑为VPP内各个分布式电源的出力计划、储能的充放电功率、可响应负荷的响应功率、风光等清洁能源的输出功率;VPP根据目前已经掌握的有限的市场信息和其余VPP竞标历史数据得到系统负荷需求量以及竞争者可能的竞标策略,来实现最大化自身效用。
进一步的,步骤S2中不同交易目标包括:
保守型VPP的博弈目标函数考虑为最大化自身效益,充分整合内部资源,包括分布式电源发电成本、储能成本、可响应负荷补偿成本;决策出最优竞标策略;VPP竞标目标函数为:
式中:λD-V,sell、λD-V,buy分别为VPP向DSO的售电电价和购电电价;分别为VPP向DSO购买的电量和出售的电量;N为VPP总数集合,VPPi∈N;Ni为第i个VPP中含有的DER集合;为DER中MT的发电成本,aj、bj、cj为DERj中MT机组的成本系数,为MT的输出功率;为ESS运行成本,为ESS的调度成本系数,为ESS的充放电功率,其值大于0表示放电,小于0表示充电;为可响应负荷成本,为响应负荷补偿电价,为响应负荷的响应功率;θ为政府对新能源发电的补贴电价;分别为光伏、风电出力功率;
激进型VPP的博弈目标函数考虑为市场占有率规模最大化,在整合内部资源满足各单元约束出力的前提下,最大化在市场中进行电量竞标,以尽可能多的在DSO层面得到交换电量,获得最大功率购买量;VPP竞标目标函数为:
竞标约束:
在获得DSO统一购售电价和交易电量后,在满足激进型VPP目标函数前提下,对VPP内部其他可调度资源进行协调,满足剩余的负荷电量缺额,电量缺额调度原则满足成本最低原则;
稳定型VPP的博弈目标函数考虑为可再生能源的最大化利用,可再生能源考虑为风电、光伏,同时充分考虑可再生能源发电的不确定性,利用ESS平抑可再生能源出力的波动性对系统供电稳定性的影响;其次再考虑VPP内部其他各单元出力情况和向上级DSO进行电量竞标,VPP竞标目标函数为:
进一步的,所述VPP在进行决策响应时,无论执行哪一种交易目标,都需要满足的约束条件包括:
功率平衡约束:
燃气轮机组约束:
可响应负荷约束:
可再生能源出力约束:
储能约束:
式中:为ESS放电和充电状态的0-1变量;NESS为ESS放电和充电状态转换次数上限;为放电功率上下限;为充电功率上下限;为ESS在第t时段的容量状态;为考虑ESS寿命等因素的容量上下限值;分别为充电和放电效率系数,其取值范围为
进一步的,所述交易在DSO层面考虑为:
DSO为VPP制定售电电价λD-V,sell和购电电价λD-V,buy;DSO以最大化自身效益作为博弈目标,包括与上级电网和下层多VPP之间的购售电成本和收益;
maxuDSO=FVPP-FDSO
式中:FVPP为DSO与VPP之间进行交易获得的收益;FDSO为DSO与上级电网进行交易所花费的成本;各项的具体表达式如下:
式中:T为总调度时段;NVPP为VPP数量,为由DSO制定的实时电价,为第i个VPP在第t调度时段与DSO交互的功率值;为DSO与上级电网之间的边界节点电价,为DSO与上级电网之间的传输功率,其值大于0表示从上级电网购电,小于0表示向上级电网售电。
进一步的,获得竞争者竞标数据后该竞争者策略判断指标系数如下:
式中:ξ1,ξ2,ξ3为对竞争者策略判断的指标系数,当ξ1≥ξlim1时判断为“激进型”,当ξ2≤ξlim2时判断为“稳定型”,当ξ3≤ξlim3时判断为“保守型”,其中ξlim1,ξlim2,ξlim3为指标系数判断常量,通过各VPP竞标历史数据获得;若在一轮博弈结束策略指标判断不在指标系数判断常量限值内,则不会输出策略结果,继续开始下一轮迭代;
各VPP通过进化学习对博弈策略的选取,决定该VPP的总效益;每种策略对应于该策略下的交易行为方式,用ρmn表示VPPi交易行为方式x(u)从策略m到策略n的条件切换状态;在获取竞争对手竞标信息后,VPPi学习进化自身策略;如果策略n比策略m具有更高的效用,那么策略m将学习策略n并作出改变;用方程表示为:
式中:xm(u)和xn(u)分别为VPPi在策略m和n下的交易行为方式,当大于0时,则根据条件切换状态更新VPP策略,当始终小于0时,则达到稳定均衡策略;um和un分别表示策略m和n对应的交易行为方式效用,[un-um]+为max{0,un-um}。
与现有技术相比,本发明所具有的有益效果为:
本方案的一个创新点在于,首先,各虚拟电厂竞标主体充分考虑自身发展阶段的目标需求和研究上层运营商动态定价行为;其次,采用不同交易目标建模对多虚拟电厂进行交易行为分析,采用粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而改进自身目标使其逐步变优;最后,研究了多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解;所提动态博弈模型的求解具有良好的收敛性,为虚拟电厂制定不同交易目标参与市场交易提供新的思路和参考。
附图说明
图1为基于有限理性的多虚拟电厂动态博弈交易行为分析方法的流程图。
图2为多虚拟电厂的系统结构图。
图3为多虚拟电厂不同交易目标示意图。
图4为多虚拟电厂动态博弈求解流程。
图5为算法收敛过程示意图。
图6为交易电价优化示意图。
图7为虚拟电厂1交易结果示意图。
图8为虚拟电厂2交易结果示意图。
图9为虚拟电厂3交易结果示意图。
具体实施方式
下面结合本发明的附图1-9,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例:
一种基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,包括以下步骤:
各虚拟电厂竞标主体充分考虑自身发展阶段的目标需求和研究上层运营商动态定价行为;
采用不同交易目标建模对多虚拟电厂进行交易行为分析;
基于粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而改进自身目标使其逐步变优;
研究多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解。
考虑VPP发展的不同阶段和竞争对手的决策行为对自身优化决策的影响,因此建立多VPP有不同目标的决策方案。本文重点考虑包括VPP内的购售电电量和分布式电源等出力为策略,设计了多种不同目标的VPP竞标决策方案,一是考虑到不同区域VPP的用户用电差异和VPP联盟目标追求不同,VPP可能选择不同的交易目标和决策方案;二是多VPP在进行博弈过程中,各主体会综合其余竞争者的市场交易行为后可能改变自身的行为规则,进行目标的转换。在进行博弈的过程中,通过进化学习会不断获取竞争者的信息。随着自身获取信息的不断增多所作出的策略也会逐步变优。最终直至任何主体都不能从单方面从改变竞标策略中获益,既达到稳定策略。
多利益主体在有限理性条件下的稳定均衡策略可以定义为:在N个利益主体进行动态博弈中,存在稳定策略组合X=(X1,X2,...,XN),对于任何其他策略组合Y=(Y1,Y2,...,YN)|Y≠X,存在有主体i和ε0∈(0,1)使得:
式中:I(Yi,S-i)和I(Xi,S-i)分别为主体i采用Yi和Xi策略时的目标效用情况;S-i为除去主体i以外系统其他主体的策略组合。
VPP的博弈策略考虑为VPP内各个分布式电源的出力计划、储能的充放电功率、可响应负荷的响应功率、风光等清洁能源的输出功率。VPP根据目前已经掌握的有限的市场信息和其余VPP竞标历史数据得到系统负荷需求量以及竞争者可能的竞标策略,来实现最大化自身效用。需要说明的是,这里的“效用”不同于传统优化调度的效用函数,而是指制定符合自身发展的最优策略。
保守型VPP的博弈目标函数考虑为最大化自身效益,充分整合内部资源,包括分布式电源发电成本、储能成本、可响应负荷补偿等成本。决策出最优竞标策略。VPP竞标目标函数为:
式中:λD-V,sell、λD-V,buy分别为VPP向DSO的售电电价和购电电价;分别为VPP向DSO购买的电量和出售的电量;N为VPP总数集合,VPPi∈N;Ni为第i个VPP中含有的DER集合;为DER中MT的发电成本,aj、bj、cj为DERj中MT机组的成本系数,为MT的输出功率;为ESS运行成本,为ESS的调度成本系数,为ESS的充放电功率,其值大于0表示放电,小于0表示充电;为可响应负荷成本,为响应负荷补偿电价,为响应负荷的响应功率;θ为政府对新能源发电的补贴电价;分别为光伏、风电出力功率。本文中主要各方主体的考虑运营效益,忽略了相对较低的新能源发电的运行和维护成本。
激进型VPP的博弈目标函数考虑为市场占有率规模最大化,在整合内部资源满足各单元约束出力的前提下,最大化在市场中进行电量竞标,以尽可能多的在DSO层面得到交换电量,获得最大功率购买量。VPP竞标目标函数为:
竞标约束:
在获得DSO统一购售电价和交易电量后,在满足目标函数式(5)前提下,对VPP内部其他可调度资源进行协调,满足剩余的负荷电量缺额,电量缺额调度原则满足成本最低原则。
稳定型VPP的博弈目标函数考虑为可再生能源的最大化利用,本文可再生能源考虑为风电、光伏,同时充分考虑可再生能源发电的不确定性,利用ESS平抑可再生能源出力的波动性对系统供电稳定性的影响。其次再考虑VPP内部其他各单元出力情况和向上级DSO进行电量竞标,VPP竞标目标函数为:
1)功率平衡约束:
2)燃气轮机组约束:
3)可响应负荷约束:
4)可再生能源出力约束:
5)储能约束:
式中:为ESS放电和充电状态的0-1变量;NESS为ESS放电和充电状态转换次数上限;为放电功率上下限;为充电功率上下限;为ESS在第t时段的容量状态;为考虑ESS寿命等因素的容量上下限值;分别为充电和放电效率系数,其取值范围为
DSO为VPP制定售电电价λD-V,sell和购电电价λD-V,buy。DSO以最大化自身效益作为博弈目标,包括与上级电网和下层多VPP之间的购售电成本和收益。
maxuDSO=FVPP-FDSO
式中:FVPP为DSO与VPP之间进行交易获得的收益;FDSO为DSO与上级电网进行交易所花费的成本。各项的具体表达式如下:
式中:T为总调度时段;NVPP为VPP数量,为由DSO制定的实时电价,为第i个VPP在第t调度时段与DSO交互的功率值;为DSO与上级电网之间的边界节点电价,为DSO与上级电网之间的传输功率,其值大于0表示从上级电网购电,小于0表示向上级电网售电。
获得竞争者竞标数据后该竞争者策略判断指标系数如下:
式中:ξ1,ξ2,ξ3为对竞争者策略判断的指标系数,当ξ1≥ξlim1时判断为“激进型”,当ξ2≤ξlim2时判断为“稳定型”,当ξ3≤ξlim3时判断为“保守型”,其中ξlim1,ξlim2,ξlim3为指标系数判断常量,通过各VPP竞标历史数据获得。若在一轮博弈结束策略指标判断不在指标系数判断常量限值内,则不会输出策略结果,继续开始下一轮迭代。
各VPP通过进化学习对博弈策略的选取,决定了该VPP的总效益。每种策略对应于该策略下的交易行为方式,用ρmn表示VPPi交易行为方式x(u)从策略m到策略n的条件切换状态。在获取竞争对手竞标信息后,VPPi学习进化自身策略。如果策略n比策略m具有更高的效用,那么策略m将学习策略n并作出改变。用方程表示为:
式中:xm(u)和xn(u)分别为VPPi在策略m和n下的交易行为方式,当大于0时,则根据条件切换状态更新VPP策略,当始终小于0时,则达到稳定均衡策略;um和un分别表示策略m和n对应的交易行为方式效用,[un-um]+为max{0,un-um}。
参考图1,图1为基于有限理性的多虚拟电厂动态博弈交易行为分析方法的流程图;该流程图包括步骤101到步骤104;
在步骤101中,各虚拟电厂竞标主体充分考虑自身发展阶段的目标需求和研究上层运营商动态定价行为;
在步骤102中,采用不同交易目标建模对多虚拟电厂进行交易行为分析;
在本发明的一个实施例中,整个配电系统运营的过程中,每个结成联盟的VPP可能会出现不同类型的购售电状态:
电量富余(售电VPP)、电量缺额(购电VPP)、和电量刚好自给自足三种类型。因此,在多VPP系统中存在着购售电需求。
VPP在满足自身内部供电的前提下,综合考虑各机组出力、可再生能源发电、电价信息和其他运营商的策略信息等,对其内部资源进行协调,制定最符合自身发展的目标,在满足各种约束情况下,与其余竞争者进行博弈。
DSO层调度中心综合考虑下层各VPP上传的电量信息,考虑最大化社会效益,制定与多VPP之间的统一市场购售电价和各VPP的交易电量,依据DSO对各VPP制定的电价和交易电量,VPP合理调整内部资源出力,制定与DSO之间的交易电量。
在步骤103中,基于粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而改进自身目标使其逐步变优;
当博弈开始后,参与者可以根据双方每一轮策略进行学习和改进,是一个连续的、重复不断的交易过程。
市场中各VPP作为个体参与博弈,分别考虑其余竞争者的决策信息制定最大化自身效用的目标函数策略,VPPi在获得竞争者的竞标数据后求解自身最优决策,重复i∈N得到各主体第一轮博弈策略。
各主体更新上一轮博弈策略直至各主体都不能从单方面改变竞标策略中获益,达到最后稳定均衡状态。
在步骤104中,研究多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解。
在本发明的一个实施例中,图4示出了多虚拟电厂动态博弈的求解流程。
在模型求解中,DSO给出购售电价。各VPP在博弈过程中,对不断获取的电价信息和其余竞争者制定的策略信息进行学习,来选择最适合自身交易目标优化求解。
在DSO层面,以电价策略作为粒子,以此电价下的收益作为粒子的适应度,使用智能优化算法寻求最优电价策略。本文整个求解过程可以看作DSO电价寻优和多VPP策略学习求解两部分。采用了改进粒子群算法结合Yalmip平台调用CPLEX求解器对多VPP不同交易目标动态博弈进化学习模型进行求解。
参考图2,图2为多虚拟电厂的系统结构图。在VPP内部主要包括微型燃汽轮机组、风电机组、光伏系统、储能系统、可响应负荷等部分组成。其中,微型燃汽轮机组所具有稳定的功率输出、能够快速响应、能够灵活应对负荷变化系统所要求满足的要求,因此将微型燃汽轮机组作为VPP内部的可调度分布式电源。
VPP作为一个整体对外加入DSO以及电力市场的调度,在VPP内部主要负责对多个DER进行协调运行控制,同时对用户负荷、设备维护管理商等多方进行交涉互动。VPP通过能量管理平台对内部各资源进行调控以及对上层DSO进行能量和信息交互。
在能量管理平台具备日前预测、决策优化、实时监控等功能,在其数据库中可存储和处理DER内部负荷历史用电情况和各电源的出力情况、用户信息、天气数据、上层DSO价格信息,这些可使VPP进行更加精准的决策。
参考图3,一是考虑到不同区域VPP的用户用电差异和VPP联盟目标追求不同,VPP可能选择不同的交易目标和决策方案;二是多VPP在进行博弈过程中,各主体会综合其余竞争者的市场交易行为后可能改变自身的行为规则,进行目标的转换。因此为VPP制定不同交易目标类型的目标函数。
在本发明的一个实施例中,各VPP内部组成元素如表1所示。选择对应的成本函数、约束条件和初始目标函数,构成包含3个VPP的有限理性博弈模型。基于上文所述的算例系统,DSO收到各VPP运营商提交的竞标策略后,确定市场交易电价,最大化社会效益。经过19轮博弈之后各VPP经进化学习后选择策略趋于稳定,获得收益也趋于稳定。各VPP均不能从单方面改变自身策略来获得更多利益,此时各主体策略目标达到最优均衡解。
表1 VPP内部资源
在本发明的一个实施例中,设置以下两种运行策略:
Case1:将电网电价和上网电价作为购售电价,各VPP之间采用非合作博弈独自优化,分别选择响应的目标策略进行优化。
Case2:各VPP之间进化学习竞争者竞标策略,获得理性程度的提高,同时进行电价优化,既采用本发明所提出的博弈模型。
由图5仿真结果可知,通过本发明的算法求解多虚拟电厂不同交易目标动态博弈均衡解是有效且收敛的,能够达到良好的求解精度和计算效率。
在两种场景下,DSO制定购售电价如图6所示,每个时刻VPP存在3中状态,若交易电量大于0,表示VPP向DSO售电,为“售电VPP”;若交易电量小于0,表示VPP向DSO购电,为“购电VPP”;若交易电量等于0,则VPP内部供需平衡,为“自平衡VPP”(图8中为画图方便,交易电量正表示购电,负表示售电)。此时,若其余竞争者策略的变化导致购售电价的变化就可能会导致VPP购售电策略的变化。
从交易行为结果图7-9可以看出,VPP1采用“保守型”以自身最大化收益为目标,VPP1在时段4-19内供大于求,因此向上级运营商售电,其余时段存在电量缺额,因此需向上级运营商购电;VPP2采用“激进型”以最大化向上级运营商购电提升VPP自身运行规模为目标,除时段8-10和16以外,其余时段均向上级运营商进行购电。相比较于VPP1和VPP3,VPP2在负荷高峰时段需求量最大。同时通过RL的调控制定最有利于自身规模扩张的负荷增加量与削减量,加大与上级运营商电量交易占有比率。VPP3采用“稳定型”考虑可再生能源最大化发电,同时保证其供电稳定性,因此建立风光储联合发电单元。从图9可以看出,在经过联合单元出力调度利用储能单元平抑可再生能源出力波动性后,呈现的功率曲线稳定性大大提升。在时段15-18内,由于储能单元容量和每时段充放电功率限制,因此为了使波动方差不越限,产生了一定量的弃风弃光量。在时段9-12和16-18内,因自身电量缺额,需向上级运营商进行购电满足功率平衡。
表2给出了DSO采用Case1电价策略和Case2电价策略的收益对比,表3给出了各VPP采用Case1电价策略和Case2电价策略的收益对比(单位/美元)。
表2 DSO不同电价策略效益对比
表3 VPP不同电价策略收益对比
通过对比可知,DSO通过电价优化后,使各VPP购售电量增加并达到平衡。与采用电网电价交易相比,DSO与VPP之间交易收益提高了7.8%。各VPP在博弈的最优目标运行条件下,可降低运行成本,并达到交易的新平衡。
本发明建立的考虑在有限理性条件下的多虚拟电厂参与的不同交易目标动态博弈模型,通过对整合不同分布式能源的多虚拟电厂系统进行算例分析得出以下结论:1)相较于传统以运行成本最小的单一目标函数,可有效满足未来各虚拟电厂在发展的不同阶段的多元化发展需求;2)建立的有限理性交易行为学习,使VPP通过对内部资源的整合和对竞争者竞标策略的学习,达到理性程度的提高,对于竞标结果也能产生较大影响。通过交易使得各VPP能够在市场运行中获得更大效用;3)将智能算法和下层混合整数规划问题相结合求解得到均衡解,与传统的单独使用智能算法相比大大的节约了求解时间并能够快速收敛。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (8)
1.基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,包括以下步骤:
S1:各虚拟电厂竞标主体根据自身发展阶段的目标需求和研究上层运营商动态定价行为;
S2:采用不同交易目标建模对多虚拟电厂进行交易行为分析;
S3:基于粒子群算法对有限理性交易行为信息进行进化学习,通过学习竞争者策略进而使自身目标逐步变优;
S4:研究多虚拟电厂动态博弈计算流程,提出动态博弈粒子群优化算法结合优化工具箱对博弈模型进行求解。
2.如权利要求1所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,根据VPP发展的不同阶段和竞争对手的决策行为对自身优化决策的影响,建立多VPP有不同目标的决策方案;考虑包括VPP内的购售电电量和分布式电源出力为策略,设计多种不同目标的VPP竞标决策方案,一是考虑到不同区域VPP的用户用电差异和VPP联盟目标追求不同,VPP可选择不同的交易目标和决策方案;二是多VPP在进行博弈过程中,各主体会综合其余竞争者的市场交易行为后可改变自身的行为规则,进行目标的转换;在进行博弈的过程中,通过进化学习会不断获取竞争者的信息;随着自身获取信息的不断增多所作出的策略也会逐步变优;最终直至任何主体都不能从单方面从改变竞标策略中获益,既达到稳定策略。
4.如权利要求3所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,VPP的博弈策略考虑为VPP内各个分布式电源的出力计划、储能的充放电功率、可响应负荷的响应功率、风光等清洁能源的输出功率;VPP根据目前已经掌握的有限的市场信息和其余VPP竞标历史数据得到系统负荷需求量以及竞争者可能的竞标策略,来实现最大化自身效用。
5.如权利要求1所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,步骤S2中不同交易目标包括:
保守型VPP的博弈目标函数考虑为最大化自身效益,充分整合内部资源,包括分布式电源发电成本、储能成本、可响应负荷补偿成本;决策出最优竞标策略;VPP竞标目标函数为:
式中:λD-V,sell、λD-V,buy分别为VPP向DSO的售电电价和购电电价;分别为VPP向DSO购买的电量和出售的电量;N为VPP总数集合,VPPi∈N;Ni为第i个VPP中含有的DER集合;为DER中MT的发电成本,aj、bj、cj为DERj中MT机组的成本系数,为MT的输出功率;为ESS运行成本,为ESS的调度成本系数,为ESS的充放电功率,其值大于0表示放电,小于0表示充电;为可响应负荷成本,为响应负荷补偿电价,为响应负荷的响应功率;θ为政府对新能源发电的补贴电价;分别为光伏、风电出力功率;
激进型VPP的博弈目标函数考虑为市场占有率规模最大化,在整合内部资源满足各单元约束出力的前提下,最大化在市场中进行电量竞标,以尽可能多的在DSO层面得到交换电量,获得最大功率购买量;VPP竞标目标函数为:
竞标约束:
在获得DSO统一购售电价和交易电量后,在满足激进型VPP目标函数前提下,对VPP内部其他可调度资源进行协调,满足剩余的负荷电量缺额,电量缺额调度原则满足成本最低原则;
稳定型VPP的博弈目标函数考虑为可再生能源的最大化利用,可再生能源考虑为风电、光伏,同时充分考虑可再生能源发电的不确定性,利用ESS平抑可再生能源出力的波动性对系统供电稳定性的影响;其次再考虑VPP内部其他各单元出力情况和向上级DSO进行电量竞标,VPP竞标目标函数为:
6.如权利要求5所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,所述VPP在进行决策响应时,无论执行哪一种交易目标,都需要满足的约束条件包括:
功率平衡约束:
燃气轮机组约束:
可响应负荷约束:
可再生能源出力约束:
储能约束:
7.如权利要求6所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,所述交易在DSO层面考虑为:
DSO为VPP制定售电电价λD-V,sell和购电电价λD-V,buy;DSO以最大化自身效益作为博弈目标,包括与上级电网和下层多VPP之间的购售电成本和收益;
maxuDSO=FVPP-FDSO
式中:FVPP为DSO与VPP之间进行交易获得的收益;FDSO为DSO与上级电网进行交易所花费的成本;各项的具体表达式如下:
8.如权利要求7所述的基于有限理性的多虚拟电厂动态博弈交易行为分析方法,其特征在于,获得竞争者竞标数据后该竞争者策略判断指标系数如下:
式中:ξ1,ξ2,ξ3为对竞争者策略判断的指标系数,当ξ1≥ξlim1时判断为“激进型”,当ξ2≤ξlim2时判断为“稳定型”,当ξ3≤ξlim3时判断为“保守型”,其中ξlim1,ξlim2,ξlim3为指标系数判断常量,通过各VPP竞标历史数据获得;若在一轮博弈结束策略指标判断不在指标系数判断常量限值内,则不会输出策略结果,继续开始下一轮迭代;
各VPP通过进化学习对博弈策略的选取,决定该VPP的总效益;每种策略对应于该策略下的交易行为方式,用ρmn表示VPPi交易行为方式x(u)从策略m到策略n的条件切换状态;在获取竞争对手竞标信息后,VPPi学习进化自身策略;如果策略n比策略m具有更高的效用,那么策略m将学习策略n并作出改变;用方程表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010832963.3A CN112001752A (zh) | 2020-08-18 | 2020-08-18 | 基于有限理性的多虚拟电厂动态博弈交易行为分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010832963.3A CN112001752A (zh) | 2020-08-18 | 2020-08-18 | 基于有限理性的多虚拟电厂动态博弈交易行为分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112001752A true CN112001752A (zh) | 2020-11-27 |
Family
ID=73473888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010832963.3A Pending CN112001752A (zh) | 2020-08-18 | 2020-08-18 | 基于有限理性的多虚拟电厂动态博弈交易行为分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112001752A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113112287A (zh) * | 2021-03-18 | 2021-07-13 | 国电南瑞南京控制系统有限公司 | 一种基于微电网调度的虚拟电厂能源管理方法及装置 |
CN113610357A (zh) * | 2021-07-16 | 2021-11-05 | 远景智能国际私人投资有限公司 | 虚拟电厂的交易流程管理方法、装置、设备及介质 |
CN115796406B (zh) * | 2023-02-13 | 2023-04-18 | 浙江浙能能源服务有限公司 | 一种用于虚拟电厂的优化调节方法及系统 |
CN115995850A (zh) * | 2023-03-06 | 2023-04-21 | 华北电力大学 | 一种虚拟电厂群的协同调度优化方法和装置 |
CN117318166A (zh) * | 2023-08-09 | 2023-12-29 | 河海大学 | 考虑故障风险下基于合作博弈的多虚拟电厂低碳调度方法 |
CN117541300A (zh) * | 2024-01-08 | 2024-02-09 | 国网浙江省电力有限公司宁波供电公司 | 虚拟电厂交易管理方法、系统、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006158189A (ja) * | 2004-11-04 | 2006-06-15 | Tokyo Electric Power Co Inc:The | 電力取引と分散電源を組み合わせたコスト最小化電力管理システム |
CN108734396A (zh) * | 2018-05-16 | 2018-11-02 | 国网山东省电力公司电力科学研究院 | 一种基于需求侧竞价和多能互补的虚拟电厂调度优化方法 |
CN108921368A (zh) * | 2018-05-03 | 2018-11-30 | 东南大学 | 基于虚拟发电厂的均衡合作博弈控制器 |
CN108960510A (zh) * | 2018-07-04 | 2018-12-07 | 四川大学 | 一种基于两阶段随机规划的虚拟电厂优化交易策略模型 |
CN109902884A (zh) * | 2019-03-27 | 2019-06-18 | 合肥工业大学 | 一种基于主从博弈策略的虚拟电厂优化调度方法 |
CN110135613A (zh) * | 2018-10-23 | 2019-08-16 | 上海交通大学 | 一种基于纳什谈判的多虚拟电厂协同优化运行方案 |
CN111402015A (zh) * | 2020-03-17 | 2020-07-10 | 南京工业大学 | 一种基于购售风险的虚拟电厂双层竞标方法及系统 |
-
2020
- 2020-08-18 CN CN202010832963.3A patent/CN112001752A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006158189A (ja) * | 2004-11-04 | 2006-06-15 | Tokyo Electric Power Co Inc:The | 電力取引と分散電源を組み合わせたコスト最小化電力管理システム |
CN108921368A (zh) * | 2018-05-03 | 2018-11-30 | 东南大学 | 基于虚拟发电厂的均衡合作博弈控制器 |
CN108734396A (zh) * | 2018-05-16 | 2018-11-02 | 国网山东省电力公司电力科学研究院 | 一种基于需求侧竞价和多能互补的虚拟电厂调度优化方法 |
CN108960510A (zh) * | 2018-07-04 | 2018-12-07 | 四川大学 | 一种基于两阶段随机规划的虚拟电厂优化交易策略模型 |
CN110135613A (zh) * | 2018-10-23 | 2019-08-16 | 上海交通大学 | 一种基于纳什谈判的多虚拟电厂协同优化运行方案 |
CN109902884A (zh) * | 2019-03-27 | 2019-06-18 | 合肥工业大学 | 一种基于主从博弈策略的虚拟电厂优化调度方法 |
CN111402015A (zh) * | 2020-03-17 | 2020-07-10 | 南京工业大学 | 一种基于购售风险的虚拟电厂双层竞标方法及系统 |
Non-Patent Citations (4)
Title |
---|
HONGJUN GAO等: "Bounded Rationality Based Multi-VPP Trading in Local Energy Markets: A Dynamic Game Approach with Different Trading Targets" * |
刘继春等: "考虑多个虚拟发电厂参与的多阶段市场交易方法" * |
周博;吕林;高红均;刘俊勇;陈庆攀;谭心怡;: "多虚拟电厂日前鲁棒交易策略研究" * |
周博等: "考虑热电联合调度的虚拟电厂交易策略研究" * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113112287A (zh) * | 2021-03-18 | 2021-07-13 | 国电南瑞南京控制系统有限公司 | 一种基于微电网调度的虚拟电厂能源管理方法及装置 |
CN113610357A (zh) * | 2021-07-16 | 2021-11-05 | 远景智能国际私人投资有限公司 | 虚拟电厂的交易流程管理方法、装置、设备及介质 |
CN115796406B (zh) * | 2023-02-13 | 2023-04-18 | 浙江浙能能源服务有限公司 | 一种用于虚拟电厂的优化调节方法及系统 |
CN115995850A (zh) * | 2023-03-06 | 2023-04-21 | 华北电力大学 | 一种虚拟电厂群的协同调度优化方法和装置 |
CN117318166A (zh) * | 2023-08-09 | 2023-12-29 | 河海大学 | 考虑故障风险下基于合作博弈的多虚拟电厂低碳调度方法 |
CN117541300A (zh) * | 2024-01-08 | 2024-02-09 | 国网浙江省电力有限公司宁波供电公司 | 虚拟电厂交易管理方法、系统、设备及存储介质 |
CN117541300B (zh) * | 2024-01-08 | 2024-06-04 | 国网浙江省电力有限公司宁波供电公司 | 虚拟电厂交易管理方法、系统、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112001752A (zh) | 基于有限理性的多虚拟电厂动态博弈交易行为分析方法 | |
CN110276698B (zh) | 基于多智能体双层协同强化学习的分布式可再生能源交易决策方法 | |
CN110147907B (zh) | 一种考虑网络约束的虚拟电厂双层优化模型的构建方法 | |
Li et al. | Trading strategy and benefit optimization of load aggregators in integrated energy systems considering integrated demand response: A hierarchical Stackelberg game | |
CN112467722B (zh) | 一种考虑电动汽车充电站的主动配电网源-网-荷-储协调规划方法 | |
Maity et al. | Simulation and pricing mechanism analysis of a solar-powered electrical microgrid | |
Chen et al. | Research on day-ahead transactions between multi-microgrid based on cooperative game model | |
CN113592648B (zh) | 一种综合能源系统的多主体交易方法及系统 | |
CN111008739B (zh) | 一种热电联产虚拟电厂优化调控及收益分配方法及系统 | |
CN109919450B (zh) | 解决综合智慧能源系统调度的博弈优化方法 | |
CN112217195B (zh) | 一种基于gru多步预测技术的云储能充放电策略形成方法 | |
CN112564102B (zh) | 多微网负荷优化调度方法和系统 | |
Gao et al. | Bounded rationality based multi-VPP trading in local energy markets: a dynamic game approach with different trading targets | |
CN111311012A (zh) | 基于多智能体的微电网电力市场双层竞价优化方法 | |
CN114256836A (zh) | 一种新能源电站共享储能的容量优化配置方法 | |
CN115147245B (zh) | 一种工业负荷参与调峰辅助服务的虚拟电厂优化调度方法 | |
CN115313520A (zh) | 分布式能源系统博弈优化调度方法、系统、设备及介质 | |
CN113690875B (zh) | 一种微网实时互动等值模型建立方法 | |
CN117332937A (zh) | 一种考虑需求响应的多能互补虚拟电厂经济调度方法 | |
CN110322092B (zh) | 一种源网荷多层博弈经济调度方法 | |
CN116961008A (zh) | 计及电力弹簧与负荷需求响应的微电网容量双层优化方法 | |
Wu et al. | Peer-to-peer energy trading optimization for community prosumers considering carbon cap-and-trade | |
CN114188940B (zh) | 一种考虑多主体利益的主动配电网调度方法 | |
CN116307029A (zh) | 促进多虚拟网格间源储荷协调的双层优化调度方法及系统 | |
CN116432862A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201127 |