Prediction of Metallic Conductor Voltage Owing to Electromagnetic Coupling Via a Hybrid ANFIS and Backtracking Search Algorithm

: The electromagnetic interference (EMI) generated by high voltage power systems can cause a serious problem for nearby electrically conductive structures, such as railroads, communication lines, or pipelines, that would place a system’s integrity and the operational safety of the structure at high level of risk. According to the IEEE standard-80, by implementing a well-designed mitigation system, the induced voltage on neighboring electrically conductive structure can reach a harmless level. The mitigation system can enhance the overall integrity of pipelines and provide higher operation safety for personal during working on the exposed parts of metallic pipelines or conductive appurtenances. An accurate prediction about the level of induced voltage is absolutely necessary to design a suitable mitigation system for metallic pipelines. Thus, in this work a hybrid prediction methodology composed of an adaptive neuro-fuzzy inference system (ANFIS) and a backtracking search algorithm (BSA) is developed to accurately predict the electromagnetic inference’s e ﬀ ects on metallic pipelines with shared right-of-way (RoW) and high voltage overhead lines (OHLs). Through the combination of BSA as a robust and e ﬃ cient optimization algorithm in the learning process of an ANFIS approach, a hybrid data mining algorithm has been developed to predict the induced voltage on mitigated and unmitigated pipelines more accurately and reliably. The simulation results are validated by data sets observed from the Current Distribution, Electromagnetic Interference, Grounding and Soil Structure Analysis (CDEGS) software. From the simulation results it was conﬁrmed that the proposed hybrid method is e ﬀ ective in accurately predicting the induced voltage on pipelines with changing system parameters. Furthermore, to evaluate the precision and applicability of the developed approach in this paper, its estimates are compared with the results obtained from an artiﬁcial neural network (ANN), a support vector regression (SVR) and an ANFIS optimized by other well-known optimization algorithms. The obtained results indicate higher accuracy of the developed hybrid method over other artiﬁcial intelligence based approaches.


Introduction
Due to the quick expansion of the economy in different countries, the demand for primary energy resources, raw materials and electrical energy is serially growing.Thus, to effectively supply such continuously increasing demands, it is essential to expand the existing high voltage overhead lines (OHLs), water, gas, and oil supply pipelines or construct new ones.To reduce both construction costs and environmental damage, a set of government regulations has been generated that limit the access to new distribution and transmission corridors [1][2][3].Therefore, high voltage OHLs are constructed with the shared transmission and distribution corridors for water, gas, and oil pipelines.When pipelines are located in shared right-of-way (RoW) with high voltage OHLs, the neighboring pipelines would suffer from very high induced currents and voltages, due to the electromagnetic interference (EMI) effects generated by high voltage OHLs [4].
The EMI can be transferred from a high voltage AC power system to nearby metallic constructions without any electrical connection.The AC interference is created in the neighboring metallic water, gas, and oil pipelines due to the electromagnetic fields produced by high voltage OHLs.Thus, in the most conditions above or underground, metallic pipelines are vulnerable to be effected by induced high AC voltages and currents [5].In the most severe conditions the metallic pipelines are located in close adjacency to high voltage OHLs without any electrical connection to a mitigation system for reduction in the level of induced voltage.These situations can be more severe during an OHL fault while the level of induced voltage on unmitigated metallic pipelines can extend to thousands of volts [6].
According to the circuit configurations, capacitive, inductive, and conductive couplings are considered as three different categories of AC interference coupling mechanisms [7].
Conductive coupling is generated through the ground voltage rise caused by discharging a very high amount of current in to the ground at power system structures, especially at the grounding systems of high voltage OHLs, high voltage substations, and power plants.The conductive coupling is one of the main concerns during power system faults, particularly in the places where the metallic pipelines are located in close proximity to high voltage OHLs [8].
Generally, inductive coupling is generated via the magnetic fields.The OHLs carrying a high level of current by their conductors generate strong magnetic fields in the neighboring area.The generated magnetic field around the high voltage OHLs induces voltage in the surrounding metallic structures coupled by this magnetic field.This coupling is influenced by different factors, such as the level of the OHLs' current, the parallelism length, and the distance between OHLs and metallic structures [9].
The capacitive coupling is generated by the potential difference between two systems.The voltage difference between high voltage OHLs and any nearby conductive structure, such as a metallic pipeline, forms an electric field between the two systems.
To avoid the electrochemical corrosion and minimize the densities of the induced current, an appropriate protection is essential for the underground metallic pipelines that are in permanent touch with the electrolyte solution of the ground [10].Uncontrolled corrosion of the buried metallic supply pipelines can cause gas/oil spills with severe economic and ecological implications.
The induced voltage on the metallic pipeline would put the safety of persons who will touch the exposed metallic parts of gas/oil pipeline at risk, due to a potentially powerful electric shock; and could also damage the cathodic protection (CP) of the pipelines that is used for protecting the pipeline from corrosion, and place the integrity of the pipeline at serious risk [11].
Additionally, extreme coating stress potentials (the potential difference between the metallic pipeline and local ground) would damage the coating, because of sped-up corrosion.If the level of voltage is sufficiently strong it can harm the gas/oil pipeline walls.This voltage can also result in damage to both insulating flanges and CP equipment [12].
The level of induced voltage on the pipelines can be reduce to a safe range, according to the IEEE standard-80, by implementing an appropriate mitigation system [10].The gradient control wire is considered one of the greatest mitigation systems of those that are widely applied [13,14].The AC interference has to be well studied to design a suitable protection for metallic pipelines [15].Especially for designing an optimized mitigation system, a proper methodology is required to accurately predict the induced voltage on metallic pipelines.

Literate Review
Multitudinous interrelations exist in the calculation of AC induced currents and voltages on metallic gas, oil, or water supply pipelines placed in the close proximity of high voltage OHLs that turn this computation as a complex task.Generally, the generated coupling mechanism and electromagnetic fields are modeled through the differential equations.The dedicated assessment approaches, such as finite difference method (FDM) [1] or the finite element method (FEM) [16] are commonly used to solve the differential equations.The EMI problems are transferred to a pure numerical simulation by FEM [17].However, the computation burden will exponentially increase by expanding dimensions of simulation geometry and simulation complexity as a huge number of extended meshes and parameters have to be assessed [18].A new mesh discretization and new evaluation has to be considered for each new simulation of the problem geometry.Nonetheless, immense computational power and time are required to profoundly evaluate the EMI between high voltage OHLs and metallic pipelines for various system arrangements via FEM [19].Hence, any prediction methodologies capable of accurately estimating the requested information from a specified set of problem configurations recently got more attention [20].
The artificial intelligence (AI) based techniques have been considered a superior solution for the computation of AC induced voltages [17,20].AI based approaches often warranty a sufficient degree of prediction accuracy for complex system modeling [21].
Among different AI-based techniques artificial neural networks (ANNs) are the most broadly used approaches for modeling the induced voltages [22].ANNs are of interest for predicting the induced voltages on metallic pipelines placed in high voltage OHLs RoW due to their capabilities such as precise pattern learning, a memory to recall information from past experience, determine relations between dependent (output) and independent (input) variables and discover diverse discriminators in the complex system [1,20,[23][24][25].The ANN can be applied for predicting induced voltage on metallic pipelines on account of its ability for handling noisy information, having memory and parallel calculation architecture.The accuracy and effectiveness of ANNs approaches are strongly associated with convergence speed, neural network architecture, and the weight updating algorithm [26].The multilayer perceptron (MLP) with the error back propagation training method as an effective class of feedforward ANN was implemented in [27], to predict the level of induced voltages on a buried metallic pipeline, located in the electromagnetic field caused by a high voltage OHL during single phase to ground fault conditions.The resistivity of the ground, the distance between the OHL and buried metallic pipeline, the magnitude of fault current, and the connection of the pipeline to the mitigation system are considered as input for ANN in this study.Additionally, the induced voltage on the buried metallic pipeline for a range of fault current and separation distances is directly estimated by the applied method.
Support vector regression (SVR) is one more AI-based method, which has been widely used for complex system modeling as an enhanced predictive method because of its capability to easily learn and adapt to the complex patterns [28].The main superiority of SVR in comparison with ANN is that the global optimum in the training phase of SVR is always funded.Moreover, in comparison with ANN, SVR has a lesser tendency for over fitting and it has plainer geometric interpretation.It also provides more sparse results [29].
Due to the application feasibility and simplicity for the development of the hardware, the fuzzy logic system (FLS) has been applied in various power system projects, and industrial processes.The FLS was effectively used in [28] for the calculation of the level of magnetic fields nearby the an OHL during a phase to ground fault condition.A major drawback of this approach is high complexity of the gradient approach used to find the optimum parameters of fuzzy logic rules.Then, the genetic algorithm (GA) optimization method is used to determine the optimum parameters of the requested fuzzy logic rules [30].Still, the FLSs are based on predetermined "if/then" rules that disrupt the capability of fuzzy systems to adapt and learn from different conditions.To overcome this shortcoming and offer a universal estimator approach with a high capability for the computation of induced voltages, hybridization of ANN with FLS as a neuro-fuzzy system has been proposed by the authors in [26].The adaptive neuro-fuzzy inference system (ANFIS) has a higher degree of generalization and consistency.It provides great prediction accuracy throughout various range of input data.
The optimized neuro-fuzzy based method is developed in this work as an enhanced version of ANFIS to further improve the prediction precision of the induced voltage on mitigated and unmitigated pipelines.The backtracking search algorithm (BSA) as an efficient optimization algorithm is applied in the learning process of ANFIS to promote the prediction accuracy by tuning the membership functions for achieving a lower prediction error.
BSA is a recently developed optimization algorithm that not only delivers highly accurate solutions but also has quite a simple mechanism with only one control parameter [31].The performance of the proposed hybrid approach was examined attentively in comparison with the results from ANN, and SVR and ANFIS models optimized by two other robust and efficient metaheuristic optimization algorithms; namely, the cuckoo search optimization algorithm (CSA) and the particle swarm optimization (PSO) algorithm.Moreover, the obtained results are further validated by the observed data from the Current Distribution, Electromagnetic Interference, Grounding and Soil Structure Analysis (CDEGS) program developed by Safe Engineering Services & technologies (SES) company [32].
This paper is organized as follows: Section 3 provides a brief description of ANFIS.Then, the main principle of backtracking search algorithm (BSA) is presented in Section 4. Section 6 presents all the simulation results together with the performance of developed hybrid ANFIS-BSA method in comparison with other applied AI-based approaches.The statistical analysis is also provided in this section to evaluate the robustness of the proposed hybrid ANFIS-BSA method for predicting the induced AC voltage in mitigated and unmitigated pipelines placed in electromagnetic fields created by an OHL in the case of a fault.Finally, the conclusion was drawn by consolidating the important features of this study.

Adaptive Neuro-Fuzzy Inference System (ANFIS)
The fuzzy logic method is based on the predetermined "if/then" rules that lack the capability of this system to adapt and learn from new conditions.Thus, to overcome this drawback, authors in [33] hybridized a ANN with a fuzzy inference system (FIS) to form the ANFIS approach.The ANFIS approach is categorized as an enhanced adaptive system like ANN, that, through training, adapts the antecedent (fuzzy membership functions) parameters and the consequent (fuzzy system output function) parameters.
The ANFIS gains the advantages of both ANN and FIS and it does not suffer from the drawbacks associated with both methods.The complicated procedures of ANNs are solved by using linguistic variables of FIS system, and the disadvantage of FIS is bypassed by using the neural inference system, which provides the capability of ANFIS to adapt and learn from new conditions.Therefore, the ANFIS approach has the ability to simulate complex systems using ANN learning with FIS, and it has been classified as a universal estimator capable of predicting the output of complex systems.
ANFIS has been developed as an enhanced adaptive system with group of "if/then" fuzzy rules and tunable membership function (MF) parameters in the training phase.The antecedent (fuzzy membership functions), and the consequent (fuzzy system output function) parameters are two different sets that are optimized during the training phase to provide the learning procedures for ANFIS approach.
ANFIS comprises five different consecutive layers: the first layer is the if-part (fuzzification); the second layer is the rules (production) part; the third layer is norm part; the fourth layer is the then-part (defuzzification); and the fifth layer is the output part [34].The main ANFIS structure with two inputs (i.e., x and y) as independent variables and one output (i.e., f out ) as the dependent variable is illustrated by Figure 1.Sugeno and Mamdani are two different types of FISs.These two types of FISs have different defuzzification procedures and consequences in "if/then" fuzzy rule sets.
ANFIS architectures represent both the Mamdani and Sugeno methods.In contrast to Sugenotype FIS, Mamdani has less flexibility to be integrated with the ANFIS approach to precisely model the complex systems [35].
The "if/then" rules of ANFIS approach integrated with first order Sugeno FIS are defined as: where fi (x, y; pi, qi, ri) is a polynomial function providing the output of fist order Sugeno FIS.In this function, the two inputs of the ANFIS model are represented by x and y, and the output of the ANFIS model is indicated by z.The fuzzy sets are represented by Ai and Bi.Generally, different node functions constrict the layers of ANFIS approach.As illustrated in Figure 1, the adaptive node, the adjustable variable, is indicated by a square and the fixed node, fixed parameter, is indicated by a circle.

• First layer (if-part):
The first layer consists only of adaptive nodes, as represented by the flowing function: , =  (),  = 3,4 The inputs to the node i are represented by x and y.Ai and Bi are different linguistic labels.The MF for Ai and Bi fuzzy sets are represented by  and  , respectively.The membership degree of a fuzzy set is represented by Q1,i .The node i output determines the grade to which each x or y input complies with the quantifiers.
Generally, any form of typical MF can be implemented in ANFIS approach.
The Gaussian membership function is the most widely used MF that is specified as follows: where  and c determined the width and center of the Guassian membership function, respectively.
• Second layer (rules): The second layer consists of only fixed nodes.The products of all connected signals to the node is considered the output of each node in this layer.This layer determines the firing strength for each fuzzy rule via multiplication of incoming signals as follows: where the firing strength for fuzzy rule is represented by output signal of wi.
• Third layer (normalization): The third layer consists of only fixed nodes.This layer normalizes the firing strength calculated in the second layer, through computing the ratio of the i th firing strength of the fuzzy rule to sum of all the firing strengths of fuzzy rules.Sugeno and Mamdani are two different types of FISs.These two types of FISs have different defuzzification procedures and consequences in "if/then" fuzzy rule sets.
ANFIS architectures represent both the Mamdani and Sugeno methods.In contrast to Sugeno-type FIS, Mamdani has less flexibility to be integrated with the ANFIS approach to precisely model the complex systems [35].
The "if/then" rules of ANFIS approach integrated with first order Sugeno FIS are defined as: where f i (x, y; p i , q i , r i ) is a polynomial function providing the output of fist order Sugeno FIS.In this function, the two inputs of the ANFIS model are represented by x and y, and the output of the ANFIS model is indicated by z.The fuzzy sets are represented by A i and B i .Generally, different node functions constrict the layers of ANFIS approach.As illustrated in Figure 1, the adaptive node, the adjustable variable, is indicated by a square and the fixed node, fixed parameter, is indicated by a circle.

• First layer (if-part):
The first layer consists only of adaptive nodes, as represented by the flowing function: The inputs to the node i are represented by x and y.A i and B i are different linguistic labels.The MF for A i and B i fuzzy sets are represented by µ Ai and µ Bi , respectively.The membership degree of a fuzzy set is represented by Q 1,i.The node i output determines the grade to which each x or y input complies with the quantifiers.
Generally, any form of typical MF can be implemented in ANFIS approach.
The Gaussian membership function is the most widely used MF that is specified as follows: where σ and c determined the width and center of the Guassian membership function, respectively.
• Second layer (rules): The second layer consists of only fixed nodes.The products of all connected signals to the node is considered the output of each node in this layer.This layer determines the firing strength for each fuzzy rule via multiplication of incoming signals as follows: Energies 2019, 12, 3651 6 of 18 where the firing strength for fuzzy rule is represented by output signal of w i .
• Third layer (normalization): The third layer consists of only fixed nodes.This layer normalizes the firing strength calculated in the second layer, through computing the ratio of the ith firing strength of the fuzzy rule to sum of all the firing strengths of fuzzy rules.
where the normalized firing strength of a fuzzy rule as output signal is represented by w.
• Fourth layer (then part): In the fourth layer a node function is adaptive with each node i.
where the "if/then" fuzzy rules are defined by the following f 1 and f 2 functions: Rule1: If x is A 1 and y is B 1 then z=f 1 (x, y; p 1 , q 1 , r 1 ) Rule2: If x is A 2 and y is B 2 then z=f 2 (x, y; p 2 , q 2 , r 2 ) where the parameter set known as the consequent parameters are represented by r i , q i , and p i .
• Fifth layer (output): The fifth layer consists of only a single fixed node that calculates the ANFIS output through adding all the arriving signals.
The total output of the ANFIS approach is a summation of all consequent signals.Hence, the total output of the ANFIS approach is defended as follows: Finally, a hybrid learning algorithm was applied with the ANFIS approach to accurately tune the parameters.

Backtracking Search Algorithm (BSA)
BSA, as one of the newly developed metaheuristic optimization algorithms, has relatively simple constriction, but is an efficient algorithm for finding the optimum solutions in non-convex and multimodal optimization problems.Due to the quite simple structure of BSA, it has been applied to solve various engineering optimization problems.During the evaluation of BSA, a trial population was generated through two enhanced crossover and mutation operators.The structures of enhanced mutation and crossover in BSA are extremely different, with the mutation and crossover operators specified in former optimization algorithms, such as differential evolution (DE) and GA.Those two enhanced operators are applied in this algorithm to provide a balance between exploitation of optimal result and exploration of the optimization search space.BSA has efficient exploitation and exploration capabilities because of several effective mechanisms defined in this optimization algorithm to generate a trial population, control the boundaries of optimization search space and adapt the magnitude of the search direction [34].
Energies 2019, 12, 3651 7 of 18 BSA has only one control parameter and the performance of this optimization in not excessively sensitive to the preliminary value of this parameter.It does not suffer from drawbacks associated with other metaheuristic methods, such as being computationally expensive, being trapped in local optima, having various control parameters, and being overly sensitive to the initial values of these parameters.Additionally, a track of previous generations is stored in memory to share the past experiences for generating a trial population.It stores a randomly chosen population from previous generation in its memory for the generation of the matrix of the search-direction.The superiority of BSA for finding optimum solutions in numerical optimization problems over different widely used metaheuristic optimization algorithms is verified through the statistical analysis reported in [31].Generally, six different steps construct the BSA structure as coded in Figure S1.The steps of BSA are as follows: Step 1: Initialization Scattering the population members in the solution space (Equation ( 11)) where: nPop is population size.nVar signifies the optimization variable.Uniform distribution function is U. low j and up j are upper and lower search space limits of j th variable.y i is productivity of i th individual.g is generation number.
Step 2: Selection-I (1) Initializing a historical population (old P) to determine the search-direction matrix (Equation ( 12)); (2) Redefining the historical population at each iteration based on (if/then) rule by comparing two random numbers, a and b.Subsequently, population (P) pursues old P until it is changed to provide a memorization process and facilitates the exploration search space (Equation ( 13)); (3) At the end of step 2, a hierarchical sequence has been permuted by shuffling random function (Equation ( 14)); where OldP is a historical population where, := is the updated operation.a and b are randomly generated numbers.

System Modeling
The system modeled in this research is provided in Figure 2. It is comprised of 132 kV OHL and a neighboring, well coated 16" pipeline.The total length of the metallic pipeline is 10 km, whereas the high voltage OHL length is 20 km.The length of the parallelism of OHL with the pipeline is 10 km; the metallic pipeline is located at the central site with a burial depth of half meter, connected to the gradient control wires as mitigation system.The gradient control wire, as an effective type of mitigation system, is a horizontally-buried zinc anode wire placed at the bottom of the trench beside the metallic pipeline [13], and regularly connected to the metallic pipeline, as illustrated in Figure 2. The gradient control wires typically connect to the pipeline at intervals varying from 150 to 600 m.The DC decoupling devices are connected between the pipeline and gradient control wire, which will decouple the DC while remaining the AC coupled.The function of this device is providing a path for the AC current to flow from the pipeline to the zinc ribbon grounding system, while blocking the DC cathodic protection current from flowing to the zinc ribbon conductor.

System Modeling
The system modeled in this research is provided in Figure 2. It is comprised of 132 kV OHL and a neighboring, well coated 16″ pipeline.The total length of the metallic pipeline is 10 km, whereas the high voltage OHL length is 20 km.The length of the parallelism of OHL with the pipeline is 10 km; the metallic pipeline is located at the central site with a burial depth of half meter, connected to the gradient control wires as mitigation system.The gradient control wire, as an effective type of mitigation system, is a horizontally-buried zinc anode wire placed at the bottom of the trench beside the metallic pipeline [13], and regularly connected to the metallic pipeline, as illustrated in Figure 2. The gradient control wires typically connect to the pipeline at intervals varying from 150 to 600 m.The DC decoupling devices are connected between the pipeline and gradient control wire, which will decouple the DC while remaining the AC coupled.The function of this device is providing a path for the AC current to flow from the pipeline to the zinc ribbon grounding system, while blocking the DC cathodic protection current from flowing to the zinc ribbon conductor.

Simulation Results and Discussion
In this paper, ANFIS-BSA is applied to improve the accuracy of induced voltage prediction on mitigated and unmitigated pipelines.The level of the induced pipelines voltage, depends on diverse factors, such as the level of OHL voltage, parallelism length, the separation distance, resistivity of the ground, the load current magnitude, OHL configuration, and the material of pipeline coating and the mitigation system.The level of fault current, the separation distance, mitigation system, and the resistivity of the ground have the higher impacts on the induced pipeline voltage, as reported in [26].Generally, the single phase to ground fault condition is considered to evaluate the induced voltage on Energies 2019, 12, 3651 9 of 18 neighboring pipelines as it is the most frequent fault in the distribution and transmission power system and it produces very high level of induced voltage on the metallic pipelines.This can be explained by the fact that within the two or three phase-to-ground fault occurrences, cancellation of the magnetic field will occur as a result of a poorer induced voltage generated on the metallic pipeline.Hence, in this study, the maximum induced voltage on the pipeline is calculated for mitigated and unmitigated pipelines with varying separation distances, various resistivities of ground, and a diverse level of single line to ground fault currents.The parameters of the system used in this study are tabulated in Figure 2.This data set of input variables covers diverse sets of possible situations.
The validity of the obtained models is verified by comparing the simulation results with the corresponding results obtained from the CDEGS software.Furthermore, to assess the effectiveness of ANFIS-BSA for predicting the induced pipeline voltage, its performances are compared with the following methods: MLP, SVR, ANFIS, MLP-PSO, MLP-CSA, MLP-BSA, SVR-PSO, SVR-CSA, SVR-BSA, ANFIS-PSO, and ANFIS-CSA.
The procedure of the proposed methodology for predicting the induced voltage on both mitigated and unmitigated pipelines is illustrated in Figure 3. Commonly, for finding the practical models generated by the machine learning algorithms, a similar process is carried out.Hence, for obtaining the optimal AI based models to predict the induced voltage on the metallic pipelines, the following serial steps are followed for all applied approaches in this study.

1.
The independent variables consisted of four inputs representing fault current, soil resistivity, separation distance, and mitigation system, while the dependent variable represented the total pipeline's maximum voltage.

2.
Both dependent and independent variables were randomly distributed into two different phases: fifty of the total sixty-five systems with different configurations as the training phase, and the reminder as the testing phase.Since the range of dependent and independent variables varies widely, both variables were normalized by Equation (19).To speed up the learning process, the observed data were normalized prior to data processing.The main purpose of raw data normalization was unifying the observed data into a common scale. __ where Z is the data to be normalized and __ Z is the normalized data, and t is the number of observations.

3.
The learning process occurred during the training phase.The computer programs that link the input variables to the output were developed during learning process.The data required for the training of the AI-based methods was obtained via the CDEGS program.This program is especially designed to automate and simplify the modeling of complex RoW arrangements involving power transmission lines and other utilities, such as water, oil, or gas pipelines.Its results were strongly validated by analytical equations and by an experimental test rig reported in [8,36].Although the testing phase does not have any role in developing the models, it was employed to assess the performance of the models obtained by AI-based methods.To measure the predictive accuracy of the generated models, several evaluation criteria were used, such as Thiel's inequality coefficient (U-statistic), root mean square error (RMSE), absolute error, and mean absolute percentage error (MAPE).The mathematical equations of those criteria are as follows: where the U-statistic provides a measure of how well fitted a time series of predicted values to a corresponding time series of observed data.The U-statistic is always in the range of zero to one, with a value closer to one indicating the estimation is no better than a naive estimate and the value closer to zero demonstrates higher prediction accuracy with a great fit.4.
The Durbin-Watson (whiteness) test was calculated to guarantee that the generated models sufficiently describe given data sets [37].The whiteness test is calculated via a confirmatory analysis.The main purpose of confirmation analysis is to guarantee the whiteness of estimated residuals.The whiteness of estimated residuals (e(t)) indicates that they are uncorrelated.
The residuals autocorrelation function (RACF) is used to study the correlation of the whiteness of estimated residuals through the following equation: Since there is no conformity in the optimal values of the machine learning algorithms' parameters, setting the control parameters of all applied AI-based approaches was done according to the successful approaches in the literature.All parameter settings of applied methods are summarized in Table 1.The optimization methods used the learning process of MLP according to [38] and the optimized SVR approaches were developed based on [39].
The performances of applied machine learning methods for predicting the total induced voltage The RACF values come into the range of zero to one; if the value of RACF is meaningfully diverse from zero, it will fall outside a confidence level.This specifies that the residuals are not correlated (white) and an important input (independent) variable has been missed in the tested model.
Since there is no conformity in the optimal values of the machine learning algorithms' parameters, setting the control parameters of all applied AI-based approaches was done according to the successful approaches in the literature.All parameter settings of applied methods are summarized in Table 1.The optimization methods used the learning process of MLP according to [38] and the optimized SVR approaches were developed based on [39].The performances of applied machine learning methods for predicting the total induced voltage on unmitigated pipeline are tabulated in Table 2.The RACF values in these tables verify the whiteness of estimated residuals within a confidence interval for all developed models.The prediction accuracy of the studied methods in terms of multi-criteria decisions, using the mean rank of the methods for each indicator (MAPE, U-statistic, absolute error, and RMSE) in whole set is ranked: ANFIS The comparison between the accuracy of methods applied for the prediction the induced voltage on an unmitigated pipeline reveals that optimized ANFIS approaches outperform other studied methods.Furthermore, it was found that the most efficient optimization algorithm applied for training ANFIS was BSA, as the superior MAPE (1.1581%), U-statistic (0.0072), RMSE (0.0197), and absolute error (1.0213) values reported in Table 2 belong to ANFIS-BSA.The performance of proposed method (ANFIS-BSA) in testing and training phases is depicted in Figures 4 and 5.
For further examination of solution methodology, the performance of optimized ANFIS for predicting the total induced voltage on mitigated pipeline is compared with those from other machine learning methods, as shown in Table 3.The computed RACF values in this table specify that the estimated residuals of all developed models are uncorrelated and attained models satisfactorily describe the given data series.Table 3.Comparison between the accuracy of the applied approached for predicting the total induced voltage on a mitigated pipeline.
of the electromagnetic field will be reduced; consequently, the pipeline's voltage level will be decreased.It also can be seen that in the small separation between the OHL and the pipeline, the effect of fault current on the total pipeline potential is more severe than the effect of soil resistivity.Moreover, the obtained results confirmed that mitigation systems can be used in order to significantly reduce the induced pipeline voltage to acceptable level for various fault conditions.

Conclusion
In this work, a hybrid approach composed of ANFIS and BSA has been developed for predicting the voltage on mitigated and unmitigated metallic pipelines built in OHLs' right-of-way.The actual parameters of the system, the mitigation system, fault current, resistivity of the ground, and separation distance, are taken as the inputs (independent variables) and the pipeline voltages predicted are the output (dependent variable).Regardless of which practical system is investigated, the obtained results indicated that the developed models by ANFIS-BSA provide enhanced estimations over other solutions from AI-based methodologies.Total induced voltage on the pipeline

Conclusions
In this work, a hybrid approach composed of ANFIS and BSA has been developed for predicting the voltage on mitigated and unmitigated metallic pipelines built in OHLs' right-of-way.The actual parameters of the system, the mitigation system, fault current, resistivity of the ground, and separation distance, are taken as the inputs (independent variables) and the pipeline voltages predicted are the output (dependent variable).Regardless of which practical system is investigated, the obtained results indicated that the developed models by ANFIS-BSA provide enhanced estimations over other solutions from AI-based methodologies.Total induced voltage on the pipeline has been modeled by ANFIS, SVR, ANN, and an optimized version of these approaches.According to the obtained results, the ANFIS-BSA models consistently give superior predictions to other studied AI-based methods.
Excellent agreement between obtained data from CDGES program and ANFIS optimized by BSA has been obtained.The relative importance of the independent variables has been also studied.This paper clearly specified the potential of the applied methodology for capturing the interactions between the levels of pipeline voltage and the independent variables, and even for the investigation of the relative importance of these input variables.Additionally, the results demonstrate that developed approach can accurately predict the total pipeline voltage for a wide range of input variables, with and without a mitigation system.The accuracy of the predicted pipeline potential is important for designing mitigation systems that will increase the overall pipeline integrity and safety.The developed methodology can be plainly applied for diverse environmental, dimensional, and power system conditions.

Figure 2 .
Figure 2. Cross section of 132-kV overhead lines (OHL) and the buried pipeline with a mitigation system.

Figure 2 .
Figure 2. Cross section of 132-kV overhead lines (OHL) and the buried pipeline with a mitigation system.

Figure 4 .
Figure 4. Performance of proposed method (ANFIS-BSA) during testing and training phases for predicting the total induced voltage on unmitigated pipeline.

Figure 4 .
Figure 4. Performance of proposed method (ANFIS-BSA) during testing and training phases for predicting the total induced voltage on unmitigated pipeline.

Figure 4 .
Figure 4. Performance of proposed method (ANFIS-BSA) during testing and training phases for predicting the total induced voltage on unmitigated pipeline.

Table 1 .
Parameter settings of applied methods.

Table 2 .
Comparison between the accuracy of the applied approached for predicting the total induced voltage on an unmitigated pipeline.