CN110555264B - Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method - Google Patents
Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method Download PDFInfo
- Publication number
- CN110555264B CN110555264B CN201910816039.3A CN201910816039A CN110555264B CN 110555264 B CN110555264 B CN 110555264B CN 201910816039 A CN201910816039 A CN 201910816039A CN 110555264 B CN110555264 B CN 110555264B
- Authority
- CN
- China
- Prior art keywords
- branch
- thermodynamic
- equation
- simulation
- thermal
- 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
Images
Abstract
The invention discloses a dynamic simulation method and a dynamic simulation device of a comprehensive energy heating system based on an extended node method, which are applied to the comprehensive energy heating system, and the method comprises the following steps: the method comprises the steps of inputting system parameters, obtaining a hydraulic part dynamic simulation model, obtaining a pipeline branch thermodynamic equation set, obtaining a heat user branch thermodynamic equation set, obtaining a heat source branch thermodynamic equation set, obtaining a system thermodynamic topology constraint equation set, obtaining a thermodynamic part dynamic simulation model, simulating an initial step, a simulation processing step, a simulation judging step, a simulation result determining step and a simulation circulating step.
Description
Technical Field
The invention relates to the field of data simulation, in particular to a dynamic simulation method and device of a comprehensive energy heating system based on an extended node method.
Background
With the increasing awareness of national energy conservation and emission reduction, the comprehensive energy heating system for improving the energy utilization rate through comprehensive utilization of energy becomes an important development direction of the future smart grid.
Because the comprehensive energy heating system relates to hydraulic dynamic and thermal dynamic at the same time, the simulation difficulty of the comprehensive energy heating system is high, and therefore a method capable of accurately simulating the comprehensive energy heating system is urgently needed by technical personnel in the field.
Disclosure of Invention
In view of the above problems, the present invention provides a method and an apparatus for dynamically simulating an integrated energy heating system based on an extended node method, which overcome or at least partially solve the above problems, and the technical solution is as follows:
a dynamic simulation method of an integrated energy heating system based on an extended node method is applied to the integrated energy heating system, and comprises the following steps: the method comprises the steps of inputting system parameters, obtaining a hydraulic part dynamic simulation model, obtaining a pipeline branch thermodynamic equation set, obtaining a heat user branch thermodynamic equation set, obtaining a heat source branch thermodynamic equation set, obtaining a system thermodynamic topology constraint equation set, obtaining a thermodynamic part dynamic simulation model, simulating an initial step, a simulation processing step, a simulation judging step, a simulation result determining step and a simulation circulating step,
the system parameter inputting step comprises: inputting system parameters of the comprehensive energy heat supply system, wherein the system parameters comprise at least one of a branch length, a branch inner wall sectional area, a branch thermal resistance, a branch resistance coefficient, a pressure increment caused by a power source in a branch and a node equivalent flow capacity, the thermal working medium parameters can comprise at least one of a thermal working medium density and a thermal working medium specific heat capacity, the heat source parameters can comprise at least one of a heat source internal heat working medium mass, a heat source heating power and a heat source heating efficiency, the heat user parameters can comprise at least one of a radiator area, a radiator heat dissipation coefficient, a heat source internal heat working medium mass, a heat supply building heat dissipation coefficient, an air density and an air capacity, and the simulation calculation parameters can comprise at least one of a simulation specific heat length, a simulation step length and a pipeline micro element length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system;
the step of obtaining the dynamic simulation model of the hydraulic part comprises the following steps: according to the system parameters, a system node flow conservation equation and a system branch momentum conservation equation are obtained; according to an implicit Euler method, carrying out differencing processing on the system node flow conservation equation and the system branch momentum conservation linearized equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation, and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system;
the step of obtaining the pipeline branch thermodynamic equation set comprises the following steps: according to the system parameters, acquiring a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch;
the step of obtaining the thermal equation set of the hot user branch comprises the following steps: according to the system parameters, obtaining thermal user branch thermodynamic equations of all thermal user branches in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch;
the step of obtaining the heat source branch thermodynamic equation set comprises the following steps: according to the system parameters, obtaining a heat source branch thermodynamic equation of each heat source branch in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the heat source branch thermodynamic equations of the heat source branches according to an implicit Euler method to obtain heat source branch discretization thermodynamic equations of the heat source branches; obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch;
the step of obtaining the system thermodynamic topological constraint equation set comprises the following steps: according to the system parameters, acquiring a thermodynamic temperature mixing equation of each node in the comprehensive energy heating system; obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node;
the step of obtaining the dynamic simulation model of the thermodynamic part comprises the following steps: obtaining a dynamic simulation model of a thermal part of the comprehensive energy heating system according to the pipeline branch thermal equation set, the heat user branch thermal equation set, the heat source branch thermal equation set and the system thermal topological constraint equation set;
the simulation initial step comprises: setting the starting time of the simulation duration to be t =0, setting initial pressure, initial flow and initial temperature, taking the initial pressure as the current pressure, the initial flow as the current flow, and the initial temperature as the current temperature;
the simulation processing step comprises: inputting the current pressure and the current flow into the hydraulic part dynamic simulation model to obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part;
the simulation judging step comprises: after waiting for the simulation step length delta t, judging whether the current time reaches the simulation termination time of the simulation duration, if so, executing the simulation result determining step, and if not, executing the simulation circulating step;
the simulation result determining step comprises: determining the pressure, flow and temperature obtained by each simulation as a simulation result;
the simulation loop step includes: and taking the first pressure as the current pressure, taking the first flow as the current flow, taking the first temperature as the current temperature, and returning to execute the simulation processing step.
Optionally, the hydraulic part dynamic simulation model is as follows:
wherein, the first and the second end of the pipe are connected with each other,
in the formula, H (t) is a column vector formed by pressure increment caused by a power source in a branch, is a known input quantity, and A is a node-branch incidence matrix of the comprehensive energy heating system: that is, if the branch k flows from the node i to the node j, a (i, k) =1,A (j, k) = -1, the remaining elements in the k-th column are 0, and the node with the node number of 0 is used as a reference node to obtain a (a-1) × b-dimensional node-branch reduced price incidence matrix a'; p (t) is a node pressure column vector at the time t; g (t) is a branch flow column vector at the time t; s is a resistance coefficient matrix which is a diagonal matrix, and main diagonal elements of the matrix are branch resistance coefficients of corresponding branches; h is a pressure increment column vector, and the element of the pressure increment column vector is the pressure increment caused by the power source in the corresponding branch; l is a radical of an alcohol m The main diagonal element of (A) is the inertia coefficient of the corresponding branch, which can be determined by the branch length l k And the cross-sectional area A of the inner wall of the branch c，k And calculating to obtain the following formula:
L m is a diagonal matrix; the main diagonal element of C is the equivalent flow capacity of the corresponding node and is a diagonal matrix; f (g) is a square function of branch flow, and J (t-delta t) is a Jacobian matrix of f (g) at t-delta t.
Optionally, the pipeline branch thermodynamic equation set is as follows:
wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
wherein q is the number of pipeline branches I q Is a q-dimensional identity matrix; t is e (t) is the temperature column vector of the micro-element section of all the pipeline branches; t is s (t) column vectors of water supply temperatures for all the branches; t is r (t) is the array vector of return water temperature of all branches; a. The se And A re Is composed of
A dimension matrix; u. of gd Number of micro-elements, T, of a branch of a pipeline s，k The temperature of the water supply for branch k; t is r，k The return water temperature of the branch k; t is e Is the temperature of the infinitesimal stage; Δ x is the length of the pipeline infinitesimal section; delta t is a simulation step length; a. The c，gd The sectional area of the inner wall of the pipeline branch gd; rho is the density of the hot working medium; c. C p The specific heat capacity of the thermal working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g gd The hot working medium flow rate of the pipeline branch gd; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
Optionally, the thermal equation set of the hot user branch is as follows:
wherein the content of the first and second substances,
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T b (t)＝[T b，1 (t) T b，2 (t) ... T b，v (t)] T
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A rl 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
in the formula, A sl 、A rl 、A bl 、A sz 、A rz 、A bz Are v-dimensional diagonal matrices, the main diagonal element pair i =1,2.
In the formula, T s Is a water supply temperature column vector; t is r Is a backwater temperature column vector; t is b Is a heating building temperature column vector; m is h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; g i+p Numbering the hot working medium flow of the hot user branch with the number of i + p for the hot user branch; k r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; delta t is a simulation step length; rho a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is o Is the outdoor temperature; q. q.s V Is a heat dissipation coefficient; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
Optionally, the thermal equation set of the heat source branch is as follows:
wherein the content of the first and second substances,
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
A sh 、A rh for an f-dimensional diagonal matrix, the main diagonal element pair i =1,2,.. Multidot.f, has:
in the formula, delta t is simulation step length; m is a unit of w Is the mass of hot working medium in the heat source; c. C p The specific heat capacity of the thermal working medium; t is s Is the temperature of the supplied water; t is a unit of r The temperature of the return water is; p W And eta are the heating power and the heating efficiency of the heat source respectively; g i+p+v The hot working medium flow of the heat source branch with the serial number of i + p + v is numbered for the heat source branch; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
Optionally, the system thermodynamic topological constraint equation set is as follows:
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T nod (t)＝[T 1 (t) T 2 (t) ... T a (t)] T
in the formula I a Is a dimension unit matrix; i is b Is a b-dimensional identity matrix; a. The rn Is a matrix of a x b dimension; a is the number of nodes, b is the number of branches and for i =1,2,.. A, k =1,2 s Is the temperature of the supplied water; t is r The temperature of the return water is; a. The sn Is a b × a dimensional matrix, and for i =1,2,.. A, k =1,2.
Optionally, the dynamic simulation model of the thermal part is as follows:
wherein, I a 、A rn 、A sn 、I b A coefficient matrix which is a system thermodynamic topological constraint equation set; a. The sh0 、A rh0 A coefficient matrix of the heat source branch thermodynamic equation set is obtained; a. The sl0 、A rl0 、A sz0 、A rz0 、A bl 、A bz A coefficient matrix of the thermal user branch thermodynamic equation set; a. The e 、I q0 、A se 、A re And the coefficient matrix of the pipeline branch thermodynamic equation set.
The utility model provides a comprehensive energy heating system dynamic simulation device based on extended node method, is applied to comprehensive energy heating system, the device includes: a system parameter input unit, a hydraulic part dynamic simulation model obtaining unit, a pipeline branch thermodynamic equation set obtaining unit, a heat user branch thermodynamic equation set obtaining unit, a heat source branch thermodynamic equation set obtaining unit, a system thermodynamic topology constraint equation set obtaining unit, a thermodynamic part dynamic simulation model obtaining unit, a simulation initial unit, a simulation processing unit, a simulation judging unit, a simulation result determining unit and a simulation circulating unit,
the system parameter input unit is used for inputting system parameters of the comprehensive energy heat supply system, wherein the system parameters comprise a network topology connection relation of a heat supply pipe network of the comprehensive energy heat supply system, network parameters, thermal working medium parameters, heat source parameters, heat user parameters and simulation calculation parameters, the network parameters comprise at least one of branch length, a branch inner wall sectional area, branch thermal resistance, branch resistance coefficients, pressure increment caused by a power source in a branch and node equivalent flow capacity, the thermal working medium parameters can comprise at least one of thermal working medium density and thermal working medium specific heat capacity, the heat source parameters can comprise at least one of heat working medium quality, heat source heating power and heat source heating efficiency in a heat source, the heat user parameters can comprise at least one of heat radiator area, heat radiator heat radiation coefficient, heat working medium quality, volume, heat supply building heat radiation coefficient, air density and air specific heat capacity in the heat source, and the simulation calculation parameters can comprise at least one of simulation duration, simulation step length and pipeline infinitesimal length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system;
the hydraulic part dynamic simulation model obtaining unit is used for obtaining a system node flow conservation equation and a system branch momentum conservation equation according to the system parameters; according to an implicit Euler method, carrying out differencing processing on the system node flow conservation equation and the system branch momentum conservation linearized equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation, and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system;
the pipeline branch thermodynamic equation set obtaining unit is used for obtaining a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system according to the system parameters; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch;
the system comprises a heat consumer branch thermodynamic equation set obtaining unit, a heat capacity particle method and a heat capacity particle method, wherein the heat consumer branch thermodynamic equation set obtaining unit is used for obtaining heat consumer branch thermodynamic equations of all heat consumer branches in the comprehensive energy heating system according to the system parameters; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch;
the heat source branch thermodynamic equation set obtaining unit is used for obtaining a heat source branch thermodynamic equation of each heat source branch in the comprehensive energy heating system based on a thermal capacity particle method according to the system parameters; carrying out differential processing on the heat source branch thermodynamic equation of each heat source branch according to an implicit Euler method to obtain the heat source branch discretization thermodynamic equation of each heat source branch; obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch;
the system thermodynamic topological constraint equation set obtaining unit is used for obtaining a thermodynamic temperature mixed equation of each node in the comprehensive energy source heat supply system according to the system parameters; obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node;
the thermodynamic part dynamic simulation model obtaining unit is used for obtaining a thermodynamic part dynamic simulation model of the comprehensive energy heating system according to the pipeline branch thermodynamic equation set, the heat user branch thermodynamic equation set, the heat source branch thermodynamic equation set and the system thermodynamic topology constraint equation set;
the simulation initialization unit is configured to set a starting time of the simulation duration to be t =0, set an initial pressure, an initial flow and an initial temperature, use the initial pressure as a current pressure, use the initial flow as a current flow, and use the initial temperature as a current temperature;
the simulation processing unit is used for inputting the current pressure and the current flow into the hydraulic part dynamic simulation model to obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part;
the simulation judging unit is used for judging whether the current time reaches the simulation termination time of the simulation duration after waiting for the simulation step length delta t, if so, triggering the simulation result determining unit, and if not, triggering the simulation circulating unit;
the simulation result determining unit is used for determining the pressure, flow and temperature obtained by each simulation as a simulation result;
and the simulation circulating unit is used for triggering the simulation processing unit by taking the first pressure as the current pressure, taking the first flow as the current flow and taking the first temperature as the current temperature.
By means of the technical scheme, the invention provides a dynamic simulation method and a dynamic simulation device of an integrated energy heating system based on an extended node method, which are applied to the integrated energy heating system, and the method comprises the following steps: the method comprises the steps of inputting system parameters, obtaining a hydraulic part dynamic simulation model, obtaining a pipeline branch thermodynamic equation set, obtaining a heat user branch thermodynamic equation set, obtaining a heat source branch thermodynamic equation set, obtaining a system thermodynamic topology constraint equation set, obtaining a thermodynamic part dynamic simulation model, performing simulation initial, performing simulation processing, performing simulation judgment, determining a simulation result and performing simulation circulation.
The above description is only an overview of the technical solutions of the present invention, and the present invention can be implemented in accordance with the content of the description so as to make the technical means of the present invention more clearly understood, and the above and other objects, features, and advantages of the present invention will be more clearly understood.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
fig. 1 is a schematic flow chart of a dynamic simulation method of an integrated energy heating system based on an extended node method according to an embodiment of the present invention;
fig. 2 is a schematic structural diagram illustrating a network topology of an integrated energy heating system according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating a head-end water supply temperature and a tail-end water return temperature of a pipeline branch according to an embodiment of the present invention;
FIG. 4 is a diagram illustrating simulation results provided by an embodiment of the present invention;
FIG. 5 is a diagram illustrating another simulation result provided by an embodiment of the present invention;
FIG. 6 is a diagram illustrating another simulation result provided by an embodiment of the present invention;
FIG. 7 is a diagram illustrating another simulation result provided by an embodiment of the present invention;
fig. 8 shows a schematic structural diagram of a dynamic simulation apparatus of an integrated energy heating system based on an extended node method according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited by the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
As shown in fig. 1, a method for dynamically simulating an integrated energy heating system based on an extended node method according to an embodiment of the present invention may include: the method comprises the following steps of inputting system parameters S100, obtaining a hydraulic part dynamic simulation model S200, obtaining a pipeline branch thermodynamic equation set S300, obtaining a hot user branch thermodynamic equation set S400, obtaining a heat source branch thermodynamic equation set S500, obtaining a system thermodynamic topological constraint equation set S600, obtaining a thermodynamic part dynamic simulation model S700, simulating an initial step S800, simulating a processing step S900 and simulating a judging step S1000.
The system parameter input step S100 includes: inputting system parameters of the comprehensive energy heat supply system, wherein the system parameters comprise at least one of a branch length, a branch inner wall sectional area, a branch thermal resistance, a branch resistance coefficient, a pressure increment caused by a power source in a branch and a node equivalent flow capacity, the thermal working medium parameters can comprise at least one of a thermal working medium density and a thermal working medium specific heat capacity, the heat source parameters can comprise at least one of a heat source internal heat working medium mass, a heat source heating power and a heat source heating efficiency, the heat user parameters can comprise at least one of a radiator area, a radiator heat dissipation coefficient, a heat source internal heat working medium mass, a heat supply building heat dissipation coefficient, an air density and an air capacity, and the simulation calculation parameters can comprise at least one of a simulation specific heat length, a simulation step length and a pipeline micro element length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system.
Optionally, the process for obtaining the node-branch incidence matrix and the node-branch price reduction incidence matrix of the integrated energy heating system in this embodiment may be:
selecting any node in a pipe network as a pressure reference node, numbering other nodes except the reference node from 1, and recording the total number of the node numbers as a-1; all branches are numbered in sequence from 1, and the number of pipeline branches is q, the number of hot user branches is v, the number of hot source branches is f, and the total number of branches is = q + v + f. Wherein: the number of the pipeline branch is from 1 to q; the hot user branch is numbered from q +1 to q + v; the heat source legs are numbered from q + v +1 to b. For example, the number of the pipe branch in this embodiment may be 1 to 13, the number of the hot user branch may be 14 to 26, and the number of the hot source branch may be 27 to 36.
According to the topological connection relation of the pipe network, an a multiplied by b dimension node-branch incidence matrix A can be established: that is, if a branch k flows from the node jdn to the node jdm, a (jdn, k) =1,A (jdm, k) = -1, the remaining elements in the k-th column are 0, and the node with the node number of 0 is used as a reference node, so as to obtain the (a-1) × b-dimensional node-branch reduced price correlation matrix a.
The type of the branch in the present embodiment may include at least one type of a pipe branch, a hot user branch, and a hot source branch.
The step S200 of obtaining a dynamic simulation model of the hydraulic part includes: according to the system parameters, a system node flow conservation equation and a system branch momentum conservation equation are obtained; according to an implicit Euler method, carrying out differential processing on the system node flow conservation equation and the system branch momentum conservation equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation; and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system.
Specifically, the process of step S200 of obtaining the hydraulic portion dynamic simulation model may include:
and S210, obtaining a system node flow conservation equation according to the system parameters.
Specifically, the system node flow conservation equation is given by formula 1:
wherein, C is a node equivalent flow capacity matrix, the main diagonal element of the node equivalent flow capacity matrix is the equivalent flow capacity corresponding to the node, and the non-main diagonal element is 0.p is a node pressure column vector whose elements are the pressure corresponding to the node relative to a reference point; g is the branch flow column vector. Differentiating the formula 1 by using an implicit Euler method, and obtaining a system node flow conservation discretization equation after sorting:
p (t) is a node pressure column vector at a time t to be solved, g (t) is a branch flow column vector at the time t to be solved, and p (t-delta t) is a node pressure historical quantity which is obtained by calculation at a time before the time t to be solved and is separated by a simulation step length delta t from the time t to be solved.
And S220, obtaining a system branch momentum conservation equation according to the system parameters.
Specifically, the system branch momentum conservation equation is formula 3:
in the formula, L m The branch inertia coefficient matrix is obtained by taking the main diagonal element as the inertia coefficient of the corresponding branch and taking the non-main diagonal element as 0. Taking the branch k as an example, the inertia coefficient thereof can be determined by the branch length lk and the branch inner wall cross-sectional area A c，k And calculating to obtain the following formula:
A′ T transposing the node-branch reduced order incidence matrix obtained in the step 1; h is a pressure increment column vector, and the element of the pressure increment column vector is the pressure increment caused by the power source in the corresponding branch; s is a resistance coefficient matrix, the main diagonal element of the matrix is the resistance coefficient of the corresponding branch, and the non-main diagonal element of the matrix is 0; f (g) is a square function of the branch flow and can be written as:
taylor expansion is carried out on the state g (t-delta t) at the time t-delta t by the nonlinear term f (g), the first two terms are reserved, and a linearized expression of the nonlinear function f (g) at the time t-delta t can be obtained and further used for approximately calculating the next time step f (g (t)):
f (g (t)) ≈ f (g (t- Δ t)) + J (t- Δ t) [ g (t) -g (t- Δ t) ] equation 6
J (t-delta t) is a Jacobian matrix of f (g) in the state at the time t-delta t. Differentiating the formula 3 by using an implicit Euler method, substituting the formula 6, and sorting to obtain a system branch momentum conservation discretization equation:
and S230, combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system.
Specifically, the hydraulic part dynamic simulation model is formula 8:
wherein the content of the first and second substances,
in the formula, H (t) is a column vector formed by pressure increment caused by a power source in a branch, and is a known input quantity, and a is a node-branch incidence matrix of the comprehensive energy heating system: that is, if a branch k flows from the node jdn to the node jdm, a (jdn, k) =1,A (jdm, k) = -1, the remaining elements in the k-th column are 0, the node with the node number 0 is used as a reference node, and the (a-1) × b-dimensional node-branch reduced price is obtainedAn incidence matrix A'; p (t) is a node pressure column vector at time t; g (t) is a branch flow column vector at the time t; s is a resistance coefficient matrix which is a diagonal matrix, and main diagonal elements of the matrix are branch resistance coefficients of corresponding branches; h is a pressure increment column vector, and the element of the pressure increment column vector is the pressure increment caused by the power source in the corresponding branch; l is m The main diagonal element of (A) is the inertia coefficient of the corresponding branch, which can be determined by the branch length l k And the cross-sectional area A of the inner wall of the branch c，k And calculating to obtain the following formula:
lm is a diagonal matrix; the main diagonal element of C is the equivalent flow capacity of the corresponding node and is a diagonal matrix; f (g) is a square function of branch flow, and J (t-delta t) is a Jacobian matrix of f (g) at t-delta t.
And solving the formula 8 to obtain the pressure column vector of each node and the flow column vector of each branch of the comprehensive energy heating system at the time t.
The step S300 of obtaining the pipeline branch thermodynamic equation set includes: according to the system parameters, acquiring a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; and obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch.
Specifically, the step S300 of obtaining the pipeline branch thermodynamic equation set may include:
and S310, acquiring a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system according to the system parameters.
For any pipeline branch gd in each pipeline branch: according to system parameters, the pipeline branch thermodynamic equation for pipeline branch k may be:
wherein, T gd Is a pipeline temperature variable.
And S320, performing space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and performing time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain the pipeline branch discretization thermodynamic equation of each pipeline branch.
For ease of understanding, the pipeline branch gd is still used as an example here: to eliminate T in equation 9 gd Spatial partial derivative of
The pipeline branch gd is divided into a plurality of infinitesimal sections. According to the length delta x of a single pipeline infinitesimal and the length l of a pipeline branch gd gd The number of infinitesimal segments can be obtained
Is a pair of
And rounding up. For the wyd infinitesimal sections of the pipeline branch gd, the thermodynamic equation form is formula 10:
the boundary conditions are as follows:
wherein the content of the first and second substances,
the temperature of wyd infinitesimal sections of the pipeline branch gd; a. The c，gd The cross section area of the inner wall of the branch of the pipeline gd; rho is the density of the hot working medium; c. C p The specific heat capacity of the thermal working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g gd The hot working medium flow rate of the pipeline branch gd; in the dynamic simulation model of the thermal power part, pipeline branches, thermal user branches and heat source branches need to be considered respectively, and each pipeline branch gd defines the water supply temperature T at the head end according to the water flow direction s，gd And tail end return water temperature T r，gd As shown in fig. 3. T is a unit of s，gd Supply water temperature, T, to the head end of the pipeline branch gd r，gd The tail end return water temperature of the pipeline branch gd.
Differentiating equation 10 using the implicit Euler method, g gd (t) is the flow of the pipeline branch gd at the time t, and the discretized thermodynamic dynamic equation of wyd infinitesimal sections can be obtained by sorting:
wherein the content of the first and second substances,
writing the discretization dynamic equations and boundary conditions of all infinitesimal sections of the pipeline branch gd according to the sequence to obtain a discretization thermodynamic equation set of the pipeline branch gd:
wherein the content of the first and second substances,
the column vector is composed of all micro-element section temperatures of the pipeline branch gd and comprises:
s330, obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch.
Specifically, the pipeline branch discretization thermodynamic equations of all pipeline branches in the comprehensive energy heating system are arranged according to the sequence of the serial numbers of the pipeline branches, and a pipeline branch thermodynamic equation set is obtained.
The thermodynamic equation set of the pipeline branch is as follows:
wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
wherein, T e (t) is the temperature column vector of the micro-element section of all the pipeline branches; i is q Is a q-dimensional identity matrix; a. The se And A re Is composed of
A dimension matrix incorporating:
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
wherein q is the number of pipeline branches I q Is a q-dimensional identity matrix; t is e (t) is the temperature column vector of the micro-element section of all the pipeline branches; t is s (t) supplying water temperature column vectors for all the branches; t is r (t) is the array vector of return water temperature of all branches; a. The se And A re Is composed of
A dimension matrix; u. of gd Number of micro-elements, T, of a branch of a pipeline s，k The temperature of the water supply for branch k; t is r，k The return water temperature of the branch k; t is e Is the temperature of the infinitesimal section; Δ x is the length of the pipeline infinitesimal section; delta t is a simulation step length; a. The c，gd The sectional area of the inner wall of the pipeline branch gd; rho is the density of the hot working medium; c. C p Is the specific heat capacity of the hot working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g gd The hot working medium flow rate of the pipeline branch gd; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The step S400 of obtaining the thermal equation set of the hot user branch includes: according to the system parameters, obtaining thermal user branch thermodynamic equations of all thermal user branches in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; and obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch.
Specifically, the step S400 of obtaining the thermal equation set of the hot user branch may include:
and S410, obtaining thermal user branch thermodynamic equations of all thermal user branches in the comprehensive energy heating system based on a thermal capacity particle method according to system parameters.
For any hot user branch ryh in each hot user branch: according to system parameters, a thermal balance equation of the hot user branch rh is obtained based on a thermal capacity particle method and is as follows:
in the formula, m h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; t is s，ryh The water supply temperature of the hot user branch ryh; t is r，ryh The return water temperature of the hot user branch ryh; for each branch of heat user, the corresponding temperature of heat supply building can be defined, and the temperature vector T of heat supply building of heat user can be defined b ＝[T b，1 ，T b，2 ，...，T b，v ]The corresponding hot user branch is numbered from q +1 to q + v, T b，ryh-q The temperature of the heating building for the hot user branch ryh; q net Supplying heat to the heat supply network; q rad Heat dissipation for the radiator to the heating building; g ryh The hot working medium flow rate of a hot user branch ryh; k r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; rho a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is a unit of o Is the outdoor temperature; q heat Heat dissipation capacity to the outside for the heat supply building; q. q.s V The heat dissipation coefficient of the building is the heat dissipation coefficient of the heat supply.
And S420, carrying out differential processing on the thermal equation of the thermal user branch of each thermal user branch according to the implicit Euler method to obtain the discretization thermal equation of the thermal user branch of each thermal user branch.
For ease of understanding, the hot subscriber leg ryh is still used as an example for explanation here:
the implicit Euler method is used to differentiate the equation 15, and the equations in lines 3, 4 and 5 in the equation 15 are substituted into the equations in lines 1 and 2 to eliminate Q net 、Q rad And Q heat And after sorting, obtaining a discretization thermodynamic equation of the hot user branch ryh:
in the formula (I), the compound is shown in the specification,
and S430, obtaining a thermal equation set of the thermal user branches according to the discretization thermal equation of the thermal user branches.
Specifically, the discretization thermodynamic equations of the thermal user branches are arranged according to the serial number of the thermal user branches, and a thermodynamic equation set of the thermal user branches is obtained.
The thermal equation set of the hot user branch is as follows:
wherein the content of the first and second substances,
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A rl 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T b (t)＝[T b，1 (t) T b，2 (t) ... T b，v (t)] T
wherein A is sl 、A rl 、A bl 、A sz 、A rz 、A bz Are v-dimensional diagonal matrices, i.e. all non-principal diagonal elements are 0, and the principal diagonal element pair i =1,2.
In the formula, T s Is a water supply temperature column vector; t is r Is a backwater temperature column vector; t is b Is a heating building temperature column vector; m is h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; g i+p Thermal working medium flow of hot user branch with serial number of i + p for hot user branch；K r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; delta t is a simulation step length; rho a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is o Is the outdoor temperature; q. q of V Is a heat dissipation coefficient; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The step S500 of obtaining the thermal equation set of the heat source branch includes: according to the system parameters, obtaining a heat source branch thermodynamic equation of each heat source branch in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the heat source branch thermodynamic equations of the heat source branches according to an implicit Euler method to obtain heat source branch discretization thermodynamic equations of the heat source branches; and obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch.
Specifically, the step S500 of obtaining the thermal equation set of the heat source branch may include:
and S510, obtaining a heat source branch thermodynamic equation of each heat source branch in the comprehensive energy heating system based on a thermal capacity particle method according to system parameters.
Specifically, the heat source is regarded as a heat capacity temperature particle, and the obtained net heat is regarded as the difference between the heat brought by the heating power and the heat taken away by the supply water and the return water. For any heat source branch ry in each heat source branch: according to system parameters, a thermal capacity particle method is used for obtaining a thermal equation of a heat source branch k, wherein the thermal equation of the heat source branch k is as follows:
in the formula, m w Is the mass of hot working medium in the heat source; c. C p The specific heat capacity of the thermal working medium; t is s，ry The water supply temperature of the heat source branch ry; t is r，ry The return water temperature of the heat source branch ry; p W And eta are respectively the heating power and heating efficiency of the heat source, g ry The flow rate of the hot working medium in the heat source of the heat source branch ry.
S520, performing differential processing on the heat source branch thermodynamic equations of the heat source branches according to the implicit Euler method to obtain the heat source branch discretization thermodynamic equations of the heat source branches.
For the sake of understanding, the heat source branch k is still taken as an example for explanation:
differentiating the formula 18 by using an implicit Euler method, and sorting to obtain:
wherein, K sh ＝(-g ry (t)c p )，
S530, obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch.
Specifically, the heat source branch discretization thermodynamic equations of each heat source branch are arranged according to the serial number sequence of the heat source branches, and a heat source branch thermodynamic equation set is obtained.
The thermodynamic equation set of the heat source branch is as follows:
wherein the content of the first and second substances,
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
A sh0 、A rh0 is an f-dimensional diagonal matrix with 0 as the non-principal diagonal element, i =1,2, a principal diagonal element pair, f has:
in the formula, delta t is a simulation step length; m is w Is the mass of hot working medium in the heat source; c. C p The specific heat capacity of the thermal working medium; t is a unit of s，k The water supply temperature of the heat source branch k; t is r The temperature of the return water is; p is W And eta are respectively the heating power and the heating efficiency of the heat source; g i+p+v The hot working medium flow of the heat source branch with the serial number of i + p + v is numbered for the heat source branch; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The step S600 of obtaining a system thermodynamic topological constraint equation set includes: according to the system parameters, acquiring a thermodynamic temperature mixing equation of each node in the comprehensive energy heating system; and obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node.
Specifically, the step S600 of obtaining the system thermodynamic topology constraint equation set may include:
s610, according to the system parameters, a thermodynamic temperature mixing equation of each node in the comprehensive energy heating system is obtained.
Specifically, the temperature of the node is calculated by weighting and averaging the branch return water temperature of the injection node according to the injection flow, and the branch water supply temperature of the outflow node is equal to the temperature of the node. According to the node temperature hybrid constraint, for any node jdn: according to the system parameters, obtaining a thermodynamic temperature mixing equation of a node jdn in the comprehensive energy heat supply system:
wherein, T nod，jdn Is the temperature, g, of node jdn k Flow for branch k, T r，k Return water temperature, R, of branch k jdn ＝{r 1 ，r 2 ，...r m The branch number set injected into the node is obtained; t is s，jdn Temperature of water supply to branch j, S i ＝{s 1 ，s 2 ，...s n And is the number set of branches flowing out of the node.
Writing equation 21 in time-stamped form:
defining a node temperature column vector T nod ＝[T nod，1 ，T nod，2 ，...，T nod，a ] T 。
And S620, obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixing equation of each node.
Specifically, the thermodynamic temperature mixed equations of the nodes are arranged according to the sequence of node numbers, and a system thermodynamic topology constraint equation set is obtained.
The system thermodynamic topological constraint equation set is as follows:
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T nod (t)＝[T 1 (t) T 2 (t) ... T a (t)] T
in the formula I a Is a dimension unit matrix; i is b Is a b-dimensional identity matrix; a. The rn Is a matrix of a x b dimension; a is the number of nodes, b is the number of branches and for i =1,2, a, k =1,2, a s Is the temperature of the supplied water; t is a unit of r The temperature of the return water; a. The sn Is a b × a dimensional matrix and for i =1,2, a, k =1,2.
The step S700 of obtaining a dynamic simulation model of the thermal part includes: and obtaining a dynamic simulation model of the thermodynamic part of the comprehensive energy heating system according to the pipeline branch thermodynamic equation set, the heat user branch thermodynamic equation set, the heat source branch thermodynamic equation set and the system thermodynamic topological constraint equation set.
Specifically, formula 14, formula 17, formula 20 and formula 23 are combined to obtain a dynamic simulation model of the thermal part. The dynamic simulation model of the thermal part is as follows:
I a 、A rn 、A sn 、I b a coefficient matrix which is a system thermodynamic topological constraint equation set; a. The sh0 、A rh0 A coefficient matrix of the heat source branch thermodynamic equation set is obtained; a. The sl0 、A rl0 、A sz0 、A rz0 、A bl 、A bz A coefficient matrix of the thermal user branch thermodynamic equation set; a. The e 、I q0 、A se 、A re And the coefficient matrix of the pipeline branch thermodynamic equation set.
Equation 24 can be abbreviated as:
A t T(t)＝b t equation 25
Wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A rl 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
the simulation initial step S800 includes: setting the starting time of the simulation duration to be t =0, setting initial pressure, initial flow and initial temperature, taking the initial pressure as the current pressure, the initial flow as the current flow, and the initial temperature as the current temperature.
The simulation processing step S900 includes: inputting the current pressure and the current flow into the hydraulic part dynamic simulation model to obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; and inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part.
During the execution of the simulation processing step of this embodiment, a certain parameter of the system parameters may be changed. Optionally, in this case, the embodiment may return to the system parameter input step to modify the corresponding changed parameter. For example, if the power of the circulation pump in the integrated energy heating system is reduced to half of the original power, the pressure increment caused by the power source in the original branch in the system parameter needs to be divided by 2 to obtain the pressure increment caused by the power source in the new branch, and the pressure increment caused by the power source in the new branch is used to replace the pressure increment caused by the power source in the original branch in the system parameter.
The simulation judging step S1000 includes: and after waiting for the simulation step length delta t, judging whether the current time reaches the simulation termination time of the simulation duration, if so, executing the simulation result determining step S1100, and if not, executing the simulation circulating step S1200.
The simulation result determining step S1100 includes: and determining the pressure, flow and temperature obtained by each simulation as a simulation result.
It is to be understood that the simulation result of the present embodiment may be a simulation result obtained from any one of the pressure, the flow rate, and the temperature obtained from each simulation. The representation of the simulation result may be a graph. It is understood that any value obtained by the simulation process may be included in the simulation results. For example, as shown in fig. 4, the simulation result may be a flow trend graph of the branch 36, the branch 5 and the branch 33 during the simulation process. As shown in fig. 5, the simulation result may be a temperature trend graph of the temperatures of the node 2, the node 12, the node 25 and the node 1 during the simulation process. As shown in fig. 6, the simulation result may be a power variation trend graph of the heat source power and the heat load power in the simulation process. As shown in fig. 7, the simulation result may be a graph of the building temperature variation trend of the branch 32, the branch 33 and the branch 34 in the simulation process.
The simulation loop step S1200 includes: and taking the first pressure as the current pressure, the first flow as the current flow, and the first temperature as the current temperature, and returning to execute the simulation processing step.
The invention provides a dynamic simulation method of a comprehensive energy heating system based on an extended node method, which is applied to the comprehensive energy heating system and comprises the following steps: the method comprises the steps of inputting system parameters, obtaining a hydraulic part dynamic simulation model, obtaining a pipeline branch thermodynamic equation set, obtaining a heat user branch thermodynamic equation set, obtaining a heat source branch thermodynamic equation set, obtaining a system thermodynamic topology constraint equation set, obtaining a thermodynamic part dynamic simulation model, performing simulation initial, performing simulation processing, performing simulation judgment, determining a simulation result and performing simulation circulation.
The comprehensive energy heat supply system in the embodiment is divided into a water supply network and a water return network, the water supply network and the water return network are completely symmetrical, corresponding pipeline branch parameters are the same, a heat source branch and a heat user branch are connected with the water supply network and the water return network to form a closed hot working medium flowing loop, and the water supply network and the water return network are composed of pipeline branches and nodes. The heat source and the circulating pump are both positioned in the branch with the branch number of 36, the circulating pump is a power source in the branch and provides power for the flow of the hot working medium, and the power is provided and described by the pressure increment caused by the power source. The topological structures of the water supply net and the water return net are completely symmetrical. The network topology of the integrated energy heating system in this embodiment may be as shown in fig. 2, and the pipe numbers may be as shown in table 1. The system parameters may be as shown in table 2 and the heat source power may be as shown in table 3.
TABLE 1
Note: the return pipes (numbered 14-26) corresponding to the water supply pipes have the same parameters and are not listed.
TABLE 2
| Parameter name | Numerical value |
| Pressure increment caused by No. 36 pipeline power source | 4000MPa |
| Density of hot working medium | 1000kg/m 3 |
| Specific heat capacity of thermal medium | 4186J/(kg·℃) |
| Mass of hot working medium in heat source | 1000kg |
| Efficiency of heat source heating | 90％ |
| Heat radiation area of radiator | 250m 2 |
| Heat dissipation coefficient of heat sink | 127.2J/(kg·m 2 ) |
| Specific heat capacity of air | 1005.6J/(kg·℃) |
| Density of air | 1.239kg/m 3 |
| Heating building volume | 3240m 3 |
| Heat radiation coefficient of heat supply building | 23.2J/(kg·m 3 ) |
| Quality of hot working medium in radiator | 200kg |
| Length of micro-element segment of pipeline | 10m |
| Node equivalent flow capacity | 0.02％ |
| Simulation step length | 5s |
| Simulated termination time | 80000s |
| Number of nodes | 28 |
| Number of branches | 36 |
| Initial value of pressure | 0Mpa |
| Initial value of flow | 0kg/s |
| Initial value of |
20℃ |
TABLE 3
| Time (10) 3 s) | ＜5 | 5-40 | 40-45 | 45-80 |
| Power (kW) | 8×10 3 | 5×10 3 | 1.5×10 4 | 1×10 4 |
Corresponding to the above method embodiment, this embodiment further provides an extended node method-based dynamic simulation apparatus for an integrated energy heating system, which is applied to the integrated energy heating system, and the structure of the apparatus is shown in fig. 8, and the apparatus may include: the system comprises a system parameter input unit 100, a hydraulic part dynamic simulation model obtaining unit 200, a pipeline branch thermodynamic equation set obtaining unit 300, a heat user branch thermodynamic equation set obtaining unit 400, a heat source branch thermodynamic equation set obtaining unit 500, a system thermodynamic topology constraint equation set obtaining unit 600, a thermodynamic part dynamic simulation model obtaining unit 700, a simulation initial unit 800, a simulation processing unit 900, a simulation judging unit 1000, a simulation result determining unit 1100 and a simulation circulating unit 1200.
The system parameter input unit 100 is configured to input system parameters of the integrated energy heating system, where the system parameters include a network topology connection relationship of a heating pipe network of the integrated energy heating system, a network parameter, a thermal working medium parameter, a heat source parameter, a heat user parameter, and a simulation calculation parameter, where the network parameter includes at least one of a branch length, a branch inner wall cross-sectional area, a branch thermal resistance, a branch resistance coefficient, a pressure increment caused by a power source in the branch, and a node equivalent flow capacity, the thermal working medium parameter may include at least one of a thermal working medium density and a thermal working medium specific heat capacity, the heat source parameter may include at least one of a mass of a thermal working medium in the heat source, a heat source heating power, and a heat source heating efficiency, the heat user parameter may include at least one of a radiator area, a radiator heat dissipation coefficient, a mass of a thermal working medium in the radiator, a volume, a heat supply building heat dissipation coefficient, an air density, and an air specific heat capacity, and the simulation calculation parameter may include at least one of a simulation step length, a pipeline infinitesimal length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system.
Optionally, in this embodiment, the process of obtaining the node-branch correlation matrix and the node-branch price reduction correlation matrix of the integrated energy heating system by the system parameter input unit 100 may be as follows:
selecting any node in the pipe network as a pressure reference node, numbering other nodes except the reference node from 1, and recording the total number of the node numbers as a-1; all branches are numbered in sequence from 1, the number of pipeline branches is q, the number of hot user branches is v, the number of hot source branches is f, and the total number of branches b = q + v + f. Wherein: the number of the pipeline branch is from 1 to q; the hot user branch is numbered from q +1 to q + v; the heat source legs are numbered from q + v +1 to b. For example, the number of the pipe branch in this embodiment may be 1 to 13, the number of the hot user branch may be 14 to 26, and the number of the hot source branch may be 27 to 36.
According to the topological connection relation of the pipe network, an a multiplied by b dimension node-branch incidence matrix A can be established: that is, if a branch k flows from the node jdn to the node jdm, a (jdn, k) =1,A (jdm, k) = -1, the remaining elements in the k-th column are 0, and the node with the node number of 0 is used as a reference node, so as to obtain a (a-1) × b-dimensional node-branch reduced price correlation matrix a'.
The type of the branch in the present embodiment may include at least one type of a pipe branch, a hot user branch, and a hot source branch.
The hydraulic part dynamic simulation model obtaining unit 200 is used for obtaining a system node flow conservation equation and a system branch momentum conservation equation according to the system parameters; and according to an implicit Euler method, carrying out differencing processing on the system node flow conservation equation and the system branch momentum conservation linearized equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation, and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system.
Specifically, the hydraulic part dynamic simulation model obtaining unit 200 obtains a system node flow conservation equation according to the system parameters.
Specifically, the system node flow conservation equation is given by formula 1:
wherein, C is a node equivalent flow capacity matrix, the main diagonal element of the node equivalent flow capacity matrix is the equivalent flow capacity corresponding to the node, and the non-main diagonal element is 0.p is a node pressure column vector whose elements are the pressure corresponding to the node relative to a reference point; g is the branch flow column vector. Differentiating the formula 1 by using an implicit Euler method, and obtaining a system node flow conservation discretization equation after sorting:
p (t) is a node pressure column vector at a time t to be solved, g (t) is a branch flow column vector at the time t to be solved, and p (t-delta t) is a node pressure historical quantity which is obtained by calculation at a time before the time t to be solved and is separated by a simulation step length delta t from the time t to be solved.
The hydraulic part dynamic simulation model obtaining unit 200 obtains a system branch momentum conservation equation according to the system parameters.
Specifically, the system branch momentum conservation equation is formula 3:
in the formula, L m The branch inertia coefficient matrix is obtained by taking the main diagonal element as the inertia coefficient of the corresponding branch and taking the non-main diagonal element as 0. Taking the branch k as an example, the inertia coefficient can be determined by the branch length l k And the cross-sectional area A of the inner wall of the branch c，k And calculating to obtain the following formula:
A′ T transposing the node-branch reduced order incidence matrix obtained in the step 1; h is a pressure increment column vector, and the element of the pressure increment column vector is the pressure increment caused by the power source in the corresponding branch; s is a resistance coefficient matrix, the main diagonal element of the matrix is the resistance coefficient of the corresponding branch, and the non-main diagonal element of the matrix is 0; f (g) is a square function of the branch flow and can be written as:
and performing Taylor expansion on the state g (t-delta t) at the time t-delta t by using the nonlinear term f (g) and reserving the first two terms to obtain a linearized expression of the nonlinear function f (g) at the time t-delta t, and further using the linearized expression to approximately calculate the next time step f (g (t)):
f (g (t)) ≈ f (g (t- Δ t)) + J (t- Δ t) [ g (t) -g (t- Δ t) ] equation 6
Wherein J (t-delta t) is a Jacobian matrix of f (g) in the state at the time t-delta t. Differentiating the formula 3 by using an implicit Euler method, substituting the formula 6, and sorting to obtain a system branch momentum conservation discretization equation:
the hydraulic part dynamic simulation model obtaining unit 200 combines the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain the hydraulic part dynamic simulation model of the comprehensive energy heating system.
Specifically, the hydraulic part dynamic simulation model is formula 8:
wherein, the first and the second end of the pipe are connected with each other,
in the formula, H (t) is a column vector formed by pressure increment caused by a power source in a branch, and is a known input quantity, and a is a node-branch incidence matrix of the comprehensive energy heating system: that is, if a branch k flows from the node jdn to the node jdm, a (jdn, k) =1,A (jdm, k) = -1, the remaining elements in the k-th column are 0, and the node with the node number of 0 is used as a reference node, so as to obtain a (a-1) × b-dimensional node-branch reduced price correlation matrix a'; p (t) is a node pressure column vector at time t; g (t) is a branch flow column vector at the time t; s is a resistance coefficient matrix which is a diagonal matrix, and main diagonal elements of the matrix are branch resistance coefficients of corresponding branches; h is a pressure increment column vector, and elements of the column vector are pressure increments caused by power sources in corresponding branches; l is m Main diagonal element ofThe element is the inertia coefficient of the corresponding branch, which can be determined by the branch length l k And the cross-sectional area A of the inner wall of the branch c，k And calculating to obtain the following formula:
L m is a diagonal matrix; the main diagonal element of C is the equivalent flow capacity of the corresponding node and is a diagonal matrix; f (g) is a square function of branch flow, and J (t-delta t) is a Jacobian matrix of f (g) at t-delta t.
The hydraulic part dynamic simulation model obtaining unit 200 solves the formula 8 to obtain pressure column vectors of each node and flow column vectors of each branch of the comprehensive energy heating system at the time t.
The pipeline branch thermodynamic equation set obtaining unit 300 is configured to obtain a pipeline branch thermodynamic equation of each pipeline branch in the integrated energy heating system according to the system parameters; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; and obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch.
Specifically, the pipeline branch thermodynamic equation set obtaining unit 300 obtains the pipeline branch thermodynamic equation of each pipeline branch in the integrated energy heating system according to the system parameters.
For any pipeline branch gd in each pipeline branch: according to system parameters, the pipeline branch thermodynamic equation for pipeline branch k may be:
wherein, T gd Is a pipeline temperature variable.
The pipeline branch thermodynamic equation set obtaining unit 300 performs spatial partial derivative term differentiation on the pipeline branch thermodynamic equations of the pipeline branches according to a infinitesimal method, and performs time partial derivative term differentiation on the pipeline branch thermodynamic equations of the pipeline branches according to an implicit euler method to obtain the pipeline branch discretization thermodynamic equations of the pipeline branches.
For ease of understanding, the conduit branch gd is still used here as an example: to eliminate T in equation 9 gd Spatial partial derivative of
The pipeline branch gd is divided into a plurality of infinitesimal segments. According to the length delta x of a single pipeline infinitesimal and the length l of a pipeline branch gd gd The number of infinitesimal segments can be obtained
Is a pair of
And rounding up. For the wyd infinitesimal sections of the pipeline branch gd, the thermodynamic equation form is formula 10:
the boundary conditions are as follows:
wherein, the first and the second end of the pipe are connected with each other,
the temperature of wyd infinitesimal sections of the pipeline branch gd; a. The c，gd The cross section area of the inner wall of the branch of the pipeline gd; ρ is the density of the hot working medium; c. C p The specific heat capacity of the thermal working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g gd The hot working medium flow rate of the pipeline branch gd; in the dynamic simulation model of the heating power part, a pipeline branch and a hot user branch are required to be connectedAnd the heat source branch is considered respectively, and each pipeline branch gd defines the head end water supply temperature T according to the water flow direction s，gd And tail end return water temperature T r，gd As shown in fig. 3. T is s，gd Supply water temperature, T, to head end of pipeline branch gd r，gd The tail end return water temperature of the pipeline branch gd.
Differentiating equation 10 using the implicit Euler method, g gd (t) is the flow of the pipeline branch gd at the time t, and the discretized thermodynamic dynamic equation of wyd infinitesimal sections can be obtained by sorting:
wherein the content of the first and second substances,
writing the discretization dynamic equations and boundary conditions of all infinitesimal sections of the pipeline branch gd according to the sequence to obtain a discretization thermodynamic equation set of the pipeline branch gd:
wherein the content of the first and second substances,
the column vector is composed of all micro-element section temperatures of the pipeline branch gd and comprises:
the pipeline branch thermal equation set obtaining unit 300 obtains a pipeline branch thermal equation set according to the pipeline branch discretization thermal equation of each pipeline branch.
Specifically, the pipeline branch discretization thermodynamic equations of all pipeline branches in the comprehensive energy heating system are arranged according to the sequence of the serial numbers of the pipeline branches, and a pipeline branch thermodynamic equation set is obtained.
The thermodynamic equation set of the pipeline branch is as follows:
wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
wherein, T e (t) is the temperature column vector of the micro-element section of all the pipeline branches; i is q Is a q-dimensional identity matrix; a. The se And A re Is composed of
A dimension matrix incorporating:
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
wherein q is the number of pipeline branches I q Is a q-dimensional identity matrix; t is e (t) is the temperature column vector of the micro-element section of all the pipeline branches; t is s (t) column vectors of water supply temperatures for all the branches; t is a unit of r (t) is the array vector of return water temperature of all branches; a. The se And A re Is composed of
A dimension matrix; u. u gd Number of micro-elements, T, of a branch of a pipeline s，k The temperature of the water supply for branch k; t is a unit of r，k The return water temperature of the branch k; t is e Is the temperature of the infinitesimal section; Δ x is the length of the pipeline infinitesimal section; delta t is a simulation step length; a. The c，gd The sectional area of the inner wall of the pipeline branch gd; rho is the density of the hot working medium; c. C p The specific heat capacity of the thermal working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g gd The hot working medium flow rate of the pipeline branch gd; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The thermal user branch thermodynamic equation set obtaining unit 400 is configured to obtain a thermal user branch thermodynamic equation of each thermal user branch in the integrated energy heating system based on a thermal capacity particle method according to the system parameter; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; and obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch.
Specifically, the thermal user branch thermodynamic equation set obtaining unit 400 obtains the thermal user branch thermodynamic equation of each thermal user branch in the integrated energy heating system based on the thermal capacity particle method according to the system parameters.
For any hot user branch ryh in each hot user branch: according to system parameters, a heat balance equation of a hot user branch rh is obtained based on a heat capacity particle method and is as follows:
in the formula, m h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; t is s，ryh The water supply temperature of the hot user branch ryh; t is r，ryh The return water temperature of the hot user branch ryh; for each branch of heat user, the corresponding temperature of heat supply building can be defined, and the temperature vector T of heat supply building of heat user can be defined b ＝[T b，1 ，T b，2 ，...，T b，v ]The corresponding hot subscriber branch is numbered from q +1 to q + v, T b，ryh-q The temperature of the heating building for the hot user branch ryh; q net Supplying heat to a heat supply network; q rad Heat dissipation for the radiator to the heating building; g ryh The hot working medium flow rate of a hot user branch ryh; k r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; rho a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is o Is the outdoor temperature; q heat Heat dissipation capacity to the outside for the heat supply building; q. q.s V The heat dissipation coefficient of the building is the heat dissipation coefficient of the heat supply.
The hot subscriber branch thermal equation set obtaining unit 400 performs differentiation processing on the hot subscriber branch thermal equations of each hot subscriber branch according to the implicit euler method to obtain the hot subscriber branch discretization thermal equations of each hot subscriber branch.
For ease of understanding, the hot subscriber leg ryh is still used as an example for explanation here:
differentiating the formula 15 using the implicit Euler method, and differentiating the one in the formula 15Substituting equations in lines 3, 4 and 5 into equations in lines 1 and 2 to eliminate Q net 、Q rad And Q heat And after sorting, obtaining a discretization thermodynamic equation of the hot user branch ryh:
in the formula (I), the compound is shown in the specification,
the thermal user branch thermodynamic equation set obtaining unit 400 obtains a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch.
Specifically, the discretization thermodynamic equations of the thermal user branches are arranged according to the serial number of the thermal user branches, and a thermodynamic equation set of the thermal user branches is obtained.
The thermal user branch thermodynamic equation system is as follows:
wherein the content of the first and second substances,
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A rl 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T b (t)＝[T b，1 (t) T b，2 (t) ... T b， v(t)] T
wherein A is sl 、A rl 、A bl 、A sz 、A rz 、A bz Is a v-dimensional diagonal matrix, i.e. all non-principal diagonal elements are 0, the principal diagonal element pair i =1,2.
In the formula, T s Is a water supply temperature column vector; t is r Is a backwater temperature column vector; t is b Is a heating building temperature column vector; m is h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; g i+p Numbering the hot working medium flow of the hot user branch with the number of i + p for the hot user branch; k r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; delta t is a simulation step length; ρ is a unit of a gradient a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is a unit of o Is the outdoor temperature；q V Is a heat dissipation coefficient; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The heat source branch thermodynamic equation set obtaining unit 500 is configured to obtain a heat source branch thermodynamic equation of each heat source branch in the integrated energy heating system based on a thermal capacity particle method according to the system parameters; carrying out differential processing on the heat source branch thermodynamic equations of the heat source branches according to an implicit Euler method to obtain heat source branch discretization thermodynamic equations of the heat source branches; and obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch.
Specifically, the heat source branch thermodynamic equation set obtaining unit 500 obtains the heat source branch thermodynamic equation of each heat source branch in the integrated energy heating system based on the thermal capacity particle method according to the system parameters.
Specifically, the heat source is regarded as a heat capacity temperature particle, and the obtained net heat is regarded as the difference between the heat brought by the heating power and the heat taken away by the supply water and the return water. For any heat source branch ry in each heat source branch: according to system parameters, a thermal capacity particle method is used for obtaining a thermal equation of a heat source branch k, wherein the thermal equation of the heat source branch k is as follows:
in the formula, m w Is the mass of hot working medium in the heat source; c. C p Is the specific heat capacity of the hot working medium; t is a unit of s，ry The temperature of the water supply for the heat source branch ry; t is r，ry The return water temperature of the heat source branch ry; p is W And eta are respectively the heating power and heating efficiency of the heat source, g ry The flow of the hot working medium in the heat source of the heat source branch ry.
The heat source branch thermal equation set obtaining unit 500 performs differentiation processing on the heat source branch thermal equations of the heat source branches according to the implicit euler method to obtain the heat source branch discretization thermal equations of the heat source branches.
For the sake of understanding, the heat source branch k is still taken as an example for explanation:
differentiating the formula 18 by using an implicit Euler method, and sorting to obtain:
wherein, K sh ＝(-g ry (t)c p )，
The heat source branch thermal equation set obtaining unit 500 obtains a heat source branch thermal equation set according to the heat source branch discretization thermal equation of each heat source branch.
Specifically, the discretization thermodynamic equations of the heat source branches of each heat source branch are arranged according to the serial number sequence of the heat source branches to obtain a thermodynamic equation set of the heat source branches.
The thermodynamic equation set of the heat source branch is as follows:
wherein the content of the first and second substances,
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
A sh0 、A rh0 is an f-dimensional diagonal matrix with 0 as the non-main diagonal element, i =1,2 as the main diagonal element pair, f has:
in the formula, delta t is a simulation step length; m is w Is the mass of hot working medium in the heat source; c. C p The specific heat capacity of the thermal working medium; t is a unit of s，k The water supply temperature of the heat source branch k; t is a unit of r The temperature of the return water is; p W And eta are the heating power and the heating efficiency of the heat source respectively; g is a radical of formula i+p+v The hot working medium flow of the heat source branch with the serial number of i + p + v is numbered for the heat source branch; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches.
The system thermodynamic topological constraint equation set obtaining unit 600 is configured to obtain a thermodynamic temperature mixing equation of each node in the integrated energy heating system according to the system parameters; and obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node.
Specifically, the system thermodynamic topological constraint equation set obtaining unit 600 obtains a thermodynamic temperature mixing equation of each node in the integrated energy heating system according to the system parameters.
Specifically, the temperature of the node is calculated by weighting and averaging the branch return water temperature of the injection node according to the injection flow, and the branch water supply temperature of the outflow node is equal to the temperature of the node. According to the node temperature mixing constraint, for any node jdn in each node: according to the system parameters, obtaining a thermodynamic temperature mixing equation of a node jdn in the comprehensive energy heat supply system:
wherein, T nod，jdn Is the temperature, g, of node jdn k Flow for branch k, T r，k Return water temperature, R, of branch k jdn ＝{r 1 ，r 2 ，...r m The branch number set injected into the node is obtained; t is s，jdn Temperature of water supply, S, for branch j i ＝{s 1 ，s 2 ，...s n And is the number set of branches flowing out of the node.
Writing equation 21 in time-stamped form:
defining a node temperature column vector T nod ＝[T nod，1 ，T nod，2 ，...，T nod，a ] T 。
The system thermodynamic topology constraint equation set obtaining unit 600 obtains a system thermodynamic topology constraint equation set according to the thermodynamic temperature mixing equation of each node.
Specifically, the thermodynamic temperature mixed equations of the nodes are arranged according to the serial number of the nodes, and a system thermodynamic topology constraint equation set is obtained.
The system thermodynamic topological constraint equation set is as follows:
T s (t)＝[T s，1 (t) T s，2 (t) ... T s，b (t)] T
T r (t)＝[T r，1 (t) T r，2 (t) ... T r，b (t)] T
T nod (t)＝[T 1 (t) T 2 (t) ... T a (t)] T
in the formula I a Is a dimension unit matrix; I.C. A b Is a b-dimensional identity matrix; a. The rn Is a x b dimensionA matrix; a is the number of nodes, b is the number of branches and for i =1,2, a, k =1,2, a s Is the temperature of the supplied water; t is a unit of r The temperature of the return water is; a. The sn Is a b × a dimensional matrix and for i =1,2, a, k =1,2.
The thermodynamic part dynamic simulation model obtaining unit 700 is configured to obtain a thermodynamic part dynamic simulation model of the integrated energy heating system according to the pipeline branch thermodynamic equation set, the heat user branch thermodynamic equation set, the heat source branch thermodynamic equation set, and the system thermodynamic topology constraint equation set.
Specifically, the thermal part dynamic simulation model obtaining unit 700 combines formula 14, formula 17, formula 20, and formula 23 to obtain the thermal part dynamic simulation model. The dynamic simulation model of the thermal part is as follows:
I a 、A rn 、A sn 、I b a coefficient matrix which is a system thermodynamic topological constraint equation set; a. The sh0 、A rh0 A coefficient matrix of the heat source branch thermodynamic equation set; a. The sl0 、A rl0 、A sz0 、A rz0 、A bl 、A bz A coefficient matrix of the thermal user branch thermodynamic equation set; a. The e 、I q0 、A se 、A re And the coefficient matrix of the pipeline branch thermodynamic equation set.
Equation 24 can be abbreviated as:
A t T(t)＝b t equation 25
Wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A r1 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
the simulation initialization unit 800 is configured to set a starting time of the simulation duration to be t =0, set an initial pressure, an initial flow and an initial temperature, use the initial pressure as a current pressure, use the initial flow as a current flow, and use the initial temperature as a current temperature.
The simulation processing unit 900 is configured to input the current pressure and the current flow into the hydraulic part dynamic simulation model, and obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; and inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part.
During the execution of the simulation processing step of this embodiment, a certain parameter of the system parameters may be changed. Optionally, in this case, the simulation processing unit 900 of the present embodiment may trigger the system parameter input unit 100 to modify the corresponding changed parameter. For example, if the power of the circulation pump in the integrated energy heating system is reduced to half of the original power, the pressure increment caused by the power source in the original branch in the system parameters needs to be divided by 2 to obtain the pressure increment caused by the power source in the new branch, and the pressure increment caused by the power source in the new branch is used to replace the pressure increment caused by the power source in the original branch in the system parameters.
The simulation judging unit 1000 is configured to judge whether the current time reaches the simulation termination time of the simulation duration after waiting for the simulation step Δ t, and if so, trigger the simulation result determining unit 1100, and if not, trigger the simulation circulating unit 1200.
The simulation result determining unit 1100 is configured to determine the pressure, flow and temperature obtained by each simulation as a simulation result.
It is to be understood that the simulation result of the present embodiment may be a simulation result obtained from any one of the pressure, the flow rate, and the temperature obtained from each simulation. The representation of the simulation result may be a graph. It is understood that any value obtained by the simulation process may be included in the simulation results. For example, as shown in fig. 4, the simulation result may be a flow trend graph of the branch 36, the branch 5 and the branch 33 during the simulation process. As shown in fig. 5, the simulation result may be a temperature trend graph of the temperatures of the node 2, the node 12, the node 25 and the node 1 during the simulation process. As shown in fig. 6, the simulation result may be a power variation trend graph of the heat source power and the heat load power in the simulation process. As shown in fig. 7, the simulation result may be a building temperature trend graph of the branch 32, the branch 33, and the branch 34 during the simulation process.
The simulation circulation unit 1200 is configured to trigger the simulation processing unit 900 by using the first pressure as a current pressure, using the first flow as a current flow, and using the first temperature as a current temperature.
The invention provides a dynamic simulation device of a comprehensive energy heating system based on an extended node method, which is applied to the comprehensive energy heating system, and comprises the following components: the system comprises a system parameter input unit 100, a hydraulic part dynamic simulation model obtaining unit 200, a pipeline branch thermodynamic equation set obtaining unit 300, a heat user branch thermodynamic equation set obtaining unit 400, a heat source branch thermodynamic equation set obtaining unit 500, a system thermodynamic topology constraint equation set obtaining unit 600, a thermodynamic part dynamic simulation model obtaining unit 700, a simulation initial unit 800, a simulation processing unit 900, a simulation judging unit 1000, a simulation result determining unit 1100 and a simulation circulating unit 1200.
The comprehensive energy heat supply system in the embodiment is divided into a water supply network and a water return network, the water supply network and the water return network are completely symmetrical, corresponding pipeline branch parameters are the same, a heat source branch and a heat user branch are connected with the water supply network and the water return network to form a closed hot working medium flowing loop, and the water supply network and the water return network are composed of pipeline branches and nodes. The heat source and the circulating pump are both positioned in the branch with the branch number of 36, the circulating pump is a power source in the branch and provides power for the flow of the hot working medium, and the size of the provided power is described by the pressure increment caused by the power source. The topological structures of the water supply net and the water return net are completely symmetrical.
In this application, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. The term "comprising", without further limitation, means that the element so defined is not excluded from the group consisting of additional identical elements in the process, method, article, or apparatus that comprises the element.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (8)
1. A dynamic simulation method of an integrated energy heating system based on an extended node method is characterized by being applied to the integrated energy heating system, and comprises the following steps: inputting system parameters, obtaining a hydraulic part dynamic simulation model, obtaining a pipeline branch thermodynamic equation set, obtaining a heat user branch thermodynamic equation set, obtaining a heat source branch thermodynamic equation set, obtaining a system thermodynamic topology constraint equation set, obtaining a thermodynamic part dynamic simulation model, simulating an initial step, a simulation processing step, a simulation judging step, a simulation result determining step and a simulation circulating step,
the system parameter inputting step comprises: inputting system parameters of the comprehensive energy heat supply system, wherein the system parameters comprise at least one of a branch length, a branch inner wall sectional area, a branch thermal resistance, a branch resistance coefficient, a pressure increment caused by a power source in a branch and a node equivalent flow capacity, the thermal working medium parameters comprise at least one of a thermal working medium density and a thermal working medium specific heat capacity, the heat source parameters comprise at least one of a mass of a heat source internal heat working medium, a heat source heating power and a heat source heating efficiency, the heat user parameters comprise at least one of a radiator area, a radiator heat dissipation coefficient, a mass and a volume of the heat working medium in the radiator, a heat supply building heat dissipation coefficient, an air density and an air specific heat capacity, and the simulation calculation parameters comprise at least one of a simulation duration, a simulation step length and a pipeline infinitesimal length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system;
the step of obtaining the dynamic simulation model of the hydraulic part comprises the following steps: according to the system parameters, a system node flow conservation equation and a system branch momentum conservation equation are obtained; according to an implicit Euler method, carrying out differential processing on the system node flow conservation equation and the system branch momentum conservation equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation, and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system;
the step of obtaining the pipeline branch thermodynamic equation set comprises the following steps: according to the system parameters, acquiring a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch;
the step of obtaining the thermal equation set of the hot user branch comprises the following steps: according to the system parameters, obtaining thermal user branch thermodynamic equations of all thermal user branches in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch;
the step of obtaining the heat source branch thermodynamic equation set comprises the following steps: according to the system parameters, obtaining a heat source branch thermodynamic equation of each heat source branch in the comprehensive energy heating system based on a thermal capacity particle method; carrying out differential processing on the heat source branch thermodynamic equations of the heat source branches according to an implicit Euler method to obtain heat source branch discretization thermodynamic equations of the heat source branches; obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch;
the step of obtaining the system thermodynamic topological constraint equation system comprises the following steps: according to the system parameters, acquiring a thermodynamic temperature mixing equation of each node in the comprehensive energy heating system; obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node;
the step of obtaining the dynamic simulation model of the thermal part comprises the following steps: obtaining a dynamic simulation model of a thermal part of the comprehensive energy heating system according to the pipeline branch thermal equation set, the heat user branch thermal equation set, the heat source branch thermal equation set and the system thermal topological constraint equation set;
the simulation initial step comprises: setting the starting time of the simulation duration to be t =0, setting initial pressure, initial flow and initial temperature, taking the initial pressure as the current pressure, the initial flow as the current flow, and the initial temperature as the current temperature;
the simulation processing step comprises: inputting the current pressure and the current flow into the hydraulic part dynamic simulation model to obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part;
the simulation judging step comprises: after waiting for the simulation step length delta t, judging whether the current time reaches the simulation termination time of the simulation duration, if so, executing the simulation result determining step, and if not, executing the simulation circulating step;
the simulation result determining step comprises: determining the pressure, flow and temperature obtained by each simulation as a simulation result;
the simulation loop step includes: and taking the first pressure as the current pressure, the first flow as the current flow, and the first temperature as the current temperature, and returning to execute the simulation processing step.
2. The method of claim 1, wherein the hydraulic segment dynamic simulation model is:
wherein the content of the first and second substances,
in the formula, H (t) is a column vector formed by pressure increment caused by a power source in a branch, and is a known input quantity, and a is a node-branch incidence matrix of the comprehensive energy heating system: that is, if the branch k flows from the node i to the node j, a (i, k) =1,A (j, k) = -1, the remaining elements in the k-th column are 0, and the node with the node number of 0 is used as a reference node to obtain a (a-1) × b-dimensional node-branch reduced price incidence matrix a'; p (t) is a node pressure column vector at time t; g (t) is a branch flow column vector at the moment t; s is a resistance coefficient matrix which is a diagonal matrix, and main diagonal elements of the matrix are branch resistance coefficients of corresponding branches; h is a pressure increment column vector, and the element of the pressure increment column vector is the pressure increment caused by the power source in the corresponding branch; l is m The main diagonal element of (A) is the inertia coefficient of the corresponding branch, which can be determined by the branch length l k And the cross-sectional area A of the inner wall of the branch c,k And calculating to obtain the following formula:
L m is a diagonal matrix; the main diagonal element of C is the equivalent flow capacity of the corresponding node and is a diagonal matrix; f (g) is the square function of the branch flow, J (t)- Δ t) is the Jacobian matrix of f (g) at time t- Δ t; a. The T Is a transposed matrix of the A matrix; delta t is a simulation step length; p (t-delta t) is a node pressure column vector at the time t-delta t; g (t-delta t) is a branch flow column vector at the time of t-delta t; f (g (t- Δ t)) is a square function of the branch flow at time t- Δ t.
3. The method of claim 1, wherein the system of pipeline branch thermodynamic equations is:
wherein the content of the first and second substances,
I q0 ＝[I q 0 q×v 0 q×f ]
T s (t)＝[T s,1 (t) T s,2 (t) … T s,b (t)] T
T r (t)＝[T r,1 (t) T r,2 (t) … T r,b (t)] T
wherein q is the number of pipeline branches I q Is a q-dimensional identity matrix; t is the time t; t is e (t) is the temperature column vector of the micro-element section of all the pipeline branches at the moment t; t is a unit of s (t) all branch water supply temperature column vectors at the moment t; t is r (t) all branch return water temperature column vectors at the moment t; a. The se And A re Is composed of
A dimension matrix; u. of gd Number of micro-elements, T, of a branch of a pipeline s,k The temperature of the water supply for branch k; t is r,k The return water temperature of the branch k; t is e,i Is the temperature of the infinitesimal section i; Δ x is the length of the pipeline infinitesimal section; delta t is a simulation step length; a. The c,gd The sectional area of the inner wall of the pipeline branch gd; rho is the density of the hot working medium; c. C p The specific heat capacity of the thermal working medium;
is the reciprocal of the thermal resistance of the pipeline branch gd; t is o Is ambient temperature; g is a radical of formula gd (t) is the hot working medium flow of the pipeline branch gd at the time t; f is the number of heat source branches; v is the number of hot subscriber legs; q is the number of pipeline branches; a is the number of nodes; b is the number of branches; 0 q×v A zero matrix of dimension qxv; 0 q×f A zero matrix with dimension of q multiplied by f; wyd is the number of infiniband segments.
4. The method of claim 1, wherein the thermal user leg thermodynamic equation set is:
wherein the content of the first and second substances,
T s (t)＝[T s,1 (t) T s,2 (t) … T s,b (t)] T
T r (t)＝[T r,1 (t) T r,2 (t) … T r,b (t)] T
T b (t)＝[T b,1 (t) T b,2 (t) … T b,v (t)] T
A sl0 ＝[0 v×q A sl 0 v×f ]，A rl0 ＝[0 v×q A rl 0 v×f ]
A sz0 ＝[0 v×q A sz 0 v×f ]，A rz0 ＝[0 v×q A rz 0 v×f ]
in the formula, A sl 、A rl 、A bl 、A sz 、A rz 、A bz All are v-dimensional diagonal matrices, and the main diagonal element pair i =1,2, …, v has:
in the formula, t is t time; t is a unit of s Is a water supply temperature column vector; t is r Is a backwater temperature column vector; t is b A heating building temperature column vector; t is b,i Building temperature column vector T for heat supply b The ith element of (1); m is h Mass of hot working medium in the radiator; c. C p The specific heat capacity of the thermal working medium; g i+p (t) the hot working medium flow of the hot user branch with the hot user branch number of i + p at the moment t; k r And A r The heat dissipation coefficient and the heat dissipation area of the radiator are respectively; delta t is a simulation step length; rho a And c a Density and specific heat capacity of air, respectively; v is the heat supply building volume; t is o Is the outdoor temperature; q. q.s V Is the heat dissipation coefficient; f is the number of heat source branches; v is the number of hot subscriber legs; q is the number of pipeline branches; a is the number of nodes; b is the number of branches; 0 v×q A zero matrix with v × q dimensions; 0 v×f Is a zero matrix with dimension v x f.
5. The method of claim 1, wherein the set of heat source branch thermal equations is:
wherein, the first and the second end of the pipe are connected with each other,
A sh0 ＝[0 f×q 0 f×v A sh ]，A rh0 ＝[0 f×q 0 f×v A rh ]
T s (t)＝[T s,1 (t) T s,2 (t) … T s,b (t)] T
T r (t)＝[T r,1 (t) T r,2 (t) … T r,b (t)] T
A sh 、A rh for an f-dimensional diagonal matrix, the main diagonal element pair i =1,2, …, f, has:
in the formula, delta t is a simulation step length; m is a unit of w Is the mass of hot working medium in the heat source; c. C p The specific heat capacity of the thermal working medium; t is s (t) the water supply temperature at time t; t is a unit of r (t) is the backwater temperature at the moment t; p w (t) the heating power of the heat source at time t; eta is the heating efficiency of the heat source; g i+p+v (t) the thermal working medium flow of the heat source branch with the serial number of i + p + v at the time t; f is the number of heat source branches; v is the number of hot user branches; q is the number of pipeline branches; a is the number of nodes; b is the number of branches; 0 f×q Is a zero matrix with f multiplied by q dimensions; 0 f×v A zero matrix with f x v dimensions; t is r,i+p+v And (t-delta t) is the return water temperature of the heat source branch with the heat source branch number of i + p + v at the time of t-delta t.
6. The method of claim 1, wherein the system thermodynamic topology constraint equation set is:
T s (t)＝[T s,1 (t) T s,2 (t) … T s,b (t)] T
T r (t)＝[T r,1 (t) T r,2 (t) … T r,b (t)] T
T nod (t)＝[T 1 (t) T 2 (t) … T a (t)] T
in the formula I a Is a dimension unit matrix; i is b Is a b-dimensional identity matrix; a. The rn Is a matrix of a x b dimension; a is the number of nodes, b is the number of branches and the pairs i =1,2, …, a, k =1,2, …, b, T s Is the temperature of the supplied water; t is r The temperature of the return water is; a. The sn Is a b x a dimensional matrix and is for i =1,2, …, a, k =1,2, …, b.
7. The method of claim 1, wherein the thermodynamic part dynamic simulation model is:
I a 、A rn 、A sn 、I b a coefficient matrix which is a system thermodynamic topological constraint equation set; a. The sh0 、A rh0 A coefficient matrix of the heat source branch thermodynamic equation set is obtained; a. The sl0 、A rl0 、A sz0 、A rz0 、A bl 、A bz A coefficient matrix of the thermal user branch thermodynamic equation set; a. The e 、I q0 、A se 、A re A coefficient matrix of the pipeline branch thermodynamic equation set; t is nod (t) is a node temperature column vector at time t; t is a unit of s (t) is a water supply temperature column vector at time t; t is r (t) is a backwater temperature column vector at the moment t; t is e (t) the temperature column vectors of the micro-element sections of all the pipeline branches at the moment t; b is a mixture of h The deviation value of the heat source branch thermodynamic equation set is obtained; b 1 And b z By-passing heat for hot usersDeviation value of force equation set.
8. The utility model provides a comprehensive energy heating system dynamic simulation device based on extended node method which characterized in that is applied to comprehensive energy heating system, the device includes: a system parameter input unit, a hydraulic part dynamic simulation model obtaining unit, a pipeline branch thermodynamic equation set obtaining unit, a heat user branch thermodynamic equation set obtaining unit, a heat source branch thermodynamic equation set obtaining unit, a system thermodynamic topology constraint equation set obtaining unit, a thermodynamic part dynamic simulation model obtaining unit, a simulation initial unit, a simulation processing unit, a simulation judging unit, a simulation result determining unit and a simulation circulating unit,
the system parameter input unit is used for inputting system parameters of the comprehensive energy heat supply system, wherein the system parameters comprise a network topology connection relation of a heat supply pipe network of the comprehensive energy heat supply system, network parameters, thermal working medium parameters, heat source parameters, heat user parameters and simulation calculation parameters, the network parameters comprise at least one of branch length, a branch inner wall sectional area, branch thermal resistance, branch resistance coefficients, pressure increment caused by a power source in a branch and node equivalent flow capacity, the thermal working medium parameters comprise at least one of thermal working medium density and thermal working medium specific heat capacity, the heat source parameters comprise at least one of heat source internal heat working medium quality, heat source heating power and heat source heating efficiency, the heat user parameters comprise at least one of heat radiator area, heat radiator heat radiation coefficient, heat working medium quality and volume in a heat radiator, heat supply building heat radiation coefficient, air density and air capacity, and the simulation calculation parameters comprise at least one of simulation duration, simulation step length and pipeline infinitesimal length; the network topology connection relation is a node-branch incidence matrix and a node-branch price reduction incidence matrix of the comprehensive energy heating system;
the hydraulic part dynamic simulation model obtaining unit is used for obtaining a system node flow conservation equation and a system branch momentum conservation equation according to the system parameters; according to an implicit Euler method, carrying out differential processing on the system node flow conservation equation and the system branch momentum conservation equation to obtain a system node flow conservation discretization equation and a system branch momentum conservation discretization equation, and combining the system node flow conservation discretization equation and the system branch momentum conservation discretization equation to obtain a hydraulic part dynamic simulation model of the comprehensive energy heating system;
the pipeline branch thermodynamic equation set obtaining unit is used for obtaining a pipeline branch thermodynamic equation of each pipeline branch in the comprehensive energy heating system according to the system parameters; carrying out space partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to a infinitesimal method, and carrying out time partial derivative term differentiation processing on the pipeline branch thermodynamic equation of each pipeline branch according to an implicit Euler method to obtain a pipeline branch discretization thermodynamic equation of each pipeline branch; obtaining a pipeline branch thermodynamic equation set according to the pipeline branch discretization thermodynamic equation of each pipeline branch;
the system comprises a heat consumer branch thermodynamic equation set obtaining unit, a heat capacity particle method and a heat capacity particle method, wherein the heat consumer branch thermodynamic equation set obtaining unit is used for obtaining heat consumer branch thermodynamic equations of all heat consumer branches in the comprehensive energy heating system according to the system parameters; carrying out differential processing on the thermal equation of each thermal user branch according to an implicit Euler method to obtain the discretization thermal equation of each thermal user branch; obtaining a thermal user branch thermodynamic equation set according to the thermal user branch discretization thermodynamic equation of each thermal user branch;
the heat source branch thermodynamic equation set obtaining unit is used for obtaining heat source branch thermodynamic equations of all heat source branches in the comprehensive energy heating system based on a thermal capacity particle method according to the system parameters; carrying out differential processing on the heat source branch thermodynamic equations of the heat source branches according to an implicit Euler method to obtain heat source branch discretization thermodynamic equations of the heat source branches; obtaining a heat source branch thermodynamic equation set according to the heat source branch discretization thermodynamic equation of each heat source branch;
the system thermodynamic topological constraint equation set obtaining unit is used for obtaining a thermodynamic temperature mixed equation of each node in the comprehensive energy source heat supply system according to the system parameters; obtaining a system thermodynamic topological constraint equation set according to the thermodynamic temperature mixed equation of each node;
the thermodynamic part dynamic simulation model obtaining unit is used for obtaining a thermodynamic part dynamic simulation model of the comprehensive energy heating system according to the pipeline branch thermodynamic equation set, the heat user branch thermodynamic equation set, the heat source branch thermodynamic equation set and the system thermodynamic topology constraint equation set;
the simulation initialization unit is configured to set a starting time of the simulation duration to be t =0, set an initial pressure, an initial flow and an initial temperature, use the initial pressure as a current pressure, use the initial flow as a current flow, and use the initial temperature as a current temperature;
the simulation processing unit is used for inputting the current pressure and the current flow into the hydraulic part dynamic simulation model to obtain a first pressure and a first flow output by the hydraulic part dynamic simulation model; inputting the first flow and the current temperature into the dynamic simulation model of the thermal part to obtain a first temperature output by the dynamic simulation model of the thermal part;
the simulation judging unit is used for judging whether the current time reaches the simulation termination time of the simulation duration after waiting for the simulation step length delta t, if so, triggering the simulation result determining unit, and if not, triggering the simulation circulating unit;
the simulation result determining unit is used for determining the pressure, flow and temperature obtained by each simulation as a simulation result;
and the simulation circulating unit is used for triggering the simulation processing unit by taking the first pressure as the current pressure, taking the first flow as the current flow and taking the first temperature as the current temperature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910816039.3A CN110555264B (en) | 2019-08-30 | 2019-08-30 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910816039.3A CN110555264B (en) | 2019-08-30 | 2019-08-30 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110555264A CN110555264A (en) | 2019-12-10 |
| CN110555264B true CN110555264B (en) | 2023-02-28 |
Family
ID=68738490
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910816039.3A Active CN110555264B (en) | 2019-08-30 | 2019-08-30 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110555264B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111046594B (en) * | 2020-01-09 | 2023-05-02 | 东南大学 | Hot water heating network dynamic simulation method based on cross iteration principle |
| CN111222213B (en) * | 2020-01-15 | 2021-08-03 | 许继集团有限公司 | Dynamic simulation method and device for thermodynamic network |
| CN111310343B (en) * | 2020-02-22 | 2021-10-15 | 清华大学 | Heat supply network hot circuit modeling method for comprehensive energy system scheduling |
| CN112163329A (en) * | 2020-09-21 | 2021-01-01 | 清华大学 | Thermodynamic system matrix model construction method and electronic equipment |
| CN112182905B (en) * | 2020-10-16 | 2023-08-29 | 北京科东电力控制系统有限责任公司 | Heat supply pipe network simulation method and device for comprehensive energy system |
| CN113221428B (en) * | 2021-06-01 | 2022-05-31 | 昆明理工大学 | Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system |
| CN114114958B (en) * | 2021-11-17 | 2023-09-29 | 天津大学 | Discrete simulation method for state of heating system by considering control characteristics of source load temperature |
| CN114004110B (en) * | 2021-11-24 | 2024-04-12 | 天津大学 | Quantized event driven simulation method for electric-thermal integrated energy system |
| CN117332725A (en) * | 2023-12-01 | 2024-01-02 | 中国科学院广州能源研究所 | Steam network dynamic calculation method and device of comprehensive energy system |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
| CN108920866A (en) * | 2018-07-20 | 2018-11-30 | 清华大学 | Heat supply network dynamic regulation operating parameter estimation method based on moving horizon estimation theory |
| CN109255194A (en) * | 2018-09-27 | 2019-01-22 | 河海大学 | A kind of electric-thermal interconnection integrated energy system photovoltaic consumption analysis method containing P2H |
| CN109255466A (en) * | 2018-07-20 | 2019-01-22 | 清华大学 | A kind of heat supply network steady-state operation method for parameter estimation measured based on multi-state |
| CN109284939A (en) * | 2018-10-26 | 2019-01-29 | 南方电网科学研究院有限责任公司 | Combined heat and power random production analog method, device and the equipment of integrated energy system |
| CN109344436A (en) * | 2018-08-28 | 2019-02-15 | 中国石油化工股份有限公司天然气分公司 | A kind of large complicated gas distributing system system in-circuit emulation method |
| CN109376428A (en) * | 2018-10-24 | 2019-02-22 | 南方电网科学研究院有限责任公司 | Reliability estimation method, device, equipment and the storage medium of integrated energy system |
| CN109978625A (en) * | 2019-03-28 | 2019-07-05 | 河海大学 | It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method |
-
2019
- 2019-08-30 CN CN201910816039.3A patent/CN110555264B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
| CN108920866A (en) * | 2018-07-20 | 2018-11-30 | 清华大学 | Heat supply network dynamic regulation operating parameter estimation method based on moving horizon estimation theory |
| CN109255466A (en) * | 2018-07-20 | 2019-01-22 | 清华大学 | A kind of heat supply network steady-state operation method for parameter estimation measured based on multi-state |
| CN109344436A (en) * | 2018-08-28 | 2019-02-15 | 中国石油化工股份有限公司天然气分公司 | A kind of large complicated gas distributing system system in-circuit emulation method |
| CN109255194A (en) * | 2018-09-27 | 2019-01-22 | 河海大学 | A kind of electric-thermal interconnection integrated energy system photovoltaic consumption analysis method containing P2H |
| CN109376428A (en) * | 2018-10-24 | 2019-02-22 | 南方电网科学研究院有限责任公司 | Reliability estimation method, device, equipment and the storage medium of integrated energy system |
| CN109284939A (en) * | 2018-10-26 | 2019-01-29 | 南方电网科学研究院有限责任公司 | Combined heat and power random production analog method, device and the equipment of integrated energy system |
| CN109978625A (en) * | 2019-03-28 | 2019-07-05 | 河海大学 | It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method |
Non-Patent Citations (1)
| Title |
|---|
| 基于键合图的区域综合能源系统动态仿真研究;葛海麟等;《陕西电力》;20170220(第02期);第43-48页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110555264A (en) | 2019-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110555264B (en) | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method | |
| CN109242365A (en) | A kind of Interval Power Flow calculation method of electric-thermal interconnection integrated energy system | |
| CN106647329B (en) | Equivalent modeling method for heat supply pipe network | |
| CN103995548A (en) | Indoor thermal environment control method based on model reduction and multi-model prediction control | |
| CN107947182B (en) | Dynamic power flow analysis method and dynamic power flow analysis system | |
| Cimmino | Semi-Analytical Method for g-Function Calculation of bore fields with series-and parallel-connected boreholes | |
| CN103605837A (en) | Multi-heat-source looped pipe network accessibility analysis and determination method, and multi-heat-source looped pipe network accessibility improvement method | |
| Lan et al. | Modeling of the enthalpy transfer using electric circuit equivalents: Theory and application to transients of multi-carrier energy systems | |
| Wang et al. | Entransy definition and its balance equation for heat transfer with vaporization processes | |
| Hirsch et al. | An efficient numerical solution method for detailed modelling of large 5th generation district heating and cooling networks | |
| CN110309591A (en) | It exchanges heat under a kind of flat finned heat exchanger air side laminar condition and drag computation method | |
| Capone et al. | Accounting for pipeline thermal capacity in district heating simulations | |
| Castorani et al. | Determination of the optimal configuration of energy recovery ventilator through virtual prototyping and DoE techniques | |
| CN112362096B (en) | Method and device for monitoring running state of multi-energy flow, terminal equipment and storage medium | |
| Bäumelt et al. | Distributed agent-based building grey-box model identification | |
| CN113111515B (en) | Unified modeling method of comprehensive energy system | |
| CN105160116B (en) | The uniform Optimization Design of heat exchanger liquid separation | |
| Li et al. | Study of distributed model predictive control for radiant floor heating systems with hydraulic coupling | |
| Quintã et al. | Virtual test bench for the design of control strategies for water heaters | |
| CN113111555B (en) | Quality control thermodynamic system energy flow rapid calculation method based on superposition decoupling method | |
| CN110705066B (en) | Dynamic simulation method for gas-electricity coupling park comprehensive energy system based on projection integration | |
| CN114048699A (en) | Method for establishing and calculating electric-heating combined electromagnetic transient simulation model of comprehensive energy system | |
| CN115289518A (en) | Heating system thermal and hydraulic balance control method and system | |
| Blizard et al. | A Graph-Based Technique for the Automated Control-Oriented Modeling of District Heating Networks | |
| CN111563332A (en) | Urban water supply pipe network leakage positioning method based on night pressure regulation actual measurement |
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 |