Renewable Energy: S. Rajanna, R.P. Saini
Renewable Energy: S. Rajanna, R.P. Saini
Renewable Energy: S. Rajanna, R.P. Saini
Renewable Energy
journal homepage: www.elsevier.com/locate/renene
a r t i c l e i n f o
a b s t r a c t
Article history:
Received 1 May 2016
Received in revised form
30 June 2016
Accepted 8 August 2016
The present study focuses on optimal sizing of an integrated renewable energy (IRE) system with battery
bank to meet the load demand of a cluster of four zones in Karnataka state of India. Hourly load demand
has been simulated using HOMER software based on three scenarios - high investment with low (HILR)
rating appliances, medium investment with moderate (MIMR) rating appliances, and low investment
with high (LIHR) rating appliances. Simulation results of optimal size, total net present cost (TNPC) and
cost of energy (COE) have been obtained using genetic algorithm (GA) considering different scenarios
with and without DSM strategy by keeping energy index at 1. Results of three scenarios with DSM
strategy were compared with those without DSM strategy. Based on the comparison of the results,
scenario of HILR with DSM strategy is found to have minimum TNPC and COE for zone 1 and zone 2,
while scenario of MIMR with DSM for zone 3 and LIHR with DSM for zone 4 are found to have least TNPC
and COE. Further, obtained results of TNPC of GA have been validated by comparing PSO technique. It has
been observed that, results of GA and PSO were found to be almost same.
2016 Elsevier Ltd. All rights reserved.
Keywords:
GA
IRES
Demand side management
Net present cost
Remote rural area
1. Introduction
Energy is generally accepted as a vital component to a nation's
growth and development of its citizens. Power generation across
the world is currently dependent to a large extent, on conventional
sources of energy such as coal, gas, oil, etc. However, the use of such
power sources has a damaging effect on environment and society,
and gives rise to various environmental and health problems.
Renewable energy resources such as biomass, micro hydro, wind,
solar, etc. offer a way of generating power with considerably less
damage to the environment. Various social and economic benets
also accompany by generating power from renewable energy resources [1e4]. It is noteworthy that economical and efcient utilization of renewable energy resources warrants development of
optimal models through which demand for energy can be met.
Further, IRE system modeling is a complex process that involves
formulation of mathematical models pertaining to each component
[5e7]. Efcient IRE system modeling depends on effective load
demand management. Demand-side management (DSM) offers
* Corresponding author.
E-mail address: srajannamce@gmail.com (S. Rajanna).
http://dx.doi.org/10.1016/j.renene.2016.08.024
0960-1481/ 2016 Elsevier Ltd. All rights reserved.
1162
Nomenclature
AC
ACnpv
ASPV
Aw
BGG
BMG
BATT/BS
BC
CFL
COE
CN
alternating current
net present value of appliances cost ($)
area of SPV (m2)
swept area of wind turbine rotor (m2)
biogas generator
biomass generator
battery system
bi-directional converter
compact uorescent
cost of energy ($/kWh)
runoff curve number for hydrological sail (CN 40e58
for dense forest)
Cb
initial capacity of battery system
CD
cow dung (ton/yr)
CP
power co-efcient (0.45)
COA
cost of appliances
CNPV
net present capital cost ($)
CRF
capital recovery factor
CVBG
caloric value of biogas (4700 kcal/kg)
CVBG
caloric value of biomass (4015 kcal/kg)
D
duration of unavailability of load (hour)
DC
direct current
DOD
depth of discharge
DSM
Demand side management
Egen
total annual generation (kWh)
EENS
expected energy not supplied (kWh)
EIR
energy index ratio
ESL
shiftable load (kWh)
ENSL
non-shiftable load (kWh)
FCNPV
net present fuel cost ($/ton)
g
acceleration due to gravity (m/s2)
GA
Genetic algorithm
hnet
net head (m)
HILR
high investment with low rating
HT
solar radiation (kW/m2)
I
monthly rainfall (mm)
IL
incandescent lamp
IRE
integrated renewable energy
L
average annual power load (kW)
LED
Light emitting diode
LIHR
low investment with high rating
Lw
length of the watershed (m)
MHG
micro hydro generator
MIMR
medium investment with moderate rating
NSPV
solar panels (Nos.)
NWTG
wind turbine (Nos.)
NBS
battery system (Nos.)
NREL
national renewable energy limited
OM
Operation and maintenance
OM npv Net present value of O&M ($)
PBMG(t) power output of BMG (kWh)
PBGG(t)
PWTG(t)
PSPV(t)
Pr
POA
PSO
Q
Qd
RNPV
S
SONPT
SOFFPT
SOC
Tc
Tp
TNOA
TNPC
Vr
Vci
Vco
WTG
hMHP
hBGG
hBMG
hSPV
hWTG
hINV
hRECT
hBC
hBD
s
aBMG
aBGG
aMHG
aSPV
aWTG
aBS
aCONV
bBMG
bBGG
bMHG
bSPV
bWTG
bBS
bCS
rW
ra
5
m
lAFC
t
$
[19] proposed IRE model for Uttarakhand state of India. Ramakumar et al. [20] proposed a linear program based IRE model
developed for developing country.
Randa Kallel et al. [21] highlighted advantages of the suggested
strategy for residential hybrid system functions under different
conditions and conducted a comparative study between energy
management strategies with and without DSM strategy. Rajanna S
and Saini RP [22] proposed a genetic algorithm based optimal IRE
1163
investment scenarios without and with DSM strategy is emphasized to develop optimal IRE models. Further, present study also
presents for a comparative analysis between the two strategies
investigated.
Table 1
The general information of the study area.
Features
Zone 1
Zone 2
Zone 3
Zone 4
11
5
302
898
12.05
77 400 4800
23
6
395
1141
11 520 58.800
77 160 58.800
24
8
763
2718
11 580 58.800
77 70 58.800
10
7
780
1906
11 450 57.600
77 50 6000
Table 2
Technology key with three investments based appliances for residential, community and commercial sectors.
Scenarios
Residential
Community
Commercial
Lighting
Entrainment TV
Cooling fan
Lighting
Street lighting
Cooling fan
Refrigerator
Lighting
HILR
MIMR
LIHR
LED
CFL
IL
TV-3
TV-2
TV-1
R3
R2
R1
LED
CFL
IL
ST3
ST2
ST1
R3
R2
R1
Tier-3
Tier-2
Tier -1
LED
CFL
IL
Technology key and investment for residential, community and commercial sectors
Appliances
Prices ($)
Appliances
Prices ($)
0.38
5.71
7.04
22.53
28.38
29.73
84.10
112.63
127.65
3.00
11.26
22.53
300.35
435.50
585.67
1164
medium investment appliances - CFL, LCD TV, CFL street light, high
efciency refrigerator (Tier-2), BEE 4 star ceiling fan (R2); and low
investment appliances - incandescent bulb, regular TV, uorescent
lamp, regular refrigerator (Tier-1), BEE 3 star ceiling fan (R1).
Detailed description regarding technology key, ratings and electrical appliances prices is given in Table 2.
Hourly load demand has been calculated as hourly segments in a
day based on operating hours and considered investment scenarios
as given in Table 2. An hourly load has been estimated for different
investment scenarios of LIHR, MIMR and HILR and hourly load
variation of three investment scenarios is shown in Fig. 6(aec). The
load assessment for all four zones has been done on the basis of
four different load sectors (domestic, community, commercial and
small industrial). Each load sector considers different electrical
appliances in three distinct scenarios. For each zone, total load
demand has been computed based on power ratings and operating
hours of appliances in a day. Detailed description of load sectors,
considered scenarios and total electrical appliances for all four
zones are given in Table 3.
The data have been synthesized through HOMER software
considering estimated values of three investment scenarios based
hourly load with and without DSM are given in Tables 4 and 5. The
total estimated annual loads of zone 1, zone 2, zone 3 and zone 4 are
shown in Fig. 1.
2.3. Resources assessment
Mathematical models of system components of the IRE system
are essential before optimizing size because they give information
about the performance of system components in different operating conditions. The mathematical modeling of solar system, micro hydro generator, wind turbine generator, biogas generator,
biomass generator and battery storage system are developed based
on the available energy potential in the study area are discussed as
follows:
2.3.1. Micro hydro generator (MHG)
Micro hydro generator is designed considering low volumetric
ow rate and low head conditions. The potential of micro hydro
power (PMHG) at a particular site area depends on net head and
dependability of the ow rate (Q). The discharge was estimated for
the considered zones on the basis of topography and monthly
average rain fall data using curve number method.
Qd
I 0:2S2
I 0:8S
25400
254
S
CN
The peak rate of runoff or discharge of a watershed is calculated
by using following Equations [28,29];
p
TP 0:6TC TC
r
0:77
.
TC 0:0195 Lw 3 Hnet
(2)
0:0208 A Qd
TP
(3)
(4)
Table 3
Categorization of load sectors with power rating of considered investment scenario for study area.
Load sectors
Parameters
Domestic
Household
Community
Hospital
Primary School
Community
Community hall
Commercial
Industrial
Appliances
Lighting
Ceiling fan
TV
Radio
Refrigerator
Lighting
Ceiling fan
Lighting
Ceiling fan
pumping water
street Lighting
Lighting
Ceiling fan
Lighting
Flour mill
Saw mills
2 point
1 point
1 point
1 point
1 point
3 point
3 point
1 2point
12 point
1 motor (5HP) set at 300 HH
1pole at cluster 6 HH
1 point
1 point
1 shop at cluster of 125 HH
1 our mill at cluster of 250
1 saw mill at cluster of 250 HH
Power rating
Nos. used
HILR
MIMR
LIHR
Zone 1
Zone 2
Zone 3
Zone 4
9
52
40
10
417
9
52
9
52
3675
20
9
52
9
3750
5000
18
55
90
10
434
18
55
18
55
3675
36
18
55
18
3750
5000
30
56
140
10
475
30
56
30
56
3675
60
30
56
30
3750
5000
604
302
302
302
1
3
3
12
12
3
30
1
1
5
3
1
790
395
395
395
1
3
3
12
12
3
39
1
1
6
3
2
1526
763
763
763
1
3
3
12
12
4
76
1
1
8
4
3
1560
780
780
780
1
3
3
12
12
4
78
1
1
8
4
3
1165
Table 4
Zone wise hourly estimated load prole in different investment scenarios without DSM strategy.
HILR
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
MIMR
LIHR
Zone 1
Zone 2
Zone 3
Zone 4
Zone 1
Zone 2
Zone 3
Zone 4
Zone 1
Zone 2
Zone 3
Zone 4
16.75
16.75
16.75
16.75
19.45
19.45
25.03
25.03
11.29
11.29
12.24
40.57
40.57
40.57
12.24
12.24
12.24
20.35
45.34
45.34
45.344
34.31
16.75
16.75
21.76
21.76
21.76
21.76
25.32
25.32
29.22
29.22
14.64
14.64
15.59
52.64
52.64
52.64
15.59
15.59
15.59
26.49
55.75
55.75
55.75
44.73
21.77
21.76
42.56
42.56
42.56
42.56
49.58
49.58
50.22
50.22
28.5
28.5
29.45
90.65
90.65
90.65
29.45
29.45
29.45
51.9
102.58
102.58
102.58
87.88
42.56
42.56
41.64
41.64
41.64
41.64
48.51
48.51
49.45
49.45
27.89
27.89
28.83
89.35
89.35
89.35
28.83
28.83
28.83
50.78
100.67
100.67
100.67
85.97
41.64
41.64
18.18
18.18
18.18
18.18
23.62
23.62
28.22
28.22
11.76
11.76
12.87
56.3
56.3
56.3
12.87
12.87
12.87
35.92
67.35
67.35
67.35
56.32
18.19
18.18
23.62
23.62
23.62
23.62
30.73
30.73
33.38
33.38
15.25
15.25
16.36
73.16
73.16
73.16
16.36
16.36
16.36
46.85
84.48
84.49
84.49
73.46
23.62
23.62
45.19
45.19
45.19
45.19
58.92
58.92
57.48
57.48
29.04
29.04
30.16
123.83
123.83
123.83
30.16
30.16
30.16
90.08
156.14
156.14
156.14
141.44
45.19
45.19
46.2
46.2
46.2
46.2
60.24
60.24
50.62
58.42
29.68
29.68
30.8
126
126
126
30.8
30.8
30.8
92.08
159.28
159.28
159.29
144.59
46.2
46.2
19.28
19.28
19.28
19.28
28.34
28.34
32.04
32.04
11.95
11.95
13.24
71.77
71.77
71.77
13.24
13.24
13.24
51.21
90.85
90.85
90.85
79.83
19.28
19.28
25.03
25.03
25.03
25.03
36.88
36.88
38.36
38.36
15.49
15.49
16.77
93.32
93.32
93.32
16.77
16.77
16.77
66.84
115.2
115.2
115.2
104.16
25.03
25.03
47.85
47.85
47.85
47.85
70.74
70.74
67.06
67.06
29.47
29.47
30.76
158.83
158.83
158.83
30.76
30.75
30.75
128.66
215.33
215.33
215.33
200.63
47.85
47.85
48.93
48.93
48.93
48.93
72.33
72.33
68.23
68.22
30.12
30.12
31.4
161.85
161.85
161.85
31.4
31.4
31.4
131.52
219.81
219.81
219.815
205.11
48.93
48.93
Table 5
Zone wise hourly estimated load prole in different investment scenarios with DSM strategy.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
HILR
MIMR
LIHR
Zone 1
Zone 2
Zone 3
Zone 4
Zone 1
Zone 2
Zone 3
Zone 4
Zone 1
Zone 2
Zone 3
Zone 4
16.75
16.75
16.75
16.75
19.47
19.47
25.03
25.03
27.54
27.54
28.49
24.32
24.32
24.32
23.26
23.26
23.26
20.35
34.31
34.31
34.31
34.31
16.75
16.75
21.76
21.76
21.76
21.76
25.32
25.32
29.22
29.22
35.89
35.89
36.84
31.39
31.39
31.39
26.61
26.61
26.61
26.49
44.73
44.73
44.73
44.73
21.77
21.77
42.56
42.56
42.56
42.56
49.58
49.58
50.22
50.22
58.5
58.5
59.45
60.65
60.65
60.65
44.15
44.15
44.15
51.9
87.88
87.88
87.88
87.88
42.56
42.56
41.64
41.64
41.64
41.64
48.51
48.51
49.45
49.45
57.89
57.89
58.83
59.35
59.35
59.35
43.53
43.53
43.53
50.78
85.97
85.97
85.97
85.97
41.64
41.64
18.18
18.18
18.18
18.18
23.65
23.61
28.2
28.22
28.01
28.01
29.12
40.05
40.05
40.05
23.9
23.9
23.9
35.92
56.32
56.32
56.32
56.32
18.18
18.18
23.62
23.62
23.627
23.62
30.73
30.73
33.38
33.38
36.5
36.5
37.61
51.91
51.91
51.91
27.39
27.39
27.39
46.85
73.46
73.46
73.46
73.46
23.62
23.62
45.19
45.19
45.19
45.19
58.92
58.92
57.48
57.48
54.04
54.04
55.16
98.83
98.83
98.83
44.86
44.86
44.86
90.08
141.44
141.44
141.44
141.44
45.19
45.19
46.2
46.2
46.2
46.206
60.24
60.24
50.634
58.43
54.68
54.68
55.8
101
101
101
45.5
45.5
45.5
92.08
144.58
144.58
144.58
144.58
46.2
46.2
19.277
19.277
19.277
19.277
28.337
28.337
32.036
32.036
28.201
28.201
29.487
55.517
55.517
55.517
24.262
24.262
24.262
51.211
79.827
79.827
79.827
79.827
19.277
19.277
25.025
25.025
25.025
25.025
36.875
36.875
38.36
38.36
36.74
36.74
38.02
72.07
72.07
72.07
27.80
27.80
27.80
66.84
104.18
104.18
104.18
104.18
25.03
25.03
47.85
47.85
47.85
47.85
70.74
70.74
67.06
67.06
50.72
50.72
52.01
137.58
137.58
137.58
45.46
45.46
45.46
128.66
200.64
200.64
200.64
200.64
47.85
47.85
48.93
48.93
48.93
48.93
72.335
72.33
68.22
68.22
51.37
51.37
52.65
140.60
140.60
140.60
46.10
46.10
46.10
131.52
205.11
205.11
205.12
205.12
48.93
48.93
Annual Load
demand (kWh)
Hours
1050000
550000
50000
HILR MIMR LIHR
Zone 1
Fig. 1. Total annual load demand considering three investment scenarios for all zones.
Appliance cost
(millions)
1166
10
8
6
4
2
0
HILR MIMR LIHR
zone 1
zone 2
zone 3
zone 4
Total gas yield m3 =day CVBG hBGS
PBGG t
860 operating hours per day
(5)
(7)
where ASPV is the active area of SPV panel (1.63 m2), HT is the solar
radiation in kW/m2 and hspv denotes the efciency of the SPV
panels (14.4%).
8
<0
PW t Pr
:
Pr
v3 t v3ci
v3r
v3ci
vt vci or vt vco
vci < vt < vr
vr < vt < vco
(8)
1
Pr $Cp $ra $hg $Aw $Vr3
2
(9)
Cmin DOD Cb
(10)
(11)
As maximum energy stored in the battery system cannot surpass maximum state of charge (SOC max), the following condition
must be met during optimization.
(12)
SOC max takes the value of nominal capacity of the battery bank
(Cb).
The discharge process of the SOC is expressed by Eq. (12) [10,38]
as:
(13)
SOCt SOCmin
Fig. 3. Flow duration curve [27,28].
(14)
Further, during discharge of the battery, if SOC (t) < SOC min, a
percentage of energy demand to the load is not met. In this case,
1167
30
20
(a)
10
11-12pm
10-11pm
8-9pm
9-10pm
7-8pm
6-7pm
5-6pm
4-5pm
3-4pm
2-3pm
12pm
121pm
1112am
910am
1011am
78am
89am
56am
67am
45am
34am
12am
23am
0
0-1am
Hour (h)
Out put power of BGG in MIMR scenario
50
40
(b)
30
20
10
10-11pm
7-8pm
7-8pm
11-12pm
6-7pm
6-7pm
9-10pm
5-6pm
5-6pm
8-9pm
3-4pm
4-5pm
4-5pm
2-3pm
2-3pm
3-4pm
12pm
12pm
121pm
1112am
1011am
89am
910am
78am
67am
56am
34am
45am
23am
0-1am
12am
Hour (h)
60
50
40
(c)
30
20
10
11-12pm
10-11pm
8-9pm
9-10pm
121pm
1112am
1011am
89am
910am
78am
67am
56am
45am
34am
23am
0-1am
12am
Hour(h)
Fig. 4. (aec). Output power of BMG and BGG in different investment scenarios for zone 1.
SOC (t) is set to SOC min and the loss of power supply is considered
as unmet load.
3. Problem formulation
Under the present study, problem formulation involves objective function, reliability and economic constraints, which are
described under this section.
3.1. Objective function
Net present cost (NPC) has been found as a useful indicator as it
includes all cost that occurs within the project lifetime. The total
net present cost (TNPC) comprises initial capital cost, replacement
costs, maintenance costs, fuel costs and investment cost of appliances. Therefore minimization of TNPC is considered as an economic indicator used for the tness function for the proposed
system. It can be expressed as [32e34];
(15)
1168
the lifetime of the project is considered for 20 years while the life of
the battery and inverter system is considered as 4 and 10 years
respectively. Therefore, 3 times replacement for battery and 1 for
inverter have been considered during the projects lifetime. As results, the net present value and number of times within the N year's
horizon that need to replacement of batteries and inverters are
calculated based on the following Eq. (17) [10,38] as:
Xb
N
1
Lb
(17)
X 1 mj
1 m 10
aINV $PINV
1i
1i
j5;10;15
(18)
n
X
1m j
j1
1i
(22)
CRF
g1 gt
1 gt 1
(23)
where Egen is the total annual energy generated by the IRE system,
TNPC is total net present cost and CRF is capacity recovery factor.
EENS
8760
X
L:D
(24)
i1
where, EIR is dened as one minus the ratio of EENS to total energy
demand of the system (EO). It can be expressed as;
EENS
EO
(25)
4. Methodology
(20)
COA POA*TNOA
TNPC CRFg; t
P8760
t1 Egen t
EIR 1
(19)
FCNPV
COE Rs=kWh
(21)
1169
Table 6
Technical and cost parameters of system components.
Parameters
Units
MHP
BGG
BMG
SPV
WTG
BATT
BC
Capital cost
O&M cost
Power rating
Area of SPV panel
Panel efciency
Rated speed
Swept area
Cut-in speed
Cut-out speed
DOD
Voltage
Efciency
Nominal capacity
$/kW
$/yr
kW
m2
%
m/s
m2
m/s
m/s
Volt
%
ah
1201.38
30.03
1
-
750.86
18.77
1
-
901.04
27.03
1
-
590.03
11.80
0.230 Wp
1.63
0.144
-
3075.54
92.27
1
10
4.26
1.5
25
-
187.72
7.51
50
24
-
170.60
3.41
1.5
90
150
AC Bus
DC Bus
DC Renewable
generators
Load
AC Renewable
generators
Charge controoler
Bidirectional
converter
Battery bank
The low cost with high power rating based appliances are
considered in this scenario.
(b) Medium investment scenario
The moderate costs with power rating based appliances are
considered in this case.
(c) High investment scenario
In this scenario, high expensive with low power rating based
appliances are considered.
Table 7
Load type and load identication.
Sl. no.
Load type
1
2
Non-shiftable
Shiftable
1170
Table 8
Peak shifting strategy based DSM strategy.
Time segment
Water pump
0:0e1:0
1:0e2:0
2:0e3:0
3:0e4:0
4:0e5:0
5:0e6:0
6:0e7:0
7:0e8:0
8:0e9:0
9:0e10:0
10:0e11:0
11:0e12:0
12:0e13:0
13:0e14:0
14:0e15:0
15:0e16:0
16:0e17:0
17:0e18:0
18:0e19:0
19:0e20:0
20:0e21:0
21:0e22:0
22:0e23:0
23:0e0:0
Duration of time length
Flourmills
Saw mills
SONPT
ROFFPT
SONPT
ROFFPT
SONPT
ROFFPT
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
5 h/day
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
3 h/day
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
3 h/day
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
3 h/day
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
3 h/day
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
ON
ON
ON
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
OFF
3 h/day
(a)
50
40
30
20
10
0
1
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour (h)
(b)
80
70
60
50
40
30
20
10
0
1
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour (h)
(c)
100
80
60
40
20
0
1
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour (h)
Fig. 6. (aec). Hourly load prole with and without DSM strategy for zone 1; (a) HILR (b) MIMR (c) LIHR.
1171
Resources assessment
DC Renewable energy
potentials
(SPV panels and
WTG)
AC Renewable energy
potentials
(MHG,BMG and BGG)
d=d+1
t= t+1
A
Yes
t=24
Yes
Simulation
Failed
d=d+1
d=365
No
Yes
Simulation successful
Finish
Fig. 7. Methodology for optimizing IRE system using with and without DSM strategy.
1172
ELoadWITHOUT:DSM t ENSL t
n
X
ESL t SONPT t
(26)
i1
Load demand with DSM strategy can be expressed by Eq. [18] as;
ELoadWITH:DSM t ENSL t
n
X
ESL t SOFFPT t
(27)
i1
Finally, hourly annual load demand of the study area with and
without DSM strategy has been simulated using HOMER software
based on low, medium and high investment scenarios, the ow
Without DSM
strategy
No
Dump Load
Yes
SOC < SOC max
Battery charging
Calculaion of new SOC
No
Yes
Yes
P
(t) +P
< P (t)
No
Yes
To
With DSM
strategy
No
Dump Load
Yes
SOC < SOC max
1173
Battery charging
Calculaion of new SOC
No
Yes
Yes
P
(t) +P
< P (t)
No
Yes
To
a) In without DSM strategy (Fig. 8), hourly electrical load has been
computed based on load modeling equation as given in Eq. (26).
b) In with DSM strategy (Fig. 9), hourly electrical load have been
estimated based on modeling equation as given in Eq. (27).
1174
Table 9
Zone wise optimization results considering investment scenarios without DSM strategy.
Zones
Zone 1
Zone 2
Zone 3
Zone 4
Scenarios
HILR
MIMR
LIHR
HILR
MIMR
LIHR
HILR
MIMR
LIHR
HILR
MIMR
LIHR
MHG
BGG
BMG
SPV
WTG
BATT
6
6
6
0
0
0
15
15
15
16
16
16
11
15
27
27
38
46
26
44
60
34
45
62
18
28
30
36
50
58
45
63
75
41
68
77
43
64
75
32
49
40
65
53
71
61
76
73
20
24
31
13
6
20
15
24
33
25
28
46
109
146
190
190
238
283
244
308
380
273
340
397
TNPC ($)
COE ($/kWh)
265355
332633
408620
370176
414777
537318
531762
649497
758672
587370
691095
790509
0.105
0.108
0.108
0.124
0.125
0.130
0.118
0.114
0.114
0.129
0.113
0.119
TNPC ($)
COE ($/kWh)
238624
298243
368974
343145
388347
528308
504730
622316
725634
565550
632678
768584
0.103
0.104
0.105
0.122
0.123
0.127
0.116
0.112
0.112
0.124
0.107
0.105
TNPC ($)
COE ($/kWh)
265355
370176
649497
790509
0.105
0.124
0.114
0.119
Table 10
Zone wise optimization results considering investment scenarios with DSM strategy.
Zones
Zone 1
Zone 2
Zone 3
Zone 4
Scenarios
HILR
MIMR
LIHR
HILR
MIMR
LIHR
HILR
MIMR
LIHR
HILR
MIMR
LIHR
MHG
BGG
BMG
SPV
WTG
BATT
6
6
6
0
0
0
15
15
15
16
16
16
11
15
27
27
38
46
26
44
60
34
45
62
18
28
30
36
50
58
45
63
75
41
68
77
27
46
69
30
40
36
42
42
68
47
52
45
18
22
27
9
4
18
12
21
25
19
27
31
101
139
183
185
215
256
228
306
378
243
315
348
Table 11
Possible conguration of IRE model without DSM strategy.
Zones
Zone
Zone
Zone
Zone
1
2
3
4
Scenarios
HILR
HILR
MIMR
LIHR
MHG
BGG
BMG
SPV
WTG
BATT
6
0
15
16
11
27
44
62
18
36
63
77
43
32
53
73
20
13
24
46
109
190
308
397
zone 1 without and with DSM are shown in Fig. 6(aec). The hourly
load demand is further synthesized to hourly annual data (8760 h)
using HOMER software. Total annual load demand for all the three
scenarios of all the zones are shown in Fig. 1. Cost parameter
comprises capital cost, maintenance cost, replacement cost, fuel
cost of the system components and appliances cost of various investment scenarios used in GA based approach. Also, project lifetime and annual real interest rate are 20 years and 0.11 respectively
considered in the present study.
Optimization results in terms of net present cost and cost of
energy under with and without DSM strategies and different values
of parameters have been considered under the present study. Three
scenarios based results of IRE system with and without DSM
strategy has been obtained for all the four zones considering three
investment scenarios of low, medium and high cost of appliances
are presented in Tables 9 and 10. Out of these three, the best
possible congurations of IRE system with and without DSM
strategy are presented in Tables 11 and 12. Further, results of
possible conguration without DSM were compared with results
with DSM strategy. It has been observed that results with DSM
Table 12
Possible conguration of IRE model with DSM strategy.
Zones
Zone
Zone
Zone
Zone
1
2
3
4
Scenarios
HILR
HILR
MIMR
LIHR
MHG
BGG
BMG
SPV
WTG
BATT
6
0
15
16
11
27
44
62
18
36
63
77
27
30
42
45
18
9
21
31
101
185
306
348
TNPC ($)
COE ($/kWh)
238624
343145
622316
768584
0.103
0.122
0.112
Parameters
Zone 1
ScenarioeI: HILR
SPV
Nos.
WTG
Nos.
BATT
Nos.
TNPC
$
COE
$/kWh
ScenarioeI: HILR
SPV
Nos.
WTG
Nos.
BATT
Nos.
TNPC
$
COE
$/kWh
ScenarioeII: MIMR
SPV
Nos.
WTG
Nos.
BATT
Nos.
TNPC
$
COE
$/kWh
ScenarioeIII: LIHR
SPV
Nos.
WTG
Nos.
BATT
Nos.
TNPC
$
COE
$/kWh
Zone 2
Zone 3
Zone 4
Units
1175
achieved. Similar line, saving for zone 2, zone 3 and zone 4 are given
in Table 13.
Without DSM
With DSM
Saving
43
20
109
265355
0.105
27
18
101
238624
0.103
16
2
8
26731
0.002
32
13
190
370176
0.124
30
9
185
343145
0.122
2
4
5
27031
0.002
53
24
308
649497
0.114
42
21
306
622316
0.112
11
3
2
27181
0.002
73
46
397
790509
0.119
45
31
348
768584
0.105
28
15
49
21925
0.014
Based on the comparison of the results of possible congurations with and without DSM strategy, results of possible conguration with peak load shifting based DSM strategy are found to be
better results with ensuring EIR of 1. It involves optimal size of 6 kW
of MHG, 18 kW of BMG, 11 kW of BGG, 27 solar modules with a area
of 1.63m2, 18 wind turbine with a swept area of 4.26 m2 and 101
batteries with a rating of 3.6 kVA. The contribution of electricity
generation from AC and DC renewable generators are found to be as
171550 kWh/yr and 111735 kWh/yr respectively. Battery contributes total discharge energy of 25205 kWh/yr during a decit in the
energy generation from AC and DC generators. Total surplus power
generates from the generators was found to be 10694 kWh/yr.
Further, the minimum NPC and COE for HILR with DSM strategy are
found to be as 238624 $ and $ 0.103 per kWh respectively. Similarly,
HILR based scenario is found to be cost effective for zone 2, while
MIMR based scenario and LIHR based scenario and LIHR based
scenario are found cost effective for zone 3 and zone 4 respectively.
Summary of results for zone 2, zone 3 and zone 4 obtained and are
presented in Tables 12e14.
6.3. Hourly simulated results with and without DSM strategy
Table 14
Possible conguration of IRE model with and without DSM strategy.
Zones
Scenario
Parameters
Units
Without DSM
With DSM
Zone 1
ScenarioeI: HILR
Zone 2
ScenarioeI: HILR
Zone 3
ScenarioeII: MIMR
Zone 4
ScenarioeIII: LIHR
AC power
INV-DC power
Discharge
Excess power
AC power
INV-DC power
Discharge
Excess power
AC power
INV-DC power
Discharge
Excess power
AC power
INV-DC power
Discharge
Excess power
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
kWh/yr
171550
130153
28805
32529
256230
52244
69341
32989.04
567940
95334
104563
27829.71
703720
155920
135147
37295
171550
111735
25205
10694
256230
41431
67394
15304
567940
81368
93278
13107
703720
124785
111502
26505
1176
Power (kW)
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour (h)
Load
AC-Power
INV-DC Power
Discharge
Excess power
Power (kW)
Discharge
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour(h)
SOC
SOC-MIN
Excess power
SOC-MAX
Fig. 10. (a) Hourly electricity production versus load and (b) SOC of a battery system for a day of 13 of Feb for a zone 1 (without DSM).
Power (kW)
Load
AC-Power
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
INV-DC
Hour (h)
Discharge
Excess power
SOC
120
100
80
60
40
20
0
1
Discharge
9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Hour(h)
Excess power
SOC
SOC-MIN
SOC-MAX
Fig. 11. (a) Hourly electricity production versus load and (b) SOC of a battery system for a day of 13 Feb (with DSM) of zone 1.
1177
Table 15
Hourly simulation results of HILR scenario with and without DSM strategy for a day on 13 Feb of zone 1.
Hours
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Load
AC power
INV-DC power
Discharge
SOC
Excess power
Load
AC power
INV-DC power
Discharge
SOC
Excess power
14.43
10.18
14.1
13.75
18.61
16.35
18.23
28.28
9.14
9.66
11.66
24.09
38.16
33.79
6.33
9.44
12
13.64
56.82
47.88
36.66
19.82
10.29
11.27
6
6
6
6
6
6
6
24
24
24
35
35
35
17
17
17
24
24
35
35
35
35
6
6
4.88
0.23
0.63
0.98
0.24
0.03
0
0.4
1.77
2.56
5.82
8.74
21.19
17.48
15.97
25.39
20.83
10.92
5.94
1.87
0.12
0
0
0
3.81
3.96
7.5
6.82
12.38
10.32
12.23
3.9
0
0
0
0
0
0.18
0
0
0
0
16.19
11.1
1.54
0
4.29
5.27
95
95
90
80
70
50
40
25
20
25
30
35
65
85
80
95
99
99
99
95
95
99
99
95
0
0
0
0
0
0
0
0
0
0
0
0
0
0
17.46
29.31
29
18.7
0
0
0
0
0
0
14.43
10.17
14.08
13.73
18.58
16.33
18.2
28.25
22.26
23.54
27.12
14.42
22.85
20.23
12.02
17.93
22.79
13.62
42.95
36.2
27.71
19.8
10.28
11.26
6
6
6
6
6
6
6
24
24
24
35
35
35
17
17
17
24
24
35
35
35
35
6
6
4.39
0.21
0.56
0.88
0.22
0.03
0
0.36
1.11
1.73
3.82
6.12
16.88
13.35
12.02
20.82
17.05
9.12
5.09
1.69
0.11
0
0
0
4.26
3.97
7.55
6.89
12.38
10.3
12.2
3.9
0
0
0
0
0
0
0
0
0
0
3.11
0
0
0
4.28
5.26
95
95
90
85
65
55
40
25
20
35
55
85
95
99
99
99
99
99
99
90
85
85
95
90
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10.38
16.38
17.08
0
0
2.25
13
0
0
Table 16
Hourly simulation results of HILR scenario with and without DSM strategy for a day on 16 Jan of zone 2.
Without DSM strategy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Load
AC-power
INV-DC
Discharge
SOC
Excess power
Load
AC-power
INV-DC
Discharge
SOC
Excess power
16.67
17.45
11.62
9.84
24.86
17.57
8.88
22.27
7.51
10.15
10.48
25.37
36.52
12.08
9.92
6.86
6.26
20.03
40.75
46.24
40.23
20.5
13.79
11.53
0
0
0
0
0
0
0
36
36
36
63
63
63
27
27
27
36
36
63
63
63
63
0
0
16.51
17.02
11.37
9.84
24.86
17.55
8.44
1.2
3.31
7.44
9.94
10.31
16.99
15.14
5.73
2.79
3.22
3.24
1.13
1.51
1.15
0.19
13.75
11.53
16.51
17.02
11.37
9.84
24.86
17.55
8.44
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
13.75
11.53
65
60
55
50
45
35
25
20
25
35
45
55
65
70
75
85
99
99
99
99
99
99
99
98
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0.4
27.38
15.53
19.69
15.69
20.5
36.51
0
0
16.67
17.45
11.63
9.84
24.86
17.57
8.88
22.28
18.41
24.89
24.77
15.13
21.78
7.21
16.93
11.72
10.69
20.03
32.7
37.1
32.28
20.5
13.79
11.53
0
0
0
0
0
0
0
36
36
36
63
63
63
27
27
27
36
36
63
63
63
63
0
0
0.11
0.32
0.18
0
0
0.01
0.32
0.83
1.9
4.28
6.09
6.42
10.34
8.96
3.52
1.69
1.89
1.87
0.66
1.05
0.79
0.13
0.03
0
16.56
17.15
11.45
9.84
24.86
17.56
8.58
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
13.76
11.53
60
55
50
45
40
30
25
20
25
30
40
50
60
65
70
80
90
99
99
99
99
99
99
95
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2.67
15.34
26.51
23.09
26.98
36.46
0
0
1178
Table 17
Hourly simulation results of MIMR scenario with and without DSM strategy for a day on 10 March of zone 3.
Without DSM strategy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Load
AC-power
INV-DC power
Discharge
SOC
Excess power
Load
AC-power
INV-DC power
Discharge
SOC
Excess power
36.18
47.9
36.36
33.55
49.89
58.93
31.13
40.83
53.49
44.87
58.57
104.13
79.77
68.42
42.05
28.63
37.36
58.74
97.01
103.85
96.81
82.06
25.05
40.11
15
15
15
15
15
15
15
78
78
78
122
122
122
59
59
59
78
78
122
122
122
122
15
15
0
1.19
0.39
5.2
3.45
2.32
1.28
3.82
1.94
5.96
5.14
7.93
8.03
10.69
10.22
9.89
8.2
3.91
1.86
0.31
1.35
7.45
6.99
7.05
21.18
31.77
20.99
13.61
31.61
41.73
14.91
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
3.41
18.41
85
80
70
60
55
40
25
20
30
40
50
70
75
80
85
90
95
99
99
99
99
99
99
99
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0.97
19.99
23.03
15.8
22.75
40.85
0
0
36.18
47.9
36.36
33.56
49.89
58.94
31.13
40.83
28.75
24.11
32.03
130.47
99.95
85.73
28.27
19.25
25.12
58.74
107.1
114.65
106.87
82.06
25.05
40.11
15
15
15
15
15
15
15
78
78
78
122
122
122
59
59
59
78
78
122
122
122
122
15
15
0
1.36
0.45
5.95
3.94
2.65
1.47
4.36
2.41
7.31
6.45
9.87
9.87
13.38
12.81
12.46
10.35
4.94
2.31
0.36
1.54
8.51
7.99
7.05
21.18
31.61
20.94
12.91
31.15
41.42
14.73
0
0
0
0
0
0
14.02
0
0
0
0
0
0
0
0
2.46
18.41
70
65
60
55
50
40
25
20
30
45
60
65
65
75
70
95
99
99
99
99
99
99
99
99
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
41.38
20.91
14.81
6.6
14.32
41.81
0
0
Table 18
Hourly simulation results of LIHR scenario with and without DSM strategy for a day on 28 Jan of zone 4.
Without DSM strategy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Load
AC-power
INV-DC
Discharge
SOC
Excess energy
Load
AC-power
INV-DC
Discharge
SOC
Excess power
39.73
43.78
44.64
38.09
56.58
54.68
52.49
54.42
24.48
28.55
17.92
152.56
134.6
103.85
21.35
22.22
26.22
116.34
142.6
144.94
242.51
130.98
41.12
33.81
16
16
16
16
16
16
16
93
93
93
155
155
155
78
78
78
93
93
155
155
155
155
16
16
0
0.22
0
0
0.13
0
0.26
0.71
3.96
5.78
13.1
15.54
10.72
18.89
18.96
15.97
9.52
6.27
11.58
15.26
8.94
0.88
0.48
0.59
23.73
27.58
28.64
22.09
40.46
38.68
36.24
0
0
0
0
0
0
7.9
0
0
0
17.39
0
0
79.01
0
24.67
17.25
80
75
65
55
50
40
30
20
30
40
50
55
60
85
80
90
95
95
95
100
100
80
90
85
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
17.31
65.67
0
0.04
22.33
0
0
0
0
39.71
43.77
44.62
38.07
56.56
54.66
52.46
54.4
41.73
48.67
30.04
132.48
116.88
90.18
31.33
32.61
38.48
116.29
133.01
135.19
226.2
130.93
41.11
33.8
16
16
16
16
16
16
16
93
93
93
155
155
155
78
78
78
93
93
155
155
155
155
16
16
0
0.32
0
0
0.2
0
0.39
1.05
4.18
4.15
10.69
10.76
7.88
15.92
16.69
12.43
7.81
6.25
14.97
22.64
13.27
1.31
0.71
0.88
23.71
27.46
28.62
22.07
40.37
38.66
36.09
0
0
0
0
0
0
0
0
0
0
17.35
0
0
58.6
0
24.44
16.96
85
80
70
60
50
40
30
20
30
40
50
70
85
90
90
95
100
100
95
100
100
90
95
90
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
20.93
53.65
0
9.26
37.31
0
0
0
0
Table 19
Comparative results of GA & PSO considering three investment scenarios with DSM
strategy for zone 1.
Scenarios
Parameters
GA
PSO
HILR
TNPC ($)
Elapsed time (seconds)
No. of Iteration
TNPC($)
Elapsed time (seconds)
No. of Iteration
TNPC($)
Elapsed time (seconds)
No. of Iteration
238624
27.38
96
298243
49.94
176
368974
123
35.874
238624
28.258
114
298093
28.258
113
368974
119
29.796
MIMR
LIHR
supplying power to four zones in Chamarajanagar district of Karnataka state of India. TNPC and COE of the IRE system were determined by considering peak load shifting based DSM strategy.
It is observed that system cost, system size and load density are
different for all considered scenarios (HILR, MIMR and LIHR). Load
density has been estimated by considering different appliances and
is found to be increased from HILR to LIHR scenarios which results
in the increase of system size and system cost.
Cost of energy has been computed as per total net present cost
of the system (including system cost, appliances cost and saving
cost) divided by total annual energy generation.
Fig.12. Convergence curve of GA and PSO for HILR scenarios for zone 1.
Fig.13. Convergence curve of GA & PSO (combined) for HILR, MIMR and LIHR scenarios for zone 1.
1179
1180