A multi-objective optimization model for the coordinated operation of hydropower and renewable energy

Abstract

The rapid growth of wind and solar energy sources in recent years has brought challenges to power systems. One challenge is surging wind and solar electric generation, understanding how to consume such generation is important. Achieving the complementarity of hydropower and renewable energies such as wind and solar power by utilizing the flexible regulation performance of hydropower is helpful to provide firm power to help renewable energy consumption. However, the multi-energy complementary operation mode will change the traditional hydropower operation mode, causing challenges to the comprehensive utilization of hydropower. In this paper, a multi-objective optimal scheduling model is built by considering coordinated hydro-wind-solar system peak shaving and downstream navigation. First, the Gaussian mixture model is adopted to quantify the uncertainty of wind and solar power. Then, a hydro-wind-solar coordinated model was built to obtain the standard deviation of the residual load and the standard deviation of the downstream water level. Finally, the ε-constraint method is used to solve for the Pareto optimality. The results demonstrate the following: 1) The proposed model can effectively determine hydropower output schemes that can coordinate wind and solar power output to reconcile peak shaving and navigation; 2) The downstream hydropower stations’ reverse regulation of the upstream hydropower station is a positive factor in reconciling conflicts; and 3) Reasonable planning of wind power and solar power is helpful for hydro-wind solar power complement operation.

1 Introduction

With the explosive growth of wind and solar power sources, grid-integrated variable renewable energy (VRE) has become a key part of achieving the “Carbon Peaking and Carbon Neutrality Goals” in China (Department of Resource Conservation and Environmental Protection, 2021). However, the inherent intermittent and random characteristics of wind and solar power seriously challenge the safety and reliability operation of power systems (Albadi and El-Saadany, 2010; Shivashankar et al., 2016; Jabir et al., 2017; Asiaban et al., 2021). Using flexible power sources to mitigate renewable intermittency is a key solution (Shivashankar et al., 2016). Hydropower is a flexible power source that has developed technology, a large scale, and a low cost. Thus, hydropower that complements VREs is important to support the wide-ranging integration of VREs such as wind and solar. However, hydropower that operates in unison with reservoirs usually serves multiple purposes. It not only needs to meet the demand for electricity generation, but also serves navigation, flood control, irrigation, ecological protection, and other comprehensive utilization tasks. Using hydropower to complement VREs will significantly change the operation mode of hydropower, and directly affect the power output, water level, and discharge flow of the reservoir. This may conflict with comprehensive utilization tasks (Jian et al., 2012; Shang et al., 2017). Therefore, exploring a reasonable operation method that considers the hydro-wind-solar complementation mode and comprehensive hydropower utilization is very significant for multi-energy complementation and the large-scale consumption of wind and solar.

Research on the operation of hybrid hydro-wind-solar systems, include studies (Zhu et al., 2018; Hu et al., 2021) that have evaluated the complementarity and united operation reliability of hydro-wind-solar power; moreover, various models suitable for different needs are established. In the literature (Liu et al., 2019a), a united operation strategy of hydro-wind-solar power, which can effectively suppress the intermittency of wind and solar output, is proposed by considering the natural complementarity of wind and solar power. The literature (Xie et al., 2021) shows that the introduction of spinning reserve and regulating reserve in the model can effectively deal with intermittent wind and solar power and improve the practical application ability of the hybrid hydro-wind-solar system. The literature (Jin et al., 2022) shows that properly developed wind and solar power can effectively utilize complementarity and reducing hydropower output fluctuations on multiple time scales. The literature (Wang et al., 2019b) simplifies the system by dividing the subsystems that can be coordinated to operate, and obtains the long-term scheduling scheme of the hydro-wind-solar power complementary system based on the whole basin. The literature (Chen et al., 2019) predicts wind and solar output through environmental conditions such as wind speed, radiation intensity, and temperature. On this basis, a hydro-wind-solar short-term complementary model that minimizes residual load fluctuations and maximizes VRE output is established.

However, with the rapid growth of integrated VREs, problems such as the risk of power output curtailment and the safe operation of high-proportion clean energy systems have become increasingly evident (Ding et al., 2016; Cheng et al., 2017; Luo et al., 2018), and scholars are constantly looking for solutions. The literature (Zhang et al., 2018) aimed to minimize the amount of curtailed wind and solar power and to maximize the stored energy in cascade hydropower stations and then established a short-term optimal scheduling model of hydro-wind-solar power. The literature (Liu et al., 2020) controlled the amount of curtailed wind and solar power by quantifying the uncertainty. The literature (Li et al., 2018) established an expansion planning method for a large-scale hybrid wind-solar multi-objective transmission network, which effectively reduced the amount of abandoned wind and solar energy in the power grid. The literature (Xi’an Jiaotong University et al., 2019) proposed a system for testing the reliability of large-scale renewable electricity integration and long-distance transmission. The literature (Ming et al., 2019) established a water-solar complementary long-term optimization model for large-scale solar power participation, and deduced a scheduling scenario that is superior to traditional rules in terms of energy production and power supply reliability.

Research on multi-energy complementarity with hydropower mostly regards hydropower stations as power production units. However, hydropower stations with reservoirs can also provide multiple benefits, such as ecological protection, irrigation, and navigation. Moreover, existing studies have shown that the participation of hydropower in peak regulation has adverse effects on the comprehensive benefits of reservoirs (Pérez-Díaz and Wilhelmi, 2010; Wan et al., 2020; Halleraker et al., 2022). Therefore, the reasonable operation of hydropower stations must also consider the conflict between different tasks. The literature (Li et al., 2020) used non-dominated sorting genetic algorithm to find the hydropower station scheduling scenario that considered the ecological needs of fish and hydropower generation. The literature (Niu et al., 2016) established a dual-objective model that takes both peak shaving and navigation under the background of demand conflict. The literature (Liu et al., 2019b) established multi-objective optimal dispatching model that considers the power generation of hydropower stations, power generation stability and ecological requirements of downstream rivers. The literature (Wang et al., 2019a) proposes a hydropower station dispatching plan that meets the requirements of hydropower compensation for wind power generation and daily water storage targets during the reservoir storage period.

The rest of this paper is structured as follows. In Section 2, the analysis method of wind-solar uncertainty and the composition of the multi-objective hydro-wind-solar collaborative scheduling model are described. In Section 3, the Gaussian mixture model, which is used to address uncertainty, and the ε-constraint method, which is used to solve the dual-objective model, are introduced. In Section 4, the model is applied to the Jinghong-Ganlanba cascade hydropower system that complements wind and solar power. In Section 5, the influence of wind and solar factors on the model results are analyzed. The main conclusions of the model analysis are given in Section 6.

2 Methodology

2.1 Uncertainty analysis of wind and solar power

It is challenging to forecast the power generation of wind and solar power due to the inherent intermittency, which is affected by environmental factors such as wind speed and light intensity. The literatures (Lingfors and Widén, 2016; Gholami et al., 2017) shows that the wind and solar output will tend to smooth out with the increase in the cluster scale of wind and solar stations. In this paper, since power stations in the same area usually send power via the same transmission channels, all wind stations in the same area are combined into a virtual wind power plant (VWP) and all solar power stations are combined into a virtual solar power plant (VSP).where is the number of wind power stations. is the number of solar power stations. is the period number, . is the number of time periods. is the wind power station number, . is the solar power station number, . Moreover, is the output of the virtual wind plant in period , MW; is the output of the virtual solar power plant in period , MW; is the output of wind power stations in period , MW; and is the output of solar power station in period , MW.

The forecast output error of the wind power plant and solar power plant is converted into the output error coefficient.where and are the forecast output error coefficients of the wind power and solar power at time , respectively; and are the actual output of wind power and solar power at time , respectively; and and are the forecast output of wind power and solar power at time , respectively.

Describing the uncertainty of wind and solar output reasonably is important for building a hybrid hydro-wind-solar complementary model. The Gaussian mixture model (GMM) is a kind of non-parametric estimation method. Compared with the parameter estimation method, it can theoretically describe any distribution. Compared with the kernel density estimation method, the GMM can avoid the bandwidth setting effect on the accuracy of the results. Therefore, in this paper, the GMM is adopted to describe the output uncertainty of wind and solar power. Given a random variable , the expression for the Gaussian distribution is:

The essence of the GMM is a simple linear combination of multiple Gaussian distributions, and its expression can be described as:where is the probability distribution of . is the weight coefficient, . Using the maximum likelihood estimation method, the parameters , , and can be determined.

2.2 Objective

With large-scale integration, the intermittent wind and solar power and the anti-peaking characteristics of wind power will challenge the peak shaving of the system. To meet the peak shaving requirements and balance the influence of wind and solar uncertainty, the flexible adjustment of hydropower output will cause frequent changes in downstream outflows and water level, which will directly affect navigation. Therefore, in this paper, peak shaving is considered the target of the hybrid hydro-wind-solar hybrid system:

represents the peak shaving target. Under certain conditions of incoming water, the goal of the hydro-wind-solar coordinate peak-shaving operation is to minimize the fluctuation of the residual load so that the residual load can be as smooth as possible. Therefore, the mathematical expression of the peak shaving target can be set to minimize the residual load variance:where is the residual load in the period, MW. is the original load of the grid. is the average residual load. is the sum of the output of the hybrid hydro-wind-solar system at Period . is the predicted output of the virtual wind plant at period ., MW. is the predicted output of the virtual solar power plant at period . is the output of the hydropower station at period . is the serial number of hydropower station, . is the water head of hydropower station at period . is the power generation flow of hydropower station at period . is the output calculation function of the hydropower station.

is the navigation target. The goal of navigation focuses on the elevation and variation in the downstream water level. Therefore, the minimum of the downstream water level change variance is taken as the objective function:where is the downstream water level of the downstream hydropower station at period , m. is the mean value of the downstream water level of the downstream hydropower station.

2.3 Constraints

2.3.1 System constraints

a) Constraints of hydro-wind-solar complementarity

Since wind and solar power are uncertain and difficult to forecast accurately, chance constraints are introduced to use hydropower flexibility to complement wind and solar power output uncertainty.

Eqs 14, 15 represent positive and negative output complement constraints, respectively. In these equations, represents the probability, and is the upper output limit of hydropower station at period . is the lower output of hydropower station at period . and are confidence levels.

b) Hybrid system reserve capacity requirements are satisfied with hydropower and can be expressed as a percentage of the power demand:where is the reserve capacity percentage factor.

2.3.2 Hydropower constraints

a) Mass balance constraints:where is the storage of hydropower station at period . , , , and are the natural incremental inflow, total discharge, turbine discharge and spillage of hydropower station at period respectively. is the water transportation time from reservoir to . is the time period duration.

b) Simultaneous regulation constraints:

For all hydropower stations in the same cascade hydropower system, the output in each period usually has the same trend of change.

c) Reservoir storage constraints:where and are the lower and upper storage bounds of hydropower station at period , respectively.

d) Total discharge constraints:where and are the lower and upper bounds of the total discharge of hydropower station at period .

e) Turbine discharge constraints:where and are the lower and upper bounds of the turbine discharge of hydropower station at period .

f) Initial storage and expected final storage and constraints:where is the initial storage volume of hydropower station . is the initial storage volume of hydropower station .

g) Navigation constraints:where is the downstream water level of the hydropower station at period , m. is the minimum downstream water level of the downstream hydropower station specified by the navigation requirements, m; and are fixed values selected according to the downstream waterway standard.

h) Water level-storage capacity function:where represents the water level-storage capacity calculation function. is the water level of hydropower station at period .

i) Downstream water level-discharge function:where represents the downstream water level-discharge function.

3 Solution method

3.1 Hydro-wind-solar complementary constraint processing

In this paper, the GMM is used to analyze the distribution of output error coefficients of wind power and solar power plants. The Gaussian mixture model (GMM) is a kind of non-parametric estimation method, it can theoretically describe any distribution. Literature (Pöthkow and Hege, 2013) compared the extract uncertainty contour line effects of the non-parametric models with Gaussian, and observed that Non-parametric models have good feasibility for various data sets. Therefore, it has widely used in uncertainty description for diverse applications (Potter et al., 2009; Mihai and Westermann, 2014), and compactly model relatively complex distributions (Liu et al., 2012).With the historical forecast and actual output coefficient data of the solar and wind power plant, the deviation value between the actual and the forecasted output coefficient can be calculated. The GMM is used to fit the deviation value data of the output coefficient, and the probability density function (PDF) of the output error coefficient is obtained as shown in Figure 1.

FIGURE 1

From the PDF graph of the output errors of the solar power plant and the wind power plant, the output coefficient deviation value can be obtained given the confidence interval [0.05, 0.95]. , , , and correspond to the endpoints of the confidence interval [0.05, 0.95] at each moment, and the upper/lower limit of the wind and solar output at each moment can be calculated:where , , , and are the lower limit of the wind power output range, the upper limit of the wind power output range, the lower limit of the solar output range and the upper limit of the solar output range during period , respectively.

Confidence in chance constraints (14) and (15) represents the probability that the power gird can consume all the wind and solar power, and whether the constraints are established is affected by the wind and solar forecast error. The chance constraint (14) means that when the hydropower is at its maximum value, the probability that the actual output of the hydro-wind-solar system is greater than the forecast output should be greater than the confidence. If the actual output of the wind and solar power is too small, the constraint may be destroyed. Therefore, the minimum output of wind power and the minimum output of solar power within the confidence level can be found through the CDF of wind and solar power forecast error. The opportunity constraint (14) can be transformed into that when the sum of wind power and solar power output is , the output of the hydro-wind-solar system meets the requirements:

The chance constraint (15) means that when the hydropower is at the minimum value, the probability that the actual output of the hydro-wind-solar system is less than the forecast output should be greater than the confidence. If the actual output of scenery is too large, the constraint may be destroyed. Similarly, the opportunity constraint (15) can be transformed into that when the sum of wind power and photovoltaic output is , the output of the hydro-wind-solar system meets the requirements:

3.2 method

The method simplifies the problem by retaining only one primary objective function and transforming the remaining objective functions into constraints (Mavrotas, 2009; De Santis et al., 2022); moreover, this method has been widely used in the solution of multi-objective models (Biswas et al., 2018; Fathipour and Saidi-Mehrabad, 2018). Compared with the weighted method, the method has a smaller search range and better performance in efficiently finding the Pareto front solution. Regarding the dual-objective optimization problem of finding the minimum value, its mathematical description is as follows:where is a matrix of appropriate dimensions; and is a column vector of appropriate dimensions.

For the dual objective, if Objective can be regarded as the primary objective function, then Objective is transformed into the constraints of . The specific operation steps are as follows:

Take as the objective function, find its minimum value under the original constraints, and find the value of at this time, which is regarded as the feasible maximum value of . The solution model required for this step is:

Take as the objective function, and find the minimum value of under the original constraints. The range of is . The solution model required for this step is:

Take points in , and generate new constraints: .

Take as the objective function and add the new constraints to the constraints. After the calculation is performed, result arrays can be obtained, that is, the approximate Pareto maximum excellent Frontier. The solution model required for this step is:

4 Case study

4.1 Background information

Yunnan Province is located in southwest China. In this region, with sufficient rainfall occurs, and several large rivers travels through, can be found. Yunnan province is endowed with abundant hydropower resources. In late 2022, the hydropower installed capacity in this region surged to 78.02 GW. This value is close to the hydropower installed capacity of Canada, which is the world’s third largest hydropower capacity. In addition, Yunnan Province is rich in renewable energy. In late 2022, wind power and solar power has surged to 8.83 and 4.17 GW respectively. However, there is very large potential for expanding wind and solar power resources, since the developed capacity accounts for only 5% of the developable capacity. According to the “New Energy Construction Plan of Yunnan Province in 2022” issued by Yunnan Province, Yunnan will accelerate the development of wind and solar power in the next few years. Since wind and solar power generation is depend on the weather, Yunnan is facing challenges of broad-scale renewable energy consumption and power grid peak shaving brought by the broad-scale integration of VREs. Therefore, using the flexibility of hydropower to achieve the multi-energy complementation of hydro-wind-solar power is one of the best ways for Yunnan to achieve large-scale consumption of VREs.

As shown in Figure 2, the Jinghong hydropower station and the Ganlanba hydropower station are both located in Jinghong City, Xishuangbanna Prefecture, Yunnan Province. They are the last two hydropower stations in the Lancang/Mekong River before flow out of China. Due to Jinghong’s rated flow and the restriction of downstream channel conditions, Jinghong’s adjustment of power generation flow will cause large fluctuations in downstream water level and flow, which will adversely affect navigation. To solve this problem, the main development task of the yet-to-be-constructed Ganlanba hydropower station is to provide reverse regulation for the Jinghong hydropower station to meet the downstream navigation requirements. Therefore, the Jinghong-Ganlanba cascade hydropwer station must undertake complex comprehensive utilization tasks. Considering Yunnan’s ever-increasing renewable energy construction plan, this cascade system will also undertake the severe task of VRES complementarity. Therefore, it is required to find a reasonable operation mode considering the conflicting objectives of generation and navigation. In this paper, the Jinghong-Ganlanba cascade hydropower stations are taken as the background. The load curve, interval runoff, initial and final water levels, and other conditions involved in the calculation process refer to the actual operating data of the power grid and power station, as shown in Supplementary Table S1.

FIGURE 2

In this paper, a non-linear model is built, where a quadratic polynomial fitting function is used to express the water level-storage function and tail water level-discharge function. The model was solved using Lingo18.0.

According to Section 2.1, all wind power plants in the same region are considered as one wind power plant with the installed capacity of 208.5 MW and a solar power plant with the installed capacity of 300 MW. Figure 3 illustrates the typical wind and solar power daily output during the dry season (January, February, and March) of the lower reaches of the Lancang River, which were chosen to evaluate the scheduling situation of cascade hydropower stations under difficult conditions.

FIGURE 3

4.2 Analysis of the results

Since Ganlanba has not been constructed, three cases are considered to analyze the impact of the construction of the Ganlanba on the overall results and the impact of the navigation requirements on the separate operation of the Jinghong hydropower station: 1) Case 1: only the Jinghong hydropower plant is considered in the hybrid hydro-wind-solar complement operation model, but does not take navigation constraints into account; 2) Case 2: only the Jinghong hydropower plant is considered and takes navigation constraints into account and 3) Case 3: both Jinghong and Ganlanba are considered, and navigational constraints are also taken into account.

Figure 4 and Supplementary Table S2 show the distribution and range of the Pareto front for the three cases which reflect the conflict between the two objective functions. The rectangular area in Figure 4 marks solutions at the turning part of the Pareto front, which are hydropower output processes that can better meet the requirements of peak shaving and navigation at the same time. Under the experimental conditions, the Pareto front turning of the cascade solution is more obvious, indicating that it can moderate the conflict more effectively. Because the navigation constraints limit the output of Jinghong, the range of feasible solutions is reduced, and even some solutions with good peak-shaving effects turn into infeasible solutions.

FIGURE 4

Figure 5 illustrate the differences of solutions on the Pareto front for Case1 and show the changes in the output and water abandonment under Case1, respectively, by selecting different downstream water level variance limits ( = 1.4, 1.0, 0.6, 0.2, 0.01, 0.0001). When the target requirements for tighten, Jinghong’s discharge adjustment ability is restricted, which result in Jinghong abandoning water for peak shaving. For Case1, the peak shaving and navigation task requirements cannot be achieved at the same time.

FIGURE 5

In the Pareto front of Case 1, three scenarios of , , and are selected as the solution that prefers a better peak shaving effect, the neutral case, and the solution that prefers a better navigation effect, respectively. These three points can effectively compare the gaps in peak regulation and navigation between these cases, as shown in Supplementary Table S3.

Supplementary Table S3 summarizes the eigenvalues of each case. It is obvious that under the same , Case 3’s peak shaving impact is superior based on the peak-to-valley difference of the residual load. Since the target tends to reduce the hydropower output at the load valley, the hydroelectric plant can only supply reserve capacity at the load valley leaving a constant minimum load at 6178.285 MW. Figure 6 show the output and water abandonment of Case 3 under three conditions. From Figure 5, 6, it can be concluded that: 1) in Figure 6, the output of Ganlanba slightly changes under different conditions, which indicates that Ganlanba has a weak adjustment ability and is mainly responsible for navigation tasks. Therefore, Figure 5, 6 show that with the addition of Ganlanba, Jinghong’s adjustment capacity can be unleashed, and a better result can be obtained. 2) Figure 5, 6 also show that Ganlanba’s participation can significantly reduce the abandoned water for peak shaving, which may successfully address conflicting objectives of peak shaving and navigation.

FIGURE 6

Figure 7 shows the downstream water levels before and after the addition of Ganlanba and displays the affected navigation due to Ganlanba’s involvement. When the downstream water level variance is set to 0.01, the corresponding solution of Case 3 is located at the turning point of the Pareto front. In Figure 7, the Pareto front of Case 2 also has a corresponding solution, which can be compared. Without Ganlanba, the navigable water level is maintained at 537 m for 1/3 of the time, which indicates that the downstream channel barely meets the navigation requirements. For Case 3, the Ganlanba’s addition lowers the expense of maintaining navigation, and it maintains the navigation target’s ideal value between 529.45 m and 529.78 m. The safety of downstream navigation has been further assured compared to the minimum water level requirement of 525 m downstream of the Ganlanba.

FIGURE 7

In conclusion, the Ganlanba hydropower plant can successfully reduce the tension between peak shaving and navigation, address the issue of water abandonment for peak shaving, and raise the level of safety for downstream navigation.

5 Impact analysis of the wind and solar output

Due to the randomness and uncontrollability of the wind and solar output, in actual operation, the real output of wind and solar power usually deviates from the predicted value. In this section, to analyze the influence of wind and solar power uncertainty on scheduling, 1) the residual load under various wind and solar output conditions is compared to analyze the influence of with/without considering hydro-wind-solar complementary constraints on the calculation results of the dual-objective model. 2) The hydropower output under various wind and solar capacity scenarios is evaluated to examine the impact of wind and solar power scale on the calculations of the dual-objective model. Therefore, 1,000 wind and solar output scenarios were generated by using the Monte Carlo method. Among the seven cases selected in Section 4.2, Case 3 () works best. Thus, in the cascade system is taken as an example.

5.1 Influence of the wind and solar uncertainty on the peak shaving effect

Literature (Shen et al., 2021) proposes that the load fluctuation coefficient can be used to evaluate the peak shaving effect, and the revised index mainly reflects the overall smoothness of the load curve, which is better than the visualization effect of the residual load mean square error. Therefore, in this paper, this index is used to reflect the peak shaving effect in different scenarios:where, is the residual load of the scenario at the period ; and is the average residual load of the scenario , with .

In Case 3, the actual operation of the cascade is simulated under two conditions: one condition with complementary constraints for the solar and wind output (Eqs 33, 34) and the other condition without such constraints. When there is no hydro-wind-solar complement constraint, the residual load in each scenario is:where, is the residual load of the th scenario at period without hydro-wind-solar complement constraints; is the output of hydropower station at period ; is the wind and solar power output of the scenario at period ; , represent the wind and solar power output of scenario at period respectively.

When considering the hydro-wind-solar complementary constraint, here are the residual loads for each scenario:where is the upper limit of the range for the sum output of wind and solar; and is the lower limit of the range for the sum output of wind and solar.

Through Eqs 40–45, the residual load under 1,000 simulation scenarios was carried out, and the distribution of the load fluctuation coefficients is shown in Figure 8. It is obvious that considering hydro-wind-solar complementarity can effectively control the residual load fluctuation in real operation where actual wind and solar power deviate from the forecast value. In summary, hydro-wind-solar complement operation can effectively reduce the adverse impact of wind and solar uncertainty on the peak shaving results.

FIGURE 8

5.2 Influence of the wind and solar uncertainty on the navigation

Under the condition of fixed hydropower output and without providing wind and solar reserve capacity, the navigation effect will not be affected by VRE when actual wind and solar output deviate from the forecast value. In the case of hydro-wind-solar complement operation, the part of the output deviation between the upper and lower limits of the reserve capacity needs to be absorbed by changing the hydropower output. Due to the small adjustment capacity of Ganlanba, in this situation, Jinghong is required to undertake the adjustment tasks. Since the navigation is determined by the discharge of Ganlanba, the picture A of Figure 9 shows the maximum and minimum water levels of the Ganlanba reservoir and the upper limit and lower limit water levels of the Ganlanba reservoir under 1,000 scenarios. In all scenarios, Ganlanba reservoir’s water level will not exceed its limit, which means that when wind and solar forecast errors cause the upper reservoir outflow to change, Ganlanba has sufficient adjustment capacity to address it. The picture B of Figure 9 shows the downstream water level variance of Ganlanba Reservoir under 1,000 scenarios. The mean square deviations of the downstream water level in most scenarios are smaller than the predicted situation, and a few are slightly larger than the predicted situation. However, the difference is small and negligible. This demonstrates that Ganlanba can maintain navigational conditions despite VREs’ forecast errors.

FIGURE 9

That is, wind and solar uncertainties have little effect on navigation.

5.3 The influence of the expanding wind power and solar power scale on hydro-wind-solar complementarity

In this section, different wind and solar power installed capacity scales are selected to explore their impact on the simulation results. There shows four wind and solar power scales. Scenario 1 elects solar power scales of 300 MW and wind power scales of 208.5 MW, which are the wind and solar installed capacity in 2020. Referring to the “List of New Energy Projects in the 14th Five-Year Plan of Yunnan Province,” Xishuangbanna will add 1.43 GW of solar power during the “14th Five-Year Plan” period. However, no wind power projects are planned. Scenario 2 sets the estimated installed capacity of solar and wind power in 2025 as 1730 and 208.5 MW. Scenario 3 sets the projected installed solar and wind capacity in 2030 as 3160 and 208.5 MW, and the growth rate is consistent with that in Scenario 2. Scenario 4 doubles the wind-power capacity installed based on Scenario 3, and sets the estimated installed capacity of solar and wind power as 3160 and 417 MW.

Figure 10 shows the residual load curves of each scenario, and Figure 11 shows the hydropower output gap between Scenario 1 and Scenario 2, 3 and 4, in which a negative value indicates that hydropower output will decrease as solar/wind power output increases, while a positive value indicates the opposite. As shown in these two figures, for Scenario 2 and 3, when solar power output increases, the residual load reduce from 8:30 to 20:00, and hydropower will reduce output during the load saddle periods and uses the saved water to increase the output during the load peak periods. This is because solar power is mainly generated during the daytime which is the load saddle, and hydropower adapts to achieve hydro-solar complementarity and improve the peak shaving effect. The simulation results of Scenarios 2 and 3 in Supplementary Table S4 also prove that increasing the output of solar power has a positive effect on the peak shaving of the model.

FIGURE 10

FIGURE 11

As shown in Table 4, the peak-to-valley difference of the residual load increases from Scenario 3 to Scenario 4. This is because hydropower has reached the best peak shaving performance in Scenario 3. When the wind power capacity is increased from Scenario 3 to Scenario 4, the hydropower output changes little with the change of reserve capacity, as shown in Figure 11, and helps little in further peak shaving. Therefore, with the increase in wind power capacity, although both the residual peak load and valley load decrease, as the anti-peak characteristic of wind power shows in Figure 3., the peak-valley difference will increase.

Figure 11 shows that the hydropower output changes little, which means that the outflow of hydropower changes little. From Scenario 1 to Scenario 4, the increase in wind power and solar power scale will not have a large negative impact on navigation.

In summary, during the dry season, the increase in the scale of solar power output has a positive impact on the peak shaving of the hydro-wind-solar complementarity system; in contrast, wind power tends to have a negative impact on peak shaving.

6 Conclusion

References

• 1

  AlbadiM. H.El-SaadanyE. F. (2010). Overview of wind power intermittency impacts on power systems. Electr. Power Syst. Res.80, 627–632. 10.1016/j.epsr.2009.10.035

  • CrossRef
  • Google Scholar

• 2

  AsiabanS.KayedpourN.SamaniA. E.BozalakovD.De KooningJ. D. M.CrevecoeurG.et al (2021). Wind and solar intermittency and the associated integration challenges: A comprehensive review including the status in the Belgian power system. Energies14, 2630. 10.3390/en14092630

  • CrossRef
  • Google Scholar

• 3

  BiswasP. P.SuganthanP. N.QuB. Y.AmaratungaG. A. J. (2018). Multiobjective economic-environmental power dispatch with stochastic wind-solar-small hydro power. Energy150, 1039–1057. 10.1016/j.energy.2018.03.002

  • CrossRef
  • Google Scholar

• 4

  ChenS.FangG.HuangX.YanM. (2019). A joint optimal dispatching method of wind-solar-hydro generation system. IOP Conf. Ser. Earth Environ. Sci.227, 032004. 10.1088/1755-1315/227/3/032004

  • CrossRef
  • Google Scholar

• 5

  ChengC.YanL.MirchiA.MadaniK. (2017). China’s booming hydropower: Systems modeling challenges and opportunities. J. Water Resour. Plann. Manage.143, 02516002. 10.1061/(ASCE)WR.1943-5452.0000723

  • CrossRef
  • Google Scholar

• 6

  De SantisM.EichfelderG.PatriaD. (2022). On the exactness of the ε-constraint method for biobjective nonlinear integer programming. Operations Res. Lett.50, 356–361. 10.1016/j.orl.2022.04.007

  • CrossRef
  • Google Scholar

• 7

  Department of Resource Conservation and Environmental Protection (2021). Action plan for carbon dioxide peaking before 2030. Available at: https://en.ndrc.gov.cn/policies/202110/t20211027_1301020.html.

  • Google Scholar

• 8

  DingM.XuZ.WangW.WangX.SongY.ChenD. (2016). A review on China׳s large-scale PV integration: Progress, challenges and recommendations. Renew. Sustain. Energy Rev.53, 639–652. 10.1016/j.rser.2015.09.009

  • CrossRef
  • Google Scholar

• 9

  FathipourF.Saidi-MehrabadM. (2018). A multi-objective energy planning considering sustainable development by a TOPSIS-based augmented e-constraint. J. Renew. Sustain. Energy10, 034901. 10.1063/1.5008545

  • CrossRef
  • Google Scholar

• 10

  GholamiM.FathiS. H.MilimonfaredJ.ChenZ.DengF. “The effect of turbulence and wake on the power fluctuation in the wind farms,” in 2017 Iranian Conference on Electrical Engineering (ICEE), Tehran, Iran, March, 2017 (IEEE), 1148–1153. 10.1109/IranianCEE.2017.7985214

  • CrossRef
  • Google Scholar

• 11

  HallerakerJ. H.KenawiM. S.L’Abée-LundJ. H.BakkenT. H.AlfredsenK. (2022). Assessment of flow ramping in water bodies impacted by hydropower operation in Norway – is hydropower with environmental restrictions more sustainable?Sci. Total Environ.832, 154776. 10.1016/j.scitotenv.2022.154776

  • CrossRef
  • Google Scholar

• 12

  HuW.WangY.SunY.NieQ.DingR.ZhangX. “Research on comprehensive complementary characteristics evaluation technology of wind-solar-hydro combined power generation system,” in 2021 International Conference on Power System Technology (POWERCON), Haikou, China, September, 2021 (IEEE), 1017–1022. 10.1109/POWERCON53785.2021.9697711

  • CrossRef
  • Google Scholar

• 13

  JabirM.Azil IlliasH.RazaS.MokhlisH. (2017). Intermittent smoothing approaches for wind power output: A review. Energies10, 1572. 10.3390/en10101572

  • CrossRef
  • Google Scholar

• 14

  JianZ.ShiqiangZ.XizhangC. (2012). Impact upon navigation conditions of river reach between the two dams by peak shaving at three gorges hydropower station. Procedia Eng.28, 152–160. 10.1016/j.proeng.2012.01.698

  • CrossRef
  • Google Scholar

• 15

  JinX.LiuB.LiaoS.ChengC.LiG.LiuL. (2022). Impacts of different wind and solar power penetrations on cascade hydroplants operation. Renew. Energy182, 227–244. 10.1016/j.renene.2021.10.022

  • CrossRef
  • Google Scholar

• 16

  LiF.-F.WangH.-R.WuZ.-G.QiuJ. (2020). Maximizing both the firm power and power generation of hydropower station considering the ecological requirement in fish spawning season. Energy Strategy Rev.30, 100496. 10.1016/j.esr.2020.100496

  • CrossRef
  • Google Scholar

• 17

  LiH.GaoH.ZhangY.AiX. (2018). Expansion planning of large scale hybrid wind\-photovoltaic transmission network considering correlation. Power Syst. Technol.42, 2120–2127. 10.13335/j.1000-3673.pst.2018.0144

  • CrossRef
  • Google Scholar

• 18

  LingforsD.WidénJ. (2016). Development and validation of a wide-area model of hourly aggregate solar power generation. Energy102, 559–566. 10.1016/j.energy.2016.02.085

  • CrossRef
  • Google Scholar

• 19

  LiuB.LundJ. R.LiaoS.JinX.LiuL.ChengC. (2020). Optimal power peak shaving using hydropower to complement wind and solar power uncertainty. Energy Convers. Manag.209, 112628. 10.1016/j.enconman.2020.112628

  • CrossRef
  • Google Scholar

• 20

  LiuJ.ShenW.YangH.SuK.ZhouH. “Operation strategy of a combined system consists of multi-type hydropower and wind-solar,” in 2019 IEEE 3rd Conference on Energy Internet and Energy System Integration (EI2), Changsha, China, November, 2019a (IEEE), 369–374. 10.1109/EI247390.2019.9062166

  • CrossRef
  • Google Scholar

• 21

  LiuS.LevineJ. A.BremerP.PascucciV. (2012). “Gaussian mixture model based volume visualization,” in IEEE symposium on large data analysis and visualization (LDAV) (Seattle, WA: IEEE), 196, 73–77. 10.1109/LDAV.2012.6378978

  • CrossRef
  • Google Scholar

• 22

  LiuW.ZhuF.ChenJ.WangH.XuB.SongP.et al (2019b). Multi-objective optimization scheduling of wind–photovoltaic–hydropower systems considering riverine ecosystem. Energy Convers. Manag.196, 32–43. 10.1016/j.enconman.2019.05.104

  • CrossRef
  • Google Scholar

• 23

  LuoG.DanE.ZhangX.GuoY. (2018). Why the wind curtailment of northwest China remains high. Sustainability10, 570. 10.3390/su10020570

  • CrossRef
  • Google Scholar

• 24

  MavrotasG. (2009). Effective implementation of the ε-constraint method in multi-objective mathematical programming problems. Appl. Math. Comput.213, 455–465. 10.1016/j.amc.2009.03.037

  • CrossRef
  • Google Scholar

• 25

  MihaiM.WestermannR. (2014). Visualizing the stability of critical points in uncertain scalar fields. Computers & Graphics41, 13–25. 10.1016/j.cag.2014.01.007

  • CrossRef
  • Google Scholar

• 26

  MingB.LiuP.GuoS.ChengL.ZhangJ. (2019). Hydropower reservoir reoperation to adapt to large-scale photovoltaic power generation. Energy179, 268–279. 10.1016/j.energy.2019.04.209

  • CrossRef
  • Google Scholar

• 27

  NiuW.ShenJ.ChengC.WuX.FengZ.GuoY. (2016). A hybrid search method for multi-objective optimization operations of cascaded hydropower plants with peak load regulation and navigation demands. Proc. Chin. Soc. Electr. Eng.36, 2331–2341. 10.13334/j.0258-8013.pcsee.2016.09.003

  • CrossRef
  • Google Scholar

• 28

  Pérez-DíazJ. I.WilhelmiJ. R. (2010). Assessment of the economic impact of environmental constraints on short-term hydropower plant operation. Energy Policy38, 7960–7970. 10.1016/j.enpol.2010.09.020

  • CrossRef
  • Google Scholar

• 29

  PöthkowK.HegeH.-C. (2023). Nonparametric models for uncertainty visualization. Computer Graphics Forum32, 131–140. 10.1111/cgf.12100

  • CrossRef
  • Google Scholar

• 30

  PotterK.WilsonA.BremerP.-T.WilliamsD.DoutriauxC.PasV.et al (2009). Ensemble-vis: a framework for the statistical visualization of ensemble data," in. 2009 IEEE international conference on data mining workshops (Miami, FL: IEEE), 233–240. 10.1109/ICDMW.2009.55

  • CrossRef
  • Google Scholar

• 31

  ShangY.LiX.GaoX.GuoY.YeY.ShangL. (2017). Influence of daily regulation of a reservoir on downstream navigation. J. Hydrol. Eng.22, 05017010. 10.1061/(ASCE)HE.1943-5584.0001522

  • CrossRef
  • Google Scholar

• 32

  ShenJ.ChenG.WeiW.ChengC.WangL. (2021). Short-term optimal scheduling model with multiple peak shaving indexes based on linkage analysis. Journal of Hydraulic Engineering52, 936–947.

  • Google Scholar

• 33

  ShivashankarS.MekhilefS.MokhlisH.KarimiM. (2016). Mitigating methods of power fluctuation of photovoltaic (PV) sources – a review. Renew. Sustain. Energy Rev.59, 1170–1184. 10.1016/j.rser.2016.01.059

  • CrossRef
  • Google Scholar

• 34

  WanZ.LiY.WangX.AnJ.DongB.LiaoY. (2020). Influence of unsteady flow induced by a large-scale hydropower station on the water level fluctuation of multi-approach channels: A case study of the three gorges project, China. Water12, 2922. 10.3390/w12102922

  • CrossRef
  • Google Scholar

• 35

  WangX.VirguezE.KernJ.ChenL.MeiY.Patiño-EcheverriD.et al (2019a). Integrating wind, photovoltaic, and large hydropower during the reservoir refilling period. Energy Convers. Manag.198, 111778. 10.1016/j.enconman.2019.111778

  • CrossRef
  • Google Scholar

• 36

  WangX.VirguezE.XiaoW.MeiY.Patiño-EcheverriD.WangH. (2019b). Clustering and dispatching hydro, wind, and photovoltaic power resources with multiobjective optimization of power generation fluctuations: A case study in southwestern China. Energy189, 116250. 10.1016/j.energy.2019.116250

  • CrossRef
  • Google Scholar

• 37

  Xi’an Jiaotong UniversityWangJ.WeiJ.Xi’an Jiaotong UniversityZhuY.Xi’an Jiaotong Universityet al (2019). The reliability and operation test system of power grid with large-scale renewable integration. CSEE JPES. 10.17775/CSEEJPES.2017.01300

  • CrossRef
  • Google Scholar

• 38

  XieJ.ZhengY.PanX.ZhengY.ZhangL.ZhanY. (2021). A short-term optimal scheduling model for wind-solar-hydro hybrid generation system with cascade hydropower considering regulation reserve and spinning reserve requirements. IEEE Access9, 10765–10777. 10.1109/ACCESS.2021.3049280

  • CrossRef
  • Google Scholar

• 39

  ZhangX.MaG.HuangW.ChenS.ZhangS. (2018). Short-term optimal operation of a wind-PV-hydro complementary installation: Yalong river, sichuan province, China. Energies11, 868. 10.3390/en11040868

  • CrossRef
  • Google Scholar

• 40

  ZhuY.ChenS.HuangW.WangL.MaG. (2018). Complementary operational research for a hydro-wind-solar hybrid power system on the upper Jinsha River. J. Renew. Sustain. Energy10, 043309. 10.1063/1.5027758

  • CrossRef
  • Google Scholar