# Dynamic Operation Management of a Renewable Microgrid including Battery Energy Storage

Authors: Shaozhen Jin; Zhizhong Mao; Hongru Li; Wenhai Qi
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In this paper, a novel dynamic programming technique is presented for optimal operation of a typical renewable microgrid including battery energy storage. The main idea is to use the scenarios analysis technique to proceed the uncertainties related to the available output power of wind and photovoltaic units and dynamic programming technique to obtain the optimal control strategy for a renewable microgrid system in a finite time period. First, to properly model the system, a mathematical model including power losses of the renewable microgrid is established, where the uncertainties due to the fluctuating generation from renewable energy sources are considered. Next, considering the dynamic power constraints of the battery, a new performance index function is established, where the Lagrange multipliers and interior point method will be presented for the equality and inequality operation constraints. Then, a feedback control scheme based on the dynamic programming is proposed to solve the model and obtain the optimal solution. Finally, simulation and comparison results are given to illustrate the performance of the presented method.


---

## Body

## 1. Introduction

Nowadays, renewable energy sources (RESs) such as wind or photovoltaics have become more wide spread due to needs for satisfying the environment concerns. On the other hand, distributed generators (DGs) like diesel engines, microturbines, and fuel cells can be used to enhance the resiliency of power system and yield other social economic benefits. Therefore, renewable microgrid is expected to play an important role in future power systems [1, 2]. As a key enabling element of renewable microgrids, battery energy storages make microgrid become a strong coupling system in the time domain. In this regard, the methodologies applied to operation management of a renewable microgrid are getting more complicated and challengeable; therefore there is a strong need for more reliable scheduling of energy sources in renewable microgrid including battery energy storage.So far, many researchers have dealt with the optimal operation scheduling of energy sources in microgrids [3–7]. Previously conventional mathematical programming such as Lagrange relaxation [8, 9], lambda iteration [10, 11], Newton-Raphson [12], interior point method [13], weighted minimax [14], and quadratic programming [15] have been used to determine the least cost solution. However, the conventional mathematical programming methods have major disadvantages such that they can be trapped in local optimal, exhibited sensitivity to the initial starting points. And many of the methods cannot solve the nonsmooth, convex, and nonmonotonically increasing cost functions. Recently, computational intelligence [16, 17] and artificial intelligence based nonconventional methods [18] have been used to solve the optimal operation scheduling of energy sources in microgrids. Artificial intelligence based methods such as artificial neural network and computational intelligence methods such as genetic algorithm, particle swarm optimization, harmony search, simulated annealing, differential evolution, gravitational search algorithm, biogeography based optimization, bacterial foraging algorithm, ant colony optimization, cuckoo search, bat algorithm, artificial bee colony, firefly algorithm, and flower pollination algorithm have been used to solve the problem. These methods can enable us to solve the nonlinear and no-convex cost functions and can obtain nearly the global solutions. However, these methods have major disadvantages such as the evolutionary algorithms greatly depending on their parameters and having high computational time. Besides, hybrid methods which combine different algorithms have been used to solve the optimal operation scheduling of energy sources in microgrid. But, these methods usually have long computational time. Also, all previous researchers search the optimal solutions in each separated time periods. In other words, they are static optimal operation scheduling which do not consider the relationships among different time periods. Dynamic optimal operation scheduling problem is another fundamental part of the renewable microgrid operation to maintain the power balance [19–24]. Among these methods, Cheng et al. [23] have used an enhanced quorum sensing based particle swarm optimization to deal with the dynamic operation optimization problem. However, it is hard to solve the complex optimization problem with high dimension variables and multiple constraints by such metaheuristic methods. Liu et al. [19, 22] have used the dynamic programming algorithm for handing the dynamic economic dispatch problem. Generally, dynamic programming can be used to solve nonlinear dynamics when the problems can be discredited with time, state, and decision variable [25]. In other words, dynamic programming can be used to solve the problem by dividing the decision process into different stages. Meanwhile, these stages interact and interconnect with decision variables. And the target of the dynamic programming is used to find the decision set over the whole searching trajectory. However, such methods may have some deficiencies handling large scale problems. And a challenging aspect of this method is determining which of the inequality constraints are binding at the solution. Another relevant aspect in renewable microgrid operation management should be coping with uncertainty in the renewable energy sources, load demand, and market prices [24–26]. Similarly, power line losses should also be considered in the renewable microgrid system. There have been some researches including the power losses in a renewable microgrid [5, 27, 28]. For example, the output power losses of distributed generators are modeled in articles [5, 29]. The power line losses which are calculated using B-coefficients are considered in the economic dispatch problems [30, 31].This paper presents several contributions as follows:( 1 ) The power losses of the renewable microgrid are included. It is noteworthy that, from an optimal operation planning point of view, such method can be optimized a renewable microgrid globally.( 2 ) Uncertainty in the renewable energy resources is considered. A methodology based on scenario analysis is used to proceed the uncertainties related to the available output power of wind and photovoltaic units.( 3 ) The dynamic performance of the battery storage system is considered. Therefore, this dynamic operation constraint makes renewable microgrid to become a strong coupling system in the time domain. Then, the Lagrange multipliers and interior point method will be presented for these operation constraints.( 4 ) A feedback control scheme based on the dynamic programming is proposed to solve the model and obtain the optimal solution. The advantage of the scheme is that the number of numerical operations is linear in the number of optimization parameters, which enables one to solve for a smooth input/control history. Meantime, this scheme enables one to initialize the optimization problem with a zero-initial guess. Finally, this scheme is a direct method to obtain the same results.


## 2. Description of the Renewable Microgrid

### 2.1. Renewable Microgrid

In this paper, the renewable microgrid system considers diesel generators, microturbines, fuel cells, wind turbines, photovoltaics, and a battery energy storage system in Figure1. The renewable microgrid can operate in parallel to the main grid or as an island.Figure 1
The smart microgrid structure diagram.


### 2.2. Active Power Balance of the Renewable Microgrid

In general, the active power balance between electrical energy production and consumption must be met at each time interval expressed as follows [28, 33]:(1)∑i=1Nwpwt+∑i=1Nspst+∑i=1Ndpdt+∑i=1Nmpmt+∑i=1Nfpft+∑i=1Nbpbt+pgt,=∑j=1NDpLt+Lpdgt,t=1,…,T,where pw(t),ps(t),pd(t),pm(t),pf(t),pb(t),pg(t),pdg(t), pL(t) are the power set-point for the wind turbines, photovoltaics, diesel generators, microturbines, fuel cells, battery energy storage, power exchange between the main grid and renewable microgrid, controllable distribution generators, and load levels. L(pdg(t)) is the power loss of the renewable microgrid. Nw,Ns,Nd,Nm,Nf,Nb,ND are the total number of wind turbines, photovoltaics, diesel generators, microturbines, fuel cells, battery energy storage, and load levels. In this paper, the controllable distribution generators (DGs) included diesel generators, microturbines, and fuel cells.


#### 2.2.1. Estimation of Power Losses Coefficients

When the power losses of the renewable microgrid are considered, these power losses coefficients can be calculated as follows [28, 34]:(2)Lpdgt=∑m=1Ndg∑n=1Ndgpdg,mt.Bmnpdg,nt+∑n=1NdgBm0.pdg,m+B00,where B=[Bmn] and Bm0 and B00 are the power losses coefficients.


### 2.3. Subsystem Model of the Renewable Microgrid

This paper considers a renewable microgrid containing several types of devices, such as battery energy storages, distributed generators, wind turbines, and photovoltaics. A block diagram is given in Figure1.


#### 2.3.1. The Model of Battery Energy Storage System

Battery energy storage can store energy when charging. On the other hand, it can supply energy to the load demand when discharging.The degradation cost of battery energy storage can be calculated as follows [33]:(3)Cb,ct=∑j=1Nb∑i=1trb.pbi,where rb is the cost coefficient which can be determined by battery energy storage charging or discharging power.p b , c ( t ) is the permitted rate of charge during a definite period of time and pb,d(t) is the permitted rate of discharge. The following equation can be expressed for a battery energy storage(4)pbt=uct.pb,ct+udt.pb,dt,where uc(t) or ud(t) indicate the working mode of the battery energy storage, uc(t),ud(t)∈0,1. Define uc(t)=1 or ud(t)=1 if the battery energy storage is charge or discharge. These variables satisfy(5)uct+udt≤1,t=1,…,T.The permitted rate of charge and discharge are limited by the following constrains:(6)0≤pb,ct≤pb,cmax,0≤pb,dt≤pb,dmax,where pb,cmax and pb,dmax are the maximum charging and discharging rates of battery energy storage.After charge or discharge, the state of charge of the battery energy storage can be calculated as follows [33]:(7)SOCbt=SOCbt-1+ηb,c.pb,ct-pb,dt/ηb,dQb,c.Δt,where ηb,c,ηb,d are the energy conversion coefficients. Qb,c is the energy capacity of the battery energy storage. On the other hand, the state of charge of the battery energy storage must be maintained in the following range. Δt is the operating time of the battery energy storage.(8)SOCb,min≤SOCbt≤SOCb,max,where SOCb,min,SOCb,max are required limits of the state of charge.


#### 2.3.2. The Model of Controllable Distributed Generators

In this paper, the distribution generators (DGs) are included diesel generators, microturbines, and fuel cells.(A) Diesel Generators. In general, diesel generators serve as a backup power source. Meanwhile, the operating cost of diesel generators is formulated as follows:(9)Cdt=∑j=1Nd∑i=1tCd,fi+Cd,omi+Rdi.1-δdi-1.δdi,where Cd,f is the fuel cost of diesel generator; Cd,om is the operating and maintenance cost; Rd is the start-up cost of diesel generator.The fuel cost of diesel generators is modeled as a function of their actual output power(10)Cd,ft=a.pdt2+b.pdt+c,where pd(t) is the actual output power. a,b,c are the fuel consumption curve fitting coefficients.The operating and maintenance cost can be expressed as(11)Cd,omt=kd,om.pdt.Δt,where kd,om is the operating and maintenance cost coefficient. Δt is the operating time of diesel generator.The start-up cost of diesel generator is relative to its operating state(12)Rdt=Rd,h,Td,off≤Td,d+Td,c,Rd,c,Td,off>Td,d+Td,c,where Rd,h,Rd,c are the hot start-up cost and cold start-up cost. Td,c is the cooling time of diesel generator. Td,off,Td,d are the off time and shut down time of diesel generator.For a stable operation, the actual output power of diesel generator is limited by lower and upper bounds as follows:(13)pd,min≤pdt≤pd,max,where pd,min is the minimum active power of diesel generator. pd,max is the maximum power of diesel generator.The following constraint ensures diesel generators do not exceed their ramp limits:(14)-rd.pd,max≤pdt-pdt-1≤rd.pd,max,where rd is the ramp constraint coefficient.Due to the physical constraints that state shifting can only take place after a fixed time interval, the state variable of two adjacent times is as follows:(15)Td,ont-1-Td,u.δdt-1-δdt≥0,(16)Td,offt-1-Td,d.δdt-δdt-1≥0,where Td,on,Td,u are the cumulative uptime and the minimum turn on time.(B) Fuel Cells. Fuel cell directly converts chemical energy to electrical energy by electrochemical reactions, which is one of the most promising energy conversion technologies. The operating cost of fuel cells are formulated as(17)Cft=∑j=1Nf∑i=1tCf,fi+Cf,omi+Rfi.1-δfi-1.δfi,where Cf,f is the fuel cost of fuel cell; Cf,om is the operating and maintenance cost; Rf is the start-up cost of fuel cell.The fuel cost of fuel cells is modeled as(18)Cf,ft=cf.Tf.∑i=1tpfi+paηfi,where cf is the fuel price; Tf is the operating time of changing two adjacent states; pa is the consumption with active power for the auxiliary equipment. ηf is the efficiency of fuel cells.The operating and maintenance cost can be obtained by(19)Cf,omt=kf,om.Tf.∑i=1tpfi,where kf,om is the operating and maintenance cost coefficient.The start-up cost of fuel cells can be expressed as(20)Rft=Rf,h+βf.1-e-Tfofft/τf,where Rf,h is the hot start-up cost of fuel cells; Rf,h+βf is the cold start-up cost of fuel cells. Tfoff is the off time of fuel cells; τf is the cooling time of fuel cells.The actual output power of fuel cell is limited by lower and upper bounds as(21)pf,min≤pft≤pf,max,where pf,min,pf,max are the minimum and maximum active power of fuel cells.The following constraint ensures fuel cell does not exceed their ramp limits(22)-Δpf,D≤pft-pft-1≤Δpf,U,where Δpf,D,Δpf,U are the minimum and maximum ramp rate of fuel cells.The minimum on or minimum down time constraints for fuel cells can be expressed as(23)Tf,ont-1-Tf,u.δft-1-δft≥0,(24)Tf,offt-1-Tf,d.δft-δft-1≥0,where Tf,on is the cumulative uptime. Tf,u,Tf,d are the minimum turn on time and shut down time of fuel cells.The number of turn on and turn off can be expressed as(25)nf,st-spt≤Nf,max,where nf,st-sp(t) is the number of turn on and turn off in time t; Nf,max is the maximum number of turn on and turn off.(C) Microturbine. Microturbines have the higher power density, produce less noise, and emit much less pollutants, especially NOx. So microturbines are effective devices of converting the fuel energy into electrical energy. Meanwhile, the operating cost of microturbines are formulated as follows:(26)Cmt=∑j=1Nm∑i=1tCm,fi+Cm,omi+Rmi.1-δmi-1.δmi,where Cm,f is the fuel cost of microturbine; Cm,om is the operating and maintenance cost; Rm is the start-up cost of microturbine.The fuel cost of microturbines is modeled as(27)Cm,ft=cmpmtηmt,where cm is the fuel price of microturbines; ηm is the efficiency of microturbines.The operating and maintenance cost of microturbines can be expressed as(28)Cm,omt=km,om.pmt.Δt,where km,om is the operating and maintenance cost coefficient of microturbines.The start-up cost of microturbines can be expressed as(29)Rmt=Rm,h+βm.1-e-Tmofft/τm,where Rm,h is the hot start-up cost of microturbines; Rm,h+βm is the cold start-up cost of microturbines. Tmoff is the off time of microturbines; τm is the cooling time of microturbines.The actual output power of microturbines is limited by lower and upper bounds as follows:(30)pm,min≤pmt≤pm,max,where pm,min,pm,max are the minimum and maximum active power of microturbines.The following constraint ensures microturbines do not exceed their ramp limits(31)-Δpm,D≤pmt-pmt-1≤Δpm,U,where Δpm,D,Δpm,U are the minimum and maximum ramp rate of microturbines.The minimum on or minimum down time constraints for microturbines can be expressed as(32)Tm,ont-1-Tm,u.δmt-1-δmt≥0,(33)Tm,offt-1-Tm,d.δmt-δmt-1≥0,where Tm,on is the cumulative uptime; Tm,u,Tm,d are the minimum turn on time and shut down time of microturbines.


#### 2.3.3. The Probability Model of Wind Turbines

Wind powers as renewable energy sources (RESs) are dependent on numerous factors such as wind velocity and efficiency. The approximate relationship between the wind power and wind speed can be expressed as(34)pwv=0,v<vi,v≥v0,v-vi.prvr-vi,vi≤v<vr,pr,vr≤v≤v0,where vi,vo,vr are the cut-in wind speed, cut-out wind speed, and rated wind speed. pr is the rated output power of wind turbines.Prior research has shown that the wind speed profile follows the Weibull distribution over time, which is [29, 32](35)Fvv=1-exp⁡-vck,where positive variables c and k are the scale parameter and shape parameter, respectively.Based on the characteristic of the wind power (34) and the probability distribution function (35), the probability distribution function of wind power is obtained by(36)FPwv=0,pwv<0,1-exp⁡-vi+vr-vipwv/prc,+exp⁡-v0ck,0≤pwv≤pr,1,pwv≥pr.


#### 2.3.4. The Probability Model of Photovoltaics

The output power of photovoltaic can be calculated by its rated output power at the standard test condition and the operating ambient temperature [35](37)ps=Ac.η.Iβ=Ac.η.T.kt-T′.kt2,where Ac is the array surface area; η is the efficiency of the photovoltaic in realistic reporting conditions; Iβ is the irradiance on a surface with inclination β to the horizontal plane; T,T′ are the parameters that depend on inclination β, reflectance of the ground ρ, etc.(38)T=Rb+ρ.1-cos⁡β2+1+cos⁡β2-Rb.p.rd.H03600,(39)T′=1+cos⁡β2-Rb.q.rd.H03600,where Rb is the ratio of beam radiation on the tilted surface to that on a horizontal surface at any time and H0 is the solar radiation, for that day, both referring to a horizontal surface; rd is ratio, diffuse radiation in hour or diffuse in day; p,q are the parameters of the daily clearness index kt.Many researches have proved that cloudiness is the main factor affecting the difference between the values of solar radiation measured outside the atmosphere and on earthly surface. The daily clearness index can be obtained by [36](40)kt=ItI0,where It,I0 are the ratio of the irradiance on horizontal plane and the extraterrestrial total solar irradiance.The effect of clouds on terrestrial irradiance is the daily clearness index;kt can not be predicted with complete confidence; it must be treated as random variable(41)fkt,k¯t=Cktu.λ.γ1eλ.kt1-γ1.kt-1,(42)k¯t=νs.ktu-1.exp⁡νs.ktu-νs.ktl-1.exp⁡νs.ktlνs.exp⁡νs.ktu-exp⁡νs.ktl,where ktl,ktu are the lower and upper bounds of the observed range for kt; νs is the parameter of daily clearness index.(43)C=λ2.ktueλ.ktu-1-λ.ktu,(44)γ1=λ1+λ.ktu,(45)λ=2Γ-17.519.exp⁡-1.3118.Γ-1062.exp⁡-5.0426.Γktu,(46)Γ=ktuktu-k¯t,where C,γ1,λ are the parameters of daily clearness index.In particular, ifT>0,T′<0, the probability density function fps can be expressed as(47)fps=C.ku-α+α′/2-ktu.Ac.η.T′.α′.eλ.α+α′/2,ps∈0,psktu,0,ps∉0,psktu,while if T>0,T′>0, the probability density function can be obtained by(48)fps=C.ku-α-α′/2-ktu.Ac.η.T′.α′.eλ.α-α′/2,ps∈0,psktu,0,ps∉0,psktu,where α=T/T′,α′=α2-4.ps/η.T′.Ac.


#### 2.3.5. The Model of Interaction with External Main Grid

In this paper, the renewable microgrid is connected to the external main grid and can trade energy with the main grid. The transaction incurs the following cost to the renewable microgrid:(49)Cgt=∑i=1tρg.pgi.Let the amount of energy be bounded by(50)-pg,max≤pgt≤pg,max(51)ρg=ρs,pg≤0.ρb,pg>0.where pg,max is the maximum transaction limits; ρs,ρb are the price coefficients to purchase and sell energy.


### 2.4. Dynamic Operation Model of the Renewable Microgrid

Dynamic economic operation management of the renewable microgrid is to determine output power of distribution generators, in order to minimize total operation cost of the renewable microgrid and meet the dynamic operation constrains.The total operating cost of the renewable microgrid can be defined as(52)min⁡Ft=∑t=1TCb,ct+Cdt+Cft+Cmt+Cgt=minΠ=u1,…,uT⁡E∑t=1TCtxt,ut,wt,where u(t) is the decision or control variables (specifically, action at state x(t) in period t); x(t) is the state variables of the renewable microgrid; w(t) is the random variable such as the output power of wind turbine or photovoltaics.u ( t ) , x ( t ) , w ( t ) can be defined as follows:(53)ut=Δpbt,Δpdt,Δpmt,Δpft,Δpgt,(54)xt=SOCbt,pbt,pdt,pmt,pft,pgt,(55)wt=pwt,pst.The state variables of renewable microgrid int+1 period can be defined as follows:(56)xt+1=fxt,ut,wt=SOCbt+1,pbt+Δpbt,pdt+Δpdt,pmt+Δpmt,pft+Δpft,pgt+Δpgt.The operation constraints of the renewable microgrid can be expressed as power balance (1)-(2), battery energy storage limits (5)-(8), distributed generations limits (13)-(16), (21)-(25), (30)-(33), and interaction with external main grid limits (50).


#### 2.4.1. The Standard Formulation of the Dynamic Operation Model

According to dynamic programming formulation, the dynamic operation model of the renewable microgrid can be formulated as(57)Jtxt=minut∈Ut⁡Ctxt,ut,wt+EJt+1xt+1.The operation constraints can be expressed as(58)xt+1=fxt,ut,wt,(59)ctExt,ut,wt=0,cTExT=0,(60)ctIxt,ut,wt≤0,cTIxT≤0,where ctE,cTE,ctI,cTI are, respectively, the equality and inequality constraints.Meanwhile, the constraint of condition should be satisfied as follows:(61)JT+1xT+1≡0.


## 2.1. Renewable Microgrid

In this paper, the renewable microgrid system considers diesel generators, microturbines, fuel cells, wind turbines, photovoltaics, and a battery energy storage system in Figure1. The renewable microgrid can operate in parallel to the main grid or as an island.Figure 1
The smart microgrid structure diagram.


## 2.2. Active Power Balance of the Renewable Microgrid

In general, the active power balance between electrical energy production and consumption must be met at each time interval expressed as follows [28, 33]:(1)∑i=1Nwpwt+∑i=1Nspst+∑i=1Ndpdt+∑i=1Nmpmt+∑i=1Nfpft+∑i=1Nbpbt+pgt,=∑j=1NDpLt+Lpdgt,t=1,…,T,where pw(t),ps(t),pd(t),pm(t),pf(t),pb(t),pg(t),pdg(t), pL(t) are the power set-point for the wind turbines, photovoltaics, diesel generators, microturbines, fuel cells, battery energy storage, power exchange between the main grid and renewable microgrid, controllable distribution generators, and load levels. L(pdg(t)) is the power loss of the renewable microgrid. Nw,Ns,Nd,Nm,Nf,Nb,ND are the total number of wind turbines, photovoltaics, diesel generators, microturbines, fuel cells, battery energy storage, and load levels. In this paper, the controllable distribution generators (DGs) included diesel generators, microturbines, and fuel cells.


### 2.2.1. Estimation of Power Losses Coefficients

When the power losses of the renewable microgrid are considered, these power losses coefficients can be calculated as follows [28, 34]:(2)Lpdgt=∑m=1Ndg∑n=1Ndgpdg,mt.Bmnpdg,nt+∑n=1NdgBm0.pdg,m+B00,where B=[Bmn] and Bm0 and B00 are the power losses coefficients.


## 2.2.1. Estimation of Power Losses Coefficients

When the power losses of the renewable microgrid are considered, these power losses coefficients can be calculated as follows [28, 34]:(2)Lpdgt=∑m=1Ndg∑n=1Ndgpdg,mt.Bmnpdg,nt+∑n=1NdgBm0.pdg,m+B00,where B=[Bmn] and Bm0 and B00 are the power losses coefficients.


## 2.3. Subsystem Model of the Renewable Microgrid

This paper considers a renewable microgrid containing several types of devices, such as battery energy storages, distributed generators, wind turbines, and photovoltaics. A block diagram is given in Figure1.


### 2.3.1. The Model of Battery Energy Storage System

Battery energy storage can store energy when charging. On the other hand, it can supply energy to the load demand when discharging.The degradation cost of battery energy storage can be calculated as follows [33]:(3)Cb,ct=∑j=1Nb∑i=1trb.pbi,where rb is the cost coefficient which can be determined by battery energy storage charging or discharging power.p b , c ( t ) is the permitted rate of charge during a definite period of time and pb,d(t) is the permitted rate of discharge. The following equation can be expressed for a battery energy storage(4)pbt=uct.pb,ct+udt.pb,dt,where uc(t) or ud(t) indicate the working mode of the battery energy storage, uc(t),ud(t)∈0,1. Define uc(t)=1 or ud(t)=1 if the battery energy storage is charge or discharge. These variables satisfy(5)uct+udt≤1,t=1,…,T.The permitted rate of charge and discharge are limited by the following constrains:(6)0≤pb,ct≤pb,cmax,0≤pb,dt≤pb,dmax,where pb,cmax and pb,dmax are the maximum charging and discharging rates of battery energy storage.After charge or discharge, the state of charge of the battery energy storage can be calculated as follows [33]:(7)SOCbt=SOCbt-1+ηb,c.pb,ct-pb,dt/ηb,dQb,c.Δt,where ηb,c,ηb,d are the energy conversion coefficients. Qb,c is the energy capacity of the battery energy storage. On the other hand, the state of charge of the battery energy storage must be maintained in the following range. Δt is the operating time of the battery energy storage.(8)SOCb,min≤SOCbt≤SOCb,max,where SOCb,min,SOCb,max are required limits of the state of charge.


### 2.3.2. The Model of Controllable Distributed Generators

In this paper, the distribution generators (DGs) are included diesel generators, microturbines, and fuel cells.(A) Diesel Generators. In general, diesel generators serve as a backup power source. Meanwhile, the operating cost of diesel generators is formulated as follows:(9)Cdt=∑j=1Nd∑i=1tCd,fi+Cd,omi+Rdi.1-δdi-1.δdi,where Cd,f is the fuel cost of diesel generator; Cd,om is the operating and maintenance cost; Rd is the start-up cost of diesel generator.The fuel cost of diesel generators is modeled as a function of their actual output power(10)Cd,ft=a.pdt2+b.pdt+c,where pd(t) is the actual output power. a,b,c are the fuel consumption curve fitting coefficients.The operating and maintenance cost can be expressed as(11)Cd,omt=kd,om.pdt.Δt,where kd,om is the operating and maintenance cost coefficient. Δt is the operating time of diesel generator.The start-up cost of diesel generator is relative to its operating state(12)Rdt=Rd,h,Td,off≤Td,d+Td,c,Rd,c,Td,off>Td,d+Td,c,where Rd,h,Rd,c are the hot start-up cost and cold start-up cost. Td,c is the cooling time of diesel generator. Td,off,Td,d are the off time and shut down time of diesel generator.For a stable operation, the actual output power of diesel generator is limited by lower and upper bounds as follows:(13)pd,min≤pdt≤pd,max,where pd,min is the minimum active power of diesel generator. pd,max is the maximum power of diesel generator.The following constraint ensures diesel generators do not exceed their ramp limits:(14)-rd.pd,max≤pdt-pdt-1≤rd.pd,max,where rd is the ramp constraint coefficient.Due to the physical constraints that state shifting can only take place after a fixed time interval, the state variable of two adjacent times is as follows:(15)Td,ont-1-Td,u.δdt-1-δdt≥0,(16)Td,offt-1-Td,d.δdt-δdt-1≥0,where Td,on,Td,u are the cumulative uptime and the minimum turn on time.(B) Fuel Cells. Fuel cell directly converts chemical energy to electrical energy by electrochemical reactions, which is one of the most promising energy conversion technologies. The operating cost of fuel cells are formulated as(17)Cft=∑j=1Nf∑i=1tCf,fi+Cf,omi+Rfi.1-δfi-1.δfi,where Cf,f is the fuel cost of fuel cell; Cf,om is the operating and maintenance cost; Rf is the start-up cost of fuel cell.The fuel cost of fuel cells is modeled as(18)Cf,ft=cf.Tf.∑i=1tpfi+paηfi,where cf is the fuel price; Tf is the operating time of changing two adjacent states; pa is the consumption with active power for the auxiliary equipment. ηf is the efficiency of fuel cells.The operating and maintenance cost can be obtained by(19)Cf,omt=kf,om.Tf.∑i=1tpfi,where kf,om is the operating and maintenance cost coefficient.The start-up cost of fuel cells can be expressed as(20)Rft=Rf,h+βf.1-e-Tfofft/τf,where Rf,h is the hot start-up cost of fuel cells; Rf,h+βf is the cold start-up cost of fuel cells. Tfoff is the off time of fuel cells; τf is the cooling time of fuel cells.The actual output power of fuel cell is limited by lower and upper bounds as(21)pf,min≤pft≤pf,max,where pf,min,pf,max are the minimum and maximum active power of fuel cells.The following constraint ensures fuel cell does not exceed their ramp limits(22)-Δpf,D≤pft-pft-1≤Δpf,U,where Δpf,D,Δpf,U are the minimum and maximum ramp rate of fuel cells.The minimum on or minimum down time constraints for fuel cells can be expressed as(23)Tf,ont-1-Tf,u.δft-1-δft≥0,(24)Tf,offt-1-Tf,d.δft-δft-1≥0,where Tf,on is the cumulative uptime. Tf,u,Tf,d are the minimum turn on time and shut down time of fuel cells.The number of turn on and turn off can be expressed as(25)nf,st-spt≤Nf,max,where nf,st-sp(t) is the number of turn on and turn off in time t; Nf,max is the maximum number of turn on and turn off.(C) Microturbine. Microturbines have the higher power density, produce less noise, and emit much less pollutants, especially NOx. So microturbines are effective devices of converting the fuel energy into electrical energy. Meanwhile, the operating cost of microturbines are formulated as follows:(26)Cmt=∑j=1Nm∑i=1tCm,fi+Cm,omi+Rmi.1-δmi-1.δmi,where Cm,f is the fuel cost of microturbine; Cm,om is the operating and maintenance cost; Rm is the start-up cost of microturbine.The fuel cost of microturbines is modeled as(27)Cm,ft=cmpmtηmt,where cm is the fuel price of microturbines; ηm is the efficiency of microturbines.The operating and maintenance cost of microturbines can be expressed as(28)Cm,omt=km,om.pmt.Δt,where km,om is the operating and maintenance cost coefficient of microturbines.The start-up cost of microturbines can be expressed as(29)Rmt=Rm,h+βm.1-e-Tmofft/τm,where Rm,h is the hot start-up cost of microturbines; Rm,h+βm is the cold start-up cost of microturbines. Tmoff is the off time of microturbines; τm is the cooling time of microturbines.The actual output power of microturbines is limited by lower and upper bounds as follows:(30)pm,min≤pmt≤pm,max,where pm,min,pm,max are the minimum and maximum active power of microturbines.The following constraint ensures microturbines do not exceed their ramp limits(31)-Δpm,D≤pmt-pmt-1≤Δpm,U,where Δpm,D,Δpm,U are the minimum and maximum ramp rate of microturbines.The minimum on or minimum down time constraints for microturbines can be expressed as(32)Tm,ont-1-Tm,u.δmt-1-δmt≥0,(33)Tm,offt-1-Tm,d.δmt-δmt-1≥0,where Tm,on is the cumulative uptime; Tm,u,Tm,d are the minimum turn on time and shut down time of microturbines.


### 2.3.3. The Probability Model of Wind Turbines

Wind powers as renewable energy sources (RESs) are dependent on numerous factors such as wind velocity and efficiency. The approximate relationship between the wind power and wind speed can be expressed as(34)pwv=0,v<vi,v≥v0,v-vi.prvr-vi,vi≤v<vr,pr,vr≤v≤v0,where vi,vo,vr are the cut-in wind speed, cut-out wind speed, and rated wind speed. pr is the rated output power of wind turbines.Prior research has shown that the wind speed profile follows the Weibull distribution over time, which is [29, 32](35)Fvv=1-exp⁡-vck,where positive variables c and k are the scale parameter and shape parameter, respectively.Based on the characteristic of the wind power (34) and the probability distribution function (35), the probability distribution function of wind power is obtained by(36)FPwv=0,pwv<0,1-exp⁡-vi+vr-vipwv/prc,+exp⁡-v0ck,0≤pwv≤pr,1,pwv≥pr.


### 2.3.4. The Probability Model of Photovoltaics

The output power of photovoltaic can be calculated by its rated output power at the standard test condition and the operating ambient temperature [35](37)ps=Ac.η.Iβ=Ac.η.T.kt-T′.kt2,where Ac is the array surface area; η is the efficiency of the photovoltaic in realistic reporting conditions; Iβ is the irradiance on a surface with inclination β to the horizontal plane; T,T′ are the parameters that depend on inclination β, reflectance of the ground ρ, etc.(38)T=Rb+ρ.1-cos⁡β2+1+cos⁡β2-Rb.p.rd.H03600,(39)T′=1+cos⁡β2-Rb.q.rd.H03600,where Rb is the ratio of beam radiation on the tilted surface to that on a horizontal surface at any time and H0 is the solar radiation, for that day, both referring to a horizontal surface; rd is ratio, diffuse radiation in hour or diffuse in day; p,q are the parameters of the daily clearness index kt.Many researches have proved that cloudiness is the main factor affecting the difference between the values of solar radiation measured outside the atmosphere and on earthly surface. The daily clearness index can be obtained by [36](40)kt=ItI0,where It,I0 are the ratio of the irradiance on horizontal plane and the extraterrestrial total solar irradiance.The effect of clouds on terrestrial irradiance is the daily clearness index;kt can not be predicted with complete confidence; it must be treated as random variable(41)fkt,k¯t=Cktu.λ.γ1eλ.kt1-γ1.kt-1,(42)k¯t=νs.ktu-1.exp⁡νs.ktu-νs.ktl-1.exp⁡νs.ktlνs.exp⁡νs.ktu-exp⁡νs.ktl,where ktl,ktu are the lower and upper bounds of the observed range for kt; νs is the parameter of daily clearness index.(43)C=λ2.ktueλ.ktu-1-λ.ktu,(44)γ1=λ1+λ.ktu,(45)λ=2Γ-17.519.exp⁡-1.3118.Γ-1062.exp⁡-5.0426.Γktu,(46)Γ=ktuktu-k¯t,where C,γ1,λ are the parameters of daily clearness index.In particular, ifT>0,T′<0, the probability density function fps can be expressed as(47)fps=C.ku-α+α′/2-ktu.Ac.η.T′.α′.eλ.α+α′/2,ps∈0,psktu,0,ps∉0,psktu,while if T>0,T′>0, the probability density function can be obtained by(48)fps=C.ku-α-α′/2-ktu.Ac.η.T′.α′.eλ.α-α′/2,ps∈0,psktu,0,ps∉0,psktu,where α=T/T′,α′=α2-4.ps/η.T′.Ac.


### 2.3.5. The Model of Interaction with External Main Grid

In this paper, the renewable microgrid is connected to the external main grid and can trade energy with the main grid. The transaction incurs the following cost to the renewable microgrid:(49)Cgt=∑i=1tρg.pgi.Let the amount of energy be bounded by(50)-pg,max≤pgt≤pg,max(51)ρg=ρs,pg≤0.ρb,pg>0.where pg,max is the maximum transaction limits; ρs,ρb are the price coefficients to purchase and sell energy.


## 2.3.1. The Model of Battery Energy Storage System

Battery energy storage can store energy when charging. On the other hand, it can supply energy to the load demand when discharging.The degradation cost of battery energy storage can be calculated as follows [33]:(3)Cb,ct=∑j=1Nb∑i=1trb.pbi,where rb is the cost coefficient which can be determined by battery energy storage charging or discharging power.p b , c ( t ) is the permitted rate of charge during a definite period of time and pb,d(t) is the permitted rate of discharge. The following equation can be expressed for a battery energy storage(4)pbt=uct.pb,ct+udt.pb,dt,where uc(t) or ud(t) indicate the working mode of the battery energy storage, uc(t),ud(t)∈0,1. Define uc(t)=1 or ud(t)=1 if the battery energy storage is charge or discharge. These variables satisfy(5)uct+udt≤1,t=1,…,T.The permitted rate of charge and discharge are limited by the following constrains:(6)0≤pb,ct≤pb,cmax,0≤pb,dt≤pb,dmax,where pb,cmax and pb,dmax are the maximum charging and discharging rates of battery energy storage.After charge or discharge, the state of charge of the battery energy storage can be calculated as follows [33]:(7)SOCbt=SOCbt-1+ηb,c.pb,ct-pb,dt/ηb,dQb,c.Δt,where ηb,c,ηb,d are the energy conversion coefficients. Qb,c is the energy capacity of the battery energy storage. On the other hand, the state of charge of the battery energy storage must be maintained in the following range. Δt is the operating time of the battery energy storage.(8)SOCb,min≤SOCbt≤SOCb,max,where SOCb,min,SOCb,max are required limits of the state of charge.


## 2.3.2. The Model of Controllable Distributed Generators

In this paper, the distribution generators (DGs) are included diesel generators, microturbines, and fuel cells.(A) Diesel Generators. In general, diesel generators serve as a backup power source. Meanwhile, the operating cost of diesel generators is formulated as follows:(9)Cdt=∑j=1Nd∑i=1tCd,fi+Cd,omi+Rdi.1-δdi-1.δdi,where Cd,f is the fuel cost of diesel generator; Cd,om is the operating and maintenance cost; Rd is the start-up cost of diesel generator.The fuel cost of diesel generators is modeled as a function of their actual output power(10)Cd,ft=a.pdt2+b.pdt+c,where pd(t) is the actual output power. a,b,c are the fuel consumption curve fitting coefficients.The operating and maintenance cost can be expressed as(11)Cd,omt=kd,om.pdt.Δt,where kd,om is the operating and maintenance cost coefficient. Δt is the operating time of diesel generator.The start-up cost of diesel generator is relative to its operating state(12)Rdt=Rd,h,Td,off≤Td,d+Td,c,Rd,c,Td,off>Td,d+Td,c,where Rd,h,Rd,c are the hot start-up cost and cold start-up cost. Td,c is the cooling time of diesel generator. Td,off,Td,d are the off time and shut down time of diesel generator.For a stable operation, the actual output power of diesel generator is limited by lower and upper bounds as follows:(13)pd,min≤pdt≤pd,max,where pd,min is the minimum active power of diesel generator. pd,max is the maximum power of diesel generator.The following constraint ensures diesel generators do not exceed their ramp limits:(14)-rd.pd,max≤pdt-pdt-1≤rd.pd,max,where rd is the ramp constraint coefficient.Due to the physical constraints that state shifting can only take place after a fixed time interval, the state variable of two adjacent times is as follows:(15)Td,ont-1-Td,u.δdt-1-δdt≥0,(16)Td,offt-1-Td,d.δdt-δdt-1≥0,where Td,on,Td,u are the cumulative uptime and the minimum turn on time.(B) Fuel Cells. Fuel cell directly converts chemical energy to electrical energy by electrochemical reactions, which is one of the most promising energy conversion technologies. The operating cost of fuel cells are formulated as(17)Cft=∑j=1Nf∑i=1tCf,fi+Cf,omi+Rfi.1-δfi-1.δfi,where Cf,f is the fuel cost of fuel cell; Cf,om is the operating and maintenance cost; Rf is the start-up cost of fuel cell.The fuel cost of fuel cells is modeled as(18)Cf,ft=cf.Tf.∑i=1tpfi+paηfi,where cf is the fuel price; Tf is the operating time of changing two adjacent states; pa is the consumption with active power for the auxiliary equipment. ηf is the efficiency of fuel cells.The operating and maintenance cost can be obtained by(19)Cf,omt=kf,om.Tf.∑i=1tpfi,where kf,om is the operating and maintenance cost coefficient.The start-up cost of fuel cells can be expressed as(20)Rft=Rf,h+βf.1-e-Tfofft/τf,where Rf,h is the hot start-up cost of fuel cells; Rf,h+βf is the cold start-up cost of fuel cells. Tfoff is the off time of fuel cells; τf is the cooling time of fuel cells.The actual output power of fuel cell is limited by lower and upper bounds as(21)pf,min≤pft≤pf,max,where pf,min,pf,max are the minimum and maximum active power of fuel cells.The following constraint ensures fuel cell does not exceed their ramp limits(22)-Δpf,D≤pft-pft-1≤Δpf,U,where Δpf,D,Δpf,U are the minimum and maximum ramp rate of fuel cells.The minimum on or minimum down time constraints for fuel cells can be expressed as(23)Tf,ont-1-Tf,u.δft-1-δft≥0,(24)Tf,offt-1-Tf,d.δft-δft-1≥0,where Tf,on is the cumulative uptime. Tf,u,Tf,d are the minimum turn on time and shut down time of fuel cells.The number of turn on and turn off can be expressed as(25)nf,st-spt≤Nf,max,where nf,st-sp(t) is the number of turn on and turn off in time t; Nf,max is the maximum number of turn on and turn off.(C) Microturbine. Microturbines have the higher power density, produce less noise, and emit much less pollutants, especially NOx. So microturbines are effective devices of converting the fuel energy into electrical energy. Meanwhile, the operating cost of microturbines are formulated as follows:(26)Cmt=∑j=1Nm∑i=1tCm,fi+Cm,omi+Rmi.1-δmi-1.δmi,where Cm,f is the fuel cost of microturbine; Cm,om is the operating and maintenance cost; Rm is the start-up cost of microturbine.The fuel cost of microturbines is modeled as(27)Cm,ft=cmpmtηmt,where cm is the fuel price of microturbines; ηm is the efficiency of microturbines.The operating and maintenance cost of microturbines can be expressed as(28)Cm,omt=km,om.pmt.Δt,where km,om is the operating and maintenance cost coefficient of microturbines.The start-up cost of microturbines can be expressed as(29)Rmt=Rm,h+βm.1-e-Tmofft/τm,where Rm,h is the hot start-up cost of microturbines; Rm,h+βm is the cold start-up cost of microturbines. Tmoff is the off time of microturbines; τm is the cooling time of microturbines.The actual output power of microturbines is limited by lower and upper bounds as follows:(30)pm,min≤pmt≤pm,max,where pm,min,pm,max are the minimum and maximum active power of microturbines.The following constraint ensures microturbines do not exceed their ramp limits(31)-Δpm,D≤pmt-pmt-1≤Δpm,U,where Δpm,D,Δpm,U are the minimum and maximum ramp rate of microturbines.The minimum on or minimum down time constraints for microturbines can be expressed as(32)Tm,ont-1-Tm,u.δmt-1-δmt≥0,(33)Tm,offt-1-Tm,d.δmt-δmt-1≥0,where Tm,on is the cumulative uptime; Tm,u,Tm,d are the minimum turn on time and shut down time of microturbines.


## 2.3.3. The Probability Model of Wind Turbines

Wind powers as renewable energy sources (RESs) are dependent on numerous factors such as wind velocity and efficiency. The approximate relationship between the wind power and wind speed can be expressed as(34)pwv=0,v<vi,v≥v0,v-vi.prvr-vi,vi≤v<vr,pr,vr≤v≤v0,where vi,vo,vr are the cut-in wind speed, cut-out wind speed, and rated wind speed. pr is the rated output power of wind turbines.Prior research has shown that the wind speed profile follows the Weibull distribution over time, which is [29, 32](35)Fvv=1-exp⁡-vck,where positive variables c and k are the scale parameter and shape parameter, respectively.Based on the characteristic of the wind power (34) and the probability distribution function (35), the probability distribution function of wind power is obtained by(36)FPwv=0,pwv<0,1-exp⁡-vi+vr-vipwv/prc,+exp⁡-v0ck,0≤pwv≤pr,1,pwv≥pr.


## 2.3.4. The Probability Model of Photovoltaics

The output power of photovoltaic can be calculated by its rated output power at the standard test condition and the operating ambient temperature [35](37)ps=Ac.η.Iβ=Ac.η.T.kt-T′.kt2,where Ac is the array surface area; η is the efficiency of the photovoltaic in realistic reporting conditions; Iβ is the irradiance on a surface with inclination β to the horizontal plane; T,T′ are the parameters that depend on inclination β, reflectance of the ground ρ, etc.(38)T=Rb+ρ.1-cos⁡β2+1+cos⁡β2-Rb.p.rd.H03600,(39)T′=1+cos⁡β2-Rb.q.rd.H03600,where Rb is the ratio of beam radiation on the tilted surface to that on a horizontal surface at any time and H0 is the solar radiation, for that day, both referring to a horizontal surface; rd is ratio, diffuse radiation in hour or diffuse in day; p,q are the parameters of the daily clearness index kt.Many researches have proved that cloudiness is the main factor affecting the difference between the values of solar radiation measured outside the atmosphere and on earthly surface. The daily clearness index can be obtained by [36](40)kt=ItI0,where It,I0 are the ratio of the irradiance on horizontal plane and the extraterrestrial total solar irradiance.The effect of clouds on terrestrial irradiance is the daily clearness index;kt can not be predicted with complete confidence; it must be treated as random variable(41)fkt,k¯t=Cktu.λ.γ1eλ.kt1-γ1.kt-1,(42)k¯t=νs.ktu-1.exp⁡νs.ktu-νs.ktl-1.exp⁡νs.ktlνs.exp⁡νs.ktu-exp⁡νs.ktl,where ktl,ktu are the lower and upper bounds of the observed range for kt; νs is the parameter of daily clearness index.(43)C=λ2.ktueλ.ktu-1-λ.ktu,(44)γ1=λ1+λ.ktu,(45)λ=2Γ-17.519.exp⁡-1.3118.Γ-1062.exp⁡-5.0426.Γktu,(46)Γ=ktuktu-k¯t,where C,γ1,λ are the parameters of daily clearness index.In particular, ifT>0,T′<0, the probability density function fps can be expressed as(47)fps=C.ku-α+α′/2-ktu.Ac.η.T′.α′.eλ.α+α′/2,ps∈0,psktu,0,ps∉0,psktu,while if T>0,T′>0, the probability density function can be obtained by(48)fps=C.ku-α-α′/2-ktu.Ac.η.T′.α′.eλ.α-α′/2,ps∈0,psktu,0,ps∉0,psktu,where α=T/T′,α′=α2-4.ps/η.T′.Ac.


## 2.3.5. The Model of Interaction with External Main Grid

In this paper, the renewable microgrid is connected to the external main grid and can trade energy with the main grid. The transaction incurs the following cost to the renewable microgrid:(49)Cgt=∑i=1tρg.pgi.Let the amount of energy be bounded by(50)-pg,max≤pgt≤pg,max(51)ρg=ρs,pg≤0.ρb,pg>0.where pg,max is the maximum transaction limits; ρs,ρb are the price coefficients to purchase and sell energy.


## 2.4. Dynamic Operation Model of the Renewable Microgrid

Dynamic economic operation management of the renewable microgrid is to determine output power of distribution generators, in order to minimize total operation cost of the renewable microgrid and meet the dynamic operation constrains.The total operating cost of the renewable microgrid can be defined as(52)min⁡Ft=∑t=1TCb,ct+Cdt+Cft+Cmt+Cgt=minΠ=u1,…,uT⁡E∑t=1TCtxt,ut,wt,where u(t) is the decision or control variables (specifically, action at state x(t) in period t); x(t) is the state variables of the renewable microgrid; w(t) is the random variable such as the output power of wind turbine or photovoltaics.u ( t ) , x ( t ) , w ( t ) can be defined as follows:(53)ut=Δpbt,Δpdt,Δpmt,Δpft,Δpgt,(54)xt=SOCbt,pbt,pdt,pmt,pft,pgt,(55)wt=pwt,pst.The state variables of renewable microgrid int+1 period can be defined as follows:(56)xt+1=fxt,ut,wt=SOCbt+1,pbt+Δpbt,pdt+Δpdt,pmt+Δpmt,pft+Δpft,pgt+Δpgt.The operation constraints of the renewable microgrid can be expressed as power balance (1)-(2), battery energy storage limits (5)-(8), distributed generations limits (13)-(16), (21)-(25), (30)-(33), and interaction with external main grid limits (50).


### 2.4.1. The Standard Formulation of the Dynamic Operation Model

According to dynamic programming formulation, the dynamic operation model of the renewable microgrid can be formulated as(57)Jtxt=minut∈Ut⁡Ctxt,ut,wt+EJt+1xt+1.The operation constraints can be expressed as(58)xt+1=fxt,ut,wt,(59)ctExt,ut,wt=0,cTExT=0,(60)ctIxt,ut,wt≤0,cTIxT≤0,where ctE,cTE,ctI,cTI are, respectively, the equality and inequality constraints.Meanwhile, the constraint of condition should be satisfied as follows:(61)JT+1xT+1≡0.


## 2.4.1. The Standard Formulation of the Dynamic Operation Model

According to dynamic programming formulation, the dynamic operation model of the renewable microgrid can be formulated as(57)Jtxt=minut∈Ut⁡Ctxt,ut,wt+EJt+1xt+1.The operation constraints can be expressed as(58)xt+1=fxt,ut,wt,(59)ctExt,ut,wt=0,cTExT=0,(60)ctIxt,ut,wt≤0,cTIxT≤0,where ctE,cTE,ctI,cTI are, respectively, the equality and inequality constraints.Meanwhile, the constraint of condition should be satisfied as follows:(61)JT+1xT+1≡0.


## 3. Solution Methodology

Based on Section2, dynamic operation management of the renewable microgrid can be regarded as a discrete time system under uncertainty. In order to solve the proposed problem, the random variables can be realized by the scenario analysis technique. Then, a feedback control scheme based on the dynamic programming is proposed to solve the model and obtain the optimal solution. The flowchart of the whole process is given in Figure 2.Figure 2
The flowchart of the proposed improved dynamic programming algorithm.


### 3.1. The Scenario Analysis Technique

In this paper, a discrete set of scenarios can be used to represent the probability realization of the output power of wind turbines or photovoltaics. On the other hand, these scenarios are generated using the Roulette wheel mechanism and Monte Carlo simulation method [19]. And the probability distribution function of the random variables can be obtained by (36) and (47)-(48).The Lattice Monte Carlo simulation(LMCS) can be used to generate the random numbers [29](62)∑j=1rkjnj∑i=1NLvimod1,kj=1,…,NL;j=1,…,r,where r is the number of random variables; NL is the number of random sampling; vi is vector with dimension d.According to the desired preciseness, the probability distribution functions (36) and (47)-(48) are divided into nl class intervals. Each class interval determines mean value Inti,i=1,…,nl; and each interval is associated with a probability denoted by βi=1,…,nl. Meanwhile, the probabilities of different intervals are normalized in which their summations become equal to unity.Therefore, each scenario comprises a vector identifying the output power of wind turbine or photovoltaics:(63)St,sc=W1,t,scw,…,Wnl,t,scw,W1,t,scs,…,Wnl,t,scs,t=1,…,T,sc=1,…,Nsc,where Wnl,t,scw,Wnl,t,scs are binary parameters indicating where nl wind power interval or photovoltaic power output whether are selected in scenarios. On the other hand, comparing the random number which follows the LMCS strategy and the probability βi, the binary parameter can be selected.Thus, the output power of wind turbine or photovoltaics for each scenario can be obtained by(64)pwt=pw,sct=∑uw=1nlWuw,t,scw.Intuw,tw,t=1,…,T,sc=1,…,Nsc,(65)pst=ps,sct=∑us=1nlWus,t,scs.Intus,ts,t=1,…,T,sc=1,…,Nsc.The normalized probability of each scenario can be expressed as follows:(66)πsc=∏t=1T∑uw=1nlWuw,t,scwβuw,t∑us=1nlWus,t,scsβus,t∑sc=1Ns∏t=1T∑uw=1nlWuw,t,scwβuw,t∑us=1nlWus,t,scwβus,t,t=1,…,T,sc=1,…,Nsc,where Nsc is the number of scenarios.


### 3.2. The Feedback Control Scheme Based on the Dynamic Programming

In the following section, we formulate the dynamic operation management of renewable microgrid which satisfies dynamic operations constraints. Firstly, the Lagrange multipliers and interior point method will be presented for the equality and inequality operation constraint. Then, a feedback control scheme based on the dynamic programming is proposed to solve the model and obtain the optimal solution.


#### 3.2.1. The Sequential Quadratic Programming Subproblem Formulation

In this paper, the sequential quadratic programming will be presented for the formulation (52)-(61). The purpose of sequential quadratic programming is to leave the model with only linear and quadratic terms. The equivalent formulation can be obtained as follows:(67)min⁡Γq=∑t=1TΓt+x¯tT.yt+u¯tT.zt+x¯tT.Qt.x¯t+2.x¯tT.Rt.u¯t+u¯tT.St.u¯t2+ΓT+x¯TT.yT+x¯TT.QT.x¯T2,where v(t) is a constant. u¯(t),x¯(t) are the perturbation variables. u^(t),x^(t) are the nominal values of x(t),u(t). And these variables can be introduced and defined as(68)Γt=xtT.yt+xtT.Qt.xt+vt,(69)ΓT=xTT.yT+xTT.QT.xT+vT,(70)u¯t=ut-u^t,(71)x¯t=xt-x^t.The operation constraints can be expressed as(72)x¯t+1=At.x¯t+Bt.u¯t,(73)AEt.x¯t+BEt.u¯t+c^Et=0,AET.x¯T+c^ET=0,(74)AIt.x¯t+BIt.u¯t+c^It≤0,AIT.xT+c^IT≤0,where y(t),z(t),Q(t),R(t),S(t) and A(t),B(t),AE(t),BE(t),c^E(t),AI(t),BI(t),c^I(t) can be obtained in Appendix A.


#### 3.2.2. Interior Point Methods

According to the sequential quadratic programming formulations (67), by applying the Lagrange multipliers and slacking variables to proceed the equalities and inequalities operation constraints, a new formulation is shown as(75)Lqk=∑t=1Tv¯kt+x~tT.y¯t+u~tT.z¯t+x~tT.Q¯t.x~t+2.x~tT.R¯t.u~t+u~tT.S¯t.u~t2+A¯Et.x~t+B¯Et.u~t+c¯EtT.λ~Et+A¯It.x~t+B¯It.u~t+c¯ItT.λ~It-λ~ItT.G¯t.λ~It2+x~TT.y¯T+x~TT.QT.x~T2+A¯ET.x~T+c¯ETT.λ~ET+A¯IT.x~T+c¯ITT.λ~IT-λ~ITT.G¯T.λ~IT2,where k is considered as the kth iteration. v¯k(t) is a constant.Meanwhile, the parameters can be defined as follows:(76)z~=u~1…u~T-1x~1T,ξ~=λ~E1…λ~Etλ~ETT,(77)λ~=λ~I1…λ~Itλ~ITT,s~=s~I1…s~Its~ITT.The approach in solving (z~,ξ~,λ~,s~) is given in Appendix B.


#### 3.2.3. Total Algorithm Procedures

In this paper, an algorithm is presented based on an adaptation of Mehrotra’s predictor-corrector algorithm [37]. Given the initial parameters δ>0,k=0 and (z0,ξ0,λ0,s0) the residual terms are defined as(78)r1k2=∑t=1T∂Lqk∂u~tT∂Lqk∂u~t,(79)r2k2=∑t=1TAEt.xtk+BEt.utk+c^EtT.AEt.xtk+BEt.utk+c^Et,(80)r3k2=∑t=1TAIt.xtk+BIt.utk+c^It+sIktT.AIt.xtk+BIt.utk+c^It+sIkt,(81)r4k2=∑t=1Tmin⁡λIkt,sIktT.min⁡λIkt,sIkt.As the solution set approaches optimum point, the summation of all residual term should approach zero.The total algorithm can be expressed as follows.Step 1.
Calculate the residual term as(82)rk=r1k2+r2k2+r3k2+r4k2,where the expressions for r1k,r2k,r3k,r4k are given in (78)-(81). If rk<δ, then go to Step 5.Step 2.
Calculate the affine-scaling direction{z~a,ξ~a,λ~a,s~a} by using the method introduced in Appendix B.d k , c I k ( t ) can be defined as(83)dk=-diagλk.diagsk.e,(84)cIkt=AIt.xkt+BIt.ukt+c~It.The scalarsαa,μa, and σ can be calculated as follows:(85)μk=λkT.skn2,where n2 is the dimension of sk.(86)αa=argmax⁡α∈0,1∣λk+α.λ~a≥0,sk+αs~a≥0,(87)μa=λk+α.λ~aT.1n2.sk+αa.s~a.Meanwhile, the centering parameter is defined as follows:(88)σ=μaμk3.Step 3.
Calculate the combined predictor-centering-corrector direction{z~p,ξ~p,λ~p,s~p} and the scalar αmax as the method introduced in Appendix B.Let(89)dk=-diagλk.diagsk.e-diagλ~a.diags~a.e+σ.μk.e,(90)cIkt=AIt.xkt+BIt.ukt+c^It+diagλIkt-1.σ.μk.e-diagλ~at.diags~at.e,where e=1…1.Calculate the scalar as follows:(91)αmax=argmax⁡α∈0,1∣λk+α.λ~p≥0,sk+α.s~p≥0.Step 4.
Calculate(zk+1,ξk+1,λk+1,sk+1) from the equation(92)zk+1,ξk+1,λk+1,sk+1=zk,ξk,λk,sk+αz~p,ξ~p,λ~p,s~p,where α=min⁡(0.995.αmax,1). Then set k←k+1 and return to Step 1.Step 5.
Calculate the search directionu~(t),x~(1),λ~E(t),λ~I(t),s~I(t) and the optimal solution(93)u∗t=ukt+u~t,x∗1=xk1+x~1.
Calculate all state variables by using the equation(94)x∗t+1=At.x∗t+Bt.u∗t,t=1,…,T.


## 3.1. The Scenario Analysis Technique

In this paper, a discrete set of scenarios can be used to represent the probability realization of the output power of wind turbines or photovoltaics. On the other hand, these scenarios are generated using the Roulette wheel mechanism and Monte Carlo simulation method [19]. And the probability distribution function of the random variables can be obtained by (36) and (47)-(48).The Lattice Monte Carlo simulation(LMCS) can be used to generate the random numbers [29](62)∑j=1rkjnj∑i=1NLvimod1,kj=1,…,NL;j=1,…,r,where r is the number of random variables; NL is the number of random sampling; vi is vector with dimension d.According to the desired preciseness, the probability distribution functions (36) and (47)-(48) are divided into nl class intervals. Each class interval determines mean value Inti,i=1,…,nl; and each interval is associated with a probability denoted by βi=1,…,nl. Meanwhile, the probabilities of different intervals are normalized in which their summations become equal to unity.Therefore, each scenario comprises a vector identifying the output power of wind turbine or photovoltaics:(63)St,sc=W1,t,scw,…,Wnl,t,scw,W1,t,scs,…,Wnl,t,scs,t=1,…,T,sc=1,…,Nsc,where Wnl,t,scw,Wnl,t,scs are binary parameters indicating where nl wind power interval or photovoltaic power output whether are selected in scenarios. On the other hand, comparing the random number which follows the LMCS strategy and the probability βi, the binary parameter can be selected.Thus, the output power of wind turbine or photovoltaics for each scenario can be obtained by(64)pwt=pw,sct=∑uw=1nlWuw,t,scw.Intuw,tw,t=1,…,T,sc=1,…,Nsc,(65)pst=ps,sct=∑us=1nlWus,t,scs.Intus,ts,t=1,…,T,sc=1,…,Nsc.The normalized probability of each scenario can be expressed as follows:(66)πsc=∏t=1T∑uw=1nlWuw,t,scwβuw,t∑us=1nlWus,t,scsβus,t∑sc=1Ns∏t=1T∑uw=1nlWuw,t,scwβuw,t∑us=1nlWus,t,scwβus,t,t=1,…,T,sc=1,…,Nsc,where Nsc is the number of scenarios.


## 3.2. The Feedback Control Scheme Based on the Dynamic Programming

In the following section, we formulate the dynamic operation management of renewable microgrid which satisfies dynamic operations constraints. Firstly, the Lagrange multipliers and interior point method will be presented for the equality and inequality operation constraint. Then, a feedback control scheme based on the dynamic programming is proposed to solve the model and obtain the optimal solution.


### 3.2.1. The Sequential Quadratic Programming Subproblem Formulation

In this paper, the sequential quadratic programming will be presented for the formulation (52)-(61). The purpose of sequential quadratic programming is to leave the model with only linear and quadratic terms. The equivalent formulation can be obtained as follows:(67)min⁡Γq=∑t=1TΓt+x¯tT.yt+u¯tT.zt+x¯tT.Qt.x¯t+2.x¯tT.Rt.u¯t+u¯tT.St.u¯t2+ΓT+x¯TT.yT+x¯TT.QT.x¯T2,where v(t) is a constant. u¯(t),x¯(t) are the perturbation variables. u^(t),x^(t) are the nominal values of x(t),u(t). And these variables can be introduced and defined as(68)Γt=xtT.yt+xtT.Qt.xt+vt,(69)ΓT=xTT.yT+xTT.QT.xT+vT,(70)u¯t=ut-u^t,(71)x¯t=xt-x^t.The operation constraints can be expressed as(72)x¯t+1=At.x¯t+Bt.u¯t,(73)AEt.x¯t+BEt.u¯t+c^Et=0,AET.x¯T+c^ET=0,(74)AIt.x¯t+BIt.u¯t+c^It≤0,AIT.xT+c^IT≤0,where y(t),z(t),Q(t),R(t),S(t) and A(t),B(t),AE(t),BE(t),c^E(t),AI(t),BI(t),c^I(t) can be obtained in Appendix A.


### 3.2.2. Interior Point Methods

According to the sequential quadratic programming formulations (67), by applying the Lagrange multipliers and slacking variables to proceed the equalities and inequalities operation constraints, a new formulation is shown as(75)Lqk=∑t=1Tv¯kt+x~tT.y¯t+u~tT.z¯t+x~tT.Q¯t.x~t+2.x~tT.R¯t.u~t+u~tT.S¯t.u~t2+A¯Et.x~t+B¯Et.u~t+c¯EtT.λ~Et+A¯It.x~t+B¯It.u~t+c¯ItT.λ~It-λ~ItT.G¯t.λ~It2+x~TT.y¯T+x~TT.QT.x~T2+A¯ET.x~T+c¯ETT.λ~ET+A¯IT.x~T+c¯ITT.λ~IT-λ~ITT.G¯T.λ~IT2,where k is considered as the kth iteration. v¯k(t) is a constant.Meanwhile, the parameters can be defined as follows:(76)z~=u~1…u~T-1x~1T,ξ~=λ~E1…λ~Etλ~ETT,(77)λ~=λ~I1…λ~Itλ~ITT,s~=s~I1…s~Its~ITT.The approach in solving (z~,ξ~,λ~,s~) is given in Appendix B.


### 3.2.3. Total Algorithm Procedures

In this paper, an algorithm is presented based on an adaptation of Mehrotra’s predictor-corrector algorithm [37]. Given the initial parameters δ>0,k=0 and (z0,ξ0,λ0,s0) the residual terms are defined as(78)r1k2=∑t=1T∂Lqk∂u~tT∂Lqk∂u~t,(79)r2k2=∑t=1TAEt.xtk+BEt.utk+c^EtT.AEt.xtk+BEt.utk+c^Et,(80)r3k2=∑t=1TAIt.xtk+BIt.utk+c^It+sIktT.AIt.xtk+BIt.utk+c^It+sIkt,(81)r4k2=∑t=1Tmin⁡λIkt,sIktT.min⁡λIkt,sIkt.As the solution set approaches optimum point, the summation of all residual term should approach zero.The total algorithm can be expressed as follows.Step 1.
Calculate the residual term as(82)rk=r1k2+r2k2+r3k2+r4k2,where the expressions for r1k,r2k,r3k,r4k are given in (78)-(81). If rk<δ, then go to Step 5.Step 2.
Calculate the affine-scaling direction{z~a,ξ~a,λ~a,s~a} by using the method introduced in Appendix B.d k , c I k ( t ) can be defined as(83)dk=-diagλk.diagsk.e,(84)cIkt=AIt.xkt+BIt.ukt+c~It.The scalarsαa,μa, and σ can be calculated as follows:(85)μk=λkT.skn2,where n2 is the dimension of sk.(86)αa=argmax⁡α∈0,1∣λk+α.λ~a≥0,sk+αs~a≥0,(87)μa=λk+α.λ~aT.1n2.sk+αa.s~a.Meanwhile, the centering parameter is defined as follows:(88)σ=μaμk3.Step 3.
Calculate the combined predictor-centering-corrector direction{z~p,ξ~p,λ~p,s~p} and the scalar αmax as the method introduced in Appendix B.Let(89)dk=-diagλk.diagsk.e-diagλ~a.diags~a.e+σ.μk.e,(90)cIkt=AIt.xkt+BIt.ukt+c^It+diagλIkt-1.σ.μk.e-diagλ~at.diags~at.e,where e=1…1.Calculate the scalar as follows:(91)αmax=argmax⁡α∈0,1∣λk+α.λ~p≥0,sk+α.s~p≥0.Step 4.
Calculate(zk+1,ξk+1,λk+1,sk+1) from the equation(92)zk+1,ξk+1,λk+1,sk+1=zk,ξk,λk,sk+αz~p,ξ~p,λ~p,s~p,where α=min⁡(0.995.αmax,1). Then set k←k+1 and return to Step 1.Step 5.
Calculate the search directionu~(t),x~(1),λ~E(t),λ~I(t),s~I(t) and the optimal solution(93)u∗t=ukt+u~t,x∗1=xk1+x~1.
Calculate all state variables by using the equation(94)x∗t+1=At.x∗t+Bt.u∗t,t=1,…,T.


## 3.2.1. The Sequential Quadratic Programming Subproblem Formulation

In this paper, the sequential quadratic programming will be presented for the formulation (52)-(61). The purpose of sequential quadratic programming is to leave the model with only linear and quadratic terms. The equivalent formulation can be obtained as follows:(67)min⁡Γq=∑t=1TΓt+x¯tT.yt+u¯tT.zt+x¯tT.Qt.x¯t+2.x¯tT.Rt.u¯t+u¯tT.St.u¯t2+ΓT+x¯TT.yT+x¯TT.QT.x¯T2,where v(t) is a constant. u¯(t),x¯(t) are the perturbation variables. u^(t),x^(t) are the nominal values of x(t),u(t). And these variables can be introduced and defined as(68)Γt=xtT.yt+xtT.Qt.xt+vt,(69)ΓT=xTT.yT+xTT.QT.xT+vT,(70)u¯t=ut-u^t,(71)x¯t=xt-x^t.The operation constraints can be expressed as(72)x¯t+1=At.x¯t+Bt.u¯t,(73)AEt.x¯t+BEt.u¯t+c^Et=0,AET.x¯T+c^ET=0,(74)AIt.x¯t+BIt.u¯t+c^It≤0,AIT.xT+c^IT≤0,where y(t),z(t),Q(t),R(t),S(t) and A(t),B(t),AE(t),BE(t),c^E(t),AI(t),BI(t),c^I(t) can be obtained in Appendix A.


## 3.2.2. Interior Point Methods

According to the sequential quadratic programming formulations (67), by applying the Lagrange multipliers and slacking variables to proceed the equalities and inequalities operation constraints, a new formulation is shown as(75)Lqk=∑t=1Tv¯kt+x~tT.y¯t+u~tT.z¯t+x~tT.Q¯t.x~t+2.x~tT.R¯t.u~t+u~tT.S¯t.u~t2+A¯Et.x~t+B¯Et.u~t+c¯EtT.λ~Et+A¯It.x~t+B¯It.u~t+c¯ItT.λ~It-λ~ItT.G¯t.λ~It2+x~TT.y¯T+x~TT.QT.x~T2+A¯ET.x~T+c¯ETT.λ~ET+A¯IT.x~T+c¯ITT.λ~IT-λ~ITT.G¯T.λ~IT2,where k is considered as the kth iteration. v¯k(t) is a constant.Meanwhile, the parameters can be defined as follows:(76)z~=u~1…u~T-1x~1T,ξ~=λ~E1…λ~Etλ~ETT,(77)λ~=λ~I1…λ~Itλ~ITT,s~=s~I1…s~Its~ITT.The approach in solving (z~,ξ~,λ~,s~) is given in Appendix B.


## 3.2.3. Total Algorithm Procedures

In this paper, an algorithm is presented based on an adaptation of Mehrotra’s predictor-corrector algorithm [37]. Given the initial parameters δ>0,k=0 and (z0,ξ0,λ0,s0) the residual terms are defined as(78)r1k2=∑t=1T∂Lqk∂u~tT∂Lqk∂u~t,(79)r2k2=∑t=1TAEt.xtk+BEt.utk+c^EtT.AEt.xtk+BEt.utk+c^Et,(80)r3k2=∑t=1TAIt.xtk+BIt.utk+c^It+sIktT.AIt.xtk+BIt.utk+c^It+sIkt,(81)r4k2=∑t=1Tmin⁡λIkt,sIktT.min⁡λIkt,sIkt.As the solution set approaches optimum point, the summation of all residual term should approach zero.The total algorithm can be expressed as follows.Step 1.
Calculate the residual term as(82)rk=r1k2+r2k2+r3k2+r4k2,where the expressions for r1k,r2k,r3k,r4k are given in (78)-(81). If rk<δ, then go to Step 5.Step 2.
Calculate the affine-scaling direction{z~a,ξ~a,λ~a,s~a} by using the method introduced in Appendix B.d k , c I k ( t ) can be defined as(83)dk=-diagλk.diagsk.e,(84)cIkt=AIt.xkt+BIt.ukt+c~It.The scalarsαa,μa, and σ can be calculated as follows:(85)μk=λkT.skn2,where n2 is the dimension of sk.(86)αa=argmax⁡α∈0,1∣λk+α.λ~a≥0,sk+αs~a≥0,(87)μa=λk+α.λ~aT.1n2.sk+αa.s~a.Meanwhile, the centering parameter is defined as follows:(88)σ=μaμk3.Step 3.
Calculate the combined predictor-centering-corrector direction{z~p,ξ~p,λ~p,s~p} and the scalar αmax as the method introduced in Appendix B.Let(89)dk=-diagλk.diagsk.e-diagλ~a.diags~a.e+σ.μk.e,(90)cIkt=AIt.xkt+BIt.ukt+c^It+diagλIkt-1.σ.μk.e-diagλ~at.diags~at.e,where e=1…1.Calculate the scalar as follows:(91)αmax=argmax⁡α∈0,1∣λk+α.λ~p≥0,sk+α.s~p≥0.Step 4.
Calculate(zk+1,ξk+1,λk+1,sk+1) from the equation(92)zk+1,ξk+1,λk+1,sk+1=zk,ξk,λk,sk+αz~p,ξ~p,λ~p,s~p,where α=min⁡(0.995.αmax,1). Then set k←k+1 and return to Step 1.Step 5.
Calculate the search directionu~(t),x~(1),λ~E(t),λ~I(t),s~I(t) and the optimal solution(93)u∗t=ukt+u~t,x∗1=xk1+x~1.
Calculate all state variables by using the equation(94)x∗t+1=At.x∗t+Bt.u∗t,t=1,…,T.


## 4. Simulation

In this section, the renewable microgrid in the simulation is shown in Figure1; it operates in parallel to the main power grid or as an island and comprises wind turbines with maximum power of 13kw, PV panels with maximum power of 15kw, and distributed generators such as microturbines, diesel generators, and full cells. A battery energy storage is included, bounded between 20 and 60 kWh and with maximal charge and discharge rates, respectively, 30 and -30kw. The charge and discharging efficiencies are both equal to 0.85. Tables 1–4 describe the distributed generators units parameters, based on data provided in [29, 32]. The simulations are carried out in the MATLAB environment on an Intel Core 2 Duo 3.00GHz running Windows 7. In the simulations study, we chose the sampling time of 1 hour. And simulations are performed over a horizon of 24 hours.Table 1
The capacities parameters of microgenerators (MGs) [29, 32].
The capacities of MGs Wind power photovoltaics Battery Diesel generator microturbine full cell Unite(kw) 13 15 60 30 30 50Table 2
The technologies parameters of distributed generators [29, 32].
Parameters Parameters a Parameters b Parameters c Diesel generator 0.4333 0.2333 0.0074  Parameters Fuel price parameter Efficiency parameter Auxiliary equipment  Full cell 0.0164 0.6 1.228  Microturbine 0.0164 0.56 -Table 3
The technologies parameters of battery storage system [29, 32].
Upper energy limits Lower energy limits Charging efficiency Discharging efficiency 20 (kw) 60(kw) 0.85 0.85Table 4
The electricity interaction parameters with dispatching networks [10, 19].
The upper energy limits The lower energy limits 40 (kw) -40(kw)The electricity usage, wind power, and solar generation data in the simulations have been provided by National Renewable Energy Laboratory. The electricity usage data is from the utility operator customers with peak usage over 90kw; the maximum power demand of load is 90.1kw. The daily electricity usage (from Miami, FI, on a certain day) is shown in Figure3.Figure 3
Load profiles for a period of one day.


### 4.1. Scenario Generation and Reduction Results

The approach used for generation of scenarios may have an impact on the determination of operational cost of renewable microgrid. In the study case, the dates from wind power profiles of the 24 hours are used as scenarios in the dynamic operation management of renewable microgrid. The scenarios of wind power and photostatic chosen for this study are shown in Figures4, 5, and 6.Figure 4
Renewable energy (wind power) scenarios for a period of one day.Figure 5
Renewable energy (photovoltaics) scenarios for a period of one day.Figure 6
The probability of renewable energy scenarios for a period of one day.


### 4.2. The Base Case of the Dynamic Operation Management of Renewable Microgrid

In this part, its assumed that the renewable microgrid is disconnected from the main power grid. The active power of load is supplied by the microgenerators such as microturbines, diesel generators, fuel cells, wind turbines and photovoltaics, and a battery storage system. The optimal results of the dynamic operation management of renewable microgrid are given by Figures7 and 8.Figure 7
Optimal results of economic operation optimization of renewable microgrid.Figure 8
Optimal results of battery of renewable microgrid.Figures7 and 8 show that microturbine and fuel cell have to generate more electricity power than the diesel generators. Because the diesel generators are expensive, they are restricted to their minimum value during most hours.


### 4.3. The Case of the Renewable Microgrid Operates in Parallel to the Main Power Grid

In the second case, its assumed that the renewable microgrid can operate in parallel to the main power grid. When the renewable microgrid is connected to the main power grid, its loads receive power from both the main power gird and the local distributed energy resources such as solar and wind power, distributed energy generators, and a battery storage system. The optimal results of the dynamic operation management of renewable microgrid are shown in Figures9 and 10.Figure 9
Optimal results of economic operation optimization of renewable microgrid (the main power grid).Figure 10
Optimal results of battery of renewable microgrid (the main power grid).It should be noted that, in Figures9 and 10, the main power grid has to generate more electricity power. Fuel cell and diesel generators are more expensive.In Table5, it is worthwhile to note that the total operation cost of second case is the lowest. In other word, the renewable microgrid which operates in parallel to the main power grid decreases the operation cost because of its loads receiving power from both the main power gird and the local distributed energy resources.Table 5
The comparison of the results obtained by different modes for smart microgrid.
Modes The base case The second case Cost ($) 1.8636e+03 1.83677e+03


### 4.4. The Case of Dynamic Operation Management of Renewable Microgrid under Uncertainties

In this case, a methodology based on scenarios analysis is used to proceed the uncertainty of the power output of distributed renewable generators, such as wind power and photovoltaics power and it is assumed that the renewable microgrid is disconnected from the main grid. The optimal results of the dynamic economic operation optimization of renewable microgrid are shown in Figures11 and 12.Figure 11
Optimal results of economic operation optimization of renewable microgrid under uncertainty.Figure 12
Optimal results of battery of renewable microgrid under uncertainty.In Table6, it is worthwhile to note that this scheme is a direct method to obtain the optimal solution without probability. And the total operation cost of the method of proposed by this paper is the lowest. In other words, modeling the renewable microgrid system under uncertainty decreases the operation cost. This is because of considering different scenarios in the stochastic model instead of single scenario in the deterministic scheme.Table 6
The comparison of the results obtained by dynamic power dispatching under uncertainty.
Dynamic power dispatching strategies under uncertainty The probability of scenario The cost of scenario The cost of scenario 1 0.003 1866.060  The cost of scenario 2 3.542121e-05 1866.166  The cost of scenario 3 3.517950e-09 1866.154  The cost of scenario 4 9.450705e-04 1866.152  The cost of scenario 5 0.232 1866.129  The cost of scenario 6 0.103 1866.166  The cost of scenario 7 0.096 1866.143  The cost of scenario 8 9.103001e-07 1866.140  The cost of scenario 9 0.437 1866.152  The cost of scenario 10 0.127 1866.165  The method of proposed by this paper - 1733.055


## 4.1. Scenario Generation and Reduction Results

The approach used for generation of scenarios may have an impact on the determination of operational cost of renewable microgrid. In the study case, the dates from wind power profiles of the 24 hours are used as scenarios in the dynamic operation management of renewable microgrid. The scenarios of wind power and photostatic chosen for this study are shown in Figures4, 5, and 6.Figure 4
Renewable energy (wind power) scenarios for a period of one day.Figure 5
Renewable energy (photovoltaics) scenarios for a period of one day.Figure 6
The probability of renewable energy scenarios for a period of one day.


## 4.2. The Base Case of the Dynamic Operation Management of Renewable Microgrid

In this part, its assumed that the renewable microgrid is disconnected from the main power grid. The active power of load is supplied by the microgenerators such as microturbines, diesel generators, fuel cells, wind turbines and photovoltaics, and a battery storage system. The optimal results of the dynamic operation management of renewable microgrid are given by Figures7 and 8.Figure 7
Optimal results of economic operation optimization of renewable microgrid.Figure 8
Optimal results of battery of renewable microgrid.Figures7 and 8 show that microturbine and fuel cell have to generate more electricity power than the diesel generators. Because the diesel generators are expensive, they are restricted to their minimum value during most hours.


## 4.3. The Case of the Renewable Microgrid Operates in Parallel to the Main Power Grid

In the second case, its assumed that the renewable microgrid can operate in parallel to the main power grid. When the renewable microgrid is connected to the main power grid, its loads receive power from both the main power gird and the local distributed energy resources such as solar and wind power, distributed energy generators, and a battery storage system. The optimal results of the dynamic operation management of renewable microgrid are shown in Figures9 and 10.Figure 9
Optimal results of economic operation optimization of renewable microgrid (the main power grid).Figure 10
Optimal results of battery of renewable microgrid (the main power grid).It should be noted that, in Figures9 and 10, the main power grid has to generate more electricity power. Fuel cell and diesel generators are more expensive.In Table5, it is worthwhile to note that the total operation cost of second case is the lowest. In other word, the renewable microgrid which operates in parallel to the main power grid decreases the operation cost because of its loads receiving power from both the main power gird and the local distributed energy resources.Table 5
The comparison of the results obtained by different modes for smart microgrid.
Modes The base case The second case Cost ($) 1.8636e+03 1.83677e+03


## 4.4. The Case of Dynamic Operation Management of Renewable Microgrid under Uncertainties

In this case, a methodology based on scenarios analysis is used to proceed the uncertainty of the power output of distributed renewable generators, such as wind power and photovoltaics power and it is assumed that the renewable microgrid is disconnected from the main grid. The optimal results of the dynamic economic operation optimization of renewable microgrid are shown in Figures11 and 12.Figure 11
Optimal results of economic operation optimization of renewable microgrid under uncertainty.Figure 12
Optimal results of battery of renewable microgrid under uncertainty.In Table6, it is worthwhile to note that this scheme is a direct method to obtain the optimal solution without probability. And the total operation cost of the method of proposed by this paper is the lowest. In other words, modeling the renewable microgrid system under uncertainty decreases the operation cost. This is because of considering different scenarios in the stochastic model instead of single scenario in the deterministic scheme.Table 6
The comparison of the results obtained by dynamic power dispatching under uncertainty.
Dynamic power dispatching strategies under uncertainty The probability of scenario The cost of scenario The cost of scenario 1 0.003 1866.060  The cost of scenario 2 3.542121e-05 1866.166  The cost of scenario 3 3.517950e-09 1866.154  The cost of scenario 4 9.450705e-04 1866.152  The cost of scenario 5 0.232 1866.129  The cost of scenario 6 0.103 1866.166  The cost of scenario 7 0.096 1866.143  The cost of scenario 8 9.103001e-07 1866.140  The cost of scenario 9 0.437 1866.152  The cost of scenario 10 0.127 1866.165  The method of proposed by this paper - 1733.055


## 5. Conclusions

In this paper, we introduce a novel dynamic operation management of renewable microgrid. First, to properly model the system, a mathematical model including power losses of the renewable microgrid is established. Then, the uncertainty in the distributed renewable generators, such as wind power and photovoltaics power, are modeled by choosing a set of possible scenarios. And each scenario is assigned a probability that reflects the scenario to be occurred. Finally, the problem of economic operation of renewable microgrid is solved by a feedback control scheme based on the dynamic programming. To solved the deterministic model efficiently, we use the Lagrange multipliers and the interior point method to solve the equality and inequality operation constraints. The numerical results indicate the validity of this optimization methodology.



---

*Source: 5852309-2018-09-20.xml*