CN113256009B - Two-stage multi-standby configuration method, system and device based on aerothermal inertia - Google Patents
Two-stage multi-standby configuration method, system and device based on aerothermal inertia Download PDFInfo
- Publication number
- CN113256009B CN113256009B CN202110600740.9A CN202110600740A CN113256009B CN 113256009 B CN113256009 B CN 113256009B CN 202110600740 A CN202110600740 A CN 202110600740A CN 113256009 B CN113256009 B CN 113256009B
- Authority
- CN
- China
- Prior art keywords
- standby
- inertia
- time period
- time
- constraint
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012937 correction Methods 0.000 claims abstract description 34
- 238000010248 power generation Methods 0.000 claims abstract description 32
- 239000007789 gas Substances 0.000 claims description 31
- 230000009194 climbing Effects 0.000 claims description 11
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 3
- 239000003345 natural gas Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 230000008878 coupling Effects 0.000 abstract description 7
- 238000010168 coupling process Methods 0.000 abstract description 7
- 238000005859 coupling reaction Methods 0.000 abstract description 7
- 238000005457 optimization Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
Abstract
The application discloses a two-stage multi-standby configuration method, a system and a device based on gas-thermal inertia, and belongs to the field of comprehensive energy systems. According to the two-stage multi-standby configuration method based on the gas-heat inertia, a complex coupling relation among the gas-heat inertia power support size, the continuous input time and the continuous correction time of the comprehensive energy system is considered, a gas-heat inertia standby unified model of the comprehensive energy system is built in combination with the actual situation of the standby configuration of the comprehensive energy system, the gas-heat inertia standby, the power generation side standby and the demand side standby are comprehensively considered, the total cost of the standby configuration of the system in a two-stage scheduling period is minimum as an objective function of the standby configuration of the comprehensive energy system, and meanwhile daily constraints, daily constraints and daily-day-related constraints are included, the standby configuration optimization is carried out on the actual park-level comprehensive energy system, and the operation economy of the system is improved on the premise of ensuring the reliability level of the system.
Description
Technical Field
The application relates to the field of comprehensive energy systems, in particular to a two-stage multi-standby configuration method, system and device based on gas-heat inertia.
Background
In recent years, the world energy crisis and the environmental pressure promote the continuous development of comprehensive energy systems, and the multi-energy complementary characteristics of the comprehensive energy systems can help to improve the energy utilization efficiency and promote the new energy consumption. At the same time, however, the comprehensive energy multipotency coupling increases the complexity of the system operation, so that the problem of optimizing the standby configuration of the comprehensive energy system is necessary to be studied in depth so as to cope with various source load uncertainties of the system operation.
However, current technology focuses on studying large grid backup configurations, while studying comprehensive energy backup configurations involves less. Considering that the large power grid and the comprehensive energy system have the risk of unbalanced supply and demand under actual operation, the reliability of the system operation is reasonably ensured, the economy of the system scheduling operation cost is considered, and the large power grid standby configuration method is necessarily extended to the comprehensive energy system standby configuration. Accordingly, a two-stage multi-standby configuration method, system and apparatus based on aerothermal inertia are presented.
Disclosure of Invention
Aiming at the defects of the prior art, the application provides a two-stage multi-standby configuration method, a system and a device based on gas-heat inertia.
The aim of the application can be achieved by the following technical scheme:
a two-stage multi-standby configuration method based on aerothermal inertia comprises the following steps:
constructing a standby model based on standby configuration and aerothermal inertia factors of the comprehensive energy system; constructing an objective function with the minimum total cost as a target based on the standby configuration cost of the comprehensive energy system;
constructing a day-ahead constraint based on the operational constraint of the integrated energy system; constructing an intra-day constraint based on the prediction error; based on the practical condition of daily endogenous load output and the correlation constraint that the daily transaction reserve capacity is called, constructing the correlation constraint of daily and daily;
and optimizing the standby configuration of the comprehensive energy system according to the standby model, the objective function, the day-ahead constraint, the day-in constraint and the association constraint.
Further, the aero-thermal inertial factors comprise aero-thermal inertial power support size, continuous input time and continuous correction time;
the standby model is as follows:wherein->Is->The gas-heat inertia standby output is in a period of time,for the upper limit of the reserve output of the aero-thermal inertia, +.>Is->Time period of continuous start-up of air-heat inertia standby time, < + >>For the maximum continuous time of the aero-thermal inertia standby, < > for>Representing the functional relation of the maximum continuous input time of the aero-thermal inertia reserve and the reserve output, and +.>Is->Time period air-starting thermal inertia standby continuous correction time, < ->Standby shortest continuous correction time for aero-thermal inertia, < >>And the function relation between the minimum continuous correction time of the aero-thermal inertia standby and the standby output and continuous input time is shown.
Further, the method comprises the steps of,、/>the expression formula of (2) is:
；
wherein,、/>、/>、/>、/>、/>、/>、/>、/>、/>、、/>are constant values determined by the gas-heat inertia characteristics of the comprehensive energy system.
Further, the standby configuration of the integrated energy system comprises a gas inertia standby, a thermal inertia standby, a power generation side standby and a demand side standby;
the objective function is:
；
wherein,is->Time period gas inertia standby cost,/->Is->Time period air inertia reserve price->Is->Time period air inertia standby day-ahead traffic capacity, < >>Is->Natural gas price in period->Is->The actual input capacity in the period of time of the gas inertia standby day, < + >>Is->Time period thermal inertial backup cost,/->Is->Time period thermal inertial reserve price,/->Is->Time period thermal inertia standby day-ahead capacity, < >>Is->Time period hot standby actual call price, +.>Is->Time period thermal inertia standby actual input capacity, +.>Is->Time period power generation side standby cost, < >>Is->Spare price on the generating side of time period->Is->Time period power generation side standby day-ahead traffic capacity, < >>Is->Time period electric quantity price @ and @>Is->Time period generation side standby actual input capacity, +.>Is->Time slot demand side standby cost,/->Is->Time period demand side reserve price,/->Is->Time period demand side spare day before transaction capacity, < ->Is->Time period demand side standby actual call price, +.>Is->The period demand side is spare to actually put into capacity.
Further, the day-ahead constraints include a power balance constraint, an external grid input constraint, an external air grid input constraint, an air inertia standby constraint, a thermal inertia standby constraint, a power generation side standby constraint and a demand side standby constraint;
the power balance constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is->Electric energy distribution coefficient of time period comprehensive energy system, +.>For transformer efficiency, < >>Is->Period external grid input power, < >>Is thatTime period photovoltaic output predictive value,/->For the electric energy production efficiency of the cogeneration unit, < >>Is->Time period external air network input power, < >>Is->Time period electrical load predictive value +.>For electric boiler efficiency>For the heat energy generation efficiency of the cogeneration unit,is->A time period thermal load prediction value;
the external grid input constraints are:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Inputting a lower power limit for an external power grid, +.>Inputting an upper power limit for an external power grid, +.>Inputting power down ramp rate for external power grid, < +.>Time of one period, +.>Inputting power ramp rate for an external power grid;
the external air network input constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Inputting power lower limit for external air network, +.>An upper power limit is input for an external air network, < >>The climbing rate under the input power of the external air network is +.>Inputting power climbing rate for an external air network;
the gas inertia reserve constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is 0-1 variable, ">Representation->Time period for inertial input of qi for standby->Representation->During the period, the air inertia is not put into standby, and the air inertia is not put into standby>Upper limit of reserve output for gas inertia, +.>Is->Time period of starting gas inertia standby continuous input time, +.>For the maximum continuous input time of gas inertia standby +.>Representing the functional relation of the maximum continuous input time of the gas inertia reserve and the reserve output, and +.>Is->Time interval gas-starting inertia standby continuous correction time, +.>Standby shortest continuous correction time for gas inertia, +.>Representing the functional relation between the gas inertia standby shortest continuous correction time and standby output and continuous input time;
the thermal inertia reserve constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is 0-1 variable, ">Representation->Time period thermal inertia put into standby,/->Representation->The thermal inertia is not put into standby in a period of time>For the upper limit of the thermal inertia reserve output, +.>Is->Time period starting thermal inertia standby continuous input time, +.>Is thermal inertiaStandby maximum continuous input time->Representing the function of the maximum continuous input time of the thermal inertia reserve and the reserve output, +.>Is->Time period starting thermal inertia standby continuous correction time, +.>Standby for thermal inertia for the shortest continuous correction time, +.>Representing the functional relation between the thermal inertia standby shortest continuous correction time and standby output and continuous input time;
the power generation side backup constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>For the lower limit of the backup output of the power generation side, +.>For the upper limit of the backup output of the power generation side, +.>For the backup ramp down rate of the power generation side, +.>The backup climbing rate is the power generation side;
the demand side standby constraint is:wherein, the->Is 0 to 1Quantity (S)>Representation->Time period thermal inertia put into standby,/->Representation->The thermal inertia is not put into standby in a period of time>For the lower limit of the demand-side reserve output, +.>For the upper limit of the demand-side reserve output, +.>The longest continuous input time is reserved for the demand side in a period of time.
Further, the prediction error includes a photovoltaic prediction error, an electrical load prediction error, and a thermal load prediction error; the intra-day constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is->Time period system actual photovoltaic prediction error, +.>Is->Time period system actual electrical load prediction error, +.>Is->Time period system actual thermal load prediction error.
Further, the association constraint is:。
in another aspect of the present application, a two-stage multi-standby configuration system based on aerothermal inertia is provided, the system comprising:
the standby configuration module is used for calculating the standby configuration of the comprehensive energy system;
the cost regulation and control module is used for analyzing and regulating and controlling the operation cost of the comprehensive energy system;
the day-ahead constraint module is used for calculating power balance constraint, external power grid input constraint, external air grid input constraint, air inertia standby constraint, thermal inertia standby constraint, power generation side standby constraint and demand side standby constraint;
the solar restraint module is used for calculating photovoltaic prediction errors, electrical load prediction errors and thermal load prediction errors;
the daily constraint module is used for calculating the correlation constraint of the actual condition of daily endogenous load output and the daily transaction reserve capacity for calling;
and the execution module optimizes the standby configuration of the comprehensive energy system based on the standby configuration module, the cost regulation module, the day-ahead constraint module and the day-ahead constraint module.
In a third aspect of the present application, a storage device is provided, in which a plurality of programs are stored, the programs being used for being loaded and executed by a processor to implement the two-stage multi-standby configuration method based on aerothermal inertia according to the first aspect.
A fourth aspect of the present application provides a processing apparatus comprising a processor adapted to execute programs loaded by the processor and executing a two-stage multi-standby configuration method implementing the aerothermal inertia-based method of the first aspect
The application has the beneficial effects that:
the application provides a two-stage multi-standby configuration method based on gas-heat inertia, which fully excavates the inertia characteristics of a gas-heat system, builds a unified model of gas-heat inertia standby of the comprehensive energy system based on complex coupling relations among the gas-heat inertia power support size, continuous input time and continuous correction time of the comprehensive energy system, utilizes the abundant flexibility of the multi-energy coupling of the comprehensive energy system, adds gas-heat inertia standby on the basis of the standby form of the traditional power generation side and the demand side, takes the lowest total cost of the standby configuration of the system in a two-stage scheduling period as an objective function of the standby configuration of the comprehensive energy system, simultaneously brings in daily constraint, daily constraint and daily-day-related constraint, optimizes the standby configuration of the actual park-level comprehensive energy system, and improves the operation economy of the system on the premise of ensuring the reliability level of the system.
Drawings
The application is further described below with reference to the accompanying drawings.
FIG. 1 is a block diagram of a comprehensive energy system of the present application;
FIG. 2 is a graph showing the configuration results of the spare capacity (day before) of the integrated energy system of the present application;
FIG. 3 shows the result of the standby actual output (daily) configuration of the integrated energy system of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
In the description of the present specification, the descriptions of the terms "one embodiment," "example," "specific example," and the like, mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
A two-stage multi-standby configuration method of a comprehensive energy system considering aero-thermal inertia standby comprises the following steps:
taking the complex coupling relation among the size of the gas-heat inertia power support, the continuous input time and the continuous correction time of the comprehensive energy system into consideration, and constructing a comprehensive energy system gas-heat inertia standby unified model by combining the actual situation of standby configuration of the comprehensive energy system;
the longest continuous input time is related to the input standby power, and the shortest continuous correction time is related to the input standby power and the continuous input time. Combining the actual situation of standby configuration of a comprehensive energy system, taking 24 hours as a complete scheduling period, and each time period is 1 hour, constructing a unified model of aero-thermal inertia standby as follows:
wherein,is->Time period of gas-heat inertia standby force, < >>For the upper limit of the reserve output of the aero-thermal inertia, +.>Is->Time period of continuous start-up of air-heat inertia standby time, < + >>Is the airMaximum continuous time of thermal inertia standby, +.>Representing the functional relation of the maximum continuous input time of the aero-thermal inertia reserve and the reserve output, and +.>Is->Time period air-starting thermal inertia standby continuous correction time, < ->Standby shortest continuous correction time for aero-thermal inertia, < >>And the function relation between the minimum continuous correction time of the aero-thermal inertia standby and the standby output and continuous input time is shown. />、/>The description is as follows:
；
wherein,、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>、/>are constant values determined by the gas-heat inertia characteristics of the comprehensive energy system.
The complex coupling relation among the size of the gas-heat inertia power support, the continuous input time and the continuous correction time of the comprehensive energy system is considered, and the gas-heat inertia standby unified model of the comprehensive energy system is constructed, so that the abundant flexibility of the comprehensive energy system can be fully excavated.
Comprehensively considering gas inertia standby, thermal inertia standby, power generation side standby and demand side standby, and taking the lowest total cost of system standby configuration in a two-stage scheduling period as an objective function of comprehensive energy system standby configuration;
objective functionThe method comprises the following steps: />
Wherein,is->Time period gas inertia standby cost,/->Is->Time period air inertia reserve price->Is->Time period air inertia standby day-ahead traffic capacity, < >>Is->Natural gas price in period->Is->The actual input capacity in the period of time of the gas inertia standby day, < + >>Is->Time period thermal inertial backup cost,/->Is->Time period thermal inertial reserve price,/->Is->Time period thermal inertia standby day-ahead capacity, < >>Is->Time period hot standby actual call price, +.>Is->Time period thermal inertia standby actual input capacity, +.>Is->Time period power generation side standby cost, < >>Is->Spare price on the generating side of time period->Is->Time period power generation side standby day-ahead traffic capacity, < >>Is->Time period electric quantity price @ and @>Is->Time period generation side standby actual input capacity, +.>Is->Time slot demand side standby cost,/->Is->Time period demand side reserve price,/->Is->Time period demand side spare day before transaction capacity, < ->Is->Time period demand side standby actual call price, +.>Is->The period demand side is spare to actually put into capacity. And on the basis of the standby forms of the traditional power generation side and the demand side, the gas-heat inertia standby is added, and the rich flexibility of the multi-energy coupling of the comprehensive energy system is fully utilized.
In the day-ahead dispatching stage, respectively considering power balance constraint, external power grid input constraint, external air grid input constraint, air inertia standby constraint, thermal inertia standby constraint, power generation side standby constraint and demand side standby constraint to jointly form day-ahead constraint;
(1) The power balance constraint is:
wherein,is->Electric energy distribution coefficient of time period comprehensive energy system, +.>For transformer efficiency, < >>Is->Period external grid input power, < >>Is->Time period photovoltaic output predictive value,/->For the electric energy production efficiency of the cogeneration unit, < >>Is->Time period external air network input power, < >>Is->Time period electrical load predictive value +.>For electric boiler efficiency>For the heat energy production efficiency of the cogeneration unit, +.>Is->Time period thermal load predictions.
(2) The external grid input constraints are:
wherein,inputting a lower power limit for an external power grid, +.>Inputting an upper power limit for an external power grid, +.>Inputting power down ramp rate for external power grid, < +.>Time of one period, +.>And inputting a power ramp rate for an external power grid.
(3) The external air network input constraint is as follows:
wherein,inputting power lower limit for external air network, +.>An upper power limit is input for an external air network, < >>The climbing rate under the input power of the external air network is +.>And (5) inputting power climbing rate for an external air network.
(4) The gas inertia reserve constraint is:
wherein,is 0-1 variable, ">Representation->Time period for inertial input of qi for standby->Representation->During the period, the air inertia is not put into standby, and the air inertia is not put into standby>Upper limit of reserve output for gas inertia, +.>Is->Time period of starting gas inertia standby continuous input time, +.>For the maximum continuous input time of gas inertia standby +.>Representing the functional relation of the maximum continuous input time of the gas inertia reserve and the reserve output, and +.>Is->Time interval gas-starting inertia standby continuous correction time, +.>Standby shortest continuous correction time for gas inertia, +.>And the functional relation between the gas inertia standby shortest continuous correction time and standby output and continuous input time is shown.
(5) The thermal inertia reserve constraints are:
wherein,is 0-1 variable, ">Representation->Time period thermal inertia put into standby,/->Representation->The thermal inertia is not put into standby in a period of time>For the upper limit of the thermal inertia reserve output, +.>Is->Time period starting thermal inertia standby continuous input time, +.>For maximum continuous time of engagement for thermal inertia, < > for standby>Representing the function of the maximum continuous input time of the thermal inertia reserve and the reserve output, +.>Is->Time period starting thermal inertia standby continuous correction time, +.>Standby for thermal inertia for the shortest continuous correction time, +.>The minimum continuous correction time of thermal inertia backup is expressed as a function of backup output and continuous input time.
(6) The power generation side standby constraint is:
wherein,for the lower limit of the backup output of the power generation side, +.>For the upper limit of the backup output of the power generation side, +.>For the backup ramp down rate of the power generation side, +.>And the power generation side is provided with a standby climbing rate.
(7) Demand side standby constraint:
wherein,is 0-1 variable, ">Representation->Time period thermal inertia put into standby,/->Representation->The thermal inertia is not put into standby in a period of time>For the lower limit of the demand-side reserve output, +.>For the upper limit of the demand-side reserve output, +.>The longest continuous input time is reserved for the demand side in a period of time.
In the daily scheduling stage, the reliability constraint of the actual photovoltaic prediction error, the electrical load prediction error and the thermal load prediction error of the system is considered to be satisfied through the standby configuration of a plurality of standby forms, so as to form the daily constraint;
the daily constraints are:
wherein,is->Time period system actual photovoltaic prediction error, +.>Is->Time period system actual electrical load prediction error, +.>Is->Time period system actual thermal load prediction error.
In the day-ahead and day-ahead associated scheduling stage, the correlation constraint of the comprehensive energy operators for calling the day-ahead bargain capacity according to the actual condition of the day-ahead bargain, namely the bargain capacity of bargain in the day-ahead bargain market, is considered.
Constitutes a pre-day-intra-day association constraint:
comprehensively considering a two-stage multi-standby configuration objective function, a day-ahead constraint and a day-ahead and day-ahead association constraint of the comprehensive energy system, and performing standby configuration optimization on the actual park-level comprehensive energy system. And taking the lowest total cost of the standby configuration of the system in the two-stage scheduling period as an objective function of the standby configuration of the comprehensive energy system, simultaneously taking into the daily constraint, the daily constraint and the daily association constraint, carrying out the standby configuration optimization on the actual park-level comprehensive energy system, and improving the running economy of the system on the premise of ensuring the reliability level of the system.
Examples: taking a park-level integrated energy system as an example, the system structure diagram of the integrated energy system is shown in figure 1. The gas inertia reserve output limit is 1000kw, the output limit is 0kw, the thermal inertia reserve output limit is 800kw, the output limit is 0kw, the power generation side reserve output limit is 1000kw, the output limit is 0kw, the climbing rate is 500 kW/h, the climbing rate is 500 kW/h, the demand side reserve output limit is 800kw, and the output limit is 0kw; the transformer efficiency is 0.98, the electricity generation efficiency of the cogeneration unit is 0.30, the heat energy generation efficiency of the cogeneration unit is 0.40, and the electric boiler efficiency is 0.98. And carrying out 1000 random samples on the actual photovoltaic prediction error, the actual electric load prediction error and the actual thermal load prediction error of the system, and requiring the standby configuration result to meet all random sample uncertainties. The two-stage multi-standby configuration method of the comprehensive energy system taking the aero-thermal inertia standby into consideration is utilized, and the configuration result of the standby transaction capacity (day before) and the configuration result of the standby actual output (day in) are respectively shown in figures 2 and 3.
The result shows that the system can fully excavate the inertia characteristics of the gas-heat system, and under the premise of meeting different constraints of various standby forms, the two-stage coordination optimization configuration of the multi-standby form is carried out so as to fully excavate the flexibility contained in the standby configuration of the comprehensive energy system, improve the running economy of the comprehensive energy system and exert the running advantage of the comprehensive energy system.
The foregoing has shown and described the basic principles, principal features and advantages of the application. It will be understood by those skilled in the art that the present application is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present application, and various changes and modifications may be made without departing from the spirit and scope of the application, which is defined in the appended claims.
Claims (5)
1. The two-stage multi-standby configuration method based on the aerothermal inertia is characterized by comprising the following steps of:
constructing a standby model based on standby configuration and aerothermal inertia factors of the comprehensive energy system; constructing an objective function with the minimum total cost as a target based on the standby configuration cost of the comprehensive energy system;
constructing a day-ahead constraint based on the operational constraint of the integrated energy system; constructing an intra-day constraint based on the prediction error; based on the practical condition of daily endogenous load output and the correlation constraint that the daily transaction reserve capacity is called, constructing the correlation constraint of daily and daily;
optimizing standby configuration of the comprehensive energy system according to the standby model, the objective function, the day-ahead constraint, the day-in constraint and the association constraint;
the standby configuration of the comprehensive energy system comprises gas inertia standby, thermal inertia standby, power generation side standby and demand side standby;
the objective function is:
；
wherein,is->Time period gas inertia standby cost,/->Is->Time period air inertia reserve price->Is->Time period air inertia standby day-ahead traffic capacity, < >>Is->Natural gas price in period->Is->The actual input capacity in the period of time of the gas inertia standby day, < + >>Is->Time period thermal inertial backup cost,/->Is->Time period thermal inertial reserve price,/->Is->Time period thermal inertia standby day-ahead capacity, < >>Is->Time period hot standby actual call price, +.>Is->Time period thermal inertia standby actual input capacity, +.>Is->Time period power generation side standby cost, < >>Is->Spare price on the generating side of time period->Is->Time period power generation side standby day-ahead traffic capacity, < >>Is->Time period electric quantity price @ and @>Is->Time period generation side standby actual input capacity, +.>Is->Time slot demand side standby cost,/->Is->Time period demand side reserve price,/->Is->Time period demand side stand-byCapacity for daily use>Is->Time period demand side standby actual call price, +.>Is->The time period demand side is used for standby actual input capacity;
the day-ahead constraint comprises a power balance constraint, an external power grid input constraint, an external air grid input constraint, an air inertia standby constraint, a thermal inertia standby constraint, a power generation side standby constraint and a demand side standby constraint;
the power balance constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is->Electric energy distribution coefficient of time period comprehensive energy system, +.>For transformer efficiency, < >>Is->Period external grid input power, < >>Is->Time period photovoltaic output predictive value,/->For the electric energy production efficiency of the cogeneration unit, < >>Is->Time period external air network input power, < >>Is thatTime period electrical load predictive value +.>For electric boiler efficiency>For the heat energy production efficiency of the cogeneration unit, +.>Is->A time period thermal load prediction value;
the external grid input constraints are:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Inputting a lower power limit for an external power grid, +.>Inputting an upper power limit for an external power grid, +.>Inputting power down ramp rate for external power grid, < +.>Time of one period, +.>Inputting power ramp rate for an external power grid;
the external air network input constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Inputting power lower limit for external air network, +.>An upper power limit is input for an external air network, < >>The climbing rate under the input power of the external air network is +.>Inputting power climbing rate for an external air network;
the gas inertia reserve constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is 0-1 variable, ">Representation->Time period for inertial input of qi for standby->Representation->During the period, the air inertia is not put into standby, and the air inertia is not put into standby>Upper limit of reserve output for gas inertia, +.>Is->Time period of starting gas inertia standby continuous input time, +.>For the maximum continuous input time of gas inertia standby +.>Representing the functional relation of the maximum continuous input time of the gas inertia reserve and the reserve output, and +.>Is->Time interval gas-starting inertia standby continuous correction time, +.>Standby shortest continuous correction time for gas inertia, +.>Representing the functional relation between the gas inertia standby shortest continuous correction time and standby output and continuous input time;
the thermal inertia reserve constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is 0-1 variable, ">Representation->Time period thermal inertia put into standby,/->Representation->The thermal inertia is not put into standby in a period of time>For the upper limit of the thermal inertia reserve output,is->Time period starting thermal inertia standby continuous input time, +.>The longest continuous input time is reserved for thermal inertia,representing the function of the maximum continuous input time of the thermal inertia reserve and the reserve output, +.>Is->Time period starting thermal inertia standby continuous correction time, +.>Standby for thermal inertia for the shortest continuous correction time, +.>Representing the functional relation between the thermal inertia standby shortest continuous correction time and standby output and continuous input time;
the power generation side backup constraint is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>For the lower limit of the backup output of the power generation side, +.>For the upper limit of the backup output of the power generation side, +.>For the backup ramp down rate of the power generation side, +.>The backup climbing rate is the power generation side;
the prediction errors comprise photovoltaic prediction errors, electrical load prediction errors and thermal load prediction errors; the intra-day constraint is:the method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>Is->Time period system actual photovoltaic prediction error, +.>Is->Time period system actual electrical load prediction error, +.>Is->Prediction error of actual thermal load of a time period system;
the association constraint is:。
2. the two-stage multi-standby configuration method based on aero-thermal inertia of claim 1, wherein the aero-thermal inertia factors include aero-thermal inertia power support size, continuous input time, and continuous correction time;
the standby model is as follows:wherein->Is->Time period of gas-heat inertia standby force, < >>For the upper limit of the reserve output of the aero-thermal inertia, +.>Is->Time period of continuous start-up of air-heat inertia standby time, < + >>For the maximum continuous time of the aero-thermal inertia standby, < > for>Representing the functional relation of the maximum continuous input time of the aero-thermal inertia reserve and the reserve output, and +.>Is->Time period air-starting thermal inertia standby continuous correction time, < ->Standby shortest continuous correction time for aero-thermal inertia, < >>And the function relation between the minimum continuous correction time of the aero-thermal inertia standby and the standby output and continuous input time is shown.
3. The two-stage multi-standby configuration method based on aerothermal inertia according to claim 2, wherein,、the expression formula of (2) is:
；
wherein,、/>、/>、/>、/>、/>、/>、/>、/>、/>、、/>are constant values determined by the gas-heat inertia characteristics of the comprehensive energy system.
4. A storage device in which a plurality of programs are stored, characterized in that the programs are for loading and execution by a processor to implement the two-stage multi-standby configuration method based on aerothermal inertia according to any of claims 1-3.
5. A processing device comprising a processor adapted to execute individual programs, characterized in that the programs are loaded and executed by the processor to implement the two-stage multi-standby configuration method based on aerothermal inertia according to any of claims 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110600740.9A CN113256009B (en) | 2021-05-31 | 2021-05-31 | Two-stage multi-standby configuration method, system and device based on aerothermal inertia |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110600740.9A CN113256009B (en) | 2021-05-31 | 2021-05-31 | Two-stage multi-standby configuration method, system and device based on aerothermal inertia |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113256009A CN113256009A (en) | 2021-08-13 |
| CN113256009B true CN113256009B (en) | 2023-11-28 |
Family
ID=77185461
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110600740.9A Active CN113256009B (en) | 2021-05-31 | 2021-05-31 | Two-stage multi-standby configuration method, system and device based on aerothermal inertia |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113256009B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114723332A (en) * | 2022-05-10 | 2022-07-08 | 广东电网有限责任公司 | Source-load cooperative scheduling method and device under wind power uncertainty |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110197312A (en) * | 2019-06-28 | 2019-09-03 | 东南大学 | A kind of user class integrated energy system Optimization Scheduling based on Multiple Time Scales |
| CN110247392A (en) * | 2019-06-14 | 2019-09-17 | 浙江大学 | More standby resources robust Optimal methods of meter and wind-powered electricity generation marginal capacity and Demand Side Response |
| CN110503271A (en) * | 2019-08-30 | 2019-11-26 | 南京工业大学 | Multi-type energy storage configuration method of comprehensive energy system |
| CN111291942A (en) * | 2020-03-05 | 2020-06-16 | 燕山大学 | Two-stage scheduling method and system of comprehensive energy system considering demand response |
| CN111950171A (en) * | 2020-07-03 | 2020-11-17 | 南京东博智慧能源研究院有限公司 | Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system |
| CN112418488A (en) * | 2020-11-04 | 2021-02-26 | 南方电网科学研究院有限责任公司 | Comprehensive energy system scheduling method and device based on two-stage energy optimization |
| CN112700066A (en) * | 2021-01-14 | 2021-04-23 | 国网山东省电力公司青岛供电公司 | Optimal time scale coordination method for scheduling of electric-thermal integrated energy system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10234886B2 (en) * | 2015-08-04 | 2019-03-19 | Nec Corporation | Management of grid-scale energy storage systems |
-
2021
- 2021-05-31 CN CN202110600740.9A patent/CN113256009B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110247392A (en) * | 2019-06-14 | 2019-09-17 | 浙江大学 | More standby resources robust Optimal methods of meter and wind-powered electricity generation marginal capacity and Demand Side Response |
| CN110197312A (en) * | 2019-06-28 | 2019-09-03 | 东南大学 | A kind of user class integrated energy system Optimization Scheduling based on Multiple Time Scales |
| CN110503271A (en) * | 2019-08-30 | 2019-11-26 | 南京工业大学 | Multi-type energy storage configuration method of comprehensive energy system |
| CN111291942A (en) * | 2020-03-05 | 2020-06-16 | 燕山大学 | Two-stage scheduling method and system of comprehensive energy system considering demand response |
| CN111950171A (en) * | 2020-07-03 | 2020-11-17 | 南京东博智慧能源研究院有限公司 | Backup configuration method for gas-thermal inertia backup participation park comprehensive energy system |
| CN112418488A (en) * | 2020-11-04 | 2021-02-26 | 南方电网科学研究院有限责任公司 | Comprehensive energy system scheduling method and device based on two-stage energy optimization |
| CN112700066A (en) * | 2021-01-14 | 2021-04-23 | 国网山东省电力公司青岛供电公司 | Optimal time scale coordination method for scheduling of electric-thermal integrated energy system |
Non-Patent Citations (2)
| Title |
|---|
| 考虑气热惯性的综合能源系统备用配置方案;孙维佳等;电力系统自动化;第1-12页 * |
| 考虑热惯性的居民楼宇综合能源系统日前运行优化;刘向向;卢婕;严勤;周琪;李昊翔;严强;;可再生能源(第03期);第129-135页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113256009A (en) | 2021-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108154309B (en) | Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity | |
| Datta et al. | A frequency-control approach by photovoltaic generator in a PV–diesel hybrid power system | |
| CN112104007A (en) | Centralized control method for generalized source storage system scheduling | |
| CN105140936B (en) | A kind of micro-capacitance sensor frequency modulation control method for minimizing cost of electricity-generating | |
| CN110991000B (en) | Modeling method for energy hub considering solid oxide fuel cell and electric conversion gas | |
| CN109242350A (en) | The cooling heating and power generation system capacity configuration optimizing method of meter and translatable load | |
| CN111144620A (en) | Electricity-hydrogen comprehensive energy system considering seasonal hydrogen storage and robust planning method thereof | |
| CN107039975A (en) | A kind of distributed energy resource system energy management method | |
| CN109919399A (en) | A kind of integrated energy system economic load dispatching method and system a few days ago | |
| CN108487994A (en) | A kind of micro- energy net composite energy storage system | |
| CN112186755A (en) | Flexible load energy storage modeling method for regional comprehensive energy system | |
| CN110165715B (en) | Method for connecting electric vehicle energy storage type charging station into virtual power plant | |
| JP2010259303A (en) | Distributed power generation system | |
| CN105958537A (en) | Energy conversion system facing energy Internet and optimal control method thereof | |
| CN103226734A (en) | Equivalent rescheduling determining method of wind power acceptation cost | |
| CN111899125B (en) | Optimized modeling operation method, device and medium of comprehensive energy system | |
| CN109286208A (en) | A kind of integrated energy system dispatching method and system | |
| CN116957294A (en) | Scheduling method for virtual power plant to participate in electric power market transaction based on digital twin | |
| CN113256009B (en) | Two-stage multi-standby configuration method, system and device based on aerothermal inertia | |
| CN113659598B (en) | Hybrid energy supply control method and system based on power and gas complementary control | |
| CN117436773B (en) | Independent micro-grid resource planning method and system containing interruptible load | |
| CN108960564A (en) | Supply of cooling, heating and electrical powers type microgrid distribution system dispatching method based on centralization modeling | |
| CN117767401A (en) | Coordinated optimization operation strategy and coordinated control device of wind-solar-hydrogen storage off-grid system | |
| CN117833273A (en) | Power system look-ahead scheduling model solving method considering frequency constraint based on neural network assistance | |
| CN117277444A (en) | New energy base power capacity optimal configuration method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |