An Optimized Protection Coordination Scheme for the Optimal Coordination of Overcurrent Relays Using a Nature-Inspired Root Tree Algorithm

In electrical engineering problems, bioand nature-inspired optimization techniques are valuable ways to minimize or maximize an objective function. We use the root tree algorithm (RTO), inspired by the random movement of roots, to search for the global optimum, in order to best solve the problem of overcurrent relays (OCRs). It is a complex and highly linear constrained optimization problem. In this problem, we have one type of design variable, time multiplier settings (TMSs), for each relay in the circuit. The objective function is to minimize the total operating time of all the primary relays to avoid excessive interruptions. In this paper, three case studies have been considered. From the simulation results, it has been observed that the RTO with certain parameter settings operates better compared to the other up-to-date algorithms.


Introduction
The power system consists of three main parts: generation, transmission, and distribution of electrical energy and works at voltage levels ranging from 415 V to 400 KV or greater.Moreover, the supply lines, which transmit electrical power, are not insulated.These lines show irregularities more often than other parts of the power system due to several issues, such as lightning, which lead to overcurrents.The imbalance due to these irregularities interrupts the behavior of power and, moreover, leads to the impairment of the other accessories linked to the power system.In order to eliminate these problems, protective measures should be taken into account.To improve this difficulty, overcurrent relays (OCRs) have been commonly used for years as a safety measure in the power system to avoid mal-operation in the power supply.Thus, OCRs are easy to use for the safety of the sub-transmission power system and a secondary layer of backup protection in transmission systems.During the designing of the electrical power system, it is important to consider the coordination of these OCRs.In case any fault occurs in the power system, OCRs help the breaker to play the logical elements and these relays are placed at both ends of the line.The OCRs decide which part of the power system has to come into action during a fault so that the faulty section is isolated and does not interrupt the power system, with certain constraints like adequate coordination tolerance and without extra disruptions.This procedure is mainly based on the nature of relays and other protection measures.The design variable of time dial settings (TDSs) for each relay in the circuit needs to be optimized.systems, which are introduced in this paper.We have actualized RTO to take care of the issue of OCR settings, and the results of our simulation are contrasted with other up-to-date algorithms.
The rest of the paper is organized as follows.In Section 2, we elaborate on the mathematical formulation of the problem.The RTO is recalled in Section 3. Parameter settings and statistical and graphical results are discussed in Section 4. Section 5 concludes the present study.

Problem Formulation of the Overcurrent Relay
The coordination of the overcurrent relay is defined as an optimization issue in a multi-or single-source loop system.Nonetheless, the coordination issue has an objective function and limitations that should fulfill the distinct constraints: where the parameters w j and T j,k are the weight and operating time of the relay, respectively.For all the relays, w j = 1 [33].Therefore, the characteristic curve for operating relay R i can be chosen from a portion of the selectable decision of IEC norms and could be characterized as: where α and k are steady parameters which characterize the relay characteristics and are expected to be α = 0.14 and k = 0.02 for normal inverse type relay.The factors TMS i and I p j are the time multiplier setting and pickup current of the ith relay, respectively, while I f i is the fault current flowing through relay R i .
where PSM remains for the plug setting multiplier and Ip j is the primary or main pickup current.
The above issue mentioned in condition ( 4) is a nonlinear issue in nature.The coordination can be planned as linear programming by considering the plug setting of the relay and the working time of the relays, a linear function of TMS.In linear programming, the TMS is ceaseless while rest of the parameters are steady, so condition (4) moves toward becoming: where: Hence, the objective function can be formulated as:

Constraints
The total operating time could be minimized under two kinds of constraints, including the constraints of the relay parameter and constraints of coordination.The first constraints contain the limits of TMS while the other constraints are appurtenant to the coordination of the main and secondary relays.The limit on the relay setting parameters necessitates constraints in light of the decisions of relay parameter and design, and the points of confinement can be communicated as follows: However, for a reasonable security margin, the pickup current should be lower than the short-circuit current and should be greater than the highest load current simultaneously.The other constraints are convenient for the acclimation of primary and secondary relay operating time.The coordination should be kept in a proper way where the secondary operating time of the relay should be greater than that of the primary relay.If the primary relay fails to clear the fault, the secondary relay should initiate its operation by opening the circuit breaker within a certain interval known as the coordination time interval (CTI).The coordination topology is shown in Figure 1.When a fault F takes place, it is sensed by both relays R i and R j .The primary relay responds first, having less operating time than that of the backup relay.For fault F, beyond bus bar i, relay R i should respond first to clear the fault.The operating time of R i is set to 0.1 s.The backup relay should wait for 0.1 s plus the operating time of the circuit breaker at bus bar i and overshoot time of the relay R j [34].To maintain selectivity of primary and backup protection and to keep the grading of relay coordination accurate, the operating time of the backup relay (R j ) should always be greater than that of primary relay (R i ) by an amount in which R i was supposed to have cleared the fault.This time is usually the sum of operating time of R i and the breaker operating time.The current and voltage wave for the primary/back up pair are shown in Figure 2 for the typical single-end ring main feeder used in the distribution system.Hence, the coordination constraints could then be defined as: where the parameter T i and T j are primary and backup relay operating times, respectively.

Root Tree Optimization Algorithm
The rooted tree optimization technique is a natural and biological algorithm inspired by the random-oriented movement of roots developed by Labbi Yacine et al. in 2016 [35].Roots work in a group to find the global optima and the best place to get water, instead of working individually.However, an individual root has a limited capacity and a greater number of roots are based around the place which links the plant with the source of water.To design the algorithm, an imaginary nature of roots should be taken into consideration for their combined decision related to the wetness degree, where the head of the root is located.To find one or more wetness locations by random movement, these roots call other roots to strengthen their presence around that position to become a new starting point for majority of root groups to get to the original place of water that will be the optimal solution.Figure 3 shows the behavior of the roots of a plant that how they search for water to find the best location.The roots which are far or have less wetness degree are replaced by new roots which are oriented randomly.However, the roots which have a greater amount of wetness ratio will preserve their orientation, whereas roots (i.e., solutions) that are far from the water location could be replaced by roots near the best roots of the previous generation.However, the completion of a maximum number of cycles or generations is the stopping criterion at the end of generations.The solution with best fitness value will be the desired solution.The evaluation of population members is based on a given objective function which is assigned with its fitness value.The candidate with the best solutions is forwarded to next generation while the others are neglected and reimbursed by new group of random solutions in each iteration/generation.The proposed algorithm starts its work by creating an initial population randomly.However, for the RTO algorithm, there are some important parameters that need to be defined regarding how roots start the random movement from the initial population to new population.Those parameters are roots and wetness degree (wd).These two parameters can give the suggested solution and fitness value to the rest of the population.In this work, the RTO algorithm is used to find the optimal solution for the relay coordination problem within a power system.The suggested procedure has an extraordinary exploration ability and merging promptness in contrast with other methodology.This distinguishing property makes the populace individual from RTO more discriminative in locating the ideal result compared to that of another developmental method.The primary point of this paper is to locate the ideal estimations of TMS, keeping in mind the end goal to limit the aggregate working time of OCRs under a few requirements like relay settings and backup constraints.

The Rate of the Root (R n ) Nearest to the Water
The term rate here represents those roots which are a member of the total population that gather around the wetter place.This root will be next in line to that root where wetness is weaker than that of the previous generation.The fresh population of the nearest root according to wetness condition can be defined as: where it is the current step of the iteration, x n (i, it + 1) is the fresh member for the next iteration, i.e., (it + 1).The best solution achieved from the previous generation can be represented by x it b , while the parameters k 1 , N, i, and l represent the adjustable parameter, population scale number, and upper limit, respectively, and randn is a normal random number having a value between −1 and 1.

The Rate of the Continuous Root (R c ) in Its District
It is the root which has greater wetness ratio and moves forward from the previous generation.The new population of the random root is expressed as follows: whereas k 2 is the adjustable parameter while x(it) is the iteration for the previous candidate at iteration it and randn is a random number that has a range between 0 and 1.

Random Root Rate (R r )
In this case, the roots spread randomly to find the best location of water to achieve the maximum rate of getting the global solution.The roots with less ratio of wetness degree are also replaced from the last generation.In this step, a new population calculated for a random root could be expressed as: where the parameter x r is the candidate which is selected randomly from the previous generation with an adjustable parameter k 3 .

The Step Root Tree Algorithm
The steps of this algorithm can be compiled as follows: Step 1: In this step, the initial generation is created randomly and is composed of N members with the variable limits in the research location with induce numerical values of R n , R r , and R c rates.
Step 2: In the second step, all population members are measured based on their respective wd for the maximum and minimum objective function: For the maximum objective Step 3: The new population is created and replacement of the member is done according to the wetness ratio of R n , R r , and R c .Equations ( 10)-( 12) are used for a candidate having the smallest value until the one with the same wetness ratio is achieved.
Step 4: In this step, the meeting criteria are satisfied to display an optimal result with the best fitness value.Otherwise there is a return to Step 2.

Implementation of the RTO Algorithm for the OCR Coordination Problem
Step 1: Define the input parameter that includes TMS, the maximum number of iterations, population size, and the adjustable parameter (i.e., different rate values of R n , R r , and R c ).
Step 2: In this step, the population is initially generated with respect to the relay bounding.These individuals must be a feasible candidate solution that satisfies the relay operating constraints.
Step 3: Measure the fitness value of each candidate x b in the population, with the help of Equation (10).Additionally, an evaluation of wetness ratio is done in this step for each candidate.
Step 4: In this step, the comparison is done between the best fitness values of candidates x b best of the individual population.
where X bi,r is the individual that is randomly selected from the current population, n is the population size, X bi,max , and X bi,max are upper and lower parameter limits, respectively, and i is the number of iterations.
Step 6: If the number of iterations reaches the maximum, then go to Step 7. Otherwise, go to Step 3.
Step 7: The candidates that generate the latest x b best is the optimal TMS of each relay with the minimum total operating time that is the objective function.
The implementation and pseudocode of RTO for solving the coordination problem of OCR is depicted in Figure 4 and Algorithm 1.
Algorithm 1.The structure of the new rooted tree algorithm [35].
Algorithm RTO Begin // Initialization: Set the rates R n , R r and R c parameters Give the maximum number of iterations: MaxIte, and the population scale: the RTO size Set iteration counter it = 1 Generate the initial population X(1) randomly within the search range of (Xmin, Xmax) // Loop Repeat Evaluate the wdi for each root // wd: Fitness, root: Individual Reorder the population according to the witness degree Identify the candidate according to the wetness place Xbest // the global best in the whole population

Results and Discussion
In this segment, a legitimate program has been created in MATLAB software to locate the optimum value of OCR in a single-and multi-loop distribution system utilizing RTO.The productivity and execution of RTO were tried for the distinctive single-and multi-loop system, and it was discovered that the RTO gives the most tasteful and best solution in all the contextual analyses.Three contextual analyses have been viewed.The system details of all contextual investigations could be found in [36][37][38][39].
In each case study, the following RTO parameters were used.The comprehensive explanation of the problem formulation and application of RTO to find the optimal solution is presented for all case studies.For the purpose of this study, population size = 50 and the maximum number of iteration = 200.

Case Study 1
Figure 5 shows a multi-loop system having eight overcurrent relays.The configuration and combustion of the primary and secondary relay pair models depend upon the location of the fault current in different feeders.In this combination, six different fault locations are taken into consideration.The aggregate fault current and the essential reinforcement relationship of the relay for the six fault focuses are given in Table 1.All relays have a plug setting of 1 and a CT ratio of 100:1.Table 2 shows the current seen by relays and the a ρ constant for various faults.In this case, eight variables had been taken into consideration (i.e., the TMS of eight relays), eight constraints due to bounds on relay operating time (or bounds on the TMS of relays), and nine constraints due to coordination criteria.The duplex of the proposed problem will have 17 variables and 8 constraints.In order to validate the efficiency of RTO, the CTI ratio of 0.6 s and minimum operating time (MOT) of the relay was taken as 0.1 s.The TMS of all the eight relays are x1-x8.
The problem for optimization can be demonstrated as: minz = 28.975x 1 + 28.975x 2 + 37.736x 3 + 11.502x 4 + 3.297x 5 + 3.297x 6 + 4.9807x 7 + 30.7994x 8 (15) The constraints owing to the MOT of the relays are: Hence, the constraints mentioned in Equations ( 18), ( 19), (22), and ( 23) violate the constraints of the minimum value of the (TMS).Hence, these constraints are reconstructed as: x 4 ≥ 0.025 (25) x 7 ≥ 0.025 (26) x 8 ≥ 0.025 (27) The constraints owing to the coordination of relays with CTI taken as 0.6 are: As discussed in Section 2, to apply RTO to this problem, the objective function was first converted into an unconstrained optimization problem by integrating the relay constraint into the objective function.The constraints pertaining to the MOT of the relay have the ability to alarm the lower bounds of the RTO.The TMS values of all the relays varies from 0.025 to 1.2.The constraint owing to the MOT of the relay was taken care of defining the lower bounds of the variable in RTO.The lower and upper limit of all TMS of the relay is considered as 0.025 to 1.2.The constraints pertaining to coordination of the relay depicted by Equations ( 28)- (35) were incorporated in the objective function by means of a penalty scheme.In order to justify the RTO, a population size of 50 (i.e., candidate solutions, eight designed variables x 1 to x 8 , and 200 iterations as meeting criterion) has been assigned to the coordination problem.The population was conveyed to the fitness function and the proceeding values of the objective function are perceived.As the objective function is of the minimization type, the lowest fitness value of the objective function is considered is the best solution among the population size.After running 200 iterations, the optimum values are shown in Table 3. Table 3 shows the optimal TMS values achieved by the proposed RTO method for this case and has been compared with the literature.Figure 6 shows the objective function values achieved during the course of the simulation for the best candidate solution in each iteration, and this figure shows that the convergence rate is faster and achieves satisfactory results in less iterations compared to other techniques cited in the literature.This figure also shows that R n = 0.4, R r = 0.3, and R c = 0.3 are the feasible values for the RTO technique.These parameter values are used for the rest of the case studies presented in this paper as it gives an optimal solution compared with other adjustable parameters of the RTO algorithm.Figure 7 shows the graphical representation of the optimized total operating time, Top (z), and demonstrates that the total operating time is minimized up to the optimum value.The optimal value of objective is found to be 26.681s by considering a CTI of 0.6 s, which is obtained in a fewer number of simulations.The values depicted in Table 3 illustrate that the RTO technique gives an ideal and best-streamlined aggregate working time (i.e., total operating time) up to the ideal optimum value.The optimum values guarantee that the relay will operate in the minimum time to detect a fault at any location.All the optimal values obtained by RTO satisfy all the coordination constraints.

Comparison of RTO with the GA Algorithm
The results obtained by using the RTO algorithm are compared with the genetic algorithm, as given in Table 3.The RTO outperforms the genetic algorithm for obtaining the optimum TMS values and gives the advantage of 5.20198 s over the GA with the same initial conditions as supposed for this case study.The proposed RTO performs outstandingly over GA, gives an advantage in total net time gain, minimizes the total operating time up to the optimum value, and maintains proper coordination during a fault condition.In this condition, the results obtained by RTO for all the relays will satisfy the coordination constraints.Furthermore, no violation has been found regarding the coordination constraints.

Case Study 2
A multi-loop system, as appeared in Figure 8, with six overcurrent relays and with inconsequential line charging admittances are considered.Four faults are taken into consideration in the midst of the lines, i.e., A, B, C, and D.

Mathematical Modelling of the Problem Formulation
Table 4 demonstrates the line information of the given system.The primary/backup pair of relays for the system are shown in Table 5.For this delineation, four fault areas are taken, as shown in Figure 8.The CT proportions and plug setting are given in Table 6.For various fault areas, the current seen by relays and the a ρ constant are given in Table 7.In this representation, the entire number of constraints that arise are 11.Six constraints are attributable to limits on relay working time, and five constraints are by reason of coordination criteria.The MOT of each relay is 0.1, and the CTI is considered as 0.3 s.The TMSs of all six relays are x 1 -x 6 .
The problem for optimization can be demonstrated as: minz = 102.404x 1 + 6.0651x 2 + 98.758x 3 + 24.403x 4 + 35.319x 5 + 11.539x 6 (36) The constraints owing to the MOT of the relays are: 6.055x 2 ≥ 0.1 (38) 8.844x 4 ≥ 0.1 4.044x 5 ≥ 0.1 (41) Hence, the constraints mentioned in Equations ( 38)-( 42) violate the constraints of the minimum value of the TMS.Hence, these constraints are reconstructed as: x 3 ≥ 0.025 (44) x 5 ≥ 0.025 (46) The constraints owing to the coordination of relays with the CTI taken as 0.3 are: The objective function was evaluated using the proposed RTO calculation with same parameters elucidated from delineation 1.The optimal values of TMS gained are given in Table 8, which shows that the RTO algorithm gives ideal and best values and that it optimized the total operating time to optimal values.The optimal values ensure that the relay will take the least possible time in the system for a fault at any area.The time taken by Relay 1 to work is the base for a fault at point A and will set aside the most noteworthy time for a relay at point C.This is alluring in light of the fact that, for a fault at point A, relay 1 is first to work while, for a fault at point C, Relay 6 should get first chance to work if it fails.Relay 3 should take over tripping action and, if Relay 3 in a similar manner fails to work, then Relay 1 ought to expect control over the tripping action.The objective function value achieved over the span of the simulation for the best applicant arrangement in every cycle is shown in Figure 9.This demonstrates that the convergence is speedier and obtained a satisfactory value in fewer iterations.The optimum value of the objective function is observed to be 11.93 by considering CTI being taken as 0.3 s. Figure 10 demonstrates the optimized total operating time accomplished by the proposed algorithm with other methods mentioned in the literature

Comparison of RTO with Other Algorithms
To evaluate the execution of the proposed RTO algorithm, this technique was compared with the other mathematical and evolutionary techniques such as the CGA, FA, CFA, and CPSO technique accessible in the literature, as shown in Table 8.The RTO outmatched the CGA, FA, and CFA given the total net gain of times and gives advantages of 3.95, 4.32, and 2.46 s over these methods, respectively.However, in the case of CPSO, this technique almost converges to the same optimum solution as obtained by using the RTO with more computational effort.Furthermore, the RTO algorithm is better than the other techniques mentioned in the literature in terms of the nature of the results.Convergence representation and less number of iterations is required to get the ideal and best result.

Case Study 3
A parallel distribution system that is sustained from a singular end with seven overcurrent relays is shown in Figure 11.The CT ratio and PS of the relays are shown in Table 9.Four faults streams are constrained in the midst of the lines, i.e., A, B, C, and D. The primary/backup pair of relays for the given system are shown in Table 10.For this case, the total number of constraints is 12.Seven imperatives arise as a result of the points of confinement of the relay action, and five prerequisites arise due to the coordination condition.The MOT of each relay is 0.1 s.The CTI is considered as 0.2 s.The TMS values of all the relays is x 1 -x 7 .For different fault regions, the current seen by relays and the a ρ constant are given in Table 11.
The constraints owing to the MOT of the relays are: Hence, the constraints mentioned in Equations ( 55)-(60) violate the constraints of the minimum value of the TMS.Hence, these constraints are reconstructed as: x 3 ≥ 0.025 (62) x 4 ≥ 0.025 (63) x 5 ≥ 0.025 (64) x 6 ≥ 0.025 (65) The constraints owing to the coordination of relays with the CTI taken as 0. The objective function was solved by the proposed RTO with vague parameters.The ideal estimations of TMS and aggregate working time found are shown in Table 12.They similarly show the comparative delayed consequence of RTO with other metaheuristic and numerical procedures in the literature.In this depiction, no encroachment and miscoordination of the relays have been found.The time taken by relay R1 to work is the least for a fault at point C and is most extreme for a fault at point A. This is attractive on the grounds that, for a fault at point C, relay 1 is the first to work.Meanwhile, for a fault at point A, relay 7 ought to get the principal opportunity to work.On the off chance that it neglects to work, relay 4 should assume a control tripping activity.If relay 4 likewise neglects to work, relay 1 should assume control over the tripping activity at that point.The values are given in Table 12 ensure that the relays will work at any rate possible time for a fault whenever in the system.The objective value acquired over the span of the simulation for the best candidate course of action in each iteration is shown in Figure 12, which exhibits that the convergence is faster and get the best qualities in a lesser number of iterations.Additionally, the RTO outmatched the simplex method for this case and obtained an agreeable and better outcome compared to the simplex technique.The graphical image of the improved aggregate working time appears in Figure 13 and shows that the total operating time is streamlined up to the ideal and optimum value when compared to the simplex method, as mentioned in the literature.

Comparison of RTO with the Simplex Method
The results found using the RTO are differentiated, and the results procured by the simplex method are shown in Table 12.The RTO outmatch the simplex method in getting the optimal values of TMS and in an aggregate net gain of times.It gives a benefit of 10.5312 s over the simplex method, with the same starting parameters gathered for this contextual investigation.This gain in all out net time gain is adequate given that it is a small system.For all the ideal and optimal values for TMS, the RTO performed exceptionally compared to the simplex method, in terms of complete net time gain.Furthermore, it limits the aggregate working time up to optimal value and will keep up appropriate coordination even in a fault condition.For all the coordination conditions, the optimal value acquired by RTO for all the relays will fulfill the coordination limitations.Additionally, no infringement has been discovered with respect to the coordination requirements.

Conclusions
This paper suggested an RTO algorithm which simulates plant roots searching for water under the ground.The aforementioned RTO technique involved some parameters for tuning and, thus, is easy to implement.The coordination problem of the overcurrent relay utilizes the RTO algorithm for different test systems to assess the execution of the proposed RTO algorithm.The competence of the proposed RTO algorithm has been confirmed and tried through its application on various single-and multi-loop systems by analyzing its superiority and contrasting its execution with the CFA, FA, CPSO, CGA, GA, and the simplex method published algorithms.The simulation results of the RTO algorithm efficiently minimize all three models of the problem.The performance of the RTO can be seen from the optimized minimum objective function values and TMS, of each relay in the systems achieved by RTO for each case studies.In case 1 the objective function value is minimized up to optimum value by RTO and gives an advantage in total net gain in time of 5.20198 s over GA. while in case 2 and 3 the RTO gives a total net gain in time 3.95 s, 4.32 s, 2.46 s, and 10.5312 s over CGA, FA, CFA, and simplex method.The RTO algorithm gives a new seeking approach for elucidation, as one of its qualities is the extensive field of research attributable to the activities of the roots.The simulation results uncover the predominance of the proposed RTO algorithm in taking care of the overcurrent relay coordination problem.In future work, RTO will be applied to solve the problem of overcurrent relay coordination of higher and complex buses in the electrical power system.

Figure 1 .
Figure 1.A single-end radial distribution system.

Figure 2 .
Figure 2. Voltage and current waveform of the primary/backup pair of the relay for the single-end radial system.

Figure 3 .
Figure 3. Searching mechanism of the RTO algorithm for searching for the water location [35].

Figure 4 .
Figure 4.The flow chart of the RTO algorithm.

Figure 6 .
Figure 6.Convergence characteristic of Case 1 for different adjustable parameters to the RTO algorithm.

Figure 7 .
Figure 7. Optimized total operating time with the literature for Case 1.

Figure 8 .
Figure 8.A single end-loop distribution system.

Figure 10 .
Figure 10.Optimized total operating time with the literature for Case 2.

Figure 13 .
Figure 13.Optimized total operating time with the literature for Case 3.

Table 1 .
Total fault current and primary/backup relationships of relays for Case 1.
-Indicates no backup relay.

Table 2 .
a ρ constants and relay currents for Case 1.

Table 4 .
Line data for Case 2.

Table 5 .
Primary and backup relationships of the relays for Case 2.
-Indicates no backup relay.

Table 6 .
CT ratios and plug settings of the relays for Case 2.

Table 7 .
a ρ constants and relay currents for Case 2.

Table 9 .
CT ratios and plug settings of the relays for Case 3.

Table 10 .
Primary and backup relationships of the relays and maximum fault current for Case 3.