CN112330493B - 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
CN112330493B
CN112330493B CN202011334339.7A CN202011334339A CN112330493B CN 112330493 B CN112330493 B CN 112330493B CN 202011334339 A CN202011334339 A CN 202011334339A CN 112330493 B CN112330493 B CN 112330493B
Authority
CN
China
Prior art keywords
energy
equation
loss
network
heat
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.)
Active
Application number
CN202011334339.7A
Other languages
Chinese (zh)
Other versions
CN112330493A (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.)
CSG Electric Power Research Institute
China Southern Power Grid Co Ltd
Original Assignee
CSG Electric Power Research Institute
China 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 CSG Electric Power Research Institute, China Southern Power Grid Co Ltd filed Critical CSG Electric Power Research Institute
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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Educational Administration (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Development Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

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 and deducing
Figure DDA0002796718630000012
While transferring in transfer lines
Figure DDA0002796718630000013
Generalized expression of loss, and establishing 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 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 conservation and the economy of the system, evaluate the energy-saving potential of the system and provide a reasonable improvement scheme for the relevant 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
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
Figure BDA00027967186100000212
Loss calculation formula, said pressure
Figure BDA00027967186100000213
Is
Figure BDA00027967186100000214
Loss calculation formula and the heat
Figure BDA00027967186100000215
Is
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:
based on the axiom of the transfer of energy flow in the energy system
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, R i Is the 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, J i For extended fluence, J i Obtained from the transfer axiom; in the transfer common, F i Is the driving force for pushing the extensive transmission, K i Is a wide spread amount x i The coefficient of transmission of (a) is,
Figure BDA0002796718610000027
is the conjugate intensity magnitude gradient.
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 g x The source intensity, χ, of the extensive amount χ in the unit volume of medium 0 A silence value that is an intensity amount x,
Figure BDA00027967186100000210
for extensive 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 flow 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
Into such a form
Figure BDA00027967186100000314
The above-mentioned
Figure BDA00027967186100000315
General damageThe calculation formula deltap represents the energy loss,
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
Building electricity
Figure BDA00027967186100000316
Is
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, R i Is the thermal resistance, l is the length of the cylindrical transfer conduit, in the electrical network, the value of voltage dead x e0 =0, 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
Figure BDA00027967186100000320
In the formula of loss calculation
Figure BDA0002796718610000038
L is the length of the cylindrical transfer line, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows, χ eA ,χ eE Voltages 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
Wherein, in the fluid network, the silent state value of the pressure energy is zero, chi pA ,χ pE Pressure of a section A and a section E respectively, wherein R p Is the flow resistance according to the Navier-Stokes equation
Figure BDA00027967186100000310
Calculation of ρ is fluid density, k p For 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 formula h0 The 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
λ h is 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
Figure BDA00027967186100000412
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure BDA0002796718610000042
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure BDA0002796718610000043
is the flow rate of the heat transfer medium,% hA ,χ hE Temperatures 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 B f Is a matrix of the basic loop,
Figure BDA0002796718610000047
is a wide-spread flux matrix, Δ χ i Is 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, R i Is 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 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, K e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows.
Preferably, said combining said electricity according to said all-node state quantities
Figure BDA00027967186100000413
Is
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
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 value, obtain unit economic cost, and unit non-energy cost
Figure BDA0002796718610000059
Flow conversion, calculation
Figure BDA00027967186100000510
The cost is reduced;
calculating the technical economic coefficient
Figure BDA0002796718610000051
Wherein D is x Representing energy transfer and conversion processes
Figure BDA00027967186100000511
Loss value, C Dx Is composed of
Figure BDA00027967186100000512
The cost is taken as input
Figure BDA00027967186100000513
Is averaged
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 and economic coefficients so that the technical and economic coefficients are equal to a preset threshold value N,
Figure BDA00027967186100000516
and the proportion 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 established, 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
Energy transfer in loss calculation equationThe influence of the parameters in the network on the technical economic coefficient, and the equipment parameters in the energy transmission network are adjusted to enable the f ex Equal to a preset threshold value N, is set,
Figure BDA00027967186100000527
the proportion of the loss cost to the non-energy cost achieves the most reasonable distribution of energy conservation and economy, and reduces the waste of energy and cost.
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
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
Figure BDA00027967186100000613
Loss calculation formula, said pressure
Figure BDA00027967186100000614
Is/are as follows
Figure BDA00027967186100000615
Loss calculation formula and the heat
Figure BDA00027967186100000616
Is
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:
based on the axiom of the transfer of energy flow in the energy system
Figure BDA0002796718610000061
Set up theBranch characteristic equation corresponding to energy subnet in energy transmission network of energy system
Figure BDA0002796718610000062
Wherein when the branch characteristic equation represents an electric network equation, x is voltage,
Figure BDA0002796718610000063
is current, K is conductivity, R i Is 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, J i For extended fluence, obtained from the transport axiom, F i Is the driving force for pushing the extensive amount to be transmitted, K i Is a wide spread amount x i The coefficient of transmission of (a) is,
Figure BDA0002796718610000066
is the 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 andheat generation
Figure BDA00027967186100000623
A loss calculation unit;
the above-mentioned
Figure BDA00027967186100000624
The loss common calculating unit is used for calculating the loss common in the energy transmission network according to the energy
Figure BDA00027967186100000625
Kinetic equations of transfer and conversion
Figure BDA0002796718610000067
Establishing
Figure BDA00027967186100000626
Loss common formula
Figure BDA0002796718610000068
Wherein rho in the kinetic equation is medium density g x The source intensity, χ, of the extensive amount χ in the unit volume of medium 0 A 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 described
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
Building 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 χ e0 =0, 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, R i Is the thermal resistance, l is the length of the cylindrical transfer line; the branch characteristicsWhen the equation of properties represents a fluid network equation, χ is the pressure,
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, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χ eA ,χ eE Voltages for section a and section E, respectively;
said pressure
Figure BDA00027967186100000724
Loss calculating unit for calculating loss according to electricity
Figure BDA00027967186100000725
Loss calculation process build-up pressure
Figure BDA00027967186100000726
Is
Figure BDA00027967186100000727
Formula for calculating loss
Figure BDA0002796718610000079
Wherein, in the fluid network, the silence state value of the pressure energy is also zero, R p Is the flow resistance, chi pA ,χ pE Pressure of the cross section A and the cross section E respectively, whereinVieie-stokes equation
Figure BDA00027967186100000710
Calculation of ρ is fluid density, k p For 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 a thermal network, an incompressible fluid is used as a heat transfer medium, the heat being
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
for transferring heatFlow rate of medium, chi hA ,χ hE Temperatures 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 formula h0 The 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
λ h is 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
Establishing an energy network equation set of the energy transmission network in combination with the topological constraint equation set, 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 B f Is a matrix of the basic loop,
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, R i Is the thermal resistance, l is the length of the cylindrical transfer line; said branch characteristic equation tableWhen showing the fluid network equation, χ is the pressure,
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
Figure BDA00027967186100000814
In the formula of loss calculation
Figure BDA0002796718610000086
L is the length of the cylindrical transfer line, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows, χ eA ,χ eE Voltages for section a and section E, respectively;
the parameter optimization module comprises
Figure BDA00027967186100000815
The device comprises a loss cost calculation unit, a technical economic coefficient calculation unit and a 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
Figure BDA00027967186100000822
Calculating 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 is x Representing energy transfer and conversion processes
Figure BDA0002796718610000092
Loss value, C Dx Is composed of
Figure BDA0002796718610000093
The cost is taken as input
Figure BDA0002796718610000094
Is averaged
Figure BDA0002796718610000095
Cost, Z stands for
Figure BDA0002796718610000096
Non-energy cost of stream value translation;
for the parameter adjustment unitOptimizing the energy transmission network according to the technical-economic coefficient so that the technical-economic coefficient is equal to a preset threshold value N,
Figure BDA0002796718610000097
and the proportion 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 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 and can also establish a branch characteristic equation of each energy subnet in the energy transmission network according to the transmission axiom of energy flow 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
The influence of parameters in the energy transmission network in the loss calculation equation on the technical economic coefficient is used for adjusting the equipment parameters in the energy transmission network so as to enable the f ex Equal to a preset threshold value N, is set,
Figure BDA00027967186100000913
the proportion of the loss cost to the non-energy cost achieves the most reasonable distribution of energy conservation and economy, and reduces the waste of energy and cost.
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, according to all the node shapesIn combination with said electricity
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 energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result.
According to the modeling and comprehensive analysis method of the energy system, provided by the embodiment of the invention, a branch characteristic equation in each energy flow subnet in an energy transmission network is established based on the transmission characteristics of different energy types according to the transmission axiom of the energy flow transmission rule in the comprehensive energy system; 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 and the economy of the system, evaluate the energy conservation potential of the system, and propose an improved scheme for related 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 transferring energy flows, the transferring 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, J i For extended fluence, F i Is the driving force for pushing the extensive amount to be transmitted, K i Is a wide extension chi i The 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 expression 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, under the pushing of corresponding intensity quantity difference, the flow transportation of extensive quantity is realized, so that various forms of energy flow transfer are completed.
The energy transfer is represented by the formulaIt can be seen that the propagation law of the spread is determined by two factors, namely, the conjugate intensity gradient, which is the root cause for the generation of the basic spread, and the propagation coefficient K i It 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., J iA Flow direction J iE Combined with transmission line pipe, can further derive extensive flow
Figure BDA0002796718610000125
The expression of (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 K i May or may not be changed, and thus will be discussed in two cases.
Referring to fig. 3, it is a schematic diagram of a physical model of energy transfer when the coefficient of transfer of the extensive quantity is constant, where the coefficient of transfer K of the extensive quantity is constant i For 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 is i Is 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 line iA -χ iE And a constant value R i And (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 constant i When not of constant value, but of coefficient of transfer of energy λ i ＝K i χ i For a fixed value, 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 is i For 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; p h Is the heat flow, i.e. the thermal power, R h Representing the thermal resistance of the radial heat conduction. If it is usedConsidering the heat convection process between the fluid inside and outside the pipe and the pipe wall, 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, λ hd Is the coefficient of thermal conductivity; lambda [ alpha ] hc1 Is the convective heat transfer coefficient of the inner wall, λ hc2 Is the convective heat transfer coefficient of the outer wall.
Through the analysis, the extensive equivalent transfer equations of the three energies of the electric energy, the pressure energy and the heat energy in the steady-state laminar flow state in the transfer pipeline have similar expression forms, namely the extensive 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 of transfer and conversion are:
Figure BDA0002796718610000142
where ρ is the density of the medium, g x Is the source intensity of the etendue χ in a unit volume of medium,
Figure BDA0002796718610000143
is a broad transfer law.
The left side of the equation is
Figure BDA00027967186100001413
Rate of change over time, right firstThe term representation flowing in through the voxel boundary
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
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
Damage calculation formula, heat
Figure BDA00027967186100001422
Is
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, χ 0 Is a silence value of the intensity amount χ.
In an electrical network, the value of voltage dead x e0 =0, so electricity
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
k e Represents 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, χ h0 The temperature value is a temperature silent value, an environment temperature value is usually taken, and the temperature value at the tail end of the transfer pipeline can be obtained by a thref temperature drop formula:
Figure BDA0002796718610000156
wherein λ is h Is the heat transfer coefficient of the transfer line.
In the incompressible steady-state laminar flow network, the reference power is zero because the silence 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 is p Is the flow resistance, which can be obtained according to the Navier-Stokes equation
Figure BDA0002796718610000161
In another preferred embodiment, step S103 is specifically:
in order to describe the inherent relation between the node intensity quantity 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, B f In the form of a matrix of elementary loops,
Figure BDA00027967186100001620
is a wide-spread flux matrix, Δ χ i Is 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
Pressing and pressing
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
Conversion of flow, to calculate the quantity of electricity
Figure BDA00027967186100001610
Cost and cold capacity
Figure BDA00027967186100001611
Cost;
by calculating the technical economic coefficient f ex ：
Figure BDA0002796718610000163
D x Representing energy transfer and conversion processes
Figure BDA00027967186100001612
Loss value, C Dx Is composed of
Figure BDA00027967186100001613
The cost is taken as input
Figure BDA00027967186100001614
Is averaged
Figure BDA00027967186100001615
Cost, Z stands for
Figure BDA00027967186100001616
Non-energy cost of stream value translation;
f is ex Can 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 coefficient ex Equal 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, B f Is a matrix of the basic loop,
Figure BDA00027967186100001717
is a wide-spread flux matrix, Δ χ i Is 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:
c gas E x,gas +Z HP ＝c e E xe +c h E xh
c h E xh +Z c ＝c c E xc
wherein c represents
Figure BDA0002796718610000174
Cost per unit of heat economy of flow, E x Represent
Figure BDA0002796718610000175
Flow number, Z HP Non-energy costs (equipment costs and labor costs of gas turbines, etc.) representing the cogeneration link, Z c Representing the non-energy costs (equipment costs and operating costs, etc.) of the waste heat recovery refrigeration equipment. Subscript g as E, 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 in units
Figure BDA00027967186100001710
Conversion of flow, to calculate the quantity of electricity
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:
c e,grid E x,grid +Z ec ＝c c2 E xc2
wherein Z is ec Representing 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.
Power grid purchase cost c e,grid Can be priced by local electricity and generate electricity from natural gas
Figure BDA00027967186100001716
Per unit economic cost c e The cost c of purchasing electricity from the power grid e,grid By contrast, the two modes of electricity production can be compared
Figure BDA0002796718610000182
The economy of the process.
By a technical economic factor f ex Can reflect
Figure BDA0002796718610000183
The proportional relationship between the loss cost and the non-energy cost is shown as the following formula:
Figure BDA0002796718610000181
in the formula, D x Representing processes of energy transfer, conversion
Figure BDA0002796718610000184
The value of the loss is reduced,
Figure BDA00027967186100001821
is composed of
Figure BDA0002796718610000185
To lose cost, take as input
Figure BDA0002796718610000186
Is averaged
Figure BDA0002796718610000187
Cost, Z stands for
Figure BDA0002796718610000188
Flow value converted non-energy cost. f. of ex When the temperature is too high, the temperature is high,
Figure BDA0002796718610000189
the loss cost is low, but the energy cost investment is too large, so that the economy is not high enough; f. of ex If 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 is ex When the pressure is not greater than 0.5%,
Figure BDA00027967186100001811
the loss and non-energy costs reach 1, which is considered the most reasonable distribution, and by this parameter, equipment replacement and capital allocation adjustments in the network can be considered.
When the waste heat recovery refrigeration equipment is used for cooling, cold water pipelines are laid at cold load positions to carry out 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 f ex (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
A kinetic equation of transfer and conversion is established 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 of the strands 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 waste of energy and cost is reduced.
Referring to fig. 6, it is a schematic diagram of a device for modeling and comprehensive analysis of an energy system according to an embodiment of the present invention, and as shown in the drawing, an energy system according to an embodiment of the present inventionThe device for modeling and comprehensive analysis of the source system comprises 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 apparatus implements a method for modeling and comprehensive analysis of an energy system as described in any of the above embodiments.
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 devices. 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 comprehensive 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 an independent 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 described above 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 embodiments of the apparatus are merely illustrative, where the units described as separate parts may or may not be physically separate, and 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 method, a device and a storage medium for modeling and comprehensive analysis of an energy system. 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 theory knowledge of heat economy, each strand in the energy network is supplied with energy
Figure BDA00027967186100002112
Analyzing and calculating the economic cost of the flow, comprehensively considering the energy saving performance and the economic performance, evaluating the energy saving potential of the system, and adjusting the equipment parameters in the energy transmission network to enable the f ex Equal to a preset threshold value N, is set,
Figure BDA00027967186100002113
the proportion of the loss cost to 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 the energy and the 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 (8)

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 FDA0003750007150000017
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA0003750007150000018
Is
Figure FDA0003750007150000019
Loss calculation formula, pressure
Figure FDA00037500071500000110
Is/are as follows
Figure FDA00037500071500000111
Damage calculation formula and heat
Figure FDA00037500071500000112
Is
Figure FDA00037500071500000113
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 FDA00037500071500000114
Is
Figure FDA00037500071500000115
Loss calculation formula, said pressure
Figure FDA00037500071500000116
Is/are as follows
Figure FDA00037500071500000117
Loss calculation formula and the heat
Figure FDA00037500071500000118
Is/are as follows
Figure FDA00037500071500000119
Loss calculation formula for said energy transmission network
Figure FDA00037500071500000120
Analyzing the system economy and energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result;
according to the energy transmission network
Figure FDA00037500071500000121
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA00037500071500000122
Is/are as follows
Figure FDA00037500071500000123
Loss calculation formula, pressure
Figure FDA00037500071500000124
Is/are as follows
Figure FDA00037500071500000125
Damage calculation formula and heat
Figure FDA00037500071500000126
Is/are as follows
Figure FDA00037500071500000127
The loss calculation formula is specifically as follows:
according to the energy transmission network
Figure FDA00037500071500000128
Equation of transfer and conversion kinetics
Figure FDA0003750007150000011
Establishing
Figure FDA00037500071500000129
General formula of calculation
Figure FDA0003750007150000012
Wherein rho in the kinetic equation is medium density g x Source intensity, χ, which is the extensive amount χ in a unit volume of medium 0 A silence value that is an intensity amount x,
Figure FDA0003750007150000013
the left side of the equation of the transfer and conversion dynamics is 15794
Figure FDA00037500071500000130
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure FDA00037500071500000138
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure FDA00037500071500000132
The third term on the right represents the other forms of the transformation between
Figure FDA00037500071500000133
Conversion to be calculated
Figure FDA00037500071500000139
The above-mentioned
Figure FDA00037500071500000135
The loss is represented by ap in the generalized equation,
Figure FDA0003750007150000014
to increase the amount of spread in the delivery process,
Figure FDA0003750007150000015
represents an extensive flow;
according to the branch characteristic equation
Figure FDA0003750007150000016
Build up electricity
Figure FDA00037500071500000136
Is/are as follows
Figure FDA00037500071500000137
Formula for calculating loss
Figure FDA0003750007150000021
Wherein, when the branch characteristic equation represents the electric network equation, χ is voltage,
Figure FDA0003750007150000022
is current, K is conductivity, R i For thermal resistance, l is the length of the cylindrical transfer line, in an electrical network, of voltageValue of silence χ e0 =0, electricity
Figure FDA00037500071500000213
The loss is equal to the electric energy loss; when the branch characteristic equation represents a fluid network equation, x is pressure intensity,
Figure FDA0003750007150000023
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA0003750007150000024
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA00037500071500000214
Is/are as follows
Figure FDA00037500071500000215
In the formula of loss calculation
Figure FDA0003750007150000025
L is the length of the cylindrical transfer line, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χ eA ,χ eE Voltages for section a and section E, respectively;
according to electricity
Figure FDA00037500071500000216
Loss calculation process build-up pressure
Figure FDA00037500071500000217
Is/are as follows
Figure FDA00037500071500000218
Formula for calculating loss
Figure FDA0003750007150000026
Wherein, in the fluid network, the silent state value of the pressure energy is zero, chi pA ,χ pE Pressure of a section A and a section E respectively, wherein R p Is flow resistance according to the Navier-Stokes equation
Figure FDA0003750007150000027
Calculation of ρ is fluid density, k p For an extensive amount of transfer coefficient in the fluid network,
Figure FDA0003750007150000028
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
heat of formation
Figure FDA00037500071500000219
Is/are as follows
Figure FDA00037500071500000220
Formula for calculating loss
Figure FDA0003750007150000029
Wherein: the heat
Figure FDA00037500071500000221
Is/are as follows
Figure FDA00037500071500000222
X in the damage calculation formula h0 The 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 FDA00037500071500000210
λ h is 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 FDA00037500071500000223
Is
Figure FDA00037500071500000224
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure FDA00037500071500000211
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure FDA00037500071500000212
as flow rate of heat transfer medium, χ hA ,χ hE Temperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure FDA0003750007150000031
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:
based on the axiom of the transfer of energy flow in the energy system
Figure FDA0003750007150000032
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure FDA0003750007150000033
Wherein when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA0003750007150000034
is current, K is conductivity, R i Is thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation expresses a fluid network equation, x is the pressure intensity,
Figure FDA0003750007150000035
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA0003750007150000036
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 FDA0003750007150000037
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, J i For extended fluence, J i Obtained from the transfer axiom; in the delivery common, F i Is the driving force for pushing the extensive amount to be transmitted, K i Is a wide extension chi i The coefficient of transmission of (a) is,
Figure FDA0003750007150000038
is the strength magnitude gradient of the conjugate.
3. The method according to claim 1, wherein the establishing of the energy network equation set of the energy transmission network and the solving of the state quantities of all the 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 FDA0003750007150000039
Equation of the branch characteristic
Figure FDA00037500071500000310
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 the energy in the energy transmission networkA node state quantity;
wherein, A in the topological constraint equation system is a correlation matrix, B f In the form of a matrix of elementary loops,
Figure FDA0003750007150000041
is a wide-spread flux matrix, Δ χ i Is an intensity quantity difference matrix; when the branch characteristic equation expresses the electrical network equation, x is voltage,
Figure FDA0003750007150000042
is current, K is conductivity, R i Is 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 FDA0003750007150000043
is volume flow, K is volume transfer coefficient rate, when the branch characteristic equation expresses a heat network equation, x is temperature,
Figure FDA0003750007150000044
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, K e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows.
4. The method of claim 1, wherein said method of modeling and analyzing comprises combining said electricity with said state quantities of all nodes
Figure FDA0003750007150000046
Is/are as follows
Figure FDA0003750007150000047
Loss calculation formula, said pressure
Figure FDA0003750007150000048
Is/are as follows
Figure FDA0003750007150000049
Loss calculation formula and the heat
Figure FDA00037500071500000410
Is/are as follows
Figure FDA00037500071500000411
Loss calculation formula for said energy transmission network
Figure FDA00037500071500000412
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 as follows:
according to the state quantity of all nodes and the electricity
Figure FDA00037500071500000413
Is/are as follows
Figure FDA00037500071500000414
Loss calculation formula, said pressure
Figure FDA00037500071500000415
Is/are as follows
Figure FDA00037500071500000416
Loss calculation formula and the heat
Figure FDA00037500071500000417
Is/are as follows
Figure FDA00037500071500000418
Calculating each node in network by loss calculation formula
Figure FDA00037500071500000419
Flow number, obtain unit economic cost, and unit non-energy cost
Figure FDA00037500071500000420
Flow conversion, calculating
Figure FDA00037500071500000421
The cost is reduced;
calculating the technical economic coefficient
Figure FDA0003750007150000045
Wherein D is x Representing energy transfer and conversion processes
Figure FDA00037500071500000422
Loss value, C Dx Is composed of
Figure FDA00037500071500000423
The cost is taken as input
Figure FDA00037500071500000424
Is averaged
Figure FDA00037500071500000425
Cost, Z stands for
Figure FDA00037500071500000426
A flow numerical reduced non-energy cost, the non-energy cost including 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 FDA00037500071500000427
the proportion of the loss cost and the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to an optimization result.
5. Modeling and comprehensive analysis of energy systemThe apparatus of (a), comprising: a branch characteristic equation calculation module,
Figure FDA0003750007150000057
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 described
Figure FDA0003750007150000058
The loss calculation module is used for calculating the loss according to the energy transmission network
Figure FDA0003750007150000059
Transfer and conversion kinetic equation and branch characteristic equation establishing electricity in transfer process
Figure FDA00037500071500000510
Is
Figure FDA00037500071500000511
Loss calculation formula, pressure
Figure FDA00037500071500000512
Is/are as follows
Figure FDA00037500071500000513
Damage calculation formula and heat
Figure FDA00037500071500000514
Is/are as follows
Figure FDA00037500071500000515
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 FDA00037500071500000516
Is/are as follows
Figure FDA00037500071500000517
Loss calculation formula, said pressure
Figure FDA00037500071500000518
Is
Figure FDA00037500071500000519
Loss calculation formula and the heat
Figure FDA00037500071500000520
Is
Figure FDA00037500071500000521
Loss calculation formula for said energy transmission network
Figure FDA00037500071500000522
Analyzing the system economy and energy conservation of the flow, and adjusting the parameters of the energy transmission network according to the analysis result;
the described
Figure FDA00037500071500000523
The loss calculation module comprises
Figure FDA00037500071500000524
Loss of generality computing unit, electricity
Figure FDA00037500071500000525
Loss calculation unit, voltage
Figure FDA00037500071500000526
Damage calculation unit and heat
Figure FDA00037500071500000527
A loss calculation unit;
the above-mentioned
Figure FDA00037500071500000528
The loss common calculating unit is used for calculating the loss common in the energy transmission network according to the energy
Figure FDA00037500071500000529
Equation of transfer and conversion kinetics
Figure FDA0003750007150000051
Establishing
Figure FDA00037500071500000530
General formula of calculation
Figure FDA0003750007150000052
Wherein rho in the kinetic equation is medium density g x The source intensity, χ, of the extensive amount χ in the unit volume of medium 0 A silence value that is an intensity amount x,
Figure FDA0003750007150000053
the transfer rule is extensive quantity; the left side of the equation of transmission and conversion kinetics is 15794
Figure FDA00037500071500000531
Rate of change over time, the first term on the right representing inflow through voxel boundaries
Figure FDA00037500071500000538
The second term on the right represents the other forms driven by the intensity magnitude gradient
Figure FDA00037500071500000533
The third term on the right represents the other forms of the transformation between
Figure FDA00037500071500000534
Conversion to be calculated
Figure FDA00037500071500000539
The above-mentioned
Figure FDA00037500071500000536
Loss in the general calculation formula deltap represents the energy loss,
Figure FDA0003750007150000054
to increase the amount of spread in the delivery process,
Figure FDA0003750007150000055
represents an extensive flow;
the electricity
Figure FDA00037500071500000537
The loss calculation unit is used for calculating the loss according to the branch characteristic equation
Figure FDA0003750007150000056
Build up electricity
Figure FDA00037500071500000613
Is/are as follows
Figure FDA00037500071500000614
Formula for calculating loss
Figure FDA0003750007150000061
Wherein, in the electric network, the voltage has a silent state value χ e0 =0, electricity
Figure FDA00037500071500000615
The loss is equal to the electric energy loss, and when the branch characteristic equation expresses an electric network equation, x is voltage,
Figure FDA0003750007150000062
is current, K is conductivity, R i Is thermal resistance, l is the length of the cylindrical transfer line; when the branch characteristic equation expresses a fluid network equation, x is the pressure intensity,
Figure FDA0003750007150000063
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA0003750007150000064
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA00037500071500000616
Is
Figure FDA00037500071500000617
In the formula of loss calculation
Figure FDA0003750007150000065
L is the length of the cylindrical transfer line, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the extensive volume flows, χ eA ,χ eE Voltages for section a and section E, respectively;
said pressure
Figure FDA00037500071500000618
Loss calculating unit for calculating loss according to electricity
Figure FDA00037500071500000619
Loss calculation process build-up pressure
Figure FDA00037500071500000620
Is/are as follows
Figure FDA00037500071500000621
Formula for calculating loss
Figure FDA0003750007150000066
Wherein, in the fluid network, the silence state value of the pressure energy is also zero, R p Is the flow resistance, chi pA ,χ pE Pressure of section A and section E, respectively, according to the Navier-Stokes equation
Figure FDA0003750007150000067
Calculation of ρ is fluid density, k p For an extensive amount of transfer coefficient in the fluid network,
Figure FDA0003750007150000068
is the volume flow, f is the fluid friction coefficient; d is the diameter of the transfer pipe;
said heat
Figure FDA00037500071500000622
Loss calculation unit for establishing heat
Figure FDA00037500071500000623
Is
Figure FDA00037500071500000624
Loss calculation formula heat
Figure FDA00037500071500000625
The loss is calculated by
Figure FDA0003750007150000069
Wherein: in a thermal network, an incompressible fluid is used as a heat transfer medium, the heat being
Figure FDA00037500071500000626
Is/are as follows
Figure FDA00037500071500000627
The calculation formula of the loss is a combination of an entropy increase calculation formula and a heat energy loss calculation formula
Figure FDA00037500071500000610
Rho and c are the density and specific heat capacity of the heat transfer medium in turn,
Figure FDA00037500071500000611
is the flow rate of the heat transfer medium,% hA ,χ hE Temperatures for section a and section E, respectively; the heat energy loss calculation formula is as follows:
Figure FDA00037500071500000612
the heat
Figure FDA00037500071500000628
Is
Figure FDA00037500071500000629
Loss calculation formula χ h0 The value is a silent value of temperature, an environmental temperature value is usually taken, and the temperature value at the tail end of a transfer pipeline is obtained by a thref temperature drop formula:
Figure FDA0003750007150000071
λ h is the heat transfer coefficient of the transfer line.
6. The apparatus for modeling and comprehensive analysis of an energy system according to claim 5, wherein the branch characteristic equation calculation module functions specifically as:
based on the axiom of the transfer of energy flow in the energy system
Figure FDA0003750007150000072
Establishing a branch characteristic equation corresponding to an energy subnet in an energy transmission network of the energy system
Figure FDA0003750007150000073
Wherein when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA0003750007150000074
is current, K is conductivity, R i Is 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 FDA0003750007150000075
is the volume flow, K is the volume transfer coefficient rate; branch characteristic equation by extensive flux
Figure FDA0003750007150000076
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, J i For extended fluence, derived from the transport axiom, F i Is the driving force for pushing the extensive amount to be transmitted, K i Is a wide extension chi i The coefficient of transmission of (a) is,
Figure FDA0003750007150000077
is the conjugate intensity magnitude gradient;
the node state calculating unit has the functions of: establishing a topological constraint equation set of the energy transmission network:
Figure FDA0003750007150000078
equation of said branch characteristic
Figure FDA0003750007150000079
Establishing the energy of the energy transmission network in combination with the topological constraint equation setThe network equation set is solved, and the state quantities of all the nodes of the energy transmission network are obtained by solving the energy network equation set; wherein A in the topological constraint equation set is a correlation matrix, and B f In the form of a matrix of elementary loops,
Figure FDA00037500071500000710
is a wide-spread flux matrix, Δ χ i Is an intensity quantity difference matrix; when the branch characteristic equation represents an electric network equation, x is voltage,
Figure FDA00037500071500000711
is current, K is conductivity, R i Is 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 FDA00037500071500000712
is volume flow, K is volume transfer coefficient rate, chi is temperature when the branch characteristic equation expresses a heat network equation,
Figure FDA00037500071500000713
the flow rate of the heat transfer medium, K is the heat transfer coefficient of the transfer line; the electricity
Figure FDA00037500071500000714
Is/are as follows
Figure FDA00037500071500000715
In the formula of loss calculation
Figure FDA0003750007150000081
L is the length of the cylindrical transfer line, k e For conductivity, A is the cross-sectional area of the cylindrical transfer line through which the elongation flows, χ eA ,χ eE Voltages for section a and section E, respectively;
the parameter optimization module comprises
Figure FDA0003750007150000083
The device comprises a loss cost calculation unit, a technical economic coefficient calculation unit and a parameter adjustment unit:
the above-mentioned
Figure FDA0003750007150000084
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 FDA0003750007150000085
Is
Figure FDA0003750007150000086
Loss calculation formula, said pressure
Figure FDA0003750007150000087
Is/are as follows
Figure FDA0003750007150000088
Loss calculation formula and the heat
Figure FDA0003750007150000089
Is/are as follows
Figure FDA00037500071500000810
Calculating each node in network by loss calculation formula
Figure FDA00037500071500000811
Flow number, non-energy cost in units
Figure FDA00037500071500000812
Flow conversion, calculating
Figure FDA00037500071500000813
The cost is reduced;
the technical economic coefficient calculating unit is used for calculating the technical economic coefficient
Figure FDA0003750007150000082
Wherein D is x Representing energy transfer and conversion processes
Figure FDA00037500071500000814
Loss value, C Dx Is composed of
Figure FDA00037500071500000815
The cost is taken as input
Figure FDA00037500071500000816
Average of (2)
Figure FDA00037500071500000817
Cost, Z stands for
Figure FDA00037500071500000818
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 FDA00037500071500000819
the proportion of the loss cost and the non-energy cost achieves reasonable distribution of energy conservation and economy, and parameter setting in the energy transmission network is adjusted according to an optimization result.
7. An apparatus for modeling and analysis-by-synthesis 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 when executing the computer program implementing a method for modeling and analysis-by-synthesis of an energy system according to any of claims 1-4.
8. 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-4.
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 CN112330493A (en) 2021-02-05
CN112330493B true 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)

Families Citing this family (1)

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

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777708B (en) * 2016-12-21 2021-02-05 天津大学 Steady state analysis method of electric power-natural gas regional comprehensive energy system
CN107563674B (en) * 2017-10-09 2020-07-10 清华大学 Electric-thermal coupling system state estimation method considering pipeline dynamic characteristics
CN108510122B (en) * 2018-03-30 2021-08-31 中国电建集团福建省电力勘测设计院有限公司 Comprehensive energy system optimization method based on optimal double-constraint condition
CN108494021B (en) * 2018-04-20 2021-06-01 东北大学 Stability evaluation and static control method of electricity-heat-gas comprehensive energy system
CN109376428B (en) * 2018-10-24 2022-10-14 南方电网科学研究院有限责任公司 Reliability evaluation method, device, equipment and storage medium of comprehensive 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
CN111815042B (en) * 2020-06-30 2022-02-18 天津大学 Electric heating comprehensive energy system optimization scheduling method considering refinement heat supply network model

Also Published As

Publication number Publication date
CN112330493A (en) 2021-02-05

Similar Documents

Publication Publication Date Title
CN110245878B (en) Distributed comprehensive energy demand response collaborative optimization method for intelligent building group
Qin et al. A generalized quasi-dynamic model for electric-heat coupling integrated energy system with distributed energy resources
Larsen et al. Aggregated dynamic simulation model of district heating networks
CN111061986B (en) Thermoelectric comprehensive energy system tide calculation method with multiple operation modes
CN109830955B (en) Electric-gas distribution network flexible planning method considering flexible constraint and full-period cost
CN109102104B (en) Multi-energy system collaborative low-carbon planning method based on carbon emission flow
CN109492325B (en) Flow analysis method of multi-energy coupling system based on extended energy concentrator
CN112330493B (en) Energy system modeling and comprehensive analysis method, device and storage medium
Nord et al. Data informed physical models for district heating grids with distributed heat sources to understand thermal and hydraulic aspects
CN110991092A (en) Complementary energy integrated utilization calculation method and system
CN111046594A (en) Hot water heating network dynamic simulation method based on cross iteration principle
Zhu et al. Integrated dispatch for combined heat and power with thermal energy storage considering heat transfer delay
CN108596526A (en) Demand-side comprehensive energy station synthetic operation Optimal Scheduling and method
CN115062555A (en) Integrated energy system exergy flow direct calculation method based on unbalanced node exergy
CN113806972B (en) Comprehensive energy system reliability analysis method considering supply and demand bilateral flexibility
Tian et al. Modeling and simulation for multi energy flow coupled network computing
CN211044253U (en) Complementary energy integrated utilization computing system
CN116739404B (en) Construction method of networking feature module in comprehensive energy system entropy state model
CN111783309A (en) Dynamic simulation method of steam heating network based on internal conservation
CN114004047B (en) Modeling method for matrixing operation model of electric heating gas interconnection multi-energy system
CN111062598A (en) Distributed optimal scheduling method and system for comprehensive energy system
CN114462163B (en) exergy hub model considering energy quality
CN112182905B (en) Heat supply pipe network simulation method and device for comprehensive energy system
CN115481858A (en) Electric-gas system optimal power flow dispersion calculation method and system considering demand response
CN115234973A (en) Method for determining dynamic heat supply strategy of heat supply network

Legal Events

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