# Dual-Layer Optimal Dispatching Strategy for Microgrid Energy Management Systems considering Demand Response

Authors: Tinglong Pan; Hui Liu; Dinghui Wu; Zeliang Hao
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

The continuous development of microgrid’s technology creates favorable conditions for the access of distributed energy. Firstly, in order to consider the interests of the demand side and the power side, this paper presents a dual-layer optimal dispatching model of microgrid based on demand response. The objective of the first-layer optimization is to obtain the maximum load satisfaction and to optimize the load curve. The objective of the second-layer optimization is to make the microgrid system economical and environmentally friendly and to optimize the power utilization ratio. And a microsource control strategy based on the isolated microgrid is proposed, which can optimize the operation state of battery and improve the economy of the system. Finally, the Nondominated Sorting Genetic Algorithm-II (NSGA-II) is adopted to solve the optimal scheduling problem of the isolated microgrid. The simulation results indicate that the microsource scheduling strategy proposed in this paper can improve the operation economy and environmental conservation of the system. It can improve the reliability of microgrid power supply and reduce energy waste.


---

## Body

## 1. Introduction

The rational allocation of microsource output plays an important role in ensuring the safe and stable operation of microgrids and improving economic efficiency. It is significant to research on optimization in microgrid energy management. The microgrid energy management system (MGEMS) can coordinate distributed generators, storage battery, and load in the microgrid through the information such as the load demand forecasting. MGEMS not only ensures the security and stability of the microgrid system, but also realizes the economic optimal operation of the microgrid [1]. Therefore, the microgrid modelling, optimization algorithms, and scheduling strategies should be deeply researched to achieve the optimal economic dispatching for energy management and improve the operation efficiency of the microgrid.For microgrid operation, many studies focus on minimizing the operation cost by scheduling different power generation units, which regard the objective function as the prerequisite for the economic operation of microgrid. At present, the study of optimization of the microgrid operation mainly considers its economy and environmental protection including initial investment cost, operation and maintenance cost, sales income of power, government subsidy, and environment benefit. Guo et al. [2] proposed a multiobjective optimization model for isolated microgrid system, which aimed at the confliction of interests between the distribution company and the distributed generation owners in the isolated microgrid system. The economic scheduling objective function in [3, 4] considered the operation costs, environmental pollution of the power system, and the equipment type, while Chen et al. [5] also considered the extra cost of battery. Chen et al. [6] proposed that the economic scheduling problem of isolated microgrid for reducing the total cost of distributed generators and satisfying the demand and supply constraints. Trivedi et al. [7] adopted the whale optimization algorithm to obtain the minimum operation cost and total pollution emissions. Alanazi studied the influence of high penetrability wind energy on the economical and reliable operation of microgrid in reducing the power transmission fluctuation of microgrid [8]. A multiobjective optimal scheduling model is proposed based on the opportunity constraint planning of microgrid, and the uncertainty of renewable energy and load is taken into account [9].The optimal scheduling strategies of microgrid play a crucial role in the scheduling optimization. In [10], a two-stage optimal scheduling method is proposed to solve the problem of optimal scheduling of the microgrid system with the combination of cooling, heat, and power. Duong Tran et al. proposed an economic dispatch method of the microgrid based on time-of-use pricing with coordinating the electricity price and battery charge-discharge management, which adopted the batteries for energy storage to achieve peak shaving during peak period [11]. Duong Tran proposed an energy management system for the battery; the stochastic dynamic programming method is employed to deal with the forecast data to improve the battery life cycle [12]. For the operation characteristics of microgrid, Bie et al. [13] proposed the reliability model of virtual power supply to realize the reliability evaluation of the distribution system with multiple microgrids.Some intelligent optimization algorithms help us better optimizing the economic operation cost of the microgrid, such as genetic algorithm, particle swarm optimization algorithm, ant colony algorithm, and differential evolution algorithm [14–16]. Mohammadi proposed an improved adaptive firefly algorithm to solve the problem of microgrid economy optimization with establishing the operation mathematical model of microgrid [17]. In [4], the K-means clustering method is used to minimize the objectives of economic and environmental during the three stages, and the use of the K-means algorithm reduced the computational time. In this paper, the NSGA-II algorithm with elite retention strategy is used to solve the multiobjective optimization model of microgrid; the algorithm can deal with the multiobjective problem efficiently.At present, most of the researches on microgrid energy management system are single preday scheduling, which can not completely reflect the randomness of new energy generation, and the optimization results may not accord with the actual operation status of microgrid. Most of the researches based on the economic optimization of demand side microgrid only consider the charge and discharge problem of storage battery or the load but do not consider the influence of both. In order to solve the above problems, this paper puts forward a dual-layer optimal dispatching strategy based on the demand side for the independent microgrid system. The simulation results show that under the premise of considering the load satisfaction on the demand side and the influence of battery, using the dual-layer optimal dispatching strategy of microgrid proposed in this paper, the economic, environmental protection and electric energy utilization ratio of the microgrid system can reach an optimal state. The major contributions and innovations of this paper are displayed as follows:(1) A dual-layer optimal dispatching strategy of microgrid is proposed. The first-layer adjusts the load according to the microsource forecast power, adjusts the load output curve according to the load satisfaction, and inputs the corresponding load curve to the second layer.(2) The optimal operation of the microgrid is a multiobjective optimization problem, and the operation of microgrid cannot be completely described by only single objective. Therefore, the multiobjective optimization algorithm is adopted in the second-layer optimization, and the three-target model of the economy, environmental protection, and power utilization rate of the microgrid system is established.(3) The battery charge and discharge strategy is introduced in the algorithm, which can reasonably dispatch batteries and microgas turbines and can improve the utilization ratio of new energy.(4) The current research neglects the diversity of user side and does not fully exploit the potential of user participation in scheduling. In this paper, based on the demand response, the microgrid energy optimization strategy is proposed to realize the economic and environmental optimization under the premise of user satisfaction and make the renewable energy to be used optimally, thus realizing the global optimization of the microgrid.The remainder of this paper is organized as follows. Section2 presents the system description. The optimization strategy of the microgrid is introduced in Section 3. The optimization model solving algorithm and the simulation results are introduced in Sections 4 and 5, respectively. The conclusions are presented in Section 6.


## 2. System Description

In this paper, the microsources of the microgrid system contain wind generators, photovoltaics, microgas turbines, and batteries. Due to the randomness of wind power and photovoltaic power generation, the microsource cannot provide stable output power; in order to solve the problem, we add the microgas turbine power generation as a supplementary power with less emission. The reasonable dispatching of microgas turbine and battery can make the utilization of new energy reach the highest and adopting the reasonable charging-discharging strategy to prolong the life cycle of battery is also the key to solve the reliability of microgrid. The microgrid structure diagram is shown in Figure1.Figure 1
The microgrid structure diagram.


### 2.1. Wind Generator Model

The wind speed-power characteristic curves of wind generator can be obtained by polynomial fitting, which can be expressed as [17](1)PWT=0v<vciav3+bv2+cv+dvci≤v≤vrPrvr<v<vco0v≥vco,where PWT is wind generator output active power; Pr is the rated power of wind generator; vci, vr, and vco are wind generator cut-in wind speed, rated wind speed, and cut-out wind speed, respectively; a, b, c, and d are wind speed parameters, respectively.The environmental parameter correction model is shown as below.The wind speed varies with the height, and the wind speed-power curve needs to provide the wind speed at the hub height of the wind generator, which needs to be converted to the wind speed at the wind generator hub. The conversion formula is shown as(2)v1=v2H1H2α,where v1 and v2 are the wind speed at different heights; H1 and H2 represent different heights, respectively; α is a power-law exponent, which is related to the roughness of the ground.


### 2.2. Photovoltaic Model

In the optimal scheduling of microgrid, it needs to forecast the power of photovoltaic cell. The power output model of the photovoltaic cell is shown as [18](3)PPV=YPVfPVGT¯G¯T,STC1+αPTc-Tc,STC,where PPV is the output active power of the photovoltaic cell; YPV is the output power of photovoltaic cells in the standard test conditions; fPV is PV derating factor, usually 0.8; G¯T is the solar irradiation intensity of the current step in actual environment; G¯T,STC is the solar radiation intensity under standard test conditions; αP is the power temperature coefficient of the PV panel; Tc is the photovoltaic cell temperature of the current time step; Tc,STC is the photovoltaic cell temperature under the standard test condition, usually 25.The PV cell temperature is the temperature of the PV array surface, which is the same as the ambient temperature and is higher than the ambient temperature during the day. According to the principle of energy conservation, the photovoltaic cell temperature model can be expressed as(4)Tc=Ta+GT¯ταUL1-ηcτα,where Ta is the ambient temperature; α is the solar absorption of photovoltaic panels; τ is the solar transmittance; UL is the heat dissipation factor; ηc is the photovoltaic conversion efficiency of photovoltaic cells.


### 2.3. Battery Model

As the energy storage equipment, the battery will release the stored energy during the lack of renewable energy, and it provides stable power supply to users to guarantee the safe and stable operation of the system.In the optimal scheduling of microgrid, the charging-discharging power of the battery is limited by the total residual energy, and in the time step ofΔt, the maximum charging-discharging power of the battery determined by the total remaining energy of the battery can be expressed as(5)Pbatt,dmax=kQ1e-kΔt+Qkc1-e-kΔt1-e-kΔt+ckΔt-1+e-kΔt,(6)Pbatt,cmax=-kcQmax+kQ1e-kΔt+Qkc1-e-kΔt1-e-kΔt+ckΔt-1+e-kΔt,where k is the charging-discharging ratio constant; it reflects the energy of the binding energy and available energy in the battery; Δt is the simulation step; c is the capacity ratio constant; Pbatt,dmax and Pbatt,cmax are the maximum charge power and discharge power of battery determined by the total remaining energy of the battery respectively; Q is the total energy at any time; Q1 is the available energy; and Q2 is the bound energy.The relationship between the state of charge (SOC) and the charge-discharge power of the battery (Pbatt) is expressed as follows [19]:(7)Soct=1-ηbattSoct-1+PbattΔtηbatt,cEbattPbatt<01-ηbattSoct-1-PbattΔtEbattηbatt,dPbatt≥0,where ηbatt is the self-discharge efficiency of the battery; Pbatt is the charge-discharge power of the battery; ηbatt,c and ηbatt,d are the charge and discharge efficiency of the battery, respectively (90%); Δt is the charge-discharge time; Ebatt is the rated capacity of the battery.


### 2.4. Microgas Turbine Model

The gas turbine generates electricity by consuming fuel. The output power of microgas turbine can be freely controlled with fast response speed; this paper mainly dispatches the power of microgas turbine. The output power of the microgas turbine is related to the fuel supply and the low calorific value of the fuel. The operating efficiency of the microgas turbine is related to the output power of its output [20]:(8)ηMT=0.0753PMT653-0.3095PMT652+0.4174PMT65+0.1068,where PMT is the active output power of the microgas turbine; ηMT is the operational efficiency of microgas turbines.


### 2.5. Load Model

For the isolated microgrid system, the load can be divided into the important load power and the interruptible load and the transferable load according to the importance degree, in which the noninterruptible loads are some basic electrical equipment, and the interruptible loads are some alternative electric or recreational equipment. Due to the limitation of power capacity and the intermittence of the renewable energy power generation for the isolated microgrid, it is difficult to ensure the real-time power balance in the microgrid only through the power generation side scheduling. There will be insufficient power supply at some periods, so the supply of part of the loads has to be interrupted. By determining the interruptible load in advance to ensure the continuous and stable operation of the microgrid, it is advantageous to ensure the safety and rationality of the load cutting operation.The load model is expressed as follows:(9)PL=PLim+PLin+Ptran,(10)Pnet-load=PL-PPV-PWT,where PL, PLim, PLin, and Ptran are the daily demand load power of users, the important load power, the load power, and the transferable load power, respectively; Pnet-load is the power of the net load.


## 2.1. Wind Generator Model

The wind speed-power characteristic curves of wind generator can be obtained by polynomial fitting, which can be expressed as [17](1)PWT=0v<vciav3+bv2+cv+dvci≤v≤vrPrvr<v<vco0v≥vco,where PWT is wind generator output active power; Pr is the rated power of wind generator; vci, vr, and vco are wind generator cut-in wind speed, rated wind speed, and cut-out wind speed, respectively; a, b, c, and d are wind speed parameters, respectively.The environmental parameter correction model is shown as below.The wind speed varies with the height, and the wind speed-power curve needs to provide the wind speed at the hub height of the wind generator, which needs to be converted to the wind speed at the wind generator hub. The conversion formula is shown as(2)v1=v2H1H2α,where v1 and v2 are the wind speed at different heights; H1 and H2 represent different heights, respectively; α is a power-law exponent, which is related to the roughness of the ground.


## 2.2. Photovoltaic Model

In the optimal scheduling of microgrid, it needs to forecast the power of photovoltaic cell. The power output model of the photovoltaic cell is shown as [18](3)PPV=YPVfPVGT¯G¯T,STC1+αPTc-Tc,STC,where PPV is the output active power of the photovoltaic cell; YPV is the output power of photovoltaic cells in the standard test conditions; fPV is PV derating factor, usually 0.8; G¯T is the solar irradiation intensity of the current step in actual environment; G¯T,STC is the solar radiation intensity under standard test conditions; αP is the power temperature coefficient of the PV panel; Tc is the photovoltaic cell temperature of the current time step; Tc,STC is the photovoltaic cell temperature under the standard test condition, usually 25.The PV cell temperature is the temperature of the PV array surface, which is the same as the ambient temperature and is higher than the ambient temperature during the day. According to the principle of energy conservation, the photovoltaic cell temperature model can be expressed as(4)Tc=Ta+GT¯ταUL1-ηcτα,where Ta is the ambient temperature; α is the solar absorption of photovoltaic panels; τ is the solar transmittance; UL is the heat dissipation factor; ηc is the photovoltaic conversion efficiency of photovoltaic cells.


## 2.3. Battery Model

As the energy storage equipment, the battery will release the stored energy during the lack of renewable energy, and it provides stable power supply to users to guarantee the safe and stable operation of the system.In the optimal scheduling of microgrid, the charging-discharging power of the battery is limited by the total residual energy, and in the time step ofΔt, the maximum charging-discharging power of the battery determined by the total remaining energy of the battery can be expressed as(5)Pbatt,dmax=kQ1e-kΔt+Qkc1-e-kΔt1-e-kΔt+ckΔt-1+e-kΔt,(6)Pbatt,cmax=-kcQmax+kQ1e-kΔt+Qkc1-e-kΔt1-e-kΔt+ckΔt-1+e-kΔt,where k is the charging-discharging ratio constant; it reflects the energy of the binding energy and available energy in the battery; Δt is the simulation step; c is the capacity ratio constant; Pbatt,dmax and Pbatt,cmax are the maximum charge power and discharge power of battery determined by the total remaining energy of the battery respectively; Q is the total energy at any time; Q1 is the available energy; and Q2 is the bound energy.The relationship between the state of charge (SOC) and the charge-discharge power of the battery (Pbatt) is expressed as follows [19]:(7)Soct=1-ηbattSoct-1+PbattΔtηbatt,cEbattPbatt<01-ηbattSoct-1-PbattΔtEbattηbatt,dPbatt≥0,where ηbatt is the self-discharge efficiency of the battery; Pbatt is the charge-discharge power of the battery; ηbatt,c and ηbatt,d are the charge and discharge efficiency of the battery, respectively (90%); Δt is the charge-discharge time; Ebatt is the rated capacity of the battery.


## 2.4. Microgas Turbine Model

The gas turbine generates electricity by consuming fuel. The output power of microgas turbine can be freely controlled with fast response speed; this paper mainly dispatches the power of microgas turbine. The output power of the microgas turbine is related to the fuel supply and the low calorific value of the fuel. The operating efficiency of the microgas turbine is related to the output power of its output [20]:(8)ηMT=0.0753PMT653-0.3095PMT652+0.4174PMT65+0.1068,where PMT is the active output power of the microgas turbine; ηMT is the operational efficiency of microgas turbines.


## 2.5. Load Model

For the isolated microgrid system, the load can be divided into the important load power and the interruptible load and the transferable load according to the importance degree, in which the noninterruptible loads are some basic electrical equipment, and the interruptible loads are some alternative electric or recreational equipment. Due to the limitation of power capacity and the intermittence of the renewable energy power generation for the isolated microgrid, it is difficult to ensure the real-time power balance in the microgrid only through the power generation side scheduling. There will be insufficient power supply at some periods, so the supply of part of the loads has to be interrupted. By determining the interruptible load in advance to ensure the continuous and stable operation of the microgrid, it is advantageous to ensure the safety and rationality of the load cutting operation.The load model is expressed as follows:(9)PL=PLim+PLin+Ptran,(10)Pnet-load=PL-PPV-PWT,where PL, PLim, PLin, and Ptran are the daily demand load power of users, the important load power, the load power, and the transferable load power, respectively; Pnet-load is the power of the net load.


## 3. Optimization Strategy

Modern power system should pay more attention to the user’s feelings, due to the difference between the optimal objective and the optimal variable; the solution speed, accuracy, and the effect for current optimization of microgrid are not very ideal. In order to consider the comfort of user side comprehensively, the dual-layer scheduling strategy is designed in this paper. The objective of first-layer optimization is to obtain the maximum load satisfaction by considering from the load side, and then the user adjusts the controllable load according to the optimization results; the new load curve will be as the input of the second layer. On the premise of guaranteeing the reliability of the microgrid system, the second-layer optimization dispatches the microgas turbine and battery according to the new load curve and the predictive power of the wind generator, the photovoltaic power, and the objective of the second layer is to optimize the economy, environment-protective, and the utilization of electricity of microgrid. The optimal solution set of Pareto is obtained by adopting the optimization algorithm NSGA-II for multiobjective and introducing method of fuzzy model recognition and then selecting an optimal compromise solution from the Pareto optimal solution set according to the actual situation.The dual-layer scheduling strategy model can better coordinate the interests of different stakeholders in the microgrid and formulate a reasonable optimal operation scheme for microgrid to coordinate the user side and power side. The scheduling strategy of the microgrid system is given in Figure2.Figure 2
The system structure diagram of the dual-layer optimization strategy for microgrid.For electricity load as the management object of user side, it is necessary to adopt different scheduling strategy in the process of management according to the characteristics of load power consumption. Some of the loads will cause a great economic impact and even endanger personal safety after power failure and therefore microgrid system should give priority to ensure the sustainable supply of such loads.The optimal scheduling strategy for each microsource is as follows: the flowchart is given in Figure3:Figure 3
The flowchart of optimal scheduling strategy for microsource.(1) To determine the power relationship of power load and the battery.(2) When the battery is sufficient to supply power to the net load, the microgas turbine is not running.(3) When the battery is not enough to supply power to the net load, the microgas turbine is turned on, and when the microgas turbine and the battery cannot meet the load demand, then cut off the unimportant load.(4) When the load demand is not very large, cut off part of the wind generator.In this paper, the fluctuation of the output of the renewable energy unit can be alleviated by the transfer of electric energy through the cycle charge and discharge of the battery. By placing user side and load on the same layer, the two can be fully coordinated to improve the net load shape. The flowchart of the battery charge-discharge management is shown in Figure4.Figure 4
The flowchart of the battery charge-discharge management.


## 4. Optimization Model Solving Algorithm

### 4.1. The First-Layer Optimal Model of Microgrid System

Adopting certain measures of demand side management in microgrid is conducive to improving the reliability of microgrid operation and the utilization of new energy power generation [21]. In order to improve the economy of microgrid operation from the perspective of load management, this paper introduces the concept of load satisfaction into the energy management of isolated microgrid to satisfy the load demand as much as possible. The original load curve is adjusted and the some periods of unimportant load are transferred by guiding users to use electricity rationally, so as to improve the utilization ratio of renewable energy. The optimization index is as follows: (11)α=PtotalPL,where Ptotal is the total power supply of the system.


### 4.2. The Second-Layer Optimal Model of Microgrid System

#### 4.2.1. Objective Function

(i) Economy. In islanding mode, the objective of microgrid optimal scheduling is to obtain the minimum operation and maintenance cost. The cost of operation and maintenance include four parts: maintenance cost, fuel cost, battery wastage cost, and interruption cost. Maintenance costs are mainly concerned with the wind turbines, photovoltaic cells, and microgas turbines. When the power supply is insufficient, part load must be cut off and the objective function needs to consider the interruption cost.(12)f1=∑i=124K1∗PWT-i+K2∗PPV-i+K3∗PMT-i+∑i=124FMT-i+∑i=124WQlifetimeηrt+εPloadloss-i,where K is the maintenance factor of generating unit; FMT-i is the fuel cost of the microgas turbine at the ith hour; Qlifetime is the output energy of battery life cycle; Ploadloss-i is the power of the resected load, due to insufficient power supply; ε is the economic losses of resected unit load due to insufficient supply of electricity; W is the battery purchase cost.(ii) Pollution Emission (13) f 2 = ∑ i = 1 24 K ′ P G - i C ′ ,where K′ is the pollution emission coefficient; C′ is the cost of sewage disposal; PG-i is output power of generating unit.(iii) The Utilization of Electricity. The highest utilization rate of electric energy is the lowest energy waste rate.(14)f3=∑i=124ΔPi∑i=124PL-i,(15)ΔPi=PPV-i+PWT-i+PMT-i+Pbatt-i-PL-i,where ΔPi>0. The optimization index of the microgrid is as follows:(16)β=∑i=124ΔPi∑i=124PL-i,where ΔPi≤0 and β is the load loss rate.


#### 4.2.2. Optimal Variables

In the island operating mode, the microgrid is not connected with the power grid. The optimal variables in this paper include the charging-discharging power and microgas turbine output power of 24 hours a day.(17)x=Pbat-1,Pbat-2,Pbat-3,…,Pbat-23,Pbat-24,PMT-1,PMT-2,PMT-3,…,PMT-23,PMT-24.


### 4.3. Constraints of the Microgrid System

(i) The Energy Balance Constraints. Microgrid system needs to ensure the power balance between power supply and power consumption. It is defined as (18)PPV-i+PWT-i+Pbatt-i+PMT-i=PL-i,where PPV-i, PWT-i, Pbatt-i, PMT-i, and PL-i are the power of photovoltaic, wind generator, battery, microgas turbine, and the load, respectively.(ii) Generation Capacity Constraints. In order to maintain the stability of the operation, the actual output power of each generator has a strict upper and lower bound:(19)Pmin≤PG-i≤Pmax.(iii) Charge Constraint of Energy Storage Battery Charge-Discharge Power. In the battery work process, with the increase of the depth of the discharge, the loss of the battery increases and the life of the battery is shortened accordingly. So the battery should run in a certain range of charge: (20)SOCmin≤SOC≤SOCmax.(iv) Transferable Load Constraints. The constraints required for transferable load include load power range constraints and load required power invariant constraints:(21)Ptran,min≤Ptran-i≤Ptran,max,(22)∑i=124Ptran-i≤Wtran,where Ptran-i is the power of transferable load; Ptran,min and Ptran,max are the upper and lower limits of power for the allowable transfer load; Wtran is the total amount of electricity needed in the scheduling cycle.


### 4.4. Optimization Algorithm

The dispatching optimization of microgrid is an optimization problem of multiobjective function and multiconstraint conditions, and it often has a conflict between subobjectives. The traditional algorithm is difficult to take into account the relationship between multiple objectives when there are huge differences in target types. Therefore, faced with the multiobjective optimization problem, this paper adopts the NSGA-II intelligent algorithms with elite strategy for solving, to find the Pareto solutions. The NSGA-II algorithm has such advantage as good robustness and can obtain nondominant sets with uniform distribution and good diversity. The multiobjective optimization problem is solved by NSGA-II, which avoids the problem of weight assignment, and it can deal with a variety of multidimensional complex planning problems; the application scope and space of the algorithm are wide. As long as the problem is accurately described in the algorithm, the algorithm can be optimized according to the evolutionary rule, and the optimal solution is obtained, no matter how complex the solution is. And how complicated the problem is, it will not affect the convergence direction and robustness of the algorithm. The strategy of fast nondominating sorting accelerates the speed of the algorithm, while the elitist strategy avoids the loss of the optimal solution and also extends the distribution range of the solution set in the Pareto optimal frontier to ensure the diversity of the population [22–24].The algorithm used in this paper is as follows.Step 1.
Establish the initial population according to the wind turbine, photovoltaic and the predictive power.Step 2 (code).
Encode the optimal variables.The real number coding is used, and the battery powerPbat-i and the microgas turbine power PMT-i are taken as the optimal variables. The chromosome coding diagram is shown in Figure 5, in which the value of P1-P24 is controlled between -21.1973 and 159.727; the value of P25-P48 is controlled between 0 and 2000.Figure 5
Chromosome coding schematic diagram.Step 3.
Set the population size and the number of objective functions and generate the initial population P0 according to the microsource scheduling strategy.Step 4.
Sort out the population based on the nondomination and calculate the crowding distance of the population, in which the nondominated sorting which is based on the three objective function values includes the cost of power system, the pollution emission penalty cost, and the reliability index; the crowding distance is obtained from the distance information of the individual vector variable in the variable space.Step 5 (genetic operation).
Select cross and mutation, and then get the subpopulation.Step 6 (elitist reservation strategy).
Combine the offspring with the parent, select the N optimal individuals according to the nondominance sorting and crowding distance, and create a new population.Step 7.
For the number of iterations plus 1, return to Step4, until the maximum number of iterations is reached.The flowchart of optimization algorithm for the optimization model of microgrid system in this paper is shown as Figure6.Figure 6
The flowchart of optimization algorithm for isolated microgrid system.The above constraint tournament criterion enables the NSGA-II algorithm to solve the interval constraints. NSGA-II uses binary competition selection processing constraints; that is, two solutions in the population are selected and the better solution is selected.The binding tournament criteria used in this paper are defined as follows:(1) When x1 is a feasible solution and x2 is an infeasible solution, then x1 is considered to be the dominant individual.(2) When x1 and x2 are both infeasible solutions, the degree of constraint violation is compared, and the individual with small degree of constraint violation is dominant.(3) When the volumes x1 and x2 are both feasible solutions, the order value and interval crowding distance between them are compared. The individual with small order value is dominant, and the individual with large interval crowding distance is superior if the order value is the same.The optimal solution set of Pareto is obtained by using NSGA-II algorithm. In this paper, a fuzzy model recognition method is used to select a suitable compromise solution which conforms to the actual operation situation from the optimal solution set in Pareto [25]. Firstly, the membership vector of noninferior solution is obtained and then combined with the important degree of each objective in the actual situation, according to the fuzzy tone operator to get the object weight vector. Finally, calculate the proximity of each noninferiority membership degree vector to the object weight vector, select the best closeness degree of membership vector, and the corresponding noninferior solution is the optimal compromise solution to be chosen.The objective membership degree formula of each particle is described as(23)Ai,m=fm,max-fi,mfm,max-fm,min,where Ai,m and fi,m are the membership degree and the actual value of the mth target of the ith particle, respectively; fm,maxand fm,min are the maximum and minimum values of the mth target, respectively.


## 4.1. The First-Layer Optimal Model of Microgrid System

Adopting certain measures of demand side management in microgrid is conducive to improving the reliability of microgrid operation and the utilization of new energy power generation [21]. In order to improve the economy of microgrid operation from the perspective of load management, this paper introduces the concept of load satisfaction into the energy management of isolated microgrid to satisfy the load demand as much as possible. The original load curve is adjusted and the some periods of unimportant load are transferred by guiding users to use electricity rationally, so as to improve the utilization ratio of renewable energy. The optimization index is as follows: (11)α=PtotalPL,where Ptotal is the total power supply of the system.


## 4.2. The Second-Layer Optimal Model of Microgrid System

### 4.2.1. Objective Function

(i) Economy. In islanding mode, the objective of microgrid optimal scheduling is to obtain the minimum operation and maintenance cost. The cost of operation and maintenance include four parts: maintenance cost, fuel cost, battery wastage cost, and interruption cost. Maintenance costs are mainly concerned with the wind turbines, photovoltaic cells, and microgas turbines. When the power supply is insufficient, part load must be cut off and the objective function needs to consider the interruption cost.(12)f1=∑i=124K1∗PWT-i+K2∗PPV-i+K3∗PMT-i+∑i=124FMT-i+∑i=124WQlifetimeηrt+εPloadloss-i,where K is the maintenance factor of generating unit; FMT-i is the fuel cost of the microgas turbine at the ith hour; Qlifetime is the output energy of battery life cycle; Ploadloss-i is the power of the resected load, due to insufficient power supply; ε is the economic losses of resected unit load due to insufficient supply of electricity; W is the battery purchase cost.(ii) Pollution Emission (13) f 2 = ∑ i = 1 24 K ′ P G - i C ′ ,where K′ is the pollution emission coefficient; C′ is the cost of sewage disposal; PG-i is output power of generating unit.(iii) The Utilization of Electricity. The highest utilization rate of electric energy is the lowest energy waste rate.(14)f3=∑i=124ΔPi∑i=124PL-i,(15)ΔPi=PPV-i+PWT-i+PMT-i+Pbatt-i-PL-i,where ΔPi>0. The optimization index of the microgrid is as follows:(16)β=∑i=124ΔPi∑i=124PL-i,where ΔPi≤0 and β is the load loss rate.


### 4.2.2. Optimal Variables

In the island operating mode, the microgrid is not connected with the power grid. The optimal variables in this paper include the charging-discharging power and microgas turbine output power of 24 hours a day.(17)x=Pbat-1,Pbat-2,Pbat-3,…,Pbat-23,Pbat-24,PMT-1,PMT-2,PMT-3,…,PMT-23,PMT-24.


## 4.2.1. Objective Function

(i) Economy. In islanding mode, the objective of microgrid optimal scheduling is to obtain the minimum operation and maintenance cost. The cost of operation and maintenance include four parts: maintenance cost, fuel cost, battery wastage cost, and interruption cost. Maintenance costs are mainly concerned with the wind turbines, photovoltaic cells, and microgas turbines. When the power supply is insufficient, part load must be cut off and the objective function needs to consider the interruption cost.(12)f1=∑i=124K1∗PWT-i+K2∗PPV-i+K3∗PMT-i+∑i=124FMT-i+∑i=124WQlifetimeηrt+εPloadloss-i,where K is the maintenance factor of generating unit; FMT-i is the fuel cost of the microgas turbine at the ith hour; Qlifetime is the output energy of battery life cycle; Ploadloss-i is the power of the resected load, due to insufficient power supply; ε is the economic losses of resected unit load due to insufficient supply of electricity; W is the battery purchase cost.(ii) Pollution Emission (13) f 2 = ∑ i = 1 24 K ′ P G - i C ′ ,where K′ is the pollution emission coefficient; C′ is the cost of sewage disposal; PG-i is output power of generating unit.(iii) The Utilization of Electricity. The highest utilization rate of electric energy is the lowest energy waste rate.(14)f3=∑i=124ΔPi∑i=124PL-i,(15)ΔPi=PPV-i+PWT-i+PMT-i+Pbatt-i-PL-i,where ΔPi>0. The optimization index of the microgrid is as follows:(16)β=∑i=124ΔPi∑i=124PL-i,where ΔPi≤0 and β is the load loss rate.


## 4.2.2. Optimal Variables

In the island operating mode, the microgrid is not connected with the power grid. The optimal variables in this paper include the charging-discharging power and microgas turbine output power of 24 hours a day.(17)x=Pbat-1,Pbat-2,Pbat-3,…,Pbat-23,Pbat-24,PMT-1,PMT-2,PMT-3,…,PMT-23,PMT-24.


## 4.3. Constraints of the Microgrid System

(i) The Energy Balance Constraints. Microgrid system needs to ensure the power balance between power supply and power consumption. It is defined as (18)PPV-i+PWT-i+Pbatt-i+PMT-i=PL-i,where PPV-i, PWT-i, Pbatt-i, PMT-i, and PL-i are the power of photovoltaic, wind generator, battery, microgas turbine, and the load, respectively.(ii) Generation Capacity Constraints. In order to maintain the stability of the operation, the actual output power of each generator has a strict upper and lower bound:(19)Pmin≤PG-i≤Pmax.(iii) Charge Constraint of Energy Storage Battery Charge-Discharge Power. In the battery work process, with the increase of the depth of the discharge, the loss of the battery increases and the life of the battery is shortened accordingly. So the battery should run in a certain range of charge: (20)SOCmin≤SOC≤SOCmax.(iv) Transferable Load Constraints. The constraints required for transferable load include load power range constraints and load required power invariant constraints:(21)Ptran,min≤Ptran-i≤Ptran,max,(22)∑i=124Ptran-i≤Wtran,where Ptran-i is the power of transferable load; Ptran,min and Ptran,max are the upper and lower limits of power for the allowable transfer load; Wtran is the total amount of electricity needed in the scheduling cycle.


## 4.4. Optimization Algorithm

The dispatching optimization of microgrid is an optimization problem of multiobjective function and multiconstraint conditions, and it often has a conflict between subobjectives. The traditional algorithm is difficult to take into account the relationship between multiple objectives when there are huge differences in target types. Therefore, faced with the multiobjective optimization problem, this paper adopts the NSGA-II intelligent algorithms with elite strategy for solving, to find the Pareto solutions. The NSGA-II algorithm has such advantage as good robustness and can obtain nondominant sets with uniform distribution and good diversity. The multiobjective optimization problem is solved by NSGA-II, which avoids the problem of weight assignment, and it can deal with a variety of multidimensional complex planning problems; the application scope and space of the algorithm are wide. As long as the problem is accurately described in the algorithm, the algorithm can be optimized according to the evolutionary rule, and the optimal solution is obtained, no matter how complex the solution is. And how complicated the problem is, it will not affect the convergence direction and robustness of the algorithm. The strategy of fast nondominating sorting accelerates the speed of the algorithm, while the elitist strategy avoids the loss of the optimal solution and also extends the distribution range of the solution set in the Pareto optimal frontier to ensure the diversity of the population [22–24].The algorithm used in this paper is as follows.Step 1.
Establish the initial population according to the wind turbine, photovoltaic and the predictive power.Step 2 (code).
Encode the optimal variables.The real number coding is used, and the battery powerPbat-i and the microgas turbine power PMT-i are taken as the optimal variables. The chromosome coding diagram is shown in Figure 5, in which the value of P1-P24 is controlled between -21.1973 and 159.727; the value of P25-P48 is controlled between 0 and 2000.Figure 5
Chromosome coding schematic diagram.Step 3.
Set the population size and the number of objective functions and generate the initial population P0 according to the microsource scheduling strategy.Step 4.
Sort out the population based on the nondomination and calculate the crowding distance of the population, in which the nondominated sorting which is based on the three objective function values includes the cost of power system, the pollution emission penalty cost, and the reliability index; the crowding distance is obtained from the distance information of the individual vector variable in the variable space.Step 5 (genetic operation).
Select cross and mutation, and then get the subpopulation.Step 6 (elitist reservation strategy).
Combine the offspring with the parent, select the N optimal individuals according to the nondominance sorting and crowding distance, and create a new population.Step 7.
For the number of iterations plus 1, return to Step4, until the maximum number of iterations is reached.The flowchart of optimization algorithm for the optimization model of microgrid system in this paper is shown as Figure6.Figure 6
The flowchart of optimization algorithm for isolated microgrid system.The above constraint tournament criterion enables the NSGA-II algorithm to solve the interval constraints. NSGA-II uses binary competition selection processing constraints; that is, two solutions in the population are selected and the better solution is selected.The binding tournament criteria used in this paper are defined as follows:(1) When x1 is a feasible solution and x2 is an infeasible solution, then x1 is considered to be the dominant individual.(2) When x1 and x2 are both infeasible solutions, the degree of constraint violation is compared, and the individual with small degree of constraint violation is dominant.(3) When the volumes x1 and x2 are both feasible solutions, the order value and interval crowding distance between them are compared. The individual with small order value is dominant, and the individual with large interval crowding distance is superior if the order value is the same.The optimal solution set of Pareto is obtained by using NSGA-II algorithm. In this paper, a fuzzy model recognition method is used to select a suitable compromise solution which conforms to the actual operation situation from the optimal solution set in Pareto [25]. Firstly, the membership vector of noninferior solution is obtained and then combined with the important degree of each objective in the actual situation, according to the fuzzy tone operator to get the object weight vector. Finally, calculate the proximity of each noninferiority membership degree vector to the object weight vector, select the best closeness degree of membership vector, and the corresponding noninferior solution is the optimal compromise solution to be chosen.The objective membership degree formula of each particle is described as(23)Ai,m=fm,max-fi,mfm,max-fm,min,where Ai,m and fi,m are the membership degree and the actual value of the mth target of the ith particle, respectively; fm,maxand fm,min are the maximum and minimum values of the mth target, respectively.


## 5. Simulation

This study takes a microgrid system including wind turbine, photovoltaic, microgas turbine, and battery as an example. The installed capacity is 1000kW photovoltaic cell, 1650kW wind generator, 530kW lithium battery, and 10 × 250kW microgas turbine. The parameters of each distributed energy resource (DER) are shown in Table1. The emission parameters of the unit are shown in Table 2, and the simulation parameters of the battery are shown in Table 3. On the Matlab 2014 platform, NSGA-II is applied to solve the model of dual-layer optimal dispatching strategy for microgrid energy management system. The population size is 50, and the cycle number is 200. Genetic manipulation parameters: the crossover rate is 0.8, and the variation rate is 0.05. The curve of radiation intensity and the curve of wind speed at the wheel hub of wind turbine are shown in Figures 7 and 8, respectively. The three-dimensional Pareto solution set as shown in Figure 9 is obtained through simulation; the number of nondominated solutions is 50. The load curve before and after demand response as shown in Figure 10 shows that the load transfer is an important management measure of load, which has fast response speed and can effectively mobilize the interactive enthusiasm of the user to participate in the optimal dispatching for microgrid. The participation of users can improve the economy of microgrid operation effectively.Table 1
DER parameters.
Type lower limit of power (kW) Upper limit of power(kW) Maintenance coefficient ($/kW.h) WT 0 1650 0.00446 PV 0 1000 0.00145 MT 0 2000 0.00604 BAT -21 159.8 0.00414Table 2
Pollution conversion costs and emission factors.
Type Emission coefficient(kg/kW.h) Emission cost($/kg) CO2 0.0016 0.0131 SO2 0.000008 0.9487 NOx 0.00044 3.9904Table 3
Simulation parameters of storage battery.
Name V n η b a t t , n α c k c I m a x Q m a x numerical value 4 V 0.80 1 0.5281 hr−1 0.254 67.5 A 1887 AhFigure 7
The curve of radiation intensity.Figure 8
The curve of wind speed at the wheel hub of wind turbine.Figure 9
The Pareto frontier of optimal dispatching results.Figure 10
The load curve of demand response.Figures11 and 12 are Pareto front diagram for each two objective functions of the three objective function, and the output of each microsource is shown in Figures 13 and 14. Table 4 lists three groups of extreme solutions in the Pareto solution sets, including the minimum operating cost, the least pollution emission, and the highest reliability, and it also lists the optimal compromise solution by using the fuzzy model recognition method. The analysis of the simulation results show that the relationship between the various goals is mutually binding. The optimality of a goal must be at the expense of the optimization of the other two goals. For instance, to achieve the optimal economy, the corresponding operation and maintenance cost is $265.89, but the pollution emissions and energy waste rate reached 49.76kg and 0.532. It can also be seen from Table 4 that the generation cost and pollution emission in the day after the optimal choice of compromise solution is slightly higher than the minimum generation cost and the lowest environmental cost, while the environmental costs that are aiming at the lowest environmental costs and the lowest ones are lower than the minimum and the lowest environmental costs. The choice of the optimal compromise solution is well coordinated economic, environmental, and energy wastage among the three balances.Table 4
Typical solution of Pareto solution.
Type Maintenance cost($) Pollution emission (kg) Power waste rate Minimum maintenance cost 265.89 49.76 0.532 Minimum pollution emission 1804.90 17.04 0.216 Minimum energy waste rate 1781.50 17.28 0.214 Optimal compromise solution 417.07 35.35 0.335Figure 11
Pareto diagram of f1 and f2.Figure 12
Pareto diagram of f1 and f3.Figure 13
Output power of renewable energy.Figure 14
Optimal dispatching results of microsource.For further processing of the data, Table5 gives the comparison of the indicators before and after the demand response. According to Table 5 the operation and maintenance cost of microgrid is $724.91 before the application of dual-layer optimal dispatching strategy for microgrid energy management system based on demand side response, which is much higher than the operation and maintenance cost after demand response. At the same time, using the optimal strategy proposed in this paper, the environment-protective and reliability of the microgrid are all improved, and the optimization effect is good. The dual-layer optimal strategy proposed in this paper takes into account the benefits of microgrid operators and users. The power supply system of microgrid has better power supply qualities and takes account of environmental protection and renewable energy utilization. It realizes the economic dispatch and rational utilization of distributed generation.Table 5
Each index value of microgrid.
Economic indicators Environmental indicators reliability index maintenance costs ($) Renewable energy output ratio Environmental cost ($) Loss rate of load Power waste rate Before the response 724.91 0.761 10.92 0.061 0.320 After the response 417.07 0.803 10.25 0.031 0.335In addition, the optimal allocation of microgrid energy management system in different seasons is shown in Figure15. Take typical days as an example; run the program twenty times to get statistical analysis results as shown in Table 6. The NSGA-II has good robustness and convergence performance in dealing with large-scale nonlinear complex programming problems. The study shows that the NSGA-II is of great significance for improving the reliability of power supply, reducing the cost, and reducing the energy waste. In summary, the multiobjective optimization model proposed in this paper can accurately reflect the actual operation of microgrid. The model can achieve better environmental protection effect and good renewable energy utilization rate with small economic cost as small as possible and realize energy saving and environmental protection power generation scheduling.Table 6
Statistical results obtained by NSGA-II for microgrid system of typical day.
Best Worst Average Std. Operation and maintenance cost($) 403.05 425.10 411.24 7.54e-00 Pollution emission(kg) 31.89 37.46 35.42 1.85e-00 Power waste rate 0.29 0.37 0.33 2.48e-02Figure 15
Optimal allocation of energy system for microgrids in different seasons.
(a)
Load demand curve (b)
Spring and autumn (c)
Summer (d)
Winter


## 6. Conclusions

In this paper, a microsource model of wind power, photovoltaic, microgas turbines, and energy storage devices is firstly established, and then the dual-layer optimization model of microgrid energy considering demand response is set up. With the objective of minimizing the operation cost of microgrid and the pollution emission and maximizing the load satisfaction and the reliability, the NSGA-II is used to obtain the optimal solution of Pareto. The methods and steps of the solution are discussed in detail. Finally, the fuzzy model recognition method is adopted to make the decision of the optimal compromise solution.The simulation results show that the dual-layer optimal dispatching model of microgrid energy considering the load satisfaction can take full account of the users. On the one hand, it solves the problem of multiple optimal variables in the optimal operation of microgrid. On the other hand, it is beneficial to coordinate the user side and the power side and greatly reduce the operation loss of the system, thereby reducing the operation cost and improving the reliability of the isolated microgrid. The dual-layer optimal dispatching model is suitable for coordinating the interests of different individuals with different objectives and different decision variables, fully exploiting the potential of user side to participate in the optimization of microgrid operation, and achieving global optimization. The system reliability index measures the power supply quality and ensures the power supply reliability, and the energy management and control strategy of microgrid for the isolated microgrid can make full use of microgrid power capacity, coordinate the start, stop and output between distributed power sources, and realize the economic operation of the system.In future work, the optimal allocation of the microgrid energy in different seasons, weather, and regions can be considered.



---

*Source: 2695025-2018-07-08.xml*