Next Article in Journal
Seismic Performance of a Green Roof Structure
Next Article in Special Issue
A Stability Improvement Method of DC Microgrid System Using Passive Damping and Proportional-Resonance (PR) Control
Previous Article in Journal
Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region
Previous Article in Special Issue
An Advanced Control Technique for Power Quality Improvement of Grid-Tied Multilevel Inverter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem

by
Cui Zheyuan
1,2,†,
Ali Thaeer Hammid
3,†,
Ali Noori Kareem
4,
Mingxin Jiang
5,
Muamer N. Mohammed
6 and
Nallapaneni Manoj Kumar
7,*
1
School of Computer Science, Baoji University of Arts and Sciences, Baoji 721007, China
2
Faculty of Management and Economics, Universiti Pendidikan Sultan Idris, Tanjong Malim 35900, Malaysia
3
Computer Engineering Techniques Department, Faculty of Information Technology, Imam Ja’afar Al-Sadiq University, Baghdad 10012, Iraq
4
Computer Engineering Department, Bilad Alrafidain University College, Ba’aqubah, Diyala 32001, Iraq
5
Faculty of Electronic Information Engineering, Huaiyin Institute of Technology, Huaian 223003, China
6
Information Technology Department, Community College of Qatar (CCQ), Doha 00974, Qatar
7
School of Energy and Environment, City University of Hong Kong, Kowloon, Hong Kong, China
*
Author to whom correspondence should be addressed.
Both have contributed equally (first and co-first author).
Sustainability 2021, 13(8), 4277; https://doi.org/10.3390/su13084277
Submission received: 27 February 2021 / Revised: 9 April 2021 / Accepted: 10 April 2021 / Published: 12 April 2021

Abstract

:
The key criteria of the short-term hydrothermal scheduling (StHS) problem is to minimize the gross fuel cost for electricity production by scheduling the hydrothermal power generators considering the constraints related to power balance; the gross release of water, and storage limitations of the reservoir, and the operating limitations of the thermal generators and hydropower plants. For addressing the same problem, numerous algorithms were being used, and related studies exist in the literature; however, they possess limitations concerning the solution state and the number of iterations it takes to reach the solution state. Hence, this article proposes using an enhanced cuckoo search algorithm (CSA) called the rigid cuckoo search algorithm (RCSA), a modified version of the traditional CSA for solving the StHS problem. The proposed RCSA improves the solution state and decreases the iteration numbers related to the CSA with a modified Lévy flight. Here, the movement distances are divided into multiple possible steps, which has infinite diversity. The effectiveness of RCSA has been validated by considering the hydrothermal power system. The observed results reveal the superior performance of RCSA among all other compared algorithms that recently have been used for the StHS problem. It is also observed that the RCSA approach has achieved minimum gross costs than other techniques. Thus, the proposed RCSA proves to be a highly effective and convenient approach for addressing the StHS problems

1. Introduction

The short-term hydrothermal scheduling (StHS) aims to minimize thermal units gross fuel cost by optimally scheduling the thermal generators and hydropower units’ operation for a given time. The support of numerous algorithms does the optimal scheduling, and there exist many studies in the literature. Owing to the StHS objective function’s nonlinear nature, Lagrange multipliers and gradient techniques must be employed. On the other side, considering the linearized behavior produces non-optimal solutions that appear in immense revenue generation losses, which was also done through scheduled operations [1,2]. In the StHS, the electrical load sharing requests have been optimally specified between thermal generators and hydropower units by appropriately planning the hydrothermal system’s operation, meeting the StHS key criteria, i.e., minimum fuel cost [3].
As mentioned earlier, many studies used optimization techniques to address the StHS problem. In ref. [3], a detailed summary of those studies along with a brief discussion on various optimization techniques and algorithms can be found. Apart from those, there exist few recent studies on the StHS problem that used a genetic algorithm (GA) [4], enhanced GA (EGA), particle swarm optimization (PSO), and enhanced PSO (EPSO) [5], classical evolutionary programming (CEP), fast evolutionary programming (FEP), and improved FEP (IFEP) [6], grasshopper optimization algorithm (GOA) [7], adaptive particle swarm optimization (APSO), modified APSO (MAPSO) [7], modified differential evolution (MDE) [8], improved PSO (IPSO) [9], teaching learning-based optimization (TLBO) [10], one rank cuckoo search algorithm (ORCSA) [11], running IFEP (RIFEP) [12], gradient search techniques (GS) [13], simulated annealing approach (SA) [14], clonal selection roy (CSR) [15], krill herd algorithm (KHA) [16], and sequential quadratic programming (SQP) [17].
All the above-suggested algorithms have contributed significantly in addressing the StHS problem; however, they possess limitations concerning the solution state and the number of iteration it takes to reach the solution state. Lately, the cuckoo search algorithm (CSA) has been proposed for economic dispatch problems; however, Yang and Deb, in 2009, used the CSA for the first time to solve optimization problems [18]. CSA is one of the metaheuristic algorithms supporting several rule parameters. It mimics the parasitism of several cuckoo kinds by depositing their eggs in the host bird nests for other kinds. Upon seeing the advantage of CSA in solving the optimization problem-solving ability, it was later used to solve non-convex and economic dispatch problems [19]. More recently, the CSA was even used for addressing the StHS problem; see in ref [20] for details. The results of ref. [20], suggested that CSA is a practical approach with improved performance, which is better than all other cases containing high-rate nonlinearity like valve point loads. However, the traditional CSA algorithm suffers from a limitation of step-length variation, which is crucial in attaining the solution. Hence, this article proposes using an enhanced CSA called rigid cuckoo search algorithm (RCSA), a modified version of the traditional CSA for solving the StHS problem.
In the proposed RCSA technique, we consider a modified Lévy flight in which the step-lengths have been classified based on its new specific possibility configuration, which produces unlimited variation. The sequential jumps or steps of a cuckoo necessarily create a stochastic step method that adopts a power law for step length configuration resulting in large final step-sizes. Hence, it is wise to make a stochastic step in a biased process with many stochastic step sizes. As a result, RCSA for the StHS optimization problems enhances the optimization convergence accuracy and rate. Additionally, the effectiveness of RCSA has been examined and validated on complicated optimization problems of hydrothermal scheduling.
The article is structured in three sections, where Section 2 provided the detailed mathematical modeling of the hydrothermal power system and the RCSA implementation. In Section 3, the results are discussed and compared with other algorithms. Finally, the concluding remarks were provided in Section 4.

2. Methods

2.1. Mathematical Model of the Hydrothermal System

This section provides a mathematical model of the hydrothermal power system that we used for optimization. Between the two generators (hydro and thermal), the input fuel cost is negligible in hydropower generating units due to water availability as a free source. However, it is different from thermal power generators; hence, our main aim was to minimize the gross input fuel cost of thermal power generating units while producing the energy by using hydro resources significantly. Considering the above condition, we selected the objective function that is shown in Equation (1). Additionally, the considered constraints for solving the StHS problem are provided in Equations (2) to (9) [21,22,23].
Objective function:
m i n   F   =   j = 1 J f i ( P T j )
Constraints:
The constraint related to the balance of power generation–load is given by:
P T j + i = 1 I P H ( i , j ) = P d j + P l o s s   j
The hydropower generation (PH(i,j)) is a rate function of water release and is given by:
P H ( i , j ) =   ϕ ( q )
The water volume stored in the reservoir can be given by:
X i ( j + 1 ) = X i j q i ( j + 1 ) s i ( j + 1 ) + r i ( j + 1 )
Operational periods of thermal power generators have been restrained according to their capacities are given by:
P T , m i n P T P T , m a x
Operational periods of hydropower generators have been restrained according to their capacities are given by:
P H , m i n P H P H , m a x
The constraints related to water release rate limitations are represented by:
q i , m i n q i , j q i , m a x
The constraints related to initial and final reservoir water volume are given by
X i 0 = X i ( 0 ) ,   X j 0 = X j ( 0 )
The constraints related to reservoir limitations of water storage are given by
X i , m i n X i , j X i , m a x

2.2. Cuckoo Search Algorithm and Lévy Flights

The CSA is a metaheuristic search algorithm developed by Yang and Deb [18]. Compared to other algorithms, the CSA is a unique technique based on a heuristic evolutionary algorithm’s population to solve optimization problems like the StHS problem with easy execution and few adjustable parameters. The CSA algorithm rule is related to the compulsory brood parasitic performance and Lévy flight performance of several birds and fruit flies. Below are three essential concepts that are helpful throughout the rule exploration.
  • Concept 1. Individually all cuckoos produce one single egg at a time that regards a proposed solution and randomly throws its egg up into the wanted nest among the set number of possible host nests.
  • Concept 2. The egg of high quality thrown up in the best nests regards a better solution transferred to the subsequent generation.
  • Concept 3. The possible number of host nests is constant, and the probability that the host bird can find a nest is indicated by the probability constant, Pa, with range [0,1]. Hence, it may either discard the egg or leave this nest and then build a new nest entirely in a different place.
Based on the above-listed principles in the order of steps, cuckoos naturally explore food options in a random or semi-random way. A cuckoo forgoing route is adequately an unexpected step due to the subsequent progress based on both the present position and the passage possibility to the next position. These unexpected steps can be modeled mathematically, and in the literature, it was explained that the flight performance could illustrate the ideal aspects by the Lévy flights concept [24]. A Lévy flight is a random step in which every step’s length is classified based on a heavy-tailed possibility configuration. Afterward, the length from the start point of the random steps leads to a steady configuration in some steps.

2.3. Rigid Cuckoo Search Algorithm

According to the three rules mentioned previously in Section 2.2, for finding the generation of new solutions, i.e., S i ( t + 1 )   for the i th cuckoo, the Lévy flight shown in Equation (10) is achieved.
S i ( t + 1 ) = S i ( t ) + α L e v y ( λ )
where α must be more than zero 0 (for all of the step size), according to the range of the interest problem.
The product indicates elementwise multiplication. This study considers Lévy flight in which the step-lengths have been distributed due to the subsequent possibility allocation, which has an unlimited variation, see in Equation (11)
L e v y   u   =   t λ ,   0.999   <   λ     3.079
Here, the steps made by any single cuckoos are typically represented in a stochastic manner that follows a power-law of allocation for the step-length with a thick tailpiece. It is worth noting that, if egg cuckoos were quite similar to the host eggs, then these egg cuckoos’ have few possibilities, which can be identified. Therefore, the fitness must relate to the variance solutions. Hence, the RCSA does entire stochastic steps in a biased approach with several unplanned step dimensions for obtaining the solutions. There are three parameters P a , λ , and α , which are included in the CSA to support this technique to obtain universally and regionally enhanced solutions. The parameters P a and α are quite significant in fine-tuning vectors of the solution and can be employed to set the algorithm convergence rate. The conventional CSA algorithm uses established values for both P a and α , defined in the initialization stage and cannot be replaced later. Additionally, the iterations number to obtain an optimal solution is another disadvantage of conventional CSA. If there were the low-value of P a and the high-value of α , this technique’s execution would be weak and manages to a significant and result in the iterations number increment. Moreover, if there were the high-value of P a and the low-value of α , the convergence speed would be immense, but it probably cannot obtain the best solutions. Hence, we used the RCSA approach to solve the StHS problem; the flow chart presented in Figure 1 illustrated solving methodology.
The essential disparity between the improved CSA (i.e., RCSA) and CSA is modifying P a and α . To enhance the CSA technique’s achievement and minimize the disadvantages that lie on established values of P a and α , the improved CSA technique employs variable values for P a and α . Hence, in the new steps, the values of P a and α obligate to be large and execute the technique to improve the solution vectors. Nevertheless, these values must be reduced in definitive productions to succeed in a more regular fine-tuning of solution vectors. The values of Pa and α are physically replaced with the production number, see in Equations (12)–(14)
P a ( I T ) = P a , m a x P a , m a x P a , m i n I c × I T
α ( I T ) = α m a x . e ( c . I T )
c = L n ( α m a x α m i n ) I c
However, it becomes difficult to find the new solutions with the restrained optimization problem. Hence, there is a need for converting the restrained optimization problem to unrestrained and is as follows:
The optimization problem is expressed as follows in Equation (15)
M i n   f i   ( P T j ) ,   s . t .   :   R j 0   j   =   1 ,   2 ,   ,   N P
So, in such cases, the impossible solution must be adjusted to be possible solutions, which prevents the severity of restrictions. For this matter, a discipline function has been used to change the restrained optimization problem in Equation (16) to the unrestrained one in Equation (17)
M a x   F *   =   f i   ( P T j ) λ i = 1 N P m i n ( 0 ,   R j )
Based on the discipline coefficient, the maximization problem of M a x   F * is made equal to M i n   F * . It is represented as follows, see in Equation (17).
M i n   F *   =   f i   ( P T j ) + λ i = 1 N P m i n ( 0 ,   R j )

2.4. Implementation of RCSA on a Hydrothermal System

Step-1.
In the StHS problem, the influential variables such as the release rate of water for the whole plants for several hours and thermal unit production for the entire period are chosen irregularly within the operating limitations. The storage capacity of every reservoir has been estimated using Equation (4), the generation of hydro plants has been calculated using Equation (3). Subsequently, the thermal power generation has been computed by applying Equation (2). The population of the host nest (NE) has been explained as:
Y = [ Y 1 ,   Y 2 ,   ,   Y N E ] L
where every nest Yi is expressed as:
Y = [ q 1 , 1 i q 1 , j i q 1 , N h i P T 1 , 1 i P T 1 , j i P T 1 , N h i q k , 1 i q k , j i q k , N h i P T k , 1 i P T k , j i P T 1 , N h i q N h , 1 i q N h , j i q N h , n h i P T N h , 1 i P T N h , j i P T N h , n h i ]
Step-2.
Set the production number.
Step-3.
Compute the objective function using Equation (17). With the equation of restraints, many restrictions irregularly have been limited. Then, enhanced fuel cost has been calculated as Equation (19).
F * * = F * + k = 1 N C ( λ k × V i o k 2 )
Step-4.
The modern solution has been created by using Levy flights. The new solution’s computation has been built in the preceding best nest by using Levy flights. For this technique, the optimal way for levy flights has been computed by Yang XS’s contribution, Deb S [18]. The new solution has been presented in Equation (20)
Y i , n e w = Y i , b e s t + ( α × r a n d 2 × Δ Y i , n e w )
where α > 0 is the updated step size, rand2 is a typical number of allocated stochastic and ΔYi,new has been computed, see Equations (21) and (22):
Δ Y i , n e w = ε × ρ 1 ( ψ ) ρ 2 ( ψ ) × ( Y i , b e s t G b e s t )
ε = r a n d A | r a n d B | 1 ψ
where randA and randB are a couple of commonly allocated stochastic variables with a standard deviation ρ 1 ( ψ ) and ρ 2 ( ψ ) that has been determined by Equations (23) and (24), respectively.
ρ 1 ( ψ ) = [ Φ ( 1 + ψ ) × s i n ( π ψ 2 ) Φ ( 1 + ψ 2 ) × ψ × 2 ( ψ 1 2 ) ] 1 ψ
ρ 2 ( ψ ) = 1
where ψ is between the range from 0.29 to 2.01 and the obtained new solution must satisfy all associated constraints for both of ψ and Φ as well.
Step-5.
The effect of the detection of an alien egg in a nest of a host bird with Pa’s possibility produces a new solution for the problem comparable with the Levy flights. The new solution has been computed as following Equations (25)–(27):
Δ Y i ,   d i s = Y i , b e s t + ( k × Δ Y i ,   d i s )
where k was the modernized coefficient defined built in the possibility of a host bird to find out an alien egg in its nest:
k = { 1 i f   r a n d 3   >   P a 0 o t h e r w i s e
The increment value of Δ Y i ,   d i s has been defined by
Δ Y i ,   d i s = r a n d 3 × [ r a n d p 1 ( Y d ,   b e s t ) r a n d p 2 ( Y d ,   b e s t ) ]
where r a n d 3 is the allocated random number within [0,1], r a n d p 1 ( Y d ,   b e s t ) and r a n d p 2 ( Y d ,   b e s t )   are the random disorder for locations of nests in ( Y d ,   b e s t ) . Once more, for the recently produced solution, its lower and upper limits should please the unit’s limitation. The best value has been modernized for every nest ( Y d ,   b e s t ) and the nest identical to the best fitness function has been defined by Gbest.
Step-6.
The technique ends if modern production gives the maximum production number.

3. Results and Discussion

This section provides the results of an RCSA application in addressing the StHS problem; hence, a test system of hydrothermal power generating units was considered and is based on the refs [7,10]. It involves a complex of four hydro plants and some thermal units regarded as a single equivalent thermal plant. To understand the RCSA technique’s feasibility for a more extensive hydrothermal power system, it has been applied on a second test system that involved a cascade of three thermal units and four hydro plants. The real data of this system was collected based on refs [10,25]. The schedule period of 24 h, and every interval equal to 1 h time, has been considered. The simulation was carried out using MATLAB 9.8, and the computer used is a core i7-8th Gen processor with 2.00 GHz and 8.00 GB RAM.

3.1. Parameter Selection

In the RCSA, there were just five parameters that can be harmonized, involving three major items from the authentic CSA and a couple of other adjustments. Firstly, a couple of parameters that involved the three major items, which influence every recent solution produced from the exploration and exploitation are considered. These aspects are the nests’ number (NE) and the possibility of an alien egg to be found out, Pa. In contrast, the number of maximum iterations should have undeviatingly an effect on the optimal solution. Secondly, a couple of other parameters that influence joining the exploration aspect and exploitation aspects are considered. These should be satisfied with the upper and lower limitations and can be adjusted via the Lévy flights power. The obliged with the best solution provides the RCSA technique to enhance its performance and convergence speed. On the other side, the three major parameters from the authentic CSA technique, a couple of others in the justification, were effortless to be selected because they had been clarified by the previously limit equations. After several number of runs with various values of RCSA control parameter, the key control parameter chosen are population (Np) = 100, maximum iteration = 500 and value of probability (Pa) = 0.7.

3.2. Obtained Results

The proposed RCSA was executed more than ten confident times within limit range values of Pa from 0.1 up to 0.9, and a specific version of RCSA was achieved more than 100 confident times. In contrast, the nests’ number and maximum iterations number are previously groups to particular values of 10 and 300, respectively. The results contain minimal total cost, average gross cost, maximum gross cost, average calculation time, and standard deviation collected by RCSA, respectively illustrated in Table 1 and Table 2.
Based on the solutions present in Table 1 and Table 2, the CSA obtained optimal solutions at Pa from 0.1 up to 0.9, and RCSA obtained optimal solutions at Pa precisely equal to 0.7. Moreover, RCSA could get less maximum gross cost, less average gross cost, and less standard deviation.
The specific optimum points of the solution involving water release and thermal power production have been illustrated in Table 3. Thus, it is proved that the proposed RCSA technique effectively attains the solution for the StHS problem with cascaded hydropower plants. The load demand and the thermal and hydropower plants’ power for every interval of time throughout the scheduling horizon correlating to the best solution for test system 1 is illustrated in Figure 2. For the same aspect, Figure 3 displays the reservoir storage volumes of entire hydro plants and the cost convergence characteristic of the suggested technique.

3.3. Proposed System Validation

The obtained results of the proposed RCSA technique were compared with several other techniques (see Table 4). To maintain the correctness of this comparison, studies conducted for test system 1 were only considered. The gross cost obtained with RCSA was relatively equal to that produced by ORCSA–Lévy flights and ORCSA–Cauchy, and less than that produced by any other techniques. Nonetheless, the comparison results proved that the proposed RCSA technique is faster and more accurate in attaining solutions for StHS problems than the compared techniques.
Thus, the proposed technique is very efficient in making a solution for the StHS problem with reservoir volume constraints. Min/max cost is considered per day in USD, representing the best/worst fuel cost. The best/worst values are arranged by sorting from a more significant value to smallest depending on the minimum (best value) of fuel cost. Moreover, Table 4 represents the property of cost convergence for this suggested version of RCSA. The convergence comparison is shown that the RCSA technique was a more suitable optimal solution and more rapid than others.
In Table 5, the proposed RCSA technique’s comparative results with other techniques for test system 2 were provided. Based on the results, it is understood that a decrease in cost is observed with RCSA. Complete results, involved to hourly water release, thermal and hydro generation, of the best solution achieved by the suggested RCSA please whole types of constraints of StHS problem while decreasing the gross fuel cost-efficiently.
Besides the above comparison made in Table 4 and Table 5 for the test systems 1 and 2, we provided a comparative discussion with recent works available in the literature. Firstly, it is challenging to sketch them in a single unique figure because there are many output costs with various X–Y axis points. Thus, we provided the most suitable algorithms that are in comparison with the proposed one. H.M Dubey et al. [20] applied conventional CSA. It was found that the cost convergence characteristic is $2.5 × 107, which is then starting to reduce until it gets to $2.5 × 107. Additionally, the conventional CSA was employed in the same system but with various parameters and different specialized system constraints and fitness function (for example, the fixed-head). Here, the cost convergence characteristic was restricted from $2 × 1011 to 0.1 × 1011 [27]. T.T Nguyen used the modified CSA and D. N Vo in ref [29]; it was observed that the fitness function’s cost convergence started from $3.9 × 105 and got settled at $3.7 × 105.
Finally, Figure 4 describes the cost convergence characteristic for some contribution that has the same range of our proposed algorithm and this was done mainly to see the validation. Then again, T.T Nguyen has applied the same with one rank on CSA [11]; and an adaptive selective approach on CSA [30].

4. Conclusions

In this study, the RCSA technique was implemented to solve the StHS problems with various complex constraints. The same was tested with four cascaded hydropower plants and one thermal plant for the scheduled operation of 24 h with subintervals of 1 h. The results proved that the proposed RCSA technique was more effective than the conventional CSA for the StHS problem. Simulation results of the cascaded hydrothermal systems have shown that the RCSA technique achieved a satisfactory optimal solution where the computing time was lesser than other recent optimization techniques.

Author Contributions

Conceptualization, A.T.H., A.N.K., and N.M.K.; methodology, C.Z., A.T.H., and, N.M.K.; software, C.Z., and M.J.; validation, A.T.H., and N.M.K.; formal analysis, A.T.H., A.N.K., and N.M.K.; investigation, C.Z., A.T.H., A.N.K., and M.N.M.; data curation, A.T.H.; writing—original draft preparation, A.T.H., C.Z., and N.M.K.; writing—review and editing, N.M.K., A.N.K., M.J., and M.N.M.; visualization, N.M.K.; supervision, N.M.K.; project administration, N.M.K., M.J., M.N.M., funding acquisition, C.Z., M.J., M.N.M., and N.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was sponsored by Qing Lan Project of JiangSu Province (No.20200420), Opening Project of Henan Engineering Laboratory of Photoelectric Sensor and Intelligent Measurement and Control of Henan Polytechnic University (No.HELPSIMC-2020-002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

F Gross cost of production
Fi(PTj)Production cost for PTj
PTj Production of power for the thermal unit at period j
PH(i, j) Production of power for hydro unit i at period j
PdjSystem load demand at period j
qijWater release rate of hydro unit i at period j
rij Stream rate into the storage reservoir of the hydro plant at a period
sij Spillage of the reservoir at a period
XijStorage volume of hydro plant i at period j
Xi0 Initial reservoir storage of hydro plant
iNumber of units
nh, npMaximum number of unit hydro/thermal
jNumber of scheduling intervals
NH, NPMaximum number of scheduling periods of hydro/thermal
NCTotal number of constraints
RjResource constraint
ITTotal iteration numbers
IC Current iterations
λk Penalty value for kth
ViokViolation amount of kth constraint
ψDistribution factor
Φ Gamma distribution function
m Time period

References

  1. Hammid, A.T.; Bin Sulaiman, M.H.; Abdalla, A.N. Prediction of small hydropower plant power production in Himreen Lake dam (HLD) using artificial neural network. Alex. Eng. J. 2018, 57, 211–221. [Google Scholar] [CrossRef]
  2. Hammid, A.T.; Kadhim, M.B.H.; Sulaiman, A.A.; Kareem, A.N.; Jadaa, K.J. Priority of kaplan turbine and small hydropower plants over other resources: An overview. J. Eng. Appl. Sci. 2017, 12, 8695–8705. [Google Scholar]
  3. Thaeer Hammid, A.; Awad, O.I.; Sulaiman, M.H.; Gunasekaran, S.S.; Mostafa, S.A.; Manoj Kumar, N.; Khalaf, B.A.; Al-Jawhar, Y.A.; Abdulhasan, R.A. A review of optimization algorithms in solving hydro generation scheduling problems. Energies 2020, 13, 2787. [Google Scholar] [CrossRef]
  4. Orero, S.; Irving, M. A genetic algorithm modelling framework and solution technique for short term optimal hydrothermal scheduling. IEEE Trans. Power Syst. 1998, 13, 501–518. [Google Scholar] [CrossRef]
  5. Yuan, X.; Wang, L.; Yuan, Y. Application of enhanced PSO approach to optimal scheduling of hydro system. Energy Convers. Manag. 2008, 49, 2966–2972. [Google Scholar] [CrossRef]
  6. Sinha, N.; Chakrabarti, R.; Chattopadhyay, P. Fast evolutionary programming techniques for short-term hydrothermal scheduling. IEEE Trans. Power Syst. 2003, 18, 214–220. [Google Scholar] [CrossRef]
  7. Amjady, N.; Soleymanpour, H.R. Daily hydrothermal generation scheduling by a new modified adaptive particle swarm optimization technique. Electr. Power Syst. Res. 2010, 80, 723–732. [Google Scholar] [CrossRef]
  8. Lakshminarasimman, L.; Subramanian, S. Short-term scheduling of hydrothermal power system with cascaded reservoirs by using modified differential evolution. IEE Proc. Gener. Transm. Distrib. 2006, 153, 693–700. [Google Scholar] [CrossRef]
  9. Hota, P.; Barisal, A.; Chakrabarti, R. An improved PSO technique for short-term optimal hydrothermal scheduling. Electr. Power Syst. Res. 2009, 79, 1047–1053. [Google Scholar] [CrossRef]
  10. Roy, P.K. Teaching learning based optimization for short-term hydrothermal scheduling problem considering valve point effect and prohibited discharge constraint. Int. J. Electr. Power Energy Syst. 2013, 53, 10–19. [Google Scholar] [CrossRef]
  11. Nguyen, T.T.; Vo, D.N.; Ongsakul, W.; Trung, N.T. One rank cuckoo search algorithm for short-term hydrothermal scheduling with reservoir constraint. In Proceedings of the IEEE Eindhoven PowerTech, Eindhoven, The Netherlands, 29 June–2 July 2015. [Google Scholar]
  12. Türkay, B.; Mecitoğlu, F.; Baran, S. Application of a fast evolutionary algorithm to short-term hydro-thermal generation scheduling. Energy Sources Part B Econ. Planning, Policy 2011, 6, 395–405. [Google Scholar] [CrossRef]
  13. Wood, A.J.; Wollenberg, B.F.; Sheblé, G.B. Power Generation, Operation, and Control; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  14. Wong, K.; Wong, Y. Short-term hydrothermal scheduling part. I. Simulated annealing approach. IEE Proc. Gener. Transm. Distrib. 1994, 141, 497–501. [Google Scholar] [CrossRef]
  15. Swain, R.; Barisal, A.; Hota, P.; Chakrabarti, R. Short-term hydrothermal scheduling using clonal selection algorithm. Int. J. Electr. Power Energy Syst. 2011, 33, 647–656. [Google Scholar] [CrossRef]
  16. Das, S.; Bhattacharya, A. Symbiotic organisms search algorithm for short-term hydrothermal scheduling. Ain Shams Eng. J. 2018, 9, 499–516. [Google Scholar] [CrossRef] [Green Version]
  17. Sivasubramani, S.; Swarup, K.S. Hybrid DE–SQP algorithm for non-convex short term hydrothermal scheduling problem. Energy Convers. Manag. 2011, 52, 757–761. [Google Scholar] [CrossRef]
  18. Yang, X.-S.; Deb, S. Cuckoo search via Lévy flights. In Proceedings of the 2009 World Congress on Nature & Biologically Inspired Computing (NaBIC), Coimbatore, India, 9–11 December 2009. [Google Scholar]
  19. Vo, D.N.; Schegner, P.; Ongsakul, W. Cuckoo search algorithm for non-convex economic dispatch. IET Gener. Transm. Distrib. 2013, 7, 645–654. [Google Scholar] [CrossRef]
  20. Dubey, H.M.; Pandit, M.; Panigrahi, B. Cuckoo Search Algorithm for Short Term Hydrothermal Scheduling. In Power Electronics and Renewable Energy Systems; Springer: Berlin/Heidelberg, Germany, 2015; pp. 573–589. [Google Scholar]
  21. Xinchun, L.; Yongde, K.; Hongna, C.; Hui, L. Hydrothermal Effects of Freeze-Thaw in the Taklimakan Desert. Sustainability 2021, 13, 1292. [Google Scholar] [CrossRef]
  22. Zeng, X.; Hammid, A.T.; Kumar, N.M.; Subramaniam, U.; Almakhles, D.J. A grasshopper optimization algorithm for optimal short-term hydrothermal scheduling. Energy Rep. 2021, 7, 314–323. [Google Scholar] [CrossRef]
  23. Zhang, Z.; Fan, W.; Bao, W.; Chen, C.-T.; Liu, S.; Cai, Y. Recent Developments of Exploration and Detection of Shallow-Water Hydrothermal Systems. Sustainability 2020, 12, 9109. [Google Scholar] [CrossRef]
  24. Kanazawa, K.; Sano, T.G.; Cairoli, A.; Baule, A. Loopy Lévy flights enhance tracer diffusion in active suspensions. Nature 2020, 579, 364–367. [Google Scholar] [CrossRef] [Green Version]
  25. Nguyen, T.T.; Vo, D.N.; Dinh, B.H. A cuckoo bird-inspired meta-heuristic algorithm for optimal short-term hydrothermal generation cooperation. Cogent Eng. 2016, 3, 1266863. [Google Scholar] [CrossRef]
  26. Nguyen, T.T.; Vo, D.N.; Dao, T.T.; Dieu, V.N. Cuckoo search algorithm using different distributions for short-term hydrothermal scheduling with cascaded hydropower plants. In Proceedings of the TENCON 2014—2014 IEEE Region 10 Conference, Bangkok, Thailand, 22–25 October 2014. [Google Scholar]
  27. Nguyen, T.T.; Vo, D.N.; Truong, A.V. Cuckoo search algorithm for short-term hydrothermal scheduling. Appl. Energy 2014, 132, 276–287. [Google Scholar] [CrossRef]
  28. Basu, M. Economic environmental dispatch of hydrothermal power system. Int. J. Electr. Power Energy Syst. 2010, 32, 711–720. [Google Scholar] [CrossRef]
  29. Nguyen, T.T.; Vo, D.N. Modified cuckoo search algorithm for short-term hydrothermal scheduling. Int. J. Electr. Power Energy Syst. 2015, 65, 271–281. [Google Scholar] [CrossRef]
  30. Nguyen, T.T.; Vo, D.N.; Dinh, B.H. An effectively adaptive selective cuckoo search algorithm for solving three complicated short-term hydrothermal scheduling problems. Energy 2018, 155, 930–956. [Google Scholar] [CrossRef]
Figure 1. Flow chart showing the steps used in the rigid cuckoo search algorithm (RCSA) approach for solving the short-term hydrothermal scheduling (StHS) problem.
Figure 1. Flow chart showing the steps used in the rigid cuckoo search algorithm (RCSA) approach for solving the short-term hydrothermal scheduling (StHS) problem.
Sustainability 13 04277 g001
Figure 2. Hydro generation, thermal generation, and load demand for test system 1.
Figure 2. Hydro generation, thermal generation, and load demand for test system 1.
Sustainability 13 04277 g002
Figure 3. The reservoir storage volume of the hydro plants for test system 1.
Figure 3. The reservoir storage volume of the hydro plants for test system 1.
Sustainability 13 04277 g003
Figure 4. Convergence characteristic of the five various algorithms with RCSA methods.
Figure 4. Convergence characteristic of the five various algorithms with RCSA methods.
Sustainability 13 04277 g004
Table 1. Brief result from proposed RCSA with various values of Pa.
Table 1. Brief result from proposed RCSA with various values of Pa.
PaMin Cost ($)Avg. Cost ($)Max Cost ($)Std. dev. ($)CPU (s)
0.1709,932.115709,995.586710,655.745327.18917.9
0.2709,922.44710,768.874711,384.563599.38217.9
0.3709,911.445709,994.412710,989.321489.73317.8
0.4709,866.727709,745.236709,999.741103.93718.1
0.5709,886.651709,887.698709,988.25847.65318.6
0.6709,862.129709,899.987709,991.95154.51218.5
0.7709,862.027709,878.852709,996.15959.66118.2
0.8709,901.478709,910.357710,920.357478.22518.3
0.9709,902.685709,901.753710,901.753471.18518.7
Table 2. The optimal solutions achieved by the proposed RCSA technique.
Table 2. The optimal solutions achieved by the proposed RCSA technique.
mPDmVm
(acre-ft)
qm
(acre-ft/hr)
Psm
(MW)
Phm
(MW)
11199101,8971832892300
2149785,9593328892600
3109893,8471340892200
4179559,9984817892900
594870,4281124783158
6128959,9982863783509
Table 3. Optimal solution obtained by RCSA for test system 1.
Table 3. Optimal solution obtained by RCSA for test system 1.
HoursWater Discharge (×104 m3)Ps (MW)
Plant 1Plant 2Plant 3Plant 4
19.48846.137726.347713.10091036.249
29.40256.310425.208813.10241066.243
39.1976.25624.50413.10351049.837
48.80966.303823.883213.1003996.409
58.52696.255222.728713.1009975.264
68.33146.427821.923813.10071076.828
78.41016.82120622613.10051294.982
88.60817.230919.242813.10231623.787
98.80727.338718.577513.10771849.655
108.4657.578317.800613.11631917.3778
118.34297.718617.0931131816.169
128.41217.803216.996913.43021885.493
138.24447.862416.314414.67421786.317
148.12248.000115.43915.93711737.617
157.94178.150114.359617.29561652.815
167.81028.551812.906618.42441580.968
177.81538.934411.435320.00651630.883
187.66269.319310.100821.41571639.432
197.6729.965510.131623.02051739.29
207.594510.525210.104324.28431787.795
217.390811.187512.161825.09851761.32
227.58612.232712.74725.06361659.735
237.551313.138613.240325.09941411.135
247.410314.253313.666225.09721177.636
Table 4. Comparison of the results achieved by the suggested RCSA technique with others for test system 1.
Table 4. Comparison of the results achieved by the suggested RCSA technique with others for test system 1.
TechniquesMin Cost ($)Average Cost ($)Max Cost ($)CPU Time (s)
GA [4]942,600946,609.1951,0871920
EGA [5]934,727.00936,058.00937,339.00---
FEP [6]930,267.92930,897.44931,396.811911.2
CEP [6]930,166.25930,373.23930,927.012292.1
IFEP [6]930,129.82930,290.13930,881.921033.2
PSO [5]928,878.00933,085.00938,012.00---
CSA-Lévy [26]927,934.23927,980.45928,000.6679.08
CSA Cauchy [26]927,967.66927,981.49927,992.5381.30
CSA Gauss [26]927,957.26927,978.911928,003.2385.75
APSO [7]926,151.54---------
EPSO [5]922,904.00923,527.00924,808.00---
MDE [8]922,555.44---------
IPSO [9]922,553.49---------
MAPSO [7]922,421.66922,544.00923,508.00---
TLBO [10]922,373.39922,462.24922,873.81---
CSA [20]913,945.87917,624.024921,994.25---
RIFEP [12]709,862.05---------
GS [13]709,877.38---------
SA [14]709,874.36------901
ORCSA–Lévy flight [11]709,862.048------18
ORCSA–Cauchy [11]709,862.048------18
Proposed RCSA709,862.027------17
Table 5. Comparison of the results achieved by the suggested RCSA technique with others for test system 2.
Table 5. Comparison of the results achieved by the suggested RCSA technique with others for test system 2.
Techniques Min Cost ($)Average Cost ($)Max Cost ($)CPU Time (s)
SA [27]47,306------
CEP [27]45,466------
CEP-IFS [25]45,036.00------
PSO [27]44,740------
DE [28]44,526.1----200
MDE [28]42,611.14----125
CSR [15]42,440----109
TLBO [10]42,385.8842,407.2342,441.36--
HDE [28]42,337.3----48
SQP [17]42,120.02----625.07
KHA [16]41,926------
MHDE [28]41,856.5----31
CSA [25]41,046.897----94.4
Proposed RCSA41,013.0941,401.541,789.917
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheyuan, C.; Hammid, A.T.; Kareem, A.N.; Jiang, M.; Mohammed, M.N.; Kumar, N.M. A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem. Sustainability 2021, 13, 4277. https://doi.org/10.3390/su13084277

AMA Style

Zheyuan C, Hammid AT, Kareem AN, Jiang M, Mohammed MN, Kumar NM. A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem. Sustainability. 2021; 13(8):4277. https://doi.org/10.3390/su13084277

Chicago/Turabian Style

Zheyuan, Cui, Ali Thaeer Hammid, Ali Noori Kareem, Mingxin Jiang, Muamer N. Mohammed, and Nallapaneni Manoj Kumar. 2021. "A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem" Sustainability 13, no. 8: 4277. https://doi.org/10.3390/su13084277

APA Style

Zheyuan, C., Hammid, A. T., Kareem, A. N., Jiang, M., Mohammed, M. N., & Kumar, N. M. (2021). A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem. Sustainability, 13(8), 4277. https://doi.org/10.3390/su13084277

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