A Compact Pigeon-Inspired Optimization for Maximum Short-Term Generation Mode in Cascade Hydroelectric Power Station

: Pigeon-inspired optimization (PIO) is a new type of intelligent algorithm. It is proposed that the algorithm simulates the movement of pigeons going home. In this paper, a new pigeon herding algorithm called compact pigeon-inspired optimization (CPIO) is proposed. The challenging task for multiple algorithms is not only combining operations, but also constraining existing devices. The proposed algorithm aims to solve complex scientiﬁc and industrial problems with many data packets, including the use of classical optimization problems and the ability to ﬁnd optimal solutions in many solution spaces with limited hardware resources. A real-valued prototype vector performs probability and statistical calculations, and then generates optimal candidate solutions for CPIO optimization algorithms. The CPIO algorithm was used to evaluate a variety of continuous multi-model functions and the largest model of hydropower short-term generation. The experimental results show that the proposed algorithm is a more effective way to produce competitive results in the case of limited memory devices.


Introduction
The metaheuristic algorithm [1] has emerged as a very promising tool to solve complex optimization problems. Original pigeon-inspired optimization (OPIO) is a new type of metaheuristic search algorithm [2]. The algorithm simulates the behavior of pigeons going home. Preliminary studies indicate that it is a very promising optimization algorithm and can outperform excellent existing algorithms [3]. OPIO exploits a population of pigeons as candidate solutions by setting boundaries and optimizing the problem by moving the candidate solutions to approach the best solutions based on a given measure of quality. The general steps of the algorithm are described below.
OPIO can solve continuous solution space problems. In addition, many versions of OPIO in the literature are proposed to solve the problem of continuous and discrete solution spaces in recent years. An improved Gaussian pigeon inspired optimization algorithm preserves the diversity of early evolution to avoid premature convergence. The entire algorithm shows excellent performance in global optimization and is effective for solving multimodal and non-convex problems with higher dimensions. Multi-objective pigeon-inspired optimization (MPIO) is used for multi-objective optimization in designing the parameters of brushless direct current motors [4]. The multimodal multi-objective

Principle of Electricity Generation of Cascade Hydropower Station
The hydropower process is actually a process of energy conversion. By constructing a hydraulic structure on a natural river, concentrating the water head, and then guiding the high water to the low-position turbine through the water channel, the water energy is converted into rotational mechanical energy, and the generator coaxial with the turbine is used to generate electricity. It is pivotal to note the conversion from water energy to electricity. The electricity generated by the generator is sent to the user through the transmission line to form the entire process of for generating electricity, as listed in Figure 1. The water body in the high-altitude reservoir has a large potential energy. When the water body flows into the downstream of the hydropower station through the hydraulic pipe installed in the hydropower station, the water flow drives the runner of the water turbine to rotate, so that the hydrodynamic energy is converted into the rotating mechanical energy. The turbine drives the coaxial generator rotor to cut the magnetic lines of force, and generates an induced electromotive force on the stator winding of the generator. When the stator winding is connected to the external circuit, the generator supplies power to the outside. This way, the selected mechanical energy of the turbine is converted into electrical energy by the generator. Water energy resources are potential energy and kinetic energy existing in rivers and are a part of renewable resources, but the reserves of water energy are related to factors such as river flow, evaporation, precipitation, etc. Rivers vary greatly from region to region and climate varies. There is also a large difference in the amount of water energy resources in the region. Due to the current technical conditions, the water volume and the drop of the river part will not be utilized, and the mutual transition between energy also has a certain loss. Therefore, the technically developable hydropower resources are usually lower than their theoretical reserves. Taking the amount of technology developable resources as the basis, the economically available hydropower resources obtained by considering factors such as transmission distance, cost, and flooding loss are smaller than the technically exploitable amount.
Hydropower station reservoirs are generally divided into two categories, one is conventional scheduling and the other is optimized scheduling. Conventional scheduling is a commonly used scheduling method. It can also be called traditional scheduling. It is the most basic adjustment method. It only bases on historical data, no longer considers any other relevant factors, and then uses classical hydraulics and runoff regulation. The scheduling diagram and scheduling rules are used to guide the operation of the reservoir, and the water level of the reservoir is calculated and then expressed as the objective function, that is, the most basic scheduling method is used to ensure the operation of the hydropower station. These schedules are drawn based on past hydrological data and tasks when the reservoir is at different water storage levels, that is, the reservoir scheduling rules when the reservoir is in different states.
The optimal dispatching model of the hydropower station reservoir is based on the optimal theory to establish a mathematical model based on the actual conditions of the hydropower station reservoir, and then using modern computer technology to find the optimal scheduling method that meets the scheduling principle in the process of establishing the optimal scheduling model of the hydropower station. In addition, there are many factors to be considered, including the connection between water and electricity as well as numerous constraints on cascade water inventory, This will minimize pollution in terms of ecological environment, and maximize profits and social benefits in terms of economic benefits. In order to meet the above requirements to the maximum extent and to minimize the water abandonment of hydropower stations, this paper establishes a short-term power generation model suitable for the optimal operation of cascade hydropower stations.

Maximum Short-Term Generation Model
With the completion and operation of a large number of hydropower stations, the effective solution of large-scale reservoir group optimization scheduling models has become an urgent problem to be solved. The reservoir optimization problem proposes a higher solution quality and running speed for the meta heuristic algorithm. Therefore, this paper proposes a new algorithm called CPIO. This aspect of research not only improves the speed of the operation, but also ensures that the quality of the solution is not worse than the original algorithm and can suppress the premature phenomenon.
Under the given inner diameter flow of the control period, the objective function of the long-term optimal scheduling model of the cascade hydropower station group is defined as: The maximum amount of cascade power generation under the condition of ensuring the output of the cascade is considered during the control period. Under the premise of satisfying the actual situation, this paper selects the maximum benefit of cascade power generation as one of the objective functions. At the same time, it is necessary in the medium and long-term optimization scheduling, it is necessary to consider the output of the period with the least output during the year as much as possible. Medium and long-term optimized dispatching provides the largest possible uniform and reliable output for the power grid, giving full play to the capacity benefits of hydropower generation which can replace thermal power.
In the formula, A, q, ∆h, ∆t represents the output coefficient, the outflow rate, the upstream and downstream water level difference respectively. In addition, unit time E is the maximum annual power generation benefit of the cascade hydropower station, N is the total number of cascade hydropower stations, and T is within one year. Calculate the total number of time slots (T = 12).
In the process of optimizing the power generation of cascade hydropower stations, it is necessary to understand the water reservoir data in order to first calculate the reservoir upstream capacity and the outflow flow value, and then calculate the downstream water level value based on the outflow flow value.
In the formula c 1 , c 2 , c 3 , c 4 , c 5 , c 6 are a set of variable constants. This set of constants has different values according to the water level in each interval, and h 1 i,j is the upstream water level value of the j-th time period of the i-th hydropower station.
In the Equation (3), q i,j is the outflow of the j-th time period of the i-th hydropower station, h o j is the initial flow of the j-th hydropower station, and t j is the number of hours of the j-th time period.
In this Equation (4), c 7 , c 8 , c 9 , c 10 , c 11 , c 12 are a set of variable constants. This set of constants has different values according to the water level in each interval, and h 2 i,j is the downstream water level value of the j-th time period of the i-th hydropower station.
In the Equation (5), ∆h i is the upstream and downstream water level difference of the i-th hydropower station. When the upstream and downstream water level difference is obtained, since the upstream water level is greatly affected, the average value of the upstream water level is selected for calculation in this paper.
Let us introduce the constraints of the objective function: The level of the water level needs to be limited between h min and h max .
The outflow of the reservoir must fluctuate between q min and q max . In order to ensure the stable operation of the power generation of the entire cascade hydropower station, the output of the power has to be relatively stable, so the outflow of the reservoir cannot be lower than the minimum flow. In addition, in order to stabilize the life of the turbine and generator, the outflow of the reservoir cannot be higher than the maximum flow.
The capacity of the reservoir should fluctuate between the V min and V max . In order to ensure that the reservoir will continue to work under special circumstances, the reservoir's capacity cannot be lower than the originally set value to ensure the safe operation and that the downstream organisms are safe, so the capacity cannot be higher than the maximum reservoir capacity.

Pigeon-Inspired Optimization
Without prejudice, the minimization problem of the objective function f (x) is discussed in this paper, where x is the vector that defines the n design variables in the domain D in the decision space.
The pigeon-inspired optimization is a meta-heuristic algorithm that is inspired by the behavior of the pigeons returning home and is widely used in most continuous or discrete optimization problems. This article mainly introduces continuity problems. Referring to extensive literature reviews, a group of pigeons move in decision space D according to the update rules to find the optimal value when looking for the solution of the problem. More formally, in order to gain the satisfactory value of the objective function f (x), the population of the pigeons is randomly sprinkled in the previously set search space. The objective function judges the equivalent quality of solution based on the position information of each pigeon. At any stage t, the i-th pigeon has its own position vector x t k and velocity vector v t k . For each pigeon, the best solution is the value of the objective function. The best position of the position where the pigeon has passed will be stored. The global optimal solution is continuously updated. To transition from the t step to the t + 1 step, a more competitive solution will be taken, and each particle is perturbed according to the following formula: and: As the formula above suggests, x t k refers to the current position of the k-th pigeon, and x gbest is the best position ever found in the entire herd, and the vector v t k is a perturbation vector, namely velocity. Finally, φ 1 is a variable constant, is a variable amount limited to 0-1, and φ 2 , φ 3 are two weight factors can be constants or variables. This stage belongs to the map and the compass operator. When the pigeon approaches the destination, the dependence on the sun and the magnetic object is reduced, and then the landmark operator is entered.
From here on, the landmark operator is entered. In this operator, the pigeons continue to iterate according to the pigeons or landmarks of the roads understood by the population. In the above formula, the purpose of this operation is to find out the pigeons with a high fitness value in the flock. This pigeon is then considered to be the pigeon that knows the road, and the pigeons are iterated according to the pigeon. N t is the population number at the t-th iteration, and F(X t k ) is the fitness function value of the k-th pigeon position.
The significance of this operation is to halve the pigeons and discard the pigeons that do not have the way to know, to prevent such pigeons from misleading the population into local optimum.
In the formula, φ 4 is a variable constant that value is a randomly generated value from (0, 1). In this operation, all pigeons that do not know the road will be iterated according to the pigeons that know the road.
For maximizing the problem, OPIO uses Equation (14) to calculate the value of F(x t k ) to find the pigeon with the ability to identify the function. For the minimization problem, OPIO uses Equation (15) to calculate the value of F(x t k ) to find the pigeon with the function of identifying. In Equation (14), χ is a non-zero constant whose purpose is to prevent the denominator from being zero.

Compact Pigeon-Inspired Optimization
The compact approach replicates the operation of the population-based algorithm by building the probability of a total solution. The optimal process encodes the probability representation of the actual population as a virtual counterpart. Compact pigeon-inspired optimization is a model built on a pigeon-inspired optimization-based framework. In the OPIO algorithm, the concept and design of the CPIO algorithm will be explored in more detail.
The purpose of CPIO is to simulate the operation of OPIO underlying overall algorithm in a smaller version of memory variable memory. By constructing a distributed data structure, the actual solution of the OPIO is transformed into a compact algorithm, the perturbation vector. The PV vector is a probabilistic model for the solution of the population.
In the formula, µ, δ are two parameters of the standard deviation and the average of the vector PV, and t is the current number of iterations. The value of µ, δ is limited to probability density functions (PDF) [38] and is changed within [-1, 1]. The magnitude of the PDF is normalized by keeping the area to 1, because by obtaining approximately sufficient in well it is the uniform distribution with a full shape.
The initialization of the virtual population is performed as follows. For each design variable µ = 0 and δ = λ , where λ is a large constant (λ = 10). This value is initialized to initially obtain a normal distribution of truncated wide shapes. The sampling mechanism of the design variable x t k associated with the generic candidate solution x in PV is not a simple process and requires extensive interpretation. For each design variable indexed by k, a truncated Gaussian PDF with the mean µ and standard deviation δ is associated, The PDF is described by the following formula: ))) (17) PDF is the probability distribution function of PV, and a truncated Gaussian PDF-related µ, and δ are formulated. A new candidate solution is generated by iteratively biasing towards a promising region of the optimal solution. Every component of the probability vector may be acquired by learning the previous generations. er f is the error function established by [39]. PDF corresponds to the cumulative distribution function (CDF) by constructing a Chebyshev polynomial [39], and the upper domain of the CDF is randomly changed between 0 and 1. CDF can be described as a real-valued random variable x with a probability distribution, and the value that can be obtained can be less than or equal to x t k . The relationship between CDF and PDF can be defined as CDF = 1 0 PDF(x)dx, and PV operations can sample the design variable x t k by randomly generating values within the range of (0, 1). In the iterative process of the compact algorithm, in order to find a better individual, a function that can be compared by two parameters is proposed in this paper. The two variable pigeon parameters are two sample individuals of the PV operation. The vector represented by the winner is the value of the fitness function. This value is higher than other virtual members, and the vector represented by the loser is that the individual fitness value is lower than the fitness evaluation standard. Two variables with return values, the winner and the loser are obtained from the calculation of the objective function, and a new candidate solution is generated to compare with the original global optimal solution to generate new winners and losers. For updating PV operations, µ and σ can be considered for updates according to the rules below. If the mean value of µ is 1, Then the update rule becomes µ t and σ t for each of its elements µ t+1 and σ t+1 [40] as described in the following: where N is the virtual population size and the value of δ is described below. The update rule for each element is given in the formula below.
Mathematical details about construction Equations (18) and (19) have been given. The persistent and non-persistent structures of rcGA have been tested and can be seen in [41]. Seeing the virtual population size N as a parameter of a compression algorithm is not a true population-based algorithm. The virtual population size, in the real-valued compression algorithm, is an algorithm that depends on the convergence speed.
In general, a probabilistic model for compact OPIO is hired to represent all of the set of solutions for the pigeon group, neither storing location information nor storing speed information; however, storing newly generated candidate solutions. Therefore, the limited storage space is required to achieve the algorithm requirements which saves a lot of time and hardware resources for the cascade hydropower station to optimize the short-term power generation model.
CPIO uses a perturbation vector PV that has the same structure as the one shown in Equation (16), at the beginning of the optimization algorithm, just like the process described in Equation (17). The PV initialization is designed as (∀k, µ[k] = 0, σ[k] = λ, λ = 10) and the variables of each design are limited to one continuous space [−1, 1], and in addition, the position x and the velocity v are randomly initialized within a certain range.
Update velocity vector and position vector by slightly revised pigeon-inspired algorithm: and: ω 1 is an inertia weight, ω 2 is a random variable between 0 and 1, and ξ 1 and ξ 2 are weighting factors that control the position update of the pigeon. It can be seen that the equation updating of speed (Equation (20)) and position (Equation (21)) is similar to the OPIO algorithm. In the original version, pigeon k was closely related to pigeon group N. In the compact version, there was no real population, but the relevance of a virtual population pigeon to the virtual population was not that great. It is easy to see that compact OPIO is just a pigeon that uses the update formula to update it, so updating it once produces a solution that saves a lot of memory.
In the landmark operator entering the second stage of CPIO, the original algorithm uses the Equation (11) to determine the pigeon with the function of identifying the function based on the fitness solution of each pigeon position, so that it becomes the center point and continues to update. Since the CPIO has only one particle to update, it is not suitable when selecting a pigeon with a path function. In this paper, a center point suitable for CPIO is proposed. By setting a virtual center position point, the guiding pigeon is updated. When the virtual center position is established, it is based on the historical fitness value of the pigeon. The number selected is also based on the size of the virtual population.
l is the number of iterations until now, N is the number of virtual populations. According to the Equation (22), the historical virtual center points of the pigeons can be selected, and it is known that they continue to iterate.
x t+1 It is easy to see that CPIO saves a lot of memory space, so this approach can be applied to other variations of OPIO.

Numerical Results
The test results of the CPIO that have been tested by 29 test functions, and these test functions come from [42]. Each test function has a very detailed introduction in Tables 1 and 2. Among these groups of questions, they have different search range and different expressions.
In Equations (20)-(23) and Algorithm 1, the parameters of the CPIO proposed herein are: N = 120, R = 0.2. The values of these parameters are referred to [43] and have a slight change. More specifically, in order to make CPIO work better, we modeled the virtual population size proposed by OPIO. In this article, CPIO is compared to the OPIO. In all test functions, CPIO is run 30 times and averaged. Take the minimum value of CPIO in all test functions. x t k = Generate from PV operation 13: // Update position and velocity 14: v t+1 end for 23: end for 24: for t = MaxDt 1 + 1 to MaxDt 1 + MaxDt 2 do 25: // Enter the second stage 26: // Select the virtual center pigeon from the historical best points 27: x t+1

Name Test Functions Range Global Minimum
Sphere Holder When initializing the two algorithms CPIO and original pigeon-inspired optimization (OPIO) , the map and compass factor R are set to 0.2, and the result is to compare CPIO and OPIO. The quality of solution and the number of runs of CPIO and OPIO optimal solutions are compared as below described. The CPIO and OPIO data results are the average of 30 runs. All algorithms operate 500 times, including 300 in the first phase and 200 in the second phase.
In Table 3, CPIO performs better than OPIO in many test functions, and most of the values perform well. In terms of the time cost comparison, it is easy to see that CPIO time spent is much better than PIO, especially in several of them, and the time spent is more than a hundred times more.
According to the comparison of the two algorithms, it can be concluded that the running time of CPIO is much lower than that of the original algorithm. This is because the number of population used in the process of iteration is different. In the new algorithm, it uses an example to keep iterating, constantly adjusting the probability distribution according to the path that has been iterated, and the greater the possibility of generating particles where the function values are superior. However, this method also has a big problem, since in the search process of a single particle, randomness is often large, and it is thus easy to fall into the local optimal. It is also relatively simple to achieve the optimal, in the case of small dimension settings, the advantages of the algorithm are not obvious. Because of this characteristic of the new algorithm, it is easy to save time and reduce the time complexity of the algorithm.

Name Test Functions Range Global Minimum
Schaffer4 Figure 2 shows the convergence trend of CPIO and OPIO. Best score obtained so far refers to the optimal value obtained by the algorithm during the iteration process. While the convergence speed of OPIO and the algebra needed to achieve optimal are small, the optimal value of CPIO is better or nearly equal to the value of OPIO. Here, CPIO uses one particle for updating and iteration, while OPIO uses the entire population for optimization. CPIO is far less than OPIO search capability, but CPIO can save a lot of memory and time to find excellence.  Among the four selected functions, Figure 3 shows the time trend of the four functions running 30 times. In general, the time spent by the CPIO and PIO algorithms does not change much, but the two algorithms compare. It is easy to see that CPIO runs much faster than the PIO. Table 4 shows the comparison of CPIO and PIO mentioned above in the memory variables, which makes it very convenient to implement the calculation algorithm. The number of variables of the two algorithms of CPIO and PIO proposed in this paper is calculated by the equation used in the computational optimization. In Table 4, it is easy to see that in the same computing situation, CPIO uses less memory than PIO. For example, during an iteration, CPIO uses an iteration Equations (16)-(23); the formula for PIO update iteration is Equations (9)-(13). Table 4. The space complexity of the two algorithms.

Algorithm Particle Memory Size Computing Complexity
Use Equations 16), (17), (18), (19), (20), (21), (22), (23)  As can be seen from Table 4, the actual population size of the PIO is N, but the actual population size in the CPIO is 1, and the virtual population number is N. In the case where the number of iterations l and the running time t are the same, the memory usage of the variables of OPIO and CPIO is iterated by 4 × t × N and 8 × t, respectively. Here, it is seen that the memory occupancy of the PIO is larger than the memory usage of the CPIO.
In Figure 4, the consequence of the presented algorithm and the else three meta-heuristics are shown. According to Table 5, the trend and optimal value of CPIO are fundamentally better than the other three algorithms, and have a superior performance. Table 5 shows the comparison of CPIO, OPIO and other algorithms, such as CPSO and PSO algorithms. Among the four meta-heuristic algorithms, the performance is as follows in 29 test functions. In the process of algorithm simulation, as part of the images are not so obvious, four relatively obvious images are extracted for display.

Experiments of Short-Term Power Generation Model for Cascade Hydropower Stations
Wanjiazhai Water Conservancy Project: The Wanjiazhai Water Conservancy Project is located in the canyon of the Tuoketuo to Longkou section of the Yellow River in the north of the Yellow River. It is the first of the eight cascades planned for the development of the middle reaches of the Yellow River. and also the Shanxi Yellow River Diversion Project. The starting point of the project the left bank is affiliated to the Pianguan County of Shanxi Province, and the right bank is subordinate to the Zhungeer Banner of Inner Mongolia Autonomous Region. The dam site controls a drainage area of 395,000 square kilometers, with a total storage capacity of 896 million cubic meters and a storage capacity of 445 million cubic meters. It has comprehensive benefits such as water supply, power generation, flood control and anti-icing.
Longkou Hydropower Station is located at the junction of two provinces, Hequ County, Shanxi Province and Zhungeer Banner, Inner Mongolia. It is 25.6 km from the upstream Wanjiazhai Water Control Project and 70 km from the downstream Tianqiao Hydropower Station. It is the regional center of energy and chemical bases in Shanxi Province and Inner Mongolia Autonomous Region, and controls the drainage area of 397, 406 square kilometers. Table 6 shows the monthly inflow values of the two cascade hydropower stations in the wet years, the flat water years and the dry years. ASP is Annual scheduling period. The short-term power generation model of cascade hydropower stations has been introduced above. Figure 5 showcases the main flow of the algorithm. In this paper, the three periods of the two cascade hydropower stations are scheduled and modeled by Equations (1)- (8) and the sum of the power generation of the two cascade hydropower stations is the largest. As shown in Figure 6, at any stage, CPIO has the largest scheduling capacity for the two cascade hydropower stations, and the total power generation is also relatively huge. CPIO dispatched the two cascade hydropower stations. The final result has the power generation at 3.968 × 10 17 KWH in the high flow year, and the total power generation at 3.108 × 10 17 KWH in the year of the median water. The power generation at 2.396 × 10 17 KWH in the low year.

Conclusions
A novel optimization approach called compact pigeon-inspired optimization (CPIO) is proposed. The proposed CPIO was tested on 29 classical test functions to demonstrate the usefulness of the proposed optimization method. A compact method is successfully used in the pigeon-inspired optimization algorithm to reduce the usage of the memory size. The proposed CPIO was also applied to cascade hydroelectric power generation. Simulation results show the CPIO may reach better results compared with some existing algorithms for the cascade hydroelectric power station.

Conflicts of Interest:
We wish to confirm that there are no known conflicts of interest and there has been no significant financial support for this work that could have influenced its outcome. We confirm that the manuscript has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed.