CN112766695B - 一种负荷聚集主体参与下的主辅联合***的均衡运行方法 - Google Patents
一种负荷聚集主体参与下的主辅联合***的均衡运行方法 Download PDFInfo
- Publication number
- CN112766695B CN112766695B CN202110038800.2A CN202110038800A CN112766695B CN 112766695 B CN112766695 B CN 112766695B CN 202110038800 A CN202110038800 A CN 202110038800A CN 112766695 B CN112766695 B CN 112766695B
- Authority
- CN
- China
- Prior art keywords
- market
- main body
- frequency modulation
- load
- capacity
- 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.)
- Active
Links
- 230000002776 aggregation Effects 0.000 title claims abstract description 123
- 238000004220 aggregation Methods 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003828 downregulation Effects 0.000 claims abstract description 16
- 238000010248 power generation Methods 0.000 claims description 78
- 230000005611 electricity Effects 0.000 claims description 26
- 230000003827 upregulation Effects 0.000 claims description 21
- 230000000295 complement effect Effects 0.000 claims description 15
- 238000007599 discharging Methods 0.000 claims description 11
- 238000012887 quadratic function Methods 0.000 claims description 9
- 230000032683 aging Effects 0.000 claims description 6
- 238000004146 energy storage Methods 0.000 claims description 6
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000003169 complementation method Methods 0.000 claims description 3
- 230000002222 downregulating effect Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 8
- 230000009467 reduction Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000000630 rising effect Effects 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
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- 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/0202—Market predictions or forecasting for commercial 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
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Finance (AREA)
- Accounting & Taxation (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Educational Administration (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开的负荷聚集主体参与下的主辅联合***的均衡运行方法，主辅联合***接入可调负荷资源，各市场主体依据市场均衡解来调节各自的机组出力和可调负荷资源使得主辅联合***均衡运行；以负荷聚集主体为对象，通过主辅联合市场模拟分析，对其市场运营能力进行探索；先建立负荷聚集主体参与下的电能量市场和调频辅助服务市场交易框架，分析负荷聚集主体参与市场方式及各市场主体的成本，基于博弈均衡理论建立负荷聚集主体参与下的市场多时段均衡模型求得衡模型的均衡解；令负荷聚集主体参与主辅市场竞争，对负荷及有功电价均有“削峰填谷”的作用，使***调频容量电价均有所下降，且下调容量电价下降更为明显；降低***调频服务的购买成本。
Description
技术领域
本发明涉及市场能源分配技术领域，具体涉及一种负荷聚集主体参与下的主辅联合***的均衡运行方法。
背景技术
随着环境问题的加剧以及需求侧响应技术的发展，未来大量需求侧负荷资源将加入市场进行交易，大规模中小负荷的参与将对电网运行提出更大的挑战。利用负荷侧资源分布广泛、快速响应的特点，负荷聚集主体(Load Aggregator，LA)可以对包括中小型用户资源、储能资源、电动汽车等负荷侧资源进行统一控制，并参与现货市场竞争。通过一定的调控策略，LA不仅能够抑制、消除不确定性负荷对电网带来的不利影响，还能利用负荷侧可调资源提供调频、调峰等辅助服务。
LA作为运营主体，在现货市场竞标中主要包括两方面策略，其一是各类型负荷每个时段的用电需求及报价，不同的市场策略将对市场出清电价及自身收益造成影响；另一方面，每个时段的电量和电价确定以后，各类负荷的可调容量也将确定，如何利用此部分可调容量参与辅助服务市场将决定LA的另一部分收益。LA作为新的市场主体参与到现货市场竞争中，对主、辅市场出清结果及各方市场主体收益带来的影响值得深入研究。
目前国内外在这方面有一些初步的研究工作。在资源控制策略研究方面，有学者从负荷聚集主体角度、投资者角度提出了响应的调度、收益模型，并证明可调负荷资源可以很好的满足调频需求，并且对电网安全稳定有着积极的作用；在电力市场竞争方面，有学者运用寡头博弈均衡理论研究各类市场行为对各市场参与者收益及市场经济效率的影响，分析了负荷聚集主体作为价格接受者和价格影响者参与电力市场的情况，以及风电主体和负荷聚集主体以VPP合作模式参与投标竞争的电力市场多时段随机博弈模型。但目前的研究中，负荷聚集主体仅作为价格接受方，并未考虑大量可调负荷资源接入后对其他市场主体策略及整个现货市场带来的影响，市场框架只包含电能量市场，而实际现货市场必定是主辅联合市场，因此无法全面体现负荷侧对现货市场带来的影响。
发明内容
为解决上述技术问题，本发明基于寡头竞争的博弈均衡理论，研究负荷聚集主体作为市场主体同时参加主辅现货市场下的均衡结果，首先建立主辅联合现货市场交易框架并对各市场主体参与调频市场的成本建模分析，在此基础上建立负荷聚集主体参与下的能量和调频市场多时段均衡模型，通过算例仿真比较有、无负荷聚集主体的参与对现货市场出清结果及各方收益的影响，并对可调负荷的调节性能进行定量灵敏度分析。
本发明通过下述技术方案实现：
本方案提供一种负荷聚集主体参与下的主辅联合***的均衡运行方法，所述主辅联合***为电能量市场和调频辅助服务市场联合***，所述主辅联合***接入可调负荷资源，电能量市场和调频辅助服务市场的市场主体均包括负荷聚集主体和传统发电主体，电能量市场和调频辅助服务市场的市场主体依据市场均衡解来调节可调负荷资源使得主辅联合***均衡运行，
所述电能量市场和调频辅助服务市场的市场主体依据市场均衡解调节可调负荷资源具体方法包括：
S1.建立负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架；
S2.在S1的交易框架下确定负荷聚集主体参与主辅联合***的方式，并依据市场主体交易成本模型对各市场主体的成本进行分析；
S3.基于S1和S2以及寡头竞争的博弈均衡理论，建立负荷聚集主体参与下的主辅联合***的均衡模型；
S4.应用非线性互补方法求解均衡模型得到市场均衡解；
S5.各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力。
所述可调负荷资源包括：用户负荷资源、储能负荷资源及电动汽车负荷资源。
进一步优化方案为，负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架，具体为：
电能量市场和调频辅助服务市场的各市场主体以古诺竞争方式参与；
参与主辅联合***中的各市场主体按自身策略分别向交易中心上报其在电能量市场、调频辅助服务市场中的竞争策略；
交易中心对上报的竞争策略进行检验，并给出电能量市场、调频辅助服务市场的出清结果和中标容量；
电能量市场、调频辅助服务市场均以出清结果及各主体中标容量进行结算。
市场均衡解包括：中标市场主体，各中标市场主体在电能市场的中标容量，各中标市场主体在调频辅助服务市场的上调或下调容量；
假定某初期日前现货市场有电能量市场和调频辅助服务市场两部分。市场中存在n个传统发电主体和1个负荷聚集主体，各市场主体以古诺竞争方式参与现货市场，按自身策略分别上报在主、辅市场中的竞争容量及报价，所上报的两部分容量需满足基本物理约束。交易中心将对上报容量进行检验并对主、辅市场进行出清。
各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力具体过程为：
调频辅助服务市场中的中标市场主体，在运行时段接收到上调(下调)指令后，提供上调(下调)容量的传统发电主体应增大(减少)机组出力，提供上调(下调)容量的负荷聚集主体应控制负荷资源进行出力增加或削减。
进一步优化方案为，电能量市场和调频辅助服务市场均采用独立出清方式。
调频市场已在欧美等成熟电力市场开展多年，市场出清规则主要以辅助服务购买成本最低为主，但出清方式主要存在独立出清和主-辅市场统一出清两种方式。就日前现货市场而言，美国PJM电力市场中，调频市场与电能量市场独立出清；美国加州电力市场中，主-辅市场实行统一出清。主-辅统一出清方式下，***总购电成本最低，但市场规则复杂、出清结果不透明，对市场运营要求较高。独立出清方式下，市场规则简单、透明度高，适合于建设初期的日前现货市场，目前我国各区域辅助服务市场试运行规则以独立出清为主，因此本方案采用独立出清方式。
主(电能量市场)、辅(调频辅助服务市场)市场均以出清结果电价及各主体中标容量进行结算，且调频容量结算时，会考虑中标机组的调频性能。本发明调频市场结算不考虑里程收益，因参与调频而产生的能量交换按中标时段能量市场电价进行结算。
进一步优化方案为，市场主体交易成本模型包括：传统发电主体的发电成本可用出力模型，传统发电主体参与调频辅助服务市场带来的附属成本模型、调频辅助服务市场中负荷聚集主体调频容量的成本模型、电能量市场中负荷聚集主体的出力模型。
进一步优化方案为，传统发电主体的发电成本可出力模型具体为：
传统发电主体i(i＝1,2,…,n)在t时段的发电成本可用出力的二次函数表示为：
式中：Pi,t为t时段第i个传统发电主体的有功出力，ai、bi为该传统发电主体的成本系数；
当传统发电主体参与调频时，机组运行过程中将不断接收调频信号并调整出力，由此造成的热效率损失成本较难核算[19]，在一个长期稳定的电力***中，这部分热效率损失成本可在实测数据的基础上对成本曲线进行拟合。因此在不考虑投资成本的情况下，传统发电主体参与调频辅助服务市场带来的附属成本也可拟合为调频容量的二次函数。
传统发电主体参与调频辅助服务市场带来的附属成本模型具体为：
传统发电主体参与调频辅助服务市场带来的附属成本拟合为调频容量的二次函数：
式中：
为第i个传统发电主体的上、下调频容量，
为第i个传统发电主体的上、下调频成本系数；
负荷聚集主体参与现货市场的成本主要包括负荷用电带来的电量成本、负荷削减带来的机会成本、资源调用带来的损耗成本以及储能、电动汽车中SOC偏移带来的电池老化成本[21]，由于本文模型中不考虑单量电动汽车的充放电模型，且模型中电动汽车调频容量相较于电池容量较小，SOC偏移成本对出清结果分析和最终结论不造成影响，因此不考虑SOC偏移成本调频辅助服务市场中负荷聚集主体调频容量的成本模型表示为：
负荷聚集主体上调频容量的成本模型：
负荷聚集主体下调频容量的成本模型：
式(3)中，第一、二项分别为调用上调容量时产生的损耗成本及电量结算成本；式(4)为调用下调容量产生的电量结算成本；其中，
为t时段***上、下调频调用率，c为资源调用的老化成本。
在电能量市场中，负荷聚集主体须通过购电满足所辖负荷次日的用电需求，同时可利用其代理的电动汽车、储能等***的剩余容量在负荷低谷时段(电价较低)进行储能，在负荷高峰时段(电价较高)进行放电，在完成充放电套利的同时，对***负荷也有削峰填谷的作用。
电能量市场中负荷聚集主体的出力模型表示为：
式中：
Lt分别为第t个时段负荷聚集主体的套利充、放电量及为满足负荷的需求电量。
进一步优化方案为，所述多时段均衡模型考虑将一天分为24个时段，主要包括市场出清规则模型和市场主体决策模型。
由市场供需平衡条件可得，t时段电能量市场出清规则如式(6)所示，根据每个时段传统发电主体和负荷聚集主体的决策功率决定能量市场出清电价及出清电量。
进一步优化方案为，市场出清规则模型有：
t时段电能量市场出清规则模型表示为：
其中有：
Dt＝αt-βtλt (7)
式中Dt为第t个时段的电力需求，αt、βt为大于零的常数，λt为t时段有功电价；
调频辅助服务市场以购买成本最小进行出清，根据每个传统发电主体的决策报价确定调频市场出清价格和出清电量，t时段调频辅助服务市场出清规则模型为：
式(8)对应的约束条件为：
式中：下标t均表示在t时段的变量；
为负荷聚集主体的上、下调频容量，
分别为第i个传统发电主体及负荷聚集主体的上、下调频容量策略报价，Pi max、Pi min分别为第i个传统发电主体的出力上下限，
为负荷聚集主体的最大充放电功率，δ为调频需求容量占负荷量的比例系数[15]；式(9)、(10)的拉格朗日乘子
分别为上、下调频容量的出清电价；
市场主体决策模型有：
传统发电主体i(i＝1,2,…,n)决策模型表示为：
式(15)约束条件为式(7)～(17)
式(15)中，第一项为能量市场及调频辅助服务市场的电量结算收益，第二项为调频辅助服务市场的容量收益，第三、四项分别为发电成本及调频成本，其中，Ai为传统发电主体i的调频性能指标，式(16)、(17)为传统发电主体i的爬坡约束，
分别为机组上、下爬坡速率；
负荷聚集主体的决策模型表示为：
式(18)的约束条件为：式(7)～(14)及
式(18)中，第一项为能量市场及调频辅助服务市场的电量结算收益，第二项为套利放电引起的损耗成本，第三项为调频辅助服务市场的容量收益，第四项、第五项为上、下调频成本；其中：ALA为负荷聚集主体的调频性能指标；式(19)表示负荷用电需求平衡；式(20)表示用于套利及辅助服务的充放电量应始终保持在负荷资源有效容量之内；其中：η为充放电转换系数，
为负荷资源有效容量。
虽然主辅市场独立出清，但是传统发电主体和负荷聚集主体需要对两个市场的投标决策进行取舍，两个市场的出清结果将同时影响市场主体的最终收益，因此两类市场主体的决策模型中，两个市场的出清规则都需考虑在约束条件中。
进一步优化方案为，应用非线性互补方法求解多时段均衡模型得到市场均衡解具体步骤为：
S41.列出电能量市场出清规则模型和调频辅助服务市场出清规则模型的KKT条件；
S42.将S41所得KKT条件作为约束条件加入到各个市场主体决策模型中，再次列出各个市场主体决策模型的KKT条件；
S43.联立所有市场主体的KKT条件并用非线性互补函数处理其中的互补条件，得到一个非线性方程组；
S44.运用Levenberg-Marquardt算法对非线性方程组进行求解得到市场均衡解。
本方案运用改进的Levenberg-Marquardt算法对非线性方程组进行求解，该算法为基于次梯度的半光滑的Newton算法，只需求解线性***解就可确定搜索方向，在大规模***求解中可快速收敛。
本发明与现有技术相比，具有如下的优点和有益效果：
本发明提出的负荷聚集主体参与下的主辅联合***的均衡运行方法，以负荷聚集主体为对象，通过主辅联合市场模拟分析，对其市场运营能力进行探索；首先建立负荷聚集主体参与下的电能量市场和调频辅助服务市场交易框架，对负荷聚集主体参与市场方式及各市场主体的成本进行分析，然后基于寡头竞争的博弈均衡理论，建立负荷聚集主体参与下的电能量市场和调频市场多时段均衡模型；最后，应用非线性互补方法对均衡模型进行求解，使负荷聚集主体参与主辅市场竞争，对负荷及有功电价均有“削峰填谷”的作用，使得***上、下调频容量电价均有所下降，且下调容量电价下降更为明显；***调频服务购买成本有所降低，且整日用户利润和社会福利均将增加。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。
图1为负荷聚集主体参与下的市场框架示意图；
图2为求解多时段均衡模型流程图；
图3为有、无负荷聚集主体参与下的电能量市场均衡价格示意图；
图4为有、无负荷聚集主体参与下的调频辅助服务市场均衡价格示意图；
图5为传统发电主体和负荷聚集主体中标有功出力示意图；
图6为各市场主体中标的上调频容量示意图；
图7为传统发电主体和负荷聚集主体下调容量示意图；
图8为社会福利随可调容量的变化示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1
一种负荷聚集主体参与下的主辅联合***的均衡运行方法，所述主辅联合***为电能量市场和调频辅助服务市场联合***，所述主辅联合***接入可调负荷资源，电能量市场和调频辅助服务市场的市场主体均包括负荷聚集主体和传统发电主体，电能量市场和调频辅助服务市场的市场主体依据市场均衡解来调节可调负荷资源使得主辅联合***均衡运行，
所述电能量市场和调频辅助服务市场的市场主体依据市场均衡解调节可调负荷资源具体方法包括：
S1.建立负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架；
S2.在S1的交易框架下确定负荷聚集主体参与主辅联合***的方式，并依据市场主体交易成本模型对各市场主体的成本进行分析；
S3.基于S1和S2以及寡头竞争的博弈均衡理论，建立负荷聚集主体参与下的主辅联合***的均衡模型；
S4.应用非线性互补方法求解均衡模型得到市场均衡解；
S5.各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力。
所述可调负荷资源包括：用户负荷资源、储能负荷资源及电动汽车负荷资源。
负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架，具体为：
电能量市场和调频辅助服务市场的各市场主体以古诺竞争方式参与；
参与主辅联合***中的各市场主体按自身策略分别向交易中心上报其在电能量市场、调频辅助服务市场中的竞争策略；
交易中心对上报的竞争策略进行检验，并给出电能量市场、调频辅助服务市场的出清结果和中标容量；
电能量市场、调频辅助服务市场均以出清结果及各主体中标容量进行结算。
市场均衡解包括：中标市场主体，各中标市场主体在电能市场的中标容量，各中标市场主体在调频辅助服务市场的上调或下调容量；
各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力具体过程为：
调频辅助服务市场中的中标市场主体在运行时段接收到上调或下调指令后，提供上调或下调容量的传统发电主体增大或减少机组出力，提供上调或下调容量的负荷聚集主体调度负荷资源进行出力增加或削减。
电能量市场和调频辅助服务市场均采用独立出清方式。
市场主体交易成本模型包括：传统发电主体的发电成本可用出力模型，传统发电主体参与调频辅助服务市场带来的附属成本模型、调频辅助服务市场中负荷聚集主体调频容量的成本模型、电能量市场中负荷聚集主体的出力模型。
传统发电主体的发电成本可出力模型具体为：
传统发电主体i(i＝1,2,…,n)在t时段的发电成本可用出力的二次函数表示为：
式中：Pi,t为t时段第i个传统发电主体的有功出力，ai、bi为该传统发电主体的成本系数；
传统发电主体参与调频辅助服务市场带来的附属成本模型具体为：
传统发电主体参与调频辅助服务市场带来的附属成本拟合为调频容量的二次函数：
式中：
为第i个传统发电主体的上、下调频容量，
为第i个传统发电主体的上、下调频成本系数；
调频辅助服务市场中负荷聚集主体调频容量的成本模型表示为：
负荷聚集主体上调频容量的成本模型：
负荷聚集主体下调频容量的成本模型：
式(3)中，第一、二项分别为调用上调容量时产生的损耗成本及电量结算成本；式(4)为调用下调容量产生的电量结算成本；其中，
为t时段***上、下调频调用率，c为资源调用的老化成本；
电能量市场中负荷聚集主体的出力模型表示为：
式中：
Lt分别为第t个时段负荷聚集主体的套利充、放电量及为满足负荷的需求电量。
所述多时段均衡模型考虑将一天分为24个时段，主要包括市场出清规则模型和市场主体决策模型。
市场出清规则模型有：
t时段电能量市场出清规则模型表示为：
其中有：
Dt＝αt-βtλt (7)
式(6)中Dt为第t个时段的电力需求，αt、βt为大于零的常数，λt为t时段有功电价；t时段调频辅助服务市场出清规则模型为：
式(8)对应的约束条件为：
式中：下标t均表示在t时段的变量；
为负荷聚集主体的上、下调频容量，
分别为第i个传统发电主体及负荷聚集主体的上、下调频容量策略报价，Pi max、Pi min分别为第i个传统发电主体的出力上下限，
为负荷聚集主体的最大充放电功率，δ为调频需求容量占负荷量的比例系数[15]；式(9)、(10)的拉格朗日乘子
分别为上、下调频容量的出清电价；
市场主体决策模型有：
传统发电主体i(i＝1,2,…,n)决策模型表示为：
式(15)约束条件为式(7)～(17)
式(15)中，第一项为能量市场及调频辅助服务市场的电量结算收益，第二项为调频辅助服务市场的容量收益，第三、四项分别为发电成本及调频成本，其中，Ai为传统发电主体i的调频性能指标，式(16)、(17)为传统发电主体i的爬坡约束，
分别为机组上、下爬坡速率；
负荷聚集主体的决策模型表示为：
式(18)的约束条件为：式(7)～(14)及
式(18)中，第一项为能量市场及调频辅助服务市场的电量结算收益，第二项为套利放电引起的损耗成本，第三项为调频辅助服务市场的容量收益，第四项、第五项为上、下调频成本；其中：ALA为负荷聚集主体的调频性能指标；式(19)表示负荷用电需求平衡；式(20)表示用于套利及辅助服务的充放电量应始终保持在负荷资源有效容量之内；其中：η为充放电转换系数，
为负荷资源有效容量。
应用非线性互补方法求解多时段均衡模型得到市场均衡解具体步骤为：
S41.列出电能量市场出清规则模型和调频辅助服务市场出清规则模型的KKT条件；
S42.将S41所得KKT条件作为约束条件加入到各个市场主体决策模型中，再次列出各个市场主体决策模型的KKT条件；
S43.联立所有市场主体的KKT条件并用非线性互补函数处理其中的互补条件，得到一个非线性方程组；
S44.运用Levenberg-Marquardt算法对非线性方程组进行求解得到市场均衡解。
本实施例使用非线性互补方法求解，求解流程如图2所示。首先列出主、辅助市场出清模型(式(7)～(14))的KKT条件，由于KKT条件中有形如：a≥0，b≥0，ab＝0的互补条件，利用非线性互补函数ψ(a,b)等效替换该互补条件。ψ(a,b)表达式为：
然后将所得KKT条件作为约束条件，加入到各个市场主体决策模型中，再次列出所有传统发电主体决策模型的KKT条件，联立所有传统发电主体的KKT条件并用同样的非线性互补函数处理其中的互补条件，得到一个非线性方程组，通过求解该非线性方程组得到市场的均衡解。
本文运用改进的Levenberg-Marquardt算法对非线性方程组进行求解，该算法为基于次梯度的半光滑的Newton算法，只需求解线性***解就可确定搜索方向，在大规模***求解中可快速收敛。
另外，由于在理论上很难保证电力市场博弈问题均衡解的唯一性，本文算例所求得的均衡可能只是所有均衡解中的一个。
实施例2
为验证实施例1所给出的针对负荷聚集主体参与下的主辅联合市场的均衡分析方法的有效性，本实施例采用以下算例进行分析。
一、场景描述
考虑某日前现货市场中存在2个传统发电主体(G1，G2)和1个负荷聚集主体，且均同时参与现货市场中的电能量市场和调频辅助服务市场。传统发电主体的参数如表1所示。假设某日中长期电量分解后，现货市场需求约为总负荷的十分之一，假设能量市场线性需求函数Dt＝αt-βtλt中，各时段βt取恒定值0.6MW/USD，αt如表2所示。负荷总需求Ld为20MW·h，负荷资源有效容量
为80MW·h，每个时段的充放电限度为20MW，充放电转换系数η为0.85，单位放电老化成本系数c为60USD/(MW·h)。
辅助服务上、下调频容量需求量占负荷量的比例系数δ取4％，***上、下调频容量调用率
均取0.2。G1、G2、LA的调频性能指标均为1。
表1传统发电主体的参数
表2各时段需求函数参数αt
二、优化结果
1.不同情况下电能量市场均衡价格
图3给出了有、无负荷聚集主体参与下的电能量市场均衡价格。可以看出：电能量市场均衡价格变化趋势与负荷变化趋势一致；负荷聚集主体的参与对有功电价具有削峰填谷的作用，这是由于负荷聚集主体在夜间充电增加了负荷量，在日间放电增加了供电量。
2.不同情况下调频辅助服务市场均衡价格
图4给出了有、无负荷聚集主体参与下的调频辅助服务市场均衡价格。首先，由于传统发电主体的上、下调频成本参数相同，在没有负荷聚集主体参与现货市场时，调频市场上、下调频容量出清电价相同；当负荷聚集主体参与现货市场时，***上、下调频容量电价在各时段均有所下降，这是由于负荷聚集主体的加入削减了传统发电主体在调频辅助服务市场中的市场力；由于电量结算对调频成本的影响，负荷聚集主体上调容量成本更低，市场竞争力更强，对市场价格的降低效果更为明显，因而上调容量电价低于下调容量电价。
3.各市场主体中标容量
图5给出了传统发电主体和负荷聚集主体中标有功出力。可以看出：有功出力方面，负荷聚集主体在夜间充电，日间放电，较无序充电情况，对总负荷具有削峰填谷的作用。
图6给出了各市场主体中标的上调频容量。夜间阶段，负荷聚集主体中标上调容量高于传统发电主体，且负荷聚集主体中标容量的上升速度高于传统发电主体，这是由于上调容量带来的电量结算降低了负荷聚集主体的上调成本，而传统发电主体随着上调容量的增加，热效率损失更大，竞争力更低。高峰负荷阶段，由于LA套利放电，受物理约束限制，所能提供的上调容量有限，因此传统发电主体的中标容量上升明显。
图7给出了各市场主体中标的下调频容量。由于电量结算导致下调频容量成本较高，负荷聚集主体在下调频容量市场竞争力不如上调频容量市场，夜间阶段，由于套利充电，负荷聚集主体下调容量为满足物理约束的最大值；日间阶段，由于负荷聚集主体竞争力不足，中标上调容量少于传统发电主体。
4.各市场主体收益分析
表3给出了市场成员利润及社会福利结果。可以看出：负荷聚集主体“底充高发”，并利用可调容量参与调频，可以抑制传统发电主体在现货市场中的市场力，使其利润下降，并降低了***调频服务购买成本；另一方面，负荷聚集主体对有功电价具有削峰填谷的作用，使用户利润在低谷时段有所下降、高峰时段有所增加，但整日用户利润和社会福利均有所增加。
表3市场成员利润及社会福利结果
5.负荷可调容量灵敏度分析
图8给出了负荷可调容量变化对社会福利影响的曲线。可以看出，随着负荷可调容量的增加，社会福利逐渐增加，但速度逐渐放缓；当容量大于160MW·h后，增长几乎停止。这是因为，***充放电套利机会有限，当容量增大到一定值后，市场出清结果将不再变化，社会福利也将停止增长。
本实施例基于寡头竞争的博弈均衡理论，建立了负荷聚集主体参与下的能量和调频市场多时段均衡模型，通过算例分析，比较研究了负荷聚集主体的参与对现货市场均衡结果的影响，所得结论如下：
1)负荷聚集主体“低充高放”的套利充放电行为对***负荷有削峰填谷的作用，同时也将使低谷时段电价有所上升，高峰时段电价有所下降。
2)负荷聚集主体利用剩余容量参与调频辅助服务市场，使***上、下调频容量电价均有所下降。由于负荷聚集主体上调容量成本低于下调容量成本，负荷聚集主体中标的上调容量整体高于下调容量，对***上调容量电价的降低作用也更为明显。
3)负荷聚集主体参与现货市场将有效削弱其他市场成员的市场力，使传统发电主体的利润有所下降，同时***的调频服务购买成本也将降低；用户利润在低谷时段有所下降、高峰时段有所增加，且整日用户利润和社会福利均将增加。
4)负荷可调容量的增加将促进社会福利进一步增长，但增长趋势逐步放缓，最终趋于稳定。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种负荷聚集主体参与下的主辅联合***的均衡运行方法，所述主辅联合***为电能量市场和调频辅助服务市场联合***，其特征在于，所述主辅联合***接入可调负荷资源，电能量市场和调频辅助服务市场的市场主体均包括负荷聚集主体和传统发电主体，电能量市场和调频辅助服务市场的市场主体依据市场均衡解来调节可调负荷资源使得主辅联合***均衡运行，
所述电能量市场和调频辅助服务市场的市场主体依据市场均衡解调节可调负荷资源具体方法包括：S1.建立负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架；S2.在S1的交易框架下确定负荷聚集主体参与主辅联合***的方式，并依据市场主体交易成本模型对各市场主体的成本进行分析；S3.基于S1和S2以及寡头竞争的博弈均衡理论，建立负荷聚集主体参与下的主辅联合***的均衡模型；S4.应用非线性互补方法求解均衡模型得到市场均衡解；S5.各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力；
市场主体交易成本模型包括：传统发电主体的发电成本可用出力模型，传统发电主体参与调频辅助服务市场带来的附属成本模型、调频辅助服务市场中负荷聚集主体调频容量的成本模型、电能量市场中负荷聚集主体的出力模型；
传统发电主体的发电成本可出力模型具体为：
传统发电主体i；i=1,2,…,n；在t时段的发电成本可用出力的二次函数表示为：
式中：为t时段第i个传统发电主体的有功出力，a i 、b i 为该传统发电主体的成本系数；
传统发电主体参与调频辅助服务市场带来的附属成本模型具体为：
传统发电主体参与调频辅助服务市场带来的附属成本拟合为调频容量的二次函数：
式中：和为第i个传统发电主体的上、下调频容量， 表示第i个传统发电主体的上调频成本系数；表示第i个传统发电主体的下调频成本系数；
调频辅助服务市场中负荷聚集主体调频容量的成本模型表示为：
负荷聚集主体上调频容量的成本模型：
负荷聚集主体下调频容量的成本模型：
式中为负荷聚集主体上调频容量的成本；为负荷聚集主体的调用下调容量产生的电量结算成本； 为t时段***上、下调频调用率，c为资源调用的老化成本；为t时段有功电价；为负荷聚集主体的上调频容量；为负荷聚集主体的下调频容量； 电能量市场中负荷聚集主体的出力模型表示为：
式中：为负荷聚集主体的出力；分别为第t个时段负荷聚集主体的套利充、放电量及为满足负荷的需求电量。
2.根据权利要求1所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，负荷聚集主体参与下的电能量市场和调频辅助服务市场的交易框架，具体为：电能量市场和调频辅助服务市场的各市场主体以古诺竞争方式参与；参与主辅联合***中的各市场主体按自身策略分别向交易中心上报其在电能量市场、调频辅助服务市场中的竞争策略；交易中心对上报的竞争策略进行检验，并给出电能量市场、调频辅助服务市场的出清结果和中标容量；电能量市场、调频辅助服务市场均以出清结果及各主体中标容量进行结算。
3.根据权利要求2所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，电能量市场和调频辅助服务市场均采用独立出清方式。
4.根据权利要求1所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，均衡模型考虑将一天分为24个时段，主要包括市场出清规则模型和市场主体决策模型。
5.根据权利要求4所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，市场出清规则模型有： t时段电能量市场出清规则模型表示为：
其中有：
式中Dt为第t个时段的电力需求，为大于零的常数； t时段调频辅助服务市场出清规则模型为：
式（8）对应的约束条件为式（9）~（14）
式中：分别为第i个传统发电主体的上、下调频容量策略报价；分别为负荷聚集主体的上、下调频容量策略报价；
和为第i个传统发电主体的上、下调频容量；为负荷聚集主体的上调频容量；为负荷聚集主体的下调频容量；
为t时段第i个传统发电主体的有功出力；分别为第i个传统发电主体的出力上下限，为负荷聚集主体的最大充放电功率，为调频需求容量占负荷量的比例系数；
式（9）、（10）中的拉格朗日乘子分别为上、下调频容量的出清电价；
市场主体决策模型有： 传统发电主体i（i=1,2,…,n）决策模型表示为：
式（15）约束条件为式（16）~（17）
其中， A i 为传统发电主体i的调频性能指标，分别为机组上、下爬坡速率；为t时段***上、下调频调用率；分别为上、下调频容量的出清电价；为t时段的发电成本可用出力的二次函数表示；为传统发电主体参与调频辅助服务市场带来的附属成本拟合为调频容量的二次函数；
负荷聚集主体的决策模型表示为：
式（18）的约束条件为： 式（19）~（20）
其中：c为资源调用的老化成本；ALA为负荷聚集主体的调频性能指标； 为充放电转换系数，为负荷资源有效容量；分别为第t个时段负荷聚集主体的套利充、放电量及为满足负荷的需求电量；为负荷聚集主体上调频容量的成本；为负荷聚集主体的调用下调容量产生的电量结算成本。
6.根据权利要求1所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，应用非线性互补方法求解均衡模型得到市场均衡解具体步骤为：S41. 列出电能量市场出清规则模型和调频辅助服务市场出清规则模型的KKT条件；S42. 将S41所得KKT条件作为约束条件加入到各个市场主体决策模型中，再次列出各个市场主体决策模型的KKT条件；S43.联立所有市场主体的KKT条件并用非线性互补函数处理其中的互补条件，得到一个非线性方程组；S44.运用 Levenberg-Marquardt算法对非线性方程组进行求解得到市场均衡解。
7.根据权利要求2所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，市场均衡解包括：中标市场主体，各中标市场主体在电能市场的中标容量，各中标市场主体在调频辅助服务市场的上调或下调容量；各市场主体按照市场均衡解统一调度可调负荷资源出力和各市场主体的机组出力具体过程为：调频辅助服务市场中的中标市场主体在运行时段接收到上调或下调指令后，提供上调或下调容量的传统发电主体增大或减少机组出力，提供上调或下调容量的负荷聚集主体调度负荷资源进行出力增加或削减。
8.根据权利要求1所述的一种负荷聚集主体参与下的主辅联合***的均衡运行方法，其特征在于，所述可调负荷资源包括：用户负荷资源、储能负荷资源及电动汽车负荷资源。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110038800.2A CN112766695B (zh) | 2021-01-12 | 2021-01-12 | 一种负荷聚集主体参与下的主辅联合***的均衡运行方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110038800.2A CN112766695B (zh) | 2021-01-12 | 2021-01-12 | 一种负荷聚集主体参与下的主辅联合***的均衡运行方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112766695A CN112766695A (zh) | 2021-05-07 |
| CN112766695B true CN112766695B (zh) | 2023-04-07 |
Family
ID=75701681
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110038800.2A Active CN112766695B (zh) | 2021-01-12 | 2021-01-12 | 一种负荷聚集主体参与下的主辅联合***的均衡运行方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112766695B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113937788B (zh) * | 2021-10-11 | 2024-02-02 | 国网山东省电力公司营销服务中心(计量中心) | 一种电力负荷侧资源调度方法及*** |
| CN114140176B (zh) * | 2022-01-30 | 2022-04-26 | 国网浙江电动汽车服务有限公司 | 一种负荷聚集平台的可调容量预测方法及装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111740452A (zh) * | 2020-07-02 | 2020-10-02 | 华能澜沧江水电股份有限公司 | 一种单建议运行区水电机组有功功率控制方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108537405B (zh) * | 2018-03-07 | 2020-10-30 | 南方电网科学研究院有限责任公司 | 一种电力市场均衡优化方法、装置、设备及介质 |
| CN112163687B (zh) * | 2020-06-28 | 2024-03-22 | 国家电网有限公司西北分部 | 一种基于负荷聚合商的用户侧两阶段调峰响应方法 |
| CN112001744A (zh) * | 2020-07-28 | 2020-11-27 | 安徽电力交易中心有限公司 | 电力现货市场基于前景理论的发电商辅助报价***和方法 |
-
2021
- 2021-01-12 CN CN202110038800.2A patent/CN112766695B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111740452A (zh) * | 2020-07-02 | 2020-10-02 | 华能澜沧江水电股份有限公司 | 一种单建议运行区水电机组有功功率控制方法 |
Non-Patent Citations (1)
| Title |
|---|
| 风电并网后电力***可靠性评估和备用优化研究;余民等;《电力***保护与控制》(第12期);106-110 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112766695A (zh) | 2021-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110188950B (zh) | 基于多代理技术的虚拟电厂供电侧和需求侧优化调度建模方法 | |
| CN112054513A (zh) | 基于混合博弈的多微电网双层协调优化调度方法 | |
| CN112766695B (zh) | 一种负荷聚集主体参与下的主辅联合***的均衡运行方法 | |
| CN112348591B (zh) | 一种电动汽车的有序充电控制方法和装置 | |
| CN112186768B (zh) | Mg、la和dno共同参与的交直流配电网协同调度方法与*** | |
| CN114662759B (zh) | 多主体双层博弈的规模化电动汽车充放电优化调度方法 | |
| CN105631553A (zh) | 一种电动汽车充放电自动需求响应优化方法 | |
| CN114155103A (zh) | 基于区块链合作博弈的能量共享联盟灵活***易方法 | |
| CN110110895A (zh) | 电动汽车优化调度的方法及终端设备 | |
| CN115271438B (zh) | 一种可计及碳排放的多主体博弈协同调度方法及电子设备 | |
| CN113888204A (zh) | 多主体博弈的虚拟电厂容量优化配置法 | |
| WO2024092954A1 (zh) | 基于深度强化学习的电力***调控方法 | |
| CN110829424A (zh) | 一种考虑风电消纳的区域多微网配电网互动方法 | |
| CN114862460A (zh) | 一种基于主从博弈的电动汽车虚拟电厂运营商的定价方法 | |
| CN110556821B (zh) | 考虑交互功率控制和双边竞价交易的多微电网双层优化调度方法 | |
| Li et al. | Energy management method for microgrids based on improved Stackelberg game real-time pricing model | |
| CN114708071A (zh) | 考虑虚拟电厂对现货市场结算影响的双层竞价***及方法 | |
| CN114362139A (zh) | 一种以消纳新能源为目标的源-荷-储多边议价日前响应控制方法 | |
| CN113870030A (zh) | 基于改进纳什议价方法的多微网能源交易机制设计方法 | |
| CN111091265A (zh) | 一种考虑电动汽车参与平衡市场的出清模型及结算方法 | |
| CN111899122B (zh) | 一种基于储能控制的用户分散出清方法 | |
| Ye et al. | Equilibrium Analysis on the Main Energy and Ancillary Service Joint Market Considering Load Aggregators | |
| Wang et al. | Research on Electric Vehicle Aggregators Participating in Demand Response Pricing Strategy and Benefit Distribution | |
| CN118350878A (zh) | 电动汽车聚合商参与需求响应定价及收益分配方法 | |
| Huang et al. | A Novel Coordination Bidding Mechanism of Virtual Power Plant Based on Stackelberg Game |
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 |