Next Article in Journal
Review on Gassing Tendency of Different Insulating Fluids towards Transformer Applications
Next Article in Special Issue
Review of Voltage-Bucking/Boosting Techniques, Topologies, and Applications
Previous Article in Journal
Public Policies for Renewable Energy: A Review of the Perspectives for a Circular Economy
Previous Article in Special Issue
Modeling of Memristors under Periodic Signals of Different Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of a 2DOF-PID Control Scheme for Frequency/Power Regulation in a Two-Area Power System Using Dragonfly Algorithm with Integral-Based Weighted Goal Objective

by
Alaa M. Abdel-hamed
1,
Almoataz Y. Abdelaziz
2 and
Adel El-Shahat
3,*
1
Electrical Power and Machines Department, High Institute of Engineering, El-Shorouk Academy, Cairo 11837, Egypt
2
Faculty of Engineering and Technology, Future University in Egypt, Cairo 11835, Egypt
3
Energy Technology Program, School of Engineering Technology, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
Energies 2023, 16(1), 486; https://doi.org/10.3390/en16010486
Submission received: 17 November 2022 / Revised: 19 December 2022 / Accepted: 28 December 2022 / Published: 2 January 2023
(This article belongs to the Special Issue Trends and Prospects in Analysis and Control of Power Electronics)

Abstract

:
The increase in power demand, nonlinearity, complexity, varying structure, and other important causes has necessitated the implementation of artificial intelligent control methodologies for safe and acceptable operation of the electric power systems. Therefore, in this article, an improved two-degrees-of-freedom (2DOF-PID) control scheme is proposed for power/frequency control of a two-area interconnected electric power system. The parameters of the 2-DOF-PID control scheme are optimized using the Dragonfly Algorithm (DA) via a new integral-based weighted goal fitness function (IB-WGFF) (i.e., DF-2DOF-PID-IB-WGFF). The superiority of the suggested scheme is proved by comparing the results obtained using the proposed IB-WGFF with those obtained using the conventional controllers, and the 2DOF-PID controllers optimized using the DA and Genetic Algorithm (GA) via the frequently published performance criterion. To verify the stability, efficacy, and robustness of the proposed control scheme, a load disturbances and parameters perturbations with various percentages are implemented in the controlled system under the same controllers. Finally, verification results proved that the proposed 2DOF-PID optimized using DA via the IB-WGFF is more stable, efficient, and robust than the other controllers recently used in the literature.

1. Introduction

A modern Power System (PS) is a dynamic behavior system where interactive dynamic elements connected via power lines. Its function is to supply power to various loads or ensemble of consumers as a one from its critical various parts [1,2].
For a power system to operate in stability region, both the tie-line power deviation (ΔP) and system frequency (Δf) had to be constant. The Frequency variation effects badly the speed of electric devices in power system such as induction and synchronous motors. High magnetizing currents in induction motors and transformers may be appear from the drop in frequency. Controlling the frequency and tie-line power, to be at their normal rated values, is a complex and difficult task for the variable nature of the power demand. The frequency/power deviation causes various undesired harms such as instability problems. The power system stability and its dynamic performance can be enhanced successfully by controlling the frequency/power via accurate design of Automatic Generation Control (AGC) or Load Frequency Control (LFC) scheme. The AGC is needed to ensure synchronization and for maintaining the frequency and tie-line power of a power system within their nominal values during load variations [3,4].
Many control methodologies have been recommended for addressing the concerns of AGC such as optimal, adaptive, and Artificial Intelligent (AI) methods [5]. Standard PID controllers are frequently used to dress the power/frequency control for their simplicity, reliability, and the suitable ratio between cost and implementation performance [6]. The nonlinearity, uncertainties, and realistic aspects such as load disturbances of an interconnected power system increase instability problems in the power system [7]. The classical PID-controllers are implemented at the explicit or normal operating conditions. The nonlinearity characteristics in the interconnected real power systems makes the Conventional PID control not convenient in the fluctuating conditions, hence the PID controllers performance must be improved [8]. The artificial optimization search techniques are utilized efficiently for load frequency Control [9], deregulated LFC of multi-area power system via mine blast algorithm [10], and PID/FOPID-based Frequency control of zero-carbon multi-sources based interconnected power systems under deregulated scenarios [11]. Various improved control strategies for LFC have been suggested in the literature.
P. Praveena et al. in [12] presented an approach for selecting the parameters of the PID controller in LFC of a single area power system. They used the Ziegler-Nichols (ZN) and Genetic Algorithm (GA) for tuning the parameters of the suggested controller. Comparative study is presented using MATLAB Simulink. The simulation results indicated that the PID tuned via GA offers lower settling time and a zero stead-state error. No verification for the suggested controller are checked under load disturbances or parameters variation in the power system.
N. El Yakine et al. in [13] used a Particle Swarm Optimization (PSO) method for optimal tuning of the PID controllers used to control the LFC and the Automatic Voltage Regulator (AVR) in a two-area power system. The proposed control scheme controls separately, the reactive and active powers. The objective of their strategy is to suppress the system fluctuations resulting from disturbances and restore system’s voltage and frequency to their nominal values. The obtained results using the suggested control are compared to the results obtained using Z-N method, GA, and Bacterial Foraging (BF) optimization methods. The authors, in their paper, used a step load disturbance, but a verification using various disturbances and parameters perturbations are not verified.
Yogendra Arya in [14] introduced a new fuzyy-PID with filter and double integral controller (FPID-II) for AGC of a two-area power system. The scaling factors of the proposed controller is tuned implementing the imperialist competitive algorithm using ISE performance criterion. The authors compared the proposed method’s results with the results obtained using various Artificial Intelligent (AI) methods-based controllers. The scalability and efficacy of the proposed approach is proved by extending the study using a non-reheat thermal power system with nonlinearity and other different power system plants. The authors in this paper carried out a sensitivity analysis, under system parameters perturbations, to prove the robustness of the proposed approach, but no verification is carried out under variable disturbances and multiple load variations.
Devbrat Gupta et al. in [15] carried out a comparative study in integrated power system using Fractional Order (FO) ID, PID, and FOPI to control the power system under sudden disturbances in generations and loads. A GA is used to tune the controller gains. Simulation results proved that the FOPID controller is better in terms of the robustness and the frequency deviation.
Ruchika Lamba et al. in [16] proposed an interval FO-PID controller methodology for controlling the load frequency in a two-area interconnected power system. The stability boundary locus methodology is used to design the proposed controller. System parameters perturbation of (±50%) is checked using Kharitonov theorem instead of using fixed values. The turbine time-constant, governor constant, and load disturbance is changed while performing the analysis of the proposed controller. Simulation results illustrated that the proposed controller gives better performance than the conventional methods.
Rajveer Singh et al. in [17] used a 2DOF-PID controllers to control and analyze the frequency of two-area power system. The Month Moth Flame Optimization method (MFO) have been used to tune the parameters of the proposed controllers for minimizing the frequency deviation and the deviation in the tie-line power. The results obtained using the 2DOF-PID are compared with the results obtained using the recent meta-heuristic technique. Results proved the superiority of the suggested controllers over the other controllers in terms of the frequency deviation and tie-line power deviation. The sensitivity and efficiency of the proposed controllers against parameter variation and varying load disturbance still need to be verified.
Appala Naidu et al. in [18] introduced a new 2DOF-PID control scheme to control the frequency of an inter-connected power system integrated with wind power generation. The used controller is optimized using the Crow Optimization (CO) technique. In the proposed method, several reference inputs for a secondary control of frequency loop are used. Results showed that the supposed controller gives better control than other controllers.
EL Yakine, N. [19] introduced a new fuzzy-PID controller to solve the AGC problem in a two-area power system. The DA is implemented to tune the proposed controller parameters including the gains of the PID and the scaling factors of the fuzzy logic. The suggested methodology is extended to be applicable for AC/DC multisource power system linked with AC tie-line and a HVDC transmission link. The result of the proposed approach is compared with the results obtained using the AGC used in the literature. The authors extended their approach to be used in a non-linear power system of three areas of AGC. The comparative study proved the effectiveness of the suggested strategy in enhancing the frequency, hoverer the proposed scheme still needs tuning for the part associated with the shaping factor of the Fuzzy logic.
Dipayan Guha in [20] discussed the LFC associated with an interconnected power system implementing a new Quasi-Oppositional-Based (QO) JAYA. The used approach enhances the speed of convergence and select the best optimal solutions of the LFC problem by housing the knowledge of QO-based learning. The effectiveness of the suggested QOJAYA scheme is verified by implementing it to a plant of two area with multi-sources and multi areas. The PI controller is tuned and applied initially as secondary control. A 2DOF-PID controller is tuned and applied to enhance the dynamic stability of the power system under study. The proposed approaches used the single valued-time based objective as the fitness function to be minimized for searching the optimal setting for the proposed controllers. Simulation results demonstrated the superiority of the proposed QOJAYA controller over the base JAYA and Teaching Learning (TL) based tuning method. Sensitivity analysis is investigated to verify the suggested controller robustness.
Watanabe, M. et al. in [21] presented an improved method for regulating the frequency of a two-area power system. The proposed method implements the cascade configuration related to tilt integral-derivative and the FO-PID with a filter controller. Integrated squared Error (ISE), Integrated Absolute Error (IAE), Integrated Time Weighted Absolute Error (ITAE), and Integrated Time weighted Squared Error (ITSE) performance criteria are implemented in the optimization method as the objective function used to track the frequency deviation. The results are validated at changing loads and the robustness is checked against various uncertainties in the power system. The results proved stable operation under the proposed controller without no need for complex design methodologies.
Based on the above discussion, this paper proposes a 2DOF-PID control scheme for regulating the frequency and power fluctuations in a two-area power system. The Dragonfly Algorithm (DA) is used to optimally design the proposed scheme. A novel Integral-Based-weighted Goal Fitness Function (IB-WGFF) is implemented as the objective function to be minimized for controlling the frequency/power of the two-area power system. The proposed 2DOF-PID control scheme using the IB-WGFF is compared with the conventional- 2DOF-PID controller and the 2DOF-PID optimized using Dragonfly Algorithm (DA) and Genetic Algorithm (GA) via the other recently implemented fitness functions in the literature for interconnected power systems. The robustness for the designed 2DOF-PID is verified by investigating sensitivity analysis implementing loading and parameters uncertainties.
The paper contributes the following:
i.
Design and implementation a 2-DOF-PID control scheme for power/frequency control of two-area interconnected electric power system.
ii.
Design and implementation of a novel IB-WGFF to tune the 2DOF-PID controller parameters using the DA and GA.
iii.
Demonstrate the superiority of the proposed scheme by comparing the results with those obtained using the frequently published fitness-based controllers.
iv.
Verifying the stability, efficacy, and robustness of the proposed approach under load disturbances and parameter perturbations under the same proposed control scheme.
The remaining sections of this article are organized as follows: The model of the interconnected power system is presented in Section 2. Section 3 presents the mathematical modelling of the DA used for tuning the parameters of the proposed 2DOF-PID control scheme. Section 4 provides the design of the proposed fitness function and implementation of the DA and GA for tuning the parameters of the 2DOF-PID controllers. Section 5 introduces the simulation results and discussions. Verification of the proposed 2DOF-PID control scheme is introduced in Section 6. Finally, the conclusion is provided in Section 7.

2. Complete System Modeling

This section describes a complete system model of the interconnected power system presented in this paper. The two-area power system investigated here is commonly applied in the literature for designing and analyzing the LFC of the inter-connected electric power system. The system model of the two-area electric power system, used to discuss and analyze the limits and issues from the LFC, is shown in Figure 1a–c. The two areas are shown in Figure 1a,b shows a LFC implemented in a single area power system. In this system, each area encompasses a speed governor, steam turbine, and generator of 2000 MW rating, and a nominal step load of 2000 MW applied at the time instant 5 s [22,23].
The non-linear equations used to describe the system-dynamic behavior can be linearized during load changes. In this case, a linearized model can be applied in the investigations [24]. The steam turbine can be modelled using the transfer function given by Equation (1).
G t ( s ) = K t T t S + 1
where T t represents turbine time constant, K t represents the turbine gain.
When the controlling system identifies frequency fluctuations, which may cause system instability, it will adjust the speed of the turbine for achieving the optimum speed that eliminates the frequency oscillations. The governor is modelled using the transfer function indicated by Equation (2) [25].
G g ( s ) = K g T g S + 1
where T g represents governor time constant, K g represents the governor gain.
The control action donated by the governor mechanism is given by Equation (3).
Δ P c a = Δ P r = 1 R Δ f
where Δ P r is the set point change in power, Δ f is the frequency change, R is the parameter of the speed governor regulator.
The load and the generator can be modelled using the transfer function described by Equation (4) [26,27].
G S y s ( s ) = K p s T p s S + 1
where K p s = 1 D is the power system gain, T p s =   2 H f D is the time constant of the power system, D represents the damping, and H represents the system inertia.
For each area of the selected power system, there are three inputs (i.e., the tie-line power deviation; Δ P t i e , load perturbation; Δ P L , and controller input; Δ P r ) and two controlled outputs (i.e., the frequency deviation of the area 1; Δ P f 1 , and the frequency deviation of area 2; Δ P f 2 ). The summation of the deviation in the tie-line power and frequency is defined by the control error (CE) of each area. The CEs are considered to be the input to the secondary controllers (i.e., the designed 2DOF-PID controllers) corresponding to their areas. CEs, used for the studied system shown in Figure 1, are indicated by Equations (5) and (6).
C E 1 = Δ P t i e + B 1 Δ f 1
C E 2 = Δ P t i e + B 2 Δ f 2
At steady state, CEs are the controlling signals used to reduce the frequency and power fluctuation to zero. Therefore, the 2DOF-PID controller must be designed for enhancing the transient of the dynamic system.
The complete system model with the controllers is shown in Figure 1c. The parameters used for the different components of the system is given in Table 1. As mentioned in Section 1, the variation of voltage level for the investigated power system is negligible with variation in load, as a result, filter impedance and the tie-line impedance between the system sources are neglected.

2.1. DOF-PID Controller Modelling

The PID controller is still one of the most commonly used controllers in recent years due to its considered various advantages. These controllers work well even with uncertainty in the controlled systems. The mathematical model of the PID controllers is used as Equation (7) demonstrates: [14].
C ( s ) = Y ( s ) R ( s ) = K P + K I S + K D S
The two-Degree of Freedom (2DOF-PID) controllers are commonly used to detect and control disturbance and rapid shock in electric power systems devoid of symbolic boost in set point tracking overshoot. Implementing the 2DOF-PID controllers can reduce the reference signal effect on the controlled electric power systems [28].
Figure 2 demonstrates the control architecture for the 2-DOF-PID controller in the controlled systems. The mathematical model used for the 2-DOF-PID controller is indicated by Equation (8).
U ( s ) = K p ( P W . R ( s ) Y ( s ) ) + K i s ( R ( s ) Y ( s ) + K d s N . s + 1 ( D W . R ( s ) Y ( s ) )
where U(s) represents the controller output, R(s) and Y(s) represent the controller inputs, Kp is the proportional gain, Ki is the integral gain, Kd is the derivative gain, PW represents the reference weight of proportional gain, DW is the weight of the derivative, and N is the filter coefficient. In this paper, a FOPID control scheme is proposed for controlling the frequency/power of a nonlinear interconnected power system, the FOPID control scheme is designed and validated.

3. Dragonfly Algorithm (DA)

DA is proposed and designed by Mirjalili [29]. Dragonfly is mainly an optimization technique inspired by the dragonflies’ warming behavior which is static and dynamic. Its model is designed on the basis of the action, swarming dynamically or statically to find the foods and bypass enemies [30].
The static and dynamic swarm functions are considered to be two optimization stages which are called exploration and exploitation. In the exploration stage, phase dragonflies will build sub-swarms and will fly over several areas. This stage is the static swarms. In the exploitation phase, dragonflies will move in one strait direction for flying in larger swarms [31].
To model the flies’ motion, five characteristics of dragonflies are described here. These features are Separation (SI), Alignment (AI), Cohesion (CI), Attraction of Food (AFI), and Enemy (EI). The five steps or features can be described and modelled mathematically as follows:

3.1. Separation (SI)

For maximizing the population search space and avoiding collision, it is necessary for the distance between adjacent dragonflies to be within the specified or given locality. Assume I represents the number of individuals in a cluster and N represents the neighbors. The Separation (SI) is modelled using Equation (9).
S I = K = 1 K = N ( X X K )
where X is the dragonfly’s current positions, XK is Kth neighboring’s position of dragonfly

3.2. Alignment (AI)

The Alignment part (AI) matches the individual‘s velocity (VK) with the other dragonflies with parallel velocity. Equation (10) models for the Alignment feature of the Dragonfly optimization algorithm.
A I = K = 1 K = N V K N  
where VK is the velocity for the Kth neighboring dragonfly.

3.3. Cohesion (CI)

All dragonfly‘s individuals in a cluster are motivated for moving in the path of center-mass of the neighboring dragonflies. The cohesion feature step (CI) of dragonfly is modelled using Equation (11).
C I = K = 1 K = N V K N

3.4. Attraction of Food (AFI)

In this feature, all individuals have a tendency for moving in food direction. The attraction for food (AFI) feature at a location (XF) is developed using Equation (12).
A F I = X F X

3.5. Enemy (EI)

In this step, all individuals in cluster have a tendency for moving away from the enemy. The mathematical model of the enemy feature (EI) at the enemy location (XE) can be designed using Equation (13).
E I = X E + X
The dragonflies’ behavior in a cluster is affected by relating all the five features. The new updated position of the individual dragonflies indicated in Equation (14) is designed by step ( Δ X I ) indicated by Equation (15).
X I = X I + Δ X I
Δ X I = ( A . S I + B . A I + C . C I + D . A F I + E . E I ) + ω Δ X I
where A, B, and C represent the weighted coefficients of the separation, alignment, and cohesion, respectively, while D represents the food coefficients and E is the enemy factor. ω Represents the inertial weighted coefficient. The exploitative and explorative variation performances of the dragonflies can be realized and modelled utilizing different values of coefficients.

4. Proposed Control Objective and Optimal Tuning of the 2DOF-PID Controller

4.1. Proposed Objective Function Formulation

The frequency and power fluctuations of interconnected power system are greatly affected by the consumer requests and the system parameter perturbations. These fluctuations make the frequency ( Δ f a ) and the tie line power ( Δ P t i e ) deviate from the nominal performance values accepted for stable power system.
A satisfactory level of stability in the controlled system and a fast behavior with speedily dampened fluctuations are desired in LFC. These requirements make the system have the capability of returning the frequency deviations ( Δ f a ) of each area and tie line power deviation ( Δ P t i e ) to the original or predetermined values.
A proper fitness/objective (FF) function should be designed and implemented before using recent optimization-based control schemes for achieving the best power/frequency performance [32]. There are many fitness/objective functions which have been implemented in the LFC problems found in the literature review, such as ITSE, ISE, IAE, and ITAE. All these fitness/objective functions have been used for restoring the normal operation or performance.
The ITAE and ISE performance criteria are applied in literature due to better performance in comparison with the ITSE and IAE performance criteria. The ISE criterion integrates the squared error over the simulation time. As the squared error has much bigger values, the ISE criterion may penalize great errors beyond smaller errors. The ITAE performance criterion integrates the absolute error over the simulation time. The ITAE performance criterion will penalize errors and transients of long duration that occur later in the system responses [33,34,35].
Although The ITSE-optimized controllers minimized the contribution of initial errors and also assured errors later in the system response, they results of more time for settling in comparison with the ITAE-optimized controllers [6,33]. The ITAE-optimized controllers decreased the settling time better than the ISE, IAE, and ITSE-optimized controllers [6,33,34,35].
In the literature, the used fitness/objective function is the integral performance index type. Equations (16)–(19) demonstrates the ITSE, ISE, IAE and the ITAE fitness/objective functions used in the literature, respectively.
F F = I T S E = 0 t [ ( a = 1 z Δ f a ) 2 + ( Δ P t i e ) 2 ] . t . d t
F F = I S E = 0 t [ ( a = 1 z Δ f a ) 2 + ( Δ P t i e ) 2 ] .   d t
F F = I A E = 0 t [ a = 1 z | Δ f a | + | Δ P t i e | ] . d t
F F = I T A E = 0 t [ a = 1 z | Δ f a | + | Δ P t i e | ] . t . d t
where Δ f a represents the frequency deviation for area number a, Δ P t i e represents the tie-line power deviation, t represents the total simulation time of over which the power system is operated.
From the literature review, the ITAE-optimized controllers are desired in LFC systems. Thus, in this paper, a novel weighted goal method based on the ITAE performance criterion was designed and implemented. This method is used as the fitness function for tuning the parameters of the 2DOF-ID controller.

4.2. Proposed Integral-Based Weighted Goal Fitness Function (IB-WGFF)

In this paper, two objective performance criteria are proposed to be the fitness function used with the DA explained in Section 3 to design the 2DOF-PID modelled in Section 2.1 and indicated by Equation 8. The first objective is to minimize the integrated time-weighted absolute frequency deviations of the two-area power system ( Δ f 1 and Δ f 1 ). The second objective is to minimize the integrated time-weighted absolute power deviation (ΔP). The two objectives are combined as a single objective with weighting coefficients C1 and C2. This control objective will make the dynamic response of the power system more robust under variable loads and uncertainties in power system parameters.
Equation (20) shows the mathematical design of the proposed Integral-Based-weighted Goal Fitness Function (IB-WGFF) to be optimized using the DA algorithm for regulating the power and frequency fluctuations of a two-area interconnected power system. The fitness/objective function (FF) to be optimized is a weighted goal type of integral performance over the simulation time (50 s). To approve the performance of the proposed IB-WGFF under various performance criteria, the corresponding values of ITSE, ISE, ITAE and the IAE will also be implemented and calculated.
F F = I B W G F F = C 1 [ 0 t ( a = 1 z | Δ f a | ) t . d t ] + C 2 [ 0 t | Δ P t i e | . t . d t ]
where a is the number of interconnected area in the electrical power system; a = 1 z | Δ f a |   =   Δ f 1 + Δ f 2 ;   Δ f 1 is the frequency deviation of area1;   Δ f 2 is the frequency deviation for are2 of two-area interconnected power system; t is the total simulation time; C1 and C2 represents the weighting coefficients of the proposed fitness function. The parameters C1 and C2 are assumed to be in the range between [0, 1] such that C1 + C2 = 1. Increasing C1 and decreasing C2 in at same time makes the frequency deviation ( Δ f a ) responses better while increasing C1 and decreasing C2 makes the tie-line power dynamic’s response smoother and better. Varying C1 and C2 will affect the power system dynamic, and this effect will be investigated later in Section 5.2.
The main procedure of the DA using the proposed IB-WGFF, modelled in Equation (20) for optimizing the parameters of the suggested 2-DOF-PID control scheme, is shown in Figure 3. The upper (U) and the lower (L) control limits used for the control parameters are defined as shown in Equation (19).
K p L K p 1 ; K p 2   K p U K i L K i 1 ; K i 2   K i U K d L K d 1 ; K d 2   K d U P W L P W 1 ; P W 2   P W U D W L D W 1 ; D W 2   D W U
The lower values of the control parameters ( K p L , K i L , and K d L ) are chosen to be a value of zero, and their upper ( K p U , K i U , and K d U ) values are chosen to be a value of 1. The lower limits of ( P W L , D W L ) are chosen to be zero whereas their upper limits ( P W U , D W U ) are set to a value of 2.

4.3. Implementation of the Proposed DA-2DOF-PID Controller

The proposed 2DOF-PID control scheme is implemented in MATLAB software using 2017b operated on Intel core i5 2.3 GHz processor PC. Equation (20) represents the proposed IB-WGFF that is solved using the DA to control and improve the power/frequency of the two-area interconnected power systems. Equations (16)–(19) represent the Fitness/Objective functions in the literature used to compare and prove the effectiveness of the proposed 2DOF-PID control scheme. The implemented DA’s flow chart is displayed in Figure 4. According to many trials, the DA’s parameters implemented in the optimization process are indicated in Table 2. The DA has been run for 15 independent trials using different settings till the solutions are tangible close to each other.

4.4. Implementation of the GA-2DOF-PID Controller

In this paper, the GA-2DOF-PID controller is assumed for comparison. The parameters of the 2DOF-PID-controller are tuned by GA using the proposed IB-WGFF indicated in Equation (18) and the performance criteria in the literature which indicated in Equations (14)–(17). The initial population size of chromosomes that signifies the 2DOF-PIDb controller parameter is randomly generated. Then, the chromosomes are calculated according to the selected fitness/objective function defined by Equations (14)–(18) is Section 4.1 and Section 4.2. A new population that is based on the fitness/objective of individuals is selected. Genetic mutation and crossover are implemented in members of the population for creating the new solutions. The process of creation and evaluation the new population is continued until termination criteria has been satisfied [36,37]. According to experience, mutation should be completed with a low probability that ranges from 0.1 to 2%, whereas crossover should be between 60 to 90%. The GA has been run for 15 independent trials using different settings till the solutions are tangible close to each other. Again, the DA-2DOF-PID controllers are compared with the 2DOF-PID controllers tuned by GA as discussed in [37]. According to trials, the basic parameters of GA are given in Table 3. The DA-2DOF-PID controllers are also compared to the classical 2DOF-PID controller.

5. Results and Discussion

This section is dedicated for investigating the performance of the proposed control scheme in controlling and suppressing the frequency/power in a two- area interconnected power system. In this section, the parameters of the proposed 2DOF-PID are designed under the conditions of a step dynamic load implemented at the time instant 5 s. Five parts are designed and investigated in this section of the paper. The first part introduces the design of the 2DOF-PID controller using the classic method (i.e., the Z-N tuning method). The second parts introduce design and investigation of the 2DOF-PID scheme using DA via the new IB-WGFF indicated in Equation (20). Section three is designed to investigate the controlled system dynamic using the DA via the Fitness/objective functions used in the literature and indicated by Equations (16)–(19). The fourth section introduces a comparison between the proposed control scheme and other methods used in the literature. Finally, the fifth section introduces a comparison with the results obtained using the GA. The complete model using MALAB Simulink of the two-area power system with the 2DOF PID controllers is shown in Figure 5.

5.1. System Dynamic using 2-DOF-PID Controllers Tuned Conventionally

In this sub-section, the two-area power system shown in Figure 5 is used to investigate the frequency and power deviations. In this system, each area encompasses a speed governor, steam turbine, and generator of 2000 MW rating, and a nominal step load of 2000 MW applied at the time instant 5 s. The ZN tuning method used to design the parameters of the 2DOF-PID-controller is completed using the MATLAB® toolbox software (R2021a). The classical design method is supposed for optimizing the 2DOF-PID controller under a normal step load of 1 pu. is applied at the time instant of 5 s. The The parameters of the 2DOF-PID controllers, using Z-N method, and the dynamic response parameters of the controlled two-area interconnected power system are demonstrated in Table 4. Figure 6 and Figure 7 indicate the frequency deviation (Δf1) for area 1, frequency deviation (Δf2) for area 2 of the controlled interconnected power system, respectively. Figure 8 demonstrates the power deviation (ΔP) of the system using conventional 2FOPID controller.
It is evident from Figure 6 and Figure 7 that the two area have oscillating frequency over a large period of time besides the multiple and large values of undershoots, thus the parameters of the 2DOF-PID controllers still need to be accurately tuned using high performance objectives and approaches. It is evident also from Figure 8 that the power deviation has multiple and very large overshoots and long settling time and therefore, the parameters of the controllers have to be optimized using powerful approach.

5.2. Results using DA with the Proposed IB-WGFF

In this sub-section of the simulation results, the parameters of the 2DOF-PID control scheme, under step load disturbance at the time instant 5 s, are designed using the DA via the proposed IB-WGFF indicated in Equation (20). The proposed IB-WGFF is implemented by trying various weighting coefficient C1 and C2 such that C1 and C2 are between [0, 1] and C1 + C2 = 1. The efforts of the two-area inter-connected power system, shown in Figure 5, optimized using the DA via the proposed IB-WGFF, are shown in Table 4. The time responses of the frequency deviations (i.e., Δf2 and Δf2) for the two areas and the tie-line power deviation (ΔP) are shown in Figure 9, Figure 10 and Figure 11, respectively.
From Figure 9 and Figure 10 and Table 4, it is evident that the frequency deviation responses (i.e., Δf2 and Δf2) of the controlled power system are improved with varying the weighting coefficients C1 and C2. In addition, it is concluded from Table 5 and Figure 11 that the time-response of the tie-line power deviation (ΔP) is improved. It is evident that the dynamic response of the controlled power system is improved in all cases denoted in Table 5. Figure 9, Figure 10 and Figure 11 and Table 5 show that the frequency and power deviations have optimum values using the proposed IB-WGFF in case of (C1 = 0.5 and C2 = 0.5). According to the priority, the designer engineer will select between different cases demonstrated in Table 5.

5.3. Results Implementing Different Performance Indices in Literature

In this sub-section of the article, the parameters of the 2FOF-PID controller is designed using the DA via the performance indices in the literature indicated by Equations (16)–(19). The parameters of the upper and lower boundaries of the 2DOF-PID controller are the same as mentioned previously in Section 4.2. The main variable parameters of the DA are set the same as indicated in Table 2. The DA has been run for 15 independent trials using different settings till the solutions are tangible close to each other. Table 6 demonstrates the efforts of the 2DOF-PID controllers optimized using DA via the ITAE, ITSE, IAE, and ISE used in the literature. The frequency deviations for the two-areas (i.e., Δf2 and Δf2) obtained by DA-2DOF-PID via the various fitness are indicated in Figure 12 and Figure 13. The tie-line power deviation (ΔP) is indicated in Figure 14.
It is evident from Figure 12 and Figure 13 that the frequency deviation of the two-area (i.e., Δf2 and Δf2) implementing the ITAE is the smoothest and the best in overshot, and settling time but, it has lager undershot compared with the other performance indices. It is evident from Figure 14 that the 2DOF-PID optimized using DA via the ITAE shows the lowest overshot, undershot, and settling time in case of the tie-line power.

5.4. Comparison with Different Controllers

To emphasize the superiority of the proposed 2DOF-PID optimized using DA via the proposed IB-WGFF indicated in Equation (20), the obtained results are compared with the results obtained using the conventional 2DOF-PID and the 2DOF-PID optimized using DA via the best performance criterion implemented in the literature and indicated in Equations (16)–(19). For a fair comparison, the DA has been run for 15 independent trials using different settings till the solutions are tangible and close to each other in each case. The same parameters and constraints indicated and discussed in Section 4.2 and Section 4.4 are used to design the 2DOF-PID control scheme via the different fitness functions. For the conventional 2-DOF-PID, the classical design method is applied for optimizing its parameters as discussed in Section 5.1. The same trials are implemented in the tuning process of the conventional 2DOF-PID controller.
The parameters of the 2DOF-PID control scheme optimized conventionally, using the DA via the ITAE, using the DA via the proposed IB-WGFF are indicated in Table 7. In addition, the dynamic responses of the frequency deviations for the two-areas (i.e., Δf2 and Δf2) and the tie-line power deviation (ΔP) for the conventional 2DOF-PID, the DA-2DOF-PID using the ITAE, and the DA-2-DOF-PID controller via the proposed IB-WGFF (for C1 = 0.5 and C2 = 0.5) are shown in Figure 15, Figure 16 and Figure 17.
It is obvious, from Table 7 and Figure 15, Figure 16 and Figure 17, that the frequency deviation of the two-area (i.e., Δf2 and Δf2) and the tie-line power deviation (ΔP) of the proposed 2DOF-PID controller is the smoothest and the best in overshot, undershot, and settling time. The conventional 2DOF-PID shows larger undershoots and overshoots, and the performance of the conventional 2DOF-PID controller is worse than the DA-2DOF-PID-ITAE and the proposed 2DOF-PID controller. It is concluded, from the comparison, that the 2DOF-PID controller tuned using DA via proposed IB-WGFF (for C1 = 0.5 and C2 = 0.5) has been verified to be the best control scheme.

5.5. Comparison with GA

In this sub-section of the results, the advantage of the proposed 2DOF-PID control scheme optimized using DA via the IB-WGFF indicated in Equation (20) is emphasized in terms of performance. The GA in [37] with the parameters indicated and discussed in Section 4.2 and Section 4.4 is used to design the 2DOF-PID control scheme. The obtained results are compared with those obtained by the proposed 2DOF-PID control scheme and the classical 2DOF-PID controller.
The efforts of the system under study and the parameters of the 2DOF-PID controllers are shown in Table 8. Figure 18, Figure 19 and Figure 20 show the frequency deviations Δf1, Δf2, and power deviation (ΔP) of the system, respectively, using GA-2DOF-PID-ITAE, GA-2DOF-PID-IB-WGFF, and the classical 2DOF-PID. Figure 21, Figure 22 and Figure 23 show the frequency deviations (Δf1, Δf2) and power deviation (ΔP) of the system, respectively, using GA with the ITAE. Figure 24, Figure 25 and Figure 26 show the frequency deviations (Δf1, Δf2) and power deviation (ΔP) of the system, respectively, using the GA and DA with the ITAE. Figure 27, Figure 28 and Figure 29 show the frequency deviations (Δf1, Δf2) and power deviation (ΔP) of the system, respectively, using GA with the proposed IB-WGFF.
From Table 8 and Figure 18, Figure 19 and Figure 20, it is found that the dynamic response of Δf2, Δf2, and ΔP, under GA-2DOF-PID-IB-WGF controller is better than the conventional 2DOF-PID and GA-DOF-PID-IB-ITAE. It has smaller overshot, undershot, settling time, and fitness function value.
It is found, from Figure 21, Figure 22 and Figure 23, that the dynamic response of Δf2, Δf2 and ΔP, using the DA-2DOF-PID-ITAE controller is better than the classical 2DOF-PID and GA-DOF-PID-IB-ITAE controllers.
Figure 24, Figure 25 and Figure 26 and Table 8 indicated that the dynamic performance of the system, implementing the proposed DA-2DOF-PID-IB-WGFF, is better than the classical 2DOF-PID and the GA-2DOF-PID-IB-WGFF.
Finally, it is concluded that the dynamic performance of the system, implementing the proposed DA-2DOF-PID-IB-WGFF, is the best in terms of overshot (zero), undershot (0.0244), settling time (1.543), and fitness function value (0.399303).
Figure 27 indicates a graph between the fitness function value versus the generations number using GA and DA in case of the ITAE fitness function. Figure 28 shows a graph using GA and DA in case of the proposed IB-WGFF. It is found, from Figure 27 and Figure 28, that the speed of convergence of DA is higher than GA in case of ITAE and the Proposed IB-WGFF. In addition, the DA has a value of fitness function of (0.399303) which is better than the value (0.657643) obtained in the case of GA.

6. Verification of the Proposed Control Scheme

This section of the article is dedicated for demonstrating the performance and verifying the superiority and robustness of the proposed control scheme under the parameters perturbations and load disturbances. Verification is done under the same parameters of controllers designed in Section 5.5 and indicated in Table 8. The following subsections demonstrates the dynamic response of the two-area power system (Δf2, Δf2, and ΔP) under parameters perturbations and load disturbances.

6.1. Parameters Perturbations

The parameters of the interconnected power system are constantly varying, this may seriously affect the system performance. Robustness against dynamic faults, uncertainties, and environmental changes is one of the essential gains of using the AI control approaches. To verify the stability, efficacy, and robustness of the proposed DA-2DOF-PID control scheme, the parameters (R, B, H, D, Tg, Tt, and T12) of the studied system given in Figure 5, were significantly changed according to Table 9. The dynamic responses of Δf2, Δf2, and ΔP of the system, after implementing these changes under the application of the same conventional 2DOF-PID, DA-2DOF-PID-ITAE, GA-2DOF-PID-ITAE, GA-2DOF-PID-IB-WGFF, and the proposed DA-2-DOF-PID-IB-WGFF are demonstrated in Figure 29, Figure 30 and Figure 31, respectively.
It can be concluded from Figure 29, Figure 30 and Figure 31 that the conventional-2DOF-PID controller cannot handle the implemented system parameters perturbation given in Table 9. It is evident that the DA-2DOF-PID-ITAE, GA-2DOF-PID-ITAE, and GA-2DOF-PID-IB-WGFF controllers can follow and handle the implemented perturbations on the parameters but the dynamic response of the proposed 2DOF-PID control scheme tuned using the IB-WGFF (for C1 = 0.5 and C2 = 0.5) is most robust and the smoothest one. It has the smallest undershot, overshoot, and settling time.

6.2. The Impact of the Parameters Perturbations

This case tested the impact of the parameters perturbations on the studied two-area power system. This case was checked because it can generate high changes in the power system dynamic. In this case, the same load was applied allowing the uncertainties shown in Table 9. The frequency performance of the system against these changes was tested. Figure 32 shows the frequency response of the nominal and perturbed system under the proposed 2DOF-PID optimized using DA via the IB-WGFF. It is concluded, from Figure 32, that the proposed control scheme effectively enhances the dynamic performance. These results prove the significant role of the proposed scheme in supporting the power system frequency that occurs with highly disrupting conditions.

6.3. Effect of Parameter Variation on System Response

Figure 33 indicates the dynamic response of frequency for the first area with 50% increase in T12. The dynamic response of the system under the proposed DA-2DOF-PID is more robust and provides better performance compared to the other controllers.

6.4. Disturbance Load

Load variation results in frequency and power variation which may affect the power quality and system stability. These effects may reduce the damping performance of the controlled power system and could even lead to instability in the power system if these changes exceed the allowable limits. To verify the robustness of the proposed control scheme, a varying load disturbance shown in Figure 34 under the same parameters obtained in Section 5.5 is implemented to the two-area power system shown in Figure 5.
In this scenario, the ability of the proposed 2DOF-PID tuned using the designed IB-WGFF, under the effect of radical multi-step load disturbance, was evaluated and proved. The step load is applied to the studied system with varying step pattern increasing at the time instants t= 5 s, 35 s, and 45 s and decreasing at t = 15 and 25 s. Figure 35, Figure 36 and Figure 37 indicate the dynamic of the two area power system implementing the conventional 2DOF-PID, the DA-2DOF-PID- ITAE, the GA-2DOF-PID-ITAE, GA-2DOF-PID-IB-WGFF, and the DA-2DOF-PID-IB-WGFF (for C1 = 0.5 and C2 = 0.5).
It is evident that the dynamic behaviors of the power and frequency, under the proposed DA-2DOF-PID-IB-WGFF, have minimum deviations in comparison with the dynamics of all other implemented controllers. The proposed DA-2DOF-PID-IB-WGFF also have the fastest control action.
It is obvious from Figure 35 and Figure 36 that the novel proposed DA-2DOF-PID PID control scheme tuned via the IB-WGFF (for C1 = 0.5 and C2 = 0.5) can robustly keep the deviations in frequency within values between −0.028 and +0.012 for area1 and within values between −0.038 and 0.017 for area2 with the smoothest and fastest settling time. Contrariwise, udder the DA-2-DOF-PID controller tuned via the ITAE comes in the second performance with deviations in frequency within values between −0.032 and 0.013 for area1 and within values between −0.044 and 0.024 for area2. The GA-2DOF-PID tuned via the IB-WGFF comes in the third order of performance with deviations in frequency within values between −0.041 and 0.021 for area1 and within values between −0.052 and 0.027 for area2. The GA-2DOF-PID tuned via the ITAE comes in the fourth order of performance with deviations in frequency within values between −0.043 and 0.023 for area1 and within values between −0.058 and 0.031 for area2. The conventional-2DOF-PID comes in the fifth order of performance with deviations in frequency within values between −0.056 and 0.028 for area1 and within values between −0.068 and 0.038 for area2.
It is evident from Figure 37 that that the new proposed DA-2DOF-PID PID control scheme tuned via the IB-WGFF (for C1 = 0.5 and C2 = 0.5) can robustly keep the power deviation (ΔP) within values of about −0.005 and 0.008 pu. The other controllers can restore the tie-line power dynamics to its steady state of zero value but with high undershoots and overshoots and with longer settling times.
Again, it is concluded From Figure 35, Figure 36 and Figure 37 that when a varying load is applied to the two-area power system, the DA-2DOF-PID controller tuned via the proposed IB-WGFF (for C1 = 0.5 and C2 = 0.5) is the best controller in terms of the overshot, undershot, and settling time. Therefore, it evident, from the above discussion, that the dynamic performance of the system is improved implementing the proposed DA-2-DOF-PID controller tuned via the designed IB-WGFF.

7. Conclusions

In this article, an improved 2DOF-PID controller is designed for frequency/power control of a two-area interconnected power system. The implemented 2DOF-PID controller is optimized using Dragonfly Algorithm via the proposed IB-WGFF. For comparison, the result of the proposed controller is compared with the results of the classical-2DOF-PID controller and the 2DOF-PID tuned using Dragonfly Algorithm and Genetic Algorithm (GA) via the most common performance criterion used in the literature (i.e., ISE, ITAE, ITSE, and IAE). The simulation results proved that the proposed control scheme is better than the classical 2DOF-PID and the 2DOF-PID controllers tuned using DA and GA via the ITAE, ISE, IAE, and the ITSE. The proposed 2-DOF-PID controller is the smoothest and the best in overshot, undershot, and settling time, and steady-state error. For proving the robustness and adaptive property of the proposed controller, the frequency deviation and the tie-line power of the two-area system are verified under the implementation of the same controllers for load disturbances and parameters perturbations. It is concluded that the 2DOF-PID control scheme tuned using the Dragonfly Algorithm via the proposed IB-WGFF (for C1 = 0.5 and C2 = 0.5) is the best controller in suppressing the frequency oscillations and tie-line power fluctuations under the load disturbance and system uncertainties.

Author Contributions

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

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Essiet, I.O.; Sun, Y.; Wang, Z. Optimized energy consumption model for smart home using improved differential evolution algorithm. Energy 2019, 172, 354–3651. [Google Scholar] [CrossRef]
  2. Tungadio, D.H.; Bansal, R.C.; Siti, M.W. Energy flow estimation-control of two interconnected microgrid. J. Energy S. Afr. 2008, 29, 69–80. [Google Scholar] [CrossRef]
  3. Barakat, M. Novel chaos game optimization tuned-fractional-order PID fractional-order PI controller for load–frequency control of interconnected power systems. Prot. Control Mod. Power Syst. 2022, 7, 16. [Google Scholar] [CrossRef]
  4. Kumar, N.; Malik, H.; Singh, A.; Alotaibi, M.A.; Nassar, M.E. Novel Neural Network-Based Load Frequency Control Scheme: A Case Study of Restructured Power System. IEEE Access 2021, 9, 162231–162242. [Google Scholar] [CrossRef]
  5. Saxena, S. Load frequency control strategy via fractional order controller and reduced-order modeling. Int. J. Electr. Power Energy Syst. 2019, 104, 603–6014. [Google Scholar] [CrossRef]
  6. Barakat, M.; Donkol, A.; Hamed, H.F.A.; Salama, G.M. Harris Hawks-Based optimization algorithm for automatic LFC of the interconnected power system using PD-PI cascade control. J. Electr. Eng. Technol. 2021, 16, 1845–1865. [Google Scholar] [CrossRef]
  7. Lu, K.-D.; Zeng, G.-Q.; Zhou, W. Adaptive constrained population external optimisation-based robust proportional-integral-derivation frequency control method for an islanded microgrid. IET Cyber Syst. Robot. 2021, 3, 210–227. [Google Scholar] [CrossRef]
  8. Hossain, S.A.; Roy, S.; Karmaker, A.; Islam, R. Performance improvement of PID controller for AVR system using Particle Swarm Optimization. In Proceedings of the 2015 International Conference on Advances in Electrical Engineering (ICAEE), Dhaka, Bangladesh, 17–19 December 2015. [Google Scholar]
  9. Abdelaziz, A.Y.; Ali, E.S. Load Frequency Controller Design via Artificial Cuckoo Search. Electr. Power Compon. Syst. 2016, 44, 90–98. [Google Scholar] [CrossRef]
  10. Fathy, A.; Kassem, A.M.; Abdelaziz, A.Y. Optimal Design of Fuzzy-PID Controller for Deregulated LFC of Multi-Area Power System via Mine Blast Algorithm. Neural Comput. Appl. 2020, 32, 4531–4551. [Google Scholar] [CrossRef]
  11. Ebrahim, M.A.; Becherif, M.; Abdelaziz, A.Y. PID/FOPID-based Frequency control of zero-carbon multi-sources based interconnected power systems under deregulated scenarios. Int. Trans. Electr. Energy Syst. 2021, 31, e12712. [Google Scholar] [CrossRef]
  12. Praveena, P.; Shubhanka, H.S.; Neha, V.; Prasath, P.; Veda, G.S. Load Frequency Control of Single Area System Using Ziegler-Nichols and Genetic Algorithm. Asian J. Electr. Sci. 2019, 8, 62–64. [Google Scholar] [CrossRef]
  13. Kouba, N.E.Y.; Menaa, M.; Hasni, M.; Boudour, M. Optimal control of frequency and voltage variations using PID controller based on Particle Swarm Optimization. In Proceedings of the 2015 4th International Conference on Systems and Control (ICSC), Sousse, Tunisia, 27–30 April 2015; pp. 424–429. [Google Scholar]
  14. Arya, Y. AGC of two-area electric power systems using optimized fuzzy PID with filter plus double integral controller. J. Frankl. Inst. 2018, 355, 4583–4617. [Google Scholar] [CrossRef]
  15. Gupta, D.; Goyal, G.; Kumar, J. An Optimized Fractional Order PID Controller for Integrated Power System. In ICICCT 2019—System Reliability, Quality Control, Safety, Maintenance and Management, Proceedings of the ICICCT 2019 International Conference on Intelligent Computing and Communication Technologies, Hyderabad, India, 9–11 January 1019; Gunjan, V.K., Diaz, V.G., Cardona, M., Solanki, V.K., Sunitha, K.V.N., Eds.; Springer: Singapore, 2019; pp. 663–672. [Google Scholar]
  16. Lamba, R.; Singla, S.K.; Sondhi, S. Design of Fractional Order PID Controller for Load Frequency Control in Perturbed Two Area Interconnected System. Electr. Power Compon. Syst. 2019, 47, 998–1011. [Google Scholar] [CrossRef]
  17. Singh, R.; Kesarwani, S.K.; Gupta, N.K.; Ashfaq, H. Design of 2-DOF PID Controller for Load Frequency Control of Two Area Power System using MFO Algorithm. Int. J. Eng. Adv. Technol. 2020, 9, 158–161. [Google Scholar] [CrossRef]
  18. Karanam, A.N.; Shaw, B. A new two-degree of freedom combined PID controller for automatic generation control of a wind integrated interconnected power system. Prot. Control Mod. Power Syst. 2022, 7, 20. [Google Scholar] [CrossRef]
  19. Kouba, N.E.Y.; Menaa, M.; Hasni, M.; Boudour, M. A Novel Optimal Combined Fuzzy PID Controller Employing Dragonfly Algorithm for Solving Automatic Generation Control Problem. Electr. Power Compon. Syst. 2018, 46, 2054–2070. [Google Scholar] [CrossRef]
  20. Guha, D.; Roy, P.K.; Banerjee, S. Quasi-oppositional JAYA optimized 2-degree-of-freedom PID controller for load-frequency control of interconnected power systems. Int. J. Model. Simul. 2022, 42, 63–85. [Google Scholar] [CrossRef]
  21. Hassan, A.; Aly, M.; Elmelegi, A.; Nasrat, L.; Watanabe, M.; Mohamed, E.A. Optimal Frequency Control of Multi-Area Hybrid Power System Using New Cascaded TID-PIλDμN Controller Incorporating Electric Vehicles. Fractal Fract. 2022, 6, 548. [Google Scholar] [CrossRef]
  22. Gheisarnejad, M. An effective hybrid harmony search and cuckoo optimization algorithm based fuzzy PID controller for load frequency control. Appl. Soft Comput. 2018, 65, 121–138. [Google Scholar] [CrossRef]
  23. Barakat, M. Optimal design of fuzzy-PID controller for automatic generation control of multi-source interconnected power system. Neural Comput. Appl. 2022, 34, 18859–18880. [Google Scholar] [CrossRef]
  24. Alhalabi, M.; Rashed, A.; Barham, Y.; Salim, R.; Ghazal, M. Two-Area Load Frequency Control for Power System Dynamic Performance Enhancement with a Graphical User Interface Integration. In Proceedings of the 2022 International Conference on Electrical, Computer and Energy Technologies (ICECET), Prague, Czech Republic, 20–22 July 2022. [Google Scholar]
  25. Sahu, R.K.; Panda, S.; Padhanm, S. A novel hybrid gravitational search and pattern search algorithm for load frequency control of nonlinear power system. Appl. Soft Comput. 2015, 29, 310–327. [Google Scholar] [CrossRef]
  26. Gbadega, P.A.; Saha, A.K. Load Frequency Control of a Two-Area Power System with a Stand-Alone Microgrid Based on Adaptive Model Predictive Control. IEEE J. Emerg. Sel. Top. Power Electron. 2021, 9, 7253–7263. [Google Scholar] [CrossRef]
  27. Alaa Abdel-hamed, M.; Abou El-Eyoun, K.M.; Amged El-Wakeel, S. Optimized Control Scheme for Frequency/Power Regulation of Microgrid for Fault Tolerant operation. Electr. Power Compon. Syst. 2016, 44, 1429–1440. [Google Scholar] [CrossRef]
  28. Tripathy, D.; Choudhury, N.D.; Sahu, B. Comparative performance assessment of several Fractional Order -Two Degree Freedom controllers tuned using GOA for LFC. Procedia Comput. Sci. 2020, 167, 2022–2032. [Google Scholar] [CrossRef]
  29. Azzini, I.; Muresano, R.; Ratto, M. Dragonfly: A multi-platform parallel toolbox for MATLAB/Octave computer Languages. Syst. Struct. 2018, 52, 21–42. [Google Scholar]
  30. Mirjalili, S. Dragonfly algorithm: A new meta-heuristic optimization technique for solving single-objective, discrete, and multiobjective problems. Neural Comput. Appl. 2015, 27, 1053–1073. [Google Scholar] [CrossRef]
  31. Venkatesh, M.; Sudheer, G. Optimal Load Frequency Regulation of Micro-Grid Using Dragonfly Algorithm. Int. Res. J. Eng. Technol. 2017, 4, 978–981. [Google Scholar]
  32. 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, 4, 2563–2581. [Google Scholar] [CrossRef]
  33. Fathy, A.; Yousri, D.; Rezk, H.; Thanikanti, S.B.; Hasanien, H.M. A Robust Fractional-Order PID Controller Based Load Frequency Control Using Modified Hunger Games Search Optimizer. Energies 2022, 15, 361. [Google Scholar] [CrossRef]
  34. Rabindra, K.S.; Tulasichandra, S.G.; Sidhartha, P. Automatic generation control of multi-area power systems with diverse energy sources using Teaching Learning Based Optimization algorithm. Eng. Sci. Technol. Int. J. 2016, 19, 113–134. [Google Scholar]
  35. Ali, T.; Malik, S.A.; Hameed, I.A.; Daraz, A.; Mujlid, H.; Azar, A.T. Load Frequency Control and Automatic Voltage Regulation in a Multi-Area Interconnected Power System Using Nature-Inspired Computation-Based Control Methodology. Sustainability 2022, 14, 12162. [Google Scholar] [CrossRef]
  36. Elwakel, A.S.; Ellissy, A.; Abdelhamed, A.M. A Hybrid Bacterial Foraging-Particle Swarm Optimization Technique for Optimal Tuning of Proportional-Integral-Derivative Controller of a Permanent Magnet Brushless DC Motor. Electr. Power Compon. Syst. 2015, 3, 309–319. [Google Scholar] [CrossRef]
  37. Cam, E.; Gorel, G.; Mamur, H. Use of the Genetic Algorithm-Based Fuzzy Logic Controller for Load-Frequency Control in a Two Area Interconnected Power System. Appl. Sci. 2017, 7, 308. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Circuit diagrams for the studied system (a) Circuit topology (b) LFC implemented a single area (c) Complete block diagram.
Figure 1. Circuit diagrams for the studied system (a) Circuit topology (b) LFC implemented a single area (c) Complete block diagram.
Energies 16 00486 g001aEnergies 16 00486 g001b
Figure 2. 2DOF-PID control architecture.
Figure 2. 2DOF-PID control architecture.
Energies 16 00486 g002
Figure 3. Tuning process of the parameters using IB-WGFF.
Figure 3. Tuning process of the parameters using IB-WGFF.
Energies 16 00486 g003
Figure 4. Dragonfly Algorithm (DA) Implementation.
Figure 4. Dragonfly Algorithm (DA) Implementation.
Energies 16 00486 g004
Figure 5. Complete SIMULINK model for the two-area power system.
Figure 5. Complete SIMULINK model for the two-area power system.
Energies 16 00486 g005
Figure 6. Frequency deviation ( Δ f 1 ) using conventional 2DOF-OPID controller.
Figure 6. Frequency deviation ( Δ f 1 ) using conventional 2DOF-OPID controller.
Energies 16 00486 g006
Figure 7. Frequency deviation ( Δ f 2 ) using conventional 2DOF-OPID controller.
Figure 7. Frequency deviation ( Δ f 2 ) using conventional 2DOF-OPID controller.
Energies 16 00486 g007
Figure 8. Tie-line Power deviation (ΔP) using conventional 2DOF-PID controller.
Figure 8. Tie-line Power deviation (ΔP) using conventional 2DOF-PID controller.
Energies 16 00486 g008
Figure 9. Frequency deviation Δ f 1 for area1 using the proposed IB-WGFF.
Figure 9. Frequency deviation Δ f 1 for area1 using the proposed IB-WGFF.
Energies 16 00486 g009
Figure 10. Frequency deviation Δ f 2 for area 2 using the proposed IB-WGFF.
Figure 10. Frequency deviation Δ f 2 for area 2 using the proposed IB-WGFF.
Energies 16 00486 g010
Figure 11. Power deviation Δ P using the proposed IB-WGFF.
Figure 11. Power deviation Δ P using the proposed IB-WGFF.
Energies 16 00486 g011
Figure 12. Frequency deviation Δ f 1 for area 1 using performance objectives in the literature.
Figure 12. Frequency deviation Δ f 1 for area 1 using performance objectives in the literature.
Energies 16 00486 g012
Figure 13. Frequency deviation Δ f 2 implementing performance objectives used in the literature.
Figure 13. Frequency deviation Δ f 2 implementing performance objectives used in the literature.
Energies 16 00486 g013
Figure 14. Power deviation Δ P using performance indices used in the literature.
Figure 14. Power deviation Δ P using performance indices used in the literature.
Energies 16 00486 g014
Figure 15. Δ f 1 using the proposed 2DOF-PID controller compared with other controllers.
Figure 15. Δ f 1 using the proposed 2DOF-PID controller compared with other controllers.
Energies 16 00486 g015
Figure 16. Δ f 2 using the proposed 2DOF-PID controller compared with other controllers.
Figure 16. Δ f 2 using the proposed 2DOF-PID controller compared with other controllers.
Energies 16 00486 g016
Figure 17. ΔP using the proposed 2DOF-PID controller compared with other controllers.
Figure 17. ΔP using the proposed 2DOF-PID controller compared with other controllers.
Energies 16 00486 g017
Figure 18. Δ f 1 , using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Figure 18. Δ f 1 , using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Energies 16 00486 g018
Figure 19. Δ f 2 , using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Figure 19. Δ f 2 , using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Energies 16 00486 g019
Figure 20. ΔP using GA using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Figure 20. ΔP using GA using GA implementing ITAE, IB-WGFF and classical 2DOF-PID.
Energies 16 00486 g020
Figure 21. Δ f 1 , using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Figure 21. Δ f 1 , using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Energies 16 00486 g021
Figure 22. Δ f 2 , using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Figure 22. Δ f 2 , using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Energies 16 00486 g022
Figure 23. ΔP using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Figure 23. ΔP using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with ITAE.
Energies 16 00486 g023
Figure 24. Δ f 1 using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Figure 24. Δ f 1 using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Energies 16 00486 g024
Figure 25. Δ f 2 using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Figure 25. Δ f 2 using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Energies 16 00486 g025
Figure 26. ΔP using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Figure 26. ΔP using GA-2DOF-PID, DA-2DOF-PID, and classical-2DOF-PID with IB-WGFF.
Energies 16 00486 g026
Figure 27. Comparative graph for the fitness function value using ITAE.
Figure 27. Comparative graph for the fitness function value using ITAE.
Energies 16 00486 g027
Figure 28. Comparative graph for the fitness function value using proposed IB-WGFF.
Figure 28. Comparative graph for the fitness function value using proposed IB-WGFF.
Energies 16 00486 g028
Figure 29. Frequency deviation ( Δ f 1 ) of area 1 under parameter perturbations.
Figure 29. Frequency deviation ( Δ f 1 ) of area 1 under parameter perturbations.
Energies 16 00486 g029
Figure 30. Frequency deviation ( Δ f 2 ) of area 2 under parameter perturbations.
Figure 30. Frequency deviation ( Δ f 2 ) of area 2 under parameter perturbations.
Energies 16 00486 g030
Figure 31. Tie-line power deviation (ΔP) after implementing parameter perturbations.
Figure 31. Tie-line power deviation (ΔP) after implementing parameter perturbations.
Energies 16 00486 g031
Figure 32. Dynamic response under normal and perturbed system using the proposed scheme.
Figure 32. Dynamic response under normal and perturbed system using the proposed scheme.
Energies 16 00486 g032
Figure 33. Dynamic Response due to the increase in T12 by a value of 50%.
Figure 33. Dynamic Response due to the increase in T12 by a value of 50%.
Energies 16 00486 g033
Figure 34. Implemented varying load disturbance (PL).
Figure 34. Implemented varying load disturbance (PL).
Energies 16 00486 g034
Figure 35. Δ f 1 for area 1 under load disturbances.
Figure 35. Δ f 1 for area 1 under load disturbances.
Energies 16 00486 g035
Figure 36. Δ f 2 for Area 2 under load disturbances.
Figure 36. Δ f 2 for Area 2 under load disturbances.
Energies 16 00486 g036
Figure 37. ΔP under load disturbances using different controllers.
Figure 37. ΔP under load disturbances using different controllers.
Energies 16 00486 g037
Table 1. Studied power system parameters.
Table 1. Studied power system parameters.
ParametersValuesParametersValuesParametersValues
R12.4D20.90B120.1
D1−0.60H24.0T120.545
H15.0Tg20.30R20.0625
Tg10.20Tt20.60Kg1.0
Tt10.50B216.90Kt1.0
Table 2. DA’s parameters.
Table 2. DA’s parameters.
ParametersValues
Population Size50
Number of iterations30
A, Separation weighted coefficients,0.25
B, Alignment weighted coefficients0.15
C, Cohesion weighted coefficients0.96
D, Food weighted coefficients0.6
E, is the enemy factor1
ω , inertial weighted coefficient[0.25, 95]
Table 3. GA Parameters.
Table 3. GA Parameters.
ParametersValues
population size50
Number of generations30
Acceleration constant, C10.5
Acceleration constant, C21.5
Initial inertia weight, wmax0.9
Final inertia weight, wmin0.4
Table 4. Efforts implementing conventional-FOPID controller.
Table 4. Efforts implementing conventional-FOPID controller.
FO-PID
Parameters
K P K I K D P W DWSeeking Time
Area 10.92130.70240.83120.90340.4434Very large time
Area 20.83240.70340.53540.60150.6415
Table 5. System efforts using DA-FOPID controller via the proposed IB-WGFF.
Table 5. System efforts using DA-FOPID controller via the proposed IB-WGFF.
Parameters2-DOF PID Controller (AREA 1)2-DOF PID Controller (AREA 2)Fitness Function
C1 and C2KpKiKdPWDWKpKiKdPWDWValue
0 and 11.00000.86240.63171.73201.18000.75051.00000.68401.66481.73150.106046
0.1 and 0.91.00000.85280.41631.55771.05540.69070.96100.55881.40081.29630.250180
0.2 and 0.81.00000.82910.26982.00002.00001.00001.00000.63491.74902.00000.262111
0.3 and 0.71.00000.86140.81692.00001.35340.91151.00000.831411.707592.00000.278407
0.4 and 0.61.00000.87120.41381.91431.51600.78191.00000.97551.82521.00900.352105
0.5 and 0.51.00000.86800.95551.80340.83081.00001.00000.61341.29112.00000.399303
0.6 and 0.40.99100.87590.42851.79241.49680.83161.00000.87261.49251.09820.465029
0.7 and 0.30.77330.87650.29562.00001.6655 0.6831 1.0000 1.0000 1.43740.81610.517983
0.8 and 0.21.00001.00000.805872.0000 1.11010.886871.00000.55451.1164 1.42860.625108
0.9 and 0.10.84480.89200.78571.88410.69180.84341.00000.66761.12991.17320.609655
1 and 01.00001.0000.6936 2.00001.17170.87811.00000.53421.19301.49380.641395
Table 6. Efforts of different performance criteria in the literature.
Table 6. Efforts of different performance criteria in the literature.
Response and PID
Parameters
ITAEITSEIAEISE
Area 1 K P 0.966421.000001.000001.00000
K I 0.937010.967731.000000.94818
K D 0.797960.828920.986771.00000
PW2.000001.940392.000002.00000
DW0.956021.864941.689532.00000
Area 2 K P 0.706930.878030.732501.00000
K I 0.60001.000001.000000.95282
K D 0.786721.000001.000000.81901
PW1.643581.153381.492481.06404
DW1.095071.339941.349482.00000
FF Val.0.352860.007510.046150.00100
Seeking Time345.105347.8026325.40310.12
Table 7. Efforts using conventional 2-DOF-PID, DA-2-DOF-PID-ITAE, and DA-2-DOF-PID-IB-WGFF.
Table 7. Efforts using conventional 2-DOF-PID, DA-2-DOF-PID-ITAE, and DA-2-DOF-PID-IB-WGFF.
Response and PID
Parameters
ITAEConventionalProposed IB-WGFF
Area 1Area 2Area 1Area 2Area 1Area 2
K P 0.966420.706930.921340.822451.00001.0000
K I 0.937011.000000.700240.733450.86801.0000
K D 0.797960.786720.831230.505430.95550.6134
PW2.000001.643580.90340.601561.80341.2911
DW0.956021.095070.443450.641560.83082.0000
FF. Value0.35286 0.399303
Seeking Time345.105Very large time343.261
OvershootΔf10.00000.02350.0
Δf20.00000.0200.0
ΔP0.03280.0920.0153
UndershootΔf10.03550.05450.0241
Δf20.04530.06650.0321
ΔP0.00000.0740.0000
S.S. ErrorΔf16.57 × 10−6344 × 10−61.931 × 10−6
Δf229.99 × 10−6131 × 10−63.143 × 10−6
ΔP18.34 × 10−61059 × 10−64.5682 × 10−6
Settling Time (Sec.)Δf17.4211.526.571
Δf26.839.436.731
ΔP3.72412.451.543
Table 8. Efforts using conventional 2DOF-PID, DA-2DOF-PID-ITAE, and GA-2DOF-PID-ITAE, and DA-2DOF-PID-IB-WGFF.
Table 8. Efforts using conventional 2DOF-PID, DA-2DOF-PID-ITAE, and GA-2DOF-PID-ITAE, and DA-2DOF-PID-IB-WGFF.
Response and PID
Parameters
TechniqueITAEProposed IB-WGOConventional
Area 1Area 2Area 1Area 2Area 1Area 2
K P GA0.726990.278710.839460.7622060.921340.82245
DA0.966420.706931.00001.0000
K I GA0.830200.947350.807000.928900.700240.73345
DA0.937011.000000.86801.0000
K D GA0.600580.393740.723540.806950.831230.50543
DA0.797960.786720.95550.6134
PWGA1.452311.1680941.460780.710850.90340.60156
DA2.000001.643581.80341.2911
DWGA0.903141.630510.727960.880710.443450.64156
DA0.956021.095070.83082.0000
FF. ValueGA0.7387610.657643
DA0.3528650.399303
Seeking TimeGA445.234442.343Very large time
DA345.105343.261
OvershootΔf1GA0.01530.00830.0235
DA0.00000.000
Δf2GA0.01780.00730.020
DA0.00000.000
ΔPGA0.04350.02730.092
DA0.03280.0153
UndershootΔf1GA0.04930.04120.0545
DA0.03550.0244
Δf2GA0.05450.0480.0665
DA0.04530.0321
ΔPGA0.03210.02340.074
DA0.00000.000
S.S. ErrorΔf1GA45.372 × 10−69.16 × 10−6344 × 10−6
DA6.57 × 10−61.931 × 10−6
Δf2GA86.31 × 10−632.78 × 10−6131 × 10−6
DA29.99 × 10−63.143 × 10−6
ΔPGA35.34 × 10−615.856 × 10−61059 × 10−6
DA18.34 × 10−64.5682 × 10−6
Settling Time (Sec.)Δf1GA9.3436.31211.52
DA7.4236.571
Δf2GA9.8436.2339.43
DA6.8326.731
ΔPGA7.5767.32212.45
DA3.7241.543
Table 9. Power system parameter variations.
Table 9. Power system parameter variations.
ParameterVariation Range (%)ParameterVariation Range (%)
R1 = 2.4;+40% (3.36)D2 = 0.90;+10% (0.99)
D1 = 0.60;−40% (0.36)H2 = 4.0;+10% (4.4)
H1 = 5.0;+50% (7.5)Tg2 = 0.30;−10% (0.33)
Tg1 = 0.20;+10% (0.22)Tt2 = 0.60;+10% (0.66)
Tt1 = 0.50;+25% (0.625)B2 = 16.90;+20% (20.28)
B1 =20.1;10% (22.11)T12=0.545;+10% (0.5995)
R2 = 0.0625;−10% (0.05625)
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

Abdel-hamed, A.M.; Abdelaziz, A.Y.; El-Shahat, A. Design of a 2DOF-PID Control Scheme for Frequency/Power Regulation in a Two-Area Power System Using Dragonfly Algorithm with Integral-Based Weighted Goal Objective. Energies 2023, 16, 486. https://doi.org/10.3390/en16010486

AMA Style

Abdel-hamed AM, Abdelaziz AY, El-Shahat A. Design of a 2DOF-PID Control Scheme for Frequency/Power Regulation in a Two-Area Power System Using Dragonfly Algorithm with Integral-Based Weighted Goal Objective. Energies. 2023; 16(1):486. https://doi.org/10.3390/en16010486

Chicago/Turabian Style

Abdel-hamed, Alaa M., Almoataz Y. Abdelaziz, and Adel El-Shahat. 2023. "Design of a 2DOF-PID Control Scheme for Frequency/Power Regulation in a Two-Area Power System Using Dragonfly Algorithm with Integral-Based Weighted Goal Objective" Energies 16, no. 1: 486. https://doi.org/10.3390/en16010486

APA Style

Abdel-hamed, A. M., Abdelaziz, A. Y., & El-Shahat, A. (2023). Design of a 2DOF-PID Control Scheme for Frequency/Power Regulation in a Two-Area Power System Using Dragonfly Algorithm with Integral-Based Weighted Goal Objective. Energies, 16(1), 486. https://doi.org/10.3390/en16010486

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