skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Bayesian Update and Method of Distributions: Application to Leak Detection in Transmission Mains

Water‐hammer equations are used to describe transient flow in pipe networks. Uncertainty in model parameters, initial and boundary conditions, and location and strength of a possible leak renders deterministic predictions of this system untenable. When deployed in conjunction with pressure measurements, probabilistic solutions of the water‐hammer equations serve as a tool for detecting leaks in pipes. We use the method of distributions to obtain a probability density function (PDF) for pressure head, whose dynamics are described by the stochastic water‐hammer equations. This PDF provides a prior distribution for subsequent Bayesian data assimilation, in which data collected by pressure sensors are combined with this prior to obtain a posterior PDF of the leak location and size. We conduct a series of numerical experiments with uncertain initial velocity and measurement noise to ascertain the robustness and accuracy of the proposed approach. The results show the method's ability to identify the leak location and strength in a water transmission main.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Leakages in water distribution networks (WDNs) are estimated to globally cost 39 billion USD per year and cause water and revenue losses, infrastructure degradation, and other cascading effects. Their impacts can be prevented and mitigated with prompt identification and accurate leak localization. In this work, we propose the leakage identification and localization algorithm (LILA), a pressure-based algorithm for data-driven leakage identification and model-based localization in WDNs. First, LILA identifies potential leakages via semi-supervised linear regression of pairwise sensor pressure data and provides the location of their nearest sensors. Second, LILA locates leaky pipes relying on an initial set of candidate pipes and a simulation-based optimization framework with iterative linear and mixed-integer linear programming. LILA is tested on data from the L-Town network devised for the Battle of Leakage Detection and Isolation Methods. Results show that LILA can identify all leakages included in the data set and locate them within a maximum distance of 374 m from their real location. Abrupt leakages are identified immediately or within 2 h, while more time is required to raise alarms on incipient leakages. 
    more » « less
  2. Abstract

    Fluid injection stimulates seismicity far from active tectonic regions. However, the details of how fluids modify on‐fault stresses and initiate seismic events remain poorly understood. We conducted laboratory experiments using a biaxial loading apparatus with a 3 m saw‐cut granite fault and compared events induced at different levels of background shear stress. Water was injected at 10 mL/min and normal stress was constant at 4 MPa. In all experiments, aseismic slip initiated on the fault near the location of fluid injection and dynamic rupture eventually initiated from within the aseismic slipping patch. When the fault was near critically stressed, seismic slip initiated only seconds after MPa‐level injection pressures were reached and the dynamic rupture propagated beyond the fluid pressure perturbed region. At lower stress levels, dynamic rupture initiated hundreds of seconds later and was limited to regions where aseismic slip had significantly redistributed stress from within the pressurized region to neighboring locked patches. We found that the initiation of slow slip was broadly consistent with a Coulomb failure stress, but that initiation of dynamic rupture required additional criteria to be met. Even high background stress levels required aseismic slip to modify on‐fault stress to meet initiation criteria. We also observed slow slip events prior to dynamic rupture. Overall, our experiments suggest that initial fault stress, relative to fault strength, is a critical factor in determining whether a fluid‐induced rupture will “runaway” or whether a fluid‐induced rupture will remain localized to the fluid pressurized region.

    more » « less
  3. Abstract

    Interaction of magma with groundwater or surface water can lead to explosive phreatomagmatic eruptions. Questions of this process center on effects of system geometry and length and time scales, and these necessitate experiments at larger scale than previously conducted in order to investigate the thermohydraulic escalation behavior of rapid heat transfer. Previous experimental work either realized melt‐water interaction at similar meter scales, using a thermite‐based magma analog in a confining vessel, or on smaller scale using ≃0.4 kg remelted volcanic rock in an open crucible, with controlled premix and a ≃5 J kinetic energy trigger event. The new setup uses 55 kg melt for interaction, and the timing and location of the magma‐water premix can be controlled on a scale up to 1 m. A trigger mechanism is a falling hammer that drives a plunger into the melt (≃28 J kinetic energy). Results show intense interaction () at relatively low magma/water mass ratio. A video analysis quantifies rate and amount of melt ejection and compares results with those using the same melt in the smaller scale setup. Experiments show that on the meter scale intense interaction can start spontaneously without an external trigger if the melt column above the initial mixing location is larger than 0.3 m. Experiments suggest that buoyant rise of water domains in a melt column could promote explosive interaction. We assess interaction scenarios between introduced water domains and a magma column, some of which could result in eruptions of Strombolian style, promoting brecciation and incorporation of wall rock debris into a magma column.

    more » « less
  4. Water distribution systems (WDSs) face a significant challenge in the form of pipe leaks. Pipe leaks can cause loss of a large amount of treated water, leading to pressure loss, increased energy costs, and contamination risks. Locating pipe leaks has been a constant challenge for water utilities and stakeholders due to the underground location of the pipes. Physical methods to detect leaks are expensive, intrusive, and heavily localized. Computational approaches provide an economical alternative to physical methods. Data-driven machine learning-based computational approaches have garnered growing interest in recent years to address the challenge of detecting pipe leaks in WDSs. While several studies have applied machine learning models for leak detection on single pipes and small test networks, their applicability to the real-world WDSs is unclear. Most of these studies simplify the leak characteristics and ignore modeling and measuring device uncertainties, which makes the scalability of their approaches questionable to real-world WDSs. Our study addresses this issue by devising four study cases that account for the realistic leak characteristics (multiple, multi-size, and randomly located leaks) and incorporating noise in the input data to account for the model- and measuring device- related uncertainties. A machine learning-based approach that uses simulated pressure as input to predict both location and size of leaks is proposed. Two different machine learning models: Multilayer Perceptron (MLP) and Convolutional Neural Network (CNN), are trained and tested for the four study cases, and their performances are compared. The precision and recall results for the L-Town network indicate good accuracies for both the models for all study cases, with CNN generally outperforming MLP. 
    more » « less
  5. Abstract

    To better quantify how injection, prior seismicity, and fault properties control rupture growth and propagation of induced earthquakes, we perform a finite‐fault slip inversion on aMw4.0 earthquake that occurred in April 2015, the largest earthquake in an induced sequence near Guthrie, Oklahoma. The slip inversion reveals a rupture with slip patches that are anti‐correlated to the locations of prior seismicity. The prior seismicity driven by low pore pressure changes and static stress changes occurred on weaker portions of the fault, while theMw4.0 earthquake likely ruptured relatively stronger portions of the fault. To resolve if pore pressure changes or the initial underlying stress distribution and fault strength controlled the final slip distribution of the GuthrieMw4.0 earthquake, we compare strike‐slip events of similar magnitude from tectonically active regions and previously inactive regions. Earthquakes on reactivated faults exhibit different slip distributions than active regions, they have more prominent and well separated slip patches, a behavior often associated with faults of lower fault maturity. Pore pressure shows little effect on the distributions. These observations suggest that the initial underlying stress distribution and fault strength of reactivated faults in low deformation regions is the primary controlling factor of the slip distribution with pore pressure perturbations and earthquake interactions being secondary. Therefore, GuthrieMw4.0 earthquakes slip distribution was enhanced by pore‐pressure perturbations and earthquake interactions by creating an optimal stress state for its failure, but the slip distribution itself is controlled by its fault's initial stress and strength state.

    more » « less