CN113515853A - Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation - Google Patents

Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation Download PDF

Info

Publication number
CN113515853A
CN113515853A CN202110636908.1A CN202110636908A CN113515853A CN 113515853 A CN113515853 A CN 113515853A CN 202110636908 A CN202110636908 A CN 202110636908A CN 113515853 A CN113515853 A CN 113515853A
Authority
CN
China
Prior art keywords
heat
electric
energy system
comprehensive energy
node
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.)
Pending
Application number
CN202110636908.1A
Other languages
Chinese (zh)
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.)
Lianyungang Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Lianyungang Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
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 Lianyungang Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical Lianyungang Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Priority to CN202110636908.1A priority Critical patent/CN113515853A/en
Publication of CN113515853A publication Critical patent/CN113515853A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an electric-thermal interconnection comprehensive energy system optimal scheduling method based on a linear equation, which comprises the steps of decoupling a complex annular heat supply network into a plurality of radiation type heat supply networks, converting ring network power flows into a plurality of radiation type heat supply network power flows, then obtaining parameter information of an electric power system and a thermodynamic system, constructing an electric-thermal interconnection comprehensive energy system nonlinear steady-state model based on the parameter information, constructing an electric-thermal interconnection comprehensive energy system rapid power flow calculation model by utilizing the model, obtaining an electric-thermal interconnection comprehensive energy system optimal scheduling model, and completing scheduling optimization of an electric-thermal interconnection comprehensive energy system. The invention provides an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation by using a least square method, the model is simple, the calculated amount is small, the convergence problem does not exist, and the calculation efficiency is improved while the calculation precision is ensured.

Description

Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation
Technical Field
The invention relates to the technical field of operation scheduling and control of an integrated energy system, in particular to an optimal scheduling method of an electric heating interconnection integrated energy system based on a linear equation.
Background
Energy is the basis of human survival and development and is a key factor influencing industrial production and human life quality, so that the reduction of environmental pollution caused by using conventional energy while ensuring sustainable energy supply has become a common concern worldwide. The regional comprehensive energy system adopts the technologies of cogeneration, electric heating technology, gas power generation and the like, takes an electric power grid as a coupling object to fuse a heating power grid and a gas power grid, fully utilizes natural complementary characteristics of various energy forms such as electricity, heat, gas and the like, and realizes the cooperative supply and comprehensive cascade utilization of multi-grid flow.
The current regional comprehensive energy system load flow calculation has many defects, for an electric-heat interconnection comprehensive energy system, a Newton method is usually adopted for solving, but the result is sensitive to an iteration initial value and has the problem of numerical stability, meanwhile, the calculation complexity of a system model is sharply increased along with the increase of the system scale and even possibly exceeds the inherent capability of the existing algorithm, and a new algorithm suitable for a complex model of a large-scale comprehensive energy system is urgently needed to be developed.
For example, a radial heat supply network model of an electric-thermal interconnection comprehensive energy system and a system thereof with the patent number of 201811642220.9 establish a radial heat supply network model of the electric-thermal interconnection comprehensive energy system which can be solved quickly, however, the radial heat supply network model cannot calculate the multi-heat-source heat supply network flow or the annular heat supply network flow and cannot be applied to a large-scale comprehensive energy system. Meanwhile, the radial heat supply network model can only carry out load flow calculation and cannot carry out optimized scheduling.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide an optimal scheduling method of an electric-heat interconnection comprehensive energy system based on a linear equation, which improves the calculation efficiency while ensuring the calculation precision.
The technical solution for realizing the purpose of the invention is as follows:
an electric-heat interconnection comprehensive energy system optimization scheduling method based on a linear equation is characterized by comprising the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
Compared with the prior art, the invention has the remarkable advantages that:
(1) the invention provides an electric-thermal interconnection comprehensive energy system optimization scheduling method based on a linear equation, which comprises the steps of firstly deducing a rapid power flow calculation method of a radiation type heat supply network, then providing a method for decoupling an annular heat supply network into the radiation type heat supply network, and converting the solved annular power flow into a plurality of solved radiation type heat supply network power flows;
(2) the invention has simple model, extremely small calculated amount and no convergence problem
(3) The invention discloses an electric-thermal interconnection comprehensive energy system optimization scheduling method based on a linear equation, which is characterized in that a radiation type heat supply network formula is expanded and deformed, and a least square method is utilized, so that the calculation efficiency is improved while the calculation precision is ensured.
The present invention will be further described with reference to the following detailed description and accompanying drawings.
Drawings
Fig. 1 is a schematic structural diagram of an ideal radial heat network model in an embodiment of the invention.
Fig. 2 is a schematic structural diagram of an actual radial heat network model in an embodiment of the present invention.
FIG. 3 is a schematic flow diagram of adjacent conduits in an embodiment of the present invention.
Fig. 4 is a schematic structural diagram of a two-heat-source radiation type heat supply network model in an embodiment of the invention.
Fig. 5 is a schematic structural diagram of two single-heat-source radiation type heat supply network models in the embodiment of the invention.
Fig. 6 is a schematic structural diagram of an electric-thermal interconnection comprehensive energy system in the embodiment of the invention.
Fig. 7 is a schematic diagram of the voltage amplitude of the 33-node distribution network in the embodiment of the present invention.
Fig. 8 is a schematic diagram of a voltage phase angle of a 33-node distribution network in the embodiment of the present invention.
Fig. 9 is a schematic structural diagram of an integrated energy system of a Bali island according to an embodiment of the present invention.
Detailed Description
An electric-heat interconnection comprehensive energy system optimal scheduling method based on a linear equation comprises the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: according to the parameter information obtained in the step 1 and the step 2, a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system is constructed, and the method specifically comprises the following steps:
Figure BDA0003105535280000031
Figure BDA0003105535280000032
Am＝mq (3)
Bhf＝0 (4)
hf＝Km|m| (5)
Figure BDA0003105535280000033
Figure BDA0003105535280000034
Figure BDA0003105535280000035
Φ＝Cpmq(Ts-To) (9)
Figure BDA0003105535280000036
(∑mout)Tout＝∑(minTin) (11)
Cm＝ΦCHP/PCHP (12)
Cz＝ΔΦ/ΔP＝ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiInjected active power and injected reactive power, U, respectively, for node iiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij＝θi-θj，θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
The equations (1) - (2) are a power grid steady-state model, the equations (3) - (8) are a heat supply network hydraulic model, the equation (3) is a node flow balance equation, the equation (4) is a loop pressure equation, the equation (5) is a pressure head loss equation, the joint equations (6) - (8) can obtain a pipeline resistance coefficient K, the equations (9) - (11) are heat supply network thermal models, the equation (9) is a heat load power equation, the equation (10) is a pipeline temperature drop equation, and the equation (11) is a node power conservation equation; the fixed hot spot ratio in the coupling element is described by equation (12), and the variable heat power ratio is described by equation (13).
And 4, step 4: the method comprises the following steps of constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system, wherein the rapid power flow calculation model comprises a power grid power flow calculation linear model and a heat supply network power flow rapid calculation model, and the method specifically comprises the following steps:
a: the power grid load flow calculation linear model specifically comprises the following steps:
Figure BDA0003105535280000041
order Sij＝Pij+jQijThen, there are:
Figure BDA0003105535280000042
Figure BDA0003105535280000043
let Pij＝Pij_1+Pij_2，Qij＝Qij_1+Qij_2Then, there are:
Figure BDA0003105535280000044
Figure BDA0003105535280000045
Figure BDA0003105535280000046
Figure BDA0003105535280000051
because the branch between two nodes is shorter in the power distribution networkTherefore, the phase angle difference is small, and sin delta is takenij≈δij＝δi-δj，cosδ ij1, and take | Vi|≈1，|V j1, and the equations (17) - (20) include:
Figure BDA0003105535280000052
Figure BDA0003105535280000053
Figure BDA0003105535280000054
Figure BDA0003105535280000055
wherein
Figure BDA0003105535280000056
Thereby obtaining a power grid load flow calculation linear model:
Figure BDA0003105535280000057
Figure BDA0003105535280000058
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1＝rijxij/(rij 2+xij 2)，kij_2＝xij 2/(rij 2+xij 2)，rijAnd xijResistance and electricity of line ij respectivelyAnti, deltaiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjThe voltages at nodes i and j, respectively;
from the equations (25) to (26), the injection power equation of the node is a linear equation with respect to the node voltage V and the phase angle δ, and a set of linear equations can be obtained by using the equations (25) to (26) for the nodes other than the balanced node, and the voltage and the phase angle of each node can be obtained by solving the set of linear equations.
B: the fast calculation model of the heat supply network load flow is as follows:
the pipeline temperature drop equation (10) is Taylor expanded, and a 2-order term is reserved, so that:
Figure BDA0003105535280000059
taking a pipeline with two nodes as an example, the following formula (27) is respectively applied, and the simplification is achieved:
Figure BDA0003105535280000061
applying the mathematical induction method to generalize equation (28) to multiple pipelines is:
Figure BDA0003105535280000062
but the heat load exists in the actual heat supply network pipeline, the pipeline flow m1≠m2≠…≠miTherefore, equation (29) is an approximate equation for solving the flow rate miLet us order
Figure BDA0003105535280000063
Then there are:
Figure BDA0003105535280000064
calculating n1,n2…ni-1When the value is the same, neglecting the heat loss of the pipe network, setting the heat quantity flowing through the pipeline i as
Figure BDA0003105535280000065
When calculating the pipe flow m flowing into the load nodeiIn time, because the temperature change of the head end and the tail end of the pipeline is small, the pipeline can be used for cooling
Figure BDA0003105535280000066
Figure BDA0003105535280000067
Wherein k represents the flow miThe number of the pipeline flowing through, when k is i, niAs for formula (30), 1 may be:
Figure BDA0003105535280000068
wherein, TsSupply water temperature to the node, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiIn order to obtain the heat transfer coefficient of the pipe i,
Figure BDA0003105535280000069
Figure BDA00031055352800000610
n i1, k is the flow miThe number of the pipe through which the fluid flows,
Figure BDA00031055352800000611
for heat energy flowing through the conduit k, CpIs the specific heat capacity of water and phi is the thermal load.
Equation (30) can be abbreviated as: ami 2+Bmi+C＝0 (31)
Since the value of C is generally small, equation (31) can be abbreviated as: ami+B＝0 (32)
By using the formula (32), only one linear equation of unity needs to be solved, and meanwhile, the formula (32) decouples the temperature and the flow, decouples the water supply network and the water return network, so that iteration is not needed, and the numerical stability and the calculation efficiency are improved.
And 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides, and the method specifically comprises the following steps:
step 5-1: splitting the annular heat supply network into a multi-heat-source radiation type heat supply network;
equation of pressure drop in the pipeline is hf＝Km2To loop pipeline
Figure BDA0003105535280000071
Wherein N is the number of loop pipelines, introducing variable x, and using equation
Figure BDA0003105535280000072
Equivalent transformation into
Figure BDA0003105535280000073
Mixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, and the loop can be disconnected at the convergence position of the pipeline flow rates;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in a multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied by a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, and the equivalent method comprises the following steps:
Figure BDA0003105535280000074
Figure BDA0003105535280000075
in the formula, M represents the number of heat sources in the multi-heat-source radiation type heat supply network and is respectively H1,...,HM，
Figure BDA0003105535280000076
Respectively represent a heat source H1,...,HMN represents the number of pipes to be branched from the node, and is respectively denoted by k1,…,kN，
Figure BDA0003105535280000077
Respectively the flow rates of the corresponding pipelines are,
Figure BDA0003105535280000078
respectively the flow of the pipeline
Figure BDA0003105535280000079
Corresponding pipeline flow in a split 1 st, … th and M single heat source radiation type heat network model;
Figure BDA00031055352800000710
for the first heat load in a multi-source radiant network,
Figure BDA00031055352800000711
are respectively as
Figure BDA00031055352800000712
At split 1 st …, M single heat source radiative heat network models correspond to the heat load at the load nodes.
Step 6: determining a heat supply network power flow linear optimization model according to the decoupling result in the step 5, thereby obtaining an electric-thermal interconnection comprehensive energy system optimization scheduling model, and outputting state quantity information of the electric-thermal interconnection comprehensive energy system, wherein the method specifically comprises the following steps:
electric-thermal interconnection based rapid comprehensive energy systemAm in load flow calculation modeli 2+Bmi+ C ═ 0:
Figure BDA0003105535280000081
omitting the quadratic term to obtain:
Figure BDA0003105535280000082
the method is simplified to obtain:
Figure BDA0003105535280000083
in the above formula, the variable is generally the heat source temperature THAnd heat load return water temperature ToFor practical systems, when determining the temperature T of the heat sourceHAnd heat load return water temperature ToBy least square method, can
Figure BDA0003105535280000084
To be K1TH+K2To+K3Will be
Figure BDA0003105535280000085
To be K4TH+K5To+K6The fitting coefficient can be generally up to 0.99, and the above equation can be:
Figure BDA0003105535280000086
wherein, K1，K2，K3，K4，K5，K6Is a constant.
And obtaining a heat supply network power flow linear optimization model in the electric-thermal interconnection comprehensive energy system optimization scheduling model.
And (5) combining the formula (25) and the formula (26) to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation.
An electricity-heat interconnection comprehensive energy system optimization scheduling system based on a linear equation is characterized by comprising the following modules:
the electric-thermal interconnection comprehensive energy system nonlinear steady-state model building module comprises: the method is used for constructing a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system based on parameter information of the electric power and thermodynamic system;
the electric-thermal interconnection comprehensive energy system rapid load flow calculation model construction module comprises: the method is used for constructing a rapid load flow calculation model of the electric-thermal interconnection comprehensive energy system;
a heat supply network decoupling module: for decoupling the looped heat network into a plurality of radiant heat networks;
the electric-thermal interconnection comprehensive energy system optimization scheduling model building module comprises: and constructing an optimized dispatching model of the electric-thermal interconnection comprehensive energy system based on the decoupling result and the quick load flow calculation model of the electric-thermal interconnection comprehensive energy system.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the computer program:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
A computer-storable medium having stored thereon a computer program which, when executed by a processor, performs the steps of:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system. Examples
An electric-heat interconnection comprehensive energy system optimal scheduling method based on a linear equation comprises the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: according to the parameter information obtained in the step 1 and the step 2, a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system is constructed, and the method specifically comprises the following steps:
Figure BDA0003105535280000101
Figure BDA0003105535280000102
Am＝mq (3)
Bhf＝0 (4)
hf＝Km|m| (5)
Figure BDA0003105535280000103
Figure BDA0003105535280000104
Figure BDA0003105535280000105
Φ＝Cpmq(Ts-To) (9)
Figure BDA0003105535280000106
(∑mout)Tout＝∑(minTin) (11)
Cm＝ΦCHP/PCHP (12)
Cz＝ΔΦ/ΔP＝ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiRespectively being nodes iInjecting active power and injecting reactive power, UiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij＝θi-θj，θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
The equations (1) - (2) are a power grid steady-state model, the equations (3) - (8) are a heat supply network hydraulic model, the equation (3) is a node flow balance equation, the equation (4) is a loop pressure equation, the equation (5) is a pressure head loss equation, the joint equations (6) - (8) can obtain a pipeline resistance coefficient K, the equations (9) - (11) are heat supply network thermal models, the equation (9) is a heat load power equation, the equation (10) is a pipeline temperature drop equation, and the equation (11) is a node power conservation equation; the fixed hot spot ratio in the coupling element is described by equation (12), and the variable heat power ratio is described by equation (13).
And 4, step 4: the method comprises the following steps of constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system, wherein the rapid power flow calculation model comprises a power grid power flow calculation linear model and a heat supply network power flow rapid calculation model, and the method specifically comprises the following steps:
a: the power grid load flow calculation linear model specifically comprises the following steps:
Figure BDA0003105535280000111
order Sij＝Pij+jQijThen, there are:
Figure BDA0003105535280000112
Figure BDA0003105535280000113
let Pij＝Pij_1+Pij_2，Qij＝Qij_1+Qij_2Then, there are:
Figure BDA0003105535280000114
Figure BDA0003105535280000121
Figure BDA0003105535280000122
Figure BDA0003105535280000123
in the power distribution network, because the branch between two nodes is short, the phase angle difference is small, and sin delta is takenij≈δij＝δi-δj，cosδij≈1，And take | Vi|≈1，|V j1, and the equations (17) - (20) include:
Figure BDA0003105535280000124
Figure BDA0003105535280000125
Figure BDA0003105535280000126
Figure BDA0003105535280000127
wherein
Figure BDA0003105535280000128
Thereby obtaining a power grid load flow calculation linear model:
Figure BDA0003105535280000129
Figure BDA00031055352800001210
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1＝rijxij/(rij 2+xij 2)，kij_2＝xij 2/(rij 2+xij 2)，rijAnd xijRespectively the resistance and reactance, delta, of the line ijiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjElectricity at nodes i and j, respectivelyPressing;
from the equations (25) to (26), the injection power equation of the node is a linear equation with respect to the node voltage V and the phase angle δ, and a set of linear equations can be obtained by using the equations (25) to (26) for the nodes other than the balanced node, and the voltage and the phase angle of each node can be obtained by solving the set of linear equations.
B: the fast calculation model of the heat supply network load flow is as follows:
the pipeline temperature drop equation (10) is Taylor expanded, and a 2-order term is reserved, so that:
Figure BDA0003105535280000131
taking the two-node pipeline shown in fig. 1 as an example, the following equations (27) are applied to node 1 and node 2 in fig. 1, respectively, and are simplified:
Figure BDA0003105535280000132
applying the mathematical induction method to generalize equation (28) to multiple pipelines is:
Figure BDA0003105535280000133
but the heat load exists in the actual heat supply network pipeline, the pipeline flow m1≠m2≠…≠miAs shown in FIG. 2, therefore, equation (29) is an approximate equation for calculating the flow rate miLet us order
Figure BDA0003105535280000134
Then there are:
Figure BDA0003105535280000135
calculating n1,n2…ni-1When the value is the same, neglecting the heat loss of the pipe network, setting the heat quantity flowing through the pipeline i as
Figure BDA0003105535280000136
When calculating the pipe flow m flowing into the load nodeiIn time, the temperature change of the head end and the tail end of the pipeline is small, so that the pipeline can be made
Figure BDA0003105535280000137
Figure BDA0003105535280000138
Wherein k represents the flow miThe number of the pipeline flowing through, when k is i, niAs for formula (30), 1 may be:
Figure BDA0003105535280000139
wherein, TsSupply water temperature to the node, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiIn order to obtain the heat transfer coefficient of the pipe i,
Figure BDA00031055352800001310
Figure BDA00031055352800001311
n i1, k is the flow miThe number of the pipe through which the fluid flows,
Figure BDA00031055352800001312
for heat energy flowing through the conduit k, CpIs the specific heat capacity of water and phi is the thermal load.
Equation (30) can be abbreviated as: ami 2+Bmi+C＝0 (31)
Since the value of C is generally small, equation (31) can be abbreviated as: ami+B＝0 (32)
By using the formula (32), only one linear equation of unity needs to be solved, and meanwhile, the formula (32) decouples the temperature and the flow, decouples the water supply network and the water return network, so that iteration is not needed, and the numerical stability and the calculation efficiency are improved.
And 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides, and the method specifically comprises the following steps:
equation of pressure drop in the pipeline is hf＝Km2To loop pipeline
Figure BDA0003105535280000141
Wherein N is the number of loop pipelines, introducing variable x, and using equation
Figure BDA0003105535280000142
Equivalent transformation into
Figure BDA0003105535280000143
Mixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, the loop can be separated at the convergence position of the pipeline flow rates, and the annular heat supply network is split into a multi-heat-source radiation type heat supply network;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
the flow direction of the adjacent pipes is shown schematically in FIG. 3, when mi+mj＝mqOr C1+C2When is equal to C, indicates miAnd mjThe flow directions of (1) are opposite;
when m isi-mj＝mqOr C1-C2When is equal to C, indicates miAnd mjThe flow direction of (a) is the same, the flow direction is from top to bottom. When m isj-mi＝mqOr C2-C1When is equal to C, indicates miAnd mjThe flow direction of (a) is the same, the flow direction is from bottom to top.
Step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in the multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied with heat from a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, for example, when the heat load is simultaneously supplied with heat from 2 heat sources, as shown in fig. 4, the heat load is
Figure BDA0003105535280000144
Can be equivalent to a heat source H1And H2And simultaneously, heat is respectively supplied, namely the heat exchanger can be split into the diagram 5, and the equivalent method comprises the following steps:
Figure BDA0003105535280000145
Figure BDA0003105535280000146
in the formula, mp、mtAnd mdPipe flow, m, for a multi-heat source networkp1、mp2And mt1、mt2And md1、md2Flow at the corresponding pipeline when the multi-heat-source network is split into 2 single-heat-source networks;
Figure BDA0003105535280000151
for the heat load of the multi-heat source network,
Figure BDA0003105535280000152
and
Figure BDA0003105535280000153
and splitting the multi-heat source network into 2 single heat source networks to correspond to the heat loads at the load nodes. Step 2 can be used to obtain m in FIG. 5d1And md2Because m isd＝md1+md2Then the flow m can be obtainedd。
Step 6: determining a heat supply network power flow linear optimization model according to the decoupling result in the step 5, thereby obtaining an electric-thermal interconnection comprehensive energy system optimization scheduling model, and outputting state quantity information of the electric-thermal interconnection comprehensive energy system, wherein the method specifically comprises the following steps:
as shown in FIG. 1, a heat source node Ts is based on Am in an electric-thermal interconnection integrated energy system fast power flow calculation modeli 2+Bmi+ C ═ 0:
Figure BDA0003105535280000154
omitting the quadratic term to obtain:
Figure BDA0003105535280000155
the method is simplified to obtain:
Figure BDA0003105535280000156
in the above formula, the ambient temperature T is generally setaIs a constant, and the variable is the heat source temperature THAnd heat load return water temperature ToFor practical systems, when determining the temperature T of the heat sourceHAnd heat load return water temperature ToBy least square method, can
Figure BDA0003105535280000157
To be K1TH+K2To+K3Will be
Figure BDA0003105535280000158
To be K4TH+K5To+K6The fitting coefficient can be generally up to 0.99, and the above equation can be:
Figure BDA0003105535280000159
wherein K1，K2，K3，K4，K5，K6Is a constant.
And obtaining a heat supply network power flow linear optimization model in the electric-thermal interconnection comprehensive energy system optimization scheduling model.
And (5) combining the formula (25) and the formula (26) to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation.
The present invention is verified by simulation as follows.
1) EXAMPLES test 1
A test example as shown in fig. 6 was constructed by CHP coupling based on a modified IEEE33 node distribution network and a 23-node heat supply network, where the grid node 1 is a balanced node, the voltage is 1.05p.u., the node 2 is a PV node, the voltage amplitude is 1.049p.u., and the others are PQ nodes. The length of each pipeline of the heat supply network is 1000 meters, the temperature of the CHP source is 100 ℃, the heat load power is 0.5MW, the return water temperature of the load node is 30 ℃, and the external environment temperature T isa10 degrees. M for flow rate is obtained by the formula (32) of the present inventionPIndicates that the flow obtained in the original document is mRThe results of the example test are shown in tables 1 and 2, and the results of the grid example test are shown in fig. 7 and 8.
TABLE 1 different model flow values
Numbering mR(kg/s) mP(kg/s) Numbering mR(kg/s) mP(kg/s)
1 22.1165 22.1672 12 1.8029 1.8056
2 20.3372 20.3864 13 3.6987 3.7076
3 16.6761 16.7184 14 1.8539 1.8584
4 14.8732 14.9128 15 1.8448 1.8492
5 11.1744 11.2052 16 1.8191 1.8229
6 3.7447 3.7556 17 5.6106 5.6267
7 1.8906 1.8964 18 3.7729 3.7845
8 1.7794 1.7808 19 1.8793 1.8850
9 3.6611 3.6680 20 1.8377 1.8423
10 1.8240 1.8273 21 1.8936 1.8995
11 1.8371 1.8408 22 1.8541 1.8592
TABLE 2 Water supply temperature and Return Water temperature for different models
Figure BDA0003105535280000161
Figure BDA0003105535280000171
As can be seen from tables 1 and 2, when there are many system nodes and the heat supply network pipes are long, the flow and temperature errors found in the original document and the present invention become larger as the pipe length increases, but the errors are within the engineering tolerance.
As can be seen from fig. 7 and 8, compared with the nonlinear alternating current model, the power distribution network linearization model does not need iteration and has smaller error.
2) EXAMPLES test 2
Selecting a Bali island comprehensive energy testing system, as shown in FIG. 9, wherein a 9-node power grid is adopted, the total active load is 1.6MW, a 32-node heat supply network is adopted, and the total active power is 2.164 MW; the power grid and the heat supply network are coupled through 3 CHP units, the CHP1 is a gas turbine with a constant heat-electricity ratio, the CHP2 is a steam extraction type turbine with a variable heat-electricity ratio, and the CHP3 is a reciprocating internal combustion engine with a constant heat-electricity ratio. The CHP water supply temperature is constant at 70 ℃, and the heat load return water temperature is constant at 30 ℃. And selecting a power grid node 9 as a power grid balance node, selecting nodes 7 and 8 as PV nodes, and selecting the other nodes as PQ nodes, wherein the heat supply network node 1 is a balance node of a heat supply network.
The flow, temperature, voltage and phase angle obtained by the method of the invention are m respectivelyP、TP、UPAnd thetaPThe flow, temperature, voltage and phase angle obtained in the original document are respectively expressed by mR、TR、URAnd thetaRAnd (4) showing. The percentage of flow error, the percentage of temperature error and the percentage of voltage error are respectively defined as deltam＝(|mP–mR|/mR)×100％、δT＝(|TP-TR|/TR)×100％、δU＝(|UP-UR|/UR) X 100%, and the phase angle error is defined as deltaθ＝|θP–θRThe results of the comparison are shown in tables 3 to 5. As can be seen from the table, the percentage of flow error is 0.5226% at the maximum, and the average of the percentage of flow error is 0.0770%; the maximum percentage of temperature error is 0.0028 percent, and the average value of the percentage of temperature error is 0.0010 percent; the maximum percentage of voltage error is 0.0057%, and the average value of the percentage of voltage error is 0.0042%; the maximum phase angle error is 0.0031 degrees, the average phase angle error value is 0.0017 degrees, and the method provided by the invention is high in precision.
TABLE 3 Heat flow network test results
Figure BDA0003105535280000172
Figure BDA0003105535280000181
Table 4 heat supply network temperature test results
Figure BDA0003105535280000182
TABLE 5 grid test results
Node numbering UP(p.u.) θP(°) UR(p.u.) θR(°) δU(％) δθ(°)
1 1.04882 -0.63068 1.04876 -0.62921 0.0057 0.0015
2 1.04888 -0.62885 1.04883 -0.62742 0.0048 0.0014
3 1.04903 -0.66235 1.04897 -0.66190 0.0057 0.0005
4 1.04937 -0.70360 1.04931 -0.70457 0.0057 0.0010
5 1.04999 -0.74550 1.04994 -0.74839 0.0048 0.0029
6 1.04998 -0.73656 1.04994 -0.73931 0.0038 0.0027
7 1.05003 -0.72020 1.05000 -0.72276 0.0029 0.0026
8 1.05005 -0.75527 1.05000 -0.75837 0.0048 0.0031
9 1.02 0 1.02 0 0 0
In summary, the invention provides an optimal scheduling method of an electric-thermal interconnection comprehensive energy system based on a linear equation, firstly, a rapid power flow calculation method of a radiation type heat supply network is deduced, then, a method for decoupling an annular heat supply network into the radiation type heat supply network is proposed, and the solved ring power flow is converted into the solution of a plurality of radiation type heat supply network power flows; and finally, the radiant heat network formula is popularized and deformed, and a least square method is utilized to provide an optimal scheduling method of the electric-heat interconnection comprehensive energy system based on the linear equation, so that the calculation efficiency is improved while the calculation precision is ensured.

Claims (10)

1. An electric-heat interconnection comprehensive energy system optimization scheduling method based on a linear equation is characterized by comprising the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
2. The optimal scheduling method for the electric-thermal interconnection comprehensive energy system based on the linear equation as claimed in claim 1, wherein the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system in the step 3 is specifically:
Figure FDA0003105535270000011
Figure FDA0003105535270000012
Am＝mq (3)
Bhf＝0 (4)
hf＝Km|m| (5)
Figure FDA0003105535270000013
Figure FDA0003105535270000014
Figure FDA0003105535270000015
Φ＝Cpmq(Ts-To) (9)
Figure FDA0003105535270000021
(∑mout)Tout＝∑(minTin) (11)
Cm＝ΦCHP/PCHP (12)
Cz＝ΔΦ/ΔP＝ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiInjected active power and injected reactive power, U, respectively, for node iiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij＝θi-θj，θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
3. The linear equation-based optimal scheduling method for the electric-thermal interconnection energy system according to claim 1, wherein the electric-thermal interconnection energy system fast load flow calculation model in the step 4 comprises a power grid load flow calculation linear model and a thermal power grid load flow fast calculation model.
4. The linear equation-based optimal scheduling method for the electricity-heat interconnection comprehensive energy system, according to claim 3, wherein the power grid load flow calculation linear model is specifically:
Figure FDA0003105535270000022
Figure FDA0003105535270000023
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1＝rijxij/(rij 2+xij 2)，kij_2＝xij 2/(rij 2+xij 2)，rijAnd xijRespectively the resistance and reactance, delta, of the line ijiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjThe voltages at nodes i and j, respectively.
5. The linear equation-based optimal scheduling method for the electric-thermal interconnection comprehensive energy system, according to claim 3, wherein the thermal network load flow fast calculation model is specifically:
Ami 2+Bmi+C＝0 (31)
wherein:
Figure FDA0003105535270000031
wherein, THTemperature of water supplied to heat source, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiFor heat transfer coefficient of conduit i, CpSpecific heat capacity of water, phi thermal load, n1＝m1/mi，n2＝m2/mi，ni-1＝mi-1/mi，ni＝1，。
Figure FDA0003105535270000032
Figure FDA0003105535270000033
ni1, k is the flow miThe number of the pipe through which the fluid flows,
Figure FDA0003105535270000034
is the heat energy flowing through the pipe k.
6. The linear equation-based optimal scheduling method for the electric-thermal interconnection comprehensive energy system according to claim 2, wherein the decoupling of the annular heat supply network in the step 5 is specifically as follows:
step 5-1: splitting the annular heat supply network into a multi-heat-source radiation type heat supply network;
equation of pressure drop in the pipeline is hf＝Km2To loop pipeline
Figure FDA0003105535270000035
Wherein N is the number of loop pipelines; introducing variablesx, will equation
Figure FDA0003105535270000036
Equivalent transformation into
Figure FDA0003105535270000037
Mixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, and the loop can be disconnected at the convergence position of the pipeline flow rates;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in a multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied by a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, and the equivalent method comprises the following steps:
Figure FDA0003105535270000041
Figure FDA0003105535270000042
in the formula, M represents the number of heat sources in the multi-heat-source radiation type heat supply network and is respectively H1,...,HM，
Figure FDA0003105535270000043
Respectively represent a heat source H1,...,HMN represents the number of pipes to be branched from the node, and is respectively denoted by k1,…,kN，
Figure FDA0003105535270000044
Respectively the flow rates of the corresponding pipelines are,
Figure FDA0003105535270000045
respectively the flow of the pipeline
Figure FDA0003105535270000046
Corresponding pipeline flow in a split 1 st, … th and M single heat source radiation type heat network model;
Figure FDA0003105535270000047
for the first heat load in a multi-source radiant network,
Figure FDA0003105535270000048
are respectively as
Figure FDA0003105535270000049
At split 1 st …, M single heat source radiative heat network models correspond to the heat load at the load nodes.
7. The optimal scheduling method for the electricity-heat interconnection comprehensive energy system based on the linear equation as claimed in claim 5, wherein the step 6 of constructing the optimal scheduling model for the electricity-heat interconnection comprehensive energy system comprises the following specific steps:
am in rapid load flow calculation model based on electricity-heat interconnection integrated energy systemi 2+Bmi+ C ═ 0:
Figure FDA00031055352700000410
omitting a quadratic term, and obtaining a heat supply network power flow linear optimization model in the electric-heat interconnection comprehensive energy system optimization scheduling model by using a least square method after simplification:
Figure FDA00031055352700000411
wherein K1，K2，K3，K4，K5，K6Is a constant.
8. An electricity-heat interconnection comprehensive energy system optimization scheduling system based on a linear equation is characterized by comprising the following modules:
the electric-thermal interconnection comprehensive energy system nonlinear steady-state model building module comprises: the method is used for constructing a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system based on parameter information of the electric power and thermodynamic system;
the electric-thermal interconnection comprehensive energy system rapid load flow calculation model construction module comprises: the method is used for constructing a rapid load flow calculation model of the electric-thermal interconnection comprehensive energy system;
a heat supply network decoupling module: for decoupling the looped heat network into a plurality of radiant heat networks;
the electric-thermal interconnection comprehensive energy system optimization scheduling model building module comprises: and constructing an optimized dispatching model of the electric-thermal interconnection comprehensive energy system based on the decoupling result and the quick load flow calculation model of the electric-thermal interconnection comprehensive energy system.
9. A computer arrangement comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1-7 are implemented by the processor when executing the computer program.
10. A computer-storable medium having a computer program stored thereon, wherein the computer program is adapted to carry out the steps of the method according to any one of claims 1-7 when executed by a processor.
CN202110636908.1A 2021-06-08 2021-06-08 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation Pending CN113515853A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110636908.1A CN113515853A (en) 2021-06-08 2021-06-08 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110636908.1A CN113515853A (en) 2021-06-08 2021-06-08 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation

Publications (1)

Publication Number Publication Date
CN113515853A true CN113515853A (en) 2021-10-19

Family

ID=78065567

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110636908.1A Pending CN113515853A (en) 2021-06-08 2021-06-08 Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation

Country Status (1)

Country Link
CN (1) CN113515853A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117332725A (en) * 2023-12-01 2024-01-02 中国科学院广州能源研究所 Steam network dynamic calculation method and device of comprehensive energy system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241662A (en) * 2018-09-29 2019-01-18 河海大学 A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies
CN109412164A (en) * 2018-12-29 2019-03-01 浙江大学 A kind of electric-thermal association system trend processing method
CN109726483A (en) * 2018-12-29 2019-05-07 国网江苏省电力有限公司南京供电分公司 A kind of radial heat supply network model of electric heating interconnection integrated energy system and its system
CN110728032A (en) * 2019-09-20 2020-01-24 河海大学 Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
WO2020155515A1 (en) * 2019-01-30 2020-08-06 广东电网有限责任公司电力调度控制中心 Blockchain-based dual-source energy internet transaction method and device
CN112615367A (en) * 2020-12-09 2021-04-06 国网湖北省电力有限公司电力科学研究院 Optimized scheduling method for comprehensive energy system in power Internet of things environment

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241662A (en) * 2018-09-29 2019-01-18 河海大学 A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies
WO2020093296A1 (en) * 2018-11-06 2020-05-14 国电南瑞科技股份有限公司 Interval power flow calculation method for power-heat integrated energy system
CN109412164A (en) * 2018-12-29 2019-03-01 浙江大学 A kind of electric-thermal association system trend processing method
CN109726483A (en) * 2018-12-29 2019-05-07 国网江苏省电力有限公司南京供电分公司 A kind of radial heat supply network model of electric heating interconnection integrated energy system and its system
WO2020155515A1 (en) * 2019-01-30 2020-08-06 广东电网有限责任公司电力调度控制中心 Blockchain-based dual-source energy internet transaction method and device
CN110728032A (en) * 2019-09-20 2020-01-24 河海大学 Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network
CN112615367A (en) * 2020-12-09 2021-04-06 国网湖北省电力有限公司电力科学研究院 Optimized scheduling method for comprehensive energy system in power Internet of things environment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张勇;马洲俊;张明;许洪华;路晓敏;: "基于线性规划的电-热互联系统参数优化方法", 电力需求侧管理, no. 02 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117332725A (en) * 2023-12-01 2024-01-02 中国科学院广州能源研究所 Steam network dynamic calculation method and device of comprehensive energy system

Similar Documents

Publication Publication Date Title
Zhang et al. Partitional decoupling method for fast calculation of energy flow in a large-scale heat and electricity integrated energy system
WO2020093296A1 (en) Interval power flow calculation method for power-heat integrated energy system
CN110728032B (en) Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network
CN109726483B (en) Radial heat supply network model of electric-heat interconnection comprehensive energy system and system thereof
CN111428351B (en) Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method
Chen et al. Fast decoupled multi-energy flow calculation for integrated energy system
Zhang et al. A fixed-point based distributed method for energy flow calculation in multi-energy systems
Sheng et al. Two-stage state estimation approach for combined heat and electric networks considering the dynamic property of pipelines
Li et al. Classification and location scheme selection of coupling components in integrated electrical and heating systems with renewable energy
CN113255105B (en) Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure
Wang et al. Decoupled optimization of integrated energy system considering CHP plant based on energy hub model
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN113515853A (en) Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation
Sun et al. Multi-energy flow calculation method for we-energy based energy internet
CN113690891B (en) Analysis-method-based probability power flow determination method for electric-thermal interconnection comprehensive energy system
CN113221428B (en) Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system
CN115081193A (en) Power grid-heat grid model construction method, fusion simulation method and system
CN114549232A (en) Hybrid energy flow calculation method for electricity-heat comprehensive energy system
Yao et al. A novel cross iteration method for dynamic energy flow calculation of the hot-water heating network in the integrated energy system
Lei et al. A Novel Methodology for Electric-Thermal Mixed Power Flow Simulation and Transmission Loss Analysis in Multi-Energy Micro-Grids
Zheng et al. Integrated planning of coupling components in rural combined heat and power systems considering biomass energy
CN115906411B (en) Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics
Shen et al. Optimal dispatch model for district heating network based on interior-point method
Markensteijn et al. A graph-based framework for steady-state load flow analysis of multi-carrier energy networks
Zhang et al. Construction and Optimization of a Thermo-economic Cost Analysis Model for the Complex Energy Network.

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