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 PDF

Info

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
Application number
CN202110692151.8A
Other languages
Chinese (zh)
Other versions
CN113505976A (en
Inventor
刘玉成
曹春诚
黄伟
罗永利
王鹏浩
彭勃
聂伟峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hulun Buir Power Supply Company State Grid Inner Mongolia Eastern Electric Power Co ltd
State Grid Corp of China SGCC
North China Electric Power University
Original Assignee
Hulun Buir Power Supply Company State Grid Inner Mongolia Eastern Electric Power Co ltd
State Grid Corp of China SGCC
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hulun Buir Power Supply Company State Grid Inner Mongolia Eastern Electric Power Co ltd, State Grid Corp of China SGCC, North China Electric Power University filed Critical Hulun Buir Power Supply Company State Grid Inner Mongolia Eastern Electric Power Co ltd
Priority to CN202110692151.8A priority Critical patent/CN113505976B/en
Publication of CN113505976A publication Critical patent/CN113505976A/en
Application granted granted Critical
Publication of CN113505976B publication Critical patent/CN113505976B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy 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

Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system
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;
step 2, 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, wherein the maximum energy supply capacity is calculated 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 respectively;
step 3, obtaining the maximum energy supply capacity under the two conditions of system splitting and non-splitting, and specifically comprising the following steps:
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):
Figure GDA0003696022250000031
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,
Figure GDA0003696022250000032
and
Figure GDA0003696022250000033
is 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):
Figure GDA0003696022250000034
Figure GDA0003696022250000035
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):
Figure GDA0003696022250000036
the generator is constrained as shown in equation (6):
Figure GDA0003696022250000037
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):
Figure GDA0003696022250000041
Figure GDA0003696022250000042
Figure GDA0003696022250000043
Figure GDA0003696022250000044
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,
Figure GDA0003696022250000045
is an auxiliary 0-1 binary variable of the pipeline flow direction;
the node pressure constraint is expressed as shown in equation (12):
Figure GDA0003696022250000046
wherein p is i Is the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
Figure GDA0003696022250000047
step S240, constructing thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψ CHPsource =ψ 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),
Figure GDA0003696022250000051
Figure GDA0003696022250000052
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):
Figure GDA0003696022250000053
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):
Figure GDA0003696022250000054
Figure GDA0003696022250000061
Figure GDA0003696022250000062
Figure GDA0003696022250000063
S ES,min ≤S ES ≤S ES,max (29),
wherein, the first and the second end of the pipe are connected with each other,
Figure GDA0003696022250000064
the charging and discharging states of the energy storage element are respectively 0-1 variable;
Figure GDA0003696022250000065
the upper limit and the lower limit of the discharge power of the energy storage element;
Figure GDA0003696022250000066
the upper limit and the lower limit of the charging power of the energy storage element;
Figure GDA0003696022250000067
Figure GDA0003696022250000068
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):
Figure GDA0003696022250000069
Figure GDA00036960222500000610
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,
Figure GDA00036960222500000611
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.
Figure GDA0003696022250000091
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,
Figure GDA0003696022250000092
and
Figure GDA0003696022250000093
is 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.
Figure GDA0003696022250000094
Figure GDA0003696022250000095
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.
Figure GDA0003696022250000096
The generator constraint is as follows.
Figure GDA0003696022250000097
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.
Figure GDA0003696022250000101
Figure GDA0003696022250000102
Figure GDA0003696022250000103
Figure GDA0003696022250000104
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,
Figure GDA0003696022250000105
is an auxiliary 0-1 binary variable of the flow direction of the pipeline.
The nodal pressure constraint is shown below.
Figure GDA0003696022250000106
Wherein p is i Is the pressure at node i.
The natural gas source constraints are shown below.
Figure GDA0003696022250000107
Step S240 is implemented as follows. The thermal equilibrium constraint is shown below.
ψ CHPsource =ψ 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)
Figure GDA0003696022250000111
Figure GDA0003696022250000112
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.
Figure GDA0003696022250000113
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:
Figure GDA0003696022250000121
Figure GDA0003696022250000122
Figure GDA0003696022250000123
Figure GDA0003696022250000124
S ES,min ≤S ES ≤S ES,max (29)
wherein the content of the first and second substances,
Figure GDA0003696022250000125
the variables of the charging state and the discharging state of the energy storage element are respectively 0-1;
Figure GDA0003696022250000126
being energy storage elementsUpper and lower limits of discharge power;
Figure GDA0003696022250000127
the upper limit and the lower limit of the charging power of the energy storage element;
Figure GDA0003696022250000128
Figure GDA0003696022250000129
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:
Figure GDA00036960222500001210
Figure GDA00036960222500001211
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,
Figure GDA00036960222500001212
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
Figure GDA0003696022250000151
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):
Figure FDA0003167092660000021
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,
Figure FDA0003167092660000022
and
Figure FDA0003167092660000023
is 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):
Figure FDA0003167092660000024
Figure FDA0003167092660000025
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):
Figure FDA0003167092660000031
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):
Figure FDA0003167092660000032
Figure FDA0003167092660000033
Figure FDA0003167092660000034
Figure FDA0003167092660000035
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,
Figure FDA0003167092660000036
Figure FDA0003167092660000037
is an auxiliary 0-1 binary variable of the pipeline flow direction;
the node pressure constraint is expressed as shown in equation (12):
Figure FDA0003167092660000041
wherein p is i Is the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
Figure FDA0003167092660000042
step S240, building thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψ CHPsource =ψ 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),
Figure FDA0003167092660000051
Figure FDA0003167092660000052
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):
Figure FDA0003167092660000053
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):
Figure FDA0003167092660000054
Figure FDA0003167092660000055
Figure FDA0003167092660000056
Figure FDA0003167092660000057
S ES,min ≤S ES ≤S ES,max (29),
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003167092660000058
the variables of the charging state and the discharging state of the energy storage element are respectively 0-1;
Figure FDA0003167092660000059
Figure FDA00031670926600000510
the upper limit and the lower limit of the discharge power of the energy storage element;
Figure FDA00031670926600000511
the upper limit and the lower limit of the charging power of the energy storage element;
Figure FDA00031670926600000512
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):
Figure FDA0003167092660000061
Figure FDA0003167092660000062
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,
Figure FDA0003167092660000063
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.
CN202110692151.8A 2021-06-22 2021-06-22 Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system Active CN113505976B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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