Search Results

Now showing 1 - 2 of 2
  • Item
    pyGIMLi: An open-source library for modelling and inversion in geophysics
    (Amsterdam [u.a.] : Elsevier Science, 2017) Rücker, Carsten; Günther, Thomas; Wagner, Florian M.
    Many tasks in applied geosciences cannot be solved by single measurements, but require the integration of geophysical, geotechnical and hydrological methods. Numerical simulation techniques are essential both for planning and interpretation, as well as for the process understanding of modern geophysical methods. These trends encourage open, simple, and modern software architectures aiming at a uniform interface for interdisciplinary and flexible modelling and inversion approaches. We present pyGIMLi (Python Library for Inversion and Modelling in Geophysics), an open-source framework that provides tools for modelling and inversion of various geophysical but also hydrological methods. The modelling component supplies discretization management and the numerical basis for finite-element and finite-volume solvers in 1D, 2D and 3D on arbitrarily structured meshes. The generalized inversion framework solves the minimization problem with a Gauss-Newton algorithm for any physical forward operator and provides opportunities for uncertainty and resolution analyses. More general requirements, such as flexible regularization strategies, time-lapse processing and different sorts of coupling individual methods are provided independently of the actual methods used. The usage of pyGIMLi is first demonstrated by solving the steady-state heat equation, followed by a demonstration of more complex capabilities for the combination of different geophysical data sets. A fully coupled hydrogeophysical inversion of electrical resistivity tomography (ERT) data of a simulated tracer experiment is presented that allows to directly reconstruct the underlying hydraulic conductivity distribution of the aquifer. Another example demonstrates the improvement of jointly inverting ERT and ultrasonic data with respect to saturation by a new approach that incorporates petrophysical relations in the inversion. Potential applications of the presented framework are manifold and include time-lapse, constrained, joint, and coupled inversions of various geophysical and hydrological data sets.
  • Item
    Noninvasive Estimation of Water Retention Parameters by Observing the Capillary Fringe with Magnetic Resonance Sounding
    (Hoboken, NJ : Wiley, 2014) Costabel, Stephan; Günther, Thomas
    The magnetic resonance sounding (MRS) method is usually applied for delineation and characterization of aquifer system stratification. Its unique property, distinct from other hydrogeophysical methods, is the direct sensitivity to water content in the subsurface. The inversion of MRS data yields the subsurface water content distribution without need of a petrophysical model. Recent developments in instrumentation, i.e., decreased instrumental dead times and advanced noise cancellation strategies, enable the use of this method for investigating the vadose zone. A possible way to interpret MRS measurements with focus on water retention (WR) parameters is an inversion approach that directly provides WR parameters by modeling the capillary fringe (CF inversion). We have developed this kind of inversion further to account for different WR models and present a sensitivity study based on both synthetic and real field data. To assess the general applicability of the CF inversion, we analyzed the resolution properties for different measurement layouts and the parameter uncertainties for different realistic scenarios. Under moderate noise conditions and if the water table position is known, all WR parameters except the residual water content can be reliably estimated. The relative accuracy of the estimated pore distribution index estimation is better for larger CF. Small measurement loops of 5-m diameter achieve the best resolution for shallow investigation depths of <10 m.