CN113515853A - Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation - Google Patents
Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation Download PDFInfo
- Publication number
- CN113515853A CN113515853A CN202110636908.1A CN202110636908A CN113515853A CN 113515853 A CN113515853 A CN 113515853A CN 202110636908 A CN202110636908 A CN 202110636908A CN 113515853 A CN113515853 A CN 113515853A
- Authority
- CN
- China
- Prior art keywords
- heat
- electric
- energy system
- comprehensive energy
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 48
- 230000005855 radiation Effects 0.000 claims abstract description 45
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 7
- RZGZMLICFFEUIQ-UHFFFAOYSA-N 5-[(1-phenylcyclohexyl)amino]pentanoic acid Chemical compound C=1C=CC=CC=1C1(NCCCCC(=O)O)CCCCC1 RZGZMLICFFEUIQ-UHFFFAOYSA-N 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000002156 mixing Methods 0.000 claims description 6
- 239000008400 supply water Substances 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000007142 ring opening reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000012360 testing method Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 6
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005485 electric heating Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 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 1
- 102100031277 Calcineurin B homologous protein 1 Human genes 0.000 description 1
- 102100031272 Calcineurin B homologous protein 2 Human genes 0.000 description 1
- 102100031077 Calcineurin B homologous protein 3 Human genes 0.000 description 1
- 241001510512 Chlamydia phage 2 Species 0.000 description 1
- 241000839426 Chlamydia virus Chp1 Species 0.000 description 1
- 101000777252 Homo sapiens Calcineurin B homologous protein 1 Proteins 0.000 description 1
- 101000777239 Homo sapiens Calcineurin B homologous protein 2 Proteins 0.000 description 1
- 101000777270 Homo sapiens Calcineurin B homologous protein 3 Proteins 0.000 description 1
- 101000943802 Homo sapiens Cysteine and histidine-rich domain-containing protein 1 Proteins 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000004083 survival effect Effects 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an electric-thermal interconnection comprehensive energy system optimal scheduling method based on a linear equation, which comprises the steps of decoupling a complex annular heat supply network into a plurality of radiation type heat supply networks, converting ring network power flows into a plurality of radiation type heat supply network power flows, then obtaining parameter information of an electric power system and a thermodynamic system, constructing an electric-thermal interconnection comprehensive energy system nonlinear steady-state model based on the parameter information, constructing an electric-thermal interconnection comprehensive energy system rapid power flow calculation model by utilizing the model, obtaining an electric-thermal interconnection comprehensive energy system optimal scheduling model, and completing scheduling optimization of an electric-thermal interconnection comprehensive energy system. The invention provides an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation by using a least square method, the model is simple, the calculated amount is small, the convergence problem does not exist, and the calculation efficiency is improved while the calculation precision is ensured.
Description
Technical Field
The invention relates to the technical field of operation scheduling and control of an integrated energy system, in particular to an optimal scheduling method of an electric heating interconnection integrated energy system based on a linear equation.
Background
Energy is the basis of human survival and development and is a key factor influencing industrial production and human life quality, so that the reduction of environmental pollution caused by using conventional energy while ensuring sustainable energy supply has become a common concern worldwide. The regional comprehensive energy system adopts the technologies of cogeneration, electric heating technology, gas power generation and the like, takes an electric power grid as a coupling object to fuse a heating power grid and a gas power grid, fully utilizes natural complementary characteristics of various energy forms such as electricity, heat, gas and the like, and realizes the cooperative supply and comprehensive cascade utilization of multi-grid flow.
The current regional comprehensive energy system load flow calculation has many defects, for an electric-heat interconnection comprehensive energy system, a Newton method is usually adopted for solving, but the result is sensitive to an iteration initial value and has the problem of numerical stability, meanwhile, the calculation complexity of a system model is sharply increased along with the increase of the system scale and even possibly exceeds the inherent capability of the existing algorithm, and a new algorithm suitable for a complex model of a large-scale comprehensive energy system is urgently needed to be developed.
For example, a radial heat supply network model of an electric-thermal interconnection comprehensive energy system and a system thereof with the patent number of 201811642220.9 establish a radial heat supply network model of the electric-thermal interconnection comprehensive energy system which can be solved quickly, however, the radial heat supply network model cannot calculate the multi-heat-source heat supply network flow or the annular heat supply network flow and cannot be applied to a large-scale comprehensive energy system. Meanwhile, the radial heat supply network model can only carry out load flow calculation and cannot carry out optimized scheduling.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide an optimal scheduling method of an electric-heat interconnection comprehensive energy system based on a linear equation, which improves the calculation efficiency while ensuring the calculation precision.
The technical solution for realizing the purpose of the invention is as follows:
an electric-heat interconnection comprehensive energy system optimization scheduling method based on a linear equation is characterized by comprising the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
Compared with the prior art, the invention has the remarkable advantages that:
(1) the invention provides an electric-thermal interconnection comprehensive energy system optimization scheduling method based on a linear equation, which comprises the steps of firstly deducing a rapid power flow calculation method of a radiation type heat supply network, then providing a method for decoupling an annular heat supply network into the radiation type heat supply network, and converting the solved annular power flow into a plurality of solved radiation type heat supply network power flows;
(2) the invention has simple model, extremely small calculated amount and no convergence problem
(3) The invention discloses an electric-thermal interconnection comprehensive energy system optimization scheduling method based on a linear equation, which is characterized in that a radiation type heat supply network formula is expanded and deformed, and a least square method is utilized, so that the calculation efficiency is improved while the calculation precision is ensured.
The present invention will be further described with reference to the following detailed description and accompanying drawings.
Drawings
Fig. 1 is a schematic structural diagram of an ideal radial heat network model in an embodiment of the invention.
Fig. 2 is a schematic structural diagram of an actual radial heat network model in an embodiment of the present invention.
FIG. 3 is a schematic flow diagram of adjacent conduits in an embodiment of the present invention.
Fig. 4 is a schematic structural diagram of a two-heat-source radiation type heat supply network model in an embodiment of the invention.
Fig. 5 is a schematic structural diagram of two single-heat-source radiation type heat supply network models in the embodiment of the invention.
Fig. 6 is a schematic structural diagram of an electric-thermal interconnection comprehensive energy system in the embodiment of the invention.
Fig. 7 is a schematic diagram of the voltage amplitude of the 33-node distribution network in the embodiment of the present invention.
Fig. 8 is a schematic diagram of a voltage phase angle of a 33-node distribution network in the embodiment of the present invention.
Fig. 9 is a schematic structural diagram of an integrated energy system of a Bali island according to an embodiment of the present invention.
Detailed Description
An electric-heat interconnection comprehensive energy system optimal scheduling method based on a linear equation comprises the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: according to the parameter information obtained in the step 1 and the step 2, a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system is constructed, and the method specifically comprises the following steps:
Am=mq (3)
Bhf=0 (4)
hf=Km|m| (5)
Φ=Cpmq(Ts-To) (9)
(∑mout)Tout=∑(minTin) (11)
Cm=ΦCHP/PCHP (12)
Cz=ΔΦ/ΔP=ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiInjected active power and injected reactive power, U, respectively, for node iiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij=θi-θj,θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
The equations (1) - (2) are a power grid steady-state model, the equations (3) - (8) are a heat supply network hydraulic model, the equation (3) is a node flow balance equation, the equation (4) is a loop pressure equation, the equation (5) is a pressure head loss equation, the joint equations (6) - (8) can obtain a pipeline resistance coefficient K, the equations (9) - (11) are heat supply network thermal models, the equation (9) is a heat load power equation, the equation (10) is a pipeline temperature drop equation, and the equation (11) is a node power conservation equation; the fixed hot spot ratio in the coupling element is described by equation (12), and the variable heat power ratio is described by equation (13).
And 4, step 4: the method comprises the following steps of constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system, wherein the rapid power flow calculation model comprises a power grid power flow calculation linear model and a heat supply network power flow rapid calculation model, and the method specifically comprises the following steps:
a: the power grid load flow calculation linear model specifically comprises the following steps:
order Sij=Pij+jQijThen, there are:
let Pij=Pij_1+Pij_2,Qij=Qij_1+Qij_2Then, there are:
because the branch between two nodes is shorter in the power distribution networkTherefore, the phase angle difference is small, and sin delta is takenij≈δij=δi-δj,cosδ ij1, and take | Vi|≈1,|V j1, and the equations (17) - (20) include:
Thereby obtaining a power grid load flow calculation linear model:
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1=rijxij/(rij 2+xij 2),kij_2=xij 2/(rij 2+xij 2),rijAnd xijResistance and electricity of line ij respectivelyAnti, deltaiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjThe voltages at nodes i and j, respectively;
from the equations (25) to (26), the injection power equation of the node is a linear equation with respect to the node voltage V and the phase angle δ, and a set of linear equations can be obtained by using the equations (25) to (26) for the nodes other than the balanced node, and the voltage and the phase angle of each node can be obtained by solving the set of linear equations.
B: the fast calculation model of the heat supply network load flow is as follows:
the pipeline temperature drop equation (10) is Taylor expanded, and a 2-order term is reserved, so that:
taking a pipeline with two nodes as an example, the following formula (27) is respectively applied, and the simplification is achieved:
applying the mathematical induction method to generalize equation (28) to multiple pipelines is:
but the heat load exists in the actual heat supply network pipeline, the pipeline flow m1≠m2≠…≠miTherefore, equation (29) is an approximate equation for solving the flow rate miLet us orderThen there are:
calculating n1,n2…ni-1When the value is the same, neglecting the heat loss of the pipe network, setting the heat quantity flowing through the pipeline i asWhen calculating the pipe flow m flowing into the load nodeiIn time, because the temperature change of the head end and the tail end of the pipeline is small, the pipeline can be used for cooling Wherein k represents the flow miThe number of the pipeline flowing through, when k is i, niAs for formula (30), 1 may be:
wherein, TsSupply water temperature to the node, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiIn order to obtain the heat transfer coefficient of the pipe i, n i1, k is the flow miThe number of the pipe through which the fluid flows,for heat energy flowing through the conduit k, CpIs the specific heat capacity of water and phi is the thermal load.
Equation (30) can be abbreviated as: ami 2+Bmi+C=0 (31)
Since the value of C is generally small, equation (31) can be abbreviated as: ami+B=0 (32)
By using the formula (32), only one linear equation of unity needs to be solved, and meanwhile, the formula (32) decouples the temperature and the flow, decouples the water supply network and the water return network, so that iteration is not needed, and the numerical stability and the calculation efficiency are improved.
And 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides, and the method specifically comprises the following steps:
step 5-1: splitting the annular heat supply network into a multi-heat-source radiation type heat supply network;
equation of pressure drop in the pipeline is hf=Km2To loop pipelineWherein N is the number of loop pipelines, introducing variable x, and using equationEquivalent transformation intoMixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, and the loop can be disconnected at the convergence position of the pipeline flow rates;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in a multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied by a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, and the equivalent method comprises the following steps:
in the formula, M represents the number of heat sources in the multi-heat-source radiation type heat supply network and is respectively H1,...,HM,Respectively represent a heat source H1,...,HMN represents the number of pipes to be branched from the node, and is respectively denoted by k1,…,kN,Respectively the flow rates of the corresponding pipelines are,respectively the flow of the pipelineCorresponding pipeline flow in a split 1 st, … th and M single heat source radiation type heat network model;
for the first heat load in a multi-source radiant network,are respectively asAt split 1 st …, M single heat source radiative heat network models correspond to the heat load at the load nodes.
Step 6: determining a heat supply network power flow linear optimization model according to the decoupling result in the step 5, thereby obtaining an electric-thermal interconnection comprehensive energy system optimization scheduling model, and outputting state quantity information of the electric-thermal interconnection comprehensive energy system, wherein the method specifically comprises the following steps:
electric-thermal interconnection based rapid comprehensive energy systemAm in load flow calculation modeli 2+Bmi+ C ═ 0:
omitting the quadratic term to obtain:
the method is simplified to obtain:
in the above formula, the variable is generally the heat source temperature THAnd heat load return water temperature ToFor practical systems, when determining the temperature T of the heat sourceHAnd heat load return water temperature ToBy least square method, canTo be K1TH+K2To+K3Will beTo be K4TH+K5To+K6The fitting coefficient can be generally up to 0.99, and the above equation can be:
wherein, K1,K2,K3,K4,K5,K6Is a constant.
And obtaining a heat supply network power flow linear optimization model in the electric-thermal interconnection comprehensive energy system optimization scheduling model.
And (5) combining the formula (25) and the formula (26) to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation.
An electricity-heat interconnection comprehensive energy system optimization scheduling system based on a linear equation is characterized by comprising the following modules:
the electric-thermal interconnection comprehensive energy system nonlinear steady-state model building module comprises: the method is used for constructing a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system based on parameter information of the electric power and thermodynamic system;
the electric-thermal interconnection comprehensive energy system rapid load flow calculation model construction module comprises: the method is used for constructing a rapid load flow calculation model of the electric-thermal interconnection comprehensive energy system;
a heat supply network decoupling module: for decoupling the looped heat network into a plurality of radiant heat networks;
the electric-thermal interconnection comprehensive energy system optimization scheduling model building module comprises: and constructing an optimized dispatching model of the electric-thermal interconnection comprehensive energy system based on the decoupling result and the quick load flow calculation model of the electric-thermal interconnection comprehensive energy system.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the computer program:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
A computer-storable medium having stored thereon a computer program which, when executed by a processor, performs the steps of:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system. Examples
An electric-heat interconnection comprehensive energy system optimal scheduling method based on a linear equation comprises the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: according to the parameter information obtained in the step 1 and the step 2, a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system is constructed, and the method specifically comprises the following steps:
Am=mq (3)
Bhf=0 (4)
hf=Km|m| (5)
Φ=Cpmq(Ts-To) (9)
(∑mout)Tout=∑(minTin) (11)
Cm=ΦCHP/PCHP (12)
Cz=ΔΦ/ΔP=ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiRespectively being nodes iInjecting active power and injecting reactive power, UiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij=θi-θj,θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
The equations (1) - (2) are a power grid steady-state model, the equations (3) - (8) are a heat supply network hydraulic model, the equation (3) is a node flow balance equation, the equation (4) is a loop pressure equation, the equation (5) is a pressure head loss equation, the joint equations (6) - (8) can obtain a pipeline resistance coefficient K, the equations (9) - (11) are heat supply network thermal models, the equation (9) is a heat load power equation, the equation (10) is a pipeline temperature drop equation, and the equation (11) is a node power conservation equation; the fixed hot spot ratio in the coupling element is described by equation (12), and the variable heat power ratio is described by equation (13).
And 4, step 4: the method comprises the following steps of constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system, wherein the rapid power flow calculation model comprises a power grid power flow calculation linear model and a heat supply network power flow rapid calculation model, and the method specifically comprises the following steps:
a: the power grid load flow calculation linear model specifically comprises the following steps:
order Sij=Pij+jQijThen, there are:
let Pij=Pij_1+Pij_2,Qij=Qij_1+Qij_2Then, there are:
in the power distribution network, because the branch between two nodes is short, the phase angle difference is small, and sin delta is takenij≈δij=δi-δj,cosδij≈1,And take | Vi|≈1,|V j1, and the equations (17) - (20) include:
Thereby obtaining a power grid load flow calculation linear model:
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1=rijxij/(rij 2+xij 2),kij_2=xij 2/(rij 2+xij 2),rijAnd xijRespectively the resistance and reactance, delta, of the line ijiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjElectricity at nodes i and j, respectivelyPressing;
from the equations (25) to (26), the injection power equation of the node is a linear equation with respect to the node voltage V and the phase angle δ, and a set of linear equations can be obtained by using the equations (25) to (26) for the nodes other than the balanced node, and the voltage and the phase angle of each node can be obtained by solving the set of linear equations.
B: the fast calculation model of the heat supply network load flow is as follows:
the pipeline temperature drop equation (10) is Taylor expanded, and a 2-order term is reserved, so that:
taking the two-node pipeline shown in fig. 1 as an example, the following equations (27) are applied to node 1 and node 2 in fig. 1, respectively, and are simplified:
applying the mathematical induction method to generalize equation (28) to multiple pipelines is:
but the heat load exists in the actual heat supply network pipeline, the pipeline flow m1≠m2≠…≠miAs shown in FIG. 2, therefore, equation (29) is an approximate equation for calculating the flow rate miLet us orderThen there are:
calculating n1,n2…ni-1When the value is the same, neglecting the heat loss of the pipe network, setting the heat quantity flowing through the pipeline i asWhen calculating the pipe flow m flowing into the load nodeiIn time, the temperature change of the head end and the tail end of the pipeline is small, so that the pipeline can be made Wherein k represents the flow miThe number of the pipeline flowing through, when k is i, niAs for formula (30), 1 may be:
wherein, TsSupply water temperature to the node, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiIn order to obtain the heat transfer coefficient of the pipe i, n i1, k is the flow miThe number of the pipe through which the fluid flows,for heat energy flowing through the conduit k, CpIs the specific heat capacity of water and phi is the thermal load.
Equation (30) can be abbreviated as: ami 2+Bmi+C=0 (31)
Since the value of C is generally small, equation (31) can be abbreviated as: ami+B=0 (32)
By using the formula (32), only one linear equation of unity needs to be solved, and meanwhile, the formula (32) decouples the temperature and the flow, decouples the water supply network and the water return network, so that iteration is not needed, and the numerical stability and the calculation efficiency are improved.
And 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides, and the method specifically comprises the following steps:
equation of pressure drop in the pipeline is hf=Km2To loop pipelineWherein N is the number of loop pipelines, introducing variable x, and using equationEquivalent transformation intoMixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, the loop can be separated at the convergence position of the pipeline flow rates, and the annular heat supply network is split into a multi-heat-source radiation type heat supply network;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
the flow direction of the adjacent pipes is shown schematically in FIG. 3, when mi+mj=mqOr C1+C2When is equal to C, indicates miAnd mjThe flow directions of (1) are opposite;
when m isi-mj=mqOr C1-C2When is equal to C, indicates miAnd mjThe flow direction of (a) is the same, the flow direction is from top to bottom. When m isj-mi=mqOr C2-C1When is equal to C, indicates miAnd mjThe flow direction of (a) is the same, the flow direction is from bottom to top.
Step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in the multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied with heat from a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, for example, when the heat load is simultaneously supplied with heat from 2 heat sources, as shown in fig. 4, the heat load isCan be equivalent to a heat source H1And H2And simultaneously, heat is respectively supplied, namely the heat exchanger can be split into the diagram 5, and the equivalent method comprises the following steps:
in the formula, mp、mtAnd mdPipe flow, m, for a multi-heat source networkp1、mp2And mt1、mt2And md1、md2Flow at the corresponding pipeline when the multi-heat-source network is split into 2 single-heat-source networks;for the heat load of the multi-heat source network,andand splitting the multi-heat source network into 2 single heat source networks to correspond to the heat loads at the load nodes. Step 2 can be used to obtain m in FIG. 5d1And md2Because m isd=md1+md2Then the flow m can be obtainedd。
Step 6: determining a heat supply network power flow linear optimization model according to the decoupling result in the step 5, thereby obtaining an electric-thermal interconnection comprehensive energy system optimization scheduling model, and outputting state quantity information of the electric-thermal interconnection comprehensive energy system, wherein the method specifically comprises the following steps:
as shown in FIG. 1, a heat source node Ts is based on Am in an electric-thermal interconnection integrated energy system fast power flow calculation modeli 2+Bmi+ C ═ 0:
omitting the quadratic term to obtain:
the method is simplified to obtain:
in the above formula, the ambient temperature T is generally setaIs a constant, and the variable is the heat source temperature THAnd heat load return water temperature ToFor practical systems, when determining the temperature T of the heat sourceHAnd heat load return water temperature ToBy least square method, canTo be K1TH+K2To+K3Will beTo be K4TH+K5To+K6The fitting coefficient can be generally up to 0.99, and the above equation can be:
wherein K1,K2,K3,K4,K5,K6Is a constant.
And obtaining a heat supply network power flow linear optimization model in the electric-thermal interconnection comprehensive energy system optimization scheduling model.
And (5) combining the formula (25) and the formula (26) to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model based on a linear equation.
The present invention is verified by simulation as follows.
1) EXAMPLES test 1
A test example as shown in fig. 6 was constructed by CHP coupling based on a modified IEEE33 node distribution network and a 23-node heat supply network, where the grid node 1 is a balanced node, the voltage is 1.05p.u., the node 2 is a PV node, the voltage amplitude is 1.049p.u., and the others are PQ nodes. The length of each pipeline of the heat supply network is 1000 meters, the temperature of the CHP source is 100 ℃, the heat load power is 0.5MW, the return water temperature of the load node is 30 ℃, and the external environment temperature T isa10 degrees. M for flow rate is obtained by the formula (32) of the present inventionPIndicates that the flow obtained in the original document is mRThe results of the example test are shown in tables 1 and 2, and the results of the grid example test are shown in fig. 7 and 8.
TABLE 1 different model flow values
Numbering | mR(kg/s) | mP(kg/s) | Numbering | mR(kg/s) | mP(kg/s) |
1 | 22.1165 | 22.1672 | 12 | 1.8029 | 1.8056 |
2 | 20.3372 | 20.3864 | 13 | 3.6987 | 3.7076 |
3 | 16.6761 | 16.7184 | 14 | 1.8539 | 1.8584 |
4 | 14.8732 | 14.9128 | 15 | 1.8448 | 1.8492 |
5 | 11.1744 | 11.2052 | 16 | 1.8191 | 1.8229 |
6 | 3.7447 | 3.7556 | 17 | 5.6106 | 5.6267 |
7 | 1.8906 | 1.8964 | 18 | 3.7729 | 3.7845 |
8 | 1.7794 | 1.7808 | 19 | 1.8793 | 1.8850 |
9 | 3.6611 | 3.6680 | 20 | 1.8377 | 1.8423 |
10 | 1.8240 | 1.8273 | 21 | 1.8936 | 1.8995 |
11 | 1.8371 | 1.8408 | 22 | 1.8541 | 1.8592 |
TABLE 2 Water supply temperature and Return Water temperature for different models
As can be seen from tables 1 and 2, when there are many system nodes and the heat supply network pipes are long, the flow and temperature errors found in the original document and the present invention become larger as the pipe length increases, but the errors are within the engineering tolerance.
As can be seen from fig. 7 and 8, compared with the nonlinear alternating current model, the power distribution network linearization model does not need iteration and has smaller error.
2) EXAMPLES test 2
Selecting a Bali island comprehensive energy testing system, as shown in FIG. 9, wherein a 9-node power grid is adopted, the total active load is 1.6MW, a 32-node heat supply network is adopted, and the total active power is 2.164 MW; the power grid and the heat supply network are coupled through 3 CHP units, the CHP1 is a gas turbine with a constant heat-electricity ratio, the CHP2 is a steam extraction type turbine with a variable heat-electricity ratio, and the CHP3 is a reciprocating internal combustion engine with a constant heat-electricity ratio. The CHP water supply temperature is constant at 70 ℃, and the heat load return water temperature is constant at 30 ℃. And selecting a power grid node 9 as a power grid balance node, selecting nodes 7 and 8 as PV nodes, and selecting the other nodes as PQ nodes, wherein the heat supply network node 1 is a balance node of a heat supply network.
The flow, temperature, voltage and phase angle obtained by the method of the invention are m respectivelyP、TP、UPAnd thetaPThe flow, temperature, voltage and phase angle obtained in the original document are respectively expressed by mR、TR、URAnd thetaRAnd (4) showing. The percentage of flow error, the percentage of temperature error and the percentage of voltage error are respectively defined as deltam=(|mP–mR|/mR)×100%、δT=(|TP-TR|/TR)×100%、δU=(|UP-UR|/UR) X 100%, and the phase angle error is defined as deltaθ=|θP–θRThe results of the comparison are shown in tables 3 to 5. As can be seen from the table, the percentage of flow error is 0.5226% at the maximum, and the average of the percentage of flow error is 0.0770%; the maximum percentage of temperature error is 0.0028 percent, and the average value of the percentage of temperature error is 0.0010 percent; the maximum percentage of voltage error is 0.0057%, and the average value of the percentage of voltage error is 0.0042%; the maximum phase angle error is 0.0031 degrees, the average phase angle error value is 0.0017 degrees, and the method provided by the invention is high in precision.
TABLE 3 Heat flow network test results
Table 4 heat supply network temperature test results
TABLE 5 grid test results
Node numbering | UP(p.u.) | θP(°) | UR(p.u.) | θR(°) | δU(%) | δθ(°) |
1 | 1.04882 | -0.63068 | 1.04876 | -0.62921 | 0.0057 | 0.0015 |
2 | 1.04888 | -0.62885 | 1.04883 | -0.62742 | 0.0048 | 0.0014 |
3 | 1.04903 | -0.66235 | 1.04897 | -0.66190 | 0.0057 | 0.0005 |
4 | 1.04937 | -0.70360 | 1.04931 | -0.70457 | 0.0057 | 0.0010 |
5 | 1.04999 | -0.74550 | 1.04994 | -0.74839 | 0.0048 | 0.0029 |
6 | 1.04998 | -0.73656 | 1.04994 | -0.73931 | 0.0038 | 0.0027 |
7 | 1.05003 | -0.72020 | 1.05000 | -0.72276 | 0.0029 | 0.0026 |
8 | 1.05005 | -0.75527 | 1.05000 | -0.75837 | 0.0048 | 0.0031 |
9 | 1.02 | 0 | 1.02 | 0 | 0 | 0 |
In summary, the invention provides an optimal scheduling method of an electric-thermal interconnection comprehensive energy system based on a linear equation, firstly, a rapid power flow calculation method of a radiation type heat supply network is deduced, then, a method for decoupling an annular heat supply network into the radiation type heat supply network is proposed, and the solved ring power flow is converted into the solution of a plurality of radiation type heat supply network power flows; and finally, the radiant heat network formula is popularized and deformed, and a least square method is utilized to provide an optimal scheduling method of the electric-heat interconnection comprehensive energy system based on the linear equation, so that the calculation efficiency is improved while the calculation precision is ensured.
Claims (10)
1. An electric-heat interconnection comprehensive energy system optimization scheduling method based on a linear equation is characterized by comprising the following steps:
step 1: acquiring parameter information of an electric power system, wherein the parameter information comprises power grid topology, branch parameter information, generator parameter information and electric load information;
step 2: acquiring parameter information of a thermodynamic system, wherein the parameter information comprises heat supply network topology, pipeline parameter information, heat source parameter information and heat load information;
and step 3: constructing a nonlinear steady-state model of the electricity-heat interconnection comprehensive energy system according to the parameter information obtained in the step 1 and the step 2;
and 4, step 4: constructing a rapid power flow calculation model of the electric-thermal interconnection comprehensive energy system by using the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system;
and 5: the flow direction of an internal pipeline of the annular heat supply network is judged and decoupled into a plurality of radiation type heat supply networks, so that the ring network tide is converted into a plurality of radiation type heat supply network tides;
step 6: and (5) determining a heat supply network power flow linear optimization model according to the decoupling result in the step (5), so as to obtain an electric-thermal interconnection comprehensive energy system optimization scheduling model and output the state quantity information of the electric-thermal interconnection comprehensive energy system.
2. The optimal scheduling method for the electric-thermal interconnection comprehensive energy system based on the linear equation as claimed in claim 1, wherein the nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system in the step 3 is specifically:
Am=mq (3)
Bhf=0 (4)
hf=Km|m| (5)
Φ=Cpmq(Ts-To) (9)
(∑mout)Tout=∑(minTin) (11)
Cm=ΦCHP/PCHP (12)
Cz=ΔΦ/ΔP=ΦCHP/(ηeFin-PCHP) (13)
in the formula, PiAnd QiInjected active power and injected reactive power, U, respectively, for node iiIs the voltage of node i, UjIs the voltage at node j, n is the number of branches connected to node i, θij=θi-θj,θijIs the voltage phase angle difference between node i and node j, θiIs the phase angle of node i, thetajIs the node j phase angle; gij、BijRespectively the conductance and susceptance of a pi-type equivalent circuit, A is a network node-branch pipeline incidence matrix, m is the heat supply network pipeline flow, m is the power supply network pipeline flowqFor node incoming load traffic, B is the loop correlation matrix, hfFor the pressure drop of the pipeline caused by friction loss, K is the drag coefficient of the pipeline, L is the length of the pipeline, f is the friction coefficient of the pipeline, D is the diameter of the pipeline, ρ is the water density, g is the acceleration of gravity, Re is the Reynolds number, μ is the kinematic viscosity of water, ε is the roughness of the pipeline, Φ is the thermal load, C is the coefficient of friction of the pipelinepIs the specific heat capacity of water, mqFor node inflow load traffic, TsSupply water temperature to the node, ToIs the node return water temperature, TendIs the temperature at the end of the pipe, TstartFor the head end temperature of the pipeline, TaIs the ambient temperature, λ is the heat transfer coefficient, moutIs the pipe flow of the outflow node, ToutIs the node mixing temperature, minIs the pipe flow into the node, TinIs the temperature at the end of the feed line, CmDetermining the thermoelectric ratio phi for cogeneration unitsCHPIs the heat output of the CHP unit, PCHPIs the electrical output of the CHP unit, CzFor converting the heat-electricity ratio, eta, of cogeneration unitseFor CHP unit condensing efficiency, FinIs the fuel input rate.
3. The linear equation-based optimal scheduling method for the electric-thermal interconnection energy system according to claim 1, wherein the electric-thermal interconnection energy system fast load flow calculation model in the step 4 comprises a power grid load flow calculation linear model and a thermal power grid load flow fast calculation model.
4. The linear equation-based optimal scheduling method for the electricity-heat interconnection comprehensive energy system, according to claim 3, wherein the power grid load flow calculation linear model is specifically:
wherein P isiAnd QiRespectively, the injected active power and the injected reactive power of the node i, NB is the branch number connected with the node i, kij_1=rijxij/(rij 2+xij 2),kij_2=xij 2/(rij 2+xij 2),rijAnd xijRespectively the resistance and reactance, delta, of the line ijiAnd deltajPhase angles, V, of nodes i and j, respectivelyiAnd VjThe voltages at nodes i and j, respectively.
5. The linear equation-based optimal scheduling method for the electric-thermal interconnection comprehensive energy system, according to claim 3, wherein the thermal network load flow fast calculation model is specifically:
Ami 2+Bmi+C=0 (31)
wherein:
wherein, THTemperature of water supplied to heat source, ToIs the node return water temperature, TaIs the ambient temperature, miIs the heat load flow in the pipeline i, LiFor the length of pipe i, λiFor heat transfer coefficient of conduit i, CpSpecific heat capacity of water, phi thermal load, n1=m1/mi,n2=m2/mi,ni-1=mi-1/mi,ni=1,。
6. The linear equation-based optimal scheduling method for the electric-thermal interconnection comprehensive energy system according to claim 2, wherein the decoupling of the annular heat supply network in the step 5 is specifically as follows:
step 5-1: splitting the annular heat supply network into a multi-heat-source radiation type heat supply network;
equation of pressure drop in the pipeline is hf=Km2To loop pipelineWherein N is the number of loop pipelines; introducing variablesx, will equationEquivalent transformation intoMixing it with Am ═ mqObtaining the flow m of each pipeline in the annular heat network simultaneouslyrA univariate linear function with respect to the variable x;
if the sum of the flow rates of the adjacent pipelines is constant, the flow direction of the adjacent loop pipelines is opposite, and the loop can be disconnected at the convergence position of the pipeline flow rates;
if the sum of the flow rates of the adjacent pipelines is still a unitary linear function of the variable x, indicating that the flow directions of the adjacent loop pipelines are the same, and not needing to be subjected to ring opening at the node;
step 5-2: decomposing the multi-heat-source radiation type heat supply network into a plurality of single-heat-source radiation type heat supply networks;
in a multi-heat-source radiation type heat supply network, if a certain heat load is simultaneously supplied by a plurality of heat sources, the heat load is equivalent to a plurality of equivalent heat sources to simultaneously supply heat, so that the multi-heat-source radiation type heat supply network is decomposed into a plurality of single-heat-source radiation type heat supply networks, and the equivalent method comprises the following steps:
in the formula, M represents the number of heat sources in the multi-heat-source radiation type heat supply network and is respectively H1,...,HM,Respectively represent a heat source H1,...,HMN represents the number of pipes to be branched from the node, and is respectively denoted by k1,…,kN,Respectively the flow rates of the corresponding pipelines are,respectively the flow of the pipelineCorresponding pipeline flow in a split 1 st, … th and M single heat source radiation type heat network model;
7. The optimal scheduling method for the electricity-heat interconnection comprehensive energy system based on the linear equation as claimed in claim 5, wherein the step 6 of constructing the optimal scheduling model for the electricity-heat interconnection comprehensive energy system comprises the following specific steps:
am in rapid load flow calculation model based on electricity-heat interconnection integrated energy systemi 2+Bmi+ C ═ 0:
omitting a quadratic term, and obtaining a heat supply network power flow linear optimization model in the electric-heat interconnection comprehensive energy system optimization scheduling model by using a least square method after simplification:
wherein K1,K2,K3,K4,K5,K6Is a constant.
8. An electricity-heat interconnection comprehensive energy system optimization scheduling system based on a linear equation is characterized by comprising the following modules:
the electric-thermal interconnection comprehensive energy system nonlinear steady-state model building module comprises: the method is used for constructing a nonlinear steady-state model of the electric-thermal interconnection comprehensive energy system based on parameter information of the electric power and thermodynamic system;
the electric-thermal interconnection comprehensive energy system rapid load flow calculation model construction module comprises: the method is used for constructing a rapid load flow calculation model of the electric-thermal interconnection comprehensive energy system;
a heat supply network decoupling module: for decoupling the looped heat network into a plurality of radiant heat networks;
the electric-thermal interconnection comprehensive energy system optimization scheduling model building module comprises: and constructing an optimized dispatching model of the electric-thermal interconnection comprehensive energy system based on the decoupling result and the quick load flow calculation model of the electric-thermal interconnection comprehensive energy system.
9. A computer arrangement comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1-7 are implemented by the processor when executing the computer program.
10. A computer-storable medium having a computer program stored thereon, wherein the computer program is adapted to carry out the steps of the method according to any one of claims 1-7 when executed by a processor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110636908.1A CN113515853B (en) | 2021-06-08 | 2021-06-08 | Optimal scheduling method of electrothermal interconnection comprehensive energy system based on linear equation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110636908.1A CN113515853B (en) | 2021-06-08 | 2021-06-08 | Optimal scheduling method of electrothermal interconnection comprehensive energy system based on linear equation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113515853A true CN113515853A (en) | 2021-10-19 |
CN113515853B CN113515853B (en) | 2024-06-18 |
Family
ID=78065567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110636908.1A Active CN113515853B (en) | 2021-06-08 | 2021-06-08 | Optimal scheduling method of electrothermal interconnection comprehensive energy system based on linear equation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113515853B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117332725A (en) * | 2023-12-01 | 2024-01-02 | 中国科学院广州能源研究所 | Steam network dynamic calculation method and device of comprehensive energy system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241662A (en) * | 2018-09-29 | 2019-01-18 | 河海大学 | A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies |
CN109412164A (en) * | 2018-12-29 | 2019-03-01 | 浙江大学 | A kind of electric-thermal association system trend processing method |
CN109726483A (en) * | 2018-12-29 | 2019-05-07 | 国网江苏省电力有限公司南京供电分公司 | A kind of radial heat supply network model of electric heating interconnection integrated energy system and its system |
CN110728032A (en) * | 2019-09-20 | 2020-01-24 | 河海大学 | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network |
WO2020093296A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Interval power flow calculation method for power-heat integrated energy system |
WO2020155515A1 (en) * | 2019-01-30 | 2020-08-06 | 广东电网有限责任公司电力调度控制中心 | Blockchain-based dual-source energy internet transaction method and device |
CN112615367A (en) * | 2020-12-09 | 2021-04-06 | 国网湖北省电力有限公司电力科学研究院 | Optimized scheduling method for comprehensive energy system in power Internet of things environment |
-
2021
- 2021-06-08 CN CN202110636908.1A patent/CN113515853B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241662A (en) * | 2018-09-29 | 2019-01-18 | 河海大学 | A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies |
WO2020093296A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Interval power flow calculation method for power-heat integrated energy system |
CN109412164A (en) * | 2018-12-29 | 2019-03-01 | 浙江大学 | A kind of electric-thermal association system trend processing method |
CN109726483A (en) * | 2018-12-29 | 2019-05-07 | 国网江苏省电力有限公司南京供电分公司 | A kind of radial heat supply network model of electric heating interconnection integrated energy system and its system |
WO2020155515A1 (en) * | 2019-01-30 | 2020-08-06 | 广东电网有限责任公司电力调度控制中心 | Blockchain-based dual-source energy internet transaction method and device |
CN110728032A (en) * | 2019-09-20 | 2020-01-24 | 河海大学 | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network |
CN112615367A (en) * | 2020-12-09 | 2021-04-06 | 国网湖北省电力有限公司电力科学研究院 | Optimized scheduling method for comprehensive energy system in power Internet of things environment |
Non-Patent Citations (1)
Title |
---|
张勇;马洲俊;张明;许洪华;路晓敏;: "基于线性规划的电-热互联系统参数优化方法", 电力需求侧管理, no. 02 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117332725A (en) * | 2023-12-01 | 2024-01-02 | 中国科学院广州能源研究所 | Steam network dynamic calculation method and device of comprehensive energy system |
Also Published As
Publication number | Publication date |
---|---|
CN113515853B (en) | 2024-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Partitional decoupling method for fast calculation of energy flow in a large-scale heat and electricity integrated energy system | |
WO2020093296A1 (en) | Interval power flow calculation method for power-heat integrated energy system | |
CN110728032B (en) | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network | |
CN109726483B (en) | Radial heat supply network model of electric-heat interconnection comprehensive energy system and system thereof | |
CN111428351B (en) | Electric-thermal comprehensive energy system tide calculation method based on forward-push back substitution method | |
Chen et al. | Fast decoupled multi-energy flow calculation for integrated energy system | |
Zhang et al. | A fixed-point based distributed method for energy flow calculation in multi-energy systems | |
CN113255105B (en) | Load flow calculation method of electric and thermal comprehensive energy system with bidirectional coupling network structure | |
Sheng et al. | Two-stage state estimation approach for combined heat and electric networks considering the dynamic property of pipelines | |
Wang et al. | Decoupled optimization of integrated energy system considering CHP plant based on energy hub model | |
Man et al. | State estimation for integrated energy system containing electricity, heat and gas | |
Li et al. | Classification and location scheme selection of coupling components in integrated electrical and heating systems with renewable energy | |
CN113515853A (en) | Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation | |
Sun et al. | Multi-energy flow calculation method for we-energy based energy internet | |
CN113690891B (en) | Analysis-method-based probability power flow determination method for electric-thermal interconnection comprehensive energy system | |
CN113221428B (en) | Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system | |
CN115081193A (en) | Power grid-heat grid model construction method, fusion simulation method and system | |
CN114549232A (en) | Hybrid energy flow calculation method for electricity-heat comprehensive energy system | |
Lei et al. | A Novel Methodology for Electric-Thermal Mixed Power Flow Simulation and Transmission Loss Analysis in Multi-Energy Micro-Grids | |
Yao et al. | A novel cross iteration method for dynamic energy flow calculation of the hot-water heating network in the integrated energy system | |
CN115906411B (en) | Optimal energy flow modeling method and system for electric heating comprehensive energy system considering full dynamics | |
Zheng et al. | Integrated planning of coupling components in rural combined heat and power systems considering biomass energy | |
Shen et al. | Optimal dispatch model for district heating network based on interior-point method | |
Wang et al. | Operational optimization on large‐scale combined heat and power units with low‐pressure cylinder near‐zero output | |
Markensteijn et al. | A graph-based framework for steady-state load flow analysis of multi-carrier energy networks |
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 |