A Novel Application of Improved Marine Predators Algorithm and Particle Swarm Optimization for Solving the ORPD Problem

: The appropriate planning of electric power systems has a significant effect on the economic situation of countries. For the protection and reliability of the power system, the optimal reactive power dispatch (ORPD) problem is an essential issue. The ORPD is a non-linear, non-convex, and continuous or non-continuous optimization problem. Therefore, introducing a reliable optimizer is a challenging task to solve this optimization problem. This study proposes a robust and flexible optimization algorithm with the minimum adjustable parameters named Improved Marine Predators Algorithm and Particle Swarm Optimization (IMPAPSO) algorithm, for dealing with the non-linearity of ORPD. The IMPAPSO is evaluated using various test cases, including IEEE 30 bus, IEEE 57 bus, and IEEE 118 bus systems. An effectiveness of the proposed optimization algorithm was verified through a rigorous comparative study with other optimization methods. There was a noticeable enhancement in the electric power networks behavior when using the IMPAPSO method. Moreover, the IMPAPSO high convergence speed was an observed feature in a comparison with its peers. & A.A.: Funding D.Y.: and


Introduction
The electrical power system is a complex system consisting of subsystems, such as generation plants and transmission networks. Therefore, the performance of electric power systems needs to be controlled and evaluated in order to ensure proper and safe system operation. An important problem in power planning and operation is the Optimum Reactive Power Dispatch (ORPD). It has a great effect when studying how to operate the power system economically. Complex power is provided by the generation side, whereas active power is consumed by the demand side, and the reactive power keeps circulating through the network. However, the reactive power is important to keep the voltage within its stable limits. Additionally, reactive power is important for the transmission of the active power through the grid. Hence, the ORPD problem is an important aspect of power system operation.
Many parameters in the power systems are to be kept within certain limits, such as the voltage on the bus or power transmitted through a transmission line. Otherwise, any violation in these limits can lead to a system failure. As a result, these problems were solved by simulating and analyzing power systems, such as the optimal power flow (OPF) problem and the ORPD problem. In general, there are many objectives for the ORPD problem such as transmission cost minimization [1], active power loss minimization [2], improving the voltage profile [3], and voltage stability enhancement. Generally speaking, there are many variables that could be considered as control variables for the ORPD optimization problem, such as the voltages at the generating units, transformer tap settings, and the buses at which the shunt reactive power compensators could be placed. After solving this optimization problem, its constraints must be rechecked to assure that they are already kept within limits and there is no violation in the voltage or power values [4][5][6][7][8][9]. The optimal solution corresponds to the best fit control variables.
The ORPD problem was described through the use of many conventional optimization methods as non-linear and non-convex optimization problems [10][11][12][13][14][15]. Samples of these traditional tools are linear programming (LP) [10], interior point method [11], quadratic programming (QP) [12], and Newton-Raphson [13]. The traditional methods encountered many points that were considered to be weak points. They need many numerical iterations to finish the calculations. This, accordingly, leads to high computational time consumption. Lately, development is taking place and is continuing in the scope of the optimization tools. This resulted in the existence of many optimization techniques that were tested widely on various applications in the scope of optimization, especially in the field of power systems.
There are many meta-heuristic optimization methods used to optimally solve the ORPD. Sunflower optimization [16,17], hybrid firefly and particle swarm optimization [18], gravitational based search optimizer [19], brain storm optimization algorithm [20], tree seed algorithm [21], salp swarm optimization [22][23][24], sine-cosine algorithm [25], coyote optimization algorithm [26], seeker optimization algorithm [27], bacteria foraging optimization [28], cuckoo search algorithm [29][30][31], grasshopper optimization [32], Harris hawks optimization (HHO) [33], artificial bee colony (ABC) [34], gravitational search algorithm (GSA) [35], and the firefly algorithm [36], represent samples of the lately employed methods for the optimization problems. The genetic algorithm (GA) was utilized in [37], and the dynamic PSO was used in [38,39]. GA is fundamentally based on genetic dynamics concepts. The efficiency of such an approach depends mainly on adequate coding of the population. The first population is generated randomly at the beginning. Individuals in this population are seen as likely solutions to the problem of optimization. The progress of the simulation converts the population dynamics partly into future generations. In [40][41][42], GA was applied to the ORPD optimization problem. Differential evolution (DE) optimizer is also one of the optimization methods used in solving the ORPD. The advantage of the DE is that, it can be implemented easily and effectively [43]. Despite the effectiveness in solving the optimization problems, the optimal solutions might not be guaranteed every time they are used. One of the negative points of these strategies is that, they could be restricted to a minimal local. The fine tuning of the tool parameter also affects the convergence speed directly. The no-free-lunch theorem of the mathematical optimization proved that none-of them could be the best option for all possible optimization problems. However, there can be a better one for one specific task. Moreover, the firefly algorithm has a similar weakness as the marine predators and other such algorithm, which uses many iterations in a single step. These two things are summarized in [44]. Therefore, the hybrid evolutionary optimization methods were then introduced to overcome such disadvantages [45][46][47][48][49].
In this circumstance, improved meta-heuristic techniques and hybrid optimization techniques were then introduced and applied to the ORPD problem to find better solutions to the problem. In [50], a new modification was inserted to the sine-cosine algorithm and applied to the ORPD problem [45], which also introduced a hybrid modified imperialist competitive algorithm (MICA) and invasive weed optimization (IWO) for the same target. An enhancement was made to the firefly optimization algorithm to solve a multi-objective ORPD in [51]. More hybrid optimization techniques were proposed in [52][53][54][55][56][57][58][59]. In [52], the ORPD problem was solved using the hybrid GA and interior point method [53], which introduced a solution to the ORPD problem using hybrid artificial bee colony and differential evolution algorithm tested with the IEEE 14-bus, 30-bus, and 57-bus systems. Deferential Evolution itself needs a large number of agents to avoid premature settlement, which takes a long time. Moreover, the results were compared with DE itself. In [54], hybrid teaching learning and a double DE algorithm was proposed and applied to a single objective function, which is power losses minimization in the IEEE 14-, 30-, and 118-bus systems. Moreover, a comparison between the results were provided. In [55], modifications were made to the social spider optimization (SSO) method on how the algorithm generates new solutions. Then, this was used for solving the ORPD problem with the IEEE 30-and 118-bus systems. The weak point of the SSO itself was that it contains many control parameters, which increases the implementation complexity and execution time. Two single objective functions were studied. In [56], hybrid multi-swarm PSO was introduced to solve a multi-objective ORPD. The idea was that, each swarm was divided into sub-swarms, then a PSO was used as a search engine for the sub-swarms. Additionally, the DE method was used to enhance the best of the particles themselves. The studied systems were the IEEE 30-bus and practical 75-bus Indian systems. In [57], ORPD was solved using a hybrid Nelder-Mead simplex and a firefly optimization method (FA). The weak point of the FA was the trapping into local optima. Two single objectives were studied; power loss minimization and voltage deviation minimization. The research was applied to the IEEE 30-, and 118-bus systems. In [59], ORPD was solved by the proposed modified NSGA-II algorithm. The problem was a multi-objective function, which minimized power loss, besides maximizing the voltage stability. The systems tested in that work were the IEEE 30-and 118-bus systems. The single objective function optimization problem focused only on getting the best solution for one target at a time. On the other hand, in the multi-objective optimization problems, more than one objective was to be targeted simultaneously. Accordingly, introducing a more flexible algorithm with minimal adjustable paragliders was the motivated key for developing a new hybrid algorithm.
The marine predators algorithm (MPA) is a population-based technique [60] that was inspired mainly by the foraging strategies of the ocean predators and their interactions with the prey. The MPA was applied to solve several optimization problems in the field of solar energy. The MPA was applied to large scale photovoltaic (PV) array reconfiguration [61] and for PV modeling [62,63]. Furthermore, the MPA confirmed an efficient performance for forecasting COVID19 cases [64]. The simplicity of MPA in the implementation attracted the researchers to examine it with their applications, whereas the MPAs main disadvantage was that, the number of iterations was divided between the main phases of the algorithm. Therefore, it required huge number of iterations, especially with nonlinear optimization problems. For solving this point, an improved version of the MPA was proposed, based on distribution and particle swarm optimization. As a result, an improved marine predators algorithm and particle swarm optimization (IMPAPSO) method was introduced in this study. This was a salient and new contribution of this research work.
This study introduces an inspection to the newly developed MPA with some modifications and hybridization between the MPA and PSO, in order to optimize the decision variables of the ORPD. The proposed algorithm was examined with three systems of the IEEE 30 bus, IEEE 57 bus, and IEEE 118 bus, and compared with recent reported techniques in literatures. The effectiveness of the IMPAPSO was verified as an optimization tool to be utilized in solving the problem of ORPD for three bus systems of the IEEE 30 bus, IEEE 57 bus, and IEEE 118 bus. Moreover, the proposed algorithm was performed competitively, in comparison with the MPA itself and other optimization tools. Therefore, the research study focused on selecting the best possible decision variables to achieve the lowest possible active power loss in the grid. This paper also demonstrated how the voltage of the PQ buses could be optimized. In this analysis, the competence of the proposed optimization methodology was verified by different test systems. The objective functions studied in this research are single ones.
The remaining parts of the manuscript are organized as follows. Section 2 describes the mathematical formula of the proposed objective function. The overview representation of the MPA and IMPAPSO was documented in Section 3. Section 4 presents the simulation and results, and finally, Section 5 summarizes the main outcomes and the future work.

Problem Formulation
Generally speaking, the ORPD target was to optimize a non-linear function [63]. The constraints must be guided within their maximum and minimum boundaries, through the process of optimization [65]. It could be represented mathematically by Equations (1) and (2): Subject to: g(x,u) = 0, The symbol F is used to refer to the objective to be minimized. Meanwhile u and x are selected to refer to the preset control variables as well as the dependent variables.
There are three variables defined as control variables in this study. They are-(1) the voltage of the generator bus, (2) the tapping of the transformers, and finally, (3) the reactive power entering the studied electric power system from the VAR compensators. They are represented by Equations (3) and (4): Additionally, the dependent variables are mathematically represented in Equation (4): The ORPD objective function seeks the minimum active power losses. Additionally, the constraints must be considered during the simulation process. The optimal solution was found by varying the decision variables mentioned in Equation (3). There are three defined dependent variables in this problem, which are expressed as shown in Equation (4). These dependent variables were-(1) the active power generation at the slack bus and (2) the reactive output power at the voltage bus. Finally, the third one was the voltage at every load bus. The ORPD is illustrated in [34] as:

Objectives
This study focused on two single objective functions, which were-(1) active power losses minimization, and (2) voltage deviation minimization at the load buses. The first objective is expressed mathematically in Equation (5): where V and V are the voltages of the buses i and j, respectively, G is the mutual conductance between bus i and j , and θ is the voltage angle difference between buses i and j . The second aim was to reduce the voltage deviation at the load buses to a minimum value. This is mathematically represented by Equation (6): where V is the reference of the buses voltages, and N is the number of PQ buses.

The Equality Limitations
The equality constraints are represented in Equations (7) and (8): , iϵN pq (8) where is the susceptance between buses and ; and are the real and reactive generated powers at bus ; , represent the active as well as the reactive load of the bus , represents the reactive power of the shunt compensator at bus i , N represents the bus number, and stands for the PQ buses number.

The Inequality Constraints
The defined inequality constraints of the ORPD problem could be found in Equations (9)- (14): where is the transformer 'm' tapping variation, defines the number of PV buses in the electric power system that was under study, defines the number of on-load tap changing transformers. Moreover, the symbol defines the number of shunt VAR compensator in the system, is the branches number, is the real power at the slack bus, besides represent the minimum and maximum boundaries of the power at the slack bus. Additionally, defines the apparent power transmitted over a branch .

The Proposed IMPAPSO Approach
In this section, a novel version of the marine predators algorithm (MPA) is introduced to enhance the algorithm's main cores. The MPA is a novel meta-heuristic algorithm developed to emulate the foraging strategies of the ocean predators and their interactions with the prey. The predators determine the most convenient strategy for maximizing the rate of encounter with their victims according to the theory of the "survival of the fittest". The MPA could be performed via three main phases across the number of iterations. Each phase's control equation was varied based on the relative velocity between the predator and the prey. The details of the proposed approach of IMPAPSO was described as:

MPA Formulation
The MPA is mainly inspired by the ocean predators' drilling strategy. Depredators determine the most convenient strategy for maximizing the rate of encounter with their victims according to the theory of "survival of the fittest". Generally, a hazardous walking strategy, is the forging style of many predators.
Like many metaheuristic optimization techniques, the MPA optimization technique is considered to be a population-based optimization method. The initial distribution of the first population over the exploration field was consistent: where and are the minimum and maximum limits of the variables; rand stands for a random vector whose values were between 0 and 1.
According to the concept that the fittest ones survive, the top predators have the highest skills in foraging. Hence, the fittest population is chosen to be the most skilled predator when constructing the Elite matrix. Elements of the matrix supervise exploration and catching the victims, according to their locations.
where ⃗ is the vector of the most skilled predator. This was simulated n times to form the Elite matrix, where n is the population size while d is the problem dimension. Both predator and prey were considered to be population. At the end of every iteration, the Elite would be modified if the fittest predator was replaced by a better one. Prey is a matrix of the same dimension as Elite. The predators change their locations accordingly. The prey is represented as follows [61]: The process of the simulations by the MPA was strongly related to the Prey and Elite matrices.

MPA Phases
The MPA simulation consisted of three major phases. Each phase was defined by a specific range of velocity ratio-(1) the first phase was the one in which the velocity ratio was high, during this phase, the victim was considered to be faster than the predator; (2) in the second phase the velocity ratio was assumed to be in units and (3) in the third phase, the velocity ratio was assumed to be lower. The predator's velocity was higher than the victim. There is a preset period of iterations for each phase.
Phase 1: This was considered in the initial iterations of the simulation process in which the exploration matters. In this phase, the prey moved with Brownian distribution, while the predator did not move at all, accordingly the velocity ratio between them was high. This phase was mathematically represented as: While where ⃗ represents a random-numbers vector. The multiplication of ⃗ by prey mimics its motion. P = 0.5 is an arbitrary constant. R represents a random numbers vector between 0 and 1. The iterations were assumed to be divided into three groups. Phase one was applied to the first group. The step size was assumed to be high in this phase. Iter stands for the ongoing iteration and Max_iter stands for the maximum iteration. Phase 2: In the second phase, the exploration as well as the exploitation were considered. The population was divided into two parts. The first part was selected for the exploration and the second one was selected for the exploitation process. Prey is responsible for the exploitation process. On the other hand, the predator is responsible for the exploration process.
While _ < Iter < _ Equation (19) represents the rule applied to the first half of the population where ⃗ is a vector of random numbers. The multiplication of ⃗ and the Prey mimics the prey motion. The rule applied to the second half of the populations is mathematically represented in Equation (20): While the adaptive parameter is The step size for the predator is controlled by this adaptive parameter. The Multiplication of ⃗ by the Elite matrix mimics the motion of the predator. The prey also changes its position according to the motion of predators. Phase 3: This is the final stage of the simulation. It is the phase of high exploitation. The last stage is mathematically represented as follows: While Multiplication of ⃗ and Elite mimics the motion of the predator.

FADs' Effect
The environmental issues affect the behavior of the marine predators. Fish aggregating devices (FADs) are an example of the environmental issues. The FADs can represent local optima in the exploration field. Assumption of longer jumps through the optimization process could avoid trapping in the local optima. Equation (22) represents the FADs effect in a mathematical form: where = 0.2 stands for how probable FADs affect the simulation progress, ⃗ is a binary vector. It changes the array to zero in case that the array is less than 0.2 and changes to one in case the array is higher than 0.2. The symbol 'r' defines a random number between [0,1]. ⃗ and ⃗ is the vector of the minimum and maximum limits of the dimensions. 1 and 2 are random indices in the prey matrix.

Marine Predator Memory
Marine predators are known by their good memory. They remember the position in which they succeeded in foraging. The ability to remember the positions is represented by memory saving. First, the prey is updated as well as the FADs. The objective function is then measured and the Elite is revised. The flowchart of the MPA is shown in Figure 1.

Particle Swarm Optimizer (PSO)
The PSO is a very commonly utilized swarm optimization technique, which was inspired by the behavior patterns of swarms, such as their communication strategy while exchanging information, chasing victims, traveling, and gathering, as introduced in [66]. The particles are the population seeking a solution by using experience to update themselves [66]. These particles are initially randomly diversified through the exploration field having positions ( ) and velocities ( ) in the dimension. After that through the simulation process, they update positions according to the individual and global best locations [67]. It is mathematically represented as follows: where and are the location and the velocity of the particle in the dimension, respectively. The symbol refers to the ongoing iteration, represents the inertia weight. The symbols and denote acceleration coefficients. ( ) is the best location of the agent, while the ( ) is the global best agent. The values and are random numbers within 0 and 1. The fitness function is evaluated iteratively then; the best and global best candidates are consequently updated. The aforementioned sequence is repeated until the stopping criteria are reached. The best solution is then obtained, as shown in Figure 2.

The Proposed IMPAPSO Structure
The profile that the agents follow while modifying their positions is one of the essential factors that reflect the algorithm response. Therefore, in this work two strategies are applied to enhance the MPA behavior, where the random walk of the predictor is boosted by the high-tailed Weibull distribution and its exploitation phase is strengthened by the PSO operator, as described in the following points and described in Algorithm 1.
In the first modification, the random walk of the agents in MPA that is controlled by Brownian movement to modify the predator movement is replaced by the high-tailed Weibull distribution ( ) [68], to boost the exploration stage, based on the following formula: where and are the scale and shape parameters. In the current work, the tuned values of and are 2 and 1, respectively. In the second modification, the PSO operator is utilized to enhance the exploitation phase of the basic MPA as a trial, to enhance the convergence behavior of the algorithm. In phase 3 of the MPA, the agents can follow MPA or PSO operators, based on a probability value, as explained below: where is a generated value on basis of the fitness function = ∑ , is a random vector drawn from uniform distribution [0,1]. If the value of prob < rand , then the operators of MPA of Equation (20) will be implemented; otherwise the operator of PSO of Equations (23) and (24) is executed.

Simulation Results
In this study, the three systems are tested, these are IEEE 30-bus, 57-bus, and 118-bus standard test networks. The main reason for choosing these systems is that, they are the most popular systems used for the ORPD problem. Moreover, the MPA quality can be tested with increasing dimension systems. The standard IEEE 30-bus system contains 5 voltage buses, 21 loads, 41 branches, 4 changers, as well as 3 shunt capacitors [68]. Regarding the IEEE 57-bus test system, it contains 6 voltage buses, 41 loads, 80 branches, 17 tap changers, and 3 shunt VAR compensators. Additionally, in the IEEE 118bus system, there are 54 voltage buses, 99 loads, 186 branches, 9 tap changers, and 14 shunt VAR compensators. The proposed IMPAPSO is inspected when applied to the aforementioned test systems for solving the ORPD Problem. Each system is studied from two points of view. These are named-'case 1' and 'case 2'. Case 1, meanwhile, decreases the active power losses, and case 2 aims to reduce the deviation in voltage. To put constraints into their limits, penalty functions are inserted with constants of high value to the objective function.

Results of the IEEE 30-Bus System
The system data are stated in [35,67]. In this network, nineteen control variables are considered. These can be stated as follows-the voltage of the generator at buses 1, 2, 5, 8, 11, and 13, transformer tap settings between the buses 6-9, 6-10, 4-12, and 28-27, in addition to the reactive power of the VAR compensators at buses 10, 12, 15, 17, 20, 21, 23, 24, and 29. The maximum and minimum values of the voltages are 1.06 and 0.94 p.u., respectively. The maximum and minimum limits of the transformer taps are 1.1 and 0.9, respectively. The minimum limit of the VAR compensators is 0 while the maximum limits are-6, 20, 6, 7, 6, 10, 6, 10, and 6 MVARs respectively. In Table 1, the results of the objective targeted in case 1 are presented in comparison with the other methods' results. Moreover, Table 2 represents the results of the tuned decision parameters in case one, when handled by the proposed IMPAPSO. The results of these variables are compared with the results of other methods such as the HHO [66] as well as the PSO.
In the second case, the IMPAPSO is compared with other optimization methods, it reaches the minimum voltage deviation. Numerical results confirming this superiority are presented in Tables 3  and 4. To affirm that the IMPAPSO is dependable on solving the ORPD problem, many trials are done for the two objectives under study. The convergence curves of the objective functions handled by the proposed IMPAPSO versus PSO and MPA for both the first and the second cases during 600 iterations are shown in Figures 3 and 4. It could be observed that, the solution convergence of the objective functions was efficient and competitive. The values of the two objectives optimized by the IMPAPSO were the best among the methods compared.   The results of power losses in the 30 bus system obtained by different optimization methods are presented in bar chart, as shown in Figure 5. A comparison of these results showed that the value of power losses obtained by the proposed IMPAPSO optimization algorithm was the best among the results obtained by the other algorithms, especially the BFOA.    The value of the voltage at every bus in the 30-bus system is shown in Figure 6. A comparison was carried out between the bus voltage, before and after performing the ORPD problem, using the proposed IMPAPSO optimization algorithm. It was clear that the voltage profiles improved and came closer to 1 p.u., after using the IMPAPSO in the simulation process.
In case 2, the proposed IMPAPSO was compared with other optimization methods, it reached a minimum voltage deviation. Numerical results confirming this superiority are presented in Tables 7  and 8. To affirm that the IMPAPSO depended on solving the ORPD problem, many trials were done for the two objectives under study. The convergence curves of the objective functions handled by the proposed IMPAPSO, MPA, and PSO for both the first and the second cases during 600 iterations are shown in Figures 7 and 8. The convergence of the objective functions solutions was competitive. The values of the two objectives optimized by the IMPAPSO were the best, compared to the MPA and PSO methods.   The results of power losses in the 57 bus system obtained by different optimization methods are presented in the bar chart, as shown in Figure 9. A comparison of these results showed that the value of power losses obtained by the proposed IMPAPSO optimization algorithm was the best among the results obtained by other algorithms especially, the PSO.  Similarly, as provided in the 30 bus system, the value of the voltage at every bus in the 57 bus system is shown in Figure 10. The bus voltages were compared, before and after performing the ORPD problem, using the proposed IMPAPSO optimization algorithm. It could be noticed that the voltage values improved and came closer to 1 p.u., after using the IMPAPSO in the simulation process. Figure 10. Bus voltages of 57 bus system before and after optimal reactive power dispatch (ORPD).

Results of the IEEE 118-Bus System
The system data are stated in [35], the number of control variables considered in this system was seventy-seven. The maximum and minimum values of the voltages were 1.06 and 0.94 p.u., respectively. The maximum and minimum limits of the transformer taps were 1.1 and 0.9, respectively. The limits of the VAR compensators varied between −20 and 20 MVARs. In Table 9, the results of the objective function in case 1 are presented in comparison with the other methods' results. Moreover, Table 10 represents the results of the adjusted control variables in case 1 when handled by the IMPAPSO. The findings are compared with the results of other approaches, including FA, ABC, and PSO. In case 2, the IMPAPSO was compared with other optimization methods, where it reached a minimum voltage deviation. Numerical results confirming this superiority are presented in Tables 11 and 12. To affirm that the IMPAPSO was dependent on solving the ORPD problem, many trials were done for the two objectives under study. The convergence curves in case 1 and case 2 of the objective functions handled by the proposed IMPAPSO, MPA, and PSO are shown in Figures 11 and  12. It could be observed that, the convergence of the solution of the objective functions was efficient and competitive. The values of the two objectives optimized by the IMPAPSO were the best among the methods compared.  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35  36  37   The results of power losses in the 118 bus system obtained by different optimization methods are presented in the bar chart, as shown Figure 13. A comparison of these results showed that the value of power losses obtained by the IMPAPSO algorithm was the best among the results obtained by other algorithms, especially the ABC and the FA.    The voltage at every bus in the 118-bus system is shown in Figure 14. The bus voltages were compared before and after performing the ORPD problem, using the IMPAPSO algorithm. Notably, the voltage values were improved after using the IMPAPSO in the simulation.

Conclusions
The paper presented a new application of the IMPAPSO to solve the ORPD problem. The objective functions were also mathematically formulated and introduced. The optimization technique was inspected with 30, 57, and 118 bus test systems. The review observed that, it was mostly studied to minimize real power losses in power systems, especially in transmission networks for the ORPD problem. Therefore, in this article, two case studies were evaluated where the ORPD problem was solved by IMPAPSO, MPA, and PSO. Two points of view optimized the objective functions of the problem. First, the active power losses must be minimized. In addition, the second goal was to minimize the deviation in voltage. The performance of IMPAPSO was compared with other well-known methods incorporated optimizers. The tool in this paper, as shown by the simulation results, was able to find improved solutions to the ORPD problem. The advantages and superiority were investigated on various test systems. The employment of the proposed method resulted in an improvement in power loss minimization up to 96% in the 30-bus system, 10% in the 57-bus system, and 9% in the 118-bus system. Additionally, the improvement in voltage deviation minimization due to the usage of IMPAPSO was up to 18% in the 30-bus system, 33% in the 57-bus system, and 12% in the 118-bus system. Furthermore, the results of the simulation confirmed that, the IMPAPSO could be applied to additional electricity grid problems. Last, the results obtained by other optimization algorithms were statistically presented to compare the performance of the discussed algorithms. It is expected that the proposed IMPAPSO would be a magnet for the researchers who are interested in optimization of different power engineering problems.