Next Article in Journal
Synchronization of Incommensurate Fractional-Order Chaotic Systems Based on Linear Feedback Control
Next Article in Special Issue
Advances in Optimization and Nonlinear Analysis
Previous Article in Journal
Mitochondrial DNA Profiling by Fractal Lacunarity to Characterize the Senescent Phenotype as Normal Aging or Pathological Aging
Previous Article in Special Issue
Non-Dominated Sorting Manta Ray Foraging Optimization for Multi-Objective Optimal Power Flow with Wind/Solar/Small- Hydro Energy Sources
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Design of TD-TI Controller for LFC Considering Renewables Penetration by an Improved Chaos Game Optimizer

1
Department of Electrical Engineering, Faculty of Engineering, Aswan University, Aswan 81542, Egypt
2
Department of Electrical Engineering, College of Engineering, Northern Border University, Arar 1321, Saudi Arabia
3
Prince Faisal bin Khalid bin Sultan Research Chair in Renewable Energy Studies and Applications (PFCRE), Northern Border University, Arar 1321, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(4), 220; https://doi.org/10.3390/fractalfract6040220
Submission received: 8 February 2022 / Revised: 30 March 2022 / Accepted: 6 April 2022 / Published: 13 April 2022
(This article belongs to the Special Issue Advances in Optimization and Nonlinear Analysis)

Abstract

:
This study presents an innovative strategy for load frequency control (LFC) using a combination structure of tilt-derivative and tilt-integral gains to form a TD-TI controller. Furthermore, a new improved optimization technique, namely the quantum chaos game optimizer (QCGO) is applied to tune the gains of the proposed combination TD-TI controller in two-area interconnected hybrid power systems, while the effectiveness of the proposed QCGO is validated via a comparison of its performance with the traditional CGO and other optimizers when considering 23 bench functions. Correspondingly, the effectiveness of the proposed controller is validated by comparing its performance with other controllers, such as the proportional-integral-derivative (PID) controller based on different optimizers, the tilt-integral-derivative (TID) controller based on a CGO algorithm, and the TID controller based on a QCGO algorithm, where the effectiveness of the proposed TD-TI controller based on the QCGO algorithm is ensured using different load patterns (i.e., step load perturbation (SLP), series SLP, and random load variation (RLV)). Furthermore, the challenges of renewable energy penetration and communication time delay are considered to test the robustness of the proposed controller in achieving more system stability. In addition, the integration of electric vehicles as dispersed energy storage units in both areas has been considered to test their effectiveness in achieving power grid stability. The simulation results elucidate that the proposed TD-TI controller based on the QCGO controller can achieve more system stability under the different aforementioned challenges.

1. Introduction

Recently, the world has become voracious in utilizing electrical power due to the growth of industrial and residential loads. Therefore, it was necessary to establish new electrical power grids to accommodate the load demands. As a result, energy planners were directed to penetrate the renewable energy sources (RESs) with the traditional power grids in the electrical power system to reduce the demerits of these traditional units. In addition, the penetration of RESs with newly established power systems is considered to have an economically good and positive rate that saves in the utilization of the oil, coal, and gas that operate traditional power plants, whereas the resulting flames from burning oil and coal lead to the release of carbon dioxide gas, causing an increase in the ozone hole and an increase in the global warming phenomenon [1]. Although the presence of RESs in electrical power grids reduces the severity of the resulting pollution from the traditional units, these renewable sources suffer from a lack of system inertia. As a result of the reduction in power system inertia caused by renewable sources, the stability and security of the system (i.e., more fluctuations in system frequency) will be affected [2,3]. Moreover, several reasons lead to more frequency fluctuations, such as a mismatching between the generated power and the demand power, system parameter variations, and different sorts of load variations. Hence, the fluctuations in system frequency can be tackled by the LFC [4]. Researchers have done their best to develop several control techniques for achieving reliability in power systems by attaining system frequency and tie-line power flow within tolerable limits.
Many interests have been prompted by researchers to address the issue of LFC in different structures of the power system; (i.e., the single-area power system [5,6], the multi-area interconnected power system [7,8,9,10], and the deregulated power system [11,12]). In addition, several different control techniques have been implemented to overcome the system frequency fluctuations, such as the intelligent control techniques (i.e., fuzzy logic controllers [13], artificial neural networks [14], and adaptive neuro-fuzzy controllers [15]). Moreover, several robust control techniques have been utilized to enhance the power system performance, such as the H-infinite technique [16] and µ-synthesis [17]. Furthermore, optimal control techniques, such as the linear quadratic Gaussian [18] and linear quadratic regulator [19], are implemented to attain the frequency within tolerable limits. In this regard, the majority of the industrial control loop is the proportional-integral-derivative (PID), due to its reputable merits (i.e., simpleness in construction, applicability, functionality, comfort, and inexpensiveness) [20]. Even so, it suffers from a bulky, complicated process when selecting its parameters using trial and error methods. Thus, researchers have been striving to accomplish the optimal PID controller, according to the different optimization techniques utilized in getting the optimal controller parameters. This design of the optimal PID controller leads to ensuring a reliable system performance in comparison to the conventional PID controller when facing the uncertainties in a studied power grid. Accordingly, several optimization techniques have been utilized to fine-tune the optimal PID controller parameters meticulously, including the grasshopper optimization algorithm [21], the ant colony optimization technique [22], the Jaya algorithm [23], and the class topper optimization algorithm [24].
On the other side, the fractional order controllers (FOCs) have become a distinct candidate in power system stabilizing due to their merits (i.e., flexibility in configuration and a higher degree of freedom). The FOCs have several types of poles, such as the hyper-damped poles, that need to be fine-tuned. Accordingly, this leads to an expansion in the stable region, giving more flexibility in the controller design process [25]. Furthermore, there are several types of controllers belonging to the FOC family; the fractional –order-proportional-integral-derivative (FOPID) is one member of this family that has been presented in [26,27]. The FOPID controller has been utilized in several electrical power systems [28,29]. Moreover, the TID controller represents one of the FOCs; it looks exactly like the PID controller in construction except for one difference, which is that the proportional parameter is tilted with a ( 1 / s 1 n ) transfer function. This additional transfer function provides the optimization process with better feedback and good tracking performance. Lately, the TID controller has been implemented for solving the LFC problem due to its good merits (i.e., it can change the parameters of the closed-loop system; it has a tremendous ability in disturbance rejection; and it has more reliability with robustness) [30,31]. There is no doubt that fractional calculus provides several options to researchers for creativity and diversity in controller designing. As a result, different engineering problems have been solved by utilizing the amalgamation of the FOPID and TID properties as a hybrid controller [32]. In addition, the researchers’ minds are destined to implement another strategy in control design, which is the cascaded controllers (CCs) form that includes one controller followed by another one; the CCs have more tuning knobs that give better results than in the utilization of non-crude CCs. Thus, many scientific studies have been presented using the different CCs to solve the LFC problem [33,34]. Another construction has been applied while designing different controllers for studying the LFC issue, which depends on the combination of two different controllers to take the benefits of both controllers. There are examples of the combination of different proposed controllers from literature, such as the combination of the model predictive control (MPC) controller with the linear quadratic Gaussian controller [35] and the combination of an adaptive MPC with the recursive polynomial model estimator [36]. Furthermore, a new controller structure, labeled as a feed-forward/feed-backward controller, has been presented to reduce the disadvantages of the PID and TID controllers during system uncertainties that affect the input of the control signal. Thus, many studies have been presented to elucidate the robustness of the feed-forward/feed-backward controller structure in achieving system stability. The integral-proportional-derivative (I-PD) controller and the integral-tilt-derivative (I-TD) controller have been proposed to cope with the LFC problem, achieving more system stability compared to the PID and TID controllers, respectively [37,38].
The achievement of system stability is not dependent on the controller design only, but the utilized optimization technique represents a critical issue that must be selected carefully to attain the optimal controller parameters. Previously, the traditional optimization methods such as the tracking approach [39] and the aggregation methods [40] were applied for regulating the system frequency. In fact, the traditional optimization methods suffer from several drawbacks, such as slump, deathtrap in local minimums, the need for more iterations, and dependence on their initial conditions to attain the optimal solution. So, meta-heuristic optimization techniques such as the artificial bee colony [41], salp swarm algorithm (SSA) [42], and whale optimization algorithm (WOA) [43] have been proposed to overcome all of the previous drawbacks. Though the meta-heuristic optimization algorithms are not usually guaranteed to find the optimal global solution, they can often find a sufficiently good solution in a reasonable time. So, they are an alternative to exhaustive search, which would take exponential time. Moreover, these techniques have several demerits, such as slowing in the rate of convergence, poor local search capability, and local optimum convergence. In this regard, algorithmic scientists have improved these techniques to diminish all of their previous drawbacks. Examples of improved algorithms utilized to achieve system stability are presented as the improved stochastic fractal search algorithm [44] and the sine augmented scaled sine cosine [45]. In this regard, the authors in this work proposed an improved algorithm known as QCGO to select the suggested combining TD-TI controller parameters to attain the optimal studied power grid performance.
Referring to the aforementioned literature related to the LFC issue, there are several control strategies that depend on the designer experience, such as the MPC, the H-infinite techniques, and the fuzzy logic control, that can attain the desired performance, but their parameter-selecting strategies take a long time. In addition, the conventional PID controller has some difficulties when facing system uncertainties. Moreover, several studies have been presented utilizing conventional algorithms and meta-heuristic optimization techniques that have many demerits in comparison to the improved techniques that develop the searching process and obtain the global solution with a few search agents. Furthermore, several previous studies did not consider the different challenges that face power systems (i.e., different types of load variations such as series SLP and RLV, the high penetration of RESs, and communication time delay). According to the above salient observations, this study proposed a new control construction labeled as a combining TD-TI controller that is derived from the form of a TID controller to enhance the studied system stability. The parameters of the proposed combining TD-TI controller can be selected by utilizing the improved algorithm QCGO when considering the challenges of high RESs penetration, different load perturbation types, and communication time delay.
The studied work in this paper is presented to overcome the limitations of the previously published works in the literature. Table 1 elucidates the differences between this work and the other published works related to the LFC issue.
The main contributions of this work can be elucidated in detail as follows:
i.
The proposal of a control structure combining TD-TI controllers for LFC of the hybrid two-area interconnected power systems.
ii.
The proposal of a novel technique known as QCGO via improving the quantum mechanics of the CGO algorithm based on the particle swarm optimizer (PSO) to improve the exploration and exploitation strategies of the main CGO algorithm.
iii.
The application of the improved CGO to select the optimal parameters of the proposed controller structure.
iv.
The validation of the performance of the proposed algorithm through a fair-maiden comparison between the proposed QCGO algorithm and other previous techniques (i.e., Supply-demand-based optimization (SDO), WOA, butterfly optimization algorithm (BOA), and the conventional CGO), based on applying 23 bench functions, as well as a fair comparison between the proposed algorithm and other previous algorithms (i.e., CGO, SSA), considering the proposed controller in the multi-area power grid for frequency stability analysis.
v.
The consideration of several challenges, such as the high RESs penetration in both areas, different load perturbation types, and communication time delay to study the system stability state.
vi.
The comparison of the performance of the proposed control TD-TI structure based on QCGO with other available controllers, such as the PID-based teaching learning-based optimization (TLBO) [46]; the PID-based arithmetic optimization algorithm (AOA) [47]; the proposed TD-TI control structure based on CGO; the proposed TD-TI control structure based on SSA; the TID controller based on CGO; and the TID controller based on QCGO, is presented to ensure the effectiveness and robustness of the proposed control structure based on the QCGO algorithm in achieving more system reliability and stability.
vii.
The consideration of the integration of electrical vehicles (EVs) in both areas to support the proposed controller in overcoming the system frequency excursions during high renewables penetration.
The remainder of this article is organized into several sections that are clarified as follows: the studied system topology which considers the high penetration of RESs and EVs is illustrated in Section 2. Section 3 discusses the proposed control approach and the formulation of the studied problem. Then, the procedure of the improved QCGO technique is given in Section 4. Moreover, the simulation results according to the different scenarios are clarified in Section 5. Finally, Section 6 summarizes the conclusions of the current work.

2. The Studied System Topology

2.1. Two-Area Interconnected Hybrid Power Grid Configuration

In this article, the issue of LFC related to electrical power grids has been addressed by conducting a study on two-area interconnected hybrid power systems. The studied power grid encompasses two interconnected areas, which include several conventional generation power plants, such as the thermal unit, hydropower unit, and gas unit. The capacity of each area in the studied power grid that includes the three traditional units (i.e., thermal, hydro, and gas) is 2000 MW of rated power [48], of which the largest percentage of electrical power sharing went to the thermal power plant, which contributes 1087 MW, then the hydropower plant, which contributes 653 MW, and the gas turbine, sharing the generated power with 262 MW. The investigated power grid is presented as a simplified model shown in Figure 1.
Figure 2 shows the block diagram of the studied two-area interconnected hybrid power grid. The transfer functions in the studied power grid are listed in Table 2. The amalgamation of the TD-TI controller is proposed to be equipped in both areas for each generation unit to minimize the oscillations in the frequencies of both areas and the tie-line power flow between them. The attitude of the input signal of the proposed combining TD-TI controller can be represented as the ACE, while the attitude of the output signal can be represented as the action of the secondary/supplementary control on each generation power plant, in order to obtain extra active power for enhancing the power grid performance. Table 3 elucidates all the parameters included in the studied power grid with their nominal values. The ACEs in both areas can be obtained according to the formulas that follow in Equations (1) and (2) [47]:
A C E 1 = Δ C tie 1 2 + B 1 Δ i 1
A C E 2 = Δ C tie 2 1 + B 2 Δ i 2

2.2. The Installation of Wind Farm Model

This work presents the high penetration of RESs, including wind power in the investigated hybrid power grid. The professional software MATLAB/SIMULINK program (R2015a) (The MathWorks, Inc., Natick, MA, USA) is used in implementing the simplified model of wind power in order to share its energy in the first area of the studied power grid. The aforementioned wind power model generates power in the same way as the real behavior of the generated power from real wind farms. This is achieved using a white-noise block that is utilized in getting a random speed, which is multiplied by the wind speed, as shown in Figure 3 [47]. The captured output power from the wind model can be formulated in the following equations [47].
P wt = 1 2 ρ A T v w 3 C p λ , β
C p λ , β = C 1 C 2 λ i C 3 β C 4 β 2 C 5 × e C 6 λ i + C 7 λ T
λ T =   λ T OP = ω T r T V W
1 λ i = 1 λ T + 0.08 β 0.035 β 3 + 1
All of these mentioned parameter values for the utilized wind farm are presented in [47]. Figure 4 shows the random output power of 257 wind turbine units of 750 KW for each wind power unit. The value of the generated power from the studied wind farm is about 192 MW.

2.3. The Installation of the PV Model

The Photovoltaic (PV) model can be built by utilizing the professional software MATLAB/SIMULINK program (R2015a) described in Figure 5. The generated output power from the model is similar to the real generated output power from a real PV plant. In addition, the output energy of the PV model is penetrated in the second area of the studied power grid at about 116 MW. Here, the white-noise block in the MATLAB program (R2015a) is used for obtaining random output oscillations that are multiplied by the standard output power generated from a real PV plant. The generated energy from the presented PV model can be obtained as formulated in Equation (7) [6]. Figure 6 clarifies the random output power generated from the PV model.
Δ P solar   = 0.6 × P solar  

2.4. The Installation of EV Model

EVs can participate in frequency regulation effectively due to the receiving of the LFC order and pass this signal to the EV to control the power during the charging and discharging process. Moreover, the response of the LFC signal can be limited through the availability of the numbers of controllable EVs in the studied power grid and by the state of the charge related to their capacity, whereas the model of the EV is similar to the model of the battery energy storage system, due to the included batteries that supply extra energy to the power grid during fluctuations for regulating the frequency excursions. However, the batteries in EVs may not be in full charging capacity due to the nature of EVs being of mobility and load, which affects the amount of extra energy to tackle the LFC problem. Thus, it is important to check the level of the EV charging to ensure more system enhancement under different system fluctuations. The output power from an EV can be obtained by the first-order transfer function, including the electrical vehicle time constant T E V , which equals 0.28 s in series with the electrical vehicle controllers’ gain, K E V , which equals 1, where K E V is represented as the ratio of the exchange in charging power of the EV’s batteries to the change of system frequency. The transfer function that represents the EV model is formulated in Equation (8) [49]. Figure 7 describes the EV model that was built in the MATLAB/SIMULINK program (R2015a).
K E V   1 + s   T E V

3. Control Methodology and Problem Formulation

Due to the high RES penetration, communication time delay, and various types of load perturbations, it is essential to implement a robust controller to enhance the system performance during abnormal conditions. Hence, this study proposes a newly developed controller construction known as a combining TD-TI controller to overcome any fluctuations resulting from the previous considerations/challenges. Moreover, the proposed controller parameters have been selected based on an improved algorithm labeled as QCGO.

3.1. The Proposed Control Strategy

This paper presents an efficient controller labeled as the combining TD-TI controller, which represents an improved modified structure of the TID controller that is shown in Figure 8. The TID controller is a sort of fractional order controller (FOC) that depends on the fractional-order calculus in its design. The TID controller construction is similar to the PID controller construction except for one difference, which is that the proportional parameter is tilted with a ( 1 / s 1 n ) transfer function. In this regard, this paper proposed a combining TD-TI controller, as derived from TID controller, due to the merits of the TID, such as the ability to tune easily, superior fluctuations rejection, and better sensitivity due to variations of the system parametric [50]. The proposed combining TD-TI controller is utilized to enhance the studied power grid performance, such as by damping frequency oscillations in both areas and overcoming fluctuations related to the tie-line power flow. Furthermore, the proposed combining TD-TI controller parameters are selected utilizing an improved QCGO algorithm. In general, the transfer function of the combining TD-TI controller is formulated as follows [50]:
G i 1 , TD s = K t i S 1 n + K d i S
G i 2 , TI s = K t i S 1 n + K i i s
G i , total s = G i 1 , TD s +   G i 2 , TI s
where i refers to the specified proposed controller of the (thermal, hydro, and gas) turbines; thus, ( i = 1, 2, 3). The gain values ( K t i ,   K i i , and K d i ) are selected within the range of [0, 10], and n is tuned in the range of [1, 10]. The control signal of the i t h area can be expressed as follows [38]:
U i s = G i , total ( s )   ×   A C E i ( s )
According to the process of controller designing, there are several sorts of performance criteria, such as the integral time absolute error (ITAE), the integral of squared error (ISE), the integral time squared error (ITSE), and the integral of absolute error (IAE). The criteria of ITAE and ISE are often utilized in the literature for minimizing the objective function due to their merits in comparison to ITSE and IAE, whereas the strategy of the ISE criteria in minimizing the objective function is the integrating of the square of error signal over simulation time. For ease, the ISE criteria can effectively dampen the large errors compared to the small errors as the square of the large errors is larger than the square of the small errors. It can be said that the ISE criteria can penalize the large errors with tolerance for the presence of continuous small errors along with time simulation. Thus, the authors of this work do not hesitate in putting in the ITAE criteria utilized in minimizing the objective function because of the multiplication of the time term by the integral of the absolute error. The multiplied time term in ITAE criteria makes the optimization process more fast which achieves more system stability than utilizing the ISE criteria [51]. The ITAE criteria can be formulated as follows [47]:
J = I T A E = 0 T s i m t .   [ Δ f 1   +   Δ f 2     +   Δ P t i e   ]   d t
where d t is represented as the time interval for taking the error signals’ samples over the simulation process.

3.2. The Proposed Optimization Technique

In this subsection, the CGO method is briefly described; then, the process of the QCGO technique is presented.

3.2.1. Chaos Game Optimization (CGO) Algorithm

This algorithm is based on certain rules of the chaos theory, where the arrangement of fractals is by the chaos game idea. Firstly, an initialization procedure is configured by determining the initial positions of the solution candidates from the following equations [52]:
X = [ X 1 X 2 . . X i . . X n ] = [ x 1 1   x 1 2   . .   x 1 j . . x 1 d x 2 1   x 2 2   . .   x 2 j . . x 2 d . . . . . . . . x i 1   x i 2   . .   x i j . . x i d . . . . . . . . x n 1   x n 2   . .   x n j . . x n d ] , i = 1 , 2 , . , m j = 1 , 2 , . . , d
x i j 0 = x i , min j + rand . x i , max j x i , min j , i = 1 , 2 , . , m j = 1 , 2 , . . , d
where d denotes the dimension of the problem and m refers to the total number of initialized candidates inside the search space. x i , min j , x i , max j are the lower and upper bounds of the decision variables. The position updating process for the temporary triangles is presented in Figure 9. The mathematical representation of the seed i 1 , as shown in Figure 9a, is as follows [52]:
seed i 1 = X i + α i × β i GB γ i × MG i , i = 1 , 2 , . , m
where GB is the global best, α i represents the movement limitation factor, and β i and γ i denote vectors randomly created by numbers in the range of [0, 1]. MG i is the mean group. From Figure 9b, seed i 2 can be calculated as follows [52]:
seed i 2 = GB + α i × β i × X i γ i × MG i , i = 1 , 2 , . , m
While seed i 3 , which is displayed in Figure 9c, is mathematically computed as below [52]:
seed i 3 = MG i + α i × β i × X i γ i × GB , i = 1 , 2 , . , m
Finally, seed i 4 , which is shown in Figure 9d, can be mathematically represented as follows [52]:
seed i 4 = X i x i k = x i k + R ,   k = 1 , 2 , . , d
where R refers to a vector with random numbers in the range of [0, 1].

3.2.2. The Proposed Quantum Chaos Game Optimization (QCGO) Algorithm

In this subsection, quantum mechanics is used to develop the original CGO algorithm. This quantum model of a CGO algorithm is called here QCGO algorithm. Quantum mechanics was employed to develop the PSO in [54]. In the quantum model, by employing the Monte Carlo method, the solution x new 4 is calculated from this equation [54]:
If h ≥ 0.5
x new 1 = p + α · Mbest i X i   · ln 1 / u
Else
x new 1 = p α · Mbest i X i   · ln 1 / u
End.
where α refers to a design parameter, u and h denote uniform probability distribution in the range [0, 1], and Mbest is the mean best of the population and is defined as the mean of the GB positions. It can be calculated as follows [54]:
Mbest = 1 N l = 1 N p g , l i
where g is the index of the best solution among all the solutions.

4. The Procedure of the Improved QCGO Algorithm

The Performance of QCGO

The proposed QCGO algorithm competency and performance are evaluated on the numerous benchmark functions, using the statistical measurements, such as best values, mean values, median values, worst values, and standard deviation (STD), for the best solutions achieved using the proposed technique and the other well-known algorithms. The results attained by the QCGO technique are compared with three recent meta-heuristic techniques, including SDO [55], WOA [56], and BOA [57], in addition to the conventional CGO. All of the mentioned techniques were executed for the maximum number of iterations of the function of 200 and a population size of 50 for 20 independent runs, using Matlab R2016a working on Windows 8.1, 64 bit (Microsoft, Albuquerque, NM, USA). All computations were performed on a Core i5-4210U CPU@ 2.40 GHz of speed (Intel Corporation, Santa Clara, CA, USA) and 8 GB of RAM. Figure 10 shows the qualitative metrics on F1, F2, F3, F5, F6, F8, F10, F12, F15, F18, and F22, with 2D views of the functions, convergence curve, average fitness history, and search history.
Table 4, Table 5 and Table 6 show the statistical results of the proposed QCGO technique and other algorithms when applied for the three types of benchmark functions (unimodal, multimodal, and composite, respectively). The best-obtained values using the QCGO, CGO, SDO, WOA, and BOA algorithms are displayed in bold. It is clearly seen that the QCGO algorithm achieves the optimal solution for most of those benchmark functions. The convergence curves of these techniques for those functions are illustrated in Figure 11, and the boxplots for each algorithm for these functions are displayed in Figure 12. From those figures, it is seen that the QCGO technique reached a stable point for all functions, and the boxplots of the proposed QCGO technique are very narrow and stable for most functions compared to the other techniques.

5. Simulation Results and Discussions

In this study, the proposed control strategy is implemented in the secondary control loop with the high integration of RESs, considering different load variation types to restore the studied system frequency at the pre-defined value, where the presented control strategy relied on the combining TD-TI controller, which is optimally designed by an improved QCGO algorithm to obtain the minimum value of the frequency fluctuations for the studied power grid. Moreover, the performance of the suggested control strategy is compared with other control strategies (i.e., TID and PID). All of the simulation results for the studied two-area, multi-unit power grid are implemented using the professional software MATLAB/SIMULINK® program (R2015a) to ensure the efficacy of the proposed controller in enhancing the studied system performance. The code of the proposed QCGO algorithm is an m-file linked to the studied model for the optimization process. The simulation results are performed on a PC with Intel Core i5-2.60 GHz with 4.00 GB of RAM. The frequency stability has been assessed by applying different operating conditions through the following scenarios.
  • Scenario A: evaluation of the studied power grid performance considering various load variation types (i.e., SLP, series SLP, and RLV).
  • Scenario B: evaluation of the studied power grid performance considering high penetration of RESs in both areas with series SLP and RLV.
  • Scenario C: evaluation of the studied power grid performance considering communication time delay.
  • Scenario D: evaluation of the studied power grid performance considering EV integration in both areas.
The studied power grid performance can be evaluated by measuring the value of the best objective function that is represented by the ITAE value over the iterations. For most, several initial considerations must be addressed while optimizing the proposed TD-TI controller using the proposed improved QCGO algorithm, such as the search agent number that equals 30 and the total iterations/attempts that equal 100. The convergence curve that is shown in Figure 13 clarifies the performance of the proposed combining TD-TI controller based on QCGO compared to the combining TD-TI controller based on CGO and SSA and compared with the TID controller based on QCGO and CGO. The demonstrated convergence curve can be obtained considering a 1% SLP at 10 s in the first area of the studied power grid, without any RESs penetration in both areas. It is clear that the proposed combining TD-TI controller based on QCGO attained the lowest value of the objective function compared to the other mentioned controllers that relied on various optimization techniques. As a result, the convergence curve elucidates the effectiveness of the proposed QCGO algorithm. It can be seen that the curve behavior of the proposed TD-TI based on QCGO starts with a 0.1098 objective function value; then, this value drops along the iterations to end up at the final iteration with a 0.0729 objective function value, whereas the behavior of the proposed controller/proposed algorithm can be described as it reaches the best objective function value quickly compared to the other utilized controllers via different techniques. Moreover, it can be said, the rest curve behaviors are far from the optimum goal achieved by the suggested controller using QCGO, demonstrating its robustness in damping the oscillations effectively.
Scenario A: evaluation of the studied system performance considering different load variation types (i.e., SLP, series SLP, and random load).
This scenario included a fair-maiden comparison between the proposed combining TD-TI controller utilizing the QCGO algorithm and the other published controllers, such as the PID controller based on TLBO and AOA. Moreover, the proposed combining TD-TI controller based on the improved QCGO technique was compared with different mentioned controllers, such as the TID controller based on QCGO and CGO and the combining TD-TI controller based on CGO and SSA, to test the stability of the studied power grid performance.
Case A.1: The SLP was selected as a challenge by applying it in the first area of the studied power grid to test the efficacy of the proposed combining TD-TI controller in enhancing the system performance. The applicable SLP occurred at 10 s with a 1% value, whereas the SLP can occur in the electrical power grids through disconnecting some generators from all the generation stations that may lead to blackouts with the shutdown of all the stations’ generators. In addition, SLP may be represented as an unexpected switch of the connected electrical loads that may lead to instability in the system performance by increasing the wear and tear on the generators in the power grid.
Case A.1.1: This case presents a comparison between the performance of the proposed combining TD-TI controller in this work and the other published performances of the PID controller, to prove the efficacy of the proposed controller in attaining the main target (damping frequency oscillations). Table 7 indicates all of the aforementioned controller parameters that are utilized in diminishing the fluctuations in the system frequency and power flow in the tie line. In addition, Figure 14 clarifies a comparison between the different dynamic studied system responses (i.e., Δ f 1 , Δ f 2 , and Δ P t i e ) of the proposed combining TD-TI controller, using QCGO and the PID controller based on TLBO and AOA, and considering a 1% SLP in the first area.
Table 8 illustrates the various specifications of the system performance, such as overshoot ( O s h ), undershoot ( U s h ), and the objective function values related to fluctuations in both the area frequencies and the power flow within the tie line. Table 8 clarifies the superiority of the proposed combining TD-TI controller-based, improved QCGO algorithm to achieve stability in the studied power grid. For ease, Table 9 denotes the percentage improvements in U s h and O s h for combining TD-TI/QCGO and PID/AOA, based on the PID/TLBO.
As can be seen, the improved QCGO algorithm utilized in fine-tuning the proposed combining TD-TI controller obtains the optimal controller parameters, which leads to attaining the optimal solution with a 0.075 objective function value. The obtained objective function value related to the proposed controller using an improved QCGO algorithm is the best compared to those attained from the published PID controller based on TLBO and AOA, which equal 0.402 and 0.189, respectively. It can be seen that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/QCGO is 60.01% and 52.43%, respectively. In contrast, the percentage improvement in U s h and O s h of Δ f 1 related to PID/AOA is 42.11% and 32.70%, respectively.
Case A.1.2: This case presents a suggestion of utilizing the TID controller based on CGO and QCGO to compare it with the proposed combining TD-TI controller based on QCGO to test the robustness of the proposed one in regulating the studied system frequency. All of the previously mentioned controller parameters are presented in Table 7. Moreover, Figure 15 describes a fair comparison between all of the dynamic system responses related to the proposed combining TD-TI controller based on QCGO and all those responses of the TID controller based on CGO and QCGO.
Table 10 illustrates the different specifications of the system performance, such as O s h , U s h , and the objective function values related to excursions in both the area frequencies and the power flow within the tie line. Table 10 clarifies the superiority of the proposed combining TD-TI controller-based improved QCGO algorithm in achieving system reliability. In addition, Table 11 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and TID/(CGO, QCGO) based on the PID/TLBO.
Table 10 clarifies that the obtained objective function value related to the proposed controller using an improved QCGO algorithm that equals 0.075 is the best compared to those attained from the TID controller based on CGO and QCGO, which equal 0.1381 and 0.1351, respectively. Moreover, Table 11 denotes that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ f 2 related to combining TD-TI/QCGO is 86.4% and 99.36%, respectively. In contrast, the percentage improvement in U s h and O s h of Δ f 2 related to TID/QCGO is 73.04% and 25.35%, respectively.
Case A.1.3: This case presents the SSA algorithm as a meta-heuristic optimization technique to tune the proposed combining TD-TI controller and make a comparison between it and the CGO and QCGO techniques in selecting the optimal controller parameters to prove that the improved QCGO algorithm can achieve more system stability compared to utilizing the different mentioned algorithms. Table 7 presents the aforementioned controller parameters that were utilized in overcoming the LFC problem in the studied power grid. Moreover, Figure 16 describes a fair comparison between all of the dynamic system responses related to the proposed combining TD-TI controller based on QCGO and all those responses of the combining TD-TI controller based on SSA and CGO.
Table 12 illustrates the different specifications of the system performance, such as O s h , U s h , and the objective function values related to the oscillations in both the area frequencies and the power flow within the tie line. Table 12 clarifies the superiority of the proposed combining TD-TI controller-based, improved QCGO algorithm in achieving system reliability. In addition, Table 13 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and combining TD-TI/(CGO, SSA), based on the PID/TLBO.
Table 12 clarifies that the obtained objective function value related to the suggested controller using an improved QCGO algorithm that equals 0.075 is the best compared to those attained from the combining TD-TI controller based on CGO and SSA, which equal 0.078 and 0.087, respectively. Moreover, Table 13 denotes that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ P t i e related to combining TD-TI/QCGO is 82.6% and 99.12%, respectively. In contrast, the percentage improvement in U s h and O s h of Δ P t i e related to combining TD-TI/SSA is 76.85% and 92.76%, respectively.
Case A.2: In this case, the performance of the proposed combining TD-TI controller optimized with an improved QCGO algorithm has been tested and assessed by subjecting a series SLP in the first area of the studied power grid. The series SLP is represented as an emulation of the series changing in the realistic connected loads. It can be said that the series SLP is considered as a series-forced switch of generators or series interrupts of the connected loads. Figure 17 describes the applied form of the series SLP. In addition, the different dynamic system responses are indicated in Figure 18 to elucidate the superiority of the suggested combining TD-TI controller based on QCGO compared to those of the other controllers optimized with different algorithms (i.e., combining TD-TI based on CGO and SSA) in the presence of the series SLP in the first area.
Table 14 illustrates the values of O s h and U s h related to the different system dynamic responses (i.e., Δ f 1 , Δ f 2 , and Δ P t i e ) according to oscillations in both the area frequencies and the power flow within the tie line. Table 14 clarifies the superiority of the proposed combining TD-TI controller-based improved QCGO algorithm in achieving system stability. In addition, Table 15 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and combining TD-TI/CGO based on the combining TD-TI/SSA.
Table 14 clarifies that the suggested controller using an improved QCGO algorithm achieves more system stability after looking at the obtained results of the O s h and U s h values. Moreover, Table 15 denotes that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/QCGO is 26.13% and 25.71%, respectively. In contrast, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/CGO is 15.81% and 14.29%, respectively.
Case A.3: In this case, the studied power grid has been subjected to RLVs in the first area. The RLVs are a diverse combination of series perturbations in industrial connected loads to the grid that cause the same effects on the grid (i.e., unbalance in electrical power grid and the occurrence of blackout). The RLV is formed in Figure 19. In addition, Figure 20 describes the different dynamic power system responses explaining the efficacy of the proposed combining TD-TI controller based on QCGO in achieving more of a reduction in the system frequency fluctuations and the power flow in the tie line compared to the other ones.
Table 16 illustrates the values of O s h and U s h related to the different system dynamic responses (i.e., Δ f 1 , Δ f 2 , and Δ P t i e ) according to the oscillations in both the area frequencies and the power flow within the tie line. Table 16 presents the robustness of the proposed combining TD-TI controller-based improved QCGO algorithm in achieving system stability. In addition, Table 17 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and combining TD-TI/CGO based on the combining TD-TI/SSA.
Table 16 clarifies that the proposed controller via an improved QCGO algorithm achieves more system stability after looking at the obtained results of the O s h and U s h values. Additionally, Table 17 denotes that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/QCGO is 20.67% and 26.00%, respectively. However, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/CGO is 10.00% and 16.00%, respectively.
Scenario B: evaluation of the studied system performance considering high penetration of RESs in both areas with series SLP and RLV.
Another challenge of high penetrating of RESs (i.e., wind energy in the first area and PV energy in the second area) is addressed in this study to test the robustness of the proposed combining TD-TI controller in reducing the studied system fluctuations. The series SLP and RLV are applied in the first area as well as integration of the RESs in the power grid. The penetration of RESs represents a burden on the studied power grid due to their demerits (i.e., lack of system inertia).
Case B.1: robustness test for the proposed combining TD-TI controller optimized by improved QCGO considering high RES penetration as well as series SLP challenge.
This section clarifies the dynamic system performance of the investigated power grid, taking into consideration a series SLP, high penetration of wind energy at t = 100 s in the first area and PV at t = 200 s in the second area. These mentioned challenges have been presented to ensure the reliability and effectiveness of the proposed combining TD-TI controller based on an improved QCGO algorithm in enhancing the studied power grid performance. Figure 21 clarifies the applicable series SLP form in the first area. Moreover, all the dynamic power grid responses represented in Δ f 1 , Δ f 2 and Δ P t i e are shown in Figure 22.
Table 18 illustrates the values of O s h and U s h related to the aforementioned system dynamic responses due to deviations in both the area frequencies and the power flow within the tie line. Table 18 presents the robustness of the proposed combining TD-TI controller-based improved QCGO algorithm in achieving system reliability. In addition, Table 19 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and combining TD-TI/CGO based on the combining TD-TI/SSA.
It can be summarized that Table 18 clarifies that the proposed controller/proposed algorithm achieves more system stability after showing the obtained results of the O s h and U s h values. In this regard, Table 19 clarifies that the proposed combining TD-TI controller-based QCGO achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ P t i e related to combining TD-TI/QCGO is 30.41% and 26.15%, respectively. However, the percentage improvement in U s h and O s h of Δ P t i e related to combining TD-TI/CGO is 7.22% and 6.15%, respectively.
Case B.2: robustness test for the proposed combining TD-TI controller optimized by improved QCGO considering high RES penetration as well as RLV.
This section includes a robustness test by the penetrating of RESs at both areas of the studied power grid with the applicable RLV in the first area. This test summarized the superiority of the proposed combining TD-TI controller based on an improved QCGO algorithm in overcoming the frequency excursions for the studied power grid. The applicable RLV is shown in Figure 23. Moreover, the behavior of both the area frequencies and the power flow in the tie line is clarified in Figure 24.
Table 20 elucidates the values of O s h and U s h related to all the mentioned system dynamic responses due to the deviations in both the area frequencies and the power flow within the tie line. Table 20 proves the robustness of the proposed controller/proposed algorithm in achieving system reliability. In addition, Table 21 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO and combining TD-TI/CGO based on the combining TD-TI/SSA.
It can be said that Table 20 clarifies that the proposed controller/proposed algorithm achieves more system stability after knowing the obtained results of the O s h and U s h values. In this regard, Table 21 clarifies that the proposed controller/proposed algorithm achieves a higher percentage in improving all of the system dynamic performance. For example, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/QCGO is 37.50% and 22.58%, respectively. However, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/CGO is 25.00% and 12.9%, respectively.
Scenario C: evaluation of the studied system performance considering communication time delay, high penetration of RESs in both areas, and RLV.
This scenario presents the suggestion of the communication time delay challenge that is applied before and after the control action with a 0.01 s time delay value and also considers the applicable random load with high RES penetration to test the robustness of the suggested combining TD-TI controller in system stabilizing. The RLV behavior is described in Figure 25. Moreover, the different dynamic system responses represented in Δ f 1 , Δ f 2 , and Δ P t i e are shown in Figure 26.
Figure 26 summarizes and elucidates the effectiveness of the proposed controller via the proposed technique in achieving system stability and reliability after testing the effect of the time delay in the controller action and in receiving the error signal. The proposed QCGO/combining TD-TI scheme shows excellent results in overcoming all the challenges and gaining more system stability.
Scenario D: evaluation of the studied system performance, considering the effect of EV integration, high penetration of RESs in both areas, and RLV.
This scenario presents the integration of EVs in both areas of the studied power grid to test the effectiveness of EVs in regulating the studied system frequency and the power flow between both areas. Figure 27 shows the applicable RLV in the first area. Figure 28 illustrates the charging/discharging power of both the EVs that are integrated into both areas of the studied power grid. Moreover, the various dynamic system responses represented in Δ f 1 , Δ f 2 and Δ P t i e are described in Figure 29.
Table 22 presents the values of O s h and U s h related to all the different mentioned system dynamic responses due to the deviations in the both area frequencies and the power flow within the tie line. Table 22 proves that the proposed controller/proposed algorithm considering EV penetration in the studied system achieves more system stability compared to not utilizing these EVs. In addition, Table 23 clarifies the percentage improvements in U s h and O s h for combining TD-TI/QCGO with and without penetration of the EVs based on the combining TD-TI/SSA.
It can be observed that Table 22 clarifies that the proposed controller/proposed algorithm achieves more system stability after presenting the values of the obtained O s h and U s h . In this regard, Table 23 clarifies that the proposed controller/proposed algorithm achieves a higher percentage in improving all system dynamic performance, whereas the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/QCGO considering EV penetration is 47.50% and 33.23%, respectively. In contrast, the percentage improvement in U s h and O s h of Δ f 1 related to combining TD-TI/CGO without EV penetration is 37.50% and 22.58%, respectively. In brief, the integration of EVs in the studied power grid can aid in dampening the frequency fluctuations due to their energy storage power which feeds the system with the extra power at abnormal conditions to obtain all the system dynamic responses within the tolerable limits.

6. Conclusions

This paper includes main points that are clarified as mentioned below:
  • A new control structure was proposed based on the TID controller labeled as a combining TD-TI controller for frequency stabilizing in the power grid.
  • A multi-area interconnected hybrid power system that includes several traditional units (i.e., thermal, hydro, and gas) has been presented in this work to test the efficacy of the combining TD-TI controller.
  • An improved algorithm was proposed named QCGO to develop the searching strategy of the main CGO algorithm to attain the optimum solution.
  • Twenty-three bench functions were applied to prove the effectiveness of the improved QCGO algorithm compared to other different techniques (i.e., SDO, WOA, BOA, and the conventional CGO).
  • The robustness of the QCGO-TD-TI controller has been validated by a fair comparison between its performance and other performances of TD-TI controllers based on the algorithms from the literature (i.e., SSA, TLBO, and AOA).
  • The CGO-TD-TI controller performance was compared with the QCGO-TD-TI controller to ensure that the improved QCGO algorithm attains more optimal results than the main CGO algorithm.
  • The efficacy of the suggested combining TD-TI controller has been ensured through a fair-maiden comparison between its performance and the performances of other mentioned controllers (i.e., TID and PID).
  • Several scenarios have been presented in this work to study the effectiveness of the suggested controller in tackling the problem of LFC, such as applying different load variation types, the high penetration of RESs in both areas, and applying the communication time delay.
  • EV integration was proposed in both areas to test its performance in enhancing the studied power grid frequency.
  • All previous simulation results have confirmed the ability of the proposed combining TD-TI controller to effectively handle the LFC problem. Moreover, the improved QCGO algorithm proved its robustness in selecting the optimal controller parameters, which led to achieving more system stability.

Author Contributions

Conceptualization, A.H.A.E., M.K., M.H.H. and S.K.; data curation, A.M.A.; formal analysis, A.H.A.E., M.K. and M.H.H.; funding acquisition, A.M.A. and S.K.; investigation, A.H.A.E., M.K. and M.H.H.; methodology, A.M.A. and S.K.; project administration, A.H.A.E., M.K. and M.H.H.; resources, A.M.A. and S.K.; supervision, S.K. and A.M.A.; validation, A.H.A.E., M.K. and M.H.H.; visualization, A.H.A.E., M.K. and M.H.H.; writing-original draft, A.H.A.E., M.K. and M.H.H.; writing-review and editing, A.M.A. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia through the project number “IF_2020_NBU_416”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number “IF_2020_NBU_416”. The authors gratefully thank the Prince Faisal bin Khalid bin Sultan Research Chair in Renewable Energy Studies and Applications (PFCRE) at Northern Border University for their support and assistance.

Conflicts of Interest

The authors declare that there is no conflict of interest.

Nomenclature

SymbolsParameters
SLPStep load perturbation
RLVRandom load variation
TIDTilt-Integral-Derivative
TI-TDCombining Tilt-Integral Tilt-Derivative
PIDProportional-Integral-Derivative
FOCsFractional-Order Controllers
FOPIDFractional-Order PID
CCsCascaded Controllers
MPCModel predictive control
I-PDIntegral-Proportional Derivative
I-TDIntegral-Tilt Derivative
PSOParticle swarm optimization
SDOSupply-demand-based optimization
WOAWhale optimization algorithm
AOAArithmetic optimization algorithm
TLBOTeaching learning-based optimization
SSASalp swarm algorithm
BOAButterfly optimization algorithm
CGOChaos game optimization
QCGOImproved chaos game optimization
LFCLoad frequency control
ACEArea control error
p.uPer unit
i t h Subscript refers to the specified area
EVsElectrical vehicles
RESsRenewable energy sources
O s h overshoot
U s h undershoot
P wt Wind turbine output power
ρThe air density
A T The area swept by the blades of a turbine
V W The wind speed
C p The coefficient of the rotor blades
C 1 - C 7 The turbine coefficients
βThe pitch angle
r T The radius of the rotor
ω T The rotor speed
λ T The optimum tip-speed ratio
λ i The intermittent tip-speed ratio
B 1 Frequency bias factor of Area 1
B 2 Frequency bias factor of Area 2
Δ f 1 Frequency deviation in area 1
Δ f 2 Frequency deviation in area 2
Δ P t i e 1 2 Tie-line power flow from area 1 to area 2
Δ P t i e 2 1 Tie-line power flow from area 2 to area 1
T 12 Coefficient of synchronizing
R 1 Regulation constant of thermal turbine
R 2 Regulation constant of hydropower plant
R 3 Regulation constant of gas turbine
a 12 Control area capacity ratio
K T Participation factor for thermal unit
K H Participation factor for hydro unit
K G Participation factor for a gas unit
K p s Gain constant of power system
T p s The time constant of the power system
T s g Governor time constant
T t Turbine time constant
K r Gain of reheater steam turbine
T r Time constant of reheater steam turbine
T g h Speed governor time constant of hydro turbine
T r s Speed governor reset time of the hydro turbine
T r h The transient droop time constant of hydro turbine speed governor
T w Nominal string time of water in penstock
b g Gas turbine constant of valve positioner
c g Valve positioner of gas turbine
Y c The lag time constant of the gas turbine speed governor
X c The lead time constant of the gas turbine speed governor
T c r Gas turbine combustion reaction time delay
T f c Gas turbine fuel time constant
T c d Gas turbine compressor discharge volume–time constant
K E V Gain of electrical vehicle
T E V The time constant of electrical vehicle
ITAEIntegral time absolute error
ISEIntegral square error
IAEIntegral absolute error
ITSEIntegral time squared error
K t The tilted gain
K i The integral gain
K d The derivative gain
n The tilt fractional component n   0
K p The proportional gain
d t The time interval for taking error signals’ samples
T s i m Total time of simulation process
J The objective function

References

  1. Ellabban, O.; Abu-Rub, H.; Blaabjerg, F. Renewable energy resources: Current status, future prospects and their enabling technology. Renew. Sustain. Energy Rev. 2014, 39, 748–764. [Google Scholar] [CrossRef]
  2. Abazari, A.; Soleymani, M.M.; Babaei, M.; Ghafouri, M.; Monsef, H.; Beheshti, M.T.H. High penetrated renewable energy sources-based AOMPC for microgrid’s frequency regulation during weather changes, time-varying parameters and generation unit collapse. IET Gener. Transm. Distrib. 2020, 14, 5164–5182. [Google Scholar] [CrossRef]
  3. Hamouda, N.; Babes, B.; Kahla, S.; Soufi, Y.; Petzoldt, J.; Ellinger, T. Predictive Control of a Grid Connected PV System Incorporating Active Power Filter functionalities. In Proceedings of the 2019 1st International Conference on Sustainable Renewable Energy Systems and Applications (ICSRESA), Tebessa, Algeria, 4–5 December 2019; pp. 1–6. [Google Scholar]
  4. Balu, N.J.; Lauby, M.G.; Kundur, P. Power System Stability and Control; Electrical Power Research Institute, McGraw-Hill Professional: Washington, DC, USA, 1994. [Google Scholar]
  5. Khamies, M.; Magdy, G.; Ebeed, M.; Kamel, S. A robust PID controller based on linear quadratic gaussian approach for improving frequency stability of power systems considering renewables. ISA Trans. 2021, 117, 118–138. [Google Scholar] [CrossRef] [PubMed]
  6. Khamari, D.; Kumbhakar, B.; Patra, S.; Laxmi, D.A.; Panigrahi, S. Load Frequency Control of a Single Area Power System using Firefly Algorithm. Int. J. Eng. Res. 2020, 9. [Google Scholar] [CrossRef]
  7. Chen, M.-R.; Zeng, G.-Q.; Xie, X.-Q. Population extremal optimization-based extended distributed model predictive load frequency control of multi-area interconnected power systems. J. Frankl. Inst. 2018, 355, 8266–8295. [Google Scholar] [CrossRef]
  8. Jagatheesan, K.; Anand, B.; Samanta, S.; Dey, N.; Santhi, V.; Ashour, A.S.; Balas, V.E. Application of flower pollination algorithm in load frequency control of multi-area interconnected power system with nonlinearity. Neural Comput. Appl. 2016, 28, 475–488. [Google Scholar] [CrossRef]
  9. Guha, D.; Roy, P.K.; Banerjee, S. Application of backtracking search algorithm in load frequency control of multi-area interconnected power system. Ain Shams Eng. J. 2018, 9, 257–276. [Google Scholar] [CrossRef] [Green Version]
  10. Guha, D.; Roy, P.K.; Banerjee, S. Load frequency control of interconnected power system using grey wolf optimization. Swarm Evol. Comput. 2016, 27, 97–115. [Google Scholar] [CrossRef]
  11. Tasnin, W.; Saikia, L.C.; Raju, M. Deregulated AGC of multi-area system incorporating dish-Stirling solar thermal and geothermal power plants using fractional order cascade controller. Int. J. Electr. Power Energy Syst. 2018, 101, 60–74. [Google Scholar] [CrossRef]
  12. Sharma, M.; Dhundhara, S.; Arya, Y.; Prakash, S. Frequency excursion mitigation strategy using a novel COA optimised fuzzy controller in wind integrated power systems. IET Renew. Power Gener. 2020, 14, 4071–4085. [Google Scholar] [CrossRef]
  13. Chen, G.; Li, Z.; Zhang, Z.; Li, S. An Improved ACO Algorithm Optimized Fuzzy PID Controller for Load Frequency Control in Multi Area Interconnected Power Systems. IEEE Access 2020, 8, 6429–6447. [Google Scholar] [CrossRef]
  14. Akula, S.K.; Salehfar, H. Frequency Control in Microgrid Communities Using Neural Networks. In Proceedings of the 2019 North American Power Symposium (NAPS), Wichita, KS, USA, 13–15 October 2019; pp. 1–6. [Google Scholar]
  15. Yousef, H. Adaptive fuzzy logic load frequency control of multi-area power system. Int. J. Electr. Power Energy Syst. 2015, 68, 384–395. [Google Scholar] [CrossRef]
  16. Zhang, H.; Liu, J.; Xu, S. H-Infinity Load Frequency Control of Networked Power Systems via an Event-Triggered Scheme. IEEE Trans. Ind. Electron. 2020, 67, 7104–7113. [Google Scholar] [CrossRef]
  17. Bevrani, H.; Feizi, M.R.; Ataee, S. Robust Frequency Control in an Islanded Microgrid: H∞ and μ-Synthesis Approaches. IEEE Trans. Smart Grid 2015, 7, 706–717. [Google Scholar] [CrossRef] [Green Version]
  18. Rahman, M.; Sarkar, S.K.; Das, S.K.; Miao, Y. A comparative study of LQR, LQG, and integral LQG controller for frequency control of interconnected smart grid. In Proceedings of the 2017 3rd International Conference on Electrical Information and Communication Technology (EICT), Khulna, Bangladesh, 7–9 December 2017. [Google Scholar]
  19. Das, S.K.; Rahman, M.; Paul, S.K.; Armin, M.; Roy, P.N.; Paul, N. High-Performance Robust Controller Design of Plug-In Hybrid Electric Vehicle for Frequency Regulation of Smart Grid Using Linear Matrix Inequality Approach. IEEE Access 2019, 7, 116911–116924. [Google Scholar] [CrossRef]
  20. Singh, V.P.; Kishor, N.; Samuel, P. Improved load frequency control of power system using LMI based PID approach. J. Frankl. Inst. 2017, 354, 6805–6830. [Google Scholar] [CrossRef]
  21. Lal, D.K.; Barisal, A.K.; Tripathy, M. Load Frequency Control of Multi Area Interconnected Microgrid Power System using Grasshopper Optimization Algorithm Optimized Fuzzy PID Controller. In Proceedings of the 2018 Recent Advances on Engineering, Technology and Computational Sciences (RAETCS), Allahabad, India, 6–8 February 2018. [Google Scholar]
  22. Dhanasekaran, B.; Siddhan, S.; Kaliannan, J. Ant colony optimization technique tuned controller for frequency regulation of single area nuclear power generating system. Microprocess. Microsyst. 2020, 73, 102953. [Google Scholar] [CrossRef]
  23. Annamraju, A.; Nandiraju, S. Coordinated control of conventional power sources and PHEVs using jaya algorithm optimized PID controller for frequency control of a renewable penetrated power system. Prot. Control. Mod. Power Syst. 2019, 4, 28. [Google Scholar] [CrossRef] [Green Version]
  24. Rai, A.; Das, D.K. Optimal PID Controller Design by Enhanced Class Topper Optimization Algorithm for Load Frequency Control of Interconnected Power Systems. Smart Sci. 2020, 8, 125–151. [Google Scholar] [CrossRef]
  25. Tepljakov, A.; Gonzalez, E.A.; Petlenkov, E.; Belikov, J.; Monje, C.A.; Petráš, I. Incorporation of fractional-order dynamics into an existing PI/PID DC motor control loop. ISA Trans. 2016, 60, 262–273. [Google Scholar] [CrossRef]
  26. Podlubny, I.; Dorcak, L.; Kostial, I. On fractional derivatives, fractional-order dynamic systems and PI/sup λ/D/sup μ/-controllers. In Proceedings of the 36th IEEE Conference on Decision and Control, San Diego, CA, USA, 12 December 1997; Volume 5, pp. 4985–4990. [Google Scholar]
  27. Podlubny, I. Fractional-order systems and PI/sup/spl lambda//D/sup/spl mu//-controllers. IEEE Trans. Autom. Control. 1999, 44, 208–214. [Google Scholar] [CrossRef]
  28. Morsali, J.; Zare, K.; Hagh, M.T. Comparative performance evaluation of fractional order controllers in LFC of two-area diverse-unit power system with considering GDB and GRC effects. J. Electr. Syst. Inf. Technol. 2018, 5, 708–722. [Google Scholar] [CrossRef]
  29. Gheisarnejad, M.; Khooban, M.H. Design an optimal fuzzy fractional proportional integral derivative controller with derivative filter for load frequency control in power systems. Trans. Inst. Meas. Control. 2019, 41, 2563–2581. [Google Scholar] [CrossRef]
  30. Topno, P.N.; Chanana, S. Differential evolution algorithm based tilt integral derivative control for LFC problem of an interconnected hydro-thermal power system. J. Vib. Control. 2017, 24, 3952–3973. [Google Scholar] [CrossRef]
  31. 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]
  32. 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]
  33. Saha, A.; Saikia, L.C. Load frequency control of a wind-thermal-split shaft gas turbine-based restructured power system integrating FACTS and energy storage devices. Int. Trans. Electr. Energy Syst. 2018, 29, e2756. [Google Scholar] [CrossRef]
  34. Prakash, A.; Murali, S.; Shankar, R.; Bhushan, R. HVDC tie-link modeling for restructured AGC using a novel fractional order cascade controller. Electr. Power Syst. Res. 2019, 170, 244–258. [Google Scholar] [CrossRef]
  35. Mohamed, T.H.; Shabib, G.; Abdelhameed, E.H.; Khamies, M.; Qudaih, Y. Load Frequency Control in Single Area System Using Model Predictive Control and Linear Quadratic Gaussian Techniques. Int. J. Electr. Energy 2015, 3, 141–143. [Google Scholar] [CrossRef]
  36. Mohamed, M.A.; Diab, A.A.Z.; Rezk, H.; Jin, T. A novel adaptive model predictive controller for load frequency control of power systems integrated with DFIG wind turbines. Neural Comput. Appl. 2019, 32, 7171–7181. [Google Scholar] [CrossRef]
  37. Daraz, A.; Malik, S.A.; Mokhlis, H.; Haq, I.U.; Laghari, G.F.; Mansor, N.N. Fitness Dependent Optimizer-Based Automatic Generation Control of Multi-Source Interconnected Power System with Non-Linearities. IEEE Access 2020, 8, 100989–101003. [Google Scholar] [CrossRef]
  38. Kumari, S.; Shankar, G. Novel application of integral-tilt-derivative controller for performance evaluation of load frequency control of interconnected power system. IET Gener. Transm. Distrib. 2018, 12, 3550–3560. [Google Scholar] [CrossRef]
  39. Moon, Y.H.; Ryu, H.S.; Kim, B.; Song, K.B. Optimal tracking approach to load frequency control in power systems. In Proceedings of the 2000 IEEE Power Engineering Society Winter Meeting, Conference Proceedings (Cat. No. 00CH37077). Singapore, 23–27 January 2000. [Google Scholar]
  40. Aoki, M. Control of large-scale dynamic systems by aggregation. IEEE Trans. Autom. Control. 1968, 13, 246–253. [Google Scholar] [CrossRef]
  41. Gozde, H.; Taplamacioglu, M.C.; Kocaarslan, İ. Comparative performance analysis of Artificial Bee Colony algorithm in automatic generation control for interconnected reheat thermal power system. Int. J. Electr. Power Energy Syst. 2012, 42, 167–178. [Google Scholar] [CrossRef]
  42. Hasanien, H.M.; El-Fergany, A.A. Salp swarm algorithm-based optimal load frequency control of hybrid renewable power systems with communication delay and excitation cross-coupling effect. Electr. Power Syst. Res. 2019, 176, 105938. [Google Scholar] [CrossRef]
  43. Hasanien, H.M. Whale optimisation algorithm for automatic generation control of interconnected modern power systems including renewable energy sources. IET Gener. Transm. Distrib. 2017, 12, 607–614. [Google Scholar] [CrossRef]
  44. Nguyen, T.T.; Nguyen, T.T.; Duong, M.Q.; Doan, A.T. Optimal operation of transmission power networks by using improved stochastic fractal search algorithm. Neural Comput. Appl. 2019, 32, 9129–9164. [Google Scholar] [CrossRef]
  45. Khadanga, R.K.; Kumar, A.; Panda, S. A novel sine augmented scaled sine cosine algorithm for frequency control issues of a hybrid distributed two-area power system. Neural Comput. Appl. 2021, 33, 12791–12804. [Google Scholar] [CrossRef]
  46. Sahu, B.K.; Pati, T.K.; Nayak, J.R.; Panda, S.; Kar, S.K. A novel hybrid LUS–TLBO optimized fuzzy-PID controller for load frequency control of multi-source power system. Int. J. Electr. Power Energy Syst. 2016, 74, 58–69. [Google Scholar] [CrossRef]
  47. Elkasem, A.H.A.; Khamies, M.; Magdy, G.; Taha, I.B.M.; Kamel, S. Frequency Stability of AC/DC Interconnected Power Systems with Wind Energy Using Arithmetic Optimization Algorithm-Based Fuzzy-PID Controller. Sustainability 2021, 13, 12095. [Google Scholar] [CrossRef]
  48. Khamies, M.; Magdy, G.; Hussein, M.E.; Banakhr, F.A.; Kamel, S. An Efficient Control Strategy for Enhancing Frequency Stability of Multi-Area Power System Considering High Wind Energy Penetration. IEEE Access 2020, 8, 140062–140078. [Google Scholar] [CrossRef]
  49. Khan, M.; Sun, H.; Xiang, Y.; Shi, D. Electric vehicles participation in load frequency control based on mixed H2/H∞. Int. J. Electr. Power Energy Syst. 2021, 125, 106420. [Google Scholar] [CrossRef]
  50. Sahu, R.K.; Panda, S.; Biswal, A.; Sekhar, G.T.C. Design and analysis of tilt integral derivative controller with filter for load frequency control of multi-area interconnected power systems. ISA Trans. 2016, 61, 251–264. [Google Scholar] [CrossRef]
  51. Mohanty, B.; Panda, S.; Hota, P.K. 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]
  52. Talatahari, S.; Azizi, M. Chaos Game Optimization: A novel metaheuristic algorithm. Artif. Intell. Rev. 2021, 54, 917–1004. [Google Scholar] [CrossRef]
  53. Azizi, M.; Aickelin, U.; Khorshidi, H.A.; Shishehgarkhaneh, M. Shape and Size Optimization of Truss Structures by Chaos Game Optimization Considering Frequency Constraints. J. Adv. Res. 2022, in press. [Google Scholar] [CrossRef]
  54. Coelho, L.d.S. A quantum particle swarm optimizer with chaotic mutation operator. Chaos Solitons Fractals 2008, 37, 1409–1418. [Google Scholar] [CrossRef]
  55. Zhao, W.; Wang, L.; Zhang, Z. Supply-demand-based optimization: A novel economics-inspired algorithm for global optimization. IEEE Access 2019, 7, 73182–73206. [Google Scholar] [CrossRef]
  56. Mirjalili, S.; Lewis, A.J. The whale optimization algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  57. Arora, S.; Singh, S. Butterfly optimization algorithm: A novel approach for global optimization. Soft Comput. 2019, 23, 715–734. [Google Scholar] [CrossRef]
Figure 1. The studied power grid schematic diagram.
Figure 1. The studied power grid schematic diagram.
Fractalfract 06 00220 g001
Figure 2. The transfer function model of the studied power grid.
Figure 2. The transfer function model of the studied power grid.
Fractalfract 06 00220 g002
Figure 3. The implemented model of wind power using MATLAB/Simulink program (R2015a).
Figure 3. The implemented model of wind power using MATLAB/Simulink program (R2015a).
Fractalfract 06 00220 g003
Figure 4. The output power of the wind model.
Figure 4. The output power of the wind model.
Fractalfract 06 00220 g004
Figure 5. The implemented model of the solar power plant using MATLAB/Simulink (R2015a) program.
Figure 5. The implemented model of the solar power plant using MATLAB/Simulink (R2015a) program.
Fractalfract 06 00220 g005
Figure 6. The output power of the photovoltaic model.
Figure 6. The output power of the photovoltaic model.
Fractalfract 06 00220 g006
Figure 7. The implemented model of the electrical vehicle using MATLAB/Simulink (R2015a).
Figure 7. The implemented model of the electrical vehicle using MATLAB/Simulink (R2015a).
Fractalfract 06 00220 g007
Figure 8. The construction of the proposed combining tilt-derivative and tilt-integral controller.
Figure 8. The construction of the proposed combining tilt-derivative and tilt-integral controller.
Fractalfract 06 00220 g008
Figure 9. Position updating process for the temporary triangles [53].
Figure 9. Position updating process for the temporary triangles [53].
Fractalfract 06 00220 g009
Figure 10. Qualitative metrics of nine benchmark functions using the proposed quantum chaos game optimizer algorithm: 2D views of the functions, search history, average fitness history, and convergence curve.
Figure 10. Qualitative metrics of nine benchmark functions using the proposed quantum chaos game optimizer algorithm: 2D views of the functions, search history, average fitness history, and convergence curve.
Fractalfract 06 00220 g010aFractalfract 06 00220 g010b
Figure 11. The convergence curves of the proposed QCGO algorithm and four other algorithms for 23 benchmark functions.
Figure 11. The convergence curves of the proposed QCGO algorithm and four other algorithms for 23 benchmark functions.
Fractalfract 06 00220 g011aFractalfract 06 00220 g011b
Figure 12. Boxplots of the proposed QCGO algorithm and four other algorithms for 23 benchmark functions.
Figure 12. Boxplots of the proposed QCGO algorithm and four other algorithms for 23 benchmark functions.
Fractalfract 06 00220 g012aFractalfract 06 00220 g012b
Figure 13. The convergence curve characteristics of QCGO, CGO, and salp swarm algorithm.
Figure 13. The convergence curve characteristics of QCGO, CGO, and salp swarm algorithm.
Fractalfract 06 00220 g013
Figure 14. Dynamic power grid responses in case A.1.1: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 14. Dynamic power grid responses in case A.1.1: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g014
Figure 15. Dynamic power grid responses in case A.1.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 15. Dynamic power grid responses in case A.1.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g015
Figure 16. Dynamic power grid responses in case A.1.3: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 16. Dynamic power grid responses in case A.1.3: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g016
Figure 17. The form of the applied series step load perturbation.
Figure 17. The form of the applied series step load perturbation.
Fractalfract 06 00220 g017
Figure 18. Dynamic power grid responses in case A.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 18. Dynamic power grid responses in case A.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g018aFractalfract 06 00220 g018b
Figure 19. The form of the applied random load variation.
Figure 19. The form of the applied random load variation.
Fractalfract 06 00220 g019
Figure 20. Dynamic power grid responses in case A.3: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 20. Dynamic power grid responses in case A.3: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g020aFractalfract 06 00220 g020b
Figure 21. The form of the applied series SLP.
Figure 21. The form of the applied series SLP.
Fractalfract 06 00220 g021
Figure 22. Dynamic power grid responses in case B.1: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 22. Dynamic power grid responses in case B.1: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g022aFractalfract 06 00220 g022b
Figure 23. The form of the applied RLV.
Figure 23. The form of the applied RLV.
Fractalfract 06 00220 g023
Figure 24. Dynamic power grid responses in case B.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 24. Dynamic power grid responses in case B.2: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g024aFractalfract 06 00220 g024b
Figure 25. The form of the applied RLV.
Figure 25. The form of the applied RLV.
Fractalfract 06 00220 g025
Figure 26. Dynamic power grid responses in case C: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 26. Dynamic power grid responses in case C: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g026
Figure 27. The form of the applied RLV.
Figure 27. The form of the applied RLV.
Fractalfract 06 00220 g027
Figure 28. The charging/discharging power of the applicable EVs in both areas.
Figure 28. The charging/discharging power of the applicable EVs in both areas.
Fractalfract 06 00220 g028
Figure 29. Dynamic power grid responses in case D: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Figure 29. Dynamic power grid responses in case D: (a) Δ f 1 (b) Δ f 2 (c) Δ P t i e .
Fractalfract 06 00220 g029
Table 1. The motivation of current work compared with other published works.
Table 1. The motivation of current work compared with other published works.
References[6][9][28][32][37][38]This Study
Controller
structure
PI/PID
controller
PI/PD
controller
FOPID/TID controllerCombining of FOPID-TID controllerI-PD
controller
I-TD
controller
Combining TD-TI
controller
Controller
design
adoption
Firefly
algorithm
Backtracking search
algorithm
Improved PSOManta ray
foraging
optimization
algorithm
Fitness
dependent
optimizer
Water cycle
algorithm
QCGO
Load
perturbation
challenge
SLPSLP/RLVSLP/RLVSLP/series
SLP
SLPSLP/RLVSLP/series SLP/RLV
Sort of studied systemSingle-area power
system
Multi-area power
system
Multi-area power systemMulti-area power systemMulti-area power systemMulti-area power systemMulti-area power
system
RESs
Penetration
Not
considered
Not
considered
Not
considered
consideredNot consideredNot
considered
Considered with high penetration
Effect of
communication time
delay
Not
considered
Not
considered
Not
considered
Not consideredConsidered
before the
action of one control unit only
Not
considered
Considered before and after the
control action
Effect of EVsNot
considered
Not
considered
Not
considered
Not consideredNot consideredNot
considered
Considered
Table 2. The transfer functions that are presented in the studied power grid.
Table 2. The transfer functions that are presented in the studied power grid.
Control BlockTransfer Functions
Thermal Governor 1 T s g · s + 1
Reheater of Thermal Turbine K r · T r · s + 1 T r · s + 1
Thermal Turbine 1 T t · s + 1
Hydro Governor 1 T g h · s + 1
Transient Droop Compensation T r s · s + 1 T r h · s + 1
Hydro Turbine T w · s + 1 0.5 T w · s + 1
Valve Positioner of Gas Turbine 1 b g · s + c g
Speed Governor of Gas Turbine x c · s + 1 Y c · s + 1
Fuel System and Combustor T c r · s + 1 T f c · s + 1
Gas Turbine Dynamics 1 T c d · s + 1
Power System 1 K p s 1 T p s 1 · s + 1
Power System 2 K p s 2 T p s 2 · s + 1
Electrical Vehicle 1 K E V 1 T E V 1 · s + 1
Electrical Vehicle 2 K E V 2 T E V 2 · s + 1
Table 3. The standard parameter values of the two interconnected identical areas [47].
Table 3. The standard parameter values of the two interconnected identical areas [47].
Parameter DescriptionsSymbolStandard Values
Frequency bias factor B i 0.4312 MW/Hz
Coefficient of synchronizing T 12 0.0433 MW
The regulation constant of thermal turbine
The regulation constant of hydropower plant
The regulation constant of gas turbine
R 1
R 2
R 3
2.4 HZ/MW
2.4 HZ/MW
2.4 HZ/MW
Control area capacity ratio a 12 −1
Participation factor for a thermal unit K T 0.543478
Participation factor for a hydro unit K H 0.326084
Participation factor for a gas unit K G 0.130438
Gain constant of power system K p s 68.9566
The time constant of the power system T p s 11.49 s
Governor time constant T s g 0.08 s
Turbine time constant T t 0.3 s
Gain of reheater steam turbine K r 0.3
The time constant of reheater steam turbine T r 10 s
Speed governor time constant of hydro turbine T g h 0.2 s
Speed governor reset time of the hydro turbine T r s 5 s
The transient droop time constant of hydro turbine speed governor T r h 28.75 s
Nominal string time of water in penstock T w 1 s
Gas turbine constant of valve positioner b g 0.05
Valve positioner of gas turbine c g 1
The lag time constant of the gas turbine speed governor Y c 1 s
The lead time constant of the gas turbine speed governor X c 0.6 s
Gas turbine combustion reaction time delay T c r 0.01 s
Gas turbine fuel time constant T f c 0.23 s
Gas turbine compressor discharge volume–time constant T c d 0.2 s
Gain of electrical vehicle K E V 1
The time constant of electrical vehicle T E V 0.28 s
Table 4. Results of unimodal benchmark functions.
Table 4. Results of unimodal benchmark functions.
FunctionQCGOCGOSDOWOABOA
F1Best2.4 × 10 126 1.52 × 10 58 1.39 × 10 55 1.92 × 10 40 3.87 × 10 08
Mean1.4 × 10 122 4.97 × 10 55 1.37 × 10 51 7.2 × 10 34 4.96 × 10 08
Median4.8 × 10 124 3.86 × 10 56 3.74 × 10 54 2.28 × 10 35 4.95 × 10 08
Worst1.2 × 10 121 3.9 × 10 54 8.43 × 10 51 4.34 × 10 33 6 × 10 08
Std3.7 × 10 122 9.85 × 10 55 2.74 × 10 51 1.34 × 10 33 4.94 × 10 09
F2Best4.2 × 10 65 3.64 × 10 31 1.83 × 10 29 4.41 × 10 24 4.26 × 10 06
Mean1.85 × 10 63 9.17 × 10 29 3.76 × 10 25 5.82 × 10 21 5.71 × 10 06
Median6.63 × 10 64 1.96 × 10 29 1.13 × 10 26 1.34 × 10 21 5.77 × 10 06
Worst7.99 × 10 63 9.73 × 10 28 3.98 × 10 24 5.99 × 10 20 7.58 × 10 06
Std2.41 × 10 63 2.23 × 10 28 9.1 × 10 25 1.34 × 10 20 9.92 × 10 07
F3Best2.68 × 10 42 2.41 × 10 40 6.27 × 10 46 0.0276083.85 × 10 08
Mean1.66 × 10 36 6.69 × 10 37 6.91 × 10 34 1.5183354.67 × 10 08
Median4.45 × 10 39 1.39 × 10 38 1.4 × 10 39 1.0113914.61 × 10 08
Worst1.82 × 10 35 7.13 × 10 36 1.38 × 10 32 3.9146955.57 × 10 08
Std4.41 × 10 36 1.68 × 10 36 3.09 × 10 33 1.184355.02 × 10 09
F4Best5.12 × 10 53 3.76 × 10 37 1.11 × 10 26 0.995288.45 × 10 06
Mean6.71 × 10 51 3.7 × 10 23 4.52 × 10 23 53.183951.02 × 10 05
Median2.13 × 10 51 1.4 × 10 23 1.14 × 10 23 60.931681.02 × 10 05
Worst3.32 × 10 50 1.81 × 10 22 1.94 × 10 22 89.099691.15 × 10 05
Std9.43 × 10 51 5.38 × 10 23 6.34 × 10 23 29.695438.51 × 10 07
F5Best18.1158217.1184527.9096727.8848328.89058
Mean19.5786119.6102628.6509628.2741928.92369
Median19.3562219.2926528.7472628.4364728.91978
Worst22.217521.5946328.9869928.722728.96927
Std1.1496091.2248820.2950260.289250.021273
F6Best1.75 × 10 14 6.75 × 10 14 0.0399570.3035424.311051
Mean2.86 × 10 12 2.63 × 10 12 2.5685410.6559075.211726
Median7.7 × 10 14 6.23 × 10 13 2.0387790.622035.06303
Worst4.89 × 10 11 2.57 × 10 11 7.2502511.164086.168001
Std1.09 × 10 11 6.11 × 10 12 1.8527010.2108110.509499
F7Best1.02 × 10 05 0.0001978.66 × 10 05 0.00040.000983
Mean0.0002630.000920.0023560.005420.002696
Median0.0002310.000850.0011360.0037630.002776
Worst0.0007680.0019750.0138130.0190690.005116
Std0.0001770.0005830.0033310.0050110.001104
The best values obtained are in bold.
Table 5. Results of multimodal benchmark functions.
Table 5. Results of multimodal benchmark functions.
FunctionQCGOCGOSDOWOABOA
F8Best−1671.01−1770.26−1655−1909.05−921.028
Mean−1465.24−1490.19−1312.83−1786.9−766.513
Median−1453.48−1483.32−1385.86−1907.06−778.594
Worst−1313.6−1235.22−598.802−1632.06−647.792
Std108.2831123.7418294.008138.075961.76107
F9Best0.000.004.33 × 10 30 0.005.17 × 10 09
Mean0.000.001.75 × 10 22 1.14 × 10 14 0.003376
Median0.000.004.17 × 10 25 0.003.86 × 10 06
Worst0.000.003.02 × 10 21 1.14 × 10 13 0.047754
Std0.000.006.75 × 10 22 2.97 × 10 14 0.010836
F10Best8.88 × 10 16 8.88 × 10 16 8.88 × 10 16 4.44 × 10 15 1.67 × 10 05
Mean2.49 × 10 15 3.2 × 10 15 8.88 × 10 16 1.33 × 10 14 4.77 × 10 05
Median8.88 × 10 16 4.44 × 10 15 8.88 × 10 16 1.15 × 10 14 4.55 × 10 05
Worst4.44 × 10 15 4.44 × 10 15 8.88 × 10 16 3.29 × 10 14 7.94 × 10 05
Std1.81 × 10 15 1.74 × 10 15 0.008.11 × 10 15 1.69 × 10 05
F11Best0.000.000.000.003.23 × 10 08
Mean0.000.000.000.0218324.29 × 10 08
Median0.000.000.000.004.22 × 10 08
Worst0.000.000.000.266265.81 × 10 08
Std0.000.000.000.0689736.29 × 10 09
F12Best3.66 × 10 16 1.34 × 10 15 0.0011520.0060520.33315
Mean5.69 × 10 15 8.04 × 10 14 0.234670.0222390.565424
Median2.26 × 10 15 1.93 × 10 14 0.0678050.0155290.562862
Worst3.32 × 10 14 5.01 × 10 13 1.4928210.0879470.754521
Std8.06 × 10 15 1.36 × 10 13 0.3520630.0187740.108748
F13Best6.36 × 10 14 7.4 × 10 13 0.0462160.4002812.497296
Mean0.0071420.0367331.8675520.6875222.894224
Median0.0054940.0109871.9342460.5980542.982946
Worst0.0439490.2334142.9999241.3213523.109356
Std0.0102540.0659780.9612840.2485230.153028
The best values obtained are in bold.
Table 6. Results of composite benchmark functions.
Table 6. Results of composite benchmark functions.
FunctionQCGOCGOSDOWOABOA
F14Best0.9980040.9980040.9980040.9980040.998004
Mean0.9980040.9980043.4946962.2302041.301281
Median0.9980040.9980041.4950171.4950171.024436
Worst0.9980040.99800412.6705110.763182.983027
Std0.005.09 × 10 17 3.9532032.2413670.534994
F15Best0.0003070.0003070.0003070.0003110.000315
Mean0.0003070.0003530.000670.0006260.000487
Median0.0003070.0003070.0005270.0005780.000405
Worst0.0003070.0012230.0021210.0015280.000917
Std1.68 × 10 19 0.0002050.0004730.0003420.000173
F16Best−1.03163−1.03163−1.03163−1.03163−1.40747
Mean−1.03163−1.03163−1.03005−1.03163−1.18199
Median−1.03163−1.03163−1.03163−1.03163−1.18517
Worst−1.03163−1.03163−1.00046−1.03163−1.07213
Std2.22 × 10 16 2.28 × 10 16 0.0069661.94 × 10 08 0.088213
F17Best0.3978870.3978870.3978870.3978870.398293
Mean0.3978870.3978870.3979870.3978960.409332
Median0.3978870.3978870.3978870.397890.406611
Worst0.3978870.3978870.3997950.3979670.461881
Std0.000.000.0004261.78 × 10 05 0.014049
F18Best3.003.003.003.0000013.000586
Mean3.003.003.0011853.0000693.092676
Median3.003.003.003.0000263.054728
Worst3.003.003.0235373.0006683.425476
Std2.7 × 10 16 6.03 × 10 16 0.0052610.0001470.108993
F19Best−0.30048−0.30048−0.30048−0.30048−0.30048
Mean−0.30048−0.30048−0.2893−0.30048−0.30048
Median−0.30048−0.30048−0.30038−0.30048−0.30048
Worst−0.30048−0.30048−0.19165−0.30048−0.30048
Std1.14 × 10 16 1.14 × 10 16 0.0265311.14 × 10 16 3.74 × 10 06
F20Best−3.322−3.322−3.322−3.31923 3.3 × 10 05
Mean−3.26849−3.28038−3.09697−2.98949 1.6 × 10 06
Median−3.322−3.322−3.2031−3.15019 1.5 × 10 40
Worst−3.2031−3.2031−0.89904−1.57922 2 × 10 134
Std0.0606850.0581820.5509860.4797957.35 × 10 06
F21Best−10.1532−10.1532−10.1532−10.1528−4.61081
Mean−10.1532−9.90058−8.703−7.35262−4.0759
Median−10.1532−10.1532−10.1532−10.0113−4.12522
Worst−10.1532−5.10077−4.99677−2.59723−3.18003
Std3.21 × 10 15 1.1297572.239523.2454450.379957
F22Best−10.4029−10.4029−10.4029−10.4008−4.76031
Mean−10.4029−10.4029−8.45822−7.90953−3.74931
Median−10.4029−10.4029−10.4029−10.2376−3.64889
Worst−10.4029−10.4029−1.0677−3.69711−2.93305
Std3.05 × 10 15 3.36 × 10 15 3.1286892.7797440.479377
F23Best−10.5364−10.5364−10.5364−10.5297−4.51577
Mean−9.99562−9.93332−7.90449−7.3919−3.38426
Median−10.5364−10.5364−10.5357−7.79854−3.60414
Worst−5.12848−3.83543−3.79083−1.67334−1.95854
Std1.6645251.8689523.0153193.339090.720921
The best values obtained are in bold.
Table 7. The optimum parameters of the different controllers.
Table 7. The optimum parameters of the different controllers.
Controller PropertiesThermalHydroGas
Combining TD-TI-based QCGO k t 1 = 9.9999 ,   k d 1 = 9.9988 ,   n 1 = 3.5626
k t 2 = 9.9991 ,   k i 2 = 5.4425 ,   n 2 = 3.5311
k t 1 = 9.9834 ,   k d 1 = 3.8871 ,   n 1 = 9.9468
k t 2 = 9.5835 ,   k i 2 = 1.0016 ,   n 2 = 9.9508
k t 1 = 9.998 ,   k d 1 = 9.9973 ,   n 1 = 3.7621
k t 2 = 9.9951 ,   k i 2 = 9.9704 ,   n 2 = 1.2938
Combining TD-TI-based CGO k t 1 = 9.9998 ,   k d 1 = 6.9628 ,   n 1 = 3.5715
k t 2 = 9.9977 ,   k i 2 = 5.033   n 2 = 3.4737
k t 1 = 9.98 ,   k d 1 = 2.7245 ,   n 1 = 9.9129
k t 2 = 7.2945 ,   k i 2 = 1.052   n 2 = 9.9827
k t 1 = 9.9998 ,   k d 1 = 8.4098 ,   n 1 = 1.2782
k t 2 = 9.9966 ,   k i 2 = 9.9989   n 2 = 6.9549
Combining TD-TI-based SSA k t 1 = 9.9998 ,   k d 1 = 8.985 ,   n 1 = 2.9819
k t 2 = 9.1794 ,   k i 2 = 9.3854   n 2 = 2.8288
k t 1 = 5.3557 ,   k d 1 = 4.68 ,   n 1 = 2.1217
k t 2 = 8.5211 ,   k i 2 = 1.0925 ,   n 2 = 5.1176
k t 1 = 9.9998 ,   k d 1 = 1.0849 ,   n 1 = 9.6003
k t 2 = 9.9628 ,   k i 2 = 7.6555 ,   n 2 = 1.4599
TID-based QCGO k t 1 = 9.8753 ,   k i 1 = 9.9302 ,   k d 1 = 7.9837,
n 1 = 2.6219
k t 1 = 9.7665 ,   k i 1 = 1.0797 ,   k d 1 = 4.9139,
n 1 = 8.0894
k t 1 = 9.9041 ,   k i 1 = 9.9922 ,   k d 1 = 1.6516,
n 1 = 9.2214
TID-based CGO k t 1 = 9.9993 ,   k i 1 = 9.7827 ,   k d 1 = 8.7199,
n 1 = 3.5979
k t 1 = 9.9525 ,   k i 1 = 1.4282 ,   k d 1 = 5.1353,
n 1 = 7.5851
k t 1 = 9.9486 ,   k i 1 = 9.9844 ,   k d 1 = 4.0435,
n 1 = 3.3106
PID-based TLBO [46] k p 1 = 4.1468 ,   k i 1 = 4.0771 ,   k d 1 = 2.0157 k p 1 = 1.0431 ,   k i 1 = 0.6030 ,   k d 1 = 2.2866 k p 1 = 4.7678 ,   k i 1 = 3.7644 ,   k d 1 = 4.9498
PID-based AOA [47] k p 1 = 10 ,   k i 1 = 1.5975 ,   k d 1 = 2.7449 k p 1 = 1.5975 ,   k i 1 = 0.0837 ,   k d 1 = 0.0875 k p 1 = 10 ,   k i 1 = 10 ,   k d 1 = 1.2779
Table 8. The transient response specifications of the presented system for case A.1.1.
Table 8. The transient response specifications of the presented system for case A.1.1.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e ) Objective Function Value (ITAE)
Combining TD-TI based on QCGO
O s h and U s h
  × ( 10 3 )
O s h = 0.819
U s h = −7.875
O s h = 0.0028
U s h = −1.744
O s h = 0.0015
U s h = −0.5361

J = 0.075
PID based on AOA
O s h and U s h
× ( 10 3 ) [47]
O s h = 1.158
U s h = −11.42
O s h = 0.02096
U s h = −4.443
O s h = 0.01107
U s h = −1.249

J = 0.189
PID based on TLBO
O s h and U s h
× ( 10 3 ) [46]
O s h = 1.7217
U s h = −19.7259
O s h = 0.4363
U s h = −12.7986
O s h = 0.1712
U s h = −3.0782
J = 0.402
Table 9. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.1.
Table 9. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.1.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO60.01  52.43 86.4  99.3682.6  99.12
PID based on AOA42.11  32.7065.29  95.259.42  93.53
The optimum values are bolded.
Table 10. The transient response specifications of the presented system for case A.1.2.
Table 10. The transient response specifications of the presented system for case A.1.2.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e ) Objective Function Value (ITAE)
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 0.819
U s h = −7.875
O s h = 0.0028
U s h = −1.744
O s h = 0.0015
U s h = −0.5361

J = 0.075
TID based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 1.893
U s h = −11.468
O s h = 0.3257
U s h = −3.45
O s h = 0.0424
U s h = −0.8862
J = 0.1351
TID based on CGO
O s h and U s h
× ( 10 3 )
O s h = 1.705
U s h = −10.341
O s h = 0.3784
U s h = −2.763
O s h = 0.0381
U s h = −0.7397

J = 0.1381
The optimum values are bolded.
Table 11. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.2.
Table 11. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.2.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO60.01   52.4386.4   99.3682.6   99.12
TID based on QCGO41.86   −9.9573.04   25.3571.21   75.23
TID based on CGO47.6   0.9778.4   13.2775.97   77.75
The optimum values are bolded.
Table 12. The transient response specifications of the presented system for case A.1.3.
Table 12. The transient response specifications of the presented system for case A.1.3.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e ) Objective Function Value (ITAE)
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 0.819
U s h = −7.875
O s h = 0.0028
U s h = −1.744
O s h = 0.0015
U s h = −0.5361

J = 0.075
Combining TD-TI based on CGO
O s h and U s h
× ( 10 3 )
O s h = 1.097
U s h = −8.95
O s h = 0.0025
U s h = −2.383
O s h = 0.00136
U s h = −0.665

J = 0.078
Combining TD-TI based on SSA
O s h and U s h
× ( 10 3 )
O s h = 1.763
U s h = −9.978
O s h = 0.0896
U s h = −2.713
O s h = 0.0124
U s h = −0.7125

J = 0.087
Table 13. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.3.
Table 13. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and PID/AOA based on PID controller via TLBO for scenario A.1.3.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO60.01   52.4386.4   99.3682.6   99.12
Combining TD-TI based on CGO54.63   36.2881.38   99.4378.4   99.21
Combining TD-TI based on SSA49.42   −2.478.8   79.4676.85   92.76
The optimum values are bolded.
Table 14. The transient response specifications of the presented system for case A.2.
Table 14. The transient response specifications of the presented system for case A.2.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e )
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 15.6
U s h = −22.9
O s h = 3.5
U s h = −5.1
O s h = 1.000
U s h = −1.67
Combining TD-TI based on CGO
O s h and U s h
× ( 10 3 )
O s h = 18.00
U s h = −26.1
O s h = 4.85
U s h = −7.3
O s h = 1.3
U s h = −1.9
Combining TD-TI based on SSA
O s h and U s h
× ( 10 3 )
O s h = 21.000
U s h = −31.000
O s h = 5.510
U s h = −8.6
O s h = 1.40
U s h = −2.15
Table 15. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario A.2.
Table 15. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario A.2.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO26.13    25.7140.7   36.4822.33   28.6
Combining TD-TI based on CGO15.81   14.2915.12   11.9811.63   6.43
The optimum values are bolded.
Table 16. The transient response specifications of the presented system for case A.3.
Table 16. The transient response specifications of the presented system for case A.3.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e )
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 7.4
U s h = −11.9
O s h = 1.4
U s h = −2.2
O s h = 0.51
U s h = −0.76
Combining TD-TI based on CGO
O s h and U s h
× ( 10 3 )
O s h = 8.40
U s h = −13.5
O s h = 2.3
U s h = −3.6
O s h = 0.65
U s h = −1.000
Combining TD-TI based on SSA
O s h and U s h
× ( 10 3 )
O s h = 10.000
U s h = −15.000
O s h = 2.60
U s h = −4.08
O s h = 0.72
U s h = −1.14
Table 17. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario A.3.
Table 17. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario A.3.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO20.67   26.0046.08   46.1533.33   29.17
Combining TD-TI based on CGO10.00   16.0011.76   11.5412.28   9.72
The optimum values are bolded.
Table 18. The transient response specifications of the presented system for case B.1.
Table 18. The transient response specifications of the presented system for case B.1.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e )
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 71.0
U s h   = −22.0
O s h = 40.3
U s h = −7.4
O s h = 4.8
U s h = −2.7
Combining TD-TI based on CGO
O s h and U s h
× ( 10 3 )
O s h = 81.0
U s h = −27.0
O s h = 46.0
U s h = −8.1
O s h = 6.1
U s h = −3.6
Combining TD-TI based on SSA
O s h and U s h
× ( 10 3 )
O s h = 96.000
U s h = −30.000
O s h = 51.1
U s h = −9.5
O s h = 6.5
U s h = −3.88
Table 19. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario B.1.
Table 19. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario B.1.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO26.67   26.0422.11   21.1430.41   26.15
Combining TD-TI based on CGO10.00   15.6314.74   9.987.22   6.15
The optimum values are bolded.
Table 20. The transient response specifications of the presented system for case B.2.
Table 20. The transient response specifications of the presented system for case B.2.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e )
Combining TD-TI based on QCGO
O s h and U s h
× ( 10 3 )
O s h = 72.0
U s h = −10.0
O s h = 40.0
U s h = −4.0
O s h = 4.7
U s h = −2.5
Combining TD-TI based on CGO
O s h and U s h
× ( 10 3 )
O s h = 81.0
U s h = −12.0
O s h = 46.0
U s h = −5.1
O s h = 6.0
U s h = −3.69
Combining TD-TI based on SSA
O s h and U s h
× ( 10 3 )
O s h = 93.000
U s h = −16.000
O s h = 51.4
U s h = −9.4
O s h = 6.4
U s h = −3.83
Table 21. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario B.2.
Table 21. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario B.2.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO37.50   22.5857.45   22.1834.73   26.56
Combining TD-TI based on CGO25.00   12.945.74   10.513.66   6.25
The optimum values are bolded.
Table 22. The transient response specifications of the presented system for case D.
Table 22. The transient response specifications of the presented system for case D.
Controller PropertiesDynamic Response of ( Δ f 1 ) Dynamic Response of ( Δ f 2 ) Dynamic Response of ( Δ P t i e )
Combining TD-TI based on QCGO with EVs
O s h and U s h
× ( 10 3 )
O s h = 62.1
U s h = −8.4
O s h = 36.0
U s h = −1.9
O s h = 4.16
U s h = −2.2
Combining TD-TI based on QCGO without EVs
O s h and U s h
× ( 10 3 )
O s h = 72.0
U s h = −10.0
O s h = 40.0
U s h = −4.0
O s h = 4.7
U s h = −2.5
Table 23. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario D.
Table 23. Percentage improvement in U s h and O s h values for combining TD-TI/QCGO and combining TD-TI/CGO based on combining TD-TI/SSA for scenario D.
Controller
Δ f 1 U s h O s h
Δ f 2 U s h O s h
Δ P t i e U s h O s h
Combining TD-TI based on QCGO with EVs47.5   33.2379.79   29.9642.56   35
Combining TD-TI based on QCGO without EVs37.5   22.5857.45   22.1834.73   26.56
The optimum values are bolded.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/fractalfract6040220

AMA Style

Elkasem AHA, Khamies M, Hassan MH, Agwa AM, Kamel S. Optimal Design of TD-TI Controller for LFC Considering Renewables Penetration by an Improved Chaos Game Optimizer. Fractal and Fractional. 2022; 6(4):220. https://doi.org/10.3390/fractalfract6040220

Chicago/Turabian Style

Elkasem, Ahmed H. A., Mohamed Khamies, Mohamed H. Hassan, Ahmed M. Agwa, and Salah Kamel. 2022. "Optimal Design of TD-TI Controller for LFC Considering Renewables Penetration by an Improved Chaos Game Optimizer" Fractal and Fractional 6, no. 4: 220. https://doi.org/10.3390/fractalfract6040220

APA Style

Elkasem, A. H. A., Khamies, M., Hassan, M. H., Agwa, A. M., & Kamel, S. (2022). Optimal Design of TD-TI Controller for LFC Considering Renewables Penetration by an Improved Chaos Game Optimizer. Fractal and Fractional, 6(4), 220. https://doi.org/10.3390/fractalfract6040220

Article Metrics

Back to TopTop