scispace - formally typeset
Search or ask a question
Journal ArticleDOI

2.5D resistivity modeling of embankment dams to assess influence from geometry and material properties

02 Jun 2006-Geophysics (Society of Exploration Geophysicists)-Vol. 71, Iss: 3, pp 107-114
TL;DR: In this article, the influence from 3D effects created by specific dam geometry and effects of water level fluctuations in the reservoir was evaluated by modeling two rockfill embankment dams with central till cores in the north of Sweden.
Abstract: Repeated resistivity measurement is a potentially powerful method for monitoring development of internal erosion and anomalous seepage in earth embankment dams. This study is part of a project to improve current longterm monitoring routines and data interpretation and increasing the understanding when interpreting existing data. This is accomplished by modeling various occurrences typical of embankment structures using properties from two rockfill embankment dams with central till cores in the north of Sweden. The study evaluates the influence from 3D effects created by specific dam geometry and effects of water level fluctuations in the reservoir. Moreover, a comparison between different layout locations is carried out, and detectability of internal erosion scenarios is estimated through modeling of simulated damage situations. Software was especially developed to model apparent resistivity for geometries and material distributions for embankment dams. The model shows that the 3D effect from the embankment geometry is clearly significant when measuring along dam crests. For dams constructed with a conductive core of fine-grained soil and high-resistive rockfill, the effect becomes greatly enhanced. Also, water level fluctuations have a clear effect on apparent resistivities. Only small differences were found between the investigated arrays. A layout along the top of the crest is optimal for monitoring on existing dams, where intrusive investigations are normally avoided, because it is important to pass the current through the conductive core, which is often the main target of investigation. The investigation technique has proven beneficial for improving monitoring routines and increasing the understanding of results from the ongoing monitoring programs. Although the technique and software are developed for dam modeling, it could be used for estimation of 3D influence on any elongated structure with a 2D cross section.

Summary (3 min read)

INTRODUCTION

  • Internal erosion is one of the major causes of embankment dam failures.
  • This method has been shown to be effective in revealing information about conditions in the core itself.
  • This means that application of standard 2D techniques on embankment dams with measurement layouts along the crest of the dam cannot be used without cau-tion because of the obvious 3D effects from the dam geometry.
  • The study covered several situations and scenarios essential for interpreting and evaluating data from resistivity measurements on embankment dams.

Software description

  • Software written for 2D resistivity/IP modeling was modified to simulate a dam-monitoring survey by allowing dam geometries in the 2D-model parameterization and a 3D measurement, which means that the current injection and potential pickup may be at any point in the dam.
  • Assumed resistivities must be constant in the electrode-layout direction, i.e., along the dam, and variable in the dam cross section, whereas the electrodes can be placed anywhere in all three dimensions.
  • Such 2.5D modeling is simply accomplished by involving the inverse Fourier transform for an electrode array parallel to the strike direction ͑Dey and Morrison, 1979a, b; Queralt et al., 1991͒.
  • The software uses the finite-element method because this method makes it easier to deal with the dam geometry, compared to the finite-difference method.
  • The authors compared the results with different element sizes and wavenumber sampling schemes.

Model geometry, material properties, and damage types

  • The dam model is a zoned embankment dam with a central till core, surrounding filter zones, and support rockfill ͑Figure 1͒.
  • Because of difficulties in estimating electrical properties of involved materials and lack of appropriate data in literature, some uncertainties are connected to these parameters.
  • For this study, the core resistivity was estimated from existing monitoring data from two Swedish dams ͑Johansson et al., 2000͒ together with laboratory resistivity measurements of similar till samples ͑Bergström, 1998͒ -even though an unsatisfying variation was found in this data.
  • Damaged zones often have this kind of extended shape because the dam is constructed in layers.
  • A resistivity increase of five times in the core was assumed because of internal erosion.

Modeling strategies

  • To evaluate responses from different electrode arrays, four arrays were selected for all modeling situations.
  • An electrode spacing of 5 m was selected for the dam model because that gives a reasonable relation between electrode spacing and dam height similar to what could be expected in an actual in situ situation.
  • All combinations, including a-spacings from one to seven ͑multiples of five͒ and n-factors ͑one to six͒, were used for the calculations.
  • Of the four examined arrays, dipole-dipole is by its nature most different from the others, and in some situations, it gave responses that were different than the others.
  • Only when examining special cases, such as cylindrical damages or elongated damage zones with lim-.

3D effects

  • The 3D effects and their dependency on material parameters were examined for a dam with the model cross section described in Figure 1 .
  • The effects were estimated by comparing the responses from two models: a 2.5D model and a 1D model with the properties of the model midsection, i.e., the section with the electrode layout extended to horizontal layers.
  • Sample results for the dipole-dipole and the Schlumberger arrays are shown in Figure 2 .
  • Next, the dependency of input-material parameters was similarly evaluated using a model with constant resistivity for the whole dam cross section, including the reservoir water.
  • It is obvious that most of the huge 3D effect arises from the contrast between the relatively conductive core and the high resistivity of the main part of the dam cross section.

Reservoir-level fluctuations

  • The effect of lowering the reservoir was examined, using the dam model in Figure 1 .
  • 3D effects estimated as relation between 1D and 2.5D models with assumed material properties for the modeled cross section and reservoir.
  • For both arrays, a-spacing is the spacing between potential electrodes, and n-factor is the shortest distance between potential and current electrode divided by the a-spacing.
  • The calculations were made once for each depth.
  • For the large lowering of the reservoir, the same effect was estimated to be moving toward approximately 40% ͑1.40 times͒ for the largest electrode distances.

Detectability of internal erosion zones

  • When internal erosion occurs, the material properties of the eroded zone will change as porosity increases and fines are washed away.
  • A permanent or possibly semipermanent change ͑because it may heal by itself͒ in the resistivity characteristics of the dam core will occur.
  • To estimate the imaging potential of the damages, standard 1D, multilayer, smooth inversion ͑Auken et al., 2004͒ was carried out on the forward model responses.
  • The anomaly effect is enhanced through inversion, but effects from the dam geometry cause the damage to localize at a shallower level than the real case ͑Figure 9͒.
  • It is not likely that the damages would be detected by a single survey, but with repeated measurements the possibilities would be fair.

Comparison of different layout locations

  • Modeling of different layout placements is helpful for interpreting data from Swedish dam monitoring, especially at the Hällby Dam, where layouts are not only placed along the crest but also on a line along the upstream and the downstream side ͑Johansson et al., 2000͒.
  • All of them are placed directly beneath the surface of the dam.
  • For the layouts along the upstream toe and the mid-upstream slope, the upstream electrodes are placed below the water table.
  • The calculated-anomaly effects are less than 1% ͑Ͻ1.01 times͒ for all different placements of the layouts, re- gardless of the damage location.
  • Obviously, the channeling effect that concentrates current flow within the conductive dam core is an important factor.

DISCUSSION AND CONCLUSIONS

  • Resistivity measurements on embankment dam geometries are influenced by many factors, such as effects caused by the geometry and variation in material properties across the dam cross section, impact of water-level changes, and electrode-layout location.
  • The influence is similar for all of the examined arrays, ranging from three to seven times the value of the standard 1D model for the geometry and material properties assumed.
  • Resistivities measured along the dam crest were shown to be significantly influenced by fluctuations in the reservoir level.
  • It is unlikely that such damages could be detected by a single resistivity survey using surface electrodes.
  • Also note that all damage types were shaped as extended layers and that the results may not be fully applicable, for instance, to a cylindrically shaped damage and other damage zones with limited extent along the dam.

Did you find this useful? Give us your feedback

Content maybe subject to copyright    Report

LUND UNIVERSITY
PO Box 117
221 00 Lund
+46 46-222 00 00
2.5D resistivity modeling of embankment dams to assess influence from geometry and
material properties
Sjödahl, Pontus; Dahlin, Torleif; Zhou, Bing
Published in:
Geophysics
DOI:
10.1190/1.2198217
2006
Link to publication
Citation for published version (APA):
Sjödahl, P., Dahlin, T., & Zhou, B. (2006). 2.5D resistivity modeling of embankment dams to assess influence
from geometry and material properties.
Geophysics
,
71
(3), G107-G114. https://doi.org/10.1190/1.2198217
Total number of authors:
3
General rights
Unless other specific re-use rights are stated the following general rights apply:
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors
and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the
legal requirements associated with these rights.
• Users may download and print one copy of any publication from the public portal for the purpose of private study
or research.
• You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal
Read more about Creative commons licenses: https://creativecommons.org/licenses/
Take down policy
If you believe that this document breaches copyright please contact us providing details, and we will remove
access to the work immediately and investigate your claim.

2.5D resistivity modeling of embankment dams to assess
influence from geometry and material properties
Pontus Sjödahl
1
, Torleif Dahlin
1
, and Bing Zhou
2
ABSTRACT
Repeated resistivity measurement is a potentially powerful
method for monitoring development of internal erosion and
anomalous seepage in earth embankment dams. This study is
part of a project to improve current longterm monitoring rou-
tines and data interpretation and increasing the understanding
when interpreting existing data. This is accomplished by mod-
eling various occurrences typical of embankment structures us-
ing properties from two rockfill embankment dams with central
till cores in the north of Sweden. The study evaluates the influ-
ence from 3D effects created by specific dam geometry and ef-
fects of water level fluctuations in the reservoir. Moreover, a
comparison between different layout locations is carried out,
and detectability of internal erosion scenarios is estimated
through modeling of simulated damage situations. Software
was especially developed to model apparent resistivity for ge-
ometries and material distributions for embankment dams. The
model shows that the 3D effect from the embankment geometry
is clearly significant when measuring along dam crests. For
dams constructed with a conductive core of fine-grained soil
and high-resistive rockfill, the effect becomes greatly enhanced.
Also, water level fluctuations have a clear effect on apparent re-
sistivities. Only small differences were found between the in-
vestigated arrays. A layout along the top of the crest is optimal
for monitoring on existing dams, where intrusive investigations
are normally avoided, because it is important to pass the current
through the conductive core, which is often the main target of
investigation. The investigation technique has proven beneficial
for improving monitoring routines and increasing the under-
standing of results from the ongoing monitoring programs. Al-
though the technique and software are developed for dam mod-
eling, it could be used for estimation of 3D influence on any
elongated structure with a 2D cross section.
INTRODUCTION
Internal erosion is one of the major causes of embankment dam
failures. Monitoring systems can significantly improve the safety
of such dams. However, to detect erosion early, monitoring sys-
tems must be highly sensitive and, at the same time, sufficiently
cover the embankment area. In addition, it should be possible to in-
stall such monitoring systems in existing dams, and these systems
should be capable of identifying small seepage changes, as well as
leakage. Experience from research and field installations carried
out in Sweden since 1993 indicates that monitoring systems based
on resistivity measurements may be able to meet this need Johans-
son and Dahlin, 1996; Johansson and Dahlin, 1998; Johansson et
al. 2000. In addition, using a resistivity monitoring technique is
essentially nondestructive. This is particularly important when
working with embankment dams, where drilling and other pen-
etrating investigations are normally avoided.
An electrode layout along the top of the dam core is the most
practical and favorable method of installing resistivity monitoring
systems on existing dams. This will be shown later in the paper.
This method has been shown to be effective in revealing informa-
tion about conditions in the core itself. In addition, good electrode
grounding conditions can be provided in the fine-grained environ-
ment commonly found in the dam core Dahlin et al., 2001. Stan-
dard 2D-inversion schemes are a common technique for processing
data from resistivity profiling Smith and Vozoff, 1984; Tripp et al.,
1984; Li and Oldenburg, 1992; Loke and Barker, 1995; LaBrecque
et al., 1996.
When doing 2D inversion, it is assumed that the properties of
the ground are constant in the third dimension, i.e., the direc-
tion perpendicular to the electrode layout. Deviations from this are
commonly referred to as 3D effects. This means that application of
standard 2D techniques on embankment dams with measurement
layouts along the crest of the dam cannot be used without cau-
Manuscript received by the Editor March 25, 2004; revised manuscript received August 13, 2005; published online May 24, 2006.
1
Lund University, Engineering Geology, Box 118, 221 00 Lund, Sweden. E-mail: pontus.sjodahl@tg.lth.se; torleif.dahlin@tg.lth.se.
2
University of Adelaide, Department of Physics, School of Chemistry & Physics, South Australia 5005, Australis. E-mail: bing.zhou@adelaide.edu.au.
© 2006 Society of Exploration Geophysicists. All rights reserved.
GEOPHYSICS, VOL. 71, NO. 3 MAY-JUNE 2006; P. G107–G114, 9 FIGS., 3 TABLES.
10.1190/1.2198217
G107

tion because of the obvious 3D effects from the dam geometry. It is
possible to use 3D inversion techniques Park and Van, 1991;
Sasaki, 1994; Zhang et al., 1995; Loke and Barker, 1996. How-
ever, they still may not be convenient for repeated measurements,
mainly because of limitations in computational resources and be-
cause data sets are 2D if only measured along the crest. Therefore,
a reasonable approach is to use common 2D techniques and then
estimate the distortions and errors that are induced in the process.
Heretofore, the terms 3D effect refer to the errors received when
measuring along an embankment, assuming standard 2D condi-
tions. The most obvious effect is the embankment topography;
however, the most significant effect might come from the variation
in electrical properties of the construction materials in the zoned
embankment dam.
The aim of this study was to improve current, longterm monitor-
ing routines on two embankment dams in the north of Sweden. The
study covered several situations and scenarios essential for inter-
preting and evaluating data from resistivity measurements on em-
bankment dams. Investigations of these different situations were
carried out through numerical calculations. The influence of the
specific dam geometry and zoned construction materials was in-
vestigated via dedicated, 2.5D software. Effects of reservoir water
level and natural, seasonal resistivity variation in the water were
examined as well. Moreover, a comparison was carried out to de-
termine the differences in the efficiency in detecting seepage zones
for four different electrode arrays.
Much work has been done on resistivity forward modeling in 2D
and 3D using the finite-difference method Mufti, 1976; Dey and
Morrison, 1979a, b; Fox et al., 1980 and the finite-element method
Pridmore et al., 1981; Queralt et al., 1991; Sasaki, 1994; Zhou and
Greenhalgh, 2001. Investigative resistivity surveys on embank-
ments to detect structural defects or anomalous seepage are fairly
widespread Abuzeid, 1994; Engelbert et al., 1997; Titov et al.,
2000; Van Tuyen et al., 2000; Buselli and Lu, 2001; Panthulu et al.,
2001; Voronkov et al., 2004. However, modeling studies to find
out more about typical effects from dam geometries is less com-
mon.
If 3D modeling were to be used for our study, large and compu-
tationally heavy models would have been needed to assess the 3D
effects without influence from the finite length of the model.
Therefore, software capable of handling typical dam geometries
was developed for the numerical calculations. This software is a
useful tool for optimizing the monitoring program design and to
improve the interpretation of collected data. It uses forward model-
ing to find the apparent resistivity distribution in earth embank-
ment dams for a given geometry and measurement layout. Addi-
tionally, it is general and may be utilized for many types of
elongated structures, as long as they can be described with an arbi-
trary although constant geometry in the plane perpendicular to
the electrode layout direction.
NUMERICAL MODELING
Software description
Software written for 2D resistivity/IP modeling was modified to
simulate a dam-monitoring survey by allowing dam geometries in
the 2D-model parameterization and a 3D measurement, which
means that the current injection and potential pickup may be at any
point in the dam. The original 2D software was written for 2D-
resistivity tomography and used the common practical situation,
where resistivity tomographic-imaging surveying is conducted in
the plane perpendicular to the strike direction, allowing arbitrary
variation of resistivity in that plane.
More precisely, the modification of the software was done in two
parts. The first considered adjustments of the elements to fit best
the outline and the inner structure of the dam Figure 1, which was
done by applying the finite-element method Zhou, 1998; Zhou
and Greenhalgh, 1999. The second part regarded the calculation
of the potentials parallel to the strike direction. This was accom-
plished by performing the inverse, Fourier-cosine transform with
nonzero y-coordinate of the potential position, according to the
method described by Queralt et al. 1991.
Hence, the modified software is applicable for modeling of the
resistivity structure with surface profile or crosshole survey. How-
ever, because the current electrodes and the potential measure-
ments must be modeled in 3D for the dam survey, we refer to it as
2.5D modeling. Assumed resistivities must be constant in the
electrode-layout direction, i.e., along the dam, and variable in the
dam cross section, whereas the electrodes can be placed anywhere
in all three dimensions. Such 2.5D modeling is simply accom-
plished by involving the inverse Fourier transform for an electrode
array parallel to the strike direction Dey and Morrison, 1979a, b;
Queralt et al., 1991. The approach is more efficient than a full 3D
model, and for an elongated embankment with constant cross sec-
tion, the drawbacks are moderate. Hence, it is an efficient tool for
assessing 3D effects on 1D and 2D resistivity surveying.
The software uses the finite-element method because this
method makes it easier to deal with the dam geometry, compared to
the finite-difference method. It is valid for calculating potential,
apparent resistivity, or IP responses for a model with arbitrary re-
sistivity distribution in the plane perpendicular to the electrode-
layout direction and for any electrode configurations, e.g., surface,
crosshole, or mise-a-la-masse, off-line and in-line measurements
with pole-pole, pole-dipole, dipole-dipole, Schlumberger, and
mixed arrays Zhou and Greenhalgh, 1999.
The accuracy of 2.5D modeling has been checked by comparing
it with some known analytic solutions Zhou, 1998. It has been
shown that the modeling accuracy mainly depends on the element
size, electrode spacings that give different ranges of the wave-
number, and the wavenumber sampling for accurate inverse-
Fourier transform. To obtain satisfactory results for the dam mod-
eling, we determined the accuracy-control parameters by applying
the dam geometry and the electrode layouts employed in the fol-
lowing simulations. We compared the results with different ele-
ment sizes and wavenumber sampling schemes. We found that the
results showed relative errors less than 1%, using element sizes of
about 1 m and 40 wavenumber sampling points.
Model geometry, material properties,
and damage types
The dam model is a zoned embankment dam with a central till
core, surrounding filter zones, and support rockfill Figure 1. This
is the most common design of large Swedish embankment dams.
Geometry and design values are given in Table 1. The electrode
layout is buried 1 m into the top of the core at the midpoint of the
cross section.
Because of difficulties in estimating electrical properties of in-
volved materials and lack of appropriate data in literature, some
uncertainties are connected to these parameters. Here, the rockfill
G108 Sjödahl et al.

was treated as an insulated matrix with all electrical conduction
concentrated to the pore spaces. Thus, Archie’s law was used using
porosity estimates. However, the porosity estimates are to some
extent uncertain in themselves. Regarding the core, the matrix can
no longer be considered an insulator, and other material models
must be used. For this study, the core resistivity was estimated
from existing monitoring data from two Swedish dams Johansson
et al., 2000 together with laboratory resistivity measurements of
similar till samples Bergström, 1998 even though an unsatis-
fying variation was found in this data.
The resistivity of the filter zones has less influence on the mod-
eling results and was assumed to be somewhere between the resis-
tivity of the core and the rockfill. The resistivity of the reservoir
water was taken from monitoring data Johansson et al., 2000.
Electrical material properties are listed in Table 2. In an interna-
tional perspective, these values are quite high, mainly because of
the high resistivity of the water. Assuming a porosity of approxi-
mately 25% may lead to resistivities of several thousand ohmme-
ters in the saturated rockfill. Keep in mind that the main factor in-
fluencing the results is the relative differences in resistivities for
the involved materials.
The simulated damages were studied for two different depths
Table 3. They could be physically interpreted as damaged layers,
possibly resulting from less compaction at initial construction and
possibly worsened as a consequence of regional piping causing a
transport of fines from the core to the filter and fill. The damages
were extended along the full length of the dam. Damaged zones of-
ten have this kind of extended shape because the dam is con-
structed in layers. Even though an extension along the full length
of the dam is not realistic, simulating these kinds of scenarios still
yields useful information. Furthermore, because of software re-
strictions, the modeled-dam cross section must be identical along
the whole length of the dam. Therefore, for example, it was impos-
sible to simulate a concentrated, cylindrical, damage zone through
the dam.
A resistivity increase of five times in the core was assumed be-
cause of internal erosion. Experiments on similar tills have shown
that resistivity can increase up to 10 times because of removal of
fines under water-saturated conditions Bergström, 1998. How-
ever, this should be handled with care because internal erosion in-
creases porosity, affecting the resistivity in the opposite direction.
The resistivity of the filter and fill was assumed not to change be-
cause of the simulated damages.
Modeling strategies
To evaluate responses from different electrode arrays, four ar-
rays were selected for all modeling situations. The dipole-dipole,
pole-dipole, Wenner-Schlumberger, and gradient arrays were cho-
sen because they have shown robust imaging quality in prior mod-
eling studies Dahlin and Zhou, 2004. An electrode spacing of
5 m was selected for the dam model because that gives a reason-
able relation between electrode spacing and dam height similar to
what could be expected in an actual in situ situation. All combina-
tions, including a-spacings from one to seven multiples of five
and n-factors one to six, were used for the calculations. The total
was 42 individual measurements for each array. Generally, the four
different arrays demonstrated similar responses for the different
modeled situations. This was particularly true for the pole-dipole,
Wenner-Schlumberger, and gradient, which are all geometrically
associated. Of the four examined arrays, dipole-dipole is by its na-
ture most different from the others, and in some situations, it gave
responses that were different than the others. Thus, only results
from dipole-dipole and Wenner-Schlumberger arrays will be pre-
sented.
Certainly, when investigating constant cross sections, i.e., no lat-
eral changes, the differences in the design of the arrays will not
show up fully in the results. Only when examining special cases,
such as cylindrical damages or elongated damage zones with lim-
Figure 1. The modeled cross section geometry. A zoned, rockfill
embankment dam with a central till core and surrounding filter
zones. Electrode layouts and damage zones that are used in the
study are marked out.
Table 1. Dam geometry design parameters (see also
Figure 1).
Dam height 60 m
Crest width 8 m
Upstream and downstream slopes 0.55:1
Distance: Top of core crest 3 m
Distance: Max reservoir level crest 6 m
Core width at top/bottom 4 m/20 m
Filter thickness outside core/top core 4 m/1 m
Table 2. Electrical material properties.
Material Resistivity
m
Core 300
Filter 2000
Upstream fill 4000
Downstream fill 20 000
Reservoir water 550
Damaged core 1500
Table 3. Damage types.
Damage type
Thickness of
damaged layer
Depth from crest to
center of damaged layer
Type 1: Thin seepage
zone layer
2 m 20 m
Type 2: Thin seepage
zone layer
2 m 50 m
2.5D Resistivity modeling of embankment dams G109

ited length, can a full verification of the performance of the differ-
ent arrays be obtained.
RESULTS
3D effects
The 3D effects and their dependency on material parameters
were examined for a dam with the model cross section described in
Figure 1. The effects were estimated by comparing the responses
from two models: a 2.5D model and a 1D model with the proper-
ties of the model midsection, i.e., the section with the electrode
layout extended to horizontal layers. The 2.5D model generated
three to seven times higher responses than the 1D model. Sam-
ple results for the dipole-dipole and the Schlumberger arrays are
shown in Figure 2.
Next, the dependency of input-material parameters was simi-
larly evaluated using a model with constant resistivity for the
whole dam cross section, including the reservoir water. The result-
ing effect, caused by the topography for a homogeneous embank-
ment, gave an increase in resistivity of about 30% 1.30 times for
the 2.5D model Figure 3. It is obvious that most of the huge 3D
effect arises from the contrast between the relatively conductive
core and the high resistivity of the main part of the dam cross sec-
tion. Most of the current flow is concentrated in the core that geo-
metrically constitutes a rather thin sheet Figure 4.
Reservoir-level fluctuations
The effect of lowering the reservoir was examined, using the
dam model in Figure 1. This was done because the reservoir water
Figure 2. 3D effects estimated as relation between 1D and 2.5D models with assumed material properties for the modeled cross section and
reservoir. a Dipole-dipole and b Wenner-Schlumberger arrays with a-spacing of 5–35 m in steps of 5 m and n-factors 1-6. For both ar-
rays, a-spacing is the spacing between potential electrodes, and n-factor is the shortest distance between potential and current electrode di-
vided by the a-spacing.
Figure 3. Purely geometrical 3D effects estimated as relation between 1D and 2.5D models with equal material properties in the whole cross
section and reservoir. a Dipole-dipole and b Wenner-Schlumberger arrays with a-spacing of 5–35 m in steps of 5 m and n-factors 1-6.
G110 Sjödahl et al.

Citations
More filters
Journal ArticleDOI
TL;DR: In this paper, a new approach based on redefining the normalization principle is proposed to derive apparent resistivities from the measured data, which are not contaminated by external 3D effects and yield more reliable results when processed within a 2D conventional inversion scheme.

32 citations


Cites methods from "2.5D resistivity modeling of embank..."

  • ...Thanks to recent improvements in data43 acquisition, ERI is now widely used for dike survey and monitoring (Hennig44 et al., 2005; Sjödahl et al., 2006, 2008; Tsourlos et al., 1999)....

    [...]

DOI
01 Jan 2013
TL;DR: In this paper, the performance of electrical resistivity tomography (ERT) applied to the investigation of seepage in earth dams through a case study was evaluated using three two-dimensional (2D) ERT survey lines deployed on the left abutments, dam crest and downstream shell, respectively.
Abstract: This study was aimed at assessing the performance of electrical resistivity tomography (ERT) applied to the investigation of seepage in earth dams through a case study. Several abnormal leaks appeared on the downstream face of an earth dam after the dam was reconstructed to raise the maximum reservoir water level. A study was conducted to investigate the mechanism of the abnormal leakage with the assistance of ERT. Three two-dimensional (2D) ERT survey lines were deployed on the left abutments, dam crest, and downstream shell, respectively. Periodic measurements were additionally collected on the downstream shell for time-lapse measurements. To gain confidence and avoid over interpretation, the results of 2D ERT were appraised by forward modeling and synthetic inversion. Combining ERT results with geotechnical monitoring data clearly indicates the likely mechanism of abnormal seepage. Time-lapse measurements further support the inducted mechanism. Integration of ERT exploration with time-lapse ERT measurements and geotechnical monitoring data was demonstrated to better understand the possible mechanism of the abnormal seepage.

32 citations


Cites background or methods from "2.5D resistivity modeling of embank..."

  • ...Furthermore, dams often presents difficult site conditions for ERT surveys and their geometry associated with the topography and zoned material property can be a complicating factor in the evaluation process (Sjödahl et al. 2006; Kim et al. 2007)....

    [...]

  • ...In a more similar context, Johansson and Dahlin 1996 and Sjödahl et al. 2006 utilized timelapse ERT to monitor seepage in earth embankments....

    [...]

Journal ArticleDOI
TL;DR: In this paper, a finite-element-based method was proposed for the fast nonlinear inversion of large 3D geoelectric data sets acquired around isolated objects or across the earth's surface.
Abstract: We introduce a new finite-element-based scheme for the fast nonlinear inversion of large 3D geoelectric data sets acquired around isolated objects or across the earth’s surface. The principal novelty of this scheme is the combination of a versatile finite-element approach with (1) a method involving minimization of an objective function using a conjugate-gradient algorithm that includes an adjoint-field technique for efficiently establishing the objective-function gradient and (2) parabolic interpolation for estimating suitable inversion step lengths. This scheme is capable of handling large volumes of data acquired using diverse electrode configurations located around or across 3D structures. Only three solutions to the forward problem are required for each iteration. Computation of the Jacobian matrix, which might require computers with a large amount of memory, is not necessary. To minimize artificial irregularities in the inverted models, particularly near the electrodes, we smooth the model parameter...

31 citations

Journal ArticleDOI
TL;DR: In this paper, electrical resistivity tomography (ERT) was used to quantify soil water content in a full-scale dike model and the results indicated that ERT may be a feasible method to quantify water content changes under certain conditions, provided an initial water content estimate is available.
Abstract: We explore the viability of electrical resistivity tomography (ERT) to quantify soil water content in a full-scale dike model. While ERT is an established method to detect changes in soil water content, the retrieval of two-dimensional water content distributions from surface-based surveys is a much more challenging application. We conducted a series of experiments on a full-scale dike model comprising rainfall and flooding events. The dike was surveyed with an ERT line down the land-side slope and monitored with permanently installed spatial time-domain reflectometry (TDR) cables for comparison. Water content changes were calculated by applying Archie’s equation and agreed well with comparison measurements from spatial TDR measurements. The saturation exponent in Archie’s equation was estimated by combining spatial TDR and ERT measurements and agreed well with literature data of comparable sands. The choice of the inversion type, e.g., time-lapse inversion or least-squares versus robust constraint, as well as the influence of inversion parameters, such as regularization and the choice of convergence criterion, is discussed in the context of an apparent inversion artefact for the rainfall experiment, using forward modelling of synthetic data sets. The results indicate that ERT may be a feasible method to quantify water content changes under certain conditions, provided an initial water content estimate is available. © 2011 EAGE

31 citations

Journal ArticleDOI
TL;DR: In this paper, an overview of research into geoelectrical monitoring at the British Geological Survey is presented, highlighting the technology and algorithms that have been developed to investigate the dynamics of landslide processes, the physical integrity of vulnerable earth structures, and the state of permafrost.
Abstract: In this overview of research into geoelectrical monitoring at the British Geological Survey, Paul Wilkinson,* Jonathan Chambers, Oliver Kuras, Philip Meldrum and David Gunn highlight the technology and algorithms that have been developed to investigate the dynamics of landslide processes, the physical integrity of vulnerable earth structures, and the thermal state of permafrost.

29 citations


Cites background from "2.5D resistivity modeling of embank..."

  • ...…recent research has shown that volumetric ERT imaging can be successfully applied to the investigation of these structures (Jackson et al., 2002; Sjodahl et al., 2006; Chambers et al., 2008), and that measures of saturation can be derived from temperature-corrected tomographic images (Brunet et…...

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this article, numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays, including pole-pole (PP), pole-dipole (PD), half-Wenner (HW), Wenner-α (WN), Schlumberger (SC), dipole-dipsole (DD), WenNER-β (WB), γ -array (GM), multiple or moving gradient array (GD) and midpoint-potential-referred measurement (MPR) arrays.
Abstract: Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include polepole (PP), pole-dipole (PD), half-Wenner (HW), Wenner-α (WN), Schlumberger (SC), dipole-dipole (DD), Wenner-β (WB), γ -array (GM), multiple or moving gradient array (GD) and midpoint-potential-referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal-to-noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L1-norm) inversion and smoothness-constrained least-squares (L2-norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high-resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi-electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L2-norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.

731 citations

Journal ArticleDOI
TL;DR: In this article, a smoothing-constrained least-squares inversion method is used for the data interpretation and the computing time required by this technique is greatly reduced by using a homogeneous half-space as the starting model so that the Jacobian matrix of partial derivatives can be calculated analytically.
Abstract: Techniques to reduce the time needed to carry out 3D resistivity surveys with a moderate number (25 to 100) of electrodes and the computing time required to interpret the data have been developed. The electrodes in a 3D survey are normally arranged in a square grid and the pole-pole array is used to make the potential measurements. The number of measurements required can be reduced to about one-third of the maximum possible number without seriously degrading the resolution of the resulting inversion model by making measurements along the horizontal, vertical and 45° diagonal rows of electrodes passing through the current electrode. The smoothness-constrained least-squares inversion method is used for the data interpretation. The computing time required by this technique can be greatly reduced by using a homogeneous half-space as the starting model so that the Jacobian matrix of partial derivatives can be calculated analytically. A quasi-Newton updating method is then used to estimate the partial derivatives for subsequent iterations. This inversion technique has been tested on synthetic and field data where a satisfactory model is obtained using a modest amount of computer time. On an 80486DX2/66 microcomputer, it takes about 20 minutes to invert the data from a 7 by 7 electrode survey grid. using the techniques described below, 3D resistivity surveys and data inversion can be carried out using commercially available field equipment and an inexpensive microcomputer.

678 citations


"2.5D resistivity modeling of embank..." refers methods in this paper

  • ...It is ossible to use 3D inversion techniques Park and Van, 1991; asaki, 1994; Zhang et al., 1995; Loke and Barker, 1996 ....

    [...]

Journal ArticleDOI
TL;DR: In this article, a smoothness-constrained least square method is used to produce a 2D subsurface model free of distortions in the apparent resistivity pseudosection caused by the electrode array geometry used.
Abstract: A fast technique for the inversion of data from resistivity tomography surveys has been developed. This technique is based on the smoothness-constrained, least-squares method, and it produces a 2-D subsurface model that is free of distortions in the apparent resistivity pseudosection caused by the electrode array geometry used. A homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. Tests with a variety of models and data from field surveys show that this technique is insensitive to random noise, provided a sufficiently large damping factor is used, and that it can resolve structures that cause overlapping anomalies in the pseudosection. On a 33 MHz 80486DX microcomputer, it takes about 5 s to process a single data set.

568 citations


"2.5D resistivity modeling of embank..." refers methods in this paper

  • ...…Geophysics, 61, 538–548. i, Y. G., and D. W. Oldenburg, 1992, Approximate inverse mapping in DC resistivity problems: Geophysical Journal International, 109, 343–362. oke, M. H., and R. D. Barker, 1995, Least-squares deconvolution of ap- parent resistivity pseudosections, Geophysics, 60, 1682–1690....

    [...]

  • ...Stanard 2D-inversion schemes are a common technique for processing ata from resistivity profiling Smith and Vozoff, 1984; Tripp et al., 984; Li and Oldenburg, 1992; Loke and Barker, 1995; LaBrecque t al., 1996 ....

    [...]

Journal ArticleDOI
TL;DR: In this paper, an Occam's inversion algorithm for crosshole resistivity data that uses a finite-element method forward solution is discussed, where the earth is discretized into a series of parameter blocks, each containing one or more elements.
Abstract: An Occam's inversion algorithm for crosshole resistivity data that uses a finite-element method forward solution is discussed. For the inverse algorithm, the earth is discretized into a series of parameter blocks, each containing one or more elements. The Occam's inversion finds the smoothest 2-D model for which the Chi-squared statistic equals an a priori value. Synthetic model data are used to show the effects of noise and noise estimates on the resulting 2-D resistivity images. Resolution of the images decreases with increasing noise. The reconstructions are underdetermined so that at low noise levels the images converge to an asymptotic image, not the true geoelectrical section. If the estimated standard deviation is too low, the algorithm cannot achieve an adequate data fit, the resulting image becomes rough, and irregular artifacts start to appear. When the estimated standard deviation is larger than the correct value, the resolution decreases substantially (the image is too smooth). The same effects are demonstrated for field data from a site near Livermore, California. However, when the correct noise values are known, the Occam's results are independent of the discretization used. A case history of monitoring at an enhanced oil recovery site is used to illustrate problems in comparing successive images over time from a site where the noise level changes. In this case, changes in image resolution can be misinterpreted as actual geoelectrical changes. One solution to this problem is to perform smoothest, but non-Occam's, inversion on later data sets using parameters found from the background data set.

459 citations

Journal ArticleDOI
TL;DR: In this article, a numerical technique is developed to solve the three-dimensional potential distribution about a point source of current located in or on the surface of a half-space containing arbitrary two-dimensional conductivity distribution.
Abstract: A numerical technique is developed to solve the three-dimensional potential distribution about a point source of current located in or on the surface of a half-space containing arbitrary two-dimensional conductivity distribution. Finite difference equations are obtained for Poisson's equations by using point- as well as area-discretization of the subsurface. Potential distributions at all points in the set defining the half-space are simultaneously obtained for multiple point sources of current injection. The solution is obtained with direct explicit matrix inversion techniques. An empirical mixed boundary condition is used at the “infinitely distant” edges of the lower half-space. Accurate solutions using area-discretization method are obtained with significantly less attendant computational costs than with the relaxation, finite-element, or network solution techniques for models of comparable dimensions.

407 citations


"2.5D resistivity modeling of embank..." refers methods in this paper

  • ...Such 2.5D modeling is simply accomlished by involving the inverse Fourier transform for an electrode rray parallel to the strike direction Dey and Morrison, 1979a, b; ueralt et al., 1991 ....

    [...]

Frequently Asked Questions (1)
Q1. What contributions have the authors mentioned in the paper ".5d resistivity modeling of embankment dams to assess nfluence from geometry and material properties" ?

This study is part of a project to improve current longterm monitoring routines and data interpretation and increasing the understanding when interpreting existing data. The study evaluates the influence from 3D effects created by specific dam geometry and effects of water level fluctuations in the reservoir.