基于组合双向拍卖的共享储能机制研究matlab程序

基于组合双向拍卖的共享储能机制研究matlab程序
参考文献：《基于组合双向拍卖的共享储能机制研究》
摘要：为满足共享储能中储能用户的互补性和替代性需求、解决常规单向拍卖中可能存在的垄断竞争问题，提出了基于组合双向拍卖的共享储能机制。在该机制中，储能运营商和储能用户在拍卖商的组织下进行组合双向拍卖。拍卖包含四个阶段：投标、竞胜标决定问题、资源匹配与定价、资金结算。相较于现有拍卖机制，该机制支持“与”投标，使得买家的互补性需求得以准确而简洁地表达。在竞胜标决定问题方面，建立了以社会福利最大化为目标的混合整数非线性规划模型。进一步，将该模型等价转化为混合整数线性规划模型，大幅降低了计算时间复杂度。在定价方面，提出了在买卖双方之间均分社会福利的定价机制，打破了传统单向拍卖的垄断弊端。最后，通过算例验证了所提机制的有效性。关键词：共享储能；组合双向拍卖；互补性；替代性储能是促进可再生能源高效利用、支撑能源互联网建设的关键技术之一。随着可再生能源发电的快速发展，储能在发电侧、用户侧的应用不断增长[4]。当前，储能仍然面临成本偏高、商业模式不成熟的发展瓶颈，而共享思维有可能帮助储能突破此瓶颈。共享储能是指利用不同储能用户储能需求的互补性提升储能设备利用率、进行价值创造的理念。共享的对象包括储能电量和储能容量。一般而言，共享储能的定价机制大致可分为4 类：1）固定价格，即价格与时间、供求关系等无关。文献[8]提出了云储能服务可按容量定价、按流量定价和按套餐定价。固定价格可稳定用户对价格的预期，但无法反映用户对储能效用的评价。2）峰谷电价，即价格与时间或供求关系相关。文献[9]中共享储能的定价策略为储能供大于求时取谷时电价，储能供不应求时取峰时电价。3）利益或成本分摊。此时储能用户实际面对的储能价格取决于利益或成本的分配策略。文献[10]建立了包含共享储能的社区综合能源系统模型并提出了相应的用户协同优化策略，即先最小化总体用能成本再按照能耗量分摊成本。文献[11]建立了基于合作博弈的分布式产消者共享储能模型，并根据 Shapely值分配共享储能收益。4）拍卖，即储能用户或储能卖家进行投标，由其中一方或拍卖商决定中标者与交易价格。文献[12]建立了拥有储能的家庭型用户与共享储能控制器之间基于 Vickery 拍卖的非合作主从博弈模型，其中共享储能的价格由 Vickery拍卖的结果确定。然而，上述共享储能机制均无法反映储能用户对储能的互补性与替代性需求。而在组合拍卖中，买家可对其感兴趣的物品组合投标，可有效解决该问题并降低交易次数与交易成本。互补性与替代性的解释见文献[14]。文献[15]将不同时段的储能视为不同的物品，建立了基于组合单向拍卖的社区级共享储能机制。然而，该机制不支持“与”投标，导致储能用户无法在投标中表达其互补性需求。另一方面，储能卖家对各个时段的投标价格必须相同，导致卖家无法充分表达对不同时段的价格要求。此外，组合单向拍卖仍然属于“一对多”的市场结构，单一卖家具有垄断优势。为克服组合单向拍卖中的垄断优势问题，研究者们提出了组合双向拍卖的思路。组合双向拍卖是组合拍卖和双向拍卖的结合，具有“多对多”的市场结构。文献[18]建立了基于组合双向拍卖的微电网间电力交易的机制，但该机制亦不支持“与”投标，这使得买家无法表达其互补性需求，从而导致组合拍卖的意义大打折扣。

1 基于组合双向拍卖的共享储能问题建模
本文研究的共享储能系统模型架构如图 1 所示。储能卖家和储能买家在拍卖商的组织下对各个时段的储能电量进行组合拍卖。不同时段的储能电量视为不同的物品。一天中有 M 个时段为储能电量拍卖的时段，其余时段为储能装置充电的时段。

1.1 参与者
1.1.1 储能买家
储能买家出于降低自身用能费用等需求购买储能服务。储能买家需要向拍卖商上报对各个时段储能的购买意愿与相应的投标价格。
1.1.2 储能卖家
储能卖家提供储能服务，需要向拍卖商上报对各个时段供给储能的意愿与相应的投标价格。储能卖家可能是独立的储能运营商，也可能是安装有储能的用户等。因此，储能卖家与买家的角色不是固定不变的，本次拍卖中的卖家在下次拍卖中亦可能变为储能买家。
1.1.3 拍卖商
拍卖商负责组织储能卖家与储能买家的整个拍卖过程，主要包括竞标信息的收集、决定中标者、资源匹配与定价、资金结算。
1.2 拍卖流程
首先由拍卖商发布投标开始信号，各储能买家和卖家向拍卖商密封提交投标信息，之后拍卖商发布投标结束信号，投标过程结束；然后拍卖商根据各参与者的竞标信息，以社会福利最大化为目标决定中标者，并对各时段中标的储能资源进行匹配并确定相应的交易价格，进而计算各买家的总支出与各卖家的总收入；最后拍卖商将中标结果、资源匹配结果、定价结果、参与者的收支情况等信息通知各参与者；拍卖结果以负荷曲线(储能买家的负荷为正值，储能卖家的负荷为负值)的形式上报给微电网控制中心或配电网调度中心，而后由微电网控制中心下达控制指令或配电网调度中心下达调度指令执行拍卖结果。拍卖结果执行后，买家将应付款转给拍卖商，之后拍卖商将应收款转给卖家，拍卖过程结束。

2 投标与竞胜标决定机制
2.1 投标
2.1.1 卖家投标信息

2.1.2 买家投标信息

2.2 竞胜标决定问题
2.2.1 目标函数
买卖双方投标完成之后，拍卖商需要确定中标者。在组合双向拍卖中，“多对多”的市场结构破除了单向拍卖中一方存在的垄断优势。在单向拍卖中往往以垄断一方的福利最大化为目标决定获胜者，而此处拍卖商以社会福利最大化为目标确定获胜者。社会福利Wso 为卖家剩余与买家剩余之和，亦等于买家出价与卖家要价之差：

2.2.2 约束条件

3 算例分析
假设一天中共有 6 个时段的储能可拍卖，不同时段的储能视为不同物品。算例包含 7 个储能买家和 6 个储能卖家。惩罚系数 的取值为 1×106。本场景中各个储能买家和储能卖家的投标策略是固定的且不考虑输配电价。储能买家的投标数据如表 1 所示，其中，(10,1.2)表示买家对该时段储能的投标电量为 10 MW·h、投标价格为 1.2 元/(kW·h)。买家 1 与买家 4 对多个时段的储能进行“与”投标，表达了其互补性储能需求。买家 2 与买家 6 对多个时段的储能进行“异或”投标，表达了其替代性储能需求。而买家 3 与买家 5 为多多益善型买家，其投标为“或”投标。买家 7 仅对第 6 时段的储能进行了原子投标。​

卖家的投标数据如表 2 所示，其中，(5, 0.9)的含义为储能卖家对该时段储能的投标电量为5 MW·h、投标价格为 0.9 元/(kW·h)。储能卖家可对若干时段进行投标，并且各个时段的投标电量和投标价格不必相同，使得各卖家可以充分表达其对不同时段的储能的价格要求。

4 算例结果程序运行结果
1）买家中标

2）卖家中标程序运行结果

文献结果
1）买家中标

2）卖家中标

4 matlab程序
clc
clear
close all
%% 买家信息
%0，原子；买家只对一个时段的储能投标
% 1或；不可超过投标时段的个数
% 2与；对各个时段的投标同时中标或同时不中标
% 3异或；若干时段储能的投标最多有 1 个时段中
%类型+容量
S_b1=[
2, 10, 10, 0, 0, 0, 0;
3, 0, 0, 5, 3, 0, 0;
1, 0, 0, 0, 0, 3, 3;
2, 4, 4, 7, 0, 0, 0;
1, 0, 0, 7, 9, 5, 6;
3, 2, 0, 0, 2, 0, 2;
0, 0, 0, 0, 0, 0, 8;
];
%序号+价格
S_b2=[
1, 1.2, 1.2, 0, 0, 0, 0;
2, 0, 0, 1.2, 1.1, 0, 0;
3, 0, 0, 0, 0, 0.9, 0.9;
4, 1, 1, 1, 0, 0, 0;
5, 0, 0, 0.5, 1.1, 1.2, 0.5;
6, 0.7, 0, 0, 0.6, 0, 0.4;
7, 0, 0, 0, 0, 0, 1;
];
%% 卖家信息
%序号+容量
S_s1=[
1, 5, 10, 10, 0, 0, 0;
2, 0, 0, 5, 5, 2, 5;
3, 2, 0, 0, 0, 4, 3;
4, 5, 5, 0, 0, 0, 0;
5, 0, 8, 8, 8, 0, 5;
6, 5, 0, 0, 8, 4, 4;
];
%序号+价格
S_s2=[
1, 0.9, 0.9, 0.9, 0, 0, 0;
2, 0, 0, 0.8, 0.7, 0.5, 0.5;
3, 0.4, 0, 0, 0, 0.6, 0.7;
4, 0.6, 0.6, 0, 0, 0, 0;
5, 0, 0.7, 0.7, 0.7, 0, 0.5;
6, 0.5, 0, 0, 1, 0.4, 0.4;
];

。。。。。。。略