Digital Distance Relay Modeling and Testing Using LabVIEW and MATLAB Simulink
Digital Distance Relay Modeling and Testing Using LabVIEW and MATLAB Simulink
discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/306106848
CITATIONS
READS
12
2 authors, including:
Ayache Mati
University M'Hamed Bougara of Boumerdes
2 PUBLICATIONS 0 CITATIONS
SEE PROFILE
MASTER
In Electrical and Electronic Engineering
Option: Power Engineering
Title:
Registration Number:..../2015
Dedication
Every challenging work needs self-efforts as well as
guidance of Elders those who were very close to our heart.
My humble effort I dedicate to my sweet and loving
Mother and my family members,
Whose affection, love, encouragement and prays of day and
night make me able to get such success and honor.
Along with all my friends, hardworking and respected
Teachers
Merouane BEGBAGUI
II
II
Dedication
I have a great pleasure to dedicate this modest work
To my Beloved Mother, my Dear Father
To my Dear Sisters, Brothers, Uncles, Aunts and Cousins
To all my Friends
To all my Teachers from primary school to my last year of
university
And to all with whom I spent wonderful moments
Ayache MATI
ACKNOWLEDGEMENT
In the name of Allah, the Most Gracious and the Most Merciful
Alhamdulillah, all praises to Allah for the strengths and His blessing in completing
this project.
We would like to express our deepest and sincere gratitude to our project
Supervisor Pr. H.BENTARZI. It was a great privilege and honor to work and study
under your supervision, we would like also to thank the other teachers for their
precious help during our work. Thank you very much.
Last but not least, we are infinitely grateful to our family members, particularly
our parents for their patience, unwavering support, continuous encouragement,
and belief in us throughout our whole life. We would have never made it this far
without them beside us every step of the way.
III
Abstract
IV
Table of Contents
Dedication
Acknowledgment
III
Abstract
IV
Table of Contents
List of Figures
VII
List of tables
VIII
Introduction
Chapter One
Chapter Two
Chapter Three
Chapter Four
Conclusion....................................................................................................... ..38
References
List of Figures
Figure 1.1
Figure 1.2
Figure 1.3
Figure 1.4
Figure 2.1
Figure 2.2
Figure 2.3
Figure 2.4
Figure 2.5
Figure 3.1
Figure 3.2
Figure 3.3
Figure 3.4
Figure 3.5
Figure 3.6
Figure 3.7
Figure 3.8
Figure 3.9
Figure 3.10
Figure 3.11
Figure 3.12
Figure 3.13
Figure 4.1
Figure 4.2
Figure 4.3
Figure 4.4
Figure 4.5
Figure 4.6
Figure 4.7
Figure 4.8
Figure 4.9
List of tables
Table 1.1
Table 3.1
Table 4.1
VII
Introduction
INTRODUCTION
Today, the major challenging task for an electrical engineer is ensuring a high
level of continuity of service to customers even under system disturbance. However, a
number of undesirable but unavoidable nature events or human-error incidents may
occur and disrupt this condition. The cause of accident includes lightings, wind
damage, ice loading, tree falling, bird shorting, aircraft colliding, vehicles hitting,
people contacting, digging into underground cable, and so on.
To avoid damage to the equipment of the utilities, long interruption service to the
customers and possible personal hazards, proper protective relays are necessary to
take suitable corrective actions during these abnormal conditions. Originally, all
protective relays were electromechanical type, which are still being widely used in
many systems. Solid state relays were introduced in the 1950s and are commonly
used today for their relative accuracy, sensitivity, ease of testing and maintaining.
Recently, researchers have been trying to develop a more reliable, secure and fast
acting relay with small space and power consumption by using microprocessor
technology.
Distance protection is the most widely used method to protect transmission lines. The
fundamental principle of distance Relying is based on the local measurement of voltages and
currents, where the Relay responds to the impedance between the relay terminal and the fault
location. There are many types of distance relay characteristic such as mho, reactance,
admittance, quadrilateral polarized-mho, offset mho etc. Every type of characteristics has
different intended function and theories behind.
In order to understanding the function of relays, software relay models must be realized,
modeling of protective relays offer an economic and feasible alternative to studying the
performance of protective relays. Relay models have been long used in a variety of tasks,
such as designing new relaying algorithms, optimizing relay settings. Electric power utilities
use computer-based relay models to confirm how the relay would perform during systems
disturbances and normal operating conditions and to make the necessary corrective
adjustment on the relay settings.
LabVIEW software for National Instruments has been used as interfacing software.
This makes the modeling process and analysis easier because LabVIEW has many
features and functions that can be used together with data acquisition card from
National Instruments. LabVIEW is a graphical programming language that uses icons
instead of lines of text to create applications. In contrast to the text-based
Page 1
Introduction
programming languages, where instructions determine program execution, LabVIEW
uses data flow programming, where the flow of data determines execution.
The goal of this paper is to explain the building process of LabVIEW model for
distance relay. Inside the modeling, fault detection, apparent impedance calculation
for all types of faults, and zone coordination were designed and implemented. A Mho
type distance characteristic was chosen to be as the protection scheme. For this relay,
the developed model can be included in one block set only by creating the subsystem
for the developed model.
This study is divided into four main chapters:
-Chapter one: presents an overview of distance protection philosophy, and its
important aspects.
-Chapter two: deals with generalities of digital and numerical relays.
-chapter three: illustrates the design model in LabVIEW and tools required for the
proposed distance relay.
-Chapter four: deals with the implementation, simulation, and testing the mho
distance relay developed in LabVIEW.
Page 2
CHAPTER 1
Distance protection philosophy
Introduction
Distance relays
Distance relay zones
Inputs signals to the relay
Distance protection comparators
Distance protection characteristics
Chapter I
Page 3
Chapter I
Chapter I
fault. Once this is achieved, the zone settings of all relays can be based upon the total positive
sequence impedance of the line, regardless of the type of the fault. We will now consider various
types of fault, and give the appropriate voltage and current inputs to be used for the distance
relays responsible for each of these fault types[1,3,5].In the following table we are going to see
fault impedance algorithm for various fault types:
Table 1.1 Fault impedance Algorithm for various fault types
Fault Type
Algorithm
AB or ABG
(VA-VB)/(IA-IB)
(1)
AC or ACG
(VA-VC)/(IA-IC)
(2)
BC or BCG
(VB-VC)/(IB-IC)
(3)
AG
VA/(IA+3K0I0)
(4)
BG
VB/(IB+3K0I0)
(5)
CG
VC/(IC+3K0I0)
(6)
ABC or ABCG
(1)=(2)=(3)
Where:
A, B and C indicates faulty phases, G indicates ground fault.
VA, VB and VC indicate voltage phases
IA, IB and IC indicate current phases
Z0 = line zero-sequence impedance
Z1 = line positive-sequence impedance
K1 = residual compensation factor where k0 = (Z0-Z1)/KZ1. K can be 1 or3 depend on the relay
design.
I0 =1/3( IA + IB + IC)
Page 5
Chapter I
The composite signals in a phase comparator are denoted by S1 and S2. An angular displacement
is considered positive if S1 leads S2. A phase comparator operates if the following condition is
(1.1)
satisfied
| |
(1.2)
(1.4)
Dividing these equations by the line current Ir -r, give the following equations.
(1.5)
(1.6)
Page 6
Chapter I
As seen in Figure 1.3, the impedances S'1 and S'2 are placed in the extremes of the constant ZR z
impedance. When the system impedance Zr r is inside the operating characteristic, as shown on
Figure (1.3(a)), the angle between S'1 and S'2 fulfills equation (1.1) and the relay operates. In
Figure (1.3(b)) is shown the case of Ir r lying outside the operating characteristic. Now, the
angle between S'1and S'2 is outside the range specified in Equation (1.1), and the relay does not
operate. The constant parameter ZR z marks the diameter of the circular characteristic that
passes through the origin.
(1.7)
(1.8)
Dividing these equations by Ir -r leads to the following equations.
(1.9)
(1.10)
When the system impedance Zr r is inside the characteristic, the absolute value of the
impedance SR is less than the absolute value of the radius SO as shown in Figure 1.4(a).
Page 7
Chapter I
the condition specified in Equation (2.2) is satisfied and the relay operates. When the system
impedance Zr r is outside the characteristic, the absolute value of SR is larger than the absolute
value of SO and the relay does not operate as shown in Figure 1.4(b).
Page 8
CHAPTER 2
Digital and numerical relays
Introduction
Relay performance
Relay technology
Generalized Numerical relay structure
Chapter II
2.1 Introduction
Modern digital and numerical relays are widely employed in protection systems
nowadays. Designing and modeling of numerical relay require establishing a generalized
numerical relay structure, which is composed by the more relevant and common internal modules
employed by typical numerical relays. The present chapter discusses the functionality of each of
the internal modules of the generalized numerical relay, namely signal conditioning and scaling
module, analog anti-aliasing filtering module, analog-to-digital conversion, phasor estimation
algorithm and relay logic. The most common techniques and methods employed in each of these
internal modules are enumerated and reviewed.
Reliability
The reliability of a relay is directly in correspondence with the concepts of dependability
and security. A relay is said to be dependable when it operates in the occurrence of a fault
relevant to its protection zone. Security is reached either when the relay will not operate for a
fault outside its operating zone, or when the system is in a healthy state.
Selectivity
Selectivity is the ability that a relay has to only open those breakers that isolate the faulted
element. Selectivity discrimination can be achieved by time grading or by unit protection.
Selectivity by time grading means that different zones of operation are graded by time and
that in the occurrence of a fault, although a number of protections equipment respond, only
those relevant to the faulty zone complete the tripping function. Selectivity by unit protection
means that the relay will only operate under certain fault conditions occurring within a clearly
defined zone.
Speed
In the occurrence of a fault, the greater the time in which the fault is affecting the power
system, the greater is the risk that the power system falls into an unstable operation point.
Relays are the greater is the risk that the power system falls into an unstable operation point.
Relays are therefore required to clear the fault as quickly as possible.
Page 9
Chapter II
Sensitivity
The relay is said to be sensitive if the relay operates to the minimum value of faulted
input signals.
Discrimination
This property allows the relay to distinguish between faults and some transient
phenomena like: an overload or transient over current in the case of transformers which is
caused by the inrush current. In addition, power swings in interconnected systems must not be
considered as faults.
Electromechanical relays
The first relays employed in the electric industry were electromechanical devices. These
relays worked based on creating a mechanical force to operate the relay contacts in response
to a fault situation. The mechanical force was established by the flow of a current that
reflected the fault current through windings mounted in magnetic cores. Due to the nature of
its principle of operation, electromechanical relays are relatively heavier and bulkier than
relays constructed with other technologies. Besides, the burden of these relays can be
extremely high, affecting protection purposes. However, electromechanical relays were so
largely employed, tested and known that even modern relays employ their principle of
operation, and still represent a good choice for certain conditions of application.
Solid-state relays
With the advances on electronics, the electromechanical technology presented in the
relays of the first generation started to be replaced by static relays in the early 60s. Static
relays defined the operating characteristic based in analog circuitry rather than in the action of
windings and coils. The advantages that static relays showed over electromechanical relays
were a reduced size, weight and electrical burden.
However, static relays showed some disadvantages since analog circuitry is extremely
affected by electromagnetic interference and the ranges of current and voltages values are
strongly restricted in analog circuits, affecting the sensitivity of the relay.
Page 10
Chapter II
Digital relays
Incorporating microprocessor into the architecture of relay to implement relay and logic
functions started happening in the 80s. Digital relays incorporated analog-to-digital converter
(ADC) to sample the analog signals incoming from instrument transformers, and used
microprocessor to define the logic of the relay. Digital relays presented an improvement in
accuracy and control over incoming signals, and the use of more complexes relay algorithms,
extra relay functions and complementary task.
Numerical relays
The difference between numerical relays and digital relays lies in the kind of
microprocessor used. Numerical relays use digital signal processors (DSP) cards, which
contain dedicated microprocessors especially designed to perform digital signal processing.
Chapter II
finally provides outputs for controlling the circuit breakers. The processor may also
support communications, self-testing, target display, time clocks, and other tasks.
In Figure 2.1 the schematic of a generalized numerical relay structure is shown. The
functionalities of each module of the generalized relay model are developed in next sections.
Page 12
Chapter II
Page 13
Chapter II
needed to perform the distance relay functions is retained and errors due to aliasing are avoided.
In practice, the sampling rate must be at least four samples per cycle (240 Hz).
Chapter II
achieve this conversion. The full-input voltage ranges for an ADC are typically 0 to +5 or 0 to
+10 volts for unipolar operations, and 5 to +5 or 10 to +10 volts for bipolar operation.
( )
(2.1)
being the frequency of the signal in radian per second, is the phase angle in radian and Xm
is the peak amplitude.
Page 15
Chapter II
( )
. The
(2.2)
(2.3)
To extract the fundamental frequency component from a signal corrupted with other frequency
components, set k=1 in equation to obtain the phasor estimate for M data samples.
(2.4)
Page 16
Chapter II
recursive algorithm. The phasor estimate for N samples from n=1 to n=N by non-recursive
algorithm is given by:
)]
(2.5)
A modified algorithm which saves computation taking into account data from previous window
is called as recursive algorithm which can be given by the equation below:
( )
(2.6)
(2.7)
Page 17
CHAPTER 3
Introduction
Proposed Relay model and tools required
Hardware Part
Software Part
Chapter III
Circuit breaker
Transmission line
CT
Load
PT
Signal Cond.Circuit
Acquisition card +
PC based program
Numerical distance relay
Figure 3.1 the general block diagram of the proposed protection scheme.
Page 18
Chapter III
Current Transformers
Smicro ModuleCS100-VP current
transformer Figure (3.2).
Where:
o Ip: primary current (input).
o Im: Secondary current (output).
Page 19
Chapter III
Features:
o Noncontact measure the high current.
o Measures DC, AC and impulse currents
o Current sensing up to 200A peak
o Very fast response and high accuracy
o High overload capacity
o Temperature range -25C to +85C
In this circuit, three current transformers (CS100-VP modules) are used.
Potential Transformers
As it is mentioned in chapter 2 the output of a voltage transformer in transmission line is
applied to an auxiliary transformer that reduces the voltage level and provides electrical
isolation to the rest of the relay equipment. We have used 220/12V transformer.as shown in
figure 3.4.
Page 20
Chapter III
| (
)|
(3.1)
| (
Then
)|
=>
( )
(
(
(3.2)
( )
( )
( )
(3.3)
(3.4)
axis.
Finding H(s) from H(s) H (-s):
H(s) is assigned all RHS poles and H (-s) is assigned all LHS poles
Following this procedure, the Butterworth LPF H(s) (H0=1, wc=1rad/sec) can be found for
various filters of order n.
( )
(3.5)
( )
(3.6)
( )
(3.7)
MATLAB is used to get this denominator polynomial (Butterworth polynomial) of the fifth
order filter.
( )
( )
(3.8)
)(
)(
(3.9)
Page 21
Chapter III
To design of a fifth order Butterworth low-pass filter with a cutoff frequency of 80Hz and
gain of K=1 a Sallen & Key Topology is used.
The general form of the transfer function of a Sallen& Key Topology is:
( )
(3.10)
If the gain k=1, the transfer function of the Sallen & Key will be:
( )
(3.11)
Computation
To determine the resistance and the capacitance of the first stage shown in figure 3.7 we have:
Chapter III
( )
(3.12)
To find the values of the resistors and the capacitors of the second and third stage:
From (3.9) and (3.12) we have:
( )
(3.13)
( )
(3.14)
C21=2*Q1 =>C21=1.236
C22=1/2Q1 =>C22=0.809
Page 23
Chapter III
C21
C22
C31
C32
Calculated value
2 F
2.46 f
1.6 F
6.5 F
618 F
New value
2 F
2.2F
1F || 470F
6.8F
470F
Page 24
Chapter III
There are five components to be considered when building a basic DAQ system
o Transducers and sensors
o Signals
o Signal conditioning
o DAQ hardware
o Driver and application software
In this chapter we focus on Signals. The appropriate transducers convert physical phenomena
into measurable signals. However, different signals need to be measured in different ways.
For this reason, it is important to understand the different types of signals and their
corresponding attributes. Signals can be categorized into two groups:
o Analog
o Digital
Analog Signals
Analog input is the process of measuring an analog signal and transferring the
measurement to a computer for analysis, display, or storage. Figure 3.9. An analog signal is a
signal that varies continuously. Analog input is most commonly used to measure voltage
or current. You can use many types of devices to perform analog input, such as multifunction
DAQ (MIO) devices, highspeed digitizers, digital multimeters, and Dynamic Signal
Acquisition (DSA) devices.
Page 25
Chapter III
Digital Signals
A digital signal cannot take on any value with respect to time. Instead, a digital
signal has two
possible levels:
high
and low.
certain specifications that define characteristics of the signal. Digital signals are commonly
referred
to as
transistortotransistor
logic
(TTL).
TTL
specifications
indicate
a digital signal to be low when the level falls within 0 to 0.8 V, and the signal is high between
2 to 5 V. The useful information that can be measured from a digital signal includes the state
and the rate.
Chapter III
As we know any Virtual Instrument (VI) has two Windows Block Diagram and Front Panel.
Block Diagram
In this window we simulate all the algorithms needed in our design, including:
o Phasor estimation algorithms: recursive algorithm is used, this stage analyze the
signals coming from DAQ Assistant and gives us the amplitude and phase of each
voltage and current in the three phase system. SeeFigure (3.10).
Figure 3.10 Block Diagram for Phasor Estimation Algorithm using Recursive DFT.
o Input signals (SubVI): this SubVI calculates all the impedances that are mentioned in
table 1.1 and gives six impedances.
o Point location (SubVI): this SubVI draws the mho characteristic with three zones and
gives graphically the location of the six impedances mentioned before.
o Comparator (SubVI): this SubVI compare the six impedances calculated from input
signals (SubVI) with the transmission line impedance for the three zones.
o Fault loc (SubVI): this SubVI calculate the fault location for each type of fault.
The remaining blocks includes the zones time delay and other mathematical calculations. As
shown in Figure (3.11).
Page 27
Chapter III
Chapter III
Front Panel:
This window is the HMI of the numerical distance relay. It is used to enter the
transmission line parameters, zones setting and time delay. It displays the amplitude and
phase of the three phase voltages and currents, fault type and location, and circuit breaker trip.
Besides there is a graph used to display the three zones of mho characteristic and the six
impedances as points. Figure (3.12).
Page 29
Chapter III
As a summary for the whole function of the developped Mho distance relay. The following
flow-chart explains the procedure step-by-step:
Start
Set relay setting (line parameters, zones length, time delay)
NO
yes
Zone 1
NO
Yes
Zone 2
NO
Delay
yes
Zone 3
Delay
Trip signal
Figure 3.13 Mho distance relay flow-chart.
Page 30
CHAPTER 4
Introduction
Power System Simulink Model
Testing procedures
Testing results
Testing results discussion
Implementation
Chapter IV
4.1 Introduction
This chapter will present the simulation of the mho distance relay developed in
LabVIEW, including testing procedures, results and discussion, and will end up with the
implementation of the proposed prototype.
Parameters
Value
Line length
100 km
Voltage
220 KV
Frequency
50 Hz
Line resistance R1
0.01165 Ohm/Km
Line resistance R0
0.2676 Ohm/Km
Line reactance X1
0.2725 Ohm/Km
Line reactance X0
0.9495 Ohm/Km
Relay settings
Zones
Settings
One
80%
Two
120%
Three
160%
Page 31
Chapter IV
Page 32
Chapter IV
Chapter IV
Page 34
Chapter IV
Page 35
Chapter IV
4.5 Implementation
The implementation of the proposed mho distance relay is shown in figure 4.9.
It consists of:
o Three potential transformer (ASTONIA220/12V).
o Three current transformers (Smicro Module CS100-VP)
o Low-Pass Filters.
o Acquisition CardNI USB-6009.
o PC ( LabVIEW model for Mho distance relay)
Page 36
Chapter IV
Chapter IV
References:
[1] M.P.Thakre, V.S.Kale, distance protection for long transmission line using pscad, International
Journal of Advances in Engineering & Technology, Jan. 2014.
[2] Dr. Hamid H. Sherwali and Eng. Abdlmnam A. Abdlrahem, Simulation of numerical distance
relays, Al-Fatah University Tripoli- Libya, 2010.
[3] Omar G. Mrehel Hassan B. Elfetori AbdAllah O. Hawal , Implementation and Evaluation a
SIMULINK Model of a Distance Relay in MATLAB/SIMULINK,2013
[4] Auday A.H. Mohamad, Essar Gafar Ahmed Design a Fast Digital Protective Relay Algorithm
For High Voltage Transmission line Received Feb. 2014, accepted after revision May 2014
[5] Stanley H . H or owitz and A r un G . Phadke, Power System Relaying, Third Edition.2008
[6] Sandro Gianny Aquiles Perez, Modeling Relays for Power System Protection Studies, Ph.D
Research, University of Saskatchewan, Saskatchewan, Canada, July 2006.
[7] H.Bentarzi, A. Ouadi, M. Chafai and A. Zitouni, Distance Protective System Performance
Enhancement Using Optimized Digital Filter, in Proc. CSECS '11 The 10th WSEAS International
Conference on CIRCUITS, SYSTEMS, ELECTRONICS, CONTROL & SIGNAL PROCESSING,
Montreux University, Switzerland, December, 29-31, 2011.
[8] Alstom, network protection & automation guide, edition may 2011
[9] Gerhard Ziegler, Numerical Distance Protection Principles and Application, Publicis
Corporate Publishing, Erlangen, Siemens, third edition, 2008.
[10] Sourav Mondal, Ch. Murthy, D. S. Roy, D. K. Mohanta, Simulation of Phasor Measurement Unit
(PMU)Using Labview
[11] http://www.ece.uic.edu/~jmorisak/blpf.html
[12] Vinicius J. & Osvaldo S. Using LabVIEW in a Mini Power System Model Allowing Remote
Access and New Implementation. International Conference on Engineering Education, 2007.
Conclusion
Conclusion
In this work, a new protection scheme that is based on the Mho distance relay algorithm
has been implemented using LabVIEW. After that the implemented Mho distance relay has been
tested. Then, the relay prototype has been realized using acquisition card NI USB-6009.
At beginning, we have presented operation principles of distance relay. We also discussed some
key aspects of distance relay protection, such as protection zones, comparators, the input signals
to the relay and Mho operating characteristic.
Consequently, the structure of a generalized numerical relay is introduced. The major internal
modules of the generalized relay model are described. These modules are the analog signal
scaling module, analog anti-aliasing filtering module, analog-to-digital conversion module,
phasor estimation algorithm module and relay logic module. The most common techniques and
methods employed in each module of the generalized numerical relay have enumerated and
reviewed.
We have then proceeded to design a new distance relay model. We have presented the proposed
relay model and tools required. Therefore, we have divided our design in two major parts,
hardware part and software part. We ended up with a flowchart that summaries the whole
function of the developped Mho distance relay.
Finally, a Mho type distance relay has been successfully developed based on LabVIEW
software. By testing the behavior of the developed relay model under different fault conditions,
the relay model has been able to recognize the appropriate fault type. From perspective
impedance calculations, the relay model has the ability of indicating the correct zone of
operation in all cases. The relay identifiers the fault locations as expected, as the fault location is
changed, the measured impedance changes consequently.
After the test, it can be noticed that the obtained results satisfy the principle operation of
numerical distance relay and its characteristics using this new frame work. Moreover, it can be
concluded that this proposed scheme has the following advantages:
1. The Mho distance protection can rapidly and reliably operate during power faults.
2. This complex protection scheme can easily be implemented on PC.
Page 39
Conclusion
3. This project is suitable for education for showing to the power engineering students the
distance relay principle of function and how to adjust it for protecting the different zones
of the transmission line.
For enhancing more the performance of the relay, power swing blocking protection function may
be implemented that can be considered as further work. Besides, quadrilateral distance relay
which is very suitable for avoiding some power swing situation, can be implemented as extended
work.
For more sophisticated distance relay that may be proposed is to design an adaptive relay which
can choose which type of distance relay may be used (Mho or quadrilateral) for such situation
and it will decide if it will use power swing blocking protection function or not. Because, the use
of the last protection function has certain limitations.
Page 40