Simultaneous Allocation of Multi-DG With Capacitor Banks in Distribution System Using IGWO
Simultaneous Allocation of Multi-DG With Capacitor Banks in Distribution System Using IGWO
net/publication/362808454
CITATION READS
1 13
4 authors, including:
All content following this page was uploaded by Sarib Malik on 15 July 2024.
Luigi Martirano
Department of Astronautical, Electrical and Energy Engineering (DIAEE)
Sapienza University of Roma
Roma, Italy
[email protected]
Abstract— Integration of Renewable based distributed inexhaustible nature of RERs. According to the U.S
generations (DGs) along with capacitor banks (CBs) aim to Department of energy, Solar PV is the fastest-growing power-
improve the performance of the distribution system (DS). This generating RER around the world [3]. 1n 2020 almost 139
paper presents an Improved grey wolf optimizer (IGWO) for GW of overall global capacity was integrated from PV into
simultaneous allocation of multi-DG with CBs in DS. The
the power system bringing it to almost 760 GW which is about
proposed technique aims to enhance technical, environmental,
and economical benefits. The multi-objective examined in this 3% of the total world’s electricity [3]. While wind power
paper are reduction of power losses, voltage profile achieved a record-breaking 93 GW new installation, to sum
improvement, cost minimization, and control of harmful up to a total of 743 GW [3]. A global total of 2799 GW of
environmental gases. Widely used and most common renewable energy has been generated from RERs in 2020 and it
energy resources (RERs) are considered for integration in DS. integrated more than 290 GW in 2021 [3].
Solar Photovoltaic (PV) and wind-based DGs are integrated into
DS with probabilistic generation to fulfill load demand in order The integration of renewable-based DG sources has the
to minimize multi-objective function (MOF). The IGWO main role of power injection along with a few different
algorithm is implemented on radial DS namely IEEE 33-bus and benefits in DS [4,5]. Mostly DG units are injected into DS in
69-bus system. The proposed IGWO provides better results as order to reduce power losses (active (P) and reactive (Q)),
compared to optimization techniques available in the literature. improve system voltage profile, enhance the power quality,
and reliability of DS. [6,7]. DG units are used to enhance Q
Keywords—Distributed generations, capacitor banks
power control, minimization of losses, and frequency control
renewable-based DG, distribution system, improved grey wolf
optimizer, photovoltaic, and wind.
[8,9]. In [10-12] DG sources are used to reduce losses,
improve load factor, system reinforcement, and enhancement
INTRODUCTION of frequency. Evolutionary programming is considered for the
Most of the power supply all around the world is generated integration of DG to reduce losses and minimize cost [13].
from conventional power sources. These conventional sources CBs can generate reactive power that can lead to
are mostly based on fossil fuels such as coal, oil, gas, thermal, improvement of voltages and minimize power losses. Several
and nuclear power plants. There are a lot of hazards from these studies discussed the integration of DGs in DS and several
sources such as high cost and emission of non-environmental deals integration of CBs in DS but there are few which
friendly harmful greenhouse gases. Moreover, these fossil considered them both. Saonerkar and Bagde presented the
fuels are depleting with the passage of time, and the demand optimal location and sizing of DGs and CBs in DS using a
for power supply is increasing every day [1]. To overcome genetic algorithm [14]. In [15] optimal siting and sizing of
these problems and meet the power demand, renewable-based DGs along with CBs while considering the uncertainty of
DG sources are considered the best alternative to traditional RERs in order to overcome the problems. In [16] optimal
energy sources for power generation [2]. These sources sizing and location of DGs along with CBs are considered in
consist of PV, wind, geothermal, biomass, etc. are inevitable DS by applying the water cycle algorithm (WCA). Optimal
and help to mitigate environmental hazards and provide siting and sizing of DGs along with CBs are examined in [17]
economic power supply. by using constriction factor particle swarm optimization
Electricity generation from renewable energy resources (PSO).
(RERs) is rapidly growing nowadays due to the depletion of
fossil fuels, increase in power demand, and everlasting or This paper presents an IGWO to investigate the
simultaneous optimal allocation of multi-DG units along with
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
CBs in DS for the sake of minimization of MOF. The MOF PROBLEM FORMULATION
includes reduction of losses, voltage profile improvement,
Objective Functions
cost minimization and environmental benefits. The IGWO is
an improved version of GWO published by Syed Ali Mir Jalili An IGWO algorithm is used for the simultaneous
in 2021 [18]. The results are compared with the other meta- allocation of multi-DG units in DS to enhance technical (P
heuristic techniques in the literature and found better. and Q power losses, V profile improvement and V stability),
economical, and environmental objective OFs. These OFs are
The section organized in this paper are as follow: section discussed in detail in this section.
1 portrayed introduction, section 2 presents DG modelling,
section 3 depicts problem formulation, section 4 deals with 1. Technical OF:
proposed technique. Section 5 determines case study and
The technical OF comprises of power losses reduction,
section 6 discusses conclusions.
voltage profile improvement and voltage stability index. The
I. DG MODELLING first OF aims to reduce power losses f1 of DS and can be
presented as[16]:
Solar PV DG Modelling
|Ii |2 * Ri
nL
A Beta probability distribution function (PDF) is
f1 x =min (5)
considered for PV power generation [19]. A beta PDF is
i=1
utilized in order to define the variability of solar irradiance (s)
The voltage profile improvement is considered as second
of each state for PV panel and is presented as.
technical OF that aims to enhance the voltage profile (f2 ) of
Γ αt + βt αt -1 βt -1 DS and this OF is depicted as[16]:
ft s = Γ(αt )* Γ(β ) *s * 1-s 0≤ s≥1 (1)
N 2
t vi -vsi
0 otherwise f2 x =min (6)
Where α and β are solar pattern parameter and are vmax
i - vi
min
i=0
represented in [19]. The last but not the least voltage stability index (VSI) is
From beta PDF the power generated by PV with installed which is one of the most significant technical OF . This third
r
rated capacity Psolar of PV at each state t with solar irradiance technical OF (f3 ) is represented as given in equation below
(s) is given by [20]. [16]:
Psolar (s)= r
s r
*Psolar 0 ≤ s ≥ sr 1
t s (2) f3 x =min (7)
r VSI b2
Psolar otherwise
ABS(|V n1 |4 -4*[P n2 *Pij -Q f1 *Rij ]2
Where
8
t
Where Psolar is the power of PV at each state t, sr is rated
r = 2
solar irradiance and Psolar is rated solar capacity. The value
of rated capacity and solar irradiance is given in [20]. -4* P n2 *P11+Q n2 *Pij *|V b1 |)
Where n1 is slack bus and n is the selected bus number.
Wind DG Modeling 2. Economic OF:
Wind speed (v) act as the deciding factor for the power
generation by wind DG. To find the wind speed generally One of the most significant is economical OF (f4 ) that
Weibull distribution is considered practically suitable [21]. A aims to reduce the generation cost and is presented as [16]
NDG
Weibull distribution function (WDF) for generation of v at
any point in time interval (t) can be calculated as [22]. f4 x =min CDGi +Csub (9)
Kt v Kt-1 -v Kt i=1
fw v = * *exp ;for Ct ≥1& Kt ≥0 (3) Where
Ct Ct Ct
Where fw v is the WDF of wind speed, Kt is the shape CDGi =a+b*PG i
while Ct is the scale factor and values can be obtained using $
capital cost * capacity kW * Gr
wind data available in article [22]. a= kW
life years * LF* 8760
The power generated from wind turbine (WT) with
probabilities of wind v, rated capacity, at time interval t is b= O&Mcost ($/kW·h) + fuel cost ($/kW·h)
calculated by using equation represented in [22].
⎧ ∑Ni=1
0 if vci ≥v≥vco Csub =PgGrid *PrGrid
c c
v-vi ei +Cci |Qci |
⎨
Pstw v = Prw * c c if vci ≤v≤vco (4) CCB =
⎩ Pw
v r -vo life years * 8760
r
otherwise Where eI is CBs installation cost, Cci is CBs purchase cost
Where Pstw is the power generated by WT at speed v, Prw is and Qci is CBs initial cost. The values of these parameters are
the rated power of WT, vci is the cut-in speed, vco is the cut-out given in Table 1.
speed and vcr is the rated speed. The cut-in wind v is 4 m/s, 3. Environmental OF:
while rated wind v is 14 m/s and cut-off wind v is 24 m/s [23].
The most hazardous pollutants produced as a result of
power generation through fossil fuels are CO2, SO2 and NOx.
The reduction of these emission during power generation is a
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
challenging task and mathematical model to achieve this OF P t+1 =Pp t -A*D (19)
is represented by (f5 ) and is as follows [16]: Where Pp indicates position vector, P indicates prey’s
NDG
position, t is current iteration, and coefficient vectors A and
b = 1, 2, . . . , n-bus.
using the equation (23).
0.95 ≤ V ≤ 1.05, (15) P1 t +P2 t +P3 t
P(t+1)= (25)
Installed CBs limit 3
a) Attacking
QCB
T
<QD (16) The attacking process starts when the hunting process
DGs Power Factor (PF) Limits terminates. According to [18] half iterations are mentioned to
the exploration stage while half iterations are assigned to the
0.8≤ PF≤1 (17) exploitation. Exploration deals with smooth adjustment while
exploitation refers to the random change of position and
PROPOSED TECHNIQUE current position of prey.
The IGWO algorithm is used for the simultaneous 2. Improved Grey Wolf Optimizer (IGWO)
allocation of multi-DG units with CBs. GWO algorithm is The best solution α, β and δ may entrap in local best
used to find the optimal solution first, then the solution obtain optima instead of finding optimal solution. To overcome this
by this algorithm is updated by IGWO in order to reduce problem an IGWO [18] is used instead of GWO. The
stuck in local optima. The detailed working philosophy of the improvement made in GWO are considered by selecting and
IGWO is presented in this section. updating the following steps [18].
a) Initializing
Working of Improved Grey Wolf Optimizer
In initializing phase number of wolves (N) are randomly
1. Grey Wolf Optimizer (GWO) allocated all over the search space with given range in
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
Then after this neighbor N of P respected to R is CASE STUDY
calculated by using the equation (27). Where D is the
Test System
Euclidean radius.
The proposed IGWO is applied in two IEEE distribution
system. These DS are IEEE 33-bus and IEEE 69-bus system.
The total real and reactive power for IEEE 33-bus system is
3715 KW and 2300 KVar. While the real and reactive power
for IEEE 69-bus system is 3802 KW and 2694 KVar [16]. The
power flow for the base case is carried out by using forward-
backward sweeping method and is represented in paper [26].
The base values for active power loss for IEEE 33-bus system
is 201.97 KW while active power loss for IEEE 69 bus system
is 225 KW.
Multi-Objective Function
The MOI is considered for the optimal allocation of CBs
and DGs in DS. DGs are used with controllable power factor
which supply active and reactive power. Two cases for OF
with weighing factors are considered.
Table. 1. Characteristics of Different types of DGs
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
proposed technique is more reliable in determining
simultaneous allocation of multi-DG with CBs in DS. The 1.04 Base Case
Table. 2. Comparison of Results for Simultaneous Allocation of DGs along With CBs
with CBs in IEEE 33-Bus System 1.02 With DGs
With DGs and CBs
PL DGs CBs Vmin 1.00
(KW (p.u)
Voltage (p.u)
Algo.
) Site Size Site Size 0.98
(MW) (MVar)
GA 71.25 16,22 0.25, 15,18, 0.3, 0.3, 0.971 0.96
[16] ,30 0.25, 0.5 29,30, 0.3, 0.6,0.3
31 0.94
BFO 41.41 17,18 0.54, 18,30,
0.16, 0.33, 0.978
A ,33 0.16, 0.89 33 0.54 0.92
[16]
WCA 24.68 11,25 0.97, 14,23, 0.46, 0.56, 0.98 0.90
[16] ,29 1.04, 0.56 30 0.53 (33) 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32
IGW 11.13 14,24 0.36, 13,25, 0.33, 0.24, 0.99 Bus No. (#)
O ,30 0.42, 0.53 31 0.41 (33)
Table. 3. Results Comparison for Multi-objective Simultaneous Allocation Fig. 2. Voltage Profiles for IEEE 33-Bus System
of DGs along with CBs in 33-Bus System
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
DGs and CBs minimizes the losses and enhance the [10] V. Pandi, H. Zeineldin, W. Xiao, ''Determining optimal location and
size of distributed generation resources considering harmonic and
utilization of energy in DS. Simultaneous allocation of multi- protection coordination limits,'' IEEE Trans. Power Syst. 28 (2) (2013)
DG units along with CBs provides better performance as 1245–1254.
compared to DGs and CBs integrated separately. The [11] S. Ghosh, S. Ghoshal, S. Ghosh, ''Optimal sizing and placement of
controllable power factor DGs provide better results as distributed generation in a network system,'' Int. J. Electr. Power
compared to unity power factor DGs. Energy Syst. 32 (8) (2010) 849–856.
The results of IEEE 69 bus system for Case 2 are presented [12] R. Al Abri, E. El-Saadany, Y. Atwa, ''Optimal placement and sizing
method to improve the voltage stability margin in a distribution system
in Table. 6. A significant amount of power losses reduction using distributed generation,'' IEEE Trans. Power Syst. 28 (1) (2013)
is obtained which is almost 214.88 KW. Along with that an 326–334.
effective minimization of generation cost and emission are [13] D. Khatod, V. Pant, J. Sharma, ''Evolutionary programming based
achieved. The power losses are minimized from 225 KW to optimal placement of renewable distributed generators,'' IEEE Trans.
10.22 KW and the generation cost for IEEE 69 bus system is Power Syst. 28 (2) (2013) 683–695.
minimized from 309.71 $/h to 272.24 $/h. While the [14] Saonerkar, A. K., and B. Y. Bagde. "Optimized DG placement in radial
distribution system with reconfiguration and capacitor placement using
environmental emission is controlled and minimized from genetic algorithm." 2014 IEEE international conference on advanced
82.50e5 to 4.13e6 lb/h. The voltage profiles for simultaneous communications, control and computing technologies. IEEE, 2014.
allocation of DGs along with CBs in IEEE 69-bus system are [15] Jannat, Mohamed B., and Aleksandar S. Savić. "Optimal capacitor
depicted in Fig. 3. placement in distribution networks regarding uncertainty in active
power load and distributed generation units production." IET
CONCLUSION Generation, Transmission & Distribution 10.12 (2016): 3060-3067.
[16] Abou El-Ela, Adel A., Ragab A. El-Sehiemy, and Ahmed Samir
This paper presents a simultaneous allocation of multi-DG Abbas. "Optimal placement and sizing of distributed generation and
with CBs in DS using IGWO. In addition, this study focusses capacitor banks in distribution systems using water cycle
on the MOF in order to enhance performance of DS. The algorithm." IEEE Systems Journal 12.4 (2018): 3629-3636.
MOF comprises of reduction of power losses, enhance [17] Balu, Korra, and Vivekananda Mukherjee. "Siting and sizing of
distributed generation and shunt capacitor banks in radial distribution
voltage profile, maintain voltage stability, minimization of system using constriction factor particle swarm optimization." Electric
generation cost and control of harmful gases. Renewable Power Components and Systems 48.6-7 (2020): 697-710.
based PV and wind along with diesel generator are [18] Nadimi-Shahraki, Mohammad H., Shokooh Taghian, and Seyedali
considered for integration in DS. Two different cases for Mirjalili. "An improved grey wolf optimizer for solving engineering
problems." Expert Systems with Applications 166 (2021): 113917.
IEEE 33-bus IEEE-69 bus DS are examined in order to find
[19] Hung, Duong Quoc, Nadarajah Mithulananthan, and Kwang Y. Lee.
MOF. The results obtained from IGWO are compared with "Determining PV penetration for distribution systems with time-
other optimization techniques. It depicts that the proposed varying load models." IEEE Transactions on Power Systems 29, no. 6
technique significantly provides better results than other (2014): 3048-3057.
optimization algorithms available in literature and also [20] Zakaria, A., Firas B. Ismail, MS Hossain Lipu, and Mahammad Abdul
improve system performance. Hannan. "Uncertainty models for stochastic optimization in renewable
energy applications." Renewable Energy 145 (2020): 1543-1571.
REFERENCES [21] Soroudi, Alireza, Morteza Aien, and Mehdi Ehsan. "A probabilistic
modeling of photo voltaic modules and wind power generation impact
on distribution networks." IEEE Systems Journal 6, no. 2 (2011): 254-
[1] L. Chiari and A. J. E. P. Zecca, "Constraints of fossil fuels depletion on 259.
global warming projections," vol. 39, no. 9, pp. 5026-5034, 2011. [22] Jithendranath, J., and Debapriya Das. "Scenario-based multi-objective
[2] Demirtas, Ozgur. "Evaluating the best renewable energy technology for optimisation with loadability in islanded microgrids considering load
sustainable energy planning." International Journal of Energy and renewable generation uncertainties." IET Renewable Power
Economics and Policy 3 (2013): 23. Generation 13, no. 5 (2019): 785-800.
[3] Renewables 2021 Global Status Report, [Online]. Available: http: [23] Atwa, Y. M., E. F. El-Saadany, M. M. A. Salama, and R. Seethapathy.
www.ren21.net "Optimal renewable resources mix for distribution system energy loss
[4] D. Hung, N. Mithulananthan, ''Multiple distributed generator minimization." IEEE Transactions on Power Systems 25, no. 1 (2009):
placement in primary distribution networks for loss reduction,'' IEEE 360-370.
Trans. Ind. Electron. 60 (4) (2013) 1700–1708. [24] Emary, Eid, Hossam M. Zawbaa, and Crina Grosan. "Experienced gray
[5] R. Rao, K. Ravindra, K. Satish, S. Narasimham, ''Power loss wolf optimization through reinforcement learning and neural
minimization in distribution system using network reconfiguration in networks." IEEE transactions on neural networks and learning
the presence of distributed generation,'' IEEE Trans. Power Syst. 28 (1) systems 29, no. 3 (2017): 681-694.
(2013) 317–325. [25] Mirjalili, Seyedali, Seyed Mohammad Mirjalili, and Andrew Lewis.
[6] D. Hung, N. Mithulananthan, R. Bansal, ''Analytical expressions for "Grey wolf optimizer." Advances in engineering software 69 (2014):
DG allocation in primary distribution networks,'' IEEE Trans. Energy 46-61.
Convers. 25 (3) (2010) 814–820. [26] Abou El-Ela, Adel Ali, Ragab A. El-Sehiemy, Abdel-Mohsen Kinawy,
[7] D. Popović, J. Greatbanks, M. Begović, A. Pregelj, ''Placement of and Mohamed Taha Mouwafi. "Optimal capacitor placement in
distributed generators and reclosers for distribution network security distribution systems for power loss reduction and voltage profile
and reliability,'' Int. J. Electr. Power Energy Syst. 27 (5–6) (2005) 398– improvement." IET Generation, Transmission & Distribution 10, no. 5
408. (2016): 1209-1221.
[8] A. Mohamed Imran, M. Kowsalya, ''Optimal size and sitting of [27] Shaheen, Abdullah Mohammed, Ragab A. El-Sehiemy, and Sobhy M.
multiple distributed generators in distribution system using bacterial Farrag. "Adequate planning of shunt power capacitors involving
foraging optimization,'' Swarm Evol. Comput. 15 (2014) 58–65. transformer capacity release benefit." IEEE Systems Journal 12, no. 1
(2015): 373-382.
[9] S. H, Lee, J.W. Park, ''Selection of optimal location and size of multiple
distributed generations by using kalman filters,'' IEEE Trans. Power
Syst. 24 (3) (2009) 1393–1400.
Authorized licensed use limited to: UNIV OF ENGINEERING AND TECHNOLOGY TAXILA. Downloaded on August 22,2022 at 18:26:12 UTC from IEEE Xplore. Restrictions apply.
View publication stats