CN120368231B - A method, device and system for detecting and locating leakage in multiple pipes of a water supply network - Google Patents
A method, device and system for detecting and locating leakage in multiple pipes of a water supply networkInfo
- Publication number
- CN120368231B CN120368231B CN202510833601.9A CN202510833601A CN120368231B CN 120368231 B CN120368231 B CN 120368231B CN 202510833601 A CN202510833601 A CN 202510833601A CN 120368231 B CN120368231 B CN 120368231B
- Authority
- CN
- China
- Prior art keywords
- leakage
- data
- node
- pipeline
- pipe
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Examining Or Testing Airtightness (AREA)
Abstract
A method for detecting and positioning multi-pipeline leakage of a water supply network includes the steps of constructing a network model and sensor deployment to create a complex hydraulic network model, screening optimal monitoring points to generate a pressure and flow sensor arrangement scheme, outputting a network topological graph and a sensor deployment result, performing scene approximation, constructing an optimization model to enable hydraulic data under original normal working conditions to fit hydraulic data under a leakage scene, finally introducing hydraulic data correction to output possible leakage nodes and corresponding hydraulic parameters, reducing the range of the leakage network, continuously updating and screening the pipeline range in which leakage is possible based on preliminary deduction of the leakage nodes, outputting a local pipeline area in which leakage is most likely to occur after iterative optimization, building an observer and positioning specific leakage points, dynamically estimating leakage parameters in suspected leakage pipeline sections, comprehensively analyzing filtering results, and outputting specific leakage pipeline sections and leakage position coordinates.
Description
Technical Field
The invention belongs to the field of multi-position leakage detection of urban water supply networks, and particularly relates to a detection and positioning method, device and system for coping with multi-pipeline leakage of a water supply network.
Background
At present, a water supply network still strives to solve the problems of operation length and water leakage at a plurality of positions. The multi-pipeline leakage of the complex pipe network is difficult to detect and position, affects the daily water consumption of urban residents, and causes the waste of natural resources and the loss of water enterprises. The positioning accuracy of multiple leakage losses of the complex water supply network is improved, and the method becomes a social problem to be solved urgently.
In the field of water supply network leakage detection, the equipment method is an advanced technology. The technology is mainly used for detecting leakage of the pipeline network by using various detection instruments. However, the accuracy of the detection result thereof depends largely on the experience of the operator. Moreover, this approach often requires significant time and effort to handle large pipe network leak detection. With the development of artificial intelligence and sensors, a data driving method has been developed, and the chinese patent publication CN115994487a "a water supply pipe network leakage positioning method based on correlation coefficients" proposes a water supply pipe network leakage positioning method based on an improved gray wolf optimization algorithm, which requires a large number of flowmeter arrangements under the conditions of annular and complex pipe networks, and has a large economic cost. The Chinese patent publication CN119123342A provides a method and a system for monitoring leakage of a town water supply network, which predict water consumption through a neural network, but the model requires a large amount of historical hydraulic data for training, has weak applicability to complex network in reality and costs a large amount of calculation cost.
Disclosure of Invention
In order to solve the defects in the prior art, the invention provides a method, a device and a system for detecting and positioning multi-pipeline leakage of a water supply network.
The first aspect of the invention relates to a detection and positioning method for coping with multi-pipeline leakage of a water supply network, which comprises the following specific steps:
s1, constructing a pipe network model and deploying sensors;
The system receives hydraulic data input by a user, analyzes the hydraulic data through a mass and energy conservation equation, combines node and pipeline information to create a complex hydraulic pipe network model, screens optimal monitoring points through calculating a correlation matrix among the nodes and combining a greedy algorithm and information redundancy degree to generate a pressure and flow sensor arrangement scheme, and outputs a pipe network topological graph and a sensor deployment result.
S2, performing scene approximation;
The method comprises the steps of judging whether a pipe network is leaked or not basically, if a sensor measured value and a pipe network model simulation value have large differences, considering that the pipe network is leaked, determining flow and pressure sensors sensitive to the change of measured data according to the relative fluctuation values of arranged sensors, determining optimization variables by combining a correlation matrix and information redundancy, constructing an optimization model to enable hydraulic data under original normal working conditions to fit the hydraulic data under a leakage scene, finally introducing hydraulic data correction, and outputting possible leakage nodes and corresponding hydraulic parameters.
S3, reducing the range of a leakage pipe network;
Adopting a node expansion strategy and a depth-first search algorithm, and continuously updating and screening a pipeline range with possible leakage based on preliminary inference of the leakage node; after iterative optimization, the local pipe network area where leakage is most likely to occur is output.
S4, establishing an observer and positioning a specific leakage point;
And constructing an observer model by using an extended Kalman filter, fusing multi-source observation information, dynamically estimating leakage parameters in the suspected leakage pipe section, comprehensively analyzing the filtering result, and outputting specific leakage pipe section and leakage position coordinates.
Further, in step S1, constructing the pipe network model and the sensor deployment specifically includes:
S11, constructing a node continuous equation and a pipe section energy equation in a pipe network according to the hydraulic data input by the system and based on the mass conservation and energy conservation principles so as to accurately describe the hydraulic behavior of the system;
S12, establishing a finite dimension mathematical model for each pipe section in the pipe network, discretizing a continuous hydraulic process to support subsequent sensitivity analysis and monitoring point optimization calculation;
s13, calculating a pressure sensitivity matrix among nodes based on a perturbation analysis method, and generating a binary correlation matrix through normalization and threshold processing The pressure correlation degree measuring device is used for measuring the pressure correlation degree between the nodes;
s14, initializing a sensor selection process, setting a selected node set as null, and setting a node set to be covered as a correlation matrix All columns in (a);
s15, under the condition that the column set to be covered is not empty, traversing all rows of the correlation matrix, searching for rows capable of covering the most unselected columns, adding the corresponding nodes of the rows into the selected set, and updating the column set to be covered according to the corresponding nodes until all columns are covered or no feasible selection exists;
And S16, returning the finally selected row index, wherein the corresponding node is the recommended position for installing the pressure sensor, and arranging a flow sensor on the pipe section connected with the nodes so as to realize accurate sensing of leakage or abnormal working conditions. The data measured by the sensor is combined into a joint observation vector WhereinIndicating a flow meter reading,Indicating the reading of the pressure gauge,And (3) withThe number of flow and pressure sensors respectively, and the total observation dimension is。
Further, in step S2, determining the leakage area and performing the scene approximation specifically includes:
S21, after the sensor arrangement is determined by utilizing a sensor deployment algorithm, residual errors between the hydraulic characteristics of the simulated pipe network and the actual measured values of the sensors can be calculated, if the large difference exists, the pipe network is considered to have leakage, and the following scene approximation method is carried out;
s22, calculating a relative fluctuation value according to the residual error in S21, selecting two pressure sensors and flow sensors with larger relative fluctuation value, calculating information redundancy of related nodes of the sensors, changing the selected sensors, and finally, according to the related nodes and the related matrix of the sensors Determining an optimization variable;
and S23, initializing an additional flow vector. Initializing an extra traffic vector of each node represented by the optimization variable to be a 0 vector;
S24, calculating pipe network characteristics, and calculating observation vector estimation in the current state WhereinThe time step in which the model is run is represented,Representing the observation vector atThe model estimate after a time step,The function represents a hydraulic calculation function based on Todini gradient algorithm,Representing the node's best list of additional traffic. Estimating an observation vectorNode indexMatching to obtain updated pressure estimate as;
S25, defining an optimization objective function;
S26, solving an optimization objective function by using a CSA color-changing Long Qun optimization algorithm to obtain an optimal extra-flow list ;
S27, using the optimumThe simulation is updated again for a plurality of steps to obtain the observation vector estimation approaching the leakage stateOptimal additional water demand list for each relevant node;
S28, hydraulic characteristic adjustment, namely after the leaking pipe section is identified, weight distribution is carried out on upstream and downstream flow so as to fully approximate the leaking scene;
s29, outputting final observation vector estimation approaching the leakage state And the optimal additional water demand for each relevant node;
Further, in step S3, the narrowing the range of the leakage pipe network specifically includes:
S31 initializing variable, leakage flow list For the calculation result in the S2 method, the head-tail node list of leakage of the leakage pipe sectionInitializing to empty list, initializing to tag access node list。
S32, traversing the extra flow listA traffic value of each node in the plurality if the traffic value of the node is greater than a threshold valueAdding the traffic value of the node into a leakage traffic listIndicating that there is a potential for leakage at the node;
s33, traversing the leakage flow list Performing a depth-first search (DFS) of the non-accessed leaky nodes to determine possible leaky nodes connected thereto;
And S34, returning to a list of nodes at two ends of all the leakage pipe sections.
Further, in step S4, establishing the observer and locating the specific leakage point specifically includes:
s41, determining a state space expression with unique leakage under a single pipe section;
firstly, determining the hydraulic boundary condition of a target pipe section, taking the pressure of nodes at two ends of the pipe section as input signals, and respectively recording And (3) withThe flow rates at the two ends are taken as system observation output signals and respectively recorded asAnd (3) with. The system state variables are set as follows:
Wherein, the Representing the inlet flow of the pipe section, wherein the unit is m3/s; Representing the pressure of the intermediate node in m (meter water column); Represents the outlet flow of the pipe section, and the unit is m3/s; the relative position of leakage in the pipe section is represented, the value range is 0 to 1, and no unit exists; the leakage flow rate per unit time is expressed in m3/s.
Defining system input variables as:
defining the system observation output as:
the continuous time state equation form of the system is as follows:
Wherein, the The function is expressed in the current stateAnd inputThe derivative of the state variable.
S42, discretizing system variables;
The continuous time state equation was discretized using Heun's method (a modified trapezoidal method). Set the time step as First, theThe status updates for the individual time steps are:
the formula may approximately describe the state change of a nonlinear system over a short period of time. Thereby defining a discrete state transfer function of the system:
s43, initializing state variable estimated values and covariance matrixes;
Setting an initial state estimate as:
Wherein the method comprises the steps of Indicating the desire. The error covariance matrix is initialized to:
s44, state prediction;
For each time instant Estimating a current state using a predictive model:
simultaneous prediction error covariance matrix:
Wherein, the A jacobian matrix representing a system state function; Representing a process noise covariance matrix.
S45, measuring and updating, and calculating residual errors;
Calculating an observation residual error:
Wherein the method comprises the steps of Is an observation matrix, here:
S46, state estimation correction;
Correcting state estimation based on residual errors:
Wherein the Kalman gain matrix is:
and updating the error covariance:
Wherein the method comprises the steps of Is the observed noise covariance matrix and,Is an identity matrix.
S47, outputting an estimation result;
through the continuous iteration of the state prediction and correction, the state variable is finally obtained And (3) withI.e. the location and leakage intensity of the only leakage point in the pipe section.
A second aspect of the present invention relates to a leakage detecting and positioning device for a water supply pipe, configured to a data processing and analyzing apparatus, comprising:
the pipeline monitoring data acquisition module is used for acquiring pipeline related data sent by monitoring equipment, wherein the monitoring data comprises water level data, water flow data and water pressure data of each node unit basin in the current monitoring area;
The pipeline leakage analysis result determining module is used for carrying out forecast evolution calculation and pipeline leakage analysis on the monitoring point data to obtain a leakage analysis result of the current monitoring area;
And the early warning module is used for generating early warning information and related data of the leakage pipeline, wherein the data comprise leakage point pipeline data, water flow data, water pressure data and position data. And sending the early warning information to early warning equipment so that the early warning equipment informs relevant maintenance personnel according to the early warning information.
The third aspect of the invention relates to a leakage detection and positioning system of a water supply pipe, which is characterized by comprising pipe network monitoring equipment, data acquisition equipment, data processing and analysis equipment, terminal equipment and early warning equipment, wherein the data acquisition equipment, the early warning equipment and the terminal equipment are connected with the data processing and analysis equipment through communication equipment;
The pipe network monitoring equipment is used for monitoring the flow and pressure states of the pipe network;
The data acquisition equipment is used for acquiring water pressure data, water level data and water flow data of each node unit basin in the area;
The data processing analysis equipment is used for acquiring data sent by the pipe network data acquisition equipment, estimating water demand based on a scene approximation method, reducing a pipe network leakage area by using a depth-first algorithm, and then carrying out specific leakage positioning analysis by using an extended Kalman observer so as to obtain a leakage detection result in a corresponding area;
The early warning equipment is used for generating early warning information including leakage point pipeline data, water flow data, water pressure data and position data if the leakage detection analysis result obtained by the data processing analysis equipment is confirmed to be in a leakage state, and sending the early warning information to the early warning equipment so that the early warning equipment informs related maintenance personnel according to the early warning information;
And the terminal equipment is used for receiving the analysis data and the pipeline leakage analysis result sent by the data processing and analysis equipment and visualizing the analysis data and the pipeline leakage analysis result.
The innovation points of the invention are as follows:
The invention is based on the observer technology (extended Kalman filter), and is carried out in three steps, firstly, the sensor arrangement is carried out on a pipe network by utilizing a pressure sensitivity matrix, then, a leakage area and a leakage pipe section are respectively determined by using a threshold method and a scene approximation method, the size of a possible leakage pipe network is reduced by a depth-first search algorithm, finally, an observer is designed on the leakage pipe section to determine leakage parameters, the leakage range is gradually reduced to determine a specific leakage pipe section, and the scale of the observer problem is obviously reduced. The detection method improves the leak detection precision and the leak detection efficiency, reduces the labor cost and helps staff find and overhaul in time.
Drawings
Fig. 1 is a schematic general flow diagram of a method for detecting and locating multiple leaks in a water supply network according to a first embodiment of the present invention.
Fig. 2 is a schematic diagram of specific implementation steps of a method for detecting and locating multiple leaks in a water supply network according to a first embodiment of the present invention.
FIG. 3 is a schematic diagram of a simulation of pipe leak detection and sensor arrangement in accordance with a first embodiment of the present invention.
Fig. 4 is a schematic structural diagram of a leakage detecting and positioning device for a water supply pipe according to a second embodiment of the present invention.
Fig. 5 is a schematic structural diagram of a water supply pipe leakage detecting and positioning system according to a third embodiment of the present invention.
Detailed Description
The following describes the detailed implementation of the embodiments of the present invention with reference to the drawings. It should be understood that the detailed description and specific examples, while indicating and illustrating the invention, are not intended to limit the invention.
The invention will be described in detail below with reference to the drawings in connection with exemplary embodiments.
Example 1
Referring to fig. 1, 2 and 3, the method for detecting and positioning the multi-pipeline leakage of the water supply network comprises the following steps:
The method comprises the steps of S1, constructing a pipe network model and sensor deployment, receiving hydraulic data input by a user, analyzing through a mass and energy conservation equation, combining node and pipeline information to create a complex hydraulic pipe network model, calculating a correlation matrix among the nodes, combining a greedy algorithm and information redundancy degree, screening optimal monitoring points to generate a pressure and flow sensor deployment scheme, and outputting a pipe network topological diagram and a sensor deployment result.
Specifically, step S1 includes:
S11, constructing a node continuous equation and a pipe section energy equation in a pipe network according to the hydraulic data input by the system and based on the mass conservation and energy conservation principles so as to accurately describe the hydraulic behavior of the system;
The present example uses the Hazen-Williams (H-W) equation to calculate head loss. The effect of the pipe parameters in this formula is summarized at the H-W drag coefficient for each pipe In which the coefficient is oneIs defined as the pipeline parameter function of. Then, calculate the pipeline through H-W formulaThe head loss is as follows:
Wherein the method comprises the steps of Is a constant index of the values of the three,Is a pipelineIs a flow rate of (a).
Pumps are another important element in water supply networks and they are characterized by pressure flow curves that are used to correlate pump flow to pump head gain according to the specifications of each pump. Described by the following equation:
Wherein the method comprises the steps of ,AndIs the coefficient of the pump pressure flow curve, and the negative sign indicates that the head of the pump is gain rather than loss.
Thus, the energy conservation equation in the water supply network can be written as follows:
Wherein the method comprises the steps of Is the water head of an unknown node,,Matrix representing water flow direction relationship between nodes, ifIllustrating that the water flow is a slave nodeFlow direction nodeOtherwise ifThen representAndThe water flow direction relationship between the nodes is fromFlow direction. If it isThen describe the nodeAndNo direct connection is available.Is in each equationThe sum of the known heads that occur. In the vector representation, the vector of the known head is represented byGiven.
The mass conservation equation in the water supply network can be written as follows:
S12, establishing a finite dimension mathematical model for each pipe section in the pipe network, discretizing a continuous hydraulic process to support subsequent sensitivity analysis and monitoring point optimization calculation;
Will be Pipeline node at momentFlow rate in (a)Is marked as. Similarly, it willPipeline node at momentPressure in (a)Is marked as. By usingAndRepresentative ofAnd. The difference between two adjacent pipeline nodes is recorded as。Is the local acceleration caused by gravity.Is the cross-sectional area of the conduit.Is the inner diameter of the pipeline. For any number of pipe segments, a finite dimensional model may be obtained:
Wherein the method comprises the steps of 、、AndRespectively corresponding to boundary conditions、、And. Coefficient of frictionWhich itself depends on the flow rate, which will make sense when considering the subsequent leakage situation.The model of (2) is:
Wherein the method comprises the steps of Indicating the roughness of the pipe and,For Reynolds number, the calculation formula is as follows: Wherein Is the kinematic viscosity of water and is used for the preparation of a liquid,Is the firstFlow of the segment pipe section.
S13, calculating a pressure sensitivity matrix among nodes based on a perturbation analysis method, and generating a binary correlation matrix through normalization and threshold processingThe pressure correlation degree measuring device is used for measuring the pressure correlation degree between the nodes;
based on pipe network microscopic analysis method, calculation Node pairThe sensitivity of the nodes is as follows:
Wherein the method comprises the steps of ,Node for requiring water under reference working conditionAndIs used for controlling the water pressure of the water,,Is thatAfter the node flow is changed,The water pressure of the node.
For the firstThe element values of the columns are subjected to minimum and maximum normalization, so that the node pairs of all the nodes can be obtainedDimensionless sensitivity matrix for node in need of water change,Is the first of (2)Line 1The column elements are:
Wherein the method comprises the steps of AndRespectively represent sensitivity matrixIs the first of (2)Minimum and maximum values of columns. In order to express the strength of the correlation, a correlation threshold is set
S14, initializing a sensor selection process, setting a selected node set as null, and setting a node set to be covered as a correlation matrixAll columns in (a);
s15, under the condition that the column set to be covered is not empty, traversing all rows of the correlation matrix, searching for rows capable of covering the most unselected columns, adding the corresponding nodes of the rows into the selected set, and updating the column set to be covered according to the corresponding nodes until all columns are covered or no feasible selection exists;
And S16, returning the finally selected row index, wherein the corresponding node is the recommended position for installing the pressure sensor, and arranging a flow sensor on the pipe section connected with the nodes so as to realize accurate sensing of leakage or abnormal working conditions. The data measured by the sensor is combined into a joint observation vector WhereinIndicating a flow meter reading,Indicating the reading of the pressure gauge,And (3) withThe number of flow and pressure sensors respectively, and the total observation dimension is。
S2, performing scene approximation, namely performing basic judgment on whether the pipe network is leaked or not, and if the measured value of the sensor arranged in the S1 and the simulation value of the pipe network model have larger differences, considering that the pipe network is leaked, introducing a scene approximation method, and determining two flow and pressure sensors according to the relative fluctuation value of the arranged sensors. And determining an optimization variable by combining the correlation matrix and the information redundancy, constructing an optimization model to enable the hydraulic data under the original normal working condition to fit the hydraulic data under the leakage scene, and finally introducing hydraulic data correction, and outputting possible leakage nodes and corresponding hydraulic parameters.
Specifically, step S2 includes:
S21, after the sensor arrangement is determined by utilizing a sensor deployment algorithm, residual errors between the hydraulic characteristics of the simulated pipe network and the actual measured values of the sensors can be calculated, if the large difference exists, the pipe network is considered to have leakage, and the following scene approximation method is carried out;
s22, calculating a relative fluctuation value according to the residual error in S21, selecting two pressure sensors and flow sensors with larger relative fluctuation value, calculating information redundancy of related nodes of the sensors, changing the selected sensors, and finally, according to the related nodes and the related matrix of the sensors Determining an optimization variable;
and S23, initializing an additional flow vector. Initializing an extra traffic vector of each node represented by the optimization variable to be a 0 vector;
S24, calculating pipe network characteristics, and calculating observation vector estimation in the current state WhereinThe time step in which the model is run is represented,Representing the observation vector atThe model estimate after a time step,The function represents a hydraulic calculation function based on Todini gradient algorithm,Representing the node's best list of additional traffic. Estimating an observation vectorNode indexMatching to obtain updated pressure estimate as;
S25, defining an optimization objective function, wherein the scene of the water supply network under the normal working condition is defined asThe scene when the pipe network is leaked is defined as a leakage sceneThe purpose of the scene approximation is to enable the node pressure and the pipe section flow characteristics in the pipe network to be equal to those of the leakage scene by changing the water demand of the related water-demand nodes under normal conditionsApproximately equal.
The sensor is noted asThese are provided withThe node with strong sensor correlation is,. Is provided withThe additional water demand at the node isAnd define. The scene approximation problem can be described by the following optimization model:
Wherein the method comprises the steps of Is thatThe measured values of the flow and pressure sensors in the pipe network at the moment,Is thatApplying at the corresponding pipe network node at any timeThen, solving the obtained pipe network flow pressure characteristics by using a pipe network characteristic calculation methodPipe section flow and node pressure with sensors installed therein are noted as。Is a weight matrix of sensors, and if the measurements between each sensor are independent, then the matrixIs a diagonal matrix and a larger value of a diagonal element indicates a higher accuracy of the corresponding sensor.Is the total water leakage of the water supply network area becauseFluctuation value with uncertainty so that leak amount is set。
S26, solving an optimization objective function by using CSA (color change Long Qun optimization algorithm) to obtain an optimal additional flow list;
S27, using the optimumThe simulation is updated again for a plurality of steps to obtain the observation vector estimation approaching the leakage stateAnd the optimal additional water demand for each relevant node;
S28, hydraulic characteristic adjustment;
The flow pressure characteristics of the leakage scene are approximated only by changing the water demand of the nodes, but the upstream and downstream flows cannot be approximated in the leakage pipe section, so that the upstream and downstream flows are required to be weighted after the leakage pipe section is identified, so that the leakage scene is approximated fully.
Is provided withIs thatThe pressure difference between the nodes at the two ends of the pipe section,Is thatThe diameter of the tube section,The multiple leaks relate to the same node but do not include the number of pipe sections with flow meters installed. According to Darcy-Wei Siba hz equation (Darcy-Weisbach equation), the flow in the pipe is related to the pressure difference at the nodes at the ends of the pipe section, and the pipe diameter is set with the following distribution weights:
related pipe section after weight distribution The flow calculation form of (a) is as follows:
s29, outputting final observation vector estimation approaching the leakage state And the optimal additional water demand for each relevant node;
S3, reducing the range of the leakage pipe network, adopting a node expansion strategy and a depth-first search algorithm, continuously updating and screening the range of the pipeline with the possible leakage based on the preliminary inference of the leakage nodes, and outputting the local pipe network area with the most possible leakage after iterative optimization.
Specifically, step S3 includes:
S31 initializing variable, leakage flow list For the calculation result in the S2 method, the head-tail node list of leakage of the leakage pipe sectionInitializing to empty list, initializing to tag access node list。
S32, traversing the extra flow listA traffic value of each node in the plurality if the traffic value of the node is greater than a threshold valueAdding the traffic value of the node into a leakage traffic listIndicating that there is a potential for leakage at the node;
s33, traversing the leakage flow list Performing depth-first search DFS on the unaccessed leaky nodes to determine possible leaky nodes connected to the DFS;
And S34, returning to a list of nodes at two ends of all the leakage pipe sections.
S4, establishing an observer and positioning specific leakage points, utilizing an extended Kalman filter to establish an observer model, fusing multi-source observation information, dynamically estimating leakage parameters in suspected leakage pipe sections, comprehensively analyzing filtering results, and outputting specific leakage pipe sections and leakage position coordinates.
Specifically, step S4 includes:
s41, determining a state space expression with unique leakage under a single pipe section;
firstly, determining the hydraulic boundary condition of a target pipe section, taking the pressure of nodes at two ends of the pipe section as input signals, and respectively recording And (3) withThe flow rates at the two ends are taken as system observation output signals and respectively recorded asAnd (3) with. The system state variables are set as follows:
Wherein, the Representing the inlet flow (unit: m 3/s) of the pipe section; representing the pressure of the intermediate node (unit: m); Represents the pipe section outlet flow (unit: m 3/s); indicating the relative position of the leak in the pipe section (value range 0 to 1, no unit); the leakage flow rate per unit time (unit: m 3/s) is shown.
Defining system input variables as:
defining the system observation output as:
the continuous time state equation form of the system is as follows:
Wherein, the The function is expressed in the current stateAnd inputThe derivative of the state variable.
S42, discretizing system variables;
The continuous time state equation was discretized using Heun's method (a modified trapezoidal method). Set the time step as First, theThe status updates for the individual time steps are:
the formula may approximately describe the state change of a nonlinear system over a short period of time. Thereby defining a discrete state transfer function of the system:
s43, initializing state variable estimated values and covariance matrixes;
Setting an initial state estimate as:
Wherein the method comprises the steps of Indicating the desire. The error covariance matrix is initialized to:
s43, initializing state variable estimated values and covariance matrixes;
Setting an initial state estimate as:
Wherein the method comprises the steps of Indicating the desire. The error covariance matrix is initialized to:
s44, state prediction;
For each time instant Estimating a current state using a predictive model:
simultaneous prediction error covariance matrix:
Wherein, the A jacobian matrix representing a system state function; Representing a process noise covariance matrix.
S45, measuring and updating, and calculating residual errors;
Calculating an observation residual error:
Wherein the method comprises the steps of Is an observation matrix, here:
S46, state estimation correction;
Correcting state estimation based on residual errors:
Wherein the Kalman gain matrix is:
and updating the error covariance:
Wherein the method comprises the steps of Is the observed noise covariance matrix and,Is an identity matrix.
S47, outputting an estimation result;
through the continuous iteration of the state prediction and correction, the state variable is finally obtained And (3) withI.e. the location and leakage intensity of the only leakage point in the pipe section.
Example 2
As shown in fig. 4, the present embodiment provides a device for detecting and positioning leakage of multiple pipelines of a water supply network, for implementing the method of embodiment 1, where the device is configured in a data processing and analyzing apparatus. The device comprises a pipeline monitoring data acquisition module, a pipeline leakage analysis result determination module and an early warning module.
The pipeline monitoring data acquisition module is used for acquiring pipeline related data sent by monitoring equipment, wherein the monitoring data comprises water level data, water flow data and water pressure data of each node unit basin in the current monitoring area;
The pipeline leakage analysis result determining module is used for carrying out forecast evolution calculation and pipeline leakage analysis on the monitoring point data to obtain a leakage analysis result of the current monitoring area;
And the early warning module is used for generating early warning information and related data of the leakage pipeline, wherein the data comprise leakage point pipeline data, water flow data, water pressure data and position data. And sending the early warning information to early warning equipment so that the early warning equipment informs relevant maintenance personnel according to the early warning information.
Example 3
As shown in fig. 5, the present embodiment provides a detection and positioning system for coping with multi-pipeline leakage of a water supply network, for implementing the method of embodiment 1, including a network monitoring device, a data acquisition device, a data processing analysis device, a terminal device and an early warning device, where the data acquisition device, the early warning device and the terminal device are connected with the data processing analysis device through a communication device;
The pipe network monitoring equipment is used for monitoring whether each node of the pipe network normally operates;
The data acquisition equipment is used for acquiring water pressure data, water level data and water flow data of each node unit basin in the area;
The data processing analysis equipment is used for acquiring data sent by the pipe network data acquisition equipment, estimating water demand based on a scene approximation method, reducing a pipe network leakage area by using a depth-first algorithm, and then carrying out specific leakage positioning analysis by using an extended Kalman observer so as to obtain a leakage detection result in a corresponding area;
The early warning equipment is used for generating early warning information including leakage point pipeline data, water flow data, water pressure data and position data if the leakage detection analysis result obtained by the data processing analysis equipment is confirmed to be in a leakage state, and sending the early warning information to the early warning equipment so that the early warning equipment informs related maintenance personnel according to the early warning information;
the terminal equipment is used for receiving the analysis data and the pipeline leakage analysis result sent by the data processing and analysis equipment and visualizing the analysis data and the pipeline leakage analysis result;
It should be noted that in the above embodiment of the pipeline leakage detection domain positioning device, the included modules are only divided according to the functional logic, and are not limited to the above division, as long as the corresponding functions can be implemented, and the specific names of the functional units are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present invention.
Finally, the above is merely an example of the invention and the technical principles applied. Those skilled in the art will appreciate that the invention is not limited to the specific embodiments described herein, as various obvious changes, rearrangements, and substitutions are possible without departing from the scope of the invention. Therefore, although the above embodiments have described the present invention in detail, the present invention is not limited to the above embodiments. Many other equivalent embodiments may be made without departing from the spirit and scope of the present invention. The scope of the invention is defined by the claims.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202510833601.9A CN120368231B (en) | 2025-06-20 | 2025-06-20 | A method, device and system for detecting and locating leakage in multiple pipes of a water supply network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202510833601.9A CN120368231B (en) | 2025-06-20 | 2025-06-20 | A method, device and system for detecting and locating leakage in multiple pipes of a water supply network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN120368231A CN120368231A (en) | 2025-07-25 |
CN120368231B true CN120368231B (en) | 2025-09-19 |
Family
ID=96441153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202510833601.9A Active CN120368231B (en) | 2025-06-20 | 2025-06-20 | A method, device and system for detecting and locating leakage in multiple pipes of a water supply network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN120368231B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101761780A (en) * | 2010-01-11 | 2010-06-30 | 中国石油大学(华东) | Gas pipeline leakage detecting and positioning device and method thereof |
CN116010800A (en) * | 2022-11-30 | 2023-04-25 | 南昌大学 | Leakage point analysis method for quantifying candidate leakage node detection priority |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101699359B1 (en) * | 2010-03-05 | 2017-01-24 | 인터디지탈 패튼 홀딩스, 인크 | Method and apparatus for providing security to devices |
CN116384277A (en) * | 2023-03-31 | 2023-07-04 | 杭州莱宸科技有限公司 | A filter-based multi-branch water supply pipe leakage detection and location method, device and system |
CN118998624A (en) * | 2024-07-03 | 2024-11-22 | 重庆大学 | Urban water supply network leakage monitoring method |
-
2025
- 2025-06-20 CN CN202510833601.9A patent/CN120368231B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101761780A (en) * | 2010-01-11 | 2010-06-30 | 中国石油大学(华东) | Gas pipeline leakage detecting and positioning device and method thereof |
CN116010800A (en) * | 2022-11-30 | 2023-04-25 | 南昌大学 | Leakage point analysis method for quantifying candidate leakage node detection priority |
Also Published As
Publication number | Publication date |
---|---|
CN120368231A (en) | 2025-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Goulet et al. | Model falsification diagnosis and sensor placement for leak detection in pressurized pipe networks | |
CN106195646B (en) | Leak is distributed estimating device, deduction system and presumption method | |
KR101105192B1 (en) | Water Supply Leak Segment Selection Method by Pipe Network Analysis and Integrated Waterworks Management System | |
JP6273125B2 (en) | Leakage investigation planning device, leakage investigation planning system, and leakage investigation planning method | |
Moser et al. | Performance comparison of reduced models for leak detection in water distribution networks | |
US20220196512A1 (en) | Detection of a leakage in a supply grid | |
CN109716090B (en) | Improved system for estimating water flow at subnet boundaries of water distribution network | |
CN107013812B (en) | A kind of THM coupling line leakage method | |
CN116906839B (en) | Intelligent monitoring and early warning method for thermal pipeline safety integrating physical measurement and soft measurement | |
WO2014115399A1 (en) | Water leak estimating device, system, and method | |
WO2020154468A1 (en) | Flow network model analysis | |
Moors et al. | Automated leak localization performance without detailed demand distribution data | |
Kutyłowska | Prediction of failure frequency of water-pipe network in the selected city | |
JP2006183274A (en) | Method and apparatus for predicting sewer pipe damage in sewer network | |
Hu et al. | Multi-objective and risk-based optimal sensor placement for leak detection in a water distribution system | |
CN115493093A (en) | Steam heating pipe network leakage positioning method and system based on mechanical simulation | |
JP7700108B2 (en) | Leak Detection Methods | |
Mazaev et al. | Probabilistic leak localization in water distribution networks using a hybrid data-driven and model-based approach | |
CN120368231B (en) | A method, device and system for detecting and locating leakage in multiple pipes of a water supply network | |
CN115075341A (en) | Regional water leakage pipeline detection method and system, storage medium and intelligent terminal | |
Rathi | A critical review of leakage detection strategies including pressure and water quality sensor placement in water distribution systems–sole and integrated approaches for leakage and contamination intrusion | |
Pandian et al. | Using computational fluid dynamics and deep learning for leak detection and localization in a smart water management system | |
Zhang et al. | Burst detection based on multi-time monitoring data from multiple pressure sensors in district metering areas | |
Takahashi et al. | Water distribution system model calibration under uncertainty environments | |
CN112840362A (en) | Computer-implemented method for modeling hydraulic networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |