ELECTRICAL & ELECTROMAGNETIC

EXPLORATION

GPH 231

Abdullah M. Al-Amri
Dept. of Geology & Geophysics
King Saud University, Riyadh
[email protected]
www.a-alamri.com
 ELECTRICAL AND ELECTROMAGNETIC EXPLORATION

PRINCIPLES OF ELECTRICAL & EM ( 2 weeks )

 Electrical properties of rocks
 Mechanism of electrical conduction in materials
 Representative resistivity values
 Conductivity mechanism

FUNDAMENTALS OF CURRENT FLOW ( 2 weeks )

 Fundamentals of the current flow in the earth.
 Potential distribution in a Homogeneous Medium
 Apparent and true resistivity
 Potential and current distribution across boundary

D.C. RESISTIVITY METHOD ( 4 weeks )

 Electrode configurations
 Electric sounding & Electric profiling
 field procedures
 Applications & Ambiguities
 Qualitative & Quantitative Interpretation
 Mise – A- la- Masse Method

ELECTROCHEMICAL METHODS ( 2 weeks )

 self-potential method
 induced polarization method .

MIDTERM EXAM

ELECTROMAGNETIC METHODS ( 4 weeks )

 Classification of electromagnetic systems

 Principles of electromagnetics
 Magnetotelluric Methods
 Vertical loop (VLEM)
 Slingram & Turam Systems
 Very Low Frequency (VLF)
 Audio Frequency Magnetics (AFMAG)
 Time-Domain systems ( TDEM )
 Airborne Method
 Ground Penetrating Radar

FINAL EXAM
HOMEWORK ASSIGNMENTS IN PAGE 97: 1-2-3-5-7-9-13-14-19-21-23-24

2
GRADING :

Midterm exam. 25 %
Lab. 20 %
Homework Assignments 15 %
Final exam. 40 %

TEXT :

Robinson & Coruh (1988 ) . Basic Exploration Geophysics. John Wiley & Sons

Lowrie, W. ( 1997). Fundamentals of Geophysics. Cambridge University Press.

INSTRUCTOR : ABDULLAH M. S. AL-AMRI

OFFICE HOURS : Sun & Tues 1 - 2

Useful Web Pages

 Introduction to Geophysical Exploration: Colorado School of Mines

 World Data Center A
 The Environmental and Engineering Geophysical Society
 Society of Exploration Geophysicists
 American Geophysical Union
U. S. Geological Survey: Geophysics Products Page

3
 ELECTRICAL RESISTIVITY TECHNIQUES

Geophysical resistivity techniques are based on the response of the earth to

the flow of electrical current. In these methods, an electrical current is passed
through the ground and two potential electrodes allow us to record the resultant
potential difference between them, giving us a way to measure the electrical
impedance of the subsurface material. The apparent resistivity is then a function of
the measured impedance (ratio of potential to current) and the geometry of the
electrode array. Depending upon the survey geometry, the apparent resistivity data
are plotted as 1-D soundings, 1-D profiles, or in 2-D cross-sections in order to look
for anomalous regions.
In the shallow subsurface, the presence of water controls much of the
conductivity variation. Measurement of resistivity (inverse of conductivity) is, in
general, a measure of water saturation and connectivity of pore space. This is
because water has a low resistivity and electric current will follow the path of least
resistance. Increasing saturation, increasing salinity of the underground water,
increasing porosity of rock (water-filled voids) and increasing number of
fractures (water-filled) all tend to decrease measured resistivity. Increasing
compaction of soils or rock units will expel water and effectively increase
resistivity. Air, with naturally high resistivity, results in the opposite response
compared to water when filling voids. Whereas the presence of water will reduce
resistivity, the presence of air in voids should increase subsurface resistivity.
Resistivity measurements are associated with varying depths depending on
the separation of the current and potential electrodes in the survey, and can be
interpreted in terms of a lithologic and/or geohydrologic model of the subsurface.
Data are termed apparent resistivity because the resistivity values measured are

4
actually averages over the total current path length but are plotted at one depth
point for each potential electrode pair. Two dimensional images of the subsurface
apparent resistivity variation are called pseudosections. Data plotted in cross-
section is a simplistic representation of actual, complex current flow paths.
Computer modeling can help interpret geoelectric data in terms of more accurate
earth models.

Geophysical methods are divided into two types : Active and Passive

Passive methods (Natural Sources): Incorporate measurements of natural

occurring fields or properties of the earth. Ex. SP, Magnetotelluric (MT), Telluric,
Gravity, Magnetic.

Active Methods (Induced Sources) : A signal is injected into the earth and then
measure how the earth respond to the signal. Ex. DC. Resistivity, Seismic
Refraction, IP, EM, Mise-A-LA-Masse, GPR.

 DC Resistivity - This is an active method that employs measurements of

electrical potential associated with subsurface electrical current flow
generated by a DC, or slowly varying AC, source. Factors that affect the
measured potential, and thus can be mapped using this method include the
presence and quality of pore fluids and clays. Our discussions will focus
solely on this method.
 Induced Polarization (IP) - This is an active method that is commonly
done in conjunction with DC Resistivity. It employs measurements of the
transient (short-term) variations in potential as the current is initially
applied or removed from the ground. It has been observed that when a
current is applied to the ground, the ground behaves much like a capicitor,
storing some of the applied current as a charge that is dissipated upon

5
 removal of the current. In this process, both capacity and electrochemical
effects are responsible. IP is commonly used to detect concentrations of
clay and electrically conductive metallic mineral grains.
 Self Potential (SP) - This is a passive method that employs measurements
of naturally occurring electrical potentials commonly associated with the
weathering of sulfide ore bodies. Measurable electrical potentials have
also been observed in association with ground-water flow and certain
biologic processes. The only equipment needed for conducting an SP
survey is a high-impedance voltmeter and some means of making good
electrical contact to the ground.
 Electromagnetic (EM) - This is an active method that employs
measurements of a time-varying magnetic field generated by induction
through current flow within the earth. In this technique, a time-varying
magnetic field is generated at the surface of the earth that produces a
time-varying electrical current in the earth through induction. A receiver
is deployed that compares the magnetic field produced by the current-
flow in the earth to that generated at the source. EM is used for locating
conductive base-metal deposits, for locating buried pipes and cables, for
the detection of unexploded ordinance, and for near-surface geophysical
mapping.
 Magnetotelluric (MT) - This is a passive method that employs
measurements of naturally occurring electrical currents, telluric currents,
generated by magnetic induction of electrical currents in the ionosphere.
This method can be used to determine electrical properties of materials at
relatively great depths (down to and including the mantle) inside the
Earth. In this technique, a time variation in electrical potential is
measured at a base station and at survey stations. Differences in the

6
 recorded signal are used to estimate subsurface distribution of electrical
resistivity.

Position of Electrical Methods in:

(1) Petroleum Exploration.

The most prominent applications of electrical techniques in petroleum expl.
Are in well logging. Resistivity and SP are standard Logging techniques.
The magnetotelluric method has found important application for pet.
Exploration. In structurally complex region (EX. Rocky Mountains).

(2) Engineering & Groundwater.

D C. Resistivity and EM have found broad use in civil Engineering and
groundwater studies. Saturated / Unsaturated, Saltwater / freshwater

(3) Mineral Exploration.

Electrical methods interpretation difficult below 1000 to 1500 ft. Electrical
exploration methods are the dominant geophysical tools in Mineral Expl.

7
 Ohm’s Law

Ohm’s Law describes the electrical properties of any medium. Ohm’s Law,
V = I R, relates the voltage of a circuit to the product of the current and the
resistance. This relationship holds for earth materials as well as simple circuits.
Resistance( R), however, is not a material constant. Instead, resistivity is an
intrinsic property of the medium describing the resistance of the medium to the
flow of electric current.
Resistivity ρ is defined as a unit change in resistance scaled by the ratio of a
unit cross-sectional area and a unit length of the material through which the current
is passing (Figure 1). Resistivity is measured in ohm-m or ohm-ft, and is the
reciprocal of the conductivity of the material. Table 1 displays some typical
resistivities.

8
 Note that, in Table 1, the resistivity ranges of different earth materials
overlap. Thus, resistivity measurements cannot be directly related to the type of soil
or rock in the subsurface without direct sampling or some other geophysical or
geotechnical information. Porosity is the major controlling factor for changing
resistivity because electricity flows in the near surface by the passage of ions
through pore space in the subsurface materials. The porosity (amount of pore
space), the permeability (connectivity of pores), the water (or other fluid) content
of the pores, and the presence of salts all become contributing factors to changing
resistivity. Because most minerals are insulators and rock composition tends to
increase resistivity, it is easier to measure conductive anomalies than resistive ones
in the subsurface. However, air, with a theoretical infinite resistivity, will produce
large resistive anomalies when filling subsurface voids.

Electric circuit has three main properties:

o Resistance (R): resistance to movement of charge

o Capacitance (C): ability to store charge

o Inductance (L): ability to generate current from changing magnetic

field arising from moving charges in circuit

9
Resistance is NOT a fundamental characteristic of the metal in the wire.

10
 MECHANISM OF ELECTRICAL CONDUCTION

Mechanism of electrical conduction in Materials the conduction of

electricity through materials can be accomplished by three means :
a) The flow of electrons Ex. In Metal
b) The flow of ions Ex. Salt water .
c) Polarization in which ions or electrons move only a short distance under
the influence of an electric field and then stop.

1 Metals :
Conduction by the flow of electrons depends upon the availability of free
electrons. If there is a large number of free electrons available, then the
material is called a metal, the number of free electrons in a metal is roughly
equal to the number of atoms.
The number of conduction electrons is proportional to a factor

n ≈ ε E/KT E ∞ 1/n T∞n

ε : Dielectric constant
K: Boltzman’s constant
T: Absolute Temperature.
E Activation Energy.

Metals may be considered a special class of electron semi conductor for

which E approaches zero.
Among earth materials native gold and copper are true metals. Most sulfide
ore minerals are electron semi conductors with such a low activation energy.

b) The flow of ions, is best exemplified by conduction through water,

especially water with appreciable salinity. So that there is an abundance of
free ions.
Most earth materials conduct electricity by the motion of ions contained in
the water within the pore spaces .

11
There are three exceptions :

1) The sulfide ores which are electron semi conductors.

2) Completely frozen rock or completely dry rock.
3) Rock with negligible pore spaces ( Massive lgneous rooks like gabbro . It
also include all rocks at depths greater than a few kilometers, where pore
spaces have been closed by high pressure, thus studies involving
conductivity of the deep crust and mantle require other mechanisms than
ion flow through connate water.

c) Polarization of ions or sometimes electrons under the influence of an

electrical field, they move a short distance then stop. Ex. Polarization of the
dielectric in a condenser polarization ( electrical moment / unit volume)

12
 Conductivity mechanism in non-water-bearing rocks

1) Extrinsic conductivity for low temperatures below 600-750o k.

2) Intrinsic conductivity for high temperatures.

Most electrical exploration will be concerned only with temperatures well below
600-750o . The extrinsic is due to weakly bonded impurities or defects in the
crystal . This is therefore sensitive to the structure of the sample and to its
thermal history .

Both of these types of conductivity present the same functional form, hence
conductivity vs. temperature for semi conductors can be written :

σ = Ai ε – Ei/RT + Ae ε – Ee/RT

Ai and Ae : Numbers of ions available . Ai is 105 times Ae

Ei and Ee are the activation energies . Ei is 2 times as large as Ee .

R: Boltzman’s constant

13
 ELECTRICAL PROPERTIES OF ROCKS :

 Resistivity (or conductivity), which governs the amount of current that passes
when a potential difference is created.

 Electrochemical activity or polarizability, the response of certain minerals to

electrolytes in the ground, the bases for SP and IP.

 Dielectric constant or permittivity. A measure of the capacity of a material to

store charge when an electric field is applied . It measure the polarizability of a
material in an electric field = 1 + 4 π X
X : electrical susceptibility .

Electrical methods utilize direct current or Low frequency alternating current to

investigate electrical properties of the subsurface.
Electromagnetic methods use alternating electromagnetic field of high frequencies.
Two properties are of primary concern in the Application of electrical methods.
(1) The ability of Rocks to conduct an electrical current.
(2) The polarization which occurs when an electrical current is passed through
them (IP).

Resistivity

For a uniform wire or cube, resistance is proportional to length and inversely

proportional to cross-sectional area. Resistivity is related to resistance but it not
identical to it. The resistance R depends an length, Area and properties of the
material which we term resistivity (ohm.m) .

Constant of proportionality is called Resistivity :

Resistivity is the fundamental physical property of the metal in the wire

14
Resistivity is measured in ohm-m

Conductivity is defined as 1/ρ and is measured in Siemens per meter (S/m),

equivalent to ohm-1m-1.

EX. 1 Copper has ρ =1.7 X 10-8 ohm.m. What is the resistance of 20 m of

copper with a cross-sectional radius of 0.005m .

EX. Quartz has ρ = 1 X 1016 ohm.m. What is the resistance at a quartz wire at
the same dimension.

15
Anisotropy : is a characteristic of stratified rocks which is generally more
conducive in the bedding plane. The anisotropy might be find in a schist (micro
anisotropic) or in a large scale as in layered sequence of shale (macro
anisotropic) .

‫ هذا يعني انه لو طبق‬1 ‫و‬1-1 ‫هو النسبة بين الحد القأصى للمقاومية إلى الحد الدأنى ويصل ما بين‬
‫ التيار في اتجاه واحد فان هذا المعامل يقوم بتغير الصفات الخواص الكهربائية للتجاه الخآر‬.

Coefficient of anisotropy λ = ρt / ρl
ρl : Longitudinal Resistivety .
ρt : Transverse Resistivity.

The effective Resistivity depends on whether the current is flowing parallel to

the layering or perpendicular to it .

R1 = ρ1 h1

The total Resistance for the unit column ( T )

T = ∑ ρ1 h1 Transverse unit resistance

The transverse resistivity ρt is defined by .

ρt = T/H H is the total thickness

For current flowing horizontally, we have a parallel circuit. The reciprocal

resistance is S = 1/ R = ∑ hi / ρi Longitudinal unit conductance

Longitudinal resistivity ρl = H / S

A geoelectric unit is characterized by two Parameters :

1) Layer Resistivity ( ρi )
2) Lager Thickness( ti )

16
 Four electrical parameters can be derived for each layer from the respective
resistivity and thickness. There are :

1) Longitudinal conductance SL= h/ρ = h.σ

2) Transverse resistance T = h.ρ
3) Longitudinal resistivity ρl = h/S
4) Transverse resistivity ρt = T/h

Anisotropy = A = Transverse resistivity ρt / Longitudinal resistivity ρl

The sums of all SL ( ∑ hi / ρi ) are called Dar Zarrouk functions.

The sums of all T ( ∑ hi . ρi ) are called Dar Zarrouk variables.

17
 Classification of Materials according to Resistivities Values

a) Materials which lack pore spaces will show high resistivity such as
- massive limestone
- most igneous and metamorphic (granite, basalt)

b) Materials whose pore space lacks water will show high resistivity such
as : - dry sand and gravel
- Ice .

c) Materials whose connate water is clean (free from salinity ) will show
high resistivity such as :
- clean sand or gravel , even if water saturated.

d) most other materials will show medium or low resistivity, especially if

clay is present such as :
- clay soil
- weathered rock.

The presence of clay minerals tends to decrease the Resistivity because :

2) The clay minerals can combine with water .
3) The clay minerals can absorb cations in an exchangeable state on the surface.
4) The clay minerals tend to ionize and contribute to the supply of free ions.

18
Factors which control the Resistivity

(1) Geologic Age

(2) Salinity.
(3) Free-ion content of the connate water.
(4) Interconnection of the pore spaces (Permeability).
(5) Temperature.
(6) Porosity.
(7) Pressure
(8) Depth

Archie’s Law
Empirical relationship defining bulk resistivity of a saturated porous rock. In
sedimentary rocks, resistivity of pore fluid is probably single most important factor
controlling resistivity of whole rock.

Archie (1942) developed empirical formula for effective resistivity of rock:

ρ0 = bulk rock resistivity

ρw = pore-water resistivity
a = empirical constant (0.6 < a < 1)
m = cementation factor (1.3 poor, unconsolidated) < m < 2.2
(good, cemented or crystalline)
φ = fractional porosity (vol liq. / vol rock)

Formation Factor:

Effects of Partial Saturation:

Sw is the volumetric saturation.

19
 n is the saturation coefficient (1.5 < n < 2.5).
 Archie’s Law ignores the effect of pore geometry, but is a reasonable
approximation in many sedimentary rocks

Resistivity survey instruments:

a- High tension battery pack (source of current).

b- Four metal stakes.
c- Milliammeter.
d- Voltmeter.
e- Four reels of insulated cable.

AC is preferred over DC as source of current. The advantage of using AC is that

unwanted potential can be avoided.

Field considerations for DC Resistivity

1- Good electrode contact with the earth

- Wet electrode location.
- Add Nacl solution or bentonite

2- Surveys should be conducted along a straight line whenever possible .

3- Try to stay away from cultural features whenever possible .

- Power lines
- Pipes
- Ground metal fences
- Pumps

20
 Sources of Noise

There are a number of sources of noise that can effect our measurements of voltage
and current.

1- Electrode polarization.
A metallic electrode like a copper or steel rod in contact with an electrolyte
groundwater other than a saturated solution of one of its own salt will
generate a measurable contact potential. For DC Resistivity, use
nonpolarizing electrodes. Copper and copper sulfate solutions are commonly
used.

2- Telluric currents.
Naturally existing current flow within the earth. By periodically reversing
the current from the current electrodes or by employing a slowly varying AC
current, the affects of telluric can be cancelled.

3- Presence of nearby conductors. (Pipes, fences)

Act as electrical shorts in the system and current will flow along these
structures rather than flowing through the earth.

4- Low resistivity at the near surface.

If the near surface has a low resistivity, it is difficult to get current to flow
more deeply within the earth.

5- Near- electrode Geology and Topography

Rugged topography will act to concentrate current flow in valleys and
disperse current flow on hills.

6- Electrical Anisotropy.
Different resistivity if measured parallel to the bedding plane compared to
perpendicular to it .

7- Instrumental Noise .

8- Cultural Feature .

21
Current Flow in Uniform Earth with Two Electrodes

Current injected by electrode at S1 and exits by electrode at S2:

Lines of constant potential (equipotential) are no longer spherical shells, but can be
calculated from expression derived previously.

Current flow is always perpendicular to equipotential lines.

 Where ground is uniform, measured resistivity should not change with

electrode configuration and surface location.

 Where inhomogeneity present, resistivity varies with electrode position.

Computed value is called apparent resistivity ρA.

22
 Current Flow in A Homogeneous Earth

1. Point current Source :

If we define a very thin shell of thickness dr we can define the potential
different dv

dv = I ( R ) = I ( ρ L / A ) = I ( ρ dr / 2π r2 )

To determine V a t a point . We integrate the above eq. over its distance D to

to infinity :

V = I ρ / 2π D

C: current density per unit of cross sectional area :

23
2. Two current electrodes

To determine the current flow in a homogeneous, isotropic earth when we

have two current electrodes. The current must flow from the positive
(source ) to the resistive ( sink ).

The effect of the source at C1 (+) and the sink at C2 (-)

Vp1 = i ρ / 2π r1 + ( - iρ / 2π r2 )

Vp1 = iρ / 2π { 1/ [ (d/2 + x )2 + Z2 ]0.5 - 1 / [ (d/2 - x )2 + Z2 ]0.5 }

24
3. Two potential Electrodes

Vp1 = i ρ / 2π r1 - iρ / 2π r2

Vp2= i ρ / 2π r3 - iρ / 2π r4

Δ V = Vp1 –Vp2 = i ρ / 2π ( 1/r1 – 1 / r2 – 1 / r3 + 1 / r4 )

25
Depth of Current Penetration

Current flow tends to occur close to the surface. Current penetration can be
increased by increasing separation of current electrodes. Proportion of current
flowing beneath depth z as a function of current electrode separation AB:

Example

If target depth equals electrode separation, only 30% of current flows beneath that
level.

 To energize a target, electrode separation typically needs to be 2-3 times its

depth.

 High electrode separations limited by practicality of working with long cable

lengths. Separations usually less than 1 km.

The fraction of the total current (if) penetrating to depth Z for an electrode
separation of d is given by :

if = 2 / π tan -1 ( 2 Z / d )

26
 ELECTRODE CONFIGURATIONS

The value of the apparent resistivity depends on the geometry of the electrode
array used (K factor)

1- Wenner Arrangement
Named after wenner (1916) .
The four electrodes A , M , N , B are equally spaced along a straight line. The
distance between adjacent electrode is called “a” spacing . So AM=MN=NB=
⅓ AB = a.

Ρa= 2 π a V /I

The wenner array is widely used in the western Hemisphere. This array is
sensitive to horizontal variations.

2- Lee- Partitioning Array .

This array is the same as the wenner array, except that an additional potential
electrode O is placed at the center of the array between the Potential electrodes
M and N. Measurements of the potential difference are made between O and
M and between O and N .
Ρa= 4 π a V / I

This array has been used extensively in the past .

3) Schlumberger Arrangement .
This array is the most widely used in the electrical prospecting . Four electrodes
are placed along a straight line in the same order AMNB , but with AB ≥ 5 MN

  AB  2  MN  2 
    
V 2   2  
a      
I  MN 
 

 


This array is less sensitive to lateral variations and faster to use as only
the current electrodes are moved.

27
4. Dipole – Dipole Array .
The use of the dipole-dipole arrays has become common since the 1950’s ,
Particularly in Russia. In a dipole-dipole, the distance between the current
electrode A and B (current dipole) and the distance between the potential
electrodes M and N (measuring dipole) are significantly smaller than the
distance r , between the centers of the two dipoles.

ρa = π [ ( r2 / a ) – r ] v/i

Or . if the separations a and b are equal and the distance between the centers
is (n+1) a then

ρa = n (n+1) (n+2) . π a. v/i

This array is used for deep penetration ≈ 1 km.

Four basic dipole- dipole arrays .
1) Azimuthal

28
 2) Radial
3) Parallel
4) Perpendicular

When the azimuth angle (Ө ) formed by the line r and the current dipole AB = π /2
, The Azimuthal array and parallel array reduce to the equatorial Array.

When Ө = O , the parallel and radial arrays reduce to the polar or axial array .

If MN only is small is small with respect to R in the equatorial array, the system is
called Bipole-Dipole (AB is the bipole and MN is the dipole ), where AB is large
and MN is small.

If AB and MN are both small with respect to R , the system is dipole- dipole

5) Pole-Dipole Array .
The second current electrode is assumed to be a great distance from the
measurement location ( infinite electrode)

ρa = 2 π a n (n+1) v/i

6) Pole – Pole.
If one of the potential electrodes , N is also at a great distance.

29
Ρa= 2 π a V /I

30
31
 REFRACTION OF ELECTRICAL RESISTIVITY

A. Distortion of Current flow

At the boundary between two media of different resistivities the potential remains
continuous and the current lines are refracted according to the law of tangents.

Ρ1 tan Ө2 = Ρ2 tan Ө1

If ρ2 < p1 , The current lines will be refracted away from the Normal. The line of
flow are moved downward because the lower resistivity below the interface results
in an easier path for the current within the deeper zone.

B. Distortion of Potential

Consider a source of current I at the point S in the first layers P1 of Semi infinite
extent. The potential at any point P would be that from S plus the amount reflected
by the layer P2 as if the reflected amount were coming from the image S/

V1 (P) = i ρ1 / 2π [ (1 / r1) + ( K / r2 ) ]

K = Reflection coefficient = ρ2 – ρ1 / ρ2 + ρ1

32
In the case where P lies in the second medium ρ2, Then transmitting light coming
from S. Since only 1 – K is transmitted through the boundary.
The Potential in the second medium is

V2(P) = I ρ2/ 2π [ (1 / r1) – (K / r1) ]

Continuity of the potential requires that the boundary where r1 = r2 , V1(p)

must be equal to V2 ( P).
At the interface r1 = r2 , V1= V2

Method of Images

Potential at point close to a boundary can be found using "Method of Images" from
optics.

In optics:

Two media separated by semi transparent mirror of reflection and

transmission coefficients k and 1-k, with light source in medium 1.
Intensity at a point in medium 1 is due to source and its reflection,
considered as image source in second medium, i.e source scaled by
reflection coefficient k. Intensity at point in medium 2 is due only to
source scaled by transmission coefficient 1-k as light passed through
boundary.

33
 Electrical Reflection Coefficient
Consider point current source and find expression for current potentials in medium
1 and medium 2: Use potential from point source, but 4π as shell is spherical:

Potential at point P in medium 1:

Potential at point Q in medium2:

At point on boundary mid-way between source and its image:

r1=r2=r3=r say. Setting Vp = Vq, and canceling we get:

k is electrical reflection coefficient and used in interpretation

34
The value of the dimming factor , K always lies between ±1

If the second layer is a pure insulator

( ρ2 = ω ) then K=+1

If the second layer is a perfect conductor

( ρ2 = O ) then K=-1

When ρ1 = ρ2 then No electrical boundary Exists and K = O .

35
 SURVEY DESIGN

Two categories of field techniques exist for conventional resistivity analysis

of the subsurface. These techniques are vertical electric sounding (VES), and
Horizontal Electrical Profiling (HEP).

1- Vertical Electrical Sounding (VES) .

The object of VES is to deduce the variation of resistivity with depth below a
given point on the ground surface and to correlate it with the available
geological information in order to infer the depths and resistivities of the
layers present.

In VES, with wenner configuration, the array spacing “a” is increased by

steps, keeping the midpoint fixed (a = 2 , 6, 18, 54…….) .

In VES, with schlumberger, The potential electrodes are moved only

occasionally, and current electrode are systematically moved outwards in
steps
AB > 5 MN.

2- Horizontal Electrical profiling (HEP) .

The object of HEP is to detect lateral variations in the resistivity of the
ground, such as lithological changes, near- surface faults…… .

In the wenner procedurec of HEP , the four electrodes with a definite array
spacing “a” is moved as a whole in suitable steps, say 10-20 m. four
electrodes are moving after each measurement.

36
In the schlumberger method of HEP, the current electrodes remain fixed at
a relatively large distance, for instance, a few hundred meters , and the
potential electrode with a small constant separation (MN) are moved
between A and B .

Multiple Horizontal Interfaces

For Three layers resistivities in two interface case , four possible curve types
exist.

1- Q – type ρ1> ρ2> ρ3

2- H – Type ρ1> ρ2< ρ3
3- K – Type ρ1< ρ2> ρ3
4- A – Type ρ1< ρ2< ρ3

37
In four- Layer geoelectric sections, There are 8 possible relations :

ρ1> ρ2< ρ3< ρ4 HA Type

ρ1> ρ2< ρ3> ρ4 HK Type
ρ1< ρ2< ρ3< ρ4 AA Type
ρ1< ρ2< ρ3> ρ4 AK Type
ρ1< ρ2> ρ3< ρ4 KH Type
ρ1< ρ2> ρ3> ρ4 KQ Type
ρ1> ρ2> ρ3< ρ4 QH Type
ρ1> ρ2> ρ3> ρ4 QQ Type

38
 Quantitative VES Interpretation: Master Curves

Layer resistivity values can be estimated by matching to a set of master curves

calculated assuming a layered Earth, in which layer thickness increases with depth.
(seems to work well). For two layers, master curves can be represented on a single
plot.

Master curves: log-log plot with ρa / ρ1 on vertical axis and a / h on

horizontal (h is depth to interface)

 Plot smoothed field data on log-log graph transparency.

 Overlay transparency on master curves keeping axes parallel.
 Note electrode spacing on transparency at which (a / h=1) to get interface
depth.
 Note electrode spacing on transparency at which (ρa / ρ1 =1) to get resistivity
of layer 1.
 Read off value of k to calculate resistivity of layer 2 from:

39
 Quantitative VES Interpretation: Inversion

Curve matching is also used for three layer models, but book of many more curves.

Recently, computer-based methods have become common:

 forward modeling with layer thicknesses and resistivities provided by user

 inversion methods where model parameters iteratively estimated from data
subject to user supplied constraints

Example (Barker, 1992)

Start with model of as many layers as data points and resistivity equal to measured
apparent resistivity value.

Calculated curve does not match data, but can be perturbed to improve fit.

40
 Applications of Resistivity Techniques

1. Bedrock Depth Determination

Both VES and CST are useful in determining bedrock depth. Bedrock usually more
resistive than overburden. HEP profiling with Wenner array at 10 m spacing and 10
m station interval used to map bedrock highs.

2. Location of Permafrost

Permafrost represents significant difficulty to construction projects due to

excavation problems and thawing after construction.

 Ice has high resistivity of 1-120 ohm-m

3. Landfill Mapping

Resistivity increasingly used to investigate landfills:

 Leachates often conductive due to dissolved salts

 Landfills can be resistive or conductive, depends on contents

41
 Limitations of Resistivity Interpretation

1- Principle of Equivalence.
If we consider three-lager curves of K (ρ1< ρ2> ρ3 ) or Q type (ρ1> ρ2> ρ3)
we find the possible range of values for the product T2= ρ2 h2 Turns out to
be much smaller. This is called T-equivalence. H = thickness, T :
Transverse resistance it implies that we can determine T2 more reliably
than ρ2 and h2 separately. If we can estimate either ρ2 or h2 independently
we can narrow the ambiguity. Equivalence: several models produce the
same results. Ambiguity in physics of 1D interpretation such that
different layered models basically yield the same response.
Different Scenarios: Conductive layers between two resistors, where
lateral conductance (σh) is the same. Resistive layer between two
conductors with same transverse resistance (ρh).

2- Principle of Suppression.
This states that a thin layer may sometimes not be detectable on the field
graph within the errors of field measurements. The thin layer will then be
averaged into on overlying or underlying layer in the interpretation. Thin
layers of small resistivity contrast with respect to background will be missed.
Thin layers of greater resistivity contrast will be detectable, but equivalence
limits resolution of boundary depths, etc.
The detectibility of a layer of given resistivity depends on its relative
thickness which is defined as the ratio of Thickness/Depth.

42
 Comparison of Wenner and Schlumberger

(1) In Sch. MN ≤ 1/5 AB

Wenner MN = 1/3 AB

(2) In Sch. Sounding, MN are moved only occasionally.

In Wenner Soundings, MN and AB are moved after each measurement.

(3) The manpower and time required for making Schlumberger soundings are
less than that required for Wenner soundings.

(4) Stray currents that are measured with long spreads effect
measurements with Wenner more easily than Sch.

(5) The effect of lateral variations in resistivity are recognized and corrected
more easily on Schlumberger than Wenner.

(6) Sch. Sounding is discontinuous resulting from enlarging MN.

Disadvantages of Wenner Array

All electrodes must be moved for each reading .1

Required more field time .2
More sensitive to local and near surface lateral variations .3
Interpretations are limited to simple, horizontally layered structures .4

Advantages of Schlumberger Array

Less sensitive to lateral variations in resistivity .1

Slightly faster in field operation .2
Small corrections to the field data .3

43
 Disadvantages of Schlumberger Array

Interpretations are limited to simple, horizontally layered structures .1

.For large current electrodes spacing, very sensitive voltmeters are required .2

Advantages of Resistivity Methods

Flexible .1
Relatively rapid. Field time increases with depth .2
Minimal field expenses other than personnel .3
Equipment is light and portable .4
Qualitative interpretation is straightforward .5
Respond to different material properties than do seismic and other methods, .6
specifically to the water content and water salinity

Disadvantages of Resistivity Methods

1- Interpretations are ambiguous, consequently, independent geophysical

and geological controls are necessary to discriminate between valid
alternative interpretation of the resistivity data ( Principles of Suppression
& Equivalence)

2- Interpretation is limited to simple structural configurations.

3- Topography and the effects of near surface resistivity variations can mask
the effects of deeper variations.

4- The depth of penetration of the method is limited by the maximum

electrical power that can be introduced into the ground and by the
practical difficulties of laying out long length of cable. The practical
depth limit of most surveys is about 1 Km.

Accuracy of depth determination is substantially lower than with .5

.seismic methods or with drilling

44
 Lateral inhomogeneities in the ground affect resistivity measurements in
different ways: The effect depends on

The size of inhomogeneity with respect to its depth

2. The size of inhomogeneities with respect to the size of electrode array
3. The resistivity contrast between the inhomogeneity and the surrounding
media
4. The type of electrode array used
5. The geometric form of the inhomogeneity
6. The orientation of the electrode array with respect to the strike of the
inhomogeneity

45
 Mise-A-LA-Masse Method

This is a charged-body potential method is a development of HEP technique

but involves placing one current electrode within a conducting body and the
other current electrode at a semi- infinite distance away on the surface .

This method is useful in checking whether a particular conductive mineral-

show forms an isolated mass or is part of a larger electrically connected ore
body.

There are two approaches in interpretation

1- One uses the potential only and uses the maximum values a being
indicative of the conductive body.

2- The other converts the potential data to apparent resistivity and thus a
high surface voltage manifests itself in a high apparent resistivity

ρa = 4Л X V/I :

Where X is the distance between C1 and P1.

46
 SELF- POTENTIAL (SP)

SP is called also spontaneous polarization and is a naturally occurring potential

difference between points in the ground. SP depends on small potentials or voltages
being naturally produced by some massive ores.
It associate with sulphide and some other types of ores. It works strongly on pyrite,
pyrrohotite, chalcopyrite, graphite.
SP is the cheapest of geophysical methods.

Conditions for SP anomalies

1) Shallow ore body
2) Continuous extension from a zone of oxidizing conditions to one of reducing
conditions, such as above and below water table.

n flow

either or both
pos/neg ions

Note that it is not necessary that an individual ion travel the entire path. Charges can
be exchanged.

47
 The implications of this for potential distribution would be

When we come to consider more specifically the mechanism, we see that it must
be consistent with
 - electron flow in the ore body
 - ion flow in surrounding rock
 - no transfer of ions across ore boundary, although electrons are free to cross

That is we must have

negative ion positive ion

neutral ion neutral ion

reverse process

48
 When we consider the possible ion species, the criteria would be
- common enough
- reversible couple under normal ground conditions
- mobile enough
Sato and Mooney proposed ferric/ferrous couples to satisfy these criteria.

made continuous by O2 – H2 O2
reaction with O2 supplied from
atmosphere

made continuous by reactions

involving ferrous and ferric
hydroxide with presence of H+

49
 Proposed electrochemical mechanism for self-potentials

This proposed mechanism have two geologic implications:

1) The ore body must be an electronic conductor with high conductivity.

This would seem to eliminate sphalerite (zinc sulfide) which has low
conductivity.

2) The ore body must be electrically continuous between a region of oxidizing

conditions and a region of reducing conditions. While water table contact
would not be the only possibility have, it would seem to be a favorable one.

50
 Instrumentation and Field Procedure

Since we wish to detect currents, a natural approach is to measure current.

However, the process of measurement alters the current. Therefore, we arrive at
it though measuring potentials.

Principle, and occasional practice:

More usual practice

Instruments

Equipment:
- potentiometer or high impedance voltmeter
- 2 non-polarizing electrodes
- wire and reel

Non-polarizing electrodes were described in connection with resistivity exploration

although they are not usually required there. Here, they are essential. The use of simple
metal electrodes would generate huge contact or corrosion potentials which would

51
 mask the desired effect. non-polarizing electrodes consist of a metal in contact with a
saturated solution of a salt of the metal . Contact with the earth can be made through a
porous ceramic pot.

The instrument which measures potential difference between the electrodes must have the
following characteristics:
a) capable of measuring +0.1 millivolt,
b) capable of measuring up to ±1000 millivolts (±1 volt)
c) input impedance greater than 10 megaohms, preferably more.

The high input impedance is required in order to avoid drawing current through the
electrodes, whose resistance is usually less than 100 kilohms. In very dry conditions (dry
rock, ice, snow, frozen soil), the electrode resistance may exceed 100 kilohms, in which case
the instrument input impedance should also be increased.

SP are produced by a number of mechanisms :

1) Mineral potential (ores that conduct electronically ) such as most sulphide ores
,Not sphalerite (zinc sulphide) magnetite, graphite. Potential anomaly over sulfide
or graphite body is negative The ore body being a good conductor. Curries current
from oxidizing electrolytes above water – table to reducing one below it .

2)Diffusion potential
RT( Ia – Ic)
Ed == = Ln (C1 / C2 )
n F(Ia+Ic)

Where
Ia , Ic Mobilities of the anions (+ve) and cations( -ve )
R= universal Gas constant ( 8.314JK-1 mol-1 )
T : absolute temperature ( K)
N : is ionic valence
F: Farady’s constant 96487 C mol-1 )
C1 , C2 Solution concentrations .

52
3) Nernst Potential

EN = - ( RT / nF ) Ln ( C1 / C2 )

Where Ia = Ic in the diffusion potential Equation .

4. Streaming potentials due to subsurface water flow are the source of many SP
anomalies. The potential E per unit of pressure drop P (The streaming potentials
coupling cocfficent) is given by :
ε ρ CE δP
EK = 4 π η

ρ Electrical Resistivity of the pore Fluld.

Ek Electro-kinetic potential as a result from an electrolyte flowing through a
porous media.
ε Dielectric constant of the pole fluid.
η Viscosity of the pore fluid
δP pressure difference
CE electro filteration coupling coefficient.

53
 Interpretation
Usually, interpretation consists of looking for anomalies.

The order of magnitude of anomalies is

0-20 mv normal variation

20-50 mv possibly of interest, especially if observed over a fairly large area
over 50 mv definite anomaly
400-1000 mv very large anomalies

Applications

Groundwater applications rely principally upon potential differences produced by

pressure gradients in the groundwater. Applications have included
detection of leaks in dams and reservoirs
location of faults, voids, and rubble zones which affect groundwater flow
delineation of water flow patterns around landslides, wells, drainage structures, and
springs, studies of regional groundwater flow

Other groundwater applications rely upon potential differences produced by

gradients in chemical concentration, Applications have included
outline hazardous waste contaminant plumes

Thermal applications rely upon potential differences produced by temperature

gradients.
Applications have included
geothermal prospecting
map burn zones for coal mine fires
monitor high-temperature areas of in-situ coal gasification processes and oil field
steam and fire floods.

54
 Induced Polarization ( IP)

IP depends on a small amount of electric charge being stored in an ore when a

current is passed Through it , to be released and measured when the current is
switched off .
The main application is in the search for disseminated metallic ores and to a
lesser extent, ground water and geothermal exploration .
Measurements of IP using 2 current electrodes and 2 non-polarizable potential
electrodes. When the current is switched off , the voltage between the potential
electrodes takes a finite to decay to zero because the ground temporarily stores
charge ( become Polarized)

Four systems of IP .
1- Time domain
2- Frequency domain < 10 HZ
3- Phase domain
4- Spectral IP 10-3 to 4000 HZ

55
 Sources of IP Effects

1) Normal IP
 Membrane Polarization
 Most Pronounced with clays
 Decreases with very high (> 10%) clay content due to few pores, low
conductivity.

2) Electrode polarization
 Most metallic minerals have EP
 Decreases with increased porosity.
 Over-voltage effect

3) IP is A bulk effect.

Grain (electrode) polarization. (A) Unrestricted electrolytic flow in an open channel.

(B) Polarization of an electronically conductive grain, blocking a channel

56
1. Time – domain measurements.

One measure of the IP effects is the ratio Vp / Vo which is known chargeability

which expressed in terms of millivolts per volt or percent.

Vp : overvoltage
Vo : observed voltage

M= Vp / Vo ( mv /v or %)

Apparent chargeability
t2
Ma = ( 1 / V0 ) ∫t1 Vp (t) dt = A / V0

Vp ( t) is the over-voltage at time t .

10 – 20 % sulphides 1000-3000 msec .

Sand stones 100-200 msec.
Shale 50-100
Water 0

2) Frequency- Domain measurements.

Frequency effect FE= (Pao –Pa1) / Pa1 ( unitless )

Pao : apparent resistivity at low frequencies

Pal : appatent resistivity at high frequencies

Pao > Pa1

Percentage frequency affect PFE = 100(Pao –Pa1) / Pa1 = 100 FE

The frequency effect in the frequency domain is equivalent to the

chargeability in the time domain for a weakly polarisable medium where FE
<1.

Metal Factor MF= A (ρa0 – ρa1) / (ρa0 ρa1)

= A ( δa1 – δa0 ) siemens / m

57
 ρa0 & ρa1 apparent resistivity.

δa0 and δa1 are apparent conductivities (1/ ρa ) at low and higher frequencies
respectively where

ρa0 > ρa1 and δa0 < δa1 A = 2 π x 105

MF = A x FE / ρa0 = A x FE / ρa0

= FE / ρa0 = A x FE x δa0

The above methods do not give a good indication of the relative amount of the
metallic mineralization within the source of the IP. It is necessary to go with
spectral IP.

3. Spectral IP and Complex Resistivity.

Is the measurement of the dielectric properties of materials

Ө is the phase lag between the applied current and the polarization voltage
measured.

| z(w) | = P0 [ 1 – M ( 1 – 1/ ( 1+(iwτ)c ) ]

Z(w) : complex resistivity

P0 : D.c. resistivity

M : IP chargeability

58
W : Angular frequency.

τ : Time constant. (relaxation time) is the behaviour between the lower

and upper frequency limits.

i : √ -1

c : frequency exponent

Critical Frequency (Fc) : Which is the specific frequency at which the

maximum phase shift is measured. This frequency is completely independent of
resistivity.

Phase angle and the critical frequency increase with increasing chargeability.

Fc = [ 2 π τ ( 1 – M)1/2c ]-1
τ Time constant

M IP chargeability .

This is call cole – cole relalaxation

59
 IP Survey Design
1- Profiling : Later contrasts in electrical properties such as lithologic contacts.
(wenner + Dipole – Dipole) .

2- Sounding : to map the depths and thickness of stratigraphic units

(Schlumberger + wenner).

3- Profiling – Sounding : in contaminant plume mapping , where subsurfae

electrical propertios are expected to vary vertically and horizontally (wenner
+ Dipole – Dipole) .

Limitations of IP
1- IP is more susceptible to sources of cultural interference (metal fences, pipe
lines , power lines) than electrical resistivity.

2- IP equipment requires more power than resistivity alone . This translates

into heavier field instruments

3- The cost of IP much greater than resistivity – alone system.

4- IP requires experience.

5- Complexity in data interpretation.

6- Intensive field work requires more than 3 crew members.

7- IP requires a fairly large area far removed from power lines , fences,
pipelines .

60
 Advantages of IP

1- IP data can be collected during an electrical resistivity survey

2- IP data and resistivity together improves the resolution of the analysis of

Resistivity data in three ways:

a. some of the ambiguities in resistivity data can be redueed by IP

analysis.

b. IP can be used to distinguish geologic layers which do not respond

well to an electrical resistivity .

c. Measurements of chargeability can be used to discriminate equally

electrically conductive target such as saline, electrolytic or metallic-
ion contaminant plumes from clay Layers.

61
 ELECTROMAGNETIC METHODS

Introduction

Electromagnetic methods in geophysics are distinguished by:

1. Use of differing frequencies as a means of probing the Earth (and other planets),
more so than source-receiver separation. Think “skin depth”. Sometimes the
techniques are carried out in the frequency domain, using the spectrum of natural
frequencies or, with a controlled source, several fixed frequencies (FDEM method
---“frequency domain electromagnetic”). Sometimes the wonders of Fourier theory
are involved and a single transient signal (such as a step function) containing, of
course, many frequencies, is employed (TDEM method - “time domain
electromagnetic”). The latter technique has become very popular.

2. Operate in a low frequency range, where conduction currents predominate over

displacement currents. The opposite is true (i.e., has to be true for the method to
work) in Ground Penetrating Radar (GPR). GPR is a wave propagation
phenomenon most easily understood in terms of geometrical optics. Low frequency
EM solves the diffusion equation.

3. Relies on both controlled sources (transmitter as part of instrumentation) and

uncontrolled sources. Mostly the latter is supplied by nature, but also can be
supplied by the Department of Defense.

 EM does not require direct Contact with the ground. So, the speed with EM can
be made is much greater than electrical methods.

 EM can be used from aircraft and ships as well as down

boreholes.

62
Advantages
 lightweight & easily portable.
 Measurement can be collected rapidly with a minimum number of
field personnel .
 Accurate
 Good for groundwater pollution investigations.

Limitations :
 Cultural Noise

Applications
1. Mineral Exploration
2. Mineral Resource Evaluation
3. Ground water Surveys
4. Mapping Contaminant Plumes
5. Geothermal Resource investigation
6. Contaminated Land Mapping
7. Landfill surveys
8. Detection of Natural and Artificial Cavities
9. Location of geological faults
10.Geological Mapping

63
 Type of EM Systems

- EM can be classified as either :

1. Time – Domain (TEM) or
2. Frequency- Domain (FEM)

- FEM use either one or more frequencies.

- TEM makes measurements as a function of time .
- EM can be either :
a- Passive, utilizing natural ground signals
(magnetotellurics)
b- Active , where an artificial transmitter is used
either in the near-field (As in ground conductivity meters) or in the far-field
(using remote high-powered military transmitters as in the case of VLF
Mapping 15-24 KHZ ).

Factors Affecting EM Signal

The signal at the Receiver depends on :
1) the material 2) Shape 3) Depth of the Target
4) Design and positions of the transmitter and receiver coils .

The size of the current induced in the target by the transmitter depends on
1) Number of lines of magnetic field through the Loop (magnetic flux ) 2) Rate of
change of this number 3) The material of the loop.

Magnetic flux Depends on :

1) The Strength of the magnetic field at the Loop 2) Area of the Loop 3)
Angle of the loop to the field
Flux Ø = Magnetic field X cos Ө X area X number of turns .

64
 Principle of EM surveying

- EM field can be generated by passing an alternating

current through either a small coil comprising many turns of wire or a large loop
of wire .
- The frequency range of EM radiation is very wide, from
< 15 HZ ( atmospheric micropulsations) , Through radar bands (108 – 1011 HZ)
up to X-ray and gamma >1016 HZ .
- For geophysical Applications less than few thousand
hertz, the wavelength of order 15-100 km , typical source- receiver separation is
much smaller ( 4-10 m )

The primary EM field travels from the transmitter coil to the receiver coil via paths
both above and below the surface.
In the presence of conducting body, the magnetic component of the EM field
penetrating the ground induces alternating currents or eddy currents to flow in the
conductor.
The eddy currents generate their own secondary EM field which travels to the
receiver. Differences between TX and RX fields reveal the presence of the
conductor and provide information on its geometry and electrical properties.

65
 Depth of Penetration of EM

Skin Depth : is the depth at which the amplitude of a plane wave has decreased to
1/e or 37% relative to its initial amplitude Ao .
Amplitude decreasing with depth due to absorption at two frequencies

Az = Ao e-1

The skin depth S in meters = √ 2 / ωσ μ = 503 √ f σ

ω = 2π f = 503 √ ρ / f = 503 √ρ λ / v
σ : conductivity in s/m
μ : magnetic permeability (usually ≈ 1)
λ : wavelength , f : frequency , v : velocity , p : Resistivity thus, the
depth increases as both frequency of EM field and conductivity decrease.
Ex. In dry glacial clays with conductivity 5x 10-4 sm-1 , S is about 225 m at a
frequency of 10 KHZ .

Skin depths are shallower for both higher frequencies and higher conductivities
(Lower resistivities ).

66
 Magnetotelluric Methods ( MT )

Telluric methods: Faraday's Law of Induction: changing magnetic fields

produce alternating currents. Changes in the Earth's magnetic field produce
alternating electric currents just below the Earth's surface called Telluric currents.
The lower the frequency of the current, the greater the depth of penetration.
Telluric methods use these natural currents to detect resistivity differences which
are then interpreted using procedures similar to resistivity methods.

MT uses measurements of both electric and magnetic components of The Natural

Time-Variant Fields generated.
Major advantages of MT is its unique Capability for exploration to very great
depths (hundreds of kilometers) as well as in shallow Investigations without using
of an artificial power source .
Natural – Source MT uses the frequency range 10-3-10 HZ , while audio –
frequency MT (AMT or AFMAG) operates within 10-104 HZ .

The main Application of MT in hydrocarbon Expl. and recently in meteoric impact,

Environmental and geotechnical Applications.

Pa = 0.2 / f │Ex / By │2 = 0.2 / f │Ex / Hy │2 = 0.2 / f │Z│2

Ex (nv/km) , By , orthogonal electric and magnetic components.

By : magnetic flux density in nT .
Hy : magnetizing force (A/m) .
Z : cagniard impedance.

The changing magnetic fields of the Earth and the telluric currents they produce
have different amplitudes. The ratio of the amplitudes can be used to determine the
apparent resistivity to the greatest depth in the Earth to which energy of that
frequency penetrates.

Typical equation:

apparent resistivity =

67
where Ex is the strength of the electric field in the x direction in millivolts
Hy is the strength of the magnetic field in the y direction in gammas
f is the frequency of the currents

Depth of penetration =

This methods is commonly used in determining the thickness of sedimentary

basins. Depths are in kilometers

68
 Field Procedure

MT Comprises two orthogonal electric dipoles to measure the two horizontal

electric components and two magnetic sensors parallel to the electric dipoles to
measure the corresponding magnetic components .

1. Two orthogonal grounded dipoles to measure electric components

2. Three orthogonal magnetic sensors to measure magnetic components.

Thus, at each location, five parameters are measured simultaneously as a function

of frequency. By measuring the changes in magnetic (H) and electric (E) fields over
a range of frequencies an apparent resistivity curve can be produced. The lower the
frequency, the grater is the depth penetration.

69
 Survey Design

EM data can be acquired in two configurations

1) Rectangular grid pattern
2) Along a traverse or profile .

EM equipment Operates in frequency domain. It allows measurement of both the .

1) in-phase (or real ) component .
2) 900 out – of – phase (or quadrature ) component.

Very Low Frequency (VLF) Method

VLF : uses navigation signal as Transmitter .

Measures tilt & phase

Main field is horizontal .

VLF detects electrical conductors by utilizing radio signal in the 15 to 30 KHZ

range that are used for military communications.
VLF is useful for detecting long, straight electrical conductors
VLF compares the magnetic field of the primary signal (Transmitted ) to that of
the secondary signal ( induced current flow with in the subsurface electrical
conductor).

Advantages of VLF
1) Very effective for locating zones of high electrical conductivity
2) fast
3) inexpensive
4) Requires one or two people .

70
 Tilt Angle Method
Tilt angle systems have no reference link between Tx and Rx coils . Rx measures
the total field irrespective of phase and the receiver coil tilted to direction of
maximum or minimum magnetic field strength .

The response parameter of a conductor is defined are the product of conductivity –

thickness ( T) , permeability (μ ) an angular frequency
ω = 2π f and the square of the target a2 .

Poor conductors have response parameter < 1

Excellent conductor have response parameter greater than 1000

A Good conductor having a higher ratio AR / Ai

AR : Amplitude of Read (in – phase )
Ai : Amplitude of imaginary ( out – of phase)

In the left side of the above figure and poor conductor having a lower ratio of AR/
Ai .

71
 Slingram System
 slingram is limited in the size of TX coil. This system has the Transmitter and
Receiver connected by a cable and their separation kept constant as they are
moved together along a traverse.

Magnetic field Through The receiver has two sources :

a) The primary field of The Transmitter .
b) The secondary field produced by The Target .

Turam system
More powerful system than Slingram. It uses a very large stationary Transmitter
coil or wire laid out on the ground, and only The receiver is moved . TX 1-2 km
long, loop over 10 km long. The receiver consists of two coils and kept a fixed
distance between 10-50 m apart.

72
 Ground Surveys of EM

A. Amplitude measurement
1- Long wire

 Receiver pick up horizontal component of field parallel to wire .

 Distortions of Normal field pattern are related to changes in
subsurface conductivity.

B. Dip-Angle
Measures combined effect of primary and secondary fields at the receiver.

AFMAG : Dip-angle method that uses Naturally occurring ELF signals

generated by Thunder storms.

Phase Component Methods

1) Work by comparing secondary & primary fields .

2) Compensator & Turam (long wire) .
3) Slingram Moving Trans / receiver
- Penetration ≈ ½ spacing of coil .
- Coil spacing critical .

 Over barren ground Null is at zero coil dip-Angle.

 Near conductor, dip angle ≠ 0
 Dip – Angle is zero over Narrow conductor, and changes sign.

73
 Dip-Angle method

1) easy , cheap
2) Quick
3) Sensitive to vertical
4) Difficult to distinguish between depth & conductivity

TDEM Method

A significant problem with many EM surveying techniques is that a small

secondary field must be measured in the presence of a much larger primary field,
with a consequent decrease in accuracy. This is circumvented to some extent in the
FDEM method described about by measuring the out-of-phase component.
In a TDEM approach, the signal is not a continuous frequency but instead consists
by a series of step-like pulses separated by periods where there is no signal
generated, but the decay of the secondary field from the ground is measured. The
induction
currents induced in a subsurface conductor diffuse outward when the inducing
energy is suddenly switched off. The measurement of the field at a number of time
steps is equivalent of measuring at several frequencies in an FDEM system.
Usually two-coil systems are used and the results can be stacked to reduce noise.
Modeling of the decay for layered systems, and more complicated conductivity
geometry, can be carried out.

74
The figure above shows the behavior of the field for two different conductivities. If
samples at different times were taken, then on could distinguish between
the two conductivities. This is the principle of Time Domain
Electromagnetic (TDEM) methods.

τ = σ μ L2 = σ μ A
where τ is a characteristic time constant and L and A correspond to a characteristic
length scale and characteristic area, respectively.

The EM61 has a single time sample at t  0.5 ms. Using a cylinder of radius 2 cm
and a conductivity of steel of 107 S m-1, then exp(-t/ τ ) = 0.56. where t is time
constant.
On the other hand, assume a plastic drum of seawater of conductivity 10 S m-1 and
radius 40 cm, then we obtain exp(-t/ τ ) = 0.

75
 Airborne Electromagnetic Surveys
The general objective of AEM (Airborne ElectroMagnetic) surveys is to conduct a
rapid and relatively low-cost search for metallic conductors, e.g. massive sulphides,
located in bed-rock and often under a cover of overburden and/or fresh water. This
method can be applied in most geological environments except where the country
rock is highly conductive or where overburden is both thick and conductive. It is
equally well suited and applied to general geologic mapping, as well as to a variety
of engineering problems (e.g., fresh water exploration.) Semi-arid areas,
particularly with internal drainage, are usually poor AEM environments.

Conductivities of geological materials range over seven orders of magnitude, with

the strongest EM responses coming from massive sulphides, followed in decreasing
order of intensity by graphite, unconsolidated sediments (clay, tills, and
gravel/sand), and igneous and metamorphic rocks. Consolidated sedimentary rocks
can range in conductivity from the level of graphite (e.g. shales) down to less than
the most resistive igneous materials (e.g. dolomites and limestones). Fresh water is
highly resistive. However, when contaminated by decay material, such lake bottom
sediments, swamps, etc., it may display conductivity roughly equivalent to clay and
salt water to graphite and sulphides.

Typically, graphite, pyrite and or pyrrhotite are responsible for the observed
bedrock AEM responses. The following examples suggest possible target types and
we have indicate the grade of the AEM response that can be expected from these
targets.

 Massive volcano-sedimentary stratabound sulphide ores of Cu, Pb, Zn, (and

precious metals), usually with pyrite and/or pyrrhotite. Fair to good AEM
targets accounting for the majority of AEM surveys.
 Carbonate-hosted Pb-Zn, often with marcasite, pyrite, or pyrrhotite, and
sometimes associated with graphitic horizons. Fair to poor AEM targets.
 Massive pyrrhotite-pentlandite bodies containing Ni and sometimes Cu and
precious metals associated with noritic or other mafic/ultramafic intrusive
rocks. Fair to good AEM targets.
 Vein deposits of Ag, often with Sb, Cu, Co, Ni, and pyrite in volcanic and
sedimentary rocks. Generally poor AEM targets.
 Quartz veins containing Au with pyrite, sometimes also with Sb, Ag, Bi, etc.,
in volcanic or sedimentary (and possibly intrusive) rocks. Poor AEM targets.

76
 Basic Principles of Airborne
Electromagnetic-induction prospecting methods, both airborne and (most) ground
techniques, make use of man-made primary electromagnetic fields in, roughly, the
following way: An alternating magnetic field is established by passing a current
through a coil, (or along a long wire). The field is measured with a receiver
consisting of a sensitive electronic amplifier and meter or potentiometer bridge.
The frequency of the alternating current is chosen such that an insignificant eddy-
current field is induced in the ground if it has an average electrical conductivity,

If the source and receiver are brought near a more conductive zone, stronger eddy
currents may be caused to circulate within it and an appreciable secondary
magnetic field will thereby be created. Close to the conductor, this secondary or
anomalous field may be compared in magnitude to the primary or normal field
(which prevails in the absence of conductors), in which case it can be detected by
the receiver. The secondary field strength, Hs, is usually measured as a proportion
of the primary field strength, Hp, at the receiver in percent or ppm (parts per
million).

Anomaly = Hs / Hp.

Increasing the primary field strength increases the secondary field strength
proportionally but the "anomaly" measured in ppm or percent remains the same.
prospecting for anomalous zones is carried out by systematically traversing the
ground either with the receiver alone or with the source and receiver in
combination, depending on the system in use. In the case of airborne systems, the
receiver coils are usually in a towed bird and the transmitter may be a large coil
encircling a fixed wing aircraft, e.g. INPUT systems, or one or more small coils in
the same bird that houses the transmitting coils, e.g. most HEM (Helicopter EM)
systems.
There are two different basic systems commonly used to generate and receive the
electromagnetic field: transient or "time domain" systems like INPUT,GEOTEM
and MEGATEM and a/c. "frequency domain" systems like most HEM systems.

77
 Transient Airborne Electromagnetics
Historically, the most commonly encountered system of this type was the INPUT
system. The newer systems GEOTEM and MEGATEM (Fugro Airborne Surveys)
function in a similar way to INPUT

In the INPUT system the transmitting coil, usually encircling a fixed wing aircraft,
is energized by what is, essentially, a step current. In the absence of conductors, a
sharp transient pulse proportional to the time derivative of the magnetic field is
induced in the receiver. When a conductor is present, however, a sudden change in
magnetic field intensity will induce in it a flow of current in the conductor which
will tend to slow the decay of the field.

The receiver "listens" only while the transmitter is "quiet" so that problems arising
out of relative motion between transmitter and receiver, because the receiver is
towed in a bird behind the aircraft, are virtually eliminated. Moreover, if the entire
decay of the secondary field could be observed, the response would be equivalent
to AC measurements made over the whole of the frequency spectrum. It is
important to note in this connection, however, that not the decay function itself but
only its time derivative can be recorded if a coil is used as the detector. This means
that the anomalous fields which decay very slowly are suppressed in amplitude
more than the others, and since these are the very ones generally associated with
good conductors, there would seem to be an inherent weakness in this system.
Because it is difficult to precisely synchronize the instant when the transmitter
becomes "quiet" with the instant that the receiver begins to "listen", it is nearly
impossible to record the entire function. This is equivalent to being unable to
record many of the lower frequencies in the a-c spectrum.

Typically, the time derivative of the decay function is measured using from six to
twelve different time delays from the instant that transmitter stops transmitting
before recording the signal received.

78
A sketch of the INPUT transient airborne EM system operation. The primary field
is a step function and the receiver records the decay of the field after the transmitter
stops transmitting. (Grant and West 1965)

79
 Frequency Domain Airborne Electromagnetics
In the typical frequency domain helicopter EM system (HEM) both the transmitting coil set and
the receiver coil set are housed in a rigid boom or "bird" that is towed beneath the helicopter.
Commonly, this boom is from three to five meters long and contains from two to six coil pairs.
Usually, half of the coils in each of the transmitter set and the receiver set are "co-axial", i.e. an
axis normal to the plane of the coils passes through the centre of both coils. The second half of
the coil sets are normally "co-planar", being equivalent to both the transmitting and receiving coil
lying flat on the ground. Other coplanar orientations have been used occasionally.

The receiver measures the in-phase and out-of-phase, or quadrature, of the

secondary field, expressed in ppm of the primary field. The two different coil
orientations provide data that is useful in discriminating between dike like
conductors that have considerable vertical extent and may be ore bodies, and
horizontal sheet like conductors that are simply conductive overburden.

Factors Affecting Detectability

80
1. Signal-to-noise ratio:

In practice, because of "system noise" (Ns) and "geological noise" (Ng), the ability
of a system to recognize and measure an anomaly is limited by the "signal-to-
noise" ratio: Signal-to-noise = Hs / (Ns + Ng)
Because Hs and Ng are proportional to the primary field strength Hp, and Ns, in
frequency-domain systems, usually contains elements proportional to Hp, there is
little to be gained by increasing the primary field power. In time domain systems
Ns is not greatly affected by Hp, so extra power does result in increased signal-to-
noise. Attempts to increase the signal-to-noise are sometimes made by increasing
the distance between the transmitter and receiver. This results in roughly the same
Hs and Ng but often a lower system noise Ns.

2. Penetration

The penetration of an AEM system is its effective depth of exploration. Commonly,

this is taken to include the elevation of the system above ground, as this is also
affected by local environment and flying conditions.
In general, systems with large transmitter-receiver coil separation, usually referred
to as Tx-Rx, have greater penetration than those with small separations. Penetration
is closely related to signal-to-noise, as the system that produces a larger anomaly
from a given conductor can, of course, look further into the ground.

3. Discrimination

The discrimination of an AEM system is the ability of the system to differentiate

between conductors of different physical properties or geometric shapes.
Discrimination, particularly between flat lying surficial conductors and steeply
dipping conductors, is vitally important. Good discrimination can be achieved in
HEM systems by using several frequencies and both co-axial and co-planar coil
pairs.

4. Resolution

Resolution refers to the ability of an AEM system to recognize and separate the
interfering effects of nearby conductors. A system that does this well also produces
sharp anomalies over isolated or discrete conductors. Resolution generally
increases with decreasing flight elevation and coil separation. Typically the HEM
systems have better resolution than the fixed wing time domain systems.

Typical Electrical Properties of Earth Materials.

81
 Conductivity Resistivity (ohm-
Rock, Mineral, etc.
(mohs/meter) meters)
Bornite 330 3 x 10-3
Chalcocite 104 10-4
Chalcopyrite 250 4 x 10-3
Galena 500 2 x 10-3
Graphite 103 10-3
Marcasite 20 5 x 10-2
Magnetite 17 x 10-4 - 2 x 104 5 x 10-5 - 6 x 10-3
Pyrite 3 0.3
Phrrhotite 104 10-4
Sphalerite 10-2 102
Igneous and
10-7 - 10-2 100 - 107
Metamorphic Rocks
Sediments 10-5 - 5 x 10-2 20 - 105
Soils 10-3 - 0.5 2 - 103
Fresh Water 5 x 10-3 - 0.1 10 - 200
Saline Overburden 0.1 - 5 0.2 - 1
Salt Water 5 - 20 0.05 - 2
Sulphide Ores 10-2 - 10 0.1 - 100
Granite Beds and Slates 10-2 - 1 1 - 100
Altered Ultramafics 10-3 - 0.8 1.25 - 103
Water-filled
10-3 - 1 1 - 103
faults/shears

82
 Ground penetrating radar (GPR )
Ground penetrating radar is a nondestructive geophysical method that produces a
continuous cross-sectional profile or record of subsurface features, without drilling,
probing, or digging. Ground penetrating radar (GPR) profiles are used for
evaluating the location and depth of buried objects and to investigate the presence
and continuity of natural subsurface conditions and features.

Ground penetrating radar operates by transmitting pulses of ultra high frequency

radio waves (microwave electromagnetic energy) down into the ground through a
transducer or antenna. The transmitted energy is reflected from various buried
objects or distinct contacts between different earth materials. The antenna then
receives the reflected waves and stores them in the digital control unit.

The ground penetrating radar antenna (transducer) is pulled along the ground by
hand or behind a vehicle.

When the transmitted signal enters the ground, it contacts objects or subsurface
strata with different electrical conductivities and dielectric constants. Part of the
ground penetrating radar waves reflect off of the object or interface; while the rest
of the waves pass through to the next interface.

The reflected signals return to the antenna, pass through the antenna, and are
received by the digital control unit. The control unit registers the reflections against
two-way travel time in nanoseconds and then amplifies the signals. The output
signal voltage peaks are plotted on the ground penetrating radar profile as different
color bands by the digital control unit.

For each reflected wave, the radar signal changes polarity twice. These polarity
changes produce three bands on the radar profile for each interface contacted by the
radar wave.

Ground penetrating radar waves can reach depths up to 100 feet (30 meters) in low
conductivity materials such as dry sand or granite. Clays, shale, and other high
conductivity materials, may attenuate or absorb GPR signals, greatly decreasing the
depth of penetration to 3 feet (1 meter) or less.

The depth of penetration is also determined by the GPR antenna used. Antennas
with low frequencies of from 25 to 200 MHz obtain subsurface reflections from
deeper depths (about 30 to 100 feet or more), but have low resolution. These low

83
frequency antennas are used for investigating the geology of a site, such as for
locating sinkholes or fractures, and to locate large, deep buried objects.

Antennas with higher frequencies of from 300 to 1,000 MHz obtain reflections
from shallow depths (0 to about 30 feet), and have high resolution. These high
frequency antennas are used to investigate surface soils and to locate small or large,
shallow buried objects and rebar in concrete.

GPR is a high-resolution technique of imaging shallow soil and ground structures

using electro-magnetic (EM) waves in the frequency band of 10-1000 MHz.
Advantage of EM waves is that signals of relatively short wavelength can be
generated and radiated into the ground to detect anomalous variations in dielectric
properties of the geological material

GPR – technical summary

  f ~ 50 MHz - 5 GHz
  v ~ 0.05 - 0.15 m/ns
  λ ~ 1 m - 1cm
  Reach : 0 - 50 λ

Used on surface with one sensor for transmitting and receiving or double system of
separate transmitter/receiver Used in boreholes (3D borehole radar, or cross
imaging)

Principles of operation
A radar system comprises a:
signal generator
transmitting / receiving antennae
recording unit

84
85
‫الشكل يوضح النظام الرادأاري حيث يحتوي على هوائي للرسال ‪TX‬وآخآر للستقبال ‪RX‬‬

‫‪86‬‬
‫الشكل يوضح المعاملت التي تسبب فقد للطاقأة وهي )‪ (1‬النعكاس وفقد الطاقأة بين الحدودأ‬
‫الطبقية‪ (2) ،‬تشتت للطاقأة بطريقة عشوائية‪ ،‬ويعرف ذلك بالتشتت أو التبعثر الكمي ‪(3) ،‬‬
‫المتصاص‪ (4) .‬النتشار الجيومتري للطاقأة‪.‬‬

‫‪87‬‬
GPR - application
  Geological
  Environmental
  Glaciological
  Engineering and construction
  Archaeology
  Forensic science

 Geological
  Detection natural cavities and fissures
  Subsidence mapping
  Mapping of superficial deposits
  Soil stratigraphy mapping
  Geological structure mapping
  Mapping of faults, dykes, coal seams
  Lake and riverbed sediment mapping
  Mineral exploration and resource evaluation

Engineering and construction

  Road pavement analysis
  Void detection
  Location of reinforcement (rebars) in concrete
  Location of public utilities (pipes, cables, etc)
  Testing integrity of building materials
  Concrete testing

88
 

What does GPR measure?

GPR systems essentially measure the signal travel time, the time between sending
the pulse at the transmitter antenna and the moment the (distorted) pulse is received
back at the receiver antenna.
In general one measures not only one peak but a series of peaks related to various
objects in the sub-surface. The inhomogeneity of the sub-surface is the main factor
that controls the number of peaks that will be returned.

Resolution
 In practice, it is often better to accept lower spatial resolution in favour of range
where there are many thin layers or scattering targets that are not of primary
interest.

 Low frequency GPR

1.  Deep penetration
2.  Low resolution

 High frequency GPR

1.  Shallow penetration
2.  Very high resolution

89
Data processing
Data processing for GPR is very similar to that used in seismics. Most seismic
software packages can deal with radar measurements.
Time – depth conversions can be made through estimation of velocity. There
are various ways to do this. CMP measurements is just one possibility.

90
91
92
93
 SOLVED PROBLEMS

1. Determine the potential difference between the two potential electrodes for
cases (a) and (b). Assume a current of 0.6 ampere.

(a)
i   i  i  1 1  i   1 1 
VP1       
2  r1  2  r2  2   r1 r2  2  12.17  m 21.63  m 

VP1  0.686  v
i  1 1 
VP2  
2  12.65  m 20  m 

VP2  0.555  v
VP1  VP2  0.131  v

(b)
VP1  0.418  v
VP2  0.418  v
VP1  VP2  0.836  v

94
2. Construct the current-flow lines beneath the interface in (a) and (b).

(a)

(b)
 2  66.6Þ

95
3. Calculate the potential at ρ1, due to a current at C, of 0.6 ampere. The material in this
section view extends to infinity in all directions. The bold line represents an interface
between p1- and p2-material.

3.

i1 ik 1 2   1
VP1   , k  0.6
4 r1 4 r2 2   1
0.6amp 200 ohm m  0.6amp 0.6 200 ohm m 
VP1  
4 10 m  4 16.67m 
VP1  0.611  v

96
 4. The following data were acquired using a constant-spread, Wenner traverse (a-
spacing = 3m). Interpret the data as completely as possible.

Horizontal Horizontal Horizontal

Position pa Position (m) pa Position (m) Pa
5 20.05 19.00 22.55 33.00 858.94
6 20.06 20.00 25.32 34.00 872.33
7.00 20.07 21.00 28.45 35.00 880.40
8.00 20.08 22.00 27.97 36.00 885.58
9.00 20.10 23.00 27.11 37.00 889.07
10.00 20.12 24.00 242.39 38.00 891.51
11.00 20.15 25.00 460.00 39.00 893.26
12.00 20.19 26.00 572.39 40.00 894.57
13.00 20.24 27.00 580.25 41.00 895.55
14.00 20.32 28.00 541.30 42.00 896.31
15.00 20.44 29.00 519.88 43.00 896.90
16.00 20.61 30.00 660.54 44.00 897.38
17.00 20.91 31.00 785.22 45.00 897.76
18.00 21.45 32.00 834.74

This data suggests a vertical discontinuity at 25 m horizontal position. Resistivity

of the material to the left of the contact is 20 ohm·m and that of the material to the
right is 900 ohm·m.

97
5. The following data were gathered with a Wenner, expanding-spread traverse in an
area with thick alluvial deposits at the surface. What is the likely depth to the water
table?
Electrode Spacing pa Electrode Spacing (m) pa
0.47 198 6..81 84
0.69 160 10..00 82
1..00 140 14.68 92
1.47 112 21.54 101
2.15 95 31.62 100
3.16 84 46.42 102
4.64 79

A consistent model shows that layer 1 is 1.3 m, 207 ohm·m; layer 2 is 15.7 m, 77
ohm·m; layer 3 is 107 ohm·m. A variety of input models produce similar results.
The water table, therefore, is judged to be at a depth of approximately 1.3 m.

98
 HOMEWORK ASSIGNMENTS

1. CONSIDER A VERTICAL CONTACT BETWEEN TWO GEOLOGIC UNITS OF GREATLY DIFFERENT

RESISTIVITIES P1 AND P2 .

A. CALCULA TE THE DIRECTION THAT ELECTRICAL CURRENT WILL LEAVE THE BOUNDARY INTO
P2 IF IT APPROACHES THE BOUNDARY AT AN ANGLE OF 45° TO THE NORMAL.

B. CALCULATE THE ELECTRICAL POTENTIAL AT POINT A FROM A SINGLE SOURGE OF CURRENT (

+I ) A T POINT B. REMEMBER THAT THE BOUNDARY WILL CREATE A REFLECTED IMAGE OF
THE CURRENT SOURCE WITH MAGNITUDE KI. FOR P1 = 100 Ωm AND P2 = 0.1 Ωm

C. REPEAT (A) AND (B) FOR P1 = 100 Ωm AND P2 = 10000 Ωm

99
2. Determine the potential difference between the two potential electrodes for cases (a) and (b). Assume a
current of 0.5 ampere for (a) and 0.8 ampere for (b).

3. Construct the current-flow lines beneath the interface in (a) and (b).

100
4 - Interpret the following data, which were obtained with a Schlumberger traverse.

Electrode  ( Ω . m) Electrode  ( Ω . m) Electrode  ( Ω . m)

Spacing (m) Spacing (m) Spacing (m)
I.00 108 14.68 307 215.44 293
1.47 121 21.54 245 316.23 381
2.15 I4X 31.62 168 464.16 479
3.16 191 46.42 122 681.29 580
4.64 244 68.13 115 1000.0 675
6.81 295 100 162
10.00 323 146.78 220

5- The following data were gathered with a Wenner, expanding-spread traverse in an area of thick
deltaic sands. Bedrock depths are greater than 30 m. What is your best estimate of the depth to the
water table in this area?

Electrode Spacing (m)  ( Ω . m) Electrode Spacing (m)  ( Ω . m)

0.47 2590 6.81 8753
0.69 3288 10.00 7630
1.00 4421 14.68 4805
1.47 5198 21.54 2160
2.15 6055 31.62 995
3.16 6686 46.42 584
4.64 7782

6. The following data were gathered with a Wenner. expanding-spread traverse in an area of dune sands
underlain by lake clays which in turn are underlain by Triassic sedimentary rocks. Estimate a
value for clay thickness.

Electrode Spacing (m)  ( Ω . m) Electrode Spacing (m)  ( Ω . m)

0.69 1298.90 10.00 360.67
1.00 1398.06 14.68 240. 1 1
1.47 1306.98 21.54 191.51
2.15 1153.02 31.62 153.81
3.16 925.27 46.42 116.38
4.64 762.40 68.13 98.03
6.81 554.13 100.00 86.08

101
7. Interpret the following data, which were obtained with a Wenner traverse.

Electrode pa (Ω • m) Electrode Spacing pa (Ω • m) Electrode pa (Ω • m)

Spacing (m) (m) Spacing (m)
1 .00 984 14.68 72 215.44 327
1.47 955 21.54 53 316.23 432
2.15 883 3 1 .62 63 464.16 554
3.16 742 46.42 87 68 1 .29 686
4.64 533 68.13 124 1000.00 818
6.81 311 100.00 174
10.00 150 146.78 241

8. Assume a homogeneous medium of resistivity 120 ohm-m. Using the wenner electrode system
with a 60-m spacing, assume a current of 0.628 ampere. What is the measured potential differ-
ence? What will be the potential difference if we place the sink (negative-current electrode) at
inf i ni t y ?

9. Suppose that the potential difference is measured with an electrode system for which one of the
current electrodes and one of the potential electrodes are at infinity. Using the Figure below
and a current of 0.5 ampere, compute the potential difference between the electrodes at PA and
infinity for d1 = 50 m, d2 = 100 m, R1 = 30 ohm-m, R2 = 350 ohm-m.

102
10. Plot resistivity data as a function of electrode spacing, and determine the particular electrode
spacing corresponding to the inflection point for the data given in Exercise 9. Compare the
electrode spacing at the inflection point with the depth of the boundary from Exercise 9.

11. If the Schlumberger electrode system with AB/2 : MN/2 = 5. is used to conduct the resistivity
survey explained in Exercise 9, what will the potential readings be? Use resistivity values found
in Exercise 9 to compute potential differences for each of the electrode spacings AB/2 = 1, 2, 4,
6, 8, 10, 15, 20, 25, 30, 40, 50 m w h i l e a constant cu rrent of 0.250 ampere is applied.

12. Suppose that an electrical resistivity survey was done using an expanding Wenner electrode-
configuration. The current of 0.25 ampere was kept the same for all the readings. Potential d i f -
ferences measured wit h different electrode spacings are given in the following table.
Interpretations of these measurements indicate that a layer of resistivity R1 lies above another
layer of resistivity R.2. Determine the depth of the boundary between these two layers. Estimate
the resistivities of the layers.

ELECTRODE SPACING (m) POTENTIAL DIFFERENCE (v)

1 0.8

2 0.42

4 0.28

6 0.155

8 0.125

10 0.120

15 0.105

20 0.10

25 0.098

30 0.086

40 0.076

50 0.064

13. For purposes of an IP survey, resistivity values are determined from both direct and alternating
current using the same electrode arrangement. If the resistivities for direct and alternating cur-
rent are Rdc = 50 ohm-m and Rac = 40 ohm-m, respectively, what w i l l th e frequency effect and
the metal factor values be.

103
14. Assume- that a telluric current survey is to be carried out to outline large-scale features
of a sedimentary basin 5 km deep. A resistivity of 50 ohm-m is supposed to represent
the sedimentary section. What is the maximum frequency of the telluric current that
will penetrate below the basin?

15. Suppose That a magnetotelluric survey indicates an apparent resistivity of 5 ohm-

m at a frequency of 1 Hz, What is the thickness of the layer?

16. Using the method of electrical images, drive therelationship between apparent resistivity,
electrode spacing, layer thicknesses and resistivities for a VES performed with a
Schlumberger spread over a single horizontal interface between media with resistivities ρ1
and p2.

17. Calculate the variation in apparent resistivity along a HEP profile at right angles to a
vertically faulted contact between sandstone and limestone, with apparent resistivities
of 50 ohm m and 600 ohm m. respectively, for a Wenner configuration. What
would be the effect on the profiles if the contact dipped at a shallower angle.

18. Why are the electrical methods of exploration particularly suited to hydrogeological
investigations? Describe other geophysical methods which could be used in this context,
stating the reasons why they are applicable.

19. Let A and B represent two different geologic sections, and let H1 and H2 represent the thicknesses of
the first and second layers in a three-layer sequence. It is well known that one type of equivalence
occurs for a three-layer case when p1 < p2 > p3 . p 1 A = p1B , p3A= p3B , H1A = H 1B , and P2A. H2A = P2B.
H2B . Demonstrate that this is true. Does equivalence exist if p1 > p2 < p3 ?. Explain.

104
20. At locations A, B, C , D along the gravity profile shown below, VES were performed with a
Wenner array with the spread laid perpendicular to the profile.

It was found that the sounding curves, were similar for locations A ,and B and for C
and D. A borehole close to A penetrated 3m of drift, 42 m of limestone and
bottomed in sandstone. Downhole geophysical surveys provided the following values
of density ( D ) and resistivity (ρR) for the lithologies encountered.

Unit R (Ω m) D ( Mg m-3)
Drift 40 2.00
Limestone 2000 2.75
Sandstone 200 2.40

A seismic refraction line near to D revealed 15m of drift, although the nature of the
underlying basement could not be assessed from the seismic velocity.

(a) Interpret the geophysical data so as to provide a geological section along the
profile.
(b) What further techniques might be used to confirm your interpretation?
(c) If a HEP were to be performed along the profile, select, giving reasons, a suitable
electrode spacing to map the basement. Sketch the expected form of the HEP for both
longitudinal and transverse traverses.

105
21. The following table represents the results of a frequency domain IP survey of a
Precambrian shield area. A dipole-dipole array was used with the separation ( x) of
both the current electrodes and the potential electrodes kept constant at 60m. n refers to
the number of separations between the current and potential electrode pairs and c to the
distance of the center of the array from the origin of the profile, where the results are
plotted (Figure below). Measurements were taken using direct current and an
alternating current of 10Hz. These provided the apparent resistivities p dc and pac
respectively,

(a) For each measurement point, calculate the percentage frequency effect (PFE) and
metal factor parameter (MF).

The dipole-dipole electrode configuration.

(b) For both the PFF, and MF plot four profiles for n = 1,2 , 3 and 4.

106
 (c) The area is covered by highly-conductive glacial deposits 30-60m thick, It is possible that
massive sulphide mineralization is present within the bedrock. Bearing this information in
mind, comment upon and interpret the profiles.
22. For each of the following subsurface models sketch an appropriate apparent resistivity curve on the
designated graph. The general shape of the curve is what is important. Base your curves on what you
know about current penetration, current density, and measured resistivities.

(a) Layer Thickness (m) Resistivity ( Ω .m)

1 1 100
2 5 10
3 5 100
4 1
Infinite

(b) Layer Thickness (m) Resistivity ( Ω .m)

1 10 100
2 10 10
3 Infinite 100

107
108
23. The figure below shows a half-Schlumberger resistivity array in which the second
current electrode is situated at a great distance from the other electrodes. Derive an
expression for the apparent resistivity of this array in terms of the electrode spacings and
the measured resistance.

The half-Schlumberger electrode configuration.

24. Calculate the potential at P1 due to a current at C1 of 0.6 ampere. The material in this section view
extends to infinity in all directions. The bold line represents an interface between p1 and p2 –material.

109
 ‫مصطلحات السإتكشاف الكهربائي‬
‫والكهرومغناطيسي‬

‫المصطلح‬ ‫الترجمة العربية‬

Self Potential ‫جهد ذاتي‬
Galvanometer ‫مقياس الجهد الكهربائي‬
Apparent Resistivity ‫مقاومية ظاهرية‬
Equipotential line ‫خآط متساوي الجهد‬
Principle of Equivalence ‫مبدأ التعادأل‬
Principle of Suppression ‫مبدأ الغطس‬
Conductivity ‫التوصيلية‬
Vertical Electrical Sounding ‫الجس الكهربائي العمودأي‬
Horizontal Electrical Profiling ‫المقطع الكهربائي الفقي‬
Dipole ‫ثنائي القطب‬
Schlumberger Arrangement ‫ترتيب شلمبرجير‬
Induced Polarization ‫الستقطابية المستحثة‬
Passive Methods ‫الطرق الخامدة‬
Electrical Conduction ‫التوصيل الكهربائي‬
Archie's Law ‫قأانون أرشي‬
Electrical Reflection Coefficient ‫معامل النعكاسية الكهربائية‬
Ground Penetrating Radar (GPR) ‫رادأار الخآتراق الرضي‬
Airborne Electromagnetic ‫المسح الكهرومغناطيسي الجوي‬
Very Low Frequency (VLF) ‫التردأدأ المنخفض جدا‬

110
‫مصطلحات السإتكشاف الجاذبي والمغناطيسي‬
‫‪Absolute gravity‬‬ ‫جاذبية مطلقة‬
‫‪Magnetic Field Strength‬‬ ‫قأوة المجال المغناطيسي‬
‫‪Magnetic Induction‬‬ ‫الحث المغناطيسي‬
‫‪Centrifugal Force‬‬ ‫قأوه طردأ مركزية‬
‫‪Intensity of Magnetization‬‬ ‫شدة التمغنط‬
‫‪Magnetic Declination‬‬ ‫انحراف مغناطيسي‬
‫‪Diamagnetic‬‬ ‫ضعيف النفاذية المغناطيسية‬
‫‪Diurnal Correction‬‬ ‫تصحيح يومي‬
‫‪Elevation Correction‬‬ ‫تصحيح الرتفاع‬
‫‪Equator‬‬ ‫خآط الستواء‬
‫‪Magnetic Permeability‬‬ ‫النفاذية المغناطيسية‬
‫‪Ferro magnetic‬‬ ‫مغناطيس حديدي‬
‫‪Gravitational Acceleration‬‬ ‫التسارع الجاذبي‬
‫‪Gravity Anomaly‬‬ ‫شاذة الجاذبية‬
‫‪Magnetic Inclination‬‬ ‫الميل المغناطيسي‬
‫‪Isostatic Correction‬‬ ‫تصحيح ايزوستاتي‬
‫‪Lunar Variations‬‬ ‫تغيرات قأمرية‬
‫‪Magnetic Moment‬‬ ‫العزم المغناطيسي‬
‫‪Magnetic Storms‬‬ ‫عواصف مغناطيسية‬
‫‪Magnetometer‬‬ ‫جهاز قأياس المغناطيسية‬
‫‪Magnetic Susceptibility‬‬ ‫قأابلية مغناطيسية )التأثرية المغناطيسية(‬
‫‪Observed Gravity‬‬ ‫جاذبية مقاسة‬
‫‪Paleomagnetism‬‬ ‫مغناطيسية قأديمة‬
‫‪Residual Magnetism‬‬ ‫مغناطيسية متخلفة‬
‫‪Remanent Magnetism‬‬ ‫مغناطيسية متبقية‬
‫‪Secular Variations‬‬ ‫تغيرات متناهية البطء‬
‫‪Geoid‬‬ ‫الجيوئد ) سطح متساوي الجهد (‬
‫‪Gravimeter‬‬ ‫جهاز قأياس الجاذبية‬
‫‪Bouguer Anomaly‬‬ ‫شاذة بوجير‬
‫‪Latitude Correction‬‬ ‫تصحيح خآط العرض‬

‫‪111‬‬
‫مصطلحات السإتكشاف السيزمي‬

‫المصطلح‬ ‫الترجمة العربية‬

Elasticity ‫المرونة‬
Stress System ‫نظام الجهادأ‬
Poisson’s ratio ‫معامل بوايسون‬
Tangential Stress ‫الجهادأ المماسي‬
Transverse Stress ‫الجهادأ المستعرض‬
Transverse Strain ‫التشوه )النفعال( المستعرض‬
Normal Stress ‫الجهادأ العمودأي‬
Rigidity Modulus ‫معامل الصلبة‬
Shear Modulus ‫معامل القص‬
Hooke’s Law ‫قأانون هوك‬
Elastic Limit ‫حد المرونة‬
Plastic Point ‫نقطة اللدونة‬
Anelastic Materials ‫الموادأ اللمرنة‬
Shear resistance ‫مقاومة القص‬
Young’s modulus ‫معامل يونج‬
Compressibility ‫النضغاطية‬
Dilatation ‫تمددأ حجمي‬
Wave Propagation ‫النتشار الموجي‬
Body Waves ‫الموجات الباطنية‬
Surface Waves ‫الموجات السطحية‬
Longitudinal Waves ‫الموجات الطولية‬
Primary Waves ‫الموجات الولية‬
Compressional Waves ‫الموجات التضاغطية‬
Shear Waves ‫موجات القص‬
Transverse Waves ‫موجات مستعرضة‬
Secondary Waves ‫موجات ثانوية‬
Birch’s Law ‫قأانون بيرش‬
Seismic Velocities ‫سرع سيزميه‬
Rayleigh Waves (LR) ‫موجات رايلي‬
Love Waves (LQ) ‫موجات لوف‬
Dispersion ‫تشتت‬
Amplitude ‫سعة الموجة‬
Wavelength ‫طول الموجة‬
Frequency ‫تردأدأ‬

112
 Seismic Refraction ‫النكسار السيزمي‬
Seismic Reflection ‫النعكاس السيزمي‬
Critical Distance ‫المسافة الحرجة‬
Thickness ‫سماكة‬
Depth ‫عمق‬
Seismic Source ‫مصدر سيزمي‬
Transmitter ‫مرسل‬
Receiver ‫مستقبل‬
Geophones ‫سماعات أرضية‬
Fermat’s Principle ‫مبدأ فيرمات‬
Huygen’s Principle ‫مبدأ هايجن‬
Reflection Coefficient (Rc) ‫معامل النعكاس‬
Transmission Coefficient (Tc) ‫معامل الخآتراق‬
Acoustic Impedance ‫العائق الصوتي‬
Wavefront ‫مقدمة الموجة‬
Raypath ‫مسار الموجة‬
Snell’s Law ‫قأانون سنيل‬
Critical Refraction ‫النكسار الحرج‬
Low- Velocity – Layer ‫طبقة منخفضة السرعة‬
Hidden Layer ‫طبقة مختبئة‬
Blind Layer ‫طبق عمياء‬
Thin Layer ‫طبقة رقأيقة‬
Diffraction ‫الحيودأ‬
Delay Time ‫زمن التأخآير‬
Dipping Layers ‫طبقات مائلة‬
Green Equation ‫معادألة جرين‬
Dynamic Correction ‫التصحيح الديناميكي‬
Multiple Reflection ‫انعكاس متعددأ‬
Time- Average Equation ‫معادألة معدل الزمن‬
Faust Equation ‫معادألة فوست‬
Apparent Velocity ‫سرعة ظاهرية‬
Average Velocity (VA) ‫معدل سرعة‬
Interval Velocity (VI) ‫سرعة بينية‬
Root Mean Square Velocity ‫سرعة تربيع متوسط الجذر‬
Dix Equation ‫معادألة دأيكس‬
Data Processing ‫معالجة المعلومات‬
Cross Over Distance ‫مسافة العبور‬
Seismic Attenuation ‫تعتيم سيزمي‬

113