Applying of Marine Predators Algorithm Linked with Reservoir Simulation Model considering Sedimentation for Reservoir Operation

This paper proposes the Marine Predators algorithm (MPA) linked with the reservoir simulation model and considering sedimentation, in order to improve reservoir rule curves. The release criteria of the hedging rule (HR) and standard operating policy (SOP) were investigated in this study. The results showed that the patterns of the new optimal rule curves from the MPA technique considering sedimentation and using HR were practically useful in this study, as were the patterns of the new optimal rule curves using SOP and those of the existing rule curves. Furthermore, the new optimal rule curves using HR criteria were able to alleviate both water scarcity and excess water situations better than both existing rule curves and optimal rule curves using SOP in terms of minimal average water shortage. The new curves reduced minimal average water shortage by 53% and excess release water deficit by 19%, whereas the frequency of water shortage term was increased by 3%. The results of rule curves efficiency from MPA were higher than GA and FPA techniques in terms of providing solutions. There was a significant difference in the efficiency of water problem alleviation between considering and not considering sedimentation. It can be concluded that the MPA linked with reservoir simulation using HR criteria and considering sedimentation can be used to find optimal rule curve solutions effectively.


Introduction
Water is an essential resource for life on Earth from past to present. Nowadays, there are problems of water quantity and quality and con icts of water use. e cause of these issues is due to economic and social expansion and due to the ine cient management of water resources [1][2][3]. ese causes lead to a crisis in water supply. In addition, agriculture and industry, which provide the basis for national security, also have e ect. For these reasons, suitable water resource management is an essential issue for all areas.
Water resource management involves planning, implementing, monitoring, analyzing, and correcting water to achieve consistency with the current situation which focuses on sustainable development on economics, society, and environment to maximize bene ts consistently with available water. Generally, water resource management can be divided into construction and nonconstruction types. Nowadays, nonconstruction sites are popular as they are relatively more e ective and take a short time and are recognized.
Reservoir management is one of the most e ective measures for the development and management of integrated water resources in nonconstruction sites. e reservoir is responsible for storing, allocating, and mitigating oods and drought, which are managed under the reservoir operation process. e mission of reservoir operation is to decide how to store and deliver water for various purposes, planning how much water should be collected and delivered from the reservoir at di erent intervals through an essential and fundamental tool called reservoir rule curves [4][5][6].
Reservoir rule curves, also known as rule curves, are composed of two curves; the upper rule curve (URC) and the lower rule curve (LRC). ere are daily, monthly, and yearly curves. It can be seen that the water level control in the reservoir is most feasible between the upper and lower levels. However, the use of rule curves is long-term management and this may impair the efficiency of reservoir operations [7,8]. erefore, a wide variety of rule curves have been developed. As an alternative or decision-making tool, the rule curves used in reservoir management have release criteria for determining water discharge as an integral part of reservoir management.
Often, the critical situation occurs in reservoir operation, when the inflow to reservoir is greater than the remaining reservoir capacity during the rainy season, but the water shortage downstream still occurs in next dry season. erefore, the optimal reservoir operation is required for solving this critical problem.
ere are many reservoirs facing this problem [9,10]. For this reason, suitable release criteria and completed physical reservoir information are required in simulating reservoir operation.
A release criterion is a condition for controlling the release or storage of water for reservoir operation. Standard operating rule (SOP) is to release as much water as the reservoir can provide to meet the target delivery [11][12][13]. Linear decision-making and hedging rules (HR) are applied to address the risks and damages caused by severe water shortages in the future [14,15]. e water release criteria of the SOP have attracted attention and have been developed and implemented in a wide range of applications. However, such SOP can cause a single high range of water shortages. is is because the amount of water flowing into the reservoir has changed. erefore, to mitigate current and future impacts, the HR policy was developed for the operation of reservoirs during the dry season under different conditions, which use the principle of high discharge [16][17][18][19][20]. In order to distribute predehydration, the HR water discharge threshold was found as a way to effectively resolve the single-period water shortage and reduce the water shortage. It can alleviate drought as well. Also, annual HR release criteria in reservoirs used to allocate agricultural water needs under the effects of climate change can very well mitigate current and future drought. Moreover, HR release criteria are suitable to be used in conjunction with a reservoir level control curve for managing reservoirs with runoff volumes greater than the crisis reservoir for managing both flood and drought situations.
Reservoir sedimentation is a major issue in many parts of the world which may be exacerbated by changes in catchment land use. e reservoir is designed to have enough capacity for the required amount of water. In addition, the reservoir function is determined by considering the correlation curve in capacity-area-elevation of the reservoir. However, sedimentation in the reservoir decreases the storage capacity of the reservoir over time [21][22][23][24][25].
In the past, finding the reservoir rule curves early was a trial and error method.
is is suitable for less complex reservoir systems and is based on the experience of calculators. It is, therefore, uncertain whether it is the optimal rule curve if the reservoir system is more complex [26,27]. Later, optimization methods were applied and developed to find rule curves, e.g., using simulation, dynamic programming [28,29], genetic algorithms [30][31][32][33], genetic programming [34], Tabu search [35,36], Harris hawks optimization [37], wind driven optimization [38], firefly algorithm [39], flower pollination algorithm [40,41], grey wolf optimizer [42], and fast orthogonal search (FOS) [43]. Optimization techniques have been developed and applied in a wide variety of applications in solving numerical and engineering problems. Because optimization can find quite global best optimal using many techniques, they are smart and can find different answers. Regardless, these optimization techniques are the ultimate optimization techniques and are inspired by evolution, certain behaviors, and randomness mimicking natural phenomena. However, there is a constant need to adopt new techniques to solve both complexity and application problems quickly and easily.
A compelling new algorithm in the metaheuristic group is the MPA inspired by the foraging of the great and intelligent sea predators [44]. e MPA has been applied to solve engineering problems [45][46][47][48][49] such as designing a spring for compression tension, welded beam, and pressure vessel. However, it is not commonly used to find the optimal rule curves. Hence, it is an interesting technique that can be applied to the reservoir simulation model for solving rule curves problem.
According to the literature study above, the MPA approach is very successful when compared to other procedures under the same conditions and it is quite valuable when applied to other issues. erefore, this research aimed to find optimal reservoir rule curves using the MPA linked with the HR release criteria considering sedimentation of the Ubolratana reservoir, Khon Kaen province in the northeast area of ailand. e results of the study were divided into three main parts: (1) the efficiency of MPA rule curves considering HR and SOP in terms of maximum water shortage, (2) the efficiency of MPA rule curves considering sedimentation using HR and SOP, and (3) the comparison of optimal rule curves of MPA, GA, and FPA algorithms and their performances in terms of water shortage situations.

Study Area.
e Ubolratana reservoir is located at longitude 102°37′06.0"E and latitude 16°46′31.4"N in Khon Kaen province in the northeast of ailand as shown in Figure 1. e normal storage capacity and dead storage capacity are 2,431 MCM (10 6 m) and 581.67 MCM, respectively. e water surface area at normal storage is 137.90 km 2 . A schematic diagram of the reservoir is presented in Figure 2 which indicates that the downstream water demands from the reservoir are electricity generation, irrigation, flood control, industrial demand, domestic water supply, and environmental conservation. e monthly water demands from the reservoir are shown in Table 1 indicating that the largest requirement is for irrigation and the least is industrial demand.

Inflow Data.
e upper watershed area of the Ubolratana reservoir is 11,960 km 2 covering three provinces of 2

Advances in Civil Engineering
Nong Bua Lamphu, Chaiyaphum, and Khon Kaen. Historic in ow data into the Ubolratana reservoir has been recorded for 52 years from 1971 to 2020; it was average when annual in ow was 2,465 MCM as shown in Figure 3. e data indicated that the maximum yearly in ow was 5,884 MCM in 1978 whereas the minimal yearly in ow was 387 MCM in 2019.

Reservoir Simulation Model.
e reservoir operation performs under the reservoir simulation model using water balance equation considering reservoir rule curves and release criteria. e operation starts from calculation of the available water using the water balance concept considering monthly in ow and water demands of downstream sites. e monthly release of water is estimated by considering the monthly available water using release criteria and reservoir rule curves. For this study, the reservoir operation model was created following the concept of the water balance. e release criteria were considered in this study to consist of the HR and the SOP which evaluated their performance. e one-point hedging rule and standard operating policy are expressed in Figure 4. e one-point hedging rule and the standard operating policy are presented in the following equations, respectively. e HR constraints are as follows:

Advances in Civil Engineering
Here, R ],τ is the total release of the aggregated reservoir at time τ; SWA τ and EWA τ are the starting and ending water availability of the aggregated reservoir at time τ; and D τ is the water demand for the water-supply system at time τ.
e SOP constraints are as follows: Here, R ] , τ is the release of water during year ] and month τ (τ �1 to 12 representing January to December); D τ is the net water demand during month τ; D t is the lower rule curve of month τ; Here, D t + C is the upper rule curve of month τ and W ] , τ is the available water by calculating the water balance concept during year v and month τ, as described in where S ] , τ is the stored water at the end of month τ; Q ] , τ is the monthly inflow to the reservoir; Eτ is the average value of the evaporation loss. e operating policy usually reserves the available water (W ],τ ) for mitigating the risk of water shortage in the future, when 0 ≤ W ],τ < D −D τ under longterm operation.

e Objective Functions for Searching Optimal Rule Curves.
e objective functions of search procedure in this study were the minimal average water shortage as described in the first following equation, the minimal frequency of water shortage is shown in the second following equation, the minimal average excess water per year is shown in the third following equation, and the minimal frequency of excess water is revealed in the final following equation. ese objective functions were used for calculations in the MPA. e minimal average water shortage per year is e minimal frequency of water shortage is e minimal average excess water per year is e minimal frequency of excess water is where H (avr) is average water shortage per year, Fre (i) is frequency of water shortage, n is the whole magnitude of examined years, and Sh V is minimization average water shortage in year v (year in which releases are less than the target demand) and Spv is the excess release water during year v (year in which releases are more than the target demand).

Application of MPA and Reservoir Simulation Model for
Searching Optimal Rule Curves. MPA is a novel metaheuristic algorithm developed to emulate the foraging strategies of the ocean predators and their interactions with the prey. MPA uses the widespread foraging strategy called the Brownian and Lévy [44]. If the concentration of prey in the hunting area is high, predators use the Brownian method and when the prey is low, they use the Lévy method. e MPA search process is divided into three phases based on different speed ratios: (1) a high-speed phase, where the prey speed is faster than the predator speed; (2) a unit speed ratio phase, where the prey speed and the predator speed are similar; and (3) a low-speed phase, where the prey speed is slower than the predator speed. In each stage, the movement of the predator and prey in nature is imitated separately [48,49]. Connect the MPA with the reservoir simulation model: (1) Start with input data and set the MPA parameters such as number of populations (N), boundaries (Xu, Xl), and total number of iterations (T). e total number of populations will take part to optimize the formulated objective functions in the search space.
(2) For this study, each decision variable represents the monthly rule curves of the reservoirs, which are defined as the upper rule curves and the lower rule curves. After the first set of fitness values of Prey Matrix in the initial population have been calculated (24 decision variables that consist of 12 values from the upper rule curves and 12 values from lower rule curves), the monthly release of water will be calculated in a reservoir simulation model using the HR and SOP considering those rule curves. (3) Next, the released water is used to evaluate the objective functions that were described in (4)- (7) of the previous section. After that, the objective functions will be used to determine the fitness value of Prey Matrix and then construct the Elite Matrix. (4) ese rule curve parameters are evaluated under three phases of MPA [44]. e process will create new rule curve values in the next iteration. is procedure is repeated until the stop criteria condition is met and the optimal 24 values of the rule curves have been obtained. e diagram flowchart of the proposed method is shown in Figure 5.

Assessment of Sediment Load in Ubolratana Reservoir.
e Ubolratana reservoir was built in 1966, with creating water surface area and capacity curves. e water surface area and storage capacity curve have been used to calculate storage and sedimentation. is curve was updated for estimation of storage capacity and sedimentation Advances in Civil Engineering accumulation using remote sensing data in 2019 and an actual survey in 2009, which calculated the decrease in capacity at each water elevation. After 10 years, the capacity of the Ubolratana reservoir had decreased [50]. is study used their data for the cases that considered sedimentation in the reservoir, as shown in Figure 6.

e E ciency of MPA Rule Curves considering HR and SOP.
e optimal rule curves using MPA linked with the reservoir simulation model considering HR and SOP are shown in Figures 7 and 8. e patterns of optimal rule curves    Advances in Civil Engineering from MPA technique considering HR (RC1-HR-Avs, RC2-HR-Fqs, RC3-HR-Exr, and RC4-HR-Fqex) in reservoir condition were higher than the patterns of MPA technique considering SOP (RC5-SOP-Avs, RC6-SOP-Fqs, RC7-SOP-Exr, and RC8-SOP-Fqex) and current rule curves (existing). Moreover, the lower rule curves from using HR criteria were higher than lower rule curves when using SOP criteria especially in the dry season (April-May). ese mean that optimal rule curves from using HR attempt to retain water by limiting the water discharge during dry season according to the concept of using HR [15,18]. It also indicates that the upper rule curves from using HR criteria are higher than upper rule curves when using SOP criteria in the end of rainy season (Oct.-Nov.). As a result, late rainy season storage capacity using HR is higher than that when using SOP and the existing rule curves for reducing severe water shortages in next dry season. is is the main purpose of applying HR criteria with rule curves for reservoir operation [14,16]. e situations of water shortage and excess release that arise from using the new rule curves generated from the MPA with the HR criteria and SOP criteria are shown in Tables 2 and 3. It is seen that the situations of water shortage when using the historic in ow under the HR with the objective functions of the minimal average water shortage rule curves (RC1-HR-Avs) were the least as 115.769 MCM/year and 742.00 MCM/year for the average water shortage and the maximum water shortage, respectively, whereas the frequency of water shortage was the highest at 0.654 times/ year as shown in Table 2. It is also clear that the situations of   Table 3.
It is concluded that the situations of water shortage and excess water when using the obtained rule curves from MPA considering HR were fewer than situations of water shortage and excess water of using the obtained rule curves from MPA considering SOP. erefore, the HR criteria control water release limitedly for saving water in order to alleviate water de cit in the next dry season [15,18]. However, the SOP criteria control release of water in order to meet target demand for all considered duration times according to many previous studies [11,30,32,33]. Hence, the SOP criteria are less suitable than HR criteria for reservoirs with high frequency of drought problems.

e E ciency of MPA Rule Curves considering Sedimentation Using HR and SOP.
e optimal rule curves from MPA technique considering sedimentation with using HR and SOP are presented in Figures 9 and 10. It can be seen that the patterns of new optimal rule curves from MPA considering sedimentation and using HR (RC9s-HR-Avs, RC10s-HR-Fq, RC11s-HR-Exr, and RC12s-HR-Fqex) were higher than optimal rule curves using SOP (RC13s-SOP-Avs, RC14s-SOP-Fq, RC15s-SOP-Exr, and RC16s-SOP-Fqex) as well as current rule curves (existing). In addition, the lower rule curves from MPA considering sedimentation using HR criteria were higher than lower rule curves of using SOP criteria for the same condition. Furthermore, the lower rule curves from MPA considering sedimentation using HR criteria were higher than lower rule curves of using SOP criteria especially in the dry season (Mar.-May) as shown in both Figures 9 and 10. e gures also indicate that the upper rule curves from MPA considering sedimentation  Advances in Civil Engineering using HR criteria were higher than upper rule curves of using SOP criteria at the end of rainy season (November). e optimal rule curves from MPA considering sedimentation using HR attempt to retain water by limiting the water discharge during the dry season (Mar.-May). Hence, the releases of water during March-May controlled by HR were smaller than releases controlled by SOP for all years.
is means that optimal rule curves from MPA considering sedimentation using HR attempt to retain water by limiting the water discharge during the dry season according to the concept of using HR [15,18]. For this reason, the storage capacity at the end of rainy season using HR considering sedimentation was higher than that when using SOP and the existing rule curves for storing water in order to reduce severe water shortage in next dry season. is is the main purpose of applying HR criteria with rule curves for reservoir operation [14,16]. e situations of water shortage and excess release when using the optimal rule curves from MPA technique and considering sedimentation and using HR and SOP are shown in Tables 4 and 5. ey indicate that the circumstances of water shortage when evaluated by reservoir simulation under historic inflow using rule curves of HR considering sedimentation and objective functions of the minimal average water shortage (RC9s-HR-Avs) were the least at 95.558 MCM/year and 693.000 MCM/year for the average water shortage and the maximum water shortage respectively, whereas the frequency water shortage was the highest at 0.654 times/year as shown in Tables 4 and 5. e tables present that the situations of excess water when  Advances in Civil Engineering evaluated by means of reservoir simulation under historic in ow using rule curves of HR considering sedimentation and objective functions of the minimal average excess water (RC9s-HR-Avs) were the least at 1,087.807 MCM/year and 4,105.658 MCM/year for the average excess water and the maximum excess water, respectively. It is concluded that the situations of water shortage and excess water when using the obtained rule curves from MPA considering sedimentation and using HR were smaller than situations of water shortage and excess water of using the obtained rule curves from MPA considering sedimentation and using SOP. e HR criteria control water release limiting and saving water in order to alleviate water de cit in next dry season [15,18]. erefore, the rule curves from the MPA considering sedimentation with the HR can be applied to reduce the risk of unacceptably large damage from water shortage during the dry season.
In addition, the results from Tables 4 and 5 also indicate that the situations of water shortage and water excess were quite di erent, when sedimentation was or was not considered. e considering and not considering sedimentation cases were evaluated via the water surface area and the capacity curve [49]. e sediment accumulation from rainfall owing into reservoir resulted in a decrease in reservoir capacity, so it is necessary to consider the reservoir's long-term sediment accumulation. erefore, longterm reservoir operation considering sediment accumulation is an important operational parameter in order to ensure future accuracy and sustainability [21,23,24]. It can be concluded that the MPA linked with reservoir simulation using HR criteria and considering sedimentation can be used to e ectively nd optimal rule curves solution.

Comparison of Optimal Rule Curves Performance of MPA, GA, and FPA.
e optimal rule curves from MPA, GA, and FPA techniques linked with reservoir simulation model considering sedimentation and using HR were plotted in Figure 11. ey indicate that the patterns from the new rule curves obtained from the MPA, the GA, and the FPA are similar because of the seasonal in ow e ect and the same conditions. e results also show that the upper rule curves of all techniques considering sedimentation (RC17-MPAs-HR, the RC18-GAs-HR, and the RC19-FPAs-HR for MPA, GA, and FPA, resp.) were higher than the existing upper rule curves. ese patterns can promote reduction of spill water and maintain full storage capacity full at the end of the rainy season. is will help prevent water shortages in the following dry season. However, lower rule curves of all techniques considering sedimentation during dry season (Jan.-May) were higher than the existing upper rule curves. ey can control water discharge by reducing water release at lower levels than target demand according to the HR concept [14][15][16]18]. e optimal rule curves from MPA, GA, and FPA techniques linked with reservoir simulation model considering sedimentation and using HR were used to evaluate the performance of mitigation water shortage and water excess situations by reservoir simulation considering historic inflow of 52 years; the results are shown in Table 6. ey indicated that the situations of water shortage and water excess when using optimal rule curves from MPA, GA, and FPA techniques were slightly different because their patterns were similar and they had the same conditions. It can be concluded that the MPA linked with reservoir simulation model considering sedimentation and using HR can be used to find optimal rule curve solution effectively like GA and FPA techniques. e efficacy of searching for rule curves solutions from MPA, GA, and FPA techniques was investigated by comparison of iteration number for all techniques. e results of searching iteration number are present in Figure 12 indicating that the optimal rule curves of MPA were obtained at 350 iteration number whereas the iteration number for GA and FPA techniques was 630 and 450, respectively.
It can be concluded that MPA technique had higher performance than GA and FPA techniques in reservoir rule curves searching. It may be inferred that, like the GA and FPA processes, the MPA methodology is effective in locating reservoir rule curves. However, while the outcomes are equivalent, the speed of search or the complexity of the system is also an essential factor, which MPA can handle better than other strategies.

Conclusions
is paper proposes the Marine Predators algorithm (MPA) linked to the reservoir simulation model considering hedging rule criteria and sedimentation in the reservoir for improving reservoir rule curves. e Ubolratana reservoir, located in Khon Kaen province, ailand, was considered for this study. e release criteria of the standard operating rule (SOP) and hedging rule (HR) were applied to solve rule curves solution. e e ciencies of rule curves from three techniques (MPA, GA, and FPA) were evaluated. e results revealed that the proposed model with four objective functions provided the new optimal rule curves. e patterns of optimal rule curves from MPA technique considering sedimentation in reservoir condition were higher than the patterns of existing rule curves for all other cases. e situations of water shortage from using optimal rule curves of HR criteria in terms of frequency were higher than those when using SOP criteria, whereas the average water shortage term of using HR criteria was less than that when using SOP criteria. is is the main objective of using HR criteria for determining release conditions. e results also showed that the optimal rule curves from considering sedimentation in the reservoir were more reasonable simulations than those not considering all cases of using historic in ow samples. In addition, the situations of water shortage and water excess were quite di erent from those concerning   considering and not considering sedimentation cases. It can be concluded that sediment accumulation from rainfall flowing into reservoir must be considered to understand the reservoir's long-term sediment accumulation. e results of the comparison for the rule curves search efficiency of the MPA technique and GA and FPA techniques showed that the optimal rule curves of MPA technique were similar to optimal rule curves of GA and FPA techniques. It can be concluded that the MPA with HR considering sedimentation can be used to find optimal reservoir rule curves effectively with both mitigating flood and drought situations. In addition, the MPA technique is faster in producing optimal rule curves compared with GA and FPA techniques.

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

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.