Cost Savings Oriented Microgrid Control Strategy Considering Battery Degradation
Cost Savings Oriented Microgrid Control Strategy Considering Battery Degradation
Cost Savings Oriented Microgrid Control Strategy Considering Battery Degradation
Abstract—In this paper a cost savings oriented centralised employing appropriate energy management system (EMS) [5],
control strategy for grid connected hybrid microgrid consisting [6]. Apart from power balancing, the other important aspects
of renewable energy sources and energy storage systems is of energy management system according to latest research are
proposed. In this work, cost savings are achieved by a power
management algorithm (PMA) which takes battery degradation to operate batteries efficiently and to minimize microgrids
cost and grid pricing into account. In-addition the proposed operating costs such as fuel costs, maintenance costs, and
PMA always operates battery within SOC limits and also ensures purchased cost of electricity from the main grid [7].
microgrid reliability. The cost savings are computed as the Grid-connected microgrids can either send/receive power
difference between battery degradation cost and electricity cost. from the grid. This electric power trading with the grid has tra-
However, the battery degradation cost is calculated using depth
of discharge (DOD) which is computed on a real time basis. In ditionally been based on a fixed tariff. However, with the incor-
this work, hybrid energy storage system (HESS) is introduced to poration of smart meters, time-varying tariffs are introduced.
absorb the sudden transients and oscillations in battery current Due to rapid development in ESS technologies and power
thereby increasing battery performance and life span. Apart from electronics, these days microgrids are equipped with wide
cost savings analysis, power quality aspects such as harmonic variety of RES powered DG units and ESSs, which allows
mitigation, reactive power support and unity power factor are
achieved. The proposed control strategy is simulated using a the consumers to sell power during surplus generation/peak
MATLAB based simulink environment during both excess and grid pricing periods, and buy power in case of generation
deficit modes of microgrid by taking PV variations into account. scarcity/off-peak grid pricing periods. To maximize the cost
savings [8] an effective EMS is required which includes the
Index Terms—Battery, supercapacitor, grid side converter, electricity price of the grid and the battery degradation cost.
state of charge, power management, cost savings, battery degra- In [9], a supervisory power management system for a
dation cost, depth of discharge, microgrid. grid interactive microgrid with HESS is proposed and the
objectives of PMA formulated is to minimize/maximize power
I. I NTRODUCTION drawn/supplied from/to grid during peak pricing. In [10], a
Renewable energy sources such as PV and wind are pen- grid adaptive power management strategy for an integrated
etrating more into present day power system in the name microgrid with HESS is proposed and the objectives of this
of microgrids, which are considered as the preferred way control strategy are grid frequency based scheduling of HESS.
to promote integration of renewable energy sources to grid. In [11], a low complexity energy management for DC-coupled
Production of electric power from renewable energy sources hybrid microgrid is proposed and the features of the proposed
results in volatile and intermittent power, to address these energy management scheme are use of sharing coefficient
problems energy storage systems (ESSs) are incorporated in for effective power sharing between energy storages and
microgrids [1], [2]. The primary application of ESSs consisting grid. In [12] a grid connected hybrid microgrid system with
of different energy storage devices (ESDs) is to coordinate dual voltage source inverter is proposed with reduced control
with RES to ensure microgrid generation adequacy. These complexity and overall system cost. From the above literature
ESDs have to be operated within upper and lower SOC [9]–[12], it is clear that none of the EMS/PMA proposed did
limits otherwise, it results in rising of internal impedance not consider battery degradation.
and decomposing of electrolyte in the battery which leads to Moreover, in [13] an EMS is proposed to minimize the
capacity loss and power fade [3]. global energy cost in a microgrid, where the global cost
Microgrids in which both AC and DC voltages are used, includes the cost of energy purchased from grid, energy
called as Hybrid AC/DC Microgrids and these are often generation cost, and the negative cost obtained by selling
implemented [4]. The microgrid can operate both in grid energy to the grid. In [14] an optimal energy management
connected and in islanded mode. The power balancing among strategy for the operation of multiple energy storage units
different sources and sinks is one of the most important (batteries) in grid-connected industrial microgrids is presented.
tasks of microgrid management, which can be achieved by The approach is based on a regrouping particle swarm opti-
,&5(5$
,(((
WK,QWHUQDWLRQDO&RQIHUHQFHRQ5HQHZDEOH(QHUJ\5HVHDUFKDQG$SSOLFDWLRQV 3DULV)5$1&(2FW
Pin(t) Pout(t)
linearAC Loads
BC
DC Loads
Pdc(t)
GRID
DC Link
HESS
PV
Supercapacitor Battery
mization. In [15] a control strategy is proposed to minimize impedances and the ac bus consists of grid and unbalanced
the industrial facility electricity bill costs and this paper deals non-linear ac loads.
with the energy management of battery storage systems in the
frame of real time pricing. From the above literature [13]–[15], III. C OST A NALYSIS
it is clear that none of the papers has considered both battery
degradation cost as well as battery SOC simultaneously for Batteries are one of the most common forms of energy
cost savings. Moreover, in the above works cost savings are storage used in microgrids. In general the battery life cycle
computed by considering constant DOD which is impractical is defined as the charging/discharging cycles performed by
because the batteries are always cycled at various levels of a battery before its capacity fades to 80% of the original
DOD in real-time operation. In-addition, these works do not capacity, the depth of discharge (DOD) is a key variable
consider a combined energy storage system such as HESS that relates to the battery life cycle. As batteries have limited
which will especially improve the transient response of battery cycle life under different DOD scenarios, the battery capacity
as well as DC link voltage control is also not explored. degradation caused by charging and discharging should be a
Considering the above challenging aspects, in this work constituent part of the energy management algorithms. In this
authors have proposed a cost savings oriented centralized work, cost savings are calculated by considering electricity
EMS for grid connected hybrid microgrids considering bat- purchase costs/selling income from/to the grid and battery
tery degradation cost, grid electricity cost and battery SOC. degradation cost. Battery degradation cost is of three types.
Moreover, the proposed energy management system considers However, degradation cost due to SOC and temperature are
the following vital aspects: a) supercapacitor is introduced to not considered in this work and only DOD related degradation
absorb the sudden transients and oscillations in battery current cost is considered.
thereby increasing battery performance and life span [16], The equivalent full cycles are used for assessing the actual
b) compensation of unbalanced and non-linear load [17] is battery cycle life. But to evaluate the degradation cost of a
achieved to improve the power quality, and c) better dc link single time charge and discharge cycle on a real time basis,
voltage regulation. in this work DOD is determined for a duration of every one
hour. To determine DOD, initial and final SOCs of battery are
II. S YSTEM D ESCRIPTION required. If the initial SOC i.e., SOCi of the battery is known
the final SOC i.e., SOCf at the end of every one hour can
The proposed hybrid microgrid structure considered for the be calculated using coulomb counting method as given below
study is shown in Fig. 1. Renewable energy sources such as [3].
PV modules are interfaced to common dc bus (dc link) using
a high gain converter (HGC). Due to the intermittent nature 1
SOCf = SOCi − ib dt. (1)
of renewable energy sources, storage systems have become 3600Cb
an integral part of the microgrid. In this work, a combination
Where ib , SOCi , and Cb are current, initial state of charge
of battery and supercapacitor called as hybrid energy storage
and nominal capacity of battery respectively. After computing
system (HESS) is considered and it is interfaced to common
the final SOC of battery, the battery degradation cost for a
dc bus using a bidirectional converter (BC), which enables
particular duration i is calculated using the (2), where the
power flow into and from the energy storage elements. Various
factor 2 in the denominator represents both charging and
topologies of grid side converters are available for interfacing
discharging energy throughputs are considered in the cost [18],
ac and dc sides. However, a neutral point clamped (NPC),
[19].
voltage source converter (VSC) topology is used for both
bidirectional power flow and for load compensation. The DOD(i)
grid side converter (GSC) is interfaced to ac bus using filter BDODcost (i) = Cbp (2)
2Ntotal (i)
,&5(5$
WK,QWHUQDWLRQDO&RQIHUHQFHRQ5HQHZDEOH(QHUJ\5HVHDUFKDQG$SSOLFDWLRQV 3DULV)5$1&(2FW
18000 0.3
where, pavg (t) is the average power difference between PV
16000 0.26
and ac, dc loads. The average power difference has to be
12000 0.2
10000 0.16 tions in the dc link voltage. The transient power ptr (t) and the
8000 0.13 oscillatory power pos (t) are provided by the supercapacitor.
6000 0.1
4000 0.06 B. Reference current generation
2000 0.03
0 0 The reference currents for various power converters can be
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Depth of discharge derived from the following effective current (ief (t)) equation
using PI controller as given below
Fig. 2. DOD vs Battery lifetime vs Degradation cost.
ief (t) = Kp ve + Ki ve dt (11)
Where DOD = |SOCi − SOCf |, Ntotal is number of cycles = iavg (t) + itr (t) + ios (t) (12)
and Cbp is battery purchase price. Fig. 2 shows the variations
in battery lifetime (number of cycles) and battery degradation where, ve = vdcref − vdc is the difference between reference
cost with respective to DOD where the battery purchase price and actual dc link voltage, which arises due to input and
is taken as 500 $ by considering 10 packs of 12 V, 26 Ahr output power difference in the dc link. In (11) Kp and Ki are
lead acid EXIDE batteries. The purchase cost/selling income proportional and integral gains of PI controller. The average
from/to grid is calculated by considering two types of tariffs component of effective current is obtained by passing through
namely time of use tariff and feed in tariff. The electricity a low pass filter whose cut-off frequency is chosen to be 5
purchase cost during time period i can be given as Hz.
1
Cp (i) = Pμgrequired (i).λ(i) (3) iavg (t) = ief (t) (13)
(1 + sτ )
Where Pμgrequired and λ are power required in microgrid and The average component of current has to be deliv-
unit electricity price in time period i. Similarly, the income ered/absorbed by grid/battery. Similarly, the transient and
from selling electricity to the grid during time period i can be oscillating components supplied by supercapacitor are given
written as follows below.
Cs (i) = Pμgavailable (i).ϑ(i) (4) iosc (t) = ios (t) + itr (t) = ief (t) − iavg (t) (14)
Where Pμgvailable and ϑ are power available in microgrid and
C. Grid Side Converter Control
feed in tariff during time period i. Considering (2), (3), and
(4) the total cost over a time period i can be given as Instantaneous symmetrical component theory (ISCT) [20]
is one of the time domain control theories used to generate
Ct (i) = BDODcost (i) + Cp (i) − Cs (i) (5) reference currents for active power filters, which are used
IV. DC L INK C ONTROL , R EFERENCE C URRENT for power quality improvement. According to this theory, the
G ENERATION , & GSC C ONTROL compensator/filter reference currents are as given below.
A. DC Link Control ∗ vg(abc) + β(vg(bca) − vg(cab) )
if (abc) = il(abc) − 2 Plavg (15)
A dc storage capacitor is used for implementing dc link and i=a,b,c vgi
the instantaneous power flow relationship across the dc link
+
)
can be given as where β, φ+ , i∗f (abc) , il(abc) , vg(abc) , and Plavg are tan(φ
√
3
,
power factor angle between voltage and current positive se-
pin (t) = pdc (t) + pout (t) (6)
quence components, filter reference currents, load currents,
pin (t) = ppv (t) ± pb (t) ± psc (t) − pdcl (t) (7) grid voltages, and average load power on ac side. For unity
pout (t) = ±pg (t) − pacl (t) (8) power factor operation β = 0. Therefore, the above equation
is written as given below.
where, pacl (t), pg (t), pdcl (t), psc (t), pb (t), and ppv (t) are
instantaneous ac load, grid, dc load, supercapacitor, battery, vg(abc)
∗
and PV powers. The bidirectional power flow is denoted by if (abc) = il(abc) − 2 Plavg (16)
i=a,b,c vgi
the notation ± . The dc link power is divided into three
components i.e., average power component pavg (t), transient To incorporate bidirectional power flow feature of grid side
power component ptr (t), and oscillatory power component converter in a microgrid environment (16) is modified as
pos (t) as shown below. follows [21].
pdc (t) = pin (t) − pout (t) (9) ∗ vg(abc)
if (abc) = il(abc) − 2 (Plavg − Pμg )(17)
= pavg (t) + ptr (t) + pos (t) (10) i=a,b,c vgi
,&5(5$
WK,QWHUQDWLRQDO&RQIHUHQFHRQ5HQHZDEOH(QHUJ\5HVHDUFKDQG$SSOLFDWLRQV 3DULV)5$1&(2FW
where Pμg is the net power available in the microgrid and it is Start
the difference between source and sink powers in microgrid.
Read iavg(t), iosc(t),vsa, il(abc), SOCb, SOCbl, SOCbu, immp
Hence, using this formulation of ISCT (17) for microgrid
applications requires measurement as well as processing of
all source and sink power signals in a microgrid. This results Microgrid Deficit Mode
iavg(t) 0
Microgrid Excess Mode
N Y
in increased sensor requirement, increased control complexity
Compute Energy Compute Energy
and high communication bandwidth for centralized control. Required in µg Available in µg
ib = ibrated, isc = iosc, ipv = immp, SOCbl = 0.5
To address the above-mentioned issues, authors have pro- ifj(abc) = il(abc) -(vsa/Ʃi = abc v2si)(1.5Vmiavg)
posed an improved version of (17) for microgrid applications N SOCb < Y
as given below. SOCbu
SOCb <
SOCbl
Y Peak
N
Pricing
∗ vg(abc)
if (abc) = il(abc) − 2 (1.5iavg Vm ) (18) Compute Battery
i=a,b,c vgi
N Y Degradation Cost
ib = 0, isc = iosc, ipv = immp, SOCbl = 0.2
For regulating DC link voltage, the effective current average Compute Battery ifj(abc) = il(abc) -(vsa/Ʃi = abc v2si)(1.5Vmiavg)
Y BDODcost
Degradation Cost
component has to be delivered/absorbed by grid/battery. But < Cs
using (18) the effective current average component is deliv- ib = -ibrated, isc = iosc, ipv = immp, SOCbu = 0.5
ifj(abc) = il(abc) -(vsa/Ʃi = abc v2si)(1.5Vmiavg) N
ered/ absorbed by grid only. However, in some cases like when
grid is in peak pricing mode battery has to supply the average BDODcost N
< Cp ib = 0, isc = iosc, ipv = immp, SOCbu = 0.9
component of effective current and grid current should be zero. ifj(abc) = il(abc) -(vsa/Ʃi = abc v2si)(1.5Vmiavg)
To incorporate this feature (18) is modified as given below. Y
Y ib = iavg, isc = iosc, ipv = immp, SOCbl = 0.2
ifj(abc) = il(abc) Peak
∗ vg(abc)
if (abc) = il(abc) − 2 (1.5Kiavg Vm )(19) ib = iavg, isc = iosc, ipv = immp, SOCbu = 0.6
Pricing
i=a,b,c vgi ifj(abc) = il(abc) N
,&5(5$
WK,QWHUQDWLRQDO&RQIHUHQFHRQ5HQHZDEOH(QHUJ\5HVHDUFKDQG$SSOLFDWLRQV 3DULV)5$1&(2FW
grid will take the entire excess power. In this case there will 8 Ipv
250
be significant amount of savings as power is fed back to grid. 6 Defict Mode Defict Mode
240
Current (A)
Voltage (V)
Subcase-I: BDODcost < Cs : Here, rather than keeping bat- 4 Excess 230 Vdc
Mode
tery idle, it is discharged till SOCb = 0.5 to maximize savings. 2
220
grid will take the entire excess power available in microgrid. 0.4 Pg Pacl
Current (A)
Current (A)
2
5
Subcase-II: Peak Pricing Disabled: Here, grid is in off-peak 0 0
pricing mode which means, the income earned by supplying -2 Isc -5
Current (A)
VI. S IMULATION R ESULTS 0.3 2
Battery SOC
0
The proposed system is simulated using MATLAB based 0.2
-2
simulink environment. The system parameters considered for 0.1 (g)
-4
(h)
simulation studies are mentioned in Table I. Tests are con- 0 4 8 12 16 Time (hrs) 24 4 Time (hrs)
ducted to verify the steady state and transient state perfor- 80 Vsa
mance of the proposed system and the detailed simulation 5 40
Current (A)
Current (A)
1
results are as follows. 0 0
The variations in PV current for a duration of 24 hrs is -40
-5
shown in Fig. 4(a), where the PV output increases till midday
(i) -80 (j) Isa
and decreases thereafter. The variations in dc link voltage are 4 Time (hrs) 9 Time (hrs)
shown in Fig. 4(b), where the dc link voltage settles at rated
value of 220 V. The pricing scheme considered in this work Fig. 4. Simulation results: (a) PV current, (b) dc link voltage, (c) Time
is time-of-use pricing, where the electricity price is set as of use tariff, (d) Various active powers as per labelling, (e) variations in
supercapacitor and battery currents, (f) variations in Grid current, (g) Battery
0.1$/kWh during off-peak hours and 0.25$/kWh at peak hours SOC, (h) Unbalanced and non-linear load currents, (i) Zoomed view of GSC
as shown in Fig. 4(c) and the selling price is considered as currents, (j) unity power factor (UPF) operation.
80% of the electricity price [23]. The variations in grid, GSC,
AC and DC loads powers are shown in Fig. 4(d). The battery
and supercapacitor currents are shown in Fig. 4(e), where the
supercapacitor absorbs the sudden transients and oscillations
TABLE I thereby smoothening the battery current.
S YSTEM PARAMETERS
The grid current variations for the entire duration of op-
Supercapacitor specifications Values eration are shown in Fig. 4(f), after t = 13 hrs grid currents
Terminal voltage (Vsc ) 16 V become zero since battery supplies the deficit power . Battery
No. of packs in series 11
Battery specifications Values SOC is shown in Fig. 4(g), where battery is charged till t =
Ah Capacity 26 Ah 9 hrs from grid. Battery kept idle from t = 9-13 hrs, since
Terminal voltage (VB ) 12 V the grid is in peak pricing mode. From t = 13-17, battery
No. of batteries in series 10
DC load parameters Rdcl =320 Ω
is discharged to supply microgrid deficit power as the battery
AC load parameters Values degradation cost is less than the price required to supply deficit
Unbalanced linear load Ra = 20 Ω, La = 65 mH power from grid and from t = 17-24, battery is charged as the
Rb = 10 Ω, Lb = 50 mH grid is in off-peak mode. The zoomed version of unbalanced,
Rc = 30 Ω, Lc = 80 mH
Nonlinear load 3-Φ, diode bridge rectifier with nonlinear load currents are shown in Fig. 4(h). The zoomed
Rl = 100 Ω, version of grid currents shown in Fig. 4(i) are balanced and
Grid side converter parameters Values sinusoidal, the in phase operation of grid currents with grid
DC link voltage Vdc =220 V
DC link capacitance Cdc = 1650 μF voltages are shown in Fig. 4(j) thus demonstrating unbalanced
Hysteresis band h = ±0.1A and nonlinear load compensation.
Grid/source voltages 3-Φ, 50 V (P-N) The microgrid operational cost analysis is shown in table-II,
,&5(5$
WK,QWHUQDWLRQDO&RQIHUHQFHRQ5HQHZDEOH(QHUJ\5HVHDUFKDQG$SSOLFDWLRQV 3DULV)5$1&(2FW
,&5(5$