CN118572703B

CN118572703B - Step water storage wind-solar-fire planning operation method, system, equipment and medium

Info

Publication number: CN118572703B
Application number: CN202411061411.1A
Authority: CN
Inventors: 但扬清; 田雪沁; 陈晴悦; 王龙; 何英静; 孙飞飞; 王智冬; 沈志恒; 吕学; 黄锦华; 黄怡; 许恩超; 王晨轩; 丁一凡; 张夏辉
Original assignee: State Grid Zhejiang Electric Power Co Ltd; Economic and Technological Research Institute of State Grid Zhejiang Electric Power Co Ltd; State Grid Economic and Technological Research Institute
Current assignee: State Grid Zhejiang Electric Power Co Ltd; Economic and Technological Research Institute of State Grid Zhejiang Electric Power Co Ltd; State Grid Economic and Technological Research Institute
Priority date: 2024-08-05
Filing date: 2024-08-05
Publication date: 2024-10-18
Anticipated expiration: 2044-08-05
Also published as: CN118572703A

Abstract

The invention relates to the technical field of multi-energy complementary power generation, in particular to a step water storage wind-solar-fire planning operation method, a system, equipment and a medium, which comprise the steps of taking the minimum annual comprehensive cost and the electricity rejection rate of a system as optimization targets, establishing a step multi-energy collaborative multi-target optimization function, linearizing the step multi-energy collaborative multi-target optimization function and a nonlinear part of step water storage wind-solar-fire constraint conditions, and constructing a mixed integer linear programming model of step water storage wind-solar-fire capacity optimization configuration; converting the cascade multipotency cooperative multi-objective optimization function into a single-objective optimization problem, solving a mixed integer linear programming model, and obtaining a capacity configuration operation scheduling strategy of a cascade multipotency complementary system; and controlling the running state of the cascade water storage wind-light-fire system according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system. Under the synergistic effect of pumping storage and thermal power, the invention realizes the optimized operation of the cascade water storage wind-solar-fire system, and ensures that the system can operate efficiently and stably.

Description

Step water storage wind-solar-fire planning operation method, system, equipment and medium

Technical Field

The invention relates to the technical field of multi-energy complementary power generation, in particular to a step water storage wind-solar-fire planning operation method, a system, equipment and a medium.

Background

Wind power and photovoltaic energy have space-time randomness and fluctuation in output, which leads to instability of electric energy output, and in order to overcome the challenge, the clean energy base usually adopts a multi-energy complementary strategy to fuse multiple energy forms, and aims to effectively stabilize the fluctuation of wind power and photovoltaic power generation through mutual coordination and complementation of the multiple energy forms, so that the sustainability and stability of energy supply are ensured.

The thermal power and the pumped storage are taken as the energy regulation sources, the renewable energy sources can be stabilized, the thermal power can reduce the output when the demands are low by virtue of the deep peak regulation capability, and the thermal power rapidly responds in the peak period to provide necessary support for the power grid; the pumped storage realizes the space-time allocation of energy through the flexible conversion of the pumped storage and the water discharge power generation, so that the adjustment elasticity of the power system is further enhanced, however, a plurality of researches at present often separate the two, and the synergistic effect caused by the combination of the two is ignored.

In addition, the development potential of the pure pumped storage power station gradually tends to be saturated in the face of the vigorous development of new energy development, but the pure pumped storage power station also has rich water and electricity resources in China, which provides space for innovative development of the pumped storage technology, existing water and electricity facilities are fused and modified, particularly existing infrastructure of a cascade hydropower station is utilized, the water and electricity resources can be converted into the efficient pumped storage power station without additionally constructing a reservoir, and therefore the flexibility and the reliability of an electric power system are greatly improved on the premise of not increasing environmental burden, and therefore, how to construct a cascade water, wind, light and fire system planning and operation method with efficient cooperation and optimal configuration between renewable energy sources such as thermal power, fused and modified water storage and wind power, photovoltaic and the like becomes a great topic to be solved in the current energy field.

Disclosure of Invention

The invention aims to provide a step water storage wind-solar-fire planning operation method, a system, equipment and a medium, which optimize the operation states of step water storage electricity, wind power and photovoltaic power generation under the synergistic effect of a pumped storage power station and a thermal power generation so as to reduce the electric quantity and improve the operation efficiency and the stability of the system.

In order to solve the technical problems, the invention provides a step water storage wind-solar-fire planning operation method, a system, equipment and a medium.

In a first aspect, the invention provides a step water storage wind-solar-fire planning operation method, which comprises the following steps:

Taking the minimum annual comprehensive cost and the electricity rejection rate of the system as optimization targets, establishing a cascade multi-energy cooperative multi-target optimization function, and establishing a cascade water storage wind-solar-fire constraint condition considering the upstream and downstream cascade hydropower constraint;

linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions, and constructing a mixed integer linear programming model of cascade water storage wind-light-fire capacity optimization configuration;

Converting the cascade multipotency cooperative multi-objective optimization function into a single-objective optimization problem by using a normal boundary crossing method, solving the mixed integer linear programming model, and obtaining a capacity configuration operation scheduling strategy of the cascade multipotency complementary system;

and controlling the running state of the cascade water storage wind-light-fire system according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system.

In a further embodiment, the step-wise multi-energy collaborative multi-objective optimization function is specifically:

Wherein,

Wherein F is an objective function; the comprehensive cost is the annual cost of the system; The power rejection rate is the electricity rejection rate; Investment cost for the system year; the annual running cost; the actual output of wind power at the time t of the s-th scene; The actual output of photovoltaic power generation at the time t of the s-th scene; predicting a force coefficient for the force of wind power at the time t of the s-th scene; Predicting a force coefficient for the force of photovoltaic power generation at the time t of the s-th scene; planning and constructing capacity for the wind farm; planning and constructing capacity for the photovoltaic power station; Is the scene number; t is the number of operating hours per scene; Is annual interest rate; the operation life of the pumping and accumulating unit is planned; the unit investment cost of the pumping and accumulating unit is set; planning and constructing capacity for the pumping and accumulating unit; planning operation years for the wind farm; The unit investment cost of the wind farm; planning and constructing capacity for the wind farm; planning the operation life of the photovoltaic power station; the unit investment cost of the photovoltaic power station; planning and constructing capacity for the photovoltaic power station; Probability of occurrence for the s-th scene; the number of the thermal power generating units; the output of the g-th thermal power generating unit is the s-th scene t moment; starting up an operation variable for a g-th thermal power generating unit at a s-th scene t moment; starting up cost for the g-th thermal power generating unit; the shutdown operation variable of the g-th thermal power generating unit is the s-th scene t moment; the shutdown cost of the g-th thermal power generating unit; the water level is the step hydropower level; the number of hydroelectric generating sets of the ith hydropower station; starting operation variables of a jth hydropower unit of an ith hydropower station at a time t of an s-th scene; the j-th station of the i-th hydropower station the starting cost of the hydroelectric generating set; the j-th water discarding cost of the i-th hydropower station; The method comprises the steps that the j-th waste water flow of the i-th hydropower station at the t moment of the s-th scene is obtained; Is a time interval; starting operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; The starting cost of the pumping working condition of the pumping and accumulating unit; the method comprises the steps of drawing a starting operation variable of a power generation working condition of a storage unit at the time t of an s-th scene; The starting cost of the power generation working condition of the pumping and accumulating unit is set; a. b and c are the coal consumption cost correlation coefficients of the thermal power unit; is a peak shaving cost function; p is the output variable of the thermal power unit; The coal consumption cost function of the thermal power generating unit is as follows; the cost of purchasing the thermal power generating unit is; Rotor fracturing cycle times at the moment t; The price of the fuel oil; The oil quantity is calculated; the price of the unit coal; The maximum output value of the thermal power unit is obtained; a first medium limit value for the output of the thermal power generating unit; a second medium limit value for the thermal power unit output; The lower limit value of the output of the thermal power generating unit is set; Is the resource consumption rate; is a time integral variable; Is a risk rate; is the reference time point.

In a further embodiment, the step water storage wind-solar-fire constraint condition comprises an upstream step water electricity constraint; the upstream and downstream cascade hydropower constraints comprise water balance constraints, hydropower output constraints, flow constraints, water head characteristic constraints, water level-reservoir capacity relation constraints, tail water level-lower drainage flow relation constraints, pumping storage capacity constraints, pumping storage start-stop constraints, pumping storage power upper and lower limit constraints and pumping storage output constraints;

The water balance constraint is specifically as follows:

In the method, in the process of the invention, The capacity of the ith level reservoir is the time t of the s-th scene; the interval flow of the ith level reservoir at the time t of the s-th scene; The outlet flow of the ith grade reservoir at the time t of the s-th scene; pumping and accumulating unit flow for the s-th scene t moment; Is a time interval; And The minimum storage capacity and the maximum storage capacity of the ith reservoir are respectively; And Respectively regulating the initial and final reservoir capacity coefficients of the reservoir in year; the capacity of the ith level reservoir is the initial moment of the s-th scene; The capacity of the ith level reservoir is the time T of the s-th scene;

the water power constraint is specifically as follows:

In the method, in the process of the invention, Is water electric power; Is the density of water; g is gravity acceleration; Generating efficiency of a j-th hydroelectric generating set of the i-th hydropower station; The water head height of the ith hydropower station at the time t of the s-th scene; Generating flow of a jth hydropower unit of the ith hydropower station at the time t of the jth scene;

The flow constraint is specifically:

In the method, in the process of the invention, The pumping flow of the pumping and accumulating unit at the s-th scene t moment is used; generating flow of the pumping and accumulating unit at the time t of the s-th scene; the method comprises the steps of (1) setting the drainage flow of an ith hydropower station as an ith hydropower station; generating flow for a jth hydroelectric generating set of an ith hydropower station at a time t of an ith scene; the method comprises the steps of (1) discarding water flow for an ith hydropower station of an ith hydropower station; The state variable of the jth hydropower unit of the ith hydropower station at the time t of the s scene; And The minimum power generation flow and the maximum power generation flow of the jth hydroelectric generating set of the ith hydropower station are respectively; Maximum water reject flow for the i-th hydropower station;

The head characteristic constraint is specifically:

In the method, in the process of the invention, The water level of the ith level reservoir at the time t of the s-th scene; The tail water level of the ith reservoir at the time t of the s-th scene; a minimum water head for the ith reservoir; Maximum water head for the ith reservoir;

the water level-reservoir capacity relation constraint is specifically:

In the method, in the process of the invention, The relation function of the level of the ith reservoir and the reservoir capacity is obtained;

The relation constraint of tail water level-downdraft flow is specifically as follows:

In the method, in the process of the invention, The relation function of the tail water level of the ith reservoir and the downward discharge flow is used;

The pumping and storage capacity constraint is specifically as follows:

In the method, in the process of the invention, The capacity of the pumping and accumulating unit; The capacity of the pumping and accumulating unit is the maximum; the capacity of the pumping and accumulating unit is the minimum;

The pumping, storage, starting and stopping constraint is specifically as follows:

In the method, in the process of the invention, The method comprises the steps of drawing and accumulating state variables under the power generation working condition of a unit at the time t of the s-th scene; Drawing and accumulating unit for the s-th scene t moment state variables under pumping conditions; The method comprises the steps of drawing a stopping operation variable of the power generation working condition of the storage unit at the time t of the s-th scene; Stopping operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; starting operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; the method comprises the steps of drawing a starting operation variable of a power generation working condition of a storage unit at the time t of an s-th scene;

The upper limit and the lower limit of the pumping and storage power are specifically defined as follows:

In the method, in the process of the invention, The lower limit power coefficient of the power generation working condition of the pumping and accumulating unit; The actual output of the power generation working condition of the pumping and accumulating unit is the s-th scene t moment; the lower limit power coefficient of the pumping working condition of the pumping and accumulating unit; the actual output of the pumping working condition of the pumping and accumulating unit is set at the time t of the s-th scene;

The extraction and storage force constraint is specifically as follows:

In the method, in the process of the invention, Generating efficiency is the generating working condition of the pumping and accumulating unit; h is the average water head of the pumping and accumulating unit; and the pumping efficiency of the pumping working condition of the pumping and accumulating unit.

In a further embodiment, the step water storage wind-solar-fire constraint conditions further comprise thermal power operation constraint, wind power and photovoltaic power generation constraint and system power balance constraint;

the thermal power operation constraint is specifically as follows:

In the method, in the process of the invention, The state constraint of the g-th thermal power generating unit is the s-th scene t moment; the minimum output of the g-th thermal power unit; maximum output of the g-th thermal power unit; the output of the g-th thermal power generating unit is the s-th scene t moment; the minimum start-up time of the g-th thermal power generating unit; the minimum shutdown time of the g-th thermal power unit;

the wind power and photovoltaic power generation constraint is specifically as follows:

In the method, in the process of the invention, The new energy installation capacity lower limit proportionality coefficient; the total assembly quantity of the multi-energy complementary system;

the system power balance constraint is specifically:

In the method, in the process of the invention, Is the load at time t of the s-th scene.

In a further embodiment, the step of linearizing the step multi-energy cooperative multi-objective optimization function and the nonlinear part of the step water storage wind-light-fire constraint condition to construct a mixed integer linear programming model of the step water storage wind-light-fire capacity optimization configuration includes:

linearizing a thermal power unit coal consumption cost function in a water level-reservoir capacity relation constraint, a tail water level-lower leakage flow relation constraint and a cascade multi-energy cooperative multi-objective optimization function by utilizing a piecewise linear method to obtain a corresponding water level-reservoir capacity relation linear constraint, tail water level-lower leakage flow relation linear constraint and thermal power unit coal consumption cost linear function;

Linearizing the hydropower output constraint and the upper and lower limit constraints of the pumping power by utilizing a convex hull relaxation method to obtain corresponding hydropower output linear constraint and upper and lower limit linear constraint of the pumping power;

constructing a linear constraint of a water level-reservoir capacity relationship, a linear constraint of a tail water level-lower leakage flow relationship, a linear function of coal consumption cost of the thermal power generating unit, a linear constraint of hydroelectric power output and an upper limit and a lower limit linear constraint of pumping and storing power into a linear set;

And constructing a mixed integer linear programming model of the cascade water storage wind-light-fire capacity optimization configuration according to the linearization set, the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint condition.

In a further embodiment, the step of converting the step of multi-energy collaborative multi-objective optimization function into a single objective optimization problem using a normal boundary crossing method comprises:

solving the minimum value of each single objective function in the cascade multi-energy collaborative multi-objective optimization function to obtain a single objective optimal solution;

Substituting the single-target optimal solution into each single-target function of the cascade multi-energy collaborative multi-target optimization function respectively, constructing a payment matrix, and connecting diagonal points of the payment matrix to obtain a Uotoban line;

carrying out normalization processing on a single objective function of the cascade multi-energy collaborative multi-objective optimization function to obtain a normalized objective function value of an arbitrary solution;

Normalizing the payment matrix to obtain a normalized payment matrix and normalized Utobang wires;

and randomly selecting a point on the normalized Utobang line by using a normal boundary crossing method as a target point, calculating the distance from the target point to each point on the Pareto front according to the normalized objective function value, and searching the Pareto front point with the maximized distance by changing the target point to obtain a Pareto solution set of the step multi-energy cooperative multi-objective optimization function so as to convert the step multi-energy cooperative multi-objective optimization function into a single-objective optimization problem.

In a further embodiment, the conversion of the step-wise multi-energy collaborative multi-objective optimization function into a single objective optimization problem is specifically expressed as:

Wherein,

In the method, in the process of the invention,The distance from the target point to each point on the Pareto front; Normalizing the point coordinates on the Utobang line; 、 respectively as a single objective function 、Is used for normalizing the objective function value; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution at that time.

In a second aspect, the invention provides a step water storage wind-solar-fire planning operation system, which comprises:

The optimization target building module is used for building a cascade multi-energy cooperative multi-target optimization function by taking the minimum system annual comprehensive cost and the electricity rejection rate as optimization targets and building a cascade water storage wind-solar-fire constraint condition considering the upstream and downstream cascade hydropower constraint;

the planning model construction module is used for linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions to construct a mixed integer linear planning model of cascade water storage wind-light-fire capacity optimization configuration;

The planning model solving module is used for converting the cascade multipotency cooperative multi-objective optimization function into a single-objective optimization problem by utilizing a normal boundary crossing method, solving the mixed integer linear planning model and obtaining a capacity configuration operation scheduling strategy of the cascade multipotency complementary system;

The operation scheduling control module is used for controlling the operation state of the cascade water storage wind-light-fire system according to the capacity configuration operation scheduling strategy of the cascade multi-energy complementary system.

In a third aspect, the present invention also provides a computer device, including a processor and a memory, where the processor is connected to the memory, the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the computer device performs steps for implementing the method.

In a fourth aspect, the present invention also provides a computer readable storage medium having stored therein a computer program which when executed by a processor performs the steps of the above method.

The invention provides a cascade water storage wind-solar-fire planning operation method, a system, equipment and a medium, wherein the method takes the minimum annual comprehensive cost and the electricity rejection rate of a system as optimization targets, establishes a cascade multi-energy cooperative multi-target optimization function, and establishes a cascade water storage wind-solar-fire constraint condition considering the upstream and downstream cascade water storage constraint; linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions, and constructing a mixed integer linear programming model of cascade water storage wind-light-fire capacity optimization configuration; converting the cascade multi-energy collaborative multi-objective optimization function into a single-objective optimization problem by using a normal boundary crossing method, solving a mixed integer linear programming model, and obtaining a capacity configuration operation scheduling strategy of the cascade multi-energy complementary system; and controlling the running state of the cascade water storage wind-light-fire system according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system. Compared with the prior art, the method realizes the optimized operation of the cascade water storage wind-solar-fire system under the influence of pumping storage and thermal power, effectively reduces the waste electric quantity in various power generation modes, improves the operation efficiency and stability of the system, realizes the full utilization and complementation of energy sources, and provides powerful support for the optimization and sustainable development of energy source structures.

Drawings

FIG. 1 is a schematic flow chart of a step water storage wind-solar-fire planning operation method provided by the embodiment of the invention;

FIG. 2 is a block diagram of a step water storage wind-solar-fire planning operation system provided by the embodiment of the invention;

fig. 3 is a schematic structural diagram of a computer device according to an embodiment of the present invention.

Detailed Description

The following examples are given for the purpose of illustration only and are not to be construed as limiting the invention, including the drawings for reference and description only, and are not to be construed as limiting the scope of the invention as many variations thereof are possible without departing from the spirit and scope of the invention.

Referring to fig. 1, an embodiment of the invention provides a step water storage wind-solar-fire planning operation method, as shown in fig. 1, comprising the following steps:

S1, taking the minimum annual comprehensive cost and the electricity rejection rate of the system as optimization targets, establishing a cascade multi-energy cooperative multi-target optimization function, and establishing a cascade water storage wind-light-fire constraint condition considering the upstream and downstream cascade hydropower constraint.

The cascade water storage system is formed by carrying out pumped storage unit transformation on the traditional cascade water storage system, the transformation not only adds the pumping capacity of the system on the original basis, but also enlarges the power generation capacity, in particular, when the capacity of the pumped storage unit arranged in the cascade water storage system is larger, the power rejection rate can be effectively reduced, namely the power waste caused by power grid dispatching or demand fluctuation is reduced, however, the problem of cost rise is caused by blind capacity increase, therefore, the capacity of the pumped storage unit is not larger and better, the embodiment builds a cascade multi-energy cooperative multi-objective optimization function on the basis of comprehensively considering two key factors of the system cost and the power rejection rate, the step water storage wind-solar-fire system comprises three step hydropower stations and a thermal power plant, a wind power plant, a photovoltaic power plant and a step pumping power storage station are newly built, the total assembly machine capacity of the step hydropower station is 436MW, the total assembly machine capacity of the thermal power plant is 432MW, and the step water storage wind-solar-fire system is specifically:

Wherein,

Middle functionThe method comprises the following steps:

wherein F is an objective function; the comprehensive cost is the annual cost of the system; The power rejection rate is the electricity rejection rate; Investment cost for the system year; the annual running cost; the actual output of wind power at the time t of the s-th scene; The actual output of photovoltaic power generation at the time t of the s-th scene; predicting a force coefficient for the force of wind power at the time t of the s-th scene; Predicting a force coefficient for the force of photovoltaic power generation at the time t of the s-th scene; planning and constructing capacity for the wind farm; planning and constructing capacity for the photovoltaic power station; Is the scene number; t is the number of operating hours per scene, e.g., T is 24 hours; Is annual interest rate; the operation life of the pumping and accumulating unit is planned; the unit investment cost of the pumping and accumulating unit is set; planning and constructing capacity for the pumping and accumulating unit; planning operation years for the wind farm; The unit investment cost of the wind farm; planning and constructing capacity for the wind farm; planning the operation life of the photovoltaic power station; the unit investment cost of the photovoltaic power station; planning and constructing capacity for the photovoltaic power station; Probability of occurrence for the s-th scene; the number of the thermal power generating units; the output of the g-th thermal power generating unit is the s-th scene t moment; starting up an operation variable for a g-th thermal power generating unit at a s-th scene t moment; starting up cost for the g-th thermal power generating unit; the shutdown operation variable of the g-th thermal power generating unit is the s-th scene t moment; the shutdown cost of the g-th thermal power generating unit; the water level is the step hydropower level; the number of hydroelectric generating sets of the ith hydropower station; starting operation variables of a jth hydropower unit of an ith hydropower station at a time t of an s-th scene; the j-th station of the i-th hydropower station the starting cost of the hydroelectric generating set; the j-th water discarding cost of the i-th hydropower station; The method comprises the steps that the j-th waste water flow of the i-th hydropower station at the t moment of the s-th scene is obtained; for the time interval, since the unit of flow is m ³/s, it is 3600 seconds; starting operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; The starting cost of the pumping working condition of the pumping and accumulating unit; the method comprises the steps of drawing a starting operation variable of a power generation working condition of a storage unit at the time t of an s-th scene; The starting cost of the power generation working condition of the pumping and accumulating unit is set; a. b and c are coal consumption cost coefficients of the thermal power unit; is a peak shaving cost function; p is the output variable of the thermal power unit; The coal consumption cost function of the thermal power generating unit is as follows; the cost of purchasing the thermal power generating unit is; the rotor fracturing cycle is the time t, and is related to the output P of the thermal power unit; The price of the fuel oil; The oil quantity is calculated; the price of the unit coal; The maximum output value of the thermal power unit is obtained; a first medium limit value for the output of the thermal power generating unit; a second medium limit value for the thermal power unit output; The lower limit value of the output of the thermal power generating unit is set; As a resource consumption rate, it may represent the change in availability of resources (e.g., coal) over time; is a time integral variable; Is a reference time point; d is a risk rate, in the embodiment, D can represent the rate of increase of equipment faults along with time, so that the model can consider the influence of risk factors in the whole time period, the thermal power unit coal consumption cost function adopted in the embodiment is enabled to consider not only static factors of cost but also dynamic factors of time and risk through the resource consumption rate and the risk rate, the influence of uncertainty on the running cost of the thermal power unit is better assessed and dealt with while the power demand is met, the power system is helped to adapt to the changes better, and meanwhile, the complementary relation between thermal power and other energy forms (such as wind power) can be recognized more effectively through the dynamic thermal power unit coal consumption cost function, and the whole energy combination is optimized.

In this embodiment, the step water storage wind-solar-fire constraint conditions include upstream and downstream step water electricity constraint, thermal power operation constraint, wind power and photovoltaic power generation constraint and system power balance constraint; the upstream and downstream cascade hydropower constraints include water balance constraints, hydropower output constraints, flow constraints, water head characteristic constraints, water level-reservoir capacity relation constraints and tail water level-lower drainage flow relation constraints, and as the cascade hydropower constraints after being modified by the pumping and storage unit are different from the system without modification, new pumping and storage constraints are added in the embodiment, namely, the upstream and downstream cascade hydropower constraints further include pumping and storage capacity constraints, pumping and storage start-stop constraints, pumping and storage power upper and lower limit constraints and pumping and storage output constraints, and each constraint condition will be specifically explained below:

The water balance constraint is specifically as follows:

In the method, in the process of the invention, The capacity of the ith level reservoir is the time t of the s-th scene; the interval flow of the ith level reservoir at the time t of the s-th scene; The outlet flow of the ith grade reservoir at the time t of the s-th scene; pumping and accumulating unit flow for the s-th scene t moment; Is a time interval; And The minimum storage capacity and the maximum storage capacity of the ith reservoir are respectively; And Respectively regulating the initial and final reservoir capacity coefficients of the reservoir in year; the capacity of the ith level reservoir is the initial moment of the s-th scene; And the capacity of the ith level reservoir is the time T of the s-th scene.

The water power constraint is specifically as follows:

In the method, in the process of the invention, Is water electric power; Is the density of water; g is gravity acceleration; Generating efficiency of a j-th hydroelectric generating set of the i-th hydropower station; The water head height of the ith hydropower station at the time t of the s-th scene; The power generation flow of the jth hydroelectric generating set of the ith hydropower station at the time t of the jth scene is obtained.

The flow constraint is specifically:

In the method, in the process of the invention, The pumping flow of the pumping and accumulating unit at the s-th scene t moment is used; generating flow of the pumping and accumulating unit at the time t of the s-th scene; the method comprises the steps of (1) setting the drainage flow of an ith hydropower station as an ith hydropower station; generating flow for a jth hydroelectric generating set of an ith hydropower station at a time t of an ith scene; the method comprises the steps of (1) discarding water flow for an ith hydropower station of an ith hydropower station; The state variable of the jth hydropower unit of the ith hydropower station at the time t of the s scene; And The minimum power generation flow and the maximum power generation flow of the jth hydroelectric generating set of the ith hydropower station are respectively; The maximum water discharge rate of the i-th hydropower station.

The head characteristic constraint is specifically:

In the method, in the process of the invention, The water level of the ith level reservoir at the time t of the s-th scene; The tail water level of the ith reservoir at the time t of the s-th scene; a minimum water head for the ith reservoir; is the maximum water head of the ith reservoir.

The water level-reservoir capacity relation constraint is specifically:

In the method, in the process of the invention, Is the relation function of the level of the ith reservoir and the reservoir capacity.

In the method, in the process of the invention, Is the relation function of the tail water level of the ith reservoir and the discharging flow.

The pumping and storage capacity constraint is specifically as follows:

In the method, in the process of the invention, The capacity of the pumping and accumulating unit; The capacity of the pumping and accumulating unit is the maximum; is the minimum pumping and accumulating unit capacity.

In the method, in the process of the invention, The method comprises the steps of drawing and accumulating state variables under the power generation working condition of a unit at the time t of the s-th scene; Drawing and accumulating unit for the s-th scene t moment state variables under pumping conditions; The method comprises the steps of drawing a stopping operation variable of the power generation working condition of the storage unit at the time t of the s-th scene; Stopping operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; starting operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; And (5) drawing a starting operation variable of the power generation working condition of the storage unit at the time t of the s-th scene.

In the method, in the process of the invention, The lower limit power coefficient of the power generation working condition of the pumping and accumulating unit; The actual output of the power generation working condition of the pumping and accumulating unit is the s-th scene t moment; the lower limit power coefficient of the pumping working condition of the pumping and accumulating unit; and (5) the actual output of the pumping working condition of the pumping and accumulating unit at the time t of the s-th scene.

The extraction and storage force constraint is specifically as follows:

The thermal power operation constraint is specifically as follows:

In the method, in the process of the invention, The state constraint of the g-th thermal power generating unit is the s-th scene t moment; the minimum output of the g-th thermal power unit; maximum output of the g-th thermal power unit; the output of the g-th thermal power generating unit is the s-th scene t moment; the minimum start-up time of the g-th thermal power generating unit; the minimum shutdown time of the g-th thermal power generating unit.

In the method, in the process of the invention, The new energy installation capacity lower limit proportionality coefficient; is the total assembly machine quantity of the multi-energy complementary system.

The system power balance constraint is specifically:

S2, linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions, and constructing a mixed integer linear programming model of cascade water storage wind-light-fire capacity optimization configuration.

In this embodiment, the step of linearizing the step multi-energy cooperative multi-objective optimization function and the nonlinear part of the step water storage wind-light-fire constraint condition to construct the mixed integer linear programming model of the step water storage wind-light-fire capacity optimization configuration includes:

Specifically, the portion of the special sequence 2 (Special Order Set of type 2, SOS 2) constraint which can be used for linearization through a piecewise linear method comprises a water level-reservoir capacity relation constraint, a tail water level-lower leakage flow relation constraint and a coal consumption cost function of the thermal power unit, wherein the coal consumption cost function of the thermal power unit can also be used for other piecewise linear methods, and the integral term of the coal consumption cost function of the thermal power unit can be approximately calculated through a numerical integration method such as a trapezoidal method or a Simpson method in the linearization process; for the water power output constraint and the pumping power upper and lower limit constraint, the water power output constraint and the pumping power upper and lower limit constraint can be linearized by adopting a McCormick convex hull relaxation method, and for the convenience of understanding the linearization process, for the portion of linearization by adopting a piecewise linearization method based on an SOS2 constraint, the water level-reservoir capacity relation constraint linearization will be specifically described, and the other two constraints are linearized by adopting the same method, and will not be repeated:

For the constraint of the water level-reservoir capacity relationship, the embodiment adopts a piecewise linear method based on SOS2 constraint for linearization, and the steps are as follows:

Let the water level be And storage capacityThe relation between the water level and the storage capacity can be approximated as a mapping function of K segments, and the linearization range of the water level and the storage capacity is determined, and the storage capacity is assumed to beThe linearization range of (2) is divided into K sections, and the storage capacity isIs of the linearization range ofStorage capacityEach segment of the linearization range ends atCorresponding toThe segment points of (a) can be written asIntroducing auxiliary variablesThe water level and the reservoir capacity may be expressed as:

Introducing 0-1 variable for determining water level and reservoir capacity in the kth section If (if)The water level and the storage capacity are indicated to run in the kth section, and the method is as follows:

。

For the water power output constraint and the pumping power upper and lower limit constraint, the embodiment adopts the McCormick convex envelope relaxation method for linearization, and for the convenience of understanding, the embodiment specifically describes the water power output constraint linearization, and the pumping power upper and lower limit constraint adopts the same method for linearization, and is not repeated:

in the context of the water electrical force constraint, AndFor multiplication of two variables, linearization is needed for solving, and the McCormick convex hull relaxation method is adopted, wherein the linearization process is as follows:

Through the steps, the nonlinear parts such as nonlinear water level-reservoir capacity relation constraint and the like can be linearized, so that the complexity of a model is simplified, the model can be suitable for an efficient mixed integer linear programming solver, the solving efficiency of a subsequent target optimization problem is remarkably improved, the accuracy and reliability of a solving result are ensured, meanwhile, the mixed integer linear programming model is built, integer variables (such as equipment capacity, unit start-stop state and the like) existing in a system are fully considered, and the optimizing result is closer to the actual running condition.

S3, converting the cascade multi-energy cooperative multi-objective optimization function into a single-objective optimization problem by using a normal boundary crossing method, and solving the mixed integer linear programming model to obtain a capacity configuration operation scheduling strategy of the cascade multi-energy complementary system.

In this embodiment, the step of converting the step multi-energy collaborative multi-objective optimization function into the single-objective optimization problem by using the normal boundary crossing method includes:

And randomly selecting a point on the normalized Utobang line by using a normal boundary intersection method (normal boundary intersection, NBI) as a target point, calculating the distance from the target point to each point on the Pareto front according to the normalized objective function value, and searching the Pareto front point with the maximized distance by changing the target point to obtain a Pareto solution set of the step multi-energy collaborative multi-objective optimization function so as to convert the step multi-energy collaborative multi-objective optimization function into a single-objective optimization problem.

In particular, for two single objective functions of a cascade of multiple-energy collaborative multiple-objective optimization functions、Respectively solving two single objective function optimal solutions, wherein the single objective functionTo minimize the annual composite cost of the system, a single objective functionTo minimize the power rejection rate, x is a decision variable; the embodiment can find the single objective function through an optimization algorithm (such as a genetic algorithm, particle swarm optimization, and the like)As a single objective functionThe optimal solution is recorded asAt this time, the corresponding x isLikewise, find a single objective functionAs a single objective functionThe optimal solution is recorded asAt this time, the corresponding x is。

Will be、Respectively substituting two single objective functions, solving to obtainAndObtaining a payment matrix:

connecting two diagonal points in the payment matrix ，) And%，) The resulting Utropina line represents the ideal case where both targets reach the optimal boundary at the same time.

In this embodiment, the objective function values that are separately solved are normalized, so that two single objective functions are compared on the same scale, and for any solution x, the normalized objective functions are as follows:

The normalized payment matrix becomes:

The normalized payment matrix represents the relative position of the two objective functions between the optimal solution and the suboptimal solution, and the embodiment selects any point on the normalized Utoban line by using the NBI algorithm And calculate the point to each point on the Pareto Front (Pareto Front, i.e. the set of all Pareto optimal solutions)，) Distance of (2)To find a Pareto solution set, wherein,+=1 And，By changingAndThe value (hold)+=1) Solving for the maximization distanceFind the distance to makeThe largest point is the Pareto solution under the current NBI algorithm, so that the original multi-objective optimization problem is converted into a single-objective optimization problem through the steps, and each problem tries to find a Pareto solution closest to the Utobang line, namely the maximum distance lambda:

In the method, in the process of the invention, The distance from the target point to each point on the Pareto front; To normalize the point coordinates on the urotobant line, AndThe normalized Uotoban line is a line segment formed by connecting two points (0, 1) and (1, 0) in a first interval;、 respectively as a single objective function 、Is used for normalizing the objective function value; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution; As a single objective function At decision variable ofSingle target optimal solution at that time.

The embodiment is realized by continuously adjustingAndAnd solving for the corresponding maximizationThe method comprises the steps of obtaining a series of Pareto solutions, converting the problems into a series of single-objective optimization problems through an NBI algorithm to solve the original multi-objective optimization problems, establishing a multi-objective function of cascade water storage wind-solar-fire planning and operation by taking the system year comprehensive cost and the electricity rejection rate as targets, linearizing nonlinear parts in a model, calling a Gurobi solver (GurobiOptimization) by adopting Yalmip (Yet Another LMI PARSER) MATLAB tool kit for modeling and solving the convex optimization problems to solve the mixed integer linear programming model, and obtaining a cascade multi-energy complementary system capacity configuration operation scheduling strategy which can comprise specific capacities of wind, light and storage and a cascade water storage wind-fire operation method, optimizing operation modes of cascade water, wind power, photovoltaic power generation, reducing the electricity rejection quantity and realizing the optimized operation of the system.

According to the embodiment, the multi-objective optimization problem can be converted into a series of single-objective optimization problems which are easy to solve through an NBI algorithm or other multi-objective optimization solving strategies, the solving difficulty is reduced, conflicts among different objectives are balanced, the overall optimization of the system is realized, one or more optimal solutions can be selected from Pareto fronts according to different system requirements, specific capacities of wind, light, storage and other energy sources and a running method of step water storage wind, light and fire are extracted, and the system can be ensured to run in an optimized mode.

S4, controlling the running state of the cascade water storage wind-light-fire system according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system.

According to the embodiment, the capacity configuration operation scheduling strategy of the cascade multipotency complementary system can be converted into specific scheduling instructions, such as starting/stopping a specific generator set, adjusting the generated energy and the like, the operation state of the cascade water storage wind-solar-fire system is dynamically adjusted through the scheduling instructions, for example, after each energy subsystem receives the scheduling instructions, corresponding adjustment operations are executed, such as adjusting the delivery flow of a hydropower station, increasing the output of photovoltaic power generation when wind power generation is insufficient, and the charging and discharging strategy of an energy storage power station, so that the relation among multiple targets is balanced to adapt to the change of the power grid demand and the fluctuation of energy supply, the accurate control of the power system is realized, the operation efficiency and the stability of the system are improved, the optimal balance is achieved while the cascade water storage wind-fire system meets the energy demand, the electric quantity is further reduced, the utilization rate of renewable energy sources is improved, and the integral optimization of the system is realized.

In the embodiment, the pumped storage power station is used as a high-quality peak regulation and frequency modulation resource, necessary flexibility and reliability are provided for the system together with thermal power, under the conditions of regulation of the pumped storage power station and support of the thermal power, a mixed integer linear programming model is constructed and solved, and an optimal set with minimum economic cost and minimum electric quantity is searched by combining NBI algorithm and other technologies on the premise of meeting safe and stable operation of the system, so that the operation cost of the system is reduced, the overall operation efficiency and stability of the system are improved, the pumped storage power station can rapidly respond when the output of intermittent renewable energy sources such as wind power, photovoltaic and the like fluctuates, the power fluctuation is stabilized, the stability of power grid frequency and voltage is ensured, and meanwhile, the phenomenon that step water and electricity are forced to be discarded due to regulation requirements is reduced.

The embodiment of the invention provides a step water storage wind-solar-fire planning operation method, which takes the minimum annual comprehensive cost and the electricity rejection rate of a system as optimization targets, establishes a step multi-energy cooperative multi-target optimization function, and establishes a step water storage wind-solar-fire constraint condition considering the upstream and downstream step water storage electricity constraint; linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions, and constructing a mixed integer linear programming model of cascade water storage wind-light-fire capacity optimization configuration; converting the cascade multi-energy collaborative multi-objective optimization function into a single-objective optimization problem by using a normal boundary crossing method, solving a mixed integer linear programming model, and obtaining a capacity configuration operation scheduling strategy of the cascade multi-energy complementary system; and controlling the running state of the cascade water storage wind-light-fire system according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system. Compared with the prior art, the method has the advantages that on the basis of comprehensively considering the pumped storage capacity and the thermal power regulating capacity, the running states of cascade hydroelectric power, wind power and photovoltaic power generation are optimized, the waste electric quantity of a renewable energy power generation mode is effectively reduced, the optimized regulation not only improves the utilization rate of electric power resources, but also promotes the running efficiency of the whole energy system, the running optimization of the system is realized through a comprehensive dispatching strategy, the energy utilization efficiency is remarkably improved, and the running cost is reduced.

It should be noted that, the sequence number of each process does not mean that the execution sequence of each process is determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present application.

In one embodiment, as shown in fig. 2, an embodiment of the present invention provides a step water storage wind-solar-fire planning operation system, which includes:

the optimization target building module 101 is used for building a cascade multi-energy cooperative multi-target optimization function by taking the minimum system annual comprehensive cost and the electricity rejection rate as optimization targets and building a cascade water storage wind-solar-fire constraint condition considering the upstream and downstream cascade hydropower constraint;

The planning model construction module 102 is used for linearizing nonlinear parts of the cascade multi-energy cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions to construct a mixed integer linear planning model of the cascade water storage wind-light-fire capacity optimization configuration;

The planning model solving module 103 is configured to convert the cascade multipotency collaborative multi-objective optimization function into a single-objective optimization problem by using a normal boundary intersection method, and solve the mixed integer linear planning model to obtain a cascade multipotency complementary system capacity configuration operation scheduling strategy;

the operation scheduling control module 104 is configured to control an operation state of the cascade water storage wind-solar-fire system according to the cascade multi-energy complementary system capacity configuration operation scheduling strategy.

The specific limitation of the step water storage wind-light-fire planning operation system can be referred to the limitation of the step water storage wind-light-fire planning operation method, and the description is omitted here. Those of ordinary skill in the art will appreciate that the various modules and steps described in connection with the disclosed embodiments of the application may be implemented in hardware, software, or a combination of both. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.

The embodiment of the invention provides a cascade water storage wind-solar-fire planning operation system, which establishes a cascade multi-energy cooperative multi-objective optimization function through an optimization objective establishment module and establishes a cascade water storage wind-solar-fire constraint condition considering upstream and downstream cascade water storage electricity constraint; linearizing nonlinear parts of the cascade multipotency cooperative multi-objective optimization function and the cascade water storage wind-light-fire constraint conditions through a planning model construction module, and constructing a mixed integer linear planning model of cascade water storage wind-light-fire capacity optimization configuration; converting the cascade multi-energy collaborative multi-objective optimization function into a single-objective optimization problem through a planning model solving module, and solving a mixed integer linear planning model to obtain a capacity configuration operation scheduling strategy of a cascade multi-energy complementary system; the operation scheduling control module is used for controlling the operation state of the cascade water storage wind-light-fire system according to the capacity configuration operation scheduling strategy of the cascade multi-energy complementary system. Compared with the prior art, the system has the advantages that on the basis of comprehensively considering pumped storage capacity and thermal power regulation capacity, the running states of cascade hydroelectric power, wind power and photovoltaic power generation are optimized, the waste electric quantity of a renewable energy power generation mode is effectively reduced, the utilization rate of electric power resources is improved through optimized regulation, the running efficiency of the whole energy system is improved, the running optimization of the system is realized through a comprehensive scheduling strategy, the energy utilization efficiency is remarkably improved, and the running cost is reduced.

FIG. 3 is a diagram of a computer device including a memory, a processor, and a transceiver connected by a bus, according to an embodiment of the present invention; the memory is used to store a set of computer program instructions and data and the stored data may be transferred to the processor, which may execute the program instructions stored by the memory to perform the steps of the above-described method.

Wherein the memory may comprise volatile memory or nonvolatile memory, or may comprise both volatile and nonvolatile memory; the processor may be a central processing unit, a microprocessor, an application specific integrated circuit, a programmable logic device, or a combination thereof. By way of example and not limitation, the programmable logic device described above may be a complex programmable logic device, a field programmable gate array, general purpose array logic, or any combination thereof.

In addition, the memory may be a physically separate unit or may be integrated with the processor.

It will be appreciated by those of ordinary skill in the art that the structure shown in FIG. 3 is merely a block diagram of some of the structures associated with the present inventive arrangements and is not limiting of the computer device to which the present inventive arrangements may be implemented, and that a particular computer device may include more or fewer components than those shown, or may combine some of the components, or have the same arrangement of components.

In one embodiment, an embodiment of the present invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the above-described method.

According to the step water storage wind-light-fire planning operation method, the system, the equipment and the medium, provided by the embodiment of the invention, the step water storage wind-light-fire planning operation method linearizes a step multi-energy cooperative multi-objective optimization function and a nonlinear part of a step water storage wind-light-fire constraint condition, and a mixed integer linear planning model of step water storage wind-light-fire capacity optimization configuration is constructed; converting the cascade multipotency cooperative multi-objective optimization function into a single-objective optimization problem, solving a mixed integer linear programming model, and obtaining a capacity configuration operation scheduling strategy of a cascade multipotency complementary system; the running state of the cascade water storage wind-light-fire system is controlled according to the capacity configuration running scheduling strategy of the cascade multi-energy complementary system, and the running state of cascade water storage wind-light-fire planning running method optimizes the running state of cascade water power, wind power and photovoltaic power generation under the synergistic effect of the pumped storage power station and the thermal power generation so as to reduce the electric quantity and realize the economic, efficient and stable optimized running of the system.

In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present invention, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital subscriber line), or wireless (e.g., infrared, wireless, microwave, etc.). The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., SSD), etc.

Those skilled in the art will appreciate that implementing all or part of the above described embodiment methods may be accomplished by way of a computer program stored on a computer readable storage medium, which when executed, may comprise the steps of embodiments of the methods described above.

The foregoing examples represent only a few preferred embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the application. It should be noted that modifications and substitutions can be made by those skilled in the art without departing from the technical principles of the present application, and such modifications and substitutions should also be considered to be within the scope of the present application. Therefore, the protection scope of the patent of the application is subject to the protection scope of the claims.

Claims

1. The step water storage wind-solar-fire planning operation method is characterized by comprising the following steps of:

the operation scheduling strategy is configured according to the capacity of the cascade multipotent complementary system to control the operation state of the cascade water storage wind-light-fire system;

The step multi-energy cooperative multi-objective optimization function specifically comprises the following steps:

Wherein,

Wherein F is an objective function; the comprehensive cost is the annual cost of the system; The power rejection rate is the electricity rejection rate; Investment cost for the system year; the annual running cost; the actual output of wind power at the time t of the s-th scene; The actual output of photovoltaic power generation at the time t of the s-th scene; predicting a force coefficient for the force of wind power at the time t of the s-th scene; Predicting a force coefficient for the force of photovoltaic power generation at the time t of the s-th scene; planning and constructing capacity for the wind farm; planning and constructing capacity for the photovoltaic power station; Is the scene number; t is the number of operating hours per scene; Is annual interest rate; the operation life of the pumping and accumulating unit is planned; the unit investment cost of the pumping and accumulating unit is set; planning and constructing capacity for the pumping and accumulating unit; planning operation years for the wind farm; The unit investment cost of the wind farm; planning and constructing capacity for the wind farm; planning the operation life of the photovoltaic power station; the unit investment cost of the photovoltaic power station; planning and constructing capacity for the photovoltaic power station; Probability of occurrence for the s-th scene; the number of the thermal power generating units; the output of the g-th thermal power generating unit is the s-th scene t moment; starting up an operation variable for a g-th thermal power generating unit at a s-th scene t moment; starting up cost for the g-th thermal power generating unit; the shutdown operation variable of the g-th thermal power generating unit is the s-th scene t moment; the shutdown cost of the g-th thermal power generating unit; the water level is the step hydropower level; the number of hydroelectric generating sets of the ith hydropower station; starting operation variables of a jth hydropower unit of an ith hydropower station at a time t of an s-th scene; the j-th station of the i-th hydropower station the starting cost of the hydroelectric generating set; the j-th water discarding cost of the i-th hydropower station; The method comprises the steps that the j-th waste water flow of the i-th hydropower station at the t moment of the s-th scene is obtained; Is a time interval; starting operation variables of pumping working conditions of the pumping and accumulating unit at the time t of the s-th scene; The starting cost of the pumping working condition of the pumping and accumulating unit; the method comprises the steps of drawing a starting operation variable of a power generation working condition of a storage unit at the time t of an s-th scene; The starting cost of the power generation working condition of the pumping and accumulating unit is set; a. b and c are the coal consumption cost correlation coefficients of the thermal power unit; is a peak shaving cost function; p is the output variable of the thermal power unit; The coal consumption cost function of the thermal power generating unit is as follows; the cost of purchasing the thermal power generating unit is; Rotor fracturing cycle times at the moment t; The price of the fuel oil; The oil quantity is calculated; the price of the unit coal; The maximum output value of the thermal power unit is obtained; a first medium limit value for the output of the thermal power generating unit; a second medium limit value for the thermal power unit output; The lower limit value of the output of the thermal power generating unit is set; Is the resource consumption rate; is a time integral variable; Is a risk rate; is a reference time point;

The step of converting the cascade multi-energy collaborative multi-objective optimization function into a single-objective optimization problem by using a normal boundary crossing method comprises the following steps of:

2. The step water storage wind-solar-fire planning operation method as claimed in claim 1, wherein: the step water storage wind-solar-fire constraint conditions comprise upstream and downstream step water electricity constraint; the upstream and downstream cascade hydropower constraints comprise water balance constraints, hydropower output constraints, flow constraints, water head characteristic constraints, water level-reservoir capacity relation constraints, tail water level-lower drainage flow relation constraints, pumping storage capacity constraints, pumping storage start-stop constraints, pumping storage power upper and lower limit constraints and pumping storage output constraints;

The water balance constraint is specifically as follows:

the water power constraint is specifically as follows:

The flow constraint is specifically:

The head characteristic constraint is specifically:

the water level-reservoir capacity relation constraint is specifically:

The pumping and storage capacity constraint is specifically as follows:

The extraction and storage force constraint is specifically as follows:

3. The step water storage wind-solar-fire planning operation method as claimed in claim 2, wherein: the step water storage wind-light-fire constraint conditions also comprise thermal power operation constraint, wind power and photovoltaic power generation constraint and system power balance constraint;

the thermal power operation constraint is specifically as follows:

the system power balance constraint is specifically:

4. The running method for step water storage wind-light-fire planning according to claim 3, wherein the step of linearizing the step multi-energy cooperative multi-objective optimization function and the nonlinear part of the step water storage wind-light-fire constraint condition to construct a mixed integer linear programming model of step water storage wind-light-fire capacity optimization configuration comprises the following steps:

5. The running method for the cascade water storage wind-solar-fire planning according to claim 1, wherein the step multi-energy cooperative multi-objective optimization function is converted into a single-objective optimization problem specifically expressed as:

Wherein,

6. A step water storage wind-light-fire planning operation system, characterized in that the step water storage wind-light-fire planning operation method according to any one of claims 1 to 5 is applied, and the system comprises:

7. A computer device, characterized by: comprising a processor and a memory, the processor being connected to the memory, the memory being for storing a computer program, the processor being for executing the computer program stored in the memory to cause the computer device to perform the method of any one of claims 1 to 5.

8. A computer-readable storage medium, characterized by: the computer readable storage medium has stored therein a computer program which, when executed, implements the method of any of claims 1 to 5.