Quantities of interest for surface based resistivity geophysical measurements

The objective of traditional goal-oriented strategies is to construct an optimal mesh that minimizes the problem size needed to achieve a user prescribed tolerance error for a given quantity of interest (QoI). Typical geophysical resistivity measurement acquisition systems can easily record electromagnetic (EM) fields. However, depending upon the application, EM fields are sometimes loosely related to the quantity that is to be inverted (conductivity or resistivity), and therefore they become inadequate for inversion. In the present work, we study the impact of the selection of the QoI in our inverse problem. We focus on two different acquisition systems: marine controlled source electromagnetic (CSEM), and magnetotellurics (MT). For both applications, numerical results illustrate the benefits of employing adequate QoI. Specifically, the use as QoI of the impedance matrix on MT measurements provides significant computational savings, since one can replace the existing absorbing boundary conditions (BCs) by a homogeneous Dirichlet BC to truncate the computational domain, something that is not possible when considering EM fields as QoI

Heat as a proxy to image dynamic processes with 4D electrical resistivity tomography

Since salt cannot always be used as a geophysical tracer (because it may pollute the aquifer with the mass that is necessary to induce a geophysical contrast), and since in many contaminated aquifer salts (e.g., chloride) already constitute the main contaminants, another geophysical tracer is needed to force a contrast in the subsurface that can be detected from surface geophysical measurements. In this context, we used heat as a proxy to image and monitor groundwater flow and solute transport in a shallow alluvial aquifer (< 10 m deep) with the help of electrical resistivity tomography (ERT). The goal of our study is to demonstrate the feasibility of such methodology in the context of the validation of the efficiency of a hydraulic barrier that confines a chloride contamination to its source. To do so, we combined a heat tracer push/pull test with time-lapse 3D ERT and classical hydrogeological measurements in wells and piezometers. Our results show that heat can be an excellent salt substitution tracer for geophysical monitoring studies, both qualitatively and semi-quantitatively. Our methodology, based on 3D surface ERT, allows to visually prove that a hydraulic barrier works efficiently and could be used as an assessment of such installations

Imaging of a fluid injection process using geophysical data - A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts

3-D Magnetotelluric Investigations for geothermal exploration in Martinique (Lesser Antilles). Characteristic Deep Resistivity Structures, and Shallow Resistivity Distribution Matching Heliborne TEM Results

Within the framework of a global French program oriented towards the development of renewable energies, Martinique Island (Lesser Antilles, France) has been extensively investigated (from 2012 to 2013) through an integrated multi-methods approach, with the aim to define precisely the potential geothermal ressources, previously highlighted (Sanjuan et al., 2003). Amongst the common investigation methods deployed, we carried out three magnetotelluric (MT) surveys located above three of the most promising geothermal fields of Martinique, namely the Anses d'Arlet, the Montagne Pel{\'e}e and the Pitons du Carbet prospects. A total of about 100 MT stations were acquired showing single or multi-dimensional behaviors and static shift effects. After processing data with remote reference, 3-D MT inversions of the four complex elements of MT impedance tensor without pre-static-shift correction, have been performed for each sector, providing three 3-D resistivity models down to about 12 to 30 km depth. The sea coast effect has been taken into account in the 3-D inversion through generation of a 3-D resistivity model including the bathymetry around Martinique from the coast up to a distance of 200 km. The forward response of the model is used to calculate coast effect coefficients that are applied to the calculated MT response during the 3-D inversion process for comparison with the observed data. 3-D resistivity models of each sector, which are inherited from different geological history, show 3-D resistivity distribution and specificities related to its volcanological history. In particular, the geothermal field related to the Montagne Pel{\'e}e strato-volcano, is characterized by a quasi ubiquitous conductive layer and quite monotonic typical resistivity distribution making interpretation difficult in terms of geothermal targets. At the opposite, the resistivity distribution of Anse d'Arlet area is radically different and geothermal target is thought to be connected to a not so deep resistive intrusion elongated along a main structural axis. Beside these interesting deep structures, we demonstrate, after analyzing the results of the recent heliborne TEM survey covering the whole Martinique, that surface resistivity distribution obtained from 3-D inversion reproduce faithfully the resistivity distribution observed by TEM. In spite of a very different sampling scale, this comparison illustrates the ability of 3-D MT inversion to take into account and reproduce static shift effects in the sub-surface resistivity distribution.Comment: Wordl Geothermal Congress 2015, Apr 2015, Melbourne, Australi

An evaluation of combined geophysical and geotechnical methods to characterize beach thickness

Beaches provide sediment stores and have an important role in the development of the coastline in response to climate change. Quantification of beach thickness and volume is required to assess coastal sediment transport budgets. Therefore, portable, rapid, non-invasive techniques are required to evaluate thickness where environmental sensitivities exclude invasive methods. Site methods and data are described for a toolbox of electrical, electromagnetic, seismic and mechanical based techniques that were evaluated at a coastal site at Easington, Yorkshire. Geophysical and geotechnical properties are shown to be dependent upon moisture content, porosity and lithology of the beach and the morphology of the beach–platform interface. Thickness interpretation, using an inexpensive geographic information system to integrate data, allowed these controls and relationships to be understood. Guidelines for efficient site practices, based upon this case history including procedures and techniques, are presented using a systematic approach. Field results indicated that a mixed sand and gravel beach is highly variable and cannot be represented in models as a homogeneous layer of variable thickness overlying a bedrock half-space

Estimation of subsurface porosities and thermal conductivities of polygonal tundra by coupled inversion of electrical resistivity, temperature, and moisture content data

Studies indicate greenhouse gas emissions following permafrost thaw will amplify current rates of atmospheric warming, a process referred to as the permafrost carbon feedback. However, large uncertainties exist regarding the timing and magnitude of the permafrost carbon feedback, in part due to uncertainties associated with subsurface permafrost parameterization and structure. Development of robust parameter estimation methods for permafrost-rich soils is becoming urgent under accelerated warming of the Arctic. Improved parameterization of the subsurface properties in land system models would lead to improved predictions and a reduction of modeling uncertainty. In this work we set the groundwork for future parameter estimation (PE) studies by developing and evaluating a joint PE algorithm that estimates soil porosities and thermal conductivities from time series of soil temperature and moisture measurements and discrete in-time electrical resistivity measurements. The algorithm utilizes the Model-Independent Parameter Estimation and Uncertainty Analysis toolbox and coupled hydrological-thermal-geophysical modeling. We test the PE algorithm against synthetic data, providing a proof of concept for the approach. We use specified subsurface porosities and thermal conductivities and coupled models to set up a synthetic state, perturb the parameters, and then verify that our PE method is able to recover the parameters and synthetic state. To evaluate the accuracy and robustness of the approach we perform multiple tests for a perturbed set of initial starting parameter combinations. In addition, we varied types and quantities of data to better understand the optimal dataset needed to improve the PE method. The results of the PE tests suggest that using multiple types of data improve the overall robustness of the method. Our numerical experiments indicate that special care needs to be taken during the field experiment setup so that (1) the vertical distance between adjacent measurement sensors allows the signal variability in space to be resolved and (2) the longer time interval between resistivity snapshots allows signal variability in time to be resolved

Simple methods for assessing groundwater resources in low permeability areas of Africa

In many areas throughout Africa, a staggering proportion of wells and boreholes fail. Failure can occur for a number of reasons – inadequate maintenance and community involvement, poor engineering or a lack of water. Often it can be difficult to work out the exact reason after the event. However, in many geological environments the impacts of poorly sited and designed boreholes and wells are a major concern to funding agencies, implementing institutions and local communities. In such areas, good supplies of groundwater cannot be found everywhere, and boreholes and wells must be sited and designed carefully to make use of the available groundwater. To appropriately site and design water sources, the groundwater resources of an area need first to be investigated to understand how water occurs in the ground. In this manual we present some techniques that allow a quick assessment of groundwater resources without requiring much expertise or expense. Some of the techniques are old and established while others are new. However, all techniques have been tested by BGS (and others) in assessing groundwater resources in Africa. This manual does not claim to be a detailed textbook for hydrogeologists – there are enough already (see reading list at the end of the chapter). Rather it is meant as a practicalaid for those involved in the practice of rural water supply, particularly in Africa. Little training or equipment is required for the tests and they can all be carried out in a short space of time. The manual is divided into six sections. The first gives an overview of the groundwater resources of sub-Saharan Africa (SSA) and discusses the scope and detail of investigations required in different geological environments. The remaining chapters describe simple techniques for assessing groundwater resources, from basic reconnaissance to assessing the yield of a borehole. In the appendix are summary sheets of the most common techniques which can be photocopied and used in the field

The α\alpha Dynamo Effects in Laboratory Plasmas

A concise review of observations of the $\alpha$ dynamo effect in laboratory plasmas is given. Unlike many astrophysical systems, the laboratory pinch plasmas are driven magnetically. When the system is overdriven, the resultant instabilities cause magnetic and flow fields to fluctuate, and their correlation induces electromotive forces along the mean magnetic field. This $\alpha$-effect drives mean parallel electric current, which, in turn, modifies the initial background mean magnetic structure towards the stable regime. This drive-and-relax cycle, or the so-called self-organization process, happens in magnetized plasmas in a time scale much shorter than resistive diffusion time, thus it is a fast and unquenched dynamo process. The observed $\alpha$-effect redistributes magnetic helicity (a measure of twistedness and knottedness of magnetic field lines) but conserves its total value. It can be shown that fast and unquenched dynamos are natural consequences of a driven system where fluctuations are statistically either not stationary in time or not homogeneous in space, or both. Implications to astrophysical phenomena will be discussed.Comment: 21 pages, 15 figures, submitted to Magnetohydrodynamic