CN112035991A - 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
- CN112035991A CN112035991A CN202011012525.9A CN202011012525A CN112035991A CN 112035991 A CN112035991 A CN 112035991A CN 202011012525 A CN202011012525 A CN 202011012525A CN 112035991 A CN112035991 A CN 112035991A
- Authority
- CN
- China
- Prior art keywords
- pipe network
- steam
- path
- matrix
- pipe
- 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 43
- 238000005457 optimization Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 51
- 239000012530 fluid Substances 0.000 claims abstract description 48
- 238000009833 condensation Methods 0.000 claims abstract description 17
- 230000008901 benefit Effects 0.000 claims abstract description 16
- 230000005494 condensation Effects 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 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
- 238000010586 diagram Methods 0.000 description 8
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000013498 data listing Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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: the physical information of the steam pipe network is displayed in a topological graph mode, and then the pipe network information is displayed in a matrix mode through a connection mode between pipelines and nodes. In the pipe network connection matrix, all conveying paths from a gas source point to a 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 calculating the condensation loss of each pipeline through heat transfer. The benefits of steam utilization are maximized by adjusting the amount of steam at different sources and users in combination with the benefits of steam usage and the losses of steam condensation. The method can calculate each conveying path and loss between the energy medium gas source point and the user point more quickly and accurately, and provides 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 utilization of steam can not be left in the conveying process of a pipe network, and the energy medium pipe network in the steel enterprise is usually a complex pipe network with multiple steam sources, multiple users and long span. Due to the lack of metering inside the pipe network, it is often difficult to find out the transport process and the transport path of the fluid in the pipe network. Especially for a steam pipe network, a large amount of condensation loss is often accompanied in the steam conveying process, and the condensation loss caused by different conveying paths is also different, which makes optimal scheduling of steam very difficult. Therefore, the calculation and steam optimization method and system based on the calculation and the calculation which can quickly and accurately find the fluid conveying path between each gas source and the user are very important for energy management and steam optimization of the iron and steel enterprises.
In the prior art, a method for finding a path is adopted by traversing all nodes of a pipe network by taking a certain member of the pipe network as a starting point and then extracting all non-coincident paths. In the existing traversing method, when the component is traversed, the current path is stopped from being traversed, and a plurality of repeated paths which are inconsistent with the actual flow exist in the path, and the path needs to be further screened and rejected; after finishing the preliminary screening and elimination, the node attributes in each path are required to be: and (4) further judging the starting point (Start) or the End point (End) to finally obtain the effective path formed by the effective nodes. In such a traversal method, a large number of invalid paths are often found, which not only wastes a large amount of computing resources, but also lengthens the running time of the whole program.
The following describes the disadvantages of this traversal method with reference to specific diagrams:
as shown in fig. 1, in the above-mentioned steam pipe network, the dots represent steam source points, and the squares represent steam user points, when the path finding is performed by traversing, the following situations may occur:
from the steam source point (Start) to the steam source point (Start) in the path found in fig. 2, such a path is practically non-existent because steam cannot enter the steam source point. The method for traversing can be used for judging the attribute of each node in the path: start point (Start) and user point (End) to cull such invalid paths. But the process of traversing searching and judging rejection consumes more operation time.
Alternatively, it is possible to find the nodes from (i) to (i) as shown in FIG. 3 by traversing
Also, only one source point (Start) and one user point (End) exist in the path. If a path is found by traversal, this is a valid path, but in practical situations, the arrows in the path may be in the opposite direction to the fluid flow in the pipe, resulting in the path being virtually non-existent or invalid.
In some existing steam pipe network calculation methods, the amount of condensed water in the steam pipe network can also be obtained through hydraulic-thermal coupling calculation, for example, after an incidence matrix is constructed, an admittance matrix between flow and impedance is written by using a flow node equation set and the relation between pipe section flow and pressure drop, and hydraulic calculation is performed through a method of matrix multiplication and equation solution; writing the relation between the flow and the temperature drop of the pipe section by the same method, writing an admittance matrix between the flow and the thermal resistance, and performing thermodynamic calculation by a method of matrix multiplication and equation solution; and continuously iterating the process to obtain the steam condensate amount by forming hydraulic-thermal coupling calculation, wherein the calculation amount of the hydraulic-thermal coupling calculation is very large, and the calculation time required to be consumed is relatively long.
In the steam optimization scheduling process, the flow rates of a steam source and a user point need to be adjusted continuously, and if hydraulic-thermal iteration coupling calculation needs to be performed after each adjustment, the whole optimization process is very difficult due to the fact that a large amount of calculation resources and calculation time are consumed, and real-time change of a steam pipe network cannot be responded.
Disclosure of Invention
In view of the above, the present invention provides a method and a system for calculating steam optimization based on a pipe network transportation path. Physical information of the energy medium pipe network is displayed in a topological graph mode, and then the pipe network information is displayed in a matrix mode through a connection mode between the pipelines and the nodes. In the pipe network connection matrix, all conveying paths from a gas source point to a 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 calculating the condensation loss of each pipeline through heat transfer. The benefits of steam utilization and the losses of steam condensation are combined, and the steam quantity of different steam sources and user points is adjusted to maximize the benefits of steam utilization. The method can calculate each conveying path and loss between the energy medium gas source point and the user point more quickly and accurately, provides reference for steam scheduling, and greatly improves the efficiency and the fineness of energy management of iron and steel enterprises.
In order to achieve the purpose, the invention provides the following technical scheme:
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 geographical position of the steam pipe network, and expressing pipe network information in the pipe network topological graph in a matrix mode according to the connection mode of pipelines and nodes, namely a pipe network information matrix;
s2: obtaining the flowing condition of the fluid in the pipeline through hydraulic calculation, and judging the flowing direction of the fluid medium of each pipeline in the pipe network;
s3: forming a new pipe network incidence matrix by combining the pipe network information matrix with the flow direction;
s4: finding all paths from the air source point to the user point by utilizing the actual fluid flow direction through a pipe network incidence 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, and the resistance coefficient of each path is calculated and used for calculating the pressure drop and the flow of the fluid on the path.
Optionally, in S1, when the pipe network information matrix is constructed:
firstly, numbering each node and each pipe section point;
the node and the pipe section point numbers form a two-dimensional matrix, and the number is represented by 1 if the pipe section and the node have a connection relation, and is represented by 0 if the pipe section and the node have no connection relation.
Optionally, in S3, when the pipe network information matrix is combined with the actual fluid flow direction to form a new pipe network association matrix:
performing 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 segment information matrix are labeled "1" and "-1" according to their relationship to flow in and out.
Optionally, in S4, when finding the flow path by using the correlation matrix:
each pipe section can only correspond to one inlet and one outlet, but one node can be communicated with the multiple pipe sections; by utilizing the property and combining the principle of 'one vertical and one horizontal and different sign connection', the whole path from the steam source to the user point can be quickly and accurately found in the two-dimensional incidence matrix by fully utilizing the actual fluid flow direction and the expression form in the incidence matrix;
optionally, in S5, when the pressure drop and the flow rate are calculated by using the resistance coefficient of each path:
by a physical formula that the multiplication of the flow rate and the resistance coefficient is equal to the pressure drop, and the pressure difference between the air source and the user is constant, the pressure drop equation on each path can be listed and the flow rate of each path can be solved:
ΔP＝∑ζi1Q1＝∑ζi2Q2＝…＝∑ζi3Qn
Q＝Q1+Q2+…+Qn
and calculating the condensation loss in the steam conveying process on each path, and optimizing the steam conveying scheme by combining the benefit of the steam, so that the overall use benefit of the steam is maximized.
An objective function: max { BenefitSteam generating device-CostCondensation of}。
A steam optimization computing system based on a pipe network transportation path, comprising:
the acquisition module is used for acquiring real-time state data of fluid states of all gas sources and user nodes in the fluid pipe network and physical parameter data for expressing pipe section structures, physical properties of the fluid pipe network and connection modes between pipe sections;
a data processing module for processing the acquired data,
the output module is used for outputting a data processing result;
the acquisition module, the data processing module and the output module are connected in sequence.
Optionally, the system further includes:
the display module is used for visualizing all paths between all air sources and user points and displaying all flow information on all paths in real time;
a data list module;
wherein, the data list module includes:
a real-time data list for storing real-time status data;
a physical parameter list for storing physical parameter data;
the result data list is used for storing all calculated results, including flow information on each path from the air source point to the user point and on the path;
and displaying the flow and pressure information measured by the metering instruments at all the air sources and the users, all the paths between the air sources and the users and the flow corresponding to each path through a display module.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the method.
An electronic terminal, comprising: a processor and a memory;
the memory is for storing a computer program and the processor is for executing the computer program stored by the memory to cause the terminal to perform 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 of the user point is adjusted each time, the change of the flow and the condensation amount on the corresponding path is only required to be calculated, and the calculation of the hydraulic-thermal coupling of the whole pipe network is not required. Therefore, the calculation resources required by the optimization process can be greatly reduced, and the calculation program can better respond to the real-time change of the steam pipe network.
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 objectives and other advantages of the invention may be realized and attained by the means of the instrumentalities and combinations particularly pointed out hereinafter.
Drawings
For the purposes of promoting a better understanding of the objects, aspects and advantages of the invention, reference will now be made to the following detailed description taken in conjunction with the accompanying drawings in which:
FIG. 1 is a diagram of an example pipe network;
FIG. 2 is a diagram of example paths one;
FIG. 3 is a second example of a path;
FIG. 4 is a schematic view of a multi-gas source multi-user complex fluid pipe network topology;
FIG. 5 is a schematic diagram of a pipe network information matrix;
FIG. 6 is a pipe network incidence matrix;
FIG. 7 is a flow path from a pipe network incidence matrix source point (first) to a user point (nine);
FIG. 8 is a schematic diagram of a path from a gas source point (c) to a user point (c);
FIG. 9 is a schematic diagram of a pipe network topology graph path;
FIG. 10 is an illustration of "connected by different sign" pipe network topology;
FIG. 11 is a diagram illustrating a "different sign connected" pipe network correlation matrix.
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention in a schematic way, and the features in the following embodiments and examples may be combined with each other without conflict.
Wherein the showings are for the purpose of illustrating the invention only and not for the purpose of limiting the same, and in which there is shown by way of illustration only and not in the drawings in which there is no intention to limit the invention thereto; to better illustrate the embodiments of the present invention, some parts of the drawings may be omitted, enlarged or reduced, and do not represent the size of an actual product; it will be understood 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 numerals in the drawings of the embodiments of the present invention correspond to the same or similar components; in the description of the present invention, it should be understood that if there is an orientation or positional relationship indicated by terms such as "upper", "lower", "left", "right", "front", "rear", etc., based on the orientation or positional relationship shown in the drawings, it is only for convenience of description and simplification of description, but it is not an indication or suggestion that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and therefore, the terms describing the positional relationship in the drawings are only used for illustrative purposes, and are not to be construed as limiting the present invention, and the specific meaning of the terms may be understood by those skilled in the art according to specific situations.
The invention will now be described with reference to a certain fluid network shown in fig. 4:
the dots shown in the above figures represent the fluid air supply, the squares represent the fluid user, and the black lines represent the fluid conduits. Each node and each segment of pipe are first numbered according to the complex network comprising a plurality of rings on such a multi-source, multi-user network. The numbers in the circles are the numbers of the nodes, such as 'r'; the individual numbers indicate the number of pipe segments, e.g. '1', and the numbering of nodes and pipe segments can be arranged in any order, but it must be ensured that each number is independent and consecutive.
In an actual large-scale pipe network, only an air source point or a user is often provided with a metering instrument, and the inside of the pipe network is not provided with the metering instrument, such as: (nine) in the figure
Node has meter and (r) () r (r)
"where there is no meter. Due to the lack of metering inside the pipe network, the flow paths of the fluid in the pipe network and the flow rate of each path are difficult to be known, which brings great difficulty to the energy fine management of enterprises.
According to the method, firstly, all paths from a steam source point to a user point and the flow of each path can be rapidly calculated through a method for calculating the conveying path based on the actual fluid flow direction, so that the flow condition of fluid in a pipe network is known, the refinement degree of enterprise energy management is improved, and the specific implementation mode is as follows:
according to the serial numbers and the 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 serial numbers and the connection modes, the ordinate represents the serial numbers of the nodes, the abscissa represents the serial numbers of the pipe sections, all data in a program are generated by program calculation, so that the number in the matrix needs to be added with 1 to correspond to the serial number in the pipe network, for example, the number '0' of the ordinate of the matrix represents the node with the serial number 'r'. The number "1" inside the matrix represents that the node is connected to the pipe, and the number "0" represents that the node is not connected to the pipe. For example, two numbers "1" in column 1 of the figure indicate that the node "r" and the node "r" are connected by the pipe segment "1".
The direction of flow of the fluid in all pipe sections can be derived from the hydraulic calculation of the pipe, i.e. the fluid always flows from the point where the pressure is high to the point where the pressure is low, as indicated by the arrows in fig. 4. Combining the fluid flow direction with the pipe network information matrix, a pipe network correlation matrix can be obtained, as shown in fig. 6, the correlation matrix is represented by the numbers "-1" and "1" to indicate the flow direction of the fluid in the pipe, "1" indicates inflow, "1" indicates outflow, and, for example, the number "1" at (0,0) indicates that the fluid flows from the node "r" into the pipe "1".
Then, all paths from the air source point to the user point can be found in the incidence matrix by using a traversal 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. This property can be exploited to quickly and accurately find all paths from the steam source point to the user point in a two-dimensional correlation matrix by path computation methods that take into account the actual fluid flow. The following is a detailed description:
if we need to find all the flow paths from the source point 'r' to the user point 'ninu' in fig. 4, the expression in the correlation matrix is as shown in fig. 7, and the flow path is found by connecting the number '1' of the source point 'r' with all other numbers according to the principle of 'connecting one vertical to one horizontal and connecting different signs' from the number '1' of the source point until the principle of 'connecting one vertical to one horizontal and connecting different signs' is not satisfied. For example, starting from the number "1" starting at the coordinate (0,0), all the odd-numbered numbers "-1" of the column where the column is vertically connected first, and then "-1" at the coordinate (0,3) can be reached; and then, all the different-sign numbers of the row in which the "-1" at the position (0,3) is transversely connected are '1', three '1' at the positions of coordinates (2,3), (3,3) and (4,3) can be respectively found, and three paths outwards from the node 'r' are represented. Repeating the steps, wherein the longitudinal connection and the transverse connection are alternately performed every time, each number is connected with all the different-sign numbers in the same row or column, and the process is ended when a certain path cannot be connected any more. If the ' 1 ' at the coordinates (2,3) can be longitudinally connected with the ' -1 ' at the positions (2,2), when the ' -1 ' at the positions (2,2) needs to be transversely connected again, all numbers of the row are ' 0 ', the connection is finished if the connection requirements are not met, and all nodes passing under the path are recorded, namely, the sequential nodes ' are ' first, fourth and third, and the ' path from the air source point ' first ' to the user point ' third '.
The above is an explanation of the principle of "connecting one vertical line to one horizontal line and connecting different signs", taking all the paths from the gas source searching point "i" to the user point "ninu" as an example, two feasible routes can be found by the above method, as shown in fig. 8:
as shown in fig. 8, two paths can be found, different places of the two paths are respectively indicated by a dotted arrow and an implementation arrow, a path one passes through a correlation matrix with coordinates of [ (0,0), (0,3), (3,3), (3,4), (5,4), (5,7), (8,7), (8,8) ], and a corresponding path node is numbered [ "first, fourth, fifth, eight, and nine ]; coordinates in the correlation matrix where another path passes through in sequence are [ (0,0), (0,3), (4,3), (4,6), (7,6), (7,7), (8,7), (8,8) ], and then the corresponding path node is numbered [ "phi, and ninc ]. The representation of the two paths on the pipe network topology is shown in fig. 9.
Two flow paths found in the incidence matrix (i, r, c, b, c) and (i, r, c, b) correspond to the path on the pipe network topology diagram.
Finally, all the paths from the gas source finding point "c" to the user point "c" are taken as an example to explain why the principle of "connected by different signs" must be satisfied. Then a path from gas source point to 'r' to user point 'is [' r ',' c, r 'and r' respectively,
⑧、⑨”]The route of the method on the pipe network topological graph is shown in figure 10.
In the path (r, c, r, c,
⑧、⑨”]In the above, the flow direction of the pipe section "11" is calculated by the hydraulic power of the pipe network, so that the fluid in the pipe section can only flow from the pipe sectionNode 'r' to node
And not in the opposite direction. The expression of such a topology map in a pipe network incidence matrix is shown in fig. 11.
After the path finds "-1" of the coordinates (11,10), the number at the coordinates (10,10) is "-1", so that the principle of "different sign connection" is not satisfied, that is, the path does not exist, and the physical meaning of the principle of "different sign connection" in the matrix is the pipe section flow direction in the corresponding pipe network topological graph, and any path finding process cannot violate the actual flow direction of the fluid in the pipe section, otherwise, the physical meaning corresponding to the fluid in the pipe section is lost.
Therefore, it can be obtained that all paths from the air source point "c to the user point" c "have only two paths, and the corresponding node sequences are respectively: (iii) phi, phi "]And [' r, c, b, c [ ] "]. Since the path calculation based on the actual fluid flow direction can be performed starting from any desired steam source point in the path calculation process, all flow paths from all desired steam source points to the user points can be found. For example, firstly, the traversal finds out the air source point (r) and sequentially finds out all the user points (c, c,
"all paths; if necessary, the source point of the gas is changed, and all the user points are found out sequentially by the same method,
"all paths. By analogy, each required flow path from the steam source point to the user point can be found. The method is also suitable for more complicated large-scale fluid pipe networks.
And finally, calculating the flow of each path by using the relationship between the resistance and the flow through the following formula:
ΔP＝∑ζi1Q1＝∑ζi2Q2＝…＝∑ζi3Qn
Q＝Q1+Q2+…+Qn
where Δ P is the pressure difference between the source point and the user point, ζi1Represents the drag coefficient, Σ ζ, of the "i" th pipe section on path 1i1The drag coefficients of all the pipe sections making up the path 1 are summed, where the drag coefficient of each pipe section is a fixed parameter determined by the physical parameters of the pipe section, such as diameter, roughness, length, etc. Q1Represents the fluid flow of path 1; sigma zetai1Q1＝∑ζi2Q2＝…＝∑ζi2QnIdentity means: since the pressure difference ap between the source point and the user point is constant (measured by the pressure meter), the multipliers of the flow rates of each path and the sum of the resistance coefficients are equal, so that the flow rate of the path with a large resistance coefficient is small, and the flow rate of the path with a small resistance coefficient is large. (for easy understanding, the relationship between voltage, current and resistance can be analogized in that the voltage between two points is constant, the current on the circuit with small resistance is larger, and the current on the circuit with large resistance is smaller); q represents the total flow from the source to the user point, i.e., equal to the sum of all flows for each path.
Finally, after all steam paths from the steam source point to the user point and the flow are obtained, the corresponding amount of condensed water on each path can be obtained through heat transfer calculation, and the formula is as follows:
wherein P isiThe heat exchange amount with the outside in 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), and r is the latent heat of the steam. The calculated Condensate is the sum of the amount of condensed water on all paths from a certain steam source to the user point.
After the amount of condensed water in the conveying process is calculated, the steam conveying scheme is optimized in combination with the benefit of steam, so that the overall use benefit of the steam is maximized, and the objective function is as follows: max { BenefitSteam generating device-CostCondensation of}. Wherein BenefitSteam generating deviceCost for the benefit of steam for heating or production after it reaches the user's siteCondensation ofRefers to the loss of steam due to condensation during transport. By continuously adjusting the steam amount of the steam source and the steam user and calculating the condensation loss on each conveying path, new Benefit can be obtainedSteam generating deviceAnd CostCondensation ofAnd ending the steam optimization process until a larger overall profit value meeting the demand is obtained. Therefore, 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 calculation of steam optimization of a pipe network conveying path, and after the flow of a steam source and the flow of a user point is adjusted each time, the change of the flow and the condensation on the corresponding path is calculated, and the calculation of the hydraulic-thermal coupling 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 better respond to the real-time change of the steam pipe network.
Hardware device
In order to accomplish the above inventive functions, a set of computing system needs to be formed by matching corresponding hardware devices, including:
the acquisition module is used for respectively acquiring real-time state data of fluid states of all gas sources and user nodes (namely nodes with metering) in the fluid pipe network and physical parameter data for expressing pipe section structures, physical properties (length, diameter and the like) of the fluid pipe network and connection modes between pipe sections;
a data processing module for processing the acquired data,
the output module is used for outputting a data processing result;
optionally, the method further includes:
and the display module is used for visualizing all paths between all the air sources and the user points and displaying all the flow information on all the paths in real time.
A data listing module, the data listing module comprising:
a real-time data list for storing real-time status data;
a physical parameter list for storing physical parameter data;
a result data list for storing all calculated results, such as information of each path from the air source point to the user point, flow rate on the path, and the like;
and displaying the information such as the flow and the pressure measured by the metering instruments at all the air sources and the users, all the paths between the air sources and the users and the flow corresponding to each path through a display module.
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 one of the above.
The present invention also provides an electronic terminal, comprising: a processor and a memory;
the memory is adapted to store a computer program and the processor is adapted to execute the computer program stored by the memory to cause the terminal to perform the method as defined in any one of the above.
The computer-readable storage medium in the present embodiment can be understood by those skilled in the art as follows: all or part of the steps for implementing the above method embodiments may be performed by hardware associated with a computer program. The aforementioned computer program may be stored in a computer readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
The electronic terminal provided by the embodiment comprises a processor, a memory, a transceiver and a communication interface, wherein the memory and the communication interface are connected with the processor and the transceiver and are used for completing mutual communication, the memory is used for storing a computer program, the communication interface is used for carrying out communication, and the processor and the transceiver are used for operating the computer program so that the electronic terminal can execute the steps of the method.
In this embodiment, the Memory may include a Random Access Memory (RAM) or may further include a non-volatile Memory (non-volatile Memory), such as at least one disk Memory.
The processor may be a general-purpose processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the Integrated Circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, or a discrete hardware component.
Finally, the above embodiments are only intended to illustrate the technical solutions of the present invention and not to limit the present invention, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions, and all of them should be covered by the claims of the present invention.
Claims (9)
1. A steam optimization calculation method based on a pipe network conveying path is characterized by comprising the following steps: the method comprises the following steps:
s1: constructing a pipe network topological graph through the geographical position of the steam pipe network, and expressing pipe network information in the pipe network topological graph in a matrix mode according to the connection mode of pipelines and nodes, namely a pipe network information matrix;
s2: obtaining the flowing condition of the fluid in the pipeline through hydraulic calculation, and judging the flowing direction of the fluid medium of each pipeline in the pipe network;
s3: forming a new pipe network incidence matrix by combining the pipe network information matrix with the flow direction;
s4: finding all paths from the air source point to the user point by utilizing the actual fluid flow direction through a pipe network incidence 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, and the resistance coefficient of each path is calculated and used for 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 as claimed in claim 1, wherein: in S1, when the pipe network information matrix is constructed:
firstly, numbering each node and each pipe section point;
the node and the pipe section point numbers form a two-dimensional matrix, and the number is represented by 1 if the pipe section and the node have a connection relation, and is represented by 0 if the pipe section and the node have no connection relation.
3. The steam optimization calculation method based on the pipe network conveying path as claimed in claim 1, wherein: in S3, when the pipe network information matrix is combined with the actual fluid flow direction to form a new pipe network association matrix:
performing 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 segment information matrix are labeled "1" and "-1" according to their relationship to flow in and out.
4. The steam optimization calculation method based on the pipe network conveying path as claimed in claim 1, wherein: in S4, when finding a flow path using the correlation matrix:
each pipe section can only correspond to one inlet and one outlet, but one node can be communicated with the multiple pipe sections; by utilizing the property and combining the principle of 'one vertical and one horizontal and different sign connection', the whole path from the steam source to the user point can be quickly and accurately found in the two-dimensional incidence matrix by fully utilizing the actual fluid flow direction and the expression form in the incidence matrix.
5. The steam optimization calculation method based on the pipe network conveying path as claimed in claim 1, wherein: in S5, when calculating the pressure drop and the flow rate by using the resistance coefficient of each path:
by a physical formula that the multiplication of the flow rate and the resistance coefficient is equal to the pressure drop, and the pressure difference between the air source and the user is constant, the pressure drop equation on each path can be listed and the flow rate of each path can be solved:
ΔP＝∑ζi1Q1＝∑ζi2Q2＝…＝∑ζi3Qn
Q＝Q1+Q2+…+Qn
calculating the condensation loss in the steam conveying process on each path, and optimizing the steam conveying scheme by combining the benefit of the steam, so that the overall use benefit of the steam is maximized;
an objective function: max { BenefitSteam generating device-CostCondensation of}。
6. A steam optimization computing system based on a pipe network conveying path is characterized in that: the system comprises:
the acquisition module is used for acquiring real-time state data of fluid states of all gas sources and user nodes in the fluid pipe network and physical parameter data for expressing pipe section structures, physical properties of the fluid pipe network and connection modes between pipe sections;
a data processing module for processing the acquired data,
the output module is used for outputting a data processing result;
the acquisition module, the data processing module and the output module are connected in sequence.
7. The system of claim 6, wherein the computing system is configured to optimize steam based on a pipe network transportation path: the system further comprises:
the display module is used for visualizing all paths between all air sources and user points and displaying all flow information on all paths in real time;
a data list module;
wherein, the data list module includes:
a real-time data list for storing real-time status data;
a physical parameter list for storing physical parameter data;
the result data list is used for storing all calculated results, including flow information on each path from the air source point to the user point and on the path;
and displaying the flow and pressure information measured by the metering instruments at all the air sources and the users, all the paths between the air sources and the users and the flow corresponding to each path through a display module.
8. 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 one of claims 1 to 5.
9. An electronic terminal, characterized by: the method comprises the following steps: a processor and a memory;
the memory is used for storing computer programs, and the processor is used for executing the computer programs stored by the memory so as to enable the terminal to execute the method of any one of 1-5.
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 true CN112035991A (en) | 2020-12-04 |
| CN112035991B 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115231504A (en) * | 2022-06-25 | 2022-10-25 | 山东鲁睿智能科技有限公司 | Double-seat valve array control method and system for reducing residues |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020108097A1 (en) * | 2000-06-27 | 2002-08-08 | Fluidigm Corporation | Object oriented microfluidic design method and system |
| 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 |
-
2020
- 2020-09-23 CN CN202011012525.9A patent/CN112035991B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020108097A1 (en) * | 2000-06-27 | 2002-08-08 | Fluidigm Corporation | Object oriented microfluidic design method and system |
| 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 |
|---|
| 吴荣炜;陈振宇;尤天蓉;王宽心;吴玉成;古勇;: "蒸汽管网运行评价技术及应用", 计算机与应用化学, no. 09, 28 September 2017 (2017-09-28) * |
| 张增刚;李继志;李永安;: "城市蒸汽管网水力热力耦合系统及其应用", 微计算机信息, no. 18, 25 June 2008 (2008-06-25) * |
| 高鲁锋;郑海村;朱启振;孙德锋;: "蒸汽管网水力热力耦合计算方法", 中国科技论文, no. 08, 15 August 2013 (2013-08-15) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115231504A (en) * | 2022-06-25 | 2022-10-25 | 山东鲁睿智能科技有限公司 | Double-seat valve array control method and system for reducing residues |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112035991B (en) | 2024-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Heat-source solutions to heat conduction in anisotropic media with application to pile and borehole ground heat exchangers | |
| Hajabdollahi et al. | Exergetic optimization of shell-and-tube heat exchangers using NSGA-II | |
| 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 | |
| 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 | |
| Li et al. | Study on thermal-hydraulic performance of printed circuit heat exchangers with supercritical methane based on machine learning methods | |
| CN104749999A (en) | Accurate guidance system for optimizing operation of turbo generator group cold end system of assembly wet-type cooling tower | |
| CN105024886B (en) | A kind of Fast W eb service QoS Forecasting Methodologies based on user metadata | |
| Yang et al. | Improvements on flow distribution and heat transfer performance of plate-fin heat exchangers by Qusai-S type header configuration | |
| CN112035991A (en) | Steam optimization calculation method and system based on pipe network conveying path | |
| CN107871058A (en) | Tidal current computing method, device, equipment and the storage medium of electric heating association system | |
| Macchitella et al. | Performance prediction of plate-finned tube heat exchangers for refrigeration: a review on modeling and optimization methods | |
| Qian et al. | An ant colony system based on the Physarum network | |
| Sridharan | Performance optimization of counter flow double pipe heat exchanger using grey relational analysis | |
| Kaya et al. | A new method to multi–objective optimization of shell and tube heat exchanger for waste heat recovery | |
| CN116227113A (en) | Solving method and device for thermodynamic system simulation, electronic equipment and medium | |
| JP7485404B2 (en) | Method and system for estimating hydraulic conditions in dynamic operation of steam heat supply networks | |
| CN114239199B (en) | Steam pipe network dynamic simulation method considering condensed water | |
| CN107016459A (en) | A kind of point-to-point shortest path computational methods based on network community message | |
| Pacheco-Vega | Simulation of compact heat exchangers using global regression and soft computing | |
| CN115081327A (en) | Partition wall type heat exchanger dynamic simulation model modeling method based on reinforcement learning | |
| Nasr et al. | Heat transfer coefficient and friction factor prediction of corrugated tubes combined with twisted tape inserts using artificial neural network | |
| Garcia et al. | Design of a numerical simulator for finned-tube heat exchangers with arbitrary circuitry | |
| Dong et al. | A novel balanced teaching-learning-based optimization algorithm for optimal design of high efficiency plate-fin heat exchanger | |
| Conté et al. | Numerical investigation of forced fluid flow and heat transfer from conically coiled pipes |
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 |