Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3243 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh-Thanh My Transmission Line Kim Hung Le The University of Danang, University of Science and Technology, Da Nang, Viet Nam lekimhung@dut.udn.vn Phan Huan Vu Central Power Corporation, Center Electrical Testing Company Limited, Da Nang, Viet Nam vuphanhuan@gmail.com Abstract—This paper presents the traveling wave based fault location methods of SEL-400L, and SFL-2000 available on the market for a 66.9km, 220kV Hoa Khanh-Thanh My transmission line in Central Viet Nam, such as single-ended, and double- ended, all of which rely on measurements from inductive CTs and capacitive VTs. Focus was given on the building process of a Matlab Simulink model to evaluate these methods. Current and voltage signals were sent to an analog Chebyshev type II filter which passes higher frequency signals at 3kHz and rejects low frequencies signal at 50Hz. After that, these output signals are used in Clarke's transformation for getting 0 and α components. The detail coefficient of the selected components after DWT using Db4 wavelet at decomposition level 1 can be used to determine the fault types, the direction of fault and propose a crest-wave comparison solution to identify exactly the adjacent bus' reflected wave from the fault point's reflected wave for the fault location. Finally, the accuracy of fault location on the transmission line is reviewed by varying various parameters like fault type, fault location and fault resistance on a given power system model. Keywords-transmission line; traveling wave fault locator; single ended method; double ended method; Matlab/Simulink I. INTRODUCTION Numerical relays are the most popular devices used for transmission line protection. They include fault location estimation based on impedance methods, which uses the voltage and current data measured in 50Hz at one or more points along the power networks after the occurrence of a fault. Typical impedance method’s error ranges from 2% to 5% depending on the relay model according to industry standards. However, the actual error is usually larger than 5% in practice operation manager at EVN of Viet Nam [1]. It can be influenced by weather, high resistance ground faults, measurement errors, line impedance errors, mutual coupling, compensated lines, and other factors [2]. Therefore, finding the accurate location of a fault constitutes a challenge in the power operator. Selecting a new appropriate technology for a fault location application can be a daunting task so a performance evaluation is required when the EVN wants to compare accuracy and finds the optimum and the most cost-effective one. In 2017, EVN’s project installed a traveling wave fault locator (TWFL) equipment Kinkei SFL-2000, SEL-400L for 220kV, 500kV transmission lines in substations such as Son Ha, Thanh My, Hoa Khanh, Hue, Dong Ha, Tam Ky, and Doc Soi in central Viet Nam. The project implementation plan followed 3 steps. The first step was setting up the TWFL to work with single-ended method. The second step was synchronizing the TWFLs by the Furuno GPS/GNSS clock receiver and using a configured communication channel. The TWFLs would send the traveling wave arrival information to F/L server that calculates the distance to fault base on double- ended method, display results and sends email to the operator. The third step, EVN installs an NPT communication network link between TWFLs as shown in Figure 1. The TWFLs receive the remote traveling wave information that is necessary to provide an automatic fault location. SFL-2000’s accuracy is reported to be200m [3], whereas the accuracy for SEL-400L is reported to be 2% [2]. Although this issue was recognized by PTC2 as a good way to overcome the shortfalls of impedance-based methods, the SFL-2000 is not able to be used in step 1. So, if it is impossible to have communication between ends, the fault location cannot be implemented automatically by the F/L server. Note that, the problem faced by the power system operator is manually collected traveling waves event reports from each substation. They lack skills to discriminate the reflection wave of fault point, the time of arrival of the wave for manual fault location estimation using the single-ended method because the reflected surge from the adjacent bus has the same polarity as the real fault point's reflected wave, so the confusion is inevitable. Almost as many respondents said all the implementation and support work would be done by experienced professionals. Consequently, EVN needs to conduct a training session about the knowledge of the problem of fault location using traveling wave signals for communication technicians, field personnel, relay engineers, and analysis an in-depth review of all TWFL system operations. In order to solve this issue, the paper focuses on TWFL based double-ended and single-ended methods implemented in a 66.9km, 220kV Hoa Khanh–Thanh My transmission line by Matlab Simulink. The proposed model has been assessed through several scenarios. Results show that methods consistently and significantly yielded the accurate Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3244 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … location of the actual fault. II. TRAVELING WAVE BASED FAULT LOCATION METHODS The TWFL can work based on current or voltage signals. It has two kinds of the traveling wave (TW) methods that are used for most fault locator systems. One is the single-ended, which captures data from the initial traveling wave and subsequent reflections at one terminal to calculate the fault location without requiring any information from the relay at the remote terminal. The other is double-ended, which requires data information from two terminals which are both equipped with a GPS receiver to time tag the exact moment the traveling wave reaches each end of the line. Fig. 1. Traveling wave fault locator SFL-2000 on transmission line A. Single Ended Method Figure 2 for a fault at location F on a line of length L=66.9km (the Bewley diagram is shown in Figure 14). The fault is m (km) away from the Hoa Khanh terminal, that is suspected to be in the second half of the line Hoa Khanh-Thanh My. A current TW is launched from the F at t0=40ms and arrives at the Hoa Khanh terminal at t1B. To discriminate the reflection wave of fault point, the number of samples in the selected window is limited to the interval from t1B to 40.825ms after the fault occurs. Part of the wave transmission travels toward the bus Da Nang and then returns to the Hoa Khanh terminal at t’1A. It has opposite polarity with the first traveling wave (t’1A-t1B=constant). So we identify and eliminate false peak due to the effect of this wave. Another part of the wave reflects, travels back toward the fault, reflects back from the fault, and then returns to the Hoa Khanh terminal at t2B. It has larger crest-current wave than the reflected surge from the Thanh My terminal at t’1C with the same polarity. Similarly, in the previously examined case, the fault is assumed to occur at the first half of the line as shown in Figure 13 then the number of samples in the chosen window has limited the interval from t1B to 40.4125ms. Now, the distance to the fault location from TWFL B is [2]: 2 1( ) 2 B Bv t tm    (1) Wave propagation velocity is 1 11 243250 km/sv L C   (where L1=0.0013H/km is the inductance and C1=0.013μF/km is the capacitance of the propagation medium). When an external fault behind the Hoa Khanh terminal launches a traveling wave as shown in Figure 3, the TWFL B sees an initial wave behind it with t1B=t1, which travels across the transmission line to the Thanh My terminal and is reflected back to the Hoa Khanh terminal after the known TW line propagation time (tL=0.275ms) with t2B=2×tL+t1. The TWFL B displays L. Fig. 2. Single end TWFL with internal fault Fig. 3. Single end TWFL with external fault Review: The single-ended method estimates the accurate fault location of an internal fault and does not require data synchronization. However, the accuracy of this method depends on the accuracy of L, v, sampling frequency, and errors in wave detection. If the impulse wave cannot be captured successfully or the impulse wave does not exist at all at the fault occurrence, the fault location will fail. For instance, strong buses on power system network influencing the voltage and current waveforms due to line impedances can reduce the amplitude of voltage waves making them harder to detect, and thus reducing the TWFL accuracy [4]. B. Double Ended Method When an internal fault F occurs at t0, waves generated at F run towards stations Hoa Khanh and Thanh My (Figure 4). The double-ended method determines the arrival time based on the rising point of surge waveform recorded at both end terminals (t1B at the Hoa Khanh, t1C at the Thanh My) and then locates the fault point from the equation shown below [2, 3]. ( )1 1 2 2 2 v t tL L v tB Cm         (2) We can calculate v using the line length L=66.9km and the arrival time of the surge generated by manually closing a circuit breaker at Thanh My terminal as shown in Figure 5 [3]: v=L/(t1B–t1C)=243270km/s. When an external fault behind the Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3245 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … Hoa Khanh terminal launches a traveling wave, the TWFL B sees an initial wave with either positive or negative polarity with t1B=t1 which travels across the transmission line to the Thanh My terminal. The TWFL C sees the same initial wave tL=0.275ms later with opposite polarity (t1C=tL+t1). The distance to the fault location from TWFL B is L. Fig. 4. Double ended TWFL with internal fault Fig. 5. Calculate propagation velocity Review: Double-ended method requires data from both terminal ends to be synchronized. It estimates the accurate fault location of an internal fault. The accuracy of this method is effected by communication and precise timing stamp GPS. This method is more expensive than the single-ended method. III. POWER SYSTEM UNDER STUDY Recorded data from a real system are not available to evaluate the performance of the TWFL. Instead, the power system supplied from both ends can be modeled by Matlab Simulink software (Figure 6). The overhead line Hoa Khanh– Thanh My is 66.9km long, and the system nominal operating voltage is 220kV, 50Hz. This model consists of: 1. The transmission line: three phase section lines are used to represent the distributed parameter transmission line. Line sequence impedance is: RL1 =0.07(Ω/km), RL0 =0.2164(Ω/km). LL1 =0.0013(H/km), LL0 =0.0044(H/km). CL1 =0.013(μF/km), CL0 =0.0085(μF/km). 2. A load of 220kV, 56MW, and 34kVAR is connected to the bus Hoa Khanh and Thanh My. 3. Three phase fault block to deduce fault types and fault resistance varies from 1 to 35 ohms. 4. Three-phase measuring blocks to measure the three phase line and load current and voltage values. 5. A TWFL model is located at Hoa Khanh bus. It has been developed for fault detection, fault classification, fault direction and fault location which will be presented in section IV. IV. TRAVELING WAVE FAULT LOCATOR MODEL This subsection helps to understand how the TWFL works. The simulation model of the TWFL has been designed with six functional flow block diagram (Figure 7) by using Simulink. Fig. 6. Power system model Fig. 7. Flowchart of travelling wave fault locator operation A. Analog Filter A fault generates current and voltage traveling waves that propagate along the overhead line. Most of them contain a significant amount of high-frequency components. The TWFL collects TWs from conventional CT class 5P20, CVT class 3P. Then it uses an analog Chebyshev type II filter to remove the fundamental component 50Hz and a high-pass filter with a cutoff frequency of 3kHz for phase currents and voltages. B. Clarke’s Tranformation To reduce the effect of mutual coupling between phases, this paper utilizes the Clarke transformation to convert the three-phase currents and voltages into the α, ß and 0 mode components. For example, we use three sets of Clarke components with reference to the A-phase, B-phase, and C- phase of current signals, as follows [5]: 0 0 0 2 1 1 1 0 3 3 3 1 1 1 1 2 1 1 3 0 3 3 1 1 1 1 1 2 1 3 3 0 3 1 1 1 A A A B A C B A B B B C C A C B C C I I I I I I I I I I I I I I I I I I                                                                                        (3) C. Discrete Wavelet Transform This function is developed by using an α mode and 0 (ground) mode component which capture voltages and current with a sampling rate of 10MHz at TWFL A (UαA, U α B, U α C, I α A, IαB, I α C, IG), and TWFL B (Bus_B_I α A, Bus_B_I α B, Bus_B_I α C, Bus_B_IG). By using DWT Daubechies 4 mother wavelet (db4) Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3246 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … during the fault surges from the state steady sampled currents and voltages, we can rapidly extract the first level detail Wavelet Coefficients (Cd1), frequencies up to 5MHz which is enough for the TWFL transient frequency. Now it is easy to determine the times of traveling wave occurrences (peaks can be observed on waveform) and to reveal their travel times between the point fault F and TWFL. D. Fault Detection and Fault Classification Under normal conditions, Cd1 of phase A, B, C, and Ground are zero (Figure 8). Under fault conditions, if the Cd1 of ground is zero, the fault is identified as an ungrounded fault or as a grounded fault if it is nonzero. The Cd1 of phase A, B, and C are available for all fault types. Consider the AG fault occurring on the transmission line at time 40ms (Figure 9). Both reflection and refraction of the Cd1 of current phase A and ground with the large amplitude occur at 40.0 to 40.8ms (or 2.0×105 to 2.04×105 samples). They are approximately twice as large as the Cd1 of current phase B, C (healthy phases), and they have opposite polarity. Consider the BC fault (Figure 10), the Cd1 of current phase B and C are also considerably larger than the Cd1 in the healthy phase (A, ground), and these faulted phases have opposite polarities. Based on information relationship between the squared Cd1 magnitude of the first current traveling waves in each phase, the TWFL can make fault type decisions that can be summarized in Table I. Fig. 8. The Cd1 of phase A, B, C and ground at normal condition Fig. 9. The Cd1 of phase A, B, C and ground at AG fault E. Direction Fault The single-ended method can make a directional decision based on the polarity relationship between the first voltage and current traveling waves. For a fault on the transmission line in the forward direction, the voltage and current traveling waves observed by the relay have opposite polarity (Figure 11). For a fault on the transmission line in the reverse direction, the voltage and current traveling waves observed by the relay have the same polarity as in Figure 12. The double end method compares time-aligned current first TWs at both ends of the protected line. For an external fault, a TW that entered one terminal with a given polarity leaves the other terminal with the opposite polarity exactly after Δt=tL. For an internal fault, a TW that entered one terminal with a given polarity leaves the other terminal with the same polarity, Δt1.5 and Mag_A/Mag_C>1.5 and Mag_G/Mag_A>0.2 Cd1_A BG Mag_B/Mag_A>1.5 and Mag_B/Mag_C>1.5 and Mag_G/Mag_B>0.2 Cd1_B CG Mag_C/Mag _A >1.5 and Mag_C/Mag_B>1.5 and Mag_G/Mag_C>0.2 Cd1_C AB Mag_A/Mag_C>1.5 and Mag_B/Mag_C>1.5 and Mag_G/Mag_A<0.2 and Mag_G/Mag_B<0.2 Cd1_A or Cd1_B BC Mag_B/Mag_A>1.5 and Mag_C/Mag_A>1.5 and Mag_G/Mag_B<0.2 and Mag_G/Mag_C<0.2 Cd1_B or Cd1_C CA Mag_C/Mag_B>1.5 and Mag_A/Mag_B>1.5 and Mag_G/Mag_C<0.2 and Mag_G/Mag_A<0.2 Cd1_C or Cd1_A ABG Mag_A/Mag_C>1.5 and Mag_B/Mag_C>1.5 and Mag_G/Mag_A>0.2 and Mag_G/Mag_B>0.2 Cd1_A or Cd1_B BCG Mag_B/Mag_A>1.5 and Mag_C/Mag_A>1.5 and Mag_G/Mag_B>0.2 and Mag_G/Mag_C>0.2 Cd1_B or Cd1_C CAG Mag_C/Mag_B>1.5 and Mag_A/Mag_B>1.5 and Mag_G/Mag_C>0.2 and Mag_G/Mag_A>0.2 Cd1_C or Cd1_A ABC Mag_A>1.5 and Mag_B>1.5 and Mag_C>1.5 Cd1_A or Cd1_B or Cd1_C F. Distance Calculations Because the frequency response for CTs is better than for CCVTs [4], we chose to use current signals for fault location. In the case when TWFL identifies a forward or internal fault, both single end and double end methods are activated. We use Matlab's find peaks function to find values and locations of Cd1 maxima in a time period set according to the travel time of the line. The fault distance is given by (1) and (2). Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3247 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … Fig. 11. Curent and voltage polarities of phase A in fault forward direction Fig. 12. Curent and voltage polarities of phase A in fault reverse direction V. SIMULATION RESULTS After building, the proposed model is ready to analyze the operation of TWFL applied with variation in fault parameters such as fault type, fault location (from -5km to 65km), fault resistance from 1Ω to 35Ω (increasing the fault resistance seen at relay point as though moving away from the transmission line). The fault creation time is t0=40ms. Figure 13 shows the phase currents captured at both terminals and a Beley diagram for a BG fault with RF=10Ω that is assumed to occur at a distance of 20km from Hoa Khanh bus (forward direction). TWFL B calculates the fault location by the single-ended method with t1B=40.0836ms, t2B=40.2476ms, due to the traveling time taken by the fault to appear at the TWFL B. Based on the measured TW arrival times, an estimated from (1) fault location 19.9523km from the Hoa Khanh terminal emerged. In Figure 14, results are shown where a BCG fault with RF=25Ω at a distance of 45km from Hoa Khanh bus (forward direction). According to the traveling time taken by the fault to appear at the point t1B=40.1866ms, t2B=40.5566ms. The TWFL B also uses the single-ended method to calculate fault location is 45.0016km. The double-end method’s (Figure 15) performance is simpler than the single-ended one. In the case an AB fault occurs on the transmission line with RF=15Ω at 10km, it calculates the distance to fault based on time tags of traveling wave records acquired at both ends of the faulty line. The first peak of Hoa Khanh terminal occurs at t1B=40.0426ms, the first peak of Thanh My terminal occurs at t1C=40.2356ms, so the distance to the fault point F of TWFL B estimate from (2) is 9.9762km. As shown in Figure 16, an ABC fault occurs on the transmission line at 35km, RF=15Ω. According to the change in current waveforms during the fault, TWFL B captured point t1B=40.1456s, TWFL C captured t1C=40.1326s, thus the double- end method of TWFL B estimate is 35.0313km. Some of the results for various fault cases are given in Table II respectively for BG, AB, BCG, and ABC. Fig. 13. Phase current waves and Bewley diagram explaining single ended method for a BG fault at 20km from the Hoa Khanh terminal Fig. 14. Phase current waves and Bewley diagram explaining single ended mehtod for a BCG fault at 45km from the Hoa Khanh terminal Fig. 15. Phase current waves and Bewley diagram explaining double ended method for AB fault at 10km from the Hoa Khanh terminal Review: The simulation results show that single ended method locates faults with accuracy less than ±114m and double ended method locates faults with accuracy less than ±112m. Engineering, Technology & Applied Science Research Vol. 8, No. 4, 2018, 3243-3248 3248 www.etasr.com Le & Vu: Performance Evaluation of Traveling Wave Fault Locator for a 220kV Hoa Khanh … Fig. 16. Phase current waves and Bewley diagram explaining double ended method for ABC fault at 35km from the Hoa Khanh terminal TABLE II. TEST RESULTS OF FAULTS ON TRANSMISSION LINE Phase Actual fault location [km] Fault resistance [Ω] Single end Method Estimated Distance Double end Method Estimated Distance m [km] Difference [m] m [km] Difference [m] BG -5 1 / / 10 5 9.9733 -26.7 9.9762 -23.8 20 10 19.9523 -53.3 19.9495 -50.5 30 15 30.0416 41.6 29.9228 -77.2 40 20 40.015 15 39.967 -33 50 25 49.9883 -11.7 50.112 112 60 30 60.889 88.9 60.0859 85.9 AB -5 1 / / 10 5 9.9733 -26.7 9.9744 -25.6 20 10 19.9467 -53.3 19.9495 -50.5 30 15 30.0416 41.6 29.9228 -77.2 40 20 40.015 15 39.967 -33 50 25 49.9883 -11.7 50.112 112 60 30 59.9616 -38.4 60.0861 86.1 BCG -5 1 / / 5 5 5.114 114 5.109 109 15 10 14.960 -40 15. 0831 83.1 25 15 25.055 55 25.0578 57.8 35 20 35.028 28 35.031 31 45 25 45.016 16 45.0045 4.5 55 30 54.932 -68 54.9778 -22.2 65 35 64.9053 -94.7 64.9512 -48.8 ABC -5 1 / / 5 5 4.9867 -13.3 5.109 109 15 10 14.96 -40 15.0831 83.1 25 15 25.055 55 25.0578 57.8 35 20 35.0283 28.3 35.0313 31.3 45 25 45.0073 7.3 45.0045 4.5 55 30 54.9806 -19.4 54.9778 -22.2 65 35 64.954 -46 64.9535 -46.5 VI. CONCLUSIONS Performance evaluation of TWFL systems has become an increasingly important issue given its design, manufacturing, sale/purchase, use, upgrade, tuning, etc. In this study, a TWFL model built into a 220kV transmission line that can be easily and reliably simulated on the Matlab Simulink software with 10MHz sampling frequency. The proposed model determines exactly fault types, fault direction and fault location. According to the obtained results, it has been shown that TWFL is more accurate than the traditional impedance-based methods in relay protection. Double-ended method is more accurate than single- ended method. Furthermore, the paper demonstrates a Bewley diagram with the times of the current impulse wave at each terminal and results which are intended to help the operator transition from beginner to experienced professionals. This can assist in the research to eliminate factors of misoperation, contribute substantially to the safe, reliable, economical operation and maintenance of overhead transmission lines. ACKNOWLEDGMENT Authors would like to thank Power Transmission Company No. 2, Viet Nam for allowing the use of fault location equipment for 220kV and 500kV transmission lines used in this study. REFERENCES [1] K. H. Le, P. H. Vu, “A Studying Of Single Ended Fault Locator On SEL Relay”, IETEC’13 Conference, Ho Chi Minh City, Vietnam, December 4-6, 2013 [2] SEL, SEL T400L Ultra high speed transmission line relay traveling wave fault locator high resolution event recorder Instruction Manual, 2018 [3] Kinkei System Corporation, Surge type fault locator system specifications SFL 2000, 2016. [4] S. Parmar, Fault Location Algorithms for Electrical Power Transmission Lines Methodology, Design, and Testing, MSc Thesis, Delft University of Technology, 2015 [5] B. Kasztenny, A. Guzmán, N. Fischer, M. V. Mynam, D. Taylor, “Practical setting considerations for protective relays that use incremental quantities and traveling waves”, 43rd Annual Western Protective Relay Conference, Washington, USA, October 18-20, 2016