Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2975 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay Kim Hung Le University of Science and Technology - The University of Da Nang, Da Nang, Vietnam lekimhung@dut.udn.vn Phan Huan Vu Center Electrical Testing Company Limited Da Nang, Vietnam vuphanhuan@gmail.com Abstract—This paper presents an application of a certain distance protection relay with a quadrilateral characteristic approach for the protection of the 110kV Duy Xuyen - Thang Binh transmission line in Vietnam using measured data from one terminal line. We propose the building process of a Matlab Simulink model for this relay that combines fault detection and classification block, apparent impedance calculation block for all types of faults and a trip logic block of three zone protection coordination. The proposed relay model is further tested using various fault scenarios on the transmission line. It is important to assess what happened, the actual conditions, the causes of mal- operation etc. Detailed explanation and results indicate that the proposed model behavior will help users to perform tests which correctly simulate real-world conditions besides that it properly interprets test results and troubleshoot distance function problems when results are not as expected. Keywords-transmission line; distance protection relay; quadrilateral characteristic; logic trip; Matlab/Simulink I. INTRODUCTION Transmission lines are one of the main components in high voltage and extra high voltage power system. The main protection used to protect Vietnam's transmission lines is a numerical distance protection relay (F21) as shown in Figure 1 because of its suitability, simplicity, economy, and reliability. There are now various distance relays that have been implemented by different vendors. These relays help to protect safeguard transmission line from exposure to any type of faults (internal or external) which could intentionally or inadvertently damage electrical equipment and avoid a great effect on the stability of the entire power system. However, the distance protection also faces a number of factors which influence negatively its operation, and thus affect its accurate determination of the measured impedance. This then results in inaccurate fault location determination [1]. A relay device manufacturer, or a power system operator, needs to conduct a simulated use of all test scenarios due to factors such as human factor, fault resistance and load current, coupling of parallel lines, multiple feeder sections, power swing or multi-end sources testing to evaluate the operating behavior of relay’s characteristics. Many operators are uncertain about the advantages and disadvantages of mho or quadrilateral characteristic related to fault type, system conditions, and look up the backgrounds of the mal-operation to remove errors. In practice, the relay provides satisfactory results for the operation condition of low resistance faults. But, in the case of high resistance faults, particularly for single phase to earth faults, it measures an incorrect value of fault impedance. Consequently, the relay may over-reach or under-reach depending upon the forward/backward flow direction and the magnitude of fault resistance [2]. In order to prove this issue, the current paper focuses on the algorithm of Siemens 7SA522 relay that uses a quadrilateral characteristic. Thereafter, the 110kV Duy Xuyen - Thang Binh transmission line model is selected as an example of fault simulation and the algorithm relay is tested in Simulink. The proposed model has been assessed under several scenarios. Conclusions based on the results are discussed along with the resistance fault factor affecting the relay’s accuracy. Fig. 1. Distance protection relay on transmission line II. RELAY PROTECTION ZONES In [3, 4], the disadvantage of the Mho relay is that, when it is used on long lines and the reach does not cover the section sufficiently along the resistance axis, it is incapable of detecting faults with high arc or fault resistances. The problem is aggravated in the case of short lines since the setting is low and the amount of the R axis covered by the mho circle is small in relation to the expected values of arc resistance. So that, the quadrilateral characteristic is a good choice to overcome the disadvantages of mho element when it is properly designed. In the case of medium lines, earth faults can, for instance, be covered with the quadrilateral characteristic and phase faults with the mho characteristic. In Vietnam, the 110kV Duy Xuyen–Thang Binh transmission line has length less than 30 km and, as shown in Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2976 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay Figure 2, is protected by a Siemens 7SA522 distance protection relay with three forward quadrilateral zones (Z1, Z2, Z3) which provide fault impedance coverage for both phase to phase and phase to earth faults. Setting calculation is important in protection zone setup, so human error can affect the accuracy of the relay protection besides the large amounts of faults to be collected and processed. The transmission line parameter (Line length=29 km, Line angle=68o, line reactance=0.3304 Ω/km), Zero sequence compensation factor (RE/RL=0.29, XE/XL=0.67), and fault resistance must also be taken into account. The zone settings of this transmission line are calculated by the central region load dispatch center in Vietnam ensuring that they are properly coordinated with the adjacent system and comply with the existing guidelines. The three zones of quadrilateral characteristic used to protect transmission line are explained in Table I [5]. We used the DIGSI operating software to enter configuration settings in the relay, it can automatically draw a relay operating characteristic graph as shown in Figure 3. The quadrilateral characteristic permits separate setting of the reactance X and the resistance R. The resistance section R can be set separately for faults with and without earth involvement, therefore, this characteristic has an optimal performance in case of faults with fault resistance [6]. The effect of fault resistance on the reach of distance relays is better discussed with the simulation results. Fig. 2. The 110kV power system TABLE I. PARAMETER SETTING ON SIEMENS 7SA522 ON A 110KV DUY XUYEN–THANG BINH OVERHEAD LINE Zone Relay Setting Setting Value Zone 1 Covers 85%AB, forward direction (0.85ZAB) Operating Mode Z1: Forward X(Z1), Rectance: 8.14Ω R(Z1), Resistance for Ph-Ph faults: 30Ω RE(Z1), Resistance for Ph-E faults: 40Ω tZ1, Delay for fault: 0s Zone 2 Covers 100%AB + 60%BC, forward direction (ZAB + 0.6ZBC) Operating Mode Z2: Forward X(Z2), Rectance: 11.5Ω R(Z2), Resistance for Ph-Ph faults: 35Ω RE(Z2), Resistance for Ph-E faults: 50Ω tZ2, Delay for fault: 0.3s Zone 3 Covers 100%AB + 100%BC + 60%CD, forward direction (ZAB + ZBC + 0.6ZCD) Operating Mode Z3: Forward X(Z3), Rectance: 19.53Ω R(Z3), Resistance for Ph-Ph faults: 40Ω RE(Z3), Resistance for Ph-E faults: 55Ω tZ3, Delay for fault: 1.5s III. POWER SYSTEM UNDER STUDY We model the power system that is supplied from both ends as shown in Figure 4. It uses Simulink to evaluate the distance protection function. The overhead line is 89.65 km long, and the system nominal operating voltage is 110 kV, 50Hz. Fig. 3. Quadrilateral characteristic of Siemens 7SA522 at the 171 overhead line of Duy Xuyen Substation Fig. 4. Power system model This model consists of: 1. The transmission line: three phase section line is used to represent the transmission line. Line sequence impedance is: RL1=0.1335Ω, RL0=0.2496Ω. LL1=1.0517mH, LL0=3.1657mH. CL1=0.038μF, CL0=0.038μF. Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2977 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay 2. A load of 110 kV, 35 MW, and 50 kVAR is connected to the bus Duy Xuyen, Thang Binh, Thang Binh 2, Tam Thang. 3. Three phase fault block to deduce fault types and fault resistance varying from 1 to 40 ohms. 4. Three-phase measuring blocks to measure the three phase line, load current and voltage values. 5. The Scopes already support the display of five inputs (Voltages, Currents, Zones 1-3) of a simulated system, which allows a better signal comparison. 6. A Siemens 7SA522 relay is located at bus Duy Xuyen. It has been developed for fault detection, fault location, and distance protection which will be presented in section IV. IV. DISTANCE PROTECTION RELAY SIEMENS 7SA522 MODEL The simulation model of the relay has been designed with four functional flow block diagrams (Figure 5) by using Matlab’s Simulink library. Firstly, measurement block reads the values of voltage and current of each phase. Next, fault detection and classification are done by FD block and Fault_Type block. Thirdly, impedance calculation block computes the value of RCAL and XCAL. Afterwards, Trip_logic block checks operating logics. If all four conditions are satisfied then they issue a trip signal otherwise all steps are repeated from start to end. This helps us monitor and analyze data to avoid severe misunderstanding in assessing any use issues that arose during testing. Fig. 5. Distance protection model A. Current and Voltage Measurement The measurement block receives two input signals corresponding to the three phase voltages and currents from CT, VT and then it provides to the multiplexer to split the voltage and current signals into three single phase voltages and currents. After that, these quantities are passed through the Fourier transform block to avoid the higher order harmonics as shown in Figure 6. B. Fault Detection and Fault Classification According to the outputs of the measurement block (IA, IB, IC, and IN), the fault detection in the transmission line is defined by the FD block. For earth current measurement, the fundamental component of the sum of the numerically filtered phase currents is supervised to detect if it exceeds the set value (parameter 3I0>0.5IN) then output “N” equals 1. For phase current measurement, a fault occurrence is detected if IA, IB, or IC is greater than the threshold 1.3INorm. Then output “A”, or “B”, or “C” is equal to 1 [7]. Based on the information above, fault_type block performs fault classifying on the transmission line which can be summarized as in Table II. Fig. 6. Measurement block TABLE II. FAULT TYPE CLASSIFICATION No Type fault A B C N 1 AG 1 0 0 1 2 BG 0 1 0 1 3 CG 0 0 1 1 4 AB 1 1 0 0 5 BC 0 1 1 0 6 CA 1 0 1 0 7 ABG 1 1 0 1 8 BCG 0 1 1 1 9 CAG 1 0 1 1 10 ABC 1 1 1 C. Impedance Calculation To calculate the RCAL, XCAL of relay protection, it uses one of the six loop measurements (AG, BG, CG, AB, BC, and CA) listed in Table III. The measured fault loop is necessary to know the fault type, local voltages and currents along with the zero sequence compensation factors. The following equations are employed to evaluate the measured values [7]. For the calculation of the phase-to-earth loop, for example during an AG short-circuit (Figure 7) it must be noted that the impedance of the earth return path does not correspond to the impedance of the phase. In the faulted loop the voltage AG, the phase current IA and the earth current ΙE are measured. The impedance to the fault location results from: 2 cos( ) cos( ) 1 ( ) cos( ) E E U A U E A A L A A E E E E E E E A L L A L L A I X U I X R I X R I R X I X R I R X I                      Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2978 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay 2 sin( ) sin( ) 1 ( ) cos( ) E E U A U E A A L A A E E E E E E E A L L A L L A I R U I R X I X R I R X I X R I R X I                       TABLE III. SHORT CIRCUIT TYPES AND FAULT LOOPS FOR THE DISTANCE MEASUREMENT No Type fault Phases involved Fault loop for distance measurement 1 Single phase to ground fault AG BG CG AG BG CG 2 Two phase short circuit without ground AB BC CA AB BC CA 3 Two phase short circuit with ground ABG BCG CAG AG, or BG, or AB BG, or CG, or BC AG, or CG, or CA 4 Three phase short circuit ABC AB, or BC, or CA Fig. 7. Single-phase earth fault, fault loop In (1) and (2), UA is the short-circuit voltage of phase A, ΙA is the short-circuit current of phase A, ΙE is the earth short- circuit current, φU is phase angle of the short-circuit voltage, φA is phase angle of the phase short-circuit current, and φE is phase angle of the earth short-circuit current. To calculate the phase- to-phase loop for instance during a two-phase short circuit BC (Figure 8), the loop equation is: 2 2 cos( ) cos( ) 2 cos( ) B B C C B C B B U I C C U I BC B B C I I C U I U I R I I I I                 2 2 cos( ) cos( ) 2 cos( ) B C C B B C B C U I C B U I B B C I I C U I U I I I I I               2 2 sin( ) sin( ) 2 cos( ) B B C C B C B B U I C C U I BC B B C I I C U I U I X I I I I                 2 2 sin( ) sin( ) 2 cos( ) B C C B B C B C U I C B U I B B C I I C U I U I I I I I               Fig. 8. Two-phase fault clear of earth, fault loop In (3) and (4), IB, and IC represents short circuit phase current of phase B, and C. UB, and UC represents short circuit voltage of phase B, and C. φUB, φUC are considered as phase angle of UB, UC. φIB, φIC are considered as phase angle of IB, IC. D. Logic Decision As mentioned in the description of the calculation impedance, the calculated RCAL and XCAL coordinate values define a point on the complex impedance plane. The logic decision compares the point with the quadrilateral characteristics of the relay, shown in Figure 3. For example, the distance characteristic has been programmed as seen in Figure 9 [7]. The point R3 and X3 are the setting points for Zone 3. The time settings tZ3=1.5s. If a measured impedance point is inside the Zone 3, then the delay time can be started and when the impedance stays longer than tZ3, the logic decision generates the true value of the related output binary signal. It means that relay will trip with Zone 3 if the following conditions are met:  Condition 1:   3 30 0 03 3 3 30 0 0 0 3 22 68 – 68 22 68 6 < 8 6 < 8 CAL CAL CAL CAL R X cotan cotan R R R X tan cotan cotan tan R R X tan X R                  Condition 2:   3 30 0 3 3 0 0 0 0 22 68 0 – 22 68 68 22 CAL CAL CAL CAL R X X cotan cotan R R R cotan cotan tan R X tan                    Condition 3: 3 0 0 30 0 CAL CAL CAL X X X tan R        Fig. 9. Zone 3 protection V. SIMULATION RESULTS After the building of the proposed model has been completed, it is ready to analyze the secure operation of the quadrilateral characteristics of the three protection zones applied with variation in fault parameters like fault types, fault location from -5 km to 50 km, fault resistance from 1 Ω to 40 Ω. During a fault, the current in the faulted phases increases. The voltage in the faulted phases decreases beginning at 0.3 s. Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2979 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay Distance protection zone outputs such as Z1, Z2, Z3, operate when the corresponding pre-set times are reached. After 10 ms the line is tripped and the voltage begins to recover. In case of an AG fault with RF=30 Ω is created on the line Duy Xuyen- Thang Binh at a distance of 35 km from Duy Xuyen bus (forward direction). The distance represents 120.7% of the line length. Figure 10 shows the current and voltage waveform and a trip signal of protection zones before, during and after a fault. The fault occurs at 0.3 s, the current of phase A is increased and voltage of phase A is decreased. As a result, impedance trajectory seen by the relay entered into the characteristic second zone. Fault location displays 28.38 km. Relay takes 0.31 s to send trip signal Z2 to the circuit breaker. Figure 11 shows the operation of relay in zone 1 with the presence of AB fault at distance 15 km with RF=10 Ω. The distance represents 51.7% of the line length. According to the change in voltage and current waveforms during the fault, the impedance trajectory moves quickly into Zone 1, fault location displays 14.29 km which results in a tripping at 0.322 s. In Figure 12, the result is shown where BCG fault is applied on the transmission line at distance 45 km with RF=40 Ω. According to the change in voltage and current waveforms during the fault, the apparent impedance seen at relay point is reduced but the trip signals are deactivated due to the movement of impedance trajectory out of Zone 3 area. Fault location displays 37.65 km. As shown in Figure 13, at 0.3 s an ABC fault occurs on the transmission line at 50 km, RF=1 Ω, and the fault duration is 1.5 s. According to the change in voltage and current waveforms during fault, relay impedance trajectory into Zone 3, fault location displays 49.47 km, thus trip signal Z3 is activated at 1.81 s. Some of the results for various fault cases (AG, BC, ABG, and ABC) are given in Table IV. Fig. 10. Voltage, current waves, trip sinals and impedance trajectory for faults at AG fault with RF=30Ω at 35km from relay point Fig. 11. Voltage, current waves, trip sinals and and impedance trajectory for AB fault with RF=10Ω at 15km from relay point Fig. 12. Voltage, current waves, trip sinals and impedance trajectory for faults at BCG fault with RF=40Ω at 45km from relay point. The power system is supplied from both ends, its setting for Zone 1 is 85 percent of line length or 24.65 km. Similarly, Zone 2 is 35.18 km, and Zone 3 is 57.9 km. In the tests, the fault resistance influence on the measured distance to fault location, the highest fault resistance that allowed for reliable distance protection relay is 20 Ω. In the case of the high Engineering, Technology & Applied Science Research Vol. 8, No. 3, 2018, 2975-2980 2980 www.etasr.com Le and Vu: Effect Evaluation of Fault Resistance on the Operating Behavior of a Distance Relay resistance fault, this value changes from 20 to 25 Ω, the final point of fault impedance locus falls under Zone 1 characteristic or under-reaching tendency. If it is higher than 40 Ω, fault impedance locus out of zone protection with corresponding error increases. TABLE IV. TEST RESULTS Phase Actual fault location (km) Fault resistance (Ω) Measured fault location (km) Zone trip AG -5 10 -12.2 / 5 1 4.718 Z1 10 5 9.244 Z1 15 10 13.7 Z1 20 15 18.21 Z1 25 20 22.78 Z1 30 25 24.75 Z1 35 30 31.42 Z2 40 35 35.25 Z2 45 40 38.35 / 50 1 49.67 Z3 AB -5 10 -9.475 / 5 1 4.923 Z1 10 5 9.633 Z1 15 10 14.29 Z1 20 15 18.96 Z1 25 20 23.63 Z1 30 25 28.15 Z1 35 30 31.92 Z2 40 35 35.31 Z2 45 40 37.65 / 50 1 49.47 Z3 BCG -5 10 -9.455 / 5 1 4.919 Z1 10 5 9.919 Z1 15 10 14.92 Z1 20 15 19.93 Z1 25 20 23.62 Z1 30 25 28.14 Z1 35 30 31.91 Z2 40 35 35.31 Z2 45 40 37.65 / 50 1 49.47 Z3 ABC -5 10 -9.48 / 5 1 4.919 Z1 10 5 9.633 Z1 15 10 14.29 Z1 20 15 18.95 Z1 25 20 23.62 Z1 30 25 28.15 Z1 35 30 31.92 Z2 40 35 35.31 Z2 45 40 37.65 / 50 1 49.47 Z3 VI. CONCLUSIONS In this study, the distance protection function with three forward zones has been implemented in a Siemens 7SA522 distance protection relay on the 110kV Duy Xuyen–Thang Binh transmission line that was simulated easily and reliably in Matlab Simulink. The proposed model allows us to include the detailed analysis of the impedance (RCAL, XCAL) measured by the distance relay on one single end, fault detection, fault classification, and tripping protection zones for fault, considering fault location, fault type existing on the transmission line, and fault resistance. According to the obtained results, it has been shown that the accuracy of the relay depends on a fault resistance value. The high fault resistance causes the distance relay to over-reach or under- reach. Therefore, a correctly set distance protection relay can correctly compute the fault location, mitigate the problem of mal-operation, thus ensuring effective and reliable protection to the transmission line. Fig. 13. Voltage, current waves, trip sinals and impedance trajectory for faults at ABC fault with RF = 1Ω at 50km from relay point ACKNOWLEDGMENT Authors would like to thank the "Electrical Testing Company Limited, Viet Nam" for allowing the use of the project of relay protection used in this study. REFERENCES [1] O. Muhayimana, Distance Protection for Parallel and Double-Circuit HV Lines, MSc Thesis, Brno University Of Technology, 2016 [2] V. H. Makwana, B. R. Bhalja, Transmission Line Protection using digital technology, Springer, 2016 [3] T. Saengsuwan, “Modelling of distance relays in EMTP”, IPST'99 International Conference on Power Systems Transients, Budapest Hungary, pp. 213-217, June 20-24, 1999 [4] L. K. Hung, N. H. Viet, V. P. Huan, “The effect of fault resistance on the performance of distance relay protection”, The 2011 International Symposium on Electrical-Electronics Engineering, Ho Chi Minh City, Vietnam, pp. 238-244, November 8-9, 2011 [5] Central Region Load Dispatch Center, Setting calculation of distance protection relay Siemens 7SA522 at 171 overhead line of Duy Xuyen Substation, Decision No A3-03-2015/DUX110, Siemens, 2015 [6] G. Ziegler, Numerical Distance Protection: Principles and Applications, John Wiley & Sons, 2011