Next Article in Journal
Enhancing Barite Carbothermal Reduction with Brown Coal by Compaction of the Charge
Next Article in Special Issue
The Mixture of Probability Distribution Functions for Wind and Photovoltaic Power Systems Using a Metaheuristic Method
Previous Article in Journal
The Fuzzy DEA-Based Manufacturing Service Efficiency Evaluation and Ranking Approach for a Parallel Two-Stage Structure of a Complex Product System on the Example of Solid Waste Recycling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving the Frequency Response of Hybrid Microgrid under Renewable Sources’ Uncertainties Using a Robust LFC-Based African Vulture Optimization Algorithm

1
Department of Electrical Power and Machines, Faculty of Engineering, Alexandria University, Alexandria 21544, Egypt
2
Electrical Engineering Department, Faculty of Engineering, Aswan University, Aswan 81542, Egypt
*
Authors to whom correspondence should be addressed.
Processes 2022, 10(11), 2320; https://doi.org/10.3390/pr10112320
Submission received: 27 September 2022 / Revised: 17 October 2022 / Accepted: 25 October 2022 / Published: 8 November 2022
(This article belongs to the Special Issue Sustainable Microgrid Systems: Technologies, Applications and Trends)

Abstract

:
Power systems have recently faced significant challenges due to the increased penetration of renewable energy sources (RES) such as frequency deviation due to fluctuations, unpredictable nature, and uncertainty of this RES. In this paper, a cascaded controller called (1+PD)-PID is proposed to reduce the influence of RES uncertainties on the system and to maintain the system’s reliability during fluctuations. The proposed controller is a combination of (1+PD) and PID controllers in order. The output signal of the (1+PD) controller along with the frequency deviation and the power difference between adjacent areas are used as inputs to the PID controller to create the load reference signal. The parameters of the suggested controller are optimally tuned using the African Vulture Optimization Algorithm (AVOA) to ensure the best performance of the controller. A two-area interconnected system with non-reheat thermal power units combined with RES such as solar and wind energy is modeled using MATLAB/Simulink to evaluate the system response. The controller effectiveness is verified by subjecting the studied system to various types of fluctuations such as step load disturbance, variable load perturbation and RES penetration. The obtained simulation results prove that the proposed (1+PD)-PID controller in integration with AVOA offers a significant improvement in the system performance specifications. Moreover, the proposed AVOA-based (1+PD)-PID controller has proven its superiority over other comparable controllers having the least fitness function of 6.01 × 10−5.

1. Introduction

Recently, renewable energy sources (RES) have become of great importance due to several reasons such as global warming caused by emissions from conventional sources and increasing demand for more electric power. Although renewable sources are considered environmentally friendly, there are various limitations to their use, including low total inertia and fault ride, uncertainty due to climatic conditions, voltage and frequency fluctuations, and poor stability [1]. Modern electric power systems consist of several interconnected areas characterized by renewable and non-renewable (conventional) energy systems that create a hybrid power system known as a Microgrid (MG). It is a feasible approach in this direction to support the integration of sustainable energy sources with conventional sources to completely reduce the dependence on conventional energy. The MG consists of small generating and load units placed at different locations in the network. The MG can work either in integration mode or in isolation mode. The maximum load size for any MG is 10 MVA [2].
To ensure the security, stability and sustainability of the integrated power system equipment, the total output power generated in each controlled area should meet the total load needs while also accounting for system losses. When the system is subjected to any unexpected variations in load demand, the system must also be able to keep the frequency- and the interconnected area tie-line power within permissible limits [3]. As the load fluctuates over time, the generated power must be automatically updated to meet the demand. This configuration is known as automatic load frequency control (LFC). Thus, the LFC provides frequency balancing by regulating the speed of the prime mover, as well as stabilizing the tie-line power flow during normal and turbulent conditions [4].
Based on the aforementioned issues, the integration of LFC with electric power systems has gained much interest, especially in hybrid systems that feature RES alongside conventional sources due to power fluctuations and uncertainties. Several researchers have proposed different ways to mimic LFC embedded in different power systems in order to deal with this problem. In this field, traditional control techniques, artificial controllers, sliding mode control, and robust control methods have been used. Ease of implementation, low cost and minimal computation time are the reasons behind using traditional controllers such as proportional-integral (PI) controllers and proportional integral derivative (PID) controllers.
In order to obtain the optimized parameters of these conventional controllers, various optimization techniques have been applied such as Genetic Algorithm (GA) [5], Ant Colony Optimization [6], Bacterial Foraging Optimization Algorithm (BFOA) [7], Differential Evolution (DE) [8], particle swarm optimization (PSO) [9], and grasshopper optimization (GOA) [10]. A model predictive control (MPC) optimized using the Multi-Purpose Practical Swarm Optimization (MOPSO) algorithm and adaptive control techniques was introduced as a frequency control mechanism for a renewable energy-powered smart grid [11]. In [12], the Firefly algorithm has been used in a two-zone hybrid system that incorporates electric vehicles (EVs) as controllable using an integrated controller, but it does not cover the development of a PID controller or even a PI controller.
The cascaded controller configurations have recently been used instead of the conventional controller configurations due to their high efficiency and performance. This has led to the use of a variety of types of cascaded controllers to improve the frequency stability of the power system [13,14]. Furthermore, fractional order (FO) based controllers have also been implemented including a FOPID tuned moving inhibitory wave algorithm [15], a modified FOPID controller optimized via the Artificial Ecosystem Optimizer (AEO) [16], a combination of FOPID and a derivative Tilt integration controller [17,18], advanced high-order differential feedback control devices (HODFC, FHODFC) tuned by PSO [19]. Artificial-intelligence-based controls such as fuzzy logic controllers are integrated with the traditional PID controller and optimally regulated by different algorithms such as the Marine Predator Algorithm (MPA) [20] and Particle Swarm Optimization (PSO) technique [21]. A new control method based on the FO Type-2 Fuzzy Logic system was introduced to solve the problem of frequency management associated with renewable energy systems so that, in case of fluctuation of the load, the systems could return to their initial state. The parameters of the controller used were adjusted online depending on the extended Kalman filter [22] or using an optimization algorithm such as PSO [23]. Due to the benefits of the TID controller, which include its excellent ability to reject disturbances and increase the reliability and durability of the power system, it has recently been used to solve the problem of LFC [2,24,25,26].
In addition, another controller architecture was recently examined while designing different types of controllers for the LFC problem which focus on combining two distinct controllers in order to take advantage of both controllers. In [27], electric vehicles (EVs) were used as distributed energy storage units, and an aggregated topology of integral tilt and tilt derivative parameters was introduced to create an optimized TD-TI controller using a new optimized optimization approach, called the Quantum Chaos Game Optimizer (QCGO). Another example of this architecture is the combination of the fuzzy proportional derivative (FPD) and tilt integral derivative (TID) to present a hybrid fuzzy-based controller by combining the fuzzy proportional derivative (FPD) controller with tilt integrated derivative (TID) sequentially to make benefit of their advantages [28].
Thus, a well-designed control unit gained great importance for its integration into the power system. Moreover, system frequency and link line capacity should be kept within permissible limits and system balance should be restored as soon as possible. In this paper, a cascaded controller called (1+PD)-PID is presented as a reliable alternative approach to enhance the stability, sustainability and reliability of a hybrid power system containing conventional and RES power systems including PV and wind energy sources. The proposed cascaded controller parameters are optimized by AVOA as a novel metaheuristic optimization technique.
The contributions, inspiration, and significance of this study can be summarized as follows:
  • Design and simulation of a robust cascaded controller called (1+PD)-PID in order to regulate the system response in terms of frequency and tie-line power deviations;
  • Using a novel AVOA optimization algorithm to find the optimal controller parameters to ensure an optimal behavior of the controller;
  • Testing the effectiveness and validation of the (1+PD)-PID controller by subjecting the microgrid to various types of fluctuations and uncertainties such as distinct step load disturbances, variable load variations, and RES fluctuations;
  • Verifying the superiority of the (1+PD)-PID controller by comparing its performance against that of other controllers such as the conventional PID controller, FOPID controller and TID controller.
The rest of the paper is organized as follows: Section 2 describes the topology of the two-area system where each component is briefly explained. The structure of the cascaded controller (1+PD)-PID is presented in Section 3, as well as the formulation of the problem. Section 4 provides an overview of the optimization algorithm (AVOA). Section 5 is devoted to the detailed results of simulating the system in different scenarios, and the conclusion is presented in Section 6.

2. Structure of The Two-Area Hybrid Power System

To evaluate the performance and efficiency of the proposed (1+PD)-PID controller, a two-area power system including non-reheat thermal power plants was constructed and studied. This system was adopted because it presents a variety of challenges to the suggested controller, providing a complete evaluation of its effectiveness and contribution to the regulation of frequency and tie-line power deviations. The two-area interconnected system under consideration, according to Figure 1, consists of two regions connected by an AC connection line. A non-reheat thermal power plant, a photovoltaic power plant, and associated loads are all found in Area-1, while Area-2 includes a non-reheating thermal power plant, a wind power plant, and loads. Each component in the system is represented by a transfer function of the first order composed of gain and time constant. Both areas of the system include generators with a rated capacity of 2000 MW. The system components are represented by linear mathematical models given using Equations (1)–(4).
  • Non-reheat Thermal System:
The model representing the non-reheat thermal power plant consists of a governor stage followed by a turbine stage. The overall transfer function can be expressed as follows:
G T S ( s ) = [ G g ( s ) ] [ G T ( s ) ] = [ K g 1 + s T g ] [ K T 1 + s T g T ]
where K T and K g are the gains of the turbine and the governor, respectively, and T T and T g represent the time constants of the turbine and the governor, respectively.
  • Power System:
The transfer function representing the power system can be written as follows:
G P S ( s ) = [ K P S 1 + s T P S ]
where K P S and T P S refer to the power system gain and time constant, respectively.
  • PV System:
A first-order transfer function is used to represent the PV system in a linear model as shown in Equation (3). Hence, at constant air temperature, the output power of the PV unit changes with solar radiance linearly.
G P V ( s ) = [ K P V 1 + s T P V ]
where K P V and T P V are the PV power plant gain and time constant, respectively.
In this paper, the input of the PV model is real solar irradiance data obtained from a PV station located in Aswan, Egypt with a capacity of 1.5 GW. Figure 2 shows the input solar irradiance during the day with a maximum of 1000 W/ m 2 and the output PV power with a maximum of 0.25 p.u.
  • Wind Turbine Generator (WTG):
In this study, simplified wind power plant modeling is used to simulate the changing nature of the output power produced by the wind production system and it is expressed as follows:
G W T ( s ) = [ K W T 1 + s T W T ]
where K W T and T W T are the wind power plant gain and time constant, respectively. In order to represent the fluctuations in wind power, the real wind speed data obtained from the wind power plant located in Zafarana, Egypt, are used as input into the wind turbine generator model. Figure 3 shows the fluctuation in wind speed and the change in wind turbine output power.
Table 1 shows the parameter values for the power system under investigation in this study, where R 1 and R 2 refer to the thermal power plant’s regulation parameters, B 1 and B 2 represent bias parameters of frequency, a a 12 is the rate of area capacity, and T 12 represent the coefficient of synchronization.

3. Structure of the (1+PD)-PID Cascaded Controller

The (1+PD)-PID controller, a cascaded set of (1+PD) controller and PID controller in sequence, was established and applied as a solution to the LFC issue in this study. In comparison to the complex control techniques found in the literature, the proposed controller employs the (1+PD) controller in sequence with the conventional PID controller. Hence, it has the features and benefits of the classic PID controller including strong performance, comprehensible architecture, and easy implementation in real time. Furthermore, since the two controllers are connected via serial communication, they have the propensity to reject disturbance causes immediately before they spread throughout the system. Figure 4 shows the proposed (1+PD)-PID controller structure and its parameters. Moreover, this controller, unlike other cascaded controller configurations, uses the deviation in frequency (ΔF) and the deviation in interconnected power between nearby areas (Δ P t i e ) as additional signals in addition to the corresponding area control error (ACE) [14].
The controller creates the necessary load reference signal(Δ P r e f ), which is applied as an input signal for the controlled system, obtained by operating on the ACE signals (ΔF and Δ P t i e ). For a two-area interconnected power system, the ACE is written as in Equation (5).
A C E 1 = B 1 Δ F 1 Δ P t i e A C E 2 = B 2 Δ F 2 Δ P t i e
where B1 = B2 = B is the frequency bias parameter, Δ F i is the deviation in frequency in the ith interconnected area and Δ P t i e is the deviation in the tie-line power of nearby control areas. The transfer functions in the s-domain of the proposed controller can be expressed by:
G 1 + P D ( s ) = 1 + K p 1 + K d 1 N 1 s s + N 1
G P I D ( s ) = K p 2 + K i s + K d 2 N 2 s s + N 2
where K p 1 , K d 1 , K p 2 , K i and K d 2 refer to proportional, derivative, proportional, integral, and derivative parameters, respectively. While N 1 and N 2 are the derivative filter coefficients whose values were estimated using the trial-and-error method and found to be equal to 600 for the optimal performance of the controller. Thus, the proposed controller has five degrees of freedom in its design, i.e., there are five parameters to be optimally tuned. These five parameters must all be optimized at the same time to achieve the best performance out of this controller, and AVOA is applied in this research to accomplish this task. The input to the first controller(1+PD) is the ACE and its output with the subtraction of ΔF and Δ P t i e are applied to the second controller (PID) to generate the controlling signal ( Δ P r e f ) to the system. Δ P r e f can be written mathematically in terms of the controller transfer function and controller inputs as follows:
Δ P r e f = ( { A C E × ( 1 + K p 1 + K d 1 N 1 s s + N 1 ) } Δ f Δ P t i e ) × ( K p 2 + K i s + K d 2 N 2 s s + N 2 )
In this work, the AVOA is applied to tune the (1+PD)-PID cascaded controller gains optimally. The integral absolute error (IAE) criterion, which is the integration of the deviations in areas frequencies ( Δ F 1 and Δ F 2 ) and in tie-line power ( Δ P t i e ), was used as the objective function to be minimized. Consequently, the applied objective function can be written as:
J o b j = 0 t s i m ( | Δ F 1 | + | Δ F 2 | + | Δ P t i e | ) . d t
where t s i m is the simulation time which was set to a value fair enough for the system to restore its balance, i.e., the response to settle. The fundamental interest of the current work is to optimally minimize the objective function ( J o b j ) using AVOA since doing so will result in fewer oscillations, reduced settling time, and no or little peak overshoot/ undershoot in the Δ f 1 , Δ f 2 , and Δ P t i e responses to a given system fluctuation or disturbance.
In order to guarantee that the optimization equation converges to an optimum solution, the following restrictions must be considered:
K p 1 m i n m a x K p 1 K p 1 m i n m a x K d 1 m i n m a x K d 1 K d 1 m i n m a x K p 2 m i n m a x K p 2 K p 2 m i n m a x K i m i n m a x K i K i m i n m a x K d 2 m i n m a x K d 2 K d 2 m i n m a x
where min is the minimum limit and max refers to the maximum limit of Kp1, Kd1, Kp2, Ki and Kd2, respectively. Finally, using AVOA’s effective search behavior, the optimum or nearly optimum combination of the presented controller parameters is obtained within the specified constraints for the minimal value of the objective function ( J o b j ) .

4. African Vulture Optimization Algorithm (AVOA)

The authors of [29] have introduced the African Vultures Optimization Algorithm (AVOA), a metaheuristic algorithm depending on the feeding and navigational habits of African vultures. The AVOA is also more flexible and adaptable than other metaheuristic techniques and has a low computing complexity. Correspondingly, Figure 5 provides a summary of the exploration and exploitation stages of AVOA, which are briefly explained as follows:

4.1. Exploration Phase

The probability of picking the chosen vultures to guide the others to one of the optimal solutions in every group is calculated in this phase according to the following equation:
P ( i + 1 ) = { R ( i ) | X × R ( i ) P ( i ) | × F                         i f   P 1 r a n d p 1 R ( i ) F + r 1 × ( ( U L ) × r 2 + L )   i f   P 1 < r a n d p 1
where P ( i ) is the vulture position in the current iteration and P ( i + 1 ) is its position in the next iteration, respectively. Additionally, F stands for the satiation rate of the vulture, U stands for the upper bound of agents, and L is the lower bound of them, while r 1 and r 2 are random variables, and X stands for vectors that depict the vultures’ random movement. Moreover, r a n d p 1 stands for a random number with a value generated in the middle of 0 and 1 which is generated for picking the appropriate strategy in this phase, and R ( i ) is provided by:
R ( i ) = { B e s t   V u l t u r e 1                                 i f   P i = L 1 B e s t   V u l t u r e 2                               i f   P i = L 2
where B e s t   V u l t u r e 1 is the first group optimal solution and B e s t   V u l t u r e 2 is the second group optimal one in the present iteration. Before the optimization search begins, L 1 and L 2 are initialized where their values lie between 1 and 0 and their sum equals 1.

4.2. Exploitation Phase

Based on the vulture’s satiation rate, two strategies are provided at this phase. When the satiation rate of a vulture (F) is greater than or equal to 0.5, the vultures battle for food in a circular movement, which may be estimated using Equation (13):
P ( i + 1 ) = { | X × R ( i ) P ( i ) | × ( F + r 3 ) ( R ( i )   P ( i ) )                       i f   P 2 r a n d p 2 R ( i ) ( S 1 + S 2 )                                                                                                                               i f   P 2 < r a n d p 2
where S 1 and S 2 are the spiral flight motion which may be provided using Equations (14) and (15), respectively as follows:
S 1 = R ( i ) × ( r 4 × P ( i ) 2 π ) × c o s ( P ( i ) )
S 2 = R ( i ) × ( r 5 × P ( i ) 2 π ) × s i n ( P ( i ) )
where R ( i ) is estimated in Equation (12) and r 3 , r 4 , a n d   r 5 represent random variables with a value between 0 and 1. Moreover, r a n d p 2 and r a n d p 3 are random integers with a value between 0 and 1 that are generated to obtain the appropriate strategy during the exploitation phase. Moreover, other vultures become aggressive during foraging when F < 0.5 and can be estimated using Equation (16):
P ( i + 1 ) = { A 1 + A 2 2                                                                                                                                                                         i f   P 3 r a n d p 3 R ( i ) | R ( i ) P ( i ) | × F × L e v y ( X × R ( i ) )                               i f   P 3 < r a n d p 3
where A 1 and A 2 are the motion of vultures and can be represented by Equations (17) and (18), respectively:
A 1 = B e s t   V u l t u r e 1 ( i ) B e s t   V u l t u r e 1 ( i ) × P ( i ) B e s t   V u l t u r e 1 ( i ) P ( i ) 2 × F
A 2 = B e s t   V u l t u r e 2 ( i ) B e s t   V u l t u r e 2 ( i ) × P ( i ) B e s t   V u l t u r e 2 ( i ) P ( i ) 2 × F
Furthermore, the Levy motion may be employed to improve the AVOA algorithm’s effectiveness. The controlling parameters of the AVOA used in this paper were set to 0.6, 0.4, 0.6, 0.8, and 0.2 for P1, P2, P3, L1 and L2, respectively. Finally, the AVOA algorithm has been shown to be efficient in solving various types of optimization problems [29]. The AVOA major steps can be outlined as follows:
  • Determine the iterations maximum number and the size of the population;
  • Compute the vulture fitness value;
  • Select R ( i ) using Equation (12) for all vultures;
  • Use Equation (11) to compute the position of the best vulture;
  • Depending on the vulture satiation rate, use Equation (13) or Equation (16) to update the position of the vulture;
  • Save the position of the optimal vulture then compute the value of the fitness function as long as the iteration maximum number is not reached.

5. Results and Discussions

In this section, the validity, effectiveness, and robustness of the proposed controller and the chosen optimization technique were evaluated by referring to the simulation results obtained under various scenarios. In order to do this, a MATLAB/SIMULINK model was established to simulate the proposed (1+PD)-PID controller integrated into a two-area interconnected hybrid power system containing thermal units and RES such as solar (PV) and wind energy units. The simulated model was tested by applying different cases of uncertainties such as step load perturbation, variable load profile and the penetration of renewable energy resources. In each case, the proposed controller performance was presented in comparison to different types of controllers such as classic PID, TID and FOPID controllers with their parameters obtained by the same optimization algorithm (AVOA). The frequency deviation in each area and the deviation of tie-line power between nearby areas are the system parameters to be shown for each type of controller. The results were used to assess the reliability, robustness and superiority of the presented controller as will be shown later.

5.1. Scenario 1: System Performance under 10% SLP in Area-1

5.1.1. Comparison of AVOA Performance with Other Optimization Algorithms

In this part of the study, various optimization algorithms such as Ant Lion Optimizer (ALO), Genetic Algorithm (GA), Grasshopper Optimization Algorithm (GOA), Harris Hawks Optimization (HHO) and African Vultures Optimization Algorithm (AVOA) were applied to find the optimal parameters of the proposed (1+PD)-PID controller. To do this, 10% SLP was applied to the system in Area-1 without incorporating RES into the system and each optimization algorithm was applied to the studied system to obtain the presented controller parameters and the fitness function was calculated. Figure 6 shows the response of the system in terms of frequency deviation in the two areas (ΔF1, ΔF2) and the change in the tie-line power (Δ P t i e ). Obviously, AVOA has the best performance among other optimization algorithms, and it has the lowest fitness function of 6.01 × 10−5 against 12.64 × 10−5, 26.28 × 10−5, 27.61 × 10−5, and 35.65 × 10−5 for GOA, HHO, ALO and GA, respectively.

5.1.2. Applying the Proposed AVOA Algorithm to Different Controllers

In this scenario, the performances of different types of controllers whose parameters were tuned by the optimization algorithm (AVOA) were examined. In order to achieve this, a 10% SLP in Area-1 was applied to the simulated system in the absence of RES generation. Figure 7 shows the frequency deviation in the two areas (ΔF1, ΔF2) and the change in the tie-line power (Δ P t i e ). It can be obviously noted that the proposed controller has a better response than the other controllers.
The parameters of each controller obtained by 100 iterations of the optimization algorithm (AVOA) with 30 search agents are shown in Table 2. It can be noticed that the proposed (1+PD)-PID controller offers the least minimum fitness function of 6.01 × 10−5 against 8.14 × 10−4, 15.72 × 10−4, and 19.42 × 10−4 for the TID, FOPID and PID controllers, respectively. Table 3 introduces the system performance in terms of maximum over-shoot (MO), maximum under-shoot (MU) and settling time (TS) values of the deviation of frequency in each area and tie-line power deviation for 10% SLP in area 1. The obtained simulation results clearly show that the (1+PD)-PID controller defeats the other competing controllers in the system performance specifications, especially the settling time. Settling time refers to how quickly the system can regain its equilibrium and the time it takes for the system to dampen the turbulence. For example, the settling time of the frequency deviation in Area-1 is 0.122246 S for the proposed controller versus 0.4792 S, 0.5353 S and 0.8649 S for TID, FOPID and PID controllers, respectively. Additionally, the maximum under-shoot of the frequency deviation in Area-1 is 0.0026 for (1+PD)-PID controller against 0.0051, 0.0059 and 0.0076 for TID, FOPID and PID controllers, respectively. Consequently, the proposed (1+PD)-PID controller proves its validity and efficacy in improving system stability and its superiority over the other competing controllers.

5.2. Scenario 2: The System Performance for Random Load Profile

In order to examine the robustness and effectiveness of the system, a successive random step load with different values is applied to the system. Figure 8 represents the applied random load profile. The positive value indicates an additional load to the grid, while the negative value shows an outgoing load from the grid. By adjusting the parameters of the controllers according to the values obtained in Table 1, the system responses for ΔF1, ΔF2 and Δ P t i e , due to the applied load profile, are shown in Figure 9.

5.3. Scenario 3: The Effect of Installing Solar PV Unit in Area-1

A solar PV unit was installed in the system in Area-1 in order to evaluate the system response under this type of power fluctuation due to the solar irradiance uncertainty. Figure 10 shows the system responses for ΔF1, ΔF2 and Δ P t i e using AVOA tuned (1+PD)-PID, TID, FOPID and PID controllers under the penetration of PV solar power generation. It is evident that the proposed controller offers the best performance and the least fluctuations in frequency and tie-line power in this case compared to the other controllers with which it is being compared.

5.4. Scenario 4: The Effect of Installing a Wind Farm Unit in Area-2

In this scenario, a wind farm unit was installed in Area-2 to evaluate the system’s robustness under the influence of another type of power fluctuation (wind speed variation). Figure 11 shows the frequency deviation in area 1, frequency deviation in area 2 and tie-line power under the effect of wind power generation. It can be easily noticed the superiority of the (1+PD)-PID controller over the TID, FOPID and PID controllers under the effect of wind speed fluctuations.

5.5. Scenario 5: The Effect of Inserting RES into the System

In this scenario, the system and the controller under study were subjected to hybrid renewable energy sources where a PV unit was installed in Area-1 and a wind farm unit was integrated into Area-2 at the same time. The efficiency, reliability, and robustness of the presented controller were evaluated under the impact of solar irradiance uncertainties and wind speed fluctuations simultaneously. The system responses in terms of frequency deviation in Area-1 (ΔF1) and in Area-2 (ΔF2) and tie-line power (Δ P t i e ) are presented in Figure 12 for different controllers under the circumstances of RES variation. It is clearly shown that the suggested controller has the best performance and the least variation in system frequency and interconnected tie-line power of nearby areas compared to other comparable controllers. Hence, the effectiveness, validity, reliability, and robustness of the proposed controller ((1+PD)-PID) tuned by the selected optimization algorithm (AVOA) were strictly proven in controlling hybrid power systems in terms of frequency and tie-line power in case of multi-area systems under various types of uncertainties.

6. Results and Discussions

In conclusion, this paper presented a robust cascaded (1+PD)-PID controller for providing a reliable, efficient, and simple solution to the LFC problem in power systems. The ACE was treated in the first controller (1+PD), then the output of it besides ΔF and Δ P t i e were applied as inputs to the latter controller (PID). The controller parameters were optimally obtained using AVOA depending on the IAE criterion as the objective function to ensure the best performance of the proposed controller. The validity, effectiveness and robustness of the controller were tested by subjecting the studied system to several scenarios such as SLP, random load profile and RES penetration. The superiority of the proposed controller was evaluated by comparing its performance, in terms of settling time, maximum overshoot and under-shoot, to other controllers such as PID controller, FOPID controller and TID controller which were also optimized by AVOA. Finally, the proposed controller effectively damped the tie-line power fluctuation and the frequency deviation in both areas and outperformed the other comparable controllers in all scenarios performed.
In the future, the proposed controller can be applied to interconnected systems containing more than two areas with different types of power generation units and energy storage systems. Additionally, it can be applied to real modern networks and tested under severe disturbances.

Author Contributions

Conceptualization, H.M., H.K. and K.M.A.; methodology, A.H.-E. and H.M.; software, H.M., A.S. and S.K.; validation, A.H.-E., H.K., K.M.A. and H.M.; writing—review and editing, H.M., A.S. and S.K.; supervision, A.H.-E., H.K., K.M.A. and S.K.; funding acquisition, A.S. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Science, Technology & Innovation Funding Authority (STDF) under grant (43180).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data sources employed for analysis are presented in the text.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abumeteir, H.A.; Vural, A.M. Design and Optimization of Fractional Order PID Controller to Enhance Energy Storage System Contribution for Damping Low-Frequency Oscillation in Power Systems Integrated with High Penetration of Renew-able Sources. Sustainability 2022, 14, 5095. [Google Scholar] [CrossRef]
  2. Singh, K.; Amir, M.; Ahmad, F.; Khan, M.A. An Integral Tilt Derivative Control Strategy for Frequency Control in Multimicrogrid System. IEEE Syst. J. 2021, 15, 1477–1488. [Google Scholar] [CrossRef]
  3. Singh, C.; Padhy, P.K. Fractional Order Controller Design for interconnected Power System using BAT optimization Al-gorithm. In Proceedings of the 2022 Second International Conference on Artificial Intelligence and Smart Energy (ICAIS), Coimbatore, India, 23–25 February 2022; pp. 1634–1639. [Google Scholar]
  4. Mohanty, P.; Sahu, R.K. Analysis of Hybrid Tilt Integral Derivative Controller for Multi area Power System with Real Time Simulation. In Proceedings of the 2020 International Conference on Computational Intelligence for Smart Power System and Sustainable Energy (CISPSSE), Keonjhar, Odisha, India, 29–31 July 2020. [Google Scholar]
  5. Rerkpreedapong, D.; Hasanovic, A.; Feliachi, A. Robust load frequency control using genetic algorithms and linear matrix inequalities. IEEE Trans. Power Syst. 2003, 18, 855–861. [Google Scholar] [CrossRef]
  6. Kaliannan, J.; Baskaran, A.; Dey, N.; Ashour, A.S. Ant colony optimization algorithm based PID controller for LFC of single area power system with non-linearity and boiler dynamics. World J. Model. Simul. 2016, 12, 3–14. [Google Scholar]
  7. Ali, E.; Abd-Elazim, S. Bacteria foraging optimization algorithm based load frequency controller for interconnected power system. Int. J. Electr. Power Energy Syst. 2011, 33, 633–638. [Google Scholar] [CrossRef]
  8. Mohanty, B.; Panda, S.; Hota, P. Controller parameters tuning of differential evolution algorithm and its application to load frequency control of multi-source power system. Int. J. Electr. Power Energy Syst. 2014, 54, 77–85. [Google Scholar] [CrossRef]
  9. Singh, K.; Shankar, G. PID parameters tuning using modified particle swarm optimization and its application in load frequency control. In Proceedings of the 2016 IEEE 6th International Conference on Power Systems (ICPS), New Delhi, India, 4–6 March 2016; pp. 1–6. [Google Scholar] [CrossRef]
  10. Maqbool, H.; Yousaf, A.; Asif, R.M.; Rehman, A.U.; Eldin, E.T.; Shafiq, M.; Hamam, H. An Optimized Fuzzy Based Control Solution for Frequency Oscillation Reduction in Electric Grids. Energies 2022, 15, 6981. [Google Scholar] [CrossRef]
  11. Ramlal, T.; Elgammal, A. Optimal Frequency Control Management of Grid Integration PV/Wind/FC/Storage Battery Based Smart Grid Using Multi Objective Particle Swarm Optimization MOPSO and Model Predictive Control (MPC). Eur. J. Eng. Technol. Res. 2021, 6, 50–56. [Google Scholar]
  12. Roshin, K.R.; Bindumol, E.K. Firefly Algorithm Based Tuning of Integral Controller for Frequency Regulation of Hybrid Two Area Power System with Nonlinearities and Electric Vehicles. In Proceedings of the 2020 International Conference on Power Electronics and Renewable Energy Applications (PEREA), Kannur, India, 27–28 November 2020; pp. 1–6. [Google Scholar] [CrossRef]
  13. Ali, M.; Kotb, H.; Aboras, K.M.; Abbasy, N.H. Design of Cascaded PI-Fractional Order PID Controller for Improving the Frequency Response of Hybrid Mi-crogrid System Using Gorilla Troops Optimizer. IEEE Access 2021, 9, 150715–150732. [Google Scholar] [CrossRef]
  14. Celik, E.; Öztürk, N.; Arya, Y.; Ocak, C. (1 + PD)-PID cascade controller design for performance betterment of load frequency control in diverse electric power systems. Neural Comput. Appl. 2021, 33, 15433–15456. [Google Scholar] [CrossRef]
  15. Fathy, A.; Alharbi, A.G. Recent Approach Based Movable Damped Wave Algorithm for Designing Fractional-Order PID Load Frequency Control Installed in Multi-Interconnected Plants With Renewable Energy. IEEE Access 2021, 9, 71072–71089. [Google Scholar] [CrossRef]
  16. Ahmed, E.M.; Mohamed, E.A.; Elmelegi, A.; Aly, M.; Elbaksawi, O. Optimum Modified Fractional Order Controller for Future Electric Vehicles and Renewable Energy-Based Interconnected Power Systems. IEEE Access 2021, 9, 29993–30010. [Google Scholar] [CrossRef]
  17. Ahmed, E.M.; Elmelegi, A.; Shawky, A.; Aly, M.; Alhosaini, W.; Mohamed, E.A. Frequency Regulation of Electric Vehicle-Penetrated Power System Using MPA-Tuned New Combined Fractional Order Controllers. IEEE Access 2021, 9, 107548–107565. [Google Scholar] [CrossRef]
  18. Mohamed, E.A.; Ahmed, E.M.; Elmelegi, A.; Aly, M.; Elbaksawi, O.; Mohamed, A.-A.A. An Optimized Hybrid Fractional Order Controller for Frequency Regulation in Multi-Area Power Systems. IEEE Access 2020, 8, 213899–213915. [Google Scholar] [CrossRef]
  19. Sahin, E. Design of an Optimized Fractional High Order Differential Feedback Controller for Load Frequency Control of a Multi-Area Multi-Source Power System With Nonlinearity. IEEE Access 2020, 8, 12327–12342. [Google Scholar] [CrossRef]
  20. Yakout, A.H.; Kotb, H.; Hasanien, H.M.; Aboras, K.M. Optimal Fuzzy PIDF Load Frequency Controller for Hybrid Microgrid System Using Marine Predator Algorithm. IEEE Access 2021, 9, 54220–54232. [Google Scholar] [CrossRef]
  21. Shafei, M.A.R.; Ibrahim, D.K.; Bahaa, M. Application of PSO tuned fuzzy logic controller for LFC of two-area power system with redox flow battery and PV solar park. Ain Shams Eng. J. 2022, 13, 101710. [Google Scholar] [CrossRef]
  22. Shakibjoo, A.D.; Moradzadeh, M.; Din, S.U.; Mohammadzadeh, A.; Mosavi, A.H.; Vandevelde, L. Optimized Type-2 Fuzzy Frequency Control for Multi-Area Power Systems. IEEE Access 2022, 10, 6989–7002. [Google Scholar] [CrossRef]
  23. Kader, M.O.A.; Akindeji, K.T.; Sharma, G. A Novel Solution for Solving the Frequency Regulation Problem of Renewable Interlinked Power System Using Fusion of AI. Energies 2022, 15, 3376. [Google Scholar] [CrossRef]
  24. Kumar Sahu, R.; Panda, S.; Biswal, A.; Sekhar, G.T.C. Design and analysis of tilt integral derivative controller with filter for load frequency control of mul-ti-area interconnected power systems. ISA Trans. 2016, 61, 251–264. [Google Scholar] [CrossRef]
  25. Elmelegi, A.; Mohamed, E.A.; Aly, M.; Ahmed, E.M.; Mohamed, A.-A.A.; Elbaksawi, O. Optimized Tilt Fractional Order Cooperative Controllers for Preserving Frequency Stability in Renewable Energy-Based Power Systems. IEEE Access 2021, 9, 8261–8277. [Google Scholar] [CrossRef]
  26. Ahmed, M.; Magdy, G.; Khamies, M.; Kamel, S. Modified TID controller for load frequency control of a two-area interconnected diverse-unit power system. Int. J. Electr. Power Energy Syst. 2022, 135, 107528. [Google Scholar] [CrossRef]
  27. Elkasem, A.H.A.; Khamies, M.; Hassan, M.H.; Agwa, A.M.; Kamel, S. Optimal Design of TD-TI Controller for LFC Considering Renewables Penetration by an Improved Chaos Game Optimizer. Fractal Fract. 2022, 6, 220. [Google Scholar] [CrossRef]
  28. Khokhar, B.; Dahiya, S.; Parmar, K.P.S. A Novel Hybrid Fuzzy PD-TID Controller for Load Frequency Control of a Standalone Microgrid. Arab. J. Sci. Eng. 2021, 46, 1053–1065. [Google Scholar] [CrossRef]
  29. Abdollahzadeh, B.; Gharehchopogh, F.S.; Mirjalili, S. African vultures optimization algorithm: A new nature-inspired metaheuristic algorithm for global optimization problems. Comput. Ind. Eng. 2021, 158, 107408. [Google Scholar] [CrossRef]
Figure 1. The components of the studied interconnected two-area non-reheat thermal power system.
Figure 1. The components of the studied interconnected two-area non-reheat thermal power system.
Processes 10 02320 g001
Figure 2. (a) Daily solar radiation and (b) PV output power.
Figure 2. (a) Daily solar radiation and (b) PV output power.
Processes 10 02320 g002
Figure 3. (a) The wind speed fluctuation and (b) the wind turbine out-power variation.
Figure 3. (a) The wind speed fluctuation and (b) the wind turbine out-power variation.
Processes 10 02320 g003
Figure 4. The configuration of the proposed (1+PD)-PID cascaded controller.
Figure 4. The configuration of the proposed (1+PD)-PID cascaded controller.
Processes 10 02320 g004
Figure 5. Flowchart of the AVOA algorithm.
Figure 5. Flowchart of the AVOA algorithm.
Processes 10 02320 g005
Figure 6. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power (Δ P t i e ) for 10% load change in Area-1 under the influence of (1+PD)-PID controller and different optimization algorithms.
Figure 6. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power (Δ P t i e ) for 10% load change in Area-1 under the influence of (1+PD)-PID controller and different optimization algorithms.
Processes 10 02320 g006aProcesses 10 02320 g006b
Figure 7. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power(Δ P t i e ) for 10% load change in Area-1.
Figure 7. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power(Δ P t i e ) for 10% load change in Area-1.
Processes 10 02320 g007aProcesses 10 02320 g007b
Figure 8. Random load profile.
Figure 8. Random load profile.
Processes 10 02320 g008
Figure 9. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power (Δ P t i e ) for random load profile in Area-1.
Figure 9. (a) The deviation of frequency in Area-1 (ΔF1); (b) the deviation of frequency in Area-2 (ΔF2); (c) the deviation of tie-line power (Δ P t i e ) for random load profile in Area-1.
Processes 10 02320 g009aProcesses 10 02320 g009b
Figure 10. System responses under solar PV fluctuations at Area-1: (a) ΔF1, (b) ΔF2 and (c) ΔPtie.
Figure 10. System responses under solar PV fluctuations at Area-1: (a) ΔF1, (b) ΔF2 and (c) ΔPtie.
Processes 10 02320 g010aProcesses 10 02320 g010b
Figure 11. System responses under wind power generation at Area-2: (a) ΔF1, (b) ΔF2 and (c) ΔPtie.
Figure 11. System responses under wind power generation at Area-2: (a) ΔF1, (b) ΔF2 and (c) ΔPtie.
Processes 10 02320 g011aProcesses 10 02320 g011b
Figure 12. System responses under RES insertion for: (a) ΔF1, (b) ΔF2 and (c) Δ P t i e .
Figure 12. System responses under RES insertion for: (a) ΔF1, (b) ΔF2 and (c) Δ P t i e .
Processes 10 02320 g012aProcesses 10 02320 g012b
Table 1. System parameters.
Table 1. System parameters.
ParameterValueParameterValue
K T 1 T T 0.3 s
K g 1 T g 0.03 s
K P S 120 Hz/pu.MW T P S 20 s
K P V 1 T P V 1.3 s
K W T 1 T W T 1.5 s
B 1 , B 2 0.425 pu.MW/Hz R 1 , R 2 2.4 Hz/pu.MW
a 12 −1 T 12 0.545 pu.MW/Hz
Table 2. Controller gains and fitness functions obtained by AVOA.
Table 2. Controller gains and fitness functions obtained by AVOA.
(1+PD)PIDTIDFOPIDPID
Area 1Kp11
Kd11
N1
Kp1
Ki1
Kd1
N2
127.0878
1.3142
600
85.536
149.906
0.88764
600
Kt1
N1
Ki1
Kd1
App1
375
7.805
375
53.75
9.5
Kp1
Ki1
λ1
Kd1
μ1
180
180
1
17
1.23
Kp1
Ki1
Kd1
148.58
150
25.67
Area 2Kp22
Kd22
N3
Kp2
Ki2
Kd2
N4
5.2932
0.3625
600
1.9495
2.7302
0.1812
600
Kt2
N2
Ki2
Kd2
App2
282.177
17.6448
14.027
209.509
16.574
Kp2
Ki2
λ2
Kd2
μ2
60.37
44.26
1.25
16.15
1.22
Kp2
Ki2
Kd2
139.08
147.71
59.76
Fitness Function6.01 × 10−58.14 × 10−415.72 × 10−419.42 × 10−4
Table 3. The system parameters for 10% load change.
Table 3. The system parameters for 10% load change.
ControllerΔF1ΔF2ΔPtie
MOMUTSMOMUTSMOMUTS
(1+PD)PID0.000660.00260.12461.95 × 10−51.65 × 10−50.98271.16 × 10−71.2 × 10−51.08
TID0.00190.00510.479202.36 × 10−44.9209.97 × 10−54.93
FOPID4.8 × 10−70.00590.53534.46 × 10−75.4 × 10−43.901.9 × 10−72.28 × 10−43.91
PID0.00130.00760.864905.4 × 10−43.9402.29 × 10−43.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hossam-Eldin, A.; Mostafa, H.; Kotb, H.; AboRas, K.M.; Selim, A.; Kamel, S. Improving the Frequency Response of Hybrid Microgrid under Renewable Sources’ Uncertainties Using a Robust LFC-Based African Vulture Optimization Algorithm. Processes 2022, 10, 2320. https://doi.org/10.3390/pr10112320

AMA Style

Hossam-Eldin A, Mostafa H, Kotb H, AboRas KM, Selim A, Kamel S. Improving the Frequency Response of Hybrid Microgrid under Renewable Sources’ Uncertainties Using a Robust LFC-Based African Vulture Optimization Algorithm. Processes. 2022; 10(11):2320. https://doi.org/10.3390/pr10112320

Chicago/Turabian Style

Hossam-Eldin, Ahmed, Hamada Mostafa, Hossam Kotb, Kareem M. AboRas, Ali Selim, and Salah Kamel. 2022. "Improving the Frequency Response of Hybrid Microgrid under Renewable Sources’ Uncertainties Using a Robust LFC-Based African Vulture Optimization Algorithm" Processes 10, no. 11: 2320. https://doi.org/10.3390/pr10112320

APA Style

Hossam-Eldin, A., Mostafa, H., Kotb, H., AboRas, K. M., Selim, A., & Kamel, S. (2022). Improving the Frequency Response of Hybrid Microgrid under Renewable Sources’ Uncertainties Using a Robust LFC-Based African Vulture Optimization Algorithm. Processes, 10(11), 2320. https://doi.org/10.3390/pr10112320

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