The Multi-Objective Optimal Scheduling of the Water–Wind–Light Complementary System Based on an Improved Pigeon Flock Algorithm

: The output of wind power and photovoltaic power is random, ﬂuctuating and intermittent, and a direct grid connection will result in the reduction of power generation income and a great ﬂuctuation in the power grid’s connection. The addition of hydropower stations can reduce the above problems. Therefore, this paper ﬁrst introduces and analyzes a typical application scenario of a water–wind–light combined power generation system. Then, a multi-objective optimization model is established, considering the two objectives of maximizing the joint generation and minimizing the system’s power ﬂuctuation. Third, the adaptive fractional order calculus strategy is introduced, and a multi-objective pigeon swarm algorithm, which can adaptively adjust the fractional order according to the location information of a ﬂock, is proposed. Finally, an optimization simulation is carried out. The simulation results show that the improved multi-objective pigeon swarm algorithm has better optimization accuracy. It provides a reference for the future implementation of hydropower stations, and the surrounding wind and photoelectric ﬁeld joint dispatching strategy.


Introduction
"Double carbon" is a strategy responding to the problems of rapid global climate change, energy scarcity and ecological and environmental damage, including efforts to reduce the consumption of fossil energy and vigorously develop new types of clean energy [1][2][3].Carbon dioxide emissions are on track to peak by 2030, and work is currently underway towards achieving carbon neutrality by 2060 [4].The proposal of "carbon peak and carbon neutral" clean energy, mainly hydropower, wind power and photovoltaic power generation, has welcomed unprecedented development opportunities.Hence, water-windlight complementary dispatching has become an important direction for the future of energy development [5][6][7].
The existing research on hydropower and wind-solar complementary power generation systems mainly focuses on capacity allocation and optimal dispatching [8].The capacity allocation of the new energy complementary power generation system is mostly a multi-objective problem due to more factors being involved, such as system stability, the economy, environmental benefits, and so on [9][10][11].At the same time, the capacity of a system often needs to consider the problem of a system optimization scheduling strategy.Optimization scheduling also typically involves multiple high-dimension decision variables; time and space coupling; the constraints of complex and various large-scale, nonlinear planning problems; and is also affected by uncertain factors such as wind, light and runoff stochasticity [12][13][14].Therefore, how to coordinate a system's capacity configuration, Energies 2023, 16, 6787 2 of 18 optimize the scheduling operation and formulate a reasonable capacity allocation scheme are the key and difficult points of such problems.
Wang et al. [15] proposed a multi-objective evolution algorithm based on Pareto optimal space, NDWA-GA and PCA for the capacity optimization allocation of a multi-energy complementary system.The correctness and effectiveness of the proposed method can be verified via a comparison with the traditional multi-objective optimization algorithm.Zhang et al. [16] put forward the three screening principles of capacity configuration to reveal techno-economic interactions.It uses the bundled system capacity configuration model of the northwest, which involves a practical engineering case exploration, and verifies the bundled system's technical and economic interactions.And it provides important theoretical support for selecting technical and economic indicators, as well as for capacity configuration.Lei et al. [17] put forward the establishment of the whole life cycle cost and unit power cost as a goal of the water-wind-light storage system capacity optimization method in the northwest region.For example, using the HOMER software to solve the optimization problem allowed the optimal capacity configuration results to be obtained.The results showed that the complementary system helped to reduce the cost of power generation and improve the power supply stability.Dong et al. [18] established investment cost, environmental protection and power supply quality as the indexes of a multi-objective optimization model.And the optimization effect of the multi-objective sparrow search algorithm was improved by introducing the Levy flight strategy to enhance the multi-objective sparrow search algorithm based on the local optimal ability.The results verify that disturbance and niche optimization technology improved the effectiveness of the multi-objective sparrow search algorithm.Yang et al. [19] proposed an optimization model of a two-layer landscape capacity configuration (including the planning layer and the operation layer) based on the dynamic regulation constraints of the system.It minimizes the investment and operation costs of the hybrid system, as well as minimizing the dynamic source and load tracking coefficient, and constructs a coordinated and optimized scheduling scheme.The results show that the proposed method, considering the above factors, has better performance in adapting to power fluctuations.Plus, it improves the capacity configuration accuracy of renewable energy, and reduces the abandoned wind volume of renewable energy.Yalcinoz et al. [20] proposed an improved multi-objective particle swarm optimization (IMOPSO) algorithm to solve the multi-objective environmental economic load scheduling problem, with a functional shrinkage factor and a functional inertial weight.It is applied to MEED problems considering wind energy, ramp limits, the embargo area, emission pollution, transmission loss and valve point effects.The test results show that IMOPSO produces better running outcomes than the comparative multiobjective optimization algorithm.Omar et al. [21] used a new algorithm based on a novel hybrid ant colony-particle swarm optimization technique, which shows that multi-target particle swarm algorithms have better performance in improving STATCOM tuning for grid code compliance.Meng et al. [22] established an optimal configuration model of wind/light/hydrogen/storage capacity, and introduced the Tent chaos mapping initialization population to improve the multi-target particle swarm algorithm.The location of leaders and followers is updated based on the adaptive spiral search strategy.The effectiveness of the algorithm is verified using the test function, and the improved multiobjective salp swarm algorithm (IMOSSA) is used to obtain the capacity configuration scheme.Kang et al. [23] established a scenery-complementary hydrogen energy storage system operation optimization model, using an adaptive simulation annealing particle swarm algorithm to solve the problem.The model can reduce the waste of wind and solar energy, and realize the system to minimize the cost of running by adjusting the power purchased from the grid, as well as the battery charge and discharge power, under the time-sharing price mechanism.Behera et al. [24] compared constriction factor-based particle swarm optimization (CFBPSO) with PSO, improved PSO and the red deer algorithm (RDA).This obtained the optimal economic and emission scheduling level, as well as an optimal solution for the renewable energy hybrid power generation system, which is composed of thermal, wind and solar generators, when the dynamic load changes within a day.He et al. [25] proposed a kind of wind-complementary power generation system, composed of heat storage, a battery, hydrogen storage and a pumped storage quantitative technical and economic comparison method.The results show that the strength Pareto evolution algorithm and the multi-target particle group optimization algorithm have good performance for capacity optimization problems.Hu et al. [26] proposed a new, multimodal, multi-objective pigeon-inspired optimization (MMOPIO) algorithm, combining self-organization mapping (SOM) and improved PIO.It establishes a good neighborhood relationship for improved PIO, and simplifies the structure of standard PIO by merging the parameters.The test results show the superiority of MMOPIO over five state-of-the-art multi-objective optimization algorithms in solving the multi-modal, multi-objective optimization problem.Wang et al. [27] proposed an improved high-level target navigation PIO (HTNPIO) and an LMS-ELS probability to balance the exploration and utilization.It compares HTNPIO with 15 other PIO and advanced heuristic algorithms, based on IEEE CEC2017 benchmark problems and five real-world optimization problems.The results show that the HTNPIO algorithm beats all competitors at CEC2017 benchmark problems, including the super-competitive LSHADE.And they are also extremely competitive in dealing with real-world problems.Most of the methods used to study the water-wind-light multi-objective capacity allocation and complementarity problems in the above literature are multi-objective particle swarm optimization algorithms, although there are problems, such as the narrow search range and low optimization accuracy.The pigeon swarm algorithm can improve the problem of the above particle swarm optimization algorithm.The memory attribute of fractional calculus can change the position state information of a current individual in a flock using the historical trajectory information in the iterative calculation process, thereby improving the optimization performance of the algorithm.At present, there is no report on the multi-objective optimization of the fractional order pigeon swarm algorithm in the water-wind-light complementary system.
The main contributions of this paper are as follows: The output characteristics of a typical water-wind-light combined system are analyzed.A multi-objective optimization model considering the maximum power generation as well as minimizing the system power fluctuation is established.According to the positional information of pigeons, the fractional order is adaptively adjusted, and an adaptive fractional order multi-objective pigeon swarm optimization algorithm is proposed.Finally, the improved multi-objective pigeon swarm algorithm, the multi-objective pigeon swarm algorithm and the multi-objective particle swarm algorithm are simulated.The effectiveness and superiority of the improved pigeon swarm algorithm are verified.

Analysis of the Water-Wind-Light Complementary Principle
Wind, light and water, as the most easily developed and utilized energy sources in renewable energy, have developed rapidly in recent years, and are not only widely used in microgrids, but the proportion of wind and photovoltaic power generation in large power grids is also increasing year by year.However, due to the random, intermittent and fluctuating characteristics of wind and photovoltaic power generation, large-scale independent grid connections often cannot meet the requirements of a large grid, and it is necessary to have a power supply with a rapid adjustment ability to compensate and adjust in order to meet the requirements of the grid connection.At present, the most widely used is the joint operation of hydroelectric power generation.The use of hydroelectric units has the ability to quickly start and stop, adjust its output and complement wind power and photovoltaic output, which can not only effectively solve the problem of new energy grid absorption but also ensure the safe and stable operation of the grid.In addition, wind power, photovoltaic power and hydropower are usually connected to the power grid through the same grid-connected connection line.The operation control strategy is unified and optimized for a hybrid system, and even has an integrated system control strategy, which shows new characteristics that are completely different from hydropower stations, photovoltaic power stations or wind farms.A power generation system operating in this way has many advantages, including better power output characteristics, maximizing the overall benefits of the integrated power generation system through the use of novel control strategies, etc.The principles of the water-wind-light integrated power generation system are shown in Figure 1.
new energy grid absorption but also ensure the safe and stable oper addition, wind power, photovoltaic power and hydropower are usua power grid through the same grid-connected connection line.Th strategy is unified and optimized for a hybrid system, and even has a control strategy, which shows new characteristics that are completel dropower stations, photovoltaic power stations or wind farms.A po tem operating in this way has many advantages, including better po teristics, maximizing the overall benefits of the integrated power through the use of novel control strategies, etc.The principles of th integrated power generation system are shown in Figure 1.

Forecast of Renewable Energy Generation Capacity
In this paper, the convolutional neural network is used to predi of the landscape [28].It uses a small number of parameters to extra between the data, and then combines them to form high-level data fea data feature extraction ability.A typical convolutional neural netwo an input layer, a hidden layer (convolutional layer, pooling layer, fu and an output layer.Its specific structure is shown in Figure 2.

Forecast of Renewable Energy Generation Capacity
In this paper, the convolutional neural network is used to predict the output power of the landscape [28].It uses a small number of parameters to extract the relationships between the data, and then combines them to form high-level data features, with a strong data feature extraction ability.A typical convolutional neural network mainly includes an input layer, a hidden layer (convolutional layer, pooling layer, full connection layer) and an output layer.Its specific structure is shown in Figure 2.
The input layer inputs the irradiance and wind speed.The convolution layer mainly performs the convolution operation, which is the operation between the feature graph of the previous layer and the convolution kernel.Then, the new feature matrix is obtained using the specific activation function, and the expression ability of the model is enhanced through the convolution layer.The mathematical expression for the convolution operation is as follows [29]: where N is the number of layers in the network; x N j is the j-th feature graph of the N-th convolutional layer; f (x) is an activation function; M is the total number of neurons in the previous layer; X N−1 i is the i-th feature diagram of the network transmitted from the The input layer inputs the irradiance and wind speed.Th performs the convolution operation, which is the operation be the previous layer and the convolution kernel.Then, the new using the specific activation function, and the expression abilit through the convolution layer.The mathematical expression tion is as follows [29]: where N is the number of layers in the network; x is the N-th convolutional layer; ( ) f x is an activation function; neurons in the previous layer; In this paper, two factors with strong correlations between wind speed and irradiance are considered in the solar power generation prediction due to the magnitude of difference in the parameters and the value difference being relatively large.Therefore, the wind speed, irradiance and historical landscape power values are normalized to map between [-1, 1].The normalization formula is as follows [29]: where x min is the minimum value of the variable and x max is the maximum value of the variable.

Constraint Condition
(1) Storage capacity constraints: where V minH and V maxH are, respectively, the maximum and minimum values of the reservoir capacity; V H t is the water storage capacity of the reservoir at the end of period t-th.
(2) Output constraints of hydropower station units: Energies 2023, 16, 6787 6 of 18 where P min h and P max h represent the minimum and maximum output of the hydropower station, respectively; P h (t) is the output of the hydropower station at the end of period t-th.
(3) Water balance constraint: where U H t is the reservoir storage capacity of the t-th period; L t is the flow of the reservoir of the t-th period; F t is the generated flow of the reservoir in the t-th period; Q t is the abandoned water quantity of the reservoir in the t period; and ∆t is a time interval.

Net Residual System Load
The difference between the load demand of the water and wind complementary system, and the photovoltaic and wind power output, is the net load of the system.The description of the overall characteristics of the net load uses the daily average net load according to the characteristics of the system net load.The average daily net load can reflect the overall level of the daily net load, according to the following formulas: where P s i is the net load of the mixed system at the i-th time; P load i is the mixed system load requirement of the i-th time; P solar i and P wind i are, respectively, the i-th time of the photovoltaic and wind output powers; and P average is the average daily net load of the hybrid system.

Objective Function
This paper involves two competing objective functions: the maximum power generation and the minimum output power fluctuation. (

1) Maximum power generation
The day is divided into 96 periods, with every 15 min being a time point; the maximum power generation of the water-wind-light complementary power generation system is the following target function: where P s t , P W t and P H t are the photoelectric output, wind power output and hydroelectricgenerated set output of the i-th time; T is the total number of intervals in the time range; M is the number of hydraulic generator sets; and ∆t is a time interval.
(2) Minimum output power fluctuation In this paper, the Wang-Mei index is used to reflect the volatility of the output power [30].This index quantifies the power output fluctuation of the system using the rotation angle and standard deviation.As shown in Figure 3, for each moment, there is a rotation angle reflecting the degree of power change at that moment.According to Energies 2023, 16, 6787 7 of 18 the Wang-Mei index, the minimum output power fluctuation of the water-wind-light complementary power generation system is taken as the following objective functions: where T is the total number of time intervals; P T t and P T t+1 is the total output power of the water-wind-light complementary power generation system at time t and time t+ 1; P T is the standard deviation of the P T ; θ t is the rotation angle of the P T t ; and k t is the gradient between the P T t and the P T t+1 .

3.1.Multi-Objective Pigeon-Inspired Optimization (MOPIO)Algorithm
The pigeon-inspired optimization (PIO) algorithm simulates the behavior-designed intelligent optimization algorithm.It designs a map

Multi-Objective Pigeon-Inspired Optimization (MOPIO) Algorithm
The pigeon-inspired optimization (PIO) algorithm simulates the pigeon homing behavior-designed intelligent optimization algorithm.It designs a map and compass operator and landmark operator phased guide on pigeons' individual flight by simulating the pigeon homing process of special navigation behavior.PIO has the characteristics of concise principles and being easy to implement.Compared with other algorithms, it has the obvious advantages of relatively simple calculation and relatively strong robustness.To solve the multi-target optimization problem, a multi-target pigeon swarm optimization algorithm based on the Pareto sorting mechanism and merging operators has emerged.At present, MOPIO has been successfully applied to solving multi-target optimization problems such as multi-UAV target searching and the fuzzy production scheduling problem [31][32][33].has emerged.At present, MOPIO has been successfully applied to solving multi-target optimization problems such as multi-UAV target searching and the fuzzy production scheduling problem [31][32][33].

Sorting Mechanism of Pareto
The ranking mechanism of Pareto consists of a non-dominant ranking operator and a crowding distance operator.First, individuals from a pigeon flock are divided into unequal non-dominance ranks based on inter-individual dominance relationships and fast non-dominance ordering.The set of non-dominated solutions for the current population consists of individuals that are located in the first rank.The crowding distance operator is reused to calculate the crowding distance between individuals in different non-dominance classes; individuals are then sorted again based on the size of the crowding distance.After the above two operators are completed, the individual pigeon is divided into different sets according to their non-dominant level, and the individuals in the set are arranged in descending order according to the size of the crowding distance.The flow of the ranking mechanism of Pareto is shown in Figure 4.

Crowd distance comparison
Save in the external archive set A The last N d was removed After the Pareto sorting is finished, the non-dominant solution set of the current population can be deposited into the external archive set.The last d N individual with the maximum non-dominance rank and minimum crowding distance will be removed.

Consolidation Operator
MOPIO uses a map operator, a compass operator and a landmark operator to obtain a merged operator, in order to update the individual speed and position of the pigeon group.The formula for the merged operator is as follows:  After the Pareto sorting is finished, the non-dominant solution set of the current population can be deposited into the external archive set.The last N d individual with the maximum non-dominance rank and minimum crowding distance will be removed.

Consolidation Operator
MOPIO uses a map operator, a compass operator and a landmark operator to obtain a merged operator, in order to update the individual speed and position of the pigeon group.The formula for the merged operator is as follows: where N p and N d respectively represent the number of populations after the i-th iteration and the number of individuals to be removed for each generation; t-max represents the maximum number of iterations; V t i represents the speed of X t i ; tr is a transition factor; X gbest is the global optimal position of the current population; and X t−1 c is the central position of the non-dominant solution set of the previous generation; the calculation formula is provided as follows: Energies 2023, 16, 6787 where S 1i is the position of the i-th individual in the previous generation of the nondominant solution concentration; n is the size of the non-dominant solution set.

Multi-Objective Fractional Pigeon-Inspired Optimization (MOFPIO) Algorithm
The Adaptive Fractional Calculus Strategy (1) Fractional calculus Fractional calculus, namely, the integral and differentiation of order as fractions, is an extension of traditional calculus and is widely used in digital image processing and digital signal processing due to its memory.There are many definitions of fractional calculus, and the commonly used G-L (Grunwald-Letnikov) definition is used here [34].Assuming the existence of continuously derivable functions of f (x) within the interval of where δ is the increment on the interval of [a, b].The n derivative of the available function is further derived from the above equation.
where n ∈ N. The υ-order fractional expression defined by G-L is as follows: where Γ(N) = ∞ 0 e −t t n−1 dt = (n − 1)! is the gamma function.According to Equation ( 17), an approximate differential expression for the υ-order integral of signal f (x) is defined by G-L as the following: (2) Pigeon flock algorithm for fractional optimization The original pigeon flock algorithm has the disadvantage of a slow convergence speed and it easily falls into the local optimum.In order to improve the global search ability, the fractional strategy is introduced into the pigeon flock algorithm.The algorithm substitutes the fractional calculus strategy into the speed update formula of the original pigeon flock algorithm.It uses the heritability and memory of the fractional calculus strategy to improve the original pigeon flock to fall into the local optimum.It also can effectively improve the convergence speed of the pigeon flock algorithm.Taking only the first four terms of Formula ( 18), the updated formula of the fractional pigeon flocking algorithm is as follows: where X t−1 inew is the location of the updated former cluster and X t i is the current location before the update.
Since the fractional order has a large effect on the algorithm, the fixed fractional order cannot be achieved, so the fractional order υ can be adaptively adjusted according to the velocity and position information of the pigeon flock.The average distance between pigeon i and the other pigeons is calculated as follows: where t is the total number of pigeons, and D is a space dimension.The evolutionary factor R is where d aver is the average of the global optimal distance of the flock from the other pigeons, the maximum mean distance between the pigeons is d max and the minimum value is d min when R∈[0, 1].When the fractional order is υ ∈ [0.5, 0.8], the algorithm works the best, so this υ can be dynamically adjusted according to the following equation:

Population Initialization
For the multi-objective problem, the power output of the unit in each scheduling cycle should be regarded as the decision variable of the individual population.Therefore, in the multi-target pigeon group algorithm, the population size is n, where each individual contains a D-dimensional decision variable, and the initialized population form is as follows: . . . . . . . . . . . .
where n is the number of generator sets; d is the number of scheduling cycles; and P ij is generated randomly between the upper and lower limits of its output.

Algorithm Process
A flowchart of the improved pigeon flock algorithm is shown in Figure 5.The process for solving multi-objective problems is shown as follows: Step 1: Initializing the population, including the population size n, position X i , speed V i , individual history optimal X p , and the number of iterations i ter .
Step 2: Initializing an external archive set A, finding the non-dominated solution set of the initial population using the Pareto sorting mechanism, and depositing it in A.
Step 3: Repeating steps 1-2 until the number of non-dominant solutions in A reaches the set point, and jumping out of the loop.
Step 4: Using the Pareto ranking mechanism to rank the external archive set A, selecting the appropriate global optimal X g for each individual separately from A, and calculating the central value X c of the initial population.
Step 5: When n ≤ 3, the speed V i and position X i are updated according to Equation ( 13).When n > 3, the fractional parameter of the current population is calculated, and the new X i (t − 1) is calculated according to Equation (19), and its speed V i and position X i are updated via substitution using Equation (13).
Step 6: Finding the non-dominant solution set of the current population and combining it into the external archive set A, sorting the external archive set A using the Pareto sorting mechanism, and keeping the first A individuals as the Pareto-optimal solution set of the current generation.
Step 7: Making n = n + 1; if n < T, return to step 4; otherwise, stop running and output the Pareto-optimal solution set.

Scenario Setting
In this paper, the landscape power generation system of a river basin in Xinjiang is taken as an example to verify and analyze the optimal dispatching model established above.The complementary system includes one hydropower station, a wind farm and a photovoltaic power station.The maximum output of the hydropower station is 720 MW and the minimum output is 60 MW.The hydropower station has a daily regulation capacity, a rated head of 68 m and a rated flow of 405 m 3 /s.The maximum capacity of the wind farm is 800 MW.The rated power of a single fan is 2 MW, and the maximum capacity of the photovoltaic power station is 1000 MW.Two typical days (a typical summer scene day and a typical winter scene day) were selected for analysis.The input and output data sizes are all 1 × 96, that is, the daily data with a time scale of 15 min.The wind power, photovoltaic power and load curves are shown in Figure 6.The correlation coefficient of the prediction results is 0.9954, the root mean square error (RMSE) is 1.6378, the average absolute error (MAE) is 1.3585 and the prediction error is small.

Scenario Setting
In this paper, the landscape power generation system of a river basin in Xinjiang is taken as an example to verify and analyze the optimal dispatching model established above.The complementary system includes one hydropower station, a wind farm and a photovoltaic power station.The maximum output of the hydropower station is 720 MW and the minimum output is 60 MW.The hydropower station has a daily regulation capacity, a rated head of 68 m and a rated flow of 405 m 3 /s.The maximum capacity of the wind farm is 800 MW.The rated power of a single fan is 2 MW, and the maximum capacity of the photovoltaic power station is 1000 MW.Two typical days (a typical summer scene day and a typical winter scene day) were selected for analysis.The input and output data sizes are all 1 × 96, that is, the daily data with a time scale of 15 min.The wind power, photovoltaic power and load curves are shown in Figure 6.The correlation coefficient of the prediction results is 0.9954, the root mean square error (RMSE) is 1.6378, the average absolute error (MAE) is 1.3585 and the prediction error is small.As can be seen in Figure 6, in one day in this region, the load is smallest in the early morning, and the output of wind power is large in the night and early morning.The peak of the photovoltaic power generation in both summer and winter is at noon, and the photovoltaic output peaks at night and in the early morning.The output of the two kinds of renewable energy has certain complementary characteristics in terms of time.

Pareto Solution Set Analysis
The MOFPIO, the MOPIO, and the MOPSO are used to solve the multi-objective model of the water-wind-light complementary power generation system.The Pareto-optimal solution set that was obtained consists of three objective in two scenarios.In each scenario, the Pareto set of all three algorithms contains 50 solutions.The resulting Pareto-optimal solution set is shown in Figure 7.In the figure, f1 is the As can be seen in Figure 6, in one day in this region, the load is smallest in the early morning, and the output of wind power is large in the night and early morning.The peak of the photovoltaic power generation in both summer and winter is at noon, and the photovoltaic output peaks at night and in the early morning.The output of the two kinds of renewable energy has certain complementary characteristics in terms of time.

Experimental Verification 6.1. Pareto Solution Set Analysis
The MOFPIO, the MOPIO, and the MOPSO are used to solve the multi-objective model of the water-wind-light complementary power generation system.The Pareto-optimal solution set that was obtained consists of three objective functions in two scenarios.In each scenario, the Pareto set of all three algorithms contains 50 solutions.The resulting Pareto-optimal solution set is shown in Figure 7.In the figure, f 1 is the power generation and f 2 is the output fluctuation.As can be seen in Figure 7, the improved pigeon swarm algorithm is closer to the objective function value, and the optimal solution set in each scenario has the same changing trend.One of the two targets increases and the other decreases, which reflects the conflicting relationship between the objective functions.The difference between each scene is found to be in the form of the value range of the objective functions; the adaptive values of the target functions are counted in different scenarios.The results of the three strategies in the different scenarios are shown in Table 1.As can be seen in Figure 7, the improved pigeon swarm algorithm is closer to the objective function value, and the optimal solution set in each scenario has the same changing trend.One of the two targets increases and the other decreases, which reflects the conflicting relationship between the objective functions.The difference between each scene is found to be in the form of the value range of the objective functions; the adaptive values of the target functions are counted in different scenarios.The results of the three strategies in the different scenarios are shown in Table 1.The results are less than the other two algorithms under the same conditions, showing the better optimization accuracy of the MOFPIO algorithm.So, MOFPIO has excellent optimization performance.

Results of the Optimization Calculation
Figure 8 presents the results of the water-wind-light independent power generation and joint scheduling.As can be seen in the figure, after the complementary optimization scheduling of water-wind-light, the total power generation is 8.4722 million kWh.The hydraulic power generation is 4.5668 million kWh, accounting for 53.90% of the total power generation.The hydraulic power generation and photovoltaic power generation are 1.9242 million kWh and 1.9812 million kWh, accounting for 22.71% and 23.38% of the total power generation.It can be seen from the comparison of the separate output and the combined output that the output of the former hydropower fluctuates greatly with the time scale, while the output of the latter fluctuates less and greatly absorbs the solar energy.Figure 8 presents the results of the water-wind-light independent power generation and joint scheduling.As can be seen in the figure, after the complementary optimization scheduling of water-wind-light, the total power generation is 8.4722 million kWh.The hydraulic power generation is 4.5668 million kWh, accounting for 53.90% of the total power generation.The hydraulic power generation and photovoltaic power generation are 1.9242 million kWh and 1.9812 million kWh, accounting for 22.71% and 23.38% of the total power generation.It can be seen from the comparison of the separate output and the combined output that the output of the former hydropower fluctuates greatly with the time scale, while the output of the latter fluctuates less and greatly absorbs the solar energy.In order to study the complementary characteristics of the water-wind-light complementary optimization system, the individual output of each power supply is compared with the complementary output of the combined hydropower station.The degree of the intraday output fluctuation is also calculated.The complementarity index is given as the following equation: In order to study the complementary characteristics of the water-wind-light complementary optimization system, the individual output of each power supply is compared with the complementary output of the combined hydropower station.The degree of the Energies 2023, 16, 6787 15 of 18 intraday output fluctuation is also calculated.The complementarity index is given as the following equation: where P t is the average value of the total output of hydropower; P t is the total output of hydropower at time t; T is the number of daily scheduling cycles, and the value is 96.The calculation results of Equation ( 24) are shown in Table 2.It can be seen in Table 2 that the independent operation of hydropower, wind power and photovoltaic power generation fluctuates greatly, and the fluctuation decreases significantly after complementary operation.The output standard deviations of MOPSO, MOPIO and MOFPIO in winter after water-wind-light complementary operation are reduced by 31.8%,35.6% and 50.2%, respectively.The output standard deviations of MOPSO, MOPIO and MOFPIO in summer after operation are reduced by 14.7%, 34.1% and 84.5%, respectively, which are better than that of a single operation, indicating that the fluctuation of the system is significantly reduced after participating in the optimal scheduling.Among them, the water-wind-light output fluctuation suppression effect based on MOFPIO is the best, which can effectively suppress the fluctuation value of wind power.Figure 9 compares the total output value of hydropower on different typical days under the three methods.As can be seen in Figure 9, the MOPSO method has the most volatility, followed by the MOPIO method.The hydropower output fluctuation of the MOFPIO method is significantly lower than that of the other two methods.A small output fluctuation is conducive to the regulation of hydropower units.The typical hydropower output fluctuation in winter is larger than that of the typical days in summer, and the typical hydropower output in summer is relatively stable.The comparison results of the different scheduling strategies are shown in Figure 10.As can be seen in Figure 10, the total output fluctuation of the MOFPIO method is the smallest, followed by the MOPIO method and the MOPSO method.The three algorithms have the same effect over the two typical days because the hydropower output fluctuation of the typical winter day is larger than that of the typical summer day.Therefore, the total output fluctuation of the typical winter system is larger than that of the typical summer system, and MOFPIO method optimization scheduling can increase the stationarity of the load.

1 NFigure 3 .
Figure 3. Schematic diagram of the rotation angle, with the Wang−Mei index as th

Figure 3 .
Figure 3. Schematic diagram of the rotation angle, with the Wang−Mei index as the indicator.

1 .
Sorting Mechanism of ParetoThe ranking mechanism of Pareto consists of a non-dominant ranking operator and a crowding distance operator.First, individuals from a pigeon flock are divided into unequal non-dominance ranks based on inter-individual dominance relationships and fast non-dominance ordering.The set of non-dominated solutions for the current population consists of individuals that are located in the first rank.The crowding distance operator is reused to calculate the crowding distance between individuals in different non-dominance classes; individuals are then sorted again based on the size of the crowding distance.After the above two operators are completed, the individual pigeon is divided into different sets according to their non-dominant level, and the individuals in the set are arranged in descending order according to the size of the crowding distance.The flow of the ranking mechanism of Pareto is shown in Figure4.

Figure 5 .
Figure 5. Flowchart of the improved pigeon flock algorithm.

Figure 6 .
Figure 6.Wind power, photovoltaic power and load curve.

Figure 6 .
Figure 6.Wind power, photovoltaic power and load curve.

Figure 7 .
Figure 7. Pareto−optimal solution set under two objective functions.

Figure 7 .
Figure 7. Pareto−optimal solution set under two objective functions.

Figure 8 .
Figure 8.Comparison before and after joint water-wind-light scheduling.

Figure 8 .
Figure 8.Comparison before and after joint water-wind-light scheduling.

Table 1 .
Target function adaptation values in different scenarios.As can be seen in Table1, the search range for the f 1 of MOFPIO is about 12,000-25,000 MW.The search range of MOPIO is about 13,000-23,000 MW.And the search range of MOPSO is about 18,000-22,000 MW.It shows that the search range of MOFPIO is similar to that of MOPIO, which is significantly larger than that of MOPSO.The MOPIO algorithm has a better global search capability.For the f 2 , the minimum values obtained by the MOFPIO algorithm are 34,676.793MW and 17,148.577MW, respectively, in summer and winter, and the maximum values are 141,018.920MW and 176,687.155MW, respectively.

Table 1 .
Target function adaptation values in different scenarios.

Table 2 .
Comparison of the fluctuation value of hydropower output.