IET Research Journals
This paper is a preprint of a paper accepted by IET Renewable Power Generation and is subject
to Institution of Engineering and Technology Copyright. When the final version is published, the
copy of record will be available at the IET Digital Library.

Heuristic Optimization for Automated
Distribution System Planning in Network
Integration Studies
Alexander Scheidler1
1

arXiv:1711.03331v2 [cs.CE] 16 Feb 2018

2

ISSN 1751-8644
doi: 0000000000
www.ietdl.org

, Leon Thurner2 , Martin Braun1,2

Fraunhofer Institute for Energy Economics and Energy System Technology (IEE), 34119 Kassel, Germany
Department of Energy Management and Power System Operation, University of Kassel, 34119 Kassel, Germany
E-mail: alexander.scheidler@iee.fraunhofer.de

Abstract:
Network integration studies try to assess the impact of future developments, such as the increase of Renewable Energy Sources
or the introduction of Smart Grid Technologies, on large-scale network areas. Goals can be to support strategic alignment in the
regulatory framework or to adapt the network planning principles of Distribution System Operators. This study outlines an approach
for the automated distribution system planning that can calculate network reconfiguration, reinforcement and extension plans in
a fully automated fashion. This allows the estimation of the expected cost in massive probabilistic simulations of large numbers
of real networks and constitutes a core component of a framework for large-scale network integration studies. Exemplary case
study results are presented that were performed in cooperation with different major distribution system operators. The case studies
cover the estimation of expected network reinforcement costs, technical and economical assessment of smart grid technologies
and structural network optimisation.

1

Introduction

1.1

Network Integration Studies

In recent years various system-wide network integration studies (e.g.
nation-wide, state-wide, DSO-wide) have been conducted in Germany [1–6]. The goal of these studies is to estimate the expected
costs for the integration of Renewable Energy Sources (RES) in the
distribution system based on different RES installation scenarios and
to assess the technical and economic impact of innovative equipment and control strategies. Such system-wide network integration
studies can support the strategic alignment in the regulatory framework as well as help to adapt the network planning principles of the
Distribution System Operators (DSOs). A focus in these studies is:
• Determination of the maximum RES hosting capacity of the
distribution networks in the considered area.
• Determination of the expected costs for RES integration in the
distribution system.
• Techno-economic assessment of smart grid applications to
increase the RES hosting capacity and to decrease the RES integration costs.
In the medium and low voltage level such studies are often performed on representative networks, which are chosen on the basis of
structural parameters like residents per square kilometer, expected
RES installation or length of feeders. The representative networks
are either built as generic network models [2, 4] or representative
real networks are chosen based on a cluster analysis of the network
parameters [1, 4]. A recent network integration study considers 130
network models of real HV, MV and LV networks [6]. The results
are then extrapolated on a system-wide perspective.
Several studies at the Fraunhofer IEE and the University of Kassel
that investigated a large number of real distribution networks show
a considerable variation in the results and the obtained findings [7].
Especially the technical and economic assessment of different technical measures in the network planning process show a high local

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

or regional diversity which can hardly be covered by characteristic network models. A high sample size is therefore crucial for the
accuracy of the results in RES network integration studies. However,
investigating a large number of networks is only possible with a high
degree of automation. Data import, analysis of the research objective
and export of the results all have to be automated, so that the whole
process can be carried out with minimum manual engagement.
To estimate the cost of RES integration and to be able to economically compare different technologies and planning principles, it
is necessary to calculate the cost of the classic network planning
measures. These measures include network reconfiguration (e.g.
reconfiguration of feeders), reinforcement (e.g. the increase of diameters of lines) and extension (e.g. founding new substations). Some
network integration studies rely on the help of experienced network
planners who plan the necessary network development for selected
representative networks. However, in the scope of automation of network integration studies the calculation of the expected cost must be
carried out fully automatized while it should still generate reasonable realistic solutions. The purpose of this paper is to present our
automated distribution system planning framework that covers both
requirements.
1.2

Automated Distribution System Planing

A novel property of the presented approach for automated distribution system planning is its abstract formulation that allows to apply it
in different studies without great need for adaptation. Its architecture
is based on three parts:
1. A simulation model that defines constraints and allows to evaluate violations of these constraints.
2. A number of reconfiguration, reinforcement and expansion measures that potentially solve these problems.
3. A heuristic optimization that searches for the best combination of
single measures to solve the defined problems.
The heuristic optimization does not rely on specific details of the
considered measures. Instead, the heuristic is only based on the

1

assumption that each measure can either be applied or not applied.
It iteratively searches for better combinations of measures by adding
and removing measures from the current best solution. While this
heuristic approach does not guarantee to find globally optimal solutions, the found solutions are reasonable realistic in the scope of
network studies and fulfil the purpose of estimating the expected
grid reconfiguration, reinforcement and expansion costs. Additionally, the approach is much more flexible and easier to customize than
e.g. analytic optimization methods. For example, it is very easy to
implement and integrate new types of measures, since a mathematical formulation of the measures is not necessary. In that way, even
complex measures, such as splitting a network in two parts or adding
new substations, can be implemented with relatively low effort.
1.3

Structure of the Paper

Section 2 briefly sketches the usual process of conducting automated
network integration studies. Section 3 introduces the network planning problem and Section 4 presents our approach for its solution. In
Section 5 we present findings from four different network integration
studies performed in collaboration with 4 different Distribution System Operators (DSO). Finally, in Section 6 we conclude and discuss
possible future research questions.

2

Automated Network Integration Studies

Network integration studies usually involve one or more DSO(s),
who provide the basic data that is necessary to conduct the study.
Typically, the data comprises a (possibly large) number of real
distribution network models. Moreover, the DSO’s basic network
planing rules and specific constraints are specified. This includes
the definition of different worst case load cases (e.g. high feed-in
case, high load case, n-1 case), the allowed network development
measures and the associated cost assumptions. Additionally, different RES scenarios, i.e. possible future development of RES in the
considered area, and load development scenarios, e.g. possible future
development of electric vehicles, are either provided by the DSO or
developed as part of the study. Finally, the specific research questions
of the study are defined, for example which Smart Grid Technologies
should be investigated.
The network data is mostly provided by the DSO(s) in a commercial network calculation software format like PowerFactoy, Sincal,
NEPLAN, CYME, etc. As a first step we convert all network models
into the pandapower format. pandapower∗ is a Python based
open source framework that is aimed at automation of power system
analysis and optimization in distribution and sub-transmission networks [8]. The software is a joint development of the University of
Kassel and the Fraunhofer IEE. It was hitherto the basis of numerous
network integration studies [7, 9].
The possibly large number of networks in conjunction with the
combinatorial complexity of the input parameters and the simulation
of different probabilistic distributions of RES often leads to several
thousands independent simulation runs. In order to cope with this
high number mostly a HPC cluster is utilized.
A single simulation run typically investigates a specific network
and a specific combination of parameters. That is, a certain random distribution of RES and loads is installed in the network model
and active elements that model Smart Grid components are added.
For the resulting network model it is then checked if violations of
constraints in any of the considered load cases occur. In case of
violations the automated network planning is started to find a combination of reconfiguration, reinforcement and extension measures that
ensures safe operation of the network while being as cost-efficient as
possible. The main result is an estimation of the expected network
reinforcement cost for the specific network.
Finally, the results of all simulation runs are collected and statistically evaluated over the different parameter sets (RES scenarios,
Smart Grid technologies, ect.).

∗ www.uni-kassel.de/go/pandapower

2

3

Automated Network Planning

The automated network planning is a core component of our
approach to conduct large-scale network studies that involve estimating RES integration cost. In the remainder of this paper, we focus on
this component. We consider a single simulation run out of possibly
thousands and assume a readily parameterized network model that is
in a state that violates one or more given constraints. The goal is then
to find a cost-optimal set of network reconfiguration, reinforcement
and expansion measures that would bring the network into a valid
state. This section introduces the network planning problem and its
formalization as a combinatorial optimization problem. The solution
of the problem is discussed in Section 4.
In the following we introduce our approach with the help of two
examples. Both examples address necessity for grid development
after the installation of additional PV plants. However, the examples are only meant as an introduction in the overall framework
and are by no means exhaustive. The versatility of the presented
approach is showcased with the case studies presented in Section
5, which, among others, also address reactive power compensation,
transformer control strategies, expansion of wind power plants, and
decommissioning of lines due to old age. Moreover, the assumed
restrictions (e.g. allowed voltage band or the restriction to radial network structures) only serve illustration purposes and may be very
different in other applications of the approach.
3.1

An Introductory Example

The proposed approach is introduced with the help of an example.
Figure 1a shows a medium voltage ring with 6 lines, 5 MV/LV substations modelled as load points and 10 load-break-switches. The
voltage at the slack node is assumed to be 1.02 p.u. and the upper
voltage limit is 1.05 p.u. Due to an increase in RES the bus voltage
in the considered high feed-in scenario exceeds the voltage limit at
several stations. Figure 1b shows the voltage profile of the network
in this case. The vertical red line in this graph represents the voltage
rise over the transformer. Note that line 4 is not depicted because
the open switch is located on this line. From the network planning
point of view, the network is not valid in its current state. To ensure
its safe operation, the network has to be reconfigured, reinforced or
extended so that all constraints are complied with.
We call every possible action that could be taken by a human network planner an abstract measure. A measure is a single action that
can be applied to the network model and that changes one or several properties of the model. A measure is atomic in the sense that
it can either be applied or not. In case of the example we consider
replacing existing lines by new lines with higher diameter as possible
measures. The set of measures M is then given as:
M = {REPLACE_LINE_1, . . . , REPLACE_LINE_6 }

(1)

We can now treat the planning problem as a combinatorial optimization problem. The overall optimization goal is to find a set
of measures with minimal cost, that leads to a network state that
complies with all technical constraints.
Formally, a solution s is a subset of the available measures s ⊆
M . A solution s corresponds to a network state that is the result of
the application of all measures in s on the initial network state. If
the resulting network complies with all given constraints we call s a
feasible solution.
A solution has associated costs c(s). In our example we define
the cost of a solution as the sum of the costs cm of all individual
measures in the solution:
c(s) =

X

cm

(2)

m∈s

Here we assume that the costs of a replacement measure is equal to
the length of the replaced line.
In order to find the cost optimal network plan for the example we
need to find a feasible solution s∗ ⊆ M that minimizes c(s):

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

Switch 1
(closed)

Line 2

Line 1
8 km

Line 3
5 km

7 km

HV/MV
Line 6

Line 5

Line 4

4 km

9 km

6 km

Switch 6
(open)

with cost c(s∗ ) = 7. As can be seen in Figure 2 b) all voltages are
within the limits after applying s∗ .
This solution does however not consider the possibility of moving the sectioning point to improve the voltage profile, which might
lead to a more efficient solution. We therefore extend the available
measures to include switching operations. We add a new degree of
freedom by allowing to open a different switch than Switch 6. To this
end, we simply use the network model with all switches closed as the
base case and define the opening of switches as additional measures.
The set of possible measures is now given as:

(a) Medium voltage line ring with 6 lines, 5 load points and 10 load-breakswitches

M=

{OPEN_SWITCH_1, . . . , OPEN_SWITCH_10,
REPLACE_LINE_1, . . . , REPLACE_LINE_6 }

(b) Voltage profile; installation of PV plants leads to voltage band violation

Fig. 1: Example network

For these 16 measures there are 216 possible solutions. However, many of these solutions would not lead to feasible network states. For example, the solution s = {OPEN_SWITCH_3,
OPEN_SWITCH_7} would leave 2 nodes unsupplied. Opening no
switch would lead to a non-radial network. To avoid this, radiality
and supply constraints are considered additionally. The radiality constraint lpmf (s) = 0 ensures that the network is not meshed (number
of load points on meshed feeders equals zero). And the supply constraints lpus (s) = 0 ensures that load points are not cut from power
supply (the number of unsupplied load points equals zero). Moreover, we also add the constraint lnol (s) = 0 to avoid overloading
of lines (length of overloaded lines equals zero) and trol (s) = 0 to
avoid overloading of transformers (sum of transformer overloading
equals zero). The optimization goal becomes:
minimize
s∈2M

minimize
s⊆M

subject to

c(s)

subject to

(3)
lpvv (s) = 0,

where lpvv (s) gives the number of load points (lp) where the
voltage limit is violated after the application of s.
replaced

Line 1
8 km

Line 2
7 km

Line 3
5 km

Line 6

Line 5

Line 4

4 km

9 km

6 km

(a) Line 2 is reinforced with a higher diameter cable

(b) Voltage Constraints are complied with after reinforcement

Fig. 2: Optimal reinforcement of the example network
The optimal solution for the example network is found by
testing all 26 possible combinations. This analysis yields s∗ =
{REPLACE_LINE_2} as the most cost efficient feasible solution

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

c(s)
lpus (s)

= 0,

lpmf (s)

= 0,

trol (s)

= 0,

lnol (s)

= 0,

lpvv (s)

= 0,

(5)

The optimal solution for the example that complies with all constraints is s∗ = {OPEN_SWITCH_4, REPLACE_LINE_6}. This
solution has cost c(s∗ ) = 4 and moves the sectioning point from
switch 6 to switch 4 (see Figure 3).
3.2

HV/ MV

(4)

A Realistic Example

In the introductory example the available measures were limited to
exchanging lines and opening switches. In general, the available
measures are more complex and depend on the planning principles of
the DSOs. In the following we present a more realistic example from
one of our case studies. Figure 4a shows the example low voltage
(LV) network after the installation of a large number of additional
PV plants. As can be seen, several buses in the LV network are subject to violations of the voltage criteria in the considered high feed-in
load case. In a first step the automated network planning component identifies the following measures that can potentially solve the
voltage problems:
• Changing the fixed tap position of the transformer. With this measure the voltage of all feeders is lowered by approximately 4 %.
However, it might lead to (more) voltage violations in the high load
case.
• Replacing cables by cables with a higher diameter. The mechanism that discovers line replacement measures only considers lines
on the path from buses with voltage violations to the MV/LV transformer, since only those can effectively influence the voltage drop.
Customer access lines are only considered for reinforcement if they
are longer than 50 m. The resulting options for the example network
are given in Figure 4b.
• Introducing additional parallel lines between two existing switching cabinets. A requirement of the involved DSO was that parallel

3

Switch 4
(open)

Line 1
8 km

Line 2
7 km

Line 3
5 km

replaced

HV/ MV
Line 6

Line 5

Line 4

4 km

9 km

6 km

Switch 6
(closed)

(a) Sectioning point is moved to Switch 3 and Line 6 is reinforced

(b) Voltage Constraints are complied with after reconfiguration and reinforcement

Fig. 3: Optimal reconfiguration and reinforcement of the example
network

lines are only allowed between switching cabinets. Figure 4c shows
the position of the switching cabinets and the parallel line measures.
• Installing new switching cabinets plus a parallel line. In this
case, the additional costs for the switching cabinet are also considered. The possible additional lines with new switch cabinets for the
example network are shown in Figure 4d.
• Splitting the LV network into two networks. The position of all
current switching cabinets further than 50 m from the current substation are considered as possible locations for new substations (see
Figure 4e). The assumed costs include the costs for the substation
itself as well as an estimate of the costs for the MV connection.
• Replacing a transformer by a transformer of larger size. Due to
the changed impedance of the transformer this measure also has
influence on the voltage profile in the network.
Clearly, these kind of detailed measures can only be applied in
models of real networks, since the necessary information is generally
not available in generic network models. When working on realistic
solutions, it is however crucial to consider all these measures and
constraints in order to obtain realistic results.
Even if only measures are considered that are potentially able to
mitigate the voltage problems still a total of 108 possible measures
are identified. In contrast to the first example, it is not possible to
find optimal solutions by brute force testing of all 2108 possible
solutions. Due to the exponential growth of the solution space, an
exhaustive search is generally not a feasible solution technique for
realistic network planning problems.
Figure 4f depicts an example solution that consists of four new
parallel line measures, one new switch cabin and replacing the
transformer. This solution was found using the heuristic solution
approach that is presented in the next section.

4

Heuristic Solution

The outlined network planning problem is a highly constrained,
high-dimensional, mixed integer, non-linear combinatorial optimisation problem. There are numerous studies dedicated to solving
network planning problems with different problem definitions and

4

optimisation approaches [10–14]. One approach is to use deterministic algorithms such as linear programming, non-linear programming
and dynamic programming [15]. While deterministic solvers guarantee to find the global optimum, they usually require simplifications
in the constraints, reduction of the solution space or linearisation
of the problem to create a problem that can be solved analytically.
The realistic example presented in Section 3.2 demonstrates that a
realistic network planning problem is too complex and non-linear
for analytical approaches. Metaheuristics, on the other hand, neither require differentiability, continuity, nor convexity of objective
functions and are efficient in handling constrained, discrete and
multi-modal optimisation problems. They are therefore very popular
in solving distribution system planning problems. Solving the grid
planning problem with a metaheuristic approach was pioneered by
V. Miranda in 1994 [16]. Since then, Different metaheuristics were
successfully applied to solve network planning problems, such as
genetic algorithms [16–21], particle swarm optimisation [22, 23],
tabu search [24, 25], artificial immune systems [26, 27], Iterated
Local Search (ILS) [7] and evolutionary algorithms [28–30]. Heuristic algorithms have also been successfully applied to line routing
problems that take the grid topology and geography into account
[21, 25, 26, 28, 31].
In this paper, we focus on so-called stochastic local search metaheuristics. These algorithms follow a single chain of solutions
through the solution space searching for good solutions [32]. They
usually consist of three parts: a method to determine initial solutions
s0 , a neighbourhood function N (s) that defines how new candidate
solutions are generated from the actual solution s, and an acceptance
criteria that is used to decide if a (randomly) selected neighbour
s0 ∈ N (s) is accepted as the new actual solution.
As initial solution s0 we usually choose the empty set that represents the initial network. Neighbourhood function, cost function and
different algorithms are outlined in the following.
4.1

Neighbourhood Function

Different variants exist for the neighbourhood function. We mostly
define N (s) as all solutions that can be derived from s by either
removing a single measure m ∈ s, adding a single measure m 6∈ s
or exchanging a measure m1 ∈ s with a measure m2 6∈ s. However,
depending on the heuristic and problem it can also be beneficial to
only allow adding and removing but not exchanging measures.
The neighbourhood can be further restricted by defining dependencies between measures. This allows to enforce additional constraints or to reduce the size of the neighbourhood. Possible dependencies are:
• a measure can exclude other measures (e.g., opening only one
switch on the same line section)
• a measure can require other measures (e.g., a parallel line to a
given line is only allowed if this line is already replaced by a line
with maximal diameter)
• at least one of a set of measures must be included in the solution
(e.g., a line has to be exchanged because of its age)
4.2

Cost Function

Most solutions in the solution space are not feasible since they violate one or more constraints. One method to cope with non-feasible
solutions would be to restrict the search process to feasible regions in
the solution space. However, this implies several difficulties. First of
all, one would require means to generate feasible starting solutions.
Yet, this is often a hard task in itself. Moreover, regions of feasible solutions in the solution space do not need to be connected with
respect to the neighbourhood function and consequently the search
might never reach certain regions in the solution space.
Instead of excluding non-feasible solutions, we extend the cost
function with the intend to guide the search process towards feasible
solutions. To this end, the degree of constraint violation is included
in the cost function. More precisely, the cost function is defined as a
tuple of two values, where the first value of this tuple is a number that

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

PV plants
Additional PV plants
Transformer station

Voltage violation

(a) installation of PV plants leads to voltage vio-

Switching cabinet
Possible parallel
lines

possible line
replacement

(b) possible line replacement measures

(c) possible parallel line measures

lation

Possible new
switch cabinets
Possible parallel
lines

(d) possible switch cabin + parallel line measures

New Transformer
Station
Parallel line
New switch cabin

Transformer station
Possible new
transformer stations

(e) possible new station measures

(f) example solution

Fig. 4: Example Low Voltage network; depicted are different discovered measures
represents the violated constraint and the second value represents
the strength of the violation. In case no constraint is violated the
first value is zero and the second gives the cost c(s) of the solution.
The cost of two solutions c0 (s1 ) = (p1 , v1 ) and c0 (s2 ) = (p2 , v2 )
is then compared lexicographically, that is, (p1 , v1 ) < (p2 , v2 ) iff
(p1 < p2 ) ∨ (p1 = p2 ∧ v1 < v2 ).
The extended cost function for the example would then be:

(5, lpus (s)),



(4, lpmf (s)),


(3, trol (s)),
c0 (s) =
(2, lnol (s)),




(1, lpvv (s)),

(0, c(s)),
4.3

if lpus (s) > 0
if lpmf (s) > 0
if trol (s) > 0
if lnol (s) > 0
if lpvv (s) > 0
otherwise

(6)

Stochastic Local Search

Different stochastic local search algorithms can be implemented
based on the same cost and neighbourhood functions.
4.3.1 Hill Climbing (HC): One of the simplest stochastic local
search algorithms is Hill Climbing. Hill Climbing starts with an
initial solution s0 and iteratively moves to a random neighbouring
solution if this step decreases the cost function. This is repeated
until no improving step is available any more. As pseudo-code, the
algorithm can be written as:
1: procedure H ILL C LIMBING(s0 )
2:
s∗ ← s0
3:
repeat
4:
Choose s0 ∈ N (s∗ )
5:
if c(s0 ) < c(s∗ ) then
6:
s∗ = s0
7:
end if
8:
until c(s) ≥ c(s∗ ), ∀s ∈ N (s∗ )
9:
return s∗
10: end procedure

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

4.3.2 Iterated Local Search: Clearly, Hill Climbing will
quickly end up in a local optimum. A valid approach for finding
better solutions is then to just restart Hill Climbing several times.
However, this would "forget" all information collected in the search
process so far. A different approach is to have means for a limited
acceptance of worsening moves. For example, to escape local optima
by slightly perturbing the actual solution. Iterated Local Search (ILS)
does exactly this. In general, the search for better solutions in ILS
occurs in the reduced solution space defined by an arbitrary blackbox heuristic. That is, ILS moves through the solution space from a
local optimum to neighbouring (improved) local optima. We use Hill
Climbing algorithm as local search algorithm for the ILS.
1: procedure I TERATED L OCAL S EARCH(s0 )
2:
s∗ ← HillClimbing(s0 )
3:
while not stopping criteria met do
4:
s0 ← P erturbate(s∗ )
0
5:
s∗ ← HillClimbing(s0 )
0
6:
if c(s∗ ) < c(s∗ ) then
0
7:
s∗ = s∗
8:
end if
9:
end while
10:
return s∗
11: end procedure

To perturbate a solution we simply move strength of arbitrary steps
in the neighbourhood of the solution:
1: procedure PERTURBATE(s0 , strength)
2:
for strength steps do
3:
s0 ← Choose solution ∈ N (s0 )
4:
end for
5:
return s0
6: end procedure

For the stopping criteria different variants can be used. In network
studies we usually stop ILS after a fixed number of iteration without
improvement.

5

(3, 210)

lnol [km]

(2, 0.10)

(1, 0.02)

c [EUR]

ILS_4_HC_AE
ILS_4_HC
LAHC_50

(3, 205)

lpvv [p.u.]

trol [%]

4.3.3 Late Acceptance Hill Climbing: Late Acceptance HillClimbing (LAHC) is a local search algorithm, which accepts nonimproving moves when a candidate cost function is better than it
was a number of iterations before [33]. It has been shown that LAHC
performs comparably to other well established heuristics like Simulated Annealing (SA) [34]. While SA and similar techniques follow a
cooling schedule that includes to calculate the cost difference in case
of worsening moves, this is not necessary in LAHC. This fact plays
well with our definition of the cost function (it is not meaningful to
quantify the difference between the cost values of two solutions that
lead to different constraint violations).

(2, 0.05)

(1, 0.01)
(1, 0)

(0, 10000)
(0, 5000)
0

1000

2000

3000

number of cost function evaluations

4000

(a) Evolution of cost function values of the actual solution for different
Metaheuristics

36

cost [TEUR]

cost [TEUR]

60
55

32
30
28

50
26

HC

_50

_A

E

AE
LA

HC
_
_5_
ILS

_A
E

HC

HC

_50

_A

E

AE

HC
LA

HC
_
ILS

_5_

_A
E
HC

HC

45

(b) Comparison of the distributions of the best found solutions (50 runs,
5000 fitness evaluations) for two different networks

Fig. 5: Comparison of different Metaheuristics

Comparison of Algorithms

Figure 5a depicts how the extended cost function of the actual
solution evolves over time for different metaheuristics. One of our
real network case studies is used to compare the following three
algorithms:
• ILS_4_HC: ILS with HC as local search, perturbation is done with
4 random steps in the neighbourhood. The neighbourhood consist of
adding and removing measures.
• ILS_4_HC_AE: same as ILS_4_HC, only that the neighbourhood
additionally allows exchange (AE) of measures.
• LAHC_50: LAHC with a threshold of 50 previous solutions
The y-axis is divided into 4 areas that represent the 4 different values at the first position in the cost tuple: transformer overloading,
line overloading, voltage violations and cost. For the ILS algorithms,
perturbations are carried out repeatedly until 5000 cost function
evaluations are reached. Between two ILS iterations, the fitness

6

Case Studies

34

65

4.4

5

The automated network planning framework allows flexible combination of measures, constraints and cost functions to reflect different objects of investigation. Within different studies we developed
measure implementations for, e.g.:

75
70

decreases monotonic due to the Hill Climbing local search. The perturbating step on the other hand usually leads to a worsening in the
cost function. In the depicted 5000 fitness evaluations only 3 perturbation steps occur for ILS_4_HC_AE whereas for ILS_4_HC more
than 20 can be observed. This is due to the fact that the neighbourhood of ILS_4_HC_AE is significantly larger than in ILS_4_HC
due to the allowed exchange of measures. The LAHC heuristic
always allows worsening steps and therefore takes more than 3000
cost function evaluations for LAHC to reach a local optimum. It is
restarted until 5000 fitness evaluations are reached.
Which metaheuristic is best suited to solve a problem can only be
decided based on the specific use case. Figure 5b shows results from
two real networks, using the same type of considered measures and
cost function. The only difference is the network, which also has an
effect on the number of measures and violated constraints. It can be
seen that the different heuristics perform differently in each network.
While in one network a restarted Hill Climbing with allowed measure exchange performs on par with the ILS variant, for the other
network the LAHC variant is the best. This shows that there is not
one metaheuristic that clearly outperforms all other algorithms. It
is however possible to draw conclusions from comparing different
algorithms. For example, it seems like the AE algorithms generally outperform the non AE algorithm. At least for this problem,
allowing an exchange of measures seems to pay off, even though it
significantly increases the neighbourhood. Larger sample sizes with
more heterogeneous problems are needed to deduce more general
conclusions about different algorithms.

• replacing existing lines and transformers
• adding parallel lines to existing line trails
• changing the switching configuration
• finding new line trails
• deploying advanced control functions to transformers and PV
systems
• replacing conventional transformers with On-Load-Tap-Changing
(OLTC) transformers
Moreover, in different studies we considered the following constraints:
• radiality, supply, n-1 and other topological constraints
• load flow constraints for bus voltage, line loading, transformer
loading
• several worst-case scenarios, e.g. high load or high generation,
simultaneously
• load flow constraints for n-1 operation with optimal resupply
• reliability constraints for outage times (ASIDI / SAIDI)
In the following we present exemplary results from four network
studies that were performed in cooperation with different major distribution system operators. Due to lack of space many details have
to be neglected. However, we think the case studies allow a general
understanding of what kind questions can be answered with the presented approach. The metaheuristics that are used in the studies are
adapted and tuned versions of the basic variants presented in Section
4. Specific cost values are omitted for confidentiality reasons.
5.1

Expected Network Reinforcement Cost for LV Networks

One goal of this study was to estimate the expected network reinforcement cost for different RES scenarios. The study was done in
cooperation with the German DSO EnergieNetz Mitte.

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

Figure 6 shows results for two different PV scenarios. The distribution of the cost values for the single LV networks is the result
of the simulation of 50 different probabilistic distributions of PV
plants within each network. As can be seen the different networks
have different sensitivities to the specific distributions of PV plants.
The results show moreover that for the conservative scenario half
of the investigated 67 LV networks (with a total of 570 feeders)
need reinforcement measures. In the progressive scenario more than
two thirds of the networks are affected and the expected costs per
network are higher.

PV Scenario

Reinforcement Cost in EUR

Conservative
Progressive
(a) Comparison of the expected network reinforcement cost using different
Smart Grid technologies for a MV network of the swiss DSO Romande
Energy

0

67 Low Voltage Networks

Fig. 6: Example results of the expected reinforcement cost for 2
different RES scenarios

(b) Comparison of avoided network reinforcement costs and OLTC invest-

5.2

Inverter Control Strategies in MV Networks

Within the scope of a PV integration study carried out for the Swiss
distribution system operator Romande Energie, the hosting capacity
of 111 LV networks and two MV networks was evaluated and compared to different PV forecast scenarios provided by the DSO (for
a definition of the term hosting capacity see e.g. [7]). The expected
PV capacity in 2035 exceeds the determined PV hosting capacity of
the LV networks only in few cases. Therefore, the LV networks are
mostly well dimensioned for the expected additional PV installations
and additional measures for PV integration are just expected in a few
LV networks.
The analysis of the PV hosting capacity of one MV network
showed need for network reinforcement. It was therefore further
investigated how control strategies can mitigate the expected reinforcement cost. The considered control strategies are:
• Constant CosPhi: reactive power provision by PV plants with a
constant power factor
• CosPhi(P): reactive power provision by PV plants with a power
factor depending on active power provision
• Peak Shaving: active power curtailment of PV plants
• AOLTC: advanced on-load tap changer control of the HV/MV
transformer, where the voltage set point of the transformer control
is adapted to the active power flow over the transformer
Exemplary results are shown in Figure 7a. The distribution of cost
values is the result of the simulation of different distributions of PV
systems in the networks.
It is clearly visible that all studied autonomous voltage control
strategies are capable of reducing the expected reinforcement costs
considerably. Due to the specific structure of the MV network mainly
voltage rise is the limiting factor. Therefore, increasing the allowed
voltage band by using an AOLTC has a significant effect on the
expected reinforcement costs. It can reduce the PV integration cost
by 90% in 2035 compared with no strategy on average.

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

ment costs for 84 LV networks of the German DSO Bayernwerk

Fig. 7: Example results for the impact of Smart Grid technologies
on reinforcement cost

5.3 Technical and Economical Assessment of MV/LV OLTC
transformer
In LV networks the application of an OLTC for the MV/LV transformer is a promising technical solution to increase the PV hosting
capacity. In a network integration study in cooperation with the German DSO Bayernwerk the technical and economic potential of the
OLTC was analysed for 84 real LV networks. The required network
reinforcement costs for PV integration are determined with and without OLTC. Figure 7b shows a comparison of the avoided network
reinforcement costs by OLTC application and the OLTC investment
costs for the LV networks. The OLTC is a cost-effective measure for
PV integration for more than two thirds of the LV networks. It should
be noted, that a high PV penetration scenario was investigated (PV
rooftop potential) and that the operational costs are not considered
in this example.
5.4

Network Topology Optimisation

A goal in strategic network planning can be to renew specific network elements, for example elements of a particular type (e.g.
overhead lines or cables with a specific insulation) or elements that
will have exceeded their life expectancy at the planning horizon.
Consequently, sometimes large parts of a network area are to be
renewed. In such a case it can be more cost-efficient to find a network
structure with new line trails than simply maintaining the existing
structure by renewing old line trails.
Figure 8 shows an example of the automated topology optimisation in a 10 kV network of the DSO Westnetz GmbH that has

7

been developed in the scope of the project ANaPlan. The lines that
are to be renewed in the target network are shown as dashed lines
in Figure 8a. Figure 8b shows the possible new line trails that are
considered in the optimization. Each line trail is modelled as one
network reinforcement measure as explained in Section 3. The line
length of a new trail is assumed to be equal to the airline distance
multiplied with a factor of 1.5 to account for obstacles. Additionally, switch measures are considered in the optimisation to allow
for a reconfiguration of the feeder partitioning. In this example, 232
line trail measures and 605 switching measures have been considered. With this kind of complexity, it is important to use problem
specific knowledge (like switches on stubs can never be opened,
two switches on the same feeder sections can never be opened at
the same time etc.) to restrict the problem. The automated network
optimisation is then used to find a solution that requires a minimum
amount of cabling while complying with all defined constraints. In
this example, topological and operational constraints are used. The
topological constraints ensure the radial structure of the network and
possibilities of resupply in case of line faults. The operational constraints ensure the compliance with voltage band and maximum line
loading in worst-case situations. The worst-case is modelled by one
low load and high RES scenario and one high load and low RES
scenario. By applying a prognosis for load development and RES
installation at the planning horizon, the methodology can be used to
find cost-efficient network structures for future power systems. The
best solution that is found by the optimisation can be seen in Figure
8c. The feeder configuration in Figure 8d shows that the radiality of
the network is maintained. The optimized network structure leads to
about 25 km of new line trails, while a renewal of the old line trails
would have resulted in about 31 km. The optimisation was therefore
able to find a more efficient structure that complies with topological
and operational constraints considering load and RES development.

6

This paper introduces an approach for automatic network planning
that is a core component of our framework for large-scale network
studies. It allows the calculation of network reconfiguration, reinforcement and extension with detailed network models. Since no
data reduction or simplification is necessary, results of the automated network planning can be directly compared to solutions of
experienced network planners to validate the results and improve the
algorithm. This increases the transparency of RES network integration studies and permits direct conclusions regarding the network
planning principles of the DSO. A further benefit of the automated approach is the possible application of probabilistic load and
generation scenarios in the network integration study, which further increases the robustness and informative values. The network
integration studies can be performed for a large number of real
distribution networks, which avoids inaccuracies caused by the simplification of characteristic network models and their projection on
a system-wide perspective. Therefore, the approach covers the complete diversity of the investigated distribution networks and allows
conclusions for the network planning process of the individual network sections. The case studies presented in this paper highlight the
need for detailed investigations, since results vary greatly in different
case studies with different DSOs.
There are two main challenges for automated network studies.
The first is the availability of data: quality and the informative value
of the RES network integration studies are strongly dependent on
the provided data base. Since a detailed knowledge of network conditions will become increasingly important in active distribution
networks, most DSOs are actively working on improving data maintenance and standardisation to facilitate automated data analysis. An
improvement in data quality and availability is thus to be expected.
Another challenge is the complexity of the network planning optimisation problem, which rises exponentially with the number of
considered measures. Metaheuristic optimisation has been shown to
find good solutions even for heterogeneous problems with several
hundreds of measures. The presented algorithms have been validated
and improved with the feedback of network planning experts of several DSOs. In future, the algorithm could not only be used for studies
but also applied as a direct supporting tool in the network planning
process.

7

(a) Initial network with lines

(b) Line trails considered for new net-

to remove (dotted)

work structure

1
2
3
4
5
6
7

(d) Feeder sectioning in optimized net-

structure with new line trails

work structure

8

Fig. 8: Case study for topological network optimization
9

8

Acknowledgment

This research was partly supported by the German Federal Ministry
for Economic Affairs and Energy and the Projektträger Jülich GmbH
(PTJ) within the framework of the projects SmartGridModels (FKZ
0325616) and ANaPlan (FKZ 0325923B); and by the Federal Ministry of Education and Research (BMBF) within the framework of
the project ENSURE (FKZ 03SFK1N0).

8

(c) Optimized network

Conclusion

References
Deutsche Energie-Agentur GmbH (dena): ‘Ausbau- und Innovationsbedarf der
Stromverteilnetze in Deutschland bis 2030’, , 2012,
Bundesministerium für Wirtschaft und Energie (BMWi): ‘Moderne Verteilernetze
für Deutschland (Verteilernetzstudie)’, , September 2014,
Witzmann, R., et al.: ‘Verteilnetzstudie Bayern 2013 (Zwischenbericht)’, , August
2013,
Ackermann, T., Martensen, N., Brown, T., Untsch, S., Tröster, E., Geidel, S., et al.:
‘Verteilnetzstudie Rheinland-Pfalz’, , Januar 2014,
Rehtanz, C., Moser, A., et al.: ‘Leistungsfähigkeit und Ausbaubedarf der Verteilnetze in Nordrhein-Westfalen (Gutachten)’, , August 2014,
Rehtanz, C., Greve, M., et al.: ‘Verteilnetzstudie für das Land BadenWürttemberg’, , April 2017,
Scheidler, A., Thurner, L., Kraiczy, M., Braun, M. ‘Automated grid planning for
distribution grids with increasing pv penetration’. In: 6th International Workshop
on Integration of Solar Power into Power Systems. (Vienna, Austria, 2016.
Thurner, L., Scheidler, A., Schäfer, F., Menke, J.H., Dollichon, J., Meier, F., et al..
‘pandapower - an Open Source Python Tool for Convenient Modeling, Analysis
and Optimization of Electric Power Systems’. (, 2017. preprint. Available from:
https://arxiv.org/abs/1709.06743
Thurner, L., Scheidler, A., Probst, A., Braun, M.: ‘Heuristic optimization for network restoration and expansion in compliance with the single contingency policy’,
IET Generation, Transmission & Distribution, 2017,

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

Georgilakis, P.S., Hatziargyriou, N.D.: ‘A review of power distribution planning
in the modern power systems era: Models, methods and future research’,
Electric Power Systems Research, 2015, 121, pp. 89 – 100.
Available
from: http://www.sciencedirect.com/science/article/pii/
S0378779614004490
Jordehi, A.R.: ‘Optimisation of electric distribution systems: A review’, Renewable and Sustainable Energy Reviews, 2015, 51, pp. 1088 – 1100. Available
from: http://www.sciencedirect.com/science/article/pii/
S1364032115006516
Sedghi, M., Ahmadian, A., Aliakbar.Golkar, M.: ‘Assessment of optimization
algorithms capability in distribution network planning: Review, comparison and
modification techniques’, Renewable and Sustainable Energy Reviews, 2016,
66, pp. 415 – 434. Available from: http://www.sciencedirect.com/
science/article/pii/S1364032116304592
Zubo, R.H.A., Mokryani, G., Rajamani, H.S., Aghaei, J., Niknam, T.,
Pillai, P.: ‘Operation and planning of distribution networks with integration of renewable distributed generators considering uncertainties: A
review’, Renewable and Sustainable Energy Reviews, 2017, 72, pp. 1177 –
1198. Available from: http://www.sciencedirect.com/science/
article/pii/S1364032116306840
Ganguly, S., Sahoo, N.C., Das, D.: ‘Recent advances on power distribution system
planning: a state-of-the-art survey’, Energy Systems, 2013, 4, (2), pp. 165–193.
Available from: http://dx.doi.org/10.1007/s12667-012-0073-x
Franco, J.F., Rider, M.J., Romero, R.: ‘A mixed-integer quadratically-constrained
programming model for the distribution system expansion planning’, International Journal of Electrical Power & Energy Systems, 2014, 62, pp. 265
– 272. Available from: http://www.sciencedirect.com/science/
article/pii/S0142061514002439
Miranda, V., Ranito, J.V., Proenca, L.M.: ‘Genetic algorithms in optimal multistage
distribution network planning’, IEEE Transactions on Power Systems, 1994, 9, (4),
pp. 1927–1933
Chen, T.H., Cherng, J.T.: ‘Optimal phase arrangement of distribution transformers
connected to a primary feeder for system unbalance improvement and loss reduction using a genetic algorithm’, IEEE Transactions on Power Systems, 2000, 15,
(3), pp. 994–1000
Camargo, V., Lavorato, M., Romero, R. ‘Specialized genetic algorithm to solve the
electrical distribution system expansion planning’. In: 2013 IEEE Power Energy
Society General Meeting. (Vancouver, Canada, 2013. pp. 1–5
Fletcher, J., Fernando, T., Iu, H., Reynolds, M., Fani, S. ‘A case study on optimizing an electrical distribution network using a genetic algorithm’. In: 2015 IEEE
24th International Symposium on Industrial Electronics (ISIE). (, 2015. pp. 20–25
Falaghi, H., Singh, C., Haghifam, M.R., Ramezani, M.: ‘Dg integrated multistage
distribution system expansion planning’, International Journal of Electrical Power
& Energy Systems, 2011, 33, (8), pp. 1489 – 1497
Kong, T., Cheng, H., Hu, Z., Yao, L.: ‘Multiobjective planning of openloop mv distribution networks using comgis network analysis and moga’,
Electric Power Systems Research, 2009, 79, (2), pp. 390 – 398. Available
from: http://www.sciencedirect.com/science/article/pii/
S0378779608002204
Sahoo, N.C., Ganguly, S., Das, D.: ‘Multi-objective planning of electrical distribution systems incorporating sectionalizing switches and tie-lines using particle
swarm optimization’, Swarm and Evolutionary Computation, 2012, 3, pp. 15 – 32
Ganguly, S., Sahoo, N., Das, D.: ‘A novel multi-objective pso for electrical distribution system planning incorporating distributed generation’, Energy Systems,
2010, 1, (3), pp. 291–337
Cossi, A.M., Silva, L.G.W.D., Lazaro, R.A.R., Mantovani, J.R.S.: ‘Primary power
distribution systems planning taking into account reliability, operation and expansion costs’, IET Generation, Transmission Distribution, 2012, 6, (3), pp. 274–284
Navarro, A., Rudnick, H.: ‘Large-scale distribution planning - part ii: Macrooptimization with voronoi’s diagram and tabu search’, IEEE Transactions on
Power Systems, 2009, 24, (2), pp. 752–758
Keko, H., Skok, M., Skrlec, D. ‘Solving the distribution network routing problem
with artificial immune systems’. In: Proceedings of the 12th IEEE Mediterranean Electrotechnical Conference (MELECON). vol. 3. (Dubrovnik, Croatia,
2004. pp. 959–962 Vol.3
Lin, C.H., Chen, C.S., Huang, M.Y., Chuang, H.J., Kang, M.S., Ho, C.Y., et al.
‘Optimal phase arrangement of distribution feeders using immune algorithm’.
In: 2007 International Conference on Intelligent Systems Applications to Power
Systems. (, 2007. pp. 1–6
Zmijarević, Z., Skok, M., Keko, H., Škrlec, D., Krajcar, S., Lang.Kosić, N., et al. ‘A
comprehensive methodology for long-term planning of distribution networks with
intrinsic contingency support’. In: 18th International Conference and Exhibition
on Electricity Distribution, CIRED. (Turin, Italy, 2005. pp. 1–5
Cossi, A.M., Romero, R., Mantovani, J.R.S.: ‘Planning of secondary distribution
circuits through evolutionary algorithms’, IEEE Transactions on Power Delivery,
2005, 20, (1), pp. 205–213
Gitizadeh, M., Vahed, A.A., Aghaei, J.: ‘Multistage distribution system
expansion planning considering distributed generation using hybrid evolutionary algorithms’, Applied Energy, 2013, 101, (Supplement C), pp. 655
– 666. sustainable Development of Energy, Water and Environment Systems. Available from: http://www.sciencedirect.com/science/
article/pii/S0306261912005223
Domingo, C.M., an A..Sanchez.Miralles, T.G.S.R., Gonzalez, J.P.P., Martinez,
A.C.: ‘A reference network model for large-scale distribution planning with automatic street map generation’, IEEE Transactions on Power Systems, 2011, 26, (1),
pp. 190–197
Lourenço, H.R., Martin, O.C., Stützle, T. In: Glover, F., Kochenberger, G.A.,
editors. ‘Iterated local search’. (Boston, USA: Springer US, 2003. pp. 320–353.
Available from: http://dx.doi.org/10.1007/0-306-48056-5_11

IET Research Journals, pp. 1–9
c The Institution of Engineering and Technology 2015

33

34

Burke, E.K., Bykov, Y.: ‘The late acceptance hill-climbing heuristic’,
European Journal of Operational Research, 2017, 258, (1), pp. 70 –
78.
Available from: http://www.sciencedirect.com/science/
article/pii/S0377221716305495
Kirkpatrick, S., Gelatt, C.D., Vecchi, M.P.: ‘Optimization by simulated annealing’,
SCIENCE, 1983, 220, (4598), pp. 671–680

9

