Reservoir rule curves are essential rules for store activity. This investigation connected the Genetic Algorithm, Firefly Algorithm, Bat Algorithm, Flower Pollination Algorithm, and Tabu Search Algorithm associated with the store reproduction model to look through the ideal supply standard bends, utilizing the Huay Ling Jone and Huay Sabag supplies situated in Yasothorn Province, Thailand, as the contextual investigation. Memorable inflow information of the two repositories were utilized in this investigation, and 1,000 examples of engineered inflows of stores were utilized to recreate the repository activity framework for assessing the acquired principle bends as displayed as far as water circumstances. Circumstances of water lack and abundance water appeared as far as the recurrence extent and length. The outcomes demonstrated that GA, FA, BA, FPA, and TS associated with the reservoir simulation model could give the ideal principle bends which better moderate the drought and flood circumstances contrasted and current guideline bends.

1. Introduction

Flood and drought situations are expanding in numerous districts in view of high populace development, monetary extension, land use change, and environmental change; therefore, to avoid the circumstances described above, water resource management must be improved. Similarly, the demand and supply management sites of water resource management are considered a solution to these serious problems; however, in order to save time and money, a nonconstruction method is required; for example, to escalate the effectiveness of irrigation, to engage public participation, and to reform reservoir operation management. In particular, reservoir operation was based on reservoir rule curves. This can be considered a widely used practice for the management and allocation of water resources from reservoirs [1]. It is necessary as guidelines for water release and storage through the long-term operation conducted monthly under the inflow conditions for the reservoir and downstream demand criteria [24].

In general, water allocation from a reservoir can be considered a complex and nonlinear problem due to a large number of conditions and variables. However, from the past ten years, improvements have been made to reservoir rule curves to address complex problems using the optimization technique, which is widely studied [5, 6], because it provides quantitative results consistent with the hydrological change situation, including its application in water allocation to reduce economic losses caused by pollution [7, 8]. It can also be used as a tool for decision support systems (DSS) in finding the optimal land use in relation to groundwater vulnerability [9].

At first, numerous improvement systems have been connected to locate the ideal reservoir rule curves; for instance, Dynamic Programming (DP), Simulated Annealing Algorithm (SA), Particle Swam Optimization (PSO), Cuckoo Search (CS), and Tabu Search (TS) [1023]. Furthermore, those obtained principles’ rule curves were viably connected by looking through conditions, for example, historic inflow data, water request in downstream, standard working, and smoothing capacity rule just as the target elements of looking through methodology. In any case, most past strategies have been connected to locate the ideal reservoir rule curves of just the enormous repositories in Thailand (storage volume higher than 100 MCM). These days, there are numerous new viable improvement methods that have not yet been connected to locate the ideal reservoir rule curves with little supply. The mentioned techniques include Genetic Algorithm (GA), Firefly Algorithm (FA), Bat Algorithm (BA), Flower Pollination Algorithm (FPA) and Tabu Search (TS) that are effective alternative optimization techniques [2436]. These optimization techniques have evolved from an evolutionary algorithm into the metaheuristic optimization algorithm and swarm intelligence, which are inspired by biological evolution and natural systems. For example, GA uses the concept of solution operators such as reproduction and mutation. FA has a concept from the firefly’s flashing behavior to attract other fireflies to fly around them and form small groups. The concept of BA mimics the behavior of bats that use sonar for navigation and echolocation. Similarly, FPA draws inspiration from self-pollination and cross pollination by other plants. In addition, for the TS algorithm, a data structure is used to remember solved problems and know the accepted answers, where the model will not resolve the same problem again or the problem has a limited time to resolve.

This study thus proposed the GA, FA, BA, FPA, and TS associated with the reservoir simulation model to find the optimal reservoir rule curves of small reservoirs in Thailand. The Huay Ling Jone and Huay Sabag reservoirs in Yasothorn Province were selected as the case study.

2. Methodology

2.1. Reservoir Operation Model

The system of reservoir operation typically comprises of the accessible water regulated by utilizing the water balance ideas and water requests from downstream; then, the month-to-month discharged water is evaluated by thinking about the month-to-month accessible water with discharge criteria, working arrangements, and repository guideline bends. For this examination, the store task model was made after the idea of the water equalization to be worked under the standard working strategy as communicated in Figure 1 and following equation [37]:where Wτ is the accessible water between month τ (τ = 1 to 12 represents January to December), Sτ 1 is the maintained water at the end of month t1, Qτ is the reservoir’s monthly inflow, Pτ is the precipitation over the reservoir, Eτ is the average value of the evaporation loss, Sτ is the maintained water at the end of month τ, Rτ is the water release during month τ, and Dτ is the water demand from the reservoir during month τ. Released water from the reservoir simulation model (Rτ) was used for a calculation of the water shortage and excessive water release circumstances in which these circumstances were used for calculating the objective function of the searching process (average annual shortage per year).

2.2. The Objective Functions for Searching Optimal Rule Curves

Currently, the hydrological system has been impacted by climate change, inflow to reservoirs, and water allocation to downstream areas differently than in the past. For example, in 2012-2013, there was a severe drought throughout the country, while in 2016, repeated floods caused great damage, especially in the northeastern region of Thailand. Therefore, water resource management of reservoirs in Thailand focuses on two main tasks: (i) providing sufficient water during the drought period and (ii) preventing flooding during the rainy season. Consequently, the evolution of optimal rule curves for reservoirs in this study would like to focus on two objective functions: (1) minimum water shortage, as shown in equations (3) and (4), and (2) minimum water exceeding release, as shown in equations (5) and (6) as follows:where Z is the minimum average volume of insufficient water per year, n is the whole magnitude of examined years, is the magnitude of water shortage during years ν (year in which overflows do not meet 100% of the objective condition), and i is the iteration rounds. Later, the finding system is extended on the circumstance that the completion criterion is conceded. This completion norm is the optimum objective function with a negligible change (slightly less than 0.10 MCM).

The minimum average volume of release excess water per year () is located as the second objective function of the considering procedure subject to the restrictions in the simulation model as illustrated in the following equations:where n is the total magnitude of examined years, is the volume of discharge excess water during years ν (year in which discharge excess water are higher than the objective condition), and i is the iteration rounds.

2.3. Applying GA, FA, BA, FPA, and TS for Searching Optimal Rule Curves

The main idea of the paper was to find the optimal monthly storage volume of the reservoir in relation to input factors (e.g., inflow, precipitation, and hydrologic conditions) and output factors (e.g., water allocation to downstream areas) over a 12-month period in a year. The search for the most appropriate response is based on the optimization techniques of each type, integrating the link with the reservoir simulation model, with the aim of comparing the results, finding the strengths, observations of each technique for the development of optimal reservoir rule curves, and their applications. They demonstrated that the choice factors (ideal standard bends, 24 choice factors comprising of 12 esteems from the upper bound and another 12 esteems from lower bound circumstances) of each inquiry were made from the beginning advance. At that point, these factors were utilized in the store recreation model for activities of the framework under the conditions relating to every repository, for example, recorded inflow information, water request information, and other physical trademark supply. In the wake of running the model, the month-to-month discharge waters were acquired, and these were utilized to compute the target capacity of the scanning process for assessing and creating new choice factors. Next, the supply reproduction model was run again until the stop condition was met. This system was rehashed until 24 estimations of the standard bends were proper. The utilization of these streamlining strategies (GA, FA, BA, FPA, and TS) to associate with the repository recreation model for looking through ideal guideline bends is descripted as follows.

2.3.1. GA Optimization for Searching Rule Curves

In the GA technique used to find for optimal rule curves, it relies on an encoding scheme to convert a decision variable into a chromosome, performing selection, crossover, and mutations [38]. After this procedure, the genetic action generates a new chromosome. Each decision variable is defined as an average monthly water retention value, which is defined as an upper margin and lower margin, subsequently calculating the first set of chromosomes in the initial population (24 decision variables comprising 12 values from the upper margin and 12 values from the lower margin situation). The amount of released water is recomputed by the reservoir model associated with the GA algorithm. Then, the discharged water is used to define the objective function to evaluate the suitability of GA. Through the reproduction process, the point values of the control curve are regenerated in the next generation. At this stage, the calculation is duplicated until the value of 24 points of rule curves’ answers is appropriate. The GA associated with the reservoir simulation model can be illustrated, as shown in Figure 2.

2.3.2. FA Optimization for Searching Rule Curves

The principle FA is based on the communication behavior of individual firefly populations that seek to find the light through the attraction of the brightest fireflies [39]. The intensity of the emitted light determines the attraction. This is usually determined by the fitness value obtained from the objective function (minimum water shortage and minimal overflowing, respectively). Figure 3 illustrates the process of finding the optimal reservoir rule curves using the FA technique. It can be explained that the initial part defines the firefly population boundaries from the distance and dimension between two fireflies (Xi and Xj and xid and xjd, respectively) [40] as 12 upper bound points and 12 lower bound points, which totaled 24 points. Then, the reservoir simulation model was associated to the population conditions with the required objective function, and a computational process was started to find out. Next, the results from the first set of computations (the best 12 rule curve points in this computational cycle) are considered as a more concentrated light that attracted the new firefly population to converge towards the response region of the new rule curves’ point series. The answer obtained from each calculation is the optimal rule curve point, which is updated until the required number of rounds of calculations are completed and satisfied with the answer.

2.3.3. BA Optimization for Searching Rule Curves

To compute the optimal rule curves using the BA technique, there is a first step consisting of three condition definitions: (i) defining a sample of the bat population; here, the possible amount of retained water between the upper and lower bounds in each month, (ii) determining the initial response set derived from the echolocation of the bats to the target to reduce the search distance as much as possible to the default, and (iii) in the search-answer phase, an algorithm was assigned to mimic the behavior of the bat’s random flight at speed () at (xi) position with different frequencies from minimum (fmin) to maximum (fmax) or using different wavelengths (λ) and loudness (A), and the wavelength (or frequency) of the emitted pulses can be adjusted automatically. In addition, the impulse delivery rate (r) is related to the target proximity range [41], which is the best rule curves’ point of 12 months in this round of calculations. Then, the first calculated answer set is considered in the context of the reservoir simulation model process and the objective function. However, if the answer in each round does not meet the criteria of the objective function, the algorithm process returns to the initial stage and finds and updates the answer until satisfactory and consistent results are obtained. As described previously, a schematic for the application of BA techniques to find the optimal rule curves is shown in Figure 4.

2.3.4. FPA Optimization for Searching Rule Curves

The concept of FPA was inspired by the nature of flower pollination, which has four rules of conduct: (i) global pollination consists of live pollination and cross pollination in which the pollen has a movement pattern that follows Levy flight, (ii) local pollination is considered as self-pollination and abiotic pollination, (iii) the stability of flowers by pollinators is equivalent to the reproduction probability, which is proportional to the similarity of the two flowers involved in the pollination process, and, (iv) the interaction of global and local pollination can be controlled by the alternation probability [42].

Figure 5 shows a process diagram of applying FPA to find the optimal rule curves for a reservoir. Consistent with the previous explanation, the first process is to determine the initial population based on the physical factor data of the reservoir (e.g., reservoir size, water surface area, precipitation, amount of water from the hydrologic system), determining the maximum number of rounds required for the calculations. This includes the extent of the search for the answer (the upper and lower bounds of each month’s optimal rule curve point) and the initial population generated by the random solution. Then, the algorithm explores the answer corresponding to the random solution and selects the best results for this population. Next, the results obtained from the previous steps are treated as a path to the computations to find answers in the form of global pollination with the directional movement of Levy Flight. The answers obtained from this round of calculations for the best control curve points for each month were linked to the reservoir simulation model. Then, a process is performed to evaluate the answer obtained according to the constraints with the created objective function. Finally, the algorithm terminates the calculation only when it obtains the optimal rule curves result with satisfactory evaluation results.

2.3.5. TS Optimization for Searching Rule Curves

The following steps describe how to use the TS algorithm in conjunction with the reservoir simulation model. The TS generates a random preliminary set (x1, x2,…, xn) within the free space as the possible value between the dead maintain capacity and the usual water level of the reservoir under consideration in the first phase. The potential distance is the difference between the dead maintain capacity and the ordinary water level of the reservoir. The method then considers 24 decision factors representing the 12 values from the upper margin and 12 values from the lower margin, respectively, and the potential outcomes of ith repetitions, as xi = [xi1, xi2, …, xi24]T. In this round of computation, the processed and agreed upon answer set is stored in the tabu list [43]. To complete the process, the newly created rule curve solution set is connected to the reservoir model. Finally, the 24 points of rule curves’ solution sets for each round of computation are connected to the monthly reservoir water released from the reservoir model. Figure 6 shows the algorithm TS in conjunction with the reservoir model.

2.4. Illustrative Application of Proposed Models

In this study, we selected the irrigation project in Yasothon Province, located in the northeast of Thailand to represent the study area. For the representative of the reservoirs in the study area, two reservoirs were selected, the Hauy Ling Jone Reservoir and the Hauy Sabag Reservoir, a small reservoir (collected less than 100 MCM) in the Moon River basin’s tributaries. The study area and schematic diagram are described in Figures 7 and 8. Geographically, the summary area of the Huay Sabag is approximately 52.00 km2 with a total cumulative annual average runoff of about 22.4 MCM during 1996 to 2017 (22 years), and its normal storage capacity is 30.03 MCM.

The Huay Ling Jone Reservoir has an irrigation area of 19.88 km2, full storage capacity of 21.06 MCM, and dead storage capacity of 0.40 MCM. The historic inflow data records of the Huay Ling Jone Reservoir were recorded from 1994 to 2017 (24 years). In essence, the demand for downstream areas includes consumption, irrigation, flood prevention, water supply, and ecological conservation. The historic inflow data records of both reservoirs are presented in Figure 9.

3. Results and Discussion

3.1. Optimal Rule Curves

The month-to-month noteworthy inflow information from the most recent 24 years (1994–2017), month-to-month information of dissipation, precipitation, water system, household and modern water requests, and month-to-month ecological protection of the Huay Ling Jone were imported for handling in the GA, FA, BA, FPA, and TA methods associated with the repository reenactment model, and after each one of those models were run, the ideal guideline bends were at long last acquired. The information from Huay Sabag repository were additionally brought into the proposed models before they were kept running along these lines as the Huay Ling Jone supply to get the bends. The ideal guideline bends of the Huay Ling Jone repository and the Huay Sabag store now are exhibited in Figures 10 and 11 separately. They demonstrate that the examples of all standard bends were normally comparable because of the inflow impact and the equivalent looking conditions [44]. Notwithstanding, the lower guideline bends existing in the Huay Ling Jone supply were not exactly the bends of GA, FA, BA, FPA, and TA principle bends, particularly during May to August. Although the acquired lower standard bends of GA, FA, BA, FPA, and TA from the Huay Sabag store were higher than the current principle bends, these standard bends could be utilized for an assessment of the model’s effectiveness on the reservoir.

3.2. Efficiency of the Obtained Rule Curves

Tables 1 and 2 demonstrate the circumstances of water deficiency and too much discharged water of the Huay Ling Jone repository utilizing the notable inflow of the most recent 24 years and the engineered inflow of 1,000 examples individually. They demonstrate that the normal water lack and the recurrence of water deficiency by utilizing RC3-FA principle bends were the lowest (1.917 MCM/year and 2.217 ± 0.932 MCM/year) as compared with the other standard bends that considered both historic inflow and synthetic inflow cases.. These were corresponding to the objective function of the searching. Moreover, the rule curves of RC3-FA could mitigate excessive water better than other rule curves in both of those two scenarios (6.214 MCM/year and 7.299 ± 1.810 MCM/year).

Tables 3 and 4 present the situations of water shortage and excessively released water of the Huay Sabag reservoir using the historical inflow of the last 22 years and the synthetic inflow of 1,000 samples, respectively. They demonstrate that the average water shortages of all rule curves were zero in both of the historical inflow and synthetic inflow cases. Thus, these situations due to the total inflow were too much demanding than the total water demand. Additionally, the study results expressed that the situation of excessive water by using the rule curves from GA, FA, BA, FPA, and TA techniques was lower than their existing rule curves in both of historic inflow and synthetic inflow scenarios.

3.3. Efficiency of the Obtained Rule Curves of the Huay Sabag Reservoir under 30% Increased Water Demand

The productivity of the acquired standard bends from the Hual Sabag repository was assessed again by 30% expanded water request in the downstream sites. Similar states of the repository reproduction models were utilized to recreate the supply framework with the goal that these assessments were performed in both noteworthy inflow situation and manufactured inflow situations, and its last outcomes are displayed in Tables 5 and 6.

The outcomes demonstrated that the water deficiency circumstances at Huay Sabag utilizing the current standard bends were not as much as utilizing the acquired guideline bends of GA, FA, BA, FPA, and TA, as the complete inflow was higher more than the all requirement. Moreover, over the top water circumstance utilizing the standard bends structure GA, FA, BA, FPA, and TA strategies was still lower than utilizing the current guideline bends in both historic inflow and synthetic inflow situations.

From the previously mentioned, the GA, FA, BA, FPA, and TA methods could be compelling and material for finding the ideal principle bends for the supply of an extreme water circumstance.

4. Conclusion

This examination proposed the improvement systems associated with the store reenactment model for the inquiry of the ideal guideline bends in the little supply in Thailand. The examination result confirms that the examples of these standard bends were comparable since the regular inflow impact and looking through conditions were equivalent. The examination additionally discovered that the obtained standard bends from the proposed GA, FA, BA, FPA, and TS procedures at the Huay Ling Jone store can relieve the water lack and extreme water circumstances superior to anything utilizing the current principle bends in the two situations. It tends to be presumed that the proposed procedures are reasonable for looking through the ideal guideline bends of the Huay Ling Jone repository where the complete inflow and the all-out-interest are near each other. Moreover, the acquired standard bends of the Huay Sabag utilizing the proposed strategy can all more likely moderate just over the top water circumstance than the current guideline bends in the noteworthy inflow situation and manufactured inflow situation just as with 30% expansion of water request. The circumstances of water lack of the Huay Sabag utilizing the current principle bends were not as much as utilizing the acquired standard bends of GA, FA, BA, FPA, and TS.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This research project was financially supported by Mahasarakham University.