Next Article in Journal
A Simple Method for Assessing Diversity and Dynamics of Microbial Community: Comparison of Dairy Phages from Industrial and Spontaneous Fermentation
Next Article in Special Issue
Performance Analysis of a Stand-Alone PV/WT/Biomass/Bat System in Alrashda Village in Egypt
Previous Article in Journal
Agricultural Waste-Based Biochar for Agronomic Applications
Previous Article in Special Issue
Optimal Incorporation of Photovoltaic Energy and Battery Energy Storage Systems in Distribution Networks Considering Uncertainties of Demand and Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cooperative Multi-Objective Optimization of DC Multi-Microgrid Systems in Distribution Networks

State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(19), 8916; https://doi.org/10.3390/app11198916
Submission received: 7 August 2021 / Revised: 14 September 2021 / Accepted: 21 September 2021 / Published: 24 September 2021
(This article belongs to the Special Issue Renewable-Based Microgrids: Design, Control and Optimization)

Abstract

:

Featured Application

A cooperative multi-objective optimization model of a DC multi-microgrid that considers across-time-and-space energy transmission of EVs is established to improve the economy of the system, decrease the loss of the distribution network, and reduce carbon emissions.

Abstract

By constructing a DC multi-microgrid system (MMGS) including renewable energy sources (RESs) and electric vehicles (EVs) to coordinate with the distribution network, the utilization rate of RESs can be effectively improved and carbon emissions can be reduced. To improve the economy of MMGS and reduce the network loss of the distribution network, a cooperative double-loop optimization strategy is proposed. The inner-loop economic dispatching reduces the daily operating cost of MMGS by optimizing the active power output of RESs, EVs, and DC/AC converters in MMGS. The outer-loop reactive power optimization reduces the network loss of the distribution network by optimizing the reactive power of the bidirectional DC/AC converters. The double-loop, which synergistically optimizes the economic cost and carbon emissions of MMGS, not only improves the economy of MMGS and operational effectiveness of the distribution network but also realizes the low-carbon emissions. The Across-time-and-space energy transmission (ATSET) of the EVs is considered, whose impact on economic dispatching is analyzed. Particle Swarm Optimization (PSO) is applied to iterative solutions. Finally, the rationality and feasibility of the cooperative multi-objective optimization model are proved by a revised IEEE 33-node system.

1. Introduction

Since the national carbon neutrality and carbon peak requirements have been put forward [1], low carbon emissions and new energy have become hot research topics [2,3]. It is a trend to replace petrol vehicles with electric vehicles (EVs) and replace regional large-scale power grids with microgrids (MGs) containing renewable energy sources (RESs) [4,5]. As the output of RESs is intermittent and uncertain, the MGs need to coordinate with the distribution network to centrally regulate the RESs, which is a challenge to the operation mode of the traditional power system. With the popularity of EVs, the burden of the distribution network will greatly be increased. Additionally, the safe operation of the distribution network will be threatened if EVs are charged in the distribution network without control.
The research on the charging and discharging dispatching strategy of EVs is mainly from the view of the economy [6,7]. Many studies have considered charging/discharging strategies of EVs but overlooked the energy storage characteristics of EVs. Through the bidirectional Vehicle-to-grid (V2G) technology, EVs can also deliver electrical energy to the grid by discharging, and improve the operation of the grid [8,9,10].
Some research combines EVs with distributed RESs in the MG. In [11], an optimization method for the operation route and charging/discharging time of EVs is proposed, which uses the timely charging/discharging of EVs to consume the output of RESs and reduce the volatility of the equivalent load. In [12], the MG energy management strategy is discussed from the perspective of system operating cost and the consumption efficiency of RESs, and V2G technology has been applied. In [13], the structure and parameter design of the system have been discussed, and the actual MG system using V2G technology has been studied. However, most of the energy scheduling in MGs and the distribution network is to adjust the output of active power.
Some research focuses on the economic dispatch of the multi-microgrid system. In [14], an interconnected multi-microgrids (IMMGs) system using various complementary power sources effectively coordinates the energy sharing/trading among the MGs and the main grid to improve energy efficiency. In [15], A probabilistic modeling of both small-scale energy resources (SSERs) and load demand at each microgrid (MGs) is performed to determine the optimal economic operation of each MG with minimum cost based on the power transaction between the MGs and the main grid. The above does not consider the reactive power exchange between MMG and the main grid.
In the current research on reactive power exchange and network loss, most studies focus on the reactive power of a single distribution network. In [16], the trend of reactive power demand in the distribution network is evaluated. Reactive power demand management plays an important role in the cost-effectiveness and stable operation of the distribution network. A multi-objective planning algorithm for reactive power compensation of radial distribution networks is proposed in [17], which uses unified power quality conditioner (UPQC) compensation for load reactive power to reduce network loss. In [18], the solid-state transformer (SST) is used to supply the load reactive power demand and inject reactive power into the grid, which reduces network losses in a radial distribution network.
Some research focuses on the impact of reactive power optimization on the loss of MG. In [19], a distributed, leaderless and randomized algorithm is proposed, which controls the microgenerators in the island-operated MG system to compensate for reactive power and reduce power distribution loss in MG. A generalized approach for probabilistic optimal reactive power planning is proposed in [20], which can reduce the annual energy losses of the grid-connected MG system.
These papers mentioned above give less consideration to the collaborative optimization of MG clusters and the distribution network. To solve the above problems, a cooperative multi-objective optimization model of a DC multi-microgrid system (MMGS) including RESs, EVs, and DC/AC converters is established. The goal of the model is to obtain the optimal MMGS economic cost and the network loss of the distribution network. The main contributions of this paper are as follows:
  • A grid-connected MMGS containing RESs and EVs is constructed, where RESs, EVs, MGs and distribution networks are combined, bidirectional V2G technology is used and the across-time-and-space energy transmission (ATSET) of EVs is thoroughly considered. The effect of the across-time-and-space energy transmission on MMGS economic operation is analyzed to state the potential benefits of cooperative multi-objective optimization.
  • A cooperative multi-objective optimization model is established, including the dynamic economic dispatch of RESs, EVs, DC/AC converters, and the reactive power optimization of DC/AC converters in MMGS. The cooperative multi-objective optimization model consists of two loops. The inner-loop model uses the active power output of RESs, EVs, DC/AC converters as variables, and the daily operating cost of MMGS is used as the optimization objective. The outer-loop model uses the reactive power output of the DC/AC converters as the variable to optimize the network loss of the distribution network, thereby reducing network loss cost and carbon emissions cost. The ultimate goal of the cooperative multi-objective is to obtain the optimal daily economic cost.
  • The concepts of carbon neutrality and carbon peaking are combined. Through the cooperative multi-objective optimization model, the carbon emissions generated by the operation of the MMGS and the distribution network are effectively reduced. The cooperative multi-objective optimization model not only improves the economy but also reduces the total carbon emissions of MMGS and the distribution network.

2. System Structure

2.1. Structure of the DC Multi-Microgrid System

The MMGS discussed in this paper includes multiple relatively independent MGs in space. The DC multi-microgrid energy management system (MMGEMS) manages all energy transactions in MMGS. Each MG is integrated into the distribution network through power electronic devices and exchanges energy with the distribution network. Each MG contains RESs and EVs charging/discharging infrastructures (EVCDIs). There are two main types of MGs in the MMGS: MGs located in residential areas (RMG) and MGs located in office areas (OBMG). The structure of the MMGS is shown in Figure 1.
The control of the system is mainly conducted by the collaboration of the MG energy management system (MEMS) and the EVs management system (EVMS). The MEMS is responsible for the energy dispatching of photovoltaics (PVs), wind turbines (WTs), and EVs in MGs, and the EVMS manages the charging and discharging behaviors of EVs. A DC multi-microgrid control system is shown in Figure 2.

2.2. DC Microgrid

The basic structure of the DC microgrid is shown in Figure 3. Each MG is connected to the distribution network through a transformer and a DC/AC converter, which can exchange energy with the distribution network. A connection switch is installed in the grid-connected circuit, which can switch the MG between island operation mode and grid-connected operation mode.

2.3. Bidirectional DC/AC Converter

The bidirectional DC/AC converters are used to connect the MMGS and the distribution network, which can output active and reactive power with the distribution network. DC/AC converters use power factor correction (PFC) to obtain the unity power factor [21]. Therefore, the DC/AC converters are set to the unity power factor in this paper when the reactive power is not optimized [22]. However, by using the appropriate pulse-width modulation (PWM) switching technique, the power factor is adjusted to control the reactive power output of the DC/AC converters to the distribution network [22]. This is the basis for reactive power optimization.

3. Mathematical Model

3.1. Renewable Power Generation

3.1.1. Photovoltaic Module

In this paper, the power prediction module based on artificial neural networks (ANNs) [23] is applied to the economic dispatching of MMGS. The weather data are from the numerical weather forecast (NWP).

3.1.2. Wind Turbine

The output of WTs is mainly affected by wind speed [24]. The ANN is still used to predict wind power [25]. The inputs are the wind speed and wind direction from NWP.

3.2. Electric Vehicles

3.2.1. EVs Model

The EVs in the MMGS are all commuter vehicles, and the residents of the residential area are the workers in the office area. As 77.95% of EVs’ users will reach the working area at 7:30–9:30 [26], the standard parking time slots in OBMG and RMG are assumed to be 9:00–17:00 and 19:00–7:00 [27,28]. The capacity of EV at t-th is
S O C E V m , n , t = S O C E V m , n , t 1 ( 1 σ ) + P m , n , t E V Δ t η C E V , i f P m , n , t E V 0
S O C E V m , n , t = S O C E V m , n , t 1 ( 1 σ ) + P m , n , t E V Δ t / η D E V , i f P m , n , t E V < 0
The power output of EV in (1) and (2) is measured on the MMGS side. Where SOCEVm,n,t is the remaining power capacity of the n-th EV in the m-th MG in the t-th hour, σ is the self-discharge coefficient. P m , n , t E V is the charging or discharging power in the t-th hour of the n-th EV in the m-th MG. If P m , n , t E V ≥ 0, EVs are charged. If P m , n , t E V < 0, EVs release energy; ∆t = 1 h. ηDEV and ηCEV are the discharging and charging efficiency of EVs to calculate the power actual charging or discharging power of EVs.

3.2.2. Across-Time-and-Space Energy Transmission of EV

In the same MG, the EV is used as an energy storage unit, and its charging/discharging power can be dispatched for the operation of the MG. When the EV is connected to the MG and the power is sufficient, MMGS controls the EV to charge during the low charging price or when the system has excess energy, and discharge during the peak discharging price or when the system is short of power. The EV is charged and discharged in the same MG to realize energy transfer over time, thereby reducing the cost of MMGS purchasing electricity directly from the distribution network. At the same time, it also allows the user of the EV to profit by selling part of the electricity, which enables both parties to obtain a certain amount of economic benefit.
On the other hand, EVs not only have energy storage characteristics but also can move between different locations. In the case of differences in the electricity price of the distribution network within a region, benefits can be obtained through the cross-space transfer of energy. For example, the electricity prices of RMG and OBMG for electricity trading with the distribution network are quite different. Most of the time, the electricity price of OBMG purchasing electricity from the distribution network is higher than RMG. Therefore, the electric energy charged in the RMG at a low charging price is sold to the MMGS at a high discharging price in the OBMG, and the electric energy is transferred between different spaces and times through charging and discharging.
EVs realize the across-time energy transmission in the same MG and realize the across-time-and-space energy transmission in different MGs, which can transfer the lower-priced electric energy in RMG to OBMG at a higher price. Under the right conditions, both MMGS and EV users can benefit. This characteristic of EVs for energy transfer between different times and different spaces is called the across-time-and-space energy transmission.

3.3. EV Charging/Discharging Infrastructures

The charging/discharging behaviors of EVs are carried out through the EVCDIs.
P m , n , t E V C D I s = P m , n , t E V η C E V , i f P m , n , t E V 0
P m , n , t E V C D I s = P m , n , t E V / η D E V , i f P m , n , t E V < 0
where P m , n , t E V C D I s is the power of the EVCDIs of the n-th EV in the m-th MG in the t-th hour.

4. Cooperative Multi-Objective Optimization Model

4.1. Description of the Optimization Model

The EVMS collects the dispatchable capacity forecast data of EVs and outputs the dispatching plan of the EVs. The MEMS collects the output of predicted RESs, the predicted load data, and the energy price of the distribution network. Based on this information, MEMS outputs the active power of RESs, EVs, and DC/AC converters in MMGS, and transmits it to the reactive power optimization module in MMGEMS. The reactive power optimization module outputs the optimal reactive power of the DC/AC converters according to the data. The two modules coordinate and output the optimal result.

4.2. Double-Loop Optimization Process

The process is shown in Figure 4. The inner-loop is the dynamic economic dispatch which is used to optimize the active power output of RESs, EVs, and DC/AC converters to obtain the optimal total operating cost of the MMGS. The outer-loop optimizes the reactive power output of the DC/AC converters according to the active power output of the inner-loop, to make the network loss of the distribution network minimum, thereby reducing the network loss cost and carbon emissions of MMGS and the distribution network. The inner-loop and the outer-loop work together to obtain the optimal active power output plan in MMGS and reactive power output of the DC/AC converters, which makes the economic cost of MMGS minimum.

4.3. Cooperative Multi-Objective Optimization Objective Function

The main goal of optimization is to reduce the daily economic total cost of MMGS. MMGS discussed in this model consists of multiple MGs, which are assumed to be owned by a single operator. Another goal of the model is the lowest network loss of the distribution network, which can be obtained through the outer-loop model. Therefore, the objective function to minimize the total economic cost of the entire system can be expressed as:
f = C E T C = C O T C + C W T C
where f is the main goal of the cooperative optimization, CETC is the economic total cost of MMGS. COTC is the operating total cost of MMGS, CWTC is the energy loss cost of the MMGS that is obtained from the outer-loop model, where
C W T C = C i l + C c o
C c o = E C O k c
E C O = ( W S G W S B ) Δ t e c
C i l = ( W S G W S B ) Δ t k i l
E C = m = 1 M E C I m + E C O
Cil and Cco are the network loss cost and carbon emissions cost caused by the increase in the distribution network loss in the outer-loop model, respectively. ECO is the carbon emissions generated by the distribution network. W S G is the total daily operating network loss of the distribution network when MMGS is integrated into the distribution network and runs.   W S B is the total daily operating network loss when there is no MMGS access, which is a fixed value also called the original baseline loss. kil, ec, kc are fixed factors, kil is the loss cost coefficient of the distribution network, ec is the carbon emissions factor, kc is the carbon cost factor. ∆t = 1 h. EC is the total carbon emissions of MMGS and the distribution network, ECIm is the carbon emissions generated by m-th MG in the inner-loop model, M is the number of MGs in the MMGS.
Since COTC and W S G are the optimization targets of the inner-loop model and the outer-loop model, respectively, the objective functions of the inner-loop model and the outer-loop model are set as follows:
f 1 = min C O T C
f 2 = min W S G
where f1 and f2 are the objective functions of the inner-loop model and the outer-loop model, respectively.
Through (5)–(12), f can be expressed as:
f = f 1 + ( f 2 W S B ) Δ t [ e c k c + k i l ]

4.4. Inner-Loop Optimization

The goal of the inner-loop optimization model is to minimize the daily operating cost of the MMGS. The daily operating cost is mainly composed of system energy transaction cost and carbon emissions cost. The objective function is as follows:
f 1 = min C O T C
C O T C = m = 1 M C O C m
C O C m = C e x m + C c i m
C c i m = E C I m k c
COCm is the operating cost of the m-th MG that is obtained from the inner-loop model. Cexm is the energy transaction cost in the m-th MG. Ccim is the carbon emissions cost due to energy exchange in the inner-loop model.

4.4.1. Energy Transaction Cost

The energy transaction cost is the sum of RESs cost, energy exchange cost between MMGS and EVs, MMGS and distribution network, and the additional cycle cost of EV batteries. P m , t E V , P m , t P V , P m , t W T , and P m , t G are the optimization variables.
C e x m = C r e s m C e v m + C g m + C c y m
C r e s m = C P V m + C W T m
C P V m + C W T m = t = 1 T P m , t P V C m , t P V Δ t + t = 1 T P m , t W T C m , t W T Δ t
C e v m = t = 1 T n = 1 N P m , n , t E V η C E V C m , t C E V Δ t , i f P m , n , t E V 0
C e v m = t = 1 T n = 1 N P m , n , t E V C m , t D E V Δ t , i f P m , n , t E V < 0
C g m = t = 1 T P m , t G C m , t G Δ t
C c y m = n = 1 N k c y C c y n E V
where Cresm is the cost of RESs of the m-th MG in a day, CPVm, and CWTm are the cost of PVs and WTs. P m , t P V is the power output of PVs in the m-th MG, at t-th hour, C m , t P V is the PV power generation cost, P m , t W T is the power output of WTs, C m , t W T is the WT power generation cost. Cevm is the cost of energy exchange between MMGS and EVs, C m , t C E V is the charging price of EVs in m-th MG, C m , t D E V is the discharging price, ∆t = 1 h, T = 24 h. Cgm is the energy exchange cost between the MG and the distribution network through the DC/AC converters, P m , t G is the active power output between the MG and the distribution network through the DC/AC converters. If P m , t G ≥ 0, MG purchases electricity from the distribution network. If P m , t G < 0, MMGS sells electricity to the distribution network. C m , t G is the electricity price that MG purchases/sells to the distribution network. Ccym is the additional cycle cost of EV batteries, C c y n E V is the additional battery charging/discharging cycle cost of n-th EV, kcy is the number of additional charging/discharging cycles, N is the number of EVs.

4.4.2. Carbon Emissions and Cost

The electricity of the distribution network mainly comes from thermal power generation. When MG exchanges energy with the distribution network, the distribution network emits more CO2. To reduce carbon emissions as much as possible and increase the use of RESs, in this paper, the cost of carbon emissions is used as the penalty cost of CO2 generated by the energy exchange between the MMGS and the distribution network [29].
C c i m = E C I m k c
E C I m = t = 1 T P m , t G Δ t e c

4.5. Constraints of the Inner-Loop Model

4.5.1. EVs Power Constraint

The charging/discharging power of EVs cannot exceed the rated power of EVCDIs.
| P m , n , t E V | P m , n , R E V C D I s
where P m , n , R E V C D I s is the rated power of the EVCDI serving the n-th EV.

4.5.2. EVs Capacity Constraint

The remaining power of EVs must meet the constraints of rated capacity.
S O C E V m , n , min S O C E V m , n , t S O C E V m , n , max
where SOCEVm,n,min, and SOCEVm,n,max are the minima and maximum capacity, respectively, of n-th EV in m-th MG.

4.5.3. RESs Output Constraint

Considering the performance limitations of renewable energy, the output of RESs in m-th MG has a certain upper limit.
0 P m , t W T P m , max W T
0 P m , t P V P m , max P V

4.5.4. System Power Balance Constraint

For MMGS, the active power output should meet the power balance constraint.
P m , t E V + P m , t G + P m , t W T + P m , t P V = P m , t L
where P m , t L is the total load of the m-th MG at time t.

4.6. Outer-Loop Optimization

The outer-loop optimization model takes the network loss as the optimization goal. By optimizing the reactive power output of the DC/AC converters Q m , t G , the daily network loss of the distribution network is minimized, thereby reducing network loss cost and carbon emissions of MMGS and the distribution network [30]. This paper assumes that the m-th MG is connected to node i of the distribution network.
f 2 = min W S G
W S G = t = 1 T i , j = 1 N b r k i R i j P i j , t 2 + Q i j , t 2 V i j , t 2
W S I = W S G W S B
P i j , t = P i j , t 0 + P m , t G
Q i j , t = Q i j , t 0 + Q m , t G
W S I = f 2 W S B
C i l + C c o = W S I Δ t k i l + W S I Δ t e c k c
W S I is the daily operating increased network loss of the distribution network. Nbr is the number of branches. i, j are the nodes, ki is the state variable of the i-th branch switch, 1 means closed, 0 means open; Rij is the resistance of branch ij, Pij,t, Qij,t are the active and reactive power of branch ij in the t-th hour, Vij,t is the voltage,   P i j , t 0 ,   Q i j , t 0 are initially active, reactive power when connected without MMGS. P m , t G , Q m , t G are the active and reactive power through the DC/AC converters injected into node i by the m-th MG connected to node i. To facilitate the calculation of network loss, a day is divided into 12 small periods, with a time interval of 2 h.

4.6.1. Network Loss Cost

The operation of MMGS connected to the distribution network will cause increased network loss in the distribution network. Therefore, the distribution network will sign a contract with the operator of MMGS, and the operator needs to pay a certain network loss fee for the daily operating increased network loss in the distribution network.
C i l = W S I Δ t k i l

4.6.2. Carbon Emissions and Cost

When the network loss of the distribution network increases by the operation of MMGS, more CO2 will be emitted. MMGS will still incur a penalty cost for carbon emissions by the increasing network loss, which differs from the carbon emissions cost due to energy exchange in the inner-loop model.
C c o = W S I Δ t e c k c

4.7. Constraints of the Outer-Loop Model

The model takes the actual power flow of the power grid as the constraints.

4.7.1. Node Power Flow Constraint

P m , t G + P i , t 0 = P L i , t + V i , t j = 1 N n V j , t ( G i j cos δ i j + B i j sin δ i j )
Q m , t G + Q i , t 0 = Q L i , t + V i , t j = 1 N n V j , t ( G i j sin δ i j + B i j cos δ i j )
where P i , t 0 and Q i , t 0 are the initial input active and reactive power of node i in the t-th hour, PLi,t and QLi,t are the active and reactive load, Vi,t and Vj,t are the voltage of node i and j, Gij, Bij, and δij are the conductance, susceptance, and phase angle difference of branch ij.

4.7.2. Node Voltage Constraint

V i min V i , t V i max
V i min and V i max are lower and upper limits of the node i voltage amplitude.

4.7.3. Branch Power Constraint

| P i j , t | P i j , max
| Q i j , t | Q i j , max
Pij,max, Qij,max are the maximum active and reactive power of the branch ij.

4.7.4. Branch Current Constraint

I i j I i j max
where I i j max is the upper limit of branch ij current carrying capacity.

4.7.5. Reactive Output Constraint of DC/AC Converter

The reactive power output of the DC/AC converters must satisfy the constraint:
| Q m , t G | S m 2 ( P m , t G ) 2
where Sm is the rated power of the DC/AC converter in m-th MG, Q m , t G is the reactive power that the DC/AC converter can output to the distribution network, P m , t G is the active power output by DC/AC converter.

4.8. Particle Swarm Algorithm

4.8.1. Procedure of PSO

The steps of particle swarm optimization (PSO) are shown in Figure 5 [31].

4.8.2. Coding

In the inner-loop, the coding about the economic dispatch of RESs, EVs, and DC/AC converters can be represented by a real-valued matrix. k is the index of the particle of the inner-loop. M is the number of MG. T is the dispatching cycle.
I MG k = [ I MG 1 I MG 2 I MGm I MGM ]
I MGm = [ P m , 1 P V P m , 2 P V P m , t P V P m , T P V P m , 1 W T P m , 2 W T P m , t W T P m , T W T P m , 1 G P m , 2 G P m , t G P m , T G P m , 1 E V P m , 2 E V P m , t E V P m , T E V ]
P m , t E V , P m , t P V , P m , t W T , and P m , t G are the power outputs of EVs, PVs, WTs, and DC/AC converters in the m-th MG in the t-th hour, respectively.
In the outer-loop, the coding about the reactive power output by DC/AC converters can be represented by a real-valued matrix. s is the index of the particle of the outer-loop.
O MG s = [ O MG 1 O MG 2 O MGm O MGM ]
O M G m s = [ Q m , 1 G Q m , 2 G Q m , t G Q m , T G ]
Q m , t G is the reactive power output of DC/AC converters in the m-th MG at t-th hour.
However, the dispatch range of the outer-loop variable also changes when the variable of the DC/AC converters changes in the inner-loop. Therefore, a dynamic range adjustment algorithm is added to the outer-loop model.
| Q m , t G | max = S m 2 ( P m , t G ) 2
The inner-loop and outer-loop cooperate to generate the optimal optimization results.

5. Case Study and Discussion

5.1. Case Description

There are 30 EVs concentrated in OBMG/RMG for the charging/discharging service [31]. The dispatching cycle is 24 h. This paper sets up four cases to analyze the optimization model. By using NWPs from Wuhan City, Hubei Province, China in June 2020, a day’s renewable power generation in summer is predicted as the input of the model.

5.1.1. Case 1

In this case, the EVs do not participate in the energy dispatch of the MG. Once they reach the MG, the EVs will be charged until the batteries are fully charged. MMGS does not optimize the reactive power output by DC/AC converters.

5.1.2. Case 2

In this case, after EVs are connected to the MG, they participate in the energy management system of each MG. Once they reach the MG, the energy in the EV battery will be dispatched by the MG’s energy management system until they leave. When the EVs leave the MG at the end of the dispatching, the energy of EVs should be fully charged. This case takes advantage of the across-time energy transmission of EVs in each independent MG, and the optimization of the reactive power output of DC/AC converters is not considered.

5.1.3. Case 3

In this case, only the inner-loop economic dispatch model is used to minimize the total cost of MMGS by optimizing the active power output of RESs, EVs, and DC/AC converters. The ATSET of EV between RMGs and OBMGs is used. However, the reactive power output of DC/AC converters is also not optimized. Case 3 can be used as a reference.

5.1.4. Case 4

In this case, cooperative multi-objective optimization combines the inner-loop economic dispatch model and the outer-loop reactive power optimization model. The ATSET of EV between RMG and OBMG is considered. By optimizing the active power output of RES, EVs, and DC/AC converters, the total daily operating cost of MMGS is reduced. By optimizing the reactive power output of DC/AC converters, the loss of the distribution network is reduced, and the total economic cost of MMGS is reduced synergistically.

5.2. Simulation System Construction

5.2.1. System Introduction

The modified IEEE 33-node system is used to prove the model, whose structure is shown in Figure 6, and its parameters can be obtained from [32]. According to the principle of distribution [32], OBMG and RMG are set at node 19 and node 20, respectively.

5.2.2. Parameters of RESs

According to the principle of renewable energy consumption [27], the RESs installed in each MG and the power generation cost are given in Table 1. The optimization time interval is 1 h, and the optimization cycle is 1 day (24 h).
The daily wind speed, radiation intensity, temperature, and load data are adopted in this area. The renewable energy output and load curves of each MG come from [27].

5.2.3. Parameters of DC/AC Converter

Considering the performance of the DC/AC converters of MMGS, Sm = 1000 kW, the power limit is set as [33]:
0 | P m , t G | 1000 kW
0 | Q m , t G | 1000 kVar

5.2.4. Parameters of EVs

Take a BYD E6 electric vehicle as an example, whose parameters are from [34]. The battery capacity is 80 kWh, and the upper limit of charging and discharging power of EVCDI is 7 kW. The charging and discharging efficiency are all 90% [34]. An EV consumes an average of 8% of electricity per way between RMG and OBMG [31]. The additional battery charging/discharging cycle cost of EV is CNY 50 each time [35]. The minimum power of the battery of EV is not less than 20% [36]. Considering the needs of users, the upper and lower limits for the battery are 100% and 35% [27].

5.2.5. Other Parameters

The time-of-use (TOU) energy prices in RMG/OBMG from [31] are shown in Figure 7. The carbon emissions factor ec is 86.47 g/kWh [29], and the carbon cost factor kc is 0.21 CNY/kg [37]. The loss cost coefficient of the distribution network kil is 0.74 CNY/kWh [38].

5.3. Simulation Results

5.3.1. Inner-Loop Optimization Results

  • Case 1
In this case, when EVs are connected to the MG, they are charged immediately. In case 1, the 24 h curve of RESs, EVs, load, and DC/AC converter active power output in OBMG/RMG is shown in Figure 8. EVs are charged as soon as they reach RMG/OBMG. The active power curve of the DC/AC converter represents the active power output curve of the MG to the distribution network. When it is below the X-axis, it means that the MG sells electric energy to the distribution network. When it is above the X-axis, it means that the MG purchases electric energy from the distribution network. The power curve of EVs has a similar definition. In Figure 8, RMG will allow EVs to be charged at maximum power from 19:00–20:00, and when RESs are insufficient, MEMS will purchase electricity from the distribution network. OBMG is also charging EVs at 9:00–10:00. The total daily operating cost of RMG is CNY 2776.3, and the total daily operating cost of OBMG is CNY 5732.6. Therefore, the total daily operating cost of MMGS is CNY 8508.9.
2.
Case 2
In this case, since EVs can participate in the energy dispatching of independent MGs, their across-time energy transmission is used. When the total generated power of the RESs in the MGs is greater than the load, the MEMS will sell the remaining energy to the distribution network or charge the EVs according to the energy prices. When the total power generation of RES is less than the load, the MEMS will purchase electricity from the distribution network or EVs according to the energy prices. In Figure 9, the active power output of RES, EV and DC/AC converters in OBMG and MG are optimized. In OBMG, due to the high energy prices of the distribution network and EVs from 9:00 to 12:00, MEMS choose to let EVs release electric energy. OBMG lowers costs by selling energy to the distribution network. When energy prices are low between 12:00 and 15:00, MEMS fully charges EVs. In RMG, MEMS chooses to charge EVs at 23:00 when energy prices are low. This is to avoid additional battery charge–discharge cycle costs due to discharge, so EVs are only charged. The across-time energy transmission of the EV in the independent MG is fully utilized. Through optimization model calculation, the total daily operating cost of RMG is CNY 2644.1, and the total daily operating cost of OBMG is CNY 5642.1. Therefore, the total daily operating cost of MMGS is CNY 8286.2. Compared with Case 1, the across-time energy transmission of EVs can reduce the overall operating cost of MMGS.
3.
Case 3
In this case, EVs can transfer energy among multiple MGs, and EVs participate in MMGS energy dispatching. For OBMG, the energy price of the distribution network and the discharging price of EVs are both high, and the difference between the energy price of the distribution network and the discharging price of EV is much higher than that of RMG. Therefore, MMGS’s energy management system will discharge almost all EVs as much as possible when EVs are connected to OBMG, and earn more profits. For RMG, its advantage is that the charging price is lower, so MMGEMS will try its best to allow almost all EVs to be charged during the low energy price of RMG to reduce the charging cost of EVs. In Figure 10, EVs are discharged as much as possible in OBMG and then charged as much as possible in RMG. After optimization model calculation, the total daily operating cost of RMG is CNY 2391.8, and the total daily operating cost of OBMG is CNY 5404.1. Therefore, the total daily operating cost of MMGS is CNY 7795.9. However, compared with case 1 case 2, by using the across-time-and-space energy transmission of EVs, the total daily operating cost of the MMGS is the lowest in this case.
Table 2 is the comparison of the results of the inner-loop economic dispatch in the three cases. Since case 4 and case 3 use the same inner-loop economic dispatch model, their inner-loop output conditions are the same. Here, the effectiveness of the inner-loop economic dispatch model is mainly discussed, so there is no need to show the results in case 4.
In case 1, EVs do not participate in the energy dispatching of the MG, and MMGS has the highest total operating cost. In case 2, the across-time energy transmission of EVs in the independent MG is used to reduce the cost. In case 3 and case 4, the across-time-and-space energy transmission of EVs is considered to further reduce the total daily operating cost of MMGS, which achieves the lowest daily operating cost COTC.
Figure 11 is the remaining capacity curve of EVs, which proves that EVs meet the power constraint in the four cases. It is also verified that the charging and discharging behaviors analyses of EVs in the three cases are correct.
Table 3 is the cost of EVs’ users. Among the three cases, the user cost of case 3 is the lowest. In case 1, EVs do not participate in dispatching, and the cost of users is the highest. In case 2, the cost of users is reduced by the across-time energy transmission. In case 3, the inner-loop economic dispatch is adopted, which makes full use of the across-time-and-space energy transmission of EVs. Additionally, the cost of users is further reduced. Combining with the lowest daily operating cost of MMGS, the inner-loop economic dispatch model using ATSET of EVs achieved a win–win situation for MMGS and EVs’ users.

5.3.2. Outer-Loop Optimization Results

In case 1, case 2, and case 3, the reactive power output of the DC/AC converters is not optimized. In case 4, the outer-loop reactive power optimization model is used to optimize the reactive power output of the DC/AC converters. The optimized reactive power output of the DC/AC converters of RMG and OBMG in case 4 is shown in Figure 12. The converters will absorb or output a certain amount of reactive power to the distribution network at every moment, which is used to optimize the operating network loss of the distribution network, thereby reducing the energy loss cost CWTC and total carbon emissions EC of MMGS and the distribution network, and cooperating with the inner-loop model to reduce the total economic cost CETC of MMGS. The comparison of the results under the four cases is shown in Table 4.
By analyzing the distribution network loss under the above different cases, it can be concluded that the reactive power output of the DC/AC converters to the distribution network will affect the distribution network loss. When MMGS is not integrated into the distribution network to work, the original baseline loss W S B of the distribution network is 14,787.9 kW. The distribution network loss under the first three cases is all greater than W S B , while the distribution network loss under case 4 is less than W S B and lower than the first three cases. Case 3 and case 4 are a set of comparisons. Under the common premise of using the inner-loop optimization model, case 4 that uses reactive power optimization has lower network loss. Figure 13 is the increased network loss diagram for each period of the distribution network which further proves that intelligently optimizing the reactive power output of DC/AC converters through the outer-loop model can effectively reduce the daily network loss of the distribution network.
Table 5 is the network loss cost and energy loss cost in four cases. Among the four cases, the network loss cost Cil and the carbon emissions cost Cco derived from the optimization of the outer-loop model are the lowest, which proves that the outer-loop optimization model plays a role in the cooperative optimization of the economic cost of MMGS.

5.3.3. Cooperative Multi-objective Optimization Results

It can be concluded from Table 6 that, under the cooperative multi-objective optimization model, the total daily economic cost CETC of MMGS is the lowest. The cost of case 4 adopting the cooperative multi-objective model is 16.3% lower than that for case 1, 13.9% lower than that for case 2, 8.6% lower than that for case 3 which only uses the economic dispatch model of the inner-loop without optimizing reactive power output of DC/AC converters. It is proved that the cooperative multi-objective optimization model improves the economy of MMGS.
It can be concluded by analyzing the carbon emissions data in Table 7 that the total carbon emissions of the MMGS and distribution network with cooperative multi-objective optimization are the lowest among the four cases, which is 24.0% lower than that for case 1, 24.6% lower than that for case 2, and 29.8% lower than that for case 3, which does not optimize the reactive power. The economic cost of MMGS, the network loss of the distribution network, and the total carbon emission of MMGS and the distribution network were all optimized, which fully proves that the cooperative multi-objective optimization achieved the effect.

5.3.4. Further Verification

To further verify the correctness and validity of the model, the weather type of a certain day in the winter of December 2020 in Wuhan City, Hubei Province, China was used as the input of the model. The results of the economic dispatch of the inner-loop model are shown in Table 8, and the results of the network loss optimization of the outer-loop model are shown in Table 9. The network loss cost and energy loss cost in four cases are shown in Table 10. The final economic cost of MMGS in the four cases on another day is shown in Table 11.
It can be concluded from the above tables that the multi-objective optimization of the model is still achieved after using the weather data of one day in winter. The optimal total economic cost of MMGS CETC and the lowest distribution network loss W S G are obtained, which further proves the correctness and effectiveness of the model.

6. Conclusions

A cooperative multi-objective optimization strategy for MMGS containing EVs and RESs is proposed, including dynamic economic dispatch and optimization of reactive power output by DC/AC converters. Dynamic economic dispatch optimizes the active power output of RESs, EVs, and DC/AC converters in MMGS to obtain the optimal daily operating cost of MMGS. Reactive power optimization reduces the daily operating network loss of the distribution network by optimizing the reactive power output of the DC/AC converters to the distribution network. By comparing the results of the four cases, the following conclusions are drawn:
  • According to the simulation results, the economic dispatch model of the inner-loop in the cooperative multi-objective optimization can reduce the operating cost of MMGS, which makes full use of the ATSET of EVs. Additionally, the optimization of the output reactive power output of the DC/AC converters of the outer-loop can reduce network loss cost and carbon emissions cost of the distribution network. The two cooperate to realize the improvement of the economy of MMGS and the efficient operation of the distribution network.
  • The cooperative multi-objective optimization model not only realizes the optimization of the economic cost of MMGS and the network loss of the distribution network, but also reduces the total carbon emissions of MMGS and the distribution network, which greatly responds to the calls for national carbon neutrality and carbon peak.

Author Contributions

Conceptualization, Z.X. and C.C.; methodology, Z.X.; software, Z.X.; validation, Z.X., M.D. and J.Z.; formal analysis, Z.X.; investigation, D.H.; resources, H.C.; data curation, Z.X.; writing—original draft preparation, Z.X.; writing—review and editing, Z.X. and C.C.; visualization, Z.X.; supervision, C.C.; project administration, C.C.; funding acquisition, C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51977086.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

BijSusceptance of branch ij in the distribution network.
CETCEconomic total cost of MMGS.
COTCOperating total cost of MMGS from the inner-loop model.
CWTCEnergy loss cost of the MMGS from the outer-loop model.
CilNetwork loss cost.
CcoCarbon emissions cost.
COCmOperating cost of the m-th MG from the inner-loop model.
CexmEnergy transaction cost in the m-th MG.
CcimCarbon emissions cost in the m-th MG.
CresmCost of RESs of the m-th MG.
CPVmCost of PVs in the m-th MG.
CWTmCost of WTs in the m-th MG.
C m , t P V PV power generation cost in the m-th MG in the t-th hour.
C m , t W T WT power generation cost in the m-th MG in the t-th hour.
CevmCost of energy exchange between the m-th MG and EVs.
C m , t C E V Charging price of EVs in m-th MG in the t-th hour.
C m , t D E V Discharging price of EVs in m-th MG in the t-th hour.
CgmEnergy exchange cost between the m-th MG and the distribution network.
C m , t G Electricity price that m-th MG purchases/sells to the distribution network in the t-th hour.
CcymAdditional cycle cost of EV batteries in m-th MG.
C c y n E V Additional battery charging/discharging cycle cost of n-th EV.
ECTotal carbon emissions of MMGS and the distribution network.
ECimCarbon emissions generated by m-th MG in the inner-loop model.
ECOCarbon emissions in the outer-loop model.
ecCarbon emissions factor.
fThe main objective function of the cooperative optimization model.
f1The objective functions of the inner-loop model.
f2The objective functions of the outer-loop model.
GijConductance of branch ij in the distribution network.
I i j max Upper limit of branch ij current carrying capacity in the distribution network.
i,jNodes of the distribution network.
kilLoss cost coefficient.
kcCarbon cost factor.
kcyNumber of additional charging/discharging cycles.
kiThe state variable of the i-th branch switch.
MNumber of MGs in the MMGS.
NNumber of EVs in m-th MG.
NbrNumber of branches in the distribution network.
P m , n , t E V Exchanging power in the t-th hour of the n-th EV in the m-th MG.
P m , t E V Exchanging power in the t-th hour of the EVs in the m-th MG.
P m , n , t E V C D I s Power of the EVCDIs of the n-th EV in the m-th MG in the t-th hour.
P m , t P V Power output of PVs in the m-th MG in the t-th hour.
P m , t W T power output of WTs in the m-th MG in the t-th hour.
P m , t G Active power output between the m-th MG and the distribution network in the t-th hour through the DC/AC converters.
P m , n , R E V C D I s Rated power of the EVCDI serving the n-th EV in the m-th MG.
P m , t L Total load of the m-th MG in the t-th hour.
Pij,tActive power of branch ij in the t-th hour.
P i j , t 0 Initially active power of branch ij when connected without MG in the t-th hour.
P i , t 0 Initial input active power of node i in the t-th hour.
PLi,tActive load of node i in the t-th hour.
Pij,maxMaximum active power of the branch ij.
Qij,tReactive power of branch ij in the t-th hour.
Q i j , t 0 Initially reactive power of branch ij when connected without MMGS in the t-th hour.
Q m , t G Reactive power output between the m-th MG and the distribution network in the t-th hour through the DC/AC converters.
Q i , t 0 Initial input reactive power of node i in the t-th hour.
QLi,tReactive load of node i in the t-th hour.
Qij,maxMaximum reactive power of the branch ij.
RijThe resistance of branch ij.
SOCEVm,n,tRemaining power capacity of the n-th EV in the m-th MG in the t-th hour.
SOCEVm,n,minMinima capacity, respectively, of the n-th EV in the m-th MG.
SOCEVm,n,maxMaximum capacity, respectively, of the n-th EV in the m-th MG.
SmRated power of the DC/AC converter in the m-th MG.
TScheduling cycle, one day, 24 h.
Vij,tVoltage of branch ij in the t-th hour.
V i min Lower limits of the node i voltage amplitude.
V i max Upper limits of the node i voltage amplitude.
W S B Original baseline network loss.
W S G Total daily operating network loss of the distribution network.
W S I Daily operating increased network loss of the distribution network.
σ Self-discharge coefficient of EV’s battery.
tLength of the time slot set for the optimization.
ηDEVEfficiency for EV discharging.
ηCEVEfficiency for EV charging.
δijPhase angle difference of branch ij.

References

  1. Chan, C.; Zhou, G.Y.; Zhang, D. Intelligent energy ecosystem based on carbon neutrality. In Proceedings of the 2017 IEEE Conference on Energy Internet and Energy System Integration (EI2), Beijing, China, 26–28 November 2017; pp. 1–5. [Google Scholar]
  2. Wang, X.; Barnett, A. The evolving value of photovoltaic module efficiency. Appl. Sci. 2019, 9, 1227. [Google Scholar] [CrossRef] [Green Version]
  3. Li, H.; Ye, Y.; Lin, L. Low-Carbon Economic Bi-Level Optimal Dispatching of an Integrated Power and Natural Gas Energy System Considering Carbon Trading. Appl. Sci. 2021, 11, 6968. [Google Scholar] [CrossRef]
  4. Izquierdo-Monge, O.; Peña-Carro, P.; Villafafila-Robles, R.; Duque-Perez, O.; Zorita-Lamadrid, A.; Hernandez-Callejo, L. Conversion of a Network Section with Loads, Storage Systems and Renewable Generation Sources into a Smart Microgrid. Appl. Sci. 2021, 11, 5012. [Google Scholar] [CrossRef]
  5. Saponara, S.; Saletti, R.; Mihet-Popa, L. Hybrid micro-grids exploiting renewables sources, battery energy storages, and bi-directional converters. Appl. Sci. 2019, 9, 4973. [Google Scholar] [CrossRef] [Green Version]
  6. Van Der Meer, D.; Mouli, G.R.C.; Mouli, G.M.-E.; Elizondo, L.R.; Bauer, P. Energy management system with PV power forecast to optimally charge EVs at the workplace. IEEE Trans. Ind. Inform. 2018, 14, 311–320. [Google Scholar] [CrossRef] [Green Version]
  7. Puma-Benavides, D.S.; Izquierdo-Reyes, J.; Calderon-Najera, J.d.D.; Ramirez-Mendoza, R.A. A Systematic Review of Technologies, Control Methods, and Optimization for Extended-Range Electric Vehicles. Appl. Sci. 2021, 11, 7095. [Google Scholar] [CrossRef]
  8. Ko, K.S.; Han, S.; Sung, D.K. Performance-based settlement of frequency regulation for electric vehicle aggregators. IEEE Trans. Smart Grid 2018, 9, 866–875. [Google Scholar] [CrossRef]
  9. Alam, M.J.E.; Muttaqi, K.M.; Sutanto, D. Effective utilization of available PEV battery capacity for mitigation of solar PV impact and grid support with integrated V2G functionality. IEEE Trans. Smart Grid 2016, 7, 1562–1571. [Google Scholar] [CrossRef] [Green Version]
  10. Kisacikoglu, M.C.; Erden, F.; Erdogan, N. Distributed control of PEV charging based on energy demand forecast. IEEE Trans. Ind. Inform. 2018, 14, 332–341. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, H.; Pan, H.; Luo, F.; Qiu, J.; Deng, Y.; Lai, M.; Dong, Z.Y. Operational planning of electric vehicles for balancing wind power and load fluctuations in a microgrid. IEEE Trans. Sustain. Energy 2017, 8, 592–604. [Google Scholar] [CrossRef]
  12. Panwar, L.K.; Konda, S.R.; Verma, A.; Panigrahi, B.K.; Kumar, R. Operation window constrained strategic energy management of microgrid with electric vehicle and distributed resources. IET Gener. Transm. Distrib. 2017, 11, 615–626. [Google Scholar] [CrossRef]
  13. Masrur, M.A.; Skowronska, A.G.; Hancock, J.; Kolhoff, S.W.; McGrew, D.Z.; Vandiver, J.C.; Gatherer, J. Military-based vehicle-to-grid and vehicle-to-vehicle microgrid—System architecture and implementation. IEEE Trans. Transp. Electrif. 2017, 4, 157–171. [Google Scholar] [CrossRef]
  14. Zou, H.; Mao, S.; Wang, Y.; Zhang, F.; Chen, X.; Cheng, L. A survey of energy management in interconnected multi-microgrids. IEEE Access 2019, 7, 72158–72169. [Google Scholar] [CrossRef]
  15. Nikmehr, N.; Ravadanegh, S.N. Optimal power dispatch of multi-microgrids at future smart distribution grids. IEEE Trans. Smart Grid 2015, 6, 1648–1657. [Google Scholar] [CrossRef]
  16. Kaloudas, C.G.; Ochoa, L.F.; Marshall, B.; Majithia, S.; Fletcher, I. Assessing the future trends of reactive power demand of distribution networks. IEEE Trans. Power Syst. 2017, 32, 4278–4288. [Google Scholar] [CrossRef] [Green Version]
  17. Ganguly, S. Multi-objective planning for reactive power compensation of radial distribution networks with unified power quality conditioner allocation using particle swarm optimization. IEEE Trans. Power Syst. 2014, 29, 1801–1810. [Google Scholar] [CrossRef]
  18. Syed, I.; Khadkikar, V.; Zeineldin, H.H. Loss reduction in radial distribution networks using a solid-state transformer. IEEE Trans. Ind. Appl. 2018, 54, 5474–5482. [Google Scholar] [CrossRef]
  19. Bolognani, S.; Zampieri, S. A distributed control strategy for reactive power compensation in smart microgrids. IEEE Trans. Autom. Control 2013, 58, 2818–2833. [Google Scholar] [CrossRef] [Green Version]
  20. Arefifar, S.A.; Mohamed, Y.A.-R.I. Probabilistic optimal reactive power planning in distribution systems with renewable resources in grid-connected and islanded modes. IEEE Trans. Ind. Electron. 2014, 61, 5830–5839. [Google Scholar] [CrossRef]
  21. Cho, Y.-W.; Kwon, J.-M.; Kwon, B.-H. Single power-conversion AC--DC converter with high power factor and high efficiency. IEEE Trans. Power Electron. 2014, 29, 4797–4806. [Google Scholar] [CrossRef]
  22. Sayed, M.A.; Takeshita, T.; Kitagawa, W. Advanced PWM switching technique for accurate unity power factor of bidirectional three-phase grid-tied DC–AC converters. IEEE Trans. Ind. Appl. 2019, 55, 7614–7627. [Google Scholar] [CrossRef]
  23. Chen, C.; Duan, S.; Cai, T.; Liu, B. Online 24-h solar power forecasting based on weather type classification using artificial neural network. Sol. Energy 2011, 85, 2856–2870. [Google Scholar] [CrossRef]
  24. Aien, M.; Fotuhi-Firuzabad, M.; Rashidinejad, M. Probabilistic optimal power flow in correlated hybrid wind–photovoltaic power systems. IEEE Trans. Smart Grid 2014, 5, 130–138. [Google Scholar] [CrossRef]
  25. Quan, H.; Srinivasan, D.; Khosravi, A. Short-term load and wind power forecasting using neural network-based prediction intervals. IEEE Trans. Neural Netw. Learn. Syst. 2014, 25, 303–315. [Google Scholar] [CrossRef]
  26. Luo, Z.; Song, Y.; Hu, Z.; Xu, Z.; Yang, X.; Zhan, K. Forecasting charging load of plug-in electric vehicles in China. In Proceedings of the 2011 IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011; pp. 1–8. [Google Scholar]
  27. Chen, C.; Xiao, L.; Duan, S.; Chen, J. Cooperative optimization of electric vehicles in microgrids considering across-time-and-space energy transmission. IEEE Trans. Ind. Electron. 2019, 66, 1532–1542. [Google Scholar] [CrossRef]
  28. Chen, C.; Duan, S. Optimal integration of plug-in hybrid electric vehicles in microgrids. IEEE Trans. Ind. Inform. 2014, 10, 1917–1926. [Google Scholar] [CrossRef]
  29. Li, Y.; Zou, Y.; Tan, Y.; Cao, Y.; Liu, X.; Shahidehpour, M.; Tian, S.; Bu, F. Optimal stochastic operation of integrated low-carbon electric power, natural gas, and heat delivery system. IEEE Trans. Sustain. Energy 2018, 9, 273–283. [Google Scholar] [CrossRef]
  30. Zhang, T.; Yu, L. Optimal allocation of DSTATCOM considering the uncertainty of photovoltaic systems. IEEJ Trans. Electr. Electron. Eng. 2020, 15, 355–363. [Google Scholar] [CrossRef]
  31. Chen, C.; Chen, J.; Wang, Y.; Duan, S.; Cai, T.; Jia, S. A price optimization method for microgrid economic operation considering across-time-and-space energy transmission of electric vehicles. IEEE Trans. Ind. Inform. 2020, 16, 1873–1884. [Google Scholar] [CrossRef]
  32. Saha, S.; Banerjee, S.; Maity, D.; Chanda, C. Optimal sizing and location determination of distributed generation in distribution networks. In Proceedings of the 2015 International Conference on Energy, Power and Environment: Towards Sustainable Growth (ICEPE), Shillong, India, 12–13 June 2015; pp. 1–5. [Google Scholar]
  33. Yang, J.; He, Z.; Ke, J.; Xie, M. A new hybrid multilevel DC–AC converter with reduced energy storage requirement and power losses for HVDC applications. IEEE Trans. Power Electron. 2018, 34, 2082–2096. [Google Scholar] [CrossRef]
  34. Chen, J.; Chen, C.; Duan, S. Cooperative optimization of electric vehicles and renewable energy resources in a regional multi-microgrid system. Appl. Sci. 2019, 9, 2267. [Google Scholar] [CrossRef] [Green Version]
  35. Bocca, A.; Baek, D. Optimal Life-Cycle Costs of Batteries for Different Electric Cars. In Proceedings of the 2020 AEIT International Conference of Electrical and Electronic Technologies for Automotive (AEIT AUTOMOTIVE), Turin, Italy, 18–20 November 2020; pp. 1–6. [Google Scholar]
  36. Yang, H.; Zhou, W.; Lu, L.; Fang, Z. Optimal sizing method for stand-alone hybrid solar–wind system with LPSP technology by using genetic algorithm. Sol. Energy 2008, 82, 354–367. [Google Scholar] [CrossRef]
  37. Xu, K.; Wang, S.; Zhao, J.; Dou, X.; Zhang, P. Research on Microgrid Dispatch Considering the Charge and Discharge Strategy of Electric Vehicles. In Proceedings of the 2020 12th IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC), Nanjing, China, 20–23 September 2020; pp. 1–5. [Google Scholar]
  38. Li, Z.; Lu, Q.; Fu, Y. State-splitting multi-objective dynamic programming algorithm for dynamic reconfiguration of active distribution network. Proc. CSEE 2019, 39, 5025–5036. [Google Scholar]
Figure 1. Structure of a DC multi-microgrid system.
Figure 1. Structure of a DC multi-microgrid system.
Applsci 11 08916 g001
Figure 2. Structure of a DC multi-microgrid control system.
Figure 2. Structure of a DC multi-microgrid control system.
Applsci 11 08916 g002
Figure 3. Structure of a DC microgrid.
Figure 3. Structure of a DC microgrid.
Applsci 11 08916 g003
Figure 4. Process of the cooperative multi-objective optimization.
Figure 4. Process of the cooperative multi-objective optimization.
Applsci 11 08916 g004
Figure 5. Process of the PSO algorithm.
Figure 5. Process of the PSO algorithm.
Applsci 11 08916 g005
Figure 6. Topological diagram of the modified IEEE 33-node system.
Figure 6. Topological diagram of the modified IEEE 33-node system.
Applsci 11 08916 g006
Figure 7. (a) Prices of energy exchanging in OBMG; (b) prices of energy exchanging in RMG.
Figure 7. (a) Prices of energy exchanging in OBMG; (b) prices of energy exchanging in RMG.
Applsci 11 08916 g007
Figure 8. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Figure 8. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Applsci 11 08916 g008
Figure 9. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Figure 9. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Applsci 11 08916 g009
Figure 10. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Figure 10. (a) Power output of RESs, EVs, and the DC/AC converter in OBMG; (b) power output of RESs, EVs, and the DC/AC converter in RMG.
Applsci 11 08916 g010
Figure 11. Remaining capacity of EVs in three cases.
Figure 11. Remaining capacity of EVs in three cases.
Applsci 11 08916 g011
Figure 12. (a) Reactive power output by the DC/AC converter of OBMG in case 4; (b) reactive power output by the DC/AC converter of RMG in case 4.
Figure 12. (a) Reactive power output by the DC/AC converter of OBMG in case 4; (b) reactive power output by the DC/AC converter of RMG in case 4.
Applsci 11 08916 g012
Figure 13. Increased network loss in four cases.
Figure 13. Increased network loss in four cases.
Applsci 11 08916 g013
Table 1. Installed RESs and the power generation cost.
Table 1. Installed RESs and the power generation cost.
MG TypeRES TypeInstalled Capacity/kWPower Generation Cost/¥·kWh−1
OBMGPV18000.24
WT18000.38
RMGPV24000.24
WT24000.38
Table 2. Comparison of operating cost in three cases.
Table 2. Comparison of operating cost in three cases.
CaseMGsCex
Energy Transaction Cost/CNY
Cci
Carbon Emissions Cost/CNY
COTC
Total Operating Cost/CNY
Case 1RMG2760.016.32776.3
OBMG5700.831.85732.6
MMGS8460.848.18508.9
Case 2RMG2627.816.32644.1
OBMG5609.632.55642.1
MMGS8237.448.88286.2
Case 3RMG2348.942.92391.8
OBMG5394.69.55404.1
MMGS7743.552.47795.9
Table 3. The cost of EVs’ users in three cases in one day.
Table 3. The cost of EVs’ users in three cases in one day.
CaseCase 1Case 2Case 3
The Cost of EVs’ Users/CNY211.2173.4−410.8
Table 4. Comparison of network loss in the four cases.
Table 4. Comparison of network loss in the four cases.
CaseCase 1Case 2Case 3Case 4
W S G
Total Network Loss/kW
14,889.214,872.514,876.013,987.2
W S I
Increased Network Loss/kW
101.384.688.1−800.7
Table 5. Network loss cost and energy loss cost in the four cases.
Table 5. Network loss cost and energy loss cost in the four cases.
CaseCil
Network Loss Cost/CNY
Cco
Carbon Emissions Cost/CNY
CWTC
Energy Loss Cost/CNY
Case 175.01.876.8
Case 262.61.564.1
Case 365.21.666.8
Case 4−592.5−14.5−607.0
Table 6. The final economic cost of MMGS in four cases.
Table 6. The final economic cost of MMGS in four cases.
CaseCase 1Case 2Case 3Case 4
CETC
Total Economic Cost of MMGS/CNY
8585.78350.37862.77188.9
Table 7. Total carbon emissions in the four cases.
Table 7. Total carbon emissions in the four cases.
CaseCase 1Case 2Case 3Case 4
EC
Total Carbon Emissions/kg
237.6239.5257.1180.5
Table 8. Comparison of operating cost in three cases on another day.
Table 8. Comparison of operating cost in three cases on another day.
CaseMGsCex
Energy Transaction Cost/CNY
Cci
Carbon Emissions Cost/CNY
COTC
Total Operating Cost/CNY
Case 1RMG2968.923.42992.3
OBMG6446.448.66495.0
MMGS9415.372.09487.3
Case 2RMG2836.723.42860.1
OBMG6304.544.66349.1
MMGS9141.268.09209.2
Case 3RMG2574.448.52622.9
OBMG6154.227.66181.8
MMGS8728.676.18804.7
Table 9. Comparison of network loss in the four cases on another day.
Table 9. Comparison of network loss in the four cases on another day.
CaseCase 1Case 2Case 3Case 4
W S G
Total Network Loss/kW
14,933.114,915.814,919.414,110.8
W S I
Increased Network Loss/kW
145.2127.9131.5−677.1
Table 10. Network loss cost and energy loss cost on four cases on another day.
Table 10. Network loss cost and energy loss cost on four cases on another day.
CaseCil
Network Loss Cost/CNY
Cco
Carbon Emissions Cost/CNY
CWTC
Energy Loss Cost/CNY
Case 1107.42.6110.0
Case 294.62.396.9
Case 397.32.499.7
Case 4−501.1−12.3−513.4
Table 11. The final economic cost of MMGS in the four cases on another day.
Table 11. The final economic cost of MMGS in the four cases on another day.
CaseCase 1Case 2Case 3Case 4
CETC
Total Economic Cost of MMGS/CNY
9597.39306.18904.48291.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, Z.; Chen, C.; Dong, M.; Zhang, J.; Han, D.; Chen, H. Cooperative Multi-Objective Optimization of DC Multi-Microgrid Systems in Distribution Networks. Appl. Sci. 2021, 11, 8916. https://doi.org/10.3390/app11198916

AMA Style

Xu Z, Chen C, Dong M, Zhang J, Han D, Chen H. Cooperative Multi-Objective Optimization of DC Multi-Microgrid Systems in Distribution Networks. Applied Sciences. 2021; 11(19):8916. https://doi.org/10.3390/app11198916

Chicago/Turabian Style

Xu, Zhiwen, Changsong Chen, Mingyang Dong, Jingyue Zhang, Dongtong Han, and Haowen Chen. 2021. "Cooperative Multi-Objective Optimization of DC Multi-Microgrid Systems in Distribution Networks" Applied Sciences 11, no. 19: 8916. https://doi.org/10.3390/app11198916

APA Style

Xu, Z., Chen, C., Dong, M., Zhang, J., Han, D., & Chen, H. (2021). Cooperative Multi-Objective Optimization of DC Multi-Microgrid Systems in Distribution Networks. Applied Sciences, 11(19), 8916. https://doi.org/10.3390/app11198916

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