Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Energy Res., 18 August 2025

Sec. Smart Grids

Volume 13 - 2025 | https://doi.org/10.3389/fenrg.2025.1612065

Coordinated multi-level scheduling method considering uncertainty of renewable energy and load

Bingbing SongBingbing Song1Kangping QinKangping Qin1Min WenMin Wen1Moyan Zhu
Moyan Zhu2*Kaiming ZouKaiming Zou2Guangyu HeGuangyu He2
  • 1East China Branch of State Grid Corporation of China, Shanghai, China
  • 2Key Laboratory of Control of Power Transmission and Conversion, Ministry of Education Shanghai Jiao Tong University, Shanghai, China

As renewable energy continues to be widely integrated, the energy structure is gradually transforming. The increasing grid connection of wind and photovoltaic power signifies a major shift in the energy mix. This change is particularly evident in heavy load areas at the regional grid and provincial dispatch levels, where uncertainties on both the supply and demand sides impact the daily operation of power systems. New dispatch strategies are urgently needed to address these uncertainties. This paper introduces a two-stage day-ahead and intra-day coordinated multi-level dispatch method that considers both the regional-level and provincial-level power systems, addressing supply-demand uncertainties from the perspective of regional grid-level and unmet load peak shaving. Unmet load refers to the load that cannot be met solely by the output of regional grid units. At the regional grid level, a unit dispatch model for unmet load peak shaving is developed. We introduce the concept of unmet load and, based on peak-valley weighting, propose a multi-province load peak shaving method, improving the approach to unmet load considerations. At the provincial level, a two-stage robust optimization dispatch model is constructed based on regional grid dispatch, and it is solved using the Karush-Kuhn-Tucker conditions and the Column-and-Constraint Generation (C&CG) algorithm. Finally, case study results validate the proposed model’s effectiveness, demonstrating its ability to provide an optimized coordinated grid-provincial dispatch strategy under supply-demand uncertainty.

1 Introduction

In the current context of widespread integration of renewable energy, considering the high degree of new energy integration in large-scale grid dispatch, the system is prone to fluctuations due to load and the variability in wind and solar power outputs, and is more susceptible to peak load situations, increasing the risk to power safety. It is particularly crucial to consider dispatch strategies at both regional grid and provincial levels to ensure power balance during peak load periods and to coordinate regional and provincial dispatch strategies. This aims to ensure stable power supply at the regional grid level and promote the consumption of new energy at the provincial level, thereby maximizing the safety of large-scale power grids.

At the regional grid level, the impact of new energy is significantly amplified due to the integration of multi-provincial loads, making it highly likely to encounter peak load situations during daily operations. Taking the East China Grid as an example, which mainly supplies power to Shanghai, Jiangsu, Zhejiang, Fujian, and Anhui, the daily power supply pressure is substantial, and it frequently participates in inter-provincial power dispatch (Li et al., 2024). Considering the limitations of line capacity and the cost of power adjustments in conventional units, using energy storage systems for cross-time scale power dispatch is an effective solution. Energy storage offers flexibility, efficiency, adjustability, rapid response, and environmental friendliness. Liu and Peng (2024) proposes a peak load transfer optimization model for wind-power-energy hybrid energy system based on situational awareness theory. Uddin et al. (2018) discusses on possible challenges and future research directions for each type of the strategy. Chua et al. (2016) provide an effective sizing method and an optimal peak shaving strategy for an energy storage system to reduce the electrical peak demand of the customers. Luo et al. (2024), Zhao et al. (2024), Cheng et al. (2018); Liao et al. (2024) construct different short-term peak-shaving frameworks to address the modeling challenge and optimization difficulty. Wallberg et al. (2024) presents a control algorithm that uses a negative correlation scheme, adjusted to the local grid load, to effectively manage the battery energy storage. Jin et al. (2022) uses flexible hydropower to buffer the volatility and the randomness of RE sources and aid peak shaving in response to the transition towards sustainability. Wang et al. (2021) proposes a nonlinear programming model to solve this problem. The research on energy storage for peak shaving is well-established, yet there is limited focus on regional grid-level dispatch scenarios. This paper takes into account the common issue of insufficient unit output at the regional grid level and considers the resulting unmet load situations. By utilizing energy storage systems to perform peak shaving and valley filling across multiple provinces, this study extends the objectives of peak shaving and valley filling to include reducing the fluctuations in unmet load.

At the provincial dispatch level, the main responsibility lies in accommodating the consumption of new energy after implementing unmet load peak shaving at the regional grid level. Dong et al. (2024) studies the operation and scheduling problem of virtual power plant with the collaborative optimization of multiple flexible loads and new energy, and improves the mismatch between power supply and demand through the efficient aggregation and optimal control of new energy and demand-side resources. Yang et al. (2024) designs a two-stage scheduling optimization framework to minimize the operating cost in the day-ahead phase and the system deviation cost in the intra-day phase. However, given the increasingly complex operational environment and various uncertainties, traditional optimization methods often fall short of meeting practical operational needs. The two-stage robust optimization method, an emerging optimization technology, has garnered considerable attention in recent years. This method has broad applicability in the field of power systems, whether in power market design, generation dispatch, or grid planning. By incorporating the two-stage robust optimization approach, complex issues in actual operations can be effectively addressed. Especially in the areas of new energy grid integration and inter-regional power trading, this method can effectively overcome challenges posed by uncertainties, ensuring efficient operation of the power system and optimal resource allocation. Zhao and Guan (2015) develops stochastic optimization models and solution methods to improve reliability unit commitment run practice. Büsing and Schmitz (2024), Zhu et al. (2024), Niu et al. (2022), Wang et al. (2022), Kong et al. (2022) study different two-stage robust optimization method under uncertainty of different variables. Zeng and Zhao (2013) present a column-and-constraint generation algorithm to solve two-stage robust optimization problems. Yang et al. (2023) proposes a distributed robust optimal scheduling method for microgrid based on discrete scenarios. Bendotti et al. (2023) proposes the anchor-robust approach as a middle ground between guaranteeing starting times and guaranteeing the thought.

This paper primarily introduces a two-stage robust optimization model under the coordination of regional grid-level dispatch. Based on regional grid dispatch strategies, it considers provincial dispatch under the influence of new energy sources. It also takes into account constraints such as tie-line power and inter-provincial electricity trading, aiming to minimize economic costs. The solution is obtained through iterative calculations.

2 Materials and methods

Figure 1 shows the model architecture in this document.

Figure 1
Flowchart illustrating a two-level dispatch model. The top section shows the Regional Grid-Level Dispatch Model, focusing on minimizing the variance of unmet load with variables for peak periods and penalty functions. The bottom section depicts the Provincial-Level Dispatch Model, including a master problem with variables for energy storage and power trading plans. It features iterative processes to update parameters, constraints, and bounds, utilizing KKT conditions and addressing worst-case scenarios with variables for power and frequency. Arrows indicate the flow between components and stages.

Figure 1. Model framework diagram.

2.1 Model formulation

2.1.1 Regional grid-level dispatch model

The dispatch optimization problem at the regional grid level studied in this paper involves two types of units: conventional generation units and pumped-storage power stations. The main constraints can be categorized as follows:

Controllable Conventional Unit Constraints:

Ramping Constraint Equation 1:

RUgPg,p,t+1Pg,p,tRUg(1)

In the formula, t denotes the time period index, and p represents the province index. Pg,p,t+1,Pg,p,t are the generation power at time t+1 and t respectively. RUg represents the ramp rate limit of the thermal power unit, indicating the maximum increase or decrease in output per unit time.

Unit Output Boundary Constraint Equation 2:

0Pg,p,tPgmax(2)

In the formula, Pgmax denotes the maximum output of the unit, and the constraint is defined with a zero lower bound, which applies to the generation stage of all units.

Pumped-Storage Unit Constraints Equations 37:

Pumped-Storage Unit Equality Constraints, including water level conversion power constraint and power balance constraint:

ηchargeΔE=Pcharget=ηchargemgΔh(3)
ΔE=ηdischargePdischarget=ηdischarge×mg+12ρV2S×Δh(4)
Et+1=Et+Ptchargeηcharge+PtdischargeηdischargeΔt(5)

The boundary constraints of pumped -storage units include energy boundary constraints as well as pumping and generating power constraints:

EminEtEmax(6)
0PtchargextchargePmaxchargextdischargePmaxdischargePtdischarge0xtcharge+xtdischarge1(7)

In the formula, Et represents the equivalent stored energy of the pumped-storage unit at time period t; ΔE denotes the equivalent energy change corresponding to the variation in water level; m,g,h,ρ,V,S respectively refer to the mass of water in the reservoir, gravitational acceleration, water level height, water density, volume, and bottom area; Ptcharge,Ptdischarge are the pumping and generating power of the pumped-storage unit at time period t; ηcharge,ηdischarge are the efficiency of pumping and power generation for the pumped-storage station; Emin,Emax denote the minimum and maximum stored energy of the pumped-storage station; Pmaxcharge,Pmaxdischarge are the maximum output for both pumping and generating; and xtcharge,xtdischarge are the pumping status of the pumped-storage station at time period t, and they are Boolean variables.

Unit Output Proportion Constraint Equation 8:

The output of all types of units needs to be allocated to each provincial grid according to the agreed proportions, and during the pumping period, the corresponding proportion of grid power is also consumed accordingly.

Pg,p,t=pprovincesPg,p,t×Rg,p(8)

In the formula, Pg,p,t represents the output power of the gth unit to province p at time period t; Rg,p denotes the power transmission ratio of the gth unit to the province p.

2.1.2 Provincial-level dispatch model

The provincial two-stage robust optimization model presented in this paper is based on the unit dispatch data obtained from the regional grid-level dispatch strategy after load peak shaving. It addresses the source-load uncertainty problem under the integration of renewable energy.

Controllable Generation Units:

Controllable generation units include adjustable gas turbines, diesel generators, and others. After linearization, their cost function can be expressed as a linear function of their generation power Equation 9:

Cg=tTaPg,t+bΔt(9)

In the formula, Pg,t represents the generation power of the unit at time period t; a, b are the cost coefficients; and Δt denotes the duration of each time period.

Output Power Constraint Equation 10:

PgminPg,tPgmax(10)

In the formula, Pgmin,Pgmax denote the minimum and maximum generation power of the unit, respectively.

Energy Storage Components.

The cost of energy storage components is the cost associated with the charging and discharging processes, which can be expressed as Equation 11:

Cs=tTKsPtdisη+PtchηΔt(11)

In the formula, Ptdis,Ptch represent the charging and discharging power of the energy storage at time period t; η is the charging and discharging efficiency; and Ks denotes the unit charging and discharging cost.

The constraints that the energy storage components must satisfy are as follows Equation 12:

0PtdisutsPmax0Ptch1utsPmaxηtTPtchΔttTPtdisΔtη=0EminE0+ηitPichΔtitPidisΔtηEmax(12)

In the formula, Pmax represents the maximum allowable output of the energy storage; uts is a Boolean variable indicating the charging or discharging state of the energy storage; and Emin,Emax denote the lower and upper bounds of the energy storage capacity.

Demand Response Load:

The flexible scheduling process of demand response load is considered in Zhang et al. (2025a), Zhang et al. (2025b), Wang et al. (2025); Li et al. (2025); Shao et al. (2024). Under the condition that the electricity usage characteristics meet the requirements for providing demand response services, the grid can adjust users’ electricity consumption plans while providing appropriate compensation to the users. Thus, the adjustment cost of the demand response load is expressed as Equation 13:

CDR=tTKDRPtDRPtDR0Δt(13)

The constraints that need to be satisfied are as follows Equation 14:

tTPtDRΔt=DDRDtminPtDRΔtDtmax(14)

In the formula, PtDR represents the actual dispatch power of the demand response load at time period t; DDR is the total electricity demand across all time periods; Dtmin,Dtmax denote the minimum and maximum demand response load at time period t, respectively.

Since Equation 13 is a nonlinear function, auxiliary variables are introduced to linearize it. The final cost function and constraints are expressed as Equations 15, 16:

CDR=tTKDRPtDR1+PtDR2Δt(15)
PtDRPtDR0+PtDR1PtDR2=0PtDR10PtDR20(16)

In the formula, PtDR1,PtDR2 are the linearized auxiliary variables.

Power Trading Interaction:

When the generation output cannot meet the actual load demand, power needs to be purchased from other sectors, and when there is excess power, it can be sold to others.

The power trading process needs to satisfy the basic power balance Equation 17:

PtbuyPtsell=Pg,t+PtchPtdis+PtDR+PtLoadPtPV(17)

In the formula, Ptbuy,Ptsell represent the power purchased or sold at time period t; PtLoad,Ptpv denote the load data and photovoltaic generation data at time period t, respectively.

The cost function and constraints of power trading are expressed as follows Equations 18, 19:

Cm=tTptPtbuyPtsellΔt(18)
0PtbuyutmPmax0Ptsell1utmPmax(19)

In the formula, pt represents the day-ahead transaction price; utm is a Boolean variable indicating the purchase or sale state; and Pmax is the maximum exchange power.

In summary, the objective function of this model is to minimize the dispatch cost, expressed as Equation 20:

minobj=Cg+Cs+CDR+Cm(20)

Uncertain Parameters:

To address the conservativeness of robust optimization results, uncertain parameter variables and the maximum values of these parameters are incorporated into the model to limit the frequency of photovoltaic and load data reaching the worst-case scenarios, ensuring that the overall model result is not overly conservative Equation 21:

Btpv,BtLoad(21)

The uncertain parameter is a Boolean variable indicating whether the worst-case scenario is reached, used to control the actual values of photovoltaic and load data in the optimization model. Typically, the worst-case scenario occurs when photovoltaic generation is insufficient while load demand is high. Thus, the model considers the lower bound deviation of photovoltaic generation and the upper bound deviation of load demand.

Ptpv=Ptpv0BtpvΔPtpvPtLoad=PtLoad0+BtpvΔPtLoad(22)

In the formula Equation 22, Ptpv0,PtLoad0 represent the predicted values of photovoltaic and load data, respectively, while ΔPtpv,ΔPtLoad denote the deviation values of photovoltaic and load power at time period t obtained through fuzzy clustering.

Additionally, the uncertain parameters need to satisfy boundary constraints Equation 23:

BtpvBmaxpvBtLoadBmaxLoad(23)

In the formula, Bmaxpv,BmaxLoad denote the upper limit of the uncertain parameter. The larger the upper limit, the more conservative the model results will be.

2.2 Solution method

2.2.1 Regional grid-level dispatch objective: Minimization of unmet load fluctuation

Based on the aforementioned variables and constraints, the optimization model is established. The primary objective at the regional grid-level is to use the output of generation units and pumped-storage power stations to smooth the unmet load curves of each province, thereby reducing the impact of peak loads on provincial-level dispatch. To achieve this goal, referring to the generation dispatch function of the Chinese power system, the minimization of the variance of unmet loads is selected as the objective function (Meng et al., 2023; Wang et al., 2023).

Thus, the objective function of this optimization problem is Equation 24:

minobjp=tTPp,tiTPp,i/T2(24)

Since there is a corresponding load variance for each province, this model is a multi-objective optimization problem. The weighted sum method is used to combine multiple objectives into an equivalent scalar objective.

To avoid unreasonable results caused by large differences in load magnitudes among provincial grids, the total accepted generation output is normalized. Then, the objective functions of different provinces are integrated into a single objective using appropriate weighting coefficients.

Therefore, the composite objective function can be expressed as:

minobj=pprovincesweightp×objp(25)

In the formula Equation 25, objp represents the normalized form of the objective function for province p.

objp=tTPp,tiTPp,i/TmaxLoadp2(26)

In the formula Equation 26, weightp is the weight of province p in the objective function, and maxLoadp is the maximum value in the load curve of province p.

By substituting the sum of generation outputs into Equation 25, we obtain Equation 27:

minobj=pprovincesweightpmaxLoadp2×tTgGPg,p,ttTgGPg,p,tT2(27)

Considering the separate optimization of peak and valley period variances, the two optimization objectives are weighted to transform the multi-objective problem into a single-objective problem again, as expressed by Equation 28:

minobj=weightpeak×objpeak+weightvalley×objvalley(28)

In the formula, weightpeak,weightvalley are the weights for peak and valley periods set in this study, and objpeak,objvalley correspond to the objective functions for peak and valley periods, respectively.

Thus, through the above formula transformations, the multi-objective problem is converted into a single-objective problem. In this model, the objective function is a quadratic function of the decision variables.

The power balance constraints in this model cannot be satisfied. However, research tests show that merely removing the power balance constraints leads to passive output from conventional generation units and non-unique solutions. The primary reason is that the relaxation of constraints is too excessive, requiring additional constraints or changes to the objective function.

This study adopts the penalty function method, adding penalty terms to the objective function in Equation 28, as shown in Equation 29:

ϕ=pprovincestTmax0,Loadp,tgGPgmax×Rg,p×Loadp,tgGPg,p,t(29)

The construction of the penalty function is derived from the weighting method. In the formula, the max function is used to determine whether the maximum generation output of units at time period ttt can meet the load demand. If it cannot be met, unmet load is generated. The difference is used as a weight and multiplied by the difference between the actual unit output and the load. By adding this penalty term to the original objective function in Equation 28, the issue of passive output from conventional generation units can be resolved.

The objective function after adding the penalty term is shown in Equation 30:

minobj=weightpeak×objpeak+weightvalley×objvalley+ϕ(30)

where:

Pp,t=gGPg,p,t(31)

In the formula Equation 31, Pp,t denotes the total received generation output of grid p at time period t; Pg,p,t is the power transmission of unit g to grid p at time period t; and T is the total number of time periods.

The highest degree term of the objective function in this model is quadratic. Therefore, this optimization problem is a quadratic programming problem. To prove that the model can converge to the global optimum, the compact form of the model is presented as follows Equation 32:

minxTax+bxs.t.Keqx=kDxd(32)

In the compact form, Keq,D represent the constraint coefficient matrix in the model, and k, d are the constant matrixes corresponding to the constraints. Since the coefficient matrix a of the quadratic term in the objective function is a symmetric matrix, its eigenvalues are real numbers. By solving the eigenvalue equation, It can be concluded that all eigenvalues are non-negative, indicating that a is a positive semi-definite matrix. Thus, this problem is a convex quadratic optimization problem. For convex optimization, under the conditions that are satisfied, any local optimal solution is a global optimal solution, and this problem is solvable within polynomial time.

2.2.2 Provincial-level dispatch objective: minimization of economic cost

The provincial two-stage robust dispatch model constructed in this paper is mainly solved in two stages: the master problem and the sub-problem.

The compact form of the provincial-level dispatch problem is as follows:

minxmaxuUminyFx,ucys.t.GyhExMuGeqy=heqEeqxMequ(33)

In the formula Equation 33, the variables represented by x Equation 34, are:

x=uts,utmT(34)

The variables represented by y Equation 35 are:

y=Pg,t,Ptch,Ptdis,PtDR,PtDR1,PtDR2,Ptbuy,Ptsell,Ptpv,PtLoadT(35)
Equation 36 indicates the frequency of reaching the worst-case scenarios.
u=Btpv,BtLoadT(36)

G,h,E,M and Geq,heq,Eeq,Meq respectively denote the coefficient matrices corresponding to inequality and equality constraints.

For the objective function, it is divided into two layers, the inner and the outer. The outer layer minx* represents finding the minimum value of the inner function and making decisions on the variables x, so that x can be used as a constant term in the inner layer’s solution. The inner function maxuUminyFx,ucy represents the cost minimization result under the worst-case scenario. The feasible region yFx,u is defined by the feasible values of variables x and uncertain parameters u. The function minyFx,ucy represents the optimization result within this feasible region, while maxuU* identifies the scenario with the worst-case photovoltaic and load values among all cost minimization results.

Based on the above analysis, the master problem and sub-problem can be organized as follows:

The master problem focuses on making decisions for the state variables of energy storage and power trading in the day-ahead stage, using the predicted photovoltaic and load data to optimize and obtain the lower bound of the final solution, as shown in Equation 37:

minas.t.acyGyhMuExGeqy=heqMequEeqx(37)

The sub-problem is solved in the intraday stage, where, based on the already determined values of x, the goal is to find the cost-minimizing decision under the worst-case scenario, as shown in Equation 38.

max min cys.t.GyhMuExvalueGeqy=heqMequEeqxvalue(38)

For the two-stage robust optimization model above, this study uses the Column-and- Constraint Generation (C&CG) algorithm to solve it. The C&CG algorithm achieves the optimal solution of the original problem by decomposing it into a master problem and a sub-problem and solving them iteratively in an alternating manner. The main solution logic is illustrated in Figure 2.

Figure 2
Flowchart outlining a computational procedure. It starts with initializing variables LB and UB. The master problem is solved, and the scheduling result updates LB. The optimal scheduling result is fixed to solve the subproblem. The target value and scenario probability update UB. If the difference between LB and UB is less than a threshold, the final results are outputted. Otherwise, variables and constraints are updated, and the process repeats. The procedure ends once the condition is satisfied.

Figure 2. Flowchart of the C&CG algorithm.

In the process of solving the master problem using the C&CG algorithm, variables and constraints related to the sub-problem are continuously introduced, allowing for a tighter lower bound of the original objective function value, thereby effectively reducing the number of iterations.

The master problem and sub-problem of the provincial-level dispatch are given by Equations 37 and 38 above. During the iterative process of the C&CG algorithm, the uncertain parameters and their corresponding constraints are dynamically updated. Thus, the formal representations of the master problem and sub-problem should be expressed as follows:

Master Problem Equation 39:

minas.t.acyGyhMukExGeqy=heqMequkEeqx(39)

Subproblem Equation 40:

max min cys.t.GyhMukExvalueGeqy=heqMequkEeqxvalue(40)

fundinWhere uk denotes the value of the uncertain parameter selected in the kth iteration. In the sub-problem, the variables x determined by the master problem are treated as constants.

Since the sub-problem is a bilevel problem, this study employs the KKT algorithm to reformulate the inner cost minimization problem into KKT conditions, transforming the bilevel problem into a single-level problem, as shown in Equation 41:

maxcys.t.GyhMukExvalueGeqy=heqMequkEeqxvalueGTπ1+GeqTπ2=chMukExvalueGyπ1=0Geqyheq+Eeqxvalue+Mequπ2=0cGTπ1GeqTπ2y=0y0π10(41)

It should be noted that the problem contains nonlinear constraints. By using the Big-M method to linearize the constraints, the sub-problem is ultimately converted into a single-level optimization problem, as shown in Equation 42:

maxcys.t.GyhMukExvalueGeqy=heqMequkEeqxvalueGTπ1+GeqTπ2=cπ1MvGyh+Exvalue+MuM1vπ2MlGeqyheq+Eeqxvalue+MequM1lyMwcGTπ1Geqπ2M1wy0π10(42)

In the formula, v,l,w represent auxiliary operators in the Big-M method, all of which are Boolean variables; M is a constant significantly larger than all parameters; π1,π2 denote the Lagrange multipliers in the KKT transformation.

3 Results and discussion

3.1 Regional grid-level dispatch: minimizing the fluctuation of unmet load in each province

Given that the focus of this chapter is on multi-province load peak shaving, a dispatch model for multiple provinces is constructed.

In the inter-provincial dispatch model, four conventional generation units and two pumped-storage units supplying power to two provinces are considered, using 24-period load demand data. The weights for receiving generation output among different provinces are determined based on the ratio of the average load values of each province, and weighting for peak and valley periods is also taken into account.

The average load values of each province and the corresponding weights, generation allocation ratios, and unit data for each generation unit are shown in Tables 13.

Table 1
www.frontiersin.org

Table 1. Average load values and weights of each province.

Table 2
www.frontiersin.org

Table 2. Generation unit allocation ratios.

Table 3
www.frontiersin.org

Table 3. Generation unit data.

From the data in Table 1, it can be observed that the total maximum output of all units in this study is less than the average load of each province, indicating that the load demand cannot be met in every time period.

Based on the data from Tables 1 and 3, a model is constructed, and the results are shown in Figures 38.

Figure 3
Chart (a) displays a bar graph of power fluctuations over a time period, showing both positive and negative values. Chart (b) shows a bar graph of energy levels over the same period, generally increasing with some fluctuations. Both charts are labeled

Figure 3. Pumped storage 1 status. (a) Pumped storage 1 Power. (b) Pumped storage 1 Energy.

Figure 4
Panel (a) is a bar chart showing

Figure 4. Pumped storage 2 status. (a) Pumped storage 2 Power. (b) Pumped storage 2 Energy.

Figure 5
Line graph titled

Figure 5. Peak shaving results of unmet load in province 1.

Figure 6
Line graph titled

Figure 6. Peak shaving results of unmet load in province 2.

Figure 7
Bar chart showing energy levels consistently at 210 across time periods from 0 to 24. The legend indicates

Figure 7. Total power received from conventional units in province 1.

Figure 8
Bar chart titled

Figure 8. Total power received from conventional units in province 2.

As shown in Figures 7, 8, the introduction of the penalty function allows conventional generation units to maintain a high power state even without power balance constraints, indicating that the penalty function construction in this study is effective. Additionally, as seen in Figures 5, 6, the fluctuation degree of the unmet load stabilizes after peak shaving. For Province 1, the valley and peak periods stabilize at 100 MW and 200 MW, respectively, while for Province 2, they stabilize at 300 MW and 440 MW, respectively. This is beneficial for determining the output of backup generation units or the power output from electricity sellers. Stable power delivery can reduce the probability of power system incidents and increase economic benefits.

3.2 Provincial dispatch: robust optimization for uncertainty

The grid operation parameters are shown in Table 4.

Table 4
www.frontiersin.org

Table 4. Grid operation parameters.

After 24 iterations, the model converged, and the upper and lower bound variation curves are shown in Figure 9.

Figure 9
Line chart titled

Figure 9. Upper and lower bound variation curves.

The final results are as follows: the upper bound value is approximately 5,068.342, and the lower bound value is approximately 5,057.569. The result using a standard economic dispatch program is 5,055.049, indicating that the model incurs additional costs to account for extreme scenarios while not deviating significantly from the results of general dispatch methods. This suggests that the model’s solution is reliable.

The load demand curve and photovoltaic output curve are shown in Figures 10, 11.

Figure 10
Line graph titled

Figure 10. Load demand curve.

Figure 11
Line graph titled

Figure 11. Photovoltaic output curve.

The power trading curve and day-ahead transaction price curve are shown in Figures 12, 13.

Figure 12
Bar chart titled

Figure 12. Power trading curve.

Figure 13
Bar chart showing day-ahead trading electricity prices over time periods. Prices range from zero to over two hundred. Peaks occur around periods ten and eighteen with prices reaching approximately one hundred seventy-five.

Figure 13. Day-ahead transaction price curve.

The controllable generation unit output curve and energy storage output curve are presented in Figures 14, 15.

Figure 14
Bar chart titled

Figure 14. Controllable generation unit output curve.

Figure 15
Bar chart titled

Figure 15. Energy storage output curve.

The expected demand response plan and actual demand response load are shown in Figures 16, 17.

Figure 16
Bar chart titled

Figure 16. Expected demand response plan.

Figure 17
Bar chart titled

Figure 17. Demand response load curve.

As shown in Figure 11, during periods 1 to 4 and 19 to 24, photovoltaic output is zero, indicating that the load demand is entirely met by controllable generation units, energy storage components, and power purchases. During these periods, when the day-ahead transaction price is lower than the generation cost of controllable units, the output power of controllable units is reduced to the minimum value, as shown in periods 0 to 6 in Figures 12, 14. In the remaining periods, the output power of controllable units is increased to the maximum value to increase the power available for external sales, as seen in periods 7 to 16 and period 20 in Figure 12, or to reduce the power purchase quantity, as shown in periods 17 to 19 and periods 21 and 22 in Figure 12, thereby reducing operational costs.

As observed in Figure 15, under the given pricing mechanism, the energy storage components are charged during periods 6 and 23 and discharged during periods 8, 19, and 20, thereby storing energy during low-price periods and selling it during high-price periods. In Figure 16, the peak electricity price period corresponds to the peak of the expected demand response load plan. By redistributing the electricity demand from periods 18 to 22 to periods 1 to 6 and period 24, the total electricity demand and time-specific electricity constraints are satisfied, reducing the amount of energy that needs to be purchased during peak price periods.

4 Conclusion

This paper primarily investigates the collaborative optimization strategy for regional grid-level and provincial-level dispatch under the influence of renewable energy, focusing on two main issues: peak load shaving in regional grid-level dispatch and uncertainty optimization in provincial-level dispatch.

At the regional grid-level dispatch stage, this paper delves into strategies for effective management of peak shaving and unmet load in a multi-province grid system using pumped-storage units. Traditional peak shaving and valley filling strategies primarily adjust grid operating modes to balance power supply and demand, reducing power consumption during peak periods and increasing power reserves during valley periods. By introducing power variance as the objective function of the optimization problem, this study not only aims to smooth the traditional load curve but also addresses the management of unmet load. Through analyzing the actual power dispatch scenarios between the East China grid and other provincial grids, the paper identifies potential issues of insufficient unit output in cross-regional power trading. To tackle this problem, the study proposes further exploration into the peak shaving demand caused by unmet load and how to manage it through optimized dispatch strategies. To effectively manage load and its residual components, a comprehensive planning method considering both load peak shaving and multi-province dispatch is proposed. By flexibly scheduling energy storage systems such as pumped-storage units, effective management of both load and unmet load can be achieved. This research provides a new perspective and approach for power system operation, enhancing the overall stability and reliability of the grid by addressing unmet load management in addition to traditional peak shaving and valley filling.

At the provincial-level dispatch stage, a two-stage robust optimization model is proposed to minimize economic costs. This model aims to enhance the operational efficiency and reliability of the power system by closely integrating day-ahead dispatch and real-time dispatch. Based on existing modeling methods considering source-load uncertainty, a two-stage robust optimization model is developed. The first stage, day-ahead dispatch, primarily determines the unit status for the following day, while the second stage, real-time dispatch, deals with uncertainties caused by forecast deviations. The model effectively decomposes the complex optimization problem into a master problem and sub-problems. Using the KKT condition algorithm and Column-and-Constraint Generation (C&CG) algorithm, the model can be solved, improving decision accuracy and execution speed. The effectiveness of the proposed two-stage robust optimization model is validated through specific case studies. The results show that this model can maintain power system stability and efficiency by appropriately increasing costs in the face of high uncertainty.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

BS: Funding acquisition, Writing – original draft, Conceptualization. KQ: Methodology, Funding acquisition, Writing – review and editing. MW: Writing – review and editing, Conceptualization, Funding acquisition. MZ: Validation, Writing – original draft. KZ: Validation, Writing – review and editing, Writing – original draft. GH: Methodology, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Science and Technology Project of East China Branch of State Grid under Grant No. 529924240008. We would like to extend our sincere gratitude to all project collaborators and contributors for their invaluable insights and support throughout the research.

Conflict of interest

Authors BS, KQ, and MW were employed by East China Branch of State Grid Corporation of China.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The authors declare that this study received funding from East China Branch of State Grid. The funder had the following involvement in the study: providing access to power grid operational data and technical expertise.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Bendotti, P., Chrétienne, P., Fouilhoux, P., and Pass-Lanneau, A. (2023). The anchor-robust project scheduling problem. Operations Res. 71 (6), 2267–2290. doi:10.1287/opre.2022.2315

CrossRef Full Text | Google Scholar

Büsing, C., and Schmitz, S. (2024). Robust two-stage combinatorial optimization problems under discrete demand uncertainties and consistent selection constraints. Discrete Appl. Math. 347, 187–213. doi:10.1016/j.dam.2023.12.028

CrossRef Full Text | Google Scholar

Cheng, C., Su, C., Wang, P., Shen, J., Lu, J., and Wu, X. (2018). An MILP-based model for short-term peak shaving operation of pumped-storage hydropower plants serving multiple power grids. Energy 163, 722–733. doi:10.1016/j.energy.2018.08.077

CrossRef Full Text | Google Scholar

Chua, K. H., Lim, Y. S., and Morris, S. (2016). Energy storage system for peak shaving. Int. J. Energy Sect. Manag. 10 (1), 3–18. doi:10.1108/ijesm-01-2015-0003

CrossRef Full Text | Google Scholar

Dong, Z., Zhang, Z., Huang, M., Yang, S., Zhu, J., Zhang, M., et al. (2024). Research on day-ahead optimal dispatching of virtual power plants considering the coordinated operation of diverse flexible loads and new energy. Energy 297, 131235. doi:10.1016/j.energy.2024.131235

CrossRef Full Text | Google Scholar

Jin, X., Liu, B., Liao, S., Cheng, C., Zhang, Y., Zhao, Z., et al. (2022). Wasserstein metric-based two-stage distributionally robust optimization model for optimal daily peak shaving dispatch of cascade hydro-plants under renewable energy uncertainties. Energy 260, 125107. doi:10.1016/j.energy.2022.125107

CrossRef Full Text | Google Scholar

Kong, F., Mi, J., and Wang, Y. (2022). A two-stage distributionally robust optimization model for optimizing water-hydrogen complementary operation under multiple uncertainties. J. Clean. Prod. 378, 134538. doi:10.1016/j.jclepro.2022.134538

CrossRef Full Text | Google Scholar

Li, L., Fan, S., Xiao, J., Zhang, Y., Huang, R., and He, G. (2025). Energy management strategy for community prosumers aggregated VPP participation in the ancillary services market based on P2P trading. Appl. Energy 384, 125472. doi:10.1016/j.apenergy.2025.125472

CrossRef Full Text | Google Scholar

Li, L., Fan, S., Xiao, J., Zhou, H., Shen, Y., and He, G. (2024). Fair trading strategy in multi-energy systems considering design optimization and demand response based on consumer psychology. Energy 306, 132393. doi:10.1016/j.energy.2024.132393

CrossRef Full Text | Google Scholar

Liao, S., Xiong, J., Liu, B., Cheng, C., Zhou, B., and Wu, Y. (2024). MILP model for short-term peak shaving of multi-grids using cascade hydropower considering discrete HVDC constraints. Renew. Energy 235, 121341. doi:10.1016/j.renene.2024.121341

CrossRef Full Text | Google Scholar

Liu, Y., and Peng, M. (2024). Research on peak load shifting for hybrid energy system with wind power and energy storage based on situation awareness. J. Energy Storage 82, 110472. doi:10.1016/j.est.2024.110472

CrossRef Full Text | Google Scholar

Luo, B., Liu, X., Chen, Y., Zhou, C., and Long, X. (2024). Short-term peak shaving model of cascade hybrid pumped storage hydropower station retrofitted from conventional hydropower. Int. J. Electr. Power Energy Syst. 162, 110226. doi:10.1016/j.ijepes.2024.110226

CrossRef Full Text | Google Scholar

Meng, Y., Fan, S., Shen, Y., Xiao, J., He, G., and Li, Z. (2023). Transmission and distribution network-constrained large-scale demand response based on locational customer directrix load for accommodating renewable energy. Appl. Energy 350, 121681. doi:10.1016/j.apenergy.2023.121681

CrossRef Full Text | Google Scholar

Niu, T., Yin, H., and Feng, E. (2022). An interval two-stage robust stochastic programming approach for steam power systems design and operation optimization under complex uncertainties. Chem. Eng. Sci. 253, 117533. doi:10.1016/j.ces.2022.117533

CrossRef Full Text | Google Scholar

Shao, Y., Fan, S., Meng, Y., Jia, K., and He, G. (2024). Personalized demand response based on sub-CDL considering energy consumption characteristics of customers. Appl. Energy 374, 123964. doi:10.1016/j.apenergy.2024.123964

CrossRef Full Text | Google Scholar

Uddin, M., Romlie, M. F., Abdullah, M. F., Abd Halim, S., Kwang, T. C., and Chia Kwang, T. (2018). A review on peak load shaving strategies. Renew. Sustain. Energy Rev. 82, 3323–3332. doi:10.1016/j.rser.2017.10.056

CrossRef Full Text | Google Scholar

Wallberg, A., Castellucci, V., Flygare, C., Lind, E., Schultz, E., Mattos, M. M., et al. (2024). Negative correlation peak shaving control in a parking garage in Uppsala, Sweden. Appl. Energy 375, 124082. doi:10.1016/j.apenergy.2024.124082

CrossRef Full Text | Google Scholar

Wang, J., Xu, X., Li, H., and Chen, H. (2022). Two-stage robust optimization of thermal-ESS units scheduling under wind uncertainty. Energy Rep. 8, 1147–1155. doi:10.1016/j.egyr.2022.02.277

CrossRef Full Text | Google Scholar

Wang, J., Zhang, S., Huo, J., Zhou, Y., Li, L., and Han, T. (2021). Dispatch optimization of thermal power unit flexibility transformation under the deep peak shaving demand based on invasive weed optimization. J. Clean. Prod. 315, 128047. doi:10.1016/j.jclepro.2021.128047

CrossRef Full Text | Google Scholar

Wang, M., Xu, X., Fan, S., Yan, Z., Yang, B., and Guan, X. (2025). Multi-period electrical array reconfiguration for multi-PV systems considering net power fluctuation mitigation. J. Mod. Power Syst. Clean Energy 203, 68–80. doi:10.35833/MPCE.2024.000779

CrossRef Full Text | Google Scholar

Wang, M., Xu, X., and Yan, Z. (2023). Online fault diagnosis of PV array considering label errors based on distributionally robust logistic regression. Renew. Energy 203, 68–80. doi:10.1016/j.renene.2022.11.126

CrossRef Full Text | Google Scholar

Yang, S., Fang, J., Zhang, Z., Lv, S., Lin, H., and Ju, L. (2024). Two-stage coordinated optimal dispatching model and benefit allocation strategy for rural new energy microgrid. Energy 292, 130274. doi:10.1016/j.energy.2024.130274

CrossRef Full Text | Google Scholar

Yang, Z., Zheng, H., Du, Y., Guo, L., Xiong, X., and Yao, G. (2023). “Two-stage distributed robust optimization based on multi-discrete scenarios,” in 2023 10th International Forum on Electrical Engineering and Automation (IFEEA), Nanjing, China, 1275–1280.

Google Scholar

Zeng, B., and Zhao, L. (2013). Solving two-stage robust optimization problems using a column-and-constraint generation method. Operations Res. Lett. 41 (5), 457–461. doi:10.1016/j.orl.2013.05.003

CrossRef Full Text | Google Scholar

Zhang, Y., Fan, S., Meng, Y., and He, G. (2025a). Payment and incentive allocation in demand response based on cost causation principle. IEEE Trans. Industry Appl. 1–14. doi:10.1109/tia.2025.3575746

CrossRef Full Text | Google Scholar

Zhang, Y., Meng, Y., Fan, S., Xiao, J., Li, L., and He, G. (2025b). Multi-time scale customer directrix load-based demand response under renewable energy and customer uncertainties. Appl. Energy 383, 125334. doi:10.1016/j.apenergy.2025.125334

CrossRef Full Text | Google Scholar

Zhao, C., and Guan, Y. (2015). Data-driven stochastic unit commitment for integrating wind generation. IEEE Trans. Power Syst. 31 (4), 2587–2596. doi:10.1109/tpwrs.2015.2477311

CrossRef Full Text | Google Scholar

Zhao, H., Liao, S., Ma, X., Fang, Z., Cheng, C., and Zhang, Z. (2024). Short-term peak-shaving scheduling of a hydropower-dominated hydro-wind-solar photovoltaic hybrid system considering a shared multi-energy coupling transmission channel. Appl. Energy 372, 123786. doi:10.1016/j.apenergy.2024.123786

CrossRef Full Text | Google Scholar

Zhu, X., Guo, Y., Li, A., Li, S., Zhang, J., Gu, B., et al. (2024). Two-stage robust optimization of unified power quality conditioner (UPQC) siting and sizing in active distribution networks considering uncertainty of loads and renewable generators. Renew. Energy 224, 120197. doi:10.1016/j.renene.2024.120197

CrossRef Full Text | Google Scholar

Glossary

Sets

T Set of time periods, indexed by t

P Set of provinces, indexed by p

G Set of generation units, indexed by g

Variables and parameters

Regional grid-level variables

Pg,p,t Output power of unit g to province p at time period t

Et Equivalent stored energy of pumped-storage unit at time period t

Ptcharge Pumping power of pumped-storage unit at time period t

Ptdischarge Generating power of pumped-storage unit at time period t

xtcharge Boolean variable for pumping status at time period t

xtdischarge Boolean variable for generating status at time period t

Provincial-level variables

Pg,t Generation power of controllable unit g at time period t

Ptch Charging power of energy storage at time period t

Ptdis Discharging power of energy storage at time period t

uts Boolean variable for energy storage charging/discharging state

PtDR Actual dispatch power of demand response load at time period t

Ptbuy Power purchased at time period t

Ptsell Power sold at time period t

utm Boolean variable for power trading state (purchase/sale)

Ptpv Actual photovoltaic power at time period t

PtLoad Actual load power at time period t

Btpv Boolean uncertain parameter for photovoltaic power

BtLoad Boolean uncertain parameter for load power

Key parameters

RUg Ramp rate limit of thermal power unit g

Pgmax Maximum output of unit g

Pgmin Minimum generation power of unit g

ηcharge Efficiency of pumping for pumped-storage station

ηdischarge Efficiency of power generation for pumped-storage station

Emin Minimum stored energy of pumped-storage station

Emax Maximum stored energy of pumped-storage station

Rg,p Power transmission ratio of unit g to province p

a,b Cost coefficients for controllable generation units

Ks Unit charging and discharging cost of energy storage

KDR Cost coefficient for demand response load

η Charging and discharging efficiency of energy storage

DDR Total electricity demand across all time periods

pt Day-ahead transaction price at time period t

Keywords: uncertainty, coordinated scheduling, unmet load, peak shaving, robust optimization, C&CG

Citation: Song B, Qin K, Wen M, Zhu M, Zou K and He G (2025) Coordinated multi-level scheduling method considering uncertainty of renewable energy and load. Front. Energy Res. 13:1612065. doi: 10.3389/fenrg.2025.1612065

Received: 15 April 2025; Accepted: 04 July 2025;
Published: 18 August 2025.

Edited by:

Morteza Nazari-Heris, East Carolina University, United States

Reviewed by:

Haitham A. Mahmoud, King Saud University, Saudi Arabia
Chengwei Lou, China Agricultural University, China

Copyright © 2025 Song, Qin, Wen, Zhu, Zou and He. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Moyan Zhu, MTM0MDE4Mzk5MzJAc2p0dS5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.