Next Article in Journal
A Generic Framework for Multiscale Simulation of High and Low Enthalpy Fractured Geothermal Reservoirs under Varying Thermodynamic Conditions
Previous Article in Journal
Part-I: State-of-the-Art Technologies of Solar Powered DC Microgrid with Hybrid Energy Storage Systems-Architecture Topologies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Participation of Renewable Energy Sources in the Frequency Regulation Issues of a Five-Area Hybrid Power System Utilizing a Sine Cosine-Adopted African Vulture Optimization Algorithm

by
Smruti Ranjan Nayak
1,
Rajendra Kumar Khadanga
1,*,
Sidhartha Panda
2,
Preeti Ranjan Sahu
3,
Sasmita Padhy
3 and
Taha Selim Ustun
4,*
1
Department of Electrical and Electronics Engineering, Centurion University of Technology and Management, Bhubaneswar 752050, Odisha, India
2
Department of Electrical Engineering, Veer Surendra Sai University of Technology, Burla 768018, Odisha, India
3
Department of Electrical and Electronics Engineering, NIST Institute of Science and Technology, Berhampur 7601008, Odisha, India
4
Fukushima Renewable Energy Institute, AIST (FREA), Fukushima 963-0298, Japan
*
Authors to whom correspondence should be addressed.
Energies 2023, 16(2), 926; https://doi.org/10.3390/en16020926
Submission received: 5 December 2022 / Revised: 8 January 2023 / Accepted: 11 January 2023 / Published: 13 January 2023

Abstract

:
In this article, a novel methodology is proposed by utilizing a technique which, in light of the change in the African vulture optimization known as Sine Cosine, adopted an African vulture optimization algorithm (SCaAVOA)-based tilt integral derivative (TID) regulator for the load frequency control (LFC) of a five-area power system with multi-type generations. At first, the execution of the Sine Cosine-adopted calculation is tried by contrasting it with the standard AVOA calculation while considering different standard benchmark functions. To demonstrate the superiority of the proposed SCaAVOA algorithm, the results are contrasted using different standard approaches. In the next stage, the proposed method is used in a five-area thermal power system and is likewise applied to a five-area, ten-unit system comprising different conventional sources as well as some renewable energy sources. The performance analysis of the planned regulator is completed for various system boundaries and loading conditions. It is seen that the said regulator is more viable in comparison to the other standard controllers.

1. Introduction

The standard undertaking of LFC is to affirm the frequency variation inside a theorized reach that is past what many would consider conceivable [1]. A fundamental issue related to LFC is how to maintain the system frequency within a certain limit [2]. The frequency regulation of a two-region interconnected power system is the focal point of research on LFC. Not many of these studies have assessed conventional energy and, incredibly, fewer have considered the effect of distributed energy sources on the LFC plan [3,4,5]. A three-area system is utilized to check the expansion of a two-area system [6], and a five-area power system is used to verify the extension of a three-area power system [7,8]. The application of different renewable energy sources in the AGC application is also explained in [9,10,11]. In any case, no exploration has considered the effect of coordinating different environmentally friendly renewable energy sources in a five-region power system. Therefore, in this work, a five-region power framework with various distributed energy sources, such as wind generators, solar generators, plug-in electric vehicles, and other sources including fuel cells, microturbines, and diesel engine generators, are taken into consideration in different areas of the system.
It is obvious from the literature that the Proportional-Integral (PI) regulators are generally utilized in AGC because of their basic executive strategy, while simultaneously giving satisfactory system responses. In any case, the presentation of a basic PI regulator decays as far as the overshoots and the settling time of the system response [12]. The PID regulator is an exceptionally famous construction because of its simple execution. However, it is not efficient to control a non-linear system that leads to unbounded responses [13]. Therefore, researchers have focused on the improvement of the PID controller such that the stability region and the application can be extended to highly nonlinear, dynamic, and complex systems [14,15,16]. In the TID regulator, the comparable part to the PID regulator is displaced with a shifted part having a transfer function s-1/n, thus achieving an upgraded feedback controller [17]. This paper uses a TID regulator associated with the frequency control issue of the five-region power system. In the proposed TID regulator, the comparable part to the PID regulator is displaced with a shifted part having a transfer function s-1/n, thus achieving an upgraded feedback controller.
Involving an evolutionary algorithm as an answer for the LFC issue is one choice (EA). The essential objective of involving an EA is to have the option of overseeing nonlinear capabilities [18]. EA applications incorporate the genetic algorithm [19], artificial bee colony algorithm [20], particle swarm optimization [21], butterfly optimization algorithm [22], biogeography-based optimization [23], satin bowerbird algorithm [24], grasshopper optimization algorithm [25], gravitational search algorithm [26], coyote optimization algorithm [27], whale optimization algorithm (WOA) [28], novel COVID-19-based optimization [29], equilibrium optimization [30], etc.
Although these strategies give a strong execution and deal with a reasonable LFC structure, their convergence rate is very sluggish and they are oftentimes caught in local optima. In recent years, another calculation known as the African vulture optimization algorithm (AVOA) has grown in prevalence and has been comprehensively used in various optimization issues. The AVOA is inspired by the foraging and navigational behaviours of African vultures and aids in improvement measures through an optimization procedure [31]. Although the said calculation has been, by and large, used, it has a couple of limitations. The most concerning issue related to the AVOA is that it effectively falls into local optimum [32]. In light of the previously mentioned issue, this paper proposes a clever methodology by bringing some scaling factors into the first AVOA calculation, known as the sine cosine-adopted AVOA algorithm, for the LFC of five-area hybrid multi-type generations.

2. Research Gap, Inspiration, and Paper Organization

2.1. Research Gap, Inspiration, and Paper Organization

The principal perceptions acquired from a review of the above studies are as follows:
  • Almost no exploration has considered the effect of coordinating different renewable power sources into an AGC.
  • Again, the significance of different important disturbances, such as solar irradiation, wind power change, and load change, on an AGC system has not been examined in previous studies.
  • To the best of the authors’ knowledge, an execution of a tilt-integral-derivative control in an AGC five-area renewable energy integrated hybrid power system for LFC applications is missing.

2.2. Inspiration and Paper Organization

Under typical operating conditions, the governor, which is in primary control, may attenuate small frequency deviations, while the LFC, which is in charge of the secondary control, is in charge of big frequency deviations. Control techniques, such as optimization, and intelligent controllers, such as fuzzy, have been applied in power system control to improve the function and performance of a system during normal and abnormal conditions. The HVDC components and their placement in the power system can also improve the operation of the AGC. The key contributions of this study are as follows:
  • The SCaAVOA algorithm’s advantages over a few alternative tactics are analyzed in terms of execution time and objective function.
  • A five-area power system’s response to the coordination of various distributed energy sources is examined.
  • For the aforementioned hybrid five-area system, a TID controller is set up using the newly introduced SCaAVOA computation, and by differentiating it from some other existing regulators/controllers, it is shown to provide a better frequency regulation.

3. Proposed System under Study

Figure 1 shows a one-area interconnected thermal power system. The single-line diagram of an interconnected five-area hybrid power system comprising some conventional, non-conventional, and other distributed generation (DG) sources is shown in Figure 2. Table 1 shows the gains (K) and time constant (T) of the system parameters for the abovementioned five-area system.

3.1. Power System Component Modelling

3.1.1. Photovoltaic (PV) System

The basic components of a PV system are the PV module, the MPPT tracker, the converter, and the filter circuit, which can be represented as follows [12]:
G P V ( s ) = K P V 1 + s T P V

3.1.2. WTG System

The wind turbine system is a non-linear system [33], and it can be addressed in low-frequency oscillation as follows:
G W T G ( s ) = K W T G 1 + s T W T G

3.1.3. Plug-in Electric Vehicle

A PEV can be characterized as a vehicle that draws power from a battery and is good for being charged from an external source. PEVs are, at present, expected to account for the reasonable power control intended to deal with the issues resulting from the significant emergence of renewable energy. In this manner, PEV can be addressed as follows [34]:
G P E V ( s ) = K P E V 1 + s T P E V = Δ P P E V Δ U

3.1.4. Hydrogen Aqua Electrolyzer

A HAE is one of the renewal energy sources which can be represented as follows:
G H A E ( s ) = K H A E 1 + s T H A E

3.1.5. Fuel Cell

A fuel cell utilizes the compound energy of hydrogen or different fuels to produce electric power neatly and productively. Its power device is comparably a critical part, considering its diminished contamination level and extended efficiency, which can be expressed as follows [35]:
G F C ( s ) = K F C 1 + s T F C

3.1.6. DEG System

Diesel motor generators are fit for addressing the shortfall of power and can limit power awkwardness, among other generators, in the power market. The transfer function for a DEG system can be expressed as follows [34]:
G D E G ( s ) = K D E G 1 + s T D E G

3.1.7. BESS System

As a prerequisite, a BESS goes about as a source or a load according to the necessity of the system and can be expressed as follows [36]:
G B E S S ( s ) = K B E S S 1 + s T B E S S

3.1.8. Thermal Power System

A non-reheat thermal system is considered for the proposed study where the value of the GRC is taken as 12% p.u. MW/minute, which can be communicated as follows:
T u r b i n e G t ( s ) = K t 1 + s T t
G e n e r a t o r G P ( s ) = K P 1 + s T P
G o v e r n o r G g ( s ) = K g 1 + s T g
R e h e a t e r G r ( s ) = 1 + s K r T r 1 + s T r

3.1.9. System Modelling

To represent the system with load, the following transfer function can be used:
G p ( s ) = K p 1 + s T p

4. Tilt-Integral-Derivative (TID) Controller

A TID regulator is fundamentally a tunable compensator having KP, KI, and KD as the three control boundaries with the tuning boundary as n. The construction of a TID controller is like that of a PID regulator: the proportional behaviour is replaced by a shifted comparable way of behaving having a transfer function 1/S1/n or s−1/n. The tilted conduct gives an input gain as a part of the frequency, which is shifted, concerning the gain of the common compensator which is insinuated as the TID compensator, as shown in Figure 3. Mathematically, it can be stated as follows [16]:
U TID = G TID s , θ E TID ( s )
E TID = R TID s , θ Y TID ( s )
where G TID s , θ shows the transfer for the TID controller, S C and θ R 4 . The G TID s , θ can be expressed as follows [17]:
G TID s , θ = KT s 1 n + KI s + KDs
where θ T = [ K T K I K D n ] and θ R 4 , n R , and n 0

5. Optimization Problem

In the proposed work, the goal is the minimization of frequency deviation by considering an ITAE objective function, which can be expressed as follows [35]:
J = ITAE = 0 t sin Δ F i + Δ P Tie i k t dt
where Δ F i , Δ P Tie i k , and tsim show the ith area frequency deviations, tie-line power deviation between the ith and kth areas, and simulation time, respectively.
The fitness function and constraints for this formulated optimization problem are given as follows:
Minimize J
Subject to K P Y K x K P Z
K I Y K x K I Z K D Y K x K D Z
where KX is the constant for controller gains, KXY denotes the smallest estimation, and KXZ is the highest value of the controller parameters. To improve control performance, four controller parameters, or twelve (20) parameters for five sections of TID, should be obtained or correctly tuned. The controller parameters’ range is taken to be between −2 and +2, and n is kept between 1 and 10 [25].

6. Proposed Sine Cosine-Adopted African Vultures Optimization Algorithm

6.1. African Vultures Optimization Algorithm

The AVOA is inspired by the foraging and navigational behaviours of African vultures. The biological nature of vultures regarding searching and competing for food is described in four different steps here [31].

6.1.1. Findings of the Best Solution

Let there be N vultures in the environment which indicates the number of population s.t n= {1, 2,…, N}. Then, the fitness function of each position is calculated. Let pn be the probability of selecting the first or second group, which is calculated as follows:
p n = F n n = 1 N F n
where Fn is the fitness function of the nth position. Then, the formation of the first and second groups in each iteration is obtained by the following equation:
R ( it ) = first   group , p n = L 1 second   group , p n = L 2
where the ranges of L1 and L2 are 0 L 1 , L 2 1 and L1 + L2 = 1.

6.1.2. The Rate of Starvation

The satiated vultures with enough energy can move a long distance to search for food while the hungry one cannot fly longer. The rate of being satiated or hungry describes the movement from the exploration phase to the exploitation phase, which is given as follows [32]:
A = 2 × r a n d 1 + 1 × x × 1 i t I T max + y
y = h × sin z Π 2 × i t I T max + cos Π 2 × i t I T max 1
where A denotes the vultures with high energy; it and ITmax denote the current and maximum iteration, respectively; x, h, and rand1 are the random numbers varying between −1 to 1, −2 to 2, and 0 to 1, respectively; and z defines the probability of entering the exploration stage.

6.1.3. Exploration

The process of finding food by African vultures describes the exploration phase, in which the parameter p1, 0 p 1 1 decides the selection of strategy. Hence,
P ( i t + 1 ) = 6 , p 1 r a n d 2 ( 8 ) , p 1 < r a n d 2
P ( it + 1 ) = R ( it ) D ( it ) × A
D ( it ) = q × R ( it ) P ( it )
where P ( it + 1 ) represent the vultures’ position vector in the next iteration. A and R(it) are obtained from (3) and (2), respectively. q = 2 × r a n d 3 , where rand3 is a random number between 0 and 1.
P ( it + 1 ) = R ( it ) A + r a n d 4 × u b l b × r a n d 5 + l b
ub and lb are the upper and lower bounds of the variable, respectively. Rand4 and Rand5 are the random numbers between 0 and 1.

6.1.4. Exploitation

The exploitation stage has two phases with two different strategies. The selection of any strategy depends on the parameters p2 and p3. The values of p2 and p3 are between 0 and 1. When F lies between 0.5 and 1, the exploitation stage enters the first stage which describes two different strategies, such as rotating flight and siege fight. Hence, the next position is updated as follows:
P ( it + 1 ) = 10 , p 2 r a n d 6 ( 11 ) , p 2 < r a n d 6
where rand6 is a random number between 0 and 1. The solutions to (27) are given as follows:
P ( it + 1 ) = D ( it ) × A + r a n d 7 d ( it )
D ( it ) = R ( it ) P ( it )
Then, the rotational flight of vultures is modelled as follows:
P ( i t + 1 ) = R ( i t ) M 1 + M 2
M 1 = R ( i t ) × r a n d 8 × P ( i t ) 2 π × cos P ( i t )
M 2 = R ( i t ) × r a n d 9 × P ( i t ) 2 π × sin P ( i t )
where rand8 and rand9 are two random numbers between 0 and 1.
The exploitation phase enters the second phase when F > 0.5 , which explains the two different strategies of vultures, such as accumulation and aggressive siege and fight to find food. The selection of any strategy is based on the following condition:
P ( i t + 1 ) = 16 , p 3 r a n d 7 ( 19 ) , p 3 < r a n d 7
where
P ( i t + 1 ) = B 1 + B 2 2
B 1 = B e s t v u l t u r e 1 ( i t ) B e s t v u l t u r e 1 ( i t ) × P ( i t ) B e s t v u l t u r e 1 ( i t ) × P ( i t ) 2 × A
B 2 = B e s t v u l t u r e 2 ( i t ) B e s t v u l t u r e 2 ( i t ) × P ( i t ) B e s t v u l t u r e 2 ( i t ) × P ( i t ) 2 × A
B e s t v u l t u r e 1 ( i t ) and B e s t v u l t u r e 2 ( i t ) represent the best vultures of the first and second groups, respectively, in the current iteration. The aggressive competition among the vultures is represented as follows:
P ( i t + 1 ) = R ( i t ) d ( i t ) × A × Levy ( d 1 )
where d1 represents the dimension of the optimization problem. Levy(d1) is calculated as follows:
Levy ( x ) = 0.01 × u × σ v 1 / β , σ = ϒ ( 1 + β ) × sin Π β 2 ϒ ( 1 + 2 β ) × β × 2 β 1 2 1 / β

6.2. Sine Cosine-Adopted African Vultures Optimization Algorithm (SCaAVOA)

In the beginning phases of the AVOA calculation, the best solution is not known. Consequently, utilizing large steps at first might bring about the moving of the calculation far away from the optimal position. Subsequently, scaling variables can be utilized to change the situations during the beginning phases of the calculation. In the proposed SCaAVOA strategies, at the end of each iteration, the vultures are repositioned using the sine and cosine adopted scaling factors, as given below:
P ( i t + 1 ) N e w = P ( i t + 1 ) × SCaSF
where SCaSF shows the scaling factors as
SCaSF = sin ( W 1 W 2 i t M a x _ i t ) i f R N D 1 < 0.5 cos ( W 1 W 2 i t M a x _ i t ) i f R N D 1 0.5
where RND1 shows a random number; W1 and W2 are the weighting factors; and it and Max_it are the present and maximum iteration. The scaling factors are utilized to control the development of the vultures during the beginning phases of the calculation. For the appropriate choice of W1 and W2, various upsides of W are appointed and tried. It is seen that the best outcomes are acquired when W1 and W2 are chosen as 10 and 9, respectively. The consideration of the scaling factor stage changes the vultures’ positions during the underlying periods of the pursuit interaction in this manner, further developing the hunting capacity of the algorithm. For the legitimate exploitation of search space, the calculation ought to have the option to track down a better position between two positions. For the exploration of the search space, the arrangement can look past the space. This cycle can ensure a better exploitation and exploration capacity of the calculations.

7. Discussion on the Simulation Results

7.1. Performance Evaluation of the Proposed Sine Cosine-Adopted AVOA

As a first step, the performance of the proposed SCaAVOA technique has been investigated with some benchmark test functions [35]. The following algorithm parameters are chosen while executing the hybrid algorithms: the search agents are 30, the maximum iterations are 500, and the iteration is 30. Table 2 shows the statistical results for several unimodal (f1 to f7) and multimodal (f8 to f10) functions. To authenticate the supremacy of the said technique, GSA [26], MWOA [28], ASO [29], AVOA [31] ,and a published hybrid hGGSA-PS algorithm [37] results are also given in Table 2. It can be seen that, for almost all test functions, the proposed hybrid technique is very competitive compared to other algorithms.
The computing time for 30 runs for all the benchmark capabilities for the SCaAVOA and AVOA computations is also shown in Table 3. Table 3 shows that the proposed SCaAVOA tends to take less time to compute than the initial AVO calculation for each of the capabilities. According to the study mentioned above, the modified algorithm performs better than the normal calculation; hence, we will use it for our proposed study.

7.2. Statistical Test of the Proposed Sine Cosine-Adopted AVOA

A nonparametric statistical test known as Wilcoxon’s test is used here to test the significance of the proposed SCaAVOA procedure with that of the standard AVOA methods. The outcomes (p esteem) are displayed in Table 4 for pairwise correlation of the SCaAVOA vs. AVOA. An algorithm is significant when the p esteem returns a worth under 0.05. It can be seen from Table 3 that, for most of the functions, the p esteem is noticeably under 0.05. This suggests that the SCaAVOA is significant in accomplishing the best encouraging answers for the standard benchmark functions [38].

7.3. Implementation of the Proposed SCaAVOA Algorithm

When the simulation is complete, the limits of the TID, PID, and PI regulators can be found by using Equation (18) to calculate the goal function for the five-area thermal system. Table 4 displays the regulator boundaries that contrast the current tactic with other conventional methods. Table 4 also shows that, when compared to the AVOA-based TID and the AVOA-based PI regulator, the Sine Cosine-adopted AVOA-based TID controller produces the best results. In a similar manner, it shows that, overall, the proposed plan with the proposed Sine Cosine-adopted AVOA-based TID system gives improved results when compared to the AVOA.

7.4. Testing on Five-Area Thermal Power System

The system response is shown in Table 5 for a 1% step load expansion in area 1 at t = 0 s. The results using the PID controllers are also shown in Table 5 and Figure 4 for comparison (A–C). It is clear from Figure 4 that the proposed SCaAVOA-based TID structure achieves improved system response when compared to the simple PID structure. Additionally, Table 5 provides the performance list for the said scenarios. The table shows that the proposed TID structure performs better than the traditional PID structure in terms of lower ITAE values and faster settling times for Δf and Ptie.

7.5. Extension to Five-Area, Ten-Unit Hybrid System

To show the capacity of the proposed plan to deal with systems using various sources of generations, a five-area, ten-unit system is examined. The hybrid system is composed of a non-reheat type of thermal power system in each area and a photovoltaic system in area 1, a wind system in area 2, an aqua electrolyzer and fuel cell in area 3, a diesel engine generator and microturbine in area 4, and a plug-in electric-vehicle in area 5, respectively. The TID structures are accepted for every area and the proposed Sine Cosine-adopted AVOA calculation is utilized to tune the regulator boundary. The following disturbance is taken into consideration:

7.5.1. Case 1: Solar Variations

For the primary occurrence, a test is performed by thinking about adjustment to solar irradiation, as shown in Figure 5A. Table 6 shows the regulator boundaries for the said disturbance. Figure 5B–D outline the frequency response bend of the hybrid power system for the TID-based planned regulators. A conclusion can be made that the TID-based hybrid power system offers better framework execution when it is contrasted with other standard regulators.

7.5.2. Case 2: Wind Penetration Variations

At the same time, one more certified test is performed by an adjustment of wind power variation, as shown in Figure 6A. Figure 6B–D outline the system reaction showing the examination of the proposed and the regular regulators. A similar conclusion can be drawn as in case 1.

7.5.3. Case 3: Simultaneous Variation of All Disturbance

For this situation, a simultaneous change in the different system parameters, including sun irradiation change, wind power change, and a step change of 1% load in locale 1 at t = 0 s, is considered. Figure 7A–F exhibit the frequency reaction of the AGC five-area power system and a conclusion can be made that, for each disturbance, the system frequency maintains a steady state of equilibrium.

7.5.4. Case 4: Stability Test of the System in Time Domain and Frequency Domain

The stability test of the system can be achieved by analyzing the circuit in time domain as well as frequency domain. The time response and the frequency response are calculated by taking the bode plot and the root locus of the hybrid power system, as shown in Figure 8A,B. Even though the system is experiencing all the disturbances, Figure 8A,B demonstrate a stable system using the SCaAVOA-adjusted TID controller. As a result, the model has been linearized, and the response has been obtained by inputting data at the controller’s output and measuring the result at the power system’s output. From the figure, it can be seen that the system is stable with the proposed TID controller.

8. Conclusions

This study proposes a novel Sine Cosine-adopted African vulture optimization algorithm (SCaAVOA) for a tilt-integral-derivative (TID) regulator structure for frequency control of a hybrid five-region power system. The hybrid system considers different conventional and non-conventional sources along with some non-linearities, such as GRC and dead band. A comparison of the proposed calculation over the standard calculation in terms of simulation time and fitness function is examined initially. Furthermore, the SCaAVOA approach is then used to additionally investigate TID regulator limits for a five-area thermal system as well as a five-area, ten-unit hybrid power system for the frequency guideline. It can be seen from the comparison table that the rate of improvement in error (J) by the application of the SCA-AVOA-based TID controller in contrast to the AVOA-based TID, the AVOA-based PID, and the AVOA-based PI are 3.79%, 28.38%, and 50%, respectively.
It can be seen from the correlation table that there is a high rate of decrease in the performance indices, i.e., in the J values. The simulation results show that the proposed SCaAVOA-based TID controller is useful for frequency regulation, when compared to the standard controller. In the present work, for each area, only thermal system modeling is included and some other system can be tested in the AGC, which can be considered in a future study of the proposed work. Future work will be focused on developing more precise control in the RES-dominated AGC. Furthermore, the authors will try to incorporate more detailed models of power system devices and different modes of operation.

Author Contributions

Conceptualization, S.R.N., R.K.K. and T.S.U.; methodology, S.R.N., R.K.K. and T.S.U.; validation, S.R.N., R.K.K.; investigation, S.R.N., R.K.K., S.P. (Sidhartha Panda) and T.S.U.; writing—original draft preparation, S.P. (Sasmita Padhy) and P.R.S.; writing—review and editing, R.K.K. and T.S.U.; visualization, S.R.N. and R.K.K.; Project administration, R.K.K. and T.S.U.; funding acquisition, T.S.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elgerd, O.I. Electric Energy Systems Theory; Tata McGraw Hill: New Delhi, India, 2006. [Google Scholar]
  2. Ranjan, M.; Shankar, R. A literature survey on load frequency control considering renewable energy integration in power system, Recent trends and future prospects. J. Energy Storage 2022, 45, 103717. [Google Scholar] [CrossRef]
  3. Khadanga, R.K.; Kumar, A. Hybrid adaptive ‘gbest’-guided gravitational search and pattern search algorithm for automatic generation control of multi-area power system. IET Gener. Transm. Distrib. 2016, 11, 3257–3267. [Google Scholar] [CrossRef]
  4. Sahu, B.K.; Mohanty, P.K. Design and implementation of Fuzzy-PID Controller with Derivative Filter for AGC of two-area interconnected Hybrid Power System. Int. J. Innov. Technol. Explor. Eng. 2019, 8, 4198–4212. [Google Scholar] [CrossRef]
  5. Khadanga, R.K.; Kumar, A.; Panda, S. A hybrid shuffled frog-leaping and pattern search algorithm for load frequency controller design of a two-area system composing of PV grid and thermal generator. Int. J. Numer. Model. Electron. Netw. Devices Fields 2020, 33, 2694. [Google Scholar] [CrossRef]
  6. Demiroren, H.L.; Zeynelgil, G.A. application to optimization of AGC in three-area power system after deregulation. Int. J. Electr. Power Energy Syst. 2007, 29, 230–240. [Google Scholar] [CrossRef]
  7. Sinha, S.K.; Patel, R.N.; Prasad, R. Application of GA and PSO tuned fuzzy controller for AGC of three area thermal-thermal-hydropower system. Int. J. Comput. Theory Eng. 2010, 2, 238. [Google Scholar]
  8. Latif, A.; Das, D.C.; Barik, A.K.; Ranjan, S. Maiden coordinated load frequency control strategy for ST-AWEC-GEC-BDDG-based independent three-area interconnected microgrid system with the combined effect of diverse energy storage and DC link using BOA-optimised PFOID controller. IET Ren. Power Gener. 2019, 13, 2634–2646. [Google Scholar] [CrossRef]
  9. Pan, I.; Das, S. Fractional order AGC for distributed energy resources using robust optimization. IEEE Trans. Smart Grid 2015, 7, 2175–2186. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, D.J.; Wang, L. Small-signal stability analysis of an autonomous hybrid renewable energy power generation/energy storage system part I: Time-domain simulations. IEEE Trans. Energy Convers. 2008, 23, 311–320. [Google Scholar] [CrossRef]
  11. Das, D.C.; Roy, A.K.; Sinha, N. GA based frequency controller for solar thermal–diesel–wind hybrid energy generation/energy storage system. Int. J. Electr. Power Energy Syst. 2012, 43, 262–279. [Google Scholar] [CrossRef]
  12. Rajesh, K.S.; Dash, S.S.; Rajagopal, R. Hybrid improved firefly-pattern search optimized fuzzy aided PID controller for automatic generation control of power systems with multi-type generations. Swarm Evol. Comput. 2019, 44, 200–211. [Google Scholar] [CrossRef]
  13. 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]
  14. Feng, H.; Ma, W.; Yin, C.; Cao, D. Trajectory control of electro-hydraulic position servo system using improved PSO-PID controller. Autom. Constr. 2021, 127, 103722. [Google Scholar] [CrossRef]
  15. Guo, B.; Zhongjie, Z.; Pan, J.S.; Chu, S.C. Optimal design and simulation for PID controller using fractional-order fish migration optimization algorithm. IEEE Access 2021, 9, 8808–8819. [Google Scholar] [CrossRef]
  16. Rai, A.; Das, D.K. The development of a fuzzy tilt integral derivative controller based on the sailfish optimizer to solve load frequency control in a microgrid. incorporating energy storage systems. J. Energy Storage 2022, 48, 103887. [Google Scholar] [CrossRef]
  17. 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]
  18. Katoch, S.; Chauhan, S.S.; Kumar, V. A review on genetic algorithm: Past, present, and future. Multimed. Tools Appl. 2021, 80, 8091–8126. [Google Scholar] [CrossRef]
  19. Latif, A.; Paul, M.; Das, D.C.; Hussain, S.M.S.; Ustun, T.S. Price Based Demand Response for Optimal Frequency Stabilization in ORC Solar Thermal Based Isolated Hybrid Microgrid under Salp Swarm Technique. Electronics 2020, 9, 2209. [Google Scholar] [CrossRef]
  20. Singh, S.; Chauhan, P.; Aftab, M.A.; Ali, I.; Hussain, S.M.S.; Ustun, T.S. Cost Optimization of a Stand-Alone Hybrid Energy System with Fuel Cell and PV. Energies 2020, 13, 1295. [Google Scholar] [CrossRef] [Green Version]
  21. Abdolrasol, M.G.; Hannan, M.; Hussain, S.S.; Ustun, T.S. Optimal PI controller based PSO optimization for PV inverter using SPWM techniques. Energy Rep. 2021, 8, 1003–1011. [Google Scholar] [CrossRef]
  22. Dey, P.P.; Das, D.C.; Latif, A.; Hussain, S.M.S.; Ustun, T.S. Active Power Management of Virtual Power Plant under Penetration of Central Receiver Solar Thermal-Wind Using Butterfly Optimization Technique. Sustainability 2020, 12, 6979. [Google Scholar] [CrossRef]
  23. Chauhan, A.; Upadhyay, S.; Khan, M.T.; Hussain, S.M.S.; Ustun, T.S. Performance Investigation of a Solar Photovoltaic/Diesel Generator Based Hybrid System with Cycle Charging Strategy Using BBO Algorithm. Sustainability 2021, 13, 8048. [Google Scholar] [CrossRef]
  24. Farooq, Z.; Rahman, A.; Hussain, S.M.S.; Ustun, T.S. Power Generation Control of Renewable Energy Based Hybrid Deregulated Power System. Energies 2022, 15, 517. [Google Scholar] [CrossRef]
  25. Latif, A.; Hussain, S.M.S.; Das, D.C.; Ustun, T.S. Double stage controller optimization for load frequency stabilization in hybrid wind-ocean wave energy based maritime microgrid system. Appl. Energy 2021, 282, 116–171. [Google Scholar] [CrossRef]
  26. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A gravitational search algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  27. Tabak, A.; Duman, S. Levy Flight and Fitness Distance Balance-Based Coyote Optimization Algorithm for Effective Automatic Generation Control of PV-Based Multi-Area Power Systems. Arab. J. Sci. Eng. 2022, 47, 14757–14788. [Google Scholar] [CrossRef]
  28. Khadanga, R.K.; Kumar, A.; Panda, S. A novel modified whale optimization algorithm for load frequency controller design of a two-area power system composing of PV grid and thermal generator. Neural Comput. Appl. 2020, 32, 8205–8216. [Google Scholar] [CrossRef]
  29. Safiullah, S.; Rahman, A.; Lone, S.A.; Hussain, S.M.S.; Ustun, T.S. Novel COVID-19 Based Optimization Algorithm (C-19BOA) for Performance Improvement of Power Systems. Sustainability 2022, 14, 14287. [Google Scholar] [CrossRef]
  30. Khadanga, R.K.; Kumar, A.; Panda, S. Frequency control in hybrid distributed power systems via type-2 fuzzy PID controller. IET Renew. Power Gener. 2021, 15, 1706–1723. [Google Scholar] [CrossRef]
  31. Abdollahzadeh, B.; Gharehchopogh, F.S.; Mirjalili, S. African vultures optimization algorithm: A new nature-inspired metaheuristic algorithm for global optimization problems. Comput. Ind. Eng. 2021, 158, 107408. [Google Scholar] [CrossRef]
  32. Kumar, C.; Mary, D.M. Parameter estimation of three-diode solar photovoltaic model using an Improved-African Vultures optimization algorithm with Newton–Raphson method. J. Comput. Electron. 2021, 20, 2563–2593. [Google Scholar] [CrossRef]
  33. Das, A.; Dawn, S.; Gope, S.; Ustun, T.S. A Strategy for System Risk Mitigation Using FACTS Devices in a Wind Incorporated Competitive Power System. Sustainability 2022, 14, 8069. [Google Scholar] [CrossRef]
  34. Hakimuddin, N.; Nasiruddin, I.; Bhatti, T.S.; Arya, Y. Optimal automatic generation control with hydro, thermal, gas, and wind power plants in 2-area interconnected power system. Electr. Power Compon. Syst. 2020, 48, 558–571. [Google Scholar] [CrossRef]
  35. Khadanga, R.K.; Kumar, A.; Panda, S. A modified Grey Wolf Optimization with Cuckoo Search Algorithm for load frequency controller design of hybrid power system. Appl. Soft Comput. 2022, 124, 109011. [Google Scholar] [CrossRef]
  36. Deng, Z.; Xu, C.; Huo, Z.; Han, X.; Xue, F. Sliding Mode Based Load Frequency Control and Power Smoothing of Power Systems with Wind and BESS Penetration. Machines 2022, 10, 1225. [Google Scholar] [CrossRef]
  37. Khadanga, R.K.; Nayak, S.R.; Panda, S.; Das, D.; Prusty, B.R.; Sahu, P.R. A Novel Optimal Robust Design Method for Frequency Regulation of Three-Area Hybrid Power System Utilizing Honey Badger Algorithm. Int. Trans. Electr. Energy Syst. 2022, 2022, 6017066. [Google Scholar] [CrossRef]
  38. Mohanty, D.; Panda, S. Modified salp swarm algorithm-optimized fractional-order adaptive fuzzy PID controller for frequency regulation of hybrid power system with electric vehicle. J. Control. Autom. Electr. Syst. 2021, 32, 416–438. [Google Scholar] [CrossRef]
Figure 1. One area of an interconnected power system.
Figure 1. One area of an interconnected power system.
Energies 16 00926 g001
Figure 2. Proposed structure of a five-area hybrid power system.
Figure 2. Proposed structure of a five-area hybrid power system.
Energies 16 00926 g002
Figure 3. Structure of a TID controller.
Figure 3. Structure of a TID controller.
Energies 16 00926 g003
Figure 4. (AC) System response for the five-area power system.
Figure 4. (AC) System response for the five-area power system.
Energies 16 00926 g004
Figure 5. (AD) System response for the five-area, ten-unit system with solar disturbance.
Figure 5. (AD) System response for the five-area, ten-unit system with solar disturbance.
Energies 16 00926 g005aEnergies 16 00926 g005b
Figure 6. (AD) System response for the five-area, ten-unit system with wind disturbance.
Figure 6. (AD) System response for the five-area, ten-unit system with wind disturbance.
Energies 16 00926 g006aEnergies 16 00926 g006b
Figure 7. (AF) System response for the five-area, ten-unit power system with all disturbances.
Figure 7. (AF) System response for the five-area, ten-unit power system with all disturbances.
Energies 16 00926 g007aEnergies 16 00926 g007b
Figure 8. (A,B) Bode and root locus responses of the hybrid power system.
Figure 8. (A,B) Bode and root locus responses of the hybrid power system.
Energies 16 00926 g008
Table 1. Parameters used for the simulation of the power plant model under study.
Table 1. Parameters used for the simulation of the power plant model under study.
ComponentsGain (K)Time Constant (T)
Wind turbine Generator (WTG)KWTG = 1TWTG = 1.5
Hydro-Aqua Electrolyzer (AE)KHAE = 0.002THAE = 0.5
Fuel Cell (FC)KFC = 0.01TFC = 4
Battery Energy Storage System (BESS)KBESS = −0.01TBESS = −0.1
Diesel Energy Storage System (DEG)KDEG = 0.003TDEG = 2
Micro-TurbineKMTG = 1TMTG = 1.5
Thermal Power SystemTg = 0.08, Tt = 0.3, T12 = 0.0866, B = 0.425, Kr = 0.5, Tr = 10.0
Damping CoefficientD = 0.03
Inertia ConstantM = 0.4
Table 2. Cost function comparison of the proposed SCaAVOA considering standard benchmark functions.
Table 2. Cost function comparison of the proposed SCaAVOA considering standard benchmark functions.
FunctionDE [20]PSO [21]GSA [26] MWOA [28]ASO [29]hGGSA-PS [34]AVOA [31]Proposed SCaAVOA
Avg.Std. DevAvg.Std. DevAvg.St. DevMeanSt. DevMeanSt. DevMeanSt. DevMeanSt. DevMeanSt. Dev
f18 × 10−145.9 × 10−140.000140.0002023.985 × 10−86.671 × 10−182.986 × 10−106.21 × 10−392.687 × 10−213.652 × 10−212.795 × 10−99.63 × 10−100.696 × 10−293000
f21.5 × 10−99.9 × 10−100.042140.0454214.196 × 1032.601 × 10−92.123 × 1032.98 × 10−253.332 × 10−101.892 × 10−103.978 × 1033.69 × 1020.028 × 10−1520.153 × 10−1520.113 × 10−1640.445 × 10−164
f36.8 × 10−117.4 × 10−1170.125622.119243.567 × 10−1441.966434.653 × 10−153.78 × 10−101.983 × 10−279.70243.567 × 10−146.39 × 10−190.421 × 10−22300.859 × 10−2490
f4001.086480.3170392.964× 10257.531 × 10−96.523 × 10−280.025143.245 × 10−96.142 × 10−98.567 × 10−273.25 × 10−320.145 × 10−1360.794 × 10−1360.406 × 10−1590.221 × 10−159
f50096.718360.115591.67411.298360.0591.23 × 10−1024.83880.5158530.9010.00250.86744.75080.2891× 10−30.447 × 10−3
f6000.00018.28 × 10−59.865 × 10−405.023 × 10−60.00219006.758 × 10−64.36 × 10−190.526 × 10−60.328 × 10−60.332 × 10−50.3101 × 10−5
f70.004630.00120.122850.0449572.4860.0077161.9562.12 × 10−50.0356410.0194982.4860.0077160.182 × 10−30.256 × 10−30.153 × 10−30.192 × 10−3
f8−11080.1574.74841.291152.8148.456 × 10−27341.60064.565 × 10−283.65 × 10−30−7428.17422.39776.787 × 10−272.65 × 10−351.239 × 10−40.036 × 10+41.251 × 10−40.016 × 10+4
f969.238.846.704211.629388.986 × 10−114.4398620.021 × 10−120006.787 × 10−114.25 × 10−170000
f109.7 × 10−84.2 × 10−80.276020.509012.5463.961 × 10−100.896 × 10−108.36 × 10−153.001 × 10−112.151 × 10−110.9010.002560.888× 10−1500.888× 10−150
Table 3. Comparison of computational time between the SCaAVOA and AVOA algorithms for the considered standard benchmark functions (f1f10).
Table 3. Comparison of computational time between the SCaAVOA and AVOA algorithms for the considered standard benchmark functions (f1f10).
Functions% Reduction in Execution TimeSCaAVOA
Elapsed Time (s)
AVOA
Elapsed Time (s)
f110.0313.4593.845
f26.2353.5493.785
f37.9457.6128.269
f46.8353.4483.701
f52.2473.9154.005
f68.0253.3123.601
f76.1015.4945.851
f85.0864.0124.227
f97.0163.5253.791
f107.7633.5643.864
Table 4. p test for pairwise comparison of the SCaAVOA vs. AVOA algorithm.
Table 4. p test for pairwise comparison of the SCaAVOA vs. AVOA algorithm.
SCA-AOA vs. AOA
FunctionsValuesFunctionsValues
f12.213 × 10−6f61.2118 × 10−12
f21.2118 × 10−12f71.2118 × 10−12
f31.2118 × 10−12f81.2118 × 10−12
f41.2118 × 10−12f9NA
f51.2118 × 10−12f101.6853 × 10−14
Table 5. Controller parameters for the five-area power system.
Table 5. Controller parameters for the five-area power system.
Component ParameterProposed SCaAVOA Tuned
TID Controller
AVOA Tuned
TID
AVOA Tuned
PID
AVOA Tuned
PI
Controller-1KP1 = 1.7995, KI1 = 1.9958, KD1 = 1.0127, n1 = 2.5109KP1 = 1.9954, KI1 = 1.9011, KD1 = 0.9092, n1 = 2.6816;KP1 = 1.9372, KI1 = 1.9958, KD1 = 1.0350KP1 = 0.0089, KI1 = 0.2184
Controller-2KP2 = 1.9952, KI2 = 1.9956, KD2 = 1.9942, n2 = 9.9192KP2 = 1.9317, KI2 = 1.9030, KD2 = 1.9793, n2 = 9.9800;KP2 = 1.9952, KI2 = 1.6493, KD2 = 1.9878KP2 = 0.0011, KI2 = 0.2516
Controller-3KP3 = 1.9946, KI3 = 1.9938, KD3 = 1.9858, n3 = 9.5630KP3 = 1.4228, KI3 = 1.9938, KD3 = 1.4069, n3 = 5.1143;KP3 = 1.9946, KI3 = 1.8576, KD3 = 1.9858KP3 = 0.1462, KI3 = 0.2022
Controller-4KP4 = 1.9946, KI4 = 1.9938, KD4 = 1.9558, n4 = 4.2015KP4 = 1.9946, KI4 = 1.9938, KD4 = 1.9558, n4 = 8.0693;KP4 = 1.9186, KI4 = 1.9938, KD4 = 1.9858KP4 = 0.0921, KI4 = 0.2291
Controller-5KP5 = 1.9946, KI5 = 1.9938, KD5 = 1.9858, n5 = 9.9192KP5 = 1.9575, KI5 = 1.9938, KD5 = 1.2533, n5 = 9.9800;KP5 = 1.9104, KI5 = 1.9938, KD5 = 1.6359KP5 = 0.2236, KI5 = 0.2227
ITAE0.10010.10450.16512.795
Table 6. Controller parameters for the five-area, ten-unit hybrid power system.
Table 6. Controller parameters for the five-area, ten-unit hybrid power system.
Component ParameterProposed SCaAVOA Tuned
TID Controller
AVOA Tuned
TID
AVOA Tuned
PID
AVOA Tuned
PI
Controller-1KP1 = 1.1381, KI1 = 0.3380, KD1 = 0.8155, n1 = 3.0131 KP1 = 0.7316, KI1 = 0.3577, KD1 = 0.2553, n1 = 2.3491KP1 = 0.7935, KI1 = 0.5885, KD1 = 0.4888KP1 = 0.0089, KI1 = 0.2184
Controller-2KP2 = 1.8385, KI2 = 0.1870, KD2 = 1.3267, n2 = 5.1948 KP2 = 0.8843, KI2 = 0.2229, KD2 = 0.3843, n2 = 3.6644KP2 = 1.9952, KI2 = 0.1775, KD2 = 0.4528KP2 = 0.0011, KI2 = 0.2516
Controller-3KP3 = 1.2227, KI3 = 0.9752, KD3 = 1.3722, n3 = 8.6644KP3 = 0.9090, KI3 = 0.575, KD3 = 0.7078, n3 = 1.4077KP3 = 0.2463, KI3 = 1.9931, KD3 = 1.9858KP3 = 0.01462, KI3 = 0.2022
Controller-4KP4 = 1.5446, KI4 = 0.5303, KD4 = 0.9427, n4 = 9.9800 KP4 = 0.7742, KI4 = 0.3660, KD4 = 0.6315, n4 = 3.5428KP4 = 0.5474, KI4 = 1.8711, KD4 = 1.9815KP4 = 0.0092, KI4 = 0.0229
Controller-5KP5 = 1.1889, KI5 = 0.6183, KD5 = 0.001, n5 = 3.0732KP5 = 1.1007, KI5 = 0.1775, KD5 = 0.7676, n5 = 3.2820KP5 = 0.5022, KI5 = 1.9576, KD5 = 0.3243KP5 = 0.0023, KI5 = 0.0022
ITAE18.9919.7224.3837.69
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nayak, S.R.; Khadanga, R.K.; Panda, S.; Sahu, P.R.; Padhy, S.; Ustun, T.S. Participation of Renewable Energy Sources in the Frequency Regulation Issues of a Five-Area Hybrid Power System Utilizing a Sine Cosine-Adopted African Vulture Optimization Algorithm. Energies 2023, 16, 926. https://doi.org/10.3390/en16020926

AMA Style

Nayak SR, Khadanga RK, Panda S, Sahu PR, Padhy S, Ustun TS. Participation of Renewable Energy Sources in the Frequency Regulation Issues of a Five-Area Hybrid Power System Utilizing a Sine Cosine-Adopted African Vulture Optimization Algorithm. Energies. 2023; 16(2):926. https://doi.org/10.3390/en16020926

Chicago/Turabian Style

Nayak, Smruti Ranjan, Rajendra Kumar Khadanga, Sidhartha Panda, Preeti Ranjan Sahu, Sasmita Padhy, and Taha Selim Ustun. 2023. "Participation of Renewable Energy Sources in the Frequency Regulation Issues of a Five-Area Hybrid Power System Utilizing a Sine Cosine-Adopted African Vulture Optimization Algorithm" Energies 16, no. 2: 926. https://doi.org/10.3390/en16020926

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

Article Metrics

Back to TopTop