Next Article in Journal
Electrolyte Leakage in Cylindrical Lithium-Ion Batteries Subjected to Temperature Cycling
Next Article in Special Issue
A Grouping and Aggregation Modeling Method of Induction Motors for Transient Voltage Stability Analysis
Previous Article in Journal
Development of Correlations Based on CFD Study for Microchannel Condensation Flow of Environmentally Friendly Hydrocarbon Refrigerants
Previous Article in Special Issue
Optimization of Ampacity in High-Voltage Underground Cables with Thermal Backfill Using Dynamic PSO and Adaptive Strategies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ground Fault in Medium-Voltage Power Networks with an Isolated Neutral Point: Spectral and Wavelet Analysis of Selected Cases in an Example Industrial Network Modeled in the ATP-EMTP Package

Department of Electrical Engineering and Automation in Industry, Faculty of Mining, Safety Engineering and Industrial Automation, Silesian University of Technology, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Energies 2024, 17(7), 1532; https://doi.org/10.3390/en17071532
Submission received: 27 February 2024 / Revised: 18 March 2024 / Accepted: 20 March 2024 / Published: 22 March 2024
(This article belongs to the Special Issue Modeling, Simulation and Optimization of Power System)

Abstract

:
The paper presents some case spectral analysis of zero-sequence voltages and currents in an example industrial power distribution network. The network layout is based on typical power delivery networks in underground coal mines. Ground fault simulations have been made using an ATP/EMTP program. Due to the high environmental risks, the reliability of the protection relay operation related to their selectivity plays an important role. This paper tries to find the reasons for nonselective operation and unnecessary tripping in extensive mine cable networks, particularly with large power sources of higher-order harmonics. It was found that in transient states—due to the decaying oscillations occurring in complex RLC circuits—the results of short time measurements of the criterion values for ground fault protective relays can be overestimated (particularly for small values of ground resistance) and lead to nonselective tripping of a healthy cable line. Therefore, it might be advisable to increase the integration time used for measuring rms values. Also, if there is a significant level of higher harmonics in the industrial network generated by high-power converters, it should be noted that the higher harmonics of the ground fault current and currents measured by ground fault protection relays assume much higher values, which may also cause nonselective tripping. In this case, it may be advisable to use higher harmonic filters in the measuring circuits and to select a sufficiently high sampling frequency in the digital protective relays.

1. Introduction

The method of grounding the neutral point of the transformer significantly affects the operational properties of the power distribution network. This applies to both the effects of possible ground fault phenomena and to the operation of ground fault protection relays. The most commonly used networks are those with a neutral point grounded directly or indirectly (through resistance or impedance) and with an insulated neutral point. The grounding method used in a particular network often results from historical events and requirements regarding the values of ground fault currents and ground fault overvoltages; therefore, different solutions can be found in individual countries and in different application areas (industry branches). Mining electrical equipment usually operates in very harsh environmental conditions (dust, humidity, methane). This leads to increased fire, shock, and explosion hazards in underground mines [1,2,3]. Therefore, power distribution networks in underground mines often operate with an isolated neutral point, as this configuration significantly reduces the current and energy values at the fault location compared to those occurring in networks with a directly or indirectly grounded neutral point. The reduction in the current value at the fault location occurs due to the fact that this current is determined not so much by the small series impedances of the network elements, as in the case of two-phase and three-phase short circuits, but by the much greater shunt impedances resulting from the earth capacitance of all galvanically connected cable lines. However, a negative aspect of this reduction in ground fault current is the difficulty of detecting these faults quickly and reliably and of selectively (and possibly quickly) tripping only the faulty cable line. This enforces the need for more sophisticated zero-sequence voltage and current measurement systems (like Ferranti transducers for zero-sequence current component measurement, and open delta connected voltage transformers for zero-sequence voltage component measurement) and protective relays. Due to the relatively low values of the zero-sequence currents and voltages in an isolated neutral point network, these measuring systems and protection devices are much more sensitive to disturbances and require more complex calculations and simulations taking into account transients and deformations of the zero-sequence voltages and currents caused by power electronic converters. The complexity of ground fault phenomena (especially transients and nonlinear arc resistances) drives the use of complex multiresolution digital signal processing methods based on wavelet transforms [4,5,6] and empirical mode decomposition (EMD) [7]. It should be noted, however, that despite obvious advances in the use of microprocessor technology and digital signal processing in protective relaying, practically applied protection in hazardous areas is still based on quite simple, strict, and conservative standards and regulations (like directional or nondirectional overcurrent relays or admittance relays). The discrete wavelet transform (DWT) in the applications mentioned above can be considered a very efficient method for digital signal processing (e.g., in the implementation of protection algorithms); however, in this paper, we use the continuous wavelet transform (CWT) as a method to analyze and visualize transient ground fault current waveforms jointly in the time and frequency domain in the first period after ground fault initiation. A comprehensive literature review on ground fault phenomena and their detection and classification is available in [8]. The analytical description of transient states in ground faults is usually focused on the issue of ground fault overvoltages and is based on simple second-order LC models described by a single transient frequency. In our paper, we present a novel two-frequency model described by the fourth-order differential equation. Also, the results of numerical analyses carried out with the use of simulation software were innovatively visualized as scalograms with the use of the continuous wavelet transform—which enabled the presentation of frequency fluctuations and waveform in the first period after the ground fault initiation.
Many research works are devoted to the analysis of ground fault phenomena in MV networks with resonant-grounded [9,10,11] or resistor-grounded [12] neutral and containing overhead lines [4,13,14,15]. However, the operating conditions of these networks differ significantly from those of distribution networks installed in underground mines [1]. In this paper, the focus is on networks with an isolated neutral point containing only cable lines. It is noteworthy that, in order to mitigate the risks associated with high-current value phase-to-phase short circuits, the MV cables used in underground mines are armored and shielded cables. This significantly increases their line-to-ground capacitance [16] and thus their specific ground current value per unit length [1].

2. Materials and Methods

2.1. Theoretical Analysis and Simulation Software

A simplified representation of an unbranched network with an isolated neutral point is depicted in Figure 1. Given the significantly smaller values of the series line impedances compared to the capacitive line-to-ground impedances, the former can be neglected. Similarly, with the network insulation resistance assuming considerably larger values in comparison to the (connected in parallel) capacitive line-to-ground impedances, the latter can also be omitted, as shown in Figure 1b. This simplification facilitates the derivation of the impedance of a Thevenin generator (Figure 1c) and a straightforward equivalent circuit (Figure 1d).
Consequently, this enables a simplified evaluation of ground fault current values across various frequencies, encompassing higher-order harmonics of the fundamental frequency waveform.
Given that the ground current value (IF) can be expressed according to the simplified Equation (1) derived from Figure 1d, it is evident that the resulting circuit impedance diminishes with increasing frequencies, including higher-order harmonics.
I F ( ω ) = V Z ( ω ) ¯ = V R F 2 + 1 3 ω C G 2 = 3 V ω C G 1 + 9 ω 2 R F 2 C G 2 ,
where IF—line-to-ground fault current, V—line-to-ground source voltage, RF—line-to-ground fault resistance, CG—line-to-ground capacitance, RG—line-to-ground insulation resistance, and ω—angular frequency of the source voltage.
This phenomenon is especially noticeable for zero or negligible values of the ground resistance (RF) when Equation (1) can be simplified to Equation (2).
I F ( ω ) = ω 3 V C G
Under such conditions, even modest magnitudes of source voltage harmonics can induce relatively significant values of ground fault current harmonics.
In real networks with a more complex (branched) structure, the flow of currents is much more complicated, as shown in Figure 2. The line-to-ground voltage in the faulty phase is reduced, and those in the other undamaged phases are elevated. Therefore, the damaged phase capacity is discharged (red arrows in Figure 2), and the remaining (healthy) phases’ capacities are additionally charged (blue arrows in Figure 2). The magnitudes of the zero-sequence currents measured by the protection relays are determined by the capacitances of the entire galvanically connected network (Figure 3).
Due to the large number of energy storage elements (including capacitances and inductances), transient dynamics is described by intricate systems of high-order differential equations posing challenges for analytical solutions. Therefore, the adoption of specialized numerical simulation software emerges as a viable methodological approach to dissect and comprehend the intricate phenomena within such a network. We used an ATP/EMTP software package [17]. This package consists of three tools: the interactive graphical user interface (ATPDraw v. 6.3), the simulation tool (TPBIG), and a postprocessor for further analysis and graphical output (PlotXY or Matlab v. R2020a in our case). The whole toolchain is shown in Figure 4.
The Π cable model [17] was used to represent the cable lines. The network design was carried out in the ATP/EMTP package, allowing precise simulation of system behavior. The ATP/EMTP package is widely used for the numerical modeling of transients in power networks [17], including the modeling and testing of the operation of protection relay systems [18], harmonics propagation [19], and power quality problems [20,21]. The ATP/EMTP software is often used to generate data used in the development of new algorithms for the detection and location of ground faults [22,23,24,25]. Other software such as PSCAD [6,26] or Matlab [6,11] can be also used for the modeling and analysis of ground fault phenomena.

2.2. Transient Components of the Ground Fault Current and Their Frequencies and Decay

When a ground fault is initiated, a transient occurs. As the network has a complex layout with a large number of energy-accumulating (capacitances and inductances) and energy-dissipating (resistances) elements, the transient waveforms are described by very high-order differential equations, which can only be solved numerically. An approximate analytical solution is only possible if certain simplifying assumptions are made. In this study, the network model shown in Figure 5a is adopted, with the inductance of the power source (transformer), the ground capacitance of the faulted line, the ground capacitance of the rest of the network, and the resistances and inductances of the faulted line as the most important elements. This model structure allows the earth capacitances of the faulty line and the rest of the network to be connected in parallel, as shown in Figure 5b.
The transient analysis of this model using the Laplace transform is possible based on a superposition method allowing the calculation of the steady-state (easily determined) component waveform to be separated from the calculation of the transient component waveform. The equivalent operator impedance circuit for the transient component and its subsequent simplifications are shown in Figure 6a–d.
The equivalent operator reactance Xz(s) of the parallel LC branch is determined as
X Z s = 1 s C 0 Σ 1 2 s C 0 Σ + 3 2 s L T 1 2 s C 0 Σ + 1 s C 0 Σ + 3 2 s L T = 1 2 s C 0 Σ + 3 2 s L T s C 0 Σ 3 2 s C 0 Σ + 3 2 s L T = 1 s C 0 Σ + 3 s L T s C 0 Σ 3 s C 0 Σ + 3 s L T ,
where LT—inductance of power transformer windings, and C—total capacitance of the entire galvanically connected network.
Therefore,
X Z s = 1 + 3 s 2 C 0 Σ L T 3 s C 0 Σ + 3 s 3 C 0 Σ 2 L T ,
and to simplify the notation, the following RΣ symbol will be adopted:
R Σ = R L + R F .
Then, the total operator impedance of the entire circuit for the transient component ZΣ(s) can be determined as
Z Σ s = R Σ + s L L + X Z s = R Σ + s L L + 1 + 3 s 2 C 0 Σ L T 3 s C 0 Σ + 3 s 3 C 0 Σ 2 L T .
So
Z Σ s = 3 C 0 Σ 2 L T L L s 4 + 3 C 0 Σ 2 L T R Σ s 3 + 3 C 0 Σ L T + L L s 2 + 3 C 0 Σ R Σ s + 1 3 C 0 Σ 2 L T s 3 + 3 C 0 Σ s ,
where LL—faulted line inductance, and RΣ—total resistance of the whole circuit (including both cable line cores and line-to-ground fault resistance).
The waveform of the transient component of the earth current I0(s) can therefore be determined from the following equation:
I 0 T s = V 0 s 3 C 0 Σ 2 L T s 3 + 3 C 0 Σ s 3 C 0 Σ 2 L T L L s 4 + 3 C 0 Σ 2 L T R Σ s 3 + 3 C 0 Σ L T + L L s 2 + 3 C 0 Σ R Σ s + 1 .
The nature of the waveforms contained in the transient component will be determined by the values of the roots of the characteristic polynomial, i.e., the solutions of the equation
Q s = 0 ,
where the characteristic polynomial Q(s) is
Q s = 3 C 0 Σ 2 L T L L s 4 + 3 C 0 Σ 2 L T R Σ s 3 + 3 C 0 Σ L T + L L s 2 + 3 C 0 Σ R Σ s + 1 .
Equation (10) describes a polynomial of order 4 that is difficult to solve analytically. However, it can be noticed that for negligibly small values of RΣ, Equation (10) takes the form of a biquadratic Equation (11), which is much simpler to solve:
3 C 0 Σ 2 L T L L s 4 + 3 C 0 Σ L T + L L s 2 + 1 = 0 .
Assuming an auxiliary variable
q = s 2 ,
Equation (11) can be written as
3 C 0 Σ 2 L T L L q 2 + 3 C 0 Σ L T + L L q + 1 = 0 .
The solution of Equation (13) are two negative-valued roots q1 and q2
q 1 = 3 C 0 Σ L T + L L + 9 C 0 Σ 2 L T + L L 2 12 C 0 Σ 2 L T L L 6 C 0 Σ 2 L T L L ,
q 2 = 3 C 0 Σ L T + L L 9 C 0 Σ 2 L T + L L 2 12 C 0 Σ 2 L T L L 6 C 0 Σ 2 L T L L ,
and after some simplification,
q 1 = 2 3 C 0 Σ L T + L L 3 C 0 Σ 3 L T + L L 2 4 L T L L ,
q 2 = 2 3 C 0 Σ L T + L L + 3 C 0 Σ 3 L T + L L 2 4 L T L L .
Equation (11) will therefore have two pairs of conjugate imaginary roots corresponding to angular frequencies:
ω 1 = 2 3 C 0 Σ L T + L L 3 C 0 Σ 3 L T + L L 2 4 L T L L ,
ω 2 = 2 3 C 0 Σ L T + L L + 3 C 0 Σ 3 L T + L L 2 4 L T L L .
This corresponds to the frequencies of
f 1 = 1 2 π 2 3 C 0 Σ L T + L L 3 C 0 Σ 3 L T + L L 2 4 L T L L = 1 2 π 3 C 0 Σ L T + L L 3 C 0 Σ 3 L T + L L 2 4 L T L L ,
f 2 = 1 2 π 2 3 C 0 Σ L T + L L + 3 C 0 Σ 3 L T 2 + L L 2 4 L T L L = 1 2 π 3 C 0 Σ L T + L L + 3 C 0 Σ 3 L T + L L 2 4 L T L L .
If an additional simplification is adopted to neglect the faulted feeder line inductance LL, then Equation (11) takes the very simple form:
3 C 0 Σ L T s 2 + 1 = 0 ,
whose solution is a single-frequency waveform with a frequency equal to (corresponding in this case to Equation (21))
f = 1 2 π 3 C 0 Σ L T .
Neglecting the value of the RΣ resistance significantly facilitates the solution of Equation (10) but makes it impossible to determine directly the coefficients describing the decay of the oscillatory transient waveform components. However, these decay coefficients can be taken into account as a small perturbation [27,28] of Equation (11).
Once the angular frequencies ω1 and ω2 have been determined in accordance with Equations (18) and (19), Equation (11) can be written as a product of linear factors:
α s j ω 1 s + j ω 1 s j ω 2 s + j ω 2 = 0 ,
where the value of the α factor (introduced to shorten the notation) will be equal to
α = 3 C 0 Σ 2 L T L L ,
expressed in s4. Equation (24) can be there simplified to the form
α s 2 + ω 1 2 s 2 + ω 2 2 = 0 .
If the RΣ resistance was included in Equation (10), the transient components of the oscillatory waveforms with frequencies f1 and f2, being its solution, would be damped sinusoidal waveforms, and the pairs of conjugated complex roots of Equation (10) would take the general form of
s 1 , 2 = σ 1 ± j ω 1 s 3 , 4 = σ 2 ± j ω 2 .
The polynomial Equation (24) can therefore be written as
α s σ 1 + j ω 1 s σ 1 j ω 1 s σ 2 + j ω 2 s σ 2 j ω 2 = 0 ,
so
α s + σ 1 j ω 1 s + σ 1 + j ω 1 s + σ 2 j ω 2 s + σ 2 + j ω 2 = 0 ,
and
α s + σ 1 2 + ω 1 2 s + σ 2 2 + ω 2 2 = 0 .
Therefore,
α s 2 + 2 σ 1 s + σ 1 2 + ω 1 2 s 2 + 2 σ 2 s + σ 2 2 + ω 2 2 = 0 ,
and, finally,
α s 4 + 2 s 3 σ 1 + σ 2 + s 2 σ 1 2 + σ 2 2 + 4 σ 1 σ 2 + ω 1 2 + ω 2 2 + + 2 s σ 1 σ 2 2 + σ 1 ω 2 2 + σ 2 σ 1 2 + σ 2 ω 1 2 + σ 1 2 σ 2 2 + σ 1 2 ω 2 2 + σ 2 2 ω 1 2 + ω 1 2 ω 2 2 = 0 .
At the same time, the polynomial Equation (26) can be reduced to the form
α s 4 + s 2 ω 1 2 + ω 2 2 + ω 1 2 ω 2 2 = 0 .
Subtracting Equation (33) from Equation (32) yields
α 2 s 3 σ 1 + σ 2 + s 2 σ 1 2 + σ 2 2 + 4 σ 1 σ 2 + + 2 s σ 1 σ 2 2 + σ 1 ω 2 2 + σ 2 σ 1 2 + σ 2 ω 1 2 + σ 1 2 σ 2 2 + σ 1 2 ω 2 2 + σ 2 2 ω 1 2 = 0 .
By expanding the polynomial Q(s) into a first-order Taylor series around the points ω1, ω2 (i.e., omitting the second- and higher-order terms involving σ12, σ22, σ1, σ2), we obtain
α 2 s 3 σ 1 + σ 2 + 2 s σ 1 ω 2 2 + σ 2 ω 1 2 = 0 .
By comparing the polynomial coefficients appearing in Equations (35) and (10) at the respective powers (s and s3), the following system of equations is obtained:
2 α σ 1 + 2 α σ 2 = 3 C 0 Σ 2 L T R Σ 2 α σ 2 ω 1 2 + 2 α σ 1 ω 2 2 = 3 C 0 Σ R Σ .
Therefore,
σ 1 + σ 2 = 3 2 α C 0 Σ 2 L T R Σ σ 2 ω 1 2 + σ 1 ω 2 2 = 3 2 α C 0 Σ R Σ .
The solution to the system of Equation (37) will be of the form
σ 1 = 3 C 0 Σ R Σ C 0 Σ L T ω 1 2 1 2 α ω 1 2 ω 2 2 ,
σ 2 = 3 C 0 Σ R Σ 1 C 0 Σ L T ω 2 2 2 α ω 1 2 ω 2 2 .
By inserting the α value described by Formula (25) and the ω1 and ω2 values described by Formulas (18) and (19) into the solution (38) and (39), we obtain, after some transformations,
σ 1 = R Σ 3 L T 3 L L + 9 L T 2 + 6 L T L L + 9 L L 2 4 L L 9 L T 2 + 6 L T L L + 9 L L 2 ,
σ 2 = R Σ 3 L L 3 L T + 9 L T 2 + 6 L T L L + 9 L L 2 4 L L 9 L T 2 + 6 L T L L + 9 L L 2 .
For the industrial network considered in the further simulation (ATP/EMTP) studies with the aggregated parameters C = 6.048 μF, LL = 2.9808 mH, LT = 1.109 mH, and RΣ = 1.07028 Ω, the above calculations were performed to show the following results:
  • f1 = 2194 Hz and σ1 = 39∙RΣ = 41.74 s−1 for the 1st transient component.
  • f2 = 606 Hz and σ2 = 128∙RΣ = 137.0 s−1 for the 2nd transient component.
The lower-frequency (f2) transient component therefore decays more than three times faster than the higher-frequency (f1) transient component. This phenomenon is shown in Figure 7. The decay of both transient components is, of course, even faster if a higher line-to-ground (RF) resistance value is present at the ground fault location.
The model considered above is a simplified model, which does not take into account the cable line frequency-dependent parameters (skin-effect phenomenon), and the obtained results are conditioned by the specific nature of the analyzed network—i.e., power supply from an HV/MV transformer with large rated power (and therefore the associated low impedance) and high inductance of the considered feeder cable line (due to its long length and the use of shielded cables). Due to the use of individually shielded cable cores, line-to-line capacitances were neglected in the model under consideration. High values of transient waveform frequencies may affect the operation of modern protection devices using digital signal processing, and knowledge of the decay rate of transient waveforms can facilitate the correct selection of the time delay.

2.3. Illustrative Case Study—An Example Network

A simplified MV network model was used to analyze ground fault disturbances in medium-voltage industrial distribution power grids with an isolated neutral point. The network layout was designed on the basis of several typical MV network layouts of Polish coal mines where ground faults have been recorded in the period 2019–2021. The key aspects of this model and the simulation procedures that were used in the study are presented here.
In the examined model, all the network feeders were assumed to be cable lines, consisting of shielded mining cable type YHKGXSFoyn with a cross-section of 3 × 185 mm2. The total length of all cables in the galvanically connected network was 10.8 km, and the ground fault current value was 19.76 A. Fault current simulations were carried out considering a fault in a cable leaving from panel 1 of the DS-1 departmental switchgear.
During the experiments, zero-sequence current i0(t) and zero-sequence voltage u0(t) waveform data have been recorded in different switchgear panels in the analyzed network model:
  • Main switchgear MS-1, panel 1;
  • Main switchgear MS-1, panel 2;
  • Departmental switchgear DS-1, panel 1;
  • Departmental switchgear DS-1, panel 2;
  • Departmental switchgear DS-1, panel 3.
The following assumptions have been made for all the simulations performed in the experiment:
  • The single line-to-ground fault was forced for 30 ms after the start of the simulation;
  • The phase voltage of the faulty line reached the maximum possible value when the fault occurred;
  • The switchgear panel of the line in which the single line-to-ground fault occurred was not disconnected until the end of the simulation;
  • In order to analyze the steady-state values of zero-sequence current i0(t) and zero-sequence voltage u0(t) waveforms, the simulation lasted 500 ms;
  • Regardless of the initial value of the line-to-ground fault resistance, this resistance value was assumed to remain constant throughout the whole duration of the fault simulation;
  • A sampling rate of 109 samples per second was assumed in the program; this enabled accurate analysis of higher harmonics in the analyzed network;
  • The simulation starts from a pre-calculated steady state, and the 30 ms interval is used only to make the signal waveform plots exhibit a more clear difference between the pre-fault state and the after-fault waveforms. As for the moment of initiation of a short circuit, we assumed the most likely case: the insulation is damaged at the moment when the line-to-ground voltage reaches its maximum value.
For the network layout shown in Figure 8, a set of zero-sequence network voltage V0 and current I0 time waveforms have been recorded for faults occurring at selected points in panel 1 of DS-1 switchgear (0, 300, 600, 900, 1200, and 1500 m counting from the beginning of the feeder cable) with fault resistance values equal to 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000, 2000, 5000, and 10,000 Ω for each point, respectively.
The network diagram utilized in the simulation, generated using the ATPDraw graphical preprocessor, is depicted in Figure 9. The subsequent section of this paper provides illustrative results corresponding to the selected length values. These results are representative of typical underground MV distribution networks.
Next, the authors directed their attention toward analyzing the rms values of the selected current and voltage waveforms. An important parameter in this analysis was the duration of the measurement window (integration time). The following 3 specific values of the measurement time have been selected:
  • 20 ms—This is the duration of one period from the start of the line-to-ground fault, calculated from the moment the line-to-ground fault occurs. The timing duration is not selected by chance. In the context of a ground fault, most of the transient waveform peaks disappear after the initial half-period of the fundamental frequency waveform. Overcurrent protection relays measure the rms value of the current/voltage waveform after a certain duration. To ensure the correct measurement of the rms value, it must be recalculated after a full period, which amounts to 20 ms for a 50 Hz fundamental frequency.
  • 80 ms—This is the most common in the Polish industry; it is the duration of a ground fault before the field is disconnected, counted from the moment of line-to-ground fault occurrence. This stems from the safety imperative in the Polish mining industry, driven by factors such as methane hazards; the safety priority is to disconnect all faults, not only line-to-line but also ground faults, without any delay. It is accepted that the fault should be disconnected in less than 100 ms. Since the fault has to be disconnected faster than 100 ms and the protection algorithm has to assess an integral multiple of the fundamental frequency period, the maximum time window duration under these conditions is 80 ms.
  • 100 ms—The time required to determine the rms values of steady-state voltages and currents, taken at 400 ms after the line-to-ground fault initiation. The steady state denotes a voltage and current waveform sufficiently temporally distant from the ground fault initiation, allowing transient waveforms to decay by that time. The selected temporal reference is a trade-off, considering factors such as simulation time and the mitigation of fault-induced noise during transients.
  • In summary, for the computation of root mean square (rms) values of steady-state voltages and currents, it was decided to select a time instant 20 periods after the ground fault initiation. The choice of a measurement window duration of 100 milliseconds (equivalent to 5 periods of the fundamental frequency) aims to mitigate the potential effects of an unplanned spike resulting from the numerical simulation method. Integrating data over 5 periods effectively “smooths out” such possible signal disturbances, rendering them negligible for the subsequent analysis.
Making the above assumptions, further studies (including spectral analysis) were carried out for two options:
  • Transient state in the network without higher harmonics;
  • Steady state in a network with higher harmonics.
In the contemporary era of industrial networks using diverse types of power electronic equipment generating higher harmonics, it may seem useless to analyze the phenomena occurring in the network without taking into account higher harmonics. However, in the case of MV distribution networks, in which only the first harmonic is often present, such an approach simplifies further consideration of ground fault phenomena in this type of network and therefore can be useful as a benchmark for comparative analysis.
In the analysis of the case of an industrial network with sources of higher harmonics, the international standard [29] defining the maximum permissible level of higher harmonics in an industrial network was used (Table 1).
The magnitudes of the above harmonics have been set in accordance with the upper limits specified by the standard. In the aforementioned standard, an additional restriction appears, stating that the average value of higher harmonics of voltages measured during any 10 min interval should not exceed 8%. For the data adopted in the above table, this value is 9.27%; however, it is assumed that in the perspective of the 10 min measurement, the contribution of higher harmonics will not exceed the aforementioned 8%.

3. Results

3.1. Analysis of a Ground Fault Transient in the Model Network without Higher Harmonics

In order to obtain data for the analysis and classification of the types of line-to-ground faults, simulations were conducted using the ATP-EMTP software. The simulations involved inducing a single line-to-ground fault at selected locations in panel 1 of the DS-1 departmental switchgear, according to the network configuration shown in Figure 4. Several sets of network voltage and current waveforms for different fault types (with different resistance values) have been recorded.
The obtained results of measurements of the rms values of currents and voltages in the first 20 ms after the fault occurrence were compared with the rms values for the steady state, which made it possible to compare the phenomena occurring in the network during a ground fault transient and in the steady state. Two single line-to-ground fault locations were selected for analysis: at the beginning and at the end of the outgoing cable line in panel 1 of departmental switchgear DS-1 (for x = 0 m—blue graph and for x = 1500 m—red graph). It is noteworthy that all the other characteristics fall in between these two cases.
Analyzing the results obtained, we focused on the characteristic that determines the ratio of the rms value of the zero-sequence component of the V0 network voltage after the fault occurrence to the rms value of the voltage of any of the phases of VPh before the fault occurrence. This ratio is described by Equation (42).
V 0 = β V p h
As is the case in reference [28], it is a common practice to depict the dependence of the β coefficient value on fault resistance value in the form of linear characteristics, but here, we depict our results on a semi-logarithmic scale. The authors hope that this enhances the clarity of these plots and makes it more convenient to analyze and understand the presented phenomena (such as the classification of ground fault types based on the inflection points of the curves visible on the graphs). Within the framework of this study—covering a wide range of fault resistance values 1 Ω … 10 kΩ—these relationships were presented on a semi-logarithmic scale. On the basis of Equation (43),
V 0 V P h = f R F
there was introduced a classification of ground faults into three categories:
  • Low-resistance faults: Those characterized by a relatively low parameter β, ranging from 0.9 to 1 (for steady state). In the analyzed network, these are assumed to be faults whose line-to-ground resistance is between 0 and 80 Ω.
  • Medium-resistance faults: Those characterized by a β parameter in the range of 0.2 to 0.9 (for steady state). In the network, they are assumed to include faults whose line-to-ground resistance is between 80 and 800 Ω.
  • High-resistance faults: These are faults that are more difficult to detect and are characterized by a β parameter in the range of 0 to 0.2. For the analyzed network, it is assumed that these are faults whose line-to-ground resistance exceeds 800 Ω.
It is important that the above classification, although found in the reference [30], was modified and extended in this work to suit the needs of the analysis. In broad terms, ground faults in an IT network in a hazardous environment can be categorized into three distinct groups:
  • Low-resistive (easily detectable but with a high ground fault current magnitude and associated hazard level);
  • High resistive (with a low earth fault current magnitude and a reduced danger level but much more difficult to detect);
  • An intermediate group between the two defined above (moderate fault current magnitude, moderate detection difficulty, and moderate hazard level).
The scope of the analysis extends to high resistance values (comparable even to the resistance value of the human body in direct or indirect contact scenarios). Line-to-ground currents then assume fairly small magnitudes, theoretically comparable to leakage currents. However, it should be noted that owing to the specific environmental conditions described in the introduction, the regulations mandate stringent high insulation resistance value requirements (at least 1 kΩ/1 V, depending on the rated network voltage), so earth leakage currents tend to be significantly smaller and, above all, (like balanced load currents) are symmetrical, i.e., there is no resulting zero-sequence component. In contrast, ground fault currents exhibit pronounced asymmetry and are therefore more easily detectable by suitable filters (core balance current transformers) [31].
As evident from the characteristics shown in Figure 10, it can sometimes be problematic to select the settings for zero-sequence component current relay protection efficiently in some special cases. Setting the integration time to 20 ms may result in the healthy line being unnecessarily disconnected for low values of fault resistances in addition to the line in which the ground fault occurred. For instance, in the above example, setting the current–overcurrent protection setting to 4 A in bay 3 of the DS-1 switchgear may result in disconnection of this bay even though the short circuit occurred at the beginning of bay 1 of the same switchgear. In the discussed example, the fault resistance would be only a few ohms.
A partial remedy to this issue involves extending the integration time for the calculation of the root mean square (rms) value. Figure 11 shows the same results as Figure 10, however, with the difference that, in this case, the integration time is not 20 ms but 80 ms. The slight modification reduces the significance of the initial distortions for the individual measurements. With this change in integration time, the possibility of selectively disconnecting only the faulty line is increased. A consequence of this adjustment is a prolonged relay response time. However, as demonstrated in the preceding sections, the protective relays in the vast majority of cases successfully disconnect the faulty line within this extended time frame.
Through numerical simulation experiments and a review of the existing literature, it becomes evident that the root mean square (rms) value of individual parameters within the initial 20 ms of a short-circuit duration is notably affected by the distortion in currents and voltages. Example transient waveforms of zero-sequence currents and voltages are presented in Figure 12. Despite the short time duration, high-magnitude pulses and slowly decaying high-frequency zero-sequence current oscillations are distinctly observable within the first period after the occurrence of a ground fault.
According to the references [25,32,33], a network with an isolated neutral point can be treated as a more complex RLC circuit. In RLC circuits, the occurrence of decaying oscillations is possible after a sudden commutation (stepwise change in structure). These decaying oscillations should be regarded as nonperiodic time functions with finite energy.
Thus, it is established that the transient magnitude spectrum is continuous. To enhance clarity, a decision was made to depict (Figure 13) the amplitude spectrum in discrete form for harmonics from 2 to 41. For ease of comparison, the percentage of the higher harmonics in relation to the first harmonic is presented. To streamline the visual analysis of higher harmonics against the first harmonic, the latter has been intentionally omitted from the individual plots.
Transient ground fault waveforms contain, in addition to the step-changing fundamental frequency component, slowly decaying oscillating waveforms with a much higher frequency—clearly visible in Figure 12 (their frequencies are determined not by the source frequency but by the inductances and capacitances occurring in the network—as described in Section 2.2). Therefore, the continuous wavelet transform (CWT) may become a more convenient method for their analysis than the Fourier transform [34,35,36]. This is related to the fact that the Fourier transform is based on sine and cosine functions characterized by infinite (temporally unlimited) support; hence, it cannot provide sufficient localization of fast-changing components. The wavelet transform, on the other hand, is based on functions with a rapidly decaying (time-constrained) support and therefore allows a much more accurate analysis of the variation in the individual components of the relay protection criterion signal. Hence, in order to delve deeper into the electromagnetic transient processes described here, a continuous wavelet transform based on analytical Morse wavelets [37,38] was used in our case study. The continuous wavelet transform makes it possible to watch the fluctuation of the individual components over shorter periods of time.
Instead of using a single (but complex) function of one variable (frequency) ω based on complex exponential support, like in a classical Fourier analysis—or alternatively two orthogonal real-valued functions based on sine and cosine support functions of one variable (frequency) defined over the whole analyzed signal range—the continuous wavelet transform (CWT) allows us to describe a signal with a function of two variables t, a:
W ( t , a ) = 1 a S ( t ) ψ τ t a   d τ ,
where ψ(t) denotes the so-called mother wavelet, i.e., basis function fulfilling certain mathematical conditions [39] (e.g., finite energy, limited both in time and in frequency space), a is a scale (basis function dilation or constriction—corresponding to frequency change) parameter, and t is a translation (time shift of the analysis window) parameter. This signal description has a local character (valid within the area defined by the time-constrained support of the wavelet applied), thus enabling more accurate analysis of the variability of individual frequency components over shorter time segments. Local matching of wavelet and signal leads to a large transform value.
A wavelet analysis of the transient and steady state was performed using the Matlab Wavelet Toolbox [40]. The outcome of this analysis (scalogram) is depicted in Figure 14.
It can be seen that in the initial (transient) time interval, the contribution from higher harmonics becomes noticeable, and they decay to a negligible magnitude after the period of approx. 40 ms. On the other hand, the 50 Hz fundamental component reaches its final value after approx. 50 ms, which, taking into account the additional time required to measure its rms value, gives a tripping time of 70–90 ms.
It is obvious that at relatively high fault resistance values (RF > 100 Ω), the percentage of higher harmonics in the current waveform is small; hence, it was consciously decided to omit spectral analysis for high-resistance faults. The first noticeable fact that can be observed is the correspondence of the percentage spectral characteristics of the currents both qualitatively and quantitatively. This makes it possible to determine the spectral character of the fault current on the basis of the selected recorded zero-sequence current I0 waveform. It was known that for small fault resistances, the higher harmonics would have a significant effect on the IF current waveform during the initial stage of the fault transient. However, for low RF values, it is possible (in some cases) that the magnitude of higher harmonics may reach values that are even twice as high as the magnitude of the first harmonic of the zero-sequence current.
What is also important is the fact that the maximum of this value is achieved for the 21st harmonic. Typically, modern digital protection relays are used in mining switchgear samples and record voltage and current waveforms with a sampling frequency of between 1000 and 1600 samples/s. At a fundamental network frequency of 50 Hz, the 21st harmonic frequency is 1050 Hz. This means that it is impossible to efficiently take into account the 21st harmonic value when measuring I0 and V0 in the network, which in some extreme situations may lead to disconnection of the healthy cable line in addition to the damaged cable line. As the fault resistance value increases, the share of higher harmonics decreases, and already, when the fault resistance value RF = 100 Ω, this share is much smaller than 50% of I1h.
It is noteworthy that the optimal sampling frequency for transient-based digital ground fault protection relays is recommended to be 100 kHz (or at least 20 kHz), as suggested in reference [25].
In the case of zero-sequence voltage V0 for a ground fault resistance value in the range from 1 to 80 Ω, the content of higher harmonics only once exceeds 10% of V0—namely for the 21st harmonic when the line-to-ground resistance value RF = 1 Ω. In the ground fault resistance range between 100 and 500 Ω, the proportion of higher harmonics increases sharply (between the second and tenth harmonics). The voltage magnitude spectrum envelope appears as a homographic function. For higher fault resistances, it can be assumed that the envelope of the zero-sequence voltage V0 spectrum has the shape of a hyperbola, but the inflection point of this hyperbolic envelope function occurs between the first and second harmonics, so the contribution of higher harmonics in this situation is not significant.

3.2. Steady-State Analysis of a Ground Fault in the Model Network Containing Higher Harmonics

This section of the paper is dedicated to the possible presence of higher harmonics in the network. It was assumed that, according to the standard [29], higher harmonics of voltages can occur in the industrial network. The maximum permissible contribution of higher harmonic voltage values was assumed in accordance with Table 1. The relevant standard imposes two conditions: the first states that none of the harmonics presented may, at any moment in time, exceed 5, 6, and 5% of the contribution of the first harmonic, respectively. The second condition implies that, in any 10 min interval, the total contribution of the higher harmonics must not exceed 8%. The second condition is put this way because in some cases (e.g., the start-up of a hoisting machine fed from a thyristor converter), the 8% condition in a 1 min interval is impossible to meet, but extending it to a 10 min interval is already within the norm. Therefore, in order to consider a worst-case scenario, it was considered that a ground fault occurred at a moment in a short interval when the contribution of higher harmonics was more than 9%, but in a 10 min interval, their contribution would be less than the required 8%. In mining facilities, thyristor-converter-powered hoisting machines are the dominant source of harmonic interference in the mine’s medium-voltage network [41,42]. They are equipped with high-power DC motors powered by controlled (usually 6-pulse or 12-pulse) thyristor rectifiers [43]. Machine speed control is achieved by varying the firing angle of the thyristors, according to a so-called travel diagram, which takes into account the allowable limits of acceleration and deceleration and the speed stabilization of the steady-state travel section. Hence, the greatest disturbances are generated during the machine’s starting and braking [44]. It is the relatively high motor power and the short, but frequently repeated, high-intensity load cycles that cause the hoisting machinery converter drives to have such a high impact on the operating conditions of the mine’s power distribution network [45,46]. Similar problems (however, without such an impact on the environmental safety of the network operation, as described in the introduction) exist in rectifier substations, supplying electrical traction with diverse and rapidly varying loads [47,48,49].
Ground fault simulations have been conducted for the same cases as outlined in Section 3.1, considering the influence of higher-order harmonics in the supply voltage. Example transient waveforms of zero-sequence currents and voltages are depicted (for the case of RF = 20 Ω) in Figure 15. Heavy waveform distortion is evident not only during the initial period (like in Figure 12) but even in the steady state. This phenomenon may disturb the correct operation of ground fault protection relays and should therefore be taken into account when selecting their settings. An overshoot of the criterion value in the first 20 ms may lead to nonselective tripping of the ground fault overcurrent protective devices, especially when their settings are selected on the basis of the analysis of the earth fault current determined in the steady state.
Comparisons of the rms values of selected parameters for the transient state (assumed as the first 20 ms of the fault duration) and the steady state are presented in Figure 16 and Figure 17. In order to ensure greater clarity, the relationship I0 = f(RF) is included in the charts below only for two cable lines outgoing from MS-1 panel 1 and DS-1 panel 1. These lines were chosen because they constitute a direct connection between the transformer supplying the network and the fault site. In addition, the relationship describing zero-sequence current and voltage components IF = f(RF) and V0 = f(RF) is provided. It can be seen that the rms value of the zero-sequence voltage V0 depends on the higher harmonics to an almost negligible extent. Therefore, it can be clearly stated that higher harmonics should not affect the operation of zero-sequence voltage protection. However, in the case of zero-sequence currents, it is slightly more important.
At low ground fault resistance values, individual currents in the network with a high content of higher harmonics can reach values up to 15% greater than in the situation without higher harmonics. When ground fault protection is set according to the I0 > criterion and a small value integration time is used, situations may arise in which the selectivity of ground fault clearing is not ensured—an undamaged cable can be switched off unnecessarily. In order to avoid such a situation, it is recommended in the literature to extend the measurement time over which the waveform will be integrated to values longer than 1 period (e.g., 40 or 60 ms) and to use a higher harmonics filter.
Both methods have their drawbacks. Increasing the integration time increases the tripping time of the protective relay. On the other hand, the higher harmonic filter does not take into account the fact that, in transient time, the distortions are nonperiodic waveforms; hence, the frequency spectrum is not discrete but continuous. In some extreme situations, even the application of the above solutions may still not be sufficient.
Spectral analysis of the same signals as in the previous graphs (Figure 18) for steady state confirms the observations made in those plots. Nevertheless, it is noteworthy to examine the influence of fault resistance on the higher harmonics’ magnitudes. For low-resistance fault values, a certain relationship can be observed, which is best described by Equation (45):
I 0 k h I 0 1 h V 0 k h V 0 1 h k
where k—harmonic order, I0kh—kth-order harmonic of the zero-sequence current at the given location, I01h—first-order harmonic of the zero-sequence current at the given location, V0kh—kth-order harmonic of the zero-sequence voltage in the galvanically connected network, V01h—first-order harmonic of the zero-sequence voltage in the galvanically connected network.
For example, it was assumed that the third harmonic of the phase voltage Vph3h before the occurrence of a ground fault was 5%. It turns out that for small ground fault resistance values, the third harmonic of the zero-sequence current (I03h) constituted about 15% of the first harmonic’s value regardless of the measurement point. As the value of RF increases, the percentage of higher harmonics in the fault current waveform decreases.
In order to highlight steady-state and transient spectral differences between the line-to-ground fault in the network with and without higher-order harmonics in the supply voltage, a differential scalogram has been calculated (difference between scalograms from Figure 14 and Figure 19), as depicted in Figure 20. It can be seen that the difference in the transient period is quite small and quickly decaying, but it still exists in the steady-state period and covers a fairly wide frequency range (several hundred Hz); however, the maximum value of this difference is several percent of the fundamental harmonic value. Thus, analysis of a model MV industrial network (designed to emulate a real coal mine network) has revealed that transients and the contribution of higher harmonics can be important factors influencing the operating conditions of ground fault protection relays in this network. A comparison of the zero-sequence current and voltage waveforms of the criterion values for the ground fault protection relays for the case of a network without higher-order harmonics and for a network with a significant proportion of harmonics in the supply voltage is shown in Figure 21.

4. Conclusions

The spectral (for steady state) and continuous wavelet (for transients) analyses indicate that increased levels of higher harmonics in the zero-sequence current waveform can cause nonselective operation of the ground fault protections, leading to unnecessary tripping of undamaged cable feeder lines and related severe disturbances in the plant operation. Due to the fact (described in Section 2) that, for zero-order sequences, the impedance of the insulated neutral point MV network is capacitive in nature (defined mainly by the values of the line-to-ground capacitances, with a relatively minor contribution from series resistances and inductances), the value of this impedance decreases significantly for higher harmonics of the network voltage. Owing to this phenomenon, even a minor contribution of higher harmonics to the line voltage can result in a large proportion of higher harmonics in the zero-sequence current measured and analyzed by the ground fault protection relays. This problem can become even more significant when there is a temporarily increased contribution of higher harmonics, e.g., in the frequently occurring dynamic states of high power drives supplied by thyristor converters. There has been also presented an innovative two-frequency analytical model of the transient component for a ground fault in an insulated neutral network. A novel derivation of the frequencies of the transient components and their decay coefficients has been described. This model is better suited to the specificities of mining industrial networks with long cable lines, shielded cables, and high-power and low-impedance power transformers. The problem of current waveform distortion poses a challenge to the functionality of ground fault protective relays based on both overcurrent and admittance criteria, so it would be advisable to increase the integration period for calculating the rms values or to incorporate harmonic filters into measurement circuits. The wavelet methods (particularly CWT used in our paper) are not applied directly to the design or setting of the relay protection, but they can help in understanding and analyzing the phenomena occurring during a ground fault and in the description and visualization of criterion signal waveforms.

Author Contributions

Conceptualization, K.K. and A.H.; methodology, K.K. and A.H.; software, A.H.; validation, K.K. and A.H.; formal analysis, A.H.; investigation, K.K.; resources, A.H.; data curation, A.H.; writing—original draft preparation, K.K.; writing—review and editing, A.H.; visualization, K.K. and A.H.; supervision, A.H.; project administration, K.K.; funding acquisition, K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by the Silesian University of Technology using resources from the Faculty of Mining, Safety Engineering, and Industrial Automation.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the large size of the data files.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Marinovic, N. Electrotechnology in Mining; Elsevier Science Publishers: Amsterdam, The Netherlands, 1990. [Google Scholar]
  2. Towsey, C. Electrical Safety in Mining—One flash and you’re ash. Australas. Mine Saf. J. 2011, 24, 42–47. [Google Scholar]
  3. Cawley, J.C. Electrical Accidents in the Mining Industry, 1990–1999. IEEE Trans. IA 2003, 39, 1570–1577. [Google Scholar] [CrossRef]
  4. Wang, D.; Ning, Y.; Zhang, C. An Effective Ground Fault Location Scheme Using Unsynchronized Data for Multi-Terminal Lines. Energies 2018, 11, 2957. [Google Scholar] [CrossRef]
  5. Azuara Grande, L.S.; Granizo, R.; Arnaltes, S. Wavelet Analysis to Detect Ground Faults in Electrical Power Systems with Full Penetration of Converter Interface Generation. Electronics 2023, 12, 1085. [Google Scholar] [CrossRef]
  6. Lee, K.M.; Park, C.W. Ground Fault Detection Using Hybrid Method in IT System LVDC Microgrid. Energies 2020, 13, 2606. [Google Scholar] [CrossRef]
  7. Li, B.; Du, X.; Miao, J.; Wang, H.; Ma, Y.; Li, Z. Arc Grounding Fault Monitoring and Fire Prediction Method Based on EEMD and Reconstruction. Electronics 2022, 11, 2159. [Google Scholar] [CrossRef]
  8. Pirmani, S.K.; Mahmud, M.A. Advances on fault detection techniques for resonant grounded power distribution networks in bushfire prone areas: Identification of faulty feeders, faulty phases, faulty sections, and fault locations. Electr. Power Syst. Res. 2023, 220, 109265. [Google Scholar] [CrossRef]
  9. Şolea, C.; Toader, D.; Vintan, M.; Greconici, M.; Vesa, D.; Tatai, I. The influence of insulation active power losses on the single line-to-ground fault in medium voltage electrical networks with resonant grounding. In Proceedings of the 10th International Conference on Modern Power Systems (MPS), Cluj-Napoca, Romania, 21–23 June 2023. [Google Scholar]
  10. Toader, D.; Greconici, M.; Vesa, D.; Vintan, M.; Şolea, 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]
  11. Yang, D.; Lu, B.; Lu, H. High-Resistance Grounding Fault Detection and Line Selection in Resonant Grounding Distribution Network. Electronics 2023, 12, 4066. [Google Scholar] [CrossRef]
  12. Toader, D.; Greconici, M.; Vesa, D.; Vintan, M.; Solea, C.; Maghet, A.; Tatai, I. The Influence of the Characteristics of the Medium Voltage Network on the Single Line-to-Ground Fault Current in the Resistor Grounded Neutral Networks. Designs 2021, 5, 53. [Google Scholar] [CrossRef]
  13. Xue, S.; Lian, J.; Qi, J.; Fan, B. Pole-to-Ground Fault Analysis and Fast Protection Scheme for HVDC Based on Overhead Transmission Lines. Energies 2017, 10, 1059. [Google Scholar] [CrossRef]
  14. Łowczowski, K.; Olejnik, B. Monitoring, Detection and Locating of Transient Earth Fault Using Zero-Sequence Current and Cable Screen Earthing Current in Medium Voltage Cable and Mixed Feeders. Energies 2022, 15, 1066. [Google Scholar] [CrossRef]
  15. Grakowski, Ł.; Gębczyk, K. Analysis of Earth Faults in the MV Grid Using the EMTP-ATP Program. Przegląd Elektrotechniczny 2019, 95, 146–149. [Google Scholar] [CrossRef]
  16. Doepken, H.C.; Schoolcraft, R.; Willis, O. Medium Voltage Cable Shielding and Grounding. In Proceedings of the IEEE Conference Record of Annual Pulp and Paper Industry Technical Conference, Atlanta, GA, USA, 19–23 June 1989; pp. 43–46. [Google Scholar]
  17. Haginomori, E.; Koshiduka, T.; Arai, J.; Ikeda, H. Power System Transient Analysis: Theory and Practice Using Simulation Programs (ATP-EMTP); John Wiley & Sons: Chichester, UK, 2008. [Google Scholar]
  18. Martinez Velasco, J.A. Transient Analysis of Power Systems: Solution Techniques, Tools, and Applications; John Wiley & Sons: Chichester, UK, 2015. [Google Scholar]
  19. Ali, S.A. Modeling of Power Networks by ATP-Draw for Harmonics Propagation Study. Trans. Electr. Electron. Mater. 2013, 13, 283–290. [Google Scholar] [CrossRef]
  20. Martinez, J.A. Power Quality Studies Using the ATP Package. Caribbean Colloquium on Power Quality (CCPQ). Dorado, 24–27 June 2003. Available online: https://inec.uprm.edu/ccpq-docs/Individual-Papers/j_martinez.pdf (accessed on 1 March 2024).
  21. Martinez, J.A. The ATP package. An environment for power quality analysis. In Proceedings of the International Conference on Harmonics and Quality of Power, Orlando, FL, USA, 1–4 August 2002. [Google Scholar]
  22. Granizo, R.; Blánquez, F.R.; Rebollo, E.; Platero, C.A. A Novel Ground Fault Non-Directional Selective Protection Method for Ungrounded Distribution Networks. Energies 2015, 8, 1291–1316. [Google Scholar] [CrossRef]
  23. Granizo, R.; Blánquez, F.R.; Rebollo, E.; Platero, C.A. New selective earth faults only current directional method for isolated neutral systems. In Proceedings of the Environment and Electrical Engineering (EEEIC), Venice, Italy, 18–25 May 2012; pp. 18–25. [Google Scholar]
  24. Ekici, S. Multi-class Support Vector Machines for Classification of Transmission Line Faults. Energy Educ. Sci. Technol. Part A Energy Sci. Res. 2012, 28, 1015–1026. [Google Scholar]
  25. Abdel-Fattah, M.L.; Lehtonen, M. Transient-Based Protection as a Solution for Earth-Fault Detection in Unearthed and Compensated Neutral Medium Voltage Distribution Networks. In Proceedings of the 2010 Electric Power Quality and Supply Reliability Conference, Kuressaare, Estonia, 16–18 June 2010. [Google Scholar]
  26. Liu, J.; Sun, K.; Ma, Z.; Liu, Z.; Li, K. Grounding Fault Model of Low Voltage Direct Current Supply and Utilization System for Analyzing the System Grounding Fault Characteristics. Symmetry 2021, 13, 1795. [Google Scholar] [CrossRef]
  27. Bender, C.M.; Orszag, S.A. Advanced Mathematical Methods for Scientists and Engineers. In Asymptotic Methods and Perturbation Theory; Springer: New York, NY, USA, 1999. [Google Scholar]
  28. Bush, A.W. Perturbation Methods for Engineers and Scientists; Taylor and Francis: Boca Raton, FL, USA, 1992. [Google Scholar]
  29. IEC Standard 61000-2-4:2002; Electromagnetic Compatibility (EMC)—Part 2–4: Environment—Compatibility Levels in Industrial Plants for Low-Frequency Conducted Disturbances. International Electrotechnical Commission: Geneva, Switzerland, 2002.
  30. Ram, B.; Vishwakarma, D.N. Power System Protection and Switchgear; McGraw Hill Education: London, UK, 2011. [Google Scholar]
  31. Kuliński, K.; Heyduk, A. Frequency Response Testing of Zero-Sequence Current Transformers for Mining Ground Fault Protection Relays. IJET 2020, 66, 701–705. [Google Scholar] [CrossRef]
  32. Ferreira, R.S.; Favoreto, H.H. Transient Overvoltages Due to Intermittent-Ground Faults in an Industrial Power System Grounded by a Resistance Connected to the Secondary of a Grounding Transformer. Electr. Power Syst. Res. 2023, 220, 109327. [Google Scholar] [CrossRef]
  33. Treider, T.; Høidalen, H.K. Estimating Distance to Transient and Restriking Earth Faults in High-Impedance Grounded, Ring-Operated Distribution Networks Using Current Ratios. Electr. Power Syst. Res. 2023, 224, 109765. [Google Scholar] [CrossRef]
  34. Addison, P.S. The Illustrated Wavelet Transform Handbook, 2nd ed.; Introductory Theory and Applications in Science, Engineering, Medicine and Finance; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  35. Addison, P.S. Introduction to redundancy rules: The continuous wavelet transform comes of age. Phil. Trans. R. Soc. A 2018, 376, 20170258. [Google Scholar] [CrossRef]
  36. He, Z. Wavelet Analysis and Transient Signal Processing Applications for Power Systems; John Wiley & Sons: Singapore, 2016. [Google Scholar]
  37. Olhede, S.C.; Walden, A.T. Generalized Morse Wavelets. IEEE Trans SP 2002, 50, 2661–2670. [Google Scholar] [CrossRef]
  38. Martinez-Riots, E.A.; Bustamante-Bello, R.; Navarro-Tuch, S. Applications of the Generalized Morse Wavelets: A Review. IEEE Access 2022, 11, 667–688. [Google Scholar] [CrossRef]
  39. Arfaoui, S.; Mabrouk, B.A.; Cattani, C. Wavelet Analysis. Basic Concepts and Applications; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  40. Misiti, M.; Misiti, Y.; Oppenheim, G.; Poggi, J.M. Wavelet ToolboxTM User’s Guide R2023b; The MathWorks Inc.: Natick, MA, USA, 2023. [Google Scholar]
  41. Ferreira, V.; Mendonça, G.; Rocha, A.; Resende, R.; Cardoso Filho, B. Medium voltage IGBT-based converters in mine hoist systems. In Proceedings of the IEEE Industry Applications Society Annual Meeting, Portland, OR, USA, 2–6 October 2016; pp. 1–8. [Google Scholar]
  42. Hyla, M. Higher harmonics filtration in the power supply system of thyristor hoisting machine of shaft transport in a mining plant. Prz. Elektrotechniczny 2022, 98, 43–48. [Google Scholar] [CrossRef]
  43. Siostrzonek, T.; Chmielowiec, K.; Piatek, K.; Dutka, M.; Firlit, A. The use of multi-pulse systems in the power supply of hoisting machine drives to improve voltage parameters in mining plants. In Proceedings of the 12th International Conference and Exhibition on Electrical Power Quality and Utilisation—(EPQU), Cracow, Poland, 14–15 September 2020. [Google Scholar]
  44. Mu Longhua, J. Measure and Suppression of Harmonic Currents for Heavy-duty Hoisting Machines. In Proceedings of the 2005 IEEE/PES Transmission & Distribution Conference & Exposition: Asia and Pacific, Dalian, China, 18 August 2005. [Google Scholar]
  45. Bebikov, Y.V.; Egorov, A.N.; Semenov, A.S. How Higher Harmonics Affect the Electrical Facilities in Mining Power Systems. In Proceedings of the 2020 International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), Sochi, Russia, 18–20 May 2020. [Google Scholar]
  46. Reshetnyak, S.; Bondarenko, A. Analysis of Technological Performance of the Extraction Area of the Coal Mine. In Proceedings of the E3S Web of Conferences, Semarang, Indonesia, 14–15 August 2018; Volume 41, p. 01014. [Google Scholar]
  47. Capasso, A.; Lamedica, R.; Manigrasso, R.; Sani, G. Reference Power Network for the Harmonic Propagation Analysis. Eur. Trans. Electr. Power 1992, 3, 167–188. [Google Scholar] [CrossRef]
  48. Kowalak, R.; Czapp, S.; Dobrzyński, K.; Klucznik, J.; Lubosny, Z. Harmonics produced by traction substations—Computer modelling and experimental verification. Przegląd Elektrotechniczny 2017, 93, 13–18. [Google Scholar] [CrossRef]
  49. Assefa, S.A.; Kebede, A.B.; Legese, D. Harmonic analysis of traction power supply system: Case study of Addis Ababa light rail transit. IET Electr. Syst. Transp. 2021, 11, 391–404. [Google Scholar] [CrossRef]
Figure 1. Simplified analysis of a single line-to-ground fault in a simple (not branched) network with an isolated neutral point: (a) network with series line impedances (ZL) and parallel insulation line-to-ground resistances (RG); (b) a simplified network with ZL and RG neglected; (c) internal impedance of the Thevenin generator; (d) straightforward equivalent circuit for the nth-order harmonic.
Figure 1. Simplified analysis of a single line-to-ground fault in a simple (not branched) network with an isolated neutral point: (a) network with series line impedances (ZL) and parallel insulation line-to-ground resistances (RG); (b) a simplified network with ZL and RG neglected; (c) internal impedance of the Thevenin generator; (d) straightforward equivalent circuit for the nth-order harmonic.
Energies 17 01532 g001
Figure 2. Ground fault current flow in a more complex network.
Figure 2. Ground fault current flow in a more complex network.
Energies 17 01532 g002
Figure 3. Simplified analysis of a single line-to-ground fault in a complex branched network with an isolated neutral point: (a) simplified network with ZL and RG neglected; (b) straightforward equivalent circuit for the nth-order harmonic.
Figure 3. Simplified analysis of a single line-to-ground fault in a complex branched network with an isolated neutral point: (a) simplified network with ZL and RG neglected; (b) straightforward equivalent circuit for the nth-order harmonic.
Energies 17 01532 g003
Figure 4. A toolchain of ATP/EMTP software used in our simulation research (*. means wildcards in the input and output filenames).
Figure 4. A toolchain of ATP/EMTP software used in our simulation research (*. means wildcards in the input and output filenames).
Energies 17 01532 g004
Figure 5. Layout of the industrial distribution network adopted for further analysis. (a) Model with separated capacitances of the faulty line and of the remaining part of the network; (b) simplified model with aggregated capacitances.
Figure 5. Layout of the industrial distribution network adopted for further analysis. (a) Model with separated capacitances of the faulty line and of the remaining part of the network; (b) simplified model with aggregated capacitances.
Energies 17 01532 g005
Figure 6. Equivalent circuit for the transient component (seen from the line-to-ground fault side) (a) network model from Figure 5. (bd) Successive simplification steps of the circuit.
Figure 6. Equivalent circuit for the transient component (seen from the line-to-ground fault side) (a) network model from Figure 5. (bd) Successive simplification steps of the circuit.
Energies 17 01532 g006
Figure 7. Comparison of the frequency and decay rate of the two transient components.
Figure 7. Comparison of the frequency and decay rate of the two transient components.
Energies 17 01532 g007
Figure 8. Layout of an example mining network analyzed in the paper. Green is the colour for the switchgear, red for the panels, red circle for the zero-sequence current i0(t).
Figure 8. Layout of an example mining network analyzed in the paper. Green is the colour for the switchgear, red for the panels, red circle for the zero-sequence current i0(t).
Energies 17 01532 g008
Figure 9. The MV network diagram used during the simulation, created in the ATPDraw graphical preprocessor as input data for the ATP package.
Figure 9. The MV network diagram used during the simulation, created in the ATPDraw graphical preprocessor as input data for the ATP package.
Energies 17 01532 g009
Figure 10. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF from the range 1 Ω … 10 kΩ for various fault locations: (a) relative value of zero-sequence voltage component V0/VPh, (b) magnitude of zero-sequence current component I0 for MS-1 panel 1, (c) magnitude of zero-sequence current component I0 for MS-1 panel 2, (d) magnitude of zero-sequence current component I0 for DS-1 panel 1, (e) magnitude of zero-sequence current component I0 for DS-1 panel 2, (f) magnitude of zero-sequence current component I0 for DS-1 panel 3. Green graph—for steady-state values, blue graph—for Δt = 20 ms and for fault location distance from departmental switchgear DS-1 x = 0 m, red graph—for Δt = 20 ms and for fault location distance from departmental switchgear DS-1 x = 1500 m.
Figure 10. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF from the range 1 Ω … 10 kΩ for various fault locations: (a) relative value of zero-sequence voltage component V0/VPh, (b) magnitude of zero-sequence current component I0 for MS-1 panel 1, (c) magnitude of zero-sequence current component I0 for MS-1 panel 2, (d) magnitude of zero-sequence current component I0 for DS-1 panel 1, (e) magnitude of zero-sequence current component I0 for DS-1 panel 2, (f) magnitude of zero-sequence current component I0 for DS-1 panel 3. Green graph—for steady-state values, blue graph—for Δt = 20 ms and for fault location distance from departmental switchgear DS-1 x = 0 m, red graph—for Δt = 20 ms and for fault location distance from departmental switchgear DS-1 x = 1500 m.
Energies 17 01532 g010
Figure 11. Zero-sequence voltage and current components as functions of line-to-ground fault resistance RF for various fault locations: (a) relative value of zero-sequence voltage component V0/VPh, (b) magnitude of zero-sequence current component I0 for MS-1 panel 1, (c) magnitude of zero-sequence current component I0 for MS-1 panel 2, (d) magnitude of zero-sequence current component I0 for DS-1 panel 1, (e) magnitude of zero-sequence current component I0 for DS-1 panel 2, (f) magnitude of zero-sequence current component I0 for DS-1 panel 3. Green curve—for steady state, blue color—for Δt = 80 ms and x = 0 m, red color—for Δt = 80 ms and x = 1500 m.
Figure 11. Zero-sequence voltage and current components as functions of line-to-ground fault resistance RF for various fault locations: (a) relative value of zero-sequence voltage component V0/VPh, (b) magnitude of zero-sequence current component I0 for MS-1 panel 1, (c) magnitude of zero-sequence current component I0 for MS-1 panel 2, (d) magnitude of zero-sequence current component I0 for DS-1 panel 1, (e) magnitude of zero-sequence current component I0 for DS-1 panel 2, (f) magnitude of zero-sequence current component I0 for DS-1 panel 3. Green curve—for steady state, blue color—for Δt = 80 ms and x = 0 m, red color—for Δt = 80 ms and x = 1500 m.
Energies 17 01532 g011
Figure 12. Example transient waveforms in the case of 20 Ω line-to-ground resistance value: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, (d) fault current.
Figure 12. Example transient waveforms in the case of 20 Ω line-to-ground resistance value: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, (d) fault current.
Energies 17 01532 g012
Figure 13. Percentage of higher harmonics relative to the first harmonic magnitude for selected values of line-to-ground resistance. Spectra shown: (a) zero-sequence current I0 for main switchgear MS-1 panel 1, (b) zero-sequence current I0 for departmental switchgear DS-1 panel 1, (c) zero-sequence voltage V0, and (d) fault current IF.
Figure 13. Percentage of higher harmonics relative to the first harmonic magnitude for selected values of line-to-ground resistance. Spectra shown: (a) zero-sequence current I0 for main switchgear MS-1 panel 1, (b) zero-sequence current I0 for departmental switchgear DS-1 panel 1, (c) zero-sequence voltage V0, and (d) fault current IF.
Energies 17 01532 g013
Figure 14. CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network without higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Figure 14. CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network without higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Energies 17 01532 g014
Figure 15. Example transient waveforms in the case of 20 Ω line-to-ground resistance value: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, and (d) fault current.
Figure 15. Example transient waveforms in the case of 20 Ω line-to-ground resistance value: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, and (d) fault current.
Energies 17 01532 g015
Figure 16. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF in the range 1 Ω … 10 kΩ. (a) Zero-sequence current component I0 for main switchgear MS-1 panel 1, (b) zero-sequence current component I0 for departmental switchgear DS-1 panel 1, (c) fault current IF, (d) zero-sequence voltage component V0. All graphs shown for Δt = 20 ms and fault distance x = 600 m. Blue graph—for network conditions without higher harmonics, red graph—for network conditions with higher harmonics.
Figure 16. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF in the range 1 Ω … 10 kΩ. (a) Zero-sequence current component I0 for main switchgear MS-1 panel 1, (b) zero-sequence current component I0 for departmental switchgear DS-1 panel 1, (c) fault current IF, (d) zero-sequence voltage component V0. All graphs shown for Δt = 20 ms and fault distance x = 600 m. Blue graph—for network conditions without higher harmonics, red graph—for network conditions with higher harmonics.
Energies 17 01532 g016
Figure 17. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF in the range 1 Ω … 10 kΩ. (a) Zero-sequence current component I0 for main switchgear MS-1 panel 1, (b) zero-sequence current component I0 for departmental switchgear DS-1 panel 1, (c) fault current IF, (d) zero-sequence voltage component V0. All graphs shown for a steady state and a fault distance x = 600 m. Blue graph—for network conditions without higher harmonics, red graph—for network conditions with higher harmonics.
Figure 17. Zero-sequence current and voltage components as functions of line-to-ground fault resistance RF in the range 1 Ω … 10 kΩ. (a) Zero-sequence current component I0 for main switchgear MS-1 panel 1, (b) zero-sequence current component I0 for departmental switchgear DS-1 panel 1, (c) fault current IF, (d) zero-sequence voltage component V0. All graphs shown for a steady state and a fault distance x = 600 m. Blue graph—for network conditions without higher harmonics, red graph—for network conditions with higher harmonics.
Energies 17 01532 g017
Figure 18. Percentage of higher harmonics relative to the first harmonic for selected values of the short circuit resistance calculated in steady state for a model network with higher harmonics. Spectra shown: (a) zero-sequence current I0 for main switchgear MS-1 panel 1, (b) zero-sequence current I0 for departmental switchgear DS-1 panel 1, (c) zero-sequence voltage V0, and (d) fault current IF.
Figure 18. Percentage of higher harmonics relative to the first harmonic for selected values of the short circuit resistance calculated in steady state for a model network with higher harmonics. Spectra shown: (a) zero-sequence current I0 for main switchgear MS-1 panel 1, (b) zero-sequence current I0 for departmental switchgear DS-1 panel 1, (c) zero-sequence voltage V0, and (d) fault current IF.
Energies 17 01532 g018
Figure 19. CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network with higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Figure 19. CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network with higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Energies 17 01532 g019
Figure 20. Differential CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network with and without higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Figure 20. Differential CWT scalogram of the zero-sequence filter output (IF = 3I0) current for the network with and without higher-order harmonics in source voltage: (a) 2D view, (b) 3D view.
Energies 17 01532 g020
Figure 21. Comparison of transient waveforms in the case of 20 Ω line-to-ground resistance value with and without influence of higher-order harmonics: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, and (d) fault current.
Figure 21. Comparison of transient waveforms in the case of 20 Ω line-to-ground resistance value with and without influence of higher-order harmonics: (a) zero-sequence current for main switchgear, (b) zero-sequence current for departmental switchgear, (c) zero-sequence voltage, and (d) fault current.
Energies 17 01532 g021
Table 1. Percentage of higher harmonic components adopted for the simulation analysis of ground faults in an example mining network.
Table 1. Percentage of higher harmonic components adopted for the simulation analysis of ground faults in an example mining network.
Harmonics Order hHarmonics Value Relative to VPh1h, %
35
56
75
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kuliński, K.; Heyduk, A. Ground Fault in Medium-Voltage Power Networks with an Isolated Neutral Point: Spectral and Wavelet Analysis of Selected Cases in an Example Industrial Network Modeled in the ATP-EMTP Package. Energies 2024, 17, 1532. https://doi.org/10.3390/en17071532

AMA Style

Kuliński K, Heyduk A. Ground Fault in Medium-Voltage Power Networks with an Isolated Neutral Point: Spectral and Wavelet Analysis of Selected Cases in an Example Industrial Network Modeled in the ATP-EMTP Package. Energies. 2024; 17(7):1532. https://doi.org/10.3390/en17071532

Chicago/Turabian Style

Kuliński, Krzysztof, and Adam Heyduk. 2024. "Ground Fault in Medium-Voltage Power Networks with an Isolated Neutral Point: Spectral and Wavelet Analysis of Selected Cases in an Example Industrial Network Modeled in the ATP-EMTP Package" Energies 17, no. 7: 1532. https://doi.org/10.3390/en17071532

APA Style

Kuliński, K., & Heyduk, A. (2024). Ground Fault in Medium-Voltage Power Networks with an Isolated Neutral Point: Spectral and Wavelet Analysis of Selected Cases in an Example Industrial Network Modeled in the ATP-EMTP Package. Energies, 17(7), 1532. https://doi.org/10.3390/en17071532

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