A Fault Location Method for Medium Voltage Distribution Network Based on Ground Fault Transfer Device

: The arc suppression device based on ground fault transfer (GFT) has been preliminarily applied in the medium voltage distribution network (MVDN). An accurate travelling wave (TW) fault location method is proposed to extend the use of the ground fault transfer device. D-PMU is used as a travelling wave detection tool to record the transient voltage travelling waves of fault grounding and bus active grounding during arc suppression. Then, the faulty section is identiﬁed through the time difference of travelling wave arrival at the upstream and downstream measurement points. On this basis, the fault location equations of the arrival time and distance of the upstream travelling wave are established, and an accurate fault location method based on the arrival time difference of the travelling wave is proposed. The simulation model is established by PSCAD/EMTDC, and the results show that the method has high location accuracy, and the absolute error is less than 30 m. It is not affected by the TW velocity, the fault conditions, or the distributed power sources.


Introduction
The GFT device was first used in the chemical industry, metallurgy, and other fields [1].Later, because the fault current was too large and accompanied by overvoltage when a single-phase grounding fault occurred in the wind farm, the arc suppression coil alone could not provide a sufficient compensation current, so the neutral grounding mode was changed and the GFT device was introduced.With the increase in MVDN capacity and the popularization and use of cable lines, MVDNs such as wind farms generally encounter the above problems.The voltage arc extinguishing scheme is adopted due to its good arc extinguishing performance, and the arc extinguishing effect of related devices has also been verified in field tests [2].
In Reference [3], the field fault treatment test of the GFT device was carried out, and the fault treatment data were compared with the treatment results of the automatic tracking arc suppression compensation device, which verified the effectiveness of the fault treatment of the voltage arc extinguishing scheme.Reference [4] analyzed the dynamic characteristics of the grounding switch by the equivalent circuit method and optimized the design and parameter analysis of the fast arc suppression device.On this basis, the parameters of grounding intermediate resistance and grounding transformer were optimized in reference [5] and verified by the online operation data of the product.In Reference [6], the current transfer performance and overvoltage suppression ability of the transfer arc extinguishing device were studied, and the influence of load on the transfer arc extinguishing technology was analyzed.Reference [7] identified whether the fault point disappeared based on the difference in electrical quantities in the station for the light load and heavy load of the line, which provided a theoretical basis for the improvement of the control method of the fault transfer device.Correspondingly, the GFT device manufactured by the equipment manufacturer has been applied in the engineering field.For example, the intelligent disposal scheme of a single-phase grounding fault in an MVDN developed by the Shaanxi Electric Power Research Institute of China can effectively solve the problem of error correction and suppress the transient process, and it has good adaptability to high resistance grounding.
At present, the research work on GFT mostly focuses on the feasibility of the voltage arc extinguishing scheme and the transfer grounding device itself, and the fault location method matching with it is less involved.Making full use of fault information to achieve the accurate location of fault points can reduce the scope of manual line patrol and shorten the power outage times of end users [8], which is also considered to be one of the key measures when building an elastic power grid [9].
According to the different principles, the current fault location methods mainly include the impedance method and the travelling wave method.The impedance method establishes a location equation based on line parameters and measurement information and determines the exact distance of the fault by solving the equation [10].Based on the distributed parameter model of the MVDN, Reference [11] constructs the location function with the current measurement information before and after the fault, which further improves the location accuracy.In Reference [12], the electrical parameters of the opposite end and the fault distance are taken as the quantities to be solved, and the parameter identification is carried out in combination with the spectrum information of the transient voltage and current.It is simple and reliable and overcomes the disadvantage that the single-ended location is easily affected by the system impedance in the traditional scheme.Reference [13] proposes an improved algorithm based on phasor analysis and a voltage and current calculation method suitable for an active MVDN, which can withstand certain load fluctuations.
Fault location based on TW information mainly uses the TW velocity and the arrival time of the wave head to determine the fault location.Theoretically, it is not affected by factors such as the system operation mode, fault type, or transition resistance [14,15], and the location accuracy is high.According to the location principle, it can be divided into the single-ended method [16], double-ended method [17], and modulus time difference method [18].The single-ended method uses the time difference between the initial TW of the fault and the reflected TW of the fault point to construct the distance equation.The principle is simple and easy to implement, but the reflected wave head identification is more difficult [15].The double-ended method calculates the fault distance according to the time difference between the initial TW and the MPs at both ends of the line.It only needs to capture the initial TW head but needs to arrange MPs at all branch nodes, and it needs to strictly meet the requirements of clock synchronization [19].
A new TW fault location method is proposed for the MVDN with a GFT device.Combined with the working characteristics of the GFT device, the arrival time of the TW head of the whole network is integrated when the fault is grounded and the bus is grounded.The fault location equation is constructed according to the arrival time difference of the travelling wave in the faulty section, the fault information of the upstream segment, and the line parameters.Because the location equation is independent of the wave velocity, the influence of the TW velocity on the location result is avoided.The simulation results show that the proposed method has high location accuracy, the absolute error is less than 30 m, and it has good location performance in different fault positions and fault conditions.

The Working Principle of the GFT Device
The GFT device [20] is generally implemented in the neutral point ungrounded system, which is connected to the substation bus through the outlet circuit breaker.The main components are the ground transformer, the control center, and the ground switch, which can be separated.The structure diagram is shown in Figure 1.In particular, different from the traditional circuit breaker using mechanical or electromagnetic suction to design the action switch, due to the need for rapid arc extinguishing, the device generally uses a Thomson coil as the action switch and uses electromagnetic repulsion to meet the requirements of the switch execution speed [21].

The Working Principle of the GFT Device
The GFT device [20] is generally implemented in the neutral point ungrounded system, which is connected to the substation bus through the outlet circuit breaker.The main components are the ground transformer, the control center, and the ground switch, which can be separated.The structure diagram is shown in Figure 1.In particular, different from the traditional circuit breaker using mechanical or electromagnetic suction to design the action switch, due to the need for rapid arc extinguishing, the device generally uses a Thomson coil as the action switch and uses electromagnetic repulsion to meet the requirements of the switch execution speed [21].After determining the grounding fault of the system, the closing pulse is sent to the grounding switch of the faulty phase bus bar, so that the faulty phase bus bar is grounded quickly; the fault point is bypassed, so that the voltage of the fault point is clamped to close to 0, and the grounding arc can be effectively extinguished.After a short delay, the switch is disconnected, and whether the fault disappears is judged according to the zero sequence voltage.If it is a transient grounding fault, the fault disappears and the system returns to normal operation.Otherwise, it is judged to be a permanent fault, and the bus grounding switch is closed again.The scheme does not need to consider the problem of system-to-ground capacitive current and compensation capacity, and the arc extinguishing effect is remarkable.At the same time, after the action of the fault transfer device, the arc grounding that originally occurred on the line can be transformed into a metal grounding in the substation to avoid the formation of step voltage outside the station and reduce the risk of personnel accidents.However, the insulation of the line may be broken after a fault occurs.To avoid the recurrence of the fault, fault location should be carried out as soon as possible to realize the isolation of the nearest fault point.
In the process of fault arc suppression, there are two transient processes of fault grounding and faulty phase bus grounding, so the same line also has two types of transient information, before and after, which provides a new idea for the research of fault location.Taking the TW fault location method studied in this paper as an example, when the GFT device acts, it is equivalent to creating a second grounding fault at the busbar.The arrival time of the line grounding TW can be combined with the arrival time of the busbar grounding TW to construct a new location equation.

Accurate Fault Location Method Based on TWATD
Considering the serious attenuation of zero-mode TW, aerial-mode TW is often used in the calculation of the TW location, but the acquisition of the aerial-mode wave velocity also needs to be matched with accurate line parameters.For the MVDN with the same voltage level, when the frequency is high enough, the aerial-mode wave velocity at the same frequency can be regarded as a certain value [22].Therefore, the same high-frequency modulus TW can be used for location calculation to avoid the error caused by inconsistent wave velocity.To achieve this goal, the extracted TWs should be at the same After determining the grounding fault of the system, the closing pulse is sent to the grounding switch of the faulty phase bus bar, so that the faulty phase bus bar is grounded quickly; the fault point is bypassed, so that the voltage of the fault point is clamped to close to 0, and the grounding arc can be effectively extinguished.After a short delay, the switch is disconnected, and whether the fault disappears is judged according to the zero sequence voltage.If it is a transient grounding fault, the fault disappears and the system returns to normal operation.Otherwise, it is judged to be a permanent fault, and the bus grounding switch is closed again.The scheme does not need to consider the problem of system-to-ground capacitive current and compensation capacity, and the arc extinguishing effect is remarkable.At the same time, after the action of the fault transfer device, the arc grounding that originally occurred on the line can be transformed into a metal grounding in the substation to avoid the formation of step voltage outside the station and reduce the risk of personnel accidents.However, the insulation of the line may be broken after a fault occurs.To avoid the recurrence of the fault, fault location should be carried out as soon as possible to realize the isolation of the nearest fault point.
In the process of fault arc suppression, there are two transient processes of fault grounding and faulty phase bus grounding, so the same line also has two types of transient information, before and after, which provides a new idea for the research of fault location.Taking the TW fault location method studied in this paper as an example, when the GFT device acts, it is equivalent to creating a second grounding fault at the busbar.The arrival time of the line grounding TW can be combined with the arrival time of the busbar grounding TW to construct a new location equation.

Accurate Fault Location Method Based on TWATD
Considering the serious attenuation of zero-mode TW, aerial-mode TW is often used in the calculation of the TW location, but the acquisition of the aerial-mode wave velocity also needs to be matched with accurate line parameters.For the MVDN with the same voltage level, when the frequency is high enough, the aerial-mode wave velocity at the same frequency can be regarded as a certain value [22].Therefore, the same high-frequency modulus TW can be used for location calculation to avoid the error caused by inconsistent wave velocity.To achieve this goal, the extracted TWs should be at the same frequency.At present, there are many research results on the detection of the TW head, among which the wavelet transform is widely used.However, the analysis scale and wavelet basis function need to be selected according to the fault characteristics, and the extraction range is limited to a certain frequency band, so it is impossible to extract the TW signal with a single frequency.Based on the Morlet wavelet, the S transform solves the problem that the frequency of the STFT window cannot be adjusted, and it has the advantage of multiresolution of the wavelet transform [23].Therefore, this paper uses the S transform to identify the fault TW head.There are many papers that introduce the algorithm [24,25], so the details are not repeated in this paper.

Faulty Section Identification Method
For the situation in which the line experiences a grounding fault, the MVDN shown in Figure 2 can be taken as an example; the power supply, transformer, and load are omitted in the figure.Among them, MN is the main trunk line, and the branch line BR is connected to the trunk line with branch point B. D k represents the MPs, which are arranged at the end of the network; B m represents the branch point; and R m represents the network endpoint.
among which the wavelet transform is widely used.However, the analysis scale and wavelet basis function need to be selected according to the fault characteristics, and the extraction range is limited to a certain frequency band, so it is impossible to extract the TW signal with a single frequency.Based on the Morlet wavelet, the S transform solves the problem that the frequency of the STFT window cannot be adjusted, and it has the advantage of multi-resolution of the wavelet transform [23].Therefore, this paper uses the S transform to identify the fault TW head.There are many papers that introduce the algorithm [24,25], so the details are not repeated in this paper.

Faulty Section Identification Method
For the situation in which the line experiences a grounding fault, the MVDN shown in Figure 2 can be taken as an example; the power supply, transformer, and load are omitted in the figure.Among them, MN is the main trunk line, and the branch line BR is connected to the trunk line with branch point B. Dk represents the MPs, which are arranged at the end of the network; Bm represents the branch point; and Rm represents the network endpoint.B m+ surrounding area are called the downstream section.Assuming that the fault initial TW detected at each MP is taken from the same higher frequency, the TW velocity of each line is consistent, and the TW transmission equation can be written for the upstream MP: where  Assuming that the branch B m B m+1 experiences a ground fault, the entire network can be divided into two parts according to the location of the fault point, where the end point M and the branch point B m surrounding area are called the upstream section, and the end point N and the branch point B m+1 surrounding area are called the downstream section.Assuming that the fault initial TW detected at each MP is taken from the same higher frequency, the TW velocity of each line is consistent, and the TW transmission equation can be written for the upstream MP: where t D i is the time at which the initial TW of the fault arrives at the MP D i , l D i B m is the distance between the MP D i and the branch point B m , v is the speed at which the travelling wave travels through the line, l fB m is the distance between the branch point B m and the fault point f, and t f is the time at which the fault occurs.By subtracting the adjacent equations in Equation ( 1) and eliminating the time of the fault, we have Similarly, for each MP in the downstream region, the following results can be obtained: where t D j is the time at which the fault initial TW reaches the MP D j , l D j B m+1 is the distance between the MP D j and the branch point B m+1 , and l fB m+1 is the distance between the branch point B m+1 and the fault point f.
The above is the TW transmission situation when a line fault occurs.When the GFT device acts, it is equivalent to artificially creating a second grounding fault at the bus f .As shown in Figure 3, the corresponding measurement points arranged at the end of the network will detect the arrival of the second fault TW.Because fault line selection will be carried out before fault location, only the outlet of the fault line on the bus is drawn in Figure 3.
Similarly, for each MP in the downstream region, the following results can be obtained: where D j t is the time at which the fault initial TW reaches the MP Dj, The above is the TW transmission situation when a line fault occurs.When the GFT device acts, it is equivalent to artificially creating a second grounding fault at the bus f ′ .As shown in Figure 3, the corresponding measurement points arranged at the end of the network will detect the arrival of the second fault TW.Because fault line selection will be carried out before fault location, only the outlet of the fault line on the bus is drawn in Figure 3. Referring to the situation of a line fault, the upstream MPs are where D i t ′ is the time at which the initial TW of the fault reaches the MP Di when the bus is grounded, and Similarly, for each MP in the downstream region, Referring to the situation of a line fault, the upstream MPs are where t D i is the time at which the initial TW of the fault reaches the MP D i when the bus is grounded, and l D 1 D i is the distance between the MP D i and the M-terminal MP D 1 .
Similarly, for each MP in the downstream region, The two TW transmission processes of line grounding and bus grounding are discussed, and the TW transmission in the upstream and downstream sections is discussed, respectively.For each MP in the upstream section, the source directions of the two TWs are different.The first one is from the direction of the fault point, and the second one is from the head end of the line.The TWATD detected by the adjacent MPs is also different.By subtracting the corresponding equations in Equations ( 2) and ( 4), it can be obtained that where ) is the difference between the TWATD at the adjacent MPs in the two grounding processes, and l B i−1 B i is the distance between the two adjacent branch points.
Different from the upstream section, the fault TW detected in the downstream section at the two grounding moments is transmitted from the bus, so the time difference before and after the adjacent MPs is equal.
Comparing Equations ( 6) and ( 7), it can be seen that for the upstream MP of the fault point, the TWATD of the adjacent MPs in the two grounding processes is different, and the difference between the front and back is proportional to the distance between the corresponding branch points of the two MPs.For the downstream section, the TWATD of the adjacent MPs in the two grounding processes is the same, so the faulty section can be identified according to the TWATD of each adjacent MP.

The Accurate Fault Location Method
Assume that the faulty section B m B m+1 has been determined, as shown in Figure 4.The precise location of the fault point is analyzed below.In the process of fault ground TW transmission, the initial TW arrival time of the adjacent MPs in the faulty section is Electronics 2023, 12, x FOR PEER REVIEW 7 of 14 Similarly, subtract the arrival time of the two MPs: Combined with Figure 4, the arrival time difference of the TWs at two adjacent MPs when the bus is grounded is By subtracting Equation ( 9) from Equation ( 10), it can be obtained that 2 Similarly, subtract the arrival time of the two MPs: Combined with Figure 4, the arrival time difference of the TWs at two adjacent MPs when the bus is grounded is By subtracting Equation ( 9) from Equation (10), it can be obtained that According to Equation ( 6) and Equation (11), It can be seen from the above formula that when the faulty section B m B m+1 is determined, the distance l fB m between the fault location f and the branch point B m can be calculated according to the arrival time of the TW at the adjacent MPs in the upstream section and the distance between the corresponding branch points.
On this basis, considering that the recording device may have problems, such as measurement error and inaccurate line lengths, only using the line length of a certain section and the arrival time of the TW at the MP may cause a large location error.Therefore, the accuracy of fault location is improved by comprehensively using the recording data of each MP in the upstream section and the line parameters: It can be seen from Equation ( 13) that the calculation result of the position of the fault point is independent of the TW velocity.Therefore, the fault location method can effectively avoid the influence of the TW velocity on the location result, and the initial TW arrival time of each MP in the upstream section of the two grounding processes is used in the calculation process.The data are highly redundant, which can avoid the influence of the TW arrival time error on the location result.
Because the location accuracy of the proposed method is related to the arrival time of the travelling wave, an advanced D-PMU measuring device is considered for the detection of travelling wave information.The D-PMU voltage travelling wave sensor module is connected to the grounding line of the capacitive device (transformer) and reflects the voltage travelling wave signal on the primary system line by detecting the capacitive current on the grounding line.Due to the derivative relationship between the capacitive current and fault current, the steepness of the fault travelling wave is amplified, so that it has better travelling wave detection ability.The device has been introduced and analyzed in detail in the existing literature [24,26,27].
Based on the above analysis, an accurate fault location method based on the TWATD is proposed.The specific contents are as follows.
Step 1: After the fault occurs, the fault recording data of each MP are uploaded to the fault location system, the fault TW head and the arrival time of each phase are detected according to the above method, and the initial TW arrival time corresponding to each MP D i is recorded in turn in the line grounding time matrix T n×1 (n is the number of MPs).
Step 2: Subtract the adjacent elements t i and t j in the line grounding time matrix T n×1 , and the difference is recorded as ∆t i−j ; then, the line grounding time difference matrix ∆T (n−1)×1 is formed.
Step 3: After the action of the GFT device, referring to Steps 1 and 2, the arrival time of the fault TW corresponding to each MP D i is recorded as t i and recorded in the bus grounding time matrix T n×1 , and then the ∆t i−j and bus grounding time difference matrix ∆T (n−1)×1 are calculated.
Step 4: Starting from the first end of the line i = 1, the corresponding elements ∆t i−j and ∆t i−j in the line grounding time difference matrix ∆T (n−1)×1 and the bus grounding time difference matrix ∆T (n−1)×1 are compared.δt i−j = ∆t i−j − ∆t i−j is compared with the given threshold t set .If δt i−j > t set , the corresponding branches of the MP D i and D j are recorded as the upstream branches, we continue to search backwards until δt i−j < t set , and the upper section of the section is determined as the faulty section.
Step 5: Combined with Equation ( 13), the fault distance is calculated according to the arrival time of the TW at each MP in the upstream section and the section distance.
As to the computational cost of the proposed method, no calculation is involved in Step 1.In addition, Step 2 and Step 3 mainly subtract the time matrix, and Step 4 is a simple logical comparison operation, so the calculation costs of Step 2 to Step 4 are very low.The calculation cost of the proposed method mainly incurred in Step 5. Step 5 mainly uses Equation ( 13) to solve the fault position.It can be found that this formula is additive, and its calculation cost depends on the value of k.The k value is related to the number of outlet branches of the distribution network, and the magnitude order of the k value for the medium voltage distribution network is usually no more than three digits.Considering the advanced computing capability of the current computer, the computational cost of the proposed method can be fully satisfied.

Simulation Analysis
To make the simulation results fit the engineering practice as much as possible, this section takes the simulation model of an MVDN in Britain as an example to verify the effectiveness of the method [28].The network consists of 15 end nodes and 12 branch points, and the network is divided into 26 sections.To meet the use conditions of the TW location method, the MPs are arranged at the end of each network, as shown in Figure 5, and the sampling rate of the MPs is 10 MHz.The simulation software is PSCAD/EMTDC V4. 6

Technical Performance Verification of Location Method
It is assumed that a phase A ground fault occurs at f1 on section B9B10 in Figure 5, and f1 is 200 m away from B10.The fault occurs at 0.1035 s, the transition resistance is 10 Ω, and the initial phase angle of the fault is 60°.After the grounding fault occurs, the 25 ms transfer device acts, i.e., the faulty phase bus switch closes at 0.1060 s.
The phase voltage travelling wave collected by each measuring point when the fault occurs is subjected to the S transform, and then the final initial travelling wave arrival time of each measuring point is obtained.The phase voltage S transform results are shown in Figure 6, and the corresponding traveling wave arrival time is shown in Table 1.By refer-

Technical Performance Verification of Location Method
It is assumed that a phase A ground fault occurs at f 1 on section B 9 B 10 in Figure 5, and f 1 is 200 m away from B 10 .The fault occurs at 0.1035 s, the transition resistance is 10 Ω, and the initial phase angle of the fault is 60 • .After the grounding fault occurs, the 25 ms transfer device acts, i.e., the faulty phase bus switch closes at 0.1060 s.
The phase voltage travelling wave collected by each measuring point when the fault occurs is subjected to the S transform, and then the final initial travelling wave arrival time of each measuring point is obtained.The phase voltage S transform results are shown in Figure 6, and the corresponding traveling wave arrival time is shown in Table 1.By referring to the above steps, the travelling wave arrival time of each measuring point when the bus is grounded can be obtained, as shown in Table 2.

Technical Performance Verification of Location Method
It is assumed that a phase A ground fault occurs at f1 on section B9B10 in Figure 5, and f1 is 200 m away from B10.The fault occurs at 0.1035 s, the transition resistance is 10 Ω, and the initial phase angle of the fault is 60°.After the grounding fault occurs, the 25 ms transfer device acts, i.e., the faulty phase bus switch closes at 0.1060 s.
The phase voltage travelling wave collected by each measuring point when the fault occurs is subjected to the S transform, and then the final initial travelling wave arrival time of each measuring point is obtained.The phase voltage S transform results are shown in Figure 6, and the corresponding traveling wave arrival time is shown in Table 1.By referring to the above steps, the travelling wave arrival time of each measuring point when the bus is grounded can be obtained, as shown in Table 2.According to Tables 1 and 2, the difference in TWATD between adjacent MPs in the two grounding processes can be obtained and compared with the given time threshold t set .The given t set is mainly related to the TW arrival time error, which can be divided into synchronous time error and TW arrival time detection error.The former mainly depends on the GPS/Beidou system clock synchronization accuracy; usually, the error does not exceed ±1 µs, and the value of the latter is related to the sampling rate of the TW head detection algorithm and the measurement device.The most commonly used detection algorithms are the wavelet transform and S transform.The error range is generally within 2 to 3 sampling moments.Corresponding to the 10 MHz sampling rate used in this paper, the error is within ±0.3 µs.Based on the above analysis, this paper takes t set = 1 µs and then obtains the comparison result of the time of the MP and t set , as shown in Table 3.
Faulty section judgment.Since ∆t 9-12 is greater than t set , it can be seen that the fault is located in one of the sections E 9 E 12 , B 9 E 11 , and B 8 E 10 .Since ∆t 10-11 and ∆t 11-12 are smaller than t set , we can exclude section B 9 E 11 and section B 8 E 10 .Thus, it can be determined that the fault occurs in the corresponding section of E 9 E 12 .Then, using the length of the line in the upstream section of the fault and the arrival time of the TW, the accurate distance between the fault point and the MP M 9 can be obtained.The calculation result is 1406.78m (the actual distance is 1401 m), and the absolute error with the actual position of the fault is 5.78 m.If the relative error is defined as the percentage of the absolute error to the total length of the line (19,495 m), the relative error of the location result is 0.03%.Single-phase ground faults are set at different positions of different lines to verify the effectiveness of the method.The fault setting is as follows: f 2 is located on section B 3 B 4 , 150 m away from section B 3 ; f 3 is located on section B 5 B 6 , 200 m away from section B 6 ; and f 4 is located on section B 8 E 10 , 1000 m away from B 8 .The fault condition is the same as that at f 1 , the transition resistance is 10 Ω, and the initial phase angle of the fault is 60 • .The fault location results are shown in Table 4.It can be seen from Table 4 that the absolute error of fault location is less than 10 m for grounding faults at different positions.It can be seen that this method has good applicability for grounding faults at different positions.The accurate identification of the faulty section is a prerequisite for accurate location.The proposed method can correctly identify the faulty section under different fault locations because the paths between the two measuring points overlap.The faulty section can be identified by reasoning similar to that shown in Table 3.

Location Results under Different Fault Conditions
On this basis, the fault transition resistance and the initial phase angle of the fault are changed to verify the effect of the method under different fault conditions.The simulation results are shown in Table 5.It can be seen from the following table that the error of the location result is the same under the same initial phase angle, indicating that the method is not affected by the transition resistance.If the initial phase angle of the fault is different, the fault transient component is different, and the subsequent TW head identification and the initial TW arrival time detection are also different, resulting in slightly different location results.However, in general, the absolute error of the location results does not exceed 30 m, indicating that the method has certain robustness under various fault conditions.

Performance Analysis of the Proposed Method in Overhead Wire-Cable Hybrid Line Scenario
To verify the effectiveness of the proposed method in the overhead wire-cable hybrid distribution network, the fault at f 1 is still used for analysis, and B 1 B 10 in Figure 6 (i.e., the line between branch point B 1 and branch point B 10 ) is modified into a cable in the simulation model.There have been a large number of studies on the wave velocity translation method of hybrid line faults [29].The wave speed of the overhead line and cable line are calculated by using theoretical line parameters.On this basis, the equivalent overhead network is obtained by converting the length of the cable line with the theory of wave velocity normalization.According to the line parameters, the aerial mode travelling wave speeds of the overhead wire and cable are v overhead = 2.94 × 10 8 m/s and v cable = 1.7052 × 10 8 m/s, respectively.However, considering that the actual line parameters have time-varying characteristics, the calculated wave velocity will be different from the actual situation, resulting in a certain conversion error when the length of the cable line is converted.To demonstrate the effectiveness of the proposed method under different conversion errors, three cases of no error, +2% conversion error, and −5% conversion error are simulated and verified.The results of fault section identification and fault location are shown in Table 6.It can be found that the proposed method can correctly identify the fault section under different degrees of conversion error, and the location error is less than 50 m.Further analysis shows that the faulty section can be correctly identified because the time threshold t set has a certain error tolerance, and its value is 1µs.Without considering the arrival time error of the travelling wave head, the allowable length of cable conversion error corresponding to t set = 1 µs is 170.52 m.In other words, if the line conversion error is 2%, as long as the length of the cable line in the overhead wireline-cable hybrid line does not exceed 8.526 km, the faulty section identification result will not be misjudged.Considering that the travelling wave arrival time may have certain detection errors when the cable exceeds 5 km, it can be considered to appropriately increase the value of time threshold t set or take other measures to deal with the problem of the unequal wave velocity of overhead lines and cables.This will be our next research direction.

Conclusions
In this paper, the TW fault location method of an MVDN under the input of a GFT device is studied, and an accurate fault location method based on the arrival time difference of the TW head is proposed.The main outcomes are as follows.
(1) Taking the idea of faulty section identification and fault point location, the two TW transmission processes of fault grounding and active grounding at the bus are analyzed.It is found that the two TW transmission directions of the upstream MPs of the fault point are opposite, and the time difference of the wave head between the adjacent MPs is different.The two TW transmission directions of the downstream MPs of the fault point are the same, and the time difference of the wave head between the adjacent MPs is the same.Thus, the faulty section can be identified.(2) According to the arrival time and distance of the TW in the upstream section, the location equation is constructed, and an accurate fault location method based on the arrival time difference of the TW head is proposed.Simulation results show that the method has high location accuracy, the absolute error is less than 30 m, and it is not affected by the TW velocity, the fault conditions, or the distributed power sources.(3) The proposed method is still valid for the overhead line and cable hybrid line with a low cable proportion.For the case of a high cable proportion, it can be considered to appropriately increase the value of time threshold t set or take other measures to deal with the problem of the unequal wave velocity of overhead lines and cables.This will be our next research direction.Funding: National Natural Science Foundation of China (51877084).

1
the time at which the initial TW of the fault arrives at the MP Di, D B i m l is the distance between the MP Di and the branch point Bm, v is the speed at which the travelling wave travels through the line, fB m l is the distance between the branch point Bm and the fault point f, and f t is the time at which the fault occurs.By subtracting the adjacent equations in Equation (1) and eliminating the time of the fault, we have
between the branch point Bm+1 and the fault point f.
between the MP Di and the M-terminal MP D 1 .

Figure 5 .
Figure 5. Simulation model of an MVDN in the UK.

Figure 5 .
Figure 5. Simulation model of an MVDN in the UK.

Figure 5 .
Figure 5. Simulation model of an MVDN in the UK.

Table 1 .
The arrival time of line ground TW at each MP.

Table 2 .
The arrival time of bus grounding TW at each MP.

Table 3 .
The TWATD of each adjacent MP.

Table 4 .
The influence of different fault positions on the results of TW fault location.

Table 5 .
The influence of different fault conditions on TW fault location results.

Table 6 .
Location results of the proposed method under different cable line length conversion errors.