CN113255105A - Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure - Google Patents

Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure Download PDF

Info

Publication number
CN113255105A
CN113255105A CN202110454003.2A CN202110454003A CN113255105A CN 113255105 A CN113255105 A CN 113255105A CN 202110454003 A CN202110454003 A CN 202110454003A CN 113255105 A CN113255105 A CN 113255105A
Authority
CN
China
Prior art keywords
network
power
heat
node
model
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
CN202110454003.2A
Other languages
Chinese (zh)
Other versions
CN113255105B (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.)
Shanghai Electric Power University
Original Assignee
Shanghai 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 Shanghai Electric Power University filed Critical Shanghai Electric Power University
Priority to CN202110454003.2A priority Critical patent/CN113255105B/en
Publication of CN113255105A publication Critical patent/CN113255105A/en
Application granted granted Critical
Publication of CN113255105B publication Critical patent/CN113255105B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

The invention relates to a load flow calculation method of an electric and thermal comprehensive energy system with a bidirectional coupling network structure, which comprises the following steps: (1) constructing a network topological structure model of the electric-gas-heat bidirectional coupling network; (2) acquiring network structure parameters and coupling equipment parameters in the electricity-gas-heat bidirectional coupling network; (3) determining load data of each network node in the electric-gas-heat bidirectional coupling network; (4) establishing a steady-state power flow model of the electricity-gas-heat bidirectional coupling network; (5) and performing electric-gas-heat tidal current solving by using a Newton-Raphson iteration method based on the steady-state tidal current model. Compared with the prior art, the invention provides the calculation method for quickly and accurately calculating the load flow of the electric and thermal comprehensive energy system with the bidirectional coupling network structure, and the calculation method has high solving speed and good convergence.

Description

Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure
Technical Field
The invention relates to the field of comprehensive energy systems and multi-energy flow load flow calculation, in particular to a load flow calculation method of an electric and thermal comprehensive energy system with a bidirectional coupling network structure.
Background
With the rapid development of technologies such as natural gas power generation, Combined Heat and Power (CHP), Combined Cooling Heating and Power (CCHP) and the like, the coupling among various energy sources such as electricity, gas, heat, cold and the like is continuously deepened, and the proposition of concepts such as an energy internet, a comprehensive energy system and the like opens up a new idea for the utilization of future energy sources. The CHP unit, the CCHP generator unit, the electric boiler and the gas boiler are used as coupling units in a regional comprehensive energy system, and are simultaneously connected with distributed power generation units such as photovoltaic power generation unit and wind power generation unit, so that the tight coupling of various heterogeneous energies such as electricity, gas and heat is realized, the energy utilization efficiency is greatly improved, a distributed power supply is consumed, and the problems of 'light abandonment' and 'wind abandonment' are effectively solved. Therefore, the regional comprehensive energy system draws high attention of the government of China, and immediately a batch of comprehensive energy microgrid exemplary projects are built, so that a theoretical basis is provided for large-scale popularization and application of the comprehensive energy system. Energy internet and comprehensive energy system have attracted wide attention of numerous scholars at home and abroad as emerging research fields. The trend of multi-energy coupling in the comprehensive energy system is gradually deepened, the running state of the electric-gas-heat comprehensive energy system can be determined through multi-energy flow trend research, and the multi-energy flow trend research is the basis of planning, running and controlling of the comprehensive energy system, so that the multi-energy flow trend research of the multi-energy coupling system has important significance.
The multi-energy flow system sub-network mainly comprises a power network, a natural gas network and a heat power network, and an electric heating coupling system, an electric coupling system and an electric-gas-heat coupling system are formed by coupling units. In the multi-energy flow power research, an electric power flow model is firstly proposed due to the development of natural gas power generation technology. Researchers analogize the power system power flow model to establish a natural gas network power flow model and provide a sequential solution of electric and hybrid power flow; then, a heat supply network recursive model is established for the tree-shaped pipe network, a heat supply network power flow forward-backward substitution algorithm is provided, a ring network power flow model is established for the annular heat supply network, and a unified method and a decomposition method are provided for solving the electric-heat combined power flow; research researchers have proposed an energy hub model aiming at electricity-gas-heat multi-energy flow trend, and coupling conversion of multi-energy flow is represented through the energy hub model. The invention aims at a bidirectional coupling electricity-gas-heat network structure, establishes a bidirectional coupling electricity-gas-heat network power flow model, and provides a power flow calculation method of the bidirectional coupling electricity-gas-heat network structure.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a load flow calculation method of an electric and thermal integrated energy system with a bidirectional coupling network structure.
The purpose of the invention can be realized by the following technical scheme:
a load flow calculation method of an electric and thermal integrated energy system with a bidirectional coupling network structure comprises the following steps:
(1) constructing a network topological structure model of the electric-gas-heat bidirectional coupling network;
(2) acquiring network structure parameters and coupling equipment parameters in the electricity-gas-heat bidirectional coupling network;
(3) determining load data of each network node in the electric-gas-heat bidirectional coupling network;
(4) establishing a steady-state power flow model of the electricity-gas-heat bidirectional coupling network;
(5) and performing electric-gas-heat tidal current solving by using a Newton-Raphson iteration method based on the steady-state tidal current model.
Preferably, the network structure parameters obtained in step (2) include the length of a wire and the resistance of a wire of each branch in the power network, the diameter of a pipe, the roughness of a pipe and the length of a pipe of each branch in the thermodynamic network, and the diameter of a pipe, the roughness of a pipe and the length of a pipe of each branch in the thermodynamic network.
Preferably, the coupled plant parameters obtained in step (2) include cogeneration plant parameters, electric boiler plant parameters.
Preferably, the steady-state power flow model established in step (4) includes a power network model, a natural gas network model, a thermal network model and a coupling equipment model.
Preferably, the power network model comprises:
voltage equation in polar form:
Figure BDA0003039901900000021
wherein, UiIs the voltage vector, | U, of power node iiI is the voltage amplitude of the power node i, z is the imaginary unit, deltaiThe phase angle of the power node i is shown, and e is a natural constant;
power equation in polar form:
Figure BDA0003039901900000022
Figure BDA0003039901900000023
wherein, Δ Pi、ΔQiRespectively an active power error vector, a reactive power error vector, Pi、QiFor active and reactive power injected into power node i, | Ui|、|UjI is the voltage amplitude of the power node i, j, Gij、BijIs the real and imaginary part, delta, of the admittance matrix between power nodes i and jijAnd n is the total number of power nodes in the power network.
Preferably, the natural gas network model is specifically:
Figure BDA0003039901900000031
BgΔpg＝0
Δpg＝Kgvg|vg|k-1
wherein A isgIs a node-branch incidence matrix of the natural gas network, vgVolume flow of pipeline for natural gas network, vq gFor node load flow of natural gas networks, BgIs a loop-branch correlation matrix, Δ p, of a natural gas networkgFor the pressure difference, K, across the pipeline joint in the natural gas networkgK is a constant coefficient;
Agmiddle element Ag ijThe value of (A) is as follows:
Figure BDA0003039901900000032
Ag ijis represented by AgThe ith row and the jth column in the natural gas network are divided into i nodes except balance nodes and j nodes except balance nodes;
Bgmiddle element Bg 1jThe value of (A) is as follows:
Figure BDA0003039901900000033
Bg 1jthe value of a pipeline j in 1 natural gas network loop and a specified positive direction is represented, and j is taken through all pipelines in the natural gas network loop;
Kgexpressed as:
Figure BDA0003039901900000034
in the formula: l isgFor the length of the natural gas pipeline, DgIs the natural gas pipeline diameter.
Preferably, the thermodynamic network model is specifically:
Figure BDA0003039901900000035
Bhhf＝0
hf＝Khmh|mh|
Figure BDA0003039901900000036
Figure BDA0003039901900000037
∑(mout)Tout＝∑minTin。
wherein A ishIs a node-branch incidence matrix, m, of a thermodynamic networkhIs the pipe mass flow vector, mq hAs heat load node mass flow vector, BhIs a loop-branch correlation matrix of the thermodynamic network, hfFor the vector of the head loss of the pipe caused by friction, KhIs a constant coefficient, phi is a thermal load node power vector, CpIs the specific heat capacity of water, Tstart、TendTemperatures of the beginning and end nodes of the pipeline, respectively, with T in the heating networkstart＝TsWith T in the regenerative networkstart＝Tr，TsSupply temperature, T, to the noderIs the return temperature of the node, ToTaking a fixed value for the outlet temperature of the node, λ being the heat transfer coefficient per unit length of the pipe, TaIs the ambient temperature, mout、ToutAnd min、TinMass flow and temperature of water flowing out and into the mixing node respectively;
Ahmiddle element Ah ijThe value of (A) is as follows:
Figure BDA0003039901900000041
Ah ijis represented by AhThe ith row and the jth column in the heat distribution network are divided into i nodes and j nodes, wherein the i nodes are the number of all nodes except balance nodes in the heat distribution network, and the j nodes are the number of all pipelines in the heat distribution network;
Bhmiddle element Bh ijThe value of (A) is as follows:
Figure BDA0003039901900000042
Bh 1jthe value of a pipeline j in 1 heat power network loop and a specified positive direction is shown, and j is taken through all pipelines in the heat power network loop;
Khexpressed as:
Figure BDA0003039901900000043
wherein L ishFor heat supply network pipe length, DhIs the diameter of the heat supply network pipeline, f is the pipeline friction factor, ρ is the density of water, and g is the acceleration of gravity.
Preferably, the coupling plant model comprises a cogeneration plant model and an electric boiler model.
Preferably, for a first cogeneration plant with an indefinite heat to power ratio, the coupling model for the first cogeneration plant is:
Figure BDA0003039901900000044
Figure BDA0003039901900000045
wherein Z is a ratio of a thermal power variation amount and an electric power variation amount of the first cogeneration apparatus, ΦCHP1Is the thermal power, P, of the first cogeneration unitCHP1Is the electric power of a first cogeneration plant, PconIs P in pure coagulation modeconElectric power of the unit in the pure condensing mode, Fg CHP1Is the gas consumption of the first cogeneration plant, ke0、ke1、ke2、ke3、ke4、ke5Is a polynomial coefficient of gas consumption of the first cogeneration unit;
for a second cogeneration plant with a certain heat-to-power ratio, the coupling model of the second cogeneration plant is:
Figure BDA0003039901900000051
Figure BDA0003039901900000052
wherein, CmIs the heat-to-power ratio, phi, of the second cogeneration unitCHP2Is the thermal power, P, of the second cogeneration unitCHP2Is the electric power of a second cogeneration plant, Fg CHP2For the gas consumption of the second cogeneration plant, kb0、kb1、kb2Is a polynomial coefficient of the gas consumption of the second cogeneration unit.
Preferably, the electric boiler model includes:
Figure BDA0003039901900000053
wherein eta isBEfficiency of electric boilers, phiBIs the thermal power of electric boiler, phiBIs the electrical power of an electrical boiler.
Compared with the prior art, the invention has the following advantages:
(1) the invention provides a load flow calculation method of a bidirectional coupling electricity-gas-heat network structure aiming at establishing a bidirectional coupling electricity-gas-heat network load flow model for an electricity-gas-heat bidirectional coupling network, and realizes accurate calculation of load flow.
(2) The combined decomposition alternate solving method adopted by the invention can effectively solve the problem of load flow calculation of the electric and thermal comprehensive energy system with the bidirectional coupling network structure, and has the advantages of high solving speed, good convergence and the like.
Drawings
Fig. 1 is a flow chart of an electrical and thermal integrated energy system load flow calculation method of a bidirectional coupling network structure according to the invention;
fig. 2 is a schematic network topology diagram of an electrical and thermal integrated energy system with a bidirectional coupling network structure according to an embodiment of the present invention;
fig. 3 is a power flow distribution of the power network according to the embodiment;
fig. 4 is a flow distribution of the natural gas network according to the embodiment;
fig. 5 is a flow distribution of the thermal network according to the embodiment.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. Note that the following description of the embodiments is merely a substantial example, and the present invention is not intended to be limited to the application or the use thereof, and is not limited to the following embodiments.
Examples
As shown in fig. 1, the present embodiment provides a method for calculating a power flow of an electrical and thermal integrated energy system in a bidirectional coupling network structure, where the method includes the following steps:
(1) constructing a network topological structure model of the electric-gas-heat bidirectional coupling network;
(2) acquiring network structure parameters and coupling equipment parameters in the electricity-gas-heat bidirectional coupling network;
(3) determining load data of each network node in the electric-gas-heat bidirectional coupling network;
(4) establishing a steady-state power flow model of the electricity-gas-heat bidirectional coupling network;
(5) and performing electric-gas-heat tidal current solving by using a Newton-Raphson iteration method based on the steady-state tidal current model.
The network topology of the electric-gas-heat bidirectional coupling network in the embodiment is shown in fig. 2, and includes 14 power nodes: e1, e2, … …, e14, 14 thermal nodes: h1, h2, … …, h14, 9 natural gas nodes g1, g2, … … and g9, and comprises two cogeneration units CHP1, CHP2 and an electric boiler B.
Therefore, the network structure parameters obtained in step (2) include the length of the wire and the resistance value of the wire of each branch in the power network, the pipe diameter, the roughness of the pipeline and the length of the pipeline of each branch in the thermodynamic network, and the pipe diameter, the roughness of the pipeline and the length of the pipeline of each branch in the thermodynamic network. The coupling equipment parameters obtained in the step (2) comprise parameters of cogeneration equipment and parameters of electric boiler equipment.
The steady-state power flow model established in the step (4) comprises a power network model, a natural gas network model, a thermal network model and a coupling equipment model.
The power network model includes:
voltage equation in polar form:
Figure BDA0003039901900000061
wherein, UiIs the voltage vector, | U, of power node iiI is the voltage amplitude of the power node i, z is the imaginary unit, deltaiThe phase angle of the power node i is shown, and e is a natural constant;
power equation in polar form:
Figure BDA0003039901900000062
Figure BDA0003039901900000063
wherein, Δ Pi、ΔQiAre error vectors of active power and reactive power, Pi、QiFor active and reactive power injected into power node i, | Ui|、|UjI is the voltage amplitude of the power node i, j, Gij、BijIs the real and imaginary part, delta, of the admittance matrix between power nodes i and jijAnd n is the total number of power nodes in the power network.
The natural gas network model specifically comprises the following steps:
Figure BDA0003039901900000071
BgΔpg＝0
Δpg＝Kgvg|vg|k-1
wherein A isgIs a node-branch incidence matrix of the natural gas network, vgVolume flow of pipeline for natural gas network, vq gFor node load flow of natural gas networks, BgIs a loop-branch correlation matrix, Δ p, of a natural gas networkgFor the pressure difference, K, across the pipeline joint in the natural gas networkgK is a constant coefficient;
Agmiddle element Ag ijThe value of (A) is as follows:
Figure BDA0003039901900000072
Ag ijis represented by AgThe ith row and the jth column in the natural gas network are divided into i nodes except balance nodes and j nodes except balance nodes;
Bgmiddle element Bg 1jThe value of (A) is as follows:
Figure BDA0003039901900000073
Bg 1jthe value of a pipeline j in 1 natural gas network loop and a specified positive direction is represented, and j is taken through all pipelines in the natural gas network loop;
Kgexpressed as:
Figure BDA0003039901900000074
in the formula: l isgFor the length of the natural gas pipeline, DgIs the natural gas pipeline diameter.
The thermodynamic network model specifically comprises:
Figure BDA0003039901900000075
Bhhf＝0
hf＝Khmh|mh|
Figure BDA0003039901900000076
Figure BDA0003039901900000077
∑(mout)Tout＝∑minTin。
wherein A ishIs a node-branch incidence matrix, m, of a thermodynamic networkhIs the pipe mass flow vector, mq hAs heat load node mass flow vector, BhIs a loop-branch correlation matrix of the thermodynamic network, hfFor the vector of the head loss of the pipe caused by friction, KhIs a constant coefficient, phi is a thermal load node power vector, CpIs the specific heat capacity of water, Tstart、TendTemperatures of the beginning and end nodes of the pipeline, respectively, with T in the heating networkstart＝TsWith T in the regenerative networkstart＝Tr，TsSupply temperature, T, to the noderIs the return temperature of the node, ToTaking a fixed value for the outlet temperature of the node, λ being the heat transfer coefficient per unit length of the pipe, TaIs the ambient temperature, mout、ToutAnd min、TinMass flow and temperature of water flowing out and into the mixing node respectively;
Ahmiddle element Ah ijThe value of (A) is as follows:
Figure BDA0003039901900000081
Ah ijis represented by AhRow i and j ofThe elements of the column, i are taken as the number of all nodes except the balance node in the thermal network, and j is taken as the number of all pipelines in the thermal network;
Bhmiddle element Bh ijThe value of (A) is as follows:
Figure BDA0003039901900000082
Bh 1jthe value of a pipeline j in 1 heat power network loop and a specified positive direction is shown, and j is taken through all pipelines in the heat power network loop;
Khexpressed as:
Figure BDA0003039901900000083
wherein L ishFor heat supply network pipe length, DhIs the diameter of the heat supply network pipeline, f is the pipeline friction factor, ρ is the density of water, and g is the acceleration of gravity.
The coupled plant model includes a cogeneration plant model and an electric boiler model.
For a first cogeneration plant with an indefinite heat to power ratio, the coupling model for the first cogeneration plant is:
Figure BDA0003039901900000084
Figure BDA0003039901900000085
wherein Z is a ratio of a thermal power variation amount and an electric power variation amount of the first cogeneration apparatus, ΦCHP1Is the thermal power, P, of the first cogeneration unitCHP1Is the electric power of a first cogeneration plant, PconIs P in pure coagulation modeconElectric power of the unit in the pure condensing mode, Fg CHP1Is the gas consumption of the first cogeneration plant, ke0、ke1、ke2、ke3、ke4、ke5Is a polynomial coefficient of gas consumption of the first cogeneration unit;
for a second cogeneration plant with a certain heat-to-power ratio, the coupling model of the second cogeneration plant is:
Figure BDA0003039901900000091
Figure BDA0003039901900000092
wherein, CmIs the heat-to-power ratio, phi, of the second cogeneration unitCHP2Is the thermal power, P, of the second cogeneration unitCHP2Is the electric power of a second cogeneration plant, Fg CHP2For the gas consumption of the second cogeneration plant, kb0、kb1、kb2Is a polynomial coefficient of the gas consumption of the second cogeneration unit.
The electric boiler model includes:
Figure BDA0003039901900000093
wherein eta isBEfficiency of electric boilers, phiBIs the thermal power of electric boiler, phiBIs the electrical power of an electrical boiler.
Solving the electricity-gas-heat current, comprising the following steps:
the system correction equation is:
X(i+1)＝X(i)-(J(i))-1F(i)
and can be represented as:
F(X(i))＝-J(i)ΔX(i)
in the formula: x is an independent variable vector corresponding to delta and U in the power grid and m in the heat supply networkh、Ts、TrCorresponding to v in the air networkgP; j is a Jacobian matrix, and the specific correspondence is provided for electricity, heat and gas in the following text; f is the error vector which has specific correspondence to electricity, heat and gas in the following.
Aiming at an electric-gas-heat bidirectional coupling network structure, firstly solving the trend of an electric-heat coupling system:
Figure BDA0003039901900000094
in the formula: (i) in the superscript represents the ith iteration in the iterative solution, the form being the same below; delta delta(i)Is delta(i+1)And delta(i)A difference of (d); delta U(i)Is U(i+1)And U(i)The difference of (a).
Figure BDA0003039901900000095
In the formula: Δ m(i)Is m(i+1)And m(i)A difference of (d); delta Ts ’(i)Is Ts ’(i+1)And Ts ’(i)A difference of (d); delta Tr ’(i)Is Tr ’(i+1)And Tr ’(i)The difference of (a).
For CHP1, phi corresponds to the ith iteration(i) CHP1And P(i) CHP1(ii) a For CHP2, phi corresponds to the ith iteration(i) CHP2And P(i) CHP2。
Phi when the power flow of the heat supply network iterates for the first time(1) CHP2The node is used as a heat supply network balance node and is firstly assigned with phi(1) CHP1(ii) a Secondly, obtaining the required quantity m of the thermodynamic system by the power flow equation of the heat supply network in the step (5)(1)、Ts(1)、Tr(1)(ii) a Calculating to obtain the thermal power phi of the balance node(1) CHP2A value of (d);
③ from phi(1) CHP2The value of (A) and the coupling model of CHP2 can result in P(1) CHP2；
Electric networkLoad flow calculation time P(1) CHP1The node is used as a power grid balancing node and is composed of P(1) CHP2The value of (5) and the power grid load flow equation in the step (5) are used for solving the delta of the power system(1)、U(1)(ii) a When the first iteration is completed, a balance node P can be obtained(1) CHP1A value of (d);
is composed of P(1) CHP1The coupling model of CHP1 and the value of phi can be updated(1)’ CHP1Calculating phi(1)’ CHP1And phi(1) CHP1If the difference value is within the error allowable range, stopping calculating; otherwise, it is updated to phi(2) CHP1And returning to the step II to continue the next iteration.
And sixthly, power flow data of the electric power system and the thermodynamic system are obtained through bidirectional coupling balance of the CHP1, the CHP2 and the electric boiler.
And (3) then obtaining gas consumption load data of CHP1 and CHP2 by the step (2) and a coupling model of CHP1 and CHP2, and further obtaining a natural gas system power flow:
Figure BDA0003039901900000101
fig. 3 shows the flow distribution of the power network, fig. 4 shows the flow distribution of the natural gas network, and fig. 5 shows the flow distribution of the thermal network.
The above embodiments are merely examples and do not limit the scope of the present invention. These embodiments may be implemented in other various manners, and various omissions, substitutions, and changes may be made without departing from the technical spirit of the present invention.

Claims (10)

1. A load flow calculation method of an electric and thermal integrated energy system with a bidirectional coupling network structure is characterized by comprising the following steps:
(1) constructing a network topological structure model of the electric-gas-heat bidirectional coupling network;
(2) acquiring network structure parameters and coupling equipment parameters in the electricity-gas-heat bidirectional coupling network;
(3) determining load data of each network node in the electric-gas-heat bidirectional coupling network;
(4) establishing a steady-state power flow model of the electricity-gas-heat bidirectional coupling network;
(5) and performing electric-gas-heat tidal current solving by using a Newton-Raphson iteration method based on the steady-state tidal current model.
2. The electrical and thermal comprehensive energy system power flow calculation method of the bidirectional coupling network structure according to claim 1, wherein the network structure parameters obtained in the step (2) include a wire length and a wire resistance value of each branch in the power network, a pipe diameter size, a pipe roughness and a pipe length of each branch in the thermodynamic network, and a pipe diameter size, a pipe roughness and a pipe length of each branch in the thermodynamic network.
3. The electrical and thermal integrated energy system power flow calculation method of the bidirectional coupling network structure as recited in claim 1, wherein the coupling equipment parameters obtained in the step (2) include cogeneration equipment parameters and electric boiler equipment parameters.
4. The electrical and thermal integrated energy system power flow calculation method of the bidirectional coupling network structure as recited in claim 1, wherein the steady-state power flow model established in the step (4) includes a power network model, a natural gas network model, a thermal network model and a coupling equipment model.
5. The electrical and thermal integrated energy system power flow calculation method of the bidirectional coupling network structure according to claim 4, wherein the power network model comprises:
voltage equation in polar form:
Figure FDA0003039901890000011
wherein, UiIs the voltage vector, | U, of power node iiI is the voltage amplitude of the power node i, z is the imaginary unit, deltaiThe phase angle of the power node i is shown, and e is a natural constant;
power equation in polar form:
Figure FDA0003039901890000012
Figure FDA0003039901890000013
wherein, Δ Pi、ΔQiRespectively an active power error vector, a reactive power error vector, Pi、QiFor active and reactive power injected into power node i, | Ui|、|UjI is the voltage amplitude of the power node i, j, Gij、BijIs the real and imaginary part, delta, of the admittance matrix between power nodes i and jijAnd n is the total number of power nodes in the power network.
6. The electrical and thermal integrated energy system power flow calculation method of the bidirectional coupling network structure according to claim 4, wherein the natural gas network model specifically comprises:
Figure FDA0003039901890000021
BgΔpg＝0
Δpg＝Kgvg|vg|k-1
wherein A isgIs a node-branch incidence matrix of the natural gas network, vgVolume flow of pipeline for natural gas network, vq gFor node load flow of natural gas networks, BgLoop-branch correlation matrix for natural gas networks，ΔpgFor the pressure difference, K, across the pipeline joint in the natural gas networkgK is a constant coefficient;
Agmiddle element Ag ijThe value of (A) is as follows:
Figure FDA0003039901890000022
Ag ijis represented by AgThe ith row and the jth column in the natural gas network are divided into i nodes except balance nodes and j nodes except balance nodes;
Bgmiddle element Bg 1jThe value of (A) is as follows:
Figure FDA0003039901890000023
Bg 1jthe value of a pipeline j in 1 natural gas network loop and a specified positive direction is represented, and j is taken through all pipelines in the natural gas network loop;
Kgexpressed as:
Figure FDA0003039901890000024
in the formula: l isgFor the length of the natural gas pipeline, DgIs the natural gas pipeline diameter.
7. The electrical and thermal integrated energy system power flow calculation method of the bidirectional coupling network structure according to claim 4, wherein the thermodynamic network model specifically comprises:
Figure FDA0003039901890000025
Bhhf＝0
hf＝Khmh|mh|
Figure FDA0003039901890000026
Figure FDA0003039901890000031
∑(mout)Tout＝∑minTin。
wherein A ishIs a node-branch incidence matrix, m, of a thermodynamic networkhIs the pipe mass flow vector, mq hAs heat load node mass flow vector, BhIs a loop-branch correlation matrix of the thermodynamic network, hfFor the vector of the head loss of the pipe caused by friction, KhIs a constant coefficient, phi is a thermal load node power vector, CpIs the specific heat capacity of water, Tstart、TendTemperatures of the beginning and end nodes of the pipeline, respectively, with T in the heating networkstart＝TsWith T in the regenerative networkstart＝Tr，TsSupply temperature, T, to the noderIs the return temperature of the node, ToTaking a fixed value for the outlet temperature of the node, λ being the heat transfer coefficient per unit length of the pipe, TaIs the ambient temperature, mout、ToutAnd min、TinMass flow and temperature of water flowing out and into the mixing node respectively;
Ahmiddle element Ah ijThe value of (A) is as follows:
Figure FDA0003039901890000032
Ah ijis represented by AhThe ith row and the jth column in the heat power network, i is taken as the number of all nodes except the balance node in the heat power network, and j is taken as the number of all tubes in the heat power networkThe number of tracks;
Bhmiddle element Bh ijThe value of (A) is as follows:
Figure FDA0003039901890000033
Bh 1jthe value of a pipeline j in 1 heat power network loop and a specified positive direction is shown, and j is taken through all pipelines in the heat power network loop;
Khexpressed as:
Figure FDA0003039901890000034
wherein L ishFor heat supply network pipe length, DhIs the diameter of the heat supply network pipeline, f is the pipeline friction factor, ρ is the density of water, and g is the acceleration of gravity.
8. The electrical heat integrated energy system flow calculation method of the bidirectional coupling network structure as recited in claim 4, wherein the coupling equipment model includes a cogeneration equipment model and an electric boiler model.
9. The electrical thermal integrated energy system power flow calculation method of the bidirectional coupling network structure as recited in claim 8,
for a first cogeneration plant with an indefinite heat to power ratio, the coupling model for the first cogeneration plant is:
Figure FDA0003039901890000041
Figure FDA0003039901890000042
wherein Z is a first cogenerationThe ratio of the variation of the thermal power to the variation of the electrical power of the device, phiCHP1Is the thermal power, P, of the first cogeneration unitCHP1Is the electric power of a first cogeneration plant, PconIs P in pure coagulation modeconElectric power of the unit in the pure condensing mode, Fg CHP1Is the gas consumption of the first cogeneration plant, ke0、ke1、ke2、ke3、ke4、ke5Is a polynomial coefficient of gas consumption of the first cogeneration unit;
for a second cogeneration plant with a certain heat-to-power ratio, the coupling model of the second cogeneration plant is:
Figure FDA0003039901890000043
Figure FDA0003039901890000044
wherein, CmIs the heat-to-power ratio, phi, of the second cogeneration unitCHP2Is the thermal power, P, of the second cogeneration unitCHP2Is the electric power of a second cogeneration plant, Fg CHP2For the gas consumption of the second cogeneration plant, kb0、kb1、kb2Is a polynomial coefficient of the gas consumption of the second cogeneration unit.
10. The electrical and thermal integrated energy system power flow calculation method of a bidirectional coupling network structure according to claim 8, wherein the electric boiler model comprises:
Figure FDA0003039901890000045
wherein eta isBEfficiency of electric boilers, phiBIs the thermal power of electric boiler, phiBIs the electrical power of an electrical boiler.
CN202110454003.2A 2021-04-26 2021-04-26 Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure Active CN113255105B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110454003.2A CN113255105B (en) 2021-04-26 2021-04-26 Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110454003.2A CN113255105B (en) 2021-04-26 2021-04-26 Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure

Publications (2)

Publication Number Publication Date
CN113255105A true CN113255105A (en) 2021-08-13
CN113255105B CN113255105B (en) 2022-10-18

Family

ID=77221715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110454003.2A Active CN113255105B (en) 2021-04-26 2021-04-26 Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure

Country Status (1)

Country Link
CN (1) CN113255105B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113836739A (en) * 2021-09-30 2021-12-24 沈阳工程学院 Power flow calculation method for electric heating combined system considering multiple losses of heat supply pipe network
CN117134418A (en) * 2023-10-27 2023-11-28 国网浙江省电力有限公司宁波供电公司 Load flow calculation method, device and equipment of multi-energy coupling system

Citations (7)

* 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
CN109347107A (en) * 2018-09-29 2019-02-15 河海大学 One kind be incorporated into the power networks electric heating interconnection integrated energy system tidal current computing method
WO2019200662A1 (en) * 2018-04-20 2019-10-24 东北大学 Stability evaluation and static control method for electricity-heat-gas integrated energy system
CN110601185A (en) * 2019-09-17 2019-12-20 武汉大学 Unified power flow model and random matrix-based comprehensive energy system weak point identification method
CN111400937A (en) * 2020-04-29 2020-07-10 上海电力大学 Load flow calculation method of comprehensive energy system
CN112016033A (en) * 2020-08-07 2020-12-01 大连理工大学 Electric-heat-gas comprehensive energy system load flow calculation method based on forward-backward substitution method
CN112421635A (en) * 2020-10-29 2021-02-26 上海电力大学 Method for calculating multi-energy flow load flow inside and outside multi-energy micro-grid cluster network

Patent Citations (7)

* 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
WO2019200662A1 (en) * 2018-04-20 2019-10-24 东北大学 Stability evaluation and static control method for electricity-heat-gas integrated energy system
CN109347107A (en) * 2018-09-29 2019-02-15 河海大学 One kind be incorporated into the power networks electric heating interconnection integrated energy system tidal current computing method
CN110601185A (en) * 2019-09-17 2019-12-20 武汉大学 Unified power flow model and random matrix-based comprehensive energy system weak point identification method
CN111400937A (en) * 2020-04-29 2020-07-10 上海电力大学 Load flow calculation method of comprehensive energy system
CN112016033A (en) * 2020-08-07 2020-12-01 大连理工大学 Electric-heat-gas comprehensive energy system load flow calculation method based on forward-backward substitution method
CN112421635A (en) * 2020-10-29 2021-02-26 上海电力大学 Method for calculating multi-energy flow load flow inside and outside multi-energy micro-grid cluster network

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
吴琼等: ""基于双向耦合的综合能源系统电热联合潮流计算"", 《热力发电》 *
张刚等: ""考虑多种耦合单元的电气热联合系统潮流分布式计算方法"", 《中国电机工程学报》 *
李杭等: ""考虑多种能源耦合运行方式的综合能源潮流计算方法"", 《分布式能源》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113836739A (en) * 2021-09-30 2021-12-24 沈阳工程学院 Power flow calculation method for electric heating combined system considering multiple losses of heat supply pipe network
CN117134418A (en) * 2023-10-27 2023-11-28 国网浙江省电力有限公司宁波供电公司 Load flow calculation method, device and equipment of multi-energy coupling system
CN117134418B (en) * 2023-10-27 2024-04-12 国网浙江省电力有限公司宁波供电公司 Load flow calculation method, device and equipment of multi-energy coupling system

Also Published As

Publication number Publication date
CN113255105B (en) 2022-10-18

Similar Documents

Publication Publication Date Title
CN111061986B (en) Thermoelectric comprehensive energy system tide calculation method with multiple operation modes
CN109242365A (en) A kind of Interval Power Flow calculation method of electric-thermal interconnection integrated energy system
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN111815042B (en) Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model
CN111428351B (en) Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method
CN113255105B (en) Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure
CN109726483B (en) Radial heat supply network model of electric-heat interconnection comprehensive energy system and system thereof
WO2020062633A1 (en) N-1 static safety analysis method for integrated energy system
CN110955954B (en) Method for reducing optimal load of layered decoupling electric heat comprehensive energy system
CN106844895B (en) Decoupling calculation method for combined cooling heating and power micro-grid energy flow
CN109347107A (en) One kind be incorporated into the power networks electric heating interconnection integrated energy system tidal current computing method
CN110728032B (en) Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network
CN112016033B (en) Electric-thermal-gas comprehensive energy system tide calculation method based on forward-push back substitution method
CN110752605B (en) Optimal power flow calculation method for electric-thermal coupling comprehensive energy system
CN111463795B (en) Electric-thermal system multi-energy-flow analysis method based on improved Suhoff cooling formula
Li et al. Dynamic modeling and operations of a heat-power station system based on renewable energy
Man et al. State estimation for integrated energy system containing electricity, heat and gas
CN112906220B (en) Method for estimating state of comprehensive energy microgrid park system
CN114549232A (en) Hybrid energy flow calculation method for electricity-heat comprehensive energy system
CN113515853A (en) Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation
CN114048699A (en) Method for establishing and calculating electric-heating combined electromagnetic transient simulation model of comprehensive energy system
CN113690891A (en) Analytic method-based power-heat interconnection comprehensive energy system probability power flow determination method
CN111275271A (en) Cold-heat-electricity comprehensive energy system power flow optimization method
Liu et al. Energy Flow Calculation Method of Combined Cooling, Heating and Power System with Terminal Cooling Network Model
Zhou et al. A Unified Energy Flow Analysis in Integrated Energy System Considering State Variables of Energy Hub

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