Next Article in Journal
Review of Recent Drilling Projects in Unconventional Geothermal Resources at Campi Flegrei Caldera, Cornubian Batholith, and Williston Sedimentary Basin
Next Article in Special Issue
Resilience Assessment of Wind Farms in the Arctic with the Application of Bayesian Networks
Previous Article in Journal
Review of Technical Design and Safety Requirements for Vehicle Chargers and Their Infrastructure According to National Swedish and Harmonized European Standards
Previous Article in Special Issue
Resilience Assessment in Electricity Critical Infrastructure from the Point of View of Converged Security
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Case Study in View of Developing Predictive Models for Water Supply System Management

by
Katarzyna Pietrucha-Urbanik
1,*,
Barbara Tchórzewska-Cieślak
1 and
Mohamed Eid
2
1
Department of Water Supply and Sewerage Systems, Faculty of Civil, Environmental Engineering and Architecture, Rzeszow University of Technology, Al. Powstancow Warszawy 6, 35-959 Rzeszow, Poland
2
National Institute of Applied Sciences of Rouen-LMN, INSA-Rouen, 685 avenue de l’Université-BP 08, 76801 St. Etienne du Rouvray, France
*
Author to whom correspondence should be addressed.
Submission received: 10 April 2021 / Revised: 23 May 2021 / Accepted: 31 May 2021 / Published: 4 June 2021
(This article belongs to the Special Issue Critical Infrastructure Resilience Assessment and Management)

Abstract

:
Initiated by a case study to assess the effectiveness of the modernisation actions undertaken in a water supply system, some R&D activities were conducted to construct a global predictive model, based on the available operational failure and recovery data. The available operational data, regarding the water supply system, are the pipes’ diameter, failure modes, materials, functional conditions, seasonality, and the number of failures and time-to-recover intervals. The operational data are provided by the water company responsible of the supply system. A predictive global model is proposed based on the output of the operational data statistical assessment. It should assess the expected effectiveness of decisions taken in support of the modernisation and the extension plan.

1. Introduction

Within the frame of the modernisation of a regional water supply system (WSS) in Poland, an assessment of the effectiveness of the related decisions and actions was demanded by the regional authorities. The assessment of this case study was accompanied by some R&D activities to develop predictive models in support of the decision making and the management of the modernisation of the WSS.
The quality of a water supply service can be measured using different metrics. The availability and the reliability of the (WSS) are the most cited amongst these metrics. The availability measures the aptitude of the WSS to provide the water at the required quality and quantity on demand at instant “t”. The reliability measures the aptitude of the WSS to provide the water at the required quality and quantity during a period “T”. Methods of assessing the availability and the reliability of the water supply networks are well developed and widely published, among others in [1,2]. It is the same situation for case studies and WSS performance assessments [3].
Throughout the paper, the abbreviation WSS will refer to the water supply system rather than to the whole water supply service. Generally, the water supply services performance is judged upon the capacity to supply the demanded water quantity with the required pressure and the required quality (chemical-physical-biological) at the convenience of the consumers [4,5]. However, the WSS performance is just one dimension of the whole water supply service performance. For this reason, water supply utilities and other stakeholders are interested in models, methods, and standards that enable them to measure, assess, and predict the WSS availability and reliability.
The availability and the reliability are probabilistic concepts. They are not directly measurable. However, they can be indirectly measured thanks to the operational experience feedback of the WSS. The operational experience feedback provides the elementary data to develop the required predictive availability and reliability models. The WSS failure and recovery records that each water supply utility has constitute valuable sources of the elementary data needed to measure the availability and the reliability. Subsequently, water utilities put strong emphases on failure and recovery data collection, storage, and analyses within the optic of developing predictive models.
The statistical analysis of the WSS failure and recovery data permits the extraction of the failure and recovery characteristic parameters. Dependant on the nature of the collected data, the characteristic parameters can be the occurrence frequencies (failures, recoveries, etc.), the time intervals (time before failure, time before recovery, etc.), or both. Thanks to these parameters, the WSS availability and reliability can be assessed. Besides, some statistical predictive models can be developed to assess the future trends of the characteristic parameters.
The available failure and recovery data are generally impacted by technical, procedural, and managerial performance effectiveness [6]. Including all these dimensions is out of the scope of our paper, as the available data do not allow the integration of such dimensions in the assessment. So, the assessment will only consider the technical dimension.
The paper presents a case study of a given regional WSS in Poland. The objective of this case study is to assess the failure and the recovery characteristic parameters, in view of building up some predictive capabilities. The predictive capabilities are intended to support decision making in the management process of the modernisation and the extension of the regional WSS. The assessment is carried out in two phases: the statistical analysis of the operational data trends and the development of a predictive global model.
The paper describes, in Section 2, the general lines of the methodology of the study.
In Section 3, the paper presents a general description of the case study. The general description of the WSS covers the pipe-lines network: pipe types, geometrical characteristics, material specifications, and lengths in km.
In Section 3, the paper introduces the statistical assessment of both the failure database and the recovery database.
In Section 4, the predictive global model developed for the case study is proposed.
In Section 5, the paper shows the principal conclusions and findings of the assessment.
In Section 6, the paper summarises a global synthesis of the work to conclude.

2. Methodology

The background of the work is to assess the failure-recovery data of a regional WSS in Poland in view of building up first order predictive availability and reliability models. These preliminary predictive models were required to serve in predicting gross control indicators to assess the effectiveness of the decisions and the actions taken within the modernisation and the extension of a regional WSS. The developed models are preliminary as the available database is of a small size, heterogeneous, and tainted with high uncertainty level.
The case study of this WSS covered two tasks: the statistical analysis of the failure-recovery data and the development of a preliminary predictive models exploiting the outputs of the statistical analysis. Our work methodology consisted of distinguishing three spheres of activities:
  • Identification and classification of the information available about the regional WSS.
  • Identification of the adequate statistical analysis methods to be exploited according to the specification of the available input data and their statistical quality.
  • Identification of the most adequate concepts to build up preliminary predictive models of availability and reliability, given the physical nature of the WSS and the nature of the outputs of the statistical analysis of the operational data.
The WSS specifications are given in Section 3.1. The specifications cover the extension of the WSS, the classification of the pipes, and the lengths and materials of each class of pipes. The specifications of the available failure and recovery data and time intervals covered by each dataset are determined.
The differences in nature between the available failure and recovery datasets are put in evidence. Subsequently, the choice of the adequate statistical analysis methods is done for each dataset. Both datasets were processed differently, but the main output has the same nature: an overall failure rate and an overall recovery rate.
These overall transition rates are the required outputs to build up the preliminary predictive models of the whole WSS. The details are given in Section 3.2.
Given the nature of the available output data resultant from the statistical assessment in Section 3, we had to develop a binary model considering globally two events: failure and recovery.
Given the limited amount of data, it was not conceivable to treat the WSS as a network even with limited number of nodes and edges. Subsequently, we had to opt for an overall predictive model. We had to focus on first order models as our immediate target was to obtain control indicators to assess the effectiveness of the decisions and the actions taken within the program of the modernisation and the extension of the regional WSS.
Although a global model does not allow for the description of individual segments of the network, it responds satisfactorily to the needs of assessing the effectiveness of the ongoing modernisation and the extension work in the WSS. The details about the proposed models are given in Section 4.

3. Input Data Description

The input data are distinguished in two classes: network descriptive characteristics and operational feedback failure-frequency and recovery-time data.
Network descriptive characteristics cover the following groups:
  • General information about the objects (locations, ages, etc.).
  • Technical data about the objects (types of objects considering their functionality, geometry, materials, technology, operational conditions, etc.).
  • Data on failure (type of event, cause, mechanisms, etc.).
  • Data on the effects and consequences of failure (type, damage severity, extensions, etc.).
  • Additional information (report date, environmental data, cost).
The operational feedback failure-frequency and recovery-time data is presented in two independent databases with two different structures: the failure-frequency database and the recovery-time database. The structural differences between the two databases may be the result of the separation between the operational surveillance task and the maintenance task due to some local organisational conditions.
The data type is different from one database to the other. In the failure database, one counts number of failures per given time-interval while one measures intervals of time-to-recover in the recovery database. Both measurables are random but different in nature.
Regarding the scope of the paper, the final outcomes we are looking for are an “overall failure rate” and an “overall recovery rate” to globally characterise the WSS. They may be called “global failure rate” and “global recovery rate” as well. They are sometimes called system equivalent parameters [7], especially if they are calculated based on the characteristic parameters of the system’s constituents. These global characteristic parameters are necessary to assess the quality of the water supply service, to support the decision making regarding the modernisation of the WSS, and to develop global predictive models of availability and reliability.
The global characteristic parameters will then be designated as:
  • The overall failure rate of the WSS, λ W S S .
  • The overall recovery rate of the WSS, μ W S S .
It is important to underline that the estimated characteristic parameters ( λ W S S , μ W S S ) are describing the global functional behaviour of the WSS. They cannot be used to characterise a zonal portion of the WSS, a single type of the pipes (main, distribution, connection), or for maintenance optimisation purposes. However, they are useful in automatic monitoring of event occurrences in the network. The predictive models that we are looking for are global. These are point-like models.

3.1. Network General Description

The analysed water supply network is a ring system formed by four mains with the total length of about 50 km. The total length of the water pipes including the mains is about 900 km. Distributional pipes represent ~60% of the total length of the pipes, the length of the water supply connections represents ~35%, and the mains represent only ~5% with a length increasing rate of about 1.3% in comparison to the total length of the network, in the last year of the analysis.
The pipes made from:
  • Polyethylene (PE) represent 32% (289.5 km) of the total length of the water pipes.
  • Polyvinyl chloride (PVC) pipes represent 28% (254.6 km).
  • Galvanized steel (GS) pipes represent 24% (213.6 km).
  • Cast iron (CI) pipes represent 11% (103 km).
  • While those from hard steel (HS) represent about 3% (30.6 km), and the ductile iron (DI) pipes represent less than one percent (~10 km).
During the last decade, the pipes made from PE and PVC are increasingly used.
The age distribution of the mains and the distributional network is:
  • Pipes aged less than 10 years represent 31.2% of the total operating water pipes.
  • Between 11 and 25 years represent 23.2%.
  • Between 26 and 50 years represent 31.8%.
  • Above 50 years represent 13.8%.
The share of the pipes aged less than 10 years is gradually increasing due to the replacement of the damaged pipes and the renewing activities following the expansion of the city districts and the creation of new settlements. To ensure better control on the network operational performance, the water company manages the network through monitoring and recording pressure and flow rate, in the first and second zones.

3.2. Failure Data Assessment

Failures are generally classified by cause, mechanism, and/or criticality. Root causes can be of different origins such as: design, construction, and/or operational conditions [8,9,10,11,12]. Failure analysis activities are very impacting on network everyday operational and maintenance decision making [13,14].
However, the failure database we have is limited. Its failure records contain only the number of failures per month in the network, over 18 years. The failures are classified by the pipe’s position, geometry, material, and type (main, distribution, connection). No direct data were available regarding the “time to fail” per month per segment of pipe. Information about failure causes or criticality are almost absent, inconsistent, or incomplete.
The assessment aims at the determination of an overall failure occurrences trend in the network considering the field operational reality and the continuous expansion of the network due to the expansion of the town itself and the modernisation of the WSS. The failure data cover eighteen years of operation.
Given the small size of samples per class of failures, we had to work out an overall failure rate, mixing failure data over all ages, failure causes, modes, mechanisms, and criticality. We kept only the functional classification: distribution, connection, and main. As it is an overall failure occurrence rate and given the heterogeneity, it cannot be used for maintenance or for safety purposes. It just counts failures occurrences.
Data quality analysis was mainly related to the correction of errors through the use of control assumptions checking arithmetic and logical relations among statistical data, their completeness, and their accuracy. In water supply practice, the collected data may be burdened with random errors, which results from mistakes in their collection and processing and systematic errors; therefore, to assess the credibility of data, it is necessary to perform substantive and formal verification of their truthfulness, which consists of random checking of reliability, comparing the consistency of related data and plotting data graphs.
Given the dataset we have, different failure analysis models and methodologies are possible: the cumulative failure distributions [15,16], the data-driven modelling techniques using neural network and neuro-fuzzy methods [17,18], the Weibull-exponential failure model and the constant probability replacement model [19], the probability of main failure method [20], or the polynomial regression methods.
The available failure database contains the number of failures per month. So, failure rates can be directly estimated using polynomial regression techniques, based on the number of failures of all pipes within a given sampling interval [1,2,6]. The choice of the unit sampling interval is important to guarantee an acceptable statistical quality of the estimated failure rate. The notion of “Guarantee” in probabilistic analyses is generally replaced by a statistical meaningful term which is “the statistical quality”. Regarding the regression and data fitting, the “statistical quality” is assessed using indicators such as: the variance, the most likelihood probability, and/or the least-square values. In our polynomial regression, we used the R-Squared indicator.
Some preliminary investigations were carried out using different sampling intervals: 1, 3, 6, and 12 months. The sampling interval of one year was selected as it showed an acceptable dispersion behaviour in terms of variance and suits the engineering use of the water sector.

3.2.1. Failure Data Regression Processing

Both linear and polynomial regressions were performed. Then, we only considered the polynomial regression of order n 2 when it gives higher statistical quality based on the R-Squared ( R 2 ) value, otherwise the average values are considered. This is the case with the distribution pipes failure rate in Figure 1a and the overall failure rate in Figure 1d where the 2nd order polynomial regression results were considered. Otherwise, the average values are considered for the connections, Figure 1b, and the mains, Figure 1c.
We recall that the overall failure rate is a global measure as it covers all pipes of all types, all ages, all geometries, and all materials, independently of the zonal variation.
The annual failure data of the distribution-pipes are traced over 18 years, Figure 1a. The average failure rate is accordingly 2.87 × 10 1   n o . o f   f a i l u r e s · k m 1 · a 1 . The polynomial regression estimates that the distribution-pipes failure rate can be determined by the following 2nd order polynomial equation, with R 2 = 91 % :
λ d i s = 0.0012 t 2 0.0452 t + 0.5777 ,
where t is given in years.
The R-squared value of 91 % indicates that the estimated polynomial is statistically of a good quality and can satisfactorily represent the variations of the distribution-pipes with the time within the vicinity of 18 year-interval.
The annual failure data of the connection-pipes are traced in Figure 1b. The best estimated polynomial regression is found to be a polynomial of 5th order with R 2 = 27.8 % , as following:
λ c o n = 7 × 10 6 t 5 + f 4 × 10 4 t 4 0.007 t 3 + 0.0675 t 2 + 0.009 t + 0.0194 ,
Given the low value of R 2 , the polynomial was rejected, and the average value of the failure rate is considered as representative of the connection-pipes failure trend in the vicinity of the 18 years, such as:
λ ¯ c o n = 3.62 × 10 1   n o . o f   f a i l u r e s · k m 1 · a 1
Similarly, the annual failure data of the main-pipes are traced in Figure 1c. The best estimated polynomial regression is found to be a polynomial of 5th order with R 2 = 49 % , as following:
λ m a i = 3 × 10 5 t 5 0.0015 t 4 + 0.0265 t 3 0.2085 t 2 + 0.6351 t + 0.7024 ,
Given the low value of R 2 , the polynomial was rejected, and the average value of the failure rate is considered as representative of the main-pipes failure trend in the vicinity of the 18 years, such as:
λ ¯ m a i = 1.06   n o .   o f   f a i l u r e s · k m 1 · a 1
It is worth mentioning that despite the low statistical quality of the failure data of the connections and the mains, their impact on the WSS overall failure is relatively minor as they represent less than 35% of the whole WSS pipes’ length. This is confirmed by the following.
The overall failure rate trend in Figure 1d shows a similar trend pattern as the distribution-pipes one in Figure 1a. The overall annual failure data of the WSS are traced over 18 years, Figure 2d. The average failure rate is accordingly 4.82 × 10 1   n o . o f   f a i l u r e s · k m 1 · a 1 . The polynomial regression estimation with the 2nd order polynomial, with R 2 = 70.3 % , is:
λ t o t = 0.0008 t 2 0.0309 t + 0.558 ,
The R-squared value R 2 = 70.3 % indicates that the estimated polynomial has an acceptable quality and can statistically represent the global failure trend of the WSS pipes around the vicinity of 18 years.
The numerical results of the statistical assessment are synthesised in Table 1. It would be useful to recall that R-squared R 2 was the only control criteria used to select the polynomial regression models. We considered that a polynomial regression model with R 2 > 70 % would be acceptable for a polynomial order n 6 . The control coupled-criteria ( R 2 , n ) is arbitrary and justified by a good applied statistics expert judgement. Investigating more advanced analysis techniques such as Fourier decomposition or stepwise regression were not justified given the small size of the dataset, the strong heterogeneity of the data, and its low statistical quality.
It can be noted that the average value of the annual failure rate in the mains ( 1.06   n o . o f   f a i l u r e s · k m 1 · a 1 ) significantly exceeds the recommended value for the mains, which should be less than 0.3   n o . o f   f a i l u r e s · k m 1 · a 1 [21]. The average values of the distribution-pipes and the connection-pipes failure rates do not exceed the recommended values given in the literature [21,22].
As we are interested in predicting the failure behaviour of the WSS after the 18 years of measures, we are then interested in the overall failure rate λ W S S in the vicinity of the 18 years.
Considering the time behaviour of the WSS overall failure rate given in Figure 1d, one can estimate the WSS characteristic overall failure rate in the vicinity of the 18th year λ w s s = 0.25 × 900 = 2.25 × 10 2   a 1 18.8   m 1 6.3 × 10 1   d 1 . This figure will be used in the global predictive models later.

3.2.2. Seasonality Analysis

The impact of the seasonality on the WSS failures has also been assessed using Statsoft [23]. Statsoft software exploits, amongst others, the autoregressive integrated moving average (ARIMA) model [24]. A preliminary analysis of the data inherent variance was conducted, so as to verify the significant differences between months (Figure 2).
It is necessary to recall that the most significant statistical seasonal impact is on the mains (Figure 2). The seasonal impact is much less significant on the overall failure rate as the contribution of the mains in terms of km is less than 5% of the total WSS 900 km.
The parameters of the ARIMA model (p = 1, d = 0, q = 0) are assumed based on the autocorrelation graphs, using the so-called Box–Jenkins approach [24], where: p is the number of autoregressive terms, d is the number of differences, and q is the number of moving averages. This is the first order ARIMA model. Using higher orders of ARIMA models or more complex Fourier decomposition analyses were justified neither by the statistical quality of the data nor by the low seasonality effect. We recall that we have only failure occurrences per month in our database, with no additional information neither about failure modes nor about failure criticality.
The basis of the implemented approach is presented among others in [24,25]. The analysis of the average failure rate (per month) indicates that the failure culmination occurred is in December (at p = 0.00251, with current effect F (11.33) = 3.5081).
Figure 3 presents main-pipes failure rate changeability with combinations of the last observations and previous predictions adjusted for seasonal fluctuations.
For the connections, the series is more stationary regarding the ARIMA model. Only the autoregression of the first order is important, the Holt-Winters model is well matched without the trend but with the seasonality. For the failure rate of the distribution pipes, a seasonal downward trend is observed. That can be particularly observed in the Holt-Winters model with a trend that seems slightly better estimated than in the ARIMA model.
Annual seasonality and autocorrelation can be worked out for the entire network. This is evidenced by the statistically significant coefficients of “normal”. The seasonal autoregression in the ARIMA model is in accord with the Holt-Winters model [25], without the trend but with the seasonal component.
The basis of the implemented approach is presented among others in [24,25]. The Holt–Winters exponential smoothing is a statistical forecasting method in time series univariate methods. Exponential smoothing models and auto regressive moving average model proved to be the most adapted to assess temperature seasonality in many comparative assessments such as in [26,27]. The statistical fundamentals are clearly exposed and applied in [28].

3.2.3. Temperature Impact on the WSS Failures Behaviour

Regarding the meteorological parameters, we observed that low temperatures increase both the failure frequency of water pipes and the recovery times. The high correlation of dependence between the cumulative failures and the months with temperatures below zero Celsius are shown in Figure 4.

3.3. Recovery Data Assessment

The recovery database records contain the component-repair time and the organisational-repair time for each recovery event. Being the sum of the component-repair time and the organisational repair time, the recovery times are classified by failure mode and pipe’s geometry, material, and type (main, distribution, connection).
The assessment aims at the determination of the global recovery trends of the network considering the field operational reality, the continuous expansion of the network, and the modernisation of the WSS. The recovery data cover ten years of operation.
Recovery can generally be classified according to: the failure mode [22], the operational conditions [29,30], the location [31], the logistic conditions [32,33], and/or the managerial context [34,35]. The recovery times and the effective repair times analysis methodologies are well developed and published as well [36,37,38,39].
Regarding the recovery database we have, it contains only time intervals in two classes: physical component-repair times and associated organisational times, over 10 years. There is no additional information in the dataset. The recovery data are measured in time-intervals, not in repair frequencies, as in the case of failure data above.
Then, the data is grouped by subclasses: annual, failure mode, pipe material, pipe diameter, and pipe type (connection, distribution, main).
We will precisely be interested in analysing the component-repair time ( T c ) and the organisation-repair time ( T o ) within the intention to estimate the WSS overall recovery-time, ( T r ). We recall that the recovery time ( T r ) is the sum of the component-repair time ( T c ) and the organizational-repair time ( T o ).
The WSS recovery time is generally impacted by many factors. The following factors are considered in the assessment: the maintenance process, pipes material and diameter, failure modes, seasonal variations, and the average atmospheric conditions.

3.3.1. Maintenance Process Impacts on the Recovery Time

In Figure 5, the meantime of component-repair ( T c ), the meantime of organizational-repair ( T o ), and the recovery meantime ( T r ) are shown.
Two observations can be stated:
  • The organisational meantime T o is significantly longer than the component-repair meantime T c by factor ~3.
  • The variability per year in the component-repair meantime T c is much more significant than in the organisational-repair meantime T o .
The resultant annual recovery meantime T r is mainly driven by the organisational and administrative characteristic of the maintenance activities and shows a weak variability over the monitored period of 10 years. The recovery meantime T r is estimated to be about 9.8   h per reparation intervention, averaged over the 10 years (Table 2).

3.3.2. The Recovery Time vs. the Pipe Type

The analysis of the factors influencing the recovery time started by classifying the pipes as: mains, distributional, and the connections. The three classes are distinguished by their size, i.e., diameters.
Figure 6 presents the annual recovery meantime per class of pipes.
Regarding the recovery meantime, we counted: for the mains (766 recovery samples within the considered period), for the distributional pipes (1928 recovery samples), and for the water connections (1576 recovery samples). The recovery meantime is 8.76 h in the mains, 7.67 h for the connections, and 7.85 h for the distributional pipes. The longest recovery meantime was 33 h for the distributional pipes. The approximate values of the recovery meantime in various water pipes are driven by the priority of repairs given to the mains and the longer repair time due to the materials and dimensions of the pipes. The analysis in function of the diameters showed shorter recovery meantime in the case of the connections compared to the pipes of larger diameter. Shorter recovery meantime is observed in the water connections due to the higher accessibility to the failure location, the smaller diameter, and the types of materials.

3.3.3. The Recovery Meantime vs. Pipe Material

The impact of the pipes’ material on the recovery meantime is shown in Figure 7.
The longest recovery meantime was recorded for the pipes made of grey cast iron. In addition, a long recovery meantime was observed for steel and asbestos-cement pipes, which were used to build water supply systems in the last century. For grey cast iron the maximum observed recovery meantime is 33 h vs. 24 h for steel.
In the case of other materials and plastics, which are currently used on a large scale in the construction of the water supply network, the recovery meantime is less than 16 h. For plastic materials, repairs are easier and do not require the mobilization of heavy special equipment.

3.3.4. The Recovery Meantime vs. Failure Mode

Another important factor determining the recovery meantime is the failure mode. Figure 8 presents the recovery meantime per failure mode over ten years.
It can then be stated that the shortest recovery time occurred for mechanically damaged pipes, amounting on average to 5.63 h. This mode of failures is due to earthworks. The shorter recovery meantime is explained by the faster failure location process. The longest recovery meantime occurred in pipes on which sealing is lost or pipe breaking occurred. It may reach in many cases the maximum value of 16 h. This can be explained by the pipe damage degree, the difficulty of failure localization, and/or the pipe’s dimensions.
Detailed analysis showed that the highest organizational-repair time was noted for the mains. The value of T o significantly exceeds the mean value of 7.15 h. The time of the total repair of the water supply system has reached the highest values for the mains. The average value of the total repair was 6.8 h.
One can then state that the characteristic WSS recovery meantime μ W S S is time constant in the vicinity of the interval of 10 years and equal to 7.9 h. That corresponds to a constant repair rate μ W S S equal to 1.2 × 10 1   h 1 2.9 × d 1 86.4   m 1 .

3.3.5. Temperature Impact on the WSS Recover Behaviour

It was also observed that longer recovery meantime was associated with the negative atmospheric temperatures. The recovery meantime in the case of the negative temperatures was generally longer than for the positive temperatures by about 1.2 h. Incidentally, longer time of recovery meantime was in the summer months, July and August; it was caused by a sudden increase in air temperature. In addition, slight correlation was observed regarding the influence of the precipitation on the recovery meantime. The failures occurring during days without rain or with little rainfall have been moved up to 8 h. Longer recovery meantime was observed in days with little snowfall and rainfall, by on average 1.5 h. The high recovery meantime around 16 h appeared in the case of heavy rainfall and snowfall.

4. Availability and Reliability Predictive Modelling

Based on the estimated WSS characteristic failure and the recovery rates ( λ W S S , μ W S S ), a preliminary point-like model can be proposed to predict the overall availability and the reliability of the WSS. The model is global and considers the whole WSS as only one localized point which is either available (on) or in recovery (off).
Within this global concept, the term “availability” refers to the supply of the water to every part in the town at the required quality and quantity on demand at time “ t ”. The unavailability refers to the lack of supplying the water at exactly the required perfect conditions somewhere in the network on demand at time “ t ”, even if the differences are minimal or insignificant. This is because critical and non-critical failures are both considered in the estimated parameters ( λ W S S , μ W S S ). Accordingly, the model is conservative as it overestimates the unavailability and underestimates the availability.
Similarly, the term “reliability” refers to the supply of the water to every part in the town at the required quality and quantity on demand over a time-interval “ T ”. While the unreliability refers to the lake of supplying the water at exactly the required perfect conditions somewhere in the network over a time-interval “ T ”, even if the differences are minimal or insignificant.
Accordingly, the model predicts the lowest level of the availability/reliability one may expect.

4.1. The Unavailability Model

The model is a binary point-like where the WSS is either available a t or unavailable u t at time t . It is a scholar model extensively explained in almost all textbooks introducing system’s reliability modelling and analysis, Appendix A. Accordingly, the unavailability u t and the availability a t of the WSS can be determined, respectively, by the following Equations (5) and (6):
u t = λ λ + μ 1 exp ( λ + μ t + u 0 . exp λ + μ t
a t = μ λ + μ 1 exp ( λ + μ t + a 0 . exp λ + μ t
where u 0 and a 0 are respectively the initial unavailability and the initial availability.
The WSS availability and unavailability have asymptotic values ( a , u ) when t equal to:
u = u t t =   λ λ + μ ,
and,
a = u t t =   μ λ + μ   ,
The estimated values of the characteristic transition rates in Section 3.2.1 and Section 3.3.4 are:
λ w s s = 18.8   m 1 μ W S S = 86.4   m 1
So, the estimated highest level of the WSS unavailability u is:
u = λ λ + μ = 18.8 18.8 + 86.4 18 %
It is worth noting that the WSS attends 99.9 % of its asymptotic unavailability u within 1.5   days. So, we can state that the WSS attends its asymptotic time condition after t = 1.5 days. The evolution with the time of the estimated unavailability u t of the WSS is traced in days, Figure 9.
In addition, it is worth noting that the estimated rates are given in “per month”. Subsequently, the failure rate value of 18.8   m 1 corresponds to 1.6   d Mean Time Before Failure (MTBF), Equation (15), while the recovery rate value of 86.4   m 1 corresponds to 0.347   d Mean Time Before Recovery (MTBR), in “day” units.

4.2. The Reliability Model

Using the same classical model as in Appendix B, one can determine the meantime before failure as:
MTBF ^ = 1 λ W S S = 1 18.8 = 0.0532   t m = 1.6   d ,
Equation (9) is expected according to the survival analysis in system reliability theory as λ W S S is time-constant parameter. This is an overall meantime, considering all failure modes over all pipe-types and materials and everywhere in the network. It does not tell if the failure is significant or not, if it is on a main-pipe or a small connection, or if it is in winter or summer.
Still, it is a significant indicator as it tells us that a failure signal will be monitored each 1.5 days on average and its significance should be controlled to take the right action.
Another similar indicator is the number of failure signals n expected per time interval T , which is the solution of:
P n T = 0 T P n 1 τ · λ d τ · E x p λ W S S T τ ,
where:
P n T : is the probability of having n failure signals, all failure types included, within an interval T .
P 0 T : is the probability to have a zero-failure within the interval T and is equal to
P 0 T   = E x p λ W S S T ,
The solution of the previous Equation (10) is:
P n T = λ W S S T n n ! E x p λ W S S T ,
The estimated distribution P n T of having n failure-occurrences with the interval T is traced in Figure 10.
The distribution P n T of having n failure-occurrences with the interval T allows us to determine two interesting measurables that are easily monitored and useful in assessing the performance of the WSS and the water supply system.
The first of these measurables is the expectation number n ^ of failure-occurrences within an interval T , such as:
n ^ T   = E n = n = 1 n · P n T ,
The expectation value n ^ is then determined as such:
n ^ T   = n = 1 n · λ W S S T n n ! E x p λ W S S T = λ W S S T · E x p λ W S S T · n = 1 λ W S S T n 1 n 1 ! ,
Considering that:
n = 1 λ W S S T n 1 n 1 ! = E x p λ W S S T ,
Then expectation value n ^ is determined by:
n ^ T = λ W S S T = T MTBF ^ ,
The 2nd measurable is the most probable interval T M P within which n failure-occurrence may be observed. This most probable interval occurs when:
d d T P n T = 0 ,
We get than:
1 λ W S S T n · λ W S S T n 1 n 1 ! · E x p λ W S S T = 0 ,
Amongst the three possible solutions, the interval T M P that corresponds to the highest probability to observe n failure-occurrences, P n T M P , is then:
T M P = n λ W S S = n · MTBF ^ ,
Substituting T M P in Equation (17) results in the probability value P n T M P for the time-interval T M P , as following:
P n T M P =   n n n ! e n ,
The most probable intervals to observe n = 1 ,   2 ,   3 ,   4 failure-occurrences and their corresponding occurrence probabilities are given in Table 3.
To synthesise the learnings from Table 2, we underline that:
  • One should expect to monitor n failure-signals, each 1.6 n days. The monitored signals do not indicate neither the mode nor the criticality of the failure in terms of water supply unavailability.
  • The highest probability of monitoring, n failure-signals within intervals equal to 1.6 n days is equal to n n n ! e n .

4.3. Measurable Indicators

The global model presented in the paper allows us to propose four measurable indicators to assess the WSS and the water supply service performance.
The 1st indicator is the expected asymptotic unavailability (Equation (5)), the second is the meantime before failure (Equation (9)), the third is the expected number of failures n within an interval T (Equation (16)), and the fourth is the expected interval T M P to observe n failures (Equation (19)). These four indicators can directly be estimated by the statistical analysis of the operational experience feedback data.
We should underline that the values estimated by the mathematical model will be called “expected values” while the term “estimated values” will be used for the values that are resulted from the statistical analyses of the operational data.

5. Predictivity and Control Indicators

In the preceding section, we have illustrated how the predictive models may help in developing indicators that are useful in controlling the WSS and the water supply services performance. The predictive quality of a mathematical model is conditioned by the founded concept of the model and the statistical quality of the basic data issued from the operational experience feedback. The predictive model proposed in the paper is based on:
  • A binary functional concept that conceives the whole network as one single point (node) where only failure and recovery occur independently.
  • Real operational experience feedback data covering 18 years of failure observations and 10 years of recovery data. Failure and recovery data are structured in different manners.
Even though both failure-recovery data have in many aspects a fair statistical quality, the resultant predictive indicators are still significant even if it is attached with some statistical uncertainties.
In Section 4, we presented four indicators: the asymptotic unavailability u , the meantime before a failure is monitored MTBF ^ , the expected number of failures n ^ T within an interval T, and the most probable interval T M P to observe n failure-occurrences. These four indicators are interested because they are field direct measurables. They can be of great interest as the modernisation of the WSS will imply the use of the ICT for in-line monitoring and control activities.
These global indicators do not distinguish between failures severity, type, or location. However, they are very indicative regarding the qualification of the modernisation actions and the shared efforts between equipment technology ( λ ) and maintenance ( μ ) logistic organisation.
The comparison between the operational estimated and the model expected values will allow us to validate and to better adjust the models. Subsequently, these predictive models present increasingly valuable support to the decision making in managing the modernization and the expansion the WSS and the water supply service performance.
As our models have been conceived to predict the indicators mentioned above based on the operational database, the validation of the model will be carried on after updating the operational databases including failure-recovery post-modernisation data.
For the time being, the conclusion of the assessment is still a predictive indicator.

6. Conclusions

Within the frame of the modernisation of a regional water supply system (WSS) in Poland, some R&D activities were conducted to develop predictive models to support decision making in managing the modernisation and the urban extension of the WSS. The modelling activities should be based on and validated by the operational experience feedback.
The operational experience feedback is generally available in the form of failure and recovery databases. The paper presents a case study to assess the failure and the recovery characteristic parameters of a given regional WSS in Poland. The assessment is carried out in two phases: the statistical analysis of the available operational data and the development of a predictive global modelling capability. The concerned databases are failure and recovery databases. In the failure database, one counts number of failures while one measures intervals of time-to-recover in the recovery database. Both measurables are random but different in nature.
Regarding the scope of the paper, the final outcomes we looked for are an “overall failure rate” and an “overall recovery rate” to globally characterise the WSS. These characteristic parameters are necessary to develop global predictive models of availability and reliability, to assess the quality of the water supply service and to support the decision making regarding the modernisation of the WSS.
Regarding the recovery database we have, it contains time intervals in two classes: physical component-repair times and associated organisational times, covering 10 years. The analyses considered the component-repair time ( T c ) and the organisation-repair time ( T o ) within the intention to estimate the WSS overall recovery-time, ( T r ). The variability per year in the component-repair meantime T c is much more significant than in the organisational-repair meantime T o , including the seasonality effect. However, the major contribution in the overall recovery-time is so far from the organisational repair time. Accordingly, one can then state that the characteristic WSS recovery repair rate μ W S S is equal to ~ 2.9 × d 1 .
Once the overall characteristic parameters ( λ W S S , μ W S S ) are determined, one can then propose a classical model to predict unavailability and reliability based on the well-known system reliability theory and survival analysis. It is a binary point-like model as it considers the whole network as one node either available or unavailable. The estimated overall failure rate is mixing all failure modes, mechanisms, criticality, locality, and pipe functionalities. It overestimates the failures’ occurrence. Accordingly, the model is conservative as it overestimates the unavailability and underestimates the reliability.
The predictive model allows us to propose four indicators: the asymptotic unavailability u , the meantime before a failure is monitored ( MTBF ^ ), the expected number of failures n ^ T within an interval T , and the most probable interval T M P to observe n failure-occurrences. These four indicators are interested because they are field direct measures. So, they are proposed to assess the effectiveness of the decisions and actions taken within the modernisation and the extension of the regional WSS.
In addition, they can be of great interest as the modernisation of the WSS will imply the use of the ICT for in-line monitoring and control activities.
The robustness of the predictive model and its control indicators will be assessed once the databases are updated by the integration of the post-modernisation failure-recovery data.
The proposed assessment, through the implementation of the failure and recovery analysis and scaled models based on the proposed global model, can be very helpful in developing the water safety plans as recommended by WHO [40] and the assessment of water leakage levels according to the revised Drinking Water Directive [41].

Author Contributions

All authors equally contributed to the development of this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The Availability Model

This model is extensively mentioned in almost all textbooks and class notes relative to system reliability modelling and analysis.
The model is a binary point-like where the WSS is either available a t or unavailable u t at time t . So, one can express the unavailability of the WSS such as:
u t + Δ t = 1 μ · Δ t u t + λ · Δ t · a t ,
where:
  • u t : is the probability that the WSS is unavailable at t
  • u t + Δ t : is the probability that the WSS is unavailable at time t + Δ t . can be developed as u t + d d t u t · Δ t + ε Δ t Δ t 0 0
  • μ . Δ t : is the probability that the WSS does not recover within the interval Δ t
  • a t : is the probability that the WSS is available at t
  • λ . Δ t : is the probability that the WSS fails within the interval Δ t
  • Δ t : is a time interval as small as possible that tends to zero.
a t + u t = 1 ,   t 0 ,
The previous equation can be developed considering the 1st order terms and taking the limits when Δ t goes to zero. We then find:
d d t u t = μ · u t + λ . a t   = λ + μ · u t + λ ,
where; u t + a t = 1 at all instants of time t .
If and only if, the WSS characteristic transition rates μ   , λ are constant with the time, the overall unavailability can be determined by:
u t =   λ λ + μ 1 exp ( λ + μ t + u 0 · exp λ + μ t ,
where, u 0 is the unavailability initial condition.
Following the same demarche above, one can deduce the expression of the availability such as:
a t =   μ λ + μ 1 exp ( λ + μ t + a 0 · exp λ + μ t ,
The WSS availability and unavailability have asymptotic values ( a , u ) when t equal to:
u = u t t =   λ λ + μ ,
and,
a = u t t =   μ λ + μ ,

Appendix B

The Reliability Model

Using the same demarche as above, one may determine the reliability R T + Δ t over a time interval T + Δ T by:
R T + Δ T = 1 λ · Δ T R T ,
where:
  • R T : is the probability that the WSS is reliable over the interval T ,
  • λ . Δ T : is the probability that the WSS does not fail within the interval Δ T ,
  • Δ T : is a time interval as small as possible that tends to zero,
  • R T + Δ T : is the probability that the WSS is reliable over the interval T + Δ T .
It can be developed to R T + d d T R T · Δ T + ε Δ T Δ T 0 0 .
The previous equation can then be developed considering the 1st order terms and taking the limits when Δ t goes to zero. We then find:
d d T R T   = λ · R T ,
whose solution is:
R T = E x p 0 T λ τ d τ ,
If and only if the WSS overall failure rate is a time constant, the reliability can then be determined by:
R T = E x p λ T ,
This reliability model is estimated to be valid in the vicinity of the 18 year-interval with a WSS overall failure rate λ w s s = 18.8   m 1 .
The model allows us to estimate the meantime before failure MTBF ^ as:
MTBF ^ = T = 0 T · E x p λ W S S T · λ W S S d T = 1 λ W S S .

References

  1. Gheisi, A.; Naser, G. Water distribution system reliability under simultaneous multicomponent failure scenario. J. Am. Water Works Assoc. 2014, 7, 83–84. [Google Scholar] [CrossRef]
  2. Thornton, J. Managing leakage by managing pressure: A practical approach. Water 2003, 6, 43–44. [Google Scholar]
  3. Sahu, N.; Panda, A.; Nayak, S.; Saini, A.; Mishra, M.; Sayama, T.; Sahu, L.; Duan, W.; Avtar, R.; Behera, S. Impact of Indo-Pacific Climate Variability on High Streamflow Events in Mahanadi River Basin, India. Water 2020, 12, 1952. [Google Scholar] [CrossRef]
  4. Boryczko, K.; Rak, J. Method for assessment of water supply diversification. Resources 2020, 9, 87. [Google Scholar] [CrossRef]
  5. EN 15975-2:2013. Security of Drinking Water Supply. Guidelines for Risk and Crisis Management. Risk Management. Available online: https://standards.iteh.ai/catalog/standards/sist/dd2df50c-59ec-40f4-845a-00b83dfdd6df/sist-en-15975-2-2013 (accessed on 9 December 2020).
  6. D’Ercole, M.; Righetti, M.; Raspati, G.S.; Bertola, P.; Ugarelli, R.M. Rehabilitation planning of water distribution network through a reliability-based risk assessment. Water 2018, 10, 277. [Google Scholar] [CrossRef] [Green Version]
  7. Coudray, R.; Eid, M.; Fiorinia, G.L. A differential model for equivalent parameters Application to two residual heat removal systems parameterization and interpretation of results. Nucl. Eng. Des. 1997, 168, 227–239. [Google Scholar] [CrossRef]
  8. Sun, C.; Parellada, B.; Puig, V.; Cembrano, G. Leak Localization in Water Distribution Networks Using Pressure and Data-Driven Classifier Approach. Water 2020, 12, 54. [Google Scholar] [CrossRef] [Green Version]
  9. Ondrejka Harbulakova, V.; Estokova, A.; Kovalcikova, M. Correlation Analysis between Different Types of Corrosion of Concrete Containing Sulfate Resisting Cement. Environments 2017, 4, 44. [Google Scholar] [CrossRef]
  10. Robles-Velasco, A.; Cortés, P.; Muñuzuri, J.; Onieva, L. Prediction of pipe failures in water supply networks using logistic regression and support vector classification. Reliab. Eng. Syst. Saf. 2020, 196. [Google Scholar] [CrossRef]
  11. Zieja, M.; Wazny, M.; Stepien, S. Outline of a method for estimating the durability of components or device assemblies while maintaining the required reliability level. Eksploat. Niezawodn. Maint. Reliab. 2018, 20, 260–266. [Google Scholar] [CrossRef]
  12. Barton, N.A.; Farewell, T.S.; Hallett, S.H.; Acland, T.F. Improving pipe failure predictions: Factors affecting pipe failure in drinking water networks. Water Res. 2019, 164. [Google Scholar] [CrossRef] [PubMed]
  13. Urbanik, M.; Tchórzewska-Cieślak, B.; Pietrucha-Urbanik, K. Analysis of the Safety of Functioning Gas Pipelines in Terms of the Occurrence of Failures. Energies 2019, 12, 3228. [Google Scholar] [CrossRef] [Green Version]
  14. Arai, Y.; Koizumi, A.; Inakazu, T.; Watanabe, H.; Fujiwara, M. Study on failure rate analysis for water distribution pipelines. J. Water Supply Res. Technol. 2010, 59, 429–435. [Google Scholar] [CrossRef]
  15. Królikowska, J. Application of PHA method for assessing risk of failure on the example of sewage system in the city of Krakow. Rocz. Ochr. Sr. 2011, 13, 693–710. [Google Scholar]
  16. Tabesh, M.; Soltani, J.; Farmani, R.; Savic, D. Assessing pipe failure rate and mechanical reliability of water distribution networks using data-driven modeling. J. Hydroinform. 2009, 11, 1–17. [Google Scholar] [CrossRef]
  17. Pietrucha-Urbanik, K. Multidimensional comparative analysis of water infrastructures differentiation. In Proceedings of the Conference on Environmental Engineering IV, Lublin, Poland, 3–5 September 2012; pp. 29–34. [Google Scholar] [CrossRef]
  18. Szpak, D. Method for Determining the Probability of a Lack of Water Supply to Consumers. Energies 2020, 13, 5361. [Google Scholar] [CrossRef]
  19. Taeho, C.; Sewan, L.; Dooil, K.; Mincheol, K.; Jayong, K. Application of Management Reliability Index for Water Distribution System Assessment. Environ. Eng. Res. 2014, 19, 117–122. [Google Scholar] [CrossRef]
  20. Martinez-Codina, A.; Castillo, M.; Gonzales-Zeas, D.; Garrote, L. Pressure as a predictor of occurrence of pipe breaks in water distribution networks. Urban. Water J. 2016, 13, 676–686. [Google Scholar] [CrossRef]
  21. Rak, J. Bases of Water Supply System Safety; Polish Academy of Science: Lublin, Poland, 2005. [Google Scholar]
  22. Kwietniewski, M.; Roman, M.; Kłos-Trębaczkiewicz, H. Reliability of Water Supply and Sewage Systems; Arkady: Warsaw, Poland, 1993. [Google Scholar]
  23. StatSoft, Inc. STATISTICA (Data Analysis Software System), Version 12. Available online: www.statsoft.com (accessed on 15 July 2020).
  24. Box, G.E.P.; Jenkins, M.G.; Reinsel, G.C.; Ljung, G.M. Time Series Analysis: Forecasting and Control; Wiley: Hoboken, NJ, USA, 2015. [Google Scholar]
  25. Sokołowski, A. Forecasting and analysis of time series. In Training Materials; StatSoft: Cracow, Poland, 2016. [Google Scholar]
  26. Guizzi, G.; Silvestri, C.; Romano, E.; Revetria, R. A comparison of forecast models to predict weather parameters. Adv. Energy Environ. Sci. Eng. 2015, 2, 88–96. Available online: https://www.researchgate.net/publication/283449890_A_comparison_of_forecast_models_to_predict_weather_parameters (accessed on 20 July 2020).
  27. Młyński, D.; Bergel, T.; Młyńska, A.; Kudlik, K. A study of the water supply system failure in terms of the seasonality: Analysis by statistical approaches. AQUA Water Infrastruct. Ecosyst. Soc. 2021, 151. [Google Scholar] [CrossRef]
  28. Rahman, H.; Salma, U.; Hossain, M.; Khan, T.F. Revenue Forecasting Using Holt–Winters Exponential Smoothing. Res. Rev. J. Stat. 2017, 5, 19–25. Available online: https://www.researchgate.net/publication/311945797_Revenue_Forecasting_using_Holt-Winters_Exponential_Smoothing (accessed on 25 July 2020).
  29. AL-Washali, T.; Sharma, S.; AL-Nozaily, F.; Haidera, M.; Kennedy, M. Modelling the Leakage Rate and Reduction Using Minimum Night Flow Analysis in an Intermittent Supply System. Water 2019, 11, 48. [Google Scholar] [CrossRef] [Green Version]
  30. Pietrucha-Urbanik, K.; Tchórzewska-Cieślak, B.; Eid, M. Water Network-Failure Data Assessment. Energies 2020, 13, 2990. [Google Scholar] [CrossRef]
  31. Vališ, D.; Hasilová, K.; Forbelská, M.; Vintr, Z. Reliability modelling and analysis of water distribution network based on backpropagation recursive processes with real field data. Measurement 2020, 149. [Google Scholar] [CrossRef]
  32. Zimoch, I. Pressure control as part of risk management for a water-pipe network in service. Ochr. Srod. 2012, 34, 57–62. [Google Scholar]
  33. Gould, S.J.F.; Boulaire, F.A.; Burn, S.; Zhao, X.L.; Kodikara, J.K. Seasonal factors influencing the failure of buried water reticulation pipes. Water Sci. Technol. 2011, 11, 2692–2699. [Google Scholar] [CrossRef]
  34. Pietrucha-Urbanik, K.; Rak, J.R. Consumers’ Perceptions of the Supply of Tap Water in Crisis Situations. Energies 2020, 13, 3617. [Google Scholar] [CrossRef]
  35. Kutyłowska, M. Comparison of two types of artificial neural networks for predicting failure frequency of water conduits. Period. Polytech. Civ. 2017, 61, 1–6. [Google Scholar] [CrossRef] [Green Version]
  36. Rak, J.R.; Tchórzewska-Cieślak, B.; Pietrucha-Urbanik, K. A Hazard Assessment Method for Waterworks Systems Operating in Self-Government Units. Int. J. Environ. Res. Public Health 2019, 16, 767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Park, S.; Kim, K. Development of new computational methods for identifying segments and estimating the risk of water supply interruption for a segment in water pipe networks. Desalination Water Treat. 2017, 99, 211–219. [Google Scholar] [CrossRef]
  38. Debón, A.; Carrión, A.; Cabrera, E.; Solano, H. Comparing risk of failure models in water supply networks using ROC curves. Reliab. Eng. Syst. Saf. 2010, 95, 43–48. [Google Scholar] [CrossRef]
  39. Parka, A.; Kuliczkowska, E.; Kuliczkowski, A.; Zwierzchowska, A. Selection of pressure linings used for trenchless renovation of water pipelines. Tunn. Undergr. Space Technol. 2019, 98, 103218. [Google Scholar] [CrossRef]
  40. World Health Organization. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011; Available online: https://apps.who.int/iris/bitstream/handle/10665/44584/9789241548151_eng.pdf;jsessionid=FB26DE4E81767BC7525DC61A1537C754?sequence=1 (accessed on 15 December 2020).
  41. Directive (EU) 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the Quality of Water Intended for Human Consumption, OJ L 435, 23.12.2020. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32020L2184&from=EN (accessed on 20 December 2020).
Figure 1. The failure rate in each year of the analysis: (a) distribution pipes, (b) water connections, (c) mains, and (d) the whole water network.
Figure 1. The failure rate in each year of the analysis: (a) distribution pipes, (b) water connections, (c) mains, and (d) the whole water network.
Energies 14 03305 g001aEnergies 14 03305 g001b
Figure 2. Whisker Plot Box: failure rate averaged per month for: (a) distribution pipes, (b) water connections, (c) mains, and (d) the whole water network.
Figure 2. Whisker Plot Box: failure rate averaged per month for: (a) distribution pipes, (b) water connections, (c) mains, and (d) the whole water network.
Energies 14 03305 g002aEnergies 14 03305 g002b
Figure 3. Predicted failures of mains by exponential smoothing, for additive seasonality during months for S0 = 1.26, T0 = 0.0017: (a) data vs. model comparison, (b) model residuals.
Figure 3. Predicted failures of mains by exponential smoothing, for additive seasonality during months for S0 = 1.26, T0 = 0.0017: (a) data vs. model comparison, (b) model residuals.
Energies 14 03305 g003
Figure 4. Cumulative failures and temperatures below 0 °C.
Figure 4. Cumulative failures and temperatures below 0 °C.
Energies 14 03305 g004
Figure 5. The component-repair meantime T c , the organizational-repair meantime T o , and the recovery meantime T r , per year.
Figure 5. The component-repair meantime T c , the organizational-repair meantime T o , and the recovery meantime T r , per year.
Energies 14 03305 g005
Figure 6. The recovery meantime in each year of the analysis.
Figure 6. The recovery meantime in each year of the analysis.
Energies 14 03305 g006
Figure 7. Recovery meantime vs. water pipe materials.
Figure 7. Recovery meantime vs. water pipe materials.
Energies 14 03305 g007
Figure 8. Failure removal time regarding failure mode of water pipe.
Figure 8. Failure removal time regarding failure mode of water pipe.
Energies 14 03305 g008
Figure 9. The time evolution of the estimated unavailability u(t) of the WSS.
Figure 9. The time evolution of the estimated unavailability u(t) of the WSS.
Energies 14 03305 g009
Figure 10. The distribution P n T of having n failure-occurrences with the interval T , for n = 0 ,   1 ,   2 ,   3 .
Figure 10. The distribution P n T of having n failure-occurrences with the interval T , for n = 0 ,   1 ,   2 ,   3 .
Energies 14 03305 g010
Table 1. Synthesis of the statistical analyses of the failure data.
Table 1. Synthesis of the statistical analyses of the failure data.
Pipes TypeAverage
N o .   o f   F a i l u r e s   · k m 1 · a 1
Estimated Trend R 2 Relative Length
Distribution 2.87 × 10 1 λ d i s = 0.0012 t 2 0.0452 t + 0.5777 91 % ~60%
Connections 3.62 × 10 1 λ c o n = 7 × 10 6 t 5 + 4 × 10 4 t 4 0.007 t 3 + 0.0675 t 2 + 0.009 t + 0.0194 27.8 % ~35%
Mains 1.06 λ m a i = 3 × 10 5 t 5 0.0015 t 4 + 0.0265 t 3 0.2085 t 2 + 0.6351 t + 0.7024 49 % ~5%
Total 4.82 × 10 1 λ t o t = 0.0008 t 2 0.0309 t + 0.558 70.3 % ~100%
Table 2. The recovery times’ mean values.
Table 2. The recovery times’ mean values.
Parameter T c T o T r
Meantime—µ (h)2.047.799.83
σ0.450.410.27
σ/µ%22.04%5.27%2.78%
Upper limit2.738.4210.27
Lower limit1.427.159.39
Table 3. The most probable intervals to observe n failure-occurrences in the WSS.
Table 3. The most probable intervals to observe n failure-occurrences in the WSS.
Metric n = 1 n = 2 n = 3 n = 4
P n T T 1.6 n n ! e T 1.6 T 1.6 e T 1.6 T 1.6 2 2 e T 1.6 T 1.6 3 6 e T 1.6 T 1.6 4 24 e T 1.6
T MP ( day ) 1.6 n 1 × 1.6 2 × 1.6 3 × 1.6 4 × 1.6
P n T MP n n n ! e n 0.368 0.271 0.2240.195
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pietrucha-Urbanik, K.; Tchórzewska-Cieślak, B.; Eid, M. A Case Study in View of Developing Predictive Models for Water Supply System Management. Energies 2021, 14, 3305. https://0-doi-org.brum.beds.ac.uk/10.3390/en14113305

AMA Style

Pietrucha-Urbanik K, Tchórzewska-Cieślak B, Eid M. A Case Study in View of Developing Predictive Models for Water Supply System Management. Energies. 2021; 14(11):3305. https://0-doi-org.brum.beds.ac.uk/10.3390/en14113305

Chicago/Turabian Style

Pietrucha-Urbanik, Katarzyna, Barbara Tchórzewska-Cieślak, and Mohamed Eid. 2021. "A Case Study in View of Developing Predictive Models for Water Supply System Management" Energies 14, no. 11: 3305. https://0-doi-org.brum.beds.ac.uk/10.3390/en14113305

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