Post-event analysis method for pipe explosion of water supply pipeline

Abstract

The invention discloses a post-analysis method for pipe explosion of a water supply pipeline. The method comprises the following steps: constructing an earth surface ponding model based on earth surface data by taking a pipe explosion point position as a center; inputting the pressure data and the flow data into an earth surface ponding model to obtain an earth surface ponding range and an earth surface ponding area; collecting a plurality of abnormal pipe points where the abnormal pressure data and the abnormal flow data are located, calculating a shortest path among the plurality of abnormal pipe points by using a graph theory algorithm, and obtaining a water cut-off range and a water cut-off area according to the shortest path; and constructing a pipe explosion influence grid and a pipe explosion influence index, obtaining pipe explosion influence degrees of a plurality of sub-grids in the pipe explosion influence grid according to the pipe explosion influence index, and performing rendering according to the pipe explosion influence degrees to obtain a visual analysis result. The method relates to the technical field of water supply pipeline emergency management, and solves the technical problems of inaccurate influence prediction, incomplete analysis indexes and low emergency response efficiency in an existing water supply pipeline explosion influence analysis method.

Description

Technical Field

[0001] The invention belongs to the field of urban emergency management, relates to fluid mechanics simulation technology, and specifically is a post-analysis method for water supply pipeline burst. Background Art

[0002] Due to the aging of urban water supply networks, over-service life, and lack of sound management, pipes often suffer structural damage, causing water leakage from the pipes to rise to the ground and burst. Pipe bursts not only directly affect users' water use, but may also trigger a series of secondary disasters, such as ground collapse, water pollution, and traffic safety issues. Therefore, accurately and quickly analyzing the impact range of water supply pipe bursts is crucial to formulating effective emergency response strategies and resource allocation plans.

[0003] In patent CN202211604720.X "A method for quickly analyzing the impact range of a pipe burst in a water supply pipeline", a method for quickly analyzing the impact range of a pipe burst in a water supply pipeline is provided, the method comprising: after determining the location of the pipe burst and the affected flow monitoring point, calculating the flow change rate of the flow monitoring point, and determining the initial impact range based on the pipeline type and the flow change rate, and then verifying the water supply impact on the initial range to determine the final impact range. This method predicts the impact range of a pipe burst in a water supply pipeline based on the acquired pipeline data, and only considers the impact of the monitoring point location information on the impact range of the pipe burst. Once a pipe burst occurs at a location far from the monitoring point, it is difficult to accurately reflect the actual impact of the pipe burst afterwards. In addition, this method renders a map based on the flow change rate of the water supply pipeline, lacks a dynamic simulation of the evolution of water accumulation, and cannot obtain the specific situation of subsequent changes in water accumulation. Its visualization effect is not conducive to decision makers quickly identifying areas that need priority treatment after a pipe burst. Summary of the invention

[0004] The present invention aims to solve at least one of the technical problems existing in the prior art; to this end, the present invention proposes a post-analysis method for a water supply pipe burst, which is used to solve the technical problems of inaccurate impact prediction, incomplete analysis indicators and low emergency response efficiency in the prior art.

[0005] To achieve the above object, the present invention provides a post-analysis method for a water supply pipeline burst, comprising:

[0006] S1, collecting the location of the pipe bursting point and the surface data, and constructing a surface water accumulation model based on the surface data and taking the location of the pipe bursting point as the center; wherein the surface data includes ground elevation data, soil type data, and soil permeability data;

[0007] S2, obtaining pressure data and flow data of the pipe burst point, inputting the pressure data and flow data into the surface water accumulation model, and obtaining the surface water accumulation range and surface water accumulation area;

[0008] S3, collecting several abnormal pipe points where abnormal pressure data and abnormal flow data are located, using a graph theory algorithm to calculate the shortest path between the several abnormal pipe points, and obtaining the water outage range and water outage area according to the shortest path;

[0009] S4, constructing a pipe burst influence grid and a pipe burst influence index, obtaining pipe burst influence degrees of several sub-grids in the pipe burst influence grid according to the pipe burst influence index, and rendering according to the pipe burst influence degrees to obtain a visual analysis result.

[0010] Based on the above technical process, by constructing a surface waterlogging model, the evolution process of road waterlogging after a pipe burst is dynamically simulated to obtain the prediction results of the waterlogging range and area, which will help to take drainage measures in a timely manner and ensure traffic safety. The graph theory algorithm is used to calculate the shortest path between the pipe burst point and other pipe points, which can quickly determine the impact range of water supply. By constructing a pipe burst impact grid and pipe burst impact indicators, the impact degree of the pipe burst is visualized, providing clearly classified impact analysis results, and improving the efficiency and effectiveness of emergency response.

[0011] Furthermore, the construction of a surface water accumulation model based on surface data and taking the pipe burst point as the center includes:

[0012] S1-1, construct the two-dimensional governing equations for the evolution of surface water based on the shallow water equation: Where E represents the conserved variable vector, F(E) represents the flux vector along the x direction, G(E) represents the flux vector along the y direction, Q(E) represents the source term vector, t represents the time of water accumulation evolution, h represents the water depth, g represents the gravitational acceleration, and q x and q y Respectively represent the single width flow in the x and y directions, and q x =hu,q y =hv, u and v are the velocity components in the x and y directions respectively, S fx and S fy They represent the slope of friction along the x and y directions, S 0x represents the bottom slope along the x direction, S 0y represents the bottom slope along the y direction, and z b Indicates the bottom elevation;

[0013] S1-2, using the local characteristics method, the shallow water equation is transformed into: The general two-dimensional control equation is obtained, where A and B represent general parameters;

[0014] S1-3, construct the Jacobian matrix of the flow, and solve the general two-dimensional control equation according to the Jacobian matrix to obtain: Among them, e A represents the eigenvalue of A, eB represents the eigenvalue of B, e represents, and f represents;

[0015] S1-4, using the difference method to solve the general two-dimensional control equation to obtain a solution result, and visualizing the solution result to obtain a surface water accumulation model.

[0016] Furthermore, the construction of the pipe burst impact grid includes:

[0017] S4-11, spatially superimpose the surface waterlogging range and water outage range to obtain the comprehensive impact range;

[0018] S4-12, using Graham scanning method to obtain the convex hull of the comprehensive influence range, and using the rotating caliper algorithm to obtain the minimum circumscribed rectangle of the convex hull;

[0019] S4-13, dividing the minimum circumscribed rectangle according to the preset length and width to obtain the pipe burst impact grid.

[0020] Furthermore, the pipe burst impact index includes human settlement density index, traffic water accumulation index and water outage water accumulation ratio;

[0021] The human settlement density index represents the ratio of the flow of people to the population density within the range of surface waterlogging during a preset time period.

[0022] The vehicle flow index represents the ratio of vehicle flow to surface waterlogging area in a preset time period within the range of surface waterlogging.

[0023] The water outage to water accumulation ratio represents the ratio of the predicted water outage area to the surface water accumulation area.

[0024] Furthermore, rendering according to the impact of the pipe explosion includes:

[0025] S4-21, collecting a number of pipe burst impact indicators of a number of sub-grids in the pipe burst impact grid in a preset time period;

[0026] S4-22, according to formula I i =0.4a i1 +0.3a i2 +0.3a i3 Get the impact of pipe burst of several sub-grids; where I i represents the impact of pipe explosion in the ith subgrid, a i1 represents the human settlement density index of the ith subgrid, a i2 represents the traffic water accumulation index of the ith subgrid, a i3 represents the water-stopping and water-accumulation ratio of the ith subgrid;

[0027] S4-23, dividing the pipe burst impact degrees of the plurality of subgrids according to a preset level threshold, to obtain a plurality of pipe burst impact degree levels; wherein the plurality of pipe burst impact degree levels include: severe impact, impact, slight impact, and no impact;

[0028] S4-24, coloring the plurality of subgrids according to the degree of influence of the pipe burst, and obtaining a visual analysis result of the influence of the pipe burst.

[0029] By comprehensively considering the scope of surface waterlogging and water outage, the pipe burst impact index and pipe burst impact grid are constructed, and the pipe burst impact degree of each sub-grid is calculated. The impact is graded and colored according to the impact degree to obtain a visual impact result. This helps the emergency management department to intuitively and quickly understand the specific scope and degree of impact of the pipe burst on residents' lives, transportation systems and water supply systems. Based on these quantitative analysis results, more scientific and reasonable emergency response strategies and resource allocation plans can be formulated, thereby effectively reducing the losses and risks caused by pipe bursts.

[0030] Compared with the prior art, the present invention has the following beneficial effects:

[0031] By constructing a surface water accumulation model and performing dynamic simulation based on the shallow water equation, the scope of road waterlogging and the area of ​​waterlogging after a pipe burst can be accurately predicted, which helps emergency management departments to take drainage measures in a timely manner and effectively ensure traffic safety. By using graph theory algorithms, the shortest path between the pipe burst point and other pipe points can be quickly calculated, and the affected range of water supply can be accurately determined, which improves the speed and efficiency of emergency response. By constructing a pipe burst impact grid and pipe burst impact indicators, the degree of impact of the pipe burst can be visualized and analyzed, providing clearly classified impact analysis results, helping emergency management departments to intuitively and quickly understand the specific scope and degree of the impact of the pipe burst on residents' lives, transportation systems and water supply systems. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0033] Figure 1 A schematic flow chart of a post-analysis method for a water supply pipe burst provided by the present invention;

[0034] Figure 2 This is a visualization effect diagram of the water accumulation evolution of the surface water accumulation model provided by the present invention. DETAILED DESCRIPTION

[0035] The technical solution of the present invention will be clearly and completely described below in conjunction with the embodiments. Obviously, the described embodiments are only 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 ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0036] See also Figure 1-Figure 2 The embodiment of the present invention provides a post-analysis method for a water supply pipeline burst, comprising:

[0037] S1, collecting the location of the pipe bursting point and the surface data, and constructing a surface water accumulation model based on the surface data with the location of the pipe bursting point as the center; wherein the surface data includes ground elevation data, soil type data, and soil permeability data;

[0038] S2, obtaining pressure data and flow data of the pipe burst point, inputting the pressure data and flow data into the surface water accumulation model, and obtaining the surface water accumulation range and surface water accumulation area;

[0039] S3, collecting several abnormal pipe points where abnormal pressure data and abnormal flow data are located, using a graph theory algorithm to calculate the shortest path between the several abnormal pipe points, and obtaining the water outage range and water outage area according to the shortest path;

[0040] S4, based on the scope of surface water accumulation and the scope of water outage, construct a pipe burst influence grid and a pipe burst influence index, obtain the pipe burst influence degree of several sub-grids in the pipe burst influence grid according to the pipe burst influence index, and visualize the pipe burst influence degrees of several sub-grids to obtain a visualized analysis result of the pipe burst influence.

[0041] In this embodiment, suppose a main water supply pipe in a city bursts near a crossroads, and there is an urgent need to quickly analyze the impact range of the burst. First, the ground elevation data, soil type data, and soil permeability data of the area should be collected to provide basic data for the subsequent construction of the ground surface water model. Then, the pressure data and flow data before and after the burst recorded by the pressure gauge and flow meter are obtained from the burst point and other pipe points, and several abnormal pipe points with abnormal pressure data and abnormal flow data are obtained;

[0042] If there is a pressure gauge and no flow meter in the burst detection, the burst point pressure is the maximum pressure difference in the pressure gauge device of the burst alarm, and the burst point flow can be calculated according to Q = (P / (ρgSL)) 0.5; where L represents the length of the pipeline, P represents the burst pressure difference of the pipeline, ρ represents the liquid density, g represents the acceleration of gravity, S represents the friction resistance of the pipeline, and the friction resistance per unit length S can be directly obtained by referring to the data in Table 1;

[0043]

[0044] Table 1

[0045] If there is a flow meter but no pressure gauge in the pipe burst detection system, the flow value at the pipe burst point should be the maximum value of the flow difference recorded by the flow meter device when the pipe burst alarm occurs. At this time, the pressure difference at the pipe burst point can be estimated by reverse deduction of the above formula for estimating flow. It should be noted that the flow difference refers to the part of the flow monitoring value that exceeds the preset dynamic threshold. When calculating, half an hour is selected as the calculation time, that is, the flow process line within half an hour is taken as the input condition for simulation. Similarly, the calculation of the pressure difference follows similar methods and principles.

[0046] Then, the two-dimensional control equation of surface water evolution is constructed based on the shallow water equation, and the conservation form of the two-dimensional control equation in the Cartesian coordinate system is obtained as follows: in E represents the conserved variable vector, F(E) represents the flux vector along the x direction, G(E) represents the flux vector along the y direction, Q(E) represents the source term vector, t represents the time of water accumulation evolution, h represents the water depth, g represents the gravitational acceleration, and q x and q y Respectively represent the single width flow in the x and y directions, and q x =hu,q y =hv, u and v are the velocity components in the x and y directions respectively, S fx and S fy They represent the slope of friction along the x and y directions, S 0x represents the bottom slope along the x direction, S 0y represents the bottom slope along the y direction, and z b Indicates the bottom elevation;

[0047] Then, we apply the local characteristics method to transform the shallow water equation into a more solvable form: The general two-dimensional governing equations are obtained, where A and B represent general parameters; then the Jacobian matrix of the flow is constructed: where is the wave speed, and c 2 =gh; In this embodiment, the eigenvalues ​​of the vector matrices A and B are: Substitute the Jacobian matrix for mathematical transformation to obtain The final form of the general two-dimensional control equation is obtained as: in,

[0048] Then, a custom script is used to divide the structured grid with the burst point as the center, and the size and shape of the grid are set according to the calculation requirements to ensure that the grid can cover the entire calculation domain; then, the numerical analysis tool is used to realize the differential approximation of the partial derivatives of the general two-dimensional control equation, and the partial differential equation is converted into a differential equation. Each node in the grid is iteratively calculated to obtain the physical quantity at each point and time in the calculation domain; then the calculation results are visualized to obtain a surface water accumulation model that can simulate the diffusion process of water accumulation on the surface; finally, the pressure data and flow data of the burst point collected above are input into the surface water accumulation model to obtain the surface water accumulation range and surface water accumulation area;

[0049] like Figure 2 As shown in the figure, based on the above process, the visualization result of the water accumulation evolution of the surface water accumulation model after a pipe burst is obtained. In the figure, x and y represent the spatial directions perpendicular to each other, and h+z represents the overall height from the reference plane to the water surface, in meters;

[0050] After a water supply pipe burst occurs, a rapid assessment of its impact on the water supply system in the surrounding area can help take timely countermeasures and minimize the inconvenience caused by the burst. The current image analysis of water supply pipe bursts relies on complex hydraulic models to estimate the pressure and flow changes after the burst, but the construction of hydraulic models is very time-consuming, and the data requirements are huge and the application scenarios are limited. In contrast, in this embodiment, a more efficient method is adopted: by collecting abnormal pressure data and abnormal flow data, a number of abnormal pipe points are identified, and then the graph theory algorithm is used to calculate the shortest path between these abnormal pipe points; based on the shortest path, the water outage range and water outage area can be quickly determined, so as to quickly obtain the impact range after the burst. This method is not only faster, but also avoids the cumbersome model construction process, providing strong support for emergency response.

[0051] In this embodiment, based on the location of the burst point, the water supply network is abstracted into a weighted graph G = (V, E, F), where V represents the node (pipe point), E represents the edge (pipe segment), and F represents the weight of the edge (such as pipe segment length, resistance, etc.). Then the Warshall-Floyd algorithm is used to solve the shortest path between any two nodes;

[0052] Specifically, first construct an n×n adjacency matrix A=(a ij ), where n represents the number of nodes, a ij represents the element in the i-th row and j-th column of the matrix A, that is, from node v i to v j The direct distance of ij =F(v i ,v j), F is a function that determines the direct distance, and can be set according to the length, resistance, flow capacity, etc. of the pipeline; and when v i With v j When there is no direct connection between ij =+∞(i≠j), indicating unreachable; when v i With v j If they are equal, then a ij = 0 (i = j), indicating the own distance; then initialize the distance matrix D = A and the path record matrix R = (r ij ), where D is used to record the shortest distance from one node to another in the graph, r ij It is used to record the next node number in the shortest path from one node to another in the graph, and r ij =j, indicating a direct connection, and setting the iteration counter k=1;

[0053] For all node pairs (i, j), check if there is a node that passes through the intermediate node v. k So that from v i to v j The path is shorter, that is, if d ik +d kj <d ij , then update d ij =d ik +d kj , and update r ij = k; and repeat this step until all possible paths have been checked;

[0054] When ij <0, it means there is a negative loop containing vertex vi, or when the number of iterations k is equal to the number of nodes n, the iteration is stopped; otherwise, let k = k + 1 and repeat the above path update process; finally, what is stored in the D matrix is ​​the shortest distance between any two points, and the R matrix provides the information of the intermediate nodes that constitute these shortest paths;

[0055] By analyzing the shortest path results, all pipe points affected by the pipe burst can be identified and marked in the geographic information system software, and then the scope and area of ​​water outage caused by the pipe burst can be obtained based on the relationship between these pipe points and specific communities.

[0056] Through the above process, not only can the impact of pipe burst on the water supply system be efficiently assessed, but also precise guidance can be provided for subsequent emergency repair work, effectively shortening the time to restore water supply.

[0057] In order to further accurately locate which areas are affected by the pipe burst incident and the degree of impact, in this embodiment, by constructing pipe burst impact indicators and pipe burst impact grids, the impact of the pipe burst incident on residents' lives and transportation systems can be quantitatively and qualitatively represented;

[0058] Specifically, firstly, the surface waterlogging range obtained by the surface waterlogging model and the water outage range obtained by the graph theory algorithm are spatially superimposed in the geographic information system software to obtain the comprehensive impact range; then the convex hull of the comprehensive impact range is obtained by the Graham scanning method, and the minimum circumscribed rectangle of the convex hull is obtained according to the rotating calotte algorithm; the minimum circumscribed rectangle is divided according to the preset length and width to obtain the pipe burst impact grid;

[0059] Then calculate the pipe burst impact index of each grid area and use formula I i =0.4a i1 +0.3a i2 +0.3a i3 Get the impact of pipe burst of several sub-grids; where I i represents the impact of pipe explosion in the ith subgrid, a i1 represents the human settlement density index of the ith subgrid, which is obtained by calculating the ratio of the flow of people in a preset time period within the range of surface water accumulation to the population settlement density; a i2 represents the traffic water accumulation index of the ith subgrid, which is obtained by calculating the ratio of the traffic volume to the surface water accumulation area in the preset time period within the water accumulation range; a i3 represents the water-out-water accumulation ratio of the ith subgrid, which is obtained by calculating the ratio of the water-out area to the surface water accumulation area; then the pipe burst impact of several subgrids is divided according to the preset level threshold, and several grids marked as severe impact, impact, slight impact, and no impact are obtained;

[0060] Assuming that there are 9 subgrids, the value of the pipe burst impact index of each subgrid is calculated, and the preset level threshold is: no impact: [0, 0.2]; slight impact: (0.2, 0.4]; impact: (0.4, 0.6]; severe impact: (0.6, 1], then the pipe burst impact degree and impact level of each subgrid are shown in Table 2:

[0061]

[0062] Table 2

[0063] Finally, the subgrids are graded and colored according to the impact level to obtain the visual analysis results of the pipe burst impact.

[0064] Part of the data in the above formula is calculated by removing the dimension and taking its numerical value. The formula is a formula closest to the actual situation obtained by software simulation of a large amount of collected data; the preset parameters and preset thresholds in the formula are set by technical personnel in this field according to actual conditions or obtained through simulation of a large amount of data.

[0065] Working principle of the present invention:

[0066] Based on the location information of the pipe burst point and related surface data, a surface water accumulation model is constructed, and the scope and area of ​​the surface water accumulation are calculated; then, the abnormal pressure and flow data of the mobile phone are used to identify the abnormal pipe points affected by the pipe burst, and the graph theory algorithm is used to calculate the shortest path between the abnormal pipe points, and the water outage range and area are determined; finally, the pipe burst impact grid is constructed, and the impact indicators of the pipe burst on residents' lives and transportation systems are comprehensively considered, and the sub-grids affected by the pipe burst are visualized using a graded coloring method to provide accurate guidance for subsequent emergency repairs and response measures.

[0067] The above embodiments are only used to illustrate the technical method of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical method of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical method of the present invention.

Claims

1. A post-analysis method for a water supply pipe burst, characterized in that: include: S1, collecting the location of the pipe bursting point and the surface data, and constructing a surface water accumulation model based on the surface data and taking the location of the pipe bursting point as the center; wherein the surface data includes ground elevation data, soil type data, and soil permeability data; S2, obtaining pressure data and flow data of the pipe burst point, inputting the pressure data and flow data into the surface water accumulation model, and obtaining the surface water accumulation range and surface water accumulation area; S3, collecting several abnormal pipe points where abnormal pressure data and abnormal flow data are located, using a graph theory algorithm to calculate the shortest path between the several abnormal pipe points, and obtaining the water outage range and water outage area according to the shortest path; S4, constructing a pipe burst influence grid and a pipe burst influence index, obtaining pipe burst influence degrees of several sub-grids in the pipe burst influence grid according to the pipe burst influence index, and rendering according to the pipe burst influence degrees to obtain a visual analysis result.

2. A post-analysis method for water supply pipe burst according to claim 1, characterized in that: The method of constructing a surface water accumulation model based on surface data and taking the pipe burst point as the center includes: S1-1, construct the two-dimensional governing equations for the evolution of surface water based on the shallow water equation: Where E represents the conserved variable vector, F(E) represents the flux vector along the x direction, G(E) represents the flux vector along the y direction, Q(E) represents the source term vector, t represents the time of water accumulation evolution, h represents the water depth, g represents the gravitational acceleration, and q x and q y Respectively represent the single width flow in the x and y directions, and q x =hu,q y =hv, u and v are the velocity components in the x and y directions respectively, S fx and S fy They represent the slope of friction along the x and y directions, S 0x represents the bottom slope along the x direction, S 0y represents the bottom slope along the y direction, and z b Indicates the bottom elevation; S1-2, using the local characteristics method, the shallow water equation is transformed into: The general two-dimensional control equation is obtained, where A and B represent general parameters; S1-3, construct the Jacobian matrix of the flow, and solve the general two-dimensional control equation according to the Jacobian matrix to obtain: Among them, e A represents the eigenvalue of A, e B represents the eigenvalue of B, e represents, and f represents; S1-4, using the difference method to solve the general two-dimensional control equation to obtain a solution result, and visualizing the solution result to obtain a surface water accumulation model.

3. A post-analysis method for water supply pipe burst according to claim 1, characterized in that: The construction of the pipe burst influence grid includes: S4-11, spatially superimpose the surface waterlogging range and water outage range to obtain the comprehensive impact range; S4-12, using Graham scanning method to obtain the convex hull of the comprehensive influence range, and using the rotating caliper algorithm to obtain the minimum circumscribed rectangle of the convex hull; S4-13, dividing the minimum circumscribed rectangle according to the preset length and width to obtain the pipe burst impact grid.

4. A post-analysis method for water supply pipe burst according to claim 1, characterized in that: The impact indicators of pipe burst include human settlement density index, traffic water accumulation index and water outage and water accumulation ratio; among them, The human settlement density index represents the ratio of the flow of people to the population density within the range of surface waterlogging during a preset time period. The vehicle flow index represents the ratio of vehicle flow to surface waterlogging area in a preset time period within the range of surface waterlogging. The water outage to water accumulation ratio represents the ratio of the predicted water outage area to the surface water accumulation area.

5. A post-analysis method for water supply pipe burst according to claim 1, characterized in that: The rendering according to the influence of the pipe explosion includes: S4-21, collecting a number of pipe burst impact indicators of a number of sub-grids in the pipe burst impact grid in a preset time period; S4-22, according to formula I i =0.4a i1 +0.3a i2 +0.3a i3 Get the impact of pipe burst of several sub-grids; where I i represents the impact of pipe explosion in the ith subgrid, a i1 represents the human settlement density index of the ith subgrid, a i2 represents the traffic water accumulation index of the ith subgrid, a i3 represents the water-stopping and water-accumulation ratio of the ith subgrid; S4-23, dividing the pipe burst influence degrees of the plurality of subgrids according to a preset level threshold, and obtaining a plurality of pipe burst influence degree levels; S4-24, coloring the plurality of subgrids according to the degree of influence of the pipe burst, and obtaining a visual analysis result of the influence of the pipe burst.

Images