Next Article in Journal
Cultural and Natural Resources in Tourism Island: Bibliometric Mapping
Next Article in Special Issue
A General Overview of the Risk-Reduction Strategies for Floods and Droughts
Previous Article in Journal
Evolution of Carbon Shadow Prices in China’s Industrial Sector during 2003–2017: A By-Production Approach
Previous Article in Special Issue
Precipitation Forecasting Using Multilayer Neural Network and Support Vector Machine Optimization Based on Flow Regime Algorithm Taking into Account Uncertainties of Soft Computing Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Sensor Placement for Disaster Prevention in Water Distribution Networks for Important Users

1
College of Architecture and Civil Engineering, Beijing University of Technology, Beijing 100124, China
2
Institute of Earthquake Resistance and Disaster Reduction, Beijing University of Technology, Beijing 100124, China
3
College of Architecture and Urban Planning, Beijing University of Technology, Beijing 100124, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(2), 723; https://0-doi-org.brum.beds.ac.uk/10.3390/su12020723
Submission received: 2 December 2019 / Revised: 10 January 2020 / Accepted: 13 January 2020 / Published: 19 January 2020
(This article belongs to the Special Issue Sustainable Water Resources Management for Disaster Risk Reduction)

Abstract

:
Sensor placement for disaster prevention for important users in urban water distribution networks is essential for post-earthquake monitoring and repair. Herein, we proposed a sensor placement approach for disaster prevention monitoring for important users, to (a) improve the fault diagnosis ability of the water distribution network and to (b) guarantee the function of emergency services for key nodes after an earthquake. First, an evaluation system of node users’ disaster prevention impact factors was presented to evaluate the node influence degree from three aspects: post-earthquake leakage, emergency support and topology structure; and the weight values of node users’ disaster prevention impact factors were obtained. Second, a post-earthquake hydraulic analysis model based on the pressure-driven demand was used to calculate the water shortage ratio of nodes. Third, using the three-way clustering integration method, the results of four clustering techniques were integrated to divide the monitoring domain in the water distribution network based on sensitivity analysis. Finally, on basis of the sensitivity matrix, the division of the monitoring area and the impact factors of node users’ disaster prevention were combined to place sensors for post-earthquake disaster prevention in the water distribution network. Detailed computational experiments for a real urban water network in China were performed and compared with the results of other traditional techniques to evaluate the performance of the proposed approach. The results show that the approach is better than traditional methods. It not only considers the actual hydraulic information of the water distribution network, but also the important user nodes after an earthquake, and is of great significance for emergency command and rescue and disaster relief after an earthquake in the city.

1. Introduction

As a key emergency infrastructure in the event of an earthquake, the urban water supply system is an important basic condition for earthquake relief and life safety guarantee [1]. As a part of the lifeline system, a water distribution network (WDN) not only plays a momentous role in the normal operation of the city, but is also particularly important for water use in earthquake shelters and the rescue of secondary disasters. However, WDNs often suffer from varying degrees of earthquake damage, resulting in loss or weakening of functions, which cannot meet the post-earthquake emergency rescue service function. If the key user nodes cannot be guaranteed to use water after an earthquake, this will have a negative impact on the rescue. Therefore, it is crucial to place sensors for disaster prevention considering important users for strengthening the efficiency of post-earthquake monitoring of WDNs, ensuring the water demand of important user nodes, and improving the emergency repair ability of the lifeline system after a disaster.
The protection and monitoring of important nodes or pipelines are important aspects of improving the resilience of the lifeline network. Jie Li et al. [2] introduced an element investment importance analysis in the seismic topology optimization research of the lifeline network, and analyzed the nodes and lines of high importance in the lifeline network in order to speed up the convergence of optimization processes. Considering node importance and power generation recovery capability, Can Zhang et al. [3] proposed a two-stage grid reconfiguration strategy, which involves a node importance evaluation method based on the concept of regret. Wenli Fan et al. [4] put forward a method for evaluating the importance of multi-attribute nodes in complex grids based on the Gini coefficient. They introduced the Gini coefficient, used in economics, to the importance evaluation of multi-attribute nodes as well as to determine the weights of indicators. Considering the node importance, Fei Wang et al. [5] proposed a classification of emergency repair scope for urban WDNs, and grouped them from three aspects using graph theory: node importance, function importance and structure importance. Anne-Marie Kermarrec et al. [6] introduced a novel form of centrality: the second order centrality, which can be used to accurately identify critical nodes as well as to globally characterize graphs’ topology in a distributed way. Through the application of the neighborhood rule, Ahmad Zareie et al. [7] proposed an index to determine node centrality using the notions of Shannon entropy and Jensen–Shannnon divergence for ranking influential nodes. Michalis Fragiadakis et al. [8] presented a methodology for the seismic assessment of the reliability of urban WDNs based on general seismic assessment standards, which incorporates data of past non-seismic damage, the vulnerabilities of the network components against seismic loading, and the topology of an urban WDN. The above research works evaluated the importance of network nodes or lines in terms of network functions, topological structure and post-disaster reliability. If the importance assessment is introduced into the monitoring sensor placement study, the priority of protection and monitoring of important nodes or pipelines in the network analysis can be improved.
The research studies on sensor placement are currently concentrated in the field of health monitoring. It is the focus of the problem to study the optimal placement of sensor from the aspects of the monitoring domain and sensitivity based on monitoring parameters. The existing research works mostly use optimization algorithms and clustering algorithms to solve these two requirements, such as sensitivity analysis, neural network algorithms, genetic algorithms [9] and fuzzy clustering algorithms. Different clustering algorithms actually reflect the different levels of clustering parameters and different degrees of consideration. Lingli Chen [10] used the effective independence method (EfI), the Fisher information matrix (FIM) maximization criterion, and the correlation coefficients method to achieve optimal sensor placement. It was concluded that the calculation process of the node correlation coefficient method is relatively simple and stable, which increases the consideration of the parameters’ difference in essence. D. Kowalski et al. [11] proposed a fractal geometry method for water quality and pressure monitoring systems in WDNs, which incorporates the impact of network structure and user characteristics into the analysis. Jonas Kjeld Kirstein et al. [12] used topological clustering as a tool for water supply utilities for monitoring preparation and contamination contingency plans. Using a pressure sensitivity matrix and an exhaustive search strategy, Fatiha Nejjarih et al. [13] presented an optimal sensor placement strategy. In their work, the evidence C-means clustering method was applied to reduce the scale and complexity of the placement problem, and the effectiveness of this strategy was demonstrated by an example of a WDN in Barcelona. Adria Soldevila et al. [14] proposed an approach for sensor placement in WDNs by using a classifier to locate leakage, which combines correlation and redundancy with a genetic algorithm by using a hybrid feature selection algorithm. These studies show the degree to which different methods consider the feature factors. In the research of disaster prevention monitoring layout, the introduction of the feature of important nodes that need to be protected and monitored after an earthquake can be used as the basic research idea of this paper.
In terms of post-earthquake hydraulic calculation, the University of Cornell has developed GIRAFFE software, which assumes that there are only two cases of node water demand after an earthquake in a WDN: fully supplied or completely not supplied [15,16,17,18]. In fact, WDN is in a low-pressure water supply state after an earthquake, and the nodes are partially supplied rather than fully supplied [19,20]. Zhao Han et al. [21] presented a post-earthquake hydraulic analysis model of WDN based on the pressure-driven demand (PDD), which overcomes the above shortcomings.
In addition, most clustering algorithms are unstable, and mostly depend on parameters and initialization to a great extent, which may produce different results for the same dataset [22,23]. Therefore, the clustering integration method has been proposed and widely considered, and it shows higher robustness and stability than single clustering algorithms. LingChao Hu [24] proposed a voting-based three-way decision clustering integration method, which is more robust than a single clustering method.
In this paper, we proposed a sensor placement strategy for disaster prevention for important users. Based on the hydraulic analysis, an evaluation system of node users’ disaster prevention impact factors was established, and the impact factor weight values of node users were calculated from the three aspects of post-earthquake water shortage, emergency support and topological structure. The post-earthquake hydraulic analysis model based on PDD was used to calculate the water shortage ratio of nodes, which was combined with the consideration of the importance features as the placement optimization parameter. On this basis, using the three-way decision cluster integration method and integrating four clustering methods to divide the monitoring domain, the decision-making value of different clustering methods can be reflected comprehensively. The structure of this paper is displayed in Figure 1.

2. Disaster Prevention Impact Factors of Node Users

Sensor placement for health monitoring considers the sensor’s optimal placement during daily leakage. After an earthquake, the key node user of the WDN plays a decisive role in disaster relief and water use for victims. Therefore, for protecting people’s lives, it is necessary to evaluate the important nodes of WDNs before earthquakes. Thus, we proposed an evaluation system of disaster prevention impact factors for important node users of WDNs. The nodal leakage can reflect the strength of the nodal water transport function after an earthquake. In addition, considering that the locations of the nodes belong to different levels of disaster prevention emergency support objects, such as municipal emergency command centers, fire control centers, special service fire stations, and central evacuation sites, we thought that the minimum guarantee levels of these emergency protection objects can directly reflect the important position and influence degree of the nodes in a WDN after an earthquake. In terms of topology, the betweenness of a node indicates the importance of the node in the network. Therefore, from three aspects: water shortage, earthquake emergency prevention and topology, the evaluation system for node users’ impact factors in a WDN was constructed (Figure 2).

2.1. Node Water Shortage Calculation Based on PDD

The PDD post-earthquake hydraulic analysis model [21] assumes that the nodal water demand is supplied when the node water pressure reaches or exceeds the design water pressure. The node water supply is 0 when the node water pressure is lower than or equal to the minimum water pressure. The node water demand and node water pressure should meet a certain functional relationship, when the node water pressure is between the design water pressure and the minimum water pressure. The relationship between nodal water pressure and nodal water demand proposed by Wagner et al. [25] was adopted, as shown in Equation (1):
Q i = { 0 , H i < H i min Q i H i H i min H i d e s H i min , H i min H i < H i d e s Q i , H i d e s H i
where Qi* is the water supply of node i considering the relationship between node water pressure (m) and water demand (L/s), Qi is the initial water demand of node i (L/s), Hi is the calculated water pressure of node i (m), Himin is the minimum node water pressure required for node i water supply greater than 0 (m), and Hides is the design water pressure required for node i to meet the initial water demand Qi (m).
To determine the number and location of the damage points, we utilized Monte Carlo simulation to generate a series of random numbers of pipeline seismic damage that obey uniform distribution along the length of the pipeline. For each damage point, two damage states (break and leakage) were generated by a Monte Carlo simulation random process again to form the post-earthquake hydraulic model. Finally, epanet2.dll was used to calculate the hydraulic adjustment of WDN based on PDD. The calculations of the average repair rate, RR, pipeline leakage form and leakage area are based on the statistical data of pipeline damage in the Northridge earthquake [26]. The basic calculation flow of the PDD model is shown in Figure 3.
The nodal demand pressure was set to 20 m; after simulating n times, the simulation results of the required water pressure were selected when the post-earthquake water pressure of the node was lower than the normal state. Note that when the water pressure of all nodes meets the minimum required water pressure in the simulation, the water demand of the simulation is the same as the water demand under normal conditions, so the calculation of the water demand in this situation was neglected.
Calculate the average water demand of the nodes after an earthquake, see Equation (2):
Q i ¯ = 1 k j = 1 k Q i k
where Q i ¯ is the average post-earthquake water demand of node i, and Qik is the water demand of the node I obtained by the k-th simulation. The post-earthquake water shortage rate of the node is calculated by Equation (3):
Δ q i = Q i Q i ¯ Q i
where Δqi is the post-earthquake water shortage rate of node i, and Qi is the water demand of node i under the normal condition. The greater the post-earthquake water shortage rate of the node, the greater the degree of influence of the node.

2.2. Emergency Support Level

The safety function guarantee level of an emergency support infrastructure, such as urban emergency traffic, water supply, power supply, communication, etc., specified in Article 6.1.2 of the “Urban Comprehensive Disaster Prevention Planning Standards” [27] in China is divided into Ⅰ, Ⅱ and Ⅲ. The classification of emergency support level reflects the importance of urban emergency support objects after an earthquake. We associated the key emergency support objects of WDN with node users, which can reflect the post-earthquake significance of nodes. According to the emergency support objects corresponding to the nodes, the nodes’ levels were determined, and the emergency support levels of the nodes were scored (Table 1).

2.3. Node Betweenness

A graph can be represented by Equation (4) [28]:
G ( V , E )
where V is a non-empty vertex set, and E is the edge set.
A graph is uniquely determined by the association between its vertices and edges, and also by the adjacency between its vertices. The matrix A = A(G) = (aij)n×n is called the adjacency matrix of graph G, where aij is used to represent the number of edges between vertex vi and vertex vj, which may be 0,1.
According to the adjacency matrix of a graph, the betweenness of a graph can be calculated. The betweenness is divided into two types: node betweenness and edge betweenness. The minimum path between node i and node j in the network will pass through some nodes; if a node k is passed by many shortest paths, it means that the node k is very important in the network. The importance or influence of node k can be expressed by the node betweenness, which is defined as:
B k = i , j C d ( i , j ) k d ( i , j )
where Bk is the betweenness of node k, d(i,j) is the shortest path numbers between node i and node j, and d(i,j)k is the minimum number of paths between node i and node j passing through node k.

2.4. Evaluation System of Disaster Prevention Impact Factors for Node Users

Factor index set z = {z1, z2, …, zk}, k = {1, 2, 3} is selected, where z1 is the node water shortage rate, z2 is the earthquake disaster prevention emergency support level score, and z3 is the node betweenness. The influence value { s i z k } under each factor indicator, indicating the influence value of node i under the factor indicator zk, is calculated. According to the influence value, the relative influence matrix F z k of the node under the factor index zk is calculated by the formula (6), which measures the relative influence of each node under one index.
F z k = { f 11 z k f 12 z k f 1 n z k f 21 z k f 22 z k f 2 n z k f n 1 z k f n 2 z k f n n z k }
where f i j z k is the relative influence value of node i and node j under the factor indicator zk, and n is the total number of nodes. When i = j, f i j z k has no practical meaning and takes a value of 0; when ij, it is calculated by Equation (7):
f i j z k = { 1 0.5 0 s i z k / s j z k > 1 s i z k / s j z k = 1 s i z k / s j z k < 1
In order to obtain the relative influence vector of the node under the factor index zk, the row vector elements of the relative influence matrix F z k of the node are summed by Equation (8):
F z k ( i ) = i = 1 n f i j z k
The relative influence vector of the node only describes the relative influence of the node under one evaluation index. To synthesize the relative influence value of the node under different indicators, the node’s comprehensive influence matrix F is constructed by using the relative influence vector of the node under different indicators (9):
F = { F 1 z 1 F 1 z 2 F 1 z k F 2 z 1 F 2 z 2 F 2 z k F n z 1 F n z 2 F n z k }
The sum of each row of the node’s comprehensive influence matrix F is summed to obtain the comprehensive relative influence values F1, F2, …, Fn of the nodes, which reflect the comprehensive relative influence of the nodes under different factors. Finally, the normalization method based on degree of membership is used to process the node’s comprehensive relative influence values to obtain the node users’ impact factor values.
The Gaussian membership function expression is:
F i ¯ = e ( F i c ) 2 2 σ 2
c = ( 1 + ε ) F m a x
σ = ( 0.5 + ε ) F m i n m a x
where F i ¯ is the node i user’s disaster prevention impact factor value of WDN, Fmax is the maximum value of the comprehensive relative influence set {F1, F2, …, Fn}, Fmin is the minimum value of the integrated relative influence set {F1, F2, …, Fn}, ε is the normalized parameter, and the value range is ε ( 0 , 1 ) .

3. Calculation of Disaster Prevention Impact Factors for Node Users

3.1. Sensitivity Matrix

The pressure sensitivity of nodes reflects the influence of the pressure on other nodes when an abnormality occurs in one node. Assuming that there are n nodes in the WDN, the node pressure vector {P0} of WDN under the normal working condition is obtained by hydraulic analysis of EPANET software:
{ P 0 } = { P 1 0 , P 2 0 , , P n 0 } T
where Pn0 is the water pressure of the n-th node(m). Then, the initial matrix [P0] of node pressure is constructed by n-column {P0}. Comparing the node pressure changes calculated by 10%, 30% and 50% of the minimum user flow, the leakage of the minimum user flow that makes the clustering effect the best is selected. After calculation, we took 50% of the minimum node user flow as the node leakage, and the node pressure matrix [P] is:
[ P ] = [ P 11 P 12 P 1 n P 21 P 22 P 2 n P n 1 P n 2 P n n ]
where Pij is the pressure value of node i when the node j leaks.
The pressure sensitivity matrix [ΔP] is obtained by subtracting the node pressure matrix [P] from the initial pressure matrix [P0] of the node pressure.
[ Δ P ] = [ P 0 ] [ P ]
[ Δ P ] = [ Δ P 11 Δ P 12 Δ P 1 n Δ P 21 Δ P 22 Δ P 2 n Δ P n 1 Δ P n 2 Δ P n n ]
where Δpij is the pressure change value of the node i when the node j leaks.
The sensitivity matrix of nodal pressure was obtained by dimensionless normalization based on the extremum method.
l i j = Δ p i j m i n ( Δ p ) m a x ( Δ p ) m i n ( Δ p )
The comprehensive sensitivity of the node is calculated from Equation (18):
L i = j = 1 n | l i j | 2
where Li is the comprehensive sensitivity of node i. The comprehensive sensitivity of a node indicates the pressure sensitivity of the node to the flow change.
On the basis of clustering integration, according to the monitoring domain node classification, the intra-class normalization of the node comprehensive sensitivity is performed by Equation (19), which indicates the sensitivity of the node in the corresponding class.
L i ¯ = L i L m i n L m i n m a x , i C
where C is the monitoring domain category, L i ¯ is the comprehensive sensitivity of node i in class C after normalization, also called the intra-class sensitivity of node i, and Lmin is the minimum comprehensive sensitivity in class C.

3.2. Strategy for Disaster Prevention Sensor Placement

The principle of the disaster prevention sensor placement is not only to ensure the pressure sensitivity of the monitoring nodes, but also to consider the needs of the post-earthquake disaster prevention function of the node users. The placement should be arranged on the important user nodes or the adjacent nodes as far as possible. Thus, we considered the principle of the disaster prevention sensor placement and proposed a layout strategy for disaster prevention monitoring points.
The node comprehensive sensitivity of each class of the monitoring domain and the corresponding impact factor weights can be obtained from Equations (19) and (10). Next, the appropriate threshold was selected to consider the degree of sensitivity at the monitoring point:
L i ¯ > ψ , i C
where ψ is the given threshold. Then, the weights of node disaster prevention impact factors corresponding to the node of Equation (20) are ranked from large to small, and the node corresponding to the maximum disaster prevention impact factor weight is selected as the node for installation of monitoring equipment.
It should be noted that the choice of threshold directly affects the sensor placement, so a suitable threshold could balance the sensitivity and the influence of the node. In actual calculations, the threshold should be adjusted in conjunction with the scale of the WDN.

4. Simulation and Comparison

Figure 4 shows a water distribution model based on actual data from a region in China [29], which has 53 nodes and 78 pipelines. Node 50 is a reservoir; node 51, node 52, and node 53 are water source nodes; and the total hydraulic heads are 76.95 m, 74.31 m, 83.58 m and 81.98 m, respectively. We assumed that pipes with a diameter greater than 600 mm are riveted steel pipes, while pipes with a diameter less than 600 mm are rubber sealed cast iron. It needs to be mentioned that all methods were based on Matlab software, and the hydraulic calculation was called epanet2.dll. For the basic data of the WDN, see [29]:
Based on the hydraulic analysis, the nodal pressure sensitivity model was used to obtain the sensitivity matrix, and the extreme value processing method was used to normalize the node pressure sensitivity matrix. Xunjian Wang, Zengyi Wang [30] suggested that the setting of the WDN monitoring points should be about 1 7 ~ 1 6 of the total number of nodes, and the number of monitoring points of a large pipe network is smaller than that of a small pipe network. Therefore, the ratio of the monitoring points number in this paper to the total number of nodes was selected as 1/7, that is, the number of monitoring points is 7.

4.1. Disaster Prevention Impact Factors of Node Users

4.1.1. Node Water Shortage Rate

According to the table of seismic intensity in China (GB/t 17742-2008) [31], we chose the peak ground velocity (PGV) of the IX seismic intensity to calculate the average repair rate, RR, in which the average repair rate of pipelines with diameter greater than 600 mm is 0.0690, and that of pipelines with diameter less than 600 mm is 0.1254. The PDD post-earthquake hydraulic analysis model was used to simulate a total of 25 times, and all nodes satisfied the water requirements in the sixth, 10th, 12th, 19th and 20th simulations. Therefore, when calculating the water shortage rate, we did not consider the results of the above five simulations. The remaining 20 simulation results of the node average shortage rate value are shown in Figure 5.

4.1.2. Node Emergency Support Rating

We assumed the emergency guarantee level of the node in this paper (Figure 6), and we analyzed the node user impact factor values of the WDN according to the node users’ impact factor evaluation system proposed in the above section.
The nodes with an emergency protection level score of 3 or more are mainly concentrated in the eastern, northeastern, western and northwestern regions of the city (Figure 6). The monitoring sensors should be placed as much as possible in these areas for better monitoring.

4.1.3. Node Betweenness

After calculation, the betweenness values of each node were obtained, as shown in Table 2.

4.1.4. Node User Disaster Prevention Impact Factor Weight

The node water shortage rate, the node emergency support level score and the node betweenness value were used as the impact factor indicators. The parameters of each node and the relative influence value were calculated by Matlab software programming, and the relative influence value of each node is normalized. The weight values of the node impact factor are shown in Figure 7.
Figure 7 shows that the darker the node color, the higher the weight of the node impact factor, and the more important the node. The weight values of the node users’ impact factor reflect the post-earthquake functional importance of the node. The nodes with relatively high impact factor weights were node 19, node 25, node 24, node 23, node 22, node 8, node 34 and node 41. These nodes are very crucial for the shelter of residents after an earthquake and the rescue of secondary disasters, which should be protected.

4.2. Sensor Placement for Disaster Prevention in WDNs

4.2.1. Monitoring Domain Division Based on Three-Way Decision Cluster Integration

The K-means clustering method, fuzzy C-means clustering method, pedigree clustering method and dynamic classification method are common clustering methods for sensor placement. Each method has its own decision-making value. Based on the three-way decision cluster integration method, we integrated the K-means clustering method, fuzzy C-means clustering method, pedigree clustering method and dynamic classification method to divide the monitoring area of the WDN.
The three-way clustering integration method integrated four clustering methods to divide the monitoring domain of the WDN. As the three-way decision cluster integration method is not the focus of this article, the details can be found in [32]. The clustering results obtained by the K-means clustering method were taken as a reference to carry out cluster label matching, and the nodes that clearly belong to the corresponding category were divided into the positive domain of the corresponding category. After the voting matrix was constructed, the monitoring domain classification of the nodes was obtained according to the previous dividing rules, as shown in Table 3 and Table 4. The nodes in the boundary domain belong to the boundary of the class, i.e., the nodes divided into the boundary domain can be effectively monitored by different monitoring points at the same time.
The three-way decision integration clustering method integrates two-way decision clustering results into one three-way decision clustering result, which makes each category more uniform (Table 3 and Table 4). When it is impossible to determine that an object belongs to a certain class of clusters, it is divided into boundary regions, indicating that the object may belong to multiple clusters, thus describing its fuzziness.

4.2.2. Sensor Placement for Disaster Prevention in WDNs for Important Users

According to the classification of the monitoring domain and the comprehensive sensitivity of nodes, the sensitivity threshold is set as 0.85, so that we obtain the sensor placement nodes, which are node 44, node 14, node 23, node 10, node 32, node 20 and node 6. MapInfo software was used to draw the final placement diagram, as shown in Figure 8.

4.3. Sensor Placement in WDNs Based on Traditional Models

In traditional sensor placement, clustering algorithms are usually used to classify monitoring points based on the sensitivity matrix. Clustering centers are used as monitoring points where sensors are placed and corresponding categories are monitoring areas. The sensitivity curves of some nodes are shown in Figure 9, and the pedigree clustering tree is shown in Figure 10. It should be noted that the fuzzifier parameter of the fuzzy C-means clustering method was 2, and the membership values are shown in Table 5. The membership degree of the dynamic classification method was determined by the entropy weight method.
In this paper, we used eight methods for sensor placement in WDNs, and the sensor placements by the eight methods are shown in Table 6. In addition to the four clustering algorithms used before, the node correlation coefficient method, the effective independence method, the improved Fisher information matrix maximization criterion and the impact ranking method were also used. The idea of voting was used to integrate the monitoring information as well as to fuse the consistency of various algorithms. The node numbers in parentheses indicate the nodes strongly related to the front nodes, which can also be used as monitoring points for the sensor installation.
Finally, based on the idea of voting, we carried out a statistical analysis of sensor placement nodes in Table 6, selected the top 10 nodes as the alternative monitoring nodes, and the placement on the basis of the division of monitoring areas is shown in Figure 11.
According to the hydraulic data of the WDN, it can be seen that pressure value of node 32 is the lowest, which is 19.32 m under the standard working condition, and the second lowest nodal pressure value is node 10, which is 20.13 m. The next lowest is node 41, for which the pressure value is 20.15 m. These three nodes are in a relatively low water pressure. The water demand of node 11 is the second lowest, at 97.2 m3/h. Comparing the hydraulic information, we found that the nodes with lower pressure and water requirement are more sensitive to the changes of pressure and water requirement caused by leakage. The monitoring sensors should be placed at the end of pipework, which is more sensitive to pressure, as far as possible, which is consistent with that described in Jianwen Liang [31].

5. Discussion

The sensor placement of the two models is as follows:
  • Sensor placement based on traditional models: node 40, node 11, node 41, node 10, node 32, node 20 and node 5.
  • Sensor placement for disaster prevention in WDNs for important users: node 44, node 14, node 23, node 10, node 32, node 20 and node 6.
In terms of the node location, Figure 7 and Figure 10 show that the sensors placed in node 40 and node 11 in the traditional model are close to node 44 and node 14, respectively, in the model presented herein. Node 10, node 32 and node 20 are selected as monitoring placements in both the two models. This means that the above three nodes are still suitable for monitoring the pressure changes of pipework as sensor placements, considering both sensitivity and node users’ impact factors.
In terms of node importance, the monitoring points node 40, node 11, node 41 and node 5 in the traditional model become node 44, node 14, node 23 and node 5 in the model presented herein. From Figure 6, it can be seen that the sensor placements based on the model presented herein are located at or near the nodes with high impact factors, indicating that the placements are more important in post-earthquake emergency support, network topology and water shortage rate.

6. Conclusions

In this paper, the post-earthquake function importance and network importance of nodes are introduced into the placement of monitoring points to realize the layout of disaster prevention monitoring points in WDNs, with good compatibility with health monitoring layout. Three-way decision clustering methods are introduced to comprehensively consider the decision value of integrating multiple clustering methods, which can be used as a way to consider the influence of different factors in sensor placement. The work in this paper has formed a technical tool for sensor placement for disaster prevention, which is of great significance for urban post-earthquake emergencies.

Author Contributions

Conceptualization, J.W., W.W. and D.M.; methodology, J.W. and W.W.; software, J.W. and Z.H.; validation, J.W.; formal analysis, D.M.; investigation, J.W.; resources, D.M. and W.W.; data curation, J.W.; writing—original draft preparation, J.W.; writing—review and editing, J.W., D.M. and W.W.; visualization, J.W.; supervision, D.M.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51678017), National Key R&D Program of China (Grant No. 2018YFD1100902-1), Research Project on Major Policy Theory and Practice of China Earthquake Administration (CEAZY2019JZ14), Key Project of Earthquake Engineering and Engineering Vibration Key Laboratory of China Earthquake Administration (2019EEEVL0501).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Liu, S.; Zhang, X. Investigation Report on Earthquake Disaster and Relief of Water Supply System Based on Wenchuan Earthquake; Tongji University Press: Shanghai, China, 2013; pp. 4–10. ISBN 978-7-5608-5060-3. [Google Scholar]
  2. Li, J.; Liu, W.; Bao, Y.F. Genetic algorithm for seismic topology optimization of lifeline network systems. Earthq. Eng. Struct. Dyn. 2008, 37, 1295–1312. [Google Scholar] [CrossRef]
  3. Zhang, C.; Lin, Z.Z.; Wen, F.S.; Ledwich, G.; Xue, Y.S. Two-stage power network reconfiguration strategy considering node importance and restored generation capacity. IET Gener. Transm. Dis. 2014, 8, 91–103. [Google Scholar] [CrossRef]
  4. Fan, W.L.; Hu, P.; Liu, Z.G. Multi-attribute node importance evaluation method based on Gini-coefficient in complex power grids. IET Gener. Transm. Dis. 2016, 10, 2027–2034. [Google Scholar]
  5. Wang, F.; Zheng, X.Z.; Chen, S.; Zhou, J.L. Emergency Repair Scope Partition of City Water Distribution Network: A Novel Approach Considering the Node Importance. Water Resour. Manag. 2017, 31, 3779–3794. [Google Scholar] [CrossRef]
  6. Kermarrec, A.M.; Merrer, E.L.; Sericola, B.; Trédan, G. Second order centrality: Distributed assessment of nodes criticity in complex networks. Comput. Commun. 2011, 34, 619–628. [Google Scholar] [CrossRef]
  7. Zareie, A.; Sheikhahmadi, A.; Jalili, M. Influential node ranking in social networks based on neighborhood diversity. Future Gener. Comput. Syst. 2019, 94, 120–129. [Google Scholar] [CrossRef]
  8. Fragiadakis, M.; Christodoulou, S.E.; Vamvatsikos, D. Reliability Assessment of Urban Water Distribution Networks Under Seismic Loads. Water Resour. Manag. 2013, 27, 3739–3764. [Google Scholar] [CrossRef]
  9. Cao, X.X.; Qie, Z.H.; Wang, W.Z.; Wu, X.M. Genetic optimization of flow monitoring points for water supply network fault diagnosis. China Rural Water Hydropower 2013, 1, 22. [Google Scholar]
  10. Chen, L.L.; Zhuang, W.T.; He, X. Pressure monitoring point layout of water supply network based on information maximization criterio. J. Shanghai Univ. Nat. Sci. 2015, 21, 640–647. [Google Scholar]
  11. Kowalski, D.; Kowalska, B.; Kwietniewski, M. Monitoring of water distribution system effectiveness using fractal geometry. Bull. Polsh Acad. Sci. Tech. 2015, 63, 155–161. [Google Scholar] [CrossRef] [Green Version]
  12. Kirstein, J.K.; Albrechtsen, H.J.; Rygaard, M. Topological clustering as a tool for planning water quality monitoring in water distribution networks. Water Sci. Tech. Water Supply 2015, 15, 1433–1435. [Google Scholar] [CrossRef] [Green Version]
  13. Nejjari, F.; Sarrate, R.; Blesa, J. Optimal pressure sensor placement in water distribution networks minimizing leak location uncertainty. Procedia Eng. 2015, 119, 953–962. [Google Scholar] [CrossRef] [Green Version]
  14. Soldevila, A.; Blesa, J.; Tornil-Sin, S.; Fernandez-Canti, R.M.; Puig, V. Sensor placement for classifier-based leak localization in water distribution networks using hybrid feature selection. Comput. Chem. Eng. 2018, 108, 152–162. [Google Scholar] [CrossRef] [Green Version]
  15. Cornell University. GIRAFFE User’s Manual Version 4.2; Cornell University: Ithaca, NY, USA, 2008. [Google Scholar]
  16. Shi, P.; O’rourke, T.D.; Wang, Y. Simulation of earthquake water supply performance. In Proceedings of the 8th National Conference on Earthquake Engineering, Oakland, CA, USA, 18–22 April 2006. [Google Scholar]
  17. Shi, P.; O’rourke, T.D. Seismic Response Modeling of Water Supply Systems. Available online: http://www.eng.buffalo.edu/mceer-reports/08/08-0016.pdf (accessed on 13 April 2019).
  18. Wang, Y.; O’rourke, T.D. Seismic Performance Evaluation of Water Supply Systems; Technical Report MCEER-08-0016; Multidisciplinary Center for Earthquake Engineering Research: Buffalo, NY, USA, 2006. [Google Scholar]
  19. Laucelli, D.; Giustolisi, O. Vulnerability Assessment of Water Distribution Networks under Seismic Actions. J. Water Resour. Plan. Manag. 2015, 141, 04014082. [Google Scholar] [CrossRef]
  20. Yoo, D.G.; Kang, D.S.; Kim, J.H. Seismic reliability assessment model of water supply networks. In World Environmental and Water Resources Congress 2013: Showcasing the Future; Patterson, C.L., Struck, S.D., Murray, D.J., Eds.; American Society of Civil Engineers: Reston, VA, USA, 2013; pp. 955–966. [Google Scholar]
  21. Han, Z.; Ma, D.H.; Hou, B.W.; Wang, W. Post-earthquake hydraulic analyses of urban water supply network based on pressure drive demand model. Sci. Sin. Technol. 2019, 49, 351–362. [Google Scholar] [CrossRef]
  22. Duda, R.; Hart, P.; Stork, D. Pattern Classification, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2001. [Google Scholar]
  23. Lingras, P.; Yan, R. Interval clustering using fuzzy and rough set theory. In IEEE Annual Meeting of the Fuzzy Information; IEEE: Piscataway, NJ, USA, 2004; pp. 780–784. [Google Scholar]
  24. Hu, L.C.; Yu, H. A voting-based clustering method for three decision-making branches. J. Chin. Comput. Syst. 2016, 37, 1741–1745. [Google Scholar]
  25. Wagner, J.M.; Shamir, U.; Marks, D.H. Water distribution reliability simulation methods. J. Water Resour. Plan. Manag. 1988, 114, 276–294. [Google Scholar] [CrossRef] [Green Version]
  26. Jeon, S.S.; O’Rourke, T.D. Northridge earthquake effects on pipelines and residential buildings. Bull. Seismol. Soc. Am. 2005, 95, 294–318. [Google Scholar] [CrossRef]
  27. GB/T 51327-2018. Standard for Urban. Planning on Comprehensive Disaster Resistance and Prevention; China Construction Industry Press: Beijing, China, 2018. [Google Scholar]
  28. Wang, H.Y.; Huang, Q.; Li, C.T. Graph. Theory Algorithm and Matlab Implementation; Beijing University of Aeronautics and Astronautics Press: Beijing, China, 2010; ISBN 978-7-8112-4940-8. [Google Scholar]
  29. Lin, Y.Y. Design of Real-Time Monitoring System for Urban. Water Supply Network and Optimal Layout of Monitoring Points; South China University of Technology: Guangzhou, China, 2013. [Google Scholar]
  30. Wang, X.J.; Wangm, Z.Y. On the selection of pressure monitoring points in water supply network. China Water Supply Drain. 1989, 5, 9–12. [Google Scholar]
  31. GB/T 17742-2008. The Chinese Seismic Intensity Scale; China Standard Press: Beijing, China, 2008. [Google Scholar]
  32. Wu, J.J. Study on Monitoring Placement of Urban. Water Supply Network Based on Disaster Prevention Function; Beijing University of Technology: Beijing, China, 2019. [Google Scholar]
  33. Jain, A.K. Data clustering: 50 years beyond K-means. Pattern Recogn. Lett. 2010, 31, 651–666. [Google Scholar] [CrossRef]
Figure 1. The structure of this paper.
Figure 1. The structure of this paper.
Sustainability 12 00723 g001
Figure 2. The evaluation system structure for disaster prevention impact factors of node users.
Figure 2. The evaluation system structure for disaster prevention impact factors of node users.
Sustainability 12 00723 g002
Figure 3. Post-earthquake hydraulic model based on pressure-driven demand in a water distribution network.
Figure 3. Post-earthquake hydraulic model based on pressure-driven demand in a water distribution network.
Sustainability 12 00723 g003
Figure 4. Water distribution network model.
Figure 4. Water distribution network model.
Sustainability 12 00723 g004
Figure 5. Node water shortage rate.
Figure 5. Node water shortage rate.
Sustainability 12 00723 g005
Figure 6. The score map of node emergency support level.
Figure 6. The score map of node emergency support level.
Sustainability 12 00723 g006
Figure 7. The ordering of node user impact factor weights.
Figure 7. The ordering of node user impact factor weights.
Sustainability 12 00723 g007
Figure 8. Sensor placement for disaster prevention in water distribution networks for important users.
Figure 8. Sensor placement for disaster prevention in water distribution networks for important users.
Sustainability 12 00723 g008
Figure 9. Sensitivity curves of some nodes.
Figure 9. Sensitivity curves of some nodes.
Sustainability 12 00723 g009
Figure 10. Pedigree clustering tree.
Figure 10. Pedigree clustering tree.
Sustainability 12 00723 g010
Figure 11. Sensor placement based on traditional models.
Figure 11. Sensor placement based on traditional models.
Sustainability 12 00723 g011
Table 1. The corresponding score of emergency support levels of the nodes.
Table 1. The corresponding score of emergency support levels of the nodes.
Level of Emergency SupportScore
4
3
2
1
Table 2. Node betweenness.
Table 2. Node betweenness.
Node NumberBetweennessNode NumberBetweennessNode NumberBetweenness
1102.001861.9235197.16
2565.8919283.1336289.28
3378.122078.9937218.11
4221.1621486.1538192.14
5130.4022137.153998.22
6306.8823411.8340132.90
7250.4624723.6241334.51
8299.6925500.1242121.37
9188.8726189.1543225.03
104.0027416.1644270.03
11102.0028288.114540.35
12230.4029398.5546199.00
1380.2030375.984744.71
14196.963134.0848194.53
15238.063234.0849102.00
16355.5833102.00
1762.1634396.83
Table 3. Positive domain of three-branch clustering integration.
Table 3. Positive domain of three-branch clustering integration.
CategoryPositive Domain
1424724304348
2111214171925
3492328293334353641
4102631
53239
61520212227
73451612
Table 4. Boundary domain of three-branch clustering integration.
Table 4. Boundary domain of three-branch clustering integration.
CategoryPositive Domain
1,5373840444546
2,4,57891318
2,4,76
Table 5. Membership values.
Table 5. Membership values.
ClassificationN1N2N3N4N5N6N7N8N9N10
10.1180.1160.0300.1720.0680.1060.0420.0150.0110.132
20.0290.0270.0060.0820.0140.0240.0120.0050.0030.116
30.0440.0460.0090.0770.0190.0650.0510.0180.0150.146
40.1010.1110.0180.1120.0350.2810.4860.8240.8500.163
50.0910.0960.0220.1350.0460.1490.0880.0270.0220.145
60.0720.0770.0140.0980.0280.1540.2660.0900.0820.161
70.5450.5270.9020.3240.7910.2210.0540.0210.0160.136
ClassificationN11N12N13N14N15N16N17N18N19N20
10.1170.0840.0380.0210.0280.1380.0270.0160.0130.015
20.0620.0360.0150.0080.8760.0310.0090.0050.0040.934
30.1140.0940.0610.0430.0170.0330.0540.0210.0300.009
40.2320.3320.4750.3740.0170.0560.3110.7760.1470.009
50.1400.1190.0700.0470.0230.0830.0710.0310.0340.012
60.1890.2230.2960.4830.0170.0460.4970.1300.7560.009
70.1460.1120.0450.0240.0230.6130.0310.0220.0150.012
ClassificationN21N22N23N24N25N26N27N28N29N30
10.1590.2190.5270.0670.0220.0470.0290.6590.4640.066
20.4360.2480.0260.0130.0070.0190.8860.0180.0290.009
30.0590.0800.0470.0790.0700.1640.0140.0310.0580.050
40.0650.0860.0460.0720.1260.1720.0140.0380.0510.036
50.1020.1460.2150.6290.0700.1060.0210.1240.2560.758
60.0640.0840.0470.0950.6830.4450.0140.0360.0540.046
70.1150.1380.0910.0450.0220.0470.0210.0950.0870.034
ClassificationN31N32N33N34N35N36N37N38N39N40
10.0580.0580.2920.5030.7710.7720.0470.0230.0560.019
20.0280.0270.0970.0580.0140.0150.0160.0090.0270.007
30.2180.3420.0880.0520.0230.0250.6390.7580.4440.799
40.1730.1330.0990.0620.0270.0250.0520.0430.1060.035
50.1210.1260.1500.1190.0710.0910.1320.0660.1230.057
60.3470.2620.0960.0580.0250.0250.0820.0820.1970.068
70.0560.0530.1790.1500.0690.0470.0320.0190.0480.016
ClassificationN41N42N43N44N45N46N47N48N49
10.8210.1360.0740.0480.0420.0260.1870.0890.373
20.0150.0280.0240.0170.0160.0100.0620.0220.048
30.0190.2430.4860.6310.6500.7420.1770.3810.060
40.0210.0780.0680.0530.0560.0450.1100.0770.078
50.0540.3430.2000.1360.1120.0700.2220.2640.129
60.0200.1010.1020.0840.0920.0860.1220.1090.070
70.0500.0710.0460.0320.0300.0210.1190.0590.242
Table 6. Comparison of various methods.
Table 6. Comparison of various methods.
MethodsSensor Placements
K-means clustering [33]48, 8, 41, 10, 32, 27, 3
Pedigree clustering10, 15 (20), 23 (29), 1 (2), 11 (12), 8 (9), 40 (46)
Fuzzy C-means clustering41 (36), 20 (27), 40 (38), 9 (8), 30 (24), 19 (25), 3 (5)
Dynamic classification10, 40, 42, 8, 41, 5, 1
Node correlation coefficient method10 (13), 11 (12), 15 (20), 47 (42), 33 (34), 4 (5), 31 (26)
Effective independence method39 (32), 10, 15 (20), 13 (14), 22, 29 (23), 3 (5)
Improved Fisher information matrix maximization criteria [32]10, 33, 4, 47, 11, 22, 15
Influence ranking method10, 15, 20, 39, 31, 32, 45

Share and Cite

MDPI and ACS Style

Wu, J.; Ma, D.; Wang, W.; Han, Z. Research on Sensor Placement for Disaster Prevention in Water Distribution Networks for Important Users. Sustainability 2020, 12, 723. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020723

AMA Style

Wu J, Ma D, Wang W, Han Z. Research on Sensor Placement for Disaster Prevention in Water Distribution Networks for Important Users. Sustainability. 2020; 12(2):723. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020723

Chicago/Turabian Style

Wu, Jiajia, Donghui Ma, Wei Wang, and Zhao Han. 2020. "Research on Sensor Placement for Disaster Prevention in Water Distribution Networks for Important Users" Sustainability 12, no. 2: 723. https://0-doi-org.brum.beds.ac.uk/10.3390/su12020723

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop