CH 3 Prospection Electromag.
CH 3 Prospection Electromag.
CH 3 Prospection Electromag.
La prospection par champs électromagnétiques (EM) artificiels ou naturels, générés par des
courants variables dans le temps, porte le nom de prospection électromagnétique. Comme on va
voir, les techniques électromagnétiques peuvent être variées quasiment à l’infini, de sorte qu’il ne
peut être question ici que d’en faire ressortir les principes généraux, les avantages, les inconvénients,
les limitations.
La complexité des phénomènes, et par conséquent les difficultés mathématiques d’interprétation
y sont bien plus grandes encore que lorsqu’on se sert de courants continus (DC). C’est pourquoi, si
les techniques électromagnétiques sont parfois intéressantes, très intéressantes même quand il s’agit
d’une reconnaissance rapide, d’une « détection » sommaire ou de la simple découverte de zones
d’anomalies, l’interprétation quantitative en 1D, 2D et 3D peut devenir très compliquée et est loin
d’être « intuitive ». Elles ont d’autre part un vice rédhibitoire : leur profondeur d’investigation est
limitée, d’autant plus limitée que la fréquence du champ EM est plus élevée.
Tous les appareils de prospection électromagnétique répondent à une grande variété de conduc-
teurs tant naturels qu’artificiels, qui peuvent se classer comme suit :
1. Conducteurs superficiels
– mort-terrain (terrain marécageux, argileux)
– fonds de lacs et lits de cours d’eau
– formations conductrices (argiles)
– topographie (relief).
2. Conducteurs dans la roche en place
– Graphite
– Sulfures massifs
– Magnétite massive
– Zones de cisaillement et failles
– Péridotite serpentinisée
3. Conducteurs artificiels
– Réservoirs métalliques
– Conduites et déchets métalliques
– Pipe-lines
– Voies ferrées
– Lignes à haute tension
1
1.L’induction électromagnétique
1.1 Généralités
En régime variable, champ électrique et champ magnétique sont étroitement liés, et même
littéralement indissociables. C’est de là d’ailleurs que vient le nom d’électromagnétisme. Un champ
magnétique variable engendre un champ électrique : en cela consiste l’induction électromagnétique
de Faraday. Dans un conducteur, un champ électrique crée un courant, lequel crée un champ
magnétique, c’est la loi d’Ampère.
La méthode de prospection EM fait intervenir simultanément ces trois processus physiques dis-
tincts (figure 1.1) :
1. le premier consiste en la production d’un champ magnétique primaire qui varie avec le temps ;
2. le second est la naissance de courants induits (courants de Foucault) dans tous les conducteurs
sur lesquels agit ce champ primaire ;
3. le dernier est la détection de ces conducteurs par la mesure des champs magnétiques secondaires
créés par les courants de Foucault.
T r a
R é
N i
C o u
F C o o u
C h
C h a
L’amplitude des courants induits dans un corps conducteur dépend de plusieurs facteurs, à peu
près équivalents, qui sont :
– les propriétés électriques du conducteur ;
– les dimensions et la forme du conducteur ;
– la fréquence du champ primaire ;
– l’emplacement du conducteur par rapport aux instruments géophysiques.
Cette discussion est résumée à la figure 1.2 où les éléments clés y sont présentés de façon
schématique. En vertu de la loi de Lenz (règle de la main droite), les courants de Foucault cir-
culent dans le conducteur de telle sorte que le champ magnétique ainsi créé (champ secondaire)
s’oppose, sur les surfaces du conducteur, au champ inducteur (champ primaire).
2
A) B)
A
H secondaire
H = Hp exp(iωt)
H = Hp exp(iωt)
Courant de Foucault
Figure 1.2: Induction électromagnétique. A) Vue en perspective. B) Vue suivant la coupe A-B.
Figure 1.3: Modèle simple de la boucle
conductrice.
Le phénomène de l’induction d’un champ magnétique variable avec un conducteur peut être
le mieux visualisé à l’aide d’un modèle simple, comme une boucle conductrice fermée (figure 1.3).
Même si de tels gisements n’existent pas dans la nature, nous savons que les réponses générales de
notre « gisement » ultra simple et du conducteur réel seront quantitativement semblables.
Si nous prenons cette boucle conductrice de rayon a et l’exposons à un champ magnétique
primaire (variable et uniforme, et dirigé le long de son axe) donné par
H = Hp cos(ωt)
2a
H = Hp cos (ωt)
dΦ
V = − ,
Z dt
Φ = B ~
~ · ds,
s
I
V = E ~
~ · dl
~
B ~
= µH,
dI
= RI + L .
dt
La solution d’une telle équation différentielle est
0 1 ωτ
I= sin ωt − cos ωt
R 1 + ω2 τ 2 1 + ω2 τ 2
avec τ = L/R.
Au centre du conducteur, on peut calculer le champ magnétique secondaire créé par ce courant
à l’aide de la loi de Biot-Savart
I
Hs |0 = .
2a
Donc
0 1 ωτ
Hs |0 = sin ωt − cos ωt
2aR 1 + ω 2 τ 2 1 + ω2 τ 2
Avec 0 = µωSHp , on obtient que
Hs µωS 1 ωτ
= sin ωt − cos ωt .
Hp 0 2aR 1 + ω 2 τ 2 1 + ω2 τ 2
4
Or τ = L/R et en sortant 1/(1 + ω 2 τ 2 )1/2 de la grande parenthèse, on obtient
Hs µS ωτ ωτ 1
=− √ √ cos ωt − √ sin ωt
Hp 0 2aL 1 + ω2 τ 2 1 + ω2 τ 2 1 + ω2 τ 2
√ 1 √ ωτ et φ = tan−1 1
Posons sin φ = 1+ω 2 τ 2
, cos φ = 1+ω 2 τ 2 ωτ , alors
Hs
µS ωτ
=− √ cos φ cos ωt − sin φ sin ωt
Hp 0 2aL 1 + ω2 τ 2
| {z } | {z }
Partie en phase Partie en quadrature
D’où
Hs µS ωτ
=− √ cos(ωt + φ)
Hp 0 2aL 1 + ω2 τ 2
µS √ ωτ
Posons maintenant G = 2aL et F (ωτ ) = 1+ω 2 τ 2
. Alors
Hs
= −G F (ωτ ) cos(ωt + φ)
Hp 0
Les divers aspects de l’exposition mathématique présentée ci-haut sont résumés sous forme gra-
phique à la figure 1.4.
H = Hp cos (ωt)
Champ primaire
0 t
φH
H = Hs cos (ωt+φ)
Champ secondaire
0
t
π−φH
Champ secondaire
(en phase)
0 t
Champ secondaire
(en quadrature)
0 t
5
Remarques :
1. Au centre du conducteur, la direction du champ secondaire est opposée à celle du champ
primaire (signe négatif) ;
2. L’amplitude du champ secondaire est proportionnelle à l’amplitude du champ primaire qui
recoupe le conducteur ;
F
3. L’amplitude du champ secondaire est fonction de :
(a) la géométrie du conducteur G = µS/2aL,
(b) du paramètre d’induction du conducteur α = ωτ = ωL/R. Il est à noter que ce dernier est
intimement lié à la fois aux paramètres électriques du conducteur (L/R) et à la fréquence
d’excitation (ω)
4. Le champ magnétique secondaire, dans le cas général, n’est pas en phase avec le champ primaire
(déphasage de φ).
Avant de continuer la discussion des caractéristiques du champ magnétique secondaire, il est
utile de mettre en graphique les relations mathématiques entre celui-ci et les propriétés électriques
du conducteur.
Ainsi, la figure 1.5 illustre la relation du déphasage et de l’amplitude du champ secondaire
en fonction du paramètre d’induction α = ωτ . Par contre, la figure 1.6 présente la variation de
l’amplitude des deux composantes phasorielles (composante en phase P et composante en quadrature
Q) du champ magnétique secondaire en fonction, toujours, du paramètre d’induction.
0.50 π 1 .
0 .
0 .
0 .
0.00 10−2 10−1 100 101 10 2
1
1
α
Figure 1.5: Relation du déphasage et de l’amplitude du champ secondaire
en fonction du paramètre d’induction α.
1.00
Partie en phase
Partie en quadrature
0.75
0.50
0.25
0.00
10−2 10−1 100 101 102
α
Figure 1.6: Variation de l’amplitude des deux composantes phasorielles P et
Q du champ magnétique secondaire en fonction du paramètre d’induction
α.
6
Finalement, la figure , combine tous les éléments des deux illustrations précédentes. Un tel
diagramme est habituellement appelé diagramme d’Argand et est très utile dans l’interprétation
quantitative de données provenant de la technique électromagnétique à cadres horizontaux
Q = H sq
GHp
0.5 (P-1/2)2 + Q2 = 1/4
!" = 1
0.4
0.3
!" = 1/3
2 !" = 3
2 +Q
P
0.2 =
!"
2 "2
)= 1 + !
F(!"
0.1
1
#(!") = tan-1( ) = tan-1( Q )
!" P
0
0 0.1 0.2 0.3 0.4 0.6 0.7 0.8 0.9 1.0
P = -Hsp
Hp = Amplitude du champ primaire
GHp
µS
G=( ), Facteur géométrique
2aL
7
2 Les équations de Maxwell
En unités SI, et dans des milieux conducteurs, les équations de Maxwell s’écrivent
∂H~
~
∇×E = −µ ,
∂t
~
∇ ×H ~
= σ E. )
~
Ces deux équations, en impliquent deux autres car on sait que pour un vecteur quelconque A
~ = 0.
∇·∇×A
~
∇·H = 0,
~
∇·E = 0.
~ ~
La solution des équations simultanées de Maxwell s’obtient à l’aide de ces deux équations auxi-
liaires pour la divergence, et à l’aide de l’identité vectorielle suivante ~
∇× ∇× A ~ =∇ ∇·A ~ − ∇2 A. ~
L’effet de peau
Dans un milieu homogène semi infini dont la surface est plane, un champ magnétique alternatif
peut se distribuer uniformément. L’équation de Helmholtz nous indique que le champ magnétique
doit varier le long d’au moins une direction spatiale, car ∂H/∂t est différent de zéro. Si nous admet-
tons l’uniformité le long des axes x et y (sol stratiforme), nous aurons alors
∂2H ∂H
− µσ =0
∂z 2 ∂t
Si on suppose que le champ magnétique est donné par
(
H = Hx (z, t) = Ho eiωt+mz),
∂ 2 Hx
= − m2 H x ,
∂z 2
7
et la dérivée selon t est
∂Hx
= iωHx .
∂t
Ainsi, nous avons p
− m2 = iωµσ ou m = ±(1−i) ωµσ/2.
Puisque le champ ne peut croı̂tre à l’infini, le signe positif de la racine est éliminé, et la solution est
alors
Hx = Ho e−z/P ei(ωt−z/P )
` r
2
P = .
µωσ
Le champ diminue donc selon une exponentielle décroissante en fonction de la profondeur. En unités
pratiques et dans un milieu non-magnétique (µ0 = 4π × 10−7 ), on utilise la relation
r
ρ
P ≈ 500
f
3 .1 Introduction
Le champ magnétique inducteur, soit le champ primaire, nécessaire à toute prospection électroma-
gnétique est le plus souvent créé par un courant alternatif circulant dans une boucle. La dimension
de la boucle varie selon les besoins de la méthode de prospection utilisée. Néanmoins, parfois on
`
profite d’autres sources conventionnelles, comme la méthode T.B.F. (V.L.F. en anglais) qui utilise
des émetteurs radio dont la fréquence d’émission se situe dans la gamme de 10 à 20 kHz. Une
autre méthode, la méthode AFMAG est basée sur les champs magnétiques de fréquences audio
engendrés par des éclairs produits lors des orages électriques. Finalement, certaines méthodes, comme
le TURAM, peuvent utiliser le champ primaire engendré par le courant alternatif qu’on fait circuler
dans un câble dont les deux extrémités sont reliées au sol.
3 .2 La loi de Biot-Savart
~ = I ~
dH dl × ~
r
4πr3
Idl sin θ
|dH| =
4πr2
~ et ~r. La direction du champ est déterminée directement de l’équation
~ est perpendiculaire à dl
où dH
vectorielle ou encore par « la règle de la main droite » appliquée à la forme scalaire. Dans l’illustration,
le vecteur du champ magnétique entre dans le plan du dessin.
En unités SI, le champ magnétique se chiffre en Ampère/mètre (A/m) avec l’équivalence suivante
en unité utilisée en prospection magnétique
3 .3 Les émetteurs
Iz l(2π/λ) Iz l
Hx = = λ
4πr 2r
où l est la longueur de l’antenne. Cette équation implique que le champ ainsi produit est raisonnable-
ment uniforme (décroissance en 1/r), mais le plus important est que l’on trouve le champ primaire
partout horizontal et dirigé dans une direction perpendiculaire à la direction de propagation.
Il est à noter que les champs magnétiques naturels, utilisés dans la méthode AFMAG, démontrent
le même comportement.
10
x = x1 x=0 dx x = x2
a
r
Le champ magnétique d’un fil droit de longueur arbitraire s’obtient par intégration à partir de
loi de Biot-Savard. Ainsi, pour le fil dessiné à la figure .
Idx sin θ
dH = ,
4πr2
mais sin θ = cos β, avec θ = π − (π − π/2 − β) = π/2 + β, et x = a tan β.
Soit dx = a sec2 β dβ avec sec2 β = 1 + tan2 β, et r2 = a2 (1 + tan2 β). Ainsi,
I a cos β sec2 β dβ
dH =
4π a2 (1 + tan2 β)
I sec2 β
= cos β 2 dβ,
4πa sec β
soit
Z β2
I
H = cos β dβ
4πa β1
I β
= |sin β|β21 ,
4πa
I
H= [sin β2 − sin β1 ]
4πa
soit " #
I x2 x1
H= p −p .
4πa a2 + x22 a2 + x21
Il est à noter que si le fil est long, soit (x1 , x2 ) a, la formule se simplifie à H = I/2πa,
démontrant ainsi la décroissance uniforme comme 1/r. Par contre, si le fil est court, soit (x1 , x2 ) a,
on retrouve la loi de Biot-Savart, c’est-à-dire
I (x2 − x1)
H= .
4πa2
On remarque aussi le fait que le champ magnétique est toujours perpendiculaire au plan qui
contient le fil et le point d’observation. La variation spatiale du champ magnétique d’un fil long peut
être mieux appréciée à l’aide d’un diagramme qui présente le fil en section (figure ).
On note que sur le plan z = h, le champ magnétique sera donné par
I I
H= = 1/2
2πa 2π (y + h2 )
2
11
y
Y z=h
z Z
a
H
ψ z=0
y
x
1.0
Y/Ymax
0.5
-4 -2 2 4
` Z/Ymax
12
L’émetteur dipolaire
Nous noterons que dans le plan de la boucle (z = 0, s = (x2 + y 2 )1/2 ), le champ magnétique a
partout la même direction que l’axe du cadre (les autres composantes sont égales à zéro) et décroı̂t
comme 1/r3 . La même chose est vraie sur l’axe du cadre sauf que cette fois-ci, le champ magnétique
est orienté dans le sens inverse et est deux fois plus fort qu’en un point également distant du centre,
mais situé dans le plan du cadre.
De plus, on remarque la dépendance de l’amplitude du champ par rapport à l’aire du cadre
(a2 , voir figure ). De fait, le champ magnétique d’un petit cadre ne dépend pas de sa forme
géométrique exacte mais, en général, il n’est lié qu’à l’aire du cadre. De tels émetteurs portent
l’appellation dipolaire parce que la distribution de leur champ magnétique ressemble fortement à
celle d’un petit aimant dont l’orientation est la même que celle de l’axe du cadre.
1 Classification
Comme on pouvait s’en douter, les méthodes EM sont extrêmement variées, suivant qu’on choi-
sit tel ou tel type de source, suivant que les sources sont fixes ou mobiles, suivant qu’on mesure
l’orientation du champ résultant ou seulement les parties en phase et quadrature d’une de ses com-
posantes et enfin suivant la fréquence adoptée.
Il est néanmoins utile, avant de procéder à des sous-classifications de séparer les méthodes cou-
rantes en deux catégories principales. Les deux catégories seront liées à la nature du champ primaire
qui peut être considéré soit comme uniforme soit comme dipolaire.
13
2. Dans la catégorie D (champ dipolaire) nous mettrons
– toutes les méthodes à émetteur mobile comme
– EM Vertical
– EM Horizontal
– toutes les méthodes aéroportées.
La base physique de notre choix de ces deux catégories principales est directement liée à la
distribution des courants de Foucault dans la cible usuelle de la prospection EM. Celle-ci prend la
forme d’un conducteur tabulaire dont les dimensions latérales excèdent sa dimension transversale
(épaisseur) de quelques ordres de grandeur.
Dans le cas où on utilise un dispositif d’excitation à champ magnétique uniforme, les courants
induits suivent, à peu près, le périmètre du conducteur (figure ). La phase du courant est plus
ou moins constante tout le long de son parcours. Dans ce cas, la phase du champ secondaire ne varie
que peu sur le plan d’observation. Le paramètre d’induction du conducteur n’est donc lié qu’à ses
dimensions physiques et se chiffre par
2el
αu = µo σωel =
p2
q
2
où p = µo ωσ et e est l’épaisseur du conducteur.
Hp
Hp
Hp I
Excitation uniforme
Par contre, dans le cas où l’excitation se fait par un champ dipolaire, la distribution des courants
induits est très complex. On note surtout que son amplitude et sa phase varient d’un point à l’autre.
Il s’en suit que les paramètres du champ secondaire, et par conséquent, le paramètre d’induction
applicable à ce type de situation, seront fonction de l’écart émetteur-récepteur.
De fait, l’expérience démontre que dans les cas où l’étendue latérale du conducteur est plus grande
que la distance récepteur-émetteur, le paramètre d’induction se chiffre par
2el
αp = µo σωel =
p2
14
2 Les méthodes de la catégorie U
La méthode T.B.F.
Comme nous l’avons déjà indiqué, cette méthode utilise un émetteur lointain qui fournit un
champ primaire horizontal et uniforme. À cause de sa fréquence élevée (10 – 20 kHz), le champ
secondaire produit par des conducteurs intéressants est d’habitude en phase avec le champ primaire.
Le champ résultant démontre par conséquent une polarisation linéaire. La détection d’un conducteur
se fait à l’aide des mesures d’angle d’inclinaison du champ résultant. Puisque le champ primaire est
horizontal, toute inclinaison du champ résultant indique la présence d’un corps anomal.
Comme on le démontre au dessin , cette méthode se pratique avec des cadres horizontaux et
coplanaire. L’un sert comme émetteur et l’autre comme récepteur.
On fait des mesures (relatives par rapport à la composante verticale du champ primaire au
point d’observation) de l’amplitude et de la phase de la composante verticale du champ secondaire
en traversant le terrain sur des lignes orientées perpendiculairement à la direction de la structure
géologique. Les mesures se font d’une façon directe à l’aide d’un dispositif électronique spécialement
destiné à cette fin. Celui-ci, en même temps, compense la présence du champ primaire au point
d’observation.
déplacement
Séparation constante
émetteur récepteur
Conducteur
Alimentation
et source
Câble de
3’ référence
3’ Amplificateur
200’
Émetteur Récepteur
´
En pratique, on utilise des fréquences d’excitation autour de 1 kHz avec des séparations des
éléments entre 25 à 200 m. L’équipement dont les éléments sont illustrés à la figure 3est léger,
et alimenté par batterie.
15