1. Introduction
In recent years, China has seen a steady rise in its new energy installed capacity. According to the National Energy Administration’s Demand Side Management for Peaking Coal power generation [
1], China is striving to balance a secure and stable energy supply with sustainable, green, and low-carbon growth. This approach is in line with the country’s methodical drive toward achieving carbon peak and carbon neutrality. The energy infrastructure in China has been progressively refined, boasting over 1.05 billion kilowatts of ultra-low emission coal-fired power generation units. Additionally, the share of clean energy consumption has climbed from 20.8% to surpass 25%, reflecting the nation’s commitment to a more environmentally considerate energy mix.
China’s energy resources are marked by a scarcity of oil and gas, with a notable abundance of coal. Wind (WD), photovoltaic (PV), and other clean energy sources are catching up at a relatively fast speed. China’s installed coal power capacity will increase from 1.01 billion kilowatts in 2018 to 1.12 billion kilowatts in 2022, a net increase of only 110 million kilowatts. Renewable energy capacity surged from 728 billion kW in 2018 to 1.213 billion kW in 2022, constituting 47.3% of the total installed capacity. In 2022, China saw a historic moment as the total installed capacity of renewable energy surpassed that of coal power [
1]. In the near future, new energy generation will become the main electricity supply. However, the high proportion of clean energy access for the stability of the power grid is a huge challenge. At the same time, we must also attach great importance to the consumption of new energy. The randomness, volatility, and uncertainty in renewable energy output pose inherent challenges to its consumption [
2]; the system in the new energy large-scale output and intermittent output peak–valley difference is large, and the peak load pressure continues to increase.
However, by the end of 2022, although the proportion of coal power capacity in China’s total installed capacity dropped to about 43.8%, the proportion of electricity generation was still as high as 58.4% [
1]. Currently, China’s power system relies heavily on thermal power in its supply structure. The flexible adaptation of traditional thermal power units for deep peak regulation can greatly enhance the adjustment capacity on the power side [
3]. Existing studies have shown that deep peak shaving is one of the most effective ways to boost renewable energy consumption in the flexible transformation of thermal power units [
4,
5,
6]. At present, there have been many studies on deep peak shaving of thermal power units. Aiming at the phenomenon of power grid frequency fluctuation caused by the increasing proportion of renewable energy, Reference [
7] proves the importance of deep peak regulation and on–off regulation for stabilizing system frequency through a series of derivations. In Reference [
8], taking into account the extra coal consumption loss and unit life reduction of thermal power units operating in deep peak regulation without oil (DPR) mode, a robust optimization scheduling model of scale WD power grid-connected and DPR cost was proposed, and, finally, a robust day-ahead scheduling scheme with optimal economy was obtained. On the basis of thermal power storage, Reference [
9] established a two-tier model involving hierarchical utilization of energy storage and DPR in the assistant service market, which can effectively improve the phenomenon of WD’s abandonment and enhance the competitiveness of fossil power in the assistant service market of peak regulation. In Reference [
10], proposing a unit commitment comprehensive optimal model aimed to minimize total costs by optimizing wind power curtailment in the context of the expensive deep peak regulation of thermal units. However, the above reference did not fully consider the probability distribution information of uncertain factors when optimizing the scheduling of uncertain outputs, such as WD and PV, which led to conservative optimization decision results and a poor economy. In order to address this issue, this paper incorporates the uncertainty of the probability distribution of WD and PV power output through the construction of fuzzy sets in day-ahead scheduling. To mitigate renewable energy abandonment, we leverage peak load balancing advantages, enhance clean energy absorption rates, and lower system operating costs.
Currently, two predominant research methodologies are employed to address the uncertainty inherent in integrating clean energy into the power system. These are stochastic optimization (SO) [
11,
12,
13,
14] and robust optimization (RO) [
15,
16,
17,
18]. Each method provides strategic frameworks for managing the unpredictable nature of clean energy outputs. Generally, SO necessitates a distribution model for uncertain parameters, which involves establishing numerous variables and constraints [
14]. Given that the coefficient probability distributions are known, these constraints can be converted into deterministic ones for resolution. However, this often results in issues, such as an extensive computational scale and reduced efficiency in solving the model. The concept of RO is comparatively conservative, as it does not demand an accurate distribution model for uncertain parameters. Instead, the stochastic variability of the variables is characterized by a specified fluctuation boundary [
19]. If the value of the variable remains within this boundary, an optimal solution can be derived using the robust optimization model [
20]. Distributionally Robust Optimization (DRO), a data-driven approach, amalgamates the strengths of both SO and RO. It offers a novel solution paradigm that addresses the low precision of the SO model and the inherent conservatism of the RO model [
21]. DRO constructs a set of uncertain probability distributions grounded in historical data and employs 1-norm and ∞-norm constraints over scenario probability distribution fuzzy sets. The aim is to ascertain the optimal solution under the premise that the prediction error of the uncertain factors adheres to the worst-case probability distribution [
22]. This methodology adeptly circumvents the issue of nonlinear relationships inherent in indeterminate polynomials and the product of dual variables.
Similarly to the traditional two-stage robust optimization (TRO) solution, DRO also obtains the optimal solution of the problem through the iteration of the main problem and subproblem, but the TRO usually only considers the robustness, and the economy is poor; DRO takes robustness and economy into consideration through the iterative solution of the main subproblem and finally achieves the unity of robustness and economy.
Currently, distributionally robust optimization has found widespread application in the energy sector and integrated energy systems to accommodate renewable energy sources and minimize their impact on the energy grid [
22,
23]. However, its utilization in power systems remains limited. With the integration of large-scale renewable energy sources into power grids, conventional thermal power units are often compelled to shift from their traditional role as primary power suppliers to serving as auxiliary power sources to balance the fluctuating demand of the grid. Consequently, thermal power units may operate at peak loads for extended durations, leading to significant wear and tear on the units. Previous research indicates that energy storage systems can help alleviate the uncertainties associated with renewable energy [
24].
Therefore, to alleviate the peak load on thermal power units and enhance the integration of renewable energy, this paper presents a distributionally robust optimization operation strategy of a WD–PV fire storage power system considering the deep peak shaving of thermal power units. Under the constraint of the comprehensive norm and considering the application of energy storage, to improve the absorption rate of renewable energy, a coordinated operation strategy of system robustness and a model economy is constructed. The main contributions of this paper are as follows:
(1) A fine modeling of a thermal power unit is carried out. By coupling energy storage equipment, the absorption rate of renewable energy is improved, and a distributionally robust optimal scheduling model of a thermal power unit coupled with WD and PV storage is established.
(2) The probability of uncertain scene distribution is constrained by combining 1-norm and ∞-norm, and the probability distribution of the worst scene is determined by updating the scene probability value.
(3) In order to maximize the consumption of new energy and reduce the abandonment of renewable energy, a collaborative optimization strategy of power system robustness and economy is established through comprehensive norm constraints. In addition, through the comparison of a variety of cases, the comprehensive effect of the model is verified.
The remainder of this paper is organized as follows:
Section 2 introduces the peak-load model of the thermal power unit;
Section 3 introduces the structure and operation strategy of the distributionally robust optimization model; and
Section 4 introduces the optimization framework, data, and simulation results of the model.
Section 5 provides a summary of the main discoveries, acknowledges limitations, and proposes future research directions.
2. Peak-Load Model of the Thermal Power Unit
In the power system with multiple thermal power units, especially when there are significant changes in power load or the integration of new energy sources, two typical modes of peak regulation operations emerge: deep peak regulation and on–off peak regulation.
2.1. Depth Peak Shaving of Thermal Power Units
Deep peak shaving in thermal power units involves modifying their output to match fluctuations in the integration of new energy sources. This adjustment aims to align with the output of wind and photovoltaic power generation, ensuring the real-time balance of system power [
25].
Based on the operational status and energy consumption traits of thermal power units, the peak regulation process can be categorized into three modes: regular peak regulation (RPR), deep peak regulation without oil (DPR), and deep peak regulation with oil (DPRO). Relevant modes of depth peak shaving are shown in
Figure 1.
Among them, Pb refers to the minimum operating output of the unit, Pa refers to the limit output of stable combustion of the unit during deep peak regulation, Pmin refers to the limit output of stable combustion of the unit during DPRO, and Pmax refers to the maximum operating output of the unit.
With the integration of large-scale, new energy sources into the power grid, the design considerations for 600 MW and 1000 MW high-capacity thermal power units emphasize the necessity of possessing a certain peak load capacity. Earlier investments in smaller-capacity thermal power units below 300 MW in China were constrained by technological limitations, enabling only base load capacity. At present, these smaller-capacity units have yet to undergo comprehensive transformation, resulting in slower responsiveness during depth peak regulation. Following the transformation of thermal power units, grid companies have significantly enhanced the depth of peak regulation, gradually augmenting the flexibility of these units to engage in deep peak regulation during periods of high demand for grid regulation [
26]. However, the deep peak shaving operation of modified coal-fired power units deviates from their rated optimal operation, leading to performance deterioration in subsystems and auxiliary equipment [
27]. Consequently, this gives rise to issues, such as a significant reduction in unit lifespan, elevated coal consumption costs, and increased carbon emissions within the power supply. The costs associated with peak regulation for thermal power units can be broadly categorized into coal consumption cost, shaft life cost, and oil injection cost.
- (1)
Coal consumption cost of the thermal power unit.
The coal consumption cost of thermal power units is usually expressed by consumption characteristics. A fixed set of abc coefficients [
28] is usually adopted for specific units participating in generation optimization scheduling. The coal consumption characteristic parameters of thermal power units can be obtained through function fitting. The coal consumption cost of units participating in peak load balancing can be shown in Equation (1):
where
represents the coal consumption cost of the thermal power unit
at
moment,
represents the actual output power of the thermal power unit,
represents the price (RMB/t) of coal purchased from the thermal power plant, and
,
,
represent the characteristic parameters for coal consumption in the thermal power unit
.
- (2)
Life cost of the rotating shaft
The steam turbine stands as the primary power generation equipment within a thermal power station, serving as the core component in most thermal power units. Its operation encompasses variable and fixed operating conditions. During transitions like start-ups and shut-downs, significant changes occur in steam temperature, pressure, and other parameters, leading to uneven heating of the rotor and primarily causing low-cycle fatigue damage. In contrast, when the steam turbine operates under fixed conditions, the metallic material comprising the rotor tends to gradually deteriorate due to thermal stress in high-temperature environments, leading to damage primarily induced by high-temperature creep [
29]. Throughout operation, both low-cycle fatigue damage and high-temperature creep damage typically manifest simultaneously. Coupled with varying capacities across different thermal power units, the life loss of the rotor shaft is exacerbated. This paper calculates the rotor shaft’s lifespan cost using the Manson–Coffin formula as a reference, expressing its life cost through Equations (2) and (3):
where
denotes the thermal power unit
in the depth peak load, the life cost of the rotating shaft at
time,
represents the function of the cracking cycle of the rotor on the actual output power [
6], and
represents the construction cost (RMB/MW) of the thermal power unit.
- (3)
Cost of fuel injection
During the oil injection deep peak shaving phase of a thermal power unit, the combustion of coal within the unit experiences instability. Injecting oil into the boiler becomes necessary to facilitate static coal burning, thereby ensuring the stable operation of both the boiler and the unit’s water cycle. The oil injection cost for the thermal power unit can be expressed using Formula (4):
In Formula (4), represents the oil cost of the thermal power unit in the depth peak load at time, represents the fuel consumption of the thermal power unit in the DPRO stage, and represents the price of unit oil.
In conclusion, the cost of deep peak shaving for thermal power units can be represented using a piecewise function (5):
- (4)
Analysis of carbon emissions in the peak regulation process
The calculation of carbon emissions of thermal power units in the process of peak load balancing is the basis of realizing low-carbon scheduling of a power system. At present, the carbon emission calculation methods for power system scheduling mainly include three methods.
It is approximately considered that the carbon emission of thermal power units is proportional to its output, and the proportional coefficient is carbon emission intensity [
30]. This method is simple and clear, and it has been adopted in most references.
Polynomial functions similar to unit consumption characteristics are used for fitting, such as quadratic function [
31], cubic function [
32], etc.
The power supply unit’s coal consumption is multiplied by the carbon dioxide emission coefficient [
33].
In this study, the quadratic function in the second method is employed to model carbon emissions in the peak shaving process, so the process can be expressed by Equations (6) and (7):
where
,
, and
, respectively, represent the total carbon emissions of the power system, the carbon emissions of thermal power units in the peak load balancing process, and the equivalent carbon emissions of the system to the power grid,
,
and
, respectively, represent the carbon emission characteristic function parameters of unit
, and
denotes the equivalent carbon emission characteristic function parameters of power purchase.
In the actual peak shaving process, the involved devices include electric energy storage and WD–PV power generation. Given that wind power and photovoltaic systems operate on clean energy sources, their carbon emissions are not factored into the actual power generation process. Additionally, the charge and discharge process of energy storage primarily involves internal chemical reactions devoid of carbon emissions. Hence, this paper disregards the consideration of carbon emissions from electric energy storage, wind power, and photovoltaic sources.
2.2. Thermal Power Unit On–Off Peak Regulation
On–off regulation in thermal power units refers to a process where, in the event of a significant increase or decrease in renewable energy penetration or a heightened load peak–valley difference within the power system, merely employing depth peak load regulation becomes insufficient for the thermal power unit to maintain the system’s power balance. Consequently, expanding the peak regulation scope necessitates the adjustment of the number of connected thermal power units to the grid, either reducing or increasing them.
On–off peak regulation stands as a primary method for thermal power units to engage in peak regulation. Units operating under variable conditions not only impact the lifespan of the rotor shaft but also induce irreversible changes in other metallic components of the unit. These alterations affect the unit’s overall service life, resulting in on–off costs. The cost associated with unit participation in on–off peak shaving can be represented by Equations (8) and (9):
In the formula, denotes the on–off cost of the unit, denotes the binary state variable of the unit in the peak load phase, denotes whether the unit is in the 0–1 state variable of RPR, denotes whether the unit is in the 0–1 state variable of DPR, denotes whether the unit is in the 0–1 state variable of DPRO, denotes the cost of the unit when it is started, and denotes the cost of the unit when it stops running.
Due to the high requirements of the reaction speed and the on–off time of the unit, the unit with small capacity and a short downtime is generally selected as the thermal power unit with on–off peak regulation [
6].
3. Data-Driven Distributionally Robust Optimal Scheduling Model
With the integration of large-scale new energy sources like wind power into the grid, acknowledging their volatility and uncertainty becomes crucial due to their impact on the power system. Previously, many attributed the abandonment of wind power to technical challenges, citing issues with volatility, uncontrollability, and a lack of grid flexibility for WD and PV power generation. However, in the current context, the abandonment rates, reaching up to 30% to 40%, cannot be solely ascribed to technical limitations. The primary root cause lies in the scenario of excessive installed power capacity and oversupply, posing the challenge of determining priority usage among available resources. Although the Renewable Energy Law [
34] mandates that renewable energy sources have priority access to the grid, the de facto priority in the current power system remains with thermal power generation. This priority status stems from the government’s annual issuance of planned electricity quantities, which has constrained the developmental space for renewable energy. Consequently, to enhance the absorption capacity of wind and photovoltaic power stations and address the uncertainties in their output alongside the adverse effects of forecast errors on the power system, this paper endeavors to construct a novel distributionally robust planning and scheduling strategy for WD–PV fire storage based on comprehensive norms. The schematic diagram of the system is shown in
Figure 2.
3.1. Two-Stage Distributionally Robust Optimization Objective Function
In short-term day-ahead scheduling, the uncertainty associated with WD and PV power generation can be characterized by their output prediction errors. However, this method has limitations when applied to long-term scheduling. This paper aims to devise an optimal scheduling scheme considering the given prediction intervals for renewable energy. It seeks to achieve a cooperative optimization aligning system economy and robustness. Following the framework outlined in Reference [
22], the objective function for the distributionally robust optimization is constructed as follows:
where
represents the probability value of the generation scenario
,
,
represents the power purchase cost of the system, and
represents the operating cost of the unit during peak load balancing.
As evident from Equation (10), the objective function takes the form of a two-stage min–max–min robust optimization. Different from TRO, DRO solves the one-stage scheduling scheme in the worst case by optimizing the lower-bound and estimating its economy by . In the additional stage of the system, the “worst” scenario that the system may encounter is determined by and the probability value is updated; the unity of robustness and economy is realized through continuous iteration in the first and second stages.
The sub-function comprises power purchase costs and operating expenses. The operating costs encompass various elements: coal consumption costs during thermal power unit peak shaving, on–off costs, shaft life expenses, oil injection expenses, deep peak shaving compensation expenses, electric energy storage utilization costs, and wind and photovoltaic abandonment expenses. The expression is as follows:
In Formula (11), is the power purchase cost of the thermal power unit during deep peak shaving, represents the purchase price of a large power grid (RMB/kWh), and represents the power purchased by the system from a large power grid during peak regulation. In Formula (12), , , and , respectively, represent the power generation cost of thermal power unit 1, unit 2, and unit 3, represents the compensation for the unit involved in deep peak regulation, represents the use cost of electric energy storage, and represents the penalty cost of abandoning renewable energy.
- (1)
Cost of power generation for the thermal power unit.
- (2)
Depth peak load compensation cost
where
indicates the electricity compensation price.
- (3)
The cost of storing electric energy
where
represents the charge power of the energy storage,
represents the discharge power of the energy storage, and
represents the unit power of the electric storage cost.
- (4)
The curtailment of abandoning WD and PV
where
denotes the amount of PV curtailment,
denotes the amount of WD curtailment, and
denotes the cost of abandoning WD and PV per unit power.
3.2. Power Constraints
- (1)
System constraints
where
indicates the electrical load of the system.
- (2)
Unit constraints
The upper and lower limits of unit output constraints are
where
and
, respectively, represent the minimum and maximum technical outputs of thermal power units.
The unit climbing power constraint is
where
and
, respectively, represent the maximum upward and downward climbing rates of thermal power units.
Unit on–off constraint
where
denotes the min continuous on–off time for units.
- (3)
Constraints of electric energy storage equipment
where
represents the real-time capacity of electric energy storage,
and
represent the upper and lower limits of electric energy storage capacity,
represents the charging state parameters of electric energy storage,
represents the discharge state parameters of electric energy storage,
and
, respectively, represent the upper and lower limits of the energy storage charging power,
and
, respectively, represent the upper and lower limits of the energy storage discharging power,
represents the energy storage charging efficiency,
represents the energy storage discharging efficiency, and
means to avoid charging and discharging electric energy storage at the same time.
- (4)
WD and PV power constraint
Wind power output constraint
where
denotes the predicted value of wind power output.
Photovoltaic power output constraint
where
denotes the predicted value of photovoltaic output.
3.3. Comprehensive Norm Constraint
In practical operations, the output of renewable energy is subject to uncertainties due to varying weather conditions. Although historical experimental data can provide a probability distribution for its output fluctuations, limitations exist in obtaining an accurate scenario of probability distribution due to data constraints. Therefore, it becomes essential to establish a suitable uncertainty set that aligns the scenario probability closer to actual conditions while fluctuating within a reasonable range. This paper proposes constructing a confidence set based on comprehensive norms to confine the fluctuation range of the probability distribution, thereby forming an uncertainty set for renewable energy output. By creating a confidence interval rooted in comprehensive norms, it constrains the fluctuation scope of the scenario distribution probability, aiming to seek the optimal solution within the worst-case uncertainty set in this scenario.
- (1)
Constraining the probability distribution values of discrete scenarios with the initial probability distribution as the center and using the 1-norm and ∞-norm as constraint conditions result in respective feasible domains and .
where
denotes the scenario
probability value that needs to be updated,
represents
(1-norm),
represents
(∞-norm),
and
represent the maximum deviation value of the probability, and
and
represent the confidence of the probability distribution value.
- (2)
It can be seen from Reference [
35] that the following confidence levels
are satisfied:
where
and
represent the number of discrete scenes and the number of sample scenes.
- (3)
Obtained from Equations (24) and (25)
Update probability scenario values using Equations (14)–(16) and in Equation (10).
3.4. Model Solving Method
Generally speaking, it is difficult to solve a model that combines a min–max–min structure and constraint with uncertain variables [
36]. In handling uncertain variables, this paper employs a method that involves generating 10 scenarios and their corresponding probability values based on 180 days of data [
37]. This is achieved through copula joint WD–PV generation and the application of the k-means clustering algorithm. The resulting probability values serve as the initial probabilities for the analysis. Current solving methods for such problems primarily include Benders’ dual decomposition (BD) and the column and constraint generation algorithm (C&CG) [
38]. In each iteration of BD, the optimal solution of the main problem is passed to the subproblem as a parameter. After solving the subproblem, a new optimal cut or feasibility cut is added to the main problem until the iteration process converges completely. However, the C&CG algorithm will add a new set of constraints and variables to the main problem after each iteration and approximate the optimal solution of the original problem by constantly cutting space. In contrast, the C&CG algorithm retains the second-stage continuous variables in the optimization calculation of the main problem, making the lower bound of the winner problem more compact. Therefore, the C&CG algorithm is often easier to converge than the BD method and requires fewer iterations. So, it is widely used.
5. Conclusions
In this study, we propose a distributionally robust optimal scheduling strategy for a WD–PV thermal storage power system while considering the deep peak shaving of thermal power units. This strategy harnesses the full potential of both conventional thermal power and deep peak shaving units, alongside the flexible utilization of electrical energy storage systems to integrate renewable energy. We introduce a data-driven, distributionally robust optimization approach to solve the uncertainties inherent in renewable energy sources, achieving a synergistic optimization of system robustness and economic efficiency. The efficacy of the proposed strategy is demonstrated through computational examples, leading to the following conclusions:
- (1)
When implementing the strategy delineated in this paper, it is observed that although the economy of the total operating costs for the coal-fired power units in Case 1 is marginally less favorable than those in Cases 2 and 3, the overall economic efficiency of the system is markedly enhanced when taking into account the fuzzy set based on the comprehensive norm distance. Concurrently, the carbon emissions in Case 1 are substantially lower in comparison to Cases 2 and 3. This underscores the profound importance of the cooperative operation of coal-fired units with peak load balancing capabilities and energy storage not only for the economic and low-carbon functioning of the power system but also for the tiered utilization of energy.
- (2)
The deployment of DRO focusing on the comprehensive norm can strike a superior equilibrium between robustness and economic efficiency. Consequently, DRO proves to be more adept at addressing uncertainties associated with variable forces, such as those from WD and PV sources. Moreover, the dual confidence intervals of the comprehensive norm serve as key parameters that reflect the risk preferences of the decision maker. A lower degree of conservatism correlates with reduced operational costs, yet this also diminishes the robustness of the strategy. Thus, decision makers are required to select different confidence intervals that align with their varying risk appetites.
- (3)
It is advisable to prioritize coal-fired units with substantial capacity and a high frequency of significant load fluctuations as the focal point for emission reduction initiatives. In tandem, regular inspections of the turbine rotor are critical to promptly detect any initial cracks and address them effectively. Furthermore, it is essential to bolster investment in the lifecycle maintenance costs of the unit’s rotating shaft to forestall any potential accidents.
The findings of this study offer valuable theoretical insights for the integrated peak-load scheduling of conventional thermal power units, deep peak shaving units, and energy storage, particularly when integrating new energy sources into the power system. However, this paper does not account for the investment costs associated with energy storage capacity or the carbon emissions costs related to low-carbon dispatch strategies, and the estimation of carbon emissions from coal-fired units is somewhat idealized. The forthcoming research will extend its focus to include the capacity configuration and investment cost of energy storage, aiming to maximize scheduling efficiency and economic benefits. Additionally, the impact of carbon emissions from thermal power units during peak load balancing cannot be overlooked, necessitating a thorough examination of carbon emission costs in future analyses. Furthermore, the study will account for demand response variability and uncertainty in load and explore the flexible conversion capabilities of thermal power units. Moreover, it will persist in investigating the feasibility of strategies under the integration of renewable energy sources.