Calculation Method for Power System Flexibility Resources

Abstract

The present disclosure provides a method for calculating flexibility resources of a power system, belonging to the technical field of power system planning. In response to receiving a request for calculating flexibility resources of a target power system, resource configuration data of the target power system is obtained from a database; the calculation of flexibility resources of the target power system is divided into M stages according to the construction time sequence; for the i-th stage, according to the calculation results of flexibility resources in the 1st stage to the (i-1)-th stage and the i-th set of sub-data, the gain of flexibility retrofit resources of thermal power units, the peak shaving resource gain of energy storage systems, the consumption of flexibility retrofit resources of thermal power units, the resource consumption of energy storage systems, and the curtailment rate of wind and light in the i-th stage are calculated; an objective function for the i-th stage is established to obtain M objective functions; and under predetermined constraint conditions, the M objective functions are solved to obtain M calculation results of flexibility resources corresponding to the M stages; the M calculation results of flexibility resources are visually displayed.

Description

Technical Field

[0001] The present disclosure relates to the technical field of power system planning, and more particularly, to a method for calculating flexibility resources of a power system. Background Art

[0002] The development and utilization of renewable energy play an irreplaceable role and significance in promoting the sustainable development of the economy and society, protecting the ecological environment, and addressing climate change. However, with the large-scale grid connection of renewable energy mainly based on wind power and photovoltaic power generation, its inherent randomness and volatility require the introduction of corresponding flexibility resources for balancing.

[0003] In order to cope with the current complex and changeable peak shaving pressure, in related technologies, the flexibility resources for absorbing the grid connection fluctuations of renewable energy are currently mainly divided into two types. One is an energy storage system with a relatively fast response speed, but due to the too high price per unit capacity, it is not suitable for large-scale configuration. The other is to transform traditional thermal power units to improve their peak shaving capacity, but thermal power units have their peak shaving limits and are not suitable for occasions with too high renewable energy penetration. Therefore, how to utilize the characteristics of various flexibility resources for complementarity and construct a reasonable configuration combination has become a difficult problem that urgently needs to be solved in the current flexibility resource configuration. Summary of the Invention

[0004] In view of this, the present disclosure provides a method for calculating flexibility resources of a power system.

[0005] One aspect of the present disclosure provides a method for calculating flexibility resources of a power system, including: in response to receiving a flexibility resource calculation request for a target power system, obtaining resource configuration data of the target power system from a database, where the resource configuration data includes M groups of sub-data divided according to the construction time sequence of the target power system, and M is a positive integer greater than or equal to 2; dividing the calculation of flexibility resources of the target power system into M stages according to the construction time sequence, where the M stages correspond to the M groups of sub-data; for the i-th stage, calculating the flexibility retrofit resource gain of thermal power units in the i-th stage, the peak shaving resource gain of energy storage systems in the i-th stage, the flexibility retrofit resource consumption of thermal power units in the i-th stage, the resource consumption of energy storage systems in the i-th stage, and the wind and light curtailment rate in the i-th stage according to the flexibility resource calculation results of the 1st stage to the (i - 1)-th stage and the i-th group of sub-data, where i is a positive integer and 1 ≤ i ≤ M; establishing an objective function for the i-th stage according to the flexibility retrofit resource gain of thermal power units in the i-th stage, the peak shaving resource gain of energy storage systems in the i-th stage, the flexibility retrofit resource consumption of thermal power units in the i-th stage, the resource consumption of energy storage systems in the i-th stage, the wind and light curtailment rate in the i-th stage, and the objective function of the (i - 1)-th stage, to obtain M objective functions; and solving the M objective functions under predetermined constraints to obtain M flexibility resource calculation results corresponding to the M stages, where the flexibility resource calculation results include the newly added energy storage capacity, the newly added converter power capacity, and the newly added thermal power unit capacity;

[0006] Another aspect of the present disclosure provides a power system flexibility resource calculation device, including: an acquisition module, configured to, in response to receiving a flexibility resource calculation request for a target power system, obtain resource configuration data of the target power system from a database, where the resource configuration data includes M groups of sub-data divided according to the construction time sequence of the target power system, and M is a positive integer greater than or equal to 2; a division module, configured to divide the flexibility resource calculation of the target power system into M stages according to the construction time sequence, where the M stages correspond to the M groups of sub-data; a calculation module, configured to, for the i-th stage, calculate the flexibility retrofit resource gain of thermal power units, the peak shaving resource gain of energy storage systems, the flexibility retrofit resource consumption of thermal power units, the resource consumption of energy storage systems, and the wind and light curtailment rate in the i-th stage according to the flexibility resource calculation results of the 1st stage to the (i - 1)-th stage and the i-th group of sub-data, where i is a positive integer and 1 ≤ i ≤ M; a construction module, configured to establish an objective function for the i-th stage according to the flexibility retrofit resource gain of thermal power units, the peak shaving resource gain of energy storage systems, the flexibility retrofit resource consumption of thermal power units, the resource consumption of energy storage systems, the wind and light curtailment rate in the i-th stage, and the objective function of the (i - 1)-th stage, and obtain M objective functions; a solution module, configured to solve the M objective functions under predetermined constraint conditions to obtain M flexibility resource calculation results corresponding to the M stages, where the flexibility resource calculation results include the newly added energy storage capacity, the newly added converter power capacity, and the newly added thermal power unit capacity; a display module, configured to visually display the M flexibility resource calculation results.

[0007] Another aspect of the present disclosure provides an electronic device, including:

[0008] One or more processors;

[0009] A memory for storing one or more programs,

[0010] where, when the one or more programs are executed by the one or more processors, the one or more processors are caused to implement the method as described above.

[0011] Another aspect of the present disclosure provides a computer-readable storage medium storing computer-executable instructions, and the instructions are used to implement the method as described above when executed.

[0012] Another aspect of the present disclosure provides a computer program product, where the computer program product includes computer-executable instructions, and the instructions are used to implement the method as described above when executed.

[0013] According to the embodiments of the present disclosure, resource configuration data of the power system is obtained from a database, and the calculation of power system flexibility resources is divided into different stages according to the construction sequence of the power system and taking into account the technical and economic changes in the process of energy storage and thermal power flexibility transformation; according to the resource configuration data of the current stage and the flexibility resource calculation results of the historical stage, the resource gain of the thermal power unit flexibility transformation, the peak-shaving resource gain of the energy storage system, the resource consumption of the thermal power unit flexibility transformation and the resource consumption of the energy storage system of the current stage resource configuration are analyzed and calculated, and the resource gain and resource consumption of the thermal power unit and the energy storage system are comprehensively considered; on the basis of considering the wind and solar power abandonment rate in the current stage, the overall economic benefits of the current stage and the previous stage are used as the objective function of the current stage, and the flexibility resource calculation results of the current stage are used as the initial conditions of the next stage for iterative optimization to obtain the objective functions of all stages; under predetermined constraints, the objective functions of all stages are solved to obtain the flexibility resource calculation results of all stages, thereby improving the accuracy and rationality of the flexibility resource calculation of the power system. This method is suitable for large-scale configuration of energy storage systems and also for situations where the penetration rate of renewable energy is too high. It explores the intrinsic connection between the resource configuration data of the power system and the flexibility resource calculation results, that is, the intrinsic connection between the resource configuration data of the power system and the newly added energy storage capacity, the newly added converter power capacity and the newly added thermal power unit capacity. At the same time, it uses computer technology to automatically generate and display the flexibility resource calculation results, thereby improving the efficiency of the flexibility resource calculation of the power system. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] The above and other objects, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which:

[0015] Figure 1 An exemplary system architecture to which the method and apparatus for calculating flexibility resources in a power system of the present disclosure can be applied is schematically shown;

[0016] Figure 2 A flowchart of a method for calculating flexibility resources of a power system according to an embodiment of the present disclosure is schematically shown;

[0017] Figure 3 The peak-shaving process of the thermal power unit flexibility transformation according to the embodiment of the present disclosure is schematically shown;

[0018] Figure 4 A flowchart of a method for calculating flexibility resources of a power system according to another embodiment of the present disclosure is schematically shown;

[0019] Figure 5 The schematic diagram shows the power generation of a renewable energy power generation base after it is put into operation in a prospective year according to an embodiment of the present disclosure;

[0020] Figure 6 Schematically shows the grid-connected load of a certain area in 2022 according to an embodiment of the present disclosure;

[0021] Figure 7 Schematically shows a block diagram of a power system flexibility resource calculation device according to an embodiment of the present disclosure; and

[0022] Figure 8 Schematically shows a block diagram of an electronic device suitable for implementing a power system flexibility resource calculation method according to an embodiment of the present disclosure. Detailed implementation manners

[0023] Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are merely exemplary and are not intended to limit the scope of the present disclosure. In the following detailed description, for the sake of explanation, many specific details are set forth to provide a comprehensive understanding of the embodiments of the present disclosure. However, obviously, one or more embodiments can also be implemented without these specific details. In addition, in the following description, descriptions of well-known structures and technologies are omitted to avoid unnecessarily confusing the concepts of the present disclosure.

[0024] The terms used herein are merely for describing specific embodiments and are not intended to limit the present disclosure. The terms "including", "comprising" and the like used herein indicate the presence of the described features, steps, operations and / or components, but do not exclude the presence or addition of one or more other features, steps, operations or components.

[0025] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification and should not be interpreted in an idealized or overly rigid manner.

[0026] In the case of using expressions such as "at least one of A, B, and C", generally, it should be interpreted according to the meaning commonly understood by those skilled in the art (for example, "a system having at least one of A, B, and C" should include, but is not limited to, a system having only A, only B, only C, having A and B, having A and C, having B and C, and / or having A, B, and C, etc.).

[0027] In the embodiments of the present disclosure, in terms of the collection, update, analysis, processing, use, transmission, provision, disclosure, storage, etc. of the data involved (for example, including but not limited to user personal information), they all comply with the provisions of relevant laws and regulations, are used for legal purposes, and do not violate public order and good customs. In particular, necessary measures are taken for user personal information to prevent illegal access to user personal information data and to maintain the security of user personal information and network security.

[0028] The development and utilization of renewable energy play an irreplaceable role and significance in promoting the sustainable development of the economy and society, protecting the ecological environment, and addressing climate change. However, with the large-scale grid connection of renewable energy mainly based on wind power and photovoltaic power generation, its inherent randomness and volatility require the introduction of corresponding flexibility resources for balancing.

[0029] To cope with the current complex and changeable peak shaving pressure, in related technologies, the flexibility resources for absorbing the grid connection fluctuations of renewable energy are currently mainly divided into two types. One is the energy storage system with a relatively fast response speed, but due to the too high price per unit capacity, it is not suitable for large-scale configuration. The other is to transform traditional thermal power units to improve their peak shaving capacity, but thermal power units have their peak shaving limits and are not suitable for occasions with too high penetration of renewable energy. Therefore, how to utilize the characteristics of various flexibility resources for complementarity and construct a reasonable configuration combination has become a difficult problem that urgently needs to be solved in the current flexibility resource configuration.

[0030] In addition, in related technologies, the configuration of flexibility resources mostly adopts static planning, that is, a one-step planning from the starting year to the target year, which cannot consider the technical and economic changes during the process of energy storage and thermal power flexibility transformation, such as the decrease in configuration cost and the innovation of energy storage technology, resulting in relatively high configuration costs and poor economy.

[0031] At present, the scale of renewable energy is growing rapidly, thermal power units are gradually undergoing transformation, and the cost of energy storage is also continuously decreasing. Therefore, when configuring flexibility resources, the embodiments of the present disclosure consider the development and changes in terms of technology and economy during the construction process.

[0032] Based on this, embodiments of the present disclosure provide a method for calculating flexibility resources of a power system, including: in response to receiving a flexibility resource calculation request for a target power system, obtaining resource configuration data of the target power system from a database, where the resource configuration data includes M groups of sub-data divided according to the construction time sequence of the target power system, and M is a positive integer greater than or equal to 2; dividing the calculation of the flexibility resources of the target power system into M stages according to the construction time sequence, where the M stages correspond to the M groups of sub-data; for the i-th stage, calculating the gain of the flexibility transformation resources of the thermal power unit in the i-th stage, the peak shaving resource gain of the energy storage system in the i-th stage, the consumption of the flexibility transformation resources of the thermal power unit in the i-th stage, the resource consumption of the energy storage system in the i-th stage, and the wind and light curtailment rate in the i-th stage according to the flexibility resource calculation results of the first stage to the (i - 1)-th stage and the i-th group of sub-data, where i is a positive integer and 1 ≤ i ≤ M; establishing an objective function for the i-th stage according to the gain of the flexibility transformation resources of the thermal power unit in the i-th stage, the peak shaving resource gain of the energy storage system in the i-th stage, the consumption of the flexibility transformation resources of the thermal power unit in the i-th stage, the resource consumption of the energy storage system in the i-th stage, the wind and light curtailment rate in the i-th stage, and the objective function of the (i - 1)-th stage, to obtain M objective functions; solving the M objective functions under predetermined constraint conditions to obtain M flexibility resource calculation results corresponding to the M stages, where the flexibility resource calculation results include the newly added energy storage capacity, the newly added converter power capacity, and the newly added thermal power unit capacity; and visually displaying the M flexibility resource calculation results.

[0033] Figure 1 Schematically shows an exemplary system architecture 100 to which the method and apparatus for calculating flexibility resources of a power system according to embodiments of the present disclosure can be applied. It should be noted that, Figure 1 The shown is only an example of the system architecture to which embodiments of the present disclosure can be applied, to help those skilled in the art understand the technical content of the present disclosure, but it does not mean that embodiments of the present disclosure cannot be used in other devices, systems, environments or scenarios.

[0034] As Figure 1 shown, the system architecture 100 according to this embodiment may include a first terminal device 101, a second terminal device 102, a third terminal device 103, a network 104, and a server 105. The network 104 is used as a medium to provide a communication link between the first terminal device 101, the second terminal device 102, the third terminal device 103, and the server 105. The network 104 may include various connection types, such as wired and / or wireless communication links, etc.

[0035] Users can use the first terminal device 101, the second terminal device 102, and the third terminal device 103 to interact with the server 105 via the network 104 to receive or send messages, etc. Various communication client applications can be installed on the first terminal device 101, the second terminal device 102, and the third terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, and / or social platform software, etc. (for example only).

[0036] The first terminal device 101, the second terminal device 102, and the third terminal device 103 can be various electronic devices with a display screen and supporting web browsing, including but not limited to smartphones, tablets, laptop portable computers, and desktop computers, etc.

[0037] The server 105 can be a server providing various services, such as a background management server that supports the websites browsed by users using the first terminal device 101, the second terminal device 102, and the third terminal device 103 (for example only). The background management server can analyze and process data such as received user requests, etc., and feedback the processing results (such as web pages, information, or data, etc. obtained or generated according to user requests) to the terminal device.

[0038] It should be noted that the power system flexibility resource calculation method provided by the embodiments of the present disclosure can generally be executed by the server 105. Correspondingly, the power system flexibility resource calculation device provided by the embodiments of the present disclosure can generally be set in the server 105. The power system flexibility resource calculation method provided by the embodiments of the present disclosure can also be executed by a server or a server cluster different from the server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or the server 105. Correspondingly, the power system flexibility resource calculation device provided by the embodiments of the present disclosure can also be set in a server or a server cluster different from the server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or the server 105. Or, the power system flexibility resource calculation method provided by the embodiments of the present disclosure can also be executed by the first terminal device 101, the second terminal device 102, or the third terminal device 103, or can also be executed by other terminal devices different from the first terminal device 101, the second terminal device 102, or the third terminal device 103. Correspondingly, the power system flexibility resource calculation device provided by the embodiments of the present disclosure can also be set in the first terminal device 101, the second terminal device 102, or the third terminal device 103, or set in other terminal devices different from the first terminal device 101, the second terminal device 102, or the third terminal device 103.

[0039] For example, the database can originally be stored in any one of the first terminal device 101, the second terminal device 102, or the third terminal device 103 (e.g., the first terminal device 101, but not limited thereto), or stored on an external storage device and can be imported into the first terminal device 101. Then, the first terminal device 101 can execute the power system flexibility resource calculation method provided by the embodiments of the present disclosure locally, or send the database to other terminal devices, servers, or server clusters, and the other terminal devices, servers, or server clusters that receive the database execute the power system flexibility resource calculation method provided by the embodiments of the present disclosure.

[0040] It should be understood that Figure 1 the numbers of terminal devices, networks, and servers in

[0041] Figure 2 is merely illustrative. According to the implementation requirements, there can be any number of terminal devices, networks, and servers.

[0042] As Figure 2 shown, the method 200 includes operations S210 to S260.

[0043] In operation S210, in response to receiving a flexibility resource calculation request for a target power system, obtain the resource configuration data of the target power system from the database, where the resource configuration data includes M groups of sub-data divided according to the construction time sequence of the target power system, and M is a positive integer greater than or equal to 2.

[0044] In operation S220, divide the flexibility resource calculation of the target power system into M stages according to the construction time sequence, where the M stages correspond to the M groups of sub-data.

[0045] In operation S230, for the i-th stage, calculate the flexibility retrofit resource gain of thermal power units in the i-th stage, the peak shaving resource gain of energy storage systems in the i-th stage, the flexibility retrofit resource consumption of thermal power units in the i-th stage, the resource consumption of energy storage systems in the i-th stage, and the curtailment rate of wind and light in the i-th stage according to the flexibility resource calculation results of the first stage to the (i - 1)-th stage and the i-th group of sub-data, where i is a positive integer and 1 ≤ i ≤ M.

[0046] In operation S240, establish the objective function for the i-th stage according to the flexibility retrofit resource gain of thermal power units in the i-th stage, the peak shaving resource gain of energy storage systems in the i-th stage, the flexibility retrofit resource consumption of thermal power units in the i-th stage, the resource consumption of energy storage systems in the i-th stage, the curtailment rate of wind and light in the i-th stage, and the objective function of the (i - 1)-th stage, and obtain M objective functions.

[0047] In operation S250, under predetermined constraints, M objective functions are solved to obtain M flexibility resource calculation results corresponding to M stages, where the flexibility resource calculation results include newly added energy storage capacity, newly added converter power capacity, and newly added thermal power unit capacity.

[0048] In operation S260, the M flexibility resource calculation results are visually displayed.

[0049] According to an embodiment of the present disclosure, the resource configuration data of the target power system may include the data required for the flexibility transformation of thermal power units and the data required for the transformation of energy storage systems. The data required for the flexibility transformation of thermal power units and the data required for the transformation of energy storage systems may be respectively divided into M groups of sub-data according to the construction time sequence of the target power system.

[0050] According to an embodiment of the present disclosure, for the i-th stage, the i-th group of sub-data may include the annual grid-connected load data of the planned year of the i-th stage, the renewable energy power generation data of the planned year of the i-th stage, the maximum number of units to be transformed in the i-th stage, the maximum scale of energy storage configuration in the i-th stage, the discount rate, the theoretical operation life of the transformed unit, the investment cost per unit capacity for the flexibility transformation of thermal power units, the lifetime maintenance cost per unit capacity of the thermal power flexibility transformation unit, the minimum technical output of the thermal power unit during three peak shaving processes, the maximum output of the thermal power unit, the quadratic term, linear term, and constant term coefficients of the consumption characteristic function of the motor unit, the unit prices of coal and fuel oil, the rotor loss life consumed when deeply reducing the output of the unit, the purchase price of the thermal power unit rotor, the average hourly fuel injection volume when the unit uses oil for peak shaving, the pollutant discharge fee generated by burning unit coal and unit fuel oil, the investment cost per unit power of the converter, the investment cost per unit capacity of the energy storage system, the theoretical operation life of the energy storage system, the annual operation and maintenance cost per unit capacity of the energy storage system, the on-grid peak-valley electricity price, the deep peak shaving compensation electricity price for the flexibility transformation unit, the charge-discharge efficiency of the energy storage system, the ramp rate limit before and after the unit transformation, the flexibility transformation upper limit of the unit, the charge-discharge power limit of the energy storage system in the i-th stage, the capacity limit of the energy storage system in the i-th stage, etc.

[0051] According to an embodiment of the present disclosure, the wind and light abandonment rate may be the sum of the wind abandonment rate and the light abandonment rate.

[0052] According to an embodiment of the present disclosure, for example, the resource allocation data may include five groups of sub-data divided according to the construction time sequence of the target power system. The calculation of the flexibility resources of the target power system can be divided into five stages according to the construction time sequence, where the five stages correspond to the five groups of sub-data. For the fourth stage, based on the calculation results of the flexibility resources in the first to third stages and the fourth group of sub-data, the gain of the flexibility transformation resources of the thermal power units in the fourth stage, the peak shaving resource gain of the energy storage system in the fourth stage, the consumption of the flexibility transformation resources of the thermal power units in the fourth stage, the resource consumption of the energy storage system in the fourth stage, and the wind and light abandonment rate in the fourth stage are calculated. Based on the gain of the flexibility transformation resources of the thermal power units in the fourth stage, the peak shaving resource gain of the energy storage system in the fourth stage, the consumption of the flexibility transformation resources of the thermal power units in the fourth stage, the resource consumption of the energy storage system in the fourth stage, the wind and light abandonment rate in the fourth stage, and the objective function of the third stage, an objective function for the fourth stage is established, and five objective functions are obtained. Under the predetermined constraint conditions, the five objective functions are solved to obtain five flexibility resource calculation results corresponding to the five stages.

[0053] According to an embodiment of the present disclosure, according to an embodiment of the present disclosure, the resource allocation data of the power system is obtained from the database. According to the construction time sequence of the power system, considering the technical and economic changes in the process of energy storage and thermal power flexibility transformation, the calculation of the flexibility resources of the power system is divided into different stages; based on the resource allocation data of the current stage and the calculation results of the flexibility resources in the historical stage, the gain of the flexibility transformation resources of the thermal power units configured in the current stage, the peak shaving resource gain of the energy storage system, the consumption of the flexibility transformation resources of the thermal power units, and the resource consumption of the energy storage system are analyzed and calculated, comprehensively considering the resource gain and resource consumption of the thermal power units and the energy storage system; based on the consideration of the wind and light abandonment rate in the current stage, with the overall economic benefit of the current stage and the previous stage as the objective function of the current stage, and the calculation result of the flexibility resources in the current stage as the initial condition for the next stage for iterative optimization, the objective functions of all stages are obtained; under the predetermined constraint conditions, the objective functions of all stages are solved to obtain the flexibility resource calculation results of each stage, improving the accuracy and rationality of the calculation of the flexibility resources of the power system. This method is applicable to both the large-scale configuration of the energy storage system and the occasion with too high penetration rate of renewable energy; it explores the internal relationship between the resource allocation data of the power system and the calculation results of the flexibility resources, that is, the internal relationship between the resource allocation data of the power system and the newly added energy storage capacity, the newly added converter power capacity, and the newly added thermal power unit capacity; at the same time, using computer technology to automatically generate and display the calculation results of the flexibility resources improves the efficiency of the calculation of the flexibility resources of the power system.

[0054] The power system flexibility resource calculation method provided by the embodiments of the present disclosure takes into account the technical and economic changes in the energy storage and thermal power flexibility transformation processes, and can achieve the overall economic optimum in each planning stage.

[0055] According to an embodiment of the present disclosure, the resource allocation planning method provided by the present disclosure is more suitable for the optimal allocation of flexibility resources in a high-proportion renewable energy power system. The obtained flexibility resource calculation results can not only provide technical support for the flexibility resource planning decision-making of a high-proportion renewable energy power system, but also effectively implement the management requirements of "cost reduction and efficiency improvement" of the power grid company, and also provide theoretical and practical basis for power enterprises to further improve the planning refinement level and market competitiveness.

[0056] According to an embodiment of the present disclosure, solving M objective functions to obtain M flexibility resource calculation results corresponding to M stages includes: for the i-th stage, solving the objective function of the i-th stage to obtain the flexibility resource calculation result of the i-th stage, and obtaining M flexibility resource calculation results corresponding to M stages.

[0057] According to an embodiment of the present disclosure, the flexibility resource calculation result of the current stage can be obtained by solving according to the objective function of the current stage. For example, the calculation of the flexibility resources of the target power system can be divided into 5 stages according to the construction sequence, where the 5 stages correspond to 5 groups of sub-data. For the 4th stage, solve the objective function of the 4th stage to obtain the flexibility resource calculation result of the 4th stage, and obtain 5 flexibility resource calculation results corresponding to 5 stages.

[0058] According to an embodiment of the present disclosure, according to the method of claim 1, wherein solving M objective functions to obtain M flexibility resource calculation results corresponding to M stages includes: summing the M objective functions to obtain a total objective function; solving the total objective function to obtain M flexibility resource calculation results corresponding to M stages.

[0059] According to an embodiment of the present disclosure, the total objective function of the power system flexibility resource calculation can be calculated according to the objective functions of all stages, and the total objective function is solved to obtain the flexibility resource calculation results of each stage. For example, the calculation of the flexibility resources of the target power system can be divided into 5 stages according to the construction sequence, where the 5 stages correspond to 5 groups of sub-data. Sum the 5 objective functions to obtain a total objective function; solve the total objective function to obtain 5 flexibility resource calculation results corresponding to 5 stages.

[0060] According to an embodiment of the present disclosure, the flexibility transformation resource consumption of thermal power units in the i-th stage includes the flexibility transformation resource consumption of thermal power units in the i-th stage, the operation resource consumption of thermal power units in the i-th stage, and the environmental resource consumption of thermal power units in the i-th stage.

[0061] According to an embodiment of the present disclosure, the resource consumption of the flexibility transformation of thermal power units in the i-th stage can be measured by the annual value of the investment cost of the flexibility transformation of thermal power units in the i-th stage.

[0062] According to an embodiment of the present disclosure, the annual value of the investment cost of the flexibility transformation of thermal power units in the i-th stage, C the-inv,i can be calculated according to formula (1):

[0063] (1)

[0064] In the formula, C the-inv,i is the annual value of the investment cost of the flexibility transformation of thermal power units in the i-th stage; C b is the investment cost of the flexibility transformation of thermal power units per unit capacity; M i is the maximum number of units to be transformed in the i-th stage; P the-min0,n , P the-min,n are the minimum outputs of the n-th thermal power unit before and after transformation; r is the discount rate; L is the theoretical operation life of the transformed unit.

[0065] According to an embodiment of the present disclosure, the operation resource consumption of thermal power units in the i-th stage includes the operation fuel resource consumption in the i-th stage.

[0066] According to an embodiment of the present disclosure, the operation fuel resource consumption in the i-th stage is calculated according to the flexibility resource calculation results from the 1st stage to the (i - 1)-th stage and the i-th group of sub-data by using the following formula (2):

[0067] (2)

[0068] In the formula, C the-op,i is the average annual operation fuel cost, N i is the sum of the total number of units to be transformed from the 1st stage to the i-th stage, P the-n,t is the output power of the n-th thermal power flexibility transformation unit at time t, C(P the-n,t ) is the functional relationship between the output power of the n-th thermal power flexibility transformation unit at time t and its fuel resource consumption, and the functional relationship is established by using formula (3):

[0069] (3)

[0070] In the formula, a, b, and c are the quadratic term, linear term, and constant term coefficients of the coal consumption characteristic function of the thermal power unit respectively, c coal is the unit coal price, c oil is the unit fuel price, L lost is the rotor loss life consumed when deeply depressing the output of the unit, c unit is the purchase price of the unit rotor, Z oilis the average hourly oil injection volume when the unit conducts oil injection for peak shaving, P a,n is the minimum technical output during the conventional peak shaving process of the nth thermal power unit, P b,n is the minimum technical output during the non - oil - injection peak shaving process of the nth thermal power unit, P c,n is the minimum technical output during the oil - injection peak shaving process of the nth thermal power unit, P max,n is the maximum output of the nth thermal power unit.

[0071] According to the embodiments of the present disclosure, based on the regulation characteristics of the flexibility transformation of thermal power units, the peak shaving process can be divided into conventional peak shaving, non - oil - injection peak shaving, and oil - injection peak shaving. Among them, both non - oil - injection peak shaving and oil - injection peak shaving involve deeply depressing the unit output, which will cause rotor fatigue loss.

[0072] Figure 3 Schematically shows the peak shaving process of the flexibility transformation of thermal power units according to the embodiments of the present disclosure.

[0073] As Figure 3 shown, P a , P b , P c are respectively the minimum technical outputs of the three peak shaving processes, and P max is the maximum output of the unit.

[0074] According to the embodiments of the present disclosure, the operation resource consumption of thermal power units in the i - th stage may further include the operation and maintenance resource consumption in the i - th stage.

[0075] According to the embodiments of the present disclosure, the operation resource consumption of thermal power units in the i - th stage can be measured by the annual average operation cost of thermal power units in the i - th stage.

[0076] According to the embodiments of the present disclosure, the annual average operation cost C the-om,i of thermal power units in the i - th stage may include the annual average maintenance cost C the-ma,i and the annual average operation fuel cost C the-op,i , and is calculated using the following formula (4):

[0077] (4)

[0078] In the formula, C ma is the lifetime maintenance cost per unit capacity of the thermal power flexibility transformation unit, N i is the sum of the total number of retrofitted units in the i - th stage and the previous stages, P the-min0,n is the minimum output of the nth thermal power unit before retrofit, and C(P the-n,t ) is the functional relationship between the output power of the thermal power flexibility transformation unit and its fuel cost.

[0079] According to an embodiment of the present disclosure, the environmental resource consumption of a thermal power unit in the i-th stage can be measured by the annual value of the environmental cost of the thermal power unit in the i-th stage.

[0080] According to an embodiment of the present disclosure, the annual value of the environmental cost C of the thermal power unit in the i-th stage the-env,i can be calculated according to formula (5):

[0081] (5)

[0082] In the formula, K coal and K oil are the pollutant discharge fees generated per unit of coal and per unit of fuel burned less.

[0083] According to an embodiment of the present disclosure, the resource gain of the flexibility transformation of the thermal power unit in the i-th stage includes the power generation grid-connected resource gain of the flexible transformation unit of the thermal power in the i-th stage and the peak shaving resource gain of the flexible transformation unit in the i-th stage.

[0084] According to an embodiment of the present disclosure, the power generation grid-connected resource gain of the flexible transformation unit of the thermal power in the i-th stage can be calculated using formula (6):

[0085] (6)

[0086] In the formula, B the-gc,i is the power generation grid-connected resource gain of the flexible transformation unit of the thermal power in the i-th stage, and c pri,t is the on-grid peak-valley electricity price of the flexible transformation unit at time t.

[0087] According to an embodiment of the present disclosure, according to relevant regulations, a certain degree of compensation can be given for the deep peak shaving of the flexible transformation unit. The peak shaving resource gain of the flexible transformation unit in the i-th stage can be calculated using formula (7):

[0088] (7)

[0089] In the formula, B the-pr,i is the peak shaving resource gain of the flexible transformation unit in the i-th stage, and c the-gp,t is the deep peak shaving compensation electricity price of the flexible transformation unit at time t.

[0090] According to an embodiment of the present disclosure, the peak shaving resource gain of the energy storage system can be that the energy storage system uses the real-time peak-valley electricity price difference to charge the energy storage device when the electricity price is low and discharge the energy storage device when the electricity price is high, thereby reducing the electricity cost of users and obtaining corresponding benefits.

[0091] According to an embodiment of the present disclosure, the peak shaving resource gain of the energy storage system in the i-th stage is calculated using formula (8) based on the calculation results of the flexibility resources in the 1st stage to the (i-1)-th stage and the i-th set of sub-data:

[0092] (8)

[0093] In the formula, B ess-pr,i is the peak shaving resource gain of the energy storage system in the i-th stage, c pri,t is the on-grid peak-valley electricity price of the flexible retrofit unit at time t, P ess-c,m,t is the charging power of the energy storage system at time t in the m-th stage, P ess-d,m,t is the discharging power of the energy storage system at time t in the m-th stage, η c is the charging efficiency of the energy storage system, η d is the discharging efficiency of the energy storage system.

[0094] According to an embodiment of the present disclosure, the resource consumption of the energy storage system in the i-th stage is calculated using formula (9) based on the calculation results of the flexibility resources in the 1st stage to the (i-1)-th stage and the i-th set of sub-data:

[0095] (9)

[0096] In the formula, C ESS,i is the resource consumption of the energy storage system in the i-th stage, P pc,i is the newly added converter power in the i-th stage, K pc,i is the investment cost per unit power of the converter in the i-th stage, E ess,i is the newly added energy storage system capacity in the i-th stage, K ess,i is the investment cost per unit capacity of the energy storage system in the i-th stage, L ess is the theoretical operation life of the energy storage system, r is the discount rate, K om is the annual operation and maintenance cost per unit capacity of the energy storage system, E ess,m is the newly added energy storage system capacity in the m-th stage.

[0097] According to an embodiment of the present disclosure, the predetermined constraint conditions include the flexible retrofit constraint of thermal power units.

[0098] According to an embodiment of the present disclosure, the flexible retrofit constraint of thermal power units includes the retrofit depth constraint of thermal power units, the output constraint of thermal power units after flexible retrofit, and the ramp rate constraint of thermal power units after flexible retrofit.

[0099] According to an embodiment of the present disclosure, the retrofit depth constraint of thermal power units is determined using formula (10):

[0100] (10)

[0101] The output constraint after the flexibility transformation of the thermal power unit is determined using Equation (11):

[0102] (11)

[0103] The ramp constraint after the flexibility transformation of the thermal power unit is determined using Equation (12):

[0104] (12)

[0105] In the formula, P the-min0,n is the minimum output of the nth thermal power unit before transformation; P the-min,n is the minimum output of the nth thermal power unit after transformation; P the-gzmax,n is the upper limit of the flexibility transformation of the nth unit; P c,n is the minimum technical output during the oil - injection peak - shaving process of the nth thermal power unit; P max,n is the maximum output of the nth thermal power unit; R n0 is the ramp - rate limit of the nth unit before transformation; R n is the ramp - rate limit of the nth unit after transformation; x n is a 0 - 1 state variable, which is 1 when the flexibility transformation is carried out, otherwise 0; P the-n,t is the output power of the nth thermal power unit with flexibility transformation at time t.

[0106] According to the embodiments of the present disclosure, the predetermined constraint conditions further include the state - of - charge constraint of the energy storage system, the power - to - capacity constraint of the energy storage system, and the system power - balance constraint.

[0107] According to the embodiments of the present disclosure, to ensure the continuous operation of the energy storage system during the day, it is necessary to consider that its initial and final states are the same during the dispatching process.

[0108] According to the embodiments of the present disclosure, the state - of - charge constraint of the energy storage system can be determined using Equation (13):

[0109] (13)

[0110] In the formula, SOC is the state - of - charge of the energy storage system; SOC(0) and SOC(24) represent the state - of - charge at 0:00 and 24:00 every day, respectively

[0111] According to the embodiments of the present disclosure, the power - to - capacity constraint of the energy storage system can be determined using Equation (14):

[0112] (14)

[0113] In the formula, P ess-c,i,min 、P ess-c,i,max 、P ess-d,i,min 、Pess-d,i,max are the charging and discharging power limits of the energy storage system in the i-th stage; E ess,i,min and E ess,i,max are the capacity limits of the energy storage system in the i-th stage.

[0114] According to an embodiment of the present disclosure, the system power balance constraint can be determined using formula (15):

[0115] (15)

[0116] In the formula, P abon,i,t is the wind and light curtailment power at time t in the i-th stage; P Load,i,t is the grid-connected load power at time t in the planned year of the i-th stage; P Wind,i and P Solar,i are the predicted wind and light power generations at time t in the planned year of the i-th stage according to their construction timings.

[0117] According to an embodiment of the present disclosure, when performing resource allocation for the power system, it should be ensured that as much renewable energy power generation as possible is consumed. However, if the objective function only includes the renewable energy utilization rate, it will lead to the configuration of a large number of energy storage and thermal power retrofit units, resulting in a sharp increase in the peak shaving cost. On the contrary, if only the principle of optimal economic benefit is considered, it will lead to a large amount of wind and light curtailment. Therefore, to balance the benefits and the new energy utilization rate, the embodiments of the present disclosure consider the impact of wind and light curtailment and optimize with the goal of maximizing the total comprehensive benefit of the current stage and the previous stage.

[0118] According to an embodiment of the present disclosure, the objective function of the i-th stage is determined using formula (16) based on the flexibility retrofit resource gain of the thermal power units in the i-th stage, the peak shaving resource gain of the energy storage system in the i-th stage, the flexibility retrofit resource consumption of the thermal power units in the i-th stage, the resource consumption of the energy storage system in the i-th stage, the wind and light curtailment rate in the i-th stage, and the objective function of the (i - 1)-th stage:

[0119] (16)

[0120] In the formula, F i is the objective function of the i-th stage, F i-1 is the objective function of the (i - 1)-th stage, N i is the sum of the total number of retrofitted units from the 1st stage to the i-th stage, B the-gc,i is the power generation grid-connected resource gain of the flexible retrofit thermal power units in the i-th stage, B the-pr,i is the peak shaving resource gain of the flexible retrofit units in the i-th stage, B ess-pr,i is the peak shaving resource gain of the energy storage system in the i-th stage, C ESS,i is the resource consumption of the energy storage system in the i-th stage, C the-inv,iThe annual value of the investment cost for the flexibility retrofit of thermal power units in the $i$-th stage, $C$ the-om,i The average annual operating cost of thermal power units in the $i$-th stage, $C$ the-env,i The annual value of the environmental cost of thermal power units in the $i$-th stage, $P$ Wind,i,t The wind power generation power at time $t$ in the $i$-th stage planning year, $P$ Solar,i,t The photovoltaic power generation power at time $t$ in the $i$-th stage planning year, $P$ the-n,m,t The output power of the $n$-th thermal power unit with flexibility retrofit at time $t$ in the $m$-th stage, $P$ ess-c,m,t The charging power of the energy storage system at time $t$ in the $m$-th stage, $P$ ess-d,m,t The discharging power of the energy storage system at time $t$ in the $m$-th stage, $P$ Load,i,t The grid-connected load power at time $t$ in the $i$-th stage planning year.

[0121] According to an embodiment of the present disclosure, the denominator represents the curtailment rate of wind and photovoltaic power in the $i$-th stage.

[0122] Figure 4 Schematically shows a flowchart of a method for calculating flexibility resources of a power system according to another embodiment of the present disclosure.

[0123] As Figure 4 shown, the method includes operations S401 to S417.

[0124] In operation S401, according to the construction schedule of renewable energy, the optimal allocation of flexibility resources is divided into different stages, and relevant data required for the first stage is investigated according to Table 1.

[0125] Table 1

[0126]

[0127] In operation S402, the flexibility configuration result of the $(i - 1)$-th stage is used as a historical condition and substituted into the flexibility configuration of the $i$-th stage, and relevant data required in the $i$-th stage is investigated.

[0128] In operation S403, the initial value of the newly added flexibility retrofitted thermal power units in the $i$-th stage is set.

[0129] In operation S404, the benefits and costs of the thermal power units after flexibility retrofit in the $i$-th stage are calculated.

[0130] In operation S405, the initial value of the newly added energy storage in the $i$-th stage is set.

[0131] In operation S406, the benefits and costs of the energy storage configuration in the $i$-th stage are calculated.

[0132] In operation S407, according to the current flexibility configuration, the curtailment rate of wind and photovoltaic power is calculated.

[0133] In operation S408, under the set constraint conditions, the objective function is obtained.

[0134] In operation S409, it is judged whether it is the maximum value of the objective function under the current flexible retrofitted thermal power unit.

[0135] If the judgment in operation S409 is yes, operation S410 is executed.

[0136] In operation S410, the maximum value of the objective function is updated, and the calculation result of the current flexible resource is recorded.

[0137] If the judgment in operation S409 is no, operation S411 is executed.

[0138] In operation S411, according to the particle swarm algorithm, the energy storage configuration variable is corrected. Repeat operations S406 - S409.

[0139] In operation S412, it is judged whether all possibilities of the flexible retrofit of thermal power units are traversed.

[0140] If the judgment in operation S412 is yes, operation S413 is executed.

[0141] In operation S413, the calculation result of the flexible resource in the i-th stage is output.

[0142] If the judgment in operation S412 is yes, operation S414 is executed.

[0143] In operation S414, the variable of the newly added flexible retrofit unit in the i-th stage is corrected. Repeat operations S404 - S412.

[0144] In operation S415, it is judged whether the calculation of flexible resources for all stages is planned.

[0145] If the judgment in operation S415 is yes, operation S416 is executed.

[0146] In operation S416, the optimal configuration scheme of flexible resources considering the construction time sequence of renewable energy is output.

[0147] If the judgment in operation S415 is yes, operation S417 is executed.

[0148] In operation S417, the stage is updated to make i = i + 1. Repeat operations S402 - S415.

[0149] According to an embodiment of the present disclosure, taking a renewable energy power generation base mainly based on wind power in a certain province as an example, the method herein is analyzed and verified. First, according to the relevant plan, the commissioning plan of this base is 180 MW in the first phase (2023, the first stage), 300 MW in the second phase (2025, the second stage), and 320 MW in the third phase (2028, the third stage), with a total of 800 MW. The power generation power after commissioning in the long-term year (2028) is as Figure 5 shown. The power generation power in each previous stage is reduced according to the installed capacity ratio.

[0150] According to an embodiment of the present disclosure, there are 16 thermal power units in this area. Among them, 12 units support flexibility transformation, and their capacities are 4 units of 50 MW and 8 units of 30 MW respectively. In addition, due to land use reasons in this area, the land for the energy storage system is at most 1200 MWh, so the maximum scale of energy storage configuration is 1200 MWh. The grid-connected load in this area in 2022 is as Figure 6 shown, and the grid-connected load in each predicted year of each stage can be predicted at an increasing rate of 8% per year.

[0151] According to an embodiment of the present disclosure, the research results of the parameters required for the optimal allocation of the remaining flexibility resources are shown in Table 2

[0152] Table 2

[0153]

[0154] According to an embodiment of the present disclosure, the particle swarm optimization algorithm can be used for correction, the population number is set to 30, and the maximum number of iterations is 100 times.

[0155] According to an embodiment of the present disclosure, to reflect the effectiveness of the power system flexibility resource calculation method considering the construction time sequence of renewable energy proposed in the present disclosure, the present disclosure sets 3 scenarios for comparison.

[0156] Scenario 1: Traditional static configuration method, that is, perform "one-step" optimization on the flexibility resource configuration.

[0157] Scenario 2: Multi-stage configuration method considering the construction time sequence of renewable energy, that is, use the static configuration method to optimize stage by stage considering the construction time sequence of renewable energy, and ensure the optimal in each stage.

[0158] Scenario 3: The flexibility resource dynamic configuration method considering the construction time sequence of renewable energy proposed in the present invention, that is, on the basis of Scenario 2, consider the influence of the configuration result of the previous stage on the future configuration, and ensure the sum of the target values of all stages is optimal.

[0159] According to an embodiment of the present disclosure, the static configuration result of the flexibility resources in Scenario 1 is shown in Table 3.

[0160] Table 3

[0161]

[0162] According to the embodiments of the present disclosure, the multi-stage configuration results of the flexibility resources in Scenario 2 are shown in Table 4 as follows.

[0163] Table 4

[0164]

[0165] According to the embodiments of the present disclosure, the costs of the dynamic configuration of the flexibility resources in Scenario 3 are shown in Table 5, and the results of the benefit indicators and the dynamic configuration are shown in Table 6.

[0166] Table 5

[0167]

[0168] Table 6

[0169]

[0170] According to the embodiments of the present disclosure, by comparing Tables 3, 4, and 6, it can be found that the profit of the static configuration result of the flexibility resources is the lowest and the objective function value is the worst. This is because the configuration of the energy storage and the retrofit of the thermal power units are both carried out in the starting year, and the technical and economic changes during the process of the energy storage and the flexibility retrofit of the thermal power (such as the decrease in the configuration cost, the innovation of the energy storage technology, etc.) cannot be considered, resulting in relatively high configuration costs for each item and the worst economy. The methods of Scenarios 2 and 3 can take into account the changes in the economic parameters of various resources during the configuration process. Therefore, generally speaking, the methods of Scenarios 2 and 3 are superior to the static configuration method of Scenario 1.

[0171] According to the embodiments of the present disclosure, the multi-stage configuration method of Scenario 2 and the dynamic configuration method of Scenario 3 are mainly compared below.

[0172] Among them, the multi-stage configuration method of Scenario 2 mainly focuses on the technical and economic optimization of each stage, while the dynamic configuration method of Scenario 3 pursues the overall technical and economic optimization during the entire planning period on the basis of the multi-stage configuration method. Because the dynamic configuration method takes into account the influence of the flexibility configuration result of the previous stage on the subsequent configuration, it can perform global optimization from an overall perspective.

[0173] According to the embodiments of the present disclosure, taking the second stage of Tables 4 and 6 as an example for explanation, when configuring flexibility resources with the same capacity in the second stage, the economy of the flexibility transformation of thermal power units is superior to that of energy storage configuration and curtailment of wind and solar power. Therefore, the multi-stage configuration method in Scenario 2 mainly focuses on the flexibility transformation of thermal power units. Compared with Scenario 3, three more thermal power units are flexibly transformed, and its profit is also 655.8 million yuan more than that of Scenario 3. However, in the third stage, due to the improvement of energy storage technology and the significant reduction of energy storage costs, the economy of energy storage is superior to that of the flexibility transformation of thermal power units. In addition, thermal power units have their peak shaving limits. In some occasions with high penetration of renewable energy, the flexibly transformed units cannot perform peak shaving normally, resulting in the phenomenon of curtailment of wind and solar power. The above two reasons lead to a significant reduction in the effectiveness of the previous thermal power transformation. Therefore, in the third stage, when overall measuring the total benefits during the construction period of renewable energy, the dynamic configuration method in Scenario 3 is superior to the multi-stage configuration method in Scenario 2.

[0174] According to the embodiments of the present disclosure, through the comparative analysis of the examples of the three scenarios of the present disclosure, the dynamic configuration method of flexibility resources considering the construction time sequence of renewable energy proposed by the present disclosure takes into account the construction time sequence of renewable energy and the improvement of the technology and economy of flexibility resources, and is more applicable to the optimal configuration of flexibility resources in a high-proportion renewable energy power system than the existing methods. Its final configuration scheme is shown in Table 6.

[0175] Figure 7 The block diagram of the power system flexibility resource calculation device according to the embodiments of the present disclosure is schematically shown.

[0176] As Figure 7 shown, the power system flexibility resource calculation device 700 includes an acquisition module 710, a division module 720, a calculation module 730, an establishment module 740, a solution module 750, and a display module 760.

[0177] The acquisition module 710 is configured to, in response to receiving a flexibility resource calculation request for a target power system, obtain resource configuration data of the target power system from a database, where the resource configuration data includes M groups of sub-data divided according to the construction time sequence of the target power system, and M is a positive integer greater than or equal to 2.

[0178] The division module 720 is configured to divide the flexibility resource calculation of the target power system into M stages according to the construction time sequence, where the M stages correspond to the M groups of sub-data.

[0179] A calculation module 730, configured to calculate, for the i-th stage, a gain of flexibility transformation resources of a thermal power unit in the i-th stage, a peak shaving resource gain of an energy storage system in the i-th stage, a consumption of flexibility transformation resources of a thermal power unit in the i-th stage, a resource consumption of an energy storage system in the i-th stage, and a wind and light abandonment rate in the i-th stage according to the calculation results of flexibility resources in the first stage to the (i-1)-th stage and the i-th group of sub-data, where i is a positive integer and 1≤i≤M.

[0180] A building module 740, configured to build an objective function for the i-th stage and obtain M objective functions according to the gain of flexibility transformation resources of a thermal power unit in the i-th stage, the peak shaving resource gain of an energy storage system in the i-th stage, the consumption of flexibility transformation resources of a thermal power unit in the i-th stage, the resource consumption of an energy storage system in the i-th stage, the wind and light abandonment rate in the i-th stage, and the objective function in the (i-1)-th stage.

[0181] A solving module 750, configured to solve the M objective functions under predetermined constraint conditions to obtain M calculation results of flexibility resources corresponding to the M stages, where the calculation results of flexibility resources include a newly added energy storage capacity, a newly added converter power capacity, and a newly added thermal power unit capacity.

[0182] A display module 760, configured to visually display the M calculation results of flexibility resources.

[0183] Any multiple of the modules, sub-modules, units, and sub-units according to the embodiments of the present disclosure, or at least part of the functions of any of them can be implemented in one module. Any one or more of the modules, sub-modules, units, and sub-units according to the embodiments of the present disclosure can be split into multiple modules for implementation. Any one or more of the modules, sub-modules, units, and sub-units according to the embodiments of the present disclosure can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on chip, a system on substrate, a system on package, an application specific integrated circuit (ASIC), or can be implemented by hardware or firmware in any other reasonable way of integrating or packaging circuits, or can be implemented in any one of the three implementation manners of software, hardware, and firmware, or in a suitable combination of any several of them. Alternatively, one or more of the modules, sub-modules, units, and sub-units according to the embodiments of the present disclosure can be at least partially implemented as a computer program module, and when the computer program module runs, it can execute the corresponding functions.

[0184] For example, any number of the acquisition module 710, the division module 720, the calculation module 730, the establishment module 740, the solution module 750, and the display module 760 can be combined and implemented in one module / unit / sub-unit, or any one of the modules / units / sub-units can be split into multiple modules / units / sub-units. Alternatively, at least part of the functions of one or more of these modules / units / sub-units can be combined with at least part of the functions of other modules / units / sub-units and implemented in one module / unit / sub-unit. According to an embodiment of the present disclosure, at least one of the acquisition module 710, the division module 720, the calculation module 730, the establishment module 740, the solution module 750, and the display module 760 can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on chip, a system on substrate, a system on package, an application specific integrated circuit (ASIC), or can be implemented by any other reasonable means such as hardware or firmware through circuit integration or packaging, or can be implemented in any one of the three implementation manners of software, hardware, and firmware, or in an appropriate combination of any several of them. Alternatively, at least one of the acquisition module 710, the division module 720, the calculation module 730, the establishment module 740, the solution module 750, and the display module 760 can be at least partially implemented as a computer program module, and when the computer program module is run, the corresponding functions can be executed.

[0185] It should be noted that the data processing system part in the embodiments of the present disclosure corresponds to the data processing method part in the embodiments of the present disclosure. For the description of the data processing system part, please refer to the data processing method part specifically, and details are not described herein again.

[0186] Figure 8 A block diagram of an electronic device suitable for implementing the method described above according to an embodiment of the present disclosure is schematically shown. Figure 8 The shown electronic device is only an example and should not impose any limitation on the functions and usage scope of the embodiments of the present disclosure.

[0187] As Figure 8As shown, the electronic device 800 according to an embodiment of the present disclosure includes a processor 801, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 802 or a program loaded from a storage section 808 into a random access memory (RAM) 803. The processor 801 can include, for example, a general microprocessor (such as a CPU), an instruction set processor, and / or a related chipset, and / or a dedicated microprocessor (such as an application specific integrated circuit (ASIC)), and so on. The processor 801 can also include on-board memory for caching purposes. The processor 801 can include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of the present disclosure.

[0188] In the RAM 803, various programs and data required for the operation of the electronic device 800 are stored. The processor 801, the ROM 802, and the RAM 803 are connected to each other via a bus 804. The processor 801 performs various operations of the method flow according to an embodiment of the present disclosure by executing the program in the ROM 802 and / or the RAM 803. It should be noted that the program can also be stored in one or more memories other than the ROM 802 and the RAM 803. The processor 801 can also perform various operations of the method flow according to an embodiment of the present disclosure by executing the program stored in the one or more memories.

[0189] According to an embodiment of the present disclosure, the electronic device 800 can further include an input / output (I / O) interface 805, and the input / output (I / O) interface 805 is also connected to the bus 804. The electronic device 800 can further include one or more of the following components connected to the input / output (I / O) interface 805: an input section 806 including a keyboard, a mouse, etc.; an output section 807 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc. and a speaker, etc.; a storage section 808 including a hard disk, etc.; and a communication section 809 including a network interface card such as a LAN card, a modem, etc. The communication section 809 performs communication processing via a network such as the Internet. A drive 810 is also connected to the input / output (I / O) interface 805 as needed. A removable medium 811, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 810 as needed so that a computer program read from it can be installed into the storage section 808 as needed.

[0190] According to an embodiment of the present disclosure, the method flow according to the embodiment of the present disclosure can be implemented as a computer software program. For example, an embodiment of the present disclosure includes a computer program product, which includes a computer program carried on a computer-readable storage medium, and the computer program includes program codes for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through the communication part 809, and / or installed from the removable medium 811. When the computer program is executed by the processor 801, the above functions defined in the system of the embodiment of the present disclosure are executed. According to an embodiment of the present disclosure, the above-described systems, devices, apparatuses, modules, units, etc. can be implemented by computer program modules.

[0191] The present disclosure also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiment; or may exist separately without being assembled into the device / apparatus / system. The above computer-readable storage medium carries one or more programs, and when the above one or more programs are executed, the method according to the embodiment of the present disclosure is implemented.

[0192] According to an embodiment of the present disclosure, the computer-readable storage medium may be a non-volatile computer-readable storage medium. For example, it may include but is not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the above. In the present disclosure, the computer-readable storage medium may be any tangible medium that contains or stores a program, and the program can be used by or in combination with an instruction execution system, apparatus, or device.

[0193] For example, according to an embodiment of the present disclosure, the computer-readable storage medium may include the above-described ROM 802 and / or RAM 803 and / or one or more memories other than ROM 802 and RAM 803.

[0194] An embodiment of the present disclosure also includes a computer program product, which includes a computer program, and the computer program includes program codes for executing the method provided by the embodiment of the present disclosure. When the computer program product runs on an electronic device, the program codes are used to enable the electronic device to implement the power system flexibility resource calculation method provided by the embodiment of the present disclosure.

[0195] When the computer program is executed by the processor 801, the above functions defined in the system / apparatus of the embodiment of the present disclosure are executed. According to an embodiment of the present disclosure, the above-described systems, apparatuses, modules, units, etc. can be implemented by computer program modules.

[0196] In one embodiment, the computer program may rely on tangible storage media such as optical storage devices, magnetic storage devices, etc. In another embodiment, the computer program may also be transmitted and distributed in the form of signals on a network medium, and downloaded and installed through the communication part 809, and / or installed from the removable medium 811. The program code included in the computer program may be transmitted using any suitable network medium, including but not limited to: wireless, wired, etc., or any suitable combination of the above.

[0197] According to an embodiment of the present disclosure, the program code for executing the computer program provided by the embodiments of the present disclosure may be written in any combination of one or more programming languages. Specifically, these computing programs may be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages include but are not limited to, such as Java, C++, Python, the "C" language, or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user device, partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device may be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computing device (e.g., by connecting through the Internet using an Internet service provider).

[0198] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagram may represent a module, a segment of a program, or a portion of code that contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions noted in the blocks may occur in a different order than noted in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, or they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram or flowchart, and combinations of blocks in the block diagram or flowchart, can be implemented by a dedicated hardware-based system that performs the specified functions or operations, or by a combination of dedicated hardware and computer instructions. Those skilled in the art will appreciate that the features described in the various embodiments of the present disclosure can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in the present disclosure. In particular, without departing from the spirit and teachings of the present disclosure, the features described in the various embodiments of the present disclosure can be combined and / or combined in various ways. All such combinations and / or combinations fall within the scope of the present disclosure.

[0199] The embodiments of the present disclosure have been described above. However, these embodiments are for illustrative purposes only and are not intended to limit the scope of the present disclosure. Although the embodiments have been described separately above, this does not mean that the measures in the various embodiments cannot be used advantageously in combination. Without departing from the scope of the present disclosure, those skilled in the art can make various substitutions and modifications, and all such substitutions and modifications should fall within the scope of the present disclosure.

Claims

1. A method for calculating power system flexibility resources, comprising: In response to receiving a flexibility resource calculation request of a target power system, obtaining resource configuration data of the target power system from a database, wherein the resource configuration data includes M groups of sub-data divided according to a construction sequence of the target power system, where M is a positive integer ≥ 2; Dividing the target power system flexibility resource calculation into M stages according to the construction sequence, wherein the M stages correspond to M groups of sub-data; For the i-th stage, according to the calculation results of flexibility resources from the 1st stage to the i-1th stage and the i-th group of sub-data, the flexibility transformation resource gain of the thermal power unit in the i-th stage, the peak-shaving resource gain of the energy storage system in the i-th stage, the flexibility transformation resource consumption of the thermal power unit in the i-th stage, the resource consumption of the energy storage system in the i-th stage and the wind and solar power abandonment rate in the i-th stage are calculated, where i is a positive integer, 1≤i≤M; According to the resource gain of the flexibility transformation of the thermal power unit in the i-th stage, the peak-shaving resource gain of the energy storage system in the i-th stage, the resource consumption of the flexibility transformation of the thermal power unit in the i-th stage, the resource consumption of the energy storage system in the i-th stage, the wind and solar power abandonment rate in the i-th stage and the objective function of the i-1-th stage, an objective function for the i-th stage is established to obtain M objective functions; and Under predetermined constraints, solving the M objective functions to obtain M flexibility resource calculation results corresponding to the M stages, wherein the flexibility resource calculation results include newly added energy storage capacity, newly added converter power capacity, and newly added thermal power unit capacity; Visually display the calculation results of the M flexibility resources; The peak-shaving resource gain of the energy storage system in the i-th stage is calculated based on the calculation results of the flexibility resources from the 1st stage to the i-1th stage and the i-th group of sub-data using formula (8): （8） In the formula, B ess-pr,i is the peak-shaving resource gain of the energy storage system in stage i, c pri,t To flexibly transform the peak and valley electricity price of the unit at time t, P ess-c,m,t is the charging power of the energy storage system at time t in the mth stage, P ess-d,m,t is the energy storage system discharge power at time t in the mth stage, η c is the charging efficiency of the energy storage system, η d is the discharge efficiency of the energy storage system; The energy storage system resource consumption in the i-th stage is calculated based on the flexibility resource calculation results from the 1st stage to the i-1th stage and the i-th group of sub-data using formula (9): （9） In the formula, C ESS,i is the energy storage system resource consumption in the i-th stage, P pc,i is the converter power added in stage i, K pc,i is the unit power converter investment cost in stage i, E ess,i is the newly added energy storage system capacity in stage i, K ess,i is the investment cost of energy storage system per unit capacity in stage i, L ess is the theoretical operating life of the energy storage system, r is the discount rate, K om is the annual operation and maintenance cost per unit capacity of the energy storage system, E ess,m It is the capacity of the energy storage system newly added in the mth phase.

2. The method according to claim 1, wherein: Solving the M objective functions to obtain M flexibility resource calculation results corresponding to the M stages includes: For the i-th stage, the objective function of the i-th stage is solved to obtain the flexibility resource calculation result of the i-th stage, and M flexibility resource calculation results corresponding to the M stages are obtained.

3. The method according to claim 1, wherein: Solving the M objective functions to obtain M flexibility resource calculation results corresponding to the M stages includes: Adding the M objective functions to obtain a total objective function; The overall objective function is solved to obtain M flexibility resource calculation results corresponding to the M stages.

4. The method according to any one of claims 1 to 3, wherein The resource consumption for flexibility transformation of thermal power units in the i-th phase includes the resource consumption for flexibility transformation of thermal power units in the i-th phase, the resource consumption for operation of thermal power units in the i-th phase and the environmental resource consumption of thermal power units in the i-th phase.

5. The method according to claim 4, wherein: The operation resource consumption of the thermal power unit in the i-th stage includes the operation fuel resource consumption in the i-th stage; The operating fuel resource consumption of the i-th stage is calculated based on the flexibility resource calculation results from the 1st stage to the i-1th stage and the i-th group of sub-data using the following formula (2): （2） In the formula, C the-op,i is the average annual operating fuel cost, N i is the sum of the total number of modified units from the first stage to the i-th stage, P the-n,t is the output power of the nth thermal power unit with flexible transformation at time t, C(P the-n,t ) is the functional relationship between the output power of the nth thermal power flexible transformation unit and its fuel resource consumption at time t, and the functional relationship is established using formula (3): （3） Where a, b, and c are the coefficients of the quadratic term, linear term, and constant term of the thermal power unit consumption characteristic function, respectively. coal is the unit coal price, c oil is the unit fuel price, L lost is the rotor loss life consumed when the unit output is deeply reduced, c unit is the purchase price of the unit rotor, Z oil is the average oil injection per hour during the peak load regulation of the unit, P a,n is the minimum technical output of the nth thermal power unit in the conventional peak load regulation process, P b,n is the minimum technical output of the nth thermal power unit during the non-oil peak load regulation process, P c,n is the minimum technical output of the nth thermal power unit during the oil-loading peak load regulation process, P max,n It is the maximum output of the nth thermal power unit.

6. The method according to any one of claims 1 to 3, wherein: The objective function of the i-th stage is determined by formula (16) based on the resource gain of the flexibility transformation of the thermal power units in the i-th stage, the peak-shaving resource gain of the energy storage system in the i-th stage, the resource consumption of the flexibility transformation of the thermal power units in the i-th stage, the resource consumption of the energy storage system in the i-th stage, the wind and solar power abandonment rate in the i-th stage and the objective function of the i-1-th stage: （16） In the formula, F i is the objective function of the i-th stage, F i-1 is the objective function of the i-1th stage, N i is the sum of the total number of modified units from the first stage to the i-th stage, B the-gc,i is the resource gain of the thermal power flexibility transformation unit in the i-th phase, B the-pr,i is the peak load resource gain of the flexibility transformation unit in stage i, B ess-pr,i is the peak load resource gain of the energy storage system in stage i, C ESS,i is the energy storage system resource consumption in the i-th stage, C the-inv,i is the annual value of the investment cost of the flexibility transformation of thermal power units in the i-th phase, C the-om,i is the average annual operating cost of the thermal power unit in the i-th phase, C the-env,i is the annual environmental cost of the thermal power unit in the i-th phase, P Wind,i,t is the wind power generation capacity at time t in the planning year of stage i, P Solar,i,t is the solar power generation power at time t in the planning year of stage i, P the-n,m,t is the output power of the nth thermal power unit with flexible transformation at time t in the mth stage, P ess-c,m,t is the charging power of the energy storage system at time t in the mth stage, P ess-d,m,t is the energy storage system discharge power at time t in the mth stage, P Load,i,t is the grid-connected load power at time t in the planning year of stage i.

7. The method according to any one of claims 1 to 3, wherein: The predetermined constraint conditions include the flexibility transformation constraint of the thermal power unit, and the flexibility transformation constraint of the thermal power unit includes the transformation depth constraint of the thermal power unit, the output constraint after the flexibility transformation of the thermal power unit, and the climbing constraint after the flexibility transformation of the thermal power unit; The transformation depth constraint of the thermal power unit is determined by formula (10): （10） The output constraint of the thermal power unit after flexibility transformation is determined by formula (11): （11） The ramp constraint after the flexibility transformation of the thermal power unit is determined by formula (12): （12） Where P the-min0,n is the minimum output of the nth thermal power unit before transformation; P the-min,n is the minimum output of the nth thermal power unit after transformation; P the-gzmax,n is the upper limit of the flexibility transformation of the nth unit; P c,n is the minimum technical output of the nth thermal power unit during the oil-loading peak-shaving process; P max,n is the maximum output of the nth thermal power unit; R n0 is the ramp rate limit of the nth unit before transformation; R n is the ramp rate limit of the nth unit after transformation; x n It is a 0-1 state variable, which is 1 when flexibility is modified, otherwise it is 0; P the-n,t The output power of the nth thermal power unit with flexibility transformation at time t.

8. The method according to claim 7, wherein: The predetermined constraints also include energy storage system state of charge constraints, energy storage system power and capacity constraints, and system power balance constraints.

9. The method according to any one of claims 1 to 3, wherein: The resource gain of the flexibility transformation of the thermal power units in the i-th stage includes the resource gain of the online power generation of the thermal power flexibility transformation units in the i-th stage and the peak-shaving resource gain of the flexibility transformation units in the i-th stage.

Images