HT2DINV: A 2D forward and inverse code for steady-state and transient hydraulic tomography problems

Abstract Hydraulic tomography is a technique used to characterize the spatial heterogeneities of storativity and transmissivity fields. The responses of an aquifer to a source of hydraulic stimulations are used to recover the features of the estimated fields using inverse techniques. We developed a 2D free source Matlab package for performing hydraulic tomography analysis in steady state and transient regimes. The package uses the finite elements method to solve the ground water flow equation for simple or complex geometries accounting for the anisotropy of the material properties. The inverse problem is based on implementing the geostatistical quasi-linear approach of Kitanidis combined with the adjoint-state method to compute the required sensitivity matrices. For undetermined inverse problems, the adjoint-state method provides a faster and more accurate approach for the evaluation of sensitivity matrices compared with the finite differences method. Our methodology is organized in a way that permits the end-user to activate parallel computing in order to reduce the computational burden. Three case studies are investigated demonstrating the robustness and efficiency of our approach for inverting hydraulic parameters.

[1]  Walter A. Illman,et al.  Capturing aquifer heterogeneity: Comparison of approaches through controlled sandbox experiments , 2011 .

[2]  Yonghong Hao,et al.  Hydraulic Tomography for Detecting Fracture Zone Connectivity , 2008, Ground water.

[3]  John Doherty,et al.  Ground Water Model Calibration Using Pilot Points and Regularization , 2003, Ground water.

[4]  W. Li,et al.  Two‐dimensional characterization of hydraulic heterogeneity by multiple pumping tests , 2007 .

[5]  William W.-G. Yeh,et al.  Coupled inverse problems in groundwater modeling: 2. Identifiability and experimental design , 1990 .

[6]  Wei Li,et al.  Geostatistical inverse modeling of transient pumping tests using temporal moments of drawdown , 2005 .

[7]  C. V. Theis.  The relation between the lowering of the Piezometric surface and the rate and duration of discharge of a well using ground‐water storage , 1935 .

[8]  Peter Dietrich,et al.  Identification of the permeability distribution in soil by hydraulic tomography , 1995 .

[9]  Tian-Chyi J. Yeh,et al.  Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones , 2007 .

[10]  Hiromitsu Saegusa,et al.  Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan , 2009 .

[11]  Jet-Chau Wen,et al.  Necessary conditions for inverse modeling of flow through variably saturated porous media , 2013 .

[12]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[13]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part II: implementation and tests in DUNE , 2008, Computing.

[14]  Tian-Chyi J. Yeh,et al.  Characterization of aquifer heterogeneity using transient hydraulic tomography , 2004 .

[15]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[16]  T. Yeh,et al.  Analysis of hydraulic tomography using temporal moments of drawdown recovery data , 2006 .

[17]  Peter K. Kitanidis,et al.  A field proof‐of‐concept of aquifer imaging using 3‐D transient hydraulic tomography with modular, temporarily‐emplaced equipment , 2012 .

[18]  M. Cardiff,et al.  3‐D transient hydraulic tomography in unconfined aquifers with fast drainage response , 2011 .

[19]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[20]  F. B. Hildebrand,et al.  Methods of applied mathematics , 1953 .

[21]  Andreas Dedner,et al.  A generic grid interface for parallel and adaptive scientific computing. Part I: abstract framework , 2008, Computing.

[22]  Peter K. Kitanidis,et al.  Efficient solution of nonlinear, underdetermined inverse problems with a generalized PDE model , 2008, Comput. Geosci..

[23]  André Revil,et al.  Hydraulic conductivity field characterization from the joint inversion of hydraulic heads and self‐potential data , 2014 .

[24]  Shuyun Liu,et al.  Effectiveness of hydraulic tomography: Sandbox experiments , 2001 .

[25]  Walter A. Illman,et al.  Three‐dimensional transient hydraulic tomography in a highly heterogeneous glaciofluvial aquifer‐aquitard system , 2011 .

[26]  Junfeng Zhu,et al.  Sequential aquifer tests at a well field, Montalto Uffugo Scalo, Italy , 2007 .

[27]  Wei Li,et al.  Efficient parallelization of geostatistical inversion using the quasi-linear approach , 2012, Comput. Geosci..

[28]  Michael A. Saunders,et al.  USER’S GUIDE FOR SNOPT 5.3: A FORTRAN PACKAGE FOR LARGE-SCALE NONLINEAR PROGRAMMING , 2002 .

[29]  P. Kitanidis.  Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[30]  Tian-Chyi J. Yeh,et al.  Robustness of joint interpretation of sequential pumping tests: Numerical and field experiments , 2011 .

[31]  T. Yeh,et al.  Hydraulic tomography: Development of a new aquifer test method , 2000 .