Use of the Genetic Algorithm-Based Fuzzy Logic Controller for Load-Frequency Control in a Two Area Interconnected Power System

The use of renewable energy resources has created some problems for power systems. One of the most important of these is load frequency control (LFC). In this study, in order to solve the LFC problem, modern control methods were applied to a two area multi source interconnected power system. A photovoltaic solar power plant (PV-SPP) was also connected, in order to identify the harmful effects on the frequency of the system. A new Genetic-based Fuzzy Logic (GA-FL) controller was designed to control the frequency of the system. For comparison, conventional proportional-integral-derivative (PID), fuzzy logic (FL), and Genetic Algorithm (GA)-PID controllers were also designed. The new control method exhibited a better performance than the conventional and other modern control methods, because of the low overshoot and short settling time. All simulations were realized with the Matlab-Simulink program.


Introduction
Electricity generated from renewable energy resources (RESs) is used as grid-on or grid-off power.Sudden climatic changes can cause large changes in their output parameters.This situation is also seen in PV-SPP systems, because of shading.Their output voltages and frequencies are affected by this [1].In the case of large frequency or voltage changes, the network is in danger of collapse.Moreover, the quality of the energy also deteriorates.These situations can be harmful to both the grid and the consumers [2].One way to eradicate this kind of problem is through load frequency control (LFC).LFC is one of the most important processes during the operation of a power system.The energy balance between the generation and loads is controlled by means of a numerical value obtained from sensors that measure the grid frequency [3,4].For this reason, to solve the load frequency problem, several control methods have been proposed by previous researchers.The aim of these studies has been to maintain the frequency at a desired value.Another problem is that multi-area power systems have a large number of parameters that must be controlled and have enormous, complex parts.Therefore, the control procedures and operation of the systems are more complex.In particular, it is quite difficult to stabilize the frequency in each area, because of the increase of the multi area.Consequently, more intelligent controllers have to be designed for LFC [5,6].The most common control strategy for such systems is the PID controller.Many system instabilities can be managed by the conventional PID controller.However, the conventional PID controller is inadequate, owing to the expansion of grids and the rapidly increasing energy demands [7].The utilization of different modern optimization methods embedded in the conventional PID controller is recommended for the load frequency control, by some researchers [8].The cuckoo search algorithm (CSA) and intelligent bee colony optimization (IBCO) are just two of the many examples of these methods [6,9,10].Khalghani et al. compared an emotional controller with the conventional PID controller [11].A sugeno fuzzy logic controller (S-FLC), optimized by particle swarm optimisation (PSO), was applied by Sangawong et al. [12].For micro grid networks with photovoltaic (PV) systems, battery units, micro wind turbines (MWT), and fuel cells, the most suitable control method is a self-adjusting control method [11,13,14].Hence, Shankar et al. realized a genetic algorithm (GA) for LFC in a two area interconnected hydro-thermal power system [15].Additionally, another comparative study presented a controller employing linear matrix inequalities (LMI) and a controller optimized GA (GALMI) for a three area interconnected power system, simulated by Rerkpreedapong et al. [16].They achieved a robust LFC performance.In addition to these, there are also studies that adjust the membership functions of the fuzzy logic controller (FLC), using a GA [17].Daneshfar and Bevrani actualized a FLC process with a GA that was applied to a grid network which had two or three areas, in a classic or hybrid power system [18].For a three area interconnected power system, a controller using an artificial neural network (ANN) was designed to regulate the load frequency imbalance, by Demiroren et al. [19].Moreover, optimization methods embedding flexible alternative current transmission systems (FACTS) were executed on multi area interconnected power systems, in some studies.In the systems, the effects of the FACTSs on the frequency of the line were monitored [20,21].
As the world rapidly changes, the development of electric grids increases the importance of choosing the best control method.Since the frequency value continuously changes, depending on the energy consumption, it is more difficult to maintain the desired frequency value [22,23].In this paper, in order to solve the LFC problem, a new GA-based FL (GA-FLC) controller is suggested, simulated, and verified, for a two area multi source interconnected power system.For this purpose, the system was first installed on Matlab-Simulink and then controlled with the help of various controllers.For comparison, conventional PID, FL, and GA-PID controllers have also been implemented in the system.The values of the overshoot and settling time of the system have been observed and examined.The results show that the proposed controller has better overshoot and settling time values than the others.

The Proposed Power System
In this study, a two area three source interconnected renewable energy power system was modelled, which is shown in Figure 1.A thermal unit (TU), a thermal generator unit with a reheat unit (RTU), and a PV-SPP were connected to the system as energy resources.The TUs and the RTUs were practically employed in some literature for the LFC [24].Since the use of PV-SPPs in grids has gradually increased in the last few years, a PV-SPP was incorporated into the system.The utilization of different modern optimization methods embedded in the conventional PID controller is recommended for the load frequency control, by some researchers [8].The cuckoo search algorithm (CSA) and intelligent bee colony optimization (IBCO) are just two of the many examples of these methods [6,9,10].Khalghani et al. compared an emotional controller with the conventional PID controller [11].A sugeno fuzzy logic controller (S-FLC), optimized by particle swarm optimisation (PSO), was applied by Sangawong et al. [12].For micro grid networks with photovoltaic (PV) systems, battery units, micro wind turbines (MWT), and fuel cells, the most suitable control method is a self-adjusting control method [11,13,14].Hence, Shankar et al. realized a genetic algorithm (GA) for LFC in a two area interconnected hydro-thermal power system [15].Additionally, another comparative study presented a controller employing linear matrix inequalities (LMI) and a controller optimized GA (GALMI) for a three area interconnected power system, simulated by Rerkpreedapong et al. [16].They achieved a robust LFC performance.In addition to these, there are also studies that adjust the membership functions of the fuzzy logic controller (FLC), using a GA [17].
Daneshfar and Bevrani actualized a FLC process with a GA that was applied to a grid network which had two or three areas, in a classic or hybrid power system [18].For a three area interconnected power system, a controller using an artificial neural network (ANN) was designed to regulate the load frequency imbalance, by Demiroren et al. [19].Moreover, optimization methods embedding flexible alternative current transmission systems (FACTS) were executed on multi area interconnected power systems, in some studies.In the systems, the effects of the FACTSs on the frequency of the line were monitored [20,21].
As the world rapidly changes, the development of electric grids increases the importance of choosing the best control method.Since the frequency value continuously changes, depending on the energy consumption, it is more difficult to maintain the desired frequency value [22,23].In this paper, in order to solve the LFC problem, a new GA-based FL (GA-FLC) controller is suggested, simulated, and verified, for a two area multi source interconnected power system.For this purpose, the system was first installed on Matlab-Simulink and then controlled with the help of various controllers.For comparison, conventional PID, FL, and GA-PID controllers have also been implemented in the system.The values of the overshoot and settling time of the system have been observed and examined.The results show that the proposed controller has better overshoot and settling time values than the others.

The Proposed Power System
In this study, a two area three source interconnected renewable energy power system was modelled, which is shown in Figure 1.A thermal unit (TU), a thermal generator unit with a reheat unit (RTU), and a PV-SPP were connected to the system as energy resources.The TUs and the RTUs were practically employed in some literature for the LFC [24].Since the use of PV-SPPs in grids has gradually increased in the last few years, a PV-SPP was incorporated into the system.

PV Power System
In order to determine the electrical properties of the PV-SPP, that provides the power to the interconnected system, it is necessary to observe how the current and voltage of the PV cell are affected by a variable load.When a PV cell is modelled, a sufficient number of PV cell groups (N S ) are connected in series, in order to generate the required output voltage.Also, a sufficient number of PV cell groups (N P ) are connected in parallel for generating the required output current.The PV cell variables of the equivalent circuit depend on the light intensity and temperature.Thus, the light and temperature values should be known, in order to calculate the output power [25].The PV cell is examined, by means of Figure 2. The voltage of the PV cell can be obtained, depending on its current of it by using Equation (1)."I 0 " is defined as the reverse saturation current of the diode.Definitions of the other symbols are shown in Abbreviations and Nomenclatures.

PV Power System
In order to determine the electrical properties of the PV-SPP, that provides the power to the interconnected system, it is necessary to observe how the current and voltage of the PV cell are affected by a variable load.When a PV cell is modelled, a sufficient number of PV cell groups ( ) are connected in series, in order to generate the required output voltage.Also, a sufficient number of PV cell groups ( ) are connected in parallel for generating the required output current.The PV cell variables of the equivalent circuit depend on the light intensity and temperature.Thus, the light and temperature values should be known, in order to calculate the output power [25].The PV cell is examined, by means of Figure 2. The voltage of the PV cell can be obtained, depending on its current of it by using Equation (1)." " is defined as the reverse saturation current of the diode.Definitions of the other symbols are shown in Abbreviations and Nomenclatures.

1
(2) where Ta and Tx represent the known reference ambient temperature and the different ambient temperature at different points during the tests, respectively.For the PV-SPP model in this study, the Tx value is taken as 28 °C.The βT and γT coefficients represent the slopes of changes in the cell voltage and current due to temperature, respectively.The βT and γT coefficients vary, depending on the PV cell type, and are determined experimentally.βT is calculated between 0.004 and 0.006, and γT is calculated between 0.02 and 0.1 [25].
Although the ambient temperature does not significantly change throughout the day, the solar radiation effectively changes because of the amount of sunlight and the cloudiness of the sky.As a result of this, the photo current, operating temperature, and output voltage of the PV cell are affected.These effects are given for the output voltage by CSV and for the photo-current by CSI, and are defined by Equations ( 4) and (5), respectively.

1
(4) where Scell is the daylight intensity used as the reference solar radiation level.Sx represents different solar radiation levels at different times, depending on the atmospheric changes.αs is the coefficient of variation in the solar radiation level, and is defined by Equation ( 6).  2) and (3), respectively.
where T a and T x represent the known reference ambient temperature and the different ambient temperature at different points during the tests, respectively.For the PV-SPP model in this study, the T x value is taken as 28 • C. The β T and γ T coefficients represent the slopes of changes in the cell voltage and current due to temperature, respectively.The β T and γ T coefficients vary, depending on the PV cell type, and are determined experimentally.β T is calculated between 0.004 and 0.006, and γ T is calculated between 0.02 and 0.1 [25].
Although the ambient temperature does not significantly change throughout the day, the solar radiation effectively changes because of the amount of sunlight and the cloudiness of the sky.As a result of this, the photo current, operating temperature, and output voltage of the PV cell are affected.These effects are given for the output voltage by C SV and for the photo-current by C SI , and are defined by Equations ( 4) and (5), respectively.
where S cell is the daylight intensity used as the reference solar radiation level.S x represents different solar radiation levels at different times, depending on the atmospheric changes.α s is the coefficient of variation in the solar radiation level, and is defined by Equation (6).
The value of α s is different for PV cells that do not have the same properties and this value is determined experimentally.However, this value is usually between 0.3 and 0.4 • C cm 2 /mW.By using the correction coefficients C TV , C TI , C SV , and C SI given above, the new output voltage V Xcell of the PV cell and the new photo current I XPV are calculated, as shown in Equations ( 7) and ( 8), for the new working temperature T X and the new solar radiation level S X .In these expressions, V cell and I PV are the values of the cell output voltage and the photo current at the reference cell operating temperature and solar radiation levels.A MATLAB-Simulink model of PV-SPP was explained in detail in Ref. [25].
The PV cell Simulink model used in the paper is shown in Figure 3a,b.
Appl.Sci.2017, 7, 308 4 of 21 The value of αs is different for PV cells that do not have the same properties and this value is determined experimentally.However, this value is usually between 0.3 and 0.4 °C cm 2 /mW.By using the correction coefficients CTV, CTI, CSV, and CSI given above, the new output voltage VXcell of the PV cell and the new photo current IXPV are calculated, as shown in Equations ( 7) and ( 8

Single Area Interconnected Power System
The blocks used in a single area interconnected power system are generally nonlinear systems.The loads at the operating point and the PV-SPP affect the frequency of the system.A LFC applied to a single area interconnected power system, which is shown in Figure 4, has a turbine, a generator, and other system blocks.Additionally, a PV plant generated through sunlight is demonstrated.For the system, a control responds as quickly as possible.The control process has been realized in two phases, as a primer and a seconder frequency control.The primer frequency control is achieved by adjusting the power balance with the turbine speed regulator, when there is a decrease or increase in electrical power between two different areas.The seconder frequency control is achieved by controlling the power flow between the areas, which reduces the permanent frequency error in the interconnected power system to zero.

Single Area Interconnected Power System
The blocks used in a single area interconnected power system are generally nonlinear systems.The loads at the operating point and the PV-SPP affect the frequency of the system.A LFC applied to a single area interconnected power system, which is shown in Figure 4, has a turbine, a generator, and other system blocks.Additionally, a PV plant generated through sunlight is demonstrated.For the system, a control responds as quickly as possible.The control process has been realized in two phases, as a primer and a seconder frequency control.The primer frequency control is achieved by adjusting the power balance with the turbine speed regulator, when there is a decrease or increase in electrical power between two different areas.The seconder frequency control is achieved by controlling the power flow between the areas, which reduces the permanent frequency error in the interconnected power system to zero.
For the system, a control responds as quickly as possible.The control process has been realized in two phases, as a primer and a seconder frequency control.The primer frequency control is achieved by adjusting the power balance with the turbine speed regulator, when there is a decrease or increase in electrical power between two different areas.The seconder frequency control is achieved by controlling the power flow between the areas, which reduces the permanent frequency error in the interconnected power system to zero.

Two Area Interconnected Power System
A two area interconnected power system was combined with a connection link.Its block diagram, to be modelled by MATLAB-Simulink, is given in Figure 5, without a controller.In the figure, U 1 and U 2 are the control signals, ∆P tie1 and ∆P tie2 are the line power changes, and ∆ f 1 and ∆ f 2 are the input frequency deviations.A MATLAB-Simulink model of the two area interconnected renewable energy power system is shown in Figure 6, where, ∆ is the frequency of the system, is the regulation constant, ∆ is the load variations, ∆ is the power variation depending on the sunlight of the PV cell, and is the field controller output.The whole system with multi variable parameters is described in terms of Equations ( 9)- (12).Parameters of two-area power system with PV-SPP are shown in Appendix A.
where A is the system matrix; B and L are the input and disturbance matrices; and x(t), u(t), and d(t) are the state, control, and load change disturbance vectors, respectively.The system output depends on the area control error (ACE), shown in Figure 6.The ACE is symbolized by Equation ( 14).
where is the frequency bias constant, ∆ is the frequency deviation, ∆ -is the change in the tie-line power, and C is the output matrix.
is the expression of the frequency change in the power A MATLAB-Simulink model of the two area interconnected renewable energy power system is shown in Figure 6, where, ∆ f i is the frequency of the system, R i is the regulation constant, ∆P Li is the load variations, ∆P PVi is the power variation depending on the sunlight of the PV cell, and U i is the field controller output.The whole system with multi variable parameters is described in terms of Equations ( 9)- (12).Parameters of two-area power system with PV-SPP are shown in Appendix A.
where A is the system matrix; B and L are the input and disturbance matrices; and x(t), u(t), and d(t) are the state, control, and load change disturbance vectors, respectively.
The system output depends on the area control error (ACE), shown in Figure 6.The ACE is symbolized by Equation (14).
Appl.Sci.2017, 7, 308 6 of 22 where b i is the frequency bias constant, ∆ f i is the frequency deviation, ∆P tie1−2 is the change in the tie-line power, and C is the output matrix.b i is the expression of the frequency change in the power unit.The transfer functions used for the modelling of a turbine and a generator are given as follows [10].
Appl.Sci.2017, 7, 308 6 of 21 In the system, a RTU type turbine was employed for increasing the efficiency of these turbines.They are made up of several turbine blades, in which the steam is reenergized between the two turbines by intermediate heaters reheating the steam.They are usually operated in the large plants [24].The RTU transfer function Equation ( 17) is explained as follows [8]: The frequency stability is defined as the ability to maintain a stable frequency after an imbalance between the generated and demanded power.The power balance fails when any difference occurs on the load side or on the generation side [26].Here, represents the load variation of the first area.When ∆ is non-zero, the system frequency changes.The frequency change is defined by the equations given below.In the system, a RTU type turbine was employed for increasing the efficiency of these turbines.They are made up of several turbine blades, in which the steam is reenergized between the two turbines by intermediate heaters reheating the steam.They are usually operated in the large plants [24].The RTU transfer function Equation ( 17) is explained as follows [8]: The frequency stability is defined as the ability to maintain a stable frequency after an imbalance between the generated and demanded power.The power balance fails when any difference occurs on the load side or on the generation side [26].Here, P L1 represents the load variation of the first area.When P ∆1 is non-zero, the system frequency changes.The frequency change is defined by the equations given below. ) where K s is defined as the characteristic constant of the system frequency in the first area.Also, the transfer function expressing the load changes affecting the system frequency is clarified as follows [8]: To identify the ACE value, Equation ( 22) is used.This gives the frequency changes when a load increase of P L1 in the first area occurs for two area interconnected power systems.Equation (23) expresses the power change in the tie line [27].The ACE used in the FLC is calculated by means of Equations ( 22) and ( 23), as follows [24]:

Control Methods
In this study, conventional PID, Fuzzy Logic, GA-PID, and GA-FL controllers have been used to perform the simulations.The sample time was taken as 0.001 s in all simulations.All controllers will be explained below, respectively.

Conventional Controllers
In the literature, proportional (P), proportional-integral (PI), and proportional-integral-derivative (PID) control methods were preferred for load-frequency control in power systems.These studies experienced some difficulties because of the sudden climatic changes and the load fluctuations [28].
In this study, a PID control has been applied to the system and is simulated by MATLAB-Simulink.The output power of PV-SPP was modelled to take into account climatic conditions.The LFC was achieved when the power obtained from the PV-SPP changed.The expression of the control signal (u (t) ) required for the PID controller is given as follows [29]: The Ziegler-Nichols (ZN) method is used for calculating the optimum values of the parameters of the PID controller [30]."K u " is defined as the constant value of oscillation of the system and "T u " is the oscillation frequency.Their formulas are given in Table 1.Many terms used randomly in daily life have a fuzzy structure.The verbal or numerical expressions used when we are describing a thing or an event, commanding, and in many other cases, are also fuzzy.Some of these terms can be said to be young, long, short, hot, cold, hot, cloudy, partly cloudy, sunny, fast, slow, very, small, little, little, too much, etc.When an event is explained and a decision is reached in a situation, terms that do not express such certainty are used.Thanks to the fuzzy logic controller, the examination of systems gained a new perspective.
A FLC fundamentally consists of three parts.These are fuzzification, rule base inferencing, and defuzzification blocks.When modelling the FLC, error and error change signals are employed.A FLC block diagram is illustrated in Figure 7.In this study, in order not to exceed the frequency limits, a designed fuzzy rule table is given in Table 2. Here, a seven stage fuzzy rule table which belongs to triangular membership functions, was formed.The regions in Table 2 were formed as three main regions: positive, zero, and negative.They were placed as two axes adjacent to each other, and their names have been given as ACE and ΔACE.Moreover, they were divided into three parts: large, medium, and small.The fuzzy logic controller has the ability to make decisions with reference to the rules in Table 2.For instance, if ACE is larger than the set value and ΔACE is rapidly increased, the output of the FLC is expected to be large.A Mamdani fuzzy interface was used when arranging the fuzzy logic blocks.Ultimately, the output of the system approached the desired value.All other rules were obtained in a similar way.The weights of the employed triangular membership functions are demonstrated in Figure 8.These weights were found by the trial and error method, to minimize the system error.In this study, in order not to exceed the frequency limits, a designed fuzzy rule table is given in Table 2. Here, a seven stage fuzzy rule table which belongs to triangular membership functions, was formed.The regions in Table 2 were formed as three main regions: positive, zero, and negative.They were placed as two axes adjacent to each other, and their names have been given as ACE and ∆ACE.Moreover, they were divided into three parts: large, medium, and small.The fuzzy logic controller has the ability to make decisions with reference to the rules in Table 2.For instance, if ACE is larger than the set value and ∆ACE is rapidly increased, the output of the FLC is expected to be large.positive, PL: large positive.
A Mamdani fuzzy interface was used when arranging the fuzzy logic blocks.Ultimately, the output of the system approached the desired value.All other rules were obtained in a similar way.The weights of the employed triangular membership functions are demonstrated in Figure 8.These weights were found by the trial and error method, to minimize the system error.

Genetic Optimization
GAs, being part of Evolutionary Algorithms (EA), are based on the creation of new and better solutions as a result of processes such as multiplication, mating, and mutation.They help obtain solutions in a much shorter time than conventional methods.In the same way, with GA, the best parameters of FLC and PID controllers in a sampling space are determined, and the system frequency is more rapidly set to nominal value [18].
Genetic optimization is based on a natural selection process.The algorithm repeatedly develops a population of individual solutions.This method randomly selects individuals from the current population and produces the children for the next generation.The population moves to an optimal solution over successive generations.A genetic algorithm uses three operators to generate the next generation from the current population.These operators are selection, crossover, and mutation.The individuals named as parents are selected.Crossover combines two parents to form children for the next generation.Random changes are applied to individual parents to form children by mutation.Genetic optimization processes are continued via a fitness function, defined and formulated by a GA.The correct determination of this fitness function is the most important step in solving the problem, since the frequency is unbalanced due to the insolation.It emphasizes that the determination stage of each problem has its own solution focus [31,32].

GA-PID Controller
Setting the parameters of a PID controller is the most important part of a control system.In this study, a GA was used to determine the parameters K p , K i , and K d of the PID controller.Then, the GA-PID controller, whose parameters were defined by the GA, was modelled by MATLAB-Simulink.The modelled GA-PID controller flow diagram is shown in Figure 9.

GA-PID Controller
Setting the parameters of a PID controller is the most important part of a control system.In this study, a GA was used to determine the parameters Kp, Ki, and Kd of the PID controller.Then, the GA-PID controller, whose parameters were defined by the GA, was modelled by MATLAB-Simulink.The modelled GA-PID controller flow diagram is shown in Figure 9.The first step includes the forming of the initial population when the PID parameters are adjusted for the LFC.The K p , K i , and K d gain parameters were calculated as applicable to the control system and as individual characteristic forms.The initial operating values were chosen to be zero, before starting the iteration.Since the LFC expresses a balance between the power in the production lines and the demanded power in the load side, the ACE value must be zero [33].For this, the fitness function was defined as in Equation (26).The ACE value obtained at the end of each iteration was measured and sent to the GA.The GA is operated with reference to this value and new K p , K i , and K d parameters were obtained to reduce the error value.At the end of each iteration, the model operated with new K p , K i , and K d parameters attempts to reach a better ACE value than before.This loop was operated for more than 800 iterations.As a result of these iterations, the obtained changes in the K p , K i and K d parameters are given in Figure 10.The GA parameters chosen for PID optimization are shown in Table 3.
As seen in Figure 10, these parameters start to stabilize after 800 iterations.After obtaining the best K p , K i , and K d parameters, the system was operated by using these values and the system results were obtained.The fitness function for the GA-PID controller is given as follows: parameters were obtained to reduce the error value.At the end of each iteration, the model operated with new Kp, Ki, and Kd parameters attempts to reach a better ACE value than before.This loop was operated for more than 800 iterations.As a result of these iterations, the obtained changes in the Kp, Ki and Kd parameters are given in Figure 10.The GA parameters chosen for PID optimization are shown in Table 3.As seen in Figure 10, these parameters start to stabilize after 800 iterations.After obtaining the best Kp, Ki, and Kd parameters, the system was operated by using these values and the system results were obtained.The fitness function for the GA-PID controller is given as follows:  In this study, a new FLC approach (GA-FLC) which has GA optimized weights, has been proposed for the LFC in two area multi source interconnected renewable energy power resources.A block model of the developed GA-FLC could be easily adapted to different systems and could save time when designing these systems.Moreover, for users, the developed GA-FLC could provide the possibility to make changes, through entering the model content into the controller.The control error E is the fundamental element of the two dimensional space DE, which is the variation of this error in a sampling process.The GA parameters chosen for the FLC optimization are shown in Table 3.
The input and output blocks of the GA flow diagram and the FLC model of the new proposed GA-FLC are presented in Figure 11.Here, the first and most important step is to compose these fuzzy membership functions for the input and output sections of the fuzzy space.To carry this out, the equation of the triangular membership function is described, as in Equation (27).
a sampling process.The GA parameters chosen for the FLC optimization are shown in Table 3.
The input and output blocks of the GA flow diagram and the FLC model of the new proposed GA-FLC are presented in Figure 11.Here, the first and most important step is to compose these fuzzy membership functions for the input and output sections of the fuzzy space.To carry this out, the equation of the triangular membership function is described, as in Equation ( 27).

,
, 0 In addition, a MATLAB-Simulink model of the triangular membership function is illustrated in Figure 12.In addition, a MATLAB-Simulink model of the triangular membership function is illustrated in Figure 12.Here, , , and are the precise parameters used to determine the location and region of the triangle [34].This triangle membership function has been operated in both the fuzzification and defuzzification processes.Altas and Sharaf [34] provided detailed information for composing the FLC algorithm via MATLAB-Simulink.
In this study, for the new proposed GA-FLC, the required weights were individually determined for error E, error change DE, and the defuzzification membership functions.For this purpose, the GA was used and the process was repeated in every simulation.The triangular membership functions used these weights and are depicted in Figure 13.Moreover, these changes in the triangular membership function weights and the fitness value of the new proposed GA-FLC are given in Figure 14.Here, X 1 , X 2 , and X 3 are the precise parameters used to determine the location and region of the triangle [34].This triangle membership function has been operated in both the fuzzification and defuzzification processes.Altas and Sharaf [34] provided detailed information for composing the FLC algorithm via MATLAB-Simulink.
In this study, for the new proposed GA-FLC, the required weights were individually determined for error E, error change DE, and the defuzzification membership functions.For this purpose, the GA was used and the process was repeated in every simulation.The triangular membership functions used these weights and are depicted in Figure 13.Moreover, these changes in the triangular membership function weights and the fitness value of the new proposed GA-FLC are given in Figure 14.
Here, , , and are the precise parameters used to determine the location and region of the triangle [34].This triangle membership function has been operated in both the fuzzification and defuzzification processes.Altas and Sharaf [34] provided detailed information for composing the FLC algorithm via MATLAB-Simulink.
In this study, for the new proposed GA-FLC, the required weights were individually determined for error E, error change DE, and the defuzzification membership functions.For this purpose, the GA was used and the process was repeated in every simulation.The triangular membership functions used these weights and are depicted in Figure 13.Moreover, these changes in the triangular membership function weights and the fitness value of the new proposed GA-FLC are given in Figure 14.The new weights , , and have been determined, depending on the error value derived from each iteration and the GA.For the next iteration, the system frequency changes have been observed using these values.The fitness function required for this solution has been governed in a similar manner, by Equation ( 26).The fuzzy logic rules, including the five rules GA-FLC results.The results obtained from these comparisons are given below.When the control process is performed, a model that gives a power increase value of 0.001 Pu in area 2 for the demanded power at the 5th s. and a power increase value of 0.01 Pu in area 1 for the demanded power at the 140th s., were formed.In addition, the fluctuations of the generated power from the solar energy plants inserted in both areas depending on the climatic conditions, starting from the 20th s., were modelled.In Figure 15, these fluctuations of solar energy plants are depicted.

Results and Discussion for Area 1
The modeled system has been operated for 200 s.The changes of ACE1 for area 1 are shown in Figure 16.In Figure 16b, the effect of increasing the demand power in area 2 on the frequency in area 1, is demonstrated.According to the results obtained from the four different controllers, the GA-FLC demonstrates a better performance than the others, in terms of the sudden oscillations.Additively, Figure 16c illustrates the effect of a long-term power drop of the solar power plant due to climatic conditions on the ACE1 value.According to the results of Figure 16, the proposed GA-FLC could set the desired value within certain limits, such that there is hardly any oscillation.In Table 6, the obtained results are presented.The proposed GA-FLC also exhibited the best performance, with a settlement time of 5.5 s. from the 70th s, as indicated in Table 6.It is also seen that the GA-FLC has a better value as 0.06 10 , when compared to other controllers at the 70th s.In Figure 17, variations of the absolute value of the ACE1 value for area 1; (a) PID; (b) GA-FLC; (c) GA-PID; and (d) FLC, are demonstrated.The reason why the absolute value of the ACE1 is displayed is that the reference error variance of the system used in the GA processes could be seen from zero, to how much it changes in the positive direction.In addition, the ACE1 signal, which is controlled by four different controllers, was converted to the form, , in Figure 18, and was then transferred to the governor and the turbines to balance the frequency value.The purpose of this is to observe which signals sent to the turbines and the generators are more effective.The signals obtained from the controllers are shown in Figure 18.As shown in Figure 18b, GA-FLC produced a wider and higher frequency value.This is due to the use of more mathematical functions while modelling the GA-FLC in

Results and Discussion for Area 1
The modeled system has been operated for 200 s.The changes of ACE 1 for area 1 are shown in Figure 16.In Figure 16b, the effect of increasing the demand power in area 2 on the frequency in area 1, is demonstrated.According to the results obtained from the four different controllers, the GA-FLC demonstrates a better performance than the others, in terms of the sudden oscillations.Additively, Figure 16c illustrates the effect of a long-term power drop of the solar power plant due to climatic conditions on the ACE 1 value.According to the results of Figure 16, the proposed GA-FLC could set the desired value within certain limits, such that there is hardly any oscillation.In Table 6, the obtained results are presented.The proposed GA-FLC also exhibited the best performance, with a settlement time of 5.5 s. from the 70th s, as indicated in Table 6.It is also seen that the GA-FLC has a better value as 0.06 × 10 3 Pu, when compared to other controllers at the 70th s.In Figure 17, variations of the absolute value of the ACE 1 value for area 1; (a) PID; (b) GA-FLC; (c) GA-PID; and (d) FLC, are demonstrated.The reason why the absolute value of the ACE 1 is displayed is that the reference error variance of the system used in the GA processes could be seen from zero, to how much it changes in the positive direction.In addition, the ACE 1 signal, which is controlled by four different controllers, was converted to the form, U (t)1 , in Figure 18, and was then transferred to the governor and the turbines to balance the frequency value.The purpose of this is to observe which signals sent to the turbines and the generators are more effective.The U (t)1 signals obtained from the controllers are shown in Figure 18.As shown in Figure 18b, GA-FLC produced a wider and higher frequency value.This is due to the use of more mathematical functions while modelling the GA-FLC in a Simulink environment.With frequency damping blocks, this high frequency signal can be prevented from damaging the turbine valves.Moreover, Figure 19 shows that the frequency value in the interconnected network controlled by the proposed GA-FLC is better than the other control methods.The effect of a load increase in area 2 on area 1 was shown to be minimized.Furthermore, the frequency oscillations of the proposed GA-FLC were shown in this way, which shows the frequency value in the interconnected network, to be kept at the minimum value.The best results obtained by GA-FLC have also been reached by a settling time of 5.1 s and an overshoot value of −0.00018 P u , as indicated in Table 6.
the interconnected network controlled by the proposed GA-FLC is better than the other control methods.The effect of a load increase in area 2 on area 1 was shown to be minimized.Furthermore, the frequency oscillations of the proposed GA-FLC were shown in this way, which shows the frequency value in the interconnected network, to be kept at the minimum value.The best results obtained by GA-FLC have also been reached by a settling time of 5.1 s and an overshoot value of −0.00018 Pu, as indicated in Table 6.   a Simulink environment.With frequency damping blocks, this high frequency signal can be prevented from damaging the turbine valves.Moreover, Figure 19 shows that the frequency value in the interconnected network controlled by the proposed GA-FLC is better than the other control methods.The effect of a load increase in area 2 on area 1 was shown to be minimized.Furthermore, the frequency oscillations of the proposed GA-FLC were shown in this way, which shows the frequency value in the interconnected network, to be kept at the minimum value.The best results obtained by GA-FLC have also been reached by a settling time of 5.1 s and an overshoot value of −0.00018 Pu, as indicated in Table 6.

Results and Discussion for Area 2
Detailed simulation results of the interconnected network in area 2 are given in Figures 20-22.The variations of ACE2 values throughout the entire simulation are illustrated in Figure 20a and the special cases are shown in Figure 20b,c.Here, the case of both sudden load fluctuations and the responses of the system variables for long-term power reduction in solar power plants, are observed.20a and the special cases are shown in Figure 20b,c.Here, the case of both sudden load fluctuations and the responses of the system variables for long-term power reduction in solar power plants, are observed.As noted in Table 6, the proposed GA-FLC reached a settling time value of 5.1 s and then quickly minimized the error.In addition, the proposed GA-FLC provided at least three times more overshoot than the GA-PID controller, as shown in Figure 20b and Table 6.In Figure 21, U (t)2 , as the new status after controlling the ACE 2 value, is demonstrated.The results of the PID, GA-PID, FLC, and the proposed GA-FLC controllers modelled to minimize the frequency variation, ∆F 2 , in area 2, are revealed in Figure 22.In Figure 22b, it is shown that the GA-FLC maintains a balanced frequency in the second region.The power change over the connection line between area 1 and area 2 is shown in Figure 23.As shown in Figure 23 and Table 6, the frequency deviation in the second area is minimized.When the situation at the 5th s is examined, the settling time value of the GA-FLC is set to 5.5 s, GA-PID is set to 5.9 s, FL is set to 8.3 s, and PID is set to 9.6 s.
When the responses given by the controllers to the situation caused by the decrease in power at the 70th s in the solar power plant are examined, it is seen that the settling time value of 70.1 s obtained by the proposed controller is the fastest and best result.
region.The power change over the connection line between area 1 and area 2 is shown in Figure 23.As shown in Figure 23 and Table 6, the frequency deviation in the second area is minimized.When the situation at the 5th s is examined, the settling time value of the GA-FLC is set to 5.5 s, GA-PID is set to 5.9 s, FL is set to 8.3 s, and PID is set to 9.6 s.
When the responses given by the controllers to the situation caused by the decrease in power at the 70th s in the solar power plant are examined, it is seen that the settling time value of 70.1 s obtained by the proposed controller is the fastest and best result.the situation at the 5th s is examined, the settling time value of the GA-FLC is set to 5.5 s, GA-PID is set to 5.9 s, FL is set to 8.3 s, and PID is set to 9.6 s.
When the responses given by the controllers to the situation caused by the decrease in power at the 70th s in the solar power plant are examined, it is seen that the settling time value of 70.1 s obtained by the proposed controller is the fastest and best result.

Conclusions
In this paper, a new GA-FLC controller was proposed for frequency stability problems in a two area multi source interconnected power system.The power system to be controlled was consisted of a thermal, a reheated thermal, and a PV unit.In addition to the proposed controller, conventional PID, FL, and GA-based PID controllers were designed.All systems and controllers were designed and simulated with the Matlab-Simulink program.In the simulations, a change in the frequencies and the tie-line power were observed by implementing each controller.Since there was a two-area power system and they were connected by a tie-line, the outputs of ∆ , ∆ , and ∆ of the system were individually plotted.As can be seen from the data in Table 6, the proposed controller produces better results than the other controllers, both in terms of overshoot and settling time values.These values were taken during the load increase and the solar power decrease periods.Therefore, this demonstrates the reaction of the system at the moment of operation and a more realistic simulation was realized.The proposed GA-FLC has three membership functions with an optimized GA, and twenty-five rules.Owing to the proposed controller, the system frequency is affected at a minimum rate when the power changes on the PV-SPP because of climatic conditions.Moreover, the ability of the system to adapt to sudden power fluctuations at the operation point was also improved.In the meantime, the optimum values of the , , and parameters for the PID controller were determined by the Zigler-Nichols method.

Conclusions
In this paper, a new GA-FLC controller was proposed for frequency stability problems in a two area multi source interconnected power system.The power system to be controlled was consisted of a thermal, a reheated thermal, and a PV unit.In addition to the proposed controller, conventional PID, FL, and GA-based PID controllers were designed.All systems and controllers were designed and simulated with the Matlab-Simulink program.In the simulations, a change in the frequencies and the tie-line power were observed by implementing each controller.Since there was a two-area power system and they were connected by a tie-line, the outputs of ∆F 1 , ∆F 2 , and ∆P tie of the system were individually plotted.As can be seen from the data in Table 6, the proposed controller produces better results than the other controllers, both in terms of overshoot and settling time values.These values were taken during the load increase and the solar power decrease periods.Therefore, this demonstrates the reaction of the system at the moment of operation and a more realistic simulation was realized.The proposed GA-FLC has three membership functions with an optimized GA, and twenty-five rules.Owing to the proposed controller, the system frequency is affected at a minimum rate when the power changes on the PV-SPP because of climatic conditions.Moreover, the ability of the system to adapt to sudden power fluctuations at the operation point was also improved.In the meantime, the optimum values of the K p , K i , and K d parameters for the PID controller were determined by the Zigler-Nichols method.
In conclusion, the frequency stability problem caused by solar energy sources in the interconnected power systems was abolished.According to these results, for the frequency stability, the proposed GA-FLC yielded an efficiency of about 90% compared to the conventional PID controllers.Moreover, it was observed that the overshoot value of the proposed control technique was 75% better than the FLCs, and approximately 40% better than the GA-PID controllers.Ultimately, since the the overshoot values and the settling time values directly affect the operating lifetimes, operating costs, and efficiencies of the grid, the proposed GA-FLC would be beneficial and advisable for an interconnected power system with solar power sources.

Figure 1 .
Figure 1.A two area three source interconnected system block diagram.

Figure 1 .
Figure 1.A two area three source interconnected system block diagram.

Figure 2 .
Figure 2. The simplified equivalent circuit of a PV cell.

Figure 2 .
Figure 2. The simplified equivalent circuit of a PV cell.

Figure 4 .
Figure 4.A LFC applied a single area interconnected power system.Figure 4. A LFC applied a single area interconnected power system.

Figure 4 .
Figure 4.A LFC applied a single area interconnected power system.Figure 4. A LFC applied a single area interconnected power system.

3 .
Two Area Interconnected Power SystemA two area interconnected power system was combined with a connection link.Its block diagram, to be modelled by MATLAB-Simulink, is given in Figure5, without a controller.In the figure, U1 and U2 are the control signals, ∆ and ∆ are the line power changes, and ∆ and ∆ are the input frequency deviations.

Figure 5 .
Figure 5. Block diagram of an uncontrolled two area interconnected power system.

Figure 5 .
Figure 5. Block diagram of an uncontrolled two area interconnected power system. .

Figure 6 .
Figure 6.Two area multi source interconnected with PV power system grid in MATLAB-Simulink block diagram.

Figure 6 .
Figure 6.Two area multi source interconnected with PV power system grid in MATLAB-Simulink block diagram.

Figure 8 .
Figure 8.The membership functions for FLC.N L : large negative, N M : medium negative, N S : small negative, Z: zero, P S : small positive, P M : medium positive, P L : large positive.

Figure 9 .
Figure 9. Flow diagram for determining the PID parameters of GA.Figure 9. Flow diagram for determining the PID parameters of GA.

Figure 9 .
Figure 9. Flow diagram for determining the PID parameters of GA.Figure 9. Flow diagram for determining the PID parameters of GA.

Figure 10 .
Figure 10.Changes in the GA-PID parameters during the iteration process; (a) Kp, (b) Ki and (c) Kd.

Figure 10 .Table 3 .
Figure 10.Changes in the GA-PID parameters during the iteration process; (a) K p , (b) K i and (c) K d .

Figure 11 .
Figure 11.The flow diagram of the new GA-FLC optimized its weights with a GA.

Figure 11 .
Figure 11.The flow diagram of the new GA-FLC optimized its weights with a GA.

Figure 12 .
Figure 12.A MATLAB-Simulink model of the triangular membership function.

Figure 12 .
Figure 12.A MATLAB-Simulink model of the triangular membership function.

Figure 13 .
Figure 13.The triangular membership functions employed the new proposed GA-FLC.Figure 13.The triangular membership functions employed the new proposed GA-FLC.

Figure 14 .
Figure 14.Changes in the triangular membership function weights and the fitness value of the new proposed GA-FLC; (a) W re f −DE (b) W re f −E ; (c) fitness value.

Figure 15 .
Figure 15.Power variations in the solar power plant and areas and: (a) solar insolation; (b) solar power; (c) area 1 demand power; (d) area 2 demand power.

Figure 15 .
Figure 15.Power variations in the solar power plant and areas and: (a) solar insolation; (b) solar power; (c) area 1 demand power; (d) area 2 demand power.

Figure 19 .
Figure 19.The detailed display of frequency, ΔF1 change in area 1.

Figure 19 .
Figure 19.The detailed display of frequency, ∆F 1 change in area 1.

Figure 20 .
Figure 20.Changes in the ACE 2 in the area 2: (a) during all simulation time; (b) 5th change; (c) 18th change.

Figure 23 .
Figure 23.The power change over the connection line between area 1 and area 2.

Figure 23 .
Figure 23.The power change over the connection line between area 1 and area 2.
d(t): load changes disturbance vector R s : Series resistance(Heat loses) G L : The Loads transfer function R G : Regulation constant G T (s): The turbines transfer function R R : Regulation constant with re-heat G R : Re-heat turbine transfer function S c : Benchmark reference solar irradiation G T : Generator transfer function T 1−2 : Synchronization constant G G (s): The generators transfer function T G : Speed governor time constant G G : Turbine transfer function T P : Control area time constant I 0 : Reverse saturation current of diode (A) T R : Re-heat time constant I D : PV solar cell formed in the reverse saturation current T T : Turbine time constant I ph : Photocurrent (A) T cell T c : Reference cell operating temperatureI cell : Cell output current (A) U 1 , U 2 :1 and 2 areas controller outputs k: Boltzmann constant (J/ • K) u(t): control vector K S : Characteristic constant of the system frequency V cell : Cell output voltage (V) K P : Control area gain x(t): state vector K R : Re-heat gain y(t): system output vector P L1 , P L2 : Changing demand power of area 1 and 2 ∆P L : Load increase P 1 , P 2 : The total power in the first region ∆ f 1 , ∆ f 2 : 1 and 2 area frequency changes P G1 : The power of generator power system ∆P tie1−2 : Change in tie-line power between two areas P L1 : Refers to the change in load first region ∆ω: Frequency changes

Table 1 .
ZN method of detection parameters.

Table 2 .
Fuzzy logic rules for FLC.

Fuzzy Logic Rules for FLC ACE ∆
NL: large negative, NM: medium negative, NS: small negative, Z: zero, PS: small positive, PM: medium positive, PL: large positive.

Table 2 .
Fuzzy logic rules for FLC.

Table 3 .
Control Parameters of GA-PID and GA-FLC.

Table 6 .
The obtained results.

Table 6 .
The obtained results.Detailed simulation results of the interconnected network in area 2 are given in Figures20-22.The variations of ACE 2 values throughout the entire simulation are illustrated in Figure