CN117113862A - Pipe network leakage identification model construction method and pipe network leakage identification method - Google Patents
Pipe network leakage identification model construction method and pipe network leakage identification method Download PDFInfo
- Publication number
- CN117113862A CN117113862A CN202310882367.XA CN202310882367A CN117113862A CN 117113862 A CN117113862 A CN 117113862A CN 202310882367 A CN202310882367 A CN 202310882367A CN 117113862 A CN117113862 A CN 117113862A
- Authority
- CN
- China
- Prior art keywords
- pipe network
- matrix
- leakage
- node
- pressure
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000010276 construction Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 143
- 238000011156 evaluation Methods 0.000 claims abstract description 48
- 238000004088 simulation Methods 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000012423 maintenance Methods 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 20
- 230000000295 complement effect Effects 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 12
- 230000000903 blocking effect Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 2
- 238000010438 heat treatment Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000005431 greenhouse gas Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
The application relates to a pipe network leakage identification model construction method and a pipe network leakage identification method, and belongs to the field of heat supply pipe network leakage identification. The pipe network leakage identification model construction method and the pipe network leakage identification method comprise the following steps: determining a basic incidence matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic incidence matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and topological connection relations among the loops; determining the simulation pressure of the node through a simulation algorithm based on the basic association matrix, the independent loop matrix and the flow parameters of the node; iteratively updating the flow parameter based on the gradient of the simulated pressure with respect to the flow parameter; and taking the flow parameters acquired when the preset iteration conditions are met as leakage evaluation parameters to obtain a pipe network leakage identification model.
Description
Technical Field
The application relates to the field of heat supply pipe network leakage identification, in particular to a pipe network leakage identification model construction method and a pipe network leakage identification method.
Background
With the increasing world population and the increasing urbanization process, the energy demands of urban areas are increasing. In order to cope with global climate change and reduce greenhouse gas emission, the aim of carbon neutralization is fulfilled, and the development and utilization of clean energy are accelerated in various countries, and meanwhile, the energy conservation and efficient utilization are encouraged. Under the background, the district heating system is used as an energy system for central heating, and generated hot water is distributed to all end users through a directly buried heating pipe network, so that the district heating system has the advantages of improving energy efficiency, reducing greenhouse gas emission, improving air quality and having better scale economy.
But the leakage condition can appear in the terminal branch pipe of heating supply network, leads to heating supply network heat loss, and the pipe network seepage mainly takes place on the branch pipe section before the heat exchange station, not only extravagant energy, can influence the income of enterprise moreover. There is therefore an urgent need to effectively identify heating network leaks.
At present, no effective solution is proposed to effectively distinguish the leakage of the heating network.
Disclosure of Invention
The embodiment of the application provides a pipe network leakage identification model construction method and a pipe network leakage identification method, which are used for at least solving the problem that the leakage of a heating pipe network cannot be effectively identified in the related technology.
In a first aspect, an embodiment of the present application provides a method for constructing a pipe network leakage identification model, including:
determining a basic incidence matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic incidence matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and topological connection relations among the loops;
determining the simulation pressure of the node through a simulation algorithm based on the basic association matrix, the independent loop matrix and the flow parameters of the node;
iteratively updating the flow parameter based on the gradient of the simulated pressure with respect to the flow parameter;
and taking the flow parameters acquired when the preset iteration conditions are met as leakage evaluation parameters to obtain a pipe network leakage identification model.
In one embodiment, the determining the basic association matrix and the independent loop matrix based on the actual pipe network topology includes:
constructing a virtual node on each tail end branch pipe in the actual pipe network topological structure, and constructing a virtual pipe section taking the virtual node as a starting point to obtain a complement pipe network topological structure;
and determining the basic association matrix and the independent loop matrix according to the layout relation of the pipe sections and the nodes in the full pipe network topological structure.
In one embodiment, the determining, by a simulation algorithm, the simulation pressure of the node based on the basic association matrix, the independent loop matrix, and the flow parameters of the node includes:
obtaining a blocking matrix according to a preset blocking rule by the basic association matrix and the independent loop matrix;
and determining the simulation pressure of the node based on the blocking matrix and the flow parameters of the node.
In one embodiment, the iteratively updating the flow parameter based on the gradient of the simulated pressure with respect to the flow parameter comprises:
determining a pressure error value according to the monitoring pressure of the node and the simulation pressure of the node;
and iteratively updating the flow parameter based on a gradient of the simulated pressure with respect to the flow parameter in response to the pressure error value being greater than a preset threshold.
In one embodiment, the obtaining the pipe network leakage identification model using the flow parameter obtained when the preset iteration condition is satisfied as the leakage evaluation parameter includes:
the meeting the preset iteration condition comprises the following steps: the error value of the monitoring pressure and the simulation pressure is smaller than or equal to the preset threshold value, or the iteration number is equal to the preset iteration number.
In a second aspect, an embodiment of the present application provides a pipe network leakage identification method, including:
acquiring a current pipe network topological structure and monitoring pressure of each node in a pipe network, and constructing a complement pipe network topological structure based on the current pipe network topological structure;
determining leakage evaluation parameters of each node based on the complement network management topological structure and the monitoring pressure by using the pipe network leakage identification model constructed by the method of any one of claims 1 to 5;
and identifying the leakage point according to the leakage evaluation parameter and a preset evaluation standard.
In one embodiment, the preset evaluation criteria includes a first interval and a second interval; the method further comprises the steps of:
determining a first maintenance strategy in response to the leakage evaluation parameter of the node being located in the first interval;
and determining a second maintenance strategy in response to the leakage evaluation parameter of the node being located in the second interval, wherein the priorities of the first maintenance strategy and the second maintenance strategy are different.
In a third aspect, an embodiment of the present application provides a pipe network leakage identification model construction apparatus, including:
the matrix determining module is used for determining a basic incidence matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic incidence matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and topological connection relation among the loops;
the simulation pressure determining module is used for determining the simulation pressure of the node through a simulation algorithm based on the basic incidence matrix, the independent loop matrix and the flow parameters of the node;
a flow parameter updating module for iteratively updating the flow parameter based on the gradient of the simulation pressure with respect to the flow parameter;
and the identification module is used for taking the flow parameter acquired when the preset iteration condition is met as a leakage evaluation parameter to acquire a pipe network leakage identification result.
In a fourth aspect, an embodiment of the present application provides a pipe network leakage identification device, including:
the construction module is used for acquiring the current pipe network topological structure and the monitoring pressure of each node in the pipe network, and constructing a complement pipe network topological structure based on the current pipe network topological structure;
the determining module is used for determining leakage evaluation parameters of each node based on the complement network management topological structure and the monitoring pressure through the pipe network leakage identification model constructed by the method of any one of claims 1-5;
and the identification module is used for identifying the leakage points according to the leakage evaluation parameters and preset evaluation standards.
In a fifth aspect, an embodiment of the present application provides a computer readable storage medium having a program stored thereon, where the program, when executed by a processor, implements the method for constructing a pipe network leakage identification model and the method for identifying pipe network leakage according to any one of the above.
The pipe network leakage identification model construction method and the pipe network leakage identification method provided by the embodiment of the application have at least the following technical effects.
The application constructs a basic incidence matrix and an independent loop matrix based on a pipe network topological structure, and calculates simulation pressure through a heating pipe network hydraulic simulation algorithm based on the basic incidence matrix, the independent loop matrix and node leakage flow parameters. And iteratively updating the leakage flow of each branch pipe section according to the gradient descent algorithm. The simulation pressure at each heat exchange station is consistent with the measured pressure, leakage flow rate at the end of iteration is used as leakage evaluation parameters, and parameters such as an impedance matrix of leakage nodes are output, so that the problem of leakage identification of the tail end branch pipes of the heat supply pipe network is solved, the leakage identification result is more accurate and comprehensive, and decisions and supports are provided for relevant heat supply pipe network management staff to maintain the heat supply pipe network.
The details of one or more embodiments of the application are set forth in the accompanying drawings and the description below to provide a more thorough understanding of the other features, objects, and advantages of the application.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
FIG. 1 is a flow chart illustrating a method of constructing a pipe network leak identification model according to an exemplary embodiment;
FIG. 2 is an association matrix and independent loop matrix validation flow diagram shown in accordance with an exemplary embodiment;
FIG. 3 is a diagram of a completed pipe network topology according to an exemplary embodiment;
FIG. 4 is a simulated pressure calculation flow chart shown according to an exemplary embodiment;
FIG. 5 is a flowchart illustrating iteratively updating flow parameters according to an exemplary embodiment;
FIG. 6 is a block diagram of a pipe network leak identification method according to an embodiment of the application;
FIG. 7 is a block diagram of a pipe network leakage identification model building apparatus according to an embodiment of the present application;
FIG. 8 is a block diagram of a pipe network leak identification device according to an embodiment of the application.
Detailed Description
The present application will be described and illustrated with reference to the accompanying drawings and examples in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application. All other embodiments, which can be made by a person of ordinary skill in the art based on the embodiments provided by the present application without making any inventive effort, are intended to fall within the scope of the present application.
It is apparent that the drawings in the following description are only some examples or embodiments of the present application, and it is possible for those of ordinary skill in the art to apply the present application to other similar situations according to these drawings without inventive effort. Moreover, it should be appreciated that while such a development effort might be complex and lengthy, it would nevertheless be a routine undertaking of design, fabrication, or manufacture for those of ordinary skill having the benefit of this disclosure, and thus should not be construed as having the benefit of this disclosure.
Reference in the specification to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment may be included in at least one embodiment of the application. The appearances of such phrases in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. It is to be expressly and implicitly understood by those of ordinary skill in the art that the described embodiments of the application can be combined with other embodiments without conflict.
Unless defined otherwise, technical or scientific terms used herein should be given the ordinary meaning as understood by one of ordinary skill in the art to which this application belongs. The terms "a," "an," "the," and similar referents in the context of the application are not to be construed as limiting the quantity, but rather as singular or plural. The terms "comprising," "including," "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, article, or apparatus that comprises a list of steps or modules (elements) is not limited to only those steps or elements but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus. The terms "connected," "coupled," and the like in connection with the present application are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The term "plurality" as used herein means two or more. "and/or" describes an association relationship of an association object, meaning that there may be three relationships, e.g., "a and/or B" may mean: a exists alone, A and B exist together, and B exists alone. The character "/" generally indicates that the context-dependent object is an "or" relationship. The terms "first," "second," "third," and the like, as used herein, are merely distinguishing between similar objects and not representing a particular ordering of objects.
Based on the above situation, the embodiment of the application provides a pipe network leakage identification model construction method and a pipe network leakage identification method.
First aspect
The embodiment of the application provides a construction method of a pipe network leakage identification model.
FIG. 1 is a flow chart illustrating a method of constructing a pipe network leak identification model according to an exemplary embodiment. As shown in fig. 1, the method includes:
step S101, determining a basic association matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic association matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and the topological connection relation among the loops. Wherein the topology of the heating network comprises leaky nodes.
In one example, the determination of the association matrix and the independent loop matrix takes place in the following manner. Fig. 2 is a flowchart illustrating an association matrix and an independent loop matrix acknowledgement procedure according to an exemplary embodiment, and as shown in fig. 2, step S101 specifically includes:
and S1011, constructing a virtual node on each end branch pipe in the actual pipe network topology structure, and constructing a virtual pipe section taking the virtual node as a starting point to obtain the completed pipe network topology structure. Fig. 3 is a diagram of a topology of a completed pipe network according to an exemplary embodiment, and reference is made to fig. 3, in which pipe segments (8) and (9) are virtual pipe segments. In this example, it is assumed that the virtual nodes on all the end branch pipes are leakage points, and further, whether each node (including the real node and the virtual node) belongs to the real leakage point is identified through the subsequent steps, so that omission is avoided.
Step S1012, determining a basic association matrix and an independent loop matrix according to the layout relation of pipe sections and nodes in the full pipe network topology structure. Table 1 is a basic association matrix structure shown according to an exemplary embodiment, as shown in table 1.
TABLE 1 basic incidence matrix Structure
Table 2 is an exemplary diagram of an independent loop matrix shown in accordance with an exemplary embodiment. As shown in table 2.
TABLE 2 independent loop matrix structure
Wherein 0, 1, -1 represents the connection and direction relation between each node and the pipe section, and the reference diagram is a numerical value of 0 when the pipe section and the node have no connection relation; when the starting point of the pipe section is the current node, the numerical value is-1; when the end point of the pipe section is the current node, the value is 1.
In this embodiment, all end nodes are assumed to be leaky nodes, added to the topology of the heating network, via step S1011. In step S1012, the basic association matrix and the independent loop matrix about the nodes and the pipe sections in the heating pipe network are resolved through the topology structure, and the matrix is partitioned, so that the linearization calculation of the subsequent equation set is facilitated. The identification accuracy is improved, and missed detection is avoided.
With continued reference to fig. 1, step S102 is performed after step S101.
Step S102, based on the basic association matrix, the independent loop matrix and the flow parameters of the nodes, the simulation pressure of the nodes is determined through a simulation algorithm. Wherein the simulation algorithm used is a basic loop method.
In one example, fig. 4 is a simulation pressure calculation flowchart, which is shown in an exemplary embodiment, and as shown in fig. 4, determining the simulation pressure of the node in step S102 specifically includes:
step S1021, obtaining a block matrix according to a preset block dividing rule by the basic association matrix and the independent loop matrix.
Wherein, a basic association matrix A is constructed, and the blocking matrix form is thatWherein O is 0 matrix, I is identity matrix, matrix A 1 The block matrix form is A 1 ＝[A 1，t A 1，r ]. Where t represents the branch portion in the heating pipe network, r represents the remaining branch portion, and referring to fig. 3, the pipe section 7 is the remaining branch portion, and the other pipe sections are branch portions.
Constructing an independent loop matrix C, wherein the block matrix form of the independent loop matrix C is C= [ C f O]Matrix C f The block matrix form is C f ＝[C f，1 C f，2 ]。
Taking Table 3 as an example, assume that a heating network comprising "virtual pipe sections" comprises n pipe sections, where n is 1 Strip true pipe section, n 2 The bar is a virtual pipe segment. The heat supply network comprises m nodes, wherein m-n is included 2 True nodes, n 2 And virtual nodes.
TABLE 3 matrix Specification Table
Step S1022, determining the simulation pressure of the node based on the block matrix and the flow parameters of the node.
Wherein update-based leaky node traffic Q 2 The method comprises the following steps of calculating a pipe network node simulation pressure matrix P by utilizing a heat supply pipe network hydraulic simulation algorithm:
(1) According to the flow matrix Q of the nodes of the heat supply pipe network 1 (leaky node traffic is 0) and leaky node traffic default value matrix Q 2 Node flow matrix Q for constructing heating network (only including leakage nodes and no virtual pipe sections) 3 . Wherein Q is 1 Is the actual measurement data; q (Q) 2 For the preset value, the values can be preset to be all 0.
(2) Assuming an initial residual pipe section flow matrix G 1，x Determining branch pipe section flow matrix based on node mass balance equation
(3) Iteratively updating the heating pipe network pipe section flow G according to the loop pressure balance equation 1 Wherein M is k 、h k Is defined as follows (M k And h k Only as an intermediate value for the calculation).
Wherein S is 1 For impedance matrix of each pipe section of heat supply pipe network, V 0 Constant term fitting coefficient matrix of flow lift quadratic polynomial of each relay pump of heat supply pipe network and V 1 First term fitting coefficient matrix of flow lift quadratic polynomial of each relay pump of heat supply pipe network and V 2 And a quadratic term fitting coefficient matrix of a flow lift quadratic polynomial of each relay pump of the heat supply pipe network. S is S 1 、V 0 、V 1 And V 2 The values of (2) are automatically resolved by the program according to the drawn CAD picture.
According to M k 、h k Calculating the flow change of the residual branch pipe section
(4) Calculating a simulation pressure matrix P of the heat supply pipe network node:
after determining the simulation pressure of the heat supply pipe network node, continuing to calculate a heat supply pipe network pipe section flow matrix G and leakage position impedance (i.e. virtual pipe section impedance) S 2 。
(5) And calculating a heat supply pipe network pipe section flow matrix G.
G 2 ＝Q 2
G＝[G 1 G 2 ]
(6) Updating leak impedance (i.e. virtual pipe section impedance) S 2 . Wherein B is a pressure matrix of the heat supply pressure of the leakage node relative to the main heat source, and B is actual measurement data.
In this embodiment, step S1021 obtains a block matrix from the basic association matrix and the independent loop matrix according to a preset partitioning rule, and divides the matrix into a branch part, a remaining branch part, a 0 matrix and a unit matrix, so that the mass balance equation and the loop energy balance equation can be conveniently used for calculation in the following steps. Step S1022 updates the simulation pressure of each pipe network node through an iteration method, and aims at enabling the simulation pressure to be consistent with the measured pressure, so that the accuracy of the leakage identification result is improved.
With continued reference to fig. 1, step S103 is performed after step S102.
Step S103, iteratively updating the flow parameters based on the gradient of the simulation pressure with respect to the flow parameters.
In one example, the flow parameters are iteratively updated. Fig. 5 is a flowchart illustrating iterative updating of flow parameters according to an exemplary embodiment, and as shown in fig. 5, step S103 specifically includes:
step S1031, determining a pressure error value according to the monitored pressure of the node and the simulation pressure of the node.
Calculating the error J between the monitored pressure and the simulated pressure new 。
J last ＝J new
Wherein J last The initial value of the pressure error is a preset value;for monitoring the pressure matrix at the heat exchange station, +.>Is the actual measurement data.
In step S1032, in response to the pressure error value being greater than the preset threshold, the flow parameter is iteratively updated based on the gradient of the simulated pressure with respect to the flow parameter.
If |J new -J las t| > epsilon (epsilon is an allowable error and is a preset value), and iteratively updating the flow parameters based on the gradient of the simulation pressure with respect to the flow parameters.
Calculating the gradient of node pressure with respect to leakage node flow
Gradient of leakage node flow based on node pressureUpdating leaky node traffic Q 2 ；
Wherein eta is the learning rate and is a preset value.
In this embodiment, step S1031 determines a pressure error value based on the monitored pressure of the node and the simulated pressure of the node, with the aim of matching the simulated pressure with the measured pressure. Step S1032 iteratively updates the leakage flow of each pipe network node based on the gradient descent algorithm, and based on the method, more accurate leakage flow data of the pipe network nodes can be obtained, thereby being beneficial to accurately judging whether the pipe network nodes belong to real leakage points or not.
With continued reference to fig. 1, step S104 is performed after step S103.
And step S104, taking the flow parameters acquired when the preset iteration conditions are met as leakage evaluation parameters to obtain a pipe network leakage identification model. And assuming that all nodes are leakage points, and determining whether the nodes belong to real leakage points according to the leakage flow of each node and a preset threshold value. And when the leakage node flow is greater than or equal to a preset threshold value, judging that the node belongs to a real leakage point.
In one example, a determination is made as to whether an iteration termination condition is met. In this example, step S104 specifically includes:
judging whether a preset iteration condition is met or not, including: the error value of the monitoring pressure and the simulation pressure is smaller than or equal to a preset threshold value, or the iteration number is equal to the preset iteration number.
If |J new -J last And stopping iteration if the I is less than or equal to epsilon or the iteration number is equal to the preset iteration number n. Outputting the node pressure matrix P, the pipe section flow matrix G and the virtual pipe section impedance matrix S calculated in the step S102 2 And the leaky node traffic matrix Q calculated in step S103 2 . Comparing the leakage node flow matrix with a preset threshold value, and when the leakage node flow is greater than or equal to the preset threshold value, determining that the node belongs to realityAnd (5) leaking points.
In the embodiment, an iteration method is adopted to iteratively calculate the error values of the simulation pressure and the actually measured pressure of the pipe network node, and finally, a more accurate leakage flow parameter is obtained, and based on the parameter, whether the node belongs to a real leakage point can be judged more accurately. And meanwhile, the simulation pressure, the pipe section flow and the virtual pipe section impedance of the node are output, more comprehensive leakage point information is obtained, and decisions and supports are provided for relevant heat supply pipe network maintenance personnel to maintain the heat supply pipe network.
In summary, according to the method for constructing the pipe network leakage identification model provided by the embodiment of the application, the basic association matrix and the independent loop matrix are constructed, and the equation set formed by the heat supply pipe network node mass balance equation, the pipe section pressure drop equation and the loop energy balance equation is linearized, so that the leakage flow of each branch pipe section is iteratively updated based on the gradient descent algorithm. The simulation pressure at each heat exchange station is consistent with the measured pressure, leakage flow rate at the end of iteration is used as leakage evaluation parameters, and parameters such as an impedance matrix of leakage nodes are output, so that the problem of leakage identification of the tail end branch pipes of the heat supply pipe network is solved, the leakage identification result is more accurate and comprehensive, and decisions and supports are provided for relevant heat supply pipe network management staff to maintain the heat supply pipe network.
Second aspect
An embodiment of the application provides a pipe network leakage identification method, and fig. 6 is a block diagram of the pipe network leakage identification method according to the embodiment of the application. As shown in fig. 6, the method includes:
step S105, acquiring a current pipe network topological structure and monitoring pressure of each node in the pipe network, and constructing a complement pipe network topological structure based on the current pipe network topological structure;
step S106, determining leakage evaluation parameters of each node based on the complement network management topological structure and the monitoring pressure through the pipe network leakage identification model constructed by the embodiment of the first aspect;
step S107, identifying the leakage points according to the leaked evaluation parameters and preset evaluation criteria.
In one example, step S1071 includes: and determining a first maintenance strategy in response to the leakage evaluation parameter of the node being located in the first interval. And determining a second maintenance strategy in response to the leakage evaluation parameter of the node being located in the second interval, wherein the priorities of the first maintenance strategy and the second maintenance strategy are different.
By adopting the mode, based on different leakage evaluation parameter results, the leakage degree corresponding to different evaluation intervals is different, and corresponding maintenance strategies are determined according to different evaluation intervals, so that data guidance and support are provided for heat supply pipe network maintenance personnel.
In summary, the pipe network leakage identification method provided by the embodiment of the present application determines the leakage evaluation parameter of each node based on the completed pipe network topology structure and the monitoring pressure and through the pipe network leakage identification model constructed by the embodiment of the first aspect. Identifying leakage points according to the leakage evaluation parameters and preset evaluation criteria, determining the interval where the leakage points are located, and determining the corresponding maintenance strategy. The pipe network leakage identification method provided by the embodiment of the application enables the leakage identification result to be more accurate and comprehensive, divides the leakage interval, corresponds to different leakage degrees and maintenance strategies, and provides decisions and supports for relevant heat supply pipe network management personnel to maintain the heat supply pipe network.
Third aspect of the application
The embodiment of the application provides a pipe network leakage identification model construction device, the system is used for realizing the embodiment and implementation mode, and fig. 7 is a block diagram of the pipe network leakage identification model construction device according to the embodiment of the application. As shown in fig. 7, the pipe network leakage identification system includes: a matrix determination module 100, a simulation pressure determination module 200, a flow parameter update module 300, an identification module 400.
The matrix determining module 100 is configured to determine a basic association matrix and an independent loop matrix based on an actual pipe network topology structure, where the basic association matrix characterizes a layout relationship of pipe segments and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and a topological connection relationship between the loops;
the simulation pressure determining module 200 is configured to determine a simulation pressure of the node through a simulation algorithm based on the basic association matrix, the independent loop matrix and the flow parameters of the node;
a flow parameter updating module 300, configured to iteratively update the flow parameter based on a gradient of the simulation pressure with respect to the flow parameter;
the identification module 400 is configured to obtain a pipe network leakage identification result by using the flow parameter obtained when the preset iteration condition is satisfied as a leakage evaluation parameter.
In one example, the matrix determination module 100 includes:
the construction unit: constructing a virtual node on each end branch pipe in an actual pipe network topological structure, and constructing a virtual pipe section taking the virtual node as a starting point to obtain a completed pipe network topological structure;
a determination unit: and determining a basic association matrix and an independent loop matrix according to the layout relation of the pipe sections and the nodes in the full pipe network topological structure.
In one example, the simulated pressure determination module 200 includes:
a blocking unit: obtaining a blocking matrix according to a preset blocking rule by the basic association matrix and the independent loop matrix;
a determination unit: and determining the simulation pressure of the node based on the block matrix and the flow parameters of the node.
In one example, the flow parameter update module 300 includes:
a determination unit: determining a pressure error value according to the monitoring pressure of the node and the simulation pressure of the node;
an updating unit: and in response to the pressure error value being greater than a preset threshold, iteratively updating the flow parameter based on a gradient of the simulated pressure with respect to the flow parameter.
In one example, the recognition module 400 includes: and judging that the error value of the monitored pressure and the simulated pressure is smaller than or equal to a preset threshold value, or the iteration number is equal to the preset iteration number.
In summary, the device for constructing the pipe network leakage identification model according to the embodiment of the present application constructs the basic association matrix and the independent loop matrix based on the pipe network topology structure, and blocks the matrix by using the preset rule. Based on the obtained blocking matrix and node leakage flow parameters, calculating simulation pressure through a heating pipe network hydraulic simulation algorithm. And iteratively updating the leakage flow of each branch pipe section according to the gradient descent algorithm. The simulation pressure at each heat exchange station is consistent with the measured pressure, leakage flow rate at the end of iteration is used as leakage evaluation parameters, and parameters such as an impedance matrix of leakage nodes are output, so that the problem of leakage identification of the tail end branch pipes of the heat supply pipe network is solved, the leakage identification result is more accurate and comprehensive, and decisions and supports are provided for relevant heat supply pipe network management staff to maintain the heat supply pipe network.
Fourth aspect of
The embodiment of the application provides a pipe network leakage identification device. FIG. 8 is a block diagram of a pipe network leak identification device according to an embodiment of the application. As shown in fig. 8, the pipe network leakage identification system includes: a construction module 500, a determination module 600, an identification module 700.
The construction module 500 is configured to obtain a current pipe network topology structure and monitoring pressures of nodes in a pipe network, and construct a complement pipe network topology structure based on the current pipe network topology structure;
a determining module 600, configured to determine, based on the completed network management topology structure and the monitoring pressure, a leakage evaluation parameter of each node by using the pipe network leakage identification model constructed by the method of any one of claims 1 to 5;
the identifying module 700 is configured to identify the leakage point according to the leakage evaluation parameter and a preset evaluation standard.
In one embodiment, the identifying module 700 is further configured to determine a first maintenance policy in response to the leakage assessment parameter of the node being located in a first interval; and determining a second maintenance strategy in response to the leakage evaluation parameter of the node being located in the second interval, wherein the priorities of the first maintenance strategy and the second maintenance strategy are different.
In summary, the pipe network leakage identification device provided by the embodiment of the present application determines the leakage evaluation parameter of each node based on the completed pipe network topology structure and the monitored pressure and through the pipe network leakage identification model constructed by the embodiment of the first aspect. Identifying leakage points according to the leakage evaluation parameters and preset evaluation criteria, determining the interval where the leakage points are located, and determining the corresponding maintenance strategy. In this way, the leakage identification result is more accurate and comprehensive, the leakage interval is divided, different leakage degrees and maintenance strategies are corresponding, and decisions and supports are provided for relevant heat supply pipe network management personnel to maintain the heat supply pipe network.
Fifth aspect of
An embodiment of the present application provides a computer readable storage medium, on which a computer program is stored, where the program when executed by a processor implements the steps of the method for constructing a pipe network leakage identification model provided in the first aspect and the method for identifying pipe network leakage provided in the second aspect.
More specifically, among others, readable storage media may be employed including, but not limited to: portable disk, hard disk, random access memory, read only memory, erasable programmable read only memory, optical storage device, magnetic storage device, or any suitable combination of the foregoing.
In a possible embodiment, the application may also be implemented in the form of a program product comprising program code for causing a terminal device to carry out the steps of the construction method implementing the pipe network leakage identification model of example 1, when said program product is run on the terminal device.
Wherein the program code for carrying out the application may be written in any combination of one or more programming languages, which program code may execute entirely on the user device, partly on the user device, as a stand-alone software package, partly on the user device and partly on the remote device or entirely on the remote device.
In a possible implementation manner, the present application may also be implemented in the form of a program product, which includes program code for causing a terminal device to perform the steps of implementing the method for constructing a pipe network leakage identification model provided in the first aspect and the method for identifying pipe network leakage provided in the second aspect, when the program product is run on the terminal device.
Wherein the program code for carrying out the application may be written in any combination of one or more programming languages, which program code may execute entirely on the user device, partly on the user device, as a stand-alone software package, partly on the user device and partly on the remote device or entirely on the remote device.
The technical features of the above-described embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above-described embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples illustrate only a few embodiments of the application, which are described in detail and are not to be construed as limiting the scope of the application. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the application, which are all within the scope of the application. Accordingly, the scope of protection of the present application is to be determined by the appended claims.
Claims (10)
1. The construction method of the pipe network leakage identification model is characterized by comprising the following steps of:
determining a basic incidence matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic incidence matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and topological connection relations among the loops;
determining the simulation pressure of the node through a simulation algorithm based on the basic association matrix, the independent loop matrix and the flow parameters of the node;
iteratively updating the flow parameter based on the gradient of the simulated pressure with respect to the flow parameter;
and taking the flow parameters acquired when the preset iteration conditions are met as leakage evaluation parameters to obtain a pipe network leakage identification model.
2. The method for constructing a pipe network leakage identification model according to claim 1, wherein the determining a basic association matrix and an independent loop matrix based on an actual pipe network topology structure comprises:
constructing a virtual node on each tail end branch pipe in the actual pipe network topological structure, and constructing a virtual pipe section taking the virtual node as a starting point to obtain a complement pipe network topological structure;
and determining the basic association matrix and the independent loop matrix according to the layout relation of the pipe sections and the nodes in the full pipe network topological structure.
3. The method for constructing a pipe network leakage identification model according to claim 2, wherein determining the simulation pressure of the node by a simulation algorithm based on the basic correlation matrix, the independent loop matrix and the flow parameters of the node comprises:
obtaining a blocking matrix according to a preset blocking rule by the basic association matrix and the independent loop matrix;
and determining the simulation pressure of the node based on the blocking matrix and the flow parameters of the node.
4. The method for constructing a pipe network leakage identification model according to claim 3, wherein iteratively updating the flow parameters based on the gradient of the simulation pressure with respect to the flow parameters comprises:
determining a pressure error value according to the monitoring pressure of the node and the simulation pressure of the node;
and iteratively updating the flow parameter based on a gradient of the simulated pressure with respect to the flow parameter in response to the pressure error value being greater than a preset threshold.
5. The method for constructing a pipe network leakage identification model according to claim 4, wherein the step of obtaining the pipe network leakage identification model by using the flow parameter obtained when the preset iteration condition is satisfied as a leakage evaluation parameter comprises the steps of:
the meeting the preset iteration condition comprises the following steps: the error value of the monitoring pressure and the simulation pressure is smaller than or equal to the preset threshold value, or the iteration number is equal to the preset iteration number.
6. A method for identifying pipe network leakage, comprising:
acquiring a current pipe network topological structure and monitoring pressure of each node in a pipe network, and constructing a complement pipe network topological structure based on the current pipe network topological structure;
determining leakage evaluation parameters of each node based on the complement pipe network topology structure and the monitoring pressure by using the pipe network leakage identification model constructed by the method of any one of claims 1 to 5;
and identifying the leakage point according to the leakage evaluation parameter and a preset evaluation standard.
7. The pipe network leakage identification method according to claim 6, wherein the preset evaluation criteria comprises a first interval and a second interval; the method further comprises the steps of:
determining a first maintenance strategy in response to the leakage evaluation parameter of the node being located in the first interval;
and determining a second maintenance strategy in response to the leakage evaluation parameter of the node being located in the second interval, wherein the priorities of the first maintenance strategy and the second maintenance strategy are different.
8. The utility model provides a pipe network seepage discernment model construction device which characterized in that includes:
the matrix determining module is used for determining a basic incidence matrix and an independent loop matrix based on an actual pipe network topological structure, wherein the basic incidence matrix represents the layout relation of pipe sections and nodes in a pipe network, and the independent loop matrix reflects independent loops in a heat supply pipe network and topological connection relation among the loops;
the simulation pressure determining module is used for determining the simulation pressure of the node through a simulation algorithm based on the basic incidence matrix, the independent loop matrix and the flow parameters of the node;
a flow parameter updating module for iteratively updating the flow parameter based on the gradient of the simulation pressure with respect to the flow parameter;
and the identification module is used for taking the flow parameter acquired when the preset iteration condition is met as a leakage evaluation parameter to acquire a pipe network leakage identification result.
9. A pipe network leakage identification device is characterized by comprising
The construction module is used for acquiring the current pipe network topological structure and the monitoring pressure of each node in the pipe network, and constructing a complement pipe network topological structure based on the current pipe network topological structure;
the determining module is used for determining leakage evaluation parameters of each node based on the complement network management topological structure and the monitoring pressure through the pipe network leakage identification model constructed by the method of any one of claims 1-5;
and the identification module is used for identifying the leakage points according to the leakage evaluation parameters and preset evaluation standards.
10. A computer-readable storage medium having stored thereon a computer program, wherein the program when executed by a processor implements the pipe network leak identification model construction method according to any one of claims 1 to 5 and the pipe network leak identification method according to any one of claims 6 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310882367.XA CN117113862A (en) | 2023-07-18 | 2023-07-18 | Pipe network leakage identification model construction method and pipe network leakage identification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310882367.XA CN117113862A (en) | 2023-07-18 | 2023-07-18 | Pipe network leakage identification model construction method and pipe network leakage identification method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117113862A true CN117113862A (en) | 2023-11-24 |
Family
ID=88802888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310882367.XA Pending CN117113862A (en) | 2023-07-18 | 2023-07-18 | Pipe network leakage identification model construction method and pipe network leakage identification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117113862A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117408007A (en) * | 2023-12-15 | 2024-01-16 | 华东交通大学 | Three-dimensional pipe network initial flow distribution method, system, storage medium and electronic equipment |
-
2023
- 2023-07-18 CN CN202310882367.XA patent/CN117113862A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117408007A (en) * | 2023-12-15 | 2024-01-16 | 华东交通大学 | Three-dimensional pipe network initial flow distribution method, system, storage medium and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117113862A (en) | Pipe network leakage identification model construction method and pipe network leakage identification method | |
| Cavazzini et al. | Optimal assets management of a water distribution network for leakage minimization based on an innovative index | |
| Marzola et al. | Leakage detection and localization in a water distribution network through comparison of observed and simulated pressure data | |
| CN116108604B (en) | Water supply network abnormality detection method, system, equipment and storage medium | |
| CN114757589B (en) | Data processing method, server and storage medium | |
| CN117114369B (en) | Intelligent gas emergency gas supply disposal method, internet of things system, device and medium | |
| CN117235950B (en) | Natural gas pipe network steady-state simulation method, medium and equipment based on Newton iteration method | |
| CN105654151B (en) | A kind of workpiece localization method and positioning system | |
| CN114418025A (en) | Real-time interaction method, device and equipment based on BIM data and Internet of things data | |
| CN112686971A (en) | Method and system for undirected graph orientation of building system relation topology | |
| CN116257559A (en) | Method and device for associating person with address based on space grid and electronic equipment | |
| CN115100378A (en) | Method and device for quickly constructing building group | |
| CN114676490A (en) | Mine intelligent ventilation network resolving method | |
| CN108304511A (en) | A kind of river network storage expression based on XML data format | |
| WO2024207867A1 (en) | Water supply pipeline topology verification method and system, storage medium, and intelligent terminal | |
| CN106575312B (en) | Method for detecting anomalies in a distribution network, in particular a water distribution network | |
| CN111553041B (en) | Gas pipe network parameter calibration method, system, terminal and storage medium | |
| CN112906121A (en) | Heating demand design method and system for floor heating system | |
| CN115982918B (en) | Resistance characteristic identification method, device, equipment and storage medium | |
| CN108133116A (en) | CAD floor plan wall thickness automatic identifying methods and system | |
| CN112528384A (en) | BIM-based electromechanical professional component comprehensive arrangement method, system and medium | |
| CN108599976B (en) | Network quality evaluation model construction method and device and computer equipment | |
| CN115809531A (en) | Simulation method and system for natural gas pipe network, storage medium and electronic equipment | |
| CN108733725B (en) | Method and system for establishing network semantic model of target system | |
| CN112612784B (en) | River basin calculation unit automatic dividing method and device and computer equipment |
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 |