Hydraulic properties at the North Sea island of Borkum derived from joint inversion of magnetic resonance and electrical resistivity soundings

For reliably predicting the impact of climate changes on salt/freshwater systems below barrier islands, a long-term hydraulic modelling is inevitable. As input we need the parameters porosity, salinity and hydraulic conductivity at the catchment scale, preferably non-invasively acquired with geophysical methods. We present a methodology to retrieve the searched parameters and a lithological interpretation by the joint analysis of magnetic resonance soundings (MRS) and vertical electric soundings (VES). Both data sets are jointly inverted for resistivity, water content and decay time using a joint inversion scheme. Coupling is accomplished by common layer thicknesses. We show the results of three soundings measured on the eastern part of the North Sea island of Borkum. Pumping test data is used to calibrate the petrophysical relationship for the local conditions in order to estimate permeability from nuclear magnetic resonance (NMR) data. Salinity is retrieved from water content and resistivity using a modified Archie equation calibrated by local samples. As a result we are able to predict porosity, salinity and hydraulic conductivities of the aquifers, including their uncertainties. The joint inversion significantly improves the reliability of the results. Verification is given by comparison with a borehole. A sounding in the flooding area demonstrates that only the combined inversion provides a correct subsurface model. Thanks to the joint application, we are able to distinguish fluid conductivity from lithology and provide reliable hydraulic parameters as shown by uncertainty analysis. These findings can finally be used to build groundwater flow models for simulating climate changes. This includes the improved geometry and lithological attribution, and also the parameters and their uncertainties.

[1]  Ugur Yaramanci,et al.  Joint inversion of Surface Nuclear Magnetic Resonance and Vertical Electrical Sounding , 2002 .

[2]  Alan G. Green,et al.  Off-resonance effects in surface nuclear magnetic resonance , 2011 .

[3]  Weichman,et al.  Theory of surface nuclear magnetic resonance with applications to geophysical imaging problems , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  D. O. Seevers,et al.  A Nuclear Magnetic Method For Determining The Permeability Of Sandstones , 1966 .

[5]  A. Christiansen,et al.  Improvement in MRS parameter estimation by joint and laterally constrained inversion of MRS and TEM data , 2012 .

[6]  F. Rubio,et al.  The use of MRS in the determination of hydraulic transmissivity: The case of alluvial aquifers , 2008 .

[7]  Anatoly Legchenko,et al.  A review of the basic principles for proton magnetic resonance sounding measurements , 2002 .

[8]  Viezzoli Andrea,et al.  Surface Water-Groundwater Exchange in Transitional Coastal Environments by Airborne Electromagnetics: The Venice Lagoon Example , 2010 .

[9]  Viezzoli Andrea,et al.  Surface Water-Groundwater Exchange in Transitional Coastal Environments by Airborne Electromagnetics: The Venice Lagoon Example , 2010 .

[10]  A. Legchenko,et al.  A new direct non-invasive groundwater detection technology for Australia , 1991 .

[11]  H. Kim,et al.  Inequality constraint in least-squares inversion of geophysical data , 1999 .

[12]  A. Christiansen,et al.  A review of helicopter‐borne electromagnetic methods for groundwater exploration , 2009 .

[13]  Jochen Kamm,et al.  Simultaneous inversion of magnetic resonance sounding in terms of water content, resistivity and decay times , 2009 .

[15]  Ugur Yaramanci,et al.  Pore size distributions and hydraulic conductivities of rocks derived from Magnetic Resonance Sounding relaxation data using multi-exponential decay time inversion , 2008 .

[16]  Walter L. Anderson,et al.  A hybrid fast Hankel transform algorithm for electromagnetic modeling , 1989 .

[17]  A. Green,et al.  Estimating the longitudinal relaxation time T1 in surface NMR , 2011 .

[18]  Ugur Yaramanci,et al.  QT inversion — Comprehensive use of the complete surface NMR data set , 2010 .

[19]  J. Roy,et al.  Hydrogeological interpretation and potential of the new magnetic resonance sounding, MRS, method , 2003 .

[20]  Rosemary Knight,et al.  The effect of pore size and magnetic susceptibility on the surface NMR relaxation parameter T 2 , 2011 .

[21]  J. Roy,et al.  MRS contribution to hydrogeological system parametrization , 2005 .

[22]  Hee Joon Kim,et al.  A unified transformation function for lower and upper bounding constraints on model parameters in electrical and electromagnetic inversion , 2011 .

[23]  A. Christiansen,et al.  A global measure for depth of investigation , 2012 .

[24]  T. Günther,et al.  Numerical modelling of climate change impacts on freshwater lenses on the North Sea Island of Borkum using hydrological and geophysical methods , 2012 .

[25]  Esben Auken,et al.  Efficient full decay inversion of MRS data with a stretched-exponential approximation of the T2* distribution , 2012 .

[26]  M. Grinat,et al.  An Automated Electrical Resistivity Tomography System to Monitor the Freshwater/saltwater Zone on a North Sea Island , 2010 .

[27]  Ugur Yaramanci,et al.  Inversion of resistivity in Magnetic Resonance Sounding , 2008 .

[28]  Thomas Günther,et al.  Assessment of the potential of a new generation of surface nuclear magnetic resonance instruments , 2011 .

[29]  B. Chatenoux,et al.  Efficiency of joint use of MRS and VES to characterize coastal aquifer in Myanmar , 2007 .

[30]  M. H. Waxman,et al.  Electrical Conductivities in Oil-Bearing Shaly Sands , 1968 .

[31]  W. E. Kenyon,et al.  Petrophysical Principles of Applications of NMR Logging , 1997 .

[32]  Marian Hertrich,et al.  Imaging of groundwater with nuclear magnetic resonance , 2008 .

[33]  Jean-Michel Vouillamoz,et al.  Magnetic Resonance Sounding Applied to Aquifer Characterization , 2004, Ground water.

[34]  Ugur Yaramanci,et al.  Evaluation of surface nuclear magnetic resonance-estimated subsurface water content , 2011 .

[35]  J. Roy,et al.  Magnetic Resonance Sounding: New Method for Ground Water Assessment , 2004, Ground water.

[36]  Joseph Robert Inman,et al.  RESISTIVITY INVERSION WITH RIDGE REGRESSION , 1975 .

[37]  T. Günther Inversion methods and resolution analysis for the 2D/3D reconstruction of resistivity structures from DC measurements , 2004 .

[38]  T. Günther,et al.  Numerical modelling of climate change impacts on freshwater lenses on the North Sea Island of Borkum , 2012 .

[39]  P. Lachassagne,et al.  The links between MRS parameters and the hydrogeological parameters , 2005 .