1. Introduction
For economic operation considerations, the non-effectively grounded neutral points have been adopted in the medium-voltage (MV) distribution network design in China [
1]. With widespread application of such system in the grid, fault line identification and fault point location have great significance for maintaining power system reliability and ensuring network equipment safety [
2,
3]. Single-phase-to-ground fault is the most common fault in distribution system. When it occurs, the grounding current of distribution system goes beyond an upper limit [
4], overvoltage is usually generated which endangers the insulation of distribution network and equipment [
5]. More seriously, when the system operates in such a fault state, shock accident or casualties might happen [
6,
7], which represents a serious public safety hazard. Specifically, when a single-phase-to-ground fault occurs in non-effectively grounded system, it is difficult to be detected due to unobvious characteristics of fault current [
8].
In last three decades, scholars have proposed various fault detection methods, which can be summarized as steady-state method, transient method, and signal injection method according to the different electrical signals used. Steady-state method [
9,
10,
11,
12,
13,
14,
15,
16] includes zero-sequence admittance method [
9,
10], zero-sequence current amplitude and phase comparison method [
11,
12,
13], zero-sequence active component direction method [
14], power angle method [
15], negative-sequence current method [
16], etc. With stable signal collection, the steady-state method has been widely used. Transient method [
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29] includes traditional energy method [
17], centralized comparison method [
18], wavelet analysis method [
19], correlation analysis method [
20], etc. It uses transient signals thus obtains abundant fault transient information and features [
21]. The signal injection method includes parallel resistance method and zero-sequence current disturbance method [
30,
31,
32]. This method demonstrates good ability in determining the line with persistent faults by injecting a large current using a dedicated or other primary equipment. Nevertheless, the above-mentioned methods all require primary feeder sensors and relevant communication capabilities for signal transmission, while the success rate of fault detection often depends on number of sensors deployed. In practice, the substantial investment required impacts the implementation of these methods, power companies often try to locate fault line through manual testing, resulting in service cut of non-fault line and down-graded power quality [
33].
In recent years, the mass deployment of low-voltage (LV) sensors, especially customer smart meters and smart distribution terminals installed at LV side of distribution transformers, provides power companies with full visibility of distribution network, which also opens a new way for detecting faults on the MV side through measurement of the multiple LV sensors on this MV feeder [
34,
35,
36,
37]. In [
34,
35], LV voltage imbalance is compared to determine the location of faults on the MV side. Although it cannot identify short circuit faults, it performs well on detecting high impedance faults such as broken wires. In [
36,
37], the author proposes a LV steady-state measurement-based method, which can accurately determine the fault location after signal classification through clustering algorithm, and it works for different neutral point grounding system, but the weak change in positive-sequence and negative-sequence steady-state LV impacts its practical applications for detecting single-phase-to-grounding fault in neutral non-effectively grounded systems.
When a single-phase-to-ground fault occurs on primary feeder or the MV side of the distribution line in neutral non-effectively grounded system, the transient aerial mode voltage (AMV) changes significantly at fault moment and demonstrates a stable pattern. On non-fault lines, the AMV magnitude decreases along the line, while on the fault line it has the highest magnitude in the fault point and decreases along the feeder in both sides. Based on this, this paper proposes a new such fault detection method by monitoring and analyzing the transient signals through LV sensors. By using the AMV for fault detection, not only can the primary fault line be accurately identified, but the primary fault location can also be determined. The main contributions are as follows:
1) A novel method for monitoring single-phase-to-ground fault in MV side of the distribution network through LV sensors is proposed. The fault transient signal collected by the LV sensor is adopted for fault characteristic analysis and decision making, which demonstrate great robustness.
2) The characteristics of transient AMV at different LV sensing points after MV single-phase-to-ground fault is analyzed and derived, which enables the fault location to a reasonable range. It is found that on the non-fault line, the magnitude decreases along the feeder to the end; while on the fault line, it has the highest value in the fault point and decreases along the feeder both sides.
3) The fault line and fault location can be determined by analyzing the magnitudes of the transient AMV sensed by multiple LV sensors.
Leveraging existing LV sensors, the proposed method can be used as independent fault detection when there is no fault detection equipment on the MV side thus avoiding MV sensor capital investment, or as supplemental measure cooperating with the existing MV detection equipment. It also works for both effectively grounded system and non-effectively grounded system.
2. Fault Characteristics Analysis
In analyzing the transient process of single-phase-to-ground fault in non-effectively grounded system, the Karenbauer phase-mode transformation method, which is also known as the 1, 2, 0 transformation method, is used. The 1-mode component and the 2-mode component are collectively referred to as the aerial-mode component, and the network composed of AB phase and AC phase, respectively. The parameters are same as those of the positive-sequence network in the fault symmetric component analysis; The 0-mode component, also known as the ground-mode component, forms a network between the three-phase conductor and the earth, so it has the same parameters as the zero-sequence network in the fault symmetric component analysis [
38]. For distribution transient fault signals within 3000Hz, the compound mode network of single-phase-to-ground fault through arc suppression coil (ASC) grounding system obtained by the mode transformation method, is shown as in the
Figure 1 below [
39].
In the figure, e(t) is the power supply voltage, Rk is the grounding resistance, and Lp is the ASC inductance; Ls1 and Rs1 in the 1-mode network are respectively the aerial-mode inductance and resistance of the power supply system behind the bus, and LL1 and RL1 are respectively the aerial-mode inductance and resistance from the fault point of the fault line to the bus; Ls2、Rs2、LL2 and RL2 in the 2-mode network have the same meaning corresponding to 1-mode network; LL0 and RL0 are the 0-mode inductance and resistance from the fault point of the fault line to the bus, and C0 is the sum of the 0-mode capacitance of the bus and its power source, the 0-mode distributed capacitance of all non-fault lines, and the 0-mode distributed capacitance from the fault point of the fault line to the bus, which is equal to the 0-mode distributed capacitance of the whole system minus the 0-mode distributed capacitance of the line after the fault point.
During the fault transient process, let
e(
t) =
Umcos
ωt,
ω is the power frequency, and the transient current in the compound mode network is:
where
is the attenuation coefficient of resonant component,
R=2(
Rs1+
RL1)+
RL0+3
Rk,
L=2(
Ls1+
LL1)+
LL0;
is the resonance angular frequency of the network,
UmωC0 is the peak value of steady-state capacitive current
Ic. Formula (1) can be simplified as:
In (2), as
,
is approximately equal to
ω0, thus
ω0 is referred to as transient dominant resonant angular frequency, and its value is generally 300~2000Hz in the distribution network [
39]. When
t equals to a quarter period of the transient dominant resonant frequency,
, and
, the formula (2) can be simplified as follows:
In (3), the value of sin
ωt is relatively small; when sin
ω0t is at its peak moment, the instantaneous value of
i(
t) is close to its maximum value, and the ratio of the peak value to the magnitude of the system steady-state capacitive current to the ground is approximately the ratio of transient dominant frequency
ω0 to power frequency
ω. Generally, the ratio can reach more than ten times of the steady-state magnitude [
33,
39], thus this transient feature has good observability in LV measurement which can be used to detect the fault occurrence.
It can be further known from
Figure 1, at the initial stage of the transient process
t = 0, the transient current is composed of capacitance charging/discharging current and inductance current. No matter how the neutral point is grounded, there is an equivalent capacitance
C0 in grounding system that affects the transient current of the compound mode network. According to formula (1), in the initial stage of the transient process, the characteristics of the transient current are mainly determined by the transient capacitive current, which is much larger than the transient inductance current. Therefore, the distribution of the transient voltage depends on the capacitance value which also determines the maximum transient voltage value. In the MV distribution system, the transient voltage value dominated by
C0 is relatively high, basically maintaining at the kilovolt level; After the transition to the steady-state, the voltage is distributed by the inductance
Lp, and transits to the hundred-volt level [
33,
39], which is far smaller than the that in initial transient process, so it is much more reliable to detect the fault using the transient signals compared to steady signals. When a single-phase-to-ground fault occurs, with the same power source, the equivalent capacitance to the ground of the non-fault line dominates its transient AMV instantaneous value. As the distribution feeders are generally short, the capacitance to the ground of each feeder has no obvious difference, resulting in small AMV instantaneous value on the non-fault line and this value decreases along the line; The AMV instantaneous value on the fault line is determined by the capacitance of all lines to the ground, which is much larger than that of the non-fault line. Therefore, the transient fault signals on the MV side are regularly distributed along the line. The transient signal within 3000Hz has little attenuation through the distribution transformer, thus it can be transmitted to LV distribution network [
39,
40] and caught in the secondary side of distribution transformer. Based on the above analysis, fault on the MV feeder can be detected with LV measurements at each point of the distribution feeder.
After identifying the fault line, the location of the fault point can be further determined through the measurement data of all LV sensors on the line. Considering aerial-mode parameters of the downstream line and load at the fault point in
Figure 1, the 1-model network is shown in the following
Figure 2.
The voltage
Ul1 of network in
Figure 2 is:
In (4),
i1 is the transient current. From
Figure 1 and
Figure 2, it can be seen that
Ul1 is the aerial-mode component of the fault transient voltage, which is the AMV value at the fault point. In order to analyze the distribution pattern of AMV along the fault line, the load side line parameters after the fault point are supplemented in
Figure 2. For the convenience of analyzing, the impact of ground capacitance is ignored [
41], the fault AMV distribution on a single line can be approximately represented by
Figure 3.
In the figure,
es is the virtual power supply at the fault point,
Rf and
I are the aerial-mode component of grounding resistance and fault current,
Is and
Il are the power source side and load side components of
I;
Zs is the system aerial-mode impedance, with value as the sum of
Ls1,
Rs1,
Ls2,
Rs2 and the equivalent capacitance of the power source;
Zkf is the aerial-mode impedance from the bus to the fault point of the fault line, with value as the sum of
LL1,
RL1,
LL2,
RL2 and the equivalent capacitance from the fault point to the bus of the fault line;
Zkl is the aerial-mode impedance of other parts of the fault line;
Zl is the aerial-mode impedance of the fault line load. For a more intuitive analysis of AMV distribution along the feeder,
Zkf and
Zkl are divided into several parts to represent the feeder segments marked by distribution transformer locations on the line, and the results are shown in the following
Figure 4.
Zkf1~
Zkfn and
Zkl1~
Zklm are the aerial-mode impedances of the lines between each distribution transformer, denoted as 1, 2, ...,
n and 1, 2, ...,
m from near to far from the fault point; The AMVs at the fault point, at the beginning and end of the line are
Uk,
Us, and
Ul.
In
Figure 4, the transient voltages
Uf1 and
Ul1 on both sides of the transformer adjacent to the fault point are obtained from Kirchhoff's voltage law (KVL) as follows:
It can be inferred that, the
Uk value at the fault point is highest, while the transient voltage value gradually decreases as the distance increases on both sides of the fault point along the feeder. The value of transition resistance
Rf at the fault point will affect the magnitude of the voltage
Uk at the fault point, but will not change the voltage distribution pattern in
Figure 4. Since the load aerial-mode impedance
Zl is much greater than the system aerial-mode impedance
Zs, the aerial-mode impedance on both sides of the fault point meets:
Therefore, the relevant magnitudes between the AMVs
Uk,
Us, and
Ul at the fault point, the beginning and end point of the line, can be illustratively shown as the following
Figure 5:
The distribution characteristics of AMV is applicable to common grounding system. According to the location of the distribution transformer on the feeder, the magnitude of the AMV signal also follows the similar distribution pattern. At the LV side, the AMV of the fault line is significantly higher than that of other non-fault lines. On the fault line, the AMV at the fault point is higher than other points on its both sides. Based on this, a single-phase grounding fault occurred on the MV side can be detected, and furthermore the fault line and location can be pinpointed.
3. Simulation Result and Analysis
Simulation is carried out using SIMULINK on a system as shown in
Figure 6. There are three 10kV feeders
L1,
L2 and
L3 in the system, with multiple distribution transformers. The transformer winding in the substation is DYn, and the neutral point grounding method in 10kV is configurable. Both 10kV lines and LV loads adopt RLC parameter models, and the distribution transformer windings are set as DYn.
Select L1 as the fault line, two single-phase grounding fault scenarios are studied with various neutral grounding modes and grounding resistances. Faults 1 and 2 are located on the primary feeder segment in the middle of distribution transformer 1 and 2, transformer 4 and 5 respectively.
When the neutral point of the system is grounded through the ASC and the fault resistance is 0Ω, assume fault 1 occurs at 0.2s.
Figure 7 shows the AMV change in the LV side of the transformers 1-5 on the fault line
L1, and
Figure 8 shows that in the LV side of the transformers 6-7 on the non-fault line
L2 and
L3.
From
Figure 7-8, it can be seen that when the neutral point is grounded through the ASC and the fault resistance is 0Ω, the AMV on the LV side of the fault line shows sharp change after the fault occurrence, and the magnitude of its instantaneous value is significant enough to be reliably detected by practical sensors. Although the AMV on the LV side of the non-fault lines changes as well at the same time, the magnitude of the AMV in the fault line is much larger than that of the non-fault line, and the difference is very noticeable. From amplitude perspective in
Table 1, The AMV has a maximum value of 699.7V at transformer 2, and 694.9V at transformer 1, while the results in transformer 3, 4, and 5 are significantly lower. Therefore, based on the comparison of amplitude results, the MV fault section can be judged as locates between transformer 1 and 2.
Table 1 also shows simulation results of fault 1 with different neutral grounding modes and grounding resistances.
Table 2 shows the results for fault 2 of all AMVs with different fault resistances and neutral grounding modes. Taking one example, when the neutral point is grounded through the ASC and the fault resistance is 0Ω, and the AMV has a maximum value of 680.6V at transformer 4, and 677.2V at transformer 5, which indicates the MV fault is located between these two transformers.
From the results shown in
Table 1 and
Table 2, it can be seen that although the AMV magnitude change significantly along with the changing the fault resistance and neutral grounding mode, it follows the distribution pattern as analyzed. Fault line
L1 can be clearly identified since its distribution transformers have noticeable higher AMV amplitude than transformers on
L2,
L3 in all cases. The MV fault location of fault line
L1 can also be pinpointed by finding its two distribution transformers with highest AMV amplitude, which are transformer 2 and 1 for fault 1, transformer 4 and 5 for fault 2. However, it can be seen that the AMV amplitude in fault segment decreases as the fault grounding resistance increases, which infers a limitation of the proposed method in high impedance fault when the AMV amplitude differences between fault line and non-fault line are not that distinguishable with practical sensor measurements.
4. Lab Experiment Result and Analysis
The proposed method is also tested in lab system shown in
Figure 9, in which power supply is provided by a 220kV transformer. There are two 10kV/400V distribution transformers M, N on the 10kV distribution feeder, with RLC adjustable loads simulating transformer loads.
The fault 3 is set on 10kV feeder primary upstream of distribution transformer M, and the fault signals are picked up through wave recorder installed at LV side of distribution transformers. When metallic single-phase-to-ground fault is triggered at phase A of the feeder, the neutral point is grounded through ASC, and AMV changes in LV side of distribution transformers are shown in
Figure 10.
Then the fault 4 is set downstream of distribution transformer N while other conditions remain unchanged, and AMV changes in LV side are shown in
Figure 11.
Table 3 shows the results of the AMV amplitude for every transformer in fault 3 and fault 4.
It can be seen that when the fault point is located upstream of the transformer M (fault 3), the amplitude of the AMV on the LV side of transformer M near the fault point is 144.98V, which is greater than 123.62V at transformer N which is further from the fault point. Similarly, after changing the location of the fault point (fault 4), the amplitude of the AMV in transformer N is 78.73V, which is greater than 69.64V in transformer M. It can also be noted that the fault location change leads to slightly change of AMV. The results conform with the characteristics proving the effectiveness of the proposed method.
Table 4 lists the lab experiment results with different combinations of fault resistances and neutral grounding methods. Constrained by the lab facilities, only two scenarios for the fault resistances (0Ω and 500Ω) are tested. In the test, changing the fault resistance, neutral grounding mode and fault point location leads to the change in AMV amplitude. The results confirm that the amplitude of the AMV near the fault point is significantly larger than that located further away, indicating that effectiveness of the proposed method in locating the MV fault section by comparing the amplitude of the AMV. However, in the case, with the neutral point grounded through ASC when the fault resistance is 500Ω, although the AMV amplitude at transformer M is still greater than that of transformer N, the amplitude of AMV is too low, which also indicates that limitation of the proposed method for high impedance fault.