Next Article in Journal
Using Modeling Tools to Better Understand Permafrost Hydrology
Previous Article in Journal
Impact of Institutional Change on Irrigation Management: A Case Study from Southern Uzbekistan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Models for the Peak Residential Water Demand

Dipartimento di Ingegneria Civile e Meccanica, University of Cassino and Southern Lazio, Via G. Di Biasio, 43, 03043 Cassino (FR), Italy
*
Author to whom correspondence should be addressed.
Water 2017, 9(6), 417; https://doi.org/10.3390/w9060417
Submission received: 31 December 2016 / Revised: 5 June 2017 / Accepted: 7 June 2017 / Published: 10 June 2017

Abstract

:
Peak water demand is one of the most stringent operative conditions for a Water Distribution System (WDS), not only for the intensity of the event itself, but also for its recurring nature. The estimation of the maximum water demand is a crucial aspect in both the design and management processes. Studies in the past have tackled this issue with deterministic approaches, even if peak phenomena are distinctly random. In this work, probabilistic models have been developed to study and forecast the daily maximum residential water demand. Some probability distributions have been tested by means of statistical inferences on different data samples related to three monitored WDS. The parameter estimations of the proposed equations have been related to the number of supplied users. Furthermore, this work investigates time scaling effects on the effectiveness of the proposed distributions and relations. Corrective factors that take into account the effect of time averaging step on the above-mentioned parameters have been proposed.

1. Introduction

The peak water demand for residential users implies one of the most onerous operative conditions for an urban Water Distribution System (WDS). The maximum water demand during the day, along with mechanical failure of one or more components of the system (e.g., pipes, pumps, and valves) or fire-fighting flow conditions, is one of the scenarios for which the WDSs may experience some deficiency.
Hence, the peak water demand is a working condition that is usually taken into account for the design and management of WDSs (e.g., [1,2]).
In the scientific literature, few papers define completely the daily maximum water demand for a considerable number of users (around one thousand inhabitants) following a probabilistic approach. This issue is often tackled by means of the more practical deterministic equations (e.g., [3,4,5]), although the random nature of the water demand is generally recognized.
On the other hand, several stochastic models describing the water demand for few end users (a dozen of dwellings) were proposed (e.g., [6,7,8,9]). These approaches model the residential water demand as it happens at the hydro-sanitary devices of the dwellings. Therefore, the demanded flow is generally composed by the overlapping of several rectangular pulses, which are generated by turning on the hydro-sanitary devices.
The rectangular pulses approaches need to model at least three random phenomena [6]: arrival, duration and intensity of the demand pulses. Hence, these models induce significant computational efforts and require a detailed knowledge of the lifestyle related to the users on which the parameters of the component distributions depend.
However, the rectangular pulse models could be overwhelming when the target of the study is the estimation of the water requirement in a specific condition, such as the maximum water demand during the day. Examples of these approaches can be found in [10,11] where a very detailed knowledge of the dwelling composition and user habits is required.
On the other hand, when the number of the users is greater than 200, the trend of the water demand during the day becomes smooth and the single pulses are no longer distinguishable [12,13]. In these conditions, the water demand can be represented by means of a single random variable that is continuous and positive.
The aim of this work is to obtain a reliable probabilistic approach that is capable of providing accurate estimations of the maximum residential water demand. For this goal, a specific field laboratory, which monitors the actual water demand of a small town in Southern Italy, was set up by the Laboratorio di Ingegneria delle Acque (LIA) of the University of Cassino and Southern Lazio. This monitoring system allowed obtaining copious data samples, where the number of the served inhabitants is known. The data sample was increased considering two further case studies whose users have very different lifestyles in respect to the monitored inhabitants of the Southern Lazio. In particular, the peak water demand for two further towns was also taken into account, one in the Netherlands (the case study of Franeker [14]) and one in Northern Italy (the case study of Castelfranco Emilia [7]).
The enlargement of the data sample showed the independence of the daily maximum demand from the local habits. Hence, the results of this study could be considered of great practical application, allowing the use of the hereby proposed probabilistic models to obtain a good estimation of the residential peak water demand.
Some effective distributions, able to model the daily maximum demand for residential users, have been proposed by means of statistical inferences on the above-mentioned observed data. In addition, practical and reliable equations have been suggested which allow estimating the parameters of the suggested Cumulative Distribution Functions (CDFs) in relation to the number of users, Nus, where Nus represents the total number of the dwelling occupants.
The temporal resolution is a fundamental choice when modeling water demand (e.g., [15,16,17,18,19]), leading to significant reductions of the peak demand when the time step increases [20,21]. Therefore, time scaling effects were also investigated and further equations that take into account such phenomenon were presented. The proposed equations allow estimating the CDFs’ parameters when the time step (∆t) varies in the interval 1 min ≤ ∆t ≤ 1 h.

2. The monitoring Systems and the Data Samples

In real WDSs monitoring systems, flow measures do not often represent the actual users demand, generally for two main issues.
The hydraulic shortages (first issue), which could be caused in different points of a WDS for under-dimensioned pipes or failure of electro-mechanical components, do not allow considering the flow measurements as an expression of the water demand. Indeed, in these working conditions, the measurements represent only the water delivery capability of the WDS under specific user requests; hence, the obtained data underestimate the actual water demand. In addition, these shortage conditions are likely when the daily maximum water demand occurs.
The monitored WDSs have to be redundant in order to avoid that the collected data are compromised by such operating conditions. In this way, the WDSs are able to totally satisfy the users’ water demand, even when this is particularly high. Hence, the monitoring system must be composed by both flow meters and pressure probes, to check that demand data obtained is reliable.
A fundamental parameter to describe the peak phenomenon is Nus (e.g., [5,22]), but for looped WSDs, when the measurements are performed on the links, it is not possible to know univocally the number of users (second issue) related to each measurement point. The served users can vary with the different working conditions.
In order to solve this issue, the measuring points have to be set on the pipes which connect the Demand Monitoring Areas (DMAs) to the rest of the WDS. A DMA is a part of the entire distribution network which is linked to the rest of the WDS by means of one or a limited number of pipes. In this way, Nus for each DMA is not dependent on the working conditions.
An accurate census must define the number of the resident population in the considered DMA, and it must verify that the users are only residential.
Based on these considerations a specific monitoring system, which involves the real WDS of a small town, Piedimonte San Germano (PSG), was realized.
In detail, the monitoring system is composed by four measurement points, each with an electromagnetic flow meter and a pressure cell, where all probes are connected to a data logger. Flows and pressures were recorded continuously, where a variable frequency data logger allowed an acquisition frequency up to 1 Hz. The location of the probes allowed the measurement of the water demand for four DMAs, where the number of users was respectively equal to: 239, 777, 981 and 1220. More details of the PSG laboratory are provided in [21,23].
In order to expand the data set, the time series of two further case studies, with similar characteristics to PSG in terms of data reliability, but characterized by a different users lifestyle, were also used: Castelfranco Emilia (CE) in Italy and Franeker (Fr) in the Netherlands.
CE monitoring system covering 596 users is part of a wider WDS. Data are collected by means of volumetric measures (Woltmann volume meter), and they are averaged over one-minute intervals. For Fr network, continuous per minute flow measurements are taken on the overall network, which supplies about 1150 users.
Alvisi et al. [7] and Blokker et al. [14] provide respectively detailed information of the CE and Fr field laboratories.
Table 1 summarizes the principal characteristic of the three in-situ laboratories.
Initially, the analysis assumed a time resolution of one minute, which was increased up to ∆t =1 h in the following elaborations in order to investigate the effects of the temporal aggregation.
The water demand presents different trends between working days and weekends which affect also the maximum water demand during the day. Therefore, only the working days have been considered in which the peak water demand is slightly more pronounced.
No significant seasonal trend of the water demands was observed for the field data examined, contrary to what happens for residential users in some countries [24,25].
On this basis, in order to obtain a more reliable data set, the time series of the three field laboratories were filtered. The daily series were removed when they presented one of the following conditions: holidays/weekends; days with missing data; days with anomalous data (e.g., fire flow, corrupted signals of the probes, and big leakages). The number of days used during the following elaborations is indicated in the last row of Table 1.
In the following sections, the maximum demand of flow during the day has been made dimensionless by means of the ratio:
C p = Q p μ Q
where μQ is the daily mean water demand, and Qp the peak flow during the day. Equation (1) gives the Peak Demand Coefficient, Cp. This work focuses on the dimensionless variable Cp.

3. Probabilistic Distributions for the Cp

The extreme events are effectively described by means of the Log-Normal distribution (LN) (e.g., [26])
f x L N = 1 x σ y 2 π exp [ 1 2 ( ln x μ y σ y ) 2 ]
and Gumbel distribution (Gu) (e.g., [26])
f x G u = α exp [ α ( x ε ) exp [ α ( x ε ) ] ]
where f x is the probability density function; x is the random variable and y = lnx; and μ and σ represent, respectively, the mean and the standard deviation. The parameters α and ε of Gu are given by the relations:
1 α = 6 π σ x ε = μ x 0.45 σ x
In the scientific literature, the peak water demand for residential users was described by means of these probabilistic models (e.g., [20,21,27]).
The LN and Gu have been tested by means of the observed data for the above described case studies, assuming a time aggregation ∆t = 1 min. In addition, the Log-Logistic (LL) model was investigated to describe the daily maximum water demand. The LL distribution presents a trend very similar to the LN, but, contrary to Gauss model, the probability density function can be analytically integrated [28,29]. The CDF of the LL distribution has the following expression:
F x L L = [ ( e μ y x ) π σ y 3 + 1 ] 1
If x = Cp (dimensionless peak demand Equation (1)), y is lnCp.
The parameters of Equation (5) can be estimated by means of statistics on the logarithmic of the original sample, or by means of μ and σ of the pristine random variable by means of the following equations [28]:
μ y = l n [ μ x σ y 3 s i n ( σ y 3 ) ] CV x = 1 σ y 3 t a n ( σ y 3 ) 1
where CVx = σxx is the variation coefficient.
However, an early and rough estimation of the standard deviation of y variable can be obtained by means of the assumption σy CVx (this assumption becomes absolutely robust when CVx is less than 0.15).
The diagrams in Figure 1 are the Quantile-Quantile (Q-Q) plots for three monitored number of users (for the sake of synthesis, the plots of PSG Nus = 239, 777 and 1220 were not here reported, but they are similar to diagrams of Figure 1). More precisely, the Q-Q plots show the comparison between the observed and the theoretical Cp quantiles with the same probability of occurrence. The theoretical values were estimated by means of the integral of Equations (2) and (3)—numerical integration for LN distribution—and Equation (5). The quantiles of the peak demand of Figure 1 fit well the bisector line, which represents the condition where the observed quantiles are equal to the theoretical quantiles.
The effectiveness of the three distributions was checked by means of the Kolmogorof–Smirnov (KS) test. Table 2 summarizes the threshold values for a 2% confidence level (D2%—column 3) and the relative KS parameters (DKS) for the LN Gu, and LL CDFs (columns 4–6). The goodness of fit tests were satisfied (DKS < D2%) for all monitored users. In addition, Table 2 shows the statistics of the data sample (µCp and σCp—columns 7–8). These probability distributions proved to be equally effective to model the peak phenomena, thus the choice is left to the decision maker.
Equations (2) and (3) allow to explicate the Peak Demand Coefficient [ C p ] F . for a predefined probability Pr[S] of not exceedance. Hence the peak coefficient with a predefined probability of not exceedance following LL is given by:
[ C p ] F = e μ y [ Pr [ S ] 1 Pr [ S ] ] σ y 3 π
while, for Gu distribution, [ C p ] F presents the following relation:
[ C p ] F = σ x [ 1 CV x 0.45 6 π ln ( ln 1 Pr [ S ] ) ]
The quantile [ C p ] F for LN model can be estimated only by means of a numerical approximation.
Figure 2 shows—as an example—the comparison between the observed data (in this case the peak demand of PSG for Nus = 1220 and Δt = 1 min) and the values of [ C p ] F estimated by means of Equations (7) and (8), assuming µx = 2.50 and σx = 0.22. The plots confirm that both the LL and Gu distributions fit well the experimental points.

4. Estimation of Parameters

All the proposed distributions (LN, Gu and LL) are bi-parametric models; hence, they require the estimation of two parameters, for instance, the mean and the variation coefficient.
The PSG demand data, jointly with data from the other monitored users (CE and Fr), shows that the dimensionless average of the peak water demand varies in relation to the number of users according to the following equation (Figure 3a):
μ C P = 10 N u s 0.2
Equation (9) resembles the relation proposed by Babbitt [30], but the coefficient 10 of the power function is significantly less than the one of the Babbitt equation, as already shown in Tricarico et al. [21]. The deterministic approaches, of which the Babbitt equation represents a classic example, imply an overestimation of the peak demand.
CV C p for peak demand can also be estimated in relation to the number of users.
CV C p = 0.1 + 20 N u s
Equation (10) fits the experimental points (Figure 3b), and it highlights a decreasing trend of the CV C p which tends asymptotically towards 0.1 when the number of the users increases. Therefore it is reasonable to assume CV C p = 0.1 when the number of residential inhabitants is greater than 1000 [21,31].
It is worth noting that the investigated distributions and the suggested relations–Equations (9) and (10)—are effective for users with completely different habits. Gargano et al. [13] give more details about the daily pattern of the water demand for Ce, Fr and PSG users.
Equations (9) and (10) adjust the relations introduced by Tricarico et al. [21], because the hereby proposed equations have been obtained on the basis of more copious and detailed data sample. In addition, the data of water demand recorded in five suburbs of Melbourne in Australia [32] show similar trends to the ones of Equations (9) and (10), respectively, for the mean and coefficient of variation.
For a predefined number of users, the parameters ( μ C P and CV C p ) can be estimated by means of the Equations (9) and (10), and hence Equations (7) and (8) allow defining the Peak Demand Coefficient with assigned probabilities of not exceedance.
The diagrams of Figure 4 show that higher values of the probability of not exceeding of Cp lead to greater peak demand coefficients, and this effect is more evident when the number of users decreases. Moreover, the comparison of Figure 4a,b denotes a close trend between Gu and LL distributions for equal values of Nus and fixed probability.
Finally, for the investigated number of users (200–1250), Figure 4 shows that the old deterministic equations (i.e., [30]) overestimate the peak phenomena with respect to the probabilistic relations, also when these latter assume high values of the success probability (i.e., Pr [ C p C ¯ p ] = Pr [ S ] = 0.999 ).
It is worth noting that the classic relation of Babbitt [30] has been recalled here only as an example of deterministic relations, in fact all of them usually lead to an overestimation of the peak phenomenon of the water demand [31]. Several researchers [3,10,22] give an interesting overview of the ancient deterministic equations for the estimation of the peak phenomena.
The Peak Coefficient curve (Figure 4) for a defined probability of non-exceedance can also be expressed by means of the product:
[ C p ] F = K F · μ C p
where μ C p is estimated by means of Equation (9) in function of the users number, and KF represents a growth factor (KF > 1), as used for flood frequency analysis (e.g., [33]); therefore, it increases when the Pr[S] increases. For LL and Gu, the growth factors are obtained, respectively, from Equations (7) and (8):
K F = [ Pr [ S ] 1 Pr [ S ] ] σ y 3 π s i n ( σ y 3 ) σ y 3 L L distribution
K F =   1 0.45 CV x ( 6 π ln ( ln 1 Pr [ S ] ) ) CV x G u distribution
KF is estimated when the probability of not exceedance Pr[S] is fixed (Equations (12) and (13)) and the probabilistic parameters are defined in relation to the number of the inhabitants.
The probability of non-exceedance usually is an engineering choice made by the decision makers and should be a consequence of a cost–benefit analysis.

5. Time Aggregation Effects

The results shown in the previous sections are related to a time step equal to ∆t = 1 min, while, in the numerical analysis of WDSs, a less refined time resolution is usually assumed (1 min < ∆t < 60 min).
The time step presents noticeable implications on the estimation of the maximum water demand [3,27], in fact when ∆t increases, the time averaging effect smooths the peak phenomena.
In the following sections, the effects of the time aggregations are investigated in order to extend the above proposed distributions and relations to other temporal averaging intervals.

5.1. Distributions in Relation to Time Step

The same data samples of the residential water demand, which were considered in previous sections for ∆t = 1 min, were resampled to define the distributions with a different time aggregation (∆t = 5, 10, 15, 30 and 60 min). The effectiveness of the LN, Gu and LL distributions was demonstrated by means of statistical inferences for different time resolutions, from 1 min up to 1 h.
For the sake of the synthesis, in Figure 5, only the diagrams for the Cp of PSG 981 with different time resolutions (∆t = 1, 10, 30 and 60 min), are reported. However, similar results have been obtained for all other monitored users with a time resolution up to one hour.
The Q-Q plots of Figure 5 show the goodness of the three distributions for all investigated time aggregations, which has been also demonstrated more rigorously by means of KS tests.
The tests with significance level 2% (the value of D2% is represented by the dashed line in Figure 6a,c,e) were all satisfied for LN, Gu and LL distributions for any value of ∆t that falls into the interval 1–60 min.
The results of the KS tests did not highlight any trends in relation to the time step, hence the effectiveness the proposed distributions seem to be independent of the time resolution. In addition to KS test, the Mean Squared Error (MSE) was investigated, where MSE is given by:
M S E = 1 n i = 1 n ( x i x i ^ ) 2
where x i and x i ^ are, respectively, the observed and the predicted Cp values, and n is total number of observations.
The MSE plots (Figure 6b,d,f) show that the time aggregation for the investigated range seems not to trigger any tendency on the effectiveness of the three probability distributions. Therefore, the estimated values of MSE vary independently by Δt and Nus.

5.2. Time Dependence of the Parameters

If the time scale does not lead to significant effects on the probabilistic models, it has to be considered in order to estimate the parameters when ∆t changes. Indeed, the statistical analysis highlighted that the proposed equations to estimate the parameters of distributions need a reduction coefficient when the time step is greater than 1 min. Therefore, Equation (9) becomes:
μ C p ( t ) = k μ C p
where μ C p ( t ) is the average value of Cp for ∆t ≥ 1 min and k represents the corrective coefficient of the mean peak demand coefficient for ∆t = 1 min ( μ C p ), which decreases when ∆t increases (Figure 7a).
Hence, k is estimated in relation to the users’ number and ∆t
k = 1 + ( 0.04 ln N u s 0.35 ) ln ( t )
where the time is counted in min. Equations (17) and (18) are valid for 1 min ≤ ∆t ≤ 60 min and 250 users < Nus < 1250 users.
The abacus of Figure 7b shows that the peak demand coefficient decreases significantly when ∆t increases from 1 min to 60 min. For instance, it decreases of about 25% for Nus = 1250, and about 50% for Nus = 250.
Likewise, Equation (10) requires a corrective coefficient β in order to estimate CV C p ( t ) when the time aggregation is greater than 1 min:
CV C p ( t ) = 0.1 + 20 β N u s
The coefficient depends on the Nus and the time step according to Equation (18)
β = e x p [ ( t 1 ) N u s 10 4 ] .
where the time is counted in minutes and Equations (17) and (18) are valid for 1 min ≤ ∆t ≤ 60 min and 250 users < Nus < 1250 users.
Figure 8 shows the effectiveness of Equations (17) and (18) to describe the effects of the time aggregation on the Cp variation coefficient. Moreover, Figure 8 confirms what has been shown in Figure 3b; in fact CV C p = 0.1 seems to be a good simplification when the number of users is greater than 250 [31,34].
Equations (9) and (10), combined with Equations (15)–(18), allow estimating the parameters of the suggested distributions, for every value of the time aggregation and users number, when 1 min ≤ ∆t ≤ 60 min and 200 users < Nus < 1250 users.

6. Conclusions

The peak water demand is a very demanding working condition in WDSs and represents an issue of great interest for both design and management problems. The results of the present research give a practical contribution to the characterization of the daily maximum water demand according to a probabilistic approach. The proposed method allows estimating the peak demand through biparametric probability distributions and relations with practical applicability.
By means of statistical inferences which took into account copious data samples related to different residential users, the effectiveness of the Log-Normal and Gumbel distributions was verified to describe the peak water demand during the day. In addition, also the Log-Logistic distribution, whose CDF presents a closed form (as Gumbel distribution), is capable of modeling the maximum daily water demand.
Therefore, some practical equations were suggested, which allow estimating the daily peak demand for predefined values of the success probability.
The investigated case studies are very significant, not for the geographic location, but for the very different users’ behavior. In fact, the three towns present daily demand patterns whose peaks are different in number, in time of occurrence and relative entity (e.g., the ratio of the morning peak in respect to the evening one). The probabilistic models proved to be independent of the daily demand pattern and of the peak occurrence time. In fact, although the lifestyles of the investigated users are completely different, the proposed models fitted very well all datasets.
In addition, statistical analysis allowed obtaining effective equations to estimate the parameters of the proposed distributions in relation to the users’ number. Even if the investigated range of the aggregated users (200–1250) is limited, it is effective to model the water demand at the nodes of a WDS, because in numerical simulations the number of served users often falls in this interval.
T time resolution has been analyzed investigating the effects of the time step, in the range 1 min–1 h, on the peak phenomena. It was observed that an aggregation of 1 min ≤ ∆t ≤ 1 h does not affect the choice of the distribution. In fact, the three distributions (LN, Gu and LL) proved to be all equally effective when modeling the peak water demand for all the investigated ∆t. On the other hand, when estimating the parameters of the distributions, corrective coefficients are needed and some effective equations have been proposed for this purpose. These latter allow estimating the parameters of the CDFs in relation to the number of users and the value of ∆t. The time resolution analysis considered the range 1 min–1 h, because it covers a wide spectrum of engineering problems (e.g., design and management problems).
In the future, it is desirable that the results herein obtained could be made more robust increasing the data sample, so that the validity field of the models could also be increased. For this end, new in-situ laboratories, whose characteristics are consistent with those described in the work, are necessary.

Acknowledgments

The Authors thank Marco Franchini, Stefano Alvisi and Mirjam Blokker for providing water demand data for Castelfranco Emilia and Franeker monitoring systems. In addition, the Authors thank the anonymous Reviewers for their valuable comments.

Author Contributions

Rudy Gargano conceived the research, and all Authors contributed to the work.

Conflicts of Interest

The authors declare no conflict of interest.

List of the Symbols and Acronyms

Cp[-] Peak demand coefficient
CDF(acronym) Cumulative Distribution Function
CE(acronym) Castelfranco Emilia (one of three case studies)
CV[-] Variation coefficient
CV C p [-] Variation coefficient for peak demand coefficient when ∆t =1 min
CV C p ( t ) [-] Variation coefficient for 1 min < ∆t ≤ 60 min
A(acronym) Demand Monitoring Area
f x [depending on variable] Probability density function
Fx[-] CDF of random variable
Fr(acronym) Franeker (one of three case studies)
Gu(acronym) Gumbel distribution
KF[-] growth factor
k [-] corrective coefficient of μ C p ( t ) , valid 1 min < ∆t ≤ 60 min
KS(acronym) Kolmogorov Smirnov test
LL(acronym) Log-Logistic distribution
LN(acronym) Log-Normal distribution
MSE[depending on variable] Mean Squared Error
n[-] total data number
Nus[-] Number of users
PSG(acronym) Piedimonte San Germano (one of three case studies)
Qp[L3/T] Peak flow of the demanded water during the day
WDS(acronym) Water Distribution System
x[depending on variable] Random variable
y[depending on variable] lnx
α[depending on variable] Gumbel parameter
β[-] valid 1 min < ∆t ≤ 60 min, valid 1 min < ∆t ≤ 60 min
t[T] Time step
ε[depending on variable] Gumbel parameter, mode
μ C p [-] mean peak demand coefficient for ∆t = 1 min
μ C p ( t ) [-] mean peak demand coefficient for 1 min < ∆t ≤ 60 min
μ Q [L3/T] mean daily water demand
σ[depending on variable] standard deviation

References

  1. Ysusi, M.A. System design: An overview. In Water Distribution Systems Handbook; Chapter 3; Mays, L.W., Ed.; McGraw-Hill: New York, NY, USA, 2000; pp. 1–49. ISBN 0-07-134213-3. [Google Scholar]
  2. Tricarico, C.; Morley, M.S.; Gargano, R.; Kapelan, Z.; de Marinis, G.; Savić, D.A. The influence of the existing network layout on water distribution system redesign analysis. J. Hydroinform. 2014, 16, 1375–1389. [Google Scholar] [CrossRef]
  3. Balacco, G.; Carbonara, A.; Gioia, A.; Iacobellis, V.; Piccinni, A.F. Evaluation of Peak Water Demand Factors in Puglia (Southern Italy). Water 2017, 9, 96. [Google Scholar] [CrossRef]
  4. Lingireddy, S.; Wood, D.J.; Nelson, A. Modified Pipe Network Model for Incorporating Peak Demand Requirements. J. Water Resour. Plan. Manag. 1998, 124, 296–299. [Google Scholar] [CrossRef]
  5. Martínez-Solano, J.; Iglesias-Rey, P.; Pérez-García, R.; López-Jiménez, P. Hydraulic analysis of peak demand in looped water distribution networks. J. Water Resour. Plan. Manag. 2008, 134, 504–510. [Google Scholar] [CrossRef]
  6. Buchberger, S.G.; Wu, L. A model for instantaneous residential water demands. J. Hydraul. Eng. 1995, 121, 232–246. [Google Scholar] [CrossRef]
  7. Alvisi, S.; Franchini, M.; Marinelli, A. A stochastic model for representing drinking water demand at residential level. Water Res. Manag. 2003, 17, 197–222. [Google Scholar] [CrossRef]
  8. Blokker, E.J.M.; Vreeburg, J.H.G.; van Dijk, J.C. Simulating residential water demand with a stochastic end-use model. J. Water Resour. Plan. Manag. 2010, 136, 19–26. [Google Scholar] [CrossRef]
  9. Gargano, R.; Di Palma, F.; de Marinis, G.; Granata, F.; Greco, R. A stochastic approach for the water demand of residential end users. Urban Water J. 2016, 13, 569–582. [Google Scholar] [CrossRef]
  10. Buchberger, S.; Blokker, M.; Cole, D.P. Estimating peak water demands in hydraulic systems I—Current practice. In Proceedings of the WDSA 2012: 14th Water Distribution Systems Analysis Conference, Adelaide, Australia, 24–27 September 2012; Engineers Australia: Barton, Australia, 2012; p. 1165. [Google Scholar]
  11. Blokker, M.; Vloerbergh, I.; Buchberger, S. Estimating peak water demands in hydraulic systems II—Future trends. In Proceedings of the WDSA 2012: 14th Water Distribution Systems Analysis Conference, Adelaide, Australia, 24–27 September 2012; Engineers Australia: Barton, Australia, 2012; p. 1165. [Google Scholar]
  12. Jankovic-Nisic, B.; Maksimovic, C.; Butler, D.; Graham, N.J.D. Use of Flow Meters for Managing Water Supply Networks. J. Water Resour. Plan. Manag. 2004, 130, 171–179. [Google Scholar] [CrossRef]
  13. Gargano, R.; Tricarico, C.; Del Giudice, G.; Granata, F. A Stochastic Model for Daily Residential Water Demand. Water Sci. Technol. Water Supply 2016, 16, 1753–1767. [Google Scholar] [CrossRef]
  14. Blokker, E.J.M.; Vreeburg, J.H.G.; Vogelaar, A.J. Combining the probabilistic demand model SIMDEUM with a network model. In Proceedings of the 8th Annual International Symposium on Water Distribution Systems Analysis, Cincinnati, OH, USA, 27–30 August 2006; p. 40. [Google Scholar] [CrossRef]
  15. Filion, Y.; Adams, B.; Karney, B. Cross correlation of demands in water distribution network design. J. Water Resour. Plan. Manag. 2007, 133, 137–144. [Google Scholar] [CrossRef]
  16. Moughton, L.J.; Buchberger, S.G.; Boccelli, D.L.; Filion, Y.R.; Karney, B.W. Effect of time step and data aggregation on cross correlation of residential demands. Proceeding of the 8th Annual International Symposium on Water Distribution Systems Analysis, Cincinnati, OH, USA, 27–30 August 2006; p. 42. [Google Scholar] [CrossRef]
  17. Magini, R.; Pallavicini, I.; Guercio, R. Spatial and temporal scaling properties of water demand. J. Water Resour. Plan. Manag. 2008, 134, 276–284. [Google Scholar] [CrossRef]
  18. Alvisi, S.; Ansaloni, N.; Franchini, M. Generation of Synthetic Water Demand Time Series at Different Temporal and Spatial Aggregation Levels. Urban Water J. 2014, 11, 297–310. [Google Scholar] [CrossRef]
  19. Vertommen, I.; Magini, R.; Cunha, M.C. Scaling Water Consumption Statistics. J. Water Resour. Plan. Manag. 2015, 141, 04014072. [Google Scholar] [CrossRef]
  20. Buchberger, S.G.; Nadimpalli, G. Leak estimation in water distribution systems by statistical analysis of flow readings. J. Water Resour. Plan. Manag. 2004, 130, 321–329. [Google Scholar] [CrossRef]
  21. Tricarico, C.; de Marinis, G.; Gargano, R.; Leopardi, A. Peak residential water demand. Water Manag. J. 2007, 160, 115–121. [Google Scholar] [CrossRef]
  22. Zhang, X.; Buchberger, S.G.; van Zyl, J.E. A theoretical explanation for peaking factors. In Proceedings of the ASCE World Water & Environmental Resources Congress, Anchorage, AK, USA, 15–19 May 2005. [Google Scholar]
  23. De Marinis, G.; Gargano, R.; Tricarico, C. Water demand models for a small number of users, ASCE. Proceeding of 8th Annual International Symposium on Water Distribution Systems Analysis, Cincinnati, OH, USA, 2007; p. 41. [Google Scholar] [CrossRef]
  24. Gato-Trinidad, S.; Jayasuriya, N.; Roberts, P. Forecasting residential water demand: Case study. J. Water Resour. Plan. Manag. 2007, 133, 309–319. [Google Scholar] [CrossRef]
  25. Rathnayaka, K.; Malano, H.; Maheepala, S.; George, B.; Nawarathna, B.; Arora, M.; Roberts, P. Seasonal demand dynamics of residential water end-uses. Water 2015, 7, 202–216. [Google Scholar] [CrossRef]
  26. Kottegoda, N.T.; Rosso, R. Applied Statistics for Civil and Environmental Engineers, 2nd ed.; Blackwell Publishing: Oxford, UK, 2008; p. 718. ISBN 978-1-4051-7917-1. [Google Scholar]
  27. Pallavicini, I.; Magini, R. Experimental analysis of residential water demand data: Probabilistic estimation of peak coefficients at small time scales. In Water Management Challenges in Global Change; Ulanicki, B., Vairavamoorthy, K., Butler, D., Bounds, P.L.M., Memon, F.A., Eds.; Taylor & Francis Group: London, UK, 2007; ISBN 978-0-415-45415-5. [Google Scholar]
  28. Swamee, P.K. Near lognormal distribution. J. Hydrol. Eng. 2002, 7, 441–444. [Google Scholar] [CrossRef]
  29. Ashkar, F.; Aucoin, F. Discriminating between the Lognormal and the Log-Logistic Distributions for Hydrological Frequency Analysis. J. Hydrol. Eng. 2012, 17, 160–167. [Google Scholar] [CrossRef]
  30. Babbitt, H.E. Sewerage and Sewage Treatment, 3rd ed.; Wiley: New York, NY, USA, 1928; pp. 28–33. [Google Scholar]
  31. Gargano, R.; Tricarico, C.; de Marinis, G. Residential water demand—Daily trends. In Proceedings of the Annual International Symposium on Water Distribution Systems Analysis, Tucson, AZ, USA, 12–15 September 2010; pp. 1314–1323. [Google Scholar] [CrossRef]
  32. Gato-Trinidad, S.; Gan, K. Characterizing maximum residential water demand. Urban Water 2012, 122, 15–24. [Google Scholar]
  33. Chow, V.T. The Log-Probability Law and its Engineering Applications; American Society of Civil Engineers: St. Louis, MO, USA, 1954; Volume 80, p. 25. [Google Scholar]
  34. Bao, Y.; Mays, L.W. Model for water distribution system reliability. J. Hydraul. Eng. 1990, 116, 1119–1137. [Google Scholar] [CrossRef]
Figure 1. Q-Q plots of the peak demand coefficient for: ∆t = 1 min; three distributions (LN, Gu and LL); and different users (Nus = 1150, 981 and 596).
Figure 1. Q-Q plots of the peak demand coefficient for: ∆t = 1 min; three distributions (LN, Gu and LL); and different users (Nus = 1150, 981 and 596).
Water 09 00417 g001
Figure 2. Peak water demand coefficient [Cp]F versus Pr[S].
Figure 2. Peak water demand coefficient [Cp]F versus Pr[S].
Water 09 00417 g002
Figure 3. Comparison of the proposed equations and of the observed data for the estimation of the: mean (a); and variation coefficients (b) of Cp versus Nus.
Figure 3. Comparison of the proposed equations and of the observed data for the estimation of the: mean (a); and variation coefficients (b) of Cp versus Nus.
Water 09 00417 g003
Figure 4. Peak Demand Coefficient versus user number for predefined success probability, and for CDFs: Log-Logistic Distribution (a); Gumbel Distribution (b).
Figure 4. Peak Demand Coefficient versus user number for predefined success probability, and for CDFs: Log-Logistic Distribution (a); Gumbel Distribution (b).
Water 09 00417 g004
Figure 5. Q-Q plots of the peak demand coefficient for Nus = 981 (PSG) for various time aggregations and distributions (LN, Gu and LL).
Figure 5. Q-Q plots of the peak demand coefficient for Nus = 981 (PSG) for various time aggregations and distributions (LN, Gu and LL).
Water 09 00417 g005
Figure 6. Effects of the time aggregation—KS parameter and MSE versus ∆t.
Figure 6. Effects of the time aggregation—KS parameter and MSE versus ∆t.
Water 09 00417 g006
Figure 7. (a) Comparison of the observed data and Equations (15) and (16); and (b) abacus of the corrective coefficient k for different values of Nus.
Figure 7. (a) Comparison of the observed data and Equations (15) and (16); and (b) abacus of the corrective coefficient k for different values of Nus.
Water 09 00417 g007
Figure 8. Effects of the time aggregation on the CV.
Figure 8. Effects of the time aggregation on the CV.
Water 09 00417 g008
Table 1. Principal characteristics of the monitoring systems.
Table 1. Principal characteristics of the monitoring systems.
PSG (I)Fr (NL)CE (I)
Nus239, 777, 981, 12201150596
Type of counterFlow meterFlow meterVolumetric meter
Frequency of sampling1 Hz1/60 Hz1/60 Hz
Period of observation24 November 2008–13 April 200921 January 2006–3 August 20066 January 2009–14 March 2009
Number of used days for the statistics8511346
Table 2. Goodness of fit and statistics for the analyzed samples.
Table 2. Goodness of fit and statistics for the analyzed samples.
Case
study
NusD2%DKSµCpσCp
LNGuLL
PSG12200.1640.11520.06420.11402.500.22
9810.1640.06160.05360.05262.640.27
7770.1640.07410.04980.08932.890.28
2390.1640.13580.12560.12384.040.64
Fr11500.1420.13220.08860.11582.840.27
CE5960.2230.04320.06730.05362.720.27

Share and Cite

MDPI and ACS Style

Gargano, R.; Tricarico, C.; Granata, F.; Santopietro, S.; De Marinis, G. Probabilistic Models for the Peak Residential Water Demand. Water 2017, 9, 417. https://doi.org/10.3390/w9060417

AMA Style

Gargano R, Tricarico C, Granata F, Santopietro S, De Marinis G. Probabilistic Models for the Peak Residential Water Demand. Water. 2017; 9(6):417. https://doi.org/10.3390/w9060417

Chicago/Turabian Style

Gargano, Rudy, Carla Tricarico, Francesco Granata, Simone Santopietro, and Giovanni De Marinis. 2017. "Probabilistic Models for the Peak Residential Water Demand" Water 9, no. 6: 417. https://doi.org/10.3390/w9060417

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