Next Article in Journal
Energy Consumption and Fume Analysis: A Comparative Analysis of the Blasting Technique and Mechanical Excavation in a Polish Gypsum Open-Pit Mine
Previous Article in Journal
Incentive Determination for Demand Response Considering Internal Rate of Return
Previous Article in Special Issue
A Fast Converter Synchronization Speed Method Based on the Frequency Difference on Both Sides of the Grid Connection Point
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nash Bargaining-Based Coordinated Frequency-Constrained Dispatch for Distribution Networks and Microgrids

1
State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300401, China
2
State Grid Tianjin Power Dongli Power Supply Branch, Tianjin 300300, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(22), 5661; https://doi.org/10.3390/en17225661
Submission received: 19 September 2024 / Revised: 31 October 2024 / Accepted: 11 November 2024 / Published: 13 November 2024

Abstract

:
As the penetration of distributed renewable energy continues to increase in distribution networks, the traditional scheduling model that the inertia and primary frequency support of distribution networks are completely dependent on the transmission grid will place enormous regulatory pressure on the transmission grid and hinder the active regulation capabilities of distribution networks. To address this issue, this paper proposes a coordinated optimization method for distribution networks and microgrid clusters considering frequency constraints. First, the confidence interval of disturbances was determined based on historical forecast deviation data. On this basis, a second-order cone programming model for distribution networks with embedded frequency security constraints was established. Then, microgrids performed economic dispatch considering the reserves requirement to provide inertia and primary frequency support, and a stochastic optimization model with conditional value-at-risk was built to address uncertainties. Finally, a cooperative game between the distribution network and microgrid clusters was established, determining the reserve capacity provided by each microgrid and the corresponding prices through Nash bargaining. The model was further transformed into two sub-problems, which were solved in a distributed manner using the ADMM algorithm. The effectiveness of the proposed method in enhancing the operational security and economic efficiency of the distribution networks is validated through simulation analysis.

1. Introduction

With the transition, the installed capacity of wind and solar power is continuously growing, gradually replacing the generation capacity of traditional coal-fired power plants [1,2]. Due to the limited ability of wind turbines and solar PV to provide inertia and primary frequency response (PFR), the inertia and frequency regulation resources of the power system are becoming increasingly constrained, posing severe challenges to the frequency security of the power system [3,4,5,6].
In response to the insufficiency of inertia and PFR capabilities in the power system, several studies have introduced frequency security constraints into the optimal dispatch models [7,8,9,10,11,12,13,14]. Reference [11] introduced the concept of frequency security margin into the dispatch model, modeled the maximum frequency deviation after a disturbance using a piecewise linear method, and incorporated it into the optimal dispatch model to optimize the base power and active reserves. Reference [12] used an analytical method to model the frequency security constraints for energy storage participation in frequency response, explored the value of inertia and PFR, and incorporated frequency response resources into the integrated market for clearing. Reference [13] proposed a convex hull approximation method to linearize the maximum frequency deviation constraint, obtained the optimal combination of inertia and droop coefficients, and provided inertia support and PFR through reserving backup capacity in energy storage. However, these studies have only focused on the optimal dispatch of transmission system, neglecting the frequency response capabilities of resources in the distribution side [15,16]. Reference [17] introduced frequency security constraints based on transmission and distribution coordination, where both the distribution and transmission networks provided bidirectional frequency support. Nonetheless, the coordination between the distribution networks and microgrids was not included.
With the development of the electricity market, numerous studies have focused on the coordination between distribution networks and microgrids to enhance the overall stability and efficiency of the power system. Reference [18] established a cooperative game model for multi-microgrids on the distribution side, enhancing energy complementarity among microgrids. This approach improved the overall energy efficiency by leveraging the diverse generation resources of different microgrids, thereby achieving optimized power sharing and reducing operational costs. Similarly, reference [19] proposed a coordination model between the distribution network and microgrids, where the distribution network guided the operational strategies of microgrids through electricity price signals. In this model, the distribution network set the price to incentivize microgrid participation in the electricity market and optimize power flow. However, these game-theoretic models mainly focus on power exchange optimization, neglecting frequency security and the microgrids’ potential to provide inertia support and PFR. This oversight limits the applicability of these models, especially in scenarios where maintaining frequency stability is crucial, such as during grid disturbances or a high penetration of renewable energy sources.
The ability of microgrids to provide inertia support and PFR has become a significant research focus in recent years, as these capabilities are vital for maintaining frequency stability in low-inertia power systems. Reference [20] integrated a dynamic frequency security model into the microgrid scheduling framework, incorporating the use of storage units within microgrids to provide inertia and PFR. This approach helped mitigate frequency fluctuations caused by variable renewable energy sources like PV and sudden load changes. By leveraging the storage systems’ fast response, microgrids could actively support the grid’s frequency, thus enhancing overall stability. On the other hand, Reference [21] addressed frequency disturbances in microgrid operations by implementing load shedding strategies during islanded operation. However, this approach primarily focused on ensuring frequency stability within isolated microgrids rather than considering the support to distribution networks during interconnected operations.
In summary, current studies have not fully considered frequency security constraints in the optimal dispatch of distribution networks. Although some studies have explored the frequency response capabilities of microgrids, they have not leveraged microgrids to provide inertia support and PFR to the distribution networks. Additionally, how to incentivize microgrids to actively provide inertia and PFR services remains an issue that requires further investigation. To excavate the capability of microgrids to provide inertia and PFR support, thereby reducing the reserve costs of distribution networks, this paper constructs a coordinated frequency-constrained dispatch model for distribution networks and microgrids, exploring a novel interactive mode between microgrids and distribution networks towards future demand for stronger frequency regulation capability of distribution networks. The main contributions are as follows:
(1) A second-order cone programming model for distribution networks with embedded frequency security constraints was established. Based on an equivalent frequency response model, constraints such as maximum frequency rate of change, maximum frequency deviation, and quasi-steady-state frequency were established. From this, the optimal combination of inertia and droop coefficients that met the safe operation requirements of the distribution network was determined, thereby defining the reserve requirements of the distribution network.
(2) A bi-level coordinated optimization model for the distribution network and microgrids was established, where the distribution network and microgrids negotiated reserve capacity and reserve prices through cooperative bargaining, thereby reducing the reserve costs for inertia and PFR resources of the distribution network and improving the revenue of microgrids. The alternating direction method of multipliers (ADMM) [22] was used to solve the problem in a distributed manner, achieving a collective optimum for the distribution network and microgrids.

2. The Interaction Framework Between Distribution Networks and Microgrids

The interaction framework between distribution networks and microgrids is shown in Figure 1, which is essentially a bi-level optimization dispatch model. Unlike the traditional model that only involves power exchange [23], this paper incorporated frequency security constraints into the distribution network model to address disturbances in the distribution network. On this basis, to explore the inertia support and PFR capabilities of microgrids, the distribution network and microgrids engage in cooperative bargaining to negotiate the active reserve capacity required for providing frequency response services and its pricing. The negotiation process is as follows:
(1)
The distribution network calculates the required inertia, droop coefficients, and corresponding reserve capacity for each period based on the optimization model considering frequency security constraints.
(2)
The distribution network interacts with microgrids for reserve energy exchange through Nash bargaining. In each round of the game, the distribution network adjusts and communicate their reserve demand based on the current supply–demand situation. The microgrids then adjust their reserve supply and electricity trading volumes for each period according to the distribution network’s demand and provide feedback. Through multiple iterations, the distribution network and microgrids gradually optimize their strategies until the reserve demand matches the reserve supply.
(3)
Following the same interaction process as in step 2, the distribution network and microgrids engage in reserve price negotiations until the price offered by the distribution network matches the price required by the microgrids.
In the model constructed in this paper, the microgrid includes distributed PV systems, wind turbines, battery energy storage systems, as well as electric boilers and absorption chillers, allowing for integrated energy utilization. By optimizing its energy interaction scheme with the distribution network, the microgrid meets various energy demands and enhances its operational economics by providing reserve services.

3. Mathematical Model

3.1. Optimization Model for Distribution Networks

3.1.1. Frequency Security Modeling

The transmission network can be modeled as an equivalent generator model [24]. The microgrid provides inertia and PFR to distribution networks using virtual synchronous generator technology, resulting in the dynamic frequency response model of distribution networks, shown in (1) [14].
2 ( H g + i m H i ) d Δ f ( t ) d t + D 0 Δ f ( t ) = Δ P g + i m Δ P i Δ P
where H g represents the inertia provided by the transmission network, H i represents the virtual inertia provided by microgrid i , and D 0 represents the system damping. Δ P g and Δ P i are the response power of the transmission network and the microgrid i , Δ P is the disturbance in distribution networks, and Δ f is the frequency deviation.
The frequency response model of the distribution network is constructed using an equivalent generator model, and its control block diagram is shown in Figure 2 [11].
In Figure 2, F H and T R represent the high-pressure turbine power generation proportion and the generator reheating time, respectively. H total is the system equivalent inertia, and D total is the system equivalent droop coefficient, which are calculated as [25]:
H total = H g P g P d + i m H i P i P d
D total = D g P g P d + i m D i P i P d
where P g is the maximum output of the transmission network, and D g is the droop coefficient it provides. P d is the maximum load demand of distribution networks. P i is the maximum output of the microgrid, and D i is the droop coefficient provided by microgrid i.
Using the equivalent system frequency response model, the time-domain model of the frequency response can be obtained [13]:
Δ f = Δ P D 0 + D total [ 1 + α e ξ ω n t sin ( ω r t + ϕ ) ]
where
ω n 2 = D total + D 0 2 H total T R ξ = D 0 T R + 2 H total + D total F H T R 2 ( D 0 + D total ) ω r = ω n 1 ξ 2 α = 1 2 T R ξ ω n + T R 2 ω n 2 1 ξ 2 ϕ = arctan ( ω r T R 1 ξ ω n T R ) arctan ( 1 ξ 2 ξ )
The natural frequency ω n , damping frequency ω r , damping ratio ξ , and coefficients α and ϕ are functions of H total and D total , as defined in (5). The physical meanings of these coefficients can be found in reference [26].
Through the time-domain model, the dynamic frequency indicators of the frequency response can be transformed into the following constraints [11]:
2 H total f rate ¯ Δ P
( D 0 + D total ) f ss ¯ Δ P
D total + D 0 1 + 1 ξ 2 α e ξ ω n t nadir Δ f nadir ¯ Δ P
Formula (6) represents the constraint of the maximum rate of frequency change, (7) represents the quasi-steady frequency constraint, and (8) represents the maximum frequency deviation constraint. f rate ¯ , f ss ¯ and Δ f nadir ¯ represent the maximum allowable rate of change of frequency (RoCoF), quasi-steady-state frequency and frequency deviation of the system, respectively. t nadir denotes the time to reach the maximum frequency deviation, which can be calculated by:
t nadir = 1 ω r arctan ( ω r T R ξ ω r T R 1 )
Let F ( H total , D total ) = D total + D 0 1 + 1 ξ 2 α e ξ ω n t nadir . Through piecewise linearization [27], a set of linear constraints on maximum frequency deviation can be obtained.
( β h C + β h H H total + β h D D total ) × Δ f nadir ¯ Δ P , h
where β h C , β h H , and β h D represent a set of hyperplanes approximating F ( H total , D total ) . For the specific derivation process, please refer to reference [11].

3.1.2. Optimization Model of Distribution Networks

(1)
Objective Function
The cost of distribution networks is divided into two parts:
C DN = min C P , DN + C R , DN
where C P , DN and C R , DN are the cost of electricity interaction and reserve service from the transmission network and microgrids, respectively. The electricity interaction cost can be expressed by (12):
C P , DN = t = 1 N T P g , t c t buy Δ t + i m ( t = 1 N T P i , t buy c t buy t = 1 N T P i . t sell c t sell ) Δ t
C R , DN = t = 1 N T ( R g , t × c t R + i m R i , t × c i , t R ) Δ t
where P g , t is the electricity purchased from the transmission network at period t , c t buy is the electricity purchase price; P i , t buy and P i . t sell are the amounts of electricity purchased from and sold to microgrid i , respectively, and c t sell denotes the electricity sale price. R g , t and R i , t are the reserve capacity provided by the transmission network and microgrid i , respectively, while c t R and c i , t R are the corresponding price for providing reserve services. N T is the scheduling cycle, and Δ t is the time step.
The relationship between the provided reserve and inertia and droop coefficients is given by [28]:
R g , t 2 H g , t P g f rate ¯ + D g , t P g Δ f nadir ¯
R i , t 2 H i , t P i f rate ¯ + D i , t P i Δ f nadir ¯ , i
where H g , t / D g , t and H i , t / D i , t are the inertia and droop coefficient from the transmission network and microgrid i at period t , respectively.
(2)
Operational Constraints
In the established distribution network model, PVs and Static Var Compensators (SVCs) are considered. Using second-order cone relaxation of the distflow equations, the steady-state operational constraints are expressed as [29]:
p j , t = k : j k P j k , t u : u j ( P u j , t r u j l u j , t ) , j N bus
q j , t = k : j k Q j k , t u : u j ( Q u j , t x u j l u j , t ) , j N bus
v j , t = v u , t 2 ( r u j P u j , t + x u j Q u j , t ) + ( r u j 2 + x u j 2 ) l u j , t , ( u , j ) N line
p u ¯ p u , t p u ¯ , u N bus
q u ¯ q u , t q u ¯ , u N bus
v u ¯ v u , t v u ¯ , u N bus
l u j l u j ¯ , ( u , j ) N line
2 P u j , t 2 Q u j , t l u j v u , t 2 l u j + v u , t , u N b u s , ( u , j ) N l i n e
p j , t = P j , g , t P j , i , t sell + P j , i , t buy + P j , t PV P j , t load
q j , t = Q j , g , t + Q j , t SVC Q j , t load
where p j , t and q j , t represent the active and reactive power flowing out of node j at period t . P j k , t and Q j k , t represent the active and reactive power flowing from node j to node k at period t , respectively, with j k indicating the direction of power flow. r u j and x u j represent the impedance and reactance parameters between nodes u and j . v j , t denotes the square of the voltage at node j at period t , and l u j , t denotes the square of the line current between nodes u and j . P j , t load and Q j , t load represent the active and reactive load power at node j at period t . P j , g , t and Q j , g , t represent the injected active and reactive power from the transmission network connected to node j at period t ; if node j is not the connection node between the distribution network and the transmission network, they are set to 0. P j , i , t sell and P j , i , t buy represent the selling and purchasing power from the microgrid i connected to node j at period t . If there is no microgrid connection, their values are set to 0. P j , t PV represents the PV generation connected to node j at period t . Q j , t SVC represents the reactive power from the SVC connected to node j at period t . N bus represents the set of nodes in the distribution network, and N line represents the set of lines in the distribution network.
The operational constraints for PVs in the distribution network at period t are:
0 P j , t PV P j , t , max PV , j N PV , t
where P j , t , max PV represents the maximum predicted active power output of the PV.
In the frequency security constraints, it is necessary to consider the uncertainty of disturbances. Disturbances in distribution networks come from PVs power fluctuations ξ PV and load power fluctuations ξ d . Assuming that the disturbance sources ξ PV and ξ d follow a normal distribution [30], the cumulative distribution function (CDF) of the disturbances is obtained based on historical data [10]:
F Δ P ( x ) = Pr ( Δ P x )
where
Δ P = ξ PV + ξ d
Let Pr represent the probability, and set the risk coefficient as α . Then, the disturbance that needs to be satisfied within the confidence interval for each period is:
F Δ P t ( K t , α ) = 1 α
where K t , α is the upper quantile of Δ P t . The dynamic frequency constraints of the distribution network for each period are transformed into:
2 H t , total f rate ¯ K t , α
( D 0 + D t , total ) f ss ¯ K t , α
( β h C + β h H H t , total + β h D D t , total ) × Δ f nadir ¯ K t , α , h
where H t , total and D t , total are the system equivalent inertia and droop coefficient. β h C , β h H and β h D are the hyperplane coefficients mentioned in Section 3.1.1.

3.2. Microgrid Operation Optimization Model

(1)
Objective Function
The optimization objective of microgrid i is to minimize scheduling costs. To account for uncertainty, a stochastic optimization approach was employed, and the impact of uncertainty risk is quantified using conditional value-at-risk (CVaR) [31,32]. The objective function is shown as follows:
C i = min C i E + σ i C i CVaR C i R
where C i E represents the expected dispatch cost for all scenarios, C i CVaR represents the risk cost, and σ i is the risk coefficient used to balance the relationship between dispatch cost and risk cost. C i R represents the reserve benefits. It can be calculated as follows:
C i R = t = 1 N T c i , t R × R i , t Δ t
The dispatch cost includes the degradation cost of electrical storage C i ESS , heat storage systems C i HS , and the cost of electricity interaction C i EX , calculated as follows:
C i E = C i ESS + C i EX + C i HS
Each cost is calculated according to the following:
C i ESS = ω = 1 S ρ i , ω t = 1 N T K ESS P i , t , ω ESS , ch + P i , t , ω ESS , dis Δ t
C i EX = t = 1 N T c t sell P i , t sell Δ t c t buy P i , t buy Δ t
C i HSS = t = 1 N T K HSS P i , t HSS , ch + P i , t HSS , dis Δ t
In (36), K ESS is the degradation coefficient of the energy storage system, and P i , t , ω ESS , ch and P i , t , ω ESS , dis represent the charging and discharging power of the energy storage system, respectively. ρ i , ω is the probability of scenario ω . K H S S represents the degradation coefficient of the heat storage electric boiler, while P i , t HSS , ch and P i , t HSS , dis denote the heat storage power and heat release power of the electric boiler at period t , respectively.
The scheduling risk cost C i CVaR is calculated according to the following:
C i CVaR = ζ + 1 1 α ω = 1 S ρ i , ω C i , ω ζ +
C i , ω = C i , t EX + C i , t HS + t = 1 N T K ESS P i , t , ω ESS , ch + P i , t , ω ESS , dis Δ t
where ζ is an auxiliary variable, whose optimal value is the Value at Risk (VaR), α is the confidence level of CVaR, C i , ω represents the scheduling cost under scenario ω . x + = max x , 0 .
(2)
Operational Constraints
(a)
Electric power Balance Constraint
P i , t PV + Δ P i , t , ω PV P i , t , ω PV , loss + P i , t Wind + Δ P i , t , ω Wind P i , t , ω Wind , loss + P i , t sell P i , t buy = P i , t , ω ESS , ch P i , t , ω ESS , diss + P i , t e , load + Δ P i , t , ω e , load + P i , t EB
where P i , t PV , P i , t Wind , and P i , t e , load represent the predicted power of PV, wind, and electrical loads at period t , respectively; Δ P i , t , ω PV , Δ P i , t , ω Wind , and Δ P i , t , ω e , load represent the power prediction errors of PV, wind, and loads in scenario ω ; P i , t , ω PV , loss and P i , t , ω Wind , loss represent the curtailed power of PV and wind; P i , t , ω ESS , ch and P i , t , ω ESS , diss represent the charging and discharging power of battery energy storage in scenario ω ; P i , t EB represents the electrical power of the heat storage electric boiler.
(b)
Thermal Power Balance Constraint
P ^ i , t he , load + P i , t AC , in + η HSS P i , t HSS , ch P i , t HSS , diss / η HSS = P i , t EB η EB
where P ^ i , t he , load represents the predicted thermal load power, P i , t AC , in represents the input power of the absorption chiller, η HSS represents the efficiency of heat storage and release, and P i , t HSS , ch and P i , t HSS , diss , respectively, represent the heat storage and release power of the heat storage electric boiler.
(c)
Cooling Power Balance Constraint
P ^ i , t c , load = P i , t AC , in η AC
where P ^ i , t c , load represents the predicted cooling load power, and η AC represents the efficiency of the absorption chiller.
(d)
Equipment Constraint
The following formulas represent the upper and lower bounds constraints of the electrical power of the absorption chiller and the heat storage electric boiler.
0 P i , t AC , in P i , max AC , in
0 P t EB P i , max EB
(e)
Battery Energy Storage System Operation Constraints
The State of Charge (SOC) in adjacent periods must satisfy the following relationship:
S O C i , t , ω S O C i , t 1 , ω = P i , t , ω ESS , ch η P i , t , ω ESS , diss / η Δ t E i ESS , cap
In (46), S O C i , t , ω represents the State of Charge (SOC) of the battery energy storage system, and η represents the charging and discharging efficiency. E i ESS , cap represents the maximum capacity of the energy storage system.
The remaining capacity of the battery energy storage system at each period during operation must satisfy the upper and lower bounds constraints as shown in the following:
S O C min S O C i , t , ω S O C max
Additionally, the charging and discharging power of the battery energy storage system must be limited within the following range:
0 P i , t , ω ESS , ch δ i , t , ω ESS λ ESS E i ESS , cap
0 P i , t , ω ESS , dis ( 1 δ i , t , ω ESS ) λ ESS E i ESS , cap
where δ i , t , ω ESS is a binary variable representing the charging and discharging state of the battery energy storage system, where δ i , t , ω ESS = 1 indicates that the system is in charging state, and δ i , t , ω ESS = 0 indicates discharging state; λ ESS is the maximum charging and discharging rate of the battery energy storage system.
Energy storage systems participate in frequency response. To ensure sufficient power and energy during frequency response, reserve constraints should be established [13]:
P i , t , ω ESS , ch + P i , t , ω ESS , dis + R i , t ESS λ ESS E i ESS , cap
S O C min + R i , t ESS Δ t f r S O C i , t , ω S O C max R i , t ESS Δ t f r
where R i , t ESS represents the reserve capacity provided by the energy storage at period t , which is equal to the reserve capacity that the microgrid provides to the distribution network. Δ t f r denotes the primary frequency response time.
(f)
Operation Constraints of the Heat Storage System in the Heat Storage Electric Boiler
The HOC in adjacent time periods must satisfy the following relationship:
H O C i , t H O C i , t 1 = P i , t HSS , ch η HSS P i , t HSS , diss / η HSS Δ t E i HSS , cap
In (52), H O C i , t represents the thermal storage state of the thermal storage system at period t . E i HSS , cap represents the maximum capacity of the thermal storage system.
H O C min H O C i , t H O C max
The remaining capacity of the heat storage system at each time period during operation must satisfy the upper and lower bounds constraints as shown in (53).
0 P i , t HSS , ch δ i , t HSS E i HSS , cap
0 P i , t HSS , diss ( 1 δ i , t HSS ) E i , t HSS , cap
In (54), δ i , t HSS is a binary variable representing the heat storage and release state of the heat storage system, where δ i , t HSS = 1 indicates that the system is in the heat storage state, and δ i , t HSS = 0 indicates the heat release state.
(g)
Electricity Interaction Constraints
The constraints on purchasing and selling electricity are as follows:
0 P i , t s e l l δ i , t EX P max sell
0 P i , t buy 1 δ i , t EX P max buy
where δ i , t EX is a binary variable representing the purchase and sale status of electricity, where δ i , t EX = 1 indicates that the microgrid purchases electricity from the distribution network, and δ i , t EX = 0 indicates that the microgrid sells electricity to the distribution network; P max sell and P max buy represent the maximum purchase and sale power, respectively.

3.3. Coordination Optimization Model of Distribution Networks and Microgrids

The distribution network and microgrids belong to different interest groups and engage in a game over reserve capacity and pricing. The following diagram represents the overall mathematical model:
This paper introduces frequency security constraints (30)–(32) into the traditional distribution network optimization model, which were derived from the frequency response model in Section 3.1.1. For the microgrid optimization model, we considered reserve capacity constraints (50) and (51) to realize inertia and PFR support for distribution networks. The coordinated model as described in Figure 3 was solved to obtain the reserve capacity interaction and reserve price interaction between distribution networks and the microgrids.

4. Cooperative Bargaining Model for Distribution Networks and Microgrids

In the cooperative model, each microgrid and the distribution network are rational entities, prioritizing their own interests when cooperating. Therefore, the cooperative model should fully consider the interests of all parties involved. The distribution network aims to reduce its reserve costs by obtaining inertia and PFR support from microgrids. Meanwhile, microgrids aim to gain reserve benefits, thereby reducing their operational costs. Hence, the cooperative model can be described as follows:
max ( C DN , 0 C DN , co ) i m ( C i 0 C i co ) s . t . C DN , 0 C DN , co 0 , C i 0 C i co 0 , i m
where C DN , co and C i co represent the costs for the distribution network and microgrid i after participating in cooperation, and C DN , 0 and C i 0 represent the costs when not participating in cooperation, which is the negotiation breakdown point. When the microgrid i does not participate in cooperation, the reserve capacity R i is 0.
By solving the equilibrium solution of problem (58), the distribution network and microgrid can derive optimal bargaining strategies R i , c i R .

4.1. Reformatting the Problem

The optimization problem (58) inherently forms a non-convex, nonlinear optimization problem considering reserve capacity and reserve price as optimization variables. To facilitate its solution, this paper transforms it into two sub-problems: minimizing the total system costs and benefit allocation. By sequentially optimizing sub-problems, the optimal solution to the original problem (58) was obtained.
(1)
Sub-problem 1: minimizing the total system costs.
The total system cost is the sum of the objective functions of all entities within the system. Let C g R = t = 1 N T ( R g , t × c g , t ) Δ t , it represents the reserve costs from the transmission network. Therefore, the overall cost is:
min C P , DN + C g R + i m ( C i E + σ i C i CVaR )
In cooperation, the distribution network does not dispatch the purchase and sale power of microgrids; thus, reducing the total system cost is transformed into a problem of reducing reserve costs from the transmission network by utilizing reserves from the microgrid.
(2)
Sub-problem 2: benefit allocation.
By solving sub-problem 1, the reserve capacity after interaction between the distribution network and each microgrid can be obtained. Substituting the reserve into (58), the reserve price can be calculated, and the objective function becomes:
max ( C DN , 0 C P , DN , * C g R , * i m t = 1 N T R i , t * c i , t R ) i m ( C i 0 C i E , * σ i C i CVaR , * + t = 1 N T R i , t * c i , t R )
where the superscript * denotes the solution obtained from solving problem 1. Then, by taking the logarithm, the product problem is transformed into a sum problem, as follows:
min ln ( C DN , 0 C P , DN , * C g R , * i m t = 1 N T R i , t * c i , t R ) i m ln ( C i 0 C i E , * σ i C i CVaR , * + t = 1 N T R i , t * c i , t R )
The derivation process is detailed in reference [33]. Problem 2 is likewise a convex problem, and it was solved using a distributed solving algorithm.

4.2. Distributed Solution Algorithm

Based on the derivation in Section 4.1, the problem was transformed into two sub-problems (59) and (61). The ADMM algorithm was used to solve subproblem 1 and subproblem 2 separately, obtaining the optimal solution to the original problem.

4.2.1. Problem 1: System Cost Minimization

In problem 1, without considering the interaction price of reserve capacity, each microgrid and the distribution network are coupled through the reserve capacity R i , t they provide. Consistency constraints are introduced to achieve decoupling.
R i , t DN = R i , t m
In (62), R i , t DN is the reserve capacity that the distribution network expects the microgrid i to provide, and R i , t m represents the reserve capacity that microgrid i is willing to provide to the distribution network. When (62) holds, problem 1 converges. The ADMM algorithm is used for decomposition and computation, resulting in the optimization model for problem 1.
The optimization model for the distribution network is:
min L pr 1 DN = C P , DN + C g R + i m t = 1 N T λ i , t , pr 1 R i , t DN R i , t m + ρ pr 1 2 R i , t DN R i , t m 2 2 s . t . ( 16 ) ( 26 ) , ( 30 ) ( 32 )
where L pr 1 DN is the augmented Lagrangian function for the distribution network; ρ pr 1 is the penalty factor; λ i , t , pr 1 is the Lagrange multiplier for the interactive power between the distribution network and microgrid i at period t .
The optimization model of microgrid i regarding subproblem 1 is:
min L i , pr 1 = C i E + σ i C i CVaR + t = 1 N T λ i , t , pr 1 R i , t DN R i , t m + ρ pr 1 2 R i , t DN R i , t m 2 2 s . t . ( 41 ) ( 57 )
where L i , pr 1 is the augmented Lagrangian function for microgrid i .
The update equations for variables during the iterative process are as follows:
R i DN ( w + 1 ) = arg min R i m ( w ) L pr 1 DN ( R i m ( w ) , λ i , pr 1 ( w ) ) R i m ( w + 1 ) = arg min R i DN ( w + 1 ) L i ( R i DN ( w + 1 ) , λ i , pr 1 ( w ) ) λ i , pr 1 ( w + 1 ) = λ i , pr 1 ( w ) + ρ pr 1 R i DN ( w + 1 ) R i m ( w + 1 )
where w represents the number of iterative steps. The convergence criterion is given by:
r pr 1 ( w ) = R i DN ( w ) R i m ( w ) 2 ε pr 1 , pri s pr 1 ( w ) = R i DN ( w ) R i DN ( w 1 ) 2 ε pr 1 , dual
where r pr 1 ( w ) and s pr 1 ( w ) represent the primal and dual residuals of subproblem 1 at the w -th iteration, respectively. ε pr 1 , pri and ε pr 1 , dual denote the convergence accuracy of the primal and dual residuals of subproblem 1 at the w -th iteration.

4.2.2. Problem 2: Benefit Allocation

In problem 2, which is the benefit allocation problem, the reserve capacities of each microgrid have already converged. The distribution network and each microgrid use the reserve price c i , t R at each period as the coupling variable. Consistency constraints were introduced to achieve decoupling.
c i , t D = c i , t m
In (67) c i , t D and c i , t m represent the bidding prices of the distribution network and the microgrid, respectively. The optimization model for Problem 2 of the distribution network is as follows:
min L pr 2 DN = ln ( C DN , 0 C P , DN , * C g R , * i m t = 1 N T R i , t * c i , t D ) + i m t = 1 N T λ i , t , pr 2 c i , t DN c i , t m + ρ pr 2 2 c i , t DN c i , t m 2 2 s . t . C DN , 0 C P , DN , * C g R , * i m t = 1 N T R i , t * c i , t D 0
In (68) L pr 2 DN represents the augmented Lagrangian function for the distribution network in problem 2; λ i , t , pr 2 represents the Lagrange multiplier, and ρ pr 2 is the penalty factor.
The optimization model for the microgrid is as follows:
min L i , pr 2 = ln ( C i 0 C i E , * σ i C i CVaR + t = 1 N T R i , t * c i , t m ) + t = 1 N T λ i , t , pr 2 c i , t DN c i , t m + ρ pr 2 2 c i , t DN c i , t m 2 2 s . t . C i 0 C i E , * σ i C i CVaR + t = 1 N T R i , t * c i , t m 0
Variable updates during the iteration process are similar to those in Problem 1 and are not reiterated here.

5. Case Study

5.1. Parameter Settings

This paper assumes that the transmission network is equivalent to a generator. Its parameters are in Table 1 [13]:
The distribution network adopted the IEEE 33-bus case. Three microgrids were distributed at nodes 24, 30, and 12. SVCs were installed at nodes 6, 17, 30, and 33. PV sources with capacities of 300 kW, 600 kW, and 400 kW were installed at nodes 2, 7, and 11, respectively. The structure is shown in Figure 4.
Microgrid 1 had 1200 kW of PV, 500 kW of wind turbines, a 1000 kW heat storage boiler, and a 500 kW absorption chiller. Its cooling, heating, and electrical loads were 120 kW, 120 kW, and 650 kW. Microgrid 2 was identical to Microgrid 1. Microgrid 3 had a 400 kW wind turbine and loads at 90% of Microgrid 1, with the rest of the same setup. Microgrids 1 and 2 had 500 kWh of energy storage, and Microgrid 3 had 400 kWh, all with a 0.5 charge/discharge rate.
The selling price of electricity from the distribution network to the microgrid is shown in Figure 5. The price remained the same throughout each hour. The distribution network purchased electricity from the microgrids and the transmission network at 80% of the selling price, while reserve was priced at 20% of the selling price [14].

5.2. Model and Algorithm Verification

This section first presents the convergence results of the interaction model and then demonstrates the feasibility of the constructed model through a sequential analysis of internal scheduling results of the microgrid, distribution networks, and the interaction results between them.
(1)
Algorithm Convergence Verification
Figure 6 shows the convergence results of the primal and dual residuals of the reserve capacity. It can be observed that the algorithm converged after 63 iterations, which met the requirement for fast convergence.
Figure 7 shows the convergence results of the primal and dual residuals for Problem 2, denoting the benefit allocation process. It can be seen that Problem 2 also converged quickly, meeting the convergence criteria after 18 iterations.
(2)
Verification of Microgrid Optimization Model
To verify the accuracy of the microgrid model, this part presents the internal scheduling results after interaction, using Microgrid 1 as an example.
Figure 8 shows the power generation of the PV and wind turbines in Microgrid 1, displaying the expected values for each scenario on the next day. The power balance results for the electric, thermal, and cooling loads are presented in Figure 9.
Figure 9 shows the power balance of heating, cooling, and electricity in the microgrid. It can be observed that the heating, cooling, and electricity demands were all met. For the electricity balance, we can see that during periods 1–32 and 66–96, the PV and wind power generation were insufficient to meet the microgrid’s electricity demand, resulting in power purchases from the distribution network. During periods 33–66, the PV generation in the microgrid was relatively high, allowing the microgrid to sell excess electricity back to the distribution network. In the heat balance chart, the microgrid showed a significant heat generation from the electric boiler at period 26, driven by a lower electricity price at this time. This allowed the microgrid to store heat initially and release it during periods 72–80, thereby reducing electricity usage in these later periods and lowering power purchase costs.
In the internal scheduling of the microgrid, it responded to time-of-use electricity pricing, which was primarily reflected in the results of the charging and discharging of energy storage systems.
From Figure 10, it can be seen that the energy storage in the microgrid charged during periods 1–32 (off-peak) and discharged during periods 71–80 (peak), thereby increasing the microgrid’s revenue. During periods 33–70, the PV generation was relatively high; due to the maximum selling limit, the charging and discharging behavior of the energy storage was closely related to the PV generation. Charging occurred during peak PV generation periods (e.g., periods 41 and 42), while discharging took place during periods of lower PV output (e.g., periods 37 and 56), aiming to sell more electricity to the distribution network and maximize revenue.
(3)
Verification of Distribution Networks Optimization Model
The frequency response results during the maximum disturbance scenario are depicted in Figure 11 to validate the effectiveness of frequency safety constraints. A disturbance risk coefficient of 0.05 was used, focusing on the period with the highest potential disturbance (0.2715 of the maximum load). Parameters f rate ¯ , Δ f ss ¯ , and Δ f nadir ¯ were set to 0.5 Hz/s, 0.25 Hz, and 0.5 Hz respectively [34].
From Figure 11, it can be observed that the maximum frequency deviation of the distribution network is 0.48 Hz, which was less than the allowable limit of 0.5 Hz for the system. The quasi-steady-state frequency deviation was 0.2 Hz, which was also below the permissible limit of 0.25 Hz. The maximum RoCoF will be provided in the next chapter, and all values remained within the allowable range for the system.
In Figure 12, the proportions of each component participating in frequency response after optimization are given.
From Figure 12, it can be seen that all three microgrids participated in the frequency response, providing part of the inertia support and PFR to the distribution network. Due to the larger energy storage capacity of Microgrid 1 and Microgrid 2, they had a greater ability to participate in the frequency response (active reserves were both 12% for Microgrid 1 and 2, while Microgrid 3 had only 10%).
Figure 13 shows the disturbance magnitude and the required reserve capacity for each period. During periods 1–26 and 77–96 (with 15 min time steps), the output from the PVs in the distribution network was low, resulting in smaller disturbances and lower reserve requirements. In periods 28–76, the distribution network experienced high load levels, and PVs contributed to significant active power fluctuations. The additional reserves provided by the microgrids were insufficient, requiring the transmission network’s participation in providing additional frequency reserves.
The interaction between the distribution network and the microgrid also involves electric energy exchange, during which the voltage should be maintained within a reasonable range. Figure 14 shows the voltage distribution of all nodes in the distribution network throughout the day. It can be observed that during the interaction between the distribution network and the microgrid cluster, the voltage remained within a normal range, ensuring the stable operation of the distribution network.

5.3. Comparison Results

To validate the effectiveness of frequency safety constraints and the economic viability of cooperative game interactions between the distribution network and the microgrid, this study compares the interaction model proposed in reference [18]. In the model presented in reference [18], the distribution network interacts solely with the microgrid in terms of electric energy. In contrast, this study includes not only electric energy interactions with the microgrid but also interactions regarding reserve capacity.
We first compared the frequency response results of the model proposed in this paper with [18].
From Figure 15, it can be seen that with a disturbance of 0.2715 times the maximum load, the frequency of the distribution network could be maintained within the specified range (with a maximum frequency deviation of 0.48 Hz and a quasi-steady-state frequency deviation of 0.2 Hz) after considering frequency safety constraints and reserve capacity. In contrast, when only electric energy interactions were considered, ensuring frequency safety in the distribution network became challenging (with a maximum frequency deviation of 1.77 Hz and a quasi-steady-state frequency deviation of 0.44 Hz).
The frequency response results of the distribution network for each time period are shown in Figure 16, with disturbances selected from the maximum disturbance within the confidence interval.
Figure 16 presents the maximum frequency deviation, maximum rate of change of frequency, and quasi-steady-state frequency deviation encountered during maximum disturbances for each time period. When frequency safety constraints and reserves were not considered, the distribution network experienced significant frequency violations during periods 27–74, primarily due to fluctuations in PV generation. However, after incorporating frequency safety constraints, the distribution network optimized its inertia and droop coefficients, and reserves were provided by the microgrid and transmission network, allowing the system frequency to remain consistently within the safe range.
In the model proposed in this study, the energy storage in the microgrid was required to provide frequency regulation reserves for the distribution network; therefore, its operational state differed from that in the model presented in reference [18]. Figure 17 provides the SOC variation curve of Microgrid 1 for analysis.
Figure 17 shows the SOC curve of Microgrid 1’s battery storage. As shown in Figure 17, in the case of electricity interaction only, the microgrid bought electricity during off-peak hours (1–32, 89–96) and sold electricity during peak hours (69–80) to obtain profits. This resulted in rapid SOC increases during charging in periods 1–32 and significant SOC drops during discharging in periods 69–80. In other periods (33–68), the storage discharged when microgrid generation was insufficient and charged when generation was ample. When considering the interaction between energy and reserve capacity, the storage system similarly charged during off-peak periods, leading to an increase in SOC, thereby providing sufficient upward frequency regulation (energy discharge capability). However, to ensure adequate downward regulation capability (energy absorption), the SOC should not be too high. Therefore, from periods 31 to 69, the SOC remained lower than it would have been in the model proposed in [18]. Conversely, from periods 81 to 86, the SOC was higher to ensure sufficient downward regulation capability (energy absorption). It is also worth noting that the significant discharging of the storage system during peak periods (65–80) was delayed until after period 74. This was because the storage system still needed to provide frequency regulation reserves during periods 65–74, requiring some power reserve and preventing large discharges.
In the mode proposed in this paper, reserve capacity and reserve pricing were determined through a cooperative game between the distribution network and the microgrid. In this cooperative game, both the distribution network and the microgrid can benefit. We also provide a cost comparison to analyze the economic differences between the model in [18] and the model presented in this paper.
Figure 18 displays the negotiated prices between each microgrid and the distribution network. During periods when additional reserves were not required, the negotiated prices were 0. The prices negotiated between the distribution network and microgrids were lower than the reserve prices from the transmission network, enabling the distribution network to achieve cost savings through microgrid reserves.
Table 2 presents the cost breakdown for each entity under two models, where negative values signify profits. In the absence of frequency regulation participation from microgrids, the distribution network’s frequency regulation demand was entirely fulfilled by the transmission network. The result indicates that each microgrid had strategically reduced its electricity trading to enhance its reserve income under the proposed model. Through this collaborative approach, the overall revenue of all three microgrids had increased. Specifically, the reserve costs incurred by the distribution network decreased by RMB 125, and the total cost was reduced by RMB 208, which underscores a successful reduction in operational expenses for distribution networks. Consequently, the collaboration between distribution networks and microgrids not only aids in diminishing the total costs for both the distribution system and the microgrid cluster but also enables all parties to derive benefits from this cooperative game.

6. Conclusions

This paper constructed a Nash bargaining-based coordinated frequency-constrained dispatch model for distribution networks and microgrids, where microgrids gained revenue through providing the inertia and PFR service, and the distribution network reduced its reserve requirement and costs from the transmission grid. The experimental results indicate that the proposed method can ensure frequency safety in the distribution network under disturbances. During the maximum disturbance, the maximum frequency deviation was reduced by 72.6%, the quasi-steady-state frequency deviation was reduced by 71.1%, and the maximum frequency change rate was reduced by 80.2%. Additionally, by utilizing reserve from the microgrid, the adjustment pressure on the transmission network was effectively reduced, with reserve capacity decreasing by 35%. Through Nash bargaining, the reserve costs for the distribution network were reduced by 9.2%, leading to a total cost reduction of 0.73%, and the total revenues for the three microgrids increased by 10.74%, 16.92%, and 15.33%, respectively.
In the future, we will further analyze the interactions and frequency support mechanisms among multiple entities, including the transmission network, distribution network, and microgrids, to fully tap into the flexible resources on the distribution side.

Author Contributions

Conceptualization, Z.Z. and Z.W.; methodology, Z.Z.; investigation, Z.Z.; data curation, Z.Z. and Y.Z.; writing—original draft, Z.Z.; writing—review and editing, Z.Z. and X.W.; funding acquisition, Z.Z. and X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China under Grant 62406099.

Data Availability Statement

The data of this study are available from the corresponding author upon reasonable request. The data of this study are not publicly available due to privacy reasons.

Conflicts of Interest

Author Yanan Zhang was employed by the company State Grid Tianjin Power Dongli Power Supply Branch. 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.

Nomenclature

Indices and Sets
N T Scheduling cycle
N bus Set of nodes in the distribution network
N line Set of lines in the distribution network
R i Reserve strategy set of microgrid i
c i R Reserve price set of microgrid i
m Microgrid cluster
Parameters
P g / P i Maximum output of the transmission network/microgrid i
P d Maximum load demand of distribution networks
D 0 System damping
Δ f Frequency deviation
f rate ¯ Maximum allowable rate of change of frequency
f ss ¯ Maximum allowable quasi-steady-state frequency deviation
Δ f nadir ¯ Maximum allowable frequency deviation
F H High-pressure turbine power generation proportion
T R Generator reheating time
β h C , β h H , β h D A set of hyperplane coefficients
c t buy / c t sell Electricity purchase/sale price of the distribution network during time period t
c t R / c i , t R Reserve service price offered by the transmission network/microgrid i at period t
r u j / x u j Impedance/reactance parameter between nodes u and j
P j , t load / Q j , t load Active and reactive load power at node j at period t
P j , t , max PV Maximum predicted active power output of the PV at node j at period t
ξ PV / ξ d Disturbances in distribution networks come from PVs/load
Δ P t Disturbance in distribution networks at period t
K t , α Upper quantile of disturbance in distribution networks
σ i Risk coefficient of microgrid i
ρ i , ω Probability of scenario ω for microgrid i
K ESS / K H S S Degradation coefficient of the energy storage system/heat storage electric boiler
P i , t PV / P i , t Wind / P i , t e , load Predicted power of PV/wind/electrical loads in microgrid i at period t
P ^ i , t he , load / P ^ i , t c , load Thermal/cooling load power in microgrid i at period t
Δ P i , t , ω PV / Δ P i , t , ω Wind / Δ P i , t , ω e , load Power prediction errors of PV, wind, and loads in scenario ω in microgrid i
η AC / η HSS Efficiency of the absorption chiller/heat storage and release
η EB Efficiency of converting electrical power to thermal energy by the electric boiler
η Efficiency of the charging/discharging efficiency of the battery energy storage system
S O C min / S O C max Upper and lower bounds of SOC
H O C min / H O C max Upper and lower bounds of the thermal storage state of the thermal storage system
P i , max AC , in / P i , max EB Upper bounds of the electrical power of the absorption chiller/the heat storage electric boiler in microgrid i
E i ESS , cap Maximum capacity of the energy storage system in microgrid i
E i HSS , cap Maximum capacity of the thermal storage system in microgrid i
L pr 1 DN / L i , pr 1 / L pr 2 DN / L i , pr 2 Lagrangian function for the distribution network/microgrid i in problem1/2
Variables
H g , t / H i , t Inertia provided by transmission network/microgrid i at period t
D g / D i Droop coefficient provided by transmission network/microgrid i at period t
H t , total Total inertia of the power system at period t
D t , total Total droop coefficient of the power system at period t
C DN Cost of distribution networks
C P , DN Cost of electricity interaction in distribution networks
C R , DN Cost of reserve services in distribution networks
C i Optimization objective of microgrid i
C i E Expected dispatch cost for all scenarios for microgrid i
C i CVaR Risk cost for microgrid i
C i R / C g R Reserve benefits for microgrid i /transmission networks
C i ESS / C i HS Degradation cost of electrical storage/heat storage systems for microgrid i
C i EX Cost of electricity interaction for microgrid i
C DN , co , C i co Costs for the distribution network and microgrid i after participating in cooperation
C DN , 0 , C i 0 Costs for the distribution network and microgrid i before participating in cooperation
P g , t Electricity purchased by the distribution network from the transmission network at period t
P i , t buy / P i . t sell Amount of electricity purchased from/sold to microgrid i at period t
p j , t / q j , t Active/reactive power flowing out of node j at period t
P j k , t / Q j k , t Active/reactive power flowing from node j to node k at period t
v j , t Square of the voltage at node j at period t
l u j , t Square of the line current between nodes u and j
P j , g , t / Q j , g , t Active/reactive power from the transmission network connected to node j at period t
P j , t PV PV generation connected to node j at period t
Q j , t SVC Reactive power from the SVC connected to node j at period t
ζ Auxiliary variable
α Confidence level of CvaR
P i , t , ω PV , loss / P i , t , ω Wind , loss Curtailed power of PV and wind in scenario ω in microgrid i
P i , t , ω ESS , ch / P i , t , ω ESS , diss charging and discharging power of battery energy storage in scenario ω in microgrid i
P i , t HSS , ch / P i , t HSS , diss Heat storage/release power of the heat storage electric boiler
P i , t EB Electrical power of the heat storage electric boiler of microgrid i at period t
P i , t AC , in Input power of the absorption chiller of microgrid i at period t
S O C i , t , ω SOC of the battery energy storage system in scenario ω in microgrid i at period t
H O C i , t Thermal storage state of the thermal storage system in microgrid i at period t
R g , t / R i , t Reserve capacity provided by the transmission network/microgrid i at period t
R i , t ESS Reserve capacity provided by the energy storage in microgrid i at period t
δ i , t , ω ESS Charging and discharging state of the battery energy storage system in scenario ω in microgrid i at period t
δ i , t HSS , δ i , t EX Charging/discharging state of thermal energy and electricity trading state in microgrid i at period t

References

  1. Li, L.L.; Pei, J.R.; Shen, Q. A Review of Research on Dynamic and Static Economic Dispatching of Hybrid Wind-Thermal Power Microgrids. Energies 2023, 16, 3985. [Google Scholar] [CrossRef]
  2. Peyerl, D.; Barbosa, M.O.; Ciotta, M.; Pelissari, M.R.; Moretto, E.M. Linkages between the Promotion of Renewable Energy Policies and Low-Carbon Transition Trends in South America’s Electricity Sector. Energies 2022, 15, 4293. [Google Scholar] [CrossRef]
  3. Ding, T.; Zeng, Z.; Qu, M.; Catalao, J.P.S.; Shahidehpour, M. Two-Stage Chance-Constrained Stochastic Thermal Unit Commitment for Optimal Provision of Virtual Inertia in Wind-Storage Systems. IEEE Trans. Power Syst. 2021, 36, 3520–3530. [Google Scholar] [CrossRef]
  4. Trovato, V. System Scheduling With Optimal Time-Varying Delivery Intervals for Frequency Response. IEEE Trans. Power Syst. 2022, 37, 4270–4285. [Google Scholar] [CrossRef]
  5. Ratnam, K.S.; Palanisamy, K.; Yang, G. Future low-inertia power systems: Requirements, issues, and solutions—A review. Renew. Sustain. Energy Rev. 2020, 124, 109773. [Google Scholar] [CrossRef]
  6. Yang, M.; Liu, Y.; Guo, L.; Wang, Z.; Zhu, J.; Zhang, Y.; Cheng, B.; Wang, C. Hierarchical Distributed Chance-constrained Voltage Control for HV and MV DNs Based on Nonlinearity-Adaptive Data-driven Method. IEEE Trans. Power Syst. 2024, 1–14. [Google Scholar] [CrossRef]
  7. Yang, Y.; Chih-Hsien, J.; Ye, Z.-S. Distributionally robust frequency dynamic constrained unit commitment considering uncertain demand-side resources. Appl. Energy 2023, 331, 120392. [Google Scholar] [CrossRef]
  8. Ahmadi, H.; Ghasemi, H. Security-Constrained Unit Commitment With Linearized System Frequency Limit Constraints. IEEE Trans. Power Syst. 2014, 29, 1536–1545. [Google Scholar] [CrossRef]
  9. Li, Z.; Guo, Q.; Sun, H.; Wang, J. Coordinated Economic Dispatch of Coupled Transmission and Distribution Systems Using Heterogeneous Decomposition. IEEE Trans. Power Syst. 2016, 31, 4817–4830. [Google Scholar] [CrossRef]
  10. Doenges, K.; Egido, I.; Sigrist, L.; Miguélez, E.L.; Rouco, L. Improving AGC Performance in Power Systems With Regulation Response Accuracy Margins Using Battery Energy Storage System (BESS). IEEE Trans. Power Syst. 2020, 35, 2816–2825. [Google Scholar] [CrossRef]
  11. Zhang, Z.; Du, E.; Teng, F.; Zhang, N.; Kang, C. Modeling Frequency Dynamics in Unit Commitment With a High Share of Renewable Energy. IEEE Trans. Power Syst. 2020, 35, 4383–4395. [Google Scholar] [CrossRef]
  12. Liang, Z.R.; Mieth, R.; Dvorkin, Y. Inertia Pricing in Stochastic Electricity Markets. IEEE Trans. Power Syst. 2023, 38, 2071–2084. [Google Scholar] [CrossRef]
  13. Shen, Y.; Wu, W.; Wang, B.; Sun, S. Optimal Allocation of Virtual Inertia and Droop Control for Renewable Energy in Stochastic Look-Ahead Power Dispatch. IEEE Trans. Sustain. Energy 2023, 14, 1881–1894. [Google Scholar] [CrossRef]
  14. Tian, Y.; Li, Z.; Wu, W.; Fan, M. Joint chance-constrained economic dispatch involving joint optimization of frequency-related inverter control and regulation reserve allocation. CSEE J. Power Energy Syst. 2024, 1–13. [Google Scholar] [CrossRef]
  15. Nakiganda, A.M.; Dehghan, S.; Markovic, U.; Hug, G.; Aristidou, P. A Stochastic-Robust Approach for Resilient Microgrid Investment Planning Under Static and Transient Islanding Security Constraints. IEEE Trans. Smart Grid 2022, 13, 1774–1788. [Google Scholar] [CrossRef]
  16. Cai, S.; Xie, Y.; Zhang, Y.; Zhang, M.; Wu, Q.; Guo, J. A Simulation-Assisted Proactive Scheduling Method for Secure Microgrid Formation Under Static and Transient Islanding Constraints. IEEE Trans. Smart Grid 2024, 15, 272–285. [Google Scholar] [CrossRef]
  17. Tian, Y.; Li, Z. Coordinated frequency-constrained stochastic economic dispatch for integrated transmission and distribution system via distributed optimization. CSEE J. Power Energy Syst. 2024, 1–13. [Google Scholar] [CrossRef]
  18. Liu, H.; Li, J.; Ge, S.; He, X.; Li, F.; Gu, C. Distributed Day-Ahead Peer-to-Peer Trading for Multi-Microgrid Systems in Active Distribution Networks. IEEE Access 2020, 8, 66961–66976. [Google Scholar] [CrossRef]
  19. Zhao, H.R.; Wang, X.J.; Wang, Y.W.; Li, B.K.; Lu, H. A dynamic decision-making method for energy transaction price of CCHP microgrids considering multiple uncertainties. Int. J. Electr. Power Energy Syst. 2021, 127, 106592. [Google Scholar] [CrossRef]
  20. Xu, H.X.; Sun, J.H.; Huang, J.G.; Lin, X.Y.; Ma, C.H. Distributed Optimization of Islanded Microgrids Integrating Multi-Type VSG Frequency Regulation and Integrated Economic Dispatch. Energies 2024, 17, 1618. [Google Scholar] [CrossRef]
  21. Hamanah, W.M.; Shafiullah, M.; Alhems, L.M.; Alam, M.S.; Abido, M.A. Realization of Robust Frequency Stability in Low-Inertia Islanded Microgrids With Optimized Virtual Inertia Control. IEEE Access 2024, 12, 58208–58221. [Google Scholar] [CrossRef]
  22. Wang, S.; Pu, L.Y.; Huang, X.D.; Yu, Y.F.; Shi, Y.W.; Wang, H.W. Online ADMM for Distributed Optimal Power Flow via Lagrangian Duality. Energies 2022, 15, 9525. [Google Scholar] [CrossRef]
  23. Xu, J.Z.; Yi, Y.Q. Multi-microgrid low-carbon economy operation strategy considering both source and load uncertainty: A Nash bargaining approach. Energy 2023, 263, 125712. [Google Scholar] [CrossRef]
  24. Shi, Q.; Li, F.; Cui, H. Analytical Method to Aggregate Multi-Machine SFR Model With Applications in Power System Dynamic Studies. IEEE Trans. Power Syst. 2018, 33, 6355–6367. [Google Scholar] [CrossRef]
  25. Guggilam, S.S.; Zhao, C.; Dall’Anese, E.; Chen, Y.C.; Dhople, S.V. Optimizing DER Participation in Inertial and Primary-Frequency Response. IEEE Trans. Power Syst. 2018, 33, 5194–5205. [Google Scholar] [CrossRef]
  26. Yang, J.W.; Zhang, N.; Kang, C.Q.; Xia, Q. A State-Independent Linear Power Flow Model With Accurate Estimation of Voltage Magnitude. IEEE Trans. Power Syst. 2017, 32, 3607–3617. [Google Scholar] [CrossRef]
  27. Paturet, M.; Markovic, U.; Delikaraoglou, S.; Vrettos, E.; Aristidou, P.; Hug, G. Stochastic Unit Commitment in Low-Inertia Grids. IEEE Trans. Power Syst. 2020, 35, 3448–3458. [Google Scholar] [CrossRef]
  28. Yang, L.; Li, Z.; Xu, Y.; Zhou, J.; Sun, H. Frequency Constrained Scheduling Under Multiple Uncertainties via Data-Driven Distributionally Robust Chance-Constrained Approach. IEEE Trans. Sustain. Energy 2023, 14, 763–776. [Google Scholar] [CrossRef]
  29. Abdelouadoud, S.Y.; Girard, R.; Neirac, F.P.; Guiot, T. Optimal power flow of a distribution system based on increasingly tight cutting planes added to a second order cone relaxation. Int. J. Electr. Power Energy Syst. 2015, 69, 9–17. [Google Scholar] [CrossRef]
  30. Shukla, A.; Singh, S.N. Clustering based unit commitment with wind power uncertainty. Energy Convers. Manag. 2016, 111, 89–102. [Google Scholar] [CrossRef]
  31. Ran, X.; Zhang, J.; Liu, K. An Interval–Probabilistic CVaR (IP-CVaR) and Modelling for Unknown Probability Distribution of Some Random Variables. IEEE Trans. Power Syst. 2023, 38, 2035–2045. [Google Scholar] [CrossRef]
  32. Wang, J.; Zhou, Y.; Zhang, Y.; Lin, F.; Wang, J. Risk-Averse Optimal Combining Forecasts for Renewable Energy Trading Under CVaR Assessment of Forecast Errors. IEEE Trans. Power Syst. 2024, 39, 2296–2309. [Google Scholar] [CrossRef]
  33. Qin, M.L.; Xu, Q.S.; Liu, W.; Xu, Z.M. Low-carbon economic optimal operation strategy of rural multi-microgrids based on asymmetric Nash bargaining. IET Gener. Transm. Distrib. 2024, 18, 24–38. [Google Scholar] [CrossRef]
  34. Yuan, Y.; Zhang, Y.; Wang, J.; Liu, Z.; Chen, Z. Enhanced Frequency-Constrained Unit Commitment Considering Variable-Droop Frequency Control From Converter-Based Generator. IEEE Trans. Power Syst. 2023, 38, 1094–1110. [Google Scholar] [CrossRef]
Figure 1. Interaction structure diagram.
Figure 1. Interaction structure diagram.
Energies 17 05661 g001
Figure 2. Equivalent system frequency response model control block diagram.
Figure 2. Equivalent system frequency response model control block diagram.
Energies 17 05661 g002
Figure 3. Mathematical model for coordination between distribution networks and microgrids.
Figure 3. Mathematical model for coordination between distribution networks and microgrids.
Energies 17 05661 g003
Figure 4. IEEE 33-bus test case structure.
Figure 4. IEEE 33-bus test case structure.
Energies 17 05661 g004
Figure 5. The selling price of electricity from the distribution network.
Figure 5. The selling price of electricity from the distribution network.
Energies 17 05661 g005
Figure 6. Primal residual and dual residual of Problem 1.
Figure 6. Primal residual and dual residual of Problem 1.
Energies 17 05661 g006
Figure 7. Primal residual and dual residual of Problem 2.
Figure 7. Primal residual and dual residual of Problem 2.
Energies 17 05661 g007
Figure 8. PV and wind power generation.
Figure 8. PV and wind power generation.
Energies 17 05661 g008
Figure 9. Power balance results of the electric, thermal, and cooling loads in Microgrid 1.
Figure 9. Power balance results of the electric, thermal, and cooling loads in Microgrid 1.
Energies 17 05661 g009aEnergies 17 05661 g009b
Figure 10. Charging and discharging results of the energy storage system in Microgrid 1.
Figure 10. Charging and discharging results of the energy storage system in Microgrid 1.
Energies 17 05661 g010
Figure 11. Frequency response results of the distribution network.
Figure 11. Frequency response results of the distribution network.
Energies 17 05661 g011
Figure 12. Frequency response coefficient ratio.
Figure 12. Frequency response coefficient ratio.
Energies 17 05661 g012
Figure 13. Reserve for each period.
Figure 13. Reserve for each period.
Energies 17 05661 g013
Figure 14. Node Voltage of the Distribution Network at Each Time Period.
Figure 14. Node Voltage of the Distribution Network at Each Time Period.
Energies 17 05661 g014
Figure 15. Frequency response process.
Figure 15. Frequency response process.
Energies 17 05661 g015
Figure 16. Frequency response results for each period.
Figure 16. Frequency response results for each period.
Energies 17 05661 g016aEnergies 17 05661 g016b
Figure 17. SOC of energy storage.
Figure 17. SOC of energy storage.
Energies 17 05661 g017
Figure 18. Negotiated reserve capacity prices.
Figure 18. Negotiated reserve capacity prices.
Energies 17 05661 g018
Table 1. Equivalent generator parameters for the transmission network.
Table 1. Equivalent generator parameters for the transmission network.
ParameterValue Range
Inertia H g ( s ) 2–20
Droop Coefficient D g ( p . u . ) 15–60
Reheat Time Constant T R ( s ) 8
Aggregated fraction F H 0.25
Table 2. Comparison of cost for each component.
Table 2. Comparison of cost for each component.
EntityModel Without Considering Reserve (RMB)The Proposed Model (RMB)
Electricity TradingReserve CostTotal CostElectricity TradingReserve CostTotal Cost
Distribution network26,839135828,19726,756123327,989
Microgrid 1−12630−1396−1235−181−1546
Microgrid 2−8160−857−791−171−1002
Microgrid 3−7260−763−704−138−880
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

Zhou, Z.; Wang, Z.; Zhang, Y.; Wang, X. Nash Bargaining-Based Coordinated Frequency-Constrained Dispatch for Distribution Networks and Microgrids. Energies 2024, 17, 5661. https://doi.org/10.3390/en17225661

AMA Style

Zhou Z, Wang Z, Zhang Y, Wang X. Nash Bargaining-Based Coordinated Frequency-Constrained Dispatch for Distribution Networks and Microgrids. Energies. 2024; 17(22):5661. https://doi.org/10.3390/en17225661

Chicago/Turabian Style

Zhou, Ziming, Zihao Wang, Yanan Zhang, and Xiaoxue Wang. 2024. "Nash Bargaining-Based Coordinated Frequency-Constrained Dispatch for Distribution Networks and Microgrids" Energies 17, no. 22: 5661. https://doi.org/10.3390/en17225661

APA Style

Zhou, Z., Wang, Z., Zhang, Y., & Wang, X. (2024). Nash Bargaining-Based Coordinated Frequency-Constrained Dispatch for Distribution Networks and Microgrids. Energies, 17(22), 5661. https://doi.org/10.3390/en17225661

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