Robust planning and operation method, device, equipment and medium for combined energy storage on the grid side

Abstract

The present invention relates to the technical field of multi-energy complementary power generation, and in particular to a method, device, equipment and medium for robust planning and operation of combined energy storage on the grid side. The method includes: considering various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, establishing a capacity configuration model of combined energy storage; taking the lowest system operation cost as the goal, considering the operation constraints and line flow constraints of cascade pumped storage and electrochemical energy storage, establishing a system operation optimization model; establishing a combined energy storage distribution robust planning-operation double-layer model including the construction of wind power and photovoltaic output probability distribution uncertainty and linearization model conversion; using a column and constraint generation algorithm to solve the double-layer model to obtain a combined energy storage planning and operation method. In the present invention, the capacity of combined energy storage under the coordination of source, grid, load and storage can be determined, and the system operation can be optimized, the amount of abandoned electricity can be reduced, the peak-shaving pressure of thermal power can be alleviated, and the frequent start and stop of thermal power and hydropower units can be reduced, and the result is robust.

Description

Technical Field

[0001] The present invention relates to the technical field of multi-energy complementary power generation, and in particular to a method, device, equipment and medium for robust planning and operation of combined energy storage on a power grid side. Background Art

[0002] Under the goal of carbon peak and carbon neutrality, wind power and photovoltaic power generation have been continuously developed. However, wind power and photovoltaic power generation are greatly affected by environmental factors, and it is difficult to accurately predict their power generation. Their uncertainty poses a challenge to the safe and stable operation of the power grid. In this case, the development of flexible resources has become an effective and necessary way, and pumped storage power stations and electrochemical energy storage have attracted attention from the industry.

[0003] However, the pumped-storage power stations currently under construction and planning have problems such as insufficient unit regulation capacity and large investment scale, making it difficult to meet the needs of new energy consumption and grid regulation. On the other hand, flexible reconstruction and expansion of pumped-storage units based on the original cascade hydropower upper and lower reservoirs to form cascade pumped-storage power stations can improve the system regulation capacity in the short term.

[0004] However, there are few studies on cascade pumped storage in the existing technology, and most of them focus on the operation field. Cascade pumped storage has not been combined with other energy storage such as electrochemical energy storage to form combined energy storage. Combined energy storage can combine the advantages of different energy storages, better absorb wind power and photovoltaic power generation, and improve the safety and stability of the power grid. Therefore, the planning and operation of combined energy storage is an issue that needs to be urgently addressed.

[0005] The information disclosed in this background technology section is only intended to deepen the understanding of the overall background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0006] The present invention provides a method, device, equipment and medium for robust planning and operation of combined energy storage on the grid side, thereby effectively solving the problems in the background technology.

[0007] In order to achieve the above object, the technical solution adopted by the present invention is: a grid-side combined energy storage robust planning and operation method, comprising the following steps:

[0008] Considering the various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, a capacity configuration model for combined energy storage is established;

[0009] With the goal of minimizing system operation cost, the operation constraints of cascade pumped storage and electrochemical energy storage, as well as line flow constraints, are considered to establish a system operation optimization model.

[0010] Establish a combined energy storage distributed robust planning-operation two-layer model that includes constructing the uncertainty of wind power and photovoltaic output probability distribution and linearizing model transformation;

[0011] The double-layer model is solved by using a column and constraint generation algorithm to obtain a combined energy storage planning and operation method.

[0012] Furthermore, the capacity configuration model of the combined energy storage includes:

[0013] minC=C in +C om -C re ;

[0014] in:

[0015]

[0016] C om =k om,ps I ps +k om,e E bat +k om,p P bat ;

[0017]

[0018] Where: C in , C om and C re are the annual equivalent investment cost, annual equivalent operation and maintenance cost, and annual equivalent residual value cost respectively; γ is the discount rate; Y ps and Y bat are the planned operating years of cascade pumped storage and electrochemical energy storage, respectively; k ps , k e and k p are the cost coefficients of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; I ps 、E bat and P bat are the capacity of cascade pumped storage, the capacity and power of electrochemical energy storage respectively; k om,ps , k om,e and k om,p are the annual operation and maintenance costs of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; C in,ps and C in,bat are the investment costs of cascade pumped storage and electrochemical energy storage respectively; δ ps and δ bat They are the residual value rates of cascade pumped storage and electrochemical energy storage, respectively.

[0019] Furthermore, the capacity configuration model includes the capacity and power constraints of cascade pumped storage and electrochemical energy storage:

[0020]

[0021] Where: and are the upper and lower limits of cascade pumped storage capacity respectively; and are the upper and lower limits of electrochemical energy storage capacity, respectively; and are the upper and lower limits of electrochemical energy storage power respectively.

[0022] Furthermore, the operation optimization model includes:

[0023] minC op =C th +C start +C cur ;

[0024] in:

[0025]

[0026]

[0027]

[0028] Where: C th , C start and C cur are the cost of thermal power units, the startup cost and the power abandonment cost of hydropower and pumped storage units respectively; N T 、N G and N H are the operating time, the number of thermal power units and the number of cascade hydropower stations respectively; a, b and c are the coefficients of the quadratic term, linear term and constant term of the thermal power coal consumption cost respectively; Provide power for thermal power units; and are the start and stop operation variables of thermal power units respectively; c G,st and c G,dn are the start-up and shutdown operation costs of thermal power units respectively; T is the carbon emission cost; Q and Q1 are the carbon emission and carbon emission quota respectively; and are the startup operation variables of the hydropower unit, pumped storage unit generating conditions and pumping conditions of the pumped storage unit respectively; c H,su 、c PSG,su and c PSP,su They are the startup costs of the hydropower unit, pumped storage unit generating conditions and pumped storage unit pumping conditions respectively; and They are the predicted output and actual output of photovoltaic power generation respectively; and are the predicted output and actual output of wind power respectively; c PV and c W are the penalty costs for curtailment of photovoltaic and wind power respectively.

[0029] Furthermore, the operation constraints and line flow constraints of cascade pumped storage and electrochemical energy storage are considered, including:

[0030] Power balance constraints:

[0031]

[0032] P t PS =P t PSG -P t PSP ;

[0033] P t bat =P t dis -P t cha ;

[0034] Where: P t PS Output of pumped storage unit; P t PSG and P t PSP are the power generation output and pumping output of the pumped storage unit respectively; P t bat Power for electrochemical energy storage; P t dis and P t cha are the discharge power and charging power of electrochemical energy storage respectively; is the load power;

[0035] Electrochemical energy storage operation constraints:

[0036]

[0037]

[0038]

[0039] Where: and are 0-1 variables for the electrochemical energy storage in the discharge state and the charge state respectively; E bat,min and E bat,max are the minimum and maximum capacity of electrochemical energy storage, respectively; It is the amount of electricity when electrochemical energy storage is working; hcha and h dis are charging power and discharging power respectively;

[0040] Hydropower unit output, start and stop, flow and storage capacity constraints:

[0041]

[0042]

[0043]

[0044]

[0045] Where: ρ is the density of water; g is the acceleration due to gravity; h i H is the power generation efficiency of the hydropower unit; H i is the average water head; is the state variable of the hydropower unit; and are the minimum and maximum outputs of the hydropower unit respectively; I i,t , and They are the generating flow of hydropower units, abandoned water flow, discharged flow, interval flow, pumped storage working flow, pumped water flow of pumped storage units and generating flow; and They are the start and stop operation variables of the hydropower unit respectively; and They are the minimum flow rate and maximum flow rate of the hydropower unit for power generation respectively; is the minimum discarded water flow; and are the marker variables of the upper reservoir and the lower reservoir of the pumped storage unit respectively; V i,t is the storage capacity; V i min and V i max are the minimum storage capacity and the maximum storage capacity respectively;

[0046] Output and start / stop constraints of pumped storage units:

[0047]

[0048]

[0049] Where: h is the average water head of the pumped storage unit; h PSG and h PSP are the power generation efficiency and pumping efficiency of the pumped storage unit respectively; and are the power generation and pumping state variables of the pumped storage unit respectively; α PSG and α PSP are the minimum coefficients of the generating power and pumping power of the pumped storage unit relative to the configured capacity, respectively; and are the stop operation variables of the pumped storage unit for power generation and water pumping, respectively;

[0050] Output, start / stop and ramp-up constraints of thermal power units:

[0051]

[0052]

[0053]

[0054] Where: and are the minimum and maximum output of thermal power units respectively; is the state variable of the thermal power unit; and It is the up and down climbing power of the thermal power unit; and are the minimum time that the thermal power unit is in the on and off state respectively;

[0055] Wind power and photovoltaic operation constraints:

[0056]

[0057] Line flow constraints:

[0058]

[0059] Where: G n-g , G n-i , G n-ps , G n-w , G n-p and G n-k are the power transfer factors of thermal power unit g, cascade hydropower station i, pumped storage unit, wind turbine unit w, photovoltaic generator unit p and load node k to line n respectively; Maximum line power limit.

[0060] Furthermore, the establishment includes constructing a combined energy storage distribution robust planning-operation two-layer model that includes building wind power and photovoltaic output probability distribution uncertainty and linearization model conversion, including:

[0061] Using clustering methods to obtain the empirical probability of typical scenarios from historical data A comprehensive norm including 1-norm and ∞-norm is used to construct the uncertainty set of wind power and photovoltaic output probability distribution;

[0062]

[0063] Where: Ψ and Ψ ∞ They are the uncertainty sets of probability distribution established based on 1-norm and ∞-norm respectively; and They are the allowable deviation ranges of probability distribution under 1-norm and ∞-norm respectively; and As shown below:

[0064]

[0065] Where: M is the total amount of historical data; α1 and α ∞ are the uncertainty confidence levels of the 1-norm and ∞-norm respectively;

[0066] The uncertainty set of wind power and photovoltaic output probability distribution is constructed as follows:

[0067]

[0068] The probability distribution uncertainty set is transformed into a linear model:

[0069]

[0070] A grid-side combined energy storage distributed robust planning-operation two-layer model considering the coordination of source, grid, load and storage is established. The outer layer of the model is the planning layer, and the decision variables are the capacity of cascade pumped storage, the configuration capacity and configuration power of electrochemical energy storage; the inner layer of the model is the operation layer, and the decision variables are the output of each unit and the power flow of the transmission network. The matrix form of the model is:

[0071]

[0072] stAy≤b;

[0073] Gx k +Hx k ≥gk=1,2,...,K;

[0074] Ux k +Vx k =wk=1,2,...,K;

[0075] Ey+Fx k ≥dk=1,2,...,K;

[0076] Where: y is the first-stage decision variable; x s is the second stage optimization variable; X is the second stage optimization variable set; x sis an uncertain parameter, which represents the typical output scenario of wind power and photovoltaic power generation after historical scenario clustering; A, b, G, H, g, U, V, w, E, F, and d are the corresponding coefficient matrices.

[0077] Further, the adopting a column and constraint generation algorithm to solve the two-layer model includes:

[0078] The two-level model of distributed robust planning and operation is decomposed to obtain the main problem and sub-problems:

[0079] Main question:

[0080]

[0081]

[0082] Sub-questions:

[0083]

[0084] Where: n is the cutting plane; J is the number of iterations; y * The outer decision variables obtained for the main problem;

[0085] The main problem obtains the optimal robust decision based on the worst scenario probability distribution set transmitted by the sub-problem, providing a lower bound for the model. The sub-problem obtains y in the main problem. * On this basis, find the worst probability distribution in the probability distribution uncertainty set and pass it to the main problem to provide an upper bound value for the model;

[0086] Yalmip is used to call the Gurobi solver to obtain the combined energy storage planning and operation method.

[0087] The present invention also includes a grid-side combined energy storage robust planning and operation device, using the above method, including:

[0088] A capacity configuration model building unit is used to consider various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, and to build a capacity configuration model for the combined energy storage;

[0089] The operation optimization model building unit is used to establish the system operation optimization model with the goal of minimizing the system operation cost, taking into account the operation constraints of cascade pumped storage and electrochemical energy storage, and the line flow constraints;

[0090] A two-layer model building unit is used to build a combined energy storage distribution robust planning-operation two-layer model including building wind power and photovoltaic output probability distribution uncertainty and linear model conversion;

[0091] A solving unit is used to solve the double-layer model by adopting a column and constraint generation algorithm to obtain a combined energy storage planning and operation method.

[0092] The present invention also includes a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method described above is implemented.

[0093] The present invention also includes a storage medium on which a computer program is stored. When the computer program is executed by a processor, the method described above is implemented.

[0094] The beneficial effects of the present invention are as follows: the present invention considers various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, and establishes a capacity configuration model for combined energy storage; considers operating constraints, line flow constraints and other conditions of cascade pumped storage and electrochemical energy storage, and establishes a system operation optimization model with the goal of minimizing system operating costs; establishes a combined energy storage distribution robust planning-operation two-layer model including constructing an uncertain set of probability distributions of wind power and photovoltaic output and linearizing model transformation; solves the model using a column and constraint generation algorithm to obtain a combined energy storage planning and operation method; the present invention can determine the capacity of combined energy storage under the coordination of source, grid, load and storage, and can optimize system operation, including: reducing the amount of abandoned electricity, alleviating the peak-shaving pressure of thermal power, and reducing the frequent start and stop of thermal power and hydropower units, and the result is robust. BRIEF DESCRIPTION OF THE DRAWINGS

[0095] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0096] Figure 1 is a flow chart of the method in Example 1;

[0097] Figure 2 It is a structural schematic diagram of the device in Example 1;

[0098] Figure 3 This is a schematic diagram of the structure of the cascade water-fire-wind-solar-storage system in Example 2;

[0099] Figure 4 This is a diagram showing the operation results of the combined power generation system on a typical day 1 in Example 2;

[0100] Figure 5 This is a diagram showing the operation results of the combined power generation system on a typical day 2 in Example 2;

[0101] Figure 6This is a diagram showing the operation results of the combined power generation system on a typical day 3 in Example 2;

[0102] Figure 7 This is a diagram showing the operation results of the 4-combined power generation system on a typical day in Example 2;

[0103] Figure 8 A schematic diagram of the structure of a computer device. DETAILED DESCRIPTION

[0104] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments.

[0105] Embodiment 1:

[0106] like Figure 1 As shown: A method for robust planning and operation of combined energy storage on the grid side includes the following steps:

[0107] Considering the various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, a capacity configuration model for combined energy storage is established;

[0108] With the goal of minimizing system operation cost, the operation constraints of cascade pumped storage and electrochemical energy storage, as well as line flow constraints, are considered to establish a system operation optimization model.

[0109] Establish a combined energy storage distributed robust planning-operation two-layer model that includes constructing the uncertainty of wind power and photovoltaic output probability distribution and linearizing model transformation;

[0110] The column and constraint generation algorithm is used to solve the two-layer model and obtain the combined energy storage planning and operation method.

[0111] By considering various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, a capacity configuration model of combined energy storage is established; considering the operating constraints, line flow constraints and other conditions of cascade pumped storage and electrochemical energy storage, and taking the minimum system operating cost as the goal, an operation optimization model of the system is established; a combined energy storage distribution robust planning-operation two-layer model is established, which includes constructing an uncertain set of probability distribution of wind power and photovoltaic output and linearizing model transformation; the model is solved using a column and constraint generation algorithm to obtain a combined energy storage planning and operation method; the present invention can determine the capacity of combined energy storage under the coordination of source, grid, load and storage, and can optimize system operation, including: reducing the amount of abandoned electricity, alleviating the peak-shaving pressure of thermal power, and reducing the frequent start and stop of thermal power and hydropower units, and the result is robust.

[0112] In this embodiment, the capacity configuration model of the combined energy storage includes:

[0113] minC=C in +Com -C re ;

[0114] in:

[0115]

[0116] C om =k om,ps I ps +k om,e E bat +k om,p P bat ;

[0117]

[0118] Where: C in , C om and C re are the annual equivalent investment cost, annual equivalent operation and maintenance cost, and annual equivalent residual value cost respectively; γ is the discount rate; Y ps and Y bat are the planned operating years of cascade pumped storage and electrochemical energy storage, respectively; k ps , k e and k p are the cost coefficients of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; I ps 、E bat and P bat are the capacity of cascade pumped storage, the capacity and power of electrochemical energy storage respectively; k om,ps , k om,e and k om,p are the annual operation and maintenance costs of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; C in,ps and C in,bat are the investment costs of cascade pumped storage and electrochemical energy storage respectively; δ ps and δ bat They are the residual value rates of cascade pumped storage and electrochemical energy storage, respectively.

[0119] The capacity configuration model includes the capacity and power constraints of cascade pumped storage and electrochemical energy storage:

[0120]

[0121] Where: and are the upper and lower limits of cascade pumped storage capacity respectively; and are the upper and lower limits of electrochemical energy storage capacity, respectively; and are the upper and lower limits of electrochemical energy storage power respectively.

[0122] Run the optimization model, including:

[0123] minC op =C th +C start +C cur ;

[0124] in:

[0125]

[0126]

[0127]

[0128] Where: C th , C start and C cur are the cost of thermal power units, the startup cost and the power abandonment cost of hydropower and pumped storage units respectively; N T 、N G and N H are the operating time, the number of thermal power units and the number of cascade hydropower stations respectively; a, b and c are the coefficients of the quadratic term, linear term and constant term of the thermal power coal consumption cost respectively; Provide power for thermal power units; and are the start and stop operation variables of thermal power units respectively; c G,st and c G,dn are the start-up and shutdown operation costs of thermal power units respectively; T is the carbon emission cost; Q and Q1 are the carbon emission and carbon emission quota respectively; and are the startup operation variables of the hydropower unit, pumped storage unit generating conditions and pumping conditions of the pumped storage unit respectively; c H,su 、c PSG,su and c PSP,su They are the startup costs of the hydropower unit, pumped storage unit generating conditions and pumped storage unit pumping conditions respectively; and They are the predicted output and actual output of photovoltaic power generation respectively; and are the predicted output and actual output of wind power respectively; c PV and c W are the penalty costs for curtailment of photovoltaic and wind power respectively.

[0129] Among them, the operation constraints of cascade pumped storage and electrochemical energy storage and line flow constraints are considered, including:

[0130] Power balance constraints:

[0131]

[0132] P t PS =P t PSG -P t PSP ;

[0133] P t bat =P t dis -P t cha ;

[0134] Where: P t PS Output of pumped storage unit; P t PSG and P t PSP are the power generation output and pumping output of the pumped storage unit respectively; P t bat Power for electrochemical energy storage; P t dis and P t cha are the discharge power and charging power of electrochemical energy storage respectively; is the load power;

[0135] Electrochemical energy storage operation constraints:

[0136]

[0137]

[0138]

[0139] Where: and are 0-1 variables for the electrochemical energy storage in the discharge state and the charge state respectively; E bat,min and E bat,max are the minimum and maximum capacity of electrochemical energy storage, respectively; It is the amount of electricity when electrochemical energy storage is working; h cha and h dis are charging power and discharging power respectively;

[0140] Hydropower unit output, start and stop, flow and storage capacity constraints:

[0141]

[0142]

[0143]

[0144]

[0145] Where: ρ is the density of water; g is the acceleration due to gravity; is the power generation efficiency of the hydropower unit; H i is the average water head; is the state variable of the hydropower unit; and are the minimum and maximum outputs of the hydropower unit respectively; I i,t , and They are the generating flow of hydropower units, abandoned water flow, discharged flow, interval flow, pumped storage working flow, pumped water flow of pumped storage units and generating flow; and They are the start and stop operation variables of the hydropower unit respectively; and They are the minimum flow rate and maximum flow rate of the hydropower unit for power generation respectively; is the minimum discarded water flow; and are the marker variables of the upper reservoir and the lower reservoir of the pumped storage unit respectively; V i,t is the storage capacity; V i min and V i max are the minimum storage capacity and the maximum storage capacity respectively;

[0146] Output and start / stop constraints of pumped storage units:

[0147]

[0148]

[0149] Where: h is the average water head of the pumped storage unit; h PSG and h PSP are the power generation efficiency and pumping efficiency of the pumped storage unit respectively; and are the power generation and pumping state variables of the pumped storage unit respectively; α PSG and α PSP are the minimum coefficients of the generating power and pumping power of the pumped storage unit relative to the configured capacity, respectively; and are the stop operation variables of the pumped storage unit for power generation and water pumping, respectively;

[0150] Output, start / stop and ramp-up constraints of thermal power units:

[0151]

[0152]

[0153]

[0154] Where: and are the minimum and maximum output of thermal power units respectively; is the state variable of the thermal power unit; and It is the up and down climbing power of the thermal power unit; and are the minimum time that the thermal power unit is in the on and off state respectively;

[0155] Wind power and photovoltaic operation constraints:

[0156]

[0157] Line flow constraints:

[0158]

[0159] Where: G n-g , G n-i , G n-ps , G n-w , G n-p and G n-k are the power transfer factors of thermal power unit g, cascade hydropower station i, pumped storage unit, wind turbine unit w, photovoltaic generator unit p and load node k to line n respectively; Maximum line power limit.

[0160] In this embodiment, a combined energy storage distribution robust planning-operation two-layer model is established, which includes constructing the uncertainty of wind power and photovoltaic output probability distribution and linearizing model conversion, including:

[0161] Using clustering methods to obtain the empirical probability of typical scenarios from historical data A comprehensive norm including 1-norm and ∞-norm is used to construct the uncertainty set of wind power and photovoltaic output probability distribution;

[0162]

[0163] Where: Ψ and Ψ ∞ They are the uncertainty sets of probability distribution established based on 1-norm and ∞-norm respectively; and They are the allowable deviation ranges of probability distribution under 1-norm and ∞-norm respectively; and As shown below:

[0164]

[0165] Where: M is the total amount of historical data; α1 and α ∞ are the uncertainty confidence levels of the 1-norm and ∞-norm respectively;

[0166] The uncertainty set of wind power and photovoltaic output probability distribution is constructed as follows:

[0167]

[0168] The linear model transformation is performed on the uncertainty set of probability distribution:

[0169]

[0170] A grid-side combined energy storage distributed robust planning-operation two-layer model considering the coordination of source, grid, load and storage is established. The outer layer of the model is the planning layer, and the decision variables are the capacity of cascade pumped storage, the configuration capacity and configuration power of electrochemical energy storage; the inner layer of the model is the operation layer, and the decision variables are the output of each unit and the power flow of the transmission network. The matrix form of the model is:

[0171]

[0172] stAy≤b;

[0173] Gx k +Hx k ≥gk=1,2,...,K;

[0174] Ux k +Vx k =wk=1,2,...,K;

[0175] Ey+Fx k ≥dk=1,2,...,K;

[0176] Where: y is the first-stage decision variable; x s is the second stage optimization variable; X is the second stage optimization variable set; x s is an uncertain parameter, which represents the typical output scenario of wind power and photovoltaic power generation after historical scenario clustering; A, b, G, H, g, U, V, w, E, F, and d are the corresponding coefficient matrices.

[0177] The two-layer model is solved using a column and constraint generation algorithm, including:

[0178] The Robust Planning-Running two-level model is decomposed into the main problem and sub-problems:

[0179] Main question:

[0180]

[0181]

[0182] Sub-questions:

[0183]

[0184] Where: n is the cutting plane; J is the number of iterations; y * The outer decision variables obtained for the main problem;

[0185] The main problem obtains the optimal robust decision based on the worst scenario probability distribution set transmitted by the sub-problem, providing a lower bound for the model. The sub-problem obtains y in the main problem. * On this basis, find the worst probability distribution in the probability distribution uncertainty set and pass it to the main problem to provide an upper bound value for the model;

[0186] Yalmip is used to call the Gurobi solver to obtain the combined energy storage planning and operation method.

[0187] like Figure 2 As shown, this embodiment also includes a grid-side combined energy storage robust planning and operation device, using the above method, including:

[0188] A capacity configuration model building unit is used to consider various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, and to build a capacity configuration model for the combined energy storage;

[0189] The operation optimization model building unit is used to establish the system operation optimization model with the goal of minimizing the system operation cost, taking into account the operation constraints of cascade pumped storage and electrochemical energy storage, and the line flow constraints;

[0190] A two-layer model building unit is used to build a combined energy storage distribution robust planning-operation two-layer model including building wind power and photovoltaic output probability distribution uncertainty and linear model conversion;

[0191] The solving unit is used to solve the double-layer model by adopting the column and constraint generation algorithm to obtain the combined energy storage planning and operation method.

[0192] Embodiment 2:

[0193] In order to demonstrate the effectiveness of the day-ahead peak-shaving robust optimization scheduling method for photovoltaic complementary power generation provided by the present invention, it is described below in conjunction with specific application scenarios.

[0194] like Figure 3As shown, the embodiment of the present invention includes three cascade hydropower stations, a wind farm, a photovoltaic power station, a thermal power plant, and a cascade pumped-storage power station and an electrochemical energy storage power station with a capacity to be determined. The total installed capacity of the cascade hydropower stations is 1900MW, the installed capacity of the wind farm is 2400MW, the installed capacity of the photovoltaic power station is 2400MW, and the installed capacity of the thermal power station is 1720MW.

[0195] First, the various costs of combined energy storage including cascade pumped storage and electrochemical energy storage are considered, and a capacity configuration model of combined energy storage is established, including:

[0196] The goal of grid-side combined energy storage planning is to minimize the annual comprehensive cost;

[0197] minC=C in +C om -C re

[0198] in:

[0199]

[0200] C om =k om,ps I ps +k om,e E bat +k om,p P bat

[0201]

[0202] Where: C in , C om and C re are the annual equivalent investment cost, annual equivalent operation and maintenance cost, and annual equivalent residual value cost respectively; γ is the discount rate; Y ps and Y bat are the planned operating years of cascade pumped storage and electrochemical energy storage, respectively; k ps , k e and k p are the cost coefficients of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; I ps 、E bat and P bat are the capacity of cascade pumped storage, the capacity and power of electrochemical energy storage respectively; k om,ps , k om,e and k om,p are the annual operation and maintenance costs of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; C in,ps and C in,bat are the investment costs of cascade pumped storage and electrochemical energy storage respectively; δ ps and δ batare the residual value rates of cascade pumped storage and electrochemical energy storage, respectively;

[0203] The constraints for grid-side combined energy storage planning include the capacity and power constraints of cascade pumped storage and electrochemical energy storage:

[0204]

[0205] Where: and are the upper and lower limits of cascade pumped storage capacity respectively; and are the upper and lower limits of electrochemical energy storage capacity, respectively; and are the upper and lower limits of electrochemical energy storage power respectively.

[0206] Then, considering the operation constraints of cascade pumped storage and electrochemical energy storage, line flow constraints and other conditions, with the goal of minimizing the system operation cost, the system operation optimization model is established, including:

[0207] The optimized operation of grid-side combined energy storage needs to take into account the operating costs of each unit, such as the coal consumption cost, start-up and shutdown cost, and carbon emission cost of thermal power units, as well as the startup cost of hydropower units and cascade pumped storage and the penalty cost of wind and solar power abandonment. On the premise of determining the cascade pumped storage capacity and electrochemical energy storage capacity and power, the output of each unit should be optimized;

[0208] The objective function of the grid-side combined energy storage optimization operation is as follows:

[0209] minC op =C th +C start +C cur

[0210] in:

[0211]

[0212]

[0213]

[0214] Where: C th , C start and C cur are the cost of thermal power units, the startup cost and the power abandonment cost of hydropower and pumped storage units respectively; N T 、N G and N H are the operating time, the number of thermal power units and the number of cascade hydropower stations respectively; a, b and c are the coefficients of the quadratic term, linear term and constant term of the thermal power coal consumption cost respectively; Provide power for thermal power units; and are the start and stop operation variables of thermal power units respectively; c G,st and c G,dn are the start-up and shutdown operation costs of thermal power units respectively; T is the carbon emission cost; Q and Q1 are the carbon emission and carbon emission quota respectively; and are the startup operation variables of the hydropower unit, pumped storage unit generating conditions and pumping conditions of the pumped storage unit respectively; c H,su 、c PSG,su and c PSP,su They are the startup costs of the hydropower unit, pumped storage unit generating conditions and pumped storage unit pumping conditions respectively; and They are the predicted output and actual output of photovoltaic power generation respectively; and are the predicted output and actual output of wind power respectively; c PV and c W are the penalty costs for curtailment of photovoltaic and wind power respectively;

[0215] The constraints for the optimal operation of combined energy storage on the grid side are as follows:

[0216] 1) Power balance constraints:

[0217]

[0218]

[0219] Where: P t PS Output of pumped storage unit; P t PSG and P t PSP are the power generation output and pumping output of the pumped storage unit respectively; P t bat Power for electrochemical energy storage; P t dis and P t cha are the discharge power and charging power of electrochemical energy storage respectively; is the load power;

[0220] 2) Electrochemical energy storage operation constraints:

[0221]

[0222]

[0223]

[0224] Where: and are 0-1 variables for the electrochemical energy storage in the discharge state and the charge state respectively; E bat,min and E bat,max are the minimum and maximum electric capacity of electrochemical energy storage respectively; E t bat It is the amount of electricity when electrochemical energy storage is working; h cha and h dis are charging power and discharging power respectively;

[0225] 3) Hydropower unit output, start and stop, flow and reservoir capacity constraints:

[0226]

[0227]

[0228]

[0229]

[0230] Where: ρ is the density of water; g is the acceleration due to gravity; h i H is the power generation efficiency of the hydropower unit; H i is the average water head; is the state variable of the hydropower unit; and P i H,max are the minimum and maximum outputs of the hydropower unit respectively; I i,t , and They are the generating flow of hydropower units, abandoned water flow, discharged flow, interval flow, pumped storage working flow, pumped water flow of pumped storage units and generating flow; and They are the start and stop operation variables of the hydropower unit respectively; and They are the minimum flow rate and maximum flow rate of the hydropower unit for power generation respectively; is the minimum discarded water flow; and are the marker variables of the upper reservoir and the lower reservoir of the pumped storage unit respectively; V i,t is the storage capacity; V i min and V i max are the minimum storage capacity and the maximum storage capacity respectively;

[0231] Cascade hydropower parameters are shown in Table 1:

[0232] Table 1 Cascade hydropower parameters

[0233]

[0234]

[0235] 4) Output and start / stop constraints of pumped storage units:

[0236]

[0237]

[0238] Where: h is the average water head of the pumped storage unit; h PSG and h PSP are the power generation efficiency and pumping efficiency of the pumped storage unit respectively; and are the power generation and pumping state variables of the pumped storage unit respectively; α PSG and α PSP are the minimum coefficients of the generating power and pumping power of the pumped storage unit relative to the configured capacity, respectively; and are the stop operation variables of the pumped storage unit for power generation and water pumping, respectively;

[0239] 5) Thermal power unit output, start / stop, and climbing constraints:

[0240]

[0241]

[0242]

[0243] Where: and are the minimum and maximum output of thermal power units respectively; is the state variable of the thermal power unit; and It is the up and down climbing power of the thermal power unit; and are the minimum time that the thermal power unit is in the on and off state respectively;

[0244] 6) Wind power and photovoltaic operation constraints:

[0245]

[0246] The wind power output forecast values ​​are shown in Table 2:

[0247] Table 2 Wind power output forecast values

[0248]

[0249] The predicted value of photovoltaic power generation output is shown in Table 3:

[0250] Table 3 Prediction of photovoltaic power generation output

[0251]

[0252]

[0253] 7) Line flow constraints:

[0254]

[0255] Where: G n-g , G n-i , G n-ps , G n-w , G n-p and G n-k are the power transfer factors of thermal power unit g, cascade hydropower station i, pumped storage unit, wind turbine unit w, photovoltaic generator unit p and load node k to line n respectively; Maximum line power limit.

[0256] Next, a combined energy storage distributed robust planning-operation two-layer model is established, which includes constructing the uncertainty set of wind power and photovoltaic output probability distribution and linearization model transformation, including:

[0257] Using clustering methods to obtain the empirical probability of typical scenarios from historical data Since the actual probability p k There is a deviation between (k=1,2,...,K) and the empirical probability. For this reason, a comprehensive norm including 1-norm and ∞-norm is used to construct the uncertainty set of wind power and photovoltaic output probability distribution;

[0258]

[0259] Where: Ψ and Ψ ∞ They are the uncertainty sets of probability distribution established based on 1-norm and ∞-norm respectively; and They are the allowable deviation ranges of probability distribution under 1-norm and ∞-norm respectively; and As shown below:

[0260]

[0261] Where: M is the total amount of historical data; α1 and α ∞ are the uncertainty confidence of 1-norm and ∞-norm respectively; thus, the uncertainty set of wind power and photovoltaic output probability distribution can be constructed as follows:

[0262]

[0263] The above probability distribution uncertainty set has absolute value constraints, and it needs to be transformed into a linear model to facilitate the solution of the model:

[0264]

[0265] At this point, considering the capacity configuration problem of combined energy storage and the optimal operation problem of cascaded water storage, thermal, wind and solar storage systems, a grid-side combined energy storage distributed robust planning-operation two-layer model considering the coordination of source, grid, load and storage can be established. The outer layer of the model is the planning layer, and the decision variables are the capacity of the cascade pumped storage, the configured capacity and configured power of the electrochemical energy storage; the inner layer of the model is the operation layer, and the decision variables are the output of each unit, the power flow of the transmission network, etc. The matrix form of the model is:

[0266]

[0267] stAy≤b

[0268] Gx k +Hx k ≥gk=1,2,...,K

[0269] Ux k +Vx k =wk=1,2,...,K

[0270] Ey+Fx k ≥dk=1,2,...,K

[0271] Where: y is the first-stage decision variable; x s is the second stage optimization variable; X is the second stage optimization variable set; x s is an uncertain parameter, which represents the typical output scenario of wind power and photovoltaic power generation after historical scenario clustering; A, b, G, H, g, U, V, w, E, F, and d are the corresponding coefficient matrices.

[0272] Finally, the model is solved using the column and constraint generation algorithm to obtain the combined energy storage planning and operation method, including:

[0273] After decomposition, the main problem and sub-problems can be obtained by running the two-level model:

[0274] Main question:

[0275]

[0276]

[0277] Sub-questions:

[0278]

[0279] Where: n is the cutting plane; J is the number of iterations; y * The outer decision variables obtained for the main problem;

[0280] The main problem obtains the optimal robust decision based on the worst scenario probability distribution set transmitted by the sub-problem, providing a lower bound for the model. * On this basis, find the worst probability distribution in the probability distribution uncertainty set and pass it to the main problem to provide an upper bound value for the model;

[0281] Yalmip is used to call the Gurobi solver to obtain the combined energy storage planning and operation method.

[0282] After solving, the combined energy storage planning and operation method is obtained, and four typical days are selected Figure 4 , Figure 5 , Figure 6 , Figure 7 and Table 4.

[0283] Table 4 Optimized operation plan

[0284]

[0285]

[0286] See also Figure 8 A computer device 400 provided in an embodiment of the present application includes: a processor 410 and a memory 420, wherein the memory 420 stores a computer program executable by the processor 410, and when the computer program is executed by the processor 410, the above method is executed.

[0287] The embodiment of the present application further provides a storage medium 430 on which a computer program is stored. When the computer program is run by the processor 410, the above method is executed.

[0288] Among them, the storage medium 430 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable red-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, disk or optical disk.

[0289] In the description of the present invention, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. "Multiple" means two or more, unless otherwise clearly and specifically defined.

[0290] In the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0291] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0292] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code that includes one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention belong.

[0293] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic devices), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0294] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0295] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0296] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present invention have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present invention. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for robust planning and operation of combined energy storage on the grid side, characterized in that: The steps include: Considering the various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, a capacity configuration model for combined energy storage is established; With the goal of minimizing system operation cost, the operation constraints of cascade pumped storage and electrochemical energy storage, as well as line flow constraints, are considered to establish a system operation optimization model. Establish a combined energy storage distributed robust planning-operation two-layer model that includes constructing the uncertainty of wind power and photovoltaic output probability distribution and linearizing model transformation; The double-layer model is solved by using a column and constraint generation algorithm to obtain a combined energy storage planning and operation method; The capacity configuration model of the combined energy storage includes: minC=C in +C om -C re ; in: C om =k om,ps I ps +k om,e E bat +k om,p P bat ; Where: C in , C om and C re are the annual equivalent investment cost, annual equivalent operation and maintenance cost, and annual equivalent residual value cost respectively; γ is the discount rate; Y ps and Y bat are the planned operating years of cascade pumped storage and electrochemical energy storage, respectively; k ps , k e and k p are the cost coefficients of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; I ps 、E bat and P bat are the capacity of cascade pumped storage, the capacity and power of electrochemical energy storage respectively; k om,ps , k om,e and k om,p are the annual operation and maintenance costs of cascade pumped storage capacity, electrochemical energy storage capacity and power respectively; C in,ps and C in,bat are the investment costs of cascade pumped storage and electrochemical energy storage respectively; δ ps and δ bat are the residual value rates of cascade pumped storage and electrochemical energy storage, respectively; The capacity configuration model includes the capacity and power constraints of cascade pumped storage and electrochemical energy storage: Where: and are the upper and lower limits of cascade pumped storage capacity respectively; and are the upper and lower limits of electrochemical energy storage capacity, respectively; and are the upper and lower limits of electrochemical energy storage power respectively; The running optimization model comprises: minC op =C th +C start +C cur ; in: Where: C th , C start and C cur are the cost of thermal power units, the startup cost and the power abandonment cost of hydropower and pumped storage units respectively; N T 、N G and N H are the operating time, the number of thermal power units and the number of cascade hydropower stations respectively; a, b and c are the coefficients of the quadratic term, linear term and constant term of the thermal power coal consumption cost respectively; Provide power for thermal power units; and are the start and stop operation variables of thermal power units respectively; c G,su and c G,sd are the start-up and shutdown operation costs of thermal power units respectively; T is the carbon emission cost; Q and Q1 are the carbon emission and carbon emission quota respectively; and are the startup operation variables of the hydropower unit, pumped storage unit generating conditions and pumping conditions of the pumped storage unit respectively; c H,su 、c PSG,su and c PSP,su They are the startup costs of the hydropower unit, pumped storage unit generating conditions and pumped storage unit pumping conditions respectively; and They are the predicted output and actual output of photovoltaic power generation respectively; and are the predicted output and actual output of wind power respectively; c PV and c W are the penalty costs for curtailment of photovoltaic and wind power respectively; The operation constraints and line flow constraints of cascade pumped storage and electrochemical energy storage include: Power balance constraints: P t PS =P t PSG -P t PSP ; P t bat =P t dis -P t cha ; Where: P t PS Output of pumped storage unit; P t PSG and P t PSP are the power generation output and pumping output of the pumped storage unit respectively; P t bat Power for electrochemical energy storage; P t dis and P t cha are the discharge power and charging power of electrochemical energy storage respectively; is the load power; Electrochemical energy storage operation constraints: Where: and are 0-1 variables for the electrochemical energy storage in the discharge state and the charge state respectively; E bat,min and E bat,max are the minimum and maximum capacity of electrochemical energy storage, respectively; The amount of electricity during electrochemical energy storage operation; η cha and η dis are charging efficiency and discharging efficiency respectively; Hydropower unit output, start and stop, flow and storage capacity constraints: Where: ρ is the density of water; g is the acceleration due to gravity; η i H is the power generation efficiency of the hydropower unit; H i is the average water head; is the state variable of the hydropower unit; and P i H,max are the minimum and maximum outputs of the hydropower unit respectively; I i,t , and They are the generating flow of hydropower units, abandoned water flow, discharged flow, interval flow, pumped storage working flow, pumped water flow of pumped storage units and generating flow; and They are the start and stop operation variables of the hydropower unit respectively; and They are the minimum flow rate and maximum flow rate of the hydropower unit for power generation respectively; is the minimum discarded water flow; and are the marker variables of the upper reservoir and the lower reservoir of the pumped storage unit respectively; V i,t is the storage capacity; V i min and V i max are the minimum storage capacity and the maximum storage capacity respectively; Output and start / stop constraints of pumped storage units: Where: h is the average water head of the pumped storage unit; η PSG and η PSP are the power generation efficiency and pumping efficiency of the pumped storage unit respectively; and are the power generation and pumping state variables of the pumped storage unit respectively; α PSG and α PSP are the minimum coefficients of the generating power and pumping power of the pumped storage unit relative to the configured capacity, respectively; and They are the start and stop operation variables of the pumped storage unit generating condition respectively; and are the start and stop operation variables of the pumping condition of the pumped storage unit respectively; Output, start / stop and ramp-up constraints of thermal power units: Where: and are the minimum and maximum output of thermal power units respectively; is the state variable of the thermal power unit; and It is the up and down climbing power of the thermal power unit; and are the minimum time that the thermal power unit is in the on and off state respectively; Wind power and photovoltaic operation constraints: Line flow constraints: Where: G n-g , G n-i , G n-ps , G n-w , G n-p and G n-k are the power transfer factors of thermal power unit g, cascade hydropower station i, pumped storage unit, wind turbine unit w, photovoltaic generator unit p and load node k to line n respectively; The maximum line limit power; The establishment includes constructing a combined energy storage distribution robust planning-operation two-layer model including building wind power and photovoltaic output probability distribution uncertainty and linearization model conversion, including: Using clustering methods to obtain the empirical probability of typical scenarios from historical data A comprehensive norm including 1-norm and ∞-norm is used to construct the uncertainty set of wind power and photovoltaic output probability distribution; Where: Ψ1 and Ψ ∞ They are the uncertainty sets of probability distribution established based on 1-norm and ∞-norm respectively; and They are the allowable deviation ranges of probability distribution under 1-norm and ∞-norm respectively; and As shown below: Where: M is the total amount of historical data; α1 and α ∞ are the uncertainty confidence levels of the 1-norm and ∞-norm respectively; The uncertainty set of wind power and photovoltaic output probability distribution is constructed as follows: The probability distribution uncertainty set is transformed into a linear model: A grid-side combined energy storage distributed robust planning-operation two-layer model considering the coordination of source, grid, load and storage is established. The outer layer of the model is the planning layer, and the decision variables are the capacity of cascade pumped storage, the configuration capacity and configuration power of electrochemical energy storage; the inner layer of the model is the operation layer, and the decision variables are the output of each unit and the power flow of the transmission network. The matrix form of the model is: stAy≤b; Gx k +Hξ k ≥gk=1,2,...,K; Ux k +Vξ k =w k=1,2,...,K; Ey+Fx k ≥d k=1,2,...,K; Where: y is the first-stage decision variable; x s is the second stage optimization variable; X is the second stage optimization variable set; ξ s is an uncertain parameter, which represents the typical output scenario of wind power and photovoltaic power generation after historical scenario clustering; A, b, G, H, g, U, V, w, E, F, and d are the corresponding coefficient matrices.

2. The grid-side combined energy storage robust planning and operation method according to claim 1 is characterized in that: The adopting a column and constraint generation algorithm to solve the two-layer model includes: The two-level model of distributed robust planning and operation is decomposed to obtain the main problem and sub-problems: Main question: Sub-questions: Where: v is the cutting plane; J is the number of iterations; y * The outer decision variables obtained for the main problem; The main problem obtains the optimal robust decision based on the worst scenario probability distribution set transmitted by the sub-problem, providing a lower bound for the model. The sub-problem obtains y in the main problem. * On this basis, find the worst probability distribution in the probability distribution uncertainty set and pass it to the main problem to provide an upper bound value for the model; Yalmip is used to call the Gurobi solver to obtain the combined energy storage planning and operation method.

3. A grid-side combined energy storage robust planning and operation device, characterized in that: Use of the method as claimed in claim 1 or 2, comprising: A capacity configuration model building unit is used to consider various costs of combined energy storage including cascade pumped storage and electrochemical energy storage, and to build a capacity configuration model for the combined energy storage; The operation optimization model building unit is used to establish the system operation optimization model with the goal of minimizing the system operation cost, taking into account the operation constraints of cascade pumped storage and electrochemical energy storage, and the line flow constraints; A two-layer model building unit is used to build a combined energy storage distribution robust planning-operation two-layer model including building wind power and photovoltaic output probability distribution uncertainty and linear model conversion; A solving unit is used to solve the double-layer model by adopting a column and constraint generation algorithm to obtain a combined energy storage planning and operation method.

4. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method as claimed in claim 1 or 2 is implemented.

5. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method as claimed in claim 1 or 2 is implemented.

Images