Hosting Capacity For Smart Power Grids
Hosting Capacity For Smart Power Grids
Hosting Capacity For Smart Power Grids
F. Zobaa
Shady H.E. Abdel Aleem
Sherif M. Ismael
Paulo F. Ribeiro Editors
Hosting
Capacity
for Smart
Power Grids
Hosting Capacity for Smart Power Grids
Ahmed F. Zobaa • Shady H. E. Abdel Aleem
Sherif M. Ismael • Paulo F. Ribeiro
Editors
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface
Today’s power systems are becoming much more complex, expandable, dynamic,
and non-predictable; therefore, accurate and up-to-date hosting capacity (HC) anal-
ysis studies are necessary to face the significant challenges associated with the
massive integration of renewable energy resources into modern electrical power
systems. This book is designed to be essentially an introduction to the concepts and
historical developments of the hosting capacity assessment. It will be beneficial to
distribution system operators, network planners, distribution generation investors,
and researchers in this field as it presents past, present, and future possibilities
regarding HC research. It is expected that the HC will play a significant role in
future power systems and smart grids. Using the HC approach to drive network
requirements could steer distributed generation (DG) toward areas of the network
where it could have the greatest positive impact on network reliability and win-win
benefits with DG owners. Finally, smarter DG integration into future electrical
systems can be met if utilities have a clear forecast of their potential network
HC. As can be seen, the chapters in each section maintain their own thematic
continuity and at the same time have significant overlaps with chapters in other
sections. Therefore, one may read the book in its entirety or focus on individual
chapters. Due to its broad scope, this book will be an ideal resource for students in
advanced graduate-level courses and special topics in the field of smart grids, hosting
capacity assessment, and enhancement in modern electrical power systems. HC
research is a key enabler for affordable, reliable, and renewable energy sources, so
it is possible to transition away from traditional high-carbon energy sources. There-
fore, it is imperative that novel solutions be sought to enable networks to cope with
future developments to realize resilient distribution networks that can host the pushy
DG penetration while enhancing the system reliability of power supplies and
controlling the over-hosted areas. Based on the editors’ experiences, we believe
that this book is an important guide to the researchers, distribution system operators,
and network planners who are interested in the integration of renewable energy
resources in the electrical power systems. Comprehensive presentation of the hosting
v
vi Preface
perspectives. It is a useful tool for the network planners, designers, operators, and
practicing engineers of modern power systems who are concerned with renewable
energy resources and their challenges, reliability, and security. Likewise, it is a key
resource for advanced students, postgraduates, academics, and researchers who had
some background in electrical power systems.
In an experiment in criticism [1], the author suggests that books allow us to see
the world through other people’s eyes and view it as someone else sees it—and
experience their view “from the inside.” This is what gives books a profound and
mysterious power. At the same time, each of us sees the world through our own
perspective. Thus, the editors wish you may find this book useful to stimulate your
imagination and go even further by applying the concept of HC, but doing it with
your own sensitive and critical eyes.
Reference
ix
x Contents
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 249
Chapter 1
An Overview of Hosting Capacity
for Modern Power Grids
1.1 Introduction
to various technical, economic, and environmental factors and political and social
initiatives [1–3].
However, the booming integration of distributed generation (DG) units, if not
properly assessed, may lead to various problems such as over- and under-voltages,
excessive line losses, overloading of transformers and feeders, protection failure, and
high harmonic distortion levels exceeding the limits of international standards.
These problems occur when the number of DG units exceeds the maximum permis-
sible penetration level, in short when the system exceeds its hosting capacity
(HC) limit. The HC concept is essential to evaluate the system capacity to host
these DG newcomers without exceeding the allowable operational performance.
Practically speaking, a conflict of interest has been found among the DG owners/
investors and distribution system operators (DSOs) in deregulated energy markets. It
is because the DG investors are looking forward to more and more DG integration
into electrical networks, while DSOs are concerned about excessive DG penetration
problems. Therefore, the HC concept was proposed to resolve this conflict, as it
provides a fair and transparent solution that clearly indicates when to accept or to
reject new DG integration requests [4–6].
Today’s power systems are becoming much more complex, expandable,
dynamic, and non-predictable; therefore, accurate to-date HC analysis studies are
necessary to face the significant challenges associated with the high penetration of
renewable energy resources into modern electrical power systems.
HC is a key enabler for affordable and reliable renewable energy sources; hence,
it is possible to transition away from traditional high-carbon energy sources. It is
expected that the HC will play a significant role in future power systems and smart
grids. Using an HC approach to drive network requirements could steer DG toward
areas of the network where it could have the greatest positive impact on network
reliability and win-win benefits with DG owners. Thus, it is imperative that novel
solutions be sought to enable networks to cope with future developments to realize
resilient distribution networks that can host the pushy DG penetration while enhanc-
ing the system reliability of power supplies and controlling the over-hosted areas,
that is, the areas that have problems because of violation of technical performance
metrics.
Not long ago, quick but conservative methods were used to limit integrated DG
capacities. Nowadays, the HC concept has been developed to define the maximum
amount of DG that can be integrated into the electrical network accurately, without
exceeding any operational performance limits. In addition, the HC is a specific,
measurable, practical, and fair concept that facilitates the ease of use of clear
performance limits as evaluation criteria for DG penetration.
Currently, enhancing the system HC is considered one of the important goals for
DSOs, worldwide [1]. In the literature, common HC enhancement techniques were
investigated, namely, reactive power control, automatic voltage control techniques
such as on-load tap changer (OLTC) transformers, active power curtailment, net-
work reconfiguration and reinforcement, and harmonic mitigation techniques.
The authors believe that along with these HC enhancement techniques, energy
storage (ES) technologies will play a significant role in the near future as a promising
1 An Overview of Hosting Capacity for Modern Power Grids 3
In [7–14], HC issues related to voltage regulation and reverse power flow are
presented. However, it is necessary to consider also the time-varying impact of
harmonics on the voltage regulation in the presence of renewables. This poses a
problem as HC conditions vary in time. The main limitation of determining the HC is
its reliability, particularly when the system is under constant harmonic injection due
to the aggregated electronic devices connected to the electrical system, thus making
it a dynamic harmonic injection issue. One question that needs to be raised is how a
dynamic harmonic injection could affect the global HC value due to its impact on the
voltage profile.
A related hypothesis maintains that a dynamic harmonic injection, originating
from external conditions, such as the background harmonic injection, as well as
distortion from the inverters (e.g., in photovoltaic (PV) systems), leads to a DHC
profile which follows the tendency of their harmonic injection. This suggests that it
presents a constructive curve in time, showing that harmonic distortion problems
should be considered when planning improvement strategies to be applied in
electrical distribution systems.
The fundamental characteristic of DHC can be described as the allowance to
increase and expand the usual approach of HC by applying a time variance, which
can designate its behavior against both external and internal variations regarding the
system, such as load conditions, background harmonic distortion, and irradiation
index, among others [15, 16]. Furthermore, in accordance with [15], the HC at the
point of common coupling (PCC) can be defined as the sum of maximum power
possible to be injected for each frequency. Furthermore, time variation is set as a
vital factor in (1.1) within this analysis.
4 T. E. Castelo de Oliveira et al.
Dynamic HC
[9] Maximum HC
Acceptable region
Minimum HC
Time
X
h Xh V max ,h V max ,h V h ðt Þ
1 g g o
g ðt Þ ¼
Pmax max ,h
PPV ðt Þ ¼ ð1:1Þ
h¼1
Rf h¼1 1 þ tan ðφÞ hX f
Rf
Distance
'V
Distance
Conventional electrical networks consider the total load flow directly from the
substation to the load. In regard to this fact, there is a voltage drop throughout the
transmission line, as shown in Fig. 1.2. This unidirectional load flow allows for easy
planning of the distribution transformers and the conductor section of the
feeders [17].
When there is a high penetration of DER in the electrical network, there may be a
reversal and/or bidirectionality of the power flow. As an example, the DER absorbs
the consumption of the load to which it is connected, and the additional power will
be injected into the distribution network, feeding other loads nearby. The phenom-
enon is shown in Fig. 1.3. Due to the bidirectional power flow, some negative
aspects can be registered, for example, a voltage rise at the connection point where
the DERs are installed, which can exceed standard indices, leading to possible
losses [17].
Alternatively, storage systems could help the system at the connection point by
absorbing excessive energy generation of renewables. Several researchers have
called this into question, claiming that storage systems can lead to improvements
into the distribution system, such as avoiding the voltage rises due to the bidirec-
tional power flow, as well as improving the HC of the system. Figure 1.4 shows the
storage system’s impact at the connection point. It is possible to understand that the
voltage would not rise in the system because there is no bidirectional load flow into
the system.
6 T. E. Castelo de Oliveira et al.
Time
This problem can be outlined in terms of the power of the system. In this case, the
power’s sum will be zero due to the relation Pg (Ps + Ps), where Pg is the power
generated by DERs, Pc is the power consumed by loads, and Ps is the power which
will be stored by the storage system. It is important to highlight that the sum will be
zero, as long as there is a generation from the DERs, for example, in regard to PVs,
only at the moment when there is solar irradiation. This allows for a formal solution
to be found, where ΔV 0. In theory, the lesser the voltage rise, the better the HC
will be.
It is easily verifiable that there are two possible options in order to bring
improvement at the connection point: either regulate the voltage through regulators
such as OLTC or connect some storage systems into the connection point. It is
necessary to analyze which option is better economically speaking. In this analysis,
the storage systems will be under investigation with the intention of bringing
improvement regarding the DHC profile for the system, as shown in Fig. 1.5.
This proposed hypothesis examines how to solve the voltage rise problem in
electrical grids using energy storage systems in order to improve the DHC approach.
The ES system is examined with the intention of finding the best contribution to
decrease the voltage rises, as well as keeping at zero the reverse power flow, thus
improving the hosting capacity. This is a feasible but expensive solution for distri-
bution systems in relation to voltage rise indices due to DERs connected at the
electrical grid.
1 An Overview of Hosting Capacity for Modern Power Grids 7
1 V o ðt ÞðV o ðt Þ 1, 0Þ
ΔPmax
storage ðt Þ ¼
ð1:2Þ
Rf 1 þ tan ðφÞ XRff
Pgmax 0 ðt Þ ¼ Pmax
g ðt Þ þ ΔPstorage ðt Þ
max
ð1:3Þ
max 0
DSTHCtotal ðt Þ ¼ HCmax
total ðt Þ þ ΔHCstorage ðt Þ
max
ð1:4Þ
The new voltage profile for the system after the energy storage system can be
verified. There will not be a voltage rise. The voltage drop can be determined by
(1.2). Thus, the new voltage value is given in (1.5).
V 0o ðt Þ ¼ V o ðt Þ ΔV storage ðt Þ ð1:5Þ
HC research is a key enabler for affordable, reliable, and renewable energy sources,
so it is possible to transition away from traditional high-carbon energy sources.
Therefore, it is imperative that novel solutions be sought to enable networks to cope
with future developments to realize resilient distribution networks that can host the
pushy DG penetration while enhancing the system reliability of power supplies and
controlling the over-hosted areas.
This chapter presents an overview of the HC concept and its developments,
defines the main HC performance limits, and its up-to-date enhancement techniques.
Then, the chapter presents the concept of dynamic HC as well as its application to
energy storage systems. The authors presented a new concept for calculating the
local HC as a contribution to the context of planning and improving DER installa-
tions into distribution systems. The definition of DHC has been discussed, as well as
its application to energy storage systems.
References
1. Ismael, S. M., Aleem, S. H. E., Abdelaziz, A. Y., & Zobaa, A. F. (2019). State-of-the-art of
hosting capacity in modern power systems with distributed generation. Renewable Energy, 130,
1002–1020. https://doi.org/10.1016/j.renene.2018.07.008.
2. Ismael, S. M., Aleem, S. H. E., Abdelaziz, A. Y., & Zobaa, A. F. (2018). Practical consider-
ations for optimal conductor reinforcement and hosting capacity enhancement in radial
1 An Overview of Hosting Capacity for Modern Power Grids 9
Hassan Al-Saadi
2.1 Introduction
H. Al-Saadi (*)
School of Electrical and Electronic Engineering, University of Adelaide, Adelaide, SA,
Australia
e-mail: hassan.al-saadi@adelaide.edu.au
specifically could occur in the mid of the feeder [2]. The problems mentioned above,
and others not mentioned here, give rise to a new expression in the taxonomy of
power system called “hosting capacity” (HC). While the general definition may
simply be referred to a maximum limit of a power network in hosting PVs without
jeopardizing connected appliances, the determination of HC is still an attractive topic
that has not reached the consensus yet. It has been found that the definitions of HC
seem to converge in their implication in the published literature. Examples of these
definitions are as follows:
The hosting capacity of the LV grid for dispersed generation is restricted by the maximum
permissible voltage rise within the grid and the maximum short-term loading of the
transformer and the cables, due to the diurnal cycle of the PV [3]
The maximum amount of new production or consumption that can be connected without
endangering the reliability or quality for other customers” and “the acceptable degree of
DER1 penetration under given circumstances [4]
The upper limit of DG before network congestion occurs [5]
The largest PV generation that can be accommodated without violating the feeder’s oper-
ational limits [6]
In this chapter, the adopted definition of hosting capacity (HC) is the maximum
capacity of DGs connected in a distribution network without disrupting the normal
operational conditions of that network.
1
DER is the abbreviation of distributed energy resource.
Table 2.1 The summary on the existing different methods for HC studies
Case study Criteria HC and penetration method Grid type
[7] Two buses with PV connected into Voltage and line current limits Analyze the steady-state bus voltage and line current Bidirectional
one bus and the another is substa- behavior when PV penetration increases with different sce-
tion bus with constant voltage in narios such as different PFs, percentages of load penetration,
addition to the assumption of loads and substation voltages
summed up and connected into the
PV bus
[8] Use an actual and large distribu- For HC purposes, ANSI C84.1 Create different scenarios wherein the PV penetration, at Unidirectional
tion feeder in California with standards as overvoltage criteria each scenario, is increased by the insertion of the additional
commercial and residential urban were taken. Other criteria taken PVs on the randomly selected locations
and rural feeders, modelled using for other analyses such as fre- The PVs inserted are selected from a set of PV sizes (43,502
OpenDSS quency distortion residential preferred PVs, 2625 commercial preferred PVs
which are based on SCADA and AMI, provided from Cal-
ifornia Solar Initiative Survey) in which that set is probabi-
listically distributed according to customer services of PV
sellers
[9] 302 low-voltage rural and subur- 3% voltage variation, no Classification approach is about classifying the HC capabil- Unidirectional
ban grids overloading electrical ity of the grids into five categories (very weak up to very
components strong) using Weibull fitting model
HC capability of a grid is probabilistically treated using
Weibull distribution as well after evaluating the single grid
with a number of DG random sizes with different configu-
rations
DG penetration is performed in piecewise increase with
discrete interval 0.1 kW
(continued)
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . .
13
Table 2.1 (continued)
14
exporting power back to the bulk grid through grid supply point. Each referenced
method is intensively reviewed and analyzed within this table.
The method that grabbed attention most and is adopted in the current work is
introduced by Bollen and Hassan [4]. They explain their concept for HC determi-
nation and introduced the general fundamentals to seize the impacts of DGs. It is
stated that “The impact of distributed generation can be quantified only by using a set
of performance indicators.” The further explanation of this concept is provided in
[13]. The summarized steps are as follows:
First stage Establish a number of performance indices.
Second stage Specify a number of standardized permissible limits for the
established indices in the 1st stage such as a limit can be defined
by EN50160 [14] or any other standards defined by the local utility.
Third stage Find the mathematical relationship between the indices and the
increases in the number of DG connections.
Fourth stage Identify network HC by streamlining the different states (acceptable
deterioration, unacceptable deterioration or critical deterioration)
according to the exceedance of the limits in second stage.
The last mentioned method is depicted in Fig. 2.1. The x-axis displays DG
penetration increases starting from a prespecified amount of power penetration.
The y-axis displays the degree of deterioration accompanied with these increases.
The region shaded in light green indicates that the system is under tolerable
deterioration. The region shaded in red is where the system is considered under
intolerable deterioration. The network HC limit is identified once the deterioration
index crosses from one region into another. There is a possibility of having critical
region that could deliver two different values of HC limits. This region is shown in
the figure with partially yellowish color. In this chapter, this method is promoted,
thereby establishing new risk indices for operational performance which is intro-
duced in the next section.
One of the deterring obstacles in the topics of risk assessment in power system is the
analytical expression of relative frequency of occurrence. The number and different
sources of uncertainties in the network make the task of finding analytical represen-
tation non-trivial, especially in linking the input variables to the response variable.
Instead of the analytical solution, the power of numerical computing techniques is
proven to provide feasible answers for many engineering problems. Therefore, in
this chapter, the nonintrusive method is employed for a set of scenarios. For each
scenario, the numerical simulation is adopted for the evaluation of the system
through probabilistic distributions and statistical analysis.
To formulate the problem mathematically, two quantities are considered: uncer-
tain quantities and fixed ones, noted by x and d, respectively. While the uncertain
quantities, x, stand for the intermittent output power of the connected PV or loads,
the fixed quantities, d, represent the level of PV increase connection and the time of
the day. In Fig. 2.2, two stages are essential in using nonintrusive method: formu-
lating the uncertainty, x, and propagating in the form of quantity of interest.
The procedures of system evaluation can be explained in three steps.
and maximum values, noted kt min and kt max, respectively. The more tendency of this
quantity to the maximum values means having more cloudless atmosphere. With the
assumption that each PV is equipped with maximum power point tracker, the output
power of PV as a function of clearness index, PPV(kt), is provided as follows:
where η is the efficiency of the entire PV unit including the aging degradation; A is
the surface area; and I βt ðkt Þ is the radiation incident on the surface with the incline
angle, β. The latest is a function of clearness index which can be expressed as
follows [15]:
8
>
> D ∙ k t þ D0
kt
D > 0&D0 0
>
> 0
< 1 þ eBðB kt Þ
β
I t ðk t Þ ¼ kt ð2:2Þ
>
> D ∙ k t D0 0 D > 0&D0 < 0
>
> 1 þ e ðB kt
B Þ
:
0 D0
where D and D0 are composed of different fixed parameters, which can be tracked
back following the same reference. B and B0 are the fixed parameters of the logistic
function that shapes the relational curve between the diffuse fraction and kt. The
probabilistic representations of this quantity have been addressed with different
probability density functions (PDFs) such as Weibull, Gaussian, double beta,
Boltzmann, single gamma, beta, bi-exponential, double Gaussian, triple Gaussian,
Weibull-logistic, etc. Single gamma [16] is adopted for this work.
The performance of the presented model is depicted in Fig. 2.3. In this figure, the
solar radiation received on the surface of the connected PV is represented in z-axis
(black striped area edged with red line) and varies according to y-axis that represents
different times during a single day. For demonstration purposes, three different times
(8:00 AM, 12:00 PM, and 1:00 PM) are considered to represent the probability
densities, which are aligned according to the axis of day hour. The x-axis, named
“power flow net,” represents the actual power import/export from or to the utility
grid. In each single hour, the probability density of solar irradiation is depicted
against the power flow net, showing the most likely occurrence of export/import
power from a PV during single hour in a day. In this step, only single hour is
considered for the uncertainty characterization.
The large number of factors considered in the evaluation of the system increases the
complexity of the evaluation task. The nonintrusive method is applicable for numer-
ical based evaluation taking into the account the recent advances in computing
power. It has been utilized in the recent studies that concern the problems of
18 H. Al-Saadi
power system such as the computation of the probabilistic power flows [17, 18]. The
nonintrusive methods govern the response variable to the uncertain input quantities
without the need to deal with the internal mathematical relation. It views the system
as a black box and allows us to analyze the stochastic performance cheaply through
the statistics such expectations, medians, quantiles, etc.
Let us consider the system to have uncertain functions represented by
n dimensional random vector, X, and parametrized by the fixed quantities
represented by d; then the function of the black box can be considered as g(X, d ),
and the response vector of variables of interest, represented by Y, can be expressed as
follows:
Y ¼ gðX, dÞ ð2:3Þ
Herein, the number of variables, n, in each vector is shared equally likely for PVs
and loads; the function g(X, d ) comprises internal computational relation that is not
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 19
checked by the author of this chapter. The internal relations are based on the used
software for computing power load flows.
The nonintrusive method simplifies the use of statistics such as sampling, quad-
rature treating, and regressions, by employing m stochastic runs of deterministic
evaluations.
Hypothesis testing is when the variables of interest are analyzed, and statistics-based
decisions are made. It should be noted that it is not possible to 100% accept or reject
the decision made out of a hypothesis. Instead, these decisions should be assessed
carefully in accordance with common traditional practices such as confidence
intervals, standard deviations, and variances. In the case of this work, the statistical
results of the sparse grid technique to the Monte Carlo technique is tested by relative
error and used for driving the statistical analysis such as the use of expectation in the
histogram and bar-stem diagram.
where Θnz denotes the index set of multi-index vectors, fℓi gni¼1 , such that Θnz ¼
n P o
ℓ 2 fℕgn : jℓ j ¼ di¼1 ℓ i ¼ z þ n &f∅ : z < 0g , and U ℓz ¼ ð1Þℓz n1 ℓz .
The extended formula of the tensor product is as follows:
Xm1 Xm2 XM n
1
Λl1 ⨂ ⨂Λln ¼ i1 ¼1 i2 ¼1
rn ¼1
w i1 ⨂ ⨂wn
in
g x1i1 , , xnin ∙ f x1i1 , , xnin : ð2:6Þ
Due to the nesting feature, the selection of Kronrod-Patterson rule is adopted for
univariate quadrature. This rule adds more nodes to the preceding level of accuracy,
li 1, thereby improving the exactness in the dimension evaluation.
The involvement of the risk assessment require clarification and definition for the
sake of approaching clear and directive framework. In this regard, the term “risk” has
been introduced by ISO/IEC [21] to be “The combination of probability of an event
and its consequence.” ISO/IEC also introduces “risk analysis” to be “systematic use
of information to identify and to estimate the risk.” Later, ISO replaced the risk’s
definition into “Effect of uncertainty on objectives” [22]. In Australia/New Zealand,
AS/NZS ISO 31000:2009 adopted the last definition [23]. In addition, the definition
of the term “risk” has been addressed, by Kaplan and Garrick [24], using triplet
questions: what are the possible events, what are the consequences of these events,
and what are the associated probabilities of an event? Afterward, these triplet
questions become the main contributor for defining the so-called quantitative risk
assessment (QRA), especially in engineering applications, such as in [25]. After
identification and description of possible hazardous phenomena, the relative fre-
quency of occurrence (occurrence probability) can be used to quantify the identified
risk [26].
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 21
The violations of the prespecified rated current in the line and standardized voltage
ranges are two problems identified previously in [27] and discussed in [28, 29] to
have a negative impact on the operational performance of distribution network when
large number of PVs is connected. These steady-state operational constraints can be
used to assess the degree of risk exposure to the network. The degree might be eased
or hardened depending on the technical standards adopted by the distribution power
utility.
During an hourly assessment of the operational performance (as discussed in
Fig. 2.3), the violations concerning bus voltages and ampacity of lines are consid-
ered for formulating the deterioration indices, signified as D.Index 1 and D.Index
2, respectively. Each risk can be quantified based on the number of violation and on
the extent of consequence according to [26]. Both can characterize the likelihood and
severity of risk components. Therefore, in each risk index, the number of violation
(s) occurring in the entire network, denoted by N , and the extended consequence of
the violation(s), denoted by D , are considered. The extended consequence, D, is
formulated here as an accumulative sum of the percentage of the overvoltage
subtracted from the maximum acceptable voltage bound, symbolized by Vmax. The
following steps counted for overvoltage violation, D.Index 1, are used to determine
the values of N and D during each stochastic pocess run:
First step Set N 0 ¼ 0 and D 0 ¼ 0.
Second step If |Vbus(i)| > Vmax, then N i ¼ N i1 þ 1 and
Third step D i ¼ D i1 þ jV bus ðViÞmax
jV max ðiÞ
ðiÞ , where i refers to the node number.
Fourth step Repeat step 2 until all buses are covered.
Fifth step Store the data so that N ) D:Index 1N and D ) D:Index 1D .
After performing m-times of evaluations, each risk index can be quantified by
computing the expectation of its likelihood of violation’s occurrence and severity of
its consequence. Following the black box in Fig. 2.2, functional relation between
joint PDFs denoted by fjoint(X, d ) and internal functions of the system denoted by g
(X, d) can be used to estimate the expectation of likelihood as follows:
Z
½D:Index 1N ¼ D:Index 1N ðgðX, dÞÞ ∙ f joint ðX, d Þ dX, ð2:7Þ
X
Z
½D:Index 1D ¼ D:Index 1D ðgðX, dÞÞ ∙ f joint ðX, dÞ dX, ð2:8Þ
X
The probability-consequence diagrams are common for risk visualization. Others are
risk bubble representations, curve approaches, uncertainty boxes, and strength-of-
evidence assessment [30]. However, a precise diagram has not been recommended
yet for the era of the current power grids. This leaves the debate open as the
informative graphical diagram for risk display may rely on the nature of the problem,
which can be ultimately decided by the risk analysts/reporters.
The expectations of the system performance when increasing multiple PV con-
nections are displayed in terms of two directions of risk index, i.e., likelihood and
severity. The likelihood takes the vertical axis, whereas the severity takes the
horizontal axis. The stepwise increases of PV connections are depicted as a third
axis using scatter points and distinguished by the size/color according to
jet-colored bar.
Using expected values, streamlining different regions, to extinguish the state of
risk from being negligible into unacceptable, is edged by a stripped red line. To
clarify further, the striped red line partitions diagram into two areas: not shaded area
where the PV connections expose no thread to system’s performance. The red
shaded area indicates that PV connection is expected to expose the system into a
thread of technical violation. In Fig. 2.4, three risk indices are drawn, and the values
of these indices when crossing among the specified areas are shown in the jet-colored
bar. Depending on the significance of each index, the limit of HC is decided
considering these three indications. Mathematically speaking, the limit of HC can
be formulated if the maximum bound of risk index is preassigned such as maximum
bound of likelihood, Lmax, and maximum bound of severity, Smax, for D.Index 1, as
follows:
ðD:Index 1N , D:Index 1D Þ : D:Index 1N Lmax
HCðD:Index 1Þ ¼ , ð2:9Þ
ε _ D:Index 1D Smax ε, ε 2 ℕ
where ε means the degree of allowance error that depends on the resolution of the
stepwise PV increase.
Similarly, other HC limits for other D.Index can be driven and concluded for final
HC determination.
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 23
Two case studies are presented in this chapter to show the effectiveness of this
approach. The work is implemented using MATLAB R2018b on a desktop com-
puter (Intel® i7-2600 4-core with 3.4GHz processor). In this section, the simulation
of superior of spare grid technique is tested against Monte Carlo technique. Then, the
approach of HC is implemented on two distribution feeders and networks named
11 bus distribution feeder and large distribution network. The operational voltage at
the main substation transformer for both cases is kept constant via On-Load Tap
Changer (OLTC) to be 1 pu all the time (the effect of OLTC on HC is out of the
scope of this chapter). The PDF of the clearness index, kt, is set to follow the time
characteristics of Adelaide city (35oS), Australia. The values of parameters of the PV
model introduced in (2.2) for three different times are shown in Table 2.2. The loads
are characterized by Gaussian PDFs using the 3-σ principle around the load maxi-
mum values. The nominal voltage is expected to be 1 pu, and voltage exceeding
1.1 pu is considered as an overvoltage violation, as specified by voltage standard
EN50160 [14]. An overcurrent is counted when the current flow at any line within
the targeted zone exceeds 150% of that line’s rated-current capacity.
In the comparison between Monte Carlo technique and sparse grid technique, two
quantities are tested: relative error and computational time. The reference result is
computed by the number of 40,000 samples using Monte Carlo technique. The aim
24 H. Al-Saadi
Fig. 2.5 Relative errors against the number of samples for fair comparison between Monte Carlo
technique and sparse grid technique
of using the relative error is to compare the modelling accuracy. The results are
shown in Fig. 2.5 displaying the values of up to 1000 samples. The relative errors
computed using sparse grid technique shows a steady tendency toward the optimal
values, i.e., without deviating arbitrarily. This is because the coverage of Monte
Carlo technique is completely random and not expected to be steadily improving,
unless a large number of samples are used, which is computationally inefficient.
Following the aforementioned analysis, the sample number used for the computation
of Eqs. (2.7) and (2.8) is 1000 that induces the relative error equal or less than the
value of 0.003.
Another feature of sparse grid technique is the time speed performance. In
Fig. 2.6, Monte Carlo technique seems to perform faster with the number of samples
less than 200. The last result is not beneficial due to the lack of stable relative errors.
Instead, sparse grid technique performs in the sample number of 1000 as five times
faster as the performance of Monte Carlo technique. Therefore, sparse grid technique
is employed only in the following simulation results. It is obvious that sparse grid
proves its superiority over Monte Carlo technique, explained by selecting the finest
number of samples through a unique quadrature rule mentioned earlier. It should be
noted that this feature is considerably important in the purpose of providing accurate-
fast hourly risk assessment that can be incorporated easily with other time-based
automation approaches in the modern power grids.
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 25
Fig. 2.7 Radial distribution network with nine points of customer connections feeding
36 PV-installed households
Table 2.3 Details of the maximum contractual loads (kW) in each house, point of customer
connection, and feeder
Level Middle feeder Sided feeder
Each house 3.3 3.3
Each point 9.9 19.8
Entire feeder 59.7 29.7
An actual distribution network with a 11/0.4 kV substation (see Fig. 2.7) is used to
implement the proposed risk-based approach. The technical details are shown in
Tables 2.3 and 2.4 (refer to [31] as the original source of this feeder). The rated
power of the substation transformer is 125 kVA. With an average of maximum
contractual load of 3.3 kW in each house and typical load profile shown in Fig. 2.8,
the performance of the network is estimated through the aforementioned risk indices.
The assessment is conducted through a 10-kW stepwise increase of PV connec-
tions at the point of common coupling which is, in this case, the same as the point of
customer connection. At each step, stochastic run is performed with the use of
26 H. Al-Saadi
2 4
8:00AM D.Index 1 D.Index 2
Likelihood & Severity
1.5 3
1 2
0.5 1
0 0
a) b)
6 6
10:00AM
D.Index 1 D.Index 2
Likelihood & Severity
5 5
4 4
3 3
2 2
1 1
c) 0 d) 0
8 8
12:00AM
D.Index 1 D.Index 2
Likelihood & Severity
6 6
4 4
2 2
0 0
e) f)
0 10 20 30 40 50 60 0 10 20 30 40 50 60
PV increase (kW) PV increase (kW)
Fig. 2.9 Three risk assessments for D.Index 1 and D.Index 2 at three different daytimes (8:00 AM,
10:00 AM, and 12:00 PM)
Table 2.5 Hosting capacity limits for 11 bus feeder according to likelihood only
Time Load percent (%) HC limit1 (kW) HC limit2 (kW)
8:00 AM 35 <50 <60
10:00 AM 60 <20 <20
12:00 PM 80 <20 <20
Table 2.6 Hosting capacity limits for 11 bus feeder according to likelihood and severity
Time Load percent (%) HC limit1 (kW) HC limit2 (kW)
8:00 AM 35 >60 >60
10:00 AM 60 <30 <30
12:00 PM 80 <20 <30
These two important indices can be used to identify two values of HC values: HC
limit1 and HC limit2. If only likelihood is considered which is the number of
technical violation(s), HC limits can be indicated as shown in Table 2.5. It is prudent
to first start by stating that the values of HC limits change according to the day hours
and loading percentages. In addition, more precise HC limit can be indicated if high
resolution of PV increasing steps is considered (10 kW is considered). In this table,
during the maximum PV production, which is at 12:00 PM, the HC limit1 and limit2
are 20 kW for both. These HC limits can be eased when considering the severity of
the violation as another factor to indicate the risk as shown in Table 2.6. It is
noticeable that some of HC limits are higher in Table 2.6 compared in Table 2.5,
indicating that not counting the severity may lead to overestimate the risk of PV
connections at some point. For instance, the violation of overvoltage could be
spotted in “1.5” buses in the entire feeder at 10:00 AM during 20 kW of PV
connection; see subfigure (d). But this violation does not pose a severe thread
which can be tolerated, subjected to the operators or policy makers, allowing more
PVs to be connected. According to Fig. 2.7, each house should be allowed to export
PV output power no more than 6.66 kW in the sided lines or 5 kW in the mid-line
during the hour centered at 12:00 PM to avoid overcurrent and overvoltage issues or
no more than 30/3 kW for the sided line and 30/4 kW for the mid-line if overcurrent
issue is tolerated but no voltage violation.
The approach is carried out on a real existing distribution network situated in South
Australia (see Fig. 2.10). The technical details of this network can be found in
[32]. The total number of 11/0.4 kV distribution transformers is 91. The network is
segmented into two zones: zone A and zone B. Zone A is assessed entirely as well as
two specified feeders: feeder 1 and feeder 2. The characteristics of these zones and
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 29
feeders are listed in Table 2.7. The cable codes mentioned in this table are 1 ¼ 61/
2.50-Flat (3 300 mm2) AL XLPE, 2 ¼ 6/1/3.75 ACSR, 3 ¼ 6/1/2.75 ACSR, and
4 ¼ 30/7/.102 ACSR. This network is modelled in OpenDSS, an advanced distri-
bution network modeller introduced by an American institute named EPRI
[33]. OpenDSS is co-simulated with MATLAB through Component Object Model
(COM) interface. The PV connections are increased in a firm stepwise and divided
by the maximum load of the transformer. The assessment for each zone or feeder is
done separately.
The use of bar and stem figures to display the risk seems to be directing toward
the HC indication, but not determination. It is because of the severity of each risk
index that cannot be judged due to the overlap with likelihood’s scale. For that
reason, the theoretical diagram shown in Fig. 2.3 is used to show the scale of both
likelihood and severity as appropriate as it can be to match the severity’s requirement
of determining HC. In this diagram, scaled axes of likelihood and severity are
displayed individually, allowing a fairer judgment for consideration. For example,
Fig. 2.10 Two zones and feeders under assessments within a large distribution network located in
South Australia
Fig. 2.11 HC
determination according to
the likelihood of violation
for two risk indices when
performing 0.5–4.5 PV
increases/load and estimated
during an hour duration
around 12:00 PM
Fig. 2.11 shows the expectations of D.Index 1 and D.Index 2, computed at 12:00 PM
for the entire zone A. For D.Index 1, the reddish shaded area in these figures
(as discussed previously in Sect. 2.4.2.1) is cropped for 1.5 severity, i.e., 150%
over specified rated ampacity, and 1 likelihood, i.e., one line expected to have
ampacity violations. For D.Index 2, the reddish shaded area in these figures
(as discussed previously in Sect. 2.4.2.1) is cropped for 0.1 severity, i.e., 1.1 pu
over nominal bus voltage, and 1 likelihood, i.e., one bus expected to have overvolt-
age violation. According to the jet-colored bar, HC limits can be determined as HC
limit1 to be 1.7 pu and HC limit2 to be 2 pu.
The selection of HC limits is meant to be conservative in Fig. 2.11. The reason for
conservativeness is that these violations are not expected to cause a severe overvolt-
age or severe ampacity exceedance. With these HC limits, the severity is just 5% of
rated ampacity (Fig. 2.11a) and 1.007 pu of over nominal voltage (Fig. 2.11b).
Following these low values of severities, the likelihood can be compromised with
higher values up to the thresholds which are 1.5 (150% over ampacity) and 0.1 (1.1
over nominal voltage). Therefore, HC limits can be further determined according to
severity; see Fig. 2.12. In this figure, the HC limit1 and HC limit2 are determined to
be 2.1 pu and 2.9 pu, respectively. In other words, the ampacity violation occurs first
and then the overvoltage violation in the scale of PV increase connections. That
would lead to a critical area that is shaded with yellow and requires more attention
from the side of policy makers, analysts, and researchers for finding ways for optimal
utilization; see Fig. 2.13.
2 Hosting Capacity Determination via Risk Analysis Approach Involving Likelihood. . . 31
Fig. 2.12 HC
determination according to
the severity of violation for
two risk indices when
performing 0.5–4.5 PV
increases/load and estimated
during an hour duration
around 12:00 PM
Critical
0.5 2.5 4.5
2.1 2.9
We assess other feeders and zones: feeder 1, feeder 2, and zone B. The outcomes
are presented in Table 2.8. These assessments reveal that the hosting capacity is
larger when the size of network or feeder is smaller in number of buses and lines.
However, these need a thorough investigation as well as other important key factors
such as the distance from the substation, the total length of the feeder, the adopted
standards for streamlining indices’ threshold, type of cables, the nominal voltage at
the substation, the percentage of loading, etc. In Table 2.8, the results that are shown
for only two HC limits indicate that the ampacity violation always occurs before the
overvoltage occurrence. The results of both risk indices according to Fig. 2.4 are
shown in Fig. 2.14 for zone B. We used the severity’s threshold to determine the
values of HC limits. It is advised that this particular early violation is to be
investigated further with more number of zones and feeders to find out whether
this is a coincidence or network related.
32 H. Al-Saadi
Fig. 2.14 HC determination for zone B when performing 1–13 PV increases/load and estimated
during an hour duration around 12:00 PM
2.6 Conclusion
residential networks being transferred to be active and feeding the nearby commer-
cial or industrial buildings or complexes around noon time. In addition, determina-
tion of HC can benefit a wide range of interested parties such technical analysts,
policy makers, distribution network operators, etc.
Future work is recommended for further development of this area. The fact that
the sparse grid technique proves to be superior over Monte Carlo technique is
explained by selecting the finest number of samples through a unique quadrature
rule mentioned earlier. The technique can be incorporated easily with other time-
based automation approaches in the modern power grids. The effect of OLTC on the
value of HC is worthy of investigation. Therefore, it is recommended to further
improve the PDFs of loads, rather than adopting Gaussian ones. Also, the critical
area is shaded with yellow and requires more attention from the side of policy
makers, analysts, and researchers for finding ways for optimal utilization.
References
1. Ismael, S., Aleem, S., Abdelaziz, A., & Zobaa, A. (2019). State-of-the-art of hosting capacity in
modern power systems with distributed generation. Renewable Energy, 130, 1002–1020.
2. Sarkar, M., Meegahapola, L., & Datta, M. (2018). Reactive power management in renewable
rich power grids: A review of grid-codes, renewable generators, support devices, control
strategies and optimization algorithms. IEEE Access, 6, 41458–41489.
3. Stetz, T., Marten, F., & Braun, M. (2013). Improved low voltage grid-integration of photovol-
taic systems in Germany. IEEE Transactions on Sustainable Energy, 4(2), 534–542.
4. Bollen, M. H., & Hassan, F. (2011). Integration of distributed generation in the power system.
Hoboken: Wiley.
5. McNab, S., Crownshaw, T., Lemon, S., Strahan, R., & Pandey, S. (2016). Determining network
congestion. Auckland: Electric Power Engineering Centre.
6. Dubey, A., & Santoso, S. (2017). On estimation and sensitivity analysis of distribution circuit’s
photovoltaic hosting capacity. IEEE Transactions on Power Systems, 324, 2779–2789.
7. Shayani, R. A., & de Oliveira, M. A. G. (2011). Photovoltaic generation penetration limits in
radial distribution systems. IEEE Transactions on Power Systems, 26(3), 1625–1631.
8. Rylander, M., & Smith, J. (2012). Stochastic analysis to determine feeder hosting capacity for
distributed solar PV. Palo Alto: EPRI technical report 1026640.
9. Breker, S., Claudi, A., & Sick, B. (2015). Capacity of low-voltage grids for distributed
generation: Classification by means of stochastic simulations. IEEE Transactions on Power
Systems, 30, 689–700.
10. Rigoni, V., Ochoa, L. F., Chicco, G., Navarro-Espinosa, A., & Gozel, T. (2016). Representative
residential LV feeders: A case study for the North West of England. EEE Transactions on
Power Systems, 31(1), 348–360.
11. Etherden, N., & Bollen, M. H. (2013). Dimensioning of energy storage for increased integration
of wind power. IEEE Transactions on Sustainable Energy, 4(3), 546–553.
12. Ochoa, L. F., Dent, C. J., & Harrison, G. (2010). Distribution network capacity assessment:
Variable DG and active networks. IEEE Transactions on Power Systems, 25(1), 87–95.
13. Etherden, N., Bollen, M., Ackeby, S., & Lennerhag, O. (2015, June). The transparent hosting-
capacity approach–overview, applications and developments. In International Conference and
Exhibition on Electricity Distribution.
14. CENELEC EN 50160. (2010). Voltage characteristics of electricity supplied by public distri-
bution systems.
34 H. Al-Saadi
15. Al-Saadi, H., Zivanovic, R., & Al-Sarawi, S. F. (2017). Probabilistic hosting capacity for active
distribution networks. IEEE Transactions on Industrial Informatics, 13(5), 1–1.
16. Hollands, K. G. T., & Huget, R. G. (1983). A probability density function for the clearness
index, with applications. Solar Energy, 30(3), 195–209.
17. Xiao, Q. (2015). Dimension reduction method for probabilistic power flow calculation. IET
Generation, Transmission & Distribution, 9(6), 540–549.
18. Ni, F., Nguyen, P. H., & Cobben, J. F. (2017). Basis-adaptive sparse polynomial chaos
expansion for probabilistic power flow. IEEE Transactions on Power Systems, 32(1), 694–704.
19. Smolyak, S. A. (1963). Quadrature and interpolation formulas for tensor products of certain
classes of functions. Soviet Mathematics Doklady, 4, 240–243.
20. Heiss, F., & Winschel, V. (2008). Likelihood approximation by numerical integration on sparse
grids. Journal of Econometrics, 144(1), 62–80.
21. British Standards Institution. (2002). Risk management - vocabulary - guidelines for use in
standards: ISO/IEC guide 73(2002). ISBN 580401782.
22. International Organization for Standardization. (2009). Risk management vocabulary: ISO
guide 73(2009).
23. Standards Australia/Standards New Zealand Standards Committee. (2009, November) AS/NZS
ISO 31000:2009: Risk management-principles and guidelines.
24. Kaplan, S. (1997). The words of risk analysis. Risk Analysis, 17(4), 407–417.
25. Nowak, J., Sarkani, S., & Mazzuchi, T. (2015). Risk assessment for a national renewable energy
target part I: Developing the model. IEEE Systems Journal, 9(3), 1045–1056.
26. Aven, T., Baraldi, P., Flage, R., & Zio, E. (2013). Uncertainty in risk assessment: The
representation and treatment of uncertainties by probabilistic and non-probabilistic methods
(pp. 1–186). Hoboken: Wiley.
27. Schwaegerl, C., Bollen, M. H. J., Karoui, K., & Yagmur, A. (2005, June). Voltage control in
distribution systems as a limitation of the hosting capacity for distributed energy resources. In
IET CIRED Conference. Turin.
28. Navarro-Espinosa, A., & Ochoa, L. F. (2016). Probabilistic impact assessment of low carbon
technologies in LV distribution systems. IEEE Transactions on Power Systems, 31(3),
2192–2203.
29. Mulenga, E., Bollen, M., & Etherden, N. (2020). A review of hosting capacity quantification
methods for photovoltaics in low-voltage distribution grids. International Journal of Electrical
Power & Energy Systems, 115, 105445.
30. Goerlandt, F., & Reniers, G. (2016). On the assessment of uncertainty in risk diagrams. Safety
Science, 84, 67–77.
31. Conti, S., & Raiti, S. (2007). Probabilistic load flow using Monte Carlo techniques for
distribution networks with photovoltaic generators. Solar Energy, 81(12), 1473–1481.
32. Al-Saadi, H., Zivanovic, R., & Al-Sarawi, S. (2017). Probabilistic hosting capacity for active
distribution networks. IEEE Transactions on Industrial Informatics, 13(5), 2519–2532.
33. Dugan, R. C., & McDermott, T. E. (2011, July). An open source platform for collaborating on
smart grid research. In 2011 PESGM IEEE. Detroit.
Chapter 3
A Methodology Proposal for Calculation
of Hosting Capacity Including Different
Power Quality Phenomena
Z m ðhÞ
E m ðhÞ ¼ E ð hÞ ð3:1Þ
Z m ð hÞ þ Z c ð hÞ c
Zm(h)
Im(h)
Fig. 3.1 Power system representing the background distortion model [1]
3 A Methodology Proposal for Calculation of Hosting Capacity Including Different. . . 37
Z c ð hÞ
E c ðhÞ ¼ E ð hÞ ð3:2Þ
Z m ð hÞ þ Z c ð hÞ m
It is important to highlight that for these models, both sides must be analyzed as
harmonic contribution, which means that the contribution of utility and the contri-
bution from the prosumer will be implemented by using a superposition method in
order to determine the background distortion at the PCC [4].
In order to understand the equivalent voltage at the PCC, the IEC voltage phasor
method cited in [4] allows us to find the total voltage harmonic distortion at the PCC
considering the background voltage and the customer installation. In this case, the
voltage harmonic distortion, before any nonlinear load on the costumer’s side is
connected, will be defined as V background. Moreover, the voltage phasor V emission is the
voltage vector that estimates the voltage drop on the costumer’s side. Thus, the
voltage vector at the PCC V PCC , which consists of the voltage background and the
customer’s harmonic voltage emission, is shown in Fig. 3.3.
Fig. 3.5 THD (Total Harmonic Distortion) of voltage at 230 kV during a Brazilian match on the
World Cup 2010 [6]
Still in accordance to [6], Fig. 3.5 in the 5th harmonic voltage at 230 kV in the
Northeast of Brazil during the Brazil versus Portugal match was measured.
According to [6], a typical 5th harmonic voltage on a typical day is under 0.5%.
The same source concludes that the impact of TV viewing must be considered on
normal days, as well as during highly popular events, as a way to analyze the impact
of the background distortion in the system. It is clear from these results that a
significant increase of harmonic pollution can take place, especially in the 5th
harmonic order. Those results show that the harmonic background can be severe
due to nonlinear loads connected at the same time.
In this section, the definition and requirements for the nonlinear contributions and
their background distortion were determined. It is well-known that harmonic distor-
tion at the PCC is caused by the harmonic sources on both sides of the system: the
utility, as well as the customer’s facilities. The theoretical aspects of the background
distortion have been reviewed in order to select the correct model. It is important to
highlight that the background distortion model has an important contribution to the
hosting capacity value regarding distribution systems.
Both concepts and applications of hosting capacity have been studied by researchers
in order to elaborate a new significant strategy to plan and improve network systems,
especially for new generation sources. For distributed energy resources, the hosting
capacity approach has been introduced as a communication tool between stake-
holders concerning the connection of DG to the electric grid [13]. In this context,
40 T. E. Castelo de Oliveira and P. F. Ribeiro
Fig. 3.6 Range curves of a generic performance index versus the amount of distributed generation
[13]
3 A Methodology Proposal for Calculation of Hosting Capacity Including Different. . . 41
Nowadays, more and more power electronics are being integrated into our electrical
system, especially at the distribution system level [9–13]. As is well-known, power
electronics are likely to inject harmonic currents and, consequently, cause a certain
level of disturbance through their background distortion. In general, the larger the
device is, the greater the impact will be, either positively or negatively [10]. In other
words, any class or type of power electronic device connected to the network will be
injecting harmonic currents, and consequently, voltage disturbance will occur as
well as voltage rise as a result of the root mean square (RMS) magnitude dependence
on voltage waveform. For example, a power electronic device connected at the PCC
in the system can be modelled as a set of controlled current sources, in parallel,
integrated with its respective impedance for each harmonic order, where each
individual source will be injecting a certain amount of harmonic current, in magni-
tude and phase, and as a consequence, a voltage RMS will appear. It is important to
highlight that the power electronic device will be inserting a harmonic current at the
PCC until order n, according to Fig. 3.7.
...
GND1 GND1 GND1 GND1
42 T. E. Castelo de Oliveira and P. F. Ribeiro
Vo PL+jQL
Vg
Load
The sum of these injected currents via the described model will impact the RMS
voltage at the PCC. These values are the product of the harmonic current magnitude
and the impedance linked to each controlled harmonic current source. Since the
RMS voltage is the square root of the sum of squares of each harmonic voltage, it can
be expressed using the following equation (3.3):
vffiffiffiffiffiffiffiffiffiffiffiffiffi
u 51
uX
V RMS t V 2h ð3:3Þ
h¼1
To consider harmonic injection by the source and distortion at the PCC, our
model needs to be rewritten based on typical parameters of an LV feeder model [10–
13]. Details of the feeder model used and corresponding voltage rise formulae are
illustrated in Fig. 3.8.
Based on the model shown in Fig. 3.8 for which the Thevenin equivalent
impedance at the PCC is given by Rf + jXf. The complex power injected by a general
DG model system and the complex power absorbed by the local load are defined as
Sg ¼ Pg + jQg and SL ¼ PL + jQL, respectively. The grid injected active and reactive
powers at the PCC are defined as P ¼ Pg PL and Q ¼ Qg QL, respectively
[10]. The fundamental frequency voltage difference between the PCC and the main
supply can be expressed by (3.4).
P þ jQ
V_ G V_ o ¼ ðRf þ jX f Þ ð3:4Þ
V_ G ejθ
The voltage at the PCC is denoted by V_ G :ejθ , and the voltage at the source is V0.
j0
e .
In order to simplify our assumptions, let us assume that cosθ 1 and sinθ 0,
and (3.4) can be approximated by (3.5) for the fundamental frequency h ¼ 1.
PRf þ QX f
ΔV 1 ð3:5Þ
Vg
Then, Eq. (3.5) can be simplified for other harmonic orders h 6¼ 1 as in the
following:
3 A Methodology Proposal for Calculation of Hosting Capacity Including Different. . . 43
Ph Rf þ Qh hX f
ΔV h ð3:6Þ
V hg
To assume the maximum voltage rise at the PCC, the load consumption (L) has to
be as minimum as possible, and the generation (G) has to be at its maximum value.
Thus, one can use (3.6) to compute an upper required for voltage rise assuming that
the power factor (PF) of net load can be controlled by the DG at the PCC system,
which means tan(φ) ¼ Q1/P1. Based on these assumptions, Eq. (3.6) can be rewritten
through (3.7) and (3.8).
ΔV 1 :V 1g P1 Rf þ Q1 X f ð3:7Þ
Going back to the assumption presented in (3.3), the active power generated in
relation to the RMS voltage can be rewritten as (3.10), using (3.3):
qffiffiffiffiffiffiffiffiffiffi
P 2ffi2
Vh P
V 2 V 2h
Pg ¼ RMS ¼ ¼ ð3:10Þ
Rf Rf Rf
In (3.9), the maximum power that can be injected by a DG was calculated, thus
defining the hosting capacity at the PCC as the sum of the maximum power possible,
based on (3.10), to be injected for each frequency as per (3.11) [10]:
X
h X51 V max ,h V max ,h V h
1 g g o
Pmax
g ¼ max ,h
PPV ¼ ð3:11Þ
h¼1
R h¼1 1 þ tan ðφÞ hX f
Rf
where:
V gmax ,h is the allowed voltage magnitude for harmonic order h
V ho is the actual voltage magnitude of harmonic order h at the PCC
φ is the arcos(PF), where the PF is the power factor at the PCC
44 T. E. Castelo de Oliveira and P. F. Ribeiro
If we agree to use the IEC and notice that, for V 1 kV, standards impose that
individual voltage distortion be below 2%, for instance, for h 2 and below 110%
for the fundamental frequency, h ¼ 1, then, the hosting capacity is given in (3.12):
1:10 1:10 V 1o Xh
0:02 0:02 V ho
Pmax
g ¼ þ ð3:12Þ
Rf 1 þ tan ðφÞ XRff h¼2 Rf 1 þ tan ðφÞ R
hX f
f
Several factors may affect the performance of calculating the hosting capacity in an
electrical system, such as the harmonic voltage, power factor, minimum load of the
distribution feeder, and voltage regulation equipment, among others [10]. We will
discuss the harmonic distortion impact on the presented Eq. (3.12).
Using (3.12), it is possible to find the value of the hosting capacity considering the
voltage harmonics from the fundamental order up to the 50th harmonic order, as well
as performance index limits, which depend on local standards [10], for instance,
IEEE 519 or IEEE 1547 or even IEC. Moreover, the injection of harmonic current
into the system depends on the power electronic equipment connected at the PCC. It
is important to highlight that the power electronics can be controlled by a prosumer,
for instance, due to the fact that the equipment is manufactured according to
international standards. According to [10], the background voltage at the PCC
cannot be controlled as it is almost impossible to know what is causing the harmonic
voltage distortion. Therefore, one extreme scenario may be imposed to analyze how
the harmonic injection can impact the hosting capacity.
Consider a scenario where the voltage for the fundamental order is set at 1.00 p.u.
and the other harmonic voltages for h 2, up to 15, are set at 0.02 p.u. For this
scenario, a high distortion is considered as it is dependent on the background and all
the electronic devices connected at the PCC. The other scenarios consider a non-high
voltage distortion for h 2 up to 15, and the voltage distortion will set in at 0.01 p.u.
In reality, a decrease in harmonic voltage would increase the hosting capacity as
shown in the equation below:
P
Pgmax 0 1:1ð1:1 1Þ þ ð0:02 0:01Þ0:02 0:113
¼ P ¼ ¼ 1:027
Pmax
g 1:1ð1:1 1Þ þ ð0:02 0:02Þ0:02 0:11
References
1. Ghijselen, J. A., Ryckaert, W. R., & Melkebeek, J. A. (2002). Required load behaviour for
power quality improvement. In IEEE Power Engineering Society summer meeting (Vol. 3, pp.
998–1003).
2. Ronnberg, S. K., Members Zavoda, F., Convenor CA Corresponding Members, et al. (2018,
March). Power quality and EMC issues with future electricity networks.
3. Ryckaert, W. R. A., Ghijselen, J. A. L., & Melkebeek, J. A. A. (2003). Harmonic mitigation
potential of shunt harmonic impedances and the influence of background distortion. In IEEE
Bologna Power Tech conference proceedings (Vol. 2, pp. 663–670).
4. Papic, I., et al. (2018). A benchmark test system to evaluate methods of harmonic contribution
determination. IEEE Transactions on Power Delivery, 8977(1), 1–8.
5. Browne, N., Perera, S., & Ribeiro, P. F. (2007). Harmonie levels and television events. In 2007
IEEE Power Engineering Society general meeting, PES (pp. 7–12).
6. Ribeiro, P. F., et al. (2011). Harmonic distortion during the 2010 FIFA World Cup. In 2011
IEEE Power and Energy Society general meeting (pp. 1–8).
7. Macedo, J. R., Martins, A. G., Carneiro, J. R. V., & Siqueira, M. J. V. (2006). The impact of
FIFA World Cup 2006 on power quality in the electric distribution systems. In 2007 9th
international conference on Electrical Power Quality and Utilisation.
8. Leitão, J. J. A. L., Fonseca, L. C. A., Lira, M. M. S., Soares, L. R., & Ribeiro, P. F. (2007).
Harmonic distortion on a transmission system during games of the Brazilian National Team in
the 2006 world cup. In 2007 IEEE Power Engineering Society general meeting, PES (pp. 1–2).
9. Testa, A., & Langella, R. (2007). Harmonic pollution in Italian distribution networks in
coincidence with important sport events. In 2007 IEEE Power Engineering Society general
meeting, PES (pp. 1–7).
10. de Oliveira, T. E. C., Carvalho, P. M. S., Ribeiro, P. F., & Bonatto, B. D. (2018). PV hosting
capacity dependence on harmonic voltage distortion in low-voltage grids: Model validation
with experimental data. Energies, 11(2), 1–13.
11. F. E. Engineering. 2015, February. PQ-SmartGrid laboratory (pp. 1–16).
12. Bollen, M. H. J., & Rönnberg, S. K. (2017). Hosting capacity of the power grid for renewable
electricity production and new large consumption equipment. Energies, 10(9), 1325.
13. Dubey, A., & Santoso, S. (2017). On estimation and sensitivity analysis of distribution circuit’s
photovoltaic hosting capacity. IEEE Transactions on Power Apparatus and Systems, 32(4),
2779–2789.
Chapter 4
Effect of Harmonics due to Distributed
Energy Resources on Hosting Capacity
of Microgrid: A Hardware in Loop-Based
Assessment
Abbreviations
CSP Concentrated solar power
DERs Distributed energy resources
DNI Direct normal irradiance (W/m2)
DNIc Clear-sky direct normal irradiance (W/m2)
EHC Enhanced hosting capacity
FPGA Field-programmable gate array
HC Hosting capacity
HIL Hardware in loop
HPR Harmonics penetration ratio
PCC Point of common coupling
SCADA Supervisory control and data acquisition
SPV Solar photovoltaic
THD Total harmonics distortion
WECS Wind energy conversion system
Nomenclature
Ai Area of the wind turbine blades
A Global irradiance on 2-axis tracking plane (W/m2)
Ac Global clear-sky irradiance on 2-axis tracking plane (W/m2)
Ad Diffuse irradiance on 2-axis tracking plane (W/m2)
4.1 Introduction
The transition of conventional power system toward smart grid has led to the
penetration of nonconventional energy sources in the form of distributed energy
resources (DERs) [1]. This is due to changes in type of energy production, that is,
shift from large-scale generation unit to small units integrated in distribution system;
type of customers, that is, the movement toward sustainable energy system is driving
a change to efficient energy system; decrement of fossil fuel; and environmental
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 49
aspects [2]. However, increased integration of DERs may hinder the grid perfor-
mance, leading to violation of voltage, thermal limits, added injection of harmonics,
and mal-operation of protection [3]. Thus, integration of these DERs should be done
as per the network hosting capacity (HC). The HC is defined as the quantity of added
generation or utilization that can be integrated to the grid without endangering the
operating limits or network reliability [4].
To maintain the operational parameters within the predefined standards and to
avoid mal-operation of the equipment, node voltage violation, ampacity of the
conductors, and uncertainty should be taken as constraints while calculating the
HC [4]. Parameter such as node voltage play a significant role in the hosting capacity
assessment due to the fact that voltage violation limit is reached far earlier than the
current limit. Upon violating the voltage limits, the distribution side may face the
uneven voltage profile which may cause local circulating currents. These circulating
currents may contribute to ohmic losses and may lead to false tripping of protecting
equipment. Ampacity should be the subsequent deciding factor, as extreme deviation
in the ampacity may result in thermal breakdown of the lines, equipment, and
conductor insulation [5]. As most of the DERs constitute solar and wind, a high
uncertainty is associated with both these sources. Any impulsive change in these
sources may result in sustainable disturbance in the grid. So to utilize these energy
resources to their optimal level, power electronic converters are used. Semiconduc-
tor switches such as IGBTs and MOSFETs which are used in these converters are
dependent on the switching frequency requirement and the uncertainty in the input
sources. The switching of these converters introduces harmonics in the system
[6, 7]. In addition to these, different nonlinear loads such as arc furnaces, saturated
magnetic core of various machines, and induction motors also contribute to the
harmonics of the utility [8].
The presence of harmonics in the distribution system causes various problems
such as degradation in voltage profile, heating of the wires and machines, and low
injection of DERs. These issues in the utility collectively hampers the HC of the
network. It also creates a problem in optimizing the size of DERs [9–12]. Although
these energy sources have negligible harmful environmental effect, they create
harmonics to a considerable extent [13, 14]. With the advancement of the maximum
power point tracking, for solar photovoltaic as well as wind power plants, the system
can get maximum power irrespective of varying sources [15, 16]. Due to the
arrangement of maximum power, the firing in the converter changes dynamically,
which creates additional distortion in the grid current and in voltage. The renewable
sources are location and time dependent. If the plant is in a geographical location
where the uncertainties of the resources are frequent, then the utility will be exposed
to more such disturbances. Due to these disturbances, the operational difficulty such
as variation in HC is a must, which is undesirable [17]. To address these operational
problems, accurate analysis of these resources is necessary to achieve a reliable and
effective operation of the utility without any voltage or ampacity limit violation.
In this chapter, effect of harmonics on HC of microgrid is tested for different
DERs. The DERs, particularly solar and wind, are uncertain in nature. Based on the
uncertainty of the different DERs, total harmonic distortion (THD) with voltage and
50 S. K. Sahu et al.
ampacity limits as the major constraint is measured. The effect of THD on the HC of
the network is monitored, which is a significant aspect for planning the size and
number of DERs in the microgrid. The entire analysis is done in MATLAB and
validated in Typhoon HIL (a real-time simulator). Typhoon HIL simulation is a
platform that includes the complexities of the plant and tests them in real time with a
controlled environment. All power system dynamics can be added in the system, and
by simulating them, the actual microgrid response can be predicted. This technique
is a realistic technique for testing, system improvement, and predicting its dynamic
behavior [13, 14]. The analysis is tested on 5-bus and IEEE 33-bus distribution
network with different DERs at a specific geographical location of India obtained
through global information system (GIS), to efficate the results. The results are also
verified at variable load conditions. Further on analyzing all the constraints, the
possible HC is proposed.
The arrangement of this chapter is as follows. The concept of HC is highlighted in
Sect. 4.2. Section 4.3 details the microgrid and its components. Section 4.4 illustrates
the effect of DER on hosting capacity. Section 4.5 details the tools used. Section 4.6
depicts the case study and results. Section 4.7 concludes the findings and highlights
the possible future work.
In the twenty-first century, energy consumption per capita is considered to be the one
of the most important indicators to know how much an individual is coping up with
the current time. In order to align with the demand of the situation, the energy
consumption pattern, quantity, and quality are drastically changing. As the demand
is increasing in never before rate, the generation should also match the increasing
demand. Due to continuous evolution in the field of science and technology,
renewable sources like solar photovoltaic, small hydro, and wind energy conversion
system are no more exorbitant. In this context, solar photovoltaic and wind energy
conversion system are most commonly used. Use of these energy sources is contin-
uously changing the dynamics of the power system, from unidirectional power flow
to bidirectional power flow. Thus, to regulate and to guide this integration process,
hosting capacity term is coined in power system. Hosting capacity can be defined as
the production or consumption of new energy considering the indicator/indicators to
lie within its allowable limits. Some of the indicators for hosting capacity are
voltage, current, power, and frequency. As frequency and power remain same as
compared to voltage and current, the indicators, that is, voltage and current are
widely used in the distribution system HC analysis. This hosting capacity can be
divided into two types:
(i) Dynamic hosting capacity
(ii) Static hosting capacity
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 51
Utility
Thevenin
Impedance
Gen
Loads
Microgrids are the modern local grids comprising of renewable and non-renewable
DERs that can operate in two different modes, namely, grid-connected mode and
standalone mode. Figure 4.2 depicts a typical microgrid structure with different
types of DERs and load integration. Some of the DERs are discussed in the
following subsection.
52 S. K. Sahu et al.
Utility
Gen
The wind energy conversion system (WECS) consists of various parts, such as rotor,
generator, gear transmission, tower, boost converter, inverter, and wind direction
sensor [18]. These are some of the most essential parts which are prerequisites for
further analysis. In addition to it, some of the operational parameters such as cut-in
wind speed, cut-out wind speed, cut brake wind speed, and Betz limit are also very
essential for the smooth operation of the WECS. Operational parameters such as
cut-in wind speed and cut-out wind speed play an important role in wind power
generation. To convert the wind current into useful electric power, there is a limit;
that is, all the energy content in the air particles cannot be converted into useful
power, which is detailed by Betz limit. The cut-in speed and cut-out speed is
normally 12.6 kmph and 90 kmph, respectively [19]. Figure 4.3 shows the cut-in
and cut-out speed graphically. The velocity of wind in between these values is
responsible for energy generation.
The wind turbine generates power by rotating the shaft with the help of the flow
of the wind through its blades. The wind which passes through the wind turbine
slows down as it gives up its energy. The kinetic energy that can be converted by the
wind turbine rotor is 59.3%, which is obtained when the residual wind, that is, the
outgoing wind after hitting the turbine blades from the wind turbine, is one-third of
the free wind. This is known as the Betz limit [20]. If the blades have ideal efficiency,
then the residual wind will result in zero value, that is, no wind condition leading to
the stop of operation of wind farms. The power content in the wind (PWind) is
depicted in Eq. (4.1).
1
PWind ¼ ρ Ai v3inf ð4:1Þ
2
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 53
Power(kW)
Rated Power
The wind power plant produces energy when the wind speed is in between cut-in
speed and cut-out speed. Equation (4.2) depicts the amount of power generated (PG, WT)
by the wind turbine.
8 9
>
<0 80 < v < vcut‐in >
=
PG,WT ¼ Jvk þ LPRated, 8vCut‐in < v < vrated ð4:2Þ
>
: >
;
PRated 8vCut‐in < v < vCutout
where
!
PRated
J¼ 3
V rated V 3 Cut‐in
!
PRated
L¼ 3
V Cut‐out V 3 Cut‐in
Solar energy is one of the widely harvested renewable sources of energy due to easy
application and installation. Solar energy can be used in two methods, namely,
concentrated solar power (CSP) and solar photovoltaic (SPV) power. In CSP, solar
thermal energy is used to make steam and run the turbine, but in SPV, plants photo
panels are used for conversion of solar energy to electric energy. Because of its
compactness, easy availability in the market, and low maintenance cost, SPV is
widely used in the modern power system. The block diagram of a typical solar
photovoltaic plant is shown in Fig. 4.4.
54 S. K. Sahu et al.
Utility Grid
IPV
RSe
RS RRSC V
IPhoto
IRSC
In the SPV plants, photovoltaic cells are the fundamental components of the
photovoltaic array. These are semiconductor devices which convert solar energy to
electric energy. This phenomenon is called the photovoltaic effect. The electric
power generated at a particular plant is a function of the geographical location,
irradiation, and optimal tilt angle. Each of these values has significant impact on
power generation. For a particular power rating of the solar photovoltaic plant, PV
arrays are connected in series to match the voltage rating and in parallel to provide
required load current. Photovoltaic cell can be modelled by a current source named
photovoltaic current (IPVC). The equivalent circuit diagram is as in Fig. 4.5, where
RSer is the series resistance and RShunt is the shunt resistor. Generally, the value of
these resistances are very less, so they are neglected to have a less complex analysis
of the circuit [21].
By solving the circuit shown in Fig. 4.5 with the above stated assumptions, the
PV array output current can be written as in Eq. (4.3) below. This current is the
resultant of NS number of parallel modules.
qV o
I PV ¼ N S I SN N S I RES exp 1 ð4:3Þ
kTAn N Se
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 55
The reverse saturation current (IRCS) in the cell is dependent on the temperature,
which is as in Eq. (4.4).
3
T Egap 1 1
I RSC ¼ I RES,I exp q ð4:4Þ
TI kAn TI T
R
I Photo ¼ Cq ðT T I Þ þ I C,SC ð4:5Þ
100
P ¼ V I PV
qV o
P ¼ N S I SN V N S I RES V exp 1 ð4:6Þ
kTAn N Se
Storage of the surplus power during the off-peak hour is one of the major aspects
of renewable integration. In this context, battery storage is very vital. Presently, the
invention of high-density Li-ion batteries is cost-effective, and these batteries have
very high energy density and can store an appreciable amount of energy. They are
used for avoiding overloading and peak saving. These batteries can be connected to
the microgrid to improve HC as well [22, 23]. Supercapacitors are also an alternative
to energy storage mainly due to their high energy density and quick charging-
discharging capacity [22].
Any fuel that is derived from the biomass is called biofuel. It is considered to be
renewable energy unlike fossil fuel as the source of biofuel is readily available like
animal waste and crop waste. Other sources such as animal fat, rapeseed oil, and
soybean oil can be used as a biofuel for generators. The net heating value can vary
between 42 MJ/kg for liquid biofuel and 39.7 MJ/kg for rapeseed oil. An example
for such system is Wartsila 9L20, which gives 1539 kW of electrical power output at
11 kV [24].
56 S. K. Sahu et al.
IG
CI ¼ ð4:7Þ
I Ext
On the other hand, the randomness for wind energy conversion system can be
determined by considering the positional and time-dependent wind flow. As the
wind pattern is totally random, various distribution functions can be implemented by
considering the best fitting of the curve. From [26], the Weibull distribution function
is taken for representing the wind speed variation as in Eq. (4.8).
k 1 k1
k 1 v 1 kv2
F Wind ðvÞjk1 ,k2 ¼ e ð4:8Þ
k2 k2
The vast use of the nonlinear types of equipment and renewable energy sources
distorts the sinusoidal current and voltage signals, which forms different types of
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 57
harmonics in the power system networks. These distorted voltage and current signal
causes various problems such as heating of equipment, false tripping of the protec-
tion equipment, and power quality issues. The decrement in the power quality and
heat formation should be kept minimal by restricting the harmonics to a safer limit as
per the standard, but the decrement in the power quality may lead to the restriction in
addition of renewables to the utility. Other than harmonics, voltage profile and
ampacity of the conductors are the key factors in the hosting capacity evaluation.
Dueto the addition of DERs, the voltage profile at point of common coupling (PCC)
and the ampacity will change. The recommended allowable deviation for different
voltage level is detailed in IEEE 519 [27, 28]. This change in the voltage or change
in the load current will cause various deformations in the network parameter profile
which may lead to decrease in the hosting capacity [29]. Keeping all the above
discussed factors in mind, the optimal sizing of the DERs can be decided for the
improvement of the HC. The optimal size of the DERs can be decided by taking
application and location into account. To optimize the sources, various methods,
parameters, and constraints can be considered [30, 31]. Various optimizing methods
such as particle swarm optimization (PSO), ant colony optimization (ACO), and
genetic algorithm (GA) can be implemented by considering the different application
of the network [32]. Mainly, the objective function and the constraints are required to
decide which optimization is applicable to the problem, and accordingly, size of
DER and types of DER along with its location can be decided. Network parameters
such as X/R ratio play an important role in optimization as the convergence of the
solution is greatly affected by this ratio [33]. Network reinforcement is one of the
methods to enhance the HC in the distribution system. The injection of the DER is
limited by the network constraints like voltage rise, thermal limit, network ampacity,
and the protection of the equipment. Along with the technical constraints, utility
rules may also restrict the DER investors to add more renewable to the utility grid.
4.5.1 MATLAB
Typhoon HIL is a scientific research software and hardware package for the real-
time testing of the power system. With its advanced features, it can simulate the
accurate condition with industry standard. Typhoon HIL simulators can be config-
ured in various modes where the user can define time according to the dynamics of
the plant. If the plant dynamics is very fast, Typhoon HIL ultra-high-fidelity
simulators can simulate the plant as low as a 1 μs time step. If the plant requires
high time steps for a specific application, then the time step can even be scaled to
10 μs. For a very large plant with multiple controllers and machines, the emulator
can simulate the plant with different time steps. These time steps are independent of
each other. Typhoon HIL real-time emulators also allow the user to optimize the
simulation so that the cores of the processor can be optimally used. Typhoon HIL is
an adaptive HIL manufacturer, which is contributing to the real-time emulator by
consistently improving the designing, testing, and verification environment. The
fast-growing microgrid library for the HIL application is revolutionary as it will help
the designers to design the microgrid with minimal approximation. Minimizing the
approximation will result in obtaining better real-time response of the system. The
FPGA-based HIL platform ensures quality performance of the emulator, and this
emulator can be used in various fields, such as:
(i) Power system and its application
(ii) Microgrid and its applications
(iii) Various FACTS devices and their application to the power network
(iv) Rapid control prototyping
To perform Typhoon HIL simulation, various series of emulators are present such
as 402, 600, 602, 602+, 603, and 604. The selection of emulator depends on various
parameters, such as number of machines, number of converters, and number of pulse
width modulation (PWM) channel requirements. Typhoon HIL emulator gives
precise control over power flow, voltage, and current at different magnitude and
phase. The various realistic grid phenomena such as sag and swell in voltage, spikes,
and rate of change of frequency can be simulated easily. Typhoon HIL can be easily
connected to the system by USB cable. It can also be connected by an ethernet port.
The interface of Typhoon HIL402 emulator is shown in Figs. 4.6, 4.7, and 4.8.
Typhoon HIL emulator works on FPGA platform with core count starting from
four in numbers. In case of higher number of core requirement, the researcher may
opt for a higher-end model, which can operate in standalone mode as well as parallel.
Each switching element in the library is assigned a weight, and each core can run in a
predefined value of weight. While simulating a model in HIL device, all the weights
and core counts are considered. The following steps are elaborated with figures for
understanding the operation through Typhoon HIL. The following steps may be
executed for connecting the HIL with the plant.
60 S. K. Sahu et al.
4.5.2.1 Step 1
Consider the core count and weight that will be used in the model. Choose a suitable
emulator, followed by a suitable configuration from the device table constituting of a
number of machine solvers as in Figs. 4.9 and 4.10, respectively.
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 61
Fig. 4.12 Selection of device (V-HIL or actual HIL) for model execution
4.5.2.2 Step 2
Create the model which is to be tested in HIL schematic editor. If multiple weights
are present in the model, then use core separator from Typhoon HIL library to assign
different weights to different physical cores of the emulator according to the
configuration table as shown in Fig. 4.10 and discussed in Step 1.
4.5.2.3 Step 3
Compile the model and run the model by clicking the button as shown in Fig. 4.11. If
there will be an error, it will appear at the bottom section, and required action needs
to be taken to fix it. After successful compilation, a command titled “compilation
finished successfully” will appear on monitor.
4.5.2.4 Step 4
Load the model into the real FPGA-based Typhoon HIL device connected to the
system by clicking on it as in Fig. 4.12.
62 S. K. Sahu et al.
4.5.2.5 Step 5
Design the supervisory control and data acquisition (SCADA) interface according to
the output requirements, that is, voltage waveform, current waveform, frequency
response, power measurement, THD calculation, etc. Data recorder to record data for
further analysis in various formats such as “.csv and .mat” is also available. For
changing the input during the execution time, there are action widgets such as knob
and slider in the SCADA panel as shown in Fig. 4.13.
4.5.2.6 Step 6
As shown in Fig. 4.13, on the right side, there are model settings. In model setting,
input all the data in the SCADA input section.
4.5.2.7 Step 7
Make sure the required contactors in the model settings are in CLOSE mode.
For analysis purposes, two test systems are considered: 5-bus distribution system
and IEEE 33-bus distribution system with variable solar irradiation and variable
wind speed, and biofuel-based power are used. The networks are also simulated with
variable load conditions. For testing, two networks, MATLAB and Typhoon
HIL402, are used.
Table 4.1 Solar irradiation data
Location: Arunachal Pradesh February
Inclination of plane: 43 (28 , 97 )
Time 06:52 07:07 07:22 07:37 07:52 08:07 08:22 08:37 08:52 09:07 09:22 09:37
Irradiation 0 0 0 0 0 0 414 429 442 452 462 469
Time 09:52 10:07 10:22 10:37 10:52 11:07 11:22 11:37 11:52 12:07 12:22 12:37
Irradiation 476 481 486 490 493 495 497 498 499 499 498 497
Time 12:52 13:07 13:22 13:37 13:52 14:07 14:22 14:37 14:52 15:07 15:22 15:37
Irradiation 495 493 490 486 481 476 469 462 452 442 429 414
Time 15:52 16:07 16:22 16:37 16:52 17:07 17:22
Irradiation 396 374 347 313 268 0 0
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . .
63
64 S. K. Sahu et al.
Table 4.3 Line and load data for 5-bus distribution system
Line data Load data Hosting capacity
Sl. From To R X P Q 70% load Full load
no. bus bus (Ohm) (Ohm) (kW) (kVA) (kVA) (kVA)
1 1 2 0.092 0.0477 0 0 605.0039 479.144
2 2 3 0.493 0.2511 100 60
3 3 4 0.366 0.1864 90 40
4 4 5 0.3811 0.1941 120 80
The solar irradiation input is taken from a specific geographical location [34]
(Table 4.1).
The irradiation from different geographical location can be extracted from geo-
graphical information system (GIS) data. These data are used as an input for
simulation to predict the exact behavior of the DER in different irradiation
conditions.
To analyze the realistic scenario of THD due to wind power plant, a realistic wind
data is taken from the east coast of India and depicted in Table 4.2.
Based on these inputs, both networks are simulated to analyze the effect of these
DERs on HC.
Pv Plant
Gen Bio-Fuel
Plant
Gen
Wind Power
Plant
Bus 5
Load
Gen Gen
Gen
Fig. 4.14 5-bus distribution system with all possible DG connections at all buses
Fig. 4.15 MATLAB-based 5-bus voltage indicator while injecting the calculated HC
66 S. K. Sahu et al.
The results from MATLAB are then imported to the real-time environment where
the THD is analyzed as shown in Figs. 4.16, 4.17, 4.18, 4.19, 4.20, 4.21, 4.22, and
4.23. Here in the plots, THD is plotted against time in seconds for different
irradiation conditions.
70 S. K. Sahu et al.
Fig. 4.24 THD at different buses at a wind speed of 6 m/s for 5- bus system
In this 5-bus system, solar photovoltaic sources are connected to all buses according
to the HC of the network. At this injection condition, THD is maximum at bus
2, which can be observed from Figs. 4.16, 4.18, 4.20, and 4.22. As bus 2 is
connected to maximum capacity line, its THD is slightly higher, but under the
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 71
Fig. 4.26 THD at different buses at a wind speed of 12 m/s for 5- bus system
standard in most of the time, which can be inferred from Figs. 4.17, 4.19, 4.21, and
4.23.
The results are also tabulated in Table 4.4 at different loading conditions. It can be
inferred that the THD violation is very minimal as per the standard.
72 S. K. Sahu et al.
Fig. 4.28 THD at different buses at a wind speed of 18 m/s for 5- bus system
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 73
Similar at all the nodes, wind power plant is connected in the considered 5-bus
network. As the uncertainty in the wind is very high, it can be observed from
Figs. 4.24, 4.25, 4.26, 4.27, 4.28, and 4.29 that the variation in the value of THD
for wind power plant is maximum among all other considered DERs.
74 S. K. Sahu et al.
The change in wind speed causes change in the THD quantity which can be
verified from Figs. 4.24, 4.25, 4.26, 4.27, 4.28, and 4.29 as well as from Table 4.5.
As wind profile is very dynamic in nature, thus, to have a stable output from wind
power plant, multiple converters are used. Due to this extensive use of converters,
the current and voltage waveforms get distorted and create THD in the utility.
Table 4.5 depicts the change in value of THD at various buses due to change in
wind profile and also shows the maximum number of violation of IEEE standard.
The biofuel-powered generators have the lowest THD contribution due to the
absence of the switching devices.
Biofuel-powered generators are the best substitute for the diesel generator due to
the low environmental impact. In these generators, normally, no power electronic
converters are used, so the distortion in the current and voltage waveform is minimal.
This can be seen from Fig. 4.30, where current and voltage THDs are shown in the Y-
axis. The values in the Y-axis are very small which resembles very low THD. As the
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 75
Table 4.7 Line and load data for IEEE 33-bus distribution system [35]
Line data Load data
Sl. no. From bus To bus R (Ohm) X (Ohm) P (kW) Q (kVA)
1 1 2 0.092 0.0477 0 0
2 2 3 0.493 0.2511 100 60
3 3 4 0.366 0.1864 90 40
4 4 5 0.3811 0.1941 120 80
5 5 6 0.819 0.7070 60 30
6 6 7 0.1872 0.6188 60 20
7 7 8 1.7114 1.2351 200 100
8 8 9 1.030 0.7400 200 100
9 9 10 1.04 0.7400 60 20
10 10 11 0.1966 0.0650 60 20
11 11 12 0.3744 0.1238 45 30
12 12 13 1.468 1.1550 60 35
13 13 14 0.5416 0.7129 60 35
14 14 15 0.591 0.5260 120 80
15 15 16 0.7463 0.5450 60 10
16 16 17 1.289 1.7210 60 20
17 17 18 0.732 0.5740 60 20
18 18 19 0.164 0.1565 90 40
19 19 20 1.5042 1.3554 90 40
20 20 21 0.4095 0.4784 90 40
21 21 22 0.7089 0.9373 90 40
22 22 23 0.4512 0.3083 90 40
23 23 24 0.898 0.7091 90 50
24 24 25 0.896 0.7011 420 200
25 25 26 0.203 0.1034 420 200
26 26 27 0.2842 0.1447 60 25
27 27 28 1.059 0.9337 60 25
28 28 29 0.8042 0.7006 60 20
29 29 30 0.5075 0.2585 120 70
30 30 31 0.9744 0.9630 200 600
31 31 32 0.3105 0.3619 150 70
32 32 33 0.341 0.5302 210 100
Hosting capacity Full load 4.8 MVA
70% load 3.5098 MVA
Gen
26 27 28 29 30 31 32 33
23 24 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
S/S
19 20 21 22
Fig. 4.32 MATLAB-based 33-bus voltage indicator while injecting the calculated HC
Assumptions:
1. All the sources are connected one by one to study the different real-time
conditions.
2. The grid harmonics are neglected to study the harmonic effect of the DERs.
3. The sources are placed only in one location, bus no. 2 since this bus is connected
to the maximum capacity line.
4. The test system is compensated.
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 77
The line data, load data, and base voltage with calculated HC are tabulated in
Table 4.7, and the corresponding single line diagram is shown in Fig. 4.31.
Similarly, to analyze the effect of THD, in case of -IEEE 33 bus distribution
system, firstly, keeping voltage limit and ampacity as constraints, the analysis is
done for HC. This can be verified from Fig. 4.32. It can also be verified that with the
increase in load, the voltage decreases.
All the variation of the voltages and its effect on THD for 33-bus distribution
system can be validated through the HIL results as in Figs. 4.33, 4.34, 4.35, 4.36,
4.37, 4.38, 4.39, and 4.40.
In this considered IEEE 33-bus distribution system, current THD at each node is
measured. This THD is mainly due to the SPV injection at bus 2 only. Here, bus 2 is
considered for the power injection as this bus is connected to the maximum power
line. The changing THD values can be observed in Figs. 4.33, 4.34, 4.35, 4.36. The
change in THD value in Fig. 4.33 for current at bus 2 can be observed; however, due
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 79
to scaling of the axes, bus 2 current plot appears to be smooth. Minute variation in
the THD values can be observed in all other buses.
Results in Figs. 4.33, 4.35, 4.35, and 4.36, is obtained for different irradation
conditions. In the foresaid Figs. 4.33, 4.35, 4.35, and 4.36, “X” axis shows the time
in seconds, and “Y” axis shows the THD magnitude. Table 4.8 summarizes all the
values of THD which are showcased in Figs. 4.33, 4.34, 4.35, and 4.36. From
Table 4.8, it can be observed that very small change in THD is taking place due to
the change in irradiation at different time of the day.
4.6.2.2 Case II: Wind Power Plant at Bus 2 in IEEE 33-Bus Distribution
System
Variation of wind speed throughout the day is a normal phenomenon that is faced by
the wind power plants. The effects of change in wind speed on THD values are
shown in Figs. 4.37, 4.38, and 4.39.
From the summary that is presented in Table 4.9, the values of THD are constant
at the lower wind speed. When the wind speed is higher, the value of THD changes
with variation in the load.
80 S. K. Sahu et al.
Fig. 4.37 THD at different buses at a wind speed of 6 m/s for 33 bus
Fig. 4.38 THD at different buses at a wind speed of 12 m/s for 33 bus
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 81
Fig. 4.39 THD at different buses at a wind speed of 18 m/s for 33 bus
Fig. 4.40 THD at different buses for 33-bus system with bio-fuel powered generator
The THD content due to the biofuel-powered generator is as shown in Fig. 4.40,
where “X” axis is the time axis and “Y” axis is the THD axis.
The summarized values in Table 4.10 show that the maximum value of THD is
0.4556. The diesel generator does not have variation in its power source, that is, fuel.
In Fig. 4.40, the THD values are stable, and at bus 2 which is connected to maximum
capacity line, the THD is slightly increasing at the initial stage due to starting of the
generator.
Although the THD values for the biofuel-powered generator are comparable to
other injections (SPV, WECS), it is of stable nature.
Thus, the overall HC as per THD analysis can be summarized as in Eq. (4.10):
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 83
(
HCCAL > HCBF PG > HCSPV > HCWECS , without compensation
¼ ð4:10Þ
HCCAL > HCBF PG > HCSPV =HCWECS , with compensation
Where,
HCCAL ¼ Calculated hosting capcity.
HCBF_PG ¼ Hosting capacity due to bio-fuel powered generator.
HCSPV ¼ Hosting capacity due to solar powered generator.
HCWECS ¼ Hosting capacity due to wind power generation system.
8
>
< HCMax , When all the sources are bio‐fuel powered genertors
HCTHD ¼ HCMin , When all the sources are wind powered generators
>
:
HCMax > HC > HCMin , When combination of sources are connected
ð4:11Þ
Based on the conclusion drawn in Eqs. (4.10) and (4.11), the HC of the network
can be planned to serve quality power without hampering the THD limits as well as
over-voltage and ampacity of the conductor so as to inject power optimally.
The sprawl of power system network has led to the injection of DERs in the form of
solar, wind, biofuel, etc. The planning of the size and quantity of these DERs in the
network is essential to obtain reliable and quality power. The DERs, especially solar
and wind, are totally nature dependent and are highly uncertain. Thus, the network
HC due to these DERs needs careful assessment. Due to the uncertain nature of these
sources, the THD, distorting the fundamental signal and violating the voltage and
ampacity limits, plays a significant role in HC assessment. The effect of THD on HC
at different irradiation conditions for SPV and different wind speed is analyzed both
in MATLAB and in real-time emulator using Typhoon HIL. The analysis is also
compared with biofuel-powered generators to monitor the varying effect of THD
on HC.
As the energy consumption per capita is exponentially growing, meeting the
energy demand is the biggest challenge. To optimally handle the emerging power
crunch problem, each bus should be injected with the maximum possible power, that
is, enhanced hosting capacity (EHC), which will not only solve the power crunch but
also encourage the individual to install DERs.
84 S. K. Sahu et al.
References
1. Harvey, R., Xu, Y., Qu, Z., & Namerikawa, T. (2017). Dissipativity-based design of local and
wide-area DER controls for large-scale power systems with high penetration of renewables. In
2017 IEEE conference on control technology and applications (CCTA), pp. 2180–2187.
2. Tsikalakis, A. G., & Hatziargyriou, N. D. (2007). Environmental benefits of distributed
generation with and without emissions trading. Energy Policy, 35(6), 3395–3409.
3. Munkhchuluun, E., & Meegahapola, L. (2017). Impact of the solar photovoltaic
(PV) generation on long-term voltage stability of a power network. In 2017 IEEE innovative
smart grid technologies – Asia (ISGT-Asia), pp. 1–6.
4. Bollen, M. H. J., & Rönnberg, S. K. (2017). Hosting capacity of the power grid for renewable
electricity production and new large consumption equipment. Energies, 10(9), 1325.
5. Zhou, Y., & Li, H. (2014). Analysis and suppression of leakage current in cascaded-multilevel-
inverter-based PV systems. IEEE Transactions on Power Electronics, 29(10), 5265–5277.
6. Stratford, R. P. (1980). Rectifier harmonics in power systems. IEEE Transactions on Industry
Applications, IA-16(2), 271–276.
7. Wang, J., Guo, Q., Yang, Z., Wang, Y., Bai, Z., & Ma, H. (2015). Modeling and research of
three-phase inverter powering nonlinear load., IEEE 24th International Symposium on Indus-
trial Electronics (ISIE), Buzios, 2015 pp, 331–336.
8. Sebaay, A. E. L., Ramadan, M., & Adma, M. A. A. (2017). Studying the effect of non-linear
loads harmonics on electric generator power rating selection. European Scientific Journal, ESJ,
13(18), 548.
9. Sakar, S., Balci, M. E., Aleem, S. H. E. A., & Zobaa, A. F. (2016). Hosting capacity assessment
and improvement for photovoltaic-based distributed generation in distorted distribution net-
works. In EEEIC 2016 – International conference on environmental and electrical engineering,
pp. 2–7.
10. Cundeva, S., Bollen, M., & Schwanz, D. (2017). Hosting capacity of the grid for wind
generators set by voltage magnitude and distortion levels. In Mediterranean conference on
power generation, transmission, distribution and energy conversion (MedPower 2016),
pp. 73–79.
11. Gupta, S. C., Kumar, Y., & Agnihotri, G. (2008). Optimal sizing of solar-wind hybrid system.
In ICTES, pp. 282–287.
12. Alturki, M., Khodaei, A., Paaso, A., & Bahramirad, S. (2018). Optimization-based distribution
grid hosting capacity calculations. Applied Energy, 219(March), 350–360.
13. Patsalides, M., et al. (2017). The effect of solar irradiance on the power quality behaviour of grid
connected photovoltaic systems. Renewable Energy and Power Quality Journal, 1(05),
323–330.
14. Al-Saadi, H., Zivanovic, R., & Al-Sarawi, S. F. (2017). Probabilistic hosting capacity for active
distribution networks. IEEE Transactions on Industrial Informatics, 13(5), 2519–2532.
15. Selvakumar, S., Madhusmita, M., Koodalsamy, C., Simon, S. P., & Sood, Y. R. (2019). High-
speed maximum power point tracking module for PV systems. IEEE Transactions on Industrial
Electronics, 66(2), 1119–1129.
16. Fathabadi, H. (2017). Novel maximum electrical and mechanical power tracking controllers for
wind energy conversion systems. IEEE Journal of Emerging and Selected Topics in Power
Electronics, 5(4), 1739–1745.
17. Sun, W., Harrison, G. P., & Djokic, S. Z. (2012). Distribution network capacity assessment:
Incorporating harmonic distortion limits. In IEEE power and energy society general meeting,
pp. 1–7.
18. Cantero, A. S. (2018). Modeling and characterization of a wind turbine emulator. In Advances
in renewable energies and power technologies (pp. 491–508). Elsevier https://www.elsevier.
com/books/advances-in-renewable-energies-and-power-technologies/yahyaoui/978-0-12-
812959-3.
4 Effect of Harmonics due to Distributed Energy Resources on Hosting Capacity. . . 85
19. Neill, S. P., Hashemi, M. R., Neill, S. P., & Hashemi, M. R. (2018). Offshore wind. In
Fundamentals of ocean renewable energy (pp. 83–106). Academic Press https://www.
elsevier.com/books/fundamentals-of-ocean-renewable-energy/neill/978-0-12-810448-4.
20. Huleihil, M., & Mazor, G. (2012). Wind turbine power: The Betz limit and beyond. In Advances
in wind power. IntechOpen.
21. Singh, G. K. (2013). Solar power generation by PV (photovoltaic) technology: A review.
Energy, 53, 1–13.
22. Poulios, V., Vrettos, E., Kienzle, F., Kaffe, E., Luternauer, H., & Andersson, G. (2014).
Optimal placement and sizing of battery storage to increase the PV hosting capacity of low
voltage grids. Master thesis. In International ETG Congress 2015; Die Energiewende – Blue-
prints for the new energy age, pp. 85–92.
23. Etherden, N., & Bollen, M. (2014). The use of battery storage for increasing the hosting
capacity of the grid for renewable electricity production. In Conference Innovations Security.
24. Kapusta, Ł. J., Sundell, J. P., & Teodorczyk, A. (2011). Liquid biofuels-promising energy
source for a small scale power plants. Silniki Spalinowe/Combustion Engines, 3(146), 1–6.
25. Shayani, R. A., Member, S., Aurélio, M., De Oliveira, G., & Member, S. (2011). Photovoltaic
Generation Penetration Limits in Radial Distribution Systems," in IEEE Transactions on Power
Systems, 3(26), 1625–1631.
26. Sharma, S. K., Chandra, A., Saad, M., Lefebvre, S., Asber, D., & Lenoir, L. (2017). Voltage
flicker mitigation employing smart loads with high penetration of renewable energy in distri-
bution systems. IEEE Transactions on Sustainable Energy, 8(1), 414–424.
27. Carnovale, D. J. (2014). Application of IEEE Std 519-2014 harmonic limits. In IEEE IAS
Atlanta Section, pp. 1–19.
28. C. Studies. (2014). IEEE_STD_519_1992vs2014.
29. Alalamat, F. (2015). Increasing the hosting capacity of radial distribution grids in Jordan.
30. Akram, U., Khalid, M., & Shafiq, S. (2018). Optimal sizing of a wind/solar/battery hybrid grid-
connected microgrid system. IET Renewable Power Generation, 12(1), 72–80.
31. Xu, L., Ruan, X., Mao, C., Zhang, B., & Luo, Y. (2013). An improved optimal sizing method
for wind-solar-battery hybrid power system. IEEE Transactions on Sustainable Energy, 4(3),
774–785.
32. Laksmi Kumari, R. V. S., Nagesh Kumar, G. V., Siva Nagaraju, S., & Babita Jain, M. (2018,
January). Optimal sizing of distributed generation using particle swarm optimization. In 2017
international conference on intelligent computing, instrumentation and control technologies
ICICICT 2017, pp. 499–505.
33. Alizadeh, S. M., Ozansoy, C., & Alpcan, T. (2016). The impact of X/R ratio on voltage stability
in a distribution network penetrated by wind farms. In Proceedings 2016 Australasian Univer-
sities Power Engineering Conference AUPEC 2016, pp. 1–6.
34. JRC photovoltaic geographical information system (PVGIS) – European Commission.
[Online]. Available: https://re.jrc.ec.europa.eu/pvg_tools/en/tools.html. Accessed 23 September
2019.
35. Rajaram, R., Sathish Kumar, K., & Rajasekar, N. (2015). Power system reconfiguration in a
radial distribution network for reducing losses and to improve voltage profile using modified
plant growth simulation algorithm with distributed generation (DG). Energy Reports, 1,
116–122.
Chapter 5
A Scenario-Based Approach for Storage
Capacity Determination to Improve
the Hosting Capacity of Distribution
Systems
5.1 Introduction
location of DERs and BESSs are given, and the unknown variables are the active and
reactive power of the BESSs. However, if BESSs are to be utilized for increasing the
HC of a system, their sizes and locations are unknown. Therefore, the methods in the
first category cannot be used directly for increasing the HC of a system. Neverthe-
less, the second category, which is based on the optimal allocation of BESSs could
have a better performance as the optimal size, location, and control action can be
identified to increase the HC.
Studies related to the allocation of BESSs have mainly focused on the following
aspects: authors in [15] allocated distributed energy storage systems to increase the
efficiency of distribution and transmission systems by using bilevel optimization. A
robust model to minimize the procurement cost by optimally allocating energy
storage systems and renewable energy sources was presented in [16]. Authors in
[17] allocated energy storage systems to maximize the benefit of the system operator.
The study exploited an accelerated Benders decomposition to solve the problem.
Energy storage systems were optimally allocated in [18] for maintaining the power
balance in active distribution networks considering the efficiency and life cycle of
the storage units. In [12], a methodology for allocating BESSs in a distribution
network is developed to minimize the wind energy curtailment as well as the annual
supply cost of the electricity. Authors in [19] defined a multi-objective optimization
to allocate the BESS by identifying a trade-off between the technical and economic
goals by minimizing the voltage deviation, network losses, feeder congestions, and
costs of supplying loads. In [13], a method for optimal placement, sizing, and control
of a community BESS to increase the HC of rooftop PVs in low-voltage
(LV) networks is proposed. Authors in [20] developed a stochastic optimization-
based method for BESS allocation and operation to ensure a desired voltage stability
margin in wind-intensive distribution systems. In [21], a convex optimization
approach has been used for BESS planning and operation in systems with high PV
penetration level. In that study, the minimum BESS required for mitigating the
technical impacts of high uptake of PVs has been estimated by using OpenDSS
simulation tool. Then, an optimization model has been used to identify the location
of BESSs. Authors in [22] proposed a heuristic strategy based on voltage sensitivity
in combination with an optimal power flow (OPF) framework to identify the optimal
number, size, and location of BESSs that prevent over-voltage and under-voltage in
an LV distribution network. To properly model the uncertainties associated with the
output of DERs and loads, authors in [6] proposed a scenario-based method to
identify the minimum required BESSs that prevents over-voltage in distribution
system with high PV penetration level. The proposed method in that study is
based on the voltage sensitivity of the scenarios associated with the customers’
injected power.
Nevertheless, to the best of the authors’ knowledge, the common drawback of the
aforementioned methods that utilize BESSs for HC improvement is that they did not
consider the uncertainty associated with the location and size of DERs. Hence, they
have only identified the required BESS capacity for a given scenario. Therefore, in
this chapter, a comprehensive method has been proposed to identify the minimum
required BESSs for increasing the HC of a system considering the uncertainties
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 89
associated with the DERs and loads. Then, an economic model has been developed
to assess the minimum required BESS from an economic perspective. Finally, the
performance of the developed models is assessed on an agricultural feeder in
Australia.
As shown in Fig. 5.1, energy storage technologies can be divided into five catego-
ries, as follows:
1. Electrochemical: These energy storage devices have been becoming more popu-
lar due to the decreasing trend of their costs as well as their high performances.
The main idea in these storages is storing electrical energy in the form of chemical
energy in batteries. The most popular electrochemical storages include lithium-
ion, lead-acid, nickel-cadmium, sodium-sulfur (NaS), and vanadium redox
(VRB) batteries.
2. Mechanical: These storages can store energy in different mechanical forms.
Mechanical storages include flywheel, pumped hydro, and compressed air energy
storages. Flywheels are electromechanical systems that can store mechanical
energy in the angular momentum of a rotating mass. Pumped hydro storages
are based on storing and releasing the potential energy of water. The charging
process consists of converting the electrical energy to potential energy by
pumping water from a lower to an upper reservoir. The discharging process
Energy storage
technologies
Sensible Pumped
NaS battery Methane
heat hydro
Lead-acid
Latent heat Flywheels
battery
Redox flow
battery
Nickel-
cadmium
consists of rotating a turbine by releasing water from the upper to lower reservoir
to generate electricity. Compressed air energy storage’s working principle is
accumulating energy by compressing and storing air in sealed tanks or under-
ground vessels. However, during the discharge process, the compressed air is
heated and released through a turbine to generate electricity.
3. Chemical: The basic idea in these storages is storing hydrogen as a compressed
gas or in the liquid form and synthesizing it to methane and finally generating
electricity by using fuel cells, combustion engines, or gas turbines.
4. Thermal: The working principle of these storages is storing energy in the form of
heat in different types of materials and generating electricity using steam turbines.
Generally, there are three types of thermal storages, namely, sensible heat, latent
heat, and thermo-mechanical storages. Please refer to [23] for more details.
5. Electrical: These storages can store electricity in the form of an electric field.
Electrical storages include supercapacitors, which have high efficiency and long
lifetime.
The electrochemical storages have been significantly progressing in recent years
[24]. The leading countries on electrochemical storages in terms of accumulated
rating power and number of installations are presented in Fig. 5.2, which is based on
all storage statuses, including operational, under construction, under repair,
announced, off-line, contracted, and decommissioned. Figure 5.3 shows the accu-
mulated rating power and the number of installed electrochemical BESSs that are
operational in 2015 [24]. As shown, the USA with a total capacity of 354 MW
1200 400
1000
350 NUMBER OF STORAGE PROJECTS
TOTAL CAPACITY (MW)
300
800
250
600 200
150
400
100
200
50
0 0
Fig. 5.2 Estimated battery storage capacity including operational, under repair, off-line, under
construction, contracted, and decommissioned by country in 2015
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 91
400 250
300
250 150
200
150 100
100
50
50
0 0
12.6
5.7
60
Share in (%)
12.3
11.1
40 6
62.6
48.7
20 39.2
0
Projects MWh Capacity MW Capacity
Lithium ion Lead acid Flow battery Sodium salt Nickel based Others
2444 MWh. As can be seen, 99.1% of all the considered projects belongs to lithium-
ion, lead-acid, flow battery, sodium-salt, and nickel-based categories.
subjected to
X
Pij,t ¼ Pjk,t þ pdj,t pgj,t 8j 2 DG, ð5:2Þ
k:j!k
X
Pij,t ¼ Pjk,t þ pdj,t 8j 2 fN ∖DG g, ð5:3Þ
k:j!k
X
Qij,t ¼ Qjk,t þ qdj,t qgj,t 8j 2 DG, ð5:4Þ
k:j!k
X
Qij,t ¼ Qjk,t þ qdj,t 8j 2 fN ∖DG g, ð5:5Þ
k:j!k
where Pij, t and Qij, t are the active and reactive power flow from bus (i) to bus ( j) at
the time step (t); pdi,t and qdi,t are the active and reactive load at bus (i) and time step
(t); pgi,t and qgi,t are the active and reactive generation at bus (i) and time step (t); ϕgi is
the power factor angle at bus (i) and time step (t); υi, t ¼ |Vi, t|2, where |Vi, t| is the
voltage magnitude at bus (i) and time step (t); Capgi and ηgi,t are the capacity and
capacity factor of the DER at the time (t) and bus (i); N gi,t represents the number of
DER at the time (t) and bus (i); N and DG are the set of all buses and the set of DG
locations, respectively; B is the set of all lines and υ represents the maximum
allowable voltage; rij + jxij is the impedance from bus (i) to ( j); (i, j) and j ! k
represents the line from bus (i) to bus ( j).
According to the third module, the HC of a distribution system could be divided
into three regions, as shown in Fig. 5.6. Region (A) shows the DER penetration
levels that do not cause any violation in technical constraints, regardless of DER
locations. Region (B) represents the DER penetration levels that do not cause any
technical violation when DERs are located at certain locations. It was demonstrated
A B C
in [25] that the probability distribution curve of the HC in the region (B) could be
approximated by a Gaussian-shape distribution. The start point of the probability
distribution curve is the border between regions (A) and (B). This start point is
defined as the minimum HC of the system. Region (C) represents all DER pene-
tration levels that would cause a violation of technical constraints, regardless of the
location of DERs. The border between regions (B) and (C) is defined as the
maximum HC. BESS is one of the options that can be used to increase both the
minimum and the maximum HC. Section 5.4 presents a methodology to identify the
minimum required BESS to increase the minimum and maximum HC up to a
certain level.
As shown in Fig. 5.6, the HC depends on the location of DERs in the region (B).
Therefore, any HC in the region (B) would cause a technical constraint violation for
some DER location scenarios. However, it is possible to resolve the constraint
violations by curtailing DERs. The optimization model for calculating the minimum
active power curtailment required to avoid technical constraint violation is as
follows:
X
maximize
|fflfflfflfflfflffl{zfflfflfflfflfflffl} pgi,t APCgi,t , ð5:11Þ
APCgi,t i2DG
subjected to
X
Pij,t ¼ Pjk,t þ pdj,t pgj,t þ APCgj,t 8j 2 DG, ð5:12Þ
k:j!k
X
Pij,t ¼ Pjk,t þ pdj,t 8j 2 fN ∖DG g, ð5:13Þ
k:j!k
X
Qij,t ¼ Qjk,t þ qdj,t qgj,t 8j 2 DG, ð5:14Þ
k:j!k
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 95
X
Qij,t ¼ Qjk,t þ qdj,t 8j 2 fN ∖DG g, ð5:15Þ
k:j!k
υ j,t ¼ υi,t 2 r ij Pij,t þ xij Qij,t 8ði, jÞ 2 B, ð5:16Þ
υi,t υ 8i 2 N , ð5:17Þ
APCgj,t pgj,t 8j 2 DG, ð5:18Þ
Solving the optimization models (5.11), (5.12), (5.13), (5.14), (5.15), (5.16), (5.17),
(5.18), (5.19), and (5.20) results in the active power curtailment required at each time
step. To identify the minimum required BESS for a scenario, the models (5.11),
(5.12), (5.13), (5.14), (5.15), (5.16), (5.17), (5.18), (5.19), and (5.20) should be
solved over the studied period. The maximum curtailed power during an hour over
the studied period is the minimum rating power required for the BESS for that
scenario. Similarly, the maximum curtailed energy during a day over the studied
period is the minimum energy rating of the required BESS. The rating power and
energy capacity of the required BESS should be identified for different scenarios that
would result in the same HC. Figure 5.7 represents the framework for identifying the
minimum BESS required to increase the minimum HC to a certain level (HCref) in
the region (B). The main idea of the framework is generating expansion scenarios
that would result in the total DER capacity of HCref. Then, the required rating power
and energy capacity of BESS should be identified. Finally, the required BESS is
estimated based on the designed BESS for the generated scenarios. In the following,
the steps of the framework are explained in detail.
• Step A: The focus of this step is addressing the uncertainty associated with the
size of DERs. Depending on the DER technology, there is a distribution curve for
the size of the units based on the historical statistics. For instance, Fig. 5.8
represents the distribution curve of PV systems in California. Generate a scenario
by sampling from the distribution curve of DERs. The sum of the samples should
be equal to the minimum HC level (HCref).
• Step B: This step is designed to address the uncertainty associated with the
location of DERs. To do so, the location of generated DER sample in Step B is
selected randomly, using a uniform distribution from the pool of potential
locations.
• Step C: After carrying out steps A and B, the optimization models (5.11), (5.12),
(5.13), (5.14), (5.15), (5.16), (5.17), (5.18), (5.19), and (5.20) should be used to
identify the required active power curtailment over the study period. Then, based
on the obtained curtailed power, determine the required BESS for the scenario. If
the curtailed power over the studied period is zero, the scenario is not a valid
scenario as it does not cause any technical constraint violation over the study
period. Thus, the following indicator for each scenario is defined to identify the
validity of each scenario.
96 M. S. S. Abad et al.
• Step A: Generate a scenario by sampling from the distribution curve of DERs. The
sum of the samples should be equal to the minimum HC level (HCref).
• Step C: Solve the optimization model (11)-(20) to identify the required active
power curtailment over the study period for the defined scenarion in steps A and B.
Capacity (kW)
Fig. 5.8 The distribution curve of PV systems based on California solar statistics [26]
1 Total APC > 0
ΘBESS ð jÞ ¼ ð5:21Þ
0 Otherwise
Z HCref
ProbCV ¼ pdf HC , ð5:22Þ
0
where pdfHC is the probability distribution function of the HC. It was also mentioned
that the active power curtailment over the study period is zero for some of the
generated scenarios in steps A and B. The probability of valid generated scenarios is
as follows:
PN BESS
j¼1 ΘBESS ð jÞ
PrAPC ¼ , ð5:23Þ
N BESS
where NBESS is the total generated scenarios. Repeat steps A, B, and C while the
following condition is held:
Estimating the required BESS using the proposed comprehensive method is based
on solving an optimization model over the study period for different expansion
scenarios, which means that the comprehensive method could have a considerable
computational burden. However, an important question is how to estimate the
required BESS without extensive computational burden? An approximated method
to estimate the BESS size is presented in this section. The proposed method only
uses the probability distribution function of the HC as well as the load and generation
profiles in the system. The approximation method is based on the maximum back-
feed in the system. To have a better understanding of the method, consider a
simplified distribution system as shown in Fig. 5.9. The back-feed can be calculated
as follows:
where sgen and PFgen are the capacity and power factor of DER and pbf represents the
active power that is injected to the upstream grid via distribution system transformer.
In this system, back-feed is the main limiting factor of the HC as it would result in an
over-voltage issue in bus (2). Thus, the HC of the system could be used to calculate
the maximum allowable back-feed as follows:
98 M. S. S. Abad et al.
1 2
pbf, qbf pgen, qgen
where smax
gen is the maximum DER capacity that could be accommodated in the
system, i.e., HC of the system. The maximum allowable back-feed can be used as
an index to estimate the minimum BESS required to increase the HC. Please refer to
[27] for more details on limitations of back-feed in real distribution systems. Any
DER with a higher capacity than smax safe
gen would result in a higher back-feed than pbf ,
hence the over-voltage problem. The excess back-feed power that would cause the
constraint violation can be calculated as follows:
time
Maximum excess
back-feed power
Fig. 5.10 Determination of excess back-feed power (red area) by comparison of maximum back-
feed power (dashed blue line) and back-feed power (solid green line)
• Step 3: Estimate the minimum required BESS based on the excess back-feed over
the study period. The maximum excess back-feed power during an hour over the
studied period is the minimum rating power required for the BESS. Similarly, the
maximum excess back-feed power during a day over the studied period is the
minimum energy rating of the required BESS.
This section explains the proposed method for performing the economic assessment
of the designed BESSs. It should be mentioned that the economic assessment in this
section is based on the annual costs and benefits of the BESSs.
Battery Inverter
Grid
Supporting structure
defined as the cost of all the auxiliary systems of a BESS including the trans-
former and supporting structures that are required to exchange the energy.
• Operation and maintenance cost: This term is defined to identify the annual cost
of maintenance and operation of the BESSs.
• Replacement cost: This term is designed to cover the cost of battery replacement
at the end of the life cycle of the battery. If the planning horizon is longer than the
life cycle of batteries, then this term should be considered.
The calculation of the total annual capital cost is as follows [28]:
EESS
BEC ¼ UCSC Total
, ð5:30Þ
Eff
where BEC and UCSC represent the total cost ($) and the unit cost ($/kWh) of
storage device, respectively; E ESS
Total is the designed energy capacity of the storage
device (kWh); and Eff represents the BESS efficiency.
where PEC and UCPE represent the total cost ($) and the unit cost ($/kW) of power
electronic device, respectively, and SESS
Total is the rated power of the BESS (kW).
E ESS
BOP ¼ UBOP Total
, ð5:32Þ
Eff
where BOP and UBOP represent the total cost ($) and the unit cost ($/kWh) of the
balance of plant, respectively. Thus, the total capital cost can be calculated as
follows:
where TCC denotes the total capital cost. Then, the total annual capital cost is as
follows:
where
IRð1 þ IRÞhzn
CAF ¼ , ð5:35Þ
ð1 þ IRÞhzn 1
where IR is the annual interest rate and hzn is the planning horizon. Further, the
operation and maintenance cost can be calculated as follows:
where AOMC and FOM denote the total annual cost ($/year) and the fixed unit cost
($/(kW.year)) of the maintenance and operation of BESS, respectively. Finally, the
replacement cost can be calculated as follows:
EESS
ARC ¼ AFactor Total
, ð5:37Þ
Eff
where
h i
AFactor ¼ F factor ð1 þ IRÞr þ ð1 þ IRÞ2r þ CAF, ð5:38Þ
CycleTotal
r¼ , ð5:39Þ
CycleDay N year
where Ffactor represents the future value of the battery replacement cost ($/kWh);
CycleTotal and CycleDay denote the number of charge/discharge cycles of the storage
during its life cycle and during a day, respectively; and Nyear is the number of
operating days in a year for the storage system. Therefore, the total annual cost
(TAC) can be calculated as follows:
Compensation of the curtailed energy highly depends on the contract between the
off-taker and the generator owner. Usually, the compensation of active power
102 M. S. S. Abad et al.
XX
8760
CostAPC ¼ T grid
t APCgi,t , ð5:41Þ
i2DG t¼1
where T grid
t is the network electricity price at time (t).
The focus of this section is identifying the benefits of the integration of BESSs in the
system for the utility. Two obvious options could benefit the utility, which are the
reduction in the network losses as well as the energy arbitrage between peak and
off-peak periods. Both of these objectives can be achieved by minimizing the costs
of supplying energy from a utility perspective. Therefore, the objective function can
be defined as follows:
X X
8760 XX
8760
Minimize costs ¼ T grid
t r ij ℓij,t i,t T t :
pESS grid
ð5:42Þ
ði, jÞ2B t¼1 i t¼1
The decision variables are the optimal size and location of BESSs. The constraints of
the optimization problem include the BESS installation and operation constraints,
network constraints, and DER constraints, which are given as follows:
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 103
X
BESS
i BESS
total , ð5:43Þ
i
BESS
i SESS,
i
min
SESS
i BESS
i SESS,
i
max
, ð5:44Þ
BESS
i E ESS,
i
min
EESS
i BESS
i E ESS,
i
max
, ð5:45Þ
X
SESS
i SESS
Total , ð5:46Þ
i
X
E ESS
i E ESS
Total , ð5:47Þ
i
BESS
i 2 f01g, ð5:48Þ
where BESStotal is the total number of locations in which BESSs could be installed;
Bi represents the binary variable to allocate BESS at bus (i); SESS
ESS
i and EESS
i
denote the power rating and energy capacity of a given BESS at bus (i),
respectively; SESS,i
max
and EESS,
i
max
represent the maximum power rating and
energy capacity of a given BESS at bus (i); and SESS, i
min
and EESS,
i
min
denote the
minimum power rating and energy capacity of a given BESS at bus (i). Con-
straints (5.46) and (5.47) limit the total power rating and energy capacity of the
BESS to the minimum required BESS, which was identified in Sect. 5.4.
2. BESS operation constraints [19]
X
t
i,t 0 Δt þ Lossi,t 0
pESS E ESS 8t 2 H , ð5:49Þ
ESS
i,t 0 ¼0
t 0 ¼0
X
t
pESS
i,t 0 Δt Lossi,t 0
ESS
EESS
i E ESS
i,t 0 ¼0 8t 2 H , ð5:50Þ
t 0 ¼0
SESS
i pESS
i,t Si
ESS
8t 2 H , ð5:52Þ
εmin
i E ESS
i E ESS
i,t εi
max
E ESS
i 8t 2 H , ð5:53Þ
i,t ¼ 1 ηi
LossESS ESS
pESS
i,t
8t 2 H , ð5:54Þ
i,t pi,t1 þ pi
pESS 8t 2 H , ð5:55Þ
ESS RURD
i,t pi,t1 pi
pESS 8t 2 H , ð5:56Þ
ESS RURD
where H and Δt are the set of time steps and the length of each time step,
i,t 0 ¼0 represents the initial stored energy in BESS at bus (i); εi ,
respectively; EESS min
104 M. S. S. Abad et al.
εmax
i , and ηi
ESS
denote the minimum and maximum allowable stored energy level
as well as the efficiency of the BESS at bus (i), respectively; and pRURD i is the
charging/discharging rate of BESS at bus (i). Constraint (5.49) is defined to make
sure that the sum of all stored and taken energy over all the steps of the study
period is less than the initial stored energy. Similarly, constraint (5.50) is defined
to keep the sum of all stored and taken energy over all the steps of the study
period below the initial available capacity of BESS. Constraint (5.51) models the
relationship between the state of charge and the output power of BESS at the time
(t) and bus (i). Constraint (5.52) keeps BESS power below their rating power;
constraint (5.53) is defined to set the minimum and maximum state of charge of
BESS at bus (i); constraint (5.54) is defined to model the energy losses in each
BESS at each time step; constraints (5.55) and (5.56) define the maximum
charging and discharging rate of BESSs. Finally, constraint (5.54) can be written
in a linear form by defining new variables as follows:
0 XX i,t SESS,
i
max
, ð5:58Þ
0 YY i,t SESS,
i
max
, ð5:59Þ
0 LossESS
i,t 1 ηi
ESS
ðXX i,t YY i,t Þ 2 SESS,
i
max
DD2,i,t , ð5:60Þ
i,t 1 ηi
0 LossESS ESS
ðYY i,t XX i,t Þ 2 SESS,
i
max
DD1,i,t , ð5:61Þ
DD1,i,t þ DD2,i,t ¼ 1, ð5:62Þ
DD1,i,t , DD2,i,t 2 f0, 1g, ð5:63Þ
where DD1,i,t and DD2,i,t are binary variables that identify the charging/
discharging status of the BESS. DD1,i,t ¼ 1 when the BESS is discharging, and
DD2,i,t ¼ 1 when the BESS is charging.
3. Network constraints
These constraints include the power flow as well as over-voltage constraints. Any
other technical constraints such as thermal capacity could be also defined as a
network constraint. The network constraints considered in this study are as
follows:
X
Pij,t ¼ Pjk,t þ pdj,t pgj,t pESS
j,t þ r ij ℓ ij,t 8j 2 DG, ð5:64Þ
k:j!k
X
Pij,t ¼ Pjk,t þ pdj,t pESS
j,t þ r ij ℓ ij,t 8j 2 fN ∖DG g, ð5:65Þ
k:j!k
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 105
X
Qij,t ¼ Qjk,t þ qdj,t qgj,t þ xij ℓij,t 8j 2 DG, ð5:66Þ
k:j!k
X
Qij,t ¼ Qjk,t þ qdj,t þ xij ℓ ij,t 8j 2 fN ∖DG g, ð5:67Þ
k:j!k
υ j,t ¼ υi,t 2 r ij Pij,t þ xij Qij,t þ r 2ij þ x2ij ℓ ij,t 8ði, jÞ 2 B, ð5:68Þ
P2ij,t þ Q2ij,t
ℓ ij,t ¼ 8ði, jÞ 2 B, ð5:69Þ
υi,t
υ υi,t υ 8i 2 N , ð5:70Þ
X
k X
k
t¼ λv t v , λv ¼ 1, λv 0 8v 2 f0, 1, 2, . . . , k g, ð5:71Þ
v¼0 v¼0
X
k
bf ðt Þ ¼ λvbf ðt v Þ, ð5:72Þ
v¼0
where t0 < t1 < t2 < . . . < tk represent the ending points of the pieces. Further,
λv 8 v 2 {0, 1, . . ., k} are special order set of type 2 variables. This means that only
two adjacent λv are non-zero. This restriction can be removed by using binary
variables as follows:
0 λ0 w 0 , ð5:73Þ
0 λv wv1 þ wv , 8v 2 f1, 2, . . . , k 1g, ð5:74Þ
0 λk wk1 , ð5:75Þ
106 M. S. S. Abad et al.
X
k1
wv ¼ 1, ð5:76Þ
v¼0
This linearization technique is only valid for single variable functions. However,
ℓ ij, t in constraint (5.69) is a function of (P, Q, v). Nevertheless, if it is supposed
that v is fixed, then constraint (5.69) turns into a separable function. A function
such as f(Var1, Var2, . . ., Varn) is called separable if it could be expressed as a
summation of n single variable function such as f1(Var1), f2(Var2), . . ., fn(Varn).
Thus, each single variable function can be approximated with a piecewise linear
function. Please refer to [27] for more information regarding approximating the
nonlinear term (5.69) with a piecewise linear function.
4. DER constraints
where ηgi,t is the capacity factor of DER at time (t) and Capgi and ϕgi are the
capacity and power factor angle of DER at bus (i).The optimization model for
allocating the BESSs to minimize the costs of supplying loads is a mixed-integer
linear programming (MILP).
The total benefits for DER owner as well as the utility are compared to the annual
cost of installing the required BESS. If the annual benefit is higher than the annual
costs, using BESSs to improve the HC is justifiable.
There are different electrical energy storages, including lead-acid (LA), sodium-
sulfur (NaS), vanadium redox (VRB), and zinc/bromine (ZnBr). The comparison of
different BESSs is presented in Table 5.1 based on the development status, energy,
and power capital costs [30]. As can be seen, there is a considerable difference
between the capital cost of these BESS technologies, which can affect the feasibility
of the required BESS capacity from an economic perspective.
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 107
Table 5.1 Comparison of the BESS technologies based on cost and development status indexes
Energy capital cost Power capital cost
Technology Development status (US $/kWh) (US $/kW)
Lead-acid (LA) Mature 200–400 200–600
Sodium-sulfur (NaS) Commercialized 275–500 1000–3000
Vanadium redox (VRB) Early commercialization 600–1500 150–1000
stage
Zinc/bromine (ZnBr) Developed 200–400 150–1000
Energy storage installation cost ($/kWh)
1000
900
800
700
600
500
400
300
200
100
0
2016 2030
Lead acid (Flooded) Lead acid (VRLA)
Sodium-sulfur Sodium-Nickel chloride
Flow (VRFB) Flow (ZBFB)
Lithium-ion (NMC/LMC) Lithium-ion (LFP)
Fig. 5.12 Installation cost estimation of battery storages in 2016 and 2030
It should be mentioned that the capital cost of BESSs has been decreasing, which
means that although a BESS technology might not be economically feasible now, it
may become feasible in the future. Figure 5.12 shows the energy capital cost
estimation for different battery technologies in 2016 and 2030 [31]. As can be
seen, it is expected to have up to 66% decrease in capital costs by 2030. Similarly,
it is expected to experience an increase in the efficiency of battery storages. Fig-
ure 5.13 presents the expected change in the energy capital cost as well as the
efficiency of battery storage technologies in 2016 and 2030 [31].
108 M. S. S. Abad et al.
Fig. 5.13 Estimated changes in the installation costs as well as the efficiency of battery storages
from 2016 to 2030
In this section, simulations are carried out to assess the performance of the proposed
methodology. Initially, the effectiveness of the proposed method for identifying the
minimum BESS capacity is evaluated. Then, the economic feasibility of the mini-
mum required BESS for different technologies is discussed. Finally, the effect of
DER technologies on the minimum required BESS is assessed.
Figure 5.15 represents the normalized load profile, which is used to model the load
variation in the test network. The load profile is derived from the data made available
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 109
1 Substation 80
166 102 91
164 165 113
163 134 33
160 161
162 167 44
14 15
159 11 55
12 13 155 153
158 10 22 77
8
157 154 152 81
147 148 149 150 151 63 62 61
146 82
16 17 71 70 69
144 83
65 64
142 143
141 84
73 72 60 59
18 19 20
140 57 85
138 139
137 58 86
21 23
136 87
24 25
135 92 90
48 50 116 52 104
117 53
by the Australian Energy Market Operator (AEMO) [33]. The normalized PV and
wind profiles are presented in Figs. 5.16 and 5.17, which are derived from [34, 35],
respectively. The distribution of PV and wind sizes are derived from [25].
110 M. S. S. Abad et al.
0.9
0.8
Load profile
0.7
0.6
0.5
0.4
0.3
300
200 20
15
100 10
5 Hour
Day
Fig. 5.15 Normalized hourly load profile for New South Wales from 1/1/2014 to 30/12/2014
0.8
0.6
PV profile
0.4
0.2
300
200 20
15
100 10
5
Day Hour
Fig. 5.16 Normalized hourly PV profiles for New South Wales from 1/1/2014 to 30/12/2014
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 111
Fig. 5.17 Wind production profile for New South Wales from 1/1/2014 to 30/12/2014
Table 5.2 represents all the parameters required for calculating the annual instal-
lation cost of BESS for different technologies [28]. Table 5.3 shows the electricity
price rate for two different schemes from Origin Energy [36]. Peak period mentioned
in Table 5.3 is from 1 to 8 pm; shoulder period is from 7 am to 1 pm and from 8 to
10 pm; and off-peak period is from 10 to 7 am.
112 M. S. S. Abad et al.
Fig. 5.18 PV HC
probability curve for the
agricultural feeder
CDF
5.7.3 Hosting Capacity Discussion
PV and wind HC of the test system are estimated based on the HC calculation
method presented in Sect. 5.3. Figure 5.18 presents PV HC curve for the test system,
which can be approximated using the Gaussian-shape distribution as (5.80) with
α ¼ 0.56, μ ¼ 19.97 MW, and σ ¼ 6.42 MW. The minimum PV HC for the test
system is 3.198 MW.
α ðxμÞ2
HC α N μ, σ 2 ¼ pffiffiffiffiffiffiffiffiffiffi e 2σ2 : ð5:80Þ
2πσ 2
Figure 5.19 presents wind HC curve for the test system, which can be approxi-
mated using the Gaussian-shape distribution as (5.80) with α ¼ 0.956,
μ ¼ 37.21 MW, and σ ¼ 0.7092 MW. The minimum wind HC for the test system
is 2.92 MW. As it can be seen, the HC distribution for wind is significantly different
from PV HC. This is mainly due to the difference between the wind and PV
generation profiles as well as the distribution of wind and PV sizes.
The performance of the proposed methods for estimating the minimum required
BESS is examined in this section. The minimum HC for PV and wind generation is
3.198 MW and 2.92 MW, respectively. Table 5.4 presents the estimated required
BESS to increase the minimum HC to 4 MW.
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 113
0.6
Histogram
Gaussian PDF
0.5
0.4
Probability
0.3
0.2
0.1
0
5 10 15 20 25 30 35
Hosting Capacity (MW)
Table 5.4 Minimum required BESS to increase the HC to 4 MW for PV and wind technologies
using the comprehensive and approximated methods
Comprehensive method Approximated method
Energy capacity Rating power Energy capacity Rating power
DER technology (MWh) (MW) (MWh) (MW)
PV 5.275 0.6697 5.515 0.75
Wind 9.879 0.5212 10.279 1.07
This section discusses the economic feasibility of the estimated BESS that is
required for different technologies. To do so, the total annual cost of the required
BESS for different technologies is compared with the annual benefits of the
corresponding BESS. Table 5.5 shows the total annual cost for the required BESS
estimated using the proposed comprehensive method for PV and wind technologies,
respectively. As can be seen, using the ZnBr battery would result in the lowest total
annual costs for both PV and wind technologies, while the valve-regulated lead-acid
(VRLA) battery has the highest total annual costs.
If the total annual profit of the designed battery is higher than the minimum
annual costs, the designed BESS is economically feasible. Tables 5.6 and 5.7 present
the optimal allocation of the required BESS in the system for both PV and wind
technologies when the electricity price follows a flat rate. As can be seen, the
required BESS capacity is not distributed equally. Further, the optimal BESS
location for PV technology is different from that of wind technology.
Tables 5.8 and 5.9 present the optimal allocation of the required BESS in the
system for both PV and wind technologies when the electricity price follows a ToU
rate. Observe that the required BESS capacity is not distributed equally. Further, the
optimal BESS location for PV technology is different from that of wind technology.
Another important factor that has a considerable impact on the allocation of the
Table 5.5 Total annual costs of the required BESS for different technologies in US $ (2014)
BESS technology
DER technology LA VRLA NaS ZnBr VRB
PV 902,323 2,124,994 1,696,292 370,760 1,056,548
Wind 1,656,803 3,953,947 3,036,043 657,629 1,964,038
Table 5.6 Optimal allocation of the required BESS in the test system for PV generation when the
electricity price follows the flat rate
BESS1 BESS2
Location 149 159
Energy capacity (MWh) 3.6 1.675
Rating power (MW) 0.3 0.3697
Table 5.7 Optimal allocation of the required BESS in the test system for wind generation when the
electricity price follows the flat rate
BESS1 BESS2
Location 151 162
Energy capacity (MWh) 5.4 4.479
Rating power (MW) 0.2643 0.2569
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 115
Table 5.8 Optimal allocation of the required BESS in the test system for PV generation when the
electricity price follows the ToU rate
BESS1 BESS2
Location 149 159
Energy capacity (MWh) 1.675 3.6
Rating power (MW) 0.297 0.3727
Table 5.9 Optimal allocation of the required BESS in the test system for wind generation when the
electricity price follows the ToU rate
BESS1 BESS2
Location 151 162
Energy capacity (MWh) 5.4 4.479
Rating power (MW) 0.2468 0.2744
200
Profit in US $ (Thousands)
180
160
140
120
100
80
60
40
20
0
PV Wind
Flat rate ToU rate
Fig. 5.20 Utility annual profit from installing the required BESS for PV and wind generation for
both flat and ToU electricity price rates
required BESS for both PV and wind technologies is the electricity price. Consid-
ering Tables 5.6 and 5.8, it can be noted that the optimal location of BESSs for both
flat and ToU electricity prices is the same. However, the energy capacity and the
rating power of the optimal BESSs are different.
Figure 5.20 demonstrates the annual benefit that utility would gain from installing
the BESSs for both PV and wind technologies. As can be seen, the annual profit of
the utility with ToU tariff is higher than that of the flat tariff. Further, the BESS
benefit for PV generation is much higher than that of wind generation in the test
system. Moreover, the total annual profit of the utility is below the total annual costs
of the minimum required BESS for different storage technologies, which might
convey that BESS is not an economically feasible solution for increasing the
116 M. S. S. Abad et al.
4500
4000
3500
3000
US $ (Thousands)
2500
2000
1500
1000
500
0
Utility DERs Utility DERs total Total LA cost VRLA NaS cost ZnBr VRB
profit profit profit profit profit Profit cost cost cost
with flat with flat with with with flat with
rate rate ToU rateToU rate rate ToU rate
PV Wind
Total annual profit of BESS for PV Total annual profit of BESS for wind
Fig. 5.21 Total annual profit compared with the total annual cost for different BESS technologies
HC. Nevertheless, the utility benefit is only a part of the BESS benefit. DER owners
also gain some benefit from avoiding active power curtailment. Figure 5.21 shows
the total annual profit of both utility and DER owners compared with the total annual
costs of different BESS technologies. As it can be seen, the total annual profit of the
BESS with both flat and ToU electricity price rates is still lower than the total annual
costs of all BESS technologies, which means that although BESSs can increase the
HC of the system, they are not a feasible option from an economic perspective.
However, it was also noted that the total annual profit of BESS for PV technology,
when the electricity price followed ToU rate, was 325,353.00$, which is consider-
ably close to the total annual cost of ZnBr battery technology, i.e., 370,760.00$.
Another considerable point in Fig. 5.21 is that the total annual profit of BESS for
wind generation is considerably lower than the total annual cost of all BESS
technologies. This implies that there are a few days that the spilled wind energy is
very high, but the spilled wind energy for the rest of the days is at a low level.
Finally, it should be mentioned that although some of BESS technologies are
economically infeasible now, it does not mean that they will continue to be infeasible
in the future. As it was mentioned in Sect. 5.6, it is predicted that the price of BESSs
would decrease and their efficiency would increase. Figure 5.22 presents the total
annual profit of both utility and DER owners compared with the expected annual
costs of different BESS technologies in 2030. As can be seen, ZnBr battery
5 A Scenario-Based Approach for Storage Capacity Determination to Improve. . . 117
3500
3000
2500
US $ (Thousands)
2000
1500
1000
500
0
Utility DERs Utility DERs total Total LA cost VRLA NaS cost ZnBr VRB
profit profit profit profit profit Profit cost cost cost
with flat with flat with with with flat with
rate rate ToU rateToU rate rate ToU rate
PV
Wind
Total annual profit of BESS for PV with ToU rate
Total annual profit of BESS for wind with ToU rate
Total annual profit of BESS for PV with flat rate
Fig. 5.22 Expected total annual profit in 2030 compared with the total annual cost for different
BESS technologies
technology is expected to be economically feasible in 2030 with both flat and ToU
electricity rates for PV technology. However, even the future decrease in the price of
energy storages would not make BESSs a feasible option for increasing the HC for
wind technology.
5.8 Summary
References
20. Jalali, A., & Aldeen, M. (2019). Risk-based stochastic allocation of ESS to ensure voltage
stability margin for distribution systems. IEEE Transactions on Power Systems, 34(2),
1264–1277.
21. Li, Q., Ayyanar, R., & Vittal, V. (2016). Convex optimization for DES planning and operation
in radial distribution systems with high penetration of photovoltaic resources. IEEE Trans-
actions on Sustainable Energy, 7(3), 985–995.
22. Giannitrapani, A., et al. (2017). Optimal allocation of energy storage systems for voltage control
in LV distribution networks. IEEE Transactions on Smart Grid, 8(6), 2859–2870.
23. Dusonchet, L., et al. (2019). Technological and legislative status point of stationary energy
storages in the EU. Renewable & Sustainable Energy Reviews, 101, 158–167.
24. Telaretti, E., & Dusonchet, L. (2017). Stationary battery technologies in the US: Development
trends and prospects. Renewable & Sustainable Energy Reviews, 75, 380–392.
25. Seydali Seyf Abad, M., et al. (2018). Probabilistic assessment of hosting capacity in radial
distribution systems. IEEE Transactions on Sustainable Energy, 9(4), 1935–1947.
26. Dubey, A., Santoso, S., & Maitra, A. (2015). Understanding photovoltaic hosting capacity of
distribution circuits. In 2015 IEEE Power & Energy Society General Meeting, Denver, CO,
pp 1–5.
27. Seydali Seyf Abad, M., Ma, J., & Han, X. (2019). Distribution systems hosting capacity
assessment: Relaxation and linearization. In Smart power distribution systems (pp 555–586).
Cambridge, Massachusetts, USA: Elsevier.
28. Poonpun, P., & Jewell, W. T. (2008). Analysis of the cost per kilowatt hour to store electricity.
IEEE Transactions on Energy Conversion, 23(2), 529–534.
29. Lori Bird, J. C., & Wang, X. (2014). Wind and solar energy curtailment: Experience and
practices in the United States. Golden: National Renewable Energy Laboratory.
30. Akinyele, D., Belikov, J., & Levron, Y. (2017). Battery storage technologies for electrical
applications: Impact in stand-alone photovoltaic systems. Energies, 10(11), 1760–1799.
31. Pablo Ralon, M. T., Ilas, A., Diaz-Bone, H., & Kairies, K.-P. (2017). Electricity storage and
renewables: Costs and markets to 2030. Australia: IRENA.
32. Berry, A., et al. (2013). National feeder taxonomy; describing a representative feeder set for
Australian networks. Newcastle: CSIRO.
33. AEMO, National transmission network development plan 2015. 2015.
34. Pfenninger, S., & Staffell, I. (2016). Long-term patterns of European PV output using 30 years
of validated hourly reanalysis and satellite data. Energy, 114, 1251–1265.
35. Staffell, I., & Pfenninger, S. (2016). Using bias-corrected reanalysis to simulate current and
future wind power output. Energy, 114, 1224–1239.
36. Energy, O. (2018). NSW residential energy price fact sheet. NSW, Australia: Origin Energy.
Chapter 6
Impact of Capacity Withholding
on Hosting Capacity Analyzing
6.1 Introduction
S. Salarkheili (*)
Faculty of Electrical Engineering, Shahid Beheshti University, A.C, Tehran, Iran
Tehran Regional Electric Company (TREC), Tehran, Iran
e-mail: s_salarkheil@sbu.ac.ir
M. SetayeshNazar
Faculty of Electrical Engineering, Shahid Beheshti University, A.C, Tehran, Iran
hosting capacity (HC) analyzing can cause to select proper capacities and locations
of DGs to mitigate capacity withholding in power markets. In this chapter, capacity
withholding analyzing of GenCos is described and different indices are proposed. In
Sect. 6.2, the effect of generation limits on the capacity withholding of GenCos is
analyzed. In the next section, the capacity withholding of GenCos is evaluated in
transmission-constrained power markets. In Sect. 6.4, the mitigation of capacity of
withholding of GenCos and the impact of DGs on capacity withholding are
described.
Noncollusive game models such as Cournot and supply function equilibrium (SFE)
models are broadly used for analyzing of market power [3–30]. The Cournot model
is used in this chapter to represent the strategic GenCos’ behavior in the presence of
elastic demands. Therefore, the theory of Cournot-NE (CNE) is used to find the
market equilibrium considering the GenCos’ capacity withholding.
Suppose in the market there are N GenCos and each GenCo has the following cost
function:
where qi is the quantity generated by GenCo i and ai and bi are the coefficients of the
GenCo’s cost function with ai 0 and bi 0.
The aggregate demand function is
λ ¼ αY þ β
ð6:2Þ
α>0
P
N
where Y ¼ yi and the transmission loss is ignorable. α and β are parameters of the
i¼1
aggregate demand function and λ is the electricity market price.
Each GenCo seeks in the market to maximize its profit (Ωi), and the optimization
problem faced by GenCo i is determined as follows:
max Ωi ¼ qi λ c i
ð6:3Þ
s:t: qmin
i qi qmax
i
where qmax
i and qmin
i are upper and lower production constraints of GenCo i,
respectively. The Lagrange function of the problem (6.3) is
Li ¼ qi λ ci þ θmax
i qi qmax
i θmin
i qi qmin
i ð6:4Þ
where θmax
i and θmin
i are Lagrange multipliers associated with lower and upper
production limits of GenCo i, respectively.
Differentiating the (6.4) with respect to qi and setting it as zero obtains
∂Li ∂λ
¼ qi þ λ ai qi bi þ θmax θmin ¼0 ð6:5Þ
∂qi ∂qi i i
In a perfectly competitive market, all GenCos are not price maker and the term
∂π/∂yi in (6.5) is equal to zero. But in an oligopoly market, GenCos are not price
taker and they are price taker and ∂π/∂yi 6¼ 0 in (6.5) can be yields. So qi of GenCo
i under perfect competition can be obtained by
λ bi θmax þ θmin
qi ¼ i i
ð6:6Þ
ai
According to (6.5) and by using (6.2), qi each GenCo under Cournot model can
be obtained by
λ bi θmax þ θmin
qi ¼ i i
ð6:7Þ
ð α þ ai Þ
124 S. Salarkheili and M. SetayeshNazar
Many indices have been proposed for market power assessment. There are different
categorizations for indices of market power. One of them categorizes indices of
market power into comparison and concentration indices. In comparison indices, the
oligopoly market’s results are compared with those of perfectly competitive market,
and concentration indices are based on the GenCos’ market shares. As another
categorization, the indices of market power can be classified into ex ante indices,
based on the simulation tools and market models, and ex post indices, in which the
results of real market are used. As a third categorization, there are global and nodal
indices. The nodal indices are used for analyzing of market power in each bus of
power system, and global indices assess the market power in the whole system.
Nodal indices are appropriate for transmission-constrained power markets [26–29].
A set of comparison indices is used in this chapter to analyze the capacity withhold-
ing of GenCos. In the proposed indices, the results of oligopoly market are compared
with those of two types of markets: perfectly competitive market with the nodal
prices of oligopoly market and perfectly competitive market. The superscripts e, p,
and p(e) are used for the values associated with the equilibriums of oligopoly market,
perfectly competitive market, and perfectly competitive market with the nodal prices
of oligopoly market, respectively. It should be noted that in this section, the
transmission constraints are ignored and nodal prices are equal. Based on (6.6) and
(6.7), the generation of GenCo i under perfectly competitive market (qpi ) can be
determined by
where θpi max and θpi min are Lagrange multipliers associated with upper and lower
production constrains of GenCo i under perfectly completive market, respectively.
The generation of GenCo under oligopoly market (qei ) can be written as
where θei min and θei max are Lagrange multipliers associated with lower and upper
production limits of GenCo i under oligopoly market, respectively. The generation
of GenCo under perfectly competitive market with the nodal prices of oligopoly
market can be also written as
Dqidistortion
e
qei q
p qp(p )i q
i
pðeÞ
Δqwithheld
i ¼ qi qei ð6:11Þ
Δqdistortion
i ¼ qpi qei ð6:12Þ
Δλdistortion ¼ λe λp ð6:13Þ
ΔQdistortion
DWI ¼ ð6:14Þ
ΔQwithheld
As shown in Fig. 6.1, λ p is the competitive market clearing price and qpi is the
GenCo i’s competitive generation. The output of GenCo i decreases to qei , and the
pðeÞ
market clearing price increases to λe due to the capacity withholding. qi is the
competitive output which would be generated at λ by GenCo i. In (6.14),
e
P
N P
N
ΔQwithheld ¼ Δqwithheld
i and ΔQdistortion ¼ Δqdistortion
i . The capacity withholding
i¼1 i¼1
indices are reported in Fig. 6.1. DWI ranges from 0 to 1 as it analyzes the market
power of GenCos from the perspective of their potential ability for capacity with-
holding. DWI indicates a decrement in the ability of GenCos for capacity withhold-
ing when it approaches 1. In contrast, when this index is lower, it is more likely that
producers have a high potential ability for capacity withholding. It can be seen that
when market power is low, the GenCos will have less incentive to raise the market
pðeÞ
price. It means that λe, qi , and Δqwithheld
i approach λ p, qpi, and Δqdistortion
i ,
respectively. It means that DWI approaches 1, as a result of the changes in both
Δqwithheld
i and Δqdistortion
i .
126 S. Salarkheili and M. SetayeshNazar
α e Δθwithheld
Δqwithheld ¼ q þ i ð6:15aÞ
i αi i ai
pðeÞ max pðeÞ min
Δθwithheld
i ¼ θi þ θei min θi θei max ð6:15bÞ
α e Δθdistortion
Δqdistortion
i ¼ qi ΔQdistortion þ i ð6:16aÞ
ai ai
X
N
Δqdistortion
i ¼ ΔQdistortion ð6:16bÞ
i¼1
Δθdistortion
i ¼ θpi max þ θei min θpi min θei max ð6:16cÞ
It can be seen that the value of capacity distortion index may be negative or
positive. It depends on slope of the marginal cost functions of GenCos and the slope
of aggregate demand function. It is to be mentioned that α/ai ratio is significant for
analyzing the ability of GenCos for capacity withholding. For GenCos with large
ratio of α/ai, capacity distortion index may be positive. This shows that the gener-
ation of these GenCos has been decreased with respect to the perfectly competitive
market. It means that these GenCos indirectly can control the market price and make
the GenCos with small α/ai ratio become the marginal producers by capacity
withholding. It may cause the producers with small α/ai to produce more than their
perfectly competitive generations (qpi), and for them qei > qpi (Δqdistortion
i < 0) can be
achieved.
Δθwithheld Δθdistortion
The terms i ai and i ai represent the effect of the capacity limits of GenCos
on capacity withholding. Then, we would have
8
pðeÞ
>
> when qmin qi qmax
<0 i i
pðeÞ
Δθi
withheld
¼ < 0 when qi ¼ qi max ð6:17Þ
>
>
: pðeÞ
> 0 when qi ¼ qmin i
8
>
<0 when qi qpi qmax
min
i
Δθdistortion ¼ < 0 when qpi ¼ qmax ð6:18Þ
i
>
:
i
> 0 when qpi ¼ qmin
i
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 127
pðeÞ
It should be noted that according to (6.8), (6.9), and (6.10), qei qpi qi can be
pðeÞ
obtained. It means that if qi is between qmin i and qmax
i , qei and qpi are also between
qmin
i and qmax
i .
By using (6.2), Δπ distortion can be obtained by
X
N
Yp ¼ qpi ¼ Qp ð6:19cÞ
i¼1
Ideally, with the increase of the total capacity distortion (ΔQdistortion), the
expected outcome would be the increment in Δqdistortion
i . (6.19) represents that
Δqdistortion
i is affected by the aggregate demand function’s slope. It indicates that in
a market with a large slope of demand function, producers will have more potential
ability to increase the market price.
By summing capacity distortions of all GenCos in (6.16), ΔQdistortion can be
obtained by
N
P
α Δθdistortion
ai qei þ i
ai
i¼1
ΔQdistortion ¼ ð6:20Þ
P
N
α
1þ ai
i¼1
N
P
α Δθdistortion
a i qi þ
e i
ai
i¼1
CF ¼ N ð6:22cÞ
P α Δθwithheld
ai qei þ i
ai
i¼1
The superscript uc shows the unconstrained production. It can be seen that in the
case of limited production, this index is composed of two parts. The first part is
unconstrained DWI (DWIuc) which is independent of the market model. The second
part is constraint factor of DWI (CF) which indicates the decrease of the potential
ability of GenCos for capacity withholding due to the capacity limits. If no GenCo is
pðπ e Þ
binding at its generation constraints (i.e., qmin i < qi < qmax
i ), that means
Δθiwithheld
¼ Δθidistortion
¼ 0 and CF ¼ 1, and then DWI ¼ DWIuc can be obtained.
We present the case study of GenCos’ capacity withholding by using a three GenCo
power system. The GenCos’ cost parameters are shown in Table 6.1. The aggregate
demand parameters are α ¼ 2 $/(MW2) and β ¼ 90 $/MW. We have used GAMS
software and CONOPT solver to solve the optimization problems of the case study.
The following simulation cases are carried out:
Case A Perfect competition
Case B Cournot-type competition
Case C Cournot-type competition with decrease in demand elasticity
Case D Cournot-type competition with decrease in demand elasticity and capacity
of GenCos
The simulation results of cases A and B are listed in Table 6.2. It is assumed that
pðeÞ
the generation upper limits of GenCos are greater than qi . Therefore, the produc-
tion constraints can be disregarded in the analyzing of capacity withholding indices.
These indices are listed in Table 6.3. It can be observed that the capacity withholding
of GenCo 1 is higher than the others due to its small slope of the marginal cost.
To analyze the impacts of capacity limits of GenCos and demand elasticity on
capacity withholding, two cases C and D are performed. In these cases, GenCos use
Cournot-type competition and the demand function’s slope is decreased from α ¼ 2
Table 6.2 Simulation results for cases A, B, C1, C2, C3, and C4
pðeÞ
Case GenCo no. qei (GW) qpi (GW) qi (GW) λ ($/MW)
A 1 – 13.375 – 25.375
2 – 10.250 –
3 – 8.687 –
B 1 9.11 – 27.330 39.34
2 8.38 – 19.553
3 7.83 – 15.660
5
4
3
DWI
2
1
0
α=1. α=1. α=1. α=0. α=0. α=0. α=0.
α=2 α=1
7 5 3 8 5 4 3
DWI(Case C) 0.187 0.214 0.235 0.262 0.315 0.365 0.488 0.532 0.606
DWI(Case D) 1 0.816 0.672 0.536 0.339 0.365 0.488 0.532 0.606
CF(Case D) 5.34 3.82 2.85 2.04 1.07 1 1 1 1
elasticity, the expected result in terms of DWI would be a decreasing pattern along
with the releasing capacities of generation until the minimum value of CF (CF ¼ 1)
is reached. The value 1 of CF indicates that the values of DWI in cases C and D are
pðeÞ
equal and for all GenCos (qi < qmax i ).
In the electricity markets, the potential ability of GenCos for capacity withholding
increases during the transmission congestion [3, 4, 15, 18, 25, 28, 29]. In Sect. 6.2, a
method for the capacity withholding analyzing of GenCos in electricity markets
based on some comparison indices has been presented. These indices can be used for
the capacity withholding assessment of GenCos. However, the impacts of transmis-
sion network congestion have been not considered in them, and they cannot be used
for nodal capacity withholding analyzing during the transmission congestion. This
section presents the development of the aforementioned method and offers some
contributions in comparison with the previous section.
1. Market model is extended to evaluate the transmission limits and other parame-
ters that affect GenCos’ capacity withholding.
2. Nodal capacity distortion index, nodal price distortion index, and nodal capacity
withholding index are formulated to assess the effect of transmission network
congestion on GenCos’ capacity withholding. In addition, DWI is reformulated to
include the transmission network limits.
3. In order to evaluate the potential ability of GenCos for exerting capacity with-
holding by considering the transmission congestion, a new index, called nodal
withholding-supply ratio (NWSR), is defined in this section.
The main goal of this subsection is to model the capacity withholding of GenCos
during transmission congestion. This will be done by modeling a game theory-based
formulation. The Cournot model is consider to model the strategic behavior of
generators in the presence of demand elasticity. Then, the model of Cournot-Nash
equilibrium (CNE) is applied to calculate the market equilibrium considering the
GenCos’ capacity withholding. For the sake of simplicity, the market equilibrium
model is first formulated with any constraints; then it is extended for the market with
transmission congestion.
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 131
Suppose that there are I GenCos in the power market, and each of them has a
generation unit. The cost function of GenCo at bus i is determined as follows:
λ j ¼ α j y j þ β j : ð6:24Þ
If there is no transmission congestion, the nodal price at bus n, λn and the price at
the reference bus λ are equal. By ignoring the transmission losses, the generation of
all GenCos and the aggregate demand are the same as described in (6.25)
X X X X
qi ¼ y j ¼ λ 1=α j þ β j =α j ¼ λ=αequal þ βequal =αequal ð6:25Þ
i2I j2J j2J j2J
where αequal and βequal are equivalent parameters of the aggregate demand. Based on
(6.25), the aggregate demand function can be expressed as an inverse demand
function as follows:
X
λ ¼ αequal y j þ βequal : ð6:26Þ
j2J
Each GenCo in the market competes with the others to maximize its own profit by
choosing its generation. The maximization problem of a GenCo at bus i can be
written as
max Ωi ¼ λi qi ci : ð6:27Þ
!
Xh i X X
max 1=2α j y2j þ β jy j 1=2ai qi þ bi qi
2
1=2αequal 2
qi
j2J i2I i2I
ð6:28Þ
X X
s:t: qi y j ¼ 0: ð6:29Þ
i2I j2J
Equation (6.28) shows the decreased social welfare in Cournot model due to the
GenCos’ capacity withholding and includes three terms:
• Demands’ cost
• GenCos’ cost
• Reduction in social welfare due to the capacity withholding of GenCos
The first two terms explain the social welfare in a perfectly competitive market.
By involving the third term, the market is cleared at the CNE point; otherwise, the
market is cleared at the perfect competition point. In fact, the third term explains the
effect of oligopolistic competition on the social welfare of perfectly competitive
market. It is to be considered that in the problem (6.28 and 6.29), the optimization
problems of each GenCo is implicit. By solving this single optimization problem, the
generations of all GenCos are determined [28, 29].
Given the strategic decisions of GenCos, the ISO would determine the
corresponding generation of GenCos, the consumption of loads, and each node’s
price by solving the market clearing problem based on DC power flow. The ISO’s
optimization problem in Cournot model can be written in the vector form as
j k j k j k
max 1=2yT αy þ yT β 1=2qT aq þ qT b 1=2qT αequal q ð6:30Þ
in which
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 133
($/MW)
a/α Diagonal matrix of the vector a=α
pmax =pmin Vector of maximum/minimum generation capacity constraints of
GenCos (MW)
ymax =ymin Vector of maximum/minimum demand capacity constraints of
consumers (MW)
F Vector of power transmission line flow constraints, F ¼
½F 1 , F 2 , . . . , F m , . . . , F M T (MW)
T Power transfer distribution factors (PTDFs) matrix
T TI =T TJ Rows of transposed T matrix for generation/consumption buses
Tmn Element of T matrix
eI =eJ Identity vector with I/J dimension
The key point of Cournot model is to solve the optimization problem of the social
welfare subject to generation and consumption balance, transmission network limits,
and capacity limits of GenCos. The Lagrangian function of the problem (6.30, 6.31,
6.32, 6.33, and 6.34) can be expressed as
h i h i
L q, y, λ, μþ , μ , θmax , θmin , ωmax , ωmin ¼ 1=2yT αy þ yT β 1=2qT aq þ qT b
h i h i h i
1=2qT αequal q þ λ eTI q eTJ y μTþ T q y F þ μT T q y þ F
θT max q qmax þ θT min q qmin ωT max y ymax þ ωT min y ymin :
ð6:35Þ
∂L=∂q ¼ aq b αequal q þ λIq T TI μþ μ θmax θmin ¼ 0
∂L=∂y ¼ αy þ β λIy þ T TJ μþ μ ωmax ωmin ¼ 0
eTI q eTJ y ¼ 0
FT qy F
qmin q qmax
ymin y ymax
h i
μþ T q y F ¼ 0
h i
μ T q y þ F ¼ 0
θmax q qmax ¼ 0
θmin qmin q ¼ 0
ωmax y ymax ¼ 0
ωmin ymin y ¼ 0
μþ , μ , θmax , θmax , ωmax , ωmin 0:
ð6:36Þ
where
μþ =μ Vector of Lagrangian multiplier associated with positive/
negative direction of transmission line flow constraints ($/MW)
μeþ =μe Vector of Lagrangian multiplier in oligopoly market associated
with positive/negative direction of transmission line flow
constraints ($/MW)
μpþ =μp Vector of Lagrangian multiplier in perfectly competitive market
associated with positive/negative direction of transmission line
flow constraints ($/MW)
θmax =θmin Vector of Lagrangian multiplier associated with maximum/
minimum generation capacity constraints of GenCos ($/MW)
θe max =θe min Vector of Lagrangian multiplier in oligopoly market associated
with maximum/minimum generation capacity constraints of
GenCos ($/MW)
θp max =θp min Vector of Lagrangian multiplier in perfectly competitive market
associated with maximum/minimum generation capacity
constraints of GenCos ($/MW)
θpðeÞ max =θpðeÞ min Vector of Lagrangian multiplier in perfectly competitive market
with the nodal prices of oligopoly market associated with
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 135
Nodal capacity withholding index and NWSR can be used for analyzing GenCos’
market power. They compare the generations at oligopoly market and perfectly
competitive market with the nodal prices of oligopoly market. Nodal price distortion
indices and nodal capacity distortion can be used to analyze the distortions of the
oligopoly market from perfectly competitive market.
pðeÞ
Δqwithheld
i ¼ qi qei : ð6:40Þ
Δqdistortion
i ¼ qpi qei : ð6:41Þ
Δλdistortion
n ¼ λen λpn : ð6:42Þ
136 S. Salarkheili and M. SetayeshNazar
NWSRi ¼ Δqwithheld
i =qei : ð6:43Þ
pðeÞ
Based on (6.46), qi can be determined by
Because of GenCos’ capacity withholding, the nodal prices and the generations of
GenCos deviate from their perfectly competitive levels to oligopoly. Then, the nodal
capacity distortion and the nodal price distortion can be formulated as Δqdistortion
i ¼
qpi qei and Δλdistortion
n ¼ λen λpn , respectively. It is to be noted that the withheld
capacity for each GenCo is not equal to qpi qei ; rather, it is formulated as
pðeÞ
Δqwithheld
i ¼ qi qei , indicating that the generation of GenCo at bus i with nodal
pðeÞ
price λei would be qi but, in fact, it was reduced to qei because of the capacity
pðeÞ
withholding of GenCo in the oligopoly market. It should be noted that qi is the
competitive generation, which could be generated at λi by GenCo at bus i. According
e
to (6.46) and (6.47), if GenCo at bus i does not hit its generation capacity constraints,
pðeÞ
qi gets larger than qpi ; then, Δqwithheld
i > Δqdistortion
i can be obtained.
By substituting (6.44) and (6.48) into (6.40), Δqwithheld
i can be determined by
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 137
Δqwithheld
i ¼ 1=ai αequal qei Δθwithheld
i ð6:49Þ
@ P P e min A
ai λ e μeþ
m T mi þ μe
m T mi bi θ i
e max
þ θi
m2M m2M
It can be seen that in addition to αequal/ai ratio, NWSRi could be influenced by μeþ m
,μe
m , Tmi, Δθwithheld
i , θei max , θei min , and bi in a transmission-constrained market.
In the same way, Δqdistortion
i can be obtained by using (6.41, 6.42, 6.44, 6.45, 6.46,
and 6.47) as follows:
Δqdistortion
i ¼ 1=ai αequal qei Δλdistortion
i Δθdistortion
i ð6:52Þ
where Δθdistortion
i ¼ θpi max þ θei min θpi min θei max .
Unlike Δqwithheld
i , which is positive, Δqdistortion
i can be negative or positive. A
positive (negative) value of Δqdistortion
i implies that the generation of GenCo has been
decreased (increased) with respect to the perfect competition model. The sign of
capacity distortion index depends on the slope of the marginal cost function of
GenCo i. For high-cost GenCos, capacity distortion index can be negative. It means
that these GenCos don’t have high market shares; however, because of the capacity
withholding of low-cost GenCos, which have high ability for the capacity
138 S. Salarkheili and M. SetayeshNazar
withholding, their market shares may be increased. Then, the high-cost GenCos in
oligopoly market may generate more than in perfectly competitive market, and
negative Δqdistortion
i may be achieved.
Based on (6.42), (6.45), and (6.47), Δλdistortion
i at the generator connected buses
can be obtained by
X
Δλdistortion
i ¼ Δλdistortion Δμdistortion
m T mi ð6:53Þ
m2M
ð6:54Þ
Another index to assess the capacity withholding may be defined by DWI, which
can be obtained by
X X
DWI ¼ Δqdistortion
i = Δqwithheld
i : ð6:55Þ
i2I i2I
As mentioned before, DWI gives the information about the potential ability of
GenCos in power markets for capacity withholding. Equation (6.56) provides
reformulation of this index by considering network and GenCos limits:
0 1
1
DWI ¼ @1 þ P αequal =a A
i
i2I
0P P P 1
αequal qei =ai Δθdistortion
i =ai Δμdistortion
m T mi =ai
Bi 2 I i2I m 2 M, i 2 I C
@ P equal e P A ¼ DWIUC CF
α qi =ai Δθwithheld
i =a i
i2I i2I
ð6:56Þ
0 1
1
where DWIUC ¼ @1 þ P αequal =a A and
i
i2I
0 P equal e P distortion P 1
α qi =ai Δθi =ai Δμdistortion
m T mi =ai
Bi 2 I i2I m 2 M, i 2 I C
CF ¼ @ P equal e P A.
α qi =ai Δθwithheld
i =a i
i2I i2I
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 139
The constraint factor (CF) shows the effect of the GenCos and transmission
network constraints on the potential ability in the market for capacity withholding.
If the GenCos and the transmission lines do not hit their constraints, therefore
Δθdistortion
i , Δθwithheld
i , and Δμdistortion
m are equal to zero. Then, CF is equal to 1 and
DWI is equal to the unconstrained DWI (DWIUC). According to (6.56), DWIUC
depends on αequal/ai ratio. CF becomes greater than 1, when the GenCos hit their
constraints. Then, DWI becomes greater than DWIUC. It means that the potential
ability in the market for capacity withholding is decreased because of the capacity
limits of GenCos. It can be seen that when the transmission line flow constraints are
exceeded, CF becomes lower than 1. The expected outcome with the increase of
transmission congestion, in terms of the capacity withholding of GenCos, could be
an increasing pattern along with the decreasing DWI.
Table 6.5 Simulation results of modified IEEE 30-bus test system for cases A to F
pðeÞ
Case Bus number of GenCo λpi ($/MW) λei ($/MW) qi (MW) qpi (MW) qei (MW)
A 1 34.559 36.028 200 193.549 145.328
2 34.559 36.028 52.937 48.740 48.088
5 34.559 36.028 20.822 19.644 20.250
8 34.559 36.028 21.151 12.291 17.456
11 34.559 36.028 12.056 10 11.261
13 34.559 36.028 12.056 12 12
B 1 27.299 29.403 125.38 97.323 85.263
2 25.189 27.648 28.994 21.968 26.338
5 38.355 38.416 22.733 22.684 22.108
8 39.210 39.311 35 35 33.702
11 38.590 38.795 17.591 17.181 16.432
13 36.990 37.464 14.930 13.981 13.945
C 1 26.106 28.842 117.886 81.415 69.116
2 24.128 27.255 27.871 20 24.211
5 35.921 36.404 21.123 20.737 20.265
8 37.271 37.795 31.749 28.604 24.1
11 36.690 37.330 14.66 13.381 13.256
13 35.190 36.127 12.254 12 12
D 1 29.742 30.996 146.620 129.896 118.694
2 27.545 29.026 32.931 28.700 31.350
5 41.540 41.558 25.246 25.232 24.895
8 42.141 42.117 35 35 35
11 41.496 41.538 23.078 22.993 22.291
13 39.831 40.044 20.090 19.662 19.405
E 1 30.724 32.067 80.444 71.492 65.125
2 29.242 30.795 37.986 33.548 34.508
5 38.553 38.799 23.039 22.842 22.406
8 39.088 39.243 35 35 33.365
11 38.653 38.869 17.740 17.306 16.570
13 37.529 37.905 15.812 15.059 14.769
F 1 23.422 26.340 151.196 112.289 102.823
2 20.575 24.020 20 20 20
5 38.355 38.416 22.733 22.684 22.108
8 39.490 39.431 35 35 34.300
11 38.654 38.750 17.502 17.308 16.348
13 36.496 36.992 13.984 12.992 13.063
Table 6.6 Capacity withholding indices of modified IEEE 30-bus test system for cases A to F
Δqwithheld
i Δqdistortion
i
Bus number of Δλdistortion
i
Case GenCo ($/MW) (MW) (MW) NWSRi DWI
A 1 1.469 54.672 48.221 0.3761 0.6472
2 1.469 4.849 0.652 0.1008
5 1.469 0.572 0.606 0.0282
8 1.469 3.695 5.165 0.2116
11 1.469 0.795 1.261 0.0705
13 1.469 0.056 0 0.0046
B 1 2.104 40.117 12.060 0.4705 0.2209
2 2.459 2.656 4.370 0.1008
5 0.061 0.625 0.576 0.0282
8 0.101 1.298 1.298 0.0385
11 0.205 1.159 0.749 0.0705
13 0.474 0.985 0.036 0.0705
C 1 2.736 48.77 12.299 0.7056 0.2106
2 3.128 3.66 4.211 0.1511
5 0.483 0.858 0.472 0.0423
8 0.524 7.649 4.504 0.3173
11 0.639 1.404 0.125 0.1059
13 0.936 0.254 0 0.0211
D 1 1.254 27.927 11.202 0.2353 0.3143
2 1.481 1.581 2.650 0.0504
5 0.018 0.352 0.337 0.0141
8 0.025 0 0 0
11 0.041 0.788 0.703 0.0353
13 0.213 0.686 0.257 0.0353
E 1 1.343 15.319 6.366 0.2352 0.3653
2 1.554 3.479 0.960 0.1008
5 0.246 0.633 0.436 0.0282
8 0.155 1.635 1.635 0.0490
11 0.216 1.170 0.736 0.0705
13 0.376 1.043 0.290 0.0705
F 1 2.918 48.373 9.466 0.4705 0.2246
2 3.446 0 0 0
5 0.061 0.625 0.575 0.0282
8 0.058 0.700 0.700 0.0204
11 0.096 1.154 0.961 0.0705
13 0.496 0.922 0.071 0.0705
When evaluating cases A and B, it can be observed that because of the transmis-
sion network congestion, NWSRi for the GenCos at buses 1 and 13 increases to
0.4705 and to 0.0705, respectively. Moreover, the NWSRi of GenCo at bus 8 reduces
to 0.0385. In addition, the values of αequal/ai ratios for GenCos at buses 1, 8, and
13 are 0.4705, 0.0705, and 0.2116, respectively. Therefore, on the one hand,NWSRi
values of GenCos at buses 1 and 13 will become lower than their αequal/ai ratios
pðeÞ
when their constraints are bounded in case A (q1 ¼ qmax 1 ¼ 200 MW and qe13 ¼
q13 ¼ 12 MW). On the other hand, in case B, the generation capacity constraints of
min
pðeÞ
GenCos at buses 1 and 13 are not limited (q1 ¼ 125:38 MW and qe13 ¼
13:945 MW), and based on (6.50) and (6.51), NWSRi of GenCos grow up to the
points where NWSR1 ¼ αequal/a1 ¼ 0.4705 and NWSR13 ¼ αequal/a13 ¼ 0.0705. The
opposite outcomes are seen for the GenCo at bus 8. In case A, the generation
pðeÞ
capacity constraints of the GenCo are not limited (qmin 8 < qe8 , q8 < qmax
8 ) and
NWSR8 ¼ α equal
/a8 ¼ 0.2116. Due to the congestion of transmission network in
case B, the GenCo hits its upper generation capacity constraint
pðeÞ
(q8 ¼ qmax 8 ¼ 35 MW ), and NWSR8 reduces to 0.0385. The constant NWSRi
values of GenCos at buses 2, 5, and 11 in cases A and B show that the constraints
pðeÞ
of qi and qei for these GenCos are not limited. Then, their NWSRi for them are not
affected by the transmission network constraints.
Based on (6.56), unconstrained DWI of the system is 0.5122. In case A, DWI and
CF are 0.6472 and 1.2635, respectively. It means that the potential ability of GenCos
for capacity withholding is decreased because of the generation limits of GenCos
1, 11, and 13. Therefore, in case B, due to the congested lines 6, 8, and 9, CF reduces
to 0.4327 causing DWI to reduce to 0.2209. It implies that the GenCos’ ability for
market power is increased. Note that the reduced CF should be ascribed to the
transmission congestion.
The effect of demand elasticity on the GenCos’ capacity withholding is also
analyzed. Cases C and D are the same as case B except for αequal ¼ 0.05292 and
αequal ¼ 0.01764, respectively. These slopes are achieved from 50% increase to 50%
decrease in the inverse demand functions’ slopes of case B, respectively. As listed in
Table III, the indices show an improvement with the decrease of αequal. In addition,
in case D, DWI increases to 0.3143 and the NWSRi decreases.
Other conditions are studied for analyzing the effect of cost parameters on the
GenCos’ capacity withholding. Cases E and F are the same as case B except for
a1 ¼ 0.15 and b1 ¼ 15, respectively. By reducing of a1 in case E, NWSR1 reduces to
0.2352, showing that the GenCos’ potential ability reduces. However, NWSR8
increases to 0.0490; the growth of DWI to 0.3653 indicates that the GenCos’ ability
for market power decreases. Therefore, GenCo at bus 1 is more effective than GenCo
at bus 8, and a decrement in its capacity withholding can affect the market consid-
erably. The values of NWSR1 for GenCos at buses 1, 5, 8, 11, and 13 in case F are
the same as those in case B. Although, NWSR2 reduces to 0, leading DWI to
increase to 0.2246.
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 143
In hosting capacity analysis, the maximum amount DERs that can be installed in
the power system is determined. In this study, the major concerns are technical
constraints such as adequacy and security limits of system. After this analysis, the
market operator or the ISO can evaluate the outcomes of the technical analysis from
market power and capacity withholding assessment. In other words, the ISO can use
the proposed global and nodal indices to measure the effects of DERs on capacity
withholding. In fact, ISO can use global index to analyze the effect of DERs on
capacity withholding in the whole system. In addition, by nodal indices the impacts
of DERs on capacity withholding at each bus of system can be analyzed. Then, the
ISO can use the optimum strategies to control the market power of GenCos. In
addition, ISO can propose new locations of DERs to decrease capacity withholding
in different buses of system by using nodal indices such as nodal price distortion
index and nodal capacity withholding index. The new locations should be analyzed
from technical constraints of system such as substations and feeder capacities. In
Fig. 6.4, the flowchart of this procedure is shown.
Yes
Change the capacity
NO or location of DERs.
Are nodal
indices
NO
acceptable?
Are Technical
Yes
constraints
exceeded?
End NO
146 S. Salarkheili and M. SetayeshNazar
In this chapter, two categories of market power indices were presented to assess the
capacity withholding of GenCos in the presence of DERs. This indices are effective
and they can be used by the ISO besides the technical assessment of hosting
capacity. Future works are needed to evaluate the defined indices when maintenance
scheduling of generators and strategic transmission system operators (TransCo) are
considered.
References
1. Gan, D., & Bourcier, D. V. (2002). Locational market power screening and congestion
management: Experience and suggestion. IEEE Transactions on Power Systems, 17(1),
180–185.
2. Wang, P., Xiao, Y., & Ding, Y. (2004). Nodal market power assessment in electricity markets.
IEEE Transactions on Power Apparatus and Systems, 19(3), 1373–1379.
3. Bompard, E., Ma, Y. C., Napoli, R., & Jiang, C. W. (2006). Assessing the market power due to
the network constraints in competitive electricity markets. Electric Power Systems Research,
76, 953–961.
4. Bompard, E., Lu, W., & Napoli, R. (2006). Network constraint impacts on the competitive
electricity markets under supply-side strategic bidding. IEEE Transactions on Power Apparatus
and Systems, 21(1), 160–170.
5. Bompard, E., Ma, Y. C., Napoli, R., & Abrate, G. (2007). The demand elasticity impacts on the
strategic bidding behavior of the electricity producers. IEEE Transactions on Power Apparatus
and Systems, 22(1), 188–197.
6. Hesamzadeh, M. R., Biggar, D. R., & Hosseinzadeh, N. (2011). The TC-PSI indicator for
forecasting the potential ability for market power in wholesale electricity markets. Energy
Policy, 39, 5988–5998.
7. Lee, Y. Y., Hur, J., Baldick, R., & Pineda, S. (2011). New indices of market power in
transmission-constrained electricity markets. IEEE Transactions on Power Apparatus and
Systems, 26(2), 681–689.
8. Gabriel, S. A., Conejo, A. J., Fuller, J. D., Hobbs, B. F., & Ruiz, C. (2013). Complementarity
modeling in energy markets. New York: Springer Press.
9. Latorre, M. D. L., & Graville, S. (2003). The Stackelberg equilibrium applied to AC power
systems-a non-interior point algorithm. IEEE Transactions on Power Apparatus and Systems,
18(2), 611–618.
10. Chen, H., Wong, K. P., Chung, C. Y., & Nguyen, D. H. M. (2006). A coevolutionary approach
to analyzing supply function equilibrium model. IEEE Transactions on Power Apparatus and
Systems, 21(3), 1019–1028.
11. Hasan, E., & Galiana, F. D. (2008). Electricity markets cleared by merit order – Part II: Strategic
offer and market power. IEEE Transactions on Power Apparatus and Systems, 23(2), 372–379.
12. Ruiz, C., Conejo, A. J., & Bertand, R. G. (2008). Some analytical results pertaining to Cournot
models for short-term electricity markets. Electric Power Systems Research, 78, 1672–1678.
13. Bushnell, J. B. (2003). Looking for trouble competition policy in the US electricity industry.
Center for the Study of Energy Markets: Univ. California Energy Institute, California.
14. Baldick, R., Grant, R., & Kahn, E. (2004). Theory and application of linear supply function
equilibrium in electricity markets. Journal of Regulatory Economics, 25, 143–167.
6 Impact of Capacity Withholding on Hosting Capacity Analyzing 147
15. Cunningham, L. B., Baldick, R., & Baughman, M. L. (2002). An empirical study of applied
game theory: Transmission constrained Cournot behavior. IEEE Transactions on Power Appa-
ratus and Systems, 17(1), 166–172.
16. Yang, Y., Zhang, Y., Li, F., & Chen, H. (2012). Computing all Nash equilibria of multiplayer
games in electricity markets by solving polynomial equations’. IEEE Transactions on Power
Apparatus and Systems, 27(1), 81–91.
17. Helman, U., & Hobbs, B. F. (2010). Large-scale market power modeling: Analyzing of the
U.S. eastern interconnection and regulatory applications. IEEE Transactions on Power Appa-
ratus and Systems, 25(3), 1434–1448.
18. Yu, C. W., Zhang, S. H., Wang, X., & Chung, T. S. (2010). Modeling and analysis of strategic
forward contracting in transmission constrained power markets. Electric Power Systems
Research, 354–361.
19. Xiaohong, G., Ho, Y. C., & Peypyne, D. L. (2001). Gaming and price spikes in electric power
markets. IEEE Transactions on Power Apparatus and Systems, 16(3), 402–408.
20. Tellidou, A. C., & Bakirtzis, A. G. (2007). Agent-based analysis of capacity withholding and
tacit collusion in electricity markets. IEEE Transactions on Power Apparatus and Systems, 22
(4), 1735–1742.
21. Hongyan, L. I., & Tesfatsion, L. (2009). Capacity withholding in restructured wholesale power
markets: An agent-based test bed study, Power Systems Conference and Exposition.
22. Mohtavipour, S. S., Haghifam, M. R., & Sheik-El-Eslami, M. K. (2011). Emergence of capacity
withholding: An agent based simulation of a double price cap electricity market. IET Genera-
tion Transmission and Distribution, 6(1), 69–78.
23. Chattopadhay, D. (2004). A game theoretic model for strategic maintenance and dispatch
decision. IEEE Transactions on Power Apparatus and Systems, 19(4), 2014–2021.
24. Wang, J., Shahidehpour, M., Li, Z., & Botterud, A. (2009). Strategic generation capacity
expansion planning with incomplete information. IEEE Transactions on Power Apparatus
and Systems, 24(2), 1002–1010.
25. Hesamzadeh, M. R., Hosseinzadeh, N., & Wolfs, P. J. (2010). Transmission system augmen-
tation based on the concepts of quantity withholding and monopoly rent for reducing market
power. IEEE Transactions on Power Apparatus and Systems, 25(1), 167–180.
26. Salarkheili, S., & Akbariforoud, A. (2012). Market power assessment in electricity markets:
Supply function equilibrium-based model. International Transactions on Electrical Energy
Systems, 23(4), 553–569.
27. Salarkheili, S., & SetayeshNazar. M. (2015). New indices of capacity withholding in power
markets. International Transaction on Electrical Energ Systems, 25(1), 180–196.
28. Salarkheili, S., & SetayeshNazar, M. (2016). Capacity withholding analysis in transmission-
constrained electricity markets. IET Generation Transmission and Distribution, 10(2),
487–495.
29. Salarkheili, S., & SetayeshNazar, M. (2017). Capacity withholding assessment in the presence
of integrated generation and transmission maintenance scheduling. IET Generation Transmis-
sion and Distribution, 11(16), 3903–3911.
30. Ameri, M., Rahimian, M., & Latifi, M. A. (2018). Capacity withholding constrained by
operational limits of generation under financial virtual divestiture in a day-ahead market.
IEEE Transactions on Power Apparatus and Systems, 33(1), 771–780.
31. Guo, H., Chen, Q., Xia, Q., & Kang, C. (2019). Market power mitigation clearing mechanism
based on constrained bidding capacities. IEEE Transactions on Power Systems, 34(6),
4817–4827.
Chapter 7
A Generalised Deterministic Approach
to Evaluate PV Hosting Capacity of LV
Distribution Networks Under Different
Operating Conditions
7.1 Introduction
D. Chathurangi (*)
Department of Electrical Engineering, University of Moratuwa, Moratuwa, Sri Lanka
Collaborated Research, School of Electrical, Computer & Telecommunications Engineering,
Faculty of Engineering and Information Science, University of Wollongong,
Wollongong, Australia
e-mail: dmcwl143@uowmail.edu.au
U. Jayatunga
Department of Electrical Engineering, University of Moratuwa, Moratuwa, Sri Lanka
e-mail: upuli@uom.lk
S. Perera · A. Agalgaonkar
School of Electrical, Computer & Telecommunications Engineering, Faculty of Engineering
and Information Science, University of Wollongong, Wollongong, Australia
e-mail: sarath@uow.edu.au; ashish@uow.edu.au
accommodate DG. Bidirectional power flows are in reality with increased solar PV
penetration levels and can degrade the performance of LV distribution networks
resulting in a number of power quality and other related issues (such as overvoltage
limit violations, voltage unbalance, overloading of transformers and feeders, exces-
sive line losses, high harmonic distortion levels exceeding the stipulated limits and
protection mal-operation [2]). Degree of the occurrence of these problems increases
with the maximum permissible penetration level, especially when the hosting capac-
ity exceeds acceptable limits [3, 4].
Key negative impacts associated with high solar PV penetration levels may be
overcome by network reinforcements such as asset upgrades (conductors and trans-
formers) and advanced monitoring and protection schemes. However, such rein-
forcements are not economically feasible in general. Given the potential obstacles, it
is desirable to develop systematic approaches for new solar PV deployment from
network planning perspectives so as to optimally utilise existing assets instead of the
network augmentation.
Distribution network service providers (DNSPs) are required to ensure that the
connection of new solar PV systems does not violate any associated technical limits.
Technical issues that limit solar PV hosting capacity of LV distribution networks
include thermal overloading of the network components, overvoltage, voltage reg-
ulation, increasing short circuit levels and power quality issues, in addition to
islanding considerations and the possibility of reverse of power flows. At present,
most DNSPs follow simple practices and rules of thumb in order to assess the solar
PV hosting capacity of existing distribution networks. However, there is a necessity
for the development of systematic approaches for hosting capacity evaluation
considering influential factors such as geographical layout of networks as well as
electrical parameters.
Solar PV hosting capacity has to be determined using a transparent approach, that
is, to define an appropriate index or indices subject to well-defined performance
criteria or limits (e.g. feeder voltage, line loading level, voltage unbalance, harmonic
levels) [5]. However, such a hosting capacity can take multiple values when different
performance indices are utilised as constraints. Moreover, a series of constraining
performance indices can be utilised in order to assess the solar PV hosting capacities
of a given distribution network by accommodating stochastic solar PV generation at
different locations through the use of Monte Carlo methodologies. Accordingly, a
single value cannot be specified as the PV hosting capacity for a given network.
The work presented in this chapter covers a generalised deterministic approach to
evaluate the solar PV hosting capacity of low-voltage (LV) distribution networks
considering overvoltage as the constraining criterion and the performance parameter.
This drive to develop a deterministic approach is justified by considering the
drawbacks of a stochastic approach. Furthermore, the generalised framework of
the solar PV hosting capacity estimation aims at capturing all network and PV
plant characteristics to maximise the solar PV penetration level in LV distribution
networks.
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 151
This section discusses the most common impacts of higher solar PV penetration
levels in distribution networks, based on a case study of an urban LV distribution
network in Sri Lanka which is reported in this section with the aid of feeders with the
highest levels of rooftop solar PV generation. This case study investigates overvolt-
age issues associated with high solar penetration levels in LV distribution networks.
The selected test network consists of a 250 kVA, 11 kV/400 V distribution
transformer with three radial feeders that supply a total of 336 customers including
253 single-phase customers and 83 three-phase customers. More than 90% of the
load on the system is shared by 315 residential customers, and the remaining is
classified as commercial consumers. The installed solar PV capacity in the system
associated with 24 solar PV systems (single-phase and three-phase) distributed along
the three feeders is 104 kW. The geographical layout of the network is shown in
Fig. 7.1, while Table 7.1 provides relevant network details. Measured voltage pro-
files of selected six customers in feeders 2 and 3 at different locations (marked as A,
B, C, D, E and F in Fig. 7.1) are analysed in Sect. 7.2.1 in relation to the overvoltage
problem.
Load profiles of more than 30% of the customers are derived based on the actual
15-minute energy data recorded remotely through smart meters. The remaining load
profiles are obtained by deriving an average load profile (using the measured data).
These average load profiles are used to generate load profiles of the remaining
customers (without smart meters) based on the monthly energy consumption data
acquired for individual customers. Further, the constant power factor of 0.8 lagging
is assumed for all the customers.
Fig. 7.2 Voltage profiles of solar/non-solar customers at different locations; (a) feeder 2 and
(b) feeder 3 [6]
154 D. Chathurangi et al.
In this section, a stochastic analysis framework followed using Monte Carlo simu-
lation (MCS) method to assess the solar PV hosting capacity of the test system is
presented. The performance indices are selected as overvoltage and overloading
criterion with regard to the selection of an appropriate solar PV hosting capacity
level. Accurate analysis of solar PV hosting capacity requires the development of
solar deployment scenarios for all possible future PV locations and corresponding
PV installation size in the network. It is not feasible to simulate all possible PV
deployment scenarios for a given network. Thus, to represent the impact of the
location of a new solar PV connection, the MCS approach is used with a limited
number of PV deployment scenarios at a particular level of penetration.
For the test network, solar PV penetration level (N ) is defined as the ratio of the
number of customers with PV units and the total number of customers. With
increasing penetration levels, a number of stochastic solar PV deployment scenarios
(M) can be considered in such a way that in each scenario, the location of the PV
customers is randomly selected from the pool of solar customers in the network.
Furthermore, each simulated scenario is unique in order that PVs are deployed. The
overall framework developed for solar deployment scenarios for each of the pene-
tration level is illustrated in Fig. 7.3.
Since this systematic approach is applied to a practical distribution network, solar
PV connections are only allocated to customers with energy consumption greater
than 120 kWh per month, assuming that the solar PV installations are financially
feasible for these customers. Capacity of the solar PV system is determined based on
the monthly energy consumption and assigned as a net metered1 customer. The total
customer demand is assumed to be constant for different solar PV penetration levels.
Following this method, several deployment scenarios are systematically developed,
until all financially feasible solar deployments are connected to the network. It is
assumed that all installed solar PV units are to provide active power at unity power
factor which is a common practice.
The influential factors governing the impact of solar PVs on the performance of
the LV distribution system include an accurate representation of the solar PV
systems, in particular the PV unit size and solar irradiance level which inherently
drives the solar PV output. Solar irradiance profiles for each PV units were devel-
oped using the System Advisor Model (SAM) software considering the irradiance
data applicable to Sri Lanka.
Maximum feeder voltage violation and line overloading are commonly reported as
the deciding and critical factors in assessing solar PV hosting capacity in distribution
networks. The considered test network is already associated with overvoltage vio-
lation issues as discussed in Sect. 7.2.1. With this in mind, solar PV hosting capacity
was evaluated for the test network by defining the maximum feeder voltage in
addition to considering the line and transformer loading violations. In each PV
deployment scenario, network power flows are monitored, and the corresponding
network voltages and line and transformer loading levels are calculated in order to
verify whether any performance standards are violated. The threshold values for
each performance index considered in this analysis are given in Table 7.2. It is
important to highlight that threshold value for a given performance index is mutually
exclusive of the rest of performance indices.
MCS results shown in Fig. 7.4a–c elaborate the resulting maximum feeder
voltage, maximum line loading and transformer loading levels, respectively, under
different solar penetration levels for a total number of 66,000 solar deployment
scenarios. Based on the MCS results, two levels of solar PV hosting capacity can be
identified as minimum hosting capacity (HCmin) and maximum hosting capacity
(HCmax) compliant with the given performance index, as shown in Fig. 7.4 and
summarised in Table 7.3. With the minimum hosting capacity level, none of the solar
penetration levels violates the criteria of the performance index, i.e. 48 kW of solar
1
This allows customers to be paid in cash for any surplus power generation from the solar PV
systems at the end of monthly billing cycles.
156 D. Chathurangi et al.
Fig. 7.4 Hosting capacity limits for performance index: (a) maximum feeder voltage criterion,
(b) maximum line loading criterion and (c) transformer loading criterion
PVs can be integrated to the test network not to violate the steady-state maximum
feeder voltage limits without any concern of the PV location. Similarly, considering
only the thermal overloading limits of line conductors and transformer, the same
network can withstand up to 244 kW solar PV units (note that the voltage violation
criterion is excluded in thermal loading analysis and vice versa). Solar penetration
levels in excess of the maximum hosting capacity limit, independent of the location
of the solar installation, will violate the relevant operational limits. Moreover, the
range in between minimum and maximum hosting capacity levels, certain solar
penetration levels which arise from random PV locations, may violate the network
constraints. Therefore, detailed studies are necessary with precise solar PV locations
to verify that a given level of penetration is safe.
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 157
Fig. 7.5 Variation of hosting capacity limit with the number of solar deployment scenarios for
overvoltage criterion; minimum hosting capacity and maximum hosting capacity
Solar PV location and its aggregated capacity along the feeder are vital factors in
assessing hosting capacity for a given distribution network. A simple case study is
presented here to show the influence of solar PV system location and aggregated
capacity on maximum penetration level. A single feeder network which is 1 km long
(conductor type: AAC-Fly type (all aluminium conductor), R ¼ 0.4505 Ω/km,
X ¼ 0.292 Ω/km) is considered with the daily peak demand of 40 kW and
30 kVAr. The total load is assumed to be divided equally between three phases
and evenly distributed along the feeder. Figure 7.6 shows the voltage profile of the
feeder for maximum solar power that can be injected at three different locations; an
aggregated capacity is considered at a given location at a time (distance are measured
from the transformer end), not to violate overvoltage conditions. The results clearly
show how power injection levels tend to reduce when solar PV is connected towards
the feeder end.
How solar PV location and aggregated capacity affect the PV hosting capacity levels
is assessed by considering a simple LV distribution network which has two feeders.
The study network comprises 11 kV/400 V, 100kVA transformer with a loading
level of 50%. Two different feeder lengths, 500 m and 1000 m (conductor type:
AAC-Fly type), are considered to analyse the effect of feeder length on the PV
hosting capacity. To assess the effect of load changes on PV hosting capacity, two
loading levels are considered—Case 1: 60% of total load is allocated to feeder 1, and
the rest (40%) is assigned to feeder 2, and Case 2: 40% of total load is allocated to
feeder 1, and the remaining (60%) is assigned to feeder 2. In addition, constant
power factor, 0.8 lagging, is assumed for all load points. For simplicity, balanced and
evenly distributed loads are considered. Each feeder length is analysed by dividing
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 159
into three equal-length segments as shown in Fig. 7.7 where the randomness of the
solar PV deployment is applied to each segment.
Accordingly, four solar deployment scenarios are considered by changing the
locations of the solar PV connections in relation to Cases 1 and 2:
Scenario 1: Solar PVs are distributed in the feeder end segment.
Scenario 2: Solar PVs are distributed in the feeder middle segment.
Scenario 3: Solar PVs are distributed in the feeder front segment.
Scenario 4: Solar PVs are distributed over all segments of the feeder.
160 D. Chathurangi et al.
Initial values of maximum feeder voltage drop and line loading levels of the network
are given in Table 7.4 for both Case 1 and Case 2 (without any solar PV
connections).
The test network is analysed to examine the maximum connectable solar PV
capacity at different segments ensuring that the voltage upper limits are not violated.
Figure 7.8 shows the maximum connectable solar PV capacities at the boundaries
(the specific case of fixed locations of PV) of the three segments for both Case 1 and
Case 2 (hosting capacity values are given in per unit using a base Sb ¼ 100 kVA
which is the step down transformer rating, thus allowing the line loading level to be
compared against). Further, solar PV units are to provide only active power operat-
ing at unity power factor. The voltage at the secondary of the transformer was
maintained to a constant value of 1 pu.
The same two-feeder network is analysed using MCS method, and for illustrative
purposes, results obtained for Case 1 scenarios are shown in Figs. 7.9 and 7.10 with
Fig. 7.8 Maximum connectable solar capacity at different locations of two-feeder network [6]
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 161
Fig. 7.9 Variation of segment-wise hosting capacity levels for overvoltage criteria; Case 1, solar
PVs distributed in (a) feeder end segment, (b) feeder middle segment, (c) feeder front segment and
(d) over all segments of the feeder [6]
regard to overvoltage and line loading criterion, respectively. The limit for overvolt-
age exceedance is 1.06 pu, and the acceptable line loading limit is considered to be
100% of the thermal loading of the conductor. In this analysis, transformer loading
level was not considered as the aim was to consider the feeder level PV hosting
capacity only. However, it is vital that transformer loading levels are considered in
the application of the above methodology.
Referring to Fig. 7.10, it can be seen that when solar PVs are distributed in the
feeder front segment, the connectable maximum solar PV capacity level is essen-
tially limited by the line conductor and transformer overloading. Table 7.5 gives the
summary for minimum and maximum hosting capacity (HC) values (considering
overvoltage and line overloading criterion) derived using MCS for both Case 1 and
Case 2.
Referring to simulation outcomes presented with regard to fixed locations, as in
Fig. 7.8, it clearly indicates that the maximum connectable PV capacity increases as
the solar PVs are distributed towards the feeder front (close to the distribution
transformer). As shown in Fig. 7.8, in Case 1, for the end segment, the lowest
hosting capacity was obtained at the feeder end as 0.36 pu (at point A2 in feeder 2).
Further, for the same case, the summation of the hosting capacities at the nearest end
to the transformer of both feeders in the end segment is 1.42 pu (hosting capacities at
points B1 and B2). These values are equal to minimum and maximum solar PV
162 D. Chathurangi et al.
Fig. 7.10 Variation of segment-wise hosting capacity levels for maximum line loading criteria;
Case 1, solar PV distributed in (a) feeder end segment, (b) feeder middle segment, (c) feeder front
segment and (d) over all segments of the feeder
Table 7.5 Solar PV deployment limits for the two-feeder network [6]
Case 1 Case 2
HCmin HCmax HCmin HCmax
Scenario 1: End segment 0.4 1.4 0.45 1.4
Scenario 2: Middle segment 0.5 2.3 0.6 2.3
Scenario 3: Front segment 0.8 3.4 0.9 3.4
Scenario 4: Over all segments 0.7 2.5 0.9 3.0
hosting capacity values obtained from MCS method. Same argument is valid for all
three segments. Hence, for a multi-feeder network, minimum hosting capacity of the
entire network is the lowest of the minimum hosting capacities of all feeders. Thus,
considering a single feeder network, minimum hosting capacity can be defined as the
maximum connectable solar PV capacity at the far end of the given feeder, while
maximum hosting capacity is the maximum connectable solar PV capacity at the
nearest end to the transformer. The same argument is valid for segment-wise analysis
as well. In other words, the minimum and maximum hosting capacities are the
connectable maximum solar PV capacity at the boundaries of each segment, far
end of the feeder segment and nearest end to the transformer, respectively.
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 163
To estimate the effect of number of parallel feeders on solar PV hosting capacity, the
analysis is carried out by increasing the number of parallel feeders. Total loading
level of the network (in terms of transformer loading level) is maintained to a
constant value (50% of the transformer capacity), and identical feeders are used to
observe the impact only of number of parallel feeders.
Analysis described in Sect. 7.3.2 is repeated for different conductor types, AAC-Fly
and ABC, in order to examine the effect of R/X ratio on solar PV hosting capacity.
The relevant conductor specifications are given in column one of Table 7.6.
164 D. Chathurangi et al.
Table 7.6 Solar deployment limits for two-feeder network per unit [6]
Number
of parallel Network 1 Network 2
Type of conductor feeders (feeder length – 600 m) kW (feeder length – 1200 m) kW
AAC-Fly 1 73 55
R ¼ 0.4505 Ω/km 2 55 36
X ¼ 0.292 Ω/km 3 50 30
R/X ¼ 1.5
4 46 28
ABC–70 mm2 1 64 45
R ¼ 0.441 Ω/km 2 51 32
X ¼ 0.08 Ω/km 3 47 27
R/X ¼ 5.5
4 45 25
ABC–50 mm2 1 51 38
R ¼ 0.641 Ω/km 2 38 26
X ¼ 0.08 Ω/km 3 34 21
R/X ¼ 8
4 32 19
Feeder loading level is one of the vital factors which govern the solar PV hosting
capacity. In the analysis of impact of multi-feeders on hosting capacity, each feeder
is applied with different loading levels while maintaining the same total transformer
loading level of 50% of the transformer rating. The corresponding results are
analysed to investigate the impact of loading level on the hosting capacity.
where HCFn, end is the maximum connectable solar PV capacity at the feeder end of
nth feeder (n ¼ 1, 2, 3,. . .. . .). Furthermore, maximum hosting capacity increases
with distribution of solar PV systems closer to the transformer. Comparison of
hosting capacity values obtained from spot simulations (hosting capacities at bound-
aries of the feeder segments) and MCS method has led to define segment-wise
minimum and maximum solar PV hosting capacities as:
• Minimum hosting capacity in a given segment is the connectable maximum PV
capacity at the furthest end (from the transformer) of that segment.
• Maximum hosting capacity in a given segment is the connectable maximum PV
capacity at the nearest end (to the transformer) of that segment.
It is evident that extensive analytical efforts are required for accurate hosting
capacity evaluation process because of the complexity associated with distribution
networks. Thus, there should be generalised approaches that are extendable to
address various aspects of concern. Thus, feeder-based hosting capacity approach
is the best as individual hosting capacity of each feeder has to be established in multi-
feeder networks. Such a methodology needs to be extended in order to better capture
the complexity of network modelling, constraints and technologies that will enable
the greatest potential in generalisation.
In essence, the development of a new deterministic approach for the evaluation of
solar PV hosting capacity in a feeder level is important from a distribution planning
perspective. Further, such an approach will facilitate the connection of solar PVs in a
cost-effective way, thus increasing the associated environmental and social benefits.
In this regard, Sect. 7.4 presents a generalised deterministic approach that can be
used to evaluate the feeder level solar PV hosting capacity levels. By implementing
such a methodology, distribution network operators can allow solar PV connections
up to the feeder minimum hosting capacity level. Further, detailed network analysis
will allow examination of whether new connections are permissible by increasing
PV capacities beyond the bounds of minimum PV hosting capacity level.
166 D. Chathurangi et al.
This section explores a novel deterministic approach which facilitates the evaluation
of maximum solar PV penetration level governed by overvoltage limits which will
enable better understanding of the associated benefits.
This section first explains the principles underlying power flow and voltage calcu-
lations using a simple feeder model with balanced and uniformly distributed loads.
The principles are then extended to deal with distribution feeders embedded with
solar PV generation using examples.
Consider the single-line diagram of a single feeder network shown in Fig. 7.11.
The resistance and reactance of the line are R and X in ohm per unit length,
respectively, and the shunt capacitance is neglected. Vs is the supply voltage at the
secondary of the transformer.
The total voltage drop VDl in per unit (pu) of the feeder can be written as [8]
where Vb is the nominal line-line voltage, Ps and Qs are the total real and reactive
power components of the total load that is uniformly distributed in the feeder, and l is
the total length of the feeder.
Active power and reactive power flow at a distance d from the transformer can be
written as (7.3).
d d
Pd ¼ PS 1 , Qd ¼ QS 1 ð7:3Þ
l l
where Pd and Qd are the active and reactive power flow at distance d from the
transformer.
Ps Qs
Pd Qd
0 d l Distance 0 d l Distance
Fig. 7.12 Active and reactive power profiles along the feeder (without solar PV)
nd
II Quadrant st Solar
Solar I Quadrant Photovoltaic
Photovoltaic
Distribution Distribution
Grid Vs POC Grid Vs R X POC
R X P Q
P Q PV , PV
PV , PV
P P
PL, QL Q PL, QL
Q
Solar Solar
Photovoltaic Photovoltaic
Distribution Distribution
Grid Vs R X POC Grid Vs POC
R X
P Q P Q
PV , PV PV , PV
P P
Q P QL Q PL, QL
L,
rd th
III Quadrant IV Quadrant
(this may include the battery). These operating scenarios of the distribution feeder
can be listed as follows:
• Ist Quadrant – POC operating at a lagging power factor while net P (PL PPV)
and net Q (QL QPV) are positive
• IInd Quadrant – POC operating at a lagging power factor while net P (PL PPV)
is negative and net Q (QL QPV) is positive
• IIIrd Quadrant – POC operating at a leading power factor while net P (PL PPV)
and net Q (QL QPV) are negative
• IVth Quadrant – POC operating at a leading power factor while net P (PL PPV)
is positive and net Q (QL QPV) is negative
In general, deterministic models should be developed in such a way that they can
represent the four-quadrant operation of the solar PV system (which may include a
battery) where possible. In specific terms, the operation at a unity power factor or at a
fixed leading or lagging power factor must be considered.
Reverse power flow due to the solar system causes a voltage rise at the POC.
Using (7.2), the new magnitude of the drop in the voltage from 0 to d0 can be written
as in (7.5):
Fig. 7.14 Active and reactive power profiles along the feeder (with unity power factor solar
inverter)
ðd d 0 Þ RðPPV Pd Þ X ð Q0 þ Qd Þ ðd d 0 Þ
VRd0 d ¼ þ j
2 V 2b 2
X ðPPV Pd Þ þ Rð Q0 þ Qd Þ
ð7:6Þ
V 2b
0
where Q0 ¼ Qs 1 dl
The total voltage rise from distribution transformer to the POC can be obtained by
subtracting (7.5) from (7.6).
dRPPV re dXPPV im
VR0d ¼ 2λ λ 2
VD l þ j 2λ λ 2
VD l ð7:7Þ
V 2b V 2b
where λ ¼ dl , VDre im
l is the real component of voltage drop and VDl is the reactive
component of voltage drop.
The voltage rise in (7.7) is a linear function of the solar PV system as well as its
location on the feeder. Hence, the maximum connectable solar capacity at a fixed
location can be derived theoretically that can give rise to the maximum allowable
voltage rise in the feeder. For unity power factor operation of solar PV inverter, the
maximum connectable solar PV capacity in three phases at a distance d from the
distribution transformer can be formulated as given in (7.8). Here, VPV, VS and VDl
are in pu, and the phase angle deviations of VPV and VS have been neglected.
V 2b
PPV ¼ ðV PV V S Þ þ 2λ λ2 VDl ð7:8Þ
Rd
170 D. Chathurangi et al.
The total voltage rise from distribution transformer to the POC can be
decomposed using active and reactive power flows.
dPPV R þ X tan cos 1 pf pv
VR0d ¼ 2λ λ2 VDre
l
V 2b
! ð7:10Þ
dPPV X R tan cos 1 pf pv
þj 2λ λ VDim
2
l
V 2b
If the PV inverter operates at a fixed power factor, then the reactive power
becomes QPV ¼ PPV tan (cos1pfpv).
Maximum amount of solar power generation that will cause the voltage to reach
the overvoltage limit at a distance d from the transformer for a leading power factor
operation can be formulated as in (7.11). Here, VPV, VS and VDl are in pu, and the
phase angle deviations of VPV and VS have been neglected.
V 2b
PPV ¼ ðV PV V S Þ þ 2λ λ2 VDl ð7:11Þ
d R þ X tan cos pf pv
1
Fig. 7.15 Active and reactive power profiles along the feeder (with leading power factor solar
inverter)
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 171
Similar to Sect. 7.4.2.2, maximum amount of solar power generation that will
cause the voltage to reach the overvoltage limit at a distance d from the transformer
for lagging power factor operation can be formulated as in (7.13). Here, VPV, VS and
VDl are in pu, and the phase angle deviations of VPV and VS have been neglected.
V 2b
PPV ¼ ðV PV V S Þ þ 2λ λ2 VDl ð7:13Þ
d R X tan cos 1 pf pv
Table 7.7 summarises the mathematical models developed for evaluation of solar
PV hosting capacity at a given point in a feeder considering three different operating
conditions of a PV inverter. The simplified equation for the feeder end hosting
capacity (i.e. when d ¼ l) is also shown in the same table.
Fig. 7.16 Active and reactive power profiles along the feeder (with lagging power factor solar
inverter)
172
Table 7.7 Generalised mathematical models for different operation conditions of PV inverters
PV inverter operation mode At a distance d from the transformer At the feeder end (d ¼ l)
Unity power factor V 2b 2
V 2b
Rd ðV PV V S Þ þ 2λ λ VDl Rd fðV PV V S Þ þ VDl g
Leading power factor V 2b V 2b
d ðRþX tan ð cos 1 pf pv ÞÞ
ðV PV V S Þ þ 2λ λ2 VDl RþX tan ð cos 1 pf pv Þ
fðV PV V S Þ þ VDl g
Lagging power factor V 2b V 2b
d ðRX tan ð cos 1 pf pv ÞÞ
ðV PV V S Þ þ 2λ λ2 VDl RX tan ð cos 1 pf pv Þ
fðV PV V S Þ þ VDl g
D. Chathurangi et al.
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 173
Table 7.8 Safe limits of solar Safe limit for HC Safe limit for HC
PV hosting capacity for the calculated from model (kW) simulated (kW)
urban LV network
Feeder 1 32 30
Feeder 2 48 41
Feeder 3 37 38
The proposed deterministic approach in Sect. 7.4.2 is used to evaluate the hosting
capacity of the practical distribution network discussed in Sect. 7.2 (named as Case
1) and of the simple distribution feeder (Case 2) with the aid of DIgSILENT
PowerFactory simulations for verification purposes.
Case 1: Practical LV distribution network
The total midday peak demand of the distribution network is (PS and QS) 82 kW
and 60 kVAr, respectively (40% by the transformer capacity). Based on the assump-
tion that all three feeders are loaded equally, the total voltage drop along each feeder
can be calculated using (7.2). These calculated total voltage drops of feeder
1 (ABC-50 mm2), feeder 2 (ABC-70 mm2) and feeder 3 (ABC-70 mm2) are about
3%, 1.7% and 2.6%, respectively, at the solar peak time. The safe limits of the solar
PV hosting capacities of each feeder are calculated using (7.8) where PV inverters
are operating at the unity power factor. Table 7.8 gives safe limits of hosting capacity
values calculated using the mathematical models and spot simulations constrained
by overvoltage for each feeder. As an example, detailed calculations relevant to
identifying safe limit of solar PV hosting capacity on feeder 1 are given below:
Example The total voltage drop in feeder 1 is approximately 3%. The secondary
line-line voltage at the distribution transformer is set at 1.04 pu to ensure that the
voltage is within the stipulated limits under maximum demand conditions at the
night peak. At this loading level (40%), supply voltage2 of the distribution trans-
former is about 1.024 pu.
From (7.8), the maximum connectable solar PV capacity at feeder end is
4002
PPV,Feeder 1 ¼ fð1:06 1:024Þ þ 0:03g 32 kW
0:641 0:51
2
Secondary voltage of the transformer depends on the transformer loading level.
174 D. Chathurangi et al.
Table 7.9 Maximum connectable solar capacity when solar PV system operates at unity power
factor [9]
HC at feeder end HC at feeder middle
Simulated Calculated Simulated Calculated
AAC-Fly 0.58 0.57 0.96 0.96
ABC – 70 mm2 0.49 0.49 0.84 0.84
ABC – 50 mm2 0.41 0.41 0.69 0.69
Table 7.10 Maximum connectable solar capacity when solar PV system operates at 0.9 leading
power factor [9]
HC at feeder end HC at feeder middle
Simulated Calculated Simulated Calculated
AAC-Fly 0.48 0.48 0.79 0.81
ABC – 70 mm2 0.50 0.50 0.84 0.85
ABC – 50 mm2 0.43 0.43 0.71 0.73
Table 7.11 Maximum connectable solar capacity when solar PV system operates at 0.9 lagging
power factor [9]
HC at feeder end HC at feeder middle
Simulated Calculated Simulated Calculated
AAC-Fly 1.41 0.92 1.98 1.55
ABC – 70 mm2 0.63 0.59 1.06 1.02
ABC – 50 mm2 0.52 0.49 0.85 0.82
The test network is analysed for maximum connectable solar PV capacity at two
different locations so that the upper voltage limit stipulated by the utility (taken as
+6%) is not violated.
The analysis has been repeated for different types of conductors considered in
Sect. 7.3 in order to investigate the effect of R/X ratio on solar PV hosting capacity.
Accordingly, three cases are presented for different operating scenarios of the PV
inverters: unity power factor, 0.9 leading power factor and 0.9 lagging power factor.
Solar PV hosting capacity values which were obtained from the simulations and
mathematical model are given in Tables 7.9, 7.10 and 7.11 for the three scenarios. In
7 A Generalised Deterministic Approach to Evaluate PV Hosting Capacity of LV. . . 175
each case, hosting capacity values relevant to feeder end and middle of the feeder
were evaluated and compared against the results from the mathematical models. All
solar PV capacity values are given in per unit with transformer ratings as base values.
As shown in Tables 7.9, 7.10 and 7.11, solar PV hosting capacities obtained from
simulations and proposed deterministic method are in close agreement verifying the
accuracy of the deterministic approach.
The minimum hosting capacity of a given radial distribution feeder is the
maximum connectable solar PV capacity at the feeder end which can be defined as
the safe limit of the hosting capacity. Thus, the proposed mathematical model
enables the evaluation of a safe limit of the PV hosting capacity of distribution
networks without resorting to complex computer simulations. From a distribution
system planning perspective, the initial use of a deterministic approach is convenient
and more practical than using extensive simulations. Therefore, utilities can make
use of the proposed approach as a guide to evaluate the maximum connectable solar
capacity of a given distribution feeder using the data which are not overly
demanding.
As stated in Sect. 7.3.3, the minimum PV hosting capacity (safe limit of the hosting
capacity) for a given radial distribution feeder is the maximum connectable solar PV
capacity at the feeder end (verified using stochastic Monte Carlo simulations). Thus,
the analysis presented in Sect. 7.4.3 (Case 2) justified that the PV hosting capacity
values obtained from both stochastic and deterministic approaches are in close
agreement. The pros and cons of stochastic and deterministic approaches in PV
hosting capacity analysis are shown in Table 7.12.
Table 7.12 Comparison of stochastic and feeder level deterministic methods in hosting capacity
analysis
Stochastic approach Deterministic approach
Pros Randomness of solar PV location and capacity of Feeder-based limitations and
PV system can be addressed easily diversity can be accommodated
Can be generalised
Improved accuracy
Easier and more practical
Whole network PV hosting
capacity can be maximised
Reduces computational effort
Cons Needs detailed network modelling Can be applied to one location at a
Increased computational burden time
Cannot be generalised Not able to capture all perfor-
Variability in accuracy depends on the number of mance criteria for hosing capacity
simulations and complexity of the network model
176 D. Chathurangi et al.
As per Table 7.12, the deterministic approach for feeder level PV hosting capacity
offers substantial benefits over a stochastic approach.
In this chapter, a feeder based, deterministic approach was proposed for the evalu-
ation of the solar PV hosting capacity in LV distribution networks. The hosting
capacity was essentially constrained by overvoltage limits of many of the LV
distribution networks. The deterministic method was developed in order to evaluate
the maximum allowable solar PV capacity at a given point of distribution feeder
constrained by overvoltage limits. Related to PV hosting capacity evaluation, Monte
Carlo-type stochastic simulations are not practical in distribution network planning
environments due to extensive network modelling efforts and the complexity. Thus,
individual feeder-based analysis approach is recommended and considered to be
more effective.
Based on the findings, a safe limit of the solar PV hosting capacity (the minimum
hosting capacity) is defined as the maximum allowable solar PV capacity at the far
end of a feeder that can be evaluated using the proposed deterministic approach.
Under this framework, distribution network operators and planners can investigate
the capability and limitations of solar PV penetration.
The proposed feeder-based hosting capacity evaluation approach can be seen to
capture all influencing factors on solar PV hosting capacity. Thus, this can be used as
an approximate guide or a rule of thumb to evaluate solar PV hosting capacity at a
given point of LV feeders for overvoltage curtailment without using complex
stochastic techniques.
Further research is needed to confirm the applicability of the developed approach
in realistic network environments. Other network performance indices such as
voltage unbalance due to PV hosting have not been considered in the present work
and require investigation of whether such parameters become dominant or not.
References
3. Bollen, M., & Rönnberg, S. (2017). Hosting capacity of the power grid for renewable electricity
production and new large consumption equipment. Energies. https://doi.org/10.3390/
en10091325.
4. Jothibasu, S., Dubey, A., & Santoso, S. (2016). Integration photovoltaic generation. Research
Gate. https://doi.org/10.13140/RG.2.2.33939.37928.
5. Chathurangi, D., Jayatunga, U., Rathnayake, M., Wickramasinghe, A., Agalgaonkar, A., &
Perera, S. (2018a). Potential power quality impacts on LV distribution networks with high
penetration levels of solar PV. Paper presented at 18th international conference on Harmonics
and Quality of Power, Ljubljana, Slovenia.
6. Chathurangi, D., Jayatunga, U., Siyambalapitiya, T., Wickramasinghe, A., Perera, S., &
Agalgaonkar, A. (2018b). Connection of solar PV to LV networks: Considerations for maximum
penetration level. Paper presented at Australasian Universities Power Engineering conference on
transition to a low-carbon energy future, University of Auckland, New Zealand.
7. Ismael, S. M., Aleem, S. H. E. A., Abdelaziz, A. Y., & Zobaa, A. F. (2019). Review: State-of-the-
art of hosting capacity in modern power system with distributed generation. Renewable Energy.
https://doi.org/10.1016/j.renene.2018.07.008.
8. Kersting, W. H. (2012). Distribution system modelling and analysis (3rd ed.). London: Taylor
and Francis Group.
9. Chathurangi, D., Jayatunga, U., Perera, S., & Agalgaonkar, A. (2019). Evaluation of maximum
solar PV penetration: Deterministic approach for overvoltage curtailments. Paper presented at
9th IEEE international conference on Innovative Smart Grid Technologies, University of
Politehnica, Bucharest, Romania.
Chapter 8
Hosting Capacity Maximization Based
on Optimal Reconfiguration of Distribution
Networks with Optimized Soft Open Point
Operation
8.1 Introduction
methodology; Sect. 8.3 illustrates the problem formulation and a brief description of
the exIWO algorithm; Sect. 8.4 presents the results of the case studies under
investigation; and Sect. 8.5 is dedicated to the conclusions of this research.
In this section, the power flow equations, DG and SOP modeling, and the proposed
DNR methodology are illustrated in detail. The distribution system model is illus-
trated in Fig. 8.1a in the presence of DGs and SOPs, and SOP-based IGBT
connection model is illustrated in Fig. 8.1b.
Fig. 8.1 Distribution system modeling with DG and SOP: (a) power injection model, (b) -
SOP-based IGBT connection model
182 I. M. Diaaeldin et al.
Power flow equations essential to solve the distribution system under study are
expressed as follows:
P2k þ Q2k
Pkþ1 ¼ Pk Pload
kþ1 r k,kþ1 ð8:1Þ
jV k j2
P2k þ Q2k
Qkþ1 ¼ Qk Qload
kþ1 xk,kþ1 ð8:2Þ
jV k j2
P2 þ Q2
jV kþ1 j2 ¼ jV k j2 2ðr k,kþ1 Pk þ xk,kþ1 Qk Þ þ r 2k,kþ1 þ x2k,kþ1 k 2 k ð8:3Þ
jV k j
where Pk and Qk are the injected active and reactive powers at the kth bus,
load
respectively; Pload
kþ1 and Qkþ1 are the active and reactive powers of loads attached to
bus k + 1 loads, respectively; |Vk| is the magnitude of the kth bus voltage; and rk, k + 1
and xk, k + 1 are the branch resistance and reactance, joining between buses k and
k + 1, respectively.
8.2.2 DG Modeling
0 PDG
k DGcap , 8k 2 N bus ð8:4Þ
where PDG
k is the injected DG active power at bus k, DGcap is the maximum DG
capacity limit, and Nbus is the number of buses in the distribution system.
Soft open points were first proposed by Bloemink and Green [16]. Three different
topologies were employed to integrate SOPs, and they are composed of back-to-back
(B2B) voltage source converter (VSC), static synchronous series compensator, and
unified power flow controller [17]. The B2BVSC is used in this study because of its
capabilities to enhance power quality. An SOP is modeled using its injected powers
(active and reactive) at its terminals as shown in Fig. 8.1. In addition, the algebraic
8 Hosting Capacity Maximization Based on Optimal Reconfiguration of Distribution. . . 183
sum of the SOP active powers injected at its terminals and its internal active losses
by its VSCs must equal to zero, as expressed in Eq. (8.5). Thus,
PSOP
I þ PSOP
J þ PSOPloss
I þ PSOPloss
J ¼0 ð8:5Þ
The reactive power limits are provided in Eq. (8.6), and the SOP size limit is
presented in Eq. (8.7). Thus,
QSOP
I
min
QSOP I QSOP I
max
, 8I, J 2 N f ð8:6Þ
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2 SOP 2 ffi
PSOP
I þ QI SSOP
I , 8I 2 N f ð8:7Þ
where Nf is the number of feeders, PSOP I is the SOP’s active power transferred to the
Ith feeder, PSOP
J is the SOP’s active power transferred to the Jth feeder, PSOPloss
I is
the internal power loss of the VSC coupled to the Ith feeder, PSOPlossJ is the internal
active loss of the VSC coupled to the Jth feeder, QSOP I is the SOP’s reactive power
transferred to the Ith feeder, QSOP
J is the SOP’s reactive power transferred to the Jth
feeder, QSOP
I
min
and Q SOP max
I are the lower and upper limits of the SOP’s reactive
power transferred to the Ith feeder, and SSOP I is the maximum size of the planned
SOP. Further, the active loss of each VSC (PSOPloss I and PSOPloss
J ) is expressed in
Eq. (8.8):
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
SOP 2 SOP 2ffi
PSOPloss
I ¼ ASOP
loss PI þ QI , 8I 2 N f ð8:8Þ
where ASOP
loss is the loss coefficient associated to converter internal power loss [18],
which signifies the percentage of leakage in the delivered power between feeders.
Fig. 8.2 DNR procedure: (a) before DNR, (b) after DNR
The overall hosting capacity (HC) is the objective function, needed to be maximized.
HC is formulated as follows:
PN DG DG
Pi
Max HC ð%Þ ¼ PN bus i¼1
loadnormal
100 ð8:9Þ
k¼1 k P
where Ploadnormal
k is the active demanded power at bus k at the normal loading
conditions.
8.3.2 Constraints
In addition to the DG capacity limit expressed in (8.4), SOP reactive power con-
straints given in Eq. (8.6), and SOP size limit given in Eq. (8.7), constraints
concerning slack bus power factor (PFslack), voltage magnitudes, and lines thermal
capacities are formulated as follows:
min PF
PFslack PFslack ð8:10Þ
slack
max
PLi α PDG
i β PSOP
I γ PSOP
J 0, 8i 2 N bus ð8:13Þ
X XN bus
QSOP
I þ QSOP
J Qload
p ð8:14Þ
p¼1
The invasive weed optimization (IWO) algorithm was designed by Mehrabian and
Lucas [20] to deal with nonlinear optimization problems (Fig. 8.3). The motivation
of the IWO is based on the colonization of weed and was further upgraded in 2014
[15] to deal with MINLP problems including discrete and continuous variables. The
effectiveness of the exIWO arises in its fast convergence toward finding the global
optimum and also its ability to solve mixed continuous-discrete problems. The
exIWO is based on the propagation of weeds and their compliance with ecological
conditions. The population is divided into individuals and seeds. First, the popula-
tion assigned by random values. Further, each individual i produces a bunch of seeds
(Si) based on its objective function value (Fi) as follows:
S Smin
Si ¼ Smin þ ðF i F min Þ max ð8:15Þ
F max F min
where Smax and Smin are the upper and lower number of seeds, respectively. Fmax and
Fmin represent the best and worst objective function value of the population. Second,
the exploration phase includes three methodologies: spreading, dispersing, and
rolling down. The spreading hinges on seed dissemination, in which a set of new
individuals is produced by each seed randomly. Besides, the dispersing is known as
the divergence among individuals and their offspring. It signifies the distance
covered by the seed from the parent plant to its final placement. The distance
formulated as a normal distribution function in which its standard deviation is
expressed Eq. (8.15).
186 I. M. Diaaeldin et al.
m
itermax iter
σ iter ¼ ðσ init σ fin Þ þ σ fin ð8:16Þ
itermax
where iter is the present iteration, itermax is the iterations’ maximum limit, σ init and
σ fin are the initial and final values of the standard deviation, and m is a nonlinear
modulation factor. Further, the rolling down resembles the placement of a seed to a
new location based on its objective function value. Finally, the procedure is then
started over again from the exploration phase until the termination criteria are met.
Termination is fulfilled if iter is higher than itermax or the difference between
individuals’ objective function values was lower than a predefined value.
Starting from the aim of this study, i.e., HC maximization resulting from the
installation of more DGs, the main obstacle facing maximization of HC is PFslack
as a result of injecting active power from DGs; the PFslack will decrease, and it will
be limited by PFslack
min . Therefore, increasing DG penetration depends on the reactive
8 Hosting Capacity Maximization Based on Optimal Reconfiguration of Distribution. . . 187
power provided by the slack bus (Qslack), where the initial injected active, reactive,
and apparent powers by the slack bus are Pslack
0 , Qslack
0 , and Sslack
0 , respectively. In
order to increase DG penetration, the following strategy is deployed as follows:
Step 1: DNR is done first to provide the resilience of the network topology.
Step 2: SOPs are installed instead of tie lines after DNR is done to inject reactive
powers to the system, and as a result, the reactive power injected by the slack bus
Qslack
0 will decrease to Qslack
1 , and also PFslack will increase until it reaches its near
maximum/maximum limit PFslack max .
Step 3: DGs are then installed, Pslack0 decreased to Pslack
2 , and PFslack will decrease
again until it reaches its near minimum/minimum limit PFslack min . The proposed
strategy will be tested in the following section.
In this study, the 83-bus distribution system [21] is the system investigated to
maximize its HC. The topology of this system consists of 83 nodes and 96 lines,
as shown in Fig. 8.4. These lines are composed of 96 lines including sectionalized
and tie lines. The voltage limits Vmin and Vmax are set to 0.95 and 1.05 p.u.,
respectively, and I rated
k,kþ1 is set to 310 A. The number of installed SOPs ranges from
one to five, and SSOP
I ¼ SSOP
J are set to 5 MVA and ASOPloss is set to 2% [18]. The results
obtained after installing DGs along with DNR with and without SOPs for the 83-bus
distribution system are shown in Tables 8.1, 8.2, 8.3, and 8.4.
Results showed that applying DNR only did not guarantee a solution due to the
violation of PFslack
min limit at different loading levels. However, installing SOPs
allowed DG installations as it injects reactive power in order to improve PFslack to
be greater than PFslack
min . Besides, it is notable that HC increases proportionally as the
number of SOPs increases. The best HC values resulted from simultaneous DNR
along with DG and SOP installations, displayed in bold in Table 8.1, where five
SOPs are installed at the three loading levels.
Optimal network configuration, capacity, and sites of SOPs and DGs for five
installed SOPs at normal loading conditions are given in Table 8.5. Branch current
and voltage profiles at normal loading conditions for five installed SOPs at normal
loading conditions are given in Figs. 8.5 and 8.6, respectively. An 83-bus compen-
sated system at normal loading condition for five installed SOPs is given in Fig. 8.7.
188 I. M. Diaaeldin et al.
Table 8.5 Optimal network configuration, capacity, and sites of SOPs and DGs for five SOPs at
normal loading condition
84
8
1 53 47
A 9 G
5 7 55
10 96
11 12 13 60 56
B 88 64 H
89 ~ ~
14 22
20 23 24 72 65
C ~ ~ I
18 21
16
16 91
86
28 29 76 73
D J
26 40
~
~ ~ 39
~
34 83 77
E 42 42 K
32 38 41
~
DG node
~
F ~ ~ SOP
43 46
Fig. 8.7 An 83-bus compensated system at normal loading condition for five installed SOPs
8.5 Conclusion
References
18. Diaaeldin, I., Abdel Aleem, S., El-Rafei, A., Abdelaziz, A., & Zobaa, A. F. (2019). Optimal
network reconfiguration in active distribution networks with soft open points and distributed
generation. Energies, 12, 4172.
19. Chiou, J.-P., Chang, C.-F., & Su, C.-T. (2005). Variable scaling hybrid differential evolution for
solving network reconfiguration of distribution systems. IEEE Transactions on Power Appa-
ratus and Systems, 20, 668–674.
20. Mehrabian, A. R., & Lucas, C. (2006). A novel numerical optimization algorithm inspired from
weed colonization. Ecological Informatics, 1, 355–366.
21. Mohamed Diaaeldin, I., Aleem, S. H. E., El-Rafei, A., Abdelaziz, A. Y., & Zobaa, A. F. (2019).
A novel graphically-based network reconfiguration for power loss minimization in large
distribution systems. Mathematics, 7, 1182.
22. Bollen, M., et al. (2017). Consequences of smart grids for power quality: Overview of the
results from CIGRE joint working group C4.24/CIRED. In 2017 IEEE PES innovative smart
grid technologies conference Europe (ISGT-Europe) (pp. 1–6). IEEE. https://doi.org/10.1109/
ISGTEurope.2017.8260116.
Chapter 9
Generation Regulation Control Systems
9.1 Introduction
Power systems are steadily growing to meet the power demand. Increase in demand
is a challenging concern in developing countries. To meet this ever-increasing
demand, enhancement of power generation is required in conventional as well as
non-conventional form. Conventional power generation has its own limitations such
as land requirement and environmental concern for future enhancement. Conven-
tional power generation has its own limitations such as land requirement, environ-
mental associated pollution and exploitation of fossil fuel. At the same time,
worldwide, distributed generations (DGs) are contributing more due to sustainabil-
ity, pollution free environmnet, ease of energy, and load management. It is expected
that the renewable energy capacity will reach 30 GW over the next two decades,
merely 10% new capacity addition of DG during this period [1]. Integration of DG
into grid is in practice for the past few decades in many countries.
Conventional power system consists of thermal, hydro, and nuclear power plants
interconnected with grid, transmission system, and distribution system. The purpose
of controlling generation is to track load change by adjusting the generation auto-
matically to restore the frequency and tie-line power flow within the prescribed limit
using automatic generation control (AGC) in a conventional power system. Under
normal conditions, generators are operated at a scheduled output, and system
frequency is maintained constant [2]. The primary objective of AGC is to minimize
K. Manickavasagam · V. Puttaraj
Department of Electrical Engineering, M. S. Ramaiah University of Applied Sciences,
Bangalore, Karnataka, India
I. Karuppasamy (*)
Department of Electrical and Electronics Engineering, Amrita School of Engineering, Chennai
Campus, Amrita Vishwa Vidyapeetham, Chennai, Tamil Nadu, India
the frequency deviation and tie-line error to zero such that operating costs of power
system are at a minimum [2].
In a power system, change in load is reflected throughout the system as a change
in frequency. In response to change in system frequency, load frequency control
(LFC) is regulating generator output within a prescribed limit. LFC controls the
frequency deviation, as well as tie-line power flow for maintaining system frequency
and scheduled interchange of power with inter-areas within prescribed limits [2].
Penetration of distributed generation (DG) into the power system is quantified by
hosting capacity approach. This approach aids in better understanding the customer
requirement on the system and requirements of the system operator on customer for
reliable operation. The maximum value of DG penetration up to which power system
operates smoothly without any disturbance is termed hosting capacity. This is
determined by performance index of DG penetration limit. Power system operates
satisfactorily if the performance index is less than hosting capacity.
The European Distributed Energy Partnership (EU-DEEP) project had been
started in 2004 for removing technical and non-technical barriers that prevent a
massive deployment of distributed energy resources (DER) in Europe. Hosting
capacity concept was introduced in a systematic approach starting with the definition
of a number of performance indices and performance objectives. Hosting capacity
calculation varies differently based on power quality events. Hosting capacity
calculation of voltage variations is different from hosting capacity calculation of
frequency variations. The hosting capacity varies depending on structure of the
network, type of DG unit, load, climate, and type of power quality event. There is
no existence of general objectives for power quality events. There are no general
accepted limits for interruptions and dips. Some countries use maximum number of
interruption and maximum duration of interruption as indices and limits to determine
the hosting capacity of their networks [28]. The main purpose of this chapter is to
present the mathematical model for incorporating hosting capacity using power flow
equations in large signal and transfer function model in small signal analysis.
In general, power system dynamics is classified as large signal (static) and small
signal (dynamic) analysis [7]. Large signal analysis is used to study major distur-
bances, in which the magnitude of voltage and power may suddenly increase up to
100% of normal operating values [6–10]. Static equations are used in steady-state
power flow analysis to control existing power system by calculating steady-state
powers and voltages at various buses [10–12]. Small signal analysis is used when
minor disturbances occur, in which the magnitude of voltage and power may increase
few percent of normal operating values. Linear differential equations and Laplace
transform analysis are used for finding the solution in small signal analysis. The
renewable energy generation is usually analyzed by small signal analysis which is
dynamic in nature, since the electrical parameters pertaining to DG such as voltage,
9 Generation Regulation Control Systems 197
current, and frequency are always continuously varying in nature and time depen-
dent. The control strategies for DG connected to microgrid are shown in Fig. 9.1.
DG operations in microgrid are classified as (i) grid-connected mode and
(ii) island or autonomous mode [3]. In grid-connected mode, main grid is absorbing
or delivering power to the microgrid similar to an infinite bus [4]. In island mode,
supply and demand have to be managed by DG connected to microgrid without
power from the main grid [5]. The control strategies for island mode are classified
further into two categories: communication based and droop based.
where the supply-demand gap always exists. The majority share of supply-demand
gap balancing is achieved by controlling the output of conventional generation
resources. In addition, with conventional energy sources, DG can be integrated
using AGC. Hence, modification is required in AGC to integrate DG in addition to
conventional energy sources. Research challenges are involved to integrate DG into
the conventional power systems [6]. In practice, frequency deviation is measured in
the power system network by using frequency meter, and the measured data such as
real power, reactive power, system frequency, and voltages at various buses are sent
to the energy control center (ECC). AGC will adjust the governor according to the
measured frequency deviation and tie-line error to generate required power to match
the load in economic way using EDC.
Solar power penetration into the grid is achieved by grid tie inverter. Wind power
generation is used to inject power into the grid when the induction motor speed is
more than synchronous speed. Thus, the DG penetration into the grid is achieved
based on the availability of natural resources. At present, the contribution of DG is
tremendously increasing worldwide due to their inherent advantages [13]. At the
same time, integration of DG into a power system is a challenging task because of
the dynamic nature of the resources [13, 14]. Setting limits to the introduction of DG
to power system network is only possible when there is agreement on performance
indices and objectives. If the power quality requirements on the network operator are
too narrow, then lesser DG can only be connected to power system network. It
necessitates the need of a wider range of performance indices and objectives [29]. In
this section, method of computing frequency deviation, voltage magnitude variation,
and tie-line error between areas from load flow equations have been explained by
mathematical model.
Generator’s prime mover responses according to the load change which has been
initiated by governor action of AGC. The generator model can be represented by the
mathematical model with the following nomenclature:
αi Participation factor of ith generation
B Bias factor setting of AGC regulator constant for area load-frequency
characteristics
F0 Scheduled system frequency in p.u.
ΔF Frequency deviation
Pgi Active power generation
Pgseti Active power generation schedule
ΔPgi Active power generation due to primary and secondary control
Pgmin i Minimum active power generation
Pgmax i Maximum active power generation
PT Actual tie-line power flow
PT0 Scheduled tie-line power flow
ri Speed-droop setting on turbine governor of ith generating unit
9 Generation Regulation Control Systems 199
Si ¼ Pi þ jQi ð9:5Þ
Pi ¼ Pgseti þ Pgi PLi ð9:6Þ
Pgseti ¼ Pgi þ PLi þ Pi ð9:7Þ
Qi ¼ Qgi QLi ; Qgi ¼ QLi Qi
200 K. Manickavasagam et al.
With bus 1 as reference, the equations for computing the voltage and phase angle
of other buses using decoupled Newton-Raphson iterative solution method are
given:
2 3 2 3
ΔP1 2 3 Δx
6 7 ∂P1 =∂x ∂P1 =∂δ2 ∂P1 =∂δN 6 7
6 ΔP 7 6 76 Δδ 7
6 2 7 ∂P =∂x ∂P =∂δ ∂P =∂δ 6 7
7¼6 N 76 2
6 6
2 2 2 2
7 7 ð9:8Þ
6 7 4 56 7
6
6 7 7
4 5 4 5
∂PN =∂x ∂PN =∂δ2 ∂PN =∂δN
ΔPN ΔδN
2 3
ΔQ1 2 3
6 7 ∂Q1 =∂V 1 ∂Q1 =∂V 2 ∂Q1 =∂δNM
6 7
ΔQ2 7 6
6
6 7¼6 ∂Q2 =∂V 1 ∂Q2 =∂V 2 ∂Q2 =∂δNM 7 7
6 7 4 6 7
6 7
5 5
4
∂QNM =∂V 1 ∂QNM =∂V 2 ∂QNM =∂δNM
ΔQN
2 3
ΔV 1
6 7
6 ΔV 7
6 2 7
66
7
7 ð9:9Þ
6 7
4 5
ΔV NM
Elements of the Jacobian matrices J1and J2 are obtained from estimated bus
voltages, x, and calculated powers. ΔQ and ΔV are obtained by the triangularization
method from the set of linear Eqs. (9.10) and (9.11) that are solved. Frequency
9 Generation Regulation Control Systems 201
deviation and tie-line error are treated as Δx is in this work. To observe range of
ΔF and ΔV from minimum to maximum, generation and load are varied accordingly,
and economic allocation of generations is stored in a lookup table.
The following steps are used to calculate the steady-state frequency deviation and
tie-line exchange deviation using decoupled NR method [8, 16, 17] which is one of
the inputs to the fuzzy logic controller (FLC). Flat frequency control (FFC), flat
tie-line control (FTC), and bias tie-line control (BTC) are different types of AGC
control strategies. This analysis is assumed in an FFC environment to integrate DG.
I. Calculation of AGC Components
1. Input power system data related to:
(a) Network, scheduled loads, and primary and secondary control
characteristics
(b) Type of control strategy used (in this case FFC)
2. Form the bus admittance matrix.
3. Evaluate steady-state frequency deviation ΔF.
Subroutine:
(i) Calculate active power mismatch vector [ΔP].
(ii) Check active power mismatch is within specified tolerance (ε) at all the
buses; if [ΔP < ε], then go to step (iv); else, continue.
(iii) Calculate the Jacobian matrix [J1], vector of unknowns [Δx], and update
voltage [Δx] and phase angles [Δδ].
(iv) Calculate the reactive power mismatch vector [ΔQ].
(v) Check reactive power mismatch is within specified tolerance (ε) at all the
buses; if [ΔQ < ε], then go to main routine; else, continue.
(vi) Compute the Jacobian matrix [J2] and vector of unknowns [ΔV], update
the estimated voltage magnitude [V], and go to step (i).
4. Check whether [ΔF < ε]. If ΔF lies within specified limits, print the results.
5. If ΔF is not less than ε, the corrective action ΔPgi is determined by conven-
tional controller, FLC, and artificial neural network controller (ANNC). ΔPgi
is added with existing generation Pgi and repeat the steps from number 3.
In large signal analysis, frequency deviation and voltage deviations are consid-
ered for performance indices. The allowable limit of frequency deviation is
0.05 Hz and voltage magnitude 0.05 p.u. as shown in Fig. 9.2. If hosting capacity
exceeds performance indices, then power system will not perform satisfactorily.
Control action of AGC is performing in between negative to positive limit for
satisfactory action of power system.
202 K. Manickavasagam et al.
αi ¼ ΔPdi =ΔPD ,
Fuzzy logic controller is used for decision-making based on the availability of wind
velocity, solar irradiation, and the calculated steady-state frequency deviation of
power system. The design procedure of FLC is discussed in the following section.
Three inputs have been considered to design FLC, i.e., ΔF calculated from power
system network data using decoupled NR method, solar irradiation, and wind speed.
The output of FLC is the available amount of power generation from solar and wind
resources.
The input ΔF varies from 0.025 to 0 Hz, solar irradiation varies from 0 to 1500
W/m2, and wind velocity varies from 0 to 20 m/s.
The input and output variables are divided into a suitable number of linguistic
variables [17]. Similarly, the output solar and wind power generated is divided
into a suitable number of linguistic variables. Based on the number of linguistic
variables, the rule base is formed. The input and output conditions have been
described by fuzzy sets. For example, “if solar irradiation is very high, wind speed
is very high, and ΔF is zero, then SPV is no power and wind power (WP) is excess.
When this rule is executed power will not be generated from DG.” These rules can be
developed from operator’s practical experience, knowledge, survey results, general
principles, etc.
9 Generation Regulation Control Systems 205
9.3.2.4 Defuzzification
The process of converting a fuzzified value into actual crisp value is known as
defuzzification. A centroid defuzzification method is normally used to represent the
center of gravity as in (9.12) of the fuzzy set.
PR 0
μr H r
u ¼ Pr¼1
R ð9:12Þ
r¼1 H r
The defuzzified values from FLC are the required amount of power generation
from DG to minimize the AGC component.
the network, and so the input vectors are equal to the number of input nodes.
Estimation of connection (synaptic) weights is done by the training set. Measure-
ment of the generalization ability of the network has been performed by the test set.
The synaptic weights are considered as free parameters. Hidden layers are using
nonlinear activation functions such as logistic sigmoid and hyperbolic tangent.
In the proposed method, the hosting capacity of the system [ΔF and ΔV] is
calculated from the steps given in flow chart Fig. 9.3 using load flow equations.
Performance indices are assumed as 0.05 Hz and 0.05 p.u. After penetration of
DG, sudden changes in input parameters such as solar irradiation wind velocity lead
to unwanted fluctuations in voltage, frequency, and power. The hosting capacity
exceeds performance indices and creates a power system stability problem. Sudden
decrease or increase in load also results in achieving higher value of hosting capacity
than performance indices and causes problem in power system. The operator guide-
line has to be fixed by company authories as a policy to choose of their own
performance indices.
All generators are installed with LFC and automatic voltage regulator (AVR) in an
interconnected power system. The controllers are designed to take care of sudden
disturbances in load to maintain frequency and voltage magnitude within the spec-
ified limits. Changes in rotor angle cause disturbance in frequency as well as real
power. Voltage magnitude affects the reactive power in the power system. When
time constant of prime mover and excitation system is compared, prime mover time
9 Generation Regulation Control Systems 207
constant is larger than excitation system time constant. Hence, comparing the effect
of LFC loop with AVR loop, AVR loop is neglected. In practice, AVR loop and LFC
loop are analyzed independently.
Turbine, generator, governor, and load are connected to feedback loop as is shown in
Fig. 9.6 to represent in block diagram an isolated power system. In block diagram,
K and T represent gain and time constant of main component governor (Kg and Tg),
turbine (Kt and Tt), generator ((Kg and Tg), and sudden load change (ΔPD).
Conventionally, PID controllers are used for small signal analysis in AGC. PID
controllers consist of the following elements: The different combinations of P, I, D,
PI, PD, and PID controllers are possible to implement in AGC. The integral
controller is an effective controller for AGC as per the literature. FLC and ANNC
can be designed for controlling AGC of small signal analysis with the same
techniques discussed in Sects. 9.3.2 and 9.3.3.
The corrective action determined by integral controller, FLC, and ANNC has to be
checked with availability of solar and wind power generation. The available DG can
be incorporated with various types of droop control and communication control
strategy as discussed in Sect. 9.2.
f y ¼ f x þ mx Px ð9:13Þ
V y ¼ V x þ nx Qx ð9:14Þ
Real and reactive power flow control is obtained from an angle of voltage in an angle
droop control.
Equations (9.19) and (9.20) give the relation of P–δ and Q–V.
δj ¼ δrat ma Pj P j,rat ð9:17Þ
V j ¼ V rat na Qj Q j,rat ð9:18Þ
where δrat is the rated voltage angle, Vrat is the rated voltage magnitude, and ma, na
are angle droop coefficients.
Angle droop control is utilized in scenario such that DG coupled converters are
used to form an AC microgrid. In angle droop control, output voltage angle of VSC
is instantaneously changed to reach steady state. Droop coefficient values are
determined by DG ratings for voltage regulation.
PB ¼ PL PðPVþWÞ ð9:19Þ
QB ¼ QL QðPVþWÞ ð9:20Þ
Frequency droop ( fdroop) is derived from active power of battery PBattery using
Eq. (9.19). Voltage droop (Vdroop) is derived from reactive power of battery QBattery
using Eq. (9.20). Control relationships for fdroop and Vdroop from Eqs. (9.19) and
(9.20) can be written as
f droop ¼ f at noload m PBðat fullloadÞ PBðat noloadÞ ð9:21Þ
h i
V droop ¼ V at noload n QBðat fullloadÞ QBðat noloadÞ ð9:22Þ
By integrating f droop in Eq. (9.23), the phase angle θ is derived. Voltage and
current controllers are used to configure the internal control loop.
Droop control can be implemented even sources are separated by long distances
with low bandwidth data lines. However, for nonlinear loads, droop control scheme
does not consider current and voltage harmonics. The presence of circulating
currents makes reactive power control difficult. The major disadvantage of using
droop control is the need for a secondary control to compensate steady-state error.
During operation of the system, droop coefficients are normally constant, and
occasionally, re-set will be done according to the system requirement.
210 K. Manickavasagam et al.
In small signal analysis, frequency and voltage deviations [ΔF and ΔV] are consid-
ered for hosting capacity, and performance indices are assumed with allowable limit
of 0.05 Hz and 0.05 p.u. Implementation of FLC by replacing integral controller
for area 1 is shown in Fig. 9.7.
In the same way, multi-area can be controlled by replacing conventional control-
ler shown in Fig. 9.8 for a multi-area system. Individual FLC and ANNC for each
area or single FLC and ANNC for three areas can be developed for a decentralized
multi-area system [30]. If hosting capacity exceeds performance indices, then power
system will not perform satisfactorily. The important task of various controllers is to
maintain performance indices within the hosting capacity.
Fig. 9.8 Incorporation of hosting capacity in small signal analysis in three area system [30]
Fig. 9.13 Single line diagram of electric circuit model of Gasa Island microgrid
Central Decentral
Central
High band Controller Local Local Local
communication link Controller Controller Controller
DC loads
DC loads
The Internet of Things (IoT) is a platform to connect physical world to digital world.
Any objects physically not connected with any type of network can be connected by
IoT. Various sensors are used to quantify physical parameters, and actuators are used
to change the condition of the Internet. Radio-frequency identification (RFID) is
encoded in digital data tags to find, distinguish, compute, and store data using cloud
technology in IoT. The term “Things” refers to the everyday physical objects, and
the term “Internet” represents the interconnection between networks of heteroge-
neous objects. There are many potential research areas of interest, and one of them is
software architecture and solution. A class is formed for a set of IoT solutions to
illustrate the solution space by investigating the similarities and differences of few
IoT solutions. At the same time, it is not possible to identify a complete set of classes
for achieving such a one-size-fits-all architecture.
In literature, most importantly, the following two qualities are discussed by
various authors:
• IoT is a network which referred to singular form with only one IoT.
• Characteristic definitions refer to the things in the IoT.
• These physical objects are installed by innovation for.
• Associate with the IoT.
• State measurement and their changes.
• Communicate with each other for greater value and service.
Since the IoT applications are enormous variety, accomplishing greater value and
service is extremely broad. Greater value and service is achieved by:
• Formation of new functionality, previously not possible without physical network
• Increasing the quality of existing processes with the help of the IoT
The requirement of IoT should have the following property:
• Self-conjuring, always adaptable for change
• Accessible everywhere
• Seamless communication infrastructure with the following:
(i). One-to-one relationship
(ii). One-to-many relationship
(iii). Many-to-many relationships
9 Generation Regulation Control Systems 217
The present state of the IoT can be categorized into connected home, connected
body, associated retail, associated transportation, smart city, industrial application,
E-Health, and smart energy.
Methodology for analyzing the solutions is performed by identifying the follow-
ing variables:
• Physical Entity being measured
• Attribute(s) of the physical element being measured
• Type of IoT connector(s) used in the solution
• Topology of diagram showing the relationship between parts
• Direction of messages passed between the components in a diagram
• Area of application with logic and data storage
• Possibilities of client interaction
• Estimation of the scalability requirement of each component with the following:
(i). Fixed scale
(ii). Potential increase
Fig. 9.18. Control, radio, and timestamp information are also carried out by Rasp-
berry Pi. Mongo data base (DB) stores data received through the LoRa modem and
works with NoSQL method to support the structure suitable for large-capacity big
data.
analysis are discussed. A mathematical model including hosting capacity with all
possible variables can be developed in the future. The requirement of bandwidth of
IoT may be increasing in the future. Multi-agent system (MAS) based control
scheme is used to monitor, coordinate, control and display data of the system
[6]. These tasks are assigned with various agents such as DER agent, control
agent, database agent, and visualizer or user agent in energy control center. Imple-
mentation of MAS will replace the burden on IoT; in turn, hosting capacity may
increase to higher value.
References
1. Basso, T., & Friedman, N. R. (2003). IEEE 1547 national standard for interconnecting
distributed generation: How could it help my facility? Preprint (No. NREL/JA-560-34875).
National Renewable Energy Lab., Golden, Co.(US).
2. Kanniah, J., Tripathy, S. C., Malik, O. P., & Hope, G. S. (1984). Microprocessor-based adaptive
load-frequency control. IEE Proceedings C Generation, Transmission and Distribution, 131(4),
121–128.
3. Radwan, A. A. A., & Mohamed, Y. A. R. I. (2014). Bidirectional power management in hybrid
AC-DC islanded microgrid system. In 2014 IEEE PES General Meeting Conference & Expo-
sition (pp. 1–5). IEEE.
4. Nejabatkhah, F., & Li, Y. W. (2014). Overview of power management strategies of hybrid
AC/DC microgrid. IEEE Transactions on Power Electronics, 30(12), 7072–7089.
5. Majumder, R., Chaudhuri, B., Ghosh, A., Majumder, R., & Ledwich, G. F. (2010). Improve-
ment of stability and load sharing in an autonomous microgrid using supplementary droop
control loop. IEEE Transactions on Power Systems, 25(2), 796–808.
6. Manickavasagam, K. (2014). Intelligent energy control center for distributed generators using
multi-agent system. IEEE Transactions on Power Systems, 30(5), 2442–2449.
7. Elgerd, O. I. (1982). Electric energy systems theory: An introduction. McGraw-Hill,
New York, US.
220 K. Manickavasagam et al.
28. Nayanatara, C., Divya, S., & Mahalakshmi, E. K. (2018). Micro-Grid management strategy
with the integration of renewable energy using IoT. International conference on Computation of
Power, Energy, Information and Communication (ICCPEIC) (pp. 160–165), IEEE.
29. Yang, Y., & Bollen, M. H. J. (2008). Power quality and reliability in distribution networks with
increased levels of distributed generation’ (Report Elforsk, Stockholm, Sweden, 2008).
30. Saha, A., Saikia, L. C., Tasnin, W., Rajbongshi, R., & Saha, D. (2018). Automatic generation
control of multi-area multisource system incorporating distributed generation units and RFB.
2nd international conference on Power, Energy and Environment: Towards Smart Technology
(ICEPE) (pp. 1–6). Shillong, India.
Chapter 10
Exploring the Concept of Hosting Capacity
from an Electricity Market Perspective
Nomenclature
In this section, the notation of the most important terms used throughout this chapter
is presented for a quick reference.
A. Index and Sets
i2I Index and set of conventional generation units
In Subset of conventional generation units connected to node n
n, m 2 M Indexes and set of nodes
Mn Set of nodes connected to node n through a line
t2T Index and set of hours
V Set of decision variables of MPEC model
VLL Set of decision variables of lower level problem
E. Valenzuela
University of Chile, Santiago, Chile
R. Moreno
University of Chile, Santiago, Chile
Imperial College London, London, UK
Instituto Sistemas Complejos de Ingeniería (ISCI), Santiago, Chile
D. Papadaskalopoulos (*) · Y. Ye
Imperial College London, London, UK
e-mail: d.papadaskalopoulos08@imperial.ac.uk
F. D. Muñoz
Instituto Sistemas Complejos de Ingeniería (ISCI), Santiago, Chile
Universidad Adolfo Ibáñez, Santiago, Chile
B. Parameters
ai Quadratic coefficient of operating cost of conventional generation unit
i ($/MW2)
bi Linear coefficient of operating cost of conventional generation unit i ($/MW)
dn, t Demand at node n and time period t (MW)
τt Weighting factor of hour t
n,m
F Capacity of line (n, m) (MW)
g Minimum output of conventional generation unit i (MW)
i
gi Maximum output of conventional generation unit i (MW)
K Annuitized investment cost of renewable generation ($/MW/year)
xn, m Reactance of line (n, m) (p.u.)
βn, t Normalized output of renewable generation at node n and hour t
γ Price of renewable energy certificate ($/MWh)
C. Decision Variables
gi, t Output of conventional generation unit i at hour t (MW)
rn, t Actual output of renewable generation at node n and hour t (MW)
Rn Capacity of renewable generation at node n (MW)
θn, t Voltage angle at node n and period t (rad)
λn, t Market clearing price at node n and hour t ($/MWh)
10.1 Introduction
10.1.1 Background
Electricity systems are currently facing fundamental challenges associated with the
necessary decarbonization measures to tackle climate change. Such measures need
to consider the deregulated setting of the electricity market, where infrastructure’s
owners, particularly in the generation sector, need clear economic incentives to
deploy future expansions.
The continuously increasing levels of greenhouse gas emissions and the associated
environmental concerns have driven governments worldwide to take significant
initiatives for the large-scale integration of renewable generation. In Europe, the
European Commission (EC) has put forward a legally binding target for renewable
energy sources to cover 20% and 27% of the total energy consumption in the European
Union (EU) by 2020 and 2030, respectively [1, 2]. In the USA, many states have opted
for different kinds of policies that promote renewable energy generation. For example,
in the state of Oregon, the Renewable Portfolio Standards (RPS) require that large
utilities supply 25% of electricity from renewable resources by 2050, whereas smaller
utilities must provide between 5% and 10% by 2025 [3]. In Latin America, there are
countries such as Chile, Peru, Colombia, Brazil, and Mexico with specific targets for
renewable generation and other regulatory policies to foster its integration. In Chile,
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 225
for instance, there is a 20% target (i.e., percentage of the electricity demand supplied
by renewables) by 2025 and a 70% target by 2050 [4].
However, this large-scale integration of renewable generation introduces signif-
icant techno-economic challenges to the operation and planning of electricity grids.
First of all, the inherent variability and limited predictability and controllability of
renewable generation imply that conventional generators need to remain in the
system and operate in an inefficient fashion (part loaded and with more frequent
start-up and shut-down cycles) to provide the required balancing services to the
system [5]. Furthermore, integration of inverter-connected renewable generation
reduces the system inertia (which is provided by the kinetic energy stored in the
rotating mass of the conventional generators), which implies that imbalances
between supply and demand change system frequency more rapidly, challenging
the stability of the system [6]. Finally, the large-scale connection of renewable
generation to transmission and distribution grids creates certain network challenges,
such as congestion (due to thermal and stability limits), increased voltage levels, and
increased short-circuit current levels, which threaten the security of these grids [7].
In order to capture and quantitatively analyze these challenges, the concept of
hosting capacity was firstly introduced to determine the maximum limits of distrib-
uted energy resources (DER) that can be connected to a distribution network without
harming power quality [8, 9]. In the same vein, references [10, 11] establish maxi-
mum quantities of distributed generation (DG) capacity, especially PV, that respect
network’s voltage limits. Other relevant papers have investigated the hosting capacity
of distribution grids beyond voltage issues, considering thermal, short-circuit current
and harmonic distortions [12–15]. Apart from the calculation of the hosting capacity,
measures to enhance it have been proposed too in [16–18], including advanced
control and operational measures, network reinforcements, and the installation of
storage facilities and reactive power equipment. Although the vast majority of the
relevant literature is focused on distribution networks [12, 19–24], at a fundamental
level, the concept can be applied to transmission networks too as in [25]. Hence, the
concept of hosting capacity is useful for governments, regulators, and central policy
makers as it expresses the maximum capacity of renewable generation that can be
integrated to electricity grids without breaching its technical limits and threatening its
secure operation and/or without surpassing certain system cost thresholds.
In combination with the abovementioned challenges of decarbonization, we need
to consider that the electricity industry has been deregulated in many countries,
allowing market forces to guide future investments. This deregulation has led to
unbundling of vertically integrated monopoly utilities and the introduction of com-
petition in the generation sector [26]. Chile and the UK are good examples of early
deregulation of electricity markets in the beginning of the 1980s and 1990s, respec-
tively. The deregulation of electricity markets, however, has not been implemented
everywhere, but more and more countries are adhering to this practice since dereg-
ulation has been key to promote cost efficiency and investments [27].
In this context, generation investment planning is not anymore carried out by a
central regulated utility aiming to minimize system costs but relies on profit-driven
decisions of self-interested generation companies, operating within a competitive
electricity market. In this setting, the large-scale integration of renewable generation
226 E. Valenzuela et al.
The motivation behind this work lies in the fact that previous work has investigated
the renewable generation hosting capacity of electricity systems from a technical
perspective, but no work has been carried out in exploring and analyzing the same
concept from a market perspective, considering the abovementioned challenges of
investment cost recovery and overall profitability.
This chapter makes the first attempt to fill this knowledge gap by proposing a
new modeling approach aiming to quantitatively determine the maximum renew-
able generation capacity that can be integrated in electricity grids, while respecting
long-term profitability constraints, which we define as the market hosting capacity
(MHC) of renewable generation. This approach involves a bi-level optimization
model. The upper level aims at optimizing renewable generation investments in
different locations of the network so as to maximize the total renewable generation
capacity integrated, subject to the non-negative long-term profit constraint imposed
by generation owners. The lower level represents endogenously the market clearing
process and determines the generation dispatch and electricity prices.
The proposed model is validated through case studies on a simple 2-node system,
the IEEE 24-node system, and a 42-node model of the Chilean electricity system.
These case studies quantitatively analyze the dependence of the MHC on a number
of factors, including the network capacity, the flexibility of the conventional gener-
ation fleet, and renewable generation subsidies. Interestingly enough, our results
demonstrate that a stronger network does not always increase the MHC.
The rest of this chapter is organized as follows. Section 10.2 details the proposed
modeling approach. Case studies and quantitative results are presented in Sects. 10.3
and 10.4. Finally, Sect. 10.5 discusses conclusions and future extensions of this work.
As mentioned earlier, the MHC concept introduced in this work aims to identify
the maximum renewable generation capacity that an electricity network can host
from a market perspective. This maximum value exists since the more the
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 227
For clarity reasons, next, we declare the main assumptions behind the proposed
model:
• The model assumes a static planning approach with a yearly operation horizon,
implying that the renewable generation capacity is maximized considering a
single, future target year. This yearly operation horizon is discretized and divided
into a number of representative hours. Both investment and operating costs and
revenues are calculated at the same yearly basis.
• The considered electricity market is a pool-based energy-only market which is
cleared by the market operator through the solution of a short-term social welfare
maximization problem. In order to account for the effect of the network, the
market clearing process incorporates a DC power flow model and yields loca-
tional marginal prices (LMP) λn, t for each node n and time period t. Note that the
LMP system can be applied on both transmission and distribution networks [31],
and thus, our results are equally applicable to both as long as energy is priced in
every node following the LMP principles.
• For simplicity reasons and without loss of generality, a single type of renewable
generation technology is considered, but its normalized output profile is generally
different in different nodes of the network in order to account for the locational
availability of renewable energy sources. Its operating costs are assumed zero,
and their output can be curtailed if required.
• Each conventional generation technology is characterized by a different operating
cost as well as different maximum output and minimum output constraints, with
the latter representing the operation of must-run generation technologies.
• The demand side is assumed inelastic.
The proposed bi-level optimization model (Fig. 10.1) aims at maximizing the total
renewable generation capacity subject to long-term profitability constraints and is
formulated as follows:
(Upper Level)
X
max Rn ð10:1Þ
n
subject to
X X
τt ðλn,t þ γ Þr n,t KRn 0 ð10:2Þ
n, t n
(Lower Level)
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 229
X
min V LL ai gi,t þ bi g2i,t ð10:3Þ
i, t
where
V LL ¼ gi,t , r n,t , θn,t ð10:4Þ
subject to
X X θn,t θm,t
d n,t gi,t r n,t þ ¼ 0 : λn,t , 8n, 8t ð10:5Þ
i2I m2M
xn,m
n n
gi gi,t gi : μ þ
i,t , μi,t , 8i, 8t ð10:6Þ
0 r n,t βn,t Rn : ν þ
n,t , νn,t , 8n, 8t ð10:7Þ
The upper level (UL) problem determines the optimal renewable generation
investments at the different nodes of the network so as to maximize the total
renewable generation capacity (10.1). This problem is subject to the renewable
generation’s long-term profitability constraint (10.2), which ensures that its revenue
from the energy market and potential renewable energy certificates [32] is higher
than or equal (i.e., it at least recovers) to its investment cost, as well as the lower level
(LL) problems (10.3), (10.4), (10.5), (10.6), (10.7), (10.8), (10.9), and (10.10). The
latter represents the market clearing process at each representative hour, minimizing
the total cost of the conventional generation units (10.3), subject to nodal demand-
supply balance constraints (10.5) (the Lagrangian multipliers of which constitute the
LMP), the operating constraints of conventional generators (10.6) and renewable
generation (10.7), and network constraints (10.8), (10.9), and (10.10).
These two problems are coupled, since the renewable generation capacity deter-
mined by the UL problem affects the constraints (10.7) of the LL problem, while the
renewable generation dispatch and the LMP determined by the LL problem affect
the constraint (10.2) of the UL problem.
enabled by the continuity and convexity of the LL problem. This converts the
bi-level problem to a single-level MPEC which is formulated as
X
max V Rn ð10:11Þ
n
where
V ¼ Rn , V LL , λn,t , μ þ þ þ þ
i,t , μi,t , νn,t , νn,t , π n,m,t , π n,m,t , ρn,t , ρn,t , φt ð10:12Þ
subject to
(10.2), (10.5), (10.10),
The objective function of the MPEC formulation (10.11) coincides with the
objective function of the UL problem. The set of decision variables of the MPEC
formulation (10.12) includes the decision variables of the UL problem and the LL
problem as well as the Lagrangian multipliers associated with the constraints of the
LL problem. The KKT optimality conditions of the LL problem correspond to
Eqs. (10.13), (10.14), (10.15), (10.16), (10.17), (10.18), (10.19), (10.20), (10.21),
(10.22), and (10.23).
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 231
X X X X θn,t θm,t
r n,t ¼ dn,t gði2I n Þ,t þ ð10:24Þ
xn,m
n, t n, t n, t n, ðm2M n Þ, t
By multiplying both sides of (10.24) by λn,t, the nonlinear term becomes equal to
By multiplying both sides of (10.13) by gi,t, summing for every i and t, and
rearranging some terms, we obtain
X X X
λðn:i2I n Þ,t gi,t ¼ ai þ bi gi,t gi,t þ μþ
i,t gi,t μi,t gi,t ¼ gði2I n Þ,t λn,t ð10:26Þ
i, t i, t n, t
By multiplying both sides of (10.15) by θn, summing for every n and t, making
use of (10.10), and rearranging some terms, we obtain
X ðθn,t θm,t Þ X
λn,t ¼ ρn,m,t þ ρþ
n,m,t F n,m ð10:29Þ
xn,m
n, ðm2M n Þ, t n, ðm2M n Þ, t
232 E. Valenzuela et al.
By substituting (10.27) and (10.29) into (10.25), this bilinear term is replaced by
the following quadratic expression:
X X X
r n,t λn,t ¼ d n,t λn,t ai þ bi gi,t gi,t þ μþ
i,t gi
n, t n, t i, t
X
ρ
n,m,t þ ρþ
n,m,t F n,m ð10:30Þ
n, ðm2M n Þ, t
The second type of nonlinearity involves the bilinear terms in constraints (10.16),
(10.17), (10.18), (10.19), (10.20), (10.21), (10.22), and (10.23), which can be written
in the generic form 0 μ ⊥ p 0, with μ and p representing generic dual and primal
terms, respectively. The linearization approach proposed by Fortuny-Amat and
McCarl in [35] replaces each of these conditions with the set of mixed-integer linear
conditions μ 0, p 0, μ ωMD, p (1 ω)MP, where ω is an auxiliary binary
variable, while MD and MP are large positive constants. The set of decision variables
of the MILP formulation includes the set (10.12) as well as the auxiliary binary
variables introduced for linearizing (10.16), (10.17), (10.18), (10.19), (10.20),
(10.21), (10.22), and (10.23).
Figure 10.2 shows the configuration of the 2-bus network with two conventional
units, two solar power units, two demands, two buses, and one line. This network is a
modified version of the 2-bus Borduria (node 1) and Syldavia (node 2) example
described in [26]. The specific data used in this chapter is shown in Table 10.1. Note
that while demand is the highest in node 2, the variable cost of thermal generation is
the lowest in node 1. Also, node 1 exhibits the highest availability of the solar power
resource. With all of the above, we can anticipate that power transfers will go from
node 1 to node 2.
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 233
For illustrative purposes, the following results are obtained assuming a single
operating hour; hence, both operating and investment costs are calculated hourly.
Overall, we study 2400 cases with different network capacities ranging from 0 to
1200 MW with a step of 0.5 MW, in which we seek to calculate, for each of them, the
maximum capacity of solar power generation that can be installed in a profitable
fashion in both nodes.
Figure 10.3 shows the market hosting capacity of the 2-bus network as a function of
the line capacity. Interestingly, for smaller line capacity values (i.e., between 0 and
400 MW), the hosting capacity decreases, while the line capacity increases. Note this
result is very counterintuitive since, in practice, we usually expand networks to
increase HC. This counterintuitive result is justified because demand in node 2 can
be supplied by higher volumes of low-cost generation in node 1 (thermal and solar
power generation) when network capacity expands, which replaces capacity of solar
power generation in node 2. Note that building more line capacity clearly hampers
business opportunities for investors in node 2, as node 1 features lower fuel costs and
higher availability factors for solar power generation, making these resources more
234 E. Valenzuela et al.
preferable. In this context, Fig. 10.4 shows how power generations change with
network capacity expansions. Note also that the lower investment levels in solar
power generation in node 2 are not replaced by more investments in solar power
generation in node 1 but are rather replaced by existing low-cost thermal generation
in node 1 that can be exported due to network expansions. This clearly reduces the
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 235
overall MHC since the utilization of existing generation capacity is favored over new
investments in renewables.
However, as existing low-cost thermal generation features a limited capacity, for
larger network expansions (above 400 MW), generation in node 2 is replaced by
investments in new solar power generation capacity in node 1 (while thermal
generation in node 1 is producing at its maximum output). Hence, it is demonstrated
that network expansions above 400 MW unlock the significant solar power gener-
ation potential in node 1, increasing the overall MHC. This is also shown in
Fig. 10.4, where, for power transfers above 400 MW, both thermal and solar
power generations in node 2 decrease due to line capacity expansions that facilitate
the installation of new solar power generation capacity in node 1.
Figure 10.5 shows the increase in the price of node 1 and the decrease in the price
of node 2, as the line capacity increases. This is because the output of the low-cost
thermal generator located in node 1 can be increased while network capacity
expands, exporting power to node 2. This low-cost energy produced in node 1, in
turn, displaces costly energy production in node 2, reducing its locational marginal
price. Therefore, the increase in line capacity constantly reduces the economic
attractiveness of node 2 for investors and increases that in node 1.
All of the abovementioned arguments demonstrate that increasing the capacity of
the network does not necessarily allow private investors to integrate higher volumes
of renewables. On the contrary, increasing network capacity can even reduce the
business opportunities for new renewable generation developments, especially in
those areas close to the load centers and with higher electricity prices.
236 E. Valenzuela et al.
This section determines and analyzes the market hosting capacity of different
network configurations, carrying out sensitivity analyses around the capacity of
the network, the flexibility of conventional generation, and the potential subsidies
that may exist to foster a larger penetration of renewable generation.
We have modified the IEEE RTS described in [36] and shown in Fig. 10.6 by
carrying out the following changes:
(a) Reduction of all lines and transformers’ capacities by 50%
(b) Increase of thermal generation capacity in nodes 1 and 2 by 10 (MW)
(c) Reduction of thermal generation capacity in node 7 by 50 (MW)
These changes (a–c), implemented in our base case, are aimed to enhance
congestion in the system and increase prices in the low voltage area. Also, we
consider that solar power generators can be installed in every node at an annuitized
investment cost of 110 k$/MW.year, and all nodes present the same availability
profile of renewable sources across different time periods. To represent the multiple
market clearing processes that occur across a year, various combinations of different
demand and solar power levels are selected through standard clustering techniques,
e.g., K-means (following [37]). For this particular case, we use demand and solar
profiles observed in Chile. The relevant operational cost data is shown in Table 10.2.
We analyze three case studies which demonstrate how the market hosting capacity
changes for:
(a) Four different network capacity levels
(b) Four different levels of flexibility in conventional generation technologies,
which is varied by changing their must-run output
(c) Four different levels of subsidies in the form of a renewables certificate
Table 10.3 shows the market hosting capacity levels of various modified versions of
the IEEE RTS system when line capacities are increased by 10%, 20%, and 30%
with respect to the base case. Interestingly, the most attractive case to maximize
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 237
several sensitivities, we can demonstrate that the efficacy of these subsidies may
depend on network capacity and generation flexibility. This is an important piece
of information for policy makers since subsidies may not present a significantly
positive impact on promoting renewables due to other inefficiencies in terms of
network configuration and level of (in)flexibility of generation plants.
The planning authority in Chile (which is within the regulatory office) has recently
approved a $1.3 billion investment in a new high voltage direct current (HVDC) line
that will interconnect the north of the country, particularly the Atacama Desert
(location with one of the highest solar irradiation levels in the world), with the
load center of the electricity network located in Santiago, where about 30% of the
national electricity demand is placed. This massive transmission network invest-
ment, with a capacity of 2000 MW and a length of 1500 km, has been justified for
shipping low-cost solar power that can be produced in the Atacama Desert to the
main load center of the country, without originating congestions and without the
need for further network investments in the rest of the AC power network corridor.
This important transmission project, however, received several objections from a
number of market participants, who remained concerned about the benefits of the
new asset in question. One of the strongest arguments used against the project was
that solar resources (and other renewable resources) in Chile are present all along the
country, not only in Atacama, and thus, renewable generation can be perfectly
connected closer to the load center without escalating network investment costs.
Furthermore, these market participants officially took a stand against this project,
bringing their arguments to the Expert Panel. The Expert Panel is a special tribunal
created in Chile to resolve disputes either among market participants or between
market participants and the regulator or planning authority [38]. The first time this
project was proposed by the regulator, the Expert Panel favored the position of those
240 E. Valenzuela et al.
market participants who opposed the project. In the second time, nevertheless, the
regulator faced no opposition, and hence, the project should go forward. Despite this,
there are currently various concerns regarding the effectiveness of this new network
investment option, especially in terms of supporting cost-effective integration of
renewable generation.
In this context, we aim to investigate the effect of the proposed HVDC line on
MHC, particularly analyzing whether (or not) more network capacity really drives
higher MHC in a realistic setting like the Chilean electricity market. Note that the
effect of network investments on the ability of the electricity grid to absorb more
renewable generation capacity is unclear due to the changes driven by network
investments on LMPs. So, we use the concept of MHC to determine whether, in
this case, more network investments may lead to a higher penetration of
renewables.
Figure 10.7 illustrates the Chilean transmission system, which is modeled by using a
40-node representative network with 56 transmission corridors. We focus our
analysis on the year 2018 with a total demand equal to 76 TWh/year and 13 GW
in terms of energy consumption and peak power, respectively. Its generation fleet
contains hydropower (28%), coal (22%), gas (19%), and variable wind and solar
resources (16%), with a total installed generation capacity of 24 GW. Regarding the
market organization, as mentioned in Sect. 10.1, Chile was deregulated in the early
1980s, unbundling generation, transmission, and distribution and establishing gen-
eration as a competitive market while networks remain as regulated monopolies.
Despite this early deregulation, the market clearing in Chile is not bid based but
rather cost based with audited costs. Under this cost-based approach, it is clear that
zero variable cost generation, such as wind and solar power, will drive energy prices
down if their penetration is too large.
We aim to investigate the effects of the HVDC line on MHC in Chile. In Fig. 10.7,
the HVDC line (indicated as a dotted line) goes from node 4 to node 29, bypassing
all existing electrical infrastructure between the Atacama Desert (node 4) and
Santiago (node 29). As the line is of HVDC type, we neglect the effect of the second
Kirchhoff’s law in the model for this particular line, considering the ability to fully
control the power flows through this asset due to the associated flexible converter-
based substations at its both ends. For this model, an additional constraint (to those
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 241
Fig. 10.7 Chilean transmission network in 2018, indicating topology, generation and demand
nodes, and number of circuits per corridor. The HVDC link is depicted as a dotted line
242 E. Valenzuela et al.
indicated in Sect. 10.2) has been added to the hydro reservoirs to make water
management more efficient, limiting the available energy generated per hydropower
plant in addition to the maximum power that each unit can generate, referring to the
limited amount of water available to generate energy.
Apart from applying the proposed model (in which the summation of the wind
and solar power generation capacity invested across the entire system is maximized)
with and without the HVDC link, we run three sensitivity scenarios with different
investment costs for wind and solar power generation, presented in Table 10.6.
These values are annuitized by using a discount rate of 10% and a life span of
25 years.
Finally, the cost data used in this study can be found in [39, 40], and the
availability profiles of renewable generation are specified in [41, 42]. As mentioned
earlier, the richest solar power resources are located in the north of Chile, specifically
in the Atacama Desert.
10.4.2.2 Results
Table 10.7 shows the MHC for the four examined cases, demonstrating that MHC
increases when the HVDC line is built. Furthermore, the MHC is clearly more
sensitive to variations in the investment cost of solar rather than wind power
generation, and this is because the new generation capacity built by the model
corresponds mostly to solar power generation as shown in Fig. 10.8. Indeed, with
the HVDC line, new solar power generation can significantly increase its market
position from 4.7% to 12.7% of the total energy produced. This significant increase
in the penetration of solar power generation due to the installation of the HVDC line
displaces mainly thermal generation from a participation of 40% in the case without
the HVDC line to 34% in the case with the HVDC line. There is also 1.5% decrease
in the participation of hydropower generation when the HVDC line is built, and this
is because the model attempts to maximize new renewable generation capacity,
regardless of the effects of doing so on existing generation capacity.
As demonstrated above, at least for this real case study, expanding network
capacity will increase MHC due to the need to export higher amounts of solar
power generation from the Atacama Desert to the rest of the country, where the
use of more expensive thermal generation will be reduced. In fact, this result is in
line with our previous finding (particularly that on the 2-bus network) that deter-
mines the value of large network capacity expansions in terms of increasing MHC
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 243
50%
45%
40%
Energy partication [%]
35%
30%
25%
20%
15%
10%
5%
0%
(a) (b) (c) (a) (b) (c) (a) (b) (c) (a) (b) (c) (a) (b) (c) (a) (b) (c)
Thermal Hydro Existing wind Existing solar New wing New solar
Fig. 10.8 Energy participation per technology in the base case and (a) without generation
expansion, (b) with generation expansion and no HVDC line, and (c) with generation expansion
and HVDC line
when there are significantly richer renewable resources in a specific area that need to
be exploited. Hence, extra network capacity will favor an increase in MHC in order
to facilitate the shipping of valuable renewable resources to other locations, partic-
ularly large load centers.
This chapter has introduced for the first time the concept of market hosting capacity,
expressing the maximum amount of renewables that can be connected to a given
power network while preserving their profitability in a deregulated environment.
This chapter has also developed a suitable mathematical model to quantify this
market hosting capacity, founded on bi-level optimization principles. This is a new
measure for hosting capacity from a market point of view, which is useful to
understand how various factors, such as network investments, subsidies, and
244 E. Valenzuela et al.
Acknowledgement The research presented in this chapter was supported by grants ANID/
FONDECYT/1181928 and FONDECYT/1190228, the Complex Engineering Systems Institute
(ANID PIA/BASAL AFB180003), SERC-CHILE (ANID/FONDAP/15110019), ANID-Basal Pro-
ject FB0008, and ANID/PIA/ACT192094.
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 245
References
1. European Commission. 2020 climate & energy package [Online]. Available: https://ec.europa.
eu/clima/policies/strategies/2020_en.
2. European Commission. 2030 climate & energy package [Online]. Available: https://ec.europa.
eu/clima/policies/strategies/2030_en.
3. Delmas, M., & Montes-Sancho, M. (2011). U.S. state policies for renewable energy: Context
and effectiveness. Energy Policy, 39(5), 2273–2288.
4. Ministerio de Energia. 2015. Energy 2050: Chile’s energy policy [Online]. Available: http://
www.energia2050.cl/en/energy-2050/energy-2050-chiles-energy-policy/.
5. Dale, L., Milborrow, D., Slark, R., & Strbac, G. (2004). Total cost estimates for large-scale wind
scenarios in UK. Energy Policy, 32(17), 1949–1956.
6. Troy, N., Denny, E., & O’Malley, M. (2010). Base-load cycling on a system with significant
wind penetration. IEEE Transactions on Power Systems, 25(2), 1088–1097.
7. von Appen, J., Braun, M., Stetz, T., Diwold, K., & Geibel, D. (2013). Time in the sun: The
challenge of high PV penetration in the German electric grid. IEEE Power and Energy
Magazine, 11(2), 55–64.
8. Bollen, M. H. J., & Häger, M. (2005). Power quality: Interactions between distributed energy
resources, the grid, and other customers and with distributed energy. Electrical Power Quality
and Utilisation, Magazine, I, 51–61.
9. Papathanassiou, S., Hatziargyriou, N., Anagnostopoulos, P., Aleixo, L., Buchholz, B., Carter-
Brown, C., et al. (2014). Capacity of distribution feeders for hosting DER (Vol. 24). Working
Group C6.
10. Bollen, M., & Hassan, F. (2011). Integration of distributed generation in the power system.
Hoboken: Wiley/IEEE Press. https://doi.org/10.1002/9781118029039.
11. Walla, T. (2012). Hosting capacity for photovoltaics in swedish distribution grids (M.Sc. thesis,
Uppsala University, Uppsala, Sweden).
12. Deuse, J., Grenard, S., Bollen, M. H. J.,Hager, M., & Sollerkvist, F. (2007, May). Effective
impact of DER on distribution system protection. Presented at the CIRED 2007-19th Interna-
tional Conference on Electricity Distribution, Vienna.
13. Etherden, N., & Bollen, M. (2011). Increasing the hosting capacity of distribution networks
by curtailment of renewable energy resources. In 2011 IEEE Trondheim PowerTech. Trond-
heim https://ieeexplore.ieee.org/abstract/document/6019292.
14. Capitanescu, F., Ochoa, L. F., Margossian, H., & Hatziargyriou, N. D. (2015). Assessing the
potential of network reconfiguration to improve distributed generation hosting capacity in active
distribution systems. IEEE Transactions on Power Systems, 30(1), 346–356.
15. Bollen, M. H. J., Sollerkvist, F., Larsson, A., & Lundmark, M. (2006). Limits to the hosting
capacity of the grid for equipment emitting high-frequency distortion. In Proceedings Nordic
distribution and asset management conference, NORDAC 2006. Stockholm, Sweden.
16. Etherden, N. (2014). Increasing the hosting capacity of distributed energy resources using
storage and communication (Ph.D. thesis, Lulea University of Technology, Lulea, Sweden).
17. Sun, W. (2015). Maximizing renewable hosting capacity in electricity networks (Ph.D. thesis,
The University of Edinburgh, Edinburgh, Scotland).
18. AlAlamat, F. (2015). Increasing the hosting capacity of radial distribution grids in Jordan
(B.Sc. thesis, Uppsala University, Uppsala, Sweden).
19. Bollen, M., & Rönnberg, S. (2017). Hosting capacity of the power grid for renewable electricity
production and new large consumption equipment. Energies, 10(9), 1325.
20. Bollen, M., & Hassan, F. (2011). IEEE press series on power engineering. In Integration of
distributed generation in the power system (pp. 88–91). Hoboken: Wiley.
21. Deuse, J., Grenard, S., & Bollen, M. (2008). EU-DEEP integrated project - technical implica-
tions of the “hosting-capacity” of the system for DER. International Journal of Distributed
Energy Resources, 4(1), 17–34.
246 E. Valenzuela et al.
22. Bollen, M., Ribeiro, P., Gu, I., Duque, Y. H., & Carlos, A. (2010). Trends, challenges and
opportunities in power quality research. European Transaction on Electrical Power, 20(1),
3–18.
23. Haesen, E., Minne, F., Driesen, J., & Bollen, M. (2005). Hosting capacity for motor starting in
weak grids. In International conference on future power systems. Amsterdam https://ieeexplore.
ieee.org/abstract/document/1600564.
24. Bollen, M., & Häger, M. (2005). Power quality: Interactions between distributed energy
resources, the grid, and other customers. Electrical Power Quality and Utilisation Magazine,
1(1), 51–61.
25. Sáiz-Marín, E., Lobato, E., & Egido, I. (2018). New challenges to wind energy voltage control.
Survey of recent practice and literature review. IET Renewable Power Generation, 12(3),
267–278.
26. Kirschen, D., & Strbac, G. (Eds.). (2004). Fundamentals of power system economics,. Univer-
sity of Manchester Institute of Science & Technology (UMIST), UK (p. 150). Hoboken: Wiley.
27. Joskow, P. (2003). Electricity sector restructuring and competition: Lessons learned. Cuadernos
de Economia, 40(121), 548–558.
28. Green, R., Vasilakos, N., & Kensington, S. (2011). The long-term impact of wind power on
electricity prices and generating capacity. In 2011 IEEE power and energy society general
meeting. Detroit https://ssrn.com/abstract=1851311.
29. Ketterer, J. (2012). The impact of wind power generation on the electricity price in Germany.
Energy Economics, 44, 270–280.
30. Pollitt, M., & Anaya, K. (2015). Can current electricity markets cope with high shares of
renewables? A comparison of approaches in Germany, the UK and the State of New York.
EPRG Working Paper 1519, University of Cambridge.
31. Sotkiewicz, P., & Vignolo, M. (2006). Nodal pricing for distribution networks: Efficient pricing
for efficiency enhancing DG. IEEE Transactions on Power Systems, 21(2), 1013–1014.
32. Johnstone, N., Haščič, I., & Popp, D. (2010). Renewable energy policies and technological
innovation: Evidence based on patent counts. Environmental and Resource Economics, 45(1),
133–155.
33. Pozo, D., Sauma, E., & Contreras, J. (2017). Basic theoretical foundations and insights on
bilevel models and their applications to power systems. Annals of Operations Research, 254
(1–2), 303–304.
34. Ye, Y., Papadaskalopoulos, D., Moreira, R., & Strbac, G. (2017, Sep). Strategic capacity
withholding by energy storage in electricity markets. In 2017 IEEE PowerTech conference.
Manchester https://ieeexplore.ieee.org/abstract/document/7981200.
35. Fortuny-Amat, J., & McCarl, B. (1981). A representation and economic interpretation of a
two-level programming problem. Journal of the Operational Research Society, 32(9), 783–792.
36. Grigg, C., Wong, P., Albrecht, P., Allan, R., Bhavaraju, M., Billinton, R., Chen, Q., Fong, C.,
Haddad, S., Kuruganty, S., Li, W., Mukerji, R., Patton, D., Rau, N., Reppen, D., Schneider, A.,
Shahidehpour, M., & Singh, C. (1999). The IEEE reliability test system-1996. A report prepared
by the reliability test system task force of the application of probability methods subcommittee.
IEEE Transactions on Power Systems, 14(3), 1010–1020.
37. Lara, C., Mallapragada, D., Papageorgiou, D., Venkatech, A., & Grossmann, I. (2017). Electric
power infrastructure planning: Mixed integer programming model and nested decomposition
algorithm. European Journal of Operational Research, 271(3), 1037–1054.
38. Rudnick, H., Araneda, J., & Mocarquer, S. (2009, Dec). Transmission planning-from a market
approach to a centralized one-the Chilean experience. In IEEE power and energy society
general meeting. Calgary https://ieeexplore.ieee.org/abstract/document/5275552.
39. Comisión Nacional de Energía. Informe de costos de tecnologías de generación [Online].
Available: https://www.cne.cl/wp-content/uploads/2017/12/Res-Ex-CNE-55-2018.pdf.
40. Diaz, G., Inzunza, A., & Moreno, R. (2019). The importance of time resolution, operational
flexibility and risk aversion in quantifying the value of energy storage in long-term energy
planning studies. Renewable and Sustainable Energy Reviews, 112, 797–812.
10 Exploring the Concept of Hosting Capacity from an Electricity Market Perspective 247
Z
U ZnBr battery technology, 116
Upper level (UL) problem, 229