Next Article in Journal
Between Fast and Sustainable Fashion: The Attitude of Young Lithuanian Designers to the Circular Economy
Next Article in Special Issue
Energy Consumption and Carbon Emissions of Nearly Zero-Energy Buildings in Hot Summer and Cold Winter Zones of China
Previous Article in Journal
Knowledge Sharing among Students in Social Media: The Mediating Role of Family and Technology Supports in the Academic Development Nexus in an Emerging Country
Previous Article in Special Issue
The Environmental Sustainability Study of an Airport Building System Based on an Integrated LCA-Embodied Energy (Emergy)-ANN Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flexibility Analysis for Multi-Energy Microgrid and Distribution System Operator under a Distributed Local Energy Market Framework

1
Center of Excellence for Power Systems Automation and Operation, School of Electrical Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran
2
College of Engineering, Lawrence Technological University, Southfield, MI 48075, USA
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(13), 9985; https://doi.org/10.3390/su15139985
Submission received: 5 May 2023 / Revised: 6 June 2023 / Accepted: 14 June 2023 / Published: 23 June 2023
(This article belongs to the Special Issue Low-Carbon Buildings and Climate Change Mitigation)

Abstract

:
Increasing the penetration rate of microgrids (MGs) for Local Energy Market (LEM) participation creates new challenges for the market-clearing process under a large number of requests for energy transactions. The market-clearing process for decentralized market frameworks is dependent on participants’ flexibility in negotiations for bilateral energy transactions. Multi-energy microgrids (MEMGs) include combined heat and power units which can be less dependent on electricity prices because of energy conversion equipment, gas infrastructure, and combined heat and power loads. In this regard, to evaluate prosumers’ flexibility role in market negotiations, a new analysis based on energy scheduling of MEMG considering a Demand-Response Program (DRP) model is executed under a distributed market structure. Moreover, two new flexibility indexes for market participants with attention to prosumers’ adaption capabilities are proposed. The results show that, under a 9.35% flexibility index improvement for the entire system, the social welfare function improved by 2.75%. Moreover, the results show that the DRP model for changeable and shiftable loads can improve the flexibility of the entire system by 35.82%. Combined heat and power load are considered as the resource of flexibility for system evaluations.

1. Introduction

1.1. Motivation and Aims

Traditionally, power systems were composed of large centralized generators usually operated in a radial form to feed large and small customers. However, these days generation can be owned by companies or individual customers such as MGs, which are not operated by Distribution System Operators (DSOs) [1]. These independent prosumers can participate in local energy markets or ancillary service markets to gain financial profits and support DSO from price fluctuations of the upstream market [2]. These new energy-trading opportunities would create new challenges for local markets because a large number of independent prosumers would increase the computation burden of the market-clearing process [3]. The flexibility of prosumers would have an intensive effect on the market-clearing process, as it would be simpler to match more flexible players to achieve an agreement on energy-trading values. In this context, by establishing reliable two-way communication systems between distributed independent prosumers, a new aspect of flexibility would be considered based on the adaptability features of each prosumer. How can this kind of flexibility, however, be measured? Generally, flexibility is calculated for energy-generation units with evaluation metrics associated with their technical characteristics. Several researchers have studied the concept, effects, and measurement metrics to evaluate a single power plant’s flexibility [4,5]. However, loss of load probability and expected unsupplied energy with an attitude of preventing the demand exceeds from the available capacity, can not be more practical to evaluate a system with the need to predict the behaviors [6]. New market structures create new requirements for the metrics of the system that can evaluate the market participants’ behavior from a market adaptation point of view. So new metrics need to be provided in response to the flexibility requirements of the market-clearing process. The flexibility index of prosumers can be used by market operators as a new instrument to improve the social welfare of the system.

1.2. Literature Review

In a power system context, the term flexibility refers to diverse things ranging from the quick response times of generation units to the degree of efficiency or robustness of a given power market setup [7]. With more attention paid to market concepts, flexibility is defined by Euro Electric as the ability to adapt the generation or consumption pattern in response to an external signal, which would be a price signal or any other incentive or inhibitory parameter [8]. This concept distinguishes a feature of flexibility in the integration of energy scheduling and market structure for signal generation. MEMGs as a new emerging part of power systems have improved the flexibility of MG operation by using Combined Heat and Power units (CHP), energy converter equipment such as Heat Pumps (HPs), and electrical and thermal energy-storage devices [9]. Moreover, demand-response programs can provide a significant opportunity for MEMGs to enhance the flexibility aspects of the systems [10]. Ref. [11] denotes the demand-flexibility definition as shifting eligible loads across the hours of a day to off-peak periods and reshaping the daily load curve to reduce generation costs. Although the traditional DRP models can control a specific kind of demand, a more integrated formulation would be required for combined heat and power loads in multi-energy systems. Ref. [12] defines a combined model of loads for changeable demands, which can be consumed in a thermal or electrical model. However, the proposed model is restricted to the changeability feature of the combined loads, and the shift ability of the DRP is not considered in the article. MEMGs with various energy resources can provide new features for a system by postponing the use of energy under energy-storage devices and changing the kind of demand with energy-converting equipment. These capabilities can create a new model of demand responses for combined heat and power loads, which have not been studied in recent papers. Many articles have been written on the optimal energy scheduling of MEMGs, each of which considered some aspects of the issue, but none of them have studied the flexibility analysis of MEMGs considering the flexibility features for MEMGs for market participation. Ref. [13] propose an optimal energy management model for multi-microgrids considering DRP with only an electric energy carrier. To consider the effects of multi-energy systems, Ref. [14] presents the optimal short-term scheduling for an MEMG considering load uncertainty. However, the role of a Distribution System (DS) is not considered in this study. DSO can play the role of an aggregator for prosumers such as MGs or energy hubs to trade energy with TSO as an intermediate platform [15,16,17]. The role of DSO as an aggregator would prevent a high-dimensional optimization problem with a large number of small controllable elements, which need to be modeled individually [18]. However, the security constraints of the DS would influence the energy-trading schedule according to the bus voltage and line current restrictions [19]. So the proposed model would consider both sides of the model from MG and DSO points of view. Ref. [20] proposes a bi-level problem structure for the operation scheduling of MG with the interaction of the distribution system. Although this study considered the role of DSO in the interaction of MGs, the problem is formulated neglecting a flexibility analysis and energy-trading negotiations. To investigate the effects of heat power, Ref. [21] uses the time-of-use tariff released by the main grid as fixed prices for energy trading of multi-energy peers. Ref. [22] proposes an energy scheduling model for multi-energy systems under two separate coalitions for heat and power energy transactions. Ref. [23] proposes uniform pricing and a pay-as-bid strategy for energy trading. However, static pricing strategies are not able to dynamically consider the effects of market-clearing prices on energy transactions. So, these methods are not good solutions for flexibility analysis. To consider the capability of negotiations for energy transactions, Ref. [24] used the Alternating Direction Method of Multipliers (ADMM) method for bilateral energy trading. Although the proposed model considers the privacy of data for market participants, the study ignored the security constraints of the distribution system, which would decrease the practicality of the model in real operation. Ref. [25] used the ADMM method for heat and power energy trading. Although the proposed model considers the energy-trading negotiation possibility under the security constraints of the system, flexibility analysis is ignored and no evaluation metric is presented in the study. Ref. [26] proposes a Flexibility Index (FI) based on the comparison of the largest variation range of uncertainty and the target range within which the system can remain feasible under a given response time horizon [27]. Ref. [28] presents a detailed flexibility study to incorporate net load forecast errors. The proposed FIs are employed to evaluate the reserve capacity requirements to cover the uncertainty of generation. Ref. [29] provides three flexibility quantification metrics: expected unserved flexible energy, expected duration of insufficient flexibility, and expected flexibility index to design an appropriate DRP. Similarly, Ref. [30] introduces the average power flexibility during the peak period index to evaluate the energy flexibility of the DRPs. To provide demand-response services, Ref. [31] presents an FI, based on a probabilistic curve of the available power from the thermostatically controlled loads. Although the proposed indexes were suitable indicators to evaluate the reliability of the system for flexibility services, none of them can be used to predict the flexibility capability of a system for market participation. In [32], with a focus on the flexible energy of building energy systems, new FIs are proposed to show a system’s capabilities in delivering power and coordination capabilities of resources. Although the different types of demands are considered in this study, the conversion of loads from one carrier to the others is not considered. Ref. [33] proposes a stochastic economical flexibility-evaluation method to scrutinize the flexibility sensitivity to the storage facility capacity. The proposed model indicates that the proposed FIs can help network operators to be aware of the available capabilities of the system. Several metrics have been proposed to measure the flexibility of a system, but no metric has been presented up to now to measure the adaptability of a system for the coordination of an independent operator with the others. This is a feature that can be used in the market-clearing process to make a margin for a market operator for easier market settlement. Moreover, several FIs have been presented to quantify the individual flexibility of an overall system, but none of them are adapted to the market requirements. This is because it is dependent on several parameters and cannot be calculated under the summation of different indicators. For example, the system flexibility cannot be calculated by using a weighted sum of the values of the flexibility indicators for all generators, because there is a relation between the operating constraints of the system and cost with flexibility that restricts the overall flexibility [34]. Table 1 reports a summary of FIs. It can be cocluded that there is no FI that could reflect the adaptability of a system according to a received signal such as prices for market participation. This is a new metric for evaluating the adoption of prosumers that can be used for the market-clearing process or ranking of prosumers in future market-clearing designs.

1.3. Contribution and Innovations

In this paper, a new model of energy scheduling for MEMGs is presented under the flexibility analysis for bilateral market energy transactions. MEMGs with electricity and gas infrastructure can provide a specific percentage of their demand from both gas and electricity energy carriers. In this regard, a new model of DRP is proposed in this paper for combined heat and power demands, which can be adapted to the market prices according to each energy carrier’s cost. According to recent studies, it can be deduced that, despite paying a lot of attention to the concept of flexibility and energy scheduling separately, no quantitative index has been presented up to now to evaluate the flexibility of prosumers according to the market negotiations. In this regard, this paper proposes a new index of flexibility for prosumers based on capabilities to change the bids/offers submitted to the local energy markets. A local energy market based on a distributed market framework is executed, which runs the market for a distribution system operator and connected multi-energy micro grid to study the effects of the flexibility of prosumers in energy-trading negotiations. By using a two-directional structure for market negotiations, the proposed indexes are evaluated from the market-clearing process point of view. As the effects of DRP on the improvement of flexibility are investigated in several types of research [36], MEMGs are employed for flexibility analysis with a new model of DRP for combined loads. In this way, the most prominent features of the current work are as follows:
  • Providing a new model of energy scheduling for multi-energy microgrids under the employment of a demand-response program for changeable and shiftable loads;
  • Providing new indexes to quantitatively evaluate the flexibility of prosumers under market negotiations;
  • Analysis of the operation schedule, market-clearing process, and energy-trading negotiations from the flexibility point of view;
  • Providing an integrated demand-response program model for combined heat and power loads.

1.4. Paper Arrangement

This article is arranged as follows:
A comprehensive mathematical formulation for the operation scheduling of MEMGs and DSO is provided in Section 2. Section 2.3 describes the centralized and proposed distributed solution approaches for bilateral energy. The results of the two structures are compared in the analysis of Section 3. Section 3.2 provides a flexibility assessment of the proposed demand-response program in MEMGs. It also defines two flexibility indexes to evaluate the prosumer adjustment capability in energy-trading negotiations. Finally, Section 4 provides a comprehensive summary of the work carried out.

2. Operation Scheduling

2.1. Multi-Energy Microgrid

MEMGs have electrical and heat demands that can be supplied using MG sources or energy trading with the distribution system. MGs in this study include CHP, Microturbines (MTs), Boilers (Bo’s), HP, Electrical Energy Systems (EESs), and Thermal Energy Systems (TESs). In this section, the mathematical model of the operation scheduling for MEMGs is presented. Figure 1 shows the different components of each MEMG.

2.1.1. Objective Function

The energy generation cost in an MG is the central part of the objective function in the problem. The results of energy management in an MG include generation scheduling for each resource, making decisions about the amount of power trading with DS, and implementation of DRP for load management. In this paper, it is assumed that a specific percent of each MG load includes combined electrical and heat demands that are changeable, and can be shifted to another time following the marginal cost of electricity and heat. Equation (1) defines the objective function of each MEMG.
Min C M G i = t T ( C i , t C H P + C i , t M T + C i , t B o R i , t )

2.1.2. Constraints

a.
CHP
The electrical and thermal output power of CHP are dependent on each other and form a closed curve called FOR. In this model, the convex form of the CHP model is used [37]. Equation (2) indicates the constraints for electrical and thermal power generation of CHP units concerning the Feasible Operating Region (FOR) zones. Equations (3) and (4) consider the minimum cost of on/off for CHP units. Finally, in accordance with the start-up/shut-down costs and the amount of electrical and thermal power generated, the total cost function of CHP is calculated in Equation (5).
{ P i , t C H P , H i , t C H P } F O R i C H P
C i , t C H P , S U ρ i C H P , S U x i , t C H P x i , t 1 C H P ; i , t
C i , t C H P , S D ρ i C H P , S D x i , t C H P x i , t 1 C H P ; i , t
C i , t T , C H P = P i , t C H P + H i , t C H P λ t N G η i C H P + C i , t C H P , S U + C i , t C H P , S D ; i , t
b.
Microturbine
The MT unit has been used to generate electric power. Equation (6) refers to the maximum power output of an MT unit. Equation (7) shows the cost function for power generation.
P i , t M T P i M T , m a x ; i , t
C i , t T , M T = P i , t M T λ t N G η i M T ; i , t
c.
Heat Pump
An MT HP works like a demand converter that can protect the MEMG from natural gas price fluctuations. Equation (8) determines the maximum value of heat generation, and Equation (9) indicates the maximum heat power generation for an HP.
H i , t H P H i , t H P , m a x ; i , t
H i , t H P = α P i , t H P ; i , t
d.
Boiler
Equation (10) presents the heat-generation range of the boiler. Equation (11) shows the boiler cost function. Boiler operation constraints and cost functions are in accordance with [38].
H i , t B o H i , t B o , m a x ; i , t
C i , t T , B o = H i , t B o λ t N G η i B o ; i , t
e.
Electrical Energy Storage
Due to the impossibility of simultaneous charging and discharging of EES units, relations (12) and (13) are formulated by considering the binary variable x i , t E E S . Equation (14) shows the amount of energy stored in the storage in accordance with the charges and discharges of each hour. Equation (15) represents the energy balance equation for an EES unit. Equation (16) indicates the state of the charge range of the storage unit.
H i , t C h , T E S H i C h , T E S , m a x x i , t T E S ; i , t
H i , t d C h , T E S H i d C h , T E S , m a x 1 x i , t T E S ; i , t
S O C i , t T E S = S O C i , t 1 T E S + H i , t C h , T E S η i C h , T E S H i , t d C h , T E S η i d C h , T E S ; i , t
S O C i , 1 T E S = S O C i , 24 T E S ; i
S O C i , t T E S , m i n S O C i , t T E S S O C i , t T E S , m a x ; i , t
f.
Thermal Energy Storage
Constraints (17)–(21), similar to the EES, represent the operational constraints for TES.
H i , t C h , T E S H i C h , T E S , m a x x i , t T E S ; i , t
H i , t d C h , T E S H i d C h , T E S , m a x 1 x i , t T E S ; i , t
S O C i , t T E S = S O C i , t 1 T E S + H i , t C h , T E S η i C h , T E S H i , t d C h , T E S η i d C h , T E S ; i , t
S O C i , 1 T E S = S O C i , 24 T E S ; i
S O C i , t T E S , m i n S O C i , t T E S S O C i , t T E S , m a x ; i , t
g.
Demand-Response Program
In this paper, by considering a specific percent of the MEMG loads which can be changed between the electrical and heat model, a new model of DRP for combined heat and power loads is proposed. It is based on the proposed model in accordance with the changeability and shiftability feature of the combined loads, which can be scheduled based on the price of electricity and gas energy carrier to reduce the operating cost. The DRP model for changeable loads considers consumption scheduling so that by maintaining consumption values for combined loads over 24 h, electricity and gas consumption on-peak price time durations are shifted to lower-priced hours of an energy carrier under the new time-of-use DRP model. Equation (22) shows the model of changeable load over 24 h. Equation (23) indicates the maximum changeable electrical load, and Equation (24) refers to the maximum changeable heat load. Although a small percentage of an MG load can be considered as combined heat and electrical demand, these values can be significant at peak times.
t = 1 24 ( e c h i , t + h c h i , t ) = 0 ; i
C h i , t M G , m a x e c h i , t C h i , t M G , m a x ; i , t
C h i , t M G , m a x h c h i , t C h i , t M G , m a x ; i , t
h.
Energy Trading with DSO
MEMGs can trade power in LEM in a connected mode of operation. Equation (25) represents the revenue of power trading with DSO. As the output power of the network on the two sides of the transformer connecting the MEMGs to the DS is different, Equation (26) makes a difference in generating power and final power value for selling. Equation (26) assumes the power output of the transformer as the sales power of the MEMG and also prevents simultaneous buy and sell bid generation for an MG.
R i , t = P i , t S e l l P i , t B u y λ i , t ; i , t
P i , t S e l l = P i , t S e l l . η i t r a n s ; i , t
i.
Power Equality
Equation (27) shows the electrical power equality, and Equation (34), similarly, states the equality constraint for heat generation and consumption in an MEMG.
P i , t C H P + P i , t M T + P i , t d C h , E E S + P i , t B u y = P i , t M G + P i , t H P + P i , t C h , E E S + P i , t S e l l + e c h i , t ; i , t
H i , t H P + H i , t C H P + H i , t B o + H i , t d C h , T E S = H i , t M G + H i , t C h , T E S + h c h i , t ; i , t

2.2. Distribution System

2.2.1. Objective Function

The DSO optimizes its objective function intending to provide its demand with the lowest cost and maximize the profit of energy trading with MGs. Equation (29) presents the revenue function for DSO as costs of supplying demand through DGs, the revenue of energy trading with MGs, and the revenue achieved via power trading with the upstream grid as TSO.
M a x R D S O = i = 1 N t = 1 T P ^ b , t b = i S e l l P ^ b , t b = i B u y × i , t + t = 1 T P b 1 , t S e l l , T S O λ t S e l l , T S O P b 1 , t B u y , T S O λ t B u y , T S O b = 1 N b t = 1 T λ b D G P b , t D G

2.2.2. Problem Constraints

The following equations are formulated according to the load flow model presented in reference [39], under the radial and active structure of the distribution network. Figure 2 shows a typical radial structure of the distribution system. Equation (30) expresses the voltage limits of the buses. Equation (31) considers the effects of transformer efficiency on energy generation for selling. Equations (32) and (33) present the active and reactive power balance constraints for each bus, respectively. Equation (34) calculates the voltage of each bus based on the previous bus voltage and the voltage drop across the two bus connections. Equation (35) defines the maximum power generation for DG at the distribution system.
V m i n V b , t V m a x ; t , b
P ^ b , t S e l l = P ^ b , t S e l l η . i t r a n s ; t , b π i
P b , b p , t + P b , t B u y , T S O + P b , t D G + P ^ b , t b π i B u y = b p π b P b p , k , t + P b p , t + P b , t S e l l , T S O + P ^ b , t b π i S e l l ; t , b
Q b , b p , t + Q b , t T = k π b Q b p , k , t + Q b p , t ; t , b
V b p , t = V b , t r b , b p , t P b p , t + x b , b p , t Q b p , t V r e f ; t , b
P b , t D G P b , t D G , m a x ; t , b π d

2.3. Solution Approach

2.3.1. Centralized Model

In this section, without considering the possibility of negotiations between the players, under a central control structure, the social welfare is modeled in accordance with Equation (36).
M i n   ψ ψ = i = 1 N C M G i R D S O = i = 1 N t T ( ( C i , t T , C H P + C i , t T , M T + C i , t T , B o ) + P i , t B u y P i , t S e l l λ i , t ) i = 1 N t = 1 T P ^ b , t b = i S e l l P ^ b , t b = i B u y × λ i , t t = 1 T P b 1 , t S e l l , T S O λ t S e l l , T S O P b 1 , t B u y , T S O λ t B u y , T S O + b = 1 N b t = 1 T λ b D G P b , t D G = i = 1 N t = 1 T ( C i , t T , C H P + C i , t T , M T + C i , t T , B o i l e r ) t = 1 T P b 1 , t S e l l , T S O λ t S e l l , T S O P b 1 , t B u y , T S O λ t B u y , T S O + b = 1 N b t = 1 T λ b D G P b , t D G
s.t. (2)–(28) and (30)–(35)
P i , t B u y = P ^ b , t b = i S e l l ; i , t , b π i
P i , t S e l l = P ^ b , t b = i B u y ; i , t , b π i
By solving the problem in a centralized manner under Equation (36), and in consideration of coupling constraints (37) and (38), the optimal value for the social welfare function of the system can be calculated. However, by eliminating the parameters related to energy trading as global variables, it will not be possible to calculate the price of power trading for MEMGs and DS. On the other hand, due to the separate ownership concept for the MGs and data privacy issues, it will not be possible to use a centralized model for a decentralized problem structure. Therefore, in the next section, using an ADMM method, the problem will be solved in a distributed manner to maintain data privacy for market participants and calculate the flexibility indexes for MEMGs.

2.3.2. Distributed Model of Market Negotiation

Using the ADMM decomposition method, the problem is broken down into two parts which, after repeating and receiving feedback from each other, finally reach convergence and an optimal solution [40]. The ADMM solution method is used in this paper as a dynamic framework to submit the updated bids to the market in response to the received signals, which can be used to determine the effects of the flexibility of prosumers on energy-trading negotiations and the market-clearing process.
a.
The upstream problem for MEMGs
Considering the penalty function of the ADMM solution approach, the objective function of MEMG is defined as in Equation (39).
M i n t = 1 T ( C i , t T , C H P + C i , t T , M T + C i , t T , B o ) + P i , t B u y P i , t S e l l λ i , t + ( ρ 2 ) t = 1 T P ^ b , t b = i B u y P i , t S e l l + P ^ b , t b = i B u y P i , t S e l l 2 2
s.t. (2)–(28).
b.
The sub-problem on the DSO side
M a x i = 1 N t = 1 T P ^ b , t b = i S e l l P ^ b , t b = i B u y × λ i , t + t = 1 T P b 1 , t S e l l , T S O λ t S e l l , T S O P b 1 , t B u y , T S O λ t B u y , T S O b = 1 N b t = 1 T λ b D G P b , t D G ( ρ 2 ) i = 1 N t = 1 T P ^ b , t b = i B u y P i , t S e l l + P ^ b , t b = i B u y P i , t S e l l 2 2
s.t. (30)–(35).
In accordance with the Lagrange coefficient equation, the price update will be applied as follows on the DSO side:
P ^ b , t b = i = P ^ b , t b = i B u y P ^ b , t b = i > 0 , P ^ b , t b = i = P ^ b , t b = i S e l l P b , t b = i < 0
λ i , t k + 1 = λ i , t k + ρ ( P ^ b , t b = i k + 1 P i , t k + 1 )
The process of solving the method is described in the ADMM algorithm in accordance with Algorithm 1. At the first iteration, MGOs optimize their cost function with an initial assumption for energy-exchange price and power-trading suggestions to update the power-trading value. After that, the DSO receives all the suggestions from MGs; it optimizes its operation planning regarding received value. The power-trading results of DSO operation planning update the dual variable as trading price and send it back to the MGs to obtain the second feedback. This process is continued until the stopping condition is satisfied for dual variable deviations. Figure 3 shows the information-exchange procedure for DSO and MEMGs. Equation (43) presents the stopping condition of the algorithm.
λ i , t k + 1 λ i , t k < ε
Algorithm 1. Proposed distributed solution approach.
Synchronous Distributed Algorithm
1Start:K = 0, P ^ i , t = P ^ i , t 0 , λ i , t = λ i , t 0 N
2Repeat
3At each individual MEMG
4Repeat
5Wait
6       Until receive update P ^ b , t from DSO
(1) solve local problem in (39) for optimal solution P i , t
7     (2) send P i , t to the DS
8At the DSO
9Repeat
10Wait
11Until receive update P i , t from all MEMGs
 (1) solve problem (40) for optimal solution P ^ b , t
 (2) update dual variables (42):
λ i , t k + 1 = λ i , t k + ρ ( P ^ b , t b = i k + 1 P i , t k + 1 )
12(3) Send P ^ b , t , λ i , t to all microgrids
13 k k + 1
14Until a stopping criterion is met
Equation (43)

2.4. Flexibility Index

According to the definition provided for flexibility as the ability to adapt a generation or consumption pattern in response to an external signal, which would be a price signal or any other incentive or inhibitory parameter [8], two different metrics of flexibility are proposed in this paper as energy- and power-flexibility indexes. With two general points of view for energy and power-trading flexibility, it would be valuable to define specific indexes of flexibility for prosumers in accordance with market requirements and valid energy-trading deviations. In this paper, the changes in energy-trading values for the first power-trading bids, which are considered without any encouraging or preventative signals, in comparison to the last power-trading values, which are determined in response to any external signals, are considered as the flexibility of a prosumer in energy-trading negotiations. This definition can be used in many different ways and would be practical in both centralized and distributed market-clearing processes. For example, the index can be calculated by a prosumer via sensitivity analyses of energy trading in accordance with price signal changes. After the prosumer executes the sensitivity analysis, the index can be sent to the market operator as complementary data, which allows the market operator to lead the bids of prosumers in the direction of improving the social welfare of the entire system to achieve a greater margin for clearing the market under security constraints of the system. Moreover, it can be used in distributed models of market negotiations to evaluate the flexibility of prosumers. For example, if a prosumer is requested to provide a specific value of flexibility in proposed energy trading of the market, it can change the required parameter in its optimization procedure to deliver the requested power. The applications of the proposed indexes are out of the scope of the current study and would be modeled in detail in future works. In fact, the flexibility index is defined as the percentage of changes in prosumers’ bids in response to the signals received from the market. Equation (44) defines the index of flexibility for the power trading of prosumers. In this equation, P i , t t r d , 0 refers to the independent power-trading scheduling results of the prosumers in the first bid submission, and P i , t t r d , c l d is the power-trading scheduling results in accordance with the signals received from the market or a new price estimation. ADMM is investigated in this paper to validate the effects of proposed flexibility indexes on market-clearing results, and P i , t t r d , c l d is devoted to the scheduled values of energy trading after market negotiations. P i t r d , b a s e is the base power-trading value to determine the per-unit values of bids. In this study, the connection transformer capacity is considered the maximum amount of power trading between the MEMGs and DS as the base power. Considering the positive values of P i , t t r d for the sales mode and the negative values for receiving power or buying, the flexibility index will show the flexibility of the MEMG against the trading of more and less power or energy. Equation (45) indicates the index of flexibility in energy trading, which is defined as the average of the flexibilities in power trading.
I P i , t F l x = P i , t t r d , 0 P i , t t r d , c l d P i t r d , b a s e
I E i F l x = 1 T t = 1 T ( P i , t t r d , 0 P i , t t r d , c l d P i , t t r d , b a s e )

3. Numerical Results

3.1. Case Study

In this study, a modified IEEE 33 bus power distribution system with four different MEMGs with integrated electrical, heat, and combined demand is adopted to demonstrate the proposed model. Figure 4 shows the case study structure for the proposed model. Table 2 displays the equipment combination in each of the MEMGs. The maximum 5% is considered for a changeable load. The electrical, heat, and combined demand profiles of each MG are shown in Figure 5. Figure 6 also shows the active and reactive load profiles of the distribution network. The characteristics of the system’s components are presented in Table 3. By examining the results of the social welfare cost function in two centralized and distributed control models, the results for the proposed model can be validated. The simulation results are provided by General Algebraic Modeling System (GAMS) software, Version 25.1.2 under a home computer with Core i7 and 8 GB RAM.

3.2. Simulation Result

3.2.1. Power Flexibility Index Analysis

Figure 7 shows the power flexibility indexes for each MG. According to the obtained results, under the proposed DRP model, I P i , t F l x in the direction of selling more power have generally increased. By shifting the load from the peak to the off-peak duration, the flexibility in the direction of providing power would increase. On the other hand, on shifting the load to the off-peak hour, the flexibility of off-peak durations in the direction of receiving more power increased. I P i , t F l x indirectly reflects the deficiency or adequacy of power for each hour. At hour 1 for MEMG1, the I P 1,1 F l x is reported as 19%, which is changed to −37% after DRP employment. For a peak hour, such as hour 16, the I P 1,16 F l x is changed from 51% to −8%, which shows a decrease in power requests in accordance with the change in demand to the thermal kind. However, the flexibility of selling power is still lesser for peak hours in comparison to off-peak hours. As an example, electrical and thermal demands were shifted from peak hours (13–14) to off-peak hours such as (5–6) and 23. This scheduling reduced the demand values in hours (13–14) and enabled MEMG1 to be more flexible in power delivery. On the other hand, with the shift of load to duration (5–6), the energy demand increased and caused MEMG1 to agree to receive more power, which creates more flexibility in receiving power. The proposed DRP model would also change the flexibility of MGs in buying energy. In MEMG2, the DRP is considered in the form of electric to thermal load change in the scheduling horizon. Thus, under the reduction of demand for power, MEMG2 showed less flexibility for receiving power in hours (1–11). However, during peak hours (15–16) and hour 23, when the price of power is low, the flexibility in receiving power increased. MEMG3 without any energy storage used the DRP as an opportunity for energy arbitrage. Results show that the flexibility of MEMG3 in buying power at off-peak hours and selling power at peak hours had already increased. MEMG4 does not include any HP unit, so the DRP model was executed in this unit to transfer the thermal and electric power from peak durations. So, with the shift of electrical and thermal demand from hour 11 to 6, the flexibility of MEMG4 in receiving power increased in the off-peak hours. An analysis of the results for the I P i , t F l x indicates that these indexes can be used in the system for demand management, and create an opportunity for market operators to enhance the system efficiency under a wider bid range of prosumers.

3.2.2. Energy-Flexibility Index Analysis

The energy-flexibility index includes a quantitative statement for the changes in the energy-trading schedule of the MEMGs compared to the initial programmed values. Without the DRPs, MEMG1 received 516.64 kWh of energy less than the initial scheduled value. Employment of DRPs changes the state to receive 206,250 kWh of energy more than the scheduled value. Employing the DRP in the MEMG1 caused the MEMG to shift the heat and power loads from peak hours to off-peak hours. This new feature of MEMG made it economic to buy more power in energy-trading negotiations. This is because under a more flexible system, MEMG has more jurisdiction to use market opportunities. So, according to the employment of DRP as flexibility resources, the energy FI changed from −7.18% to 28.65%, which indicates achieving more energy in comparison to the scheduled value. The flexibility of MEMG2 in receiving energy increased from −4.73% to −5.10% with DRP. These negative FIs mean MEMG2 agreed to achieve 23 kWh energy less underemployment of DRP (340.85 kWh in comparison to 367.23 kWh). Similarly, the energy FI increased from −5.82% to −6.20% under DRPs for MEMG3. This 0.38% difference is equivalent to an agreement on achieving 27.82 kWh of energy less than the state of no DRP employment (418.69 kWh in comparison to 446.51 kWh). Table 4 presents the energy FIs for each MEMG with and without DRPs. The results indicate that employing the proposed DRP model would enhance the energy-flexibility index of MEMGs. Although the proposed DRP model was the same for all of MEMGs, the results show that according to the different energy equipment of MGs, the results of flexibility changes would be different. In this study, MEMG1 with the greatest diversity of resources, experienced the greatest flexibility improvement.

3.2.3. Proposed Shiftable and Changeable DRP Model

Figure 8 illustrates the results of the proposed DRP model in different MEMGs. With the obtained results, it can be concluded that the employment of different DRPs is integrated with MEMG resources and price parameters. All of the results indicate that the DRP model tried to shift and change the energy demand to a more economic type, and shifted the demand from peak hours and more expensive energy carriers to the cheaper kind. To analyze the effect of the proposed DRP model on the operating costs, the related results are presented in Table 5. As can be seen in the results, the use of DRP models improved the social welfare function by 2.75%, which shows the effectiveness of the proposed model in reducing system-operating costs. Under demand-response programming, the operating costs of MEMGs decreased between 6.28 and 10.98 percent. The proposed DRP also reduced the costs of operation for DSO by 0.09% under the reduction of the market-clearing prices. Due to the increase in flexibility on the MG side, DRPs demonstrated a prominent effect on the MG side. However, the effects of increasing flexibility in social welfare improvement and reducing DSO costs have also been reported.

3.2.4. Convergency and Flexibility

This section discusses the impact of flexibility on the convergence process of a market. To analyze the results, hours 6, 14, and 16 for the negotiations between MEMG4 and DS were studied. Figure 9 shows the negotiations for the selected hours of MEMG4. Hour 6 with a flexibility index of 32.16% has the highest index in the direction of receiving energy. By examining the convergence process of this MG at 6 o’clock, it can be seen that this MG reached an agreement with DSO after only 11 iterations to receive 230.063 kW power instead of the scheduled 133.593 kWh of energy. In contrast, hour 14 with a flexibility index of 16.67 for the sale of energy could reach an agreement with DSO after 499 iterations. This indicates the role of the player’s flexibility on the speed of convergence. At hour 16, MEMG4 had no flexibility and emphasized the sale of 152.52 kWh of energy to DSO. In this regard, DSO had to employ more flexibility to agree with this power purchase and provide the balance of a transaction to make an agreement. The results show that, if a player is not flexible at all, the opposite player in the transaction needs to employ more flexibility if it is essential to reach an agreement. On the other hand, more flexibility would increase the possibility of reaching a common solution and would improve the convergence speed of the problem. This is because there are more possibilities and strategies to reach an agreement. In this study, the players are considered rational without any intrigue in abusing the flexibility while employing the ADMM method for modeling the negotiations. Moreover, in the ADMM method, the number of iterations continues until all transactions reach convergence, and the results may change until the global optimal solution is determined. So, the number of iterations mentioned in this section to consider the required iteration of negotiation is related to the first stable solution without any changes.

3.2.5. Operation Scheduling and Market-Clearing Results

Figure 10 shows the electrical and thermal power balance constraints for each MEMG. Employing the DRP in the MEMG1 causes the CHP unit to allocate more capacity to generate thermal energy for the duration of (1–6) and change the combined load to the heat kind. MEMG2 is considered without any CHP unit. In this MG, the MT generates electrical energy at maximum capacity. Due to the lack of CHP units, energy converter equipment and storage devices play a fundamental role in supporting DRP. MEMG3 does not include EES, so it cannot use the energy arbitrage revenue. Despite the excess energy generation in this MG, MEMG3 had to purchase power in the peak hours (19–22) with high power prices. The use of changeable load models reduced the amount of power purchased in peak hours (13–16) by about 18 kW, and in the period (20–23) by about 24 kW. MEMG4 does not have any energy converter equipment or a heat pump unit, and thus the role of TES devices during peak heat demand hours is very significant. As can be deduced from the results, MEMG1 with the greatest diversity of resources and the highest flexibility index, and MEMG4 with the lowest index value due to the lack of an HP, made maximum use of shifting feature of the load. In this study, DSO operation scheduling is related to determining the amount of generation for each of the DGs and scheduling the purchase and sale of energy from TSO and MGs. With Figure 11, it can be deduced that in times of peak energy prices, DGs generated power with their maximum capacity, so that in addition to buying energy from MGs and internal generation, they can trade the excess energy purchased in the upstream market. Due to the higher marginal cost of DG2 compared to 1 (8 ¢/kWh and 6 ¢/kWh), DG1 generated more energy in comparison to DG2. Figure 12 shows the market-clearing prices for power trading. The energy requirement of the DSO is mostly provided through buying energy from the upstream market, and using the DG capacity for providing power is not an effective solution with respect to the market prices. This is because, in general, the market-clearing price is less than 6 ¢/kWh, which creates a more economical solution for DSO to trade energy with MGs. It is clear why DSO shows more flexibility in energy-trading negotiations in purchasing power from MGs.

4. Conclusions

In this paper, a new flexibility analysis for MEMGs is executed in accordance with bilateral energy transactions in a local energy market. In this regard, a new model of DRP is proposed to enhance the flexibility of MEMG operation and analyze the effects of participants’ flexibility improvement in the market-clearing process and social welfare enhancement. Moreover, two new metrics for flexibility indexes are proposed. The simulation results indicate that the proposed DRP model increased the flexibility index in buying and selling energy by 0.81% and 35.82%, respectively, which improved the social welfare value of the system by 2.75%. Moreover, the results of the iteration number for the convergency process indicate that, with more flexible market participants, the market-clearing process would be easier under more adaptable negotiators. This concept can be used as an option for market participants in market payments to obtain more financial benefits for the participants who submit bids with a flexibility index for energy transactions. This study showed that using the flexibility of prosumers in market clearing would enhance the social welfare of the entire system; this can be employed in new market structures. For future studies, it is recommended that the effects of the flexibility of prosumers on the security margins of the system also be studied.

Author Contributions

Conceptualization, A.S.; Methodology, A.S., S.J. and M.N.-H.; Software, A.S.; Validation, S.J. and M.N.-H.; Formal analysis, S.J. and M.N.-H.; Investigation, A.S.; Resources, A.S. and S.J.; Writing—original draft, A.S.; Writing—review & editing, S.J. and M.N.-H.; Visualization, A.S.; Supervision, S.J. and M.N.-H.; Project administration, S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data would be available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Indexes
iIndexes of MEMGs in distribution system
tIndexes of time periods, hour
b/bpIndexes of buses
π i Sets of connections for MEMGs and buses
π d Sets of buses that are connected to DG
π Sets of all branches that are between bus bp and k
Scalars
MA big number
ρPenalty coefficient
Parameter
P i , t M G , Q i , t M G , H i , t M G Active / Reactive / Heat   demand   of   M E M G i at time t, kW, kvar, kW
C h i , t M G , m a x Maximum   changeable   load   of   M E M G i at time t, kW
P b , t , Q b , t Active/Reactive power consumed at bus b of distribution system at time t, kW, kvar
r b , b p , x b , b p Resistance/Reactance of a line connecting bus b and bp, ohm
V m a x , V m i n Maximum/Minimum Voltage limit, kV
λ t S e l l , T S O , λ t B u y , T S O Forecasted TSO price for energy trading with DSO at time t, ¢/kW·h
λ t N G Forecasted Natural Gas price at time t, ¢/Btu.h
λ b D G Cost of energy generation for DG of bus b (¢/kW·h)
π i C H P , S U , π i C H P , S D Start-Up/ Shout-down cost of CHP unit i, ¢
η i C H P Efficiency of CHP unit i,%
P i , A C H P / P i , B C H P / P i , C C H P / P i , D C H P
H i , A C H P / H i , B C H P / H i , C C H P / H i , D C H P
Feasible operation region of a CHP plant, kW
P i M T , m a x Maximum electrical capacity of MT unit i, kW
η i M T Efficiency of MT unit i, %
P i C h , E E S , m a x , P i d C h , E E S , m a x Maximum Charge/Discharge power of EES, kW
S O C i , t E E S , m a x , S O C i , t E E S , m i n Maximum/Minimum Electric quantity of EES during time t, kW·h
λ i C h , E E S , λ i d C h , E E S Charge/Discharge efficiency of EES i, %
H i B O , m a x Maximum thermal capacity power of boiler i, kW
λ i B o Efficiency of Boiler i, %
H i H P , m a x Maximum thermal power of Heat Pump i, kW
η i H P Efficiency of Heat Pump i, %
H i C h , T E S , m a x , H i D C h , T E S , m a x Maximum Charge/Discharge thermal power of TES, kW
η i C h , T E S , η i d C h , T E S Charge/Discharge efficiency of TES unit, %
I P i , t F l x Index of flexibility for power trading of MEMG i at time t
I E i F l x Index of flexibility for energy trading of MEMG i
Positive Variables
P b 1 , t S e l l , T S O , P b 1 , t B u y , T S O Electrical power of DSO for power trading with TSO at time t, kW
P i , t C H P Active power generation of CHP unit i at time t, kW
H i , t C H P Heat power generation of CHP unit i at time t, kW
C i , t C H P , S U , C i , t C H P , S D Start-Up/ Shout-down cost of CHP unit i at time t, ¢
C i , t C H P Cost of CHP unit i at time t, ¢
C i T , C H P Total cost of CHP unit i, ¢
P i , t M T Active power generation of MT unit i at time t, kW
C i , t N G Cost of MT unit i at time t, ¢
C i T , N G Total cost of MT unit i, ¢
P i , t C h , E E S , P i , t d c h , E E S Charge/Discharge power of EES unit i at time t, kW
S O C i , t E E S Electric quantity of EES i at time t, kW·h
Heat power generation of Boiler i at time t, kW
H i , t B o
C i , t B o Cost of Boiler i at time t, ¢
C i T , B O Total cost of Boiler i, ¢
H i , t H P Heat power generation of Heat Pump i at time t, kW
P i , t H P Electrical power consumed by Heat Pump i at time t, kW
H i , t C h , T E S , H i , t d c h , T E S Charge/Discharge thermal power of TES unit i at time t, kW
S O C i , t T E S Electric quantity of TES during time t, kW·h
V b , t Voltage magnitude at bus b at time t, kV
P b , t D G Power generation of DG unit of bus b at time t, kW
C D S O D G Cost of DG generations for DSO, ¢
Free Variables
e c h i , t , h c h i , t Changeable   load   of   M E M G i at time t, kW·h
Q b , t T Reactive power trading of DSO with TSO at bus b and time t, kvar
P b , b p , t , Q b , b p , t Active/Reactive power flowing between branch of bus b and bp at time t, kW/kvar
R D S O M G Revenue of energy trading with MGs for DSO, ¢
R D S O T S O Revenue of energy trading with TSO for DSO, ¢
R D S O Total revenue of DSO in operation planning, ¢
R i , t Revenue of M E M G i of energy trading with DSO at time t, ¢
C M G i Total   operation   cos t   of   M E M G i in operation planning, ¢
Binary Variable
x i , t C H P Binary variable for CHP unit i at time t (0 for off status and 1 for on status)
x i , t E E S , x i , t T E S Binary variable for EES/TES unit i at time t (1 for charging status and 0 for discharging status)
Global variable
λ i , t Electricity   Price   for   power   trading   of   DSO   and M E M G i at time t, ¢/kW·h
P i , t S e l l , P i , t B u y Electrical   power - trading   schedule   of M E M G i at time t, kW
P ^ b , t S e l l , P ^ b , t B u y Electrical power-trading schedule of DSO at time t at bus b π i , kW

References

  1. Wu, R.; Sansavini, G. Active Distribution Networks or Microgrids? Optimal Design of Resilient and Flexible Distribution Grids with Energy Service Provision. Sustain. Energy Grids Netw. 2021, 26, 100461. [Google Scholar] [CrossRef]
  2. Sahebi, A.; Jadid, S. A Robust Model of Local Energy Market under a Security Constraint-Based Approach for Distribution System Operator and Multi-Energy Microgrids. Electr. Power Syst. Res. 2023, 217, 109164. [Google Scholar] [CrossRef]
  3. Papazoglou, G.K.; Forouli, A.A.; Bakirtzis, E.A.; Biskas, P.N.; Bakirtzis, A.G. Day-Ahead Local Flexibility Market for Active and Reactive Power with Linearized Network Constraints. Electr. Power Syst. Res. 2022, 212, 108317. [Google Scholar] [CrossRef]
  4. Ma, J.; Silva, V.; Belhomme, R.; Kirschen, D.S.; Ochoa, L.F. Evaluating and Planning Flexibility in Sustainable Power Systems. IEEE Trans. Sustain. Energy 2013, 4, 200–209. [Google Scholar] [CrossRef] [Green Version]
  5. Oree, V.; Sayed Hassen, S.Z. A Composite Metric for Assessing Flexibility Available in Conventional Generators of Power Systems. Appl. Energy 2016, 177, 683–691. [Google Scholar] [CrossRef]
  6. Mouton, D.; Mararakanye, N.; Bekker, B. Assessment of Methods to Measure Power System Flexibility: A Review. In Proceedings of the 2021 56th International Universities Power Engineering Conference: Powering Net Zero Emissions, UPEC 2021, Middlesbrough, UK, 31 August–3 September 2021. [Google Scholar]
  7. Ulbig, A.; Andersson, G. Analyzing Operational Flexibility of Electric Power Systems. Int. J. Electr. Power Energy Syst. 2015, 72, 155–164. [Google Scholar] [CrossRef]
  8. Minniti, S.; Haque, N.; Nguyen, P.; Pemen, G. Local Markets for Flexibility Trading: Key Stages and Enablers. Energies 2018, 11, 3074. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, Z.; Hu, J.; Ai, X.; Wu, J.; Yang, G. Transactive Energy Supported Economic Operation for Multi-Energy Complementary Microgrids. IEEE Trans. Smart Grid 2021, 12, 4–17. [Google Scholar] [CrossRef]
  10. Chen, J.J.; Qi, B.X.; Rong, Z.K.; Peng, K.; Zhao, Y.L.; Zhang, X.H. Multi-Energy Coordinated Microgrid Scheduling with Integrated Demand Response for Flexibility Improvement. Energy 2021, 217, 119387. [Google Scholar] [CrossRef]
  11. Löschenbrand, M. A Transmission Expansion Model for Dynamic Operation of Flexible Demand. Int. J. Electr. Power Energy Syst. 2021, 124, 106252. [Google Scholar] [CrossRef]
  12. Shekari, T.; Gholami, A.; Aminifar, F. Optimal Energy Management in Multi-Carrier Microgrids: An MILP Approach. J. Mod. Power Syst. Clean Energy 2019, 7, 876–886. [Google Scholar] [CrossRef] [Green Version]
  13. Karimi, H.; Jadid, S. Optimal Energy Management for Multi-Microgrid Considering Demand Response Programs: A Stochastic Multi-Objective Framework. Energy 2020, 195, 116992. [Google Scholar] [CrossRef]
  14. Shen, Z.; Wu, L.; Sedgh, S.A.; Radmehr, F. Optimal Operation Scheduling of a Microgrid Using a Novel Scenario-Based Robust Approach. Eng. Optim. 2023, 22. [Google Scholar] [CrossRef]
  15. Bolfek, M.; Capuder, T. An Analysis of Optimal Power Flow Based Formulations Regarding DSO-TSO Flexibility Provision. Int. J. Electr. Power Energy Syst. 2021, 131, 106935. [Google Scholar] [CrossRef]
  16. Vijay, R.; Mathuria, P. Complex Power Flexibility Evaluation Using Energy Arbitrage between Transmission and Distribution. Electr. Power Syst. Res. 2022, 203, 107641. [Google Scholar] [CrossRef]
  17. Kumamoto, T.; Aki, H.; Ishida, M. Provision of Grid Flexibility by Distributed Energy Resources in Residential Dwellings Using Time-of-Use Pricing. Sustain. Energy Grids Netw. 2020, 23, 100385. [Google Scholar] [CrossRef]
  18. Yi, Z.; Xu, Y.; Wu, C. Improving Operational Flexibility of Combined Heat and Power System through Numerous Thermal Controllable Residents Aggregation. Int. J. Electr. Power Energy Syst. 2021, 130, 106841. [Google Scholar] [CrossRef]
  19. Babagheibi, M.; Jadid, S.; Kazemi, A. Distribution Locational Marginal Pricing for Congestion Management of an Active Distribution System with Renewable-Based Microgrids under a Privacy-Preserving Market Clearing Approach and Load Models. Sustain. Energy Grids Netw. 2022, 32, 100935. [Google Scholar] [CrossRef]
  20. Bahramara, S.; Parsa Moghaddam, M.; Haghifam, M.R. A Bi-Level Optimization Model for Operation of Distribution Networks with Micro-Grids. Int. J. Electr. Power Energy Syst. 2016, 82, 169–178. [Google Scholar] [CrossRef]
  21. Zhang, H.; Zhang, S.; Hu, X.; Cheng, H.; Gu, Q.; Du, M. Parametric Optimization-Based Peer-to-Peer Energy Trading among Commercial Buildings Considering Multiple Energy Conversion. Appl. Energy 2022, 306, 118040. [Google Scholar] [CrossRef]
  22. Wang, N.; Liu, Z.; Heijnen, P.; Warnier, M. A Peer-to-Peer Market Mechanism Incorporating Multi-Energy Coupling and Cooperative Behaviors. Appl. Energy 2022, 311, 118572. [Google Scholar] [CrossRef]
  23. Esmat, A.; de Vos, M.; Ghiassi-Farrokhfal, Y.; Palensky, P.; Epema, D. A Novel Decentralized Platform for Peer-to-Peer Energy Trading Market with Blockchain Technology. Appl. Energy 2021, 282, 116123. [Google Scholar] [CrossRef]
  24. Nguyen, D.H.; Ishihara, T. Distributed Peer-to-Peer Energy Trading for Residential Fuel Cell Combined Heat and Power Systems. Int. J. Electr. Power Energy Syst. 2021, 125, 106533. [Google Scholar] [CrossRef]
  25. Babagheibi, M.; Sahebi, A.; Jadid, S.; Kazemi, A. An Integrated Design of Heat and Power Market for Energy Hubs Considering the Security Constraints of the System. Sustain. Cities Soc. 2023, 96, 104616. [Google Scholar] [CrossRef]
  26. Zhao, J.; Zheng, T.; Litvinov, E. A Unified Framework for Defining and Measuring Flexibility in Power System. IEEE Trans. Power Syst. 2016, 31, 339–347. [Google Scholar] [CrossRef]
  27. Tuohy, A.; Lannoye, E. Metrics for Quantifying Flexibility in Power System Planning; Electric Power Research Institute: Palo Alto, CA, USA, 2014. [Google Scholar]
  28. Krommydas, K.F.; Dikaiakos, C.N.; Papaioannou, G.P.; Stratigakos, A.C. Flexibility Study of the Greek Power System Using a Stochastic Programming Approach for Estimating Reserve Requirements. Electr. Power Syst. Res. 2022, 213, 108620. [Google Scholar] [CrossRef]
  29. Gusain, D.; Cvetković, M.; Palensky, P. Quantification of Operational Flexibility from a Portfolio of Flexible Energy Resources. Int. J. Electr. Power Energy Syst. 2022, 141, 107466. [Google Scholar] [CrossRef]
  30. Karimi, H.; Gharehpetian, G.B.; Ahmadiahangar, R.; Rosin, A. Optimal Energy Management of Grid-Connected Multi-Microgrid Systems Considering Demand-Side Flexibility: A Two-Stage Multi-Objective Approach. Electr. Power Syst. Res. 2023, 214, 108902. [Google Scholar] [CrossRef]
  31. Gasca, M.V.; Ibáñez, F.; Pozo, D. Flexibility Quantification of Thermostatically Controlled Loads for Demand Response Applications. Electr. Power Syst. Res. 2022, 202, 107592. [Google Scholar] [CrossRef]
  32. Zhou, Y.; Cao, S. Quantification of Energy Flexibility of Residential Net-Zero-Energy Buildings Involved with Dynamic Operations of Hybrid Energy Storages and Diversified Energy Conversion Strategies. Sustain. Energy Grids Netw. 2020, 21, 100304. [Google Scholar] [CrossRef]
  33. Alipour, M.; Gharehpetian, G.B.; Ahmadiahangar, R.; Rosin, A.; Kilter, J. Energy Storage Facilities Impact on Flexibility of Active Distribution Networks: Stochastic Approach. Electr. Power Syst. Res. 2022, 213, 108645. [Google Scholar] [CrossRef]
  34. Wu, Y.K.; Tan, W.S.; Huang, S.R.; Chiang, Y.S.; Chiu, C.P.; Su, C.L. Impact of Generation Flexibility on the Operating Costs of the Taiwan Power System under a High Penetration of Renewable Power. IEEE Trans. Ind. Appl. 2020, 56, 2348–2359. [Google Scholar] [CrossRef]
  35. Mohandes, B.; El Moursi, M.S.; Hatziargyriou, N.; Khatib, S. El A Review of Power System Flexibility with High Penetration of Renewables. IEEE Trans. Power Syst. 2019, 34, 3140–3155. [Google Scholar] [CrossRef]
  36. Heydarian-Forushani, E.; Golshan, M.E.H. Quantitative Flexibility Assessment of a Comprehensive Set of Demand Response Programs. Int. J. Electr. Power Energy Syst. 2020, 116, 105562. [Google Scholar] [CrossRef]
  37. Nazari-Heris, M.; Abapour, S.; Mohammadi-Ivatloo, B. Optimal Economic Dispatch of FC-CHP Based Heat and Power Micro-Grids. Appl. Therm. Eng. 2017, 114, 756–769. [Google Scholar] [CrossRef]
  38. Hemmati, M.; Abapour, M.; Mohammadi-Ivatloo, B.; Anvari-Moghaddam, A. Optimal Operation of Integrated Electrical and Natural Gas Networks with a Focus on Distributed Energy Hub Systems. Sustainability 2020, 12, 8320. [Google Scholar] [CrossRef]
  39. Li, R.; Wei, W.; Mei, S.; Hu, Q.; Wu, Q. Participation of an Energy Hub in Electricity and Heat Distribution Markets: An MPEC Approach. IEEE Trans. Smart Grid 2019, 10, 3641–3653. [Google Scholar] [CrossRef]
  40. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers. Found. Trends Mach. Learn. 2010, 3, 1–122. [Google Scholar] [CrossRef]
Figure 1. MEMG components.
Figure 1. MEMG components.
Sustainability 15 09985 g001
Figure 2. A typical radial network.
Figure 2. A typical radial network.
Sustainability 15 09985 g002
Figure 3. The information exchange of DSO and MEMG.
Figure 3. The information exchange of DSO and MEMG.
Sustainability 15 09985 g003
Figure 4. IEEE 33 bus test system with four connected MEMGs and DGs.
Figure 4. IEEE 33 bus test system with four connected MEMGs and DGs.
Sustainability 15 09985 g004
Figure 5. Load profile of each MEMG.
Figure 5. Load profile of each MEMG.
Sustainability 15 09985 g005
Figure 6. Load profile of distribution system.
Figure 6. Load profile of distribution system.
Sustainability 15 09985 g006
Figure 7. Flexibility indexes for MEMGs.
Figure 7. Flexibility indexes for MEMGs.
Sustainability 15 09985 g007
Figure 8. DRP scheduling for MEMGs.
Figure 8. DRP scheduling for MEMGs.
Sustainability 15 09985 g008
Figure 9. Market process negotiations for MEMG4 at hours 6, 14, and 16.
Figure 9. Market process negotiations for MEMG4 at hours 6, 14, and 16.
Sustainability 15 09985 g009
Figure 10. Power balance constraint.
Figure 10. Power balance constraint.
Sustainability 15 09985 g010
Figure 11. Operation scheduling of DSO.
Figure 11. Operation scheduling of DSO.
Sustainability 15 09985 g011
Figure 12. Market-clearing prices.
Figure 12. Market-clearing prices.
Sustainability 15 09985 g012
Table 1. Literature review.
Table 1. Literature review.
Index FeaturePropertiesReasonEffects
[4]Power capacityAvailable power and ramp rateVariability and uncertaintyComparison of generators
[5]Technical weighted IndexOperational range, ramp rate, start-up and shut-down times, minimum up and down timesFrequent and larger variations in the net loadComparison of generators
[26]Uncertainty coverage capabilityAvailable powerUncertaintyMeasure the flexibility of power system
[27]Expected unserved rampingUnserved rampingVariability and uncertaintyRisk management
[28]Insufficient ramping resource expectation, loss of load expectation, flexibility residual indexesRamp rate, loss of loadLoad forecast errorsReserve requirements
[29]Unserved flexibility serviceUnserved energy, duration of insufficient flexibilityRenewable generation integrationDesign appropriate DRP, formulate operational policies
[30]Flexibility during peak periodEnergy flexibility of the demand-response programsUncoordinated DRPDetermines the best strategy for operation planning
[32]Electricity-providing capability, time durations of servicesForced flexibility factor, delayed flexibility factorAdvanced energy conversionsDemand side management, control strategy
[35]Response timeTime durationGeneration outageContingency analysis
C.PRescheduling capabilityEnergy-trading changes/power-trading changesEvaluation of prosumers’ adoptionComparison of market participants, market-clearing development
Table 2. Component equipment of each Multi-Carrier Microgrid.
Table 2. Component equipment of each Multi-Carrier Microgrid.
MEMGComponent
CHPMTHPBoilerTESEES
M E M G 1
M E M G 2
M E M G 3
M E M G 4
Table 3. Characteristics of the systems’ components.
Table 3. Characteristics of the systems’ components.
CHP P i , A C H P = 250 / P i , B C H P = 217.6 / P i , C C H P = 81.98 / P i , D C H P = 100 H i , A C H P = 0 / H i , B C H P = 180 / H i , C C H P = 104.8 / H i , D C H P = 0 η i C H P = 0.9
MT P i M T , m a x = 200 η i M T = 0.9
EES S O C i E E S , m a x = 300 P i d C h , E E S , m a x = 50 η i C h , E E S = 0.9 / η i d C h , E E S = 0.8
TES S O C i T E S , m a x = 200 P i d C h , T E S , m a x = 100 / P i C h , T E S , m a x = 100 η i C h , T E S = 0.9 / η i d C h , T E S = 0.8
Bo H i B o , m a x = 45 η i B o = 0.9
HP H i H P , m a x = 135 α = 2
DG P i D G , m a x = 900 λ b 8 D G = 6 / λ b 8 D G = 8
TR η i t r a n s = 0.9
Table 4. Energy-flexibility index.
Table 4. Energy-flexibility index.
Value ChangesWith DRPWithout DRP
MEMG135.83%28.65%−7.18%
MEMG2−0.37%−5.10%−4.73%
MEMG3−0.38%−6.20%−5.82%
MEMG4−0.80%−3.23%−2.43%
Average Index Values9.34%10.79%5.04%
Table 5. Cost analysis under DRP.
Table 5. Cost analysis under DRP.
ImprovementWithout DRPWith DRP
DSO0.09%133,305133,189
MEMG17.47%24,08722,414
MEMG26.55%13,93113,075
MEMG36.28%35,93433,810
MEMG410.98%10,7059646
Total Cost2.75%217,964212,133
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

Sahebi, A.; Jadid, S.; Nazari-Heris, M. Flexibility Analysis for Multi-Energy Microgrid and Distribution System Operator under a Distributed Local Energy Market Framework. Sustainability 2023, 15, 9985. https://doi.org/10.3390/su15139985

AMA Style

Sahebi A, Jadid S, Nazari-Heris M. Flexibility Analysis for Multi-Energy Microgrid and Distribution System Operator under a Distributed Local Energy Market Framework. Sustainability. 2023; 15(13):9985. https://doi.org/10.3390/su15139985

Chicago/Turabian Style

Sahebi, Ali, Shahram Jadid, and Morteza Nazari-Heris. 2023. "Flexibility Analysis for Multi-Energy Microgrid and Distribution System Operator under a Distributed Local Energy Market Framework" Sustainability 15, no. 13: 9985. https://doi.org/10.3390/su15139985

APA Style

Sahebi, A., Jadid, S., & Nazari-Heris, M. (2023). Flexibility Analysis for Multi-Energy Microgrid and Distribution System Operator under a Distributed Local Energy Market Framework. Sustainability, 15(13), 9985. https://doi.org/10.3390/su15139985

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