CN112182905B - Heat supply pipe network simulation method and device for comprehensive energy system - Google Patents
Heat supply pipe network simulation method and device for comprehensive energy system Download PDFInfo
- Publication number
- CN112182905B CN112182905B CN202011108790.7A CN202011108790A CN112182905B CN 112182905 B CN112182905 B CN 112182905B CN 202011108790 A CN202011108790 A CN 202011108790A CN 112182905 B CN112182905 B CN 112182905B
- Authority
- CN
- China
- Prior art keywords
- pipe network
- supply pipe
- heat supply
- module
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The application discloses a heating network simulation method and a device for a comprehensive energy system, which construct a heating network module, wherein the heating network module comprises: input quantity, output quantity and characteristic parameters; constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module; and according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the heat supply pipe network is simulated by considering the change of load nodes and/or heat source nodes, and the hydraulic-thermal balance point of the heat supply pipe network is obtained. The application realizes the development of the heat supply pipe network module based on the function of the TRNSYS rapid development module, takes external factors as input variables of a load model, fully considers the influence of external factor changes on load nodes and the influence of comprehensive energy system operation mode changes on heat supply sources, considers the load nodes and/or heat source node changes to simulate the heat supply pipe network, and improves the accuracy of comprehensive energy simulation.
Description
Technical Field
The application relates to the technical field of comprehensive energy system simulation, in particular to a heating pipe network simulation method and device for a comprehensive energy system.
Background
The comprehensive energy system organically integrates energy links such as electric power, fuel gas, heat supply/cold supply, hydrogen supply and the like with the supporting system such as traffic, information and the like, so that the cooperative coupling among different energy sources is realized, the efficient utilization of energy sources is realized, the multiple energy source demands of users are met, and the reliability and the safety of social energy supply are improved. The electric heating combined system is a typical comprehensive energy multi-energy flow network system, and the tide calculation is the basis of other various analysis calculations.
The typical integrated energy system consists of an energy supply system, an energy network and a load. The energy network mainly comprises a power grid and a heat supply pipe network, and is used for connecting a user with the energy supply system. The heat supply pipe network has complex network structure, delay in energy transmission and close correlation between the change condition of the hydraulic and thermal balance points of the heat supply pipe network and the load change. Compared with the traditional energy supply mode, the comprehensive energy system has the advantages that certain elements (such as a combined cooling, heating, power supplying system) in the comprehensive energy system can supply heat and power, and are generally called as coupling elements. In order to obtain more accurate results in the simulation process of the comprehensive energy system, the solving method of the change of the hydraulic and thermal balance points of the heating network along with the load and the heat source is a key technical problem of the simulation of the heating network of the comprehensive energy system.
The heat supply network consists of a heat source, a heat supply network and heat users. The heat supply network comprises a water supply network and a water return network which take water or steam as medium, and connects a heat source with heat users. The variables involved in the modeling simulation of the heating network are as follows: pressure and flow in the pipeline, water supply temperature, backwater temperature and thermal power of pipe network nodes.
The existing solving method aiming at the hydraulic and thermal balance points of the heating pipe network is mature, a graph theory method is generally adopted, a specific pipe network is abstracted into a graph, the pipe network can be simplified into a pipe network model comprising pipe sections and node two elements by utilizing the concept of node flow, the attributes of the nodes and the pipe sections are described by a correlation matrix and a basic loop matrix in mathematics, the continuous equation of the node flow and the pressure loss equation of the loop are combined, and the hydraulic and thermal solving of the heating pipe network is carried out by an iteration method such as Newton-Lafson.
Where matrix a= (a) ij ) J×N A is the incidence matrix (J rows, N columns) of the heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the nodes and the pipeline number of the heat supply pipeline network, wherein the nodes comprise heat source nodes and loadsNodes and general nodes; m is the flow of each pipeline, m q Is the consumed traffic of the node.
B＝(b kl ) (N-J+1，N) Is a heat supply network basic loop matrix (N-J+1 rows, N columns), b kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
j and N respectively represent the node and pipeline number of the heat supply pipe network, and N-J+1 represents the loop number of the heat supply pipe network.
h f Representing the head loss of the pipes in the basic circuit of the heating network.
At present, no simulation method or scheme exists for the change of the hydraulic and thermal balance points of the heat supply network along with the load and the heat source.
The existing mathematical modeling method for the heat supply pipe network is mature. However, the heating network model is relatively independent, and cannot provide an input/output interface for the comprehensive energy system simulation.
In the prior art, when solving a hydraulic/thermal model of a heat supply network, load nodes are often simplified, and the influence caused by the change of the load nodes along with external factors such as time, ambient temperature and the like is often not considered, so that the load nodes are inconsistent with the actual situation.
Disclosure of Invention
In order to solve the defects in the prior art, the application provides a heating network simulation method and a heating network simulation device for a comprehensive energy system, and solves the problem that the current heating network simulation method does not consider the change of the hydraulic and thermal balance points of the heating network along with the load and the heat source.
In order to achieve the above object, the present application adopts the following technical scheme: a heating pipe network simulation method for a comprehensive energy system comprises the following steps:
constructing a heating network module, the heating network module comprising: input quantity, output quantity and characteristic parameters;
constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module;
and according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the heat supply pipe network is simulated by considering the change of load nodes and/or heat source nodes, and the hydraulic-thermal balance point of the heat supply pipe network is obtained.
Further, the input amount includes: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: total number of heat supply pipe network nodes, inner diameter of the pipeline, length of the pipeline, roughness coefficient and heat dissipation coefficient.
Further, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix are as follows:
where matrix a= (a) ij ) J×N A is an incidence matrix of a heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node;
B＝(b kl ) (N-J+1，N) b, a basic loop matrix of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f representing the head loss of the pipes in the basic circuit of the heating network.
Further, the simulating the heating pipe network by considering the load node change to obtain the hydraulic-thermal balance point of the heating pipe network comprises:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to carry out simulation calculation, when the output of an internal factor module and an external factor module of a certain simulation step changes, the load module takes the changed output of the internal factor module and the external factor module as own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value;
in the step length, the load module transmits the calculated new load node power as output to the heat supply pipe network module, the heat supply pipe network module uses the new load node power as input value, and the Newton-Lafson method is adopted to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, so as to obtain new hydraulic and thermal balance points;
the method for simulating the heat supply pipe network by considering the heat source node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to perform simulation calculation, and when the output of a thermoelectric coupling element module with a certain simulation step changes, the heat source change of the thermoelectric coupling element module is transmitted to a heating pipe network module as a new output value;
and in the step length, the heat supply pipe network module adopts a Newton-Lapherson method to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, and a new hydraulic and thermal balance point is obtained.
Further, the method also comprises the steps of delaying the new hydraulic and thermal balance points by delta N steps and then transmitting the delayed hydraulic and thermal balance points to an external module,
ΔN＝floor(T/timestep)
wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation steps of the delay.
A heating network simulation device for a comprehensive energy system, comprising:
constructing a heating network module, the heating network module comprising: input quantity, output quantity and characteristic parameters;
constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module;
and according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the heat supply pipe network is simulated by considering the change of load nodes and/or heat source nodes, and the hydraulic-thermal balance point of the heat supply pipe network is obtained.
Further, the input amount includes: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: total number of heat supply pipe network nodes, inner diameter of the pipeline, length of the pipeline, roughness coefficient and heat dissipation coefficient.
Further, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix are as follows:
where matrix a= (a) ij ) J×N A is an incidence matrix of a heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node;
B＝(b kl ) (N-J+1，N) b, a basic loop matrix of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f representing the head loss of the pipes in the basic circuit of the heating network.
Further, the simulating the heating pipe network by considering the load node change to obtain the hydraulic-thermal balance point of the heating pipe network comprises:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to carry out simulation calculation, when the output of an internal factor module and an external factor module of a certain simulation step changes, the load module takes the changed output of the internal factor module and the external factor module as own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value;
and in the step length, the load module transmits the calculated new load node power as output to the heat supply pipe network module, the heat supply pipe network module uses the new load node power as input value, and the Newton-Lafson method is adopted to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, so as to obtain a new hydraulic and thermal balance point.
The method for simulating the heat supply pipe network by considering the heat source node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to perform simulation calculation, and when the output of a thermoelectric coupling element module with a certain simulation step changes, the heat source change of the thermoelectric coupling element module is transmitted to a heating pipe network module as a new output value;
and in the step length, the heat supply pipe network module adopts a Newton-Lapherson method to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, and a new hydraulic and thermal balance point is obtained.
Further, the method also comprises the steps of delaying the new hydraulic and thermal balance points by delta N steps and then transmitting the delayed hydraulic and thermal balance points to an external module,
ΔN＝floor(T/timestep)
wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation steps of the delay.
The application has the beneficial effects that:
(1) The application realizes the development of the heat supply pipe network module based on the function of the TRNSYS rapid development module, provides a standardized input/output interface and establishes the connection between the heat supply pipe network module and other modules;
(2) According to the application, external factors are used as input variables of a load model, the influence of external factor changes on load nodes and the influence of comprehensive energy system operation mode changes on a heat supply source are fully considered, and the load nodes and/or the heat source node changes are considered to simulate a heat supply pipe network, so that the accuracy of comprehensive energy simulation is improved;
(3) The application fully considers the time delay of the heat supply pipe network, so that the simulation result is more approximate to the actual operation.
Drawings
FIG. 1 is a flow chart of a heating network simulation method in a specific embodiment of the application;
fig. 2 is a schematic diagram of a solution method for hydraulic and thermal balance points of a heating network in a specific embodiment of the application.
Detailed Description
The application is further described below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present application, and are not intended to limit the scope of the present application.
The integrated energy simulation system adopts the modular idea and is constructed by a TRNSYS (transient system simulation program). In TRNSYS, a comprehensive energy system is organically combined by a plurality of small modules, each module defines input and output variables and physical parameters of the module, simulation personnel is allowed to define and assign the characteristics of the module according to actual engineering equipment parameters, and connection of the comprehensive energy system can be completed through standardized input and output interfaces of each module. After connection is completed, a simulation personnel can set a simulation time step and total simulation time, so that TRNSYS calls each module to complete simulation calculation by taking the simulation time step as an interval in the total simulation time.
Meanwhile, the TRNSYS has the function of rapidly developing a new module, and a developer can complete the development of the new module according to own requirements.
The existing load module and various thermoelectric coupling element modules in the TRNSYS have no heating pipe network simulation module, and the embodiment realizes a heating pipe network module applicable to a comprehensive energy system based on the rapid development function of the TRNSYS.
Example 1:
as shown in fig. 1, a heating network simulation method suitable for a comprehensive energy system includes the steps:
step 1, constructing a heat supply pipe network module, packaging the heat supply pipe network module and adding the heat supply pipe network module to a TRNSYS module library, wherein the heat supply pipe network module comprises: input quantity, output quantity and characteristic parameters which are all set in a built-in function interface of TRNSYS;
the input quantity includes: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: total number of heat supply pipe network nodes, inner diameter of the pipeline, length of the pipeline, roughness coefficient and heat dissipation coefficient.
Step 2, constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters of the heat supply pipe network module;
where matrix a= (a) ij ) J×N A is the incidence matrix (J rows, N columns) of the heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node.
B＝(b kl ) (N-J+1，N) B is a basic loop matrix (N-J+1 rows and N columns) of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f the pressure head loss of the pipeline in the basic loop of the heat supply pipeline network is represented, and the pressure head loss is determined by the inner diameter of the pipeline, the length of the pipeline, the roughness coefficient and the heat dissipation coefficient.
And 3, simulating the heating pipe network according to the constructed heating pipe network module, the heating pipe network association matrix and the heating pipe network basic loop matrix by considering the change of load nodes and/or heat source nodes, so as to obtain the hydraulic-thermal balance point of the heating pipe network.
The calculation method of the hydraulic-thermal balance point of the heating pipe network comprises the following steps: according to the input quantity of the heating pipe network module, the heating pipe network association matrix and the heating pipe network basic loop matrix, carrying out the hydraulic-thermal joint iterative solution of the heating pipe network by adopting a Newton-Lapherson method until the iterative result meets the convergence condition; and transmitting the solving result to external modules such as load, heat source and the like, wherein the solving result is output quantity.
As shown in fig. 2:
the method for solving the hydraulic and thermal balance points of the heating pipe network by considering the change of the load node comprises the following steps:
1) TRNSYS calls corresponding modules in a certain order according to a fixed simulation time step to perform simulation calculation, and when the output of an internal factor module and an external factor module of a certain simulation step (for example: weather, personnel flow in building (load), etc.), the load module takes the output of the change of the internal and external factor modules as its own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value.
The internal and external factor modules are used for simulating the changing conditions of weather, personnel flow in the load and the like;
the load module is used for calculating the thermal power of the new load node and transmitting the thermal power value of the new load node to each load node of the heat supply pipe network module.
2) In the step length, the load module transmits the calculated load node power as output to the heat supply pipe network module, the heat supply pipe network module takes the new load node power as input value, the heat supply pipe network hydraulic-thermal joint iteration solution is carried out by adopting the Newton-Lafson method through the constructed heat supply pipe network joint matrix and the heat supply network basic loop matrix until the iteration result meets the convergence condition, and the new hydraulic and thermal balance points (comprising: the flow of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network) are used as the output quantity of the step length and transmitted to external modules such as loads, heat sources and the like.
The method for solving the hydraulic and thermal balance points of the heating network by considering the change of the heat source comprises the following steps:
1) TRNSYS calls the corresponding modules in a certain order according to a fixed simulation time step to perform simulation calculation, and when the output of the thermoelectric coupling element module with a certain simulation step (for example: the temperature of the heat source and the flow rate of the heat source), the heat source change of the thermoelectric coupling element module is used as a new output value to be transmitted to the heat supply pipe network module.
The thermoelectric coupling element module is used for simulating and calculating the change condition of the heat source and transmitting parameters such as the temperature, the flow and the like of the heat source to the heat source node of the heat supply pipe network module.
2) In this step, the heat source of the heat supply pipe network module changes (e.g.: heat source temperature and heat source flow) as input values, and carrying out hydraulic-thermal joint iterative solution on the heat supply pipe network by adopting the Newton-Lapherson method through the constructed heat supply gateway joint matrix and the heat supply network basic loop matrix until the iterative result meets the convergence condition. And solving new hydraulic and thermal balance points (comprising the flow of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network) as output of the step length, and transmitting the output to external modules such as loads, heat sources and the like.
The application also fully considers the time delay of energy transfer of the heat supply pipe network, the calculated thermodynamic and hydraulic balance points are not immediately transferred to each load module, and a time delay module is added between the heat supply pipe network module and the load modules. The delay characteristics are as follows:
ΔN＝floor(T/timestep)
wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation step sizes of delay, namely, after the delay delta N step sizes of new hydraulic and thermal balance points, the delay delta N step sizes are transmitted to a load module.
Example 2:
a heating network simulation device for a comprehensive energy system, comprising:
constructing a heating network module, the heating network module comprising: input quantity, output quantity and characteristic parameters;
constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module;
and according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the heat supply pipe network is simulated by considering the change of load nodes and/or heat source nodes, and the hydraulic-thermal balance point of the heat supply pipe network is obtained.
Further, the input amount includes: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: total number of heat supply pipe network nodes, inner diameter of the pipeline, length of the pipeline, roughness coefficient and heat dissipation coefficient.
Further, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix are as follows:
where matrix a= (a) ij ) J×N A is an incidence matrix of a heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node;
B＝(b kl ) (N-J+1，N) b, a basic loop matrix of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f representing the head loss of the pipes in the basic circuit of the heating network.
Further, the simulating the heating pipe network by considering the load node change to obtain the hydraulic-thermal balance point of the heating pipe network comprises:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to carry out simulation calculation, when the output of an internal factor module and an external factor module of a certain simulation step changes, the load module takes the changed output of the internal factor module and the external factor module as own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value;
and in the step length, the load module transmits the calculated new load node power as output to the heat supply pipe network module, the heat supply pipe network module uses the new load node power as input value, and the Newton-Lafson method is adopted to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, so as to obtain a new hydraulic and thermal balance point.
The method for simulating the heat supply pipe network by considering the heat source node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to perform simulation calculation, and when the output of a thermoelectric coupling element module with a certain simulation step changes, the heat source change of the thermoelectric coupling element module is transmitted to a heating pipe network module as a new output value;
and in the step length, the heat supply pipe network module adopts a Newton-Lapherson method to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, and a new hydraulic and thermal balance point is obtained.
Further, the method also comprises the steps of delaying the new hydraulic and thermal balance points by delta N steps and then transmitting the delayed hydraulic and thermal balance points to an external module,
ΔN＝floor(T/timestep)
wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation steps of the delay.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing is merely a preferred embodiment of the present application, and it should be noted that modifications and variations could be made by those skilled in the art without departing from the technical principles of the present application, and such modifications and variations should also be regarded as being within the scope of the application.
Claims (4)
1. A heating pipe network simulation method for a comprehensive energy system is characterized by comprising the following steps of: comprising the following steps:
constructing a heating network module based on TRNSYS, the heating network module comprising: input quantity, output quantity and characteristic parameters;
constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module;
according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the load node and/or the heat source node change is considered to simulate the heat supply pipe network, so as to obtain a hydraulic-thermal balance point of the heat supply pipe network;
the input amounts include: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: the total number of heat supply pipe network nodes, the inner diameter of a pipeline, the length of the pipeline, the roughness coefficient and the heat dissipation coefficient;
the method for simulating the heat supply pipe network by considering the load node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to carry out simulation calculation, when the output of an internal factor module and an external factor module of a certain simulation step changes, the load module takes the changed output of the internal factor module and the external factor module as own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value;
in the step length, the load module transmits the calculated new load node power as output to the heat supply pipe network module, the heat supply pipe network module uses the new load node power as input value, and the Newton-Lafson method is adopted to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, so as to obtain new hydraulic and thermal balance points;
the method for simulating the heat supply pipe network by considering the heat source node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to perform simulation calculation, and when the output of a thermoelectric coupling element module with a certain simulation step changes, the heat source change of the thermoelectric coupling element module is transmitted to a heating pipe network module as a new output value;
in the step length, the heat supply pipe network module adopts Newton-Lapherson method to carry out hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, and a new hydraulic and thermal balance point is obtained;
the method also comprises the steps of delaying the new hydraulic and thermal balance points for delta N steps, and then transmitting the delayed hydraulic and thermal balance points to an external module, wherein delta N=floor (T/timer)
Wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation steps of the delay.
2. The heating network simulation method for the comprehensive energy system according to claim 1, wherein the heating network simulation method is characterized in that: the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix are as follows:
where matrix a= (a) ij ) J×N A is an incidence matrix of a heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node;
B＝(b kl ) (N-J+1，N) b, a basic loop matrix of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f representing the head loss of the pipes in the basic circuit of the heating network.
3. A heating network simulation device for synthesizing energy system, characterized by: comprising the following steps:
constructing a heating network module based on TRNSYS, the heating network module comprising: input quantity, output quantity and characteristic parameters;
constructing a heat supply pipe network association matrix and a heat supply pipe network basic loop matrix according to characteristic parameters in the heat supply pipe network module;
according to the constructed heat supply pipe network module, the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix, the load node and/or the heat source node change is considered to simulate the heat supply pipe network, so as to obtain a hydraulic-thermal balance point of the heat supply pipe network;
the input amounts include: the power required by the load node, the heat source node temperature and the heat source node flow;
the output quantity includes: the flow rate of each pipeline of the heat supply pipe network, the water supply and return temperature of the load node, the water supply and return temperature of the heat source node and the heat loss of the heat supply pipe network;
the characteristic parameters include: the total number of heat supply pipe network nodes, the inner diameter of a pipeline, the length of the pipeline, the roughness coefficient and the heat dissipation coefficient;
the method for simulating the heat supply pipe network by considering the load node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to carry out simulation calculation, when the output of an internal factor module and an external factor module of a certain simulation step changes, the load module takes the changed output of the internal factor module and the external factor module as own input, and calculates the load node power required by the step load by adopting a volumetric heat index method according to the changed input value;
in the step length, the load module transmits the calculated new load node power as output to the heat supply pipe network module, the heat supply pipe network module uses the new load node power as input value, and the Newton-Lafson method is adopted to carry out the hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, so as to obtain new hydraulic and thermal balance points;
the method for simulating the heat supply pipe network by considering the heat source node change to obtain the hydraulic-thermal balance point of the heat supply pipe network comprises the following steps:
the TRNSYS calls corresponding modules in a certain sequence according to a fixed simulation time step to perform simulation calculation, and when the output of a thermoelectric coupling element module with a certain simulation step changes, the heat source change of the thermoelectric coupling element module is transmitted to a heating pipe network module as a new output value;
in the step length, the heat supply pipe network module adopts Newton-Lapherson method to carry out hydraulic-thermal joint iterative solution of the heat supply pipe network through the heat supply pipe network association matrix and the basic loop matrix until the iterative result meets the convergence condition, and a new hydraulic and thermal balance point is obtained;
the method also comprises the steps of delaying the new hydraulic and thermal balance points for delta N steps, and then transmitting the delayed hydraulic and thermal balance points to an external module, wherein delta N=floor (T/timer)
Wherein floor represents a rounding down; t represents the flow time from a heat supply source to a certain load hot water under a new hydraulic power and thermal balance point; timetep represents the simulation step size; delta N represents the number of simulation steps of the delay.
4. A heating network simulation device for a comprehensive energy system according to claim 3, wherein: the heat supply pipe network association matrix and the heat supply pipe network basic loop matrix are as follows:
where matrix a= (a) ij ) J×N A is an incidence matrix of a heat supply pipe network ij An element representing the ith row and j column of the incidence matrix of the heating network, a ij The value of (2) is determined by:
j and N respectively represent the node and the pipe section number of the heat supply pipe network; m is the flow of each pipeline, m q Is the consumed traffic for each node;
B＝(b kl ) (N-J+1，N) b, a basic loop matrix of a heat supply pipe network kl Elements representing the kth row/column of the basic loop matrix of the heating network, b kl The value of (2) is determined by:
h f representing the head loss of the pipes in the basic circuit of the heating network.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011108790.7A CN112182905B (en) | 2020-10-16 | 2020-10-16 | Heat supply pipe network simulation method and device for comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011108790.7A CN112182905B (en) | 2020-10-16 | 2020-10-16 | Heat supply pipe network simulation method and device for comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112182905A CN112182905A (en) | 2021-01-05 |
| CN112182905B true CN112182905B (en) | 2023-08-29 |
Family
ID=73950634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011108790.7A Active CN112182905B (en) | 2020-10-16 | 2020-10-16 | Heat supply pipe network simulation method and device for comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112182905B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115906587B (en) * | 2022-12-28 | 2024-02-20 | 南京聚润信息科技有限公司 | Method for establishing dynamic coupling model of heat supply pipe network of comprehensive energy system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110543661A (en) * | 2019-07-18 | 2019-12-06 | 国网江苏省电力有限公司 | correlation-considered probability energy flow calculation method for electricity-heat interconnection comprehensive energy system |
| CN110555264A (en) * | 2019-08-30 | 2019-12-10 | 南方电网科学研究院有限责任公司 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
| WO2020093296A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Interval power flow calculation method for power-heat integrated energy system |
| CN111191182A (en) * | 2019-12-17 | 2020-05-22 | 东南大学 | Static combined heat and power flow hybrid solving method based on linearized thermal model |
| CN111222213A (en) * | 2020-01-15 | 2020-06-02 | 许继集团有限公司 | Dynamic simulation method and device for thermodynamic network |
-
2020
- 2020-10-16 CN CN202011108790.7A patent/CN112182905B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020093296A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Interval power flow calculation method for power-heat integrated energy system |
| CN110543661A (en) * | 2019-07-18 | 2019-12-06 | 国网江苏省电力有限公司 | correlation-considered probability energy flow calculation method for electricity-heat interconnection comprehensive energy system |
| CN110555264A (en) * | 2019-08-30 | 2019-12-10 | 南方电网科学研究院有限责任公司 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
| CN111191182A (en) * | 2019-12-17 | 2020-05-22 | 东南大学 | Static combined heat and power flow hybrid solving method based on linearized thermal model |
| CN111222213A (en) * | 2020-01-15 | 2020-06-02 | 许继集团有限公司 | Dynamic simulation method and device for thermodynamic network |
Non-Patent Citations (1)
| Title |
|---|
| 典型综合能源系统仿真与规划平台综述;王梦雪;电网技术(第12期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112182905A (en) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106647329B (en) | Equivalent modeling method for heat supply pipe network | |
| Gabrielaitiene | Numerical simulation of a district heating system with emphases on transient temperature behaviour | |
| Li et al. | Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors | |
| CN109376428A (en) | Reliability estimation method, device, equipment and the storage medium of integrated energy system | |
| Leśko et al. | Modeling of district heating networks for the purpose of operational optimization with thermal energy storage | |
| CN112182905B (en) | Heat supply pipe network simulation method and device for comprehensive energy system | |
| CN110795845A (en) | Heat supply system margin quantitative calculation method and system | |
| Dénarié et al. | Dynamical modelling and experimental validation of a fast and accurate district heating thermo-hydraulic modular simulation tool | |
| CN111783309A (en) | Dynamic simulation method of steam heating network based on internal conservation | |
| CN112100778A (en) | Scheduling method and system of comprehensive energy system | |
| CN115374527A (en) | Construction method of regional electricity-heat-cold comprehensive energy dynamic simulation system | |
| CN115935566A (en) | Simulation method and system for natural gas pipeline network, storage medium and electronic equipment | |
| CN115079592A (en) | Pipe network simulation method for thermodynamic system of ship nuclear power device | |
| CN114021490A (en) | Dynamic full-analysis method for hot water network | |
| Anderson | Coupling RELAP5-3D and Fluent to analyze a Very High Temperature Reactor (VHTR) outlet plenum | |
| Dalvi et al. | Integrated system model of district cooling for energy consumption optimization | |
| Blizard et al. | A Graph-Based Technique for the Automated Control-Oriented Modeling of District Heating Networks | |
| CN113706094B (en) | Comprehensive energy real-time collaborative simulation system and method based on message bus | |
| Sorocins et al. | District heating simulation model development to solve optimization problems in the market conditions | |
| Landman et al. | Dynamic Systems CFD Simulation Code for the Modeling of HTGR Power Plants | |
| CN109611949B (en) | Pipeline thermal calculation method and monitoring and heating system adopting same | |
| Pan et al. | Quasi-dynamic interactions between electricity systems and heating systems in normal operation | |
| CN113268824B (en) | Fluid loop numerical simulation method, system and medium based on flow field in temperature control valve | |
| CN116187586B (en) | Distributed planning solving method and system for comprehensive energy system | |
| Chen et al. | Quasi-dynamic modeling and simulation of district heating systems: A laplace transform-based approach |
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 |