Next Article in Journal
Fabrication of a Sensitive and Stable NiO Uric Acid Biosensor Using Ag Nanowires and Reduced Graphene Oxide
Previous Article in Journal
Analysis of Daily Energy Demand for Cooling in Buildings with Different Comfort Categories—Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Faulty Phase Selection Method for Single-Phase-to-Ground Fault in Distribution System Based on Transient Current Similarity Measurement

1
Wenzhou Power Supply Company, State Grid Zhejiang Electric Power Co., Ltd., Wenzhou 325000, China
2
Department of Electrical Engineering, School of Automation, Wuhan University of Technology, Wuhan 430070, China
3
Ji’an Power Supply Company, State Grid Jiangxi Electric Power Co., Ltd., Ji’an 343000, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(15), 4695; https://doi.org/10.3390/en14154695
Submission received: 29 June 2021 / Revised: 25 July 2021 / Accepted: 27 July 2021 / Published: 3 August 2021
(This article belongs to the Topic Power Distribution Systems)

Abstract

:
In modern electrical power distribution systems, the effective operation of inverter-based arc suppression devices relies on the accuracy of faulty phase selection. In the traditional methods of faulty phase selection for single-phase-to-ground faults (SPGs), power frequency-based amplitude and phase characteristics are used to identify the faulty phase. In the field, when a high-resistance SPG occurs in the system, traditional methods are difficult for accurately identifying the faulty phase because of the weak fault components and complicated process. A novel realizable and effective method of faulty phase selection based on transient current similarity measurements is presented when SPGs occur in resonantly grounded distribution systems in this paper. An optimized Hausdorff distance matrix (MOHD) is proposed and constructed by the transient currents of three phases’ similarity measurements within a certain time window of our method. This MOHD is used to select the sampling time window adaptively, which allows the proposed method to be applied to any scale of distribution systems. Firstly, when a SPG occurs, the expressions for the transient phase current mutation in the faulty and sound phases are analyzed. Then, the sampling process is segmented into several selection units (SUs) to form the MOHD-based faulty phase selection method. Additionally, the Hausdorff distance algorithm (HD) is used to calculate the waveform similarities of the transient phase current mutation among the three phases to form the HD-based faulty phase selection method. Finally, a practical resonant grounded distribution system is modeled in PSCAD/EMTDC, and the effectiveness and performance of the proposed method is compared and verified under different fault resistances, fault inception angles, system topologies, sampling time windows and rates of data missing.

1. Introduction

Arc suppression coils have been widely used in power distribution systems in China and Europe [1,2]. When a single-phase-to-ground fault (SPG) occurs in the distribution network, the arc suppression device will compensate the fault current and then prevent the fault arc from reigniting [3,4]. The effective operation of arc suppression devices relies on the accurate selection of faulty phases. Faulty phase selection failure will lead to an in-crease in the ground fault current and not only cannot achieve the fault arcing but also increases the chance of arc reignition [5,6,7]. Seriously, there may be a large-area short circuit of the system due to intermittent overvoltage.
Traditional faulty phase selection methods directly judge the phase with the lowest phase voltage as the faulty phase when a SPG occurs in a resonant grounded distribution network [8]. However, in the case of high-resistance ground faults, the phase with the lowest phase voltage is no longer the faulty phase. To solve this problem, another power frequency-based faulty phase selection method is proposed by some scholars. When the system is in the overcompensation state, the phase leading the highest voltage phase is the faulty phase, and when the system is in the under-compensation state, the phase lagging the highest voltage phase is the faulty phase. Since the above methods can only be used for symmetrical systems, Xu [9] proposed a method to select the faulty phase by constructing a “phase vector”, which used an asymmetric vector to determine the faulty phase, but its calculation process was too complicated. In addition, Meng [10] proposed a faulty phase identification method based on the zero-sequence voltage variation law, which was easy to implement in a digital processor. In reference [11], the Clark Modal Transform (CMT) was used to extract the zero-sequence components from the three-phase voltages, combined with the voltage phase angle shift to determine the faulty phase in the 33-kV distribution system. To avoid the complicated selection process, Wang et al. [12] determined the faulty phase based on the phase voltage variation after the injection of a compensation current at the neutral point in a distribution system with fully compensated arc suppression technology, while Fan et al. [13] proposed a faulty phase recognition method based on phase-to-ground voltage variations, who analyzed the effects of system asymmetry and fault resistance in ungrounded distribution systems. All of the above methods are based on steady-state signals, where the selection of the faulty phase is achieved by the power frequency-based amplitude and phase characteristics of the voltage and current before and after the fault. However, the application of these methods is susceptible to interference from factors such as fault resistance values, system grounding methods and system asymmetries.
The transient process of the power system contains abundant fault signals, and various criteria using fault transient characteristics have been applied in the field of relay protection of the distribution system, such as faulty line selection, fault diagnosis, fault section location, etc. [14]. The relay protection method based on a transient quantity is effective, unaffected by the arc suppression coil, and has a higher sensitivity and better reliability [15]. In reference [16], the gray-scale image was synthesized from the zero-sequence current signal of the system through continuous wavelet transform (CWT), and then, the neural network was used to adaptively extract the fault features to identify the faulty feeder. Compared with reference [16], in addition to using the neural network to detect the faulty feeder, Guo et al. [17] used a feature fusion method, which relieved the dependence of the neural network-based methods on the amount of data. Guo [18] also used a 1D convolutional neural network (1D-CNN) and waveform concatenation to overcome the difficulty of feature extraction and classifier selection. This method can effectively realize the fault location in the case of high-resistance SPG.
In the field, the transient faulty phase selection methods are interfered with less by the grid structure parameters and fault resistance, which have become a hot research topic in power transmission systems [19]. In references [20,21], the affiliation function was constructed from the angle difference and amplitude difference among the basic fault current sequences through WT, and fuzzy logic was used to realize the fault-type diagnosis. Salim et al. [22] extracted the characteristic quantities in a specific frequency band from the transient voltage and current information collected and selected the faulty phase by comparing the threshold values of the faulty phase and the sound phase. Hong et al. [23] proposed an intelligent faulty phase selection method by applying trained deep belief networks (DBNs), which preprocessed the original data to obtain one-dimensional data by min–max standardization and waveform splicing. These faulty phase selections can be effectively realized in the transmission systems. However, they cannot be directly used in a resonant grounded distribution network because of the different sequence impedances, especially the zero-sequence impedance.
Indeed, once a SPG fault occurred on the distribution feeder, there were abundant electromagnetic transient signals because of the energy interactions among the faulty phase and two sound phases. The waveforms of different phase currents may exist as distinguished features. Some papers have used the similarity measurement results of the voltage or current waveforms to achieve a power line protection function [24,25,26,27]. Our research also showed that there is a significant difference between the waveforms of the phase current mutation of the faulty and sound phases. Therefore, it is a viable method to analyze the waveform similarity to select the faulty phase.
In this paper, a faulty phase selection method based on the transient current similarity measurements is proposed for a power distribution system. After a SPG occurs in the distribution network, there is a difference in the waveform of the transient phase current mutation between the faulty phase and the sound phase. The HD algorithm is able to determine the degree of similarity between the waveforms and can be used in cases of high-resistance ground faults. The optimized HD matrix (MOHD) was first proposed in this paper, and it can adaptively determine the sampling time window. The MOHD is constructed by the transient currents similarity measurement in the divided faulty phase selection units (SUs). The effectiveness of the proposed faulty phase selection method is verified under various kinds of fault conditions in PSCAD/EMTDC and MATLAB. The rest of this article is organized as follows. Section 2 describes the reasons why the traditional faulty phase selection methods are not feasible when a high-resistance SPG occurs. Section 3 introduces the basic principles of the proposed faulty phase selection method. Section 4 presents the flow of the SPG faulty phase selection method based on the HD and the MOHD. In Section 5, a typical resonant grounded distribution system is modeled based on PSCAD/EMTDC to verify the accuracy of the proposed faulty phase selection methods. In Section 6, the effectiveness of the proposed method is discussed in different fault resistances, fault inception angles, system topologies and sampling time windows. The conclusions are presented in Section 7.

2. Analysis of Phase Voltage under High-Resistance Ground Fault

When SPGs occur in the noneffectively grounded distribution system, the phase with the lowest phase voltage is identified as the faulty phase. When the fault resistance is large, the steady-state characteristics of the system are weak, and the use of this method can lead to misclassifications. There is a SPG in the resonant grounded distribution system, as shown in Figure 1.
The composite sequence network diagram formed based on the symmetrical component method is shown in Figure 2. X C + , X C and X C 0 , respectively, are the system’s capacitive reactions of the positive-sequence, negative-sequence and zero-sequence. U g + , U g and U g 0 , respectively, are equivalent three-sequence voltage sources at the fault point. I f + , I f and I f 0 are the three-sequence fault currents at the fault point, respectively. Z S + , Z S and Z S 0 , respectively, are the three-sequence equivalent impedances of the distribution system feeders. Z T + , Z T and Z T 0 are the positive-sequence, negative-sequence and zero-sequence impedances of the HV/MV transformer, respectively. E φ + is the positive sequence power electromotive force of the distribution system.
The three-sequence fault current at the fault point in Figure 2 is given by
I f + = I f = I f 0 = 1 3 I f = E A 3 R f + Z + + Z + Z 0
E A is the electromotive force of the power supply. Z + , Z and Z 0 , respectively, are the three-sequence equivalent impedances of the feeder from the bus to the fault point. In a resonant grounded system, there are Z 0 > > Z + and Z 0 > > Z , and the Equation (1) becomes
I f + = I f = I f 0 = 1 3 I f = E A 3 R f + Z 0
where
Z 0 = X C 0 / / j 3 ω L = j 3 ω L 1 3 ω 2 L C
There is a large resistance between the distribution network and ground. The zero-sequence impedance of the HV/MV transformer and the zigzag transformer is very small, and the resistance of the arc suppression coil is also very small. Therefore, the r 0 , the Z T 0 and r L are ignored in Equation (3). According to Equations (2) and (3), the zero-sequence voltage U g 0 can be obtained as follows:
U g 0 = I f 0 Z 0 = E A 1 + j R f ( 3 ω C 1 / ω L )
Therefore, the three-phase voltage U A , U B and U C can be expressed by
U A = E A E A 1 + j R f ( 3 ω C 1 / ω L )
U B = E B E A 1 + j R f ( 3 ω C 1 / ω L )
U C = E C E A 1 + j R f ( 3 ω C 1 / ω L )
A medium voltage distribution system is modeled in PSCAD/EMTDC. The three-phase power supply voltage is 110 kV, connected to a HV/MV transformer with a rated capacity of 50 MVA and a rated ratio of 110/10. Seven feeders are connected to the secondary transformer. The length of all the feeders is 20 km. The arc suppression coil is grounded through the secondary side of the zigzag transformer, and the equivalent inductance value is 0.1078 H. The rated capacity of the distribution transformer is 0.63 MVA. The distribution loads are connected to the end of the feeders with capacity of 0.39 + j0.18 MVA.
The parameters of the feeder are as follows:
r 1 = 0.11   Ω / km ,   L 1 = 0.52   mH / km ,   C 1 = 0.29   μ F / km ; r 0 = 0.34   Ω / km ,   L 0 = 1.54   mH / km ,   C 0 = 0.19   μ F / km .
Assuming a SPG occurs in phase A of the distribution system, the relationship between the three-phase voltage amplitude and the fault resistance is shown in Figure 3.
It can be seen from Figure 3 that, when a SPG occurs in phase A of the distribution system, with the increase of the fault resistance, the voltage amplitude of phase A will no longer be the lowest. When the fault resistance is greater than a certain threshold, the traditional faulty phase selection method will misjudge phase C as the faulty phase. The proposed faulty phase selection method based on transient current similarity measurements in this paper can avoid the interference of the fault resistance and select the faulty phase exactly.

3. Basic Principle of the Proposed Method

The distribution system currents under a SPG are shown in Figure 4. i k p is the phase current near the bus bar on line k. i k p C is the capacitance current to the ground on line k. i k L p is the load current of line k. i f is the fault current at the grounding point (k is the serial number of the feeder; p can be the A, B or C phase).
When there is no fault in the distribution network, the feeder phase currents are expressed by
i k p = i k p C + i k L p = C k d ( e p + u 0 ) d t + i k L p
u 0 is the neutral point voltage before the fault. e p is the effective value of the phase voltage. In Figure 4, the phase currents of the sound feeder and the sound phase of the faulty feeder are expressed by
i k p = i k p C + i k L p = C k d ( e p + u 0 ) d t + i k L p
u 0 is the neutral point voltage after the fault. i k p is the phase current near the bus after the fault. When a SPG occurs in the distribution system, the line voltage of the system remains constant, so the load current i k L p before and after the fault also remains constant, assuming that the source voltages before the fault have sinusoidal variations over time. In Figure 4, the faulty phase current of fault feeder in the distribution system is expressed by
i 2 A = i 2 A C + i 2 L A             = C 2 d ( e A + u 0 ) d t + i 2 L A + ( I c m E p ω L ) cos ( ω t + φ )   + I c m ( ω f ω sin φ sin ω f t cos φ cos ω f t ) e t τ C   + E p ω L cos e t τ L
i 2 A C is the capacitance current to the ground after the fault. I c m is the capacitance current amplitude. φ is the initial phase angle of the fault voltage. ω f is the angular frequency of the free oscillation component. τ C and τ L are the equivalent time constant. ω f , τ C and τ L are related to the inductance value of the arc suppression coil, the leakage impedance to the ground and the capacitance to the ground.
i 2 A is made up of the capacitive current to the ground, the load current, the steady-state component of the fault current and the transient component of the fault current. The steady-state component of the fault current is expressed by ( I c m E p ω L ) cos ( ω t + φ ) , which is the difference between the magnitude of the steady-state capacitive current and the steady-state inductive current. The transient component of the fault current is expressed by I c m ( ω f ω sin φ sin ω f t cos φ cos ω f t ) e t τ C + E p ω L cos e t τ L , which is the sum of the transient free oscillatory component of the capacitive current and the transient DC component of the inductive current.
By subtracting Equation (8) from Equation (10), the phase current mutation of phase A of feeder 2 where a SPG occurs can be expressed by
Δ i 2 A = i 2 A i 2 A = C 2 d ( u 0 + u 0 ) d t + ( I c m E p ω L ) cos ( ω t + φ ) + I c m ( ω f ω sin φ sin ω f t cos φ cos ω f t ) e t τ C + E p ω L cos e t τ L
By subtracting Equation (8) from Equation (9), the phase current mutation of the sound feeder and the sound phase of the faulty feeder can be expressed by
Δ i k p = i k p i k p = C k d ( u 0 u 0 ) d t
In a distribution system with k feeders, when a SPG occurs in phase A, the phase current mutation in phase A on the low-voltage side of the HV/MV transformer can be expressed by
Δ i A = k = 1 n Δ i k A               = k = 1 n C k d ( u 0 + u 0 ) d t + ( I c m E P ω L ) cos ( ω t + φ )               + I c m ( ω f ω sin φ sin ω f t cos φ cos ω f t ) e t τ C               + E p ω L cos e t τ L
The phase current mutation of the sound feeder and the sound phase of the faulty feeder on the low-voltage side of the HV/MV transformer can be expressed by
Δ i B = Δ i C = k = 1 n Δ i k C = k = 1 n C k d ( u 0 u 0 ) d t
According to Equations (13) and (14), the conclusions are as follows:
  • When a SPG occurs in the distribution system, for the sound phase current on the low-voltage side of the HV/MV transformer, the phase current mutation is the sum of the phase-to-ground capacitance current changes. For the faulty phase current on the low voltage side of the HV/MV transformer, the phase current mutation is the sum of the phase-to-ground capacitance current changes and the power–frequency steady-state current components and the power–frequency transient current components at the fault point.
  • When a SPG occurs in the distribution system, one faulty phase current and two sound-phase currents are collected on the low voltage side of the HV/MV transformer. The similarity between the phase current mutation waveforms of the two sound phases is high. There is a significant difference between the phase current mutations of the faulty phase and the sound phase, and the similarity between the phase current mutation waveforms of the faulty and sound phases is low.

4. Faulty Phase Selection Method Based on HD and MOHD

At present, there are two main image similarity algorithms. One is to judge, based on the characteristic points between the images, the typical Hausdorff distance algorithm [25,26]; the other is to judge, based on the gray information between the images, a typical gray correlation analysis [27]. In this section, the HD was used to calculate the similarities among the transient three-phase current mutations when the SPGs occurred. Then, the faulty phase selection criterions suitable for the distribution system were derived, and the SPG faulty phase selection method based on HD was formed.

4.1. HD-Based Faulty Phase Selection Process

The calculation object of the Hausdorff distance algorithm is the discrete dataset of the two images to be compared. Its calculation results are the characteristics’ distances between the two datasets, which can reflect the similarity between the two datasets and the similarity of the two images. The proposed faulty phase selection method based on HD includes the following four parts.

4.1.1. Data Sampling and Normalization

The proposed faulty phase selection method performs similarity calculations and faulty phase selection based on the phase current mutation of a one-quarter cycle.
Using 20 kHz as the sampling frequency, we collected the three-phase current of the low-voltage side of the HV/MV transformer 20~15 ms before the fault and 5 ms after the fault, totaling 600 data. The phase current sequences after the fault minus the ones before the fault achieve the phase current mutation data of the three-phase current, with 100 data per phase. The sampling data sequence of the phase A current mutation is defined as X A ( n ) ; the sampling data sequence of the phase B current mutation is defined as X B ( n ) . The sampling data sequence of the phase C current mutation is defined as X C ( n ) .
Normalizing the sampling data of the current mutation obtains the three-phase current mutation data sequence Y A ( n ) , Y B ( n ) and Y C ( n ) , and the following relations are satisfied:
Y A ( n ) = X A ( n ) X A ( n ) 2 + X B ( n ) 2 + X C ( n ) 2
Y B ( n ) = X B ( n ) X A ( n ) 2 + X B ( n ) 2 + X C ( n ) 2
Y C ( n ) = X C ( n ) X A ( n ) 2 + X B ( n ) 2 + X C ( n ) 2

4.1.2. Obtaining the Similarity Coefficient among the Three Phases

The similarity coefficient of the A–B phase current mutation is P A B , the similarity coefficient of the B–C phase current mutation is P B C and the similarity coefficient of the C–A phase current mutation is P C A .
The calculation process of the similarity coefficient of the A–B phase current mutation is as follows:
• First step
We calculated the absolute value of the difference between the first data Y A ( 1 ) in Y A ( n ) and all the data in Y B ( n ) and selected the minimum value Δ Y A B ( 1 ) among the above absolute values, which is given by
Δ Y A B ( 1 ) = min | Y A ( 1 ) Y B ( n ) |
• Second step
According to the similar method, we calculated Δ Y A B ( 1 ) to Δ Y A B ( 100 ) in turn to form sequence Δ Y A B ( n ) and selected the maximum value in Δ Y A B ( n ) as the one-way distance h ( A , B ) from A to B.
• Third step
We used the same method to calculate the one-way distance h ( B , A ) from B to A and selected the larger of h ( A , B ) and h ( B , A ) as the similarity coefficient P A B of the A–B phase current mutation, which is expressed by
P A B = max [ h ( A , B ) , h ( B , A ) ]
The calculation process of the similarity coefficient P A B of the A–B phase current mutation is completed. The similarity coefficient P B C of the B–C phase current mutation and the similarity coefficient P C A of the C–A phase current mutation can be calculated by the same method.

4.1.3. Obtaining Faulty Phase Selection Parameters

We defined the average similarity of the current mutation sequence of phase A as P A , defined the average similarity of the current mutation sequence of phase B as P B and defined the average similarity of the current mutation sequence of phase C as P C . They can be expressed by
P A = P A B + P C A 2
P B = P A B + P B C 2
P C = P B C + P C A 2
We defined the phase A fault evaluation parameter as λ A , defined the phase B fault evaluation parameter as λ B and defined the phase C fault evaluation parameter as λ C . They can be expressed by
λ A = | P A P B | + | P A P C |
λ B = | P B P A | + | P B P C |
λ C = | P C P A | + | P C P B |
We set the faulty phase selection standard λ s e t and the sensitivity factor ε ( λ A , B , C / λ s e t ) to indicate the practicality of the proposed method. λ s e t is expressed by
λ s e t = 1 3 ( λ A + λ B + λ C )

4.1.4. Faulty Phase Selection Process

When a SPG occurs in a resonant grounded distribution system, the process of the proposed faulty phase selection method based on HD in this paper is as follows.
Firstly, we performed data sampling and normalization. Secondly, we calculated the phase-to-phase similarity coefficient. Then, we calculated the faulty phase selection parameters. Finally, we compared the A, B and C phase faulty evaluation parameters and the faulty phase selection standard, respectively.
If λ A > λ s e t , it is judged that there is a SPG occurring in phase A; if λ B > λ s e t , it is judged that there is a SPG occurring in phase B and if λ C > λ s e t , it is judged that there is a SPG occurring in phase C. The faulty phase selection flowchart is shown in Figure 5.

4.2. MOHD-Based Faulty Phase Selection Process

In the field, there is a short transient process when a SPG occurs in the distribution system. The duration of the transients also varies for different scales of the system. The moment of fault occurrence is determined by detecting a sudden rise in zero-sequence voltage at the substation. The distribution system operates in a harsh environment, and the signals are easily disturbed by a variety of factors. It may happen that the collection terminal online is not sampled for a very short time after the zero-sequence voltage fluctuates, resulting in the loss of transient fault data. The HD-based faulty phase selection method proposed in Section 4.1 uses a fixed 5-ms time window. In the field, the faulty phase may be incorrectly selected.
The faulty phase selection matrix (MOHD) based on an optimized HD algorithm (OHD) was first proposed in this paper. With the MOHD, the faulty phase selection time window can be determined adaptively. Suppose a SPG occurs in the distribution system, the transient data within T after the fault is missing, and the transient process ends at T’ after the fault. It therefore makes sense for transient fault data to be sampled from T to T’ after the fault. This valid data sampling time window is defined as W T T .
The sampling process is divided into n faulty phase selection units (SUs), each with a time duration of ms. For each SU, a data sequence of the three-phase current mutation will be synthesized, and then, the faulty phase selection will be completed. The above process will be repeated n times, and there will be a continuous data sampling process. The transient data sampling time window is shown in Figure 6.
In Figure 6, a SPG occurs at 0.2 s, with a data loss of 0.003 s after the fault. At T’ seconds after the fault, the fault transient process ends. The transient data sampling time window W T T is made up of n SUs.
The faulty phase selection vector α n is set to represent the faulty phase selection result of the nth SU.
When the faulty phase is phase A, α n = [ 1 , 0 , 0 ] T ;
When the faulty phase is phase B, α n = [ 0 , 1 , 0 ] T ;
When the faulty phase is phase C, α n = [ 0 , 0 , 1 ] T .
The matrix MOHD is expressed by
M OHD = [ α 2 α 1 , α 3 α 2 , , α n α n 1 ] 3 × n 1                       = [ β 1 , β 2 , , β n 1 ] 3 × n 1
The time window 1-norm μ was set to indicate the consistency of the α n and the α n 1 , which are expressed by
μ = β n 1 1 = i = 1 3 | β n 1 |
When μ = 0 , this means that the nth SU is still in W T T , and faulty phase selection continues in the n + 1th SU. When μ 0 , this means that the nth SU is already outside W T T , and the faulty phase selection process is complete. The flowchart of faulty phase selection based on MOHD is shown in Figure 7.
The time window WT−T’ is determined as from T to T + a(n − 1) after the fault. With the MOHD and WT−T’, the proposed faulty phase selection method can be adapted to different scales of the distribution system. The effect of the missing sampling data on the faulty phase selection results is also taken into account in this method. In addition, the WT−T’ can provide a reliable transient database for the rest of the transient-based relay protection methods.

5. Simulation Verification and Results Analysis

In this paper, a typical medium voltage distribution system is modeled in a Power Systems Computer-Aided Design/Electromagnetic Transients, including the DC. The three-phase current mutation is collected in PSCAD/EMTDC, and the similarity coefficient is calculated in MATLAB. The distribution system model is shown in Figure 8.
The parameters of the HV/MV transformer and the 10-kV feeder are shown in Table 1 and Table 2, respectively. The arc suppression coil is grounded through the secondary side of the zigzag transformer, and its equivalent inductance value is 0.1078 H. The distribution loads are connected to the end of the feeders with capacities of 0.39 + j0.18 MVA.
Suppose a SPG occurs in phase A of feeder 2 of the distribution system simulation model with a system sampling rate of 20 kHz. The fault inception angle of the voltage is 90°, and the fault resistance is 1 Ω, 25 Ω, 100 Ω, 500 Ω, 1000 Ω and 2000 Ω, respectively. After the SPG occurs, the waveforms of the phase current mutation in phases A, B and C are shown in Figure 9.
As can be seen from Figure 9, the phase current mutations in phase A, phase B and phase C all decrease gradually as the fault resistance increases. During this process, the differences between the phase current mutation waveforms of the faulty phase and the ones of the sound phase gradually decrease, and the difference between the phase current mutation waveforms of the two sound phases gradually increases. When the fault resistance is 1000 Ω and above, there is no situation where the difference between the faulty phase and the sound phase is more significant than the difference between the two sound phases, which increases the difficulty of faulty phase selection.
We calculated the similarity coefficient of the data sequence of the A, B and C phase current mutations through the Hausdorff distance algorithm in MATLAB. Then, we calculated the fault evaluation parameter and faulty phase selection standard. The faulty phase selection was completed, and the results are shown in Table 3.
Table 3 shows that when the fault resistance varies from 1 Ω to 2000 Ω, the fault evaluation parameters λA of the faulty phase A are all much greater than the faulty phase selection standard λset. When the fault resistance is a constant value, the proposed faulty phase selection method based on the transient current similarity measurements correctly identifies the faulty phase A. In all cases, the sensitivity factors ε (λA/λset) exceeded 1.4, so the proposed method is effective in harsh environments.
Three sets of fault signals were collected in a real power distribution system in the Wenzhou area of China. The fault signals were collected in a fault indicator. The fault indicator was installed near the output of the substation with a sampling frequency of 4096 Hz. The performance of the proposed method was verified, and the results are shown in Table 4. The faulty phase selection method based on the transient current similarity measurements proposed in this paper correctly selected the faulty phase in a real distribution system.
The three methods of faulty phase selection that are widely used in the field are listed as follows. The method proposed in this paper (Method 1) is compared with these three methods.
  • Method 1: The novel faulty phase selection method for SPG in a distribution system based on the transient current similarity measurements proposed in this paper.
  • Method 2: The phase with the lowest phase voltage is judged to be the faulty phase when a SPG occurs in a resonant grounded distribution system.
  • Method 3: When the system is in the overcompensation state, the phase leading to the highest voltage phase is the faulty phase, when the system is in the under-compensation state, the phase lagging the highest voltage phase is the faulty phase.
  • Method 4: The first half-wave method.
The first half-wave method has been widely used in faulty line selection, which was used to select the fault phase in this paper. The polarity standard is defined as ζ . When the polarity of the phase current is positive, ζ > 0 . When the polarity of the phase current is negative, ζ < 0 . ζ can be expressed by
ζ = { k = 1 M s i g n ( i ( k ) ) , s i g n ( i ( M ) · i ( M + 1 ) ) 0 , M [ 1 , N ) k = 1 N s i g n ( i ( k ) ) , s i g n ( i ( m ) · i ( m + 1 ) ) = 1 , m [ 1 , N )
i ( k ) is the kth sample value of the fault current. N is the number of half-cycle sampling points. In the simulation model in Figure 8, the fault resistance was set as 100 Ω and 2000 Ω, respectively. The distribution system was set up as a symmetrical system and an asymmetrical system, respectively, and the asymmetry of the system was set as 0% and −1.095%, respectively. The results of the faulty phase selection for Method 1 and the other methods are shown in Table 5. The method proposed in this paper (Method 1) is able to select the faulty phase accurately, without interference from fault resistance, system compensation and system symmetry. When the fault resistance is high, Methods 2 and 4 will identify the sound phase as the faulty phase by mistake. Method 3 will identify the sound phase as the faulty phase in an asymmetrical distribution network.

6. Discussion of the Anti-Interference

6.1. The Influence of the Sudden Change of the Fault Resistance Value on the Proposed Method

When SPGs occur in the distribution system, there may be an arc at the fault point, so the fault resistance is not constant. In a resonant grounded distribution system, the fault point is set at phase A of feeder 2. There is a sudden change in the fault resistance value 2.5 ms after the fault. The change of the fault resistance is shown below:
  • The fault resistance changes abruptly from 25 to 50 and 1000, respectively.
  • The fault resistance changes abruptly from 1000 to 800 and 25, respectively.
The similarity parameters and faulty phase selection results under the influence of the sudden change of the fault resistance are shown in Table 6. The results show that the proposed method is able to accurately select the faulty phase A when there is a sudden increase or a sudden decrease in the fault resistance.

6.2. The Influence of the Fault Inception Angle on the Proposed Method

The fault inception angle of the phase voltage affects the waveform of the phase current mutation. The performance of the proposed method may be disturbed by the fault inception angle. In the simulation model in Section 5, the fault resistance is set to 100 Ω, and the fault phase is phase A of feeder 2. The fault inception angle is set to 0°, 30°, 60° and 90°, respectively. The similarity parameters among the waveforms of the three-phase current mutations are shown in Table 7. The results show that the faulty phase selection method based on the transient current similarity measurements is not affected by the fault inception angle.

6.3. The Influence of the System Topology on the Proposed Method

The distribution system is directly connected to a large number of electricity consumers. There are a variety of network topologies in the power distribution system. In this paper, the effect of system topology on the proposed method is explored by varying the number of feeders. The model in Section 5 will be designed as four new scenarios that contain one feeder, three feeders, five feeders and seven feeders, respectively. The fault resistance is set to 100Ω, and the fault inception angle is set to 90°. The phase similarity coefficients are calculated by the HD, and the faulty phase selection results are listed in Table 8. The results show that the HD-based faulty phase selection method is not affected by changes in the number of feeders. This method has a certain degree of anti-disturbance when the system topology changes.

6.4. The Influence of the Sampling Time Window on the Proposed Method

The sampling data length will affect the similarity coefficient of the three-phase current mutations calculated by the HD. When the sampling time window is increased, there may be a loss of accuracy because of the inclusion of useless data from the fault steady state. At the same time, the efficiency of the faulty phase selection may be reduced because of the larger size of the sampling database. When the sampling data length is small, the calculation speed is improved, but a smaller amount of sampled data may lead to incorrect faulty phase selection. When the fault resistance is large or the fault inception angle is small, the transient decay process is also very short. We modified the sampling time length to 1 ms, 2 ms, 5 ms and 10 ms, respectively. The simulation results are shown in Table 9. As shown in the table, the method based on HD can correctly select the faulty phase when the sampling time window fluctuates within a limited range.

6.5. The Influence of the Missing Sampling Data on the Proposed Method

There is a complex operating environment for power distribution systems. Missing sampling data often occurs in rural distribution networks. Assuming that a small amount of data is missing in a short period of time after the fault, the data missing rate is set at 10%, 20%, 30% and 40%, respectively. The faulty phase current is collected in the simulation model in Section 5. The sensitivity factor ε (λA/λset) is set to 1.4. The data missing rate and faulty phase selection results are shown in Table 10. When the sampling data missing rate is less than 30%, the method proposed in this paper still accurately identifies the faulty phase.

7. Conclusions

When a high-resistance SPG occurs in a distribution system, the steady-state components of the fault signals are too weak to be used for faulty phase selection. The traditional faulty phase selection methods based on the power–frequency amplitude and phase characteristics have a complicated and unfeasible process. A novel faulty phase selection method based on the transient current similarity measurements is proposed in this paper.
The data sampling process has been divided into several faulty phase selection units (SUs). For each SU, the transient phase current mutation is synthesized from the difference in the phase currents before and after the fault. The HD algorithm is used to normalize the phase current mutations and calculate the waveform similarities of the transient phase current mutations among the three phases of the distribution system. The MOHD is synthesized based on the faulty phase selection results of the SUs. The phases where the fault evaluation parameter (λA,B,C) exceeds the selection standard (λset) are identified as faulty phases. With the MOHD, the proposed faulty phase selection method can be applied to any scale of the distribution systems, and the proposed method can be applied in cases where a small amount of sampling data is missing. After a theoretical analysis and simulation verification, the following two conclusions can be drawn:
  • When any SPG occurs in a typical medium-voltage distribution system, the proposed faulty phase selection method can accurately select the faulty phase over a wide range of fault resistances with a system sampling rate of 20 kHz and a time window of 5 ms. This method is still valid, even when the fault resistance reaches 2000 Ω or the fault resistance varies during the fault process. In all cases, the sensitivity factors ε (λA,B,C/λset), the ratio of the evaluation parameter of the faulty phase to the selection standard, all exceeded 1.4.
  • When the fault inception angle varies from 0° to 90°or the number of system feeders changes, the proposed method can select the faulty phase correctly. Although the similarity between the faulty phase and sound phase increases as the sampling time, the window changes from 1 ms to 10 ms, which raises the difficulty of faulty phase selection, and the proposed method can still accurately select the faulty phase. In addition, the proposed method is still valid when the missing rate of the sampled data is below 30%.

Author Contributions

Conceptualization, J.T. and Y.T.; methodology, J.T.; validation, Y.T., Y.C. and H.Y.; formal analysis, H.Y. and B.X.; writing—original draft preparation, Y.C. and H.Y.; writing—review and editing, J.T. and project administration, M.Y. and Y.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51707139, and by Wenzhou Power Supply Company of State Grid Zhejiang Electric Power Co., Ltd. of China, grant number 5211WZ1900RS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, B.; Yao, G.; Yu, K.; Zeng, X.; Wang, W.; Zhuo, C.; Guerrero, J.M. Principle of Flexible Ground-Fault Arc Suppression Device Based on Zero-Sequence Voltage Regulation. IEEE Access 2021, 9, 2382–2389. [Google Scholar] [CrossRef]
  2. Toader, D.; Greconici, M.; Vesa, D.; Vintan, M.; Solea, C. Analysis of the Influence of the Insulation Parameters of Medium Voltage Electrical Networks and of the Petersen Coil on the Single-Phase-to-Ground Fault Current. Energies 2021, 14, 1330. [Google Scholar] [CrossRef]
  3. Wang, W.; Yan, L.; Zeng, X.; Fan, B.; Guerrero, J.M. Principle and Design of a Single-Phase Inverter-Based Grounding System for Neutral-to-Ground Voltage Compensation in Distribution Networks. IEEE Trans. Ind. Electron. 2017, 64, 1204–1213. [Google Scholar] [CrossRef] [Green Version]
  4. Peng, W.; Baichao, C.; Cuihua, T.; Sun, B.; Zhou, M.; Yuan, J. A Novel Neutral Electromagnetic Hybrid Flexible Grounding Method in Distribution Networks. IEEE Trans. Power Deliv. 2017, 32, 1350–1358. [Google Scholar] [CrossRef]
  5. Zeng, X.; Wang, Y.; Li, J.; Xiong, T. Novel Principle of Faults Arc Extinguishing & Feeder Protection Based on Flexible Grounding Control for Distribution Networks. Proc. CSEE 2012, 32, 137–143. [Google Scholar]
  6. Guo, M.; You, J.; Zhang, W.; Yang, G.; Miao, X. Separate-Phase Flexible Arc-Suppression Method of Earth-Fault in Distribution Systems Based on Three-Phase Cascaded H-Bridge Converter. Trans. China Electrotech. Soc. 2016, 31, 11–22. [Google Scholar]
  7. Burgess, R.; Ahfock, A. Minimising the Risk of Cross-Country Faults in Systems Using Arc Suppression Coils. IET Gener. Transm. Distrib. 2011, 5, 703–711. [Google Scholar] [CrossRef] [Green Version]
  8. Zeng, X.; Zhuo, C.; Yu, K.; Xiang, G.; Li, J.; Zou, H. A Novel Method of Faults Arc Extinguishing and Feeder Protection Based on Voltage Regulating Intervention With Grounding Transformer Winding Taps for Distribution Networks. Proc. CSEE 2020, 40, 1523–1533. [Google Scholar]
  9. Xu, B.; Zhang, J.; Cai, X.; Guan, S.; Yao, Y. Grounding Phase Determination in Non-Effective Grounding Systems in Case of Unsymmetrical Voltage. Trans. China Electrotech. Soc. 2011, 26, 175–182. [Google Scholar]
  10. Meng, J.; Wang, W.; Tang, X.; Xu, X. Zero-sequence voltage trajectory analysis for unbalanced distribution networks on single-line-to-ground fault condition. Electr. Power Syst. Res. 2018, 161, 17–25. [Google Scholar] [CrossRef]
  11. Alanzi, E.A.; Younis, M.A.; Ariffin, A.M. Detection of faulted phase type in distribution systems based on one end voltage measurement. Int. J. Electr. Power 2014, 54, 288–292. [Google Scholar] [CrossRef]
  12. Wang, J.; Tang, J.; Hou, H.; Liu, B.; Yang, C. A New Phase Selection Method for Single-Phase Grounding Faults in Distribution Networks with Full Compensation Arc Suppression Technology. In Proceedings of the 2018 China International Conference on Electricity Distribution (CICED), Tianjin, China, 17–19 September 2018. [Google Scholar]
  13. Fan, B.; Yao, G.; Wang, W.; Yang, X.; Ma, H.; Yu, K.; Zhuo, C.; Zeng, X. Faulty phase recognition method based on phase-to-ground voltages variation for neutral ungrounded distribution networks. Electr. Power Syst. Res. 2021, 190, 106848. [Google Scholar] [CrossRef]
  14. Haleem, A.M.I.; Sharma, M.; Sajan, K.S.; Babu, K.N.D. A Comparative Review of Fault Location/Identification Methods in Distribution Networks. In Proceedings of the 2018 1st International Conference on Advanced Research in Engineering Sciences (ARES), Dubai, United Arab Emirates, 15–15 June 2018. [Google Scholar]
  15. Sapountzoglou, N.; Lago, J.; De Schutter, B.; Raison, B. A generalizable and sensor-independent deep learning method for fault detection and location in low-voltage distribution grids. Appl. Energy 2020, 276, 115299. [Google Scholar] [CrossRef]
  16. Guo, M.; Zeng, X.; Chen, D.; Yang, N. Deep-Learning-Based Earth Fault Detection Using Continuous Wavelet Transform and Convolutional Neural Network in Resonant Grounding Distribution Systems. IEEE Sens. J. 2018, 18, 1291–1300. [Google Scholar] [CrossRef]
  17. Du, Y.; Liu, Y.; Shao, Q.; Luo, L.; Dai, J.; Sheng, G.; Jiang, X. Single Line-to-Ground Faulted Line Detection of Distribution Systems With Resonant Grounding Based on Feature Fusion Framework. IEEE Trans. Power Deliv. 2019, 34, 1766–1775. [Google Scholar] [CrossRef]
  18. Guo, M.; Gao, J.; Shao, X.; Chen, D. Location of Single-Line-to-Ground Fault Using 1-D Convolutional Neural Network and Waveform Concatenation in Resonant Grounding Distribution Systems. IEEE Trans. Instrum. Meas. 2021, 70, 1–9. [Google Scholar] [CrossRef]
  19. Guo, M.; Yang, N.; Chen, W. Deep-Learning-Based Fault Classification Using Hilbert–Huang Transform and Convolutional Neural Network in Power Distribution Systems. IEEE Sens. J. 2019, 19, 6905–6913. [Google Scholar] [CrossRef]
  20. Zhang, J.; He, Z.Y.; Lin, S.; Zhang, Y.B.; Qian, Q.Q. An ANFIS-based fault classification approach in power distribution system. Int. J. Electr. Power 2013, 49, 243–252. [Google Scholar] [CrossRef]
  21. Das, B. Fuzzy Logic-Based Fault-Type Identification in Unbalanced Radial Power Distribution System. IEEE Trans. Power Deliv. 2006, 21, 278–285. [Google Scholar] [CrossRef]
  22. Salim, R.H.; de Oliveira, K.; Filomena, A.D.; Resener, M.; Bretas, A.S. Hybrid Fault Diagnosis Scheme Implementation for Power Distribution Systems Automation. IEEE Trans. Power Deliv. 2008, 23, 1846–1856. [Google Scholar] [CrossRef]
  23. Hong, C.; Zeng, Z.Y.; Fu, Y.Z.; Guo, M.F. Deep-Belief-Networks Based Fault Classification in Power Distribution Networks. IEEJ Trans. Electr. Electron. Eng. 2020, 15, 1428–1435. [Google Scholar] [CrossRef]
  24. Song, G.; Wang, C.; Wang, T.; Kheshti, M.; Kang, X. A Phase Selection Method for Wind Power Integration System Using Phase Voltage Waveform Correlation. IEEE Trans. Power Deliv. 2017, 32, 740–748. [Google Scholar] [CrossRef]
  25. Zhang, D.; Zou, L.; Chen, Y.; He, F. Efficient and Accurate Hausdorff Distance Computation Based on Diffusion Search. IEEE Access 2018, 6, 1350–1361. [Google Scholar] [CrossRef]
  26. Hua, T.; Dai, K.; Zhang, X.; Yao, Z.; Wang, H.; Xie, K.; Feng, T.; Zhang, H. Optimal VMD-Based Signal Denoising for Laser Radar via Hausdorff Distance and Wavelet Transform. IEEE Access 2019, 7, 167997–168010. [Google Scholar] [CrossRef]
  27. Wang, Y.; Huang, Y.; Zeng, X.; Wei, G.; Zhou, J.; Fang, T.; Chen, H. Faulty Feeder Detection of Single Phase-Earth Fault Using Grey Relation Degree in Resonant Grounding System. IEEE Trans. Power Deliv. 2017, 32, 55–61. [Google Scholar] [CrossRef]
Figure 1. A phase-to-ground fault in the resonant grounded distribution system, where T is the HV/MV transformer, and R f is the fault resistance. C is the equivalent capacitance of the system. r is the equivalent resistance of the system. L is the equivalent inductance of the arc suppression coil. r L is the equivalent inductance of the arc suppression coil.
Figure 1. A phase-to-ground fault in the resonant grounded distribution system, where T is the HV/MV transformer, and R f is the fault resistance. C is the equivalent capacitance of the system. r is the equivalent resistance of the system. L is the equivalent inductance of the arc suppression coil. r L is the equivalent inductance of the arc suppression coil.
Energies 14 04695 g001
Figure 2. Composite sequence network of a phase-to-ground fault. X C + , X C and X C 0 , respectively, are the system’s capacitive reactions of the positive-sequence, negative-sequence and zero-sequence. U g + , U g and U g 0 , respectively, are equivalent three-sequence voltage sources at the fault point. I f + , I f and I f 0 are the three-sequence fault currents at the fault point, respectively. Z S + , Z S and Z S 0 , respectively, are the three-sequence equivalent impedances of the distribution system feeders. Z T + , Z T and Z T 0 are the positive-sequence, negative-sequence and zero-sequence impedances of the HV/MV transformer, respectively. E φ + is the positive sequence power electromotive force of the distribution system.
Figure 2. Composite sequence network of a phase-to-ground fault. X C + , X C and X C 0 , respectively, are the system’s capacitive reactions of the positive-sequence, negative-sequence and zero-sequence. U g + , U g and U g 0 , respectively, are equivalent three-sequence voltage sources at the fault point. I f + , I f and I f 0 are the three-sequence fault currents at the fault point, respectively. Z S + , Z S and Z S 0 , respectively, are the three-sequence equivalent impedances of the distribution system feeders. Z T + , Z T and Z T 0 are the positive-sequence, negative-sequence and zero-sequence impedances of the HV/MV transformer, respectively. E φ + is the positive sequence power electromotive force of the distribution system.
Energies 14 04695 g002
Figure 3. Changes in the voltage amplitude of three phases: A, B and C with fault resistance.
Figure 3. Changes in the voltage amplitude of three phases: A, B and C with fault resistance.
Energies 14 04695 g003
Figure 4. Currents of a phase-to-ground fault distribution system. T is the HV/MV transformer. L is the equivalent inductance of the arc suppression coil. r L is the equivalent inductance of the arc suppression coil. u 0 is the neutral point voltage before the fault. i f is the fault current at the grounding point. e A , e B and e C , respectively, are the effective values of phase A, phase B and phase C. i 1 A , i 1 B and i 1 C , respectively, are the phase currents of phase A, phase B and phase C near the bus bar on line 1. i 1 L A , i 1 L B and i 1 L C , respectively, are the load currents of phase A, phase B and phase C on line 1. i 1 A C , i 1 B C and i 1 C C , respectively, are the capacitance currents to the ground of phase A, phase B and phase C on line 1. i 2 A , i 2 B and i 2 C are the phase currents on line 2. i 2 L A , i 2 L B and i 2 L C are the load currents on line 2. i 1 A C , i 1 B C and i 1 C C are the capacitance currents to the ground on line 2.
Figure 4. Currents of a phase-to-ground fault distribution system. T is the HV/MV transformer. L is the equivalent inductance of the arc suppression coil. r L is the equivalent inductance of the arc suppression coil. u 0 is the neutral point voltage before the fault. i f is the fault current at the grounding point. e A , e B and e C , respectively, are the effective values of phase A, phase B and phase C. i 1 A , i 1 B and i 1 C , respectively, are the phase currents of phase A, phase B and phase C near the bus bar on line 1. i 1 L A , i 1 L B and i 1 L C , respectively, are the load currents of phase A, phase B and phase C on line 1. i 1 A C , i 1 B C and i 1 C C , respectively, are the capacitance currents to the ground of phase A, phase B and phase C on line 1. i 2 A , i 2 B and i 2 C are the phase currents on line 2. i 2 L A , i 2 L B and i 2 L C are the load currents on line 2. i 1 A C , i 1 B C and i 1 C C are the capacitance currents to the ground on line 2.
Energies 14 04695 g004
Figure 5. Flowchart of the faulty phase selection method based on HD.
Figure 5. Flowchart of the faulty phase selection method based on HD.
Energies 14 04695 g005
Figure 6. The transient data sampling time window.
Figure 6. The transient data sampling time window.
Energies 14 04695 g006
Figure 7. Flowchart of faulty phase selection method based on MOHD.
Figure 7. Flowchart of faulty phase selection method based on MOHD.
Energies 14 04695 g007
Figure 8. Simulation model of the resonant grounded distribution.
Figure 8. Simulation model of the resonant grounded distribution.
Energies 14 04695 g008
Figure 9. The waveform diagram of the phase current mutation when the inception-phase angle of the fault voltage is 90°. (a) The phase current mutation waveform when the fault resistance is 1 Ω. (b) The phase current mutation waveform when the fault resistance is 25 Ω. (c) The phase current mutation waveform when the fault resistance is 100 Ω. (d) The phase current mutation waveform when the fault resistance is 500 Ω. (e) The phase current mutation waveform when the fault resistance is 1000 Ω. (f) The phase current mutation waveform when the fault resistance is 2000 Ω.
Figure 9. The waveform diagram of the phase current mutation when the inception-phase angle of the fault voltage is 90°. (a) The phase current mutation waveform when the fault resistance is 1 Ω. (b) The phase current mutation waveform when the fault resistance is 25 Ω. (c) The phase current mutation waveform when the fault resistance is 100 Ω. (d) The phase current mutation waveform when the fault resistance is 500 Ω. (e) The phase current mutation waveform when the fault resistance is 1000 Ω. (f) The phase current mutation waveform when the fault resistance is 2000 Ω.
Energies 14 04695 g009
Table 1. Parameters of the HV/MV transformer.
Table 1. Parameters of the HV/MV transformer.
Rated Capacity
(MVA)
Voltage Rate (kV)No Load Loss, P0 (kW)Copper Loss, Pk (kW)I0 (%)Uk (%)
50.0110/1025.6125.80.210.5
Table 2. Parameters of the feeders (per kilometer).
Table 2. Parameters of the feeders (per kilometer).
Feeder TypesPhase Sequencer
(Ω/km)
L
(mH/km)
C
(μF/km)
Cable feederPositive sequence0.110.520.25
Zero sequence0.341.540.19
Overhead feederPositive sequence0.0961.220.011
Zero sequence0.233.660.007
Table 3. Similarity parameters and faulty phase selection results.
Table 3. Similarity parameters and faulty phase selection results.
Fault Resistance
(Ω)
Similarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
10.530.010.520.530.270.270.520.260.270.49A
251.000.041.001.000.520.520.960.480.480.91A
1001.000.081.001.000.540.540.920.460.460.85A
5001.030.231.031.020.630.620.790.400.420.76A
10001.050.421.051.040.740.730.620.320.330.59A
20001.100.781.101.100.940.940.320.160.160.29A
Table 4. Similarity parameters and faulty phase selection results in a real power distribution system.
Table 4. Similarity parameters and faulty phase selection results in a real power distribution system.
NumberSimilarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
10.450.950.890.670.700.920.280.250.470.46C
21.071.070.220.651.070.650430.850.430.79B
32.000.451.861.931.231.161.480.780.851.45A
Table 5. The faulty phase selection results of Method 1, Method 2, Method 3 and Method 4.
Table 5. The faulty phase selection results of Method 1, Method 2, Method 3 and Method 4.
SymmetryFault ResistanceMethod 1Method 2Method 3Method 4
symmetrical systemhigh fault resistanceYesNoYesNo
low fault resistanceYesYesYesYes
asymmetric systemhigh fault resistanceYesNoNoNo
low fault resistanceYesYesNoNo
Table 6. Similarity parameters and faulty phase selection results under the influence of the sudden change of the fault resistance value.
Table 6. Similarity parameters and faulty phase selection results under the influence of the sudden change of the fault resistance value.
Fault Resistance (Ω)Similarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
25→501.000.041.001.000.520.520.960.480.480.90A
25→10001.000.041.001.000.520.520.960.480.480.90A
1000→8001.040.421.041.040.730.730.620.310.310.57A
1000→251.010.051.011.010.530.530.960.480.480.90A
Table 7. Similarity parameters and faulty phase selection results under the influence of the fault inception angle.
Table 7. Similarity parameters and faulty phase selection results under the influence of the fault inception angle.
Inception Angle (°)Similarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
1.070.121.011.040.600.570.920.480.500.90A
30°1.010.081.011.010.550.550.930.470.470.87A
60°1.010.071.011.010.540.540.940.470.470.88A
90°1.000.081.001.000.540.540.920.460.460.85A
Table 8. Similarity parameters and faulty phase selection results under the influence of the system topology.
Table 8. Similarity parameters and faulty phase selection results under the influence of the system topology.
Number of FeedersSimilarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
11.000.011.001.000.510.510.990.500.500.92A
31.000.031.001.000.520.520.970.490.490.91A
51.000.051.001.000.530.530.950.480.480.90A
71.000.061.001.000.530.530.940.470.470.88A
Table 9. Similarity parameters and faulty phase selection results under the influence of the sampling time window.
Table 9. Similarity parameters and faulty phase selection results under the influence of the sampling time window.
Time Window
(ms)
Similarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
11.030.051.011.020.540.530.970.490.500.91A
21.030.051.011.020.540.530.970.490.500.91A
51.010.081.011.010.530.530.960.480.480.90A
100.510.090.510.510.300.300.420.210.210.39A
Table 10. Similarity parameters and faulty phase selection results under the influence of the missing sampling data.
Table 10. Similarity parameters and faulty phase selection results under the influence of the missing sampling data.
Missing Rate
(%)
Similarity ParametersResult
PABPBCPCAPAPBPCλAλBλC1.4λset
101.000.081.001.000.540.540.920.460.460.86A
201.000.111.001.000.560.560.890.450.450.83A
301.000.111.001.000.560.560.890.450.450.83A
400.770.130.640.710.450.390.580.320.390.60No
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, Y.; Chang, Y.; Tang, J.; Xu, B.; Ye, M.; Yang, H. A Novel Faulty Phase Selection Method for Single-Phase-to-Ground Fault in Distribution System Based on Transient Current Similarity Measurement. Energies 2021, 14, 4695. https://doi.org/10.3390/en14154695

AMA Style

Tang Y, Chang Y, Tang J, Xu B, Ye M, Yang H. A Novel Faulty Phase Selection Method for Single-Phase-to-Ground Fault in Distribution System Based on Transient Current Similarity Measurement. Energies. 2021; 14(15):4695. https://doi.org/10.3390/en14154695

Chicago/Turabian Style

Tang, Yaojing, Yongle Chang, Jinrui Tang, Bin Xu, Mingkang Ye, and Hongbo Yang. 2021. "A Novel Faulty Phase Selection Method for Single-Phase-to-Ground Fault in Distribution System Based on Transient Current Similarity Measurement" Energies 14, no. 15: 4695. https://doi.org/10.3390/en14154695

APA Style

Tang, Y., Chang, Y., Tang, J., Xu, B., Ye, M., & Yang, H. (2021). A Novel Faulty Phase Selection Method for Single-Phase-to-Ground Fault in Distribution System Based on Transient Current Similarity Measurement. Energies, 14(15), 4695. https://doi.org/10.3390/en14154695

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop