Next Article in Journal
Critical Raw Materials and Supply Chain Disruption in the Energy Transition
Previous Article in Journal
A Review of Using IoT for Energy Efficient Buildings and Cities: A Built Environment Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Term Hydro-Thermal-Solar Scheduling with CCGT Based on Self-Adaptive Genetic Algorithm

1
Faculty of Electrical Engineering and Information Technologies, Ss. Cyril and Methodius University, Ruger Boshkovikj 18, P.O. Box 574 Skopje, North Macedonia
2
Faculty of Electronic Engineering and Technologies, Technical University of Sofia, 8 Kl. Ohridski Blvd, 1000 Sofia, Bulgaria
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(16), 5989; https://doi.org/10.3390/en15165989
Submission received: 20 July 2022 / Revised: 12 August 2022 / Accepted: 14 August 2022 / Published: 18 August 2022
(This article belongs to the Section B2: Clean Energy)

Abstract

:
This paper presents a new metaheuristic approach based on a self-adaptive genetic algorithm (SAGA) for solving the short-term hydro-thermal-solar scheduling with combined-cycle (CCGT) units. First of all, the proposed approach is applied to a test system with different characteristics, considering the valve-point effect. The simulation results obtained from the new SAGA are compared with the results obtained from some other metaheuristic methods, such as AIS, DE, and EP to reveal the validity and verify the feasibility of the proposed approach. The test results show that the proposed metaheuristic approach proves the effectiveness and superiority of the SAGA algorithm for solving the short-term hydro-thermal-solar scheduling (SHTSS) problem.

1. Introduction

In the recent period, the electric power demand has increased and fossil fuel prices have risen, which has led to growing the world energy crisis. Therefore, the world has tended to reduce the use of TPP by using renewable energy sources to reduce the emissions and harmful gases that are released from TPP. These harmful gases harm the environment, and this leads to an increase in the temperature of the planet, which causes global warming. Accordingly, the European Union has introduced stricter directives, which provide for the decommitment of TPP, and a gradual increase in the renewable penetrations in the system [1,2,3].
Therefore, the incorporation of the short-term hydrothermal generation scheduling (STHS) problem and new renewable energy integration has great importance in the power system operation [4]. The STHS problem is one of the most important optimization problems in power system planning. The primary goal of the STHS problem is to determine the optimal power-generation schedule of the TPP and HPP to minimize the total operation cost of the system. On other hand, the total operating costs and flexibility of the TPP’s operation can be favorably influenced by committing combined cycle gas turbines (CCGT) into the system and the appropriate decommitment of existing conventional TPPs, especially in the transition period, which has already begun. However, in SHTSS, one of the basic features of solar power plants (SPP) is the extreme variability and unpredictability of the output power, primarily due to the dependence on weather conditions [5].
These two attributes are especially evident in modern power systems, in which solar power plants have an increasing penetration. In addition to this, modern PV modules use bifacial technology [5], so their orientation significantly affects the output power of the solar power plant and, thus, all the parameters of the system (including costs). According to what has been said so far, it can be said that this paper deals with the impact of the orientation of the modules on the total operating costs. This analysis is of particular importance for the integration of SPP in the power system.

2. Motivation, Literature Review, and Contributions

To solve the SHTS problem, various approaches are applied based on the achievements of applied mathematics and mathematical software, and, initially, various classical optimization methods were proposed, such as the dynamic programing method [6,7], the Lagrange relaxation method [8,9], and mixed integer programing [10,11]. In this aspect, it is known that the considered optimization methods are difficult to apply in solving complex optimization problems. On the other hand, gradient methods are characterized by fast convergence but are inefficient and of limited applicability to problems that are described by nonconvex and discontinuous objective functions, such as the SHTSS problem.
In articles [12,13,14,15,16,17], techniques based on application of artificial intelligence are considered: evolutionary programing, simulated annealing, differential evolution, artificial neural network, genetic algorithm (GA), and particle swarm optimization (PSO) for solving of the SHTS optimization problem. A characteristic of these manuscripts is that only a quadratic objective function is used and, also, significantly simpler constraints. Thus, the SHTS optimization problem has become more complex by incorporating the valve point load and transmission loss into the objective function while considering various hydro, thermal, system, and security constraints.
In this regard, it is known that, when applying metaheuristic approaches to solve optimization problems such as SHTS, there is a significantly greater possibility of avoiding local minima and obtaining an optimal global approach. Due to the widespread distribution of decentralized power systems, the hydrothermal scheduling has been widely studied by researchers worldwide.
In previous studies, some of the considered algorithms failed to obtain optimal solutions when considering various complex constraints based on multi-objective tasks. An excellent real-coded genetic algorithm (RCGA) solution was obtained by introducing a self-adaptive penalty function, Laplace crossover, and Makinen, Perriot, and Toivanen (MPTM) mutation. Thus, the planning problem is more accurately described based on the formulation of constraints inherent in the operation of a real system, such as reservoir storage capacity constraints, water discharge constraints, available production constraints, generator constraints, ramp rate constraint, valve point effect, and transmission line constraint.
All optimization methods strive to obtain the global optimum, but gradient methods easily become stuck in local optima. GA has proven to be a better method of obtaining the global optimum because it works with a population which performs diversification and intensification of the search space, i.e., a group of solutions, compared to gradient methods, which work with a single solution.
In this article, a new constraint handling technique is proposed through the fitness function evaluation by a new self-adaptive penalty, Laplace crossover, and adaptive crossover and mutation strategy. A newly proposed constraint handling repair mechanism has been proposed for consideration of the constraints that are most difficult to satisfy, especially the power balance constraint and the hydro constraints. The main contributions of this paper are:
  • A new self-adaptive penalty for constraints handling, which require no tuning.
  • A new crossover technique, i.e., Laplace crossover, which has a self-adaptive tuning ability, which is important to maintain population diversity.
  • A new constraint handling repair mechanism for simultaneous satisfaction of all constraints, especially power balance and hydraulic continuity equation, which are neglected in other papers. This allows for a significantly more physically realistic solution.
  • Analysis of the influence of the orientation of the modules on the output power of the solar power plant and, thus, on the overall system parameters. This is especially important for system operators, as they receive future scenarios for operational planning. Therefore, they will be able to decide on the maximum installed power of solar power plants in the system.
  • The paper is of particular importance to the academic community, as it presents a scenario in operational planning in a “green energy transition”, with a slow departure from coal, i.e., in a period of committing CCGTs, decommissioning TPP, and solar power plants taking an increasing penetration.
  • An additional advantage of the proposed approach is that it can be used to create a graphical user interface (GUI) that can be developed and enhanced. Having a GUI is a very important tool for system operators in operational planning and real-time operational decision making.

3. Problem Formulation

3.1. Combined Cycle Gas Turbines

TPP with combined cycle gas turbines (CCGT), shown in Figure 1, has a gas and steam turbine. The main purpose of such plants is to use the heat generated at the outlet of the gas turbine. Since the exhaust gases coming out of the gas turbine have extremely high temperatures, about 600 °C can be used as a means to heat water and produce steam for the steam turbine. This increases the efficiency of the process itself, as the heat that would otherwise be lost is used for further steam production. The efficiency of such a plant reaches up to 60%. In a CCGT, the compressor compresses the air and sends it to the combustion chamber, where the combustion fuel is supplied at the same time.
Very high-temperature combustion gases are led from the combustion chamber to the gas turbine, where they expand, giving useful work on the shaft connected to the gas turbine rotor. The shaft drives the generator and produces electricity that is sent to the grid.
After expansion, the exhaust gases from the gas turbine are led to the heat recovery steam generator (HRSG). There is still a lot of unused air at the output of the gas turbine, and this excess air is used to burn additional fuel in the HRSG. In the HRSG, the feed water is heated to evaporation and overheats to the set parameters. The superheated steam goes from the HRSG to the steam turbine, where it expands and hands over the mechanical work to the electric generator. After that, the steam, now of low parameters, goes to the condenser, where it condenses. After condensation, the water is returned to the HRSG by the feed pump for reheating. It has already been mentioned that, with this principle, we increase the usability of the whole process, by changing CCGT’s operation mode, which is shown in Figure 2 [17], i.e., the optimal output power, according to the daily load diagram [18,19].

3.2. Objective Function

The main objective in solving the SHTSS optimization problem is to minimize fuel costs in the thermal power plants (TPPs). In this regard, fuel costs in HPPs are neglected, as the ongoing costs of hydropower generation systems are significantly lower compared to other production costs, such as thermal energy. In this way, the optimal fuel costs are determined for the electricity generation from the TPPS. When using an optimization procedure in modern power systems, the main objective is to minimize the total fuel cost by properly using renewable energy sources, such as SPPs, while providing the demand power for the whole optimization period, which, in SHTSS, is one day, i.e., 24 h. It should be noted that, firstly, the units are committed according to [18,19,20] (considering various constraints), and then the SHTSS problem is solved.
The fuel costs of TPPs are usually represented by a quadratic function and depend on the generator output power. Therefore, the objective function of the problem is the sum of the fuel costs of all TPPs that are subject to optimization, i.e., [21,22,23]:
min F = j = 1 J t = 1 N T [ F t , j j ] u t , j t N T ; j J ; u 0 , 1
where u represents the commitment state of TPP t at interval j (0 for a decommitted unit or 1 for a committed unit). On other hand, Ft,j is the fuel cost function of the TPP t and is represented by the nonconvex function:
F t , j = a t + b t P GT t , j + c t P GT t , j 2 + | d t sin ( e t ( P GT t min P GT t , j ) ) | t N T ; j J
where at, bt, ct, dt, and et represent cost coefficients of the TPP t, and P GT t min is the technical minimum, while j is the index of the optimization interval.

3.3. Constraints

3.3.1. Generator Constraint

The generation output power of TPPs and HPPs is enclosed in the technical minimum and the technical maximum, and is represented by:
u t , j P GT t min P GT t , j u t , j P GT t max u h , j P GH h min P GH h , j u h , j P GH h max

3.3.2. Power Balance Constraint

The total output power from all TPP (including CCGT), HPP, and SPP must be equal to the total load and losses in the system:
t = 1 N T u t , j P GT t , j + h = 1 N H u h , j P GH h , j + s = 1 N S u s , j P GS s , j = P P , j + P L , j
P L , j = i = 1 N G j = 1 N G P G i B i j P G j + j = 1 N G B i 0 P G j + B 00 i N G ; j N G ; N G = N T + N H

3.3.3. Spinning Reserve Constraint

The required level of spinning reserve, which is required for the system security, is given as:
i = 1 N T + N H u i , j P G i max P P , j + P L , j + R ;   j J ;   u 0 , 1
while the required spinning reserve is calculated according to the empirical formula of ENTSO (UCTE), i.e., [24]:
R = a res P P , max + b res 2 b res ;   a res = 10   MW ;   b res = 150   MW
where ares and bres represent empirical constants, and PP,max the system’s peak load.
Thus, the spinning reserve constraint with available spinning reserve from TPPs and HPPs is represented by:
t = 1 N T u t , j P GT t AV R T h = 1 N H u h , j P GH h AV R H
where P GT , t AV and P GH , h AV are the available spinning reserve of TPP t or HPP h; R T is the total required spinning reserve of TPPs; and R H is the totally required spinning reserve of HPPs. It is considered that R T = 0.75 R and R H = 0.25 R .

3.3.4. Ramp Rate Constraint

The sudden change in the output power of the generators is limited by the ramp rate constraint, i.e.,
P GT t , j = max ( P GT t min , ( P GT t , j D R T t ) ) P GH h , j = max ( P GH h min , ( P GH h , j D R H h ) )
P GT t , j = min ( P GT t max , ( P GT t , j + U R T t ) ) P GH h , j = min ( P GH h max , ( P GH h , j + U R H h ) )
where URTt, DRTt, URHh, and DRHh represent the allowable upper and down rates of TPP t or HPP h.

3.3.5. Transmission Line Constraint

The maximum power that can be transmitted by a transmission line is represented by the following constraint:
P line , g P line , g max ,   g = 1 , , G
where G is the total number of transmission lines in the system. The active power of the transmission line is obtained by power flow calculation by applying the DC model, i.e., DC power flow [20].

3.3.6. Water Availability Constraint

The total water discharge is physically constrained by a total available volume i.e.,
j = 1 J Q t h , j T j V h , k
where Tj is the duration of interval j, Qth is the water discharge, i.e., the input–output curve of the HPP, and is represented by a quadratic function:
Q t h P GH , h = α h + β h P GH h + γ h P GH h 2
where PGHh is the output power of HPP h, and αh, βh, and γh are constant coefficients of the input–output curve.

3.3.7. Available Production Constraint

The available production of TPPs is defined according to the following expression:
i = 1 N T + N H P G i , j T j = W max , i
The available production of HPPs is defined according to the available (initial) volume V k and the total discharge time Tdis, i.e.,
Q max , h = Q ins , h = f ( α h , β h , γ h , P GH , h max )   ( m 3 / h )
T dis , h = V k , h Q max , h   h
W max , h = P GH h max T dis , h   MWh

3.3.8. Dynamic Balance of the Reservoir Storage

Reservoir storage of HPP is compressed by spillage and inflow at the preceding event and it must track the continuity equations of the hydraulic system at each time interval j, and it is described as follows:
V h , j = V h , ( j 1 ) + I h , j Q t h , j S h , j
where V h , j is the storage volume of HPP h; I h , j is the inflow in reservoir h; and S h , j is the water spillage of the reservoir h. In this paper, water spillage is neglected.

3.3.9. Initial and Final Reservoir Storage Constraint

The volume constraints of the reservoir at the beginning and end of the scheduling period are determined by:
V h , 0 = V h initial V h , 24 = V h end

3.3.10. Water Discharge Constraint

The water discharge constraint is determined by:
Q h min Q h , j Q h max

3.3.11. Reservoir Volume Constraint

The physical limitation of reservoir volume is given by:
V h min V h , j V h max

4. Genetic Algorithm

4.1. Initialization

The initial population is modeled based on the optimal solution, previously obtained from unit commitment (as conducted in reference [20]), as follows:
P G 0 , m a = 1 ψ P GUC , m ; P G 0 , m b = 1 + ψ P GUC , m P n , m = P G 0 , m a + rand P G 0 , m b P G 0 , m a n = 1 , , N pop ; m = 1 , , N var = N T + N H J
where ψ is the diversity factor and has a value of 0.2.

4.2. Fitness Function Evaluation and Constraint Handling

The initial population contains 200 chromosomes. To increase the robustness of the algorithm but also to maintain the quality of the initial population, before the main stage of the algorithm, the best 100 chromosomes are selected. This selection is made according to the following expression [25,26]:
F x = f X if   X   is   feasible f max + i = 1 I g i X + k = 1 K | h k X | if   X   is   infeasible
where fmax is the value of the objective function of the worst feasible solution, i.e., chromosome in the population.
In this paper, a new constraints handling technique is proposed for fitness function evaluation by a new self-adaptive penalty function, which has not been applied to this type of optimization problem so far. It consists of distance function d(X) and penalty function p(X). The proposed technique works only with inequality constraints, i.e., [27,28,29,30]:
G i X = max g i X , 0 , i = 1 , , I max h k X ε , 0 , k = I + 1 , , K
and, after that, the final fitness function receives a self-adaptive value (but not the sum of the objective function and classical penalty function):
F X = d X + p X
d X = v ( X ) ,   if   r f = 0 f X 2 + v ( X ) 2 ,   otherwise
v X = i = 1 m w i ( G i ( X ) ) i = 1 m w i
f X = f X f min f max f min
p X = 1 r f M X + r f N X
M X = 0 ,   if   r f = 0 v X ,   otherwise
N X = 0 ,   if   X   is   a   feasible f X ,   if   X   is   an   infeasible
where, wi = 1/Gmax,i is a weight parameter; Gmax,i is the maximum value for violation of constraint Gi(X) recorded thus far; and v(X) is the total constraint violation. The main advantage of the proposed algorithm is that, unlike the static penalty, the self-adaptive penalty does not require parameter setting.

4.3. Selection

To prevent premature convergence, the fitness function is scaled linearly. The relationship between the original fitness function and the scaled one is given by the expression:
f s = a s f + b s a s = s p 1 f av / f max f min b s = 1 a s f av
where sp is a parameter that controls the selection pressure and is in the interval (1.2–2), fs is the scaled fitness value, F is the fitness value, fav is the mean fitness of the population, fmax and fmin are the maximum and minimum fitness value in the current population, and as and bs are scaling coefficients [31,32,33].

4.4. Crossover

For the proposed SAGA, Laplace crossover (LX) was applied. Two offspring are generated from a pair of parents, as described. First, a uniformly distributed random number αL is generated, so α L 0 , 1 . Then, a random number βL is calculated, which follows the Laplace distribution by simply inverting the distribution function of the Laplace distribution, i.e., [34]:
β L = a L b L ln α L ,   α L 0.5 a L + b L ln α L ,   α L > 0.5
The offspring are given by the equations:
y i 1 = x i 1 + β L x i 1 x i 2 y i 2 = x i 2 + β L x i 1 x i 2

4.5. Mutation

In the newly proposed self-adaptive genetic algorithm (SAGA), the Makinen, Periaux, and Toivanen mutation (MPTM) is applied (which is originally proposed by Makinen [35] and is applied for a general optimization problem). In this paper, MPTM is applied for the first time to that optimization problem. The implementation of MPTM in SAGA consists of the following. First, an evenly distributed random number r is generated, such that r 0 , 1 . Then, the new mutated chromosome x ^ is calculated as follows:
x ^ = 1 p ^ x l + p ^ x u
p ^ = p p p r p b if   r < p p if   r = p p + 1 p r p 1 p b if   r > p
p = x x l x u x
where x l and x u are the upper bound and the lower bound, respectively.

4.6. Adaptive Crossover and Mutation Strategy

The target of the newly proposed crossover and mutation strategies is to regulate the crossover and mutation probability of the chromosome, and thus maintain population diversity. These strategies are based on the difference between its fitness function and the average fitness function of the population [20], i.e.,
P c = k 1 f max f c / f max f av , f c f av k 2 , f c < f av
P m = k 3 f max f c / f max f av , f c f av k 4 , f c < f av
where k1 < 1, k2 < 1, k3 < 1, and k4 < 1. The values of the k parameters, i.e., k1, k2, k3, and k4, were obtained by the “trial and error” method. Because these parameters regulate the crossover and mutation probability, a large number of simulations were performed on a suitable test system, and then results are compared with the results of other proposed algorithms. Taking into account that these parameters significantly affect the performance of the algorithm and the CPU time, an optimal interval is determined, which represents a balance between the quality of the results (the global optimum) and the CPU time, so it is recommended 0.8 ≤ k1 ≤ 1; 0.17 ≤ k2 ≤ 0.22; and 0.94 ≤ k3 ≤ 1; 0.14 ≤ k4 ≤ 0.18. In this paper, k1 = 0.95, k2 = 0.20, and k3 = 1, k4 = 0.18.

4.7. Elitism Strategy

In this paper, an elitism strategy is proposed by creating a combined population that consists of the best 50 parents and the best 50 children. This combined population represents the population for the next generation. Figure 3 shows the simplified block diagram of the SAGA, while Algorithms 1–3 represent different parts of the newly proposed constraint handling repair mechanism.
Algorithm 1: New constraints handling for reservoir storage volume (initial and final reservoir storage)
  1:
Set V 1 end = v 10 4 ;   V 2 end = v 10 4 ;
  2:
Set ε h = 0.1 ;   i t e r max = 7 ;   j min = 1 ;   j max = N D ;   i t e r = 1 ;
  3:
while h N H
    a.
Δ V h = V h initial + j = 1 J I h , j j = 1 J Q h , j V h end
    b.
if Δ V h ε h   o r   i t e r > i t e r max then
  i.
j rand = r a n d i j min   j max , 1 , 1
  ii.
Q h , j rand = max Q h min , min Q h , j rand + Δ V h , Q h max
    c.
else
  i.
Q h , J = Q h , 1 Q h , 2 Q h , 24
  ii.
Q h , J = max Q h min , min Q h , J + Δ V h / N D , Q h max
  iii.
i t e r = i t e r + 1
    d.
end if
  4:
end while
Algorithm 2: New constraints handling for constraints handling for water discharge
  1:
for h = 1 : N H
   a.
for j = 1 : N D
    i.
if Q h , j < Q h min then
 1.
Q h , j = Q h min
    ii.
else if Q h , j > Q h max then
 1.
Q h , j = Q h max
    iii.
end if
   b.
end for
  2:
end for
Algorithm 3: New constraints handling for real power balance and ramp rate.
  1:
for i = 1 : N H
   a.
for j = 1 : J
    i.
if P GH i , j 0 then
    1.
if P GH i , j max P GH min i , P GH i , j D R H i then
    a.
P GH i , j = max P GH min i , P GH i , j D R H i
    2.
else if P GH i , j min P GH max i , P GH i , j + U R H i
    a.
P GH i , j = min P GH max i , P GH i , j + U R H i
    3.
end if
    ii.
end if
   b.
end if
  2:
end for
  3:
for j = 1 : J
   a.
W H new j = i = 1 N H P GH i , j
   b.
P PT j = P P j + P L j W H new j
  4:
end for
  5:
P L j = m = 1 N T + N H n = 1 N T + N H P G , m j B m n P G , n j + m = 1 N T + N H B m 0 P G m + B 00
  6:
for j = 1 : J
   a.
Δ P PT j = t = 1 N T P GT t , j P PT j P L j
   b.
while Δ P PT j > 10 3
    i.
Δ P PT , mean j = Δ P PT j / N T
    ii.
for t = 1 : N T
   1.
if P GT t , j 0 then
   a.
P GT t , j = P GT t , j Δ P PT , mean j
   b.
if P GT t , j max P GT min t , P GT t , j D R T t then
   i.
P GT t , j = max P GT min t , P GT t , j D R T t
   c.
else if P GT t , j min P GT max t , P GT t , j + U R T t then
   i.
P GT t , j = min P GT max t , P GT t , j + U R T t
   d.
end if
   2.
end if
   3.
Δ P PT j = t = 1 N T P GT t , j P PT j P L j
    iii.
end for
   c.
end while
  7:
end for

4.8. Computational Procedure

The computational procedure of the newly proposed SAGA is shown in Figure 3, and is briefly described as follows:
Step 1: Preparation and initialization. Determine the necessary computational parameters of the algorithm. Generate the initial population from the previous UC solution, perform fitness function evaluation using Equation (24), and select the top 100 chromosomes.
Step 2: Calculate the fitness function values and total constraint violation of each chromosome in the parent population using Equations (25)–(31). Then, determine the elite chromosome.
Step 3: Apply a linear rank selection. Each chromosome in the population is ranked in increasing order of fitness function, from 1 to n. Linear ranking assigns a selection probability to each chromosome that is proportional to the individual’s rank.
Step 4: LX and MPTM is used to generate the offspring population. To enhance the population diversity and prevent premature convergence, before the predefined crossover operator and mutation operator, an appropriate strategy will be applied. These strategies regulate the crossover and mutation probability, based on the fitness deviation of the chromosome from the mean fitness function of the entire population.
Step 5: Fitness function evaluation. Calculate the fitness function values and total constraint violation of each chromosome.
Step 6: Determine constraint violations, and then apply the new constraint handling repair mechanism, in order for the chromosome to be feasible again.
Step 7: Apply an elitism strategy. The parent population and offspring population are combined in one group. They are then sorted based on the fitness function in descending order. Then, the better, i.e., best 100 chromosomes, will be chosen as the members in the new generation.
Step 8: Repeat steps 3 to 7 until the maximum number of generations is reached, then export the optimal solution, i.e., chromosome.

5. Simulation Experiment and Result Analysis

5.1. Test System 1

The investigated system consists of two TPPs and two HPPs, and the effect of the valve point is considered in the simulation. The data for conducting the numerical experiments were taken from [36]. Table 1 shows the optimal generation of the studied system. In the performed simulations, the population size and maximum generations are 100 and 500, respectively, as shown in [37].
In the simulations, the spin reserve is chosen to be 10% of the total load demand. Table 2 presents a comparison of the results for the test system using different optimization methods, which was obtained in the authors’ previous research [37]. In the particular example, the experiment was conducted 50 times. Thus, the average fuel costs were obtained based on the data from conducting 50 experiments.

5.2. Test System 2

In this subchapter, the efficiency of the first part of the proposed algorithm, i.e., the unit commitment part, will be verified. For this purpose, the algorithm will be applied to a test system, which consists of six thermal power plants. All data, such as unit characteristics and load demand, are taken from [38,39]. On an Intel Core i7-9750H [email protected] with 16GB RAM, 50 independent simulations were made for a relevant comparison with other methods. Table 3 shows the total fuel costs for test system 2.

5.3. IEEE 30 Bus System

The performance of the proposed algorithm, after its verification, has been evaluated using the Modified IEEE 30 bus system. This system consists of 30 buses, 7 generators, i.e., four TPP (buses 1, 2, 5, and 8), of which TPP1 (in bus 1) is with CCGT, two HPP (buses 11 and 13), and one SPP (bus 10), and 41 transmission lines. The data for this system are taken from [38,40,41,42,43]. The data for inflows in reservoirs 1 and 2 are taken from [20]. The data for CCGT TPP are shown in Table 4 and Table 5 and are taken from [44]. A spinning reserve value of 10% of the system’s peak load is selected. When solving the optimization problem, 75% of the spinning reserve was chosen to be covered by TPPs, and the remaining 25% by HPPs. Figure 4 shows the daily load diagram of the system, and Figure 5 shows its single-line diagram. For an IEEE 30 bus system, the proposed approach runs for 50 independent simulations, and the population size and maximum number of iterations are set to 100 and 500, respectively.

5.3.1. Solar Power Plant Model

The SPP model was developed in the PVSYST software package. It implements 34,000 Trina Solar bifacial 500 Wp PV modules and 170 grid-tie inverters of the type Fronius International—Tauro Eco 100-3-P, with a modeled DC/AC ratio of 1.00. The SPP consists of 2125 strings and has an installed power of 17 MWp. The albedo of the substrate is 0.7, and the height above ground is 1.5 m. The data for PV modules are taken from [45].
To prove the initial hypothesis, i.e., that the orientation and location of the SPP affect the total fuel costs, two SPP (i.e., cases) are modeled. In case 1, the location of the SPP is Gaag, Netherlands (latitude: 51.96° and longitude: 4.29), while, in case 2, the location of the SPP is Bitola, Macedonia (latitude: 41.03° and longitude: 21.33). Given the HUPX’s electricity price diagram, which is shown in Figure 6 [46], each of these two cases involves three different sub-cases. In other words, the electricity price diagram has two peaks, one in the morning and the other in the afternoon. To achieve the greatest possible energy balance (largest production in the morning or afternoon), in two of the three sub-cases, the orientation of the modules, i.e., the azimuth, will be 90° and −90°.
The parameters for both cases are shown in Table 6 and Table 7. The best values of parameters, such as height above ground sheds spacing, injected energy per year (EGRID), and tilt angle, are obtained through “trial and error”.
Figure 7 shows the output power of the SPP for case 1, and Figure 8 shows the output power for case 2 for the specific optimization period, i.e., 9 June 2021.
From Table 6 and Table 7, it can be seen that the annual production for case 2 is significantly higher than the production for case 1. This should not come as a surprise, given the fact that the global horizontal radiation for case 2 is 1480.3 kWh/m2, while, for case 1, it is 1040.0 kWh/m2, meaning case 2 is 1.42 times higher compared to case 1.

5.3.2. Case 1: Gaag, Netherlands

Obtained heat rate (HR) for all TPPs, including the CCGT, is given in Table 8. The average value of total fuel cost for all sub-cases, obtained using the proposed method is shown in Table 9, and the computation time taken by the algorithm is 21.28 s. The average fuel costs are obtained from 50 independent simulations. Figure 9, Figure 10, and Figure 11, respectively, show the optimal output power of TPPs and HPPs for sub-case 1, 2, and 3, while the optimal generation scheduling for the optimization period for the lowest cost sub-case, using the proposed optimization algorithm, is shown in Figure 10 and in Table 10.

5.3.3. Case 2: Bitola, Macedonia

The average fuel cost for all sub-cases, obtained using the proposed method, is shown in Table 11, and the computation time taken by the algorithm is 22.14 s. The average fuel costs are obtained from the 50 independent simulations. Figure 12, Figure 13, and Figure 14, respectively, show the optimal output power of TPPs and HPPs for sub-case 1, 2, and 3, while the optimal generation scheduling for the optimization period for the lowest cost sub-case, obtained by newly proposed SAGA, is shown in Figure 14 and in Table 12.

5.4. Analysis of the Obtained Results

First of all, the total fuel costs of the test system 1 are EUR 57,640, compared to other previously proposed algorithms, where they are EUR 66,117, EUR 66,121, and EUR 66,198. The same was established for test system 2, which refers to unit commitment, where the total costs amount to EUR 13,126.31, compared to the previously proposed methods, i.e., EUR 13,292.28 and EUR 13,600.00. According to this, it can be concluded that the proposed SAGA provides better results (for all parts, i.e., unit commitment and hydro-thermal scheduling) compared to optimization algorithms proposed in [30,38,39].
For case 1, the sub-case with the lowest total fuels costs is sub-case 2, where they are EUR 10,255.34. For sub-cases 1 and 3, the total fuel costs are EUR 10,316.19 and EUR 10,320.99, which is a relative change, i.e., sensitivity by 0.59% and 0.63%, compared to sub-case 2, which is not a small amount, in terms of short-term planning.
For case 2, the best scenario is sub-case 3, where the total fuel costs are EUR 10,184.02. For the sub-cases 1 and 2, the total fuel costs are EUR 10,397.02 and EUR 10,345.02, which is a relative change, i.e., sensitivity by 2.05% and 1.55%, respectively, compared to sub-case 3, which is also not negligible. The different sensitivities for cases 1 and 2 are based on the fact that these two different locations have significantly different global horizontal solar radiation. The annual global horizontal solar radiation for Gaag, Netherlands is 1040.0 kWh/m2, while, for Bitola, Macedonia, it is 1480.3 kWh/m2, i.e., even 50 percent higher.
This statement can confirm the initial hypothesis, i.e., that the different orientation of the PV modules affects the total fuel costs in the system. It is important that the different orientation of the PV modules also affects other system parameters, such as active power of transmission lines and optimal output power of TPPs and HPPs, according to Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. On the other hand, the auxiliary hypothesis can also be confirmed, i.e., that the different location of the SPP installation also has an impact on the total fuel costs, even for the same topological structure of the system, which is the case in this paper. The confirmation of the hypothesis is also based on the fact that the transmission line constraint has a major impact on fuel costs because it directly affects the optimal power flow of the system. This means that if the active power of one of the transmission lines in a given time interval is equal to its maximum capacity, then the optimal output power from an HPP or TPP (regardless of whether it is the most economical TPP) will be redirected to another transmission line or distributed to several transmission lines. The consequences of this scenario are greater transmission losses and, hence, reduced overall efficiency and increased fuel costs.
In other words, the transmission line constraint limits the production from the “optimal” hydroelectric plant or thermal power plant and forces production from another thermal power plant that may have a higher HR, which significantly affects the overall fuel costs.
On the other hand, the power that the solar power plant injects into the system also has an impact on power flows, but also on TPPs and HPPs. In other words, in some intervals, the output power of the SPP is almost maximum and, in other periods, it is zero. In order to cover such a sudden change, the TPPs are exposed to a sudden change in their output power, i.e., their ramp rate limit is directly affected, and thus to their optimal power, i.e., total fuel costs. In this case, the cost sensitivity of the ramp rate limit is not so pronounced due to the implementation of CCGT in the system. The advantage of CCGT is that it has a negligibly low ramp rate limit, which significantly contributes to maintaining the stability of the system, as well as reducing overall fuel costs.
From all this, another hypothesis can be confirmed. It means that the solution for the SPP installation with an orientation that “covers” the peaks in the electricity price diagram will give the lowest total fuel costs.

6. Conclusions

The optimal generation and optimal power flow in SHTSS are of the most important goals in the planning and operation of modern power systems. The use of power electronic devices and systems expands the possibilities of power regulation and control of power plants. In this aspect, in recent years, more and more research has focused on the solution of this problem. Unfortunately, in most of the published manuscripts, a number of factors that affect the quality of the obtained results are neglected: CCGT, security constraints (valve point effect, ramp rate, spinning reserve, and transmission line constraint), and power flow calculations. In this sense, this paper focuses on solving the problem of optimal short-term hydrothermal planning in a more realistic manner, i.e., focuses on the application of CCGT and solar power plants in short-term hydrothermal planning. This was achieved by introducing the security and power flow constraints into the SHTSS mathematical model. A novel GA algorithm is proposed that uses a self-adaptive mechanism, LX operator, MPTM operator, adaptive crossover, and mutation strategy to solve the SHTSS problem under an imposed DC power flow constraint. The complex constraints of the SHTSS problem are solved by using a new approach to their handling. The optimization results obtained by using the proposed SAGA algorithm on a known hydrothermal test system are compared with those obtained by other algorithms. The comparisons convincingly proved that the proposed SAGA performs better than other algorithms in terms of both determining the optimal value and the average CPU time used to solve the problem. After verifications of traditional STHS problems, SAGA is used for SHTSS problem of IEEE 30 bus systems considering transmission networks, CCGT and SPP. These results also show the better qualities of the proposed SAGA. To better illustrate the impact of different parameters on the final solutions, the proposed constraint handling techniques are neglected and the hydrothermal test system is solved by classical GA. The numerical results obtained from this test system show that LX, MPTM, self-adaptive fitness, and adaptive crossover and mutation strategy also play an important role in the SAGA algorithm. Numerical results show that using these parameters would significantly increase the feasible solutions and reduce the total fuel costs. All this has an extraordinary impact in the current state of the energy crisis and constantly increasing prices of energy carriers.
On the other hand, taking into account the expansion of the share of produced energy from decentralized sources and that, in accordance with the EU directives, there is a tendency towards the gradual decommissioning of thermal power plants, the proposed algorithm can be successfully applied to a hydro-thermal-solar system. Therefore, the proposed metaheuristic algorithm gives a clear picture of the operational planning of the system in the phase of gradual replacement of classical TPPs with combined cycle units (CCGT), but also the growing influence of large SPPs on total fuel costs and overall system parameters.
Possible disadvantages of the proposed algorithm can be the long time required for its modeling, as well as the strong intercorrelation of the proposed techniques that could affect the output parameters. However, keeping in mind the theorem “there is no free lunch”, this should not affect the positive outcomes, such as solving a very complex optimization problem, avoiding local optima, and obtaining a better and more physically realistic solution compared to other methods.
On the other hand, the author performs further research and upgrades to the mathematical model of the proposed algorithm in order to solve the AC short-term hydro-thermal-solar scheduling, i.e., considering voltage bus constraint, reactive power constraint, and reactive power flow limit. This means that the Newton–Raphson’s method will be implemented in the algorithm for power flow solving. This model will enable considering the stochastic nature of SPPs, and objective function will be extended with a penalty function for solar irradiation.
Therefore, the developed algorithm and the presented opportunities for its development allow a much more realistic and physically feasible solution to the SHTSS optimization problem to be obtained. This is also necessary for optimal planning of the operation of the power system from the point of view of investments, energy security, and reliability.

Author Contributions

B.P., N.H., A.I., and D.D. were involved in the full process of producing this paper, including conceptualization, methodology, modeling, validation, visualization, and preparing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Bulgarian National Scientific Fund, grant number KII-06-Н57/7/16.11.2021 and the APC was funded by KII-06-Н57/7/16.11.2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was carried out within the framework of the project “Artificial Intelligence-Based modeling, design, control and operation of power electronic devices and systems”, KII-06-H57/7/16.11.2021, Bulgarian National Scientific Fund.

Conflicts of Interest

The authors declare no conflict of interest.

Glossary

Symbols
NTnumber of TPPs
NHnumber of HPPs
NSnumber of SPPs
Ftotal costs, EUR
PGT,toutput power of TPP t, MW
PGH,houtput power of HPP h, MW
PGS,soutput power of SPP s, MW
PPsystem load, MW
PLtransmission loses, MW
Tjduration of interval j, h
Jduration of optimization period, h
Rspinning reserve, MW
PP,maxsystem’s peak load, MW
PGminunit technical minimum, MW
PGmaxunit technical maximum, MW
Bij, Bi0, B00Crohn’s B coefficients for losses
Qthwater discharge of hydro unit h, m3/h
P GT , t AV available spinning reserve of TPP t, MW
P GH , h AV available spinning reserve of HPP h, MW
RTtotal required spinning reserve of TPPs, MW
RHtotal required spinning reserve of HPPs, MW
Tdisdischarge time, h
URT, URHup rate of TPP and HPP, MW
DRT, DRHdown rate of TPP and HPP, MW
PGR,gactive power of transmission line g, MW
P GR , g max maximum transmission capacity of transmission line g, MW
Gnumber of transmission lines
Vh,kavailable water volume of HPP h, 103 m3
Wmax,itotal available energy of generator i, MWH
Vh,jstorage volume of HPP h at interval j, 103 m3
Ih,jinflow in reservoir h at interval j, 103 m3
Sh,jwater spillage of reservoir h at interval j, 103 m3
Vh0, Vh24initial and final volume of reservoir h, 103 m3
Vhmin, Vhmaxminimum and maximum volume of reservoir h, 103 m3
d(X), p(X)distance value and penalty value
v(X)total constraint violation
fsscaled fitness function
favaverage fitness function of current population
fmax, fminmaximal and minimal fitness function of current population
as, bsscaling coefficients
βL, aL, bLLX parameters
x1, x2parent chromosomes
y1, y2children chromosomes
p ^ , p, bmMPTM parameters
k1, k2, k3, k4constants of crossover strategy and mutation strategy
Pc, Pmcrossover probability and mutation probability
gi, hkinequality and equality constraints
Subscripts
jhour number
tTPP index
hHPP index
gtransmission line index
mgene index
nchromosome index
iinequality constraint index
kequality constraint index
Abbreviations
SHTSSshort-term hydro-thermal-solar scheduling
SAGAself-adaptive genetic algorithm
LXLaplace crossover
MPTMMakinen, Periaux, and Toivanen mutation
UCUnit Commitment

References

  1. Braun, S. Improving Flexibility of Fossil Fired Power Plants. Encycl. Energy Storage 2022, 2, 133–140. [Google Scholar]
  2. Zhao, X.; Yin, H.; Zhao, Y. Impact of environmental regulations on the efficiency and CO2 emissions of power plants in China. Appl. Energy 2015, 149, 238–247. [Google Scholar] [CrossRef]
  3. Schreiber, A.; Zapp, P.; Kuckshinrichs, W. Environmental assessment of German electricity generation from coal-fired power plants with amine-based carbon capture. Int. J. Life Cycle Assess. 2009, 14, 547–559. [Google Scholar] [CrossRef]
  4. Singh, R.; Banerjee, R. Impact of large-scale rooftop solar PV integration: An algorithm for hydrothermal-solar scheduling (HTSS). Sol. Energy 2017, 157, 988–1004. [Google Scholar] [CrossRef]
  5. Libal, J.; Kopecek, R. Bifacial Photovoltaics: Technology, Applications and Economics, 1st ed.; IET: London, UK, 2018; pp. 153–220. [Google Scholar]
  6. Tang, J.; Luh, P.B. Hydro thermal scheduling via extended differential dynamic programming and mixed coordination. IEEE Trans. Power Syst. 1995, 10, 2021–2028. [Google Scholar] [CrossRef]
  7. Jin-Shyr, Y.; Nanming, C. Short term hydro thermal coordination using multi pass dynamic programming. IEEE Trans. Power Syst. 1989, 4, 1050–1056. [Google Scholar] [CrossRef]
  8. Guan, X.; Ni, E.; Li, R.; Luh, P.B. An optimization-based algorithm for scheduling hydro thermal power systems with cascaded reservoirs and discrete hydro constraints. IEEE Trans. Power Syst. 1997, 12, 1775–1780. [Google Scholar] [CrossRef]
  9. Al-Agtash, S. Hydro thermal scheduling by augmented lagrangian: Consideration of transmission constraints and pumped storage units. IEEE Trans. Power Syst. 2001, 16, 750–756. [Google Scholar] [CrossRef]
  10. Nilsson, O.; Sjelvgren, D. Mixed integer programming applied to short term planning of a hydro thermal system. IEEE Trans. Power Syst. 1996, 11, 281–286. [Google Scholar] [CrossRef]
  11. Maturana, J.; Riff, M.C. Solving the short-term electrical generation scheduling problem by an adaptive evolutionary approach. Eur. J. Oper. Res. 2007, 179, 677–691. [Google Scholar] [CrossRef]
  12. Wong, D.P.; Wong, Y.W. Short term hydro thermal scheduling part. I. Simulated annealing approach. IEE Proc.-Gener. Transm. Distrib. 1994, 141, 497–501. [Google Scholar] [CrossRef]
  13. Jayabarathi, T.; Chalasani, S.; Shaik, Z.A. Hybrid differential evolution and particle swarm optimization-based solutions to short term hydro thermal scheduling. WSEAS Trans. Power Syst. 2007, 2, 245–254. [Google Scholar]
  14. Diew, V.N.; Ongsakul, W. Enhanced merit order and augmented Lagrange Hopfield network for hydro thermal scheduling. Int. J. Electr. Power Energy Syst. 2008, 30, 93–101. [Google Scholar]
  15. Gil, E.; Bustos, J.; Rudnick, H. Short term hydrothermal generation scheduling model using genetic algorithm. IEEE Trans. Power Syst. 2003, 18, 1256–1264. [Google Scholar] [CrossRef]
  16. Titus, S.; Jeyakumar, A.E. Hydrothermal scheduling using an improved particle swarm optimization technique considering prohibited operating zones. Int. J. Soft Comput. 2007, 2, 313–319. [Google Scholar]
  17. Liu, C.; Shahidehpour, M.; Li, Z.; Fotuhi-Firuzabad, M. Component and Mode Models for the Short-Term Scheduling of Combined-Cycle Units. IEEE Trans. Power Syst. 2009, 24, 976–990. [Google Scholar]
  18. General Electri GE Power. Combined Cycle Power Plant. Available online: www.ge.com/power/resources/knowledge-base/combined-cycle-power-plant-how-it-works (accessed on 14 July 2022).
  19. Kehlhofer, R.; Rukes, B.; Hannemann, F.; Stirnimann, F. Combined-Cycle Gas & Steam Turbine Power Plants, 3rd ed.; PennWell: Tulsa, OK, USA, 2009; pp. 35–68. [Google Scholar]
  20. Postolov, B.; Iliev, A. New metaheuristic methodology for solving Security Constrained Hydrothermal Unit Commitment based on Adaptive Genetic Algorithm. Int. J. Electr. Power Energy Syst. 2022, 134, 107163. [Google Scholar] [CrossRef]
  21. Zhu, J. Optimization of Power System Operation, 2nd ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2015; pp. 297–346. [Google Scholar]
  22. Soliman, A.; Mantawy, A.H. Modern Optimization Techniques with Applications in Electric Power Systems, 1st ed.; Springer: New York, NY, USA, 2012; pp. 83–280. [Google Scholar]
  23. Ongsakul, W.; Vo, D.N. Artificial Intelligence in Power System Optimization, 1st ed.; CRC Press: New York, NY, USA, 2013; pp. 118–314. [Google Scholar]
  24. Postolov, B.; Iliev, A. Adaptive Genetic Algorithm for Hydro-thermal Unit Commitment Considering the Security Constraints. Int. J. Electr. Eng. Comput. 2020, 4, 61–69. [Google Scholar] [CrossRef]
  25. Michalewicz, Z. Genetic Algorithms + Data Structures = Evolution Programs, 3rd ed.; Springer: New York, NY, USA, 1996; pp. 283–314. [Google Scholar]
  26. Deb, K. An efficient constraint handling method for genetic algorithms. Comput. Meth. Appl. Mech. Eng. 2000, 186, 311–338. [Google Scholar] [CrossRef]
  27. Qin, A.K.; Huang, V.L.; Suganthan, P.N. Differential evolution algorithm with strategy adaptation for global numerical optimization. IEEE Trans. Evol. Comput. 2009, 13, 398–417. [Google Scholar] [CrossRef]
  28. Mallipeddi, R.; Suganthan, P.N. Ensemble of constraint handling techniques. IEEE Trans. Evol. Comput. 2010, 14, 561–579. [Google Scholar] [CrossRef]
  29. Farmani, R.; Wright, J.A. Self-adaptive fitness formulation for constrained optimization. IEEE Trans. Evol. Comput. 2003, 7, 445–455. [Google Scholar] [CrossRef]
  30. Tessema, B.G.; Yen, G.G. A Self Adaptive Penalty Function Based Algorithm for Constrained Optimization. In Proceedings of the IEEE International Conference on Evolutionary Computation, Vancouver, BC, Canada, 16–21 July 2006. [Google Scholar]
  31. Sivanandam, S.N.; Deepa, N. Introduction to Genetic Algorithms, 1st ed.; Springer: Berlin, Germany, 2008; pp. 83–122. [Google Scholar]
  32. Kramer, O. Genetic Algorithm Essentials, 1st ed.; Springer International Publishing: Cham, Switzerland, 2017; pp. 11–28. [Google Scholar]
  33. Bansal, J.C.; Singh, P.K.; Pal, N.R. Evolutionary and Swarm Intelligence Algorithms, 1st ed.; Springer International Publishing: Cham, Switzerland, 2019; pp. 61–85. [Google Scholar]
  34. Deep, K.; Thakur, M. A new crossover operator for real coded genetic algorithms. Appl. Math. Comput. 2007, 188, 895–911. [Google Scholar] [CrossRef]
  35. Makinen, R.A.E.; Periaux, J.; Toivanen, J. Multidisciplinary shape optimization in aerodynamics and electromagnetic using genetic algorithms. Int. J. Numer. Methods Fluids 1999, 30, 149–159. [Google Scholar] [CrossRef]
  36. Basu, M. Artificial immune system for fixed head hydrothermal power system. Energy 2011, 36, 606–612. [Google Scholar] [CrossRef]
  37. Postolov, B.; Iliev, A.; Dimitrov, D.; Dimishkoska, N. N–1 Security Constrained Short-Term Hydrothermal Scheduling by Self Adaptive Genetic Algorithm with PTDF. In Proceedings of the 2021 International Conference on Information Technologies (InfoTech), Varna, Bulgaria, 16–17 September 2021; pp. 1–7. [Google Scholar]
  38. Khunkitti, S.; Watson, N.R.; Chatthaworn, R.; Premrudeepreechacharn, S.; Siritaratiwat, A. An Improved DA-PSO Optimization Approach for Unit Commitment Problem. Energies 2019, 12, 2335. [Google Scholar] [CrossRef]
  39. Kamboj, V.K. A novel hybrid PSO–GWO approach for unit commitment problem. Neural Comput. Appl. 2016, 27, 1643–1655. [Google Scholar] [CrossRef]
  40. Bouchekara, H.R.E.H.; Chaib, A.E.; Abido, M.A.; El-Sehiemy, R.A. Optimal Power Flow Using an Improved Colliding Bodies Optimization Algorithm. Appl. Soft Comput. 2016, 42, 119–131. [Google Scholar] [CrossRef]
  41. Aghdam, F.H.; Hagh, M.T. Security Constrained Unit Commitment (SCUC) formulation and its solving with Modified Imperialist Competitive Algorithm (MICA). J. King Saud Univ.-Eng. Sci. 2019, 31, 253–261. [Google Scholar] [CrossRef]
  42. Surekha, P. Investigation of Efficient Bio Inspired Intelligent Paradigms for Solving Unique Constraint Based Optimization Problems. Ph.D. Thesis, Anna University, Chennai, India, 1 August 2014. [Google Scholar]
  43. Nguyen, T.T.; Vo, D.N.; Truong, A.V.; Ho, L.D. Meta-Heuristic Algorithms for Solving Hydrothermal System Scheduling Problem Considering Constraints in Transmission Lines. Glob. J. Technol. Optim. 2016, 7, 1000192. [Google Scholar]
  44. Lu, B. Short-Term Scheduling of Combined Cycle Units. IEEE Trans. Power Syst. 2004, 19, 1616–1625. [Google Scholar] [CrossRef]
  45. Trinasolar. PV Modules Vertex 500W TSM-DEG18MC.20(II) Bifacial, 1/3-Cut, MBB, 485-505W. Available online: https://www.trinasolar.com/en-glb/product/VERTEX-DEG18MC20II (accessed on 14 July 2022).
  46. HUPX. Historical Data. Available online: https://hupx.hu/en/market-data/dam/historical-data (accessed on 14 July 2022).
Figure 1. CCGT with two GTs and one ST.
Figure 1. CCGT with two GTs and one ST.
Energies 15 05989 g001
Figure 2. State transition diagram for CCGT with two GTs and one ST [17].
Figure 2. State transition diagram for CCGT with two GTs and one ST [17].
Energies 15 05989 g002
Figure 3. Simplified flowchart of the proposed SAGA.
Figure 3. Simplified flowchart of the proposed SAGA.
Energies 15 05989 g003
Figure 4. Daily load diagram of the IEEE 30 bus system.
Figure 4. Daily load diagram of the IEEE 30 bus system.
Energies 15 05989 g004
Figure 5. Single-line diagram of the modified IEEE 30 bus system.
Figure 5. Single-line diagram of the modified IEEE 30 bus system.
Energies 15 05989 g005
Figure 6. HUPX’s electricity price diagram for the chosen optimization period (9 June 2021).
Figure 6. HUPX’s electricity price diagram for the chosen optimization period (9 June 2021).
Energies 15 05989 g006
Figure 7. Daily output power of the SPP for case 1.
Figure 7. Daily output power of the SPP for case 1.
Energies 15 05989 g007
Figure 8. Daily output power of the SPP for case 2.
Figure 8. Daily output power of the SPP for case 2.
Energies 15 05989 g008
Figure 9. Optimal generation scheduling for case 1–sub-case 1.
Figure 9. Optimal generation scheduling for case 1–sub-case 1.
Energies 15 05989 g009
Figure 10. Optimal generation scheduling for case 1–sub-case 2.
Figure 10. Optimal generation scheduling for case 1–sub-case 2.
Energies 15 05989 g010
Figure 11. Optimal generation scheduling for case 1–sub-case 3.
Figure 11. Optimal generation scheduling for case 1–sub-case 3.
Energies 15 05989 g011
Figure 12. Optimal generation scheduling for case 2–sub-case 1.
Figure 12. Optimal generation scheduling for case 2–sub-case 1.
Energies 15 05989 g012
Figure 13. Optimal generation scheduling for case 2–sub-case 2.
Figure 13. Optimal generation scheduling for case 2–sub-case 2.
Energies 15 05989 g013
Figure 14. Optimal generation scheduling for case 2–sub-case 3.
Figure 14. Optimal generation scheduling for case 2–sub-case 3.
Energies 15 05989 g014
Table 1. Optimal generation scheduling of the test system [37].
Table 1. Optimal generation scheduling of the test system [37].
IntervalPGT1
(MW)
PGT2
(MW)
PGH1
(MW)
PGH2
(MW)
1181.14264.66400.0085.11
2300.00340.23400.00217.02
3140.76309.24400.00300.00
Table 2. Optimal generation scheduling of the test system [37].
Table 2. Optimal generation scheduling of the test system [37].
SAGAGA 1AIS [30]DE [30]EP [30]AIS [30]
FT (EUR)47,18466,34166,11766,12166,198−40.12 (%)
CPU time (s)6.0328.5353.4360.7675.48−47.4 (s)
STD (EUR)26,861-----
1 This means that generation scheduling is obtained with Classical GA (modeled by the author, neglecting all newly proposed techniques from SAGA), to represent the difference between GA and SAGA [32].
Table 3. Optimal generation scheduling of the test system 2 [38,39].
Table 3. Optimal generation scheduling of the test system 2 [38,39].
Best (EUR)Average (EUR)Worst (EUR)Change (%)
DA-PSO [38]13,292.28---
PSO-GWO [39]13,600.00---
GA13,221.5513,288.1913,321.57-
SAGA13,126.3113,171.6913,193.28−1.264
Newton’s method-13,170.68--
Difference SAGA—Newton---0.0077
- This means that generation scheduling is obtained with Classical GA (modeled by the author, neglecting all newly proposed techniques from SAGA), to represent the difference between GA and SAGA.
Table 4. Composition data for CCGT.
Table 4. Composition data for CCGT.
ModeCompositionGTST P GT , t min
(MW)
P GT , t max
(MW)
URt
(MW)
DRt
(MW)
00 + 0000000
11 + 01020806060
22 + 02040160120120
31 + 11130135105105
42 + 12150215165165
Table 5. Cost curve coefficients data for CCGT.
Table 5. Cost curve coefficients data for CCGT.
ModeCompositionatbtctdtet
00 + 000000
11 + 0169.927002.369290.00051180.037
22 + 0339.854002.369290.00025180.037
31 + 1149.396301.400330.00063180.037
42 + 1247.069161.530060.00036180.037
Table 6. Data for case 1: Gaag, Netherlands.
Table 6. Data for case 1: Gaag, Netherlands.
Sub-CaseTilt Angle
(°)
Azimuth
(°)
Height
(m)
Sheds Spacing (m)AlbedoEGRID
(MWh)
12501.590.720,457
220901.5120.718,883
320−901.5120.718,734
Table 7. Data for case 2: Bitola, Macedonia.
Table 7. Data for case 2: Bitola, Macedonia.
Sub-CaseTilt Angle
(°)
Azimuth
(°)
Height
(m)
Sheds Spacing (m)AlbedoEGRID
(MWh)
12501.590.728,620
220901.5120.726,448
320−901.5120.726,481
Table 8. Obtained HR for TPPs (EUR/MW).
Table 8. Obtained HR for TPPs (EUR/MW).
CCGT-
Mode1
CCGT-Mode2CCGT-Mode3CCGT-Mode4PGT,2PGT,3PGT,4
HR4.714.642.682.773.274.243.83
Table 9. Total fuel cost for case 1.
Table 9. Total fuel cost for case 1.
Sub-Case 1Sub-Case 2Sub-Case 3
FT (EUR)10,316.1910,255.3410,320.99
Relative change compared to sub-case 2 (%)0.59-0.63
Table 10. Optimal generation scheduling using SAGA for case 1.
Table 10. Optimal generation scheduling using SAGA for case 1.
IntervalCCGT
(MW)
PGT2
(MW)
PGT3
(MW)
PGT4
(MW)
PGH1
(MW)
PGH2
(MW)
PGS
(MW)
PL
(MW)
Q1
(m3/h)
Q2
(m3/h)
171.05 m431.2018.2412.9722.9411.600.002.0258.02228.89
278.92 m438.2021.1717.0226.5116.860.002.7290.02321.49
310.62 m439.9315.0018.4223.1733.170.004.3260.10612.33
4132.28 m452.600.0026.7428.8133.610.007.0310.80620.28
5130.26 m454.4915.0028.8028.2132.421.006.8305.38598.82
6119.02 m449.2923.3026.5428.8528.621.985.6311.17530.74
7110.97 m441.3025.6719.4327.5022.253.444.6298.97417.10
8102.65 m434.1121.3313.9123.8014.606.203.6265.75281.69
993.38 m428.3817.0312.2823.4611.408.942.9262.75225.42
1078.59 m320.0015.0010.0018.3410.0011.011.9217.03200.81
1165.18 m320.0015.0010.0015.4810.0012.781.4191.71200.81
1274.21 m320.0015.0010.0018.6410.0013.931.8219.67200.81
1378.66 m320.0015.0010.9022.5110.2014.712.0254.22204.35
1489.29 m424.4915.0011.0822.5910.0015.092.5254.96200.81
1594.16 m431.6619.6813.5524.3913.2714.363.1271.06258.16
16100.44 m438.0524.4714.4627.1517.6513.513.7295.79335.43
17107.23 m440.8824.1817.7927.8720.8511.484.3302.31392.12
18109.53 m440.2224.5817.0026.5619.048.414.3290.49360.05
19103.00 m439.7826.9119.5228.2818.793.724.0305.97355.70
20111.21 m437.0923.2314.3525.9417.320.134.3284.99329.58
2197.33 m434.0522.1014.1824.6515.040.003.3273.35289.44
2297.02 m427.5416.1010.7122.9910.640.003.0258.53212.02
2388.13 m320.0015.0010.0020.2210.000.002.3233.75200.81
2496.03 m30.000.000.0023.1514.120.002.3259.97273.16
Table 11. Total fuel cost for case 2.
Table 11. Total fuel cost for case 2.
Sub-Case 1Sub-Case 2Sub-Case 3
FT (EUR)10,397.0210,345.0210,184.02
Relative change compared to sub-case 3 (%)2.051.55-
Table 12. Optimal generation scheduling using SAGA for case 2.
Table 12. Optimal generation scheduling using SAGA for case 2.
IntervalCCGT
(MW)
PGT2
(MW)
PGT3
(MW)
PGT4
(MW)
PGH1
(MW)
PGH2
(MW)
PGS
(MW)
PL
(MW)
Q1
(m3/h)
Q2
(m3/h)
171.68 m430.8717.5113.9922.7011.270.002.0255.94223.06
271.09 m438.4924.7915.8527.2121.010.002.4296.33395.08
393.38 m439.9829.6220.8128.6020.160.003.6308.92379.87
4121.68 m445.7327.0523.0028.8526.040.175.5311.17484.57
5127.31 m446.6532.0623.5028.8526.464.546.0311.17492.11
6117.57 m446.2725.1924.0528.6325.4010.175.3309.13473.13
7111.55 m441.7616.0717.0328.0923.0813.014.6304.30431.72
897.86 m437.230.019.3328.4519.4414.463.8307.51367.21
9103.39 m428.670.012.4422.6913.1215.163.5255.79255.52
1063.04 m320.2415.0012.5223.4912.6715.511.5262.97247.72
1159.45 m320.0015.0010.0018.3010.0015.521.3216.72200.81
1267.26 m320.0016.5811.2221.3410.9414.241.6243.72217.30
1383.41 m322.1515.0010.0019.1310.0012.512.2224.02200.81
1494.72 m425.170.0010.0017.6130.4810.273.2210.51564.10
15103.06 m429.6315.0015.8024.5616.217.243.5272.52310.09
16119.02 m420.9824.3516.2627.4724.543.694.3298.69457.81
17124.51 m40.0033.4528.5628.8532.971.984.3311.17608.71
18114.34 m420.5135.2019.5028.8526.440.324.2311.17491.76
19102.08 m430.8131.3521.8028.8524.910.003.8311.17464.37
20111.04 m431.9124.5516.0126.7218.870.004.1291.92357.16
2198.91 m434.5121.2713.9824.1114.630.003.4268.57282.26
2283.69 m328.5820.3613.1324.9213.830.002.5275.81268.05
2388.07 m320.0015.0910.0020.1310.050.002.3232.98201.63
2487.82 m30.0015.000.0020.0510.000.001.9232.23200.81
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Postolov, B.; Hinov, N.; Iliev, A.; Dimitrov, D. Short-Term Hydro-Thermal-Solar Scheduling with CCGT Based on Self-Adaptive Genetic Algorithm. Energies 2022, 15, 5989. https://doi.org/10.3390/en15165989

AMA Style

Postolov B, Hinov N, Iliev A, Dimitrov D. Short-Term Hydro-Thermal-Solar Scheduling with CCGT Based on Self-Adaptive Genetic Algorithm. Energies. 2022; 15(16):5989. https://doi.org/10.3390/en15165989

Chicago/Turabian Style

Postolov, Borche, Nikolay Hinov, Atanas Iliev, and Dimitar Dimitrov. 2022. "Short-Term Hydro-Thermal-Solar Scheduling with CCGT Based on Self-Adaptive Genetic Algorithm" Energies 15, no. 16: 5989. https://doi.org/10.3390/en15165989

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

Article Metrics

Back to TopTop