Allocation of synchronized phasor measurement units for power grid observability using advanced binary accelerated particle swarm optimization approach

Large-scale power grid observability is still a challenge because of deteriorating infrastructure and the incorporation of renewable energy sources. A smart grid that makes use of cutting-edge technology, such as a phasor measurement unit (PMU), is an excellent option for monitoring and bringing networks up to speed with the latest information. Latterly, the considerable investment required for the deployment locations has slowed down the adoption of PMU. Therefore, because PMUs are expensive, it is necessary to deploy them in the best possible places on large-scale power grids. The most significant share of optimal PMU placement problems (OPPP) is defined as 0–1 knapsack problems. Considering this, the development of an effective optimization technique that can handle difficulties has emerged as an appealing topic in recent years. In this paper, a meta-heuristic algorithm based on the binary particle swarm algorithm (BPSO), a binary accelerated particle swarm optimization (BAPSO), is offered for solving OPPP. Since earlier research has shown that BPSO is likely to stick to local optima, the majority of them evaluated their suggested technique using small-scale test systems. The technique that has been suggested searches for the optimal solution by employing two topologies—one global and one local—that are analogous to BPSO. This work determines the optimal PMU position for a large network in a reasonable amount of time by fine-tuning the acceleration factor. Additionally, in order to employ fewer PMUs, an integration strategy was put into place for the radial buses. The OPPP solutions are provided by the suggested method within a reasonable period with prior solutions published in reliable publications, according to computational findings.


Introduction
The power networks of today are being run under difficult conditions in order to supply the rapidly growing demand for electrical resources and to keep the commercial activity going in the midst of a very dynamic, deregulated market.Therefore, power grid monitoring, preservation, and control become increasingly important for enhanced systems operation, maintenance, planning, and energy trading.As a result, PMU has evolved as a valuable piece of equipment for measuring phasors of voltage and current, which are synchronized with signals collected using GPS technology.They can enhance operations such as bad data detection, corrective action schemes, state estimation, stability control, and disturbance monitoring.When it comes to installing PMU on the electrical grid, one of the most important issues that have to be taken into account is the expense of doing so.It is, therefore, of the utmost importance to determine the optimum position of the PMU where reliability is maintained while minimizing the costs involved.Recently, many methods to analyze the OPPP employing various sets of optimization algorithms have been presented.Deterministic and stochastic algorithms are two categories that describe optimization strategies that can be used to solve the PMU placement problem.
The integer linear programming (ILP)-based formulation to evaluate the OPPP was initially suggested in [1].This formulation, in which linear constraints are established based on a binary bus-to-bus connection matrix, evaluates the observability of power networks considerably more simply and straightforwardly.An ILP method was suggested in [2]; this method took traditional measurement as well as zero injection bus (ZIB) into consideration.The use of a permutation matrix is included in the suggested method, which helps reduce the nonlinear limitations.There is also an explanation of the idea of partial observability in [3].In addition, a malfunctioning of PMU was incorporated into the strategy that was presented.The bus observability index (BOI) and the system observability redundancy index (SORI) were described so that the optimal PMU employment set could be obtained.
To address the effective solutions specified for the OPPP while taking into consideration the impact of ZIB, such as line and PMU outages, a mixed ILP (MILP) is presented in [4] and [5].A technique that is based on integer programming and genetic algorithms (GA) was developed by [6] to install PMU in order to obtain full monitorability of the power network.A combination of GA with a simulated annealing strategy was offered by Kerdchuen and Ongsakul in [7] as a way of obtaining a solution for the OPPP.In [8], researchers investigated a unique cellular GA-based approach for OPPP that takes into consideration the availability of channel capacity as well as single-line loss.Ahmadi et al. [9] recommended using conventional BPSO to decide on the OPPP with and without ZIBs.The measurement of redundancy is presented as a method for ranking the solutions.
In the research conducted by Chakrabarti et al. [10], an enhanced particle swarm optimization (EPSO) for power grids, as described by Valle et al. [11], was applied to the OPPP.Further velocity update rules are implemented by EPSO if the particles cannot identify a viable solution.Similar to the study conducted by Chakrabarti et al. [10], the authors of [12] suggested a novel velocity update equation to locate the OPPP using BPSO.In addition to the velocity update equation, the authors created additional observability techniques for ZIB, a PMU, and a line failure.In reference [13], the authors introduced the exponential BPSO as a novel way of controlling the inertia mass of BPSO.The authors assert that it improves the searchability of the method.Wang et al. [14] presented a hybrid technique for the OPPP that combines simulated annealing and BPSO.In order to place PMUs in power distribution systems optimally, a tri-objective strategy has been presented in [15].Its goals are to reduce the number of PMU channels, state estimation uncertainty, and sensitivity to line parameter tolerances.Observability propagation depth and probabilistic observability are taken into account in [16] to improve formulation for the best placement of PMUs in power grids.
A two-stage approach to optimize the placement of PMUs was proposed in [17,18] to achieve complete system visibility while minimizing cost, taking into account objectives such as cost minimization, redundancy, and efficiency maximization, as well as constraints such as zero injection buses, single PMU failure, single-line outage, and flow measurements.Article [19] addresses the issue of incomplete observability under single PMU loss (N − 1) contingencies and proposes an enhanced two-archive algorithm and a fuzzy decision-making method for PMU placement optimization.In ref. [20], author has presented a BPSO technique for the optimal allocation of PMUs in connected power networks, demonstrating its effectiveness and superiority compared to other methods through testing on various test systems.In addition, a technique for integration is not used for the radial buses; instead, the approach entails taking into consideration as well as ignoring ZIBs.
In this article, author has proposed a meta-heuristic algorithm, based on BPSO, a BAPSO, to solve the OPPP in large-scale power grids, aiming to find the best locations for deploying expensive PMUs and achieve grid observability while considering the challenges posed by deteriorating infrastructure and cost constraints.The algorithm combines global and local search topologies and fine-tunes the acceleration factor to efficiently determine the optimal PMU positions, and it also incorporates an integration strategy for radial buses to reduce the number of PMUs required, providing solutions within a reasonable time frame compared to the previous research.As BAPSO is a meta-heuristic algorithm, it is expected to generate multiple PMU placement sets and to determine the quality of each set with the same number of PMUs, the one with the highest SORI value is chosen as the optimal result.The PMU placement set with higher measurement redundancy is considered better than the one with lower measurement redundancy.BAPSO is proposed to determine the minimum number and optimal locations of PMUs for complete monitoring of the power grid, taking into account factors such as normal operation and zero injection measurements.

Method used for the optimal PMU placement problem
In general, the primary goal of the OPPP is to obtain the fewest number of PMUs necessary, along with the location of those PMUs, to ensure full observability of the power grid.As a result, the following is the model for the generalized objective function that is used for the identification of the OPPP in this work [21]: (1) where n is the number of buses, c i is the vector of PMU price coefficients, Y is the binary design variable vector having components y i which decide the feasibility of PMUs on ith bus, and H and B PMU are interpreted as the transformation matrix that may be modified according to the contingency cases.A PMU = [A i,k ] n×n is the binary connectivity matrix that explains the bus-to-bus connection whose entries are shown in Eq. (5).
Y provides the decision for the placement of PMU as given in Eq. (6).
] n×1 is the column vector that signifies the redundancy, that is, essential for the specific case.

Radial BUS
It is noted that installing PMU on a bus that is linked to more than one neighboring bus would have greater coverage of the connected power grid relative to the bus that has very few adjacent buses, in particular, the radial bus network [22].Hence, if the PMUequipped bus is radial, the PMU can only monitor two buses-the radial bus and its neighbor.Radial buses are excluded from prospective OPPP solutions since their PMU setup will measure the voltage phasors on that bus and one associated bus.

Modeling of ZIB
The consideration of ZIB may benefit in further reducing the PMU numbers necessary to achieve maximal observability of the power grid.Several methods for coping with ZIBs have been suggested in the previous research.The bus integration approach is one of the strategies that have been established to cope with the characteristics of the ZIB [22].The bus integration strategy requires an integration process between the ZIB and one of the neighboring buses.As a consequence of this, during the process of integration, the limits placed on both buses may be combined into a single constraint.As a result, the number of constraints that need to be satisfied to guarantee that the installed PMUs will observe each bus will be reduced.It is believed that if all observable buses except for the unobservable one are interconnected to the ZIB, then the unobservable bus can be construed as being observable.Because of this, the integration of the bus shows that if it is measurable, the bus that was picked to be integrated will also be observable.
The 14-bus system of IEEE is taken into consideration to comprehend the bus integrating strategy.This system is illustrated in Fig. 1, and it is important to note that bus "7" is a ZIB, and it is coupled with bus "4, 8, " and "9." To identify a candidate bus to integrate with the ZIB, the following process may be utilized: (i) randomly integrate the ZIB with one of the buses that are near it.In this example, bus 7 is integrated with one of its neighbors.As an illustration, bus 7 and bus 9 are combined into one, (ii) integrate the (5) ZIB with one of the surrounding buses that have the fewest number of buses attached to it.Using this technique, bus 7 is integrated with bus 8, which only has one bus connected to it; therefore, the total number of buses linked to bus 8 is reduced to one, and (iii) integrate the ZIB with its adjacent buses that have a higher number of buses linked to it-in contrast with the plan that was shown previously, bus number 7 is interconnected with bus number 4, which has a total of five buses that are connected to it.Bus number 9 only has a total of four buses that are connected to it.
When dealing with the presence of ZIB, a bus integrating approach may be used to establish the bare minimum number of PMUs that are required; despite this, there are a few drawbacks that need to be brought to the attention: (i) If a PMU is necessary to be installed on an integrated bus, this might imply that the PMU has to be installed on the ZIB, or on the bus chosen to be integrated with the ZIB, or on both buses.However, it could also mean that the PMU needs to be installed on both buses.Because of this circumstance, a further monitoring test needs to be conducted in order to determine which of the two buses should have the PMU placed on it, and (ii) every time an integration procedure has been carried out, the topology of the system has been modified.When referring to a power grid on a massive scale, this may cause the topology to become more complicated.

OPPP rules without ZIB
Rule 1 A PMU installed at a specific bus has the ability to compute not only the voltage phasors of that bus but also the current phasors of all of the lines that are related to it.In Fig. 2, bus {1} is PMU-equipped bus.Here, V 1 , I 12 , I 13 , and I 41 can be unswervingly meas- ured by the employed PMU.Rule 2 It is possible to calculate the voltage at the other end if the voltage at one end and the line currents of that end are known.Taking into consideration (Fig. 3) and assume that the values of the line current I 12 , I 13 , and I 41 are known, then Ohm's law can be used to compute the voltages at the buses {1} , {3} , and {4} .The values of V 2 and V 3 are the results of V 1 subtracting the potential drop induced by current flowing over the line.Therefore, the values of V 2 , V 3 , and V 4 are solved as follows: (8) Rule 3 If the voltages at both ends of the buses are known, then Ohm's law may be utilized to calculate the line currents that flow between the buses.Given (Fig. 4) that the values of V 1 and V 2 are well known, the line current I 12 may be determined by using Ohm's law, which is presented in the following form:

OPPP rules with ZIB
A bus is known as a ZIB when neither the load nor the generator is connected.As a result, the summation of line currents used at a ZIB is zero.If ZIB, which includes its neighbors, has N z members, then monitoring N z − 1 buses is enough to turn an unobservable bus into an observable bus.Because of this, while considering ZIB, the number of buses that need to be observed drops by one for each ZIB that is present in the power grid.This, in turn, reduces the minimum number of PMUs that are required for total observability.The following PMU observability criteria are implemented to analyze the topological observability using ZIB: Rule 4 If there is one bus that is not observable that is adjacent to a ZIB that can be observed, then the bus that cannot be observed can be deemed to be observable.Take, for instance, if the values of V 1 , V 2 , and V 3 are known, then V 4 may be determined with the use of the KCL at bus {2} which is a ZIB.Refer to Fig. 5, where bus {2} is a ZIB that is observable.Assuming for the moment that the values of V 1 , V 2 , and V 3 are identified, then the value of line currents I 12 and I 23 can be determined by using rule 3 as mentioned above.So, by using KCL at bus {2}, the value of I 12 is I 12 = I 23 + I 24 .For that reason, the value of I 24 and V 4 can be obtained as follows: (14) Fig. 4 Modeling PMU placement rule 3 Rule 5 If the observable buses are linked to ZIB which is unobservable, then the ZIB can be considered as observable.Consider Fig. 6, where bus {2} unobservable ZIB which is connected by all the observable buses such as buses {1}, {3}, and {4}, then the voltage of bus {2} can be obtained as follows:

Algorithmic perspective of BAPSO in OPPP
Refer to [23], Eqs. ( 22) and ( 23) are used to update the velocity of a particle i at each iteration m in the original PSO.
(17 where w (m) is the inertia weight at the m iteration [24].The inertia weight in PSO con- trols the dynamics of flying among particles, and a higher value of this weight leads to global exploration, while a lower value promotes local search.If the inertia weight is set too high, the algorithm may focus too much on exploring new areas and neglect local search, making it challenging to find the exact optimal point.In order to accelerate convergence to the true optimum by balancing global and local exploration, a linearly decreasing inertia weight has been employed: The values of inertia weight are of w max = 0.9 and w min = 0.4 , and M max is the max- imum number of iterations used in PSO [24,25].By introducing a virtual mass to stabilize the motion of the particles, the algorithm is anticipated to have a faster convergence rate.A velocity threshold is introduced [26]: where v ij (m + 1) is the velocity component of the ith particle along the jth direction at the (m + 1)th iteration of the algorithm, and v max j is the maximum absolute value of velocity allowed along the same jth direction in the parameter space.The adaptation of the inertia weight allows the swarm to achieve convergence with greater accuracy and efficiency as compared to the original PSO.r 1 , r 2 are random vectors from the uniform distribution in the range [0, 1] to maintain the swarm diversity.The acceleration constants are c 1 = c 2 = 2 called cognitive parameters, so that c 1 r 1 and c 2 r 2 ensure that the particles would overfly the target about half the time.
The present study presents a meta-heuristic optimization algorithm named BAPSO, which builds upon the BPSO algorithm by incorporating global and local topologies.The BPSO algorithm is known to face the issue of premature convergence and tends to get stuck in local minima.However, the newly introduced mutation strategies in BAPSO can effectively prevent agents from being quickly trapped in local optima, especially when dealing with complex combinatorial problems.BAPSO has the capability to explore the entire solution space for a global search and conduct a local search, leading to the identification of global minima [27].The evolutionary equation of BAPSO is as follows: where a is the acceleration factor.Compared to the conventional PSO, the evolution equation of BAPSO involves an additional parameter "a," while it includes one more parameter "w" than the PSO with the contraction factor.Despite this, BAPSO has demonstrated impressive results in solving complex OPPP for large-scale inter-power grids within a reasonable time.The proposed BAPSO in this study shares a similar structure with PSOCF.The equation of PSOCF is given as follows [28]: (23) .9(3.9 − 4) does not exist.The PSOCF could not be used, but BAPSO in this work is used for solutions of the OPPP of large-scale inter-power grids and obtained satisfying results.For The acceleration factor a is outside the space of the contraction factor , and the optimization of BAPSO in this article is performed as usual.Although BAPSO is similar to PSOCF, PSOCF did not have good adaptability as BAPSO [27].The pseudo-code of the BAPSO is as follows: For loop over all n particles and all j dimensions Generate new velocity v PSO was originally intended to handle unconstrained optimization, but it has the potential to solve constrained problems with modifications.To locate the global minimum while accounting for constraints, BAPSO employs a constraint-handling approach that updates both a particle's best position and the swarm's global best position.To steer the search toward the feasible area, a feasibility term is included, which determines the extent of the overall constraint violation.The choice of the global best (Gbest) topology in BAPSO depends on the dimension of the search space.In order to enable BAPSO to work with binary problems, the initial Gbest is represented as a binary column vector [25].The population size is selected according to the network size [29].The initial inertia parameter could be selected as w = (0.9 − 0.7 × rand).
The objective is evaluated with a number of moving particles at each iteration.As observed, the BAPSO starts with the iteration to find the global minimum point, whereas the velocity tends to go into v max or −v max .The value of v max is carefully selected [25].When the size is insufficient, the algorithm can get stuck in a local minimum or have to perform more iterations to arrive at the correct solution.The particle is positioned within the binary search space [27], and its current velocity and position impact its future position.The BAPSO is capable of conducting both global and local searches (26) of the solution space without being confined to local minimum points.To attain better convergence, an inertia weight is utilized to maintain a balance between global and local searches.

Particles
For the OPP problem, every particle has a promising solution.The objective of this work is to determine the optimal minimum number and strategic locations of PMUs to maximize the observability of the power grid.As a result, the configuration of each particle is designed to indicate the availability of PMUs on a particular bus.When determining the OPPP for a 7-bus system (as shown in Fig. 7), the construction of each particle is depicted in Fig. 8, which can be found below.Each dimension of the power grid is linked to a specific bus, and each particle is developed according to these dimensions.A value of {1} at bus {2} indicates that a PMU is installed at that bus, while a value of {0} denotes that there is no PMU installed at bus {2}.

Redundancy measurement
In order to determine the most effective sets of PMU placements, the BOI and SORI redundancy measurement concepts, as described in reference [30], are utilized.BOI refers to the number of times a particular PMU observes a bus, while SORI is the sum of all BOI values.The solution sets that have the least number of PMUs and the greatest sum of BOI, represented by SORI, are considered to be the most optimal.The BOI is the performance metric, which can be calculated using Eq. ( 27), while Eq. ( 28) shows how to calculate SORI.

Fitness function
The BAPSO involves particles that carry potential solutions to the OPPP, and in order to determine the best solution, a fitness function is used to evaluate each solution during the investigation.The fitness function must meet three important criteria: ensuring power grid observability, determining the minimum number of PMUs needed for full observability, and measuring redundancy.Following these guidelines, the fitness function for identifying the desired target can be expressed as shown in [22].
where w 1 (= −2) , w 2 (= 1) , and C (= 0.01) are the weight parameters, N obs is the total number of a bus which is observable, N PMU is the number of PMUs equipped bus, and R 1 is the redundancy measurement.The fitness function described in Eq. ( 29) is com- prised of three components: (i) the count of observable buses, (ii) the count of PMUs, and (iii) the redundancy measurement.It is important to highlight that the first component determines the number of buses that can be monitored through the placement of installed PMUs.The value of N obs can be given as follows: Additionally, the second component determines the quantity of PMUs, which can be interpreted as follows: Moreover, regarding the third component, the value of redundancy measurement is established by:

Results and discussion
The OPPP is solved using a modified particle swarm optimization approach in this study.The traditional BPSO method is limited by premature convergence and is prone to get stuck in local minima.The proposed BAPSO method, on the other hand, can carry out both global and local searches to locate global minima.It effectively prevents agents from quickly becoming trapped in local optima, which is particularly useful in addressing complex combinatorial problems.
In order to implement the proposed method effectively, it is necessary to determine the appropriate parameter values such as population size.To this end, various trial runs have been conducted on all the test systems studied for solving the OPPP, and the optimal results are presented here.The population size is four times the (28 number of buses, which is sufficient for solving the OPPP in the present work.The maximum number of iterations has been set to 250 for smaller systems such as IEEE 14-bus, IEEE 30-bus, New England 39-bus, and IEEE 57-bus systems, while it is 1000 for larger systems such as IEEE 118-bus, IEEE 300-bus, and NRPG 246-bus systems.MATLAB R2013a software was used to conduct the simulations, and the computer used had an Intel Core i3-5005U (2.0 GHz, 3 MB L3 Cache) processor and 8 GB of RAM.The number and location of radial and ZIBs are shown in Appendix "Number and location of a radial and zero injection buses, " while the connection of ZIBs is displayed in Appendix "Connection of ZIBs." The parameter values used in the proposed method for the PMU placement problem are listed in Table 1.
The parameter values used to unravel the OPPP were carefully selected through extensive testing to ensure feasible performance.The proposed BAPSO algorithm was found to converge faster than the standard BPSO algorithm for every bus system.The results obtained were satisfactory, and the proposed method achieved adequate computational time, which was only slightly longer than the standard BPSO algorithm.Interestingly, the computational time was found to be superior to that of existing studies.Additionally, the consideration of ZIBs and radial buses from the OPPP minimized the number of PMUs necessary for the entire power network observability.
According to Table 2, it is possible to ensure the observability of the power grid under normal operation for the standard 14-bus system by placing PMUs in the optimum locations and with the minimum number required.After considering the number of trials and with redundancy measurement, it is found that buses {2, 6, 7, and 9} are the best promising set for the OPPP.Here, solution set 2 has a maximum number of SORI, that is, 19.The entries of BOI signify that how many numbers of times the PMU-equipped bus observes each bus, and it is also explained in Sect."Redundancy measurement." Table 3 provides the details on the minimum number and optimum locations of PMUs required to achieve full observability of the power grid during normal Inertia parameter ( w max and w min ) 0.9-0.4  6, 9, 10, 13, 16, 17, 19, 20, 22, 23, 25, and 29} is the most promising optimal locations of PMU as depicted in Table 4 for NE 39-bus system to make power grid utterly observable during the normal operating condition.
The solution to the OPPP may be found in Table 6, and it is presented here for 118bus, 246-bus, and 300-bus systems accordingly.
Table 7 indicates that the number of PMUs required for achieving maximum observability increases with the expansion of the power grid.As the size of the network increases, computational time also increases.Table 8 displays the optimal position and minimum number of PMUs needed, taking into account ZIBs.It is worth mentioning  that the inclusion of ZIBs in simulations reduces the number of PMUs required for observing all buses.For instance, in a 14-bus system under normal operating conditions, four PMUs are required for maximum network observability, but with the consideration of ZIB, only three PMUs are needed.Tables 9 and 10 compare the results obtained from the BAPSO technique with those obtained from Guo [31], Chakrabarti et al. [32], Milosevic et al. [33], Manousakis et al. [34], and Sodhi et al. [35] for IEEE 14-, 30-, 57-, 118-, NE 39-, NRPG 246-, and 300-bus systems, with and without ZIBs, respectively.In this study, the BAPSO technique was employed to determine the optimal number of PMUs and their positions while maximizing redundancy measurement, ensuring full observability of the power grid.The BAPSO approach was applied to IEEE networks, and the results were compared with those obtained using various programming methods proposed in the previous literature.The comparative analysis demonstrates that the BAPSO technique provides alternative methods where the objective function takes a minimal value in full agreement with the ones defined by the current programming techniques for each case study.
Table 11 presents a comparison between the computational time of the proposed method and the results obtained from using the BPSO algorithm in recent studies.The study found that an increase in the number of buses resulted in a longer computational time.However, the proposed approach significantly outperformed the previous studies in terms of computational time.This demonstrates that the proposed    125,128,132,134,140,141,142,153,157,158,160,163,168,173,181,185,187,190,191,194,199,201,202,203  approach not only yields high-quality solutions but also operates at a faster computational pace.

Conclusions and scopes of future work
The purpose of this paper is to introduce a novel BAPSO algorithm that incorporates global and local topologies to solve OPPP, to enhance the learning and convergence procedure of classifiers.The proposed algorithms have numerous benefits, such as simplicity, ease of implementation, and the lack of need for algorithm-specific parameters.Instead, they require only common controlling parameters, such as the number of generations, population size, and tuning of the acceleration coefficient.The efficacy of the BAPSO algorithm in achieving an OPPP solution is demonstrated using IEEE bus systems.In binary PSO, the population size is a crucial factor in achieving optimal execution time and solution consistency.However, increasing the population size also increases the total execution time.The study finds that the algorithm's average execution time and performance are directly proportional to the size of the population and the maximum number of iterations.In a large-scale network, conventional BPSO can generate a set of optimum solutions, but it is not feasible within a reasonable timeframe.Conversely, the proposed BAPSO approach offers a fast OPPP solution  tations into the PMU allocation procedure.In order to obtain optimal PMU placement, take into account the effects of network structure, such as the presence of radial or meshed networks, and design a strategy that integrates topological considerations.(d) Investigate the best location for PMUs when using them for wide-area monitoring applications, taking into account local or global power grids.Create a framework that considers geographic and connectivity factors in order to improve situational awareness and system stability in massive power systems.(e) Cybersecurity Considerations: Examine the possible hazards and vulnerabilities related to the installation of PMUs in the electrical grid.To lessen the danger of cyberattacks and unauthorized access to vital power system infrastructure, look into ways to protect the security and integrity of PMU data and suggest solutions for safe PMU installation.(f ) Cost-Effectiveness Study: Conduct a thorough cost-effectiveness study to assess the financial advantages of the suggested PMU allocation strategy.Think about things such as the price of PMUs, installation, communication setup, and upkeep.Create optimization models that strive to achieve the required level of observability while minimizing the total cost.(g) Real-Time Implementation: Examine the viability and practicality of putting the suggested PMU allocation technique into use during the real-time operation of the power system.Think about the computational effectiveness, the communication needs, and the SCADA system integration.Create methods for real-time PMU allocation changes and ongoing power system monitoring.(h) Application to Renewable Energy Integration: The suggested PMU allocation methodology should be expanded to accommodate the unique difficulties involved in integrating renewable energy sources into the power grid, such as solar and wind.Develop methods for the best PMU deployment in grids with a high concentration of renewable energy sources by looking at the effects of distributed generation and intermittent power generation on observability needs.
(i) These suggestions can act as a springboard for more study, enabling the development and improvement of the suggested PMU allocation strategy and eventually advancing the observability and stability of the power grid.Velocity component of the ith particle along the jth direction at the (m + 1)th iteration of the algorithm v (m) Velocity of particle i at iteration m w 1 Weight parameter for the number of bus observed w 2 Weight parameter for the number of PMUs C Weight parameter value for the measurement redundancy

Algorithm 1 :
Pseudo-code for BAPSO Objective function f ( − → y ), − → y = y 1 , y 2 , ..., y n T Initialize locations y i and velocity v i of n particles Find Gbest from min {f (y 1 ), ...., f (y n )} at ( m = 0) while (criterion) m = m + 1(pseudo time or iteration counter) Assess the objective functions at new positions y m+1 i Determine the present optimum for every particle Pbest end ε Find the current global best Gbest end while output the final results Pbest and Gbest

Table 1
Configuration settings for the optimization method

Table 2
Optimum locations of PMU for 14-bus under normal operations

Table 3
Optimum locations of PMU for 30-bus under normal operationsBold highlighted values are results that are considered the best in terms of time efficiency.They represent optimal outcomes achieved swiftly and efficiently within a given context

Table 4
Optimum locations of PMU for 39-bus under normal operationsBold highlighted values are results that are considered the best in terms of time efficiency.They represent optimal outcomes achieved swiftly and efficiently within a given context

Table 5
Optimum locations of PMU for 57-bus under normal operations

Table 6
Optimum locations of PMU for 118-bus, 246-bus, and 300-bus under normal operations

Table 7
Optimum locations of PMUs under normal operations

Table 9
Comparison of obtained results with existing methods without a ZIB

Table 10
Comparison of obtained results with existing methods with ZIB

Table 11
Comparison of computational time with existing methods ignoring and considering ZIBs for large-scale power grids.The results indicate that the proposed algorithms outperform other meta-heuristic algorithms available in the state-of-the-art literature.The suggested method may be developed further in a number of different ways that may be researched.These include, among others things: (a) Performance Assessment: Compare the proposed BAPSO technique to other optimization algorithms utilized for PMU allocation.To illustrate the usefulness and efficiency of the suggested technique, this assessment should encompass a wide range of test systems with different sizes and complexity.(b) Resilience Analysis: Evaluate the suggested allocation method's resilience by taking into account power system parameter uncertainties including demand fluctuations, line outages, and generator failures.Look at the PMU allocation scheme's capacity to adjust to such dynamic events and provide dependable observability under challenging circumstances.(c) Network Topology Incorporation: Look into incorporating network topology limi- "-" means not reportedIg.and Cons.ignoring and considering; Comp.computational

Table 13
Buses connected with ZIB for all test systems

Table 13 (
continued)Maximum absolute value of velocity allowed along the same jth direction in the parameter space M max Maximum number of iterations used in PSO w max Personal best position for particle i discovered so far pop Population size r 1 and r 2Random numbers that are uniformly distributed between [0, 1] to maintain the swarm diversity R 1 PMUTransformation matrix that may be modified according to the contingency cases W N obs Total number of a bus which is observable c iVector of PMU price coefficients v ij (m + 1)