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.

Time-lapse geophysical investigations over known archaeological features using electrical resistivity imaging and earth resistance

Electrical methods of geophysical survey are known to produce results that are hard to predict at different times of the year, and under differing weather conditions. This is a problem which can lead to misinterpretation of archaeological features under investigation. The dynamic relationship between a ‘natural’ soil matrix and an archaeological feature is a complex one, which greatly affects the success of the feature’s detection when using active electrical methods of geophysical survey. This study has monitored the gradual variation of measured resistivity over a selection of study areas. By targeting difficult to find, and often ‘missing’ electrical anomalies of known archaeological features, this study has increased the understanding of both the detection and interpretation capabilities of such geophysical surveys. A 16 month time-lapse study over 4 archaeological features has taken place to investigate the aforementioned detection problem across different soils and environments. In addition to the commonly used Twin-Probe earth resistance survey, electrical resistivity imaging (ERI) and quadrature electro-magnetic induction (EMI) were also utilised to explore the problem. Statistical analyses have provided a novel interpretation, which has yielded new insights into how the detection of archaeological features is influenced by the relationship between the target feature and the surrounding ‘natural’ soils. The study has highlighted both the complexity and previous misconceptions around the predictability of the electrical methods. The analysis has confirmed that each site provides an individual and nuanced situation, the variation clearly relating to the composition of the soils (particularly pore size) and the local weather history. The wide range of reasons behind survey success at each specific study site has been revealed. The outcomes have shown that a simplistic model of seasonality is not universally applicable to the electrical detection of archaeological features. This has led to the development of a method for quantifying survey success, enabling a deeper understanding of the unique way in which each site is affected by the interaction of local environmental and geological conditions.

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.

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.

Electrical source imaging for presurgical focus localization in epilepsy patients with normal MRI.

PURPOSE: Patients with magnetic resonance (MR)-negative focal epilepsy (MRN-E) have less favorable surgical outcomes (between 40% and 70%) compared to those in whom an MRI lesion guides the site of surgical intervention (60-90%). Patients with extratemporal MRN-E have the worst outcome (around 50% chance of seizure freedom). We studied whether electroencephalography (EEG) source imaging (ESI) of interictal epileptic activity can contribute to the identification of the epileptic focus in patients with normal MRI. METHODS: We carried out ESI in 10 operated patients with nonlesional MRI and a postsurgical follow-up of at least 1 year. Five of the 10 patients had extratemporal lobe epilepsy. Evaluation comprised surface and intracranial EEG monitoring of ictal and interictal events, structural MRI, [(18)F]fluorodeoxyglucose positron emission tomography (FDG-PET), ictal and interictal perfusion single photon emission computed tomography (SPECT) scans. Eight of the 10 patients also underwent intracranial monitoring. RESULTS: ESI correctly localized the epileptic focus within the resection margins in 8 of 10 patients, 9 of whom experienced favorable postsurgical outcomes. DISCUSSION: The results highlight the diagnostic value of ESI and encourage broadening its application to patients with MRN-E. If the surface EEG contains fairly localized spikes, ESI contributes to the presurgical decision process.

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.

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.

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.

The effects of moisture, cation concentration, temperature, density and composition of soils on their electrical resistivity

The effects. of moisture, cation concentration, dens ity , temper~ t ure and grai n si ze on the electrical resistivity of so il s are examined using laboratory prepared soils. An i nexpen si ve method for preparing soils of different compositions was developed by mixing various size fractions i n the laboratory. Moisture and cation c oncentration are related to soil resistivity by powe r functions, whereas soil resistiv ity and temperature, density, Yo gravel, sand , sil t, and clay are related by exponential functions . A total of 1066 cases (8528 data) from all the experiments were used in a step-wise multiple linear r egression to determine the effect of each variable on soil resistivity. Six variables out of the eight variables studied account for 92.57/. of the total variance in so il resistivity with a correlation coefficient of 0.96. The other two variables (silt and gravel) did not increase the · variance. Moisture content was found to be - the most important Yo clay. variable- affecting s oil res istivi ty followed by These two variables account for 90.81Yo of the total variance in soil resistivity with a correlation ~oefficient ·.of 0 . 95. Based on these results an equation to ' ~~ed{ ct soil r esist ivi ty using moisture and Yo clay is developed . To t est the predicted equation, resistivity measurements were made on natural soils both in s i tu a nd i n the laboratory. The data show that field and laboratory measurements are comparable. The predicted regression line c losely coinciqes with resistivity data from area A and area B soils ~clayey and silty~clayey sands). Resistivity data and the predicted regression line in the case of c layey soils (clays> 40%) do not coincide, especially a t l ess than 15% moisture. The regression equation overestimates the resistivity of so i l s from area C and underestimates for area D soils. Laboratory prepared high clay soils give similar trends. The deviations are probably caused by heterogeneous distribution of mo i sture and difference in the type o f cl ays present in these soils.

Resistivity imaging survey of the Roman barrows at Bartlow, Cambridgeshire, UK

Resistivity imaging was carried out on four large Roman barrows at Bartlow in Cambridgeshire. The geophysical survey formed part of a wider research project designed to record and assess the landscape context of the largest surviving Roman burial mounds in Britain. The barrows today range in height from 6.6 m to 13.2 m and their steep profile loosed particular practical and modelling challenges. Data were obtained using a Campus Geopulse resistance meter with up to 50 electrodes spaced at 1 m intervals and lines up to 76 m long. A total of 24 lines was obtained. Topographic corrections were applied to the pseudosections, whichwere inverted using Res 2 Dinv and Res3 Dinv. Resistivity imaging was particularly successful in identifying evidence for the antiquarian explorations of the site. Central collapse features or in-filled tunnels image as high resistance features in all barrows and in one (Barrow IV) there is also a low resistance feature in the approximate position of a known antiquarian tunnel. Barrow VI had a thick covering of high-resistivity that may relate to nineteenth century landscaping and reconstruction of this monument. Resistivity imaging also revealed possible evidence for ancient revetments in all four large barrows. Copyright (c) 2006 John Wiley & Sons, Ltd.

Use of electrical resistivity to identify collapsible soils in Brazil

Three soil profiles in Ilha Solteira, Brazil were investigated to establish their potential for collapsible behavior. The soil profiles were identified using terrain evaluation techniques and simple laboratory tests. Geophysical surveys were undertaken as they are quick and relatively cheap. The results were correlated with trial pit descriptions and cone and standard penetration tests. The study has shown that electrical resistivity is a useful tool for the preliminary identification of horizons of collapsible soils, before more expensive intrusive and laboratory work is undertaken.

Relaxation of the electrical resistivity in Cr-doped Nd(0.5)Ca(0.5)MnO(3) single crystals

We have performed a systematic study of the time and temperature dependencies of the electrical resistivity (rho(T, t)) inNd(0.5)Ca(0.5)Mn(1-x)Cr(x)O(3) single crystals with x = 0.02 and 0.07 in order to examine the dynamics of the phase separation. The relaxation effects can be described by the combination of a rapid exponential increase/decrease with a slower logarithmic contribution at longer times. The experimental results suggest the existence of a large temperature window in which huge relaxation effects occur, and the relative fraction of the coexisting phases rapidly changes as a function of time, depending on the initial magnetic state of the sample. The rho(T, t) relaxation measurements were shown to be a suitable tool for probing the dynamical nature of the phase separation, in which magnetically distinct phases compete against each other in a wide temperature range. In addition, the features observed in the rho(T, t) curves were found to be in excellent agreement with both the magnetic properties and the structural transitions observed in these manganites.

Use of electrical resistivity to identify collapsible soils in Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)