Fault Model and Travelling Wave Matching Based Single Terminal Fault Location Algorithm for T-Connection Transmission Line: A Yunnan Power Grid Study

Due to the complex structure of the T-connection transmission lines, it is extremely difficult to identify the reflected travelling wave from the fault point and that from the connection point by the measurement from only one terminal. According to the characteristics of the structure of the T-connection transmission line, the reflection of the travelling wave within the line after the failure of different sections in T-connection transmission line are analyzed. Based on the lattice diagram of the travelling wave, the sequence of travelling waves detected at the measuring terminal varies with the fault distance and the faulty section. Moreover, the sequence of travelling waves detected in one terminal is unique at each faulty section. This article calculates the arrival time of travelling waves of fault points at different locations in different sections to form the collection of the travelling wave arrival time sequence. Then the sequence of travelling waves of the new added fault waveforms is extracted to compare with the sequences in the collection for the faulty section identification and fault location. This proposed method can accurately locate the fault with different fault types, fault resistances and system impedances by only single-terminal fault data. Both Power Systems Computer Aided Design/ Electromagnetic Transients including DC (PSCAD/EMTDC) and actual measurement data are implemented to verify the effectiveness of this method.


Introduction
With the global emphasis on environmental protection and the recent acceleration of urbanization, an increasing number of clean energy sources, e.g., wind turbines and photovoltaics, have been merged into the Chinese power grid through 110 kV transmission lines [1][2][3][4]. Particularly, the Yunnan power grid is the regional power grid with the highest proportion of clean energy power generation in China, which is an important power base for China's significant national strategy of power transmission from west to east [5][6][7][8]. Therefore, 110 kV grids are the significant part of the Yunnan power grid for large-scale clean energy integration and power supply for cities [9][10][11][12]. Due to the land shortage in cities, the T-connection transmission line, which requires less area coverage and economic cost than multiple two-terminal transmission lines, is usually selected in 110 kV grids in China [13][14][15].

Travelling Wave Propagation Characteristics of Faults in Section MT
We assume that the distance from the fault point F to the measurement terminal M is F x , l1 is the line distance from terminal M to connection point T, l2 is the line distance from terminal P to connection point T, while l3 is the distance from terminal N to connection point T, furthermore l1 > l3 > l2. It is supposed in this paper that t0, tF, tT, tN, tP are the times of the travelling waves arriving terminal M, which are initial travelling wave, reflected wave surge from fault point F, reflected travelling wave from connection point T, reflected travelling wave from terminal N and reflected travelling wave from terminal P. The time of travelling wave surges arriving terminal M are listed in Table 1.

Travelling Wave Propagation Characteristics of Faults in Section MT
We assume that the distance from the fault point F to the measurement terminal M is x F , l 1 is the line distance from terminal M to connection point T, l 2 is the line distance from terminal P to connection point T, while l 3 is the distance from terminal N to connection point T, furthermore l 1 > l 3 > l 2 .
It is supposed in this paper that t 0 , t F , t T , t N , t P are the times of the travelling waves arriving terminal M, which are initial travelling wave, reflected wave surge from fault point F, reflected travelling wave from connection point T, reflected travelling wave from terminal N and reflected travelling wave from terminal P. The time of travelling wave surges arriving terminal M are listed in Table 1. Table 1. Travelling wave types with their reach times.

Time
Travelling Wave Surge Types t 0 Initial travelling wave t F Reflected travelling wave from F t T Reflected travelling wave from T t N Reflected travelling wave from N t P Reflected travelling wave from P Under this assumption, it is necessary to analyze the possible quantitative relationship between x F , l 1 − x F , and l 1 + l 2 − x F . There are three situations among x F , l 1 − x F , and l 1 + l 2 − x F , which is (1) If the fault occurs in the range from M to T and the fault distance satisfies x F < l 1 − x F < l 1 + l 2 − x F , the travelling wave propagation path of the T-connection transmission line is shown in Figure 2.
The time and the propagation paths of travelling waves arriving the terminal M are listed in Table 2.
The time and the propagation paths of travelling waves arriving the terminal M are listed in Table 2.
Therefore, the reflected travelling wave from the fault point F reaches the measuring terminal M earlier than the reflected travelling wave from the connection point T. Thus, it can be clearly seen from Figure 2 that t0<tF<tT<tP<tN.
The distance between the initial travelling wave and the reflected travelling wave from the fault point F is 2x F . The relationship of t 0 , t F , t T , t P and t N can be determined as follows: Therefore, the reflected travelling wave from the fault point F reaches the measuring terminal M earlier than the reflected travelling wave from the connection point T. Thus, it can be clearly seen from Figure 2 that t 0 <t F <t T <t P <t N .
(2) If the fault occurs in the range from the terminal M to the connection point T with relationship of l 1 − x F < x F < l 1 + l 2 − x F , the travelling wave propagation path of the T-connection transmission line is shown in Figure 3. Measurement terminal The time of travelling waves arriving the terminal M and the propagation paths are listed in Table 3.
The time of travelling waves arriving the terminal M and the propagation paths are listed in Table 3.  The distance between the initial travelling wave and the reflected travelling wave from connection point T is 2(l 1 − x F ). The relationship of t 0 , t T , t F , t P and t N can be determined as follows: can be concluded that the reflected travelling wave from the connection point T reaches the measuring terminal M earlier than the reflected travelling wave from the fault point F. From Equation (2), it is easy to determine t 0 <t T <t F <t P <t N .
(3) If the fault occurs in the range from terminal M to connection point T within the relationship of l 1 − x F < l 1 + l 2 − x F < x F , the travelling wave propagation path of the T-connection transmission line is shown in Figure 4. Measurement terminal  Table 4.
The relationship of t0, tT, tP, tF and tN can be determined as follows: The times of travelling waves arriving terminal M are listed in Table 4. The distance between the initial travelling wave and the reflected travelling wave from connection point T is 2(l 1 − x F ). The relationship of t 0 , t T , t P , t F and t N can be determined as follows: the reflected travelling wave from the connection point T reaches the measuring terminal M earlier than the reflected travelling wave from the fault point F. From Equation (3), it is easy to determine t 0 <t T <t P <t F <t N .

Travelling Wave Propagation Characteristics of Faults in Section TP
When the fault occurs at the section TP, the fault travelling wave propagation path is shown in Figure 5. Measurement terminal The times of travelling waves arriving terminal M are listed in Table 5. In the condition shown in Figure 5, the distance between the initial travelling wave and the reflected travelling wave from terminal P is 1 2 . The relationship of t0, tP, tT, tN and tF can be determined as follows: Under this condition, it is not difficult to see that can be concluded that the reflected travelling wave from the terminal P reaches the measuring terminal M earlier than the reflected travelling wave from the terminal N. Besides it is easy to determine t0<tP<tN<tT<tF. The times of travelling waves arriving terminal M are listed in Table 5. In the condition shown in Figure 5, the distance between the initial travelling wave and the reflected travelling wave from terminal P is 2(l 1 + l 2 − x F ). The relationship of t 0 , t P , t T , t N and t F can be determined as follows: M earlier than the reflected travelling wave from the terminal N.

Travelling Wave Propagation Characteristics of Faults in Section TN
When the fault occurs at the section of TN, the situation is more complex than mentioned above. The fault travelling wave propagation path is shown in Figure 6. Measurement terminal The distance between the initial travelling wave and the reflected travelling wave from terminal N is The relationship of t0, tN, tP, tT, and tF can be determined as follows: In this condition, we cannot determine which of the reflected travelling waves from terminal P and the reflected travelling wave from terminal N reaches the terminal M first, because the relationship between  Table 6. In this table, the positional relationship of tN and tP may be reversed. Therefore, it is easy to determine t0<tP<tN<tT<tF, or t0<tN<tP<tT<tF.

Analysis of Connection Point T Fault Travelling Wave Propagation Characteristics
If the fault occurs at the connection point T, this case is the most special of all cases. The reflected travelling wave from the fault point T and the reflected travelling wave from the The distance between the initial travelling wave and the reflected travelling wave from terminal N is 2(l 1 + l 3 − x F ). The relationship of t 0 , t N , t P , t T , and t F can be determined as follows: In this condition, we cannot determine which of the reflected travelling waves from terminal P and the reflected travelling wave from terminal N reaches the terminal M first, because the relationship between t P − t N = 2(l 2 + x F − l 1 − l 3 ) and 0 cannot be judged. However, this does not affect the judgment of the fault section, because only when the section TN is faulty, the two travelling waves that finally arrived at terminal M are the reflected travelling wave from the connection point T and the reflected travelling wave from fault point F. Moreover, no matter whether the reflected travelling wave from terminal N reaches the terminal M first or the reflected travelling wave from terminal P first reaches the terminal M, their time relationship with the arrival time of the initial wave will not change, as shown in Equation (5). The times of travelling waves arriving terminal M are listed in Table 6. In this table, the positional relationship of t N and t P may be reversed. Therefore, it is easy to determine t 0 <t P <t N <t T <t F , or t 0 <t N <t P <t T <t F .

Analysis of Connection Point T Fault Travelling Wave Propagation Characteristics
If the fault occurs at the connection point T, this case is the most special of all cases. The reflected travelling wave from the fault point T and the reflected travelling wave from the connection point T is coincided. The fault travelling wave propagation path is shown in Figure 7.
Energies 2017, 10, xFOR PEER REVIEW 9 of 23 In this case, it is only necessary to establish the relationship between the initial travelling wave and the reflected travelling wave from the connection point T, the reflected travelling wave from the terminal P, and the reflected travelling wave from terminal N. The distance between the initial travelling wave and the reflected travelling wave from terminal P is 2 2l . The distance between the initial travelling wave and the reflected travelling wave from terminal N is 3 2l . The distance between the initial travelling wave and the reflected travelling wave from the connection point T (the reflected travelling wave from fault point F) is F 2 x . The relationship of t0, tP, tN, tT(tF) can be determined as follows: The times of travelling waves arriving terminal M are listed in Table 7.

T(F)→M→T(F)→M F(T)→M→F(T)→M
In this case, it is not difficult to see that By analyzing the above fault conditions, it can be found that when the fault occurs on any point of the T-connection transmission line, there is a specific travelling wave sequence, which corresponds to the location of the fault point. In other words, a specific travelling wave sequence corresponds to only one fault point on the T-connection transmission line. In order to identify the faulty section and locate the fault point, a collection of the sequences of travelling wave arriving In this case, it is only necessary to establish the relationship between the initial travelling wave and the reflected travelling wave from the connection point T, the reflected travelling wave from the terminal P, and the reflected travelling wave from terminal N. The distance between the initial travelling wave and the reflected travelling wave from terminal P is 2l 2 . The distance between the initial travelling wave and the reflected travelling wave from terminal N is 2l 3 . The distance between the initial travelling wave and the reflected travelling wave from the connection point T (the reflected travelling wave from fault point F) is 2x F . The relationship of t 0 , t P , t N , t T (t F ) can be determined as follows: The times of travelling waves arriving terminal M are listed in Table 7. Table 7. Travelling wave types with their reach times.

Time Travelling Wave Types Propagation Paths
In this case, it is not difficult to see that t N − t P = 2(l 3 − l 2 )/v > 0, therefore, it can be concluded that the reflected travelling wave from the terminal P reaches the measuring terminal M earlier than the reflected travelling wave from the terminal N.
By analyzing the above fault conditions, it can be found that when the fault occurs on any point of the T-connection transmission line, there is a specific travelling wave sequence, which corresponds to the location of the fault point. In other words, a specific travelling wave sequence corresponds to only one fault point on the T-connection transmission line. In order to identify the faulty section and locate Energies 2020, 13, 1506 9 of 22 the fault point, a collection of the sequences of travelling wave arriving time is established. In this collection, each data is simulated as a fault point in the transmission line with different fault location to obtain the arrival moment of reflected travelling wave from fault point F, reflected travelling wave from connection point T, reflected travelling wave from terminal N, and reflected travelling wave from terminal P. Finally, the sequence of the travelling waves of the researched data is applied to compare with the data in the collection.

Signal Singularity Detection Method of Multiple Recursive SVD
SVD is a method of numerical analysis with clear physical significance, so it has been applied in many fields [31][32][33]. SVD can represent a relatively complex matrix with a relatively simple matrix. These simple matrices describe the important characteristics of the matrix.

Principle of Singular Value Decomposition
SVD is the decomposition of a matrix H of length m × n into an orthogonal matrix U = (u 1 , u 2 , . . . , u n ) ∈ R m×m and an orthogonal matrix V = (v 1 , v 2 , . . . , v n ) ∈ R n×n . For the generally complex valued elements' matrix H, SVD is obtained by factorizing it into a product of three matrices as follows in Equation (7): where S is an m × n matrix whose elements are zero except possibly along the main diagonal (the singular values of H), S = [diag( σ 1 , σ 2 , · · ·, σ n )] or it is transposed, depending on whether m≤n or m>n, S ∈ R m×n , q = min(m, n), while σ 1 ≥ σ 2 ≥ · · · ≥ σ q > 0, σ i (i = 1, 2, · · · , q) are the singular values of matrix H [34].

n-Division Recursive SVD Algorithm
Although the dichotomous recursive SVD method has a good detection effect on singular points, but it has a poor effect on signal noise reduction. Based on this feature, this paper improves the dichotomous recursion method and proposes an n-division recursion SVD method that can effectively suppress noise. The principle of the n-recursive SVD method is as follows: A discrete signal X= (x 1 , x 2 , . . . , x N ) is sent, N is the original signal length, and an n × (N + 1 − n)-dimensional Hankel matrix H is constructed for the original signal, H ∈ R (N+1−n)×n : By performing SVD on the H matrix, n singular values can be obtained. In the matrix H, starting from the second row, each row lags one data later than the previous row, and the correlation between the two rows of data is high. H can be rewritten as a vector u i and v i according to Equation (9): in Equation (9), so H 1 ∈ R n×(N+1−n) . The structure of H 1 is similar to H, but not identical, H 1 is just a detail matrix that can best reflect the H information of the original signal, and the specific expression of H 1 shown in Equation (10): In the matrix H 1 , H 1 (n,1)=H 1 (n-1,2)= . . . = H 1 (1, n) should be satisfied. In the same way, T to H n = σ in u in v in T can be calculated as above [35].
Each layer of the n-division recursion SVD is a further fine characterization of the original signal, but at the same time the energy of the detailed component signal of each layer also decreases accordingly. It is particularly important to be able to adaptively determine an optimal decomposition layer according to the target signal. Referring to the concept of information entropy, a singular entropy is introduced to represent the uncertainty of the variable. The calculation method as shown in Equation (11): In Equation (11), m is the number of singular values, and σ i / m k=1 σ k is the weight of the i-th singular value to all singular values. The singular entropy reflects the uncertainty of the signal singularity energy distribution. The fuzzier the signal singularity is, the more concentrated the energy distribution is on a few singular values, and the smaller the signal singular entropy; on the contrary, the more obvious the signal singularity, the more dispersed the energy, and the larger the signal singular entropy. After the n-level multi-scale SVD, n detail component matrices and one approximate component matrix can be obtained. The new matrix composed of n + 1 matrices can represent all the information of the original signal. According to the SVD theory, the singular value of a matrix can reflect the entire amount of information of the matrix. When the number of decomposition layers is increased to a certain extent, if the singular entropy of the new matrix does not increase significantly, the number of decomposition layers is considered to be optimal.
We define the n-th layer singular entropy increment ∆E n as: where n≥2. If ∆E obviously approaches to 0, it is not necessary to continue the decomposition, and the number of decomposition layers reaches the optimal. By decomposing the travelling wave data with different layers, the singular entropy results are tabulated in Table 8. It can be concluded from Table 8 that the singular entropy decreases with the increase of the decomposition layer. However, the decreasing of the singular entropy is much slower if the decomposition layer is larger than 3. Therefore, the decomposition is chosen as 3.

Single-Terminal Fault Location Scheme for T-Connection Transmission Line
The selection of the fault time window plays a vital role in the realization of travelling wave fault location. When the time window is too short, the travelling wave required for fault location may not be within the selected time window. Loss of important information required for fault location will directly lead to failure. If the time window selected is too long, the reflection waves of various nodes that are not needed for fault location and more noise interference will appear in the time window. The extra travelling wave may lead to incorrect results in fault location. In this method, full consideration should be given to the time window containing the fault point reflected travelling wave, reflected travelling wave from the connection point, reflected travelling wave from the terminal N and reflected travelling wave from the terminal P. The length of the time window determined by Equation (13): In Equation (13), max represents the maximum value, l 1 , l 2 , and l 3 represent the distances from the connection point T to the three terminals, the bus at the terminal M is the measurement terminal, and v is the travelling wave velocity.
A fault point x F is set every 100 m from the terminal M to calculate the time of different reflected travelling wave arriving the terminal M. Then, all the calculated sequences are applied to form a collection as the fault cases. After extracting the sequence of travelling waves by SVD, the specific fault data are compared with the sequences of travelling waves in the collection. Owing to the uniqueness of the travelling wave arrival sequence, the faulty section and the fault point can be determined by the comparison of cases in the collection. In order to measure the differences of the sequence of the specific data and the data in the collection, δ is calculated by the sum of time differences.
In Equation (14), t i represents the time difference between the travelling wave moments calculation and the original fault signal after n-division recursive SVD, and i represents the different travelling waves that reach the measurement terminal M, as shown in Table 9.  Time difference of reflected travelling waves from N of the specific data and the data in the collection δ, which represents the similarity of the sequences, decreases with the increasing of the similarity of two sequences.
When the matching algorithm is satisfied, the matched fault distance x F is determined, while t 0 , t F are determined by x F . The fault distance of the reflected travelling wave corresponding to the fault point F of the n-recursion SVD method calculated as shown in Equation (15): Energies 2020, 13, 1506

of 22
In Equation (15) presents the fault distance calculated by the n-recursive SVD calibration wave; t 0 is the arrival time of the travelling wave from the fault point, and t F is the arrival moment of the reflected travelling wave from the fault point F. Finally, by averaging the distance calculated by the full-line matching algorithm and the fault distance calculated by the n-recursive recursion method, the final fault location result is obtained. We used the hypothesis result x F and the waveform calculation result L n to obtain the average value, which can make the fault point location result more accurate. The algorithm flowchart is shown in Figure 8.

Simulation Model
PSCAD/EMTDC of version 4.5 [36][37][38][39] with sampling rate of 1 MHz is implemented to simulate the faults in the T-connection transmission lines as demonstrated in Figure 1, and its parameters are given in Table 10. The lengths of the T-connection transmission line are set to l1 = 50 km, l2 = 40 km, and l3 = 45 km.

Simulation Model
PSCAD/EMTDC of version 4.5 [36][37][38][39] with sampling rate of 1 MHz is implemented to simulate the faults in the T-connection transmission lines as demonstrated in Figure 1, and its parameters are given in Table 10. The lengths of the T-connection transmission line are set to l 1 = 50 km, l 2 = 40 km, and l 3 = 45 km.

Single-Phase Ground Fault Simulation
A phase-to-ground (AG) fault occurs at section l 1 with fault distance x F =30 km and 20 Ω fault resistance at 0.202 s, this simulation corresponds to the l 1 − x F < x F < l 1 + l 2 − x F case in Section 2.1. The simulation current waveform and the result of SVD are illustrated in Figure 9.     It can be seen in the Figure 9 that the proposed method accurately detects and calibrates singular points in the transient waveform. Assume the first reflected travelling wave is from each feasible impedance discontinuity, the sequences of the travelling wave of the assumptions are compared with the result of SVD as shown in Figure 10. It can be seen that only the correct matched aligned all singular value points. It can be seen in the Figure 9 that the proposed method accurately detects and calibrates singular points in the transient waveform. Assume the first reflected travelling wave is from each feasible impedance discontinuity, the sequences of the travelling wave of the assumptions are compared with    The matching result is that the fault occurred in l1 section, and the fault distance is km, which is quite similar to the actual fault position.

Lightning Fault Simulation
The travelling wave arrive sequence shown in Figure 11, it can calculate the fault location

Lightning Fault Simulation
A lightning induced Phase A ground fault (AG) occurs 75 km away from the terminal M at 0.204 ms. This simulation corresponds to the l 1 − x F < x F < l 1 + l 2 − x F case in Section 2.3. The lightning current model adopts the double exponential model [40][41][42]. The fault phase current measured at the measuring terminal of the T-connection transmission line shown in Figure 12a. The fault phase data in Figure 12a subjected to n-division recursive SVD singularity calibration, and the result shown in Figure 12b. The results of the combination of the full-line matching algorithm and SVD shown in Figure 13. The result of the calculation is that the fault occurred within l3 of the fault section, and the fault distance is F x = 75.1 km.  The results of the combination of the full-line matching algorithm and SVD shown in Figure 13. The result of the calculation is that the fault occurred within l3 of the fault section, and the fault distance is F x = 75.1 km.  The travelling wave sequence shown in Figure 14 can calculate the fault location within the initial fault travelling wave and the reflected travelling wave from fault point F. The fault distance L SVD is 75 km, so the final result is 75.05 km.
Energies 2017, 10, xFOR PEER REVIEW 17 of 23 The travelling wave sequence shown in Figure 14 can calculate the fault location within the initial fault travelling wave and the reflected travelling wave from fault point F. The fault distance LSVD is 75 km, so the final result is 75.05 km. In order to fully verify the reliability of the algorithm, different types of faults (single-phase ground fault, phase-to-phase short-circuit fault, three-phase short-circuit fault, and lightning strike failure) at different fault sections (l1, l2, l3, and connection point T) were simulated for the test. Owing to the page limits, a part of the simulation results is tabulated in Table 11. In order to fully verify the reliability of the algorithm, different types of faults (single-phase ground fault, phase-to-phase short-circuit fault, three-phase short-circuit fault, and lightning strike failure) at different fault sections (l 1 , l 2 , l 3 , and connection point T) were simulated for the test. Owing to the page limits, a part of the simulation results is tabulated in Table 11.

Practical Field Data Test in Yunnan Province
In order to verify the accuracy of the proposed fault location method for the actual fault data, a travelling wave analysis and locating device was established, as illustrated in Figure 15. The device can complete the synchronous sampling of 16 channel records at a sampling rate of 1 MHz.

Practical Field Data Test in Yunnan Province
In order to verify the accuracy of the proposed fault location method for the actual fault data, a travelling wave analysis and locating device was established, as illustrated in Figure 15. The device can complete the synchronous sampling of 16 channel records at a sampling rate of 1 MHz. The geographical map of Yunnan Province is shown in Figure 16a, while a part of the topology of the 110 kV transmission line in Pu'er Power Grid is shown in Figure 16b. Ximahe substation, which is equipped with the established travelling wave analysis and locating device is and 38.83 km away from the connection point T, is assumed as the terminal M. Similarly, Jiangcheng substation, which is 64.83 km away from the connection point, and Laoshizhai substation, which is 14.4 km away from the connection point, are assumed as the terminal N and P, respectively. A fault induced by the transmission line swing is recorded on June 30, 2018. The measured current waveform is shown in Figure 17.

Practical Field Data Test in Yunnan Province
In order to verify the accuracy of the proposed fault location method for the actual fault data, a travelling wave analysis and locating device was established, as illustrated in Figure 15. The device can complete the synchronous sampling of 16 channel records at a sampling rate of 1 MHz. The geographical map of Yunnan Province is shown in Figure 16a, while a part of the topology of the 110 kV transmission line in Pu'er Power Grid is shown in Figure 16b. Ximahe substation, which is equipped with the established travelling wave analysis and locating device is and 38.83 km away from the connection point T, is assumed as the terminal M. Similarly, Jiangcheng substation, which is 64.83 km away from the connection point, and Laoshizhai substation, which is 14.4 km away from the connection point, are assumed as the terminal N and P, respectively. A fault induced by the transmission line swing is recorded on June 30, 2018. The measured current waveform is shown in Figure 17.       After selecting the time window of the fault phase data, the n-division recursive SVD transformation is applied to detect the wavefronts. The result of SVD is compared with all the assumptions as shown in Figure 18. The types of travelling waves are illustrated in Figure 19. After selecting the time window of the fault phase data, the n-division recursive SVD transformation is applied to detect the wavefronts. The result of SVD is compared with all the assumptions as shown in Figure 18. The types of travelling waves are illustrated in Figure 19.  By the n-division recursive SVD for comparison, it can be correctly identified that the faulty section is the terminal P, the fault distance is 97.1 km away from the terminal M. The actual line inspection result is 97.4 km, which is only 0.3 km different from the fault location result.

Conclusion
This paper presents a single terminal fault location algorithm for T-connection transmission lines based on the sequences of travelling waves, which is able to accurate identify the faulty section and locate the fault point. The main contributions of this paper can be summarized into the following three aspects: (1) This paper discusses the sequence of fault induced travelling waves reaching the measuring terminal to manifest the uniqueness of the sequences of different fault points. A reliable single-terminal fault location method for T-connection transmission lines is proposed, which does not need to identify the travelling wave property of the fault current waveform. After selecting the time window of the fault phase data, the n-division recursive SVD transformation is applied to detect the wavefronts. The result of SVD is compared with all the assumptions as shown in Figure 18. The types of travelling waves are illustrated in Figure 19.  By the n-division recursive SVD for comparison, it can be correctly identified that the faulty section is the terminal P, the fault distance is 97.1 km away from the terminal M. The actual line inspection result is 97.4 km, which is only 0.3 km different from the fault location result.

Conclusion
This paper presents a single terminal fault location algorithm for T-connection transmission lines based on the sequences of travelling waves, which is able to accurate identify the faulty section and locate the fault point. The main contributions of this paper can be summarized into the following three aspects: (1) This paper discusses the sequence of fault induced travelling waves reaching the measuring terminal to manifest the uniqueness of the sequences of different fault points. A reliable single-terminal fault location method for T-connection transmission lines is proposed, which does not need to identify the travelling wave property of the fault current waveform. By the n-division recursive SVD for comparison, it can be correctly identified that the faulty section is the terminal P, the fault distance is 97.1 km away from the terminal M. The actual line inspection result is 97.4 km, which is only 0.3 km different from the fault location result.

Conclusions
This paper presents a single terminal fault location algorithm for T-connection transmission lines based on the sequences of travelling waves, which is able to accurate identify the faulty section and locate the fault point. The main contributions of this paper can be summarized into the following three aspects: (1) This paper discusses the sequence of fault induced travelling waves reaching the measuring terminal to manifest the uniqueness of the sequences of different fault points. A reliable single-terminal fault location method for T-connection transmission lines is proposed, which does not need to identify the travelling wave property of the fault current waveform.
(2). This paper uses the n-recursive SVD to obtain accurate singular value information of the original fault current signal, and then finds the assumed fault location that best matches the fault signal by calculating the singular value information.
(3). Simulation results of different case studies and the practical field data have verified that the proposed method is able to accurately identify the faulty section and locate fault point.