Next Article in Journal
Tyrosinase/Chitosan/Reduced Graphene Oxide Modified Screen-Printed Carbon Electrode for Sensitive and Interference-Free Detection of Dopamine
Next Article in Special Issue
Optimal Non-Integer Sliding Mode Control for Frequency Regulation in Stand-Alone Modern Power Grids
Previous Article in Journal
High-Speed Photodetectors for Microwave Photonics
Previous Article in Special Issue
Optimal Operation of Isolated Microgrids Considering Frequency Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hierarchical Optimization Method for Energy Scheduling of Multiple Microgrids

1
School of Computer Science and Technology, Anhui University, Hefei 230601, China
2
Engineering Research Center of Power Quality, Ministry of Education, Anhui University, Hefei 230601, China
3
Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne 3122, Australia
4
State Grid Anhui Electric Power Co. Ltd. Research Institute, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(4), 624; https://doi.org/10.3390/app9040624
Submission received: 29 December 2018 / Revised: 10 February 2019 / Accepted: 11 February 2019 / Published: 13 February 2019
(This article belongs to the Special Issue Microgrids)

Abstract

:
This paper proposes a hierarchical optimization method for the energy scheduling of multiple microgrids (MMGs) in the distribution network of power grids. An energy market operator (EMO) is constructed to regulate energy storage systems (ESSs) and load demands in MMGs. The optimization process is divided into two stages. In the first stage, each MG optimizes the scheduling of its own ESS within a rolling horizon control framework based on a long-term forecast of the local photovoltaic (PV) output, the local load demand and the price sent by the EMO. In the second stage, the EMO establishes an internal price incentive mechanism to maximize its own profits based on the load demand of each MG. The optimization problems in these two stages are solved using mixed integer programming (MIP) and Stackelberg game theory, respectively. Simulation results verified the effectiveness of the proposed method in terms of the promotion of energy trading and improvement of economic benefits of MMGs.

1. Introduction

With the increasing penetration of renewable energy resources (RESs) in the distribution networks of regional power grids, intermittent RESs, which are connected to distribution networks in a distributed way with a small capacity and high density, have great influence on the stability of regional power grids. The energy management and control of distributed generation, load demand and energy storage systems (ESSs) using microgrid (MG) technology can effectively improve the stability of regional power grids [1,2,3]. MGs with different characteristics coexist in distribution networks of regional power grids and form multiple microgrids (MMGs). MGs with different characteristics can achieve energy interaction through energy management systems (EMSs), which not only can enhance the reliability of the regional power supply [4], but also promote the utilization of renewable energy, and improve the economic benefits of MGs [5].
In order to manage the energy scheduling between different MGs, centralized coordination control is commonly adopted in MMG architecture [6,7,8]. In [9], the energy market operator (EMO) acts as a centralized control system to coordinate the energy iteration between a cluster of selling MGs and a cluster of buying MGs. A two-stage robust optimization for energy transactions in MMGs is proposed in [10], which could minimize system cost under the worst realization of uncertain PV output. In [11], a multiple agent system (MAS)-based hierarchical energy management strategy for MMGs is proposed with easy implementation and low computation cost. A practical model is proposed for distribution companies to minimize the total operation cost of the system including distribution networks and MGs through coordinated operation [12]. However, the above optimization methods are all aimed at a situation in which there is no direct conflict of interest between MGs and upper manager. And these methods are not very suitable for competitive hierarchical power market structures [13].
To stimulate the energy transaction potential of MGs, a multi-market participation framework is proposed for distribution network operators (DNOs) [14]. Since market participants belong to different stakeholders, DNOs and MGs can trade energy not only through cooperation, but also through competition. The relationship among participants is suitable to be solved by game theory. In [15], an incentive mechanism based on cooperative game theory is proposed to reduce the peak ramp of distribution networks and improve the benefits of MGs. In [16], a two-level day-ahead scheduling structure based on Stackelberg game is proposed to stimulate MGs participating in power sale bidding, which could reduce the total cost of the DNO. A bi-level programming based on Stackelberg is adopted by the DNO to improve the benefits of MGs while the cost of the DNO reaches the minimum [17]. However, these methods are based on day-head energy market, and there are still some shortcomings in dealing with the uncertainties of renewable energy generation and load demand.
In this paper, a two-level optimization method is proposed for hour-ahead MMG energy scheduling in distribution network energy markets, where the EMO is the upper manager of the whole MMG system, and the MGO is lower manager of the local MG. This method can be implemented in two stages. According to the short-term forecast information, the local MGO optimizes energy scheduling of storage units by adopting rolling optimization in the first stage. The energy transactions between the EMO and MGs are optimized in the second stage, which is modeled as a Stackelberg game. The EMO is the leader of the game who determines the prices of the next hour to maximize its own utility, while the MGs are followers who respond to prices by adjusting local load demands. Through two stages of optimization, the economic utilities of the EMO and MGOs are both improved. The other significant feature of this proposed method is that energy storage and load demands are dispatched hourly, which makes it more reasonable to manage energy trading in MMGs under the uncertainties of PV output and load demands.

2. Framework

2.1. Structure of an MMG

The structure of an MMG in this paper is shown in Figure 1. The parties involved in energy trading in the MMG are MGs, the EMO and a power grid. Each MG mainly consists of a PV system, an ESS, loads, a smart meter and a local EMS. For a rational MG, the first choice of PV generation is to supply its own loads and charge its own ESS. The MG can then act as a seller when there is an energy surplus or as a buyer when there is an energy deficit. The net energy profile of the MG is optimized by the local EMS according to the benefits generated by energy consumption and the internal prices from the EMO. The EMO is responsible for stimulating energy trading in MMGs by establishing a reasonable internal buying price ( p c b h ) and selling price ( p c s h ) each hour. At the occurrence of an internal energy mismatch, the EMO trades with power grids to balance supply and demand. In order to guarantee that the profits of the MG obtained from energy sharing among MGs are better than those of the MG obtained from energy trading with power grids directly, the internal prices produced by the EMO should be between the selling price ( p g s h ) and the buying price ( p g b h ) of power girds, allowing the EMO to maximize its own profits under the constraint as follows:
p g s h p c s h < p c b h p g b h

2.2. Operation Strategy

Figure 2 shows the details of the operation strategy of an MMG, which can be described by two stages of optimization.

2.2.1. First Stage

Considering the time-coupling characteristics of a local ESS, a local EMS adopts a rolling optimization method to determine the scheduling plan of the ESS for next period, which can minimize the cost for the MG to purchase energy from the EMO and determine the role of the MG to participate in energy trading in the next period. The inputs of the rolling optimization are the long-term forecast information of local PV output, local energy demand and the price obtained from the EMO.

2.2.2. Second Stage

According to the short term forecast information collected from MGs and the prices of power grids, the EMO establishes the internal price optimization model for next period based on Stackelberg game theory. The optimization target of the EMO is to maximize its own profits while considering the demand responses of MGs. Furthermore, the EMO broadcasts the results of the internal prices to the MGs, which respond to the prices by adjusting load demands.

3. System Model

3.1. Utility Model of MGs

Each MG equipped with a PV system receives a government subsidy for clean energy plus the revenue from selling surplus energy to MMGs. The MG prefers the PV generation to satisfy its own load demands and charge its own ESS. If the PV generation is insufficient, the MG will buy energy form the MMG. The energy consumption in the MG can create revenue, especially for industrial and commercial users [18,19]. Therefore, the utility model of the MG mainly considers the government subsidy for PV generation, the benefits of energy consumption, the benefits from selling surplus energy to the EMO and the costs of purchasing energy from the EMO. As the MG may be either a buyer or a seller in different periods, the utility function in time slot h can be expressed as follows:
U i h = θ p v i h + k i h ln ( + 1 l i h ) p c b h n l s , i h p c s h n l b , i h c i e i h 2
where θ is the subsidy for each kWh generated by the PV system, p v i h is the PV generation in the time slot h . k i h ln ( 1 + l i h ) is the benefit that the MG i consumes energy l i h . p c b h and p c s h are the buying price and selling price from the MG to the MMG, respectively. n l s , i h and n l b , i h respectively represent the selling energy and buying energy of the MG i in time slot h . c i is the degradation cost coefficient of the ESS in the MG i , and e i h is the charging energy to the ESS in time slot h . The other maintenance and scheduling costs are neglected in this paper. Since each MG in the MMG can only act as a buyer or seller at a special time slot, and the net energy should satisfy the following constraints:
n l b , i h n l s , i h = l i h + e i h p v i h ,
0 n l b , i h D b , i h n l i max ,
0 n l s , i h D s , i h n l i max ,
D b , i h + D s , i h 1 .
Equation (3) is the constraint of energy balance in the MG i . n l i max is the maximum energy transaction of the MG with the power grid in the connection line, D b , i h and D s , i h are the binary variables indicating the state of buying and selling energy of the MG i .

3.2. Profit Model of EMOs

The MG may be either a seller or a buyer depending on its requirement of net energy. The MMG energy importing from the MGs or exporting to the MGs in time slot h can be expressed by
E i m h = D s , i h = 1 n l s , i h ,
E e x h = D b , i h = 1 n l b , i h .
Usually, the mismatch between imported energy and exported energy always exists, and the EMO should trade with the power grid to maintain an internal energy balance. Therefore, the profit function of the EMO can be written as
P r o E M O h = { p c s h E e x h p c b h E i m h + p g b h ( E i m h E e x h ) , E i m h < E e x h p c s h E e x h p c b h E i m h + p g s h ( E i m h E e x h ) , E i m h > E e x h .

4. Optimization Scheduling

4.1. Rolling Optimization for Local MGs

The scheduling of the ESS is optimized based on the long term forecast information on PV generation, load demands and internal prices, which can be expressed as p v i h , f , l i h , f , p c b h , f and p c s h , f , respectively. The primary target of a local MG is to maximize the operation utility. According to (2), the target of the rolling optimization is equivalent to minimize the cost of energy trading with the EMO, which can be expressed as
min C i r o l l = h h + K Δ h [ p c b h , f n l s , i h + p c s h , f n l b , i h + c i e i h 2 ] ,
where k is the length of rolling optimization, and Δ h is the rolling step. In addition, the objective function in (10) needs to satisfy not only the constraints in (3)–(6) but also the constraints of battery charging and discharging, which can be expressed as follows.
e i h = e c h , i h e d i s , i h ,
0 e c h , i h D c h , i h e c h , i max ,
0 e d i s , i h D d i s , i h e d i s , i max ,
D c h , i h + D d i s , i h 1 ,
S O C i h + 1 = S O C i h + D c h , i h e c h , i h η c h , i D d i s , i h e d i s , i h / η d i s , i ,
S O C i min S O C i h S O C i max ,
where e c h , i h and e d i s , i h are the charging and discharging energy of the ESS, respectively. D c h , i h and D d i s , i h indicate the charging and discharging state of the ESS, e c h , i max and e d i s , i max are the maximum charging and discharging energy in time slot h , respectively. Equation (15) shows the state of charge (SOC) of the ESS during its charging and discharging at the end of the time slot h , η c h , i and η d i s , i are the charging and discharging efficiency of the ESS and S O C i min and S O C i max are the lower and upper limits of the SOC. Therefore, the rolling optimal scheduling in (10) can be modeled as a mixed integer programming (MIP) problem.

4.2. Stackelberg Game for EMOs

The optimal scheduling of the ESS in the first stage is regarded as one of the inputs of the EMO in second stage. The other inputs include the short term forecast information on PV generation, schedulable and unscheduled loads and the trading roles of the MG in the following hour. After receiving these inputs, the EMO will determine the internal prices and send to local EMSs for energy consumption optimization in the following hour.

4.2.1. Formulation of a Stackelberg Game

The hour-ahead energy sharing within the MMG is formulated as a Stackelberg game. The EMO is the leader of the game, and stimulates energy sharing by setting the internal prices with the goal of maximizing profits, while the MGs act as followers that optimize their utility through properly responding to internal prices. The game between the EMO and the MGs can be defined by its strategic form as
G = { ( M G O { E M O } ) , { L i h } i N , { P c b h } , { P c s h } , { U i h } i N , P r o M G C h } ,
where { L i h } i N is the set of load strategies adopted by each MG i in the time slot h constrained by l i h , min l i h l i h , max ; { P c b h } and { P c s h } are the strategic set of the EMO, which ensures that the internal prices are constrained by p g s h p c s h < p c b h p g b h ; { U i h } i N and P r o E M O h are the utility of the MG and the profit of the EMO which are expressed by (2) and (9), respectively.
Definition: Consider the game G defined in (17) as a set of strategies ( L i h * , P c b h * , P c s h * ) constituting a Stackelberg equilibrium (SE) if (and only if) the following set of inequalities are satisfied:
U i h ( L i h * , P c b h * , P c s h * ) U i h ( l i h , L i h * , P c b h * , P c s h * )   i N , l i h L i h ,
P r o E M O h ( L i h * , P c b h * , P c s h * ) P r o E M O h ( L i h * , p c b h , p c s h * )   p c b h P c b h ,
P r o E M O h ( L i h * , P c b h * , P c s h * ) P r o E M O h ( L i h * , p c b h * , p c s h )   p c s h P c s h ,
where L i h * = [ l 1 h * , ... , l i h * , ... , l N h * ] , L i h * = [ l 1 h * , ... , l i - 1 h * , l i + 1 h * , ... , l N h * ] . When all the players in ( M G O { E M O } ) reach the SE, the EMO cannot improve its profit by adjusting the internal prices from the SE prices P c b h * and P c s h * . Likewise, no MGs can enhance their utilities by selecting different strategies from L i h * .

4.2.2. Achievement of Game Equilibrium

It is known from (2) and (3) that the utility of MGs can be modified as
U i h = { θ p v i h + k i h ln ( 1 + l i h ) p c b h ( l i h + e i h p v i h ) ,   l i h + e i h p v i h < 0 θ p v i h + k i h ln ( 1 + l i h ) p c s h ( l i h + e i h p v i h ) ,   l i h + e i h p v i h 0 .
For the given internal prices p c b h and p c s h , the optimal energy consumption l i h o can be easily obtained by making U i h / l i h = 0 , which leads to
l i h o = { k i h / p c b h 1 ,   l i h + e i h p v i h < 0 k i h / p c s h 1 ,   l i h + e i h p v i h > 0 .
Equation (22) can be substituted into (7)–(9) to solve the optimal internal price of the EMO based on the role of energy sharing that each MG wants to play. The following section provides a discussion of the relationship between the optimal prices and the roles of MGs.
For each MG acting as a buyer, the net energy n l i h satisfies l i h + e i h p v i h 0 . The range of energy consumption can be redefined as
max ( l i h , min , p v i h e i h ) l i h l i h , max .
Substituting the optimal value in (22) into the energy consumption in (23), we can obtain
k i h l i max + 1 p c s , i h k i h max ( l i h , min , p v i h e i h ) + 1 ,
where p c s , i h is the flexible selling price that the MG i as a buyer expects to trade with the EMO. If p c s h < k i h / ( l i max + 1 ) , the optimal energy consumption will be l i max ; if p c s , i h > k i h / ( max ( l i h , min , p v i h e i h ) + 1 ) , the optimal energy consumption will be max ( l i h , min , p v i h e i h ) . Similarly, the flexible buying price p c b , i h of a seller is constrained by
k i h min ( p v i h e i h , l i h , max ) + 1 p c b , i h k i h l i min + 1 .
If p c b h < k i h / ( min ( p v i h e i h , l i h , max ) + 1 ) , the optimal energy consumption will be min ( p v i h e i h , l i h , max ) ; if p c s h > k i h / ( l i h , min + 1 ) , the optimal energy consumption will be max ( l i h , min , p v i h e i h ) . As a result, we can simplify (24) and (25) into
p c s , i h [ p c s , i h , min , p c s , i h , max ] ,
p c b , i h [ p c b , i h , min , p c b , i h , max ] .
Equation (26) expresses the feasible region of the optimal price for a seller and equation (27) expresses the feasible region of the optimal price for a buyer. Therefore, the profit functions of the EMO in (7)–(9) can be updated with the optimal energy consumption as
E e x h o = p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) D b , i h + p c s h < p s , i h , min l i h , max D b , i h + p c s h > p s , i h , max max ( l i h , min , p v i h e i h ) D b , i h + D b , i h = 1 e i h p v i h
E i m h o = p b , i h , min p c b h p b , i h , max ( k i h p c b h 1 ) D s , i h + p c s h > p s , i h , max l i h , min D s , i h + p c b h < p b , i h , min min ( p v i h e i h , l i h , max ) D s , i h + D s , i h = 1 e i h p v i h
P r o M G C h = { p c s h E e x h o p c b h E i m h o + p g b h ( E i m h o E e x h o ) , E i m h o < E e x h o p c s h E e x h o p c b h E i m h o + p g s h ( E i m h o E e x h o ) , E i m h o E e x h o .
In order to optimize the objective function in (30), we define a two-dimensional coordinate system for the selling price and buying prices, as shown in Figure 3. The critical values in (26) of the sellers and the power grid prices are located in a transverse axis, such as [ p g b h , p c s , i h , min , p c s , j h , min , ... , p c s , k h , max , p c s , l h , max , p g s h ] , and the critical values in (27) of the buyers and the power grid prices are located in a longitudinal axis, such as [ p g b h , p c b , q h , min , p c b , p h , min , ... , p c s , m h , max , p c s , n h , max , p g s h ] . The feasible region S of the profit function (30) is divided into a certain number of sub-regions, such as sub-region s S . For each sub-region s S , it is easily known that the profit function in (30) is concave (the proof can be seen in Appendix A). The optimal buying and selling prices in sub-region s can be found as [ p c b , s h * , p c s , s h * ] by using CPLEX which is the commercial solver for MIP problem. Thus, the optimal prices in S can be obtained as
[ p c b h * , p c s h * ] = arg max ( p c b h , p c s h ) ( P r o M G C h ( p c b , s h * , p c s , s h * ) , s S ) .
Furthermore, the optimal prices [ p c b , s h * , p c s , s h * ] will be sent to MGOs. The MGOs can calculate the optimal load demands L i h * according to formulas (22)–(27). Thus, the SE of the proposed game is reached. The sensitivity analysis of MGs’ utilities over the internal prices adopted by the market operator are shown in Appendix B.

5. Case Studies

5.1. Basic Data

We employed MATLAB software to program the proposed model and analyze the simulation results. The MIP problem and convex optimization problem were solved by CPLEX. The model was applied to an MMG consisting of 3 MGs. All of the MGs had a PV system and an ESS installed, and the maximum schedulable loads were set to nearly 20% of the maximum load demand. The capacity of the ESS in each MG was 100kWh at the maximum charging/discharging rate of 0.5 C, and the range of SOC is from 0.2 to 1. The degradation cost coefficient c i is 0.005.
The PV generation and load demands of the MGs in a typical day are shown in Figure 4, which were collected from the operation data of different MGs located in Jinzhai, Anhui Province, China. The time-of-use tariffs of the distribution network are shown in Table 1. The subsidy for PV energy was CNY 0.42 per kWh. The length of the rolling optimization was 8 hours at a step of 1 hour. The forecast information, which included load demands and PV generation, was generated from the collected historical data by using the prediction algorithms [20,21].

5.2. Internal Prices of the MMG

For the convenience of the following analysis, we provide the optimal prices in the MMG under the condition of market equilibrium as shown in Figure 5. It can be seen that that during the periods of 0:00–8:00 and 18:00–24:00 the internal buying and selling prices were equal to those of the power grid because the PV generation of each MG was very small, or even 0, and the internal load demand was relatively large. While all MGs are in the purchasing state, the EMO cannot improve its own benefits by adjusting the internal prices. During the period of 8:00–18:00, the internal buying price was always higher than that of the power grid, and the MG could sell more electricity by adjusting load demand to maximize operation utility. During the periods of 11:00–12:00 and 13:00–16:00, the internal sale price was lower than the selling price of the power grid, and the MG could increase energy consumption to reduce the purchasing cost of the MG and create more economic benefits.
The above results show that the EMO could only adjust the internal electricity price during the period of 8:00–18:00, when PV generation is relatively strong. During the daytime, it could promote energy trading in the MMG, and realize the improvement of the operation efficiency of each MG. Next, we analyze the benefits of the MGs and the EMO in combination with the internal prices proposed in this paper.

5.3. Results of Local MGs

5.3.1. Rolling Optimization of ESSs

In the proposed method, the local ESS is charged in the case of a PV generation surplus or low electricity price, and discharged in the case of a PV energy deficit or high electricity price to reduce the cost at which the MG purchases the electricity from the power grid. Figure 6 shows the dispatching results of the ESSs in the local MGs. It can be seen that the ESS charged when the internal selling price was equal to the valley price of the power grid in the period of 0:00–6:00. In the period of 7:00–9:00, as load demands increased, the local PV energy was insufficient, the internal selling price was high and discharge of the ESS occurred to reduce the cost of buying energy from the power grids in that period. During the period of 15:00–18:00, the PV generation surplus was absorbed by the ESS. During the period of 18:00–22:00, when the electricity price was high, the absorbed energy was released to reduce the cost of power consumption in the MG.
In order to further analyze the benefits of ESS scheduling, the benefits of ESS charging and discharging in each period are shown in Figure 7. These benefits were calculated on the basis of the internal prices. The results show that the positive benefits of each ESS in a day were higher than the total negative benefits. The increased benefits of three MGs were CNY 130.08, CNY 133.66 and CNY 143.96, respectively. This validates the necessity and rationality of ESS scheduling in the first stage.

5.3.2. Demand Response of Local MGs

On the basis of optimal ESS scheduling, each MG will determine its own trading role (see Appendix C) to participate in energy trading in the MMGs by transferring the local information to the EMO. Under the incentive of internal prices, the MGs in the state of selling are encouraged to sell more energy to the EMO, and the MGs in the state of buying are encouraged to buy more energy from the EMO through the load adjustment as shown in Figure 8. During the periods of 11:00–12:00 and 13:00–16:00, MG1 and MG3 were stimulated by internal selling prices in the corresponding hours, which improved the operation utility by increasing energy consumption. On the contrary, the MG2 in the state of selling during the period of 7:00–18:00 could improve benefits by reducing load demand and selling more energy. Compared with the direct transaction with the power grid, the proposed method improves the utility of MGs, and the increased utility in each hour is shown in Figure 9.

5.4. Results of the EMO

The adjustment of energy consumption in local MGs directly affects the total energy sharing and total energy demand in the energy market, as shown in Figure 10. It can be seen that as the internal selling prices were lower than the grid selling prices during the periods of 11:00–12:00 and 13:00–16:00, the total energy demand in the energy market increased significantly. The total energy sharing in the energy market also increased as internal buying prices rose higher than the grid buying prices during the period of 7:00–18:00. The results show that MGs will actively participate in energy market transactions under the incentive of internal electricity prices.
Energy demand and sharing in the energy market is the premise for the EMO to make profits. In the periods of 0:00–7:00 and 18:00–24:00, MGs were in the buying state thus the EMO could only purchase energy from the power grids to meet the demands of MGs and could not obtain profits. Figure 11 compares the profits of the EMO under the internal price strategy and the grid price strategy during the period of 7:00 and 18:00. Under the incentive of the internal price strategy, the profit of the EMO increased significantly in these periods. The daily profit of the EMO increased from CNY 171.6 to CNY 277.7 and the rate of increase is 61.82%.
To further illustrate the advantages of the proposed method, the net energy and PV utilization of the MMG were analyzed. The comparison of net energy curves is shown in Figure 12. During the period of 0:00–6:00, the charging of the ESS led to an increase in net energy. During the period of 18:00–24:00, the discharging of the ESS led to a decline in net energy. During the period of 6:00–18:00, the change in net energy was affected by ESS scheduling and the load demand response under internal prices. Obviously, the proposed method produced smaller fluctuations of the net load than the original method, and the PV energy reversal in the MMG during the period of 10:00–15:00 was significantly suppressed.
Moreover, the comparison of peak-to-average ratio (PAR) and PV utilization ratio are shown in Table 2. Through comparisons, it is found that the proposed method can effectively reduce the peak-to-average ratio of the net load and improve the utilization ratio of PV energy.

5.5. Utility Comparisons with Other Methods

To further illustrate the advantages of the proposed hour-ahead optimization over day-ahead optimization, the methods in references [16] and [17] were applied to the proposed case, and the cost optimization problem of references [16] and [17] were transformed into the utility optimization problem to make a comparative analysis of the results. In addition, the day-ahead stochastic prediction errors of PV and the load demand were set to 10% and 12%, while the hour-ahead stochastic prediction errors were set to 5% and 6%, respectively. The results are shown in Table 3.
The results show that EMO and MGO can achieve more benefits by using the optimization method proposed in this paper.

6. Conclusions

In this paper, a hierarchical optimization method is proposed to optimize the energy transaction of a MMG in two stages. Firstly, the EMS of each MG determines the scheduling of ESS in the next hour by adopting rolling optimization, and decides its participating role in energy trading market. Secondly, according to forecast information and the energy trading roles collected from MGs, the EMO optimizes the internal prices of next hour based on Stackelberg game theory. The simulation results show that the utility of both the EMO and MGO are increased by using the proposed method. In addition, the net load curve and utilization ratio of PV energy in the whole MMG system are both improved.

Author Contributions

Formal analysis, C.H., W.S. and B.X.; Funding acquisition, Q.W.; Methodology, T.R. and G.L.; Writing—original draft, T.R.

Funding

This work was supported by the National Key R&D Program of China (Nos. 2016YFB0900400)

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

It is assumed that
C e x = p c s h < p s , i h , min l i h , max D b , i h + p c s h > p s , i h , max max ( l i h , min , p v i h e i h ) D b , i h + D b , i h = 1 e i h p v i h ,
C i m = p c s h > p s , i h , max l i h , min D s , i h + p c b h < p b , i h , min min ( p v i h e i h , l i h , max ) D s , i h + D s , i h = 1 e i h p v i h .
Thus, expressions (32) and (33) can be simplified as
E e x h o = p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) + C e x ,
E i m h o = p b , i h , min p c s h p b , i h , max ( k i h p c b h 1 ) + C i m .
By substituting (42) and (43) into (34), the objective function of the EMO in sub-region s S is described as:
If E i m h o < E e x h o
P r o E M O h = p c s h ( p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) D b , i h + C e x ) p c b h ( p b , i h , min p c b h p b , i h , max ( k i h p c b h 1 ) D s , i h + C i m ) + p g b h ( p b , i h , min p c b h p b , i h , max ( k i h p c b h 1 ) D s , i h + C i m ) p g b h ( p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) D b , i h + C e x )
If E i m h o E e x h o
P r o E M O h = p c s h ( p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) D b , i h + C e x ) p c b h ( p b , i h , min p c s h p b , i h , max ( k i h p c b h 1 ) D s , i h + C i m ) + p g s h ( p b , i h , min p c s h p b , i h , max ( k i h p c b h 1 ) D s , i h + C i m ) p g s h ( p s , i h , min p c s h p s , i h , max ( k i h p c s h 1 ) D b , i h + C e x )
where ( p c b h , p c s h ) s , and p c b h < p c s h . Therefore, the Hessian matrix of P r o E M O h is
H = { [ 2 p g b h p c s h 3 p s , i h , min p c s h p s , i h , max D b , i h k i h 0 0 2 p g b h p c b h 3 p b , i h , min p c b h p b , i h , max D s , i h k i h ] , E i m h o < E e x h o [ 2 p g s h p c s h 3 p s , i h , min p c s h p s , i h , max D b , i h k i h 0 0 2 p g s h p c b h 3 p b , i h , min p c b h p b , i h , max D s , i h k i h ] , E i m h o E e x h o .
Consider the fact that k i h > 0 , p g b h > 0 , p g s h > 0 , p c s h > 0 , p c b h > 0 i N , h H , H is thus negative definite and P r o M G C h is strictly concave with respect to p c s h and p c b h .

Appendix B

Formulas (26) and (27) show that if p c s , i h [ p c s , i h , min , p c s , i h , max ] , then U i h / p c s , i h = 0 , and if p c b , i h [ p c b , i h , min , p c b , i h , max ] , then U i h / p c b , i h = 0 . Therefore, we mainly discuss the case of p c s , i h [ p c s , i h , min , p c s , i h , max ] and p c b , i h [ p c b , i h , min , p c b , i h , max ] . For MG i , who is a buyer, if p c s , i h [ p c s , i h , min , p c s , i h , max ] , then the optimal load demand can be expressed as
l i h o = k i h / p c b h 1 = ( l i h f + 1 ) p g b h / p c b h 1 .
The corresponding utility is rewritten as
U i h = θ p v i h + ( l i h f + 1 ) p g b h ln ( ( l i h f + 1 ) p g b h p c b h ) p c b h ( ( l i h f + 1 ) p g b h p c b h 1 + e i h p v i h ) .
Therefore, the derivative of benefit U i h with respect to price p c b h is
U i h p c b , i h = ( l i h f + 1 ) p g b h p c b h 1 + e i h p v i h .
For MG i , who is a seller, if p c b , i h [ p c b , i h , min , p c b , i h , max ] , then the optimal load demand can be expressed as
l i h o = k i h / p c s h 1 = ( l i h f + 1 ) p g b h / p c s h 1 .
The corresponding utility is rewritten as
U i h = θ p v i h + ( l i h f + 1 ) p g b h ln ( ( l i h f + 1 ) p g b h p c s h ) p c s h ( ( l i h f + 1 ) p g b h p c s h 1 + e i h p v i h ) .
Therefore, the derivative of benefit U i h with respect to price p c s h is
U i h p c s , i h = ( l i h f + 1 ) p g b h p c s h 1 + e i h p v i h .
To further demonstrate the utility of MG response to the internal electricity price, the utilities of MG1 and MG3 in time slot 14 were chosen to show the relationship between utility and changing prices, which are shown in Figure A1 and Figure A2, respectively.
Figure A1. The utility of MG1 response to internal buying prices. ( k 1 14 = 46.1 , l 1 14 f = 130.7 , p v 1 14 = 141.6 , e 1 14 = 0 , p c s 14 = 0.95 , p g b 14 = 0.35 , p g s 14 = 1.189 ).
Figure A1. The utility of MG1 response to internal buying prices. ( k 1 14 = 46.1 , l 1 14 f = 130.7 , p v 1 14 = 141.6 , e 1 14 = 0 , p c s 14 = 0.95 , p g b 14 = 0.35 , p g s 14 = 1.189 ).
Applsci 09 00624 g0a1
Figure A2. The utility of MG3 response to internal selling prices. ( k 3 14 = 271.6 , l 3 14 f = 227.4 , p v 3 14 = 172.7 , e 3 14 = 30 , p c b 14 = 0.49 , p g b 14 = 0.35 , p g s 14 = 1.189 ).
Figure A2. The utility of MG3 response to internal selling prices. ( k 3 14 = 271.6 , l 3 14 f = 227.4 , p v 3 14 = 172.7 , e 3 14 = 30 , p c b 14 = 0.49 , p g b 14 = 0.35 , p g s 14 = 1.189 ).
Applsci 09 00624 g0a2

Appendix C

The optimal results of D b , i h and D b , i h determine the roles of MGs participating in energy trading, which are all listed in the follow Table A1.
Table A1. Roles of Energy Trading.
Table A1. Roles of Energy Trading.
PeriodMG1MG2>MG3
0:00–8:00buyerbuyerbuyer
8:00–9:00buyerbuyerseller
9:00–10:00buyerbuyerseller
10:00–11:00buyerbuyerseller
11:00–12:00buyerbuyerseller
12:00–13:00buyersellerbuyer
13:00–14:00sellersellerbuyer
14:00–15:00sellersellerbuyer
15:00–16:00sellersellerbuyer
16:00–17:00sellerbuyerbuyer
17:00–18:00sellersellerbuyer
18:00–24:00buyerbuyerbuyer

References

  1. Rafique, S.F.; Zhang, J. Energy management system, generation and demand predictors: A review. IET Gener. Transm. Distrib. 2018, 7, 519–530. [Google Scholar] [CrossRef]
  2. Han, Y.; Zhang, K.; Li, H.; Coelho, E.A.; Guerrero, J.M. MAS-based Distributed Coordinated Control and Optimization in Microgrid and Microgrid Clusters: A Comprehensive Overview. IEEE Trans. Power Electron. 2018, 33, 6488–6508. [Google Scholar] [CrossRef]
  3. Nguyen, T.A.; Crow, M.L. Stochastic Optimization of Renewable-Based Microgrid Operation Incorporating Battery Operating Cost. IEEE Trans. Power Syst. 2016, 31, 2289–2296. [Google Scholar] [CrossRef]
  4. Liu, W.; Gu, W.; Sheng, W.; Meng, X.; Wu, Z.; Chen, W. Decentralized Multi-Agent System-Based Cooperative Frequency Control for Autonomous Microgrids With Communication Constraints. IEEE Trans. Sustain. Energy 2014, 5, 446–456. [Google Scholar] [CrossRef]
  5. Li, Q.; Chen, F.; Chen, M.; Guerrero, J.M.; Abbott, D. Agent-Based Decentralized Control Method for Islanded Microgrids. IEEE Trans. Smart Grid 2016, 7, 637–649. [Google Scholar] [CrossRef]
  6. Parisio, A.; Wiezorek, C.; Kyntäjä, T.; Elo, J.; Strunz, K.; Johansson, K.H. Cooperative MPC-Based Energy Management for Networked Microgrids. IEEE Trans. Smart Grid 2017, 8, 3066–3074. [Google Scholar] [CrossRef]
  7. Najmeh, B.; Ahmadreza, T.; Amjad, A.; Josep, M.G. Optimal operation management of a regionalnetwork of microgrids based on chanceconstrained model predictive control. IET Gener. Transm. Distrib. 2018, 12, 3772–3779. [Google Scholar]
  8. Rafiee Sandgani, M.; Sirouspour, S. Priority-based Microgrid Energy Management in a Network Environment. IEEE Trans. Sustain. Energy 2018, 9, 980–990. [Google Scholar] [CrossRef]
  9. Jadhav, A.M.; Patne, N.R.; Guerrero, J.M. A Novel Approach to Neighborhood Fair Energy Trading in a Distribution Network of Multiple Microgrid Clusters. IEEE Trans. Ind. Electron. 2018, 66, 1520–1531. [Google Scholar] [CrossRef]
  10. Zhang, B.; Li, Q.; Wang, L.; Feng, W. Robust optimization for energy transactions in multi-microgrids under uncertainty. Appl. Energy 2018, 217, 346–360. [Google Scholar] [CrossRef]
  11. Bui, V.H.; Hussain, A.; Kim, H.M. A Multiagent-Based Hierarchical Energy Management Strategy for Multi-Microgrids Considering Adjustable Power and Demand Response. IEEE Trans. Smart Grid 2018, 9, 1323–1333. [Google Scholar] [CrossRef]
  12. Maryam, M.; Hassan, M.; Amjad, A.; Josep, M.G.; Hamid, L. A decentralized robust model for optimal operation of distribution companies with private microgrids. Electr. Power Energy Syst. 2019, 106, 105–123. [Google Scholar]
  13. Manshadi, S.D.; Khodayar, M.E. A Hierarchical Electricity Market Structure for the Smart Grid Paradigm. IEEE Trans. Smart Grid 2016, 7, 1866–1875. [Google Scholar] [CrossRef]
  14. Yue, J.; Hu, Z.; Amjad, A.M.; Josep, M.G. A Multi-Market-Driven Approach to Energy Scheduling of Smart Microgrids in Distribution Networks. Sustainability 2019, 11, 301. [Google Scholar] [CrossRef]
  15. Fan, S.; Ai, Q.; Piao, L. Bargaining-based cooperative energy trading for distribution company and demand response. Appl. Energy 2018, 226, 469–482. [Google Scholar] [CrossRef]
  16. Liu, Y.; Guo, L.; Wang, C. A robust operation-based scheduling optimization for smart distribution networks with multi-microgrids. Appl. Energy 2018, 228, 130–140. [Google Scholar] [CrossRef]
  17. Jalali, M.; Zare, K.; Seyedi, H. Strategic decision-making of distribution network operator with multi-microgrids considering demand response program. Energy 2017, 141, 1059–1071. [Google Scholar] [CrossRef]
  18. Liu, N.; Cheng, M.; Yu, X.; Zhong, J.; Lei, J. Energy Sharing Provider for PV Prosumer Clusters: A Hybrid Approach using Stochastic Programming and Stackelberg Game. IEEE Trans. Ind. Electron. 2018, 65, 6740–6750. [Google Scholar] [CrossRef]
  19. Liu, N.; Yu, X.; Wang, C.; Wang, J. Energy Sharing Management for Microgrids with PV Prosumers: A Stackelberg Game Approach. IEEE Trans. Ind. Inform. 2017, 13, 1088–1098. [Google Scholar] [CrossRef]
  20. Abedinia, O.; Amjady, N.; Zareipour, H. A New Feature Selection Technique for Load and Price Forecast of Electrical Power Systems. IEEE Trans. Power Syst. 2017, 32, 62–74. [Google Scholar] [CrossRef]
  21. Gigoni, L.; Betti, A.; Crisostomi, E. Day-Ahead Hourly Forecasting of Power Generation from Photovoltaic Plants. IEEE Trans. Sustain. Energy 2018, 9, 831–842. [Google Scholar] [CrossRef]
Figure 1. Structure of an MMG (multiple microgrid).
Figure 1. Structure of an MMG (multiple microgrid).
Applsci 09 00624 g001
Figure 2. Operation strategy of an MMG (EMO: Energy Market Operator, MGO: Microgrid Operator, PV: photovoltaic).
Figure 2. Operation strategy of an MMG (EMO: Energy Market Operator, MGO: Microgrid Operator, PV: photovoltaic).
Applsci 09 00624 g002
Figure 3. Feasible region of internal prices.
Figure 3. Feasible region of internal prices.
Applsci 09 00624 g003
Figure 4. Basic data in hour: (a) PV energy output for 3 MGs, and (b) load demands for 3 MGs.
Figure 4. Basic data in hour: (a) PV energy output for 3 MGs, and (b) load demands for 3 MGs.
Applsci 09 00624 g004
Figure 5. Comparison of internal prices and grid prices.
Figure 5. Comparison of internal prices and grid prices.
Applsci 09 00624 g005
Figure 6. Optimal scheduling of energy storage systems (ESSs). SOC: state of charge.
Figure 6. Optimal scheduling of energy storage systems (ESSs). SOC: state of charge.
Applsci 09 00624 g006
Figure 7. Impact of ESS scheduling on the benefits of MGs.
Figure 7. Impact of ESS scheduling on the benefits of MGs.
Applsci 09 00624 g007
Figure 8. Adjustment of load demands.
Figure 8. Adjustment of load demands.
Applsci 09 00624 g008
Figure 9. Increased utility of MGs based on internal prices.
Figure 9. Increased utility of MGs based on internal prices.
Applsci 09 00624 g009
Figure 10. Comparison of energy trading in the energy market.
Figure 10. Comparison of energy trading in the energy market.
Applsci 09 00624 g010
Figure 11. Comparison of the utility of the EMO.
Figure 11. Comparison of the utility of the EMO.
Applsci 09 00624 g011
Figure 12. Comparison of net energy of MMGs.
Figure 12. Comparison of net energy of MMGs.
Applsci 09 00624 g012
Table 1. Time-of-use tariff.
Table 1. Time-of-use tariff.
Distribution NetworkPrices (kWh/h)Hours
SellingPeak: 1.1898:00–11:00; 13:00–16:00; 18:00–22:00
Flat: 0.7387:00–8:00; 11:00–13:00; 16:00–18:00
Valley: 0.4230:00–7:00; 22:00–24:00
Buying0.3520:00–24:00
Table 2. Comparison of pea-to-average ratio (PAR) and PV utilization ratio.
Table 2. Comparison of pea-to-average ratio (PAR) and PV utilization ratio.
MethodPARPV Utilization Ratio
Original method3.159685.25%
Proposed method2.699298.07%
Table 3. Comparisons with other methods.
Table 3. Comparisons with other methods.
Utility/CNYReference [16]Reference [17]Proposed
EMO234.2207.6277.8
MG14792.34796.44823.0
MG25897.25952.05973.8
MG36897.86946.47011.6

Share and Cite

MDPI and ACS Style

Rui, T.; Li, G.; Wang, Q.; Hu, C.; Shen, W.; Xu, B. Hierarchical Optimization Method for Energy Scheduling of Multiple Microgrids. Appl. Sci. 2019, 9, 624. https://doi.org/10.3390/app9040624

AMA Style

Rui T, Li G, Wang Q, Hu C, Shen W, Xu B. Hierarchical Optimization Method for Energy Scheduling of Multiple Microgrids. Applied Sciences. 2019; 9(4):624. https://doi.org/10.3390/app9040624

Chicago/Turabian Style

Rui, Tao, Guoli Li, Qunjing Wang, Cungang Hu, Weixiang Shen, and Bin Xu. 2019. "Hierarchical Optimization Method for Energy Scheduling of Multiple Microgrids" Applied Sciences 9, no. 4: 624. https://doi.org/10.3390/app9040624

APA Style

Rui, T., Li, G., Wang, Q., Hu, C., Shen, W., & Xu, B. (2019). Hierarchical Optimization Method for Energy Scheduling of Multiple Microgrids. Applied Sciences, 9(4), 624. https://doi.org/10.3390/app9040624

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