CN112330493A - Energy system modeling and comprehensive analysis method, device and storage medium - Google Patents

Energy system modeling and comprehensive analysis method, device and storage medium Download PDF

Info

Publication number
CN112330493A
CN112330493A CN202011334339.7A CN202011334339A CN112330493A CN 112330493 A CN112330493 A CN 112330493A CN 202011334339 A CN202011334339 A CN 202011334339A CN 112330493 A CN112330493 A CN 112330493A
Authority
CN
China
Prior art keywords
energy
equation
loss
network
follows
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011334339.7A
Other languages
Chinese (zh)
Other versions
CN112330493B (en
Inventor
谢平平
卢斯煜
禤培正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Southern Power Grid Co Ltd
Research Institute of Southern Power Grid Co Ltd
Original Assignee
China Southern Power Grid Co Ltd
Research Institute of Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Southern Power Grid Co Ltd, Research Institute of Southern Power Grid Co Ltd filed Critical China Southern Power Grid Co Ltd
Priority to CN202011334339.7A priority Critical patent/CN112330493B/en
Publication of CN112330493A publication Critical patent/CN112330493A/en
Application granted granted Critical
Publication of CN112330493B publication Critical patent/CN112330493B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067—Enterprise or organisation modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80—Management or planning
    • Y02P90/82—Energy audits or management systems therefor
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention discloses a method, a device and a storage medium for modeling and comprehensive analysis of an energy system. According to the transmission axiom of energy flow in the comprehensive energy system, a branch characteristic equation of each energy subnet in the energy transmission network is established; according to the energy transmission network
Figure DDA0002796718630000011
Transfer and conversion kinetic equation to derive
Figure DDA0002796718630000012
While transferring in transfer lines
Figure DDA0002796718630000013
Generalized watch of damageAnd establishing the electricity in the transmission process according to the branch characteristic equation
Figure DDA0002796718630000014
Is/are as follows
Figure DDA0002796718630000015
Loss calculation formula, pressure
Figure DDA0002796718630000016
Is/are as follows
Figure DDA0002796718630000017
Damage calculation formula and heat
Figure DDA0002796718630000018
Is/are as follows
Figure DDA0002796718630000019
A loss calculation formula; establishing an energy network equation set of the energy transmission network, solving state quantities of all nodes in the energy transmission network according to the energy network equation set, and combining the state quantities
Figure DDA00027967186300000110
A loss calculation formula for each strand in the energy transmission network based on the principle of thermal economy
Figure DDA00027967186300000111
The system economy and energy conservation of the stream are evaluated, and system-related parameters are improved, so that the waste of energy and cost is reduced.

Description

Energy system modeling and comprehensive analysis method, device and storage medium
Technical Field
The invention relates to the field of energy and the field of thermodynamics, in particular to a method, a device and a storage medium for modeling and comprehensive analysis of an energy system.
Background
Under the background of energy shortage and more serious environmental problems, the comprehensive utilization of various forms of energy becomes a necessary trend, and in order to realize the comprehensive unification of various forms of energy systems, it is necessary to research the transmission and conversion rules of various energy flows from the basic axiom of energy transmission and conversion; at present, in the aspects of energy efficiency analysis, economic analysis and the like of a multi-energy flow system, the energy angle and the independent analysis of each energy subsystem are mostly limited, the energy saving performance and the economic performance of the system are not analyzed and evaluated simultaneously, so that related parameters of the energy system are set unreasonably, and energy and cost are wasted.
Disclosure of Invention
The embodiment of the invention provides a method, a device and a storage medium for modeling and comprehensive analysis of an energy system, which can simultaneously consider the energy saving performance and the economical efficiency of the system, evaluate the energy saving potential of the system and provide a reasonable improvement scheme for the related parameter setting of the system, thereby reducing the waste of energy and cost.
The embodiment of the invention provides a method for modeling and comprehensive analysis of an energy system, which comprises the following steps:
according to the transmission axiom of energy flow in an energy system, establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system;
according to the energy transmission network
Figure BDA0002796718610000011
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure BDA0002796718610000012
Is/are as follows
Figure BDA0002796718610000013
Loss calculation formula, pressure
Figure BDA0002796718610000014
Is/are as follows
Figure BDA0002796718610000015
Damage calculation formula and heat
Figure BDA0002796718610000016
Is/are as follows
Figure BDA0002796718610000017
A loss calculation formula;
establishing an energy network equation set of the energy transmission network, and solving state quantities of all nodes in the energy transmission network according to the state equation set;
according to the state quantities of all nodes and combining the electricity
Figure BDA00027967186100000211
Is/are as follows
Figure BDA00027967186100000212
Loss calculation formula, said pressure
Figure BDA00027967186100000213
Is/are as follows
Figure BDA00027967186100000214
Loss calculation formula and the heat
Figure BDA00027967186100000215
Is/are as follows
Figure BDA00027967186100000216
Loss calculation formula for said energy transmission network
Figure BDA00027967186100000217
And analyzing the system economy and the energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
Preferably, the establishing of the branch characteristic equation corresponding to the energy subnet in the energy transmission network of the energy system according to the transmission axiom of the energy flow in the energy system specifically includes:
according to the axiom of the transfer of energy flow in energy systems
Figure BDA0002796718610000021
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure BDA0002796718610000022
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure BDA0002796718610000023
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure BDA0002796718610000024
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure BDA0002796718610000025
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; branch characteristic equation by extensive flux
Figure BDA0002796718610000026
The integral in the length direction of the transfer line in the transmission network, A being the cross-sectional area of the cylindrical transfer line through which the extensive quantity flows, JiFor extended fluence, JiObtained from the transfer axiom; in the delivery common, FiIs the driving force for pushing the extensive amount to be transmitted, KiIs a wide extension chiiThe coefficient of transmission of (a) is,
Figure BDA0002796718610000027
is the strength magnitude gradient of the conjugate.
Preferably, said energy transmission network is
Figure BDA00027967186100000218
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure BDA00027967186100000219
Is/are as follows
Figure BDA00027967186100000220
Loss calculation formula, pressure
Figure BDA00027967186100000221
Is/are as follows
Figure BDA00027967186100000222
Damage calculation formula and heat
Figure BDA00027967186100000223
Is/are as follows
Figure BDA00027967186100000224
The loss calculation formula is specifically as follows:
according to the energy transmission network
Figure BDA00027967186100000225
Kinetic equations of transfer and conversion
Figure BDA0002796718610000028
Establishing
Figure BDA00027967186100000226
General formula of calculation
Figure BDA0002796718610000029
Wherein rho in the kinetic equation is medium density gxThe source intensity, χ, of the extensive amount χ in the unit volume of medium0A silence value that is an intensity amount x,
Figure BDA00027967186100000210
for the extended amount of transfer law, the left side of the equation is
Figure BDA00027967186100000227
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure BDA00027967186100000228
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure BDA00027967186100000229
The third term on the right represents the other forms of the transformation between
Figure BDA00027967186100000313
Converted into such a form
Figure BDA00027967186100000314
The above-mentioned
Figure BDA00027967186100000315
The loss is represented by ap in the generalized equation,
Figure BDA0002796718610000031
to increase the amount of spread in the delivery process,
Figure BDA0002796718610000032
represents an extensive flow;
according to the branch characteristic equation
Figure BDA0002796718610000033
Build up electricity
Figure BDA00027967186100000316
Is/are as follows
Figure BDA00027967186100000317
Formula for calculating loss
Figure BDA0002796718610000034
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure BDA0002796718610000035
is current, K is conductivity, RiIs the thermal resistance, l is the length of the cylindrical transfer conduit, in the electrical network, the value of voltage dead xe00, electricity
Figure BDA00027967186100000318
The loss is equal to the electric energy loss; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure BDA0002796718610000036
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure BDA0002796718610000037
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure BDA00027967186100000319
Is/are as follows
Figure BDA00027967186100000320
In the formula of loss calculation
Figure BDA0002796718610000038
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
according to electricity
Figure BDA00027967186100000321
Loss calculation process build-up pressure
Figure BDA00027967186100000322
Is/are as follows
Figure BDA00027967186100000323
Formula for calculating loss
Figure BDA0002796718610000039
WhereinIn a fluid network, the value of the pressure energy in the silent state is zero, χpA,χpEPressure of a section A and a section E respectively, wherein RpIs the flow resistance according to the Navier-Stokes equation
Figure BDA00027967186100000310
Calculation, p is the fluid density, kpFor an extensive amount of transfer coefficient in the fluid network,
Figure BDA00027967186100000311
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
heat of formation
Figure BDA00027967186100000324
Is/are as follows
Figure BDA00027967186100000325
Formula for calculating loss
Figure BDA00027967186100000312
Wherein: the heat
Figure BDA00027967186100000326
Is/are as follows
Figure BDA00027967186100000327
X in the damage calculation formulah0The value is a silent state value of temperature, an environment temperature value is usually taken, and the temperature value at the tail end of a transmission pipeline is obtained by a Suhoff temperature drop formula:
Figure BDA0002796718610000041
λhis the heat transfer coefficient of the transfer line; in a thermal network, an incompressible fluid is used as a heat transfer medium, the heat being
Figure BDA00027967186100000411
Is/are as follows
Figure BDA00027967186100000412
The formula of the loss is used as the formula of entropy increaseThe formula is combined with a heat energy loss calculation formula, wherein the entropy increase calculation formula is
Figure BDA0002796718610000042
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure BDA0002796718610000043
as flow rate of heat transfer medium, χhA,χhETemperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure BDA0002796718610000044
preferably, the establishing an energy network equation set of the energy transmission network, and solving the state quantities of all nodes in the energy transmission network according to the state equation set specifically include:
establishing the energy transmission network topology constraint equation set
Figure BDA0002796718610000045
Equation of the branch characteristic
Figure BDA0002796718610000046
Establishing an energy network equation set of the energy transmission network in combination with the topological constraint equation set, and solving the energy network equation set to obtain state quantities of all nodes in the energy transmission network;
wherein A in the topological constraint equation set is a correlation matrix, and BfIn the form of a matrix of elementary loops,
Figure BDA0002796718610000047
is a wide-spread flux matrix, Δ χiIs an intensity quantity difference matrix; when the branch characteristic equation represents an electric network equation, x is voltage,
Figure BDA0002796718610000048
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; the branchWhen the characteristic equation represents a fluid network equation, χ is pressure,
Figure BDA0002796718610000049
is volume flow, K is volume transfer coefficient rate, when the branch characteristic equation expresses a heat network equation, x is temperature,
Figure BDA00027967186100000410
flow rate of heat transfer medium, K is the heat transfer coefficient of the transfer line, l is the length of the cylindrical transfer line, KeFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows.
Preferably, said combining of said electricity and said according to said all-node state quantities
Figure BDA00027967186100000413
Is/are as follows
Figure BDA00027967186100000414
Loss calculation formula, said pressure
Figure BDA00027967186100000415
Is/are as follows
Figure BDA00027967186100000416
Loss calculation formula and the heat
Figure BDA00027967186100000417
Is/are as follows
Figure BDA00027967186100000418
Loss calculation formula for said energy transmission network
Figure BDA00027967186100000419
Analyzing the system economy and energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result specifically comprises the following steps:
according to the state quantity of all nodes and the electricity
Figure BDA0002796718610000052
Is/are as follows
Figure BDA0002796718610000053
Loss calculation formula, said pressure
Figure BDA0002796718610000054
Is/are as follows
Figure BDA0002796718610000055
Loss calculation formula and the heat
Figure BDA0002796718610000056
Is/are as follows
Figure BDA0002796718610000057
Calculation of each node in network by loss calculation formula
Figure BDA0002796718610000058
Flow number, obtain unit economic cost, and unit non-energy cost
Figure BDA0002796718610000059
Flow conversion, calculating
Figure BDA00027967186100000510
The cost is reduced;
calculating technical economic coefficient
Figure BDA0002796718610000051
Wherein D isxRepresenting energy transfer and conversion processes
Figure BDA00027967186100000511
Loss value, CDxIs composed of
Figure BDA00027967186100000512
The cost is taken as input
Figure BDA00027967186100000513
Average of (2)
Figure BDA00027967186100000514
Cost, Z stands for
Figure BDA00027967186100000515
A non-energy cost of stream numerical translation, the fractional energy cost comprising equipment cost, labor cost, and operational cost;
optimizing the energy transmission network according to the technical economic coefficient to make the technical economic coefficient equal to a preset threshold value N,
Figure BDA00027967186100000516
and the ratio of the loss cost to the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to the optimization result.
The invention discloses a method for modeling and comprehensive analysis of an energy system, which comprises the steps of establishing a branch characteristic equation of each energy subnet in an energy transmission network according to the transmission axiom of energy flow in the comprehensive energy system; according to the energy transmission network
Figure BDA00027967186100000517
A kinetic equation of transfer and conversion is carried out, and electricity in the transfer process is established according to the branch characteristic equation
Figure BDA00027967186100000518
Is/are as follows
Figure BDA00027967186100000519
Loss calculation formula, pressure
Figure BDA00027967186100000520
Is/are as follows
Figure BDA00027967186100000521
Damage calculation formula, heat
Figure BDA00027967186100000522
Is/are as follows
Figure BDA00027967186100000523
A loss calculation formula; the energy network equation set of the energy transmission network is established, all node state quantities in the energy transmission network are solved according to the system state equation, and the state quantities are combined
Figure BDA00027967186100000524
A loss calculation formula for each strand in the energy transmission network
Figure BDA00027967186100000525
Evaluating system economy and energy conservation of the flow according to the network topology constraint equations
Figure BDA00027967186100000526
And adjusting the equipment parameters in the energy transmission network to enable the f to be in accordance with the influence of the parameters in the energy transmission network in the loss calculation equation on the technical economic coefficientexEqual to a preset threshold value N, is set,
Figure BDA00027967186100000527
the proportion of the loss cost and the non-energy cost reaches the most reasonable distribution of energy conservation and economy, and the waste of energy and cost is reduced.
The embodiment of the invention also provides a device for modeling and comprehensive analysis of an energy system, which comprises: a branch characteristic equation calculation module,
Figure BDA00027967186100000528
The system comprises a loss calculation module, a node state calculation module and a parameter optimization module;
the branch characteristic equation calculation module is used for establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system according to the transmission axiom of the energy flow in the energy system;
the above-mentioned
Figure BDA00027967186100000529
The loss calculation module is used for calculating the loss according to the energy transmission network
Figure BDA00027967186100000530
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure BDA00027967186100000531
Is/are as follows
Figure BDA00027967186100000532
Loss calculation formula, pressure
Figure BDA00027967186100000533
Is/are as follows
Figure BDA00027967186100000534
Damage calculation formula and heat
Figure BDA00027967186100000610
Is/are as follows
Figure BDA00027967186100000611
A loss calculation formula;
the node state calculation module is used for establishing an energy network equation set of the energy transmission network and solving state quantities of all nodes in the energy transmission network according to the state equation set;
the parameter optimization module is used for combining the electricity according to the state quantities of all the nodes
Figure BDA00027967186100000612
Is/are as follows
Figure BDA00027967186100000613
Loss calculation formula, said pressure
Figure BDA00027967186100000614
Is/are as follows
Figure BDA00027967186100000615
Loss calculation formula and the heat
Figure BDA00027967186100000616
Is/are as follows
Figure BDA00027967186100000617
Loss calculation formula for said energy transmission network
Figure BDA00027967186100000618
And analyzing the system economy and the energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
Preferably, the branch characteristic equation calculation module specifically functions as:
according to the axiom of the transfer of energy flow in energy systems
Figure BDA0002796718610000061
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure BDA0002796718610000062
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure BDA0002796718610000063
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a steady-state incompressible fluid network equation, x is pressure intensity,
Figure BDA0002796718610000064
is the volume flow, K is the volume transfer coefficient rate; branch characteristic equation by extensive flux
Figure BDA0002796718610000065
The integral in the length direction of the transfer line in the transmission network, A being the cross-sectional area of the cylindrical transfer line through which the extensive quantity flows, JiFor extended fluence, obtained from the transport axiom, FiIs the driving force for pushing the extensive amount to be transmitted, KiIs a wide extension chiiThe coefficient of transmission of (a) is,
Figure BDA0002796718610000066
is a conjugate intensity magnitude gradient;
the above-mentioned
Figure BDA00027967186100000619
The loss calculation module comprises
Figure BDA00027967186100000620
Loss of generality computing unit, electricity
Figure BDA00027967186100000621
Loss calculation unit, voltage
Figure BDA00027967186100000622
Loss calculation unit and heat
Figure BDA00027967186100000623
A loss calculation unit;
the above-mentioned
Figure BDA00027967186100000624
The loss-to-prevalence calculation unit is used for calculating the loss-to-prevalence in the energy transmission network
Figure BDA00027967186100000625
Kinetic equations of transfer and conversion
Figure BDA0002796718610000067
Establishing
Figure BDA00027967186100000626
General formula of calculation
Figure BDA0002796718610000068
Wherein rho in the kinetic equation is medium density gxThe source intensity, χ, of the extensive amount χ in the unit volume of medium0A silence value that is an intensity amount x,
Figure BDA0002796718610000069
the transfer rule is extensive quantity; the left side of the equation is
Figure BDA00027967186100000627
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure BDA00027967186100000713
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure BDA00027967186100000714
The third term on the right represents the other forms of the transformation between
Figure BDA00027967186100000715
Converted into such a form
Figure BDA00027967186100000716
The above-mentioned
Figure BDA00027967186100000717
The loss is represented by ap in the generalized equation,
Figure BDA0002796718610000071
to increase the amount of spread in the delivery process,
Figure BDA0002796718610000072
represents an extensive flow;
the electricity
Figure BDA00027967186100000718
The loss calculation unit is used for calculating the loss according to the branch characteristic equation
Figure BDA0002796718610000073
Build up electricity
Figure BDA00027967186100000719
Is/are as follows
Figure BDA00027967186100000720
Formula for calculating loss
Figure BDA0002796718610000074
Wherein, in the electric network, the voltage has a silent state value χe00, electricity
Figure BDA00027967186100000721
The loss is equal to the electric energy loss, and when the branch characteristic equation represents an electric network equation, x is voltage,
Figure BDA0002796718610000075
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure BDA0002796718610000076
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure BDA0002796718610000077
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure BDA00027967186100000722
Is/are as follows
Figure BDA00027967186100000723
In the formula of loss calculation
Figure BDA0002796718610000078
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
said pressure
Figure BDA00027967186100000724
For loss calculating unitsIn accordance with electricity
Figure BDA00027967186100000725
Loss calculation process build-up pressure
Figure BDA00027967186100000726
Is/are as follows
Figure BDA00027967186100000727
Formula for calculating loss
Figure BDA0002796718610000079
Wherein, in the fluid network, the silence state value of the pressure energy is also zero, RpIs the flow resistance, chipA,χpEPressure of section A and section E, respectively, according to the Navier-Stokes equation
Figure BDA00027967186100000710
Calculation, p is the fluid density, kpFor an extensive amount of transfer coefficient in the fluid network,
Figure BDA00027967186100000711
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
the heat
Figure BDA00027967186100000728
Loss calculation unit for establishing heat
Figure BDA00027967186100000729
Is/are as follows
Figure BDA00027967186100000730
Loss calculation formula heat
Figure BDA00027967186100000731
The loss is calculated by
Figure BDA00027967186100000712
Wherein: in thermal networks, the use ofCompressed fluid as a heat transfer medium, said heat
Figure BDA00027967186100000732
Is/are as follows
Figure BDA00027967186100000733
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure BDA0002796718610000081
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure BDA0002796718610000087
as flow rate of heat transfer medium, χhA,χhETemperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure BDA0002796718610000082
the heat
Figure BDA0002796718610000088
Is/are as follows
Figure BDA0002796718610000089
X in the damage calculation formulah0The value is a silent state value of temperature, an environment temperature value is usually taken, and the temperature value at the tail end of a transmission pipeline is obtained by a Suhoff temperature drop formula:
Figure BDA0002796718610000083
λhis the heat transfer coefficient of the transfer line;
the node state calculating unit has the functions of: establishing a topological constraint equation set of the energy transmission network:
Figure BDA0002796718610000084
equation of the branch characteristic
Figure BDA0002796718610000085
With said topological constraint equationEstablishing an energy network equation set of the energy transmission network in a combined manner, and solving the equation set to obtain state quantities of all nodes of the energy transmission network; wherein A in the topological constraint equation set is a correlation matrix, and BfIn the form of a matrix of elementary loops,
Figure BDA00027967186100000826
the matrix is an extensive flow matrix, and the delta chi i is an intensity quantity difference matrix; when the branch characteristic equation represents an electric network equation, x is voltage,
Figure BDA00027967186100000810
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure BDA00027967186100000811
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure BDA00027967186100000812
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure BDA00027967186100000813
Is/are as follows
Figure BDA00027967186100000814
In the formula of loss calculation
Figure BDA0002796718610000086
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
the parameter optimization module comprises
Figure BDA00027967186100000815
Loss cost calculation unit and technical economyCoefficient calculation unit and parameter adjustment unit:
the above-mentioned
Figure BDA00027967186100000816
The loss cost calculation unit is used for calculating the loss cost according to the state quantity of all the nodes and the electricity
Figure BDA00027967186100000817
Is/are as follows
Figure BDA00027967186100000818
Loss calculation formula, said pressure
Figure BDA00027967186100000819
Is/are as follows
Figure BDA00027967186100000820
Loss calculation formula and the heat
Figure BDA00027967186100000821
Is/are as follows
Figure BDA00027967186100000822
Calculation of each node in network by loss calculation formula
Figure BDA00027967186100000823
Flow number, non-energy cost in units
Figure BDA00027967186100000824
Flow conversion, calculating
Figure BDA00027967186100000825
The cost is reduced;
the technical economic coefficient calculating unit is used for calculating the technical economic coefficient
Figure BDA0002796718610000091
Wherein D isxRepresenting energy transfer and conversion processes
Figure BDA0002796718610000092
Loss value, CDxIs composed of
Figure BDA0002796718610000093
The cost is taken as input
Figure BDA0002796718610000094
Average of (2)
Figure BDA0002796718610000095
Cost, Z stands for
Figure BDA0002796718610000096
Non-energy cost of stream value translation;
the parameter adjusting unit is used for optimizing the energy transmission network according to the technical economic coefficient to ensure that the technical economic coefficient is equal to a preset threshold value N,
Figure BDA0002796718610000097
and the ratio of the loss cost to the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to the optimization result.
The embodiment of the present invention further provides an apparatus for modeling and comprehensive analysis of an energy system, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, the method for modeling and comprehensive analysis of an energy system described in any of the above embodiments is implemented.
The embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, where when the computer program runs, a device in which the computer-readable storage medium is located is controlled to execute the method for modeling and comprehensive analysis of an energy system according to any of the above embodiments.
The invention discloses a method, a device and a storage medium for modeling and comprehensive analysis of an energy system, which can establish energy according to the transmission axiom of energy flow in the comprehensive energy systemThe branch characteristic equation of each energy sub-network in the source transmission network is determined according to the branch characteristic equation in the energy transmission network
Figure BDA0002796718610000098
A kinetic equation of transfer and conversion is established according to the branch characteristic equation in the transfer process
Figure BDA0002796718610000099
A loss calculation formula; establishing an energy network equation set of an energy transmission network, solving state quantities of all nodes in the energy transmission network according to the system state equation, and combining the state quantities
Figure BDA00027967186100000910
A loss calculation formula for each strand in the energy transmission network
Figure BDA00027967186100000911
Evaluating system economy and energy conservation of the flow according to the network topology constraint equations
Figure BDA00027967186100000912
And adjusting the equipment parameters in the energy transmission network to enable the f to be in accordance with the influence of the parameters in the energy transmission network in the loss calculation equation on the technical economic coefficientexEqual to a preset threshold value N, is set,
Figure BDA00027967186100000913
the proportion of the loss cost and the non-energy cost reaches the most reasonable distribution of energy conservation and economy, and the waste of energy and cost is reduced.
Drawings
Fig. 1 is a schematic flow chart of a method for modeling and comprehensive analysis of an energy system according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a cylindrical transfer line energy transfer provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a physical model of energy transfer with constant coefficient of extensive mass transfer provided by an embodiment of the present invention;
FIG. 4 is a schematic diagram of a physical model of energy transfer with constant coefficient of extensive mass transfer provided by an embodiment of the present invention;
fig. 5 is a schematic configuration diagram of an electricity-cold cogeneration network system according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of an apparatus for modeling and comprehensive analysis of an energy system according to an embodiment of the present invention
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a schematic flow chart of a method for modeling and comprehensive analysis of an energy system according to an embodiment of the present invention.
The embodiment of the invention provides a method for modeling and comprehensive analysis of an energy system, which comprises the following steps from S101 to S104:
s101, establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of an energy system according to the transmission axiom of energy flow in the energy system;
s102, according to the energy source in the transmission network
Figure BDA0002796718610000101
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure BDA0002796718610000102
Is/are as follows
Figure BDA0002796718610000103
Loss calculation formula, pressure
Figure BDA0002796718610000104
Is/are as follows
Figure BDA0002796718610000105
Damage calculation formula and heat
Figure BDA0002796718610000106
Is/are as follows
Figure BDA0002796718610000107
A loss calculation formula;
s103, establishing an energy network equation set of the energy transmission network, and solving state quantities of all nodes in the energy transmission network according to the state equation set;
s104, combining the electricity according to the state quantities of all the nodes
Figure BDA0002796718610000112
Is/are as follows
Figure BDA0002796718610000113
Loss calculation formula, said pressure
Figure BDA0002796718610000114
Is/are as follows
Figure BDA0002796718610000115
Loss calculation formula and the heat
Figure BDA0002796718610000116
Is/are as follows
Figure BDA0002796718610000117
Loss calculation formula for said energy transmission network
Figure BDA0002796718610000118
And analyzing the system economy and the energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
The embodiment of the invention provides a method for modeling and comprehensive analysis of an energy system, which is a transmission public according to the energy flow transmission rule in the comprehensive energy systemEstablishing branch characteristic equations in energy flow subnetworks in the energy transmission network based on the transmission characteristics of different energy types; in energy-based networks
Figure BDA0002796718610000119
Transfer and conversion kinetic equation to derive
Figure BDA00027967186100001110
While transferring in transfer lines
Figure BDA00027967186100001111
Generalized calculation formula of loss and specifically gives out electricity, voltage and heat in the transfer process
Figure BDA00027967186100001112
Is/are as follows
Figure BDA00027967186100001113
A loss calculation formula; according to the theoretical knowledge of heat economy, the energy resources of each strand in the network
Figure BDA00027967186100001114
The method can simultaneously consider the energy conservation performance and the economical efficiency of the system, evaluate the energy conservation potential of the system, and propose an improved scheme for the relevant parameters of the system, thereby reducing the waste of energy and cost.
In another preferred embodiment, step S101 is specifically:
according to the principle of energy flow transfer, the transfer law of the extensive quantity corresponding to the ith energy form in the comprehensive energy system can be represented by the following formula:
Figure BDA0002796718610000111
wherein, JiFor extended fluence, FiIs the driving force for pushing the extensive amount to be transmitted, KiIs a wide extension chiiThe coefficient of transmission of (a) is,
Figure BDA00027967186100001115
is the strength magnitude gradient of the conjugate.
In addition, in any form of energy flow transfer process, the form of the transfer is different, but the same law is followed in nature, and the transfer axiom summarizes the physical nature of various transfer processes, namely, the flow transportation of extensive amount is realized under the push of corresponding intensity amount difference, so that various forms of energy flow transfer are completed.
The energy transfer is shown by a formula, and the transfer rule of the extensive quantity is determined by two factors, namely a conjugate intensity quantity gradient which is the root cause for generating basic extensive quantity on the one hand, and a transfer coefficient K on the other handiIt is also important for the delivery process of extensive amounts.
Referring to fig. 2, a schematic energy transfer diagram of a cylindrical transfer pipe according to an embodiment of the present invention is shown, in which the various energy transmission networks are generally formed by cylindrical transfer pipes, which are radially sealed and have an extension that flows only in the axial direction from section a to section E, i.e., JiAFlow direction JiECombined with transmission line pipe, can further derive extensive flow
Figure BDA0002796718610000125
Expression (c):
Figure BDA0002796718610000121
wherein, A is the sectional area of the cylindrical transfer pipeline through which the extensive quantity flows, and the formula is integrated along the length direction of the transfer pipeline, so that the relational expression of the strength quantity and the extensive quantity can be obtained.
For different forms of energy, an extensive amount of transfer coefficient KiMay or may not be changed, and is therefore divided intoTwo cases are discussed.
Referring to fig. 3, it is a schematic diagram of an energy transfer physical model when the extensive quantity transfer coefficient is constant, provided by an embodiment of the present invention, and when the extensive quantity transfer coefficient K is constantiFor a fixed value, such as the transmission of electric energy, the transmission of pressure energy in steady laminar flow, and the like, the equivalent transmission equation of the extensive quantity in the cylindrical transmission pipeline is as follows:
Figure BDA0002796718610000122
wherein R isiIs constant and l is the length of the cylindrical transfer line.
Figure BDA0002796718610000123
Is an integral
Figure BDA0002796718610000124
The Lagrange median value of (1), i.e. the equivalent extensive flux in the transfer process, the value of which can be determined by the difference x of the strength of two ends of the lineiA-χiEAnd a constant value RiAnd (4) obtaining.
Referring to fig. 4, it is a schematic diagram of an energy transfer physical model when the extensive quantity transfer coefficient is constant, provided by an embodiment of the present invention, and when the extensive quantity transfer coefficient K is constantiWhen not of constant value, but of coefficient of transfer of energy λi＝KiχiFor fixed values, such as a one-dimensional constant heat transfer process, the equivalent transfer equation of the spread in the cylindrical transfer pipe is:
Figure BDA0002796718610000131
wherein R isiFor thermal resistance, l is the length of the cylindrical transfer line. At this time
Figure BDA0002796718610000134
The relation between the X and the X is non-linear,
Figure BDA0002796718610000135
the relationship with ln χ is linear. In the heat transfer in the actual engineering, the temperature change along the axial direction is small, the heat transfer caused by the small temperature change can be ignored, and the heat transfer along the radial direction causes main heat loss in the transfer process, and the following conditions are adopted:
Figure BDA0002796718610000132
wherein r1 is the radius of the inner wall of the cylindrical transfer pipeline, and r2 is the radius of the outer wall; phIs the heat flow, i.e. the thermal power, RhRepresenting the thermal resistance of the radial heat conduction. If the heat convection process between the fluid inside and outside the pipe and the pipe wall is considered, several corresponding thermal resistances can be connected in series to serve as the total thermal resistance of the radial heat conduction process, and the expression is as follows:
Figure BDA0002796718610000133
in the formula, λhdIs the coefficient of thermal conductivity; lambda [ alpha ]hc1Is the convective heat transfer coefficient of the inner wall, λhc2Is the convective heat transfer coefficient of the outer wall.
Through the analysis, the equivalent transfer equations of the three energies of the electric energy, the pressure energy and the heat energy in the steady laminar flow state in the transfer pipeline have similar expression forms, namely the extension amounts are all linear responses of the medium to the difference of the external strength. Because of the inconstant thermal energy transfer coefficient, the extensive quantity in the thermal energy transfer process is converted into thermal energy, and the equivalent transfer equations of the electric network and the steady-state incompressible fluid network have strong similarity, so that the branch characteristic equations of each energy sub-network in the time-invariant energy transmission network can be uniformly described as follows:
Figure BDA0002796718610000141
when the equation represents an electrical network equation, χ is voltage,
Figure BDA0002796718610000145
is current, K is conductivity;
when the equation represents a fluid network equation, χ is the pressure,
Figure BDA0002796718610000146
is the volume flow, K is the volume transfer coefficient rate;
when the formula represents a heat network equation, χ is temperature,
Figure BDA0002796718610000147
the heat transfer medium flow rate, K, is the transfer line heat transfer coefficient.
In another preferred embodiment, step S103 is specifically:
Figure BDA0002796718610000148
while transferring in transfer lines
Figure BDA0002796718610000149
The generalized expression derivation process for impairments is as follows:
Figure BDA00027967186100001410
is the significant portion of energy, i.e., the portion of energy that can be converted to useful work to the greatest extent, representing the "mass" of energy. In comparison with the amount of energy,
Figure BDA00027967186100001411
the method can further reflect the essence of the work-doing capability loss in the process of energy transmission and conversion and reflect the external work-doing capability of the energy, so the method is suitable for being used as the standard for evaluating the energy-saving benefit of the comprehensive energy system.
In the energy transmission network
Figure BDA00027967186100001412
The kinetic equations for transfer and conversion are:
Figure BDA0002796718610000142
where ρ is the density of the medium, gxIs the source intensity of the etendue χ in a unit volume of medium,
Figure BDA0002796718610000143
the transfer rule of extensive quantity.
The left side of the equation is
Figure BDA00027967186100001413
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure BDA00027967186100001414
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure BDA00027967186100001415
The third term on the right represents the other forms of the transformation between
Figure BDA00027967186100001416
Converted into such a form
Figure BDA00027967186100001417
Electricity in the process of transmission
Figure BDA00027967186100001418
Is/are as follows
Figure BDA00027967186100001419
Loss calculation formula, pressure
Figure BDA00027967186100001420
Is/are as follows
Figure BDA00027967186100001421
Loss calculationFormula, heat
Figure BDA00027967186100001422
Is/are as follows
Figure BDA00027967186100001423
The loss calculation formula is derived as follows:
for the time-invariant energy system, the specific transmission pipeline is combined, and the volume can be obtained through volume division
Figure BDA00027967186100001424
The generalized calculation of the loss is:
Figure BDA0002796718610000144
wherein, Δ P represents an energy loss,
Figure BDA00027967186100001425
for the increase of the spread during the transfer, χ0Is a silence value of the intensity amount χ.
In an electrical network, the value of voltage dead xe00, so that
Figure BDA0002796718610000158
The loss is equal to the electric energy loss, and the electricity is obtained by a branch characteristic equation
Figure BDA0002796718610000159
The loss is calculated by the formula:
Figure BDA0002796718610000151
wherein
Figure BDA0002796718610000152
keRepresents the electrical conductivity;
in a thermal network, when an incompressible fluid is used as the heat transfer medium, the formula for the entropy increase in the transfer line is:
Figure BDA0002796718610000153
rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure BDA00027967186100001510
is the heat transfer medium flow rate. The heat energy loss calculation formula is as follows:
Figure BDA0002796718610000154
will be transferred in the pipeline
Figure BDA00027967186100001511
The expression of the transmission loss, the formula of the entropy increase calculation and the formula of the heat energy loss are combined to obtain the heat
Figure BDA00027967186100001512
The loss is calculated by the formula:
Figure BDA0002796718610000155
wherein, χh0The temperature value is a temperature silent state value, an environment temperature value is usually taken, and the temperature value at the tail end of a transmission pipeline can be obtained by a Suhoff temperature drop formula:
Figure BDA0002796718610000156
wherein λ ishIs the heat transfer coefficient of the transfer line.
In the incompressible steady-state laminar flow network, the reference voltage is zero because the silence state value of the pressure energy is also zero
Figure BDA00027967186100001513
Process of calculating damage, pressure
Figure BDA00027967186100001514
The expression for the loss is:
Figure BDA0002796718610000157
wherein R ispIs the flow resistance, which can be obtained from the Navier-Stokes equation
Figure BDA0002796718610000161
In another preferred embodiment, step S103 is specifically:
in order to describe the inherent relationship between the node intensity amount and the loop extensive flow in the energy transmission network, an energy transmission network topology constraint equation set is established:
Figure BDA0002796718610000162
wherein A is a correlation matrix, BfIn the form of a matrix of elementary loops,
Figure BDA00027967186100001620
is a wide-spread flux matrix, Δ χiIs an intensity quantity difference matrix;
combining the branch characteristic equation with an energy transmission network topology constraint equation set to establish an energy network equation set of an energy transmission network, and solving the energy network equation set to obtain state quantities of all nodes of the energy system;
in another preferred embodiment, step S104 specifically includes:
according to the node state quantity and electricity
Figure BDA0002796718610000164
Press and press
Figure BDA0002796718610000165
Heat generation
Figure BDA0002796718610000166
Is/are as follows
Figure BDA0002796718610000167
Calculation of loss equation for each node in network
Figure BDA0002796718610000168
Flow number, non-energy cost in units
Figure BDA0002796718610000169
Converting the flow to calculate the electric quantity
Figure BDA00027967186100001610
Cost and cold capacity
Figure BDA00027967186100001611
Cost;
by calculating the technical economic coefficient fex：
Figure BDA0002796718610000163
DxRepresenting energy transfer and conversion processes
Figure BDA00027967186100001612
Loss value, CDxIs composed of
Figure BDA00027967186100001613
The cost is taken as input
Figure BDA00027967186100001614
Average of (2)
Figure BDA00027967186100001615
Cost, Z stands for
Figure BDA00027967186100001616
Non-energy cost of stream value translation;
f isexCan reflect
Figure BDA00027967186100001617
Proportional relation between loss cost and non-energy cost, and according to the network topology constraint equation set
Figure BDA00027967186100001618
And adjusting the equipment parameters in the energy transmission network to enable the f to be in accordance with the influence of the parameters in the energy transmission network in the loss calculation equation on the technical economic coefficientexEqual to a preset threshold value N, is set,
Figure BDA00027967186100001619
the proportion of the loss cost to the non-energy cost reaches the most reasonable distribution of economy and energy conservation.
Referring to fig. 5, which is a schematic configuration diagram of an electricity-cold cogeneration network system according to an embodiment of the present invention, as shown in fig. 5, natural gas purchased through a gas turbine is combusted in the electricity-cold cogeneration network system to generate power to supply power to a complex building, an experimental building and a new building, a part with insufficient power is purchased from a power grid, in addition, high-temperature waste heat after combustion of the natural gas is utilized, cooling is supplied to the new building and the complex building through a waste heat recovery refrigeration device and a lithium bromide absorption refrigerator, and a part with insufficient cooling capacity is provided by an electric refrigeration air conditioner.
From the perspective of network topology, the inherent relationship between the node strength amount and the loop spread flow in each energy transmission network is described, so that an energy transmission network topology constraint equation system can be established:
Figure BDA0002796718610000171
wherein A is a correlation matrix, BfIn the form of a matrix of elementary loops,
Figure BDA00027967186100001717
is a wide-spread flux matrix, Δ χiIs an intensity quantity difference matrix.
Forming an energy network equation set of the energy transmission network by the branch characteristic equation set and the energy transmission network topology constraint equation set, and solving the equation set to obtain all state quantities of the energy system;
according to the theory of heat economy, calculating the electricity for supplying power and cooling for the gas turbine and the waste heat recovery unit
Figure BDA0002796718610000172
Cost and coldness
Figure BDA0002796718610000173
Cost:
cgasEx,gas+ZHP＝ceExe+chExh
chExh+Zc＝ccExc
wherein c represents
Figure BDA0002796718610000174
Cost per unit of heat economy of flow, ExRepresents
Figure BDA0002796718610000175
Flow number, ZHPNon-energy costs (equipment costs and labor costs of gas turbines, etc.) representing the cogeneration link, ZcRepresenting the non-energy costs (equipment costs and operating costs, etc.) of the waste heat recovery refrigeration equipment. Subscript gasE, h and c respectively represent natural gas and electricity
Figure BDA0002796718610000176
Heat generation
Figure BDA0002796718610000178
And cool
Figure BDA0002796718610000177
By all state quantities of the energy system and
Figure BDA0002796718610000179
loss calculation formula, non-energy cost as a unitBit
Figure BDA00027967186100001710
Converting the flow to calculate the electric quantity
Figure BDA00027967186100001711
Cost and cold capacity
Figure BDA00027967186100001712
And (4) cost.
Similarly, the cold quantity generated by refrigeration of the electric air conditioner in the process of purchasing electricity through the power grid is calculated
Figure BDA00027967186100001713
Cost:
ce,gridEx,grid+Zec＝cc2Exc2
wherein Z isecRepresenting the equipment cost of the electric air conditioner. Subscript grid represents the power grid purchase, and subscript 2 is used for cooling two production processes
Figure BDA00027967186100001714
The difference can also be calculated from the cold produced in this production process
Figure BDA00027967186100001715
And (4) cost.
Cost of purchasing electricity from power gride,gridCan be priced by local electricity and generate electricity from natural gas
Figure BDA00027967186100001716
Per unit economic cost ceThe cost c of purchasing electricity from the power gride,gridBy contrast, the two modes of electricity production can be compared
Figure BDA0002796718610000182
The economy of (2).
By a technical economic factor fexCan reflect
Figure BDA0002796718610000183
Loss of expense and incorgulationThe proportional relationship of the amount and the cost is shown as the following formula:
Figure BDA0002796718610000181
in the formula, DxRepresenting processes of energy transfer, conversion
Figure BDA0002796718610000184
The value of the loss is reduced,
Figure BDA00027967186100001821
is composed of
Figure BDA0002796718610000185
The cost is taken as input
Figure BDA0002796718610000186
Average of (2)
Figure BDA0002796718610000187
Cost, Z stands for
Figure BDA0002796718610000188
Flow value converted non-energy cost. f. ofexWhen the temperature is too high, the temperature is high,
Figure BDA0002796718610000189
the loss cost is low, but the non-energy cost investment is too large and is not economic enough; f. ofexIf too low, the non-energy costs are low, but concomitantly
Figure BDA00027967186100001810
The loss cost is too high, and energy is not saved; when f isexWhen the content is equal to 0.5,
Figure BDA00027967186100001811
the loss and non-energy costs are 1:1, considered the most reasonable distribution at this time, and by this parameter, equipment replacement and capital allocation adjustments in the network can be considered.
Waste heat recovery refrigerating equipmentWhen in cold supply, cold water pipelines are laid at cold load positions for cooling
Figure BDA00027967186100001812
Is transmitted because
Figure BDA00027967186100001813
Irreversibility of the transfer process will inevitably occur
Figure BDA00027967186100001814
To damage this part
Figure BDA00027967186100001815
The cost is calculated, and the transmission process is
Figure BDA00027967186100001816
The damage formula is given in S102 as a function of the pipe length L, the pipe diameter d and the pipe material epsilon, and the non-energy cost of the pipe is also a function of the pipe length L, the pipe diameter d and the pipe material epsilon, so that the technical-economic coefficient of the pipe can be expressed as fex(L, d, epsilon) which can be used to guide the parameters of the laid pipeline to achieve the most rational distribution of economy and energy saving.
The invention discloses a method for modeling and comprehensive analysis of an energy system, which comprises the steps of establishing a branch characteristic equation of each energy subnet in an energy transmission network according to the transmission axiom of energy flow in the comprehensive energy system; according to the energy transmission network
Figure BDA00027967186100001817
Transfer and conversion kinetics equation, and establishing the transfer process according to the branch characteristic equation
Figure BDA00027967186100001818
Establishing an energy network equation set of the energy transmission network by using a loss calculation formula, solving state quantities of all nodes in the energy transmission network according to the state equation set, and combining the state quantities
Figure BDA00027967186100001819
A loss calculation formula for each strand in the energy transmission network
Figure BDA00027967186100001820
The system economy and energy conservation of the stream are evaluated, and system-related parameters are improved, so that the waste of energy and cost is reduced.
Referring to fig. 6, which is a schematic diagram of an apparatus for modeling and comprehensive analysis of an energy system according to an embodiment of the present invention, the apparatus for modeling and comprehensive analysis of an energy system according to the embodiment of the present invention includes a branch characteristic equation calculation module,
Figure BDA0002796718610000191
The system comprises a loss calculation module, a node state calculation module and a parameter optimization module.
In a specific implementation, the apparatus for modeling and comprehensive analysis of an energy system can complete specific functions of the method for modeling and comprehensive analysis of an energy system provided in any one of the above embodiments, and a specific implementation process is specifically described in any one of the embodiments of the method for modeling and comprehensive analysis of an energy system, which is not described in detail in this embodiment.
The embodiment of the present invention further provides an apparatus for modeling and comprehensive analysis of an energy system, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, the method for modeling and comprehensive analysis of an energy system described in any of the above embodiments is implemented.
The device for modeling and comprehensive analysis of the energy system can be a desktop computer, a notebook computer, a palm computer, a cloud server and other computing equipment. The device/terminal equipment for modeling and comprehensive analysis of the energy system can comprise, but is not limited to, a processor and a memory.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component, or the like. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like, the processor is a control center of the apparatus for modeling and analyzing integrated of one kind of energy system, and various interfaces and lines are used to connect various parts of the apparatus for modeling and analyzing integrated of the whole one kind of energy system.
The memory may be used to store the computer programs and/or modules, and the processor may be used to implement various functions of the apparatus/terminal device for modeling and comprehensive analysis of the energy system by operating or executing the computer programs and/or modules stored in the memory and calling up the data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, where when the computer program runs, a device in which the computer-readable storage medium is located is controlled to execute the method for modeling and comprehensive analysis of an energy system according to any of the above embodiments. The device-integrated module for modeling and integrated analysis of an energy system may be stored in a computer-readable storage medium if it is implemented in the form of a software functional unit and sold or used as a separate product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
It should be noted that the above-described device embodiments are merely illustrative, where the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. In addition, in the drawings of the embodiment of the apparatus provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, and may be specifically implemented as one or more communication buses or signal lines. One of ordinary skill in the art can understand and implement it without inventive effort.
The invention discloses a modeling and comprehensive analysis method of an energy systemA method, an apparatus, and a storage medium. Establishing branch characteristic equations of each energy subnet in the energy transmission network according to the energy flow transmission axiom in the comprehensive energy system, and establishing branch characteristic equations of each energy subnet in the energy transmission network according to the branch characteristic equations
Figure BDA0002796718610000211
A kinetic equation of transfer and conversion is carried out, and electricity in the transfer process is established according to the branch characteristic equation
Figure BDA0002796718610000212
Press and press
Figure BDA0002796718610000213
Heat generation
Figure BDA0002796718610000214
Is/are as follows
Figure BDA0002796718610000215
A loss calculation formula; and particularly gives the electricity in the transfer process
Figure BDA0002796718610000216
Is/are as follows
Figure BDA0002796718610000217
Loss calculation formula, pressure
Figure BDA0002796718610000218
Is/are as follows
Figure BDA0002796718610000219
Damage calculation formula, heat
Figure BDA00027967186100002110
Is/are as follows
Figure BDA00027967186100002111
A loss calculation formula; according to the theoretical knowledge of heat economy, the energy resources of each strand in the network
Figure BDA00027967186100002112
Economic cost of flowAnalyzing and calculating, comprehensively considering energy conservation and economy, evaluating the energy conservation potential of the system, and adjusting the equipment parameters in the energy transmission network to enable the fexEqual to a preset threshold value N, is set,
Figure BDA00027967186100002113
the proportion of the loss cost and the non-energy cost achieves the most reasonable distribution of energy conservation and economy, reduces the waste of energy and cost, and can reduce the waste of energy and cost.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

Claims (9)

1. A method for modeling and comprehensive analysis of an energy system, comprising:
according to the transmission axiom of energy flow in an energy system, establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system;
according to the energy transmission network
Figure FDA0002796718600000018
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA0002796718600000019
Is/are as follows
Figure FDA00027967186000000110
Loss calculation formula, pressure
Figure FDA00027967186000000111
Is/are as follows
Figure FDA00027967186000000112
Damage calculation formula and heat
Figure FDA00027967186000000113
Is/are as follows
Figure FDA00027967186000000114
A loss calculation formula;
establishing an energy network equation set of the energy transmission network, and solving state quantities of all nodes in the energy transmission network according to the energy network equation;
according to the state quantities of all nodes and combining the electricity
Figure FDA00027967186000000115
Is/are as follows
Figure FDA00027967186000000116
Loss calculation formula, said pressure
Figure FDA00027967186000000117
Is/are as follows
Figure FDA00027967186000000118
Loss calculation formula and the heat
Figure FDA00027967186000000119
Is/are as follows
Figure FDA00027967186000000120
Loss calculation formula for said energy transmission network
Figure FDA00027967186000000121
And analyzing the system economy and the energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
2. The method according to claim 1, wherein the establishing of the branch characteristic equation corresponding to the energy subnet in the energy transmission network of the energy system according to the axiom of energy flow transmission in the energy system specifically comprises:
according to the axiom of the transfer of energy flow in energy systems
Figure FDA0002796718600000011
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure FDA0002796718600000012
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure FDA0002796718600000013
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA0002796718600000014
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA0002796718600000015
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; branch characteristic equation by extensive flux
Figure FDA0002796718600000016
The integral in the length direction of the transfer line in the transmission network, A being the cross-sectional area of the cylindrical transfer line through which the extensive quantity flows, JiFor extended fluence, JiObtained from the transfer axiom; in the delivery common, FiIs the driving force for pushing the extensive amount to be transmitted, KiIs a wide extension chiiThe coefficient of transmission of (a) is,
Figure FDA0002796718600000017
is the amount of strength of the conjugateAnd (4) gradient.
3. The method of claim 1, wherein the energy transmission network is based on energy from the group consisting of
Figure FDA00027967186000000212
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA00027967186000000213
Is/are as follows
Figure FDA00027967186000000214
Loss calculation formula, pressure
Figure FDA00027967186000000215
Is/are as follows
Figure FDA00027967186000000216
Damage calculation formula and heat
Figure FDA00027967186000000217
Is/are as follows
Figure FDA00027967186000000218
The loss calculation formula is specifically as follows:
according to the energy transmission network
Figure FDA00027967186000000219
Kinetic equations of transfer and conversion
Figure FDA0002796718600000021
Establishing
Figure FDA00027967186000000220
General formula of calculation
Figure FDA0002796718600000022
Wherein rho in the kinetic equation is medium density gxThe source intensity, χ, of the extensive amount χ in the unit volume of medium0A silence value that is an intensity amount x,
Figure FDA0002796718600000023
for the extended amount of transfer law, the left side of the equation is
Figure FDA00027967186000000221
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure FDA00027967186000000231
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure FDA00027967186000000223
The third term on the right represents the other forms of the transformation between
Figure FDA00027967186000000224
Converted into such a form
Figure FDA00027967186000000230
The above-mentioned
Figure FDA00027967186000000226
The loss is represented by ap in the generalized equation,
Figure FDA0002796718600000024
to increase the amount of spread in the delivery process,
Figure FDA0002796718600000025
represents an extensive flow;
according to the branch characteristic equation
Figure FDA0002796718600000026
Build up electricity
Figure FDA00027967186000000227
Is/are as follows
Figure FDA00027967186000000228
Formula for calculating loss
Figure FDA0002796718600000027
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure FDA0002796718600000028
is current, K is conductivity, RiIs the thermal resistance, l is the length of the cylindrical transfer conduit, in the electrical network, the value of voltage dead xe00, electricity
Figure FDA00027967186000000229
The loss is equal to the electric energy loss; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA0002796718600000029
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA00027967186000000210
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA00027967186000000232
Is/are as follows
Figure FDA00027967186000000233
In the formula of loss calculation
Figure FDA00027967186000000211
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
according to electricity
Figure FDA00027967186000000310
Loss calculation process build-up pressure
Figure FDA00027967186000000311
Is/are as follows
Figure FDA00027967186000000312
Formula for calculating loss
Figure FDA0002796718600000031
Wherein, in the fluid network, the silent state value of the pressure energy is zero, chipA,χpEPressure of a section A and a section E respectively, wherein RpIs the flow resistance according to the Navier-Stokes equation
Figure FDA0002796718600000032
Calculation, p is the fluid density, kpFor an extensive amount of transfer coefficient in the fluid network,
Figure FDA0002796718600000033
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
heat of formation
Figure FDA00027967186000000313
Is/are as follows
Figure FDA00027967186000000314
Formula for calculating loss
Figure FDA0002796718600000034
Wherein: the heat
Figure FDA00027967186000000315
Is/are as follows
Figure FDA00027967186000000316
X in the damage calculation formulah0The value is a silent state value of temperature, an environment temperature value is usually taken, and the temperature value at the tail end of a transmission pipeline is obtained by a Suhoff temperature drop formula:
Figure FDA0002796718600000035
λhis the heat transfer coefficient of the transfer line; in a thermal network, an incompressible fluid is used as a heat transfer medium, the heat being
Figure FDA00027967186000000317
Is/are as follows
Figure FDA00027967186000000318
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure FDA0002796718600000036
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure FDA0002796718600000037
as flow rate of heat transfer medium, χhA,χhETemperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure FDA0002796718600000038
4. the method according to claim 1, wherein the establishing an energy network equation set of the energy transmission network and solving the state quantities of all nodes in the energy transmission network according to the energy network equation set are specifically:
establishing the energy transmission network topology constraint equation set
Figure FDA0002796718600000039
Equation of the branch characteristic
Figure FDA0002796718600000041
Establishing an energy network equation set of the energy transmission network in combination with the topological constraint equation set, and solving the energy network equation set to obtain state quantities of all nodes in the energy transmission network;
wherein A in the topological constraint equation set is a correlation matrix, and BfIn the form of a matrix of elementary loops,
Figure FDA0002796718600000042
is a wide-spread flux matrix, Δ χiIs an intensity quantity difference matrix; when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA0002796718600000043
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA0002796718600000044
is volume flow, K is volume transfer coefficient rate, when the branch characteristic equation expresses a heat network equation, x is temperature,
Figure FDA0002796718600000045
flow rate of heat transfer medium, K is the heat transfer coefficient of the transfer line, l is the length of the cylindrical transfer line, KeFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows.
5. An energy system according to claim 1The method for modeling and analyzing comprehensively, characterized in that said method is based on said state quantities of all nodes and combines said electricity
Figure FDA0002796718600000047
Is/are as follows
Figure FDA0002796718600000048
Loss calculation formula, said pressure
Figure FDA0002796718600000049
Is/are as follows
Figure FDA00027967186000000410
Loss calculation formula and the heat
Figure FDA00027967186000000411
Is/are as follows
Figure FDA00027967186000000412
Loss calculation formula for said energy transmission network
Figure FDA00027967186000000413
Analyzing the system economy and energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result specifically comprises the following steps:
according to the state quantity of all nodes and the electricity
Figure FDA00027967186000000414
Is/are as follows
Figure FDA00027967186000000415
Loss calculation formula, said pressure
Figure FDA00027967186000000416
Is/are as follows
Figure FDA00027967186000000417
Loss meterFormula of calculation and the heat
Figure FDA00027967186000000418
Is/are as follows
Figure FDA00027967186000000419
Calculation of each node in network by loss calculation formula
Figure FDA00027967186000000420
Flow number, obtain unit economic cost, and unit non-energy cost
Figure FDA00027967186000000421
Flow conversion, calculating
Figure FDA00027967186000000422
The cost is reduced;
calculating technical economic coefficient
Figure FDA0002796718600000046
Wherein D isxRepresenting energy transfer and conversion processes
Figure FDA00027967186000000423
Loss value, CDxIs composed of
Figure FDA00027967186000000424
The cost is taken as input
Figure FDA00027967186000000425
Average of (2)
Figure FDA00027967186000000426
Cost, Z stands for
Figure FDA00027967186000000427
Non-energy costs of stream numerical conversion, said fractional energy costs including equipment costs, labor costs and operating costs；
Optimizing the energy transmission network according to the technical economic coefficient to make the technical economic coefficient equal to a preset threshold value N,
Figure FDA00027967186000000428
and the ratio of the loss cost to the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to the optimization result.
6. An apparatus for modeling and integrated analysis of an energy system, comprising: a branch characteristic equation calculation module,
Figure FDA0002796718600000056
The system comprises a loss calculation module, a node state calculation module and a parameter optimization module;
the branch characteristic equation calculation module is used for establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system according to the transmission axiom of the energy flow in the energy system;
the above-mentioned
Figure FDA0002796718600000057
The loss calculation module is used for calculating the loss according to the energy transmission network
Figure FDA0002796718600000058
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA0002796718600000059
Is/are as follows
Figure FDA00027967186000000510
Loss calculation formula, pressure
Figure FDA00027967186000000511
Is/are as follows
Figure FDA00027967186000000512
Damage calculation formula and heat
Figure FDA00027967186000000513
Is/are as follows
Figure FDA00027967186000000514
A loss calculation formula;
the node state calculation module is used for establishing an energy network equation set of the energy transmission network and solving all node state quantities in the energy transmission network according to the energy network equation;
the parameter optimization module is used for combining the electricity according to the state quantities of all the nodes
Figure FDA00027967186000000515
Is/are as follows
Figure FDA00027967186000000516
Loss calculation formula, said pressure
Figure FDA00027967186000000517
Is/are as follows
Figure FDA00027967186000000518
Loss calculation formula and the heat
Figure FDA00027967186000000519
Is/are as follows
Figure FDA00027967186000000520
Loss calculation formula for said energy transmission network
Figure FDA00027967186000000521
And analyzing the system economy and the energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
7. The apparatus for modeling and comprehensive analysis of an energy system according to claim 1, wherein the branch characteristic equation calculation module functions specifically as:
according to the axiom of the transfer of energy flow in energy systems
Figure FDA0002796718600000051
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure FDA0002796718600000052
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure FDA0002796718600000053
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a steady-state incompressible fluid network equation, x is pressure intensity,
Figure FDA0002796718600000054
is the volume flow, K is the volume transfer coefficient rate; branch characteristic equation by extensive flux
Figure FDA0002796718600000055
The integral in the length direction of the transfer line in the transmission network, A being the cross-sectional area of the cylindrical transfer line through which the extensive quantity flows, JiFor extended fluence, obtained from the transport axiom, FiIs the driving force for pushing the extensive amount to be transmitted, KiIs a wide extension chiiThe coefficient of transmission of (a) is,
Figure FDA0002796718600000061
is a conjugate intensity magnitude gradient;
the above-mentioned
Figure FDA00027967186000000613
The loss calculation module comprises
Figure FDA00027967186000000614
Loss of generality computing unit, electricity
Figure FDA00027967186000000615
Loss calculation unit, voltage
Figure FDA00027967186000000616
Loss calculation unit and heat
Figure FDA00027967186000000617
A loss calculation unit;
the above-mentioned
Figure FDA00027967186000000618
The loss-to-prevalence calculation unit is used for calculating the loss-to-prevalence in the energy transmission network
Figure FDA00027967186000000619
Kinetic equations of transfer and conversion
Figure FDA0002796718600000062
Establishing
Figure FDA00027967186000000620
General formula of calculation
Figure FDA0002796718600000063
Wherein rho in the kinetic equation is medium density gxThe source intensity, χ, of the extensive amount χ in the unit volume of medium0A silence value that is an intensity amount x,
Figure FDA0002796718600000064
the transfer rule is extensive quantity; the left side of the equation is
Figure FDA00027967186000000621
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure FDA00027967186000000633
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure FDA00027967186000000623
The third term on the right represents the other forms of the transformation between
Figure FDA00027967186000000624
Converted into such a form
Figure FDA00027967186000000634
The above-mentioned
Figure FDA00027967186000000626
The loss is represented by ap in the generalized equation,
Figure FDA0002796718600000065
to increase the amount of spread in the delivery process,
Figure FDA0002796718600000066
represents an extensive flow;
the electricity
Figure FDA00027967186000000627
The loss calculation unit is used for calculating the loss according to the branch characteristic equation
Figure FDA0002796718600000067
Build up electricity
Figure FDA00027967186000000628
Is/are as follows
Figure FDA00027967186000000629
Formula for calculating loss
Figure FDA0002796718600000068
Wherein, in the electric network, the voltage has a silent state value χe00, electricity
Figure FDA00027967186000000630
The loss is equal to the electric energy loss, and when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA0002796718600000069
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA00027967186000000610
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA00027967186000000611
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA00027967186000000631
Is/are as follows
Figure FDA00027967186000000632
In the formula of loss calculation
Figure FDA00027967186000000612
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
said pressure
Figure FDA00027967186000000713
Loss calculating unit for calculating loss according to electricity
Figure FDA00027967186000000714
Loss calculation process build-up pressure
Figure FDA00027967186000000715
Is/are as follows
Figure FDA00027967186000000716
Formula for calculating loss
Figure FDA0002796718600000071
Wherein, in the fluid network, the silence state value of the pressure energy is also zero, RpIs the flow resistance, chipA,χpEPressure of section A and section E, respectively, according to the Navier-Stokes equation
Figure FDA0002796718600000072
Calculation, p is the fluid density, kpFor an extensive amount of transfer coefficient in the fluid network,
Figure FDA0002796718600000073
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
the heat
Figure FDA00027967186000000717
Loss calculation unit for establishing heat
Figure FDA00027967186000000718
Is/are as follows
Figure FDA00027967186000000719
Loss calculation formula heat
Figure FDA00027967186000000720
The loss is calculated by
Figure FDA0002796718600000074
Wherein: in a thermal network, an incompressible fluid is used as a heat transfer medium, the heat being
Figure FDA00027967186000000721
Is/are as follows
Figure FDA00027967186000000722
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure FDA0002796718600000075
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure FDA0002796718600000076
as flow rate of heat transfer medium, χhA,χhETemperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure FDA0002796718600000077
the heat
Figure FDA00027967186000000723
Is/are as follows
Figure FDA00027967186000000724
X in the damage calculation formulah0The value is a silent state value of temperature, an environment temperature value is usually taken, and the temperature value at the tail end of a transmission pipeline is obtained by a Suhoff temperature drop formula:
Figure FDA0002796718600000078
λhis the heat transfer coefficient of the transfer line;
the node state calculating unit has specific functionsComprises the following steps: establishing a topological constraint equation set of the energy transmission network:
Figure FDA0002796718600000079
equation of the branch characteristic
Figure FDA00027967186000000710
Establishing an energy network equation set of the energy transmission network in combination with the topological constraint equation set, and solving the energy network equation set to obtain state quantities of all nodes of the energy transmission network; wherein A in the topological constraint equation set is a correlation matrix, and BfIn the form of a matrix of elementary loops,
Figure FDA00027967186000000711
is a wide-spread flux matrix, Δ χiIs an intensity quantity difference matrix; when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA00027967186000000712
is current, K is conductivity, RiIs thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA0002796718600000081
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA0002796718600000082
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA0002796718600000085
Is/are as follows
Figure FDA0002796718600000086
In the formula of loss calculation
Figure FDA0002796718600000083
L is the length of the cylindrical transfer line, keFor conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χeA,χeEVoltages for section a and section E, respectively;
the parameter optimization module comprises
Figure FDA0002796718600000087
The device comprises a loss cost calculation unit, a technical economic coefficient calculation unit and a parameter adjustment unit:
the above-mentioned
Figure FDA0002796718600000088
The loss cost calculation unit is used for calculating the loss cost according to the state quantity of all the nodes and the electricity
Figure FDA0002796718600000089
Is/are as follows
Figure FDA00027967186000000810
Loss calculation formula, said pressure
Figure FDA00027967186000000811
Is/are as follows
Figure FDA00027967186000000812
Loss calculation formula and the heat
Figure FDA00027967186000000813
Is/are as follows
Figure FDA00027967186000000814
Calculation of each node in network by loss calculation formula
Figure FDA00027967186000000815
Flow number, non-energy cost in units
Figure FDA00027967186000000816
Flow conversion, calculating
Figure FDA00027967186000000817
The cost is reduced;
the technical economic coefficient calculating unit is used for calculating the technical economic coefficient
Figure FDA0002796718600000084
Wherein D isxRepresenting energy transfer and conversion processes
Figure FDA00027967186000000818
Loss value, CDxIs composed of
Figure FDA00027967186000000819
The cost is taken as input
Figure FDA00027967186000000820
Average of (2)
Figure FDA00027967186000000821
Cost, Z stands for
Figure FDA00027967186000000822
Non-energy cost of stream value translation;
the parameter adjusting unit is used for optimizing the energy transmission network according to the technical economic coefficient to ensure that the technical economic coefficient is equal to a preset threshold value N,
Figure FDA00027967186000000823
and the ratio of the loss cost to the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to the optimization result.
8. An apparatus for modeling and ensemble analysis of an energy system, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing a method for modeling and ensemble analysis of an energy system according to any of claims 1-5 when executing the computer program.
9. A computer-readable storage medium, comprising a stored computer program, wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the method of modeling and ensemble analysis of an energy system according to any of claims 1-5.
CN202011334339.7A 2020-11-24 2020-11-24 Energy system modeling and comprehensive analysis method, device and storage medium Active CN112330493B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011334339.7A CN112330493B (en) 2020-11-24 2020-11-24 Energy system modeling and comprehensive analysis method, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011334339.7A CN112330493B (en) 2020-11-24 2020-11-24 Energy system modeling and comprehensive analysis method, device and storage medium

Publications (2)

Publication Number Publication Date
CN112330493A true CN112330493A (en) 2021-02-05
CN112330493B CN112330493B (en) 2022-10-28

Family

ID=74309464

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011334339.7A Active CN112330493B (en) 2020-11-24 2020-11-24 Energy system modeling and comprehensive analysis method, device and storage medium

Country Status (1)

Country Link
CN (1) CN112330493B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116759004A (en) * 2023-08-14 2023-09-15 宁德时代新能源科技股份有限公司 Model correction method, device, computer equipment and storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777708A (en) * 2016-12-21 2017-05-31 天津大学 Steady state analysis method of electric power-natural gas regional comprehensive energy system
CN108510122A (en) * 2018-03-30 2018-09-07 中国电建集团福建省电力勘测设计院有限公司 The integrated energy system optimization method of optimal double constraints is utilized based on environmental emission, hot *
CN109376428A (en) * 2018-10-24 2019-02-22 南方电网科学研究院有限责任公司 Reliability estimation method, device, equipment and the storage medium of integrated energy system
WO2019071763A1 (en) * 2017-10-09 2019-04-18 清华大学 Method for estimating state of electro-thermal coupling system in combination with dynamic characteristics of pipeline
CN109684763A (en) * 2019-01-02 2019-04-26 华南理工大学 It is a kind of based on individual be this model integrated energy system modeling method
WO2019200662A1 (en) * 2018-04-20 2019-10-24 东北大学 Stability evaluation and static control method for electricity-heat-gas integrated energy system
CN111815042A (en) * 2020-06-30 2020-10-23 天津大学 Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777708A (en) * 2016-12-21 2017-05-31 天津大学 Steady state analysis method of electric power-natural gas regional comprehensive energy system
WO2019071763A1 (en) * 2017-10-09 2019-04-18 清华大学 Method for estimating state of electro-thermal coupling system in combination with dynamic characteristics of pipeline
US20200232886A1 (en) * 2017-10-09 2020-07-23 Tsinghua University Method for estimating state of combined heat and power system
CN108510122A (en) * 2018-03-30 2018-09-07 中国电建集团福建省电力勘测设计院有限公司 The integrated energy system optimization method of optimal double constraints is utilized based on environmental emission, hot *
WO2019200662A1 (en) * 2018-04-20 2019-10-24 东北大学 Stability evaluation and static control method for electricity-heat-gas integrated energy system
CN109376428A (en) * 2018-10-24 2019-02-22 南方电网科学研究院有限责任公司 Reliability estimation method, device, equipment and the storage medium of integrated energy system
CN109684763A (en) * 2019-01-02 2019-04-26 华南理工大学 It is a kind of based on individual be this model integrated energy system modeling method
CN111815042A (en) * 2020-06-30 2020-10-23 天津大学 Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
左剑，谢平平，李银红等: "基于智能优化算法的互联电网负荷频率控制器设计及其控制性能分析", 《电工技术学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116759004A (en) * 2023-08-14 2023-09-15 宁德时代新能源科技股份有限公司 Model correction method, device, computer equipment and storage medium
CN116759004B (en) * 2023-08-14 2024-01-12 宁德时代新能源科技股份有限公司 Model correction method, device, computer equipment and storage medium

Also Published As

Publication number Publication date
CN112330493B (en) 2022-10-28

Similar Documents

Publication Publication Date Title
Qin et al. A generalized quasi-dynamic model for electric-heat coupling integrated energy system with distributed energy resources
CN110245878B (en) Distributed comprehensive energy demand response collaborative optimization method for intelligent building group
Larsen et al. Aggregated dynamic simulation model of district heating networks
Weber et al. Design and optimization of district energy systems
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN109830955B (en) Electric-gas distribution network flexible planning method considering flexible constraint and full-period cost
Gan et al. Decentralized computation method for robust operation of multi-area joint regional-district integrated energy systems with uncertain wind power
CN110991092A (en) Complementary energy integrated utilization calculation method and system
Nord et al. Data informed physical models for district heating grids with distributed heat sources to understand thermal and hydraulic aspects
Zhu et al. Integrated dispatch for combined heat and power with thermal energy storage considering heat transfer delay
CN112330493B (en) Energy system modeling and comprehensive analysis method, device and storage medium
CN108596526A (en) Demand-side comprehensive energy station synthetic operation Optimal Scheduling and method
Khalid et al. Thermal control strategies for reliable and energy-efficient data centers
CN113806972B (en) Comprehensive energy system reliability analysis method considering supply and demand bilateral flexibility
CN107748963A (en) A kind of electric network coordinated operation method based on the optimization of OCD distributings
Zhang et al. Calculation and analysis of energy storage in heat supply nets of distributed energy
CN211044253U (en) Complementary energy integrated utilization computing system
CN113886761A (en) Energy efficiency analysis and evaluation method for comprehensive energy system
Zong et al. Multi-objective optimization of multi-energy flow coupling system with carbon emission target oriented
CN111783309A (en) Dynamic simulation method of steam heating network based on internal conservation
CN116739404A (en) Construction method of networking feature module in comprehensive energy system entropy state model
CN114462163B (en) exergy hub model considering energy quality
CN114004047A (en) Electric heating gas interconnection multi-energy system matrix operation model modeling method
CN105337316B (en) A kind of distributed power generation Optimal Configuration Method based on smart allocation
Hu et al. Multi-Energy Flow Calculation Model and Solution Method for Regional Integrated Energy System

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant