Real-time low-carbon scheduling for the wind–thermal–hydro-storage resilient power system using linear stochastic robust optimization

Abstract

With the large-scale wind power integration, power systems have to address not only the conventional power demand fluctuations but also the wind uncertainty. To improve the economical effectiveness, resilience, and environmental protection of power systems in the source-load uncertainty, a real-time low-carbon scheduling for the wind–thermal–hydro-storage integrated system is proposed. The power imbalance caused by the uncertainty is neutralized by the synergetic linear decision of multiple resources. To address the source-load uncertainty, a stochastic robust optimization is introduced, which establishes the system constraints by robust optimization for the resilience operation, while optimizing the expected operation cost in the empirical uncertainty distribution for economic efficiency. Moreover, a multi-point estimation is applied to formulate the expected operation cost precisely and quickly. By using the dual theory, the proposed real-time power scheduling is derived as a mixed integer bilinear constrained programming. A multi-step sequential convexified solution is developed to solve the complex scheduling problem, which linearizes the bilinear constraints with alternate optimization and relaxes the state variables of energy storages with an “estimation–correction” strategy. Finally, case studies show the superiority of the proposed scheduling method and convexified solution.

1 Introduction

Large-scale wind power integration is of great significance to energy conservation and emission reduction in power systems (Zhu et al., 2021; He et al., 2022). However, due to the fluctuating nature of wind power, the resilience operation of power systems faces the uncertainty challenge of both wind power and power demand. Real-time scheduling, as a part of multiple time-scale power scheduling, plays an important role in restoring the power balance and improving the system resilience (Surender et al., 2015; Zhou et al., 2018). Therefore, how to address the source-load uncertainty with the co-regulation of wind farms, thermal units, hydro units, energy storages, and other regulation resources in real-time scheduling to ensure the system security should be paid more attention.

At present, methods for the uncertainty problems mainly contain the stochastic optimization and robust optimization. The stochastic optimization (Liu et al., 2022; Li et al., 2023) formulates the source-load uncertainty with discrete scenarios. The method samples multiple scenarios according to the probability distribution of the source-load uncertainty and optimizes the expected operation cost formulated with these sampled scenarios (Li and Xu, 2019). However, it is difficult for these sampled scenarios to cover the whole uncertainty distribution, and hence, the system security cannot be guaranteed in some extreme scenarios. Moreover, the huge sample size and complex calculation also make it difficult for the stochastic optimization to meet the timeliness in real-time scheduling. By contrast, robust optimization methods are widely applied since they can cover the extreme scenarios of the source-load uncertainty. The widely applied robust optimization includes conventional robust optimization (CRO) (Zeng and Zhao, 2013), distributionally robust optimization (DRO) (Ding et al., 2022), and chance-constrained robust optimization (Qi et al., 2020; Li et al., 2022). CRO formulates the source-load uncertainty with the uncertainty set, i.e., the maximum possible space of uncertain scenarios. A two-stage robust unit commitment based on it is established (Bertsimas et al., 2013; Zhou et al., 2019), and the day-ahead scheduling is usually solved with the Benders duality method (Yang and Su, 2021) or column and constraint generation (C&CG) methods (Wang et al., 2022). Moreover, the regulation of energy storages is incorporated in the unit commitment, and a nested C&CG is developed to address the integer state variables of energy storages in the intraday regulation (Cobos et al., 2018). CRO only takes the boundary information of the uncertainty into account, while ignores the distribution information, and the methods usually lead to a conservative solution. In view of this, DRO formulates the source-load uncertainty with the ambiguous set (Wei et al., 2016; Zheng et al., 2021), i.e., the maximum possible space of uncertainty distribution, and optimizes the expected operation cost in the worst-case distributions. However, chance-constrained robust optimization (Qi et al., 2022) assumes that the constraints hold with a certain confidence level. The method improves the flexibility of the scheduling, while the system security is affected.

The aforementioned robust methods can be concluded as the adaptive optimization decision-based robust optimization. In the second stage of the scheduling (i.e., the adjustment stage), the adjustment of regulation resources adapting to the immediate source-load uncertainty is determined through an optimization procedure. However, the system power balance in the real-time scheduling is restored through the automatic generation control (AGC) system (Wu et al., 2017). According to the frequency deviation and tie-line power deviation of the regional power system, the AGC system obtains the power fluctuation from the deviation signals and allocates the adjusted power to regulation resources proportionally. The AGC system has a very high requirement for the decision time, and hence, the aforementioned adaptive optimization decision-based robust methods are inapplicable.

In addition to the adaptive optimization decision-based robust methods, the linear decision-based CRO and DRO (LCRO and LDRO) are also proposed. When the source-load uncertainty occurs, the regulation tasks are allocated to regulation resources, according to the affine factors, i.e., allocation coefficients. Therefore, the linear decision-based robust methods adapt to the real-time scheduling. The co-regulation of generation units and energy storages to wind uncertainty is realized based on LCRO (Jabr 2013). To improve the applicability of LCRO, the allowed interval of wind power is optimized in (Li X et al., 2015; Wang et al., 2019). The allowed upper bound is lower than the forecast upper bound, and hence the strategy makes LCRO applicable for systems with inadequate regulation capacity. In view of the conservativeness of LCRO, which optimizes the operation at the expected scenario, the LDRO-based scheduling model is developed for power systems (Xiong et al., 2017; Alismail et al., 2018). Considering that the worst-case distributions of the source-load uncertainty rarely happen, LDRO may lead to a non-optimal economical effectiveness in common distributions. To enhance the economic efficiency of robust methods in common distributions, a linear decision-based stochastic robust optimization (LSRO) is proposed in (Qu et al., 2020). The method establishes the system constraints using the robust optimization to guarantee the system resilience, while optimizes the operation cost in the empirical distribution of uncertainties for the economical effectiveness. The empirical distribution of uncertainties is the fitting of the existing historical scenarios, and it is closest to the real distribution of uncertainties. Therefore, LSRO shows an economical effectiveness in common distributions of uncertainties.

2 Real-time power scheduling model with linear stochastic robust optimization

2.1 Source-load uncertainty set

The available wind power of wind farms and the real power demand equal the sum of forecast expectations and forecast errors as follows:where is the forecast available wind power of the wind farm m, is the forecast error of the wind power, is the forecast power demand of the node n and, v_n,t is the forecast error of the power demand.

Considering that the regulation resources may be insufficient, a certain amount of wind curtailment is allowed to ensure the system security. Therefore, the system will determine an allowed upper bound of total wind fluctuations (). is smaller than the forecast upper bound (). The superfluous wind power exceeding the allowed threshold will be curtailed as follows:where PW m,t is the baseline power output of the wind farm m, is the real wind power output, is the total available wind fluctuations, and U_t is the lower bound of wind fluctuations.

When considering the wind curtailment, the real wind fluctuation is smaller than the available wind fluctuation, and thus the system security constraints are only associated with the allowed wind fluctuation (i.e., real wind fluctuation). The uncertainty set of the allowed wind fluctuation is established aswhere U_t is the total allowed wind fluctuation and is the allowed upper bound of the total wind fluctuation.

As for the power demand, load shedding is not allowed, and thus the uncertainty set for the forecast error of the power demand is formulated aswhere V_t is the total power demand fluctuation.

Therefore, the source-load uncertainty set is a union set of the wind uncertainty set and power demand uncertainty set, which is given as: Ξ_t = ΞU t∪ΞV t.

2.2 Overall architecture of real-time scheduling

Real-time power scheduling optimizes the baseline power output of wind farms, thermal units, hydro units, and baseline charging power of energy storages, as well as the affine factors of non-wind regulation resources for the source-load uncertainty. The real-time power scheduling is usually performed in a rolling optimization way to determine the scheduling plan for the following 1 h period with 5 min resolution. When the wind power and power demand deviate from the baseline value, non-wind regulation resources including thermal units, hydro units, and energy storages co-regulate to restore the system power balance with the affine factors.where , and are the baseline power output of the thermal unit i, hydro unit j, and baseline charging power of the energy storage k; , , and are their real power outputs or charging power under the source-load uncertainty and; are downward/upward affine factors of non-wind regulation resources for the source-load uncertainty.

To sum up, the strategy for balancing the source-load uncertainty is presented as follows:

Power demand: Load shedding is not allowed, and power demands always trace the actual demands.

Wind farms: If the total available wind fluctuation is below the allowed upper bound , each wind farm operate at the available value, otherwise, the superfluous wind power beyond is curtailed, as shown in Eq. 2.

Non-wind resources: Thermal units, hydro units, and energy storages co-regulate as Eqs 5–7 to balance the source-load uncertainty.

2.3 System security constraints

The maximum/minimum operation of wind farms, thermal units, hydro units, and energy storages occur at the extreme total source-load fluctuation, and hence their constraints are established using the extreme scenario method. As for transmission lines, their transmitted power is related to the power output of each wind farm and power demand of each node, which cannot be described with the total source-load fluctuation. Therefore, we formulate the transmission limits of lines as that concerning the source-load uncertainty as follows:where is the baseline transmitted power of line l, h_l,i, h_l,j, h_l,k, h_l,m, and h_l,n are transfer distribution factors from the thermal unit i, hydro unit j, energy storage k, wind farm m, and power demand n to the line l, and are transmitted power limits of line l. Ξ_t,1 and Ξ_t,2 are two subsets of the source-load uncertainty set. In Ξ_t,1, U_t-V_t≤0, the system adopts the downward affine factors to balance the uncertainty, while in Ξ_t,2, U_t-V_t≥0, the system adopts the upward affine factors to balance the uncertainty.

The transmission limits of lines can be described as the following compact form:where A, B, C are constant vectors, d is a constant scalar, x is the decision variable vector, α_j is the affine factor matrix, concerning the jth uncertainty subset, and π, λ, τ, and ϑ are dual variables, corresponding to constraints in the uncertainty subset Ξ_t,j.

By adopting the dual theory (Li et al., 2018), the constraint (29) with the source-load uncertainty is converted to a deterministic constraint as follows:where e is a vector whose elements are all equal 1.

2.4 Optimized objective

The scheduling optimizes the baseline operation cost (i.e., the first stage) and the incremental cost expectation in the adjustment state (i.e., the second stage) as follows:where F_base is the baseline operation cost, is the empirical distribution of the source-load uncertainty, ΔF is the incremental adjustment cost under the source-load uncertainty .

This paper introduces carbon trading (Wang et al., 2020a) to improve the environmental benefits of the resilient power system. The operation cost in the baseline state includes fuel costs of thermal units, system carbon trading costs, and penalty costs on the wind curtailment.where C_i() is the quadratic fuel cost function of the thermal unit i, f_c is the carbon trading costs of the system, and W() is the quadratic penalty cost function of the wind curtailment.

This study adopts the baseline method (Wang et al., 2020b; Guo et al., 2020) to calculate the carbon quota of the system, and the carbon trading cost of the system is formulated as follows:where p_c is the carbon trading price, E_R and E_Q are the baseline carbon emission and carbon quota of the system, respectively, β_i,t is the carbon emission coefficient of the thermal unit i, anf η_c is the carbon quota for per unit of power generation.

The incremental cost expectation in the adjustment state includes the adjustment cost of thermal units and energy storages, incremental penalty costs on the wind curtailment, and incremental carbon trading costs. The traditional random sampling methods to formulate the output expectation usually requires a huge sampling scale to achieve an accurate estimation. In this paper, the multi-point estimation (Li Z et al., 2015) can uniformly cover the source-load uncertainty with a small scale of estimation points.

The multi-point estimation actually first determines the estimate vectors in the standard Gaussian distribution, then converts these estimate vectors to that of the source-load uncertainty, and finally calculates the output expectation (Li et al., 2008). The optimal estimate point number will be determined through sensitivity analysis in numerical simulations. The source-load uncertainty is assumed to be irrelevant in the time dimension, and thus only two uncertain variables are considered in the formulation of the incremental cost expectation, i.e. the available wind fluctuation and forecast error of the power demand. The incremental cost expectation in the adjustment state is formulated aswhere r∈N is the index of uncertain variables, s∈M is the index of estimate points, ω_s is the weight of the sth estimate point, (r, s) corresponds to the estimate vector concerning the rth uncertain variable and the sth estimate point, and z_r,s is the estimate vector in the standard Gaussian distribution. In the estimate vector, the uncertain variable r equals the sth estimate point and other uncertain variables equal 0. Ψ is the cumulative distribution function of the standard Gaussian distribution. F_t⁻¹ is the inverse of the cumulative distribution function of the source-load uncertainty.

Assume that the sign variable concerning the adopted affine factors in each estimate vector is σ_t,r,s. If σ_t,r,s = 1, the source-load power fluctuation is larger than 0, the system adopts the upward affine factors for the adjustment; else if σ_t,r,s = 0, the system adopts the downward affine factors for the adjustment. In each estimate vector, the incremental cost produced by thermal units and energy storages are calculated aswhere U_t,r,s is the allowed wind fluctuation in the estimate vector (r, s), ρ_i,t, and p_i,t are the incremental cost coefficient and regulation cost coefficient of the thermal unit i, respectively, and ρ− k, t ,and ρ+ k,t are the downward and upward incremental cost coefficient of the energy storage k. Since hydro units do not produce fuel costs and carbon emissions in the operation, their incremental cost coefficients equal 0. As for energy storages, since they also play a role of peak-load shifting and standby power generation during intraday operation, the stored power in storages after the scheduling shall not deviate too far from the required value. Therefore, the incremental cost coefficients are set as follows to reduce the adjustment of energy storages: ρ− k,t = -1.1 × max_i (ρ_i,t), ρ+ k,t = -0.9 × max_i (ρ_i,t).

To sum up, the incremental cost in each estimate vector is calculated aswhere ρ_w is the penalty coefficient on the wind curtailment.

3 Solution methodology

The proposed scheduling model is a mixed integer bilinear constraint programming. The integer variables include the sign variables σ_t,r,s in each estimate vector concerning the incremental cost expectation and the charging/discharging state variables of energy storages. The bilinear constraints include Eq. 11 (15) (23) (30). Due to the timeliness required in the real-time scheduling, the integer variables and bilinear constraints shall be addressed, and such that the scheduling model can be converted to a convex quadratic programming.

3.1 Approximation for the sign variables

Therefore, the initial model does not contain the sign variables σ_t,r,s and bilinear constraints. Compare the available wind power in each estimate vector () with the approximated value for the allowed upper bound of wind fluctuation () and obtain the sign variables σ_t,r,s in each estimate vector.

3.2 Linearization for bilinear constraints

This paper adopts the AO (Li et al., 2021) to address the bilinear constraints in the real-time scheduling, which fixes partial decision variables to optimize the rest of the decision variables in the bilinear constraints. The real-time scheduling with bilinear constraints can be described as the following compact form:where A_i is a constant matrix, g_i is a convex function concerning decision variables, y corresponds to the allowed upper bound of wind fluctuation (), and x corresponds to the remaining decision variables. Detailed steps for AO to solve the real-time scheduling model are presented in Algorithm 1.

3.3 Approximation for state variables of energy storages

This paper adopts an “estimation–correction” strategy (Shi and Guo, 2022) to approximate the state variables of energy storages. The estimation stage relaxes the state variables of energy storages, and the correction stage directly determines the state variables according to the results obtained in the estimation stage.

Therefore, the scheduling problems in these two stages do not contain the state variables of energy storages. The relaxed model of energy storages in the estimation stage is formulated aswhere are the charging power and discharging power of energy storages in the baseline scenario and two extreme scenarios.

The state variables of energy storage are directly fixed in the correction stage. For the sake of convenience, the following takes the baseline scenario as an example, while the other two extreme scenarios can be disposed in a similar way.where are the charging power and discharging power of energy storages obtained in the estimation stage.

In fact, there is generally little difference between the optimal state variables of energy storages in the estimation stage and the correction stage. On the one hand, the charging and discharging efficiency of energy storages are generally higher. In the real-time scheduling, only the energy storages with the rapid adjustment ability could be involved in the AGC system, such as the lithium storages, and the efficiency of lithium storages can reach about 90% at present (Shi and Guo, 2022). On the other hand, the current capacity of energy storages is generally small compared with the total power demand, and thus the relaxed model of energy storages in the estimation stage has slight influence on the optimized result.

The overall flowchart for the convexified solution is presented in Figure 1.

FIGURE 1

4 Numerical simulation

To illustrate the effectiveness of the proposed real-time low-carbon scheduling and solving method, numerical simulations are carried out on the IEEE 39-node power system (Zimmerman et al., 2011). The system has seven thermal units, one hydro unit, and two energy storages with 50MW capacity. Moreover, four wind farms are deployed. The scheduling cycle is from 6:00 pm to 7:00 pm with 5 min time resolution. The forecast expectation and forecast interval of wind power and power demand are presented in Figure 2 and Figure 3. Finally, 5,000 scenarios of the source-load uncertainty are sampled by Monte Carlo simulation (MCS) to verify the scheduling performance. The optimizations are solved with the CPLEX solver on a computer with Intel(R) Xeon(R) Silver 4216 CPU, 2.1GHz clock speed, and 16GB RAM.

FIGURE 2

FIGURE 3

4.1 Effects of the allowed upper bound of wind fluctuation

In the proposed scheduling model, the allowed upper bound of wind fluctuation actually has two functions: 1) it reasonably determines the wind power accommodation and hence enhances the system flexibility, 2) when the system regulation capacity is insufficient, discard the excess wind power to ensure the system security. Figure 4 and Figure 5 present the allowed upper bounds (UB in the figures) of wind fluctuation in the system with sufficient and insufficient regulation capacity, where the maximum capacity of wind power accommodation is the obtained allowed upper bound of wind fluctuation when adding the wind penalty objective, i.e., the second line in Eq. 42, to the original objective. Moreover, in the system with insufficient regulation capacity, change the ramp rates of thermal units and hydro units as well as the maximum charging and discharging power of energy storages to 75% of the original.

FIGURE 4

FIGURE 5

As shown in Figure 4, although the system with sufficient regulation capacity can completely accommodate wind power, partial wind power is still discarded. This is mainly because the wind fluctuation exceeding the allowed upper bound rarely appears, and hence these wind powers are discarded to improve the system flexibility. In addition to, as shown in Figure 5, the system with insufficient regulation capacity cannot completely accommodate wind power, and hence the allowed upper bound is below the forecast upper bound to ensure the system security. To further validate the effects of the allowed upper bound in improving the system flexibility, Table 1 compares the optimization results and MCS results with two wind power accommodation strategies. The maximum accommodation strategy is to absorb wind power as far as possible, while the flexible accommodation strategy is that in the proposed scheduling model. As shown in Table 1, although the allowed upper bound of wind fluctuation in the flexible strategy is obviously smaller than the maximum strategy, the discarded wind power is very slight. This is because the wind fluctuation exceeding the allowed upper bound rarely appears. Through the flexible accommodation strategy, the system operation cost is reduced, which validates the effects of allowed upper bound in improving the system flexibility.

4.2 Baseline state and affine factors of regulation resources

Figure 6 presents the baseline outputs of thermal units and hydro units, and the baseline charging power of energy storages. Hydro units always operate at the maximum output since they do not produce fuel costs and carbon trading costs. Thermal units produce fuel costs and carbon trading costs, and hence, they arrange the outputs according to the cost-optimal principle. In the scheduling cycle of the testing case, the load peak is in the front period, while the load trough is in the latter period. Therefore, the energy storage is discharged according to the required state in the front period to realize the peak load shifting and reduce the operation cost.

FIGURE 6

Figure 7 and Figure 8 present the affine factors of thermal units, hydro units, and energy storages for the downward and upward source-load uncertainty, respectively. Affine factors of regulation resources are closely related to their incremental cost coefficients. When the source-load uncertainty fluctuates downward, the system operation cost is increased. The incremental cost coefficients of regulation resources satisfy , and hence the regulation priority for the downward source-load uncertainty satisfies: hydro units > thermal units > energy storages. Since hydro units operate at the maximum output and have no room for the upward adjustment, they do not undertake the regulation task for the downward source-load uncertainty, as shown in Figure 7. During periods t = 1–6, the source-load uncertainty is relatively small, thermal units have sufficient regulation capacity, and hence the regulation tasks for the downward source-load uncertainty are fully undertaken by the thermal units. While during periods t = 7–12, the regulation capacity of the thermal units is insufficient, and hence energy storages participate in the regulation for the downward source-load uncertainty. When the source-load uncertainty fluctuates upward, the system operation cost is decreased. The incremental cost coefficients of regulation resources satisfy , and hence the regulation priority for the upward source-load uncertainty satisfies thermal units > energy storages > hydro units. As shown in Figure 8, during periods t = 1–7, thermal units completely undertake the regulation task for the upward source-load uncertainty, while energy storages and hydro units participate in the regulation for the upward source-load uncertainty during next periods t = 8–12. As indicated in the simulation results, the LSRO reasonably arranges the affine factors, according to the incremental cost coefficients of regulation resources, and thus it can reduce the system incremental cost in the adjustment state.

FIGURE 7

FIGURE 8

Moreover, the existing AGC system usually utilizes units for the regulation. Due to the rapid adjustment speed, energy storages can participate in the AGC system to regulate the source-load uncertainty. To validate the effects of energy storages in the regulation, Table 2 compares the optimization results and MCS results without and with the regulation of energy storages.

4.3 Comparison of different robust methods

In order to further verify the economic efficiency of the LSRO, we compare it with LCRO and LDRO. The robust methods all optimize the allowed upper bound of wind fluctuation. The difference of these three robust methods is that, LSRO optimizes the expected cost under the empirical distribution of the source-load uncertainty, LCRO minimizes the baseline cost and maximizes the allowed upper bound of wind fluctuation, and LDRO minimizes the expected cost under the worst-case distributions of the source-load uncertainty. As shown in Table 3, LCRO minimizing the baseline cost ignores the impact of the adjustment process on the operation cost, and hence its adjustment is uneconomical. LDRO minimizes the expected cost under the worst-case distributions, however, the worst-case distributions rarely appear, and hence LDRO actually sacrifices the economic efficiency in general distributions for robustness in worst cases. The empirical distribution is the fitting of the historical scenarios, which is closest to the real distribution of the source-load uncertainty, and hence LSRO shows the best economic efficiency.

4.4 Impacts of estimate point numbers

This paper upgrades the traditional three-point estimation to multi-point estimation. The estimate point numbers are set at 3, 5 and 7, respectively and results with different numbers are compared. Figure 9 shows the allowed upper bound of wind fluctuation obtained by multi-point estimation with different point numbers. The maximum cumulative distribution function covered by estimate points in three-point estimation, five-point estimation, and seven-point estimation are 0.9584, 0.9979, and 0.9999, respectively. The distribution of uncertain variables covered by three-point estimation is limited, while the uncertainty distribution covered by five-point estimation and seven-point estimation is large enough. Therefore, the allowed upper bound of wind fluctuation obtained by three-point estimation is much smaller compared to that obtained by five-point estimation and seven-point estimation. In addition, Table 4 compares the results of optimization and MCS with different point numbers, and it can be seen that there is little difference between the results. Considering that the seven-point estimation can cover the maximum distribution of uncertain variables, and the approximate degree of the expected value evaluation generally enhances with the increase in the point number (Li X et al., 2015), the estimate point number is set at 7 in this paper to ensure the robustness of the multi-point estimation method.

FIGURE 9

4.5 Validation of solving methods

As indicated in subsection 3.2, AO is a sub-optimal method, whose optimization is related to initialization. This paper evenly sets the confidence level corresponding to initialization and updates the initial value of the allowed upper bound of wind fluctuation as Eq. 51. The results of optimization and MCS associated with the confidence level are presented in Table 5. Generally, the obtained allowed upper bound of wind fluctuation decreases with decrease in the confidence level. When the confidence level is high, the obtained allowed upper bound is large, and hence the wind curtailment is small. However, the system flexibility is affected due to excessive wind power accommodation, and thus the system operation cost increases. As the confidence level decreases, the obtained allowed upper bound decreases, the system flexibility improves, and the system operation cost is reduced. However, when the confidence level reaches 97.5%, the obtained allowed upper bound is too small, and hence the wind power accommodation is inadequate and the system operation cost gradually increases. Considering that the wind power prediction accuracy has a certain degree of error, the confidence level can be set at 99.7% to absorb wind power as far as possible without compromising the system’s flexibility.

In addition, the integer variables of charging/discharging state of energy storages are addressed with the estimation–correction strategy. Table 6 compares the optimization results with different strategies for the integer variables. The precise strategy takes the state variable of energy storage as variables in the optimization, and hence the scheduling problem is formulated as an iterative mixed integer quadratic constraint programming, which takes a long time to solve. The estimation stage relaxes the state variables, and the correction stage directly fixes the state variables as parameters, according to the results obtained in the estimation stage. Therefore, the scheduling problem is formulated as an iterative quadratic constraint programming with the estimation–correction strategy, and the computation time is short. On the other hand, the theoretically optimal objectives of these three strategies satisfy the following relations: F_estimation < F_precise < F_correction. Since the deviation between the optimization results of the estimation stage and the correction stage is small, it can be deduced that the estimation–correction strategy can find a good enough solution.

4.6 Impacts of carbon trading price

Carbon trading cost is an important part of the proposed real-time low-carbon scheduling. To explore the influence of the carbon trading price on the scheduling, this subsection changes the carbon trading price from 0 to 60$/t and obtains the variation trend of the operation cost and carbon emission of the system relative to the carbon trading price. As indicated in Figure 10, carbon trading has little impact on the system operation when the carbon trading price is low, and thus the operation cost and carbon emissions change slightly. With the increase in the carbon trading price, the proportion of the carbon trading cost in the operation cost increases, the system gradually utilizes the low-carbon resources for power generation, and the carbon emission show a significant downward trend. When the carbon trading price rises to 26$/t, the system can make profits by the carbon trading mechanism, and the operation cost gradually decreases. When the carbon trading price rises to 40$/t, the carbon reduction space of the system becomes smaller, and the downward trend of carbon emissions slows down.

FIGURE 10

5 Simulation in large-scale systems

5.1 Simulation results

To further verify the effectiveness of LSRO, supplementary simulations are performed in a large-scale IEEE 118-node power system. The system has seventeen thermal units, two hydro units, and four energy storages with 40MW capacity. Furthermore, six wind farms are deployed. As shown in Table 7, the results in the IEEE 118-node system are similar to those in the IEEE 39-node system. LCRO focuses on the baseline cost while ignores the influence of the adjustment on the operation cost, and the adjustment of the method is uneconomical. LDRO minimizes the expected cost under the worst-case distribution, in view of the extremely low probability of the worst-case distribution, and the method is not effective enough in common distributions. Compared to LCRO and LDRO, LSRO minimizes the expected cost under the empirical distribution of the source-load uncertainty. In view of that, the empirical distribution is closest to the real distribution of the source-load uncertainty, and LSRO shows the best economic efficiency.

5.2 Further discussion on the robust method

In fact, the performance of LSRO is related to the accuracy of the empirical distribution of the source-load uncertainty. If the scale of historical scenarios is not huge enough or the adopted fitting method to determine the empirical distribution has a poor performance, the empirical distribution may deviate far from the real distribution. Therefore, the performance of the robust method in other distributions is affected in these cases. To improve the robustness of LSRO, our future work will combine the empirical distribution and worst-case distribution together in the robust method.

6 Conclusion

This paper proposes a real-time low-carbon scheduling for the wind–thermal–hydro-storage resilient power system based on the linear decision rule-based stochastic robust optimization. A multi-step alternate optimization-based convexified solution is formulated to solve the complex scheduling problem. The allowed upper bound of wind fluctuation provides the system with more flexibility. Monte Carlo simulation shows that very slight redundant wind power is sacrificed to significantly reduce the system operation cost. By the stochastic robust optimization, the regulation task to regulation sources is reasonably arranged according to the incremental cost coefficients of resources, and thus an economically effective adjustment is realized. The comparison of the adopted stochastic robust optimization with conventional robust optimization and distributionally robust optimization in two case studies validates the economic efficiency of the method. Numerical simulations indicate that the performance of alternate optimization is related to initialization, and the confidence level concerning the initialization can be reasonably set to absorb wind power as far as possible without compromising the system flexibility.

The empirical distribution of the source-load uncertainty may not be an accurate estimation of the real distribution, when the number of historical scenarios is inadequate or the fitting method is inaccurate. Our future work will combine the empirical distribution and worst-case distribution to improve the applicability of the linear decision-based robust method.

References

• 1

  AlismailF.XiongP.SinghC. (2018). Optimal wind farm allocation in multi-area power systems using distributionally robust optimization approach. IEEE Trans. Power Syst.33 (1), 536–544. 10.1109/TPWRS.2017.2695002

  • CrossRef
  • Google Scholar

• 2

  BertsimasD.LitvinovE.SunX.ZhaoJ.ZhengT. (2013). Adaptive robust optimization for the security constrained unit commitment problem. IEEE Trans. Power Syst.28 (1), 52–63. 10.1109/TPWRS.2012.2205021

  • CrossRef
  • Google Scholar

• 3

  CobosN. G.ArroyoJ. M.AlguacilN.WangJ. (2018). Robust energy and reserve scheduling considering bulk energy storage units and wind uncertainty. IEEE Trans. Power Syst.33 (5), 5206–5216. 10.1109/TPWRS.2018.2792140

  • CrossRef
  • Google Scholar

• 4

  DingX.MaH.YanZ.XingJ.SunJ. (2022). Distributionally robust capacity configuration for energy storage in microgrid considering renewable utilization. Front. Energy Res.10, 923985. 10.3389/fenrg.2022.923985

  • CrossRef
  • Google Scholar

• 5

  GuoZ.LiG.ZhouM. (2020). Fast and dynamic robust optimization of integrated electricity-gas system operation with carbon trading. Power Syst. Technol.44 (4), 1220–1228. 10.13335/j.1000-3673.pst.2019.2313

  • CrossRef
  • Google Scholar

• 6

  HeR.DongP.LiuM.HuangS.LiJ. (2022). A review of wind energy output simulation for new power system planning. Front. Energy Res.10, 942450. 10.3389/fenrg.2022.942450

  • CrossRef
  • Google Scholar

• 7

  JabrR. A. (2013). Adjustable robust OPF with renewable energy sources. IEEE Trans. Power Syst.28 (4), 4742–4751. 10.1109/TPWRS.2013.2275013

  • CrossRef
  • Google Scholar

• 8

  LiH.LvZ.YuanX. (2008). Nataf transformation based point estimate method. Chin. Sci. Bull.53 (6), 2586–2592. 10.1007/s11434-008-0351-0

  • CrossRef
  • Google Scholar

• 9

  LiX.ZhangX.WuL.LuP.ZhangS. (2015). Transmission line over-load risk assessment for power systems with wind and load-power generation correlation. IEEE Trans. Smart Grid6 (2), 1233–1242. 10.1109/TSG.2014.2387281

  • CrossRef
  • Google Scholar

• 10

  LiZ.WuW.ZhangB.WangB. (2015). Adjustable robust real-time power dispatch with large-scale wind power integration. IEEE Trans. Sustain. Energy6 (2), 357–368. 10.1109/TSTE.2014.2377752

  • CrossRef
  • Google Scholar

• 11

  LiP.YuD.YangM.WangJ. (2018). Flexible look-ahead dispatch realized by robust optimization considering CVaR of wind power. IEEE Trans. Power Syst.33 (5), 5330–5340. 10.1109/TPWRS.2018.2809431

  • CrossRef
  • Google Scholar

• 12

  LiP.YangM.WuQ. (2021). Confidence interval based distributionally robust real-time economic dispatch approach considering wind power accommodation risk. IEEE Trans. Sustain. Energy12 (1), 58–69. 10.1109/TSTE.2020.2978634

  • CrossRef
  • Google Scholar

• 13

  LiZ.WuL.XuY.ZhengX. (2022). Stochastic-weighted robust optimization based bilayer operation of a multi-energy building microgrid considering practical thermal loads and battery degradation. IEEE Trans. Sustain. Energy13 (2), 668–682. 10.1109/TSTE.2021.3126776

  • CrossRef
  • Google Scholar

• 14

  LiZ.WuL.XuY.WangL.YangN. (2023). Distributed tri-layer risk-averse stochastic game approach for energy trading among multi-energy microgrids. Appl. Energy331, 120282. 10.1016/j.apenergy.2022.120282

  • CrossRef
  • Google Scholar

• 15

  LiZ.XuY. (2019). Temporally-coordinated optimal operation of a multi-energy microgrid under diverse uncertainties. Appl. Energy240, 719–729. 10.1016/j.apenergy.2019.02.085

  • CrossRef
  • Google Scholar

• 16

  LiuY.ShenZ.TangX.LianH.LiJ.GongJ. (2019). Worst-case conditional value-at-risk based bidding strategy for wind-hydro hybrid systems under probability distribution uncertainties. Appl. Energy256, 113918. 10.1016/j.apenergy.2019.113918

  • CrossRef
  • Google Scholar

• 17

  LiuH.MiaoZ.WangN.YangY. (2022). A two-layer optimization of design and operational management of a hybrid combined heat and power system. Front. Energy Res.10, 959774. 10.3389/fenrg.2022.959774

  • CrossRef
  • Google Scholar

• 18

  QiF.ShahidehpourM.LiZ.WenF.ShaoC. (2020). A chance-constrained decentralized operation of multi-area integrated electricity-natural gas systems with variable wind and solar energy. IEEE Trans. Sustain. Energy11 (4), 2230–2240. 10.1109/TSTE.2019.2952495

  • CrossRef
  • Google Scholar

• 19

  QuK.YuT.PanZ.ZhangX. (2020). Point estimate-based stochastic robust dispatch for electricity-gas combined system under wind uncertainty using iterative convex optimization. Energy211, 118986. 10.1016/j.energy.2020.118986

  • CrossRef
  • Google Scholar

• 20

  ShiY.GuoC. (2022). Flexibility reinforcement method for integrated electricity and heat system based on decentralized and coordinated multi-stage robust dispatching. Automation Electr. Power Syst.46 (6), 10–19. 10.7500/AEPS20210311013

  • CrossRef
  • Google Scholar

• 21

  SurenderR. S.BijweP. R.AbhyankarA. R. (2015). Real-time economic dispatch considering renewable power generation variability and uncertainty over scheduling period. IEEE Syst. J.9 (4), 1440–1451. 10.1109/JSYST.2014.2325967

  • CrossRef
  • Google Scholar

• 22

  WangC.WeiW.WangJ.BiT. (2019). Convex optimization based adjustable robust dispatch for integrated electric-gas systems considering gas delivery priority. Appl. Energy239, 70–82. 10.1016/j.apenergy.2019.01.121

  • CrossRef
  • Google Scholar

• 23

  WangY.QiuJ.TaoY.ZhaoJ. (2020a). Carbon-oriented operational planning in coupled electricity and emission trading markets. IEEE Trans. Power Syst.35 (4), 3145–3157. 10.1109/TPWRS.2020.2966663

  • CrossRef
  • Google Scholar

• 24

  WangY.QiuJ.TaoY.ZhangX.WangG. (2020b). Low-carbon oriented optimal energy dispatch in coupled natural gas and electricity systems. Appl. Energy280, 115948. 10.1016/j.apenergy.2020.115948

  • CrossRef
  • Google Scholar

• 25

  WangY.QiuJ.TaoY. (2022). Robust energy systems scheduling considering uncertainties and demand side emission impacts. Energy239, 122317. 10.1016/j.energy.2021.122317

  • CrossRef
  • Google Scholar

• 26

  WeiW.LiuF.MeiS. (2016). Distributionally robust co-optimization of energy and reserve dispatch. IEEE Trans. Sustain. Energy7 (1), 289–300. 10.1109/TSTE.2015.2494010

  • CrossRef
  • Google Scholar

• 27

  WuH.KradI.FloritaA.HodgeB.IbanezE.ZhangJ.et al (2017). Stochastic multi-timescale power system operations with variable wind Generation. IEEE Trans. Power Syst.32 (5), 3325–3337. 10.1109/TPWRS.2016.2635684

  • CrossRef
  • Google Scholar

• 28

  XiongP.JirutitijaroenP.SinghC. (2017). A distributionally robust optimization model for unit commitment considering uncertain wind power generation. IEEE Trans. Power Syst.32 (1), 39–49. 10.1109/TPWRS.2016.2544795

  • CrossRef
  • Google Scholar

• 29

  YangJ.SuC. (2021). Robust optimization of microgrid based on renewable distributed power generation and load demand uncertainty. Energy223, 120043. 10.1016/j.energy.2021.120043

  • CrossRef
  • Google Scholar

• 30

  ZengB.ZhaoL. (2013). Solving two-stage robust optimization problems using a column-and-constraint generation method. Oper. Res. Lett.41 (5), 457–461. 10.1016/j.orl.2013.05.003

  • CrossRef
  • Google Scholar

• 31

  ZhengX.QuK.LvJ.LiZ.ZengB. (2021). Addressing the conditional and correlated wind power forecast errors in unit commitment by distributionally robust optimization. IEEE Trans. Sustain. Energy12 (2), 944–954. 10.1109/TSTE.2020.3026370

  • CrossRef
  • Google Scholar

• 32

  ZhouA.YangM.WangZ.LiP. (2018). A linear solution method of generalized robust chance constrained real-time dispatch. IEEE Trans. Power Syst.33 (6), 7313–7316. 10.1109/TPWRS.2018.2865184

  • CrossRef
  • Google Scholar

• 33

  ZhouB.AiX.FangJ.YaoW.ZuoW.ChenZ.et al (2019). Data-adaptive robust unit commitment in the hybrid AC/DC power system. Appl. Energy254, 113784. 10.1016/j.apenergy.2019.113784

  • CrossRef
  • Google Scholar

• 34

  ZhuQ.WangY.SongJ.JiangL.LiY. (2021). Coordinated frequency regulation of smart grid by demand side response and variable speed wind turbines. Front. Energy Res.9, 754057. 10.3389/fenrg.2021.754057

  • CrossRef
  • Google Scholar

• 35

  ZimmermanR. D.Murillo-SanchezC. E.ThomasR. J. (2011). Matpower: Steady-state operations, planning, and analysis tools for power systems research and education. IEEE Trans. Power Syst.26 (1), 12–19. 10.1109/TPWRS.2010.2051168

  • CrossRef
  • Google Scholar