CN112035991B - Steam optimization calculation method and system based on pipe network conveying path - Google Patents
Steam optimization calculation method and system based on pipe network conveying path Download PDFInfo
- Publication number
- CN112035991B CN112035991B CN202011012525.9A CN202011012525A CN112035991B CN 112035991 B CN112035991 B CN 112035991B CN 202011012525 A CN202011012525 A CN 202011012525A CN 112035991 B CN112035991 B CN 112035991B
- Authority
- CN
- China
- Prior art keywords
- steam
- path
- pipe network
- flow
- matrix
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 46
- 238000005457 optimization Methods 0.000 title claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 58
- 238000000034 method Methods 0.000 claims abstract description 50
- 239000012530 fluid Substances 0.000 claims abstract description 48
- 230000008901 benefit Effects 0.000 claims abstract description 21
- 230000005494 condensation Effects 0.000 claims abstract description 21
- 238000009833 condensation Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 230000000704 physical effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 abstract description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 6
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000007726 management method Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening 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/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/02—CAD in a network environment, e.g. collaborative CAD or distributed simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention relates to a steam optimization calculation method and system based on a pipe network conveying path, and belongs to the field of automation. The method comprises the following steps: firstly, physical information of a steam pipe network is displayed in a topological graph mode, and then the pipe network information is displayed in a matrix mode in a connection mode between a pipeline and a node. In the network connection matrix, all the conveying paths from the gas source point to the user point are quickly and accurately found by using a method based on actual fluid flow. And finally, calculating the flow of each path by combining the resistance loss coefficient in each path, and obtaining the condensation loss of each pipeline through heat transfer calculation. The benefits of steam utilization are maximized by adjusting the amount of steam at different steam sources and users in combination with the benefits of steam use and the losses of steam condensation. The method can calculate each conveying path and loss thereof between the energy medium gas source point and the user point more quickly and accurately, and provide data support for steam optimal scheduling.
Description
Technical Field
The invention belongs to the field of automation, and relates to a steam optimization calculation method and system based on a pipe network conveying path.
Background
The steam is utilized without separating the conveying process of the pipe network, and the energy medium pipe network in the iron and steel enterprises is usually a multi-steam source, multi-user and long-span complex pipe network. Because of the lack of metering inside the pipe network, it is often difficult to ascertain the delivery process and delivery path of the fluid in the pipe network. In particular, the steam pipe network often causes a large amount of condensation loss in the steam conveying process, and the condensation loss caused by different conveying paths is also different, which makes optimal scheduling of the steam very difficult. Therefore, the calculation of the fluid conveying paths between each air source and the user and the steam optimization method and system based on the calculation are very important for energy management and steam optimization of iron and steel enterprises.
In the prior art, a method for finding paths by traversing all nodes of a pipe network by taking a certain component of the pipe network as a starting point and extracting all non-coincident paths is provided. When the existing traversing method traverses the component, the current path is stopped traversing, and a plurality of repeated paths which are not consistent with the actual flow exist in the paths, so that the paths are required to be further screened and removed; after the preliminary screening and culling are completed, node attributes in each path are also required to be: and (3) further judging the starting point (Start) or the ending point (End) to finally obtain an effective path formed by the effective nodes. Such traversal methods tend to find a large number of invalid paths, which not only wastes a large amount of computation resources, but also lengthens the overall program running time.
The following describes the disadvantages of this traversal method in conjunction with the specific illustrations:
as shown in fig. 1, in the steam pipe network, the dots represent steam source points, the squares represent steam user points, and the following cases may occur when path searching is performed in a traversing manner:
from the source point (Start) (1) back to the source point (Start) (2) in the path found in fig. 2, such a path is virtually nonexistent because steam is unlikely to enter the source point. The traversing method can judge the attribute of each node in the path: the Start point (Start) and the user point (End) cull such invalid paths. However, the process of traversing searching and judging rejection consumes more operation time.
Alternatively, it is possible to find the node (1) to node as shown in fig. 3 by a traversal methodAlso, the path satisfies that only one source point (Start) and one user point (End) exist. If the path is found by traversal, this is an effective path, but in practice the arrows in the path may be in the opposite direction to the fluid flow in the pipe, resulting in the path being virtually nonexistent or ineffective.
In some existing steam pipe network calculation methods, the condensation water quantity in the steam pipe network can be obtained through hydraulic-thermal coupling calculation, for example, after a correlation matrix is built, an admittance matrix between flow and impedance is written by utilizing a flow node equation set and the relation between pipe section flow and pressure drop, and hydraulic calculation is carried out by a method of multiplying the matrix and then solving an equation; writing out the relation between the flow and the temperature drop of the pipe section by the same method, writing into an admittance matrix between the flow and the thermal resistance, and performing thermal calculation by multiplying the matrix and then solving an equation; the process is iterated continuously, so that the hydraulic-thermal coupling calculation is formed to obtain the steam condensate water, the calculation amount of the hydraulic-thermal iterative coupling calculation is very large, and the calculation time required to be consumed is relatively long.
In the process of optimizing and dispatching steam, the flow rates of a steam source and a user point are required to be continuously adjusted, and if hydraulic-thermal iterative coupling calculation is required after each adjustment, the whole optimization process becomes very difficult due to the consumption of a large amount of calculation resources and calculation time, and the real-time change of a steam pipe network cannot be responded.
Disclosure of Invention
In view of the above, the present invention aims to provide a steam optimization calculation method and system based on a pipe network conveying path. The physical information of the energy medium pipe network is displayed in a topological graph mode, and the pipe network information is displayed in a matrix mode through a connection mode between the pipeline and the nodes. In the network connection matrix, all the conveying paths from the gas source point to the user point are quickly and accurately found by using a method based on actual fluid flow. And finally, calculating the flow of each path by combining the resistance loss coefficient in each path, and obtaining the condensation loss of each pipeline through heat transfer calculation. The benefits brought by the steam use and the losses brought by the steam condensation are combined, and the benefits of the steam utilization are maximized by adjusting the steam quantity of different steam sources and user points. The method can calculate each conveying path and loss thereof between the energy medium gas source point and the user point more quickly and accurately, provides reference for steam dispatching, and greatly improves the efficiency and the fineness of energy management of iron and steel enterprises.
In order to achieve the above purpose, the present invention provides the following technical solutions:
a steam optimization calculation method based on a pipe network conveying path comprises the following steps:
s1: constructing a pipe network topological graph through the geographic position of a steam pipe network, and representing pipe network information in the pipe network topological graph in a matrix mode according to the connection mode of the pipeline and the node, namely a pipe network information matrix;
s2: obtaining the flow condition of fluid in the pipeline through hydraulic calculation, and judging the flow direction of fluid medium of each pipeline in the pipeline network;
s3: combining the pipe network information matrix with the flowing direction to form a new pipe network association matrix;
s4: finding all paths from the gas source point to the user point by utilizing the actual fluid flow direction through the pipe network association matrix;
s5: according to the physical parameters of the pipeline on each path: the length of the pipeline, the diameter of the pipeline and the roughness of the pipeline are used for calculating the resistance coefficient of each path and calculating the pressure drop and the flow of the fluid on the path.
Optionally, in the step S1, when the pipe network information matrix is constructed:
firstly numbering each node and pipe section point;
the numbers of the nodes and the pipe section points form a two-dimensional matrix, the number is represented by a number of 1 if the pipe section and the nodes have a connection relationship, and the number is represented by a number of 0 if the pipe section and the nodes have no connection relationship.
Optionally, in the step S3, when the pipe network information matrix is combined with the actual fluid flow direction to form a new pipe network association matrix:
carrying out hydraulic calculation according to data measured by the steam source point and the user point instrument to obtain the flow direction of the fluid in each pipe section; the numbers in the pipe section information matrix are marked as "1" and "-1" according to the relation of the flow in and out.
Optionally, in the step S4, when the correlation matrix is used to find the flow path:
each pipe section can only correspond to one inlet and one outlet, but one node can be communicated with multiple pipe sections; by utilizing the property and combining the principle of 'one-to-one and different number connection', the actual fluid flow direction and the expression form in the correlation matrix are fully utilized, and all paths from a steam source to a user point can be quickly and accurately found in the two-dimensional correlation matrix;
optionally, in S5, when calculating the pressure drop and the flow rate using the resistance coefficient of each path:
the pressure drop equation on each path can be listed and the flow rate of each path can be found by multiplying the flow rate by the drag coefficient and the physical formula of the pressure drop and the pressure difference between the air source and the user is constant:
ΔP＝∑ζ i1 Q 1 ＝∑ζ i2 Q 2 ＝…＝∑ζ i3 Q n
Q＝Q 1 +Q 2 +…+Q n
and (3) calculating the condensation loss in the steam conveying process on each path, and optimizing the steam conveying scheme by combining the benefits of the steam, so that the overall use benefits of the steam are maximized.
Objective function: max { Benefit Steam generation -Cost Condensation }。
A steam optimization computing system based on a pipe network transport path, comprising:
the acquisition module is used for acquiring real-time state data of fluid states of all air sources and user nodes in the fluid pipe network and physical parameter data used for representing the pipe section structure, the physical properties of the fluid pipe network and the connection mode between the pipe sections;
a data processing module for processing the collected data,
the output module is used for outputting a data processing result;
the acquisition module, the data processing module and the output module are sequentially connected.
Optionally, the system further comprises:
the display module is used for visualizing all paths between all air sources and the user point and displaying all flow information on each path in real time;
a data list module;
wherein the data list module comprises:
the real-time data list is used for storing real-time state data;
a physical parameter list for storing physical parameter data;
the result data list is used for storing all the calculated results, including each path from the gas source point to the user point and flow information on the path;
and displaying all paths between the air source and the user and the flow corresponding to each path through the display module according to the flow and pressure information measured by the metering instruments at the air source and the user.
A computer readable storage medium having stored thereon a computer program which when executed by a processor implements the method.
An electronic terminal, comprising: a processor and a memory;
the memory is used for storing a computer program, and the processor is used for executing the computer program stored in the memory so as to enable the terminal to execute the method.
The invention has the beneficial effects that: the calculation of steam optimization is carried out based on the pipe network conveying path, and after the flow of the steam source and the flow pointed out by a user are regulated each time, the change of the flow and the condensation amount on the corresponding path is calculated, and the hydraulic-thermal coupling calculation of the whole pipe network is not needed. Therefore, the calculation resources required by the optimization process can be greatly reduced, and the calculation program can respond to the real-time change of the steam pipe network better.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objects and other advantages of the invention may be realized and obtained by means of the instrumentalities and combinations particularly pointed out in the specification.
Drawings
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in the following preferred detail with reference to the accompanying drawings, in which:
FIG. 1 is a diagram of an example pipe network;
FIG. 2 is a path example one;
FIG. 3 is a second path example;
FIG. 4 is a schematic diagram of a multi-source multi-user complex fluid pipe network topology;
FIG. 5 is a schematic diagram of a network information matrix;
FIG. 6 is a pipe network association matrix;
FIG. 7 is a flow path from a pipe network incidence matrix gas source point "(1)" to a user point "(9)";
FIG. 8 is a schematic diagram of the air source point "(1)" to user point "(9)" path;
FIG. 9 is a schematic diagram of a pipe network topology path;
FIG. 10 is an illustration of a "foreign link" pipe network topology;
fig. 11 is an explanatory diagram of the "foreign-number-connected" pipe network association matrix.
Detailed Description
Other advantages and effects of the present invention will become apparent to those skilled in the art from the following disclosure, which describes the embodiments of the present invention with reference to specific examples. The invention may be practiced or carried out in other embodiments that depart from the specific details, and the details of the present description may be modified or varied from the spirit and scope of the present invention. It should be noted that the illustrations provided in the following embodiments merely illustrate the basic idea of the present invention by way of illustration, and the following embodiments and features in the embodiments may be combined with each other without conflict.
Wherein the drawings are for illustrative purposes only and are shown in schematic, non-physical, and not intended to limit the invention; for the purpose of better illustrating embodiments of the invention, certain elements of the drawings may be omitted, enlarged or reduced and do not represent the size of the actual product; it will be appreciated by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.
The same or similar reference numbers in the drawings of embodiments of the invention correspond to the same or similar components; in the description of the present invention, it should be understood that, if there are terms such as "upper", "lower", "left", "right", "front", "rear", etc., that indicate an azimuth or a positional relationship based on the azimuth or the positional relationship shown in the drawings, it is only for convenience of describing the present invention and simplifying the description, but not for indicating or suggesting that the referred device or element must have a specific azimuth, be constructed and operated in a specific azimuth, so that the terms describing the positional relationship in the drawings are merely for exemplary illustration and should not be construed as limiting the present invention, and that the specific meaning of the above terms may be understood by those of ordinary skill in the art according to the specific circumstances.
The summary will now be described with reference to a fluid network as shown in fig. 4:
the dots as shown in the above figures represent the fluid gas source, the squares represent the fluid user, and the black lines represent the fluid conduits. Each node and each segment of pipe is first numbered according to a complex network of multiple sources and users comprising multiple loops. The numbers in the circles are the numbers of the nodes, such as "(1)"; the individual numbers represent the numbers of the pipe segments, such as '1', the numbers of the nodes and pipe segments may be arranged in any order, but it must be ensured that each number is independent and consecutive.
In a practical large pipe network, only an air source point is usually arrangedOr the user has a metering device, and the pipe network is internally provided with no metering device, as shown in the figure as "(1) (2) (3) (6) (9)"there is a meter at the node, and" (4) (5) (7) (8) is->"where there is no metering device. Because of the lack of metering inside the pipe network, it is difficult to ascertain the flow paths of the fluid in the pipe network and the flow rate of each path, which makes energy fine management of enterprises more difficult.
According to the method, firstly, through a conveying path calculation method based on the actual fluid flowing direction, all paths from a steam source point to a user point and the flow of each path can be rapidly calculated, so that the flowing condition of fluid in a pipe network can be known, the refinement degree of enterprise energy management is improved, and the method is concretely implemented as follows:
according to the numbers and connection modes of the nodes and the pipe sections, the pipe network topological graph is expressed in a matrix form, the expressed matrix is shown in fig. 5, in a two-dimensional matrix formed according to the numbers and the connection modes, the ordinate represents the node numbers, the abscissa represents the pipe section numbers, and because the matrix is generated through calculation of a program, all data in the program starts from 0, the numbers in the matrix need to be added with 1 to correspond to the numbers in the pipe network, for example, the number ' 0 ' on the ordinate of the matrix represents the node with the number (1 '). The number "1" inside the matrix indicates that the node is connected to the pipe, and the number "0" indicates that the node is not connected to the pipe. For example, in the figure, the two numbers "1" in column 1 indicate that node "(1)" and node "(4)" are connected by pipe segment "1".
From the hydraulic calculation of the pipe, i.e. the fluid always flows from the point of high pressure to the point of low pressure, the direction of flow of the fluid in all pipe sections can be derived, as indicated by the arrows in fig. 4. Combining the fluid flow direction with the network information matrix, a network correlation matrix can be obtained, as shown in fig. 6, where the correlation matrix is represented by the numbers "-1" and "1" to indicate the fluid flow direction in the pipeline, "1" indicates the inflow, "1" indicates the outflow, and the number "1" at (0, 0) indicates that the fluid flows from node "(1)" into the pipeline "1".
Then, all paths from the gas source point to the user point can be found in the incidence matrix by using a traversing method, and the rule for finding the paths is based on the characteristic that each pipe section of the incidence matrix can only correspond to one inlet and one outlet, but one node can be communicated with multiple pipe sections. By using this characteristic, the whole path from the steam source point to the user point can be quickly and accurately found in the two-dimensional correlation matrix by a path calculation method taking the actual fluid flow into consideration. The following is a detailed description:
if we need to find all the flow paths from the gas source point "(1)" to the user point "(9)" in fig. 4, the representation in the correlation matrix is as shown in fig. 7, and the number "1" of the gas source point "(1) is connected with all other numbers according to the" one-to-one-cross, different number connection "principle, namely, the flow path is found, firstly, starting from the number" 1 "of the gas source point. For example, starting from the numeral "1" at the coordinates (0, 0), all the opposite numerals "-1" of the column are connected longitudinally, and then "-1" at the coordinates (0, 3) can be reached; from the "—1" at (0, 3) cross-connect all the different numbers "1" of the row, three "1" at coordinates (2, 3), (3, 3) and (4, 3) can be found, respectively, representing three paths outward from node "(4)". The method is repeated, longitudinal connection and transverse connection are alternately performed each time, all different numbers in the same row or the same column are connected by each number, and a path is ended when no connection is performed any more. If the '1' at the coordinates (2, 3) can be longitudinally connected with the '1' at the (2, 2), but when the '1' at the (2, 2) needs to be transversely connected again, all numbers of the row are ' 0 ', the connection is not met, the connection is ended, all nodes passing under the path, namely the sequence nodes ' (1), (4) and (3), are recorded, and the path from the air source point ' (1) to the user point ' (3).
The above is an explanation of the principle of "one-to-one-horizontal, different number connection", taking all paths from the source point "(1)" to the user point "(9)" as an example, two possible routes can be found by the above method, as shown in fig. 8:
as shown in fig. 8, two paths can be found, where the two paths are different and are respectively represented by a dashed arrow and an implementation arrow, and the coordinates of the first path in the association matrix passing through in turn are [ (0, 0), (0, 3), (3, 4), (5, 7), (8, 8) ], and then the corresponding path node numbers are [ "(1), (4), (5), (8), (9)" ]; the coordinates of the other path in the association matrix sequentially pass through are [ (0, 0), (0, 3), (4, 6), (7, 7), (8, 8) ], and the corresponding path node numbers are [ "(1), (4), (7), (8), (9)" ]. The representation of the two paths on the network topology is shown in fig. 9.
The two flow paths [ "(1), (4), (5), (8), (9)" ] and [ "(1), (4), (7), (8), (9)" ] found by the correlation matrix can correspond exactly to the paths on the pipe network topology.
Finally, taking all paths from the source point (1) to the user point (9) as an example, the principle of "different number connection" must be satisfied. Then a path from the air source point to the point (1) to the user point (9) is [ "(1), (4), (7), or (2),⑧、⑨”]Its route on the network topology is shown in fig. 10.
In the paths [ "(1), (4), (7), gamma,⑧、⑨”]The flow direction of the pipe section 11 is calculated hydraulically through the pipe network, so that the fluid in the pipe section can only flow from the node 8 to the node>And it is impossible to flow in the opposite direction. The expression of such an in-network topology in the network association matrix is shown in fig. 11.
When the path finds the "-1" of the coordinates (11, 10), the number at the coordinates (10, 10) is "-1", so that the principle of "different number connection" is not satisfied, i.e. the path does not exist, and the physical meaning of the principle of "different number connection" in the matrix corresponds to the pipe section flow direction in the pipe network topological graph, and any path finding process cannot violate the actual flow direction of the fluid in the pipe section, otherwise the corresponding physical meaning of the fluid is lost.
Therefore, there are only two paths from the source point "(1)" to the user point "(9)", and the corresponding node sequences are respectively: [ "(1), (4), (5), (8), (9)".]And [ "(1), (4), (7), (8), (9)".]. Since the path calculation based on the actual fluid flow direction can be performed starting from any desired vapor source point in the path calculation process, all flow paths from all desired vapor source points to the user point can be found. For example, the gas source point (1) "is found by traversing, and all the user points (3), (9) are found in turn,"all paths; if necessary, the source point (2) is updated again, and all user points (3), (9) and +.>"all paths. Similarly, each flow path from the vapor source point to the user point may be reversed for all that is needed. The method is also suitable for a more complex large-scale fluid pipe network.
And finally, calculating the flow of each path by using the relation between the resistance and the flow through the following formula:
ΔP＝∑ζ i1 Q 1 ＝∑ζ i2 Q 2 ＝…＝∑ζ i3 Q n
Q＝Q 1 +Q 2 +…+Q n
wherein DeltaP is the pressure difference between the air source point and the user point, ζ i1 Representing the drag coefficient of the "i" th tube segment on path 1, Σζ i1 The drag coefficients of all the pipe segments that make up path 1 are summed, wherein the drag of each pipe segmentThe force coefficient is a fixed parameter determined by the physical parameters of the pipe section, such as diameter, roughness, length, etc. Q (Q) 1 Representing the fluid flow of path 1; Σζ i1 Q 1 ＝∑ζ i2 Q 2 ＝…＝∑ζ i2 Q n Identity means: since the pressure difference Δp between the source point and the user point is constant (measured by the pressure meter), the flow of each path is equal to the sum of the drag coefficients, and the path flow with a large drag coefficient is small and the path flow with a small drag coefficient is large. (for ease of understanding, the relationship between voltage, current and resistance may be analogized: the voltage between two points is constant, the current on a circuit with small resistance is larger, and the current on a circuit with large resistance is smaller); q represents the total flow of the air source to the user point, i.e. equal to the sum of all flows of each path.
Finally, after all steam paths from the steam source point to the user point and after the flow are obtained, the corresponding condensation water quantity on each path can be obtained through heat transfer calculation, and the formula is as follows:
wherein P is i For the heat exchange amount with the outside on the ith path (the calculation of the heat exchange amount is described in detail in the heat transfer science, and is not described in detail in the patent), r is the latent heat of steam. The calculated condensation is the sum of the condensation water amounts on all paths from a certain steam source to the user point.
After the condensate water amount in the conveying process is calculated, the steam conveying scheme is optimized in combination with the benefits of steam, so that the overall use benefits of the steam are maximized, and the objective function is as follows: max { Benefit Steam generation -Cost Condensation }. Wherein Benefit is Steam generation For benefits of steam to be used for heating or production after reaching customer points, cost Condensation Refers to the loss of steam due to condensation during transportation. By continuously adjusting the steam source and the steam user, and then calculating the condensation loss on each conveying path, a new Benefit can be obtained Steam generation And Cost Condensation Until a large overall benefit value is obtained that meets the demand, the steam optimization process ends. The steam optimization calculation method and the steam optimization calculation system based on the pipe network conveying path have the greatest characteristics of saving a large amount of calculation resources and calculation time and greatly improving the real-time response capability of the system.
The invention is based on the pipe network conveying path to perform steam optimization calculation, and after the flow of the steam source and the flow pointed out by a user are regulated each time, the change of the flow and the condensation quantity on the corresponding path is calculated, and the hydraulic-thermal coupling calculation of the whole pipe network is not needed. Therefore, the calculation resources required by the optimization process can be greatly reduced, and the calculation program can respond to the real-time change of the steam pipe network better.
Hardware device
To accomplish the above-mentioned inventive functions, a set of computing systems is formed by matching corresponding hardware devices, including:
the acquisition module is used for respectively acquiring real-time state data of fluid states of all air sources and user nodes (namely, metered nodes) in the fluid pipe network and physical parameter data used for representing the pipe section structure, the physical properties (length, diameter and the like) of the fluid pipe network and the connection mode between the pipe sections;
a data processing module for processing the collected data,
the output module is used for outputting a data processing result;
optionally, the method further comprises:
and the display module is used for visualizing all paths between all air sources and the user point and displaying all flow information on each path in real time.
A data list module, the data list module comprising:
the real-time data list is used for storing real-time state data;
a physical parameter list for storing physical parameter data;
the result data list is used for storing all the calculated results, such as information of each path from the gas source point to the user point, flow on the path and the like;
and displaying all paths between the air source and the user and the flow corresponding to each path through the display module according to the information such as the flow and the pressure measured by the metering instruments at the air source and the user.
The present invention also provides a computer-readable storage medium having stored thereon a computer program characterized in that: the program, when executed by a processor, implements the method of any of the above.
The invention also provides an electronic terminal, which is characterized by comprising: a processor and a memory;
the memory is configured to store a computer program, and the processor is configured to execute the computer program stored in the memory, so as to cause the terminal to execute the method according to any one of the above methods.
The computer readable storage medium in this embodiment, as will be appreciated by those of ordinary skill in the art: all or part of the steps for implementing the method embodiments described above may be performed by computer program related hardware. The aforementioned computer program may be stored in a computer readable storage medium. The program, when executed, performs steps including the method embodiments described above; and the aforementioned storage medium includes: various media that can store program code, such as ROM, RAM, magnetic or optical disks.
The electronic terminal provided in this embodiment includes a processor, a memory, a transceiver, and a communication interface, where the memory and the communication interface are connected to the processor and the transceiver and complete communication with each other, the memory is used to store a computer program, the communication interface is used to perform communication, and the processor and the transceiver are used to run the computer program, so that the electronic terminal performs each step of the above method.
In this embodiment, the memory may include a random access memory (RandomAccess Memory, abbreviated as RAM), and may further include a non-volatile memory (non-volatile memory), such as at least one magnetic disk memory.
The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a Network Processor (NP), etc.; but also digital signal processors (Digital Signal Processing, DSP for short), application specific integrated circuits (Application Specific Integrated Circuit, ASIC for short), field-programmable gate arrays (Field-Programmable Gate Array, FPGA for short) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the present invention, which is intended to be covered by the claims of the present invention.
Claims (7)
1. A steam optimization calculation method based on a pipe network conveying path is characterized by comprising the following steps of: the method comprises the following steps:
s1: constructing a pipe network topological graph through the geographic position of a steam pipe network, and representing pipe network information in the pipe network topological graph in a matrix mode according to the connection mode of the pipeline and the node, namely a pipe network information matrix;
s2: obtaining the flow condition of fluid in the pipeline through hydraulic calculation, and judging the flow direction of fluid medium of each pipeline in the pipeline network;
s3: combining the pipe network information matrix with the flowing direction to form a new pipe network association matrix;
in the step S3, when the pipe network information matrix is combined with the actual fluid flow direction to form a new pipe network association matrix:
carrying out hydraulic calculation according to data measured by the steam source point and the user point instrument to obtain the flow direction of the fluid in each pipe section; the numbers in the pipe section information matrix are marked as '1' and '-1' according to the relation of the flow in and flow out;
s4: finding all paths from the gas source point to the user point by utilizing the actual fluid flow direction through the pipe network association matrix;
in the step S4, when the flow path is found by using the correlation matrix:
each pipe section can only correspond to one inlet and one outlet, but one node can be communicated with multiple pipe sections; by utilizing the property and combining the principle of 'one-to-one and different number connection', the actual fluid flow direction and the expression form in the correlation matrix are fully utilized, and all paths from a steam source to a user point can be quickly and accurately found in the two-dimensional correlation matrix;
s5: according to the physical parameters of the pipeline on each path: the length of the pipeline, the diameter of the pipeline and the roughness of the pipeline are used for calculating the resistance coefficient of each path and calculating the pressure drop and the flow of the fluid on the path.
2. The steam optimization calculation method based on the pipe network conveying path according to claim 1, wherein the method comprises the following steps of: in the step S1, when the pipe network information matrix is constructed:
firstly numbering each node and pipe section point;
the numbers of the nodes and the pipe section points form a two-dimensional matrix, the number is represented by a number of 1 if the pipe section and the nodes have a connection relationship, and the number is represented by a number of 0 if the pipe section and the nodes have no connection relationship.
3. The steam optimization calculation method based on the pipe network conveying path according to claim 1, wherein the method comprises the following steps of: in the step S5, when the pressure drop and the flow rate are calculated using the resistance coefficient of each path:
the pressure drop equation on each path can be listed and the flow rate of each path can be found by multiplying the flow rate by the drag coefficient and the physical formula of the pressure drop and the pressure difference between the air source and the user is constant:
ΔP＝∑ζ i1 Q 1 ＝∑ζ i2 Q 2 ＝…＝∑ζ i3 Q n
Q＝Q 1 +Q 2 +…+Q n
the condensation loss in the steam conveying process on each path is calculated, and the steam conveying scheme is optimized by combining the benefits of the steam, so that the overall use benefits of the steam are maximized;
P i the heat exchange quantity with the outside on the ith path; r is the latent heat of steam;
objective function: max { Benefit Steam generation -Cost Condensation }。
4. A steam optimization computing system based on a pipe network conveying path based on the method of any one of claims 1 to 3, characterized in that: the system comprises:
the acquisition module is used for acquiring real-time state data of fluid states of all air sources and user nodes in the fluid pipe network and physical parameter data used for representing the pipe section structure, the physical properties of the fluid pipe network and the connection mode between the pipe sections;
a data processing module for processing the collected data,
the output module is used for outputting a data processing result;
the acquisition module, the data processing module and the output module are sequentially connected.
5. The steam optimization computing system based on a pipe network transport path of claim 4, wherein: the system further comprises:
the display module is used for visualizing all paths between all air sources and the user point and displaying all flow information on each path in real time;
a data list module;
wherein the data list module comprises:
the real-time data list is used for storing real-time state data;
a physical parameter list for storing physical parameter data;
the result data list is used for storing all the calculated results, including each path from the gas source point to the user point and flow information on the path;
and displaying all paths between the air source and the user and the flow corresponding to each path through the display module according to the flow and pressure information measured by the metering instruments at the air source and the user.
6. A computer-readable storage medium having stored thereon a computer program, characterized by: the program, when executed by a processor, implements the method of any one of claims 1 to 3.
7. An electronic terminal, characterized in that: comprising the following steps: a processor and a memory;
the memory is used for storing a computer program, and the processor is used for executing the computer program stored in the memory, so that the terminal executes the method of any one of 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011012525.9A CN112035991B (en) | 2020-09-23 | 2020-09-23 | Steam optimization calculation method and system based on pipe network conveying path |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011012525.9A CN112035991B (en) | 2020-09-23 | 2020-09-23 | Steam optimization calculation method and system based on pipe network conveying path |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112035991A CN112035991A (en) | 2020-12-04 |
| CN112035991B true CN112035991B (en) | 2024-02-27 |
Family
ID=73575146
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011012525.9A Active CN112035991B (en) | 2020-09-23 | 2020-09-23 | Steam optimization calculation method and system based on pipe network conveying path |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112035991B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8612189B2 (en) * | 2005-10-03 | 2013-12-17 | Tyco Fire Products Lp | System and method for evaluation of fluid flow in a piping system |
| WO2015071510A1 (en) * | 2013-11-12 | 2015-05-21 | Regasificadora Del Noreste, S.A. | System and method for optimising the management of a gas transportation network |
| CN106647329A (en) * | 2016-12-16 | 2017-05-10 | 中国电力科学研究院 | Heat supply network equivalent modeling method |
| CN109840305A (en) * | 2019-03-26 | 2019-06-04 | 中冶赛迪技术研究中心有限公司 | A kind of steam pipe network waterpower-thermal calculation method and system |
| CN110991092A (en) * | 2020-01-03 | 2020-04-10 | 中冶赛迪技术研究中心有限公司 | Complementary energy integrated utilization calculation method and system |
| CN111222211A (en) * | 2019-12-12 | 2020-06-02 | 合肥通用机械研究院有限公司 | Micro-channel unit CFD fluid-solid coupling structured hexahedral mesh division method |
| CN111680374A (en) * | 2020-05-21 | 2020-09-18 | 合肥泽众城市智能科技有限公司 | Method for checking and repairing topological relation of pipe network |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6885982B2 (en) * | 2000-06-27 | 2005-04-26 | Fluidigm Corporation | Object oriented microfluidic design method and system |
-
2020
- 2020-09-23 CN CN202011012525.9A patent/CN112035991B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8612189B2 (en) * | 2005-10-03 | 2013-12-17 | Tyco Fire Products Lp | System and method for evaluation of fluid flow in a piping system |
| WO2015071510A1 (en) * | 2013-11-12 | 2015-05-21 | Regasificadora Del Noreste, S.A. | System and method for optimising the management of a gas transportation network |
| CN106647329A (en) * | 2016-12-16 | 2017-05-10 | 中国电力科学研究院 | Heat supply network equivalent modeling method |
| CN109840305A (en) * | 2019-03-26 | 2019-06-04 | 中冶赛迪技术研究中心有限公司 | A kind of steam pipe network waterpower-thermal calculation method and system |
| CN111222211A (en) * | 2019-12-12 | 2020-06-02 | 合肥通用机械研究院有限公司 | Micro-channel unit CFD fluid-solid coupling structured hexahedral mesh division method |
| CN110991092A (en) * | 2020-01-03 | 2020-04-10 | 中冶赛迪技术研究中心有限公司 | Complementary energy integrated utilization calculation method and system |
| CN111680374A (en) * | 2020-05-21 | 2020-09-18 | 合肥泽众城市智能科技有限公司 | Method for checking and repairing topological relation of pipe network |
Non-Patent Citations (3)
| Title |
|---|
| 吴荣炜 ; 陈振宇 ; 尤天蓉 ; 王宽心 ; 吴玉成 ; 古勇 ; .蒸汽管网运行评价技术及应用.计算机与应用化学.2017,(第09期),全文. * |
| 张增刚 ; 李继志 ; 李永安 ; .城市蒸汽管网水力热力耦合系统及其应用.微计算机信息.2008,(第18期),全文. * |
| 高鲁锋 ; 郑海村 ; 朱启振 ; 孙德锋 ; .蒸汽管网水力热力耦合计算方法.中国科技论文.2013,(第08期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112035991A (en) | 2020-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hadidi et al. | Design and economic optimization of shell-and-tube heat exchangers using biogeography-based (BBO) algorithm | |
| Jung et al. | Irie: Scalable and robust influence maximization in social networks | |
| Tutar et al. | Numerical analysis of fluid flow and heat transfer characteristics in three-dimensional plate fin-and-tube heat exchangers | |
| CN109840305B (en) | Steam pipe network hydraulic-thermal calculation method and system | |
| Turgut | Thermal and economical optimization of a shell and tube evaporator using hybrid backtracking search—sine–cosine algorithm | |
| CN107871058B (en) | Load flow calculation method, device, equipment and storage medium of electric heating combined system | |
| Chandrakant et al. | Numerical simulation of flow-through heat exchanger having helical flow passage using high order accurate solution dependent weighted least square based gradient calculations | |
| Qian et al. | An ant colony system based on the Physarum network | |
| CN112035991B (en) | Steam optimization calculation method and system based on pipe network conveying path | |
| Branisavljevic et al. | Fuzzy approach in the uncertainty analysis of the water distribution network of Becej | |
| Mastrippolito et al. | Kriging metamodels-based multi-objective shape optimization applied to a multi-scale heat exchanger | |
| Peng et al. | Pressure sensor placement in water supply network based on graph neural network clustering method | |
| Wang et al. | Community discovery algorithm of complex network attention model | |
| Villemonteix et al. | Global optimization of expensive-to-evaluate functions: an empirical comparison of two sampling criteria | |
| CN115982918B (en) | Resistance characteristic identification method, device, equipment and storage medium | |
| Puliyeri et al. | Effect of cylinder arrangement on fluid flow and heat transfer characteristics past four elliptic cylinders | |
| Chen et al. | Towards Self-supervised Learning on Graphs with Heterophily | |
| Tricarico et al. | Economic level of reliability for the rehabilitation of hydraulic networks | |
| Liu et al. | Link prediction in aligned heterogeneous networks | |
| Aho et al. | Non-perennial stream networks as directed acyclic graphs: The R-package streamDAG | |
| Banasiak et al. | Generalized network transport and Euler-Hille formula | |
| Liu et al. | Signal-to-noise ratio analysis on saw-tooth vortex generator in vehicle radiators | |
| Conté et al. | Numerical investigation of forced fluid flow and heat transfer from conically coiled pipes | |
| Yongxin et al. | Pipe friction parameters identification method based on Moore-Penrose pseudo-inverse solution | |
| Zimmer et al. | The DLR ThermoFluid Stream Library |
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 |