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 PDFInfo
- 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
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 77
- 238000010168 coupling process Methods 0.000 title claims abstract description 77
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 77
- 230000002457 bidirectional effect Effects 0.000 title claims abstract description 46
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 88
- 239000003345 natural gas Substances 0.000 claims description 44
- 239000013598 vector Substances 0.000 claims description 22
- 239000007789 gas Substances 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 238000010438 heat treatment Methods 0.000 claims description 4
- 108010046685 Rho Factor Proteins 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000015271 coagulation Effects 0.000 claims description 3
- 238000005345 coagulation Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000001172 regenerating effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- JPKJQBJPBRLVTM-OSLIGDBKSA-N (2s)-2-amino-n-[(2s,3r)-3-hydroxy-1-[[(2s)-1-[[(2s)-1-[[(2s)-1-[[(2r)-1-(1h-indol-3-yl)-3-oxopropan-2-yl]amino]-1-oxo-3-phenylpropan-2-yl]amino]-1-oxo-3-phenylpropan-2-yl]amino]-1-oxo-3-phenylpropan-2-yl]amino]-1-oxobutan-2-yl]-6-iminohexanamide Chemical compound C([C@H](NC(=O)[C@@H](NC(=O)[C@@H](N)CCCC=N)[C@H](O)C)C(=O)N[C@@H](CC=1C=CC=CC=1)C(=O)N[C@@H](CC=1C=CC=CC=1)C(=O)N[C@H](CC=1C2=CC=CC=C2NC=1)C=O)C1=CC=CC=C1 JPKJQBJPBRLVTM-OSLIGDBKSA-N 0.000 description 8
- 102100031277 Calcineurin B homologous protein 1 Human genes 0.000 description 8
- 102100031272 Calcineurin B homologous protein 2 Human genes 0.000 description 8
- 241001510512 Chlamydia phage 2 Species 0.000 description 8
- 241000839426 Chlamydia virus Chp1 Species 0.000 description 8
- 101000777252 Homo sapiens Calcineurin B homologous protein 1 Proteins 0.000 description 8
- 101000777239 Homo sapiens Calcineurin B homologous protein 2 Proteins 0.000 description 8
- 101000943802 Homo sapiens Cysteine and histidine-rich domain-containing protein 1 Proteins 0.000 description 8
- 238000011160 research Methods 0.000 description 6
- 238000010248 power generation Methods 0.000 description 5
- 230000005611 electricity Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
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]
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)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
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:
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:
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:
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:
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:
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:
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:
Bhhf＝0
hf＝Khmh|mh|
∑(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:
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:
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:
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:
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:
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:
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:
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:
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:
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:
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:
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:
in the formula: l isgFor the length of the natural gas pipeline, DgIs the natural gas pipeline diameter.
The thermodynamic network model specifically comprises:
Bhhf＝0
hf＝Khmh|mh|
∑(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:
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:
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:
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:
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:
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:
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:
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).
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:
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:
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:
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:
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:
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:
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:
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:
Bhhf＝0
hf＝Khmh|mh|
∑(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:
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:
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:
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:
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:
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:
wherein eta isBEfficiency of electric boilers, phiBIs the thermal power of electric boiler, phiBIs the electrical power of an electrical boiler.
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 (4)
| 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 |
| CN114091244A (en) * | 2021-11-11 | 2022-02-25 | 广东电网有限责任公司广州供电局 | Electric-gas interconnection energy network model construction method and device |
| CN114418220A (en) * | 2022-01-21 | 2022-04-29 | 广西大学 | Construction method and system of comprehensive energy network multi-energy flow data set |
| CN117134418A (en) * | 2023-10-27 | 2023-11-28 | 国网浙江省电力有限公司宁波供电公司 | Load flow calculation method, device and equipment of multi-energy coupling system |
Citations (7)
| 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 |
-
2021
- 2021-04-26 CN CN202110454003.2A patent/CN113255105B/en active Active
Patent Citations (7)
| 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)
| Title |
|---|
| 吴琼等: ""基于双向耦合的综合能源系统电热联合潮流计算"", 《热力发电》 * |
| 张刚等: ""考虑多种耦合单元的电气热联合系统潮流分布式计算方法"", 《中国电机工程学报》 * |
| 李杭等: ""考虑多种能源耦合运行方式的综合能源潮流计算方法"", 《分布式能源》 * |
Cited By (5)
| 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 |
| CN114091244A (en) * | 2021-11-11 | 2022-02-25 | 广东电网有限责任公司广州供电局 | Electric-gas interconnection energy network model construction method and device |
| CN114418220A (en) * | 2022-01-21 | 2022-04-29 | 广西大学 | Construction method and system of comprehensive energy network multi-energy flow data set |
| 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 |
|---|---|---|
| CN113255105B (en) | Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure | |
| CN111061986B (en) | Thermoelectric comprehensive energy system tide calculation method with multiple operation modes | |
| CN111428351B (en) | Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method | |
| 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 | |
| CN112016033B (en) | Electric-thermal-gas comprehensive energy system tide calculation method based on forward-push back substitution method | |
| CN109726483B (en) | Radial heat supply network model of electric-heat interconnection comprehensive energy system and system thereof | |
| CN111815042A (en) | Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model | |
| 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 | |
| CN110752605B (en) | Optimal power flow calculation method for electric-thermal coupling comprehensive energy system | |
| CN110728032B (en) | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network | |
| Li et al. | Dynamic modeling and operations of a heat-power station system based on renewable energy | |
| CN111463795B (en) | Electric-thermal system multi-energy-flow analysis method based on improved Suhoff cooling formula | |
| Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
| CN113515853B (en) | Optimal scheduling method of electrothermal interconnection comprehensive energy system based on linear equation | |
| CN112906220B (en) | Method for estimating state of comprehensive energy microgrid park system | |
| Xiang et al. | Planning Research on Electrically Coupled Integrated Energy System Based on ARIMA-LSTM Model | |
| CN114549232A (en) | Hybrid energy flow calculation method for electricity-heat comprehensive energy system | |
| Hu et al. | Multi-Energy Flow Calculation Model and Solution Method for Regional Integrated Energy System | |
| Liu et al. | Energy Flow Calculation Method of Combined Cooling, Heating and Power System with Terminal Cooling Network Model | |
| CN111275271A (en) | Cold-heat-electricity comprehensive energy system power flow optimization method | |
| Zhou et al. | A Unified Energy Flow Analysis in Integrated Energy System Considering State Variables of Energy Hub | |
| CN115906411B (en) | Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics |
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 |