![]() |
Computational & Technology Resources
an online resource for computational,
engineering & technology publications |
Civil-Comp Proceedings
ISSN 1759-3433 CCP: 84
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 155
Calibration of the Model of an Operational Water Distribution System A. Vassiljev and T. Koppel
Department of Mechanics, Tallinn University of Technology, Estonia Full Bibliographic Reference for this paper
A. Vassiljev, T. Koppel, "Calibration of the Model of an Operational Water Distribution System", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Fifth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 155, 2006. doi:10.4203/ccp.84.155
Keywords: calibration, optimization, response surface, Levenberg-Marquardt, Epanet, water distribution system.
Summary
Any water distribution system (WDS) model needs to be calibrated. The accuracy of
the hydraulic model depends on how well it has been calibrated. Calibration
parameters usually include pipe roughness, pipe diameters and demands, where the
two first belong to flow conditions and the last one to boundary conditions.
Economic development in Estonia and the accompanying increase of water price
have led to the situation in which practically each consumer has a water flow meter.
Therefore, only roughness of pipes is considered in this paper.
Calibration of pipe roughness is a difficult task if the WDS contains pipes of
different ages. In an ideal case roughness coefficients for all pipes must be calibrated.
However, measurements are quite expensive and therefore we always have much
fewer measurements than we require for calibration. Even with an extensive data
collection effort, pipe roughnesses for all links cannot be determined exactly [1].
Thus, a major problem associated with model calibration is the need to calibrate a
large number of pipes using only a few measurements. Therefore, grouping of pipes
is widely used at calibration [1]. Clustering of pipes is a difficult task for old WDS
because great age differences of pipes lead to a large number of clusters. Therefore,
instead of grouping of pipes with different age we propose to describe dependence
of roughness on age by some flexible function with a small number of parameters. It
may be as follows for the pipe roughness
where ![]() ![]() ![]() ![]() ![]() Approximation (87) is quite flexible and describes both the linear and nonlinear dependencies with different degree of nonlinearity.
Parameters
Our calculations for a real WDS, containing thousands of pipes, showed that the Levenberg-Marquardt algorithm failed because of oscillations of the objective function. Calculations of the objective function showed that its surface looked like emery paper with a large number of small local maximums and minimums. At the same time, the surface would be very good for the Levenberg-Marquardt algorithm (one minimum with good gradients in any direction) if we could eliminate these small variations. The simplest way to do it is to select a long enough step size for the parameters, when calculating partial derivatives for the Levenberg-Marquardt method [2]. Two methods have been used to estimate the step. One of them is based on the number of reliable digits in the objective function values [2]. The other method, proposed in the paper, is based on trying different step sizes to find the step size, the Levenberg-Marquardt method can work with. A WDS representing a part of the whole WDS of the City of Tallinn was used as an example of the real WDS. This WDS includes thousands of pipes. The age of pipes varies from 0 to 41 years. Calculations showed that the Levenberg-Marquardt method worked very well with a long enough step size for calculations of the gradient of the objective function. The surface of the objective function obtained by the trial-and-error method facilitates calibration significantly. Visualization of the response surface discovers lack of data and helps to plan measurements. It must also be mentioned that we used the TOOLKIT developed by Rossman [3] for the EPANET2 to complete the calculations automatically. The TOOLKIT permits the creation of external programs, which govern the calculations and therefore is much better for calibration than commercial programs. References
purchase the full-text of this paper (price £20)
go to the previous paper |
|