Water supply network pipe burst positioning method, device and equipment and storage medium

Abstract

The invention relates to the technical field of municipal engineering, in particular to a water supply network pipe burst positioning method, device and equipment and a storage medium. According to the water supply network pipe burst positioning method provided by the embodiment of the invention, the moment when the pipe burst pressure wave reaches each monitoring point is determined through the high-frequency pressure monitoring data, the possible pipe burst position is delineated on the basis of the moment and the propagation time of the transient flow pressure wave, and then the transient flow pressure data at each monitoring point are simulated and calculated on the basis of the transient flow model. And through considering the continuous change information of the pressure wave, the simulated transient pressure data is matched with the high-frequency monitoring data, so that the pipe explosion position can be accurately positioned from the possible pipe explosion positions. Therefore, the problem that the positioning result is easily influenced by the wave velocity estimation error due to the fact that the pipe explosion position is positioned only based on the time difference of the transient pressure wave reaching each monitoring point in the prior art is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of municipal engineering, and in particular to a method, device, equipment and storage medium for locating burst pipes in a water supply network. Background Art

[0002] In the water supply system, the water supply network is responsible for delivering treated water to thousands of households, and is an important part of ensuring water supply safety. Pipe bursts will affect the water supply capacity and water quality of the water supply network. Therefore, it is particularly important to detect and locate pipe bursts in the water supply network as early as possible.

[0003] With the development of high-frequency pressure monitoring equipment, it has become possible to perform high-frequency pressure monitoring in water supply networks at a relatively low cost. Compared with traditional monitoring equipment, high-frequency pressure monitoring equipment can capture instantaneous pressure changes in the network, which is conducive to more timely and accurate positioning of pipe burst events in the network. However, the current method for locating pipe bursts in water supply networks based on high-frequency pressure monitoring data is simply based on the time difference between the transient pressure waves caused by the pipe burst reaching each monitoring point, which makes the positioning results easily affected by wave velocity estimation errors. Summary of the invention

[0004] In view of this, the present invention provides a method, device, equipment and storage medium for locating a burst pipe in a water supply network, so as to solve the problem that the burst pipe locating method in the prior art is easily affected by wave velocity estimation errors.

[0005] In a first aspect, the present invention provides a method for locating a pipe burst in a water supply network, the method comprising: obtaining high-frequency pressure monitoring data of each monitoring point in the water supply network, the high-frequency pressure monitoring data comprising a pressure value and a monitoring time corresponding to each pressure value; based on the high-frequency pressure monitoring data, using a cumulative sum algorithm to determine each time when the pipe burst pressure wave arrives at each monitoring point; determining a possible pipe burst location based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network; matching new steady-state pressure monitoring data obtained after the pipe burst with pressure estimation data corresponding to different pipe burst coefficients determined based on a steady-state model to determine the pipe burst coefficient of the possible pipe burst location; based on a transient flow model, using the pipe burst coefficient to simulate and calculate the transient flow pressure data of each possible pipe burst location, matching the transient flow pressure data with the high-frequency pressure monitoring data to determine the pipe burst location.

[0006] The method for locating a pipe burst in a water supply network provided by an embodiment of the present invention determines the time when the pipe burst pressure wave reaches each monitoring point through high-frequency pressure monitoring data, first delineates the possible pipe burst position based on the time and the propagation time of the transient flow pressure wave, then simulates and calculates the transient flow pressure data at each monitoring point based on the transient flow model, and then matches the simulated transient pressure data with the high-frequency monitoring data by considering the information of the continuous change of the pressure wave, thereby realizing accurate positioning of the pipe burst position from the possible pipe burst positions. Thus, the problem that the positioning result of the pipe burst position is easily affected by the wave velocity estimation error caused by the positioning of the pipe burst position based only on the time difference when the transient pressure wave reaches each monitoring point in the related art is solved.

[0007] In an optional embodiment, based on the high-frequency pressure monitoring data, a cumulative sum algorithm is used to determine the time when the burst pressure wave arrives at each monitoring point, including: using the cumulative sum algorithm to calculate the time when the cumulative pressure change is greater than a threshold based on the high-frequency pressure monitoring data, and determining that a burst pressure wave occurs when the cumulative pressure change is greater than the threshold; determining the time when the burst pressure wave arrives at each monitoring point according to the time when the pressure change recorded in the cumulative sum algorithm starts to accumulate.

[0008] In this embodiment, a cumulative sum algorithm is used to determine the time when the pipe burst pressure wave appears, which provides a data basis for determining the possible pipe burst position.

[0009] In an optional embodiment, before determining the possible burst location based on the relationship between the transient flow pressure wave propagation time and each time instant between any two nodes in the water supply network, the method further includes: determining the nodes of the water supply network and the pipes between the nodes based on the topological structure of the water supply network; and calculating the transient flow pressure wave propagation time in any pipe based on the length of the pipe between the two nodes and the transient flow wave velocity of the pipe.

[0010] In an optional embodiment, the possible burst location is determined based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network, including: using the shortest path algorithm to determine the shortest time for the pressure wave to propagate between any two nodes based on the transient flow pressure wave propagation time in any pipe of the water supply network; determining the objective function based on the minimum value of the difference between the time when the burst pressure wave reaches any two monitoring points and the difference between the shortest time for the pressure wave to propagate from any node to the corresponding two monitoring points, and solving the objective function to determine the possible burst location.

[0011] In an optional implementation, the objective function is expressed using the following formula:

[0012]

[0013] In the formula, k represents any node of the water supply network, M represents the set of monitoring points, and t i represents the time when the burst pressure wave reaches the monitoring point i, tj represents the time when the burst pressure wave reaches the monitoring point j, T ik represents the shortest time for transient flow pressure wave propagation between monitoring point i and node k, T jk It represents the shortest time for transient flow pressure wave to propagate between monitoring point j and node k.

[0014] In this embodiment, based on the transient flow pressure wave that occurs when a pipe burst occurs, the possible pipe burst location is first determined according to the propagation time of the transient flow pressure wave and the time when the pipe burst pressure wave reaches each monitoring point, thereby achieving preliminary screening of the pipe burst location.

[0015] In an optional implementation, new steady-state pressure monitoring data acquired after a pipe burst is matched with pressure estimation data corresponding to different pipe burst coefficients determined based on a steady-state model to determine the pipe burst coefficient of a possible pipe burst location, including: acquiring new steady-state pressure monitoring data of each monitoring point after a pipe burst occurs; calculating pressure estimation data of a possible pipe burst location under different pipe burst coefficients based on a pre-constructed steady-state model of a water supply network; determining an objective function based on the difference between the new steady-state pressure monitoring data and the pressure estimation data, solving the objective function using a genetic algorithm, and determining the pipe burst coefficient of each possible pipe burst location.

[0016] In this embodiment, a steady-state model is used to calculate the pressure estimation data of possible pipe burst positions under different pipe burst coefficients, and the pipe burst coefficient of each possible pipe burst position is determined based on the relationship between the pressure estimation data and the pressure monitoring data, thereby providing a data basis for the subsequent analysis of the transient flow model.

[0017] In an optional implementation, based on a transient flow model and using a pipe burst coefficient, transient flow pressure data of a pipe burst at each possible pipe burst position is simulated and calculated, and the transient flow pressure data is matched with high-frequency pressure monitoring data to determine the pipe burst position, including: based on a transient flow model and using a pipe burst coefficient, transient flow pressure data of a pipe burst at each possible pipe burst position is simulated and calculated; a dynamic time warping algorithm is used to match the transient flow pressure data of each monitoring point with the high-frequency pressure monitoring data, and the pipe burst position is determined according to the matching result.

[0018] In this embodiment, a transient flow model is used to calculate the pressure data of each monitoring point corresponding to the burst coefficient of the possible burst position, thereby realizing the utilization of the continuous transformation information of the pressure wave; when the dynamic time warping algorithm is used for matching, it can be achieved that the points in one sequence correspond to other time points in another sequence, and the influence of shape stretching or compression can be eliminated.

[0019] In a second aspect, the present invention provides a device for locating a pipe burst in a water supply network, the device comprising: a data acquisition module, for acquiring high-frequency pressure monitoring data of each monitoring point in the water supply network, the high-frequency pressure monitoring data comprising a pressure value and a monitoring time corresponding to each pressure value; a pipe burst monitoring module, for determining, based on the pressure monitoring data, each time when the pipe burst pressure wave reaches each monitoring point by using a cumulative sum algorithm; a preliminary positioning module, for determining a possible pipe burst location based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network; a pipe burst coefficient determination module, for matching new steady-state pressure monitoring data acquired after the pipe burst with pressure estimation data corresponding to different pipe burst coefficients determined based on a steady-state model, to determine the pipe burst coefficient of the possible pipe burst location; a positioning module, for simulating and calculating the transient flow pressure data of each possible pipe burst location by using the pipe burst coefficient based on a transient flow model, matching the transient flow pressure data with the high-frequency pressure monitoring data, to determine the pipe burst location.

[0020] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for locating burst pipes in a water supply network according to the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0021] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for locating burst pipes in a water supply network according to the first aspect or any corresponding embodiment thereof.

[0022] In a fifth aspect, the present invention provides a computer program product, including computer instructions, and the computer instructions are used to enable a computer to execute the method for locating burst pipes in a water supply network according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0024] Figure 1 is a schematic flow chart of a method for locating a burst pipe in a water supply network according to an embodiment of the present invention;

[0025] Figure 2 is a schematic diagram of the location of pipe burst in a water supply network according to an embodiment of the present invention;

[0026] Figure 3is a schematic diagram of locating a burst pipe in a water supply network according to an embodiment of the present invention;

[0027] Figure 4 is a structural block diagram of a device for locating a burst pipe in a water supply network according to an embodiment of the present invention;

[0028] Figure 5 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0029] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0030] According to an embodiment of the present invention, an embodiment of a method for locating a burst pipe in a water supply network is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0031] In this embodiment, a method for locating a burst pipe in a water supply network is provided, which can be used in electronic devices such as computers, mobile phones, tablet computers, etc. Figure 1 FIG. 4 is a flow chart of a method for locating a burst pipe in a water supply network according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0032] Step S101, obtain the pressure monitoring data of each monitoring point in the water supply network, the pressure monitoring data includes the pressure value and the monitoring time corresponding to the pressure value. Specifically, the number and location of the monitoring points can be determined according to the topological structure of the water supply network. At the same time, the actual situation such as cost can also be considered when setting the monitoring points, so as to determine the final monitoring points. After determining the monitoring points, high-frequency pressure sensors can be set at the monitoring points to monitor the pressure and obtain high-frequency pressure monitoring data. Among them, when performing pressure monitoring, high-frequency pressure monitoring data can be collected based on a pre-set sampling frequency, and the high-frequency pressure monitoring data obtained thereby includes the pressure values ​​specifically collected by the high-frequency pressure sensor and the collection time corresponding to each pressure value.

[0033] Step S102, based on the high-frequency pressure monitoring data, a cumulative sum algorithm is used to determine the time when the burst pressure wave reaches each monitoring point. Specifically, when a burst occurs in a water supply pipeline, a pressure wave will be generated and propagated in the pipe network, so whether a burst occurs and the time of the burst can be determined by processing and calculating the high-frequency pressure monitoring data. In this embodiment, a cumulative sum algorithm is used to process the pressure monitoring data. It quickly identifies when the system deviates from the normal operating state, that is, a burst occurs, by calculating the cumulative deviation of the pressure monitoring data stream, and determines the time when the burst pressure wave reaches each monitoring point according to the time when the deviation begins to accumulate.

[0034] Step S103, determining the possible burst location based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network. The nodes in the water supply network can be determined based on the topological structure of the water supply network and the graph theory algorithm. For example, the water supply network can be abstracted into a weighted graph using a graph theory algorithm. The vertices in the weighted graph represent the nodes in the network. The edges in the weighted graph represent the pipes in the network, connecting two nodes, and are undirected edges, and their weights are the propagation time of the transient flow pressure wave in the pipe corresponding to the edge.

[0035] Specifically, when a pipe burst occurs, a transient flow pressure wave appears in the pipeline, and the propagation time of the transient flow pressure wave between any two nodes can be predetermined based on the wave velocity of the transient flow pressure wave. Therefore, the location where the pipe burst may occur can be delineated based on the propagation time and the time when the burst pressure wave reaches each monitoring point. In addition, it should be noted that when determining the monitoring points, they should be determined from the nodes of the water supply network, that is, the nodes of the water supply network should be determined first, and then some nodes should be selected as monitoring points.

[0036] Step S104, matching the new steady-state pressure monitoring data obtained after the pipe burst with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model, and determining the pipe burst coefficient of the possible pipe burst position. Specifically, when a pipe burst occurs, the pressure fluctuation will dissipate within tens of seconds and reach a new steady state. Therefore, after the pipe burst is detected, the new steady-state pressure monitoring data of each monitoring point can be obtained. At the same time, for each possible pipe burst position, the steady-state model and the pipe burst coefficient can also be used to calculate the pressure estimation data of each monitoring point; then, the pipe burst coefficient corresponding to the possible pipe burst position is determined based on the relationship between the new steady-state pressure monitoring data and the pressure estimation data.

[0037] Step S105, based on the transient flow model, using the burst coefficient, simulate and calculate the transient flow pressure data of each possible burst position, match the transient flow pressure data with the high-frequency pressure monitoring data, and determine the burst position. Specifically, after determining the burst coefficient of the possible burst position, it can be substituted into the transient flow model to calculate the transient flow pressure data of the corresponding monitoring point when the burst occurs at the possible burst position. Match the calculated transient flow pressure data with the high-frequency pressure monitoring data, and accurately locate the burst position from the possible burst positions according to the matching result.

[0038] The method for locating a pipe burst in a water supply network provided by an embodiment of the present invention determines the time when the pipe burst pressure wave reaches each monitoring point through high-frequency pressure monitoring data, first delineates the possible pipe burst position based on the time and the propagation time of the transient flow pressure wave, then simulates and calculates the transient flow pressure data at each monitoring point based on the transient flow model, and then matches the simulated transient pressure data with the high-frequency monitoring data by considering the information of the continuous change of the pressure wave, thereby realizing accurate positioning of the pipe burst position from the possible pipe burst positions. Thus, the problem that the positioning result of the pipe burst position is easily affected by the wave velocity estimation error caused by the positioning of the pipe burst position based only on the time difference when the transient pressure wave reaches each monitoring point in the related art is solved.

[0039] In this embodiment, a method for locating a burst pipe in a water supply network is provided, and the process includes the following steps:

[0040] Step S201, obtaining high-frequency pressure monitoring data of each monitoring point in the water supply network, wherein the high-frequency pressure monitoring data includes pressure values ​​and monitoring times corresponding to each pressure value. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0041] Step S202: Based on the pressure monitoring data, a cumulative sum algorithm is used to determine the time when the burst pressure wave reaches each monitoring point.

[0042] Specifically, the above step S202 includes:

[0043] Step S2021, using the cumulative sum algorithm, based on the high-frequency pressure monitoring data, calculates the time when the pressure change is accumulated to be greater than the threshold value, and when the pressure change is accumulated to be greater than the threshold value, it is determined that a burst pressure wave occurs. Specifically, for any monitoring point, the monitored pressure change is H(t). In this embodiment, the high-frequency pressure monitoring data is collected at a certain sampling frequency, and thus, the collected high-frequency pressure monitoring data is a discrete data sequence, which can be represented by:

[0044] H n =H(nΔt)

[0045] Where Δt represents the time step, which is the inverse of the sampling frequency, and n represents the nth time step.

[0046] Then the pressure change at the nth time step is:

[0047] ΔH n =H n -H n-1

[0048] Using the Cumulative Sum (CUSUM) algorithm, A is iteratively calculated as follows: n , B n and S n :

[0049] A0=B0=0

[0050] A n =max(A n-1 -ΔH n -ν,0)

[0051]

[0052] Where ν is the drift parameter of the CUSUM algorithm. n is the cumulative sum of the CUSUM (cumulative sum) algorithm; S n is the sequence number of the time when the accumulation starts. The accumulation starts, reflecting that the pressure of the system begins to change. n If the pressure drop is greater than the threshold value K, a pressure drop warning is issued, indicating that the accumulated pressure change is greater than the threshold value, and it is believed that a pipe burst pressure wave may have occurred. At the same time, S is recorded. n , which is the ordinal number of the moment when the pressure change begins to accumulate.

[0053] Step S2022, according to the time when the pressure change recorded in the cumulative sum algorithm starts to accumulate, determine the time when the burst pressure wave arrives at the monitoring point. Specifically, the time when the burst pressure wave arrives at the monitoring point is S n Δt. For each monitoring point, the S determined in step S2021 can be used respectively. n , and further according to S n Δt is the time when the burst pressure wave reaches the monitoring point.

[0054] Step S203, based on the topological structure of the water supply network, determine the nodes of the water supply network and the pipes between the nodes; calculate the transient flow pressure wave propagation time in any pipe based on the length of the pipe between the two nodes and the transient flow wave velocity of the pipe. Specifically, through the graph theory algorithm, it is determined that the set of all nodes of the water supply network can be recorded as V, the set of all pipes can be recorded as E, and the pipe network can be recorded as G(V,E). The adjacency matrix of graph G can be recorded as A, and its elements are:

[0055]

[0056] Among them, τ ij represents the pipeline weight, which can be calculated as follows:

[0057]

[0058] Where, L ij is the length of the pipeline between node i and node j, a ij is the transient flow wave velocity of the pipeline between node i and node j. Pipeline weight τ ij The physical meaning of is the time it takes for the transient flow pressure wave to pass through the pipe between node i and node j.

[0059] Step S204, determining a possible pipe burst location based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network.

[0060] Specifically, the above step S204 includes:

[0061] Step S2041, using the shortest path algorithm, based on the transient flow pressure wave propagation time in any pipe of the water supply network, determine the shortest time for the pressure wave to propagate between any two nodes; specifically, using the shortest path algorithm such as the Dijkstra algorithm on the above matrix A, determine the shortest time T required for the transient flow pressure wave to propagate between any point i and another point j in the water supply network. ij Among them, the method of using Dijkstra algorithm to determine the shortest time is similar to the shortest path. When the weight between nodes is time instead of distance, one node can be used as the starting point and the other node as the end point to obtain the shortest time required from the starting point to the end point.

[0062] Step S2042, determine the objective function based on the minimum value of the difference between the time when the burst pressure wave reaches any two monitoring points and the difference between the shortest propagation time of the pressure wave from any node to the corresponding two monitoring points, and solve the objective function to determine the possible burst position. In this embodiment, the set of monitoring points is recorded as M. If the burst occurs at the kth node in the pipe network, then theoretically:

[0063]

[0064] Among them, t i and t j is the time when the pressure fluctuation caused by the pipe burst propagates to monitoring points i and j. This formula indicates that the difference between the propagation time difference of the pressure wave caused by the pipe burst and any two monitoring points and the propagation time difference of the transient flow pressure wave is theoretically 0. However, in actual situations, due to factors such as monitoring errors, even if the pipe burst occurs at k, S kIt is not 0. Therefore, minimization can be used as the objective function, and the possible burst position can be determined by solving the objective function.

[0065] Specifically, the objective function is expressed as follows:

[0066]

[0067] In the formula, k represents any node of the water supply network, M represents the set of monitoring points, and t i represents the time when the burst pressure wave reaches the monitoring point i, t j represents the time when the burst pressure wave reaches the monitoring point j, T ik represents the shortest time for transient flow pressure wave propagation between monitoring point i and node k, T jk It represents the shortest time for transient flow pressure wave to propagate between monitoring point j and node k.

[0068] When solving the objective function, an exhaustive method can be used to find the optimal solution and several suboptimal solutions as possible pipe bursting positions. The set of these possible pipe bursting positions is denoted as Y.

[0069] Step S205 , matching the new steady-state pressure monitoring data acquired after the pipe burst with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model to determine the pipe burst coefficient of the possible pipe burst position.

[0070] Specifically, the above step S205 includes:

[0071] Step S2051, obtain the new steady-state pressure monitoring data of each monitoring point after the pipe burst occurs; specifically, when the pipe burst occurs, the pressure fluctuation will dissipate within tens of seconds and reach a new steady state. Therefore, after the pipe burst is detected, the high-frequency pressure sensor can be used to re-acquire the new steady-state pressure monitoring data of each monitoring point. In this embodiment, the steady-state water head monitoring value obtained is H i,obs As new steady-state pressure monitoring data.

[0072] Step S2052, based on the pre-built steady-state model of the water supply network, calculate the pressure estimation data of the possible burst position under different burst coefficients; specifically, the steady-state model of the water supply network can be constructed using EPANET software, such as first constructing a model of the water supply network, and then setting boundary conditions to configure the model, so as to obtain a steady-state model. After obtaining the steady-state model, the steady-state model can be used to perform steady-state analysis, such as calculating the burst coefficient at the possible burst position k, such as the orifice outflow coefficient size c k (c k >0), the new steady-state water head estimate H of each monitoring point i i,cal (k,c k) as the pressure estimation data, where the burst coefficient c k The size of the tube also reflects the size of the burst tube.

[0073] Step S2053, determine the objective function based on the difference between the new steady-state pressure monitoring data and the pressure estimation data, use a genetic algorithm to solve the objective function, and determine the burst coefficient of each possible burst position. Specifically, for each possible burst position, the pressure estimation data corresponding to different burst coefficients and the new steady-state pressure monitoring data can be used for difference calculation. When the difference is the smallest, the corresponding burst coefficient is the burst coefficient corresponding to the possible burst position. Therefore, for each possible burst position, the following objective function can be used to determine the corresponding burst coefficient:

[0074]

[0075] The objective function can be solved by a genetic algorithm. The specific solution method can be implemented by referring to related technologies and will not be described in detail here.

[0076] Step S206, based on the transient flow model and using the pipe burst coefficient, simulate and calculate the transient flow pressure data of pipe burst at each possible pipe burst position, match the transient flow pressure data with the high-frequency pressure monitoring data, and determine the pipe burst position.

[0077] Specifically, the above step S206 includes:

[0078] Step S2061, based on the transient flow model, using the burst coefficient, simulate and calculate the transient flow pressure data of each possible burst position. The transient flow model can be used to simulate the instantaneous change of water flow pressure in the water supply network. For example, the transient flow model of the water supply network can be constructed using TSNet software. When constructing, the topological model of the water supply network can be constructed first, and then the transient conditions and solver parameters can be set to obtain the transient flow model. For the constructed transient flow model, transient flow simulation can be performed, such as calculating and displaying the dynamic parameters such as pressure and flow rate of the water supply network that change with time. In this embodiment, the burst coefficient c corresponding to the possible burst position k is calculated. k That is, the burst size is used as the transient condition simulated by the transient flow model, so the possible burst position k is calculated by the transient flow model and the size is c k When the pipe bursts, the pressure change time series h of each monitoring point i,cal (k) Transient flow pressure data.

[0079] Step S2063, using a dynamic time warping algorithm to match the transient flow pressure data of each monitoring point with the high-frequency pressure monitoring data, and determine the burst position based on the matching results. Among them, the dynamic time warping (Dynamic TimeWarping, DTW) algorithm is an algorithm for measuring the similarity between two sequences, and can process time series data with inconsistent lengths. When using the DTW algorithm, the two sequences are no longer one-to-one corresponding to each other at the same time as using the Euclidean distance, but points in one sequence are allowed to correspond to other time points in the other sequence to eliminate the influence of shape stretching or compression. Therefore, the dynamic time adjustment algorithm is used to match the pressure data of each monitoring point with the new pressure monitoring data, which can be expressed by the following formula:

[0080]

[0081] Therefore, through the above formula, the DTW algorithm is used to calculate the pressure change time series and steady-state head monitoring value H at each possible blasting position. i,obs The possible pipe burst position corresponding to the minimum distance is taken as the final located pipe burst position.

[0082] Compared with the method for locating a pipe burst in a water supply network provided by the embodiment of the present invention, the method utilizes high-frequency pressure monitoring data, and has better timeliness and positioning effect. The present invention not only uses the arrival time of the pressure wave to locate the pipe burst, but also couples the transient flow model of the pipe network to utilize the continuous change information of the pressure wave, thereby improving the accuracy of pipe burst positioning.

[0083] As a specific application example of the embodiment of the present invention, Figure 2 and Figure 3 The pipe burst location monitoring of the water supply network shown is used as an example for explanation. For the water supply network, its nodes are first determined, and high-frequency pressure monitoring data is collected at the locations of nodes J57, J100, J151, J193 and J365, with a sampling frequency of 250Hz. If relevant technologies such as only using the arrival time of the pressure wave for positioning are used, the result is that the pipe burst occurred at node J147. Other points where pipe bursts may occur are ranked from most to least likely: J158, J429, J124, J144... If the method provided in the embodiment of the present invention is used, the pipe burst can be accurately located at node J124 after analysis using the information of the continuous change of the pressure wave.

[0084] In this embodiment, a water supply network burst pipe locating device is also provided, which is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0085] This embodiment provides a water supply network burst pipe positioning device, such as Figure 4 As shown, including:

[0086] The data acquisition module 41 is used to acquire high-frequency pressure monitoring data of each monitoring point in the water supply network, wherein the high-frequency pressure monitoring data includes pressure values ​​and monitoring times corresponding to each pressure value;

[0087] A pipe burst monitoring module 42, for determining the time when the pipe burst pressure wave reaches each monitoring point by using a cumulative sum algorithm based on high-frequency pressure monitoring data;

[0088] A preliminary positioning module 43 is used to determine the possible burst position of the pipe based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network;

[0089] A pipe burst coefficient determination module 44 is used to match the new steady-state pressure monitoring data obtained after the pipe burst with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model to determine the pipe burst coefficient of the possible pipe burst position;

[0090] The positioning module 45 is used to simulate and calculate the transient flow pressure data of pipe burst at each possible pipe burst position based on the transient flow model and the pipe burst coefficient, match the transient flow pressure data with the high-frequency pressure monitoring data, and determine the pipe burst position.

[0091] In an optional embodiment, the pipe burst monitoring module is specifically used to: use a cumulative sum algorithm to calculate the moment when the accumulated pressure change is greater than a threshold based on the high-frequency pressure monitoring data, and determine that a pipe burst pressure wave occurs when the accumulated pressure change is greater than the threshold; determine the moment when the pipe burst pressure wave arrives at each monitoring point based on the moment when the pressure change recorded in the cumulative sum algorithm begins to accumulate.

[0092] In an optional embodiment, the device also includes: a time calculation module, which is used to determine the nodes of the water supply network and the pipes between the nodes based on the topological structure of the water supply network; and calculate the propagation time of the transient flow pressure wave in any pipe based on the length of the pipe between the two nodes and the transient flow wave velocity of the pipe.

[0093] In an optional implementation, the preliminary determination module is specifically used to: use the shortest path algorithm to determine the shortest time for pressure wave propagation between any two nodes based on the propagation time of transient flow pressure waves in any pipe of the water supply network; determine the objective function based on the minimum value of the difference between the time when the burst pressure wave reaches any two monitoring points and the difference in the shortest time for pressure wave propagation from any node to the corresponding two monitoring points, and solve the objective function to determine the possible burst location.

[0094] In an optional implementation, the objective function is expressed using the following formula:

[0095]

[0096] In the formula, k represents any node of the water supply network, M represents the set of monitoring points, and t i represents the time when the burst pressure wave reaches the monitoring point i, t j represents the time when the burst pressure wave reaches the monitoring point j, T ik represents the shortest time for transient flow pressure wave propagation between monitoring point i and node k, T jk It represents the shortest time for transient flow pressure wave to propagate between monitoring point j and node k.

[0097] In an optional implementation, the pipe burst coefficient determination module is specifically used to: obtain new steady-state pressure monitoring data of each monitoring point after a pipe burst occurs; calculate pressure estimation data of possible pipe burst locations under different pipe burst coefficients based on a pre-constructed steady-state model of the water supply network; determine an objective function based on the difference between the new steady-state pressure monitoring data and the pressure estimation data, solve the objective function using a genetic algorithm, and determine the pipe burst coefficient of each possible pipe burst location.

[0098] In an optional implementation, the positioning module is specifically used to: based on the transient flow model, use the burst coefficient to simulate and calculate the transient flow pressure data of each possible burst position; use the dynamic time warping algorithm to match the transient flow pressure data of each monitoring point with the high-frequency pressure monitoring data, and determine the burst position according to the matching result.

[0099] The further functional description of each of the above modules is the same as that of the above corresponding embodiments and will not be repeated here.

[0100] The embodiment of the present invention also provides a computer device having the above Figure 4 The water supply network burst pipe locating device shown.

[0101] See also Figure 5 , Figure 5 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Figure 5As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 5 A processor 10 is taken as an example.

[0102] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0103] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0104] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created by the use of a computer device based on the presentation of a small program landing page, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0105] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive; the memory 20 may also include a combination of the above types of memory.

[0106] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0107] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0108] A part of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the existence of the computer program instruction in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc., and accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium accessible to the computer.

[0109] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for locating a burst pipe in a water supply network, characterized in that: The method comprises: Acquire high-frequency pressure monitoring data of each monitoring point in the water supply network, wherein the high-frequency pressure monitoring data includes pressure values ​​and monitoring times corresponding to each pressure value; Based on the high-frequency pressure monitoring data, a cumulative sum algorithm is used to determine the time when the burst pipe pressure wave reaches each monitoring point; Determine the possible burst position based on the relationship between the transient flow pressure wave propagation time and the above mentioned moments between any two nodes in the water supply network; Match the new steady-state pressure monitoring data obtained after the pipe burst with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model to determine the pipe burst coefficient of the possible pipe burst location; Based on the transient flow model and using the pipe burst coefficient, transient flow pressure data of pipe bursts at each possible pipe burst position are simulated and calculated, and the transient flow pressure data is matched with the high-frequency pressure monitoring data to determine the pipe burst position.

2. The method according to claim 1, characterized in that Based on the high-frequency pressure monitoring data, a cumulative sum algorithm is used to determine the time when the burst pressure wave reaches each monitoring point, including: A cumulative sum algorithm is used to calculate the time when the pressure variation accumulation is greater than a threshold value based on the high-frequency pressure monitoring data. If the pressure variation accumulation is greater than the threshold value, it is determined that a pipe burst pressure wave occurs; The time when the pressure change recorded in the cumulative sum algorithm starts to accumulate is used to determine the time when the burst pipe pressure wave reaches each monitoring point.

3. The method according to claim 1, characterized in that Before determining the possible pipe burst position based on the relationship between the transient flow pressure wave propagation time between any two nodes in the water supply network and the time instants, the method further includes: Based on the topological structure of the water supply network, determine the nodes of the water supply network and the pipes between the nodes; The propagation time of transient flow pressure waves in any pipe is calculated based on the length of the pipe between two nodes and the transient flow wave velocity of the pipe.

4. The method according to claim 1, characterized in that: Determine the possible burst position based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network, including: The shortest path algorithm is used to determine the shortest time for pressure wave propagation between any two nodes based on the propagation time of transient flow pressure waves in any pipe of the water supply network. The objective function is determined based on the minimum value of the difference between the time when the burst pressure wave reaches any two monitoring points and the difference between the shortest propagation time of the pressure wave from any node to the corresponding two monitoring points, and the possible burst position is determined by solving the objective function.

5. The method according to claim 4, characterized in that The objective function is expressed by the following formula: In the formula, k represents any node of the water supply network, M represents the set of monitoring points, and t i represents the time when the burst pressure wave reaches the monitoring point i, t j represents the time when the burst pressure wave reaches the monitoring point j, T ik represents the shortest time for transient flow pressure wave propagation between monitoring point i and node k, T jk It represents the shortest time for transient flow pressure wave to propagate between monitoring point j and node k.

6. The method according to claim 1, characterized in that The new steady-state pressure monitoring data obtained after the pipe burst is matched with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model to determine the pipe burst coefficient of the possible pipe burst location, including: Obtain new steady-state pressure monitoring data of each monitoring point after pipe burst occurs; Calculate the pressure estimation data of possible burst locations under different burst coefficients based on the pre-built steady-state model of the water supply network; The objective function is determined based on the difference between the new steady-state pressure monitoring data and the pressure estimation data, and the objective function is solved by using a genetic algorithm to determine the pipe burst coefficient of each possible pipe burst position.

7. The method according to claim 1, characterized in that Based on the transient flow model, using the pipe burst coefficient, simulating and calculating the transient flow pressure data of pipe bursts at each possible pipe burst position, matching the transient flow pressure data with the high-frequency pressure monitoring data, and determining the pipe burst position, including: Based on the transient flow model, using the pipe burst coefficient, simulate and calculate the transient flow pressure data of pipe burst at each possible pipe burst position; The transient flow pressure data of each monitoring point is matched with the high-frequency pressure monitoring data using a dynamic time warping algorithm, and the burst pipe position is determined based on the matching result.

8. A water supply pipe burst positioning device, characterized in that: The device comprises: A data acquisition module, used to acquire high-frequency pressure monitoring data of each monitoring point in the water supply network, wherein the high-frequency pressure monitoring data includes pressure values ​​and monitoring times corresponding to each pressure value; A pipe burst monitoring module, used to determine the time when the pipe burst pressure wave reaches each monitoring point by using a cumulative sum algorithm based on the pressure monitoring data; A preliminary positioning module, used to determine the possible burst position of the pipe based on the relationship between the transient flow pressure wave propagation time and each time between any two nodes in the water supply network; A pipe burst coefficient determination module is used to match the new steady-state pressure monitoring data obtained after the pipe burst with the pressure estimation data corresponding to different pipe burst coefficients determined based on the steady-state model to determine the pipe burst coefficient of the possible pipe burst position; The positioning module is used to simulate and calculate the transient flow pressure data of pipe burst at each possible pipe burst position based on the transient flow model and the pipe burst coefficient, match the transient flow pressure data with the high-frequency pressure monitoring data, and determine the pipe burst position.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the water supply network burst pipe locating method according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for locating a burst pipe in a water supply network according to any one of claims 1 to 7.

Images