CN114881749A - 一种计及风险约束的虚拟电厂参与主辅市场的投标方法 - Google Patents
一种计及风险约束的虚拟电厂参与主辅市场的投标方法 Download PDFInfo
- Publication number
- CN114881749A CN114881749A CN202210461034.5A CN202210461034A CN114881749A CN 114881749 A CN114881749 A CN 114881749A CN 202210461034 A CN202210461034 A CN 202210461034A CN 114881749 A CN114881749 A CN 114881749A
- Authority
- CN
- China
- Prior art keywords
- power plant
- virtual power
- market
- bidding
- constraint
- 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
- 238000000034 method Methods 0.000 title claims abstract description 71
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000003860 storage Methods 0.000 claims abstract description 5
- 238000004146 energy storage Methods 0.000 claims description 33
- 238000007599 discharging Methods 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 20
- 230000005611 electricity Effects 0.000 claims description 14
- 230000008901 benefit Effects 0.000 claims description 13
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 239000000370 acceptor Substances 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 238000012725 vapour phase polymerization Methods 0.000 abstract description 6
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 241000756943 Codonopsis Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 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/06—Buying, selling or leasing transactions
- G06Q30/08—Auctions
-
- 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/0635—Risk analysis of enterprise or organisation activities
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards
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)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Game Theory and Decision Science (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (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
本发明涉及一种计及风险约束的含电动汽车虚拟电厂同时参与主辅市场的投标方法，虚拟电厂含EV，包括：确定含EV的虚拟电厂参与电力市场的交易调度流程；基于交易调度流程和风险条件价值CVaR理论建立计及风险约束的虚拟电厂同时参与主辅市场的投标模型；求解计及风险约束的虚拟电厂同时参与主辅市场的投标模型；基于投标模型的求解结果将虚拟电厂参与主辅市场获得的运营收益合理分配，从而获得虚拟电厂同时参与主辅市场的投标方案。还公开了对应的系统、终端设备以及计算机可读存储介质。为单个容量较小的需求侧资源参与电力市场，为VPP同时参与主辅市场，为VPP运营商根据自身风险偏好程度灵活衡量风险与风险收益之间的关系提供了手段。
Description
技术领域
本发明涉及含电动汽车虚拟电厂参与电力市场的优化领域，尤其是涉及一种计及风险约束的含电动汽车虚拟电厂同时参与主辅市场的投标方法。
背景技术
近年来在可持续发展的指导理念下，风电及光伏等清洁能源快速发展，但风电、光伏出力的不确定性给电力系统的安全稳定运营带来了巨大的挑战。在此背景下，需加快调峰辅助服务市场的建设，在提高电力系统调控能力的同时促进可再生能源的消纳。电动汽车(Electric Vehicle，EV)以其节能、零排放的优势被广泛使用，据国家工业和信息化部称，我国2030年EV的保有量可能约为6000万辆，大量EV的并网，既提高了我国电力市场资源的利用率，也加快了实现“双碳”目标的步伐，但单体EV因其容量小且地理位置分散无法直接参与主能量市场及辅助服务市场。而虚拟电厂(Virtual Power Plant，VPP)则是一种通过先进信息通信技术和软件系统，实现分布式电源、储能系统、可控负荷、电动汽车等分布式能源的聚合和协调调度，以作为一个特殊的虚拟市场主体参与电力市场。VPP的出现，其不仅能够有效聚合各类可再生能源及分散的EV共同参与电力市场，还能协调EV应对风电、光伏机组的“弃风、弃光”现象，这种运营思路在我国有着非常大的市场潜力。
公开于该背景技术部分的信息仅仅旨在加深对本发明总体背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成本领域技术人员所公知的现有技术。
发明内容
本发明所要解决的技术问题是：提供一种计及风险约束的含电动汽车的虚拟电厂同时参与主能量市场与调峰辅助服务市场的投标方法，即虚拟电厂运营商在风电、光伏出力风险水平可接受的前提下，追求含电动汽车的虚拟电厂在主辅电力市场的收益最大。
为了达到上述目的，本发明所采用的技术方案是：一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，包括：
S1，确定含EV的虚拟电厂参与电力市场的交易调度流程；
S2，基于所述交易调度流程和风险条件价值(Conditional Value at Risk,CVaR)理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S3，求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
优选的，所述S1确定含EV的虚拟电厂参与电力市场的交易调度流程基于所述虚拟电厂能够聚合EV、风电机组、光伏机组、柔性负荷与储能装置资源作为价格接受者参与主能量市场和调峰辅助服务市场，所述交易调度流程包括两个阶段：
第一阶段为日前运行阶段，在所述第一阶段内确定所述虚拟电厂参与日前调峰辅助服务市场的调峰时段、容量及柔性负荷的出力安排和运行基线，所述第一阶段与具体运行场景无关；
第二阶段为基于风光实际出力场景的实时运行阶段，在所述第二阶段内确定所述虚拟电厂中EV的充放电功率、储能的充放电功率、弃风弃光功率、参与调峰市场的实际功率与直接交易功率，本阶段的决策变量与具体场景相关。
优选的，所述S2计及风险约束的虚拟电厂同时参与主辅市场的投标模型包括目标函数以及约束条件，所述目标函数为所述虚拟电厂在日前主能量市场与调峰辅助服务市场中获得的总收益；所述约束条件包括虚拟电厂购售电约束、虚拟电厂外特性约束、虚拟电厂功率平衡约束、柔性负荷约束、EV约束、储能设备约束、调峰投标约束、调峰实际出力约束、CVaR相关约束。
优选的，所述目标函数如下式(1)所示：
式中，T为虚拟电厂的一个投标周期，取24h；Nω为风电、光伏机组出力的场景总数；ρω为场景ω下的概率；β为风险偏好系数，代表虚拟电厂的运营商对风险收益的偏好程度，β∈[0,1)，β值越大则代表虚拟电厂的运营商对风险收益越厌恶，其设计的定价策略就越保守；δ为CVaR值；α为置信水平；ξω式引入的辅助变量，为场景ω下虚拟电厂运营总收益与δ之间的差值；
所述目标函数主要由两部分组成，第一部分为虚拟电厂的运营收益，包括虚拟电厂与电网直接交易的收益Bgrid,tω、参与调峰市场的收益Bp,tω、调峰惩罚成本Cp,tω、柔性负荷调度成本Cload,tω与弃风弃光惩罚成本Cabon,tω；第二部分为CVaR值与风险偏好系数β的乘积；其中：
(1)与电网直接交易的收益Bgrid,tω为：
式中，
分别为t时段VPP在日前主能量市场中的购、售电价；
分别为场景ω下t时段VPP在日前著能量市场中的购、售电功率；Δt为调度步长，取1h；
(2)柔性负荷调度成本Cload,tω为：
Cload,tω＝πload,tPload,tωΔt (3)；
式中，πload,t为t时段柔性负荷调度的补偿价格；Pload,tω为场景ω下t时段柔性负荷的调度功率；Δt为调度步长，取1h；
(3)弃风弃光惩罚成本Cabon,tω为：
Cabon,tω＝πabon,t[Pwindb,tω+Ppvb,tω]Δt (4)；
式中，πabon,t为t时段虚拟电厂弃风、弃光惩罚电价；Pwindb,tω为场景ω下t时段风电机组的弃风功率；Ppvb,tω为场景ω下t时段光伏机组的弃光功率；Δt为调度步长，取1h；
(4)参与调峰辅助服务市场的收益Bp,tω为：
式中，Ppa,t为t时段VPP在调峰辅助服务市场的日前投标容量；Ppd,tω为场景ω下t时段VPP的调峰实际出力；
为t时段VPP在调峰辅助服务市场中的补偿价格；Δt为调度步长，取1h；
(5)调峰惩罚成本Cp,tω为：
当虚拟电厂调峰实际出力未达到投标容量的70％时，应受到调峰市场的惩罚Cp,tω如下，
式中，
为t时段虚拟电厂在调峰辅助服务市场的惩罚价格；Δt为调度步长，取1h。
优选的，所述约束条件包括：
(1)虚拟电厂购售电约束，表示虚拟电厂无法同时进行购售电行为：
(2)虚拟电厂外特性：
Pbase,t＝Dt-Pwind,t-Ppv,t (10)；
其中，Pbase,t为虚拟电厂根据风光预测值上报的t时段的运行基线；Pwind,t、Ppv,t分别为t时段的风、光预测值；Dt为柔性负荷调度前t时段的负荷功率；Pout,tω为场景ω下t时段VPP实时运行的外特性功率；Pwinda,tω、Ppva,tω分别为场景ω下t时段的风、光实际出力；Dt'为柔性负荷调度后t时段的负荷功率；
分别为EV在场景ω下t时段的充、放电功率；
分别为储能设备在场景ω下t时段的充、放电功率；
(3)虚拟电厂的功率平衡约束：
(4)柔性负荷约束，柔性负荷具有可转移特性，在峰时段削减的负荷总功率与平谷时段增加的负荷总功率相等：
Dt-D′t＝Pla,t-Plc,t (14)；
0≤Pla.t≤Pload,max t∈TV (15)；
0≤Plc.t≤γDt t∈TP (16)；
式中，TV为平谷时段；TP为峰时段；Pla,t为平谷时段中t时段柔性负荷增加的功率，其功率值不应大于最大值Pload,max；Plc,t为在峰时段中t时段柔性负荷削减的负荷总功率；γ是虚拟电厂中柔性负荷的比例；
(5)EV约束：
Eev,min≤Eev,tω≤Eev,max (21)；
式中，Pev,max为EV的最大充、放电功率；Eev,tω为EVt时段的荷电量；
为EV的充、放电效率；Qev为EV车载电池的最大容量；Eev,min、Eev,max为EV车载电池所允许的最小、最大荷电量；
分别为EV在场景ω下t时段的充、放电功率；
(6)储能约束：
Eess,min≤Eess,tω≤Eess,max (26)；
Eess,1ω＝Eess,24ω (27)；
式中，Pess,max为储能的最大充、放电功率；Eess,tω、Eess,(t+1)ω、Eess,1ω、Eess,24ω分别为储能t时段、t+1时段、1时、24时的荷电量；
为储能的充、放电效率；Qess为储能设备的最大容量；Eess,min、Eess,max为储能设备允许的最小、最大荷电量；
分别为储能设备在场景ω下t时段的充、放电功率；
(7)调峰投标约束：
Ppa,minu(t)≤Ppa,tω≤Mu(t)t∈TV (28)；
Ppa,tω＝0 t∈TP (29)；
式中，u(t)为布尔变量，表示该时段虚拟电厂是否参与调峰市场，值为1时表示VPP参与调峰市场，值为0时表示VPP不参与调峰市场；M表示一个很大的正数；Ppa,min为VPP日前向调峰市场投标的最小容量；
(8)调峰实际出力约束：
0≤Ppd,tω≤Pout,tω-Pbase,tω+M(1-u(t))t∈TV (30)；
Ppd,tω≤u(t)M t∈TV (31)；
Ppd,tω＝0 t∈TP (32)；
(8)CVaR相关约束：
对于一个离散的收益分布，当其置信水平为α时，CVaR近似为小概率1-α情景集合的期望损失；CVaR及其相关约束可转化为如下线性问题进行求解：
ξω≥0 (34)；
式中，ξω为引入的辅助变量；δ的最优值为对应的风险价值，即VPP运营收益小于或等于ξ的概率小于或等于1-α的最大收益。
优选的，所述S3包括：
S31，将所述目标函数及所述约束条件共同构成的非线性模型转化为混合整数线性规划问题，即模型的线性化；
S32，通过商业求解器Gurobi求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；其中：
所述S31包括：
对于式(9)、(19)、(24)中连续变量的乘积为0，引入布尔变量；对于式(5)、(8)的分段函数，线性化包括：
0≤Fp,tω≤Ppa,t (35)；
Fp,tω≤Ppd,tω (36)；
Ppd,tω-0.7Ppa,t≥(atω-1)M (37)；
Ppd,tω-0.7Ppa,t≥(1-btω)M (38)；
atω+btω＝1 (39)；
Fp,tω≤atωM (40)；
Ktω≥Ppa,t-Ppd,tω+(btω-1)M (41)；
Ktω≥0 (42)；
式中，atω、btω均为布尔变量。
优选的，所述S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案包括：基于Shapley值法将虚拟电厂参与电力市场获得的总运营收益合理分配，具体步骤如下：对于一个由J个参与者组成的联盟，将成员集合用Ω＝{1,2,…,J}表示，Ω的任意子集称为一个子联盟。利用Shapley值法对虚拟电厂的市场收益进行分配，得到各成员的收益为：
式中，Cj为虚拟电厂内成员j分配的收益；sj为包含成员j的子联盟；|sj|为子联盟sj中的成员数；C(sj)、C(sj/j)分别为包含成员j与不包含成员j的子联盟的运行收益。
本发明的目的还在于提供一种计及风险约束的虚拟电厂同时参与主辅市场的投标系统，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，包括：
交易调度流程确定模块，用于确定含EV的虚拟电厂参与电力市场的交易调度流程；
投标模型建立模块，用于基于所述交易调度流程和风险条件价值(ConditionalValue at Risk,CVaR)理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
模型求解模块，用于求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
投标模块：用于基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
本发明的目的还在于提供一种终端设备，包括：一个或多个处理器；存储器，与所述处理器耦接，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如上所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法。
本发明的目的还在于提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如上所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法。
本发明的有益效果为：
1、计及风险约束的含电动汽车虚拟电厂同时参与主辅市场的投标方法，为单个容量较小的需求侧资源参与电力市场提供了途径，也为VPP同时参与主能量市场与辅助服务市场提供了运营思路。
2、针对风电机组、光伏机组出力的不确定性引入了CVaR理论，为VPP运营商根据自身风险偏好程度，灵活衡量风险与风险收益之间的关系提供了重要手段。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中VPP参与主能量市场及辅助服务市场的交易流程图；
图2为本发明实施例的终端设备结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
需要说明的是，当元件被称为“固定于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。本文所使用的术语“垂直的”、“水平的”、“左”、“右”以及类似的表述只是为了说明的目的，并不表示是唯一的实施方式。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
实施例一
本实施例的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，包括：
S1，确定含EV的虚拟电厂参与电力市场的交易调度流程；
S2，基于所述交易调度流程和风险条件价值(Conditional Value at Risk,CVaR)理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S3，求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
作为优选的实施方式，所述S1确定含EV的虚拟电厂参与电力市场的交易调度流程基于所述虚拟电厂能够聚合EV、风电机组、光伏机组、柔性负荷与储能装置资源作为价格接受者参与主能量市场和调峰辅助服务市场，所述交易调度流程包括两个阶段：
第一阶段为日前运行阶段，在所述第一阶段内确定所述虚拟电厂参与日前调峰辅助服务市场的调峰时段、容量及柔性负荷的出力安排和运行基线，所述第一阶段与具体运行场景无关；所述确定所述虚拟电厂参与日前调峰辅助服务市场的调峰时段、容量及柔性负荷的出力安排和运行基线的方法包括：电力市场交易中心向虚拟电厂先下发日前能量市场的交易价格及调峰辅助服务市场的价格，虚拟电厂作为价格接受者在整合内部的风电机组、光伏机组、储能设备、柔性负荷、EV资源后，日前向调度机构申报参与日前调峰辅助服务市场的调峰时段、投标容量、柔性负荷的出力安排和运行基线。为保证市场公平，VPP上报的运行基线为日前预测净负荷曲线，电力市场交易中心则根据所述虚拟电厂上报的运行基线与实际运行曲线差值确定所述虚拟电厂的调峰实际出力；
第二阶段为基于风光实际出力场景的实时运行阶段，在所述第二阶段内确定所述虚拟电厂中EV的充放电功率、储能的充放电功率、弃风弃光功率、参与调峰市场的实际功率与直接交易功率，本阶段的决策变量与具体场景相关；所述确定所述虚拟电厂中EV的充放电功率、储能的充放电功率、弃风弃光功率、参与调峰市场的实际功率与直接交易功率包括：所述虚拟电厂获得当日的风电机组、光伏机组的出力曲线；基于获得的信息确定VPP中EV的充放电功率、储能的充放电功率、弃风弃光功率、参与调峰辅助服务市场的实际功率与直接交易功率；将确定的信息提交电力市场交易中心，交易中心确认信息无误后，完成实时功率出清。因风光出力的不确定性的影响，本阶段的决策变量与具体场景有关。
基于上述流程，虚拟电厂参与主能量市场及辅助服务市场的交易流程如附图1所示。
作为优选的实施方式，所述S2计及风险约束的虚拟电厂同时参与主辅市场的投标模型包括目标函数以及约束条件，所述目标函数为所述虚拟电厂在日前主能量市场与调峰辅助服务市场中获得的总收益；所述约束条件包括虚拟电厂购售电约束、虚拟电厂外特性约束、虚拟电厂功率平衡约束、柔性负荷约束、EV约束、储能设备约束、调峰投标约束、调峰实际出力约束、CVaR相关约束。
作为优选的实施方式，所述目标函数如下式(1)所示：
式中，T为虚拟电厂的一个投标周期，取24h；Nω为风电、光伏机组出力的场景总数；ρω为场景ω下的概率；β为风险偏好系数，代表虚拟电厂的运营商对风险收益的偏好程度，β∈[0,1)，β值越大则代表虚拟电厂的运营商对风险收益越厌恶，其设计的定价策略就越保守；δ为CVaR值；α为置信水平；ξω式引入的辅助变量，为场景ω下虚拟电厂运营总收益与δ之间的差值。
作为优选的实施方式，所述目标函数主要由两部分组成，第一部分为虚拟电厂的运营收益，包括虚拟电厂与电网直接交易的收益Bgrid,tω、参与调峰市场的收益Bp,tω、调峰惩罚成本Cp,tω、柔性负荷调度成本Cload,tω与弃风弃光惩罚成本Cabon,tω；第二部分为CVaR值与风险偏好系数β的乘积；其中：
(1)与电网直接交易的收益Bgrid,tω为：
式中，
分别为t时段VPP在日前主能量市场中的购、售电价；
分别为场景ω下t时段VPP在日前著能量市场中的购、售电功率；Δt为调度步长，取1h；
(2)柔性负荷调度成本Cload,tω为：
Cload,tω＝πload,tPload,tωΔt (46)；
式中，πload,t为t时段柔性负荷调度的补偿价格；Pload,tω为场景ω下t时段柔性负荷的调度功率；Δt为调度步长，取1h；
(3)弃风弃光惩罚成本Cabon,tω为：
Cabon,tω＝πabon,t[Pwindb,tω+Ppvb,tω]Δt (47)；
式中，πabon,t为t时段虚拟电厂弃风、弃光惩罚电价；Pwindb,tω为场景ω下t时段风电机组的弃风功率；Ppvb,tω为场景ω下t时段光伏机组的弃光功率；Δt为调度步长，取1h；
(4)参与调峰辅助服务市场的收益Bp,tω为：
式中，Ppa,t为t时段VPP在调峰辅助服务市场的日前投标容量；Ppd,tω为场景ω下t时段VPP的调峰实际出力；
为t时段VPP在调峰辅助服务市场中的补偿价格；Δt为调度步长，取1h；
(5)调峰惩罚成本Cp,tω为：
当虚拟电厂调峰实际出力未达到投标容量的70％时，应受到调峰市场的惩罚Cp,tω如下，
式中，
为t时段虚拟电厂在调峰辅助服务市场的惩罚价格；Δt为调度步长，取1h。
作为优选的实施方式，目标函数的约束条件有虚拟电厂购售电约束、虚拟电厂外特性约束、虚拟电厂功率平衡约束、柔性负荷约束、EV约束、储能设备约束、调峰投标约束、调峰实际出力约束、CVaR相关约束；其中：
(1)虚拟电厂购售电约束，表示虚拟电厂无法同时进行购售电行为：
(2)虚拟电厂外特性：
Pbase,t＝Dt-Pwind,t-Ppv,t (53)；
其中，Pbase,t为虚拟电厂根据风光预测值上报的t时段的运行基线；Pwind,t、Ppv,t分别为t时段的风、光预测值；Dt为柔性负荷调度前t时段的负荷功率；Pout,tω为场景ω下t时段VPP实时运行的外特性功率；Pwinda,tω、Ppva,tω分别为场景ω下t时段的风、光实际出力；Dt'为柔性负荷调度后t时段的负荷功率；
分别为EV在场景ω下t时段的充、放电功率；
分别为储能设备在场景ω下t时段的充、放电功率；
(3)虚拟电厂的功率平衡约束：
(4)柔性负荷约束，柔性负荷具有可转移特性，在峰时段削减的负荷总功率与平谷时段增加的负荷总功率相等：
Dt-D′t＝Pla,t-Plc,t (57)；
0≤Pla.t≤Pload,max t∈TV (58)；
0≤Plc.t≤γDtt∈TP (59)；
式中，TV为平谷时段；TP为峰时段；Pla,t为平谷时段中t时段柔性负荷增加的功率，其功率值不应大于最大值Pload,max；Plc,t为在峰时段中t时段柔性负荷削减的负荷总功率；γ是虚拟电厂中柔性负荷的比例；
(5)EV约束：
Eev,min≤Eev,tω≤Eev,max (64)；
式中，Pev,max为EV的最大充、放电功率；Eev,tω为EVt时段的荷电量；
为EV的充、放电效率；Qev为EV车载电池的最大容量；Eev,min、Eev,max为EV车载电池所允许的最小、最大荷电量；
分别为EV在场景ω下t时段的充、放电功率；
(6)储能约束：
Eess,min≤Eess,tω≤Eess,max (69)；
Eess,1ω＝Eess,24ω (70)；
式中，Pess,max为储能的最大充、放电功率；Eess,tω、Eess,(t+1)ω、Eess,1ω、Eess,24ω分别为储能t时段、t+1时段、1时、24时的荷电量；
为储能的充、放电效率；Qess为储能设备的最大容量；Eess,min、Eess,max为储能设备允许的最小、最大荷电量；
分别为储能设备在场景ω下t时段的充、放电功率；
(6)调峰投标约束：
Ppa,minu(t)≤Ppa,tω≤Mu(t) t∈TV (71)；
Ppa,tω＝0 t∈TP (72)；
式中，u(t)为布尔变量，表示该时段虚拟电厂是否参与调峰市场，值为1时表示VPP参与调峰市场，值为0时表示VPP不参与调峰市场；M表示一个很大的正数；Ppa,min为VPP日前向调峰市场投标的最小容量；
(7)调峰实际出力约束：
0≤Ppd,tω≤Pout,tω-Pbase,tω+M(1-u(t)) t∈TV (73)；
Ppd,tω≤u(t)M t∈TV (74)；
Ppd,tω＝0 t∈TP (75)；
(8)CVaR相关约束：
本发明利用风险条件价值CVaR度量风光出力不确定性给调峰市场收益带来的影响。条件风险价值方法不依赖于投资回报符合正态分布的假设，而电力市场中由于各种因素的影响，电价以及负荷的波动总是难以满足正态分布的要求，所以采用条件风险价值测量和管理虚拟电厂的风险，具有重要的现实意义，将模型转换为线性规划求解，避免了其他模型复杂的求解过程以及不可能收敛的情况。
对于本实施例，基于之前分析，第一阶段已确定调峰时段与投标容量，在实时运行阶段为满足调峰实际出力的要求会出现弃风、弃光的现象。同时，直接交易收益与调峰收益是相对独立的两部分。因此，调峰风险收益应考虑调峰收益、调峰惩罚与弃风弃光成本。对于一个离散的收益分布，当其置信水平为α时，CVaR近似为小概率1-α情景集合的期望损失。CVaR及其相关约束可转化为如下线性问题进行求解：
ξω≥0 (77)；
式中，ξω为引入的辅助变量；δ的最优值为对应的风险价值，即VPP运营收益小于或等于ξ的概率小于或等于1-α的最大收益。
作为优选的实施方式，所述S3包括：
S31，将所述目标函数及所述约束条件共同构成的非线性模型转化为混合整数线性规划问题，即模型的线性化；
S32，通过商业求解器Gurobi求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型。
作为优选的实施方式，所述S31包括：
需要注意的是，模型的非线性部分主要由两部分组成：式(9)、(19)、(24)中连续变量的乘积为0；目标函数式(5)、(8)的分段函数；对于连续变量乘积为0，引入布尔变量便可解决，此处不再赘述；这里主要就(5)、(8)两式的线性化进行说明。
0≤Fp,tω≤Ppa,t (78)；
Fp,tω≤Ppd,tω (79)；
Ppd,tω-0.7Ppa,t≥(atω-1)M (80)；
Ppd,tω-0.7Ppa,t≥(1-btω)M (81)；
atω+btω＝1 (82)；
Fp,tω≤atωM (83)；
Ktω≥Ppa,t-Ppd,tω+(btω-1)M (84)；
Ktω≥0 (85)；
式中，atω、btω均为布尔变量。
作为优选的实施方式，所述S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案包括：在VPP的统一管理下，各成员以追求经济效益最大化为目标组成联盟，进行合作博弈共同在日前能量市场与调峰辅助服务市场获得收益，但VPP应考虑各成员对总运营收益的贡献程度，需基于Shapley值法将VPP参与电力市场获得的总运营收益合理分配。
目前常见的利益分配法为核仁法与Shapley值法。核仁法使成员形成的所有联盟的不满意度达到最小，而Shapley值法则关注各成员对联盟的边际贡献。本实施例中，选择Shapley值法，具体步骤如下：对于一个由J个参与者组成的联盟，将成员集合用Ω＝{1,2,…,J}表示，Ω的任意子集称为一个子联盟。利用Shapley值法对虚拟电厂的市场收益进行分配，得到各成员的收益为：
式中，Cj为虚拟电厂内成员j分配的收益；sj为包含成员j的子联盟；|sj|为子联盟sj中的成员数；C(sj)、C(sj/j)分别为包含成员j与不包含成员j的子联盟的运行收益。
实施例二
一种计及风险约束的虚拟电厂同时参与主辅市场的投标系统，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，包括：
交易调度流程确定模块，用于确定含EV的虚拟电厂参与电力市场的交易调度流程；
投标模型建立模块，用于基于所述交易调度流程和风险条件价值(ConditionalValue at Risk,CVaR)理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
模型求解模块，用于求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
投标模块：用于基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
请参阅图2，某一实施例提供一种终端设备，包括：
一个或多个处理器；
存储器，与所述处理器耦接，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如上所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法。
处理器用于控制该终端设备的整体操作，以完成上述的新能源与常规能源的多目标协调控制方法的全部或部分步骤。存储器用于存储各种类型的数据以支持在该终端设备的操作，这些数据例如可以包括用于在该终端设备上操作的任何应用程序或方法的指令，以及应用程序相关的数据。该存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read-Only Memory，简称EPROM)，可编程只读存储器(Programmable Read-Only Memory，简称PROM)，只读存储器(Read-Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。
在一示例性实施例中，终端设备可以被一个或多个应用专用集成电路(Application Specific 1ntegrated Circuit，简称AS1C)、数字信号处理器(DigitalSignal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device,简称DSPD)、可编程逻辑器件(Programmable Logic Device，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行如上述任一项实施例所述的投标方法，并达到如上述方法一致的技术效果。
在另一示例性实施例中，还提供一种包括计算机程序的计算机可读存储介质，该计算机程序被处理器执行时实现如上述任一项实施例所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法的步骤。例如，该计算机可读存储介质可以为上述包括计算机程序的存储器，上述计算机程序可由终端设备的处理器执行以完成如上述任一项实施例所述的投标方法，并达到如上述方法一致的技术效果。
以上实施例计及风险约束的含电动汽车虚拟电厂同时参与主辅市场的投标方法，为单个容量较小的需求侧资源参与电力市场提供了途径，也为VPP同时参与主能量市场与辅助服务市场提供了运营思路；针对风电机组、光伏机组出力的不确定性引入了CVaR理论，为VPP运营商根据自身风险偏好程度，灵活衡量风险与风险收益之间的关系提供了重要手段。
本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (10)
1.一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，其特征在于，包括：
S1，确定含EV的虚拟电厂参与电力市场的交易调度流程；
S2，基于所述交易调度流程和风险条件价值CVaR理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S3，求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
2.根据权利要求1所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述S1确定含EV的虚拟电厂参与电力市场的交易调度流程基于所述虚拟电厂能够聚合EV、风电机组、光伏机组、柔性负荷与储能装置资源作为价格接受者参与主能量市场和调峰辅助服务市场，所述交易调度流程包括两个阶段：
第一阶段为日前运行阶段，在所述第一阶段内确定所述虚拟电厂参与日前调峰辅助服务市场的调峰时段、容量及柔性负荷的出力安排和运行基线，所述第一阶段与具体运行场景无关；
第二阶段为基于风光实际出力场景的实时运行阶段，在所述第二阶段内确定所述虚拟电厂中EV的充放电功率、储能的充放电功率、弃风弃光功率、参与调峰市场的实际功率与直接交易功率，本阶段的决策变量与具体场景相关。
3.根据权利要求2所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述S2计及风险约束的虚拟电厂同时参与主辅市场的投标模型包括目标函数以及约束条件，所述目标函数为所述虚拟电厂在日前主能量市场与调峰辅助服务市场中获得的总收益；所述约束条件包括虚拟电厂购售电约束、虚拟电厂外特性约束、虚拟电厂功率平衡约束、柔性负荷约束、EV约束、储能设备约束、调峰投标约束、调峰实际出力约束、CVaR相关约束。
4.根据权利要求3所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述目标函数如下式(1)所示：
式中，T为虚拟电厂的一个投标周期，取24h；Nω为风电、光伏机组出力的场景总数；ρω为场景ω下的概率；β为风险偏好系数，代表虚拟电厂的运营商对风险收益的偏好程度，β∈[0,1)，β值越大则代表虚拟电厂的运营商对风险收益越厌恶，其设计的定价策略就越保守；δ为CVaR值；α为置信水平；ξω式引入的辅助变量，为场景ω下虚拟电厂运营总收益与δ之间的差值；
所述目标函数主要由两部分组成，第一部分为虚拟电厂的运营收益，包括虚拟电厂与电网直接交易的收益Bgrid,tω、参与调峰市场的收益Bp,tω、调峰惩罚成本Cp,tω、柔性负荷调度成本Cload,tω与弃风弃光惩罚成本Cabon,tω；第二部分为CVaR值与风险偏好系数β的乘积；其中：
(1)与电网直接交易的收益Bgrid,tω为：
式中，
分别为t时段VPP在日前主能量市场中的购、售电价；
分别为场景ω下t时段VPP在日前著能量市场中的购、售电功率；Δt为调度步长，取1h；
(2)柔性负荷调度成本Cload,tω为：
Cload,tω＝πload,tPload,tωΔt (3)；
式中，πload,t为t时段柔性负荷调度的补偿价格；Pload,tω为场景ω下t时段柔性负荷的调度功率；Δt为调度步长，取1h；
(3)弃风弃光惩罚成本Cabon,tω为：
Cabon,tω＝πabon,t[Pwindb,tω+Ppvb,tω]Δt (4)；
式中，πabon,t为t时段虚拟电厂弃风、弃光惩罚电价；Pwindb,tω为场景ω下t时段风电机组的弃风功率；Ppvb,tω为场景ω下t时段光伏机组的弃光功率；Δt为调度步长，取1h；
(4)参与调峰辅助服务市场的收益Bp,tω为：
式中，Ppa,t为t时段VPP在调峰辅助服务市场的日前投标容量；Ppd,tω为场景ω下t时段VPP的调峰实际出力；
为t时段VPP在调峰辅助服务市场中的补偿价格；Δt为调度步长，取1h；
(5)调峰惩罚成本Cp,tω为：
当虚拟电厂调峰实际出力未达到投标容量的70％时，应受到调峰市场的惩罚Cp,tω如下，
式中，
为t时段虚拟电厂在调峰辅助服务市场的惩罚价格；Δt为调度步长，取1h。
5.根据权利要求4所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述约束条件包括：
(1)虚拟电厂购售电约束，表示虚拟电厂无法同时进行购售电行为：
(2)虚拟电厂外特性：
Pbase,t＝Dt-Pwind,t-Ppv,t (10)；
其中，Pbase,t为虚拟电厂根据风光预测值上报的t时段的运行基线；Pwind,t、Ppv,t分别为t时段的风、光预测值；Dt为柔性负荷调度前t时段的负荷功率；Pout,tω为场景ω下t时段VPP实时运行的外特性功率；Pwinda,tω、Ppva,tω分别为场景ω下t时段的风、光实际出力；D′t为柔性负荷调度后t时段的负荷功率；
分别为EV在场景ω下t时段的充、放电功率；
分别为储能设备在场景ω下t时段的充、放电功率；
(3)虚拟电厂的功率平衡约束：
(4)柔性负荷约束，柔性负荷具有可转移特性，在峰时段削减的负荷总功率与平谷时段增加的负荷总功率相等：
Dt-D′t＝Pla,t-Plc,t (14)；
0≤Pla.t≤Pload,max t∈TV (15)；
0≤Plc.t≤γDt t∈TP (16)；
式中，TV为平谷时段；TP为峰时段；Pla,t为平谷时段中t时段柔性负荷增加的功率，其功率值不应大于最大值Pload,max；Plc,t为在峰时段中t时段柔性负荷削减的负荷总功率；γ是虚拟电厂中柔性负荷的比例；
(5)EV约束：
Eev,min≤Eev,tω≤Eev,max (21)；
式中，Pev,max为EV的最大充、放电功率；Eev,tω为EVt时段的荷电量；
为EV的充、放电效率；Qev为EV车载电池的最大容量；Eev,min、Eev,max为EV车载电池所允许的最小、最大荷电量；
分别为EV在场景ω下t时段的充、放电功率；
(6)储能约束：
Eess,min≤Eess,tω≤Eess,max (26)；
Eess,1ω＝Eess,24ω (27)；
式中，Pess,max为储能的最大充、放电功率；Eess,tω、Eess,(t+1)ω、Eess,1ω、Eess,24ω分别为储能t时段、t+1时段、1时、24时的荷电量；
为储能的充、放电效率；Qess为储能设备的最大容量；Eess,min、Eess,max为储能设备允许的最小、最大荷电量；
分别为储能设备在场景ω下t时段的充、放电功率；
(7)调峰投标约束：
Ppa,minu(t)≤Ppa,tω≤Mu(t) t∈TV (28)；
Ppa,tω＝0 t∈TP (29)；
式中，u(t)为布尔变量，表示该时段虚拟电厂是否参与调峰市场，值为1时表示VPP参与调峰市场，值为0时表示VPP不参与调峰市场；M表示一个很大的正数；Ppa,min为VPP日前向调峰市场投标的最小容量；
(8)调峰实际出力约束：
0≤Ppd,tω≤Pout,tω-Pbase,tω+M(1-u(t)) t∈TV (30)；
Ppd,tω≤u(t)M t∈TV (31)；
Ppd,tω＝0 t∈TP (32)；
(9)CVaR相关约束：
对于一个离散的收益分布，当其置信水平为α时，CVaR近似为小概率1-α情景集合的期望损失；CVaR及其相关约束转化为如下线性问题进行求解：
ξω≥0 (34)；
式中，ξω为引入的辅助变量；δ的最优值为对应的风险价值，即VPP运营收益小于或等于ξ的概率小于或等于1-α的最大收益。
6.根据权利要求5所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述S3包括：
S31，将所述目标函数及所述约束条件共同构成的非线性模型转化为混合整数线性规划问题，即模型的线性化；
S32，通过商业求解器Gurobi求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；其中：
所述S31包括：
对于式(9)、(19)、(24)中连续变量的乘积为0，引入布尔变量；对于式(5)、(8)的分段函数，线性化包括：
0≤Fp,tω≤Ppa,t (35)；
Fp,tω≤Ppd,tω (36)；
Ppd,tω-0.7Ppa,t≥(atω-1)M (37)；
Ppd,tω-0.7Ppa,t≥(1-btω)M (38)；
atω+btω＝1 (39)；
Fp,tω≤atωM (40)；
Ktω≥Ppa,t-Ppd,tω+(btω-1)M (41)；
Ktω≥0 (42)；
式中，atω、btω均为布尔变量。
7.根据权利要求6所述的一种计及风险约束的虚拟电厂同时参与主辅市场的投标方法，其特征在于，所述S4，基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案包括：基于Shapley值法将虚拟电厂参与电力市场获得的总运营收益合理分配，具体步骤如下：对于一个由J个参与者组成的联盟，将成员集合用Ω＝{1,2,…,J}表示，Ω的任意子集称为一个子联盟；利用Shapley值法对虚拟电厂的市场收益进行分配，得到各成员的收益为：
式中，Cj为虚拟电厂内成员j分配的收益；sj为包含成员j的子联盟；|sj|为子联盟sj中的成员数；C(sj)、C(sj/j)分别为包含成员j与不包含成员j的子联盟的运行收益。
8.一种计及风险约束的虚拟电厂同时参与主辅市场的投标系统，所述虚拟电厂含EV，所述主辅市场为主能量市场与调峰辅助服务市场，用于实施如权利要求1-7任一所述的方法，其特征在于，包括：
交易调度流程确定模块，用于确定含EV的虚拟电厂参与电力市场的交易调度流程；
投标模型建立模块，用于基于所述交易调度流程和风险条件价值理论建立所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
模型求解模块，用于求解所述计及风险约束的虚拟电厂同时参与主辅市场的投标模型；
投标模块：用于基于所述投标模型的求解结果将所述虚拟电厂参与所述主辅市场获得的运营收益合理分配，从而获得所述虚拟电厂同时参与主辅市场的投标方案。
9.一种终端设备，其特征在于，包括：一个或多个处理器；存储器，与所述处理器耦接，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如1-7任一所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法。
10.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-7任一所述的计及风险约束的虚拟电厂同时参与主辅市场的投标方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210461034.5A CN114881749A (zh) | 2022-04-28 | 2022-04-28 | 一种计及风险约束的虚拟电厂参与主辅市场的投标方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210461034.5A CN114881749A (zh) | 2022-04-28 | 2022-04-28 | 一种计及风险约束的虚拟电厂参与主辅市场的投标方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114881749A true CN114881749A (zh) | 2022-08-09 |
Family
ID=82672594
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210461034.5A Pending CN114881749A (zh) | 2022-04-28 | 2022-04-28 | 一种计及风险约束的虚拟电厂参与主辅市场的投标方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114881749A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115953011A (zh) * | 2023-03-10 | 2023-04-11 | 中国铁塔股份有限公司 | 通信基站储能资源调度方法及设备 |
| CN117196680A (zh) * | 2023-09-21 | 2023-12-08 | 江苏方天电力技术有限公司 | 一种虚拟电厂参与市场交易的内部主体效益分配方法及装置 |
-
2022
- 2022-04-28 CN CN202210461034.5A patent/CN114881749A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115953011A (zh) * | 2023-03-10 | 2023-04-11 | 中国铁塔股份有限公司 | 通信基站储能资源调度方法及设备 |
| CN117196680A (zh) * | 2023-09-21 | 2023-12-08 | 江苏方天电力技术有限公司 | 一种虚拟电厂参与市场交易的内部主体效益分配方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fasugba et al. | Cost benefits and vehicle-to-grid regulation services of unidirectional charging of electric vehicles | |
| CN114881749A (zh) | 一种计及风险约束的虚拟电厂参与主辅市场的投标方法 | |
| CN109389327B (zh) | 基于风光不确定性的多虚拟电厂时前合作方法 | |
| Zhong et al. | Cooperative operation of battery swapping stations and charging stations with electricity and carbon trading | |
| KR20230103799A (ko) | 가상 발전소 운영 시스템 및 방법 | |
| CN114037192A (zh) | 基于大数据的虚拟电厂交易管理方法、装置、设备及介质 | |
| CN115759556A (zh) | 一种碳-电市场中的多虚拟电厂优化运行与分配方法 | |
| Li et al. | Towards double auction for assisting electric vehicles demand response in smart grid | |
| Wang et al. | Peer-to-peer energy trading for residential prosumers with photovoltaic and battery storage systems | |
| CN115423260A (zh) | 一种电力市场及政策服务新能源利用的量化分析方法 | |
| CN117439080B (zh) | 一种虚拟电厂的调度方法、系统、存储介质及设备 | |
| CN112232716B (zh) | 一种考虑调峰辅助服务的智慧园区优化决策方法 | |
| TWI725606B (zh) | 電動車充電站分散式電能管理方法 | |
| CN117391718A (zh) | 一种基于动态减排因子的绿电-ccer互认交易系统 | |
| Behboodi et al. | Integration of price-driven demand response using plug-in electric vehicles in smart grids | |
| Liu et al. | Market for multi-dimensional flexibility with parametric demand response bidding | |
| Zhang et al. | Shared energy storage market operation mechanism to promote new energy consumption | |
| Savelli et al. | A local market model for community microgrids | |
| CN112365059B (zh) | 一种基于区块链的微网系统设计方法 | |
| Esmaeilian et al. | Techno-economical evaluation of EVs parking lot construction by distribution companies | |
| CN111784540A (zh) | 一种基于多目标粒子群算法的云储能优化出清方法 | |
| CN117081169B (zh) | 一种聚合园区分布式光伏能源的运营方法 | |
| Falabretti et al. | Residential Users as Flexibility Providers: a Techno-Economic Analysis | |
| Jiang et al. | Two-stage optimization model of electric vehicle participating in balancing market | |
| CN109473981B (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 |