A Two-Stage Robust Optimization Method Based on the Expected Scenario for Islanded Microgrid Energy Management

One of the main challenges in microgrid system energy management is dealing with uncertainties such as the power output from renewable energy sources. The classic two-stage robust optimization (C-TSRO) method was proposed to cope with these uncertainties. However, this method is oriented to the worst-case scenario and is therefore somewhat conservative. In this study, focusing on the energy management of a typical islanded microgrid and considering uncertainties such as the power output of renewable energy sources and the power demand of loads, an expected-scenario-oriented two-stage robust optimization (E-TSRO) method is proposed to alleviate the conservative tendency of the C-TSRO method because the E-TSRO method chooses to optimize the system cost according to the expected scenario instead, while ensuring the feasibility of the ﬁrst-stage variables for all possible scenarios, including the worst case. According to the structural characteristics of the proposed model based on the E-TSRO method, a column-and-constraint generation (C & CG) algorithm is utilized to solve the proposed model. Finally, the eﬀectiveness of the E-TSRO model and the solution algorithm are analysed and validated through a series of experiments, thus obtaining some important conclusions, i.e., the economic eﬃciency of system operation can be improved at about 6.7% in comparison with the C-TSRO results.


Introduction
Microgrids are small power generation and distribution systems composed of distributed power supplies, energy storage systems (ESSs), user loads, energy conversion devices, and protective devices. ey are important for the effective utilization of renewable energy resources and to improve the reliability of the energy supply under various loads. Microgrids can be classified into two types, namely, islanded and grid-connected, based on whether or not they connect to the power grid. Islanded microgrids are primarily used to supply power to remote regions that have yet to be reached by the power grid (e.g., islands). Grid-connected microgrids are connected to the power grid, through which energy can be traded based on the difference between the peak and off-peak electricity prices. In the event of a large power grid malfunction, a grid-connected microgrid can switch to the islanding mode.
Energy management is an important part of the scheduling and operational management of a microgrid, wherein scheduling decisions are made based on the behaviour of its components (e.g., the power generators, ESS, and response loads) to ensure a stable and efficient operation. Various uncertainties such as the power output from renewable energy sources (e.g., wind and solar energy) and loads pose significant challenges in microgrid energy management, and their effective treatment has become an important research topic in this domain. Accordingly, researchers have proposed two main types of uncertainty modelling methods for microgrid energy management, namely, stochastic programming and robust optimization [1].
In the stochastic programming method, the probability distribution of an uncertain parameter is assumed to be known in advance. is probability distribution is then used to transform an uncertainty problem to a deterministic one.
Researchers [2][3][4][5][6] have developed stochastic programming models for the energy management of grid-connected microgrids to minimize the overall cost under uncertainties in single or multiple factors (e.g., electricity transaction price, power output from renewable energy sources, and loads). ese models differ mainly in the composition of the target cost, the uncertain factors considered, and the scenario generation method. Other researchers [7][8][9][10] have developed models for grid-connected microgrids with optimal energy management objectives that are not limited to minimization of the overall cost. e energy management models developed by Tabar et al. [8] and Gazijahani et al. [10] use multiple optimization objectives. e model built by Tabar et al. [8] has two optimization objectives, namely, to minimize the operational cost and the level of air pollution, whereas those in the model developed by Gazijahani et al. [10] aim to minimize the operational cost and to maximize the power supply reliability. e optimization objectives of the stochastic programming energy management model developed by Hu et al. [9] are divided into two stages: reduction of the investment cost of the microgrid in the first stage and minimization of the operational and management cost of the microgrid in the second stage. In addition to these studies on grid-connected microgrids, studies have also been conducted on islanded microgrids [11][12][13]. Cau et al. [11] investigated an islanded microgrid with a hydrogen production unit and hydrogen fuel cells as well as renewable (wind and solar) energy sources and energy storage batteries; the optimization objective of the proposed energy management model was to maximize the utilization of the renewable energy sources. Rezaei and Kalantar [12] considered the frequency control problem of islanded microgrids and developed a microgrid energy management model that minimized the deviation from the day-ahead frequency. Farzin et al. [13] optimized the expected operational and management cost of a microgrid during an islanding event. In addition to scenario-probability-based stochastic programming, researchers [14] have used chance constraintbased stochastic programming to handle uncertain factors in microgrid energy management. e main difficulty in applying stochastic programming-based microgrid energy management is the accurate determination of the multiple scenarios used to describe uncertainties, particularly the probability of occurrence of each scenario. Moreover, the complexity of the model solution scheme increases with the number of scenarios used to represent the uncertainties [15]. e robust optimization method does not require a priori knowledge of the probability distribution for an uncertain parameter. Instead, the fluctuation range of an uncertain parameter is represented by a specified interval, and a feasible solution can be guaranteed with a high probability as long as the realized value of the uncertain parameter lies in the given interval [16]. According to the stage division of decision variables, available studies on microgrid energy management based on robust optimization can be divided into two types, namely, those based on single-stage robust optimization [17][18][19][20][21] and those based on multistage robust optimization. In the former, all the decision variables for microgrid energy management are determined at the same time. For example, Kuznetsova et al. [20] and Zhang et al. [21] investigated microgrids that only contain electrical power whereas Wang et al. [17] and Luo et al. [18] studied microgrids that consist of CHP (combined heat and power) and CCHP (combined cooling, heating, and power), respectively. Microgrid energy management methods based on single-stage robust optimization offer the advantage of relatively simple models that can be solved easily, but these models require decision-makers to make decisions in parallel before knowing the realized values of the uncertain parameters. Typically, there are equality constraints (i.e., energy supply and demand should be balanced at all times) in microgrid energy management models, which leads to conservative decisions and is inconsistent with the fact that decisions are often made sequentially in stages in engineering practice. In contrast, in multistage optimization such as two-stage robust optimization [22][23][24][25], some decisions (e.g., the on/off status of the generators, the charge/discharge status of the ESS, the transaction status, or preagreed quantity of power transacted with the power grid) are made before knowing the realized values of the uncertain parameters whereas the remaining variables are determined after the values of the uncertain parameters have been determined. erefore, this approach can be used to deal with equality constraints in microgrid energy management models. To cope with the min-max-min optimization objectives embedded in two-stage robust optimization models, researchers typically adopt the columnand-constraint generation (C & CG) algorithm [26] to solve these models [23][24][25] but other algorithms have also been proposed. For example, the Lyapunov optimization algorithm (Hu et al. [27]) has been utilized to yield approximate solutions. In addition, Xiang et al. [22] established a scenario-robust-combined two-stage robust optimization model for energy management and proposed an orthogonal array-based search algorithm to solve the model. For an isolated microgrid on Nanji Island in China that contains a seawater pumped storage system, Tian et al. [28] established a multitime-scale islanded microgrid energy management model based on robust optimization and developed a Bender's decomposition and cutting plane algorithm to solve the proposed model. Zhang et al. [29] established a two-stage robust optimization model for multimicrogrid energy management, which was solved using the C & CG algorithm. Hussain et al. [30] studied grid resilience by investigating possible outages to ensure optimal operation and scheduling of critical loads. In the aforementioned studies, the decision variables were determined using the classic two-stage robust optimization (C-TSRO) method that is oriented to the worst-case scenario. However, in practice, the probability of the worst-case scenario occurring is very low, making the solutions obtained using the C-TSRO method relatively conservative, especially as continuous improvement by an interval prediction algorithm often decreases the difference between the actual and predicted expected values of the uncertain parameters. As a result, the solution obtained based on the worst-case scenario is always worse than that based on the expected scenario considering the optimization objectives (i.e., economic efficiency) [31]. 2 Discrete Dynamics in Nature and Society Targeting on the drawback in adopting the C-TSRO method, an expected-scenario-oriented two-stage robust optimization (E-TSRO) model is proposed for microgrid energy management of typical islanded microgrid systems composed of renewable energy equipment (wind and solar), ESSs, diesel generators, and various types of loads. e model accounts for uncertainties in factors such as the power output from renewable energy sources and the load demand.
ere are two stages in the E-TSRO model for determining the decision variables of microgrid energy management, namely, prescheduling (day-ahead) and rescheduling (intraday). Compared to the C-TSRO method oriented to the worst-case scenario, the E-TSRO method primarily focuses on the optimization function (i.e., maximizing the economic efficiency of the system) based on the expected scenario, while ensuring the feasibility of solutions under any realized scenario. e E-TSRO method ensures that there are always feasible rescheduled solutions for these second-stage decision variables under any realized scenario based on the predetermined values of the first-stage decision variables, resulting in less conservative decisions. us, safe and stable microgrid operation can be achieved in uncertain environments while maintaining a relatively high economic efficiency. Taking into consideration the nonlinear formulation of the proposed model and the complexity of solving this type of model, the model is equivalently transformed and a C & CG algorithm is developed to solve the transformed model. Experiments are used to validate the advantages of the E-TSRO model and the effectiveness of the developed algorithm.
e main contributions of this paper are summarized as follows.
(1) Instead of orienting to the worst case, an expectedscenario-oriented two-stage robust optimization (E-TSRO) method is proposed to alleviate the conservative tendency of the classic two-stage robust optimization (C-TSRO) method. (2) According to the structural characteristics of the proposed model based on the E-TSRO method, a column-and-constraint generation (C & CG) algorithm is utilized to solve the proposed model. is paper is organized as follows. In Section 1, the research background and literature review are presented. Section 2 provides a basic model for islanded microgrid energy management. In order to solve the proposed model, the E-TSRO framework and its transformation along with the C & CG solution algorithm are developed in Section 3. Section 4 depicts the experimental validation process and analytical results. Finally, summary and future work are concluded in Section 5.

Basic Model for Islanded Microgrid
Energy Management e islanded microgrid investigated in this study primarily consists of wind and photovoltaic (PV) power generators, loads (including critical and curtailable loads), one ESS, and distributed generators (DGs).
e notations used in this study and their corresponding abbreviations are listed as follows.

Notation
Description. Tables 1-3 depict system parameters, predicted parameters, and decision variables, respectively, that are used in the model. e main constraints and optimization objectives of the basic model for islanded microgrid energy management are described below.

DG Constraints.
In this study, DGs primarily refer to controllable generators and are represented by diesel generators. For a DG, constraints are placed on the maximum ramp power, the maximum and minimum operating power, and the minimum duration before startup and shutdown.
(1) Generation power and ramp power constraints: Constraint equations (1) and (2) are the ramp power and generation power constraints on the generator, respectively. (2) Duration constraints for startup and shutdown: Constraint equation (3) Discrete Dynamics in Nature and Society Constraint equations (5) and (6) are methods for calculating the startup and shutdown costs, respectively, for DG i over the time interval k, which are identical with [15].

Constraints on the ESS. An ESS plays a vital role in an
islanded microgrid and can effectively mitigate asynchronicity between power generation from renewable energy sources and power consumption by users. e main constraints on an ESS are given below: A positive P bess (k) corresponds to a charged ESS whereas a negative P bess (k) corresponds to ESS discharge. Constraint equations (8) and (9) show the charge and discharge power limit and capacity limit, respectively. Equation (10) describes the relationship between the ESS capacity and the ESS charge/discharge power. Since the primary focus of this study is the performance of two-stage robust optimization for microgrid energy management based on the expected scenario, to simplify the subsequent dual transformation process of the model, this study does not distinguish the charge power coefficient from the discharge power coefficient, and both coefficients are denoted as η bess in the remainder of this study.

Constraint on Curtailable
Load. Two types of loads are mainly considered in the proposed model, namely, critical and curtailable loads. Critical loads are closely related to the basic daily activities of users. As the normal operation of critical loads directly affects the daily life of users, the power demand of these loads (e.g., indoor illumination and  Startup and shutdown costs for DG i (€), respectively Table 2: Predicted parameters. P solar (k), P wind (k), P cri (k), and P cur (k) Predicted power of PV generators, wind generators, critical loads, and curtailable loads at time interval k (kW), respectively P solar (k), P wind (k), P cri (k), and P cur (k) Expected value of the predicted power for PV generators, wind generators, critical loads, and curtailable loads during time interval k (kW), respectively P solar (k), P wind (k), P cri (k), and P cur (k) Deviation from the expected value of predicted power for PV generators, wind generators, critical loads, and curtailable loads during time interval k, respectively Capacity of the ESS at time interval k (kW · h) P bess (k) Charging/discharging power of the ESS at time interval k (kW) β cur (k) Curtailed power percentage of the curtailable loads for time interval k (%) infrastructure) must be met unconditionally. In contrast, curtailable loads (e.g., air conditioners and ventilation systems) do not have to operate at rated powers, that is, the operating power can be adjusted within a certain range, which is shown below: Moreover, considering that load curtailment affects user comfort, the penalty or compensation cost associated with load curtailment will be considered in calculating the total operational cost.

Microgrid Energy Balance Constraint.
Stable and reliable operation at all times is the precondition for optimizing the energy scheduling of an islanded microgrid; therefore, a balance of energy supply and demand must be maintained as follows:

Optimization Objective of the Basic Model.
e optimization objective of microgrid energy management is to minimize the total operational cost of the system, which is expressed as e three terms within the brackets are the startup cost, the shutdown cost, and the maintenance and management costs, respectively. e term C total is the sum of the energy storage cost, the fuel cost of the generators, and the penalty cost of load curtailment. C total is expressed as where C bess (k) is the ESS maintenance cost for the time , and C cur (k) is the cost associated with the impact of load curtailment on user comfort, i.e., the penalty or compensation cost of load curtailment (C cur (k) � α cur (k)β cur (k)P cur (k)Δt, where α cur,i (k) reflects the preference of various users).

Description of Uncertain Parameters.
e primary focus of this study is the uncertainty in two types of factors that affect islanded microgrid energy management, namely, the power output from renewable (wind and solar) energy sources and the power of loads (critical and curtailable loads). It is assumed that the uncertainties in these factors can be represented by intervals. Here, the predicted output power of solar energy, P solar (k), for the time interval k is used as an example.
where P solar (k) and P solar (k) are the predicted expected output power and the deviation from the expected value of solar energy over time interval k, respectively. Spatial and temporal fluctuations in the uncertainties are incorporated in the model using two adjustable budget-of-uncertainty (BoU) parameters, Γ S and Γ T [33], respectively, which reflect the risk preference of the decision-maker.
As all the constraints in the proposed model are linear, the worst-case scenario will inevitably happen at one of the vertices of the polyhedral uncertain set [33]. e uncertain parameters considered in this study are given below: Considering that the uncertain parameters incorporated in this paper include the power of renewable energy (i.e., wind and solar) and the power of critical and curtailable loads, the probability that actual outputs of these uncertain parameters reach their upper or lower bound at the same time is very small. erefore, a spatial constraint on the BoU Discrete Dynamics in Nature and Society parameter of Γ S is imposed denoting by equation (16). Similarly, equation (17) represents the temporal constraint on the BoU of Γ T .

E-TSRO Model.
e aforementioned uncertain parameters and their expressions were used to classify all the decision variables into decision variables for the day-ahead (prescheduling) stage and the intraday (rescheduling) stage. Considering limitations such as the startup and shutdown costs of DGs, the required durations before startup and shutdown, and the ramp power, only the on/off status variable of the DGs is treated as a day-ahead stage decision variable, i.e., the value of the on/off status variable must be determined a day ahead based on the predicted ranges of the uncertain parameters. e remaining decision contents are intraday or rescheduling decision variables.
e E-TSRO model is expressed below [34]: ∀u ∈ U, ∃y: x ∈ X, y 0 ≥ 0, y ≥ 0, where x and y are first-and second-stage decision variables, respectively, u is an uncertain parameter, and y 0 is the decision corresponding to y when the value of the uncertain parameter u is set to its predicted expected value u 0 . Equation (19a) shows the operational cost of the system for minimization under the expected scenario. Equation (19b) is the constraint on the system under the expected scenario. Equation (19c) expresses the robust feasibility of the prescheduling solution x, i.e., for a given prescheduling strategy x, when the uncertain parameter u is set to a value within the given range, there is a feasible solution for the second-stage decision variable y. Equations (19a)-(19d) show that, in the E-TSRO model, the optimization of the cost under the worst-case scenario in the C-STRO model is replaced by the optimization of the cost d T y 0 under the expected scenario. e E-TSRO model also requires the prescheduling strategy to be feasible under any scenario. e change in the focus of the model from the worst-case scenario to the expected scenario results in a less conservative decision, while ensuring that a feasible decision can be obtained when knowledge of the uncertain parameter values becomes available. us, an optimal decision is ensured under the expected scenario and feasible under the worst-case scenario.
e correspondence between equations (19a)-(19d) and the microgrid energy management model presented in Section 2 is as follows: equation (19a) corresponds to equation (13); equation (19b) corresponds to equations (1)- (12), and the only difference being that the predicted range of the uncertain parameter is replaced by the predicted expected value; corresponds to replacing the predicted value of the uncertain parameter by its actual value in equations (1)- (12); for equation (19d), x represents δ DG (k), and y represents the remaining decision variables. However, equation (19c) makes it impossible to directly solve the model given by equations (19a)-(19d). Instead, a set of prescheduling solutions must be found, which satisfies robust feasibility, and a corresponding test model must be established. A method for establishing a test model based on equations (19a)-(19d) and a solution algorithm are presented below.

Robust Feasibility Test Model for Prescheduling Solutions.
e method reported by [33,35] for establishing a robust feasibility test model is used here to determine a set of prescheduling solutions that satisfy robust feasibility. First, the robust feasible region of the prescheduling variable x denoted by equation (19c) is rewritten as follows: where Y(x, u) � y|By ≤ b − Cu − Ax is a feasible solution to the rescheduling decision variable y based on the given prescheduling variable x and the uncertain parameter u. For the given u ∈ U, two positive slack variables, s + and s − , are introduced to test whether Y(x, u) is null. ese variables can be interpreted as emergency measures to prevent system failure. Accordingly, a linear programming problem is defined as follows: e objective function in equation (21) is to minimize the sum of the slack variables s + and s − based on the given prescheduling decision x * and knowledge of the value of the uncertain parameter u. r(x * , u) � 0 means Y(x * , u) ≠ ∅. Based on the model established in Section 2.2, the constraints corresponding to those in equation (21) are as follows: Based on the definition above, if r(x * , u) � 0, then Y(x * , u) ≠ ∅; otherwise, if r(x * , u) > 0, then Y(x * , u) � ∅, which means a feasible rescheduling result cannot be obtained. To analyse the robustness of x * , the result of equation 6 Discrete Dynamics in Nature and Society (21) must be tested under the worst-case scenario based on the uncertainty in u. us, the following equation is defined based on equation (21): e uncertainty set U is determined by equations (15)- (18). is max-min form is used to determine whether a feasible rescheduling decision can be obtained under the worst-case scenario based on the given prescheduling decision x * , which can be addressed with KKT condition or dual theory [26]. As for equation (21), if R(x * ) � 0, then Y(x * , u) ≠ ∅, and vice versa. As there is no integer variable in y, the max-min decision problem of equation (23) is dually transformed to the max bilinear programming problem shown below: where w is the dual of the variable in equation (24) and (23) shows that the prescheduling strategy is robust feasible only for R(x * ) � 0 [33]. Otherwise, new constraints must be continually added to the main problem (MP) given by equations (19a) and (19b) until R(x * ) � 0. Here, we use the C & CG algorithm to solve the model according to the main principle as follows. Based on the aforementioned two-stage classification, a set of prescheduling solutions is obtained by solving the MP given by equations (19a) and (19b); then, the robust feasibility test model given by equation (24) is solved using the obtained solutions. If R(x * ) > 0, it means that no feasible solution can be obtained during the rescheduling stage, i.e., the power balance of the system cannot be achieved, and a constraint that corresponds to the scenario determined by solving equation (24) should be added to the MP. e MP is then solved again. Cyclic interactive iterations are performed until R(x * ) � 0, which means a robust feasible prescheduling solution is found. e steps of the C & CG algorithm are outlined below, and the corresponding diagram is also illustrated as in Figure 1.

Algorithm for Solving the E-TSRO Model. Equation
Step 1. Set the number of iterations k � 0, x 0 � 0, and R 0 � 0.
Step 2. Solve the following MP: and denote the solution as x * k+1 .
Step 3. Solve the subproblem shown in equation (24) with the obtained solution x * k+1 . Denote the value of the second-stage variable y by y k+1 , the determined value of u by u * k+1 , and the obtained value of the optimization objective by R k+1 . If R k+1 � 0, take x * k+1 as a robust feasible optimal prescheduling solution and proceed to Step 4; otherwise, add the following constraint to the MP: Let k � k + 1 and return to Step 2 Step 4. Solve the rescheduling model as follows: where u r is the actual value of u and x * k+1 is the robust feasible optimal prescheduling solution obtained in Step 3.

Experimental Environment and Parameter Settings.
Simulations were performed on an islanded microgrid to analyse and validate the established E-TSRO model and the solution algorithm. e islanded microgrid contains three diesel generators, a battery ESS, PV and wind power generators, and critical and curtailable loads. e development, programming, and solution of the model, as well as the implementation of the rescheduling stage were built using IBM ILOG CPLEX V12.2. e model was used to perform day-ahead scheduling. e duration of energy scheduling was set to 24 hours that starts at 0:00 and ends at 24:00, and the unit time interval was set to 1 hour. e rated power of the wind and PV generators was set to 200 kW in order to make the total energy consumption and the energy supply of the system approximately balance within the setting 24 hours. Tables 4 and 5 summarize the relevant parameters of the DGs and ESS, respectively, and the parameter settings of diesel generators are mainly from literature [36].
Data for the wind, PV devices, and loads used in this study were obtained by integrating relevant data published by the Belgian Power Transmission Database in May 2018, as shown in Figure 2. e solid lines represent the predicted expected values of the uncertain parameters (wind power, solar power, and loads). For each time interval, the ratio of the critical loads to the curtailable loads in the system is 7 : 3. e areas enclosed by the dotted lines are the fluctuation ranges of the uncertainties. In addition, the prediction error percentage was set to increase linearly with time (from 3.5% at the 1st h to 15% at the 24th h at a rate of 0.5% per h), in order to reflect the trend that the prediction error normally increases over time.
Load curtailment can alleviate the energy supply burden when the energy supply is inadequate and is an indispensable adjustment measure for an islanded microgrid system. However, load curtailment may affect user comfort and is therefore generally not encouraged. As a result, the penalty for load curtailment should be set based on the actual situation. In this study, the proportion of curtailable Discrete Dynamics in Nature and Society loads at a single time point was set to 0.2. e power generation cost per unit of the DGs considered in the experiment is approximately from 0.4 €/kW to 0.65 €/kW. Hence, the load curtailment cost was set to 1 €/kW, which is approximately twice the power generation cost per unit of the generators.

Impact of the BoUs on the System Behaviour.
Equations (16) and (17) show that the prescheduling result is affected by the two BoUs used to describe an uncertain parameter, that is, Γ T , which is related to the number of periods, and Γ S , which is related to the number of uncertain parameters. Accordingly, the model was solved where Γ T is fixed as 24 and Γ S is set from 0 to 3, respectively, where 0 means all uncertain parameters are deterministic and 3 denotes that the output of solar energy, wind energy, and user load are simultaneously uncertain. e results of on/off status of the three DGs are summarized in Table 6. Similarly, Table 7 summarizes the status for the three DGs with different settings of Γ T and a fixed Γ S setting of 3. Tables 6 and 7 show that, as Γ S or Γ T increases, there is an overall increase in the conservativeness of the prescheduling solution obtained from the model, as reflected in the increased number of operational DGs. In addition, changing the BoU settings does not necessarily result in changes to the    (Table 6) and Γ T from 18 to 24 (Table 7) does not change the status of DGs. at is, the prescheduling results are clearly not overly sensitive to the changes in the BoUs.
To further analyse the impact of the BoUs on the rescheduling behaviour of the system, the following three special scenarios were defined according to the actual values assigned to the uncertain parameters at the second stage. (1) Expected scenario (Scenario A): each uncertain parameter is set to its predicted expected value. (2) Inadequate supply scenario (Scenario B): for each time interval, the wind and solar power are set to values to the lower bound of the predicted range, whereas the load demand is set to a value with the upper bound, denoting that a severely inadequate energy supply exists in this scenario. (3) Adequate supply scenario (Scenario C): for each time interval, the wind and solar power are set to values with the upper bound of the maximum predicted range, whereas the load demand is set with the lower bound. Tables 8 and 9 summarize the total cost of the system, the sum of the curtailed power, and the sum of the curtailed wind and PV power for the three typical scenarios under various BoU settings. e following conclusions can be drawn from the results in Tables 8 and 9. (1) In Scenario A, the total cost of the system increases slightly (not significantly) with the BoUs.
is finding is consistent with the impact of the BoUs on the system prescheduling result shown in Tables 6 and 7. (2) In Scenario B or C, at relatively low BoU values, no feasible rescheduling decision can be obtained for the system because as per the system requirements, the load curtailment ratio over each time interval cannot exceed 0.2 and wind and solar power curtailment is prohibited. us, for the given prescheduling decision, no rescheduling solution can be found that satisfies the power balance constraint on the system during the rescheduling stage, and this energy balance constraint can only be satisfied by increasing curtailment of the loads or wind and solar power. Analysing the relationship of the BoU settings to the total cost of the system, the total power of the curtailed loads, and the sum of the curtailed wind and solar power shows that increasing the

Period
Strategies in predispatch stage    BoUs moderately increases the operational cost of the system and reduces load or wind and solar power curtailment in extreme cases. us, the system stability is moderately improved at an economic cost. (3) For the BoU settings that a solution can be obtained, a comparison of the three scenarios shows that the total system cost is highest under Scenario B (i.e., an extremely inadequate energy supply) and lowest under Scenario C (i.e., adequate energy supply). is result can be attributed to the large number of loads that must be curtailed in Scenario B, which results in a relatively high penalty cost, whereas the relatively low operational power of the generators in Scenario C results in a significant reduction of the fuel cost.

Analysis of the System Behaviour for Given BoUs.
To determine the behaviour of typical devices (including the DGs and ESS) and users under given BoU settings, the BoUs were set to Γ T � 24 and Γ T � 3 using the aforementioned three scenarios as rescheduling scenarios, and the behaviour of the DGs, ESS, and curtailed loads for each scenario was recorded. Figure 3(a) shows the power of the DGs and the proportion of curtailed loads for each time interval under the three typical scenarios. e prescheduling decision results in a consistent on/off status, but there is a notable difference in the total output power of the DGs among the three scenarios. As the energy supply in Scenario B is the most inadequate of the three scenarios by comparing the power generated from the renewable energy sources with the actual load demand, the generators maintain relatively high power during their on status periods. In contrast, the operational generators in Scenario C output relatively low power for their working periods. User loads need to be curtailed only in Scenario B because of the relative large inadequacy in the energy supply of the system. Figure 3(b) shows the ESS behaviour under the three scenarios. In Scenario B, to effectively respond to the inadequacy in the energy supply resulting from the shutdown of the generators in period 15, the ESS starts charging in cycle 0 until it is fully charged. e ESS starts discharging in cycle 15 to compensate for the energy supply deficit caused by the shutdown of the generators. In contrast, in Scenario C, the energy supply is relatively adequate. us, the ESS starts discharging in cycle 0 and shifts to discharge mode in cycle 15 to effectively utilize the remaining generation power after the wind power generators have started to supply power to the loads.

Comparison with the C-TSRO Model.
Compared with the C-TSRO model, the most distinctive feature of the E-TSRO model is that it is oriented to the expected scenario instead of the worst-case scenario. erefore, experiments were conducted to compare the economic efficiencies of the two models.
First, prescheduling decisions were obtained by solving both models with the setting of Γ S � 3 and Γ T � 0, 6, 12, 18, and 24. en, equation (27) was solved by assigning values to the uncertain parameters to obtain rescheduling decisions under the different scenarios. For example, under Scenario B and Γ T � 12, P win d (k) was set to the lower bound of its predicted range for only 12 periods and to its expected value for the remaining periods. Table 10 is a comparison of the prescheduling and overall costs obtained by solving the two models for various Γ T values under the three rescheduling scenarios.  10 Discrete Dynamics in Nature and Society e following interesting phenomena can be observed from Table 10. (1) Both models yield rescheduling solutions for Γ T � 0 and 6 because the values assigned to the uncertain parameters at the rescheduling stage were based on the actual value of Γ T . In contrast, the results in Table 8 were obtained using a fixed Γ T � 24, resulting in no solutions under Scenarios B and C. (2) An analysis of the prescheduling costs in columns 2-4 of Table 10 shows that the C-TSRO model results in the generators being switched on/off a relatively large number of times and remaining operational over a relatively large number of periods. erefore, the prescheduling cost from the C-TSRO model is higher (between 6.5-8.9%) than that from the E-TSRO model.     shows that, under Scenarios A and C, the economic efficiency of the E-TSRO model is notably higher (between 4.5-8.3%) than that of the C-TSRO model. However, under Scenario B (i.e., an extremely inadequate supply), the E-TSRO model does not offer any advantages over the C-TSRO model (from −3.4% to −4.6%). is result is obtained because the C-TSRO model is oriented to the worstcase scenario, as represented by Scenario B. As a result, the generators are operational over more periods during the prescheduling stage compared with the solution of the E-TSRO model, preventing the curtailment of relevant loads at the rescheduling stage and the associated penalty effectively. e economic efficiencies of the two models were further compared by using the aforementioned BoU settings to obtain prescheduling decisions. irty scenarios were randomly generated at the rescheduling stage, and the rescheduling models were solved under these 30 scenarios to calculate the total system cost under each scenario. Table 11 summarizes the maximum, minimum, and average costs for the 30 random scenarios, showing that the E-TSRO model is moderately more effective than the C-TSRO model. In addition, comprehensive statistical results show that the average total cost is reduced by about 6.7% based on different setting of Γ T . e experiment also shows that the extreme case wherein load curtailment is required at the rescheduling stage does not occur under any of the 30 scenarios randomly generated at each Γ T value. is result further demonstrates that there is an infinitesimal probability of an extreme scenario occurring and underscores the necessity of responding to the impact of uncertain factors from an expected scenario perspective.

Conclusions
In this study, the uncertainties in factors such as power output from new (wind and solar) energy sources and user loads are addressed by dividing microgrid energy management decisions into day-ahead (prescheduling) stage decisions and real-time (rescheduling) stage decisions. ese decisions are based on the operational characteristics and constraints of typical devices (e.g., generators) in a microgrid. is approach is moderately effective in reducing the impact of uncertain factors. An E-TSRO model for energy management is developed based on a typical islanded microgrid, using the aforementioned two decision-making stages and modelling the aforementioned uncertainties by ranges. Unlike the C-TSRO model, the E-TSRO model is oriented towards the expected scenario, while ensuring the feasibility of these prescheduling decisions under any realized scenario. As a result, the E-TSRO model can effectively improve the economic efficiency of the C-TSRO model at about 6.7%. A C&CG algorithm is developed to solve the E-TSRO model considering the model characteristics. e effectiveness of the E-TSRO model is validated and compared to that of the C-TSRO model through numerous experiments.
e main conclusions are summarized as follows. (1) Although the BoU settings have some effect on the prescheduling decisions, this effect is not strictly linear, i.e., the prescheduling decision is not affected for some ranges of the BoU values. (2) Increasing the BoUs moderately increases the operational cost of the system. However, this increase also reduces load or wind and solar power curtailment in some extreme cases, thereby moderately improving the system stability. (3) In most of the investigated scenarios, the E-TSRO model is notably superior (i.e., from 4.5% to 8.9%) to the C-TSRO model in terms of the economic efficiency of system operation, with the exception of the extreme scenario of a severely inadequate energy supply.
In this study, the charge efficiency coefficient is not distinguished from the discharge efficiency coefficient to simplify the dual transformation of the proposed model. In subsequent studies, a variable that characterizes the ESS charge/discharge status should be introduced to distinguish between these two coefficients. In addition, the E-TSRO model can also be analysed and validated by application to the energy management of grid-connected microgrids.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.