Probabilistic electrical resistivity tomography of a CO2 sequestration analog

Electrical resistivity tomography (ERT) is a well-established method for geophysical characterization and has shown potential for monitoring geologic CO2 sequestration, due to its sensitivity to electrical resistivity contrasts generated by liquid/gas saturation variability. In contrast to deterministic inversion approaches, probabilistic inversion provides the full posterior probability density function of the saturation field and accounts for the uncertainties inherent in the petrophysical parameters relating the resistivity to saturation. In this study, the data are from benchtop ERT experiments conducted during gas injection into a quasi-2D brine-saturated sand chamber with a packing that mimics a simple anticlinal geological reservoir. The saturation fields are estimated by Markov chain Monte Carlo inversion of the measured data and compared to independent saturation measurements from light transmission through the chamber. Different model parameterizations are evaluated in terms of the recovered saturation and petrophysical parameter values. The saturation field is parameterized (1) in Cartesian coordinates, (2) by means of its discrete cosine transform coefficients, and (3) by fixed saturation values in structural elements whose shape and location is assumed known or represented by an arbitrary Gaussian Bell structure. Results show that the estimated saturation fields are in overall agreement with saturations measured by light transmission, but differ strongly in terms of parameter estimates, parameter uncertainties and computational intensity. Discretization in the frequency domain (as in the discrete cosine transform parameterization) provides more accurate models at a lower computational cost compared to spatially discretized (Cartesian) models. A priori knowledge about the expected geologic structures allows for non-discretized model descriptions with markedly reduced degrees of freedom. Constraining the solutions to the known injected gas volume improved estimates of saturation and parameter values of the petrophysical relationship. (C) 2014 Elsevier B.V. All rights reserved.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.

Comparative analysis between biogas flow in landfill and electrical resistivity tomography in Rio Claro city, Brazil

The biogas originated from anaerobic degradation of organic matter in landfills consists basically in CH4, CO2, and H2O. The landfills represent an important depository of organic matter with high energetic potential in Brazil, although with inexpressive use in the present. The estimation of production of the productive rate of biogas represents one of the major difficulties of technical order to the planning of capture system for rational consumption of this resource. The applied geophysics consists in a set of methods and techniques with wide use in environmental and hydrogeological studies. The DC resistivity method is largely applied in environmental diagnosis of the contamination in soil and groundwater, due to the contrast of electrical properties frequent between contaminated areas and the natural environment. This paper aims to evaluate eventual relationships between biogas flows quantified in drains located in the landfill, with characteristic patterns of electrical resistivity in depth. The drain of higher flow (117 m3 /h) in depth was characterized for values between 8000 Ω⋅m and 100.000 Ω⋅m, in contrast with values below 2000 Ω⋅m, which characterize in subsurface the drain with less flow (37 m3 /h), besides intermediary flow and electrical resistivity values, attributed to the predominance of areas with accumulation or generation of biogas.

Ground penetrating radar (GPR) and electrical resistivity tomography (ERT) survey on the neotectonic activity of the Padul-Nigüelas Fault Zone (PNFZ), Southern Spain

The Padul-Nigüelas Fault Zone (PNFZ) is situated at the south-western mountain front of the Sierra Nevada (Spain) in an extensive regime and belongs to the internal zone of the Betic Cordilleras. The aim of this study is a collection of new evidence for neotectonic activity of the fault zone with classical geological field work and modern geophysical methods, such as ground penetrating radar (GPR). Among an apparently existing bed rock fault scarp with triangular facets, other evidences, such as deeply incised valleys and faults in the colluvial wedges, are present in the PNFZ. The preliminary results of our recent field work have shown that the synsedimentary faults within the colluvial sediments seem to propagate basinwards and the bed rock fault is only exhumed due to erosion for the studied segment (west of Marchena). We will use further GPR data and geomorphologic indices to gather further evidences of neotectonic activity of the PNFZ.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.

Dynamic Imaging in Electrical Impedance Tomography of the Human Chest With Online Transition Matrix Identification

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

Three-dimensional electrical impedance, tomography: A topology optimization approach

Electrical impedance tomography is a technique to estimate the impedance distribution within a domain, based on measurements on its boundary. In other words, given the mathematical model of the domain, its geometry and boundary conditions, a nonlinear inverse problem of estimating the electric impedance distribution can be solved. Several impedance estimation algorithms have been proposed to solve this problem. In this paper, we present a three-dimensional algorithm, based on the topology optimization method, as an alternative. A sequence of linear programming problems, allowing for constraints, is solved utilizing this method. In each iteration, the finite element method provides the electric potential field within the model of the domain. An electrode model is also proposed (thus, increasing the accuracy of the finite element results). The algorithm is tested using numerically simulated data and also experimental data, and absolute resistivity values are obtained. These results, corresponding to phantoms with two different conductive materials, exhibit relatively well-defined boundaries between them, and show that this is a practical and potentially useful technique to be applied to monitor lung aeration, including the possibility of imaging a pneumothorax.

The borehole-fluid effect in electrical resistivity imaging

Fluid that fills boreholes in crosswell electrical resistivity investigations provides the necessary electrical contact between the electrodes and the rock formation but it is also the source of image artifacts in standard inversions that do not account for the effects of the boreholes. The image distortions can be severe for large resistivity contrasts between the rock formation and borehole fluid and for large borehole diameters. We have carried out 3D finite-element modeling using an unstructured-grid approach to quantify the magnitude of borehole effects for different resistivity contrasts, borehole diameters, and electrode configurations. Relatively common resistivity contrasts of 100:1 and borehole diameters of 10 and 20 cm yielded, for a bipole length of 5 m, apparent resistivity underestimates of approximately 12% and 32% when using AB-MN configurations and apparent resistivity overestimates of approximately 24% and 95% when using AM-BN configurations. Effects are generally more severe at shorter bipole spacings. We report the results obtained by either including or ignoring the boreholes in inversions of 3D field data from a test site in Switzerland, where approximately 10,000 crosswell resistivity-tomography measurements were made across six acquisition planes among four boreholes. Inversions of raw data that ignored the boreholes filled with low-resistivity fluid paradoxically produced high-resistivity artifacts around the boreholes. Including correction factors based on the modeling results fora ID model with and without the boreholes did not markedly improve the images. The only satisfactory approach was to use a 3D inversion code that explicitly incorporated the boreholes in the actual inversion. This new approach yielded an electrical resistivity image that was devoid of artifacts around the boreholes and that correlated well with coincident crosswell radar images.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

(Table T1) Permeability, electrical resistivity and density of ODP Hole 193-1188A and Site 193-1189 minicores, PACMANUS field

Magmatic fluids, heat fluxes, and fluid/rock interactions associated with hydrothermal systems along spreading centers and convergent margins have a significant impact on the genesis of major sulfide deposits and biological communities. Circulation of hydrothermal fluids is one of the most fundamental processes associated with localized mineralization and is controlled by inherent porous and permeable properties of the ocean crust. Heat from magmatic intrusions drives circulation of seawater through permeable portions of the oceanic crust and upper mantle, discharging at the seafloor as both focused high-temperature (250°-400°C) fluids and diffuse lower-temperature (<250°C) fluids. This complex interaction between the circulating hydrothermal fluids and the oceanic basement greatly influences the physical properties and the composition of the crust (Thompson, 1983; Jacobson, 1992, doi:10.1029/91RG02811; Johnson and Semyan, 1994, doi:10.1029/93JB00717). During Ocean Drilling Program (ODP) Leg 193, 13 holes were drilled in the PACMANUS hydrothermal system (Binns, Barriga, Miller, et al., 2002, doi:10.2973/odp.proc.ir.193.2002). The hydrothermal system consists of isolated hydrothermal deposits lined along the main crest of the Pual Ridge, a 500- to 700-m-high felsic neovolcanic ridge in the eastern Manus Basin. The principal drilling targets were the Snowcap (Site 1188) and Roman Ruins (Site 1189) active hydrothermal fields. Samples from these two sites were used for a series of permeability, electrical resistivity, and X-ray computed tomography measurements.

USE OF AN ELECTRICAL IMPEDANCE TOMOGRAPHY METHOD TO DETECT AND TRACK FRACTURES IN A GELATIN MEDIUM

Electrical impedance tomography is applied to the problem of detecting, locating, and tracking fractures in ballistics gelatin. The hardware developed is intended to be physically robust and based on off-the-shelf hardware. Fractures were created in two separate ways: by shooting a .22 caliber bullet into the gelatin and by injecting saline solution into the gelatin. The .22 caliber bullet created an air gap, which was seen as an increase in resistivity. The saline solution created a fluid filled gap, which was seen as a decrease in resistivity. A double linear array was used to take data for each of the fracture mechanisms and a two dimensional cross section was inverted from the data. The results were validated by visually inspecting the samples during the fracture event. It was found that although there were reconstruction errors present, it was possible to reconstruct a representation of the resistive cross section. Simulations were performed to better understand the reconstructed cross-sections and to demonstrate the ability of a ring array, which was not experimentally tested.

Logarithmic contribution to the electrical resistivity in (Ru(1-x)Ir(x))Sr(2)GdCu(2)O(8) compounds

A systematic study of magnetoresistance and dc magnetization was conducted in polycrystalline (Ru(1-x)Ir(x))Sr(2)GdCu(2)O(8) [(Ru,Ir)-1212] compounds, for 0 <= x <= 0.15. We found that a deviation from linearity in the normal-state electrical resistivity (rho) curves for temperatures below the magnetic transition temperature T(M) < 130 K can be properly described by a logarithmic term. The prefactor C(x, H) of this anomalous ln T contribution to rho(T) increases linearly with the Ir concentration, and diminishes rapidly with increasing applied magnetic field up to H approximate to 4 T, merging with the C(0,H) curve at higher magnetic fields. Correlation with magnetic susceptibility measurements supports a scenario of local perturbations in the orientation of Ru moments induced in the neighborhood of the Ir ions, therefore acting as scattering centers. The linear dependence of the prefactor C(x,H=0) and the superconducting transition temperature T(SC) on x points to a common source for the resistivity anomaly and the reduction in T(SC), suggesting that the CuO(2) and RuO(2) layers are not decoupled.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.