Magnetic Surveying

Dr. Laurent Marescot

[email protected]

1
 Introduction
Magnetic surveying…

Investigation on the basis of anomalies in the Earth´s

magnetic field resulting from the magnetic properties of
the underlying rocks (magnetic susceptibility and
remanence)

2
 Application
• Exploration of fossil fuels (oil, gas, coal)
• Exploration of ore deposits
• Regional and global tectonics
• Large scale geological structures, volcanology
• Buried conductive objects (cables, drums)
• Unexploded ordnance (UXO)
• Archaeological investigations
• Engineering/construction site investigation

3
 Structure of the lecture
1. Equations in magnetic surveying
2. Geomagnetic field (refresher)
3. Magnetic properties of rocks (refresher)
4. Survey strategies and interpretation
5. Conclusions

We will not talk about magnetic properties at an atomic scale,

paleomagnetics or the magnetic structure of the Earth. These
notions were developed last year. We will focus on magnetics
for environmental and engineering applications and
emphasize links with gravimetry.
4
1. Equations in magnetic surveying

5
 Basic magnetic theory
Let us first define the following terms…
G
F magnetic force
G
B magnetic induction field
G
H magnetic field strenght
G
M magnetic moment
G
Ji intensity of induced magnetization
k magnetic susceptibility

6
 Basic concepts

Within the vicinity of a bar

magnet a magnetic flux is
developed which flows from
one end of the magnet to the
other (the poles of the magnet).

7
This flux can be
mapped by a small
compass needle
suspended within it.
Similarly, a magnet
aligns in the flux of
the Earth´s magnetic
field.

8
 Force between two poles
The force F between two magnetic poles of strengths m1 and m2:

G μ0 m1 m2 G
F= r
4πμr r 2

r = (x2 − x1 )2 + ( y2 − y1 )2 + (z2 − z1 )2

magnetic permeability of vacuum : μ 0 = 4π 10 −7 Vs/Am

μ
relative magnetic permeability : μ r =
μ0
9
 Magnetic induction field
The magnetic induction field is expressed as the force
created by a pole m and applied on a unitary pole m1

G
G F μ0 m G
B= = r
m1 4πμr r 2

2
in Vs/m or Tesla (T)
10
Dipole

11
 Magnetic moment
The magnetic moment is the vector joining the two poles
–m and +m (at distance l apart)

G G
M = ml r

2
in Am
12
 Induced magnetization
When a magnetic material is placed in a
magnetic field, elementary dipoles in
the material align in the direction of the
field.

The resulting magnetization gives rise

to an additional magnetic field in the
region occupied by the material.

13
Intensity of induced magnetization
The intensity of induced magnetization is defined as the
magnetic moment per unit volume of material
G
G M
Ji =
v

in A/m

14
 Magnetic susceptibility
The magnetic susceptibility k is dimensionless and
determines the degree of magnetization of material

G
Ji
k= G
H
H simply describes how B is modified by the magnetic
polarization (or magnetization M).
In a non-polarizable body, H can be regarded as simply a
computational parameter proportional to B G G
B = μ0 H 15
 Magnetic induction
When a magnetic material is placed in a magnetic field, the
resulting magnetization gives rise to an additional
magnetic field in the region occupied by the material.
Within the body, the total magnetic induction is given by:
G G G G G G G
B = μ 0 H + μ 0 J i = μ 0 H + μ 0 kH = (1 + k )μ 0 H = μ r μ 0 H
G
B in Vs/m 2 or T (Tesla)

magnetic permeability of vacuum : μ 0 = 4π 10 −7 Vs/Am

μ
relative magnetic permeability: μr = 16
μ0
 Units of magnetism
The unit used in geomagnetic surveys is the Tesla

• 1 Tesla = 1 T = 1 N/Am
• 1 nT = 10-9 T=1 γ = 10-5 Oersted

• c.g.s unit:
1 gauss (G)=10-4 T
1 gamma (γ)=10-5 G
17
2. Geomagnetic field

18
 Geomagnetic elements
The geomagnetic elements are…

• Inclination
• Declination

19
20
Simplified model for the Earth
More complex than the
gravity field: irregular
variations with latitude,
longitude and time

Inclination varies
depending on the
hemisphere

Geocentric dipole is
inclined at about 11.4°
21
22
Simplified model for the Earth
• Dipole field in first approximation
• ~ 60 000 nT (poles)
• ~ 30 000 nT (equator)
• ~ 47 000 nT in Switzerland

23
Changes in the geomagnetic field
• The exchange of dominance between the cells produce the
periodic changes in polarity imaged in paleomag studies

• Slow variations in the circulation patterns within the core

produce temporal changes in the geomagnetic field
(secular variations, e.g. gradual rotation of the magnetic
pole)

24
The magnetic pole is moving in time

25
Change in the magnetic field intensity with time

26
Origin of the geomagnetic field

• 99 % from the Earth (94 % dipole field + 5 % non-

dipole field)

• 1 % current in the ionosphere (diurnal variations,

magnetic storms)

27
Not a remanent origin
(temperature too high).

Dynamo action produced by

the circulation of charged
particles in couples
convective cells within the
outer, fluid, part of the
Earth's core.

28
29
 Diurnal variations
• Variations of external origin. Results from the magnetic
field induced by the flow of charged particles within the
ionized ionosphere towards the poles
• Movements in ionosphere:
Difference in temperature in atmosphere
Sun-Moon attraction
• Varies with latitude and seasons (max. in summer, max in
polar regions)
• Smooth variations. Amplitude 20-80 nT

30
 47230

47225

47220
Champ magnétique terrestre total (nT)

47215

47210

47205

47200

47195

47190

47185

47180
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24

temps (h)

Observatory of Neuchâtel du 1 mars 2000 (1 sec) 31

 Magnetic storms
• Associated with intense solar activity, results from the arrival
in ionosphere of charged solar particles
• Less regular than diurnal variations. Amplitude up to 1000 nT!
• No magnetic surveys during storms (impossibility of
correcting the data)

32
29 octobre at 6h 14 Period from 18 octobre to 5 novembre 2003
Perturbation more than 350 nT
Observatory of Neuchâtel (Suisse)
33
 Geomagnetic Reference Field
• The International Geomagnetic Reference Field (IGRF)
defines the theoretical undisturbed magnetic field at any
point on the Earth´s surface in simulating the observed
geomagnetic field by a series of dipoles

• This formula is used to remove from the magnetic data

those magnetic variations attributable to this theoretical
field

http://www-geol.unine.ch/GEOMAGNETISME/HomePage.html

34
3. Magnetic properties of rocks

35
 Rock magnetism
The measured total magnetic field is the sum of the
geomagnetic field and the remanent magnetic field

36
 Rock magnetism
• All substances are magnetic at the atomic scale. Each atom
acts as a dipole due to both the spin of its electrons and the
orbital path of the electrons around the nucleus

• Two electrons can exist in the same state provided their

spins are in opposite directions (paired electrons). In this
case their spins cancelled. When unpaired electrons are
present, a magnetic moment at the atomic scale appears

• Paired and unpaired electrons are mainly at the origin of

the various magnetic rock properties
37
 Rock magnetism
G G G G
B = μ0 H + μ0 J i = (1 + k ) μ0 H

• Diamagnetic: k<0
• Paramagnetic: k>0
• Ferromagnetic (e.g. iron), ferrimagnetic (e.g. magnetite)
and antiferromagnetic (e.g. haematite)

38
Magnetic properties of rocks
G G
B = (1 + k ) μ0 H
Magnetic
properties of rock
depend mainly on
the concentration
size, shape and
dispersion of
magnetite

39
Magnetite content
 Remanent magnetization
The strength of the magnetization of ferro and
ferrimagnetic material decreases with temperature and
disappears at the Curie temperature (for most of the rocks
about 500 oC, i.e. to a depth of 40 to 50 km).

Origin of remanent magnetization:

• Thermoremanent magnetization
• Detrital remanent magnetization
• Chemical remanent magnetization
• Viscous remanent magnetization

These notions were developed in last year lectures… 40

4. Survey strategies and interpretation

41
 Magnetic surveys
How to proceed…

(1) Data measurements, basis measurements, data location

(2) Calculation of the theoretical field (IGRF)
(3) Calculation of the geomagnetic anomalies, reduction
(4) Removal of the regional trend
(5) Modeling, inversion
(6) Interpretation

42
 Magnetic surveying instruments
Two types of
magnetometers are
frequently used in
magnetic surveying:
• Proton magnetometer
• Optically pumped
magnetometer
• Other device: fluxgate
magnetometer

Precision required: about ± 0.1 nT (about one part in 5x106 of

the background field) 43
Proton magnetometer

44
 Proton magnetometer
ω = 2π f = γ p H
γ p is the gyromagnetic ratio of the proton (constant)
Precession frequency f ≅ 2000 Hz
2π f
H total = ≅ 23.49 f
γp

• Sensitivity about ± 0.1nT

• Frequency measurement 2-3 s
45
Optically pumped magnetometer
• A glass cell containing an evaporated alkali (e.g. Cesium)
is energized by a light of particular wavelength.

• Measurement principle based on the partition of valence

electrons into different energy levels.

• Very rapid and sensitive measurements: used in

gradiometers

Sensitivity < 0.01nT

Frequency measurement 0.1 s
46
Magnetic gradiometers
Meters are used in pairs
to measure either
horizontal or vertical
magnetic gradients.
Useful in shallow
geophysics to resolve
complex anomalies.
Regional and temporal
variations are
automatically removed.

47
 Reduction of magnetic data
The main corrections are…

• Diurnal variation correction

• Elevation and terrain correction
• Geomagnetic correction

48
Diurnal variation correction
• Loop to a reference basis
(tedious…)

• Use a fixed magnetometer

located at the basis to
correct the data collected
with a second
magnetometer

• Use the record of a

regional magnetic
observatory
49
Using a basis: some considerations

Remember: for gravity, the basis readings are taken both to

correct for the drift and the tidal effects.

In magnetic, only for the diurnal effect since magnetometers

do not drift!

50
Elevation and terrain corrections

The gradient of the magnetic field is only some 0.03 nTm-1

the pole and -0.015 nTm-1 at the equator: no elevation
correction is applied for ground surveys.

The terrain correction is very difficult to applied (generally

rarely applied) since we need to know about the magnetic
properties of the topographic features.

51
 Corrections for aeromagnetic meas.
• To be considered in case of aeromagnetic measurements

• The reduction to a datum is applied for local measurements

with steep topography. The measurement located at z=h is
reduced to a datum z=0 using:

δZ
Z ( x, y, 0) = Z ( x, y, h) − h
δh z =h

52
Drift, secular variations and storm
• Drift: fluxgate and proton magnetometers do not drift
• Secular variation are yearly variations. Too slow for a
influencing a survey.
• Magnetic storm: stop the survey!

Moreover…
…do not carry out magnetic surveys in the vicinity of
metallic objects such as railway lines, cars, fencing !!!
… as an operator, do not carry metallic objects!!!

53
 Latitude (geomagnetic) correction
• Equivalent of the latitude correction in gravimetry (reference
ellipsoid)

• We can use the International Geomagnetic Reference Field

(IGRF), updated every 5 years, which defines the theoretical
undisturbed magnetic field at any point of the Earth surface.
Warning: the IGRF is imperfect and in areas remote from
observatories can be substantially in error!

• Alterative method for small surveys: use a trend analysis, where

the regional field is approximated by a linear trend
54
 Residual magnetic anomaly

Complex shape: anomaly not only positive or negative like in

gravity surveys!
55
More complex than
gravity anomalies
(vary not only in
amplitude but also in
direction)

56
Bodies with identical
shapes and intensity
of magnetization can
give rise to very
different magnetic
anomalies depending
their latitude

57
Interpretation of magnetic surveys is mainly qualitative (maps)
58
 Interpretation
Like for gravity, we can use…

• Direct interpretation
• Indirect interpretation and automatic inversion

59
Direct interpretation

60
 Simple geological structures
• Ball: Compact bodies (salt domes, karst)
• Horizontal cylinder: paleo-valleys, tunnel, karst, cables
• Vertikale cylinder : volcanic structures, karst

Body Anomaly Depth

( )( ) z = 2.00 x1/ 2
Ball 4
Fz = 100 πR 3 ΔJ z 2 z 2 − x 2 x 2 + z 2
5/ 2
(magnetic dipole)
3
Horizontale cylinder
(line of dipoles) ( )(
Fz = 200πR 2 ΔJ z z 2 − x 2 x 2 + z 2 )
2
z = 1.75 x1/ 2
Vertical cylinder 1 z = 1.30 x1/ 2
Fz = 100πR 2 ΔJ z
(x )
(magnetic monopole) 3/ 2
2
+ z2
61
 Other interpretation techniques
Other techniques:

• Euler deconvolution: a complex but more rigorous method

of determining depth to magnetic sources

• Reduction to the pole: simplify anomaly, produce anomaly

that are axisymmetric

62
 Indirect interpretation
• Same approach than in gravimetry (improvement of a
initial model, see picture)
• Automatic inversion useful since anomalies are complex
• Model built using a series of dipoles (sum of positive and
negative poles)

63
Ambiguity in interpretation

64
 Comparison grav/mag 1

• Magnetic properties of the rocks disappear at about 20 to

40 km depth (Curie temperature)

• Variations of magnetic permeability over several orders of

magnitude, density over only a range of 20-30%

• Density is a scalar, intensity of magnetization is a vector

65
 Comparison grav/mag 2
• 2:1 length-width ratio sufficient to validate 2D
approximation in gravimetry, but 10:1 for magnetics!

• Survey faster and simpler than gravimetry, since no

leveling required

• The magnetic anomalies are asymmetric depending on the

latitude! The magnetic anomalies are more complex than
the gravity anomalies

66
Examples

67
68
69
70
71
72
73
74
75
76
77
Unexploded ordnance (UXO)

78
5. Conclusions

79
 Advantages
• Simple
• Fast
• Cost effective
• No artificial source required
• Good qualitative tool for mapping

80
 Drawbacks
• Very sensitive to non-unicity in the modeling solutions
• Mainly qualitative
• Very sensitive to metallic fences, rails (difficult to use in
urbanized regions)

81