Next Article in Journal
Cost-Risk Analysis Reconsidered—Value of Information on the Climate Sensitivity in the Integrated Assessment Model PRICE
Next Article in Special Issue
Numerical Simulation for Redundant Electro-Hydrostatic Servo-Actuators under Certain Special Conditions
Previous Article in Journal
Study on the Positioning Accuracy of the GNSS/INS System Supported by the RTK Receiver for Railway Measurements
Previous Article in Special Issue
An Impedance Matching Solution to Increase the Harvested Power and Efficiency of Nonlinear Piezoelectric Energy Harvesters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Artificial Neural Networks and Nonlinear Optimization Techniques to Fault Location in Transmission Lines—Statistical Analysis

by
Simone A. Rocha
1,*,
Thiago G. Mattos
1,
Rodrigo T. N. Cardoso
1 and
Eduardo G. Silveira
2
1
Program in Mathematical and Computational Modeling, Federal Center for Technological Education of Minas Gerais, Belo Horizonte 30180-001, Brazil
2
Department of Electrical Engineering, Federal Center for Technological Education of Minas Gerais, Belo Horizonte 30180-001, Brazil
*
Author to whom correspondence should be addressed.
Energies 2022, 15(11), 4095; https://doi.org/10.3390/en15114095
Submission received: 2 April 2022 / Revised: 21 May 2022 / Accepted: 24 May 2022 / Published: 2 June 2022

Abstract

:
This study presents applications of artificial neural networks and nonlinear optimization techniques for fault location in transmission lines using simulated data in an electromagnetic transient program and actual data occurring in transmission lines. The localization is performed by a modular structure of 4 neural networks and by the minimization of objective functions descriptive of the problem, defined according to the parameters of the line and the type of short circuit, submitted to the methods Quasi-Newton, Ellipsoidal, and Real Polarized Genetic Algorithm. The results obtained are compared statistically with those of a classical analytical method. The analysis of the variance of location errors presented by the methods revealed, with 5% significance, statistical evidence that allowed the conclusion that the type of method used affects fault location indication. In simulated scenarios, minor errors were obtained with the neural network and larger with the analytical method. For field oscillographic, the largest errors were in the neural network; there is no evidence to reject the equality between the results of the analytical method and the nonlinear optimization techniques. The Tukey test identified no differences between the nonlinear optimization methods applied to the proposed objective functions, but the low computational cost associated with the Quasi-newton method highlights it. The nonlinear optimization methods used for the localization function proved to be promising for application in companies that operate electrical systems, providing localization errors similar to those presented by the classical analytical method.

1. Introduction

Overhead transmission lines (TL) are integral components of an electric power system that enable the supply of electric energy from its generation to the distribution networks. Generally, a short circuit (fault) in a transmission line represents a phenomenon that is difficult to predict, and such faults are characterized by the instant of occurrence, the type classification: phase-to-ground (SPG), two-phase (DP), two-phase-to-ground (DPG), or three-phase (TP), the location indication d F , and the fault resistance value ( R F ) . In the event of a fault, the protection relay detects, identifies, and signals the event, commanding the circuit breakers to remove the short-circuited line from service. Following its actuation, automatic restart attempts are made. If successful, the line is reintegrated into the system; however, if the attempts fail, the line needs to be shut down until repair work is performed, which should ideally take place in the shortest time possible and with adequate reliability. Therefore, indicating the point of occurrence of the defect with the most minor possible error reduces the costs related to the shutdown.
Since the 1970s, there has been ongoing research about locating faults in transmission lines. However, it remains a considerable problem. The development of hardware has allowed access to more measurement points with time synchronization via the global positioning system (GPS), creating the possibility of applying new mathematical techniques to problems involving the electric power system. In addition to avoiding the payment of hefty fines to which electricity concessionaires are subject to regulatory agents, locating faults rapidly and accurately directly reflects the quality of the electricity supply and the stability of the system. To address the problem, according to the proposed location method, scholars on the subject establish, as shown in Table 1, the location scenario, the origin and type of signal used, and the characteristics related to the pre-processing of the data. Existing methods typically use phasors of the voltages and currents of the transmission line terminals. Some authors applied the wavelet transform to the sampled signals to extract the detail coefficients and then the energy or frequencies.
Several authors have addressed the issue of fault location in TL through artificial neural networks (ANN). In [21], the authors compare these structures to a black box. According to [7], there is no single neural network structure and the results obtained with the technique highlight the possibilities presented by researchers, a situation verified through the bibliographic survey, as shown in Table 2. Multilayer perceptron feedforward networks are predominantly used, with supervised learning backpropagation and the Levenberg Marquardt training algorithm. In most existing studies, with variations among the authors, the training data were obtained by simulating faults along the transmission line, with different values of R F , fault incidence angle, and equivalent impedance of the sources generated through electromagnetic transient programs. Modular networks, such as those used by [2,5,8,9,12,14,16,17,26,29,30,31], represent a possibility in solving the problem of locating different types of faults in LTs, contributing to the reduction of computational effort.
Considering some recent work, in the study presented in [23], instantaneous fault current measurements obtained by simulation are used in modular networks with training guided by Bayesian regularization. Most of these tests in that study indicate the short circuit distance with errors below 0.5% in relation to the simulated location. In [25], location errors in the order of 10−5 were achieved by employing the magnitudes of the frequencies of simulated fault voltage and current waves to feed the inputs of a single ANN. Moreover, simulated voltage and current signals from one of the line terminals of the electrical system in Great Britain were used in [26] as input in modular networks, obtaining errors of less than 0.7% in relation to the location of the short circuit. Current and voltage wave frequencies from one of the terminals of a 230 kV and 100 km long transmission line, categorized into 5 groups, were used in [27] as inputs to an ANN for the location of simulated SPG type faults, allowing to achieve errors of less than 0.4% in relation to the location of the fault. From the survey of academic production, there is a predominance of the use of ANN in simulated situations, and the results of an application of ANN to real cases of faults in TL are presented only in [30], where reduced errors were achieved in a transmission line of the electricity supply system in Iran. The situation is related to the difficulty of generalizing ANN, which provides reduced errors in simulated fault location cases, but presents a lack of performance when applied to real cases, indicating the difficulty of making the training data represent the systems of the concessionaires assertively.
Regarding the use of optimization techniques, the Simplex method of linear optimization, Nelder-Mead, was applied by [6] to an objective function to obtain the distance and fault resistance of simulated cases. The Nelder-Mead and Broyden–Fletcher–Goldfarb–Shanno numerical optimization methods were proposed in [10] to locate simulated faults in transmission lines. Harmony Search optimization was presented by [15] to estimate the distance from the substation to the fault point in simulated cases. Different optimization techniques are applied by [21] to estimate the fault location from the minimization of an objective function of a variable. In the authors’ analysis, although small localization errors were achieved in all situations, the Teaching Learning Based Optimization technique presented a shorter convergence time to obtain the results. A hybrid model that combines the Relief algorithm (which performs the sequential comparison of the entire database, returning a logical response associated with the operating functionality parameters) with the Transformed Wavelet for fault detection and location is proposed [22]. The implementation of genetic algorithms to the problem of fault location in TL is presented in [28] to compare waveforms from digital fault recorders and simulated waves. The phasors obtained during the fault were used to calculate the fitness value. The results of the proposal were not satisfactory when applied to real faults. In [31], the authors proposed objective functions based on symmetrical components that, applied to simulated and real faults in the Brazilian electrical system and minimized by the ellipsoidal algorithm, provided promising results for using the method in practical situations. In the context presented, this paper aims to:
  • to analyze the use of ANN and nonlinear optimization techniques (NLO) in the location of TL faults for simulated and real cases of the Brazilian electric power transmission system;
  • to compare the results obtained by the analyzed methods with those provided by the classical analytical method (AM), proposed in [32] and used by some electric utilities in practical applications;
  • apply a statistical analysis that allows evaluating the existence of differences between the responses obtained by the methods in the fault location process.
The research proposal is related to evaluating and comparing the errors achieved by the implemented methods so that the possibility of joint use can be verified to improve the fault location process in TL. The application of more than one method can give engineering and maintenance teams confidence in indicating the distance to failure. In practice, inaccurate results drive repair personnel away from the point of failure and reduce confidence in the process. The proposed methods, developed and applied in simulated and real situations of faults in the Brazilian electrical system, allow statistical analysis of the variance and application of the Tukey’s test to localization errors, validating the difference between the results achieved by the different methods. Also, considering the bibliographic contribution used, the present study made it possible to verify that:
  • in simulated scenarios, smaller localization errors are obtained using ANN and larger using the AM, leaving the NLO techniques with intermediate errors. In real scenarios, different from the situation verified in simulations, greater errors were evidenced with the use of ANN, with no statistical evidence to reject the equality between the fault location performed by the AM and the NLO methods;
  • the neural networks implemented in the proposed paper, based on techniques found in the literature, proved to be incompatible for application in companies that operate electrical systems, providing location errors significantly greater than the other methods, far from acceptable in practical situations (up to 5%);
  • the Quasi-Newton (QN), Ellipsoidal (EL), and genetic (PRGA) of NLO methods were used in the study. Among them, considering the precision and computational effort, the QN method was the most suitable to be used in field applications, being able to be used together with the AM for the defect location indication.
Notably, for the location methods presented, this study does not intend to evaluate the sensitivity of the algorithms to factors such as current and potential transformer errors, line model, errors inline parameters, synchronization, mutual inductance, transposition, and phasor estimation method. Data are presented in equal conditions to all methods. Accordingly, computer programs containing routines and mathematical techniques necessary to provide the distance to the fault in relation to one of the line terminals were implemented.
Although the paper deals with fault location methods, preventive measures such as transient earth voltage measurement to detect partial discharges [33,34], methods to detect faulty insulators [35,36], AM modulation methods to detect sources of problems in voltage parameters [37,38] and thermography [39] can be adopted to reduce or avoid future problems.

2. Materials and Methods

2.1. Simulated Test System

Consider the three-line model of Figure 1, which represents the electrical system of the faulty transmission line proposed for simulation, where L is the extension of the transmission line, and dF is the distance from the fault point to Terminal S, which occurs through a fault resistance (RF).
The process begins by defining the electrical system used and simulating short circuit cases. A pilot project was developed for a real transmission line of a Brazilian utility that is 74.4 km in length with a capacity of 345 kV. Resistance r and reactance x data from the terminal equivalent sources and the line are considered to be ideally transposed along its length with lumped parameters, as listed in Table 3.
The transmission line parameters were inserted into the Alternative Transients Program (ATP) [40] to obtain simulated data, as shown in Figure 2, where the voltmeter and ammeter represent the protection relay or digital disturbance recorder installed at each line terminal. The simulator contained models of the electrical components of the electrical power system and was used because it is a public domain software. After the simulation, ATP exported a .txt file containing voltage and current data to be used in the various steps that make up fault classification and location programs. The electric energy concessionaire provided only the parameters R, L, and C of the transmission lines used in real case studies. The authors chose to represent the transmission line by resistance and inductance to carry out the studies. More complete models for transmission lines can be used, such as the one by JMARTI [41] and LMARTI [42], as long as the line geometry and cable data are available.
From the model, variations in the short circuit location, the RF value, and the different types of faults were simulated represented in Figure 3: AG—phase A to ground, BG—phase B to ground, CG—phase C to ground, ABG—phase A to phase B to ground, BCG—phase B to phase C to ground, AB—phase A to phase B, BC—phase B to phase C, AC—phase A to phase C and ABC—phase A to phase B to phase C.
As specified in Table 4 and Table 5, the procedure resulted in the generation of 1368 failure scenarios for the composition of the ANN training database and a further 336 cases for the validation of the ANN and application of the three NLO techniques and the AM method. As indicated in [43], the number of training and validation data respectively corresponds to 80% and 20% of the data set. In the case of faults involving the ground, according to references [44,45], higher values of fault resistance can be used to generate data for training neural networks.

2.2. Process Steps

The necessary procedures involve the steps indicated in the flowchart presented in Figure 4. Data adequacy begins with reading both line terminals’ voltage and current inputs. Subsequently, fault detection occurs, which entails the identification of the short circuit start instant in the voltage and current waves, thereby allowing the separation of the data corresponding to the pre-fault and fault periods.
In the pre-processing of voltage and current signals, necessary for data preparation prior to fault location, higher frequencies were removed by digital low-pass filtering Butterworth with a cut-off frequency of 100 Hz. After that, the sampling frequency was reduced to the desired value of 16 points per cycle of the fundamental frequency (60 Hz in Brazil). The voltages V A S ,   V B S ,   V C S ,   V A R ,   V B R ,   V C R and currents I A S ,   I B S ,   I C S ,   I A R ,   I B R ,   I C R phasors associated with the fundamental frequency were estimated through the do least squares method, as described in [46]. Figure 5 shows the voltage and current waves for an AG fault at the sender terminal of the transmission line, which started at 0.049 s. Figure 6 shows the filtered waves, and Figure 7 shows the estimated phasor modules. Figure 8 demonstrates the process of selecting current phasors in phase A. Finally, the fault classification, necessary for the NLO methods, was performed, and then the location process.
Fault location occurs after the fault detection and classification steps. In this study, it was estimated through ANN and NLO techniques applied to the objective functions adapted from [31], as well as the analytical method proposed in [32], implemented through MATLAB® software.

2.3. ANN for Fault Location

For fault location, ANN was implemented using the neural network toolbox in MATLAB®, as described by [47]. In structuring the networks, the training was performed separately according to the type of short circuit identified by the classification. Feedforward networks were used with supervised learning backpropagation and the Levenberg-Marquardt training algorithm hyperbolic tangent activation function. In the feedforward model, the network was arranged in layers, and the propagation of information followed from the input to the output without feedback from previous units. The supervised learning backpropagation and Levenberg Marquardt training algorithm were used to minimize the error obtained by the difference between the response in the output layer and the desired value from the correction of the weights in all layers, starting from output to input. The organization of neurons and layers was defined experimentally, with 2 hidden layers of 30 and 20 neurons. After pre-processing the data simulated by ATP, the voltage (VAS, VBS, VCS, VAR, VBR, VCR) and current (IAS, IBS, ICS, IAR, IBR, ICR) phasors amplitude from the local and remote terminals of the line were used as the input data for training the ANN, with the function of providing the estimated fault location as an output, according to the selected modular structure, as indicated in Figure 9.
The voltage and current phasors from each validation file in Table 5 were inserted into the trained networks to verify the generalization and adaptation capacity of the location ANN, returning the distance of fault.

2.4. Nonlinear Optimization for Fault Location

The NLO problem used in this study involves minimizing a nonlinear objective function f : D f n     , overall vectors x D f . Depending on the logic employed to determine the solution to the problem, different minimization techniques can be applied to the function f to provide a particular solution x * D f n such that   x D f , f x * f x . In this context, the fault location problem is related to the minimization of one of the mono-objective functions, F d F , R F , indicated in Table 6. These functions, adapted from [31], contain fault location and fault resistance as variables, which are based on the currents and voltages of the TL terminals and depend on the types of existing faults. Therefore, the variables must be included in minimization algorithms. The subscripts 0 and 1 indicate, respectively, the values of zero and positive sequence, Z denotes the impedance of the line in Ω / km , I S ,   I R ,   V S and V R denote the respective currents and voltages at the local (S) and remote (R) terminals.
The graphical and analytical analyses of the objective functions that describe the process make it possible to place them as differentiable, convex, and unimodal representations of two real variables, represented by F d F ,   R F . As an example of the characteristics of the function, the graph of a simulated AG fault event occurring 45 km away with R F of 20 Ω is illustrated in Figure 10. NLO methods based on the search direction and the exclusion of regions and population dynamics were applied to minimize the objective functions. The NLO techniques are primarily distinguished by the logic used in the search for the solution to the problem.

2.4.1. Search Direction Methods

The Quasi-Newton search direction method used in the research indicates the minimum of the objective function. Information from previous searches is considered until the algorithm satisfies some previously established stopping criterion. From the gradient information, approximations are generated recursively to estimate the inverse of the hessian of the objective function. There are different strategies for this approximation, and, among the methods described by [48], the combination between the DFP (Davidson-Fletcher-Powell) and BFGS (Broyden-Fletcher-Goldfarb-Shanno) methods was used, as indicated in lines 14, 15, and 16 of the algorithm presented in Algorithm 1. The vector r represents the difference between the gradient vector of the objective function in the current position and the previous position; the vector v is the difference between the current position and the previous one, the coefficient γ indicates the weight of DFP and BFGS, and H represents the approximation to the hessian of the objective function.
Algorithm 1 Quasi_Newton Algorithm (y,x)
Energies 15 04095 i001
The search direction methods perform local searches and are strongly dependent on the derivative, the hessian of the objective function, or both at the points obtained at each iteration. Consequently, it may present difficulties for their use. When applied to convex functions, converge to the local minimum, which in these cases also represents the global minimum, a condition verified for the adapted objective functions of [32]. However, if the function were not convex, the local minimum cannot be considered the global minimum. If the objective function were multi-modal, there would be several basins of attraction. Thus, the problem of determining the global minimum would be complex using only local information. This is because of the inability to recognize a global minimum, even within a set of well-defined local minima. The graphic representation shown in Figure 11 exemplifies the sequence of points generated by the QN method for a simulated AG fault at 29.76 km from the transmission line’s emitting terminal, with an R F of 20 Ω.

2.4.2. Region Exclusion Methods

The ellipsoidal method is a region exclusion method. Accordingly, in a convex set with a non-empty interior, the minimum can be found in a finite number of steps through projections in half-spaces and investigation of the subgradients of the objective function. Intuitively, as described in [49], from an initial ellipsoid, which contains the point to be reached, cuts were made to always generate ellipsoids of decreasing volumes for a fixed ratio, which only depends on the dimension of space. The process is outlined by the algorithm shown in Algorithm 2.
Algorithm 2 Ellipsoidal Algorithm (x,Q)
Energies 15 04095 i002
In the algorithm, n denotes the dimension of the initial point,   Q represents the matrix of the ellipse that surrounds the optimal point to be obtained at each iteration, and the vector m k   is the subgradient of the objective function at a given point belonging to the domain of the function. A sequence of points is generated in which x k represents the center of each ellipsoid. After a finite number of iterations, the smallest ellipsoid over the sought point x * D f n is obtained such that   x D f , f x * f x is given by the center of the ellipsoid. In the Ellipsoidal algorithm, the graph in Figure 12 represents the sequence of points generated by the EL method for the same fault simulated in Figure 11. At each iteration, the EL method returns the center of the ellipsoid generated in the process. The exclusion of regions occurs from the subgradient without any direct dependence on the derivatives of the objective function. However, the method requires the objective function to be convex or quasi-convex and unimodal, which is also verified with the adapted functions in [31]. If these conditions are not met, it will lead the algorithm to an unpredictable situation.

2.4.3. Population Methods

The search with population methods seeks to reach optimal solutions from a set of solutions. The minimization of a function using a genetic algorithm represented schematically in Algorithm 3 commences from the random generation of a set of points (individuals) that form the initial population studied, representing possible solutions for a given problem. Crossover and mutation genetic operators are applied sequentially during the iterative process, followed by evaluation and selection processes to generate the surviving population. The evaluation of each results in a value called fitness. The better the fitness value, the greater the individual’s chance to remain in the surviving population.
Algorithm 3 Genetic Algorithm (x)
Energies 15 04095 i003
For the NLO desired in this study, the best fitness is associated with the set of solutions ( d F , R F ) that minimize f : D f n . The population evolution along the iterations occurs according to an established stopping criterion that leads the set of individuals to a region close to the optimal point. Population methods do not use information involving derivatives of the objective function for convergence. Instead, they operate with parameter encoding and indicate optimal solutions from a set of solutions. Owing to its non-deterministic character, the process is indicated for the treatment of functions and constraints that are difficult to model, requiring statistical analysis to study their convergence. In this study, we employed the polarized real genetic algorithm described in [50]. In PRGA, each optimization parameter is represented by a real variable, and the set of parameters is stored in a vector that represents an individual. For the simulated AG fault located 29.76 km from the emitting terminal of the transmission line, with an R F of 20 Ω, the graph presented in Figure 13 refers to the last generation of the PRGA execution.
The previously described NLO algorithms were applied to the minimization of the objective functions adapted from [31]. Based on the type of fault, the estimated distance to the point of occurrence of the short circuit in relation to the line terminal was obtained for each of the scenarios specified in Table 5. A variation of the analytical method, which disregards the capacitance of the TL, proposed in [32] was also implemented and applied to the validation cases for comparisons of its results with those provided by the ANN and NLO methods. In this method, fault location is obtained from descriptive functions of the electrical circuit, capable of representing the voltage and current variations along the transmission line from its two terminals.

3. Statistical Study

An analysis of possible statistical differences between the proposed location methods and the types of fault that can occur in a fault event was developed by analyzing the variance with two factors (Anova2) of the location errors presented by the methods using Python [51]. The procedure, as described in [52], allows us to investigate the equality of means in experiments with more than one factor. For the application of the process, the test data must be obtained randomly, and the data must be adjusted to a normal distribution, with mean μ = 0 and constant variance σ2. These requirements, used in the process of formulating the Anova2 model, allow the generalization of a population from a sample. Regarding the test, it is possible to verify the existence of evidence of variations between the levels of factors that may interfere in the process, signaling the existence of at least one level that differs from the others. In applying the method to the fault location problem, the type of method and the type of fault that occurred were considered factors to be analyzed. Internally, different levels were considered for each of these factors, as shown in Figure 14.
Mathematically, Anova2 can be expressed using Equation (1), and it is suitable for the problem to be studied. Accordingly, µ is considered as the global average effect of errors arising from the experimental location process, τ i denotes the effect of the ith level of the fault type factor, β j represents the effect of the jth level of the method type factor, τ β i j indicates the effect of the interaction between the fault type and method type factors, i is the number of levels of the fault type factor,   j indicates the number of levels of the method type factor, and k denotes the number of samples collected for each level.
y i j k = μ + τ i + β j + τ β i j + ε i j k   i = 1 ,   , a   j = 1 ,   ,   b k = 1 ,   2 ,   ,   k
The test deals with the influence of the interaction between the two factors analyzed and occurred on the hypothesis tests described in Equations (2)–(4). Each null hypothesis ( H 0 ) supports the assumption that the averages of the errors of the levels are the same for the analyzed factor. Moreover, with the alternative hypothesis ( H 1 ), the existence of at least one different average is assumed. The test results conducted and interpreted in [52] provide tools that allow us to either accept or reject a statistical hypothesis through the evidence provided by the sample.
H 0   : τ 1 = = τ i = 0 ;   H 1   :   i /   τ i 0
H 0   : β 1 = = β j = 0 ;   H 1   :   j /   τ j 0
H 0   : τ β i j = 0   i , j ;   H 1   :   i ,   j /   τ β i j 0  
When the Anova2 results reject H 0 and, consequently, accept H 1 , there is statistical evidence that at least one of the means of the levels differs from the others for the analyzed factor. This conclusion, analyzed for the fault location problem, allows us to verify whether the type of fault that occurred or the type of method used leads to a differentiation in the efficiency of the location process. The rejection of H 0 is performed by considering a certain level of significance related to the probability of rejecting this hypothesis when it is true. In statistical hypothesis tests, the rejection of H 0 is considered to be significant when the observed p-value is lower than the significance level defined for the study.
Although the Anova2 test indicates the existence of a discrepant level among the others, it does not locate the observed difference. It is necessary to apply a test of comparison of means between the levels to determine the best (or worst) level of the factor. As a complement to the Anova2 study, the Tukey test was applied at a significance level of 5%. The test involves building confidence intervals for all pairs of means in such a manner that the set of all intervals has a certain degree of confidence.

4. Results

In the fault location process, deterministic (QN, EL, AM) and non-deterministic (ANN and PRGA) procedures were implemented and applied. In the application of ANN to the problem, intelligent processes are used to recognize patterns simulated by the ATP related to the identification of the fault location. These structures characterized by adaptation by experiences, learning ability, generalization ability, data organization, etc., can provide a satisfactory solution to the problem. We attempted to obtain the minimum point of an objective function of two variables by applying the implemented NLO algorithms, f d F , R F , from the methods of QN, EL, and PRGA. Table 7 shows the number of iterations performed and the execution time required by each algorithm to reach the result referring to a given fault case.

4.1. Fault Location: Simulated Cases

The quantities expressed in Table 8 refer to the number of simulated cases for validation of ANN and the application of NLO techniques and the AM method, which depends on the type of fault and the methods implemented. A statistical study is needed to compare the methods used in this study. Owing to the stochastic characteristics of the location by the ANN and by the PRGA, the means and variances of the location values of these methods were obtained from 40 executions.
Following the implementation and execution of the location step, the balancing of the number of samples per class was applied for the adequacy of the data used. Accordingly, 24 samples were randomly selected from each of the possible combinations of the studied factors, a number corresponding to the smaller quantity of samples of the relationship between the fault event and the implemented methods, as listed in Table 8. Associated with the Anova2 test, the results achieved in the balanced experiment allow the use of the differences between the averages of the levels for the estimates of the main effects of the factors and the interaction between them. When the experiment is unbalanced, differences in factor level means may be associated with unbalanced observations rather than changes in factor levels.
From Equation (5), the percentage errors of the results indicated for each fault location were calculated, provided by the difference between the estimated value by each method and the simulated value in relation to the length of the line.
e % = e s t i m a t e d   l o c a t i o n s i m u l a t e d   l o c a t i o n l i n e   l e n g t h x 100
In simulated scenarios, regardless of the type of fault that occurred, the analysis of the Box-plots of the percentage errors, represented in Figure 15, revealed the better performance of the ANN.
The verification of the existence of differences between the types of methods was based on the Anova2 statistical test. Particularly, considering the results provided by the AM method and the NLO methods, when conducting the analysis of variance for the comparison of samples, it was assumed that the observations were independent and normally distributed, maintaining the variance constant in each treatment. According to [52], the graphic results expressed in Figure 16 made it possible to verify that these assumptions were satisfied.
In the quantile-quantile plot of Normal shown in Figure 16a, the configuration of points approached a straight line, indicating that the residuals are normally distributed. The residuals shown in Figure 16b were located approximately around a horizontal band, signaling the validity of the independence assumption. Regarding the homoscedasticity of the levels of the type of fault and type of method factors, Figure 16c,d present the residual plots from randomly arranged points without a set pattern for level change, having adjusted values both for the type of fault factor and type of method factor. In addition to graphical analysis, the verification of the premises for the application of Anova2 was based on statistical tests, as indicated in [53]. In this sense, the randomness of the residues was verified by applying the Durbin-Watson test (DW), with the result DW = 1.805907. The DW test statistic is in the range between 0 and 4; values close to 2 indicate no correlation between the residues, ensuring the randomness of the sample. The normality of the residues, with 5% significance, was verified by the Shapiro-Wilk hypothesis test (SW), with a p-value result (SW) = 6.023376 × 10−2. To test the normal distribution, it is assumed that H0 is symmetric, whereas H1 assumes that the variable distribution is asymmetrical. In the Shapiro-Wilk test, the test’s significance was greater than 0.050, p-value > 0.050, indicating that the residues follow a normal distribution. Regarding the homoscedasticity of the residues, Bartlett’s test (B) was used, which is not affected by the sample size and should be used when the residues present normal distribution, and also the Levene test (L), which, although limited to balanced samples, does not require the assumption of normality. In these tests, H0 indicates that the residue variances are homoscedastic, and H1 assumes that at least one variance differs from the others. For the type of fault, the tests provided as a result p-value (B) = 6.071991 × 10−1, p-value (L) = 3.748992 × 10−1 and for type of method p-value (B) = 2.511740 × 10−1, p-value (L) = 5.339021 × 10−2.
As presented in Table 9, the application of the Anova2 test leads to the rejection of the null hypothesis for both factors analyzed. Specifically, about the location methods evaluated, the indication is that at a 5% significance level, the effects of the type of method affect the percentage error for the fault location, with at least one of the location methods that differ from the others. The result of the p-value (PR) obtained by the test statistic was close to 5.816931 × 10−12, a value considerably lower than 5%. With the data from Table 9, it is also possible to make inferences regarding the types of faults. However, this analysis was not addressed in the proposal of this study.
As a complement to the Anova2 study in simulated cases, the Tukey test was applied at a significance level of 5%, and the results expressed in Table 10 were obtained.
The Tukey test is interpreted based on the value of the minimum significant difference (MSD), obtained from the distribution of the studentized amplitude, the mean square of the residues of Anova2, and the sample size of the groups, in the confidence interval determined. For the test, represented in Table 9 and Figure 17, the modulus of the mean difference between the pairs of methods used In the fault’s location was greater than the MSD value obtained for the pairs AM—QN, EL—AM, and PRGA—AM. The value 0 (zero) is not contained in the confidence intervals of these pairs, indicating that the average performance is significantly different between them. Furthermore, the percentage errors provided by the AM method were greater than those provided by the three NLO methods.

4.2. Fault Location: Real Cases

The procedures used for the pilot study were applied in other TL. In addition to the 74.4 km transmission line used initially, the techniques for locating faults in three other lines were implemented. Accordingly, the implemented methods were tested in seven real cases of short circuits caused by atmospheric discharge (AD) and fire in the modeled lines, as indicated in Table 11.
The results achieved by the methods for the real cases are recorded in Table 12.
For the statistical analysis of the results, the location errors were calculated from the difference between the location estimated by the method and the actual distance from the defect (provided by the inspection team) in relation to the length of the line where the bending occurred. For the ANN and the PRGA, the errors were indicated from the average error of the 40 executions.
To identify the existence of differences in the location results that occurred depending on the type of method used, the verification of the randomness, normality, and homoscedasticity supported the assumption that the observations are independent, normally distributed, and maintain the same variance in each treatment, allowing the application of the Anova2 model.
Using the results of the test application presented in Table 13, the null hypothesis of equality between the methods was discarded, and the result of the p-value obtained by the test statistic was close to 9.92 × 10−4, which is considerably less than 5%.
Tukey test results presented in Table 14 made it possible to identify a difference in the location of faults performed by the ANN, which provided greater location errors than the other methods.
In real cases, using the Tukey test, as shown in Table 14 and Figure 18, there was no statistical evidence identified for the difference between the location of faults carried out between the different optimization methods implemented, nor between these methods and the classic analytical. The module of the mean difference between the pairs of methods used was greater than the value of the MSD obtained for the pairs PRGA—ANN, EL—ANN, AM—ANN, and QN—ANN. The value 0 (zero) is not contained in the interval’s confidence of these peers, indicating that average performance is significantly different between them. Furthermore, the worst performance was presented by the ANN.
An interesting point to be highlighted for the real cases is that, although the objective functions for the optimization algorithm were developed for the short transmission line model (with only the parameters of resistance and series inductance concentrated), there was no significant difference in the errors obtained in relation to the analytical method, which considers the long-line model. It is noteworthy that only seven cases were investigated (for lines of 74.4 km, 105.58 km, 248.44 km, and 342.71 km), requiring further studies from a larger database. For the cases presented, the short line model used for the transmission line was not preponderant for the results, probably due to other errors such as those caused by current and potential transformer, line parameters, mutual inductance, transposition, and phasor estimation method.
The neural networks, as they have been conventionally used, with voltage and current phasors as input structures, provide reduced errors in simulated cases, but they own the difficulty generalizing for the real fault location function. It should be mentioned that the electrical system model in simulated cases is the same for training and validation, which leads the network to obtain accurate results. The largest errors of the ANN in the real fault location process are related to the values of the Thevenin equivalents (sources and impedance) at the local and remote ends of the transmission line. Electrical systems are dynamic, with load and generation varying 24 h a day. In order to generate the fault files used in the training of location ANN on real cases, the Thevenin voltages and impedances provided by the electric utility were inserted into the ATP, calculated from a short circuit program that considers the static system with light or heavy load. However, the exact values of these equivalents at the time of the fault cannot be determined, which compromises the results of neural networks applied to real cases, contributing to unsatisfactory results. To mitigate the problem, a neural network fault location method, in which the Thevenin equivalents of the line terminals are not needed, is already being developed by the authors. The proposal will enable the practical application of these structures to solve this type of problem in conventional lines and on lines with series compensation, in which fault location becomes more complex due to the non-linearity of the Metal Oxide Varistor, the capacitor’s protection element.

5. Conclusions

This study aimed to present and statistically compare applications of artificial neural networks, nonlinear optimization techniques, and a classical analytical method for the problem of fault location in transmission lines in simulated and real scenarios. In the simulated scenarios, the analysis of the percentage errors of the location revealed, with 5% of significance, that the type of method used affects the indication of the location of the short circuit. Smaller errors were observed with neural networks and larger ones with the analytical method, leaving the nonlinear optimization methods with intermediate errors. In real scenarios, unlike the situation verified in simulated environments, there was statistical evidence of higher percentages of errors with neural networks, not rejecting the equality between the location indicated by the analytical method and the nonlinear optimization methods. The neural networks conventionally used for the location function proved incompatible for application in companies that operate electrical systems. No differences were identified among the nonlinear optimization methods. The three methods implemented and applied to the proposed function converged to the expected local minimum. However, it should be noted that the low computational cost and fast response associated with the Quasi-Newton method make it a prominent method. Moreover, the nonlinear optimization methods used for the location function proved to be promising for application in companies that operate electrical systems, providing location errors similar to those presented by the classical analytical method. In this sense, and based on the proposal of this study, it is suggested to consider the application of the methods Quasi-Newton and analytic to improve the answer to the problem of locating faults in transmission lines.

Author Contributions

Conceptualization, S.A.R. and E.G.S.; methodology, S.A.R. and E.G.S.; software, R.T.N.C. and S.A.R.; validation, S.A.R. and E.G.S.; supervision, R.T.N.C. and T.G.M.; writing, S.A.R., E.G.S. and T.G.M. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Federal Center for Technological Education of Minas Gerais.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Z.; Maun, J.-C. Artificial neural network approach to single-ended fault locator for Transmission Lines. IEEE Trans. Power Syst. 2000, 15, 370–375. [Google Scholar] [CrossRef]
  2. Joorabian, M. Artificial neural network based fault locator for EHV transmission system. In Proceedings of the 10th Mediterranean Electrotechnical Conference. Information Technology and Electrotechnology for the Mediterranean Countries, Lemesos, Cyprus, 29–31 May 2000; Volume 3, pp. 1003–1006. [Google Scholar]
  3. Mazon, A.J.; Zamora, I.; Gracia, J.; Sagastabeutia, K.J.; Saenz, J.R. Selecting ANN structures to find transmission faults. IEEE Comput. Appl. Power 2001, 14, 44–48. [Google Scholar] [CrossRef]
  4. Martin, F.; Aguado, J.A. Wavelet-based ANN approach for transmission line protection. IEEE Trans. Power Deliv. 2003, 18, 1572–1574. [Google Scholar] [CrossRef]
  5. Amorim, H.P.; Huais, L. Faults location in transmission lines through neural networks. In Proceedings of the IEEE/PES Transmision and Distribution Conference and Exposition: Latin America, Sao Paulo, Brazil, 8–11 November 2004; pp. 691–695. [Google Scholar]
  6. Pereira, C.E.M.; Zanetta, L.C. Fault location in transmission lines using one-terminal postfault voltage data. IEEE Trans. Power Deliv. 2004, 19, 570–575. [Google Scholar] [CrossRef]
  7. Gracia, J.; Mazon, A.J.; Zamora, I. Best ANN structures for fault location in single-and double-circuit transmission lines. IEEE Trans. Power Deliv. 2005, 20, 2389–2395. [Google Scholar] [CrossRef]
  8. Lin, X.; Mao, P.; Weng, H.; Wang, B.; Bo, Z.Q.; Klimek, A. Study on Fault Location for High Voltage Overhead Transmision Lines Based on Neural Network System. In Proceedings of the International Conference on Intelligent Systems Applications to Power Systems, Kaohsiung, Taiwan, 5–8 November 2007; pp. 1–5. [Google Scholar]
  9. Hagh, M.T.; Razi, K.; Taghizadeh, H. Fault classification and location of power transmission lines using artificial neural network. In Proceedings of the International Power Engineering Conference, Singapore, 3–6 December 2007; pp. 1109–1114. [Google Scholar]
  10. Pereira, C.E.M.; Zanetta, L.C., Jr. An optimization approach for fault location in transmission lines using one terminal data. Int. J. Electr. Power Energy Syst. 2007, 29, 290–296. [Google Scholar] [CrossRef]
  11. Souza, S.M.; Silva, A.P.A.; Lima, C.S. Voltage and Current Patterns for Fault Location in Transmission Lines. In Proceedings of the International Joint Conference on Neural Networks, Orlando, FL, USA, 12–17 August 2007; pp. 737–742. [Google Scholar]
  12. Ekici, S.; Yildirim, S.; Poyraz, M. A transmission line fault locator based on Elman recurrent networks. Sci. Direct Appl. Soft Comput. 2009, 9, 341–347. [Google Scholar] [CrossRef]
  13. Ngaopitakkul, A.; Pothisarn, C. Discrete wavelet transform and back-propagation neural networks algorithm for fault location on single-circuit transmission line. In Proceedings of the IEEE International Conference on Robotics and Biomimetics, Guilin China, 19–23 December 2009; pp. 1613–1618. [Google Scholar]
  14. Sahoo, S.; Ray, P.; Panigrahi, B.K.; Senroy, N. A computational intelligence approach for fault location in transmission lines. In Proceedings of the Joint International Conference on Power Electronics, Drives and Energy Systems & 2010 Power India, New Delhi, India, 20–23 December 2010; pp. 1–6. [Google Scholar]
  15. Wei, L.; Guo, W.; Wen, F.; Ledwich, G.; Liao, Z.; Xin, J. Waveform Matching Approach for Fault Diagnosis of a High-Voltage Transmission Line Employing Harmony Search Algorithm. IET Gener. Transm. Distrib. 2010, 4, 801–809. [Google Scholar] [CrossRef]
  16. Jiang, J.A.; Chuang, C.L.; Wang, Y.C.; Hung, C.H.; Wang, J.Y.; Lee, C.H.; Hsiao, Y.T. A Hybrid Framework for Fault Detection, Classification, and Location—Part I: Concept, Structure, and Methodology. IEEE Trans. Power Deliv. 2011, 26, 1988–1998. [Google Scholar] [CrossRef]
  17. Lout, K.; Aggarwal, R.K. A feedforward Artificial Neural Network approach to fault classification and location on a 132kV transmission line using current signals only. In Proceedings of the 47th International Universities Power Engineering Conference (UPEC), London, UK, 4–7 September 2012; pp. 1–6. [Google Scholar]
  18. Teklić, L.; Filipović-Grčić, B.; Pavičić, I. Artificial neural network approach for locating faults in power transmission system. Eurocon 2013, 1425–1430. [Google Scholar]
  19. Raoofat, M.; Mahmoodian, A.; Abunasri, A. Fault location in transmission lines using neural network and wavelet transform. In Proceedings of the International Congress on Electric Industry Automation, Shiraz, Iran, 24–25 February 2015; pp. 1–6. [Google Scholar]
  20. Chaitanya, B.K.; Yadav, A.; Andanapalli, K.; Varma, B.R.K. A comparative study of different signal processing techniques for fault location on transmission lines using hybrid Generalized Regression Neural Network. In Proceedings of the 2016 International Conference on Signal Processing, Communication, Power and Embedded System (SCOPES), Odisha, India, 3–5 October 2016; pp. 1246–1250. [Google Scholar]
  21. Ahmed, A.S.; Attia, M.A.; Hamed, N.M.; Abdelaziz, A.Y. Modern optimization algorithms for fault location estimation in power systems. Eng. Sci. Technol. Int. J. 2017, 20, 1475–1485. [Google Scholar]
  22. Carrión, D.; González, J.W.; IssacI, A.; López, G.J. Optimal fault location in transmission lines using hybrid method. In Proceedings of the IEEE PES Innovative Smart Grid Technologies Conference—Latin America (ISGT Latin America), Quito, Ecuador, 20–22 September 2017; pp. 1–6. [Google Scholar]
  23. Fahim, S.R.; Sarker, Y.; Islam, O.K.; Sarker, S.K.; Ishraque, M.F.; Das, S.K. An intelligent approach of fault classification and localization of a power transmission line. In Proceedings of the IEEE International Conference on Power, Electrical, and Electronics and Industrial Applications, Dhaka, Bangladesh, 29 November–1 December 2019; pp. 53–56. [Google Scholar]
  24. Wang, C.; Yun, Z. Parameter-free fault location algorithm for distribution network T-Type transmission lines. Energies 2019, 12, 1534. [Google Scholar] [CrossRef] [Green Version]
  25. Elnozahy, A.; Sayed, K.; Bahyeldin, M. Artificial neural network based fault classification and location for transmission lines. In Proceedings of the IEEE Conference on Power Electronics and Renewable Energy, Aswan, Egypt, 23–25 October 2019; pp. 140–144. [Google Scholar]
  26. Onaolapo, A.K.; Pillay Carpanen, R.; Dorrell, D.G.; Ojo, E.E. Transmission line fault classification and location using multi-layer perceptron artificial neural network. In Proceedings of the IECON 2020 The 46th Annual Conference of the IEEE Industrial Electronics Society, Singapore, 18–21 October 2020; pp. 5182–5187. [Google Scholar]
  27. Coban, M.; Tezcan, S.S. Artificial neural network based fault location on 230 kv transmission line using voltage and current signals. In Proceedings of the 4th International Symposium on Multidisciplinary Studies and Innovative Technologies (ISMSIT), Istanbul, Turkey, 22–24 October 2020; pp. 1–4. [Google Scholar]
  28. Kezunovic, M.; Lou, S.; Galijasevic, Z.; Ristanovic, D. Accurate fault location in transmission networks using modeling, simulation and limited field recorded data. Power Syst. Eng. Res. Cent. 2002, 2–44. [Google Scholar]
  29. Ezquerra, J.; Valverde, V.; Mazón, A.J.; Zamora, I.; Zamora, J.J. Field programmable gate array implementation of a fault location system in transmission lines based on artificial neural networks. IET Gener. Transm. Distrib. 2011, 5, 191–198. [Google Scholar] [CrossRef]
  30. Joorabian, M.; Taleghani, S.M.A.; Aggarwal, R.K. Accurate fault locator for EHV transmission lines based on radial basis function neural networks. Electr. Power Syst. Res. 2015, 71, 195–202. [Google Scholar] [CrossRef]
  31. Silveira, E.G.; Paula, H.R.; Rocha, S.A.; Pereira, C.S. Hybrid fault diagnosis algorithms for transmission lines. Electr. Eng. 2017, 100, 1689–1699. [Google Scholar] [CrossRef]
  32. Johns, A.T.; Jamali, S. Accurate fault location technique for power transmission line. IEEE Proc. C Gener. 1990, 137, 395–402. [Google Scholar] [CrossRef]
  33. Zhang, X.; Pang, B.; Liu, Y.; Liu, S.; Xu, P.; Li, Y.; Xie, Q. Review on detection and analysis of partial discharge along power cables. Energies 2021, 14, 7692. [Google Scholar] [CrossRef]
  34. Yaacob, M.M.; Alsaedi, M.A.; Rashed, J.R.; Dakhil, A.M.; Atyah, S.F. Review on partial discharge detection techniques related to high voltage power equipment using different sensors. Photonic Sens. 2014, 4, 325–337. [Google Scholar] [CrossRef] [Green Version]
  35. Birlasekaran, S.; Li, H.J. Detection of faulty insulators on power transmission line. In Proceedings of the 2000 IEEE Power Engineering Society Winter Meeting, Singapore, 23–27 January 2000; Volume 4, pp. 2817–2821. [Google Scholar]
  36. Han, J.; Yang, Z.; Zhang, Q.; Chen, C.; Li, H.; Lai, S.; Hu, G.; Xu, C.; Xu, H.; Wag, D.; et al. A method of insulator faults detection in aerial images for high-voltage transmission lines inspection. Appl. Sci. 2019, 9, 2009. [Google Scholar] [CrossRef] [Green Version]
  37. Kuwalek, P. AM modulation signal estimation allowing further research on sources of voltage fluctuations. IEEE Trans. Ind. Electron. 2019, 67, 6937–6945. [Google Scholar] [CrossRef]
  38. Kuwałek, P. Selective identification and localization of voltage fluctuation sources in power grids. Energies 2021, 14, 6585. [Google Scholar] [CrossRef]
  39. Taqi, A.; Beryozkina, S. Overhead transmission line thermographic inspection using a drone. In Proceedings of the IEEE 10th GCC Conference & Exhibition (GCC), Kuwait, Kuwait, 19–23 April 2019; pp. 1–6. [Google Scholar]
  40. EMTP Center. Alternative Transients Program Rule Book; EMTP Center: Leuven, Belgium, 1990. [Google Scholar]
  41. Yu, T.C.; Martí, J.R. A robust phase-coordinates frequency-dependent underground cable model (zCable) for the EMTP. IEEE Trans. Power Deliv. 2003, 18, 189–194. [Google Scholar] [CrossRef]
  42. Marti, L. Simulation of transients in underground cables with frequency-dependent modal transformation matrices. IEEE Trans. Power Deliv. 1988, 3, 1099–1110. [Google Scholar] [CrossRef]
  43. Coury, D.V.; Oleskovicz, M.; Aggarwal, R.K. An ANN routine for fault detection, classification and location in transmission lines. Electr. Power Compon. Syst. 2002, 30, 1137–1149. [Google Scholar] [CrossRef]
  44. Santos, A.; Gaspar, C.; Barros, M.T.C.; Duarte, P. Transmission line fault resistance values based on field data. IEEE Trans. Power Deliv. 2020, 35, 1321–1329. [Google Scholar] [CrossRef]
  45. Andrade, V.D.; Sorrentino, E. Typical expected values of the fault resistance in power systems. In Proceedings of the IEEE/PES Transmission and Distribution Conference and Exposition: Latin America, Sao Paulo, Brazil, 8–10 November 2010; pp. 602–609. [Google Scholar]
  46. Phadke, A.G.; Thorp, J.S. Computer Relaying for Power Systems; John Wiley & Sons Inc.: Hoboken, NJ, USA, 1998; p. 289. [Google Scholar]
  47. Demuth, H.; Beale, M. Neural Network Toolbox User’S Guide; The MathWork: Natick, MA, USA, 2015. [Google Scholar]
  48. Luenberger, D.G.; Ye, Y. Linear and Nonlinear Programming; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  49. Bland, R.G.; Goldfarb, D.; Todd, M.J. Feature Article—The ellipsoid method: A survey. Inst. Oper. Res. Manag. Sci. 1981, 29, 1039–1091. [Google Scholar] [CrossRef] [Green Version]
  50. Takahashi, R.H.C.; Vasconcelos, J.A.; Ramirez, J.A.; Krahenbuhl, L. A multiobjective methodology for evaluating genetic operators. IEEE Trans. Magn. 2003, 39, 1321–1324. [Google Scholar] [CrossRef]
  51. Python Software Foudation. Python Language Site: Documentation. 2022. Available online: https://www.python.org/psf/ (accessed on 7 January 2022).
  52. Montgomery, D.C.; Runger, G.C. Applied Statistics and Probability for Engineers; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  53. Montgomery, D.C.; Peck, E.A.; Vining, G.G. Introduction to Linear Regression Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2021. [Google Scholar]
Figure 1. A transmission line with a phase A to ground fault.
Figure 1. A transmission line with a phase A to ground fault.
Energies 15 04095 g001
Figure 2. ATP model for short circuit simulation.
Figure 2. ATP model for short circuit simulation.
Energies 15 04095 g002
Figure 3. Types of faults: (a) SPG (AG, BG, CG); (b) DP (AB, AC, BC); (c) TP (ABC); (d) DPG (ABG, ACG, BCG).
Figure 3. Types of faults: (a) SPG (AG, BG, CG); (b) DP (AB, AC, BC); (c) TP (ABC); (d) DPG (ABG, ACG, BCG).
Energies 15 04095 g003
Figure 4. Fault location process.
Figure 4. Fault location process.
Energies 15 04095 g004
Figure 5. Voltage and current input waves—AG fault.
Figure 5. Voltage and current input waves—AG fault.
Energies 15 04095 g005
Figure 6. Filtered voltage and current waves—AG fault.
Figure 6. Filtered voltage and current waves—AG fault.
Energies 15 04095 g006
Figure 7. The amplitude of voltage and current phasors—AG fault.
Figure 7. The amplitude of voltage and current phasors—AG fault.
Energies 15 04095 g007
Figure 8. Detail of selection of phasors from the missing phase A current signal, after low-pass filtering.
Figure 8. Detail of selection of phasors from the missing phase A current signal, after low-pass filtering.
Energies 15 04095 g008
Figure 9. The Modular structure of fault location ANN.
Figure 9. The Modular structure of fault location ANN.
Energies 15 04095 g009
Figure 10. Objective function: fault AG (dF = 45 km, RF = 20 Ω).
Figure 10. Objective function: fault AG (dF = 45 km, RF = 20 Ω).
Energies 15 04095 g010
Figure 11. Fault location and resistance: Quasi-Newton.
Figure 11. Fault location and resistance: Quasi-Newton.
Energies 15 04095 g011
Figure 12. Fault location and resistance: Ellipsoidal.
Figure 12. Fault location and resistance: Ellipsoidal.
Energies 15 04095 g012
Figure 13. Fault location and resistance: PRGA.
Figure 13. Fault location and resistance: PRGA.
Energies 15 04095 g013
Figure 14. Factors and levels of the location problem.
Figure 14. Factors and levels of the location problem.
Energies 15 04095 g014
Figure 15. Location errors: simulated cases.
Figure 15. Location errors: simulated cases.
Energies 15 04095 g015
Figure 16. Assumptions for the application of the Anova2 model. (a) Data normality; (b) Data randomness; (c) Homoscedasticity—Type of fault; (d) Homoscedasticity—Type of method.
Figure 16. Assumptions for the application of the Anova2 model. (a) Data normality; (b) Data randomness; (c) Homoscedasticity—Type of fault; (d) Homoscedasticity—Type of method.
Energies 15 04095 g016
Figure 17. Tukey test: simulated cases.
Figure 17. Tukey test: simulated cases.
Energies 15 04095 g017
Figure 18. Tukey’s test: real cases.
Figure 18. Tukey’s test: real cases.
Energies 15 04095 g018
Table 1. Signal characterization and preprocessing.
Table 1. Signal characterization and preprocessing.
SignalsReferences
ScenarioSimulated[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27].
Simulated and Real[28,29,30,31].
OriginOne Terminal[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30].
Two Terminals/PMUs[3,7,13,14,18,19,20,21,31].
SignCurrent[14,17,23].
Voltage[19].
Current and Voltage[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16].
Pre-processingPhasors/Discrete Fourier Transform/Least Squares[1,3,5,6,7,8,9,10,11,18,22,24,26,28,30,31].
Frequency Spectrum/Discrete Fourier Transform/Least Squares[2,17,27,30].
Energy/detail coefficient of the Discrete Wavelet Transform[4,12,13,14,15,16,19].
Time signals[15,21,23,25,29].
Table 2. Characterization of ANN in fault location studies.
Table 2. Characterization of ANN in fault location studies.
ANNReferences
ModelMulti Layer Perceptron[1,2,3,4,5,6,7,9,11,12,13,14,18,19,23,25,26,27,29,30,31].
Others[3,12,16,20].
FeedingFeedforward[1,2,4,7,8,9,13,14,17,18,25,26,27,28,29,30,31].
LearningBackpropagation[1,2,7,9,11,13,14,18,23,25,26,27,29,31].
TrainingLevenberg-Marquardt[1,2,7,9,13,14,18,19,25,26,27,29,31].
Others[16,23].
Activation FunctionSigmoid[9,19,27,29].
Hyperbolic Tangent[1,13,14,31].
Logarithmic[19,25,27].
Gaussian[30].
Table 3. Line parameters and equivalent fonts.
Table 3. Line parameters and equivalent fonts.
ElementPositive SequenceZero Sequence
r 1   ( Ω ) x 1   ( Ω ) r 0   ( Ω ) x 0   ( Ω )
Line2.6927.9726.94106.58
Local Source4.00134.1104.06933.603
Remote Source6.33253.8452.73139.363
Table 4. Composition of fault scenarios: Training of ANN.
Table 4. Composition of fault scenarios: Training of ANN.
Training DataTotal
dF (km)Every 5% of the line19 locations
R F (Ω)Cases SPG e DPG: 0-6-12-18-24-30-36-428 values
Cases DP e TP: 0-1-2-3-4-56 values
Number of scenariosSPG e DPG: 2∙3∙8∙19 = 9121368 scenarios
DP: 3∙6∙19 = 342
TP: 1∙6∙19 = 114
Table 5. Composition of fault scenarios: validation of ANN, application in NLO techniques and in the analytical method.
Table 5. Composition of fault scenarios: validation of ANN, application in NLO techniques and in the analytical method.
Validation and Application DataTotal
dF (km)4-11-17-26-41-53-59-688 locations
R F (Ω)Cases SPG e DPG: 10-21-26-31-345 values
Cases DP e TP: 2.5-3.5-4.53 values
Number of scenariosSPG e DPG: 2∙3∙8∙5 = 240336 scenarios
DP: 3∙3∙8 = 72
TP: 1∙3∙8 = 24
Table 6. Objective functions.
Table 6. Objective functions.
Type of Fault F d F , R F
AG, BG, CG V S d F Z 1 I S + Z 0 Z 1 Z 1 I S 0 3 R F I S 1 + I R 1 2
AB, BC, CA V S 1 V S 2 d F Z 1 I S 1 + I S 2 2 R F I S 1 + I R 1 2
ABG, BCG, CAG V S 1   V S 0 d F Z 1 I S 1 + d F Z 0 I S 0 + 3 R F I S 0 + I R 0 2
ABC V S 1 d F Z 1 I S 1 R F I S 1 + I R 1 2  
Table 7. The number of iterations and execution time per method.
Table 7. The number of iterations and execution time per method.
MethodANNPRGAQNELAM
Number of iterations-42569-
Time (s)579160.70.80.6
Table 8. Validation cases by fault type and method.
Table 8. Validation cases by fault type and method.
Type of MethodANNPRGAQNELAM
Type of FaultDPG48004800120120120
SPG48004800120120120
DP28802880727272
TP960960242424
Table 9. Anova2 testes: simulated cases.
Table 9. Anova2 testes: simulated cases.
sum_sqDfFPR(>F)
C(Type_Fault)0.2914823.017.3840731.362907 × 10−10
C(Type_Method)0.3326683.019.8404635.816931 × 10−12
Residual2.107073377.0NaNNaN
Table 10. Multiple comparisons of means: Tukey HSD, FWER = 0.05—simulated cases.
Table 10. Multiple comparisons of means: Tukey HSD, FWER = 0.05—simulated cases.
Group1Group2MeandiffLowerUpperReject
PRGAEL−0.0162−0.04570.0134False
PRGAAM0.05220.02260.0818True
PRGAQN−0.0228−0.05240.0068False
ELAM0.06840.03880.0980True
ELQN−0.0067−0.03620.0229False
AMQN−0.0750−0.1046−0.0454True
Table 11. Characterization: real fault scenarios.
Table 11. Characterization: real fault scenarios.
LineVoltage (kV)Length (km)TypeLocation (km)Cause
134574.40AG60.0AD
BT54.0AD
2500105.58AG30.0Fire
3500342.71CT317.0AD
AG76.0Fire
CT55.0Fire
4500248.44ACT91.0AD
Table 12. Location results: real faults.
Table 12. Location results: real faults.
LineInspection Results
(km)
Fault Location (km)
ANNQNELPRGAAM
16055.665.465.465.464.4
5433.554.454.454.454.3
23024.023.623.623.629.7
3317306.4311.0311.0311.0324.0
76103.078.478.478.482.7
5560.555.755.755.756.6
491124.091.591.591.594.1
Average Error (%)9.22.42.42.41.7
Median Error (%)6.81.21.21.21.4
Table 13. Anova2 test: real cases.
Table 13. Anova2 test: real cases.
sum_sqdfFPR(>F)
C(Type_Fault)266.1309276.01.2396800.321484
C(Type_Method)944.3113114.06.5981270.000992
Residual858.70848824.0NaNNaN
Table 14. Multiple comparisons of means: Tukey HSD, FWER = 0.05—real cases.
Table 14. Multiple comparisons of means: Tukey HSD, FWER = 0.05—real cases.
Group1Group2MeandiffLowerUpperReject
PRGAEL−1.0258−10.51938.4678False
PRGAAM−1.6714−11.1657.8221False
PRGAQN−1.0259−10.51958.4676False
PRGAANN11.98572.492221.4793True
ELAM−0.6457−10.13928.8479False
ELQN−0.0001−9.49379.4934False
ELANN13.01153.517922.505True
AMQN0.6455−8.84810.1391False
AMANN13.65714.163623.1507True
QNANN13.01163.518122.5052True
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rocha, S.A.; Mattos, T.G.; Cardoso, R.T.N.; Silveira, E.G. Applying Artificial Neural Networks and Nonlinear Optimization Techniques to Fault Location in Transmission Lines—Statistical Analysis. Energies 2022, 15, 4095. https://doi.org/10.3390/en15114095

AMA Style

Rocha SA, Mattos TG, Cardoso RTN, Silveira EG. Applying Artificial Neural Networks and Nonlinear Optimization Techniques to Fault Location in Transmission Lines—Statistical Analysis. Energies. 2022; 15(11):4095. https://doi.org/10.3390/en15114095

Chicago/Turabian Style

Rocha, Simone A., Thiago G. Mattos, Rodrigo T. N. Cardoso, and Eduardo G. Silveira. 2022. "Applying Artificial Neural Networks and Nonlinear Optimization Techniques to Fault Location in Transmission Lines—Statistical Analysis" Energies 15, no. 11: 4095. https://doi.org/10.3390/en15114095

APA Style

Rocha, S. A., Mattos, T. G., Cardoso, R. T. N., & Silveira, E. G. (2022). Applying Artificial Neural Networks and Nonlinear Optimization Techniques to Fault Location in Transmission Lines—Statistical Analysis. Energies, 15(11), 4095. https://doi.org/10.3390/en15114095

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