Rational Application of Electric Power Production Optimization through Metaheuristics Algorithm

: The aim of this manuscript is to introduce solutions to optimize economic dispatch of loads and combined emissions (CEED) in thermal generators. We use metaheuristics, such as particle swarm optimization (PSO), ant lion optimization (ALO), dragonﬂy algorithm (DA), and differential evolution (DE), which are normally used for comparative simulations, and evaluation of CEED optimization, generated in MATLAB. For this study, we used a hybrid model composed of six (06) thermal units and thirteen (13) photovoltaic solar plants (PSP), considering emissions of contaminants into the air and the reduction in the total cost of combustibles. The implementation of a new method that identiﬁes and turns off the least efﬁcient thermal generators allows metaheuristic techniques to determine the value of the optimal power of the other generators, thereby reducing the level of pollutants in the atmosphere. The results are presented in comparative charts of the methods, where the power, emissions, and costs of the thermal plants are analyzed. Finally, the comparative results of the methods were analyzed to characterize the efﬁciency of the proposed algorithm.


Introduction
The energy crises that occur worldwide can be managed through the connection of sustainable and renewable energy systems, which are needed as we attempt to reduce the use of fossil fuels as populations increase. Interconnected grids can be divided into a power system structure based on the joint operations to generate and transmit power to the load demand as operational and technical controls [1][2][3].
The standard economic load dispatch (ELD) solution seeks to assign the total power demand among all generators used to achieve the minimum fuel cost [4]. Deterministic methods applied in the ELD solution are difficult to apply because of the non-continuous, non-convex, and nonlinear nature of the problem [5,6]. However, new rules have forced Thermoelectric Power Plant (TPP) to reduce the amount of polluting emissions into the atmosphere, expanding the ELD issue to the economic combined emissions dispatch (CEED) created to minimize fuel costs and pollutant emissions such as NO x , SO x , and CO x elements from the TPP. Although the ELD and CEED problems are nonlinear optimization problems, many heuristic methods have been implemented to solve them, CEED is a multiobjective optimization [7,8]. The concept of the price penalty factor was proposed by some researchers to transform the CEED multi-objective issue into a single-objective issue by unifying the emission cost equations with the fuel cost equations [9,10].
Many techniques have been proposed to solve the ELD problem in power systems [11][12][13], and a nonconvex ELD difficulty has been addressed by several hybrid optimization tech-

Problem Formulation
To streamline the multicriteria issue of CEED, two objective functions must be considered: one for oil consumption f 1 (P i ) and one for environmental cost f 2 (P i ) [29]. The objective function coefficients are achieved by adjusting the techniques of the curves based on the generator performance test [30,31].

Mathematical Model of Generation by Solar Power Plant (g 1 )
The model of solar power plant generation (SPP) is described in Equation (1) [8,[32][33][34]: g 1 Pgs j = P rated 1 + T re f − T amb * al pha * S i 1000 (1) where: g 1 = solar power plant Pgs j = generated power by the solar plant P rated = rated power; T ref = reference temperature; T amb = room temperature; alpha = temperature coefficient; and S i = incident solar radiation.
With the contribution of the SPP stemming system, the added solar energy is described in (2) [8,[32][33][34]: where Pgs j is the energy available at the jth SPP and Us j indicates the status of the jth SPP, which is 1 (ON) or 0 (OFF). The cost of the SPP is given by (3) [8]: where PUCost j is the unit cost of jth SPP.

Mathematical Model of Cost of Thermal Plants (f 1 )
The principal ELD characteristic is to ascertain the optimum energy demand between the generating units of the station, minimize the total operational cost, and satisfy one cluster of restrictions of equality and inequality, which is one of the principal tasks of energy system optimization. Due to environmental responsibility, ELD was devoted to a multi-criterion optimization issue of CEED, aiming to reduce the emission of pollutant gases in the atmosphere [35].
The equation for the fuel consumption of every unit is represented by the quadratic function in Equation (4), considering the generator output power P i , given as $/h, as in [4]: where, a i b i and c i represent the fuel consumption coefficients of each generating unit.
The problem of minimizing the total cost of a TPP is represented by Equation (5): where N is the total number of generating units of the TPP and P i is the output power of each generating unit.

Mathematical Model of Emissions of Thermal Plants (f 2 )
The function of the total thermal station emissions formulated in Equation (6) relates the emission to the power generated by every generator unit. This function represents the emission of SO 2 and NO x in kg/h, which can be expressed as [4,9,29,36]: where α i β i and γ i represent the emission coefficients of all generating units. Equation (7) represents the total minimizations of the TPP problem.
where N is the total number of generators in the TPP and P i is the energy output of each generator.

Equal Power Constraints
The Equation (8) is expressed as a parity constraint on the nominal power, given the lower and upper bounds of each unit, following [37,38]: where P i is the nominal power, P D is the power demand, P L is the transmission loss. Therefore, the total generation should be similar to the power demand, plus the actual loss in the transmission lines. According to Equation (9): ∑ n i=1 P i = P D + P L (9) The restriction on the total power generated in Equation (8) considers the total power generated by the SPP, that is, Equation (10) [8]. The P L dimensioning is equal to the sum of losses versus power, which have the same active and reactive power constraints at each bar according to Equation (11) [25]: Transmittance losses are defined as an output generator function arising from the Kron loss coefficients of Kron's loss formula, Equation (12) [38,39]: where B ij , B oi , and B oo are the coefficients of the system transmission loss, and n is the generator number. The B coefficients can be accurately identified when the actual operating conditions are close to the base case [38,39].

Generation Constraint
Equation (13) expresses the power of each generating unit by the upper and lower limits, according to [31,40,41]: where P i is the output power of each generating unit and P min.i and P max.i are the minimum and maximum power of the generator.
Another constraint that has to be met in a mixed system with photovoltaics is to maximize the utilization to 80% of SPP capacity installed, owing to the generator instability, presented in Equation (14) [8]:

Optimization Problem
For this issue, the energy generation is calculated by analyzing the SPP capacity installed; maximum value is being utilized to 80% of the total generation, and then the SPP energy generation is calculated, applying the function g 1 (Pgs j ) described in (1). The remaining pendent demand is solved by minimizing the multi-criteria problem in Equation (15): where f 1 (P i ) and f 2 (P i ) are the objective functions of the costs and emissions, respectively, to be optimized.

Formulation of the Incremental Cost
The incremental cost of oil ($/MWh) is given by Equation (16) [42]: where: IC i = incremental fuel cost; a i = actual incremental cost curve; b i = is an approximate (linear) incremental cost curve; P gi = total power generation.

Optimization Technique
Despite using different techniques and methods to solve the CEED, new and innovative algorithms are still being proposed for the solution. The main intrinsic motivation for this is the theorem No Free Lunch (NFL) [43]. The NFL theorem states that a single algo- rithm does not work as accurately for all optimization problems. Therefore, it is justified to search for new, more suitable algorithms and to refine existing algorithms.
To experiment with this theory and seek new results, we used the DA, DE, and ALO techniques to compare the results of the three techniques to determine the one that can obtain better solutions to the CEED problem.
A new swarm intelligence optimization technique called the DA was proposed by [44]. This technique considers the proposal of binary and multi-objective versions. Dragonflies are small predators that hunt almost every small insect in nature. An interesting fact about dragonflies is their unique and rare swarming behavior. Dragonflies group together for two purposes: the first being hunting (static) and the second being migration (dynamic swarm). The hunting is in small groups. Dynamic swarms form large groups and travel long distances [44][45][46][47]. To point out three primitive principles of swarm behavior [44,47]: Separation is about preventing static collisions of individuals with other individuals in the neighborhood. The alignment indicates the similar speed of the individuals with that of other individuals in the neighborhood. Cohesion refers to the predisposition of individuals to move toward the center of mass in the neighborhood.
For any swarm its survival is necessary, so all individuals must be attracted to its food sources to avoid confrontation with external enemies.
Each of these behaviors is mathematically modeled as follows: The separation is calculated as follows [44,47]: where X is the current individual, X j is the position of the jth neighboring individual, and N is the number of neighboring individuals. The positioning is calculated by Equation (18): where V j denotes the velocity of the jth neighboring individual. Cohesion is calculated as follows: where X is the position of the current individual, N is the number of neighborhoods, and X j is the position of the jth neighbor individual. Attraction to a food source is calculated as follows: where X is the location of the current individual and X + is the position of the food source.
The Equation (21) shows the division to the outside of an enemy individual: where X is the position of the current individual and X − shows the position of the enemy. The dragonfly assumes a behavior based on five patterns. To update the positions of the artificial dragonflies at a search location and simulate their displacements, two vectors were considered: pitch (∆X) and position (X). The pitch vector is similar to the velocity vector in PSO, and the DA algorithm is developed based on the structure of the PSO algorithm. The pitch vector shows the direction of displacements of dragonflies; however, the introduced method can be extended to higher dimensions [44,46,47].
For DE, [48,49] proposed a new encoded evolution algorithm using a fluctuating point for global optimization and was nominated as DE algorithm. DE has four main stages: initialization, mutation, crossover, and selection. Many optimization guidelines should also be adjusted. These guidelines are bonded under the control guidelines of the common name. There are only three guidelines for the real control of the algorithm: the differential constant F (or mutation), crossover constant Cr, and population length Np. The remaining guideline dimensions of issue D measure optimization task difficulty, and generation maximum number (or iteration) Gen. This can be suited as an interruption condition, and low and high limit restrictions are variables that range a viable area [41,48,50].
The ith vector of current population G with size D can be defined by: Initialization: The initialization of the DE uses Np D-dimensional real-valued parameter vectors to perform random generation of candidate solutions.
where rand ij [0, 1] is a random number, 0 ≤ rand ij [0, 1] ≤ 1 which is multiplied by the interval length, X j,max − X j,min ensures a distributed sampling of the parameter's domain interval X j,min , X j,max . Different approaches with random uniformity are the most common, so they are used to generate the initial population [50,51].
Mutation: Mutation: Differential mutation adds a random-scale vector subtraction to a third vector. The vectors of the variable → V i,G , also called donors, are obtained by this operation [50,51]: where F > 0 is a real number that controls the growth rate of a population. Vectors → X r i 1 ,G , → X r i 2 ,G and → X r i 3 ,G from the current population the samples are randomly, and r i 1 , r i 2 , r i 3 Are integers, respectively, are chosen from the given interval {1, . . . , N p}. In classical DE, small F-values are associated with the analysis, which is implied as a strategy if some of the test solutions are in the area of the global minimum. On the other hand, large F-values are associated with research, since new variable experimental solutions (donors) incorporate larger differences from the original population (targets) [50,51]. Crossover: Increases the potential disparity of the population. In the case of the crossover binomial, the vectors for evaluation U i,G are produced according to [50,51]: Selection: Selection may be understood as a form of competition, in line with many examples that are directly observable in nature. Many evolutionary optimization schemes, such as DE or GAs, use some form of selection [50].
For a selection operation, pairwise selection (elitist selection) is used in the algorithm. The maximum number of generations G max is defined as a stopping criterion [50,51].
Ant lion, the optimization technique of ant lion, is a stochastic research algorithm based on a recently developed population propounded by [52] to solve problems of restricted optimization engineering issues. ALO is inspired by the lifespan of ant lions (doodlebugs), which belong to the family Myrmeleontidae and the order Neuroptera (grid of insects with wings). This technique is a free algorithm of the grid and lacks a baseline for adjustment. Because ALO is a population-based algorithm, the avoidance of ideal places is inherently high. The ALO algorithm has a high probability of solving ideal place closeness because of the random loops and roulette swivels. The pursuit of space exploration in the ALO algorithm is guaranteed by random selection of ant lions and random loops of ants all around it, and the pursuit of space exploration is guaranteed by the adaptive shrinkage limits of ant lion traps [53]: The ants random loop Building traps Ants entanglement with traps Catching preys and Traps restoration.
Adult and larval stages are both important stages of the ALO lifespan. Ant lions prey on the larval and adult phases. The larval period was based on the ALO algorithm. It digs the sand into a cone-shaped pit, following a circular path, and displaces the sand with its jaw. After the trap is built, the larva waits for its prey. The trap size differs depending on the maggots' hungriness, ant lion, and moon size [53].
The ALO algorithm can also be optimized to issue an EELD [35,53]. This algorithm is effectively applied to solve many categories of test functions (non-modal, multi-model, and compounded). The ALO technique is converted to a fast-converted general resolution owing to the use of the roulette selection method, which is also read with optimization issues that remain discrete. With different application problems, the use of ALO technique is valid in comparison with other algorithms such as PSO, DE and DA.
As the ants transit on the pursuit space randomly to find prey, a random loop is selected to demonstrate the ant's effort in Equation (27) [35,52,53]: (27) where cumsum computes the accumulated total, n is the maximum number of iterations, t is the current iteration, and r(t n ) is the stochastic function described in Equation (28) [35,52,53]: To maintain the random walks of ants inside the search space, the positions of their walks are normalized using the following min-max normalization described in Equation (29): where a i and b i are the minimum and maximum of the ith variable of random iterations, respectively, c t i and d t i indicate the minimum and maximum ith iteration in the tth variable, respectively.
The mathematical expression of the trapping of ants in the ant lion pits is given by Equations (30) and (31) as follows: The fittest ant lions were selected using the roulette-wheel method. The sliding of ants into pits is given by Equations (32) and (33) [35,52,53]: where c t is the minimum of all variables at the tth iteration, and d t the vector including the maximum of all variables at the tth iteration, and I is a ratio according to Equation (34): where t is the current iteration, S is the maximum number of iterations, and w is a constant given by Equation (35) [54]: The ant lion captures the ant when it hits the rock bottom and consumes it. Subsequently, the ant needs to update its position to capture the new prey. This process is represented by Equation (36).
On what Antlion t i indicates the position jth to the ant lion selected in the iteration, Ant t i displays the ant position in the iteration, and t indicates the current iteration. Elitism was used to preserve the best solutions in each procedure. The lion ant result is treated as elite, where the technique is more apt. The elite should affect the ant at each step (random movement). For this, each ant is assigned a breeder by the roulette wheel and elite, which is given by Equation (37) [35]: where R t A is the random path around the ant lion selected by the roulette method in the tth iteration, R t B is the random trajectory around the elitism of tth e Ant t i indicates the positioning of the ith ant in the tth iteration.

Applied Procedures to Solve the CEED Problem
Actions used to solve the CEED problem: The first step is the formation of an objective function of the CEED based on the optimization of Equations (5) and (7); the second is a power integration system, considering the IEEE-6 bus system; the third step is to apply the DA, DE, and ALO metaheuristic techniques to implement their usage guidelines; and the fourth is the progress of the technique used to obtain results from the energy system.
The CEED impacts were managed as a simultaneous multi-objective optimization problem. To solve this problem, the techniques must reduce the total cost of fuel consumption and pollutant emissions from the generating units [35,55,56].
We used three meta-heuristic techniques to compare the results obtained and evaluate which achieved the best result for the CEED problem. In the following example, details of the implementation of the ALO are provided.
Optimization steps: Step 1. The main agent of the ALO search, are characterized by the set of ants with random values.
Step 2. The capacity value of each ant is evaluated using an objective function (Equation (15)) for each iteration.
Step 3. The ants' random paths through the search space are expected by the ant lion ant traps.
Step 4. The position of the ants are evaluated in each iteration and the ones in the best position are relocated to capture the others.
Step 5. The Lion ant is more agile, as it needs its position updated to catch the ant that becomes fitter. Step 6. An elite ant lion can affect the movement of the other ants, regardless of their displacement.
Step 7. If a lion ant becomes better than the elite, then it is replaced by the new aptitude.
Step 8. Steps 2 to 7 are repeated until the final parameter is satisfied.
Step 9. The position and fitness coordinates of the elite ant lion are replicated as the best inferences for the overall optimization.
Random trajectories of the ants using Equation (27). In addition to ants, it can be assumed that ants are hidden somewhere in the search space; to save their positions and aptitude values, the following matrices are used (38)-(41) [35,52]: where M Ant is the matrix that stores each ant lion position, AL i,j is the value of the jth dimension of the ith ant lion, n is the number of ant lions, and d is the number of variables.
To evaluate the generating units (ants), objective functions referred to in Equations (4) and (5) are used in developing the optimization and matrix (38), to keep the value of the characteristics of all generators: . . .
where M OA is the matrix that maintains the fitness of each ant, Ant i,j shows the value of the jth dimension of the ith, n is the number of ants, and f is the objective function.
To optimize the cost and power generation, matrices (40) and (41) are used: where M AL stores the position of each ant, AL i,j shows the ALO proportion value j, n is the number of ants, and d is the number of variables (i.e., generators).
where M OAL stores the properties of each ALO, AL i,j is the value of the jth dimension of the jth ALO, n is the number of ants, and f is the objective function of the optimization. The number of generations of the system involving the solution that will be optimized shows the results for minimizing the costs and emissions of the contaminant gases, where it meets the constraints of Equations (8)- (13).
Equation (15) was applied to assess the performance of the CEED, although optimal costs and emissions were achieved. For inequality restrictions analogous to other techniques, the solutions achieved for any iteration are off-limits, causing the algorithm to continue until the maximum iteration and best results are obtained [57]. A flowchart of the ALO method is shown in Figure 1. The number of generations of the system involving the solution that will be optimized shows the results for minimizing the costs and emissions of the contaminant gases, where it meets the constraints of Equations (8)- (13).
Equation (15) was applied to assess the performance of the CEED, although optimal costs and emissions were achieved. For inequality restrictions analogous to other techniques, the solutions achieved for any iteration are off-limits, causing the algorithm to continue until the maximum iteration and best results are obtained [57]. A flowchart of the ALO method is shown in Figure 1.

Case Study: IEEE 6-Units Test System and 13 Solar Plants
The proposed model satisfied 100% of the demand. For safety reasons, 80% or less of the SPPs capacity will be used; owing to the instability in energy capitation by the SPP generators, the rest of the demand will be met by TPP. In terms of the percentage of TPP to be served, optimization techniques were applied to solve the CEED problem, considering a test system with six units to meet some demand needs [8]. The selected TPP case study is composed of (six) generating units, which present fuel consumption coefficients (a, b, and c) and minimum (Pmin) and maximum (Pmax) power limits, as listed in Table  1.

Case Study: IEEE 6-Units Test System and 13 Solar Plants
The proposed model satisfied 100% of the demand. For safety reasons, 80% or less of the SPPs capacity will be used; owing to the instability in energy capitation by the SPP generators, the rest of the demand will be met by TPP. In terms of the percentage of TPP to be served, optimization techniques were applied to solve the CEED problem, considering a test system with six units to meet some demand needs [8]. The selected TPP case study is composed of (six) generating units, which present fuel consumption coefficients (a, b, and c) and minimum (Pmin) and maximum (Pmax) power limits, as listed in Table 1.    Table 3 presents the power ratings and unit costs of the various SPPs, which are estimated to be within the established range.  Table 4 covers global solar radiation, exactly like the temperature profiles and Islamabad charge for 17 July 2012. The solar radiation aggregated data were generated using the Geospatial Toolkit, the data related to the energy demand of the Islamabad region were taken off IESCO, and the temperature profile was obtained from [8]. Data on 17 July 2012 were selected, arbitrarily, from the only data available on that day.

Analysis and Discussion of Results
The results obtained from the three simulations were compared with those obtained in [8]. The proposed model was developed in a computational environment using a Core i7 machine, 8GB ram, win10 operating system with MATLAB R2015a software. The control settings used in the simulation were C1 and C2 = 2, r1 and r2 were random numbers between 0 and 1, and the maximum number of iterations was 1800, obtaining the best results for 1500, with a fixed initial population of 500 individuals considering a significant increase in computational cost for a larger number of individuals without major contributions to the results. Table 5 presents the results obtained in the simulation at 10:00 am, which reached the total energy demand of 1244 energy. When evaluating Table 5, it is noted that the demand of 1244 MW required for the schedule is achieved by the three optimization techniques, and the algorithm intelligently verifies the power of the six generators, shutting down with the less powerful ones. With the application of this technique, we were able to meet the demand efficiently, guaranteeing all the operating parameters of the system, and reducing fuel consumption and pollutants in the atmosphere.
We verified that the generator turned off in the simulations was the P2 UG. Even so, the energy generation did not suffer great variations, providing the stable energy generation to the system, in addition to enabling its predictive maintenance. On the other hand, the P1 UG generator presented lower power in the simulations, but in the DE technique, it exhibited better performance. In the P3 UG generator, the power variations were more notable in relation to the P1 UG generator when applied to the three techniques. The same occurs for the other P4, P5, and P6 UG generators.
We obtained the average fuel cost of the techniques (ALO, DA, and DE) compared to the Khan [8] technique, applied at 10:00 a.m. to meet the demand of 1244 MW, and achieved a reduction of approximately 6.58% corresponding to $3462.32 in the cost of fuel with the use of PV solar energy; it was possible to reduce pollutants in the atmosphere (ALO 21.58%, DA 21.69%, and DE 21.67%). The average reduction in pollutants in the atmosphere was approximately 21.64%. Figures 2-5 show the results of the data obtained in Table 5, which represent the values of power, emissions and fuel costs generated by the techniques (ALO, DA and DE) at 10:00 a.m. demand. phere (ALO 21.58%, DA 21.69%, and DE 21.67%). The average reduction in p the atmosphere was approximately 21.64%. Figures 2-5 show the results of the data obtained in Table 5, which re values of power, emissions and fuel costs generated by the techniques (AL DE) at 10:00 a.m. demand.             Table 6 shows the results of the comparative simulation at 11:00 a.m., reaching a total power demand of 1088 MW. When evaluating Table 6, it is noted that the demand of 1088 MW required for the hour is achieved in the three optimization techniques, and the algorithm intelligently verifies the power of the six generators by turning off the lowest power generators. With the application of this technique, we were able to meet the demand efficiently, guaranteeing all the operating parameters of the system and reducing fuel consumption and pollutants in the atmosphere.
We verified that the generators turned off in the simulations were P2 and P3 UG; even so, the energy generation did not suffer great variations, providing the stability of the energy generation to the system, in addition to enabling predictive maintenance as the generators are turned off. On the other hand, the P1 UG generator presented lower power in the simulations, but in the DE technique, it exhibited better performance. In the P4 UG generator, the power variations were more significant in relation to the P1 UG generator when applied to the three techniques. The same occurred for the other P5 and P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA, and DE) compared to the Khan [8] technique, applied at 11:00 a.m. to meet the demand of 1088 MW, we were able to obtain a reduction of approximately 6.83%, which corresponds to $2691.41 in fuel cost with the use of PV solar energy. With this optimization, it was possible to reduce pollutants in the atmosphere (ALO 33.36%, DA 33.30%, and DE 33.40%) and the average reduction in pollutants in the atmosphere was approximately 33.35%. Figures 6-9 show the results of the data obtained in Table 6, which represent the fuel cost values generated by the techniques (ALO, DA and DE) at 11:00 a.m. demand.         Table 7 presents the comparative simulation results at 12:00, reaching a total power demand of 1240 MW.   Table 7 presents the comparative simulation results at 12:00, reaching a total power demand of 1240 MW. When evaluating Table 7, it is noted that the demand of 1240 MW necessary for the schedule is reached in the three optimization techniques, and the algorithm intelligently verifies the power of the six generators by turning off the less powerful ones. With the application of this technique, we were able to efficiently meet the demand, guaranteeing all the operating parameters of the system and reducing the fuel consumption and pollutants in the atmosphere.
We verified that the generator turned off in the simulations was the P2 UG. Even so, the energy generation did not suffer great variations, providing the stability of the energy generation to the system, in addition to allowing its predictive maintenance. On the other hand, the P1 UG generator presented lower power in the simulations, but in the DE technique, it exhibited better performance in relation to ALO and DA. In the P3 UG generator, the power variations were more significant in relation to the P1 UG generator when applied to the three techniques. The same was observed for the other P4, P5, and P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA, and DE), compared to the Khan [8] technique, applied at 12:00 to meet the demand of 1240 MW, we obtained a reduction of approximately 6.82%, which corresponds to $3188.84 in the cost of fuel, with the use of PV solar energy, it was possible to reduce pollutants in the atmosphere (ALO 30.52%, DA 30.52%, and DE 30.58%), the reduction in the average number of pollutants in the atmosphere was approximately 30.54%. Figures 10-13 show the results of the data obtained in Table 7, which represent the values of power, emissions and fuel costs generated by the techniques (ALO, DA and DE) at 12:00 demand. Obtaining the average fuel cost of the techniques (ALO, DA, and DE), compared to the Khan [8] technique, applied at 12:00 to meet the demand of 1240 MW, we obtained a reduction of approximately 6.82%, which corresponds to $3188.84 in the cost of fuel, with the use of PV solar energy, it was possible to reduce pollutants in the atmosphere (ALO 30.52%, DA 30.52%, and DE 30.58%), the reduction in the average number of pollutants in the atmosphere was approximately 30.54%. Figures 10-13 show the results of the data obtained in Table 7, which represent the values of power, emissions and fuel costs generated by the techniques (ALO, DA and DE) at 12:00 demand.             When evaluating Table 8, it was observed that the demand of 1135 MW necessary for the schedule was reached in the three optimization techniques, and the algorithm intelligently checked the power of the six generators by turning off the less powerful ones. With the application of this technique, we were able to efficiently meet the demand, guaranteeing all the operating parameters of the system and reducing the fuel consumption and pollutants in the atmosphere.
We verified that the generators turned off in the simulations were P2 and P3 UG; even so, the energy generation did not suffer great variations, providing the stability of the energy generation to the system, in addition to enabling predictive maintenance. On the other hand, the P1 UG generator presented lower power in the simulations; however, in the ALO technique, it performed better than the other techniques. In the P4 UG generator, the power variations were more significant in relation to the P1 UG generator when applied to the three techniques. The same occurred for the other P5 and P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA and DE), compared to the Khan [8] technique, applied at 13:00 to meet the demand of 1135 MW, we obtained a reduction of approximately 8.30% that corresponds to $3662.63 in the cost of fuel, with the use of photovoltaic solar energy, it was possible to reduce pollutants in the atmosphere (ALO 29.81%, DA 29.60% and DE 29.82%), the average reduction in pollutants in the atmosphere was approximately 29.75%. Figures 14-17 show the results of the data obtained in Table 8, which represent the fuel cost values generated by the techniques (ALO, DA and DE) at 13:00 p.m. demand. When evaluating Table 8, it was observed that the demand of 1135 MW necessary for the schedule was reached in the three optimization techniques, and the algorithm intelligently checked the power of the six generators by turning off the less powerful ones. With the application of this technique, we were able to efficiently meet the demand, guaranteeing all the operating parameters of the system and reducing the fuel consumption and pollutants in the atmosphere.
We verified that the generators turned off in the simulations were P2 and P3 UG; even so, the energy generation did not suffer great variations, providing the stability of the energy generation to the system, in addition to enabling predictive maintenance. On the other hand, the P1 UG generator presented lower power in the simulations; however, in the ALO technique, it performed better than the other techniques. In the P4 UG generator, the power variations were more significant in relation to the P1 UG generator when applied to the three techniques. The same occurred for the other P5 and P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA and DE), compared to the Khan [8] technique, applied at 13:00 to meet the demand of 1135 MW, we obtained a reduction of approximately 8.30% that corresponds to $3662.63 in the cost of fuel, with the use of photovoltaic solar energy, it was possible to reduce pollutants in the atmosphere (ALO 29.81%, DA 29.60% and DE 29.82%), the average reduction in pollutants in the atmosphere was approximately 29.75%. Figures 14-17 show the results of the data obtained in Table 8, which represent the fuel cost values generated by the techniques (ALO, DA and DE) at 13:00 p.m. demand.       Table 9 presents the comparative results of the simulation at 14:00, reaching a total demand of 1318 MW of power.   Table 9 presents the comparative results of the simulation at 14:00, reaching a total demand of 1318 MW of power. When evaluating Table 9, it is noted that the demand of 1318 MW required for the schedule was reached in the three optimization techniques, and the algorithm intelligently checks the power of the six generators and is unable to turn off the generators owing to the required demand. With the application of this technique, we were able to efficiently meet the demand, guaranteeing all the operating parameters of the system and reducing the fuel consumption and pollutants in the atmosphere.
We verified that the P1 and P2 UG generators presented lower power in the simulations, but the ALO technique performed better than the other techniques. In the P3 UG generator, the power variations were more significant in relation to the P1 and P2 UG generators when applied to the three techniques. The same occurred for the other P4, P5, and P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA and DE), compared to the Khan [8] technique, applied at 14:00 to meet the demand of 1318 MW, we obtained a reduction of about 3.51% that corresponds to $1930.38 in the cost of fuel, with the use of PV solar energy, it was possible to reduce pollutants in the atmosphere (ALO 21.47%, DA 21.57% and DE 21.53%), the average reduction in pollutants in the atmosphere was approximately 21.52%. Figures 18-21 show the results of the data obtained in Table 5, which represent the values of power, emissions and fuel costs generated by the techniques (ALO, DA and DE) at 14:00 p.m. demand. DA 21.57% and DE 21.53%), the average reduction in pollutants in the atmo approximately 21.52%. Figures 18-21 show the results of the data obtained in Table 5, which re values of power, emissions and fuel costs generated by the techniques (AL DE) at 14:00 p.m. demand.   DA 21.57% and DE 21.53%), the average reduction in pollutants in the atmo approximately 21.52%. Figures 18-21 show the results of the data obtained in Table 5, which re values of power, emissions and fuel costs generated by the techniques (AL DE) at 14:00 p.m. demand.    The Table 10 presents the comparative results of the simulation at 15:00, reaching a total demand of 1074 MW of power.   The Table 10 presents the comparative results of the simulation at 15:00, reaching a total demand of 1074 MW of power.  The Table 10 presents the comparative results of the simulation at 15:00, reaching a total demand of 1074 MW of power.
From Table 10, it is observed that the demand of 1074 MW necessary for the schedule is reached in the three optimization techniques, and the algorithm intelligently checks the power of the six generators, regardless of whether it is possible to turn off the less powerful ones. With the application of this technique, we were able to efficiently meet the demand, guaranteeing all the operating parameters of the system and reducing the fuel consumption and pollutants in the atmosphere.
We verified that the generator turned off in the simulations was the P2 UG. Even so, the energy generation did not suffer great variations, providing the stability of the energy generation to the system, in addition to allowing its predictive maintenance. On the other hand, the P1 UG generator presented lower power in the simulations, but in the DA technique, it showed better performance in relation to ALO and DE. In the P3 UG generator, the power variations were more significant in relation to the P1 UG generator when applied to the three techniques. The same occurred for the other P4, P5, and P6 UG generators.  Table 10, it is observed that the demand of 1074 MW necessary ule is reached in the three optimization techniques, and the algorithm checks the power of the six generators, regardless of whether it is possible less powerful ones. With the application of this technique, we were able meet the demand, guaranteeing all the operating parameters of the system the fuel consumption and pollutants in the atmosphere.
We verified that the generator turned off in the simulations was the so, the energy generation did not suffer great variations, providing the s energy generation to the system, in addition to allowing its predictive ma the other hand, the P1 UG generator presented lower power in the simu the DA technique, it showed better performance in relation to ALO and UG generator, the power variations were more significant in relation to th erator when applied to the three techniques. The same occurred for the oth P6 UG generators.
Obtaining the average fuel cost of the techniques (ALO, DA, and DE) the Khan [8] technique, applied at 15:00 to meet the demand of 1074 MW, reduction of approximately 5.75%, which corresponds to $2591.25 in th With the use of PV solar energy, it was possible to reduce pollutants in th (ALO 22.65%, DA 22.45%, and DE 22.70%), the average reduction in po atmosphere was approximately 22.60%. Figures 22-25 show the results of the data obtained in Table 5, which values of power, emissions and fuel costs generated by the techniques ( DE) at 15:00 p.m. demand.        It can be seen in the graphs of the Pareto Front figures (Cost × Emissions), for example, Figure 25, that the curves are inversely proportional, therefore, a multi-objective problem is configured. Table 11 presents the simulated metaheuristic techniques from 10:00 a.m. to 15:00 p.m. and the results achieved by each technique. We observed that the DE technique achieved the best fuel cost reduction five times in the six simulations, DA obtained the second-best result, ALO was in third place, and obtained better results only on PSO. The proposed algorithms maintained the generators at their optimal power to obtain better efficiency at all times. The new proposal shows that among the techniques used, the DE technique was the best, guaranteeing the reduction in fossil fuel use by 6.67%, corresponding to $18,894.31, DA reduction by 6.09%, corresponding to $17,228.66, ALO obtained the smallest reduction 5.81%, corresponding to $16,458.44, the results were compared with the PSO technique.
During the simulations, the greatest reduction occurred at 13:00 p.m. with 8.64% and the smallest reduction occurred at 14:00 p.m. with 4.34%. Figure 26 shows the results of the data obtained in Table 11, in which the general fuel cost obtained between 10:00 a.m. and 15:00 p.m., achieves the best performance (cost reduction) with the DE technique. It can be seen in the graphs of the Pareto Front figures (Cost × Emissions), for example, Figure 25, that the curves are inversely proportional, therefore, a multi-objective problem is configured. Table 11 presents the simulated metaheuristic techniques from 10:00 a.m. to 15:00 p.m. and the results achieved by each technique. We observed that the DE technique achieved the best fuel cost reduction five times in the six simulations, DA obtained the second-best result, ALO was in third place, and obtained better results only on PSO. The proposed algorithms maintained the generators at their optimal power to obtain better efficiency at all times. The new proposal shows that among the techniques used, the DE technique was the best, guaranteeing the reduction in fossil fuel use by 6.67%, corresponding to $18,894.31, DA reduction by 6.09%, corresponding to $17,228.66, ALO obtained the smallest reduction 5.81%, corresponding to $16,458.44, the results were compared with the PSO technique.
During the simulations, the greatest reduction occurred at 13:00 p.m. with 8.64% and the smallest reduction occurred at 14:00 p.m. with 4.34%. Figure 26 shows the results of the data obtained in Table 11, in which the general fuel cost obtained between 10:00 a.m. and 15:00 p.m., achieves the best performance (cost reduction) with the DE technique.

Conclusions
Aiming to reliably show alternatives to supply the generation of energy has been one of the research objectives for several years. Thus, this work is a viable solution for sustainable clean energy, responsible for solving the problem of optimization of CEED using a hybrid system composed of TPP and photovoltaic generation. By choosing photovoltaic plants to meet the energy demand, a cleaner production with less environmental impact was achieved, contributing to a decrease in contaminants released into the atmosphere. The new method intelligently checks the capacity of the generators in use and, depending on the demand requested for the time, switches off the least efficient thermal generator. With the application of these techniques, it is possible to meet the energy demand efficiently, guarantee all the operating parameters of the system, reduce fuel costs, and reduce pollutant emissions. Thus, the application of metaheuristics to optimize CEED proved to be efficient and reliable with excellent results. The new proposal shows that among the techniques used, DE was the method that presented the best results, guaranteeing a reduction in the use of fossil fuel by 6.67%, corresponding to $18,894.31. The DA also obtained a reduction of 6.09% corresponding to $17,228.66, while ALO obtained 5.81% corresponding to $16,458.44. These results were compared with the results of Khan simulations, which used PSO. The results achieved using the ALO, DA, and DE techniques demonstrate the robustness of the application of these algorithms in optimization problems.

Conclusions
Aiming to reliably show alternatives to supply the generation of energy has been one of the research objectives for several years. Thus, this work is a viable solution for sustainable clean energy, responsible for solving the problem of optimization of CEED using a hybrid system composed of TPP and photovoltaic generation. By choosing photovoltaic plants to meet the energy demand, a cleaner production with less environmental impact was achieved, contributing to a decrease in contaminants released into the atmosphere. The new method intelligently checks the capacity of the generators in use and, depending on the demand requested for the time, switches off the least efficient thermal generator. With the application of these techniques, it is possible to meet the energy demand efficiently, guarantee all the operating parameters of the system, reduce fuel costs, and reduce pollutant emissions. Thus, the application of metaheuristics to optimize CEED proved to be efficient and reliable with excellent results. The new proposal shows that among the techniques used, DE was the method that presented the best results, guaranteeing a reduction in the use of fossil fuel by 6.67%, corresponding to $18,894.31. The DA also obtained a reduction of 6.09% corresponding to $17,228.66, while ALO obtained 5.81% corresponding to $16,458.44. These results were compared with the results of Khan simulations, which used PSO. The results achieved using the ALO, DA, and DE techniques demonstrate the robustness of the application of these algorithms in optimization problems.