CN117332725A - Steam network dynamic calculation method and device of comprehensive energy system - Google Patents
Steam network dynamic calculation method and device of comprehensive energy system Download PDFInfo
- Publication number
- CN117332725A CN117332725A CN202311629074.7A CN202311629074A CN117332725A CN 117332725 A CN117332725 A CN 117332725A CN 202311629074 A CN202311629074 A CN 202311629074A CN 117332725 A CN117332725 A CN 117332725A
- Authority
- CN
- China
- Prior art keywords
- network
- dynamic
- node
- steam
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 52
- 230000008878 coupling Effects 0.000 claims abstract description 11
- 238000010168 coupling process Methods 0.000 claims abstract description 11
- 238000005859 coupling reaction Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 55
- 230000008569 process Effects 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 13
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 108091027981 Response element Proteins 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- General Business, Economics & Management (AREA)
- Animal Husbandry (AREA)
- Mining & Mineral Resources (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a steam network dynamic calculation method and a steam network dynamic calculation device of a comprehensive energy system. When programming is carried out, the Taylor processing model is packaged into a function, and a time-based curve and a function expression of the waterpower and the heat of each node in the steam network in a scheduling period are obtained by inputting the pressure and the temperature at the heat source point and the curve of the mass flow of each user along with the time change. In a scheduling period, the real-time prediction of the comprehensive energy system on the thermodynamic load in the industrial park is realized by utilizing the time-based function expression of each node, the dynamic coupling with other energy flows is facilitated, and a control model for optimizing the dynamic operation of the whole system is simplified.
Description
Technical Field
The invention relates to the technical field of digital simulation of comprehensive energy systems, in particular to a steam network dynamic calculation method and device of a comprehensive energy system.
Background
In an industrial park, the thermal inertia characteristic of a thermodynamic network is utilized to improve the operation flexibility of the comprehensive energy system (Integrated energy system, IES), and the system heat energy supply is accurately controlled. The thermal inertia characteristics of the thermal network are embodied in the dynamic operation of the system, and the dynamic characteristics of the thermal network need to be completely described. However, the dynamic model of the thermodynamic network presents a high degree of nonlinearity that makes it impossible to solve using analytical methods, which typically require numerical methods. However, the existing numerical solution has high space-time precision and long calculation, and needs to repeatedly iterate continuously to achieve a convergence result, so that the complexity of overall calculation is increased. The existing thermodynamic network model is not suitable for collaborative optimization operation coupled to the integrated energy system. On the one hand, the dynamic nature of the heat supply network makes it difficult to select a reasonable time resolution for both the electrical and thermal subsystems when operated in combination. On the other hand, the nonlinear characteristic of the heat supply network increases the solving difficulty of the IES collaborative operation control model, and the generated non-convex feasible domain makes the original problem difficult to obtain the optimal result of the system operation. Therefore, how to accurately describe the dynamic operation process of the thermal network, a thermal network dynamic calculation method with reasonable complexity is provided, which is a technical problem to be solved in the collaborative optimization operation of the comprehensive energy system in the industrial park.
In addition, the dynamic analysis of thermodynamic networks in IES is currently under study mostly using hot water as the fluid medium, and in fact the industrial park requires steam more than hot water. The dynamic modeling of the existing steam network mostly depends on more mature software (such as FLUENT, AFT, PIPENET and the like) and cannot be applied to the joint modeling of the IES. Although steam networks have a similar topology to large-scale hot water networks, steam networks are more modeling challenging than hot water networks in terms of mode of operation, heating requirements. Therefore, an efficient method is needed for the optimization system of IES to handle the steam network to simplify the dynamic response process of the steam for integration with other components.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a steam network dynamic calculation method and a steam network dynamic calculation device for an integrated energy system so as to improve the system efficiency.
In order to achieve the above purpose, the technical scheme of the invention is as follows:
in a first aspect, the present invention provides a steam network dynamic computing method of an integrated energy system, the method comprising:
establishing a steady-state simplified model of a steam network, and calculating the hydraulic and thermal coupling steady-state process of each node in the network one by one;
and establishing a dynamic processing model of the steam network, and performing dynamic processing based on time factors on the hydraulic and thermal processes of the network.
Further, the establishing a steady-state simplified model of the steam network, and calculating the hydraulic and thermal coupling steady-state process of each node in the network one by one, includes:
s1, acquiring a topological structure of a dendritic steam network in an industrial park, wherein the total node of the network is K, and the total branch number is K-1;
s2, writing the topological structure of the steam network into a form of a directed adjacency list, and storing the topological structure by utilizing a dictionary form in computer programming;
s3, initializing the steam mass flow at the heat source point in the networkTemperature T, pressure P, density->Viscosity->Constant pressure specific heat capacity->And mass flow at the user node +.>；
And S4, calculating hydraulic and thermal parameters of each node and the connected nodes one by utilizing the adjacency list in the step S2.
Further, the establishing a dynamic processing model of the steam network, performing dynamic processing based on time factors on hydraulic and thermal processes of the network, includes:
s5, inquiring i+1 branched by any node i in the heat supply network by utilizing the calculation result of the step S4,judging whether the density and the constant pressure specific heat at the position are converged or not, if yes, turning to the step S6, if not, updating the density and the constant pressure specific heat at the positions of the nodes i+1 and i+2, and returning to the step S4;
s6, setting: the pressure, temperature, mass flow parameters at the beginning of the steam network branch j are expressed as:the terminal pressure, temperature, mass flow parameters are expressed as +.>The method comprises the steps of carrying out a first treatment on the surface of the The steady state moment of the steam network is the initial dynamic moment t 0 The method comprises the steps of carrying out a first treatment on the surface of the The time factor of dynamic calculation is t, and the time period of one-time dynamic calculation is t T The dynamic time node is m=t t The recursion step length in dynamic calculation is n, and n=0 is initialized;
s7, obtaining t from the step S5 0 Starting end parameters of each branch of time steam networkAnd end parametersUsing taylor modelsf h Calculation of t+t 0 The end parameters of the branches at the moment +.>And start parameter；
S8, judging whether the result in the step S7 meets a balance equation, if yes, turning to the step S9, and if not, turning to the step S7, wherein n=n+1;
s9, judging whether the dictionary calculates the last key value pair combination, if so, turning to the step S10, and if not, turning to the step S4;
s10, judging t+t 0 Whether or not it is smaller than t 0 +t T If yes, t= (s+1) t, s=1, 2,..m, and go to step S7, if no, finish the calculation, and output the result.
Further, in the step S2, the dictionary form in the computer programming is used for storing as follows:
wherein, i is { i+1, i+2} is a key value pair combination, which represents the i+1 and i+2 nodes connected with any node i.
Further, in the step S4, hydraulic and thermal parameters at the nodes i+1 and i+2 are calculated in a coupling manner, and the calculation method is as follows:
wherein, subscripts (i, i+1) and (i, i+2) respectively represent branches between the nodes i and i+1 and i+2;l, D is the tube length and diameter of the branch, respectively, for the coefficient of friction resistance; />Constant pressure specific heat capacity and density of steam respectively; η is a local heat dissipation correction parameter; initializing calculation +.>And->Respectively +.>。
Further, in the step S5, the density and the constant pressure specific heat at the nodes i+1, i+2 are queriedJudging that the density and specific heat at the node i+1 and the node i+2 are respectively equal toIf the error between the two reaches the convergence standard, if yes, go to step S6, if no, make，/>，/>，/>The process returns to step S4.
Further, the step S7 is:
in the method, in the process of the invention,f h the calculation model of (2) is as follows:
carrying out space dispersion on any branch with the length L, dividing the branch into N pipe sections, wherein the length of each section is as follows:the method comprises the steps of carrying out a first treatment on the surface of the The hydraulic and thermodynamic parameters for the start and end of any pipe segment k (k=1, 2, … N) are developed as taylor functions based on time factors:
wherein,for the pipeline taylor expansion coefficient, n is the function order and is also the recursion step size of the dynamic calculation.
In a second aspect, the present invention provides a dynamic computing device for a steam network in an industrial park integrated energy system, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method described above when the computer program is executed.
In a third aspect, the present invention provides a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as described above.
Compared with the prior art, the invention has the beneficial effects that:
the hydraulic dynamic process and the thermodynamic dynamic process in the steam network are a set of complex mathematical models based on partial differential equation sets, the partial differential equation is subjected to Taylor expansion processing by using a Taylor function, the equation is simplified into a series of functions based on time factors, and the problem of repeated iteration in numerical solution is effectively avoided by carrying out recursive solution in a limited recursive step length. Meanwhile, when programming is carried out, the Taylor processing model is packaged into a function, and the time-based curve and the function expression of the hydraulic power and the thermal power of each node in the steam network in a scheduling period can be obtained by inputting the pressure and the temperature at the heat source point and the curve of the change of the mass flow of each user along with time. In a scheduling period, the real-time prediction of the comprehensive energy system on the thermodynamic load in the industrial park can be realized by utilizing the time-based function expression of each node, and the dynamic coupling with other energy flows is facilitated, so that a control model for optimizing the dynamic operation of the whole system is simplified, the linearization degree of the control model is improved, and the difficulty of overall solving is reduced. The invention is suitable for the current situation that a plurality of response elements based on different time exist in the comprehensive energy system, fully ensures that the system always obtains the optimal operation effect, and achieves the purpose of improving the system efficiency.
Drawings
FIG. 1 is a flow chart of a steam network dynamic calculation method of the integrated energy system provided in embodiment 1 of the present invention;
FIG. 2 is a topology diagram of a steam network in the integrated energy system of the industrial park provided in embodiment 1 of the present invention;
fig. 3 is a schematic diagram of a dynamic computing device for a steam network in an integrated energy system for an industrial park according to embodiment 2 of the present invention.
Detailed Description
The technical scheme of the invention is further described below with reference to the accompanying drawings and examples.
Example 1:
the steam network dynamic calculation method of the comprehensive energy system mainly comprises two parts, namely, establishing a steady-state simplified model of a steam network, and calculating the hydraulic and thermal coupling steady-state processes of all nodes in the network one by one; and establishing a dynamic processing model of the steam network, and performing dynamic processing based on time factors on the hydraulic and thermal processes of the network. Specifically, referring to fig. 1, the method specifically includes the following steps:
s1: the topological structure of the dendritic steam network in the industrial park is shown in figure 2, wherein the total node of the network is K, and the total branch number is K-1;
s2, writing the topological structure of the steam network into a form of a directed adjacency list, and storing the topological structure in a dictionary form in computer programming:
wherein, i is { i+1, i+2} is a key value pair combination, which represents the i+1 and i+2 nodes connected with any node i.
S3, initializing the steam mass flow at the heat source point in the networkTemperature TPressure P, density->Viscosity->Constant pressure specific heat capacity->And mass flow at the user node +.>；
S4: and (3) calculating the hydraulic power and thermal parameters of each node and the connected nodes one by utilizing the adjacency list in the step (S2). In a key value pair of i { i+1, i+2}, the hydraulic and thermal parameters at the node i are known, and the hydraulic and thermal parameters at the nodes i+1 and i+2 need to be subjected to coupling calculation, wherein the calculation method is as follows:
in the formula, subscripts (i, i+1) and (i, i+2) respectively represent branches between the nodes i and i+1 and i+2.The coefficient of friction, L, D, is the tube length and diameter of the leg, respectively. />The specific heat capacity and the density of the constant pressure of the steam are respectively. η is a local heat dissipation correction parameter. Initializing calculation +.>And->Respectively +.>。
S5: inquiring the node by using the calculation result of the step S4Density and specific heat at constant pressure at i+1, i+2, i.e. Judging that the density and specific heat at the node i+1 and the node i+2 are respectively equal toIf the error between the two reaches the convergence standard, if yes, go to step S6, if no, make，/>，/>，/>The process returns to step S4.
S6: setting: the pressure, temperature, mass flow parameters at the beginning of the steam network branch j are expressed as:the terminal pressure, temperature, mass flow parameters are expressed as +.>The method comprises the steps of carrying out a first treatment on the surface of the The steady state moment of the steam network is the initial dynamic moment t 0 The method comprises the steps of carrying out a first treatment on the surface of the The time factor of dynamic calculation is t, and the time period of one-time dynamic calculation is t T The dynamic time node is m=t t And (3) the recursion step length in dynamic calculation is n, and n=0 is initialized.
S7: obtaining t 0 Starting end parameters of each branch of time steam network,And terminal parameters->Using taylor modelsf h Calculation of t+t 0 Each branch end at momentParameter->And start parameter->The method comprises the following steps:
in the method, in the process of the invention,f h the calculation model of (2) is as follows:
carrying out space dispersion on any branch with the length L, dividing the branch into N pipe sections, wherein the length of each section is as follows:. The hydraulic and thermodynamic parameters for the start and end of any pipe segment k (k=1, 2, … N) are developed as taylor functions based on time factors:
wherein,for the pipeline taylor expansion coefficient, n is the function order and is also the recursion step size of the dynamic calculation.
Therefore, only the Taylor expansion coefficient of the branch start-end parameter needs to be calculatedAnd obtaining a complete time-based function of each parameter. The Taylor expansion coefficient method is as follows:
the outlet pressure Taylor expansion coefficients of the N pipe sections in the nth recursive calculation step length are integrated into a one-dimensional array, and then:
in the method, in the process of the invention,representing the Taylor expansion coefficient of the outlet pressure of the kth pipe section in the nth recursive step, wherein +.>I.e. the end pressure taylor expansion coefficient of branch j,/>。
In the same way, the processing method comprises the steps of,
in the method, in the process of the invention,i.e., the end temperature Taylor expansion coefficient of branch j, ">。I.e. the start temperature taylor expansion coefficient of branch j, ">。
P{n}、Q m The calculation modes of { n } and T { n } are as follows:
in the method, in the process of the invention,
wherein,，/>，/>a is the cross-sectional area of the branch, < >>Is the inertia factor, D is the diameter of the branch, < >>. +.>，/>，/>，The calculation mode of (a) is as follows:
s8: substituting the result in step S7 into the equilibrium equationJudging->Whether or not to converge, i.e.: />If yes, the process goes to step S9, and if not, n=n+1 is entered and goes to step S7.
Wherein,
s9: judging whether the dictionary calculates the last key value pair combination, if so, turning to the step S10, and if not, turning to the step S4.
S10: judging t+t 0 Whether or not it is smaller than t 0 +t T If yes, t= (s+1) t, s=1, 2,..m, and go to step S7, if no, finish the calculation, and output the result.
The hydraulic dynamic process and the thermodynamic dynamic process in the steam network are a set of complex mathematical models based on partial differential equation sets, the partial differential equation is subjected to Taylor expansion processing by using a Taylor function, the equation is simplified into a series of functions based on time factors, and the problem of repeated iteration in numerical solution is effectively avoided by carrying out recursive solution in a limited recursive step length. Meanwhile, when programming is carried out, the Taylor processing model is packaged into a function, and the time-based curve and the function expression of the hydraulic power and the thermal power of each node in the steam network in a scheduling period can be obtained by inputting the pressure and the temperature at the heat source point and the curve of the change of the mass flow of each user along with time. In a scheduling period, the real-time prediction of the comprehensive energy system on the thermodynamic load in the industrial park can be realized by utilizing the time-based function expression of each node, and the dynamic coupling with other energy flows is facilitated, so that a control model for optimizing the dynamic operation of the whole system is simplified, the linearization degree of the control model is improved, and the difficulty of overall solving is reduced. The invention is suitable for the current situation that a plurality of response elements based on different time exist in the comprehensive energy system, fully ensures that the system always obtains the optimal operation effect, and achieves the purpose of improving the system efficiency.
Example 2:
referring to fig. 3, the dynamic computing device for a steam network in an industrial park comprehensive energy system according to the present embodiment includes a processor 21, a memory 22, and a computer program 23 stored in the memory 22 and executable on the processor 21, for example, a dynamic computing program for a steam network in an industrial park comprehensive energy system. The processor 21, when executing the computer program 23, implements the steps of embodiment 1 described above, such as the steps shown in fig. 1.
Illustratively, the computer program 23 may be partitioned into one or more modules/units that are stored in the memory 22 and executed by the processor 21 to complete the present invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions for describing the execution of the computer program 23 in the dynamic computing device for a steam network in an industrial park integrated energy system.
The dynamic computing device for the steam network in the industrial park comprehensive energy system can be computing equipment such as a desktop computer, a notebook computer, a palm computer and a cloud server. The dynamic computing device for the steam network in the industrial park integrated energy system may include, but is not limited to, a processor 21, a memory 22. It will be appreciated by those skilled in the art that fig. 3 is merely an example of a dynamic computing device for a steam network in an industrial park integrated energy system, and is not limiting of a dynamic computing device for a steam network in an industrial park integrated energy system, and may include more or fewer components than shown, or may combine certain components, or different components, e.g., the dynamic computing device for a steam network in an industrial park integrated energy system may also include input and output devices, network access devices, buses, etc.
The processor 21 may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (FieldProgrammable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 22 may be an internal storage element of the dynamic computing device for the steam network in the industrial park integrated energy system, such as a hard disk or a memory of the dynamic computing device for the steam network in the industrial park integrated energy system. The memory 22 may also be an external storage device of the dynamic computing device for the steam network in the industrial park integrated energy system, such as a plug-in hard disk, smart Media Card (SMC), secure Digital (SD) Card, flash Card (Flash Card) or the like, which is provided on the dynamic computing device for the steam network in the industrial park integrated energy system. Further, the memory 22 may also include both internal and external memory devices of the dynamic computing device for the steam network in the industrial park integrated energy system. The memory 22 is used to store the computer program and other programs and data required by the dynamic computing device for the steam network in the industrial park integrated energy system. The memory 22 may also be used to temporarily store data that has been output or is to be output.
Example 3:
the present embodiment provides a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method described in embodiment 1.
The computer readable medium can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer readable medium may even be paper or another suitable medium upon which the program is printed, such as by optically scanning the paper or other medium, then editing, interpreting, or otherwise processing as necessary, and electronically obtaining the program, which is then stored in a computer memory.
The above embodiments are only for illustrating the technical concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement the same, and are not intended to limit the scope of the present invention. All equivalent changes or modifications made in accordance with the essence of the present invention are intended to be included within the scope of the present invention.
Claims (8)
1. A steam network dynamic computing method of an integrated energy system, the method comprising:
establishing a steady-state simplified model of a steam network, and calculating the hydraulic and thermal coupling steady-state process of each node in the network one by one;
establishing a dynamic processing model of a steam network, and performing dynamic processing based on time factors on the hydraulic and thermal processes of the network;
the establishing a steady-state simplified model of the steam network, and calculating the hydraulic and thermal coupling steady-state process of each node in the network one by one comprises the following steps:
s1, acquiring a topological structure of a dendritic steam network in an industrial park, wherein the total node of the network is K, and the total branch number is K-1;
s2, writing the topological structure of the steam network into a form of a directed adjacency list, and storing the topological structure by utilizing a dictionary form in computer programming;
s3, initializing the steam mass flow at the heat source point in the networkTemperature T, pressure P, density->Viscosity->Constant pressure specific heat capacity->And mass flow at the user node +.>；
And S4, calculating hydraulic and thermal parameters of each node and the connected nodes one by utilizing the adjacency list in the step S2.
2. The method for dynamic calculation of steam network of integrated energy system according to claim 1, wherein said creating a dynamic processing model of steam network, performing dynamic processing based on time factors on hydraulic and thermal processes of network, comprises:
s5, inquiring i+1 branched by any node i in the heat supply network by utilizing the calculation result of the step S4,judging whether the density and the constant pressure specific heat of the node i+1, i+2 are converged or not, if yes, turning to the step S6, if not, updating the density and the constant pressure specific heat of the node i+1, i+2, and returning to the step S4;
s6, setting: the pressure, temperature, mass flow parameters at the beginning of the steam network branch j are expressed as:the terminal pressure, temperature, mass flow parameters are expressed as +.>The method comprises the steps of carrying out a first treatment on the surface of the The steady state moment of the steam network is the initial dynamic moment t 0 The method comprises the steps of carrying out a first treatment on the surface of the The time factor of dynamic calculation is t, and one-time dynamic calculation is carried outTime period t T The dynamic time node is m=t t The recursion step length in dynamic calculation is n, and n=0 is initialized;
s7, obtaining t from the step S5 0 Starting end parameters of each branch of time steam networkAnd terminal parameters>Using taylor modelsf h Calculation of t+t 0 The end parameters of the branches at the moment +.>And start parameter->；
S8, judging whether the result in the step S7 meets the balance equation, if yes, turning to the step S9, and if not, makingAnd goes to step S7;
s9, judging whether the dictionary calculates the last key value pair combination, if so, turning to the step S10, and if not, turning to the step S4;
s10, judging t+t 0 Whether or not it is smaller than t 0 +t T If yes, t= (s+1) t, s=1, 2,..m, and go to step S7, if no, finish the calculation, and output the result.
3. The steam network dynamic computing method of an integrated energy system according to claim 1, wherein in the step S2, the dictionary form in the computer programming is used for storing:wherein i { i+1, i+2} is a combination of key value pairs, and represents the combination with any node iAnd the i+1 and i+2 nodes are connected.
4. The method for dynamic calculation of steam network of integrated energy system according to claim 1, wherein in step S4, the hydraulic and thermal parameters at nodes i+1 and i+2 are calculated by coupling, and the calculation method is as follows:
；
wherein, subscripts (i, i+1) and (i, i+2) respectively represent branches between the nodes i and i+1 and i+2;l, D is the tube length and diameter of the branch, respectively, for the coefficient of friction resistance; />Constant pressure specific heat capacity and density of steam respectively; η is a local heat dissipation correction parameter; initializing calculation +.>And->Respectively +.>。
5. The method according to claim 4, wherein in the step S5, the density and the specific heat of constant pressure at the nodes i+1, i+2 are inquiredJudging that the density and specific heat at the node i+1 and the node i+2 are respectively equal to +.>If the error between the two reaches the convergence standard, if yes, go to step S6, if no, let +.>，/>，/>，The process returns to step S4.
6. The steam network dynamic computing method of the integrated energy system according to claim 1 or 5, wherein the step S7 is:
in the method, in the process of the invention, f h the calculation model of (2) is as follows:
carrying out space dispersion on any branch with the length L, dividing the branch into N pipe sections, wherein the length of each section is as follows:the method comprises the steps of carrying out a first treatment on the surface of the The hydraulic and thermodynamic parameters for the start and end of any pipe segment k (k=1, 2, … N) are developed as taylor functions based on time factors:
wherein,the Taylor expansion coefficient of the pipeline, n is the function order and is also the recursion step size of dynamic calculation.
7. A dynamic computing device for a steam network in an industrial park integrated energy system, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method of claim 1 when the computer program is executed.
8. A computer readable storage medium storing a computer program, which when executed by a processor performs the steps of the method according to claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311629074.7A CN117332725B (en) | 2023-12-01 | 2023-12-01 | Steam network dynamic calculation method and device of comprehensive energy system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311629074.7A CN117332725B (en) | 2023-12-01 | 2023-12-01 | Steam network dynamic calculation method and device of comprehensive energy system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117332725A true CN117332725A (en) | 2024-01-02 |
| CN117332725B CN117332725B (en) | 2024-07-02 |
Family
ID=89279684
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311629074.7A Active CN117332725B (en) | 2023-12-01 | 2023-12-01 | Steam network dynamic calculation method and device of comprehensive energy system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117332725B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102444784A (en) * | 2011-11-14 | 2012-05-09 | 上海金自天正信息技术有限公司 | Pressure control system for steel enterprise steam pipe network based on dynamic matrix control |
| CN108599137A (en) * | 2017-12-28 | 2018-09-28 | 东南大学 | A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character |
| CN110555264A (en) * | 2019-08-30 | 2019-12-10 | 南方电网科学研究院有限责任公司 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
| CN113283077A (en) * | 2021-05-21 | 2021-08-20 | 福州大学 | Correlation-considered power flow calculation method for electricity-heat-hydrogen comprehensive energy system interval |
| CN113515853A (en) * | 2021-06-08 | 2021-10-19 | 国网江苏省电力有限公司连云港供电分公司 | Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation |
-
2023
- 2023-12-01 CN CN202311629074.7A patent/CN117332725B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102444784A (en) * | 2011-11-14 | 2012-05-09 | 上海金自天正信息技术有限公司 | Pressure control system for steel enterprise steam pipe network based on dynamic matrix control |
| CN108599137A (en) * | 2017-12-28 | 2018-09-28 | 东南大学 | A kind of multipotency streaming system optimizing operation method considering region heat supply network transient state heat-transfer character |
| CN110555264A (en) * | 2019-08-30 | 2019-12-10 | 南方电网科学研究院有限责任公司 | Dynamic simulation method and device of comprehensive energy heat supply system based on extended node method |
| CN113283077A (en) * | 2021-05-21 | 2021-08-20 | 福州大学 | Correlation-considered power flow calculation method for electricity-heat-hydrogen comprehensive energy system interval |
| CN113515853A (en) * | 2021-06-08 | 2021-10-19 | 国网江苏省电力有限公司连云港供电分公司 | Electric-heat interconnection comprehensive energy system optimization scheduling method based on linear equation |
Non-Patent Citations (3)
| Title |
|---|
| WEIJIA YANG, YUPING HUANG, DAIQING ZHAO: "A coupled hydraulic –thermal dynamic model for the steam network in a heat–electricity integrated energy system", ENERGY，HTTPS://DOI.ORG/10.1016/J.ENERGY.2022. 125800, pages 1 - 18 * |
| 张增刚;李继志;李永安;: "城市蒸汽管网水力热力耦合系统及其应用", 微计算机信息, no. 18, pages 207 - 209 * |
| 陈彬彬;孙宏斌;吴文传;郭庆来;乔铮;: "综合能源系统分析的统一能路理论(三)：稳态与动态潮流计算", 中国电机工程学报, vol. 40, no. 15, pages 4820 - 4831 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117332725B (en) | 2024-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109242365B (en) | Interval power flow calculation method of electricity-heat interconnection comprehensive energy system | |
| Niu et al. | A new control method for state-constrained nonlinear switched systems with application to chemical process | |
| Michel et al. | Stability of dynamical systems | |
| Wei et al. | New approach to delay‐dependent H∞ filtering for discrete‐time Markovian jump systems with time‐varying delay and incomplete transition descriptions | |
| CN107609141B (en) | Method for performing rapid probabilistic modeling on large-scale renewable energy data | |
| Kitzhofer et al. | The new Matlab code bvpsuite for the solution of singular implicit BVPs | |
| CN107871058A (en) | Tidal current computing method, device, equipment and the storage medium of electric heating association system | |
| Poirier et al. | Efficient reduced-radial basis function-based mesh deformation within an adjoint-based aerodynamic optimization framework | |
| Chen et al. | MGNet: a novel differential mesh generation method based on unsupervised neural networks | |
| Wei et al. | Modelling and simulation of nabla fractional dynamic systems with nonzero initial conditions | |
| CN113111515A (en) | Unified modeling method of comprehensive energy system | |
| Mayes et al. | Approximation of potential-driven flow dynamics in large-scale self-similar tree networks | |
| Feng | Parameter independent model order reduction | |
| Li et al. | Generalized linear‐constrained optimal power flow for distribution networks | |
| CN109145441B (en) | Thermoelectric energy supply calculation method and device, electronic equipment and storage medium | |
| CN117332725B (en) | Steam network dynamic calculation method and device of comprehensive energy system | |
| Zahra et al. | Spline difference scheme for two-parameter singularly perturbed partial differential equations | |
| Yuan et al. | Stabilization of coupled ODE‐PDE system with intermediate point and spatially varying effects interconnection | |
| Kurth et al. | Inversion-based and optimal feedforward control for population dynamics with input constraints and self-competition in chemostat reactor applications | |
| CN113221428A (en) | Rapid decomposition method for dynamic energy flow calculation of electricity-heat comprehensive energy system | |
| Wang et al. | Optimal reduction of models using a hybrid searching strategy | |
| CN109611949B (en) | Pipeline thermal calculation method and monitoring and heating system adopting same | |
| Aziz et al. | Polynomial interpolation in matlab | |
| Monaco et al. | Input-state matching under digital control | |
| Zhang et al. | Truncation error calculation based on Richardson extrapolation for variable-step collaborative simulation |
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 |