Next Article in Journal
The Effects of Virtual Reality-Assisted Language Learning: A Meta-Analysis
Previous Article in Journal
Constructed Wetlands as Nature-Based Solutions in the Post-COVID Agri-Food Supply Chain: Challenges and Opportunities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Weak Buses for Optimal Load Shedding Using Differential Evolution

by
Olumuyiwa T. Amusan
,
Nnamdi I. Nwulu
and
Saheed Lekan Gbadamosi
*
Center for Cyber Physical Food, Energy and Water Systems, University of Johannesburg, Johannesburg 2006, South Africa
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(6), 3146; https://doi.org/10.3390/su14063146
Submission received: 25 December 2021 / Revised: 12 February 2022 / Accepted: 14 February 2022 / Published: 8 March 2022
(This article belongs to the Section Energy Sustainability)

Abstract

:
For the sustainability of power supply and operation systems, planners aim to deliver power at an optimum value to consumers, while maintaining stability in the system. The load-shedding approach has proven to be an effective means of achieving the desired stability. This paper presents a nodal analysis to establish critical bus identification in the power grid. A power simulation for load shedding was created using the power system analysis toolbox (PSAT) for identifying and isolating weak buses on the power system. A computational algorithm was developed using differential evolution (DE) for minimizing service interruptions and blackouts, and was tested against the conventional genetic algorithm (GA). The proposed algorithm was implemented on an IEEE 30-bus test system. The simulation results were analyzed before and after the application of DE. It was observed that after the application of DE, load shedding gives an efficient result of 10.6%, 8.7%, and 13.4% improvement at buses 26, 29, and 30, respectively, after being tested using a genetic algorithm (GA), with a result of 10.2%, 7.6% and 13.1% on the same respective buses. This work will further expand the reliability and availability of power systems toward a sustainable, steady power supply that is void of nodal or bus cutoffs.

1. Introduction

The rise in industrialization and technological advancements around the world has necessitated the need for power utility industries to meet the rising demand of consumers at every point in time [1]. The increased size and capacity of renewables to meet rising consumer load demand, combined with a lower inertia and more sensitive voltage profiles, are making the system operators’ roles more complex. As a matter of fact, a mismatch between the power that is demanded and supplied may be a driver of system collapse and the consequent widespread blackouts [2]. As such, the operating margins of conventional power systems are now moving close to their limits, thereby creating increased voltage profile violations, power-balancing problems, and an overall grid maintenance crisis that could cause blackouts and outages. This was demonstrated by the US and Canada 96-hour blackout in 2003, which affected over 50 million people. In order to avoid a repeat of this situation, it becomes necessary to shed some of the load at specific buses, thus stabilizing voltages at steady-state values in all the buses of the distribution network. In this way, such a mismatch in consumer load consumption value and the amount of power generated will not affect the frequency of the system, necessitating load shedding [3,4].
Voltage instability gradually begins as a fall in the voltage profile, which eventually leads to unsafe levels at sensitive and critical nodes. If this is left unmodified, nodal system instability may arise at voltage profiles close to nominal values. This is the same situation as when the combined generation and transmission system is unable to fulfill the power demands required by the loads, thus shifting the voltage stability framework. This imbalance between power generation and demand power might cascade into system outage, reactive power source restrictions, or a combination of both. In severe cases, voltage instability can lead to voltage collapse, resulting in a blackout or abnormally low voltages in a substantial area of the power system.
Various load-shedding schemes have been proposed by researchers and engineers. Traditional load-shedding methods monitor frequency and voltage using under-frequency, activating load cut-offs until the frequency returns to a stable point of operation. This happens during situations when the active power available is incapable of meeting the load and there is a need to maintain frequency within safe limits. This functionality has necessitated a new voltage stability index, taking into account the load characteristics, to improve voltage stability [5]. In one study [6], wind generators were used to maintain voltage stability using a load impedance to source impedance ratio to calculate the number of induction motors that must be dropped when there are dominant loads in the feeders [7,8].
Evolutionary algorithms or their hybrids have also been employed to study load management. In a study by [9], a hybrid mix using the firefly algorithm with particle swarm optimization (PSO) was applied to an IEEE-33 bus system, to determine the maximum load to be shed from network operation, while maintaining the voltage profile of the system within allowable limits. The result shows the capability of the method of restoring the system’s frequency without overshooting. An analytical basis for under-voltage load shedding to mitigate voltage collapse in the power system was presented in [6], using a dynamic load model as the criterion by which to evaluate the amount of load to be shed. The authors of [10] proposed the implementation of the glowworm swarm optimization (GSO) algorithm to minimize load shedding, testing the proposed algorithm on IEEE 14, 30, and other bus test systems. This technique, when compared with other conventional methods, provided better results on evaluation. The authors of [11] suggested the application of a genetic algorithm (GA) to solve the nonlinear optimization problem of load shedding on the IEEE 30-bus test system. In binary combinatorial problems, the use of metaheuristic techniques is highly favorable in dealing with a large number of variables. Among them, the genetic algorithm (GA) technique is promising, due to its execution speed and considerable accuracy, which is the situation in the current problem, as will be validated in the case study.
This paper proposes a new algorithm for optimum load shedding using differential evolution(DE). DE has been found to offer a suitable solution method, due to its ability and effectiveness when solving joint replenishment and delivery (JRD) problems in the studies by [12,13], and as reported by [14]. This is further supported here since the parameters of this intelligent solution’s algorithms play an important role in its performance. A power simulation model was developed for the identification of weak buses, using the power system analysis toolbox (PSAT), and an eigenvalue was set as an index using sensitivity analysis to provide situational awareness. The developed algorithm is of two sections, the first section is used for identifying the load buses required for appropriate load shedding, while the other section decides which of the optimum loads is required to be shed at the chosen load buses, using differential evolution. Validation of the method is conducted using a genetic algorithm, a method that has often been proven by different authors to be fast and accurate [15,16,17,18]. The remaining of this paper is structured as follows: Section 2 describes the proposed methodology, which includes nodal analysis and a DE algorithm for optimizing the load shedding. Section 3 presents a numerical simulation of the load flow analysis and of the identification of weak load buses using the test system. In Section 4, the results are discussed, and Section 5 provides our conclusions.

2. Methodology

The process of carrying out optimized load shedding that is adopted in this work is divided into two parts. The weak buses in the system are first identified, before implementing a differential evolution algorithm. This model of the IEEE 30-bus system has been simulated using PSAT, and a Newton–Rhapson power flow solution was employed to ascertain the magnitudes of voltage variation and losses in the network, owing to its fast convergence characteristics.
For a k-bus power system with a nodal current and voltages I and V, the system forms an admittance matrix, Y, which is formed based on the apparent power, S, expressed as follows [19]:
S = I V = Y V
This can also be expressed in a complex relationship, as seen in Equation (2), given real power, P, and reactive power, Q, as follows:
S = P + j Q
Hence, in a k-bus system, the real and reactive power are individually given by:
P k = V k m = 1 n V k Y k m cos ( δ k δ m Φ k m )
Q k = V k m = 1 n V k Y k m sin ( δ k δ m Φ k m )
The power mismatch at the bus, k, is given by the following equations:
Δ P k = P k k s P k k c
Δ Q k = Q k k s Q k k c
where the superscripts ks and kc signify the specified and calculated values, respectively.
The power-flow elements of the network are used to construct and populate the Jacobian matrix. The network elements are added one after the other, and the Jacobian matrix is updated in a simple way, as obtained via the Newton–Raphson load flow (NRLF) model [20].
A Jacobian matrix, J, is formed from the Newton–Rhapson load flow obtained through Equation (7), from the active and reactive power mismatch ( Δ P and Δ Q ), and is given as:
Δ P Δ Q = J Δ Φ Δ V
where Δ P and Δ Q represent the mismatch in active and reactive power vectors; Δ Φ and Δ V represent the unknown voltage angle and magnitude at the buses.

2.1. Nodal Analysis and Participation Factor

Nodal analysis is used to determine the most vulnerable location in terms of voltage stability. This will help select the best site for locating the most effective actions to alleviate changing voltage conditions. A modal analysis equation can be used to compute the smallest eigenvalues and the respective eigenvectors of the voltage and reactive power variations in the system [21]. The results of the analysis help to determine the voltage stability performance of the system. The participation factor is used to find the weakest bus in the system.
Δ P Δ Q = J 11 J 12 J 21 J 22 Δ Φ Δ V
At constant power, Equation (8) becomes:
0 Δ Q = J 11 J 12 J 21 J 22 Δ Φ Δ V
Solving (9), we obtain the following:
Δ Φ = J 11 1 J 12 Δ V
Δ Q = J 21 Δ Φ + J 22 Δ V
Substituting (10) into (11) yields:
Δ Q = J r Δ V
where Jr is the reduced Jacobian matrix of the system. The eigenvalues and eigenvectors of the reduced Jacobian matrix of the system indicate the voltage stability characteristics of the system, the magnitudes of which are proportional to how close the system is to instability.
The reduced Jacobian matrix, Jr, can be found by solving the equation below:
J r = J 22 J 21 J 11 1 J 12
which results in the following equation:
J r 1 = Φ γ 1 ξ
where Φ and ξ represent the right and left eigenvalue, γ , which gives information about the observability and controllability of the system. For non-dynamic systems, positive eigenvalues of Jr indicate that the system is stable, while negative eigenvalues indicate voltage instability. In the case of a zero eigenvalue, the system is on the verge of voltage instability or has just reached its stability limit [22].

2.2. Differential Evolution

Differential evolution, as developed by Price and Storn, is a population-based optimization tool that has been applied in various engineering scenarios to propose a feasible solution. It solves the problem using an objective function at multiple randomly chosen initial points. As a method, it offers great convergence characteristics and needs only a few control parameters that do not need to be modified during the optimization process. It is governed by four main processes, as depicted in Figure 1.
  • Initialization
The first step of the DE algorithm is to initialize the problem parameters, i.e., the initial population, which must lie within the feasible bounds of the decision variables. This can be generated by following Equation (15):
X i , j G = 0 = X j l + λ j ( X j u X j l )
where i represents 1, 2, 3, …, Np, and Np is the vector population size; j represents 1, 2, 3, …, D, and D stands for the dimension of the decision variable. X j u a n d X j l represent the upper and lower bounds of the decision variable, respectively; λ j represents the uniformly distributed random number, lying between 0 and 1. For each new j, X i , j 0 is the initial value of the jth parameter of the ith decision variable.
  • Mutation
Three random vectors are selected from the vector population, i.e., X a , X b   a n d   X c , such as a b c i , where X i ( G ) is the already chosen target vector from the initialization step. Two of these chosen random numbers (say, X a , X b ) have their weighed difference perturbed with the third random vector X c , to give a mutant vector V i , j G , such that:
V i , j G = X c + α ( X b G X a G )
where V i , j G = [ V i , 1 ( G ) , , V i , D ( G ) ] is the mutant vector and α is the scaling factor, which is usually chosen to be between 0 and 2.
  • Crossover
The crossover step is essential to increase the potential diversity of the parameter vectors. This step generates the trial vector, U i G , by mixing the target vectors, X i G with the mutant vector, V i G . Crossing over is achieved by satisfying any of the following two schemes, as given by:
U i , j ( G ) = V i , j ( G ) i f λ j C r   o r   j = q X i , j ( G ) o t h e r w i s e
where λ j is a randomly newly chosen vector between 0 and 1 for each value of j . Cr is the crossover rate (or probability), which is user-defined and usually lies between 0 and 1.
The index, q, is chosen to be between 1 and D to ensure that the trial vector gets at least one of its parameters from the mutant vector. If X j is less than or equal to C r , the trial vector is adopted from the mutant vector; otherwise, the trial vector comes from the target vector, ensuring that the trial vector does not completely duplicate the target vector.
  • Selection
In this stage, the vectors that would populate the next generation are selected from either the trial vectors or the target vectors. The fitness value of the trial vectors and target vectors are compared, and the best solution is selected according to the following equation:
X i , j ( G + 1 ) = U i , j ( G ) i f   f ( U i , j ( G ) ) f ( J ) X i , j ( G ) o t h e r w i s e
where f ( J ) is the function to be minimized, defined as:
J = i = 1 N L S l s i
l s i = min { i = 1 N L S ( P L i f i ( x min ) f i ( x max ) }
where lsi is the total load shed at the i th bus, N L S is the set of load buses selected for load shedding, P L i is the i th load-shedding bus, f i ( x min )   and   f i ( x max ) are the minimum and maximum permissible load shed at the selected i th bus. This is subjected to any voltage constraints as given in Equations (21) and (22). The complete algorithm for shedding the load is depicted in Figure 2.
0.05 V     V   0.05 V
π     Φ   π

2.3. Genetic Algorithm Solution Space

The genetic algorithm is a method used to search for the optimal solution by simulating the natural evolutionary process. It involves encoding, population initialization, fitness function, genetic operators, and elite reservation strategy [23], which can be summarized as follows:
  • Encoding. Encoding is the representation of chromosomes (different voltages at the nodes), which establishes the mapping between the solution space of the original problem and the environment of the population in GA. The form of encoding is related to the characteristics of the original problem’s decision variables.
  • Population initialization. Population initialization is used to obtain initial solutions and prepare for the subsequent evolution.
  • Fitness function. Fitness function is designed according to the original problem’s objective function, and the fitness function value of a chromosome determines the probability of its retention.
  • Genetic operators. Genetic operators include the selection operator, crossover operator, and mutation operator. The selection operator selects chromosomes from the current population to enter the next generation population.
In order to visualize the optimization process, the multivariate power problem is decomposed, and the optimization process values of each variable are saved and displayed in sequence. In the operation of the genetic algorithm, new individuals are constructed manually and are added to the chromosome population to replace the poorer-fitting individuals in the population. When the genetic algorithm makes an optimization selection, all optimization variables are coded as search objects, which enhances the scope of the optimization selection [24]. A PSAT simulation of the IEEE 30-bus network structure is shown in Figure 3. With the increase in node relationships, the shortest path between nodes and the average distance of the whole network are decreasing, thus gradually drawing the system into an unstable state.

3. Numerical Simulation

The power system analysis toolbox (PSAT) is employed in this work to model and perform load flow steady-state analysis using the Newton–Raphson model.
It is used for performing power flow analysis and to establish continuation power flow, optimal power flow, and small-signal stability, as well as simulations in the time domain. This aids in evaluating the system operating conditions, evolving due to small and gradual changes, thus making the computation of power flow analysis simpler and more explicit [25]. This method was used to model the IEEEE 30-bus system shown in Figure 3. Eigenvalue analysis is performed using PSAT, to obtain the most critical load buses and participation factors that can take the system toward a state of near-voltage collapse. With this result, a load-shedding decision can be made, based on buses with low participation factors.

4. Results and Discussion

The proposed DE algorithm was evaluated using an IEEE 30-bus test system. In this study, equal loading was applied to the test bus system, with the set indicator reduced to a small value; the system shows a significant breach of bus voltages. Under these conditions, a simulation was carried out using the power system analysis toolbox (PSAT) to model and perform load flow steady-state analysis with Newton–Raphson.
The test bus system consists of 30 buses, six generation buses, 41 transmission branches, and 20 load buses. The voltage of the load buses ranges from 0.95 to 1.05 p.u. Table 1 presents a power flow solution that contains the bus voltage, system load, and PQ buses. Table 2 shows the results of the eigenvalue for the test bus system. It can be observed that bus 30 has the smallest eigenvalue of 0.11052, which is an indication that bus 30 is the most critical bus in the system. Table 3 presents the participation factor for each bus. It can be observed that the bus with the highest participation factor is bus 30, with a participation factor value of 0.4251. This is an indication that bus 30 makes the highest contribution to voltage collapse in the system. Therefore, this is referred to as the weakest bus in the system. Similarly, buses 26 and 29 have both been identified as contributing to voltage collapse on the system, as signified by their participation factors of 0.2545 and 0.2287, respectively. Therefore, the load can be shed at these three buses, based on the load-shedding function.
Table 4 presents the minimum load at critical-load buses, based on two scenarios: (a) before and after load shedding, without the application of optimization; (b) with optimized load shedding after application of the DE optimization technique. Table 5 shows the voltage profile of the grid network after load shedding. Table 6 presents the percentage improvement of bus voltages using DE and GA, after load shedding. It can be observed that the critical buses, 26, 29, and 30, experienced an increase in the voltage profile, with 10.2%, 7.6%, and 13.1%, respectively, using GA. This was shown to be better, with an increased percentage improvement of 10.6%, 8.7%, and 13.4%, respectively, using DE after load shedding was performed on the network.

5. Conclusions

This study has presented a novel algorithm for load-shedding optimization of the required load buses in a 30-bus IEEE test network. The following has been achieved, and the conclusions herein are as follows:
  • A simulation was performed on the grid, using PSAT for nodal load analysis to identify the critical load buses.
  • A computational algorithm was developed for load shedding, based on the sensitivities of the load bus using differential evolution, and was validated using a genetic algorithm.
  • The results proved that DE significantly improved the voltage profile of the test bus system after load shedding, with a percentage improvement of 10.6, 8.7, and 13.4 on buses 26, 29, and 30, respectively, as opposed to 10.2, 7.6, and 13.1 on the same buses, using GA.
  • The system network also experienced an appropriate increment in voltage values per unit at the selected load buses and other subsequent buses.

Author Contributions

Conceptualization, O.T.A.; methodology, O.T.A. and S.L.G.; software, O.T.A.; validation, O.T.A., S.L.G. and N.I.N.; formal analysis, O.T.A. and S.L.G.; investigation, O.T.A.; resources, O.T.A.; data curation, O.T.A.; writing—original draft preparation, O.T.A.; writing—review and editing, S.L.G. and N.I.N.; visualization, S.L.G.; supervision, S.L.G. and N.I.N.; project administration, N.I.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used is available in the PSAT documentation software by Milano, F. Version 2.0.0.

Acknowledgments

The authors thank the University of Johannesburg, South Africa for providing an enabling top-class academic environment for the success of this research work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gbadamosi, S.L.; Nwulu, N.I. A multi-period composite generation and transmission expansion planning model incorporating renewable energy sources and demand response. Sustain. Energy Technol. Assess. 2020, 39, 100726. [Google Scholar] [CrossRef]
  2. Gbadamosi, S.L.; Nwulu, N.I. A comparative analysis of generation and transmission expansion planning models for power loss minimization. Sustain. Energy Grids Netw. 2021, 26, 100456. [Google Scholar] [CrossRef]
  3. Gbadamosi, S.L.; Nwulu, N.I. Reliability assessment of composite generation and transmission expansion planning incorporating renewable energy sources. J. Renew. Sustain. Energy 2020, 12, 026301. [Google Scholar] [CrossRef]
  4. Gbadamosi, S.L.; Nwulu, N.I.; Sun, Y. Multi-objective optimisation for composite generation and transmission expansion planning considering offshore wind power and feed-in tariffs. IET Renew. Power Gener. 2018, 12, 1687–1697. [Google Scholar] [CrossRef]
  5. Deng, J.; Liu, J. A Study on a Centralized Under-Voltage Load Shedding Scheme Considering the Load Characteristics. Phys. Procedia 2012, 24, 481–489. [Google Scholar] [CrossRef] [Green Version]
  6. Arief, A.; Dong, Z.; Nappu, M.B.; Gallagher, M. Under voltage load shedding in power systems with wind turbine-driven doubly fed induction generators. Electr. Power Syst. Res. 2013, 96, 91–100. [Google Scholar] [CrossRef] [Green Version]
  7. Larik, R.; Mustafa, M.W.; Aman, M.; Jumani, T.; Sajid, S.; Panjwani, M. An Improved Algorithm for Optimal Load Shedding in Power Systems. Energies 2018, 11, 1808. [Google Scholar] [CrossRef] [Green Version]
  8. Wiszniewski, A. New criteria of voltage stability margin for the purpose of load shedding. IEEE Trans. Power Deliv. 2007, 22, 1367–1371. [Google Scholar] [CrossRef] [Green Version]
  9. Jallad, J.; Mekhilef, S.; Mokhlis, H.; Laghari, J.; Badran, O. Application of Hybrid Meta-Heuristic Techniques for Optimal Load Shedding Planning and Operation in an Islanded Distribution Network Integrated with Distributed Generation. Energies 2018, 11, 1134. [Google Scholar] [CrossRef] [Green Version]
  10. Mageshvaran, R.; Jayabarathi, T. GSO based optimization of steady state load shedding in power systems to mitigate blackout during generation contingencies. Ain Shams Eng. J. 2015, 6, 145–160. [Google Scholar] [CrossRef] [Green Version]
  11. Raghu, C.; Aravind, P.M. Assessing Effectiveness of Research for Load Shedding in Power System. Int. J. Electr. Comput. Eng. 2017, 7, 3235–3245. [Google Scholar]
  12. Cui, L.; Wang, L.; Deng, J.; Zhang, J. Intelligent algorithms for a new joint replenishment and synthetical delivery problem in a warehouse centralized supply chain. Knowl.-Based Syst. 2015, 90, 185–198. [Google Scholar] [CrossRef]
  13. Wang, L.; Dun, C.X.; Bi, W.J.; Zeng, Y.R. An effective and efficient differential evolution algorithm for the integrated stochastic joint replenishment and delivery model. Knowl.-Based Syst. 2012, 36, 104–114. [Google Scholar] [CrossRef]
  14. Qu, H.; Ai, X.Y.; Wang, L. Optimizing an integrated inventory-routing system for multi-item joint replenishment and coordinated outbound delivery using differential evolution algorithm. Appl. Soft Comput. 2020, 86, 105863. [Google Scholar] [CrossRef]
  15. Ajmal, M.S.; Iqbal, Z.; Khan, F.Z.; Ahmad, M.; Ahmad, I.; Gupta, B.B. Hybrid ant genetic algorithm for efficient task scheduling in cloud data centers. Comput. Electr. Eng. 2021, 95, 107419. [Google Scholar] [CrossRef]
  16. Liu, Y.; Ćetenović, D.; Li, H.; Gryazina, E.; Terzia, V. An optimized multi-objective reactive power dispatch strategy based on improved genetic algorithm for wind power integrated systems. Int. J. Electr. Power Energy Syst. 2022, 136, 107764. [Google Scholar] [CrossRef]
  17. Yuan, H.-B.; Zou, W.-J.; Jung, S.; Kim, Y.-B. Optimized rule-based energy management for a polymer electrolyte membrane fuel cell/battery hybrid power system using a genetic algorithm. Int. J. Hydrogen Energy 2022, 47, 7932–7948. [Google Scholar] [CrossRef]
  18. Ramos-Figueroa, O.; Quiroz-Castellanos, M.; Mezura-Montes, E.; Kharel, R. Variation Operators for Grouping Genetic Algorithms: A Review. Swarm Evolut. Comput. 2021, 60, 100796. [Google Scholar] [CrossRef]
  19. Kothari, D.P.; Nagrath, I.J. Modern Power System Analysis; McGraw Hill Education: New York, NY, USA, 2003. [Google Scholar]
  20. Jegatheesan, R.; Nor, N.M.; Romlie, M.F. Newton-Raphson power flow solution employing systematically constructed jacobian matrix. In Proceedings of the 2nd IEEE International Conference on Power and Energy (PECON 08), Johor Bahru, Malaysia, 1–3 December 2008. [Google Scholar]
  21. Konoval, V.; Prytula, R. Participation Factor in Modal Analysis of Power System Stability. Pozn. Univ. Technol. Acad. J. 2016, 86, 97–104. [Google Scholar]
  22. Enemuoh, F.O.; Onuegbu, J.C.; Anazie, E.A. Modal based analysis of and Evaluation of Voltage Stability of Bulk Power System. Int. J. Eng. Res. Dev. 2013, 6, 71–79. [Google Scholar]
  23. Tao, Z.; Moncada, J.A.; Poncelet, K.; Delarue, E. Review and analysis of investment decision making algorithms in long-term agent-based electric power system simulation models. Renew. Sustain. Energy Rev. 2021, 136, 110405. [Google Scholar] [CrossRef]
  24. Xu, L. Research on computer interactive optimization design of power system based on genetic algorithm. Energy Rep. 2021, 7, 1–13. [Google Scholar] [CrossRef]
  25. Milano, F. Documentation for PSAT, Version 2.0.0; University of Castilla-La Mancha: Ciudad Real, Spain, 2008. [Google Scholar]
Figure 1. Diagram of the DE optimization procedure.
Figure 1. Diagram of the DE optimization procedure.
Sustainability 14 03146 g001
Figure 2. Flowchart algorithm for load shedding using DE.
Figure 2. Flowchart algorithm for load shedding using DE.
Sustainability 14 03146 g002
Figure 3. Scheme of the IEEE 30-bus system simulation, using PSAT.
Figure 3. Scheme of the IEEE 30-bus system simulation, using PSAT.
Sustainability 14 03146 g003
Table 1. Initial power flow solution.
Table 1. Initial power flow solution.
BusV (p.u)Phase (p.u)PG (p.u)QG (p.u)PL (p.u)QL (p.u)
1100.5695−0.17100
21−0.01230.8−0.00310.2170.127
30.994−0.0467000.0240.012
40.9918−0.0555000.0760.016
51−0.08750.50.37380.9420.19
60.9885−0.07130000
70.9838−0.0848000.2280.109
81−0.07410.350.4180.30.3
90.9781−0.08990000
100.9575−0.138000.058−0.1542
111−0.05290.30.084600
120.9849−0.0654000.1120.075
131−0.05350.40.25800
140.9321−0.1189000.0620.016
150.904−0.1213000.0820.025
160.9513−0.1045000.0350.018
170.9344−0.1484000.090.058
180.7854−0.1631000.03080.0867
190.8211−0.2052000.0950.034
200.8519−0.1918000.0220.007
210.925−0.1652000.1750.112
220.9249−0.16470000
230.8459−0.1855000.0320.016
240.835−0.2208000.0870.0391
250.886−0.19570000
260.7543−0.2577000.03110.0204
270.984−0.14950000
280.9862−0.08020000
290.818−0.2999000.0240.009
300.7301−0.4294000.08830.0158
V: Voltage; PG: Active Power; QG: Reactive Power; PL: Active Load; QL: Reactive Load.
Table 2. Eigenvalue result for an IEEE 30-bus system.
Table 2. Eigenvalue result for an IEEE 30-bus system.
EigenvalueMost AssociatedReal PartImaginary Part
14101.25150
2668.81450
31235.20220
42130.77520
51026.78560
62823.21850
7718.25070
8913.83030
9209.01770
1038.37850
11154.50320
12174.34970
13273.68120
14162.43360
15161.94240
16241.95310
17141.54990
18251.36710
19300.11050
20260.16910
21180.33760
22231.11360
23230.62990
24290.77440
25119990
26139990
2729990
2859990
2989990
30Slack Bus 19980
Table 3. Participation factors of buses in IEEE 30-bus test system.
Table 3. Participation factors of buses in IEEE 30-bus test system.
BusParticipation Factor
Slack bus0
20
30
40
50
60.0001
70.0001
80
90.0001
100.0003
110
120
130
140.0001
150.0003
160.0001
170.0003
180.0009
190.0008
200.0007
210.0007
220.0008
230.0038
240.0103
250.0514
260.2287
270.0205
280.0006
290.2545
300.4251
Table 4. Minimum total load shed, before and after applying the DE algorithm.
Table 4. Minimum total load shed, before and after applying the DE algorithm.
BusInitial Active PowerMinimum Load ShedMaximum Load ShedTotal Load ShedActive Power after Optimized Load Shedding (p.u.)Reactive Power after Optimized Load Shedding (p.u.)
260.0350.0070.02450.00350.03150.0207
290.0240.00480.01680.00240.02160.0081
300.1060.02120.07420.01060.09540.0171
Table 5. Active power (P), reactive power (Q), and phase, after load shedding using DE.
Table 5. Active power (P), reactive power (Q), and phase, after load shedding using DE.
BusPhase (p.u.)PG (p.u.)QG (p.u.)PL (p.u.)QL (p.u.)
Slack Bus00.453−0.153900
2−0.0090.8−0.03140.2170.127
3−0.0396000.0240.012
4−0.047000.0760.016
5−0.0810.50.36060.9420.19
6−0.0610000
7−0.0755000.2280.109
8−0.06190.350.36740.30.3
9−0.07220000
10−0.1175000.058−0.1579
11−0.03410.30.051800
12−0.049000.1120.075
13−0.03660.40.219600
14−0.0989000.06060.0156
15−0.1004000.0770.0235
16−0.0868000.0350.018
17−0.128000.08890.0573
18−0.1356000.02460.0692
19−0.1686000.0780.0279
20−0.159000.01910.0061
21−0.1424000.17080.1093
22−0.14150000
23−0.1548000.02730.0137
24−0.1859000.0740.0263
25−0.15960000
26−0.2014000.02430.016
27−0.12130000
28−0.06810000
29−0.231000.01890.0071
30−0.3177000.07250.013
Table 6. Comparison of the percentage improvement of bus voltages, using DE and GA.
Table 6. Comparison of the percentage improvement of bus voltages, using DE and GA.
BusV (before Load Shedding) (p.u.)V (after Load Shedding) (p.u.)Percentage Improvement Using GAPercentage Improvement Using DE
Slack Bus1.00001.000000
21.00001.000000
30.99400.99610.20050.2062
40.99180.99430.21010.2563
51.00001.000000
60.98850.99180.26890.3312
70.98380.98630.21800.2484
81.00001.000000
90.97810.98250.38970.4441
100.95750.96761.04351.0501
111.00001.000000
120.98490.98630.05640.1484
131.00001.000000
140.93210.93900.70410.7441
150.90400.92051.01181.8322
160.95130.95570.44430.4646
170.93440.94401.00001.0246
180.78540.83304.63876.0636
190.82110.86082.85864.8275
200.85190.88481.51743.8585
210.92500.93851.41561.4542
220.92490.93921.36921.5564
230.84590.87821.43973.8156
240.83500.87623.14724.9397
250.88600.92672.11684.5927
260.75430.834210.200010.5823
270.98401.00432.01242.0583
280.98620.99140.45960.5240
290.81800.88937.57068.7188
300.73010.828013.104013.4060
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Amusan, O.T.; Nwulu, N.I.; Gbadamosi, S.L. Identification of Weak Buses for Optimal Load Shedding Using Differential Evolution. Sustainability 2022, 14, 3146. https://doi.org/10.3390/su14063146

AMA Style

Amusan OT, Nwulu NI, Gbadamosi SL. Identification of Weak Buses for Optimal Load Shedding Using Differential Evolution. Sustainability. 2022; 14(6):3146. https://doi.org/10.3390/su14063146

Chicago/Turabian Style

Amusan, Olumuyiwa T., Nnamdi I. Nwulu, and Saheed Lekan Gbadamosi. 2022. "Identification of Weak Buses for Optimal Load Shedding Using Differential Evolution" Sustainability 14, no. 6: 3146. https://doi.org/10.3390/su14063146

APA Style

Amusan, O. T., Nwulu, N. I., & Gbadamosi, S. L. (2022). Identification of Weak Buses for Optimal Load Shedding Using Differential Evolution. Sustainability, 14(6), 3146. https://doi.org/10.3390/su14063146

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop