ORIGINAL RESEARCH article

Front. Energy Res., 14 July 2022

Sec. Smart Grids

Volume 10 - 2022 | https://doi.org/10.3389/fenrg.2022.902776

Robust Load Frequency Control of Hybrid Solar Power Systems Using Optimization Techniques

  • Department of Electrical and Computer Engineering, King Abdulaziz University, Jeddah, Saudi Arabia

Abstract

It is necessary to predict solar photovoltaic (PV) output and load profile to guarantee the security, stability, and reliability of hybrid solar power systems. Severe frequency fluctuations in hybrid solar systems are expected due to the intermittent nature of the solar photovoltaic (PV) output and the unexpected variation in load. This paper proposes designing a PID controller along with the integration of a battery energy storage system (BESS) and plug-in hybrid electric vehicle (PHEV) for frequency damping in the hybrid solar power system. The solar PV output is predicted with high accuracy using artificial neural networks (ANN) given that solar irradiance and cell temperature are inputs to the model. The variation in load is also forecasted considering the factors affecting the load using ANN. Optimum values of the PID controller have been found using genetic algorithm, particle swarm optimization, artificial bee colony, and firefly algorithm considering integral absolute error (IAE), integral square error (ISE), and integral time absolute error (ITAE) objective functions. IAE, ISE, ITAE, Rise time, settling time, peak overshoot and maximum frequency deviation have been measured for comparison and effectiveness. The transient behavior has been further improved by utilizing the power from BESS/PHEV to the power system. The results demonstrate the efficacy of the suggested design for frequency control using the genetic algorithm method along with ISE objective function compared with those obtained from the conventional, particle swarm optimization, artificial bee colony, and firefly algorithm techniques.

Introduction

Solar photovoltaic (PV) systems are intermittent in nature. The variation in generated power and load demand is common in solar PV connected power systems. Under normal operating situations, microgrids are powered by the integration of both solar PV systems and the power grid. However, when the output from the solar PV system is lower than expected, and the load demand is high, the power grid may not be able to support enough power because of its slow dynamic response. This will cause the power system to be unstable which leads to severe frequency fluctuations in the system ().

Hybrid solar power systems are a source for generating electrical power that use inverters coupled with batteries to store energy later use. In this way, hybrid solar systems will be able to use the stored energy during blackouts and compensate for renewable energy sources output fluctuation (). Integrating microgrids with renewable energy systems is considered an adequate solution to meet the increasing demand of electric power. Such microgrids can utilize plug-in hybrid electric vehicles and battery energy storage systems to maintain the balance between generated power and the load demand. Due to the intermittent nature of renewable energy sources, the change in generation power and load demand would lead the system to instability causing severe frequency fluctuations. Mathematical models were developed to predict the behavior of microgrids in order to dampen frequency fluctuations. Power systems with a high penetration of renewable energy sources were integrated with battery energy storage system (BESS) and plug-in hybrid electric vehicles to maintain the frequency deviations within limits ().

Mellit () uses artificial neural network methods for forecasting the electric power produced by a 20 grid-connected PV system placed on a rooftop in the municipality of Trieste, Italy. According to Huang (), improved the solar photovoltaic (PV) output by using the solar irradiance and temperature along solar zenith angle and solar azimuth angle with the help of artificial neural network methods. A data-driven ensemble approach was used by Al-Dahidi in () to predict the day ahead solar photovoltaic (PV) output (230 capacity) placed on the top of the Applied Science Private University, Amman, Jordan to overcome the intermittent nature of solar energy system.

studied the PID controller calibrating approaches for load frequency control of Power systems. The prime challenge in LFC was to design a PID controller that could maintain the frequency fluctuations rigorously within the defined limits. Fractional Order PID controller was tuned by for load frequency control for three area power system with the help of imperialist competitive algorithm (ICA). The response of the interconnected power system towards disruption was smooth and less vibratory as a result of using the proposed controller. Mosaad () presented a methodological design for an adaptive PID load frequency using artificial neural networks (ANN) and adaptive Neuro-Fuzzy Inference systems (ANFIS). The PID Controller was able to maintain the effective performance at any load point in the power system with the superiority of ANFIS over ANN with respect to integral absolute error (IAE), integral square of errors. According to , there are possible chances of frequency fluctuations after losing balance between power generated and load demand. Battery storage was incorporated because of its faster response to overcome load frequency deviations using recurrent neural networks.

The desired frequency of the power system deviates due to disturbances in the electric power generation and load fluctuations. Load frequency control has been considered the main issue for power system operation and control for effective and reliable power supply (). PI, PID and fuzzy controllers have been used to stabilize the power system when it is integrated with wind power. introduce the ANN-based PID controller to sufficiently damp frequency fluctuations produced as a result of load changes. An artificial intelligent controller has been designed using the backpropagation algorithm to reduce the area control error of a two-area hydrothermal power system for automatic generation control purposes in ().

Meta-heuristic optimization techniques have been employed to get the optimal gains of PID controller for LFC problem. According to (), genetic algorithm based PID controller was tuned for autonomous hybrid generation system. Fractional order PID controller was tuned using genetic algorithm for hybrid power system integrated with renewable energy sources (). Multi-objective artificial bee colony technique was adopted to tune the PID controller parameters with the help of integral time multiplied absolute error (ITAE) and integral of time weighted squared error (ITSE) for two area thermal power system (). Particle swarm optimization has been utilized to tune the fuzzy logic controller for the frequency restoration of multi-area power system (). PSO is used for obtaining the optimal PID parameters for two area system with integral square error as cost function and showing best convergence profile (). PID optimal parameters obtained for two and three area power system using firefly algorithm considering the generation rate constraint and governor dead band ().

Energy storage systems are also used in previous studies for frequency regulation. The role of battery energy storage system has been presented for regulating the frequency in interconnected power system. The results are compared to the conventional PID controller that describe BESS as efficient for LFC purpose (). Peak frequency deviation and settling time of two area power system is reduced using the battery energy storage system and choosing the proper integral gain (). Plug-in electric vehicles power is utilized for frequency stability of the islanded power system with the help of controllable power rate strategy (). The energy from electric vehicles that are charged during daytime can be utilized in the night-time for small residential power system ().

The purpose of this study to address the limitations in the previous studies and fulfill the gap. Many load frequency control techniques have been applied to power system for its stability purposes, i.e., robust control, decentralized aspect, linear quadratic, pole shifting, and variable structure. But there have been some drawbacks of these techniques which decrease their execution efficiency. To overcome such problem, artificial intelligence techniques, fuzzy logic (FL) and neural network (NN) have been adopted and applied to solve the non-linear models of the power system. These techniques have proved their effectiveness and efficiency to solve the problem. But they still have some disadvantages. For example, it is hard to train the neural network again and again to set the number of neurons and other parameters to get the best output response. Similarly, fuzzy logic requires hard work to get the influential signal response from the power system (). Optimization techniques such as genetic algorithm (GA) () (; ; ; ; ; ; ; ; ), particle swarm optimization (PSO) (; ) (; ; ; ; ; ; ; ; ; ; ) () ()- (), firefly algorithm (FA) (; ; ; ; ; ; ), and artificial bee colony (; ; ; ; ; ; ; ), have been used and found to be effective for load frequency control problem. But proportional, integral and derivative gain range of PID controller during optimization problem has been narrowed. Renewable energy source, photovoltaic and wind turbine, transfer function model has been adopted in the previous studies. In the same way, load frequency control has also been done in previous studies using the energy storage systems but without efficient PID controller that ultimately require large energy storage capacity to damp the frequency ()- (). Taking into account the mentioned limitations, the current study discusses and solve the problem by taking real-time data to train the network for solar photovoltaic prediction, load deviation forecasting, widening proportional, integral, and derivative controller gain range, and integrating the BESS/PHEV with efficient tuned PID controller.

The major contributions of this paper are threefold and shown below.

  • 1) Solar photovoltaic generated power and load deviation have been forecasted from the affecting factors using artificial neural networks.

  • 2) The forecasted solar photovoltaic power and load deviation have been utilized by the power system model to find the optimal PID controller parameters using genetic algorithm, particle swarm optimization, artificial bee colony and firefly algorithm optimization techniques considering the integral absolute error, integral square error, and integral time absolute error objective functions.

  • 3) BESS/PHEV model, with gain levels 10 and 50 respectively, has been integrated to the power system model to further dampen the frequency deviation.

Power System Model

Power system control concentrate on steady state operation. This study presents how to deal with active power to keep the power system in steady state. The major objective of the control strategy is to provide qualitative and reliable power to customers within an interconnected system. Changes in the active power disturb the frequency of the system. Therefore, a control strategy is designed to regulate load frequency control using control loops. Two common approaches, transfer function and state variable, are adopted to transform the power system model into a mathematical model by making some proper assumptions (; ; ; ; ; ). A single area power system model is shown in Figure 1.

FIGURE 1

Generator Model

The generator equation has been extracted from the swing equation as shown in Eq. 1.

Load Model

The power system contains resistive and inductive load that is frequency independent and dependent respectively (). Therefore, the net change in load power can be described as the combination of frequency sensitive and frequency non-sensitive load changes shown in Eq. 2.Where is the frequency independent load change and denotes the frequency-sensitive load change and represent the ratio of percent change in load to the percent change in frequency. So, the relationship for variation in load with respect to frequency fluctuation can be expressed as,

Many factors affect the electric load (

). These factors could be short term, middle term, and long-term influence factors. Some important factors used for load forecasting are categorized as follows.

  • • Meteorological factors involving temperature, wind speed, humidity, surface pressure, and precipitation factor are responsible for load forecasting.

  • • Temporal or calendar factors involving hour of the day and day of the week.

  • • Random factors such as sports activities.

Frequency dependent load change has been observed from the Atlas Power Plant. The Atlas power plant is a coal fired thermal power plant located at the geographic location (Longitude:31.5204, Latitude:74.3587). The historical load demand data has been extracted from the power data reference book produced by the National Transmission and Dispatch Company Limited (NTDCL) (). Historical daily load data from 2011 to 2016 has been taken as input along with weather and seasonal data. Seasonal and weather data consists of wind speed, surface pressure, temperature, precipitation factor, and humidity. Average load has been calculated for 365 days from load profile data between 2011 and 2016 (). Frequency dependent change in load has been extracted from the average load for consecutive 365 days as shown in Figure 2. This change in load power has been taken as an output ().

FIGURE 2

The weather data has been first normalized by taking the maximum value as base from each factor affecting the load and dividing the entirety of the weather data with base value within the range (0–1). Similarly, historical data has also been normalized by dividing with the base value 219 MW and converting the load data within the range of (0–1). After normalizing the input and target data, NNTOOL is used to train the network. A multi-layer feed forward neural network was trained using wind speed, surface pressure, temperature, precipitation factor, humidity, and historical load data as input and deviation in load power data as target as shown in Figure 1 of the load deviation model (). The number of neurons has been chosen as 10, Tan sigmoid activation function, and Levenberg Marquardt as learning algorithm as shown in Table 1. The trained model depicts its effectiveness showing the coefficient of regression as 1.0 and mean squared error (MSE) as at 525 epochs. After training model, the network model was deployed in Simulink, where the neural network can successfully forecast the deviation in load power.

TABLE 1

Network typeMulti-layer feedforward neural network
Inputs6 inputs- wind speed, surface pressure, temperature, precipitation factor, humidity, historical load 2011–2016
Outputs1 output-deviation in load
Number of layersInput, hidden, output
Number of hidden neurons10
Activation functionTan-sigmoid, linear
Learning algorithmLevenberg-Marquardt

Summary of artificial neural network design and architecture for forecasting load deviation.

Prime Mover Model

It is the origin of mechanical power whose energy is obtained from burning coal or gas, or nuclear fission. The transfer function of the turbine can be represented as the ratio of the change in mechanical output power to the change in steam valve position as represented by Eq. 4.Where is the turbine time constant.

Governor Model

The speed governor model can be expressed as Eq. 5.The speed governor operates as a comparator and its output can be written as difference between the reference set power and power , where R indicates the speed regulation. The relation between governor input and valve opening can be expressed by Eq. 6 ().Where is taken as the governor time constant in seconds.

Solar Photovoltaic Model

It is fact that photovoltaic output power is intermittent in nature. It depends on two main factors: solar irradiance and temperature (). Solar photovoltaic power is directly proportional to solar irradiance while inversely proportional to temperature (). Solar photovoltaic model of 5 MW is considered at location (Longitude:31.5204, Latitude:74.3587) near the thermal power plant. Solar PV system is integrated to the thermal power plant. The solar power plant is built by LONGI solar modules of power 540W. LONGI solar module characteristics are shown in Table 2. The solar PV module output (W) was considered at different solar irradiance (ranges 100–1100 ) and temperature (0°C–55°C). These changes were applied to the whole 5 MW power plant to see the effect of intermittency on the power plant. The solar PV generated power 5 MW has been normalized using the base power of 219 MW according to the thermal power plant. A multi-layer feed forward neural network was trained using solar irradiance and temperature as input data and solar power as target data (; ). The number of neurons is chosen to be 10, Tan sigmoid activation function, and Levenberg Marquardt as learning algorithm as shown in Table 3. The trained model depicts its effectiveness showing the coefficient of regression as 1.0 and mean squared error as at epochs 22. After successfully training the model, the network model was deployed using Simulink, where neural network can successfully forecast the solar power using the input data taken from the National Aeronautics and Space Administration (NASA). The peak value of solar irradiance and its corresponding temperature at a specific time of the day is considered to forecast the solar photovoltaic output power for that day.

TABLE 2

Module typeLR5-72HPH-540M
Testing conditionSTC:AM1.5, 1000 , 25 NOCT:AM1.5, 800 , 20
Maximum power 540403.3
Open circuit voltage 49.5046.41
Short circuit current 13.8511.20
Voltage at maximum power 41.6538.78
Current at maximum power 12.9710.40
Module efficiency (%)21.1

Electrical characteristics of LONGI Solar 540W module.

TABLE 3

Network typeMulti-layer feedforward neural network
Inputs2 inputs- solar irradiance, temperature
Outputs1 output-module power
Number of layersInput, hidden, output
Number of hidden neurons10
Activation functionTan-sigmoid, linear
Learning algorithmLevenberg-Marquardt

Summary of artificial neural network design and architecture for forecasting solar photovoltaic power.

Battery Energy Storage System/Plug-in Hybrid Electric Vehicle Model

Energy storage systems are a great source of frequency damping in the power system. As discussed in the introduction section, power from battery energy storage systems and plug-in hybrid electric vehicle play an important role for maintaining the frequency. Energy from PHEV can be utilized for small residential systems. The circle of willingness of customers for utilizing PHEV can be enhanced. Figure 1 consists of the control system model for a battery energy storage system or plug-in hybrid electric vehicle integrated to the power system. The term depicts the mechanism of auto balancing charging loop. This close-loop control will maintain the energy of the BESS/PHEV near the initial state, which is considered 50% charged. is the feedback gain from the frequency fluctuation in the power grid. represent the first order conversion delay from DC to AC of the BESS/PHEV (; ).

PID Controller Model for Optimization Problem

The PID Controller model has been adopted in the power system model as shown in Figure 1. The transfer function of the controller is as follows.

The control signal for maintaining the system frequency is given by Eq. 8 ().where is the area control error of the power system and is the input signal to the governor for controlling the valve output according to load demand of the power system.

Area Control Error (ACE) of studied single-area power system can be written by Eq. 9. is called bias factor and denoted with , and represents the frequency deviation.

Multi-objective function is used in () for load frequency control using the artificial bee colony optimization approach considering ITAE and ITSE. It has been observed the performance of the PID controller can be compromised, i.e., Either frequency deviation or settling time is compromised on the cost of each other. Therefore, it is recommended to use objective functions separately. But, the most suitable objective function should be identified based on the performance indices and convergence profile.

Three objective functions are used to check the performance indices of the PID controller, which are the integral absolute error (IAE), integral square error (ISE), and integral time absolute Error (ITAE) and given by Eqs 10–12.Where is area control error (ACE).

Frequency deviation permissible limit is calculated using the following steady-state equation ().

The maximum steady-state frequency deviation in Hz has found to be −0.7932 Hz for maximum positive load change and 0.5468 Hz for maximum negative load change for studied model. So, the frequency deviation permissible limit range should be (−0.7932, 0.5468) in Hz.

Proposed Approaches for Tuning the PID Controller

Particle Swarm Optimization

Particle swarm optimization technique was developed by Kennedy and Eberhart in 1995. They used nature inspired optimization algorithms in their technique. In this technique, particles are flown through the search space and update the position of the particle at time step of . The expression for the velocity updates is given by the Eq. 14.

This technique is adopted to find the optimal values of the PID controller parameters. Initially, the PID controller proportional gain , integral gain , and derivative gain optimal values have been searched in the range (−10, 1000). The searched controller gains produce insufficient closed-loop stability of the power system in this range. Therefore, the gain parameters range has been narrowed iteratively until closed-loop stability condition of power system model is achieved within the range (−10, 150). PID controller produces best results and minimizes the error significantly within the selected range (−10, 150) as compared to very short range (0, 5) (; ) (; ; ; ; ; ; ; ; ; ; ) () ()- (). Table 4 depicts the PSO operators used to find out the optimal PID controller parameters. The best performing PID controller parameters are mentioned in Table 5.

TABLE 4

OptionNumber/Type
Genetic algorithm
Number of variables3
Limit(−10, 150)
Population sizeDefault, 50
Creation functionUniform
Fitness scaling functionRank
Selection functionTournament
Crossover fractionDefault, 0.8
Mutation functionAdaptive feasible
Crossover functionArithmetic
Fitness functionIAE, ISE, ITAE
Stopping criteriaDefault, 100*Number of variables
Particle swarm optimization
No. of variables3
No. of particles30
No. of iterations50
Objective functionsIAE, ISE, ITAE
Limit(−10, 150)
2
0.9
0.4
Artificial bee colony
No. of populations30
No. of iterations60
Number of variables3
Objective functionsIAE, ISE, ITAE
Limit(−10, 150)
a (acceleration coefficient)1
Number of onlooker bees30
Firefly algorithm
No. of populations30
No. of iterations50
Number of variables3
Objective functionsIAE, ISE, ITAE
Limit(−10, 150)
1
1
0.2
0.98

Meta-heuristic techniques operators used to find out optimal PID gain parameters.

TABLE 5

PID gain parametersIAEISEITAE
Conventional0.81770.18930.39080.81770.18930.39080.81770.18930.3908
GA22.210567.344614.112218.79140.5523.1820.907133.76249.5543
PSO22.172266.980114.108121.446018.6322.807051.491514.4405
ABC22.195962.805013.482716.907138.07920.72922.111465.772713.6213
FA5.89338.922121.554810.0141.091514.7154.97956.45862.0571

PID controller gain parameters using IAE, ISE, and ITAE objective function.

Genetic Algorithm

Genetic Algorithm is an optimization technique based on heredity and evolution. This technique was suggested by John Holland in 1960 to search for the best solution of complex problems. It is an iterative process which maintains the constant population size of desired solution. GA initially begins with a randomly selected population of function input that is represented as a bit of strings. The Population is evaluated in each iterative step called generations to give a new population of the desired solution. That means GA uses the current population to produce a new population such that the new population is better on average if compared to the previous population. Best elements are used from the current population to form the best population. The successful process will produce a population better than the old one. Three steps selection, mating, and mutation are used to produce new population from the old population. The basic genetic algorithm cycle is shown in Figure 3. New generations are produced iteratively with the repetition of these three steps. The process continues until the stopping criteria is reached like the maximum number of iterations is achieved or no improvements () (; ; ; ; ; ; ; ; ).

FIGURE 3

This technique is adopted to find the optimal values of the PID controller parameters. The PID controller proportional gain , integral gain , and derivative gain optimal values have been searched in the range (−10, 150). Different parameters and functions for applying GA technique are used in OPTIMTOOL toolbox of MATLAB and given in Table 4. The optimal parameters obtained are depicted in Table 5.

Firefly Algorithm

The flashing light of fireflies is of great importance to study their pattern and rhythmic movement. Studies discussed the main aim of pattern of flashes is either to communicate with the mating partners, to attract the potential prey or to signal other fireflies for potential safety warning. This signal consists of rhythmic flash, rate of flash and amount of flashing. This behavior of flashing light from fireflies can be used to solve many optimization problems. This algorithm has been developed by Xin-She Yang in 2008.

The Euclidian distance between two fireflies and with respect to their positions and can be found with the flowing Eq. 15Where denotes the element of the geographical coordinates. The attractiveness between fireflies can be denoted with the following attractiveness parameter. is called the coefficient of absorption. It is used to control the flashing light concentration. Fireflies’ positions and movement can be described by the following equations. indicates the instantaneous position of the firefly whereas represents the random behavior of firefly when the firefly can’t see another brighter firefly (; ; ; ; ; ; ).

The PID controller proportional gain , integral gain , and derivative gain optimal values have been searched in the range (−10, 150) using firefly algorithm. Firefly algorithm operators used to find out the PID gain values are given in Table 4 and obtained optimal PID parameters in Table 5.

Artificial Bee Colony

Artificial bee colony is meta-heuristic optimization algorithm developed by Karaboğa in 2005. Honey bees foraging behavior has been utilized in this technique searching for nectar at different positions. Usually, three kind of honey bees are considered to be responsible for finding nectar positions: employed, onlooker and scouts. This algorithm adopts some variables and steps to proceed for finding optimal solution. The variables of ABC algorithm include population of honey bees (SN), maximum cycle number (MCN), objective functions and some functions to test the probability, fitness and optimization of the objective functions. The steps in the ABC algorithm includes initialization, employed, onlooker and scout bees. The initialization step starts by generating the number of solutions equal to total number of bees (SN). Each solution is obtained within the upper and lower limit of the selected decision variable with the following equation.Where denotes the total number of decision variables used in the objective function. The next step starts by generating the new food source by employed and onlooker bees by updating their positions as given by,Where is random number within range (−1, 1) and used to change the position around . The last step includes the probability of finding the nectar by onlooker bees, which is given by following expression,Where is the fitness value corresponding to nectar position . The fitness value is checked for the next iteration, and if it is better than previous fitness value the latest one is stored and the algorithm continues until the stopping criteria is reached (; ; ; ; ; ; ; ). Using the artificial bee colony, optimized PID controller parameters are found as shown in Table 5.

Conventional PID Controller

The Simulink model of the single-area power system containing PID controller is tuned using the classical method to enhance the performance and robustness. PID Tuner toolbox is launched from the Simulink model that automatically computes the linear power system model to find gain parameters corresponding to initial controller conditions. The obtained PID controller gains can robustly stabilize the power system. The performance of controller could be checked using frequency deviation performance indicators, i.e., settling time, zero steady state error. If the desired response is not achieved by auto tunning, PID controller gains can be adjusted manually by varying the response time (seconds) and transient behavior. The PID controller gains are updated to the Simulink model after meeting the desired requirement ()- (). The compensator formula for PID controller is given as,

In the studied model, best performing PID controller gains have been found to be 0.8177, 0.1893, and 0.3908 within the range (0, 1) corresponding to initial controller conditions. PID controller has been tuned with response time value of 13.63 s and transient behavior value of 0.6 between aggressiveness and robustness. Filter coefficient value has been chosen as 100 for this case.

Results and Discussion

Three objective functions IAE, ISE, and ITAE have been applied to the genetic algorithm, particle swarm optimization, artificial bee colony, and firefly algorithm to find out the optimal proportional gain , integral gain , and derivative gain . Table 5 illustrates the PID gain parameters for GA, PSO, ABC, and FA using three objective functions. GA, PSO, and ABC techniques have found to be competitive and efficient than FA and conventional methods. The convergence characteristics of GA, PSO, and ABC are shown in Figure 4. The graph represents the objective function vs. number of iterations. The graph depicts almost two iterations are required to reach the best fitness 0.0219, 0.00034, and 0.0166 using IAE, ISE, and ITAE objective function respectively for genetic algorithm. It can also be seen that almost 40 iterations are required to reach the best fitness 2.20, 0.0344, and 44.0306 using IAE, ISE, and ITAE cost function respectively for particle swarm optimization. No significant change has been observed after 40 iterations. The fitness value has been found 0.73, 0.0115, and 14.55 using IAE, ISE, and ITAE objective function respectively in less than 20 iterations for artificial bee colony. Convergence properties describe the effective genetic algorithm that requires very few iterations to reach its best fitness value.

FIGURE 4

The results of the studied model have been categorized into three cases.

  • • Frequency damping using integral absolute error (IAE) objective function, with and without BESS/PHEV support.

  • • Frequency damping using integral square error (ISE) objective function, with and without BESS/PHEV support.

  • • Frequency damping using integral time absolute error (ITAE) objective function, with and without BESS/PHEV support.

Several parameters are used to check the performance indices of hybrid power system for frequency deviation including integral absolute error, integral square error, and integral time absolute error. Rise time , settling time , peak overshoot , and maximum frequency deviation are also used to check the transient response of the power system. These performance indices are also checked by integrating power from BESS/PHEV for different capacities to see the frequency deviation suppression.

Frequency of the power system deviates either positive or negative depending upon the increase or decrease in load and photovoltaic power. The mismatch between load demand and photovoltaic power give rise to changes in frequency. Five key days, 1st January ( = 0.3394 pu), 30th March ( = 0.2374 pu), 13th December ( = 0.1575 pu), 24th May ( = −0.2275 pu), and 28th September ( = −0.1499 pu), showing significant load changes have been selected to check the conventional PID, GA optimized PID, PSO optimized PID, ABC optimized PID, and FA optimized PID responses against frequency deviation. Figures 5–7 depicts the frequency deviation responses of the power system in the same order as mentioned above from top to bottom.

FIGURE 5

Frequency deviation responses of the GA optimized PID, PSO optimized PID, ABC optimized PID, and FA optimized PID using the IAE objective function is shown in Table 6 and Figure 5. Their performance indices show that GA optimized PID, PSO optimized PID, and ABC optimized PID produce approximately same results as their proportional, integral, and derivative values have been found to be same. IAE, ISE, ITAE and are higher for higher mismatches between load and photovoltaic generated power and lower for lower mismatches, i.e., The values of IAE, ISE, ISE and are 0.02199, 0.0003987, 0.01991, and −0.1148 for GA, PSO, and ABC considering the 0.3394 pu load deviation and 0.01351 pu solar PV power for 1st January loading situation in single-area power system. These values of IAE, ISE, ITAE and are lowered to 0.01469, 0.000178, 0.01331, and −0.07668 for GA, PSO, and ABC for 0.2374 pu load deviation and 0.0147 pu solar PV power on 13th December loading conditions. FA optimized PID produces relatively poor result as the IAE, ISE, and ITAE values 0.1106, 0.001322, 0.61 for 1st January loading conditions indicate much higher than GA, PSO and ABC. The values of rise time , settling , and peak overshoot are same for all loading conditions. For example, the values of , settling , and peak overshoot are 0.0673, 2.87, and 78.7 for GA and PSO, 0.069, 2.77, and 78 for ABC on 1st January.

TABLE 6

Without BESS/PHEV support, KB = 10 and KB = 50
DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.7180.38936.5020.9127.978.15−0.9184
 30th Mar0.23740.019691.1480.17384.3450.9127.978.15−0.6136
 13th Dec0.15750.011920.76740.077722.9050.9127.978.15−0.4103
 24th May−0.22750.01851.2970.22194.9090.9127.978.150.6934
 28th Sep−0.14990.014510.86660.09913.2810.9127.978.150.4634
GA-optimized-PID
 1st Jan0.33940.013510.021990.00039870.019910.06732.8778.7−0.1148
 30th Mar0.23740.019690.014690.0001780.013310.06732.8778.7−0.07668
 13th Dec0.15750.011920.0098230.00007950.0088970.06732.8778.7−0.05132
 24th May−0.22750.01850.01660.00022720.015030.06732.8778.70.08664
 28th Sep−0.14990.014510.011090.00010150.010050.06732.8778.70.05792
PSO-optimized-PID
 1st Jan0.33940.013510.021990.00039870.019940.06732.8678.6−0.1148
 30th Mar0.23740.019690.014690.0001780.013320.06732.8678.6−0.07668
 13th Dec0.15750.011920.0098240.00007950.008910.06732.8678.6−0.05132
 24th May−0.22750.01850.01660.00022720.015060.06732.8678.60.08664
 28th Sep−0.14990.014510.011090.00010150.010060.06732.8678.60.05792
ABC-optimized-PID
 1st Jan0.33940.013510.0220.00041330.019130.0692.7778−0.1171
 30th Mar0.23740.019690.01470.0001840.012780.0692.7778−0.0783
 13th Dec0.15750.011920.0098310.0000820.008540.0692.7778−0.0523
 24th May−0.22750.01850.016610.0002350.014440.0692.77780.0885
 28th Sep−0.14990.014510.01110.0001050.009650.0692.77780.05918
FA-optimized-PID
 1st Jan0.33940.013510.11060.0013220.610.05483.8476.9−0.09573
 30th Mar0.23740.019690.07390.000590.40760.05483.8476.9−0.0641
 13th Dec0.15750.011920.04930.0002640.27250.05483.8476.9−0.0428
 24th May−0.22750.01850.083520.0007530.46050.05483.8476.90.0722
 28th Sep−0.14990.014510.058820.0003360.30780.05483.8476.90.0484
With BESS/PHEV support, and
 DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.721.2990.23560.0895610.2411.277.93—693—10.1—−0.6085−0.2554
 30th Mar0.23740.019691.1490.86810.10520.039986.8397.537.93—693—10.1—−0.4066−0.1705
 13th Dec0.15750.011920.78850.58050.047030.017884.5735.0357.93—693—10.1—−0.2719−0.1140
 24th May−0.22750.01851.2990.98090.13430.051057.7288.5097.93—693—10.1—0.45940.1928
 28th Sep−0.14990.014510.86780.65550.059970.02285.1645.6867.93—693—10.1—0.30700.1289
GA-Optimized-PID
 1st Jan0.33940.013510.018730.013530.00032030.00019020.016960.013990.06830.07382.31.757041−0.1103−0.09469
 30th Mar0.23740.019690.012520.009040.0001430.00008480.011330.0093440.06830.07382.31.757041−0.07368−0.06326
 13th Dec0.15750.011920.008370.0060450.00006390.00003790.0075760.0062480.06830.07382.31.757041−0.04928−0.04229
 24th May−0.22750.01850.014140.010220.00018260.00010830.01280.010560.06830.07382.31.7570410.083230.07146
 28th Sep−0.14990.014510.0094520.0068270.00008150.00004830.0085550.0070560.06830.07382.31.7570410.055610.04779
PSO-Optimized-PID
 1st Jan0.33940.013510.018740.013550.00032050.00019020.016990.014020.06830.07382.31.7569.940.9−0.1103−0.09469
 30th Mar0.23740.019690.012520.0090520.00014310.00008490.011350.009370.06830.07382.31.7569.940.9−0.07368−0.06326
 13th Dec0.15750.011920.0083730.0060530.00006390.00003790.0075890.0062650.06830.07382.31.7569.940.9−0.04928−0.04229
 24th May−0.22750.01850.014150.010230.00018270.00010840.012820.010590.06830.07382.31.7569.940.90.083230.07146
 28th Sep−0.14990.014510.0094550.0068350.00008150.00004840.008570.0070750.06830.07382.31.7569.940.90.055610.04779
ABC-optimized-PID
 1st Jan0.33940.013510.018720.013460.00032360.0001970.016190.013240.070.07562.31.7769.140−0.1122−0.0959
 30th Mar0.23740.019690.012510.008990.00014850.0000880.010820.008840.070.07562.31.7769.140−0.0749−0.064
 13th Dec0.15750.011920.008360.006110.00006640.0000390.007230.005910.070.07562.31.7769.140−0.0501−0.0429
 24th May−0.22750.01850.014130.010170.0001890.0001120.012220.009990.070.07562.31.7769.1400.08470.0724
 28th Sep−0.14990.014510.009440.006790.0000840.0000500.0081680.006670.070.07562.31.7769.1400.05660.0484
FA-optimized-PID
 1st Jan0.33940.013510.10410.08440.0011920.000870.55580.3990.05530.0592.758.9869.946.1−0.0926−0.0823
 30th Mar0.23740.019690.06950.056420.0005320.0003880.37140.2660.05530.0592.758.9869.946.1−0.06190.05504
 13th Dec0.15750.011920.04650.03770.0002370.0001730.24830.17820.05530.0592.758.9869.946.1−0.0414−0.0366
 24th May−0.22750.01850.07860.06370.0006790.0004960.41960.30120.05530.0592.758.9869.946.10.0700.06194
 28th Sep−0.14990.014510.052540.04260.0003030.0002210.28040.20130.05530.0592.758.9869.946.10.04680.0415

Frequency deviation parameter indices using integral absolute error objective function.

Considering the effect of BESS/PHEV with its different capacities, BESS/PHEV of 10 and 50 are used to check frequency deviation response. Results indicate that settling , peak overshoot are reduced to 2.3, 70 for GA, PSO, and ABC for of 10. These values are further reduced to 1.75 and 40 for of 50. The values of IAE, ISE, ISE and are reduced to 0.01873, 0.0003203, 0.01696, and −0.1103 for GA, PSO, ABC for of 10, and 0.01353, 0.0001902, 0.01399, and −0.0946 for 50 for 1st January loading conditions. Rise time remains unchanged in any case.

Frequency deviation responses of the GA optimized PID, PSO optimized PID, ABC optimized PID, and FA optimized PID using the ISE objective function is shown in Table 7 and Figure 6. Their performance indices show that GA optimized PID, PSO optimized PID, and ABC optimized PID produce efficient and competitive results. IAE, ISE, ITAE and are higher for higher mismatches between load and photovoltaic generated power and lower for lower mismatches, i.e., The values of IAE, ISE, ISE and are 0.02719, 0.0003375, 0.04641, and −0.09178 for GA, 0.02339, 0.0003489, 0.02866, and −0.1014 for PSO, 0.02629, 0.000327, 0.04707, and −0.0966 for ABC considering the 0.3394 pu load deviation and 0.01351 pu solar PV power on 1st January in a single-area power system. These values of IAE, ISE, ITAE and are lowered to 0.01817, 0.0001507, 0.03101, and −0.06102 for GA, 0.01563, 0.0001557, 0.01915, and −0.06781 for PSO, 0.01758, 0.000146, 0.03145 and −0.0644 for ABC for 0.2374 pu load deviation and 0.0147 pu solar PV power on 13th December loading conditions. The values of rise time , settling , and peak overshoot remain same for different loading conditions. For example, the values of rise time , settling , and peak overshoot are 0.0518, 4.51, and 87 for GA, 0.0581, 3.35, and 81.6 for PSO, 0.0551, 3.37, and 84 for ABC on 1st January loading conditions. FA optimized PID produces relatively poor result as the IAE, ISE, and ITAE values are much higher than GA, PSO and ABC as can be seen their values are 0.275, 0.005127, and 1.896.

TABLE 7

Without BESS/PHEV support
DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.7180.38936.5020.9127.978.15−0.9184
 30th Mar0.23740.019691.1480.17384.3450.9127.978.15−0.6136
 13th Dec0.15750.011920.76740.077722.9050.9127.978.15−0.4103
 24th May−0.22750.01851.2970.22194.9090.9127.978.150.6934
 28th Sep−0.14990.014510.86660.09913.2810.9127.978.150.4634
GA-Optimized-PID
 1st Jan0.33940.013510.027190.00033750.046410.05184.5187−0.09178
 30th Mar0.23740.019690.018170.00015070.031010.05184.5187−0.06102
 13th Dec0.15750.011920.012150.00006730.020740.05184.5187−0.04110
 24th May−0.22750.01850.020530.00019240.035040.05184.51870.06937
 28th Sep−0.14990.014510.013720.00008590.023420.05184.51870.04637
PSO-Optimized-PID
 1st Jan0.33940.013510.023390.00034890.028660.05813.3581.6−0.1014
 30th Mar0.23740.019690.015630.00015570.019150.05813.3581.6−0.06781
 13th Dec0.15750.011920.010450.000069640.01280.05813.3581.6−0.04532
 24th May−0.22750.01850.017660.00019880.021640.05813.3581.60.0766
 28th Sep−0.14990.014510.01180.00008880.014460.05813.3581.60.05118
ABC-optimized-PID
 1st Jan0.33940.013510.026290.0003270.047070.05513.3784−0.0966
 30th Mar0.23740.019690.017580.0001460.031450.05513.3784−0.0644
 13th Dec0.15750.011920.011740.0000650.021030.05513.3784−0.0431
 24th May−0.22750.01850.019850.00018680.035540.05513.37840.0731
 28th Sep−0.14990.014510.013260.00008340.023750.05513.37840.0487
FA-optimized-PID
 1st Jan0.33940.013510.27560.0051271.8960.06772.2767.9−0.1146
 30th Mar0.23740.019690.18410.002281.2670.06772.2767.9−0.0764
 13th Dec0.15750.011920.12310.001020.84690.06772.2767.9−0.0511
 24th May−0.22750.01850.20810.0029221.4320.06772.2767.90.0865
 28th Sep−0.14990.014510.1390.0013050.95640.06772.2767.90.0577
With BESS/PHEV support, and
DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.721.2990.23560.0895610.2411.277.93—693—10.1—−0.6085−0.2554
 30th Mar0.23740.019691.1490.86810.10520.039986.8397.537.93—693—10.1—−0.4066−0.1705
 13th Dec0.15750.011920.78850.58050.047030.017884.5735.0357.93—693—10.1—−0.2719−0.1140
 24th May−0.22750.01851.2990.98090.13430.051057.7288.5097.93—693—10.1—0.45940.1928
 28th Sep−0.14990.014510.86780.65550.059970.02285.1645.6867.93—693—10.1—0.30700.1289
GA-optimized-PID
 1st Jan0.33940.013510.022580.016860.00025320.00013980.039970.034580.05220.05513.422.778055.9−0.08901−0.07942
 30th Mar0.23740.019690.015090.011260.0001130.00006240.028710.023110.05220.05513.422.778055.9−0.05950−0.05293
 13th Dec0.15750.011920.010090.0075310.00005050.00002790.017880.015450.05220.05513.422.778055.9−0.03987−0.03547
 24th May−0.22750.01850.017050.012730.00014430.00007960.030180.026110.05220.05513.422.778055.90.067380.05988
 28th Sep−0.14990.014510.011390.0085050.00006440.00003550.020170.017450.05220.05513.422.778055.90.045060.04005
PSO-optimized-PID
 1st Jan0.33940.013510.019970.014750.00027660.00016460.024970.020490.05880.06282.642.27447.8−0.09794−0.08589
 30th Mar0.23740.019690.013340.0098560.00012350.000073470.016680.013690.05880.06282.642.27447.8−0.06544−0.05740
 13th Dec0.15750.011920.008920.006590.00005520.000032850.011160.009150.05880.06282.642.27447.8−0.04378−0.03838
 24th May−0.22750.01850.015080.011140.00015770.00009380.018850.01540.05880.06282.642.27447.80.074010.06494
 28th Sep−0.14990.014510.010080.0074420.00007040.00004190.01260.010340.05880.06282.642.27447.80.049440.04341
ABC-optimized-PID
 1st Jan0.33940.013510.023230.018350.000260.0001550.04370.039270.0550.0592.92.776.751.6−0.09344−0.08275
 30th Mar0.23740.019690.015520.012260.000110.0000690.02920.026240.0550.0592.92.776.751.6−0.0624−0.0554
 13th Dec0.15750.011920.010380.008190.0005190.0000310.01950.017550.0550.0592.92.776.751.6−0.041−0.037
 24th May−0.22750.01850.017540.013850.0001480.0000880.032990.029650.0550.0592.92.776.751.60.07060.0625
 28th Sep−0.14990.014510.011720.0092570.0000660.0000390.022050.019810.0550.0592.92.776.751.60.04710.0417
FA-optimized-PID
 1st Jan0.33940.013510.27260.26010.004860.004051.9251.9860.06868.248.4748859.87.85−0.1098−0.09428
 30th Mar0.23740.019690.18210.17380.002170.001811.2861.3270.06868.248.4748859.87.85−0.0733−0.0629
 13th Dec0.15750.011920.12180.11620.000970.000810.860.88720.06868.248.4748859.87.85−0.049−0.042
 24th May−0.22750.01850.20580.19640.002770.00231.4531.4990.06868.248.4748859.87.850.08280.0712
 28th Sep−0.14990.014510.13750.13120.0012370.001030.97111.0020.06868.248.4748859.87.850.05540.0476

Frequency deviation parameter indices using integral square error objective function.

FIGURE 6

Considering the effect of BESS/PHEV with its different capacities, gain for 0f 10 and 50 are taken to check the stability of power system. Results indicate that settling time , peak overshoot are reduced to 3.4, 80 for GA, 2.64, 74 for PSO, and 2.9, 76.7 for ABC for of 10. These values are further reduced to 2.77, 55.9 for GA, 2.2, 47.8 for PSO, and 2.7, 51.6 for ABC for of 50. Similarly, the values of IAE, ISE, ITAE and are reduced to 0.02258, 0.0002532, 0.03997, and −0.08901 for GA, 0.01997, 0.0002766, 0.02497, and −0.09794 for PSO, 0.02323, 0.00026, 0.0437, −0.09344 for ABC, for of 10, and further reduced to 0.01686, 0.0001398, 0.03458, and −0.07942 for GA, 0.01475, 0.0001646, 0.02049, and −0.085859 for PSO, 0.01835, 0.000155, 0.03927, and −0.08275 for ABC, for of 50 for 1st January loading conditions. Rise time remains almost unchanged in any case. FA optimized PID produces high error values although it is connected to BESS/PHEV support. Instead, the settling time reached to 488 s when BESS/PHEV gain is 50. It is concluded that GA optimized PID produces comparatively better result that PSO, ABC, and FA when ISE cost function is used for the load frequency control.

Frequency deviation responses of the GA optimized PID, PSO optimized PID, ABC optimized PID, and FA optimized PID using the ITAE objective function is shown in Table 8 and Figure 7. Their performance indices show that GA optimized PID, PSO optimized PID, and ABC optimized PID produce approximately efficient competitive results. IAE, ISE, ITAE and are higher for higher mismatches between load and photovoltaic generated power and lower for lower mismatches, i.e., The values of IAE, ISE, ITAE and are 0.02375, 0.0005794, 0.01643, and −0.1374 for GA, 0.02201, 0.0003988, 0.01991, and −0.1136 for PSO, 0.022, 0.000409, 0.01936, and −0.1165 for ABC, considering the 0.3394 pu load deviation and 0.01351 pu solar PV power on 1st January in a single-area power system. These values of IAE, ISE, ITAE and are lowered to 0.01587, 0.0002587, 0.01098, and −0.09176 for GA, 0.0147, 0.000178, 0.0133, and −0.07594 for PSO, 0.0147, 0.0001827, 0.01293 and −0.0778 for ABC, for 0.2374 pu load deviation and 0.0147 pu solar PV power on 13th December loading conditions. The values of rise time , settling , and peak overshoot remain same for different loading conditions. For example, the values of rise , settling , and peak overshoot are 0.0827, 2.55, and 73.4 for GA, 0.0665, 3.03, and 78.4 for PSO, 0.0686, 2.76, and 78.2 for ABC on 1st January loading conditions. FA optimized PID produces relatively poor result as the IAE, ISE, ITAE, and values are much higher than GA, PSO as can be seen the values are 0.07121, 0.005486, 0.04984, and −0.3024.

TABLE 8

Without BESS/PHEV support
DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.7180.38936.5020.9127.978.15−0.9184
 30th Mar0.23740.019691.1480.17384.3450.9127.978.15−0.6136
 13th Dec0.15750.011920.76740.077722.9050.9127.978.15−0.4103
 24th May−0.22750.01851.2970.22194.9090.9127.978.150.6934
 28th Sep−0.14990.014510.86660.09913.2810.9127.978.150.4634
GA-optimized-PID
 1st Jan0.33940.013510.023750.00057940.016430.08272.5573.4−0.1374
 30th Mar0.23740.019690.015870.00025870.010980.08272.5573.4−0.09176
 13th Dec0.15750.011920.010610.00011570.0073390.08272.5573.4−0.06134
 24th May−0.22750.01850.017930.00033020.01240.08272.5573.40.1037
 28th Sep−0.14990.014510.011980.00014750.0082870.08272.5573.40.06925
PSO-optimized-PID
 1st Jan0.33940.013510.022010.00039880.019910.06653.0378.4−0.1136
 30th Mar0.23740.019690.01470.0001780.01330.06653.0378.4−0.07594
 13th Dec0.15750.011920.0093880.00007250.0084930.06653.0378.4−0.04849
 24th May−0.22750.01850.016620.00022730.015030.06653.0378.40.08579
 28th Sep−0.14990.014510.01110.00010150.010040.06653.0378.40.05736
ABC-optimized-PID
 1st Jan0.33940.013510.0220.0004090.019360.06862.7678.2−0.1165
 30th Mar0.23740.019690.01470.00018270.012930.06862.7678.2−0.0778
 13th Dec0.15750.011920.009820.00008170.008640.06862.7678.2−0.052
 24th May−0.22750.01850.016610.0002330.014620.06862.7678.20.088
 28th Sep−0.14990.014510.01110.0001040.009760.06862.7678.20.0588
FA-optimized-PID
 1st Jan0.33940.013510.071210.0054680.049840.2051.8445.5−0.3024
 30th Mar0.23740.019690.047580.0024410.03330.2051.8445.5−0.202
 13th Dec0.15750.011920.031820.001090.02220.2051.8445.5−0.135
 24th May−0.22750.01850.053760.0031170.037630.2051.8445.50.228
 28th Sep−0.14990.014510.035930.001390.025150.2051.8445.50.1526
With BESS/PHEV support, and
DateIAEISEITAE
Conventional-PID
 1st Jan0.33940.013511.721.2990.23560.0895610.2411.277.93—693—10.1—−0.6085−0.2554
 30th Mar0.23740.019691.1490.86810.10520.039986.8397.537.93—693—10.1—−0.4066−0.1705
 13th Dec0.15750.011920.78850.58050.047030.017884.5735.0357.93—693—10.1—−0.2719−0.1140
 24th May−0.22750.01851.2990.98090.13430.051057.7288.5097.93—693—10.1—0.45940.1928
 28th Sep−0.14990.014510.86780.65550.059970.02285.1645.6867.93—693—10.1—0.30700.1289
GA-optimized-PID
 1st Jan0.33940.013510.019680.014120.00046870.00028010.012840.0097430.08440.0952.040.9463.129.9−0.1302−0.1080
 30th Mar0.23740.019690.013150.0094360.00020920.0001250.0085810.006510.08440.0952.040.9463.129.9−0.08699−0.07205
 13th Dec0.15750.011920.0087930.0063090.00009350.00005590.0057380.0043530.08440.0952.040.9463.129.9−0.05821−0.04821
 24th May−0.22750.01850.014860.010660.00026720.00015960.0096960.0073560.08440.0952.040.9463.129.90.098310.08150
 28th Sep−0.14990.014510.0099290.0071250.00011930.00007120.0064790.0049150.08440.0952.040.9463.129.90.065690.0544
PSO-optimized-PID
 1st Jan0.33940.013510.018630.013470.00032070.00019190.016730.013340.06470.0722.251.9769.841.4−0.1092−0.09385
 30th Mar0.23740.019690.012450.0090030.00014320.00008560.011180.0089120.06470.0722.251.9769.841.4−0.07292−0.06276
 13th Dec0.15750.011920.0079470.0057480.00005830.00003490.0071380.005690.06470.0722.251.9769.841.4−0.04656−0.04008
 24th May−0.22750.01850.014060.010170.00018280.00010940.012630.010070.06470.0722.251.9769.841.40.082400.07091
 28th Sep−0.14990.014510.0093990.0067980.00008160.00004880.0084420.006730.06470.0722.251.9769.841.40.055070.04742
ABC-optimized-PID
 1st Jan0.33940.013510.018770.013520.0003290.00019550.016530.013570.06960.0752.31.7569.440.3−0.116−0.095
 30th Mar0.23740.019690.012540.0090350.0001470.0000870.011040.0090650.06960.0752.31.7569.440.3−0.074−0.064
 13th Dec0.15750.011920.008380.006040.0000650.0000390.007380.006060.06960.0752.31.7569.440.3−0.0499−0.0427
 24th May−0.22750.01850.014170.010210.00018780.0001110.012480.010240.06960.0752.31.7569.440.30.08430.0722
 28th Sep−0.14990.014510.009470.006820.0000830.0000490.008330.0068450.06960.0752.31.7569.440.30.05640.04823
FA-optimized-PID
 1st Jan0.33940.013510.068250.05750.004250.002200.051320.051320.2271.121.891.7925.61.8−0.264−0.175
 30th Mar0.23740.019690.04560.038420.0018970.0009850.034290.034290.2271.121.891.7925.61.8−0.176−0.1169
 13th Dec0.15750.011920.03040.025690.0008480.000440.02290.022930.2271.121.891.7925.61.8−0.118−0.0781
 24th May−0.22750.01850.05150.043410.002420.0012590.038750.038740.2271.121.891.7925.61.80.20.132
 28th Sep−0.14990.014510.034430.029010.001080.0005620.025890.025890.2271.121.891.7925.61.80.1330.0882

Frequency deviation parameter indices using integral time absolute error objective function.

FIGURE 7

BESS/PHEV with its different capacities, gain 10 and 50 are taken into consideration. Results depict that settling time , peak overshoot are reduced to 2.04, 63.1 for GA, 2.25, 69.8 for PSO, and 2.3, 69.4 for ABC for of 10. These values are further reduced to 0.94, 29.9 for GA, 1.97, 41.4 for PSO, and 1.75, 40.3 for ABC for of 50. Similarly, the values of IAE, ISE, ITAE and are reduced to 0.01968, 0.0004687, 0.01284, and −0.1302 for GA, 0.01347, 0.0001919, 0.01673, and −0.1092 for PSO, 0.01877, 0.000329, 0.01653, and −0.116 for ABC, for of 10, and further reduced to 0.01412, 0.0002801, 0.009743, and −0.1080 for GA, 0.01475, 0.0001646, 0.01334, and −0.09385 for PSO, 0.01352, 0.0001955, 0.01357, and −0.095 for ABC, for of 50 for 1st January loading conditions. Rise time remains almost unchanged in any case. FA optimized PID produces high error values and frequency deviation although it is connected to BESS/PHEV support. It is concluded that GA optimized produces better results with respect to settling time and peak overshoot but PSO optimized PID produces minimum IAE, ISE, ITAE and using ITAE objective function.

The conventional PID controller is tuned using the PID tuner toolbox within MATLAB. Tables 6–8 depict that performance indices IAE, ISE, ITAE and increase with the higher load changes and decrease with lower load changes. But rise time , settling time , and peak overshoot remains same regardless of the load and photovoltaic power changes. It is depicted that performance indices IAE and ISE decreases when BESS/PHEV is integrated to the system at the gain levels 10 and 50, while ITAE increases when the BESS/PHEV is integrated. Frequency deviation decreases with the attachment of BESS/PHEV support. Rise time , settling time , and peak overshoot increase significantly for gain 10, but are infinite at gain 50 because of the unstable response of conventional PID controller. Therefore, the values are replaced by – (dash) where the conventional PID controller response is unstable. It is considered to be major drawback of the conventional PID controller. The conventional PID controller response is shown in Figures 5–7.

The maximum steady state frequency deviation permissible limit in Hz is found to be in the range (−0.7932, 0.5468) for maximum positive and negative load deviation respectively. It can be seen from Table 6, Table 7 and Table 8 that the frequency deviation in Hz has found to be within permissible range using the GA, PSO, ABC, and FA optimization methods. The frequency deviation has not found to be within steady state limit for Conventional PID controller. The frequency deviation range in Hz is (−0.9184, 0.4634) for Conventional-PID. The frequency deviation is suppressed by adding the BESS/PHEV source. But, the PID controller response becomes unstable if the large capacity of BESS/PHEV source is attached to the power system model. The power delivered to power system by BESS and/or PHEV is 0.025 pu for gain of 10, and 0.105 pu for gain of 50 taking into account the maximum positive loading condition on 1st January.

After comparison between the applied techniques and objective functions, it is concluded that PSO optimized PID controller using the ITAE objective function performs well just to mitigate frequency deviation of the power system but GA optimized PID controller using the ISE objective function gives better overall results considering all the performance indices. The convergence characteristics and parameter indices indicate the effectiveness, sustainability and consistency of proposed GA optimized PID tested within the loading range (−0.2275, 0.3394).

Conclusion

The load deviation of single area power system and solar photovoltaic power is forecasted using the state-of-the art artificial neural network. PSO, GA, ABC, and FA metaheuristic techniques have been employed to find out the optimal PID controller gain parameters using the IAE, ISE and ITAE objective functions. Some key days, showing significant load changes, have been chosen to see the performance of the single area hybrid power system. Results indicate the efficiency and competency of GA, PSO, and ABC among each other that sufficiently reduce frequency fluctuations. It has been found that PSO optimized PID along with ITAE cost function can sufficiently minimize the frequency. The performance indices and convergence characteristics show overall better results produced from GA optimized PID than others using ISE objective function. These frequency fluctuations are further improved by integrating power from BESS/PHEV source. Small BESS/PHEV capacity is required to damp the frequency when integrated with GA optimized PID.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://ntdc.gov.pk/services.

Author contributions

SM and AM conceived the conceptualization and methodology. SM conducted the normal analysis and writing. MR and AM conducted the review and editing. MR managed funding acquisition. All authors have agreed to submitted version of the manuscript.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IFPIP-1183-135-1442) and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Nomenclature

  • IAE

    Integral absolute error

  • ISE

    Integral square error

  • ITAE

    Integral time absolute error

  • GA

    Genetic algorithm

  • PSO

    Particle swarm optimization

  • ABC

    Artificial bee colony

  • FA

    Firefly algorithm

  • Rise time

  • Settling time

  • Peak overshoot

  • Maximum frequency deviation in Hertz

  • Proportional gain

  • Integral gain

  • Derivative gain

  • ACE

    Area control error

  • Change in mechanical power

  • Solar photovoltaic power

  • Net change in electrical load demand

  • Battery Energy Storage System/Plug—in Hybrid Electric Vehicle Power

  • BESS/PHEV gain

  • DC-AC converter time constant

  • BESS/PHEV charging time constant

  • Rated power plant power

  • Change in nominal frequency

  • Generator inertia constant

  • Change in frequency sensitive load power

  • Turbine time constant

  • Governor time constant

  • Ratio of percent change in load to percent change in frequency

  • Governor speed regulation per unit

  • Attractiveness parameter

  • Coefficient of absorption

References

  • 1

    Abd-ElazimS. M.AliE. S. (2018). Load Frequency Controller Design of a Two-Area System Composing of PV Grid and Thermal Generator via Firefly Algorithm. Neural Comput. Applic30 (2), 607–616. 10.1007/s00521-016-2668-y

  • 2

    Abo-ElyousrF. K. (2018). “Load Frequency Controller with Virtual Inertia Generator for Interconnected Power Systems via Artificial Bee Colony,” in 2018 Twentieth International Middle East Power Systems Conference (MEPCON) (Giza, Egypt: IEEE). 10.1109/mepcon.2018.8635143

  • 3

    AdityaS. K.DasD. (2001). Battery Energy Storage for Load Frequency Control of an Interconnected Power System. Electr. power Syst. Res.58 (3), 179–185. 10.1016/s0378-7796(01)00129-8

  • 4

    AkulaS. K. (2019). Frequency Control in Microgrid Communities Using Neural Networks. Wichita, KS, USA: IEEE.

  • 5

    Al-DahidiS.AyadiO.AlrbaiM.AdeebJ. (2019). Ensemble Approach of Optimized Artificial Neural Networks for Solar Photovoltaic Power Prediction. IEEE Access7, 2923905. 10.1109/access.2019.2923905

  • 6

    AlessandraB. A.Carlos RobertoM.Mara LuciaM. L.Fábio RobertoC.Nelson JoseP. (2011). Multinodal Load Forecasting in Power Electric Systems Using a Neural Network with Radial Basis Function. Adv. Mater. Res.217, 39. 10.4028/www.scientific.net/amr.217-218.39

  • 7

    AzeerS. A.Ramjug-BallgobinR.Sayed HassenS. Z. (2017). Intelligent Controllers for Load Frequency Control of Two-Area Power System. IFAC-PapersOnLine50 (2), 301–306. 10.1016/j.ifacol.2017.12.062

  • 8

    BoddepalliM. K.NavuriP. K. (2018). “Design and Analysis of Firefly Algorithm Based PID Controller for Automatic Load Frequency Control Problem,” in 2018 Technologies for Smart-City Energy Security and Power (ICSESP) (Bhubaneswar, India: IEEE).

  • 9

    Chandra SekharG. T.SahuR. K.BaliarsinghA. K.PandaS. (2016). Load Frequency Control of Power System under Deregulated Environment Using Optimal Firefly Algorithm. Int. J. Electr. Power & Energy Syst.74, 195–211. 10.1016/j.ijepes.2015.07.025

  • 10

    ChangC. S.FuW.Wen.F. (1998). Load Frequency Control Using Genetic-Algorithm Based Fuzzy Gain Scheduling of PI Controllers. Electr. Mach. Power Syst.26 (1), 39–52. 10.1080/07313569808955806

  • 11

    DasD. C.RoyA. K.SinhaN. (2012). GA Based Frequency Controller for Solar Thermal-Diesel-Wind Hybrid Energy Generation/energy Storage System. Int. J. Electr. Power & Energy Syst.43 (1), 262–279. 10.1016/j.ijepes.2012.05.025

  • 12

    DasD. C.RoyA. K.SinhaN. (2010). “Genetic Algorithm Based PI Controller for Frequency Control of an Autonomous Hybrid Generation System,” in World Congress on Engineering 2012, July 4-6, 2012 (London, UK: International Association of Engineers). Vol. 2189.

  • 13

    ElsisiM.SolimanM.AboelelaM. A. S.MansourW. (2015). ABC Based Design of PID Controller for Two Area Load Frequency Control with Nonlinearities. Telkomnika Indonesian J. Electr. Eng.16 (1), 58–64. 10.11591/tijee.v16i1.1588

  • 14

    GhasemiH.ShayeghiA. (2011). Market Based LFC Design Using Artificial Bee Colony. Int. J. Tech. Phys. Problems Eng.3 (6), 1–10.

  • 15

    GozdeH.Cengiz TaplamaciogluM.Kocaarslanİ. (2012). Comparative Performance Analysis of Artificial Bee Colony Algorithm in Automatic Generation Control for Interconnected Reheat Thermal Power System. Int. J. Electr. Power & Energy Syst.42 (1), 167–178. 10.1016/j.ijepes.2012.03.039

  • 16

    GözdeH.TaplamacioğluM. C.Kocaarslanİ.ÇamE. (2008). Particle Swarm Optimization Based Load Frequency Control in a Single Area Power System. Sci. Bull.2 (8).

  • 17

    GuptaD. K.SoniA. K.JhaA. V.MishraS. K.AppasaniB.SrinivasuluA.et al (2021). Hybrid Gravitational-Firefly Algorithm-Based Load Frequency Control for Hydrothermal Two-Area System. Mathematics9 (7), 712. 10.3390/math9070712

  • 18

    Hazlee AzilI.Mohamed ZahariA. F.MokhlisH. (2016). Optimisation of PID Controller for Load Frequency Control in Two-Area Power System Using Evolutionary Particle Swarm Optimisation. J. Electr. Syst.12 (2), 315–324.

  • 19

    HemeidaA.MohamedS.MahmoudM. (2020). Load Frequency Control Using Optimized Control Techniques. JES. J. Eng. Sci.48 (6), 1119–11136. 10.21608/jesaun.2020.42349.1011

  • 20

    HoteY. V. (2018). PID Controller Design for Load Frequency Control: Past, Present and Future Challenges. IFAC Pap.51, 608. 10.1016/j.ifacol.2018.06.162

  • 21

    HuangC.BensoussanA.EdesessM.TsuiK. L. (2016). Improvement in Artificial Neural Network-Based Estimation of Grid Connected Photovoltaic Power Output. Renew. Energy97, 838–848. 10.1016/j.renene.2016.06.043

  • 22

    JaberA. S.AhmadA. Z.AbdallaA. N. (2013). An Investigation of Scaled-FLC Using PSO for Multi-Area Power System Load Frequency Control. Energy Power Eng.05 (04), 458–462. 10.4236/epe.2013.54b088

  • 23

    JeyalakshmiV.SubburajP. (2016). PSO-scaled Fuzzy Logic to Load Frequency Control in Hydrothermal Power System. Soft Comput.20 (7), 2577–2594. 10.1007/s00500-015-1659-8

  • 24

    KalyaniS.NagalakshmiS.MarishaR. (2012). “Load Frequency Control Using Battery Energy Storage System in Interconnected Power System,” in 2012 Third International Conference on Computing, Communication and Networking Technologies (ICCCNT'12) (Coimbatore, India: IEEE). 10.1109/icccnt.2012.6396052

  • 25

    KhatoonS.SinghA. K. (2014). “Effects of Various Factors on Electric Load Forecasting: An Overview,” in 2014 6th IEEE Power India International Conference (PIICON) (Delhi, India: IEEE). 10.1109/poweri.2014.7117763

  • 26

    KonarG.MandalK. K.ChakrabortyN. (2014). “Two Area Load Frequency Control Using GA Tuned PID Controller in Deregulated Environment,” in Proceedings of the International Multi Conference of Engineers and Computer Scientists, 2.

  • 27

    KoturD.ŽarkovićM. (2016). “Neural Network Models for Electricity Prices and Loads Short and Long-Term Prediction,” in 2016 4th International Symposium on Environmental Friendly Energies and Applications (EFEA) (Belgrade, Serbia: IEEE). 10.1109/efea.2016.7748787

  • 28

    KoubaN. E. Y.MenaaM.HasniM.BoudourM. (2015). “Optimal Load Frequency Control Based on Artificial Bee Colony Optimization Applied to Single, Two and Multi-Area Interconnected Power Systems,” in 2015 3rd International Conference on Control, Engineering & Information Technology (CEIT) (Tlemcen, Algeria: IEEE). 10.1109/ceit.2015.7233027

  • 29

    KumarD.MishraA.ChatterjeeK. (2017). “Power and Frequency Control of a Wind Energy Power System Using Artificial Bee Colony Algorithm,” in 2017 Third International Conference on Science Technology Engineering & Management (ICONSTEM) (Chennai, India: IEEE). 10.1109/iconstem.2017.8261385

  • 30

    KumariK.ShankarG.KumariS.GuptaS. (2016). “Load Frequency Control Using ANN-PID Controller,” in 2016 IEEE 1st International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES) (Delhi, India: IEEE). 10.1109/icpeices.2016.7853516

  • 31

    KumariN.JhaA. N. (2014). Frequency Response Enhancement of Hybrid Power System by Using PI Controller Tuned with PSO Technique. Int. J. Adv. Comput. Res.4 (1), 116.

  • 32

    LiangL.JinZ.JiaoZ. (2012). “Frequency Regulation for a Power System with Wind Power and Battery Energy Storage,” in 2012 IEEE International Conference on Power System Technology (POWERCON) (Auckland, New Zealand: IEEE). 10.1109/powercon.2012.6401357

  • 33

    LoneA. H.YousufV.PrakashS.bid BazazM. A. (2018). “Load Frequency Control of Two Area Interconnected Power System Using SSSC with PID, Fuzzy and Neural Network Based Controllers,” in 2018 2nd IEEE International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES) (Delhi, India: IEEE). 10.1109/icpeices.2018.8897468

  • 34

    MalleshamG.MishraS.JhaA. N. (2012). “Automatic Generation Control of Microgrid Using Artificial Intelligence Techniques,” in 2012 IEEE Power and Energy Society General Meeting (San Diego, CA, USA: IEEE).

  • 35

  • 36

  • 37

    MellitA.PavanA. M. (2010). Performance Prediction of 20kWp Grid-Connected Photovoltaic Plant at Trieste (Italy) Using Artificial Neural Network. Energy Convers. Manag.51, 2431–2441. 10.1016/j.enconman.2010.05.007

  • 38

    MilaniA. E.MozafariB. (2009). Genetic Algorithm Based Optimal Load Frequency Control in Two‐Area Interconected Power Systems. AIP Conf. Proc.1159 (1), 43. 10.1063/1.3223953

  • 39

    ModiN.KhareM.ChaturvediK. (2013). Performance Analysis of Load Frequency Control in Single Area Power System Using GA and PSO Based PID Controller. Int. J. Electr. Electron. Comput. Eng.2 (1), 108–114.

  • 40

    MosaadM. I.Salem.F. (2014). LFC Based Adaptive PID Controller Using ANN and ANFIS Techniques. J. Electr. Syst. Inf. Technol.1 (3), 212–222. 10.1016/j.jesit.2014.12.004

  • 41

    NagarjunaN.ShankarG. (2015). “Load Frequency Control of Two Area Power System with AC-DC Tie Line Using PSO Optimized Controller,” in 2015 International Conference on Power and Advanced Control Engineering (ICPACE) (Bengaluru, India: IEEE). 10.1109/icpace.2015.7274948

  • 42

    NaiduK.MokhlisH.BakarA. H. A. (2013). “Application of Firefly Algorithm (FA) Based Optimization in Load Frequency Control for Interconnected Reheat Thermal Power System,” in 2013 IEEE Jordan Conference on Applied Electrical Engineering and Computing Technologies (AEECT) (Amman, Jordan: IEEE). 10.1109/aeect.2013.6716461

  • 43

    NaiduK.MokhlisH.BakarA. H. A. (2014). Multiobjective Optimization Using Weighted Sum Artificial Bee Colony Algorithm for Load Frequency Control. Int. J. Electr. Power & Energy Syst.55, 657–667. 10.1016/j.ijepes.2013.10.022

  • 44

    NcaneZ. P.SahaA. K. (2019). “Forecasting Solar Power Generation Using Fuzzy Logic and Artificial Neural Network,” in 2019 Southern African Universities Power Engineering Conference/Robotics and Mechatronics/Pattern Recognition Association of South Africa (SAUPEC/RobMech/PRASA) (Bloemfontein, South Africa: IEEE). 10.1109/robomech.2019.8704737

  • 45

  • 46

    OtaniT. (2017). “Cooperative Load Frequency Control of Generator and Battery Using a Recurrent Neural Network,” in Proc. of the 2017 IEEE Region 10 Conference (TENCON) (Malaysia: TENCON), 918–923. 10.1109/tencon.2017.8227989

  • 47

    PadhanS.SahuR. K.PandaS. (2014). Application of Firefly Algorithm for Load Frequency Control of Multi-Area Interconnected Power System. Electr. power components Syst.42 (13), 1419–1430. 10.1080/15325008.2014.933372

  • 48

    PainS.AcharjeeP. (2014). Multiobjective Optimization of Load Frequency Control Using Pso. Int. J. Emerg. Technol. Adv. Eng.4 (7), 16–22.

  • 49

    Power Data Reference Book (2017). Power Data Reference Book, Vol 2, 2010-11 to 2016-17. Lahore: Planning Power National Transmission and Despatch Company Limited (NTDCL). Available at: https://ntdc.gov.pk/services.

  • 50

    PrajapatiP.ParmarA. (2016). “Multi-area Load Frequency Control by Various Conventional Controller Using Battery Energy Storage System,” in 2016 International Conference on Energy Efficient Technologies for Sustainability (ICEETS) (Nagercoil, India: IEEE). 10.1109/iceets.2016.7583800

  • 51

    PrakashS.SinhaS. K. (2012). “Four Area Load Frequency Control of Interconnected Hydro-Thermal Power System by Intelligent PID Control Technique,” in 2012 Students Conference on Engineering and Systems (Allahabad, India: IEEE). 10.1109/sces.2012.6199090

  • 52

    QiX.BaiY.LuoH.ZhangY.ZhouG.WeiZ. (2018). Fully-distributed Load Frequency Control Strategy in an Islanded Microgrid Considering Plug-In Electric Vehicles. Energies11 (6), 1613. 10.3390/en11061613

  • 53

    RaoC. S. (2012). Design of Artificial Intelligent Controller for Automatic Generation Control of Two Area Hydrothermal System. Int. J. Electr. Comput. Eng.2 (2), 183.

  • 54

    RaoR. N.Rama Krishna ReddyP. (2015). PSO Based Tuning of PID Controller for a Load Frequency Control in Two Area Power System. Int. J. Eng. Res. Appl.1 (3), 1499–1505.

  • 55

    RathorS.AcharyaD. S.GudeS.MishraP. (2011). “Application of Artificial Bee Colony Optimization for Load Frequency Control,” in 2011 World Congress on Information and Communication Technologies (Mumbai, India: IEEE). 10.1109/wict.2011.6141339

  • 56

    RegadM.HelaimiM.TalebR.GabbarH. A.OthmanA. M. (2019). “Fractional Order PID Control of Hybrid Power System with Renewable Generation Using Genetic Algorithm,” in 2019 IEEE 7th International Conference on Smart Energy Grid Engineering (SEGE) (Oshawa, ON, Canada: IEEE). 10.1109/sege.2019.8859970

  • 57

    RerkpreedapongD.HasanovicA.FeliachiA. (2003). Robust Load Frequency Control Using Genetic Algorithms and Linear Matrix Inequalities. IEEE Trans. Power Syst.18 (2), 855–861. 10.1109/tpwrs.2003.811005

  • 58

    RodríguezF.FleetwoodA.GalarzaA.FontánL. (2018). Predicting Solar Energy Generation through Artificial Neural Networks Using Weather Forecasts for Microgrid Control. Renew. energy126, 855–864. 10.1016/j.renene.2018.03.070

  • 59

    SadatH. (2012). Power System Analysis and Design, 33–41.

  • 60

    SafariA.BabaeiF.FarrokhifarM. (2021). A Load Frequency Control Using a PSO-Based ANN for Micro-grids in the Presence of Electric Vehicles. Int. J. Ambient Energy42 (6), 688–700. 10.1080/01430750.2018.1563811

  • 61

    SedaghatiF.NahavandiA.Ali BadamchizadehM.GhaemiS.Abedinpour FallahM. (2012). PV Maximum Power-point Tracking by Using Artificial Neural Network. Math. Problems Eng.2012, 506709. 10.1155/2012/506709

  • 62

    SelvakumaranS.ParthasarathyS.KarthigaivelR.RajasekaranV. (2012). Optimal Decentralized Load Frequency Control in a Parallel AC-DC Interconnected Power System through HVDC Link Using PSO Algorithm. Energy Procedia14, 1849–1854. 10.1016/j.egypro.2011.12.1178

  • 63

    ShakaramiM. R.FarajiI.AsghariI.AkbariM. (2013). “Optimal PID Tuning for Load Frequency Control Using Lévy-Flight Firefly Algorithm,” in 2013 3rd International Conference on Electric Power and Energy Conversion Systems (Istanbul, Turkey: IEEE). 10.1109/EPECS.2013.6713008

  • 64

    ShankarG.LakshmiS.NagarjunaN. (2015). “Optimal Load Frequency Control of Hybrid Renewable Energy System Using PSO and LQR,” in 2015 International Conference on Power and Advanced Control Engineering (ICPACE) (Bengaluru, India: IEEE). 10.1109/icpace.2015.7274942

  • 65

    Sidi BrahimR. M.M’hamedH.TalebR.LemraboutA. (2021). Load Frequency Control of Hybrid Power System Using Classical PID Controller. iKSP J. Comput. Sci. Eng.1 (1).

  • 66

    SinghRRameshL. (2019). Comparison of Automatic Load Frequency Control in Two Area Power Systems Using Pso Algorithm Based Pid Controller and Conventional Pid Controller. J. Phys. Conf. Ser.1172, 012054. 10.1088/1742-6596/1172/1/012054

  • 67

    SrinivasanD.LiewA. C.ChenJ. S. P. (1991). “Short Term Forecasting Using Neural Network Approach,” in Proceedings of the First International Forum on Applications of Neural Networks to Power Systems (Piscataway, NJ, USA: IEEE).

  • 68

    SundaramV. S.JayabarathiT. (2011). “Load Frequency Control Using PID Tuned ANN Controller in Power System,” in 2011 1st International Conference on Electrical Energy Systems (Chennai, India: IEEE). 10.1109/icees.2011.5725341

  • 69

    TaherS. A. (2014). Fractional Order PID Controller Design for LFC in Electric Power Systems Using Imperialist Competitive Algorithm. Ain Shams Eng. J.5, 133–135. 10.1016/j.asej.2013.07.006

  • 70

    TakagiM.YamamotoH.YamajiK.OkanoK.HiwatariR.IkeyaT. (2009). Load Frequency Control Method by Charge Control for Plug-In Hybrid Electric Vehicles with LFC Signal. IEEJ Trans. Power Energy129 (11), 1342–1348. 10.1541/ieejpes.129.1342

  • 71

    UeharaA.SenjyuT.YonaA.FunabashiT. (2009). “Frequency Control by Coordination Control of WTG and Battery Using Load Estimation,” in 2009 International Conference on Power Electronics and Drive Systems (PEDS) (Taipei, Taiwan: IEEE). 10.1109/peds.2009.5385765

  • 72

    VeerasamyV.WahabN. I. A.RamachandranR.OthmanM. L.HizamH.IrudayarajA. X. R.et al (2020). A Hankel Matrix Based Reduced Order Model for Stability Analysis of Hybrid Power System Using PSO-GSA Optimized Cascade PI-PD Controller for Automatic Load Frequency Control. IEEE Access8, 71422–71446. 10.1109/access.2020.2987387

  • 73

    Yang-WuS.XunM.AoP.Yang-GuangW.TingC.DingW.et al (2019). “Load Frequency Control Strategy for Wind Power Grid-Connected Power Systems Considering Wind Power Forecast,” in 2019 IEEE 3rd Conference on Energy Internet and Energy System Integration (EI2) (Changsha, China: IEEE). 10.1109/EI247390.2019.9062084

Appendix

The typical values for the hybrid power system are given below

Summary

Keywords

load frequency control, PID controller design, load deviation forecasting, solar photovoltaic power forecasting, particle swarm optimization, genetic algorithm, artificial bee colony (ABC) algorithm

Citation

Mahboob Ul Hassan S, Ramli MAM and Milyani AH (2022) Robust Load Frequency Control of Hybrid Solar Power Systems Using Optimization Techniques. Front. Energy Res. 10:902776. doi: 10.3389/fenrg.2022.902776

Received

23 March 2022

Accepted

04 May 2022

Published

14 July 2022

Volume

10 - 2022

Edited by

Chuan-Ke Zhang, China University of Geosciences Wuhan, China

Reviewed by

D. Boopathi, Paavai Engineering College, India

Hazlie Mokhlis, University of Malaya, Malaysia

Updates

Copyright

*Correspondence: Syed Mahboob Ul Hassan, ,

This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics