CN113505976B - Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system - Google Patents
Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system Download PDFInfo
- Publication number
- CN113505976B CN113505976B CN202110692151.8A CN202110692151A CN113505976B CN 113505976 B CN113505976 B CN 113505976B CN 202110692151 A CN202110692151 A CN 202110692151A CN 113505976 B CN113505976 B CN 113505976B
- Authority
- CN
- China
- Prior art keywords
- supply capacity
- natural gas
- electric
- energy supply
- power
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 166
- 239000003345 natural gas Substances 0.000 claims abstract description 80
- 238000004146 energy storage Methods 0.000 claims abstract description 30
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 239000007789 gas Substances 0.000 claims description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 238000000034 method Methods 0.000 claims description 11
- 229910052739 hydrogen Inorganic materials 0.000 claims description 9
- 239000001257 hydrogen Substances 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 8
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 6
- 230000002457 bidirectional effect Effects 0.000 claims description 6
- 238000007599 discharging Methods 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 3
- 239000001569 carbon dioxide Substances 0.000 claims description 3
- 239000003054 catalyst Substances 0.000 claims description 3
- 238000005868 electrolysis reaction Methods 0.000 claims description 3
- 239000002918 waste heat Substances 0.000 claims 2
- 238000007596 consolidation process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000009133 cooperative interaction Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 230000009044 synergistic interaction Effects 0.000 description 1
Images
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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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/067—Enterprise or organisation modelling
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Mathematical Physics (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a multi-scene maximum energy supply capacity calculation method for an electric and thermal interconnection system considering an energy storage element and a wind turbine generator, which comprises the following steps of: calculating the maximum energy supply capacity of the normal-situation electric thermal coupling system; calculating the maximum energy supply capacity of the electric thermal coupling system under the condition of N-1 fault of the power system line; calculating the maximum energy supply capacity of the electric heat coupling system under the condition of N-1 fault of the natural gas system line; calculating the maximum energy supply capacity of the electric heat coupling system under the condition of N-1 fault of the thermodynamic system line; the maximum multi-scene energy supply capacity of the electric and heat coupling system under two conditions of considering the splitting of the electric and heat interconnection comprehensive energy system and not considering the splitting of the electric and heat interconnection comprehensive energy system.
Description
Technical Field
The invention belongs to the field of maximum energy supply capacity of comprehensive energy systems, relates to a multi-scene maximum energy supply capacity calculation method for an electric and thermal interconnection system, and particularly relates to a multi-scene maximum energy supply capacity calculation method for the electric and thermal interconnection system considering an energy storage element and a wind turbine generator.
Background
The traditional energy source is an N separation supply mode, occupies more land, is lack of sharing, is less in interconnection, is disorderly managed and has high cost. The electric-thermal interconnection comprehensive energy system is an N-in-one supply mode, realizes the cooperative interaction among transverse electric, gas and thermal energy sources and among a plurality of sections of longitudinal sources, nets, loads, storage and utilization, and realizes land consolidation, equipment consolidation, function consolidation and operation consolidation. The electric-thermal interconnection comprehensive energy system makes traditional energy and clean energy draw strong points to make up for shortages through integration and optimization of the system, and realizes synergistic interaction so as to achieve the aims of multi-energy synergistic complementation and energy coordination and mutual assistance, and has important practical significance and long-term strategic significance for promoting energy production and consumption revolution, and constructing a safe, high-efficiency and clean low-carbon energy system.
Object of the Invention
The invention aims to overcome the defects in the prior art and provides a method for calculating the maximum energy supply capacity of an electric-thermal interconnection system in multiple scenes, which considers the access of an energy storage element and a wind generating set, calculates the maximum energy supply capacity under the fault scene of a heating system line N-1 under the fault scene of a power system line N-1 and a natural gas system line N-1 under the normal scene, and finally obtains the maximum energy supply capacity calculation method under the two conditions of considering the disconnection of the electric-thermal interconnection comprehensive energy system and not considering the disconnection of the electric-thermal interconnection comprehensive energy system.
Disclosure of Invention
The invention provides a method for calculating the multi-scene maximum energy supply capacity of an electric-thermal interconnection system, which refers to an electric-thermal comprehensive energy system, and comprises the following steps of:
step 1, calculating the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal condition, and specifically comprising the following steps of: inputting parameters of a power system, natural gas system and thermodynamic system; the power system parameters comprise a unit number, a branch number, branch starting and ending nodes, a node where a generator is located, upper and lower limits of branch power, a node where a wind turbine is located, a node where an energy storage element is located, and upper and lower limits of wind turbine output; the natural gas system parameters comprise natural gas source parameter numbers, pipeline starting and ending nodes, nodes where natural gas sources are located, upper and lower limits of pipeline flow and upper and lower limits of node pressure; the thermodynamic system parameters comprise heat source numbers, pipeline starting and ending nodes, nodes where heat sources are located, upper and lower pipeline flow limits and upper and lower node temperature limits;
considering the constraint conditions of the power system, the natural gas system, the thermodynamic system and the coupling element, and aiming at the maximum energy supply capacity of the power system, the natural gas system and the thermodynamic system, the maximum energy supply capacity of the electrical and thermal interconnection system of the energy storage element and the wind turbine generator is calculated;
the maximum energy supply capacity when the N-1 fault occurs in the power system line, the N-1 fault occurs in the natural gas system line and the N-1 fault occurs in the thermodynamic system line is analyzed, the minimum value under the two conditions of system splitting and non-splitting is selected respectively and taken as the maximum energy supply capacity under the two conditions of system splitting and non-splitting of the electro-thermal interconnection system considering the energy storage element and the wind turbine generator.
Preferably, step 1 comprises in particular the following sub-steps:
step S210, determining an objective function, specifically, determining the maximum electrical, gas and thermal load value of the electrical and thermal interconnection comprehensive energy system as an objective, and representing the maximum electrical, gas and thermal load value as the minimum value of energy supply capacity under each scene, as shown in formula (1):
F=min(P electric +P gas,electric +P heat ) (1),
wherein F is the maximum energy supply capacity of the required electric-thermal interconnection comprehensive energy system, and P is electric Is the power supply capacity of the power system under each scene, P gas,electric Is the natural gas system energy supply capacity under various circumstances, P heat The energy supply capacity of the thermodynamic system under each scene;
step S220, constructing power system constraint, and representing the power system constraint as shown in formula (2):
wherein, P G Is the power generation capacity of the generator, P w Is the generated energy of the wind turbine generator, P CHP Is the electricity quantity produced by the CHP unit, P G2P Is the amount of power produced by the G2P device,andis the discharge and charge quantity of the energy storage element, P P2G Is the power consumption of the P2G device, P L Is an electrical load;
the power flow constraint is expressed as shown in formulas (3) and (4):
wherein i and j are power system nodes, P ij Is the power flow of branch ij, θ i And theta j Is the phase angle of nodes i and j;
the phase angle constraint is expressed as shown in equation (5):
the generator is constrained as shown in equation (6):
wherein, P Gi Is the output of generator i;
step S230, constructing natural gas system constraints, including representing the gas balance constraints as shown in equation (7):
Q source +Q P2G =Q CHP +Q G2P +Q L (7),
wherein Q source Is the gas output of the natural gas source, Q P2G Is P2G plant Natural gas yield, Q CHP Is the natural gas consumption of the CHP unit, Q G2P Is the natural gas consumption, Q, of the G2P plant L Is the natural gas load;
constructing a bidirectional second-order conical convex relaxation model of the natural gas pipeline flow, introducing an auxiliary 0-1 binary variable in the pipeline flow direction into the model, and considering the bidirectional flow condition of natural gas in the pipeline, wherein the model is characterized as shown in formulas (8) to (11):
wherein, C mn Is the pipe constant, f, of the pipe mn mn Is the natural gas flow rate, T, of the pipeline mn mn Is an auxiliary variable of the pressure difference between the two ends of the pipe mn, pi m 、π n Is the pressure at the two ends of the pipe mn,is an auxiliary 0-1 binary variable of the pipeline flow direction;
the node pressure constraint is expressed as shown in equation (12):
wherein p is i Is the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
step S240, constructing thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψ CHP +ψ source =ψ L (14),
wherein psi CHP Is the heat generated by the CHP unit, # source Is heat source generating heat quantity psi L Is the thermal load;
constructing a hydraulic network model, wherein the hydraulic network model is described by a flow continuity equation and a loop pressure drop equation, and the flow continuity equation and the loop pressure drop equation are shown as formulas (15) to (18):
A s m=m q (15),
B h h f =0 (16),
h f =Km|m| (17),
m min ≤m≤m max (18),
wherein, A s Is node-branch incidence matrix, m is pipe flow, m is q Is the node flow, h f Is the pipeline head loss, K is the friction coefficient matrix of each pipeline;
constructing a thermodynamic model, wherein the thermodynamic model is composed of a heat equation, a pipeline temperature drop equation and a node mixed temperature equation, the heat equation reflects the relation between the power of a heat supply network and the temperature and the flow, the pipeline temperature drop equation reflects the heat loss phenomenon of a pipeline, and the node mixed temperature equation reflects the temperature relation before and after hot water mixing and is expressed as formulas (19) to (23):
Φ=c p mT s -T o (19),
T end =(T start -T a )×e+T a (20),
∑m out T out =∑m in T in (21),
where Φ is the node thermoelectric load, c p Is the specific heat capacity of water, T s 、T o Is the temperature of the water supply and the water discharge, T start 、T end Is the temperature at the beginning and end of the pipeline, T a Is the ambient temperature, e is the pipeline temperature loss constant, m in 、m out 、T in 、T out Flow and temperature of the inflow and outflow nodes;
step S250, constructing element constraints, as shown in formula (24):
wherein, P w Is the output of the wind turbine;
the storage battery is charged in the valley-load period and discharged in the peak-load period, as an auxiliary means, the source-load matching degree is improved, and a storage battery model is constructed to include the constraints that the storage battery cannot be charged and discharged at the same time, the upper limit of the charging and discharging power and the upper limit of the capacity are restricted, as shown in formulas (25) to (29):
S ES,min ≤S ES ≤S ES,max (29),
wherein, the first and the second end of the pipe are connected with each other,the charging and discharging states of the energy storage element are respectively 0-1 variable;the upper limit and the lower limit of the discharge power of the energy storage element;the upper limit and the lower limit of the charging power of the energy storage element; the charging and discharging efficiency of the energy storage element is improved; s ES Stored electricity being an energy storage elementAn amount; s ES,max 、S ES,min The upper limit and the lower limit of the stored electric quantity of the energy storage element are set;
when the CHP unit generates power through gas, the waste heat boiler collects waste heat to supply heat energy, so that the energy utilization efficiency of the system is improved, and a thermoelectric coupling CHP model of the CHP unit is described as shown in formulas (30) to (31):
wherein Q is CHP Is the natural gas consumption of the CHP unit, P CHP Is the electric quantity, psi, produced by the CHP unit CHP Is the heat produced by the CHP unit,the power generation and heat production efficiency of the CHP unit q gas Is the natural gas calorific value;
the P2G equipment in the thermodynamic system firstly generates hydrogen through alkaline electrolysis of water, and the hydrogen and carbon dioxide react under a catalyst to generate methane and water on the basis of electricity-to-hydrogen conversion, and the model for constructing the P2G equipment is as shown in formula (32):
η P2G P P2G =Q P2G ·q Gas (32),
wherein, P P2G Is the power consumption, Q, of the P2G device P2G Is the P2G plant natural gas output, η P2G Is the energy conversion efficiency of the P2G plant;
constructing the G2P plant constraint is shown in equation (33):
P G2P =η G2P Q G2P ·q Gas (33),
wherein, P G2P Is the amount of electricity, Q, produced by the G2P plant G2P Is the natural gas consumption, eta, of the G2P plant G2P Is the energy conversion efficiency of the G2P plant;
and calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal situation.
Preferably, the step 2 further comprises the following sub-steps:
step S310, calculating the maximum energy supply capacity of the power system line under the condition of N-1 fault, specifically calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the condition of fault of each line of the power system;
step S320, calculating the maximum energy supply capacity of the natural gas system line under the condition of N-1 fault, and calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the condition of fault of each line of the power system;
and a substep S330, calculating the maximum energy supply capacity of the thermodynamic system line under the condition of N-1 fault, specifically calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnected comprehensive energy system of the electric power system under the condition of fault of each line of the electric power system.
Drawings
Fig. 1 is a flow chart of a multi-scenario maximum energy supply capability model of an electrical and thermal interconnection system considering an energy storage element and a wind turbine generator according to an embodiment of the present invention;
fig. 2 is a topological diagram of an electrical-thermal interconnection integrated energy system according to an embodiment of the present invention;
fig. 3 is a diagram of a total load value of an N-1 fault scenario of the power system according to an embodiment of the present invention;
FIG. 4 is a total load value of a natural gas system N-1 fault scenario according to an embodiment of the present invention;
FIG. 5 is a total load value of a thermodynamic system N-1 in a fault scenario according to an embodiment of the present invention;
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The large number of accesses of the coupling device leads to an improved coupling between the electrical thermal energy sources and to an increased uncertainty. The calculation of the maximum energy supply capacity of the electric and thermal interconnection comprehensive energy system considers the coupling property established among all energy sources through the coupling equipment, measures the load supply capacity of the electric and thermal interconnection comprehensive energy system, and can provide decision support for the safe operation of the comprehensive energy system.
Based on this, according to the multi-scenario maximum energy supply capacity model of the electrical-thermal interconnection system considering the energy storage element and the wind turbine generator, the maximum energy supply capacity model considering that the electrical-thermal coupling comprehensive energy system has the energy storage element and the wind turbine generator is connected is obtained, and the maximum energy supply capacity under the situation of failure of the line N-1 of the thermal power system and the line N-1 of the natural gas system under the situation of failure of the line N-1 of the natural gas system under the normal situation is calculated, and the maximum energy supply capacity model considering the disconnection of the electrical-thermal interconnection comprehensive energy system and the disconnection of the electrical-thermal interconnection comprehensive energy system is finally obtained.
To facilitate understanding of the embodiment, first, a multi-scenario maximum energy supply capability model of an electrical and thermal interconnection system considering an energy storage element and a wind turbine generator disclosed in the embodiment of the present invention is described in detail:
the first embodiment is as follows:
fig. 1 is a flow chart of a multi-scenario maximum energy supply capability model of an electrical and thermal interconnection system considering an energy storage element and a wind turbine generator according to an embodiment of the present invention.
Referring to fig. 1, the multi-scenario maximum energy supply capacity model of the electrical thermal interconnection system with the energy storage element and the wind turbine generator includes the following steps:
step S110, calculating the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal condition;
step S120, calculating the maximum energy supply capacity under the fault situations of the N-1 of the power system line, the N-1 of the natural gas system line and the N-1 of the thermodynamic system line;
and step S130, obtaining the maximum energy supply capacity under the two conditions of considering the splitting and non-splitting of the system.
In the embodiment, considering the electrical thermal interconnection system multi-scenario maximum energy supply capability model of the energy storage element and the wind turbine generator, first, step S110 may be implemented by the following steps, including:
step S210, determining a target function;
step S220, restraining the power system;
step S230, natural gas system constraint;
step S240, constraint of a thermodynamic system;
step S250, element constraint;
specifically, step S210 is implemented by using the maximum electrical, gas and thermal load values of the electrical and thermal interconnection integrated energy system as the target, which is the minimum value of the energy supply capacity in each scenario.
F=min(P electric +P gas,electric +P heat ) (1)
Wherein, F is the maximum energy supply ability of the electric-thermal interconnection comprehensive energy system of asking, and it includes the triplex: p is electric Power system energy supply capability under each scenario, P gas,electric Energy supply capability of natural gas system under various conditions, P heat Thermodynamic system energy supply capacity under each scene.
Step S220 is implemented as follows. The electrical balance constraint is shown below.
Wherein, P G Is the power generation capacity of the generator, P w Is the generated energy of the wind turbine generator, P CHP Is the electricity quantity produced by the CHP unit, P G2P Is the amount of power produced by the G2P plant,andis the energy storage element discharges andamount of charge, P P2G Is the power consumption of the P2G device, P L Is an electrical load.
The power flow constraint is shown in the following equation.
Where i and j are power system nodes, P ij Is the power flow of branch ij, θ i And theta j Is the phase angle of nodes i and j.
The phase angle constraint is shown below.
The generator constraint is as follows.
Wherein, P Gi Is the output of the generator i.
Step S230 is implemented as follows. The gas balance constraint is shown below.
Q source +Q P2G =Q CHP +Q G2P +Q L (7)
Wherein Q source Is the gas output of the natural gas source, Q P2G Is P2G plant Natural gas output, Q CHP Is the natural gas consumption of the CHP unit, Q G2P Is the natural gas consumption, Q, of the G2P plant L Is the natural gas load.
The method for constructing the bidirectional two-order conical convex relaxation model of the natural gas pipeline flow is shown as the formula (8) -11. The model introduces an auxiliary 0-1 binary variable in the flow direction of the pipeline, and can consider the bidirectional flow condition of natural gas in the pipeline.
Wherein, C mn Is the pipe constant, f, of the pipe mn mn Is the natural gas flow of the pipeline mn, T mn Is an auxiliary variable of the pressure difference between the two ends of the mn pipeline, pi m 、π n Is the pressure at the two ends of the pipe mn,is an auxiliary 0-1 binary variable of the flow direction of the pipeline.
The nodal pressure constraint is shown below.
Wherein p is i Is the pressure at node i.
The natural gas source constraints are shown below.
Step S240 is implemented as follows. The thermal equilibrium constraint is shown below.
ψ CHP +ψ source =ψ L (14)
Wherein psi CHP Is heat produced by CHP unitAmount psi source Is the heat produced by the heat source, /) L Is the thermal load.
The hydraulic network belongs to a fluid network, and a model of the hydraulic network is described by a flow continuity equation and a loop pressure drop equation.
A s m=m q (15)
B h h f =0 (16)
h f =Km|m| (17)
m min ≤m≤m max (18)
Wherein A is s Is node-branch incidence matrix, m is pipe flow, m is q Is the node flow, h f Is the pipeline head loss and K is the coefficient of friction matrix for each pipeline.
The thermodynamic model is composed of a heat equation, a pipeline temperature drop equation and a node mixing temperature equation. The heat equation reflects the relationship between the power of the heat supply network and the temperature and the flow, the pipeline temperature drop equation reflects the phenomenon of pipeline heat loss, and the node mixing temperature equation reflects the temperature relationship before and after hot water mixing.
Φ=c p mT s -T o (19)
T end =(T start -T a )×e+T a (20)
∑m out T out =∑m in T in (21)
Where Φ is the node thermoelectric load, c p Is the specific heat capacity of water, T s 、T o Is the temperature of the water supply and the water discharge, T start 、T end Is the temperature at the beginning and end of the pipeline, T a Is the ambient temperature, e is the pipeline temperature loss constant, m in 、m out 、T in 、T out Flow into and out of the node, and temperature.
Step S250 is implemented as follows. The output of the wind turbine generator is not higher than the maximum output limit value, the system does not have enough margin to consume a large amount of wind power in a windy period, and a wind abandon condition may occur, and the constraint is shown as the following formula.
Wherein, P w Is the output of the wind turbine generator.
The storage battery can be charged in the valley-load period and discharged in the peak-load period, and is used as an auxiliary means for improving the source-load matching degree. The model comprises the following constraints that the charging and discharging can not be carried out at the same time, the upper limit of the charging and discharging power and the upper limit of the capacity are restricted:
S ES,min ≤S ES ≤S ES,max (29)
wherein the content of the first and second substances,the variables of the charging state and the discharging state of the energy storage element are respectively 0-1;being energy storage elementsUpper and lower limits of discharge power;the upper limit and the lower limit of the charging power of the energy storage element; the charging and discharging efficiency of the energy storage element; s. the ES Is the stored electric quantity of the energy storage element; s ES,max 、S ES,min The upper limit and the lower limit of the stored electric quantity of the energy storage element.
When the CHP unit generates electricity through gas, the waste heat boiler collects waste heat to supply heat energy, so that the energy utilization efficiency of the system is improved. A common thermoelectric coupling CHP model is described as:
wherein Q is CHP Is the natural gas consumption of CHP unit, P CHP Is the electric quantity, psi, produced by the CHP unit CHP Is the heat produced by the CHP unit,is the power generation and heat production efficiency of the CHP unit, q gas Is the natural gas heating value.
The P2G equipment firstly generates hydrogen by alkaline electrolysis of water, and then the hydrogen and carbon dioxide react under a catalyst to generate methane and water on the basis of electro-transformation of hydrogen, and the model of the equipment is described as follows:
η P2G P P2G =Q P2G ·q Gas (32)
wherein, P P2G Is the power consumption, Q, of the P2G device P2G Is the P2G plant natural gas output, η P2G Is the energy conversion efficiency of the P2G plant.
The G2P device constraint is shown below.
P G2P =η G2P Q G2P ·q Gas (33)
Wherein, P G2P Is the amount of electricity, Q, produced by the G2P plant G2P Is the natural gas consumption, eta, of the G2P plant G2P Is the energy conversion efficiency of the G2P plant.
And calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal situation.
Further, the implementation of the electrical and thermal interconnection system multi-scenario maximum energy supply capability model with the energy storage element and the wind turbine generator is implemented by the following steps of S120, including:
step S310, calculating the maximum energy supply capacity of the power system circuit under the condition of N-1 fault;
step S320, calculating the maximum energy supply capacity of the natural gas system line under the condition of N-1 fault;
step S330, calculating the maximum energy supply capacity of the thermodynamic system line under the condition of N-1 fault;
specifically, step S310 is implemented by calling a Gurobi solver to solve the maximum energy supply capacity of the electrical-thermal interconnection comprehensive energy system in the case of a fault in each line of the power system.
Step S320 is implemented by calling a Gurobi solver to solve the maximum energy supply capacity of the electrical-thermal interconnection comprehensive energy system when each line of the power system has a fault.
Step S330 is realized by calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the condition that each line of the power system has a fault.
Further, the implementation of the electrical and thermal interconnection system multi-scenario maximum energy supply capability model with the energy storage element and the wind turbine generator is implemented by the following steps of S130, including:
due to the fact that faults occur to some lines, partial disconnection of the system can be caused, and the maximum power supply capacity under the two conditions of considering the disconnection of the electric-thermal interconnection comprehensive energy system and not considering the disconnection of the electric-thermal interconnection comprehensive energy system is discussed when the maximum power supply capacity of the final electric-thermal interconnection comprehensive energy system is calculated.
The examples are given below:
for example, based on the feasibility and effectiveness of the improved IEEE39 node power system, the 6 node natural gas system and the 6 node thermodynamic system, a topological graph is shown in the figure 2. Using MATLAB R2016a simulation, invoking a Gurobi solver, and performing simulation analysis for the following four scenarios:
(1) normal scenarios (power system, natural gas system, thermal system are all in normal condition);
the normal situation, namely the power system, the natural gas system and the thermodynamic system are all in normal condition, and each load value is shown in the following table, P load Is an electrical load, Q load Is referred to as the heat load, H load The unit of the heat load is kW.
TABLE 1 Normal situation load values
Pload1 | Pload2 | Pload3 | Pload4 | Pload5 | Pload6 |
271.9114 | 393.2046 | 183.7857 | 379.4632 | 6.2573 | 215.8110 |
Pload7 | Pload8 | Pload9 | Pload10 | Pload11 | Pload12 |
239.1199 | 128.0152 | 496.1179 | 185.0995 | 166.8385 | 224.3704 |
Pload13 | Pload14 | Pload15 | Pload16 | Pload17 | Pload18 |
151.4061 | 93.7948 | 189.6647 | 149.7293 | 191.3248 | 6.2253 |
Pload19 | Qload1 | Qload2 | Hload1 | Hload2 | Hload3 |
715.3201 | 176.2536 | 262.8035 | 26.4823 | 4.5061 | 18.0439 |
The total electric load value obtained from the table is 4387.4589kW, and the total gas load value is 439.0571kW, namely 44.4222m 3 And/h, the total heat load value is 49.0323kW, and the energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal situation is 4875.5492 kW.
(2) An N-1 fault situation occurs in a power system line;
because the electric power system circuit is complicated and various, the first five conditions with the minimum energy supply capacity are selected for analysis, including: nodes 34, 20 are disconnected, nodes 35, 22 are disconnected, nodes 30, 2 are disconnected, nodes 6, 31 are disconnected, and nodes 4, 34 are disconnected. The total electric, gas and heat load values under each condition are shown in figure 3
(3) The natural gas system line has an N-1 fault situation;
the natural gas system line N-1 fault situation comprises five conditions of disconnection of nodes 1 and 2, disconnection of nodes 2 and 4, disconnection of nodes 3 and 5, disconnection of nodes 5 and 6 and disconnection of nodes 2 and 5, wherein the total electric, gas and heat load values under each condition are shown in the attached figure 4.
(4) The thermodynamic system circuit experiences an N-1 fault scenario.
The fault situation of the thermodynamic system line N-1 comprises five conditions of disconnection of the nodes 1 and 2, disconnection of the nodes 2 and 5, disconnection of the nodes 2 and 3, disconnection of the nodes 3 and 4 and disconnection of the nodes 3 and 6, wherein the total electric, gas and heat load values in each condition are shown in the attached figure 5.
The maximum energy supply capacities under all the scenes are compared, and the minimum value is selected as the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system as shown in the following table.
As can be seen from the table, the maximum energy supply capacity of the electrical and thermal interconnection integrated energy system is 4271.3694kW when the natural gas system nodes 2 and 5 are disconnected, but since the heat supply network is disconnected at this time, the maximum energy supply capacity when the electrical and thermal interconnection is not disconnected is 4441.7197kW when the natural gas system nodes 5 and 6 are disconnected.
TABLE 2 energy supply capability for each scenario
According to the multi-scenario maximum energy supply capacity model of the electric and thermal interconnection system considering the energy storage element and the wind turbine generator, the access of the energy storage element and the wind turbine generator is considered, the maximum energy supply capacity under the fault scenario of the line N-1 of the electric power system, the fault scenario of the line N-1 of the natural gas system and the fault scenario of the line N-1 of the thermodynamic system is calculated under the normal scenario, and the maximum energy supply capacity model under the two situations of considering the disconnection of the electric and thermal interconnection comprehensive energy system and not considering the disconnection of the electric and thermal interconnection comprehensive energy system is finally obtained.
The computer program product for the multi-scenario maximum energy supply capability model of the electrical and thermal interconnection system involving the energy storage element and the wind turbine generator, provided by the embodiment of the present invention, includes a computer-readable storage medium storing a program code, where instructions included in the program code may be used to execute the method described in the foregoing method embodiment, and specific implementation may refer to the method embodiment, and details are not described herein.
It can be clearly understood by those skilled in the art that, for convenience and simplicity of description, the specific working process of the system and the apparatus described above may refer to the corresponding process in the foregoing method embodiment, and details are not described herein again.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that the following descriptions are only illustrative and not restrictive, and that the scope of the present invention is not limited to the above embodiments: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (3)
1. The method for calculating the multi-scene maximum energy supply capacity of the electric thermal interconnection system is characterized by comprising the following steps of:
step 1, calculating the maximum energy supply capacity of the electric and thermal interconnection comprehensive energy system under the normal situation, and specifically comprising the following steps: inputting parameters of a power system, natural gas system and thermodynamic system; the power system parameters comprise a unit number, a branch number, branch starting and ending nodes, a node where a generator is located, upper and lower limits of branch power, a node where a wind turbine is located, a node where an energy storage element is located, and upper and lower limits of wind turbine output; the natural gas system parameters comprise natural gas source parameter numbers, pipeline starting and ending nodes, nodes where natural gas sources are located, upper and lower limits of pipeline flow and upper and lower limits of node pressure; the thermodynamic system parameters comprise a heat source number, a pipeline number, pipeline starting and ending nodes, nodes where the heat source is located, upper and lower limits of pipeline flow and upper and lower limits of node temperature;
considering constraint conditions of an electric power system, a natural gas system, a thermodynamic system and a coupling element, and aiming at the maximum energy supply capacity of the electric power system, the natural gas system and the thermodynamic system, the maximum energy supply capacity of an electric and thermal interconnection system considering an energy storage element and a wind turbine generator set is obtained;
step 2, calculating the maximum energy supply capacity under the fault situations of the N-1 line of the power system, the N-1 line of the natural gas system and the N-1 line of the thermodynamic system, and specifically calculating the maximum energy supply capacity when the N-1 fault occurs in the line of the power system, the N-1 fault occurs in the line of the natural gas system and the N-1 fault occurs in the line of the thermodynamic system respectively;
step 3, obtaining the maximum energy supply capacity under two conditions of considering system disconnection and non-disconnection, and specifically comprising the following steps:
the maximum energy supply capacity when an N-1 fault occurs in a power system line, an N-1 fault occurs in a natural gas system line and an N-1 fault occurs in a thermodynamic system line is analyzed, the minimum value under the two conditions of considering system disconnection and non-disconnection is selected respectively and used as the maximum energy supply capacity under the two conditions of considering system disconnection and non-disconnection of an electric and thermal interconnection system of an energy storage element and a wind turbine generator.
2. The electrical thermal interconnection system multi-scenario maximum energy supply capacity calculation method according to claim 1, wherein the step 1 specifically comprises the following sub-steps:
step S210, determining an objective function, specifically, determining the maximum electrical, gas and thermal load value of the electrical and thermal interconnection comprehensive energy system as an objective, and representing the maximum electrical, gas and thermal load value as the minimum value of energy supply capacity under each scene, as shown in formula (1):
F=min(P electric +P gas,electric +P heat ) (1),
wherein F is the maximum energy supply capacity of the required electric-thermal interconnection comprehensive energy system, and P is electric Is the power supply capacity of the power system under each scene, P gas,electric Is the natural gas system energy supply capacity under various circumstances, P heat The energy supply capacity of the thermodynamic system under each scene;
step S220, constructing power system constraint, and representing the power system constraint as shown in formula (2):
wherein, P G Is hairElectric power generation amount, P w Is the generated energy of the wind turbine generator, P CHP Is the electricity quantity produced by the CHP unit, P G2P Is the amount of power produced by the G2P device,andis the discharge and charge amount of the energy storage element, P P2G Is the power consumption of the P2G device, P L Is an electrical load;
the power flow constraint is expressed as shown in formulas (3) and (4):
where i and j are power system nodes, P ij Is the power flow of branch ij, θ i And theta j Is the phase angle of nodes i and j;
the phase angle constraint is expressed as shown in equation (5):
θ i min ≤θ i ≤θ i max (5),
the generator is constrained as shown in equation (6):
wherein, P Gi Is the output of generator i;
step S230, constructing natural gas system constraints, including representing the gas balance constraints as shown in equation (7):
Q source +Q P2G =Q CHP +Q G2P +Q L (7),
wherein Q source Is the gas output of the natural gas source, Q P2G Is P2G plant Natural gas yield, Q CHP Is the natural gas consumption of the CHP unit, Q G2P Is the natural gas consumption, Q, of the G2P plant L Is the natural gas load;
constructing a natural gas pipeline flow bidirectional second-order conical convex relaxation model, introducing an auxiliary 0-1 binary variable in the pipeline flow direction into the model, and considering the bidirectional flow condition of natural gas in a pipeline, wherein the model is characterized by being shown in formulas (8) to (11):
wherein, C mn Is the pipe constant, f, of the pipe mn mn Is the natural gas flow of the pipeline mn, T mn Is an auxiliary variable of the pressure difference between the two ends of the pipe mn, pi m 、π n Is the pressure at the two ends of the pipe mn, is an auxiliary 0-1 binary variable of the pipeline flow direction;
the node pressure constraint is expressed as shown in equation (12):
wherein p is i Is the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
step S240, building thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψ CHP +ψ source =ψ L (14),
wherein psi CHP Is the heat generated by the CHP unit, # s o urce Is heat source generating heat quantity psi L Is the thermal load;
constructing a hydraulic network model, wherein the model is described by a flow continuity equation and a loop pressure drop equation, and the equations (15) to (18) are shown as follows:
A s m=m q (15),
B h h f =0 (16),
h f =Km|m| (17),
m min ≤m≤m max (18),
wherein, A s Is node-branch incidence matrix, m is pipe flow, m is q Is the node flow, h f Is the pipeline head loss, and K is the friction coefficient matrix of each pipeline;
the method comprises the following steps of constructing a thermal model, wherein the thermal model is composed of a heat equation, a pipeline temperature drop equation and a node mixed temperature equation, the heat equation reflects the relation between the power of a heat supply network and the temperature and the flow, the pipeline temperature drop equation reflects the pipeline heat loss phenomenon, and the node mixed temperature equation reflects the temperature relation before and after hot water mixing, and is expressed as formulas (19) - (23):
Φ=c p mT s -T o (19),
T end =(T start -T a )×e+T a (20),
∑m out T out =∑m in T in (21),
where Φ is the node thermoelectric load, c p Is the specific heat capacity of water, T s To is the temperature of water supply and outlet, T start 、T end Is the temperature at the beginning and end of the pipeline, T a Is the ambient temperature, e is the pipeline temperature loss constant, m in 、m out 、T in 、T out Flow and temperature of the inflow and outflow nodes;
step S250, constructing element constraints, as shown in formula (24):
wherein, P w Is the output of the wind turbine;
the storage battery is charged in the valley-load period and discharged in the peak-load period, as an auxiliary means, the source-load matching degree is improved, and a storage battery model is constructed to include the constraints that the storage battery cannot be charged and discharged at the same time, the upper limit of the charging and discharging power and the upper limit of the capacity are restricted, as shown in formulas (25) to (29):
S ES,min ≤S ES ≤S ES,max (29),
wherein, the first and the second end of the pipe are connected with each other,the variables of the charging state and the discharging state of the energy storage element are respectively 0-1; the upper limit and the lower limit of the discharge power of the energy storage element;the upper limit and the lower limit of the charging power of the energy storage element;the charging and discharging efficiency of the energy storage element; s. the ES The stored electric quantity of the energy storage element; s ES,max 、S ES,min The upper limit and the lower limit of the stored electric quantity of the energy storage element;
when the CHP unit generates power through gas, the waste heat boiler collects waste heat to supply heat energy, the energy utilization efficiency of the system is improved, and a thermoelectric coupling CHP model of the CHP unit is described as shown in formulas (30) to (31):
wherein Q is CHP Is the natural gas consumption of the CHP unit, P CHP Is the electric quantity, psi, produced by the CHP unit CHP Is the heat produced by the CHP unit,is the power generation and heat production efficiency of the CHP unit, q gas Is the natural gas calorific value;
the P2G equipment in the thermodynamic system firstly generates hydrogen through alkaline electrolysis of water, and the hydrogen and carbon dioxide react under a catalyst to generate methane and water on the basis of electricity-to-hydrogen conversion, and the model for constructing the P2G equipment is as shown in formula (32):
η P2G P P2G =Q P2G ·q Gas (32),
wherein, P P2G Is the power consumption, Q, of the P2G device P2G Is the P2G plant natural gas output, η P2G Is the energy conversion efficiency of the P2G plant;
constructing the G2P plant constraint is as shown in equation (33):
P G2P =η G2P Q G2P ·q Gas (33),
wherein, P G2P Is the amount of electricity, Q, produced by the G2P plant G2P Is the natural gas consumption, eta, of the G2P plant G2P Is the energy conversion efficiency of the G2P plant;
and calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the normal situation.
3. The electrical thermal interconnection system multi-scenario maximum energy supply capacity calculation method according to claim 1, wherein the step 2 further comprises the following sub-steps:
step S310, calculating the maximum energy supply capacity of the power system line under the condition of N-1 fault, specifically calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the condition of fault of each line of the power system;
step S320, calculating the maximum energy supply capacity of the natural gas system line under the condition of N-1 fault, and calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnection comprehensive energy system under the condition of fault of each line of the power system;
and a substep S330, calculating the maximum energy supply capacity of the thermodynamic system line under the condition of N-1 fault, specifically calling a Gurobi solver to solve the maximum energy supply capacity of the electric-thermal interconnected comprehensive energy system of the electric power system under the condition of fault of each line of the electric power system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110692151.8A CN113505976B (en) | 2021-06-22 | 2021-06-22 | Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110692151.8A CN113505976B (en) | 2021-06-22 | 2021-06-22 | Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113505976A CN113505976A (en) | 2021-10-15 |
CN113505976B true CN113505976B (en) | 2022-07-26 |
Family
ID=78010380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110692151.8A Active CN113505976B (en) | 2021-06-22 | 2021-06-22 | Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113505976B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107730129A (en) * | 2017-10-24 | 2018-02-23 | 重庆大学 | Consider the electrical heat interacted system methods of risk assessment of photo-thermal cogeneration of heat and power and electric boiler |
CN109241662A (en) * | 2018-09-29 | 2019-01-18 | 河海大学 | A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies |
CN109583141A (en) * | 2018-12-29 | 2019-04-05 | 国网江苏省电力有限公司南京供电分公司 | A kind of meter and the electric electric-gas interacted system Unit Combination linear model and system for turning gas coupling |
CN110955954A (en) * | 2019-08-19 | 2020-04-03 | 天津大学 | Optimal load reduction method for layered decoupling electric and thermal comprehensive energy system |
CN112448404A (en) * | 2020-11-19 | 2021-03-05 | 国网经济技术研究院有限公司 | Power distribution network reliability efficiency improvement calculation method under electric-gas-heat interconnection background |
-
2021
- 2021-06-22 CN CN202110692151.8A patent/CN113505976B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107730129A (en) * | 2017-10-24 | 2018-02-23 | 重庆大学 | Consider the electrical heat interacted system methods of risk assessment of photo-thermal cogeneration of heat and power and electric boiler |
CN109241662A (en) * | 2018-09-29 | 2019-01-18 | 河海大学 | A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies |
CN109583141A (en) * | 2018-12-29 | 2019-04-05 | 国网江苏省电力有限公司南京供电分公司 | A kind of meter and the electric electric-gas interacted system Unit Combination linear model and system for turning gas coupling |
CN110955954A (en) * | 2019-08-19 | 2020-04-03 | 天津大学 | Optimal load reduction method for layered decoupling electric and thermal comprehensive energy system |
CN112448404A (en) * | 2020-11-19 | 2021-03-05 | 国网经济技术研究院有限公司 | Power distribution network reliability efficiency improvement calculation method under electric-gas-heat interconnection background |
Non-Patent Citations (2)
Title |
---|
考虑运行约束的区域电力–天然气–热力综合能源系统能量流优化分析;王伟亮 等;《中国电机工程学报》;20171231;第37卷(第24期);7108-7119 * |
计及综合能源服务商参与的电力系统低碳经济调度;胡静哲 等;《电网技术》;20200229;第44卷(第2期);514-520 * |
Also Published As
Publication number | Publication date |
---|---|
CN113505976A (en) | 2021-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Optimal dispatch of zero-carbon-emission micro Energy Internet integrated with non-supplementary fired compressed air energy storage system | |
CN109711601B (en) | Distributed optimal scheduling method and device for electricity-gas-heat integrated energy system | |
CN110620402B (en) | Distributed planning operation joint optimization method and system for electricity-gas hybrid system | |
CN109376428B (en) | Reliability evaluation method, device, equipment and storage medium of comprehensive energy system | |
CN109031952B (en) | Hybrid control method for electricity-gas interconnection comprehensive energy system | |
CN107257132B (en) | Method for constructing comprehensive load reduction model of electricity-gas interconnection system considering wind power abandonment | |
CN109818347B (en) | Assessment method for wind power consumption capability of electric power system | |
CN111598289A (en) | Distributed optimization method of integrated energy system considering LSTM photovoltaic output prediction | |
CN110991857A (en) | Method for evaluating wind power consumption capability of electric heating integrated energy system | |
CN112396238A (en) | Real-time optimized scheduling method for thermoelectric combined system | |
CN113256152A (en) | Reliability evaluation method and system for gas-electricity coupling system based on Bayesian network | |
Zhang et al. | Bi-level optimization dispatch of integrated-energy systems with P2G and carbon capture | |
CN109510238B (en) | Coordinated dispatching unit combination method for efficiently solving hydroelectric power, thermal power and wind power | |
CN113268698B (en) | New energy consumption capacity optimization method and device of comprehensive energy system and terminal equipment | |
CN113011655B (en) | Two-stage random planning-based hybrid energy storage system planning method for community multi-energy system | |
CN113505976B (en) | Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system | |
CN108717598A (en) | Index calculating method temporarily drops in the network voltage of electric-gas interconnection integrated energy system | |
Zong et al. | Two-stage robust optimization of regional integrated energy systems considering uncertainties of distributed energy stations | |
Liu et al. | Influence Evaluation of Integrated Energy System on the Unit Commitment in Power System | |
CN110716429A (en) | Control method and device of combined cooling heating and power system, computer and storage medium | |
CN113131513B (en) | Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium | |
CN115081867A (en) | Energy storage planning method, device and equipment for comprehensive energy system | |
Jiang et al. | Combined heat and power dispatch using simplified district heat flow model | |
Li et al. | Scenario-based unit commitment model considering stochastic wind power | |
Xu et al. | Distributionally Robust Operating Reserve Quantification of Integrated Electricity and Heating System |
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 |