Next Article in Journal
Literature Review on Thermodynamic and Kinetic Limitations of Thermal Decomposition of Methane
Next Article in Special Issue
Research on Electric Hydrogen Hybrid Storage Operation Strategy for Wind Power Fluctuation Suppression
Previous Article in Journal
Analysis of Electricity Supply and Demand Balance in Residential Microgrids Integrated with Micro-CAES in Northern Portugal
Previous Article in Special Issue
Research on Optimal Scheduling Strategy of Differentiated Resource Microgrid with Carbon Trading Mechanism Considering Uncertainty of Wind Power and Photovoltaic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bidding Strategy for the Alliance of Prosumer Aggregators in the Distribution Market

by
Chunyi Wang
1,
Jiawei Xing
2,
Yuejiao Wang
2,
Jing Xu
3,
Zhixin Fu
3,*,
Benjie Xu
3 and
Haoming Liu
3
1
State Grid Shandong Electric Power Company, Jinan 250001, China
2
State Grid Shandong Electric Power Research Institute, Jinan 250002, China
3
School of Electrical and Power Engineering, Hohai University, Nanjing 211100, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(19), 5006; https://doi.org/10.3390/en17195006
Submission received: 12 August 2024 / Revised: 24 September 2024 / Accepted: 30 September 2024 / Published: 8 October 2024

Abstract

:
Photovoltaic energy storage system (PV-ESS) prosumer aggregators are characterized by a large number but small scale in the distribution system and are not competitive enough to participate in market transactions. For this reason, a prosumer aggregator alliance is proposed to participate in the distribution market bidding strategy. Firstly, based on the framework for prosumer aggregator alliances participating in distribution market trading, a bilevel bidding model is constructed. The upper level represents the optimal decision-making model for the prosumer aggregators, while the lower level constitutes the distribution market-clearing model. Secondly, the additional benefits obtained by the alliance are distributed more fairly using the improved Shapley value based on the PV self-consumption rate. Given the problem that the traditional diagonalization algorithm (DA) has an excessive number of iterations when solving the game equilibrium problem of multiple subjects, the DA is improved by optimizing the initial value of the inputs. Finally, case studies are conducted based on the improved IEEE-33 bus distribution system to validate the feasibility and economic viability of the proposed strategy. The case study results show that forming cooperative alliances to participate in market bidding can significantly increase overall profits. The improved DA reduces the number of bids and computation time by 75% and 80%, respectively. Additionally, the improved Shapley value facilitates compensation for some of the aggregators.

1. Introduction

In recent years, there has been a notable surge in the advancement of distributed generation, represented by distributed photovoltaics (PVs) [1]. Numerous distributed PV aggregators are now participating in electricity spot market trading [2]. Due to their small individual capacities and fluctuating output, distributed PV aggregators face challenges in market competitiveness, compounded by significant curtailment issues [3]. However, energy storage systems (ESSs) enable time-shifted reuse of energy [4]. PV systems combined with ESSs can effectively reduce curtailment issues when participating in the market as a whole [5]. However, the number of distributed PV-ESSs is large but typically small in scale [6]. Compared to traditional units, their market competitiveness remains limited. These PV and ESS resources can be managed and scheduled by prosumer aggregators or virtual power plants. This approach allows for better resource utilization and increased overall revenue [7]. This paper primarily considers the agency of prosumer aggregators for these PV and ESS resources.
As new market entities, prosumer aggregators have attracted significant attention from researchers. Reference [8] proposes a trading strategy for prosumer aggregators based on P2P technology to participate in microgrid transactions. However, due to the small scale of prosumers and user demand, they can only engage in transactions within limited areas. Reference [9] proposes a mechanism based on a novel transactive energy framework. However, their research focuses only on the energy market. In contrast, reference [10] considers the production–consumption characteristics of PV-ESSs and proposes a two-stage bidding strategy for prosumers in the energy-reserve market. Reference [11] highlights the excellent frequency modulation (FM) capabilities of ESSs and presents a two-stage bidding framework for prosumers participating in the energy–FM auxiliary service market. However, references [8,9,10,11] treat prosumer aggregators as price takers, neglecting their influence on clearing prices. With the continuous improvement of auxiliary service market rules and the increasing number of prosumer aggregators, those with strong regulation capabilities are becoming competitive in the electricity market, and their bidding strategies can impact clearing outcomes. Therefore, references [12,13] propose a bidding strategy for prosumer aggregators participating in the market, with the market operator responsible for calculating and announcing the final clearing price. However, these studies consider prosumer aggregators as independent entities directly participating in market bidding, without accounting for the potential collaboration among storage stations of different aggregators to enhance overall returns. References [14,15], respectively, study the formation of alliances by combining wind ESSs and virtual power plants for market transactions. However, no studies in the literature have examined prosumer aggregators as members of an alliance yet. Therefore, this study is highly significant. Furthermore, after forming an alliance, fair profit distribution is crucial for maintaining cooperation stability. The Shapley value is a commonly used method for profit distribution in cooperative games [16]. However, in practical scenarios, reallocation based on the contributions and sacrifices of alliance members is necessary.
Multi-agent participation in electricity markets is often modeled as an Equilibrium Problem with Equilibrium Constraint (EPEC), where traditional approaches include intelligent algorithms [17] or the combination of KKT conditions with diagonalization algorithms (DAs) [18,19,20] for the solution. Reference [17] employs a prioritized experience replay deep deterministic policy gradient method to address this issue. However, this method is based on reinforcement learning, making the training process complex. Moreover, its convergence is influenced by numerous factors. As the number of market participants and constraints increases, the number of iterations required for such intelligent algorithms will increase significantly, leading to a dramatic rise in solution time, which makes them less suitable for practical problem-solving. Conversely, leveraging KKT conditions and strong duality transforms the EPEC into solving a Mixed-Integer Linear Programming (MILP) problem for multiple agents. Subsequently applying the DA effectively mitigates these issues. The core idea of the DA involves iteratively solving each market participant’s optimal decision model while keeping other participants’ decision variables fixed, iterating until strategies stabilize or reach a maximum iteration limit. References [18,19,20] specifically address strategic bidding by multiple market participants, utilizing the DA to find Nash equilibrium solutions. This methodical approach ensures robust convergence towards market equilibrium, maintaining stability in competitive electricity markets. However, significant discrepancies between the initial operational strategies chosen by each agent and their optimal strategies can result in excessive iteration. Moreover, the assumption of knowing other participants’ strategies during iterative self-strategy optimization does not reflect real-world scenarios. Therefore, this paper proposes an improved DA to address these issues, where each agent optimizes its bidding strategy by predicting other agents’ bidding strategies, thereby significantly reducing the solution time to just one bidding round.
In summary, a bidding strategy is developed for a coalition of prosumer aggregators to engage in distribution market bidding. The key contributions of the paper are as follows:
A framework for prosumer aggregator alliances to participate in electricity market trading is proposed. This framework uses an improved DA and KKT conditions to solve the EPEC problem. It ensures that the solution’s accuracy remains within an acceptable range while significantly reducing computation time.
Additional revenues earned by the alliances are allocated based on improved Shapley values. This approach effectively compensates the ESS for the losses incurred from participating in the FM market to support PV generation.

2. Trading Framework of Prosumer Aggregator Alliance Participation in Market

With the deepening reform of the electricity market, more and more new market players are involved in the market competition. Distributed PVs and distributed ESSs cannot directly participate in electricity market transactions due to their small volume, decentralized layout, and difficult regulation and management. Therefore, distributed PVs and distributed ESSs can be aggregated by aggregators in a certain region for unified scheduling. However, the capacity of distributed PVs and distributed ESSs represented by a single aggregator remains limited, with minimal market competitiveness. Therefore, forming alliances through cooperative gaming can leverage resource complementarity and enhance overall revenue. The distribution market trading framework based on cooperative gaming with the participation of multiple prosumer aggregators is shown in Figure 1, and this paper mainly considers the multi-market joint bidding of prosumer aggregators in the energy–FM market.
The specific market trading process is as follows:
  • Each market participant is required to verify unit parameters and relevant details concerning FM auxiliary service provision capabilities.
  • On declaration day, the electricity trading center releases information to all market participants about the next day’s active power demand, frequency modulation demand, bidding limits for power, and similar details. Subsequently, the operator organizes market participants to participate in the declaration process.
  • The trading period is set at 15 min, totaling 96 periods. Microturbine (MT) units declare bids based on their respective generation costs, while prosumer aggregators search for optimal alliances. These alliances then collectively submit unified “price-quantity” bids.
  • Considering the secure functioning of the distribution network and physical constraints, the objective is to meet day-ahead market load and FM requirements. Electricity trading centers strive to optimize day-ahead market clearing, with the primary objective of minimizing energy costs, thus determining the day-ahead clearing results.
  • The electricity trading center announces the final clearing result to all market members.

3. Optimal Decision Model for Upper-Level Prosumer Aggregator

The upper-level model is a trading decision model for prosumer aggregators, with the objective of maximizing aggregator profits. This model generates the optimal bidding strategies for each aggregator. The uncertainty of PV power generation is characterized through scenario analysis.

3.1. Objective Function of Prosumer Aggregator

The objective function of the aggregated prosumer aggregator is described in Equation (1), which includes participation in energy–FM market revenues and curtailment penalty costs.
max F k = s N s ρ s t = 1 T λ t , s en P k , t , s en + λ t , s cap P k , t , s cap + λ t , s mil P k , t , s mil C k , t , s PV
C k , t , s PV = v = 1 V k c PV P k , v , t , s PV , fore P k , v , t , s PV
P k , t , s en = v = 1 V k P k , v , t , s PVe + e = 1 E k P k , e , t , s dis e = 1 E k P k , e , t , s ch P k , t , s cap = v = 1 V k P k , v , t , s Pcap + e = 1 E k P k , e , t , s Ecap P k , t , s mil = v = 1 V k P k , v , t , s Pmil + e = 1 E k P k , e , t , s Emil
P k , v , t , s PV = P k , v , t , s PVe + P k , v , t , s Pcap
Here, ρ s is the probability of scenario s; P k , t , s en , P k , t , s cap , and P k , t , s mil are the amount of awarded electricity, the FM capacity, and mileage in the market by the aggregator k at time t under scenario s, respectively. λ t , s en , λ t , s cap , and λ t , s mil are the price of electricity, the FM capacity, and the mileage clearing price, obtained through the joint clearing of the lower-level market, respectively; C k , t , s PV is the costs incurred in the process of participating in market transactions for distributed PV; V k and E k are distributed PV and distributed ESS quantities represented by aggregator k, respectively; c PV is the penalty factor for PV abandonment; P k , v , t , s PV and P k , v , t , s PV , fore are the aggregators’ planned PV output and predicted PV output. Equation (3) represents the allocation of the aggregator k allocating the awarded quantities of charging, discharging, the FM capacity, and the FM mileage to each member, where P k , v , t , s PVe , P k , e , t , s dis , and P k , e , t , s ch represent the planned output of distributed PVs and the planned charging and discharging power of distributed ESSs by aggregator k; and P k , v , t , s Pcap , P k , v , t , s Pmil , P k , e , t , s Ecap , and P k , e , t , s Emil are aggregator k planned FM capacity and FM mileage output for distributed PVs and distributed ESSs in the FM market, respectively.

3.2. Constraints

3.2.1. Constraints of Prosumer Aggregators’ Declared Power

Constraints (5) limit the declared power range of the prosumer aggregator in the energy–FM market.
e = 1 E k P k , e , max ch p k , t , s en + p k , t , s cap v = 1 V k P k , v , t PV , fore + e = 1 E k P k , e , max dis 0 p k , t , s cap v = 1 V k P k , v , max Pcap + e = 1 E k P k , e , max Ecap 0 p k , t , s mil v = 1 V k s pv , v mc P k , v , max Pcap + e = 1 E k s ess , e mc P k , e , max Ecap
Here, p k , t , s en , p k , t , s cap , and p k , t , s mil are the power declared by aggregator k in the energy market and the FM capacity and mileage declared in the FM market, respectively; P k , e , max dis and P k , e , max ch are the maximum charging and discharging power of distributed ESSs; P k , v , max Pcap and P k , e , max Ecap are the maximum FM capacity for distributed PVs and distributed ESSs, respectively; s pv , v mc and s ess , e mc are mileage-to-capacity ratios for distributed PVs and distributed ESSs, respectively.

3.2.2. Constraints of Prosumer Aggregators’ Declared Price

Constraints (6) limit the declared price range of the prosumer aggregator in the energy–FM market.
0 b k , t , s en b k , max en 0 b k , t , s cap b k , max cap 0 b k , t , s mil b k , max mil
Here, b k , t , s en , b k , t , s cap , and b k , t , s mil are the price of electric energy, the FM capacity, and the mileage price declared by the aggregator, respectively; and b k , max en , b k , max cap , and b k , max mil are the maximum prices declared by aggregator k in the energy–FM market, respectively.

3.2.3. Constraint of Distributed PVs

Constraint (7) means that the planned PV output of the aggregator cannot exceed the forecasted output.
0 P k , v , t , s PV P k , v , t PV , fore

3.2.4. Constraints of Distributed ESSs

Constraints (8) describe the range of ESSs’ charging and discharging, along with the non-simultaneous charging and discharging characteristics.
0 P k , e , t , s dis α k , e , t , s dis P k , e , max dis 0 P k , e , t , s ch α k , e , t , s ch P k , e , max ch α k , e , t , s dis + α k , e , t , s ch 1
Here, P k , e , t , s dis and P k , e , t , s ch are the planned charging and discharging power of aggregator k for the distributed ESS, respectively; α k , e , t , s dis and α k , e , t , s ch are binary variables representing the discharge and charge states of the ESS; E k , e , t , s is the state of charge of the distributed ESS; E k , e , min and E k , e , max are the minimum and maximum capacities allowed by the ESS, respectively; and η k , e ch and η k , e dis are charge and discharge efficiency, respectively.

4. Integrated Clearing Model for Lower-Level Energy–FM Market

The lower-level model is a joint clearing model for the energy–FM market. Its objective function aims to minimize operational costs while determining clearing prices and the quantities allocated to various market participants.

4.1. Objective Function of Joint Market Clearing

The objective function of the energy–FM market clearing is described in Equation (10), which includes the operational costs of MT units and prosumer aggregators.
min F L = s N S ρ s t = 1 T m = 1 N DG ( b m , t , s Ge P m , t , s Ge + b m , t , s Gcap P m , t , s Gcap + b m , t , s Gmil P m , t , s Gmil ) + t = 1 T k = 1 N AG ( b k , t , s en P k , t , s en + b k , t , s cap P k , t , s cap + b k , t , s mil P k , t , s mil )
Here, N DG and N AG are the quantities of MT units and prosumer aggregators in the system, respectively; P m , t , s Ge , P m , t , s Gcap , and P m , t , s Gmil are the winning bid in the energy market and the winning bid in the FM market for FM capacity and FM mileage for MT m, respectively; b m , t , s Ge , b m , t , s Gcap , and b m , t , s Gmil are the declared generation price, FM capacity, and mileage price for MT m, respectively.

4.2. Constraints

4.2.1. Constraints of Network Operation

Constraint (10) ensures nodal power balance, constraint (10) limits branch transmission capacity, and constraint (12) establishes the balance node.
m Φ g ( i ) P m , t , s G + e Φ ess ( i ) ( P k , e , t , s dis P k , e , t , s ch ) + v Φ pv ( i ) P k , v , t , s PVe d Φ d ( i ) P d , t L = j Φ ( i ) B i j ( θ i , t θ j , t ) , i , t : λ 1 , i , t , s
P i j , max B i j ( θ i , t θ j , t ) P i j , max , i , q Ψ ( i ) : μ 1 , i j , t , s min , μ 1 , i j , t , s max
θ 1 , t = 1 : λ 1 , 1 , t , s
Here, Φ g ( i ) , Φ ess ( i ) , Φ pv ( i ) , and Φ d ( i ) represent the aggregation of MT units, distributed ESSs, distributed PVs, and loads set at node i, respectively; Φ ( i ) is the set of nodes connected to node i; Ψ ( i ) is the set of all branches connected to node i; P i , t L is the load of node i; B i j is the line conductance between nodes i and j; θ i , t and θ j , t are the phase angles of node i and node j, respectively; P i j , max represents the maximum transmission capacity of the branch (i,j); and λ 1 , i , t , s , μ 1 , i j , t , s min , and μ 1 , i j , t , s max are the pairwise multipliers corresponding to the equality constraint and the inequality constraint, respectively, where λ 1 , i , t , s is the clearing price of node i.

4.2.2. Constraints of FM Capacity and Mileage Demand

Constraints (13) ensure that the awarded quantities of various market participants in the FM market meet the system’s FM requirements.
m = 1 N DG P m , t , s Gcap + k = 1 N AG P k , t , s cap = R t SVS : λ 2 , t , s m = 1 N DG P m , t , s Gmil + k = 1 N AG P k , t , s mil = M t SYS : λ 3 , t , s
Here, R t SVS and M t SYS are the system’s FM capacity, and mileage requirements, respectively; λ 2 , t , s and λ 3 , t , s are dual variables corresponding to the constraints of the distribution system, namely FM capacity and mileage clearing price.

4.2.3. Constraints of Market Clearing

Constraints (14)–(15) limit the clearing power of MT units and aggregators in the energy–FM market, respectively.
P m , t , s Ge P m , min G + P m , t , s Gcap : μ 2 , m , t , s min P m , t , s Ge P m , max G P m , t , s Gcap : μ 2 , m , t , s max 0 P m , t , s Gcap p m , t , s Gcap : μ 3 , m , t , s min , μ 3 , m , t , s max 0 P m , t , s Gmil p m , t , s Gmil : μ 4 , m , t , s min , μ 4 , m , t , s max 0 P m , t , s Gmil s g , m mc P m , t , s Gcap : μ 5 , m , t , s min , μ 5 , m , t , s max
0 P k , t , s en p k , t , s en : μ 6 , k , t , s min , μ 6 , k , t , s max 0 P k , t , s cap p k , t , s cap : μ 7 , k , t , s min , μ 7 , k , t , s max 0 P k , t , s mil p k , t , s mil : μ 8 , k , t , s min , μ 8 , k , t , s max 0 P k , t , s mil v = 1 V k s pv , v mc P k , v , t , s Pcap + e = 1 E k s ess , e mc P k , e , t , s Ecap : μ 9 , k , t , s min , μ 9 , k , t , s max
Here, P m , max G and P m , min G denote the maximum and minimum output for MT m, respectively; s g , m mc is the FM mileage-to-capacity ratio for MT m; and μ o , m , t , s min , μ o , m , t , s max , o ( 2 , 10 ) are dual variables.

5. Formation of Cooperative Alliances and Benefit-Sharing Methods

As a single prosumer aggregator participates in market transactions, its capacity is small compared to other units, and it does not have an advantage in market competition, while the quantity of PVs and ESSs represented by some aggregators do not match, which can cause serious light abandonment problems in the transaction process. Therefore, aggregators are encouraged to cooperate to form an alliance, which increases the market competitiveness of aggregators and reduces the problem of light abandonment at the same time. How to fairly apportion the additional benefits generated after the formation of a cooperative alliance is the key to maintaining the stability of the alliance.

5.1. Improved Shapley Value Method Based on PV Self-Consumption Rate

The Shapley value method is the average marginal contribution of a particular alliance member in the alliance, and the Shapley value of aggregator k is assumed to be the cooperative game G = N , v , with N = n denoting the number of all aggregators:
φ k ( v ) = S N n ( | S | 1 ) ! ( n | S | ) ! n ! [ v ( S ) v ( S { k } ) ]
Here, φ k ( v ) is the revenue apportioned to aggregator k from the total revenue of alliance S; set N n is the n aggregators that can participate in the alliance S; v ( S ) is the total revenue of the union S; and v ( S { k } ) is the total revenue of the affiliate after removing the aggregator k.
This paper introduces the concept of a PV self-absorption rate, which represents the rate at which prosumer aggregators strategically charge ESSs during the bidding process in the market to accommodate the PV systems. Therefore, to encourage aggregator entities with greater distributed ESS capacity to join the alliance and enhance overall PV integration rates, a certain compensation for lost profits in the FM market will be provided. This involves adjusting the Shapley value method for profit distribution based on each aggregator’s self-consumption rate of solar energy. Let the PV self-consumption rate of aggregator k when bidding independently be ω k :
ω k = t = 1 T e = 1 E k P k , e , t PV - ESS t = 1 T v = 1 V k P k , v , t PV , fore
Here, P k , e , t PV - ESS is the energy absorbed from the PV system by the energy storage in aggregator k at time t.
The new weight m k of aggregator k is shown in Equation (18) and the gain M k of the improved aggregator k is shown in Equation (19).
m k = ω k k = 1 n ω k
M k = φ k ( v ) + ( m k 1 n ) k = 1 n φ k ( v )
Here, m k 1 n is the correction factor for each aggregator k. The sum of the correction factors for each aggregator is zero, so the total return after reallocation for each aggregator is the same as before correction, as demonstrated in Appendix A.

5.2. Model Transformation and Solution

The prosumer aggregator day-ahead market bidding model presented in the preceding section can be formulated as a multi-objective bilevel optimization problem.
max x S F S x S , y S N S   s . t .   P S x S 0 Q S x S , y = 0 y = argmin y   F L ( x , y ) s . t . p ( x , y ) 0   q ( x , y ) = 0
Here, F S · is the objective function of the aggregator alliance S; x S is a decision variable for the aggregator alliance S; y is a variable that is part of the market-clearing problem; and N S is a collection of aggregators. P S · and Q S · are inequality constraints and equality constraints for the aggregator alliance S; F L ( · ) is an objective function of the market-clearing problem; x is a decision variable for all aggregator alliances in the distribution system; and p ( · ) and q ( · ) are inequality constraints and equation constraints on the market-clearing problem.
The lower-level market-clearing problem is structured as a linear programming model. Therefore, it can be reformulated into its corresponding KKT conditions, which serve as constraints in the upper-level optimal decision model for prosumer aggregators.
max x S , y , λ , μ F S x S , x S , y , λ , μ S N S   s . t . P S x S 0 Q S x S , y = 0 { y , λ , μ } C KKT
Here, λ and μ are the dyadic variables of the equilibrium and inequality constraints in the market-clearing problem; x S is the decision-making of alliances other than alliance S; and C KKT is the set of constraints corresponding to the KKT condition for the market-clearing problem.
Using KKT conditions and strong duality theory, the bilevel model is transformed into a single-level model, thereby converting the multi-agent bilevel programming problem into an EPEC problem. The DA is a commonly used method for solving such problems. Furthermore, this paper improves upon traditional DAs. The improved DA first determines the optimal bidding strategy set x k ( 0 ) when all prosumer aggregators form a cooperative alliance, using it as the initial input, and then sequentially compute the bidding strategy sets x k ( f ) for each agent. For each aggregator, while fixing the bidding strategies of the other aggregators, the strategies are subsequently updated until convergence or reaching the maximum iteration count for each agent. This process decomposes the global problem into independent subproblems, allowing each agent to independently seek its optimal strategy and thereby achieve the global optimal solution. The specific steps for solving are detailed as follows:
Step 1: Obtain initial values. Obtain the bidding strategy when all aggregation merchants form an alliance as initial values x k ( 0 ) .
Step 2: Parameter settings. Set the maximum number of iterations f m a x and convergence criteria ε .
Step 3: Iterative solving. Let the iteration count be f = 1 , and input the optimized initial values x k ( 0 ) . Update each aggregator merchant’s solution by fixing the bidding strategy of the remaining aggregator merchants.
Step 4: Convergence criteria. If for all k, three consecutive iterations satisfy x k ( f ) x k ( f 1 ) ε , terminate and report the F k , p k , t , s e n , p k , t , s c a p , p k , t , s m i l . Else if f = f m a x , terminate and report non-convergence. Otherwise, f = f + 1 and return to Step 3.

6. Case Study

6.1. Example Setup

This case study is tested on the improved IEEE-33 node system, incorporating four MT units and three prosumer aggregators, as shown in Figure 2. The parameter information for the MT units and prosumer aggregators is provided in Appendix B [21]. Distributed PVs managed by Aggregator 1 are situated at nodes 3, 4, and 20; those managed by Aggregator 2 are at nodes 11 and 13; and those managed by Aggregator 3 are at nodes 9 and 28, with equal installed capacities. The PV output scenarios are set to five, with typical PV output scenarios shown in Figure 3 and their corresponding probabilities listed in Table 1. Assuming the total system FM requirement equals 10% of the total load, with a curtailment penalty cost of 0.6 CNY/kWh, the system daily load profile is depicted in Figure 4 [22,23]. To avoid vicious competition, given the offer range of the PV-ESS alliance, the electric energy offer is [0.1, 0.3] CNY/kWh, and the FM capacity and FM mileage offer range is [0, 0.02] CNY/kWh [24]. The convergence criterion for the solving process is set such that the difference between successive iterations is less than or equal to 1.5%. Convergence is considered achieved if this criterion is met for three or more consecutive iterations [25]. This paper uses Gurobi 10.0 in MATLAB R2021b to solve the problem.

6.2. Algorithm Comparison Analysis

The upper limit values of declared prices in each market were used as initial inputs for solving the traditional DA and compared with the improved DA proposed in this paper. Taking the case of three independent aggregators participating in market bidding as an example, a comparative analysis was conducted. The analysis focused on iteration count, solution time, and overall profit, as shown in Figure 5. The improved DA reduced the number of bidding attempts by three and significantly shortened the total solution time. Due to the sensitivity of the diagonalization algorithm to initial values, the iteration count increases substantially as the initial inputs deviate further from optimal results. Furthermore, an increase in the number of target entities leads to an exponential rise in iteration count. In terms of overall profits, the improved algorithm shows a decrease of 0.37% in total profits compared to the equilibrium solution after the first round of bidding. Combined with actual solving times, utilizing the improved algorithm for a single bidding round can significantly enhance solving efficiency. Therefore, subsequent case studies in this paper employ the improved DA, which only requires one round of bidding for solving.

6.3. PV-ESS Alliance Bidding Results

Based on the merge–split rule introduced in Section 5.1, the optimal alliance combination identified is 1 , 2 , 3 . Specifically, prosumer Aggregator 1 forms Alliance 1, while prosumer Aggregators 2 and 3 form Alliance 2. The bidding outcomes of each alliance in the energy–FM market are illustrated in Figure 6 and Figure 7.
Figure 6 depicts the bidding outcomes of various MT and PV-ESS alliances in the energy market. It can be observed that the two PV-ESS alliances account for approximately 50% of the total energy demand cumulatively across different periods, primarily concentrated during hours 20–77. Due to the higher solar generation during this period, the solar-storage alliances aim to minimize curtailment penalties. They achieve this by offering prices slightly below those of MT units to capture market share while maximizing overall revenue. MT units predominantly generate power during early morning and evening hours when distributed PVs are not producing power. MT unit 1 is cleared first due to having the lowest bid. When its output reaches the upper limit, units 2 and 3 with relatively lower bids supplement the generation capacity.
The winning bids of each market player in the FM market for each time period are shown in Figure 7. As can be seen from the figures, the FM capacity and mileage provided by the two PV-ESS alliances account for about 99% of the system demand, respectively. Since the distributed PVs and distributed ESSs in the alliance have a higher FM mileage-capacity factor and can provide more FM mileage than MT units when providing the same FM capacity, the PV-ESS alliance with stronger FM responsiveness is preferentially called upon to provide FM services in the FM auxiliary service market.
Taking Aggregator 2 as an example, a comparative analysis of the scheduling plans of the distributed PVs and distributed ESSs that it represents in the case of bidding independently and bidding after allying is conducted, and the results are shown in Figure 8 and Figure 9. The dispatch of distributed PVs and distributed ESSs by the aggregator increased by 6% and 10% during the 36–67 time period in the energy market, respectively. During midday hours, PV generation is high while load demand is low. Therefore, to maximize overall revenue, excess electricity generated by PV systems represented by Aggregator 2 is partially stored in distributed ESS 3 and ESS 4. Subsequently, it is sold during the nighttime peak load periods.
As shown in Figure 9a,b, it can be observed that when participating in the FM market, the aggregator’s dispatch of PVs and ESSs is reduced in almost all time periods. Simultaneously, distributed PV outputs remain concentrated during the midday period, whereas distributed ESSs exhibit an opposite trend, with their output times displaying complementary characteristics. In the energy market, the unit power price is significantly higher than in the FM market. Aggregators prioritize participation in the energy market to pursue greater profits, thereby increasing their involvement as much as possible. This inevitably reduces the capacity available for participation in the FM market. Additionally, aggregators primarily utilize dispatchable distributed ESSs to participate in the FM market, accounting for 83% of total ESS output. During peak PV generation periods, the primary task of distributed ESSs is to absorb excess solar power and opportunistically sell it. During periods of low PV generation, the main task shifts to frequency modulation.
Figure 10 depicts the cumulative curtailment of distributed PV systems for Aggregator 2 and Aggregator 3 across different time periods, under cooperative and non-cooperative scenarios. The total curtailment without cooperation amounts to 11.13 MWh, whereas with cooperation in forming an alliance, the total curtailment decreases to 7.99 MWh, a reduction of 28.22% in curtailed energy.

6.4. Benefits of the PV-ESS Alliance

The PV-ESS alliance currently participates primarily in the energy market and FM market. To illustrate the necessity of the PV-ESS alliance’s participation in the energy–FM market, comparisons were also made regarding the daily revenue performance of the PV-ESS alliance in different market cases, as detailed in Table 2.
The comparative results are shown in Table 3. It is clear from the table that the alliance generates higher revenue by participating in the combined market. Compared to participating in only one market, the total revenue increases by CNY 8361.55 and CNY 16,556.78, respectively. Therefore, actively participating in the combined energy–FM market is beneficial for the alliance to obtain higher revenue, encouraging their engagement in market competition and promoting the development of PV-ESSs.
In existing research, there are cases where only distributed ESSs are considered as individual market entities participating in bidding. Therefore, this paper considers all ESSs as being represented by a single ESS aggregator. It compares and analyzes the revenues of the ESS aggregator participating in market bidding independently under this scenario. The comparison results are presented in Table 4. The table indicates that the total revenue of the integrated prosumer aggregator exceeds the combined revenue of the ESS aggregator and the PV generator by CNY 1682.06. When the ESS aggregator participates in the market independently, it competes with the distributed PV generator, which results in both ESS and PV resources being underutilized. The aggregator can minimize resource wastage by conducting unified scheduling and management of ESSs and PVs.
Table 5 presents the revenue outcomes for each aggregator under cooperative and non-cooperative conditions. According to the table, forming alliances through mutual cooperation increases the total revenue for prosumer aggregators by CNY 887.72. Moreover, to highlight the benefits of cooperation, the revenue includes not only market participation earnings but also curtailment penalty costs. When allocating additional revenue gained by the alliance using Shapley values, Aggregator 2 and Aggregator 3 see revenue increases of CNY 495.65 and CNY 402.07, respectively. Utilizing the improved Shapley values for revenue distribution, the new weights for Aggregator 2 and Aggregator 3 are 0.47 and 0.53, resulting in revenue increases of CNY 239.14 and CNY 648.58, respectively. Under both allocation methods, aggregators experience enhanced revenue, with Aggregator 3 showing a higher growth rate in revenue when using the improved Shapley values. As analyzed earlier, the ESS in Aggregator 3 sacrifices some revenue from FM markets for overall benefit by accommodating more PV power. Compared to non-cooperation, this cooperation results in a loss of CNY 150.62 in frequency modulation revenue. The use of improved Shapley values compensates Aggregator 3, encouraging more aggregators to join the alliance through fairer revenue distribution.

7. Conclusions

This paper addresses the bidding strategy problem of prosumer aggregators participating in multi-species electricity markets, proposing a bilevel bidding model for prosumer aggregators considering uncertainty. A revenue allocation method based on the improved Shapley value of the PV self-consumption rate is proposed and solved using the improved DA. The following conclusions can be obtained from the case study:
  • Compared to traditional algorithms, the improved DA proposed in this paper achieves market equilibrium solutions more rapidly. Considering the actual computation time, using this improved DA approach requires only one round of bidding, which greatly boosts solution efficiency.
  • Taking Aggregator 2 as an example, after forming the optimal alliance, the scheduling volume of distributed PVs increases by 6% and distributed ESSs by 10%. Furthermore, Alliance 2 demonstrates a 28% reduction in total curtailed solar energy, leading to a more efficient utilization of all resource types.
  • After forming the optimal alliance, total revenue increased by 7.3%. Considering the structural differences in the composition of prosumer aggregators, an improved Shapley value based on PV self-consumption rates was proposed. This approach allows for compensating aggregators that handle a higher proportion of distributed ESSs, thereby achieving a more equitable distribution of benefits.
In the subsequent work, further consideration will be given to the participation of prosumer aggregators in the operation of the real-time market within cooperative alliances. Additionally, a bidding strategy for the participation of prosumer aggregators in the day-ahead–real-time market will be proposed, taking into account the strong coupling between day-ahead and real-time markets. Additionally, when studying the real-time market, it is important to consider the impacts of special situations, such as the forced disconnection of photovoltaic sources under adverse operating conditions.

Author Contributions

Conceptualization, C.W., Z.F. and H.L.; methodology, Y.W., J.X. (Jing Xu) and B.X.; software, J.X. (Jing Xu) and B.X.; formal analysis, J.X. (Jing Xu); investigation, J.X. (Jiawei Xing); resources, J.X. (Jiawei Xing) and Y.W.; data curation, J.X. (Jiawei Xing), Y.W. and B.X.; writing—original draft, J.X. (Jing Xu); writing—review and editing, C.W., J.X. (Jiawei Xing), Y.W.; supervision, C.W. and Z.F.; and H.L.; project administration, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Project of State Grid Corporation of China: “The Research and Application of Distributed Photovoltaic and Energy Storage Aggregation Regulation Technology” (No. 5419-202316463A-3-2-ZN).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

Author Chunyi Wang was employed by the company State Grid Shandong Electric Power Company. Authors Jiawei Xing and Yuejiao Wang were employed by the company State Grid Shandong Electric Power Research Institute. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

The abbreviations used in this manuscript are as follows:
PVphotovoltaic
ESSenergy storage system
DAdiagonalization algorithm
FMfrequency modulation
EPECEquilibrium Problem with Equilibrium Constraint
MILPMixed-Integer Linear Programming

Appendix A

It has been shown through previous studies that when an ESS participates in the electric energy–FM market, most of its revenue comes from the participating FM market. It is assumed that the aggregator k1 acts as an agent for a larger number of distributed PV with a smaller number, while the aggregator k2 acts as an agent managing a larger number of distributed PV systems alongside a smaller number of energy storage units, and a cooperative game is played to ally the aggregators k1 and k2. Aggregator k1 allocates a portion of its energy storage capacity originally intended for FM market participation to store excess energy from distributed PV systems within aggregator k2 during peak solar generation periods. This ESS is subsequently sold during nighttime hours, generating additional revenue for the alliance in the energy market. Using the traditional Shapley value method for revenue allocation only considers the average marginal contribution of aggregators k1 and k2 and ignores the sacrifices made by aggregator k1 in the alliance. Therefore, this paper proposes an improved Shapley value.
The proof regarding the sum of correction factors for all aggregators being zero is as follows:
k = 1 n ( m k 1 n ) = k = 1 n ω k k = 1 n ω k k = 1 n 1 n = k = 1 n ω k k = 1 n ω k 1 = 0 k = 1 n M k = k = 1 n φ k ( v ) + ( m k 1 n ) k = 1 n φ k ( v ) = k = 1 n φ k ( v ) + k = 1 n ( m k 1 n ) k = 1 n φ k ( v ) = k = 1 n φ k ( v ) + k = 1 n ( m k 1 n ) k = 1 n φ k ( v ) = k = 1 n φ k ( v )

Appendix B

Table A1. The parameter information of MT units.
Table A1. The parameter information of MT units.
MT UnitsG1G2G3G4
Rated power/kW4000400035003000
Power generation quotation/CNY/kWh0.1830.1900.2080.218
Frequency modulation capacity quotation/CNY/kWh0.0120.0120.0160.013
Frequency modulation mileage quotation/CNY/kWh0.0160.0150.0150.017
Mileage and capacity ratio1210712
Table A2. The parameter information of distributed ESSs.
Table A2. The parameter information of distributed ESSs.
Aggregators123
ESS IDESS1ESS2ESS3ESS4ESS5ESS6ESS7
Charge/Discharge Power/kW400200900900800800800
Capacity/kWh200080036003600320032003200
Maximum Storage Capacity0.90.90.90.90.90.90.9
Minimum Storage Capacity0.10.10.10.10.10.10.1
Charge/Discharge Efficiency0.90.90.90.90.90.90.9
Energy Density Ratio15151515151515
Table A3. The parameter information of distributed PVs.
Table A3. The parameter information of distributed PVs.
Aggregators123
PV IDPV1PV2PV3PV4PV5PV6PV7
Capacity/kW2000200020002000200020002000
Energy Density Ratio15151515151515

References

  1. Yuan, L.; Yongxiu, H.; Shanzi, L.; Jinghan, Z.; Bo, L. Impact of grid connection cost channeling mechanisms on the development of distributed photovoltaic: The case of China. Energy Policy 2024, 187, 114041. [Google Scholar]
  2. Qiangang, J.; Yiyan, L.; Zheng, Y.; Sijie, C. Reactive Power Market Design for Distribution Networks With High Photovoltaic Penetration. IEEE Trans. Smart Grid 2023, 14, 1642–1652. [Google Scholar]
  3. Prado, J.C.D.; Wei, Q. A Stochastic Bilevel Model for an Electricity Retailer in a Liberalized Distributed Renewable Energy Market. IEEE Trans. Sustain. Energy 2020, 11, 2803–2812. [Google Scholar] [CrossRef]
  4. Li, X.; Yang, L.; Xu, Y.; Luo, X.; Yang, X.; Fang, J.; Lu, Y. Research on the Coordinated Trading Mechanism of Demand-Side Resources and Shared Energy Storage Based on a System Optimization Model. Energies 2024, 17, 3378. [Google Scholar] [CrossRef]
  5. Makowski, R.; Marcin, J.; Pawlicki, W. Operation of the Hybrid Photovoltaic-Battery System on the Electricity Market—Simulation, Real-Time Tests and Cost Analysis. Energies 2020, 13, 1402. [Google Scholar] [CrossRef]
  6. Li, H.; Jie, Z. Energy Trading in Local Electricity Markets With Behind-the-Meter Solar and Energy Storage. IEEE Trans. Energy Mark. 2023, 1, 107–117. [Google Scholar]
  7. Guoqiang, S.; Sichen, S.; Sheng, C.; Yizhou, Z.; Zhinong, W. Bidding strategy for a prosumer aggregator with stochastic renewable energy production in energy and reserve markets. Renew. Energy 2022, 191, 278–290. [Google Scholar]
  8. Zhenyuan, Z.; Haoyue, T.; Jie, R.; Qi, H.; Wei-Jen, L. Strategic Prosumers-Based Peer-to-Peer Energy Market Design for Community Microgrids. IEEE Trans. Ind. Appl. 2021, 57, 2048–2057. [Google Scholar]
  9. Peng, H.; Guangya, Y.; Junjie, H.; Philip, J.D.; Yusheng, X. A Distributed Transactive Energy Mechanism for Integrating PV and Storage Prosumers in Market Operation. Engineering 2022, 12, 171–182. [Google Scholar]
  10. Amaia, G.G.; Andoni, S.I.; Haizea, G.; Aitor, M.; Pablo, E. Annual Optimized Bidding and Operation Strategy in Energy and Secondary Reserve Markets for Solar Plants With Storage Systems. IEEE Trans. Power Syst. 2018, 34, 5115–5124. [Google Scholar]
  11. Bian, J.; Song, Y.; Ding, C.; Cheng, J.; Li, S.; Li, G. Optimal Bidding Strategy for PV and BESSs in Joint Energy and Frequency Regulation Markets Considering Carbon Reduction Benefits. J. Mod. Power Syst. Clean Energy 2024, 12, 427–439. [Google Scholar] [CrossRef]
  12. Ankun, Y.; Chaorui, Z.; Ying, J.A.Z. Optimal Bidding Strategy of Prosumers in Distribution-Level Energy Markets. IEEE Trans. Power Syst. 2020, 35, 1695–1706. [Google Scholar]
  13. José, I.; Filipe, S.; Manuel, M. Optimal bidding strategy for an aggregator of prosumers in energy and secondary reserve markets. Appl. Energy 2019, 238, 1361–1372. [Google Scholar]
  14. Li, P.; Wang, Z.; Jin, J. Market Impact of Wind-Energy Storage Alliance Strategic Bidding Under Uncertainty. IEEE Access 2021, 9, 156537–156547. [Google Scholar] [CrossRef]
  15. Zhu, Z.; Chan, K.W.; Xia, S.; Bu, S. Optimal Bi-Level Bidding and Dispatching Strategy Between Active Distribution Network and Virtual Alliances Using Distributed Robust Multi-Agent Deep Reinforcement Learning. IEEE Trans. Smart Grid 2022, 13, 2833–2843. [Google Scholar] [CrossRef]
  16. Pedrero, R.A.; Pisciella, P.; Granado, P.C.D. Fair investment strategies in large energy communities: A scalable Shapley value approach. Energy 2024, 295, 131033. [Google Scholar] [CrossRef]
  17. Yujian, Y.; Dawei, Q.; Mingyang, S.; Dimitrios, P.; Goran, S. Deep Reinforcement Learning for Strategic Bidding in Electricity Markets. IEEE Trans. Smart Grid 2020, 11, 1343–1355. [Google Scholar]
  18. Wang, C.; Wei, W.; Wang, J.; Liu, F.; Mei, S. Strategic Offering and Equilibrium in Coupled Gas and Electricity Markets. IEEE Trans. Power Syst. 2018, 33, 290–306. [Google Scholar] [CrossRef]
  19. Sheng, C.; Antonio, J.C.; Ramteen, S.; Zhinong, W. Equilibria in Electricity and Natural Gas Markets With Strategic Offers and Bids. IEEE Trans. Power Syst. 2020, 35, 1956–1966. [Google Scholar]
  20. Houhe, C.; Linbo, F.; Linquan, B.; Tao, J.; Yaosuo, X.; Rufeng, Z.; Badrul, C. Distribution Market-Clearing and Pricing Considering Coordination of DSOs and ISO: An EPEC Approach. IEEE Trans. Smart Grid 2021, 12, 3150–3162. [Google Scholar]
  21. Yang, L.; Bin, W.; Zhen, Y.; Jiazheng, L.; Guoqing, L. Optimal Scheduling of Integrated Demand Response-Enabled Community-Integrated Energy Systems in Uncertain Environments. IEEE Trans. Ind. Appl. 2022, 58, 2640–2651. [Google Scholar]
  22. Alex, D.P.; Panagiotis, E.A. Performance-Based Pricing of Frequency Regulation in Electricity Markets. IEEE Trans. Power Syst. 2014, 29, 441–449. [Google Scholar]
  23. Navid, A.R.; Krishnamurthy, D.; Niangjun, C.; Adam, W. Opportunities for Price Manipulation by Aggregators in Electricity Markets. IEEE Trans. Smart Grid 2018, 9, 5687–5698. [Google Scholar]
  24. Jing, X.; Zhixin, F.; Chunyi, W.; Haoming, L. The Trading Strategy of Distributed Energy Storage Participating in Energy and Frequency Modulation Market. In Proceedings of the 2024 7th International Conference on Energy, Electrical and Power Engineering (CEEPE), Yangzhou, China, 26–28 April 2024; pp. 1694–1700. [Google Scholar]
  25. Ting, D.; Wei, Q. Finding Equilibria in the Pool-Based Electricity Market With Strategic Wind Power Producers and Network Constraints. IEEE Trans. Power Syst. 2017, 32, 389–399. [Google Scholar]
Figure 1. Cooperative game-based distribution market trading framework with the participation of prosumer aggregator alliances.
Figure 1. Cooperative game-based distribution market trading framework with the participation of prosumer aggregator alliances.
Energies 17 05006 g001
Figure 2. Schematic diagram of improved IEEE-33 node system topology.
Figure 2. Schematic diagram of improved IEEE-33 node system topology.
Energies 17 05006 g002
Figure 3. Typical set of scenarios for PV output.
Figure 3. Typical set of scenarios for PV output.
Energies 17 05006 g003
Figure 4. System daily load profile.
Figure 4. System daily load profile.
Energies 17 05006 g004
Figure 5. Algorithm performance comparison results.
Figure 5. Algorithm performance comparison results.
Energies 17 05006 g005
Figure 6. Winning bids in the energy market for each subject.
Figure 6. Winning bids in the energy market for each subject.
Energies 17 05006 g006
Figure 7. Winning bids in the FM market for each subject. (a) FM capacity. (b) FM mileage.
Figure 7. Winning bids in the FM market for each subject. (a) FM capacity. (b) FM mileage.
Energies 17 05006 g007
Figure 8. Prosumer Aggregator 2 scheduling plan in the energy market. (a) ESS. (b) PV.
Figure 8. Prosumer Aggregator 2 scheduling plan in the energy market. (a) ESS. (b) PV.
Energies 17 05006 g008
Figure 9. Prosumer Aggregator 2 scheduling plan in the FM market. (a) FM capacity (b) FM mileage.
Figure 9. Prosumer Aggregator 2 scheduling plan in the FM market. (a) FM capacity (b) FM mileage.
Energies 17 05006 g009
Figure 10. The curtailed solar energy amount for cooperation and non-cooperation.
Figure 10. The curtailed solar energy amount for cooperation and non-cooperation.
Energies 17 05006 g010
Table 1. Probabilities corresponding to each scenario.
Table 1. Probabilities corresponding to each scenario.
Scenario12345
Probability0.190.210.210.220.17
Table 2. Comparison of cases for prosumer aggregator participation in various markets.
Table 2. Comparison of cases for prosumer aggregator participation in various markets.
CaseEnergy MarketFM Ancillary Service Market
Case 1YESNO
Case 2NOYES
Case 3YESYES
Table 3. Benefits of the alliance when participating in different markets.
Table 3. Benefits of the alliance when participating in different markets.
Benefits/CNYCase 1Case 2Case 3
Alliance 14105.991941.927681.34
Alliance 28650.392619.2313,436.58
All Benefits12,756.384561.1521,117.93
Table 4. Benefits of ESS and PVs managed by different entities.
Table 4. Benefits of ESS and PVs managed by different entities.
Benefits/CNYESSPVPV-ESS
Energy Benefits802.7112,036.5814,584.24
FM Capacity Benefits289.74149.47464.04
FM Mileage Benefits4062.112095.266069.65
All Benefits5154.5614,281.3121,117.93
Table 5. Benefits in case of cooperation and non-cooperation.
Table 5. Benefits in case of cooperation and non-cooperation.
SubjectBenefits/CNYBenefit Allocation/CNY
Non-CooperationCooperationShapley ValueImproved Shapley Value
Aggregator 14464.834464.83--
Aggregator 23776.718550.394272.364015.85
Aggregator 33885.964288.034534.54
All Benefits12,127.5013,025.22--
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.; Xing, J.; Wang, Y.; Xu, J.; Fu, Z.; Xu, B.; Liu, H. Bidding Strategy for the Alliance of Prosumer Aggregators in the Distribution Market. Energies 2024, 17, 5006. https://doi.org/10.3390/en17195006

AMA Style

Wang C, Xing J, Wang Y, Xu J, Fu Z, Xu B, Liu H. Bidding Strategy for the Alliance of Prosumer Aggregators in the Distribution Market. Energies. 2024; 17(19):5006. https://doi.org/10.3390/en17195006

Chicago/Turabian Style

Wang, Chunyi, Jiawei Xing, Yuejiao Wang, Jing Xu, Zhixin Fu, Benjie Xu, and Haoming Liu. 2024. "Bidding Strategy for the Alliance of Prosumer Aggregators in the Distribution Market" Energies 17, no. 19: 5006. https://doi.org/10.3390/en17195006

APA Style

Wang, C., Xing, J., Wang, Y., Xu, J., Fu, Z., Xu, B., & Liu, H. (2024). Bidding Strategy for the Alliance of Prosumer Aggregators in the Distribution Market. Energies, 17(19), 5006. https://doi.org/10.3390/en17195006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop