Front. Energy Res. Frontiers in Energy Research Front. Energy Res. 2296-598X Frontiers Media S.A. 1381396 10.3389/fenrg.2024.1381396 Energy Research Original Research A two-stage optimization configuration model for multi-type flexible resource considering spatiotemporal response characteristics Liang and Zhou 10.3389/fenrg.2024.1381396 Liang Yan 1 2 * Zhou Ming 1 1 School of Electrical and Electronic Engineering, North China Electric Power University, Changping, China 2 Economic and Technical Research Institute of Shanxi Electrical Power Company of SGCC, Taiyuan, China

Edited by: Sheng Wang, University of Macau, China

Reviewed by: Xianyou Pan, Shanghai University of Electric Power, China

Wei Fan, Beihang University, China

Shenbo Yang, Beijing University of Technology, China

 *Correspondence: Yan Liang, 3495966833@qq.com 03 07 2024 2024 12 1381396 03 02 2024 04 03 2024 Copyright © 2024 Liang and Zhou. 2024 Liang and Zhou

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.

As China gradually transitions towards a low-carbon energy structure, the proportion of grid-connected new energy sources like wind and solar power continues to increase. To ensure the safe and reliable operation of the power system while meeting the capacity planning for future new energy installations, there is a need for flexible resources with corresponding adjustment capabilities in the power system. In response to this situation, this paper presents an optimization model for the allocation of multiple types of flexible resources that takes into account spatiotemporal response characteristics. Firstly, a flexibility evaluation model is developed based on spatial and temporal response characteristics. Flexibility evaluation indices, such as flexibility average deficit and flexibility coverage index, are constructed. These indices are used for screening nodes with inadequate flexibility in the power system and analyzing the flexibility adequacy at various nodes. Next, the adjustment characteristics of multiple types of flexible resources are analyzed, and a model for their adjustment capabilities is established. Finally, by considering constraints based on time flexibility evaluation indices, a two-stage optimization model for flexible resource allocation is constructed. This model leverages the multiscale matching characteristics between flexibility resources and the fluctuation patterns of new energy sources to guide the allocation of flexible resources at nodes with insufficient flexibility. The effectiveness and applicability of the proposed flexible resource allocation method are validated using the IEEE 9-node system.

 flexible resource spatiotemporal response characteristics flexibility evaluation nodes with insufficient flexibility optimal allocation section-at-acceptance Process and Energy Systems Engineering 1 Introduction

With the rise of the green and low-carbon concept and the continuous deepening of energy transition, to reduce dependence on traditional fossil fuels and decrease carbon emissions (Huang et al., 2023), and in response to the challenges of climate change, the grid integration ratio of new energy sources such as wind and solar will continue to expand, posing higher demands on the flexibility of the power system due to their randomness and uncertainty (Turk et al., 2020).

Currently, there have been numerous studies on the assessment of flexibility in power systems with a high proportion of renewable energy. Existing evaluation indicators and methods on flexibility can be broadly categorized as deterministic and probabilistic. Yan et al. conducted a practical framework study on the flexibility evaluation of power systems at different time scales, considering the uncertainty of renewable energy through the Monte Carlo method (Yan et al., 2020). Lu et al. proposed a novel assessment method based on the probability distribution of flexibility abundance, capable of linearly reflecting the relationship between flexibility shortage and renewable energy reduction (Lu et al., 2018). Tang et al. proposed flexibility evaluation indicators from three perspectives: intra-area supply-demand balance, intra-area power flow distribution, and inter-area transmission capacity (Tang et al., 2020). Stephen et al. introduced a novel metric using the concept of inter-layer pipe bundling to assess the integrated flexibility of natural gas and electricity (Clegg and Mancarella, 2016). Guo et al. introduced a production simulation method using an improved generalized generating function as a flexibility evaluation tool and proposed a flexibility measurement method based on the definition of flexibility and physical mechanisms (Guo et al., 2020). The aforementioned literature collectively indicates that the application of flexibility evaluation methods and indicators can assist power system planners and operators in gaining a better understanding of system flexibility. It enables them to formulate corresponding strategies to cope with the fluctuations in renewable energy output and load. Simultaneously, the assessment of power system flexibility is a multi-layered, multidimensional issue that requires comprehensive consideration of spatial, temporal, and resource-related factors. The traditional flexibility evaluation only analyzes from a single dimension of time or space, without considering the temporal and spatial coupling characteristics of the system, and the evaluation results deviate from the reality. Therefore, it is necessary to conduct a refined evaluation starting from the spatiotemporal characteristics of flexibility resources, providing new insights for enhancing system flexibility and optimizing resource allocation.

In terms of the optimization of flexibility resource allocation, due to the large variability of flexibility resource regulation characteristics at different scales, it is often necessary to consider how to complement the allocation of different types of flexibility resources to meet the different dimensions of the flexibility needs of the new power system. Zhang et al. proposed a multi-flexibility resource collaborative configuration optimization model, considering the Stackelberg game relationship between the flexibility adjustment demands and capabilities of different segments of the power system under various fluctuation scenarios (Ting and Yunna, 2024). Ren et al. comprehensively considered the economic, security, and flexibility aspects of the system, establishing a dual-layer operational planning joint optimization model for flexibility resources (Ren et al., 2020). Ji et al. developed a mixed-integer linear programming model to optimize the design and scheduling of hybrid energy systems, utilizing rooftop photovoltaics and solid waste biomass to meet electricity demands (Ji et al., 2022). Li et al. proposed a new perspective on modeling and planning the flexibility resources at multiple time scales in power systems with high penetration of variable renewable energy. This approach transforms the operational boundaries of flexibility resources into characteristic domains (Li et al., 2022). Zhang et al. proposed a two-dimensional mixed energy storage optimization configuration model for a novel power system with the coupling of multiple flexible resources, aiming to meet the diverse flexibility adjustment requirements at various stages of the novel power system (Zhang et al., 2023). The models in the aforementioned literature often comprehensively consider the coupling of planning and operational aspects, leading to a complex computational process. Existing literature typically plans operations on a single time scale and does not consider grid flexibility at finer time scales. At the same time, existing studies have not fully explored how to optimize power system configuration while simultaneously considering the economic and flexibility resource characteristics. This optimization aims to maximize the potential of flexible resources such as solar-thermal power stations and energy storage.

In response to the aforementioned limitations, this paper proposes a multi-type flexibility resource configuration optimization model that takes into account spatiotemporal response characteristics. Firstly, starting from the net load demand and resource allocation status quo of each node, the flexibility evaluation model based on spatial response characteristics screens the flexibility-deficient nodes and determines the specific location of the allocated flexibility resources. Secondly, considering the multi-timescale fluctuation characteristics of renewable energy/load, the flexibility evaluation constraints based on the time response characteristics, and with the objective of minimizing the investment and construction cost and operation cost of flexibility resources, the optimization model for the allocation of multiple types of flexibility resources is established. Finally, the effectiveness of the configuration optimization model proposed in this paper is verified by IEEE 9-node system simulation.

 2 Spatiotemporal flexibility evaluation model for power systems 2.1 Flexibility evaluation model based on spatial response characteristics

To determine the optimal spatial locations for the configuration of flexible resources, it is necessary to construct a flexibility assessment model based on spatial response characteristics, considering the net load demands and resource allocation status at each node. Firstly, setting the evaluation time scale to 1 h, the power system network structure, resource allocation status, and predicted values of node net load are obtained. With the objective of minimizing the total flexibility deficit in the system, an optimization model for unit combination based on spatial response characteristics is formulated, resulting in the baseline output scheme. The objective function is shown in Formula (1). m i n f = ∑ t = 1 T 1 ∑ m N m P l a c k , t m

Where, P l a c k , t m represents the flexibility deficit at node m at time t , with a positive value indicating insufficient downward ramping capability and a negative value indicating insufficient upward ramping capability.

Next, the average flexibility deficit within the time period is calculated as the evaluation indicator to reflect the flexibility deficiency at each node. This provides a reference for screening nodes with insufficient flexibility. The function is shown in Formula (2). P l a c k m = 1 T 1 ∑ t = 1 T 1 P l a c k , t m

Where, P l a c k m represents the average flexibility deficit at node m .

The constraint conditions are as follows:

1) This paper assume that network losses are ignored and only power balance is considered. The unit output scheduling at each node in the region must meet the predicted net load demand. The function is shown in Formula (3).

∑ k = 1 N k P k , t m + ∑ n ∈ N m P m n , t = P n e t , t m + P l a c k , t m + P e x , t m

Where, P k , t m represents the planned output of unit k at node m at time t , the unit types include all power generation resources and flexibility resources within the system; P n e t , t m represents the net load forecast value for node m at time t . P m n , t represents the power interaction between nodes m and n at time t . P l a c k , t m is the difference between the load forecast and the wind power forecast; P e x , t m is the outgoing power demand of node m at time t .

2) To prevent all units’ reserve output from concentrating on meeting the flexibility requirements of a specific node, the flexibility deficit at each node must satisfy an upper limit constraint. The function is shown in Formula (4).

− P l a c k , t m ≤ P l a c k , max u p , P l a c k , t m ≤ 0 P l a c k , t m ≤ P l a c k , max d n , P l a c k , t m > 0

Where, P l a c k , max u p and P l a c k , max d n are the upper limit values of the upward and downward flexibility deficits at each node.

3) The reserve output provided by each unit must satisfy an upper limit constraint. The function is shown in Formula (5).

0 ≤ R k , t m , u p ≤ ∆ P k , m 0 ≤ R k , t m , d n ≤ ∆ P k , m

Where, R k , t m , u p and R k , t m , d n represent the upward and downward reserve output that unit k under node m can provide at time t; ∆ P k , m denotes the ramping rate of unit k under nod m .

4) The output of each unit needs to satisfy the constraints of output upper and lower limits and ramping constraints. The function is shown in Formula (6).

s k m P k , min m ≤ P k , t m + R k , t m , u p ≤ s k m P k , max m s k m P k , min m ≤ P k , t m − R k , t m , d n ≤ s k m P k , max m − ∆ P k m ≤ P k , t m + R k , t m , u p − P k , t − 1 m − R k , t − 1 m , d n ≤ ∆ P k m

Where, P k , ⁡ max m and P k , ⁡ min m represent the upper and lower limits of the output of unit k under node m ; s k m represents the operating state of unit k under node m .

5) Power interaction between nodes is subject to network constraints. The function is shown in Formula (7).

− P m n , max < P m n , t + R m n , t u p < P m n , max − P m n , max < P m n , t + R m n , t d o w n < P m n , max

Where, P m n , max represents the upper limit of power interaction between node m and node n ; R m n , t u p and R m n , t d o w n represent the upward and downward reserves that nodes m and n can interact with at time t .

6) Each node needs to satisfy the spinning reserve constraint at a 1-h time scale. The function is shown in Formula (8).

R u p , t m = ∑ k = 1 N k R k , t m , u p + ∑ n ∈ N m R m n , t u p ≥ α w p p u p P w p p , t m + α p v u p P p v , t m + α l o a d u p P l o a d , t m R d n , t m = ∑ k = 1 N k R k , t m , d n + ∑ n ∈ N m R m n , t d n ≥ α w p p d n P w p p , t m + α p v d n P p v , t m + α l o a d d n P l o a d , t m

Where, R u p , t m and R d n , t m are the total upward/downward spinning reserve output of node m at time t ; α w p p u p , α w p p d n , α p v u p , α p v d n , α l o a d u p , α l o a d d n are the upward/downward spinning reserve coefficients for wind power, photovoltaic, and load.

 2.2 Flexibility evaluation model based on temporal response characteristics

To determine the capacity of flexibility resources with different response characteristics, it is necessary to build a flexibility assessment model based on time response characteristics. This model should consider the multi-time scale fluctuation characteristics of new energy/load and the benchmark output schemes obtained in Section 2.2. The evaluation time scales are set at 1 h and 15 min.

Assuming the flexibility demand generated by the fluctuation of new energy/load (upward or downward) in the scheduling period for node m is represented by F R m , t , which can be obtained from the stochastic fluctuation characteristics of new energy/load. Then, based on the existing flexibility resources at node m and the benchmark output scheme of the thermal power unit, considering the short-time scale response characteristics of flexibility resources, the ability of node m to provide (upward or downward) flexibility regulation within the scheduling period, denoted as F S m , t , is calculated. The details are as follows: The function is shown in Formula (9). F R m , t u p = F R m , t d n = ξ l o a d P l o a d , t m + ξ w p p P w p p , t m + ξ p v P p v , t m F S m , t u p = ⁡ min   P max m − P t m , ∆ P t m , E max m − E min m / η d i s − P t m ∆ t ∆ t F S m , t d n = ⁡ min   P t m − P min m , ∆ P t m , E t m − E min m / η c h + P t m ∆ t ∆ t

Where, ξ l o a d , ξ w p p , and ξ p v are the random fluctuation error coefficients for load, wind power, and photovoltaic; P max m and P min m are the maximum and minimum values of the operating power; E max m and E min m are the maximum and minimum values of the energy storage capacity; η d i s and η c h are the charge and discharge efficiency; The third term in the calculation formula of F S m , t u p and F S m , t d n only considers energy storage flexibility resources.

Then, define the Flexibility Coverage Index (FCI), which is the ratio of the flexibility regulation capacity to the flexibility demand of the power system during the scheduling period. This index is used to reflect the flexibility abundance of each flexibility-deficient node at different time levels. The specific calculation method for the Flexibility Coverage Index is shown in Figure 1. The specific calculation method for the Flexibility Coverage Index is the area formed by the envelope of the system’s flexibility regulation capacity divided by the area formed by the envelope of the system’s flexibility demand, as shown in Formula (10). F C I m = ∫ t − 1 t F S m , t ∫ t − 1 t F R m , t

Illustration of flexibility coverage index

Where, F C I m represents the Flexibility Coverage Index.

Finally, the minimum value of the appropriate flexibility coverage index according to different risk preferences will be selected as the evaluation index of temporal flexibility, and it will be added as a constraint to the flexibility resource allocation optimization model.

 3 Modeling of multi-type flexibility resource adjustment capacity

The considered flexibility resources in this study include thermal power units, solar-thermal power stations, and energy storage devices such as hydrogen storage, electrochemical storage, and pumped storage.

3.1 Thermal power units

Conventional thermal power units typically exhibit stable operating characteristics, providing continuous power output and meeting power demand over long time scales. Their advantages become more apparent, especially when dealing with significant fluctuations or long-term demands within the system. However, conventional thermal power units cannot achieve frequent charging and discharging in a short time, resulting in a slower response speed and insufficient flexibility when handling high-frequency but low-amplitude random fluctuations. The ramping performance, as well as the upper and lower limits of instantaneous output and output power, of conventional thermal power units affect their flexibility regulation capacity. Therefore, the flexibility of conventional thermal power plants on a time scale of Δ t is as follows. The function is shown in Formula (11), (12). F G , ∆ t + = ⁡ min   P G , ⁡ max − P G , t , R ∙ ∆ t F G , ∆ t − = ⁡ min   P G , t − P G , ⁡ min , R ∙ ∆ t

Where, P G , ⁡ max , P G , ⁡ min , and P G , t represent the maximum, minimum, and current output of the unit at time t, respectively; R is the ramping rate of the unit.

 3.2 Solar-thermal power station

Solar-Thermal Power Stations typically have controllable thermal energy release characteristics, allowing adjustment of the energy release rate and duration according to the power system’s needs. This controllability enables them to flexibly participate in the power system’s response scheduling, providing power according to demand. Their flexibility for upward and downward adjustments on a time scale is shown in formula (13), (14). F c s p , ∆ t + = ⁡ min   P c s p , ⁡ max + , E t − E min + P S F − H T F ∙ ∆ t η R C ∆ t F c s p , ∆ t − = ⁡ min   P c s p , ⁡ max − , E max − E t − P S F − H T F ∙ ∆ t η E H ∆ t

Where, P c s p , ⁡ max + and P c s p , ⁡ max − are the maximum power generation and heat storage power of the solar-thermal power station, respectively; E t , E max , and E min are the upper and lower limits of the thermal storage capacity of the solar-thermal power station at time t ; P S F − H T F is the thermal power collected by the solar field.

 3.3 Energy storage devices

Energy storage devices not only serve as responsive power sources with excellent performance to meet large-scale, system-level applications on the grid side, but also provide bidirectional regulation flexibility by frequently converting electrical energy in a short time. This effectively handles high-frequency but low-amplitude random fluctuations.

Its upward and downward flexibility at time scale Δ t is shown in formula (15), (16): F e s , ∆ t + = ⁡ min   P e s , ⁡ max d − P e s , t , S O C t − S O C min η d i s − P e s , t ∆ t ∆ t F e s , ∆ t + = ⁡ min   P e s , ⁡ max c + P e s , t , S O C max − S O C t / η c h + P e s , t ∆ t ∆ t

Where, P e s , ⁡ max d and P e s , ⁡ max c are the maximum discharge and charge power of the energy storage device, respectively; P e s , t is the power at time t, positive for discharge and negative for charge; S O C t , S O C max , and S O C min are the upper and lower limits of the stored equivalent energy and device capacity at time t, respectively; η d i s and η c h are the efficiency of discharge and charge.

The paper considers pumped storage power stations, hydrogen energy storage systems, and electrochemical energy storage as the main components for configuring energy storage devices. Among them, pumped storage units have strong capacity benefits, fast response rates, and deep response capabilities, effectively addressing peak shaving and valley filling of the net load; hydrogen energy storage systems store electrical energy on a large scale using hydrogen as a medium, which can also achieve peak shaving and valley filling of the net load. At the same time, they can respond rapidly to fluctuations in the net load. However, the current stage has a relatively low energy conversion efficiency. Electrochemical energy storage has an extremely high response rate, can smooth out high-frequency random fluctuations in the net load, but its installed capacity is limited by technical characteristics and cost factors, leading to limited response depth.

 4 Two-stage optimization model for flexible resource configuration 4.1 Two-stage configuration optimization framework system

The paper first filters out nodes with insufficient flexibility based on spatial flexibility assessment indicators. Secondly, considering the constraints of temporal flexibility assessment indicators and the adjustment characteristics of flexibility resources at different time scales, different types of flexibility resources are phased in to meet the flexibility requirements of the power system at both long and short time scales. In the first phase, the focus is on configuring flexibility resources such as pumped storage power stations, solar-thermal power stations, and hydrogen energy storage systems with high adjustment depth and long duration to meet the requirements of long-time scale climbing or peak shaving and valley filling. In the second phase, the emphasis is on flexibility resources with extremely fast response rates, such as electrochemical energy storage, to cope with real-time high-frequency fluctuations in new energy and load. The two-stage framework for flexible resource allocation is shown in Figure 2.

Two-stage framework for flexible resource allocation.

 4.2 Two-stage configuration optimization model for flexibility resource 4.2.1 First stage configuration optimization model

In Stage 1, with a time scale of 1 h, it is mainly used to meet the forecast net load demand of each node, ensuring the full consumption of new energy. The objective function is to minimize the investment construction cost and operating cost. The function is shown in Formula (17). min ⁡ C c b , t o t a l m = C c b , i n v m + C c b , c p e m C c b , i n v m = ∑ i ∈ N c b r 1 + r y i 1 + r y i − 1 c c b , i P P c b , i , ⁡ max m + c c b , i e E c b , i , ⁡ max m C c b , c p e m = ∑ i = 1 T 1 ∑ i ∈ N c b c c b , i o p e P c b , i , t m + ∑ k ∈ N g C g , k , t m

Where, C c b , t o t a l m is the total cost of flexibility resources for node m ; C c b , i n v m and C c b , c p e m are the investment construction cost and operating cost of flexibility resources for node m ; N c b and N g are the sets of flexibility resources and conventional thermal power units configured in Stage 1; r is the discount rate; y i is the lifespan of the ith resource; c c b , i P and c c b , i e are the investment construction costs per unit rated power and unit storage capacity of the ith resource; P c b , i , ⁡ max m and E c b , i , ⁡ max m are the rated power and storage capacity of the ith resource at node m ; c c b , i o p e is the cost of per unit operating power for the ith resource; P c b , i , t m is the operating power of the ith resource at node m in the timeslot t ; C g , k , t m is the operating cost of the ith conventional thermal power unit at node m in the timeslot k .

Constraint conditions are as follows.

1) After configuring flexibility resources, it is necessary to meet the net load demands of nodes and ensure that no flexibility shortfall occurs. The function is shown in Formula (18).

∑ i ∈ N c b P c b , i , t m + ∑ k ∈ N g P g , k , t m + ∑ n ∈ N m P m n , t = P n e t , t m + P e x , t m

Where, P c b , i , t m and P g , k , t m are the operating powers of the flexibility resource i and conventional thermal power unit k at node m in timeslot t during the first stage of configuration.

2) Pumped storage power stations and hydrogen energy storage systems, as energy storage devices, also need to satisfy the energy storage capacity constraints. The function is shown in Formula (19).

E c b , i , t m = E c b , i , t − 1 m + η c b , i c h P c b , i , t m , c h ∆ t − P c b , i , t m , d i s ∆ t η c b , i d i s 0 ≤ E c b , i , t m ≤ E c b , i , ⁡ max m

Where, E c b , i , t m is the equivalent energy stored by the flexibility resource i at node m at timeslot t ; η c b , i c h and η c b , i d i s are the charging and discharging efficiencies of the flexibility resource i ; P c b , i , t m , c h and P c b , i , t m , d i s are the charging and discharging powers of the flexibility resource i at node m at timeslot t . The function is shown in Formula (20). P c b , i , t m , c h = − ⁡ min   P c b , i , t m , 0 P c b , i , t m , d i s = ⁡ max   P c b , i , t m , 0

Additionally, it is necessary to consider the output upper and lower limits, reserve output upper limits, ramping constraints, and rotating reserve constraints of flexibility resources and existing resources, which are not further elaborated here. The transmission power and reserve output between nodes P m n , t , R m n , t u p , R m n , t d o w n are the result values obtained from the optimized base output scenario in Section 2.2. Finally, E c b , i , ⁡ max m represents the configuration results of various types of flexibility resources obtained through optimization, and P c b , i , t m will be passed as a parameter to the second-stage configuration optimization model.

 4.2.2 Second stage configuration optimization model

The short-time scale configuration optimization model for stage 2, based on the initial output plan, generates a baseline output plan for a 15-min time scale. Subsequently, it optimizes the configuration of units with fast regulation rates according to the fluctuation characteristics of new energy and load. The objective function is to minimize the total investment and operating costs, and it is represented in formula (21). min ⁡ C f m , t o t a l m = C f m , i n v m + C f m , o p e m C f m , i n v m = ∑ i ∈ N f m r 1 + r y j 1 + r y j − 1 c f m , j P P f m , j , ⁡ max m + c f m , j e E f m , j , ⁡ max m C f m , o p e m = ∑ t = 1 T 2 ∑ j ∈ N f m c f m , j o p e P f m , j , t m + ∑ i ∈ N c b c c b , i o p e ∆ P c b , i , t m

Where, C f m , t o t a l m is the total cost of flexibility resources for node m ; C f m , i n v m and C f m , o p e m are the investment construction cost and operating cost of flexibility resources for node m ; N f m is the set of electrochemical energy storage units; other variables are consistent with the flexibility resource configuration optimization model in Stage 1, which is not repeated here. P f m , j , ⁡ max m and E f m , j , ⁡ max m are the optimized configuration results of electrochemical energy storage units.

Constraint conditions are as follows.

1) Some of the flexibility resources configured in Stage 1 can also provide a certain degree of regulation capability, and both need to meet the flexibility requirements at a 15-min time scale, it is represented in formula (22).

∑ j ∈ N f m P f m , j , t m + ∑ i ∈ N c b ∆ P c b , i , t m = δ m , t n e t ∆ P c b , i , t m = P c b , i , t m − P c b , i , t − 1 m

Where, P f m , j , t m represents the operating power of the flexibility resource j configured at node m in Stage 2 at time t ; ∆ P c b , i , t m represents the power variation of the flexibility resource i configured at node m in Stage 1 at time t , which still needs to satisfy the output upper limit constraint and ramping constraint of the flexibility resource. Pumped storage and hydrogen storage systems also need to meet the energy storage capacity constraint, which is not elaborated here; δ m , t n e t represents the fluctuation characteristics of the net load at a 15-min time scale, which can be derived from historical data for convenience.

2) The flexibility adjustment capability of each node at a 15-min time scale needs to be greater than the stochastic fluctuations of wind power, photovoltaics, and load at the same time scale. The function is shown in Formula (23).

F u p , t m = ∑ j ∈ N f m F f m , j , t m , u p + ∑ i ∈ N c b F c b , i , t m , u p ≥ β w p p P w p p , t m + β p v P p v , t m + β l o a d P l o a d , t m F d n , t m = ∑ j ∈ N f m F f m , j , t m , d n + ∑ i ∈ N c b F c b , i , t m , d n ≥ β w p p P w p p , t m + β p v P p v , t m + β l o a d P l o a d , t m

Where, F u p , t m and F d n , t m represent the upward/downward flexibility adjustment capabilities of node m at time t ; F f m , j , t m , u p and F f m , j , t m , d n represent the upward/downward flexibility adjustment capabilities of flexibility resource j configured in Stage 2 at time t ; F c b , i , t m , u p and F c b , i , t m , d n represent the upward/downward flexibility adjustment capabilities of flexibility resource i configured in Stage 1 at time t; β w p p , β p v and β l o a d are the random fluctuation coefficients of wind power, photovoltaics, and load at a 15-min time scale, respectively.

3) Flexibility resources need to satisfy output upper and lower limits, ramping constraints, and electrochemical energy storage also needs to consider energy storage capacity constraints; this is not reiterated here.

4) In order to enhance the system’s ability to cope with the random fluctuations of new energy/load, the optimization models for flexibility resource configuration in both stages need to satisfy flexibility indicator constraints. The function is shown in Formula (24).

 F C I m , s u p ≥ ε F I C , s u p F C I m , s d o w n ≥ ε F I C , s d o w n , s ∈ S

Where, S represents the set of considered time scales, including the 1-h time scale and the 15-min time scale; ε F I C , s u p and ε F I C , s d o w n are the thresholds set for the flexibility coverage index at scale s, which is typically chosen between one and 1.2. The specific values of the time flexibility indicator thresholds depend on the decision-maker’s requirements for system reliability and flexibility.

 4.3 Solution procedure

The method for configuring multiple types of flexibility resources, considering spatiotemporal response characteristics, is illustrated in Figure 3. The steps for solving are as follows:

1) Obtain information about the system’s network structure, current resource configuration, node net load forecasts, resource technical parameters, etc.

2) Utilize the flexibility evaluation model based on spatial response characteristics to generate baseline output plans for each node. Calculate the spatial flexibility indicators for each node to determine the node type. If the indicator value for a node is less than the threshold, it is classified as a node with sufficient flexibility; otherwise, it is classified as a node with insufficient flexibility.

3) Based on the multiscale matching characteristics of multiple types of flexibility resources and the fluctuation characteristics of new energy sources, construct an optimization model for flexibility resource capacity configuration. Use the flexibility evaluation model based on temporal response characteristics to guide the formation of flexibility resource configuration schemes that consider time flexibility constraints.

Flowchart of flexible resource configuration steps.

 5 Example analysis 5.1 Example data

In order to simplify the computational complexity and save the simulation computation time, this paper adopts the IEEE 9-node network system shown in Figure 4 for the simulation analysis, and performs the arithmetic simulation on matlab R2022a, and calls the cplex solver through the yalmip toolbox to solve the model, so as to validate the validity and applicability of the flexibility resource allocation method proposed in this paper. Theoretically, the number of nodes can be expanded as needed to accommodate larger systems. In this case, the data for the whole year is scaled by a suitable ratio to obtain the data of wind power output, photovoltaic output and load level at each node for the whole year of 365 × 24 h. Due to the extremely uneven distribution of energy resources and energy demand market in China, the system considers the situation of power transmission. In order to give full play to the renewable energy supply potential in the region, a larger regional power outflow demand is set. The power outgoing demand in one dispatching cycle is shown in Figure 5.

IEEE 9-node network diagram.

Electricity export demand chart.

In addition, the current installed capacity of new energy sources and the configuration status of flexibility resources for each node are presented in Table 1. Currently, only thermal power units and hydropower units are used for flexibility adjustment. Since Nodes 4, 8, and 9 do not have hydropower units, it is assumed in this paper that these nodes have relatively scarce water resources and cannot configure pumped storage units. The cost parameters for different types of flexibility resources are outlined in Table 2.

Current status of new energy installations and flexible resource allocation at each node.

Node Wind power installed Capacity/MW Photovoltaic installed Capacity/MW Thermal power installed Capacity/MW Hydroelectric installed Capacity/MW
1 0 50 10 120
2 0 150 25 160
3 50 100 20 150
4 200 50 30 0
5 0 200 20 0
6 0 200 10 180
7 100 100 20 0
8 200 100 15 200
9 200 100 10 200

Cost parameters for various types of flexible resources.

Type/Parameter Lifecycle (Years) Unit rated power investment cost (104 CNY/MW) Unit energy storage capacity investment cost (104 CNY/MW) Unit operating power operating Cost/MW)
Hydrogen Energy Storage System 10 200 60 200
Pumped Storage Power Station 50 600 84 120
Concentrated Solar Power Station 30 300 0 150
Electrochemical Energy Storage 5 165 127 80
5.2 Example result

This paper mainly considers four types of flexibility resources: solar-thermal, pumped storage, hydrogen storage systems, and electrochemical energy storage. Solar-thermal, pumped storage, and hydrogen storage systems are considered as the initial stage configuration flexibility resources, while electrochemical energy storage is considered as the adjustment stage configuration flexibility resource.

Firstly, using the flexibility assessment model based on spatial response characteristics in Section 1.2, the current flexibility deficits for each node are obtained, and nodes with insufficient flexibility are identified. The flexibility deficits in each time period of the scheduling cycle are shown in Figure 5. The mutual support of electrical energy and flexibility between nodes is illustrated in Figure 6.

Flexibility deficit chart for each node within the scheduling period.

Combining Figure 6 and Figure 7, it is evident that there is no upward flexibility deficit in the system. However, all nodes except Node four experience downward flexibility deficits, concentrated in the time period from 10 to 14. Nodes such as 3, 8, and 9 have abundant hydroelectric resources, resulting in surplus flexibility adjustment capabilities. These nodes will provide some upward and downward reserves through regional power grid support for other nodes. Based on this, the average flexibility deficit and spatial flexibility index for each node are calculated, with results presented in Table 3.

Electricity and flexibility transfer diagram among nodes.

Average flexibility deficit across nodes.

Indicator/Node 1 2 3 4 5 6 7 8 9
P l a c k m M W / h 1.91 2.08 1.49 0 6.04 2.08 3.42 2.08 4.17

Assuming the threshold for the average flexibility deficit at each node is set to 2 MW/h, nodes with insufficient flexibility, namely, 2, 5, 6, 7, 8, and 9, are filtered. Subsequently, a two-stage configuration optimization model is employed to allocate the required types of flexibility resources for each node. The flexibility coverage index thresholds for Stage 1 and Stage 2 of the flexibility resource allocation optimization model are set to 1.05 and 1.1, respectively. The optimized flexibility resource allocation results for each node are presented in Figure 8, and the corresponding investment and construction costs as well as operating costs are detailed in Table 4.

Optimization results for flexible resource allocation at each node.

Investment and operating costs for each node.

Node\cost Total Cost (104 CNY) Investment and construction Cost (104 CNY) Operating Cost (104 CNY)
1 0 0 0
2 16057.96 2775.19 13282.77
3 0 0 0
4 0 0 0
5 8037.12 4295.99 3741.13
6 11834.24 4689.42 7144.82
7 7350.75 3422.61 3928.14
8 22455.1 6232.61 16222.49
9 19834.91 5568.35 14266.56

From Figure 8, it can be observed that due to the ability of hydrogen storage systems and pumped storage power stations to store a large amount of electrical energy while also considering economic factors, when the flexibility deficit is significant, nodes will prioritize the configuration of pumped storage power stations. Nodes that cannot configure pumped storage power stations will choose to configure hydrogen storage systems. Although solar-thermal power stations have a certain degree of flexibility, they are still limited by sunlight conditions. Therefore, their regulation capabilities are affected during nighttime or cloudy weather, making them somewhat limited, with smaller configuration capacities at each node. Additionally, to mitigate the stochastic fluctuations of net load on shorter time scales, nodes will configure a certain capacity of electrochemical energy storage, and the configuration capacity will increase with the increase in the installed capacity of renewable energy sources in the region where the node is located.

Next, the flexibility coverage index threshold is set to 1, 1.05, 1.1, 1.15, and 1.2, respectively. The total investment and construction costs, as well as operating costs, for configuring electrochemical energy storage at each node in Phase 2 are calculated and presented in Table 5, and the configuration optimization results are shown in Figure 9. The configuration results are analyzed based on the threshold of 1.1, and the relative changes in the configured node energy storage capacity are depicted in Figure 10.

Total investment and operating costs for node system in stage 2.

Threshold/Cost Total Cost (104 CNY) Investment and construction Cost (104 CNY) Total operating Cost (104 CNY)
1 13682.12 6320.94 7361.17
1.05 13788.23 6427.06 7361.17
1.1 13894.35 6533.17 7361.17
1.15 14000.46 6639.29 7361.17
1.2 14106.58 6745.41 7361.17

Optimization results for flexible resource allocation at each node in stage 2 under different thresholds.

The relative energy storage capacity of flexibility resources allocated to nodes in stage 2 at different thresholds.

From Table 5, it can be observed that with the increase in the flexibility coverage index threshold, the total investment and construction costs continuously increase, while the total operating costs remain constant. This is because, with the increase in the threshold, the requirement for reserve resources to better cope with the uncertainty of load and wind-solar output increases, necessitating the expansion of reserve capacity, thus increasing the overall construction costs. As for operating costs, they are only related to the operating power of resources and the unit operating cost, independent of changes in the configuration of reserve capacity. Therefore, regardless of how the configuration is altered, operating costs remain constant.

 6 Conclusion

To enhance the safety and reliability of future renewable energy power systems, this paper proposes a spatiotemporal response-aware multi-type flexibility resource configuration optimization model. Considering the regulating characteristics of various flexibility resources, the paper optimizes the configuration of flexibility resources in the power system to obtain an optimal planning solution. Simulation results indicate:

1) Based on the set average flexibility deficit indicator method, this paper conducts a spatial flexibility assessment, allowing system operators or planners to concentrate resources and attention on nodes that most require additional flexibility resources. This targeted approach facilitates the allocation of appropriate types of flexibility resources to specific nodes.

2) The two-stage configuration optimization model for flexibility resource capacity proposed in this paper, by considering the different regulation characteristics and mutual cooperation of multiple types of flexibility resources, enables flexible decision-making. This phased configuration of flexibility resources helps the power system better accommodate and integrate new energy, reducing dependence on traditional fossil fuels.

3) The proposed flexibility coverage index threshold in this paper can effectively guide the flexibility resource allocation schemes for nodes with insufficient flexibility. Additionally, as the flexibility coverage index threshold increases, the requirements for the capacity of flexibility resource allocation increase, leading to an overall increase in investment and construction costs.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary materials, further inquiries can be directed to the corresponding author.

 Author contributions

YL: Conceptualization, Methodology, Validation, Writing–original draft, Conceptualization, Methodology, Validation, Writing–original draft, MZ: Visualization, Writing–review and editing, Visualization, Writing–review and editing.

 Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work is supported by Science and Technology Project of Shanxi Electric Power Company (520533230006).

We would like to thank Shanxi Electric Power Company for their help in the design of our study.

 Conflict of interest

Author YL was employed by Economic and Technical Research Institute of Shanxi Electrical Power Company of SGCC.

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 Science and Technology Project of Shanxi Electric Power Company. The funder had the following involvement in the study: design of our study.

 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 Clegg S. Mancarella P. (2016). Integrated electrical and gas network flexibility assessment in low-carbon multi-energy systems. IEEE Trans. Sustain. Energy 7 (2), 718–731. 10.1109/tste.2015.2497329 Degefa M. Z. Sperstad I. B. Sæle H. (2021). Comprehensive classifications and characterizations of power system flexibility resources. Electr. Power Syst. Res. 194, 107022. 10.1016/j.epsr.2021.107022 Guo Z. Zheng Y. Li G. (2020). Power system flexibility quantitative evaluation based on improved universal generating function method: a case study of Zhangjiakou. Energy 205, 117963. 10.1016/j.energy.2020.117963 Hong TANG Wang S. Hangxin L. I. (2021). Flexibility categorization, sources, capabilities and technologies for energy-flexible and grid-responsive buildings: state-of-the-art and future perspective. Energy 219, 119598. 10.1016/j.energy.2020.119598 Huang Y. Kuldasheva Z. Bobojanov S. Djalilov B. Salahodjaev R. Abbas S. (2023). Exploring the links between fossil fuel energy consumption, industrial value-added, and carbon emissions in G20 countries. Environ. Sci. Pollut. Res. 30 (4), 10854–10866. 10.1007/s11356-022-22605-9 Ji L. Wu Y. Liu Y. Sun L. Xie Y. Huang G. (2022). Optimizing design and performance assessment of a community-scale hybrid power system with distributed renewable energy and flexible demand response. Sustain. Cities Soc. 84, 104042. 10.1016/j.scs.2022.104042 Jiawei FENG Yang J. Wang H. Ji H. Okoye M. O. Cui J. (2020). Optimal dispatch of high-penetration renewable energy integrated power system based on flexible resources. Energies 13 (13), 3456. 10.3390/en13133456 Jie JIAN Peng L. I. Haoran J. I. Bai L. Yu H. Xi W. (2022). DLMP-based quantification and analysis method of operational flexibility in flexible distribution networks. IEEE Trans. Sustain. Energy 13 (4), 2353–2369. 10.1109/tste.2022.3197175 Li H. Zhang N. Bao W. (2022). Modeling and planning of multi-timescale flexible resources in power systems. CSEE J. Power Energy Syst., 1–10. 10.17775/CSEEJPES.2021.04850 Lu Z. Li H. Qiao Y. (2018). Probabilistic flexibility evaluation for power system planning considering its association with renewable power curtailment. IEEE Trans. Power Syst. 33 (3), 3285–3295. 10.1109/tpwrs.2018.2810091 Ren Z. Guo H. Yang P. Zuo G. Zhao Z. (2020). Bi-level optimal allocation of flexible resources for distribution network considering different energy storage operation strategies in electricity market. IEEE Access 8, 58497–58508. 10.1109/access.2020.2983042 Ren Z. Guo H. Yang P. Zuo G. Zhao Z. (2020a). Bi-level optimal allocation of flexible resources for distribution network considering different energy storage operation strategies in electricity market. IEEE Access 8, 58497–58508. 10.1109/access.2020.2983042 Ren Z. Guo H. Yang P. Zuo G. Zhao Z. (2020b). Bi-level optimal allocation of flexible resources for distribution network considering different energy storage operation strategies in electricity market. IEEE Access 8, 58497–58508. 10.1109/access.2020.2983042 Tang X. Hu Y. Chen Z. You G. (2020). Flexibility evaluation method of power systems with high proportion renewable energy based on typical operation scenarios. Electronics 9 (4), 627. 10.3390/electronics9040627 Ting Z. Yunna W. (2024). Collaborative allocation model and balanced interaction strategy of multi flexible resources in the new power system based on Stackelberg game theory. Renew. Energy 220, 119714. 10.1016/j.renene.2023.119714 Turk A. Wu Q. Zhang M. Østergaard J. (2020). Day-ahead stochastic scheduling of integrated multi-energy system for flexibility synergy and uncertainty balancing. Energy 196, 117130. 10.1016/j.energy.2020.117130 Xiaowei Y. U. Xu LING Xiaogang ZHOU Li K. Feng X. (2022). Flexibility evaluation and index analysis of distributed generation planning for grid-source coordination. Front. Energy Res. 10. 10.3389/fenrg.2022.1019352 Yan X. Jiang H. Gao Y. Li J. Abbes D. (2020). “Practical flexibility analysis on europe power system with high penetration of variable renewable energy,” in 2020 IEEE Sustainable Power and Energy Conference (iSPEC), Chengdu, China, November, 2020, 395–402. 10.1109/ispec50848.2020.9351168 Yang H. Tian X. Liu F. Liu L. Li L. Wang Q. (2023). A multi-objective dispatching model for a novel virtual power plant considering combined heat and power units, carbon recycling utilization and flexible load response. Front. Energy Res. 11, 1332474. 10.3389/fenrg.2023.1332474 Zhang T. Ma Y. Wu Y. Yi L. (2023). Optimization configuration and application value assessment modeling of hybrid energy storage in the new power system with multi-flexible resources coupling. J. Energy Storage 62, 106876. 10.1016/j.est.2023.106876 Zhijie L. I. U. Zhu S. Zhang P. Zhang Z. (2022). “Flexibility evaluation index system and calculation method of active distribution network with high proportion of renewable energy,” in 2022 Asian Conference on Frontiers of Power and Energy (ACFPE), Chengdu, China, October, 2022, 720–724. 10.1109/acfpe56003.2022.9952278