Cours Mthodesdeprospectionlectriques

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 27

Cours :

Les méthodes de prospection


électriques
• Établissement : École Nationale d’Ingénieurs de Sfax
• Département : Génie Géologique
• Classe : 2ème, tronc commun Géoressources et Environnement (GG2 TC)
• Responsable : Riadh AHMADI
• Code module : P16
• Année universitaire : 2020-2021

Plan
I. Introduction
II. Résistivité : Notions et caractéristiques
III. Principe des méthodes électriques
IV. Dispositifs de mesure de résistivité
IV.1) Trainée électrique

IV.2) Sondage électrique vertical

IV.3) Tomographie de résistivité électrique

V. Interprétation des données


V.1) Interprétation qualitative

V.2) Interprétation quantitative

VI. Exemple et étude de cas


VII. Exercices et TP (abaques et logiciels)

1
Cours :

Les méthodes de prospection


électriques

Objectif du module :

Ce module consiste à une initiation aux méthodes de prospection électrique


applicables à des problèmes géologiques, géotechnique et à l’exploration des ressources
naturelles (surtout en hydrogéologie).

Au cours de cet apprentissage, les étudiants vont acquérir progressivement les


compétences suivantes :

• Principe physique des méthodes électriques (loi dl’Ohm)


• Différents types de dispositifs d’acquisition (dipôles et quadripôles)
• Méthodes de prospection électriques (trainées, sondages électriques, coupes
géo-électriques, tomographie électrique)
• Techniques d’interprétation qualitative et quantitative

Prérequis :

• Indispensable : Géologie structurale, Sédimentologie,


• Souhaitable : Techniques de prospection, Géophysique (fondamentale),

2
I. Introduction
La plus ancienne utilisation des méthodes électriques date de Robert W. Fox (1830). Il a
remarqué qu’il y a une circulation du courant électrique dans les mines de cuivre.

En 1882, Carl Barus conduit des expériences de prospection électrique en introduisant des
électrodes non polarisables.

Les bases de cette technique ont été faites par Conrad Schlumberger, en 1912, pour
l’utilisation active des méthodes électriques dans la terre via des sources contrôlées (Fig. 1). Cette
technique devient une méthode commerciale utilisée, surtout, dans le domaine minier.

Figure 1 : Photo prise durant les premières compagnies d'acquisition électrique en 1912.

II. Résistivité : Notions et caractéristiques


La résistivité électrique d’une formation géologique est un paramètre de base qui reflète son
contenu en fluides (eau ou hydrocarbures).

La résistivité d’une roche dépend de la qualité de l’électrolyte (résistivité de l’eau interstitielle), de


sa quantité (porosité) dans la roche ainsi que de son mode de distribution (Isotropie).

La prospection n’atteint que des valeurs moyennes de la résistivité, relatives à des volumes
de terrain en place. Les valeurs de la résistivité sont très différentes des mesures obtenues sur
échantillon à cause des hétérogénéités locales que peuvent présenter les formations. Les valeurs de
la résistivité des eaux ou des roches varient de 0,2 (l'eau de mer) à 10000 Ohm*m (sable et gravier
sec et calcaire) d'un milieu à un autre (Tab. 1).

1
Tableau 1 : Intervalles de variation des valeurs de résistivité de différents milieux.

- Relation entre la Résistivité/Roche

La relation entre la résistivité et la nature du terrain est très étroite. Cependant toute
variation de faciès va influencer directement sur la variation de résistivité. Mais la valeur de
résistivité ne permet jamais la détermination du faciès ou l’âge de la couche.

La résistivité mesurée (moyenne) ne montre pas les variations lithologiques dans le cas
d’une couche mince empilée dans une formation de grande épaisseur (Fig. 2).

Figure 2: illustration schématique qui


illustre la négligence électrique d'une
intercalation de mince couche de sable
dans une couche d'argile d'épaisseur
importante.

2
- Anisotropie de terrain

Les résistivités des roches varient en fonction de la direction du courant qui les traverse, on
dit qu’elles sont anisotropes. Cette anisotropie est due à la structure-même des roches.
Cependant, les terrains sédimentaires sont généralement plus résistants dans la direction
perpendiculaire au plan de stratification, il s’agit alors de micro-anisotropie (Fig. 3).

Figure 3 : Illustration schématique de l'anisotropie de la résistivité d'une couche


en relation avec le plan de stratification.

III. Principe des méthodes électriques


Le principe de la méthode électrique repose essentiellement sur l’interprétation de la
résistance électrique du terrain. C’est par la mesure d’une différence de potentiel (méthodes des
potentiels), d’une part, et surtout d’une intensité, d’autre part, que l’on accède à la valeur de la
résistance.

Figure 4 : Schéma de la mise en œuvre par un quadripôle

Dans un milieu homogène de résistivité ρ et avec deux pôles A et B, l’action conjuguée de


ces derniers donnera :

• Le potentiel en M : VM = VM(A) –VM(B) = ρ I/2π . (1/MA - 1/MB)

• Le potentiel en N : VN = VN(A) –VN(B) = ρ I/2π . (1/NA - 1/NB)

3
Avec: ρ : résistivité (Ohm.m); V : potentiel (mV); I : intensité de courant (mA); MA, MB, NA et NB :
distances entre les électrodes (m).
La différence de potentiel entre M et N est :
• VM – VN = ρ I/2·. (1/MA – 1/MB – 1/NA + 1/NB)
-1
• La résistivité ρ = 2 π ∆V/ I. (1/MA – 1/MB – 1/NA + 1/NB)
 Si le sous-sol est homogène et isotrope, lors d’une prospection avec ce type de dispositif
(quadripôle AMNB), on obtiendra la résistivité vraie.
 Si par contre le sous-sol est hétérogène, on mesurera la résistivité apparente.
 Si on considère un sous-sol constitué de n couches planes, homogènes et isotropes, un tel
modèle où la résistivité est fonction de la profondeur uniquement, correspond au modèle
tabulaire ou modèle à une dimension (1D).
-1
La résistivité ρ = 2 π ∆V/ I. (1/MA – 1/MB – 1/NA + 1/NB) ; Que l’on peut écrire ρ = K. ∆V/ I
-1
Avec : K = 2π (1/MA – 1/MB – 1/NA + 1/NB) ; Où K est le facteur géométrique.
Pour un dispositif symétrique par rapport à O (milieu d'AB) le facteur géométrique K vaut :
• K = π . MA. NA/ MN
Le coefficient géométrique K est tributaire de la configuration géométrique des électrodes.
Selon leurs dispositions on peut distinguer plusieurs dispositifs dont les plus utilisés sont le
quadripôle de Schlumberger et le quadripôle de Wenner.

 Dispositif Schlumberger

Il s’agit d’un dispositif symétrique dont A, M, N et B sont alignées, O représente le centre de


symétrie (Fig. 5).

Figure 5 : Représentation schématique de la configuration Schlumberger.

On pose OA = OB =a OM = ON = b

2 2 2
Par conséquent ρa = 2 π (a – b ) /4b. ∆V/I  ρa ≈ π a /2b. ∆V/I

4
2 2
Le coefficient géométrique vaut donc : K = π a /2b [si (b / a) <<1]
 Les résultats de calcul pour les structures 1D ont conduit à la construction d’abaques, très
utilisés autre fois pour l’interprétation manuelle des sondages électriques (Fig. 6). Actuellement,
l’interprétation se fait systématiquement, à l’aide de programmes sur ordinateurs.

Figure 6 : Explication du mode d'utilisation des abaques pour l’interprétation manuelle


des sondages électriques.

 Dispositif Wenner

Il s’agit d’un dispositif dont les électrodes sont aussi alignées et réparties d'une manière
équidistante (AM=MN=NB=a) (Fig. 7). Avec K le coefficient géométrique est égal : K = 2 π a.

5
Figure 7 : Représentation schématique de la configuration Wenner.

 Dispositif dipôle-dipôle axiale

Ce dispositif est un peu différent des deux premiers (Schlumberger et Wenner) puisque les
électrodes M et N sont en dehors de A et B mais toujours dans la même ligne (Fig. 8). Chaque paire
d'électrode à une constante mutuelle de séparation égale à a.

Figure 8 : Représentation schématique de la configuration dipôle-dipôle axiale

Si la distance entre les deux électrodes internes B et N est un multiple de a (soit n x a), la
résistivité apparente sera égale a :
ρa = π n a (n +1) (n + 2) ∆V/I .

Le coefficient géométrique K est alors : K = π n a (n +1) (n + 2).

Les profondeurs d'investigations sont variables d'un dispositif à un autre (Roy, 1971 et
Barker, 1989) ainsi que le pouvoir de résolution (Fig.9). De ce fait, le choix de dispositif est fonction
de l'objectif de prospection.

6
Figure 9: Récapitulation et caractérisations des profondeurs d'investigation et du pouvoir
de résolution de différents dispositifs de prospection électrique.

IV. Dispositifs de mesure de résistivité


Les propriétés électriques d'un terrain peuvent être exploitées de manières différentes. Pour
ce la c'est en fonction de l'objectif d'investigation il y a une multitude de technique et de dispositifs.
IV.1) Trainée électrique

Il s’agit de déplacer tout le dispositif en gardant fixes ses dimensions sur le site à explorer ;
à chaque station on mesure une résistivité apparente (Fig.10). Le but est de mettre en évidence les
variations latérales de résistivités dues aux changements de faciès, la présence de failles ou des
filons…etc.
Le dispositif le mieux utilisé en traîné électrique est celui de Wenner.

Figure 10 : principe d'investigation d'une zone par trainée électrique en utilisant le dispositif Wenner.

7
IV.2) Sondage électrique vertical "SE ou SEV"

L’objectif de ce dispositif est de reconstituer le profil vertical du domaine prospecter en


profondeur (+/- AB/4), en déterminant le nombre de couches à partir de leurs résistivités et leurs
épaisseurs respectives (Fig. 11).

Figure 11 : Mise en œuvre du sondage électrique vertical

L’exécution d’un SEV sert à caractériser la variation de la résistivité du sous-sol à la verticale


en un point donné de la surface. Cette variation dépend de la profondeur pour les structures 1D.
On exécute, au même endroit, une succession de mesures en augmentant à chaque fois la distance
AB (en laissant fixe le point O : centre de AB et MN) : On augmente ainsi la profondeur
d’investigation qui est proportionnelle à la longueur AB (Fig. 11, Fig. 12). On considère que les
couches sont parallèle et isotrope (Fig. 12).

8
Figure 12 : Mise en œuvre du sondage électrique vertical, la profondeur d’investigation augmente
proportionnellement fonction de la longueur AB. On considère que les couches sont parallèles et isotrope.

En principe, tous les dispositifs (dipôles, tripôles, quadripôles, … etc) peuvent être utilisés
pour exécuter des sondages électriques. Mais dans la pratique on n’utilise que les quadripôles
symétriques, et le plus souvent le dispositif Schlumberger.
-1
ρ = 2 π ∆V/ I. (1/MA – 1/MB – 1/NA + 1/NB)
La distance MN est maintenue fixe et aussi petite que possible pendant un certain nombre
de mesures. Cela allège le travail manuel et évite le nombre des « à coups de prise ».
Mais à un certain moment, le voltage entre M et N devient faible et de l’ordre de la précision du
resistivimeter. Dans ce cas, Il faut augmenter la distance MN en respectant la condition MN <<AB
(généralement MN ≤ AB/5). Cette dernière étape effectuée est appelée débrayage ou aussi
embrayage (Fig. 13).

9
Figure 13 : Explication de débrayage (embrayage) et leurs représentations.

IV.3) Tomographie de résistivité électrique

C'est une acquisition directe d’une surface en 2D effectuée en utilisant un équipement multi
électrodes (Fig. 14). Avec ce dispositif on peu cartographier les variations de résistivité en fonction
de profondeur (Fig. 14).

Figure 14 : A) Mise en œuvre du dispositif de tomographie électrique. B) exemple d'une section de


tomographie électrique (Résistivité vs. Profondeur).
10
V. Interprétation des données
L'interprétation des données se fait en suivant deux manières successives et
complémentaires à savoir interprétation Qualitative et Quantitative.

V.1) Interprétation qualitative

L’analyse qualitative peut nous fournir des indications sur la nature des terrains à explorer.
C’est pour cette raison qu’on doit avoir recours à une approche qui peut nous indiquer sur les zones
d’anomalies tout en précisant leur nature et leur géométrie.

D'abord, après avoir recueilli les données brutes (SEV ou TE), on doit avoir recours à une
interprétation qualitative en vue d’isoler les zones d’anomalies.
Ensuite, l'établissement des cartes d’iso-résistivités apparentes pour des AB/2 croissantes.
Enfin, des indications sur l’évolution verticale et latérale des résistivités peuvent être fournies.

L'interprétation qualitative est fiable dans le cas de terrains homogènes et subhorizontaux


La tranche de sol concernée par cette interprétation est généralement comprise entre AB/10 et
AB/4. De telles approximations doivent être considérées avec prudence car la résistivité apparente
dépend de plusieurs facteurs tels que le quadripôle utilisé (le dispositif) et la nature de terrains.

Cette interprétation peut se faire également en dressant les pseudo-coupes suivant un profil
donné : à l’aplomb d’un SE, les valeurs de résistivités apparentes correspondant aux différentes
longueurs de AB sont portées et des courbes d’égales valeurs sont ensuite tracées. Cette manière
d’interprétation peut nous renseigner sur les variations latérales et verticales. Les changements
latéraux de faciès et les accidents sont facilement détectés (Fig. 15).

11
Figure 15 : Exemple d'interprétation qualitative et corrélation géo-électrique entre les différents SEV.

V.2) Interprétation quantitative

Le but de l’interprétation quantitative est de déterminer les résistivités vraies des différents
terrains traversés par le courant. Les résultats d’un SEV sont portés sous forme d’un diagramme où
l’on porte en abscisses les de AB/2 et en ordonnées les valeurs des résistivités apparentes ρa

correspondantes, avec une échelle bilogarithmique (Fig. 16).

12
Figure 16 : Diagramme bi-logarithmique pour la représentation des sondages électriques.

 Cas de deux couches:

1- On superposer le sondage effectué aux abaques à deux couches disponibles, les axes
doivent être de même dimension (Fig. 17).

2- On déplace la courbe sur les abaques jusqu'à ce qu'elle se superpose à une des courbes des
abaques (ou à une courbe imaginaire puisque tous les ρ2/ρ1 ne sont pas indiquées) (Fig. 1).

NB: Il faut garder les axes des deux graphiques bien parallèles.

13
Figure 17: Abaque utilisé pour calculer
manuellement la résistivité vraie des couches
dans le cas de deux couches.

Exemple 1: Calcul manuelle

Figure 18: A) tableaux des valeurs de résistivité apparentes. B) Représentation des valeurs de la résistivité vs.
AB/2 sur une courbe bi-logarithmique. C) Abaque utilisé pour l'interprétation.

14
Les étapes à suivre pour achever l'interprétation quantitative manuelle des données brutes
du cas d'étude sont on nombre de cinq (les étapes sont illustrées dans les figures 19A, B, C, D et E).

Figure 19 A :1ère étape, Représentation des données brutes sur la feuille bi-logarithmique.

Figure 19B : 2ème étape, On place le calque de l'abaque sur la courbe bi-logarithmique,

en cherchant la meilleure superposition.

15
Figure 19C : 3ème étape, En choisissant la meilleure superposition entre la courbe établie sur la feuille
bilogarithmique et celle de l'abaque.

Figure 19D : 4ème étape, On va projeter de la croix de référence sur l'axe d'ordonnée du diagramme bi-
logarithmique pour déterminer la résistivité réelle de la première couche. Projection de la croix de référence
sur l'axe des abscisses pour déterminer l'épaisseur de la première couche.

16
Figure 19 E : 5ème étape, On détermine la résistivité réelle du deuxième terrain suit au rapport atteint par
l'asymptote de la courbe.

Actuellement le traitement des données et l'interprétation se fait par des logicielles appropriés.

Exemple 1 : Application par logiciel (IP2WIN)

Figure 20 : Exemple d'un model établi suit à une application du programme IP2WIN du cas de la figure 18A.

17
Il suffit de faire entrer les données brutes (Fig. 18A) dans le logicielle IP2WIN (input). Le
programme va faire les inversions nécessaires et les interprétations et il va donner le model fini
(output) Fig.20.

NB. : L'ingénieur géologue doit contrôler la qualité de données et bien évidement contrôler
et interpréter le model donné par le logiciel en fonction du contexte géologique de la zone

 Cas de trois couches:

L'interprétation devient un peu plus compliquée, puisque aux deux courbes possibles qu'on
aurait dans le cas de deux couches, on passe à quatre cas possibles (Fig.21).

Figure 21 : Les différent cas possible dans le cas d'un model de trois couches.

Les abaques à deux couches ne dépendaient que de trois paramètres ρ1, ρ2/ρ1 et e1
(Fig.19).
Dans le cas à trois couches, il y a 5 paramètres ρ1, ρ2/ρ1, ρ3/ρ2, e1 et e2/e1. Les abaques
sont donc plus compliqués et plus nombreuses (Fig.21).
On présente les abaques selon les contrastes de résistivité et on trace les courbes pour
différent e2/e1.

L'interprétation se fait en ajustant la courbe expérimentale avec les abaques. On interprète


la partie gauche (AB/2 petit) de la courbe à partir de l'abaque à deux couches et on obtient ainsi e1,
ρ1 et ρ2/ρ1. On note alors le point (e1, ρ1) sur la courbe expérimentale.

18
Exemple 2: Calcul manuelle

Figure 22: A) tableaux des valeurs de résistivité apparentes. B) Représentation des valeurs de la résistivité vs.
AB/2 sur une courbe bi-logarithmique.

Exemple 2: Application par logiciel IP2WIN

Figure 23 : Exemple d'un model établi suit à une application du programme IP2WIN du cas de la figure 22A.

19
VI. Exemple et étude de cas
 1er cas : Tomographie de résistivité électrique et Trainée électrique.

 2ème cas : Coupes sériées de Tomographie de rés. électrique sur des alluvions fluviaux.

20
 3ème cas : Passage de 2 à 3 D sur des coupes orthogonaux de Tomographie électriques.

21
VII. Exercices et TP (abaques et logiciels)
TD- 1 (Trainé électrique)
Dans une région donnée, nous avons mené une campagne électrique avec AB = constante (AB
= 100 m) répartie sure une maille orthogonale équidistante (50 m). La maille est orientée N-S et E-
W. Les mesures enregistrées sont représenté sur le tableau ci-dessous:

1) Dressez les points selon leur coordonnés relatifs et projeter les valeurs sur chaque point.
2) Pour permettre l’évaluation spatiale de la répartition des résistivités de surface, dessiner la
carte d’iso-valeurs de résistivité avec une équirésistivité adéquate.
3) Interpréter les mesures observées.

TD- 2 (SEV, Analyse quantitative)


En considérant les données d’un SEV ci-contre :
1) Dessiner la courbe correspondante sur papier bi-logarithmique
2) Identifier le nombre de couches
3) Déterminer les résistivités et épaisseurs de chaque couche.

22
TD- 3 (SEV, Analyse quantitative)

On considérant les donnés d’un SEV ci-contre :


1) Dessiner la courbe correspondante sur papier bi-logarithmique
2) Identifier le nombre de couches
3) Déterminer les résistivités et épaisseurs de chaque couche.

TD- 4 (SEV, Analyse quantitative)


Un sondage électrique vertical réalisé, au cours de l’été, dans une région à couches
tabulaires et selon le dispositif Schlumberger, on a enregistré les données suivantes :

• Dessiner la courbe représentative de ces mesures sur le papier bi-logarithmique.


• Conclure le nombre de couches et leurs résistivités relatives.
• Calculer les épaisseurs et les résistivités de chaque couche
Le même point, mesuré avec le même dispositif, le même équipement et les mêmes
espacements AB, mais en période hivernale, a donné les mesures suivantes :

• Juste en observant la variation de la résistivité par rapport à la 1ère mesure, qu’est-ce qui a
changée ? et que pouvez-vous conclure ?

23
TD- 5
Dans une région agricole à structure géologique sub-tabulaire, on cherche à explorer un
aquifère gréseux situé à environ 50 m de profondeur. La prospection électrique par SEV a donné les
résultats suivants (voir les figures) :

1) Evaluer, en une seule phrase, la qualité du modèle et des résultats obtenues ?


2) Identifier le nombre de couches, leurs résistivités et leurs épaisseurs relatives.
3) Identifier la profondeur de l’aquifères gréseux prospecté.
4) Proposer une profondeur raisonnable pour exploiter cet aquifère par un forage.

TD- 6
Sur un terrain constitué de 3 couches : sol argileux à la surface, suivit par une couche à
dominance gréseuse et à potentialité aquifère, et enfin des argiles silteuses comme substratum, on
a exécuté la coupe tomographiques suivante :
1) Identifié, en argumentant par une phrase, les couches géologiques (remplir la colonne
lithologie sur la figure) et conclure : leurs épaisseur respectives au niveau du forage projeté.
2) Décrire l’aspect structural de cette coupe.
3) Que représente le corps triangulaire à forte résistivité (rouge) sur cette coupe ?

Lithologie

24
Références bibliographiques

• Barker R. D. (2001). Principles of electrical imaging.


• Ahmed Srdah AL-Zubedi (2015). Principles of Electrical Resistivity Techniques. Publisher:
LAP LAMBERT Academic PublishingEditor: Olga PomogaibogISBN: 978-3-659-78918-2.
• Chapellier, D. (1987). Diagraphie appliquée à l’hydrologie. Techniques et documentation.
Edition LAVOISIER.
• Chapellier D. (2000). Prospection électrique de surface. Cours en ligne de Géophysique.
Université de Lausane. Institut Français du Pétrole. https://www-
ig.unil.ch/pdf/doc_pro/pro_f.pdf. https://www-ig.unil.ch/acc5f.htm
• Torleif Dahlin (2001). The development of DC resistivity imaging techniques. Computers &
Geosciences 27 (2001) 1019–1029.
• Keller,G.V. and F.C. Frischknecht (1966). Electrical methods in geophysical prospecting.
Pergamon Edition.

25

Vous aimerez peut-être aussi