Optimal Network Reconfiguration With Distributed G
Optimal Network Reconfiguration With Distributed G
Abstract—This paper presents a method to solve electrical Network reconfiguration is a very effective and efficient
network reconfiguration problem in the presence of distributed way to ensure more even load distribution of network’s
generation (DG) with an objective of minimizing real power loss elements, improve system reliability and voltage profile, and
and energy not supplied function in distribution system. A to reduce power losses. All modes are subject to
method based on NSGA II multi-objective algorithm is used to reconfiguration: normal, critical and failure. Provided that all
simultaneously minimize two objective functions and to identify variables are within acceptable limits, network reconfiguration
the optimal distribution network topology. The constraints of will achieve optimal working conditions in normal mode.
voltage and branch current carrying capacity are included in the
evaluation of the objective function. The method has been tested Taking into consideration a large number of switches in
on radial electrical distribution network with 213 nodes, 248 lines distribution network, whose on/off switching affects the
and 72 switches. Numerical results are presented to demonstrate network topology, reconfiguration problem can be defined as
the performance and effectiveness of the proposed methodology. a complex combinatorial, non-differentiable, and constrained
multi-objective optimization problem. Radial network
Keywords—radial distribution network; distributed generation; conditions, explicit voltage constraints in all node, line
genetic algorithms; NSGA II; loss reduction capacities, etc. are viewed as some of the constraints that have
I. INTRODUCTION to be taken into consideration.
Newly formed market conditions articulate the need for In recent years, different methods and approaches have
adjusted approach in managing distribution network in order been used to solve the problem of distribution system
to meet not only the requirements imposed by technical reconfiguration with distribution generators installed. The
conditions of the system but also the requirements imposed by literature related to this problem mainly refers to application
consumers and network regulators. Significant changes in of heuristic algorithms and artificial intelligence-based
distribution system have been caused by installing distribution algorithms such as Genetic algorithms, Fuzzy logic, Particle
generation units which have considerable impact on system swarm optimization, Tabu search, etc. [2-7]. Most cases
voltage profile and power losses, both being important address reducing power losses and load balancing, taking into
quantities in the process of planning and reconfiguration of account the effect of generators distributed in the network,
electrical network. while very little attention is paid to system reliability.
However, special attention should been paid to the issue of
DGs are grid-connected or stand-alone electric generation reliability of power supply, in order to increase economic
units located within the electric distribution system at or near efficiency of distribution companies [14]. Network
the end user. The integration of DGs in distribution system reconfiguration process can be used as one possibility to
would lead to improving the voltage profile, reliability improve network reliability indicators. Furthermore, reliability
improvement such as service restoration and uninterruptible improvement by DGs is possible when intended islanding
power supply and increase in energy efficiency. The operation is allowed [8].
distribution feeder reconfiguration (DFR) is one of the most
significant control schemes in the distribution networks which In paper [15] NSGA II (Non-Dominated Genetic
can be affected by the interconnection of DGs [1]. Algorithm II) is applied to the planning of distribution
electrical network problem. This paper is focused on the
Generally, the DFR is defined as altering the topological application of NSGA II on resolving the problem of
structure of distribution feeders by changing the open/closed reconfiguration of distribution network with distributed
status of automatic and tie switches or protective devices generation. The effect of distributed generation on voltage in
located in strategic places in distribution system. By changing the network, taking into account two objective functions:
the statuses of the sectionalizing and tie switches, the power losses and reliability function is presented. Depending
configuration of distribution system is varied, and loads are on characteristics of the power distribution networks (network
transferred among the feeders while the radial configuration parameters, characteristics of power lines, failure rates, types
format of electrical supply is still maintained. of consumers, etc) simultaneous optimization of these
6|P age
www.ijarai.thesai.org
(IJARAI) International Journal of Advanced Research in Artificial Intelligence,
Vol. 5, No.2, 2016
two objective functions simultaneously, the problem is defined where Di = 1 if line i , otherwise is equal to 0.
as multi-objective problem taking into account the defined
system constraints. The effectiveness of the methodology is NV
demonstrated on real distribution network consisting of 213-
buses showing its potential of applicability to the large
min f ens P r
i 1
i i i (2)
distribution systems. where Ploss is total real power losses function, fens is the
ENS function, Ii current in the branch i; Ri resistance in the
The problem formulation is discussed in detail in section
branch i; α is set of branches connected to the distribution
II. The network reconfiguration algorithm using a multi-
generators node m, Pi real power flow through branch i; λi
objective NSGA II is described in section III. The simulation
failure rate of branch i (number of failure per year and per
results in terms of power loss and energy not supplied are
kilometer of branch i); ri failure duration of branch i; Nv
discussed in Section IV and finally the last section presents the
number of branches. The DG produces active current I DI, and
conclusion of the study.
for a radial network it changes only the active component of
II. PROBLEM FORMULATION current of branch set α.
The objective of the proposed solution model for Subject to the system constraints:
reconfiguration of distribution network problem is to
minimize two functions as follows: power losses function and - I i I maxi - the current in each branch cannot exceed the
reliability function presented by Energy not supplied index branch capacity,
(ENS). The optimal results of the defined functions do not - Vmin V j Vmax - voltage constraint.
lead to the same optimal network topology what creates trade-
off between reliability and power losses function. The voltage in each load buses in the system has to be
within the defined limits. The minimum voltage is 0.95 and
Electrical power losses are one of the most important maximum voltage is 1.05 (±5%).
factors which point to business cost-effectiveness and quality
of distribution. Energy losses in electrical distribution network n n
in the amount of 1% cause the increase in company’s business
costs of up to 2% to supply energy to cover the losses [9].
- PS
i 1
Pg ,i P P
i 1
i loss - power balance constraint
and
The purpose of distribution network reconfiguration is to T c T
Vi Vi a , Vib , Vic i , i , i
find optimal radial operating structure that minimizes two a b
i are
functions: the system power losses and ENS function within
the operating constraints. Accorindg to the literature [2], [14] voltage magnitude and angels respectively for load buses i,
thus the problem can be formulated as follows: and the bus 1 is substation.
7|P age
www.ijarai.thesai.org
(IJARAI) International Journal of Advanced Research in Artificial Intelligence,
Vol. 5, No.2, 2016
For three-phase distribution network with n buses, bus 1 with the same solution dominance. Furthermore, NSGA II
presents substation and buses 2, 3, ...n, are load buses. uses elite strategy that significantly helps in speeding up the
Equation (3) can be solved by power flow calculation solving performance of the genetic algorithm [13].
the system of (6n-6) non-linear algebraic equations.
B. Proposed Methodology
To calculate the second objective function, ENS index due Algorithm starts with randomly selected radial functional
to interruption in supply, it is necessary to consider two solution that is typical for electrical distribution network, as a
elements: failure rate and the length of interruption in power basis for a first generation of trade-offs in the part of a genetic
supply for each load point. The latter is consisted of two algorithm code. By applying NSGA II algorithm new potential
components: time necessary to locate failure and the time solutions of the network are generated. The binary alphabet
necessary to repair it. Automatic sectionalizers and switches has been used to implement the optimization model, in which
separate the part of the network where the failure occurred, every bite of chromosome represents the status of switches
reducing the risk for other consumers in the network. The time (open/closed). Every bite can have value of 0 or 1, which
needed to repair the failure is usually the time needed to identifies the status of every electric line, 0-open, 1-closed. In
isolate the failure, to connect the affected consumers to the reconfiguration network problems, only a certain number
reserve power supply (if possible) and to repair the fault itself of lines have a changeable bit in chromosome, and therefore
[10]. In order to calculate this function, load flow studies only those lines are subject to genetic operators, crossover and
should be performed to calculate not-distributed energy in all mutation, while other lines have a fixed value in chromosome
consumer nodes without supply, which are located ―under‖ the (always have the value of 1 in operation).
fault in the network.
If newly created solutions meet topological constraints
III. OPTIMIZATION METHOD FOR MULTI-OBJECTIVE (radial conditions), evaluation of the objective functions is
RECONFIGURATION NETWORK carried out, i.e. power flow and calculation of objective
The development of heuristic algorithms and computer functions are performed. Power flow calculation is done in
performances have contributed towards solving the problem of MATLAB. For that purpose, a part of the code for power flow
multi-objective optimization. While solving multi-objective calculation based on Newton-Raphson method is modified for
optimization problems, it is necessary to pay attention to the need of objective functions evaluation, transfer of
convergence to optimal set of solutions (Pareto set) and variables, storage of diverging solutions and visualization.
maintain diversity of solutions within the set of current Based on the power flow results, convergence of specific
solutions [11]. Suggested methodology for solving the defined network configuration is verified, as well as other constraints
multi-objective optimization problem is based on multi- which refer to the capacities of lines, power stations and
objective Non-Dominated Sorting Genetic Algorithm II. distributed generations. Solutions which do not satisfy defined
(NSGA II). constraints are eliminated or penalized, depending on
convergence of power flow calculation. For other solutions,
Genetic algorithms use population of solutions in every which meet defined limits, evaluation of objective function is
optimization path within optimization process. The objective done.
is to come as close as possible to the true Pareto-front and
simultaneously gain as many solutions as possible. This The procedure is repeated until stopping criteria are met.
ensures that the decision-maker will have a wider choice of The criteria for stopping calculation can be based on a
quality solutions with a better overview of all possible optimal maximum number of generations, minimum of evaluated
topologies of a distribution network [11]. solutions, time limits to simulation, average change in solution
distribution, etc.
A. NSGA II Algorithm
Suggested model uses the concept of Pareto domination in
Multi-objective evolutionary algorithms are suitable for the evaluation of the objective functions. Input data to
multi-objective optimization due to their ability to handle describe multi-objective optimization problem are system
complex problems, involving features such as discontinuities, parameters and constraints, lines, the loads, reliability
multimodality, disjoint feasible spaces and noise functions parameters, failure rates, and the repair times.
evaluation [12]. NSGA II is a multi-objective genetic
algorithm developed by Deb, 2003 [13]. Basic advantage of MATLAB functions for genetic algorithms which are used
NSGA II over other multi-objective genetic algorithms is for calculation are modified for specific discrete function for
reflected in possibilities for diversity preservation of calculating power flow, power losses and testing system’s
population, which further enables uniform distribution of constraints for network solutions.
solutions within Pareto front. The crowding distance approach CPU time spent for calculating and identifying the set of
is introduced into NSGA II as the fitness measure to make possible solutions depends on the time necessary for the
comparison of solutions in the same Front. This approach objective functions evaluation, time for verification of defined
estimates the density of solutions surrounding a particular constraints by power flow calculation and active losses. To
solution by calculating the average distance of two points on speed up the calculation, parallel processing of genetic
either side of the observed solution for all objective functions algorithm, in the part of evaluation of the objective function
defined for particular problem. The fast non-dominated sort and constraints, is done.
strategy is used to evaluate solution dominance and classify
the solution into Pareto fronts that corresponds to the cluster
8|P age
www.ijarai.thesai.org
(IJARAI) International Journal of Advanced Research in Artificial Intelligence,
Vol. 5, No.2, 2016
9|P age
www.ijarai.thesai.org
(IJARAI) International Journal of Advanced Research in Artificial Intelligence,
Vol. 5, No.2, 2016
A. Result Analysis solutions depends on all set values, where different intensities
Application of the described methodology to identify and length of fault on lines are of special importance.
optimal network configuration based on NSGA II algorithm If intensity and fault length at all lines are equal, the
resulted in a set of possible solutions, out of which 9 Pareto variability of solution would be considerably lower, with a
optimal solutions are presented in Fig. 2. Pareto set of optimal unique optimum for both functions. It is obvious that many
solutions is achieved for 138 generations, while total number searched solutions can be simultaneously improved from the
of achieved possible solutions is 49. aspect of both functions which are optimized; subsequently,
Total algorithm execution time on 18 processor units was the considered objectives are not necessarily in conflict with
15 minutes and 36 seconds. each other. This does not provide values which are
approximately optimal for either of the objective functions.
Table 1. shows values of objective functions for solutions
from Pareto optimal set, as well as the on/off change of Objective function assessing reliability, interruptions in
switches. supply, is incidental. Therefore, when descending sort order
strategy is applied, the probability for local minimum is higher
for reliability criterion that for power losses function.
Optimal minimization of losses will be achieved when the
voltage in lines is closer to the maximum allowed value U max.
Since calculations were done with assumed constant load
values (values of peak loads), maintaining voltage at lines as
closer as possible to Umax ensures considerably less values for
losses in the network. If voltage limitation is U max = 1.20,
losses values for solution 3 would be 11.256 KW.
Power losses for line 211-212 for initial solution are 0.73
KW, and for the identified optimal solution it is 0.49 KW
(solution 3). Values for losses at the same line without
distributed generation connected is 0.58 KW with the same
Fig. 2. Pareto optimal solutions for network reconfiguration with distributed switch state.
generation
Since distributed generation is of small capacity in relation
TABLE I. VALUES OF OBJECTIVE FUNCTION FOR PARETO OPTIMAL to the strength of distributive network into which it is
SOLUTIONS FROM FIG. 2 connected, there is a reduction in the losses in the line onto
Losses ENS which it is connected.
No. Switch on Switch off
KW KWh
1. 62 3 20.6344 24.4875
Impact on voltage
2. 7 72 23.7643 24.2683 On Fig. 3. is shown voltage profile of bus system for
3. 63 6 25.8799 24.2139 solution 3 from Table 1. It is clear that the voltage is within
4. 68 67 37.1474 23.5844 the allowed limits. On Fig. 4. is shown the change in voltage
5. 60 59 38.7454 23.2314 in network nodes when distributed generation is not connected
6. 31 5 40.1339 23.1442 into the network (voltage shown in blue) and when the
7. 62 5 44.4429 21.8045 distributed generation is connected (voltage shown in red).
8. 71 61 46.8351 21.6585 The shapes of voltage profiles are almost the same in both
9. 63 8 52.2096 21.2066 cases, except for minor changes in the voltage strength at end
lines, which is a consequence of connected distributed
Considering the results shown in Table 1, the best solution
generation. Lines with distributive generation connected have
for losses function is solution 1. However, this is the worst
an increase in voltage from 0.9723 p.u. to 0.9813 p.u. after
solution for ENS function. It is evident that changes in the
installing the distributed generation into the network.
value for ENS function are smaller than for losses function.
Therefore, based on practical network topology
implementation, functional and economical benefit the best
compromise solution can be solution 3, given the evident
small changes in ENS function between solutions 3, 4 and 5,
while the difference in losses function is somewhat bigger.
Solution 3 has power losses of 25.8799 KW, while ENS is
24.2139 KWh.
The achieved near-optimal solutions show traits of each
solution from the Pareto front (the fact that not a single
individual solution from Pareto front can be improved for one
function without affecting the other in the opposite manner).
This trait does not apply for all permissible searched solutions.
The character of the achieved searched and near-optimal Fig. 3. Voltage profile for solution 3
10 | P a g e
www.ijarai.thesai.org
(IJARAI) International Journal of Advanced Research in Artificial Intelligence,
Vol. 5, No.2, 2016
REFERENCES
[1] W. M. Dahalan, H. Mokhlis, “Network Reconfiguration for Loss
Reduction with Distributed Generations Using PSO‖, IEEE International
Conference on Power and Energy (PECon), Malaysia, pp. 823-828G,
December 2012.
[2] N. Rugthaicharoencheep, S. Sirisumranukul, Optimal Feeder
Reconfiguration with Distributed Generation in Three-Phase
Distribution System by Fuzzy Multiobjective and Tabu Search, Energy
Technology and Management, 1st ed, Tauseef Aized, ISBN 978-953-
307-742-0, October , 2011 under CC BY-NC-SA 3.0 license
[3] M. N. M. Nasir, N. M. Shahrin, M. F. Sulaima, M. H. Jali, M. F.
Baharoma, ―Optimum Network Reconfiguration and DGs Sizing With
Allocation Simultaneously by Using Particle Swarm Optimization
(PSO)‖, International Journal of Engineering and Technology (IJET),
Vol 6, No 2, pp. 773-780, Apr-May 2014.
[4] T.P. Sai Sree, N.P. Chandra Rao, ―Optimal Network Reconfiguration
and Loss Minimization Using Harmony Search Algorithm In The
Presence Of Distributed Generation‖, International Journal of Electrical
Fig. 4. Voltage profile with distributed generation connected and and Electronics Research, Vol. 2, Issue 3, pp. 251-265,July - September
disconnected 2014.
[5] Y. Shu-jun, Y.Zhi, W.Yan, Y. Yu-xin, S.Xiao-yan, ―Distribution
V. CONCLUSION network reconfiguration with distributed power based on genetic
algorithm‖, 4th International Conference on Electric Utility
Reconfiguration represents one of the most important Deregulation and Restructuring and Power Technologies (DRPT), ,pp.
measures which can improve performance in the operation of 811-815, 6-9 July 2011.
a distribution system. The paper shows application of multi- [6] A.Hajizadeh, E. Hajizadeh, ―PSO-Based Planning of Distribution
objective genetic algorithm NSGA II on resolving the problem Systems with Distributed Generators‖, World Academy of Science,
of reconfiguration of distributive network with distributed Engineering and Technology 45, pp.598-603, 2008.
generation, in order to identify optimal topological solution [7] B. Tomoiagă , M. Chindriş, A. Sumper, A. Sudria-Andreu, R.
taking into account the set limitations. Algorithm is tested on a Villafafila-Robles, ―Pareto Optimal Reconfiguration of Power
Distribution Systems Using a Genetic Algorithm Based on NSGA-II‖,
part of a network with 213 nodes, 248 lines and 72 switches. Energies 2013, 6, pp.1439-1455, 2013.,
Multi-objective problem is formulated in order to decrease www.mdpi.com/journal/energies
overall network losses and improve system reliability through [8] S.M. Cho, H.S. Shin, J.H. Park, J.C. Kim, ―Distribution System
minimizing ENS. The achieved results show the efficiency of Reconfiguration Considering Customer and DG Reliability Cost‖,
the proposed methodology. Identification of near-optimal Journal of Electrical Engineering & Technology Vol. 7, No. 4, pp.
network configuration is presented. It is evident in decreasing 486~492, 2012.
overall network losses and ENS index compared to network [9] K.Šimleša, Z.Vrančić, A.Matković, F.Lulić, G.Juretić, ―Gubitak
električne energije u distribucijskoj elektroenergetskoj mreži‖, CIRED
initial state. The paper also shows the effect of distributed hrvatski ogranak, 1. savjetovanje, S01-2, Šibenik, Croatia 2008.
generation on voltage profile in distributive network. The
[10] A.Hajizadeh, E.Hajizadeh : PSO-based Planning of Distribution Systems
results show that network reconfiguration in the presence of with Distributed Generators, World Academy of Science, Engineering
DG improve the voltage profile in the network. and Technology, No. 45, pp.598–603, 2008.
It is obvious that application of the proposed methodology [11] F.Rivas-Davalos, E.Moreno-Goytia, G.Gutierrez-Alacaraz, J.Tovar-
Hernandez, ―Evolutionary Multi-Objective Optimization in Power
enables a more complex approach to improving the Systems: State-of-the-Art‖, Proceedings on IEEE PowerTech
operational conditions in distributive networks, compared to Conference, pp.2093-2098, Lausanne, Switzerland, 2007..
traditional methods. The methodology proposed is also a [12] C.M.Fonseca, P.J.Fleming, ―An overview of evolutionary algorithms in
useful tool for quick identification of optimal network multiobjective optimization‖, Evolutionary Computation, pp.1-16, 3(1),
configuration in case of faults, and it can also be beneficial for Spring 1995.
planning and upgrading existing network. [13] K.Deb, A.Pratpat, S.Agarwal, T.Meyarivan, ―A Fast and Elitist Multi-
Objective Genetic Algorithm: NSGA-II‖, IEEE Transaction on
So, both the NSGA II efficiency in finding solutions and Evolutionary Computation, Vol. 6, No2, pp.182-197, 2002.
the increased efficiency of the the distribution network after [14] I.J.Ramirez-Rosado, J.L.Bernal-Agustin, ―Reliability and Costs
using NSGA II are presented in the paper. Optimization for Distribution Networks Expansion Using an
Evolutionary Algorithm‖, IEEE Power Engineering Review, Vol. 21,
The achieved results for both objective functions can be pp. 70-76, 2001
represented in financial terms as well, and these are economic [15] Muhammad Ahmadi, Ashkan Yousefi, Alireza Soroudi, Mehdi Ehsan,
indicators to improve efficiency in managing a distributive ―Multi Objective Distributed Generation Planning Using NSGA-II‖,
network. 13th International Power Electronics and Motion Control Conference
EPE-PEMC, pp.1847-1851, 2008.
11 | P a g e
www.ijarai.thesai.org