Joint hydrogeophysical inversion: state estimation for seawater intrusion models in 3D

Seawater intrusion (SWI) is a complex process, where 3D modeling is often necessary in order to monitor and manage the affected aquifers. Here, we present a synthetic study to test a joint hydrogeophysical inversion approach aimed at solving the inverse problem of estimating initial and current saltwater distribution. First, we use a 3D groundwater model for variable density flow based on discretized flow and solute mass balance equations. In addition to the groundwater model, a 3D geophysical model was developed for direct current resistivity imaging and inversion. The objective function of the coupled problem consists of data misfit and regularization terms as well as a coupling term that relates groundwater and geophysical states. We present a novel approach to solve the inverse problem using an alternating direction method of multipliers (ADMM) to minimize this coupled objective function. ADMM enables to treat the groundwater and geophysical part separately and thus use the existing software with minor changes. To further reduce the computational cost, the sensitivities are derived analytically for the discretized system of equations, which allows us to efficiently compute the gradients in the minimization procedure. The method was tested on different synthetic scenarios with groundwater and geophysical data represented by solute mass fraction data and direct current resistivity data. With the ADMM approach, we were able to obtain better estimates for the solute distribution compared to just considering each data separately, solving the problem with a simple coupled approach or by a direct substitution of the coupling constraint.

[1]  Weixing Guo,et al.  MODFLOW/MT3DMS–Based Simulation of Variable‐Density Ground Water Flow and Transport , 2006, Ground water.

[2]  J. Hem.  Study and Interpretation of the Chemical Characteristics of Natural Water , 1989 .

[3]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[4]  J. Harris,et al.  Coupled seismic and tracer test inversion for aquifer property characterization , 1993 .

[5]  Eldad Haber,et al.  Model Fusion and Joint Inversion , 2013, Surveys in Geophysics.

[6]  Frédéric Nguyen,et al.  Calibration of seawater intrusion models: Inverse parameter estimation using surface electrical resistivity tomography and borehole data , 2014 .

[7]  Stanley H. Ward,et al.  The Resistivity And Induced Polarization Methods , 1988 .

[8]  E. Haber,et al.  Joint inversion: a structural approach , 1997 .

[9]  Alyssa M. Dausman,et al.  SEAWAT Version 4: A Computer Program for Simulation of Multi-Species Solute and Heat Transport , 2008 .

[10]  Andrew Binley,et al.  Monitoring Unsaturated Flow and Transport Using Cross‐Borehole Geophysical Methods , 2008 .

[11]  J. S. Perret,et al.  Control of sea-water intrusion by salt-water pumping: Coast of Oman , 2009 .

[12]  David V. Fitterman,et al.  Mapping Saltwater Intrusion in the Biscayne Aquifer, Miami-Dade County, Florida using Transient Electromagnetic Sounding , 2014 .

[13]  P. Barlow,et al.  Saltwater intrusion in coastal regions of North America , 2010 .

[14]  A. Binley,et al.  DC Resistivity and Induced Polarization Methods , 2005 .

[15]  A. C. Hinnell,et al.  Improved extraction of hydrologic information from geophysical data through coupled hydrogeophysical inversion , 2010 .

[16]  C. Langevin,et al.  Effect of Sea‐Level Rise on Salt Water Intrusion near a Coastal Well Field in Southeastern Florida , 2013, Ground water.

[17]  Michael I. Jordan,et al.  A General Analysis of the Convergence of ADMM , 2015, ICML.

[18]  D. Oldenburg,et al.  A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems , 2004 .

[19]  Martin Rumpf,et al.  A Variational Approach to Nonrigid Morphological Image Registration , 2004, SIAM J. Appl. Math..

[20]  Y. Rubin,et al.  Hydrogeological parameter estimation using geophysical data: a review of selected techniques , 2000 .

[21]  Application and evaluation of electromagnetic methods for imaging saltwater intrusion in coastal aquifers: Seaside Groundwater Basin, California , 2013 .

[22]  Peter Bauer-Gottwein,et al.  Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data , 2012, Ground water.

[23]  A. Cheng,et al.  Seawater intrusion in coastal aquifers : concepts, methods, and practices , 1999 .

[24]  Olaf A. Cirpka,et al.  Fully coupled hydrogeophysical inversion of a laboratory salt tracer experiment monitored by electrical resistivity tomography , 2012 .

[25]  G. O. Essink,et al.  Shallow rainwater lenses in deltaic areas with saline seepage , 2011 .

[26]  O. Kolditz,et al.  Coupled groundwater flow and transport : 1. Verification of variable density flow and transport models , 1998 .

[27]  P. Hoekstra,et al.  Time Domain Electromagnetic Soundings for Mapping Sea‐Water Intrusion in Monterey County, California , 1988 .

[28]  Michael Commer,et al.  Three-dimensional controlled-source electromagnetic and magnetotelluric joint inversion , 2009 .

[29]  J. Tarhouni,et al.  Delineation of saltwater intrusion zones using the time domain electromagnetic method: the Nabeul–Hammamet coastal aquifer case study (NE Tunisia) , 2013 .

[30]  E. Auken,et al.  Sequential and joint hydrogeophysical inversion using a field-scale groundwater model with ERT and TDEM data , 2013 .

[31]  A new coupling algorithm for density‐driven flow in porous media , 2004 .

[32]  M. Antonellini,et al.  Salt water intrusion in the coastal aquifer of the southern Po Plain, Italy , 2008 .

[33]  Frédéric Nguyen,et al.  Characterization of seawater intrusion using 2D electrical imaging , 2009 .

[34]  Wotao Yin,et al.  On the Global and Linear Convergence of the Generalized Alternating Direction Method of Multipliers , 2016, J. Sci. Comput..

[35]  Mary C. Hill,et al.  UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API , 2006 .

[36]  Roland Martin,et al.  Imaging artificial salt water infiltration using electrical resistivity tomography constrained by geostatistical data , 2012 .

[37]  E. Custodio.  Coastal aquifers of Europe: an overview , 2010 .

[38]  Stefan Finsterle,et al.  Coupled modeling of hydrogeochemical and electrical resistivity data for exploring the impact of recharge on subsurface contamination , 2011 .

[39]  Max A. Meju,et al.  Structure‐coupled multiphysics imaging in geophysical sciences , 2011 .

[40]  C. Vogel.  Computational Methods for Inverse Problems , 1987 .

[41]  E. Haber,et al.  RESINVM3D: A 3D resistivity inversion package , 2007 .

[42]  Rita Deiana,et al.  A tracer test in a shallow heterogeneous aquifer monitored via time-lapse surface electrical resistivity tomography , 2010 .

[43]  G. E. Archie.  The electrical resistivity log as an aid in determining some reservoir characteristics , 1942 .

[44]  A. Revil,et al.  Stochastic joint inversion of hydrogeophysical data for salt tracer test monitoring and hydraulic conductivity imaging , 2013 .

[45]  Luis Tenorio,et al.  Prior information and uncertainty in inverse problems , 2001 .

[46]  T. F. Russell,et al.  A Finite‐Volume ELLAM for Three‐Dimensional Solute‐Transport Modeling , 2003, Ground water.

[47]  O. Banton,et al.  Cross‐validation of geo‐electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers , 2007 .

[48]  P. Ackerer,et al.  A new efficient Eulerian–Lagrangian localized adjoint method for solving the advection–dispersion equation on unstructured meshes , 2006 .

[49]  A. Englert,et al.  SOLUTE TRANSPORT PROCESSES , 2006 .

[50]  T. F. Russell,et al.  An overview of research on Eulerian-Lagrangian localized adjoint methods (ELLAM). , 2002 .

[51]  L. Slooten,et al.  Computational and conceptual issues in the calibration of seawater intrusion models , 2010 .

[52]  Stefan Finsterle,et al.  MPiTOUGH2: A parallel parameter estimation framework for hydrological and hydrogeophysical applications , 2014, Comput. Geosci..

[53]  A. Cardona,et al.  Salinization in coastal aquifers of arid zones: an example from Santo Domingo, Baja California Sur, Mexico , 2004 .

[54]  Shaul Sorek,et al.  Modified Eulerian-Lagrangian formulation for hydrodynamic modeling , 2012, J. Comput. Phys..

[55]  Eldad Haber,et al.  Computational Methods in Geophysical Electromagnetics , 2014, Mathematics in Industry.

[56]  John Doherty,et al.  Approaches to Highly Parameterized Inversion: A Guide to Using PEST for Model-Parameter and Predictive-Uncertainty Analysis , 2014 .

[57]  Jesús Carrera,et al.  Coupled estimation of flow and solute transport parameters , 1996 .

[58]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[59]  Daniel M. Tartakovsky,et al.  Linear functional minimization for inverse modeling , 2015 .

[60]  P. Hansen.  Discrete Inverse Problems: Insight and Algorithms , 2010 .

[61]  G. Z. Hornberger,et al.  A three-dimensional method-of-characteristics solute-transport model (MOC3D) , 1996 .

[62]  P. Bauer‐Gottwein,et al.  Hydrogeophysical exploration of three-dimensional salinity anomalies with the time-domain electromagnetic method (TDEM) , 2010 .

[63]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[64]  Euhanna Ghadimi,et al.  Optimal Parameter Selection for the Alternating Direction Method of Multipliers (ADMM): Quadratic Problems , 2013, IEEE Transactions on Automatic Control.

[65]  G. O. Essink.  Modeling Three-Dimensional Density Dependent Groundwater Flow at the Island of Texel, The Netherlands , 2003 .

[66]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .

[67]  C. Paniconi,et al.  Picard and Newton linearization for the coupled model for saltwater intrusion in aquifers , 1995 .

[68]  L. Smith,et al.  Saline water intrusion adjacent to the Fraser River, Richmond, British Columbia , 2001 .

[69]  E. Custodio,et al.  State of knowledge of coastal aquifer management in South America , 2010 .

[70]  W. Sanford,et al.  Current challenges using models to forecast seawater intrusion: lessons from the Eastern Shore of Virginia, USA , 2010 .

[71]  Jesús Carrera,et al.  Anisotropic dispersive Henry problem , 2007 .

[72]  K. Singha,et al.  Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities , 2010 .

[73]  Warren Barrash,et al.  Hydrological parameter estimations from a conservative tracer test with variable‐density effects at the Boise Hydrogeophysical Research Site , 2011 .

[74]  LINEARIZED FUNCTIONAL MINIMIZATION FOR INVERSE MODELING , 2012 .