Multi-Objective Optimal Scheduling for Multi-Renewable Energy Power System Considering Flexibility Constraints

: As renewable energy penetration increases, the lack of ﬂexibility in a multi-renewable power system can seriously affect its own economics and reliability. To address this issue, three objectives are considered in this study: power ﬂuctuations on tie-line, operating cost, and curtailment rate of renewable energy. Presented also is an optimal day-ahead scheduling model based on the MREPS for distributed generations with ﬂexibility constraints. The multi-objective particle swarm optimization (MOPSO) algorithm can be applied to obtain a set of Pareto non-dominated solutions for the day-ahead scheduling strategy with the proposed model. By using fuzzy comprehensive evaluation, the optimal compromise solution is determined in the set. The presented method sacriﬁces a small amount of economy and power ﬂuctuation, but it can reduce the deviation between forecast and realized power ﬂuctuations on the tie-line, while improving the utilization of renewable energy.


Introduction
In modern power systems, the scarcity of fossil fuels and increasing pollution of the environment contribute to the development of renewable energy sources, such as solar and wind. Despite this, the stochastic nature of renewable energy generation is likely to have significant effects on system reliability and economy [1][2][3].
When operating a multi-renewable energy power system (MREPS), it is necessary to develop an optimal schedule to cope with the stochasticity of renewable energy generation [4,5]. MREPS scheduling is divided into two categories: day-ahead and real-time scheduling. There are several strategies for achieving various operational objectives using day-ahead scheduling. Reference [6] used a two-stage stochastic optimization model for an MREPS for minimizing the short-term operation cost, which introduces uncertainty in renewable generation, and showed that stochastic scheduling can provide significant reliability benefits to multi-energy supply systems. Reference [7] presented a day-ahead scheduling model that considers the seasonal uncertainty of renewable energy for a microgrid equipped with multi-renewable energy units. An improved optimization algorithm was proposed to solve optimization problems that focus on minimizing the operation cost.
The results indicate that day-ahead scheduling based on the proposed algorithm can provide an efficient solution for managing MREPS energy. An MREPS is highly influenced by both economic indicators and the rate at which renewable energy is utilized, as illustrated in [8][9][10]. Furthermore, power fluctuations on the tie-line serve as relevant indicators for the main grid connected to an MREPS [11,12]. In the current MREPS scheduling process, the impact of these power fluctuations on the main grid is rarely considered. Although operating costs, curtailment rates of renewable energy, and power fluctuations on the tie-lines are all commonly considered in current studies on optimal dispatch [4][5][6][7][8][9][10][11][12][13][14][15], few studies consider these three indicators simultaneously. This can lead to a situation where, while some of these indicators are optimal, the other indicators may be poor. As a result, the present study proposes a day-ahead scheduling strategy based on operating costs, curtailment rates of renewable energy, and power fluctuations on the tie-lines.
In the case of real-time scheduling, the schedule is amended if the actual renewable energy output differs from the forecasted output. An improved particle swarm-optimization (PSO)-based strategy for managing energy over a two-time scale was presented in [16]. In day-ahead scheduling, one of the primary objectives is to achieve the most cost-effective schedule, and in real-time scheduling, the primary goal is to track day-ahead scheduling, compensate for power fluctuations, and maintain system stability. According to the experimental results, the proposed method could minimize the cost of generated electricity and maximize the efficiency of renewable energy systems. Real-time scheduling, however, relies on day-ahead scheduling. It is possible to increase the utilization rate of renewable energy if uncertainties associated with the generation of renewable energy can be properly incorporated into day-ahead scheduling. To address uncertainties in power systems [17][18][19], flexibility has been proposed. In [20], an optimal scheduling model for flexible resources was presented from both the generation and load sides. According to the results, a dynamic line rating model that incorporates optimal scheduling can maximize the utilization of flexible resources without curtailing wind power and minimize dispatch costs. Consequently, we propose a day-ahead scheduling strategy for an MREPS that accounts for flexibility constraints and concentrates on operation costs, renewable energy curtailment rates, and power fluctuations on the tie-line.
In this study, we attempt to solve a multi-objective optimization problem. It is possible to solve a multi-objective optimization problem in several ways. Among them, the linear weighted sum method and intelligent algorithms are frequently adopted. The linear weighted sum method is simple and fast, but it usually gives unclear physical results. In addition, the resulting error is typically large because several targets are of different dimensions and orders of magnitude. These factors remarkably impact the results and conclusions of practical problems [21]. Intelligent algorithms, which have clearer physical meaning, outperform the linear weighted sum method in accuracy, flexibility, and effectiveness in solving multi-objective problems. A number of intelligent algorithms have been successfully applied to engineering optimization, such as multi-objective particle swarm optimization (MOPSO) due to its unique search mechanism, excellent convergence performance, and convenient calculation capabilities [22,23]. As a result, MOPSO is adopted in this study to identify the Pareto non-dominated set of objective functions. Following the determination of the Pareto non-dominated solution set, a fuzzy comprehensive evaluation method [24] is adopted to determine the optimal compromise solution. The optimal day-ahead scheduling strategy can be determined based on the optimum compromise solution.
The contributions of this paper are summarized as follows: 1 Considering the operation cost, renewable energy curtailment rates, and power fluctuations on the tie-line, a day-ahead scheduling model for the MREPS is established. 2 MOPSO and a fuzzy comprehensive evaluation method are used to evaluate the day-ahead scheduling model, and a day-ahead scheduling strategy for the MREPS considering flexibility is proposed.
Following is the remainder of this paper. The MREPS presents the day-ahead optimal scheduling model in Section 2, along with its constraints, taking flexibility into account. In Section 3, MOPSO and fuzzy comprehensive evaluation are discussed. An analysis of the experimental results is presented in Section 4, which simulates an actual MREPS. In Section 5, the conclusions are summarized.

Model for Multi-Objective Optimal Scheduling
The MREPS, consisting of wind turbines (WTs), photovoltaic (PV) arrays, diesel generators (DGs), energy storage systems (ESSs), and loads, only purchases electricity from the main grid. Figure 1 illustrates the details of the MREPS.
Following is the remainder of this paper. The MREPS presents the day-ahead scheduling model in Section 2, along with its constraints, taking flexibility into In Section 3, MOPSO and fuzzy comprehensive evaluation are discussed. An an the experimental results is presented in Section 4, which simulates an actual MR Section 5, the conclusions are summarized.

Model for Multi-Objective Optimal Scheduling
The MREPS, consisting of wind turbines (WTs), photovoltaic (PV) array generators (DGs), energy storage systems (ESSs), and loads, only purchases el from the main grid. Figure 1 illustrates the details of the MREPS.

Operation Cost
A major consideration in the MREPS's day-ahead scheduling is the operat The operation cost mainly includes fuel cost, operation and maintenance, p control, and purchased power cost. The daily operation cost of the MREPS expressed as follows: where Cfuel, COM, and CPC are the total fuel, operation and maintenance (OM), and p control (PC) costs associated with each distributed power generation in one day is the total cost of the purchased power (PP) from the main grid in one day. The f for Cfuel, COM, CPC, and CPP are as follows:

Operation Cost
A major consideration in the MREPS's day-ahead scheduling is the operation cost. The operation cost mainly includes fuel cost, operation and maintenance, pollutant control, and purchased power cost. The daily operation cost of the MREPS can be expressed as follows: where C fuel , C OM , and C PC are the total fuel, operation and maintenance (OM), and pollutant control (PC) costs associated with each distributed power generation in one day and C PP is the total cost of the purchased power (PP) from the main grid in one day. The formulas for C fuel , C OM , C PC , and C PP are as follows: where N is the total number of DGs; T is the total scheduling time; ξ f.DG.n is the consumption cost per kW·h of the nth DG; P DG.n (t) is the output power of the nth DG in period t; ∆t is the scheduling interval; M is the total number of distributed generation types; ξ OM.m is the cost of operation and maintenance per kW·h of the mth type of distributed generation variety; P m (t) is the output power of the mth type of distributed generation during period t; J is the total number of pollutant types; ξ m.j is the emission cost factor for the jth type of pollutants generated by the mth type of distributed power generation; ξ PP (t) and P PP (t) are the power purchased by the MREPS from the main grid for the period t, respectively.

Renewable Energy Curtailment Rate
The renewable energy curtailment rate is the ratio of the power curtailed by all renewable energy units to the total power that can be generated by renewable energy units during the dispatching period. The lower the value, the higher the utilization rate of renewable energy is. Utilization of the renewable energy is ideally increased while reducing the operation cost of the MREPS [15]. Thus, this study considers the curtailment rate of renewable energy. Following is the daily curtailment rate of the MREPS's renewable energy: where H represents the total number of the renewable energy types; P h.c (t) and P h.all (t) reflect the power curtailment and the available power generation of the hth type of renewable energy in period t, respectively.

Tie-Line Power Fluctuations
The frequent power exchange between the MREPS and the main grid can be attributed to the stochastic nature of renewable energy. In general, the unit commitment and economic load dispatch of the main grid have no impact on the MREPS. This is because the MREPS can be self-sufficient. Only when the MREPS is not self-sufficient, it is necessary for the main grid to help the MREPS to achieve power balance through unit commitment or economic load dispatch. As a consequence, the significant fluctuation in power on the tie-line has a negative impact on the main grid. The MREPS day-ahead scheduling must be designed to take into account power fluctuations across the tie-line.
As the square root of the variance, the standard deviation is a metric used to measure how dispersed a dataset is relative to its mean. Consequently, the standard deviation is typically used to describe the fluctuation of a data series. However, the standard deviation lacks comparability for different objects or samples with varying means of the same object. Hence, the coefficient of variation is used in this study to avoid these problems.
As a measure of data dispersion around the mean, the coefficient of variation represents the ratio of the standard deviation to the mean for a series of data points. The degree of variation from one data series to another can be compared, although the means are remarkably different from one another. A small variation indicates a small fluctuation degree [11]. On the tie-line, the power fluctuation can be defined as follows: where P TL (t) represents the transmission power of the tie-line during period t and µ TL represents the average transmission power of the tie-line during one day, which can be expressed as follows: Based on the above three objectives, the total objective function for the multi-objective optimal scheduling model is as follows: During the operation of the MREPS, at the end of each period, the sum of the output power of the DGs, the charging and discharging power of the ESS, and the purchased power from the main grid should be equal to the net load (NL) power of the system. As a result of this relationship, we may state: where P DG (t) is the sum of the output power of the DGs in period t; P ESS (t) is the charge and discharge power of the ESS in period t and becomes negative when the ESS is being charged; P TL (t) is the transmission power of the tie-line in period t; P L (t), P WT (t), P PV (t), and P NL (t) are the load power, output of WTs, output of PV arrays, and net load power in period t, respectively.

Constraints of Output Power for DGs Considering Flexibility
In light of the uncertainty associated with renewable energy generation, it is imperative that the MREPS be flexible. The DG is a commonly used flexible resource in MREPSs because of its satisfactory controllability and fast response speed. Therefore, flexibility should be considered and exploited in the output power constraints of DGs.
The uncertainty of renewable energy generation leads to prediction errors in net load power and the MREPS's flexibility requirements. These flexibility requirements are provided by DGs in the MREPS. DGs are typically restricted in their output power range as follows [7]: where P DG.min and P DG.max represent the minimum and maximum output power of DGs, respectively. In terms of conventional constraints, DGs are only considered for their minimum and maximum output. Due to these constraints, the output power of DGs may be inflexible during operation. When large prediction errors occur, DGs output power fails to meet the flexibility requirements. Therefore, new constraints that consider flexibility based on conventional constraints can be expressed as follows: where F NL.U (t) and F NL.D (t) represent the maximum upward and downward flexibility requirements of the MREPS in each period t.
Because of the error in predicting net load power, the MREPS is required to be flexible. Variables such as wind, solar, and load power are among the factors that affect the prediction errors of the net load power. As a result, wind, solar, and load power prediction errors are assumed to follow a normal distribution with a zero mean [25]. σ PV (t), σ WT (t), and σ L (t) have the following standard deviations: where W PV.F (t), W WT.F (t), and W L.F (t) are the predicted wind, solar, and load power in period t, respectively, and W PV.C and W WT.C are the total installed capacity of PV arrays and WTs, respectively. There are no correlations among wind, solar, and load power prediction errors. According to the basic characteristics of the normal distribution, the sum of two independent normally distributed random variables is normal, with the sum of the two means as the mean and the sum of the two variances as the variance. Hence, prediction errors of the net load power also follow the normal distribution with mean zero and standard deviation σ NL (t), as shown below: The MREPS must offer flexibility equal to the forecast errors of net load power. However, in the prediction of net load power, there is a very small probability of maximum errors. Hence, the scenario where the flexibility requirement and forecast error are equal is uneconomical and wasteful. As a result, we need to consider the confidence intervals of these prediction errors. If the confidence level is 1 − α, then F NL.U (t) and F NL.D (t) can be expressed as: where u 1−α/2 is the upper (1 − α/2) critical value for the standard normal distribution. Conventional ramping-rate constraints of DGs are expressed as follows: where R D DG.max and R U DG.max are the maximum downward and upward ramping rates of DGs, respectively, with both values positive, and ∆t is the time interval.
Conventional constraints only consider the minimum and maximum ramping rates. The maximum ramping rate of DGs may occur in some periods under these constraints. The output power of DGs fails to meet the fluctuation of the net load power in this situation when DGs lack the ramping rate. Consequently, ramping-rate constraints should take flexibility into consideration.
As shown in Figure 2, A 0 , B 0 , and C 0 are the prediction output power values of DGs at different times. DGs must reduce their output power to A 2 if the maximum downward prediction errors of the net load power occur at time t − ∆t. According to Equation (15), the reduced output power should be equal to F NL.D (t − ∆t). It is imperative that the DGs increase their output power to B 1 if they experience maximum upward prediction errors at time t. According to Equation (15), the increased output power should be equal to F NL.U (t). Therefore, the margin of the ramping rate that DGs should reserve is the sum of F NL.D (t − ∆t) and F NL.U (t) from time t − ∆t to time t. Similarly, the margin of the ramping rate that DGs should reserve is the sum of F NL.U (t) and F NL.D (t + ∆t) from time t to time t + ∆t. As a result, new constraints that consider flexibility can be expressed as follows:

Constraints of the ESS Considering Flexibility in Charging and Dischargin
In this study, the main function of the ESS is to provide assistance to DGs an that DGs have adequate flexibility. Thus, the ESS and its conventional constr similar in terms of charging and discharging power: where SESS(t) represents the state of charge (SOC) of the ESS in period t; δ repre self-discharge efficiency of the ESS; EESS represents the maximum capacity of SESS.min and SESS.max represent the minimum and maximum SOC of the ESS, respect represents an hour; ηESS represents the percentage of the maximum charge-d capacity per hour to the maximum capacity of the ESS; αESS(t) represents the discharge efficiency of the ESS in period t. The following formula can be used to αESS(t): where αc and αd are the charge and discharge efficiencies of the ESS, respectively DGs, however, fail to achieve the maximum prediction error of the net loa Because of this, we should take into account flexibility in the power constraints o when charging and discharging. This is because the secondary function of the meet these maximum prediction errors. Here are the revised constraints: where FESS.U and FESS.D are the reserved SOC of the ESS and ηF.ESS is the percenta maximum charge-discharge capacity per hour to the maximum capacity of the considering flexibility.
If the maximum prediction error of the net load power occurs, DGs are ca

Constraints of the ESS Considering Flexibility in Charging and Discharging
In this study, the main function of the ESS is to provide assistance to DGs and ensure that DGs have adequate flexibility. Thus, the ESS and its conventional constraints are similar in terms of charging and discharging power: where S ESS (t) represents the state of charge (SOC) of the ESS in period t; δ represents the selfdischarge efficiency of the ESS; E ESS represents the maximum capacity of the ESS; S ESS.min and S ESS.max represent the minimum and maximum SOC of the ESS, respectively; ∆t represents an hour; η ESS represents the percentage of the maximum charge-discharge capacity per hour to the maximum capacity of the ESS; α ESS (t) represents the charge-discharge efficiency of the ESS in period t. The following formula can be used to calculate α ESS (t): where α c and α d are the charge and discharge efficiencies of the ESS, respectively. DGs, however, fail to achieve the maximum prediction error of the net load power. Because of this, we should take into account flexibility in the power constraints of the ESS when charging and discharging. This is because the secondary function of the ESS is to meet these maximum prediction errors. Here are the revised constraints: where F ESS.U and F ESS.D are the reserved SOC of the ESS and η F.ESS is the percentage of the maximum charge-discharge capacity per hour to the maximum capacity of the ESS after considering flexibility. If the maximum prediction error of the net load power occurs, DGs are capable of meeting most of these prediction errors. There is only a small portion of these prediction errors that must be met by the ESS. As a result, the constraints of F ESS.U , F ESS.D , and η F.ESS are as follows: Since the ESS operation is periodic, the initial SOC should be equal to the final SOC in a day and expressed as follows: where T represents the total scheduling time.

Tie-Line Transmission Power Constraints in Consideration of Flexibility
In the MREPS, DGs and the ESS can sometimes fail. In this case, the tie-line must be flexible to facilitate the normal functioning of the MREPS. Therefore, the power limitation for tie-line transmission considering flexibility is expressed as follows: where P TL.min and P TL.max represent the minimum and maximum transmission power of the tie line, respectively, and F TL.U (t) represents the reserved transmission power of the tie-line that can be increased in period t.
Load shedding is prohibited in the operation of the MREPS. In the event that DGs and the ESS fail, it is necessary to reserve the transmission power on the tie-line, which may be increased in case of a breakdown. However, if the load power decreases, then we may be able to achieve power balance by reducing the output power of DGs or curtailing the use of renewable energy sources. As a consequence, it is not necessary to reserve the transmission power of the tie-line.

MOPSO
The multi-objective optimization problem in the presented model is solved via MOPSO in this study. Although MOPSO is one of the commonly used intelligent optimization algorithms [22,23], it is briefly introduced since it is the core of the algorithm in this paper. MOPSO updates the position and the velocity of each particle in every iteration and searches for the local and global best positions of particles. The velocity is modified and updated as follows: where v k i and x i are the velocity and position of the ith particle in the kth iteration; p k l is the local optimal position of the ith particle; p k g is the global optimal position of particles; ω is the inertia weight that influences the local and global exploitation abilities for MOPSO; c 1 and c 2 are the cognitive and social learning factors that maintain the movement of particles to the local and global optimal positions, respectively; r 1 and r 2 are two uniform random functions in the range [0,1].
Constraints of the velocity are expressed as follows: where v min and v max are the minimum and maximum velocities of particles, respectively. The position of particles is indicated as follows: MOPSO obtains the non-dominated solution set, unlike the single-objective PSO algorithm. Therefore, an external file is needed to store the set. The flowchart of MOPSO is shown in Figure 3. Processes 2022, 10, x FOR PEER REVIEW MOPSO obtains the non-dominated solution set, unlike the single-objectiv algorithm. Therefore, an external file is needed to store the set. The flowchart of M is shown in Figure 3. The SOC of the ESS, the output power of DGs, and the tie-line transmission are unknown in this study. Therefore, the SOC of the ESS and the output power in each period are considered decision variables. The tie-line power of transmis each period can be obtained using Equation (10). Due to the fact that the set of Pare dominated solutions is determined by these decision variables, each solution in corresponds to an operation state of the MREPS. Therefore, if we wish to determ optimal day-ahead scheduling strategy, then we need to first identify the op compromise solution from the set of non-dominated Pareto solutions.

Fuzzy Comprehensive Analysis Methodology
This study employs the fuzzy comprehensive evaluation method to determ most appropriate compromise solution. An effective and widely adopted meth evaluating hierarchical and integrated problems is a fuzzy comprehensive evaluati top two critical steps of this method are determining the weight vector of each eva objective appropriately and selecting the appropriate fuzzy membership function. 4 illustrates the fuzzy comprehensive evaluation method flowchart.  Assume that X is the Pareto non-dominated solution set, which is expre follows: The SOC of the ESS, the output power of DGs, and the tie-line transmission power are unknown in this study. Therefore, the SOC of the ESS and the output power of DGs in each period are considered decision variables. The tie-line power of transmission in each period can be obtained using Equation (10). Due to the fact that the set of Pareto non-dominated solutions is determined by these decision variables, each solution in the set corresponds to an operation state of the MREPS. Therefore, if we wish to determine the optimal day-ahead scheduling strategy, then we need to first identify the optimum compromise solution from the set of non-dominated Pareto solutions.

Fuzzy Comprehensive Analysis Methodology
This study employs the fuzzy comprehensive evaluation method to determine the most appropriate compromise solution. An effective and widely adopted method for evaluating hierarchical and integrated problems is a fuzzy comprehensive evaluation. The top two critical steps of this method are determining the weight vector of each evaluation objective appropriately and selecting the appropriate fuzzy membership function. Figure 4 illustrates the fuzzy comprehensive evaluation method flowchart.
MOPSO obtains the non-dominated solution set, unlike the single-objective PSO algorithm. Therefore, an external file is needed to store the set. The flowchart of MOPSO is shown in Figure 3. The SOC of the ESS, the output power of DGs, and the tie-line transmission power are unknown in this study. Therefore, the SOC of the ESS and the output power of DGs in each period are considered decision variables. The tie-line power of transmission in each period can be obtained using Equation (10). Due to the fact that the set of Pareto nondominated solutions is determined by these decision variables, each solution in the set corresponds to an operation state of the MREPS. Therefore, if we wish to determine the optimal day-ahead scheduling strategy, then we need to first identify the optimum compromise solution from the set of non-dominated Pareto solutions.

Fuzzy Comprehensive Analysis Methodology
This study employs the fuzzy comprehensive evaluation method to determine the most appropriate compromise solution. An effective and widely adopted method for evaluating hierarchical and integrated problems is a fuzzy comprehensive evaluation. The top two critical steps of this method are determining the weight vector of each evaluation objective appropriately and selecting the appropriate fuzzy membership function. Figure  4 illustrates the fuzzy comprehensive evaluation method flowchart.  Assume that X is the Pareto non-dominated solution set, which is expressed as follows: 11 1 1 Assume that X is the Pareto non-dominated solution set, which is expressed as follows: where m is the total number of objectives, n is the total number of the Pareto non-dominated solutions, and x ij is the value of the ith objective of the jth Pareto non-dominated solution.
To determine the most effective compromise solution from the Pareto non-dominated solution set, follow the steps outlined below: (1) For determining the membership degree of each objective in each Pareto non-dominated solution, a single-factor fuzzy evaluation is adopted. The fuzzy relation matrix can then be obtained as follows: where r ij is the membership degree of the ith objective in the jth Pareto non-dominated solution. r ij can be calculated as follows: where X i and X i are the minimum and maximum expectations of the decision-maker for the ith objective, respectively. (2) The analytic hierarchy process (AHP)-entropy weight method (EWM) can be employed to determine the comprehensive weight vector for each objective. Assume that the comprehensive weight vector is: where ω i is the comprehensive weight vector of the i th objective. A can then be calculated using the following equation: where ω i and ω i are the objective weights based on AHP and EWM, respectively. The discussion on AHP [26] and EWM [27] is excluded from this paper given that both have been extensively analyzed in the literature. (3) The comprehensive fuzzy evaluation vector B can be calculated as follows: where b j is the membership degree of the jth Pareto non-dominated solution.
When b j is close to 1, we can evaluate the jth Pareto non-dominated solution comprehensively. Therefore, from the set of Pareto non-dominated solutions, the optimal compromise solution corresponds to the maximum b j .

Case Study
Practical MREPS data are presented in this study within a simulation and analysis environment. The parameters of each power supply unit in the MREPS are listed in Table 1. The cost of each distributed generation is presented in Table 2. Table 3 presents the pollutant emission coefficients generated by each distributed generation. Table 4 shows the electricity purchase price of the MREPS from the main grid [15]. This study adopts a simulation of the load and output date of a typical day in the MREPS for WTs and PV arrays. To determine an optimal scheduling strategy, the forecast values of these data are used. Utilizing the realized values of these data, the optimal scheduling strategy is evaluated. Figures 5 and 6 present the predicted and realized curves for the load and output from WTs and PV arrays in the MREPS for a typical day, respectively.    As shown in Section 4, it is possible to determine the optimal scheduling strategy of a typical day after obtaining the non-dominated solution set. On the basis of the models, constraints, data, and decision variables of this study, the non-dominated solution set of objectives can therefore be obtained using MOPSO. Figure 7 and Figure 8 illustrate nondominated solutions based on the conventional and proposed strategies, respectively. From Figures 7 and 8, we can see that the three objectives cannot reach the optimal solution at the same time for the non-dominated solution set obtained by MOPSO. The distribution of the non-dominated solution set shows a narrow shape at both ends and is wide in the middle.    As shown in Section 4, it is possible to determine the optimal scheduling strategy of a typical day after obtaining the non-dominated solution set. On the basis of the models, constraints, data, and decision variables of this study, the non-dominated solution set of objectives can therefore be obtained using MOPSO. Figure 7 and Figure 8 illustrate nondominated solutions based on the conventional and proposed strategies, respectively. From Figures 7 and 8, we can see that the three objectives cannot reach the optimal solution at the same time for the non-dominated solution set obtained by MOPSO. The distribution of the non-dominated solution set shows a narrow shape at both ends and is wide in the middle. As shown in Section 4, it is possible to determine the optimal scheduling strategy of a typical day after obtaining the non-dominated solution set. On the basis of the models, constraints, data, and decision variables of this study, the non-dominated solution set of objectives can therefore be obtained using MOPSO. Figures 7 and 8 illustrate nondominated solutions based on the conventional and proposed strategies, respectively. From Figures 7 and 8, we can see that the three objectives cannot reach the optimal solution at the same time for the non-dominated solution set obtained by MOPSO. The distribution of the non-dominated solution set shows a narrow shape at both ends and is wide in the middle.   Furthermore, Figures 7 and 8 indicate that the non-dominated solution set of the proposed strategy is larger than that of the conventional strategy. According to this finding, after evaluating the flexibility constraints, the MREPS will sacrifice a portion of the economy, the utilization rate of renewable energy, and power fluctuations.
We must then determine the optimal compromise solution from the Pareto nondominated solution set to find the optimal day-ahead scheduling strategy for a typical day. The fuzzy comprehensive evaluation method can be used to obtain optimal compromise solutions from non-dominated Pareto solution sets, based on the conventional and proposed strategies, as discussed in Section 4. By the fuzzy comprehensive evaluation method, Point A and Point B are the solutions with the highest evaluation scores in the respective non-dominated solution sets of Figure 7 and Figure 8, respectively. Therefore, Points A and B represent optimal compromise solutions in Figure  7 and Figure 8, respectively.
Finally, based on the relevant decision variables, the optimal day-ahead scheduling strategy for a typical day according to every optimal compromise solution can be determined. These are illustrated in Figures 9 and 10. Furthermore, Figures 7 and 8 indicate that the non-dominated solution set of the proposed strategy is larger than that of the conventional strategy. According to this finding, after evaluating the flexibility constraints, the MREPS will sacrifice a portion of the economy, the utilization rate of renewable energy, and power fluctuations.
We must then determine the optimal compromise solution from the Pareto nondominated solution set to find the optimal day-ahead scheduling strategy for a typical day. The fuzzy comprehensive evaluation method can be used to obtain optimal compromise solutions from non-dominated Pareto solution sets, based on the conventional and proposed strategies, as discussed in Section 4. By the fuzzy comprehensive evaluation method, Point A and Point B are the solutions with the highest evaluation scores in the respective nondominated solution sets of Figures 7 and 8, respectively. Therefore, Points A and B represent optimal compromise solutions in Figures 7 and 8, respectively.
Finally, based on the relevant decision variables, the optimal day-ahead scheduling strategy for a typical day according to every optimal compromise solution can be determined. These are illustrated in Figures 9 and 10. Using the conventional and proposed strategies, Strategies A and B provide an optimal scheduling strategy for the typical day based on what has been denoted respectively, for ease of use. Figures 9 and 10 illustrate Strategy A's preference fo purchasing electricity from the main grid. It was found that the tie-line transmission power of Strategy A is close to its average. The consequence is a relatively low level o power fluctuations at the tie-line under Strategy A, whereas Strategy B prefers to generat electricity from DGs. It is more expensive to generate electricity through DGs than to purchase electricity from the main grid; therefore, Strategy B has a higher operation cost To better compare Strategies A and B, they were both applied in the context of a realistic scenario based on a typical day. Figure 5 shows such a situation. In Figure 11 and Figure 12, the realized operation results under Strategies A and B are presented respectively. Using the conventional and proposed strategies, Strategies A and B provide an optimal scheduling strategy for the typical day based on what has been denoted, respectively, for ease of use. Figures 9 and 10 illustrate Strategy A's preference for purchasing electricity from the main grid. It was found that the tie-line transmission power of Strategy A is close to its average. The consequence is a relatively low level of power fluctuations at the tie-line under Strategy A, whereas Strategy B prefers to generate electricity from DGs. It is more expensive to generate electricity through DGs than to purchase electricity from the main grid; therefore, Strategy B has a higher operation cost.
To better compare Strategies A and B, they were both applied in the context of a realistic scenario based on a typical day. Figure 5 shows such a situation. In Figures 11 and 12  The realized output of WTs and PV arrays is larger than the predicted va depicted in Figures 5 and 6. Consequently, the realized output power of DGs Figures 11 and 12 is smaller than the predicted value. However, DGs based on Stra generally have a minimum output power. DGs fail to reduce the output power case. Hence, curtailing renewable energy is the most economical way to maintain balance after the ESS absorbs the excess power to the extent possible. By compariso based on this strategy are able to reduce the output power because Strategy B flexibility constraints into account. Therefore, the MREPS can decrease the output of DGs to increase the utilization rate of renewable energy. In addition, Strategy B r the power purchased from the main grid during peak load periods from 18:00 to 2 comparison to Strategy A. This phenomenon significantly eases the demand for po the main grid during peak load periods.
The simulation results of Strategies A and B are listed in Table 5. It can be se Strategy A has a lower forecast operation cost, forecast curtailment rate of ren energy, and forecast power fluctuations on the tie-line. All three indicators are hig Strategy B. However, the IFR is 31.47% higher, the FSR is 45.83% lower, and the 13.36 kW/h higher for Strategy A compared to Strategy B. This indicates that Stra is less capable of coping with uncertainties during the real-time operation of the M compared to Strategy B. This is corroborated by the real-time operation results of t strategies. The realized output of WTs and PV arrays is larger than the predicted value, as depicted in Figures 5 and 6. Consequently, the realized output power of DGs in both Figures 11 and 12 is smaller than the predicted value. However, DGs based on Strategy A generally have a minimum output power. DGs fail to reduce the output power in this case. Hence, curtailing renewable energy is the most economical way to maintain power balance after the ESS absorbs the excess power to the extent possible. By comparison, DGs based on this strategy are able to reduce the output power because Strategy B takes flexibility constraints into account. Therefore, the MREPS can decrease the output power of DGs to increase the utilization rate of renewable energy. In addition, Strategy B reduces the power purchased from the main grid during peak load periods from 18:00 to 22:00 in comparison to Strategy A. This phenomenon significantly eases the demand for power on the main grid during peak load periods.
The simulation results of Strategies A and B are listed in Table 5. It can be seen that Strategy A has a lower forecast operation cost, forecast curtailment rate of renewable energy, and forecast power fluctuations on the tie-line. All three indicators are higher for Strategy B. However, the IFR is 31.47% higher, the FSR is 45.83% lower, and the AIF is 13.36 kW/h higher for Strategy A compared to Strategy B. This indicates that Strategy A is less capable of coping with uncertainties during the real-time operation of the MREPS compared to Strategy B. This is corroborated by the real-time operation results of the two strategies.

Parameters and Units Strategy A Strategy B
Insufficient flexibility rate of the day-ahead scheme (IFR) (%) [15] 46.21 14.74 Flexibility sufficiency rate (FSR) (%) [28] 4.17 50.00 Average insufficiency of flexibility (AIF) (kW·h) [28] 28 What is more, the MREPS is operating at a higher cost under Strategy B, and the tie-line power fluctuations are 9.52% higher than they are under Strategy A. However, the Processes 2022, 10, 1401 16 of 17 rate of utilization of renewable energy for the MREPS under Strategy B is 18.03% higher than that under Strategy A. In addition, the deviation rate of forecast and realized power fluctuations of power under Strategy B is 14.59% lower than that under Strategy A. These findings indicate that Strategy B can reduce the deviation between forecast and realized power fluctuations on the tie-line and increase the utilization of renewable energy at the expense of a small amount of economy and power fluctuations at the same time. Moreover, the realized values of the three objectives under Strategy A increased by 16.79%, 37.58%, and 0.47% compared with their forecast values. In comparison with their forecast values, the realized values of the three objectives under Strategy B increased by 12.27%, 12.59%, and 0.05%. Based on this finding, the results formulated by Strategy B are more consistent with the realized operation results.

Conclusions
In this study, a day-ahead scheduling strategy designed to account for flexibility constraints was presented. A multi-objective problem was solved using MOPSO, and a set of the non-dominated solutions was derived for the three objectives. Using the fuzzy comprehensive evaluation method, the optimal compromise solution of the non-dominated solution set was determined. In addition, the optimal strategy proposed in this study was the day-ahead scheduling strategy corresponding to the optimal compromise solution. During peak load periods, the simulation results showed that the proposed strategy was effective in relieving the main grid's pressure on power supply. Moreover, compared with those formulated using the conventional strategy, the results obtained from day-ahead scheduling formulated using the proposed strategy were closer to the results obtained from the MREPS. Although the economy and power fluctuations on the tie-line were slightly higher under the revised strategy, renewable energy usage was significantly higher, and the differences between forecast and realized power fluctuations on the tie-line were relatively small. This finding showed that the revised strategy has the potential to significantly improve the flexibility and reliability of the MREPS's operation at the cost of a small amount of economy and fluctuations in power supply.