Moses Kavi Thesis
Moses Kavi Thesis
Thesis
Moses Kavi
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources i
Abstract
Existing distribution feeders and their integrated protection systems are not
designed for high penetration of renewable energy (RE) based distributed energy
resources (DERs). The overcurrent protection systems are designed considering the
passive, unidirectional current flow. However, integration of the RE based DERs
such as PV systems through power electronic inverter interfaces fundamentally
changes the distribution network from passive to active network with bidirectional
current flow. The increased use of inverter interfaced RE based DERs and loads will
result in increased harmonic injection affecting power quality. Moreover, increased
penetration of RE based DERs will reduce the level of fault current magnitude from
the feeder substation source. This will adversely affect the feeder protection system
to provide effective protection as the fault current could fall below the overcurrent
threshold.
Faults in power systems (both in AC and DC system) are inevitable and will
occur at one time or another. Certain fault types, such as high impedance faults (HIF)
in AC systems generate low fault current magnitude as opposed to high fault current
magnitude from common short circuit faults which renders the feeder overcurrent
(OC) protection mechanism ineffective in HIF detection. This type of faults must be
detected and isolated as they can cause fire hazards and increase the risk of
electrocution. The inherent difficulty in HIF detection using OC protection scheme in
medium- (MV) to low voltage (LV) where HIFs are a common occurrence can be
aggravated by penetration of RE based DERs. HIF detection and classification based
on feature extraction rather than simply using current magnitude as a metric for HIF
detection will fail. This is due to low fault current magnitude from HIFs and
moreover, increased penetration of RE based DERs that reduces the fault current
magnitude.
Short circuit faults on the other hand result in large fault current having
potential to cause severe damage to power system apparatus and switchgear as well
as causing instability to the unaffected portion of the power system, thus must be
speedily detected and isolated. Short-circuit fault conditions generate transients in
fault current with an exponentially decaying DC-offset. The DC-offset distorts the
iiSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
fault signal waveform and may compromise the integrity of the relay algorithms such
as those based on fast Fourier transform (FFT) and wavelet transform (WT) thereby
resulting in computational delays in the detection of the fault condition. As the
accuracy and speed of convergence of conventional FFT and WT relies on the
periodicity of the fault current and voltage, their effectiveness under DC-offset and
HIFs are limited. Moreover, most DC-offset suppression techniques utilise parameter
estimation, and can add additional computational delay.
Fault protection systems in DC distribution are at their infancy as compared to
the fault protection systems in AC distribution. Faults in DC systems including DC
side of PV system exhibit characteristics quite different from AC system generally
because of different voltage (V) and current (I) characteristics in DC systems. DC
systems generally suffer from short circuit as well as open circuit faults resulting
from mechanical separation of conductors, and in most cases resulting in sustained
arcing. An overcurrent protection strategy using current magnitude as a threshold
metric is applied for all types of faults in the DC power systems including PV
systems. However, not all fault conditions on the DC system can be adequately
protected using such a strategy. One such fault condition is the DC arc-fault
occurring on the DC systems including the PV system. DC arc-fault can either be a
parallel fault (a short-circuit fault) or a series fault (an open-circuit fault). In PV
systems, the detection mechanism relies on backfed current to detect theses faults.
The nature of the faults, especially the series fault contravenes the logic in its
detection using current as the threshold metric. The difficulty in DC arc-fault
detection is compounded in PV systems, particularly at low irradiance which also
includes night to day transition and partial shading. The fast action of the maximum
power point tracking (MPPT) algorithm to put the system at different MPP operation
also imposes additional difficulties in the task of developing accurate reliable DC
arc-fault detection techniques.
In this thesis, a fault detection and diagnostic tool call the decomposed open-
closed alternating sequence (DOCAS) morphological fault detector (MFD) has been
proposed for application in fault detection in both AC and DC systems. The DOCAS
algorithm is a multistage morphological filter constructed from two nonlinear
Mathematical Morphological (MM) filters called the Morphological Median Filter
(MMF) and the Alternating Sequential Filters (ASF). The MM based technique
analyses the topography of the input signal waveforms by means of a probing signal
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources iii
called the structuring element (SE) in complete time domain. MM has the ability to
detect seemingly insignificant changes in the topography of the signal waveform
being investigated. The DOCAS algorithm uses a decomposed weighted SE to
enhance its performance in fault detection. The designed structure of DOCAS
algorithm allows it to be seamlessly applied in fault detection in both AC and DC
systems without any structural change. The characteristics of the MM technique
make the DOCAS algorithm convenient for the detection and classification of HIFs
as well as DC Arc-Faults in PV systems.
The performance of the DOCAS algorithm has been tested in radial
distribution feeder with PV based RE sources connected as DERs for short circuit
faults studies. From these studies, a strategy for adaptive radial distribution feeder
OC protection with built-in DC offset suppression capability is proposed. The
DOCAS algorithm’s capabilities in HIF detection and classification based on feature
extraction has been tested on various contact surfaces using the IEEE 13 bus test
system. The test results showed that DOCAS is capable of extracting successfully the
two target HIF features including randomness and arc extinction and re-ignition
characteristics. A strategy for HIF detection based on the extraction of the two target
features is proposed. The DOCAS algorithm was tested in a radial distribution feeder
with PV based RE sources as DERs for DC arc-fault detection on the PV side. The
performance of the algorithm has been remarkable with all cases of DC arc-fault
detected under all simulated conditions including low irradiance and changing
maximum power point (MPP).
ivSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
Table of Contents
Keywords .................................................................................................................................. i
Abstract .................................................................................................................................... ii
Table of Contents ......................................................................................................................v
List of Figures ........................................................................................................................ vii
List of Tables ......................................................................................................................... xii
List of Abbreviations ............................................................................................................ xiii
Statement of Original Authorship ...........................................................................................xv
Acknowledgements ............................................................................................................... xvi
Publications .......................................................................................................................... xvii
Chapter 1: Introduction ...................................................................................... 1
1.1 Overview of Power System Protection ...........................................................................2
1.2 Conventional Distribution Network ...............................................................................5
1.3 Problem Statement ........................................................................................................11
1.4 Research Aims ..............................................................................................................13
1.5 Research Questions.......................................................................................................14
1.6 Research Contributions .................................................................................................15
1.7 Significance of the Research in Bushfires and Wildfires Prevention ...........................16
1.8 Thesis Outline ...............................................................................................................18
Chapter 2: Literature Review ........................................................................... 21
2.1 Introduction ..................................................................................................................21
2.2 Radial Distribution Feeder OverCurrent Protection .....................................................22
2.3 High Impedance Faults (HIF) Detection and Classification .................................................34
2.4 DC Arc-Fault Detection in Photovoltaic Systems ........................................................41
2.5 Summary and Implications ...........................................................................................49
Chapter 3: Designing the Multistage MM Arc Fault Detection Algorithm . 54
3.1 Introduction ..................................................................................................................54
3.2 Research Methodology Utilizing the MM Technique ..................................................54
3.3 Background of MM Based Techniques ........................................................................55
3.4 design of the morphological algorithm for power system fault detection ....................63
3.5 Attributes of the DOCAS Algorithm ............................................................................71
3.6 Conclusion ....................................................................................................................82
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with
RE Based DERs ................................................................................. 84
4.1 Introduction ..................................................................................................................84
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources v
4.2 Thevenin Equi valent Parameter Esti mation .................................................. 85
4.3 Effect of PV System Penetration on Feeder Substation Fault Current Level .............. 88
4.4 Effect of Fault Location on Feeder Substation Fault Current Level ............................ 91
4.5 DOCAS Algorithm in Adaptive Overcurrent Protection of Radial Distribution Feeder
with PV Penetration ..................................................................................................... 94
4.6 Application of MFD Output in Adaptive Radial Distribution Feeder OC Protection 102
4.7 Simulations and Discussion ....................................................................................... 104
4.8 Inverse-Time Over Current Relaying Using MFD Output Signal ............................. 122
4.9 Conclusion ................................................................................................................. 126
Chapter 5: HIF Detection and Classification in Distribution Feeders ............................. 129
5.1 introduction ................................................................................................................ 129
5.2 Proposed Method for HIF detection and Classification................................................................ 130
5.3 Application of the MFD Output in Detection and Classification of HIF ................... 132
5.4 Simulations and Discussions ...................................................................................... 136
5.5 Challenges in HIF Detection in the Presence of Increasing RE based DER in Radial
DistrIbution Networks................................................................................................ 156
5.6 Conclusion ................................................................................................................. 164
Chapter 6: DC Arc-Fault Detection in PV Systems ...................................... 167
6.1 Introduction ................................................................................................................ 167
6.2 The Proposed DC Arc-Fault Detection Technique .................................................... 167
6.3 Simulations and Discussions ...................................................................................... 169
6.4 Conclusion ................................................................................................................. 183
Chapter 7: Conclusions and Future Directions ............................................. 185
7.1 Summary of Conclusions ........................................................................................... 185
7.2 Future Directions........................................................................................................ 188
Bibliography ........................................................................................................... 205
viSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
List of Figures
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources vii
Figure 3.11: DOCAS response to DC arc-fault in PV systems, (a) DC arc-fault
voltage, (b) Average MMF output, (c) diff DC fault voltage, ΔV
and (c) MFD output. ............................................................................. 82
Figure 3.12: DOCAS response to DC arc-fault in PV systems, (a) DC arc-fault
current, (b) Average MMF output, (c) diff DC fault current, ΔI
and (c) MFD output. ............................................................................. 82
Figure 4.1: A typical radial distribution feeder with PV penetration ......................... 85
Figure 4.2: PV system sequence networks ................................................................ 87
Figure 4.3: Per phase circuit diagram of the typical radial feeder system with
PV ............................................................................................................ 90
Figure 4.4: Per phase circuit diagram categorising PV sources into downstream
and upstream sources with respect to point of fault. ............................... 94
Figure 4.5: The OC fault detection and diagnostic scheme incorporating the
DOCAS algorithm. .................................................................................. 95
Figure 4.6: Flowchart showing the OC fault detection process ............................... 104
Figure 4.7: Test feeder for modelled in Simulink for simulations ........................... 105
Figure 4.8: Characteristic curves for the PV strings at STC, (a) I-V and (b) P-V
curves..................................................................................................... 107
Figure 4.9: Circuit topology of a DC-DC boost converter [171] ............................. 107
Figure 4.10: Switching waveforms for the voltage and current in the DC-DC
boost converter .................................................................................... 108
Figure 4.11: Characteristic curves for the PV strings at STC, (a) I-V and (b) P-
V curves with increased temperatures ................................................ 109
Figure 4.12: Circuit topology of two-level voltage source converter with a
phase output voltage waveform [179] ................................................. 112
Figure 4.13: Fault current signals with DC-offset, (a) The fault current signal
and average MMF output (b) Difference fault current signal. ............ 113
Figure 4.14: Power spectral density plots for (a) fault current (b) average MMF
output and (c) difference fault current, ∆i. .......................................... 113
Figure 4.15: Difference fault current signal, ∆i and the EWMA filter output ......... 114
Figure 4.16: Power spectral density of the EWMA filter output ............................. 114
Figure 4.17: RLSE filter signals (a)RLSE filter input signal, (b) RLSE filter
output signal ........................................................................................ 114
Figure 4.18: RLSE filter magnitude response for SLG fault at fault location 2 at
0%, 28% and 33% PV penetration ...................................................... 117
Figure 4.19: MFD output corresponding to fault current magnitude for SLG
fault at fault location 2 at 0%, 28% and 33% PV penetration. ........... 117
Figure 4.20: MFD output corresponding to fault current magnitude for SLG
fault at fault location 2 at 0%, 28% and 33%, 44% and 50% PV
penetration levels. ............................................................................... 119
viiiSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
Figure 4.21: RLSE filter magnitude response for SLG fault at fault locations 1,
2 and 3 at 28% .................................................................................... 120
Figure 4.22: MFD output corresponding to fault current magnitude for SLG
fault at fault locations 1,2 and 3 at 28% PV penetration. ................... 120
Figure 4.23: MFD output corresponding to fault current magnitude for SLG
fault at fault location 2 at (a) 0% and 28% PV penetration. ............... 122
Figure 4.24: Standard moderately inverse ITOC relay curve with M values in
Table 4.10 ........................................................................................... 124
Figure 4.25: ITOC relay curves at various TDS values ........................................... 126
Figure 5.1: Structure of the Morphological HIF detector ........................................ 130
Figure 5.2: MFD fault windows partitions for HIF detection .................................. 131
Figure 5.3: Flowchart of Proposed HIF Detection and Declaration ........................ 136
Figure 5.4: IEEE 13 bus test system ........................................................................ 137
Figure 5.5: Emanuel Arc model in HIF simulation.................................................. 138
Figure 5.6: V-I characteristic curves of the simulated contact surfaces .................. 139
Figure 5.7: MFDvWindow spikes for SLG fault ......................................................... 141
Figure 5.8: MFDvWindow spikes for capacitor switching .......................................... 141
Figure 5.9: MFDvWindow spikes for induction motor switching ............................... 141
Figure 5.10: MFDvWindow spikes for step load increase ........................................... 141
Figure 5.11: Signals for HIF at 602 on conc. surface (a) fault voltage and HIF
current, and (b)MFDv output .............................................................. 146
Figure 5.12: MFDvTall edge spikes for HIF at 602 on conc. surface (a) MFDvTall
(1), (b) MFDvTall(1) pu increase,(c)MFDvTall(2) and (d)
MFDvTall(2) pu increase. ..................................................................... 146
Figure 5.13: MFDvShort edge spikes for HIF at 602 on conc. surface (a)
MFDvShort (1), (b) MFDvShort(1) pu increase,(c)MFDvShort(2) and
(d) MFDvShort(2) pu increase. .............................................................. 147
Figure 5.14: Random MFDvWindow spikes for HIF at 602 on conc. surface ............. 147
Figure 5.15: MFDvArc spikes for HIF at 602 on conc. surface, (a) MFDvArc
spikes and ............................................................................................ 148
Figure 5.16: MFDvTall edge spikes for HIF at 605 on dry grass (a) MFDvTall (1),
(b) MFDvTall(1) pu increase,(c)MFDvTall(2) and (d) MFDvTall(2) pu
increase. .............................................................................................. 152
Figure 5.17: MFDvShort edge spikes for HIF at 605 on dry grass (a) MFDvShort
(1), (b) MFDvShort (1) pu increase, (c)MFDvShort(2) and (d)
MFDvShort(2) pu increase.. ................................................................... 152
Figure 5.18: Random MFDvWindow spikes for HIF at 605 on dry grass .................... 152
Figure 5.19: MFDvArc spikes for HIF at 605 on dry grass, (a) MFDvArc spikes
and ....................................................................................................... 153
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources ix
Figure 5.20: Signals for HIF at 602 on conc. surface (a) fault voltage and HIF
current, and (b) MFDv output with noise. ........................................... 154
Figure 5.21: MFDvTall edge spikes for HIF at 602 on conc. surface (a) MFDvTall
(1), (b) MFDvTall(1) pu increase,(c)MFDvTall(2) and (d)
MFDvTall(2) pu increase with noise. .................................................... 155
Figure 5.22: MFDvShort edge spikes for HIF at 602 on conc. surface (a)
MFDvShort (1), (b) MFDvShort(1) pu increase,(c)MFDvShort(2) and
(d) MFDvShort(2) pu increase with noise. ............................................. 155
Figure 5.23: Random MFDvWindow spikes for HIF at 602 on conc. surface with
noise. ................................................................................................... 155
Figure 5.24: MFDvArc spikes for HIF at 602 on conc. surface, (a) MFDvArc
spikes and ............................................................................................ 156
Figure 5.25: Radial feeder with PV .......................................................................... 157
Figure 5.26: DOCAS MFDi outputs for HIF at XF1 closer to the feeder
substation (a) MFDi without PV, (b) MFDi with PV .......................... 158
Figure 5.27: DOCAS MFDi outputs for HIF at XF2 further from the feeder
substation (a) MFDi without PV, (b) MFDi with PV .......................... 158
Figure 5.28: HIF signals for HIF at XF2 on dry grass (a) fault currents for with
and without PV, (b) MFDi outputs for with and without PV .............. 161
Figure 5.29: MFDiTall spikes for HIF at XF2 on dry grass (a) MFDiTall (1) spikes
for with and without PV, (b) MFDiTall (2) spikes for with and
without PV .......................................................................................... 162
Figure 5.30: MFDiTall spikes and MFDiTall pu increases for HIF on dry grass at
XF2 with PV (a) MFDiTall (1) spikes, (b) MFDiTall (1) pu increase,
(c) MFDiTall (2) spikes, (c) MFDiTall (2) pu increase ........................... 162
Figure 5.31: HIF signals for HIF at XF2 on dry grass (a) fault voltage with PV,
(b) MFDv output .................................................................................. 162
Figure 5.32: MFDvTall spikes and MFDvTall pu increases for HIF on dry grass at
XF2 with PV (a) MFDvTall (1) spikes, (b) MFDvTall (1) pu increase,
(c) MFDvTall (2) spikes, (c) MFDvTall (2) pu increase .......................... 163
Figure 5.33: MFDvShort spikes and MFDvShort pu increases for HIF on dry grass
at XF2 with PV (a) MFDvShort (1) spikes, (b) MFDvShort (1) pu
increase, (c) MFDvShort (2) spikes, (c) MFDvShort (2) pu increase ........ 163
Figure 5.34: Window, MFDvWindow spikes ................................................................ 163
Figure 5.35: HIF arc extinction and re-ignition, MFDvArc spikes ............................ 164
Figure 6.1: Block diagram of the DC Arc-Fault detection system........................... 168
Figure 6.2: Radial distribution feeder with PV penetration used in the
simulation study .................................................................................... 169
Figure 6.3: PV array configuration in the simulation system .................................. 170
Figure 6.4: A typical PV system configuration array configuration ....................... 170
Figure 6.5: Characteristic curves for the PV String configuration (a) I-V curve,
and (b)P-V curve ................................................................................... 171
xSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
Figure 6.6: Layout of the PV modules in the PV strings for fault simulations. ....... 171
Figure 6.7: Current signal measured at DC bus with associated signals for DC
arc-fault FP1 ........................................................................................... 176
Figure 6.8: Voltage signal measured at DC bus with associated signals for DC
arc-fault FP1 ........................................................................................... 176
Figure 6.9: MFD outputs for PV string currents for fault DC arc-fault FP1 ............ 176
Figure 6.10: Current signal measured at DC bus and MFD out for DC arc-fault
FP7 ....................................................................................................... 177
Figure 6.11: Voltage signal measured at DC bus and MFD out for DC arc-fault
FP7 ....................................................................................................... 177
Figure 6.12: Current signals measured at PV strings for DC arc-fault FP7 ............. 178
Figure 6.13: MFD outputs for PV string currents for DC arc-fault FP7 .................. 178
Figure 6.14: Current signal measured at DC bus and MFD out for DC arc-fault
Fs2 ....................................................................................................... 179
Figure 6.15: Voltage signal measured at DC bus and MFD out for DC arc-fault
Fs2 ....................................................................................................... 179
Figure 6.16: Current signals measured at PV strings for DC arc-fault Fs2 .............. 179
Figure 6.17: MFD outputs for PV string currents for DC arc-fault FP2 .................. 180
Figure 6.18: DC arc-fault current and MFD output high to low transition after
fault ..................................................................................................... 181
Figure 6.19: DC arc-fault voltage and MFD output high to low transition after
fault ..................................................................................................... 181
Figure 6.20: DC arc-fault current and MFD output high to low transition before
fault ..................................................................................................... 182
Figure 6.21: DC arc-fault voltage and MFD output high to low transition
before fault .......................................................................................... 182
Figure 6.22: DC arc-fault current and MFD output high to low transition with
fault ..................................................................................................... 183
Figure 6.23: DC arc-fault voltage and MFD output high to low transition with
fault ..................................................................................................... 183
Figure 7.1: A typical DC power distribution system ............................................... 199
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources xi
List of Tables
xiiSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
List of Abbreviations
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources xiii
PV Photovoltaic
RE Renewable Energy
RLSE Recursive Least Square Error
RMPP Resistance at Maximum Power Point
RMS Root Mean Square
SE Structuring Element
SS Substation
STC Standard Test Condition
TDS Time Dial Setting
TMS Time Multiplier Setting
VMPP Voltage at maximum Power Point
VT Voltage Transformer
WT Wavelet Transform
xivSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
Statement of Original Authorship
The work contained in this thesis has not been previously submitted to meet
requirements for an award at this or any other higher education institution. To the
best of my knowledge and belief, the thesis contains no material previously
published or written by another person except where due reference is made.
Date:
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources xv
Acknowledgements
Last but certainly not the least; I would like to specially thank my wife and my
children for having the patience to put up with my long absence while I pursued my
PhD research. I would like to also thank God, my mother and my father for bringing
me into this world.
xviSmart Protection System for Future power System Distribution Networks with Increased Distributed Energy Resources
Publications
Smart Protection System for Future power System Distribution Networks with Increased Distributed Energy
Resources xvii
Chapter 1: Introduction
This chapter introduces the topic of the thesis and provides background
discussion related to the topic including general overview of distribution network
protection, its components and their functional attributes. The technical challenges
introduced due to the increased penetration of RE based DERs in distribution
network feeder protection are discussed. The discussions further highlight the
coexistence of AC and DC systems within the distribution network made possible by
the advances in power electronic converters and interfacing technology, and the need
to provide effective protection in both systems.
The current trend in the design and utilisation of modern distribution networks
includes diverse generating sources including renewable energy (RE) sources
directly connected to the distribution feeders as distributed energy resources (DERs).
The inclusion of RE based DERs changes the passive unidirectional power delivery
nature of the radial distribution networks to active networks with bidirectional
current flow. The conventional OC protection system in radial distribution feeder
relies on current magnitude as a threshold metric, thus it is imperative that sufficient
fault current magnitude above the predefined threshold limit during fault must exist
for this scheme to be effective. The direct interconnection of RE based DERs at the
distribution feeder contravenes this fundamental requirement as the DERs supply
power to distributed loads along the feeder length thereby reducing the current
supply emanating from the feeder substation. This can affect the coordination of the
protective devices in the feeder OC protection scheme hence compromising its
effectiveness in responding to fault conditions reliably.
Furthermore, other fault conditions such as high impedance faults (HIFs)
resulting from fallen conductors as well as energised conductors making unwanted
contacts with tree trunks and branches are quite common in medium (MV) and low
voltage (LV) distribution networks. Unlike OC faults, this category of faults
generates a low current magnitude rendering the feeder OC protection scheme
ineffective in detecting such faults. This inherent difficulty can be exacerbated by
increased levels of RE based DERs in the distribution feeders.
Chapter 1: Introduction 1
In recent times, the advances in power electronic converters and interfacing
technology has enabled the creation of DC subsystems within the AC distribution
networks to directly supply DC power to DC loads from DC sources such as
photovoltaic(PV) based RE sources or through AC-DC inverters from traditional AC
sources. This adds further challenges in designing protection systems for DC systems
as the protection mechanism in AC systems cannot be applied in DC systems.
These technical challenges introduced due to the increased penetration of RE
based DERs define the scope of the research problem. The aims and objectives of the
research, the research questions and the research contributions are also defined in this
chapter.
Speed: The protection system must operate rapidly to interrupt the fault to
minimise damage or possible system collapse. Intentional delay can be
introduced as part of the relaying strategy to coordinate between zones
of protection. The minimum time taken to isolate hazard is called the
clearing time.
Sensitivity: This ensures that fault (current/voltage) magnitude (however small)
should be detected by the protection system. The protection system
must be able to correctly discriminate between fault and normal
operating condition based on comparison with some predefined
inequality constraint or threshold quantity.
Selectivity: This is defined by the relaying or switching strategy where only the
portion of the network or equipment under fault is isolated. Relaying
strategy sectionalises the network into zones of protection whereby only
faulty equipment or network portion within a zone of protection are
isolated. Zones of protections are defined by proper grading of
2 Chapter 1: Introduction
protection threshold, time delay and/or operating characteristics of the
protective relaying devices.
Security: Security of the power distribution network operation means that the
protection system must be reliable. Moreover, robustness is also
anticipated as it would add to the notion of reliability of the protection
system. Reliability defines the expectation that the protection system
devices will correctly operate when expected to. Robustness suggests
that the devices can operate reliably under dynamic network condition.
The security of the protection system further is enhanced by having
backup protection. If a protection device in a zone fails to operate, then
backup protection is provided by devices in the neighbouring zones.
The protection system continuously monitors the state of the power system by
taking measurements of quantity such as voltage (V´) and current (I´) via the
secondary windings of the VT and CT respectively. The voltage (V) and current (I)
signals are the primary quantities, other quantities, such as impedance, admittance
Chapter 1: Introduction 3
and power, etc., are secondary quantities can be computed from the primary
quantities if required to execute a decision. The functional attributes of the protection
components are graphically illustrated in Figure 1.2. The occurrence of any
abnormality or disturbance is detected by comparing the measured and/or computed
quantities against a preset threshold metric. The network is said to be in normal
(healthy) state and is allowed to continue operating if the measured quantities and/or
their derivatives are within the limit defined by the preset threshold metric. If the
value of the measured quantities and/or their derivatives violates the limit of the
threshold metric, the network is said to be in abnormal (unhealthy) state. The
protection system goes into state of alert to indicate abnormal network state, and the
decision logic element is activated to issue a trip command to take out the faulty
system apparatus or isolate the faulty section of the network and allow service to the
rest of the network.
The functional block diagram in Figure 1.2 shows that, the protection relay is
the key component in defining the functional attributes of the protection system. The
protection relay executes the decision to send appropriate trip signal to the associated
circuit breaker to isolate the unhealthy part of the network and allow service
continuity to the rest of network. The relay senses the existence or occurrence of any
abnormal condition, then in collaboration with the circuit breaker to isolate the
disturbances in a fast and reliable manner to minimise any harmful consequence to
the power system as well as human beings. Thus, to ensure the security, integrity and
reliability of a power system, it is highly imperative to have a fast and effective relay
operation that is sensitive and respond rapidly to hazardous condition to isolate any
faulted equipment or section of the network under any operating conditions and
allow the continuous operation of the healthy section of the network.
4 Chapter 1: Introduction
1.2 CONVENTIONAL DISTRIBUTION NETWORK
Chapter 1: Introduction 5
Considering the system referred to in Figure 1.3, the distribution feeders exist
at the part of the distribution subsystem at 10-40 kV. The feeders are radial in nature
and provide means to distribute energy to the consumers as tapped loads through
distribution transformers. A typical example of a radial distribution feeder with step-
down transformers supplying power to the consumers as tapped loads is shown in
Figure 1.4. The Substation source represents the HV primary source in the
distribution substation which is stepped down to 10-40 kV. The consumers are
connected through distribution transformers at lower voltages of between 220-400 V.
6 Chapter 1: Introduction
modern power system design, there are proposals to address the issues of low
efficiency and, high maintenance and operational costs by integrating small scale
power generating systems at the distribution network as DERs [6] coupled to radial
distribution feeders. Direct integration of DERs reduces the requirement for long
transmission lines, which reduces maintenance and operational costs, and improve
energy efficiency by reducing transmission losses and better management of
resources. While achieving improvements in cost reduction and energy efficiency, it
is envisaged that the transformed power system will meet the obligation for reduced
carbon emission. This requirement is inevitable and mandatory, thus must be
incorporated into the design and implementation of the modern power system as
international governments develop and adopt policies to reduce global warming.
Furthermore, with the advancement in power electronic technologies there are
increasing number of load types that require DC power supply. Therefore, to supply
both AC and DC load demands, the modern power system must be more flexible,
and configurable energy system [7]. This view has seen the proliferation in
production and deployment of RE based DERs and energy storage systems (ESSs)
[8] at the distribution network. Some RE technologies that have matured and widely
used include wind turbines [9], photovoltaic (PV) arrays [10] fuel cells and micro-
turbines [11] while ESSs include batteries [12], flywheels [13] and super capacitors
[14]. A power system that exhibits such characteristics with integration of PV system
at the distribution feeder level is depicted in Figure.1.5.
The ever-increasing demand for energy, and more so the need for clean,
sustainable environmentally friendly energy sources has seen the proliferation in the
development and deployment of renewable energy harvesting technologies such as
photovoltaic modules, wind generators, fuel cells, etc., to name a few. Moreover,
with the advancements in the development of power electronic converters and
interfacing technologies, the RE sources can now be easily integrated into medium
voltage (MV) or low voltage (LV) distribution networks as DERs. Figure.1.5 shows
the integration of PV based RE sources as DERs into the radial distribution feeder.
The integration of RE based DERs such as PV systems provides opportunity for
diversification of load types that can be supplied directly or indirectly by the DERs.
The RE based DERs such as PV systems can supply DC loads directly through a DC-
Chapter 1: Introduction 7
DC converter as well as supply AC loads through DC-AC converter. There is
increasing load diversification today due to advances in power electronic devices
where more and more loads today require DC power supply as opposed to the
traditional AC power supply. Thus, the development trend in modern power system
must consider, and be capable of supporting the integration of diverse and increasing
penetration of RE based DERs, energy storage, electric vehicles while meeting the
high power quality standard required by sensitive modern digital devices and loads
[15], [16].
The existing power system, more so the distribution subsystem and its
integrated protection mechanism were not designed for high penetration of DERs,
more so renewable energy resources. The existing distribution network/feeder
overcurrent protection system was designed considering the passive, unidirectional
current flow. However, integration of the RE based DERs fundamentally changes the
distribution network from passive to active network with bidirectional current flow
[17]. Moreover, increased deployment of power electronic converters and loads
increases harmonic injection thus increasing the level of voltage and current
distortion. Furthermore, high penetration of RE based DERs contribute large
component of the total fault current thus reducing the magnitude of the fault current
at the feeder substation. This will adversely affect the main feeder protection relay at
the feeder substation to provide effective protection as the fault current could fall
below the overcurrent threshold [18] ,[19]. Further challenge to distribution network
feeder protection is that increased penetration of the RE based DERs will decrease
the reach of the protective devices within a zone of protection for which it was
configured to operate [20].
Some researchers have proposed that the future distribution network will
incorporate DC buses to connect DC sources and loads [21]. To maintain optimum
system operation, the DC generation sources as well as associated components and
loads must be protected from harmful fault conditions. The AC feeder OC protection
system cannot be used as it is not designed for the DC system protection. Thus, a
separate DC protection system is necessary for protection of DC systems. The
existing DC fault protection scheme which is contingent on detecting sufficient fault
8 Chapter 1: Introduction
current magnitude, [22] cannot guarantee protection against all types of faults in DC
systems. For instance, the DC system OC protection will fail to operate in PV system
under low irradiance as the fault current magnitude would be insufficient to trigger
the OCPD and/or GFPD.
In designing protection systems for AC and DC networks, consideration for
adaptability is a key factor. The dynamics of both the AC and DC systems including
network topology change under different operating conditions with RE penetration
must be considered. In this thesis, the study is based on PV as DER; however, it is
assumed that the analysis, observations and conclusions, and the outcome of the
research will be applicable to all RE based DERs unless specified otherwise. Thus, to
maintain generality, the use of PV based DERs will be dropped, instead RE based
DERs is used.
Chapter 1: Introduction 9
substation and the type of fault. The network topology change influenced by the
intermittency of the RE sources is another parameter that impacts on the fault current
magnitude. Thus, for the feeder OC protection scheme to be effective, the threshold
parameter must be adaptive to all the changes introduced by the RE based DERs.
High Impedance Faults: The power system network not only suffers from low
impedance faults, but also other fault conditions resulting in low fault current
magnitudes known as high impedance faults (HIFs) are prevalent in medium voltage
(MV) and low voltage (LV) networks. Because of the inherent low fault current
magnitude from HIFs, the conventional OC protection will become ineffective under
such fault conditions. HIFs are very difficult to detect because of the low fault
current magnitude as well as their highly random and nonlinear characteristics which
most often involves arcing.
In the face of the changing landscape in the distribution network brought about
by the integration of the DERs, more specifically the RE based DERs, effective
protection system, considering all fault types must be considered.
10 Chapter 1: Introduction
reliable detection of large fault current, above a predefined threshold to operate the
protection devices. The DC arc-faults, particularly those occurring in the DC bus of
the PV system as well as on the PV strings, including inside the PV cells and
modules are difficult to detect, especially under low irradiance. The longer the fault
exists; the optimal operational condition of the DC system is degraded. Moreover,
the persistent and sustained DC arc increases the risk of fire.
The challenges and issues in power system protection considering the changing
landscape in distribution network due to the integration of RE based DERs are
summarised here, and they constitute the problems that motivated this research.
Chapter 1: Introduction 11
1.3.1 Challenges in Overcurrent Protection in Distribution Network with
Increasing DER Penetration
12 Chapter 1: Introduction
▪ The HIF current has high frequency harmonic components from 2 to
10 kHz
The aim of this research is to develop a Fault Detection and Diagnostic Tool
that can be used seamlessly in AC and DC distribution systems. The attributes of this
tool are then used to propose strategies for; Adaptive OC protection in distribution
networks with increased DER penetration, HIF detection and DC Arc-Fault detection
in PV systems which are defined herein.
Chapter 1: Introduction 13
1) Adaptive to the change in landscape of power system structure at the
distribution network.
5) Must be able to detect power system fault and issue trip signal within 1
cycle of the fundamental frequency.
The attributes of the HIF detection and classification strategy would include;
2) Detect Arc fault under all conditions including, low irradiance, partial
shading, night to day transition.
To meet the research aims, the following research questions are defined;
14 Chapter 1: Introduction
1) What is the effective strategy for fault detection and diagnosis that can be
applied in both AC and DC power networks with DC-offset suppression
capability in AC power system fault detection as well as convenience of
application in adaptive overcurrent protection in radial distribution feeders
with increased RE based DER penetration?
2) What are the analytical and computational methods for developing a fault
detection and diagnostic tool that can be seamlessly utilized in both AC and
DC power systems?
3) Are the analytical and computational tools capable of performing feature
extraction for high impedance faults detection?
4) Are the analytical and computational analytical tools capable of performing
feature extraction for DC Arc-faults detection?
-Tool for Fault Detection and Diagnosis: A fault detection and diagnostic tool
based on Mathematical Morphology for time-domain analysis of the fault signal
called the decomposed open close alternating sequence (DOCAS) morphological
fault detector algorithm is proposed. This tool is a multistage filter based on two
nonlinear morphological filters namely; the morphological median filter (MMF) and
the alternating sequential filter (ASF). The MMF filter is comprised of two cascaded
stages where the output of the first stage becomes the input to the second stage while
the ASF has two layers; the open-close and the close-open alternating sequential
filters each with four stages. The two layers of the ASF operate simultaneously, and
each stage of the ASF is cascaded in a hierarchical manner where the output of the
previous stage is cascaded to the next stage. This operational sequence is achieved
through the decomposition of the filtering signal call the structuring element (SE)
into two SEs used in the two different filters. The underlying nature and computation
technique in the DOCAS algorithm makes it possible for its application in both the
AC and DC power system fault detection and diagnosis. The complete process in
developing the DOCAS algorithm is presented in section 3.4.
Chapter 1: Introduction 15
Distribution Feeders with RE based DERs: An analytical method by means of
Thevenin parameter estimation and two distance factors, distance to the RE based
DER and distance to fault, to analysing the impact of the level of DER penetration on
feeder current magnitude is proposed. Then a technique for adaptive pickup setting in
overcurrent relays in distribution feeder overcurrent protection schemes using the
DOCAS MFD output is proposed. Moreover, a technique for determining the relay
trip time using inverse time overcurrent relay based on the adaptive pickup parameter
is proposed. The techniques are presented in Chapter 4.
-Method for the Detection and Classification of HIFs: A technique for the
detection and classification of HIFs based on feature extraction using the DOCAS
algorithm was developed. The HIF features extracted from the DOCAS MFD output
include the HIF randomness and the arc extinction and re-ignition feature. The
feature extraction and HIF detection and classification are presented in Chapter 5.
16 Chapter 1: Introduction
173 lives. The Victorian Bush Fire Royal Commission identified electricity
distribution infrastructure as the cause of these fires [24]. The Royal Commission
based on its investigations made 67 recommendations of which 8 are directly to the
electricity distribution infrastructure. The Victorian government in response to the
recommendations of The Royal Commission established a Powerline Bushfire Safety
Task force to investigate the recommendations directly related to the electricity
distribution infrastructure. The Task force furnished a report on the 30th of
September 2011 to the Victorian government, and out of several recommendations,
one was for further research into improving safety by identifying and introducing
new technology and methods for reducing risks and preventing bushfires ignited by
the electricity distribution system.
Powerlines span several hundreds of kilo meters and are exposed to weather
and often come into contact with trees and vegetation. Fallen powerlines due to
support structure failure and/or powerlines coming into contact with vegetation
almost always result in arcing. Unlike short circuit faults, these conditions are
classified as high impedance faults (HIFs) and generate low fault current magnitudes
resulting in not being detected by the overcurrent protection system. The longer the
arcing fault conditions persist, the higher the risk of igniting fire and electrocution of
people. These conditions must be detected and eliminated before they escalate into
any catastrophic events like bushfire and wildfires.
Bushfires and wildfires are disastrous events that are bound to happen if the
arcing sources such as HIFs are not detected and removed. HIFs generally result in
low fault current magnitude thus making it extremely difficult for their detection
using the conventional overcurrent (OC) protection scheme. Thus, specialise
algorithm specifically for arc-fault detection including HIFs have to be developed
that can work in tandem with the conventional OC protection system. Thus, the
DOCAS algorithm satisfies this requirement. Moreover, the DOCAS algorithm can
be used as part of an online intelligent power line condition monitoring system, for
the detection of arcing events.
Fire in DC systems including PV based DERs is major cause for concern. The
DOCAS algorithm can be seamlessly utilised in DC systems to detect the
occurrences of DC arc-fault as fire prevention mechanism in DC systems including
PV systems.
Chapter 1: Introduction 17
1.8 THESIS OUTLINE
This thesis consists of seven chapters including Chapter 1 which covers the
Introduction where the problem statement, the research objectives, the research
questions and research contributions were defined. The overviews of the rest of the
chapters of the thesis are presented herein:
Chapter 2: This chapter reviews literature on the three topics defined in the
problem statement including, 1) Challenges in overcurrent protection in radial
distribution feeders with increasing penetration of RE based DERs, 2) Challenges in
HIF detection in distribution networks and 3) Challenges in DC arc-fault detection
on the DC bus in PV systems.
A background review on the existing practice in radial distribution feeder OC
current protection and relaying strategy, and the drawback in its non-adaptive nature
in the face of increasing DER penetration are discussed. Moreover, a discussion on
the prevalence of the exponentially decaying DC-offset that occurs during short-
circuits faults and its ramification on feeder OC protection is provided. Several
existing techniques proposed in the literature to overcome the drawback and
challenges including suppression of the exponentially decay DC-offset in feeder OC
protection are discussed. The discussion is concluded with justification for a new
fault detection and diagnostic tool that is adaptive to the changing distribution feeder
network in the presence of increasing DER penetration as well as having immunity to
the effect of the exponentially decaying DC-offset.
A background review on the characteristics of HIFs and the challenges in their
detection based on the OC protection scheme is discussed. The identifying
characteristics of HIFs and the proposed HIF feature extraction techniques proposed
in the literature are discussed. Moreover, since this research is based on modelling
and simulation, a discussion on the proposed HIF models used in simulation studies
is provided to determine the appropriate model used in this research. The topic
review is concluded with justification of the proposed HIF detection technique as
well as HIF simulation model used in this research.
A comprehensive review on PV systems as well as their integration to the
distribution network feeders is discussed. Furthermore, the types of the DC arc-faults
and their causes discussed. A background review on the existing DC bus protection
scheme in the PV systems and the challenges in DC arc-fault detection is discussed.
18 Chapter 1: Introduction
Then a review on the proposed techniques existing in the literature to overcome the
challenges in DC arc-fault detection is presented. Furthermore, a review on the
proposed DC arc models used in DC arc simulations is done to determine the model
for use in this research. The topic review is concluded with the justification for the
proposed DC arc-fault detection and the DC arc model used in the simulation studies
in this research.
Finally, a summary of the literature reviews and their implications in
developing a research methodology that answers the research questions, and thereby
achieving the research objectives is presented.
Chapter 1: Introduction 19
MFD output is proposed to overcome the challenges in feeder OC protection
influenced by the increasing RE based DER penetration.
The theoretical analysis and simulation result of the technique used in
suppressing the exponentially decaying DC-offset are also presented.
Chapter 7: In this chapter, the conclusions from the research finds are
provided. The significant research contributions are specified, and the benefits and
importance of the proposed methods are summarised. Finally, recommendations for
future research directions are suggested.
20 Chapter 1: Introduction
Chapter 2: Literature Review
2.1 INTRODUCTION
The issues and challenges from the paradigm shift in reshaping the landscape
of the power distribution networks due to increasing levels of penetration of RE
based DERs at the radial distribution feeders are discussed in this chapter. The issues
and challenges including: 1) Challenges in overcurrent protection in radial
distribution feeders with increasing penetration of RE based DERs such as PV
systems, 2) Challenges in HIF detection in distribution networks and 3) Challenges
in DC arc-fault detection on the DC bus in PV systems are discussed.
The chapter begins with discussions on the existing radial distribution feeder
OC current protection including relaying strategy. The drawback in the existing
feeder OC protection system is that, the pickup OC threshold setting is fixed, thus
making it non-adaptive in the face of increasing DER penetration. A review of the
proposed methods in the literature aimed at providing an adaptive solution to the
fixed pickup OC relay setting is discussed. Furthermore, a discussion on the
prevalence of the exponentially decaying DC-offset that occurs during short-circuit
faults and its ramifications on feeder OC protection, and the proposed methods for its
elimination are presented.
Moreover, in power system distribution networks, occurrences of HIF are quite
common. Unlike short-circuit faults, HIFs generally results in low fault current
magnitudes, thus rendering the OC protection system incapable of detecting HIFs.
The methods for HIF detection and classification based on HIF feature extraction
rather than fault current magnitude have spanned decades and are continuing to be an
active part of research today. The highly random and illusive nature of HIF has
proven to be very challenging in developing a universal HIF detection and diagnostic
tool. However, over the years, new methods of HIF modelling and feature extraction
techniques have evolved. A review of the characteristics and challenges in HIF
detection, and the proposed techniques for HIF detection and classification existing
in the literature are presented.
Moreover, the integration of the RE based DERs at the distribution network
feeders have introduced new opportunity for load diversification such as DC loads
The radial distribution network referred to in Figure 1.4 indicates that the
current is directional and emanating from the substation source through the main
radial feeder and distributed to all the users connected to it. The current flow must
always be maintained to guarantee service continuity to all connected users.
Fault conditions such as short-circuit faults occur in distribution network
feeders due to exposure to the elements such as lightning strike, falling trees from
bad weather as well as through unintentional human-caused accidents, etc. Thus, it is
imperative on the service provider to implement a reliable protection system to
minimize damage to network apparatus as well as the possibility of total power
The integration of the DERs requires that the feeder OC protection must be
adaptive to the changing network condition to provide reliable and effective OC
protection. One of the aims of this research defined in section 1.4 is to develop a
strategy for an Adaptive Distribution Feeder OC protection, and the vehicle for such
implementation is the digital OC protective relay.
Time setting
Pickup setting
The protective relays continuously monitor the condition of the power system
by reading voltage and/or current signals, comparing the measured quantity values
with threshold inequality constraint. The relay sends appropriate trip signal
depending on the state of the power system network. The relay determines the state
of the power system network as stated here [1] ,[27]: The operation of the relay to
determine the state of the power system can be graphically illustrated by the block
diagram of Figure 2.3.
Normal State: The power system is said to be operating at its normal state if there
is no sudden or abnormal change (increase/decrease) in the power
system current magnitude. Any normal operating condition that
resembles characteristics of abnormal state must be discriminated
by the relay. Referring to Figure 2.3, lets denote the measured
(input) current x, and constraint to test the measured current against
the pickup setting, Xm. The system is said to be in normal state if x
< Xm.
Abnormal State: The power system enters an abnormal state if x > Xm. This does
not necessarily mean a fault has occurred. It simply implies
abnormal current has been detected. The protection relay gets into
an alert mode to indicate presence of an abnormal condition.
Action State: A time delay, t is initiated and compared against an inequality time
constrain, Tm through the time setting. If conditions, t ≥ Tm, and x
Several different techniques have been proposed for the treatment of the effects
of the DC-offset. The widely proposed technique in the mitigation and elimination of
the DC-offset in the fault current signal require some form of filtering algorithm to
remove the DC-offset as highlighted in [47] and [48]. Most removal algorithms
based on, DFT, Kalman Filter, Morphological Filter, etc., require some form of
parameter estimations, particularly the inductance (L) and resistance (R) parameters
to estimate the Time-Constant, τ of the exponentially decaying DC-offset, the phase
angle and the amplitude. Depending on the algorithm, a full-cycle and two samples
or a half-cycle and one sample is necessary to estimate the required parameters [49],
and if the estimated parameter is matched with the preset value, then a total
elimination of the DC-offset occurs, otherwise the DC-offset can only be suppressed,
but not eliminated. Moreover, the computation for parameter estimation increases the
time delay in extinguishing the fault and can become unacceptable. A notable
contribution in the search for a fast and effective algorithm for the elimination of the
DC-offset is that proposal by Rahmati et al. [50]. This method is DFT-based where
the DC-offset can be reduced to an acceptable level by subtracting the odd and even
samples of the original phasor without estimating any parameter. This algorithm
greatly reduces the time delay associated with all those other algorithms. However
further analysis of this technique needs to be conducted for its real-time application
and robustness.
In summary, the conventional power system distribution network is a passive
network with unidirectional current flow supplying distributed loads along the feeder
Any power system fault conditions or events resulting in low fault current
magnitude not sufficient to overcome the threshold setting of the conventional
overcurrent protection scheme are classified as HIFs. Unlike short-circuit faults,
HIFs normally have high impedance fault path which limits the fault current
magnitude. Occurrences of HIFs are frequently experienced at medium-voltage
(MV) to low-voltage (LV) networks. HIFs can be grouped as, 1) fallen energized
conductor either broken or intact making unwanted electrical contact with high
impedance surfaces such as bitumen, concrete, sand, soil, etc., and 2) energized
overhead conductor making unwanted contact with tree trunk and branches [51].
Such contact surfaces provide a high impedance fault path to ground resulting in a
fault current magnitude of typically between 10-50 A rms with erratic waveform.
While the low HIF current magnitude may not pose immediate significant threat to
power system infrastructure; nonetheless, a fault condition exits. Moreover, if a HIF
is allowed to persist, it increases risk of fire and safety hazard for people. Moreover,
HIFs introduce harmonics, and can degrade the quality of power supply in the long
run therefore must be detected and removed within reasonable time [52].
HIF characteristics are highly random and nonlinear which adds to the
challenges in developing HIF detection and classification using pattern recognition
and other feature extraction techniques. There are several factors including, ground
where Rarc is the arc equivalent resistance, L is the arc length, and Iarc is the arc
current. In practical sense, considering the highly random nature of HIF, this model
fails to capture that. To overcome this limitation, a second category of HIF arc are
proposed where the conductance of the arc is calculated by means of differential
equations [53],[84, 85]. The randomness parameter is introduced with a randomly
varying voltage parameter per arc length. This model is defined by the (2.2) [85];
dg 1 i(t ) − g (t )
=
( )
(2.2)
dt u + R i(t ) l
0
However, it was noted that the single resistor model does not properly account
for randomness in positive and negative half cycle when in contact with a high
impedance surface. Thus, a modification to the original Emanuel arc model is done
by placing either one or two resistor variable resistors in series with the DC voltage
sources to model the randomness introduced by different contact surfaces
[68],[69],[73],[89] . One such model is shown in Figure 2.6. Another HIF arc model
structurally different from the modified Emanuel HIF arc model is shown in Figure
2.7. This model is referred to as the Transient Analysis of Control System (TACS)
model. The TACS model while structurally different, attempts to accurately capture
the characteristics of the HIFs. The TACS model incorporates a nonlinear resistance,
time-varying voltage sources connected in parallel, but only one effective in each
half cycle, and transient analysis of controlled system switch to connect and
disconnect the fault and randomly vary the fault resistance.
The third category of HIF models has found wide application in time domain
simulations of HIF conditions.
Typically, parallel faults are line-to-line short circuit faults that occur between
two points of unequal voltages of opposite polarities. This can be a bridge between
two lines (line-to-line) or two different points on the same line at different voltages
[95] . If a gap exists between the bridging conductors at the point of contact [96] a
DC arc can be ignited creating a route for the fault current through the restrictive arc
channel having resistance defined as arc resistance. This fault is classified as parallel
DC Arc-Fault. Moreover, the arc resistance limits the level of backed current from
the unfaulted PV strings which can be problematic in detecting the fault when fault
current is used as a threshold metric for fault detection. The probability of such faults
occurring is increased by exposed conductors due to insulator breakdown, exposed
dry solder joints, etc., making it possible for accidental bridging of conductors.
Parallel faults with arcing or otherwise occur in parallel to the load.
Unlike parallel faults, series faults are typically open-circuit faults; however, a
high resistive connection is maintained through the arc channel with limited backed
current through the arc resistance. Series arc faults usually result from mechanical
separation of conductor at the solder points or break in the conductor [97],[98] , and
occur in series with the load. Electrical discharge with arcing can occur between the
The research work involved in developing this thesis has involved simulation
of different test conditions. This approach was selected to have flexibility in
analysing different arcing faults under varying environmental conditions such as
changing irradiance. Therefore, to perform DC arc-fault simulation tests as
accurately as possible, developing or selecting an appropriate DC arc model is of
paramount importance.
With DC power systems gaining popularity, the momentum is on the rise to
develop adequate protection scheme for DC systems. Thus, to this effect, proper
modelling of DC arc is essential in analysing DC arc-faults in DC systems including
PV systems. DC Arc-faults when ignited create resistive arc channel with high
resistance that allow enough current to sustain the arc. Thus, there exist a relationship
between voltage and current (V-I) to determine the nonlinearity in DC arc resistance.
Hence, DC arc models have been modelled empirically based on staged tests by
curve fitting technique to determine the nonlinear V-I relationship. A summary of the
common V-I relationships based on the review by Ammarman et al [124] given by
[125], are presented in Table 2.2.
where g is the gap length given in millimetres, and Iarc is the arc current. Using (2.3)
Ammarman [124] developed a model for calculating the DC arc current using the
model for DC rms steady-state arc resistance in terms of the arc current as given in
(2.4)
Rarc =
( 20 + 0.534 g ) (2.4)
0.88
I arc
Vs
where I arc =
Rs + Rarc
C + DL
Ayrton Varc = A + BL + Carbon electrodes
I arc
C ( L + D) Carbon and magnetic
Steinmetz Varc = A + 0.5
I arc electrodes
B
Nottingham Varc = A + n n is related to the electrode
I arc material, L = 0.039 to 0.39 in
a
Paukert Varc = b L = 0.039 to 7.78 in
I arc Iarc = 0.3 to 100 kA
a + cL
Modified Paukert Varc = b + dL L = 0.04 to 0.12 in
I arc Iarc = 3 to 25 A
2.5.1 Summary
In consideration of the dissimilar voltage and current signal attributes in the
domain of AC and DC voltage signals, a computational technique based on time
domain analysis is determined to be most suitable for developing a fault detection
and diagnostic tool that can be seamlessly applied in both AC and DC power
systems. The functional characteristics of MM to detect insignificant variations in the
topography of the graph of the fault signal makes it suitable for detecting changes in
voltage and current signal due to HIFs and DC arc-faults which are both arcing faults
but in the domain of AC and DC signals.
Radial Distribution Feeder Adaptive OC Protection: The integration of the
RE based DERs at the radial distribution network feeders changes the conventional
distribution network from passive network with unidirectional current flow to active
dynamic network with bidirectional current flow. The preceding discussions on
radial distribution feeder OC protection have shown that, increased RE based DER
penetration results in increased current contribution consequently reducing the fault
current level seen by the OC protection relay thus having the potential to
compromise the effectiveness of the OC protection system coordination. Moreover,
during short circuit faults, the fault current is normally affected by the high
frequency transients as well as the prevalence of the exponentially decaying DC-
offset. While the transients can be effectively filtered by means of low pass filtering,
the DC-offset on the other hand is not so easily removed. Its presence can affect the
fault current magnitude, which the feeder OC protection system is contingent on to
mitigate the damaging effects of short circuit faults and any other conditions
resulting in abnormally large current magnitudes.
The existing relay algorithms in power system protection including distribution
feeder OC protection relays are predominantly based on integral transforms such as
Fourier or Wavelet Transforms to estimate the fault current magnitude. The integrity
2.5.2 Implications
The desired research outcome is to develop a universal fault detection and
diagnostic tool that can be utilised in distribution feeder adaptive OC protection as
well as having the ability to detect HIF based on feature extraction. Moreover, the
3.1 INTRODUCTION
The research methodology presented herein defines the overall design of the
research, and the way the research has been conducted to answer research questions
in achieving the research aims and objectives.
In designing the fault detection and diagnostic tool, various features and
functional attributes of the MM signal processing technique have been investigated
- HIF detection based on slope detection aiming two features of the HIFs
including, randomness and arc extinction and re-ignition
iDI G = ( i g ) ( x ) = max i ( x − p ) + g ( p )
p (3.1)
0 ( x − p ) x, p 0
where iDIG is the new finite set formed by the dilation operation on the fault
signal i(x) by the SE g(p). The set iDIG contains the maximum selection of a set of
sums within the neighbourhood of the points around the origin of the SE g. The
process defined by (3.1) involves spatial translation of the signal i(x) by each point of
g(p) then offset by the value of g at point p [151]. Then the dilated signal iDIG is the
pointwise maximum selection of the translated and offset version of the original
signal i(x). The sequence of operation in the dilation transform can be demonstrated
by considering an illustrative example.
Supposing the signal i(x) to be dilated by a SE has 6 points represented on the
Cartesian coordinate with values such that i = [0, 0.1, 0.2, 0.3, 0.2, 0.1]. The SE g
has 3 points with origin at the center, such that the points are p = (-1, 0, 1) with
values of g at these points being g = [0.1, 0, 0.1]. Application of (3.1) involves a
pointwise maximum of three signals.
Let h(p) be a 1 by X matrix to hold the value of i(x) offset by g at each point p.
The result of the addition in the dilation operation on the signal i(x) by values of g at
each point p is held in the h(p) matrix as illustrated below.
The results in the h(p) matrices can be placed in another matrix denoted as H. To
account for the effect of translation by each point of g(p), the signals h(-1), h(0) and h(1)
must be translated as shown in the H matrix. The negative infinity (-∞) is used in
dilation and meaning that the signal is undefined at that point.
h−1 1 X − −
H = − h0 1 X −
− − h1 1 X P X
The signal iDIG is formed by selecting the maximum value in each column in H as
shown below.
iERG = ( i g ) ( x) = min i ( x + p ) − g ( p )
p (3.2)
0 ( x + p ) x, p 0
where iERG is the new finite set formed by the erosion transform using the
minimum selection of a set of differences within the neighbourhood of the points
around the origin of the structuring element, g. The sequence of operation in the
erosion transform can be demonstrated by considering an illustrative example.
The operation of the erosion transform can be demonstrated by considering the
same illustrative example previously used. The signal i(x) used previously is to be
eroded by the same SE, g. Application of (3.2) involves a pointwise minimum of
three signals.
The h(p) matrices and their contents are defined for the erosion process herein.
The result of the subtraction in the erosion operation on the signal i(x) by values of g
at each point p is held in the h(p) matrix as illustrated below.
The results in the h(p) matrices can be placed in H as previously defined. To account
for the effect of translation by each point of g(p), the signals h(-1), h(0) and h(1) must be
translated as shown in the H matrix. The positive infinity (∞) is used in erosion to
represent points where signal is undefined.
h−1 1 X
H = h0 1 X
h1 1 X P X
The eroded value of i(x) results from the minimum selection of each column in the H
matrix. Let iERG be a 1 by X finite set containing the minimum selection of columns
of H. The minimum selection of the valid columns of the H matrix is shown below.
The reclassified iERG considering the effect of translating the SE is as given below.
( i g ) = ( i g ) g = max
p (min i ( x + p) − g( p)) ( x − p ) + g( p)
p
(3.3)
As can be seen in (3.3), opening transform is the dilation of the erosion of the
signal i by the same SE g, and in (3.4), the closing transform is the erosion of the
dilation of the signal i by the structuring element g. The opening and closing
transforms are also dual transforms, where the opening action results in the reduction
of small positive regions within the signal, whilst the closing action will result in the
reduction of small negative regions within the signal. Furthermore, the opening and
closing transforms allow for the recovery of features of the signal lost due to the
erosion and dilation operations respectively, thus providing the characteristics
necessary for feature extraction [133].
The physical effect of opening and closing transforms are graphically
illustrated in Figures 3.2(a) and (b) respectively by considering the same fault current
i(x) and the structuring element g(p). The opening and closing transforms provide the
building block for the construction of two fundamental nonlinear morphological
filters called, Closing-opening and Opening-closing filters. These filters when
connected in cascade and operated simultaneously with the same SE can be very
effective in eliminating positive and negative impulse noise in the signal. These
properties of mathematical morphological transforms and filters can be exploited by
uniquely combining them to develop morphological fault detection and identification
algorithm which is fast, robust and reliable.
Figure 3.2: Physical effect of (a) opening and (b) closing transforms
iDI G + iERG ( i g )( x ) + ( i g )( x )
AvMMF = = (3.5)
2 2
OC ASF = ((( f )
g1 ) • g1 ) g 2 • g 2 (3.6)
The fault detection and classification technique presented in this thesis proposes
to use voltage and current signals measured at the distribution feeder substation. The
measured voltage and current signals provide inputs to the HIF detection and
classification algorithm centrally located at the substation. The occurrence/existence
of faults is detected by continuously monitoring and analysing the extracted voltage
and current signals for changes in current and/or voltage magnitude.
In the proposed algorithm, two MM filters, namely MMF and ASF are
strategically combined, and with the use of the designed SE, create a MFD signal
output from the input fault signal quantities for fault detection. The mathematical
derivations for the entire signal processing technique in creating the MFD output is
described herein.
3.4.2 Sampling
A sampling frequency of 1.92 kHz suitable for real-time application [78] in the
detection of power system disturbances is applied in sampling the signal. It has been
further suggested that two times recommended sampling frequency could be used.
Thus, considering 64 samples per fundamental cycle in the design, a sampling
frequency of fs= 3200 Hz for 50 Hz system is used. All designed parameters
presented herein and thereafter will refer to 50 Hz system; however, values can be
recalculated by interchanging the frequency for 60 Hz system. Considering the
number of samples, and using (3.9) [153] , the sampling rate used in the design is
312.5 µs.
The signal can be written as sampled signal for a fundamental cycle as given in
(3.11);
f ( n ) = F0 r n + F1 sin n + − n = 0,....., N − 1 (3.11)
N
where r = e-(∆T/τ), and the first and second terms are the DC and the AC components
respectively [153]
W −1
f w (t ) = f (t ) (t ) w(t ) = f (wT ) (t − wT ) (3.12)
w=0
where fw(t) is the fault signal window. The wth sample in the data window can be
represented (3.13);
The first and the last sample in the data window are at position k = 0 and k = N-1.
These samples define the beginning and the end of the data window. Each data
window is separated by ∆T. The data window can be written as in (3.14);
f ( n ) = F0 r n + F1 sin n + − n = 0,....., N − 1 (3.14)
N
Therefore, let F{n1 n2 n3…….nN} represent a window in the domain of the fault signal
f(n) to be processed by the MFD algorithm at each natural time step.
g (m) = cos m cos(m − 1).... cos 1 cos ..... cos(m − 1) cos m (3.15)
where m is the number of points in the SE; ϕ is the phase angle (ϕ = 2ᴫf0∆T); f0
is fundamental frequency and ∆T is the sampling interval. Note that, points of the
grayscale SE denoted p in (3.1) and (3.2) and all MM operation thereafter is replaced
by m in (3.15), and the reason will become obvious in the proceeding section.
Decomposition of Structuring Element: While the general SE as given in
(3.15) could be used in its entirety, it is not computationally efficient to do so as the
same result can still be archived with a smaller SE with less number of points. In the
overall design of the algorithm, a sampling frequency of 3.2 kHz to get 64 samples at
a sampling rate of ∆T = 312.5μs per fundamental in 50 Hz was used. Considering
that the 64 samples are segmented into data windows of 16 samples, the SE can have
points (or size) up to 16. However, more sampling points will increase computational
MMF-Stage 1: The dilation and erosion of original fault signal for a data
window segment, f(n) by A1 are given by (3.19) and (3.20) respectively.
DI1 = ( f A1 ) = max f ( n − a1, m ) / a1, m m = 1, 2, 3
m
(3.19)
( )
DI 2 = ( f A1 ) A2 = max DI1 n − a 2,m / a 2,m m = 1,2,3
m
(3.21)
( )
ER2 = ( f A1 ) A2 = min ER1 n + a2, m / a1, m m = 1,2,3
m
(3.22)
where DI1 and DI2 contain maximum selections of sets of sums from dilation
operation while ER1 and ER2 contain the minimum selections of set of differences
from erosion operation within the neighbourhood A1 and A2 in stage 1 and 2 of the
decomposed MMF respectively. The processed signal can be reconstructed by taking
the average output of the decomposed MMF as defined by (3.23);
where (DI1 + ER1)/2 and (DI2 + ER2)/2 are the average outputs of MMF-Stage 1 and
MMF-Stage 2 respectively. Equation (3.23) represents the magnitude transformation
of the original fault signal. Let fmmf (n) denote the transformed version of f(n) given
by (3.24);
f mmf ( n ) = F0 mmf r n + F1mmf sin n + − n = 0,....., N − 1 (3.24)
N
where F0mmf and F1mmf are the dc and ac components of the MMF average output
respectively.
The DOCAS algorithm is based on slope detection to generate spikes in
response to variation in the topography of the fault signal waveform caused by the
fault. The distortion of the signal waveform from the transients following a fault are
emphasised by subtracting fmmf (n), (3.24) from the original fault input signal, f(n)
(3.16) to produce the initial fault detection signal given in (3.25);
N −1
f (n) = F0 r n + F1 sin n + (3.25)
n =1 N
OCASF =
(((( f ) )
B2 ( ) • B2 ( ) ) B1 (u) • B1 = 1,...,5 and u = 1,.,3
(3.26)
where OCASF is the entire sequential transformation covering the four stages. A stage
is defined by a closing (●) and (○) operator in each ASF layer. Thus, each stage of
the Open-close ASF layer can be elaborated as follows;
OCASF-Stage 1:
OC1 (n) = ( f B2 ) ( )
(
(3.27)
= max
)
min f ( n + ) / B2 ( ) ( n − ) / B2 ( ) = 1,..,5
OCASF-Stage 2:
OC2 (n) = ( ( f B2 ) • B2 ) ( )
(3.28)
= min
( max OC ( n − ) / B ( )) ( n + ) / B ( ) = 1,..,5
1 2 2
OCASF-Stage 3:
OC3 ( n) = ( ( ( f )
B2 ) • B2 ) ( ) B1 (u )
(3.29)
= max
( min OC ( n + u ) / B (u)) ( n − u ) / B (u) u = 1,., 3
2 1 1
OCASF-Stage 4:
COASF =
(((( f • B ( ) 2
) )
B2 ( ) ) • B1 (u) B1 = 1,...,5 and u = 1,.,3 (3.31)
where COASF is the entire sequential transformation covering the four stages. A
stage has been previously defined is denoted by a closing (●) and (○) operator. The
signal transformation at each stage is thus described as follows;
COASF-Stage 1:
CO1 (n) = ( f • B2 ) ( )
(3.32)
= min
( )
max f ( n − ) / B2 ( ) ( n + ) / B2 ( ) = 1,..,5
COASF-Stage 2:
CO2 (n) = ( ( f • B2 ) B2 ) ( )
(3.33)
= max
( min CO ( n + ) / B ( )) ( n − ) / B ( ) = 1,..,5
1 2 2
COASF-Stage 3:
CO3 (n) = ( (( f • B ) 2 B2 ) ( ) • B1 (u ) )
(3.34)
(
= min max CO2 ( n − u ) / B1 (u ) ( n + u ) / B1 (u ) u = 1,.,3
)
COASF-Stage 4:
CO4 (n) = ( ( (( f • B ) 2 B2 ) ( ) • B1 ) )
B1 (u )
(3.35)
= max
( min CO ( n + u ) / B (u)) ( n − u ) / B (u) u = 1,., 3
3 1 1
The entire signal transformation process to create the MFD output can be
graphically illustrated by considering a simple AC signal. Moreover, the attributes of
the MFD output and its application fault detection and classification are presented
herein.
The technique applied in this algorithm makes use of one unique characteristic
of the MM in edge detection [159], [160]. This can be observed on the average
MMF output of Figure 3.4 (b), the leading and the lagging edges of the data window
Figure 3.4: DOCAS Response at the MMF stages (a) Simple AC input signal, (b) MMF average
output and (c) Difference fault signal, ∆f
The passage of the initial fault detection signal, Δf through the DOCAS
algorithm at the ASF stages, and the extraction and combination of signal at each
stage of the ASF layers results in the MFD signal output of Figure 3.5.
The MFD output signal is created through the process described by (3.25) –
(3.34). In Figure 3.5, the MFD output signal is segmented into fault detection
windows accomplished through the edge detection characteristic of the MMF. The
process of summing the differences of closing and opening transformation at each
ASF stages results in the suppression of smaller spikes within the fault window and
giving visibility to the edge spikes. The MFD fault windows are fundamental to the
fault detection functionality of the DOCAS algorithm. Based on the design
consideration of 64 samples per fundamental period and data window size of 16
samples, a fault detection window is of a quarter cycle duration resulting four fault
detection windows per fundamental period. A fault detection window is defined by
two edge spikes, one beginning and one ending. Moreover, it can be observed that
2kNT
MFDTall = k = 0,.., − 1 (3.37)
2k (N + 1)T
kNT
MFDShort = k = 1,.., (3.38)
k (N + 1)T
Figure 3.5: Fault detection windows of the DOCAS MFD output signal
where Iinc(pu) is the per unit increase in the current magnitude. Generally, in
overcurrent protection, increase in current magnitude is sensed for abnormalities or
short circuit fault conditions. Thus, if Iinc ≤ 1, no abnormality exists, or no
disturbance has occurred, and the system is in normal state. If Iinc > 1, then
disturbance has occurred. The detection of any disturbance does not necessarily
mean a fault exists. To possibility of a fault exist though and must be tested against
overcurrent threshold parameter to determine if the disturbance is a fault or a normal
operating condition.
Figure 3.7: DOCAS response to SLG fault (a) Fault voltage waveforms for SLG fault on
phase A, (b)-(d) corresponding MFD outputs for each phase.
Figure 3.8: The MFD tall edge spikes for current, (a) MFDTall (1) and (b) MFDTall (2) The
MFD tall edge spikes for current, (a) MFDTall (1) and (b) MFDTall (2)
The DOCAS algorithm operates on the fault signal inputs to extract features for
HIF detection and classification. The two HIF characteristics targeted by the DOCAS
algorithm are 1) arc extinction and re-ignition of the HIF arc around the fundamental
period resulting in shoulder shaped unsymmetrical fault current waveform and 2)
randomness, due to randomly changing effective fault resistance, Rf resulting in
erratic fault current waveform. The DOCAS algorithm generates MFD spikes at
specific locations on the output MFD signal that correlate to the target HIF
characteristics. The MFD fault detection windows are fundamental to this concept
where changes in the edge spikes as well as appearance of random spikes within the
fault detection window correlate to transients generated by a disturbance or an event.
The HIF current flowing through the high impedance path represented R f
normally result is very low RMS current magnitude of between 10A-50A. While this
will not have any noticeable or significant impact on the fault current seen by the
feeder OC relay, it is however assumed that the HIF current will leave its footprint
on the feeder current and voltage signals. The DOCAS algorithm can detect
insignificant changes on the topography of the input signals and generates MFD
spikes at its output that correlate to these changes. While both voltage and current
signals can be considered for HIF feature extraction, in the DOCAS algorithm the
voltage signal is utilised as it was observed that voltage signal presented better
transient characteristics compared to current signal.
HIF Arc Extinction and Re-ignition Characteristic: The current and voltage
waveform intersect once every half cycle. If they are in phase, then this happens at
zero magnitude of each signal. If they are not in phase, then at the point of
intersection, the instantaneous voltage and current magnitude are equal. In a 50 Hz
The AC arc extinguishes and re-ignites twice every have cycle. The arc
remains extinguished for a short period of time while waiting for the voltage to
regain until it reaches the restriking level to overcome the gap separation break down
voltage for arcing to happen again. This phenomenon impacts the fault voltage at
specific and fixed period defined in (3.40). The arc extinction and re-ignition thus
causes the MFD spikes defined as MFDArc to grow in height when HIF occurs.
Hence, the MFDArc spikes can be separated to observe the arc extinction and re-
ignition HIF characteristic. The MFDArc spikes maintain uniform height under
normal condition as observed in Figure 3.10(c) for the period before 0.03s. A HIF is
initiated at 0.03s and the HIF arc extinction and re-ignition occurs about zero-
crossing of the fundamental period shown in Figure 3.10(b) resulting in increased
height of the MFDARC spikes above the prefault level as shown in Figure 3.10(c) for
the period after 0.03s.
HIF Randomness Characteristic: The randomness characteristic would be
captured at the edge spikes and fault windows as previously defined (refer Figure
3.5). The tall and short edge spikes, MFDTall and MFDShort respectively defined in
(3.37) and (3.38) randomly vary in height under HIF condition. Transients occurring
closer to the zero-crossing of the fundamental cycle will cause random variation in
the height of MFDTall spikes while those closed to the positive and negative peaks
will cause random variation in the height of the MFDShort Spikes.
The appearance of the elongated edges spikes is non-uniform and random, and
can occur at any point from k = 0 to ∞-1 and k = 1,..,∞ as defined by (3.37) and
(3.38) if HIF exists. Those transients occurring away from the zero-crossing and the
positive and negative peaks will cause random spikes to grow inside the fault
windows. These spikes defined as MFDWindow will appear inside the fault windows
and will persist in the presence of HIF.
The DOCAS algorithm is a multistage MM filter designed for the detection and
diagnosis of all types of faults in radial distribution feeders and DC arc-faults on the
DC bus of DC power systems including PV systems. The structure of the DOCAS
algorithm and its mathematical derivations are universally applied in all respect.
However, the input signals are different where its application in fault detection and
diagnosis in the radial distribution feeder relies on AC voltage and current signals as
input quantities while in DC arc-fault detection, the input signals are DC voltage and
current. Moreover, the SEs, A1, A2, B1 and B2 derived from the eccentrically
decreasing convex SE are utilized in the DOCAS algorithm in DC arc-fault
detection. While the SE used is designed from attributes of the AC signal, its reduced
structure for computational efficiency giving five points eccentrically decreasing
weighted convex having the following points [0.957, 0.995, 1, 0.995, 0.957] with
two unequal slopes is adopted without any change for DC arc-fault detection. The SE
as previously applied in AC system fault detection are decomposed into, A1 = [0.995,
1, 0.995], A2 = [0.957, 1, 0.957], B1 = [0.995, 1, 0.995] and B2 = [0.957, 0.995, 1,
0.995, 0.957]. A1 and A2 are used in the MMF section while B1 and B2 are used in the
ASF section.
The DOCAS algorithm as a DC Arc-Fault Detector retains its original design,
however its response to AC and DC signal inputs are significantly different. The
The difference DC current and voltage signals are transformed by the two layers
ASF, and at each stage of the ASF layers, the transformed signals are extracted and
combined to form the MFD outputs according to (3.36), and the results are shown in
Figures 3.11(d) and 3.12(d) for the DC current and voltage respectively.
It is obvious from Figures 3.11 and 3.12 that the concept of fault detection
window is not available because of the non-sinusoidal nature of the (DC) signal. The
onslaught of the arc is detected by appearance of spikes which increase in height in
response to the fast changing (rate of change of the random DC arc) noise like DC
arc characteristic that sustains when ignited. The spikes in the MFD output appear
chaotically. The DOCAS algorithm detects DC arc-fault by means of detecting the
chaotic behaviour of the random DC arc. Moreover, the MFD output only produces
spikes when arcing occurs.
Figure 3.12: DOCAS response to DC arc-fault in PV systems, (a) DC arc-fault current, (b)
Average MMF output, (c) diff DC fault current, ΔI and (c) MFD output.
3.6 CONCLUSION
A multistage fault detection and diagnostic tool called the DOCAS algorithm
utilising two classes of nonlinear morphological filters called the MMF and ASF was
developed. A SE with a sinusoidal geometrical structure was designed based on the
sampling rate and the fundamental frequency with points as many as 16 points
representing the data window. The SE was restructured in the interest of maintaining
computational efficiency to five points without compromising the integrity of the
algorithm. The restructured algorithm resulted in an eccentrically decreasing SE with
two unequal slopes. Further restructuring was made by decomposing the SE into two
4.1 INTRODUCTION
84 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
feeder OC protection. In section 4.8, a method for adaptive inverse time OC relaying
using the DOCAS algorithm is presented followed by conclusion in section 4.9.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 85
𝑃𝑉(𝑢)∀𝑢 = 1. . , 𝑈. The voltage and current measured at the point of coupling of the
uth PV can be defined as 𝑉𝑃𝑉 (𝑢) and 𝐼𝑃𝑉 (𝑢) respectively. Considering that the PV
sources are assumed to inject only positive sequence voltage and current under load
condition, then the positive-sequence voltage and current measured at PV(u) can be
defined as 𝑉 (1) 𝑃𝑉 (𝑢) and 𝐼 (1) 𝑃𝑉 (𝑢) respectively. Then an equation that relates the
per phase (positive-sequence) prefault voltage and current for the uth PV source can
be expressed in (4.1) [161]
E (u ) = VPV
(1)
(u ) + Z PV
(1) (1)
(u ) I PV (u ) (4.1)
(1)
where 𝑍𝑃𝑉 (𝑢) is the series positive-sequence Thevenin equivalent impedance at
PV(u) for all interfacing devices including the transformer. Supposing a fault occurs
at some point on the feeder, the equation in (4.1) can be modified to reflect the fault
condition as given in (4.2);
E(u) = VF(1)
PV
(u) + Z PV
(1)
(u) I F(1)PV (u) (4.2)
(1) (1)
where 𝑉𝐹𝑃𝑉 (𝑢) and 𝐼𝐹𝑃𝑉 (𝑢) are the positive-sequence fault voltage and current
respectively seen at PV(u) injected by the uth PV source.
VPV(1)
(1)
Z PV (u ) = (u ) (4.3)
I PV
(1)
The positive sequence voltage and current measured by the DFRs at the PCC
of the PV sources are phase values. The zero- and negative sequence impedance can
be obtained following the standard convention as indicated by Figure 4.2. The zero-
86 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
and negative-sequence impedances at PV(u) can be obtained by transforming the
voltage to sequence components as follows;
Vs (u ) = AVp (u ) (4.4)
I s (u ) = AI p (u ) (4.5)
where Vs and Is are sequence voltage and current respectively, and similarly Vp and Ip
are phase voltage and current respectively. A is the transformation matrix such that,
1 1 1
A = 1 a 2 a (4.6)
1 a a 2
2𝜋
where 𝑎 = 𝑒 𝑗 3 . Then using standard notations of sequence networks shown in
Figure 4.2, the negative- and zero-sequence Thevenin impedances for the uth PV
source are;
VF(0)
Z (0)
PV (u ) = − PV
(u ) (4.7)
I F(0)
PV
VF(2)
(2)
Z PV (u ) = − PV
(u ) (4.8)
I F(2)
PV
All inertial sources including the substation source are modelled similarly using the
voltage and current measurements at their point of coupling. Replacing all subscripts
in (4.1) - (4.2) from PV to SS, the positive-, negative- and zero-sequence Thevenin
impedances of the substation (SS) source as seen from the point of measurement are;
(1) (2) (0)
𝑍𝑆𝑆 , 𝑍𝑆𝑆 and 𝑍𝑆𝑆 respectively.
0
𝑍𝑃𝑉 1
(𝑢) 𝑍𝑃𝑉 2
(𝑢) 𝑍𝑃𝑉 (𝑢)
1 2 0
𝐼𝑃𝑉 𝐼𝑃𝑉 𝐼𝑃𝑉
1 2 0
E(m) 𝑉𝑃𝑉 𝑉𝑃𝑉 𝑉𝑃𝑉
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 87
on the type of the fault at the point of common coupling with the sequence networks
of the distribution feeder.
Considering the system in Figure 4.1, assume that there are U number of PV
systems connected at locations along the feeder length such that, the distance to the
point of coupling of the uth PV from the substation is denoted dPV(u). Supposing the
total feeder distance is l, then a distance factor, denoted df(u)∀ u=1,..,U can be
defined by (4.9).
dPV (u )
d f (u ) = u = 1,..., U (4.9)
l
PV1 is the nearest PV system to the substation with a distance of dPV(1). It is
assumed in this analysis, that PV1 is connected at a close proximity to the feeder
substation, thus is assumed to have a zero meter (dPV(1) ≈ 0 ) from the substation .
The furthest PV system is PV(U), where dPV(u) = l, then df can have possible values
between 0 and 1.
Supposing a fault occurs at some point on the distribution feeder of Figure 4.1.
Similarly, the distance to fault as seen from the substation can also be defined by a
fault distance factor denoted dff. Supposing the distance to fault is denoted dF in
meters; then the fault distance factor can be defined as;
dF
d ff = (4.12)
l
dff can have values between dffmin = 0 and dffmax = 1, where fault closes to the
substation, assuming at PV(1) will have dffmin value while a fault at the end of the
feeder, assuming at PV(U) will have dffmax value.
88 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
to the uth PV system. The factor ε can have absolute value between 0 and 1, such that
| ε(u) | = 0 indicates fault location is at PV(u) and | ε(u) | = 1 indicates fault location
either at the feeder substation or at the end of the feeder length with respect to PV(U)
and PV(1) respectively. The smaller the | ε(u) | is, the closer the fault is to PV(u).
Fault can be located at either downstream (on the right) or upstream (on the left) of
PV(u) with respect to the direction of the current flow. The direction of the fault
location with respect to PV(u) depends on the polarity of ε(u). If ε(u) is a positive
value the fault is downstream, and if ε(u) is a negative value the fault is upstream.
U
I FPVup =
u =U − D
I FPV ( u ) (4.14)
The total fault current injected downstream to the point of fault is;
D
I F1 = I Fss + I FPV ( u ) (4.15)
u =1
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 89
where IFss is the fault current contribution from the substation source which is
additive to the fault current injection by all downstream PV sources.
Assuming 𝐼𝐿 ≈ 0, the total fault current injected upstream to the point of fault
is;
U
I F2 =
u =U − D
I FPV ( u ) (4.16)
D U
I F = I Fss + I FPV ( u ) + I FPV ( u ) (4.17)
u =1 u =U − D
Noting that 𝐼𝐹𝑆𝑆 is the fault current contribution by the substation source, it is
obvious from (4.17) that high penetration of PV systems will increase the fault
current contribution by them, and consequently reduce the fault current magnitude
seen at the feeder substation. Moreover, 𝐼𝐹𝑆𝑆 is a component of 𝐼𝐹1 which is the total
fault current injection from all sources, including the substation source feeding the
fault downstream. Hence, a large number of downstream sources with respect to the
fault location will also reduce the fault current level seen at the feeder substation.
This implies that, fault closest to the substation will register higher fault current
contribution from the substation source as there will be less downstream PV sources.
As the fault moves further away, the fault current contribution by the substation
source will reduce as the PV sources feed the fault. However, with high penetration
of PV DERs with large capacitor, the fault current magnitude seen at the feeder
substation can be significantly reduced irrespective of the fault location thus
affecting the over current system protection system coordination.
Figure 4.3: Per phase circuit diagram of the typical radial feeder system with PV
90 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
4.4 EFFECT OF FAULT LOCATION ON FEEDER SUBSTATION FAULT
CURRENT LEVEL
1
'
ZThdn = (4.18)
1 1
+
1
+ Z L (3) Z PV (3)
1 1
+
1 Z PV (2)
+ Z L (2)
Z PV (1)
1
'
ZThup = (4.19)
1
1 1
+
1 Z PV (4)
+ Z L (5)
Z PV (5)
where ZPV(u) is the Thevenin impedance of the PV source seen at the PCC and
ZL(u) is the impedance of the line segments between PV sources along the feeder.
Supposing the impedance per meter of the feeder line can be represented by z Ω/m,
then the total feeder line impedance is Z = zl. The feeder line impedance between
PV(u) and the substation, in terms of the distance factor is then,
Z L (u ) = d f (u ) Z Ω (4.10)
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 91
Then, the total impedance of the feeder line in terms of the line segments from the
PV(1) to PV(U) can be expressed by (4.11),
U
Z = ( Z L (u ) − Z L (u − 1) ) Ω (4.11)
u =2
The impedance of the line segment between PV sources in terms of the distance
factor, df(u) is, ZL(u) = Z(df(u) – df(u – 1)). It must be noted that, the Thevenin
impedances defined in (4.18) and (4.19) do not include the impedance of the faulted
line segment. Considering that the fault occurs in the line segment between PV(3)
and PV(4), the impedance of the fault line segment is then ZL(4) such that in terms of
PV distance factors, ZL(4) = Z(df(4) – df(3)). Fault can be located at any point within
ZL(4). The distance to fault with respect to the substation is represented by the
distance factor dff. The proximity to the fault with respect to PV(3) is, ε(u) = dff –
df(u) |u = D = 3, then the impedance of the faulted line segment with respect to PV(3) is
ZL(4) ε(u)|u = D = 3. Conversely, the proximity to the fault with respect to PV(4) is, ε(u)
= dff – df(u) |u = ( D+1) = 4. However, df(4) > dff which means fault is upstream with
respect to PV(4) then the impedance of the fault line segment with respect to PV(4)
is ZL(4) (1- ε(u))|u = D = 3 . These impedances are added in series to ZThdn’ and ZThup’
respectively. The following relationship can be deduced for the downstream sources
by applying KVL to the downstream sources in Figure 4.4;
Vss = I Fss ( Zss + Z L (4) (u) ) + I FPVdn Z L (4) (u) |u =D=3 (4.20)
(
VPVdn = I Fss Z L (4) (u ) + I FPVdn ZThdn
'
)
+ Z L (4) (u ) |u = D =3 (4.21)
where Vss is the substation source Thevenin voltage, and VPVdn is the Thevenin
voltage of all the downstream sources. The Thevenin equivalent impedance for the
downstream side of the feeder under fault seen from the point of fault is;
'
Z ss ZThdn
ZThdn = + Z L (4) (u ) |u = D =3 (4.22)
Z ss + ZThdn
'
Let VThdn be the Thevenin equivalent voltage for all the downstream sources
including the substation source, then the fault current contribution from all the
downstream sources including the substation source is;
92 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
VThdn
I F1 = (4.23)
ZThdn
Let VPVup be the Thevenin voltage of all the upstream PV sources, then the upstream
sources can be similarly expressed as given in (4.24)
(
VPVup = I FPVup ZThup
'
+ Z L (4) (1 − (u ) ) |u = D =3 ) (4.24)
The Thevenin equivalent impedance for the upstream side of the feeder under fault
seen from the point of fault is;
ZThup = ZThup
'
+ Z L (4) (1 − (u) ) |u =D=3 (4.25)
Let VThup = VPVup be the Thevenin equivalent voltage for all the upstream PV sources,
the fault current contribution from all the upstream PV sources is;
VThup
IF 2 = (4.26)
ZThup
Assuming a three-phase fault such that ZF ≠ 0, and applying KVL, the following
relationships can be derived;
VThdn ZThdn
'
+ Z L (4) (3) + Z F ZF I F1
V =
+ Z L (4) (1 − (3) ) + Z F I F 2
' (4.27)
Thup ZF Z Thup
ZThF =
(Z '
Thdn + Z L (4) (3) + Z F ) ( ZThup
'
+ Z L (4) (1 − (3) ) + Z F )
(4.28)
(Z '
Thdn + Z L (4) (3) + Z F ) + ( ZThup
'
+ Z L (4) (1 − (3) ) + Z F )
In the case of bolted fault, ZF is removed from (4.28). Moreover, the sequence
components are appropriately connected subject to the type of fault.
Supposing VThF is the Thevenin equivalent voltage seen at the point of fault,
then the fault current at the point of fault is;
VThF
IF = (4.29)
ZThF
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 93
Considering equation (4.28), the Thevenin impedance seen at the point of fault F is a
function of fault distance and level of PV penetration. Increased PV penetration will
result in reduced Thevenin impedance resulting in more fault current contribution by
the PV sources consequently lowering the fault current magnitude seen at the feeder
substation.
Figure 4.4: Per phase circuit diagram categorising PV sources into downstream and
upstream sources with respect to point of fault.
In sections 4.3, 4.4 and 4.5, theoretical network analysis was done to
demonstrate the effect of PV penetration on fault current magnitude seen by the main
OC protection relay at the feeder substation. In this section, the DOCAS algorithm is
used to analyze short circuit faults, then to propose a method for making the feeder
OC protection relay threshold setting adaptive to increasing PV penetration.
Figure.4.5 illustrates the structure of the proposed scheme incorporating the
DOCAS algorithm. Mathematical derivations and the functions of the major blocks
excluding the DOCAS algorithm are presented herein. Voltage and current signals
measured at the feeder substation constitute inputs to the algorithm. Discussions
considers current signal, however is equally applicable to voltage signal as well.
94 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Figure 4.5: The OC fault detection and diagnostic scheme incorporating the DOCAS
algorithm.
Power system signals (voltage and current) are normal distorted during a fault
condition. These signals usually contain 1) Decaying DC component, 2)
Fundamental Frequency Component, and 3) Integer harmonics as well non-integer
harmonic components. The presence of the decaying dc-offset and the harmonics
causes distortion on the fundamental frequency component, and this can affect the
accuracy of the algorithm. While the harmonic components can be effectively
filtered by appropriate low pass filters, the decaying dc-offset on the other hand is
not so straight forward. The presence of the exponentially decaying dc-offset is of
concern especially in estimating the magnitude of the fault current and/or voltage
signal.
DC Offset Removal in the Proposed Overcurrent Protection Method: The
removal of the exponentially decaying DC-offset is accomplished through the natural
process of creating the initial fault detection signal, Δf, such that Δf = ∆i is the
difference fault current signal. The method does not require any parameter estimation
and is the integral part of the fault detection process. The technique implemented in
this algorithm is based on adaption of the technique proposed by A.Rahmati et al
[50]. However, in this case two signals, the original fault signal I and the average out
of the MMF filter, Immf are subtracted to eliminate the DC-offset as opposed to the
subtraction of the even and odd samples implemented by [50]. The mathematical
derivation given is intended to describe removal of the DC-offset through the
subtraction of the signals. It does not imply the implementation of any DFT based
filter for this purpose.
Based on the equation for the discretised signal in (3.11), a discretised input
current signal is expressed in (4.30);
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 95
H
i(nT ) = I0 r n +
h =1
I h sin ( 2 f h nT + h ) (4.30)
where r = e-ΔT/τ such that I0rn is the DC component with a magnitude and time
constant represented by I0 and τ respectively. The second part of (4.30) is the AC
component with Ih, fh and ϕh being the amplitude, frequency and phase angle of the
hth order harmonic with H representing the highest order harmonic in the fault
current. The phasors for the AC and DC parts of the fundamental frequency (h = 1)
component can be calculated by applying DFT as follows;
AC Part:
H
i(nT ) AC = I
h =1
h sin ( 2 f h nT + h ) (4.31)
1 2 1 2
I ac (n) = I1 1 + cos n + 1 + jI1 1 − sin n + 1 (4.32)
N N N N
1 2
Re I ac (n) = I1 1 + cos n + 1 (4.33)
N N
1 2
Im I ac ( n) = I1 1 − sin n + 1 (4.34)
N N
( Re I ac (n)) + ( Im I ac (n))
2 2
I ac (n) = (4.35)
DC Part:
nT
−
i ( nT ) DC = I 0 r = I 0 e
n (4.36)
96 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
I DC (n) = Mr n (4.37)
where,
2 2
r cos − 1 + j r sin N
M =
N
2I0 r − N (4.38)
2 2
N 2r cos − 1 − r − 1
N
2 n
2 I 0 r − N r cos − 1 r
N
Re I dc (n) = (4.39)
2 2
N 2r cos − 1 r − 1
N
2 n
2 I 0 r − N r sin r
N
Im I dc (n) = (4.40)
2 2
N 2r cos − 1 r − 1
N
The fundamental frequency component of the fault signal in (4.30) can be expressed
in terms of the real and imaginary parts of the AC and DC parts are follows;
Then the amplitude of the fundamental frequency component of the fault current
with exponentially decaying DC offset is;
The fault current signal transformation at the decomposed MMF stages and its
reconstruction resulted in the average MMF output signal given (3.22). This signal
can be expressed as given in (4.42)
H
immf (nT ) = I 0_ mmf r + n
I h _ mmf (
sin 2 f h _ mmf nT + h _ mmf ) (4.42)
h =1
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 97
where I0mmfrn is the DC component of the average MMF output and the second part
of (4.42) is the AC part with Ih_mmf, fh_mmf, ϕh_mmf are the amplitude, hth order harmonic
frequency and phases angle. It must be noted that, at the MMF stages, the fault signal
results in only magnitude transformation while it is assumed that all other signal
attribute, for the DC and the AC part remain unchanged.
Applying DFT to the AC and DC parts of the average MMF output will give
the amplitude of the fundamental frequency component of the average MMF output
with exponentially decaying DC offset as;
It must be noted that, the MMF stages emphasizes the samples at the edges of the
data window on the average MMF output. Thus, the AC signal component is also
reduced except at the edges of the data windows. The peak value of the fundamental
frequency component of Δi without the exponentially decaying DC-offset are relative
to the height of the tall edge spikes (ref Figure 3.10(c)) appearing at the regular
intervals given by (4.44).
i(2kN T )
k = 1,...., (4.44)
i(2k ( N + 1)T )
98 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
appropriate low pass filter must be designed to reduce the high frequency ripples. In
this filter design, selection of any cut-off frequency is not necessary. However, a
digital filter that is time invariant and does spectral shaping, and selective frequency
filtering without assigning any cut-off frequency is desirable. An averaging filter is
one such filter, and the one used in this design is the exponentially weighted moving
average (EWMA) filter [162],[163] . The signal needs smoothing is Δi.
Given that Δi is a sequence of sampled data, lets denote the sequence as Δis.
Supposing at any instant, t = nΔT, the average of the previous samples, denoted p
can be computed as follows, [163];
t
1
i (pt ) = is (4.45)
p s =t − p +1
Let’s consider the point t+1 = (n+1)ΔT where the average of the previous samples of
the data sequence can be computed as [163].
1
t +1 t
1
i (pt +1) = is = it +1 + is (4.46)
p + 1 s =t − p +1 p +1
s = k − p +1
t
Let’s denote
s =t − p +1
is = pi (pt ) , then (4.46) becomes;
1 1 p (t )
i (pt +1) = it +1 + pi (pt ) = it +1 + i p (4.47)
p +1
p +1 p +1
The expression for the exponentially weighted moving average filter can be obtained
by shifting the time index back one step. Thus, the expression for Δip(t) is;
1 p (t −1)
i (pt ) = it + i p (4.48)
p +1 p +1
p 1
Let’s denote, = , then = (1 − ) , such that the exponentially weighted
P +1 P +1
moving average (EWMA) filter equation is [163].
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 99
The parameter, α is a weighting constant that dictates the degree of filtering, where it
can have any value within the range 0 ≤ α < 1. Moreover, in (4.49), the EWMA Δip(t)
denotes a weighted average of all past data in the data sequence Δis at point t. When
α → 0, a large number of past data points are considered in computing the EWMA
Δip(t), and as α → 1, there is less and less averaging being performed, and Δip(t) →
Δip(t-1), and at α = 1 essentially Δip(t) = Δip(t-1). By recognising (4.49), the EWMA
equation resembles a discrete first-order low-pass filter.
y (t | ) = T (t ) (4.50)
where φ(t) is the regression vector whose elements are known and σ is the
unknown parameter vector whose elements are unknown and can be estimated by
means of LSE. Suppose y(t) is the measured data, the task is to estimate the unknown
parameters such that, y(t)-φT(t)σ is minimized. In this algorithm, the non-recursive
LSE is transformed to a recursive LSE based on the technique proposed by [165].
Referring to Figure,4.1, the input signal to LSE parameter estimation block is the
difference fault signal, ∆i, and suppose the output of the model is ŷ(t), and
considering sinusoidal nature of the power system signal, then the model output can
be described by a sum of N sinusoids as defined by (4.51)[166].
where αs,n = Ansin(δn) and αc,n = Ancos(δn) with An and δn respectively being the
amplitude and the phase angle of the nω0 frequency component. ∆T is the sampling interval
as previously defined. From (4.51), the regression vector can be defined as having elements
given by (4.52).
100 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
where φ(t) is the regression vector, and ωo is the fundamental frequency. The
elements of the regression vector are known since ω0 and ∆T are known. The
unknown parameters vector whose elements are the coefficients of the sine and
cosine terms of (4.51), to be estimated by (4.53)
= ( T ) T i
−1
(4.54)
where (φTφ)-1φT is the pseudoinverse of the regression vector. Let σN-1 denote
the non-recursive estimate of the unknown parameter, then (4.54) can be written as
(4.55). The unknown parameter can now be estimated recursively using (4.56) [165]
N −1 = ( T ) T i
−1
(4.55)
where (∆i(n) - φ(n)σN-1) is the error between the model output and the
measured data at the sample point n∆T. G(n) is an element of the time varying gain
vector G at the same sample point. The elements of G are recursively updated
according to (4.57) [165]
Using the estimated parameters, the instantaneous fault signal at any sample
point, n∆T can be calculated using (4.58); and the instantaneous magnitude and the
phase angle is determined using (4.59) and (4.60) respectively.
f (n ) = s (n ) s (n ) + c (n ) c (n ) (4.58)
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 101
s
(n ) = tan −1 (4.60)
c
I inc =
(I Fault − I Prefault )
(4.61)
I Prefault
102 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
overcurrent threshold is set about 2 to 3 times the prefault load current, the OC
threshold value based on current magnitude is IThreshold = 2IPrefault, such that;
I Fault
M Threshold = (4.62)
2 I Prefault
where MThreshold is the ratio of the fault current IFault to IThreshold and is taken as OC
threshold parameter in the proposed scheme. Using (4.62), MThreshold can have a
minimum and maximum of 0.5 and ≥ 1.0 per unit and correspond to MThresholdmin and
MThresholdmax respectively. The increase in current magnitude, Iinc is compared with
MThreshold to declare a fault. The following scaling factor is applied to Iinc for
comparison, MThresholdmin(1+Iinc) such that, Iincmin = 0.5(1+0) = 0.5 and Iincmax =
0.5(1+1) = 1.0. Considering noise and other factors including transient disturbances,
the threshold is set at 15% above prefault value such that MThresholdmin = 0.575 and
MThresholdmax = 1.15 per unit. It must be noted that the 15% above minimum and
maximum values of MThreshold for noise consideration is not based on any criteria,
however has been selected based on the value used by Gautam and Brahma [78] .
This value is subject to review for different fault conditions, network configuration,
level of DER penetration, noise level, etc. This, for initial fault detection, if 0.575 <
MThreshold < 1.15, then DD flag is activated, followed by L2 flag to initiate HIF
feature extraction and classification process discussed in Chapter 5. For OC
parameter setting, the maximum value is considered, thus a fault current increase of
greater than or equal to 1.15 pu must be detected to initiate the OC protection.
Supposing a fault occurs such that, Iinc > 1.15, then Timer and Counter
parameters are used to count the consecutive number of tall edge MFD spikes with
increased heights above the OC MThresholdmax of 1.15 pu within OC fault detection
delay (Tdmax) of 15.3125ms. This is the time delay in which a fault must be declared
and is considered as the waiting time. Moreover, within the fault detection time
delay, a minimum of two MFDTall spikes must be detected to declare OC fault. To set
the counter value, the following has been considered. Within Tdmax, there are five
edge spikes, and the edge spikes can come in the following combinations, two
MFDTall spikes with three MFDShort spikes or three MFDTall spikes with two MFDShort
spikes. Considering that fault detection is based on MFDTall spikes, a minimum of
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 103
two MFDTall spikes exist, hence, the counter value is set at two. It must be noted that,
if two tall edge spikes are detected before Tdmax lapses, the algorithm will not activate
L1 flag, both Counter and Timer conditions must be satisfied for the OC protection
procedure to proceed any further.
104 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
4.7.1 Simulation System Description
To test and verify the effectiveness of the proposed adaptive overcurrent
protection scheme for radial distribution network feeder, and test system was
modelled in MATLAB/Simulink. The layout of the testbed for simulations to analyse
the impact of PV penetration on radial distribution feeder overcurrent protection is
shown in Figure 4.7. The testbed was modelled using power system device objects
available in Simscape tool box in the Simulink library. The distribution feeder was
modelled to operate at a nominal system voltage of 25kV with system frequency of
50Hz. The synchronous generator source at the substation represents power supply
from the grid, and has total generating capacity of 5MW. The total distance of the
feeder is 30km, thus a short line model is used for each line segment between the PV
systems with positive- and zero- sequence impedances, z(1) = 0.1153 + j0.53 Ω/km
and z(0) = 0.413 + j1.043 Ω/km respectively. The system is comprised of four PV
systems connected along the feeder length. The distance to PV systems from the
feeder substation are; PV1 = 0.5 km, PV2 = 10 km, PV2 = 20 km and PV4 = 30 km.
Each tapped load is P = 1.1MW and Q = 0.15 Mvar. The complete distribution feeder
system with PV penetration incorporates measurement, interfacing as well as control
devices.
260V/25kV
Substation 50Hz
SS PCC
CT Δ/Y
FP 1 FP 2 FP 3 FP 4
VT
25kV, DOCAS Load Load Load Load Load
50Hz
V, I Measurements
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 105
fault current signal as input to the DOCAS algorithm. In the simulation, the CT is
model with turns ratio of 500:1 and operates in the linear of the saturation
characteristic curve. The efficiency of the CT has not been considered in the
modelling on the assumption that, losses are minimal, thus CT operates at near rated
efficiency.
Distribution transformers: Each PV system is connected to the main feeder
through a 600kVA, 50 Hz, 260V/25kV delta/wye distribution transformer with wye
grounded. The transformer generally provides a galvanic isolation between the PV
system and main feeder.
Photovoltaic Systems: Each PV system is modelled identically to generate
500kW at 260V DC. A PV system is comprised of five 100kW PV strings. Each 100
kW PV string has 64 parallel strings with 5 series connected PV modules. The PV
module used are SunPower SPR-315E-WHT-D found in the Simscape library in
MATLAB/Simulink [167] with following characteristics; maximum power =
315.072 W, open circuit voltage, Voc = 64.4 V, short circuit current, Isc = 6.14 A.
The voltage and current at maximum power point are Vmpp = 54.7 V and Impp = 5.76
A. The equation that models the current output of each PV module is give in (4.1)
[168], [169].
V + IRs V + IRs
I = N pp I PV − I o exp − 1 − (4.63)
Vt N ss Rp
where Nss and Npp are number of series and parallel connected cells in the PV
module. The parameters are defined as follows:
IPV is the current generated by the incident light
Io is the diode reverse saturation current
Vt is the thermal voltage of PV module, such that Vt = NsskT/q, where k
is the Boltzmann constant (1.3806503 x 10-23 J/K), q is the electron
charge (1.60217646 x 10-19 C), T is the temperature of the p-j
junction.
Rs, RP are the series resistance and parallel resistances respectively.
The current versus voltage (I-V) and power versus voltage (P-V) curves at
standard test condition (STC), 1000W/m2, 25ºC for each 100kW PV string is shown
106 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
in Figure.4.8. The system voltage (VMPP) and current (IMPP) at maximum power point
(MPP) respectively are 273.5 V DC and 368.6 A.
Figure 4.8: Characteristic curves for the PV strings at STC, (a) I-V and (b) P-V curves
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 107
To analyse the boost converter circuit in Figure 4.9, we consider the switching
intervals of Q1 and D1.
Interval 1: Q1 is in “on” state and D1 “off”, the sections of the circuit after and
before the diode become isolated. The inductor L is shorted to ground. Under this
condition we derive the following equations [171];
vL = Vin (4.64)
Vout
iC = − (4.65)
R
It is assumed that voltage and current ripples are very small, thus iL = I; then
Vout
iC = I − (4.68)
R
Using the switching intervals, and the voltage across the inductor, v L(t) and the
capacitor through current, iC(t) can be represented graphically as shown in Figures
4.10(a) and (b) respectively. The area under the curves in the plots of Figure 4.10(a)
and (b) are equated to zero in steady state as given by (4.69) and (4.70) respective.
Figure 4.10: Switching waveforms for the voltage and current in the DC-DC boost converter
108 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
0 = D(Vin ) + (1 − D)(Vin − Vout ) (4.69)
Vout V
0 = D( − ) + (1 − D)( I − out ) (4.70)
R R
The parameter D is the on/off duty cycle where D = Vout/Vin = Ton/Toff. The steady
state output voltage and the inductor current can be respectively solved from (4.69)
and (4.70) as defined by equations (4.71) and (4.72) respectively.
1
Vout = Vin (4.71)
1− D
1 Vout 1 Vin 1 Vin
I= = = (4.72)
1− D R (1 − D ) R 1 + D2 R
2
Equation (4.71) shows that, by increasing D in the range, 0 < D < 1, the output
voltage can be made larger than the input voltage. At D = 0, the output voltage is
equal to the input voltage. The output current is determined according to (4.71), and
by inspection, the output current will be lower than the input current for the same
value of D applied in both equations.
In PV systems, the voltage and current are unregulated due to the changing
solar irradiance and temperature. This can be visualised by varying the temperature
in the 100kW PV string in the test system. Plots of I-V and P-V curves at
temperatures of 25ºC, 30ºC and 40ºC for the 100kW strings are shown in Figure
4.11.
Figure 4.11: Characteristic curves for the PV strings at STC, (a) I-V and (b) P-V curves
with increased temperatures
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 109
Under test these conditions the PV strings operate at lower MPP where the
respective VMPP are; 273.5 V, 268.6 V and 259.5V. The decrease in power is related
to the temperature coefficient of material in the PV modules [172]. There is minimal
change in the current due to temperature change where the respective currents at the
new MPP due to changing temperature are; 368.6 A, 369.7 A, and 370.9 A. The
impact of current is normally significant with changing irradiance. Under different
ambient conditions (change in temperature and irradiance) which the PV systems are
subjected to, the MPP changes. It is expected that under these conditions, the DC bus
voltage must be kept constant. This can be achieved by appropriately scaling the
value of D in equation (4.8).
Maximum Power Point Tracking (MPPT): In PV systems, the MPPT
algorithm is used in conjunction with the DC-DC boost converter by appropriately
adjusting the value of D to maintain a constant DC voltage output at different MMP.
Several different MPPT algorithms have been proposed, including Perturb and
Observe, Incremental Conductance and Fuzzy logic methods [173], [174]. The
MPPT algorithm used in the simulation system is the Perturb and Observe (P&O)
algorithm. The choice of this algorithm was based on simplicity and speed of
execution based on only two inputs as compared to the other popular method, the
incremental conductance method which requires four input parameters including
incremental and instantaneous PV array conductance and VMPP and IMPP [175], [176].
Moreover, the P&O algorithm was selected on the assumption that the simulated
faults occur under steady state environmental conditions (slowly changing irradiance
and temperature). Furthermore, the algorithm is still widely used. The inherent lack
of speed in tracking the MPP and loss of power when subjected to rapidly changing
irradiance and temperature widely reported in literature regarding the P&O algorithm
can be reduced by increasing the sampling rate and the execution speed [177],[178].
The P&O algorithm works by periodical perturbation (increment or decrement)
of the array voltage and comparing the array output power with power output of the
previous perturbation cycle. If the power increases, the perturbation will continue in
the incremental direction in the next perturbation cycle, if otherwise, the perturbation
will reverse and go in the decremented direction. This process continues until the
new maximum power point is reached. To further illustrate the P&O algorithm,
supposing RMPP is the array characteristic impedance, such that RMPP = VMPP/IMPP,
and R is the array load impedance, for maximum power transfer to occur, R MPP = R.
110 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
RMPP depends on the operating environmental conditions, irradiance and temperature.
The P&O algorithm either increment or decrement the duty cycle D to control VMPP
to match RMPP to R.
DC-AC Inverter: The DC to AC converter used to interface the PV systems to
the radial feeder via the distribution transformer is the Voltage Source Converter
(VSC)[179] . The model used is the two-level average VSC model with six pulse
insulated gate bipolar transistor with inverse-parallel diodes. The VSC takes the 500
V DC from the DC-DC boost converter and converts it to 260 V AC, three-phase, 50
Hz output. A brief description of the VSC operation is provided herein. Figure 4.12
shows topology of the two-level VSC and the AC voltage output at one of the
phases. The two-level VSC model is an average model. As it is obvious from the
graph of output waveform, the positive and negative halves of the AC voltage output
are relative to the positive and negative half values of the DC input voltage taken at
the midpoint. There are three legs in the SVC topology with upper and lower IGBT
switches, S1 to S6. These legs, call phase-legs generate the phase AC voltages by
appropriate (complimentary) switching of the upper and lower IGBTs connected to
that leg. S1, S3 and S5 are the upper switches for phases A, B and C respectively
while S2, S4 and S6 are the lower switches for the respective phases. All the upper
switches are simultaneously in the “On” state for 180º (half fundamental period) with
phase shift of 60º between each other obtained through successive gating signal
between each IGBT. In the next 180º all the upper IGBT switches are in the “Off”
state and the lower IGBT switches are in the “On” state. Any simultaneous switch on
of the upper and lower IGBT switches in a phase-leg will result in short circuiting
that leg, thus this is avoided. The magnitude of the fundamental component of
generated AC output voltages at each phase is controlled by pulse width modulation
(PWM). Moreover, the PWM reduces harmonic content present in the voltage (and
current) signal.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 111
Figure 4.12:Circuit topology of two-level voltage source converter with a phase output
voltage waveform [179]
112 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Figure 4.13: Fault current signals with DC-offset, (a) The fault current signal and average
MMF output (b) Difference fault current signal.
Figure 4.14: Power spectral density plots for (a) fault current (b) average MMF output and
(c) difference fault current, ∆i.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 113
Figure 4.15: Difference fault current signal, ∆i and the EWMA filter output
Recursive Least Square Filter Output: The EWMA filtered version of the
difference fault current signal, ∆i provides input to the recursive least error square
filter for phasor estimation. The input and the output signals are presented in Figure
4.17 for comparison. It shows that the two signals are in phase and almost exact
replication of each other. Thus, it can be concluded that the magnitude estimation is
accurate.
Figure 4.17: RLSE filter signals (a)RLSE filter input signal, (b) RLSE filter output signal
114 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
lumped load at 5 MW. The substation source has a total capacity of 5 MW while the
total PV generation is 2 MW, thus decrementing the substation source contribution
allows for increased PV contribution to maintain the 5 MW supply. The PV systems
were maintained at a fixed irradiance of 300 W/m2 at 30°C for all simulations.
Prefault Condition: The substation source was decremented to allow
increasing PV penetration at 0, 28, 33, 44 and 50 percent. The MFD values for
prefault current magnitudes at these PV penetration levels for each phase of the
feeder are given in Table 4.1, and the results show that as the PV penetration level
increases, the MFD values decrease indicating decrease in the current magnitude at
the substation.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 115
Using the MFD data recorded in Tables 4.1 and 4.2, the fault current increase in each
phase can be calculated using (4.61) for different faults types at the specified location
for case. Cases of fault current increase in phase A under the simulated fault
conditions without PV penetration are given in Table 4.3. The other two phases can
be treated in the same manner. The results show that, even without any PV
penetration, the fault current magnitude decreases with increased fault distance.
Table 4.4: MFD values fault currents for faults at Fault location 2
MFD Output – Fault Current Magnitudes at
Fault
Phase PV Penetration Level (%)
Types
0 28 33 44 50
SLG A 0.0096 0.0086 0.0071 0.0060 0.0053
A 0.0088 0.0083 0.0077 0.0059 0.0047
LLG
B 0.0098 0.0082 0.0067 0.0055 0.0054
A 0.0082 0.0081 0.0065 0.0050 0.0042
3-Ph B 0.0079 0.0077 0.0062 0.0052 0.0040
C 0.0077 0.0071 0.0064 0.0040 0.0039
116 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Figure 4.18: RLSE filter magnitude response for SLG fault at fault location 2 at 0%, 28%
and 33% PV penetration
Figure 4.19: MFD output corresponding to fault current magnitude for SLG fault at fault
location 2 at 0%, 28% and 33% PV penetration.
Table 4.5: Fault current increase at various PV levels for faults at location 2
Inc. in Fault Mag. using (4.61)
Fault
at given PV Levels (%)
Types
0 28 33 44 50
SLG 2.0000 5.6154 5.4545 5.0000 4.8889
LLG 1.7500 5.3846 6.0000 4.9000 4.2222
3-Ph 1.5625 5.2308 4.9091 4.0000 3.6667
The maximum fault current under any fault condition would occur at fault
location nearest to the relay, and in this case nearest to the feeder substation, and at
zero PV penetration. The maximum fault current occurs at location 1 (FP1) at 0%
PV, thus referring back to the results in Table 4.3, the highest increase in fault
current magnitude is 2.1875 pu for SLG fault on phase A. All increase in fault
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 117
current magnitudes in phase A will be lower than this value, thus this is defined as
Iincmax. With increasing PV penetration, any fault along the feeder length will result in
lower current magnitude compared to Iincmax. To account for this, fault increases at
different PV penetration levels and fault locations are made relative to Iincmax. The
relative increase in fault current magnitude from without PV penetration to with PV
penetration can be computed according to (4.73);
I Fault ( PV )
I incRel = I inc max (4.73)
I Fault ( No _ PV )
where IincRel. is the increase in fault current magnitude relative to the maximum
fault current increase without PV penetration. IFault(PV) and IFault(No_PV) are the fault
current magnitudes with and without PV penetration respectively represented by
respective MFD values.
Using MFD data in Tables 4.1 and 4.4, the relative increase in fault current
magnitude at the given PV penetration levels for various fault conditions at fault
location 2, calculated using (4.73) are given in Table 4.6. The trend shown by the
results in Table 4.6 indicate that, the fault current magnitude decreases with
increased PV penetration level.
Table 4.6: Increase fault current magnitude at different PV level for fault a location 2
Relative pu Increase in Fault Current Magnitude
Fault
at given PV Levels using (4.73)
Types
0% 28% 33% 44% 50%
SLG 2.0000 1.7917 1.4792 1.2500 1.1042
LLG 1.7500 1.6506 1.5313 1.1733 0.9347
3-Ph 1.5625 1.5434 1.2386 0.9527 0.8003
Supposing the main feeder relay has total reach up to location 4, then a fixed
pu threshold (pickup) value using the prefault and fault MFD values for the same
fault at location 4 in Tables 4.1 and 4.2 respectively can be calculated using (4.62)
giving 1.4015625 pu. A trend can be observed by comparing the different pu fault
current increases in Table 4.6 with the fixed threshold value which suggests that the
relay will under reach for some fault conditions at 33% PV penetration and all faults
occurring at 44% and 50% PV penetration levels for faults at location 2. This trend is
graphically illustrated in Figure 4.20. In Figure 4.20(a), the MFD values for the fault
current magnitude decrease with increased PV penetration. Similarly, the pu fault
level increase using (4.73) in Figure 4.20(b) decreases with increased PV
118 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
penetration, and begins to under reach at 33% PV where the fault increase drop
below Mpickup which is a fixed threshold. In Figure 4.20(c), a different trend is
established using (4.61). The trend generally shows that the pu increase in fault
current magnitude decreases with increased PV penetration. However, the pu
increase in fault current magnitude is significantly higher than the actual increase in
fault magnitude as indicated by the trend in Figure 4.20(b). This trend helps define
the basis for the adaptive OC protection scheme developed in this thesis and
described in section 4.8.
Figure 4.20: MFD output corresponding to fault current magnitude for SLG fault at fault
location 2 at 0%, 28% and 33%, 44% and 50% PV penetration levels.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 119
Figure 4.21: RLSE filter magnitude response for SLG fault at fault locations 1, 2 and 3 at
28%
Figure 4.22: MFD output corresponding to fault current magnitude for SLG fault at fault
locations 1,2 and 3 at 28% PV penetration.
Table 4.7: Fault level increase at different fault location at 28% PV penetration
Increase in Fault Current Magnitude relative to maximum fault current
Fault
at 28% PV at given Locations using (4.73)
Types
1 2 3 4
SLG 1.9945 1.7917 1.6209 1.4520
LLG 1.8032 1.6506 1.4702 1.3529
3-Ph 1.6057 1.5434 1.4438 1.3638
The impact on the fault current can be further analysed by using (4.73) to
calculate the relative increase in fault current magnitude with increased fault distance
under PV penetration. The results from this calculation are tabulated in Table 4.9,
and graphically illustrated in Figure 4.23.
Figure 4.23(a) shows the decrease in fault current magnitude due to increased
fault distance without PV penetration using (4.73). It can be observed that the relay
120 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
can maintain its distance of coverage without under reaching. Figure 4.23(b) shows
similar trend as in Figure 4.23(a) the case of increased fault distance with 28% PV
penetration. However, in this case the reach of the relay is decreased and the relay
under reaches at location for the LLG and 3-Ph faults when a fixed pick up threshold
is applied. The reach of the relay is affected by increased PV penetration, and as the
distance to fault increases, the relay will not be able to maintain its coverage of the
feeder.
Consider the case of the SLG fault, at 0% PV penetration, a SLG fault at
location 1 registered an increase of 2.1875 pu (Iincmax), and at location 4, the same
type of fault registered an increase in fault current magnitude of 1.6875 as given in
Table 4.3. This is a decrease of 23%. Consider the same fault at 28% PV penetration.
At fault location 1, the increase is 1.9945 and at fault location 4 the increase is
1.4520 pu as shown in Table 4.9. This is a decrease of 25%. When considered in
terms of Iincmax, the SLG fault at location 4 with 28% PV penetration results in a
reduced fault magnitude of 38%. This trend shows that increased PV penetration will
have an impact on the existing feeder protection system device coordination based on
passive unidirectional current flow. Moreover, it also shows that the relay will under
reach much faster than it would have without PV penetration.
While it was shown that fault current magnitude decreases with increased PV
penetration, a different trend was observed when directly applying (4.61), where the
increase in fault current magnitude is quite significant even with PV penetration as
was shown in Tables 4.5 and 4.8, and graphically illustrated in Figure 4.23(c). The
increase in fault current magnitude determined using equation (4.61) is adaptive to
any change in current magnitude influenced by PV penetration, load change and
network topology change, etc. Moreover, the threshold parameter defined in equation
(4.62) is made adaptive due to the memory update every 9.6875ms defined in section
3.5.5. Thus, a scheme for adaptive inverse time overcurrent (ITOC) relaying based
on equations (4.61) and (4.62) and applying appropriate scaling to the fault current
increases determine using equation (4.73) given in Tables 4.6 and 4.9 can be
designed and appropriately matched to a standard ITOC relay curve.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 121
Table 4.8: Increase in fault current magnitude at different fault distance 28% PV
Inc. in Fault Mag. using (4.61)
Fault
for 28% PV at given Locations
Types
1 2 3 4
SLG 6.1538 5.6154 5.2538 4.6923
LLG 5.8446 5.3846 4.8462 4.4625
3-Ph 5.3846 5.2308 4.9231 4.6923
Figure 4.23: MFD output corresponding to fault current magnitude for SLG fault at fault
location 2 at (a) 0% and 28% PV penetration.
The inverse time overcurrent relaying (ITOC) parameters for adaptive feeder
overcurrent protection are determined as defined herein.
122 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Ipickup |. The values of M must fall within a range defined by Mmin < M < Mmax, where
these are defined as follows;
Mmax = IFault_No_PV(max)/2IPrefault_No_PV(max)
The minimum value of Mpickup is set at 15 percent above Mmin, such that, Mpickup =
1.15Mmin. The value for Mmax is determined by means of power flow algorithm and is
updated only when major network changes occur. Moreover, the Mmax is based on the
prefault and fault current values without the penetration of PV such that,
IFault_No_PV(max) is the maximum fault current without PV and IPrefault_No_PV(max) is the
maximum prefault current without PV.
Increased penetration of the PV reduces the fault current magnitude on the
feeder, thus to account for this, all values of M are calculated as a ratio of Mmax
defined by (4.74),
I Fault _ PV
M Rel = M max (4.74)
I Fault _ No _ PV (max)
where MRel. is the value of M relative to Mmax and IFault_PV is the fault current with PV
penetration. It must be noted that, for all fault conditions at any point on the feeder
under any PV penetration level, load condition or network topology, the value of
Mpickup has to be maintained. This is the minimum point on the Inverse Time
Overcurrent (ITOC) relay curve. Therefore, all other M values must be shifted to the
right on the relay curve (above Mpickup) to compensate for the reduction in fault
current magnitude with increasing PV level. The compensated value of M is then
given by (4.75);
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 123
values for prefault current recorded in Table 4.1 and MFD value for the maximum
fault current without PV penetration for the SLG fault at location 1 given in Table
4.2 is, (0.0102 / (2 x 0.0032)) = 1.59375. Following the procedure defined in (4.74)
and (4,75), the fully compensated MComp values for the SLG faults, taking into
consideration the effect of increased PV penetration and fault distance are calculated
and tabulated in Table 4.10 and Table 4.11 respectively. The MComp values are
adaptive, however they must operate on a certain relay operating curve to determine
the relay operating time.
Figure 4.24: Standard moderately inverse ITOC relay curve with M values in Table 4.10
124 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Table 4.9: Per unit increases (M) in fault current magnitude at different PV levels at FP2
MFD Trip Fault at different PV Levels at Location 2 (FP2)
Parameters 0% 28% 33% 44% 50%
Ipickup 0.0064 0.0026 0.0022 0.0022 0.0018
M 1.5000 3.3077 3.2273 3.0000 2.9444
MRel. 1.5000 1.3438 1.1094 0.9375 0.8281
MComp 2.6500 2.4938 2.2594 2.0875 1.9781
Table 4.10: Per unit increases (M) in fault current magnitude at different fault locations
MFD Trip Fault at different Locations for 28% PV
Parameters 1 2 3 4
Ipickup 0.0026 0.0026 0.0026 0.0026
M 3.5769 3.3077 3.0769 2.8462
MRel. 1.4531 1.3438 1.2500 1.1563
MComp 2.6031 2.4938 2.4000 2.3063
The corresponding operating times for the SLG fault under consideration are
shown in Tables 4.12 and 4.13 for PV penetration level and fault distance. Note that
the maximum trip time (Ttripmax) is 18.512s which corresponds to Mpickup = 1.15Mmin
irrespective of the PV penetration level and fault distance. Moreover, the minimum
trip time (Ttripmin) is also maintained as it corresponds to Mmax and is associated with
fault condition without PV penetration. In Table 4.12, it can be noted that the time to
trip (Ttrip) increases with increasing PV penetration for fault at the same location. To
obtain same trip time, the CTS trip parameters A, B and ρ have to be redefined or
redesigned as part of the future work. The trip times (also relay pick up times) in
Table 4.13 are consistent, where time to trip increases with fault distance.
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 125
By selecting appropriate TDS values faster or slower operating time can be
selected as illustrated in Figure 4.25.
4.9 CONCLUSION
126 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
fault current contribution from the PV systems, and consequently reducing the fault
current contribution from the substation source. Moreover, it was shown that with
increased PV penetration, and fault occurring at distance along the feeder length will
result in reduced Thevenin impedance seen at the point of fault allowing high fault
current contribution by the PV systems. The implication in terms of feeder OC
protection is that increased PV penetration will affect the fixed OC protection pickup
setting of the relays and ACRs thus affecting the protection coordination. The
theoretical analysis was verified through simulations. The proposed feeder OC
protection strategy using the DOCAS algorithm includes other features including
DC-offset suppression and fault current magnitude estimation. A fast DC-offset
suppression technique requiring no parameter estimation was achieved through
subtraction of the fault current and the average MMF output. This method incurred
no additional computational overhead as the process is an integrated part of the
algorithm in fault detection and imposes no additional delay in fault detection. A
feature for fault current magnitude estimation through a RLSE filter is incorporated
in the fault detection and diagnostic tool. The performance and functional attributes
of these features were verified through simulations and were seen to have performed
as expected without compromising the overall feeder OC protection using the
DOCAS algorithm.
The reduction in the fault current magnitude contribution by the feeder
substation source was analysed through simulation of fault conditions including
SLG, DLG, and 3-Ph faults under different PV penetration levels at fixed fault
location (FP2). It was observed that, with increased PV penetration, the fault current
contribution from the substation source reduced accordingly. Further simulations
involving the same fault conditions, at a fixed PV level (28%) and different fault
distance along the feeder length were simulated. It was observed that as fault
distance increased, the fault current magnitude seen at the feeder substation
decreased. Moreover, with PV penetration at 28%, it was observed that the fault
current contribution by the substation source was lower for the same fault conditions
at the same fault location for the case without PV penetration. This shows that the
relay will under reach faster with PV penetration and fault occurring at distance
further from the substation.
The maximum fault current would occur at a point nearest to the substation
without any PV penetration as demonstrated by the simulation results. All fault
Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs 127
currents; for faults along the feeder length either with or without PV will be less than
this current. The trends in the reduction of the fault current magnitude were verified
through simulation of faults at locations along the feeder and measuring the fault
current magnitudes at feeder substation. Moreover, the trends in the reduced fault
current magnitude were established through developing relationship as a ratio of the
maximum fault current measured without PV penetration.
While it was observed that fault current contribution from the feeder substation
source decreased with increased PV penetration and fault distance, a different trend
was observed when computing the actual percentage increase in percentage rather
than the actual magnitude using the concept of dynamic memory update. The MFD
value of the tall edge spike in the memory is used to calculate the percentage/per unit
increase. It must be reiterated that the MFD value in the memory before the fault is
used to determine the increase, and only updated when the fault is cleared. It was
observed through this analysis that the percentage increase in fault current magnitude
was significantly higher as compared to increase in the actual fault current magnitude
for each type of fault. Based on this observation, a strategy for declaring an OC fault
using the percentage increase in fault current magnitude in conjunction with the
adaptive feeder OC protection threshold (pickup) parameter was proposed. The
adaptiveness and scalability of the strategy was demonstrated by applying inverse-
time overcurrent (ITOC) protection strategy. A minimum time delay of ¾ cycles + 1
sample which is defined as Tdmax = 15.3125 ms is required to declare a OC fault
before the ITOC relaying is initiated to determine the trip time. The proposed method
showed that irrespective of the fault location under any PV penetration level, the
relay threshold (pickup) value was always maintained at the same pickup value.
Furthermore, it was observed that protection coordination can be achieved at
different PV penetration level by appropriately selecting relay curves at different
TDS values.
128 Chapter 4: Adaptive Overcurrent Protection in Active Radial Distribution Feeders with RE Based DERs
Chapter 5: HIF Detection and Classification in
Distribution Feeders
5.1 INTRODUCTION
The structure of the proposed algorithm for HIF detection and classification
developed from this research is shown in Figure 5.1. Notably, the block diagram
structure of Figure 5.1 is the DOCAS MFD algorithm.
The mathematical derivations and the operation attributes of this algorithm
were discussed in detail in Chapter 3, therefore no mathematical derivations are
provided in this chapter. This chapter discuss the application of the DOCAS MFD
algorithm for HIF detection and classification. The DOCAS MFD algorithm is a
multistage Morphological filter, constructed from two nonlinear MM filters called
the morphological median filter (MMF) and the Alternating Sequential Filters (ASF).
As was previously discussed, MM is a nonlinear image/signal processing technique
that analyses the topography of the input signals waveforms by means of a probing
signal called the structuring element (SE) in complete time domain. The SE is the
filtering signal that provides general functional attributes of the MM signal
processing technique. MM can detect seemingly insignificant changes in the
topography of the signal waveform being investigated, thus making it convenient for
the detection and classification of HIFs.
The DOCAS MFD algorithm extract features from the HIF current to detect
and classify HIFs based on two HIF identifying characteristics, 1) randomness, due
to randomly changing the effective fault resistance, Rf resulting in erratic fault
current waveform and, 2) arc extinction and re-ignition of the AC HIF arc around the
fundamental period resulting in a shoulder shaped unsymmetrical fault current
waveform.
B1 B2
Initial Fault Detection
Open-close ● Close
∑●
Sampled Fault
Signal Input
f ASF ○ +
Δf
MMF -+ -
f (i,v) Av Close-open ● MFD
MMF ASF ∑○ Open Output
A1 A2 ○
B1 B2
5.3.3 HIF Classification Using Feature Extraction Using MFD Output Signal
The L2 flag is activated to initiate the HIF feature extraction using the MFD
output from the voltage signal input denoted MFDv after Tdmax timer times out. The
targeted HIF feature characteristics are the randomness and HIF AC arc extinction
and re-ignition resulting in the shoulder-shaped unsymmetrical fault current
waveform. These HIF characteristics cause variation MFD value of the edges as well
as cause spikes to grow in specific regions of the MFD output signal.
Randomness: The randomness feature is extracted from the edge spikes;
MFDTall and MFDShort spikes as well as the random spikes in the fault windows
kN T
MFDvShort = k = 1,.., (5.2)
k ( N + 1) T
The appearance of the elongated edges spikes is non-uniform and random, and
can occur at any point from k = 0 to ∞ and k = 1,..∞, while HIF persists.
The third class of MFD spikes considered in HIF randomness characteristic are
the MFDvWindow spikes. The MFDvWindow spikes are generated by those transients
occurring away from the zero-crossing and the positive and negative peak. The
MFDvWindow spikes appear randomly inside the fault windows, and unlike MFD spikes
due to transients from device switching and/or short-circuit faults, the MFDvWindow
spikes from HIF will prolong and span the entire length (in time) of the MFD output
so long as the HIF persists. The HIF feature extraction/classifier for the HIF
randomness feature shown in the flowchart targets the regions defined in Figure 5.2
to extract the randomness feature by observing the MFDvTall, MFDvShort and
MFDvWindow spikes.
HIF Arc Extinction and Re-ignition: The arc-extinction and re-ignition
feature of the AC arc due to HIF is extracted by observing the MFD output around
the region defined by Arc Spikes in Figure 5.2. The HIF feature extraction/classifier
for the HIF arc extinction and re-ignition feature shown in the flowchart targets this
region to extract the arc extinction and re-ignition feature by observing the MFDvArc
spikes. The MFDvArc spikes are naturally generated by the intersection of the current
and voltage waveform once every half cycle. If the voltage and current signal
waveforms are in phase, then this happens at zero magnitude of each signal. If they
are not in phase, then at the point of intersection, the instantaneous voltage and
current values are equal. In a 50 Hz system, this repeats every 0.01s. A spike is
generated in the MFDv output at the point of intersection that correlates to this
The AC arc extinguishes and re-ignites twice every have cycle. The arc
remains extinguished for a short period of time while waiting for the voltage to
regain until it reaches the restriking level to overcome the gap separation break down
voltage for arcing to happen again. This phenomenon impacts the fault voltage at
specific and fixed period defined in (5.3). The arc extinction and re-ignition thus
causes MFDvArc to grow in height when HIF occurs. The MFDvArc spikes can be
separated to observe for the arc extinction and re-ignition HIF characteristics.
Rule 4: elseif no HIF alarm, then reset, Tr, DD flag and L2 flag.
HIF is only declared if both features are present. The condition for HIF must be
detected and declared within the time limits, thus it takes 1.0153125s from the point
of fault inception to reach a decision if HIF exists or not
Contact Surfaces: In the HIF model used, the parameters Vp and Vn model the
contact surfaces. It was experimentally shown by Emanuel et al [56] that, during
HIF, the HIF current is asymmetrical with the positive half cycle having higher value
than the negative half cycle[181]. Hence, to model this phenomenon, Vn must be
greater than Vp (Vn > Vp), and Vn – Vp = ΔV, where ΔV is unsymmetrical voltage.
Moreover, it was shown that, less densely packed contact surface (soil) yield higher
arc voltage than contact surface with high density. Using this as the guide, contact
surfaces in[182] were modelled to obtain the specified current magnitudes.
Furthermore, the values for Rp and Rn are specified in Table 5.1. These parameters
randomly vary between +10% of the specified steady state values and represent the
effective fault resistance for positive and negative half cycles respectively.
Moreover, the rate of variation of the resistance used for the respective surfaces are
40ΔT, 35ΔT, 30ΔT, 20ΔT, 15ΔT and 10ΔT going from contact surface 1 to 6. The
unequal positive and negative arc voltages and the randomly changing arc resistances
generate nonlinear V-I characteristics in each of the contact surfaces. This
phenomenon is graphically illustrated in Figure 5.6 where the V-I characteristic
curves for each contact surface are shown.
Table 5.2: Time duration for existence of the MFD spikes due to non HIF transients
MFDiTall Edge MFDvWindow
Event Node Phase Ts (s)
Value IInc (pu) Tdw (ms)
0.4066 1.8299
SLG 632 A 0.030 -
0.4279 1.9257
0.1364 0.6058
611 A 0.030 4.6875
0.1111 0.5000
0.0821 0.5782
100kVar 646 B 0.035 9.6875
0.0721 0.5077
cap
0.1110 0.6160
switching 652 C 0.037 12.031
0.0901 0.5000
0.1212 0.5455
675 ABC 0.033 7.0312
0.1111 0.5000
0.1936 0.8713
100 kW 632 ABC 0.035 9.6875
0.1942 0.8740
IM
0.1942 0.8740
switching 675 ABC 0.037 7.0313
0.1950 0.8776
0.1184 0.5329
645 AB 0.037 2.4375
50 kW step 0.1179 0.5306
Load Inc. 0.1328 0.5977
675 ABC 0.030 2.4375
0.1313 0.5910
It must be noted that, Table 5.3 only contains the time of appearance of the first
seven MFDv spikes that went above the lower threshold value within the upper time
constrain (Tr). It must be stated that, there is no Counter Classification to determine
the minimum number of MFDv spikes to determine if HIF exists. The first seven
MFDv spikes are recorded to show the random time sequence of appearance of the
spikes. Moreover, any assumptions of HIF detection for each test scenario are based
on the recorded data in Table 5.3.
It was noted during simulation that in all cases, the random MFDvWindow spikes
appeared and sustained for the duration of the HIF. As defined, the random
appearance of the MFDvTall, MFDvShort and sustained appearance of MFDvWindow
spikes indicate the HIF randomness characteristic. The randomness is signified by
the non-uniform time of appearance of the MFDvTall and MFDvShort spikes. The non-
uniform time sequence can be easily observed in the time data recorded in Table 5.3
for all cases. The HIF arc extinction and re-ignition characteristic is signified by the
MFDvArc, and it can be observed that in all cases, the arc extinction and re-ignition
characteristic is consistently present, and exists when arcing occurs.
The results in Table 5.3 must be tested case by case against the Decision Logic
to determine if the HIF condition has been detected. Referring to the Decision Logic
in Subsection 5.3.4, all cases except HIF phase C of line 605 contacting with dry
grass surface has not been detected based on the first seven random MFDv spikes.
While dispersed sustained MFDvWindow spikes were observed in this case, however,
145
Figure 5.11: Signals for HIF at 602 on conc. surface (a) fault voltage and HIF current, and
(b)MFDv output
Figure 5.12: MFDvTall edge spikes for HIF at 602 on conc. surface (a) MFDvTall (1), (b)
MFDvTall(1) pu increase,(c)MFDvTall(2) and (d) MFDvTall(2) pu increase.
Figure 5.14: Random MFDvWindow spikes for HIF at 602 on conc. surface
Figure 5.15: MFDvArc spikes for HIF at 602 on conc. surface, (a) MFDvArc spikes and
(b) MFDvArc spikes with pu increase.
where MDvInc is the measure by how much the MFDvFault deviates from the
MFDvPrefault, where both respectively are the prefault and fault MFDv values of each
target MFDv spikes used in HIF feature extraction except the random windows
spikes. The signal for HIF extraction is then, MFDvInc = 0.5(1+MFDvInc) – threshold
or if MFDvInc value ≤ Threshold, then MFDvInc value = 0.
The prefault values for the respective MFDv spikes for each phase of the
simulation test system are shown in Table 5.4.
The revised method is applied to the undetected HIF case on phase C of line
605, which is the line segment 684-611 and given as HIF case 3 in Table 5.5. The
revised threshold limits are given in the same table, and the result for this case is
graphically illustrated in Figures 5.16 – 5.19. It must be noted that the threshold
limits are lower than the level previously determined without considering the
distance to fault. While no empirical relationship is presented, the levels have been
Figure 5.17: MFDvShort edge spikes for HIF at 605 on dry grass (a) MFDvShort (1), (b)
MFDvShort (1) pu increase, (c)MFDvShort(2) and (d) MFDvShort(2) pu increase..
Figure 5.18: Random MFDvWindow spikes for HIF at 605 on dry grass
The DOCAS algorithm HIF feature extraction capability under noise condition
has been tested by corrupting the HIF fault current and faulted phase voltage signals
with different signal to noise ratio (SNR). The success of the technique as previously
observed relies on selecting threshold limits. The threshold limits for the different
MFDv spikes signals at the respective SNR are tabulated in Table 5.6. The values
given in the table were determined through simulations and appropriate adjustment
of the threshold to eliminate all prefault tall spikes associated with the noise.
The DOCAS algorithm MFDv output under noise condition can be graphically
illustrated by considering the case of HIF on phase A of line 602 contacting with
concrete represented by Figures 5.20 – 5.24. The HIF current and phase voltage
signals have been corrupted with white Gaussian noise with SNR of -10dB.
Figure 5.21(a) shows that none of the MFDvTall(1) edge spikes had increase
above the threshold limit. However, Figure 5.21(b) shows the appearance of MFDvTall
edge spikes above the threshold limit.
The presence of the randomness feature can be further observed in the
MFDvShort edge and the MFDvWindow spikes respectively shown in Figures 5.22 and
5.23 with the MFDvWindow spikes having sustained chaotic appearance with increased
Figure 5.20: Signals for HIF at 602 on conc. surface (a) fault voltage and HIF current, and
(b) MFDv output with noise.
Figure 5.22: MFDvShort edge spikes for HIF at 602 on conc. surface (a) MFDvShort (1), (b)
MFDvShort(1) pu increase,(c)MFDvShort(2) and (d) MFDvShort(2) pu increase with noise.
Figure 5.23: Random MFDvWindow spikes for HIF at 602 on conc. surface with noise.
The prefault MFDiTall edge values representing the current magnitude without
and with PV penetration are shown in Table 5.7.
The prefault MFDiTall edge values measured at the SS indicates a current
reduction of 18.02% because of PV penetration. The considered case of HIF is
initiated at 0.06s. Considering the case of HIF occurring without PV penetration, the
fault current increased by 10.66%. This is much lower than the normal OC threshold
limit of more than 2 to 3 times the prefault current to detect any fault. With PV
penetration the HIF results in a fault current magnitude increase of 7.43%. The net
effect of PV penetration on fault current magnitude measured at the SS is considered
relative to the prefault current without PV penetration resulting in the decrease in
current magnitude calculated as (1 – 0.0347/0.0394) x 100 = 11.3%. This shows that
increased penetration of PV can aggravate the difficulty in HIF detection using the
OC protection scheme. Moreover, the reduced fault current magnitude reaffirms the
dynamic behavior of the feeder due to penetration of RE based DERs such as PVs.
The plots in Figures 5.26 to 5.27 show the effect of PV penetration on the current
magnitude. Figure 5.26 (a) shows the fault current waveforms without and with PV
Figure 5.26: DOCAS MFDi outputs for HIF at XF1 closer to the feeder substation (a) MFDi
without PV, (b) MFDi with PV
Figure 5.27: DOCAS MFDi outputs for HIF at XF2 further from the feeder substation (a)
MFDi without PV, (b) MFDi with PV
The prefault MFDiTall edge values representing the current magnitude without
and with PV penetration are shown in Table 5.8. The effect on the substation current
due to HIF occurring at a distance further away from the substation in a radial
distribution network with PV integration is demonstrated. The decrease in the
prefault MFDiTall for without and with PV in Table 5.8 compared to Table 5.7 is
associated with the increase in the Thevenin equivalent impedance with increased
fault distance.
Figures 5.28 to 5.35 demonstrate the impact on fault current magnitude in
radial feeders with increased RE based DERs for HIF occurring at a distance further
away from the point of signal measurements. Figures 5.28(a) shows the fault current
signal waveforms with and without PV penetration and respective MFDi outputs
shown in Figures 5.28 (b) and (c). The net effect of PV penetration on fault current
magnitude measured at the SS for the fault at XF2 from prefault value without PV to fault
with PV results in a decrease in current magnitude calculated as (1 – 0.0323/0.0370) x 100 =
12.7%.
The fault location XF2 is towards the end of the feeder and between PV3 and
PV4. The results in Table 5.8 show that the MFDiTall values for the two PV systems
increased by 10% in response to the HIF. The increased current injection by PV3 and
PV4 caused feeder current seen at the SS to drop from its prefault value of 0.0353 to
0.0323 resulting in an 8.5% drop fault in current magnitude. This phenomenon can
cause the feeder relay to see this as a load reduction. Figures 5.29 (a) and (b)
graphically illustrate this condition where the MFDiTall (1) and MFDiTall (2) edge
spikes show the decrease in fault current magnitude seen at the feeder substation.
Such a condition will aggravate the deficiency in the feeder OC protection
system rendering it totally blind to the existence of this fault condition.
The results from these simulations reaffirm the deficiency of the OC protection
system in the detection and identification of HIFs thus requiring other methods such
as the DOCAS algorithm.
The condition defined in subsection 5.5.3, and for which the MFDiTall output for the
simulated case depicted in Figure 5.27(b) poses a serious challenge in HIF detection
using any conventional OC protection technique. HIF detection under such condition
will have to rely on feature extraction. The effectiveness of the DOCAS algorithm is
tested under this condition. For this purpose, additional HIF simulation has
performed at location XF2 under same load condition with fallen conductor on dry
grass contact surface. The fault current signals for the case of with and without PV
penetration are shown in Figure 5.28(a). The MFDi outputs of the DOCAS algorithm
for these scenarios are respectively shown in Figure 5.28(b). The trend in the
behaviour of the fault current magnitude under these scenarios is graphically
illustrated by extracting the respective MFDiTall (1) and MFDiTall (2) outputs (spikes)
as shown in Figure 5.29 (a) and (b) respectively. The HIF is initiated at 0.4s, and the
trend shows that HIF occurring without PV results in an increase in fault current
magnitude while in the case of HIF with PV penetration, the fault current magnitude
decreases. The effectiveness of the DOCAS algorithm in detecting HIF is tested
against HIF occurring at XF2 with PV penetration.
The procedure in HIF detection using the DOCAS algorithm is implemented
according to the flowchart in Figure 5.3. The first step is detection of disturbance
based on threshold classification using the MFDiTall edge spikes. The increase in
MFDiTall edge spikes is determined according to (5.5). The MFDiTall (1) and MFDiTall
(2) edge spikes with their respective pu increases above the cut-off threshold 0.575
pu are shown in Figure 5.30. The disturbance is dictated by MFDiTall (1) edge spikes
as shown in Figure 5.30(b). The MFDiTall (2) edge spikes show no increase above the
threshold as show by Figure 5.30(d). The increase in the MFDiTall (1) edge spikes fall
within 0.575 < IThreshold <1.0 which satisfies the requirement for further HIF
extraction based on DOCAS MFDv output based on the fault voltage signal input.
Figure 5.28: HIF signals for HIF at XF2 on dry grass (a) fault currents for with and without
PV, (b) MFDi outputs for with and without PV
Figure 5.30:MFDiTall spikes and MFDiTall pu increases for HIF on dry grass at XF2 with PV
(a) MFDiTall (1) spikes, (b) MFDiTall (1) pu increase, (c) MFDiTall (2) spikes, (c) MFDiTall (2)
pu increase
Figure 5.31: HIF signals for HIF at XF2 on dry grass (a) fault voltage with PV, (b) MFDv
output
Figure 5.33: MFDvShort spikes and MFDvShort pu increases for HIF on dry grass at XF2 with
PV (a) MFDvShort (1) spikes, (b) MFDvShort (1) pu increase, (c) MFDvShort (2) spikes, (c)
MFDvShort (2) pu increase
5.6 CONCLUSION
The strategy for HIF detection and classification require feature extraction
from the DOCAS MFD output. The MFD output from current signal input,
designated MFDi and voltage signal input designated MFDv are used in the
classification process. The segregation of the non HIF related features from the HIF
features are achieved through conceptualized classifiers including Threshold
Classifier, Timer Classifier, and MFD output Feature Extraction/Classifiers including
Randomness and HIF Arc Extinction and Re-ignition Classifiers. In threshold
classification, the Threshold Classifier uses the upper and lower level threshold
values to differentiate a possible HIF related disturbance from a non HIF related
disturbance. The Time Classifier is used to time out any non HIF related MFD spikes
generated by the transients from the disturbance. The Threshold Classifier uses the
MFDi tall edge spikes while the Timer Classifier observes the fault window of the
MFDv spikes existence of random window spikes.
Simulation case studies including, a case of SLG fault, 100 kVar capacitor
switching, 100 kW induction motor switching and 50 kW step load increase were
simulated using the IEEE 13 bus test system to validate the Threshold and Timer
Classifiers. It was observed that, through these classification processes, only those
disturbances meeting the characteristics of HIF are processed. The HIFs generally
result in lower fault current, thus Threshold parameter with lower and upper
threshold limits have been defined with noise tolerance of 15% above prefault MFDi
value for both upper and lower thresholds. It was further observed that, the window
spikes from MFDv output are short lived, and die out before the Timer Classifier
limit of 15.3125 ms. The Threshold and the Timer Classifiers operate simultaneously
once a disturbance is detected. If the condition of the Threshold and Timer
Classifiers are met, further HIF classification based on feature extraction is initiated.
6.1 INTRODUCTION
The DC arc fault detection technique presented herein proposes to use current
and voltage signal extracted by placement of appropriate sensors at the PV strings
and DC bus. The measured current and voltage signals are fed as input to the DC arc
fault detector which is centrally located within the vicinity of the inverter. Presence
of arc is detected by continuously monitoring and Analyzing variation in signals
extracted at the input of the inverter and on the PV array DC bus (before the DC-DC
converter). The structure of the DC arc-fault detector is shown in Figure 6.1.
The procedure for DC Arc-fault detection is the same as fault detection in the
AC power system. However, in the DC Arc-Fault detection, the concept of fault
detection window is not available because of the non-sinusoidal nature of the (DC)
signal. The onslaught of the arc is detected by appearance of spikes which increase in
height in response to the fast changing (rate of change of the random DC arc) noise
like DC arc characteristic that sustains when ignited. The spikes in the DOCAS MFD
output appear chaotically. The DOCAS MFD algorithm detects DC arc-fault by
means of detecting the chaotic behaviour of the random DC arc.
The fault voltage and current, Vdc and Idc are operated on at the MMF stage to
create the initial fault detection signals, ΔVdc and ΔIdc by subtracting the average
MMF output from the input signals to detect existence of DC arc-fault. The chaotic
behaviour of the DC arc is detected further at the ASF stages for the classification of
DC arc-fault.
Figure 6.2: Radial distribution feeder with PV penetration used in the simulation study
Figure 6.6: Layout of the PV modules in the PV strings for fault simulations.
The DC Arc model used in the simulation to generate the DC Arc-fault current
and fault voltage is the model developed by Stokes and Oppenlander [31] given in
(6.1);
Varc = (20 + 0.534 zg )I arc
1.12
(6.1)
Several cases of arc faults were simulated at points indicated on Figure 6.3.
The faults are generally defined as, parallel faults, for Fp1-Fp7 while Fs1-Fs4 are
series faults. The parallel faults can be reclassified as mismatch faults resulting in
voltage mismatch in the PV strings. Faults Fp1-Fp6 are mismatch faults while Fp7
will be considered as parallel fault.
Mismatch Arcing Faults: The mismatch faults Fp1-Fp3 and Fp6 bridge two
points at different voltages on the same PV string. Fp4 and Fp5 are mismatch faults
that bridge two points at different voltages on adjacent PV strings. The severity of
the mismatch faults can be evaluated base on their mismatch percentage. The
mismatch percentage is determined based on the number of PV modules in the string
that are affected by the fault. Consider a mismatch fault on the same string such as
Fp1. This fault has three PV modules affected, thus the mismatch percentage is (3/15)
x 100 = 20%. Now consider a mismatch fault between two adjacent strings such as
Fp4 affecting PV modules on string 2 (PVS2) and string 3 ((PVS2). There are 9
affected PV modules in PVS2 and 3 in PVS3, thus the mismatch percentage is ((6-
3)/15) x 100 = 40%. The mismatch percentages for all mismatch faults are
determined following this procedure. The fault currents obtained for each simulated
fault condition are tabulated in Table 6.2.
During fault, the unfaulted strings view the faulted string as a load, hence feed
current (negative current, Ireverse) to the faulted string, where this current is the
backfed current. The PV string load current under fault is;
where Ibackfed is the backfed current from the adjacent unfaulted PV strings. From
(6.4), the backfed current is;
The backfed current calculated according to (6.4) are tabulated in Table 6.3.
The general trend in the results in Table 6.3 show that the backfed currents increase
with increase in irradiance. However, the arcing faults simulated at the low
irradiance levels still resulted in insignificant backfed currents which indicate that,
for the simulated PV system, the backfed current at each faulting PV strings will not
reach 2.1 times load current of the strings at STC. Thus, the arcing fault will persist.
Observing the results in Tables 6.2 and 6.3, it can be noted that some mismatch
faults having equal mismatch percentage do not exhibit same fault characteristics.
Compare for instance Fp2 (or Fp3) with Fp4 having 40% mismatch. This alludes to
the fact that fault location and fault resistance also dictates the level of backfed
current. The faults Fp2 and Fp3 are in parallel to the faulted PV modules, but are in
series with unfaulted modules in the same string. The resistance of each PV module
in the string is given as, RPV = v/i, where v and i of the PV module at any irradiance
and is inversely proportional to the irradiance. The fault is in parallel to mRPV where
m is the number of faulted PV modules. In arcing fault, the arc channel is represented
by Rarc which is parallel to mRPV to give effective fault resistance as (mRPV x Rarc)/
(mRPV + Rarc). Thus, fault location influences the fault resistance, and consequently
the backfed current.
Table 6.2: Load current under simulated fault conditions at different irradiances
PV String Load Current (A) under Fault at Simulated Irradiances
Fault Mismatch
(W/m2)
Type (%)
50 75 100 150
Fp1 20 0.2691 0.4247 0.5804 0.8378
Fp2 40 0.2484 -0.1544 -0.1872 -0.3643
Fp3 40 0.2484 -0.1554 -0.1872 -0.3643
Fp4 40 0.2070 0.3089 0.4493 0.6192
Fp5 20 0.2691 0.4247 0.5804 0.8378
Fp6 60 0.1656 -0.5019 -0.7489 -0.4371
The proposed arcing fault detection method can be visualized through graphical
illustration of one of the simulated cases, where faults are initiated at 1.5s. Fault Fp1
was simulated on PV string 1 at all levels of irradiance; the graphical illustrations in
Figure 6.7 are related to irradiance of 50 W/m2. Figure 6.7 (a) shows the drop in PV
string 1 load current. Figure 6.7 (b) is the fault current transformed by the MMF
stage of the DOCAS algorithm. Figure 6.7 (c) is the difference fault current signal
for the DC arc-fault detection. In Figure 6.7 (d) is the MFD output showing the
onslaught of arcing fault at 1.5s. The algorithm detects the rate of change in the fault
current signal to generate spikes. The chaotic nature of the arc generates random
spikes that grow in height randomly to correlate to the occurrence and changes in the
arc fault current and voltage signals. Arcing is characterized by random and
persistent spikes where as in normal transients the spikes will be short lived.
Fault voltage signals in Figure 6.8 (b)-(d) illustrate the same phenomena based
on voltage signal input. Graphs in Figures 6.7 (a) 6.8 (a) are generated from signal at
the input of the inverter. It is not possible to locate the faulted string based only on
the output of the signal at the inverter input. Figure 6.9 shows the PV string fault
current signals. As it is obvious, fault, Fp1 occurs at string 1 as shown in Figure
6.9(a) where the load current drops as compared to the other unfaulted PV strings.
Figure 6.10 is the corresponding MFD out of the PV strings fault current signals.
While the MFD output of the faulted string, Figure 6.10(a) shows taller spikes
compared to the unfaulted strings, it alone cannot allude to the faulted string. Both
signals are needed to locate the faulted string.
Figure 6.8: Voltage signal measured at DC bus with associated signals for DC arc-fault FP1
Figure 6.9: MFD outputs for PV string currents for fault DC arc-fault FP1
where IPV is the sum of PV string currents and Ifault is the fault current through
the arc channel. The DC input current to the inverter drops as shown in Figure
6.11(a). Since it is an arcing fault, the presence of arc is detected by the appearance
of random and sustained spikes in the MFD output shown in Figure 6.11 (b). The
presence of the arc is also shown on the MFD output of the voltage signal in Figure
6.12 (b) as well. In Figure 6.12 (a) the initial decrease in voltage due to the fault is
shown. The MPPT quickly restore the voltage and puts the system to operate at a
lower maximum power (MPP).
The PV strings initially seem to increase their output as indicated in Figure
6.13 but due to the fast action of the MPPT, the PV strings maintain the mean load
current before the fault with sustained arc fault current generation. The
corresponding MFD outputs in Figure 6.14 also indicate the presence of the arc.
Figure 6.10:Current signal measured at DC bus and MFD out for DC arc-fault FP7
Figure 6.11: Voltage signal measured at DC bus and MFD out for DC arc-fault FP7
Figure 6.13: MFD outputs for PV string currents for DC arc-fault FP7
Series Arcing Faults: Four cases of series arc faults defined as, Fs1 to Fs4 were
simulated at the simulated low irradiance levels at the locations indicated on Figure
6.3. The general trend in the load current during all cases of series fault is graphically
illustrated by considering the case for series fault Fs2 at irradiance of 50 W/m2 shown
in Figure 6.15, and Figure 6.16 for load current and voltage respectively.
Series arcing fault introduces a restrictive arc resistance in series to the load
resulting in significant reduction in faulted PV string load current with almost no
backfed current from the unfaulted strings. This reduction is seen at the DC input of
the inverter. Figure 6.15(a) shows the drop in the load current at the inverter input.
The corresponding MFD output in Figure 6.15 (b) indicates the existence of
sustained DC arc-fault. The MPPT response by putting the PV system to operate at a
lower MPP as shown in Figure 6.16 (a) where the voltage is quickly restored to
500V. Even at the reduced MPP, the DC arc-fault is still sustained as indicated by the
respective current and voltage MFD outputs.
Figure 6.14: Current signal measured at DC bus and MFD out for DC arc-fault Fs2
Figure 6.15: Voltage signal measured at DC bus and MFD out for DC arc-fault Fs2
Figure 6.18: DC arc-fault current and MFD output high to low transition after fault
Figure 6.19: DC arc-fault voltage and MFD output high to low transition after fault
Figure 6.20: DC arc-fault current and MFD output high to low transition before fault
Figure 6.21: DC arc-fault voltage and MFD output high to low transition before fault
Transition from Low to High irradiance with existing Arc Fault: Arc fault
exists at irradiance of 50W/m2 then transitions into higher irradiance of 75W/m2.
This case demonstrates transition of the irradiance as the day progresses from
morning upwards. The arc fault occurs at 1.2s, and the arc persists as the irradiance
transitions from 50W/m2 to 75W/m2 at 1.3s as illustrated in Figure 6.23 (a) and
Figure 6.22: DC arc-fault current and MFD output high to low transition with fault
Figure 6.23: DC arc-fault voltage and MFD output high to low transition with fault
6.4 CONCLUSION
In this chapter, research outcomes from the thesis are summarised. The
significant research contributions are specified, and the benefits and importance of
the proposed methods are summarised. Finally, recommendations for future research
directions are suggested.
The main objective of this research work is to develop a fault detection and
diagnostic tool which can be applied in AC as well as in DC systems. MM which is a
nonlinear signal processing tool is used to develop the DOCAS algorithm. The MFD
output signal of the DOCAS algorithm is used to propose strategies for adaptive
radial distribution feeder OC protection, strategy for HIF detection and classification
based on feature extraction. The same algorithm is used for DC arc-fault detection in
PV systems.
The DOCAS MFD is a multistage filter algorithm that uses two nonlinear MM
called the MMF and two classes of ASFs filters called the open-close and close-open
ASFs. The performance of the algorithm was enhanced through the application of a
weighted eccentrically decreasing convex SE, g with five points and two unequal
slopes. The SE was decomposed into two smaller SEs, designated A1 and A2 each
with three points to preserve computational efficiency and applied at the decomposed
MMF filtering stages in cascade. Further two SEs designated B1 = A1 and B2 = g,
such that B2 > B1 where applied at the ASF stages.
The attributes of the DOCAS algorithm in AC system fault detection and
diagnosis has been evaluated through simulation and analysis of simple AC signal as
well as a SLG fault. The output of the algorithm provided characteristics and features
such as the fault detection windows and MFD edge spikes for application in adaptive
OC protection and HIF detection. Moreover, the attributes for HIF detection and
classification based on extraction of two HIF features; the randomness and arc
extinction and re-ignition were established. Furthermore, the MFD output provided
mechanism for the concept of memory update at a fixed time interval by the MFD
The fault detection and diagnostic tool also referred to as the decomposed
open-close alternating sequence (DOCAS) algorithm has been developed using the
mathematical morphology signal processing technique. This algorithm has been
extensively tested for performance verification and characterisation involving three
different fault scenarios, including OC fault detection in radial distribution network
with increased penetration of DERs, high impedance faults detection and DC arc-
faults detection in PV systems
• Computing the available maximum and minimum fault MVA values at every
point in the distribution network feeder where a protective device will be
located. This is equivalent to computing the maximum and minimum values
of the Thevenin equivalent impedance at those locations. Usually these data
are obtained via load flow and short circuit (fault) studies.
• The threshold limit and rating of the protective devices are determined based
on the fault MVA information obtained.
Table 7.1: Summary of OC Protection Schemes for Distribution Networks with Increased
DER Penetration
Method Main Feature Advantages Disadvantages
There are two types of protection schemes widely researched and proposed for
implementation in DC systems, and these are known as the non-unit protection and
unit protection [212]. The discussions contained herein highlight some of the
advantages and disadvantages of these schemes as an attempt at directing further
research in proposing and developing sound protection schemes for DC systems.
Non-unit Protection Scheme: The Non-unit protection scheme requires
placement of fault current interrupting devices in series with the element to be
protected or at the point of protection. This scheme operates on similar principle as
the OC protection in AC systems where current is used as the threshold metric, such
that, if the fault current exceeds the pickup setting of the protective devices which
violates their normal limit of operation defined by their thermal or magnetic
characteristics, the fault current is interrupted.
The Non-unit protection depends on accurately monitoring and/or measuring the
voltage and current signals as well as the time derivatives of voltage (dv/dt), current
(di/dt) and the impedance at the point of fault [198],[213]. There is merit in this
technique; however, there is no flexibility in the technique being adaptive
considering the fixed threshold setting and/or rating of the protective devices, and the
intermittent nature of RE based DERs as well as multisource fault current. The
threshold parameters in the protective devices in the scheme must be adaptive to
ensure secure and reliable protection. Moreover, while being sensitive, the Non-unit
protection scheme lacks the flexibility for selectivity. The protective devices will
instantaneously operate as long as their threshold limits are violated, which is an
advantage as backup protection. Furthermore, the instantaneous tripping
characteristic presents difficulties in enabling coordination with other protective
devices within the non-unit protection scheme [201],[211].
To ensure secure and robust protection of the DC system, improvement to the
non-unit protection scheme would require research into making this scheme adaptive
while enabling the flexibility for selectivity and coordination [212].
Unit Protection Scheme: The unit protection scheme relies on the use of
intelligent electronic devices (IEDs) such as digital relays as current sensors to
continuously monitor the status of the network as well as maintaining an active
communication network to transfer network status information (fault data) between
Bibliography 205
[18] G. Chen, D. Jiang, Z. Lu, and Z. Wu, "A new proposal for solid state fault
current limiter and its control strategies," in Power Engineering Society
General Meeting, 2004. IEEE, 2004, pp. 1468-1473: IEEE.
[19] Z. Lu, D. Jiang, and Z. Wu, "A new topology of fault-current limiter and its
parameters optimization," in Power Electronics Specialist Conference, 2003.
PESC'03. 2003 IEEE 34th Annual, 2003, vol. 1, pp. 462-465: IEEE.
[20] A. Girgis and S. Brahma, "Effect of distributed generation on protective
device coordination in distribution system," in Power Engineering, 2001.
LESCOPE'01. 2001 Large Engineering Systems Conference on, 2001, pp.
115-119: IEEE.
[21] K. Kaygusuz, "Energy for sustainable development: A case of developing
countries," Renewable and Sustainable Energy Reviews, vol. 16, no. 2, pp.
1116-1126, 2012.
[22] S. Xue, F. Gao, W. Sun, and B. Li, "Protection principle for a DC distribution
system with a resistive superconductive fault current limiter," Energies, vol.
8, no. 6, pp. 4839-4852, 2015.
[23] M. K. Alam, F. Khan, J. Johnson, and J. Flicker, "A comprehensive review of
catastrophic faults in PV arrays: types, detection, and mitigation techniques,"
IEEE Journal of Photovoltaics, vol. 5, no. 3, pp. 982-997, 2015.
[24] V. B. R. Commission, "Power line Bushfire Safety Taskforce, Final Report
Chapter 3.4. 1. Rapid Earth Fault Current Limiter, Melbourne," Victoria,
2011.
[25] S. M. Madani, "Analysis and design of power system protections using graph
theory," 1999.
[26] A. G. Phadke and J. S. Thorp, Computer relaying for power systems. John
Wiley & Sons, 2009.
[27] Y. G. Paithankar and S. Bhide, Fundamentals of power system protection.
PHI Learning Pvt. Ltd., 2011.
[28] M. Baran and I. El-Markabi, "Adaptive over current protection for
distribution feeders with distributed generators," in Power Systems
Conference and Exposition, 2004. IEEE PES, 2004, pp. 715-719: IEEE.
[29] S. H. Horowitz and A. G. Phadke, Power system relaying. John Wiley &
Sons, 2008.
[30] S. Zocholl, J. Akamine, A. Hughes, M. Sachdev, L. Scharf, and H. Smith,
"Computer representation of overcurrent relay characteristics: IEEE
committee report," IEEE transactions on power Delivery, vol. 4, no. 3, pp.
1659-1667, 1989.
[31] C. A. Kramer and W. A. Elmore, "Flexible inverse overcurrent relaying using
a microprocessor," IEEE transactions on power delivery, vol. 5, no. 2, pp.
915-923, 1990.
[32] M. V. Yalla, "A digital multifunction protective relay," IEEE Transactions on
Power Delivery, vol. 7, no. 1, pp. 193-201, 1992.
[33] IEEE STD C37.112-1996: Inverse-Time Characteristics Equations for
Overcurrent Relays, 1996.
[34] G. Benmouyal et al., "IEEE standard inverse-time characteristic equations for
overcurrent relays," IEEE Transactions on Power Delivery, vol. 14, no. 3, pp.
868-872, 1999.
[35] Z. N. Stojanovic and M. B. Djuric, "Table Based Algorithm for Inverse-Time
Overcurrent Relay," Journal of Electrical Engineering, vol. 65, no. 4, pp.
213-220, 2014.
206 Bibliography
[36] K. A. Wheeler, S. O. Faried, and M. Elsamahy, "Assessment of distributed
generation influences on fuse-recloser protection systems in radial
distribution networks," in Transmission and Distribution Conference and
Exposition (T&D), 2016 IEEE/PES, 2016, pp. 1-5: IEEE.
[37] P. Mohammadi and S. Mehraeen, "Challenges of PV integration in low-
voltage secondary networks," IEEE Transactions on Power Delivery, vol. 32,
no. 1, pp. 525-535, 2017.
[38] S. M. Brahma and A. A. Girgis, "Development of adaptive protection scheme
for distribution systems with high penetration of distributed generation,"
IEEE Transactions on power delivery, vol. 19, no. 1, pp. 56-63, 2004.
[39] S. Javadian and M.-R. Haghifam, "Implementation of a new protection
scheme on a real distribution system in presence of DG," in Power System
Technology and IEEE Power India Conference, 2008. POWERCON 2008.
Joint International Conference on, 2008, pp. 1-7: IEEE.
[40] P. Mahat, Z. Chen, B. Bak-Jensen, and C. L. Bak, "A simple adaptive
overcurrent protection of distribution systems with distributed generation,"
IEEE Transactions on Smart Grid, vol. 2, no. 3, pp. 428-437, 2011.
[41] C.-S. Chen, C.-T. Tsai, S.-C. Hsieh, C.-T. Hsu, and C.-H. Lin, "Adaptive
relay setting for distribution systems considering operation scenarios of wind
generators," in Industrial & Commercial Power Systems Technical Conf
(I&CPS), 2013 IEEE/IAS 49th, 2013, pp. 1-8: IEEE.
[42] H. Yazdanpanahi, Y. W. Li, and W. Xu, "A new control strategy to mitigate
the impact of inverter-based DGs on protection system," IEEE Transactions
on Smart grid, vol. 3, no. 3, pp. 1427-1436, 2012.
[43] C. Su, Z. Liu, Z. Chen, and Y. Hu, "An adaptive control strategy of converter
based DG to maintain protection coordination in distribution system," in
Innovative Smart Grid Technologies Conference Europe (ISGT-Europe),
2014 IEEE PES, 2014, pp. 1-6: IEEE.
[44] S. Shen et al., "An adaptive protection scheme for distribution systems with
DGs based on optimized Thevenin equivalent parameters estimation," IEEE
Transactions on Power Delivery, vol. 32, no. 1, pp. 411-419, 2017.
[45] D. S. Kumar, D. Srinivasan, and T. Reindl, "A fast and scalable protection
scheme for distribution networks with distributed generation," IEEE
Transactions on Power Delivery, vol. 31, no. 1, pp. 67-75, 2016.
[46] F. B. Costa, A. Monti, and S. C. Paiva, "Overcurrent protection in distribution
systems with distributed generation based on the real-time boundary wavelet
transform," IEEE Transactions on Power Delivery, vol. 32, no. 1, pp. 462-
473, 2017.
[47] J. Buse, D. Shi, T. Ji, and Q. Wu, "Decaying DC offset removal operator
using mathematical morphology for phasor measurement," in Innovative
Smart Grid Technologies Conference Europe (ISGT Europe), 2010 IEEE
PES, 2010, pp. 1-6: IEEE.
[48] Y. Guo, M. Kezunovic, and D. Chen, "Simplified algorithms for removal of
the effect of exponentially decaying DC-offset on the Fourier algorithm,"
IEEE Transactions on power delivery, vol. 18, no. 3, pp. 711-717, 2003.
[49] J.-C. Gu and S.-L. Yu, "Removal of DC offset in current and voltage signals
using a novel Fourier filter algorithm," IEEE Transactions on Power
Delivery, vol. 15, no. 1, pp. 73-79, 2000.
Bibliography 207
[50] A. Rahmati and R. Adhami, "An accurate filtering technique to mitigate
transient decaying DC offset," IEEE Transactions on Power Delivery, vol.
29, no. 2, pp. 966-968, 2014.
[51] V. Torres, J. Guardado, H. Ruiz, and S. Maximov, "Modeling and detection
of high impedance faults," International Journal of Electrical Power &
Energy Systems, vol. 61, pp. 163-172, 2014.
[52] A. Ghaderi, H. L. Ginn III, and H. A. Mohammadpour, "High impedance
fault detection: A review," Electric Power Systems Research, vol. 143, pp.
376-388, 2017.
[53] N. I. Elkalashy, M. Lehtonen, H. A. Darwish, M. A. Izzularab, and I. T.
Abdel-maksoud, "Modeling and experimental verification of high impedance
arcing fault in medium voltage networks," IEEE Transactions on Dielectrics
and Electrical Insulation, vol. 14, no. 2, 2007.
[54] C. Kim, B. D. Russell, and K. Watson, "A parameter-based process for
selecting high impedance fault detection techniques using decision making
under incomplete knowledge," IEEE Transactions on Power Delivery, vol. 5,
no. 3, pp. 1314-1320, 1990.
[55] A.-R. Sedighi, M.-R. Haghifam, O. Malik, and M.-H. Ghassemian, "High
impedance fault detection based on wavelet transform and statistical pattern
recognition," IEEE Transactions on Power Delivery, vol. 20, no. 4, pp. 2414-
2421, 2005.
[56] A. Emanuel, D. Cyganski, J. Orr, S. Shiller, and E. Gulachenski, "High
impedance fault arcing on sandy soil in 15 kV distribution feeders:
contributions to the evaluation of the low frequency spectrum," IEEE
Transactions on Power Delivery, vol. 5, no. 2, pp. 676-686, 1990.
[57] N. Zamanan and J. K. Sykulski, "Modelling arcing high impedances faults in
relation to the physical processes in the electric arc," WSEAS Transactions on
power systems, vol. 1, no. 8, pp. 1507-1512, 2006.
[58] A. Sultan, G. Swift, and D. Fedirchuk, "Detecting arcing downed-wires using
fault current flicker and half-cycle asymmetry," IEEE Transactions on Power
Delivery, vol. 9, no. 1, pp. 461-470, 1994.
[59] W. Santos, F. Lopes, N. Brito, B. Souza, D. Fernandes Jr, and W. Neves,
"High impedance fault detection and location based on electromagnetic
transient analysis," in International Conference on Power Systems Transients
(IPST2013) em Vancouver, Canadá, Julho 18, 2013, vol. 20.
[60] M. S. Thomas, N. Bhaskar, and A. Prakash, "Voltage based detection method
for high impedance fault in a distribution system," Journal of The Institution
of Engineers (India): Series B, vol. 97, no. 3, pp. 413-423, 2016.
[61] J. Chen, B. Phung, D. Zhang, T. Blackburn, and E. Ambikairajah, "Study on
high impedance fault arcing current characteristics," in Power Engineering
Conference (AUPEC), 2013 Australasian Universities, 2013, pp. 1-6: IEEE.
[62] C. G. Wester, "High impedance fault detection on distribution systems," in
Rural Electric Power Conference, 1998. Papers Presented at the 42nd
Annual Conference, 1998, pp. c5-1: IEEE.
[63] A. Sharaf and S. Abu-Azab, "A smart relaying scheme for high impedance
faults in distribution and utilization networks," in Electrical and Computer
Engineering, 2000 Canadian Conference on, 2000, vol. 2, pp. 740-744:
IEEE.
[64] H. Calhoun, M. T. Bishop, C. H. Eichler, and R. E. Lee, "Development and
testing of an electro-mechanical relay to detect fallen distribution
208 Bibliography
conductors," IEEE Transactions on Power Apparatus and Systems, no. 6, pp.
1643-1650, 1982.
[65] C.-L. Huang, H.-Y. Chu, and M.-T. Chen, "Algorithm comparison for high
impedance fault detection based on staged fault test," IEEE transactions on
power delivery, vol. 3, no. 4, pp. 1427-1435, 1988.
[66] B. Aucoin and B. Russell, "Distribution high impedance fault detection
utilizing high frequency current components," IEEE Transactions on Power
Apparatus and Systems, no. 6, pp. 1596-1606, 1982.
[67] S.-J. Huang and C.-T. Hsieh, "High-impedance fault detection utilizing a
Morlet wavelet transform approach," IEEE Transactions on Power Delivery,
vol. 14, no. 4, pp. 1401-1410, 1999.
[68] Y. Sheng and S. M. Rovnyak, "Decision tree-based methodology for high
impedance fault detection," IEEE Transactions on Power Delivery, vol. 19,
no. 2, pp. 533-536, 2004.
[69] S. Samantaray and P. Dash, "High impedance fault detection in distribution
feeders using extended kalman filter and support vector machine,"
International Transactions on Electrical Energy Systems, vol. 20, no. 3, pp.
382-393, 2010.
[70] M. Eissa, G. A. Sowilam, and A. Sharaf, "A new protection detection
technique for high impedance fault using neural network," in Power
Engineering, 2006 Large Engineering Systems Conference on, 2006, pp. 146-
151: IEEE.
[71] I. Baqui, I. Zamora, J. Mazón, and G. Buigues, "High impedance fault
detection methodology using wavelet transform and artificial neural
networks," Electric Power Systems Research, vol. 81, no. 7, pp. 1325-1333,
2011.
[72] S. Samantaray, B. Panigrahi, and P. Dash, "High impedance fault detection in
power distribution networks using time–frequency transform and
probabilistic neural network," IET generation, transmission & distribution,
vol. 2, no. 2, pp. 261-270, 2008.
[73] A. Etemadi and M. Sanaye-Pasand, "High-impedance fault detection using
multi-resolution signal decomposition and adaptive neural fuzzy inference
system," IET generation, transmission & distribution, vol. 2, no. 1, pp. 110-
118, 2008.
[74] M.-R. Haghifam, A.-R. Sedighi, and O. Malik, "Development of a fuzzy
inference system based on genetic algorithm for high-impedance fault
detection," IEE Proceedings-Generation, Transmission and Distribution, vol.
153, no. 3, pp. 359-367, 2006.
[75] A. Bakar, M. Ali, C. Tan, H. Mokhlis, H. Arof, and H. Illias, "High
impedance fault location in 11 kV underground distribution systems using
wavelet transforms," International Journal of Electrical Power & Energy
Systems, vol. 55, pp. 723-730, 2014.
[76] A. Mahari and H. Seyedi, "High impedance fault protection in transmission
lines using a WPT-based algorithm," International Journal of Electrical
Power & Energy Systems, vol. 67, pp. 537-545, 2015.
[77] F. B. Costa, B. Souza, N. Brito, J. Silva, and W. Santos, "Real-time detection
of transients induced by high-impedance faults based on the boundary
wavelet transform," IEEE Transactions on Industry Applications, vol. 51, no.
6, pp. 5312-5323, 2015.
Bibliography 209
[78] S. Gautam and S. M. Brahma, "Detection of high impedance fault in power
distribution systems using mathematical morphology," IEEE Transactions on
Power Systems, vol. 28, no. 2, pp. 1226-1234, 2013.
[79] T. Gammon and J. Matthews, "Instantaneous arcing-fault models developed
for building system analysis," IEEE Transactions on Industry Applications,
vol. 37, no. 1, pp. 197-203, 2001.
[80] M. Jannati and L. Eslami, "Precise modeling of high impedance faults in
power distribution system in emtpworks software," Journal of Electrical
Engineering, 2013.
[81] W. Santos, N. Brito, B. Souza, and F. Pereira, "High impedance fault
modeling based on statistical data," in International Conference on Power
Systems Transients (IPST2015), 2015, pp. 1-8.
[82] J. Carr, "Detection of high impedance faults on multi-grounded primary
distribution systems," IEEE Transactions on Power Apparatus and Systems,
no. 4, pp. 2008-2016, 1981.
[83] L. E. Fisher, "Resistance of low-voltage ac arcs," IEEE Transactions on
Industry and General Applications, no. 6, pp. 607-616, 1970.
[84] M. Michalik, W. Rebizant, M. Lukowicz, S.-J. Lee, and S.-H. Kang, "High-
impedance fault detection in distribution networks with use of wavelet-based
algorithm," IEEE Transactions on Power Delivery, vol. 21, no. 4, pp. 1793-
1802, 2006.
[85] V. Torres, S. Maximov, H. Ruiz, and J. Guardado, "Distributed parameters
model for high-impedance fault detection and localization in transmission
lines," Electric Power Components and Systems, vol. 41, no. 14, pp. 1311-
1333, 2013.
[86] G. Vijayachandran and B. K. Mathew, "High impedance arcing fault
detection in MV networks using discrete wavelet transform and Artificial
Neural Networks," in Green Technologies (ICGT), 2012 International
Conference on, 2012, pp. 089-098: IEEE.
[87] M.-T. Yang, J.-C. Gu, J.-L. Guan, and C.-Y. Cheng, "Evaluation of
algorithms for high impedance faults identification based on staged fault
tests," in Power Engineering Society General Meeting, 2006. IEEE, 2006, p.
8 pp.: IEEE.
[88] N. Zamanan and J. Sykulski, "The evolution of high impedance fault
modeling," in Harmonics and Quality of Power (ICHQP), 2014 IEEE 16th
International Conference on, 2014, pp. 77-81: IEEE.
[89] T. Lai, L. Snider, E. Lo, and D. Sutanto, "High-impedance fault detection
using discrete wavelet transform and frequency range and RMS conversion,"
IEEE Transactions on Power Delivery, vol. 20, no. 1, pp. 397-407, 2005.
[90] Y. Zhao, B. Lehman, J.-F. de Palma, J. Mosesian, and R. Lyons, "Challenges
to overcurrent protection devices under line-line faults in solar photovoltaic
arrays," in Energy Conversion Congress and Exposition (ECCE), 2011 IEEE,
2011, pp. 20-27: IEEE.
[91] J. C. Wiles, "Photovoltaic system grounding," Southwest Technology
Development Institute College of Engineering New Mexico State University
October, 2012.
[92] O. Hegazy, J. Van Mierlo, and P. Lataire, "Analysis, modeling, and
implementation of a multidevice interleaved DC/DC converter for fuel cell
hybrid electric vehicles," IEEE transactions on power electronics, vol. 27,
no. 11, pp. 4445-4458, 2012.
210 Bibliography
[93] J. Flicker and J. Johnson, "Electrical simulations of series and parallel PV
arc-faults," in Photovoltaic Specialists Conference (PVSC), 2013 IEEE 39th,
2013, pp. 3165-3172: IEEE.
[94] X. Yao, L. Herrera, and J. Wang, "Impact evaluation of series dc arc faults in
dc microgrids," in Applied Power Electronics Conference and Exposition
(APEC), 2015 IEEE, 2015, pp. 2953-2958: IEEE.
[95] A. Tjahjono et al., "Adaptive modified firefly algorithm for optimal
coordination of overcurrent relays," IET Generation, Transmission &
Distribution, vol. 11, no. 10, pp. 2575-2585, 2017.
[96] S. Lu, B. Phung, and D. Zhang, "Study on DC series arc fault in photovoltaic
systems for condition monitoring purpose," in Universities Power
Engineering Conference (AUPEC), 2017 Australasian, 2017, pp. 1-6: IEEE.
[97] J. Johnson et al., "Differentiating series and parallel photovoltaic arc-faults,"
in Photovoltaic Specialists Conference (PVSC), 2012 38th IEEE, 2012, pp.
000720-000726: IEEE.
[98] D. A. Dini, P. W. Brazis, and K.-H. Yen, "Development of arc-fault circuit-
interrupter requirements for photovoltaic systems," in Photovoltaic
Specialists Conference (PVSC), 2011 37th IEEE, 2011, pp. 001790-001794:
IEEE.
[99] A. Dhoke and A. Mengede, "Challenges to overcurrent protection devices in
PV array during winter and low irradiation conditions in Australia," in
Innovative Smart Grid Technologies-Asia (ISGT-Asia), 2017 IEEE, 2017, pp.
1-6: IEEE.
[100] Y. Zhao, J.-F. De Palma, J. Mosesian, R. Lyons, and B. Lehman, "Line–line
fault analysis and protection challenges in solar photovoltaic arrays," IEEE
transactions on Industrial Electronics, vol. 60, no. 9, pp. 3784-3795, 2013.
[101] F. Schimpf and L. E. Norum, "Recognition of electric arcing in the DC-
wiring of photovoltaic systems," in Telecommunications Energy Conference,
2009. INTELEC 2009. 31st International, 2009, pp. 1-6: IEEE.
[102] F. M. Uriarte et al., "A DC arc model for series faults in low voltage
microgrids," IEEE Transactions on Smart Grid, vol. 3, no. 4, pp. 2063-2070,
2012.
[103] G.-S. Seo, K. A. Kim, K.-C. Lee, K.-J. Lee, and B.-H. Cho, "A new DC arc
fault detection method using DC system component modeling and analysis in
low frequency range," in Applied Power Electronics Conference and
Exposition (APEC), 2015 IEEE, 2015, pp. 2438-2444: IEEE.
[104] M. Naidu, T. J. Schoepf, and S. Gopalakrishnan, "Arc fault detection scheme
for 42-V automotive DC networks using current shunt," IEEE Transactions
on power electronics, vol. 21, no. 3, pp. 633-639, 2006.
[105] H. Hamilton and N. N. Schulz, "DC protection on the electric ship," in
Electric Ship Technologies Symposium, 2007. ESTS'07. IEEE, 2007, pp. 294-
300: IEEE.
[106] J. A. Momoh and R. Button, "Design and analysis of aerospace DC arcing
faults using fast fourier transformation and artificial neural network," in
Power Engineering Society General Meeting, 2003, IEEE, 2003, vol. 2, pp.
788-793: IEEE.
[107] X. Yao, L. Herrera, S. Ji, K. Zou, and J. Wang, "Characteristic study and
time-domain discrete-wavelet-transform based hybrid detection of series DC
arc faults," IEEE Transactions on Power Electronics, vol. 29, no. 6, pp.
3103-3115, 2014.
Bibliography 211
[108] R. D. Telford, S. Galloway, B. Stephen, and I. Elders, "Diagnosis of Series
DC Arc Faults—A Machine Learning Approach," IEEE Transactions on
Industrial Informatics, vol. 13, no. 4, pp. 1598-1609, 2017.
[109] T. Shimakage, K. Nishioka, H. Yamane, M. Nagura, and M. Kudo,
"Development of fault detection system in PV system," in
Telecommunications Energy Conference (INTELEC), 2011 IEEE 33rd
International, 2011, pp. 1-5: IEEE.
[110] R. Platon, J. Martel, N. Woodruff, and T. Y. Chau, "Online fault detection in
PV systems," IEEE Transactions on Sustainable Energy, vol. 6, no. 4, pp.
1200-1207, 2015.
[111] P. Ducange, M. Fazzolari, B. Lazzerini, and F. Marcelloni, "An intelligent
system for detecting faults in photovoltaic fields," in Intelligent systems
design and applications (ISDA), 2011 11th international conference on,
2011, pp. 1341-1346: IEEE.
[112] A. Chouder and S. Silvestre, "Automatic supervision and fault detection of
PV systems based on power losses analysis," Energy conversion and
Management, vol. 51, no. 10, pp. 1929-1937, 2010.
[113] S. Silvestre, A. Chouder, and E. Karatepe, "Automatic fault detection in grid
connected PV systems," Solar Energy, vol. 94, pp. 119-127, 2013.
[114] S. Spataru, D. Sera, T. Kerekes, and R. Teodorescu, "Photovoltaic array
condition monitoring based on online regression of performance model," in
Photovoltaic Specialists Conference (PVSC), 2013 IEEE 39th, 2013, pp.
0815-0820: IEEE.
[115] S. Vergura, G. Acciani, V. Amoruso, G. E. Patrono, and F. Vacca,
"Descriptive and inferential statistics for supervising and monitoring the
operation of PV plants," IEEE Transactions on Industrial Electronics, vol.
56, no. 11, pp. 4456-4464, 2009.
[116] B. Andò, S. Baglio, A. Pistorio, G. M. Tina, and C. Ventura, "Sentinella:
Smart monitoring of photovoltaic systems at panel level," IEEE Transactions
on Instrumentation and Measurement, vol. 64, no. 8, pp. 2188-2199, 2015.
[117] P. Guerriero, F. Di Napoli, G. Vallone, V. d'Alessandro, and S. Daliento,
"Monitoring and diagnostics of PV plants by a wireless self-powered sensor
for individual panels," IEEE Journal of Photovoltaics, vol. 6, no. 1, pp. 286-
294, 2016.
[118] Z. Yi and A. H. Etemadi, "Fault detection for photovoltaic systems based on
multi-resolution signal decomposition and fuzzy inference systems," IEEE
Transactions on Smart Grid, vol. 8, no. 3, pp. 1274-1283, 2017.
[119] C.-L. Kuo, J.-L. Chen, S.-J. Chen, C.-C. Kao, H.-T. Yau, and C.-H. Lin,
"Photovoltaic energy conversion system fault detection using fractional-order
color relation classifier in microdistribution systems," IEEE Transactions on
Smart Grid, vol. 8, no. 3, pp. 1163-1172, 2017.
[120] S. Chae, J. Park, and S. Oh, "Series DC arc fault detection algorithm for DC
microgrids using relative magnitude comparison," IEEE Journal of Emerging
and Selected Topics in Power Electronics, vol. 4, no. 4, pp. 1270-1278, 2016.
[121] H. Haeberlin and M. Real, "Arc detector for remote detection of dangerous
arcs on the DC side of PV plants," in 22nd European Photovoltaic Solar
Energy Conference, Milano, Italy, 2007, vol. 200.
[122] C. Strobl and P. Meckler, "Arc faults in photovoltaic systems," in 2010
Proceedings of the 56th IEEE Holm Conference on Electrical Contacts,
2010, pp. 1-7.
212 Bibliography
[123] Z. Wang and R. S. Balog, "Arc fault and flash detection in DC photovoltaic
arrays using wavelets," in Photovoltaic Specialists Conference (PVSC), 2013
IEEE 39th, 2013, pp. 1619-1624: IEEE.
[124] R. F. Ammerman, T. Gammon, P. K. Sen, and J. P. Nelson, "DC arc models
and incident energy calculations," in Petroleum and Chemical Industry
Conference, 2009. PCIC 2009. 2009 Record of Conference Papers-Industry
Applications Society 56th Annual, 2009, pp. 1-13: IEEE.
[125] X. Yao, J. Wang, and D. L. Schweickart, "Review and recent developments
in DC arc fault detection," in Power Modulator and High Voltage Conference
(IPMHVC), 2016 IEEE International, 2016, pp. 467-472: IEEE.
[126] A. Stokes and W. Oppenlander, "Electric arcs in open air," Journal of Physics
D: Applied Physics, vol. 24, no. 1, p. 26, 1991.
[127] S. Samantaray, "Ensemble decision trees for high impedance fault detection
in power distribution network," International Journal of Electrical Power &
Energy Systems, vol. 43, no. 1, pp. 1048-1055, 2012.
[128] M. Banejad and H. Ijadi, "High impedance fault detection: Discrete wavelet
transform and fuzzy function approximation," Journal of AI and Data
Mining, vol. 2, no. 2, pp. 149-158, 2014.
[129] N. Varma, D. Ram, and D. Anjaneyulu, "Development of fault detection
algorithm for high impedance faults in distribution network using
multiresolution analysis," International Journal of Engineering Research and
Technology, vol. 3, no. 9, pp. 573-576, 2014.
[130] G. Matheron, G. Matheron, G. Matheron, and G. Matheron, "Random sets
and integral geometry," 1975.
[131] J. Serra, Image analysis and mathematical morphology. Academic Press, Inc.,
1983.
[132] Y. Tingfang, L. Pei, Z. Xiangjun, and K. Li, "Application of adaptive
generalized morphological filter in disturbance identification for power
system signatures," in Power System Technology, 2006. PowerCon 2006.
International Conference on, 2006, pp. 1-7: IEEE.
[133] O. Sen and R. Zhen, "Application of improved mathematical morphology
method in the power quality monitoring," in Power System Technology, 2006.
PowerCon 2006. International Conference on, 2006, pp. 1-6: IEEE.
[134] Z. Lu, D. Turner, Q. Wu, J. Fitch, and S. Mann, "Morphological transform
for detection of power quality disturbances," in Power System Technology,
2004. PowerCon 2004. 2004 International Conference on, 2004, vol. 2, pp.
1644-1649: IEEE.
[135] G. Li, M. Zhou, Y. Luo, and Y. Ni, "Power quality disturbance detection
based on mathematical morphology and fractal technique," in Transmission
and Distribution Conference and Exhibition: Asia and Pacific, 2005
IEEE/PES, 2005, pp. 1-6: IEEE.
[136] P. M. Ramos, F. M. Janeiro, and A. C. Serra, "PQ monitoring system for real-
time detection and classification of disturbances in a single-phase power
system," IEEE Transactions on Instrumentation and Measurement, vol. 57,
no. 8, pp. 1725-1733, 2008.
[137] L.-a. Chen and P. Zhang, "Detection and protection of short circuit fault
based on morphology-wavelet," in Transmission and Distribution Conference
and Exhibition: Asia and Pacific, 2005 IEEE/PES, 2005, pp. 1-5: IEEE.
Bibliography 213
[138] Z. Lu, Q. Wu, and J. Fitch, "A morphological filter for estimation of power
system harmonics," in Power System Technology, 2006. PowerCon 2006.
International Conference on, 2006, pp. 1-5: IEEE.
[139] X. Lin, L. Zou, Q. Tian, H. Weng, and P. Liu, "A series multiresolution
morphological gradient-based criterion to identify CT saturation," IEEE
transactions on power delivery, vol. 21, no. 3, pp. 1169-1175, 2006.
[140] Z. Lu, J. S. Smith, and Q. H. Wu, "Morphological lifting scheme for current
transformer saturation detection and compensation," IEEE Transactions on
Circuits and Systems I: Regular Papers, vol. 55, no. 10, pp. 3349-3357, 2008.
[141] C. Zhang, X. Lin, Z. Chen, Z. Bo, and M. Redfern, "A novel mathematical
morphology based anti-CT saturation criterion for generator differential
protection," in Universities Power Engineering Conference, 2004. UPEC
2004. 39th International, 2004, vol. 2, pp. 763-767: IEEE.
[142] M. Jing, W. Zengping, X. Yan, and W. Jie, "A novel adaptive algorithm to
identify inrush using mathematical morphology," in Power Systems
Conference and Exposition, 2006. PSCE'06. 2006 IEEE PES, 2006, pp. 1020-
1028: IEEE.
[143] M. Jing, X. Yan, W. Zengping, and L. Haofang, "A novel adaptive scheme of
discrimination between internal faults and inrush currents of transformer
using mathematical morphology," in Power Engineering Society General
Meeting, 2006. IEEE, 2006, p. 7 pp.: IEEE.
[144] L. Zou, P. Liu, and Q. Zhao, "Applications of multiresolution morphological
analysis in ultra high speed protection of transmission line," IEE
Proceedings-Generation, Transmission and Distribution, vol. 153, no. 5, pp.
515-523, 2006.
[145] Q. Wu, J. Zhang, and D. Zhang, "Ultra-high-speed directional protection of
transmission lines using mathematical morphology," IEEE transactions on
power delivery, vol. 18, no. 4, pp. 1127-1133, 2003.
[146] M. Jing, W. Zengping, X. Yan, and M. Lei, "Single-ended transient positional
protection of transmission lines using mathematical morphology," in Power
Engineering Conference, 2005. IPEC 2005. The 7th International, 2005, pp.
1-603: IEEE.
[147] L. Zou, P. Liu, and Q. Zhao, "Mathematical morphology based phase
selection scheme in digital relaying," IEE Proceedings-Generation,
Transmission and Distribution, vol. 152, no. 2, pp. 157-163, 2005.
[148] X. Lin, Z. Bo, P. Liu, and M. Redfern, "A novel directional protective relay
using mathematic morphology," in Power Engineering Society General
Meeting, 2003, IEEE, 2003, vol. 3, pp. 1854-1861: IEEE.
[149] F. Y. Shih, Image processing and mathematical morphology: fundamentals
and applications. CRC press, 2009.
[150] Q.-H. Wu, Z. Lu, and T. Ji, Protective relaying of power systems using
mathematical morphology. Springer Science & Business Media, 2009.
[151] E. R. Dougherty and R. A. Lotufo, Hands-on morphological image
processing. SPIE press, 2003.
[152] S. Gautam and S. M. Brahma, "Application of mathematical morphology
based filters to detect a power swing," in Power and Energy Society General
Meeting, 2010 IEEE, 2010, pp. 1-6: IEEE.
[153] A. G. Phadke and J. S. Thorp, Synchronized phasor measurements and their
applications. Springer, 2017.
214 Bibliography
[154] M. Yang and S. Tan, "Research on an ultra-high-speed protection based on
multiresolution morphological gradient," in Power System Technology, 2006.
PowerCon 2006. International Conference on, 2006, pp. 1-5: IEEE.
[155] L. Zhang, Y. ping Lv, and H. Yin, "A multi-resolution morphology gradient
based non-communication protection scheme for transmission lines," in
Electric Utility Deregulation and Restructuring and Power Technologies,
2008. DRPT 2008. Third International Conference on, 2008, pp. 1608-1613:
IEEE.
[156] D. Zhang, Q. Li, J. Zhang, Q. Wu, and D. Turner, "Improving the accuracy of
single-ended transient fault locators using mathematical morphology," in
Power System Technology, 2002. Proceedings. PowerCon 2002.
International Conference on, 2002, vol. 2, pp. 788-792: IEEE.
[157] S. Gautam and S. M. Brahma, "Guidelines for selection of an optimal
structuring element for mathematical morphology based tools to detect power
system disturbances," in Power and Energy Society General Meeting, 2012
IEEE, 2012, pp. 1-6: IEEE.
[158] F. Y.-C. Shih and O. R. Mitchell, "Decomposition of gray-scale
morphological structuring elements," Pattern recognition, vol. 24, no. 3, pp.
195-203, 1991.
[159] H.-x. Zhang, J.-y. Li, M.-z. Wang, and H. Li, "Image edge detection based on
fusion of wavelet transform and mathematical morphology," in Computer
Science & Education (ICCSE), 2016 11th International Conference on, 2016,
pp. 981-984: IEEE.
[160] W. Zhang, D. Shi, and X. Yang, "An improved edge detection algorithm
based on mathematical morphology and directional wavelet transform," in
Image and Signal Processing (CISP), 2015 8th International Congress on,
2015, pp. 335-339: IEEE.
[161] S. M. Brahma, "Fault location in power distribution system with penetration
of distributed generation," IEEE transactions on power delivery, vol. 26, no.
3, pp. 1545-1553, 2011.
[162] S.-H. Hwang, H.-T. Chen, and C.-T. Chang, "An exponentially weighted
moving average method for identification and monitoring of stochastic
systems," Industrial & Engineering Chemistry Research, vol. 47, no. 21, pp.
8239-8249, 2008.
[163] S. Shome, S. Vadali, U. Datta, S. Sen, and A. Mukherjee, "Performance
evaluation of different averaging based filter designs using digital signal
processor and its synthesis on FPGA," International Journal of Signal
Processing, Image Processing and Pattern Recognition, vol. 5, no. 3, pp. 75-
92, 2012.
[164] L. Ljung, "System identification," in Signal analysis and prediction:
Springer, 1998, pp. 163-173.
[165] M. Sachdev and M. Nagpal, "A recursive least error squares algorithm for
power system relaying and measurement applications," IEEE Transactions on
Power Delivery, vol. 6, no. 3, pp. 1008-1015, 1991.
[166] F. Costa, D. Formiga, R. Ferreira, T. Sousa, and F. Costa, "A recursive least-
squares aided by pre-filtering for phasor-estimation in distance protection," in
PowerTech (POWERTECH), 2013 IEEE Grenoble, 2013, pp. 1-6: IEEE.
[167] T. M. MATLAB and Simscape Toolbox Release 2016a, Inc, Natick,
Massachusetts, USA.
Bibliography 215
[168] J. P. Ram, T. S. Babu, and N. Rajasekar, "A comprehensive review on solar
PV maximum power point tracking techniques," Renewable and Sustainable
Energy Reviews, vol. 67, pp. 826-847, 2017.
[169] K. Ishaque, Z. Salam, and H. Taheri, "Simple, fast and accurate two-diode
model for photovoltaic modules," Solar energy materials and solar cells, vol.
95, no. 2, pp. 586-594, 2011.
[170] P. Sanjeevikumar, G. Grandi, P. W. Wheeler, F. Blaabjerg, and J. Loncarski,
"A simple MPPT algorithm for novel PV power generation system by high
output voltage DC-DC boost converter," in Industrial Electronics (ISIE),
2015 IEEE 24th International Symposium on, 2015, pp. 214-220: IEEE.
[171] R. W. Erickson and D. Maksimovic, Fundamentals of power electronics.
Springer Science & Business Media, 2007.
[172] H. Ibrahim and N. Anani, "Variations of PV module parameters with
irradiance and temperature," Energy Procedia, vol. 134, pp. 276-285, 2017.
[173] S. Xiao and R. S. Balog, "An improved adaptive perturb & observe maximum
power point tracking technique," in Texas Power and Energy Conference
(TPEC), 2018 IEEE, 2018, pp. 1-6: IEEE.
[174] J. Chauhan, P. Chauhan, T. Maniar, and A. Joshi, "Comparison of MPPT
algorithms for DC-DC converters based photovoltaic systems," in Energy
Efficient Technologies for Sustainability (ICEETS), 2013 International
Conference on, 2013, pp. 476-481: IEEE.
[175] C. Hua, J. Lin, and C. Shen, "Implementation of a DSP-controlled
photovoltaic system with peak power tracking," IEEE Transactions on
Industrial Electronics, vol. 45, no. 1, pp. 99-107, 1998.
[176] K. Hussein, I. Muta, T. Hoshino, and M. Osakada, "Maximum photovoltaic
power tracking: an algorithm for rapidly changing atmospheric conditions,"
IEE Proceedings-Generation, Transmission and Distribution, vol. 142, no. 1,
pp. 59-64, 1995.
[177] A. K. Abdelsalam, A. M. Massoud, S. Ahmed, and P. N. Enjeti, "High-
performance adaptive perturb and observe MPPT technique for photovoltaic-
based microgrids," IEEE Transactions on Power Electronics, vol. 26, no. 4,
pp. 1010-1021, 2011.
[178] J. Ahmed and Z. Salam, "A modified P&O maximum power point tracking
method with reduced steady-state oscillation and improved tracking
efficiency," IEEE Transactions on Sustainable Energy, vol. 7, no. 4, pp.
1506-1515, 2016.
[179] O. Anaya-Lara, D. Campos-Gaona, E. Moreno-Goytia, and G. Adam,
Offshore wind energy generation: control, protection, and integration to
electrical systems. John Wiley & Sons, 2014.
[180] W. H. Kersting, "Radial distribution test feeders," in Power Engineering
Society Winter Meeting, 2001. IEEE, 2001, vol. 2, pp. 908-912: IEEE.
[181] A. Sedighi and M. Haghifam, "Simulation of high impedance ground fault In
electrical power distribution systems," in Power System Technology
(POWERCON), 2010 International Conference on, 2010, pp. 1-7: IEEE.
[182] B. M. Aucoin and R. H. Jones, "High impedance fault detection
implementation issues," IEEE Transactions on Power Delivery, vol. 11, no.
1, pp. 139-148, 1996.
[183] H. Sharaf, H. Zeineldin, and E. El-Saadany, "Protection coordination for
microgrids with grid-connected and islanded capabilities using
216 Bibliography
communication assisted dual setting directional overcurrent relays," IEEE
Trans. Smart Grid, vol. 9, pp. 143-151, 2018.
[184] V. Telukunta, J. Pradhan, A. Agrawal, M. Singh, and S. G. Srivani,
"Protection challenges under bulk penetration of renewable energy resources
in power systems: A review," CSEE Journal of Power and Energy Systems,
vol. 3, no. 4, pp. 365-379, 2017.
[185] D. M. Vilathgamuwa, P. C. Loh, and Y. Li, "Protection of microgrids during
utility voltage sags," IEEE Transactions on Industrial Electronics, vol. 53,
no. 5, pp. 1427-1436, 2006.
[186] M. Ebrahimpour, B. Vahidi, and S. H. Hosseinian, "A hybrid
superconducting fault current controller for DG networks and microgrids,"
IEEE Transactions on Applied Superconductivity, vol. 23, no. 5, pp.
5604306-5604306, 2013.
[187] P. Naveen and P. Jena, "A Review on Issues and Coordination Strategies for
Over Current Protection in Microgrid," in 2017 14th IEEE India Council
International Conference (INDICON), 2017, pp. 1-6: IEEE.
[188] M. Y. Shih, A. Conde, Z. Leonowicz, and L. Martirano, "An adaptive
overcurrent coordination scheme to improve relay sensitivity and overcome
drawbacks due to distributed generation in smart grids," IEEE Transactions
on Industry Applications, vol. 53, no. 6, pp. 5217-5228, 2017.
[189] E. Dehghanpour, H. K. Karegar, R. Kheirollahi, and T. Soleymani, "Optimal
coordination of directional overcurrent relays in microgrids by using cuckoo-
linear optimization algorithm and fault current limiter," IEEE Transactions
on Smart Grid, vol. 9, no. 2, pp. 1365-1375, 2018.
[190] H. M. Zeineldin, H. H. Sharaf, and E. El-Saadany, "Protection coordination
for microgrids with grid-connected and islanded capabilities using dual
setting directional overcurrent relays," IEEE Transactions on Smart Grid,
2016.
[191] W. J. Park, B. C. Sung, K.-B. Song, and J.-W. Park, "Parameter optimization
of SFCL with wind-turbine generation system based on its protective
coordination," IEEE Transactions on Applied Superconductivity, vol. 21, no.
3, p. 2153, 2011.
[192] D. K. Ibrahim, E. E. D. A. El Zahab, and S. A. E. A. Mostafa, "New
coordination approach to minimize the number of re-adjusted relays when
adding DGs in interconnected power systems with a minimum value of fault
current limiter," International Journal of Electrical Power & Energy Systems,
vol. 85, pp. 32-41, 2017.
[193] V. R. Pandi, H. Zeineldin, and W. Xiao, "Determining optimal location and
size of distributed generation resources considering harmonic and protection
coordination limits," IEEE transactions on power systems, vol. 28, no. 2, pp.
1245-1254, 2013.
[194] H. Zhan et al., "Relay protection coordination integrated optimal placement
and sizing of distributed generation sources in distribution networks," IEEE
Transactions on Smart Grid, vol. 7, no. 1, pp. 55-65, 2016.
[195] P. Mahat, Z. Chen, B. Bak-Jensen, and C. L. Bak, "A Simple Adaptive
Overcurrent Protection of Distribution Systems With Distributed
Generation," IEEE Trans. Smart Grid, vol. 2, no. 3, pp. 428-437, 2011.
[196] J. Ma, X. Wang, Y. Zhang, Q. Yang, and A. Phadke, "A novel adaptive
current protection scheme for distribution systems with distributed
Bibliography 217
generation," International Journal of Electrical Power & Energy Systems,
vol. 43, no. 1, pp. 1460-1466, 2012.
[197] M. Ojaghi, Z. Sudi, and J. Faiz, "Implementation of full adaptive technique to
optimal coordination of overcurrent relays," IEEE Transactions on Power
Delivery, vol. 28, no. 1, pp. 235-244, 2013.
[198] J.-D. Park and J. Candelaria, "Fault detection and isolation in low-voltage
DC-bus microgrid system," IEEE transactions on power delivery, vol. 28, no.
2, pp. 779-787, 2013.
[199] L. ZHANG, T. Nengling, W. HUANG, L. Jian, and W. Yanhong, "A review
on protection of DC microgrids," Journal of Modern Power Systems and
Clean Energy, pp. 1-15, 2018.
[200] N. Bayati, A. Hajizadeh, and M. Soltani, "Protection in DC microgrids: a
comparative review," IET Smart Grid, vol. 1, no. 3, pp. 66-75, 2018.
[201] D. Salomonsson, L. Soder, and A. Sannino, "Protection of low-voltage DC
microgrids," IEEE Transactions on Power Delivery, vol. 24, no. 3, pp. 1045-
1053, 2009.
[202] D. M. Bui, S.-L. Chen, C.-H. Wu, K.-Y. Lien, C.-H. Huang, and K.-K. Jen,
"Review on protection coordination strategies and development of an
effective protection coordination system for DC microgrid," in Power and
Energy Engineering Conference (APPEEC), 2014 IEEE PES Asia-Pacific,
2014, pp. 1-10: IEEE.
[203] S. D. Fletcher, P. J. Norman, S. J. Galloway, P. Crolla, and G. M. Burt,
"Optimizing the roles of unit and non-unit protection methods within DC
microgrids," IEEE transactions on Smart Grid, vol. 3, no. 4, pp. 2079-2087,
2012.
[204] B. Devarajan, "Modeling of a DC fuse for protection of semiconductor
devices using PSCAD/EMTDC," 2010.
[205] N. Sadiq and T. Noor, "Selecting fuses: simple procedures to get the right
overcurrent protection for DC–DC converters," Power Electronics
Technology, 2010.
[206] J. Candelaria and J.-D. Park, "VSC-HVDC system protection: A review of
current methods," in Power Systems Conference and Exposition (PSCE),
2011 IEEE/PES, 2011, pp. 1-7: IEEE.
[207] C. M. Franck, "HVDC circuit breakers: A review identifying future research
needs," IEEE Transactions on Power Delivery, vol. 26, no. 2, pp. 998-1007,
2011.
[208] C. Meyer, S. Schroder, and R. W. De Doncker, "Solid-state circuit breakers
and current limiters for medium-voltage systems having distributed power
systems," IEEE transactions on power electronics, vol. 19, no. 5, pp. 1333-
1340, 2004.
[209] D. Lawes, L. Ran, and Z. Xu, "Design of a solid-state DC circuit breaker for
light rail transit power supply network," in Energy Conversion Congress and
Exposition (ECCE), 2014 IEEE, 2014, pp. 350-357: IEEE.
[210] Y. Zhang and Y. C. Liang, "Over-current protection scheme for SiC power
MOSFET DC circuit breaker," in 2014 IEEE Energy Conversion Congress
and Exposition (ECCE), 2014, pp. 1967-1971: IEEE.
[211] M. E. Baran and N. R. Mahajan, "Overcurrent protection on voltage-source-
converter-based multiterminal DC distribution systems," IEEE transactions
on power delivery, vol. 22, no. 1, pp. 406-412, 2007.
218 Bibliography
[212] R. M. Cuzner and G. Venkataramanan, "The status of DC micro-grid
protection," in Industry Applications Society Annual Meeting, 2008. IAS'08.
IEEE, 2008, pp. 1-8: IEEE.
[213] L. Tang and B. Ooi, "Locating and Isolating DC Faults in Multi-Terminal DC
Systems," IEEE Transactions on Power Delivery, vol. 22, no. 3, pp. 1877-
1884, 2007.
[214] M. Monadi, C. Koch-Ciobotaru, A. Luna, J. I. Candela, and P. Rodriguez, "A
protection strategy for fault detection and location for multi-terminal MVDC
distribution systems with renewable energy systems," in 2014 International
Conference on Renewable Energy Research and Application (ICRERA),
2014, pp. 496-501.
[215] J. Park, J. Candelaria, L. Ma, and K. Dunn, "DC Ring-Bus Microgrid Fault
Protection and Identification of Fault Location," IEEE Transactions on
Power Delivery, vol. 28, no. 4, pp. 2574-2584, 2013.
[216] M. Farhadi and O. A. Mohammed, "Event-Based Protection Scheme for a
Multiterminal Hybrid DC Power System," IEEE Transactions on Smart Grid,
vol. 6, no. 4, pp. 1658-1669, 2015.
[217] M. Farhadi and O. A. Mohammed, "A New Protection Scheme for Multi-Bus
DC Power Systems Using an Event Classification Approach," IEEE
Transactions on Industry Applications, vol. 52, no. 4, pp. 2834-2842, 2016.
[218] G. Venkataramanan and C. Marnay, "A larger role for microgrids," IEEE
Power and Energy Magazine, vol. 6, no. 3, pp. 78-82, 2008.
[219] V. K. Sood, D. Fischer, J. M. Eklund, and T. Brown, "Developing a
communication infrastructure for the Smart Grid," in 2009 IEEE Electrical
Power & Energy Conference (EPEC), 2009, pp. 1-7.
Bibliography 219