# An Improved Traveling-Wave-Based Fault Location Method with Compensating the Dispersion Effect of Traveling Wave in Wavelet Domain

Authors: Huibin Jia
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The fault generated transient traveling waves are wide band signals which cover the whole frequency range. When the frequency characteristic of line parameters is considered, different frequency components of traveling wave will have different attenuation values and wave velocities, which is defined as the dispersion effect of traveling wave. Because of the dispersion effect, the rise or fall time of the wavefront becomes longer, which decreases the singularity of traveling wave and makes it difficult to determine the arrival time and velocity of traveling wave. Furthermore, the dispersion effect seriously affects the accuracy and reliability of fault location. In this paper, a novel double-ended fault location method has been proposed with compensating the dispersion effect of traveling wave in wavelet domain. From the propagation theory of traveling wave, a correction function is established within a certain limit band to compensate the dispersion effect of traveling wave. Based on the determined arrival time and velocity of traveling wave, the fault distance can be calculated precisely by utilizing the proposed method. The simulation experiments have been carried out in ATP/EMTP software, and simulation results demonstrate that, compared with the traditional traveling-wave fault location methods, the proposed method can significantly improve the accuracy of fault location. Moreover, the proposed method is insensitive to different fault conditions, and it is adaptive to both transposed and untransposed transmission lines well.


---

## Body

## 1. Introduction

Power systems have grown rapidly over the last few decades, and the number and length of transmission lines increased. Transmission lines are exposed in the field, and, especially in the mountains and hilly terrains, they are prone to failure. In this scenario, a fast and accurate fault location technique is essential to reduce the restoration time of power systems, which is important with respect to technical and economic issues. Therefore, the study and development of fault location have been motivated since the 1950s [1]. The traveling-wave-based fault location method has been proposed by lots of researchers because these are insensitive to load flow, transition resistance, wiring ways, and series compensation.The traveling-wave-based fault location methods for transmission lines can generally be classified as single- and double-ended methods in terms of their different ways of obtaining the fault information. For many years, the single-ended traveling-wave-based methods were recognized by utilities as a good way to overcome the drawbacks of impedance-based approaches [2]. However, these methods commonly have problems of distinguishing between traveling waves reflected from the fault point and from power system terminals, which decreases the reliability of the fault location method [3]. Therefore, double-ended traveling-wave-based methods have been reported for overcoming these drawbacks. The double-ended method employs the data from the two ends of transmission lines; these data are synchronized by using GPS [4]. And the double-ended method has higher reliability and accuracy than the single-ended one.For the traveling-wave-based fault location method, the accuracy of fault location lies in the arrival time and velocity of traveling wave. Several methods have been proposed to determine the arrival time of traveling wave [5–10]. Wavelet transform has strong time-frequency analysis capability, which can effectively improve the accuracy of the singular value detection. Therefore, it has been firstly used for determining the arrival time of traveling wave, which is represented by the wavelet modulus maxima point [5]. Hence, several researchers have exploited the continuous wavelet transform to extract the arrival time of traveling wave [6–10]. As a matter of fact, it can be easier to extract the singular value point when mother wavelet is similar to traveling wave. Therefore, the mother wavelet has been extracted from the fault transient signal to improve the accuracy of fault location in distribution network [11, 12]. All these methods described above determine the arrival time of wavefront with singular value detection algorithm without considering different frequency components of traveling wave having different arrival times.On the other hand, the wave velocity also directly affects the accuracy of fault location. Some researchers have proposed some methods which are insensitive to the velocity of traveling wave [13–15]. The method proposed in [13] has exploited the three-terminal synchronized transient data, which increases the cost of fault location and decreases the reliability of fault location. The method proposed in [14] is excessively dependent on the first consecutive transient wavefronts of traveling wave, which makes it difficult to distinguish the wavefront of the reflected traveling wave. The method in [15] needs a dedicated communication system to locate the faults, but unfortunately the latency of the dedicated communication system is uncertain, which increases the error of fault location. Currently, the light speed is selected in most traveling-wave-based fault location methods. A few researchers have made significant attempts to determine the wave velocity. Reference [16] calculated the wave velocity with the arrival time of initial wave at nonfault phase when an out-area fault happens. Reference [17] calculated the wave velocity with high-frequency signal generated by autoreclosure. There is a certain difference between the fault signal and the analyzed high-frequency signal. Lots of experts have made discussions and attempts on the problem of determining the wave velocity, but there is still not a brilliant solution. Moreover, these methods have not also considered the different frequency components having different propagation velocities.When a fault happens, the generated traveling wave contains a lot of frequency components, which can be viewed as very important fault information. Because of the frequency-dependent parameters of transmission line, each component has different velocity and attenuation. This phenomenon is defined as traveling-wave dispersion effect [18–20], which distorts the wavefront of traveling wave and lengthens the fall or rise time of the wavefront. Therefore, it is difficult to accurately determine the arrival time and velocity of traveling wave. In [18–20], the authors just analyzed the influence of the dispersion effect on the traveling-wave-based fault location and did not provide the perfect solution. In this paper, a correction method for dispersion effect of traveling wave is proposed. The correction method can shorten the fall or rise time of the wavefront and enhance the singularity of traveling wave. Thus, the accuracy and reliability of fault location can be improved significantly.The paper is organized as follows. Section2 introduces the principle of traveling-wave-based fault location. Section 3 analyzes dispersion effect of traveling wave. Section 4 proposes a correction method for dispersion effect of traveling wave in detail. Section 5 shows the analysis of simulations and results. Finally, Section 6 is the conclusion.


## 2. Principle of Traveling-Wave-Based Fault Location

### 2.1. Traveling-Wave Theory

When a fault happens on a transmission line, the generated voltage and current surges will travel towards both ends of the line. This is equivalent to generating a virtual voltage source at the fault point, as shown in Figure1. The traveling wave is driven by the voltage source and travels towards both ends of the line. Its velocity is close to the velocity of light.Figure 1
Traveling-wave theory.Refraction and reflection will occur at the fault point, ends of lines, and other discontinuity points. Thus, the generated reflected and refracted wave will propagate along the transmission line, as shown in Figure2. The wavefront is detected when the initial traveling wave first reaches M-side and N-side, denoted as M1 and N1. Then, the signal moves to the fault point from M-side (N-side) of the line and returns to the fault point F. The refracted wavefront N1 will reach M-side, marked as M2. And the reflected wavefront M1 will go back to M-side again, marked as M3. M1 and N1 are used for double-ended fault location in this paper.Figure 2
Refection and refraction of traveling wave.


### 2.2. Double-Ended Traveling-Wave-Based Method

In the double-ended traveling-wave-based method, the distance between a fault point and measuring point atM and N terminals can be obtained according to(1)lM=12vtM-tN+L(2)lN=12vtN-tM+L,where lM and lN are the distance away from the measuring points at M and N terminals to the fault point, respectively. L is the total length of the transmission line. v is the propagation velocity of traveling wave.


### 2.3. Phase-to-Modal Transformation

There is an electromagnetic coupling effect between transmission lines. Therefore, by using a modal transformation matrix, we can decompose the traveling wave between phases into several independent modes. In the case of transposed transmission lines, the parameters of transmission line are balanced and the modal transformation can be defined as constant and real [21], such as Karrenbauer transformation and Clarke transformation. In this paper, Karrenbauer transformation is adopted for decoupling the transposed transmission lines. The Karrenbauer transformation is as follows: (3)U0UαUβ=13111-1010-11UAUBUC,where UA, UB, and UC are the phase voltages, respectively. U0 is the ground mode component, and Uα and Uβ are two independent mode components, respectively. Each independent mode component has different velocity and attenuation.For untransposed transmission line, the phase-to-modal transformation matrix is frequency-dependent and unsymmetrical [22]. Fortunately, it is shown that the frequency effect on the modal transformation matrix can be neglected when the frequency is high enough [23]. Therefore, transformation matrix can be approximated as a real and frequency-independent one. In the case of untransposed transmission line, the Karrenbauer transformation is also adopted in this paper.


### 2.4. Wavelet Transform

Wavelet transform (WT) is a common mathematical tool for digital signal processing. WT has been widely applied in lots of fields, such as time series analysis, speech processing, digital image processing, and power system transient analysis [24]. The continuous wavelet transform (CWT) of a signal f(t) is the integral of the product between f(t) and the daughter wavelet. The daughter wavelet is constructed by a mother wavelet, which is dilated with a scale parameter a and translated by b. Therefore, the mother wavelet of CWT can be defined as (5). Hence,(4)CWTft;a,b=∫-∞∞ftΨa,b∗tdt(5)Ψa,b∗t=1aΨt-ba,where Ψa,b(t) is the mother wavelet and ∗ is the complex conjugate operation. a is a scale factor and b is the transition factor. Wavelet transform has a time-frequency resolution. A short window is used for high frequencies while a long window is for low frequencies. Sharp signal transitions create large-amplitude wavelet coefficients. Large wavelet coefficients can detect and measure short high-frequency variations because they have narrow time localization at high frequency. Therefore, WT is very attractive for the analysis of transient signals [25].Wavelet transform is regarded as a filter. The filter of wavelet transform can be employed as follows:(6)gt=∫-∞∞fτht-τdτ(7)Gω=Fω·Hω,where h(t) is the pulse response of the filter and H(ω) is the frequency response of the filter. Comparing (4) and (6), h(t) can be described as(8)ht=a-1/2ψ-ta¯.H ( ω ) is also represented as(9)Hω=a·Ψaω¯.From the analysis above, wavelet transform can be calculated by using Fast Fourier Transform (FFT) and Inverse Fast Fourier Transform (IFFT). Therefore, wavelet transform of traveling wave can be employed by using the following equation:(10)CWTa=IFFTfωHω,where f(ω) is the calculation result of FFT for traveling wave and CWTa is the result of wavelet transform under the scale a and it is used for fault location. In this paper, the continuous wavelet transform is utilized to extract the wavefront of traveling wave and determine the arrival time of traveling wave.


## 2.1. Traveling-Wave Theory

When a fault happens on a transmission line, the generated voltage and current surges will travel towards both ends of the line. This is equivalent to generating a virtual voltage source at the fault point, as shown in Figure1. The traveling wave is driven by the voltage source and travels towards both ends of the line. Its velocity is close to the velocity of light.Figure 1
Traveling-wave theory.Refraction and reflection will occur at the fault point, ends of lines, and other discontinuity points. Thus, the generated reflected and refracted wave will propagate along the transmission line, as shown in Figure2. The wavefront is detected when the initial traveling wave first reaches M-side and N-side, denoted as M1 and N1. Then, the signal moves to the fault point from M-side (N-side) of the line and returns to the fault point F. The refracted wavefront N1 will reach M-side, marked as M2. And the reflected wavefront M1 will go back to M-side again, marked as M3. M1 and N1 are used for double-ended fault location in this paper.Figure 2
Refection and refraction of traveling wave.


## 2.2. Double-Ended Traveling-Wave-Based Method

In the double-ended traveling-wave-based method, the distance between a fault point and measuring point atM and N terminals can be obtained according to(1)lM=12vtM-tN+L(2)lN=12vtN-tM+L,where lM and lN are the distance away from the measuring points at M and N terminals to the fault point, respectively. L is the total length of the transmission line. v is the propagation velocity of traveling wave.


## 2.3. Phase-to-Modal Transformation

There is an electromagnetic coupling effect between transmission lines. Therefore, by using a modal transformation matrix, we can decompose the traveling wave between phases into several independent modes. In the case of transposed transmission lines, the parameters of transmission line are balanced and the modal transformation can be defined as constant and real [21], such as Karrenbauer transformation and Clarke transformation. In this paper, Karrenbauer transformation is adopted for decoupling the transposed transmission lines. The Karrenbauer transformation is as follows: (3)U0UαUβ=13111-1010-11UAUBUC,where UA, UB, and UC are the phase voltages, respectively. U0 is the ground mode component, and Uα and Uβ are two independent mode components, respectively. Each independent mode component has different velocity and attenuation.For untransposed transmission line, the phase-to-modal transformation matrix is frequency-dependent and unsymmetrical [22]. Fortunately, it is shown that the frequency effect on the modal transformation matrix can be neglected when the frequency is high enough [23]. Therefore, transformation matrix can be approximated as a real and frequency-independent one. In the case of untransposed transmission line, the Karrenbauer transformation is also adopted in this paper.


## 2.4. Wavelet Transform

Wavelet transform (WT) is a common mathematical tool for digital signal processing. WT has been widely applied in lots of fields, such as time series analysis, speech processing, digital image processing, and power system transient analysis [24]. The continuous wavelet transform (CWT) of a signal f(t) is the integral of the product between f(t) and the daughter wavelet. The daughter wavelet is constructed by a mother wavelet, which is dilated with a scale parameter a and translated by b. Therefore, the mother wavelet of CWT can be defined as (5). Hence,(4)CWTft;a,b=∫-∞∞ftΨa,b∗tdt(5)Ψa,b∗t=1aΨt-ba,where Ψa,b(t) is the mother wavelet and ∗ is the complex conjugate operation. a is a scale factor and b is the transition factor. Wavelet transform has a time-frequency resolution. A short window is used for high frequencies while a long window is for low frequencies. Sharp signal transitions create large-amplitude wavelet coefficients. Large wavelet coefficients can detect and measure short high-frequency variations because they have narrow time localization at high frequency. Therefore, WT is very attractive for the analysis of transient signals [25].Wavelet transform is regarded as a filter. The filter of wavelet transform can be employed as follows:(6)gt=∫-∞∞fτht-τdτ(7)Gω=Fω·Hω,where h(t) is the pulse response of the filter and H(ω) is the frequency response of the filter. Comparing (4) and (6), h(t) can be described as(8)ht=a-1/2ψ-ta¯.H ( ω ) is also represented as(9)Hω=a·Ψaω¯.From the analysis above, wavelet transform can be calculated by using Fast Fourier Transform (FFT) and Inverse Fast Fourier Transform (IFFT). Therefore, wavelet transform of traveling wave can be employed by using the following equation:(10)CWTa=IFFTfωHω,where f(ω) is the calculation result of FFT for traveling wave and CWTa is the result of wavelet transform under the scale a and it is used for fault location. In this paper, the continuous wavelet transform is utilized to extract the wavefront of traveling wave and determine the arrival time of traveling wave.


## 3. Analysis of Traveling-Wave Dispersion Effect

### 3.1. Analysis of Dispersion Effect in Frequency Domain

The traveling wave has a lot of frequency components. As a matter of fact, the distributed parameters of transmission line are frequency-dependent, which results in different frequency components having different velocities and attenuation values. This phenomenon is defined as the dispersion effect of traveling wave [18].There are two important parameters for the traveling-wave dispersion effect. One is the propagation coefficient, and the other is phase velocity. The propagation coefficient of them-modal component is calculated as follows:(11)γmω=Rm+jωLmGm+jωCm,where Rm, Lm, Gm, and Cm are the m-modal resistance, inductance, conductance, and capacitance, respectively, which are frequency-dependent because of the influence of skin effect [26]. These distributed parameters can be calculated by using Carson formulation according to geometrical parameters of transmission line. It can be observed from Figure 3 that the modal resistance increases with the increase of the frequency, and the modal inductance decreases with the increase of frequency.Figure 3
Frequency-dependent parameters of transmission line.
(a)
Resistance varied with frequency (b)
Reactance varied with frequency (c)
Attenuation coefficient varied with frequency (d)
Propagation velocity varied with frequencyThe propagation coefficient is also described as follows:(12)γmω=αmω+jβmω,where αm is the attenuation coefficient and βm is the phase coefficient. They are also frequency-dependent parameters. The phase velocity can be calculated by phase coefficient βm, which is defined as(13)Vpm=ωβm.The modal attenuation coefficient and phase velocity, which vary with frequencies, are shown in Figures3(c) and 3(d). The frequency ranges from 1 Hz to 1 MHz.It can be seen from Figure3 that the ground-modal parameters are more affected by the frequency than the line-modal parameters. The line-modal component is widely used for the traveling-wave-based fault location. Nevertheless, the characteristics of frequency-dependent parameters should not be ignored for improving the accuracy of fault location. For line-modal parameters, attenuation coefficients and propagation velocities increase along with the increase of frequency. For attenuation coefficients, the change of the low frequency is less than that of high frequency; for propagation velocities, the change of the low frequency is greater than that of the high frequency. The characteristic of frequency-dependent parameters makes the rise or fall time of wavefront spread out, which decreases the singularity of transient traveling wave.


### 3.2. Analysis of Dispersion Effect in Time Domain

In frequency domain, attenuation and wave velocity of traveling wave increase along with frequency. In time domain, the arrival times of different frequency components are not the same. The high-frequency component will arrive at the measuring point first of all, and the low-frequency component will be delayed much more in reaching the measuring point than the high-frequency component. Therefore, all of the detected wavefronts at different measuring points are not ideal step signals; these rise times or fall times get long, which can been seen from Figure4. In the figure, the signal at the fault point is an ideal step one. At the measuring point 150 kilometers away from the fault point, the fall time of transient traveling wave becomes larger than that at the fault point. At the measuring point 300 kilometers away from the fault point, the fall time of the transient traveling wave becomes larger than that at 150 kilometers. The fall time will become larger and larger with the increase of the distance between the measuring point and the fault point. Thus, the singularity of transient traveling wave gradually decreases. In summary, as a result of the dispersion effect of traveling wave, the wavefront is distorted, which makes it more difficult to be detected. Furthermore, the dispersion effect of traveling wave will decrease the accuracy and reliability of fault location.Figure 4
α-mode wavefronts at different measuring points.


## 3.1. Analysis of Dispersion Effect in Frequency Domain

The traveling wave has a lot of frequency components. As a matter of fact, the distributed parameters of transmission line are frequency-dependent, which results in different frequency components having different velocities and attenuation values. This phenomenon is defined as the dispersion effect of traveling wave [18].There are two important parameters for the traveling-wave dispersion effect. One is the propagation coefficient, and the other is phase velocity. The propagation coefficient of them-modal component is calculated as follows:(11)γmω=Rm+jωLmGm+jωCm,where Rm, Lm, Gm, and Cm are the m-modal resistance, inductance, conductance, and capacitance, respectively, which are frequency-dependent because of the influence of skin effect [26]. These distributed parameters can be calculated by using Carson formulation according to geometrical parameters of transmission line. It can be observed from Figure 3 that the modal resistance increases with the increase of the frequency, and the modal inductance decreases with the increase of frequency.Figure 3
Frequency-dependent parameters of transmission line.
(a)
Resistance varied with frequency (b)
Reactance varied with frequency (c)
Attenuation coefficient varied with frequency (d)
Propagation velocity varied with frequencyThe propagation coefficient is also described as follows:(12)γmω=αmω+jβmω,where αm is the attenuation coefficient and βm is the phase coefficient. They are also frequency-dependent parameters. The phase velocity can be calculated by phase coefficient βm, which is defined as(13)Vpm=ωβm.The modal attenuation coefficient and phase velocity, which vary with frequencies, are shown in Figures3(c) and 3(d). The frequency ranges from 1 Hz to 1 MHz.It can be seen from Figure3 that the ground-modal parameters are more affected by the frequency than the line-modal parameters. The line-modal component is widely used for the traveling-wave-based fault location. Nevertheless, the characteristics of frequency-dependent parameters should not be ignored for improving the accuracy of fault location. For line-modal parameters, attenuation coefficients and propagation velocities increase along with the increase of frequency. For attenuation coefficients, the change of the low frequency is less than that of high frequency; for propagation velocities, the change of the low frequency is greater than that of the high frequency. The characteristic of frequency-dependent parameters makes the rise or fall time of wavefront spread out, which decreases the singularity of transient traveling wave.


## 3.2. Analysis of Dispersion Effect in Time Domain

In frequency domain, attenuation and wave velocity of traveling wave increase along with frequency. In time domain, the arrival times of different frequency components are not the same. The high-frequency component will arrive at the measuring point first of all, and the low-frequency component will be delayed much more in reaching the measuring point than the high-frequency component. Therefore, all of the detected wavefronts at different measuring points are not ideal step signals; these rise times or fall times get long, which can been seen from Figure4. In the figure, the signal at the fault point is an ideal step one. At the measuring point 150 kilometers away from the fault point, the fall time of transient traveling wave becomes larger than that at the fault point. At the measuring point 300 kilometers away from the fault point, the fall time of the transient traveling wave becomes larger than that at 150 kilometers. The fall time will become larger and larger with the increase of the distance between the measuring point and the fault point. Thus, the singularity of transient traveling wave gradually decreases. In summary, as a result of the dispersion effect of traveling wave, the wavefront is distorted, which makes it more difficult to be detected. Furthermore, the dispersion effect of traveling wave will decrease the accuracy and reliability of fault location.Figure 4
α-mode wavefronts at different measuring points.


## 4. Correction Method of Traveling-Wave Dispersion Effect

### 4.1. Correction Method of Traveling-Wave Dispersion Effect

The propagation theory of forward and reverse traveling wave is shown in Figure5. The direction of the “forward” and “reverse” is relative. The component that leaves the terminal is defined as “forward”; the component that enters the terminal from the line is defined as “reverse.” The length of transmission line between M and N terminal is l.Figure 5
The propagation theory of forward and reverse traveling wave.According to the propagation equation of single conductor line, the voltagesUM and UN at M and N terminal can be obtained as follows:(14)UN=FN+BN=FN+FNe-γωlUM=FM+BM=FM+FNe-γωl,where e-γ(ω)l is the propagation function. Compared with (14), each reverse traveling wave can be obtained as follows:(15)BM=FNe-γωlBN=FMe-γωl.From signal and system theory, the reverse traveling waveBM is the frequency response of the forward traveling wave FN; similarly, the reverse traveling wave BM is the frequency response of the forward traveling wave FM. For lossy transmission line, γ(ω) varies with frequency. It distorts the wavefront of traveling wave and decreases its singularity. Obviously, γ(ω) causes the dispersion effect of traveling wave.Likewise, the initial surge of traveling wave, which is used for double-ended fault location method, is the frequency response of the signal at the fault point. Suppose that the distance between the fault point and the measuring point isL and the transient traveling wave at the fault point is referred to as fω in frequency domain, which is an ideal step signal in time domain as shown in Figure 4. Thus, the traveling wave at the measuring points can be described as f(ω)e-γ(ω)L in frequency domain from the analysis above, and this wave is distorted by e-γ(ω)L. The traveling wave at the measuring points can be multiplied by A(ω), and the traveling wave at the fault point with the constant delay can be obtained, and it is described as follows: (16)fωe-γωLAω=fωe-jωτ(17)τ=Lυ,where υ is set as the velocity at 50 Hz. f(ω)e-jωτ is the perfect solution for the traveling-wave-based fault location, because it is an ideal step signal and it has perfect singularity and constant delay. In theory, f(ω)e-jωτ is very suitable for fault location. Consequently, the correction function in frequency domain is obtained as follows: (18)Aω=Q-1ω=eγωL-jωτ.For fault location,A(ω) multiplies the receiving traveling wave at the measuring point in frequency domain; the perfect signal f(ω)e-jωτ can be obtained and it is used for the traveling-wave-based fault location. That is to say, all of the frequency components have the same velocity and attenuation to the frequency component at 50 Hz by using the correction function. When the fault happens on the transmission line, the forward or reverse traveling wave travels towards the measuring points. The distorted wavefront of the initial traveling wave is corrected by the correction function. All frequency components of the initial traveling wave will be postponed for some constant time to reach the measuring point and all of them have constant attenuation values, which make it easy to accurately determine the arrival time and wave velocity. By the correction method proposed in this paper, the singularity of the receiving transient traveling wave is significantly enhanced.


### 4.2. Implementation of the Correction Method

When a fault happens, the three-phase voltage atM and N measuring terminals can be obtained. Firstly, approximate fault distance can be obtained by using wavelet transform, which is described in Sections 2.2 and 2.4. The approximate fault distance is used for calculating the correction function; secondly, we select several discrete frequencies and calculate the propagation function of each frequency with R, L, and C parameters by using Carson formulation, and then the correction function A(ω) can be obtained by using (18); thirdly, the corrected traveling wave can be obtained according to (16) in frequency domain and the fault location procedure is implemented by wavelet transform again. The specific steps are shown as follows.Step 1.
The three-phase voltage is obtained when a fault happens.Step 2.
The transient three-phase voltages are first decoupled into their independent modal components by using (3), and then α-mode traveling wave, uMt and uNt, is used for fault location.Step 3.
u M ω and uN(ω) of α-mode traveling wave can be obtained by Fast Fourier Transform (FFT).Step 4.
At the scalea, the frequency response of the mother wavelet function, H(ω), can be calculated by using (9).Step 5.
At the scalea, the decomposition result of wavelet transform can be obtained by IFFT according to (10); this result can be described as follows:(19)WTauMt=IFFTuMω·HωWTauNt=IFFTuNω·Hω.Step 6.
Based on the double-ended fault location theory described in Section2.2, approximate fault distance, lM and lN, can be obtained.Step 7.
Calculate the propagation function based on Carson formulation, and then calculate the corresponding corrected functionAM(ω) and AN(ω) for uM(ω) and uN(ω), respectively.Step 8.
At the scalea of wavelet transform, the corrected decomposition result of wavelet transform can be obtained by IFFT; this result can be shown as follows:(20)WTuM′t=IFFTuMω·AMω·HωWTuNt=IFFTuNω·ANω·Hω.Step 9.
Based on the double-ended fault location theory, accurate fault distance,lM′ and lN′, can be obtained.From the analysis above, it can be observed that the wavelet-transform-based fault location method is described in frequency domain from Step3 to Step 6. The method is implemented again from Step 7 to Step 9, and the traveling wave is corrected in frequency domain.


## 4.1. Correction Method of Traveling-Wave Dispersion Effect

The propagation theory of forward and reverse traveling wave is shown in Figure5. The direction of the “forward” and “reverse” is relative. The component that leaves the terminal is defined as “forward”; the component that enters the terminal from the line is defined as “reverse.” The length of transmission line between M and N terminal is l.Figure 5
The propagation theory of forward and reverse traveling wave.According to the propagation equation of single conductor line, the voltagesUM and UN at M and N terminal can be obtained as follows:(14)UN=FN+BN=FN+FNe-γωlUM=FM+BM=FM+FNe-γωl,where e-γ(ω)l is the propagation function. Compared with (14), each reverse traveling wave can be obtained as follows:(15)BM=FNe-γωlBN=FMe-γωl.From signal and system theory, the reverse traveling waveBM is the frequency response of the forward traveling wave FN; similarly, the reverse traveling wave BM is the frequency response of the forward traveling wave FM. For lossy transmission line, γ(ω) varies with frequency. It distorts the wavefront of traveling wave and decreases its singularity. Obviously, γ(ω) causes the dispersion effect of traveling wave.Likewise, the initial surge of traveling wave, which is used for double-ended fault location method, is the frequency response of the signal at the fault point. Suppose that the distance between the fault point and the measuring point isL and the transient traveling wave at the fault point is referred to as fω in frequency domain, which is an ideal step signal in time domain as shown in Figure 4. Thus, the traveling wave at the measuring points can be described as f(ω)e-γ(ω)L in frequency domain from the analysis above, and this wave is distorted by e-γ(ω)L. The traveling wave at the measuring points can be multiplied by A(ω), and the traveling wave at the fault point with the constant delay can be obtained, and it is described as follows: (16)fωe-γωLAω=fωe-jωτ(17)τ=Lυ,where υ is set as the velocity at 50 Hz. f(ω)e-jωτ is the perfect solution for the traveling-wave-based fault location, because it is an ideal step signal and it has perfect singularity and constant delay. In theory, f(ω)e-jωτ is very suitable for fault location. Consequently, the correction function in frequency domain is obtained as follows: (18)Aω=Q-1ω=eγωL-jωτ.For fault location,A(ω) multiplies the receiving traveling wave at the measuring point in frequency domain; the perfect signal f(ω)e-jωτ can be obtained and it is used for the traveling-wave-based fault location. That is to say, all of the frequency components have the same velocity and attenuation to the frequency component at 50 Hz by using the correction function. When the fault happens on the transmission line, the forward or reverse traveling wave travels towards the measuring points. The distorted wavefront of the initial traveling wave is corrected by the correction function. All frequency components of the initial traveling wave will be postponed for some constant time to reach the measuring point and all of them have constant attenuation values, which make it easy to accurately determine the arrival time and wave velocity. By the correction method proposed in this paper, the singularity of the receiving transient traveling wave is significantly enhanced.


## 4.2. Implementation of the Correction Method

When a fault happens, the three-phase voltage atM and N measuring terminals can be obtained. Firstly, approximate fault distance can be obtained by using wavelet transform, which is described in Sections 2.2 and 2.4. The approximate fault distance is used for calculating the correction function; secondly, we select several discrete frequencies and calculate the propagation function of each frequency with R, L, and C parameters by using Carson formulation, and then the correction function A(ω) can be obtained by using (18); thirdly, the corrected traveling wave can be obtained according to (16) in frequency domain and the fault location procedure is implemented by wavelet transform again. The specific steps are shown as follows.Step 1.
The three-phase voltage is obtained when a fault happens.Step 2.
The transient three-phase voltages are first decoupled into their independent modal components by using (3), and then α-mode traveling wave, uMt and uNt, is used for fault location.Step 3.
u M ω and uN(ω) of α-mode traveling wave can be obtained by Fast Fourier Transform (FFT).Step 4.
At the scalea, the frequency response of the mother wavelet function, H(ω), can be calculated by using (9).Step 5.
At the scalea, the decomposition result of wavelet transform can be obtained by IFFT according to (10); this result can be described as follows:(19)WTauMt=IFFTuMω·HωWTauNt=IFFTuNω·Hω.Step 6.
Based on the double-ended fault location theory described in Section2.2, approximate fault distance, lM and lN, can be obtained.Step 7.
Calculate the propagation function based on Carson formulation, and then calculate the corresponding corrected functionAM(ω) and AN(ω) for uM(ω) and uN(ω), respectively.Step 8.
At the scalea of wavelet transform, the corrected decomposition result of wavelet transform can be obtained by IFFT; this result can be shown as follows:(20)WTuM′t=IFFTuMω·AMω·HωWTuNt=IFFTuNω·ANω·Hω.Step 9.
Based on the double-ended fault location theory, accurate fault distance,lM′ and lN′, can be obtained.From the analysis above, it can be observed that the wavelet-transform-based fault location method is described in frequency domain from Step3 to Step 6. The method is implemented again from Step 7 to Step 9, and the traveling wave is corrected in frequency domain.


## 5. Simulations and Results Analysis

### 5.1. The Model of Simulation

A model for 500 kV transposed transmission line is constructed in ATP/EMTP, as shown in Figure6(a). The geometry structure of transmission line is shown in Figure 6(b). The classical J. Marti model is adopted for analyzing the frequency characteristic of transmission lines. The voltage measuring terminal is at each transformer and is used to acquire the fault voltage data. The total length of the transmission line is 414 km. The fault distance is 100 km. Besides, an A-phase grounding fault happens.Figure 6
Model of transmission line.
(a)
500 kV transmission line model (b)
Geometry structure of transmission lineIn ATP/EMTP software, the sample rate is set as 1 MHz. The traveling-wave data generated in ATP software is imported into MATLAB software. In MATLAB, the correction method is implemented. In this paper,α-mode component is used for fault location. Due to the fact that the velocity is dependent on frequency, the α-mode velocity of traveling wave is calculated at 50 Hz, which is 2.9724 × 105 km/s.


### 5.2. Simulation Waveform Analysis

When we select several discrete frequencies and calculate the propagation function of each frequency withR, L, and C parameters by using Carson formulation, the propagation function A(ω) can be obtained by using (19).When we obtainA(ω), it is used for correcting the receiving traveling wave by using (21). The corrected traveling wave is shown in Figure 7. At the moment, the transmission line is grounded in A phase, the transition resistance is 10 Ω, and the fault inception angle is 106°.Figure 7
Corrected waveforms with the distance of 345 km.As can be seen from Figure7, the fall time of the corrected traveling wave becomes shorter than of the uncorrected one, which enhances the singularity of the traveling wave. That is to say, the corrected traveling wave shows more singularity than the uncorrected one. The correction method proposed in this paper makes the wavefront easier to be detected by the extraction algorithms, such as wavelet transform and HHT. In this paper, the comparison analysis is also implemented in wavelet field. As shown in Figure 8, the modulus maximum of uncorrected wavelet component is less than that of corrected wavelet component.Figure 8
Comparisons of the wavelet coefficients.


### 5.3. Results of Fault Location

In this paper, we compare the fault location accuracy between the traditional method and the proposed method. The traditional method is based on the wavelet transform (Morlet wavelet function) and it is not corrected by the correction method proposed in this paper. In order to test the applicability of the proposed method, various fault conditions are simulated, respectively. In the simulated experiments, the error of fault location is calculated as a performance index by the following equation:(21)error=l′-lS×100%,where l′ is the calculated fault distance, l is the actual fault distance, and S is the length of total transmission line.


#### 5.3.1. Performance of Different Transition Resistance Values

The transition resistance directly affects the voltage amplitude of the initial traveling wave. In this paper, we set the fault inception angle to analyze the errors of fault location, when the transition resistance varies from 10 Ω to 500 Ω. The distance between the fault point and the measuring terminal is 100 km, the fault type is A-G, and the fault inception angle is 106°. The simulated results are shown in Tables1 and 2. As can be seen from the table, the location error of the proposed method becomes smaller than that of the traditional one, no matter what the transition resistance is.Table 1
Error under different transition resistance values (transposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.3 0.10 100.23 0.07 50 99.50 0.18 99.60 0.13 100 99.34 0.22 99.62 0.12 500 98.70 0.43 99.30 0.23Table 2
Error under different transition resistance values (untransposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.57 0.19 100.30 0.10 50 101.20 0.40 100.56 0.19 100 98.75 0.42 99.18 0.27 500 101.02 0.34 100.55 0.18


#### 5.3.2. Performance of Different Inception Angles

The fault inception angle is one of the most effective parameters which influence the accuracy of the traveling-wave-based fault location. In order to test the influence of fault inception angle, the simulated cases are constructed with the fault inception angle varied from 18° to 90°. At the moment, the fault type is A-G, the transition resistance is 10 Ω, and the fault distance is 100 kilometers. The experiment results are shown in Tables3 and 4. The proposed method also shows great adaptability to the different inception angles.Table 3
Error under different inception angles (transposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.60 0.20 100.42 0.07 36° 100.31 0.10 100.22 0.10 72° 100.86 0.29 99.96 0.18 90° 100.80 0.27 100.78 0.26Table 4
Error under different inception angles (untransposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.85 0.28 100.50 0.17 36° 100.62 0.21 100.40 0.13 72° 99.33 0.22 99.60 0.13 90° 100.56 0.19 100.33 0.11


#### 5.3.3. Performance of Different Fault Types

To investigate the effect of different fault types based on the proposed method in this paper, several experiments have been implemented and experiment results are shown in Tables5 and 6. In these experiments, all transition resistances are 10 Ω, the fault inception angle is 106°, and the fault distance is 100 kilometers.Table 5
Error under different fault types (transposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.30 0.10 100.23 0.08 B-G 99.50 0.17 99.63 0.12 BC 99.21 0.26 100.45 0.15 BC-G 100.87 0.29 100.14 0.05 ABC-G 99.25 0.25 99.45 0.18Table 6
Error under different fault types (untransposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.57 0.19 100.30 0.10 B-G 99.54 0.15 99.64 0.12 BC 101.05 0.35 100.85 0.28 BC-G 100.78 0.26 100.23 0.08 ABC-G 100.50 0.16 100.26 0.09


#### 5.3.4. Performance of Antinoise

To investigate the performance of antinoise under different levels, several experiments have been implemented. In these experiments, the transition resistance is 10 Ω with BC-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Since noise has no meaningful information about fault location, it would be rational to define a threshold to reject the noise-associated WT coefficient [8]. Simulation results in Table 7 reveal that the additional noise has little influence on traveling-wave-based fault location.Table 7
Error under different noise levels (transposed circuit).
Noise (unit: dBw) The traditional method The proposed method Location results Error Location results Error 10 100.78 0.26 100.14 0.05 20 100.78 0.26 100.15 0.05 30 99.21 0.2633 100.30 0.10 40 100.76 0.2533 100.21 0.07


#### 5.3.5. Performance Analysis of Geometry Parameters

The production error of the towers supporting transmission line is very low, and it is within 1 cm. The geometrical parameters of the corner tower and the tower in substation are sometimes different from other towers along the transmission line. The number of the corner towers and the towers in substation is very small and they can be omitted. Therefore, from the point of view of theoretical calculation, we assume that the relative position between line conductors is unchanged. On the other hand, the height of line conductors will sometimes increase when the transmission line passes through a road, village, city, mountain, and so forth. The change is very complicated; it is very difficult to model the real transmission line. In the paper, for analyzing the applicability of the proposed method, the length of the transmission line which is supported by the changed tower is 20 km. In these experiments, the transition resistance is 10 Ω with A-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Simulation results in Table8 reveal that the proposed method can also improve the accuracy of fault location when the height of tower varies.Table 8
Error under different tower heights (transposed circuit).
Tower height (mile) The traditional method The proposed method Location results Error Location results Error 0 100.60 0.20 100.14 0.05 +5 100.47 0.16 100.22 0.07 +10 99.18 0.27 100.45 0.15 +20 100.80 0.27 100.55 0.18


## 5.1. The Model of Simulation

A model for 500 kV transposed transmission line is constructed in ATP/EMTP, as shown in Figure6(a). The geometry structure of transmission line is shown in Figure 6(b). The classical J. Marti model is adopted for analyzing the frequency characteristic of transmission lines. The voltage measuring terminal is at each transformer and is used to acquire the fault voltage data. The total length of the transmission line is 414 km. The fault distance is 100 km. Besides, an A-phase grounding fault happens.Figure 6
Model of transmission line.
(a)
500 kV transmission line model (b)
Geometry structure of transmission lineIn ATP/EMTP software, the sample rate is set as 1 MHz. The traveling-wave data generated in ATP software is imported into MATLAB software. In MATLAB, the correction method is implemented. In this paper,α-mode component is used for fault location. Due to the fact that the velocity is dependent on frequency, the α-mode velocity of traveling wave is calculated at 50 Hz, which is 2.9724 × 105 km/s.


## 5.2. Simulation Waveform Analysis

When we select several discrete frequencies and calculate the propagation function of each frequency withR, L, and C parameters by using Carson formulation, the propagation function A(ω) can be obtained by using (19).When we obtainA(ω), it is used for correcting the receiving traveling wave by using (21). The corrected traveling wave is shown in Figure 7. At the moment, the transmission line is grounded in A phase, the transition resistance is 10 Ω, and the fault inception angle is 106°.Figure 7
Corrected waveforms with the distance of 345 km.As can be seen from Figure7, the fall time of the corrected traveling wave becomes shorter than of the uncorrected one, which enhances the singularity of the traveling wave. That is to say, the corrected traveling wave shows more singularity than the uncorrected one. The correction method proposed in this paper makes the wavefront easier to be detected by the extraction algorithms, such as wavelet transform and HHT. In this paper, the comparison analysis is also implemented in wavelet field. As shown in Figure 8, the modulus maximum of uncorrected wavelet component is less than that of corrected wavelet component.Figure 8
Comparisons of the wavelet coefficients.


## 5.3. Results of Fault Location

In this paper, we compare the fault location accuracy between the traditional method and the proposed method. The traditional method is based on the wavelet transform (Morlet wavelet function) and it is not corrected by the correction method proposed in this paper. In order to test the applicability of the proposed method, various fault conditions are simulated, respectively. In the simulated experiments, the error of fault location is calculated as a performance index by the following equation:(21)error=l′-lS×100%,where l′ is the calculated fault distance, l is the actual fault distance, and S is the length of total transmission line.


### 5.3.1. Performance of Different Transition Resistance Values

The transition resistance directly affects the voltage amplitude of the initial traveling wave. In this paper, we set the fault inception angle to analyze the errors of fault location, when the transition resistance varies from 10 Ω to 500 Ω. The distance between the fault point and the measuring terminal is 100 km, the fault type is A-G, and the fault inception angle is 106°. The simulated results are shown in Tables1 and 2. As can be seen from the table, the location error of the proposed method becomes smaller than that of the traditional one, no matter what the transition resistance is.Table 1
Error under different transition resistance values (transposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.3 0.10 100.23 0.07 50 99.50 0.18 99.60 0.13 100 99.34 0.22 99.62 0.12 500 98.70 0.43 99.30 0.23Table 2
Error under different transition resistance values (untransposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.57 0.19 100.30 0.10 50 101.20 0.40 100.56 0.19 100 98.75 0.42 99.18 0.27 500 101.02 0.34 100.55 0.18


### 5.3.2. Performance of Different Inception Angles

The fault inception angle is one of the most effective parameters which influence the accuracy of the traveling-wave-based fault location. In order to test the influence of fault inception angle, the simulated cases are constructed with the fault inception angle varied from 18° to 90°. At the moment, the fault type is A-G, the transition resistance is 10 Ω, and the fault distance is 100 kilometers. The experiment results are shown in Tables3 and 4. The proposed method also shows great adaptability to the different inception angles.Table 3
Error under different inception angles (transposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.60 0.20 100.42 0.07 36° 100.31 0.10 100.22 0.10 72° 100.86 0.29 99.96 0.18 90° 100.80 0.27 100.78 0.26Table 4
Error under different inception angles (untransposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.85 0.28 100.50 0.17 36° 100.62 0.21 100.40 0.13 72° 99.33 0.22 99.60 0.13 90° 100.56 0.19 100.33 0.11


### 5.3.3. Performance of Different Fault Types

To investigate the effect of different fault types based on the proposed method in this paper, several experiments have been implemented and experiment results are shown in Tables5 and 6. In these experiments, all transition resistances are 10 Ω, the fault inception angle is 106°, and the fault distance is 100 kilometers.Table 5
Error under different fault types (transposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.30 0.10 100.23 0.08 B-G 99.50 0.17 99.63 0.12 BC 99.21 0.26 100.45 0.15 BC-G 100.87 0.29 100.14 0.05 ABC-G 99.25 0.25 99.45 0.18Table 6
Error under different fault types (untransposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.57 0.19 100.30 0.10 B-G 99.54 0.15 99.64 0.12 BC 101.05 0.35 100.85 0.28 BC-G 100.78 0.26 100.23 0.08 ABC-G 100.50 0.16 100.26 0.09


### 5.3.4. Performance of Antinoise

To investigate the performance of antinoise under different levels, several experiments have been implemented. In these experiments, the transition resistance is 10 Ω with BC-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Since noise has no meaningful information about fault location, it would be rational to define a threshold to reject the noise-associated WT coefficient [8]. Simulation results in Table 7 reveal that the additional noise has little influence on traveling-wave-based fault location.Table 7
Error under different noise levels (transposed circuit).
Noise (unit: dBw) The traditional method The proposed method Location results Error Location results Error 10 100.78 0.26 100.14 0.05 20 100.78 0.26 100.15 0.05 30 99.21 0.2633 100.30 0.10 40 100.76 0.2533 100.21 0.07


### 5.3.5. Performance Analysis of Geometry Parameters

The production error of the towers supporting transmission line is very low, and it is within 1 cm. The geometrical parameters of the corner tower and the tower in substation are sometimes different from other towers along the transmission line. The number of the corner towers and the towers in substation is very small and they can be omitted. Therefore, from the point of view of theoretical calculation, we assume that the relative position between line conductors is unchanged. On the other hand, the height of line conductors will sometimes increase when the transmission line passes through a road, village, city, mountain, and so forth. The change is very complicated; it is very difficult to model the real transmission line. In the paper, for analyzing the applicability of the proposed method, the length of the transmission line which is supported by the changed tower is 20 km. In these experiments, the transition resistance is 10 Ω with A-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Simulation results in Table8 reveal that the proposed method can also improve the accuracy of fault location when the height of tower varies.Table 8
Error under different tower heights (transposed circuit).
Tower height (mile) The traditional method The proposed method Location results Error Location results Error 0 100.60 0.20 100.14 0.05 +5 100.47 0.16 100.22 0.07 +10 99.18 0.27 100.45 0.15 +20 100.80 0.27 100.55 0.18


## 5.3.1. Performance of Different Transition Resistance Values

The transition resistance directly affects the voltage amplitude of the initial traveling wave. In this paper, we set the fault inception angle to analyze the errors of fault location, when the transition resistance varies from 10 Ω to 500 Ω. The distance between the fault point and the measuring terminal is 100 km, the fault type is A-G, and the fault inception angle is 106°. The simulated results are shown in Tables1 and 2. As can be seen from the table, the location error of the proposed method becomes smaller than that of the traditional one, no matter what the transition resistance is.Table 1
Error under different transition resistance values (transposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.3 0.10 100.23 0.07 50 99.50 0.18 99.60 0.13 100 99.34 0.22 99.62 0.12 500 98.70 0.43 99.30 0.23Table 2
Error under different transition resistance values (untransposed circuit).
Transition resistance The traditional method The proposed method Location results Error Location results Error 10 100.57 0.19 100.30 0.10 50 101.20 0.40 100.56 0.19 100 98.75 0.42 99.18 0.27 500 101.02 0.34 100.55 0.18


## 5.3.2. Performance of Different Inception Angles

The fault inception angle is one of the most effective parameters which influence the accuracy of the traveling-wave-based fault location. In order to test the influence of fault inception angle, the simulated cases are constructed with the fault inception angle varied from 18° to 90°. At the moment, the fault type is A-G, the transition resistance is 10 Ω, and the fault distance is 100 kilometers. The experiment results are shown in Tables3 and 4. The proposed method also shows great adaptability to the different inception angles.Table 3
Error under different inception angles (transposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.60 0.20 100.42 0.07 36° 100.31 0.10 100.22 0.10 72° 100.86 0.29 99.96 0.18 90° 100.80 0.27 100.78 0.26Table 4
Error under different inception angles (untransposed circuit).
Fault inception angle The traditional method The proposed method Location results Error Location results Error 18° 100.85 0.28 100.50 0.17 36° 100.62 0.21 100.40 0.13 72° 99.33 0.22 99.60 0.13 90° 100.56 0.19 100.33 0.11


## 5.3.3. Performance of Different Fault Types

To investigate the effect of different fault types based on the proposed method in this paper, several experiments have been implemented and experiment results are shown in Tables5 and 6. In these experiments, all transition resistances are 10 Ω, the fault inception angle is 106°, and the fault distance is 100 kilometers.Table 5
Error under different fault types (transposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.30 0.10 100.23 0.08 B-G 99.50 0.17 99.63 0.12 BC 99.21 0.26 100.45 0.15 BC-G 100.87 0.29 100.14 0.05 ABC-G 99.25 0.25 99.45 0.18Table 6
Error under different fault types (untransposed circuit).
Fault type The traditional method The proposed method Location results Error Location results Error A-G 100.57 0.19 100.30 0.10 B-G 99.54 0.15 99.64 0.12 BC 101.05 0.35 100.85 0.28 BC-G 100.78 0.26 100.23 0.08 ABC-G 100.50 0.16 100.26 0.09


## 5.3.4. Performance of Antinoise

To investigate the performance of antinoise under different levels, several experiments have been implemented. In these experiments, the transition resistance is 10 Ω with BC-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Since noise has no meaningful information about fault location, it would be rational to define a threshold to reject the noise-associated WT coefficient [8]. Simulation results in Table 7 reveal that the additional noise has little influence on traveling-wave-based fault location.Table 7
Error under different noise levels (transposed circuit).
Noise (unit: dBw) The traditional method The proposed method Location results Error Location results Error 10 100.78 0.26 100.14 0.05 20 100.78 0.26 100.15 0.05 30 99.21 0.2633 100.30 0.10 40 100.76 0.2533 100.21 0.07


## 5.3.5. Performance Analysis of Geometry Parameters

The production error of the towers supporting transmission line is very low, and it is within 1 cm. The geometrical parameters of the corner tower and the tower in substation are sometimes different from other towers along the transmission line. The number of the corner towers and the towers in substation is very small and they can be omitted. Therefore, from the point of view of theoretical calculation, we assume that the relative position between line conductors is unchanged. On the other hand, the height of line conductors will sometimes increase when the transmission line passes through a road, village, city, mountain, and so forth. The change is very complicated; it is very difficult to model the real transmission line. In the paper, for analyzing the applicability of the proposed method, the length of the transmission line which is supported by the changed tower is 20 km. In these experiments, the transition resistance is 10 Ω with A-G fault, the fault inception angle is 106°, and the fault distance is 100 kilometers. Simulation results in Table8 reveal that the proposed method can also improve the accuracy of fault location when the height of tower varies.Table 8
Error under different tower heights (transposed circuit).
Tower height (mile) The traditional method The proposed method Location results Error Location results Error 0 100.60 0.20 100.14 0.05 +5 100.47 0.16 100.22 0.07 +10 99.18 0.27 100.45 0.15 +20 100.80 0.27 100.55 0.18


## 6. Conclusions

In this paper, the dispersion characteristic of traveling wave is analyzed in time and frequency domain, respectively. When traveling wave travels along transmission line, the fall or rise time of the wavefront will become long with the increase of propagation distance. Thus, the singularity of the transient wavefront decreases. In this paper, a novel double-ended fault location method has been proposed to overcome the dispersion effect of traveling wave. In the method, a correction algorithm for overcoming the dispersion effect of traveling wave enhances the singularity of the transient traveling wave. The proposed method is tested under various experiment conditions, such as different fault distances, different transition resistances, and different fault inception angles. The simulation experiments demonstrate that the proposed method is better than the traditional traveling-wave-based fault location method. Furthermore, the novel method is suitable for both transposed and untransposed transmission lines. All of the advantages prove that the proposed method is available. However, there are several other parameters which affect the slope of traveling wave, such as the impedance characteristic of traveling-wave measurement equipment and shunt reactor connected to the line. Future work will consider those parameters.



---

*Source: 1019591-2017-02-08.xml*