Power Prioritization and Load Shedding in an Island with RESs Using ABC Algorithm

Department of Electrical Engineering, Pan African University Institute for Basic Sciences, Technology and Innovation (PAUSTI), P.O. Box 62000-00200, City Square, Nairobi, Kenya Department of Electrical Engineering, Jomo Kenyatta University of Agriculture and Technology (JKUAT), P.O. Box 62000-00200, City Square, Nairobi, Kenya Department of Electrical Engineering, Technical University of Mombasa, P. O. Box 90420-80100, Mombasa, Kenya


Introduction
Globally, energy demand has been steadily increasing as a result of increase in population and industrial and economic growth in various sectors [1]. Hence, in order to provide energy to this steady increase in power demand, clean RESs distributed generators (DGs) are being installed to operate either in grid connected or islanded modes. ese sources are preferred because they are environmentally friendly and assist in reducing the cost of transmission line construction. However, these DGs may not adequately meet the load requirements within the island in case of a contingency that can lead to loss of mains.
In order to ensure there is continuity of supply of power to important loads and maintain the system stability, utility companies resort to optimal shedding of loads within the island based on the given criteria guided by IEEE 1547 standard [2,3]. Optimal shedding of some loads within the island is used to minimize the difference between the connected active and reactive load and the power that the connected RESs can supply [4]. rough the load shedding action, the disturbed MG is made to settle to a new different equilibrium state. Undervoltage Load Shedding (UVLS) and Underfrequency Load Shedding (UFLS) are the criteria mainly used for load shedding [5]. rough the shedding of some loads, the perturbed system is forced to settle to a new state of equilibrium [6]. However, for any method of load shedding to be successful, the amount to be shed should be optimal, at the correct location and be done on the right time so as to prevent any possibility of voltage instability and also avoid losing the customer trust in the utility company [7,8]. For improved accuracy in the determination of the amount of load to be shed in a given island, real-time monitoring and control of the system is preferred. However, this is a costly undertaking especially for small grids.
A number of studies have been done on UFLS as can be seen in [9][10][11]. In this method, frequencies are measured at a point of interest and compared with the set thresholds. If the values exceed the limits, then some loads will be shed to equalize power generation and supply. It is shown that automatic UVLS performs far much better than manual UVLS during emergencies in [12]. In [13], a Comprehensive Learning Particle Swarm Optimization (CLPSO) method in the optimal subdivision of a section of the grid when there is loss of mains is suggested. Power supply and demand are then balanced by optimally shedding some loads from each islanded subsection. In [14], optimal load shedding is achieved using Artificial Neural Network (ANN) algorithm. In [15], loads were classified into noncritical and critical loads before optimally shedding them in case generation is less than loads. However, this method takes a lot of time to complete the load shedding process as the algorithm has to go through the two classes of loads. In [5], the author used a combination of Particle Swarm Optimization (PSO) and Genetic Algorithm (GA) to find an optimal amount of load to be shed in smart grids under stress. e sum of squares of the difference between the connected loads and the supplied power and ABC algorithm to arrive at the optimal amount of load to be shed during a contingency is proposed in [6]. However, in this case, the priority of important loads and buses was not considered in the study before shedding of the loads. In [16], the author evaluated an implementation of optimal UVLS in competitive electricity markets using Particle Swarm Optimization. In [17], the authors developed an optimal load shedding model using Particle Swarm Optimization in order to control voltage instability in a system. However, this approach has a long convergence time since it comprises of multiobjectives. A hybrid of GA and PSO is proposed in [18] for load shedding. However, this method leads to suboptimal load shedding since it does not consider voltage instability constraint. ere is continued research being undertaken on the usage of computational intelligence algorithms such as fuzzy logic, genetic algorithms, particle swarm optimization among others. e merits and drawbacks of these methods are briefly discussed in [9]. ere is need for additional studies that involves high penetration of RESs with the aim of minimizing the total amount of load to be shed during a contingency. ere is also need to improve the voltage profile of the system and improve the security of the power system so as to reduce the possibilities of the system being overloaded. is is because a technology that sheds more than or less than the required optimal load amount will not control the instability of the voltage but can instead lead to the collapse of the system voltage and frequency [5].
In case of a contingency that leads to a section of the grid operating as an island, the isolated generators should operate in frequency and voltage control mode to provide constant voltage to local loads so as to maintain the voltage and frequency within the island. Otherwise, this may affect the power system operation and pose safety risks to the maintenance personnel [19]. In this paper, loads connected to the different buses are assigned priority indices based on Fast Voltage Stability Index (FVSI). en, ABC algorithm is used to evaluate the amount of load to be shed to ensure there is balance between power generated and connected loads within an islanded microgrid. e rest of the paper is structured as follows. Section 2 introduces the formulation of the problem for the proposed load prioritization and shedding scheme using ABC algorithm, while Section 3 discusses the methodology used in this approach. e study simulation results and their discussions are presented in Section 4. And finally, the paper conclusion is highlighted in Section 5. Figure 1 is the structure of an islanded MG that was used in this study. It comprises RESs supplying loads locally, while the rest of the system is connected to the main grid through a transformer.

Load Prioritization.
In this study, loads and buses were assigned priority indices which determines their order of shedding in case of a contingency. Fast Voltage Stability Index (FVSI) was used to determine the weak buses and rank them in order of their priority as expressed by the following equation: where X is the reactance of the line, V 2 i is the square of the sending end voltage, Z is the impedance of the line, and Q r is the receiving end reactive power.
FVSI is an essential tool that helps in identifying weak buses and lines that are critically loaded and how close the system is to voltage collapse [20]. It is proven to be a good indicator for voltage stability in a power system. It is a sensitive index and hence an indicator of the best location to perform load shedding.
is information is useful when making decisions on which loads and buses to shed in case of a contingency. e higher the value of FVSI, the weaker the bus or line. Other indices include New Voltage Stability Index (NVSI), LQP, and Lmn. FVSI is the simplification of Lmn, while LQP is the simplification of NVSI. FVSI, Lmn, and LQP are found to be sensitive to the changes in reactive power, while NVSI is prone to the changes in real power load [21].

Load Shedding.
In this study, ABC algorithm was applied in shedding loads. is system is designed to calculate power imbalance in the island and then determine the minimum amount of load to shed and ensure as many loads as possible remain connected to the system. e main objective when microgrids operating as islands due to loss of mains is to ensure the supplied power and connected loads are equal, hence adequately supplied. is can be expressed by the following function: where P is the power to be shed, P di and Q di are the active and reactive powers that are supplied to the load, N is the total number of the buses in the system, P di and P di are active and reactive loads connected, and α i and β i are weights associated with bus priority in the system. At the same time, the power losses in the system should be minimal as expressed by the following equation: where G k is the conductance of the k th branch connected between the i th and j th buses. e system equality and inequality constraints can be expressed by (4)- (8). Equations (4) and (5) are equality constraints for active and reactive powers, whil (7) and (8) are constraints for minimum and maximum powers in the system: where G ij is the conductance between bus i and j, V gi is the bus voltage, P gi and Q gi are the active and reactive power generated at bus i, and P G−D and Q G−D are the active and reactive power flow between generation and demand points.

ABC Algorithm.
Artificial Bee Colony (ABC) algorithm is based on the behavior of the bees colony when looking for quality food sources (nectar). It is a metaheuristic optimization algorithm that was proposed by Karaboga from Erciyes University of Turkey [22]. e algorithm has been applied to solve a number of optimization problems, especially in engineering [23]. is is because it is robust and can be implemented very easily. Unlike other stochastic methods of optimization, this algorithm has excellent performance in exploration. However, it performs poorly when it comes to exploitation. is algorithm is divided into three groups as was observed from the behavior of a swarm of bees as follows.
(i) e Employed Bees Group. is group is equal to the possible solutions of the problem. ey continuously update the rest of the bees in the hive about the quantity, quality, and the direction of the food source through the performance of a waggle dance. e duration of the waggle dance depends on the quality and quantity of food source. (ii) e Onlooker Bees Group. ese are a group of bees which pick on a food source to exploit based on the information provided by the employed bees waggle dance. (iii) e Scout Bees Group. e work of this group is to look for new sources of food for exploitation.

ABC Control Parameters.
e parameters of this algorithm are set as follows: (i) NP: this is the total number of the bee population.
(ii) e total number of the available sources of food is set as NP/2, that is, half of the bee colony size. (iii) Limit: this is the number of trials for a given source of food before an employed bee abandons it and becomes a scout bee. (iv) Max cycle: this is the number of foraging cycles before an algorithm stops. e algorithm can be summarized into four steps as elaborated below.
(1) Initialization of the artificial bee stage where the bee colony size, the number of food sources, and limit of iterations parameters are set. e fitness of a food source can be expressed by the following equations:

Journal of Engineering
(2) e employee bee stage where a new fitness value is evaluated and compared with the previous one. e two values are compared and the inferior one is discarded. e new food sources are arrived at using the following equation: All food sources have the potential of being selected based on the fitness value of employed bee i, j as shown by the following equation: (3) e onlooker bee stage where the onlooker bees decide which food sources to exploit based on the probability as in stage 2 above. (4) e scout bee stage where if the food source does not improve after a set number of trials, that food source is left and the bee becomes a scout. e scout bee then goes out to search for new food sources randomly using the following equation: (5) e food with highest value of fitness is selected at this stage and printed.

Methodology
is proposed prioritization and optimal load shedding for MGs with RESs using ABC algorithm was done on the standard IEEE 30-bus system. is system has six generators located at buses 1, 2, 5, 8, 11, and 13; two static VAR sources are located at buses 10 and 24, thirty buses, forty-one lines, four tap changing transformers at lines 6-9, 6-10, 4-12, and 28-27 branches and twenty-one loads. It has a total of 283.400 MW and 126.200 MVAR base load.
First, the connected loads and buses were assigned priority indices to be used to determine which loads to shed based on their priority. is was based on FVSI values. Using ABC algorithm, the optimal amount of load to be shed within an island was determined. Both loss of generation and overload contingencies were analyzed using this algorithm. In each of the test cases, Newton-Raphson load flow was done to give the status of the bus voltage and line currents.
e simulation was carried out on MATLAB/SIMU-LINK platform on an AMD 4C + 6G 2.10 GHz processor with 4 GB RAM. e flow chart in Figure 2 was used in this study.

Simulation Results and Analysis
Using FVSI, the weak buses with their maximum loading limit were identified and arranged as shown in Table 1.

Generation Contingencies.
Loss of generation scenarios ranging from 0 to 75 MW at intervals of 5 MW were simulated while keeping the load demand constant at 328.4 MW. e line losses, convergence characteristics, and the amount of load to be shed by ABC algorithm at each instant were observed. is is as shown by Figures 3-7. Figure 3 shows the system loss for the 30-bus system for a generation deficit contingency. e active power generation is varied from 335 MW to 260 MW at intervals of 5 MW. From the observation, this approach shows a considerable reduction in power losses. However, it shows an increase in losses as the generation deficit increases.
Using this algorithm, the suggested amount of load to be disconnected so as to bring the system to a stable operating point was evaluated. is was recorded as shown in Figure 4. Figure 5 shows the convergence characteristic for a 283.4 MW load with 190 MW generation. e system converges after 23 iterations.
Keeping the load demand constant at 328.4 MW, the amount of load to be shed under various generation contingencies and their appropriate buses were monitored. is is as shown in Table 2. Figure 6 is the generation contingency characteristics when power generation is reduced from 322 to 290 MW. From the observation, the suggested load to be shed increases proportionally from 14.844 MW to 46.002 MW. e line losses decrease as the generation contingency increases. is is evidenced in Figure 7 where the losses are 8.484 MW when the power of 322 MW is being produced and reduces to 7.602 at 290 MW power generation.

Overload Contingencies.
To test the performance of the proposed approach in load shedding, the system was subjected to an overloading contingency of 85.02 MW on weak buses just like the case in [8]. is approach suggests a total of 30.491 MW of load to be shed. is is shown in Figure 8. is indicates that the approach performs better than particle swarm optimization (PSO) but not as good as genetic algorithm (GA) and GA-PSO hybrid working on the same conditions. GA-PSO hybrid exploits the advantages of both GA and PSO and overcoming of individual drawbacks of each of these algorithms applied separately. Hence, it gives the  lowest amount of load to shed in this analysis. Similarly, it is expected that the overall efficiency will be greatly improved with the GA-ABC hybridization and hence have even less amount of load to be shed when put into the same application. Figure 9 shows the voltage magnitude characteristics during the contingency and after load shedding using ABC algorithm in a standard IEEE 30-bus system. ese are compared against the voltage characteristics for normal operating condition. is shows that load shedding using ABC algorithm can reliably restore deteriorating voltages close to normal operating conditions. Figure 10 shows an observation of the line active and reactive power losses before, during the overload contingency, and after load shedding. When compared to the losses recorded by [6], this approach records lower losses, which is in line of minimizing system losses and load to be shed. From the figure, a total 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 Voltage magnitude

Bus number
Voltage magnitude during normal operating condition Voltage magnitude during overload contigency Voltage magnitude a er load shedding   In order to shed as minimal load as possible based on the priority attached, Table 3 shows the amount of load to be shed from respective buses based on their bus priority index.

Conclusion
In this study, load prioritization and optimal shedding using ABC algorithm have been presented for islands with RESs. IEEE 30-bus system was used to show the possibility of this approach and its validation. As it can be seen from these results, the proposed load shedding strategy was applied satisfactorily, and at the same time the objective function was minimized while satisfying the system's operational constraints. e results obtained also show that as minimal power as possible can be shed in case of a contingency, that can lead to islanding operation.
is is due to superior convergence characteristics of ABC algorithm in optimal load shedding in MGs with high penetration of RESs. It can, therefore, be concluded that this approach can be used in the determination of minimal amount to be shed to prevent voltage collapse and frequency instability.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper. Journal of Engineering 9