1. Introduction
Due to the expansion of industries and population growth, there has been a significant surge in the demand for electrical energy. This demand puts pressure on power systems to provide reliable and efficient electricity supply [
1]. In the management of power systems, Load Frequency Control (LFC) plays a critical role in maintaining balance between the production and consumption of electricity [
2]. In order to maintain the control of system frequency within a reasonable range, LFC regulates the output power of units to meet the power demands [
3]. LFC is a system used to prevent frequency fluctuations in electrical systems. Frequency fluctuations caused by imbalances between electricity production and consumption can lead to unwanted consequences, such as destabilizing the energy system, causing power outages, and even risking the collapse of the entire electrical grid [
4]. The importance of LFC is to bear significant importance in maintaining the stability of the energy system, preventing power outages, and even preventing the collapse of the energy system [
5].
Wind and solar power which are clean energy sources are vital for meeting global energy demands. However, their fluctuating nature poses a challenge for maintaining a stable energy system as their production relies on natural conditions that can vary throughout the day and seasons. This can cause sudden surges or drops in energy supply, which destabilizes the system [
6]. To address this, LFC systems must respond quickly to changes in energy demand and supply. Traditional energy sources are more easily controlled, while renewable sources require innovative solutions like energy storage systems such as batteries [
7].
There has been a surge in research on various types of controllers and optimization techniques aimed at improving LFC in recent times. Several classical controllers such as, Proportional Integral (PI) and Proportional Integral Derivative (PID) have been employed in LFC regulation [
8,
9]. In addition, the controllers’ effectiveness was examined in [
10] with regards to their ability to regulate LFC in a three-area thermal system that was merged with a solar thermal power plant. Although classical controllers are simple to implement, they require considerable effort adjusting controller’s parameters.
To address the limitations of conventional controllers, several researchers have utilized physics-based algorithms like Simulated Annealing (SA) [
11], Big Bang Big Crunch (BBBC) [
12], and Gravitational Search Algorithm (GSA) [
13], as well as swarm-based algorithms including Artificial Bee Colony (ABC) algorithm [
14], Ant Colony Optimization (ACO) [
15], Cuckoo Search (CS) [
16], and Particle Swarm Optimization (PSO) [
17,
18]. Evolutionary algorithms such as Genetic Algorithm (GA) [
19], Differential Evolution (DE) [
20], Firefly Algorithm (FA) [
21], Flower Pollination Algorithm (FPA) [
22], Whale Optimization Algorithm (WOA) [
23], Water Cycle Algorithm (WCA) [
24], Grasshopper Optimization Algorithm (GOA) [
25], and Grey Wolf Optimizer [
26] have also been employed for solving the frequency regulation problem. However, these meta-heuristic algorithms may suffer from certain limitations such as vulnerability to local minimum, gradual convergence, and longer computing times for highly-complex problems [
27].
On the other hand, artificial intelligence techniques such as Fuzzy Logic (FL) [
28,
29,
30,
31] and Neural Networks (NN) [
32,
33] have been utilized for solving the LFC problem owing to their ability to deal with nonlinearities in the power system. However, these algorithms suffer from various disadvantages, such as the requirement of hard work to define influential signals, determining the optimal quantity of layers and neurons for a neural network can pose a significant challenge.
In this paper, we propose a novel Sea Horse Optimizer (SHO) [
34] based tuning approach for PI/PID load frequency controllers in a two-area power system with PV and reheat thermal units. The proposed method aims to improve the performance of the LFC system in terms of frequency response to load changes, settling time, maximum overshoot and undershoot values. To evaluate the efficiency of the proposed SHO algorithm, comparisons were made with other optimization algorithms published in the literature, such as Genetic Algorithm (GA) [
35], Firefly Algorithm (FA) [
35], Whale Optimization Algorithm (WOA) [
36], and Modified Whale Optimization Algorithm (MWOA) [
36]. The results showed that the proposed SHO algorithm outperformed the other optimization methods. Therefore, the SHO algorithm can be considered as a potential alternative for optimizing PI/PID LFC systems.
The main contributions of this research work are as follows:
1. Development of a novel Sea Horse Optimizer algorithm for PI/PID load frequency controller tuning in a two-area power system with PV and reheat thermal units, which has not been previously proposed in the literature.
2. Application of the proposed SHO-based tuning approach to effectively optimize PI and PID controllers for LFC problem.
3. Demonstration of the effectiveness of the proposed SHO algorithm in improving the performance of the LFC system in terms of frequency response, settling time, and overshoot/undershoot values.
4. Identification of the potential of the SHO algorithm as an alternative method for optimizing PI/PID LFC systems, due to its ability to handle non-linear and complex systems more effectively, and provide better overall system performance.
The rest of this paper is organized as follows:
Section 2 presents the system modeling for the two-area power system with PV and reheat thermal units.
Section 3 describes the proposed Sea Horse Optimizer algorithm for tuning the PI/PID load frequency controllers. In
Section 4, the optimization problem is formulated, and the proposed SHO-based tuning approach is applied to optimize the PI and PID controllers.
Section 5 presents the simulation results and discussion on different conditions. Finally,
Section 6 provides the conclusion part and suggests further work.
2. System Modeling
The system model is presented as a two-area test system shown in
Figure 1 [
35,
36]. The model consists of two areas: Area 1 has a photovoltaic (PV) system with maximum power point tracking (MPPT), while Area 2 has a reheat thermal unit. The MATLAB/Simulink program was used to model the PV-thermal testing system. The components shown in
Figure 1, including the PV system, Governor, Turbine, Reheater, and Power System, are modeled using their respective transfer functions. The mathematical models of these components are approximated as first-order linear transfer function models. Equation (1) to (5) represent the mathematical frameworks of the system components. The utilization of the MPPT algorithm is aimed at enhancing the effectiveness of the PV system, given that the performance of a solar panel is contingent on several factors, including solar radiation intensity and ambient temperature.
The PV system’s transfer function, which includes the photovoltaic panel, MPPT, inverter, and filter, is represented by Equation (1) [
37]:
where the PV system’s transfer function is
GPV. The governor’s transfer function is given by following Equation (2).
where
GGVR is the transfer function of the governor,
is the governor’s gain parameter,
is the governor’s time constants. The turbine’s transfer function is given by following Equation (3).
where the turbine’s transfer function is
GTRB,
is the turbine’s gain parameter, the turbine’s time constant is
. The reheater’s transfer function is given by following Equation (4).
where the reheater’s transfer function is
GRHT,
is the reheater’s gain parameter, the reheater’s time constant is
. The power system’s transfer function is given by following Equation (5).
where the power system’s transfer function is
GPS,
is the power system’s gain parameter,
is the power system’s time constants.
3. Sea Horse Optimizer Algorithm
A new meta-heuristic method, the Sea Horse Optimizer (SHO) algorithm, was proposed by Zhao in 2022 [
34]. The SHO algorithm is an meta-heuristic algorithm based on population that imitates the social behavior of sea horses. The algorithm comprised in three main components: movement, hunting, and reproduction. To achieve a balance between exploration and exploitation capabilities of the SHO algorithm, local and global search abilities are designed for the movement and hunting behaviors, respectively, while the reproductive behavior complements them.
The SHO algorithm initializes by creating a population of candidate solutions.
where
represents the number of dimension in the search space while
indicates the size of the population used in the SHO algorithm. Each member of the population of seahorses corresponds to a possible solution within the problem’s search space. In a minimization optimization problem, the individual with the lowest fitness value is considered as the elite individual and is denoted by
.
can be obtained using Equation (7).
where
denotes the value of the cost function for a particular problem, which is used to evaluate the fitness of candidate solutions in the search space.
Sea horses motion behavior is composed of two states: Brownian motion and Levy flight. The brownian motion of seahorses enables better exploration in the search space, while Levy flight simulates the movement step size of seahorses, allowing them to migrate allowing them to move to different locations early on to avoid excessive local exploitation. To compute the updated position of a seahorse in iteration
t, we can formulate these two scenarios as follows:
where
is defined by Lévy flight distribution function with a randomly generated parameter
from the interval [0, 2]. The spiral movement component of SHO is represented by the coordinates
. The constant coefficient
is used to control the step size of the Lévy flight,
is Brownian motion’s random walk coefficient. The normal random number
is used to introduce stochasticity in the Brownian motion component.
The hunting behavior of seahorses can result in either success or failure. Successful hunting occurs when a seahorse catches its prey by moving faster than it, while failure means the seahorse explores the search space further. The mathematical notation representing this hunting behavior can be expressed as:
where, the new location of the sea horse after hunting at iteration
t is denoted as
,
is the randomly generater number within [0, 1],
is a directly decreasing parameter that adjusts seahorse-based step length during the hunting process.
The reproductive behavior of seahorses is split into male and female groups based on population fitness values, with male seahorses responsible for reproduction.
where,
fathers and
mothers refer to the male and female populations, respectively, while
denotes all
arranged in ascending order of their correspoding fitness values. The algorithm selects half of the best-fit individuals from the population to create a new generation of candidate solutions. The expression of the
ith offspring is as follows:
where
is the random number between [0, 1],
and
individuals chosen at random from the male and female populations.
The SHO algorithm is specifically designed to tackle optimization problems that involve continuous search spaces, and has demonstrated favorable outcomes in various applications.
Figure 2 shows the flowchart of the proposed SHO algorithm.
The SHO algorithm provides a novel approach to solving optimization problems, and its effectiveness and efficiency make it a promising technique for various applications. This study proposes the utilization of the SHO algorithm for optimizing PI/PID controllers in the LFC problem.
4. Sea Horse Optimization Based Load Frequency Control
When there is an mismatch between power production and demand in a power system, the system frequency may change. This frequency change leads to the formation of a signal called Area Control Error (ACE). The ACE signal is a critical input for the LFC controller, which is responsible for maintaining the equilibrium between the total system load and the total system generation. By monitoring the ACE signal, the LFC controller can adjust the power output of the generators to guarantee that the power system remains stable and reliable. Mathematically, the ACE signals for Area-1 and Area-2 can be expressed as:
where,
and
are frequency bias parameters of the both Areas 1 and 2, respectively.
and
are the frequency deviations for the both Areas 1 and 2, respectively.
is the tie-line power variation. The LFC controller can apply various strategies to manage power generation using the ACE signal. For example, the controller can increase or decrease generator outputs based on whether the ACE signal is positive or negative, which can help to mitigate power imbalances between supply and demand in the control area. Moreover, the controller is capable of regulating generator outputs based on the rate of change of the ACE signal, which is instrumental in preserving the stability of the power system.
PI or PID controllers are control strategies that are often used to reset the ACE signal. These controllers measure the ACE signal and generate the necessary control signals to keep it close to zero.
4.1. PI Controller
The PI controller is a type of feedback control that is commonly used in LFC applications. The PI controller uses two control terms to adapt the system frequency, which are the proportional term and the integral term. The proportional term acts to stabilize the frequency quickly, while the integral term acts to eliminate the steady-state error by adjusting the control signal over time.
The transfer function of the PI controller is formulated as:
where,
is the proportional gain and
is the integral gain.
4.2. PID Controller
The PID controller is a widely adopted feedback control strategy that encompasses three key elements: proportional, integral, and derivative. It serves as an extension of the PI controller and is developed to enhance the performance of the control system, particularly with regards to stability and transient response.
The transfer function of the PID controller is expressed as:
where,
,
and
are the proportional, integral and derivative gains, respectively.
The controller gains in this study will be optimized through the proposed SHO algorithm, with the objective function being the integral of time multiply absolute error (ITAE) of the frequency deviations in both areas and tie-line power variation.
The ITAE objective function is expressed in Equation (16) and is utilized to determine the optimal controller gains, ultimately improving the power system’s performance. The ITAE is a performance metric. This criterion helps optimize the controller for fast response and stable operation. ITAE calculates the integral of multiplying the controller by the error in the ACE signal at a given time and is used to evaluate the controller performance.
where
J is the objective function of the LFC problem. It is aimed to minimize the
J.
When optimizing PID controllers, it is necessary to constrain the gain parameters
,
and
within predefined lower and upper bounds. These bounds represent variable constraints that ensure the stability and performance of the controller. Equation (17) specifies the limit ranges for the PID controller.
The proposed SHO algorithm aims to optimize the gain parameters of a PI/PID controller by utilizing the
J objective function in Equation (16) and considering the vari-able constraints in Equation (17).
Figure 3 illustrates how the SHO algorithm optimizes PI/PID parameters for the LFC problem. The proposed SHO algorithm was employed to the LFC problem using the test system shown in
Figure 1, and the parameter values of the test model are provided in
Table 1 [
35,
36].
Table 2 presents the parameters of the proposed SHO algorithm, while
Figure 3 illustrates the block diagram of the SHO-based PI/PID controller for the LFC problem.
In the proposed SHO algorithm, the number of sea horse is set to 50, the maximum iteration number is 100, and the gain parameter range of the controller is set between -2 and 2. These parameter values were chosen to facilitate comparison with similar algorithm parameters used in other published studies, such as FA, GA, WOA, and MWOA [
35,
36]. The speed of a selected sea horse (u) is used when calculating a new position (v). The learning factor (l) of a selected sea horse is used when updating its position. These parameters can be adjusted to improve the effectiveness of the SHO algorithm and to reach optimal solutions faster.
Figure 1.
The two-area power system.
Figure 1.
The two-area power system.
Figure 2.
The flowchart of the proposed SHO algorithm.
Figure 2.
The flowchart of the proposed SHO algorithm.
Figure 3.
The proposed SHO algorithm optimizing PI/PID parameters for LFC problem.
Figure 3.
The proposed SHO algorithm optimizing PI/PID parameters for LFC problem.
Figure 4.
Frequency response of Area-1 to load changes in both regions with PI controllers.
Figure 4.
Frequency response of Area-1 to load changes in both regions with PI controllers.
Figure 5.
Frequency response of Area-2 to load changes in both regions with PI controllers.
Figure 5.
Frequency response of Area-2 to load changes in both regions with PI controllers.
Figure 6.
Power response in tie-line to load changes in both regions with PI controllers.
Figure 6.
Power response in tie-line to load changes in both regions with PI controllers.
Figure 7.
Converge curve of the proposed SHO algorithm for LFC problem.
Figure 7.
Converge curve of the proposed SHO algorithm for LFC problem.
Figure 8.
Frequency response of Area-1 to load changes in both regions with PID controllers.
Figure 8.
Frequency response of Area-1 to load changes in both regions with PID controllers.
Figure 9.
Frequency response of Area-2 to load changes in both regions with PID controllers.
Figure 9.
Frequency response of Area-2 to load changes in both regions with PID controllers.
Figure 10.
Power response in tie-line to load changes in both regions with PID controllers.
Figure 10.
Power response in tie-line to load changes in both regions with PID controllers.
Figure 11.
Frequency response of Area-1 to large load disturbances.
Figure 11.
Frequency response of Area-1 to large load disturbances.
Figure 12.
Frequency response of Area-2 to large load disturbances.
Figure 12.
Frequency response of Area-2 to large load disturbances.
Figure 13.
Power response in tie-line to large load disturbances.
Figure 13.
Power response in tie-line to large load disturbances.
Figure 14.
Change of solar radiation.
Figure 14.
Change of solar radiation.
Figure 15.
Response of LFC to change in solar radiation.
Figure 15.
Response of LFC to change in solar radiation.
Table 1.
Parameters of the test system.
Table 1.
Parameters of the test system.
|
PV system gain 1 |
18 |
|
PV system gain 2 |
900 |
|
PV system time constant 1 |
100 |
|
PV system time constant 2 |
50 |
|
Governor gain |
1 p.u. MW |
|
Governor time constant |
0.08 sec |
|
Turbine gain |
1 p.u. MW |
|
Turbine time constant |
0.3 sec |
|
Reheat gain |
0.33 p.u. MW |
|
Reheat time constant |
10 sec |
|
Power system gain of thermal area |
120 Hz/p.u. MW |
|
Power system time constant |
20 sec |
|
Regulation droop |
0.4 Hz/p.u. MW |
|
Frequency bias constant |
0.8 p.u. |
|
Tie-line power coefficient |
0.545 |
Table 2.
Parameters of SHO algorithm for solving LFC problem.
Table 2.
Parameters of SHO algorithm for solving LFC problem.
Parameter |
Value |
Number of sea horse |
50 |
Iteration number |
100 |
|
0.05 |
|
0.05 |
|
0.05 |
Lower bound for [; ; ] |
[-2; -2; -2] |
Upper bound for [; ; ] |
[2; 2; 2] |
Table 3.
Controller parameters for the PV-thermal power system.
Table 3.
Controller parameters for the PV-thermal power system.
Parameter |
Methods |
SHO-tuned PID (proposed) |
MWOA-tuned PID [36] |
SHO-tuned PI (proposed) |
WOA-tuned PI [36] |
FA-tuned PI [35] |
GA-tuned PI [35] |
|
-0.8599 |
-0.1070 |
-0.67012 |
-0.4563 |
-0.8811 |
-0.5663 |
|
-0.1290 |
-0.0906 |
-0.5371 |
-0.2254 |
-0.5765 |
-0.4024 |
|
-1.9396 |
-0.6112 |
- |
- |
- |
- |
|
-2.0000 |
-1.8938 |
-2.0000 |
-0.8967 |
-0.7626 |
-0.5127 |
|
-2.0000 |
-1.8935 |
-0.8476 |
-0.9865 |
-0.8307 |
-0.7256 |
|
-0.2614 |
-0.2505 |
- |
- |
- |
- |
ITAE |
0.8582 |
1.5602 |
2.5308 |
4.1211 |
7.4259 |
12.1244 |
Table 4.
Specifications of dynamic responses with optimized PI controller.
Table 4.
Specifications of dynamic responses with optimized PI controller.
Parameters |
GA |
FA |
WOA |
SHO (proposed) |
∆f1 |
Overshoot (M+) |
0.1638 |
0.1577 |
0.07997 |
0.001733 |
Undershoot (M-) |
-0.2966 |
-0.3154 |
-0.2015 |
-0.04374 |
Settling Time (s) |
26.73 |
26.44 |
26.30 |
12.62 |
∆f2 |
Overshoot (M+) |
0.1571 |
0.1228 |
0.09816 |
0.1012 |
Undershoot (M-) |
-0.2435 |
-0.2295 |
-0.2216 |
-0.1807 |
Settling Time (s) |
23.64 |
23.60 |
25.54 |
15.75 |
∆Ptie |
Overshoot (M+) |
0.05636 |
0.04643 |
0.0534 |
0.03823 |
Undershoot (M-) |
-0.04921 |
-0.04778 |
-0.03814 |
-0.03215 |
Settling Time (s) |
27.73 |
26.45 |
21.07 |
18.04 |
Table 5.
Performance indisces values of compared algorithms.
Table 5.
Performance indisces values of compared algorithms.
Techniques |
IAE |
ITAE |
ISE |
ITSE |
GA-tuned PI |
2.3341 |
12.1244 |
0.3202 |
0.8618 |
FA-tuned PI |
1.7207 |
7.4259 |
0.2907 |
0.4723 |
WOA-tuned PI |
1.0566 |
4.1211 |
0.1663 |
0.4262 |
SHO-tuned PI |
0.6491 |
2.5308 |
0.1021 |
0.26179 |
MWOA-tuned PID |
0.5625 |
1.5602 |
0.0815 |
0.0601 |
SHO-tuned PID |
0.3091 |
0.8582 |
0.0448 |
0.0369 |