Water supply network leakage location identification method and storage medium

Abstract

The present invention discloses a method for identifying leakage locations in a water supply network and a storage medium. The method comprises real-time monitoring of pressure and flow data at different locations in the water supply network; using lattice Boltzmann parallel methods to solve the right-hand residual term of a pressure-flow relationship equation as a real-time monitoring residual term; and determining the leakage location of the water supply network corresponding to the real-time monitoring residual term according to a pre-established mapping relationship between the leakage location and the residual term. The present invention significantly improves the efficiency of solving the equation by using lattice Boltzmann parallel methods, and at the same time avoids the complex process of feature extraction by establishing a mapping relationship between the leakage location and the residual term, thereby achieving rapid and accurate identification of the leakage location in the network.

Description

Technical Field

[0001] The present invention relates to the technical field of fire water supply network fault monitoring and diagnosis, and in particular to a method for identifying a leakage position of a water supply network and a storage medium. Background Art

[0002] The water supply network plays a crucial role in the firefighting water supply system. First, it is a vital component of the firefighting system, responsible for transporting water to various firefighting points, providing the necessary water for extinguishing fires. Second, the water supply network must maintain sufficiently high pressure to support firefighting efforts at the scene. If a water supply network leaks, sufficient water will not be available to extinguish the fire in a timely manner, causing the fire to spread rapidly, resulting in property damage and casualties. Therefore, quickly and accurately identifying the location of water supply network leaks is crucial.

[0003] Flow and pressure sensors are often deployed in water supply networks to monitor their operational status. The distribution of pressure and flow data at different locations in the network varies significantly between normal operating conditions (no leaks) and leaking conditions, as well as for different leak locations. Therefore, changes in water supply pressure and flow are the primary means of detecting leaks in water supply networks.

[0004] In related technologies, Feng Kai's master's thesis, "Research on Leak Detection and Location Methods for Ship Pipeline Networks," describes a real-time modeling approach: a precise real-time model of the pipeline is constructed using computer calculations and executed synchronously with the actual pipeline. By periodically comparing pressure and flow rates with those calculated by the theoretical model, the leak location is determined. This method relies heavily on the accuracy of the initial model and software, and its leak location mechanism is still based on the pressure gradient method, which has significant limitations.

[0005] Therefore, how to quickly and accurately identify the leakage location of the water supply network by monitoring the changes in pressure and flow data at different locations in the water supply network is an urgent problem that needs to be solved. Summary of the Invention

[0006] The technical problem to be solved by the present invention is how to quickly and accurately identify the leakage location of a water supply network.

[0007] The present invention solves the above technical problems through the following technical means:

[0008] The present invention proposes a method for identifying a leakage location in a water supply network, the method comprising:

[0009] Real-time monitoring of pressure and flow data at different locations in the water supply network;

[0010] The residual term on the right side of the pressure-flow relationship equation is solved in parallel using the lattice Boltzmann method as the real-time monitoring residual term;

[0011] According to the pre-established mapping relationship between the leakage position and the residual term, the leakage position of the water supply network corresponding to the real-time monitoring residual term is determined.

[0012] Furthermore, the process of constructing the mapping relationship between the leakage position and the residual term includes:

[0013] Obtain pressure and flow data at different locations in the water supply network under normal operating conditions, and use the lattice Boltzmann parallel method to solve the residual term on the right side of the pressure-flow relationship equation under normal operating conditions as the residual term under normal operating conditions;

[0014] Set a set number of leakage locations and collect pressure and flow data at different locations of the water supply network under different leakage locations;

[0015] For each leakage location, the lattice Boltzmann method is used to solve the residual term on the right side of the pressure-flow relationship equation under the leakage condition in parallel as the residual term under the leakage condition;

[0016] Based on all leakage locations and the residual items under the corresponding leakage conditions, a mapping relationship between different leakage locations and residual items is established to form a leakage location-residual item mapping relationship comparison table, including leakage locations, residual items and pressure and flow data at different locations of the water supply network under the working condition.

[0017] Furthermore, the leakage positions are arranged evenly or unevenly.

[0018] Furthermore, the residual term on the right side of the pressure-flow relationship equation solved in parallel using the lattice Boltzmann method is

[0019]

[0020] Where δ is the residual term on the right side, φ is the pressure, and f is the fluid velocity. f=∑ i c i h i / ∑ i h i , c i is the model discrete velocity, h i is the discrete distribution function, i is the discrete velocity index, is the gradient operator.

[0021] Furthermore, the pressure-flow relationship equation is:

[0022]

[0023] The pressure-flow relationship equation is solved in parallel using the lattice Boltzmann method, including the collision step, migration step, and statistical pressure-flow data. The collision step formula is expressed as:

[0024]

[0025] The migration step formula is expressed as:

[0026]

[0027] The statistical pressure data and flow data are:

[0028]

[0029] Where, is the discrete distribution function after collision, h i (x, t) is the discrete distribution function at position x and time t, Ω i (h) is the collision operator, x is the position index, and t is the time index.

[0030] Furthermore, after determining the water supply network leakage position corresponding to the real-time monitoring residual term according to the pre-established mapping relationship between the leakage position and the residual term, the method further includes:

[0031] Calculate the non-local cosine similarity between the pressure and flow data of different locations in the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions;

[0032] When the non-local cosine similarity is greater than a set similarity threshold, it is determined that the identified pipe network leakage location is correct;

[0033] When the non-local cosine similarity is less than or equal to a set similarity threshold, it is determined that the identified pipe network leakage location is wrong.

[0034] Furthermore, the calculation formula of the non-local cosine similarity is:

[0035]

[0036] Where cosα Φ is the non-local cosine similarity of pressure data, Φ is the pressure vector monitored at the current moment, Φ 0 is the pressure vector under historical leakage conditions.

[0037] Similarly, the calculation formula for the non-local cosine similarity of traffic data is:

[0038]

[0039] Where cosα F is the non-local cosine similarity of traffic data, F is the traffic data monitored at the current moment, and F 0 This is the flow rate data under historical leakage conditions.

[0040] Furthermore, the calculation of the non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions includes:

[0041] According to the vector length n, the first non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions is calculated;

[0042] Determining whether the first non-local cosine similarity is within a set range;

[0043] If yes, then comparing the non-local cosine similarity with the similarity threshold to determine whether the pipeline network leakage location identification result is accurate;

[0044] Otherwise, the vector length is reset to 2n, and the second non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions is calculated.

[0045] Furthermore, after calculating the second non-local cosine similarity, the method further includes:

[0046] The difference between the second non-local cosine similarity and the first non-local cosine similarity is calculated, and the first non-local cosine similarity is output when the difference is within a set range; otherwise, the second non-local cosine similarity is output for comparison with a similarity threshold.

[0047] In addition, the present invention also proposes a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the method for identifying a leakage location in a water supply network as described above is implemented.

[0048] The advantages of the present invention are:

[0049] (1) The present invention preliminarily uses the lattice Boltzmann parallel method to solve the pressure-flow relationship equation to obtain the mapping relationship between the leakage position and the residual term of the equation. When monitoring the leakage position of the water supply network, the pressure and flow data of different positions of the water supply network are monitored in real time. The residual term on the right side of the pressure-flow relationship equation is solved in parallel by the lattice Boltzmann method as the real-time monitoring residual term. The pre-constructed mapping relationship is used to identify the leakage position of the network corresponding to the real-time monitoring residual term. By using the lattice Boltzmann parallel method to solve the pressure and flow relationship equation, the efficiency of the equation solution is significantly improved. At the same time, by establishing the mapping relationship between the leakage position and the residual term, the complex process of feature extraction is avoided, and the leakage position of the network is quickly and accurately identified.

[0050] (2) Based on the established mapping relationship between the leakage location and the residual term, the accuracy of the leakage location of the water supply network is further determined by comparing the pressure and flow data using non-local cosine similarity. This reduces the number of data comparison samples and improves the accuracy of calculating the leakage location of the water supply network.

[0051] Additional aspects and advantages of the present invention will be set forth in part in the description which follows and, in part, will be obvious from the description which follows, or may be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] Figure 1 This is a flow chart of a method for identifying a leak location in a water supply network according to an embodiment of the present invention;

[0053] Figure 2 Schematic diagram of the arrangement of pressure and flow sensors in one embodiment of the present invention;

[0054] Figure 3 This is a graph showing the changing trend of pressure data at different leakage locations in one embodiment of the present invention;

[0055] Figure 4 This is a graph showing the changing trend of flow data at different leakage locations in one embodiment of the present invention. DETAILED DESCRIPTION

[0056] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0057] like Figure 1 As shown, an embodiment of the present invention provides a method for identifying a leakage location in a water supply network, the method comprising the following steps:

[0058] S10, real-time monitoring of pressure and flow data at different locations of the water supply network;

[0059] It should be noted that, in this embodiment, pressure data and flow data at different locations of the pipeline network can be measured by means of pressure sensors and flow sensors arranged at different locations of the pipeline network.

[0060] S20, using the lattice Boltzmann parallel method to solve the right-hand side residual term of the pressure-flow relationship equation as the real-time monitoring residual term;

[0061] S30: Determine the water supply network leakage position corresponding to the real-time monitoring residual item according to the pre-established mapping relationship between the leakage position and the residual item.

[0062] It should be noted that the pressure / flow coupling relationship equation is a nonlinear partial differential equation. If the traditional method is used to solve it, the differential equation needs to be discretized and then the large-scale linear equation group after discretization needs to be solved iteratively. The amount of calculation is huge and the solution time is long. As a mesoscopic method, Lattice Boltzmann does not need to discretize macroscopic equations. It performs collision and migration through the distribution function of "fluid particles" to restore the macroscopic flow behavior of the fluid. Since it does not involve solving complex linear equations, LBM is easy to implement, and local operations are very suitable for large-scale parallel computing. Using Lattice Boltzmann to solve the pressure / flow relationship equation in parallel significantly improves the efficiency of equation solution. At the same time, by establishing a mapping relationship between the leakage location and the residual term, the complex process of feature extraction is avoided, and the leakage location of the pipeline network can be quickly and accurately identified.

[0063] As a further preferred technical solution, in step S30, the process of constructing the mapping relationship between the leakage position and the residual term includes the following steps:

[0064] Obtain pressure and flow data at different locations in the water supply network under normal operating conditions, and use the lattice Boltzmann parallel method to solve the residual term on the right side of the pressure-flow relationship equation under normal operating conditions as the residual term under normal operating conditions;

[0065] Set a set number of leakage locations and collect pressure and flow data at different locations of the water supply network under different leakage locations;

[0066] For each leakage location, the lattice Boltzmann method is used to solve the residual term on the right side of the pressure-flow relationship equation under the leakage condition in parallel as the residual term under the leakage condition;

[0067] Based on all leakage locations and the residual items under the corresponding leakage conditions, a mapping relationship between different leakage locations and residual items is established to form a leakage location-residual item mapping relationship comparison table, including leakage locations, residual items and pressure and flow data at different locations of the water supply network under the working condition.

[0068] Specifically, this embodiment sets up a water supply network including 2 water tanks, 2 water pumps, 6 fire hydrants, 10 valves and 26 pipes. First, the pressure / flow data at different locations of the water supply network under normal working conditions are obtained. Then, several groups of leakage locations are set as needed. The setting of leakage locations is very flexible and can be evenly distributed or unevenly distributed. The pressure / flow data at different locations of the water supply network under different leakage locations are collected. Only one group of leakage locations is set for each test. The layout positions of pressure / flow sensors are as follows: Figure 2 As shown, the pressure data change trends and flow data change trends at different leakage locations are shown in Figure 3 and Figure 4 .

[0069] In this example, the pipeline is 400 meters long, with 20 sets of leak locations evenly spaced. In the first test, leak hole 1 was opened, and the remaining leak holes were closed. Pressure and flow data were collected at different locations in the water supply network under different leak locations and stored in a database, as shown in Table 1.

[0070] Table 1 stores pressure / flow data

[0071]

[0072] In the calculation process, truncation error, rounding error, etc. are inevitably introduced. In order to correct the right-hand residual term under the leakage condition, the pressure and flow data at different locations of the water supply network under normal conditions are obtained, and the lattice Boltzmann method is used to solve the right-hand side residual term of the pressure-flow relationship equation under normal conditions in parallel as the residual term under normal conditions. This residual term is used to correct the residual term under the leakage condition. Therefore, for each group of leakage locations, the lattice Boltzmann method is used to solve the pressure / flow relationship equation in parallel, and the residual term on the right side of the equation is calculated. The mapping relationship between different leakage locations and residual terms is established, as shown in Table 2.

[0073] Table 2 Leakage points and residual items

[0074]

[0075] The pressure / flow data at different locations of the water supply network are monitored in real time by arranged pressure / flow sensors. The pressure / flow relationship equation is solved in parallel using the lattice Boltzmann method to calculate the residual term on the right side of the pressure / flow relationship equation at the current moment, as shown in Table 3.

[0076] Table 3 Current pressure / flow and right-hand residual

[0077]

[0078] According to the mapping relationship between leakage location and residual term, different leakage locations lead to different pressure / flow distributions. Therefore, the residual term of the pressure / flow equation is generally different. The leakage location can be preliminarily determined by comparing the value of the residual term with the residual terms of different leakage locations.

[0079] Furthermore, water is an incompressible fluid, and the pressure / flow of the water supply network satisfies the following equation:

[0080]

[0081] Where φ is the pressure (MPa) and f is the fluid velocity (m / s). The flow rate can be obtained by integrating the fluid velocity along the cross-sectional area. The pressure / flow relationship equation is solved in parallel using the lattice Boltzmann method. The steps are as follows:

[0082] (1) Collision step:

[0083]

[0084] (2) Migration step:

[0085]

[0086] (3) Statistical pressure / flow:

[0087]

[0088] It should be noted that using the lattice Boltzmann method to solve the pressure / flow relationship equation is a local operation and is very suitable for parallel computing.

[0089] Furthermore, due to the existence of measurement errors and calculation errors, the right-hand side of the pressure / flow relationship equation is generally not 0, and the residual term on the right-hand side of the equation needs to be calculated. The calculation formula is as follows:

[0090]

[0091] Where c i is the model discrete velocity, h i is the discrete distribution function, i is the discrete velocity index, is the gradient operator, is the discrete distribution function after collision, h i (x, t) is the discrete distribution function at position x and time t, Ω i (h) is the collision operator, x is the position index, and t is the time index.

[0092] The following is the partial calculation code for solving the pressure / flow relationship equation using the lattice Boltzmann method:

[0093] (1) Collision step:

[0094] void COLL(void){

[0095] int X;

[0096] int Y;

[0097] int I;

[0098] #pragma omp parallel for private(Y,I)num_threads(16)

[0099] for(X=0; X<=GRIX; X++){

[0100] for(Y=0;Y<=GRIY;Y++){

[0101] for (I = 0; I < 9; I++) {

[0102] FUN1[X][Y][I] = (1.0 - 1.0 / TIME) * FUN0[X][Y][I] + 1.0 / TIME * EQUI(I, DENS[X][Y], VELX[X][Y], VELY[X][Y]);

[0103] }

[0104] }

[0105] }

[0106] }

[0107] (2) Migration step:

[0108] void STRE(void) {

[0109] int X;

[0110] int Y;

[0111] int I;

[0112] #pragma omp parallel for private(Y, I) num_threads(16)<000028​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​void PHYS(void){

[0125] int X;

[0126] int Y;

[0127] int I;

[0128] #pragma omp parallel for private(Y,I)num_threads(16)

[0129] for(X=0; X<=GRIX; X++){

[0130] for(Y=0;Y<=GRIY;Y++){

[0131] DENS[X][Y]=0.0;

[0132] VELX[X][Y]=0.0;

[0133] VELY[X][Y]=0.0;

[0134] for(I=0;I<9;I++){

[0135] DENS[X][Y]+=FUN0[X][Y][I];

[0136] VELX[X][Y]+=INT1[I]*FUN0[X][Y][I];

[0137] VELY[X][Y]+=INT2[I]*FUN0[X][Y][I];

[0138] }

[0139] VELX[X][Y] / =DENS[X][Y];

[0140] VELY[X][Y] / =DENS[X][Y];

[0141] }

[0142] }

[0143] }

[0144] Furthermore, since different leakage locations correspond to the same or very similar residual terms, in order to avoid misjudgment of the leakage location, this embodiment locates the leakage location based on the mapping relationship between the leakage location and the residual term, and combines the non-local cosine similarity to compare the pressure / flow data, thereby further improving the accuracy of calculating the leakage location of the water supply network.

[0145] Accordingly, after determining the water supply network leakage position corresponding to the real-time monitoring residual item according to the pre-established mapping relationship between the leakage position and the residual item in step S30, the method further comprises the following steps:

[0146] S31, calculating the non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions;

[0147] S32: When the non-local cosine similarity is greater than a set similarity threshold, determining that the identified pipe network leakage location is correct;

[0148] S33: When the non-local cosine similarity is less than or equal to a set similarity threshold, determine that the identified pipe network leakage location is wrong.

[0149] Specifically, assuming that the real-time monitoring pressure vector at the current moment is:

[0150] Φ=(φ -n ,φ -n+1 …,φ -1 ,φ o ,φ1…,φ n-1 ,φ n )

[0151] Where φ0 is the pressure value at the node where the leakage location is preliminarily determined.

[0152] Assuming that the leak location is the preliminary judgment leak location, according to the leak test, the historical pressure vector under the leak condition is

[0153]

[0154] Calculate the non-local cosine similarity of pressure vectors:

[0155]

[0156] Where cosα Φ is the non-local cosine similarity of pressure data, Φ is the pressure vector monitored at the current moment, Φ 0 is the pressure vector under historical leakage conditions.

[0157] Similarly, the calculation formula for the non-local cosine similarity of traffic data is:

[0158]

[0159] Where cosα F is the non-local cosine similarity of traffic data, F is the traffic data monitored at the current moment, and F 0 This is the flow rate data under historical leakage conditions.

[0160] Furthermore, if and only if cosα Φ ≥0.9 and cosα F If the value is ≥0.9, the pipeline network leakage location is correctly identified. Otherwise, it is considered that the pressure and flow data under this working condition do not match, and the pipeline network leakage location is incorrectly identified. The leakage location should be re-selected for identification.

[0161] As a further preferred technical solution, step S31: calculating the non-local cosine similarity between the real-time monitoring pressure and flow data at different locations of the water supply network and the pressure and flow data under historical leakage conditions, specifically includes the following steps:

[0162] S311, calculating the first non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions, based on a vector length of n (N / 4≤n≤N / 2), where N is the number of pressure and flow data;

[0163] S312, determining whether the first non-local cosine similarity is within a set range, if so, executing step S313, otherwise executing step S314;

[0164] S313, comparing the non-local cosine similarity with the similarity threshold to determine whether the pipeline network leakage location identification result is accurate;

[0165] S314: Reset the vector length to 2n, and calculate the second non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions.

[0166] As a further preferred technical solution, after step S314, the method further includes:

[0167] The difference between the second non-local cosine similarity and the first non-local cosine similarity is calculated, and the first non-local cosine similarity is output when the difference is within a set range; otherwise, the second non-local cosine similarity is output for comparison with a similarity threshold.

[0168] It should be noted that the non-local cosine similarity calculation method proposed in this embodiment is different from the traditional method. The traditional method only calculates the cosine similarity once and then determines the similarity between the two vectors. The innovation of this method is that the determination of the vector length n is very flexible. At the beginning, a smaller n is selected and the cosα of the two vectors is calculated. If cosα is too high or too low, it means that the vector length may be too short. The vector length is reset to 2n and the cosα of the two vectors is calculated again. If the cosα calculated twice does not change much, the cosα calculated for the first time is output, and it is considered that the vector length has no significant effect on the similarity, and the calculation is stopped. Otherwise, the cosα calculated for the second time shall prevail.

[0169] In addition, another embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the method for identifying a leakage location in a water supply network as described in the above embodiment is implemented.

[0170] It should be noted that other embodiments of the computer-readable storage medium of the present invention or implementation methods thereof can refer to the above-mentioned method embodiments and will not be described in detail here.

[0171] It should be noted that the logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, apparatus, or device (such as a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device), or in conjunction with such instruction execution system, apparatus, or device. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transmit a program for use by an instruction execution system, apparatus, or device, or in conjunction with such instruction execution system, apparatus, or device. More specific examples (non-exhaustive list) of computer-readable media include the following: an electrical connection portion having one or more wires (electronic device), a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disc read-only memory (CDROM). Furthermore, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing it in another suitable manner if necessary, and then storing it in a computer memory.

[0172] It should be understood that various parts of the present invention can be implemented using hardware, software, firmware, or a combination thereof. In the above-described embodiments, multiple steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0173] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in conjunction with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0174] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of the technical features being referred to. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.

[0175] Although the embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and are not to be construed as limitations on the present invention. A person skilled in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for identifying a leak location in a water supply network, characterized in that: The method comprises: Real-time monitoring of pressure and flow data at different locations in the water supply network; The residual term on the right side of the pressure-flow relationship equation is solved in parallel using the lattice Boltzmann method as the real-time monitoring residual term. The pressure-flow relationship equation is: Where, For pressure, is the fluid velocity, is the gradient operator; According to the pre-established mapping relationship between the leakage position and the residual term, the leakage position of the water supply network corresponding to the real-time monitoring residual term is determined.

2. The method for identifying a leak location in a water supply network according to claim 1, wherein: The process of constructing the mapping relationship between the leakage position and the residual term includes: Obtain pressure and flow data at different locations in the water supply network under normal operating conditions, and use the lattice Boltzmann parallel method to solve the residual term on the right side of the pressure-flow relationship equation under normal operating conditions as the residual term under normal operating conditions; Set a set number of leakage locations and collect pressure and flow data at different locations of the water supply network under different leakage locations; For each leakage location, the lattice Boltzmann method is used to solve the residual term on the right side of the pressure-flow relationship equation under the leakage condition in parallel as the residual term under the leakage condition; According to all leakage locations and the residual items under the corresponding leakage conditions, a mapping relationship between different leakage locations and residual items is established.

3. The method for identifying a leak location in a water supply network according to claim 2, wherein: The leakage positions are arranged evenly or unevenly.

4. The method for identifying a leak location in a water supply network according to claim 1, wherein: The residual term on the right side of the pressure-flow relationship equation solved in parallel using the lattice Boltzmann method is: Where, is the residual term on the right side.

5. The method for identifying a leak location in a water supply network according to claim 4, wherein: The method comprises: The pressure-flow relationship equation is solved in parallel using the lattice Boltzmann method, including the collision step, migration step, and statistical pressure-flow data. The collision step formula is expressed as: The migration step formula is expressed as: The statistical pressure data and flow data are: Where, is the discrete distribution function after collision, For location ,time The discrete distribution function of is the collision operator, is the position index, is the time index, is the model discrete velocity, is a discrete distribution function, is the discrete velocity index.

6. The method for identifying a leak location in a water supply network according to claim 1, wherein: After determining the water supply network leakage position corresponding to the real-time monitoring residual item according to the pre-established mapping relationship between the leakage position and the residual item, the method further includes: Calculate the non-local cosine similarity between the pressure and flow data of different locations in the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions; When the non-local cosine similarity is greater than a set similarity threshold, it is determined that the identified pipe network leakage location is correct; When the non-local cosine similarity is less than or equal to a set similarity threshold, it is determined that the identified pipe network leakage location is wrong.

7. The method for identifying a leak location in a water supply network according to claim 6, wherein: The calculation formula of the non-local cosine similarity is: Where, is the non-local cosine similarity of pressure data, is the pressure vector monitored at the current moment, is the pressure vector under historical leakage conditions; The calculation formula of non-local cosine similarity of traffic data is: Where, is the non-local cosine similarity of traffic data, The traffic data monitored at the current moment, This is the flow rate data under historical leakage conditions.

8. The method for identifying a leak location in a water supply network according to claim 6, wherein: The calculation of the non-local cosine similarity between the real-time monitoring pressure and flow data at different locations of the water supply network and the pressure and flow data under historical leakage conditions includes: According to the vector length , calculate the first non-local cosine similarity between the pressure and flow data of different locations of the real-time monitoring water supply network and the pressure and flow data under historical leakage conditions; Determining whether the first non-local cosine similarity is within a set range; If yes, then comparing the non-local cosine similarity with the similarity threshold to determine whether the pipeline network leakage location identification result is accurate; Otherwise reset the vector length to 2 , calculate the second non-local cosine similarity between the real-time monitoring pressure and flow data at different locations of the water supply network and the pressure and flow data under historical leakage conditions.

9. The method for identifying a leak location in a water supply network according to claim 8, wherein: After calculating the second non-local cosine similarity, the method further includes: The difference between the second non-local cosine similarity and the first non-local cosine similarity is calculated, and the first non-local cosine similarity is output when the difference is within a set range; otherwise, the second non-local cosine similarity is output for comparison with a similarity threshold.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for identifying a leakage location in a water supply network according to any one of claims 1 to 9 is implemented.

Images