上海交通大学学报, 2024, 58(12): 1868-1880 doi: 10.16183/j.cnki.jsjtu.2023.149

新型电力系统与综合能源

计及风电不确定性及碳交易的储能参与电力市场竞价策略

张宪文1, 殷高文2, 沈非凡,2, 黄晟2, 魏娟2

1.临朐县重大项目服务中心,山东 潍坊 262600

2.湖南大学 电气与信息工程学院,长沙 410082

Bidding Strategies for Energy Storage Participation in Electricity Market Considering Uncertainty of Wind Power and Carbon Trading

ZHANG Xianwen1, YIN Gaowen2, SHEN Feifan,2, HUANG Sheng2, WEI Juan2

1. Linqu County Major Project Service Center, Weifang 262600, Shandong, China

2. College of Electrical and Information Engineering, Hunan University, Changsha 410082, China

通讯作者: 沈非凡,助理教授;E-mail:shenfeifan@hnu.edu.cn.

责任编辑: 王历历

收稿日期: 2023-04-20   修回日期: 2023-06-29   接受日期: 2023-06-30  

基金资助: 湖南省制造业关键产品“揭榜挂帅”项目(2022GXGG019)
湖南省重点研发计划项目(2021GK2019)
湖南省科技创新计划资助(2021RC2046)

Received: 2023-04-20   Revised: 2023-06-29   Accepted: 2023-06-30  

作者简介 About authors

张宪文(1970—),高级工程师,现主要从事电力系统运行优化和电力市场研究.

摘要

电源侧储能是促进新能源消纳、践行“双碳”战略的有力抓手.针对电源侧储能作为新兴主体参与竞争性联合电力市场竞价机制尚不完善的问题,从电源侧储能作为价格制定者参与电力现货市场的博弈关系出发,考虑风电出力的不确定性以及碳交易引入下电力市场的耦合性,提出一种电源侧储能参与电能量及备用市场多场景应用的竞标博弈双层优化模型.上层电源侧储能作为领导者以自身综合收益最大为目标进行投标决策,下层电力交易中心作为追随者以社会购电成本及平衡调节成本最小完成市场出清.利用MATLAB仿真验证了所提模型在优化电源侧储能参与电力市场多场景应用竞标、提升储能综合运行收益方面的合理性和有效性.

关键词: 电源侧储能; 电力市场; 风电不确定性; 碳交易; 竞价策略; 双层优化模型

Abstract

Power side energy storage is a powerful tool for promoting the consumption of new energy and implementing the “dual carbon” strategy. In view of the imperfection of the bidding mechanism of the power side energy storage as an emerging entity participating in the competitive joint electricity market, and starting from the game relationship of the power side energy storage as a price maker participating in the spot electricity market, a bidding game two-level optimization model of the power side energy storage participating in the multi scenario application of the electric energy and reserve ancillary service market is proposed, considering the uncertainty of wind power output and the coupling of the electricity market with the introduction of carbon trading. The energy storage at the upper power supply side, as a leader, makes bidding decisions with the goal of maximizing its own comprehensive income, while the lower power trading center as a follower completes market clearing with the minimum cost of social power purchase and balance regulation. The MATLAB simulation verifies the rationality and effectiveness of the proposed model in optimizing the participation of power side energy storage in multi scenario application bidding in the electricity market and improving the comprehensive operating income of the energy storage.

Keywords: power side energy storage; electricity market; uncertainty of wind power; carbon trading; bidding strategy; two-level optimization model

PDF (2880KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张宪文, 殷高文, 沈非凡, 黄晟, 魏娟. 计及风电不确定性及碳交易的储能参与电力市场竞价策略[J]. 上海交通大学学报, 2024, 58(12): 1868-1880 doi:10.16183/j.cnki.jsjtu.2023.149

ZHANG Xianwen, YIN Gaowen, SHEN Feifan, HUANG Sheng, WEI Juan. Bidding Strategies for Energy Storage Participation in Electricity Market Considering Uncertainty of Wind Power and Carbon Trading[J]. Journal of Shanghai Jiaotong University, 2024, 58(12): 1868-1880 doi:10.16183/j.cnki.jsjtu.2023.149

随着新型电力系统的加紧构建,新能源的反调峰、波动性等不利特性日益凸显,威胁电力系统稳定可靠运行.储能因能暂存电能,摆脱了传统发电机组即发即用的束缚,可进行低储高发,有效减少弃风、弃光现象,同时平滑新能源出力,利于电力系统安全[1-2].因此,大力发展电源侧储能产业已成为全球共识[3-4].

电源侧储能产业的健康发展依赖于合理的市场政策与运营指导机制.然而,我国电力市场建设起步较晚,尚未形成完善的针对电源侧储能参与电力市场竞价的整体运行机制[5-6].总体而言,现阶段电源侧储能参与电力市场运行调度策略方面主要存在以下问题:一是由于缺乏成熟的商业模式,电源侧储能多参与单一应用场景运行,利用率低,营收情况不佳.目前我国电源侧储能多为新能源发电厂商兴建,主要作为辅助设施优化新能源出力[7].仅参与单一应用场景使得电源侧储能非工作时段功率及容量大量闲置,未充分发挥应用潜力,运营经济性差,成本回收期长.二是缺乏具体有效的指导电源侧储能同其他市场主体竞标的策略机制.与传统发电机组仅作为市场中的发电方不同,储能在新能源电力超发阶段还可充当用能角色进行充电,且储存的电力又可参与其他应用场景[8-9].如何考虑市场动态且兼顾储能特性进行最优的市场竞标是储能良好参与电力市场的关键.三是碳交易引入下传统机组在电能量市场的报价与其碳交易成本存在耦合性[10],影响传统机组的投标策略,进而增加电力系统中能量分配的难度.因此,研究电源侧尤其是新能源侧储能在电力市场环境下的经济调度策略对于其优化运营、提高营收、降低扶持政策依赖,并促进储能参与下的电力市场健康有序发展具有重要意义.

针对储能多场景运行的投标策略研究中,文献[11-12]中提出一种风储系统参与电能量市场和调频市场的报价策略,可指导风储系统在联合市场中竞价投标.文献[13-14]中综合考虑储能调频可靠性与经济性,研究储能参与电能量和一次调频联合市场的竞价与出清机制.但是,上述针对联合市场的竞价策略研究虽然考虑了多场景应用,但在市场类型上又均只考虑日前市场,未考虑与日前市场高度耦合的实时平衡市场.缺乏对实时平衡市场不确定性考量,会使得日前决策在应对电力市场实际交割时新能源出力不确定性方面存在局限性,增加储能的机会损失.

目前,针对储能参与电力现货市场日前-实时两阶段、多场景应用的竞标指导策略还未完全形成,但已有相关研究开展.文献[15]中提出一种基于风险评估的储能参与日前电能量市场、旋转备用市场以及日内平衡市场的竞标策略,计及市场价格不确定性以及备用需求不确定性,进行报价决策优化.文献[16-17]中研究储能在日前电能量和辅助服务联合市场中的运行情况,并分别采用概率分布、多场景生成等方法模拟实时市场的不确定性影响,优化竞标决策方案.

进一步地,电力市场作为一个规模庞大、各市场主体关系耦合的复杂交易体系,每个个体的决策都会对市场出清产生影响.随着我国引入碳排放权交易,出于发电成本的考量,火电机组在电力现货市场中的售电报价将直接受到影响,因而引起电能量市场出清结果变动.文献[18-20]中均基于碳排放和运行经济性考虑,对能源系统进行优化调度研究.但上述研究虽然对碳交易下能源调度做了研究,但未同时计及新能源出力的不确定性.

综上所述,针对缺乏有效的能同时计及新能源出力不确定性及电-碳市场耦合性的电源侧储能参与电力市场多场景运行调度的指导策略这一问题,以风电为例,开展电源侧储能多应用场景协同运行的竞价投标策略研究,细化储能与电力交易中心之间的博弈关系,充分考虑风电出力不确定性和碳排放权交易对储能投标策略的影响,借鉴双层模型优化思路[21-23],构建电源侧储能参与日前电能量市场、日前备用辅助服务市场以及实时平衡市场竞价投标双层优化模型.借助Karush-Kuhn-Tucker(KKT)条件[24]将上述双层模型转化为单层模型后,进一步对其线性化.最后,基于MATLAB仿真平台求解模型,仿真结果验证了所提策略在市场出清过程中指导储能优化自身决策和提升收益的有效性.

1 电源侧储能参与电力市场主从博弈竞标机制分析

1.1 主从博弈竞标机制

主从博弈是基于博弈论的一种动态博弈,参与博弈的各方在决策时存在先后顺序,即领导与跟随的主从关系[25-26].先做出决策的主体称为领导者,根据领导者的决策情况后做出决策的主体称为跟随者.主从博弈基本原理如图1所示.

图1

图1   主从博弈基本原理

Fig.1   Basic principles of master-slave game


在电力市场中,电力交易中心需依托储能等各竞标方的投标数据,组织电力市场有序运行,降低电力交易成本,而关系储能收益的中标量与出清电价又受到电力交易中心的决策影响.因此,考虑到电力市场运行过程中储能与电力交易中心之间存在利益制衡,且各自追求收益最大化的特点与博弈论研究对象特点类似,故可将电源侧储能参与日前及日内电力市场竞标及市场出清过程建模为主从博弈模型.

从电源侧储能角度出发,假设储能具有优先决策权,电力交易中心跟随决策,则电源侧储能作为主从博弈的领导者,在考虑自身运行约束的基础上充分考虑市场对其自身投标计划的反应,提前对可能出现的电力市场出清情况进行预测,优化决策出最佳竞价投标方案.

1.2 考虑碳交易的电力市场机制

碳排放权交易是我国为减少温室气体排放所进行的碳排放配额(carbon emission allowance, CEA)交易机制,主要面向高碳排放的重点单位,如燃煤电厂等,交易对象为CEA.重点单位超排即配额缺额部分可向碳交易市场购买,满足自身排放需求,反之盈余部分也可向碳排放权交易市场出售进行获利.

传统机组在碳排放权交易市场购买缺额后,相应单位发电成本增加,因而会选择在电力现货市场提高报价以补贴发电成本.若传统机组碳排放量未超出配额,可出售配额盈余,则单位发电成本降低,传统机组将在电力现货市场降低报价以提高竞标竞争力[27].因此,可根据传统机组在碳交易市场和电力现货市场中的成本转移行为构建传统机组的电-碳耦合报价机制,其形成机理如图2所示.

图2

图2   考虑碳排放权交易的传统机组电-碳耦合报价机制

Fig.2   Electricity-carbon coupling pricing mechanism for traditional units considering carbon trading


传统机组在电力市场中的电-碳耦合报价为

γg,tec,o= γeg,t+ γcg
γcg=(Crg- Cbg)Vc

式中:γg,tec,o为传统机组g在电力市场t时刻的电-碳耦合报价;γeg,t为传统机组g在电力市场t时刻基本的电报价;γcg为传统机组g的碳交易成本;Crg为传统机组g的碳排放量实际值;Cbg为传统机组g的碳排放量配额值;Vc为碳排放权交易市场中的碳单价.

2 电源侧储能电力市场多场景应用竞标博弈双层优化模型

鉴于电源侧储能竞标报价与市场出清电价之间具有相互影响的耦合关系,电源侧储能参与电力市场竞标过程可具体建模为一个决策变量能在上下层间传递的双层模型,以凸显主从博弈过程.上层为电源侧储能竞标决策层,电源侧储能根据自身功率、报价等约束,追求日前与实时平衡市场综合营收最大,规划自己的量-价投标方案并报至下层电力交易中心;下层为电力市场出清层,电力交易中心统筹电源侧储能、传统机组等各交易主体上报的竞标数据,综合考虑新能源预测出力偏差带来的净负荷偏差,以及传统机组碳交易下的电-碳耦合报价,进行日前-实时平衡市场的出清,最后把电源侧储能在各市场中的中标量和市场出清电价返回给电源侧储能决策层,进而储能可优化初始投标策略,得到更优竞标方案.双层优化模型构建如图3所示.

图3

图3   电源侧储能电力市场多场景应用竞标博弈双层优化模型

Fig.3   Two-level optimization model for participation of power side energy storage in multi-scenario application bidding games in an electricity market


2.1 上层电源侧储能竞标决策模型

(1) 目标函数.储能收益来源包括日前市场的电能量市场、备用辅助服务市场以及实时平衡市场.每个收益来源都基于储能的两种工作模式,即充电和放电.上层模型以储能在3种市场中的综合收益最大为目标建立目标函数

$ \begin{array}{c} \max \sum_{t=1}^{T} \sum_{e=1}^{E}\left\{\left[\lambda_{t}^{\mathrm{e}}\left(P_{e, t}^{\mathrm{dis}}-P_{e, t}^{\mathrm{ch}}\right)\right]+\right. \\ {\left[\lambda_{t}^{\mathrm{up}}\left(S_{e, t}^{\mathrm{ch}, \text { up }}+S_{e, t}^{\text {dis, up }}\right)+\right.} \\ \left.\lambda_{t}^{\text {down }}\left(S_{e, t}^{\mathrm{ch}, \text { down }}+S_{e, t}^{\text {dis.down }}\right)\right]+ \\ \sum_{k=1}^{K} \phi_{k}\left[\lambda_{t, k}^{\mathrm{b}}\left(P_{e, t, k}^{\mathrm{ch}, \text { up }}-P_{e, t, k}^{\mathrm{ch}, \text { down }}\right)+\right. \\ \left.\left.\lambda_{t, k}^{\mathrm{b}}\left(P_{e, t, k}^{\mathrm{dis}, \mathrm{pp}}-P_{e, t, k}^{\text {dis,down }}\right)\right]\right\} \end{array}$

式中:TE分别为储能调度周期和储能数量;Kϕk分别为实时平衡市场中净负荷偏差不确定性场景k的数量和对应场景出现概率;λett时刻电能量市场出清价格;λtupλtdown分别为t时刻备用辅助服务市场上调备用和下调备用的市场出清价格;λbt,k为场景kt时刻实时平衡市场中容量调用出清价格;Pe,tchPe,tdis分别为储能e在电能量市场t时刻中标的充、放电功率;Se,tch,upSe,tch,down分别为储能e在备用辅助服务市场t时刻中标的充电状态下所提供的上、下调备用容量;Se,tdis,upSe,tdis,down分别为储能e在备用辅助服务市场t时刻中标的放电状态下所提供的上、下调备用容量;Pe,t,kch,upPe,t,kch,down分别为储能e在场景k的实时平衡市场t时刻中标的充电状态下所提供的上、下调容量;Pe,t,kdis,upPe,t,kdis,down分别为储能e在场景k的实时平衡市场t时刻中标的放电状态下所提供的上、下调容量.

(2) 约束条件.充、放电状态约束为

Ue,tch+ Ue,tdis≤1

式中:Ue,tchUe,tdis分别为储能e充、放电动作状态0-1变量.约束式(4)的目的在于禁止储能e同时进行充、放电.

功率约束为

0≤ Pe,tch,oUe,tchPerated
0≤ Se,tch,up,oUe,tchSeup,max
0≤ Se,tch,down,oUe,tchSedown,max
Se,tch,down,o+ Pe,tch,oUe,tchPerated
Se,tch,up,o- Pe,tch,o≤0
0≤ Pe,tdis,oUe,tdisPerated
0≤ Se,tdis,up,oUe,tdisSeup,max
0≤ Se,tdis,down,oUe,tdisSedown,max
Se,tdis,up,o+ Pe,tdis,oUe,tdisPerated
Se,tdis,down,o- Pe,tdis,o≤0

式中:Pe,tch,oPe,tdis,o分别为储能e在电能量市场t时刻申报的充、放电功率;Perated为储能e的额定功率;Se,tch,up,oSe,tch,down,o分别为储能e在备用辅助服务市场t时刻申报的充电状态下所能提供的上、下调备用容量;Se,tdis,up,oSe,tdis,down,o分别为储能e在备用辅助服务市场t时刻申报的放电状态下所能提供的上、下调备用容量;Seup,maxSedown,max分别为储能e的最大上、下调节能力.约束式(5)和(10)限制储能申报的充、放电功率不超过额定值.式(6)和(7)限制储能申报的其充电模式下的上、下调备用容量不超过额定调节能力.式(8)表示储能所申报的在t时刻以充电方式提供下调备用的容量与所申报的t时刻的充电功率之和应不超过储能额定充电功率.式(9)表示储能申报的在t时刻以减小充电功率方式提供的上调备用容量不应大于此时储能申报的充电功率.式(11)和(12)限制储能申报的放电模式下的上、下调备用容量不超过其最大调节能力.式(13)表示储能所申报的在t时刻以放电模式提供上调备用的容量与申报的同一时刻的放电功率之和应不超过储能额定放电功率.式(14)表示储能申报的在t时刻以减小放电功率方式提供的下调备用容量不应大于此时储能申报的放电功率.

能量约束为

0≤Ee,tEerated
$ \begin{aligned} E_{e, t} & =E_{e}^{\mathrm{init}}-\left[P_{e, t}^{\mathrm{dis}}+\right. \\ & \left.\sum_{k=1}^{K} \phi_{k}\left(P_{e, t, k}^{\mathrm{dis}, \mathrm{pp}}-P_{e, t, k}^{\mathrm{dis}, \mathrm{down}}\right)\right] / \eta_{\mathrm{dis}} \Delta t+ \\ & {\left[P_{e, t}^{\mathrm{ch}}+\sum_{k=1}^{K} \phi_{k}\left(-P_{e, t, k}^{\mathrm{ch}, \mathrm{up}}+P_{e, t, k}^{\mathrm{ch}, \mathrm{down}}\right)\right] \eta_{\mathrm{ch}} \Delta t } \end{aligned}$
Eeinit=Ee,24

式中:Ee,t为储能et时刻的能量值;EeratedEeinit分别为储能e的额定能量值和在调度周期初始时刻能量值;Ee,24为调度周期末储能e的能量值;ηchηdis分别为储能e的充、放电效率.约束式(15)表示储能的能量值应在额定范围内.式(16)为储能在任一时刻的能量值.式(17)表示为保证储能调度的循环进行,应当令每个调度周期末尾的能量值回到初始值.

申报价格约束为

0≤ γe,te,oγemax
0≤ γe,tch,up,o, γe,tch,down,o, γe,tdis,up,o, γe,tdis,down,oγemax

式中:γe,te,o为储能e在电能量市场t时刻申报的充放电价格;γe,tch,up,oγe,tch,down,o分别为储能e在备用辅助服务市场t时刻申报的充电状态下的上、下调备用容量价格;γe,tdis,up,oγe,tdis,down,o分别为储能e在备用辅助服务市场t时刻申报的放电状态下的上、下调备用容量价格;γemax为储能e在电力市场的报价上限.约束式(18)和(19)用以限制储能报价非负且在限定范围之内.

2.2 下层电力市场出清模型

2.2.1 日前电力市场出清模型

(1) 目标函数

$ \begin{array}{l} \min \left[\sum_{e=1}^{E}\left(\gamma_{e, t}^{\mathrm{c}, 0} P_{e, t}^{\mathrm{dis}}-\gamma_{e, t}^{\mathrm{c}, 0} P_{e, t}^{\mathrm{ch}}\right)+\sum_{g=1}^{G} \gamma_{k, t}^{\mathrm{ce}, P_{k}} P_{g, t}+\right. \\ \sum_{e=1}^{E}\left(\gamma_{e, t}^{\text {ch, up }, 0} S_{e, t}^{\text {ch, up }}+\gamma_{e, t}^{\text {ch, down }, o} S_{e, t}^{\text {ch,down }}+\right. \\ \left.\gamma_{e, t}^{\text {dis, up }, o} S_{e, t}^{\text {dis. up }}+\gamma_{e, t}^{\text {dis.down } o \mathrm{o}} S_{e, t}^{\text {dis.down }}\right)+ \\ \left.\sum_{g=1}^{G} \gamma_{k, t}^{s, 0}\left(S_{g, t}^{\mathrm{up}}+S_{g, t}^{\text {down }}\right)\right] \end{array}$

式中:G为参与电力市场的传统机组数量;γg,ts,o为传统机组g在日前备用辅助服务市场t时刻的报价;Pg,t为传统机组g在电能量市场t时刻的中标电量;Sg,tupSg,tdown分别为传统机组g在备用辅助服务市场t时刻中标的上调备用容量和下调备用容量.目标函数式(20)以市场购电成本最小为目标构建目标函数.考虑到对风电的保障性消纳,将风电作为负的负荷对待,其不参与市场报价,只在约束中对其考量.

(2) 约束条件.潮流平衡约束为

e=1EPe,tdis-Pe,tch+ w=1WPw,t+ g=1GPg,t= l=1LPl,t: λet

式中:W为风电场数量;Pw,t为风电场wt时刻的出力;L为负荷数量;Pl,t为负荷lt时刻的功率需求值.

备用平衡约束为

e=1E(Se,tch,up+ Se,tdis,up)+ g=1GSg,tup= Stup: λtup
e=1E(Se,tch,down+ Se,tdis,down)+ g=1GSg,tdown= Stdown: λtdown

式中:StupStdown分别为t时刻市场上、下调备用容量总需求.

电能量市场中标功率约束为

0≤Pg,tPog,t: μ_g,t, μ-g,t
0≤ Pe,tchPe,tch,o: μ_e,tch, μ-e,tch
0≤ Pe,tdisPe,tdis,o: μ_e,tdis, μ-e,tdis

式中:Pog,t为传统机组g在电能量市场t时刻的申报功率.约束式(24)~(26)用于限制各市场主体所中标的功率不超过其申报值.

备用辅助服务市场中标备用容量约束为

0≤ Sg,tupSg,tup,o: μ_g,tup, μ-g,tup
0≤ Sg,tdownSg,tdown,o: μ_g,tdown, μ-g,tdown
0≤ Se,tch,upSe,tch,up,o: μ_e,tch,up, μ-e,tch,up
0≤ Se,tch,downSe,tch,down,o: μ_e,tch,down, μ-e,tch,down
0≤ Se,tdis,upSe,tdis,up,o: μ_e,tdis,up, μ-e,tdis,up
0≤ Se,tdis,downSe,tdis,down,o: μ_e,tdis,down, μ-e,tdis,down

式中:Sg,tup,oSg,tdown,o分别为传统机组g在备用辅助服务市场t时刻的上、下调备用容量申报量.约束式(27)~(32)用于限制各市场主体中标量不超过其申报量.

电能量与备用容量耦合约束为

Sg,tup+Pg,tPog,t: μg,tup
Sg,tdown-Pg,t≤0: μg,tdown
Se,tdis,downPe,tdis: μe,tdis,down
Se,tch,upPe,tch: μe,tch,up

约束式(33)~(36)用于限制传统机组和储能各自的中标能量及其中标容量之间的耦合关系.

2.2.2 日内实时平衡市场出清模型

(1) 目标函数

$ \begin{array}{l} \min \left[\sum_{g=1}^{G} \gamma_{g, t}^{\mathrm{cc}, \mathrm{o}}\left(P_{g, t, k}^{\mathrm{up}}-P_{g, t, k}^{\mathrm{down}}\right)+\right. \\ \sum_{e=1}^{E} \gamma_{e, t}^{\mathrm{c}, \mathrm{o}}\left(P_{e, t, k}^{\mathrm{dis}, \mathrm{pp}}-P_{e, t, k}^{\mathrm{dis} \text { down }}\right)+ \\ \left.\sum_{e=1}^{E} \gamma_{e, t}^{\mathrm{c}, \mathrm{o}}\left(P_{e, t, k}^{\mathrm{ch}, \mathrm{up}}-P_{e, t, k}^{\mathrm{ch}, \mathrm{down}}\right)+\sum_{w=1}^{\mathrm{d}} V_{w, t} l_{w, t, k}\right] \end{array}$

式中:Pg,t,kupPg,t,kdown分别为传统机组g在场景kt时刻上、下调节的中标量;Vw,t为风电场wt时刻单位弃风量损失;lw,t,k为因电力市场对净负荷偏差调节不足导致的风电场wt时刻的弃风量.目标函数式(37)以电力市场对净负荷偏差的不平衡调节成本最小为目标构建目标函数.

(2) 约束条件.净负荷偏差等式约束为

$ \begin{array}{c} \sum_{e=1}^{E}\left(P_{e, t, k}^{\mathrm{ch}, \mathrm{up}}-P_{e, t, k}^{\mathrm{ch}, \mathrm{down}}\right)+\sum_{e=1}^{E}\left(P_{e, t, k}^{\mathrm{dis}, \text { up }}-P_{e, t, k}^{\mathrm{dis}, \text { down }}\right)+ \\ \sum_{g=1}^{G}\left(P_{g, t, k}^{\mathrm{up}}-P_{g, t, k}^{\mathrm{down}}\right)+l_{w, t, k}=P_{t, k}: \lambda_{t, k}^{\mathrm{b}} \end{array}$

式中:Pt,k为系统在场景kt时刻的净负荷偏差.

传统机组中标功率约束为

0≤ Pg,t,kupSg,tup: ρ_g,t,kup, ρ-g,t,kup
0≤ Pg,t,kdownSg,tdown: ρ_g,t,kdown, ρ-g,t,kdown

约束式(39)和(40)表示为保证实时平衡市场调节容量的可调用,传统机组在实时平衡市场的中标量应小于其在该时刻中标的备用容量.

电源侧储能中标功率约束为

0≤ Pe,t,kch,upSe,tch,up: ρ_e,t,kch,up, ρ-e,t,kch,up
0≤ Pe,t,kch,downSe,tch,down: ρ_e,t,kch,down, ρ-e,t,kch,down
0≤ Pe,t,kdis,upSe,tdis,up: ρ_e,t,kdis,up, ρ-e,t,kdis,up
0≤ Pe,t,kdis,downSe,tdis,down: ρ_e,t,kdis,down, ρ-e,t,kdis,down

约束式(41)~(44)表示为保证实时平衡市场调节容量的可调用,储能的中标量应小于其在该时刻中标的备用容量.

弃风量约束为

0≤lw,t,k≤Pw,t: ρ_w,t,k, ρ-w,t,k

约束式(45)旨在限制风电场wt时刻的弃风量不超过其在该时刻的计划出力.

上述讨论的约束式(21)~(36)、(38)~(45)中各式冒号后的变量均为该约束式所对应的对偶变量.其中,对偶变量μρ无实际物理意义.

3 双层优化模型求解

所提电源侧储能电力市场多场景应用竞标博弈双层优化模型为混合整数非线性(mixed integer nonlinear programming, MINLP)问题,直接进行求解具有难度.因此,采用KKT条件先将双层模型转化为单层模型,即变为带均衡约束的数学规划(mathematical programs with equilibrium constraints, MPEC)模型.但因该模型中含有双变量乘积项,属于非线性问题,故再进一步利用大M[28]和二元展开法[29]将MPEC模型作线性处理,最终变换成混合整数线性规划(mixed integer linear programming, MILP)模型,此时即可借助Yamlip+Gurobi专业求解器在MATLAB平台上完成模型求解.所提竞标博弈双层优化模型整体求解流程如图4所示.

图4

图4   竞标博弈双层优化模型求解流程

Fig.4   Solution process of biding game two-level optimization model


以下将对KKT条件转化及线性化处理过程作进一步说明,因下层中实时平衡市场出清模型的处理过程与日前市场出清模型处理过程类似,所以仅以日前电力市场出清模型为例进行说明.

3.1 下层模型的KKT条件

(1) KKT等式约束条件

γe,te,o- λet+ μ-e,tdis- μ_e,tdis- μe,tdis,down=0
- γe,te,o+ λet+ μ-e,tch- μ_e,tch- μe,tch,up=0
γg,tec,o- λet+ μ-g,t- μ_g,t+ μg,tup- μg,tdown=0
γe,tdis,up,o- λtup+ μ-e,tdis,up- μ_e,tdis,up=0
γe,tdis,down,o- λtdown+ μ-e,tdis,down- μ_e,tdis,down+ μe,tdis,down=0
γe,tch,up,o- λtup+ μ-e,tch,up- μ_e,tch,up+ μe,tch,up=0
γe,tch,down,o- λtdown+ μ-e,tch,down- μ_e,tch,down=0
γg,ts,o- λtup+ μ-g,tup- μ_g,tup+ μg,tup=0
γg,ts,o- λtdown+ μ-g,tdown- μ_g,tdown+ μg,tdown=0

式(46)~(54)分别为对原模型中决策变量Pe,tdisPe,tchPg,tSe,tdis,upSe,tdis,downSe,tch,upSe,tch,downSg,tupSg,tdown求偏导形成的KKT等式约束条件.

(2) KKT互补松弛条件

0≤Pg,tμ_g,t≥0
0≤(Pog,t-Pg,t)⊥ μ-g,t≥0
0≤ Sg,tupμ_g,tup≤0
0≤(Sg,tup,o- Sg,tup)⊥ μ-g,tup≥0
0≤ Sg,tdownμ_g,tdown≥0
0≤(Sg,tdown,o- Sg,tdown)⊥ μ-g,tdown≥0
0≤(Pog,t-Pg,t- Sg,tup)⊥ μg,tup≥0
0≤(Pg,t- Sg,tdown)⊥ μg,tdown≥0
0≤ Pe,tchμ_e,tch≥0
0≤(Pe,tch,o- Pe,tch)⊥ μ-e,tch≥0
0≤ Pe,tdisμ_e,tdis≥0
0≤(Pe,tdis,o- Pe,tdis)⊥ μ-e,tdis≥0
0≤ Se,tch,upμ_e,tch,up≥0
0≤(Se,tch,up,o- Se,tch,up)⊥ μ-e,tch,up≥0
0≤ Se,tch,downμ_e,tch,down≥0
0≤(Se,tch,down,o- Se,tch,down)⊥ μ-e,tch,down≥0
0≤ Se,tdis,upμ_e,tdis,up≥0
0≤(Se,tdis,up,o- Se,tdis,up)⊥ μ-e,tdis,up≥0
0≤ Se,tdis,downμ_e,tdis,down≥0
0≤(Se,tdis,down,o- Se,tdis,down)⊥ μ-e,tdis,down≥0
0≤(Pe,tch- Se,tch,up)⊥ μe,tch,up≥0
0≤(Pe,tdis- Se,tdis,down)⊥ μe,tdis,down≥0

式(55)~(76)为利用KKT理论推导的式(24)~(36)的互补松弛条件.其中,“⊥”表示其左右侧两个非负值至多有一个可以严格大于0.

经过上述KKT条件数学变换,下层模型已转化为上层模型的约束形式,即原双层模型已变换为单层模型.

3.2 模型线性化

(1) 互补松弛条件的线性化.因互补松弛约束式(55)~(76)中存在决策变量与对偶变量的乘积项,故采用大M法对其进行精确化线性处理,在此仅以对式(55)线性化为例进行说明,结果如式下:

0≤Pg,tMu
0≤ μ_g,t≤M(1-u)

式中:M为一个足够大的正数;u为0-1二进制辅助变量.

(2) 目标函数的线性化.因单层化后的MPEC模型目标函数式(3)中存在大量双变量乘积的非线性项,所以采用强对偶定理对其进行处理

$ \begin{array}{l} \max \sum_{t=1}^{T} \sum_{e=1}^{E}\left\{\lambda_{t}^{\mathrm{e}}\left(P_{e, t}^{\mathrm{dis}}-P_{e, t}^{\mathrm{ch}}\right)+\right. \\ {\left[\lambda_{t}^{\mathrm{up}}\left(S_{e, t}^{\mathrm{ch}, \mathrm{up}}+S_{e, t}^{\mathrm{dis}, \text { up }}\right)+\right.} \\ \left.\lambda_{t}^{\text {down }}\left(S_{e, t}^{\text {ch,down }}+S_{e, t}^{\text {dis.down }}\right)\right]+ \\ \sum_{k=1}^{K} \phi_{k}\left[\lambda_{t, k}^{\mathrm{b}}\left(P_{e, t, k}^{\mathrm{ch}, \text { up }}-P_{e, t, k}^{\mathrm{ch}, \text { down }}\right)+\right. \\ \left.\left.\lambda_{t, k}^{\mathrm{b}}\left(P_{e, t, k}^{\mathrm{dis}, \text { up }}-P_{e, t, k}^{\text {dis, down }}\right)\right]\right\}= \\ \max \left[-\sum_{g=1}^{G} \gamma_{g, t}^{\mathrm{cc}, \mathrm{o}} P_{g, t}-\sum_{g=1}^{G} \gamma_{g, t}^{s, 0}\left(S_{g, t}^{\mathrm{up}}+S_{g, t}^{\mathrm{down}}\right)+\right. \\ \sum_{g=1}^{G} \bar{\mu}_{g, t} P_{g, t}^{\mathrm{o}}-\sum_{\mathrm{g}=1}^{G} \bar{\mu}_{g, t}^{\mathrm{up}} S_{g, t}^{\mathrm{up}, \mathrm{o}}- \\ \sum_{g=1}^{G} \bar{\mu}_{g, t}^{\text {down }} S_{g, t}^{\text {down }, o}-\sum_{g=1}^{G} \mu_{g, t}^{\text {up }} P_{g, t}^{o}+ \\ \lambda_{t}^{e}\left(P_{g, t}-P_{w, t}\right)+\left(\lambda_{t}^{\mathrm{up}} S_{t}^{\mathrm{up}}+\lambda_{t}^{\mathrm{down}} S_{t}^{\mathrm{down}}\right)- \\ \sum_{w=1}^{W} V_{w, t} l_{w, t, k}-\sum_{g=1}^{G} \gamma_{g, t}^{\mathrm{e}, \mathrm{o}}\left(P_{g, t, k}^{\mathrm{up}}-P_{g, t, k}^{\mathrm{down}}\right)- \\ \sum_{w=1}^{W} \bar{\rho}_{w, t, k} P_{w, t}-\sum_{g=1}^{G}\left(\bar{\rho}_{g, t, k}^{\text {up }} S_{g, t}^{\text {up }}+\right. \\ \left.\left.\bar{\rho}_{g, t, k}^{\text {down }} S_{g, t}^{\text {down }}\right)+\lambda_{t, k}^{b} P_{t, k}\right] \end{array}$

然而,观察式(79)可知,其中仍存在双变量乘积项ρ-g,t.kupSg,tupρ-g,t.kdownSg,tdown,此处采用二元展开法对其进行近似化线性处理.接下来以线性化ρ-g,t.kupSg,tup为例进行说明.

决策变量Sg,tup的上下限已知,为[Sg,t,minup, Sg,t,maxup],故可利用二进制变量ug,tup对其进行展开:

Sg,tup= Sg,t,minupSg,tupi=1I2i-1ug,tup
Δ Sg,tup=(Sg,t,maxup- Sg,t,minup)/2I

式中:ΔSg,tup为线性化展开步长;I为展开项数,决定线性化的精度,以增加计算量及计算时间为代价,求解精度随着I的增大而提升.

式(80)两边同乘以ρ-g,t.kup,可转化为

ρ-g,t.kupSg,tup=ρ-g,t.kupSg,t,minup+ ρ-g,t.kupΔ Sg,tupi=1I2i-1ug,tup

定义变量Xρu=ρ-g,t.kupug,tup,则式(82)可转化为

ρ-g,t.kupSg,tup= ρ-g,t.kupSg,t,minupSg,tupi=1I2i-1Xρu

Xρu=ρ-g,t.kupug,tup可进一步借助大M法进行线性化表述为

0≤ ρ-g,t.kup-Xρu≤M ug,tup
0≤Xρu≤M(1- ug,tup)

4 算例仿真及分析

4.1 算例设置

为验证所提模型的有效性,设置电力系统仿真场景进行仿真验证.系统中包含1个装机容量为200 MW的风电场W1,1个储能电站E1,5台传统机组G1~G5.其中,系统负荷需求及风电场预测出力如图5所示.图5中阴影区域为系统净负荷,即系统总负荷与风电场出力的差值.同时净负荷偏差考虑全部由风电场出力不确定性引起的预测偏差产生,系统通过调用日前备用辅助服务市场中标的储能及传统机组来对净负荷偏差进行消除,满足实时潮流平衡.

图5

图5   系统负荷需求及风电预测出力

Fig.5   System load demand and predicted wind power output


储能规模设置为40 MW/160 MW·h,最大上调备用功率Seup,max和最大下调备用功率Sedown,max均为15 MW,充电效率Pch和放电效率Pdis均为0.95,初始荷电状态Eeinit设为0.2Eerated.参照文献[27],传统机组G1~G5的CEA配额值均为0.877 t/(MW·h),其实际CEA分别设为1.087、0.727、0.887、0.861、1.083 t/(MW·h),碳交易单价定为14 美元/t,传统机组其余参数设置如表1所示.系统日前备用辅助服务市场上调备用和下调备用需求如图6所示.

表1   传统机组参数设置

Tab.1  Parameter settings of traditional unit

机组γg,te,o/
[美元·(MW·h)-1]
γg,tec,o/
[美元·(MW·h)-1]
γg,ts,o/
[美元·(MW·h)-1]
Sg,tup,o/MWSg,tdown,o/MWPog,t/MW
G12022.94061212100
G22724.90012101060
G34040.14015101040
G44544.776208840
G56062.884365560

新窗口打开| 下载CSV


图6

图6   系统日前市场总备用容量需求

Fig.6   Total reserve capacity demand of the system in current market


4.2 仿真结果及分析

(1) 各竞标主体市场收益及中标总体情况.各竞标主体在电能量市场、备用辅助服务市场及实时平衡市场的总体收益情况如表2所示.由表可见,电源侧储能与传统机组的主要收益来源均是电能量收益,这也从侧面体现电源侧发电机组作为能量供应者的基本角色定位.此外,电源侧储能综合市场与自身情况进行策略性竞标,因而在各市场中均有中标,并获得可观收益.

表2   各竞标主体市场收益情况

Tab.2  Profitability of each bidding entity in electricity market

竞标
主体
电能量
收益/
美元
上调备
用收益/
美元
下调备
用收益/
美元
实时平衡
市场收益/
美元
合计/
美元
E18 273.631 7341 2964 227.0215 530.65
G194 411.184741 044555.7796 484.95
G251 003.90120014.3551 138.25
G325 045.46120020.0725 185.53
G417 128.1400017 128.14
G58 316.820008 316.82

新窗口打开| 下载CSV


各竞标主体总体中标情况如图7所示.由图可见,电能量市场中传统机组的中标量占绝大多数,与传统机组仍是电能量提供的主力实际情况相吻合.而在备用辅助服务市场中,储能中标量占比更大,这主要得益于储能可以通过改变自身充放电状态,灵活提供上下调备用,因而能够在市场竞争中占据优势.综合表2图7可知,各主体在市场中的收益情况与其中标量表现一致,在市场统一出清下,中标量大的主体收益也更大.对于未竞标成功的机组则市场收益为0,如日前备用辅助服务市场中的机组G4、G5,结合市场出清规则来看,主要源于其高报价难以被优先安排出清,这也进一步凸显出市场主体应注重降本增效,从而可在低生产成本下得以进行较低报价,以提高中标机会.

图7

图7   各竞标主体市场总体中标情况

Fig.7   Overall winning situation of each bidding entity in electricity market


(2) 各竞标主体各时段中标情况及市场出清电价.图8~10分别展示各竞标主体在不同市场各时段具体中标情况.由图8可更清晰地看出,储能在电能量市场各时段中标量占比较少,传统机组G1的中标量则在所有竞价主体中占比最大.市场呈现此种出清表现,一方面与G1机组报价较低,得以优先出清,且其中标量已能覆盖系统的大部分能量需求有关;另一方面也源于储能容量有限且其需要对自身能量及功率进行全局规划.此外,由图8还可看出储能能够在电能量市场适时切换放能、用能角色,以对能量进行调配,如 3:00—5:00、14:00—17:00.

图8

图8   各竞标主体电能量市场中标情况

Fig.8   Winning situation of each bidding entity in electric energy market


图9

图9   各竞标主体备用辅助服务市场中标情况

Fig.9   Winning situation of each bidding entity in reserve ancillary service market


图10

图10   各竞标主体实时平衡市场中标情况

Fig.10   Winning situation of each bidding entity in real-time balanced market


图9图10则进一步体现在备用辅助服务市场中可灵活切换充放电状态的储能是调节主力.此外,对比图9及和图10还可看出,在备用辅助服务市场中标的备用容量未必在实时平衡市场会得到全部调用,具体调用情况取决于电力系统在实时运行过程中的实际需求.

图11所示为各时段电力市场出清电价,可看出在电力负荷需求较高、供电紧张的 9:00—12:00、18:00—20:00 时段,电能量市场出清电价也较高,充分体现供求关系下电力的商品价值属性.

图11

图11   电力市场出清电价

Fig.11   Electricity market clearing price


(3) 电源侧储能竞标策略.图12图13所示分别为电源侧储能在电能量市场和备用辅助服务市场中的竞标策略,即时序耦合下的功率报量和能量报价.其中,由图12可以看出,储能中标充电功率时,如 1:00、3:00—5:00、8:00、12:00—13:00,会以一个极高值进行报价,在市场对购电方报价由高到低排序作为优先出清顺序的机制下,该策略无疑可以提高充电功率中标机会.

图12

图12   电源侧储能电能量市场竞标策略

Fig.12   Bidding strategy of power side energy storage in electric energy market


图13

图13   电源侧储能备用辅助服务市场竞标策略

Fig.13   Bidding strategy of power side energy storage in reserve ancillary service market


4.3 碳单价参数灵敏度分析

碳交易成本主要附加在传统机组电能量市场报价上,故在此分别设置碳单价为0、14、28、42 美元/t,重点研究碳单价改变引起的电能量市场运行结果变化.不同碳单价下各机组中标量及电能量市场出清电价分别如图14图15所示.

图14

图14   不同碳单价下各机组电能量中标量

Fig.14   Winning bid of electric energy for each unit at different carbon unit prices


图15

图15   不同碳单价下电能量市场出清电价

Fig.15   Clearing prices in the electricity market at different carbon unit prices


随着碳单价的增加,传统机组的报价将整体发生变化,因而不同碳单价下节点边际电价呈现相近的变化趋势,电源侧储能的总中标量趋向相同,如图14所示.但随着碳单价增加,G1中标量逐渐减小,G2中标量逐渐增加,高碳排放量的G1出力逐渐向较低碳排放量的G2发生转移;其原因是高碳排放机组报价G1随碳单价的增加而增大,而G2由于其碳排放实际值低于配额值,其报价与碳单价变化呈负相关,且当碳单价升高到一定范围,G2报价会低于G1报价,此时G2取代G1,被优先出清.机组电-碳耦合报价大小关系无改变情况下,其出清顺序也不改变.G3~G5中标总量无明显变化主要源于所测试的碳单价还未增大到一定程度,以至于其引起的附加碳交易成本还未覆盖甚至超过G3~G5之间的电能量报价差,故在负荷需求及报价大小关系不变下,G3~G5的中标总量基本不变.

图15可得出,当碳单价不断升高时,传统机组报价相应升高,引起市场出清电价随之升高.因此,电源侧储能在进行投标决策时应充分了解相关市场机制,紧跟市场动态,有助于己方进行更加合理且更高收益的策略竞标.

5 结论

通过分析电力市场交易机制,从电源侧储能角度出发,建立计及风电出力不确定性及碳交易的储能参与电力市场多场景应用竞标博弈双层优化模型,指导优化电源侧储能竞标策略,提升其多场景应用综合经济效益.同时,进行算例仿真并得到以下结论.

(1) 基于博弈关系所构建的双层优化模型,借助数学变换,可将复杂的多目标、非线性市场竞标优化问题转化为单一目标的线性优化求解问题,有效反映出实际出清过程中电源侧储能与市场的博弈关系.

(2) 所提模型能以电源侧储能综合收益最大为目标,兼顾社会购电成本及风电出力不确定性,完成对电源侧储能在电力市场的竞标策略优化及市场模拟出清,实现多方主体互利共赢.

(3) 因充当的市场角色不同,电源侧储能在不同市场中的竞价策略存在差异,在电能量市场中充当用能角色时会以极高报价提高充电中标概率,而充当电源角色进行放电或在备用辅助服务市场发挥调节作用时,其报价多与优化模拟的市场出清电价相近.

(4) 计及碳交易后,碳单价的改变对电源侧储能在电能量市场中的中标量无明显影响,但会使得传统机组的中标量在不同碳排放强度机组间发生转移.同时,碳单价升高会引起电能量市场出清电价升高,进而使得电源侧储能在电能量市场的中标电价同步升高.

本文基于我国负荷侧报量不报价的现行电力市场规则,在电源侧储能参与电力市场竞标博弈建模中对负荷侧的投标方式同样采用报量不报价模式.然而,随着电力市场改革的深入推进,为充分发挥市场价格发现及调控作用,势必会向负荷侧报量报价转变.因此,下一步将考虑负荷侧报量报价的竞标机制,对目前的优化模型加以改进.

参考文献

黄强, 郭怿, 江建华, .

“双碳” 目标下中国清洁电力发展路径

[J]. 上海交通大学学报, 2021, 55(12): 1499-1509.

DOI:10.16183/j.cnki.jsjtu.2021.272      [本文引用: 1]

目前,第三次能源革命已经开始,为了减少碳排放,发达国家先后制定了清洁能源发展战略,并公布了放弃火电、核电的时间表.与此同时,中国也向世界作出承诺:2030年前碳排放达峰,2060年实现碳中和.因此,在&#x0201c;双碳&#x0201d;目标下,研究清洁电力发展路径具有重要意义.分析了中国水、风、光等清洁能源储量及其特征,预测了中长期电力需求,依据电力电量平衡原理,估算了2030和2050规划水平年电力系统结构组成,并分析了未来CO<sub>2</sub>排放趋势,提出了未来中国清洁电力的发展对策与建议.结果表明:预计2027年中国电力系统将实现&#x0201c;碳达峰&#x0201d;;2030年中国清洁电力发电量将超过总发电量的50%;2050年火电、核电将被水、风、光等清洁电力全部取代,电力行业将实现CO<sub>2</sub>的&#x0201c;零排放&#x0201d;,基本全面实现电力系统的绿色转型,以响应国家的&#x0201c;双碳&#x0201d;目标.

HUANG Qiang, GUO Yi, JIANG Jianhua, et al.

Development pathway of China’s clean electricity under carbon peaking and carbon neutrality goals

[J]. Journal of Shanghai Jiao Tong University, 2021, 55(12): 1499-1509.

[本文引用: 1]

SHI J, LEE W J, LIU X F.

Generation scheduling optimization of wind-energy storage system based on wind power output fluctuation features

[C]// 2017 IEEE/IAS 53rd Industrial and Commercial Power Systems Technical Conference. Niagara Falls, Canada: IEEE, 2017: 1-7.

[本文引用: 1]

SIOSHANSI R, DENHOLM P, ARTEAGA J, et al.

Energy-storage modeling: State-of-the-art and future research directions

[J]. IEEE Transactions on Power Systems, 2022, 37(2): 860-875.

[本文引用: 1]

李姚旺, 张宁, 张世旭, .

面向电力系统的多能源云储能模式:基本概念与研究展望

[J]. 中国电机工程学报, 2023, 43(6): 2179-2190.

[本文引用: 1]

LI Yaowang, ZHANG Ning, ZHANG Shixu, et al.

Multi-energy cloud energy storage for power systems: Basic concepts and research prospects

[J]. Proceedings of the CSEE, 2023, 43(6): 2179-2190.

[本文引用: 1]

陈国平, 梁志峰, 董昱.

基于能源转型的中国特色电力市场建设的分析与思考

[J]. 中国电机工程学报, 2020, 40(2): 369-379.

[本文引用: 1]

CHEN Guoping, LIANG Zhifeng, DONG Yu.

Analysis and reflection on the marketization construction of electric power with Chinese characteristics based on energy transformation

[J]. Proceedings of the CSEE, 2020, 40(2): 369-379.

[本文引用: 1]

陈启鑫, 房曦晨, 郭鸿业, .

储能参与电力市场机制:现状与展望

[J]. 电力系统自动化, 2021, 45(16): 14-28.

[本文引用: 1]

CHEN Qixin, FANG Xichen, GUO Hongye, et al.

Participation mechanism of energy storage in electricity market: Status quo and prospect

[J]. Automation of Electric Power Systems, 2021, 45(16): 14-28.

[本文引用: 1]

闫斯哲, 王维庆, 李笑竹, .

储能-机组联合调频的动态经济环境跨区灵活性鲁棒优化调度

[J]. 电力系统自动化, 2022, 46(9): 61-70.

[本文引用: 1]

YAN Sizhe, WANG Weiqing, LI Xiaozhu, et al.

Cross-regional flexible robust optimal scheduling in dynamic economic environment with joint frequency regulation of energy storage and units

[J]. Automation of Electric Power Systems, 2022, 46(9): 61-70.

[本文引用: 1]

ZHOU B, FANG J K, AI X M, et al.

Storage right-based hybrid discrete-time and continuous-time flexibility trading between energy storage station and renewable power plants

[J]. IEEE Transactions on Sustainable Energy, 2023, 14(1): 465-481.

[本文引用: 1]

ABDELTAWAB H, MOHAMED Y A R I.

Energy storage planning for profitability maximization by power trading and ancillary services participation

[J]. IEEE Systems Journal, 2022, 16(2): 1909-1920.

[本文引用: 1]

尚楠, 陈政, 冷媛.

电碳市场背景下典型环境权益产品衔接互认机制及关键技术

[J]. 中国电机工程学报, 2024, 44(7): 2558-2577.

[本文引用: 1]

SHANG Nan, CHEN Zheng, LENG Yuan.

Mutual recognition mechanism and key technologies of typical environmental interest products in power and carbon markets

[J]. Proceedings of the CSEE, 2024, 44(7): 2558-2577.

[本文引用: 1]

张巍, 缪辉.

基于云储能租赁服务的风储参与能量-调频市场竞价策略研究

[J]. 电网技术, 2021, 45(10): 3840-3850.

[本文引用: 1]

ZHANG Wei, MIAO Hui.

Bidding strategies of wind power and energy storage participating in energy and frequency regulation market based on cloud energy storage leasing services

[J]. Power System Technology, 2021, 45(10): 3840-3850.

[本文引用: 1]

孙鑫, 魏文荣, 李琼林, .

计及不确定性的风储联合系统多时间尺度市场参与策略

[J]. 电力自动化设备, 2024, 44(2): 18-25.

[本文引用: 1]

SUN Xin, WEI Wenrong, LI Qionglin, et al.

Multi-time scale market participation strategy of wind-energy storage combined system considering uncertainty

[J]. Electric Power Automation Equipment, 2024, 44(2): 18-25.

[本文引用: 1]

王霞. 含储能的能量与一次调频服务市场联合优化出清研究[D]. 武汉: 武汉大学, 2021.

[本文引用: 1]

WANG Xia. Research on the joint optimization clearing for eneregy and primary frequency regulation markets incorporating energy storage[D]. Wuhan: Wuhan University, 2021.

[本文引用: 1]

葛晓琳, 凡婉秋, 符杨, .

基于改进柔性策略评价的风火储多主体博弈电能-调频市场联合竞价模型

[J]. 电网技术, 2023, 47(5): 1920-1933.

[本文引用: 1]

GE Xiaolin, FAN Wanqiu, FU Yang, et al.

Joint bidding model of electricity and frequency regulation market with wind fire storage multi-agent games based on improved soft actor-critic

[J]. Power System Technology, 2023, 47(5): 1920-1933.

[本文引用: 1]

KAZEMI M, ZAREIPOUR H, AMJADY N, et al.

Operation scheduling of battery storage systems in joint energy and ancillary services markets

[J]. IEEE Transactions on Sustainable Energy, 2017, 8(4): 1726-1735.

[本文引用: 1]

MORENO R, MOREIRA R, STRBAC G.

A MILP model for optimising multi-service portfolios of distributed energy storage

[J]. Applied Energy, 2015, 137: 554-566.

[本文引用: 1]

AKHAVAN-HEJAZI H, MOHSENIAN-RAD H.

Optimal operation of independent storage systems in energy and reserve markets with high wind penetration

[J]. IEEE Transactions on Smart Grid, 2014, 5(2): 1088-1097.

[本文引用: 1]

车泉辉, 吴耀武, 祝志刚, .

基于碳交易的含大规模光伏发电系统复合储能优化调度

[J]. 电力系统自动化, 2019, 43(3): 76-82.

[本文引用: 1]

CHE Quanhui, WU Yaowu, ZHU Zhigang, et al.

Carbon trading based optimal scheduling of hybrid energy storage system in power systems with large-scale photovoltaic power generation

[J]. Automation of Electric Power Systems, 2019, 43(3): 76-82.

[本文引用: 1]

YAN N, MA G C, LI X J, et al.

Low-carbon economic dispatch method for integrated energy system considering seasonal carbon flow dynamic balance

[J]. IEEE Transactions on Sustainable Energy, 2023, 14(1): 576-586.

[本文引用: 1]

MU C G, DING T, ZHU S Y, et al.

A decentralized market model for a microgrid with carbon emission rights

[J]. IEEE Transactions on Smart Grid, 2023, 14(2): 1388-1402.

[本文引用: 1]

张钧钊, 姜欣, 段世杰, .

虚拟电厂参与电-碳联合市场运行的竞价策略研究

[J]. 电力系统保护与控制, 2023, 51(11): 108-118.

[本文引用: 1]

ZHANG Junzhao, JIANG Xin, DUAN Shijie, et al.

Bidding strategy for a virtual power plant to participate in the power-carbon joint market

[J]. Power System Protection & Control, 2023, 51(11): 108-118.

[本文引用: 1]

马云聪, 武传涛, 林湘宁, .

计及碳排放权交易的光热电站市场竞价策略研究

[J]. 电力系统保护与控制, 2023, 51(4): 82-92.

[本文引用: 1]

MA Yuncong, WU Chuantao, LIN Xiangning, et al.

Bidding strategy for a concentrated solar power plant participating in the electricity market with the background of carbon trading

[J]. Power System Protection & Control, 2023, 51(4): 82-92.

[本文引用: 1]

米阳, 赵海辉, 付起欣, .

考虑风光不确定与碳交易的区域综合能源系统双层博弈优化运行

[J]. 电网技术, 2023, 47(6): 2174-2188.

[本文引用: 1]

MI Yang, ZHAO Haihui, FU Qixin, et al.

Two-level game optimal operation of regional integrated energy system considering wind and solar uncertainty and carbon trading

[J]. Power System Technology, 2023, 47(6): 2174-2188.

[本文引用: 1]

RUIZ C, CONEJO A J.

Pool strategy of a producer with endogenous formation of locational marginal prices

[J]. IEEE Transactions on Power Systems, 2009, 24(4): 1855-1866.

[本文引用: 1]

潘虹锦, 高红均, 杨艳红, .

基于主从博弈的售电商多元零售套餐设计与多级市场购电策略

[J]. 中国电机工程学报, 2022, 42(13): 4785-4800.

[本文引用: 1]

PAN Hongjin, GAO Hongjun, YANG Yanhong, et al.

Multi-type retail packages design and multi-level market power purchase strategy for electricity retailers based on master-slave game

[J]. Proceedings of the CSEE, 2022, 42(13): 4785-4800.

[本文引用: 1]

姜欣, 刘萌, 王天梁, .

电网侧储能电站参与现货联合市场的竞价策略

[J]. 电网技术, 2021, 45(9): 3398-3408.

[本文引用: 1]

JIANG Xin, LIU Meng, WANG Tianliang, et al.

Bidding strategy for grid-side energy storage power stations to participate in the spot joint market

[J]. Power System Technology, 2021, 45(9): 3398-3408.

[本文引用: 1]

王凯, 延肖何, 蒋凯, .

考虑碳交易的风光储场站参与电力现货市场报价策略与调控方法

[J]. 中国电机工程学报, 2023, 43(18): 7091-7103.

[本文引用: 2]

WANG Kai, YAN Xiaohe, JIANG Kai, et al.

Bidding strategy and regulation method for the unified wind/photovoltaic/energy storage power stations in electricity spot market considering carbon trading

[J]. Proceedings of the CSEE, 2023, 43(18): 7091-7103.

[本文引用: 2]

NASROLAHPOUR E, KAZEMPOUR J, ZAREIPOUR H, et al.

A bilevel model for participation of a storage system in energy and reserve markets

[J]. IEEE Transactions on Sustainable Energy, 2018, 9(2): 582-598.

[本文引用: 1]

TOUBEAU J F, BOTTIEAU J, DE GRÈVE Z, et al.

Data-driven scheduling of energy storage in day-ahead energy and reserve markets with probabilistic guarantees on real-time delivery

[J]. IEEE Transactions on Power Systems, 2021, 36(4): 2815-2828.

[本文引用: 1]

/