CN113505976A - 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
CN113505976A
CN113505976A CN202110692151.8A CN202110692151A CN113505976A CN 113505976 A CN113505976 A CN 113505976A CN 202110692151 A CN202110692151 A CN 202110692151A CN 113505976 A CN113505976 A CN 113505976A
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.)
Granted
Application number
CN202110692151.8A
Other languages
Chinese (zh)
Other versions
CN113505976B (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)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Mathematical Physics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Educational Administration (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Game Theory and Decision Science (AREA)
  • Mathematical Analysis (AREA)
  • Quality & Reliability (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Algebra (AREA)
  • Public Health (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (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 scene 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 N-1 fault situation 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 line of the thermodynamic system; 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 a large area, is lack of sharing, has less 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 and thermal interconnection system, wherein the electric and thermal interconnection system refers to an electric and thermal comprehensive energy system, and the calculation method comprises the following steps:
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 parameters of the power system 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 output of the wind turbine; 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 pipeline flow limits and upper and lower node pressure limits; 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 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.
Preferably, step 1 specifically comprises the following sub-steps:
step S210, determining an objective function, specifically, the maximum electric, gas and heat load value of the electric-thermal interconnection comprehensive energy system is the target, and the maximum electric, gas and heat load value is represented as the minimum value of energy supply capacity under each scene, as shown in formula (1):
F=min(Pelectric+Pgas,electric+Pheat) (1),
wherein F is the maximum energy supply capacity of the required electric-thermal interconnection comprehensive energy system, and P iselectricIs the power supply capacity of the power system under each scene, Pgas,electricIs the natural gas system energy supply capacity under various circumstances, PheatThe 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 RE-GDA0003167092670000031
wherein, PGIs the power generation capacity of the generator, PwIs the generated energy of the wind turbine generator, PCHPIs the electricity quantity produced by the CHP unit, PG2PIs the amount of power produced by the G2P device,
Figure RE-GDA0003167092670000032
and
Figure RE-GDA0003167092670000033
is the discharge and charge quantity of the energy storage element, PP2GIs the power consumption of the P2G device, PLIs an electrical load;
the power flow constraint is expressed as shown in formulas (3) and (4):
Figure RE-GDA0003167092670000034
Figure RE-GDA0003167092670000035
where i and j are power system nodes, PijIs the power flow of branch ij, θiAnd thetajIs the phase angle of nodes i and j;
the phase angle constraint is expressed as shown in equation (5):
Figure RE-GDA0003167092670000036
the generator is constrained as shown in equation (6):
Figure RE-GDA0003167092670000037
wherein, PGiIs the output of generator i;
step S230, constructing natural gas system constraints, including representing the gas balance constraints as shown in equation (7):
Qsource+QP2G=QCHP+QG2P+QL (7),
wherein Q issourceIs the gas output of the natural gas source, QP2GIs P2G plant Natural gas output, QCHPIs the natural gas consumption of the CHP unit, QG2PIs the natural gas consumption, Q, of the G2P plantLIs 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 RE-GDA0003167092670000041
Figure RE-GDA0003167092670000042
Figure RE-GDA0003167092670000043
Figure RE-GDA0003167092670000044
wherein, CmnIs the pipe constant of the pipe mn, fmnIs the natural gas flow of the pipeline mn, TmnIs an auxiliary variable of the pressure difference between the two ends of the pipe mn, pim、πnIs the pressure at the two ends of the pipe mn,
Figure RE-GDA0003167092670000045
is an auxiliary 0-1 binary variable of the pipeline flow direction;
the node pressure constraint is expressed as shown in equation (12):
Figure RE-GDA0003167092670000046
wherein p isiIs the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
Figure RE-GDA0003167092670000047
step S240, constructing thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψCHPsource=ψL (14),
wherein psiCHPIs the heat generated by the CHP unit, #sourceIs heat source generating heat quantity psiLIs 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:
Asm=mq (15),
Bhhf=0 (16),
hf=Km|m| (17),
mmin≤m≤mmax (18),
wherein A issIs a node-branch incidence matrix, m is the pipe flow, m isqIs the node flow, hfIs the pipeline head loss, and 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):
Φ=cpmTs-To (19),
Tend=(Tstart-Ta)×e+Ta (20),
∑moutTout=∑minTin (21),
Figure RE-GDA0003167092670000051
Figure RE-GDA0003167092670000052
where Φ is the node thermoelectric load, cpIs the specific heat capacity of water, Ts、ToIs the temperature of water supply and outlet, Tstart、TendIs the temperature at the beginning and end of the pipeline, TaIs the environmentTemperature, e is the pipeline temperature loss constant, min、mout、Tin、ToutFlow and temperature of the inflow and outflow nodes;
step S250, constructing element constraints, as shown in formula (24):
Figure RE-GDA0003167092670000053
wherein, PwIs 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 for improving the source-load matching degree, and a storage battery model is constructed to include the conditions 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 RE-GDA0003167092670000054
Figure RE-GDA0003167092670000061
Figure RE-GDA0003167092670000062
Figure RE-GDA0003167092670000063
SES,min≤SES≤SES,max (29),
wherein,
Figure RE-GDA0003167092670000064
the charging and discharging states of the energy storage element are respectively 0-1 variable;
Figure RE-GDA0003167092670000065
is an energy storage elementUpper and lower limits of the discharge power of (1);
Figure RE-GDA0003167092670000066
the upper limit and the lower limit of the charging power of the energy storage element;
Figure RE-GDA0003167092670000067
Figure RE-GDA0003167092670000068
the charging and discharging efficiency of the energy storage element; sESIs the stored electric quantity of the energy storage element; sES,max、SES,minThe 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 RE-GDA0003167092670000069
Figure RE-GDA00031670926700000610
wherein Q isCHPIs the natural gas consumption of CHP unit, PCHPIs the electric quantity, psi, produced by the CHP unitCHPIs the heat produced by the CHP unit,
Figure RE-GDA00031670926700000611
is the power generation and heat production efficiency of the CHP unit, qgasIs 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 to generate methane and water under a catalyst on the basis of electricity-to-hydrogen conversion, and the model for constructing the P2G equipment is as shown in formula (32):
ηP2GPP2G=QP2G·qGas (32),
wherein, PP2GIs the power consumption, Q, of the P2G deviceP2GIs P2G plant natural gas output, ηP2GIs the energy conversion efficiency of the P2G plant;
constructing the G2P plant constraint is as shown in equation (33):
PG2P=ηG2PQG2P·qGas (33),
wherein, PG2PIs the amount of electricity, Q, produced by the G2P plantG2PIs the natural gas consumption, eta, of the G2P plantG2PIs 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 interconnection comprehensive energy system under the condition of each line of the electric power system under the fault.
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 derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within 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 heat 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 heat interconnection comprehensive energy system, and can provide decision support for the safe operation of the comprehensive energy system.
Based on this, in the electrical and thermal interconnection system multi-scenario maximum energy supply capacity model considering the energy storage element and the wind turbine generator, the maximum energy supply capacity under the situation that the electrical and thermal coupling comprehensive energy system has the energy storage element and the wind turbine generator is connected is considered, under the normal situation through calculation, under the situation that a power system line N-1 has a fault, under the situation that a natural gas system line N-1 has a fault, and under the situation that a thermodynamic system line N-1 has a fault, the maximum energy supply capacity model under the two situations that the electrical and thermal interconnection comprehensive energy system is disconnected is considered and the electrical and thermal interconnection comprehensive energy system is not considered is finally obtained.
In order 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 situation;
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 disconnection and the non-disconnection 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 constraints;
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 energy system as the target, which is the minimum value of the energy supply capacity in each scenario.
F=min(Pelectric+Pgas,electric+Pheat) (1)
Wherein, F is the maximum energy supply capacity of the required electric-thermal interconnection comprehensive energy system, and comprises three parts: pelectricPower system energy supply capability under each scenario, Pgas,electricEnergy supply capability of natural gas system under various conditions, PheatThermodynamic system energy supply capacity under each scene.
Step S220 is implemented as follows. The electrical balance constraint is shown below.
Figure RE-GDA0003167092670000092
Wherein, PGIs the power generation capacity of the generator, PwIs the generated energy of the wind turbine generator, PCHPIs the electricity quantity produced by the CHP unit, PG2PIs the amount of power produced by the G2P device,
Figure RE-GDA0003167092670000093
and
Figure RE-GDA0003167092670000094
is the discharge and charge quantity of the energy storage element, PP2GIs the power consumption of the P2G device, PLIs an electrical load.
The power flow constraint is shown as follows.
Figure RE-GDA0003167092670000101
Figure RE-GDA0003167092670000102
Where i and j are power system nodes, PijIs the power flow of branch ij, θiAnd thetajIs the phase angle of nodes i and j.
The phase angle constraint is shown below.
Figure RE-GDA0003167092670000103
The generator constraint is as follows.
Figure RE-GDA0003167092670000104
Wherein, PGiIs the output of the generator i.
Step S230 is implemented as follows. The gas balance constraint is shown below.
Qsource+QP2G=QCHP+QG2P+QL (7)
Wherein Q issourceIs the gas output of the natural gas source, QP2GIs P2G plant Natural gas output, QCHPIs the natural gas consumption of the CHP unit, QG2PIs the natural gas consumption, Q, of the G2P plantLIs 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 RE-GDA0003167092670000105
Figure RE-GDA0003167092670000106
Figure RE-GDA0003167092670000107
Figure RE-GDA0003167092670000108
Wherein, CmnIs the pipe constant of the pipe mn, fmnIs the natural gas flow of the pipeline mn, TmnIs an auxiliary variable of the pressure difference between the two ends of the pipe mn, pim、πnIs the pressure at the two ends of the pipe mn,
Figure RE-GDA0003167092670000115
is an auxiliary 0-1 binary variable of the pipeline flow direction.
The nodal pressure constraint is shown below.
Figure RE-GDA0003167092670000111
Wherein p isiIs the pressure at node i.
The natural gas source constraints are shown below.
Figure RE-GDA0003167092670000112
Step S240 is implemented as follows. The thermal equilibrium constraint is shown below.
ψCHPsource=ψL (14)
Wherein psiCHPIs the heat generated by the CHP unit, #sourceIs heat source generating heat quantity psiLIs 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.
Asm=mq (15)
Bhhf=0 (16)
hf=Km|m| (17)
mmin≤m≤mmax (18)
Wherein A issIs a node-branch incidence matrix, m is the pipe flow, m isqIs the node flow, hfIs 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 mixed 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.
Φ=cpmTs-To (19)
Tend=(Tstart-Ta)×e+Ta (20)
∑moutTout=∑minTin (21)
Figure RE-GDA0003167092670000113
Figure RE-GDA0003167092670000114
Where Φ is the node thermoelectric load, cpIs the specific heat capacity of water, Ts、ToIs the temperature of water supply and outlet, Tstart、 TendIs the temperature at the beginning and end of the pipeline, TaIs the ambient temperature, e is the pipeline temperature loss constant, min、mout、 Tin、ToutFlow 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 RE-GDA0003167092670000121
Wherein, PwIs the output of the wind turbine generator.
The storage battery can be discharged in the charging peak-load period in the valley-load period, and the source-load matching degree is improved as an auxiliary means. 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 RE-GDA0003167092670000122
Figure RE-GDA0003167092670000123
Figure RE-GDA0003167092670000124
Figure RE-GDA0003167092670000125
SES,min≤SES≤SES,max (29)
wherein,
Figure RE-GDA0003167092670000126
the charging and discharging states of the energy storage element are respectively 0-1 variable;
Figure RE-GDA0003167092670000127
the upper limit and the lower limit of the discharge power of the energy storage element;
Figure RE-GDA0003167092670000128
the upper limit and the lower limit of the charging power of the energy storage element;
Figure RE-GDA0003167092670000129
Figure RE-GDA00031670926700001210
the charging and discharging efficiency of the energy storage element; sESIs the stored electric quantity of the energy storage element; sES,max、SES,minThe 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 RE-GDA00031670926700001211
Figure RE-GDA00031670926700001212
wherein Q isCHPIs the natural gas consumption of CHP unit, PCHPIs the electric quantity, psi, produced by the CHP unitCHPIs the heat produced by the CHP unit,
Figure RE-GDA0003167092670000131
is the power generation and heat production efficiency of the CHP unit, qgasIs the natural gas heating value.
The P2G equipment firstly generates hydrogen by alkaline electrolysis of water, and on the basis of electricity-to-hydrogen, the hydrogen and carbon dioxide react under a catalyst to generate methane and water, and the model of the equipment is described as follows:
ηP2GPP2G=QP2G·qGas (32)
wherein, PP2GIs the power consumption, Q, of the P2G deviceP2GIs P2G plant natural gas output, ηP2GIs the energy conversion efficiency of the P2G plant.
The G2P device constraint is shown below.
PG2P=ηG2PQG2P·qGas (33)
Wherein, PG2PIs the amount of electricity, Q, produced by the G2P plantG2PIs the natural gas consumption, eta, of the G2P plantG2PIs 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, in the above implementation, the energy storage element and the maximum multi-scenario energy supply capability model of the electrical-thermal interconnection system of the wind turbine generator are taken into consideration, and step S120 may be implemented by the following steps, 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 N-1 fault situation;
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 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.
Further, in the above implementation, the energy storage element and the maximum multi-scenario energy supply capability model of the electrical-thermal interconnection system of the wind turbine generator are taken into consideration, and step S130 may be implemented by the following steps, 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, the feasibility and the effectiveness of the model are verified based on the improved IEEE39 node power system, the 6-node natural gas system and the 6-node thermal system, and the topological diagram is shown in the figure 2. Using MATLAB R2016a simulation, invoking the Gurobi solver, simulation analysis was performed for the following four scenarios:
(1) normal situations (the power system, the natural gas system and the thermodynamic system are all in normal conditions);
the normal situation, i.e. the power system, the natural gas system and the thermodynamic system are all in normal conditions, and the load values are shown in the following table, PloadReferred to as electrical load, QloadIs referred to as the heat load, HloadThe unit is kW, which means the heat load.
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.4222m3And/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 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 scene;
the natural gas system line N-1 fault scene comprises five conditions of node 1 and node 2 disconnection, node 2 and node 4 disconnection, node 3 and node 5 disconnection, node 5 and node 6 disconnection and node 2 and node 5 disconnection, wherein the total electric load, gas load 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 situations 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 considered to be 4441.7197kW when the natural gas system nodes 5 and 6 are disconnected.
TABLE 2 energy supply capability for each scenario
Figure RE-GDA0003167092670000151
Figure RE-GDA0003167092670000161
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 situation of N-1 fault of a power system line, N-1 fault of a natural gas system line and N-1 fault of a thermodynamic system line is calculated under the normal situation, 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 is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: 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 and 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 parameters of the power system 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 output of the wind turbine; 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 pipeline flow limits and upper and lower node pressure limits; 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 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, the maximum electric, gas and heat load value of the electric-thermal interconnection comprehensive energy system is the target, and the maximum electric, gas and heat load value is represented as the minimum value of energy supply capacity under each scene, as shown in formula (1):
F=min(Pelectric+Pgas,electric+Pheat) (1),
wherein F is the maximum energy supply capacity of the required electric-thermal interconnection comprehensive energy system, and P iselectricIs the power supply capacity of the power system under each scene, Pgas,electricIs the natural gas system energy supply capacity under various circumstances, PheatThe 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 RE-FDA0003167092660000021
wherein, PGIs the power generation capacity of the generator, PwIs the generated energy of the wind turbine generator, PCHPIs the electricity quantity produced by the CHP unit, PG2PIs the amount of power produced by the G2P device,
Figure RE-FDA0003167092660000022
and
Figure RE-FDA0003167092660000023
is the discharge and charge quantity of the energy storage element, PP2GIs the power consumption of the P2G device, PLIs an electrical load;
the power flow constraint is expressed as shown in formulas (3) and (4):
Figure RE-FDA0003167092660000024
Figure RE-FDA0003167092660000025
where i and j are power system nodes, PijIs the power flow of branch ij, θiAnd thetajIs 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 RE-FDA0003167092660000031
wherein, PGiIs the output of generator i;
step S230, constructing natural gas system constraints, including representing the gas balance constraints as shown in equation (7):
Qsource+QP2G=QCHP+QG2P+QL (7),
wherein Q issourceIs the gas output of the natural gas source, QP2GIs P2G plant Natural gas output, QCHPIs the natural gas consumption of the CHP unit, QG2PIs the natural gas consumption, Q, of the G2P plantLIs 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 RE-FDA0003167092660000032
Figure RE-FDA0003167092660000033
Figure RE-FDA0003167092660000034
Figure RE-FDA0003167092660000035
wherein, CmnIs the pipe constant of the pipe mn, fmnIs the natural gas flow of the pipeline mn, TmnIs an auxiliary variable of the pressure difference between the two ends of the pipe mn, pim、πnIs the pressure at the two ends of the pipe mn,
Figure RE-FDA0003167092660000036
Figure RE-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 RE-FDA0003167092660000041
wherein p isiIs the pressure at node i;
the natural gas source constraint is expressed as shown in equation (13):
Figure RE-FDA0003167092660000042
step S240, constructing thermodynamic system constraints, including: the thermal equilibrium constraint is expressed as shown in equation (14):
ψCHPsource=ψL (14),
wherein psiCHPIs the heat generated by the CHP unit, #sourceIs heat source generating heat quantity psiLIs 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:
Asm=mq (15),
Bhhf=0 (16),
hf=Km|m| (17),
mmin≤m≤mmax (18),
wherein A issIs a node-branch incidence matrix, m is the pipe flow, m isqIs the node flow, hfIs the pipeline head loss, and 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):
Φ=cpmTs-To (19),
Tend=(Tstart-Ta)×e+Ta (20),
∑moutTout=∑minTin (21),
Figure RE-FDA0003167092660000051
Figure RE-FDA0003167092660000052
where Φ is the node thermoelectric load, cpIs the specific heat capacity of water, TsTo is the temperature of water supply and outlet, Tstart、TendIs the temperature at the beginning and end of the pipeline, TaIs the ambient temperature, e is the pipeline temperature loss constant, min、mout、Tin、ToutFlow and temperature of the inflow and outflow nodes;
step S250, constructing element constraints, as shown in formula (24):
Figure RE-FDA0003167092660000053
wherein, PwIs 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 for improving the source-load matching degree, and a storage battery model is constructed to include the conditions 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 RE-FDA0003167092660000054
Figure RE-FDA0003167092660000055
Figure RE-FDA0003167092660000056
Figure RE-FDA0003167092660000057
SES,min≤SES≤SES,max (29),
wherein,
Figure RE-FDA0003167092660000058
the charging and discharging states of the energy storage element are respectively 0-1 variable;
Figure RE-FDA0003167092660000059
Figure RE-FDA00031670926600000510
the upper limit and the lower limit of the discharge power of the energy storage element;
Figure RE-FDA00031670926600000511
the upper limit and the lower limit of the charging power of the energy storage element;
Figure RE-FDA00031670926600000512
the charging and discharging efficiency of the energy storage element; sESIs the stored electric quantity of the energy storage element; sES,max、SES,minThe 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 RE-FDA0003167092660000061
Figure RE-FDA0003167092660000062
wherein Q isCHPIs the natural gas consumption of CHP unit, PCHPIs the electric quantity, psi, produced by the CHP unitCHPIs the heat produced by the CHP unit,
Figure RE-FDA0003167092660000063
is the power generation and heat production efficiency of the CHP unit, qgasIs 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 to generate methane and water under a catalyst on the basis of electricity-to-hydrogen conversion, and the model for constructing the P2G equipment is as shown in formula (32):
ηP2GPP2G=QP2G·qGas (32),
wherein, PP2GIs the power consumption, Q, of the P2G deviceP2GIs P2G plant natural gas output, ηP2GIs the energy conversion efficiency of the P2G plant;
constructing the G2P plant constraint is as shown in equation (33):
PG2P=ηG2PQG2P·qGas (33),
wherein, PG2PIs the amount of electricity, Q, produced by the G2P plantG2PIs the natural gas consumption, eta, of the G2P plantG2PIs 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 interconnection comprehensive energy system under the condition of each line of the electric power system under the fault.
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 true CN113505976A (en) 2021-10-15
CN113505976B 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
王伟亮 等: "考虑运行约束的区域电力–天然气–热力综合能源系统能量流优化分析", 《中国电机工程学报》 *
胡静哲 等: "计及综合能源服务商参与的电力系统低碳经济调度", 《电网技术》 *

Also Published As

Publication number Publication date
CN113505976B (en) 2022-07-26

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
Chen et al. Adaptive robust day-ahead dispatch for urban energy systems
CN110571789B (en) Electric heating air network three-stage scheduling method based on wind power uncertainty under data driving
CN109034508B (en) Comprehensive energy system robust optimization scheduling method considering electric heating double uncertainty
CN110318877B (en) Wind power consumption method of cogeneration system based on bimodal decoupling
CN109818347B (en) Assessment method for wind power consumption capability of electric power system
CN109510238B (en) Coordinated dispatching unit combination method for efficiently solving hydroelectric power, thermal power and wind power
CN110991857A (en) Method for evaluating wind power consumption capability of electric heating integrated energy system
Zhang et al. Bi-level optimization dispatch of integrated-energy systems with P2G and carbon capture
CN113256152A (en) Reliability evaluation method and system for gas-electricity coupling system based on Bayesian network
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN113268698B (en) New energy consumption capacity optimization method and device of comprehensive energy system and terminal equipment
CN113131513B (en) Method for optimizing operation of electric, thermal and gas conversion system with consideration of carbon emission and storage medium
CN110119855A (en) A kind of pneumoelectric coupled system running optimizatin method
Su et al. An interval based cost-emissions optimization strategy for gas-electricity integrated energy systems under uncertainty and demand response
CN113505976B (en) Multi-scene maximum energy supply capacity calculation method for electric-thermal interconnection system
Liu et al. Influence evaluation of integrated energy system on the unit commitment in power system
CN108717598A (en) Index calculating method temporarily drops in the network voltage of electric-gas interconnection integrated energy system
CN110716429A (en) Control method and device of combined cooling heating and power system, computer and storage medium
CN115081867A (en) Energy storage planning method, device and equipment for comprehensive energy system
CN111598320B (en) Micro-energy system scheduling method
Xiang et al. Planning Research on Electrically Coupled Integrated Energy System Based on ARIMA-LSTM Model
Yang et al. Coordinated Optimal Configuration Method of Hybrid Energy Storage Systems in Energy Internet System
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