The decision-making process of wind–photovoltaic–hydropower systems involves knowledge from many fields. Influenced by the knowledge level of the decision-maker and the attribute information of the scheme set, there exists a certain uncertainty in the indicator weights. In view of this, this paper proposes a stochastic multi-criteria decision-making framework for scheduling of wind–photovoltaic–hydropower systems, which overcomes the difficulty of uncertainty in indicator weights or even completely unknown information about indicator weights at the time of decision-making. The Stochastic Multi-criteria Acceptability Analysis (SMAA) theory and the VIKOR model are introduced, and the proposed SMAA–VIKOR model makes the indicator weight space explicit. The study shows that the proposed SMAA–VIKOR model can overcome the obstacle of decision-makers’ lack of information on indicator weights. The ranking acceptability indicators calculated by the model show a more obvious trend of advantages and disadvantages, which gives full confidence to the decision-making group to formulate a plan to be implemented. It breaks through the bottleneck of group decision-making, which is difficult to make effective decisions due to the condition of incomplete information, and enriches the library of stochastic multi-criteria decision-making methods for the scientific formulation of scheduling schemes of wind–photovoltaic–hydropower systems under uncertainty conditions.

  • A stochastic multi-criteria decision-making framework for wind–photovoltaic–hydropower systems is proposed.

  • The SMAA–VIKOR model is proposed to clarify the indicator weight space.

  • The proposed SMAA–VIKOR model can overcome the obstacle of decision-makers’ lack of information on indicator weights.

In response to global climate change, more and more countries are developing new energy sources such as wind power and photovoltaic power generation (Rigatos et al. 2019; Russo et al. 2023; Xia et al. 2023). The scientific development of wind–photovoltaic–hydropower system scheduling scheme is an important basis for guaranteeing system efficiency, system safety, downstream river health, reservoir safety, and other goals (Zhang et al. 2021). However, the wind–photovoltaic–hydropower system scheduling scheme is affected by many uncertainties in the scheduling process, such as the subjective preference of decision-makers and the importance of the indicators themselves (Liu et al. 2020).

A number of multi-criteria decision-making methods have been applied to scheduling of wind–photovoltaic–hydropower systems. In the deterministic field, Kang et al. (2011) used the Fuzzy Analytical Hierarchy Process (FAHP) to analyze the indicators of benefits, opportunities, costs and risks of wind farms to assess the expected comprehensive benefits of wind farm projects in order to select the most appropriate wind farm construction option in the comprehensive wind farm siting decision. Based on the subjective preference of the decision-maker and the improved entropy weight to determine the target weights, combined with the fuzzy set theory, Lu et al. (2011) proposed a decision-making method for multi-objective joint scheduling of a group of reservoirs. Perera et al. (2013) obtained Pareto frontiers for the objectives of energy cost, load deficit, energy wastage and fuel consumption through multi-objective optimization during the design of a hybrid energy systems (HESs), and then considered the existence of ambiguity in the level of importance between the indicators, the weights of the indicators were processed using the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS). Sanchez-Lozano et al. (2016) analyzed 10 indicators related to wind farm siting by coupling FAHP and fuzzy TOPSIS (FTOPSIS) methods in wind farm siting assessment. The study shows that the method is suitable not only for evaluating quantitative indicators, but also for qualitative indicators.

In the field of uncertainty, Liu et al. (2019) considered different decision-makers' knowledge and preference of indicators in the group decision-making process in the formulation of scheduling scheme for wind and hydropower system, and considered that the weights of indicators obeyed the uniform distribution. The second generation of SMAA (SMAA-2) was used for decision-making, and a scheme that meets the preferences of decision-makers was obtained. Qin et al. (2010) described the indicators such as maximum water level, maximum discharge flow, annual power generation and minimum output as random variables obeying interval normal distribution, and proposed a risk-based multi-criteria decision-making method based on the combination of dominant likelihood degree and comprehensive assignment. Zhu et al. (2017) proposed a stochastic multi-criteria decision-making model for reservoir flood control scheduling based on Qin's work by assuming the indicator weights to be uniformly and normally distributed, in addition to considering that the indicator values obey the normal distribution, and that there may be conflicts in the indicator weights. The results show that the model has sorted acceptability indicators, center vectors and global acceptability indicators, which can provide rich decision-making information for decision-makers when making decisions.

The development of scheduling schemes for wind–photovoltaic–hydropower systems involves knowledge from a number of fields, including wind power, photovoltaic power, hydropower, ecology and group decision-making, and requires a high level of knowledge on the part of the decision-maker in order to develop a scientifically sound scheduling scheme (Tan et al. 2021). At the same time, it is difficult to obtain a scheduling scheme that satisfies the preferences of the decision-making group due to conflicting preferences and competing objectives in the group decision-making process (Wang et al. 2022). The subjective preferences of different interest decision-makers are not consistent, and the weights of indicators derived from different objective weighting methods are not the same, leading to a certain degree of uncertainty in the weights of indicators (Chen et al. 2023). Therefore, this paper proposes a stochastic multi-criteria decision-making framework for wind–photovoltaic–hydropower systems. In response to the decision-making process, decision-makers are afraid to express their own preference information due to their own knowledge limitations, SMAA-2 and VIKOR (VlseKriterijumska Optimizacija I Kompromisno Resenje) models are reviewed, and the coupled SMAA–VIKOR model is used to analyze the inverse weight space of each completely unknown indicator weight, and to clarify the indicator weight space. In the later stages of decision-making, as information availability continues to improve, the decision-maker's understanding of the options and indicator weights becomes clearer, but there is still a corresponding ambiguity in the indicator weights, which relies on the partly clear and partly ambiguous subjective opinions of the decision-making group. For this reason, the Intuitionistic Fuzzy Analytic Hierarchy Process (IFAHP) is introduced to allow decision-making groups to express their fuzzy preference information. Finally, the IFAHP- SMAA–VIKOR model is developed to scientifically develop long-term scheduling schemes for wind–photovoltaic–hydropower systems. The proposed stochastic multi-criteria decision-making framework for wind–photovoltaic–hydropower systems overcomes the difficulty of uncertainty in indicator weights or even completely unknown information of indicator weights when making decisions. This study enriches and refines the stochastic multi-criteria decision-making methods under uncertainties.

The rest of this paper is organized as follows. Section 2 provides the details of the SMAA-2 Model, VIKOR Model, and group decision-making model for wind–photovoltaic–hydropower systems. Section 3 provides the details of the case study, which includes the analysis of sampling method, utility functions, and deterministic and stochastic multi-criteria decision-making. Section 4 presents the discussion and conclusions of this study.

The development of scheduling scheme for wind–photovoltaic–hydropower systems involves knowledge from multiple disciplines such as wind power, photovoltaic power, hydropower, ecology, and group decision-making, which requires the decision-making group to be familiar with these disciplines, and it is difficult to develop a scientifically sound scheduling scheme. Decision-attribute information accompanied by the subjective will of the decision-maker and its own attributes exists a certain uncertainty, even the indicator weight information is completely unknown. In view of this, we reviewed the SMAA-2 model to carry out inverse weight space analysis of indicator weight space to overcome the obstacle of no weight information, reviewed the VIKOR model to improve the linear summation type utility function in SMAA-2, and established a new SMAA–VIKOR model to derive a scientific and reasonable scheduling scheme for wind–photovoltaic–hydropower systems.

SMAA-2 model

Lahdelma & Salminen (2001) proposed a SMAA-2 that can take into account the weight information and the uncertainty of the alternatives. SMAA-2 is based on a decision model that performs inverse weight space analysis on indicator weights and calculates the probability of each alternative becoming the optimized one, which overcomes the problem that traditional deterministic decision models are not applicable to multi-criteria decision problems with unknown weight information.

It is assumed that the problem to be decided consists of M schemes and N attributes , where denotes the Mth solution and denotes the Nth attribute. The attribute weights can be expressed as and denotes the weight of attribute . The indicator evaluation matrix is , as shown in Equation (1):
(1)
Uncertainty in stochastic multi-criteria decision-making arises from both indicator evaluation values and indicator weights. In general, uncertainty in indicator evaluation values arises from uncertainty in the outcome of stochastic optimization. Indicator weights reflect information on both the level of importance of the indicator itself and the subjective preferences of the decision-maker, which is the source of uncertainty in indicator weights. The uncertainty in both the indicator evaluation values and the indicator weights can be described using the corresponding probability density functions. Therefore, the indicator evaluation value of the above equation can be described by a random variable obeying a probability density function , as shown in Equation (2):
(2)
where is the indicator evaluation value for the nth attribute of the mth scheme.

In the stochastic multi-attribute decision-making process, the indicator weight space W can be represented by a probability density function , and there are usually four possibilities, which are (1) the indicator weight values are uniquely determined, (2) the indicator weight values obey a uniform distribution in the specified interval, (3) the indicator weight values obey an arbitrary type of distribution in the specified interval, and (4) the information about the indicator weights is completely unknown (Barron & Barrett 1996).

  • (1)

    The indicator weight values are uniquely determined.

The feasible weight space represents the preferences of the decision-maker. When the indicator weights are deterministic values, the feasible weight space will be uniquely determined. In this paper, taking the three-dimensional space as an example, the feasible weight space is a point in the three-dimensional space, where , and represent the indicator weights, respectively, as shown in Figure 1.
  • (2)

    The indicator weight values obey a uniform distribution in the specified interval.

Figure 1

Criteria weight with deterministic information.

Figure 1

Criteria weight with deterministic information.

Close modal
The complexity of indicators and the uncertainty of people's subjective preferences may make it difficult for decision-makers to reach an agreed understanding of weights in the face of incomplete information. At this time, if we use purely determined numerical information to express the corresponding weights, it may not be able to cover the information of the indicator itself and the subjective preference of the decision-makers. The interval estimation based on weights can express such information well, avoiding the corresponding lack of information. In this paper, taking the three-dimensional spatial decision-making as an example, it is specified that the indicator weights obey the uniform distribution on the weight intervals, as shown in Figure 2, and the specific formula is shown in Equation (3):
(3)
  • (3)

    The indicator weight values obey an arbitrary type of distribution in the specified interval.

Figure 2

Criteria weight follows uniform weight distribution with interval constraint.

Figure 2

Criteria weight follows uniform weight distribution with interval constraint.

Close modal
In addition to the above indicator weight values obeying uniform distribution in the specified interval, the indicator weight values can also be described in a more general form, i.e., the indicator weight values obeying any type of distribution in the specified interval. In this paper, taking the three-dimensional decision space as an example, the indicator weight space is described as the shaded part shown in Figure 3, and the specific formula is shown in Equation (4):
(4)
  • (4)

    The information about the indicator weights is completely unknown.

Figure 3

Criteria weight follows any weight distribution with interval constraint.

Figure 3

Criteria weight follows any weight distribution with interval constraint.

Close modal
At the early stage of multi-attribute decision-making, the decision-maker may be affected by the intertwining of their own knowledge and the decision-making environment, and it is difficult to accurately give the corresponding weighting information of the indicators. The decision-making model proposed in this paper can be adapted to search the feasible domain of the whole decision space through the idea of inverse weight analysis, so as to provide the decision-maker with a corresponding reference when the decision-maker cannot give any preference information at the early stage of decision-making. Taking the three-dimensional decision space as an example, the indicator weight space is described as a triangular region as shown in Figure 4, and the specific formula is shown in Equation (5):
(5)
Figure 4

Criteria weight with no weight information.

Figure 4

Criteria weight with no weight information.

Close modal
The probability density function of the indicator evaluation and the probability density function of the indicator weights are known, and furthermore, SMAA-2 obtains the combined utility of each scheme by weighting and summing the utility values of each attribute through the linear utility function of Equation (6). The schemes are then ranked in order of their advantages and disadvantages based on , and the scheme that satisfies the decision-making requirements is selected from them.
(6)
where a random variable can be used instead of the constant to reflect the uncertainty of the indicator value, then Equation (6) becomes .
The scheme ranking function is defined as shown in Equation (7):
(7)
where the range of values of the scheme ranking function is ; and are 1 and 0, respectively.
The ranking propensity weight is then defined. SMAA-2 defines the space in which the scheme obtains a ranking of as the ranking propensity weight , as shown in Equation (8):
(8)
The ranking acceptability metric is defined as shown in Equation (9). It represents the acceptability of alternative ranked at the rth position, which can also be regarded as the probability that alternative is ranked at the rth position.
(9)
Define the global acceptability level . It is a synthesis of all rankings obtained for the scheme and describes the overall acceptable level of the scheme as a whole.
(10)
where indicator takes a value in the range of and is a secondary weight that indicates the contribution of a ranking r of the scheme to the indicator .
Commonly used secondary weights include linear weights , inverse weights , and centroid weights , where the smaller the weight r, the larger the corresponding secondary weight, indicating a greater emphasis on acceptability at the top of the rankings. Figure 6 depicts the three different secondary weights described above (with a number of schemes of 10 as an example). As can be seen in Figure 5, all three secondary weights give higher weights to the top-ranked schemes, but the linear weights give greater weight to the centrally ranked schemes compared to the inverse and centroid weights. Compared to the centroid weights, the inverse weights have a more even distribution of weights for the schemes at the bottom of the ranking. As a result, the inverse weights for aggregation are less sensitive to the schemes that are ranked lower than the centroid weights. Also, Barron & Barrett (1996) pointed out that when normalizing the secondary weights for multi-criteria decision-making, compared to other forms of weights based on ranking, the center of gravity weights contain valid information about other forms of weights and ranking, and it is more accurate and effective to use the centroid weights as secondary weights. Therefore, in this paper, the centroid weights are selected as the secondary weights in the multi-criteria decision-making of wind–photovoltaic–hydropower systems.
Figure 5

Schematic of meta-weight.

Figure 5

Schematic of meta-weight.

Close modal
Figure 6

A set of alternatives for WPHSs for multi-criteria decision-making.

Figure 6

A set of alternatives for WPHSs for multi-criteria decision-making.

Close modal
The central weight vector is the center of the set of weights corresponding to the ranking r obtained by the scheme , which can be obtained by binary integration of the probability density functions of the evaluation indicators and indicator weights, as shown in Equation (11):
(11)
It can be seen from Equation (11) that when , the center weight vector is the center of the space of optimal ranking weights obtained by the scheme , therefore, the center weight vector responds to the decision-maker's preference information for the scheme .

VIKOR model

The VIKOR method is a multi-criteria decision-making method for trade-off ranking. It starts by identifying a set of positive and negative ideal solutions, and subsequently calculates the distances between the alternatives and the positive and negative ideal solutions, and performs a trade-off ranking of the alternatives based on maximizing the group benefits and minimizing the individual regrets (Opricovic & Tzeng 2004). Both the VIKOR method and the TOPSIS method are trade-off ranking methods for near-ideal solutions. The TOPSIS method bases the trade-off ranking on the fact that the alternatives are the closest to the positive ideal solution and the farthest from the negative ideal solution. This may lead to a reverse order result, while the VIKOR method does not need to consider the problem that the closest alternative needs to be closest to the ideal point and farthest from the negative ideal point, which is a good way to avoid the reverse order problem.

Group decision-making for wind–photovoltaic–hydropower systems

A linear utility function is used in the original SMAA-2 model to evaluate the utility values of the alternatives. Due to the characteristic of non-metricity of multi-objective problems in wind–photovoltaic–hydropower systems, the decision-making directly through the linear sum-type utility function in the SMAA-2 model may not be able to obtain a fair solution. However, the SMAA-2 model has good expandability and is easy to be coupled with other multi-attribute decision-making methods (Corrente et al. 2014). And the VIKOR model is widely used in the field of multi-attribute decision-making, which has been proved to have better maneuverability, convenience and robustness by calculating the proximity of the evaluated value of each alternative to the ideal solution and avoiding the inverse order problem. Therefore, we try to review the VIKOR model into the SMAA-2 model of inverse weight space analysis, and propose a stochastic multi-criteria decision-making method, i.e., the SMAA–VIKOR model.

The SMAA–VIKOR model inherits the SMAA-2 model with the characteristic of inverse weight space analysis, and the specific operation procedure is as follows:

  • (1)

    Normalize the decision matrix of the alternatives to obtain a standardized decision matrix ;

  • (2)

    Determine the distribution of the feasible weight space and perform Latin Hypercube Sampling (LHS) based on its probability density function to randomly generate feasible weights w;

  • (3)

    Call the VIKOR model and calculate to get the corresponding ranking of each scheme;

  • (4)

    Determine whether the number of iterations is satisfied; if so, next step; if not, return to step (2);

  • (5)

    Calculate the ranked acceptability indicators and global acceptability indicators and center weight vectors based on the ranking of the alternatives obtained during the cycle;

  • (6)

    End.

The proposed methodology is applied to the Yalong River Basin. The experimental runoff data were measured at the hydrological station in 2016. The wind speed, solar radiation, and temperature were obtained from National Meteorological Administration (http://data.cma.cn/data/detail/dataCode/A.0012.0001.html). A more detailed description of the study area and data can be found in our previous study (Liu et al. 2019). In this study, on the basis of the non-inferior solution set of the wind–photovoltaic–hydropower system obtained after the multi-objective optimization by Liu et al. (2019), 20 representative schemes were selected uniformly on the Pareto frontier surface obtained in that study, as shown in Figure 6. In order to efficiently solve the multi-criteria decision model, the Monte Carlo (MC) and LHS were compared for the SMAA-2 and SMAA–VIKOR models in Section 3.1. In Section 3.2, the SMAA-2 and SMAA–VIKOR models were compared to verify the validity of the VIKOR method as a utility function for the decision model over the simple linear summation type utility function. Finally, the SMAA–VIKOR was compared with the VIKOR in Section 3.3 to confirm the superiority of SMAA–VIKOR over the deterministic VIKOR method.

Sampling method

The SMAA-2 model involves a number of multiple integration operations, and the dimension of its integrals is high, which may be difficult to solve analytically if it is performed directly. The original SMAA-2 model is solved by MC random sampling, which may occur when samples that have already been taken are taken again, reducing the sampling efficiency. The LHS is a stratified sampling method in which the sampling points are distributed throughout the sampling area to avoid drawing samples that are already present. In order to compare the efficiency difference between the two sampling methods and to validate SMAA–VIKOR, five sets of comparative experiments were conducted for the SMAA-2 and SMAA–VIKOR models under the two sampling methods, with 100, 500, 1,000, 5,000, and 10,000 simulations in each set, respectively. In order to eliminate the errors due to the randomness of the experiments, each group of experiments was independently repeated 50 times. The results are shown in Figure 7 (horizontal coordinate MC-100 refers to the number of times the SMAA-2 or SMAA–VIKOR model was simulated in the case of MC random sampling as 100, and so on, and LHS-100 refers to the number of times the SMAA-2 or SMAA–VIKOR model was simulated in the case of Latin Hypercubic Random Sampling as 100). As can be seen in Figure 7(a), the SMAA-2 model has a smaller range of global acceptability indicator distributions obtained from LHS than MC random sampling for all five sets of comparison experiments. As shown in Figure 7(b), the range of the distribution of global acceptability indicators obtained by the SMAA–VIKOR model with LHS is also all smaller than MC method. This indicates that LHS has better stability in solving both SMAA-2 and SMAA–VIKOR models.
Figure 7

Comparison of MC and LHS.

Figure 7

Comparison of MC and LHS.

Close modal

Comparative analysis of utility functions

In order to verify the effectiveness of the improved utility function and to validate SMAA–VIKOR, the effectiveness of the SMAA-2 and SMAA–VIKOR models in the unweighted preference information improvement method was compared. Both of the above models were solved using LHS. Figure 8 demonstrates the global acceptability metrics based on the SMAA-2 and SMAA–VIKOR models with unknown weight information. From the global acceptability metrics, it can be seen that the overall ranking of each scheme obtained based on the SMAA-2 and SMAA–VIKOR models is similar, with the top two ranked schemes being A2 and A3, and the last ranked scheme being A12.
Figure 8

Holistic acceptability metrics of SMAA-2 and SMAA–VIKOR models under no weight information.

Figure 8

Holistic acceptability metrics of SMAA-2 and SMAA–VIKOR models under no weight information.

Close modal
Figure 9 further illustrates the ranking acceptability metrics obtained by the SMAA-2 and SMAA–VIKOR models. It can be seen that both models yield the highest probability of ranking first for Scheme 2, and the probability of ranking first for Scheme 2 obtained by the SMAA–VIKOR model is 0.69, which is greater than the probability of ranking first for Scheme 2 obtained by the SMAA-2 model, which is 0.46. At the same time, the ranking of each scheme obtained based on the SMAA–VIKOR model is more centralized than that obtained by the SMAA-2 model. This can more effectively reduce the impact of weight uncertainty on the results of stochastic decision-making, thus obtaining more obvious probabilistic ranking results and providing more explicit information for decision-makers.
Figure 9

Ranking of acceptable metrics of SMAA-2 and SMAA–VIKOR models under no weight information.

Figure 9

Ranking of acceptable metrics of SMAA-2 and SMAA–VIKOR models under no weight information.

Close modal

Comparative analysis of deterministic and stochastic multi-criteria decision-making

The SMAA–VIKOR was compared with VIKOR in this section to validate SMAA–VIKOR and to further confirm the superiority of SMAA–VIKOR over the deterministic VIKOR. Since the deterministic VIKOR model has to be input with indicator weight information for multi-criteria decision-making, in this section, experts are invited to give the corresponding preference information for the three indicators, as shown in Table 1. The three indicators correspond to the three objectives of multi-objective optimization, with weights of 0.26, 0.36, and 0.38, respectively.

Table 1

Deterministic weight information of three indicators

IndicatorPower generationMinimum output for the time periodAPFD
Indicator weight 0.26 0.36 0.38 
IndicatorPower generationMinimum output for the time periodAPFD
Indicator weight 0.26 0.36 0.38 

By inputting the indicator weight information in Table 1 into the deterministic VIKOR model for decision-making, the values of the integrated indicators can be obtained as shown in Table 2. As can be seen from Table 2, the comprehensive ranking structure of each scheme is A2A6A3A5A7A1A9A10A8A4A11A12.

Table 2

Comprehensive indicators and their corresponding rankings under the deterministic VIKOR model for each scheme

Scheme numberA1A2A3A4A5A6A7A8A9A10A11A12
Integrated indicator 0.38 0.00 0.20 0.71 0.31 0.15 0.36 0.61 0.46 0.52 0.79 0.93 
Rank 10 11 12 
Scheme numberA1A2A3A4A5A6A7A8A9A10A11A12
Integrated indicator 0.38 0.00 0.20 0.71 0.31 0.15 0.36 0.61 0.46 0.52 0.79 0.93 
Rank 10 11 12 

Less information is available to the decision-maker at the early stage of decision-making, at which time the weighting information can be considered completely unknown. Using the proposed SMAA–VIKOR model with LHS and through inverse weight space analysis, the global acceptability indicator of each scheme can be obtained, as shown in Figure 10. As can be seen from Figure 10, the ordering of the schemes obtained based on the SMAA–VIKOR model is A2A3A6A5A1A4A7A9A10A11A8A12, which is similar to the ordering of the deterministic VIKOR model.
Figure 10

Holistic acceptability indices of the SMAA–VIKOR model under no weight information.

Figure 10

Holistic acceptability indices of the SMAA–VIKOR model under no weight information.

Close modal
Figure 11 further illustrates the acceptability metrics for the ranking of each scenario obtained based on the SMAA–VIKOR model. As can be seen from Figure 11, each scheme has a certain probability of getting a different ranking, and the probability of each ranking is not equal. Scheme A2 has the highest probability of obtaining the first place in the ranking, which is 0.69, but scheme A2 also has a certain probability of obtaining other rankings, such as the fourth place with a probability of 0.13. Similarly, scheme A12, which has the worst ranking obtained through the global ranking acceptability metrics of Figure 10, also has a probability of obtaining the first place in the ranking with a probability of 0.017. The centrally ranked schemes A1, A4, and A7 are obtained through the global ranking acceptability metrics in Figure 10, and they have probabilities of 0.227, 0.211, and 0.069, respectively, for obtaining the first place in the ranking, and the probabilities of 0.084, 0.020, and 0.117, respectively, for obtaining the sixth place in the ranking.
Figure 11

Ranking of acceptable metrics of the SMAA–VIKOR model under no weight information.

Figure 11

Ranking of acceptable metrics of the SMAA–VIKOR model under no weight information.

Close modal
Figure 12 illustrates the central weight vector of each scheme when it obtains the optimized ranking for the SMAA–VIKOR model with unknown weight information. It can be seen that the corresponding typical weight vectors are different when each scheme obtains the optimized ranking. Taking Scheme A3 as an example, from Figure 6, it can be seen that Scheme A3 has larger values of power generation and minimum output during the time period, and if Scheme A3 obtains the optimized ranking, the two indicators of power generation and minimum output during the time period need to obtain larger weights, which are 0.41 and 0.37, respectively, and at this time the corresponding ecological indicator APFD has smaller weights.
Figure 12

Center weight vectors of the SMAA–VIKOR model under no weight information.

Figure 12

Center weight vectors of the SMAA–VIKOR model under no weight information.

Close modal

With the above results, it can be seen that the SMAA–VIKOR model can obtain the probability of each scheme to achieve different ranking results, while the deterministic VIKOR model can only obtain the unique ranking results. In addition, the deterministic VIKOR model requires the decision-makers to give decision information before the multi-criteria decision-making, which creates a big difficulty for the decision-makers. which does not require information about the weights of the indicators given by the decision-making group, has a certain superiority over the deterministic VIKOR model and can provide more useful information for the decision-makers when making decisions.

In group decision-making for scheduling schemes of wind–photovoltaic–hydropower systems, there exists a certain uncertainty in the indicator weights or even no information about the indicator weights due to the knowledge level of the decision-makers and the attribute information of the scheme set itself. This study proposes a stochastic multi-criteria decision-making framework for wind–photovoltaic–hydropower systems to overcome the difficulty of uncertainty in indicator weights or even completely unknown information of indicator weights when making decisions. Under this framework, SMAA theory and VIKOR model are reviewed, and the SMAA–VIKOR model is proposed to clarify the indicator weight space. The study shows that the proposed SMAA–VIKOR model can overcome the obstacle of decision-makers' lack of information on indicator weights. The ranking acceptability indicators calculated by the model show a more obvious trend of advantages and disadvantages, which gives full confidence to the decision-making group to formulate a plan to be implemented. It breaks through the bottleneck of group decision-making, which is difficult to make effective decisions due to the condition of incomplete information, and enriches the library of stochastic multi-criteria decision-making methods for the scientific formulation of scheduling schemes of wind–photovoltaic–hydropower systems under uncertainty conditions. Several findings can be revealed as follows:

In solving the SMAA-2 and SMAA–VIKOR models, the global acceptability indicators obtained by the LHS method are generally better than those obtained by the MC method, and the stability of the LHS method is also better than that of the MC method.

The SMAA–VIKOR model has a probability of 0.69 to rank first globally for Scheme A2, which is higher than the probability of 0.46 for the SMAA-2 model to rank first for Scheme A2. Additionally, the schemes obtained from the SMAA–VIKOR model are more concentrated in their rankings compared to the SMAA-2 model. This concentration can effectively reduce the impact of weight uncertainty on random decision results in wind–photovoltaic–hydropower systems decision-making, resulting in more distinct probability ranking outcomes and providing decision-makers with clearer information.

In contrast to the deterministic VIKOR model, the SMAA–VIKOR model, which considers the uncertainty of indicator weights, provides varying probabilities for each scheme across different rankings. Scheme A2 has the highest probability of ranking first at 0.69 but also holds probabilities for other rankings, such as 0.13 for ranking fourth. The scheme with the worst ranking acceptability index, A12, has a probability of 0.017 to rank first. Schemes A1, A4, and A7, which rank in the middle, have probabilities of 0.227, 0.211, and 0.069 to rank first, and probabilities of 0.084, 0.020, and 0.117 to rank sixth, respectively. Unlike the deterministic VIKOR model, which provides only one ranking for each scheme, the SMAA–VIKOR model does not require decision-makers to provide decision information in advance before multi-attribute decision-making. This eliminates the difficulty decision-makers face in determining wind-hydro system scheduling plans at the initial stages of decision-making. The SMAA–VIKOR model, which does not require decision groups to provide indicator weight information, has certain advantages and can offer decision-makers more useful information during the decision-making process.

The complementary characteristics of wind, photovoltaic, and hydro power generation exist at different time scales, and the uncertainty of prediction presents a complex coupling relationship. Conducting research on the prediction and uncertainty analysis of basin wind, photovoltaic, and hydro power generation is of great significance for elucidating the coupling laws of complementary relationship and prediction uncertainty. This study focuses on the deterministic input used in the long-term scheduling of basin wind, photovoltaic, and hydro power systems, without involving the study of long-term prediction and uncertainty of wind and photovoltaic energy. Future research needs to start from the long-term prediction uncertainty of basin wind, photovoltaic, and hydro power, introduce prediction models for the coupling relationship of long-term prediction uncertainty in basin wind, photovoltaic, and hydro power generation, and use them to scientifically formulate long-term scheduling plans for basin wind, photovoltaic, and hydro power systems.

W.L., Y.Z., Y.L. conceptualized the study; W.L., Y.Z., X.X. performed the methodology; W.L., X.G., R.M., Y.Z. did formal analysis and investigation; W.L., Y.Z. wrote and prepared the original draft; W.L., X.G., R.M, Y.L., Y.Z. wrote, reviewed, and edited the article; Y.L., Y.Z., W.L. acquired the funds; W.L., Y.Z. collected resources; Y.L. supervised the article.

This work was supported by the National Key Research and Development Program of China (Grant No. 2022YFC3202300); the National Natural Science Foundation of China (Grant No. 52209032); the China Postdoctoral Science Foundation Funded Project (Grant No. 2021M702313).

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Barron
F. H.
&
Barrett
B. E.
1996
Decision quality using ranked attribute weights
.
Manage. Sci.
42
(
11
),
1515
1523
.
https://doi.org/10.1287/mnsc.42.11.1515
.
Chen
Y. B.
,
Wang
M. L.
,
Zhang
Y.
,
Lu
Y.
,
Xu
B.
&
Yu
L.
2023
Cascade hydropower system operation considering ecological flow based on different multi-objective genetic algorithms
.
Water Resour. Manage.
37
,
3093
3110
.
https://doi.org/10.1007/s11269-023-03491-3
.
Corrente
S.
,
Figueira
J. R.
&
Greco
S.
2014
The SMAA-PROMETHEE method
.
Eur. J. Oper. Res.
239
(
2
),
514
522
.
https://doi.org/10.1016/j.ejor.2014.05.026
.
Kang
H. Y.
,
Hung
M. C.
,
Pearn
W. L.
,
Lee
A. H. I.
&
Kang
M. S.
2011
An integrated multi-criteria decision-making model for evaluating wind farm performance
.
Energies
4
(
11
),
2002
2026
.
https://doi.org/10.3390/en4112002
.
Lahdelma, R. & Salminen, P. 2001 SMAA-2: Stochastic multicriteria acceptability analysis for group decision making. Operations Research 49 (3), 444–454. https://doi.org/10.1287/opre.49.3.444.11220
.
Liu
W. F.
,
Zhu
F. L.
,
Chen
J.
,
Wang
H.
,
Xu
B.
,
Song
P. B.
,
Zhong
P. A.
,
Lei
X. H.
,
Wang
C.
,
Yan
M. J.
,
Li
J. Y.
&
Yang
M. Z.
2019
Multi-objective optimization scheduling of wind–photovoltaic–hydropower systems considering riverine ecosystem
.
Energ. Convers. Manage.
196
,
32
43
.
https://doi.org/10.1016/j.enconman.2019.05.104
.
Liu
W. F.
,
Zhu
F. L.
,
Zhao
T. T.
,
Wang
H.
,
Lei
X. H.
,
Zhong
P. A.
&
Fthenakis
V.
2020
Optimal stochastic scheduling of hydropower-based compensation for combined wind and photovoltaic power outputs
.
Appl. Energy
276
,
115501
.
https://doi.org/10.1016/j.apenergy.2020.115501
.
Lu
Y. L.
,
Zhou
J. Z.
,
Qin
H.
,
Wang
Y.
&
Zhang
Y. C.
2011
Environmental/economic dispatch problem of power system by using an enhanced multi-objective differential evolution algorithm
.
Energ. Convers. Manage.
52
(
2
),
1175
1183
.
https://doi.org/10.1016/j.enconman.2010.09.012
.
Opricovic
S.
&
Tzeng
G. H.
2004
Compromise solution by MCDM methods: A comparative analysis of VIKOR and TOPSIS
.
Eur. J. Oper. Res.
156
(
2
),
445
455
.
https://doi.org/10.1016/s0377-2217(03)00020-1
.
Perera
A. T. D.
,
Attalage
R. A.
,
Perera
K. K. C. K.
&
Dassanayake
V. P. C.
2013
A hybrid tool to combine multi-objective optimization and multi-criterion decision-making in designing standalone hybrid energy systems
.
Appl. Energy
107
,
412
425
.
https://doi.org/10.1016/j.apenergy.2013.02.049
.
Qin
H.
,
Zhou
J. Z.
,
Lu
Y. L.
,
Wang
Y.
&
Zhang
Y. C.
2010
Multi-objective differential evolution with adaptive Cauchy mutation for short-term multi-objective optimal hydro-thermal scheduling
.
Energ. Convers. Manage.
51
(
4
),
788
794
.
https://doi.org/10.1016/j.enconman.2009.10.036
.
Rigatos
G.
,
Siano
P.
,
Abbaszadeh
M.
&
Wira
P.
2019
Nonlinear optimal control for wind power generators comprising a multi-mass drivetrain and a DFIG
.
J. Franklin Inst.
356
(
5
).
https://doi.org/10.1016/j.jfranklin.2018.12.017
.
Russo
M. A.
,
Carvalho
D.
,
Martins
N.
&
Monteiro
A.
2023
Future perspectives for wind and solar electricity production under high-resolution climate change scenarios
.
J. Cleaner Prod.
404
,
136997
.
https://doi.org/10.1016/j.jclepro.2023.136997
.
Sanchez-Lozano
J. M.
,
Garcia-Cascales
M. S.
&
Lamata
M. T.
2016
GIS-based onshore wind farm site selection using fuzzy multi-criteria decision-making methods. Evaluating the case of southeastern Spain
.
Appl. Energy
171
,
86
102
.
https://doi.org/10.1016/j.apenergy.2016.03.030
.
Tan
Q. F.
,
Wen
X.
,
Sun
Y. L.
,
Lei
X. H.
,
Wang
Z. N.
&
Qin
G. H.
2021
Evaluation of the risk and benefit of the complementary operation of the large wind–photovoltaic–hydropower system considering forecast uncertainty
.
Appl. Energy
285
,
116442
.
https://doi.org/10.1016/j.apenergy.2021.116442
.
Wang
M. L.
,
Zhang
Y.
,
Lu
Y.
,
Wan
X. Y.
,
Xu
B.
&
Yu
L.
2022
Comparison of multi-objective genetic algorithms for optimization of cascade reservoir systems
.
J. Water Clim. Change
13
(
11
),
4069
4086
.
https://doi.org/10.2166/wcc.2022.290
.
Xia
L.
,
Chen
G.
,
Wu
T.
,
Gao
Y.
,
Mohammadzadeh
A.
&
Ghaderpour
E.
2023
Optimal intelligent control for doubly Fed induction generators
.
Mathematics
11
,
20
.
https://doi.org/10.3390/math11010020
.
Zhang
Y.
,
Yu
L.
,
Wu
S. Q.
,
Wu
X. F.
,
Dai
J. Y.
,
Xue
W. Y.
&
Yang
Q. Q.
2021
A framework for adaptive control of multi-reservoir systems under changing environment
.
J. Cleaner Prod.
316
,
128304
.
https://doi.org/10.1016/j.jclepro.2021.128304
.
Zhu
F. L.
,
Zhong
P. A.
,
Wu
Y. N.
,
Sun
Y. M.
,
Chen
J.
&
Jia
B. Y.
2017
SMAA-based stochastic multi-criteria decision-making for reservoir flood control operation
.
Stochastic Environ. Res. Risk Assess.
31
(
6
),
1485
1497
.
https://doi.org/10.1007/s00477-016-1253-3
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).