Data-domain correlation approach for joint hydrogeologic inversion of time-lapse hydrogeologic and geophysical data

Inverse estimations of hydrogeologic properties often are highly uncertain because of the expense of collecting hydrogeologic data and the subsequent lack of information. Geophysical data potentially can help fill this information gap because geophysical methods can survey large areas remotely and relatively inexpensively. However, geophysical data are difficult to incorporate into hydrogeologic parameter estimations primarily because of a lack of knowledge concerning the petrophysical relationships between hydrogeologic and geophysical parameters. A method can be used that allows time-lapse geophysical data to be incorporated directly into a hydrogeologic parameter estimation when a strong correlation exists between changes in geophysical and hydrogeologic properties. This approach bypasses the need for an explicit petrophysical transform by formulating the geophysical part of the hydrogeologic inversion in terms of a data-domain correlation operator. A synthetic electrical resistivity monitoring application is used to estimate the hydraulic conductivity distribution. Including time-lapse resistivity data to supplement sparse hydrologic data appears to improve greatly the resolution of hydraulic conductivity in this case. More generally, the formulation and results suggest that geophysical monitoring data can be incorporated effectively into a hydrogeologic parameter estimation using a data-domain correlation operator, assuming a strong correlation exists between changes in hydrogeologic and geophysical properties.

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

[2]  T. Günther,et al.  Three‐dimensional modelling and inversion of dc resistivity data incorporating topography – II. Inversion , 2006 .

[3]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[4]  Partha S. Routh,et al.  Time-lapse ERT monitoring of an injection/withdrawal experiment , 2006 .

[5]  Andrew Binley,et al.  Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution‐dependent limitations , 2005 .

[6]  Arlen W. Harbaugh,et al.  User's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model , 1996 .

[7]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[8]  Kamini Singha,et al.  Effects of spatially variable resolution on field-scale estimates of tracer concentration from electrical inversions using Archie’s law , 2006 .

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

[10]  James P. McNamara,et al.  Application of time-lapse ERT imaging to watershed characterization , 2008 .

[11]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[12]  Per Christian Hansen,et al.  Regularization methods for large-scale problems , 1993 .

[13]  S. Gorelick,et al.  Hydrogeophysical tracking of three‐dimensional tracer migration: The concept and application of apparent petrophysical relations , 2006 .

[14]  Douglas LaBrecque,et al.  Difference Inversion of ERT Data: a Fast Inversion Method for 3-D In Situ Monitoring , 2001 .

[15]  J. Harris,et al.  Traveltime inversion for the geometry of aquifer lithologies , 1996 .

[16]  H. Vereecken,et al.  Imaging and characterisation of subsurface solute transport using electrical resistivity tomography (ERT) and equivalent transport models , 2002 .

[17]  Andrew Binley,et al.  Electrical imaging of saline tracer migration for the investigation of unsaturated zone transport mechanisms. , 1997 .

[18]  Catherine Certes,et al.  Application of the pilot point method to the identification of aquifer transmissivities , 1991 .

[19]  Kamini Singha,et al.  A framework for inferring field‐scale rock physics relationships through numerical simulation , 2005 .

[20]  Alan C. Tripp,et al.  Two-dimensional resistivity inversion , 1984 .

[21]  Max A. Meju,et al.  Joint two‐dimensional cross‐gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification , 2007 .